The dawn of human matrilineal diversity

Deciphering the maternal ancestral lineage of Greek Cypriots, Armenian Cypriots and Maronite Cypriots

Cyprus was conquered from several populations because of its special geographical location. In this study, 406 unrelated Cypriot samples were tested based on their mitochondrial DNA. In more detail, 185 were Greek Cypriots, 114 Armenian Cypriots and 107 Maronite Cypriots. This is the first time where the mitochondrial DNA of Greek Cypriots, Armenian Cypriots and Maronite Cypriots is compared with the aim of characterizing the maternal ancestry of Cypriots. The control region of the mtDNA is the most informative in terms of studying maternal ancestry and consists of three hypervariable regions (HVS-I, HVS-II, HVS-III). The hypervariable regions can provide important information regarding the maternal ancestor of the tested samples. The entire control region of the mtDNA was used to determine the mitotypes and subsequently the haplogroups of all the Cypriot DNA samples. Based on the aforementioned analyses, Greek Cypriots were found to be genetically closer to Armenian Cypriots, while Greek Cypriots and Armenian Cypriots showed moderate genetic differentiation with Maronite Cypriots. The most prevalent haplogroups among Cypriots were haplogroups H and U, while R0 is common but in different frequencies for Greek Cypriots, Armenian Cypriots and Maronite Cypriots. It is proposed that the maternal ancestor may have originated during the Neolithic period and/or the Bronze age.

Refining the Global Phylogeny of Mitochondrial N1a, X, and HV2 Haplogroups Based on Rare Mitogenomes from Croatian Isolates

Mitochondrial DNA (mtDNA) has been used for decades as a predominant tool in population genetics and as a valuable addition to forensic genetic research, owing to its unique maternal inheritance pattern that enables the tracing of individuals along the maternal lineage across numerous generations. The dynamic interplay between evolutionary forces, primarily genetic drift, bottlenecks, and the founder effect, can exert significant influence on genetic profiles. Consequently, the Adriatic islands have accumulated a subset of lineages that exhibits remarkable absence or rarity within other European populations. This distinctive genetic composition underscores the islands' potential as a significant resource in phylogenetic research, with implications reaching beyond regional boundaries to contribute to a global understanding. In the initial attempt to expand the mitochondrial forensic database of the Croatian population with haplotypes from small isolated communities, we sequenced mitogenomes of rare haplogroups from different Croatian island and mainland populations using next-generation sequencing (NGS). In the next step and based on the obtained results, we refined the global phylogeny of haplogroup N1a, HV2, and X by analyzing rare haplotypes, which are absent from the current phylogenetic tree. The trees were based on 16 novel and 52 previously published samples, revealing completely novel branches in the X and HV2 haplogroups and a new European cluster in the ancestral N1a variant, previously believed to be an exclusively African-Asian haplogroup. The research emphasizes the importance of investigating geographically isolated populations and their unique characteristics within a global context.

Mitochondrial DNA haplogroup, genetic ancestry, and susceptibility to Ewing sarcoma

Based on current studies, the incidence of Ewing sarcoma (ES) varies significantly by race and ethnicity, with the disease being most common in patients of European ancestry. However, race/ethnicity has generally been self-reported rather than formally evaluated at a population level using DNA evidence. Additionally, mitochondrial dysfunction is a hallmark of ES, yet there have been no reported studies of mitochondrial genetics in ES. Thus, we evaluated both the mitochondrial and nuclear ancestries of 420 pediatric ES patients in the United States using whole-genome sequencing. We found that the mitochondrial DNA (mtDNA) genomes of only six (1.4 %) patients belonged to African L haplogroups, while those of 90 % of the patients belonged to macrohaplogroup R, which includes haplogroup H, the most common maternal lineage in Europe. Compared to the general US population, European haplogroups were significantly enriched in ES patients (p < 2.2e-16) and the African haplogroups are significantly impoverished (p < 4.6e-16). Using the ancestry informative markers defined in a National Genographic study, the vast majority of patients exhibited significant nuclear ancestry originating from the Mediterranean, Northern Europe, and Southwest Asia, including all six patients with African L mtDNAs. Very few had primarily African nuclear ancestry. This is the first genomic epidemiology study to simultaneously interrogate the mitochondrial and nuclear ancestries of ES patients. While supporting previous findings of enriched European ancestry in ES patients, these results also suggest alternative hypotheses for the significant contribution of mitochondrial ancestry in ES patients, as well as the protective role of African ancestry.

Whole mitochondrial genome analysis of Tai-Kadai-speaking populations in Southwest China

As a single matrilineal gene, human mitochondrial DNA plays a very important role in the study of population genetics. The whole mitogenome sequences of 287 individuals of the Tai-Kadai-speaking population in Guizhou were obtained. It was discovered that there were 82, 104, and 94 haplotypes in 83 Bouyei individuals, 107 Dong individuals, and 97 Sui individuals, respectively; and the haplotype diversity in Bouyei, Dong, and Sui groups was 1.000 ± 0.02, 0.9993 ± 0.0015, and 0.999 ± 0.002, respectively. The result of neutrality tests of the Tai-Kadai-speaking population in Guizhou showed significant negative values, and the analysis of mismatch distribution showed an obvious unimodal distribution. The results implied that Guizhou Tai-Kadai-speaking populations had high genetic diversities and may have experienced recent population expansion. In addition, the primary haplogroups of studied populations were M*, F, B, D, and R*, implying that they may origin from Southern China. The matrilineal genetic structure of the Tai-Kadai-speaking populations in Guizhou was analyzed by merging the mitogenome data of 79 worldwide populations as reference data. The results showed that there were close relationships between studied populations and other Tai-Kadai as well as some Austronesian populations in East and Southeast Asia. Overall, the mitogenome data generated in this study will provide important data for the study of genetic structure of Tai-Kadai speaking populations.

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.

Demographic history was a formative mechanism of the genetic structure for the taste receptor TAS2R16 in human populations inhabiting Africa's Sahel/Savannah Belt

OBJECTIVES

Mode of subsistence is an important factor influencing dietary habits and the genetic structure of various populations through differential intensity of gene flow and selection pressures. Previous studies suggest that in Africa Taste 2 Receptor Member 16 (TAS2R16), which encodes the 7-transmembrane receptor protein for bitterness, might also be under positive selection pressure.

