From evolutionary bystander to master manipulator: the emerging roles for the mitochondrial genome as a modulator of nuclear gene expression

Evolutionary genetics of the mitochondrial genome: insights from Drosophila

Mitochondria are key to energy conversion in virtually all eukaryotes. Intriguingly, despite billions of years of evolution inside the eukaryote, mitochondria have retained their own small set of genes involved in the regulation of oxidative phosphorylation (OXPHOS) and protein translation. Although there was a long-standing assumption that the genetic variation found within the mitochondria would be selectively neutral, research over the past 3 decades has challenged this assumption. This research has provided novel insight into the genetic and evolutionary forces that shape mitochondrial evolution and broader implications for evolutionary ecological processes. Many of the seminal studies in this field, from the inception of the research field to current studies, have been conducted using Drosophila flies, thus establishing the species as a model system for studies in mitochondrial evolutionary biology. In this review, we comprehensively review these studies, from those focusing on genetic processes shaping evolution within the mitochondrial genome, to those examining the evolutionary implications of interactions between genes spanning mitochondrial and nuclear genomes, and to those investigating the dynamics of mitochondrial heteroplasmy. We synthesize the contribution of these studies to shaping our understanding of the evolutionary and ecological implications of mitochondrial genetic variation.

Pathogenic mitochondrial DNA 3243A>G mutation: From genetics to phenotype

The mitochondrial DNA (mtDNA) m.3243A>G mutation is one of the most common pathogenic mtDNA variants, showing complex genetics, pathogenic molecular mechanisms, and phenotypes. In recent years, the prevention of mtDNA-related diseases has trended toward precision medicine strategies, such as preimplantation genetic diagnosis (PGD) and mitochondrial replacement therapy (MRT). These techniques are set to allow the birth of healthy children, but clinical implementation relies on thorough insights into mtDNA genetics. The genotype and phenotype of m.3243A>G vary greatly from mother to offspring, which compromises genetic counseling for the disease. This review is the first to systematically elaborate on the characteristics of the m.3243A>G mutation, from genetics to phenotype and the relationship between them, as well as the related influencing factors and potential strategies for preventing disease. These perceptions will provide clarity for clinicians providing genetic counseling to m.3243A>G patients.

Interspecific Gene Flow and Mitochondrial Genome Capture During the Radiation of Jamaican Anolis Lizards (Squamata; Iguanidae)

Gene flow and reticulation are increasingly recognized as important processes in the diversification of many taxonomic groups. With the increasing ease of collecting genomic data and the development of multispecies coalescent network approaches, such reticulations can be accounted for when inferring phylogeny and diversification. Caribbean Anolis lizards are a classic example of an adaptive radiation in which species have independently radiated on the islands of the Greater Antilles into the same ecomorph classes. Within the Jamaican radiation at least one species, A. opalinus, has been documented to be polyphyletic in its mitochondrial DNA, which could be the result of an ancient reticulation event or incomplete lineage sorting. Here we generate mtDNA and genotyping-by-sequencing (GBS) data and implement gene-tree, species-tree, and multispecies coalescent network methods to infer the diversification of this group. Our mtDNA gene-tree recovers the same relationships previously inferred for this group, which is strikingly different from the species-tree inferred from our GBS data. Posterior predictive simulations suggest that our genomic data violate commonly adopted assumptions of the multispecies coalescent model, so we use network approaches to infer phylogenetic relationships. The inferred network topology contains a reticulation event but does not explain the mtDNA polyphyly observed in this group, however coalescent simulations suggest that the observed mtDNA topology is likely the result of past introgression. How common a signature of gene flow and reticulation is across the radiation of Anolis is unknown; however, the reticulation events that we demonstrate here may have allowed for adaptive evolution, as has been suggested in other, more recent adaptive radiations.

Interactions between cytoplasmic and nuclear genomes confer sex-specific effects on lifespan in Drosophila melanogaster

