Ultraviolet Raman spectroscopy revealed the existence of an unusual large frequency shift occurring to a non-soft mode of E(TO4) when BaTiO3 is strained to a SrTiO3 substrate [D. Tenne et al., Science 313, 1614]. It raised two interesting questions: (i) whether there are other non-soft modes that possess similar or even larger strain-induced frequency shifts; (ii) how the mode sequence is altered by these shifts in frequency. Note that mode sequence is also pivotal in correctly indexing and assigning the spectroscopy peaks observed in all Raman experiments. By mapping out the evolutions of individual phonon modes as a function of strain using first-principles density functional perturbation calculations, we determine the mode sequence and strain-induced phonon frequency shifts in prototypical BaTiO3. Our study reveals that the mode sequence is drastically different when BaTiO3 is strained to SrTiO3 as compared to that in the unstrained structure, caused by multiple mode crossings. Furthermore, we predict that three other non-soft modes-A1(TO2), E(LO4), and A1(TO3)-display even larger strain-induced frequency shifts than E(TO4). The strain responses of individual modes are found to be highly mode specific, and a mechanism that regulates the magnitude of the frequency shift is provided. As another key outcome of this study, we tackle a long-standing problem of LO-TO splitting in ferroelectrics. A rigorous definition for the LO-TO splitting is formulated, which allows this critical quantity to be calculated quantitatively. The definition immediately reveals a new finding, that is, a large LO-TO splitting not only exists for E(LO4), which is previously known and originates from a soft mode, it also occurs to a non-soft A1(LO3) mode. The LO-TO splitting is shown to decrease drastically with compressive strain, and this decrease cannot be explained by the Born effective charges and high-frequency dielectric constants.