Recall the following properties of the Fourier transform:
If $x(t)$ is an even function, then its Fourier transform $X(\omega)$ is purely real.
If $x(t)$ is an odd function, then its Fourier transform $X(\omega)$ is purely imaginary.
Thus, we can think of the real and imaginary parts of the Fourier transform of a zero-mean Gaussian random process $x(t)$ as the Fourier transforms of two separate inputs $x_e(t)$ and $x_o(t)$: the components of the process that have even and odd symmetry, respectively. We split the process into these components as follows:
$$
x_e(t) = \frac{x(t)+x(-t)}{2}
$$
$$
x_o(t) = \frac{x(t)-x(-t)}{2}
$$
Note that $x(t) = x_e(t) + x_o(t)$. Moving through this step by step,
Your observation was that $\text{Re}\{X(\omega)\}$ and $\text{Im}\{X(\omega)\}$ (the real and imaginary parts of the process's DFT) are uncorrelated.
Using the Fourier transform properties mentioned above, we can deduce that $\text{Re}\{X(\omega)\} = X_e(\omega)$ and $\text{Im}\{X(\omega)\} = X_o(\omega)$; the real and imaginary parts of $X(\omega)$ are none other than the Fourier transforms of $x_e(t)$ and $x_o(t)$, respectively.
Therefore, your observation is equivalent to saying that $X_e(\omega)$ and $X_o(\omega)$ are uncorrelated.
Since the Fourier transform is a one-to-one mapping between the time and frequency domains, I posit that the lack of correlation between $X_e(\omega)$ and $X_o(\omega)$ would imply a lack of correlation between $x_e(t)$ and $x_o(t)$ as well.
What is the correlation between $x_e(t)$ and $x_o(t)$? Simple:
$$
\begin{align}
\mathbb{E}(x_e(t)x_o(t)) &= \mathbb{E}\left(\left(\frac{x(t)+x(-t)}{2}\right)\left(\frac{x(t)-x(-t)}{2}\right)\right) \\
&= \mathbb{E}\left(\frac{1}{4}\left(x^2(t) - x^2(-t)\right)\right) \\
&= \frac{1}{4} \left(\mathbb{E}(x^2(t)) - \mathbb{E}(x^2(-t))\right) \\
&= \frac{1}{4} \left(\sigma^2 - \sigma^2\right) \\
&= 0
\end{align}
$$
As expected, the even and odd components are uncorrelated. So, to summarize, I would say the following:
The real and imaginary components of a Fourier transform correspond to the individual Fourier transforms even and odd components of the input function.
For a zero-mean Gaussian random process, these even and odd components are uncorrelated.
Therefore, their Fourier transforms (the real and imaginary components that you asked about) are also uncorrelated.
Edit: To address your followup:
If $x(t)$ is Gaussian, then its even and odd components $x_e(t)$ and $x_o(t)$ are as well, due to the property that any weighted sum of Gaussian random variables is also Gaussian.
If $x_e(t)$ and $x_o(t)$ are Gaussian random processes, then their Fourier transforms $X_e(\omega)$ and $X_o(\omega$) are as well. This follows from the same property as the previous statement; if you look at the transform, you're computing a weighted sum of a bunch of Gaussian random variables.
If $X_e(\omega)$ and $X_o(\omega)$ are Gaussian, and they are uncorrelated with one another (as described above), then they are also independent. This is a property of the Gaussian distribution.