Edit: Right now I will not give a complete argument in what follows (the previous version of the answer missed the part of the question asking for $S$ to be the boundary of a spacelike hypersurface in $M$), I will add the missing ingredients and/or additional hypotheses (if needed) as time allows in a future edit.
If by "future null cone" you mean the (achronal) boundary $\partial I^+(p)$ of the chronological future $I^+(p)$ of $p$, then the answerone can show that $S$ is a smooth, acausal hypersurface of (almost) yes$\partial M$ even without the nontrapping condition, thanks to your smoothness assumption. What may fail is that tangent vectors will be either spacelike or null, and the latter case cannot be excluded without additional assumptions on $\partial M$. For example, take $M$ as the Minkowski space-time minus an open cylinder around the time axis - $S$ will have a null tangent vector at any point where a null generator of $\partial I^+(p)$ meets $\partial M$ tangentially. Of course, this example has trapped null geodesics, but even with the nontrapping assumption the null generators could in principle meet $\partial M$ tangentially, causing $S$ to have null tangent vectors at the meeting point.
To rule out $S$ having null tangent vectors, a sufficient condition is that $\partial M$ is totally geodesic. More generally, one could impose a sort of "convexity" assumption on $\partial M$ with respect to null geodesics.
Edit: Let us assume from now on that $S$ is indeed spacelike. It remains to show the existence of a spacelike hypersurface $\Sigma\subset M$ such that $\partial\Sigma=S$. We will sketch below a partial argument assuming that $M_0$ is connected and $g$ is ultrastatic (i.e. $g_0(t,x)=g_0(x)$ does not depend on $t$). Let $p=(t_0,p_0)$ and let $\gamma_0:[0,1]\rightarrow M_0$ be a smooth curve connecting $p_0=\gamma_0(0)$ to $r_0=\gamma_0(1)$. Given $$\alpha>\sup_{s\in[0,1]}\sqrt{g_0(\gamma_0(s))(\dot{\gamma}_0(s),\dot{\gamma}_0(s))}\ ,$$ we have that the curve $$\gamma(s)=(t_0+\alpha s,\gamma_0(s))\ ,\quad s\in[0,1]$$ is timelike and links $p$ to the point $r=(t_0+\alpha,r_0)\in\partial M$. This shows that the timelike curve $t\mapsto (t,r_0)$ must cross $\partial I^+(p)$ at some $t\in(t_0,t_0+\alpha)$, which must be unique since $S$ is acausal. This shows that $\partial I^+(p)$ crosses the timelike curve $t\mapsto (t,q_0)$ exactly once for all $q_0\in\partial M_0$, thus implying that $S$ is the graph of a smooth function $\phi:\partial M_0\rightarrow\mathbb{R}$ - namely, if $q\in S$ then $q=(\phi(q_0),q_0)$ for a unique $q_0\in\partial M_0$. Since $S$ is spacelike, one must have $g_0^{-1}(q_0)(d\phi(q_0),d\phi(q_0))<1$. The missing piece is to extend $\phi$ to a smooth function $\Phi:M_0\rightarrow\mathbb{R}$ such that $g_0^{-1}(q_0)(d\Phi(q_0),d\Phi(q_0))<1$ - if that is accomplished, then the graph of $\Phi$ is the desired spacelike hypersurface $\Sigma$ with $\partial\Sigma=S$. It is straightforward to extend $\phi$ to a collar neighborhood of $\partial M_0$ in $M_0$ while keeping the derivative bound above. Using this to extend it further to a smooth function $\Phi$ on $M_0$ using a partition of unity is also straightforward - the problem is to do it in such a way so as to keep the derivative bound. I will complete the details in a future edit.
Connectedness of $M_0$ seems necessary to achieve the result, for otherwise it is hard to tell how to extend $S$ into the interior of $M$. As for ultrastaticity of $g$, it may be possible to weaken it to a form of global hyperbolicity suitable for Lorentzian manifolds with boundary.