Investigating mitonuclear interactions in human admixed populations

Mitochondrial haplotype and mito-nuclear matching drive somatic mutation and selection throughout ageing

Mitochondrial genomes co-evolve with the nuclear genome over evolutionary timescales and are shaped by selection in the female germline. Here we investigate how mismatching between nuclear and mitochondrial ancestry impacts the somatic evolution of the mitochondrial genome in different tissues throughout ageing. We used ultrasensitive duplex sequencing to profile ~2.5 million mitochondrial genomes across five mitochondrial haplotypes and three tissues in young and aged mice, cataloguing ~1.2 million mitochondrial somatic and ultralow-frequency inherited mutations, of which 81,097 are unique. We identify haplotype-specific mutational patterns and several mutational hotspots, including at the light strand origin of replication, which consistently exhibits the highest mutation frequency. We show that rodents exhibit a distinct mitochondrial somatic mutational spectrum compared with primates with a surfeit of reactive oxygen species-associated G > T/C > A mutations, and that somatic mutations in protein-coding genes exhibit signatures of negative selection. Lastly, we identify an extensive enrichment in somatic reversion mutations that 're-align' mito-nuclear ancestry within an organism's lifespan. Together, our findings demonstrate that mitochondrial genomes are a dynamically evolving subcellular population shaped by somatic mutation and selection throughout organismal lifetimes.

Hidden in the sand: Phylogenomics unravel an unexpected evolutionary history for the desert-adapted vipers of the genus Cerastes

The desert vipers of the genus Cerastes are a small clade of medically important venomous snakes within the family Viperidae. According to published morphological and molecular studies, the group is comprised by four species: two morphologically similar and phylogenetically sister taxa, the African horned viper (Cerastes cerastes) and the Arabian horned viper (Cerastes gasperettii); a more distantly related species, the Saharan sand viper (Cerastes vipera), and the enigmatic Böhme's sand viper (Cerastes boehmei), only known from a single specimen in captivity allegedly captured in Central Tunisia. In this study, we sequenced one mitochondrial marker (COI) as well as genome-wide data (ddRAD sequencing) from 28 and 41 samples, respectively, covering the entire distribution range of the genus to explore the population genomics, phylogenomic relationships and introgression patterns within the genus Cerastes. Additionally, and to provide insights into the mode of diversification of the group, we carried out niche overlap analyses considering climatic and habitat variables. Both nuclear phylogenomic reconstructions and population structure analyses have unveiled an unexpected evolutionary history for the genus Cerastes, which sharply contradicts the morphological similarities and previously published mitochondrial approaches. Cerastes cerastes and C. vipera are recovered as sister taxa whilst C. gasperettii is a sister taxon to the clade formed by these two species. We found a relatively high niche overlap (OI > 0.7) in both climatic and habitat variables between C. cerastes and C. vipera, contradicting a potential scenario of sympatric speciation. These results are in line with the introgression found between the northwestern African populations of C. cerastes and C. vipera. Finally, our genomic data confirms the existence of a lineage of C. cerastes in Arabia. All these results highlight the importance of genome-wide data over few genetic markers to study the evolutionary history of species.

Mitonuclear interactions modulate nutritional preference

In nature, organisms are faced with constant nutritional options which fuel key life-history traits. Studies have shown that species can actively make nutritional decisions based on internal and external cues. Metabolism itself is underpinned by complex genomic interactions involving components from both nuclear and mitochondrial genomes. Products from these two genomes must coordinate how nutrients are extracted, used and recycled. Given the complicated nature of metabolism, it is not well understood how nutritional choices are affected by mitonuclear interactions. This is under the rationale that changes in genomic interactions will affect metabolic flux and change physiological requirements. To this end we used a large Drosophila mitonuclear genetic panel, comprising nine isogenic nuclear genomes coupled to nine mitochondrial haplotypes, giving a total of 81 different mitonuclear genotypes. We use a capillary-based feeding assay to screen this panel for dietary preference between carbohydrate and protein. We find significant mitonuclear interactions modulating nutritional choices, with these epistatic interactions also being dependent on sex. Our findings support the notion that complex genomic interactions can place a constraint on metabolic flux. This work gives us deeper insights into how key metabolic interactions can have broad implications on behaviour.

Mitochondrial haplotype and mito-nuclear matching drive somatic mutation and selection throughout aging

Mitochondrial genomes co-evolve with the nuclear genome over evolutionary timescales and are shaped by selection in the female germline. Here, we investigate how mismatching between nuclear and mitochondrial ancestry impacts the somatic evolution of the mt-genome in different tissues throughout aging. We used ultra-sensitive Duplex Sequencing to profile ~2.5 million mt-genomes across five mitochondrial haplotypes and three tissues in young and aged mice, cataloging ~1.2 million mitochondrial somatic and ultra low frequency inherited mutations, of which 81,097 are unique. We identify haplotype-specific mutational patterns and several mutational hotspots, including at the Light Strand Origin of Replication, which consistently exhibits the highest mutation frequency. We show that rodents exhibit a distinct mitochondrial somatic mutational spectrum compared to primates with a surfeit of reactive oxygen species-associated G>T/C>A mutations, and that somatic mutations in protein coding genes exhibit signatures of negative selection. Lastly, we identify an extensive enrichment in somatic reversion mutations that "re-align" mito-nuclear ancestry within an organism's lifespan. Together, our findings demonstrate that mitochondrial genomes are a dynamically evolving subcellular population shaped by somatic mutation and selection throughout organismal lifetimes.

The genetic and phenotypic correlates of mtDNA copy number in a multi-ancestry cohort

