If $A$ and $B$ are positive definite, then is $B^{-1} - A^{-1}$ positive semidefinite?
Yes. Here is one way to prove it:
Since $A-B\geq0$, we have (conjugating with $B^{-1/2}$) that $B^{-1/2}AB^{-1/2}-I\geq0$. This tells you that all the eigenvalues of the positive definite matrix $B^{-1/2}AB^{-1/2}$ are $\geq1$.
Now note that $B^{-1/2}AB^{-1/2}=(B^{-1/2}A^{1/2})(A^{1/2}B^{-1/2})$. Since commuting the product of two matrices does not change the eigenvalues, we get that $A^{1/2}B^{-1/2}B^{-1/2}A^{1/2}=A^{1/2}B^{-1}A^{1/2}$ has also all eigenvalues $\geq1$. So $$ A^{1/2}B^{-1}A^{1/2}-I\geq0. $$ Now conjugating with $A^{-1/2}$ we get $B^{-1}-A^{-1}\geq0$.
I like Martin Argerami's answer better, but here is a slightly different approach which does not use that the nonzero eigenvalues of $CD$ are the same as those of $DC$.
Note: this approach applies verbatim to the more general case of a $C^*$-algebra. Note also that in view of the notion of operator monotone function, we could interpret this result as saying that $x\longmapsto x^{-1}$ is operator nonincreasing.
The rather lengthy following paragraph shows that we can assume $B=I_n$ without loss of generality. Then it is obvious.
Since $B$ is positive definite, we have its square root $B^{1/2}$ which is also positive definite with positive definite inverse $B^{-1/2}$. Now since $A-B$ is positive semidefinite, so is $B^{-1/2}(A-B)B^{-1/2}=B^{-1/2}AB^{-1/2}-I_n=A_0-I_n$. And to prove $B^{-1}-A^{-1}$ positive semidefinite is equivalent to proving $B^{1/2}(B^{-1}-A^{-1})B^{1/2}=I_n-A_0^{-1}$ is positive semidefinite. Since $A_0=B^{-1/2}AB^{-1/2}$ is still positive definite, this shows that we can assume $B=I_n$ withtout loss of generality.
Now using the usual notation $C\geq D$ for symmetric matrices to mean that $C-D$ is positive semidefinite, or equivalently that the spectrum of $C-D$ is nonnegative, we get $$ A-I_n\geq 0\Rightarrow A^{-1/2}(A-I_n)A^{-1/2}=I_n-A^{-1}\geq 0. $$ So the result you mentioned holds, together, of course, with its converse.
Let me offer a constructive proof.
Lemma 1. For invertible and symmetric matrix $A$, and vectors $x$ and $y$, we have
$$ x^{'}A^{-1}x = \underset{y}{max\, \, } 2{x}'y - {y}'Ay $$
Lemma 2. For invertable matrices $A$ and $B$, we have
$$ (A + C) ^{-1} = A^{-1} - A^{-1} (A ^{-1} + C ^{-1}) ^{-1} A ^{-1} $$
Theorem. If $A > 0$, $B > 0$, $A$ and $B$ are symmetric matricex, and $A > B$, then $ B ^{-1} > A ^{-1} $.
Proof. By Lemma 1, we have
$$x ^{'} ( B ^{-1} - A ^{-1} ) x = \underset{y}{max\, \, } 2{x}'y - {y}' ( B ^{-1} - A ^{-1} ) ^{-1} y$$
Let $ y = A ^{-1} x $, we have $$ x ^{'} ( B ^{-1} - A ^{-1} ) x \ge 2x ^{'} A ^{-1} x - x ^{'} A ^{-1} ( B ^{-1} - A ^{-1} ) ^ {-1} A ^{-1} x $$
Then by Lemma 2, we have
$$ x ^{'} ( B ^{-1} - A ^{-1} ) x \ge x ^{'} A ^{-1} x + x ^{'} ( A - B ) ^{-1} x $$
By assumption, $ A > 0 $, so $ A ^{-1} > 0 $. Also $ A > B $, i.e. $ A - B > 0 $, so $ (A - B ) ^ {-1} > 0$. Thus the LHS of above equation is greater than 0, so $ B ^{-1} - A ^{-1} > 0 $, i.e. $ B ^{-1} > A ^{-1} $.
EDIT: typos.