Life insurance contracts typically possess various embedded options. We focus on common options with early exercise features such as paid-up options, resumption options, surrender options and combinations among them. Contrary to the existing literature, showing these options value little under a deterministic-term structure, we demonstrate the situation changes dramatically whenever stochastic interest rates are introduced. With two stochastic sources (asset and interest rates), deriving optimal stopping strategies for multiple options becomes complex. Therefore, in this paper, an extension of the least squares Monte Carlo method (LSMC) is developed to allow valuation of these multi-feature options.