Mitochondrial DNA copy number (mtCN) is often treated as a proxy for mitochondrial (dys-) function and disease risk. Pathological changes in mtCN are common symptoms of rare mitochondrial disorders, but reported associations between mtCN and common diseases vary across studies. To understand the biology of mtCN, we carried out genome- and phenome-wide association studies of mtCN in 30,666 individuals from the Penn Medicine BioBank (PMBB)-a diverse cohort of largely African and European ancestry. We estimated mtCN in peripheral blood using exome sequence data, taking cell composition into account. We replicated known genetic associations of mtCN in the PMBB and found that their effects are highly correlated between individuals of European and African ancestry. However, the heritability of mtCN was much higher among individuals of largely African ancestry ( h 2 = 0.3 ) compared with European ancestry individuals( h 2 = 0.1 ) . Admixture mapping suggests that there are undiscovered variants underlying mtCN that are differentiated in frequency between individuals with African and European ancestry. We show that mtCN is associated with many health-related phenotypes. We discovered robust associations between mtDNA copy number and diseases of metabolically active tissues, such as cardiovascular disease and liver damage, that were consistent across African and European ancestry individuals. Other associations, such as epilepsy and prostate cancer, were only discovered in either individuals with European or African ancestry but not both. We show that mtCN-phenotype associations can be sensitive to blood cell composition and environmental modifiers, explaining why such associations are inconsistent across studies. Thus, mtCN-phenotype associations must be interpreted with care.

G × G × E effect on phenotype expression in a non-conventional model organism, the unicellular slime mould Physarum polycephalum

In metazoans, the expression of key phenotypic traits is sensitive to two- and three-way interactions between variation in mitochondrial DNA, nuclear DNA and the external environment. Whether gene-by-environment interactions affect phenotypes in single-celled eukaryotes is poorly studied, except in a few species of yeast and fungi. We developed a genetic panel of the unicellular slime mould, Physarum polycephalum containing strains differing in mitochondrial and nuclear DNA haplotypes. The panel also included two strains harbouring a selfishly replicating mitochondrial-fusion (mF) plasmid that could affect phenotype expression. We assayed movement and growth rate differences among the strains across two temperature regimes: 24° and 28°C. We found that the slime mould's growth rate, but not movement, is affected by G × G × E interactions. Predictably, mtDNA × nDNA interactions significantly affected both traits. The inter-trait correlation across the strains in each temperature regime was positive. Surprisingly, the mF plasmid had no negative effects on our chosen traits. Our study is the first to demonstrate genetic regulation of phenotype expression in a unicellular slime mould. The genetic effect on phenotypes manifests via epistatic interactions with the thermal environment, thus shedding new light on the role of G × G × E interactions in trait evolution in protists.

MultiWaverX: modeling latent sex-biased admixture history

Sex-biased gene flow has been common in the demographic history of modern humans. However, the lack of sophisticated methods for delineating the detailed sex-biased admixture process prevents insights into complex admixture history and thus our understanding of the evolutionary mechanisms of genetic diversity. Here, we present a novel algorithm, MultiWaverX, for modeling complex admixture history with sex-biased gene flow. Systematic simulations showed that MultiWaverX is a powerful tool for modeling complex admixture history and inferring sex-biased gene flow. Application of MultiWaverX to empirical data of 17 typical admixed populations in America, Central Asia, and the Middle East revealed sex-biased admixture histories that were largely consistent with the historical records. Notably, fine-scale admixture process reconstruction enabled us to recognize latent sex-biased gene flow in certain populations that would likely be overlooked by much of the routine analysis with commonly used methods. An outstanding example in the real world is the Kazakh population that experienced complex admixture with sex-biased gene flow but in which the overall signature has been canceled due to biased gene flow from an opposite direction.

Exploring the Effects of Mitonuclear Interactions on Mitochondrial DNA Gene Expression in Humans

Most mitochondrial protein complexes include both nuclear and mitochondrial gene products, which coevolved to work together. This coevolution can be disrupted due to disparity in genetic ancestry between the nuclear and mitochondrial genomes in recently admixed populations. Such mitonuclear DNA discordance might result in phenotypic effects. Several nuclear-encoded proteins regulate expression of mitochondrial DNA (mtDNA) genes. We hypothesized that mitonuclear DNA discordance affects expression of genes encoded by mtDNA. To test this, we utilized the data from the GTEx project, which contains expression levels for ∼100 African Americans and >600 European Americans. The varying proportion of African and European ancestry in recently admixed African Americans provides a range of mitonuclear discordance values, which can be correlated with mtDNA gene expression levels (adjusted for age and ischemic time). In contrast, European Americans did not undergo recent admixture. We demonstrated that, for most mtDNA protein-coding genes, expression levels in energetically-demanding tissues were lower in African Americans than in European Americans. Furthermore, gene expression levels were lower in individuals with higher mitonuclear discordance, independent of population. Moreover, we found a negative correlation between mtDNA gene expression and mitonuclear discordance. In African Americans, the average value of African ancestry was higher for nuclear-encoded mitochondrial than non-mitochondrial genes, facilitating a match in ancestry with the mtDNA and more optimal interactions. These results represent an example of a phenotypic effect of mitonuclear discordance on human admixed populations, and have potential biomedical applications.

Mitochondrial genome variations, mitochondrial-nuclear compatibility, and their association with metabolic diseases

Two genomes regulate the energy metabolism of eukaryotic cells: the nuclear genome, which codes for most cellular proteins, and the mitochondrial genome, which, together with the nuclear genome, coregulates cellular bioenergetics. Therefore, mitochondrial genome variations can affect, directly or indirectly, all energy-dependent cellular processes and shape the metabolic state of the organism. This review provides a current and up-to-date overview on how codependent these two genomes are, how they appear to have coevolved, and how variations within the mitochondrial genome might be associated with the manifestation of metabolic diseases. This review summarizes and structures results obtained from epidemiological studies that identified links between mitochondrial haplogroups and individual risks for developing obesity and diabetes. This is complemented by findings on the compatibility of mitochondrial and nuclear genomes and cellular bioenergetic fitness, which have been acquired from well-controlled studies in conplastic animal models. These elucidate, more mechanistically, how single-nucleotide variants can influence cellular metabolism and physiology. Overall, it seems that certain mitochondrial genome variations negatively affect mitochondrial-nuclear compatibility and are statistically linked with the onset of metabolic diseases, whereas, for others, greater uncertainty exists, and additional research into this exciting field is required.