METHODS

However, since sampling coverage of populations was limited, we created a new TAS2R16 population dataset from across the African Sahel/Savannah belt representing various local populations of differing subsistence modes, linguistic affiliations, and geographic provenience. We sequenced the TAS2R16 exon gene and analyzed 2250 haplotypes among 19 populations.

RESULTS

We found no evidence for selection as a driving force of genetic variation at this locus; instead, we discovered a highly significant correlation between TAS2R16 genetic and geographical distances based on provenience of the sampled populations, strongly suggesting that genetic drift most likely prevailed over positive selection at this specific locus. We also found significant correlations with other independent loci, mainly in sedentary farmers.

DISCUSSION

Our results do not support the notion that the genetic diversity of TAS2R16 in Sahelian populations was shaped by selective pressures. This could result from several alternative and not mutually exclusive mechanisms, of which the possibility that, due to the pleiotropic nature of TAS2R16, selective pressures on other traits could counterbalance those acting on bitter taste perception, or that the change of diet in the Neolithic generally relaxed selective pressure on this gene.

Male-biased migration from East Africa introduced pastoralism into southern Africa

BACKGROUND

Hunter-gatherer lifestyles dominated the southern African landscape up to ~ 2000 years ago, when herding and farming groups started to arrive in the area. First, herding and livestock, likely of East African origin, appeared in southern Africa, preceding the arrival of the large-scale Bantu-speaking agro-pastoralist expansion that introduced West African-related genetic ancestry into the area. Present-day Khoekhoe-speaking Namaqua (or Nama in short) pastoralists show high proportions of East African admixture, linking the East African ancestry with Khoekhoe herders. Most other historical Khoekhoe populations have, however, disappeared over the last few centuries and their contribution to the genetic structure of present-day populations is not well understood. In our study, we analyzed genome-wide autosomal and full mitochondrial data from a population who trace their ancestry to the Khoekhoe-speaking Hessequa herders from the southern Cape region of what is now South Africa.

RESULTS

We generated genome-wide data from 162 individuals and mitochondrial DNA data of a subset of 87 individuals, sampled in the Western Cape Province, South Africa, where the Hessequa population once lived. Using available comparative data from Khoe-speaking and related groups, we aligned genetic date estimates and admixture proportions to the archaeological proposed dates and routes for the arrival of the East African pastoralists in southern Africa. We identified several Afro-Asiatic-speaking pastoralist groups from Ethiopia and Tanzania who share high affinities with the East African ancestry present in southern Africa. We also found that the East African pastoralist expansion was heavily male-biased, akin to a pastoralist migration previously observed on the genetic level in ancient Europe, by which Pontic-Caspian Steppe pastoralist groups represented by the Yamnaya culture spread across the Eurasian continent during the late Neolithic/Bronze Age.

CONCLUSION

We propose that pastoralism in southern Africa arrived through male-biased migration of an East African Afro-Asiatic-related group(s) who introduced new subsistence and livestock practices to local southern African hunter-gatherers. Our results add to the understanding of historical human migration and mobility in Africa, connected to the spread of food-producing and livestock practices.

Mitochondrial genetic variation in human bioenergetics, adaptation, and adult disease

OBJECTIVES

Mitochondria are critical for the survival of eukaryotic organisms due to their ability to produce cellular energy, which drives virtually all aspects of host biology. However, the effects of mitochondrial DNA (mtDNA) variation in relation to disease etiology and adaptation within contemporary global human populations remains incompletely understood.

METHODS

To develop a more holistic understanding of the role of mtDNA diversity in human adaptation, health, and disease, we investigated mitochondrial biology and bioenergetics. More specifically, we synthesized details from studies of mitochondrial function and variation in the context of haplogroup background, climatic adaptation, and oxidative disease.

RESULTS

The majority of studies show that mtDNA variation arose during modern human dispersal around the world. Some of these variants appear to have been positively selected for their adaptiveness in colder climates, with these sequence changes having implications for tissue-specific function and thermogenic capacity. In addition, many variants modulating energy production are also associated with damaging metabolic byproducts and mitochondrial dysfunction, which, in turn, are implicated in the onset and severity of several different adult mitochondrial diseases. Thus, mtDNA variation that governs bioenergetics, metabolism, and thermoregulation may potentially have adverse consequences for human health, depending on the genetic background and context in which it occurs.

CONCLUSIONS

Our review suggests that the mitochondrial research field would benefit from independently replicating mtDNA haplogroup-phenotype associations across global populations, incorporating potentially confounding environmental, demographic, and disease covariates into studies of mtDNA variation, and extending association-based studies to include analyses of complete mitogenomes and assays of mitochondrial function.

The genomic origin of Zana of Abkhazia

Enigmatic phenomena have sparked the imagination of people around the globe into creating folkloric creatures. One prime example is Zana of Abkhazia (South Caucasus), a well-documented 19th century female who was captured living wild in the forest. Zana's appearance was sufficiently unusual, that she was referred to by locals as an Almasty-the analog of Bigfoot in the Caucasus. Although the exact location of Zana's burial site was unknown, the grave of her son, Khwit, was identified in 1971. The genomes of Khwit and the alleged Zana skeleton were sequenced to an average depth of ca. 3× using ancient DNA techniques. The identical mtDNA and parent-offspring relationship between the two indicated that the unknown woman was indeed Zana. Population genomic analyses demonstrated that Zana's immediate genetic ancestry can likely be traced to present-day East-African populations. We speculate that Zana might have had a genetic disorder such as congenital generalized hypertrichosis which could partially explain her strange behavior, lack of speech, and long body hair. Our findings elucidate Zana's unfortunate story and provide a clear example of how prejudices of the time led to notions of cryptic hominids that are still held and transmitted by some today.

The genomic history of the Aegean palatial civilizations

The Cycladic, the Minoan, and the Helladic (Mycenaean) cultures define the Bronze Age (BA) of Greece. Urbanism, complex social structures, craft and agricultural specialization, and the earliest forms of writing characterize this iconic period. We sequenced six Early to Middle BA whole genomes, along with 11 mitochondrial genomes, sampled from the three BA cultures of the Aegean Sea. The Early BA (EBA) genomes are homogeneous and derive most of their ancestry from Neolithic Aegeans, contrary to earlier hypotheses that the Neolithic-EBA cultural transition was due to massive population turnover. EBA Aegeans were shaped by relatively small-scale migration from East of the Aegean, as evidenced by the Caucasus-related ancestry also detected in Anatolians. In contrast, Middle BA (MBA) individuals of northern Greece differ from EBA populations in showing ∼50% Pontic-Caspian Steppe-related ancestry, dated at ca. 2,600-2,000 BCE. Such gene flow events during the MBA contributed toward shaping present-day Greek genomes.

