Large scale mitochondrial sequencing in Mexican Americans suggests a reappraisal of Native American origins

Bioanthropological analysis of human remains from the archaic and classic period discovered in Puyil cave, Mexico

OBJECTIVES

Determine the geographic place of origin and maternal lineage of prehistoric human skeletal remains discovered in Puyil Cave, Tabasco State, Mexico, located in a region currently populated by Olmec, Zoque and Maya populations.

MATERIALS AND METHODS

All specimens were radiocarbon (14C) dated (beta analytic), had dental modifications classified, and had an analysis of 13 homologous reference points conducted to evaluate artificial cranial deformation (ACD). Following DNA purification, hypervariable region I (HVR-1) of the mitogenome was amplified and Sanger sequenced. Finally, Next Generation Sequencing (NGS) was performed for total DNA. Mitochondrial DNA (mtDNA) variants and haplogroups were determined using BioEdit 7.2 and IGV software and confirmed with MITOMASTER and WebHome softwares.

RESULTS

Radiocarbon dating (14C) demonstrated that the inhabitants of Puyil Cave lived during the Archaic and Classic Periods and displayed tabular oblique and tabular mimetic ACD. These pre-Hispanic remains exhibited five mtDNA lineages: A, A2, C1, C1c and D4. Network analysis revealed a close genetic affinity between pre-Hispanic Puyil Cave inhabitants and contemporary Maya subpopulations from Mexico and Guatemala, as well as individuals from Bolivia, Brazil, the Dominican Republic, and China.

CONCLUSIONS

Our results elucidate the dispersal of pre-Hispanic Olmec and Maya ancestors and suggest that ACD practices are closely related to Olmec and Maya practices. Additionally, we conclude that ACD has likely been practiced in the region since the Middle-Archaic Period.

High Mitochondrial Haplotype Diversity Found in Three Pre-Hispanic Groups from Colombia

The analysis of mitochondrial DNA (mtDNA) hypervariable region (HVR) sequence data from ancient human remains provides valuable insights into the genetic structure and population dynamics of ancient populations. mtDNA is particularly useful in studying ancient populations, because it is maternally inherited and has a higher mutation rate compared to nuclear DNA. To determine the genetic structure of three Colombian pre-Hispanic populations and compare them with current populations, we determined the haplotypes from human bone remains by sequencing several mitochondrial DNA segments. A wide variety of mitochondrial polymorphisms were obtained from 33 samples. Our results support a high population heterogeneity among pre-Hispanic populations in Colombia.

Demographic history and genetic structure in pre-Hispanic Central Mexico

Aridoamerica and Mesoamerica are two distinct cultural areas in northern and central Mexico, respectively, that hosted numerous pre-Hispanic civilizations between 2500 BCE and 1521 CE. The division between these regions shifted southward because of severe droughts ~1100 years ago, which allegedly drove a population replacement in central Mexico by Aridoamerican peoples. In this study, we present shotgun genome-wide data from 12 individuals and 27 mitochondrial genomes from eight pre-Hispanic archaeological sites across Mexico, including two at the shifting border of Aridoamerica and Mesoamerica. We find population continuity that spans the climate change episode and a broad preservation of the genetic structure across present-day Mexico for the past 2300 years. Lastly, we identify a contribution to pre-Hispanic populations of northern and central Mexico from two ancient unsampled "ghost" populations.

Bioarchaeological study of ancient Teotihuacans based on complete mitochondrial genome sequences and diet isotopes

BACKGROUND

The Teotihuacan civilisation was the largest one in ancient Mesoamerica. The Teotihuacan city was born in the north-eastern Basin of Mexico around the second century BC, reached its peak in the fourth century AD, and had cultural influence throughout Mesoamerica. At its peak, the size of the city reached more than 20 km2, and the total population is estimated to have increased from 100,000 to 200,000. However, knowledge of the genetic background of the Teotihuacan people is still limited.

AIM

We aimed to determine the mitogenome sequences of the Teotihuacan human remains and compare the ancient and present Mesoamericans. In addition, we aimed to identify the food habits of ancient Teotihuacans.

SUBJECTS AND METHODS

We determined the mitogenome sequences of human remains dated to 250-636 cal AD using target enrichment-coupled next generation sequencing. We also performed stable isotope analysis.

RESULTS

We successfully obtained nearly full-length sequences newly unearthed from a civilian dwelling in the Teotihuacan site. Teotihuacan mitochondrial DNA was classified into the haplogroups in present and ancient Mesoamericans. In addition, Teotihuacan individuals had a diet dependent on C4 plants such as maize.

CONCLUSION

Genetic diversity varied among the Teotihuacans.

Fitness and Power: The Contribution of Genetics to the History of Differential Reproduction

Textual evidence from pre-modern societies supports the prediction that status differences among men translate to variance in reproductive success. In recent years, analysis of genetic data has opened up new ways of studying this relationship. By investigating cases that range over several millennia, these analyses repeatedly document the replacement of local men by newcomers and reveal instances of exceptional reproductive success of specific male lineages. These findings suggest that violent population transfers and conquests could generate considerable reproductive advantages for male dominants. At the same time, this does not always seem to have been the case. Moreover, it is difficult to link such outcomes to particular historical characters or events, or to identify status-biased reproductive inequalities within dominant groups. The proximate factors that mediated implied imbalances in reproductive success often remain unclear. A better understanding of the complex interplay between social power and genetic fitness will only arise from sustained transdisciplinary engagement.

The Mitochondrial DNA Landscape of Modern Mexico

Mexico is a rich source for anthropological and population genetic studies with high diversity in ethnic and linguistic groups. The country witnessed the rise and fall of major civilizations, including the Maya and Aztec, but resulting from European colonization, the population landscape has dramatically changed. Today, the majority of Mexicans do not identify themselves as Indigenous but as admixed, and appear to have very little in common with their pre-Columbian predecessors. However, when the maternally inherited mitochondrial (mt)DNA is investigated in the modern Mexican population, this is not the case. Control region sequences of 2021 samples deriving from all over the country revealed an overwhelming Indigenous American legacy, with almost 90% of mtDNAs belonging to the four major pan-American haplogroups A2, B2, C1, and D1. This finding supports a very low European contribution to the Mexican gene pool by female colonizers and confirms the effectiveness of employing uniparental markers as a tool to reconstruct a country’s history. In addition, the distinct frequency and dispersal patterns of Indigenous American and West Eurasian clades highlight the benefit such large and country-wide databases provide for studying the impact of colonialism from a female perspective and population stratification. The importance of geographical database subsets not only for forensic application is clearly demonstrated.

Cultivating Positive Health, Learning, and Community: The Return of Mesoamerica’s Quetzalcoatl and the Venus Star

For more than 3500 years, since Olmec times (1500–400 BC), the peoples of Mesoamerica have shared with one another a profound way of living involving a deep understanding of the human body and of land and cosmology. As it stands, healing ways of knowing that depend on medicinal plants, the Earth’s elements, and knowledge of the stars are still intact. The Indigenous Xicana/o/xs who belong to many of the mobile tribes of Mesoamerica share a long genealogical history of cultivating and sustaining their Native American rituals, which was weakened in Mexico and the United States during various periods of colonization. This special edition essay sheds light on the story of Quetzalcoatl and the Venus Star as a familial place of Xicana/o/x belonging and practice. To do so, we rely on the archaeological interpretation of these two entities as one may get to know them through artifacts, monuments, and ethnographic accounts, of which some date to Mesoamerica’s Formative period (1500–400 BC). Throughout this paper, ancestral medicine ways are shown to help cultivate positive health, learning, and community. Such cosmic knowledge is poorly understood, yet it may further culturally relevant education and the treatment of the rampant health disparities in communities of Mesoamerican ancestry living in the United States. The values of and insights into Indigenous Xicana/o/x knowledge and identity conclude this essay.

