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.

Functional analysis of ADARs in planarians supports a bilaterian ancestral role in suppressing double-stranded RNA-response.

ADARs (adenosine deaminases acting on RNA) are known for their adenosine-to-inosine RNA editing activity, and most recently, for their role in preventing aberrant dsRNA-response by activation of dsRNA sensors (i.e., RIG-I-like receptor homologs). However, it is still unclear whether suppressing spurious dsRNA-response represents the ancestral role of ADARs in bilaterians. As a first step to address this question, we identified ADAR1 and ADAR2 homologs in the planarian Schmidtea mediterranea, which is evolutionarily distant from canonical lab models (e.g., flies and nematodes). Our results indicate that knockdown of either planarian adar1 or adar2 by RNA interference (RNAi) resulted in upregulation of dsRNA-response genes, including three planarian rig-I-like receptor (prlr) homologs. Furthermore, independent knockdown of adar1 and adar2 reduced the number of infected cells with a dsRNA virus, suggesting they suppress a bona fide anti-viral dsRNA-response activity. Knockdown of adar1 also resulted in lesion formation and animal lethality, thus attesting to its essentiality. Simultaneous knockdown of adar1 and prlr1 rescued adar1(RNAi)-dependent animal lethality and rescued the dsRNA-response, suggesting that it contributes to the deleterious effect of adar1 knockdown. Finally, we found that ADAR2, but not ADAR1, mediates mRNA editing in planarians, suggesting at least in part non-redundant activities for planarians ADARs. Our results underline the essential role of ADARs in suppressing activation of harmful dsRNA-response in planarians, thus supporting it as their ancestral role in bilaterians. Our work also set the stage to study further and better understand the regulatory mechanisms governing anti-viral dsRNA-responses from an evolutionary standpoint using planarians as a model.

RNA editing in bacteria: occurrence, regulation and significance.

DNA harbors the blueprint for life. However, the instructions stored in the DNA could be altered at the RNA level before they are executed. One of these processes is RNA editing, which was shown to modify RNA sequences in many organisms. The most abundant modification is the deamination of adenosine (A) into inosine (I). In turn, inosine can be identified as a guanosine (G) by the ribosome and other cellular machineries such as reverse transcriptase. In multicellular organisms, enzymes from the ADAR (adenosine deaminase acting on RNA) family mediate RNA editing in mRNA, whereas enzymes from the ADAT family mediate A-to-I editing on tRNAs. In bacteria however, until recently, only one editing site was described, in tRNAArg, but never in mRNA. The tRNA site was shown to be modified by tadA (tRNA specific adenosine deaminase) which is believed to be the ancestral enzyme for the RNA editing family of enzymes. In our recent work, we have shown for the first time, editing on multiple sites in bacterial mRNAs and identified tadA as the enzyme responsible for this editing activity. Focusing on one of the identified targets - the self-killing toxin hokB, we found that editing is physiologically regulated and that it increases protein activity. Here we discuss possible modes of regulation on hokB editing, potential roles of RNA editing in bacteria, possible implications, and future research directions.

The m1A landscape on cytosolic and mitochondrial mRNA at single-base resolution.

Modifications on mRNA offer the potential of regulating mRNA fate post-transcriptionally. Recent studies suggested the widespread presence of N1-methyladenosine (m1A), which disrupts Watson-Crick base pairing, at internal sites of mRNAs. These studies lacked the resolution of identifying individual modified bases, and did not identify specific sequence motifs undergoing the modification or an enzymatic machinery catalysing them, rendering it challenging to validate and functionally characterize putative sites. Here we develop an approach that allows the transcriptome-wide mapping of m1A at single-nucleotide resolution. Within the cytosol, m1A is present in a low number of mRNAs, typically at low stoichiometries, and almost invariably in tRNA T-loop-like structures, where it is introduced by the TRMT6/TRMT61A complex. We identify a single m1A site in the mitochondrial ND5 mRNA, catalysed by TRMT10C, with methylation levels that are highly tissue specific and tightly developmentally controlled. m1A leads to translational repression, probably through a mechanism involving ribosomal scanning or translation. Our findings suggest that m1A on mRNA, probably because of its disruptive impact on base pairing, leads to translational repression, and is generally avoided by cells, while revealing one case in mitochondria where tight spatiotemporal control over m1A levels was adopted as a potential means of post-transcriptional regulation.

