The role of selection in the evolution of human mitochondrial genomes

The rate and nature of mitochondrial DNA mutations in human pedigrees

We examined the rate and nature of mitochondrial DNA (mtDNA) mutations in humans using sequence data from 64,806 contemporary Icelanders from 2,548 matrilines. Based on 116,663 mother-child transmissions, 8,199 mutations were detected, providing robust rate estimates by nucleotide type, functional impact, position, and different alleles at the same position. We thoroughly document the true extent of hypermutability in mtDNA, mainly affecting the control region but also some coding-region variants. The results reveal the impact of negative selection on viable deleterious mutations, including rapidly mutating disease-associated 3243A>G and 1555A>G and pre-natal selection that most likely occurs during the development of oocytes. Finally, we show that the fate of new mutations is determined by a drastic germline bottleneck, amounting to an average of 3 mtDNA units effectively transmitted from mother to child.

Mitochondria in biology and medicine - 2023

Mitochondria are an indispensable part of the cell that plays a crucial role in regulating various signaling pathways, energy metabolism, cell differentiation, proliferation, and cell death. Since mitochondria have their own genetic material, they differ from their nuclear counterparts, and dysregulation is responsible for a broad spectrum of diseases. Mitochondrial dysfunction is associated with several disorders, including neuro-muscular disorders, cancer, and premature aging, among others. The intricacy of the field is due to the cross-talk between nuclear and mitochondrial genes, which has also improved our knowledge of mitochondrial functions and their pathogenesis. Therefore, interdisciplinary research and communication are crucial for mitochondrial biology and medicine due to the challenges they pose for diagnosis and treatment. The ninth annual conference of the Society for Mitochondria Research and Medicine (SMRM)- India, titled "Mitochondria in Biology and Medicine" was organized at the Centre for DNA Fingerprinting and Diagnostics (CDFD), Hyderabad, India, on June 21-23, 2023. The latest advancements in the field of mitochondrial biology and medicine were discussed at the conference. In this article, we summarize the entire event for the benefit of researchers working in the field of mitochondrial biology and medicine.

Spatiotemporal dynamics in butterfly hybrid zones

Evaluating whether hybrid zones are stable or mobile can provide novel insights for evolution and conservation biology. Butterflies exhibit high sensitivity to environmental changes and represent an important model system for the study of hybrid zone origins and maintenance. Here, we review the literature exploring butterfly hybrid zones, with a special focus on their spatiotemporal dynamics and the potential mechanisms that could lead to their movement or stability. We then compare different lines of evidence used to investigate hybrid zone dynamics and discuss the strengths and weaknesses of each approach. Our goal with this review is to reveal general conditions associated with the stability or mobility of butterfly hybrid zones by synthesizing evidence obtained using different types of data sampled across multiple regions and spatial scales. Finally, we discuss spatiotemporal dynamics in the context of a speciation/divergence continuum, the relevance of hybrid zones for conservation biology, and recommend key topics for future investigation.

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.

Evaluating the role of selection in the evolution of mitochondrial genomes of aboriginal peoples of Siberia

Studies of the nature of mitochondrial DNA (mtDNA) variability in human populations have shown that protein-coding genes are under negative (purifying) selection, since their mutation spectra are characterized by a pronounced predominance of synonymous substitutions over non-synonymous ones (Ka/Ks < 1). Meanwhile, a number of studies have shown that the adaptation of populations to various environmental conditions may be accompanied by a relaxation of negative selection in some mtDNA genes. For example, it was previously found that in Arctic populations, negative selection is relaxed in the mitochondrial ATP6 gene, which encodes one of the subunits of ATP synthase. In this work, we performed a Ka/Ks analysis of mitochondrial genes in large samples of three regional population groups in Eurasia: Siberia (N = 803), Western Asia/Transcaucasia (N = 753), and Eastern Europe (N = 707). The main goal of this work is to search for traces of adaptive evolution in the mtDNA genes of aboriginal peoples of Siberia represented by populations of the north (Koryaks, Evens) and the south of Siberia and the adjacent territory of Northeast China (Buryats, Barghuts, Khamnigans). Using standard Ka/Ks analysis, it was found that all mtDNA genes in all studied regional population groups are subject to negative selection. The highest Ka/Ks values in different regional samples were found in almost the same set of genes encoding subunits of ATP synthase (ATP6, ATP8), NADH dehydrogenase complex (ND1, ND2, ND3), and cytochrome bc1 complex (CYB). The highest Ka/Ks value, indicating a relaxation of negative selection, was found in the ATP6 gene in the Siberian group. The results of the analysis performed using the FUBAR method (HyPhy software package) and aimed at searching for mtDNA codons under the influence of selection also showed the predominance of negative selection over positive selection in all population groups. In Siberian populations, nucleotide sites that are under positive selection and associated with mtDNA haplogroups were registered not in the north (which is expected under the assumption of adaptive evolution of mtDNA), but in the south of Siberia.

m.4216 T > C polymorphism in JT cluster determines a lower pregnancy rate in response to controlled ovarian stimulation treatment

PURPOSE

To analyze the influence of Caucasian mitochondrial haplogroups on controlled ovarian stimulation outcome (COS), embryo (E), and pregnancy success.

METHODS

In a Caucasian population (n = 517) undergoing COS, mitochondrial haplogroups and physiological parameters were determined. Patients were classified, according to Bologna criteria, as good (>3)/poor ≤3) responder, on dependence of recruited oocytes (RO), and in pregnancy/non-pregnancy groups. Haplogroups were determined by sequencing mitochondrial hypervariable sequence I and confirmed by polymerase chain reaction (PCR), followed by restriction fragment length polymorphisms (RFLP).