Y chromosome diversity in Aztlan descendants and its implications for the history of Central Mexico

Native Mexican populations are crucial for understanding the genetic ancestry of Aztec descendants and coexisting ethnolinguistic groups in the Valley of Mexico and elucidating the population dynamics of the prehistoric colonization of the Americas. Mesoamerican societies were multicultural in nature and also experienced significant admixture during Spanish colonization of the region. Despite these facts, Native Mexican Y chromosome diversity has been greatly understudied. To further elucidate their genetic history, we conducted a high-resolution Y chromosome analysis with Chichimecas, Nahuas, Otomies, Popolocas, Tepehuas, and Totonacas using 19 Y-short tandem repeat and 21 single nucleotide polymorphism loci. We detected enormous paternal genetic diversity in these groups, with haplogroups Q-MEH2, Q-M3, Q-Z768, Q-L663, Q-Z780, and Q-PV3 being identified. These data affirmed the southward colonization of the Americas via Beringia and connected Native Mexicans with indigenous populations from South-Central Siberia and Canada. They also suggested that multiple population dispersals gave rise to Y chromosome diversity in these populations.

•

Enormous Y chromosome diversity observed in Native Mexican populations.

•

Haplogroups Q-MEH2, Q-M3, Q-Z768, Q-L663, Q-Z780, and Q-PV3 were identified.

•

Patterns of Y chromosome diversity not shaped by ethnicity, geography, or language.

•

Multiple population dispersals contributed to Y chromosome diversity in Mexico.

A New Approach to Identify the Methylation Sites in the Control Region of Mitochondrial DNA

Mitochondrial DNA (mtDNA) methylation has the potential to be used as a biomarker of human development or disease. However, mtDNA methylation procedures are costly and time-consuming. Therefore, we developed a new approach based on an RT-PCR assay for the base site identification of methylated cytosine in the control region of mtDNA through a simple, fast, specific, and low-cost strategy. Total DNA was purified, and methylation was determined by RT-PCR bisulfite sequencing. This procedure included the DNA purification, bisulfite treatment and RT-PCR amplification of the control region divided into three subregions with specific primers. Sequences obtained with and without the bisulfite treatment were compared to identify the methylated cytosine dinucleotides. Furthermore, the efficiency of C to U conversion of cytosines was assessed by including a negative control. Interestingly, mtDNA methylation was observed mainly within non-Cphosphate- G (non-CpG) dinucleotides and mostly in the regions containing regulatory elements, such as OH or CSBI, CSBII, and CSBIII. This new approach will promote the generation of new information regarding mtDNA methylation patterns in samples from patients with different pathologies or that are exposed to a toxic environment in diverse human populations.

A Mitochondrial DNA Variant Elevates the Risk of Gallstone Disease by Altering Mitochondrial Function

BACKGROUND AND AIMS

Gallstone disease (cholelithiasis) is a cholesterol-related metabolic disorders with strong familial predisposition. Mitochondrial DNA (mtDNA) variants accumulated during human evolution are associated with some metabolic disorders related to modified mitochondrial function. The mechanistic links between mtDNA variants and gallstone formation need further exploration.

METHODS

In this study, we explored the possible associations of mtDNA variants with gallstone disease by comparing 104 probands and 300 controls in a Chinese population. We constructed corresponding cybrids using trans-mitochondrial technology to investigate the underlying mechanisms of these associations. Mitochondrial respiratory chain complex activity and function and cholesterol metabolism were assessed in the trans-mitochondrial cell models.

RESULTS

Here, we found a significant association of mtDNA 827A>G with an increased risk of familial gallstone disease in a Chinese population (odds ratio [OR]: 4.5, 95% confidence interval [CI]: 2.1-9.4, P=1.2×10-4). Compared with 827A cybrids (haplogroups B4a and B4c), 827G cybrids (haplogroups B4b and B4d) had impaired mitochondrial respiratory chain complex activity and function and activated JNK and AMPK signaling pathways. Additionally, the 827G cybrids showed disturbances in cholesterol transport and accelerated development of gallstones. Specifically, cholesterol transport through the transporter ABCG5/8 was increased via activation of the AMPK signaling pathway in 827G cybrids.

CONCLUSIONS

Our findings reveal that mtDNA 827A>G induces aberrant mitochondrial function and abnormal cholesterol transport, resulting in increased occurrence of gallstones. The results provide an important biological basis for the clinical diagnosis and prevention of gallstone disease in the future.

Ancient DNA Reconstructs the Genetic Legacies of Precontact Puerto Rico Communities

Indigenous peoples have occupied the island of Puerto Rico since at least 3000 BC. Due to the demographic shifts that occurred after European contact, the origin(s) of these ancient populations, and their genetic relationship to present-day islanders, are unclear. We use ancient DNA to characterize the population history and genetic legacies of precontact Indigenous communities from Puerto Rico. Bone, tooth, and dental calculus samples were collected from 124 individuals from three precontact archaeological sites: Tibes, Punta Candelero, and Paso del Indio. Despite poor DNA preservation, we used target enrichment and high-throughput sequencing to obtain complete mitochondrial genomes (mtDNA) from 45 individuals and autosomal genotypes from two individuals. We found a high proportion of Native American mtDNA haplogroups A2 and C1 in the precontact Puerto Rico sample (40% and 44%, respectively). This distribution, as well as the haplotypes represented, supports a primarily Amazonian South American origin for these populations and mirrors the Native American mtDNA diversity patterns found in present-day islanders. Three mtDNA haplotypes from precontact Puerto Rico persist among Puerto Ricans and other Caribbean islanders, indicating that present-day populations are reservoirs of precontact mtDNA diversity. Lastly, we find similarity in autosomal ancestry patterns between precontact individuals from Puerto Rico and the Bahamas, suggesting a shared component of Indigenous Caribbean ancestry with close affinity to South American populations. Our findings contribute to a more complete reconstruction of precontact Caribbean population history and explore the role of Indigenous peoples in shaping the biocultural diversity of present-day Puerto Ricans and other Caribbean islanders.

Reduced mitochondrial DNA copy number is associated with the haplogroup, and some clinical features of breast cancer in Mexican patients

Mitochondrial DNA (mtDNA) copy number and mitochondrial DNA haplogroups have been associated with different types of cancer, including breast cancer, because they alter cellular energy metabolism. However, whether mtDNA copy number or haplogroups are predictors of oxidative stress-related risks in human breast cancer tissue in Mexican patients remains to be determined. Using quantitative real-time PCR assays and sequencing of the mtDNA hypervariable region, analysis of mtDNA copy numbers in 82 breast cancer tissues (BCT) and matched normal adjacent tissues (NAT) was performed to determine if copy number correlated with clinical features and Amerindian haplogroups (A2, B2, B4, C1 and D1) . The results showed that the mtDNA copy number was significantly decreased in BCT compared with NAT (p = 0.010); it was significantly decreased in BCT and NAT in women > 50 years of age, compared with NAT in women < 50 years of age (p = 0.032 and p = 0.037, respectively); it was significantly decreased in NAT and BCT in the postmenopausal group and in BCT in the premenopausal group compared with NAT in the premenopausal group (p = 0.011, p = 0.010 and, p = 0.018; respectively); and it was also significantly decrease in members of the BCT group classified as having invasive ductal carcinoma I-III (IDC-I, IDC-II and IDC-III) and IDC-II for NAT compared to IDC-I of NAT (p = 0.025, p = 0.022 and p = 0.031 and p = 0.020; respectively). The mtDNA copy number for BCT from patients with haplogroup B2 was decreased compared to patients with haplogroup D1 (p = 0.01); for BCT from patients with haplogroup C1 was also decreased compare with their NAT counterpart (p = 0.006) and with BCT patients belonging to haplogroups A2 and D1 (p = 0.01 and p = 0.03; respectively). In addition, the mtDNA copy number was decrease in the sequences with three deletions relative to the rCRS at nucleotide positions A249del, A290del and A291del, or C16327T polymorphism with the same p = 0.019 for all four variants. Contrary, the copy number increased in sequences containing C16111T, G16319A or T16362C polymorphisms (p = 0.021, =0.048, and = 0.001; respectively). In conclusion, a decrease in the copy number of mtDNA in BCT compared with NAT was shown by the results, which suggests an imbalance in oxidative phosphorylation (OXPHOS) that can affect the apoptosis pathway and cancer progression. It was also observed an increase of the copy number in samples with specific polymorphisms, which may be a good sign of favourable prognosis.