Genome-wide local ancestry and evidence for mitonuclear coadaptation in African hybrid cattle populations (Bos taurus/indicus)

The phenotypic diversity of African cattle reflects adaptation to a wide range of agroecological conditions, human-mediated selection preferences, and complex patterns of admixture between the humpless Bos taurus (taurine) and humped Bos indicus (zebu) subspecies, which diverged 150-500 thousand years ago. Despite extensive admixture, all African cattle possess taurine mitochondrial haplotypes, even populations with significant zebu biparental and male uniparental nuclear ancestry. This has been interpreted as the result of human-mediated dispersal ultimately stemming from zebu bulls imported from South Asia during the last three millennia. Here, we assess whether ancestry at mitochondrially targeted nuclear genes in African admixed cattle is impacted by mitonuclear functional interactions. Using high-density SNP data, we find evidence for mitonuclear coevolution across hybrid African cattle populations with a significant increase of taurine ancestry at mitochondrially targeted nuclear genes. Our results, therefore, support the hypothesis of incompatibility between the taurine mitochondrial genome and the zebu nuclear genome.

Genome and cuticular hydrocarbon-based species delimitation shed light on potential drivers of speciation in a Neotropical ant species complex

Geographic separation that leads to the evolution of reproductive isolation between populations generally is considered the most common form of speciation. However, speciation may also occur in the absence of geographic barriers due to phenotypic and genotypic factors such as chemical cue divergence, mating signal divergence, and mitonuclear conflict. Here, we performed an integrative study based on two genome-wide techniques (3RAD and ultraconserved elements) coupled with cuticular hydrocarbon (CHC) and mitochondrial (mt) DNA sequence data, to assess the species limits within the Ectatomma ruidum species complex, a widespread and conspicuous group of Neotropical ants for which heteroplasmy (i.e., presence of multiple mtDNA variants in an individual) has been recently discovered in some populations from southeast Mexico. Our analyses indicate the existence of at least five distinct species in this complex: two widely distributed across the Neotropics, and three that are restricted to southeast Mexico and that apparently have high levels of heteroplasmy. We found that species boundaries in the complex did not coincide with geographic barriers. We therefore consider possible roles of alternative drivers that may have promoted the observed patterns of speciation, including mitonuclear incompatibility, CHC differentiation, and colony structure. Our study highlights the importance of simultaneously assessing different sources of evidence to disentangle the species limits of taxa with complicated evolutionary histories.

Phylogenetic inference of changes in amino acid propensities with single-position resolution

Fitness conferred by the same allele may differ between genotypes and environments, and these differences shape variation and evolution. Changes in amino acid propensities at protein sites over the course of evolution have been inferred from sequence alignments statistically, but the existing methods are data-intensive and aggregate multiple sites. Here, we develop an approach to detect individual amino acids that confer different fitness in different groups of species from combined sequence and phylogenetic data. Using the fact that the probability of a substitution to an amino acid depends on its fitness, our method looks for amino acids such that substitutions to them occur more frequently in one group of lineages than in another. We validate our method using simulated evolution of a protein site under different scenarios and show that it has high specificity for a wide range of assumptions regarding the underlying changes in selection, while its sensitivity differs between scenarios. We apply our method to the env gene of two HIV-1 subtypes, A and B, and to the HA gene of two influenza A subtypes, H1 and H3, and show that the inferred fitness changes are consistent with the fitness differences observed in deep mutational scanning experiments. We find that changes in relative fitness of different amino acid variants within a site do not always trigger episodes of positive selection and therefore may not result in an overall increase in the frequency of substitutions, but can still be detected from changes in relative frequencies of different substitutions.

Why do we pick similar mates, or do we?

Humans often mate with those resembling themselves, a phenomenon described as positive assortative mating (PAM). The causes of this attract broad interest, but there is little agreement on the topic. This may be because empirical studies and reviews sometimes focus on just a few explanations, often based on disciplinary conventions. This review presents an interdisciplinary conceptual framework on the causes of PAM in humans, drawing on human and non-human biology, the social sciences, and the humanities. Viewing causality holistically, we first discuss the proximate causes (i.e. the 'how') of PAM, considering three mechanisms: stratification, convergence and mate choice. We also outline methods to control for confounders when studying mate choice. We then discuss ultimate explanations (i.e. 'the why') for PAM, including adaptive and non-adaptive processes. We conclude by suggesting a focus on interdisciplinarity in future research.

Mitonuclear mismatch alters nuclear gene expression in naturally introgressed Rhinolophus bats

Background

Mitochondrial function involves the interplay between mitochondrial and nuclear genomes. Such mitonuclear interactions can be disrupted by the introgression of mitochondrial DNA between taxa or divergent populations. Previous studies of several model systems (e.g. Drosophila) indicate that the disruption of mitonuclear interactions, termed mitonuclear mismatch, can alter nuclear gene expression, yet few studies have focused on natural populations.

Results

Here we study a naturally introgressed population in the secondary contact zone of two subspecies of the intermediate horseshoe bat (Rhinolophus affinis), in which individuals possess either mitonuclear matched or mismatched genotypes. We generated transcriptome data for six tissue types from five mitonuclear matched and five mismatched individuals. Our results revealed strong tissue-specific effects of mitonuclear mismatch on nuclear gene expression with the largest effect seen in pectoral muscle. Moreover, consistent with the hypothesis that genes associated with the response to oxidative stress may be upregulated in mitonuclear mismatched individuals, we identified several such gene candidates, including DNASE1L3, GPx3 and HSPB6 in muscle, and ISG15 and IFI6 in heart.

Conclusion

Our study reveals how mitonuclear mismatch arising from introgression in natural populations is likely to have fitness consequences. Underlying the processes that maintain mitonuclear discordance is a step forward to understand the role of mitonuclear interactions in population divergence and speciation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12983-021-00424-x.

