A "Copernican" reassessment of the human mitochondrial DNA tree from its root

The discovery of a ten-generation m.C1494T pedigree in the east of England with probable links to King Richard III

This paper reports the discovery of a m.C1494T pedigree in the east of England made during a search for matrilineal relations of King Richard III. The mitochondrial DNA variant m.C1494T has been associated with aminoglycoside-induced deafness. This variant is very uncommon. although pedigrees with this variant have previously been found in China and Spain. The members of the newly identified pedigree all belong to the mitochondrial haplogroup J1c2c3, which is also the haplogroup of King Richard III. The presence of a few people in the USA from the same haplogroup has previously been noted, and it is now known that one of the people can show his descent from a couple who lived in Nottinghamshire, England, in the late 1700's. The mitochondrial DNA sequence of this man, at present living in the USA, and of his 4th cousin, twice removed, living in Lincoln, England, has shown they belong to haplogroup J1c2c3 and both have the variant m.C1494T; thereby, allowing the production of a multi-generational pedigree originating in the east of England. Fortunately, deafness has not been found in any living member of this large pedigree. It was also noted that the link to the family of King Richard III has not been firmly defined; however the circumstantial evidence is strong as many of his family members lived in this part of England.

CD59 gene: 143 haplotypes of 22,718 nucleotides length by computational phasing in 113 individuals from different ethnicities

BACKGROUND

CD59 deficiency due to rare germline variants in the CD59 gene causes disabilities, ischemic strokes, neuropathy, and hemolysis. CD59 deficiency due to common somatic variants in the PIG-A gene in hematopoietic stem cells causes paroxysmal nocturnal hemoglobinuria. The ISBT database lists one nonsense and three missense germline variants that are associated with the CD59-null phenotype. To analyze the genetic diversity of the CD59 gene, we determined long-range CD59 haplotypes among individuals from different ethnicities.

METHODS

We determined a 22.7 kb genomic fragment of the CD59 gene in 113 individuals using next-generation sequencing (NGS), which covered the whole NM_203330.2 mRNA transcript of 7796 base pairs. Samples came from an FDA reference repository and our Ethiopia study cohorts. The raw genotype data were computationally phased into individual haplotype sequences.

RESULTS

Nucleotide sequencing of the CD59 gene of 226 chromosomes identified 216 positions with single nucleotide variants. Only three haplotypes were observed in homozygous form, which allowed us to assign them unambiguously as experimentally verified CD59 haplotypes. They were also the most frequent haplotypes among both cohorts. An additional 140 haplotypes were imputed computationally.

DISCUSSION

We provided a large set of haplotypes and proposed three verified long-range CD59 reference sequences, based on a population approach, using a generalizable rationale for our choice. Correct long-range haplotypes are useful as template sequences for allele calling in high-throughput NGS and precision medicine approaches, thus enhancing the reliability of clinical diagnostics. Long-range haplotypes can also be used to evaluate the influence of genetic variation on the risk of transfusion reactions or diseases.

Investigating demic versus cultural diffusion and sex bias in the spread of Austronesian languages in Vietnam

Austronesian (AN) is the second-largest language family in the world, particularly widespread in Island Southeast Asia (ISEA) and Oceania. In Mainland Southeast Asia (MSEA), groups speaking these languages are concentrated in the highlands of Vietnam. However, our knowledge of the spread of AN-speaking populations in MSEA remains limited; in particular, it is not clear if AN languages were spread by demic or cultural diffusion. In this study, we present and analyze new data consisting of complete mitogenomes from 369 individuals and 847 Y-chromosomal single nucleotide polymorphisms (SNPs) from 170 individuals from all five Vietnamese Austronesian groups (VN-AN) and five neighboring Vietnamese Austroasiatic groups (VN-AA). We found genetic signals consistent with matrilocality in some, but not all, of the VN-AN groups. Population affinity analyses indicated connections between the AN-speaking Giarai and certain Taiwanese AN groups (Rukai, Paiwan, and Bunun). However, overall, there were closer genetic affinities between VN-AN groups and neighboring VN-AA groups, suggesting language shifts. Our study provides insights into the genetic structure of AN-speaking communities in MSEA, characterized by some contact with Taiwan and language shift in neighboring groups, indicating that the expansion of AN speakers in MSEA was a combination of cultural and demic diffusion.

Natural selection shaped the protective effect of the mtDNA lineage against obesity in Han Chinese populations

Mitochondria play a key role in lipid metabolism, and mitochondrial DNA (mtDNA) mutations are thus considered to affect obesity susceptibility by altering oxidative phosphorylation and mitochondrial function. In this study, we investigated mtDNA variants that may affect obesity risk in 2,877 Han Chinese individuals from three independent populations. The association analysis of 16 basal mtDNA haplogroups with body mass index (BMI), waist circumference (WC) and waist-to-hip ratio (WHR) revealed that only haplogroup M7 was significantly negatively correlated with all three adiposity-related anthropometric traits in the overall cohort (P=0.003 for BMI, P=1×10-5 for WC, P=0.005 for WHR), which was verified by the analysis of a single population, i.e., the Zhengzhou population. Furthermore, subhaplogroup analysis suggested that M7b1a1 was the most likely haplogroup associated with a decreased obesity risk, and the variation T12811C (causing Y159H in ND5) harbored in M7b1a1 may be the most likely candidate for altering mitochondrial function. Specifically, we found that proportionally more nonsynonymous mutations accumulated in M7b1a1 carriers, indicating that M7b1a1 was either under positive selection or subject to a relaxation of selective constraints. We also found that nuclear variants, especially in DACT2 and PIEZO1, may functionally interact with M7b1a1.

G-quadruplex propensity in H. neanderthalensis, H. sapiens and Denisovans mitochondrial genomes

Current methods of processing archaeological samples combined with advances in sequencing methods lead to disclosure of a large part of H. neanderthalensis and Denisovans genetic information. It is hardly surprising that the genome variability between modern humans, Denisovans and H. neanderthalensis is relatively limited. Genomic studies may provide insight on the metabolism of extinct human species or lineages. Detailed analysis of G-quadruplex sequences in H. neanderthalensis and Denisovans mitochondrial DNA showed us interesting features. Relatively similar patterns in mitochondrial DNA are found compared to modern humans, with one notable exception for H. neanderthalensis. An interesting difference between H. neanderthalensis and H. sapiens corresponds to a motif found in the D-loop region of mtDNA, which is responsible for mitochondrial DNA replication. This area is directly responsible for the number of mitochondria and consequently for the efficient energy metabolism of cell. H. neanderthalensis harbor a long uninterrupted run of guanines in this region, which may cause problems for replication, in contrast with H. sapiens, for which this run is generally shorter and interrupted. One may propose that the predominant H. sapiens motif provided a selective advantage for modern humans regarding mtDNA replication and function.

Mitochondrial DNA: Inherent Complexities Relevant to Genetic Analyses

