An evolutionary preserved intergenic spacer in gadiform mitogenomes generates a long noncoding RNA

Human mtDNA-Encoded Long ncRNAs: Knotty Molecules and Complex Functions

Until a few decades ago, most of our knowledge of RNA transcription products was focused on protein-coding sequences, which were later determined to make up the smallest portion of the mammalian genome. Since 2002, we have learnt a great deal about the intriguing world of non-coding RNAs (ncRNAs), mainly due to the rapid development of bioinformatic tools and next-generation sequencing (NGS) platforms. Moreover, interest in non-human ncRNAs and their functions has increased as a result of these technologies and the accessibility of complete genome sequences of species ranging from Archaea to primates. Despite not producing proteins, ncRNAs constitute a vast family of RNA molecules that serve a number of regulatory roles and are essential for cellular physiology and pathology. This review focuses on a subgroup of human ncRNAs, namely mtDNA-encoded long non-coding RNAs (mt-lncRNAs), which are transcribed from the mitochondrial genome and whose disparate localisations and functions are linked as much to mitochondrial metabolism as to cellular physiology and pathology.

mtR_find: A Parallel Processing Tool to Identify and Annotate RNAs Derived from the Mitochondrial Genome

RNAs originating from mitochondrial genomes are abundant in transcriptomic datasets produced by high-throughput sequencing technologies, primarily in short-read outputs. Specific features of mitochondrial small RNAs (mt-sRNAs), such as non-templated additions, presence of length variants, sequence variants, and other modifications, necessitate the need for the development of an appropriate tool for their effective identification and annotation. We have developed mtR_find, a tool to detect and annotate mitochondrial RNAs, including mt-sRNAs and mitochondria-derived long non-coding RNAs (mt-lncRNA). mtR_find uses a novel method to compute the count of RNA sequences from adapter-trimmed reads. When analyzing the published datasets with mtR_find, we identified mt-sRNAs significantly associated with the health conditions, such as hepatocellular carcinoma and obesity, and we discovered novel mt-sRNAs. Furthermore, we identified mt-lncRNAs in early development in mice. These examples show the immediate impact of miR_find in extracting a novel biological information from the existing sequencing datasets. For benchmarking, the tool has been tested on a simulated dataset and the results were concordant. For accurate annotation of mitochondria-derived RNA, particularly mt-sRNA, we developed an appropriate nomenclature. mtR_find encompasses the mt-ncRNA transcriptomes in unpreceded resolution and simplicity, allowing re-analysis of the existing transcriptomic databases and the use of mt-ncRNAs as diagnostic or prognostic markers in the field of medicine.

Emerging functions of mitochondria-encoded noncoding RNAs

Mitochondria, organelles that harbor their own circular genomes, are critical for energy production and homeostasis maintenance in eukaryotic cells. Recent studies discovered hundreds of mitochondria-encoded noncoding RNAs (mt-ncRNAs), including novel subtypes of mitochondria-encoded circular RNAs (mecciRNAs) and mitochondria-encoded double-stranded RNAs (mt-dsRNAs). Here, we discuss the emerging field of mt-ncRNAs by reviewing their expression patterns, biogenesis, metabolism, regulatory roles, and functional mechanisms. Many mt-ncRNAs have regulatory roles in cellular physiology, and some are associated with, or even act as, causal factors in human diseases. We also highlight developments in technologies and methodologies and further insights into future perspectives and challenges in studying these noncoding RNAs, as well as their potential biomedical applications.

Non-coding 7S RNA inhibits transcription via mitochondrial RNA polymerase dimerization

The mitochondrial genome encodes 13 components of the oxidative phosphorylation system, and altered mitochondrial transcription drives various human pathologies. A polyadenylated, non-coding RNA molecule known as 7S RNA is transcribed from a region immediately downstream of the light strand promoter in mammalian cells, and its levels change rapidly in response to physiological conditions. Here, we report that 7S RNA has a regulatory function, as it controls levels of mitochondrial transcription both in vitro and in cultured human cells. Using cryo-EM, we show that POLRMT dimerization is induced by interactions with 7S RNA. The resulting POLRMT dimer interface sequesters domains necessary for promoter recognition and unwinding, thereby preventing transcription initiation. We propose that the non-coding 7S RNA molecule is a component of a negative feedback loop that regulates mitochondrial transcription in mammalian cells.

Complete mitochondrial genome of Episymploce splendens (Blattodea: Ectobiidae): A large intergenic spacer and lacking of two tRNA genes