Validation of low-coverage whole-genome sequencing for mitochondrial DNA variants suggests mitochondrial DNA as a genetic cause of preterm birth

Preterm birth (PTB), or birth that occurs earlier than 37 weeks of gestational age, is a major contributor to infant mortality and neonatal hospitalization. Mutations in the mitochondrial genome (mtDNA) have been linked to various rare mitochondrial disorders and may be a contributing factor in PTB given that maternal genetic factors have been strongly linked to PTB. However, to date, no study has found a conclusive connection between a particular mtDNA variant and PTB. Given the high mtDNA copy number per cell, an automated pipeline was developed for detecting mtDNA variants using low‐coverage whole‐genome sequencing (lcWGS) data. The pipeline was first validated against samples of known heteroplasmy, and then applied to 929 samples from a PTB cohort from diverse ethnic backgrounds with an average gestational age of 27.18 weeks (range: 21–30). Our new pipeline successfully identified haplogroups and a large number of mtDNA variants in this large PTB cohort, including 8 samples carrying known pathogenic variants and 47 samples carrying rare mtDNA variants. These results confirm that lcWGS can be utilized to reliably identify mtDNA variants. These mtDNA variants may make a contribution toward preterm birth in a small proportion of live births.

Boosting the potential of cattle breeding using molecular biology, genetics, and bioinformatics approaches – a review

Cattle are among the most important farm animals that underwent an intense selection with the aim to increase milk production and to improve growth and meat properties, meanwhile reducing the generation interval allowing for a faster herd turnover. Recently, a shift from traditional breeding methods to breeding based on genetic testing has been observed. In this perspective, we review the techniques of molecular biology, genetics, and bioinformatics that are expected to further boost the agricultural potential of cattle. We discuss embryo selection based on next-generation and Nanopore sequencing and in vitro embryo production, boosting the potential of genetically superior animals. Gene editing of embryos could further speed up the selection process, essentially introducing a change in a single generation. Lastly, we discuss the host-microbiome co-evolution and adaptation. For example, cattle already adapted to low-quality low-cost fodder could be bred to achieve desired properties for the beef and dairy industry. The challenge of breeding and genetic editing is to accompany the selection on desired consumer-oriented traits with the push for sustainability and the adaptation to a changing climate while remaining economically viable. We propose that we are yet to see the limits of what is possible to achieve with modern technology for the cattle of the future; the ultimate goal will be to produce and maintain genetically elite individuals that can sustain the growing demands on the production.

The Role of Mitonuclear Incompatibility in Bipolar Disorder Susceptibility and Resilience Against Environmental Stressors

It has been postulated that mitochondrial dysfunction has a significant role in the underlying pathophysiology of bipolar disorder (BD). Mitochondrial functioning plays an important role in regulating synaptic transmission, brain function, and cognition. Neuronal activity is energy dependent and neurons are particularly sensitive to changes in bioenergetic fluctuations, suggesting that mitochondria regulate fundamental aspects of brain function. Vigorous evidence supports the role of mitochondrial dysfunction in the etiology of BD, including dysregulated oxidative phosphorylation, general decrease of energy, altered brain bioenergetics, co-morbidity with mitochondrial disorders, and association with genetic variants in mitochondrial DNA (mtDNA) or nuclear-encoded mitochondrial genes. Despite these advances, the underlying etiology of mitochondrial dysfunction in BD is unclear. A plausible evolutionary explanation is that mitochondrial-nuclear (mitonuclear) incompatibility leads to a desynchronization of machinery required for efficient electron transport and cellular energy production. Approximately 1,200 genes, encoded from both nuclear and mitochondrial genomes, are essential for mitochondrial function. Studies suggest that mitochondrial and nuclear genomes co-evolve, and the coordinated expression of these interacting gene products are essential for optimal organism function. Incompatibilities between mtDNA and nuclear-encoded mitochondrial genes results in inefficiency in electron flow down the respiratory chain, differential oxidative phosphorylation efficiency, increased release of free radicals, altered intracellular Ca²⁺ signaling, and reduction of catalytic sites and ATP production. This review explores the role of mitonuclear incompatibility in BD susceptibility and resilience against environmental stressors.

Pervasive Mitonuclear Coadaptation Underlies Fast Development in Interpopulation Hybrids of a Marine Crustacean

Cellular energy production requires coordinated interactions between genetic components from the nuclear and mitochondrial genomes. This coordination results in coadaptation of interacting elements within populations. Interbreeding between divergent gene pools can disrupt coadapted loci and result in hybrid fitness breakdown. While specific incompatible loci have been detected in multiple eukaryotic taxa, the extent of the nuclear genome that is influenced by mitonuclear coadaptation is not clear in any species. Here, we used F₂ hybrids between two divergent populations of the copepod Tigriopus californicus to examine mitonuclear coadaptation across the nuclear genome. Using developmental rate as a measure of fitness, we found that fast-developing copepods had higher ATP synthesis capacity than slow developers, suggesting variation in developmental rates is at least partly associated with mitochondrial dysfunction. Using Pool-seq, we detected strong biases for maternal alleles across 7 (of 12) chromosomes in both reciprocal crosses in high-fitness hybrids, whereas low-fitness hybrids showed shifts toward the paternal population. Comparison with previous results on a different hybrid cross revealed largely different patterns of strong mitonuclear coadaptation associated with developmental rate. Our findings suggest that functional coadaptation between interacting nuclear and mitochondrial components is reflected in strong polygenic effects on this life-history phenotype, and reveal that molecular coadaptation follows independent evolutionary trajectories among isolated populations.

Contribution of Mitochondrial DNA Heteroplasmy to the Congenital Cardiac and Palatal Phenotypic Variability in Maternally Transmitted 22q11.2 Deletion Syndrome