Mitochondrial DNA (mtDNA) exhibits distinct characteristics distinguishing it from the nuclear genome, necessitating specific analytical methods in genetic studies. This comprehensive review explores the complex role of mtDNA in a variety of genetic studies, including genome-wide, epigenome-wide, and phenome-wide association studies, with a focus on its implications for human traits and diseases. Here, we discuss the structure and gene-encoding properties of mtDNA, along with the influence of environmental factors and epigenetic modifications on its function and variability. Particularly significant are the challenges posed by mtDNA's high mutation rate, heteroplasmy, and copy number variations, and their impact on disease susceptibility and population genetic analyses. The review also highlights recent advances in methodological approaches that enhance our understanding of mtDNA associations, advocating for refined genetic research techniques that accommodate its complexities. By providing a comprehensive overview of the intricacies of mtDNA, this paper underscores the need for an integrated approach to genetic studies that considers the unique properties of mitochondrial genetics. Our findings aim to inform future research and encourage the development of innovative methodologies to better interpret the broad implications of mtDNA in human health and disease.

The Genetic Origin of the Indo-Europeans

The Yamnaya archaeological complex appeared around 3300BCE across the steppes north of the Black and Caspian Seas, and by 3000BCE reached its maximal extent from Hungary in the west to Kazakhstan in the east. To localize the ancestral and geographical origins of the Yamnaya among the diverse Eneolithic people that preceded them, we studied ancient DNA data from 428 individuals of which 299 are reported for the first time, demonstrating three previously unknown Eneolithic genetic clines. First, a "Caucasus-Lower Volga" (CLV) Cline suffused with Caucasus hunter-gatherer (CHG) ancestry extended between a Caucasus Neolithic southern end in Neolithic Armenia, and a steppe northern end in Berezhnovka in the Lower Volga. Bidirectional gene flow across the CLV cline created admixed intermediate populations in both the north Caucasus, such as the Maikop people, and on the steppe, such as those at the site of Remontnoye north of the Manych depression. CLV people also helped form two major riverine clines by admixing with distinct groups of European hunter-gatherers. A "Volga Cline" was formed as Lower Volga people mixed with upriver populations that had more Eastern hunter-gatherer (EHG) ancestry, creating genetically hyper-variable populations as at Khvalynsk in the Middle Volga. A "Dnipro Cline" was formed as CLV people bearing both Caucasus Neolithic and Lower Volga ancestry moved west and acquired Ukraine Neolithic hunter-gatherer (UNHG) ancestry to establish the population of the Serednii Stih culture from which the direct ancestors of the Yamnaya themselves were formed around 4000BCE. This population grew rapidly after 3750-3350BCE, precipitating the expansion of people of the Yamnaya culture who totally displaced previous groups on the Volga and further east, while admixing with more sedentary groups in the west. CLV cline people with Lower Volga ancestry contributed four fifths of the ancestry of the Yamnaya, but also, entering Anatolia from the east, contributed at least a tenth of the ancestry of Bronze Age Central Anatolians, where the Hittite language, related to the Indo-European languages spread by the Yamnaya, was spoken. We thus propose that the final unity of the speakers of the "Proto-Indo-Anatolian" ancestral language of both Anatolian and Indo-European languages can be traced to CLV cline people sometime between 4400-4000 BCE.

Medieval DNA from Soqotra points to Eurasian origins of an isolated population at the crossroads of Africa and Arabia

Soqotra, an island situated at the mouth of the Gulf of Aden in the northwest Indian Ocean between Africa and Arabia, is home to ~60,000 people subsisting through fishing and semi-nomadic pastoralism who speak a Modern South Arabian language. Most of what is known about Soqotri history derives from writings of foreign travellers who provided little detail about local people, and the geographic origins and genetic affinities of early Soqotri people has not yet been investigated directly. Here we report genome-wide data from 39 individuals who lived between ~650 and 1750 CE at six locations across the island and document strong genetic connections between Soqotra and the similarly isolated Hadramawt region of coastal South Arabia that likely reflects a source for the peopling of Soqotra. Medieval Soqotri can be modelled as deriving ~86% of their ancestry from a population such as that found in the Hadramawt today, with the remaining ~14% best proxied by an Iranian-related source with up to 2% ancestry from the Indian sub-continent, possibly reflecting genetic exchanges that occurred along with archaeologically documented trade from these regions. In contrast to all other genotyped populations of the Arabian Peninsula, genome-level analysis of the medieval Soqotri is consistent with no sub-Saharan African admixture dating to the Holocene. The deep ancestry of people from medieval Soqotra and the Hadramawt is also unique in deriving less from early Holocene Levantine farmers and more from groups such as Late Pleistocene hunter-gatherers from the Levant (Natufians) than other mainland Arabians. This attests to migrations by early farmers having less impact in southernmost Arabia and Soqotra and provides compelling evidence that there has not been complete population replacement between the Pleistocene and Holocene throughout the Arabian Peninsula. Medieval Soqotra harboured a small population that showed qualitatively different marriage practices from modern Soqotri, with first-cousin unions occurring significantly less frequently than today.

Association analysis of mitochondrial DNA heteroplasmic variants: methods and application

We rigorously assessed a comprehensive association testing framework for heteroplasmy, employing both simulated and real-world data. This framework employed a variant allele fraction (VAF) threshold and harnessed multiple gene-based tests for robust identification and association testing of heteroplasmy. Our simulation studies demonstrated that gene-based tests maintained an appropriate type I error rate at α=0.001. Notably, when 5% or more heteroplasmic variants within a target region were linked to an outcome, burden-extension tests (including the adaptive burden test, variable threshold burden test, and z-score weighting burden test) outperformed the sequence kernel association test (SKAT) and the original burden test. Applying this framework, we conducted association analyses on whole-blood derived heteroplasmy in 17,507 individuals of African and European ancestries (31% of African Ancestry, mean age of 62, with 58% women) with whole genome sequencing data. We performed both cohort- and ancestry-specific association analyses, followed by meta-analysis on both pooled samples and within each ancestry group. Our results suggest that mtDNA-encoded genes/regions are likely to exhibit varying rates in somatic aging, with the notably strong associations observed between heteroplasmy in the RNR1 and RNR2 genes (p<0.001) and advance aging by the Original Burden test. In contrast, SKAT identified significant associations (p<0.001) between diabetes and the aggregated effects of heteroplasmy in several protein-coding genes. Further research is warranted to validate these findings. In summary, our proposed statistical framework represents a valuable tool for facilitating association testing of heteroplasmy with disease traits in large human populations.

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.

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.

Genetic Analysis of Mingrelians Reveals Long-Term Continuity of Populations in Western Georgia (Caucasus)

