Discretizing a Controller with the Backward Difference Method
This is kind of hand-wavy, but you can look at this from two different perspectives:

One, you can look at $z^{-1}$ as a "back-step" operator; i.e. if $X(z) = \mathcal{Z}\lbrace x_n \rbrace$, then (with a few 'i's left undotted and 't's uncrossed) $\frac{X(z)}{z} = \mathcal{Z}\lbrace x_{n-1} \rbrace$.  You can also look at $s$ as a derivative operator: if $X(s) = \mathcal{Z}\lbrace x(t) \rbrace$, then $s X(s) = \mathcal{Z}\lbrace \frac{d}{dt}x(t) \rbrace$.

And -- the time derivative of $x(t)$ is approximately $\frac{d}{dt}x(t) \simeq \frac{x(t) - x(t - T)}{T}$.  If I throw caution to the wind and just sample the right side, that translates to, again roughly, $\frac{d}{dt}x(t) \simeq \frac{x_n - x_{n-1}}{T}$.  Now the Laplace of the left and the Z transform of the right gives $s X \simeq \frac{1 - z^{-1}}{T} X$.  Note that I have left off all pretense of a frequency domain variable on this -- I think that if there's an antithesis of "rigor" it must be "limpor", and we've just found it.

So you can take a system transfer function in $s$, expand it into a differential equation, then you can use the approximation $\frac{d}{dt}x(t) \simeq \frac{x(t) - x(t - T)}{T}$, then you can take the Z transform, then you turn the resulting difference equation into a transfer function in $z$ -- and you have just taken the very long way around to simply substituting $\frac{1 - z^{-1}}{T}$ for $s$.

The other way you can look at this is to observe that with a continuous-time system, the delay operator has a Laplace-domain transform: $\mathcal{L}\lbrace x(t - T)\rbrace = e^{-sT} X(s)$.  So you just set $z^{-1} = e^{-sT}$.  There's ways you can make this more rigorous, but compared to my above treatment this is already granite.

For $sT \ll 1$, $e^{-sT} \simeq 1 - sT$.  Substitute in $z^{-1} = e^{-sT}$ and solve for $s$ and you get $s \simeq \frac{1 - z^{-1}}{T}$.  Voila.  To get the first forward difference start with $z = e^{sT} \simeq 1 + sT$ and you find that $s \simeq \frac{z - 1}{T}$.  Use $z = e^{sT} = \frac{e^{\frac{sT}{2}}}{e^{-\frac{sT}{2}}}$ and you get the Tustin approxmation.  

Note, when you actually go to do these, the word approximation.  You can generally get by with control systems because you need to sample way faster than the plant dynamics anyway, but -- approximation.