There is growing evidence that the Galactic Center Excess identified in the Fermi-LAT gamma-ray data arises from a population of faint astrophysical sources. We provide compelling supporting evidence by showing that the morphology of the excess traces the stellar over-density of the Galactic bulge. By adopting a template of the bulge stars obtained from a triaxial 3D fit to the diffuse near-infrared emission, we show that it is detected at high significance. The significance deteriorates when either the position or the orientation of the template is artificially shifted, supporting the correlation of the gamma-ray data with the Galactic bulge. In deriving these results, we have used more sophisticated templates at low-latitudes for the Fermi bubbles compared to previous work and the three-dimensional Inverse Compton (IC) maps recently released by the GALPROP team. Our results provide strong constraints on Millisecond Pulsar (MSP) formation scenarios proposed to explain the excess. We find that an admixture formation scenario, in which some of the relevant binaries are primordial and the rest are formed dynamically, is preferred over a primordial-only formation scenario at 7.6σ confidence level. Our detailed morphological analysis also disfavors models of the disrupted globular clusters scenario that predict a spherically symmetric distribution of MSPs in the Galactic bulge. For the first time, we report evidence of a high energy tail in the nuclear bulge spectrum that could be the result of IC emission from electrons and positrons injected by a population of MSPs and star formation activity from the same site. pulsars, in the form of either old "recycled" millisecond pulsars [13,18], or young pulsar remnants from Galactic Center supernovae [19,20]. The gamma-ray spectra of known pulsars are similar to the GCE, and the capture of pulsars onto the Galactic Center region by Globular Cluster disruptions could explain the GCE's quasi-spherical spatial distribution [21,22]. Other explanations for the GCE discussed in the literature include outbursts of cosmic-ray production by the central supermassive black hole (e.g., [23]).While pulsars are a natural astrophysical candidate, there is ongoing debate regarding the consistency of the luminosity function required to explain the GCE with that measured for the pulsar population observed elsewhere as point sources [24][25][26][27]. Nevertheless, there are strengthening indications that the GCE may have an astrophysical origin. Cosmic-ray interactions in the Galactic Center region are notoriously complex to model. New analyses have cast doubt on some of the main properties claimed previously for the GCE. First, the energy spectrum of the GCE is subject to large systematic uncertainties arising from the incomplete understanding of cosmic-ray interactions in the Galactic Center region (e.g., [13][14][15]17]). A variety of scenarios might, therefore, describe the GCE based on the energy spectrum alone. Second, it has been argued that the photon count distribution of the GCE is ...