Adiposity and Leukocyte Telomere Length in US Adults by Sex-Specific Race/Ethnicity: National Health and Nutrition Examination Survey

Objective

Little is known about the relationship between adiposity and telomere length in the United States population. The objective of our research was to examine this relationship in a representative, socioeconomically and sex-specific, diverse racial/ethnic population in the United States.

Methods

Body mass index (BMI), % total body fat (TBF) and waist circumference (WC) with leukocyte telomere length (LTL) were examined according to sex-specific race/ethnicity using separate adjusted multivariate linear regressions on a sample of 4,919 respondents aged 20-84 years from the National Health and Nutrition Examination Survey's 1999-2002 data.

Results

LTL was shortened .41%, .44%, and .16% in African American (AA) women and was associated with increasing BMI, %TBF, and WC, (β:-.0041, 95%CI: -.0070, -.0012; P=.007; β:-.0044, 95% CI: -.0081, -.0007; P=.02; β:-.0016, 95%CI: -.0031, -.0001; P=.04, respectively). LTL was shortened .29% in White women and was associated with increasing %TBF (β:-.0029, 95%CI: -.0048, -.0009; P=.006). There were no associations among AA men, White men or Mexican American men and women.

Conclusions

LTL is associated with an obesity phenotype in AA women. Tailored intervention is needed to ameliorate the burden of excess adiposity and subsequent cellular aging.

Mitochondrial DNA, a Powerful Tool to Decipher Ancient Human Civilization from Domestication to Music, and to Uncover Historical Murder Cases

Mitochondria are unique organelles carrying their own genetic material, independent from that in the nucleus. This review will discuss the nature of mitochondrial DNA (mtDNA) and its levels in the cell, which are the key elements to consider when trying to achieve molecular identification in ancient and degraded samples. mtDNA sequence analysis has been appropriately validated and is a consistent molecular target for the examination of biological evidence encountered in forensic cases—and profiling, in certain conditions—especially for burnt bodies and degraded samples of all types. Exceptional cases and samples will be discussed in this review, such as mtDNA from leather in Beethoven’s grand piano, mtDNA in mummies, and solving famous historical criminal cases. In addition, this review will be discussing the use of ancient mtDNA to understand past human diet, to trace historical civilizations and ancient trade routes, and to uncover geographical domestication origins and lineage relationships. In each topic, we will present the power of mtDNA and how, in many cases, no nuclear DNA was left, leaving mitochondrial DNA analysis as a powerful alternative. Exploring this powerful tool further will be extremely useful to modern science and researchers, due to its capabilities in providing us with previously unattainable knowledge.

The utility of genomic public databases to mitochondrial haplotyping in contemporary Mestizo population of Mexican origin

There are different public databases and open access information that can be exploited to be reused in different research projects. With this concept in mind, we carried out a study to answer the question about the prevalence of haplogroups in human populations of modern Mexico. Since the publication of genomic and mitochondrial data in Latin American populations are very scarce and with very small samples, our work proposes to consider the availability of genomic and genetic data collections that can be reused for other purposes, different from those initially proposed in the investigations where the sequences were obtained. The objective of the present study was to explore the population structure of Mexico using available information in the public database. Through the search of information in the nucleotide database of National Center of Biotechnology Information (NCBI) of complete sequences of mitochondrial genome (16 Kb) of indigenous people, Mexican Mestizo population and Mexican-Americans living in the United States, they were classified according to the polymorphisms associated with haplogroups A, B, C and D reported in the literature as the most frequent. We obtained 283 sequences, of which 255 were selected with the criteria mentioned. The haplotyping results showed 113 different clades and subclades distributed in a general way in eight haplogroups. The most frequent groups that dominate the population were the haplogroup A with 90 individuals representing 36%, followed by haplogroup B in 65 individuals representing 26% of the sample.

Whole sequence of the mitochondrial DNA genome of Kearns Sayre Syndrome patients: Identification of deletions and variants

Mitochondria both produce the energy of the cell as ATP via respiration and regulate cellular metabolism. Accordingly, any deletion or mutation in the mitochondrial DNA (mtDNA) may result in a disease. One of these diseases is Kearns Sayre syndrome (KSS), described for the first time in 1958, where different large-scale deletions of different sizes and at different positions have been reported in the mitochondrial genome of patients with similar clinical symptoms. In this study, sequences of the mitochondrial genome of three patients with clinic features of KSS were analyzed. Our results revealed the position, heteroplasmy percentage, size of deletions, and their haplogroups. Two patients contained deletions reported previously and one patient showed a new deletion not reported previously. These results display for the first time a systematic analysis of mtDNA variants in the whole mtDNA genome of patients with KSS to help to understand their association with the disease.

Early human dispersals within the Americas

Studies of the peopling of the Americas have focused on the timing and number of initial migrations. Less attention has been paid to the subsequent spread of people within the Americas. We sequenced 15 ancient human genomes spanning from Alaska to Patagonia; six are ≥10,000 years old (up to ~18× coverage). All are most closely related to Native Americans, including those from an Ancient Beringian individual and two morphologically distinct “Paleoamericans.” We found evidence of rapid dispersal and early diversification that included previously unknown groups as people moved south. This resulted in multiple independent, geographically uneven migrations, including one that provides clues of a Late Pleistocene Australasian genetic signal, as well as a later Mesoamerican-related expansion. These led to complex and dynamic population histories from North to South America.

Cultural Innovations Influence Patterns of Genetic Diversity in Northwestern Amazonia

Human populations often exhibit contrasting patterns of genetic diversity in the mtDNA and the nonrecombining portion of the Y-chromosome (NRY), which reflect sex-specific cultural behaviors and population histories. Here, we sequenced 2.3 Mb of the NRY from 284 individuals representing more than 30 Native American groups from Northwestern Amazonia (NWA) and compared these data to previously generated mtDNA genomes from the same groups, to investigate the impact of cultural practices on genetic diversity and gain new insights about NWA population history. Relevant cultural practices in NWA include postmarital residential rules and linguistic exogamy, a marital practice in which men are required to marry women speaking a different language. We identified 2,969 SNPs in the NRY sequences, only 925 of which were previously described. The NRY and mtDNA data showed different sex-specific demographic histories: female effective population size has been larger than that of males through time, which might reflect larger variance in male reproductive success. Both markers show an increase in lineage diversification beginning ∼5,000 years ago, which may reflect the intensification of agriculture, technological innovations, and the expansion of regional trade networks documented in the archaeological evidence. Furthermore, we find similar excesses of NRY versus mtDNA between-population divergence at both the local and continental scale, suggesting long-term stability of female versus male migration. We also find evidence of the impact of sociocultural practices on diversity patterns. Finally, our study highlights the importance of analyzing high-resolution mtDNA and NRY sequences to reconstruct demographic history, since this can differ considerably between sexes.

Historical records under the genetic evidence: "Chiriguano" tribe genesis as a test case