RESULTS

The rank of total dose of FSH (TD FSH) was similar in all clusters/haplogroups, except in JT, which is narrower (950-3,650 IU), particularly in T (1,350-3,650 IU). The statistical analysis showed higher RO and E in JT when compared to U, although it was only Uk which accumulated significantly in pregnancy respect to JT. Pearson's correlations between TD FSH and RO showed negative statistical significance in all population (P = 0.001), H (P = 0.03), JT (P = 0.01), and T (P = 0.03). The percentage of contribution of TD FSH on RO was almost nine times in the JT cluster as compared to all population one.

CONCLUSIONS

JT cluster shows a different influence of TD FSH on RO. JT cluster shows higher RO and E than U, but it is Uk which exhibits a significant higher pregnancy rate than JT. The negative influence of the JT cluster on pregnancy success strongly suggests that the m.4216 T > C polymorphism could be responsible.

Organelle bottlenecks facilitate evolvability by traversing heteroplasmic fitness valleys

Bioenergetic organelles-mitochondria and plastids-retain their own genomes (mtDNA and ptDNA), and these organelle DNA (oDNA) molecules are vital for eukaryotic life. Like all genomes, oDNA must be able to evolve to suit new environmental challenges. However, mixed oDNA populations in cells can challenge cellular bioenergetics, providing a penalty to the appearance and adaptation of new mutations. Here we show that organelle "bottlenecks," mechanisms increasing cell-to-cell oDNA variability during development, can overcome this mixture penalty and facilitate the adaptation of beneficial mutations. We show that oDNA heteroplasmy and bottlenecks naturally emerge in evolutionary simulations subjected to fluctuating environments, demonstrating that this evolvability is itself evolvable. Usually thought of as a mechanism to clear damaging mutations, organelle bottlenecks therefore also resolve the tension between intracellular selection for pure cellular oDNA populations and the "bet-hedging" need for evolvability and adaptation to new environments. This general theory suggests a reason for the maintenance of organelle heteroplasmy in cells, and may explain some of the observed diversity in organelle maintenance and inheritance across taxa.

Tracing genetic connections of ancient Hungarians to the 6th-14th century populations of the Volga-Ural region

Most of the early Hungarian tribes originated from the Volga-Kama and South-Ural regions, where they were composed of a mixed population based on historical, philological and archaeological data. We present here the uniparental genetic makeup of the mediaeval era of these regions that served as a melting pot for ethnic groups with different linguistic and historical backgrounds. Representing diverse cultural contexts, the new genetic data originate from ancient proto-Ob-Ugric people from Western Siberia (6th-13th century), the pre-Conquest period and subsisting Hungarians from the Volga-Ural region (6th-14th century) and their neighbours. By examining the eastern archaeology traits of Hungarian prehistory, we also study their genetic composition and origin in an interdisciplinary framework. We analyzed 110 deep-sequenced mitogenomes and 42 Y-chromosome haplotypes from 18 archaeological sites in Russia. The results support the studied groups' genetic relationships regardless of geographical distances, suggesting large-scale mobility. We detected long-lasting genetic connections between the sites representing the Kushnarenkovo and Chiyalik cultures and the Carpathian Basin Hungarians and confirmed the Uralic transmission of several East Eurasian uniparental lineages in their gene pool. Based on phylogenetics, we demonstrate and model the connections and splits of the studied Volga-Ural and conqueror groups. Early Hungarians and their alliances conquered the Carpathian Basin around 890 AD. Re-analysis of the Hungarian conquerors' maternal gene pool reveals numerous surviving maternal relationships in both sexes; therefore, we conclude that men and women came to the Carpathian Basin together, and although they were subsequently genetically fused into the local population, certain eastern lineages survived for centuries.

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.

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 Role of Mitochondria in Carcinogenesis

The mitochondria are essential for normal cell functioning. Changes in mitochondrial DNA (mtDNA) may affect the occurrence of some chronic diseases and cancer. This process is complex and not entirely understood. The assignment to a particular mitochondrial haplogroup may be a factor that either contributes to cancer development or reduces its likelihood. Mutations in mtDNA occurring via an increase in reactive oxygen species may favour the occurrence of further changes both in mitochondrial and nuclear DNA. Mitochondrial DNA mutations in postmitotic cells are not inherited, but may play a role both in initiation and progression of cancer. One of the first discovered polymorphisms associated with cancer was in the gene NADH-ubiquinone oxidoreductase chain 3 (mt-ND3) and it was typical of haplogroup N. In prostate cancer, these mutations and polymorphisms involve a gene encoding subunit I of respiratory complex IV cytochrome c oxidase subunit 1 gene (COI). At present, a growing number of studies also address the impact of mtDNA polymorphisms on prognosis in cancer patients. Some of the mitochondrial DNA polymorphisms occur in both chronic disease and cancer, for instance polymorphism G5913A characteristic of prostate cancer and hypertension.

Mitogenomes Reveal Multiple Colonization of Mountains by Rattus in Sundaland

