Consuelo D. Quinto-Cortés

Type 2 diabetes (T2D) is a highly heritable, polygenic disease with over 600 loci identified through genome-wide association studies (GWAS). However, despite possessing unique genetic variation shaped by demographic history and admixture, Latin American populations remain markedly underrepresented in global genomic research. To address this gap, we conducted genome- and exome-wide analyses of 19,431 T2D cases and 105,611 controls from the Mexico City Prospective Study (MCPS). We identified 86 independent GWAS associations, including 21 novel signals, 15 of which replicated in external cohorts. Risk alleles at novel loci were enriched in individuals with Indigenous American ancestry. Exome analyses revealed rare and ultra-rare missense variants with substantial risk effects at HNF1A and GCK , as well as a protein-damaging variant in SLC30A8 that reduced T2D risk by 45% in carriers. Integrative analyses indicate that T2D genetic architecture in Mexico is predominantly driven by common regulatory variation acting in the endocrine pancreas. Polygenic risk scores strongly stratified T2D risk and transferred to Indigenous Mexican populations. These findings demonstrate the power of large-scale genetic discovery in diverse populations to refine disease architecture and identify loci with potential therapeutic relevance.

Papua New Guinea (PNG) hosts distinct environments mainly represented by the ecoregions of the Highlands and Lowlands that display increased altitude and a predominance of pathogens, respectively. Since its initial peopling approximately 50,000 years ago, inhabitants of these ecoregions might have differentially adapted to the environmental pressures exerted by each of them. However, the genetic basis of adaptation in populations from these areas remains understudied. Here, we investigated signals of positive selection in 62 highlanders and 43 lowlanders across 14 locations in the main island of PNG using whole-genome genotype data from the Oceanian Genome Variation Project (OGVP) and searched for signals of positive selection through population differentiation and haplotype-based selection scans. Additionally, we performed archaic ancestry estimation to detect selection signals in highlanders within introgressed regions of the genome. Among highland populations we identified candidate genes representing known biomarkers for mountain sickness (SAA4, SAA1, PRDX1, LDHA) as well as candidate genes of the Notch signaling pathway (PSEN1, NUMB, RBPJ, MAML3), a novel proposed pathway for high altitude adaptation in multiple organisms. We also identified candidate genes involved in oxidative stress, inflammation, and angiogenesis, processes inducible by hypoxia, as well as in components of the eye lens and the immune response. In contrast, candidate genes in the lowlands are mainly related to the immune response (HLA-DQB1, HLA-DQA2, TAAR6, TAAR9, TAAR8, RNASE4, RNASE6, ANG). Moreover, we find two candidate regions to be also enriched with archaic introgressed segments, suggesting that archaic admixture has played a role in the local adaptation of PNG populations.

Background: Polygenic risk scores (PRS) derived from European individuals have reduced portability across global populations, limiting their clinical implementation at worldwide scale. Here, we investigate the performance of a wide range of PRS models across four ancestry groups (Africans, Europeans, East Asians, and South Asians) for 14 conditions of high-medical interest. Methods: To select the best-performing model per trait, we first compared PRS performances for publicly available scores, and constructed new models using different methods (LDpred2, PRS-CSx and SNPnet). We used 285K European individuals from the UK Biobank (UKBB) for training and 18K, including diverse ancestries, for testing. We then evaluated PRS portability for the best models in Europeans and compared their accuracies with respect to the best PRS per ancestry. Finally, we validated the selected PRS models using an independent set of 8,417 individuals from Biobank of the Americas-Genomelink (BBofA-GL); and performed a PRS-Phewas. Results: We confirmed a decay in PRS performances relative to Europeans when the evaluation was conducted using the best-PRS model for Europeans (51.3% for South Asians, 46.6% for East Asians and 39.4% for Africans). We observed an improvement in the PRS performances when specifically selecting ancestry specific PRS models (phenotype variance increase: 1.62 for Africans, 1.40 for South Asians and 0.96 for East Asians). Additionally, when we selected the optimal model conditional on ancestry for CAD, HDL-C and LDL-C, hypertension, hypothyroidism and T2D, PRS performance for studied populations was more comparable to what was observed in Europeans. Finally, we were able to independently validate tested models for Europeans, and conducted a PRS-Phewas, identifying cross-trait interplay between cardiometabolic conditions, and between immune-mediated components. Conclusion: Our work comprehensively evaluated PRS accuracy across a wide range of phenotypes, reducing the uncertainty with respect to which PRS model to choose and in which ancestry group. This evaluation has let us identify specific conditions where implementing risk-prioritization strategies could have practical utility across diverse ancestral groups, contributing to democratizing the implementation of PRS.

Micronesia began to be peopled earlier than other parts of Remote Oceania, but the origins of its inhabitants remain unclear. We generated genome-wide data from 164 ancient and 112 modern individuals. Analysis reveals five migratory streams into Micronesia. Three are East Asian related, one is Polynesian, and a fifth is a Papuan source related to mainland New Guineans that is different from the New Britain–related Papuan source for southwest Pacific populations but is similarly derived from male migrants ~2500 to 2000 years ago. People of the Mariana Archipelago may derive all of their precolonial ancestry from East Asian sources, making them the only Remote Oceanians without Papuan ancestry. Female-inherited mitochondrial DNA was highly differentiated across early Remote Oceanian communities but homogeneous within, implying matrilocal practices whereby women almost never raised their children in communities different from the ones in which they grew up.