Historical records suggest that Chiriguano tribe is the result of a genetic admixture event. The process involved the arrival of Guaraní tribesmen descending from Amazonian region of Brazil along with groups of Arawak origin that inhabited the foothill plains of Bolivia. Later they arrived in Argentina at the beginning of the twentieth century. Aiming to test the historical records, we analysed a set of 46 samples collected at San Ramon de la Nueva Orán, Province of Salta, Argentina. A wide set of uni- and biparentally transmitted genetic markers were analysed, including 23 autosomal STRs; 46 AIM-DIPs and 24 AIM-SNPs all located at diverse autosomal chromosome locations; 23 Y-STRs and the entire mtDNA D-Loop sequence. Ancestry informative markers allowed for the detection of a strong Native American component in the genomes (> 94%), while all mtDNA haplotypes showed Native American characteristic motives, and 93% of Y-haplotypes belonged to the Q1a3a Y-haplogroup. The analysis of mitochondrial haplotypes and Y chromosome, although they did not match other populations, revealed a relationship between the Chiriguano and other groups of Guaraní and Arawak origin inhabiting Brazil and Bolivia, confirming, at least in part, the historical records describing the origins of Chiriguano tribal settlements in northwestern Argentina.

Ancient human mitochondrial DNA and radiocarbon analysis of archived quids from the Mule Spring Rockshelter, Nevada, USA

Chewed and expectorated quids, indigestible stringy fibers from the roasted inner pulp of agave or yucca root, have proven resilient over long periods of time in dry cave environments and correspondingly, although little studied, are common in archaeological archives. In the late 1960s, thousands of quids were recovered from Mule Spring Rockshelter (Nevada, USA) deposits and stored without consideration to DNA preservation in a museum collection, remaining unstudied for over fifty years. To assess the utility of these materials as repositories for genetic information about past inhabitants of the region and their movements, twenty-one quids were selected from arbitrary excavation depths for detailed analysis. Human mitochondrial DNA sequences from the quids were amplified by PCR and screened for diagnostic single nucleotide polymorphisms. Most detected single nucleotide polymorphisms were consistent with recognized Native American haplogroup subclades B2a5, B2i1, C1, C1c, C1c2, and D1; with the majority of the sample set consistent with subclades C1, C1c, and C1c2. In parallel with the DNA analysis, each quid was radiocarbon dated, revealing a time-resolved pattern of occupancy from 347 to 977 calibrated years before present. In particular, this dataset reveals strong evidence for the presence of haplogroup C1/C1c at the Southwestern edge of the US Great Basin from ~670 to 980 cal YBP, which may temporally correspond with the beginnings of the so-called Numic Spread into the region. The research described here demonstrates an approach which combines targeted DNA analysis with radiocarbon age dating; thus enabling the genetic analysis of archaeological materials of uncertain stratigraphic context. Here we present a survey of the maternal genetic profiles from people who used the Mule Spring Rockshelter and the historic timing of their utilization of a key natural resource.

The Paleo-Indian Entry into South America According to Mitogenomes

Recent and compelling archaeological evidence attests to human presence ∼14.5 ka at multiple sites in South America and a very early exploitation of extreme high-altitude Andean environments. Considering that, according to genetic evidence, human entry into North America from Beringia most likely occurred ∼16 ka, these archeological findings would imply an extremely rapid spread along the double continent. To shed light on this issue from a genetic perspective, we first completely sequenced 217 novel modern mitogenomes of Native American ancestry from the northwestern area of South America (Ecuador and Peru); we then evaluated them phylogenetically together with other available mitogenomes (430 samples, both modern and ancient) from the same geographic area and, finally, with all closely related mitogenomes from the entire double continent. We detected a large number (N = 48) of novel subhaplogroups, often branching into further subclades, belonging to two classes: those that arose in South America early after its peopling and those that instead originated in North or Central America and reached South America with the first settlers. Coalescence age estimates for these subhaplogroups provide time boundaries indicating that early Paleo-Indians probably moved from North America to the area corresponding to modern Ecuador and Peru over the short time frame of ∼1.5 ka comprised between 16.0 and 14.6 ka.

Successful reconstruction of whole mitochondrial genomes from ancient Central America and Mexico

The northern and southern peripheries of ancient Mesoamerica are poorly understood. There has been speculation over whether borderland cultures such as Greater Nicoya and Casas Grandes represent Mesoamerican outposts in the Isthmo-Colombian area and the Greater Southwest, respectively. Poor ancient DNA preservation in these regions challenged previous attempts to resolve these questions using conventional genetic techniques. We apply advanced in-solution mitogenome capture and high-throughput sequencing to fourteen dental samples obtained from the Greater Nicoya sites of Jícaro and La Cascabel in northwest Costa Rica (n = 9; A.D. 800–1250) and the Casas Grandes sites of Paquimé and Convento in northwest Mexico (n = 5; A.D. 1200–1450). Full mitogenome reconstruction was successful for three individuals from Jícaro and five individuals from Paquimé and Convento. The three Jícaro individuals belong to haplogroup B2d, a haplogroup found today only among Central American Chibchan-speakers. The five Paquimé and Convento individuals belong to haplogroups C1c1a, C1c5, B2f and B2a which, are found in contemporary populations in North America and Mesoamerica. We report the first successfully reconstructed ancient mitogenomes from Central America, and the first genetic evidence of ancestry affinity of the ancient inhabitants of Greater Nicoya and Casas Grandes with contemporary Isthmo-Columbian and Greater Southwest populations, respectively.

Enclaves of genetic diversity resisted Inca impacts on population history

The Inca Empire is claimed to have driven massive population movements in western South America, and to have spread Quechua, the most widely-spoken language family of the indigenous Americas. A test-case is the Chachapoyas region of northern Peru, reported as a focal point of Inca population displacements. Chachapoyas also spans the environmental, cultural and demographic divides between Amazonia and the Andes, and stands along the lowest-altitude corridor from the rainforest to the Pacific coast. Following a sampling strategy informed by linguistic data, we collected 119 samples, analysed for full mtDNA genomes and Y-chromosome STRs. We report a high indigenous component, which stands apart from the network of intense genetic exchange in the core central zone of Andean civilization, and is also distinct from neighbouring populations. This unique genetic profile challenges the routine assumption of large-scale population relocations by the Incas. Furthermore, speakers of Chachapoyas Quechua are found to share no particular genetic similarity or gene-flow with Quechua speakers elsewhere, suggesting that here the language spread primarily by cultural diffusion, not migration. Our results demonstrate how population genetics, when fully guided by the archaeological, historical and linguistic records, can inform multiple disciplines within anthropology.

Tri-allelic heteroplasmies, DNA-RNA differences and their polynucleotide tract associations in the mitochondrial genome

The human mitochondrial genome has been extensively studied for its function and disease associations. Utilizing five types of high-throughput sequencing data on ten breast cancer patients (total N=50), we examined several aspects of the mitochondrial genome that have not been thoroughly studied, including the occurrence of tri-allelic heteroplasmy, the difference between DNA and RNA, and the variants association with polynucleotide tracts. We validated four previously reported and identified 23 additional tri-allelic positions. Furthermore, we detected 18 single nucleotide and seven InDel differences between DNA and RNA. Previous studies have suggested that some of these differences are caused by post transcriptional methylation. The rest can be accredited to RNA editing, polyadenylation or sequencing errors. Most importantly, we found that the tri-allelic positions, and differences between DNA and RNA, are strongly associated with polynucleotide tracts in the mitochondrial genome, suggesting DNA instability or difficulty sequencing around the polynucleotide tract regions.

Characterization of complete mitochondrial genomes of indigenous Mayans in Mexico

BACKGROUND

The authors have previously published the complete mitochondrial genome (mitogenome) sequences of two indigenous Mesoamerican populations, Mazahua (n = 25) and Zapotec (n = 88).

METHODS

This study determined the complete mitogenome sequences of nine unrelated individuals from the indigenous Maya population living in Mexico.

RESULTS

