Improved calibration of the human mitochondrial clock using ancient genomes

Convergent Mutations and Single Nucleotide Variants in Mitochondrial Genomes of Modern Humans and Neanderthals

The genetic contributions of Neanderthals to the modern human genome have been evidenced by the comparison of present-day human genomes with paleogenomes. Neanderthal signatures in extant human genomes are attributed to intercrosses between Neanderthals and archaic anatomically modern humans (AMHs). Although Neanderthal signatures are well documented in the nuclear genome, it has been proposed that there is no contribution of Neanderthal mitochondrial DNA to contemporary human genomes. Here we show that modern human mitochondrial genomes contain 66 potential Neanderthal signatures, or Neanderthal single nucleotide variants (N-SNVs), of which 36 lie in coding regions and 7 result in nonsynonymous changes. Seven N-SNVs are associated with traits such as cycling vomiting syndrome, Alzheimer's disease and Parkinson's disease, and two N-SNVs are associated with intelligence quotient. Based on recombination tests, principal component analysis (PCA) and the complete absence of these N-SNVs in 41 archaic AMH mitogenomes, we conclude that convergent evolution, and not recombination, explains the presence of N-SNVs in present-day human mitogenomes.

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 (14 C) 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 (14 C) 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.

Dating ancient splits in phylogenetic trees, with application to the human-Neanderthal split

BACKGROUND

We tackle the problem of estimating species TMRCAs (Time to Most Recent Common Ancestor), given a genome sequence from each species and a large known phylogenetic tree with a known structure (typically from one of the species). The number of transitions at each site from the first sequence to the other is assumed to be Poisson distributed, and only the parity of the number of transitions is observed. The detailed phylogenetic tree contains information about the transition rates in each site. We use this formulation to develop and analyze multiple estimators of the species' TMRCA. To test our methods, we use mtDNA substitution statistics from the well-established Phylotree as a baseline for data simulation such that the substitution rate per site mimics the real-world observed rates.

RESULTS

We evaluate our methods using simulated data and compare them to the Bayesian optimizing software BEAST2, showing that our proposed estimators are accurate for a wide range of TMRCAs and significantly outperform BEAST2. We then apply the proposed estimators on Neanderthal, Denisovan, and Chimpanzee mtDNA genomes to better estimate their TMRCA with modern humans and find that their TMRCA is substantially later, compared to values cited recently in the literature.

CONCLUSIONS

Our methods utilize the transition statistics from the entire known human mtDNA phylogenetic tree (Phylotree), eliminating the requirement to reconstruct a tree encompassing the specific sequences of interest. Moreover, they demonstrate notable improvement in both running speed and accuracy compared to BEAST2, particularly for earlier TMRCAs like the human-Chimpanzee split. Our results date the human - Neanderthal TMRCA to be [Formula: see text] years ago, considerably later than values cited in other recent studies.

Relationships of Late Pleistocene giant deer as revealed by Sinomegaceros mitogenomes from East Asia

The giant deer, widespread in northern Eurasia during the Late Pleistocene, have been classified as western Megaloceros and eastern Sinomegaceros through morphological studies. While Megaloceros's evolutionary history has been unveiled through mitogenomes, Sinomegaceros remains molecularly unexplored. Herein, we generated mitogenomes of giant deer from East Asia. We find that, in contrast to the morphological differences between Megaloceros and Sinomegaceros, they are mixed in the mitochondrial phylogeny, and Siberian specimens suggest a range contact or overlap between these two groups. Meanwhile, one deep divergent clade and another surviving until 20.1 thousand years ago (ka) were detected in northeastern China, the latter implying this area as a potential refugium during the Last Glacial Maximum (LGM). Moreover, stable isotope analyses indicate correlations between climate-introduced vegetation changes and giant deer extinction. Our study demonstrates the genetic relationship between eastern and western giant deer and explores the promoters of their extirpation in northern East Asia.

North and East African mitochondrial genetic variation needs further characterization towards precision medicine

INTRODUCTION

Mitochondria are maternally inherited cell organelles with their own genome, and perform various functions in eukaryotic cells such as energy production and cellular homeostasis. Due to their inheritance and manifold biological roles in health and disease, mitochondrial genetics serves a dual purpose of tracing the history as well as disease susceptibility of human populations across the globe. This work requires a comprehensive catalogue of commonly observed genetic variations in the mitochondrial DNAs for all regions throughout the world. So far, however, certain regions, such as North and East Africa have been understudied.

OBJECTIVES

To address this shortcoming, we have created the most comprehensive quality-controlled North and East African mitochondrial data set to date and use it for characterizing mitochondrial genetic variation in this region.

METHODS

We compiled 11 published cohorts with novel data for mitochondrial genomes from 159 Sudanese individuals. We combined these 641 mitochondrial sequences with sequences from the 1000 Genomes (n = 2504) and the Human Genome Diversity Project (n = 828) and used the tool haplocheck for extensive quality control and detection of in-sample contamination, as well as Nanopore long read sequencing for haplogroup validation of 18 samples.

RESULTS

Using a subset of high-coverage mitochondrial sequences, we predict 15 potentially novel haplogroups in North and East African subjects and observe likely phylogenetic deviations from the established PhyloTree reference for haplogroups L0a1 and L2a1.

CONCLUSION

Our findings demonstrate common hitherto unexplored variants in mitochondrial genomes of North and East Africa that lead to novel phylogenetic relationships between haplogroups present in these regions. These observations call for further in-depth population genetic studies in that region to enable the prospective use of mitochondrial genetic variation for precision medicine.

Unraveling the mitochondrial phylogenetic landscape of Thailand reveals complex admixture and demographic dynamics

The evolutionary dynamics of mitochondrial DNA within the Thai population were comprehensively explored with a specific focus on the influence of South Asian admixture. A total of 166 samples were collected through randomized sampling, ensuring a diverse representation. Our findings unveil substantial genetic and haplogroup diversity within the Thai population. We have identified 164 haplotypes categorized into 97 haplogroups, with a notable inclusion of 20 novel haplogroups. The distribution of haplogroups exhibited variations across different populations and countries. The central Thai population displayed a high diversity of haplogroups from both the M and N clades. Maternal lineage affinities were discerned between several Mainland Southeast Asia (MSEA) and South Asian populations, implying ancestral genetic connections and a substantial influence of South Asian women in establishing these relationships. f4-statistics indicates the presence of a Tibeto-Burman genetic component within the Mon population from Thailand. New findings demonstrate two phases of population expansion occurring 22,000-26,000 and 2500-3800 years ago, coinciding with the Last Glacial Maximum, and Neolithic demographic transition, respectively. This research significantly enhances our understanding of the maternal genetic history of Thailand and MSEA, emphasizing the influence of South Asian admixture. Moreover, it underscores the critical role of prior information, such as mutation rates, within the Bayesian framework for accurate estimation of coalescence times and inferring demographic history.

Mitogenome analysis reveals predominantly ancient Yellow River origin of population inhabiting Datong agro-pastoral ecotone along Great Wall