Tropical mountains are cradles of biodiversity and endemism. Sundaland, tropical Southeast Asia, hosts 3 species of Rattus endemic to elevations above 2000 m with an apparent convergence in external morphology: Rattus korinchi and R. hoogerwerfi from Sumatra, and R. baluensis from Borneo. A fourth one, R. tiomanicus, is restricted to lowland elevations across the whole region. The origins of these endemics are little known due to the absence of a robust phylogenetic framework. We use complete mitochondrial genomes from the 3 high altitude Rattus, and several related species to determine their relationships, date divergences, reconstruct their history of colonization, and test for selection on the mitochondrial DNA. We show that mountain colonization happened independently in Borneo (<390 Kya) and Sumatra (~1.38 Mya), likely from lowland lineages. The origin of the Bornean endemic R. baluensis is very recent and its genetic diversity is nested within the diversity of R. tiomanicus. We found weak evidence of positive selection in the high-elevation lineages and attributed the greater nonsynonymous mutations on these branches (specially R. baluensis) to lesser purifying selection having acted on the terminal branches in the phylogeny.

Association of whole mtDNA, an NADPH G11914A variant, and haplogroups with high physical performance in an elite military troop

Physical performance is a multifactorial and complex trait influenced by environmental and hereditary factors. Environmental factors alone have been insufficient to characterize all outstanding phenotypes. Recent advances in genomic technologies have enabled the investigation of whole nuclear and mitochondrial genome sequences, increasing our ability to understand interindividual variability in physical performance. Our objective was to evaluate the association of mitochondrial polymorphic loci with physical performance in Brazilian elite military personnel. Eighty-eight male military personnel who participated in the Command Actions Course of the Army were selected. Total DNA was obtained from blood samples and a complete mitochondrial genome (mtDNA) was sequenced using Illumina MiSeq platform. Twenty-nine subjects completed the training program (FINISHED, 'F'), and fifty-nine failed to complete (NOT_FINISHED, 'NF'). The mtDNA from NF was slightly more similar to genomes from African countries frequently related to endurance level. Twenty-two distinct mtDNA haplogroups were identified corroborating the intense genetic admixture of the Brazilian population, but their distribution was similar between the two groups (FST=0.0009). Of 745 polymorphisms detected in the mtDNA, the position G11914A within the NADPH gene component of the electron transport chain, was statistically different between F and NF groups (P=0.011; OR: 4.286; 95%CI: 1.198-16.719), with a higher frequency of the G allele in group F individuals). The high performance of military personnel may be mediated by performance-related genomic traits. Thus, mitochondrial genetic markers such as the ND4 gene may play an important role on physical performance variability.

Is Global Microbial Biodiversity Increasing, Decreasing, or Staying the Same?

Animal and plant biodiversity is decreasing. In contrast, the global direction and the pace of change in microbial, including viral, biodiversity is unknown. Important niches for microbial diversity occur in highly specific associations with plants and animals, and these niches are lost as hosts become extinct. The taxonomic diversity of human gut bacteria is reported to be decreasing. On the other hand, SARS-CoV-2 variation is increasing. Where microbes are concerned, Darwin’s “tangled bank” of interdependent organisms may be composed mostly of other microbes. There is the likelihood that as some classes of microbes become extinct, others evolve and diversify. A better handle on all processes that affect microbial biodiversity and their net balance is needed. Lack of insight into the dynamics of evolution of microbial biodiversity is arguably the single most profound and consequential unknown with regard to human knowledge of the biosphere. If some or all parts of microbial diversity are relentlessly increasing, then survey approaches may be too slow to ever catch up. New approaches, including single-molecule or single-cell sequencing in populations, as well as focused attention on modulators and vectors of vertical and horizontal evolution may offer more direct insights into some aspects of the pace of microbial evolution.

Ancient and modern mitogenomes from Central Argentina: new insights into population continuity, temporal depth and migration in South America

The inverted triangle shape of South America places Argentina territory as a geographical crossroads between the two principal peopling streams that followed either the Pacific or the Atlantic coasts, which could have then merged in Central Argentina (CA). Although the genetic diversity from this region is therefore crucial to decipher past population movements in South America, its characterization has been overlooked so far. We report 92 modern and 22 ancient mitogenomes spanning a temporal range of 5000 years, which were compared with a large set of previously reported data. Leveraging this dataset representative of the mitochondrial diversity of the subcontinent, we investigate the maternal history of CA populations within a wider geographical context. We describe a large number of novel clades within the mitochondrial DNA tree, thus providing new phylogenetic interpretations for South America. We also identify several local clades of great temporal depth with continuity until the present time, which stem directly from the founder haplotypes, suggesting that they originated in the region and expanded from there. Moreover, the presence of lineages characteristic of other South American regions reveals the existence of gene flow to CA. Finally, we report some lineages with discontinuous distribution across the Americas, which suggest the persistence of relic lineages likely linked to the first population arrivals. The present study represents to date the most exhaustive attempt to elaborate a Native American genetic map from modern and ancient complete mitochondrial genomes in Argentina and provides relevant information about the general process of settlement in South America.

Impact of modular mitochondrial epistatic interactions on the evolution of human subpopulations

Investigation of human mitochondrial (mt) genome variation has been shown to provide insights to the human history and natural selection. By analyzing 24,167 human mt-genome samples, collected for five continents, we have developed a co-mutation network model to investigate characteristic human evolutionary patterns. The analysis highlighted richer co-mutating regions of the mt-genome, suggesting the presence of epistasis. Specifically, a large portion of COX genes was found to co-mutate in Asian and American populations, whereas, in African, European, and Oceanic populations, there was greater co-mutation bias in hypervariable regions. Interestingly, this study demonstrated hierarchical modularity as a crucial agent for these co-mutation networks. More profoundly, our ancestry-based co-mutation module analyses showed that mutations cluster preferentially in known mitochondrial haplogroups. Contemporary human mt-genome nucleotides most closely resembled the ancestral state, and very few of them were found to be ancestral-variants. Overall, these results demonstrated that subpopulation-based biases may favor mitochondrial gene specific epistasis.