Congenital heart disease (CHD) and palatal anomalies (PA), are among the most common characteristics of 22q11.2 deletion syndrome (22q11.2DS), but they show incomplete penetrance, suggesting the presence of additional factors. The 22q11.2 deleted region contains nuclear encoded mitochondrial genes, and since mitochondrial function is critical during development, we hypothesized that changes in the mitochondrial DNA (mtDNA) could be involved in the intrafamilial variability of CHD and PA in cases of maternally inherited 22q11.2DS. To investigate this, we studied the transmission of heteroplasmic mtDNA alleles in seventeen phenotypically concordant and discordant mother-offspring 22q11.2DS pairs. We sequenced their mtDNA and identified 26 heteroplasmic variants at >1% frequency, representing 18 transmissions. The median allele frequency change between a mother and her child was twice as much, with a wider distribution range, in PA discordant pairs, p-value = 0.039 (permutation test, 11 concordant vs. 7 discordant variants), but not in CHD discordant pairs, p-value = 0.441 (9 vs. 9). Only the variant m.9507T>C was considered to be pathogenic, but it was unrelated to the structural phenotypes. Our study is novel, yet our results are not consistent with mtDNA variation contributing to PA or CHD in 22q11.2DS. Larger cohorts and additional factors should be considered moving forward.

Mitochondrial DNA Replacement Techniques to Prevent Human Mitochondrial Diseases

Background: Mitochondrial DNA (mtDNA) diseases are a group of maternally inherited genetic disorders caused by a lack of energy production. Currently, mtDNA diseases have a poor prognosis and no known cure. The chance to have unaffected offspring with a genetic link is important for the affected families, and mitochondrial replacement techniques (MRTs) allow them to do so. MRTs consist of transferring the nuclear DNA from an oocyte with pathogenic mtDNA to an enucleated donor oocyte without pathogenic mtDNA. This paper aims to determine the efficacy, associated risks, and main ethical and legal issues related to MRTs. Methods: A bibliographic review was performed on the MEDLINE and Web of Science databases, along with searches for related clinical trials and news. Results: A total of 48 publications were included for review. Five MRT procedures were identified and their efficacy was compared. Three main risks associated with MRTs were discussed, and the ethical views and legal position of MRTs were reviewed. Conclusions: MRTs are an effective approach to minimizing the risk of transmitting mtDNA diseases, but they do not remove it entirely. Global legal regulation of MRTs is required.

Investigating Mitonuclear Genetic Interactions Through Machine Learning: A Case Study on Cold Adaptation Genes in Human Populations From Different European Climate Regions

Cold climates represent one of the major environmental challenges that anatomically modern humans faced during their dispersal out of Africa. The related adaptive traits have been achieved by modulation of thermogenesis and thermoregulation processes where nuclear (nuc) and mitochondrial (mt) genes play a major role. In human populations, mitonuclear genetic interactions are the result of both the peculiar genetic history of each human group and the different environments they have long occupied. This study aims to investigate mitonuclear genetic interactions by considering all the mitochondrial genes and 28 nuclear genes involved in brown adipose tissue metabolism, which have been previously hypothesized to be crucial for cold adaptation. For this purpose, we focused on three human populations (i.e., Finnish, British, and Central Italian people) of European ancestry from different biogeographical and climatic areas, and we used a machine learning approach to identify relevant nucDNA–mtDNA interactions that characterized each population. The obtained results are twofold: (i) at the methodological level, we demonstrated that a machine learning approach is able to detect patterns of genetic structure among human groups from different latitudes both at single genes and by considering combinations of mtDNA and nucDNA loci; (ii) at the biological level, the analysis identified population-specific nuclear genes and variants that likely play a relevant biological role in association with a mitochondrial gene (such as the “obesity gene” FTO in Finnish people). Further studies are needed to fully elucidate the evolutionary dynamics (e.g., migration, admixture, and/or local adaptation) that shaped these nucDNA–mtDNA interactions and their functional role.

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.

Disrupted architecture and fast evolution of the mitochondrial genome of Argeia pugettensis (Isopoda): implications for speciation and fitness

BACKGROUND

Argeia pugettensis is an isopod species that parasitizes other crustaceans. Its huge native geographic range spans the Pacific from China to California, but molecular data are available only for a handful of specimens from North-American populations. We sequenced and characterised the complete mitogenome of a specimen collected in the Yellow Sea.

RESULTS

It exhibited a barcode (cox1) similarity level of only 87-89% with North-American populations, which is unusually low for conspecifics. Its mitogenome is among the largest in isopods (≈16.5 Kbp), mostly due to a large duplicated palindromic genomic segment (2 Kbp) comprising three genes. However, it lost a segment comprising three genes, nad4L-trnP-nad6, and many genes exhibited highly divergent sequences in comparison to isopod orthologues, including numerous mutations, deletions and insertions. Phylogenetic and selection analyses corroborated that this is one of the handful of most rapidly evolving available isopod mitogenomes, and that it evolves under highly relaxed selection constraints (as opposed to positive selection). However, its nuclear 18S gene is highly conserved, which suggests that rapid evolution is limited to its mitochondrial genome. The cox1 sequence analysis indicates that elevated mitogenomic evolutionary rates are not shared by North-American conspecifics, which suggests a breakdown of cox1 barcoding in this species.

CONCLUSIONS

A highly architecturally disrupted mitogenome and decoupling of mitochondrial and nuclear rates would normally be expected to have strong negative impacts on the fitness of the organism, so the existence of this lineage is a puzzling evolutionary question. Additional studies are needed to assess the phylogenetic breadth of this disrupted mitochondrial architecture and its impact on fitness.

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.

Y Haplogroup Diversity of the Dominican Republic: Reconstructing the Effect of the European Colonization and the Trans-Atlantic Slave Trades