The Datong Basin was an important arena for population movement and admixture between the Yellow River Valley and Eastern Steppe. In historical materials, the region was often the setting for a tug-of-war between Han farmers and non-Han nomads. The genetic makeup and population history of this Datong population has, however, remained uncertain. In this study, we analysed 289 mitogenomes from Datong individuals. Our primary findings were: (1) population summary statistics analysis revealed a high level of genetic diversity and strong signals of population expansion in the Datong population; (2) inter-population comparisons (PCA and Fst heatmap) exhibited a close clustering between the Datong population and Northern Han, especially northern frontier groups, such as the Inner Mongolia Han, Heilongjiang Han, Liaoning Han and Tianjin Han; (3) phylogeographic analysis of complete mitogenomes revealed the presence of different components in the maternal gene pools of Datong population-the northern East Asian component was dominant (66.44%), whereas the southern East Asians were the second largest component with 31.49%. We also observed a much reduced west Eurasian (2.07%) component; (4) direct comparisons with ancient groups showed closer relationship between Datong and Yellow River farmers than Eastern Steppe nomads. Despite, therefore, centuries of Eastern Steppe nomadic control over the Datong area, Yellow River farmers had a much more significant impact on the Datong population.

Mitogenomics of the Koryaks and Evens of the northern coast of the Sea of Okhotsk

Due to the geographical proximity of the northern coast of the Sea of Okhotsk and Kamchatka Peninsula to the Beringia, the indigenous populations of these territories are of great interest for elucidating the human settlement history of northern Asia and America. Meanwhile, there is a clear shortage of genetic studies of the indigenous populations of the northern coast of the Sea of Okhotsk. Here, in order to examine their fine-scale matrilineal genetic structure, ancestry and relationships with neighboring populations, we analyzed 203 complete mitogenomes (174 of which are new) from population samples of the Koryaks and Evens of the northern coast of the Sea of Okhotsk and the Chukchi of the extreme northeast Asia. The patterns observed underscore the reduced level of genetic diversity found in the Koryak, Even, and Chukchi populations, which, along with the high degree of interpopulation differentiation, may be the result of genetic drift. Our phylogeographic analysis reveals common Paleo-Asiatic ancestry for 51.1% of the Koryaks and 17.8% of the Evens. About third of the mitogenomes found in the Koryaks and Evens might be considered as ethno-specific, as these are virtually absent elsewhere in North, Central and East Asia. Coalescence ages of most of these lineages coincide well with the emergence and development of the Tokarev and Old Koryak archaeological cultures associated with the formation of the Koryaks, as well as with the period of separation and split of the North Tungusic groups migrated northwards from the Lake Baikal or the Amur River area.

Mitochondrial DNA in Human Diversity and Health: From the Golden Age to the Omics Era

Mitochondrial DNA (mtDNA) is a small fraction of our hereditary material. However, this molecule has had an overwhelming presence in scientific research for decades until the arrival of high-throughput studies. Several appealing properties justify the application of mtDNA to understand how human populations are-from a genetic perspective-and how individuals exhibit phenotypes of biomedical importance. Here, I review the basics of mitochondrial studies with a focus on the dawn of the field, analysis methods and the connection between two sides of mitochondrial genetics: anthropological and biomedical. The particularities of mtDNA, with respect to inheritance pattern, evolutionary rate and dependence on the nuclear genome, explain the challenges of associating mtDNA composition and diseases. Finally, I consider the relevance of this single locus in the context of omics research. The present work may serve as a tribute to a tool that has provided important insights into the past and present of humankind.

Herbarium specimen sequencing allows precise dating of Xanthomonas citri pv. citri diversification history

Herbarium collections are an important source of dated, identified and preserved DNA, whose use in comparative genomics and phylogeography can shed light on the emergence and evolutionary history of plant pathogens. Here, we reconstruct 13 historical genomes of the bacterial crop pathogen Xanthomonas citri pv. citri (Xci) from infected Citrus herbarium specimens. Following authentication based on ancient DNA damage patterns, we compare them with a large set of modern genomes to estimate their phylogenetic relationships, pathogenicity-associated gene content and several evolutionary parameters. Our results indicate that Xci originated in Southern Asia ~11,500 years ago (perhaps in relation to Neolithic climate change and the development of agriculture) and diversified during the beginning of the 13th century, after Citrus diversification and before spreading to the rest of the world (probably via human-driven expansion of citriculture through early East-West trade and colonization).

Ancient mitogenomes from the Southern Pampas of Argentina reflect local differentiation and limited extra-regional linkages after rapid initial colonization

OBJECTIVE

This study aims to contribute to the recovery of Indigenous evolutionary history in the Southern Pampas region of Argentina through an analysis of ancient complete mitochondrial genomes.

MATERIALS AND METHODS

We generated DNA data for nine complete mitogenomes from the Southern Pampas, dated to between 2531 and 723 cal BP. In combination with previously published ancient mitogenomes from the region and from throughout South America, we documented instances of extra-regional lineage-sharing, and estimated coalescent ages for local lineages using a Bayesian method with tip calibrations in a phylogenetic analysis.

RESULTS

We identified a novel mitochondrial haplogroup, B2b16, and two recently defined haplogroups, A2ay and B2ak1, as well as three local haplotypes within founder haplogroups C1b and C1d. We detected lineage-sharing with ancient and contemporary individuals from Central Argentina, but not with ancient or contemporary samples from North Patagonian or Littoral regions of Argentina, despite archeological evidence of cultural interactions with the latter regions. The estimated coalescent age of these shared lineages is ~10,000 years BP.

DISCUSSION

The history of the human populations in the Southern Pampas is temporally deep, exhibiting long-term continuity of mitogenome lineages. Additionally, the identification of highly localized mtDNA clades accords with a model of relatively rapid initial colonization of South America by Indigenous communities, followed by more local patterns of limited gene flow and genetic drift in various South American regions, including the Pampas.

Museum specimens of a landlocked pinniped reveal recent loss of genetic diversity and unexpected population connections

The Saimaa ringed seal (Pusa hispida saimensis) is endemic to Lake Saimaa in Finland. The subspecies is thought to have originated when parts of the ringed seal population of the Baltic region were trapped in lakes emerging due to postglacial bedrock rebound around 9000 years ago. During the 20th century, the population experienced a drastic human-induced bottleneck. Today encompassing a little over 400 seals with extremely low genetic diversity, it is classified as endangered. We sequenced sections of the mitochondrial control region from 60 up to 125-years-old museum specimens of the Saimaa ringed seal. The generated dataset was combined with publicly available sequences. We studied how genetic variation has changed through time in this subspecies and how it is phylogenetically related to other ringed seal populations from the Baltic Sea, Lake Ladoga, North America, Svalbard, and the White Sea. We observed temporal fluctuations in haplotype frequencies and loss of haplotypes accompanied by a recent reduction in female effective population size. In apparent contrast with the traditionally held view of the Baltic origin of the population, the Saimaa ringed seal mtDNA variation also shows affinities to North American ringed seals. Our results suggest that the Saimaa ringed seal has experienced recent genetic drift associated with small population size. The results further suggest that extant Baltic ringed seal is not representative of the ancestral population of the Saimaa ringed seal, which calls for re-evaluation of the deep history of this subspecies.

Kazak mitochondrial genomes provide insights into the human population history of Central Eurasia