The deep population history in Africa

Africa is the continent with the greatest genetic diversity among humans and the level of diversity is further enhanced by incorporating non-majority groups, which are often understudied. Many of today's minority populations historically practiced foraging lifestyles, which were the only subsistence strategies prior to the rise of agriculture and pastoralism, but only a few groups practicing these strategies remain today. Genomic investigations of Holocene human remains excavated across the African continent show that the genetic landscape was vastly different compared to today's genetic landscape and that many groups that today are population isolate inhabited larger regions in the past. It is becoming clear that there are periods of isolation among groups and geographic areas, but also genetic contact over large distances throughout human history in Africa. Genomic information from minority populations and from prehistoric remains provide an invaluable source of information on the human past, in particular deep human population history, as Holocene large-scale population movements obscure past patterns of population structure. Here we revisit questions on the nature and time of the radiation of early humans in Africa, the extent of gene-flow among human populations as well as introgression from archaic and extinct lineages on the continent.

Population inference based on mitochondrial DNA control region data by the nearest neighbors algorithm

Population and geographic assignment are frequently undertaken using DNA sequences on the mitochondrial genome. Assignment to broad continental populations is common, although finer resolution to subpopulations can be less accurate due to shared genetic ancestry at a local level and members of different ancestral subpopulations cohabiting the same geographic area. This study reports on the accuracy of population and subpopulation assignment by using the sequence data obtained from the 3070 mitochondrial genomes and applying the K-nearest neighbors (KNN) algorithm. These data also included training samples used for continental and population assignment comprised of 1105 Europeans (including Austria, France, Germany, Spain, and England and Caucasian countries), 374 Africans (including North and East Africa and non-specific area (Pan-Africa)), and 1591 Asians (including Japan, Philippines, and Taiwan). Subpopulations included in this study were 1153 mitochondrial DNA (mtDNA) control region sequences from 12 subpopulations in Taiwan (including Han, Hakka, Ami, Atayal, Bunun, Paiwan, Puyuma, Rukai, Saisiyat, Tsou, Tao, and Pingpu). Additionally, control region sequence data from a further 50 samples, obtained from the Sigma Company, were included after they were amplified and sequenced. These additional 50 samples acted as the "testing samples" to verify the accuracy of the population. In this study, based on genetic distances as genetic metric, we used the KNN algorithm and the K-weighted-nearest neighbors (KWNN) algorithm weighted by genetic distance to classify individuals into continental populations, and subpopulations within the same continent. Accuracy results of ethnic inferences at the level of continental populations and of subpopulations among KNN and KWNN algorithms were obtained. The training sample set achieved an overall accuracy of 99 to 82% for assignment to their continental populations with K values from 1 to 101. Population assignment for subpopulations with K assignments from 1 to 5 reached an accuracy of 77 to 54%. Four out of 12 Taiwanese populations returned an accuracy of assignment of over 60%, Ami (66%), Atayal (67%), Saisiyat (66%), and Tao (80%). For the testing sample set, results of ethnic prediction for continental populations with recommended K values as 5, 10, and 35, based on results of the training sample set, achieved overall an accuracy of 100 to 94%. This study provided an accurate method in population assignment for not only continental populations but also subpopulations, which can be useful in forensic and anthropological studies.

Later Stone Age human hair from Vaalkrans Shelter, Cape Floristic Region of South Africa, reveals genetic affinity to Khoe groups

Previous studies show that the indigenous people of the southern Cape of South Africa were dramatically impacted by the arrival of European colonists starting ~400 years ago and their descendants are today mixed with Europeans and Asians. To gain insight on the occupants of the Vaalkrans Shelter located at the southernmost tip of Africa, we investigated the genetic make-up of an individual who lived there about 200 years ago. We further contextualize the genetic ancestry of this individual among prehistoric and current groups. From a hair sample excavated at the shelter, which was indirectly dated to about 200 years old, we sequenced the genome (1.01 times coverage) of a Later Stone Age individual. We analyzed the Vaalkrans genome together with genetic data from 10 ancient (pre-colonial) individuals from southern Africa spanning the last 2000 years. We show that the individual from Vaalkrans was a man who traced ~80% of his ancestry to local southern San hunter-gatherers and ~20% to a mixed East African-Eurasian source. This genetic make-up is similar to modern-day Khoekhoe individuals from the Northern Cape Province (South Africa) and Namibia, but in the southern Cape, the Vaalkrans man's descendants have likely been assimilated into mixed-ancestry "Coloured" groups. The Vaalkrans man's genome reveals that Khoekhoe pastoralist groups/individuals lived in the southern Cape as late as 200 years ago, without mixing with non-African colonists or Bantu-speaking farmers. Our findings are also consistent with the model of a Holocene pastoralist migration, originating in Eastern Africa, shaping the genomic landscape of historic and current southern African populations.

Origins of modern human ancestry

New finds in the palaeoanthropological and genomic records have changed our view of the origins of modern human ancestry. Here we review our current understanding of how the ancestry of modern humans around the globe can be traced into the deep past, and which ancestors it passes through during our journey back in time. We identify three key phases that are surrounded by major questions, and which will be at the frontiers of future research. The most recent phase comprises the worldwide expansion of modern humans between 40 and 60 thousand years ago (ka) and their last known contacts with archaic groups such as Neanderthals and Denisovans. The second phase is associated with a broadly construed African origin of modern human diversity between 60 and 300 ka. The oldest phase comprises the complex separation of modern human ancestors from archaic human groups from 0.3 to 1 million years ago. We argue that no specific point in time can currently be identified at which modern human ancestry was confined to a limited birthplace, and that patterns of the first appearance of anatomical or behavioural traits that are used to define Homo sapiens are consistent with a range of evolutionary histories.

A Late Pleistocene human humerus from Rusinga Island, Lake Victoria, Kenya