RNA editing in bacteria recodes multiple proteins and regulates an evolutionarily conserved toxin-antitoxin system.

Adenosine (A) to inosine (I) RNA editing is widespread in eukaryotes. In prokaryotes, however, A-to-I RNA editing was only reported to occur in tRNAs but not in protein-coding genes. By comparing DNA and RNA sequences of Escherichia coli, we show for the first time that A-to-I editing occurs also in prokaryotic mRNAs and has the potential to affect the translated proteins and cell physiology. We found 15 novel A-to-I editing events, of which 12 occurred within known protein-coding genes where they always recode a tyrosine (TAC) into a cysteine (TGC) codon. Furthermore, we identified the tRNA-specific adenosine deaminase (tadA) as the editing enzyme of all these editing sites, thus making it the first identified RNA editing enzyme that modifies both tRNAs and mRNAs. Interestingly, several of the editing targets are self-killing toxins that belong to evolutionarily conserved toxin-antitoxin pairs. We focused on hokB, a toxin that confers antibiotic tolerance by growth inhibition, as it demonstrated the highest level of such mRNA editing. We identified a correlated mutation pattern between the edited and a DNA hard-coded Cys residue positions in the toxin and demonstrated that RNA editing occurs in hokB in two additional bacterial species. Thus, not only the toxin is evolutionarily conserved but also the editing itself within the toxin is. Finally, we found that RNA editing in hokB increases as a function of cell density and enhances its toxicity. Our work thus demonstrates the occurrence, regulation, and functional consequences of RNA editing in bacteria.

Correction: Mitochondrial 16S rRNA Is Methylated by tRNA Methyltransferase TRMT61B in All Vertebrates.

[This corrects the article DOI: 10.1371/journal.pbio.1002557.].

Mitochondrial 16S rRNA Is Methylated by tRNA Methyltransferase TRMT61B in All Vertebrates.

The mitochondrial ribosome, which translates all mitochondrial DNA (mtDNA)-encoded proteins, should be tightly regulated pre- and post-transcriptionally. Recently, we found RNA-DNA differences (RDDs) at human mitochondrial 16S (large) rRNA position 947 that were indicative of post-transcriptional modification. Here, we show that these 16S rRNA RDDs result from a 1-methyladenosine (m1A) modification introduced by TRMT61B, thus being the first vertebrate methyltransferase that modifies both tRNA and rRNAs. m1A947 is conserved in humans and all vertebrates having adenine at the corresponding mtDNA position (90% of vertebrates). However, this mtDNA base is a thymine in 10% of the vertebrates and a guanine in the 23S rRNA of 95% of bacteria, suggesting alternative evolutionary solutions. m1A, uridine, or guanine may stabilize the local structure of mitochondrial and bacterial ribosomes. Experimental assessment of genome-edited Escherichia coli showed that unmodified adenine caused impaired protein synthesis and growth. Our findings revealed a conserved mechanism of rRNA modification that has been selected instead of DNA mutations to enable proper mitochondrial ribosome function.

LEMONS - A Tool for the Identification of Splice Junctions in Transcriptomes of Organisms Lacking Reference Genomes.

RNA-seq is becoming a preferred tool for genomics studies of model and non-model organisms. However, DNA-based analysis of organisms lacking sequenced genomes cannot rely on RNA-seq data alone to isolate most genes of interest, as DNA codes both exons and introns. With this in mind, we designed a novel tool, LEMONS, that exploits the evolutionary conservation of both exon/intron boundary positions and splice junction recognition signals to produce high throughput splice-junction predictions in the absence of a reference genome. When tested on multiple annotated vertebrate mRNA data, LEMONS accurately identified 87% (average) of the splice-junctions. LEMONS was then applied to our updated Mediterranean chameleon transcriptome, which lacks a reference genome, and predicted a total of 90,820 exon-exon junctions. We experimentally verified these splice-junction predictions by amplifying and sequencing twenty randomly selected genes from chameleon DNA templates. Exons and introns were detected in 19 of 20 of the positions predicted by LEMONS. To the best of our knowledge, LEMONS is currently the only experimentally verified tool that can accurately predict splice-junctions in organisms that lack a reference genome.