Extreme heterogeneity of human mitochondrial DNA from organelles to populations

Contrary to the long-held view that most humans harbour only identical mitochondrial genomes, deep resequencing has uncovered unanticipated extreme genetic variation within mitochondrial DNA (mtDNA). Most, if not all, humans contain multiple mtDNA genotypes (heteroplasmy); specific patterns of variants accumulate in different tissues, including cancers, over time; and some variants are preferentially passed down or suppressed in the maternal germ line. These findings cast light on the origin and spread of mtDNA mutations at multiple scales, from the organelle to the human population, and challenge the conventional view that high percentages of a mutation are required before a new variant has functional consequences.

Genetic hitchhiking, mitonuclear coadaptation, and the origins of mt DNA barcode gaps

DNA barcoding based on mitochondrial (mt) nucleotide sequences is an enigma. Neutral models of mt evolution predict DNA barcoding cannot work for recently diverged taxa, and yet, mt DNA barcoding accurately delimits species for many bilaterian animals. Meanwhile, mt DNA barcoding often fails for plants and fungi. I propose that because mt gene products must cofunction with nuclear gene products, the evolution of mt genomes is best understood with full consideration of the two environments that impose selective pressure on mt genes: the external environment and the internal genomic environment. Moreover, it is critical to fully consider the potential for adaptive evolution of not just protein products of mt genes but also of mt transfer RNAs and mt ribosomal RNAs. The tight linkage of genes on mt genomes that do not engage in recombination could facilitate selective sweeps whenever there is positive selection on any element in the mt genome, leading to the purging of mt genetic diversity within a population and to the rapid fixation of novel mt DNA sequences. Accordingly, the most important factor determining whether or not mt DNA sequences diagnose species boundaries may be the extent to which the mt chromosomes engage in recombination.

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.

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.

A worldwide view of matriliny: using cross-cultural analyses to shed light on human kinship systems

Although matriliny and matrilocality are relatively rare in contemporary human populations, these female-based descent and residence systems are present in different cultural contexts and across the globe. Previous research has generated numerous hypotheses about which cultural traits are associated with the stability or loss of matrilineal descent. In addition, several studies have examined matrilineal descent with phylogenetic analyses; however, the use of language phylogenies has restricted these analyses to comparisons within a single language family, often confined to a single continent. Cross-cultural comparisons are particularly informative when they account for the relationships between widely distributed populations, as opposed to treating each population as an independent sample or focusing on a single region. Here, we study the evolution of descent systems on a worldwide scale. First, we test for significant associations between matriliny and numerous cultural traits that have been theoretically associated with its stability or loss, such as subsistence strategy, animal domestication, mating system, residence pattern, wealth transfer and property succession. In addition, by combining genetic and linguistic information to build a global supertree that includes 16 matrilineal populations, we also perform phylogenetically controlled analyses to assess the patterns of correlated evolution between descent and other traits: for example, does a change in subsistence strategy generally predict a shift in the rules of descent, or do these transitions happen independently? These analyses enable a worldwide perspective on the pattern and process of the evolution of matriliny and matrilocality. This article is part of the theme issue 'The evolution of female-biased kinship in humans and other mammals'.

Discovery of mitochondrial DNA variants associated with genome-wide blood cell gene expression: a population-based mtDNA sequencing study

The effect of mitochondrial DNA (mtDNA) variation on peripheral blood transcriptomics in health and disease is not fully known. Sex-specific mitochondrially controlled gene expression patterns have been shown in Drosophila melanogaster but in humans, evidence is lacking. Functional variation in mtDNA may also have a role in the development of type 2 diabetes and its precursor state, i.e. prediabetes. We examined the associations between mitochondrial single-nucleotide polymorphisms (mtSNPs) and peripheral blood transcriptomics with a focus on sex- and prediabetes-specific effects. The genome-wide blood cell expression data of 19 637 probes, 199 deep-sequenced mtSNPs and nine haplogroups of 955 individuals from a population-based Young Finns Study cohort were used. Significant associations were identified with linear regression and analysis of covariance. The effects of sex and prediabetes on the associations between gene expression and mtSNPs were studied using random-effect meta-analysis. Our analysis identified 53 significant expression probe-mtSNP associations after Bonferroni correction, involving 7 genes and 31 mtSNPs. Eight probe-mtSNP signals remained independent after conditional analysis. In addition, five genes showed differential expression between haplogroups. The meta-analysis did not show any significant differences in linear model effect sizes between males and females but identified the association between the OASL gene and mtSNP C16294T to show prediabetes-specific effects. This study pinpoints new independent mtSNPs associated with peripheral blood transcriptomics and replicates six previously reported associations, providing further evidence of the mitochondrial genetic control of blood cell gene expression. In addition, we present evidence that prediabetes might lead to perturbations in mitochondrial control.

Mutation signatures in germline mitochondrial genome provide insights into human mitochondrial evolution and disease

Variations in mitochondrial DNA (mtDNA) have been fundamental for understanding human evolution and are causative for a plethora of inherited mitochondrial diseases, but the mutation signatures of germline mtDNA and their value in understanding mitochondrial pathogenicity remain unknown. Here, we carried out a systematic analysis of mutation patterns in germline mtDNA based on 97,566 mtDNA variants from 45,494 full-length sequences and revealed a highly non-stochastic and replication-coupled mutation signature characterized by nucleotide-specific mutation pressure (G > T>A > C) and position-specific selection pressure, suggesting the existence of an intensive mutation-selection interplay in germline mtDNA. We provide evidence that this mutation-selection interplay has strongly shaped the mtDNA sequence during evolution, which not only manifests as an oriented alteration of amino acid compositions of mitochondrial encoded proteins, but also explains the long-lasting mystery of CpG depletion in mitochondrial genome. Finally, we demonstrated that these insights may be integrated to better understand the pathogenicity of disease-implicated mitochondrial variants.