The complete mitochondrial genome of Episymploce splendens, 15,802 bp in length, was determined and annotated in this study. The mito-genome included 13 PCGs, 20 tRNAs and 2 rRNAs. Unlike most typical mito-genomes with conservative gene arrangement and exceptional economic organization, E. splendens mito-genome has two tRNAs (tRNA-Gln and tRNA-Met) absence and a long intergenic spacer sequence (93 bp) between tRNA-Val and srRNA, showing the diversified features of insect mito-genomes. This is the first report of the tRNAs deletion in blattarian mito-genomes and we supported the duplication/random loss model as the origin mechanism of the long intergenic spacer. Two Numts, Numt-1 (557 bp) and Numt-2 (975 bp) transferred to the nucleus at about 14.15 Ma to 22.34 Ma, and 19.19 Ma to 24.06 Ma respectively, were found in E. splendens. They can be used as molecular fossils in insect phylogenetic relationship inference. Our study provided useful data for further studies on the evolution of insect mito-genome.

The mitochondrial transcriptome of the anglerfish Lophius piscatorius

Objective

Analyze key features of the anglerfish Lophius piscatorius mitochondrial transcriptome based on high-throughput total RNA sequencing.

Results

We determined the complete mitochondrial DNA and corresponding transcriptome sequences of L. piscatorius. Key features include highly abundant mitochondrial ribosomal RNAs (10–100 times that of mRNAs), and that cytochrome oxidase mRNAs appeared > 5 times more abundant than both NADH dehydrogenase and ATPase mRNAs. Unusual for a vertebrate mitochondrial mRNA, the polyadenylated COI mRNA was found to harbor a 75 nucleotide 3′ untranslated region. The mitochondrial genome expressed several non-canonical genes, including the long noncoding RNAs lncCR-H, lncCR-L and lncCOI. Whereas lncCR-H and lncCR-L mapped to opposite strands in a non-overlapping organization within the control region, lncCOI appeared novel among vertebrates. We found lncCOI to be a highly abundant mitochondrial RNA in antisense to the COI mRNA. Finally, we present the coding potential of a humanin-like peptide within the large subunit ribosomal RNA.

Mitogenomics of Central American weakly-electric fishes

Electric fishes are a diverse group of freshwater organisms with the ability to generate electric organ discharges (EODs) that are used for communication and electrolocation. This group (ca. 200 species) has originated in South America, and six species colonized the Central American Isthmus. Here, we assembled the complete mitochondrial genomes (mitogenomes) for three Central American electric fishes (i.e. Sternopygus dariensis, Brachyhypopomus occidentalis, and Apteronotus rostratus), and, based on these data, explored their phylogenetic position among Gymnotiformes. The three mitogenomes show the same gene order, as reported for other fishes, with a size ranging from 16,631 to 17,093 bp. We uncovered a novel 60 bp intergenic spacer (IGS) located between the COII and tRNA^Lys genes, which appears to be unique to the Apteronotidae. Furthermore, phylogenetic relationships supported the traditional monophyly of Gymnotiformes, with the three species positioned within their respective family. In addition, the genus Apteronotus belongs to the early diverging lineage of the order. Finally, we found high sequence divergence (13%) between our B. occidentalis specimen and a sequence previously reported in GenBank, suggesting that the prior mitogenome of B. occidentalis represents a different South American species. Indeed, phylogenetic analyses using Cytochrome b gene across the genus placed the previously reported individual within B. bennetti. Our study provides novel mitogenome resources that will advance our understanding of the diversity and phylogenetic history of Neotropical fishes.

A mitochondrial long noncoding RNA in atlantic cod harbors complex heteroplasmic tandem repeat motifs

A heteroplasmic tandem repeat (HTR) array occupies 100 to 300 bp of the mitochondrial DNA control region in the Atlantic cod, and recently we noted that the repeat appeared integrated in a polyadenylated mitochondrial long noncoding RNA. Here we provide a more detailed analysis of the mitochondrial HTR in the mitochondrial genome of 134 Atlantic cod specimens. We report all specimens to harbor mitochondrial HTRs in the control region, and identified 26 distinct variants among the 402 repeat motifs assessed. Whereas most specimens contained HTR profiles of 2-5 copies consisting of the same 40-bp motif, 22 specimens showed compound HTR arrays of at least two types of motifs present in the same mitochondrial DNA molecule. We found HTR profiles to be highly conserved between different tissue types of a single individual, and strictly maternally inherited in a mating experiment between parental Atlantic cod expressing different HTR profiles and array motifs. We conclude that mitochondrial heteroplasmy in the control region is very common in Atlantic cod, and results in length heterogenity of the long noncoding RNA lncCR-H.

Mitochondrial genome variation of Atlantic cod

OBJECTIVE

The objective of this study was to analyse intraspecific sequence variation of Atlantic cod mitochondrial DNA, based on a comprehensive collection of completely sequenced mitochondrial genomes.

RESULTS

We determined the complete mitochondrial DNA sequence of 124 cod specimens from the eastern and western part of the species' distribution range in the North Atlantic Ocean. All specimens harboured a unique mitochondrial DNA haplotype. Nine hundred and fifty-two polymorphic sites were identified, including 109 non-synonymous sites within protein coding regions. Eighteen variable sites were identified as indels, exclusively distributed in structural RNA genes and non-coding regions. Phylogeographic analyses based on 156 available cod mitochondrial genomes did not reveal a clear structure. There was a lack of mitochondrial genetic differentiation between two ecotypes of cod in the eastern North Atlantic, but eastern and western cod were differentiated and mitochondrial genome diversity was higher in the eastern than the western Atlantic, suggesting deviating population histories. The geographic distribution of mitochondrial genome variation seems to be governed by demographic processes and gene flow among ecotypes that are otherwise characterized by localized genomic divergence associated with chromosomal inversions.