In 2010, a hominin right humerus fragment (KNM-RU 58330) was surface collected in a small gully at Nyamita North in the Late Pleistocene Wasiriya Beds of Rusinga Island, Kenya. A combination of stratigraphic and geochronological evidence suggests the specimen is likely between ∼49 and 36 ka in age. The associated fauna is diverse and dominated by semiarid grassland taxa. The small sample of associated Middle Stone Age artifacts includes Levallois flakes, cores, and retouched points. The 139 mm humeral fragment preserves the shaft from distal to the lesser tubercle to 14 mm below the distal end of the weakly projecting deltoid tuberosity. Key morphological features include a narrow and weakly marked pectoralis major insertion and a distinctive medial bend in the diaphysis at the deltoid insertion. This bend is unusual among recent human humeri but occurs in a few Late Pleistocene humeri. The dimensions of the distal end of the fragment predict a length of 317.9 ± 16.4 mm based on recent samples of African ancestry. A novel method of predicting humeral length from the distance between the middle of the pectoralis major and the bottom of the deltoid insertion predicts a length of 317.3 mm ± 17.6 mm. Cross-sectional geometry at the midshaft shows a relatively high percentage of cortical bone and a moderate degree of flattening of the shaft. The Nyamita humerus is anatomically modern in its morphology and adds to the small sample of hominins from the Late Pleistocene associated with Middle Stone Age artifacts known from East Africa. It may sample a population closely related to the people of the out-of-Africa migration.

Profiling of Mitochondrial DNA Heteroplasmy in a Prospective Oral Squamous Cell Carcinoma Study

While a shift in energy metabolism is essential to cancers, the knowledge about the involvement of the mitochondrial genome in tumorigenesis and progression in oral squamous cell carcinoma (OSCC) is still very limited. In this study, we evaluated 37 OSCC tumors and the corresponding benign mucosa tissue pairs by deep sequencing of the complete mitochondrial DNA (mtDNA). After extensive quality control, we identified 287 variants, 137 in tumor and 150 in benign samples exceeding the 1% threshold. Variant heteroplasmy levels were significantly increased in cancer compared to benign tissues (p = 0.0002). Furthermore, pairwise high heteroplasmy frequency difference variants (∆HF% > 20) with potential functional impact were increased in the cancer tissues (p = 0.024). Fourteen mutations were identified in the protein-coding region, out of which thirteen were detected in cancer and only one in benign tissue. After eight years of follow-up, the risk of mortality was higher for patients who harbored at least one ∆HF% > 20 variant in mtDNA protein-coding regions relative to those with no mutations (HR = 4.6, (95%CI = 1.3-17); p = 0.019 in primary tumor carriers). Haplogroup affiliation showed an impact on survival time, which however needs confirmation in a larger study. In conclusion, we observed a significantly higher accumulation of somatic mutations in the cancer tissues associated with a worse prognosis.

The emergence of emotionally modern humans: implications for language and learning

According to the Cooperative Breeding Hypothesis, apes with the life-history attributes of those in the line leading to the genus Homo could not have evolved unless male and female allomothers had begun to help mothers care for and provision offspring. As proposed elsewhere, the unusual way hominins reared their young generated novel phenotypes subsequently subjected to Darwinian social selection favouring those young apes best at monitoring the intentions, mental states and preferences of others and most motivated to attract and appeal to caretakers. Not only were youngsters acquiring information in social contexts different from those of other apes, but they would also have been emotionally and neurophysiologically different from them in ways that are relevant to how humans learn. Contingently delivered rewards to dependents who attracted and ingratiated themselves with allomothers shaped their behaviours and vocalizations and transformed the way developing youngsters learned from others and internalized their preferences. This article is part of the theme issue 'Life history and learning: how childhood, caregiving and old age shape cognition and culture in humans and other animals'.

Genetic Affinities among Southern Africa Hunter-Gatherers and the Impact of Admixing Farmer and Herder Populations

Southern African indigenous groups, traditionally hunter-gatherers (San) and herders (Khoekhoe), are commonly referred to as “Khoe-San” populations and have a long history in southern Africa. Their ancestors were largely isolated up until ∼2,000 years ago before the arrival of pastoralists and farmers in southern Africa. Assessing relationships among regional Khoe-San groups has been challenging due to admixture with immigrant populations that obscure past population affinities and gene flow among these autochthonous communities. We re-evaluate a combined genome-wide data set of previously published southern Africa Khoe-San populations in conjunction with novel data from Khoe-San individuals collected in Xade (Central Kalahari Game Reserve, Botswana) prior to their resettlement outside the reserve. After excluding regions in the genome that trace their ancestry to recent migrant groups, the genetic diversity of 20 Khoe-San groups fitted an isolation-by-distance model. Even though isolation-by-distance explained most genetic affinities between the different autochthonous groups, additional signals of contact between Khoe-San groups could be detected. For instance, we found stronger genetic affinities, than what would be explained by isolation-by-distance gene flow, between the two geographically separated Khoe-San groups, who speak branches of the Kx’a-language family (ǂHoan and Ju). We also scanned the genome-wide data for signals of adaptive gene flow from farmers/herders into Khoe-San groups and identified a number of genomic regions potentially introduced by the arrival of the new groups. This study provides a comprehensive picture of affinities among Khoe-San groups, prior to the arrival of recent migrants, and found that these affinities are primarily determined by the geographic landscape.

Mitochondrial DNA variation in Sub-Saharan Africa: Forensic data from a mixed West African sample, Côte d'Ivoire (Ivory Coast), and Rwanda

This study provides 398 novel complete mitochondrial control region sequences that augment the still underrepresented data from Africa by three datasets: a mixed West African sample set deriving from 12 countries (n = 145) and datasets from Côte d'Ivoire (Ivory Coast) (n = 100) as well as Rwanda (n = 153). The analysis of mtDNA variation and genetic comparisons with published data revealed low random match probabilities in all three datasets and typical West African and East African diversity, respectively. Genetic parameters indicate that the presented mixed West African dataset may serve as first forensic mtDNA control region database for West Africa in general. In addition, a strategy for responsible forensic application of precious mtDNA population samples potentially containing close maternal relatives is outlined. The datasets will be uploaded to the forensic mtDNA database EMPOP (https://empop.online) upon publication.

Whole mitochondrial genome analysis of highland Tibetan ethnicity using massively parallel sequencing