As a historical nomadic group in Central Asia, Kazaks have mainly inhabited the steppe zone from the Altay Mountains in the East to the Caspian Sea in the West. Fine scale characterization of the genetic profile and population structure of Kazaks would be invaluable for understanding their population history and modeling prehistoric human expansions across the Eurasian steppes. With this mind, we characterized the maternal lineages of 200 Kazaks from Jetisuu at mitochondrial genome level. Our results reveal that Jetisuu Kazaks have unique mtDNA haplotypes including those belonging to the basal branches of both West Eurasian (R0, H, HV) and East Eurasian (A, B, C, D) lineages. The great diversity observed in their maternal lineages may reflect pivotal geographic location of Kazaks in Eurasia and implies a complex history for this population. Comparative analyses of mitochondrial genomes of human populations in Central Eurasia reveal a common maternal genetic ancestry for Turko-Mongolian speakers and their expansion being responsible for the presence of East Eurasian maternal lineages in Central Eurasia. Our analyses further indicate maternal genetic affinity between the Sherpas from the Tibetan Plateau with the Turko-Mongolian speakers.

Genomic ancestry, diet and microbiomes of Upper Palaeolithic hunter-gatherers from San Teodoro cave

Recent improvements in the analysis of ancient biomolecules from human remains and associated dental calculus have provided new insights into the prehistoric diet and genetic diversity of our species. Here we present a multi-omics study, integrating metagenomic and proteomic analyses of dental calculus, and human ancient DNA analysis of the petrous bones of two post-Last Glacial Maximum (LGM) individuals from San Teodoro cave (Italy), to reconstruct their lifestyle and the post-LGM resettlement of Europe. Our analyses show genetic homogeneity in Sicily during the Palaeolithic, representing a hitherto unknown Italian genetic lineage within the previously identified Villabruna cluster. We argue that this lineage took refuge in Italy during the LGM, followed by a subsequent spread to central-western Europe. Analysis of dental calculus showed a diet rich in animal proteins which is also reflected on the oral microbiome composition. Our results demonstrate the power of this approach in the study of prehistoric humans and will enable future research to reach a more holistic understanding of the population dynamics and ecology.

Uniparental Genetic Analyses Reveal Multi-Ethnic Background of Dunhuang Foyemiaowan Population (220–907 CE) With Typical Han Chinese Archaological Culture

The relationship between archeological culture and ethnicity is invariably complex. This is especially the case for periods of national division and rapid inter-ethnic exchange, such as China’s Sixteen Kingdoms (304–439 CE) and Northern and Southern Dynasties (420–589 CE). Going by tomb shape and grave goods, the Foyemiaowan cemetery at Dunhuang exhibits a typical third–tenth century Han style. Despite this, the ethnic makeup of the Foyemiaowan population has remained unclear. We therefore analyzed 485 Y-chromosomal SNPs and entire mitochondrial genomes of 34 Foyemiaowan samples. Our study yielded the following discoveries: (1) principal component analysis revealed that the Foyemiaowan population was closely clustered with Tibeto-Burman populations on the paternal side and close to Mongolic-speaking populations on the maternal side; (2) lineage comparisons at the individual level showed that the Foyemiaowan population consisted of primarily Tibeto-Burman and Han Chinese related lineages (Oα-M117, 25%;Oβ-F46, 18.75%), partially Altaic speaking North Eurasian lineages (N-F1206, 18.75%) and a slight admixture of southern East Asian lineages (O1b1a2-Page59, 6.25%; O1b1a1-PK4, 3.13%). Similarly, the maternal gene pool of Foyemiaowan contained northern East Asian (A, 4.17%; CZ, 16.67%; D, 20.83%; G, 4.17%; M9, 4.17%), southern East Asian (B, 12.51%; F, 20.83%) and western Eurasian (H, 4.17%; J, 4.17%) related lineages; (3) we discovered a relatively high genetic diversity among the Foyemiaowan population (0.891) in our ancient reference populations, indicating a complex history of population admixture. Archeological findings, stable isotope analysis and historical documents further corroborated our results. Although in this period China’s central government had relinquished control of the Hexi Corridor and regional non-Han regimes became the dominant regional power, Foyemiaowan’s inhabitants remained strongly influenced by Han culture.

African mitochondrial haplogroup L7: a 100,000-year-old maternal human lineage discovered through reassessment and new sequencing

Archaeological and genomic evidence suggest that modern Homo sapiens have roamed the planet for some 300–500 thousand years. In contrast, global human mitochondrial (mtDNA) diversity coalesces to one African female ancestor (“Mitochondrial Eve”) some 145 thousand years ago, owing to the ¼ gene pool size of our matrilineally inherited haploid genome. Therefore, most of human prehistory was spent in Africa where early ancestors of Southern African Khoisan and Central African rainforest hunter-gatherers (RFHGs) segregated into smaller groups. Their subdivisions followed climatic oscillations, new modes of subsistence, local adaptations, and cultural-linguistic differences, all prior to their exodus out of Africa. Seven African mtDNA haplogroups (L0–L6) traditionally captured this ancient structure—these L haplogroups have formed the backbone of the mtDNA tree for nearly two decades. Here we describe L7, an eighth haplogroup that we estimate to be ~ 100 thousand years old and which has been previously misclassified in the literature. In addition, L7 has a phylogenetic sublineage L7a*, the oldest singleton branch in the human mtDNA tree (~ 80 thousand years). We found that L7 and its sister group L5 are both low-frequency relics centered around East Africa, but in different populations (L7: Sandawe; L5: Mbuti). Although three small subclades of African foragers hint at the population origins of L5'7, the majority of subclades are divided into Afro-Asiatic and eastern Bantu groups, indicative of more recent admixture. A regular re-estimation of the entire mtDNA haplotype tree is needed to ensure correct cladistic placement of new samples in the future.

Variation in the substitution rates among the human mitochondrial haplogroup U sublineages

Resolving the absolute timescale of phylogenetic trees stipulates reliable estimates for the rate of DNA sequence evolution. For this end, various calibration methods have been developed and studied intensively. Intraspecific rate variation among distinct genetic lineages, however, has gained less attention. Here, we have assessed lineage-specific molecular rates of human mitochondrial DNA (mtDNA) by performing tip-calibrated Bayesian phylogenetic analyses. Tip-calibration, as opposed to traditional nodal time stamps from dated fossil evidence or geological events, is based on sample ages and becoming ever more feasible as ancient DNA data from radiocarbon-dated samples accumulate. We focus on subhaplogroups U2, U4, U5a, and U5b, the data including ancient mtDNA genomes from ¹⁴C-dated samples (n = 234), contemporary genomes (n = 301), and two outgroup sequences from haplogroup R. The obtained molecular rates depended on the data sets (with or without contemporary sequences), suggesting time-dependency. More notable was the rate variation between haplogroups: U4 and U5a stand out having a substantially higher rate than U5b. This is also reflected in the divergence times obtained (U5a: 17,700 years and U5b: 29,700 years), a disparity not reported previously. After ruling out various alternative causes (e.g., selection, sampling, and sequence quality), we propose that the substitution rates have been influenced by demographic histories, widely different among populations where U4/U5a or U5b are frequent. As with the Y-chromosomal subhaplogroup R1b, the mitochondrial U4 and U5a have been associated with remarkable range extensions of the Yamnaya culture in the Bronze Age.

Anatomical and molecular analyses of the deltoid muscle in chimpanzees (Pan troglodytes) and modern humans (Homo sapiens): Similarities and differences due to the uses of the upper extremity

