The metazoan landscape of mitochondrial DNA gene order and content is shaped by selection and affects mitochondrial transcription.

Mitochondrial DNA (mtDNA) harbors essential genes in most metazoans, yet the regulatory impact of the multiple evolutionary mtDNA rearrangements has been overlooked. Here, by analyzing mtDNAs from ~8000 metazoans we found high gene content conservation (especially of protein and rRNA genes), and codon preferences for mtDNA-encoded tRNAs across most metazoans. In contrast, mtDNA gene order (MGO) was selectively constrained within but not between phyla, yet certain gene stretches (ATP8-ATP6, ND4-ND4L) were highly conserved across metazoans. Since certain metazoans with different MGOs diverge in mtDNA transcription, we hypothesized that evolutionary mtDNA rearrangements affected mtDNA transcriptional patterns. As a first step to test this hypothesis, we analyzed available RNA-seq data from 53 metazoans. Since polycistron mtDNA transcripts constitute a small fraction of the steady-state RNA, we enriched for polycistronic boundaries by calculating RNA-seq read densities across junctions between gene couples encoded either by the same strand (SSJ) or by different strands (DSJ). We found that organisms whose mtDNA is organized in alternating reverse-strand/forward-strand gene blocks (mostly arthropods), displayed significantly reduced DSJ read counts, in contrast to organisms whose mtDNA genes are preferentially encoded by one strand (all chordates). Our findings suggest that mtDNA rearrangements are selectively constrained and likely impact mtDNA regulation.

Human mitochondrial RNA modifications associate with tissue-specific changes in gene expression, and are affected by sunlight and UV exposure.

RNA-DNA differences (RDD) have previously been identified in the human mitochondrial RNA (mt-RNA) transcripts, yet their functional impact is poorly understood. By analyzing 4928 RNA-seq samples from 23 body sites, we found that mtDNA gene expression negatively correlated with the levels of both m1A 947 16 S rRNA modification (mtDNA position 2617) and the m1A 1812 ND5 mRNA modification (mtDNA position 13,710) in 15 and 14 body sites, respectively. Such correlation was not evident in all tested brain tissues, thus suggesting a tissue-specific impact of these modifications on mtDNA gene expression. To assess the response of the tested modifications to environmental cues, we analyzed pairs of skin samples that were either exposed to the sun or not. We found that the correlations of mtDNA gene expression with both mt-RNA modifications were compromised upon sun exposure. As a first step to explore the underlying mechanism, we analyzed RNA-seq data from keratinocytes that were exposed to increasing doses of UV irradiation. Similar to sun exposure, we found a significant decrease in mtDNA gene expression upon increase in UV dosage. In contrast, there was a significant increase in the m1A 947 16 S rRNA modification levels upon elevation in UV dose. Finally, we identified candidate modulators of such responses. Taken together, our results indicate that mt-RNA modifications functionally correlate with mtDNA gene expression, and responds to environmental cues, hence supporting their physiological importance.

Mutant C. elegans mitofusin leads to selective removal of mtDNA heteroplasmic deletions across generations to maintain fitness.