Mitochondrial DNA (mtDNA) is a key player in numerous multifaceted and intricate biological processes and plays a pivotal role in dissecting the peopling of different populations, due to its maternally inherited property and comparatively high mutation rate. In this study, 119 Tibetan individuals from the Muli Tibetan Autonomous County of China (average altitude above 3,000 m) were employed in mitochondrial genome (mitogenome) sequencing by massively parallel sequencing (MPS) techniques using the Precision ID mtDNA Whole Genome Panel on an Ion S5XL system. The dataset presented 88 distinct haplotypes, resulting in the haplotype diversity of 0.9909. The majority of haplotypes were assigned to East Asian lineages and the distribution of haplogroups of Muli Tibetan significantly differed from reference Tibetan populations. The maximum parsimony phylogeny reconstructed by 119 newly generated mitogenomes revealed 12 major Muli Tibetan lineages. Intriguingly, a Sherpa-specific sub-haplogroup A15c1 with the lack of mutations at 4216 and 15,924 was discerned in our dataset, which suggested that the maternal gene pool of Sherpas may derive from Tibetan populations. The shared haplogroups between Muli Tibetan and lowland Han Chinese hinted that these lineages may derive from non-Tibetans and have already differentiated before their arrival on the Tibetan Plateau. Furthermore, extensive pairwise population comparisons displayed that Muli Tibetan had a closer genetic relationship with ethnically or linguistically close Nyingtri Tibetan, Nyingtri Lhoba and Chamdo Tibetan populations. Genetic affinity was also observed between the Muli Tibetan and North Han Chinese. Collectively, the results generated in this study enriched the existing forensic mtDNA database and raised additional interest in the application of whole mitogenome sequencing in forensic investigations.

Human origins in a southern African palaeo-wetland and first migrations

Anatomically modern humans originated in Africa around 200 thousand years ago (ka)^1-4. Although some of the oldest skeletal remains suggest an eastern African origin², southern Africa is home to contemporary populations that represent the earliest branch of human genetic phylogeny^5,6. Here we generate, to our knowledge, the largest resource for the poorly represented and deepest-rooting maternal L0 mitochondrial DNA branch (198 new mitogenomes for a total of 1,217 mitogenomes) from contemporary southern Africans and show the geographical isolation of L0d1'2, L0k and L0g KhoeSan descendants south of the Zambezi river in Africa. By establishing mitogenomic timelines, frequencies and dispersals, we show that the L0 lineage emerged within the residual Makgadikgadi-Okavango palaeo-wetland of southern Africa⁷, approximately 200 ka (95% confidence interval, 240-165 ka). Genetic divergence points to a sustained 70,000-year-long existence of the L0 lineage before an out-of-homeland northeast-southwest dispersal between 130 and 110 ka. Palaeo-climate proxy and model data suggest that increased humidity opened green corridors, first to the northeast then to the southwest. Subsequent drying of the homeland corresponds to a sustained effective population size (L0k), whereas wet-dry cycles and probable adaptation to marine foraging allowed the southwestern migrants to achieve population growth (L0d1'2), as supported by extensive south-coastal archaeological evidence^8-10. Taken together, we propose a southern African origin of anatomically modern humans with sustained homeland occupation before the first migrations of people that appear to have been driven by regional climate changes.

Ancient DNA and bioarchaeological perspectives on European and African diversity and relationships on the colonial Delaware frontier

OBJECTIVES

Ancient DNA (aDNA) and standard osteological analyses applied to 11 skeletons at a late 17th to early 18th century farmstead site in Delaware to investigate the biological and social factors of settlement and slavery in colonial America.

MATERIALS AND METHODS

Osteological analysis and mitochondrial DNA (mtDNA) sequencing were conducted for all individuals and the resulting data contextualized with archaeological and documentary evidence.

RESULTS

Individuals of European and African descent were spatially separated in this colonial cemetery. The skeletal remains exhibited differences in osteological features and maternal genetic ancestry. A specific mtDNA haplotype appeared in a subset of the European-descended individuals suggesting they were maternally related. Individuals of African descent were not maternally related, and instead showed a diversity of haplotypes affiliated with present-day Western, Central, and Eastern regions of Africa.

DISCUSSION

Along with the bioarchaeological and documentary evidence, the aDNA findings contribute to our understanding of life on the colonial Delaware frontier. Evidence of maternal relatedness among European-descended individuals at the site demonstrates kin-based settlements in 17th century Delaware and provides preliminary identifications of individuals. The maternal genetic diversity of the individuals with African descent aligns with the routes of the trans-Atlantic slave trade but broadens our understanding of the ancestries of persons involved in it. Burial positioning, osteological pathology, and lack of maternal kinship among individuals of African descent provide tangible evidence for the emergence of racialized labor and society in Delaware during the late 17th century.

Population Structure and Implications on the Genetic Architecture of HIV-1 Phenotypes Within Southern Africa

The interesting history of Southern Africa has put the region in the spotlight for population medical genetics. Major events including the Bantu expansion and European colonialism have imprinted unique genetic signatures within autochthonous populations of Southern Africa, this resulting in differential allele frequencies across the region. This genetic structure has potential implications on susceptibility and resistance to infectious diseases such as human immunodeficiency virus (HIV) infection. Southern Africa is the region affected worst by HIV. Here, we discuss advances made in genome-wide association studies (GWAS) of HIV-1 in the past 12 years and dissect population diversity within Southern Africa. Our findings accentuate that a plethora of factors such as migration, language and culture, admixture, and natural selection have profiled the genetics of the people of Southern Africa. Genetic structure has been observed among the Khoe-San, among Bantu speakers, and between the Khoe-San, Coloureds, and Bantu speakers. Moreover, Southern African populations have complex admixture scenarios. Few GWAS of HIV-1 have been conducted in Southern Africa, with only one of these identifying two novel variants (HCG22rs2535307 and CCNG1kgp22385164) significantly associated with HIV-1 acquisition and progression. High genetic diversity, multi-wave genetic mixture and low linkage disequilibrium of Southern African populations constitute a challenge in identifying genetic variants with modest risk or protective effect against HIV-1. We therefore posit that it is compelling to assess genome-wide contribution of ancestry to HIV-1 infection. We further suggest robust methods that can pin-point population-specific variants that may contribute to the control of HIV-1 in Southern Africa.

