The genetic characterization of Native American groups provides insights into their history and demographic events. We sequenced the mitochondrial D-loop region (control region) of 520 samples from eight Mexican indigenous groups. In addition to an analysis of the genetic diversity, structure and genetic relationship between 28 Native American populations, we applied Bayesian skyline methodology for a deeper insight into the history of Mesoamerica. AMOVA tests applying cultural, linguistic and geographic criteria were performed. MDS plots showed a central cluster of Oaxaca and Maya populations, whereas those from the North and West were located on the periphery. Demographic reconstruction indicates higher values of the effective number of breeding females (Nef) in Central Mesoamerica during the Preclassic period, whereas this pattern moves toward the Classic period for groups in the North and West. Conversely, Nef minimum values are distributed either in the Lithic period (i.e. founde)