Globally, autosomal recessive IFNAR1 deficiency is a rare inborn error of immunity underlying susceptibility to live attenuated vaccine and wild-type viruses. We report seven children from five unrelated kindreds of western Polynesian ancestry who suffered from severe viral diseases. All the patients are homozygous for the same nonsense IFNAR1 variant (p.Glu386*). This allele encodes a truncated protein that is absent from the cell surface and is loss-of-function. The fibroblasts of the patients do not respond to type I IFNs (IFN-α2, IFN-ω, or IFN-β). Remarkably, this IFNAR1 variant has a minor allele frequency >1% in Samoa and is also observed in the Cook, Society, Marquesas, and Austral islands, as well as Fiji, whereas it is extremely rare or absent in the other populations tested, including those of the Pacific region. Inherited IFNAR1 deficiency should be considered in individuals of Polynesian ancestry with severe viral illnesses.

The population of Mexico has a considerable genetic substructure due to both its pre-Columbian diversity and due to genetic admixture from post-Columbian trans-oceanic migrations. The latter primarily originated in Europe and Africa, but also, to a lesser extent, in Asia. We analyze previously understudied genetic connections between Asia and Mexico to infer the timing and source of this genetic ancestry in Mexico. We identify the predominant origin within Southeast Asia—specifically western Indonesian and non-Negrito Filipino sources—and we date its arrival in Mexico to approximately 13 generations ago (1620 CE). This points to a genetic legacy from the seventeenth century Manila galleon trade between the colonial Spanish Philippines and the Pacific port of Acapulco. Indeed, within Mexico we observe the highest level of this trans-Pacific ancestry in Acapulco, located in the state of Guerrero. This colonial Spanish trade route from East Asia to Europe was centred on Mexico and appears in historical records, but its legacy has been largely ignored. Identities and stories were suppressed due to slavery, assimilation of the immigrants as ‘Indios’ and incomplete historical records. Here we characterize this understudied Mexican ancestry. This article is part of the theme issue ‘Celebrating 50 years since Lewontin's apportionment of human diversity’.

Stochastic simulation is a key tool in population genetics, since the models involved are often analytically intractable and simulation is usually the only way of obtaining ground-truth data to evaluate inferences. Because of this, a large number of specialized simulation programs have been developed, each filling a particular niche, but with largely overlapping functionality and a substantial duplication of effort. Here, we introduce msprime version 1.0, which efficiently implements ancestry and mutation simulations based on the succinct tree sequence data structure and the tskit library. We summarize msprime’s many features, and show that its performance is excellent, often many times faster and more memory efficient than specialized alternatives. These high-performance features have been thoroughly tested and validated, and built using a collaborative, open source development model, which reduces duplication of effort and promotes software quality via community engagement.

Stochastic simulation is a key tool in population genetics, since the models involved are often analytically intractable and simulation is usually the only way of obtaining ground-truth data to evaluate inferences. Because of this necessity, a large number of specialised simulation programs have been developed, each filling a particular niche, but with largely overlapping functionality and a substantial duplication of effort. Here, we introduce msprime version 1.0, which efficiently implements ancestry and mutation simulations based on the succinct tree sequence data structure and tskit library. We summarise msprime ’s many features, and show that its performance is excellent, often many times faster and more memory efficient than specialised alternatives. These high-performance features have been thoroughly tested and validated, and built using a collaborative, open source development model, which reduces duplication of effort and promotes software quality via community engagement.

We used 15 short tandem repeat (STR) loci (D8S1179, D21S11, D7S820, CSF1PO, D3S1358, TH01, D13S317, D16S539, D2S1338, D19S433, VWA, TPOX, D18S51, D5S818, and FGA) to genetically characterize 361 individuals from 11 indigenous populations (Amuzgo, Chinanteco, Chontal, Huave, Mazateco, Mixe, Mixteco, Triqui, Zapoteco del Istmo, Zapoteco del Valle, and Zoque) from Oaxaca, Mexico. We also used previously published data from other Mexican peoples (Maya, Chol, Tepehua, Otomi, and Mestizos from northern and central Mexico) to delineate genetic relations, for a total of 541 individuals. Average heterozygosity (H) was lower in most populations from Oaxaca (range 0.687 in Zoque to 0.756 in Chontal) than values observed in Mestizo populations from Mexico (0.758 and 0.793 in central and northern Mestizo, respectively) but higher than values observed in other Amerindian populations from South America; the same relation was true for the number of alleles (n(a)). We tested (using the software Structure) whether major geographic or linguistic barriers to gene flow existed among the populations of Oaxaca and found that the populations appeared to constitute one or two genetic groups, suggesting that neither geographic location nor linguistics had an effect on the genetic structure of these culturally and linguistically highly diverse indigenous peoples. Moreover, we found a low but statistically significant between-population differentiation. In addition, the genetic structure of Oaxacan populations did not fit an isolation-by-distance model. Finally, using AMOVA and a Bayesian clustering approach, we did not detect significant geographic or linguistic barriers to gene flow within Oaxaca. These results suggest that the indigenous communities of Oaxaca, although culturally isolated, can be genetically defined as a large, nearly panmictic population in which migration could be a more important population mechanism than genetic drift. Finally, compared with outgroups in Mexico (both indigenous peoples and Mestizos), three groups were apparent. Among them, only the Otomi population from Hidalgo has a different culture and language.