Mexico harbors great cultural and ethnic diversity, yet fine-scale patterns of human genome-wide variation from this region remain largely uncharacterized. We studied genomic variation within Mexico from over 1,000 individuals representing 20 indigenous and 11 mestizo populations. We found striking genetic stratification among indigenous populations within Mexico at varying degrees of geographic isolation. Some groups were as differentiated as Europeans are from East Asians. Pre-Columbian genetic substructure is recapitulated in the indigenous ancestry of admixed mestizo individuals across the country. Furthermore, two independently phenotyped cohorts of Mexicans and Mexican Americans showed a significant association between sub-continental ancestry and lung function. Thus, accounting for fine-scale ancestry patterns is critical for medical and population genetic studies within Mexico, in Mexican-descent populations, and likely in many other populations worldwide.

Performing genetic studies in multiple human populations can identify disease risk alleles that are common in one population but rare in others1, with the potential to illuminate pathophysiology, health disparities, and the population genetic origins of disease alleles. We analyzed 9.2 million single nucleotide polymorphisms (SNPs) in each of 8,214 Mexicans and Latin Americans: 3,848 with type 2 diabetes (T2D) and 4,366 non-diabetic controls. In addition to replicating previous findings2–4, we identified a novel locus associated with T2D at genome-wide significance spanning the solute carriers SLC16A11 and SLC16A13 (P=3.9×10−13; odds ratio (OR)=1.29). The association was stronger in younger, leaner people with T2D, and replicated in independent samples (P=1.1×10−4; OR=1.20). The risk haplotype carries four amino acid substitutions, all in SLC16A11; it is present at ≈50% frequency in Native American samples and ≈10% in East Asian, but rare in European and African samples. Analysis of an archaic genome sequence indicated the risk haplotype introgressed into modern humans via admixture with Neandertals. The SLC16A11 mRNA is expressed in liver, and V5-tagged SLC16A11 protein localizes to the endoplasmic reticulum. Expression of SLC16A11 in heterologous cells alters lipid metabolism, most notably causing an increase in intracellular triacylglycerol levels. Despite T2D having been well studied by genome-wide association studies (GWAS) in other populations, analysis in Mexican and Latin American individuals identified SLC16A11 as a novel candidate gene for T2D with a possible role in triacylglycerol metabolism.

Common polymorphisms in the fat mass and obesity-associated gene (FTO) have shown strong association with obesity in several populations. In the present study, we explored the association of FTO gene polymorphisms with obesity and other biochemical parameters in the Mexican population. We also assessed FTO gene expression levels in adipose tissue of obese and nonobese individuals. The study comprised 788 unrelated Mexican-Mestizo individuals and 31 subcutaneous fat tissue biopsies from lean and obese women. FTO single-nucleotide polymorphisms (SNPs) rs9939609, rs1421085, and rs17817449 were associated with obesity, particularly with class III obesity, under both additive and dominant models (P = 0.0000004 and 0.000008, respectively). These associations remained significant after adjusting for admixture (P = 0.000003 and 0.00009, respectively). Moreover, risk alleles showed a nominal association with lower insulin levels and homeostasis model assessment of B-cell function (HOMA-B), and with higher homeostasis model assessment of insulin sensitivity (HOMA-S) only in nonobese individuals (P dom = 0.031, 0.023, and 0.049, respectively). FTO mRNA levels were significantly higher in subcutaneous fat tissue of class III obese individuals than in lean individuals (P = 0.043). Risk alleles were significantly associated with higher FTO expression in the class III obesity group (P = 0.047). In conclusion, FTO is a major risk factor for obesity (particularly class III) in the MexicanMestizo population, and is upregulated in subcutaneous fat tissue of obese individuals.

It has been suggested that the higher susceptibility of Hispanics to metabolic disease is related to their Native American heritage. A frequent cholesterol transporter ABCA1 (ATP-binding cassette transporter A1) gene variant (R230C, rs9282541) apparently exclusive to Native American individuals was associated with low high-density lipoprotein cholesterol (HDL-C) levels, obesity and type 2 diabetes in Mexican Mestizos. We performed a more extensive analysis of this variant in 4405 Native Americans and 863 individuals from other ethnic groups to investigate genetic evidence of positive selection, to assess its functional effect in vitro and to explore associations with HDL-C levels and other metabolic traits. The C230 allele was found in 29 of 36 Native American groups, but not in European, Asian or African individuals. C230 was observed on a single haplotype, and C230-bearing chromosomes showed longer relative haplotype extension compared with other haplotypes in the Americas. Additionally, single-nucleotide polymorphism data from the Human Genome Diversity Panel Native American populations were enriched in significant integrated haplotype score values in the region upstream of the ABCA1 gene. Cells expressing the C230 allele showed a 27% cholesterol efflux reduction (P< 0.001), confirming this variant has a functional effect in vitro. Moreover, the C230 allele was associated with lower HDL-C levels (P = 1.77 x 10(-11)) and with higher body mass index (P = 0.0001) in the combined analysis of Native American populations. This is the first report of a common functional variant exclusive to Native American and descent populations, which is a major determinant of HDL-C levels and may have contributed to the adaptive evolution of Native American populations.

Although ATP-binding cassette transporter A1 (ABCA1) is well known for its role in cholesterol efflux and HDL formation, it is expressed in various tissues, where it may have different functions. Because hypoalphalipoproteinemia is highly prevalent in Mexico, we screened the ABCA1 coding sequence in Mexican individuals with low and high HDL cholesterol levels to seek functional variants. A highly frequent nonsynonymous variant (R230C) was identified in low-HDL cholesterol but not in high-HDL cholesterol individuals (P ‫؍‬ 0.00006). We thus assessed its frequency in the Mexican-Mestizo general population, seeking possible associations with several metabolic traits. R230C was screened in 429 Mexican Mestizos using Taqman assays, and it was found in 20.1% of these individuals. The variant was significantly associated not only with decreased HDL cholesterol and apolipoprotein A-I levels but also with obesity (odds ratio 2.527, P ‫؍‬ 0.005), the metabolic syndrome (1.893, P ‫؍‬ 0.0007), and type 2 diabetes (4.527, P ‫؍‬ 0.003). All of these associations remained significant after adjusting for admixture (P ‫؍‬ 0.011, P ‫؍‬ 0.001, and P ‫؍‬ 0.006, respectively). This is the first study reporting the association of an ABCA1 variant with obesity and obesityrelated comorbidities as being epidemiologically relevant in the Mexican population.