Genetic variation outside of the cell nucleus can affect the phenotype. The cytoplasm is home to the mitochondria, and in arthropods often hosts intracellular bacteria such as Wolbachia. Although numerous studies have implicated epistatic interactions between cytoplasmic and nuclear genetic variation as mediators of phenotypic expression, two questions remain. Firstly, it remains unclear whether outcomes of cyto-nuclear interactions will manifest differently across the sexes, as might be predicted given that cytoplasmic genomes are screened by natural selection only through females as a consequence of their maternal inheritance. Secondly, the relative contribution of mitochondrial genetic variation to other cytoplasmic sources of variation, such as Wolbachia infection, in shaping phenotypic outcomes of cyto-nuclear interactions remains unknown. Here, we address these questions, creating a fully crossed set of replicated cyto-nuclear populations derived from three geographically distinct populations of Drosophila melanogaster, measuring the lifespan of males and females from each population. We observed that cyto-nuclear interactions shape lifespan and that the outcomes of these interactions differ across the sexes. Yet, we found no evidence that placing the cytoplasms from one population alongside the nuclear background of others (generating putative cyto-nuclear mismatches) leads to decreased lifespan in either sex. Although it was difficult to partition mitochondrial from Wolbachia effects, our results suggest at least some of the cytoplasmic genotypic contribution to lifespan was directly mediated by an effect of sequence variation in the mtDNA. Future work should explore the degree to which cyto-nuclear interactions result in sex differences in the expression of other components of organismal life history.

Resolving complexity in mitochondrial disease: Towards precision medicine

Mitochondrial diseases, caused by mutations in either the nuclear or mitochondrial genomes (mtDNA), are the most common form of inherited neurometabolic disorders. They are remarkably heterogeneous, both in their clinical presentation and genetic etiology, presenting challenges for diagnosis, clinical management and elucidation of molecular mechanism. The multifaceted nature of these diseases, compounded by the unique characteristics of mitochondrial genetics, cement their space in the field of complex disease. In this review we examine the m.3243A>G variant, one of the most prevalent mitochondrial DNA mutations, using it as an exemplar to demonstrate the challenges presented by these complex disorders. Disease caused by m.3243A>G is one of the most phenotypically diverse of all mitochondrial diseases; we outline known causes of this heterogeneity including mtDNA heteroplasmy, mtDNA copy number and nuclear genetic factors. We consider the impact that this has in the clinic, discussing the personalized management of common manifestations attributed to this pathogenic mtDNA variant, including hearing impairment, diabetes mellitus, myopathy, cardiac disease, stroke-like episodes and gastrointestinal disturbances. Future research into this complex disorder must account for this heterogeneity, benefitting from the use of large patient cohorts to build upon current clinical expertise. Through multi-disciplinary collaboration, the complexities of this mitochondrial disease can be addressed with the variety of diagnostic, prognostic, and treatment approaches that are moulded to best fit the needs of each individual patient.

Age Dependent Dysfunction of Mitochondrial and ROS Metabolism Induced by Mitonuclear Mismatch

Mitochondrial and nuclear genomes have to coevolve to ensure the proper functioning of the different mitochondrial complexes that are assembled from peptides encoded by both genomes. Mismatch between these genomes is believed to be strongly selected against due to the consequent impairments of mitochondrial functions and induction of oxidative stress. Here, we used a Drosophila model harboring an incompatibility between a mitochondrial tRNA^tyr and its nuclear-encoded mitochondrial tyrosine synthetase to assess the cellular mechanisms affected by this incompatibility and to test the relative contribution of mitonuclear interactions and aging on the expression of impaired phenotypes. Our results show that the mitochondrial tRNA mutation caused a decrease in mitochondrial oxygen consumption in the incompatible nuclear background but no effect with the compatible nuclear background. Mitochondrial DNA copy number increased in the incompatible genotype but that increase failed to rescue mitochondrial functions. The flies harboring mismatch between nuclear and mitochondrial genomes had almost three times the relative mtDNA copy number and fifty percent higher rate of hydrogen peroxide production compared to other genome combinations at 25 days of age. We also found that aging exacerbated the mitochondrial dysfunctions. Our results reveal the tight interactions linking mitonuclear mismatch to mitochondrial dysfunction, mitochondrial DNA regulation, ROS production and aging.

Evolved genetic and phenotypic differences due to mitochondrial-nuclear interactions