Their mitogenome sequences could be classified into either of the haplogroups A2 and C1. Surprisingly, there were no mitogenome sequences (haplotypes) that the Maya, Mazahua, and Zapotec people share in common.

CONCLUSIONS

This indicates that no genetic exchange, at least matrilineally, has occurred among them.

Archaeogenomic evidence reveals prehistoric matrilineal dynasty

For societies with writing systems, hereditary leadership is documented as one of the hallmarks of early political complexity and governance. In contrast, it is unknown whether hereditary succession played a role in the early formation of prehistoric complex societies that lacked writing. Here we use an archaeogenomic approach to identify an elite matriline that persisted between 800 and 1130 CE in Chaco Canyon, the centre of an expansive prehistoric complex society in the Southwestern United States. We show that nine individuals buried in an elite crypt at Pueblo Bonito, the largest structure in the canyon, have identical mitochondrial genomes. Analyses of nuclear genome data from six samples with the highest DNA preservation demonstrate mother–daughter and grandmother–grandson relationships, evidence for a multigenerational matrilineal descent group. Together, these results demonstrate the persistence of an elite matriline in Chaco for ∼330 years.

In ancient cultures without a writing system, it is difficult to infer the basis of status and rank. Here the authors analyse ancient DNA from nine presumed elite individuals buried successively over a 300-year period at Chaco Canyon, and show evidence of matrilineal relationships.

Imputation approach for deducing a complete mitogenome sequence from low-depth-coverage next-generation sequencing data: application to ancient remains from the Moon Pyramid, Mexico

It is considered that more than 15 depths of coverage are necessary for next-generation sequencing (NGS) data to obtain reliable complete nucleotide sequences of the mitogenome. However, it is difficult to satisfy this requirement for all nucleotide positions because of problems obtaining a uniform depth of coverage for poorly preserved materials. Thus, we propose an imputation approach that allows a complete mitogenome sequence to be deduced from low-depth-coverage NGS data. We used different types of mitogenome data files as panels for imputation: a worldwide panel comprising all the major haplogroups, a worldwide panel comprising sequences belonging to the estimated haplogroup alone, a panel comprising sequences from the population most closely related to an individual under investigation, and a panel comprising sequences belonging to the estimated haplogroup from the population most closely related to an individual under investigation. The number of missing nucleotides was drastically reduced in all the panels, but the contents obtained by imputation were quite different among the panels. The efficiency of the imputation method differed according to the panels used. The missing nucleotides were most credibly imputed using sequences of the estimated haplogroup from the population most closely related to the individual under investigation as a panel.

Ancient mitochondrial DNA provides high-resolution time scale of the peopling of the Americas

The exact timing, route, and process of the initial peopling of the Americas remains uncertain despite much research. Archaeological evidence indicates the presence of humans as far as southern Chile by 14.6 thousand years ago (ka), shortly after the Pleistocene ice sheets blocking access from eastern Beringia began to retreat. Genetic estimates of the timing and route of entry have been constrained by the lack of suitable calibration points and low genetic diversity of Native Americans. We sequenced 92 whole mitochondrial genomes from pre-Columbian South American skeletons dating from 8.6 to 0.5 ka, allowing a detailed, temporally calibrated reconstruction of the peopling of the Americas in a Bayesian coalescent analysis. The data suggest that a small population entered the Americas via a coastal route around 16.0 ka, following previous isolation in eastern Beringia for ~2.4 to 9 thousand years after separation from eastern Siberian populations. Following a rapid movement throughout the Americas, limited gene flow in South America resulted in a marked phylogeographic structure of populations, which persisted through time. All of the ancient mitochondrial lineages detected in this study were absent from modern data sets, suggesting a high extinction rate. To investigate this further, we applied a novel principal components multiple logistic regression test to Bayesian serial coalescent simulations. The analysis supported a scenario in which European colonization caused a substantial loss of pre-Columbian lineages.

Why are there race/ethnic differences in adult body mass index-adiposity relationships? A quantitative critical review

Body mass index (BMI) is now the most widely used measure of adiposity on a global scale. Nevertheless, intense discussion centers on the appropriateness of BMI as a phenotypic marker of adiposity across populations differing in race and ethnicity. BMI-adiposity relations appear to vary significantly across race/ethnic groups, but a collective critical analysis of these effects establishing their magnitude and underlying body shape/composition basis is lacking. Accordingly, we systematically review the magnitude of these race-ethnic differences across non-Hispanic (NH) white, NH black and Mexican American adults, their anatomic body composition basis and potential biologically linked mechanisms, using both earlier publications and new analyses from the US National Health and Nutrition Examination Survey. Our collective observations provide a new framework for critically evaluating the quantitative relations between BMI and adiposity across groups differing in race and ethnicity; reveal new insights into BMI as a measure of adiposity across the adult age-span; identify knowledge gaps that can form the basis of future research and create a quantitative foundation for developing BMI-related public health recommendations.

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.

Mitochondrial genome diversity at the Bering Strait area highlights prehistoric human migrations from Siberia to northern North America

The patterns of prehistoric migrations across the Bering Land Bridge are far from being completely understood: there still exists a significant gap in our knowledge of the population history of former Beringia. Here, through comprehensive survey of mitochondrial DNA genomes retained in 'relic' populations, the Maritime Chukchi, Siberian Eskimos, and Commander Aleuts, we explore genetic contribution of prehistoric Siberians/Asians to northwestern Native Americans. Overall, 201 complete mitochondrial sequences (52 new and 149 published) were selected in the reconstruction of trees encompassing mtDNA lineages that are restricted to Coastal Chukotka and Alaska, the Canadian Arctic, Greenland, and the Aleutian chain. Phylogeography of the resulting mtDNA genomes (mitogenomes) considerably extends the range and intrinsic diversity of haplogroups (eg, A2a, A2b, D2a, and D4b1a2a1) that emerged and diversified in postglacial central Beringia, defining independent origins of Neo-Eskimos versus Paleo-Eskimos, Aleuts, and Tlingit (Na-Dene). Specifically, Neo-Eskimos, ancestral to modern Inuit, not only appear to be of the High Arctic origin but also to harbor Altai/Sayan-related ancestry. The occurrence of the haplogroup D2a1b haplotypes in Chukotka (Sireniki) introduces the possibility that the traces of Paleo-Eskimos have not been fully erased by spread of the Neo-Eskimos or their descendants. Our findings are consistent with the recurrent gene flow model of multiple streams of expansions to northern North America from northeastern Eurasia in late Pleistocene-early Holocene.

Genetic diversity in Puerto Rico and its implications for the peopling of the Island and the West Indies

Puerto Rico and the surrounding islands rest on the eastern fringe of the Caribbean's Greater Antilles, located less than 100 miles northwest of the Lesser Antilles. Puerto Ricans are genetic descendants of pre-Columbian peoples, as well as peoples of European and African descent through 500 years of migration to the island. To infer these patterns of pre-Columbian and historic peopling of the Caribbean, we characterized genetic diversity in 326 individuals from the southeastern region of Puerto Rico and the island municipality of Vieques. We sequenced the mitochondrial DNA (mtDNA) control region of all of the samples and the complete mitogenomes of 12 of them to infer their putative place of origin. In addition, we genotyped 121 male samples for 25 Y-chromosome single nucleotide polymorphism and 17 STR loci. Approximately 60% of the participants had indigenous mtDNA haplotypes (mostly from haplogroups A2 and C1), while 25% had African and 15% European haplotypes. Three A2 sublineages were unique to the Greater Antilles, one of which was similar to Mesoamerican types, while C1b haplogroups showed links to South America, suggesting that people reached the island from the two distinct continental source areas. However, none of the male participants had indigenous Y-chromosomes, with 85% of them instead being European/Mediterranean and 15% sub-Saharan African in origin. West Eurasian Y-chromosome short tandem repeat haplotypes were quite diverse and showed similarities to those observed in southern Europe, North Africa and the Middle East. These results attest to the distinct, yet equally complex, pasts for the male and female ancestors of modern day Puerto Ricans.