The Dominican Republic is one of the two countries on the Hispaniola island, which is part of the Antilles. Hispaniola was affected by the European colonization and massive deportation of African slaves since the XVI century and these events heavily shaped the genetic composition of the present-day population. To shed light about the effect of the European rules, we analyzed 92 single nucleotide polymorphisms on the Y chromosome in 182 Dominican individuals from three different locations. The Dominican Y haplogroup composition was characterized by an excess of northern African/European lineages (59%), followed by the African clades (38%), whereas the Native-American lineages were rare (3%). The comparison with the mitochondrial DNA variability, dominated by African clades, revealed a sex-biased admixture pattern, in line with the colonial society dominated by European men. When other Caribbean and non-Caribbean former colonies were also considered, we noted a difference between territories under a Spanish rule (like the Dominican Republic) and British/French rule, with the former characterized by an excess of European Y lineages reflecting the more permissive Iberian legislation about mixed people and slavery. Finally, we analyzed the distribution in Africa of the Dominican lineages with a putative African origin, mainly focusing on central and western Africa, which were the main sources of African slaves. We found that most (83%) of the African lineages observed in Santo Domingo have a central African ancestry, suggesting that most of the slaves were deported from regions.

Multi-locus interactions and the build-up of reproductive isolation

All genes interact with other genes, and their additive effects and epistatic interactions affect an organism's phenotype and fitness. Recent theoretical and empirical work has advanced our understanding of the role of multi-locus interactions in speciation. However, relating different models to one another and to empirical observations is challenging. This review focuses on multi-locus interactions that lead to reproductive isolation (RI) through reduced hybrid fitness. We first review theoretical approaches and show how recent work incorporating a mechanistic understanding of multi-locus interactions recapitulates earlier models, but also makes novel predictions concerning the build-up of RI. These include high variance in the build-up rate of RI among taxa, the emergence of strong incompatibilities producing localized barriers to introgression, and an effect of population size on the build-up of RI. We then review recent experimental approaches to detect multi-locus interactions underlying RI using genomic data. We argue that future studies would benefit from overlapping methods like ancestry disequilibrium scans, genome scans of differentiation and analyses of hybrid gene expression. Finally, we highlight a need for further overlap between theoretical and empirical work, and approaches that predict what kind of patterns multi-locus interactions resulting in incompatibilities will leave in genome-wide polymorphism data. This article is part of the theme issue 'Towards the completion of speciation: the evolution of reproductive isolation beyond the first barriers'.

The combined effect of nuclear and mitochondrial genomes on the risk of developing multiple sclerosis

Multiple sclerosis (MS) is a severe chronic CNS disease characterized by autoimmune inflammation, demyelination, and neurodegeneration. The interaction of mitochondrial and nuclear genomes is shown to be important in the formation of a predisposition to many diseases.

Objective: to analyze the association of MS with the carriage of biallelic combinations, including as components the polymorphisms of three genes of mitochondrial DNA (mtDNA) and those of 16 nuclear genes, the products of which are involved in the functioning of the immune system and may participate in the development of autoimmune inflammation in MS; and, if these combinations are identified, to determine the nature of an interaction between their components.

Patients and methods. The investigation enrolled 540 MS patients and 406 control group individuals; all were Russians. The mitochondrial genome was genotyped by polymerase chain reaction-restriction fragment length polymorphism analysis. APSampler software was used for multilocus association analysis. 

Results and discussion. The investigators identified five biallelic combinations that were associated with MS (p=0.0036–0.022) and possessed protective properties (odds ratio (OR) 0.67–0.75). The mitochondrial component of the identified combinations was the polymorphisms m.4580 (rs28357975), m.13368 (rs3899498), and m.13708 (rs28359178) mtDNA; the nuclear component was CXCR5 (rs523604), TNFRSF1A (rs1800693), and CD86 (rs2255214) gene polymorphisms. The interaction between the components of the identified combinations was additive.

Conclusion. The data obtained in the Russian population suggest that the combined contribution of the mitochondrial and nuclear genomes may affect the risk of developing MS.

The history behind the mosaic of the Americas

Focusing on literature published in 2018-2020, we review inferences about: (i) how ancient DNA is contributing to clarify the peopling of the Americas and the dispersal of its first inhabitants, (ii) how the interplay between environmental diversity and culture has influenced the genetic structure and adaptation of Andean and Amazon populations, (iii) how genetics has contributed to our understanding of the Pre-Columbian Tupi expansion in Eastern South America, (iv) the subcontinental origins and dynamics of Post-Columbian admixture in the Americas, and finally, (v) episodes of adaptive natural selection in the American continent, particularly in the high altitudes of the Andes.

Genetic and phenotypic landscape of the mitochondrial genome in the Japanese population

The genetic landscape of mitochondrial DNA (mtDNA) has been elusive. By analyzing mtDNA using the whole genome sequence (WGS) of Japanese individuals (n = 1928), we identified 2023 mtDNA variants and high-resolution haplogroups. Frequency spectra of the haplogroups were population-specific and were heterogeneous among geographic regions within Japan. Application of machine learning methods could finely classify the subjects corresponding to the high-digit mtDNA sub-haplogroups. mtDNA had distinct genetic structures from that of nuclear DNA (nDNA), characterized by no distance-dependent linkage disequilibrium decay, sparse tagging of common variants, and the existence of common haplotypes spanning the entire mtDNA. We did not detect any evidence of mtDNA–nDNA (or mtDNA copy number–nDNA) genotype associations. Together with WGS-based mtDNA variant imputation, we conducted a phenome-wide association study of 147,437 Japanese individuals with 99 clinical phenotypes. We observed pleiotropy of mtDNA genetic risk on the five late-onset human complex traits including creatine kinase (P = 1.7 × 10⁻¹²).

Kenichi Yamamoto et al. report a genetic analysis of mitochondrial DNA (mtDNA) and a phenome-wide association study in Japanese individuals from the BioBank Japan Project. They describe the genetic landscape of the mitochondria and identify pleiotropic mtDNA variants associated with 5 late-onset complex traits.

Mitonuclear genomics and aging

