Polaron spectral functions are computed for highly doped graphene-on-substrate and other atomically thin graphitic systems using the diagrammatic Monte Carlo technique. The specific aim is to investigate the effects of interaction on spectral functions when the symmetry between sub-lattices of a honeycomb lattice has been broken by the substrate or ionicity, inducing a band gap. Introduction of electron-phonon coupling leads to several polaronic features, such as band-flattening and changes in particle lifetimes. At the K point, differences between energies on each sub-lattice increase with electron-phonon coupling, indicating an augmented transport gap, while the spectral gap decreases slightly. Effects of phonon dispersion and long-range interactions are investigated, and found to lead to only quantitative changes in spectra.