MitoIMP: A Computational Framework for Imputation of Missing Data in Low-Coverage Human Mitochondrial Genome

The incompleteness of partial human mitochondrial genome sequences makes it difficult to perform relevant comparisons among multiple resources. To deal with this issue, we propose a computational framework for deducing missing nucleotides in the human mitochondrial genome. We applied it to worldwide mitochondrial haplogroup lineages and assessed its performance. Our approach can deduce the missing nucleotides with a precision of 0.99 or higher in most human mitochondrial DNA lineages. Furthermore, although low-coverage mitochondrial genome sequences often lead to a blurred relationship in the multidimensional scaling analysis, our approach can correct this positional arrangement according to the corresponding mitochondrial DNA lineages. Therefore, our framework will provide a practical solution to compensate for the lack of genome coverage in partial and fragmented human mitochondrial genome sequences. In this study, we developed an open-source computer program, MitoIMP, implementing our imputation procedure. MitoIMP is freely available from https://github.com/omics-tools/mitoimp.

Mitochondrial DNA in human identification: a review

Mitochondrial DNA (mtDNA) presents several characteristics useful for forensic studies, especially related to the lack of recombination, to a high copy number, and to matrilineal inheritance. mtDNA typing based on sequences of the control region or full genomic sequences analysis is used to analyze a variety of forensic samples such as old bones, teeth and hair, as well as other biological samples where the DNA content is low. Evaluation and reporting of the results requires careful consideration of biological issues as well as other issues such as nomenclature and reference population databases. In this work we review mitochondrial DNA profiling methods used for human identification and present their use in the main cases of humanidentification focusing on the most relevant issues for forensics.

No Country for Oldowan Men: Emerging Factors in Language Evolution

Language evolution has long been researched. I will review a number of broad, emerging research directions which arguably have the potential to contribute to our understanding of language evolution. Emerging topics in genomics and neurolinguistics are explored, and human-specific levels of braincase globularity - and the broader process of self-domestication within which globularity seems capable of being encapsulated - will be argued to be the central pillars of any satisfactory and interdisciplinary model of language evolution.

Whole-genome sequence analysis of a Pan African set of samples reveals archaic gene flow from an extinct basal population of modern humans into sub-Saharan populations

Background

Population demography and gene flow among African groups, as well as the putative archaic introgression of ancient hominins, have been poorly explored at the genome level.

Results

Here, we examine 15 African populations covering all major continental linguistic groups, ecosystems, and lifestyles within Africa through analysis of whole-genome sequence data of 21 individuals sequenced at deep coverage. We observe a remarkable correlation among genetic diversity and geographic distance, with the hunter-gatherer groups being more genetically differentiated and having larger effective population sizes throughout most modern-human history. Admixture signals are found between neighbor populations from both hunter-gatherer and agriculturalists groups, whereas North African individuals are closely related to Eurasian populations. Regarding archaic gene flow, we test six complex demographic models that consider recent admixture as well as archaic introgression. We identify the fingerprint of an archaic introgression event in the sub-Saharan populations included in the models (~ 4.0% in Khoisan, ~ 4.3% in Mbuti Pygmies, and ~ 5.8% in Mandenka) from an early divergent and currently extinct ghost modern human lineage.

Conclusion

The present study represents an in-depth genomic analysis of a Pan African set of individuals, which emphasizes their complex relationships and demographic history at population level.

Electronic supplementary material

The online version of this article (10.1186/s13059-019-1684-5) contains supplementary material, which is available to authorized users.

African evolutionary history inferred from whole genome sequence data of 44 indigenous African populations

BACKGROUND

Africa is the origin of modern humans within the past 300 thousand years. To infer the complex demographic history of African populations and adaptation to diverse environments, we sequenced the genomes of 92 individuals from 44 indigenous African populations.

RESULTS

Genetic structure analyses indicate that among Africans, genetic ancestry is largely partitioned by geography and language, though we observe mixed ancestry in many individuals, consistent with both short- and long-range migration events followed by admixture. Phylogenetic analysis indicates that the San genetic lineage is basal to all modern human lineages. The San and Niger-Congo, Afroasiatic, and Nilo-Saharan lineages were substantially diverged by 160 kya (thousand years ago). In contrast, the San and Central African rainforest hunter-gatherer (CRHG), Hadza hunter-gatherer, and Sandawe hunter-gatherer lineages were diverged by ~ 120-100 kya. Niger-Congo, Nilo-Saharan, and Afroasiatic lineages diverged more recently by ~ 54-16 kya. Eastern and western CRHG lineages diverged by ~ 50-31 kya, and the western CRHG lineages diverged by ~ 18-12 kya. The San and CRHG populations maintained the largest effective population size compared to other populations prior to 60 kya. Further, we observed signatures of positive selection at genes involved in muscle development, bone synthesis, reproduction, immune function, energy metabolism, and cell signaling, which may contribute to local adaptation of African populations.

CONCLUSIONS

We observe high levels of genomic variation between ethnically diverse Africans which is largely correlated with geography and language. Our study indicates ancient population substructure and local adaptation of Africans.

Tales of Human Migration, Admixture, and Selection in Africa

In the last three decades, genetic studies have played an increasingly important role in exploring human history. They have helped to conclusively establish that anatomically modern humans first appeared in Africa roughly 250,000–350,000 years before present and subsequently migrated to other parts of the world. The history of humans in Africa is complex and includes demographic events that influenced patterns of genetic variation across the continent. Through genetic studies, it has become evident that deep African population history is captured by relationships among African hunter–gatherers, as the world's deepest population divergences occur among these groups, and that the deepest population divergence dates to 300,000 years before present. However, the spread of pastoralism and agriculture in the last few thousand years has shaped the geographic distribution of present-day Africans and their genetic diversity. With today's sequencing technologies, we can obtain full genome sequences from diverse sets of extant and prehistoric Africans. The coming years will contribute exciting new insights toward deciphering human evolutionary history in Africa.

Carriers of mitochondrial DNA macrohaplogroup L3 basal lineages migrated back to Africa from Asia around 70,000 years ago

Background