Complete mitogenome analysis of indigenous populations in Mexico: its relevance for the origin of Mesoamericans

Mesoamerica has an important role in the expansion of Paleoamericans as the route to South America. In this study, we determined complete mitogenome sequences of 113 unrelated individuals from two indigenous populations of Mesoamerica, Mazahua and Zapotec. All newly sequenced mitogenomes could be classified into haplogroups A2, B2, C1 and D1, but one sequence in Mazahua was D4h3a, a subclade of haplogroup D4. This haplogroup has been mostly found in South America along the Pacific coast. Haplogroup X2a was not found in either population. Genetic similarity obtained using phylogenetic tree construction and principal component analysis showed that these two populations are distantly related to each other. Actually, the Mazahua and the Zapotec shared no sequences (haplotypes) in common, while each also showed a number of unique subclades. Surprisingly, Zapotec formed a cluster with indigenous populations living in an area from central Mesoamerica to Central America. By contrast, the Mazahua formed a group with indigenous populations living in external areas, including southwestern North America and South America. This intriguing genetic relationship suggests the presence of two paleo-Mesoamerican groups, invoking a scenario in which one group had expanded into South America and the other resided in Mesoamerica.

The impact of modelling rate heterogeneity among sites on phylogenetic estimates of intraspecific evolutionary rates and timescales

Phylogenetic analyses of DNA sequence data can provide estimates of evolutionary rates and timescales. Nearly all phylogenetic methods rely on accurate models of nucleotide substitution. A key feature of molecular evolution is the heterogeneity of substitution rates among sites, which is often modelled using a discrete gamma distribution. A widely used derivative of this is the gamma-invariable mixture model, which assumes that a proportion of sites in the sequence are completely resistant to change, while substitution rates at the remaining sites are gamma-distributed. For data sampled at the intraspecific level, however, biological assumptions involved in the invariable-sites model are commonly violated. We examined the use of these models in analyses of five intraspecific data sets. We show that using 6-10 rate categories for the discrete gamma distribution of rates among sites is sufficient to provide a good approximation of the marginal likelihood. Increasing the number of gamma rate categories did not have a substantial effect on estimates of the substitution rate or coalescence time, unless rates varied strongly among sites in a non-gamma-distributed manner. The assumption of a proportion of invariable sites provided a better approximation of the asymptotic marginal likelihood when the number of gamma categories was small, but had minimal impact on estimates of rates and coalescence times. However, the estimated proportion of invariable sites was highly susceptible to changes in the number of gamma rate categories. The concurrent use of gamma and invariable-site models for intraspecific data is not biologically meaningful and has been challenged on statistical grounds; here we have found that the assumption of a proportion of invariable sites has no obvious impact on Bayesian estimates of rates and timescales from intraspecific data.

AmericaPlex26: a SNaPshot multiplex system for genotyping the main human mitochondrial founder lineages of the Americas

Phylogeographic studies have described a reduced genetic diversity in Native American populations, indicative of one or more bottleneck events during the peopling and prehistory of the Americas. Classical sequencing approaches targeting the mitochondrial diversity have reported the presence of five major haplogroups, namely A, B, C, D and X, whereas the advent of complete mitochondrial genome sequencing has recently refined the number of founder lineages within the given diversity to 15 sub-haplogroups. We developed and optimized a SNaPshot assay to study the mitochondrial diversity in pre-Columbian Native American populations by simultaneous typing of 26 single nucleotide polymorphisms (SNPs) characterising Native American sub-haplogroups. Our assay proved to be highly sensitive with respect to starting concentrations of target DNA and could be applied successfully to a range of ancient human skeletal material from South America from various time periods. The AmericaPlex26 is a powerful assay with enhanced phylogenetic resolution that allows time- and cost-efficient mitochondrial DNA sub-typing from valuable ancient specimens. It can be applied in addition or alternative to standard sequencing of the D-loop region in forensics, ancestry testing, and population studies, or where full-resolution mitochondrial genome sequencing is not feasible.

Characterization of mitochondrial haplogroups in a large population-based sample from the United States

Mitochondrial DNA (mtDNA) haplogroups are valuable for investigations in forensic science, molecular anthropology, and human genetics. In this study, we developed a custom panel of 61 mtDNA markers for high-throughput classification of European, African, and Native American/Asian mitochondrial haplogroup lineages. Using these mtDNA markers, we constructed a mitochondrial haplogroup classification tree and classified 18,832 participants from the National Health and Nutrition Examination Surveys (NHANES). To our knowledge, this is the largest study to date characterizing mitochondrial haplogroups in a population-based sample from the United States, and the first study characterizing mitochondrial haplogroup distributions in self-identified Mexican Americans separately from Hispanic Americans of other descent. We observed clear differences in the distribution of maternal genetic ancestry consistent with proposed admixture models for these subpopulations, underscoring the genetic heterogeneity of the United States Hispanic population. The mitochondrial haplogroup distributions in the other self-identified racial/ethnic groups within NHANES were largely comparable to previous studies. Mitochondrial haplogroup classification was highly concordant with self-identified race/ethnicity (SIRE) in non-Hispanic whites (94.8 %), but was considerably lower in admixed populations including non-Hispanic blacks (88.3 %), Mexican Americans (81.8 %), and other Hispanics (61.6 %), suggesting SIRE does not accurately reflect maternal genetic ancestry, particularly in populations with greater proportions of admixture. Thus, it is important to consider inconsistencies between SIRE and genetic ancestry when performing genetic association studies. The mitochondrial haplogroup data that we have generated, coupled with the epidemiologic variables in NHANES, is a valuable resource for future studies investigating the contribution of mtDNA variation to human health and disease.

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.

Populations of the Americas have a rich and heterogeneous genetic and cultural heritage that draws from a diversity of pre-Columbian Native American, European, and African populations. Characterizing this diversity facilitates the development of medical genetics research in diverse populations and the transfer of medical knowledge across populations. It also represents an opportunity to better understand the peopling of the Americas, from the crossing of Beringia to the post-Columbian era. Here, we take advantage sequencing of individuals of Colombian (CLM), Mexican (MXL), and Puerto Rican (PUR) origin by the 1000 Genomes project to improve our demographic models for the peopling of the Americas. The divergence among African, European, and Native American ancestors to these populations enables us to infer the continent of origin at each locus in the sampled genomes. The resulting patterns of ancestry suggest complex post-Columbian migration histories, starting later in CLM than in MXL and PUR. Whereas European ancestral segments show evidence of relatedness, a demographic model of synonymous variation suggests that the Native American Ancestors to MXL, PUR, and CLM panels split within a few hundred years over 12 thousand years ago. Together with early archeological sites in South America, these results support rapid divergence during the initial peopling of the Americas.

Complexities of CYP2D6 gene analysis and interpretation

Cytochrome P450 2D6 (CYP2D6) plays an important role in the metabolism and bioactivation of about 25% of clinically used drugs including many antidepressants, antipsychotics and opioids. CYP2D6 activity is highly variably ranging from no activity in so-called poor metabolizers to ultrarapid metabolism at the other end of the extreme of the activity distribution. A large portion of this variability can be explained by the highly polymorphic nature of the CYP2D6 gene locus for which > 100 variants and subvariants identified to date. Allele frequencies vary markedly between ethnic groups; some have exclusively or predominantly only been observed in certain populations. Pharmacogenetic testing holds the promise of individualizing drug therapy by identifying patients with CYP2D6 diplotypes that puts them at an increased risk of experiencing dose-related adverse events or therapeutic failure. Inferring a patient's CYP2D6 metabolic capacity, or phenotype, however, is a challenging task due to the complexity of the CYP2D6 gene locus. Allelic variation includes SNPs, small insertions and deletions, gene copy number variation and rearrangements with CYP2D7, a highly related non-functional gene. This review provides a summary of the intricacies of CYP2D6 variation and genotype analysis, knowledge that is invaluable for the translation of genotype into clinically useful information.

