What constraints are imposed on a wavefunction by the symmetry of the system?
I use quantum chemistry about as much as a lumberjack uses math (that is, just enough not to fell a tree on one's own head), so my answer will be short and hopefully useful as a tl;dr to someone else's more comprehensive take.

In quantum chemistry, if two operators commute, their values are simultaneously observable. In other words, there is a basis of functions which are eigenfunctions to both operators at once. For an opposite example, think of coordinate and momentum; they do not commute, hence the Heisenberg uncertainty and stuff.

When our problem has some symmetry (including, but not limited to, geometrical symmetry like $T_d$), this means that symmetry operators commute with its Hamiltonian. Therefore the eigenfunctions of the latter must also be eigenfunctions of the former. In other words, the application of symmetry operator to any of our states is equivalent to a multiplication by some number. That's the constraint you were after. Since a symmetry operator of a finite point group must have finite order (by applying it a few times, you get back where you started), the said number must be a root of unity.

That's about the size of it. The numbers for different symmetry operators can be chosen in a multitude of ways, each corresponding to some representation. The character tables tell you how many functions of each representation do you have. Since the functions of different representations are orthogonal to each other and don't interact, this helps us to split a huge problem into smaller chunks.



With methane, the qualitative result could have been envisioned in a simpler way, almost without any calculations. Think of a carbon atom with four $sp^3$ orbitals. They are all equivalent and orthogonal to each other. Now bring in the four H atoms. This will produce four $\sigma$ orbitals (and $\sigma^*$ too, we just don't care about those at the moment). But wait, those $\sigma$ orbitals are no longer orthogonal, so they must interact, if only a little. And because of symmetry, each two of them overlap by exactly the same amount. So it is as simple as Hückel can be: find the eigenvalues of the matrix $$\begin{pmatrix}
 0&amp;-1&amp;-1&amp;-1\\
-1&amp; 0&amp;-1&amp;-1\\
-1&amp;-1&amp; 0&amp;-1\\
-1&amp;-1&amp;-1&amp; 0
\end{pmatrix}$$
The symmetry approach I outlined above still applies; it is just that the system became so small that one can manage without it. The thing can be solved by hand. Or it can be noted that we're dealing with a complete graph, and their spectra are well known (just remember that math people have it upside down, for they use $\mathbf1$s instead of $\mathbf{-1}$s). Either way, the eigenvalues are $\{-3,1,1,1\}$, which kinda explains the 1:3 thing.

So it goes.