To elucidate the population history of the Caucasus, we conducted a survey of genetic diversity in Samegrelo (Mingrelia), western Georgia. We collected DNA samples and genealogical information from 485 individuals residing in 30 different locations, the vast majority of whom being Mingrelian speaking. From these DNA samples, we generated mitochondrial DNA (mtDNA) control region sequences for all 485 participants (female and male), Y-short tandem repeat haplotypes for the 372 male participants, and analyzed all samples at nearly 590,000 autosomal single nucleotide polymorphisms (SNPs) plus around 33,000 on the sex chromosomes, with 27,000 SNP removed for missingness, using the GenoChip 2.0+ microarray. The resulting data were compared with those from populations from Anatolia, the Caucasus, the Near East, and Europe. Overall, Mingrelians exhibited considerable mtDNA haplogroup diversity, having high frequencies of common West Eurasian haplogroups (H, HV, I, J, K, N1, R1, R2, T, U, and W. X2) and low frequencies of East Eurasian haplogroups (A, C, D, F, and G). From a Y-chromosome standpoint, Mingrelians possessed a variety of haplogroups, including E1b1b, G2a, I2, J1, J2, L, Q, R1a, and R1b. Analysis of autosomal SNP data further revealed that Mingrelians are genetically homogeneous and cluster with other modern-day South Caucasus populations. When compared with ancient DNA samples from Bronze Age archaeological contexts in the broader region, these data indicate that the Mingrelian gene pool began taking its current form at least by this period, probably in conjunction with the formation of a distinct linguistic community.

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.

South Asian maternal and paternal lineages in southern Thailand and the role of sex-biased admixture

Previous genome-wide studies have reported South Asian (SA) ancestry in several Mainland Southeast Asian (MSEA) populations; however, additional details concerning population history, in particular the role of sex-specific aspects of the SA admixture in MSEA populations can be addressed with uniparental markers. Here, we generated ∼2.3 mB sequences of the male-specific portions of the Y chromosome (MSY) of a Tai-Kadai (TK)-speaking Southern Thai group (SouthernThai_TK), and complete mitochondrial (mtDNA) genomes of the SouthernThai_TK and an Austronesian (AN)-speaking Southern Thai (SouthernThai_AN) group. We identified new mtDNA haplogroups, e.g. Q3, E1a1a1, B4a1a and M7c1c3 that have not previously reported in Thai populations, but are frequent in Island Southeast Asia and Oceania, suggesting interactions between MSEA and these regions. SA prevalent mtDNA haplogroups were observed at frequencies of ~35-45% in the Southern Thai groups; both of them showed more genetic relatedness to Austroasiatic (AA) speaking Mon than to any other group. For MSY, SouthernThai_TK had ~35% SA prevalent haplogroups and exhibited closer genetic affinity to Central Thais. We also analyzed published data from other MSEA populations and observed SA ancestry in some additional MSEA populations that also reflects sex-biased admixture; in general, most AA- and AN-speaking groups in MSEA were closer to SA than to TK groups based on mtDNA, but the opposite pattern was observed for the MSY. Overall, our results of new genetic lineages and sex-biased admixture from SA to MSEA groups attest to the additional value that uniparental markers can add to studies of genome-wide variation.

A systematic comparison of human mitochondrial genome assembly tools

BACKGROUND

Mitochondria are the cell organelles that produce most of the chemical energy required to power the cell's biochemical reactions. Despite being a part of a eukaryotic host cell, the mitochondria contain a separate genome whose origin is linked with the endosymbiosis of a prokaryotic cell by the host cell and encode independent genomic information throughout their genomes. Mitochondrial genomes accommodate essential genes and are regularly utilized in biotechnology and phylogenetics. Various assemblers capable of generating complete mitochondrial genomes are being continuously developed. These tools often use whole-genome sequencing data as an input containing reads from the mitochondrial genome. Till now, no published work has explored the systematic comparison of all the available tools for assembling human mitochondrial genomes using short-read sequencing data. This evaluation is required to identify the best tool that can be well-optimized for small-scale projects or even national-level research.

RESULTS

In this study, we have tested the mitochondrial genome assemblers for both simulated datasets and whole genome sequencing (WGS) datasets of humans. For the highest computational setting of 16 computational threads with the simulated dataset having 1000X read depth, MitoFlex took the least execution time of 69 s, and IOGA took the longest execution time of 1278 s. NOVOPlasty utilized the least computational memory of approximately 0.098 GB for the same setting, whereas IOGA utilized the highest computational memory of 11.858 GB. In the case of WGS datasets for humans, GetOrganelle and MitoFlex performed the best in capturing the SNPs information with a mean F1-score of 0.919 at the sequencing depth of 10X. MToolBox and NOVOPlasty performed consistently across all sequencing depths with a mean F1 score of 0.897 and 0.890, respectively.

CONCLUSIONS

Based on the overall performance metrics and consistency in assembly quality for all sequencing data, MToolBox performed the best. However, NOVOPlasty was the second fastest tool in execution time despite being single-threaded, and it utilized the least computational resources among all the assemblers when tested on simulated datasets. Therefore, NOVOPlasty may be more practical when there is a significant sample size and a lack of computational resources. Besides, as long-read sequencing gains popularity, mitochondrial genome assemblers must be developed to use long-read sequencing data.

The ancestry and geographical origins of St Helena's liberated Africans

The island of St Helena played a crucial role in the suppression of the transatlantic slave trade. Strategically located in the middle of the South Atlantic, it served as a staging post for the Royal Navy and reception point for enslaved Africans who had been "liberated" from slave ships intercepted by the British. In total, St Helena received approximately 27,000 liberated Africans between 1840 and 1867. Written sources suggest that the majority of these individuals came from West Central Africa, but their precise origins are unknown. Here, we report the results of ancient DNA analyses that we conducted as part of a wider effort to commemorate St Helena's liberated Africans and to restore knowledge of their lives and experiences. We generated partial genomes (0.1-0.5×) for 20 individuals whose remains had been recovered during archaeological excavations on the island. We compared their genomes with genotype data for over 3,000 present-day individuals from 90 populations across sub-Saharan Africa and conclude that the individuals most likely originated from different source populations within the general area between northern Angola and Gabon. We also find that the majority (17/20) of the individuals were male, supporting a well-documented sex bias in the latter phase of the transatlantic slave trade. The study expands our understanding of St Helena's liberated African community and illustrates how ancient DNA analyses can be used to investigate the origins and identities of individuals whose lives were bound up in the story of slavery and its abolition.

Descent, marriage, and residence practices of a 3,800-year-old pastoral community in Central Eurasia

Our understanding of prehistoric societal organization at the family level is still limited. Here, we generated genome data from 32 individuals from an approximately 3,800-y-old burial mound attributed to the Bronze Age Srubnaya-Alakul cultural tradition at the site of Nepluyevsky, located in the Southern Ural region of Central Eurasia. We found that life expectancy was generally very low, with adult males living on average 8 y longer than females. A total of 35 first-degree, 40 second-degree, and 48 third-degree biological relationships connected 23 of the studied individuals, allowing us to propose a family tree spanning three generations with six brothers at its center. The oldest of these brothers had eight children with two women and the most children overall, whereas the other relationships were monogamous. Notably, related female children above the age of five were completely absent from the site, and adult females were more genetically diverse than males. These results suggest that biological relationships between male siblings played a structural role in society and that descent group membership was based on patrilineality. Women originated from a larger mating network and moved to join the men, with whom they were buried. Finally, the oldest brother likely held a higher social position, which was expressed in terms of fertility.