In the deltoid muscles of Pan troglodytes and Homo sapiens, we have analyzed the muscle architecture and the expression of the myosin heavy chain (MHC) isoforms. Our aim was to identify differences between the two species that could be related to their different uses of the upper limb. The deltoid muscle of six adult Pan troglodytes and six adult Homo sapiens were dissected. The muscle fascicle length (MFL) and the physiological cross-sectional area (PCSA) of each muscle were calculated in absolute and normalized values. The expression pattern of the MHC-I, MHC-IIa and MHC-IIx isoforms was analyzed in the same muscles by real-time polymerase chain reaction. Only the acromial deltoid (AD) presented significant architectural differences between the two species, with higher MFL values in humans and higher PCSA values in chimpanzees. No significant differences in the expression pattern of the MHC isoforms were identified. The higher PCSA values in the AD of Pan troglodytes indicate a greater capacity of force generation in chimpanzees than in humans, which may be related to a greater use of the upper limb in locomotion, specifically in arboreal locomotion like vertical climbing. The functional differences between chimpanzees and humans in the deltoid muscle are more related to muscle architecture than to a differential expression of MHC isoforms.

Pedigree derived mutation rate across the entire mitochondrial genome of the Norfolk Island population

Estimates of mutation rates for various regions of the human mitochondrial genome (mtGenome) vary widely, depending on whether they are inferred using a phylogenetic approach or obtained directly from pedigrees. Traditionally, only the control region, or small portions of the coding region have been targeted for analysis due to the cost and effort required to produce whole mtGenome Sanger profiles. Here, we report one of the first pedigree derived mutation rates for the entire human mtGenome. The entire mtGenome from 225 individuals originating from Norfolk Island was analysed to estimate the pedigree derived mutation rate and compared against published mutation rates. These individuals were from 45 maternal lineages spanning 345 generational events. Mutation rates for various portions of the mtGenome were calculated. Nine mutations (including two transitions and seven cases of heteroplasmy) were observed, resulting in a rate of 0.058 mutations/site/million years (95% CI 0.031-0.108). These mutation rates are approximately 16 times higher than estimates derived from phylogenetic analysis with heteroplasmy detected in 13 samples (n = 225, 5.8% individuals). Providing one of the first pedigree derived estimates for the entire mtGenome, this study provides a better understanding of human mtGenome evolution and has relevance to many research fields, including medicine, anthropology and forensics.

Simulated patterns of mitochondrial diversity are consistent with partial population turnover in Bronze Age Central Europe

OBJECTIVES

The analysis of ancient mitochondrial DNA from osteological remains has challenged previous conclusions drawn from the analysis of mitochondrial DNA from present populations, notably by revealing an absence of genetic continuity between the Neolithic and modern populations in Central Europe. Our study investigates how to reconcile these contradictions at the mitochondrial level using a modeling approach.

MATERIALS AND METHODS

We used a spatially explicit computational framework to simulate ancient and modern DNA sequences under various evolutionary scenarios of post Neolithic demographic events and compared the genetic diversity of the simulated and observed mitochondrial sequences. We investigated which-if any-scenarios were able to reproduce statistics of genetic diversity similar to those observed, with a focus on the haplogroup N1a, associated with the spread of early Neolithic farmers.

RESULTS

Demographic fluctuations during the Neolithic transition or subsequent demographic collapses after this period, that is, due to epidemics such as plague, are not sufficient to explain the signal of population discontinuity detected on the mitochondrial DNA in Central Europe. Only a scenario involving a substantial genetic input due to the arrival of migrants after the Neolithic transition, possibly during the Bronze Age, is compatible with observed patterns of genetic diversity.

DISCUSSION

Our results corroborate paleogenomic studies, since out of the alternative hypotheses tested, the best one that was able to recover observed patterns of mitochondrial diversity in modern and ancient Central European populations was one were immigration of populations from the Pontic steppes during the Bronze Age was explicitly simulated.

Weaving Mitochondrial DNA and Y-Chromosome Variation in the Panamanian Genetic Canvas

The Isthmus of Panama was a crossroads between North and South America during the continent’s first peopling (and subsequent movements) also playing a pivotal role during European colonization and the African slave trade. Previous analyses of uniparental systems revealed significant sex biases in the genetic history of Panamanians, as testified by the high proportions of Indigenous and sub-Saharan mitochondrial DNAs (mtDNAs) and by the prevalence of Western European/northern African Y chromosomes. Those studies were conducted on the general population without considering any self-reported ethnic affiliations. Here, we compared the mtDNA and Y-chromosome lineages of a new sample collection from 431 individuals (301 males and 130 females) belonging to either the general population, mixed groups, or one of five Indigenous groups currently living in Panama. We found different proportions of paternal and maternal lineages in the Indigenous groups testifying to pre-contact demographic events and genetic inputs (some dated to Pleistocene times) that created genetic structure. Then, while the local mitochondrial gene pool was marginally involved in post-contact admixtures, the Indigenous Y chromosomes were differentially replaced, mostly by lineages of western Eurasian origin. Finally, our new estimates of the sub-Saharan contribution, on a more accurately defined general population, reduce an apparent divergence between genetic and historical data.

Assessing the Forensic Value of DNA Evidence from Y Chromosomes and Mitogenomes

Y chromosome and mitochondrial DNA profiles have been used as evidence in courts for decades, yet the problem of evaluating the weight of evidence has not been adequately resolved. Both are lineage markers (inherited from just one parent), which presents different interpretation challenges compared with standard autosomal DNA profiles (inherited from both parents). We review approaches to the evaluation of lineage marker profiles for forensic identification, focussing on the key roles of profile mutation rate and relatedness (extending beyond known relatives). Higher mutation rates imply fewer individuals matching the profile of an alleged contributor, but they will be more closely related. This makes it challenging to evaluate the possibility that one of these matching individuals could be the true source, because relatives may be plausible alternative contributors, and may not be well mixed in the population. These issues reduce the usefulness of profile databases drawn from a broad population: larger populations can have a lower profile relative frequency because of lower relatedness with the alleged contributor. Many evaluation methods do not adequately take account of distant relatedness, but its effects have become more pronounced with the latest generation of high-mutation-rate Y profiles.

Genome-scale sequencing and analysis of human, wolf, and bison DNA from 25,000-year-old sediment

Cave sediments have been shown to preserve ancient DNA but so far have not yielded the genome-scale information of skeletal remains. We retrieved and analyzed human and mammalian nuclear and mitochondrial environmental "shotgun" genomes from a single 25,000-year-old Upper Paleolithic sediment sample from Satsurblia cave, western Georgia:first, a human environmental genome with substantial basal Eurasian ancestry, which was an ancestral component of the majority of post-Ice Age people in the Near East, North Africa, and parts of Europe; second, a wolf environmental genome that is basal to extant Eurasian wolves and dogs and represents a previously unknown, likely extinct, Caucasian lineage; and third, a European bison environmental genome that is basal to present-day populations, suggesting that population structure has been substantially reshaped since the Last Glacial Maximum. Our results provide new insights into the Late Pleistocene genetic histories of these three species and demonstrate that direct shotgun sequencing of sediment DNA, without target enrichment methods, can yield genome-wide data informative of ancestry and phylogenetic relationships.

Ancient mitochondrial genomes from the Argentinian Pampas inform the early peopling of the Southern Cone of South America