BACKGROUND: Mitochondrial DNA (mtDNA) is present at high copy numbers in animal cells, and though characterized by a single haplotype in each individual due to maternal germline inheritance, deleterious mutations and intact mtDNA molecules frequently co-exist (heteroplasmy). A number of factors, such as replicative segregation, mitochondrial bottlenecks, and selection, may modulate the exitance of heteroplasmic mutations. Since such mutations may have pathological consequences, they likely survive and are inherited due to functional complementation via the intracellular mitochondrial network. Here, we hypothesized that compromised mitochondrial fusion would hamper such complementation, thereby affecting heteroplasmy inheritance. RESULTS: We assessed heteroplasmy levels in three Caenorhabditis elegans strains carrying different heteroplasmic mtDNA deletions (ΔmtDNA) in the background of mutant mitofusin (fzo-1). Animals displayed severe embryonic lethality and developmental delay. Strikingly, observed phenotypes were relieved during subsequent generations in association with complete loss of ΔmtDNA molecules. Moreover, deletion loss rates were negatively correlated with the size of mtDNA deletions, suggesting that mitochondrial fusion is essential and sensitive to the nature of the heteroplasmic mtDNA mutations. Introducing the ΔmtDNA into a fzo-1;pdr-1;+/ΔmtDNA (PARKIN ortholog) double mutant resulted in a skewed Mendelian progeny distribution, in contrast to the normal distribution in the fzo-1;+/ΔmtDNA mutant, and severely reduced brood size. Notably, the ΔmtDNA was lost across generations in association with improved phenotypes. CONCLUSIONS: Taken together, our findings show that when mitochondrial fusion is compromised, deleterious heteroplasmic mutations cannot evade natural selection while inherited through generations. Moreover, our findings underline the importance of cross-talk between mitochondrial fusion and mitophagy in modulating the inheritance of mtDNA heteroplasmy.

Immune system cells from COVID-19 patients display compromised mitochondrial-nuclear expression co-regulation and rewiring toward glycolysis.

Mitochondria are pivotal for bioenergetics, as well as in cellular response to viral infections. Nevertheless, their role in COVID-19 was largely overlooked. Here, we analyzed available bulk RNA-seq datasets from COVID-19 patients and corresponding healthy controls (three blood datasets, N = 48 healthy, 119 patients; two respiratory tract datasets, N = 157 healthy, 524 patients). We found significantly reduced mtDNA gene expression in blood, but not in respiratory tract samples from patients. Next, analysis of eight single-cells RNA-seq datasets from peripheral blood mononuclear cells, nasopharyngeal samples, and Bronchoalveolar lavage fluid (N = 1,192,243 cells), revealed significantly reduced mtDNA gene expression especially in immune system cells from patients. This is associated with elevated expression of nuclear DNA-encoded OXPHOS subunits, suggesting compromised mitochondrial-nuclear co-regulation. This, together with elevated expression of ROS-response genes and glycolysis enzymes in patients, suggest rewiring toward glycolysis, thus generating beneficial conditions for SARS-CoV-2 replication. Our findings underline the centrality of mitochondrial dysfunction in COVID-19.

Two Homogametic Genotypes - One Crayfish: On the Consequences of Intersexuality.

In the Australian redclaw crayfish, Cherax quadricarinatus (WZ/ZZ system), intersexuals, although exhibiting both male and female gonopores, are functional males bearing a female genotype (WZ males). Therefore, the occurrence of the unusual homogametic WW females in nature is plausible. We developed W/Z genomic sex markers and used them to investigate the genotypic structure of experimental and native C. quadricarinatus populations in Australia. We discovered, for the first time, the natural occurrence of WW females in crustacean populations. By modeling population dynamics, we found that intersexuals contribute to the growth rate of crayfish populations in the short term. Given the vastly fragmented C. quadricarinatus habitat, which is characterized by drought-flood cycles, we speculate that intersexuals contribute to the fitness of this species since they lead to occasional increment in the population growth rate which potentially supports crayfish population restoration and establishment under extinction threats or colonization events.

Mitochondria Are Fundamental for the Emergence of Metazoans: On Metabolism, Genomic Regulation, and the Birth of Complex Organisms.

Out of many intracellular bacteria, only the mitochondria and chloroplasts abandoned their independence billions of years ago and became endosymbionts within the host eukaryotic cell. Consequently, one cannot grow eukaryotic cells without their mitochondria, and the mitochondria cannot divide outside of the cell, thus reflecting interdependence. Here, we argue that such interdependence underlies the fundamental role of mitochondrial activities in the emergence of metazoans. Several lines of evidence support our hypothesis: (a) Differentiation and embryogenesis rely on mitochondrial function; (b) mitochondrial metabolites are primary precursors for epigenetic modifications (such as methyl and acetyl), which are critical for chromatin remodeling and gene expression, particularly during differentiation and embryogenesis; and (c) mitonuclear coregulation adapted to accommodate both housekeeping and tissue-dependent metabolic needs. We discuss the evolution of the unique mitochondrial genetic system, mitochondrial metabolites, mitonuclear coregulation, and their critical roles in the emergence of metazoans and in human disorders.

