We find a large solution to a semilinear Dirichlet problem in a bounded $C^{1,1}$ domain for a non-local operator $φ(-Δ\vert_{D})$, an extension of the infinitesimal generator of a subordinate killed Brownian motion. The setting covers and extends the case of the spectral fractional Laplacian. The upper bound for the explosion rate of the large solution is obtained, and is given in terms of the renewal function, distance to the boundary, and the Keller-Osserman-type transformation of the nonlinearity. Additionally, we prove interior higher regularity results for this operator.