A Customized Human Mitochondrial DNA Database (hMITO DB v1.0) for Rapid Sequence Analysis, Haplotyping and Geo-Mapping

The field of mitochondrial genomics has advanced rapidly and has revolutionized disciplines such as molecular anthropology, population genetics, and medical genetics/oncogenetics. However, mtDNA next-generation sequencing (NGS) analysis for matrilineal haplotyping and phylogeographic inference remains hindered by the lack of a consolidated mitogenome database and an efficient bioinformatics pipeline. To address this, we developed a customized human mitogenome database (hMITO DB) embedded in a CLC Genomics workflow for read mapping, variant analysis, haplotyping, and geo-mapping. The database was constructed from 4286 mitogenomes. The macro-haplogroup (A to Z) distribution and representative phylogenetic tree were found to be consistent with published literature. The hMITO DB automated workflow was tested using mtDNA-NGS sequences derived from Pap smears and cervical cancer cell lines. The auto-generated read mapping, variants track, and table of haplotypes and geo-origins were completed in 15 min for 47 samples. The mtDNA workflow proved to be a rapid, efficient, and accurate means of sequence analysis for translational mitogenomics.

Novel B2 mitogenomes from Continental southern Patagonia's Late Holocene: New insights into the peopling of the Southern Cone

OBJECTIVES

The main aim of this study is to discuss the migratory processes and peopling dynamics that shaped the genetic variability of populations during the settlement of the Southern Cone, through the analysis of complete mitogenomes of individuals from southern Patagonia.

MATERIALS AND METHODS

Complete mitogenomes were sequenced through massively parallel sequencing from two late Holocene individuals (SAC 1-1-3 and SAC 1-1-4) buried in the same chenque at Salitroso Lake Basin (Santa Cruz province, Argentina). To evaluate matrilineal phylogenetic affinities with other haplotypes, maximum likelihood and Bayesian phylogenetic reconstructions were performed, as well as a haplotype median-joining network.

RESULTS

The mitogenomes were assigned to haplogroups B2 and B2b, exhibiting an average depth of 54X and 89X (≥1X coverage of 98.6% and 100%), and a high number of nucleotide differences among them. The phylogenetic analyses showed a relatively close relationship between the haplotype found in SAC 1-1-4 and those retrieved from a Middle Holocene individual from Laguna Chica (Buenos Aires province), and from a group of individuals from the Peruvian coast. For the SAC 1-1-3, no clear affiliations to any other haplotype were established.

DISCUSSION

The large divergence between the haplotypes presented in this study suggests either a highly variable founder gene pool, or a later enrichment by frequent biological contact with other populations. Our results underline the persistence of genetic signals related to the first waves of peopling in South America, suggesting that the regional settlement of the southern end of the continent has been much more complex than initially thought.

The genetic legacy of African Americans from Catoctin Furnace

Few African Americans have been able to trace family lineages back to ancestors who died before the 1870 United States Census, the first in which all Black people were listed by name. We analyzed 27 individuals from Maryland's Catoctin Furnace African American Cemetery (1774-1850), identifying 41,799 genetic relatives among consenting research participants in 23andMe, Inc.'s genetic database. One of the highest concentrations of close relatives is in Maryland, suggesting that descendants of the Catoctin individuals remain in the area. We find that many of the Catoctin individuals derived African ancestry from the Wolof or Kongo groups and European ancestry from Great Britain and Ireland. This study demonstrates the power of joint analysis of historical DNA and large datasets generated through direct-to-consumer ancestry testing.

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.

The first maternal genetic study of hunter-gatherers from Vietnam

The current limitation of ancient DNA data from Vietnam led to the controversy surrounding the prehistory of people in this region. The combination of high heat and humidity damaged ancient bones that challenged the study of human evolution, especially when using DNA as study materials. So far, only 4 k years of history have been recorded despite the 65 k years of history of anatomically modern human occupations in Vietnam. Here we report, to our knowledge, the oldest mitogenomes of two hunter-gatherers from Vietnam. We extracted DNA from the femurs of two individuals aged 6.2 k cal BP from the Con Co Ngua (CCN) site in Thanh Hoa, Vietnam. This archeological site is the largest cemetery of the hunter-gatherer population in Southeast Asia (SEA) that was discovered, but their genetics have not been explored until the present. We indicated that the CCN haplotype belongs to a rare haplogroup that was not detected in any present-day Vietnamese individuals. Further matrilineal analysis on CCN mitogenomes showed a close relationship with ancient farmers and present-day populations in SEA. The mitogenomes of hunter-gatherers from Vietnam debate the "two layers" model of peopling history in SEA and provide an alternative solution for studying challenging ancient human samples from Vietnam.

The role of Beringia in human adaptation to Arctic conditions based on results of genomic studies of modern and ancient populations

The results of studies in Quaternary geology, archeology, paleoanthropology and human genetics demonstrate that the ancestors of Native Americans arrived in mid-latitude North America mainly along the Pacific Northwest Coast, but had previously inhabited the Arctic and during the last glacial maximum were in a refugium in Beringia, a land bridge connecting Eurasia and North America. The gene pool of Native Americans is represented by unique haplogroups of mitochondrial DNA and the Y chromosome, the evolutionary age of which ranges from 13 to 22 thousand years. The results of a paleogenomic analysis also show that during the last glacial maximum Beringia was populated by human groups that had arisen as a result of interaction between the most ancient Upper Paleolithic populations of Northern Eurasia and newcomer groups from East Asia. Approximately 20 thousand years ago the Beringian populations began to form, and the duration of their existence in relative isolation is estimated at about 5 thousand years. Thus, the adaptation of the Beringians to the Arctic conditions could have taken several millennia. The adaptation of Amerindian ancestors to high latitudes and cold climates is supported by genomic data showing that adaptive genetic variants in Native Americans are associated with various metabolic pathways: melanin production processes in the skin, hair and eyes, the functioning of the cardiovascular system, energy metabolism and immune response characteristics. Meanwhile, the analysis of the existing hypotheses about the selection of some genetic variants in the Beringian ancestors of the Amerindians in connection with adaptation to the Arctic conditions (for example, in the FADS, ACTN3, EDAR genes) shows the ambiguity of the testing results, which may be due to the loss of some traces of the "Beringian" adaptation in the gene pools of modern Native Americans. The most optimal strategy for further research seems to be the search for adaptive variant.

Maternal Ancestry of First Parsi Settlers of India Using Ancient Mitogenome