Disease-causing mutations in subunits of OXPHOS complex I affect certain physical interactions.

Mitochondrial complex I (CI) is the largest multi-subunit oxidative phosphorylation (OXPHOS) protein complex. Recent availability of a high-resolution human CI structure, and from two non-human mammals, enabled predicting the impact of mutations on interactions involving each of the 44 CI subunits. However, experimentally assessing the impact of the predicted interactions requires an easy and high-throughput method. Here, we created such a platform by cloning all 37 nuclear DNA (nDNA) and 7 mitochondrial DNA (mtDNA)-encoded human CI subunits into yeast expression vectors to serve as both 'prey' and 'bait' in the split murine dihydrofolate reductase (mDHFR) protein complementation assay (PCA). We first demonstrated the capacity of this approach and then used it to examine reported pathological OXPHOS CI mutations that occur at subunit interaction interfaces. Our results indicate that a pathological frame-shift mutation in the MT-ND2 gene, causing the replacement of 126 C-terminal residues by a stretch of only 30 amino acids, resulted in loss of specificity in ND2-based interactions involving these residues. Hence, the split mDHFR PCA is a powerful assay for assessing the impact of disease-causing mutations on pairwise protein-protein interactions in the context of a large protein complex, thus offering a possible mechanistic explanation for the underlying pathogenicity.

mtDNA Chromatin-like Organization Is Gradually Established during Mammalian Embryogenesis.

Unlike the nuclear genome, the mammalian mitochondrial genome (mtDNA) is thought to be coated solely by mitochondrial transcription factor A (TFAM), whose binding sequence preferences are debated. Therefore, higher-order mtDNA organization is considered much less regulated than both the bacterial nucleoid and the nuclear chromatin. However, our recently identified conserved DNase footprinting pattern in human mtDNA, which co-localizes with regulatory elements and responds to physiological conditions, likely reflects a structured higher-order mtDNA organization. We hypothesized that this pattern emerges during embryogenesis. To test this hypothesis, we analyzed assay for transposase-accessible chromatin sequencing (ATAC-seq) results collected during the course of mouse and human early embryogenesis. Our results reveal, for the first time, a gradual and dynamic emergence of the adult mtDNA footprinting pattern during embryogenesis of both mammals. Taken together, our findings suggest that the structured adult chromatin-like mtDNA organization is gradually formed during mammalian embryogenesis.

Higher Order Organization of the mtDNA: Beyond Mitochondrial Transcription Factor A.

The higher order organization of eukaryotic and prokaryotic genomes is pivotal in the regulation of gene expression. Specifically, chromatin accessibility in eukaryotes and nucleoid accessibility in bacteria are regulated by a cohort of proteins to alter gene expression in response to diverse physiological conditions. By contrast, prior studies have suggested that the mitochondrial genome (mtDNA) is coated solely by mitochondrial transcription factor A (TFAM), whose increased cellular concentration was proposed to be the major determinant of mtDNA packaging in the mitochondrial nucleoid. Nevertheless, recent analysis of DNase-seq and ATAC-seq experiments from multiple human and mouse samples suggest gradual increase in mtDNA occupancy during the course of embryonic development to generate a conserved footprinting pattern which correlate with sites that have low TFAM occupancy in vivo (ChIP-seq) and tend to adopt G-quadruplex structures. These findings, along with recent identification of mtDNA binding by known modulators of chromatin accessibility such as MOF, suggest that mtDNA higher order organization is generated by cross talk with the nuclear regulatory system, may have a role in mtDNA regulation, and is more complex than once thought.

A common pattern of DNase I footprinting throughout the human mtDNA unveils clues for a chromatin-like organization.