Mitochondrial Involvement in Vertebrate Speciation? The Case of Mito-nuclear Genetic Divergence in Chameleons.

Compatibility between the nuclear (nDNA) and mitochondrial (mtDNA) genomes is important for organismal health. However, its significance for major evolutionary processes such as speciation is unclear, especially in vertebrates. We previously identified a sharp mtDNA-specific sequence divergence between morphologically indistinguishable chameleon populations (Chamaeleo chamaeleon recticrista) across an ancient Israeli marine barrier (Jezreel Valley). Because mtDNA introgression and gender-based dispersal were ruled out, we hypothesized that mtDNA spatial division was maintained by mito-nuclear functional compensation. Here, we studied RNA-seq generated from each of ten chameleons representing the north and south populations and identified candidate nonsynonymous substitutions (NSSs) matching the mtDNA spatial distribution. The most prominent NSS occurred in 14 nDNA-encoded mitochondrial proteins. Increased chameleon sample size (N = 70) confirmed the geographic differentiation in POLRMT, NDUFA5, ACO1, LYRM4, MARS2, and ACAD9. Structural and functionality evaluation of these NSSs revealed high functionality. Mathematical modeling suggested that this mito-nuclear spatial divergence is consistent with hybrid breakdown. We conclude that our presented evidence and mathematical model underline mito-nuclear interactions as a likely role player in incipient speciation in vertebrates.

Parkin modulates heteroplasmy of truncated mtDNA in Caenorhabditis elegans.

Parkin, which is mutated in most recessive Parkinsonism, is a key player in the selective removal of damaged mitochondria via mitophagy. Damaged mitochondria may carry mitochondrial DNA (mtDNA) mutations, thus creating a mixed mtDNA population within cells (heteroplasmy). It was previously shown that Parkin over-expression reduced the level of heteroplasmic mutations that alter mitochondrial membrane potential in human cytoplasmic hybrids. However, it remained unclear whether Parkin serves a similar role at the entire living organism, and whether this role is evolutionarily conserved. Here, we show that mutation in the Caenorhabditis elegans orthologue of Parkin (pdr-1) modulates the level of a large heteroplasmic mtDNA truncation. Massive parallel sequencing revealed that the mtDNAs of C. elegans wild type and pdr-1(gk448) mutant strains were virtually deprived of heteroplasmy, thus reflecting strong negative selection against dysfunctional mitochondria. Therefore, our findings show that the role of Parkin in the modulation of heteroplasmy is conserved between human and worm and raise the interesting possibility that mitophagy modulates the striking lack of heteroplasmy in C. elegans.

The first Chameleon transcriptome: comparative genomic analysis of the OXPHOS system reveals loss of COX8 in Iguanian lizards.

Recently, we found dramatic mitochondrial DNA divergence of Israeli Chamaeleo chamaeleon populations into two geographically distinct groups. We aimed to examine whether the same pattern of divergence could be found in nuclear genes. However, no genomic resource is available for any chameleon species. Here we present the first chameleon transcriptome, obtained using deep sequencing (SOLiD). Our analysis identified 164,000 sequence contigs of which 19,000 yielded unique BlastX hits. To test the efficacy of our sequencing effort, we examined whether the chameleon and other available reptilian transcriptomes harbored complete sets of genes comprising known biochemical pathways, focusing on the nDNA-encoded oxidative phosphorylation (OXPHOS) genes as a model. As a reference for the screen, we used the human 86 (including isoforms) known structural nDNA-encoded OXPHOS subunits. Analysis of 34 publicly available vertebrate transcriptomes revealed orthologs for most human OXPHOS genes. However, OXPHOS subunit COX8 (Cytochrome C oxidase subunit 8), including all its known isoforms, was consistently absent in transcriptomes of iguanian lizards, implying loss of this subunit during the radiation of this suborder. The lack of COX8 in the suborder Iguania is intriguing, since it is important for cellular respiration and ATP production. Our sequencing effort added a new resource for comparative genomic studies, and shed new light on the evolutionary dynamics of the OXPHOS system.

RNA-DNA differences in human mitochondria restore ancestral form of 16S ribosomal RNA.

