Karla Sandoval

The genetics of Mexico recapitulates Native American substructure and affects biomedical traits

The population structure of Native Mexicans The genetics of indigenous Mexicans exhibit substantial geographical structure, some as divergent from each other as are existing populations of Europeans and Asians. By performing genome-wide analyses on Native Mexicans from differing populations, Moreno-Estrada et al. successfully recapitulated the pre-Columbian substructure of Mexico. This ancestral structure is evident among cosmopolitan Mexicans and is correlated with subcontinental origins and medically relevant aspects of lung function. These findings exemplify the importance of understanding the genetic contributions of admixed individuals. Science, this issue p. 1280 Indigenous and cosmopolitan Mexican populations are highly structured with genetic variation of medical relevance. 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 1000 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 subcontinental 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.

Pulling out the 1%: whole-Genome capture for the targeted enrichment of ancient dna sequencing libraries

Most ancient specimens contain very low levels of endogenous DNA, precluding the shotgun sequencing of many interesting samples because of cost. Ancient DNA (aDNA) libraries often contain <1% endogenous DNA, with the majority of sequencing capacity taken up by environmental DNA. Here we present a capture-based method for enriching the endogenous component of aDNA sequencing libraries. By using biotinylated RNA baits transcribed from genomic DNA libraries, we are able to capture DNA fragments from across the human genome. We demonstrate this method on libraries created from four Iron Age and Bronze Age human teeth from Bulgaria, as well as bone samples from seven Peruvian mummies and a Bronze Age hair sample from Denmark. Prior to capture, shotgun sequencing of these libraries yielded an average of 1.2% of reads mapping to the human genome (including duplicates). After capture, this fraction increased substantially, with up to 59% of reads mapped to human and enrichment ranging from 6- to 159-fold. Furthermore, we maintained coverage of the majority of regions sequenced in the precapture library. Intersection with the 1000 Genomes Project reference panel yielded an average of 50,723 SNPs (range 3,062-147,243) for the postcapture libraries sequenced with 1 million reads, compared with 13,280 SNPs (range 217-73,266) for the precapture libraries, increasing resolution in population genetic analyses. Our whole-genome capture approach makes it less costly to sequence aDNA from specimens containing very low levels of endogenous DNA, enabling the analysis of larger numbers of samples.

Reconstructing the population genetic history of the Caribbean.

The Caribbean basin is home to some of the most complex interactions in recent history among previously diverged human populations. Here, we investigate the population genetic history of this region by characterizing patterns of genome-wide variation among 330 individuals from three of the Greater Antilles (Cuba, Puerto Rico, Hispaniola), two mainland (Honduras, Colombia), and three Native South American (Yukpa, Bari, and Warao) populations. We combine these data with a unique database of genomic variation in over 3,000 individuals from diverse European, African, and Native American populations. We use local ancestry inference and tract length distributions to test different demographic scenarios for the pre- and post-colonial history of the region. We develop a novel ancestry-specific PCA (ASPCA) method to reconstruct the sub-continental origin of Native American, European, and African haplotypes from admixed genomes. We find that the most likely source of the indigenous ancestry in Caribbean islanders is a Native South American component shared among inland Amazonian tribes, Central America, and the Yucatan peninsula, suggesting extensive gene flow across the Caribbean in pre-Columbian times. We find evidence of two pulses of African migration. The first pulse—which today is reflected by shorter, older ancestry tracts—consists of a genetic component more similar to coastal West African regions involved in early stages of the trans-Atlantic slave trade. The second pulse—reflected by longer, younger tracts—is more similar to present-day West-Central African populations, supporting historical records of later transatlantic deportation. Surprisingly, we also identify a Latino-specific European component that has significantly diverged from its parental Iberian source populations, presumably as a result of small European founder population size. We demonstrate that the ancestral components in admixed genomes can be traced back to distinct sub-continental source populations with far greater resolution than previously thought, even when limited pre-Columbian Caribbean haplotypes have survived.

Reconstructing native American migrations from whole-genome and whole-exome data.

