= Solution
Expand $X^a=x^a+\theta\psi^a$, $DX^a=\psi^a-i\theta\dot x^a$, and $g_{ab}(X)=g_{ab}(x)+\theta\psi^c\partial_cg_{ab}$. Performing the <Berezin integral> gives
$$
S=\int dt\left[\frac12g_{ab}\dot x^a\dot x^b
+\frac i2g_{ab}\dot\psi^a\psi^b
+\frac i2(\partial_cg_{ab})\psi^c\dot x^a\psi^b\right].
$$
After a fermionic integration by parts, this is the manifestly covariant expression
$$
\boxed{S=\frac12\int dt\left[g_{ab}\dot x^a\dot x^b-i g_{ab}\psi^a\frac{D\psi^b}{dt}\right]},
\qquad
\frac{D\psi^a}{dt}=\dot\psi^a+\Gamma^a{}_{bc}\dot x^b\psi^c.
$$
Under a <diffeomorphism>, $x$ is a point of $M$, $\psi\in T_xM$ is a tangent vector, $D_t\psi$ is its <covariant derivative along a curve>, and every index is contracted with the <Riemannian metric>; the action is therefore invariant.
Back to article page