The rich cultural and genetic diversity of South Asia emerged from multiple migrations and cultural assimilation of multiple waves of migrants. The Parsi community of North-western India were one of those who migrated from West Eurasia in the aftermath of 7^th century CE and assimilated into the local cultural framework. Earlier genetic studies further strengthened this notion with the finding that they harbour both Middle Eastern and South Asian genetic components. Although these studies covered both autosomal and uniparental markers, still maternal ancestry was not covered in depth and with good resolution of mitochondrial markers. Hence in our current study, we have first time generated a complete mitogenome of 19 ancient samples of the first Parsi settlers excavated from the archaeological site of Sanjan and performed detailed phylogenetic analysis to infer their maternal genetic affinity. In our analysis, we found that the Parsi mitogenome with mtDNA haplogroup M3a1+204 shares clade with both Middle Eastern and South Asian modern individuals in both the Maximum Likelihood tree and Bayesian phylogenetic tree. This haplogroup was also prevalent among the medieval Swat valley population of present-day Northern Pakistan and was also observed in two Roopkund A individuals. In the phylogenetic network this sample share haplotype with both South Asian and Middle Eastern samples. So conclusively, the first Parsi settlers' maternal ancestry encompasses both South Asian and Middle Eastern genetic composition.

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.

Mitogenome evidence shows two radiation events and dispersals of matrilineal ancestry from northern coastal China to the Americas and Japan

Although it is widely recognized that the ancestors of Native Americans (NAs) primarily came from Siberia, the link between mitochondrial DNA (mtDNA) lineage D4h3a (typical of NAs) and D4h3b (found so far only in East China and Thailand) raises the possibility that the ancestral sources for early NAs were more variegated than hypothesized. Here, we analyze 216 contemporary (including 106 newly sequenced) D4h mitogenomes and 39 previously reported ancient D4h data. The results reveal two radiation events of D4h in northern coastal China, one during the Last Glacial Maximum and the other within the last deglaciation, which facilitated the dispersals of D4h sub-branches to different areas including the Americas and the Japanese archipelago. The coastal distributions of the NA (D4h3a) and Japanese lineages (D4h1a and D4h2), in combination with the Paleolithic archaeological similarities among Northern China, the Americas, and Japan, lend support to the coastal dispersal scenario of early NAs.

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.

Entwined African and Asian genetic roots of medieval peoples of the Swahili coast

The urban peoples of the Swahili coast traded across eastern Africa and the Indian Ocean and were among the first practitioners of Islam among sub-Saharan people1,2. The extent to which these early interactions between Africans and non-Africans were accompanied by genetic exchange remains unknown. Here we report ancient DNA data for 80 individuals from 6 medieval and early modern (AD 1250-1800) coastal towns and an inland town after AD 1650. More than half of the DNA of many of the individuals from coastal towns originates from primarily female ancestors from Africa, with a large proportion-and occasionally more than half-of the DNA coming from Asian ancestors. The Asian ancestry includes components associated with Persia and India, with 80-90% of the Asian DNA originating from Persian men. Peoples of African and Asian origins began to mix by about AD 1000, coinciding with the large-scale adoption of Islam. Before about AD 1500, the Southwest Asian ancestry was mainly Persian-related, consistent with the narrative of the Kilwa Chronicle, the oldest history told by people of the Swahili coast3. After this time, the sources of DNA became increasingly Arabian, consistent with evidence of growing interactions with southern Arabia4. Subsequent interactions with Asian and African people further changed the ancestry of present-day people of the Swahili coast in relation to the medieval individuals whose DNA we sequenced.

Phylogeographic history of mitochondrial haplogroup J in Scandinavia

BACKGROUND

Mitochondrial DNA haplogroup J is the third most frequent haplogroup in modern-day Scandinavia, although it did not originate there. To infer the genetic history of haplogroup J in Scandinavia, we examined worldwide mitogenome sequences using a maximum-likelihood phylogenetic approach.

METHODS

Haplogroup J mitogenome sequences were gathered from GenBank (n = 2245) and aligned against the ancestral Reconstructed Sapiens Reference Sequence. We also analyzed haplogroup J Viking Age sequences from the European Nucleotide Archive (n = 54). Genetic distances were estimated from these data and projected onto a maximum likelihood rooted phylogenetic tree to analyze clustering and branching dates.

RESULTS

Haplogroup J originated approximately 42.6 kya (95% CI: 30.0-64.7), with several of its earliest branches being found within the Arabian Peninsula and Northern Africa. J1b was found most frequently in the Near East and Arabian Peninsula, while J1c occurred most frequently in Europe. Based on phylogenetic dating, subhaplogroup J1c has its early roots in the Mediterranean and Western Balkans. Otherwise, the majority of the branches found in Scandinavia are younger than those seen elsewhere, indicating that haplogroup J dispersed relatively recently into Northern Europe, most plausibly with Neolithic farmers.

CONCLUSIONS

Haplogroup J appeared when Scandinavia was transitioning to agriculture over 6 kya, with J1c being the most common lineage there today. Changes in the distribution of haplogroup J mtDNAs were likely driven by the expansion of farming from West Asia into Southern Europe, followed by a later expansion into Scandinavia, with other J subhaplogroups appearing among Scandinavian groups as early as the Viking Age.

Complete mitochondrial DNA sequence tries to settle hitherto putative history of Kayastha population of India

OBJECTIVES

Indian caste system is unique as it has an inimitable type of class system where the social ordering is done based on birth. Within the caste system, there is a distinct endogamous population known as the Kayastha, who have had inconclusive stratification records due to unidentified historical records.

METHODS

To gain a more inclusive view on the history and genetic affinities of Kayastha people, complete mitochondrial genomes from 15 individuals of a Kayastha population from North-western India have been sequenced.

RESULTS

Interestingly, three novel sub-clades (U2b2a, M3d2, and M33a3b) have been identified that represent unique Kayastha motifs.

CONCLUSION

The haplotype-based analysis suggests that the Kayastha population shares genetic affinities with the Indo-European and Sino-Tibetan populations found in the trans-Himalayan region. The F_ST based population comparison and the MDS plot indicates that Kayastha people have close maternal genetic affinity with the available genetic database of Brahmins, Kashmiris, and Tharus. The maternal genetic lineages among Kayastha population shows deep in situ origin that emerged much before settled life developed on this sub-continent. Both mtDNA and Y-chromosome markers, trace the genetic lineages of Kayastha population with Tharus, who regard themselves Kshatriya, corroborated by the oral history of the Kayasthas for their Kshatriya affiliation. It also validates genetic heritage of earliest settlers of India in both indigenous tribal and caste populations.

Genome-wide data from medieval German Jews show that the Ashkenazi founder event pre-dated the 14th century

We report genome-wide data from 33 Ashkenazi Jews (AJ), dated to the 14^th century, obtained following a salvage excavation at the medieval Jewish cemetery of Erfurt, Germany. The Erfurt individuals are genetically similar to modern AJ, but they show more variability in Eastern European-related ancestry than modern AJ. A third of the Erfurt individuals carried a mitochondrial lineage common in modern AJ and eight carried pathogenic variants known to affect AJ today. These observations, together with high levels of runs of homozygosity, suggest that the Erfurt community had already experienced the major reduction in size that affected modern AJ. The Erfurt bottleneck was more severe, implying substructure in medieval AJ. Overall, our results suggest that the AJ founder event and the acquisition of the main sources of ancestry pre-dated the 14^th century and highlight late medieval genetic heterogeneity no longer present in modern AJ.

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.