Human mitochondrial DNA (mtDNA) is believed to lack chromatin and histones. Instead, it is coated solely by the transcription factor TFAM. We asked whether mtDNA packaging is more regulated than once thought. To address this, we analyzed DNase-seq experiments in 324 human cell types and found, for the first time, a pattern of 29 mtDNA Genomic footprinting (mt-DGF) sites shared by ∼90% of the samples. Their syntenic conservation in mouse DNase-seq experiments reflect selective constraints. Colocalization with known mtDNA regulatory elements, with G-quadruplex structures, in TFAM-poor sites (in HeLa cells) and with transcription pausing sites, suggest a functional regulatory role for such mt-DGFs. Altered mt-DGF pattern in interleukin 3-treated CD34+ cells, certain tissue differences, and significant prevalence change in fetal versus nonfetal samples, offer first clues to their physiological importance. Taken together, human mtDNA has a conserved protein-DNA organization, which is likely involved in mtDNA regulation.

Mitochondrial DNA associations with East Asian metabolic syndrome.

Mitochondrial dysfunction has repeatedly been reported associated with type 2 diabetes mellitus (T2DM) and metabolic syndrome (MS), as have mitochondrial DNA (mtDNA) tRNA and duplication mutations and mtDNA haplogroup lineages. We identified 19 Taiwanese T2DM and MS pedigrees from Taiwan, with putative matrilineal transmission, one of which harbored the pathogenic mtDNA tRNALeu(UUR) nucleotide (nt) 3243A>G mutation on the N9a3 haplogroup background. We then recruited three independent Taiwanese cohorts, two from Taipei (N = 498, mean age 52 and N = 1002, mean age 44) and one from a non-urban environment (N = 501, mean age 57). All three cohorts were assessed for an array of metabolic parameters, their mtDNA haplogroups determined, and the haplogroups correlated with T2DM/MS phenotypes. Logistic regression analysis revealed that mtDNA haplogroups D5, F4, and N9a conferred T2DM protection, while haplogroups F4 and N9a were risk factors for hypertension (HTN), and F4 was a risk factor for obesity (OB). Additionally, the 5263C>T (ND2 A165V) variant commonly associated with F4 was associated with hypertension (HTN). Cybrids were prepared with macro-haplogroup N (defined by variants m.ND3 10398A (114T) and m.ATP6 8701A (59T)) haplogroups B4 and F1 mtDNAs and from macro-haplogroup M (variants m.ND3 10398G (114A) and m.ATP6 8701G (59A)) haplogroup M9 mtDNAs. Additionally, haplogroup B4 and F1 cybrids were prepared with and without the mtDNA variant in ND1 3394T>C (Y30H) reported to be associated with T2DM. Assay of mitochondria complex I in these cybrids revealed that macro-haplogroup N cybrids had lower activity than M cybrids, that haplogroup F cybrids had lower activity than B4 cybrids, and that the ND1 3394T>C (Y30H) variant reduced complex I on both the B4 and F1 background but with very different cumulative effects. These data support the hypothesis that functional mtDNA variants may contribute to the risk of developing T2DM and MS.

Mitochondrial DNA Transcription and Its Regulation: An Evolutionary Perspective.

The bacterial heritage of mitochondria, as well as its independent genome [mitochondrial DNA (mtDNA)] and polycistronic transcripts, led to the view that mitochondrial transcriptional regulation relies on an evolutionarily conserved, prokaryotic-like system that is separated from the rest of the cell. Indeed, mtDNA transcription was previously thought to be governed by a few dedicated direct regulators, namely, the mitochondrial RNA polymerase (POLRMT), two transcription factors (TFAM and TF2BM), one transcription elongation (TEFM), and one known transcription termination factor (mTERF1). Recent findings have, however, revealed that known nuclear gene expression regulators are also involved in mtDNA transcription and have identified novel transcriptional features consistent with adaptation of the mitochondria to the regulatory environment of the precursor of the eukaryotic cell. Finally, whereas mammals follow the human mtDNA transcription pattern, other organisms notably diverge in terms of mtDNA transcriptional regulation. Hence, mtDNA transcriptional regulation is likely more evolutionary diverse than once thought.

