This is not really rigorous (the other answer covers that) and might be trivial, but I wanted to add this answer from a more intuitive perspective. For me, $\int f(W_s)\,dW_s$ always makes me do a double take, since $dW_s$ is such an odd creature, i.e. not an independent variable in the sense one is used to seeing.
Consider the Brownian increment $\Delta W_{t_{i+1}} = W_{t_{i+1}} - W_{t_i}$. By the definition of a Wiener process, if $\Delta t = t_{i+1}-t_i$, then (1) $\Delta W_{t_{i+1}} \sim \mathcal{N}(0,\Delta t) $ and (2) $W_{t_i}$ is independent from $\Delta W_{t_{i+1}}$.
Thus, expanded as a sort of Ito-Riemann sum:
\begin{align}
\mathbb{E}\left[ \int_0^t f(W_s)\;dW_s \right]
&\approx \mathbb{E}\left[ \sum_i f(W_{t_i})[W_{t_{i+1}} - W_{t_i}] \right] \\
&= \sum_i \mathbb{E}\left[ f(W_{t_i})\Delta W_{t_{i+1}} \right] \\
&= \sum_i \mathbb{E}\left[ f(W_{t_i})\right] \mathbb{E}\left[\Delta W_{t_{i+1}} \right]\;\;\;\;\;\;\;\;\;\;\;\text{ By (2)} \\
&= \sum_i \mathbb{E}\left[ f(W_{t_i})\right] 0 \;\;\;\;\;\;\;\;\;\;\;\;\,\;\;\;\;\;\;\;\;\;\;\;\text{ By (1)} \\
&=0
\end{align}