Genomes from a medieval mass burial show Ashkenazi-associated hereditary diseases pre-date the 12th century

We report genome sequence data from six individuals excavated from the base of a medieval well at a site in Norwich, UK. A revised radiocarbon analysis of the assemblage is consistent with these individuals being part of a historically attested episode of antisemitic violence on 6 February 1190 CE. We find that four of these individuals were closely related and all six have strong genetic affinities with modern Ashkenazi Jews. We identify four alleles associated with genetic disease in Ashkenazi Jewish populations and infer variation in pigmentation traits, including the presence of red hair. Simulations indicate that Ashkenazi-associated genetic disease alleles were already at appreciable frequencies, centuries earlier than previously hypothesized. These findings provide new insights into a significant historical crime, into Ashkenazi population history, and into the origins of genetic diseases associated with modern Jewish populations.

The diverse genetic origins of a Classical period Greek army

Trade and colonization caused an unprecedented increase in Mediterranean human mobility in the first millennium BCE. Often seen as a dividing force, warfare is in fact another catalyst of culture contact. We provide insight into the demographic dynamics of ancient warfare by reporting genome-wide data from fifth-century soldiers who fought for the army of the Greek Sicilian colony of Himera, along with representatives of the civilian population, nearby indigenous settlements, and 96 present-day individuals from Italy and Greece. Unlike the rest of the sample, many soldiers had ancestral origins in northern Europe, the Steppe, and the Caucasus. Integrating genetic, archaeological, isotopic, and historical data, these results illustrate the significant role mercenaries played in ancient Greek armies and highlight how participation in war contributed to continental-scale human mobility in the Classical world.

The Anglo-Saxon migration and the formation of the early English gene pool

The history of the British Isles and Ireland is characterized by multiple periods of major cultural change, including the influential transformation after the end of Roman rule, which precipitated shifts in language, settlement patterns and material culture1. The extent to which migration from continental Europe mediated these transitions is a matter of long-standing debate2-4. Here we study genome-wide ancient DNA from 460 medieval northwestern Europeans-including 278 individuals from England-alongside archaeological data, to infer contemporary population dynamics. We identify a substantial increase of continental northern European ancestry in early medieval England, which is closely related to the early medieval and present-day inhabitants of Germany and Denmark, implying large-scale substantial migration across the North Sea into Britain during the Early Middle Ages. As a result, the individuals who we analysed from eastern England derived up to 76% of their ancestry from the continental North Sea zone, albeit with substantial regional variation and heterogeneity within sites. We show that women with immigrant ancestry were more often furnished with grave goods than women with local ancestry, whereas men with weapons were as likely not to be of immigrant ancestry. A comparison with present-day Britain indicates that subsequent demographic events reduced the fraction of continental northern European ancestry while introducing further ancestry components into the English gene pool, including substantial southwestern European ancestry most closely related to that seen in Iron Age France5,6.

Genetic structure and differentiation from early bronze age in the mediterranean island of sicily: Insights from ancient mitochondrial genomes

Sicily is one of the main islands of the Mediterranean Sea, and it is characterized by a variety of archaeological records, material culture and traditions, reflecting the history of migrations and populations’ interaction since its first colonization, during the Paleolithic. These deep and complex demographic and cultural dynamics should have affected the genomic landscape of Sicily at different levels; however, the relative impact of these migrations on the genomic structure and differentiation within the island remains largely unknown. The available Sicilian modern genetic data gave a picture of the current genetic structure, but the paucity of ancient data did not allow so far to make predictions about the level of historical variation. In this work, we sequenced and analyzed the complete mitochondrial genomes of 36 individuals from five different locations in Sicily, spanning from Early Bronze Age to Iron Age, and with different cultural backgrounds. The comparison with coeval groups from the Mediterranean Basin highlighted structured genetic variation in Sicily since Early Bronze Age, thus supporting a demic impact of the cultural transitions within the Island. Explicit model testing through Approximate Bayesian Computation allowed us to make predictions about the origin of Sicanians, one of the three indigenous peoples of Sicily, whose foreign origin from Spain, historically attributed, was not confirmed by our analysis of genetic data. Sicilian modern mitochondrial data show a different, more homogeneous, genetic composition, calling for a recent genetic replacement in the Island of pre-Iron Age populations, that should be further investigated.

The genetic history of the Southern Arc: A bridge between West Asia and Europe

By sequencing 727 ancient individuals from the Southern Arc (Anatolia and its neighbors in Southeastern Europe and West Asia) over 10,000 years, we contextualize its Chalcolithic period and Bronze Age (about 5000 to 1000 BCE), when extensive gene flow entangled it with the Eurasian steppe. Two streams of migration transmitted Caucasus and Anatolian/Levantine ancestry northward, and the Yamnaya pastoralists, formed on the steppe, then spread southward into the Balkans and across the Caucasus into Armenia, where they left numerous patrilineal descendants. Anatolia was transformed by intra-West Asian gene flow, with negligible impact of the later Yamnaya migrations. This contrasts with all other regions where Indo-European languages were spoken, suggesting that the homeland of the Indo-Anatolian language family was in West Asia, with only secondary dispersals of non-Anatolian Indo-Europeans from the steppe.

Sequence analyses of Malaysian Indigenous communities reveal historical admixture between Hoabinhian hunter-gatherers and Neolithic farmers

Southeast Asia comprises 11 countries that span mainland Asia across to numerous islands that stretch from the Andaman Sea to the South China Sea and Indian Ocean. This region harbors an impressive diversity of history, culture, religion and biology. Indigenous people of Malaysia display substantial phenotypic, linguistic, and anthropological diversity. Despite this remarkable diversity which has been documented for centuries, the genetic history and structure of indigenous Malaysians remain under-studied. To have a better understanding about the genetic history of these people, especially Malaysian Negritos, we sequenced whole genomes of 15 individuals belonging to five indigenous groups from Peninsular Malaysia and one from North Borneo to high coverage (30X). Our results demonstrate that indigenous populations of Malaysia are genetically close to East Asian populations. We show that present-day Malaysian Negritos can be modeled as an admixture of ancient Hoabinhian hunter-gatherers and Neolithic farmers. We observe gene flow from South Asian populations into the Malaysian indigenous groups, but not into Dusun of North Borneo. Our study proposes that Malaysian indigenous people originated from at least three distinct ancestral populations related to the Hoabinhian hunter-gatherers, Neolithic farmers and Austronesian speakers.

Ancient DNA reveals five streams of migration into Micronesia and matrilocality in early Pacific seafarers

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

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.

The Mitogenome Relationships and Phylogeography of Barn Swallows (Hirundo rustica)

