Here is a simple proof I thought, tell me if anything is wrong.
First claim. Let $k$ be a field, $V$ a vector space of dimension at least the cardinality of $k$ and infinite. Then $\operatorname{dim}V^{*} >\operatorname{dim}V$.
Indeed let $E$ be a basis for $V$. Elements of V* correspond bijectively to functions from $E$ to $k$, while elements of $V$ correspond to such functions with finite support. So the cardinality of $V^{*}$ is $k^E$, while that of $V$ is, if I'm not wrong, equal to that of $E$ (in this first step I am assuming $\operatorname{card} k \le \operatorname{card} E$).
Indeed $V$ is a union parametrized by $\mathbb{N}$ of sets of cardinality equal to $E$. In particular $\operatorname{card} V < \operatorname{card} V^{*}$, so the same inequality holds for the dimensions.
Second claim. Let $h \subset k$ two fields. If the thesis holds for vector spaces on $h$, then it holds for vector spaces on $k$.
Indeed let $V$ be a vector space over $k$, $E$ a basis. Functions with finite support from $E$ to $h$ form a vector space $W$ over $h$ such that $V$ is isomorphic to the extension of $W$, i.e. to $W\otimes_h k$. Every functional from $W$ to $h$ extends to a functional from $V$ to $k$, hence
$$\operatorname{dim}_k V = \operatorname{dim}_h W < \operatorname{dim}_h W^* \leq \operatorname{dim}_k V^*.$$
Putting the two claims together and using the fact that every field contains a field at most denumerable yields the thesis.
Let me give an alternative proof of the result in Ron Maimon's answer, that the canonical map $V\to V^{**}$ can consistently, in the absence of the axiom of choice, be surjective. The only technical advantage of my proof is that it uses only the consistency of ZF plus countable choice plus "all sets (in Polish spaces) have the Baire property". That consistency was proved (by Shelah) relative to just ZF, whereas the consistency of "all sets are Lebesgue measurable" needs an inaccessible cardinal. I think my proof is also a bit simpler than the one using measure. (Also, I don't need to mention ultrafilters, which some people might consider an advantage.)
Let $f:V^*\to\mathbb R$ be a linear map, and let me identify $V^*$ with the space of infinite sequences of reals. Topologize $V^*$ with the product topology, where each factor $\mathbb R$ has the usual topology of the reals. That makes $V^*$ a Polish space, so I can use the assumption about Baire category. In particular, if I partition $\mathbb R$ into intervals of length 1, then the inverse images of these intervals under $f$ have the Baire property, and they can't all be meager, by the Baire category theorem. So at least one of them, call it $f^{-1}(I)$, differs by a meager set from a nonempty open set. Inside that nonempty open set, I can find a basic open set of the form `B=$\prod_iU_i$ where, for some $n$, the first $n$ of the factors $U_i$ are intervals of some length $\delta$ and the later factors $U_i$ are $\mathbb R$. For the first $n$ indices $i$, let $U'_i$ be the interval with the same midpoint as $U_i$ but only half the length, and let $B'$ be the product that is like $B$ except using the $U'_i$ instead of the $U_i$ for the first $n$ factors. Consider an arbitrary $z\in V^*$ whose first $n$ components are smaller than $\delta/2$ in absolute value. Then translation by $z$ in $V^*$ maps $B'$ homeomorphically to another subset of $B$. So the two sets $f^{-1}(I)\cap B'$ and $\{x\in B':z+x\in f^{-1}(I)\}$ are both comeager in $B'$ and therefore must intersect. Let $x$ be in their intersection. Both $f(x)$ and $f(z+x)$ are in the interval $I$ of length 1. Subtracting (and remembering that $f$ is linear), we get that $|f(z)|\leq 2$.
Summarizing, we have that $|f(z)|\leq 2$ for all $z\in V^*$ whose first $n$ components are smaller than $\delta/2$ in absolute value. By linearity, if the first $n$ components of $z$ are smaller than $\alpha\delta/2$ in absolute value, for some positive $\alpha$, then $|f(z)|\leq2\alpha$. In particular, if the first $n$ components of $z$ are zero, then so is $f(z)$. That is, the kernel of $f$ includes the subspace $N$ consisting of those $z\in V^*$ whose first $n$ components vanish. So $f$ factors through the quotient $V^*/N$, which is finite-dimensional (in fact, $n$-dimensional). Knowing what linear functionals on finite-dimensional spaces look like, we immediately conclude that $f$ is given by inner product with a member of $V$ (having non-zero components in at most the first $n$ positions).
Best Answer
Search for Dirac structures or Courant algebroids in MathSciNet: These are common generalizations of symplectic and Poisson structures and use the symmetric bilinear form on $TM\times_M T^*M$ on a manifold: Namely, the graph of a symplectic structure as well as the graph of a Poisson structure are maximal isotropic subbundles, with further properties.
There is a lot of literature on them now.