Human primitive brain displays negative mitochondrial-nuclear expression correlation of respiratory genes.

Oxidative phosphorylation (OXPHOS), a fundamental energy source in all human tissues, requires interactions between mitochondrial (mtDNA)- and nuclear (nDNA)-encoded protein subunits. Although such interactions are fundamental to OXPHOS, bi-genomic coregulation is poorly understood. To address this question, we analyzed ∼8500 RNA-seq experiments from 48 human body sites. Despite well-known variation in mitochondrial activity, quantity, and morphology, we found overall positive mtDNA-nDNA OXPHOS genes' co-expression across human tissues. Nevertheless, negative mtDNA-nDNA gene expression correlation was identified in the hypothalamus, basal ganglia, and amygdala (subcortical brain regions, collectively termed the "primitive" brain). Single-cell RNA-seq analysis of mouse and human brains revealed that this phenomenon is evolutionarily conserved, and both are influenced by brain cell types (involving excitatory/inhibitory neurons and nonneuronal cells) and by their spatial brain location. As the "primitive" brain is highly oxidative, we hypothesized that such negative mtDNA-nDNA co-expression likely controls for the high mtDNA transcript levels, which enforce tight OXPHOS regulation, rather than rewiring toward glycolysis. Accordingly, we found "primitive" brain-specific up-regulation of lactate dehydrogenase B (LDHB), which associates with high OXPHOS activity, at the expense of LDHA, which promotes glycolysis. Analyses of co-expression, DNase-seq, and ChIP-seq experiments revealed candidate RNA-binding proteins and CEBPB as the best regulatory candidates to explain these phenomena. Finally, cross-tissue expression analysis unearthed tissue-dependent splice variants and OXPHOS subunit paralogs and allowed revising the list of canonical OXPHOS transcripts. Taken together, our analysis provides a comprehensive view of mito-nuclear gene co-expression across human tissues and provides overall insights into the bi-genomic regulation of mitochondrial activities.

The First Mitochondrial Genomics and Evolution SMBE-Satellite Meeting: A New Scientific Symbiosis.

The central role of the mitochondrion for cellular and organismal metabolism is well known, yet its functional role in evolution has rarely been featured in leading international conferences. Moreover, the contribution of mitochondrial genetics to complex disease phenotypes is particularly important, and although major advances have been made in the field of genomics, mitochondrial genomic data have in many cases been overlooked. Accumulating data and new knowledge support a major contribution of this maternally inherited genome, and its interactions with the nucleus, to both major evolutionary processes and diverse disease phenotypes. These advances encouraged us to assemble the first Mitochondrial Genomics and Evolution (MGE) meeting-an SMBE satellite and Israeli Science foundation international conference (Israel, September 2017). Here, we report the content and outcome of the MGE meeting (https://www.mge2017.com/; last accessed November 5, 2017).

The genomic landscape of evolutionary convergence in mammals, birds and reptiles.

Many lineage-defining (nodal) mutations possess high functionality. However, differentiating adaptive nodal mutations from those that are functionally compensated remains challenging. To address this challenge, we identified functional nodal mutations (fNMs) in ~3,400 nuclear DNA (nDNA) and 4 mitochondrial DNA (mtDNA) protein structures from 91 and 1,003 species, respectively, representing the entire mammalian, bird and reptile phylogeny. A screen for candidate compensatory mutations among co-occurring amino acid changes in close structural proximity revealed that such compensated fNMs encompass 37% and 27% of the mtDNA and nDNA datasets, respectively. Analysis of the remaining (non-compensated) mutations, which are enriched for adaptive mutations, showed that birds and mammals share most such recurrent fNMs (N = 51). Among the latter, we discovered mutations in thermoregulation-related genes. These represent the best candidates to explain the molecular basis of convergent body thermoregulation in birds and mammals. Our analysis reveals the landscape of possible mutational compensation and convergence in amniote phylogeny.