Cancer-specific SNPs originate from low-level heteroplasmic variants in human mitochondrial genomes of a matched cell line pair

Low-level mitochondrial heteroplasmy is a common phenomenon in both normal and cancer cells. Here, we investigate the link between low-level heteroplasmy and mitogenome mutations in a human breast cancer matched cell line by high-throughput sequencing. We identified 23 heteroplasmic sites, of which 15 were common between normal cells (Hs578Bst) and cancer cells (Hs578T). Most sites were clustered within the highly conserved Complex IV and ribosomal RNA genes. Two heteroplasmic variants in normal cells were found as fixed mutations in cancer cells. This indicates a positive selection of these variants in cancer cells. RNA-Seq analysis identified upregulated L-strand specific transcripts in cancer cells, which include three mitochondrial long non-coding RNA molecules. We hypothesize that this is due to two cancer cell-specific mutations in the control region.

The Evolution of Dark Matter in the Mitogenome of Seed Beetles

Animal mitogenomes are generally thought of as being economic and optimized for rapid replication and transcription. We use long-read sequencing technology to assemble the remarkable mitogenomes of four species of seed beetles. These are the largest circular mitogenomes ever assembled in insects, ranging from 24,496 to 26,613 bp in total length, and are exceptional in that some 40% consists of non-coding DNA. The size expansion is due to two very long intergenic spacers (LIGSs), rich in tandem repeats. The two LIGSs are present in all species but vary greatly in length (114-10,408 bp), show very low sequence similarity, divergent tandem repeat motifs, a very high AT content and concerted length evolution. The LIGSs have been retained for at least some 45 my but must have undergone repeated reductions and expansions, despite strong purifying selection on protein coding mtDNA genes. The LIGSs are located in two intergenic sites where a few recent studies of insects have also reported shorter LIGSs (>200 bp). These sites may represent spaces that tolerate neutral repeat array expansions or, alternatively, the LIGSs may function to allow a more economic translational machinery. Mitochondrial respiration in adult seed beetles is based almost exclusively on fatty acids, which reduces the need for building complex I of the oxidative phosphorylation pathway (NADH dehydrogenase). One possibility is thus that the LIGSs may allow depressed transcription of NAD genes. RNA sequencing showed that LIGSs are partly transcribed and transcriptional profiling suggested that all seven mtDNA NAD genes indeed show low levels of transcription and co-regulation of transcription across sexes and tissues.

Variation and evolution of polyadenylation profiles in sauropsid mitochondrial mRNAs as deduced from the high-throughput RNA sequencing

BACKGROUND

Genes encoded in vertebrate mitochondrial DNAs are transcribed as a polycistronic transcript for both strands, which is later processed into individual mRNAs, rRNAs and tRNAs, followed by modifications, such as polyadenylation at the 3' end of mRNAs. Although mechanisms of the mitochondrial transcription and RNA processing have been extensively studied using some model organisms, structural variability of mitochondrial mRNAs across different groups of vertebrates is poorly understood. We conducted the high-throughput RNA sequencing to identify major polyadenylation sites for mitochondrial mRNAs in the Japanese grass lizard, Takydromus tachydromoides and compared the polyadenylation profiles with those identified similarly for 23 tetrapod species, featuring sauropsid taxa (reptiles and birds).

RESULTS

As compared to the human, a major polyadenylation site for the NADH dehydrogenase subunit 5 mRNA of the grass lizard was located much closer to its stop codon, resulting in considerable truncation of the 3' untranslated region for the mRNA. Among the other sauropsid taxa, several distinct polyadenylation profiles from the human counterpart were found for different mRNAs. They included various truncations of the 3' untranslated region for NADH dehydrogenase subunit 5 mRNA in four taxa, bird-specific polyadenylation of the light-strand-transcribed NADH dehydrogenase subunit 6 mRNA, and the combination of the ATP synthase subunit 8/6 mRNA with a neighboring mRNA into a tricistronic mRNA in the side-necked turtle Pelusios castaneus. In the last case of P. castaneus, as well as another example for NADH dehydrogenase subunit 1 mRNAs of some birds, the association between the polyadenylation site change and the gene overlap was highlighted. The variations in the polyadenylation profile were suggested to have arisen repeatedly in diverse sauropsid lineages. Some of them likely occurred in response to gene rearrangements in the mitochondrial DNA but the others not.

CONCLUSIONS

These results demonstrate structural variability of mitochondrial mRNAs in sauropsids. The efficient and comprehensive characterization of the mitochondrial mRNAs will contribute to broaden our understanding of their structural and functional evolution.