The Southern Cone of South America (SCSA) is a key region for investigations about the peopling of the Americas. However, little is known about the eastern sector, the Argentinian Pampas. We analyzed 18 mitochondrial genomes—7 of which are novel—from human skeletal remains from 3 Early to Late Holocene archaeological sites. The Pampas present a distinctive genetic makeup compared to other Middle to Late Holocene pre-Columbian SCSA populations. We also report the earliest individuals carrying SCSA-specific mitochondrial haplogroups D1j and D1g from Early and Middle Holocene, respectively. Using these deep calibration time points in Bayesian phylogenetic reconstructions, we suggest that the first settlers of the Pampas were part of a single and rapid dispersal ∼15,600 years ago. Finally, we propose that present-day genetic differences between the Pampas and the rest of the SCSA are due to founder effects, genetic drift, and a partial population replacement ∼9,000 years ago.

•

Analysis of 18 ancient human mitochondrial genomes from the Argentinian Pampas.

•

Genetic makeup of Early-Mid Holocene Pampas distinct from later neighboring peoples.

•

Earliest individuals carrying region-specific mitochondrial haplogroups D1j and D1g.

•

First Pampean settlers were part of a single and rapid dispersal ∼15,600 years ago.

Exploring the phylogeography and population dynamics of the giant deer (Megaloceros giganteus) using Late Quaternary mitogenomes

Late Quaternary climatic fluctuations in the Northern Hemisphere had drastic effects on large mammal species, leading to the extinction of a substantial number of them. The giant deer (Megaloceros giganteus) was one of the species that became extinct in the Holocene, around 7660 calendar years before present. In the Late Pleistocene, the species ranged from western Europe to central Asia. However, during the Holocene, its range contracted to eastern Europe and western Siberia, where the last populations of the species occurred. Here, we generated 35 Late Pleistocene and Holocene giant deer mitogenomes to explore the genetics of the demise of this iconic species. Bayesian phylogenetic analyses of the mitogenomes suggested five main clades for the species: three pre-Last Glacial Maximum clades that did not appear in the post-Last Glacial Maximum genetic pool, and two clades that showed continuity into the Holocene. Our study also identified a decrease in genetic diversity starting in Marine Isotope Stage 3 and accelerating during the Last Glacial Maximum. This reduction in genetic diversity during the Last Glacial Maximum, coupled with a major contraction of fossil occurrences, suggests that climate was a major driver in the dynamics of the giant deer.

Middle Pleistocene genome calibrates a revised evolutionary history of extinct cave bears

Palaeogenomes provide the potential to study evolutionary processes in real time, but this potential is limited by our ability to recover genetic data over extended timescales.¹ As a consequence, most studies so far have focused on samples of Late Pleistocene or Holocene age, which covers only a small part of the history of many clades and species. Here, we report the recovery of a low coverage palaeogenome from the petrous bone of a ∼360,000 year old cave bear from Kudaro 1 cave in the Caucasus Mountains. Analysis of this genome alongside those of several Late Pleistocene cave bears reveals widespread mito-nuclear discordance in this group. Using the time interval between Middle and Late Pleistocene cave bear genomes, we directly estimate ursid nuclear and mitochondrial substitution rates to calibrate their respective phylogenies. This reveals post-divergence mitochondrial transfer as the dominant factor explaining their mito-nuclear discordance. Interestingly, these transfer events were not accompanied by large-scale nuclear introgression. However, we do detect additional instances of nuclear admixture among other cave bear lineages, and between cave bears and brown bears, which are not associated with mitochondrial exchange. Genomic data obtained from the Middle Pleistocene cave bear petrous bone has thus facilitated a revised evolutionary history of this extinct megafaunal group. Moreover, it suggests that petrous bones may provide a means of extending both the magnitude and time depth of palaeogenome retrieval over substantial portions of the evolutionary histories of many mammalian clades.

Human molecular evolutionary rate, time dependency and transient polymorphism effects viewed through ancient and modern mitochondrial DNA genomes

Human evolutionary genetics gives a chronological framework to interpret the human history. It is based on the molecular clock hypothesis that suppose a straightforward relationship between the mutation rate and the substitution rate with independence of other factors as demography dynamics. Analyzing ancient and modern human complete mitochondrial genomes we show here that, along the time, the substitution rate can be significantly slower or faster than the average germline mutation rate confirming a time dependence effect mainly attributable to changes in the effective population size of the human populations, with an exponential growth in recent times. We also detect that transient polymorphisms play a slowdown role in the evolutionary rate deduced from haplogroup intraspecific trees. Finally, we propose the use of the most divergent lineages within haplogroups as a practical approach to correct these molecular clock mismatches.

Sometimes hidden but always there: the assumptions underlying genetic inference of demographic histories

Demographic processes directly affect patterns of genetic variation within contemporary populations as well as future generations, allowing for demographic inference from patterns of both present-day and past genetic variation. Advances in laboratory procedures, sequencing and genotyping technologies in the past decades have resulted in massive increases in high-quality genome-wide genetic data from present-day populations and allowed retrieval of genetic data from archaeological material, also known as ancient DNA. This has resulted in an explosion of work exploring past changes in population size, structure, continuity and movement. However, as genetic processes are highly stochastic, patterns of genetic variation only indirectly reflect demographic histories. As a result, past demographic processes need to be reconstructed using an inferential approach. This usually involves comparing observed patterns of variation with model expectations from theoretical population genetics. A large number of approaches have been developed based on different population genetic models that each come with assumptions about the data and underlying demography. In this article I review some of the key models and assumptions underlying the most commonly used approaches for past demographic inference and their consequences for our ability to link the inferred demographic processes to the archaeological and climate records. This article is part of the theme issue 'Cross-disciplinary approaches to prehistoric demography'.

Utilizing auxology to understand ontogeny of extinct hominins: A case study on Homo naledi

The methods used to study human growth and development (auxology) have not previously been applied within the setting of hominin maturation (ontogeny). Ontogeny is defined here as the pattern of biological change into an adult form, both at the individual and species level. The hominin fossil record has a lack of recovered immature materials, due to such factors as taphonomic processes that destroy pre-adults; the fragility of immature compared to adult bone; and the lower mortality rates of juveniles compared to adults. The recent discovery of pre-adult hominin skeletal material from a single, homogeneous Homo naledi species from the Rising Star cave system in South Africa provides the opportunity for a broader application of auxology methods and thus the need to understand their use in a modern context. Human auxology studies benefit from a robust database, across multiple populations, and with longitudinal studies in order to assess the patterns and variations in typical growth, development and life history stages. Here, we review the approach, vocabulary, and methods of these human studies, investigate commonalities in data with the fossil record, and then advance the reconstruction of ontogeny for the extinct hominin species H. naledi. To this end, we apply an auxology model into the paleontological context to broadly predict H. naledi birthweight of the offspring at 2.06 kg with a range (±1 SD) of 1.89 to 2.24 kg, with a length at birth 45.5 cm. We estimate a H. naledi juvenile partial skeleton DH7 to be a height of 111-125 cm at death.

Complete chemical structures of human mitochondrial tRNAs