The oxidative phosphorylation (OxPhos) pathway is responsible for most aerobic ATP production and is the only pathway with both nuclear and mitochondrial encoded proteins. The importance of the interactions between these two genomes has recently received more attention because of their potential evolutionary effects and how they may affect human health and disease. In many different organisms, healthy nuclear and mitochondrial genome hybrids between species or among distant populations within a species affect fitness and OxPhos functions. However, what is less understood is whether these interactions impact individuals within a single natural population. The significance of this impact depends on the strength of selection for mito-nuclear interactions. We examined whether mito-nuclear interactions alter allele frequencies for ~11,000 nuclear SNPs within a single, natural Fundulus heteroclitus population containing two divergent mitochondrial haplotypes (mt-haplotypes). Between the two mt-haplotypes, there are significant nuclear allele frequency differences for 349 SNPs with a p-value of 1% (236 with 10% FDR). Unlike the rest of the genome, these 349 outlier SNPs form two groups associated with each mt-haplotype, with a minority of individuals having mixed ancestry. We use this mixed ancestry in combination with mt-haplotype as a polygenic factor to explain a significant fraction of the individual OxPhos variation. These data suggest that mito-nuclear interactions affect cardiac OxPhos function. The 349 outlier SNPs occur in genes involved in regulating metabolic processes but are not directly associated with the 79 nuclear OxPhos proteins. Therefore, we postulate that the evolution of mito-nuclear interactions affects OxPhos function by acting upstream of OxPhos.

MOTS-c: A novel mitochondrial-derived peptide regulating muscle and fat metabolism

Mitochondria are ancient organelles that are thought to have emerged from once free-living α-proto-bacteria. As such, they still possess several bacterial-like qualities, including a semi-autonomous genetic system, complete with an independent genome and a unique genetic code. The bacterial-like circular mitochondrial DNA (mtDNA) has been described to encode 37 genes, including 22 tRNAs, 2 rRNAs, and 13 mRNAs. Two additional peptides reported to originate from the mtDNA, namely humanin (Hashimoto et al., 2001; Ikone et al., 2003; Guo et al., 2003) [1-3] and MOTS-c (mitochondrial ORF of the twelve S c) (Lee et al., 2015) [4], indicate a larger mitochondrial genetic repertoire (Shokolenko and Alexeyev, 2015) [5]. These mitochondrial-derived peptides (MDPs) have profound and distinct biological activities and provide a paradigm-shifting concept of active mitochondrial-encoded signals that act at the cellular and organismal level (i.e. mitochondrial hormone) (da Cunha et al., 2015; Quiros et al., 2016) [6,7]. Considering that mitochondria are the single most important metabolic organelle, it is not surprising that these MDPs have metabolic actions. MOTS-c has been shown to target the skeletal muscle and enhance glucose metabolism. As such, MOTS-c has implications in the regulation of obesity, diabetes, exercise, and longevity, representing an entirely novel mitochondrial signaling mechanism to regulate metabolism within and between cells.

Spotlight on the relevance of mtDNA in cancer

The potential role of the mitochondrial genome has recently attracted interest because of its high mutation frequency in tumors. Different aspects of mtDNA make it relevant for cancer's biology, such as it encodes a limited but essential number of genes for OXPHOS biogenesis, it is particularly susceptible to mutations, and its copy number can vary. Moreover, most ROS in mitochondria are produced by the electron transport chain. These characteristics place the mtDNA in the center of multiple signaling pathways, known as mitochondrial retrograde signaling, which modifies numerous key processes in cancer. Cybrid studies support that mtDNA mutations are relevant and exert their effect through a modification of OXPHOS function and ROS production. However, there is still much controversy regarding the clinical relevance of mtDNA mutations. New studies should focus more on OXPHOS dysfunction associated with a specific mutational signature rather than the presence of mutations in the mtDNA.

Complete mitochondrial genome sequences of thirteen globally sourced strains of fruit fly (Drosophila melanogaster) form a powerful model for mitochondrial research