Rectifying long-standing misconceptions about the ρ statistic for molecular dating

When divided by a given mutation rate, the ρ (rho) statistic provides a simple estimator of the age of a clade within a phylogenetic tree by averaging the number of mutations from each sample in the clade to its root. However, a long-standing critique of the use of ρ in genetic dating has been quite often cited. Here we show that the critique is unfounded. We demonstrate by a formal mathematical argument and illustrate with a simulation study that ρ estimates are unbiased and also that ρ and maximum likelihood estimates do not differ in any systematic fashion. We also demonstrate that the claim that the associated confidence intervals commonly estimate the uncertainty inappropriately is flawed since it relies on a means of calculating standard errors that is not used by any other researchers, whereas an established expression for the standard error is largely unproblematic. We conclude that ρ dating, alongside approaches such as maximum likelihood (ML) and Bayesian inference, remains a useful tool for genetic dating.

Roles of the mitochondrial genetics in cancer metastasis: not to be ignored any longer

Mitochondrial DNA (mtDNA) encodes for only a fraction of the proteins that are encoded within the nucleus, and therefore has typically been regarded as a lesser player in cancer biology and metastasis. Accumulating evidence, however, supports an increased role for mtDNA impacting tumor progression and metastatic susceptibility. Unfortunately, due to this delay, there is a dearth of data defining the relative contributions of specific mtDNA polymorphisms (SNP), which leads to an inability to effectively use these polymorphisms to guide and enhance therapeutic strategies and diagnosis. In addition, evidence also suggests that differences in mtDNA impact not only the cancer cells but also the cells within the surrounding tumor microenvironment, suggesting a broad encompassing role for mtDNA polymorphisms in regulating the disease progression. mtDNA may have profound implications in the regulation of cancer biology and metastasis. However, there are still great lengths to go to understand fully its contributions. Thus, herein, we discuss the recent advances in our understanding of mtDNA in cancer and metastasis, providing a framework for future functional validation and discovery.

Ancient nuclear genomes enable repatriation of Indigenous human remains

After European colonization, the ancestral remains of Indigenous people were often collected for scientific research or display in museum collections. For many decades, Indigenous people, including Native Americans and Aboriginal Australians, have fought for their return. However, many of these remains have no recorded provenance, making their repatriation very difficult or impossible. To determine whether DNA-based methods could resolve this important problem, we sequenced 10 nuclear genomes and 27 mitogenomes from ancient pre-European Aboriginal Australians (up to 1540 years before the present) of known provenance and compared them to 100 high-coverage contemporary Aboriginal Australian genomes, also of known provenance. We report substantial ancient population structure showing strong genetic affinities between ancient and contemporary Aboriginal Australian individuals from the same geographic location. Our findings demonstrate the feasibility of successfully identifying the origins of unprovenanced ancestral remains using genomic methods.

Experimental evidence that thermal selection shapes mitochondrial genome evolution

Mitochondria are essential organelles, found within eukaryotic cells, which contain their own DNA. Mitochondrial DNA (mtDNA) has traditionally been used in population genetic and biogeographic studies as a maternally-inherited and evolutionary-neutral genetic marker. However, it is now clear that polymorphisms within the mtDNA sequence are routinely non-neutral, and furthermore several studies have suggested that such mtDNA polymorphisms are also sensitive to thermal selection. These observations led to the formulation of the “mitochondrial climatic adaptation” hypothesis, for which all published evidence to date is correlational. Here, we use laboratory-based experimental evolution in the fruit fly, Drosophila melanogaster, to test whether thermal selection can shift population frequencies of two mtDNA haplogroups whose natural frequencies exhibit clinal associations with latitude along the Australian east-coast. We present experimental evidence that the thermal regime in which the laboratory populations were maintained drove changes in haplogroup frequencies across generations. Our results strengthen the emerging view that intra-specific mtDNA variants are sensitive to selection, and suggest spatial distributions of mtDNA variants in natural populations of metazoans might reflect adaptation to climatic environments rather than within-population coalescence and diffusion of selectively-neutral haplotypes across populations.

Evolutionary analysis of mitochondrially encoded proteins of toad-headed lizards, Phrynocephalus, along an altitudinal gradient

Background

Animals living at high altitude must adapt to environments with hypoxia and low temperatures, but relatively little is known about underlying genetic changes. Toad-headed lizards of the genus Phrynocephalus cover a broad altitudinal gradient of over 4000 m and are useful models for studies of such adaptive responses. In one of the first studies to have considered selection on mitochondrial protein-coding regions in an ectothermic group distributed over such a wide range of environments, we analysed nineteen complete mitochondrial genomes from all Chinese Phrynocephalus (including eight genomes sequenced for the first time). Initial analyses used site and branch-site model (program: PAML) approaches to examine nonsynonymous: synonymous substitution rates across the mtDNA tree.

Results