Mitochondria generate most cellular energy via oxidative phosphorylation. Twenty-two species of mitochondrial (mt-)tRNAs encoded in mtDNA translate essential subunits of the respiratory chain complexes. mt-tRNAs contain post-transcriptional modifications introduced by nuclear-encoded tRNA-modifying enzymes. They are required for deciphering genetic code accurately, as well as stabilizing tRNA. Loss of tRNA modifications frequently results in severe pathological consequences. Here, we perform a comprehensive analysis of post-transcriptional modifications of all human mt-tRNAs, including 14 previously-uncharacterized species. In total, we find 18 kinds of RNA modifications at 137 positions (8.7% in 1575 nucleobases) in 22 species of human mt-tRNAs. An up-to-date list of 34 genes responsible for mt-tRNA modifications are provided. We identify two genes required for queuosine (Q) formation in mt-tRNAs. Our results provide insight into the molecular mechanisms underlying the decoding system and could help to elucidate the molecular pathogenesis of human mitochondrial diseases caused by aberrant tRNA modifications.

Counterbalancing the time-dependent effect on the human mitochondrial DNA molecular clock

BACKGROUND

The molecular clock is an important genetic tool for estimating evolutionary timescales. However, the detection of a time-dependent effect on substitution rate estimates complicates its application. It has been suggested that demographic processes could be the main cause of this confounding effect. In the present study, I propose a new algorithm for estimating the coalescent age of phylogenetically related sequences, taking into account the observed time-dependent effect on the molecular rate detected by others.

RESULTS

By applying this method to real human mitochondrial DNA trees with shallow and deep topologies, I obtained significantly older molecular ages for the main events of human evolution than were previously estimated. These ages are in close agreement with the most recent archaeological and paleontological records favoring the emergence of early anatomically modern humans in Africa 315 ± 34 thousand years ago (kya) and the presence of recent modern humans outside of Africa as early as 174 ± 48 thousand years ago. Furthermore, during the implementation process, I demonstrated that in a population with fluctuating sizes, the probability of fixation of a new neutral mutant depends on the effective population size, which is in better accordance with the fact that under the neutral theory of molecular evolution, the fate of a molecular mutation is mainly determined by random drift.

CONCLUSIONS

I suggest that the demographic history of populations has a more decisive effect than purifying selection and/or mutational saturation on the time-dependent effect observed for the substitution rate, and I propose a new method that corrects for this effect.

Ancient DNA suggests modern wolves trace their origin to a Late Pleistocene expansion from Beringia

Grey wolves (Canis lupus) are one of the few large terrestrial carnivores that have maintained a wide geographical distribution across the Northern Hemisphere throughout the Pleistocene and Holocene. Recent genetic studies have suggested that, despite this continuous presence, major demographic changes occurred in wolf populations between the Late Pleistocene and early Holocene, and that extant wolves trace their ancestry to a single Late Pleistocene population. Both the geographical origin of this ancestral population and how it became widespread remain unknown. Here, we used a spatially and temporally explicit modelling framework to analyse a data set of 90 modern and 45 ancient mitochondrial wolf genomes from across the Northern Hemisphere, spanning the last 50,000 years. Our results suggest that contemporary wolf populations trace their ancestry to an expansion from Beringia at the end of the Last Glacial Maximum, and that this process was most likely driven by Late Pleistocene ecological fluctuations that occurred across the Northern Hemisphere. This study provides direct ancient genetic evidence that long‐range migration has played an important role in the population history of a large carnivore, and provides insight into how wolves survived the wave of megafaunal extinctions at the end of the last glaciation. Moreover, because Late Pleistocene grey wolves were the likely source from which all modern dogs trace their origins, the demographic history described in this study has fundamental implications for understanding the geographical origin of the dog.

see also the Perspective by Rena M. Schweizer and Robert K. Wayne.

The molecular clock of Mycobacterium tuberculosis

The molecular clock and its phylogenetic applications to genomic data have changed how we study and understand one of the major human pathogens, Mycobacterium tuberculosis (MTB), the etiologic agent of tuberculosis. Genome sequences of MTB strains sampled at different times are increasingly used to infer when a particular outbreak begun, when a drug-resistant clone appeared and expanded, or when a strain was introduced into a specific region. Despite the growing importance of the molecular clock in tuberculosis research, there is a lack of consensus as to whether MTB displays a clocklike behavior and about its rate of evolution. Here we performed a systematic study of the molecular clock of MTB on a large genomic data set (6,285 strains), covering different epidemiological settings and most of the known global diversity. We found that sampling times below 15-20 years were often insufficient to calibrate the clock of MTB. For data sets where such calibration was possible, we obtained a clock rate between 1x10-8 and 5x10-7 nucleotide changes per-site-per-year (0.04-2.2 SNPs per-genome-per-year), with substantial differences between clades. These estimates were not strongly dependent on the time of the calibration points as they changed only marginally when we used epidemiological isolates (sampled in the last 40 years) or three ancient DNA samples (about 1,000 years old) to calibrate the tree. Additionally, the uncertainty and the discrepancies in the results of different methods were sometimes large, highlighting the importance of using different methods, and of considering carefully their assumptions and limitations.

Evolutionary history and palaeoecology of brown bear in North-East Siberia re-examined using ancient DNA and stable isotopes from skeletal remains

Over 60% of the modern distribution range of brown bears falls within Russia, yet palaeoecological data from the region remain scarce. Complete modern Russian brown bear mitogenomes are abundant in the published literature, yet examples of their ancient counterparts are absent. Similarly, there is only limited stable isotopic data of prehistoric brown bears from the region. We used ancient DNA and stable carbon (δ¹³C) and nitrogen (δ¹⁵N) isotopes retrieved from five Pleistocene Yakutian brown bears (one Middle Pleistocene and four Late Pleistocene), to elucidate the evolutionary history and palaeoecology of the species in the region. We were able to reconstruct the complete mitogenome of one of the Late Pleistocene specimens, but we were unable to assign it to any of the previously published brown bear mitogenome clades. A subsequent analysis of published mtDNA control region sequences, which included sequences of extinct clades from other geographic regions, assigned the ancient Yakutian bear to the extinct clade 3c; a clade previously identified from Late Quaternary specimens from Eastern Beringia and Northern Spain. Our analyses of stable isotopes showed relatively high δ¹⁵N values in the Pleistocene Yakutian brown bears, suggesting a more carnivorous diet than contemporary brown bears from Eastern Beringia.

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.

How many individuals share a mitochondrial genome?

Mitochondrial DNA (mtDNA) is useful to assist with identification of the source of a biological sample, or to confirm matrilineal relatedness. Although the autosomal genome is much larger, mtDNA has an advantage for forensic applications of multiple copy number per cell, allowing better recovery of sequence information from degraded samples. In addition, biological samples such as fingernails, old bones, teeth and hair have mtDNA but little or no autosomal DNA. The relatively low mutation rate of the mitochondrial genome (mitogenome) means that there can be large sets of matrilineal-related individuals sharing a common mitogenome. Here we present the mitolina simulation software that we use to describe the distribution of the number of mitogenomes in a population that match a given mitogenome, and investigate its dependence on population size and growth rate, and on a database count of the mitogenome. Further, we report on the distribution of the number of meioses separating pairs of individuals with matching mitogenome. Our results have important implications for assessing the weight of mtDNA profile evidence in forensic science, but mtDNA analysis has many non-human applications, for example in tracking the source of ivory. Our methods and software can also be used for simulations to help validate models of population history in human or non-human populations.