The first peopling of South America: new evidence from Y-chromosome haplogroup Q

Recent progress in the phylogenetic resolution of the Y-chromosome phylogeny permits the male demographic dynamics and migratory events that occurred in Central and Southern America after the initial human spread into the Americas to be investigated at the regional level. To delve further into this issue, we examined more than 400 Native American Y chromosomes (collected in the region ranging from Mexico to South America) belonging to haplogroup Q – virtually the only branch of the Y phylogeny observed in modern-day Amerindians of Central and South America – together with 27 from Mongolia and Kamchatka. Two main founding lineages, Q1a3a1a-M3 and Q1a3a1-L54(xM3), were detected along with novel sub-clades of younger age and more restricted geographic distributions. The first was also observed in Far East Asia while no Q1a3a1-L54(xM3) Y chromosome was found in Asia except the southern Siberian-specific sub-clade Q1a3a1c-L330. Our data not only confirm a southern Siberian origin of ancestral populations that gave rise to Paleo-Indians and the differentiation of both Native American Q founding lineages in Beringia, but support their concomitant arrival in Mesoamerica, where Mexico acted as recipient for the first wave of migration, followed by a rapid southward migration, along the Pacific coast, into the Andean region. Although Q1a3a1a-M3 and Q1a3a1-L54(xM3) display overlapping general distributions, they show different patterns of evolution in the Mexican plateau and the Andean area, which can be explained by local differentiations due to demographic events triggered by the introduction of agriculture and associated with the flourishing of the Great Empires.

Sequencing Y chromosomes resolves discrepancy in time to common ancestor of males versus females

The Y chromosome and the mitochondrial genome have been used to estimate when the common patrilineal and matrilineal ancestors of humans lived. We sequenced the genomes of 69 males from nine populations, including two in which we find basal branches of the Y-chromosome tree. We identify ancient phylogenetic structure within African haplogroups and resolve a long-standing ambiguity deep within the tree. Applying equivalent methodologies to the Y chromosome and the mitochondrial genome, we estimate the time to the most recent common ancestor (T(MRCA)) of the Y chromosome to be 120 to 156 thousand years and the mitochondrial genome T(MRCA) to be 99 to 148 thousand years. Our findings suggest that, contrary to previous claims, male lineages do not coalesce significantly more recently than female lineages.

Mitochondrial diversity in human head louse populations across the Americas

Anthropological studies suggest that the genetic makeup of human populations in the Americas is the result of diverse processes including the initial colonization of the continent by the first people plus post-1492 European migrations. Because of the recent nature of some of these events, understanding the geographical origin of American human diversity is challenging. However, human parasites have faster evolutionary rates and larger population sizes allowing them to maintain greater levels of genetic diversity than their hosts. Thus, we can use human parasites to provide insights into some aspects of human evolution that may be unclear from direct evidence. In this study, we analyzed mitochondrial DNA (mtDNA) sequences from 450 head lice in the Americas. Haplotypes clustered into two well-supported haplogroups, known as A and B. Haplogroup frequencies differ significantly among North, Central and South America. Within each haplogroup, we found evidence of demographic expansions around 16,000 and 20,000 years ago, which correspond broadly with those estimated for Native Americans. The parallel timing of demographic expansions of human lice and Native Americans plus the contrasting pattern between the distribution of haplogroups A and B through the Americas suggests that human lice can provide additional evidence about the human colonization of the New World.

The genome of Anopheles darlingi, the main neotropical malaria vector

Anopheles darlingi is the principal neotropical malaria vector, responsible for more than a million cases of malaria per year on the American continent. Anopheles darlingi diverged from the African and Asian malaria vectors ∼100 million years ago (mya) and successfully adapted to the New World environment. Here we present an annotated reference A. darlingi genome, sequenced from a wild population of males and females collected in the Brazilian Amazon. A total of 10 481 predicted protein-coding genes were annotated, 72% of which have their closest counterpart in Anopheles gambiae and 21% have highest similarity with other mosquito species. In spite of a long period of divergent evolution, conserved gene synteny was observed between A. darlingi and A. gambiae. More than 10 million single nucleotide polymorphisms and short indels with potential use as genetic markers were identified. Transposable elements correspond to 2.3% of the A. darlingi genome. Genes associated with hematophagy, immunity and insecticide resistance, directly involved in vector–human and vector–parasite interactions, were identified and discussed. This study represents the first effort to sequence the genome of a neotropical malaria vector, and opens a new window through which we can contemplate the evolutionary history of anopheline mosquitoes. It also provides valuable information that may lead to novel strategies to reduce malaria transmission on the South American continent. The A. darlingi genome is accessible at www.labinfo.lncc.br/index.php/anopheles-darlingi.

The genetic legacy of the pre-colonial period in contemporary Bolivians

Only a few genetic studies have been carried out to date in Bolivia. However, some of the most important (pre)historical enclaves of South America were located in these territories. Thus, the (sub)-Andean region of Bolivia was part of the Inca Empire, the largest state in Pre-Columbian America. We have genotyped the first hypervariable region (HVS-I) of 720 samples representing the main regions in Bolivia, and these data have been analyzed in the context of other pan-American samples (>19,000 HVS-I mtDNAs). Entire mtDNA genome sequencing was also undertaken on selected Native American lineages. Additionally, a panel of 46 Ancestry Informative Markers (AIMs) was genotyped in a sub-set of samples. The vast majority of the Bolivian mtDNAs (98.4%) were found to belong to the main Native American haplogroups (A: 14.3%, B: 52.6%, C: 21.9%, D: 9.6%), with little indication of sub-Saharan and/or European lineages; however, marked patterns of haplogroup frequencies between main regions exist (e.g. haplogroup B: Andean [71%], Sub-Andean [61%], Llanos [32%]). Analysis of entire genomes unraveled the phylogenetic characteristics of three Native haplogroups: the pan-American haplogroup B2b (originated ∼21.4 thousand years ago [kya]), A2ah (∼5.2 kya), and B2o (∼2.6 kya). The data suggest that B2b could have arisen in North California (an origin even in the north most region of the American continent cannot be disregarded), moved southward following the Pacific coastline and crossed Meso-America. Then, it most likely spread into South America following two routes: the Pacific path towards Peru and Bolivia (arriving here at about ∼15.2 kya), and the Amazonian route of Venezuela and Brazil southwards. In contrast to the mtDNA, Ancestry Informative Markers (AIMs) reveal a higher (although geographically variable) European introgression in Bolivians (25%). Bolivia shows a decreasing autosomal molecular diversity pattern along the longitudinal axis, from the Altiplano to the lowlands. Both autosomes and mtDNA revealed a low impact (1-2%) of a sub-Saharan component in Bolivians.

DNA analysis of an early modern human from Tianyuan Cave, China

Hominins with morphology similar to present-day humans appear in the fossil record across Eurasia between 40,000 and 50,000 y ago. The genetic relationships between these early modern humans and present-day human populations have not been established. We have extracted DNA from a 40,000-y-old anatomically modern human from Tianyuan Cave outside Beijing, China. Using a highly scalable hybridization enrichment strategy, we determined the DNA sequences of the mitochondrial genome, the entire nonrepetitive portion of chromosome 21 (∼30 Mbp), and over 3,000 polymorphic sites across the nuclear genome of this individual. The nuclear DNA sequences determined from this early modern human reveal that the Tianyuan individual derived from a population that was ancestral to many present-day Asians and Native Americans but postdated the divergence of Asians from Europeans. They also show that this individual carried proportions of DNA variants derived from archaic humans similar to present-day people in mainland Asia.

