In many of our classical textbook toy problems (mass moving in a certain potential), yes there is a ‘hidden’ Riemannian metric (because we’re talking about massive particles) which is how we formulate the kinetic term. But, abstractly, no it is not necessary to have a Riemannian metric in order to define a symplectic form or to formulate Hamiltonian mechanics.
You write
In the Hamiltonian formalism (working on flat space in Cartesian coordinates), the positions $r^i$ transform as vectors, while the canonical momenta $q_j := \frac{\partial \mathcal{L}}{\partial \dot{r}^j}$ transform as one-forms.
Well, that’s just confusing. The $r^i$ are not vectors or covectors or anything. They are coordinate functions. I have written about this in several answers, so you may want to read some of them and go down the linked rabbit hole starting from Help with geometric view of conjugate momenta and Legendre transformation (which has all the relevant links inside).
Next, regarding your middle paragraphs, you’re right, you can’t talk about symmetry or skew-symmetry without first fixing an isomorphism between the vector space and its dual. But this is not relevant in this discussion, and I think you’re only bringing it up out of confusion.
You write
Therefore, it isn't obvious that the symplectic form can directly input canonical momenta, which are one-forms
and
This Wikipedia page seems to imply that there is a way to define symplectic forms on smooth manifolds without a metric structure, but it's way over my head and I don't understand it…
So, it seems to me you’re misunderstanding the spaces things are defined on.
What we’re saying is you start with a smooth manifold $Q$. Then, on $T^*Q$, there is a natural 1-form $\theta$ (note this is a much more complicated beast than a 1-form on $Q$. We have $\theta\in \Gamma(T^*(T^*Q))$, i.e it is a mapping $\theta:T^*Q\to T^*(T^*Q)$ such that composing with the projection $\pi_{T^*(T^*Q)\to T^*Q}$ gives you the identity mapping on $T^*Q$). I emphasize again, $\theta$ is NOT a 1-form on $Q$, which is what your wording seems to suggest your interpretation is, instead it is a 1-form on the larger space $M:=T^*Q$. The canonical symplectic form is then given by the exterior derivative $\omega:=d\theta$. This is a 2-form on $M=T^*Q$, giving us the symplectic manifold $(M=T^*Q,\omega)$.
Now, I should mention that not every symplectic manifold arises as the contangent bundle of some manifold. For example, the complex projective spaces $\Bbb{C}P^n$ have standard symplectic structures (by using the standard one, $\sum_{j=1}^{n+1}dx^j\wedge dy^j$, on $\Bbb{C}^{n+1}=\Bbb{R}^{2n+2}$, pulling it back to the sphere $S^{2n+1}$, and then noting that it is invariant under the $U(1)$-action, and hence descends to a closed 2-form, i.e a symplectic form, on the quotient $S^{2n+1}/U(1)=\Bbb{C}P^n$). Complex projective spaces are compact, while cotangent bundles are non-compact, hence they provide examples of symplectic manifolds which are not cotangent bundles. By the way, a concrete special case is when $n=1$, in which case one has the well-known diffeomorphism $\Bbb{C}P^1\cong S^2$, meaning that the sphere $S^2$ has a symplectic structure. This can even be described explicitly as follows:
\begin{align}
\omega_{S^2}=\iota^*(x\,dy\wedge dz+y\,dz\wedge dx+z\,dx\wedge dy),
\end{align}
where $(x,y,z)$ are the usual Cartesian coordinate functions on $\Bbb{R}^3$, and $\iota:S^2\hookrightarrow \Bbb{R}^3$ is the canonical injection. In other words, if you let $(\xi,\eta,\zeta)$ be the functions obtained by restricting the functions $(x,y,z)$ from $\Bbb{R}^3$ to $S^2$, then $\omega_{S^2}=\xi\,d\eta\wedge d\zeta+\eta\,d\zeta\wedge d\xi+\zeta\,d\xi\wedge d\eta$. This is none other than the area 2-form of the sphere (so it is also not exact, i.e this $\omega$ doesn’t arise from a $\theta$ like in the cotangent bundle case).
As you can see from these concrete examples, I never used any Riemannian metric to define the $\omega$ (I didn’t assume any Riemannian metric on $Q$ or $T^*Q$, yet I defined $\omega$ on $T^*Q$. Likewise, I defined $\omega$ on $\Bbb{C}P^n$, and especially on $S^2$, all without referencing a Riemannian metric).
Next, if you have a Poisson manifold $(P,\{\cdot,\cdot\})$, for which symplectic manifolds $(M,\omega)$ serve as basic examples, this again has nothing to do with a Riemannian metric anywhere. The bracket is a skew-symmetric $C^{\infty}(P)$-bilinear mapping $C^{\infty}(P)\times C^{\infty}(P)\to C^{\infty}(P)$, and hence it gives rise to a bivector field $\pi$ on $P$, i.e $\pi$ is a section of the bundle $\bigwedge^2(TP)$, i.e $\pi:P\to\bigwedge^2(TP)$ such that it projects back to the identity on $P$ (you can think of it as a $(2,0)$ tensor field on $P$ satisfying a bunch of properties). But again, this section has been defined without referencing any Riemannian metric.