The time integration of Brownian motion is given by $X(t)=\int_0^tB(s)ds$. In order to show that $X(t)$ is a gaussian process I work with the definition that for any $t_0<t_1 < t_2<\ldots <t_{n-1}$ the random variable
$$ \langle v,(X_{t_0}, ... , X_{t_{n-1}}) \rangle $$
is again Gaussian.
I tried to write out the quantity in terms of Riemann sums i.e.
$$ X(t) = \lim_{n\to {\infty}}\sum_{k=0}^{n-1}B(t_k)(t_{k+1}-t_k) $$ But got stuck with all the indices, how to continue?
We can assume that limits of gaussian random variables are again Gaussian.
Using integration by parts we get: $$X(t) = t B(t) - \int_0^ts dB(s)$$ Recall that the Itô integral of a deterministic function is a Gaussian process; thus $X(t)$ is the linear combination of two Gaussian process, making it Gaussian itself.
Edit: We can prove this without directly appealing to results from Itô calculus (as requested in the comments). Using summation by parts, we get: $$\sum_{k=0}^{n-1}B(t_{k+1})(t_{k+1} - t_k) = B(t_n)t_n - \sum_{k=1}^{n-1} t_k(B(t_{k+1})- B(t_{k}))$$ The right hand side can be seen to be a Gaussian process; sending $n \to \infty$ we get that $X(t)$ is a Gaussian process.