What follows is not rigorous, but hopefully has the main idea. First, you probably want to justify
$$X_t-X_s \sim N\left(0, \int_s^t f(u)^2 du \right), $$
which can be done by approximating $f$ by a simple function $g = a_1 1_{(s,t_1)} + a_21_{(t_1,t_2)} + \ldots + a_n1_{(t_{n-1},t)}$ and then using
$$\int_s^t g(u)dW_u = a_1(W_{t_1} - W_s) + a_2 (W_{t_2}-W_{t_1}) + \ldots a_n (W_{t}-W_{t_{n-1}}) \sim N(0, a_1^2t_1 + a_2^2 (t_2-t_1)+\ldots a_n^2(t_n-t_{n-1})) = N\left(0, \int_s^t g(u)^2 du \right)$$
for $W_{t_i}-W_{t_{i-1}} \sim N(0,t_i-t_{i-1})$ all independent. The trick then is to show that $X_s$ and $X_t-X_s$ are uncorrelated in order to help us take the expectation $\mathbb{E}[\exp(a_1X_s + a_2 (X_t-X_s))]$. Now,
$$\mathbb{E}[(X_t-X_s)^2] = \mathbb{E}[X_t^2] - 2\mathbb{E}[X_tX_s] + \mathbb{E}[X_s^2],$$
hence,
$$\mathbb{E}[X_sX_t] = \frac{1}{2}\left(\mathbb{E}[X_t^2] + \mathbb{E}[X_s^2] - \mathbb{E}[(X_t-X_s)^2] \right) = \frac{1}{2} \left( \int_0^t f(u)du + \int_0^s f(u)du - \int_s^t f(u) du \right) = \int_0^sf(u)du = \mathbb{E}[X_s^2]$$
and so $\mathbb{E}[X_s(X_t-X_s)]=0,$ hence $X_s$ and $X_t-X_s$ are uncorrelated.
Now, since two uncorrelated normally distributed random variables $Y_1 \sim N(0,\sigma_1^2)$ and $Y_2 \sim N(0,\sigma_2^2)$ satisfy $\mathbb{E}[\exp(a_1Y_1 + a_2 Y_2)] = \exp(\frac{a_1}{2} Y_1 + \frac{a_2}{2} Y_2) = \mathbb{E}[\exp(a_1Y_1)]\mathbb{E}[\exp(a_2Y_2)]$, we have
$$\mathbb{E}[a_1X_s + a_2(X_t-X_s)] = \mathbb{E}[\exp(a_1X_s)] \mathbb{E}[\exp(a_2(X_t-X_s))]$$
for all $a_1$, $a_2$. Taking $n$ derivatives with respect to $a_1$ and $m$ derivatives with respect to $a_2$, and then setting $a_1=0=a_2$, we get
$$\mathbb{E}[X_s^n(X_t-X_s)^m] = \mathbb{E}[X_s^n]\mathbb{E}[(X_t-X_s)^m].$$
Using this, we now know that for any polynomial functions $p$ and $q$, we have $\mathbb{E}[p(X_s)q(X_t-X_s)] = \mathbb{E}[p(X_s)]\mathbb{E}[q(X_t-X_s)].$ Since polynomial functions are weakly dense, the probability distribution function $\rho_{X_s,X_t-X_s}$ of $X_s$ and $X_t-X_s$ is the product of the probability density functions $\rho_{X_s}$ of $X_s$ and $\rho_{X_t-X_s}$ of $X_t-X_s$. Hence, $X_s$ and $X_t-X_s$ are independent.