The answer is no, even for locally constant families. Let $F$ be a smooth projective variety with an automorphism $\sigma$, and let $L$ be a line bundle on $F$ such that $\sigma^* L$ is not isomorphic to $L$. For example, take $X=\mathbf{P}^1\times \mathbf{P}^1$$F=\mathbf{P}^1\times \mathbf{P}^1$ with $\sigma$ the coordinate-switching involution, and $L=\mathcal{O}(1,0)$. Take $S$ to be the nodal curve obtained by gluing $0$ to $\infty$ on $\mathbf{P}^1$, and let $X$ be obtained from $F\times \mathbf{P}^1$ by identifying $(f, 0)$ and $(\sigma(f), 1)$. The natural $X\to S$ is a locally trivial fibration with fiber $F$. Take $s$ to be the node and $L_0$ to be $L$.
Think about it this way: in general, suppose that $X\to S$ is nice enough so that there exists a relative Picard scheme $P\to S$. The assumption that $H^2(X_s, \mathcal{O})=0$ for all $s\in S$ implies that $P\to S$ is (formally) smooth, and so every section over a point (i.e. your $\mathcal{L}_0$) extends to an infinitesimal/formal/etale neighborhood. But there is no reason to have a non-zero global section of $P\to S$ at all.