Big Picture: What is the connection of Malliavin calculus with differential geometry? I know that Paul Malliavin was heavily influenced by ideas from differential geometry while developing his calculus on Wiener space. But what are the concrete analogies between both areas of mathematics? What has this to do with Hörmander's theorem (if so)?
 A: Some of Malliavin's ideas are well-explained in his own book: Stochastic Analysis. 
I think the main geometric idea behind his proof of Hormander's theorem is the idea of submersion. More precisely, if we consider a stochastic differential equation in Stratonovitch form
$dX_t=V_0(X_t) dt +\sum_{i=1}^n V_i(X_t) \circ dW^i_t
$
where $W$ is a Wiener process and the $V_i$'s vector fields on $R^n$ (or a manifold), we can see the solution $X_t$ as function of the path $\{ W_s , s \le t\}$. We thus have a map $\Pi$ from the Wiener space $\mathbb{W}$, which is the path space of $W$ and $\mathbb{R}^n$. 
A key insight of Malliavin is now to see the Wiener space $\mathbb{W}$ as an infinite-dimensional manifold with tangent space the Cameron-Martin space (space of functions with derivatives in $L^2$) and to realize that Hormander's conditions make the  map $\Pi$ a submersion. As a consequence the Wiener measure is mapped into a smooth measure, which means that if Hormander's conditions are satisfied then the distribution of $X_t$ has a smooth density. It implies the hypoellipticity of the generator of the Markov process $(X_t)_{t \ge 0}$ which is $L=V_0+\frac{1}{2} \sum_{i=1}^n V_i^2$.
Of course, this has to be made rigorous and some care has to be given because we have to use Sobolev type derivatives instead of Frechet derivatives.
This idea that there is a "natural" manifold structure on path spaces is quite fruitful and besides Malliavin's proof of Hormander's theorem explains and leads to several interesting properties of the operator $V_0+\frac{1}{2} \sum_{i=1}^n V_i^2$.
A: I can't speak for Paul Malliavin's influences, but I do know a bit about Hormander's theorem (by no means, an expert), and it is naturally suited to differentiable manifolds involving largely the idea of pullbacks of vectors. Malliavin calculus was apparently initiated to give a probabilistic proof of Hormander's theorem.
Following Rogers and Williams, Hormander's theorem concerns stochastic differential equations of the form
$$
\partial X = \sum_q U_q(X) \partial B^q + W(X)\partial t.
$$
Here, X is a stochastic process taking values in Rn, Uq and W are smooth vector fields, Bq are Brownian motions and  ∂ represents the Stratonovich integral.
As Stratonovich integration satisfies the standard change of variables formula, this SDE makes sense on an arbitrary differentiable manifold.
Next, according to the statement of Hormander's theorem, let [.,.] be the usual Lie Bracket for vector fields, and let A0,A1,... be the sequence of Lie algebras defined as follows.
$$
\begin{align}
& A_0={\rm Lie}(U_1,U_2,...),\\
& A_k={\rm Lie}([U,W]\colon U\in A_{k-1})
\end{align}
$$
Then Hormander's theorem states that if ∪nAn spans the tangent space at each point of Rn, then X has smooth transition densities.
Hormander's theorem is naturally a statement concerning diffusions on differentiable manifolds, as everything I said above makes perfect sense, and is true, if Rn is replaced by any differentiable manifold.
The idea behind the Malliavin proof is to consider differentiating with respect to perturbations of the Brownian motions Bq. The point is, that if Bq has a small bump applied at time t, then this creates a small bump in the solution X proportional to Uq at this time, which will then be propagated along the solution. In fact, solutions to the original SDE with all different points give rise to stochastically moving frames, and bumps in the solution X are transported along with these frames in a similar as vector fields give rise to transport of vectors along these fields.
The solution for X is smooth with respect to smooth bumps in the Brownian motions, which can be shown by converting bumps in the Brownian motion into changes in the probability measure, using Girsanov transforms. So, according to Malliavin calculus, you can always differentiate the solution with respect to the Brownian motions. The idea behind the proof of Hormanders theorem, is to invert this process of varying the Brownian motion→bumps in the final position of X. To do this, it is necessary to invert the process of transporting along the moving frames. That is, it must be a 1-1 map on the tangent spaces. Then, by a stochastic "pull-back" on the manifold (or Rn), you can interpret differentiating the solution with respect to its position at any time in terms of differentiating with respect to the Brownian motions.
So, this method of proving Hormander's theorem requires you to be able to differentiate with respect to the Brownian motion (i.e., Malliavin calculus) and the rest is all differential geometry (i.e., pullbacks of tangent vectors).