The barn swallow (Hirundo rustica) poses a number of fascinating scientific questions, including the taxonomic status of postulated subspecies. Here we obtained and assessed the sequence variation of 411 complete mitogenomes, mainly from the European H. r. rustica, but other subspecies as well. In almost every case, we observed subspecies-specific haplogroups, which we employed together with estimated radiation times to postulate a model for the geographical and temporal worldwide spread of the species. The female barn swallow carrying the Hirundo rustica ancestral mitogenome left Africa (or its vicinity) around 280 thousand years ago (kya), and her descendants expanded first into Eurasia and then, at least 51 kya, into the Americas, from where a relatively recent (< 20 kya) back migration to Asia took place. The exception to the haplogroup subspecies specificity is represented by the sedentary Levantine H. r. transitiva that extensively shares haplogroup A with the migratory European H. r. rustica and, to a lesser extent, haplogroup B with the Egyptian H. r. savignii. Our data indicate that rustica and transitiva most likely derive from a sedentary Levantine population source that split at the end of the Younger Dryas (11.7 kya). Since then, however, transitiva received genetic inputs from and admixed with both the closely related rustica and the adjacent savignii. Demographic analyses confirm this species' strong link with climate fluctuations and human activities making it an excellent indicator for monitoring and assessing the impact of current global changes on wildlife.

Evolution and dispersal of mitochondrial DNA haplogroup U5 in Northern Europe: insights from an unsupervised learning approach to phylogeography

Background

We combined an unsupervised learning methodology for analyzing mitogenome sequences with maximum likelihood (ML) phylogenetics to make detailed inferences about the evolution and diversification of mitochondrial DNA (mtDNA) haplogroup U5, which appears at high frequencies in northern Europe.

Methods

Haplogroup U5 mitogenome sequences were gathered from GenBank. The hierarchal Bayesian Analysis of Population Structure (hierBAPS) method was used to generate groups of sequences that were then projected onto a rooted maximum likelihood (ML) phylogenetic tree to visualize the pattern of clustering. The haplogroup statuses of the individual sequences were assessed using Haplogrep2.

Results

A total of 23 hierBAPS groups were identified, all of which corresponded to subclades defined in Phylotree, v.17. The hierBAPS groups projected onto the ML phylogeny accurately clustered all haplotypes belonging to a specific haplogroup in accordance with Haplogrep2. By incorporating the geographic source of each sequence and subclade age estimates into this framework, inferences about the diversification of U5 mtDNAs were made. Haplogroup U5 has been present in northern Europe since the Mesolithic, and spread in both eastern and western directions, undergoing significant diversification within Scandinavia. A review of historical and archeological evidence attests to some of the population interactions contributing to this pattern.

Conclusions

The hierBAPS algorithm accurately grouped mitogenome sequences into subclades in a phylogenetically robust manner. This analysis provided new insights into the phylogeographic structure of haplogroup U5 diversity in northern Europe, revealing a detailed perspective on the diversity of subclades in this region and their distribution in Scandinavian populations.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08572-y.

Neither cardiac mitochondrial DNA variation nor copy number contribute to congenital heart disease risk

The well-established manifestation of mitochondrial mutations in functional cardiac disease (e.g., mitochondrial cardiomyopathy) prompted the hypothesis that mitochondrial DNA (mtDNA) sequence and/or copy number (mtDNAcn) variation contribute to cardiac defects in congenital heart disease (CHD). MtDNAcns were calculated and rare, non-synonymous mtDNA mutations were identified in 1,837 CHD-affected proband-parent trios, 116 CHD-affected singletons, and 114 paired cardiovascular tissue/blood samples. The variant allele fraction (VAF) of heteroplasmic variants in mitochondrial RNA from 257 CHD cardiovascular tissue samples was also calculated. On average, mtDNA from blood had 0.14 rare variants and 52.9 mtDNA copies per nuclear genome per proband. No variation with parental age at proband birth or CHD-affected proband age was seen. mtDNAcns in valve/vessel tissue (320 ± 70) were lower than in atrial tissue (1,080 ± 320, p = 6.8E-21), which were lower than in ventricle tissue (1,340 ± 280, p = 1.4E-4). The frequency of rare variants in CHD-affected individual DNA was indistinguishable from the frequency in an unaffected cohort, and proband mtDNAcns did not vary from those of CHD cohort parents. In both the CHD and the comparison cohorts, mtDNAcns were significantly correlated between mother-child, father-child, and mother-father. mtDNAcns among people with European (mean = 52.0), African (53.0), and Asian haplogroups (53.5) were calculated and were significantly different for European and Asian haplogroups (p = 2.6E-3). Variant heteroplasmic fraction (HF) in blood correlated well with paired cardiovascular tissue HF (r = 0.975) and RNA VAF (r = 0.953), which suggests blood HF is a reasonable proxy for HF in heart tissue. We conclude that mtDNA mutations and mtDNAcns are unlikely to contribute significantly to CHD risk.

Mitochondrial 1555 G&gt;A variant as a potential risk factor for childhood glioblastoma

Background

Childhood glioblastoma multiforme (GBM) is a highly aggressive disease with low survival, and its etiology, especially concerning germline genetic risk, is poorly understood. Mitochondria play a key role in putative tumorigenic processes relating to cellular oxidative metabolism, and mitochondrial DNA variants were not previously assessed for association with pediatric brain tumor risk.

Methods

We conducted an analysis of 675 mitochondrial DNA variants in 90 childhood GBM cases and 2789 controls to identify enrichment of mitochondrial variant associated with GBM risk. We also performed this analysis for other glioma subtypes including pilocytic astrocytoma. Nuclear-encoded mitochondrial gene variants were also analyzed.

Results

We identified m1555 A>G was significantly associated with GBM risk (adjusted OR 29.30, 95% CI 5.25–163.4, P-value 9.5 X 10^–4). No association was detected for other subtypes. Haplotype analysis further supported the independent risk contributed by m1555 G>A, instead of a haplogroup joint effect. Nuclear-encoded mitochondrial gene variants identified significant associations in European (rs62036057 in WWOX, adjusted OR = 2.99, 95% CI 1.88–4.75, P-value = 3.42 X 10^–6) and Hispanic (rs111709726 in EFHD1, adjusted OR = 3.57, 95% CI 1.99–6.40, P-value = 1.41 X 10^–6) populations in ethnicity-stratified analyses.

Conclusion

We report for the first time a potential role played by a functional mitochondrial ribosomal RNA variant in childhood GBM risk, and a potential role for both mitochondrial and nuclear-mitochondrial DNA polymorphisms in GBM tumorigenesis. These data implicate cellular oxidative metabolic capacity as a contributor to the etiology of pediatric glioblastoma.

South-to-north migration preceded the advent of intensive farming in the Maya region

