The Weierstrass $\wp$ function satisfies the addition formula
$$\wp(z+Z)+\wp(z)+\wp(Z) = \left(\frac{\wp'(z)-\wp'(Z)}{\wp(z)-\wp(Z)}\right)^2.$$
Of course, this is just the $x$-coordinate of the sum of the points $(\wp'(z), \wp(z))$ and $(\wp'(Z), \wp(Z))$ on the Weierstrass elliptic curve $y^2=4x^3-g_2x-g_3$. If one has an a priori knowledge of this fact, the computation of the addition formula is absolutely trivial. However, it seems that, in the literature, there is no "illuminating" proof from first principles. Perhaps my standards for "illuminating" are too high, but in comparison with proofs of addition formulas for trigonometric functions, the addition formula for $\wp$ is painstaking to establish.
I would like to know how the formula above, or an equivalent formula (perhaps for the Weierstrass $\sigma$-functions?) may be deduced in the most direct possible manner. Of course, a proof such as one that involves the comparison of power series around the origin is very direct, but it requires a first-hand knowledge of the formula, so it's not quite what I'm looking for.
All the best, and thank you!
I'm not quite sure if this is the sort of proof that'll satisfy, but I'll write out (a sketch of) one proof in Akhiezer's book.
The prerequisite is the knowledge that an arbitrary elliptic function $f(u)$ can be represented in terms of Weierstrass $\sigma$:
$$f(u)=k\frac{\prod\limits_{j=1}^n\sigma(u-p_i;g_2,g_3)}{\prod\limits_{j=1}^n\sigma(u-q_i;g_2,g_3)}$$
where the $p_i$ are zeroes and $q_i$ are poles within the period parallelogram (multiple poles/zeroes being indicated according to their order), and these satisfy the equation $\sum_i p_i=\sum_i q_i$.
Consider the function $\wp(u)-\wp(v)$ where $v$ is a constant and $\wp(v)$ is finite. (I do not indicate the invariants $g_2,g_3$ here and in the sequel due to my laziness. ;)) The function has a double pole at $u=0$ and two zeroes at $u=v$ and $u=-v$ (hint: $\wp$ is even), so we have the representation
$$\wp(u)-\wp(v)=k\frac{\sigma(u-v)\sigma(u+v)}{\sigma^2(u)}$$
$k$ can be determined by multiplying both sides of the equation by $u^2$ and letting $u\to0$, which gives $k=-\frac1{\sigma^2(v)}$; thus ending up with
$$\wp(u)-\wp(v)=-\frac{\sigma(u-v)\sigma(u+v)}{\sigma^2(u)\sigma^2(v)}$$
If we logarithmically differentiate both sides, we have
$$\frac{\wp^\prime(u)}{\wp(u)-\wp(v)}=\zeta(u+v)-2\zeta(u)+\zeta(u-v)$$
where $\zeta(u)$ is Weierstrass $\zeta$. If $u$ and $v$ are swapped and we remember that $\zeta$ is odd, we have
$$-\frac{\wp^\prime(v)}{\wp(u)-\wp(v)}=\zeta(u+v)-2\zeta(v)-\zeta(u-v)$$
Adding up these two and rearranging yields the (pseudo-)addition formula for $\zeta$:
$$\zeta(u+v)=\zeta(u)+\zeta(v)+\frac12\frac{\wp^\prime(u)-\wp^\prime(v)}{\wp(u)-\wp(v)}$$
Differentiating the (pseudo-)addition formula with respect to $u$ (and remembering that $\zeta^\prime=-\wp$) gives
$$-\wp(u+v)=-\wp(u)+\frac12\frac{\wp^{\prime\prime}(u)(\wp(u)-\wp(v))-\wp^\prime(u)(\wp^\prime(u)-\wp^\prime(v))}{(\wp(u)-\wp(v))^2}$$
Adding this up with the derivative with respect to $v$, and remembering that $\wp^{\prime\prime}=6\wp^2-\frac{g_2}{2}$, we can simplify the mess and thus obtain the addition formula for $\wp$.
$$\wp(u+v)+\wp(u)+\wp(v)=\frac14\left(\frac{\wp'(u)-\wp'(v)}{\wp(u)-\wp(v)}\right)^2$$
(i.e. there's a typo in the OP... ;))