Global demographic history of human populations inferred from whole mitochondrial genomes

The Neolithic transition has led to marked increases in census population sizes across the world, as recorded by a rich archaeological record. However, previous attempts to detect such changes using genetic markers, especially mitochondrial DNA (mtDNA), have mostly been unsuccessful. We use complete mtDNA genomes from over 1700 individuals, from the 1000 Genomes Project Phase 3, to explore changes in populations sizes in five populations for each of four major geographical regions, using a sophisticated coalescent-based Bayesian method (extended Bayesian skyline plots) and mutation rates calibrated with ancient DNA. Despite the power and sophistication of our analysis, we fail to find size changes that correspond to the Neolithic transitions of the study populations. However, we do detect a number of size changes, which tend to be replicated in most populations within each region. These changes are mostly much older than the Neolithic transition and could reflect either population expansion or changes in population structure. Given the amount of migration and population mixing that occurred after these ancient signals were generated, we caution that modern populations will often carry ghost signals of demographic events that occurred far away from their current location.

Using Phylogenomic Data to Explore the Effects of Relaxed Clocks and Calibration Strategies on Divergence Time Estimation: Primates as a Test Case

Primates have long been a test case for the development of phylogenetic methods for divergence time estimation. Despite a large number of studies, however, the timing of origination of crown Primates relative to the Cretaceous-Paleogene (K-Pg) boundary and the timing of diversification of the main crown groups remain controversial. Here, we analysed a data set of 372 taxa (367 Primates and 5 outgroups, 3.4 million aligned base pairs) that includes nine primate genomes. We systematically explore the effect of different interpretations of fossil calibrations and molecular clock models on primate divergence time estimates. We find that even small differences in the construction of fossil calibrations can have a noticeable impact on estimated divergence times, especially for the oldest nodes in the tree. Notably, choice of molecular rate model (autocorrelated or independently distributed rates) has an especially strong effect on estimated times, with the independent rates model producing considerably more ancient age estimates for the deeper nodes in the phylogeny. We implement thermodynamic integration, combined with Gaussian quadrature, in the program MCMCTree, and use it to calculate Bayes factors for clock models. Bayesian model selection indicates that the autocorrelated rates model fits the primate data substantially better, and we conclude that time estimates under this model should be preferred. We show that for eight core nodes in the phylogeny, uncertainty in time estimates is close to the theoretical limit imposed by fossil uncertainties. Thus, these estimates are unlikely to be improved by collecting additional molecular sequence data. All analyses place the origin of Primates close to the K-Pg boundary, either in the Cretaceous or straddling the boundary into the Palaeogene.

A comparison of methods for estimating substitution rates from ancient DNA sequence data

Background

Phylogenetic analysis of DNA from modern and ancient samples allows the reconstruction of important demographic and evolutionary processes. A critical component of these analyses is the estimation of evolutionary rates, which can be calibrated using information about the ages of the samples. However, the reliability of these rate estimates can be negatively affected by among-lineage rate variation and non-random sampling. Using a simulation study, we compared the performance of three phylogenetic methods for inferring evolutionary rates from time-structured data sets: regression of root-to-tip distances, least-squares dating, and Bayesian inference. We also applied these three methods to time-structured mitogenomic data sets from six vertebrate species.

Results

Our results from 12 simulation scenarios show that the three methods produce reliable estimates when the substitution rate is high, rate variation is low, and samples of similar ages are not all grouped together in the tree (i.e., low phylo-temporal clustering). The interaction of these factors is particularly important for least-squares dating and Bayesian estimation of evolutionary rates. The three estimation methods produced consistent estimates of rates across most of the six mitogenomic data sets, with sequence data from horses being an exception.

Conclusions

We recommend that phylogenetic studies of ancient DNA sequences should use multiple methods of inference and test for the presence of temporal signal, among-lineage rate variation, and phylo-temporal clustering in the data.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1192-3) contains supplementary material, which is available to authorized users.

Significant loss of mitochondrial diversity within the last century due to extinction of peripheral populations in eastern gorillas

Species and populations are disappearing at an alarming rate as a direct result of human activities. Loss of genetic diversity associated with population decline directly impacts species' long-term survival. Therefore, preserving genetic diversity is of considerable conservation importance. However, to assist in conservation efforts, it is important to understand how genetic diversity is spatially distributed and how it changes due to anthropogenic pressures. In this study, we use historical museum and modern faecal samples of two critically endangered eastern gorilla taxa, Grauer's (Gorilla beringei graueri) and mountain gorillas (Gorilla beringei beringei), to directly infer temporal changes in genetic diversity within the last century. Using over 100 complete mitochondrial genomes, we observe a significant decline in haplotype and nucleotide diversity in Grauer's gorillas. By including historical samples from now extinct populations we show that this decline can be attributed to the loss of peripheral populations rather than a decrease in genetic diversity within the core range of the species. By directly quantifying genetic changes in the recent past, our study shows that human activities have severely impacted eastern gorilla genetic diversity within only four to five generations. This rapid loss calls for dedicated conservation actions, which should include preservation of the remaining peripheral populations.

Genetic Ancestry of Hadza and Sandawe Peoples Reveals Ancient Population Structure in Africa

The Hadza and Sandawe populations in present-day Tanzania speak languages containing click sounds and therefore thought to be distantly related to southern African Khoisan languages. We analyzed genome-wide genotype data for individuals sampled from the Hadza and Sandawe populations in the context of a global data set of 3,528 individuals from 163 ethno-linguistic groups. We found that Hadza and Sandawe individuals share ancestry distinct from and most closely related to Omotic ancestry; share Khoisan ancestry with populations such as ≠Khomani, Karretjie, and Ju/'hoansi in southern Africa; share Niger-Congo ancestry with populations such as Yoruba from Nigeria and Luhya from Kenya, consistent with migration associated with the Bantu Expansion; and share Cushitic ancestry with Somali, multiple Ethiopian populations, the Maasai population in Kenya, and the Nama population in Namibia. We detected evidence for low levels of Arabian, Nilo-Saharan, and Pygmy ancestries in a minority of individuals. Our results indicate that west Eurasian ancestry in eastern Africa is more precisely the Arabian parent of Cushitic ancestry. Relative to the Out-of-Africa migrations, Hadza ancestry emerged early whereas Sandawe ancestry emerged late.

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.

Investigating the origins of eastern Polynesians using genome-wide data from the Leeward Society Isles

The debate concerning the origin of the Polynesian speaking peoples has been recently reinvigorated by genetic evidence for secondary migrations to western Polynesia from the New Guinea region during the 2nd millennium BP. Using genome-wide autosomal data from the Leeward Society Islands, the ancient cultural hub of eastern Polynesia, we find that the inhabitants' genomes also demonstrate evidence of this episode of admixture, dating to 1,700-1,200 BP. This supports a late settlement chronology for eastern Polynesia, commencing ~1,000 BP, after the internal differentiation of Polynesian society. More than 70% of the autosomal ancestry of Leeward Society Islanders derives from Island Southeast Asia with the lowland populations of the Philippines as the single largest potential source. These long-distance migrants into Polynesia experienced additional admixture with northern Melanesians prior to the secondary migrations of the 2nd millennium BP. Moreover, the genetic diversity of mtDNA and Y chromosome lineages in the Leeward Society Islands is consistent with linguistic evidence for settlement of eastern Polynesia proceeding from the central northern Polynesian outliers in the Solomon Islands. These results stress the complex demographic history of the Leeward Society Islands and challenge phylogenetic models of cultural evolution predicated on eastern Polynesia being settled from Samoa.