Ten positively selected sites were discovered, nine of which corresponded to subunits ND2, ND3, ND4, ND5, and ND6 within the respiratory chain enzyme mitochondrial Complex I (NADH Coenzyme Q oxidoreductase). Four of these sites showed evidence of general long-term selection across the group while the remainder showed evidence of episodic selection across different branches of the tree. Some of these branches corresponded to increases in altitude and/or latitude. Analyses of physicochemical changes in protein structures revealed that residue changes at sites that were under selection corresponded to major functional differences. Analyses of coevolution point to coevolution of selected sites within the ND4 subunit, with key sites associated with proton translocation across the mitochondrial membrane.

Conclusions

Our results identify mitochondrial Complex I as a target for environment-mediated selection in this group of lizards, a complex that frequently appears to be under selection in other organisms. This makes these lizards good candidates for more detailed future studies of molecular evolution.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4569-1) contains supplementary material, which is available to authorized users.

Mitogenomic differences between the normal and tumor cells of colorectal cancer patients

So far, a reliable spectrum of mitochondrial DNA mutations in colorectal cancer cells is still unknown, and neither is their significance in carcinogenesis. Indeed, it remains debatable whether mtDNA mutations are "drivers" or "passengers" of colorectal carcinogenesis. Thus, we analyzed 200 mitogenomes from normal and cancer tissues of 100 colorectal cancer patients. Minority variant mutations were detected at the 1% level. We showed that somatic mutations frequently occur in colorectal cancer cells (75%) and are randomly distributed across the mitochondrial genome. Mutational signatures of somatic mitogenome mutations suggest that they might arise through nucleotide deamination due to oxidative stress. The majority of somatic mutations localized within the coding region (in positions not known from the human phylogeny) and was potentially pathogenic to cell metabolism. Further analysis suggested that the relaxation of negative selection in the mitogenomes of colorectal cancer cells may allow accumulation of somatic mutations. Thus, a shift in glucose metabolism from oxidative phosphorylation to glycolysis may create advantageous conditions for accumulation of mtDNA mutations. Considering the fact that the presence of somatic mtDNA mutations was not associated with any clinicopathological features, we suggested that mtDNA somatic mutations are "passengers" rather than the cause of colorectal carcinogenesis.

Mitochondrial DNA Heteroplasmy and Purifying Selection in the Mammalian Female Germ Line

Inherited mutations in the mitochondrial (mt)DNA are a major cause of human disease, with approximately 1 in 5000 people affected by one of the hundreds of identified pathogenic mtDNA point mutations or deletions. Due to the severe, and often untreatable, symptoms of many mitochondrial diseases, identifying how these mutations are inherited from one generation to the next has been an area of intense research in recent years. Despite large advances in our understanding of this complex process, many questions remain unanswered, with one of the most hotly debated being whether or not purifying selection acts against pathogenic mutations during germline development.

Ancient Human Migrations to and through Jammu Kashmir- India were not of Males Exclusively

Jammu and Kashmir (J&K), the Northern most State of India, has been under-represented or altogether absent in most of the phylogenetic studies carried out in literature, despite its strategic location in the Himalayan region. Nonetheless, this region may have acted as a corridor to various migrations to and from mainland India, Eurasia or northeast Asia. The belief goes that most of the migrations post-late-Pleistocene were mainly male dominated, primarily associated with population invasions, where female migration may thus have been limited. To evaluate female-centered migration patterns in the region, we sequenced 83 complete mitochondrial genomes of unrelated individuals belonging to different ethnic groups from the state. We observed a high diversity in the studied maternal lineages, identifying 19 new maternal sub-haplogroups (HGs). High maternal diversity and our phylogenetic analyses suggest that the migrations post-Pleistocene were not strictly paternal, as described in the literature. These preliminary observations highlight the need to carry out an extensive study of the endogamous populations of the region to unravel many facts and find links in the peopling of India.

Mutation in populations governed by a Galton-Watson branching process

A population genetics model based on a multitype branching process, or equivalently a Galton-Watson branching process for multiple alleles, is presented. The diffusion limit forward Kolmogorov equation is derived for the case of neutral mutations. The asymptotic stationary solution is obtained and has the property that the extant population partitions into subpopulations whose relative sizes are determined by mutation rates. An approximate time-dependent solution is obtained in the limit of low mutation rates. This solution has the property that the system undergoes a rapid transition from a drift-dominated phase to a mutation-dominated phase in which the distribution collapses onto the asymptotic stationary distribution. The changeover point of the transition is determined by the per-generation growth factor and mutation rate. The approximate solution is confirmed using numerical simulations.

Cellular and population level processes influence the rate, accumulation and observed frequency of inherited and somatic mtDNA mutations

Mitochondria are found in all animals and have the unique feature of containing multiple copies of their own small, circular DNA genome (mtDNA). The rate and pattern of mutation accumulation in the mtDNA are influenced by molecular, cellular and population level processes. We distinguish between inherited and somatic mtDNA mutations and review evidence for the often-made assumption that mutations accumulate at a higher rate in mtDNA than in nuclear DNA (nDNA). We conclude that the whole genome mutation accumulation rate is higher for mtDNA than for nDNA but include the caveat that rates overlap considerably between the individual mtDNA- and nDNA-encoded genes. Next, we discuss the postulated causal mechanisms for the high rate of mtDNA mutation accumulation in both inheritance and in somatic cells. Perhaps unexpectedly, mtDNA is resilient to many mutagens of nDNA but is prone to errors of replication. We then consider the influence of maternal inheritance, recombination and selection on the observed accumulation pattern of inherited mtDNA mutations. Finally, we discuss environmental influences of temperature and diet on the observed frequency of inherited and somatic mtDNA mutations. We conclude that it is necessary to understand the cellular processes to fully interpret the pattern of mutations and how they influence our interpretations of evolution and disease.