The main unequivocal conclusion after three decades of phylogeographic mtDNA studies is the African origin of all extant modern humans. In addition, a southern coastal route has been argued for to explain the Eurasian colonization of these African pioneers. Based on the age of macrohaplogroup L3, from which all maternal Eurasian and the majority of African lineages originated, the out-of-Africa event has been dated around 60-70 kya. On the opposite side, we have proposed a northern route through Central Asia across the Levant for that expansion and, consistent with the fossil record, we have dated it around 125 kya. To help bridge differences between the molecular and fossil record ages, in this article we assess the possibility that mtDNA macrohaplogroup L3 matured in Eurasia and returned to Africa as basal L3 lineages around 70 kya.

Results

The coalescence ages of all Eurasian (M,N) and African (L3 ) lineages, both around 71 kya, are not significantly different. The oldest M and N Eurasian clades are found in southeastern Asia instead near of Africa as expected by the southern route hypothesis. The split of the Y-chromosome composite DE haplogroup is very similar to the age of mtDNA L3. An Eurasian origin and back migration to Africa has been proposed for the African Y-chromosome haplogroup E. Inside Africa, frequency distributions of maternal L3 and paternal E lineages are positively correlated. This correlation is not fully explained by geographic or ethnic affinities. This correlation rather seems to be the result of a joint and global replacement of the old autochthonous male and female African lineages by the new Eurasian incomers.

Conclusions

These results are congruent with a model proposing an out-of-Africa migration into Asia, following a northern route, of early anatomically modern humans carrying pre-L3 mtDNA lineages around 125 kya, subsequent diversification of pre-L3 into the basal lineages of L3, a return to Africa of Eurasian fully modern humans around 70 kya carrying the basal L3 lineages and the subsequent diversification of Eurasian-remaining L3 lineages into the M and N lineages in the outside-of-Africa context, and a second Eurasian global expansion by 60 kya, most probably, out of southeast Asia. Climatic conditions and the presence of Neanderthals and other hominins might have played significant roles in these human movements. Moreover, recent studies based on ancient DNA and whole-genome sequencing are also compatible with this hypothesis.

Electronic supplementary material

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

Genomes of all known members of a Plasmodium subgenus reveal paths to virulent human malaria

Plasmodium falciparum, the most virulent agent of human malaria, shares a recent common ancestor with the gorilla parasite Plasmodium praefalciparum. Little is known about the other gorilla- and chimpanzee-infecting species in the same (Laverania) subgenus as P. falciparum, but none of them are capable of establishing repeated infection and transmission in humans. To elucidate underlying mechanisms and the evolutionary history of this subgenus, we have generated multiple genomes from all known Laverania species. The completeness of our dataset allows us to conclude that interspecific gene transfers, as well as convergent evolution, were important in the evolution of these species. Striking copy number and structural variations were observed within gene families and one, stevor, shows a host-specific sequence pattern. The complete genome sequence of the closest ancestor of P. falciparum enables us to estimate the timing of the beginning of speciation to be 40,000-60,000 years ago followed by a population bottleneck around 4,000-6,000 years ago. Our data allow us also to search in detail for the features of P. falciparum that made it the only member of the Laverania able to infect and spread in humans.

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.

Human phylogeography and diversity

Homo sapiens phylogeography begins with the species' origin nearly 200 kya in Africa. First signs of the species outside Africa (in Arabia) are from 125 kya. Earliest dates elsewhere are now 100 kya in China, 45 kya in Australia and southern Europe (maybe even 60 kya in Australia), 32 kya in northeast Siberia, and maybe 20 kya in the Americas. Humans reached arctic regions and oceanic islands last-arctic North America about 5 kya, mid- and eastern Pacific islands about 2-1 kya, and New Zealand about 700 y ago. Initial routes along coasts seem the most likely given abundant and easily harvested shellfish there as indicated by huge ancient oyster shell middens on all continents. Nevertheless, the effect of geographic barriers-mountains and oceans-is clear. The phylogeographic pattern of diasporas from several single origins-northeast Africa to Eurasia, southeast Eurasia to Australia, and northeast Siberia to the Americas-allows the equivalent of a repeat experiment on the relation between geography and phylogenetic and cultural diversity. On all continents, cultural diversity is high in productive low latitudes, presumably because such regions can support populations of sustainable size in a small area, therefore allowing a high density of cultures. Of course, other factors operate. South America has an unusually low density of cultures in its tropical latitudes. A likely factor is the phylogeographic movement of peoples from the Old World bringing novel and hence, lethal diseases to the New World, a foretaste, perhaps, of present day global transport of tropical diseases.

Southern African ancient genomes estimate modern human divergence to 350,000 to 260,000 years ago

Southern Africa is consistently placed as a potential region for the evolution of Homo sapiens We present genome sequences, up to 13x coverage, from seven ancient individuals from KwaZulu-Natal, South Africa. The remains of three Stone Age hunter-gatherers (about 2000 years old) were genetically similar to current-day southern San groups, and those of four Iron Age farmers (300 to 500 years old) were genetically similar to present-day Bantu-language speakers. We estimate that all modern-day Khoe-San groups have been influenced by 9 to 30% genetic admixture from East Africans/Eurasians. Using traditional and new approaches, we estimate the first modern human population divergence time to between 350,000 and 260,000 years ago. This estimate increases the deepest divergence among modern humans, coinciding with anatomical developments of archaic humans into modern humans, as represented in the local fossil record.

MitoSuite: a graphical tool for human mitochondrial genome profiling in massive parallel sequencing

Recent rapid advances in high-throughput, next-generation sequencing (NGS) technologies have promoted mitochondrial genome studies in the fields of human evolution, medical genetics, and forensic casework. However, scientists unfamiliar with computer programming often find it difficult to handle the massive volumes of data that are generated by NGS. To address this limitation, we developed MitoSuite, a user-friendly graphical tool for analysis of data from high-throughput sequencing of the human mitochondrial genome. MitoSuite generates a visual report on NGS data with simple mouse operations. Moreover, it analyzes high-coverage sequencing data but runs on a stand-alone computer, without the need for file upload. Therefore, MitoSuite offers outstanding usability for handling massive NGS data, and is ideal for evolutionary, clinical, and forensic studies on the human mitochondrial genome variations. It is freely available for download from the website https://mitosuite.com.

Origin and spread of human mitochondrial DNA haplogroup U7

