In this study, we compared the simulated absorption and fluorescence lineshapes computed with three combined ensemble-Franck-Condon methods that account for both explicit environmental effects and vibronic effects. The explicit solvent is modeled through QM/MM AIMD trajectories performed on the ground state and the excited state. The vibronic effects are captured through the simulation of Franck-Condon spectra for the chromophore embedded in a frozen QM solvent shell. These E-FC methods are used to generate spectral lineshapes for three different chromophore-solvent systems and are compared to the more standard ensemble and implicit solvent Franck-Condon approaches.