Our cells operate based on two distinct genomes that are enclosed in the nucleus and mitochondria. The mitochondrial genome presumably originates from endosymbiotic bacteria. With time, a large portion of the original genes in the bacterial genome is considered to have been lost or transferred to the nuclear genome, leaving a reduced 16.5 Kb circular mitochondrial DNA (mtDNA). Traditionally only 37 genes, including 13 proteins, were thought to be encoded within mtDNA, its genetic repertoire is expanding with the identification of mitochondrial-derived peptides (MDPs). The biology of aging has been largely unveiled to be regulated by genes that are encoded in the nuclear genome, whereas the mitochondrial genome remained more cryptic. However, recent studies position mitochondria and mtDNA as an important counterpart to the nuclear genome, whereby the two organelles constantly regulate each other. Thus, the genomic network that regulates lifespan and/or healthspan is likely constituted by two unique, yet co-evolved, genomes. Here, we will discuss aspects of mitochondrial biology, especially mitochondrial communication that may add substantial momentum to aging research by accounting for both mitonuclear genomes to more comprehensively and inclusively map the genetic and molecular networks that govern aging and age-related diseases.

Mitonuclear interactions influence Alzheimer's disease risk

We examined the associations between mitochondrial DNA haplogroups (MT-hgs; mitochondrial haplotype groups defined by a specific combination of single nucleotide polymorphisms labeled as letters running from A to Z) and their interactions with a polygenic risk score composed of nuclear-encoded mitochondrial genes (nMT-PRS) with risk of dementia and age of onset (AOO) of dementia. MT-hg K (Odds ratio [OR]: 2.03 [95% CI: 1.04, 3.97]) and a 1 SD larger nMT-PRS (OR: 2.2 [95% CI: 1.68, 2.86]) were associated with elevated odds of dementia. Significant antagonistic interactions between the nMT-PRS and MT-hg K (OR: 0.45 [95% CI: 0.22, 0.9]) and MT-hg T (OR: 0.22 [95% CI: 0.1, 0.49]) were observed. Individual MT-hgs were not associated with AOO; however, a significant antagonistic interactions was observed between the nMT-PRS and MT-hg T (Hazard ratio: 0.62 [95% CI: 0.42, 0.91]) and a synergistic interaction between the nMT-PRS and MT-hg V (Hazard ratio: 2.28 [95% CI: 1.19, 4.35]). These results suggest that MT-hgs influence dementia risk and that variants in the nuclear and mitochondrial genome interact to influence the AOO of dementia.

Psychiatric drugs impact mitochondrial function in brain and other tissues

Mitochondria have been linked to the etiology of schizophrenia (SZ). However, studies of mitochondria in SZ might be confounded by the effects of pharmacological treatment with antipsychotic drugs (APDs) and other common medications. This review summarizes findings on relevant mitochondria mechanisms underlying SZ, and the potential impact of psychoactive drugs including primarily APDs, but also antidepressants and anxiolytics. The summarized data suggest that APDs impair mitochondria function by decreasing Complex I activity and ATP production and dissipation of the mitochondria membrane potential. At the same time, in the brains of patients with SZ, antipsychotic drug treatment normalizes gene expression modules enriched in mitochondrial genes that are decreased in SZ. This indicates that APDs may have both positive and negative effects on mitochondria. The available evidence suggests three conclusions i) alterations in mitochondria functions in SZ exist prior to APD treatment, ii) mitochondria alterations in SZ can be reversed by APD treatment, and iii) APDs directly cause impairment of mitochondria function. Overall, the mechanisms of action of psychiatric drugs on mitochondria are both direct and indirect; we conclude the effects of APDs on mitochondria may contribute to both their therapeutic and metabolic side effects. These studies support the hypothesis that neuronal mitochondria are an etiological factor in SZ. Moreover, APDs and other drugs must be considered in the evaluation of this pathophysiological role of mitochondria in SZ. Considering these effects, pharmacological actions on mitochondria may be a worthwhile target for further APD development.

Impact of mitonuclear interactions on life-history responses to diet

Mitochondria are central to both energy metabolism and biosynthesis. Mitochondrial function could therefore influence resource allocation. Critically, mitochondrial function depends on interactions between proteins encoded by the mitochondrial and nuclear genomes. Severe incompatibilities between these genomes can have pervasive effects on both fitness and longevity. How milder deficits in mitochondrial function affect life-history trade-offs is less well understood. Here, we analyse how mitonuclear interactions affect the trade-off between fecundity and longevity in Drosophila melanogaster. We consider a panel of 10 different mitochondrial DNA haplotypes against two contrasting nuclear backgrounds (w¹¹¹⁸ (WE) and Zim53 (ZIM)) in response to high-protein versus standard diet. We report strikingly different responses between the two nuclear backgrounds. WE females have higher fecundity and decreased longevity on high protein. ZIM females have much greater fecundity and shorter lifespan than WE flies on standard diet. High protein doubled their fecundity with no effect on longevity. Mitochondrial haplotype reflected nuclear life-history trade-offs, with a negative correlation between longevity and fecundity in WE flies and no correlation in ZIM flies. Mitonuclear interactions had substantial effects but did not reflect genetic distance between mitochondrial haplotypes. We conclude that mitonuclear interactions can have significant impact on life-history trade-offs, but their effects are not predictable by relatedness. This article is part of the theme issue 'Linking the mitochondrial genotype to phenotype: a complex endeavour'.

Unlocking the origins and biology of domestic animals using ancient DNA and paleogenomics

Animal domestication has fascinated biologists since Charles Darwin first drew the parallel between evolution via natural selection and human-mediated breeding of livestock and companion animals. In this review we show how studies of ancient DNA from domestic animals and their wild progenitors and congeners have shed new light on the genetic origins of domesticates, and on the process of domestication itself. High-resolution paleogenomic data sets now provide unprecedented opportunities to explore the development of animal agriculture across the world. In addition, functional population genomics studies of domestic and wild animals can deliver comparative information useful for understanding recent human evolution.

