On a compact Riemannian manifold $M$ (we assume Dirichlet boundary condition if $\partial M \neq \emptyset$), the Laplace-Beltrami operator $-\Delta$ has a discrete spectrum $0 < \lambda_1 \leq \lambda_2 < .....$ going to $\infty$. Consider the eigenspace corresponding to the first eigenvalue $\lambda_1$. My question is, what is the dimension of this eigenspace?

It seems that the answer should be one, and the otherwise one can locate two nodal domains corresponding to one of the first eigenfunctions, and then $\lambda_1$ is the first eigenvalue of each of these nodal domains, which should violate domain monotonicity.

On the other hand, suppose we are given one (even non-sign-changing) eigenfunction $\varphi$ corresponding to $\lambda$. Suppose the manifold $M$ is such that it has a non-trivial isometry $T$. Then $T\varphi (x) := \varphi (Tx)$ is also an eigenfunction, but in general one does not expect that $T\varphi = \varphi$.

What am I missing here?