Semiflexible polymers undergo a weakly first order isotropic-to-nematic (IN) phase transition when the volume fraction φ is high enough that random alignment of the backbone segments is no longer viable. For semiflexible chains, the critical volume fraction φc is governed by the backbone stiffness Np. To locate the IN phase transition, we perform molecular dynamics (MD) simulations of bead-spring chains confined between two impenetrable parallel surfaces. We use the impenetrable surfaces to induce nematic-isotropic interfaces for semiflexible chains in the isotropic phase. By progressively increasing the backbone stiffness Np, we observe the propagation of surface-induced nematic order above a critical stiffness Ncp for a given φ. Using the simulation results Ncp(φ), we construct the IN phase boundry in the φ-Np plane, from which the scaling relation between φc and Np is obtained. For semiflexible chains with Np ≤ 5.78, our results suggest φc ∼ Np-1, consistent with prediction by Khokhlov and Semenov. For chains with Np ≥ 5.78, we observe a new scaling regime in which φc ∼ Np-2/3.