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.

To evaluate the associations between six single-nucleotide polymorphisms (SNPs) in intron 1 of FTO and body mass index (BMI), a case-control association study of 2314 unrelated Mexican-Mestizo adult subjects was performed. The association between each SNP and BMI was tested using logistic and linear regression adjusted for age, gender, and ancestry and assuming additive, recessive, and dominant effects of the minor allele. Association analysis after BMI stratification showed that all five FTO SNPs (rs1121980, rs17817449, rs3751812, rs9930506, and rs17817449), were significantly associated with obesity class II/III under an additive model (P<0.05). Interestingly, we also documented a genetic model-dependent influence of gender on the effect of FTO variants on increased BMI. Two SNPs were specifically associated in males under a dominant model, while the remainder were associated with females under additive and recessive models (P<0.05). The SNP rs9930506 showed the highest increased in obesity risk in females (odds ratio = 4.4). Linear regression using BMI as a continuous trait also revealed differential FTO SNP contributions. Homozygous individuals for the risk alleles of rs17817449, rs3751812, and rs9930506 were on average 2.18 kg/m2 heavier than homozygous for the wild-type alleles; rs1121980 and rs8044769 showed significant but less-strong effects on BMI (1.54 kg/m2 and 0.9 kg/m2, respectively). Remarkably, rs9930506 also exhibited positive interactions with age and BMI in a gender-dependent manner. Women carrying the minor allele of this variant have a significant increase in BMI by year (0.42 kg/m2, P = 1.17 x 10−10). Linear regression haplotype analysis under an additive model, confirmed that the TGTGC haplotype harboring all five minor alleles, increased the BMI of carriers by 2.36 kg/m2 (P = 1.15 x 10−5). Our data suggest that FTO SNPs make differential contributions to obesity risk and support the hypothesis that gender differences in the mechanisms involving these variants may contribute to disease development.

Methylenetetrahydrofolate reductase (MTHFR) is a key enzyme in folate metabolism. Folate deficiency has been related to several conditions, including neural tube defects (NTDs) and cardiovascular diseases. Hence, MTHFR genetic variants have been studied worldwide, particularly the C677T and A1298C. We genotyped the C677T and A1298C MTHFR polymorphisms in Mexican Amerindians (MAs), from the largest sample included in a genetic study (n = 2026, from 62 ethnic groups), and in a geographically-matched Mexican Mestizo population (MEZ, n = 638). The 677T allele was most frequent in Mexican individuals, particularly in MAs. The frequency of this allele in both MAs and MEZs was clearly enriched in the South region of the country, followed by the Central East and South East regions. In contrast, the frequency of the 1298C risk allele in Mexicans was one of the lowest in the world. Both in MAs and MEZs the variants 677T and 1298C displayed opposite allele frequency gradients from southern to northern Mexico. Our findings suggest that in Mestizos the 677T allele was derived from Amerindians while the 1298C allele was a European contribution. Some subgroups showed an allele frequency distribution that highlighted their genetic diversity. Notably, the distribution of the frequency of the 677T allele was consistent with that of the high incidence of NTDs reported in MEZ.

Extensive research has shown that aberrant expression of microRNAs (miRNAs) plays an important role in innate and adaptive immune responses. The rs2910164 polymorphism has been identified as a functional variant, which affects the transcription and expression level of miR-146a and, thereby, contributes to the pathogenesis of several inflammatory and autoimmune diseases. To investigate whether the rs2910164 G/C polymorphism was associated with asthma, systemic lupus erythematosus (SLE) or juvenile rheumatoid arthritis (JRA), we performed an association study in a pediatric Mexican cohort. We included 979 pediatric patients (asthma: 402, SLE: 367 and JRA: 210) and 531 control subjects without inflammatory or immune diseases. Genotyping was performed using the 5' exonuclease technique. The genotype distribution of the rs2910164 polymorphism was in Hardy-Weinberg equilibrium in each group. No significant differences were detected in the distribution of this polymorphism between cases and controls (P = 0.108, 0.609 and 0.553 for subjects with asthma, JRA and SLE, respectively). However, stratification by gender showed a statistically significant difference between asthmatic and control females, where the C allele was significantly associated with protection to asthma (odds ratio = 0.694, 95% confidence interval 0.519-0.929, P = 0.0138). Our results provide evidence that rs2910164 may play a role in the susceptibility to childhood-onset asthma, but not SLE or JRA in Mexicans. Further association studies may contribute to determining the role of miR-146a single-nucleotide polymorphisms in immune-mediated diseases.