On the origins, rapid expansion and genetic diversity of Native Americans from hunting-gatherers to agriculturalists

Given the importance of Y-chromosome haplogroup Q to better understand the source populations of contemporary Native Americans, we studied 8 biallelic and 17 microsatellite polymorphisms on the background of 128 Q Y-chromosomes from geographically targeted populations. The populations examined in this study include three from the Tuva Republic in Central Asia (Bai-Tai, Kungurtug, and Toora-Hem, n = 146), two from the northeastern tip of Siberia (New Chaplino and Chukchi, n = 32), and two from Mesoamerica (Mayans from Yucatan, Mexico n = 72, and Mayans from the Guatemalan Highlands, n = 43). We also see evidence of a dramatic Mesoamerican post-migration population growth in the ubiquitous and diverse Y-STR profiles of the Mayan and other Mesoamerican populations. In the case of the Mayans, this demographic growth was most likely fueled by the agricultural- and trade-based subsistence adopted during the Pre-Classic, Classic and Post-Classic periods of their empire. The limited diversity levels observed in the Altaian and Tuvinian regions of Central Asia, the lowest of all populations examined, may be the consequence of bottleneck events fostered by the spatial isolation and low effective population size characteristic of a nomadic lifestyle. Furthermore, our data illustrate how a sociocultural characteristic such as mode of subsistence may be of impact on the genetic structure of populations. We analyzed our genetic data using Multidimensional Scaling Analysis of populations, Principal Component Analysis of individuals, Median-joining networks of M242, M346, L54, and M3 individuals, age estimations based on microsatellite variation utilizing genealogical and evolutionary mutation rates/generation times and estimation of Y- STR average gene diversity indices.

Arrival of Paleo-Indians to the southern cone of South America: new clues from mitogenomes

With analyses of entire mitogenomes, studies of Native American mitochondrial DNA (mtDNA) variation have entered the final phase of phylogenetic refinement: the dissection of the founding haplogroups into clades that arose in America during and after human arrival and spread. Ages and geographic distributions of these clades could provide novel clues on the colonization processes of the different regions of the double continent. As for the Southern Cone of South America, this approach has recently allowed the identification of two local clades (D1g and D1j) whose age estimates agree with the dating of the earliest archaeological sites in South America, indicating that Paleo-Indians might have reached that region from Beringia in less than 2000 years. In this study, we sequenced 46 mitogenomes belonging to two additional clades, termed B2i2 (former B2l) and C1b13, which were recently identified on the basis of mtDNA control-region data and whose geographical distributions appear to be restricted to Chile and Argentina. We confirm that their mutational motifs most likely arose in the Southern Cone region. However, the age estimate for B2i2 and C1b13 (11–13,000 years) appears to be younger than those of other local clades. The difference could reflect the different evolutionary origins of the distinct South American-specific sub-haplogroups, with some being already present, at different times and locations, at the very front of the expansion wave in South America, and others originating later in situ, when the tribalization process had already begun. A delayed origin of a few thousand years in one of the locally derived populations, possibly in the central part of Chile, would have limited the geographical and ethnic diffusion of B2i2 and explain the present-day occurrence that appears to be mainly confined to the Tehuelche and Araucanian-speaking groups.

MtDNA analysis of global populations support that major population expansions began before Neolithic Time

Agriculture resulted in extensive population growths and human activities. However, whether major human expansions started after Neolithic Time still remained controversial. With the benefit of 1000 Genome Project, we were able to analyze a total of 910 samples from 11 populations in Africa, Europe and Americas. From these random samples, we identified the expansion lineages and reconstructed the historical demographic variations. In all the three continents, we found that most major lineage expansions (11 out of 15 star lineages in Africa, all autochthonous lineages in Europe and America) coalesced before the first appearance of agriculture. Furthermore, major population expansions were estimated after Last Glacial Maximum but before Neolithic Time, also corresponding to the result of major lineage expansions. Considering results in current and previous study, global mtDNA evidence showed that rising temperature after Last Glacial Maximum offered amiable environments and might be the most important factor for prehistorical human expansions.

Reconstructing the history of Mesoamerican populations through the study of the mitochondrial DNA control region

The study of genetic information can reveal a reconstruction of human population’s history. We sequenced the entire mtDNA control region (positions 16.024 to 576 following Cambridge Reference Sequence, CRS) of 605 individuals from seven Mesoamerican indigenous groups and one Aridoamerican from the Greater Southwest previously defined, all of them in present Mexico. Samples were collected directly from the indigenous populations, the application of an individual survey made it possible to remove related or with other origins samples. Diversity indices and demographic estimates were calculated. Also AMOVAs were calculated according to different criteria. An MDS plot, based on FST distances, was also built. We carried out the construction of individual networks for the four Amerindian haplogroups detected. Finally, barrier software was applied to detect genetic boundaries among populations. The results suggest: a common origin of the indigenous groups; a small degree of European admixture; and inter-ethnic gene flow. The process of Mesoamerica’s human settlement took place quickly influenced by the region’s orography, which development of genetic and cultural differences facilitated. We find the existence of genetic structure is related to the region’s geography, rather than to cultural parameters, such as language. The human population gradually became fragmented, though they remained relatively isolated, and differentiated due to small population sizes and different survival strategies. Genetic differences were detected between Aridoamerica and Mesoamerica, which can be subdivided into “East”, “Center”, “West” and “Southeast”. The fragmentation process occurred mainly during the Mesoamerican Pre-Classic period, with the Otomí being one of the oldest groups. With an increased number of populations studied adding previously published data, there is no change in the conclusions, although significant genetic heterogeneity can be detected in Pima and Huichol groups. This result may be explained because populations historically assigned as belonging to the same group were, in fact, different indigenous populations.

Miocene and Pliocene dominated diversification of the lichen-forming fungal genus Melanohalea (Parmeliaceae, Ascomycota) and Pleistocene population expansions

Background

Factors promoting diversification in lichen symbioses remain largely unexplored. While Pleistocene events have been important for driving diversification and affecting distributions in many groups, recent estimates suggest that major radiations within some genera in the largest clade of macrolichens (Parmeliaceae, Ascomycota) vastly predate the Pleistocene. To better understand the temporal placement and sequence of diversification events in lichens, we estimated divergence times in a common lichen-forming fungal genus, Melanohalea, in the Northern Hemisphere. Divergence times were estimated using both concatenated gene tree and coalescent-based multilocus species tree approaches to assess the temporal context of major radiation events within Melanohalea. In order to complement our understanding of processes impacting genetic differentiation, we also evaluated the effects of Pleistocene glacial cycles on population demographics of distinct Melanohalea lineages, differing in reproductive strategies.

Results

We found that divergence estimates, from both concatenated gene tree and coalescent-based multilocus species tree approaches, suggest that diversification within Melanohalea occurred predominantly during the Miocene and Pliocene, although estimated of divergence times differed by up to 8.3 million years between the two methods. These results indicate that, in some cases, taxonomically diagnostic characters may be maintained among divergent lineages for millions of years. In other cases, similar phenotypic characters among non-sister taxa, including reproductive strategies, suggest the potential for convergent evolution due to similar selective pressures among distinct lineages. Our analyses provide evidence of population expansions predating the last glacial maximum in the sampled lineages. These results suggest that Pleistocene glaciations were not inherently unfavorable or restrictive for some Melanohalea species, albeit with apparently different demographic histories between sexually and vegetatively reproducing lineages.

Conclusions

Our results contribute to the understanding of how major changes during the Miocene and Pliocene have been important in promoting diversification within common lichen-forming fungi in the northern Hemisphere. Additionally, we provide evidence that glacial oscillations have influenced current population structure of broadly distributed lichenized fungal species throughout the Holarctic.