Matriclans shape populations: Insights from the Angolan Namib Desert into the maternal genetic history of southern Africa

OBJECTIVES

Southern Angola is a poorly studied region, inhabited by populations that have been associated with different migratory movements into southern Africa. Apart from Kx'a-speaking San foragers and Bantu-speaking pastoralists, ethnographic and linguistic studies have suggested the existence of an enigmatic array of pre-Bantu communities, like the Kwepe (formerly Khoe-Kwadi speakers), Twa and Kwisi. Here, we evaluate previous peopling hypotheses by assessing the relationships between different southern Angolan populations, based on newly collected linguistic data and complete mtDNA genomes.

MATERIALS AND METHODS

We analyzed 295 complete mtDNA genomes and linguistic data from seven groups from the Namib Desert (Himba, Kuvale, Tjimba, Twa, Kwisi, Kwepe) and Kunene Province (!Xun), placing special emphasis on the evaluation of the genealogical consistency of the matriclanic system that characterizes most of these groups.

RESULTS

We found that the maternal genetic structure of all groups from the Namib Desert was strongly shaped by the consistency of their matriclanic system. The tracking of the maternal heritage enhanced population differentiation by genetic drift and is likely to have caused the divergent mtDNA profiles of the Kwepe, Twa, and Kwisi, who probably formed a single population within the spectrum of Bantu genetic variation. Model-based analyses further suggest that the dominant pastoral groups Kuvale and Himba may be grouped into a Bantu proto-population which also included the ancestors of present-day Tjimba and Herero, as well as the Khoe-Kwadi speaking Damara foragers from Namibia.

DISCUSSION

The view from southwestern Angola offers a new perspective on the populating history of southern Africa and the Bantu expansions by showing that social stratification and different subsistence patterns are not always indicative of remnant groups, but may reflect Bantu-internal variation and ethnogenesis.

Investigating population continuity with ancient DNA under a spatially explicit simulation framework

BACKGROUND

Recent advances in sequencing technologies have allowed for the retrieval of ancient DNA data (aDNA) from skeletal remains, providing direct genetic snapshots from diverse periods of human prehistory. Comparing samples taken in the same region but at different times, hereafter called "serial samples", may indicate whether there is continuity in the peopling history of that area or whether an immigration of a genetically different population has occurred between the two sampling times. However, the exploration of genetic relationships between serial samples generally ignores their geographical locations and the spatiotemporal dynamics of populations. Here, we present a new coalescent-based, spatially explicit modelling approach to investigate population continuity using aDNA, which includes two fundamental elements neglected in previous methods: population structure and migration. The approach also considers the extensive temporal and geographical variance that is commonly found in aDNA population samples.

RESULTS

We first showed that our spatially explicit approach is more conservative than the previous (panmictic) approach and should be preferred to test for population continuity, especially when small and isolated populations are considered. We then applied our method to two mitochondrial datasets from Germany and France, both including modern and ancient lineages dating from the early Neolithic. The results clearly reject population continuity for the maternal line over the last 7500 years for the German dataset but not for the French dataset, suggesting regional heterogeneity in post-Neolithic migratory processes.

CONCLUSIONS

Here, we demonstrate the benefits of using a spatially explicit method when investigating population continuity with aDNA. It constitutes an improvement over panmictic methods by considering the spatiotemporal dynamics of genetic lineages and the precise location of ancient samples. The method can be used to investigate population continuity between any pair of serial samples (ancient-ancient or ancient-modern) and to investigate more complex evolutionary scenarios. Although we based our study on mitochondrial DNA sequences, diploid molecular markers of different types (DNA, SNP, STR) can also be simulated with our approach. It thus constitutes a promising tool for the analysis of the numerous aDNA datasets being produced, including genome wide data, in humans but also in many other species.

Identification and analysis of mtDNA genomes attributed to Finns reveal long-stagnant demographic trends obscured in the total diversity

In Europe, modern mitochondrial diversity is relatively homogeneous and suggests an ubiquitous rapid population growth since the Neolithic revolution. Similar patterns also have been observed in mitochondrial control region data in Finland, which contrasts with the distinctive autosomal and Y-chromosomal diversity among Finns. A different picture emerges from the 843 whole mitochondrial genomes from modern Finns analyzed here. Up to one third of the subhaplogroups can be considered as Finn-characteristic, i.e. rather common in Finland but virtually absent or rare elsewhere in Europe. Bayesian phylogenetic analyses suggest that most of these attributed Finnish lineages date back to around 3,000–5,000 years, coinciding with the arrival of Corded Ware culture and agriculture into Finland. Bayesian estimation of past effective population sizes reveals two differing demographic histories: 1) the ‘local’ Finnish mtDNA haplotypes yielding small and dwindling size estimates for most of the past; and 2) the ‘immigrant’ haplotypes showing growth typical of most European populations. The results based on the local diversity are more in line with that known about Finns from other studies, e.g., Y-chromosome analyses and archaeology findings. The mitochondrial gene pool thus may contain signals of local population history that cannot be readily deduced from the total diversity.

Deeply divergent archaic mitochondrial genome provides lower time boundary for African gene flow into Neanderthals

Ancient DNA is revealing new insights into the genetic relationship between Pleistocene hominins and modern humans. Nuclear DNA indicated Neanderthals as a sister group of Denisovans after diverging from modern humans. However, the closer affinity of the Neanderthal mitochondrial DNA (mtDNA) to modern humans than Denisovans has recently been suggested as the result of gene flow from an African source into Neanderthals before 100,000 years ago. Here we report the complete mtDNA of an archaic femur from the Hohlenstein–Stadel (HST) cave in southwestern Germany. HST carries the deepest divergent mtDNA lineage that splits from other Neanderthals ∼270,000 years ago, providing a lower boundary for the time of the putative mtDNA introgression event. We demonstrate that a complete Neanderthal mtDNA replacement is feasible over this time interval even with minimal hominin introgression. The highly divergent HST branch is indicative of greater mtDNA diversity during the Middle Pleistocene than in later periods.

Ancient DNA keeps expanding our understanding of complex genetic relationships between Pleistocene hominins. Here, Posth and colleagues analyse the mitochondrial genome of an archaic human that diverged from other Neanderthals ∼270,000 years ago, providing the minimum age for an African introgression into Neanderthals.

tipdatingbeast: an r package to assist the implementation of phylogenetic tip-dating tests using beast

Molecular tip dating of phylogenetic trees is a growing discipline that uses DNA sequences sampled at different points in time to coestimate the timing of evolutionary events with rates of molecular evolution. In this context, beast, a program for Bayesian analysis of molecular sequences, is the most widely used phylogenetic tool. Here, we introduce tipdatingbeast, an r package built to assist the implementation of various phylogenetic tip-dating tests using beast. tipdatingbeast currently contains two main functions. The first one allows preparing date-randomization analyses, which assess the temporal signal of a data set. The second function allows performing leave-one-out analyses, which test for the consistency between independent calibration sequences and allow pinpointing those leading to potential bias. We apply those functions to an empirical data set and supply practical guidance for results interpretation.