The aetiology of cardiovascular disease: a role for mitochondrial DNA?

Cardiovascular disease (CVD) is a world-wide cause of mortality in humans and its incidence is on the rise in Africa. In this review, we discuss the putative role of mitochondrial dysfunction in the aetiology of CVD and consequently identify mitochondrial DNA (mtDNA) variation as a viable genetic risk factor to be considered. We then describe the contribution and pitfalls of several current approaches used when investigating mtDNA in relation to complex disease. We also propose an alternative approach, the adjusted mutational load hypothesis, which would have greater statistical power with cohorts of moderate size, and is less likely to be affected by population stratification. We therefore address some of the shortcomings of the current haplogroup association approach. Finally, we discuss the unique challenges faced by studies done on African populations, and recommend the most viable methods to use when investigating mtDNA variation in CVD and other common complex disease.

Characterization of complete mitochondrial genomes of indigenous Mayans in Mexico

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Aboriginal Australian mitochondrial genome variation - an increased understanding of population antiquity and diversity

Aboriginal Australians represent one of the oldest continuous cultures outside Africa, with evidence indicating that their ancestors arrived in the ancient landmass of Sahul (present-day New Guinea and Australia) ~55 thousand years ago. Genetic studies, though limited, have demonstrated both the uniqueness and antiquity of Aboriginal Australian genomes. We have further resolved known Aboriginal Australian mitochondrial haplogroups and discovered novel indigenous lineages by sequencing the mitogenomes of 127 contemporary Aboriginal Australians. In particular, the more common haplogroups observed in our dataset included M42a, M42c, S, P5 and P12, followed by rarer haplogroups M15, M16, N13, O, P3, P6 and P8. We propose some major phylogenetic rearrangements, such as in haplogroup P where we delinked P4a and P4b and redefined them as P4 (New Guinean) and P11 (Australian), respectively. Haplogroup P2b was identified as a novel clade potentially restricted to Torres Strait Islanders. Nearly all Aboriginal Australian mitochondrial haplogroups detected appear to be ancient, with no evidence of later introgression during the Holocene. Our findings greatly increase knowledge about the geographic distribution and phylogenetic structure of mitochondrial lineages that have survived in contemporary descendants of Australia’s first settlers.

Selection on the mitochondrial ATP synthase 6 and the NADH dehydrogenase 2 genes in hares (Lepus capensis L., 1758) from a steep ecological gradient in North Africa

Background

Recent studies of selection on mitochondrial (mt) OXPHOS genes suggest adaptation due mainly to environmental variation. In this context, Tunisian hares that display several external phenotypes with phylogenetically rather homogenous gene pool and shallow population structure provide a good precondition to detect positive selection on mt genes related to environmental/climatic variation, specifically ambient temperature and precipitation.

Results

We used codon-based methods along with population genetic data to test for positive selection on ATPase synthase 6 (ATP6) and NADH dehydrogenase 2 (ND2) of cape hares (Lepus capensis) collected along a steep ecological gradient in Tunisia. We found significantly higher differentiation at the ATP6 locus across Tunisia, with sub-humid Mediterranean, semi-arid, and arid Sahara climate than for fourteen unlinked supposedly neutrally evolving nuclear microsatellites and mt control region sequences. This suggested positive selection on ATP6 sequences, which was confirmed by several codon-based tests for one sequence site that together with a second site translated into four different amino acids. Positive selection on ND2 sequences was also confirmed by several codon-based tests. The corresponding frequencies of the two most prevalent variants at each locus varied significantly across climate regions, and our logistic general linear models of occurrence of those proteins indicated significant effects of mean annual temperature for ATP6 and mean minimum temperature of the coldest month of the year for ND2, independent of geographical location, annual precipitation, and the respective co-occurring protein at the second locus. Moreover, presence of the ancestral ATP6 protein, as inferred from phylogenetic networks, was positively affected by the simultaneous presence of the derived ND2 protein and vice versa, independent of temperature, precipitation, or geographic location. Finally, we obtained a significant coevolution signal for the ancestral ATP6 and derived ND2 sequences and vice versa.

Conclusions

positive selection was strongly suggested by the population genetic approach and the codon-based tests in both mtDNA genes. Moreover, the two most prevalent proteins at the ATP6 locus were distributed at significantly varying frequencies across the study area with a significant effect of mean annual temperature on the occurrence of the ATP6 proteins independent of geographical coordinates and the co-occuring ND2 protein variant. For ND2, occurrence of the two most frequent protein variants was significantly influenced by the mean minimum temperature of the coldest month, independent of the co-occurring ATP6 protein variant and geographical coordinates. This strongly suggests direct involvement of ambient temperature in the adaptation of the studied mtOXPHOS genes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-0896-0) contains supplementary material, which is available to authorized users.

Mitochondrial DNA diversity of present-day Aboriginal Australians and implications for human evolution in Oceania

Aboriginal Australians are one of the more poorly studied populations from the standpoint of human evolution and genetic diversity. Thus, to investigate their genetic diversity, the possible date of their ancestors' arrival and their relationships with neighboring populations, we analyzed mitochondrial DNA (mtDNA) diversity in a large sample of Aboriginal Australians. Selected mtDNA single-nucleotide polymorphisms and the hypervariable segment haplotypes were analyzed in 594 Aboriginal Australians drawn from locations across the continent, chiefly from regions not previously sampled. Most (~78%) samples could be assigned to mtDNA haplogroups indigenous to Australia. The indigenous haplogroups were all ancient (with estimated ages >40 000 years) and geographically widespread across the continent. The most common haplogroup was P (44%) followed by S (23%) and M42a (9%). There was some geographic structure at the haplotype level. The estimated ages of the indigenous haplogroups range from 39 000 to 55 000 years, dates that fit well with the estimated date of colonization of Australia based on archeological evidence (~47 000 years ago). The distribution of mtDNA haplogroups in Australia and New Guinea supports the hypothesis that the ancestors of Aboriginal Australians entered Sahul through at least two entry points. The mtDNA data give no support to the hypothesis of secondary gene flow into Australia during the Holocene, but instead suggest long-term isolation of the continent.