RNA transcripts are generally identical to the underlying DNA sequences. Nevertheless, RNA-DNA differences (RDDs) were found in the nuclear human genome and in plants and animals but not in human mitochondria. Here, by deep sequencing of human mitochondrial DNA (mtDNA) and RNA, we identified three RDD sites at mtDNA positions 295 (C-to-U), 13710 (A-to-U, A-to-G), and 2617 (A-to-U, A-to-G). Position 2617, within the 16S rRNA, harbored the most prevalent RDDs (>30% A-to-U and ∼15% A-to-G of the reads in all tested samples). The 2617 RDDs appeared already at the precursor polycistrone mitochondrial transcript. By using traditional Sanger sequencing, we identified the A-to-U RDD in six different cell lines and representative primates (Gorilla gorilla, Pongo pigmaeus, and Macaca mulatta), suggesting conservation of the mechanism generating such RDD. Phylogenetic analysis of more than 1700 vertebrate mtDNA sequences supported a thymine as the primate ancestral allele at position 2617, suggesting that the 2617 RDD recapitulates the ancestral 16S rRNA. Modeling U or G (the RDDs) at position 2617 stabilized the large ribosomal subunit structure in contrast to destabilization by an A (the pre-RDDs). Hence, these mitochondrial RDDs are likely functional.

Mitochondrial DNA variation, but not nuclear DNA, sharply divides morphologically identical chameleons along an ancient geographic barrier.

The Levant is an important migration bridge, harboring border-zones between Afrotropical and palearctic species. Accordingly, Chameleo chameleon, a common species throughout the Mediterranean basin, is morphologically divided in the southern Levant (Israel) into two subspecies, Chamaeleo chamaeleon recticrista (CCR) and C. c. musae (CCM). CCR mostly inhabits the Mediterranean climate (northern Israel), while CCM inhabits the sands of the north-western Negev Desert (southern Israel). AFLP analysis of 94 geographically well dispersed specimens indicated moderate genetic differentiation (PhiPT = 0.097), consistent with the classical division into the two subspecies, CCR and CCM. In contrast, sequence analysis of a 637 bp coding mitochondrial DNA (mtDNA) fragment revealed two distinct phylogenetic clusters which were not consistent with the morphological division: one mtDNA cluster consisted of CCR specimens collected in regions northern of the Jezreel Valley and another mtDNA cluster harboring specimens pertaining to both the CCR and CCM subspecies but collected southern of the Jezreel Valley. AMOVA indicated clear mtDNA differentiation between specimens collected northern and southern to the Jezreel Valley (PhiPT = 0.79), which was further supported by a very low coalescent-based estimate of effective migration rates. Whole chameleon mtDNA sequencing (∼17,400 bp) generated from 11 well dispersed geographic locations revealed 325 mutations sharply differentiating the two mtDNA clusters, suggesting a long allopatric history further supported by BEAST. This separation correlated temporally with the existence of an at least 1 million year old marine barrier at the Jezreel Valley exactly where the mtDNA clusters meet. We discuss possible involvement of gender-dependent life history differences in maintaining such mtDNA genetic differentiation and suggest that it reflects (ancient) local adaptation to mitochondrial-related traits.

Mitochondrial-nuclear co-evolution and its effects on OXPHOS activity and regulation.

Factors required for mitochondrial function are encoded both by the nuclear and mitochondrial genomes. The order of magnitude higher mutation rate of animal mitochondrial DNA (mtDNA) enforces tight co-evolution of mtDNA and nuclear DNA encoded factors. In this essay we argue that such co evolution exists at the population and inter-specific levels and affect disease susceptibility. We also argue for the existence of three modes of co-evolution in the mitochondrial genetic system, which include the interaction of mtDNA and nuclear DNA encoded proteins, nuclear protein - mtDNA-encoded RNA interaction within the mitochondrial translation machinery and nuclear DNA encoded proteins-mtDNA binging sites interaction in the frame of the mtDNA replication and transcription machineries. These modes of co evolution require co-regulation of the interacting factors encoded by the two genomes. Thus co evolution plays an important role in modulating mitochondrial activity. This article is part of a Special Issue entitled: Mitochondrial Gene Expression.