Human mitochondrial DNA haplogroup U is among the initial maternal founders in Southwest Asia and Europe and one that best indicates matrilineal genetic continuity between late Pleistocene hunter-gatherer groups and present-day populations of Europe. While most haplogroup U subclades are older than 30 thousand years, the comparatively recent coalescence time of the extant variation of haplogroup U7 (~16–19 thousand years ago) suggests that its current distribution is the consequence of more recent dispersal events, despite its wide geographical range across Europe, the Near East and South Asia. Here we report 267 new U7 mitogenomes that – analysed alongside 100 published ones – enable us to discern at least two distinct temporal phases of dispersal, both of which most likely emanated from the Near East. The earlier one began prior to the Holocene (~11.5 thousand years ago) towards South Asia, while the later dispersal took place more recently towards Mediterranean Europe during the Neolithic (~8 thousand years ago). These findings imply that the carriers of haplogroup U7 spread to South Asia and Europe before the suggested Bronze Age expansion of Indo-European languages from the Pontic-Caspian Steppe region.

The Impact of Genetic and Non-Genetic Factors on Warfarin Dose Prediction in MENA Region: A Systematic Review

Background

Warfarin is the most commonly used oral anticoagulant for the treatment and prevention of thromboembolic disorders. Pharmacogenomics studies have shown that variants in CYP2C9 and VKORC1 genes are strongly and consistently associated with warfarin dose variability. Although different populations from the Middle East and North Africa (MENA) region may share the same ancestry, it is still unclear how they compare in the genetic and non-genetic factors affecting their warfarin dosing.

Objective

To explore the prevalence of CYP2C9 and VKORC1 variants in MENA, and the effect of these variants along with other non-genetic factors in predicting warfarin dose.

Methods

In this systematic review, we included observational cross sectional and cohort studies that enrolled patients on stable warfarin dose and had the genetics and non-genetics factors associated with mean warfarin dose as the primary outcome. We searched PubMed, Medline, Scopus, PharmGKB, PHGKB, Google scholar and reference lists of relevant reviews.

Results

We identified 17 studies in eight different populations: Iranian, Israeli, Egyptian, Lebanese, Omani, Kuwaiti, Sudanese and Turkish. Most common genetic variant in all populations was the VKORC1 (-1639G>A), with a minor allele frequency ranging from 30% in Egyptians and up to 52% and 56% in Lebanese and Iranian, respectively. Variants in the CYP2C9 were less common, with the highest MAF for CYP2C9*2 among Iranians (27%). Variants in the VKORC1 and CYP2C9 were the most significant predictors of warfarin dose in all populations. Along with other genetic and non-genetic factors, they explained up to 63% of the dose variability in Omani and Israeli patients.

Conclusion

Variants of VKORC1 and CYP2C9 are the strongest predictors of warfarin dose variability among the different populations from MENA. Although many of those populations share the same ancestry and are similar in their warfarin dose predictors, a population specific dosing algorithm is needed for the prospective estimation of warfarin dose.

Worldwide genetic and cultural change in human evolution

Both genetic variation and certain culturally transmitted phenotypes show geographic signatures of human demographic history. As a result of the human cultural predisposition to migrate to new areas, humans have adapted to a large number of different environments. Migration to new environments alters genetic selection pressures, and comparative genetic studies have pinpointed numerous likely targets of this selection. However, humans also exhibit many cultural adaptations to new environments, such as practices related to clothing, shelter, and food. Human culture interacts with genes and the environment in complex ways, and studying genes and culture together can deepen our understanding of human evolution.

mtDNA diversity in human populations highlights the merit of haplotype matching in gene therapies

STUDY QUESTION

Does mitochondrial DNA (mtDNA) diversity in modern human populations potentially pose a challenge, via mtDNA segregation, to mitochondrial replacement therapies?

SUMMARY ANSWER

The magnitude of mtDNA diversity in modern human populations is as high as in mammalian model systems where strong mtDNA segregation is observed; consideration of haplotype pairs and/or haplotype matching can help avoid these potentially deleterious effects.

WHAT IS KNOWN ALREADY

In mammalian models, substantial proliferative differences are observed between different mtDNA haplotypes in cellular admixtures, with larger proliferative differences arising from more diverse haplotype pairings. If maternal mtDNA is 'carried over' in human gene therapies, these proliferative differences could lead to its amplification in the resulting offspring, potentially leading to manifestation of the disease that the therapy was designed to avoid-but existing studies have not investigated whether mtDNA diversity in modern human populations is sufficient to permit significant amplification.

STUDY DESIGN, SIZE, DURATION

This theoretical study used over 7500 human mtDNA sequences from The National Center for Biotechnology Information (NCBI), a range of international and British mtDNA surveys, and 2011 census data.

PARTICIPANTS/MATERIALS, SETTING, METHODS

A stochastic simulation approach was used to model random haplotype pairings from within different regions. In total, 1000 simulated pairings were analysed using the basic local alignment search tool (BLAST) for each region. Previous data from mouse models were used to estimate proliferative differences.

MAIN RESULTS AND THE ROLE OF CHANCE

Even within the same haplogroup, differences of around 20-80 single-nucleotide polymorphisms (SNPs) are common between mtDNAs admixed in random pairings. These values are sufficient to lead to substantial segregation in mouse models over an organismal lifetime, even given low starting heteroplasmy, inducing increases from 5% to 35% over 1 year. Substantial population mixing in modern UK cities increases the expected genetic differences. Hence, the likely genetic differences between humans randomly sampled from a population may well allow substantial amplification of a disease-carrying mtDNA haplotype over the timescale of a human lifetime. We report ranges and mean differences for all statistics to quantify uncertainty in our results.

LIMITATIONS/REASONS FOR CAUTION

The mapping from mouse and other mammalian models to the human system is challenging, as timescales and mechanisms may differ. Reporting biases in NCBI mtDNA data, if present, may affect the statistics we compute. We discuss the robustness of our findings in the light of these concerns.

WIDER IMPLICATIONS OF THE FINDINGS

Matching the mtDNA haplotypes of the mother and third-party donor in mitochondrial replacement therapies is supported as a means of ameliorating the potentially deleterious results of human mtDNA diversity. We present a chart of expected SNP differences between mtDNA haplogroups, allowing the selection of optimal partners for therapies.

LARGE SCALE DATA

N/A STUDY FUNDING/COMPETING INTERESTS: The authors report no external funding sources or conflicts of interest.