Consider exercise 34 from chapter 4 ("Hilbert Spaces: An Introduction") of [1] (p. 201):
Let K be a Hilbert-Schmidt kernel which is real and symmetric. Then, as we saw, the operator $T$ whose kernel is $K$ is compact and symmetric. Let $\{\varphi_k(x)\}$ be the eigenvectors (with eigenvalues $\lambda_k$) that diagonalize $T$. Then:
(a) $\sum_k |\lambda_k|^2 < \infty$
(b) $K(x,y) \sim \sum\lambda_k\varphi_k(x)\varphi_k(y)$ is the expansion of $K$ in the basis $\{\varphi_k(x)\varphi_k(y)\}$.
(c) Suppose $T$ is a compact operator which is symmetric. Then $T$ is of Hilbert-Schmidt type if and only if $\sum_n |\lambda_n|^2 < \infty$, where $\{\lambda_n\}$ are the eigenvalues of $T$ counted according to their multiplicities.
I'd appreciate help in solving parts (b) and (c) (I have figured out part (a) on my own). For completeness, here's how [1] defines a Hilbert-Schmidt kernel (pp. 186-187). Note that $\mathcal{H}$ stands for a Hilbert space and that only separable Hilbert spaces are considered throughout [1].
Let $\mathcal{H} = L^2(\mathbb{R}^d)$. If we can define an operator $T:\mathcal{H}\rightarrow\mathcal{H}$ by the formula $$ T(f)(x) = \int_{\mathbb{R}^d} K(x,y)f(y)\ dy,\hspace{1cm} \mbox{whenever } f \in L^2(\mathbb{R}^d), $$ we say that the operator $T$ is an integral operator and $K$ is its associated kernel. [...] Hilbert-Schmidt operators [... are integral operators ...] with a kernel $K$ that belongs to $L^2(\mathbb{R}^d\times \mathbb{R}^d)$.
References
[1] Stein, Elias M. and Shakarchi, Rami. Real Analysis: Measure Theory, Integration, and Hilbert Spaces. Princeton University Press (2005)
Let $H$ be a separable Hilbert space as desired. For simplicity assuming we are working over the real field. Then being an integral operator, we have $$ K: H\rightarrow H, (Kf)(x)=\int K(x,y)f(y)dy $$ Since $K$ is a symmetric compact operator, it is normal and can be diagonalized. Let us write its eigenvectors as $\phi_{j}$ with eigenvalue $\lambda_{j}$. Then we know $\phi_{j}$s span $H$. Thus $K$ can be defined alternatively by $$ \forall f=\sum a_{j}\phi_{j}, (K(\sum a_{j}\phi_{j}))(x)=\sum a_j \lambda_{j}\phi_{j}(x)=\int_{H} \sum\lambda_{j} \phi_{j}(x)\phi_{j}(y)(\sum a_{i}\phi_{i}(y))dy $$ where we normalized $\langle \phi_{i}, \phi_{j}\rangle=\delta_{ij}$. The equality holds because for each term we have $$ a_{i}\lambda_{j}\phi_{j}(x)\int_{H}\phi_{j}(y)\phi_{i}(y)dy=a_{i}\lambda_{j}\phi_{j}(x)\delta_{ij} $$ and this concludes the proof.