There is great scientific and popular interest in understanding the genetic history of populations in the Americas. We wish to understand when different regions of the continent were inhabited, where settlers came from, and how current inhabitants relate genetically to earlier populations. Recent studies unraveled parts of the genetic history of the continent using genotyping arrays and uniparental markers. The 1000 Genomes Project provides a unique opportunity for improving our understanding of population genetic history by providing over a hundred sequenced low coverage genomes and exomes from Colombian (CLM), Mexican-American (MXL), and Puerto Rican (PUR) populations. Here, we explore the genomic contributions of African, European, and especially Native American ancestry to these populations. Estimated Native American ancestry is in MXL, in CLM, and in PUR. Native American ancestry in PUR is most closely related to populations surrounding the Orinoco River basin, confirming the Southern America ancestry of the Taíno people of the Caribbean. We present new methods to estimate the allele frequencies in the Native American fraction of the populations, and model their distribution using a demographic model for three ancestral Native American populations. These ancestral populations likely split in close succession: the most likely scenario, based on a peopling of the Americas thousand years ago (kya), supports that the MXL Ancestors split kya, with a subsequent split of the ancestors to CLM and PUR kya. The model also features effective populations of in Mexico, in Colombia, and in Puerto Rico. Modeling Identity-by-descent (IBD) and ancestry tract length, we show that post-contact populations also differ markedly in their effective sizes and migration patterns, with Puerto Rico showing the smallest effective size and the earlier migration from Europe. Finally, we compare IBD and ancestry assignments to find evidence for relatedness among European founders to the three populations.

Genetic ancestry influences asthma susceptibility and lung function among Latinos.

BACKGROUND Childhood asthma prevalence and morbidity varies among Latinos in the United States, with Puerto Ricans having the highest and Mexicans the lowest. OBJECTIVE To determine whether genetic ancestry is associated with the odds of asthma among Latinos, and secondarily whether genetic ancestry is associated with lung function among Latino children. METHODS We analyzed 5493 Latinos with and without asthma from 3 independent studies. For each participant, we estimated the proportion of African, European, and Native American ancestry using genome-wide data. We tested whether genetic ancestry was associated with the presence of asthma and lung function among subjects with and without asthma. Odds ratios (OR) and effect sizes were assessed for every 20% increase in each ancestry. RESULTS Native American ancestry was associated with lower odds of asthma (OR = 0.72, 95% CI: 0.66-0.78, P = 8.0 × 10(-15)), while African ancestry was associated with higher odds of asthma (OR = 1.40, 95% CI: 1.14-1.72, P = .001). These associations were robust to adjustment for covariates related to early life exposures, air pollution, and socioeconomic status. Among children with asthma, African ancestry was associated with lower lung function, including both pre- and post-bronchodilator measures of FEV1 (-77 ± 19 mL; P = 5.8 × 10(-5) and -83 ± 19 mL; P = 1.1 x 10(-5), respectively) and forced vital capacity (-100 ± 21 mL; P = 2.7 × 10(-6) and -107 ± 22 mL; P = 1.0 x 10(-6), respectively). CONCLUSION Differences in the proportions of genetic ancestry can partially explain disparities in asthma susceptibility and lung function among Latinos.

Linguistic and maternal genetic diversity are not correlated in native mexicans

Mesoamerica, defined as the broad linguistic and cultural area from middle southern Mexico to Costa Rica, might have played a pivotal role during the colonization of the American continent. The Mesoamerican isthmus has constituted an important geographic barrier that has severely restricted gene flow between North and South America in pre-historical times. Although the Native American component has been already described in admixed Mexican populations, few studies have been carried out in native Mexican populations. In this study, we present mitochondrial DNA (mtDNA) sequence data for the first hypervariable region (HVR-I) in 477 unrelated individuals belonging to 11 different native populations from Mexico. Almost all of the Native Mexican mtDNAs could be classified into the four pan-Amerindian haplogroups (A2, B2, C1, and D1); only two of them could be allocated to the rare Native American lineage D4h3. Their haplogroup phylogenies are clearly star-like, as expected from relatively young populations that have experienced diverse episodes of genetic drift (e.g., extensive isolation, genetic drift, and founder effects) and posterior population expansions. In agreement with this observation, Native Mexican populations show a high degree of heterogeneity in their patterns of haplogroup frequencies. Haplogroup X2a was absent in our samples, supporting previous observations where this clade was only detected in the American northernmost areas. The search for identical sequences in the American continent shows that, although Native Mexican populations seem to show a closer relationship to North American populations, they cannot be related to a single geographical region within the continent. Finally, we did not find significant population structure in the maternal lineages when considering the four main and distinct linguistic groups represented in our Mexican samples (Oto-Manguean, Uto-Aztecan, Tarascan, and Mayan), suggesting that genetic divergence predates linguistic diversification in Mexico.

Y‐chromosome diversity in Native Mexicans reveals continental transition of genetic structure in the Americas