Phylogenetic and population-based approaches to mitogenome variation do not support association with male infertility

Infertility has a complex multifactorial etiology and a high prevalence worldwide. Several studies have pointed to variation in the mitochondrial DNA (mtDNA) molecule as a factor responsible for the different disease phenotypes related to infertility. We analyzed 53 mitogenomes of infertile males from Galicia (northwest Spain), and these haplotypes were meta-analyzed phylogenetically with 43 previously reported from Portugal. Taking advantage of the large amount of information available, we additionally carried out association tests between patient mtDNA single-nucleotide polymorphisms (mtSNPs) and haplogroups against Iberian matched controls retrieved from The 1000 Genomes Project and the literature. Phylogenetic and association analyses did not reveal evidence of association between mtSNPs/haplogroups and infertility. Ratios and patterns in patients of nonsynonymous/synonymous changes, and variation at homoplasmic, heteroplasmic and private variants, fall within expected values for healthy individuals. Moreover, the haplogroup background of patients was variable and fits well with patterns typically observed in healthy western Europeans. We did not find evidence of association of mtSNPs or haplogroups pointing to a role for mtDNA in male infertility. A thorough review of the literature on mtDNA variation and infertility revealed contradictory findings and methodological and theoretical problems that overall undermine previous positive findings.

Mitochondrial coupling and capacity of oxidative phosphorylation in skeletal muscle of Inuit and Caucasians in the arctic winter

During evolution, mitochondrial DNA haplogroups of arctic populations may have been selected for lower coupling of mitochondrial respiration to ATP production in favor of higher heat production. We show that mitochondrial coupling in skeletal muscle of traditional and westernized Inuit habituating northern Greenland is identical to Danes of western Europe haplogroups. Biochemical coupling efficiency was preserved across variations in diet, muscle fiber type, and uncoupling protein-3 content. Mitochondrial phenotype displayed plasticity in relation to lifestyle and environment. Untrained Inuit and Danes had identical capacities to oxidize fat substrate in arm muscle, which increased in Danes during the 42 days of acclimation to exercise, approaching the higher level of the Inuit hunters. A common pattern emerges of mitochondrial acclimatization and evolutionary adaptation in humans at high latitude and high altitude where economy of locomotion may be optimized by preservation of biochemical coupling efficiency at modest mitochondrial density, when submaximum performance is uncoupled from VO2max and maximum capacities of oxidative phosphorylation.

MtDNA analysis reveals enriched pathogenic mutations in Tibetan highlanders

Tibetan highlanders, including Tibetans, Monpas, Lhobas, Dengs and Sherpas, are considered highly adaptive to severe hypoxic environments. Mitochondrial DNA (mtDNA) might be important in hypoxia adaptation given its role in coding core subunits of oxidative phosphorylation. In this study, we employed 549 complete highlander mtDNA sequences (including 432 random samples) to obtain a comprehensive view of highlander mtDNA profile. In the phylogeny of a total of 36,914 sequences, we identified 21 major haplogroups representing founding events of highlanders, most of which were coalesced in 10 kya. Through founder analysis, we proposed a three-phase model of colonizing the plateau, i.e., pre-LGM Time (30 kya, 4.68%), post-LGM Paleolithic Time (16.8 kya, 29.31%) and Neolithic Time (after 8 kya, 66.01% in total). We observed that pathogenic mutations occurred far more frequently in 22 highlander-specific lineages (five lineages carrying two pathogenic mutations and six carrying one) than in the 6,857 haplogroups of all the 36,914 sequences (P = 4.87 × 10⁻⁸). Furthermore, the number of possible pathogenic mutations carried by highlanders (in average 3.18 ± 1.27) were significantly higher than that in controls (2.82 ± 1.40) (P = 1.89 × 10⁻⁴). Considering that function-altering and pathogenic mutations are enriched in highlanders, we therefore hypothesize that they may have played a role in hypoxia adaptation.

Yeast models of mutations in the mitochondrial ATP6 gene found in human cancer cells

Since the discovery of somatic mtDNA mutations in tumor cells, multiple studies have focused on establishing a causal relationship between those changes and alterations in energy metabolism, a hallmark of cancer cells. Yet the consequences of these mutations on mitochondrial function remain largely unknown. In this study, Saccharomyces cerevisiae has been used as a model to investigate the functional consequences of four cancer-associated missense mutations (8914C>A, 8932C>T, 8953A>G, 9131T>C) found in the mitochondrial MT-ATP6 gene. This gene encodes the a-subunit of F1FO-ATP synthase, which catalyzes the last steps of ATP production in mitochondria. Although the four studied mutations affected well-conserved residues of the a-subunit, only one of them (8932C>T) had a significant impact on mitochondrial function, due to a less efficient incorporation of the a-subunit into ATP synthase. Our findings indicate that these ATP6 genetic variants found in human tumors are neutral mitochondrial genome substitutions with a limited, if any, impact on the energetic function of mitochondria.