In finite dimensions, if $T$ is a linear operator and $x$ is a (centered) Gaussian random variable, then $Tx$ is again a (centered) Gaussian random variable.
Now suppose that $x$ is a (say, centered) random variable on $\mathbb{R}^\mathbb{N}$.
Lemma 2.2.8 together with Example 2.3.5 in Bogachev's monograph on Gaussian measures seem to demonstrate that $Y =\sum_{j=1}^\infty f_j x_j$ is (distributed with) a centered Gaussian measure on the real line with variance $\|f\|_{\ell^2}^2$, when $f$ is in $\ell^2$.
This is interesting because even though $x$ lies in $\ell^2$ with probability 0, one can nonetheless define the random variable $\sum_{j=1}^\infty f_j x_j$---looking through the proof it seems to be because if we define $f^{(n)} \colon x \mapsto \sum_{j \leq n} f_j x_j$, then $\|f^{(n)} - f\|_{L^2(\gamma)}^2 = \sum_{j>n} f_j^2 \to 0$, demonstrating that $f$ is the limit of elements in $(\mathbb{R}^\mathbb{N})^\ast \cong c_{00}$, the space of sequences with finitely many nonzero. One can then use this to analyze the characteristic functional of $Y$.
Question: Now suppose that $T \colon \ell^2 \to \ell^2$ is a bounded linear operator and as above that $x$ is a random variable drawn according to the Gaussian law on $\mathbb{R}^\mathbb{N}$. What additional assumptions are needed so that $Z = Tx$ is a centered Gaussian random variable?