The genetic characterization of Native Mexicans is important to understand multiethnic based features influencing the medical genetics of present Mexican populations, as well as to the reconstruct the peopling of the Americas. We describe the Y-chromosome genetic diversity of 197 Native Mexicans from 11 populations and 1,044 individuals from 44 Native American populations after combining with publicly available data. We found extensive heterogeneity among Native Mexican populations and ample segregation of Q-M242* (46%) and Q-M3 (54%) haplogroups within Mexico. The northernmost sampled populations falling outside Mesoamerica (Pima and Tarahumara) showed a clear differentiation with respect to the other populations, which is in agreement with previous results from mtDNA lineages. However, our results point toward a complex genetic makeup of Native Mexicans whose maternal and paternal lineages reveal different narratives of their population history, with sex-biased continental contributions and different admixture proportions. At a continental scale, we found that Arctic populations and the northernmost groups from North America cluster together, but we did not find a clear differentiation within Mesoamerica and the rest of the continent, which coupled with the fact that the majority of individuals from Central and South American samples are restricted to the Q-M3 branch, supports the notion that most Native Americans from Mesoamerica southwards are descendants from a single wave of migration. This observation is compatible with the idea that present day Mexico might have constituted an area of transition in the diversification of paternal lineages during the colonization of the Americas.

Demographic History of Indigenous Populations in Mesoamerica Based on mtDNA Sequence Data

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. founder effects) or in recent periods (i.e. population declines). The Mesomerican regions showed differences in population fluctuation as indicated by the maximum Inter-Generational Rate (IGRmax): i) Center-South from the lithic period until the Preclassic; ii) West from the beginning of the Preclassic period until early Classic; iii) North characterized by a wide range of temporal variation from the Lithic to the Preclassic. Our findings are consistent with the genetic variations observed between central, South and Southeast Mesoamerica and the North-West region that are related to differences in genetic drift, structure, and temporal survival strategies (agriculture versus hunter-gathering, respectively). Interestingly, although the European contact had a major negative demographic impact, we detect a previous decline in Mesoamerica that had begun a few hundred years before.

Whole Genome Sequencing of Pharmacogenetic Drug Response in Racially Diverse Children with Asthma

Rationale: Albuterol, a bronchodilator medication, is the first‐line therapy for asthma worldwide. There are significant racial/ethnic differences in albuterol drug response. Objectives: To identify genetic variants important for bronchodilator drug response (BDR) in racially diverse children. Methods: We performed the first whole‐genome sequencing pharmacogenetics study from 1,441 children with asthma from the tails of the BDR distribution to identify genetic association with BDR. Measurements and Main Results: We identified population‐specific and shared genetic variants associated with BDR, including genome‐wide significant (P < 3.53 × 10−7) and suggestive (P < 7.06 × 10−6) loci near genes previously associated with lung capacity (DNAH5), immunity (NFKB1 and PLCB1), and &bgr;‐adrenergic signaling (ADAMTS3 and COX18). Functional analyses of the BDR‐associated SNP in NFKB1 revealed potential regulatory function in bronchial smooth muscle cells. The SNP is also an expression quantitative trait locus for a neighboring gene, SLC39A8. The lack of other asthma study populations with BDR and whole‐genome sequencing data on minority children makes it impossible to perform replication of our rare variant associations. Minority underrepresentation also poses significant challenges to identify age‐matched and population‐matched cohorts of sufficient sample size for replication of our common variant findings. Conclusions: The lack of minority data, despite a collaboration of eight universities and 13 individual laboratories, highlights the urgent need for a dedicated national effort to prioritize diversity in research. Our study expands the understanding of pharmacogenetic analyses in racially/ethnically diverse populations and advances the foundation for precision medicine in at‐risk and understudied minority populations.

Whole Genome Sequencing of Pharmacogenetic Drug Response in Racially and Ethnically Diverse Children with Asthma

Albuterol, a bronchodilator medication, is the first-line therapy for asthma treatment worldwide. However it has a wide variation of drug response among different racial/ethnic groups. We performed the largest pharmacogenetics study to date, using whole genome sequencing data from 1,441 minority children with asthma from the extremes of bronchodilator drug response (BDR) to albuterol. We identified population-specific and shared pharmacogenetic variants associated with BDR, including genome-wide significant and suggestive loci near genes previously associated with lung capacity (DNAH5), immunity (NFKB1 and PLCB1), and β-adrenergic signaling pathways (ADAMTS3 and COX18). Functional assays revealed that the BDR-associated SNP within NFKB1 is in linkage disequilibrium with SNPs in a functionally active enhancer and is also associated with the expression of a neighboring gene SLC39A8. Our study expands the understanding of pharmacogenetic analyses in racially and ethnically diverse populations and advances the foundation for precision medicine in at-risk and understudied minority populations.