Mapping and editing animal mitochondrial genomes: can we overcome the challenges?

The animal mitochondrial genome, although small, can have a big impact on health and disease. Non-pathogenic sequence variation among mitochondrial DNA (mtDNA) haplotypes influences traits including fertility, healthspan and lifespan, whereas pathogenic mutations are linked to incurable mitochondrial diseases and other complex conditions like ageing, diabetes, cancer and neurodegeneration. However, we know very little about how mtDNA genetic variation contributes to phenotypic differences. Infrequent recombination, the multicopy nature and nucleic acid-impenetrable membranes present significant challenges that hamper our ability to precisely map mtDNA variants responsible for traits, and to genetically modify mtDNA so that we can isolate specific mutants and characterize their biochemical and physiological consequences. Here, we summarize the past struggles and efforts in developing systems to map and edit mtDNA. We also assess the future of performing forward and reverse genetic studies on animal mitochondrial genomes. This article is part of the theme issue 'Linking the mitochondrial genotype to phenotype: a complex endeavour'.

Mitocellular communication: Shaping health and disease

Throughout the animal kingdom, mitochondria are the only organelles that retain their own genome and the transcription and translation machineries that are all essential for energy harvesting. Mitochondria have developed a complex communication network, allowing them to stay in tune with cellular needs and nuclear transcriptional programs and to alleviate mitochondrial dysfunction. Here, we review recent findings on the wide array of mechanisms that contribute to these mitocellular communication networks, spanning from well-studied messenger molecules to mitonuclear genetic interactions. Based on these observations and developments, we advocate a broad and inclusive view on mitocellular interactions, which can have profound impacts on physiological, pathological, and evolutionary processes.

Mitochondrial Ecophysiology: Assessing the Evolutionary Forces That Shape Mitochondrial Variation

The mitonuclear species concept hypothesizes that incompatibilities between interacting gene products of the nuclear and mitochondrial genomes are a major factor establishing and maintaining species boundaries. However, most of the data available to test this concept come from studies of genetic variation in mitochondrial DNA, and clines in the mitochondrial genome across contact zones can be produced by a variety of forces. Here, we show that using a combination of population genomic analyses of the nuclear and mitochondrial genomes and studies of mitochondrial function can provide insight into the relative roles of neutral processes, adaptive evolution, and mitonuclear incompatibility in establishing and maintaining mitochondrial clines, using Atlantic killifish (Fundulus heteroclitus) as a case study. There is strong evidence for a role of secondary contact following the last glaciation in shaping a steep mitochondrial cline across a contact zone between northern and southern subspecies of killifish, but there is also evidence for a role of adaptive evolution in driving differentiation between the subspecies in a variety of traits from the level of the whole organism to the level of mitochondrial function. In addition, studies are beginning to address the potential for mitonuclear incompatibilities in admixed populations. However, population genomic studies have failed to detect evidence for a strong and pervasive influence of mitonuclear incompatibilities, and we suggest that polygenic selection may be responsible for the complex patterns observed. This case study demonstrates that multiple forces can act together in shaping mitochondrial clines, and illustrates the challenge of disentangling their relative roles.

MOTS-c: A Mitochondrial-Encoded Regulator of the Nucleus

Mitochondria are increasingly being recognized as information hubs that sense cellular changes and transmit messages to other cellular components, such as the nucleus, the endoplasmic reticulum (ER), the Golgi apparatus, and lysosomes. Nonetheless, the interaction between mitochondria and the nucleus is of special interest because they both host part of the cellular genome. Thus, the communication between genome-bearing organelles would likely include gene expression regulation. Multiple nuclear-encoded proteins have been known to regulate mitochondrial gene expression. On the contrary, no mitochondrial-encoded factors are known to actively regulate nuclear gene expression. MOTS-c (mitochondrial open reading frame of the 12S ribosomal RNA type-c) is a recently identified peptide encoded within the mitochondrial 12S ribosomal RNA gene that has metabolic functions. Notably, MOTS-c can translocate to the nucleus upon metabolic stress (e.g., glucose restriction and oxidative stress) and directly regulate adaptive nuclear gene expression to promote cellular homeostasis. It is hypothesized that cellular fitness requires the coevolved mitonuclear genomes to coordinate adaptive responses using gene-encoded factors that cross-regulate the opposite genome. This suggests that cellular gene expression requires the bipartite split genomes to operate as a unified system, rather than the nucleus being the sole master regulator.

CyMIRA: The Cytonuclear Molecular Interactions Reference for Arabidopsis

The function and evolution of eukaryotic cells depend upon direct molecular interactions between gene products encoded in nuclear and cytoplasmic genomes. Understanding how these cytonuclear interactions drive molecular evolution and generate genetic incompatibilities between isolated populations and species is of central importance to eukaryotic biology. Plants are an outstanding system to investigate such effects because of their two different genomic compartments present in the cytoplasm (mitochondria and plastids) and the extensive resources detailing subcellular targeting of nuclear-encoded proteins. However, the field lacks a consistent classification scheme for mitochondrial- and plastid-targeted proteins based on their molecular interactions with cytoplasmic genomes and gene products, which hinders efforts to standardize and compare results across studies. Here, we take advantage of detailed knowledge about the model angiosperm Arabidopsis thaliana to provide a curated database of plant cytonuclear interactions at the molecular level. CyMIRA (Cytonuclear Molecular Interactions Reference for Arabidopsis) is available at http://cymira.colostate.edu/ and https://github.com/dbsloan/cymira and will serve as a resource to aid researchers in partitioning evolutionary genomic data into functional gene classes based on organelle targeting and direct molecular interaction with cytoplasmic genomes and gene products. It includes 11 categories (and 27 subcategories) of different cytonuclear complexes and types of molecular interactions, and it reports residue-level information for cytonuclear contact sites. We hope that this framework will make it easier to standardize, interpret, and compare studies testing the functional and evolutionary consequences of cytonuclear interactions.