The genetic prehistory of human populations in Central America is largely unexplored leaving an important gap in our knowledge of the global expansion of humans. We report genome-wide ancient DNA data for a transect of twenty individuals from two Belize rock-shelters dating between 9,600-3,700 calibrated radiocarbon years before present (cal. BP). The oldest individuals (9,600-7,300 cal. BP) descend from an Early Holocene Native American lineage with only distant relatedness to present-day Mesoamericans, including Mayan-speaking populations. After ~5,600 cal. BP a previously unknown human dispersal from the south made a major demographic impact on the region, contributing more than 50% of the ancestry of all later individuals. This new ancestry derived from a source related to present-day Chibchan speakers living from Costa Rica to Colombia. Its arrival corresponds to the first clear evidence for forest clearing and maize horticulture in what later became the Maya region.

The genetic prehistory of central America has not been well explored. Here, the authors find evidence from ancient DNA from twenty individuals who lived in Belize 9,600 to 3,700 years ago of a migration from the south that coincided with the first evidence for forest clearing and the spread of maize horticulture.

Ancient DNA and deep population structure in sub-Saharan African foragers

Multiple lines of genetic and archaeological evidence suggest that there were major demographic changes in the terminal Late Pleistocene epoch and early Holocene epoch of sub-Saharan Africa^1–4. Inferences about this period are challenging to make because demographic shifts in the past 5,000 years have obscured the structures of more ancient populations^3,5. Here we present genome-wide ancient DNA data for six individuals from eastern and south-central Africa spanning the past approximately 18,000 years (doubling the time depth of sub-Saharan African ancient DNA), increase the data quality for 15 previously published ancient individuals and analyse these alongside data from 13 other published ancient individuals. The ancestry of the individuals in our study area can be modelled as a geographically structured mixture of three highly divergent source populations, probably reflecting Pleistocene interactions around 80–20 thousand years ago, including deeply diverged eastern and southern African lineages, plus a previously unappreciated ubiquitous distribution of ancestry that occurs in highest proportion today in central African rainforest hunter-gatherers. Once established, this structure remained highly stable, with limited long-range gene flow. These results provide a new line of genetic evidence in support of hypotheses that have emerged from archaeological analyses but remain contested, suggesting increasing regionalization at the end of the Pleistocene epoch.

DNA analysis of 6 individuals from eastern and south-central Africa spanning the past approximately 18,000 years, and of 28 previously published ancient individuals, provides genetic evidence supporting hypotheses of increasing regionalization at the end of the Pleistocene.

Ancient DNA at the edge of the world: Continental immigration and the persistence of Neolithic male lineages in Bronze Age Orkney

The Orcadian Neolithic has been intensively studied and celebrated as a major center of cultural innovation, whereas the Bronze Age is less well known and often regarded as a time of stagnation and insularity. Here, we analyze ancient genomes from the Orcadian Bronze Age in the context of the variation in Neolithic Orkney and Bronze Age Europe. We find clear evidence for Early Bronze Age immigration into Orkney, but with an extraordinary pattern: continuity from the Neolithic on the male line of descent but immigration from continental Europe on the female side, echoed in the genome-wide picture. This suggests that despite substantial immigration, indigenous male lineages persisted for at least a thousand years after the end of the Neolithic.

Orkney was a major cultural center during the Neolithic, 3800 to 2500 BC. Farming flourished, permanent stone settlements and chambered tombs were constructed, and long-range contacts were sustained. From ∼3200 BC, the number, density, and extravagance of settlements increased, and new ceremonial monuments and ceramic styles, possibly originating in Orkney, spread across Britain and Ireland. By ∼2800 BC, this phenomenon was waning, although Neolithic traditions persisted to at least 2500 BC. Unlike elsewhere in Britain, there is little material evidence to suggest a Beaker presence, suggesting that Orkney may have developed along an insular trajectory during the second millennium BC. We tested this by comparing new genomic evidence from 22 Bronze Age and 3 Iron Age burials in northwest Orkney with Neolithic burials from across the archipelago. We identified signals of inward migration on a scale unsuspected from the archaeological record: As elsewhere in Bronze Age Britain, much of the population displayed significant genome-wide ancestry deriving ultimately from the Pontic-Caspian Steppe. However, uniquely in northern and central Europe, most of the male lineages were inherited from the local Neolithic. This suggests that some male descendants of Neolithic Orkney may have remained distinct well into the Bronze Age, although there are signs that this had dwindled by the Iron Age. Furthermore, although the majority of mitochondrial DNA lineages evidently arrived afresh with the Bronze Age, we also find evidence for continuity in the female line of descent from Mesolithic Britain into the Bronze Age and even to the present day.

An invasive Haemophilus influenzae serotype b infection in an Anglo-Saxon plague victim

Background

The human pathogen Haemophilus influenzae was the main cause of bacterial meningitis in children and a major cause of worldwide infant mortality before the introduction of a vaccine in the 1980s. Although the occurrence of serotype b (Hib), the most virulent type of H. influenzae, has since decreased, reports of infections with other serotypes and non-typeable strains are on the rise. While non-typeable strains have been studied in-depth, very little is known of the pathogen’s evolutionary history, and no genomes dating prior to 1940 were available.

Results

We describe a Hib genome isolated from a 6-year-old Anglo-Saxon plague victim, from approximately 540 to 550 CE, Edix Hill, England, showing signs of invasive infection on its skeleton. We find that the genome clusters in phylogenetic division II with Hib strain NCTC8468, which also caused invasive disease. While the virulence profile of our genome was distinct, its genomic similarity to NCTC8468 points to mostly clonal evolution of the clade since the 6th century. We also reconstruct a partial Yersinia pestis genome, which is likely identical to a published first plague pandemic genome of Edix Hill.

Conclusions

Our study presents the earliest genomic evidence for H. influenzae, points to the potential presence of larger genomic diversity in the phylogenetic division II serotype b clade in the past, and allows the first insights into the evolutionary history of this major human pathogen. The identification of both plague and Hib opens questions on the effect of plague in immunocompromised individuals already affected by infectious diseases.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-021-02580-z.

Large-scale migration into Britain during the Middle to Late Bronze Age

Present-day people from England and Wales harbour more ancestry derived from Early European Farmers (EEF) than people of the Early Bronze Age¹. To understand this, we generated genome-wide data from 793 individuals, increasing data from the Middle to Late Bronze and Iron Age in Britain by 12-fold, and Western and Central Europe by 3.5-fold. Between 1000 and 875 BC, EEF ancestry increased in southern Britain (England and Wales) but not northern Britain (Scotland) due to incorporation of migrants who arrived at this time and over previous centuries, and who were genetically most similar to ancient individuals from France. These migrants contributed about half the ancestry of Iron Age people of England and Wales, thereby creating a plausible vector for the spread of early Celtic languages into Britain. These patterns are part of a broader trend of EEF ancestry becoming more similar across central and western Europe in the Middle to Late Bronze Age, coincident with archaeological evidence of intensified cultural exchange^2-6. There was comparatively less gene flow from continental Europe during the Iron Age, and Britain's independent genetic trajectory is also reflected in the rise of the allele conferring lactase persistence to ~50% by this time compared to ~7% in central Europe where it rose rapidly in frequency only a millennium later. This suggests that dairy products were used in qualitatively different ways in Britain and in central Europe over this period.