Unraveling Signatures of Local Adaptation among Indigenous Groups from Mexico.

Few studies have addressed how selective pressures have shaped the genetic structure of the current Native American populations, and they have mostly limited their inferences to admixed Latin American populations. Here, we searched for local adaptation signals, based on integrated haplotype scores and population branch statistics, in 325 Mexican Indigenous individuals with at least 99% Native American ancestry from five previously defined geographical regions. Although each region exhibited its own local adaptation profile, only PPARG and AJAP1, both negative regulators of the Wnt/ß catenin signaling pathway, showed significant adaptation signals in all the tested regions. Several signals were found, mainly in the genes related to the metabolic processes and immune response. A pathway enrichment analysis revealed the overrepresentation of selected genes related to several biological phenotypes/conditions, such as the immune response and metabolic pathways, in agreement with previous studies, suggesting that immunological and metabolic pressures are major drivers of human adaptation. Genes related to the gut microbiome measurements were overrepresented in all the regions, highlighting the importance of studying how humans have coevolved with the microbial communities that colonize them. Our results provide a further explanation of the human evolutionary history in response to environmental pressures in this region.

DNA methylation and gene expression analysis in adipose tissue to identify new loci associated with T2D development in obesity.

BACKGROUND: Obesity is accompanied by excess adipose fat storage, which may lead to adipose dysfunction, insulin resistance, and type 2 diabetes (T2D). Currently, the tendency to develop T2D in obesity cannot be explained by genetic variation alone-epigenetic mechanisms, such as DNA methylation, might be involved. Here, we aimed to identify changes in DNA methylation and gene expression in visceral adipose tissue (VAT) that might underlie T2D susceptibility in patients with obesity. METHODS: We investigated DNA methylation and gene expression in VAT biopsies from 19 women with obesity, without (OND = 9) or with T2D (OD = 10). Differences in genome-scale methylation (differentially methylated CpGs [DMCs], false discovery rate < 0.05; and differentially methylated regions [DMRs], p value < 0.05) and gene expression (DEGs, p value <0.05) between groups were assessed. We searched for overlap between altered methylation and expression and the impact of altered DNA methylation on gene expression, using bootstrap Pearson correlation. The relationship of altered DNA methylation to T2D-related traits was also tested. RESULTS: We identified 11 120 DMCs and 96 DMRs distributed across all chromosomes, with the greatest density of epigenomic alterations at the MHC locus. These alterations were found in newly and previously T2D-related genes. Several of these findings were supported by validation and extended multi-ethnic analyses. Of 252 DEGs in the OD group, 68 genes contained DMCs (n = 88), of which 24 demonstrated a significant relationship between gene expression and methylation (p values <0.05). Of these, 16, including ATP11A, LPL and EHD2 also showed a significant correlation with fasting glucose and HbA1c levels. CONCLUSIONS: Our results revealed novel candidate genes related to T2D pathogenesis in obesity. These genes show perturbations in DNA methylation and expression profiles in patients with obesity and diabetes. Methylation profiles were able to discriminate OND from OD individuals; DNA methylation is thus a potential biomarker.

Kazak mitochondrial genomes provide insights into the human population history of Central Eurasia.

As a historical nomadic group in Central Asia, Kazaks have mainly inhabited the steppe zone from the Altay Mountains in the East to the Caspian Sea in the West. Fine scale characterization of the genetic profile and population structure of Kazaks would be invaluable for understanding their population history and modeling prehistoric human expansions across the Eurasian steppes. With this mind, we characterized the maternal lineages of 200 Kazaks from Jetisuu at mitochondrial genome level. Our results reveal that Jetisuu Kazaks have unique mtDNA haplotypes including those belonging to the basal branches of both West Eurasian (R0, H, HV) and East Eurasian (A, B, C, D) lineages. The great diversity observed in their maternal lineages may reflect pivotal geographic location of Kazaks in Eurasia and implies a complex history for this population. Comparative analyses of mitochondrial genomes of human populations in Central Eurasia reveal a common maternal genetic ancestry for Turko-Mongolian speakers and their expansion being responsible for the presence of East Eurasian maternal lineages in Central Eurasia. Our analyses further indicate maternal genetic affinity between the Sherpas from the Tibetan Plateau with the Turko-Mongolian speakers.

Ancestry-dependent genetic structure of the Xq28 risk haplotype in the Mexican population and its association with childhood-onset systemic lupus erythematosus.

Objective: Here we aimed to investigate the association of the Xq28 risk haplotype (H1) with susceptibility to childhood-onset systemic lupus erythematosus (SLE), and to compare its frequency and genetic structure in the Mexican population with those in other continental populations. Methods: We genotyped 15 single-nucleotide variants (SNVs) that form the H1 haplotype, using TaqMan real-time PCR. The association analysis [case-control and transmission disequilibrium test (TDT)] included 376 cases and 400 adult controls, all of whom were mestizos (MEZ). To identify risk alleles in Mexican Indigenous individuals, SNVs were imputed from whole-exome sequencing data of 1,074 individuals. The allelic frequencies determined in MEZ and Indigenous individuals were compared with those of the continental populations from the 1,000 Genomes database phase 3. Linkage disequilibrium (LD) analysis of risk alleles was performed on all populations. Interleukin-1 receptor associated kinase 1 (IRAK1) and methyl CpG binding protein 2 (MECP2) mRNA levels were determined using real-time PCR. Results: Case-control analysis revealed genetic association with childhood-onset SLE for all 15 SNVs (OR = 1.49-1.75; p = 0.0095 to 1.81 × 10-4) and for the Xq28 risk haplotype (OR = 1.97, p = 4 × 10-6). Comparing with individuals of European ancestry (0.14-0.16), the frequencies of the risk alleles were significantly higher in the MEZ individuals (0.55-0.68) and even higher in Indigenous individuals (0.57-0.83). LD analysis indicated a differential haplotype structure within the Indigenous groups, which was inherited to the MEZ population as a result of genetic admixture. Individuals homozygous for the Xq28 risk haplotype exhibited decreased levels of both MECP2A and B transcripts. Conclusion: We found that the H1 risk haplotype differs in its conformation in the Mexican population. This difference could be attributed to positive selection within the Indigenous population, with its inheritance now having an autoimmune health impact in both the Mexican Indigenous and MEZ populations.

The genomic landscape of Mexican Indigenous populations brings insights into the peopling of the Americas.

The genetic makeup of Indigenous populations inhabiting Mexico has been strongly influenced by geography and demographic history. Here, we perform a genome-wide analysis of 716 newly genotyped individuals from 60 of the 68 recognized ethnic groups in Mexico. We show that the genetic structure of these populations is strongly influenced by geography, and our demographic reconstructions suggest a decline in the population size of all tested populations in the last 15-30 generations. We find evidence that Aridoamerican and Mesoamerican populations diverged roughly 4-9.9 ka, around the time when sedentary farming started in Mesoamerica. Comparisons with ancient genomes indicate that the Upward Sun River 1 (USR1) individual is an outgroup to Mexican/South American Indigenous populations, whereas Anzick-1 was more closely related to Mesoamerican/South American populations than to those from Aridoamerica, showing an even more complex history of divergence than recognized so far.