The mitochondrial transcriptome of the anglerfish Lophius piscatorius
BMC Research Notes volume 12, Article number: 800 (2019)
Analyze key features of the anglerfish Lophius piscatorius mitochondrial transcriptome based on high-throughput total RNA sequencing.
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.
The mitochondrial genome (mtDNA) gene content and organization is highly conserved among vertebrates . All species investigated to date encode the same 37 canonical gene products of 13 hydrophobic membrane proteins, 2 ribosomal RNAs (mt-rRNAs), and 22 transfer RNAs (tRNAs), as well as several non-canonical peptides and long noncoding RNAs (lncRNAs) . The corresponding mitochondrial transcriptomes are less studied and have mainly been investigated in a small number of vertebrates including some mammalian cells and tissues [3, 4] and in gadiform fishes [5, 6]. Only minor differences were noted between the mammals and fish. In general, three polycistronic transcripts initiated from two H-strand promoters (HSP1 and HSP2) and one L-strand promoter (LSP) are involved in mitochondrial gene expression. Whereas the highly abundant HSP1 transcript mainly generates mt-rRNAs, the HSP2 transcript is responsible for most messenger RNAs (mRNAs) and tRNAs. The LSP transcript generates one mRNA and eight tRNAs.
Atlantic cod mt-rRNAs are oligo-adenylated , and fold into similar secondary structures as in other fish species [7, 8]. Interestingly, several mitochondrial-derived peptides (MDP) have been proposed to be encoded on both strands of the mt-rRNA gene locus , and two MDPs (MOTS-c and Humanin) have coding potential in Atlantic cod . Mature tRNAs carry the non-template CCA at their 3′ ends and fold into the common tRNA patterns [7, 10]. Eleven mature mRNAs were found expressed in the Atlantic cod mitochondria, 10 from the HSP2 transcript and one from LSP, and two of the HSP2-specific mRNAs were bicistronic (ND4/4L and ATPase8/6) . All mRNAs, except the LSP-specific ND6 mRNA, were found polyadenylated.
Mitochondrial lncRNAs have been identified and investigated in Atlantic cod . Here, lncCR-H and lncCR-L correspond to different strands of the mitochondrial control region (CR). Both lncRNAs are clearly expressed and appear to generate small stable mitochondrial RNA (mitosRNA) [2, 6, 11, 12]. We recently reported low-level substitution heteroplasmy of the anglerfish Lophius piscatorius based on SOLiD deep sequencing . As part of a study to generate a full reference genome and transcriptome for L. piscatorius, we here present the complete mitochondrial genome and key features of the corresponding mitochondrial transcriptome.
Nucleic acid extraction and high-throughput sequencing
Lophius piscatorius tissue samples were collected from two specimens obtained by commercial fishery off the coast of Nordland County, Northern Norway, in 2015 (BF1) and 2017 (BF2). Total DNA from BF1 was extracted from muscle tissue and sequenced by the SOLiD5500 and Ion PGM platforms as described previously . Total DNA sequencing (head kidney) of BF2 using the Illumina HiSeqX platform was performed by Dovetail Genomics (Chicago, US) as a service . Total RNA from heart muscle tissue of specimen BF2 was isolated using QIAzol Lysis Reagent (QIAGEN, Hilden—Germany) according to the manufacturers protocol. Cellular rRNA was depleted from 1 μg of total RNA using the RiboMinus™ Eukaryote System v2 (Thermo Fisher Scientific, Waltham, MA—USA), and whole transcriptome library was constructed using the Ion Total RNA-seq kit v2 (Thermo Fisher Scientific) according to the manufacturers protocols. Manual template preparation on an Ion OneTouch™ 2 System (Thermo Fisher Scientific) and sequencing of two Ion 540™ chips on the Ion GeneStudio™ S5 System (Thermo Fisher Scientific) were carried out at our Genomics Platform (Nord University) according to the manufacturers protocols. The sequencing resulted in a total of 154,741,088 reads with a mean read length of 169 nt, corresponding to 26 billion nt.
RNA reads were quality trimmed with Cutadapt  using q20 as a threshold. The minimum read length was set to 50 nt. Trimmed RNA reads were then mapped to the BF2 mitochondrial genome with CLC Genomics Workbench v12 (QIAGEN). The “Length fraction” parameter was set to 0.9 and “Similarity fraction” to 0.96, requiring at least 90% of the read length to map with 96% similarity. Other parameters were set to their defaults. The resulting BAM file was coordinate sorted with SAMtools  and then processed with BEDTools  (genomecov command) to obtain a base level coverage of the mitogenome. Mean coverage for each gene and non-coding region was calculated from bed file. Alignments were visually examined to identify non-coding RNAs and polyA tails.
Canonical mitochondrial genes in L. piscatorius
Complete mitochondrial genome sequences of two L. piscatorius specimens were determined using the Ion PGM and SOLiD5500 technologies (BF1; 2532 times mean coverage; MF994812; ) and the Illumina HiSeqX pair-end reads (BF2; 7643 times mean coverage; MN240767). The circular mtDNA possesses the conventional gene content and organization typical in vertebrates (Fig. 1a). Among the nine polymorphic sites between BF1 and BF2, seven were located in protein coding genes, representing both synonymous and non-synonymous amino acid substitutions (Additional file 1: Table S1).
Mitochondrial transcripts from L. piscatorius BF2 were generated by Ion S5 sequencing. About 145.2 million quality-filtered total RNA reads were obtained, including 510,484 reads (0.35%) unambiguously identified as mitochondrial transcripts when mapped to the BF2 mitochondrial genome. Several features were noted when inspecting the mitochondrial transcripts and correlating the expression values to specific mitochondrial gene regions (Fig. 1b): (1) reads from mt-rRNA gene transcripts were 10–100 times more abundant than protein coding transcripts. This observation is likely underestimated due to rRNA depletion of input RNA. (2) Of coding transcripts, cytochrome oxidase subunits were the most abundant, with NADH dehydrogenase subunits and ATPase subunits transcripts being much less abundant. (3) Highly abundant lncRNAs mapping to opposite strands within the mitochondrial CR and cytochrome oxidase I gene (COI) were noted. (4) Most mRNAs were polyadenylated and lacked 5′ and 3′ untranslated regions (UTRs) (Additional file 2: Table S2). A notable exception was the 75 nt 3′UTR of the COI mRNA (see below). Secondary structure predictions of L. piscatorius mt-SSU rRNA (Additional file 3: Figure S1) and mt-LSU rRNA (Additional file 4: Figure S2) showed typical fish mitochondrial features [7, 8]. Secondary structure predictions of all 22 tRNAs (Additional file 5: Figure S3) followed the general pattern of fish mitochondrial tRNAs .
Non-canonical mitochondrial genes in L. piscatorius
The two CR specific lncRNAs (lncCR-H and lncCR-L), transcribed from opposite strands in a non-overlapping organization (Fig. 2a), have previously been reported in Atlantic cod [11, 12] and human . The L-strand specific lncCR-L was found to be 30 times more abundant than the L-strand specific ND6 mRNA (Fig. 1b). The vertebrate mitochondrial COI mRNA is unusual due to the presence of a structured 3′UTR. We identified a polyadenylated COI mRNA containing a 75-nt 3′UTR in L. piscatorius (Fig. 2b). RNA-Seq data revealed a highly abundant 178 nt antisense RNA to the 5′ end of COI mRNA (Figs. 1b and 2b), which appeared novel among vertebrate mitochondrial lncRNAs and named lncCOI.
MDPs have been reported in vertebrates, and the best characterized is the humanin peptide . The humanin gene is located within the mt-LSU rDNA locus. L. piscatorius contains a humanin-like open reading frame (ORF) in the mt-LSU rRNA Domain IV, at the exact same location as in Atlantic cod and human (Fig. 2c, left panel). Sequence analysis revealed the derived peptide sequence to be invariant within the Lophius genus, highly conserved among fishes, and well conserved between fish and mammals (Fig. 2c, right panel).
Here we provide the complete mitochondrial genome sequence and key features of the corresponding transcriptome of the anglerfish L. piscatorius. We found all canonical mitochondrial genes to be expressed. Mt-rRNAs were clearly more abundant than mRNAs. Two lncRNAs (lncCR-L and lncCR-H) mapped to the mitochondrial CR, a finding that corroborates recent reports of Atlantic cod and human cells [2, 18]. Interestingly, we identified a novel and highly abundant antisense RNA (lncCOI). Finally, we present feature support for the encoding of a humanin-like peptide within the mt-LSU rRNA.
Teleost fish mitochondria generate 10 mature mRNAs from a single primary transcript (HSP2) that subsequently are translated into 12 mitochondrial proteins in OxPhos complexes I, III, IV and V [2, 6]. Thus, the observed differences in transcript abundance may be explained by differential stability of individual mRNAs, and not by transcription initiation. Fish mitochondrial mRNAs contain no, or very short UTRs. A notable exception is the approximately 75-nt 3′UTR of the COI mRNA, which is conserved between fish species [2, 6] and mammals . A study in rat showed that the nuclear miR-181c was regulating COI mRNA stability in heart tissue by 3′UTR binding . A similar 75-nt 3′UTR was detected in the polyadenylated L. piscatorius COI mRNA. It is plausible, that the 3′UTR structure in L. piscatorius contributes to the COI mRNA stability.
A number of mitochondrial lncRNAs have been noted and characterized in vertebrates [reviewed in 2, 22, 23], but no lncRNA has so far been linked to COI gene sequences. Our observation of lncCOI appears novel among vertebrates. If the highly abundant lncCOI contributes to mRNA stability, translational regulation, or other mitochondrial roles is currently not known. We also detected two CR-specific lncRNAs (lncCR-L and lncCR-H) in L. piscatorius. lncCR-L corresponds to the 5′ end region of the LSP primary transcript and has been detected in Atlantic cod . lncCR-L appears homologous to the 7S RNA reported in human mitochondria more than three decades ago , that was recently shown to be aberrantly expressed in human cancer cells . Interestingly, lncCR-L was the most abundant non-ribosomal mitochondrial transcripts in L. piscatorius. lncCR-H, on the other hand, corresponds to the 3′ end region of the HSP2 primary transcript. It has been reported in Atlantic cod to be polyadenylated, to harbor a mirror tRNA, a noncoding intergenic spacer, and heteroplasmic tandem repeats [11, 12]. Similar to that of Atlantic cod, the L. piscatorius lncCR-H contains a mirror tRNA and a polyA tail. lncCR-L and lncCR-H may function as precursors for mitosRNAs , but their biological role has not been elucidated.
Reports in mammals conclude that the humanin peptide has important roles in cellular signaling [19, 25,26,27]. Previously we presented evidence supporting the encoding of humanin-like peptides in Domain IV of the mt-LSU rRNA in gadiform fishes , and similar features have recently been reported in avians . Here we show that several anglerfishes, including all Lophius species where mtDNA sequences are available, possess humanin-like ORFs. How vertebrate humanin is translated is under debate, but different scenarios may be considered; (1) The humanin ORF is recognized in mt-rRNA by mitochondrial ribosomes and translated in mitochondria. This scenario is supported by a recent study in rat . (2) Translation may also occur in cytosolic ribosomes, which would require mitochondrial export. Interestingly, a chimeric mt-LSU rRNA (lncRNA SncmtRNA) was reported to be expressed in human proliferating cells and localized in the cytoplasm and the nucleus [29, 30]. (3) Humanin may also be expressed from a nuclear copy of mt-LSU rRNA (Numt sequence). Studies from human cells provide support for the expression of nuclear-encoded humanin isoforms . The latter scenario may explain why most, but not all, fish species have intact humanin-like ORFs in Domain IV.
Our study provides a mitochondrial transcriptome resource from L. piscatorius heart muscle tissue. All mitochondrial genes were expressed, and different mRNAs had different abundances. Two lncRNAs mapped to the control region, we identified one novel lncRNA antisense to the COI mRNA, and the mt-LSU rRNA has the potential of coding a humanin-like peptide.
Mitochondrial RNA sequencing was performed in one tissue type in one individual and has to be considered as a snapshot of the mitochondrial transcriptome of L. piscatorius.
Availability of data and materials
The RNA-seq raw sequencing data accession number at NCBI’s Sequence Read Archived (SRA) is SRS5181095. Accession numbers of mitogenomes are available from GenBank under the Accession Number MF994812 (BF1) and MN240767 (BF2).
long noncoding RNA
mitochondrial small RNA
Boore JL. Animal mitochondrial genomes. Nucleic Acids Res. 1999;27:1767–80.
Jørgensen TE, Johansen SD. Expanding the coding potential of vertebrate mitochondrial genomes: lesson learned from the Atlantic cod. In: Mattila M, editor. Mitochondrial DNA—new insights. New York: InTech Open; 2018. p. 95–111.
Mercer TR, Neph S, Dinger ME, Crawford J, Smith MA, Shearwood A-MJ, Haugen E, Bracken CP, Rackham O, Stamatoyannopoulos JA, Filipovska A, Mattick JS. The human mitochondrial transcriptome. Cell. 2011;146:645–58.
Siira SJ, Spåhr H, Shearwood AMJ, Ruzzenente B, Larsson NG, Rackham O, Filipovska A. LRPPRC-mediated folding of the mitochondrial transcriptome. Nat Commun. 2017;8:1532.
Bakke I, Johansen S. Characterization of mitochondrial ribosomal RNA genes in gadiformes: sequence variation, secondary structure features, and phylogenetic implications. Mol Phylogen Evol. 2002;25:87–100.
Coucheron DH, Nymark M, Breines R, Karlsen BO, Andreassen M, Jørgensen TE, Moum T, Johansen SD. Characterization of mitochondrial mRNAs in codfish reveals unique features compared to mammals. Curr Genet. 2011;57:213–22.
Satoh TP, Miya M, Mabuchi K, Nishida M. Structure and variation of the mitochondrial genome of fishes. BMC Genomics. 2016;17:719.
Jørgensen TE, Karlsen BO, Emblem Å, Breines R, Andreassen M, Rounge TB, Nederbragt AJ, Jakobsen KS, Nymark M, Ursvik A, Coucheron DH, Jakt LM, Nordeide JT, Moum T, Johansen SD. Mitochondrial genome variation of Atlantic cod. BMC Res Notes. 2018;11:397.
Cobb LJ, Lee C, Xiao J, Yen K, Wong RG, Nakamura HK, Mehta HH, Gao Q, Ashur C, Huffman DM, Wan J, Muzumdar R, Barzilai N, Cohen P. Naturally occurring mitochondrial-derived peptides are age-dependent regulators of apoptosis, insulin sensitivity, and inflammatory markers. Aging. 2016;8:796–809.
Johansen S, Guddal PH, Johansen T. Organization of the mitochondrial genome of Atlantic cod, Gadus morhua. Nucleic Acids Res. 1990;18:411–9.
Jørgensen TE, Bakke I, Ursvik A, Andreassen M, Moum T, Johansen SD. An evolutionary preserved intergenic spacer in gadiform mitogenomes generates a long noncoding RNA. BMC Evol Biol. 2014;14:182.
Jørgensen TE, Karlsen BO, Emblem Å, Jaky LM, Nordeide JT, Moum T, Johansen SD. A mitochondrial long noncoding RNA in Atlantic cod harbors complex heteroplasmic tandem repeat motifs. Mitochondrial DNA Part A. 2019;30:307–11.
Dubin A, Jørgensen TE, Jakt LM, Moum T, Johansen SD. The mitochondrial genome of the European Anglerfish Lophius piscatorius express low-level substitution heteroplasmy. Ann Mar Biol Res. 2017;4:1019.
Dubin A, Jørgensen TE, Moum T, Johansen SD, Jakt LM. Complete loss of the MHC II pathway in an anglerfish, Lophius piscatorius. Biol Lett. 2019;15:20190594.
Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet. J. 2011;17:10–2. https://doi.org/10.14806/ej.17.1.200.
Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R, et al. The sequence alignment/map format and SAMtools. Bioinformatics. 2009;25:2078–9.
Quinlan AR, Hall IM. BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics. 2010;26:841–2.
Hedberg A, Knutsen E, Løvhaugen AS, Jørgensen TE, Perander M, Johansen SD. Cancer-specific SNPs originate from low-level heteroplasmic variants in human mitochondrial genomes of a matched cell line pair. Mitochondrial DNA Part A. 2019;30:82–91.
Guo B, Zhai D, Cabezas E, Welsh K, Nouraini S, Satterthwait AC, Reed JC. Humanin peptide suppresses apoptosis by interfering with Bax activation. Nature. 2003;423:456–61.
Slomovic S, Laufer D, Geiger D, Schuster G. Polyadenylation and degradation of human mitochondrial RNA: the prokaryotic past leaves its mark. Mol Cell Biol. 2005;25:6427–35.
Das SD, Ferlito M, Kent OA, Fox-Talbot K, Wang R, Liu D, Raghavachari N, Yang Y, Wheelan SJ, Murphy E, Steenbergen C. Nuclear miRNA regulates the mitochondrial genome in the heart. Circ Res. 2012;110:1596–603.
Dietrich A, Wallet C, Iqbal RK, Gualberto JM, Lotfi F. Organellar non-coding RNAs: emerging regulation mechanisms. Biochimie. 2015;117:48–62.
Zhao Y, Sun L, Wang RR, Hu J-F, Cui J. The effects of mitochondria-associated long noncoding RNAs in cancer mitochondria: new players in an old arena. Crit Rev Oncol Hematol. 2018;131:76–82.
Ojala D, Crews S, Montoya J, Gelfand R, Attardi G. A small polyadenylated RNA (7S RNA), containing a putative ribosome attachment site, maps near the origin of human mitochondrial DNA replication. J Mol Biol. 1981;150:303–14.
Lee C, Yen K, Cohen P. Humanin: a harbinger of mitochondrial-derived peptides? Trends Endocrin Metabol. 2013;24:222–8.
Paharkova V, Alvarez G, Nakamura H, Cohen P, Lee KW. Rat humanin is encoded and translated in mitochondria and is located to the mitochondrial compartment where it regulates ROS production. Mol Cell Endocrinol. 2015;413:96–100.
Zarate SC, Traetta ME, Codagnone MG, Seilicovich A, Reines AG. Humanin, a mitochondrial-derived peptide released by astrocytes, prevents synapse loss in hippocampal neurons. Front Aging Neurosci. 2019;11:123.
Moritz M, Degletagne C, Romestaing C, Duchamp C. Comparative genomic analysis identifies small open reading frames (sORFs) with peptide-encoding features in avian 16S rDNA. Genomics. 2019. https://doi.org/10.1016/j.ygeno.2019.06.026.
Burzio V, Villota C, Villegas J, Landerer E, Boccardo E, Villa LL, Martinez R, Lopez C, Gaete F, Toro V, Rodrigues X, Burzio LO. Expression of a family of noncoding mitochondrial RNAs distinguishes normal from cancer cells. Proc Natl Acad Sci USA. 2009;106:9430–4.
Fitzpatrick C, Bendek MF, Briones M, Farfan N, Silva VA, Nardocci G, Montecino M, Boland A, Deleuze JF, Villegas J, Villota C, Silva V, Lobos-Gonzalez L, Borgna V, Barrey E, Burzio LO, Burzio VA. Mitochondrial ncRNA targeting induces cell cycle arrest and tumor growth inhibition of MDA-MB-231 breast cancer cells through reduction of key cell cycle progression factors. Cell Death Dis. 2019;10:423.
Bodzioch M, Lapicka-Bodzioch K, Zapala B, Kamysz W, Kiec-Wilk B, Dembinska-Kiec A. Evidence for potential functionality of nuclearly-encoded humanin isoforms. Genomics. 2009;94:247–56.
We thank Dovetail Genomics for Illumina sequencing of total DNA, and Harald Lorentzen, Bofisk – Bodø, for providing L. piscatorius specimen samples.
This work was supported by an open grant from Nord University covering the study design, data collection and sequencing, analysis, and interpretation of data.
Ethics approval and consent to participate
Fish tissue samples were obtained at site of fisheries, and do not involve research on animals.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Polymorphic sites in the mitochondrial genome of L. piscatorius specimens BF1 and BF2.
5′ and 3′ sequence features of L. piscatorius mitochondrial mRNAs derived from RNA-seq reads.
Complete secondary structure diagram of L. piscatorius mitochondrial small subunit rRNA.
Complete secondary structure diagram of L. piscatorius mitochondrial large subunit rRNA. Polymorphic site between BF1 and BF2 is indicated in Domain I. Low-level heteroplasmic sites in BF1 are indicated in Domains I and VI.
Secondary structure diagram of L. piscatorius mitochondrial tRNAs. Anti-codon triplets and the non-template CCA are indicated.
About this article
Cite this article
Dubin, A., Jørgensen, T.E., Jakt, L.M. et al. The mitochondrial transcriptome of the anglerfish Lophius piscatorius. BMC Res Notes 12, 800 (2019). https://doi.org/10.1186/s13104-019-4835-6
- Antisense RNA
- Long noncoding RNA