The complete mitogenomes of 13 strains of the fruit fly Drosophila melanogaster were sequenced. Haplotypes varied between 19 532 and 19 537 bp in length, and followed standard dipteran mitogenome content and organization. We detected a total of 354 variable sites between all thirteen haplotypes, while single pairs of haplotypes were separated by an average of 123 variable sites. The sequenced fly strains form a powerful model for mitochondrial research, when it comes to elucidating the links between the mitochondrial genotype and the phenotype.

Mitochondrial genome diversity among six laboratory zebrafish (Danio rerio) strains

The mitochondrial genome (mtDNA) is distinct from the nuclear genome and is known to play a significant role in several disease phenotypes, such as longevity and fertility. Here we characterize the complete mitochondrial genomes (∼16 590 bp), and the extent of within and between strain variation for 27 adult zebrafish, representing five commonly used laboratory strains (AB, TL, HL, WIK, and SJD) and one line acquired from a local pet shop. These data were subsequently analyzed to determine the phylogenetic relationships between strains and ascertain if positive selection might be operating on any mtDNA genes. Relationships between strains are not entirely consistent with those observed previously using nuclear DNA. Further there is a substantial body of variation within current zebrafish lines, with 172 variants described across lines. Of these, 27 changes are non-synonymous and there is nominal evidence for positive selection in the mtDNA sequences at some of these sites. We further identify novel frameshift mutations in eight genes, which are all predicted to have functional consequences. Our study provides the first information on mtDNA diversity in zebrafish, identifies multiple non-synonymous substitutions and other mutations expected to have functional effects, and represents an important first step in establishing zebrafish as a model for investigating the phenotypic effects of mtDNA mutations.

Mitonuclear Ecology

Eukaryotes were born of a chimeric union between two prokaryotes--the progenitors of the mitochondrial and nuclear genomes. Early in eukaryote evolution, most mitochondrial genes were lost or transferred to the nucleus, but a core set of genes that code exclusively for products associated with the electron transport system remained in the mitochondrion. The products of these mitochondrial genes work in intimate association with the products of nuclear genes to enable oxidative phosphorylation and core energy production. The need for coadaptation, the challenge of cotransmission, and the possibility of genomic conflict between mitochondrial and nuclear genes have profound consequences for the ecology and evolution of eukaryotic life. An emerging interdisciplinary field that I call "mitonuclear ecology" is reassessing core concepts in evolutionary ecology including sexual reproduction, two sexes, sexual selection, adaptation, and speciation in light of the interactions of mitochondrial and nuclear genomes.

Mitochondrial DNA variation and virologic and immunological HIV outcomes in African Americans

OBJECTIVE

To evaluate the impact of mitochondrial DNA (mtDNA) haplogroups on virologic and immunological outcomes of HIV infection.

DESIGN

HAART-naive African American adolescent participants to the Reaching for Excellence in Adolescent Care and Health study.

METHODS

The mtDNA haplogroups were inferred from sequenced mtDNA hypervariable regions HV1 and HV2 and their predictive value on HIV outcomes were evaluated in linear mixed models, controlled for human leukocyte antigen (HLA)-B27, HLA-B57 and HLA-B35-Px alleles and other covariates.

RESULTS

We report data showing that the mtDNA L2 lineage, a group composed of L2a, L2b and L2e mtDNA haplogroups in the studied population, is significantly associated (beta  = -0.08; Bonferroni-adjusted P = 0.004) with decline of CD4 T cells (median loss of 8 ± 1 cells per month) in HAART-naive HIV-infected individuals of African American descent (n = 133). No significant association (P < 0.05) with set-point viral load was observed with any of the tested mtDNA haplogroups. The present data concur with previous findings in the AIDS Clinical Trials Group study 384, implicating the L2 lineage with slower CD4 T-cell recovery after antiretroviral therapy in African Americans.

CONCLUSIONS

Whereas the L2 lineage showed an association with unfavorable immunological outcomes of HIV infection, its phylogenetic divergence from J and U5a, two lineages associated with accelerated HIV progression in European Americans, raises the possibility that interactions with common nucleus-encoded variants drive HIV progression. Disentangling the effects of mitochondrial and nuclear gene variants on the outcomes of HIV infection is an important step to be taken toward a better understanding of HIV/AIDS pathogenesis and pharmacogenomics.