- Research note
- Open Access
Functional conservation of mitochondrial RNA levels despite divergent mtDNA organization
BMC Research Notes volume 13, Article number: 334 (2020)
Mitochondria-encoded ribosomal RNA (rRNA) genes in humans are expressed at a higher rate than protein coding genes of the mitochondria. The organization of the human mitochondrial genome (mtDNA) is amenable to differential expression of rRNAs as the rRNA encoding genes lie in tandem immediately downstream of the promoter-containing region. However, mtDNA is not organized in the same way as humans in all metazoans. In the nematode, Caenorhabditis elegans, the rRNA genes are on opposite sides of the mtDNA molecule and there are no obvious promoter sequences specific to the rRNA genes. Thus, we asked whether rRNA levels are higher relative to mRNAs in mitochondria of C. elegans as they are in humans.
Using droplet digital PCR, we discovered that steady-state mitochondrial rRNA transcript levels are approximately 120 times higher than the levels of mitochondrial mRNAs. These data demonstrate that despite the lack of conservation in mitochondrial genome organization, a high mitochondrial rRNA-to-mRNA ratio is a conserved feature of metazoans.
As some of the most ancient and essential structural components of ribosomes, rRNAs constitute the majority of total cellular RNA in both prokaryotes and eukaryotes. Mitochondria contain their own ribosomes, which are, in part, composed of two rRNAs encoded by the organellar genome, mtDNA. In humans, genes encoded by mtDNA are transcribed as polycistrons, which are subsequently processed into discrete transcripts . Higher mitochondrial rRNA levels relatively to mRNA are in part accounted for by about 2.5 times higher stability of structured rRNAs . Additionally, in humans, mitochondrial rRNAs are calculated to be synthesized at a rate 50 to 100 times higher than for mRNAs . Higher expression of mitochondrial rRNAs relative to mRNAs appears to be made possible by the organization of mtDNA [2, 3]. In humans, the two rRNA encoding genes are positioned in tandem, and lie directly downstream of a promoter-containing region  (Fig. 1a). This arrangement is thought to allow the rRNA genes to be transcribed independently from their own promoter called heavy strand promoter 1 (HSP1) [4,5,6]. Termination of the resulting polycistron immediately downstream of the second rRNA gene ensures differential expression of the rRNAs from the rest of mtDNA transcripts encoded on the heavy strand [7, 8]. A second promoter called HSP2 is believed to regulate transcription of a polycistron that corresponds to almost the entire length of mtDNA. While there may be some controversy surrounding this “two-promoter” hypothesis , it is clear that the mitochondrial rRNAs are synthesized at a higher rate than the mRNA genes. However, mtDNA organization is not the same across all metazoans. A notable example is that of the model species Caenorhabditis elegans, in which the two rRNA encoding genes in mtDNA are separated by protein coding genes and tRNAs, making it difficult to transcribe them independently from the rest of the genome  (Fig. 1b). Organization of mtDNA, similar to C. elegans, in which the rRNAs are distant from each other and not near any rRNA-specific promoters is seen in species across metazoan taxa (phylogeny adapted from  (Fig. 1c). Here, we ask whether, and to what extent, C. elegans have a high mitochondrial rRNA-to-mRNA ratio.
Materials and methods
N2 (C. elegans wild isolate) animals were maintained at 20 °C on standard nematode growth media (NGM) plates seeded with OP50 E. coli bacteria.
For each biological replicate, four 60 mm NGM plates of mixed-stage N2 worms were washed off of plates with M9 buffer into a single 15 mL conical tube. Worm suspensions were spun-down at 1300g for 5 min, supernatants were aspirated, and worm pellets were transferred to 2 ml RNase-free micro-centrifuge tube. Worm pellets were washed with 1 ml of M9 buffer two additional times to remove residual bacteria. Worm pellets (~ 100 μl) were snap frozen in liquid nitrogen. RNA was extracted using TRIzol™ (Thermo Scientific #15596026) as per manufacturer’s directions.
cDNA was synthesized using Maxima H Minus First Strand cDNA Synthesis Kit (Thermo Scientific #K1681) following manufacturer’s directions. Briefly, ~ 500 ng of RNA from each biological replicate was combined with 1 μl of dsDNase, 1 μl of 10 × dsDNase buffer and brought-up to 10ul with nuclease-free H2O. This mixture was incubated at 37 °C for 2 min and then immediately returned to ice. Then, to the RNA containing tube, 1 μl random hexamer primers, 1 μl of 10 mM dNTP, 4 μl 5x RTase buffer, 1 μl of Maxima H Minus reverse transcriptase, and 3 μl of nuclease-free H2O were added. The cDNA synthesis reaction was run in a thermocycler: 25 °C for 10 min, 55 °C for 30 min, 85 °C for 2 min, 4 °C hold. Following cDNA synthesis, cDNA was diluted to 50 μl with nuclease-free H2O and stored at − 80 °C.
Droplet digital PCR (ddPCR)
cDNA dilution: cDNA was diluted in nuclease-free water to run in the dynamic range of ddPCR. Initially, all targets were diluted 1:200. This dilution fell within the dynamic range for the protein-genes. However, at this dilution the rRNA levels saturated the capacity of ddPCR. From this initial result it was determined that a 100x further dilution would be sufficient to accurately quantify rRNA levels. Consequently, cDNA was diluted 1:20,000 for reactions that measured rRNA targets.
ddPCR reaction: In the wells of a Eppendorf™ 96-Well twin.tec™ PCR Plate (Fisher Scientific # E951020303), 2 μl of diluted cDNA was combined with 10 μl of nuclease-free water, 12.5 μl of QX200™ ddPCR™ EvaGreen Supermix (Bio-Rad #1864034), and 0.25 μl of each corresponding forward and reverse primer (10 μM working stock). Droplets were generated on a QXDx Automated Droplet Generator. Following droplet generation cDNA was amplified: 95 °C for 5 min, 39 cycles of 95 °C for 30 s then 58 °C for 1 min, 4 °C for 5 min, 90 °C for 5 min, 10 °C hold. Following cDNA amplification, droplets were read on a QXDx Droplet Reader.
Primers: 12S rRNA: CTTGTTCCAGAATAATCGGCTAGACTTG/CTAACCAGGTACTAATCTGCTTTGTTCAAC
16S rRNA: CAGTCTTAGCGTGAGGACATTAAGGTA/CTAACCAATAACTTCATTCATACTGGAACTC
Results and discussion
To determine the level of mitochondrial rRNAs and mRNAs in C. elegans we extracted RNA from mixed-stage populations of wildtype worms followed by first strand cDNA synthesis using random hexamer priming. We then used highly sensitive quantitative droplet digital PCR (ddPCR) to measure mitochondrial mRNA and rRNA steady-state levels. Upon initial quantification, the rRNA transcript counts exceeded the dynamic range of ddPCR while the mRNA levels did not. In fact, cDNA had to be diluted 100-fold further to precisely quantify rRNA levels. Using a generalized linear mixed effects model in R (package = ImerTest) we compared mitochondrial rRNA levels to mitochondrial mRNA levels after log-transformation of RNA copy number to achieve a normal distribution. We treated RNA type as a fixed effect and replicate and gene as random effects, to give a model of log (relative copy number) ~ rna type + (1|replicate) + (1|gene). Our analysis revealed that rRNA species levels are 122 times higher than mRNA species levels, after accounting for variance associated with the random effects. A Type II ANOVA on the fixed effect, using Satterthwaite’s method to approximate the degrees of freedom from the mixed model, provided an estimated F = 56.7 and p = 0.00028 (Fig. 1d).
Our finding that C. elegans mitochondrial rRNA levels are over 100 times higher than mitochondrial mRNAs suggests that a high mitochondrial rRNA-to-mRNA ratio is evolutionarily conserved. Furthermore, our finding suggests that multiple strategies have been employed across evolution to ensure a high mitochondrial rRNA-to-mRNA ratio. In humans, and many other metazoans, mtDNA organization provides the means for differential expression of the mitochondrial rRNA genes. However, the rRNA genes in C. elegans are separated by coding genes and tRNA genes and there are no apparent promoter sequences directly upstream of rRNA genes. Thus, the mechanisms to achieve high mitochondrial rRNA-to-mRNA ratio in humans are likely not available in C. elegans, forcing it to rely on potentially different mechanisms. One explanation is that C. elegans mtDNA contains cryptic promoter sequences within coding genes upstream of rRNA genes to facilitate specific expression of rRNAs. This has been proposed as a potential explanation for higher mitochondrial rRNA levels in sea urchins (S. purpuratus) . However, this model is not consistent with the report of a single transcriptional start site in C. elegans mtDNA . Alternatively, the post-transcriptional stability of rRNAs could be higher than that of mRNAs. Indeed, it has been shown that the half-lives of mitochondrial rRNAs in human cell lines are longer than that of mitochondrial mRNAs but not long enough to fully account for the difference in rRNA and mRNA levels . Such increased stability of rRNAs may be a product of their secondary structure or a result of rRNA association with ribosomal proteins or other RNA-binding proteins. Furthermore, we cannot exclude the possibility that other post-transcriptional mechanisms are involved in facilitating high mitochondrial rRNA levels. Future investigations promise to reveal insights into mitochondrial RNA-level regulation. Nevertheless, our finding that C. elegans mitochondrial rRNA levels are much higher than mitochondrial mRNAs strongly suggests that a high mitochondrial rRNA-to-mRNA ratio is evolutionarily conserved and is not restrained by mtDNA organization. We propose that C. elegans can be used as a complimentary model system to gain insights into the mechanisms employed to maintain a high mitochondrial rRNA-to-mRNA ratio.
The current study reports the observation that rRNA levels are higher than mRNA levels in mitochondria. Determining the mechanism underlying this difference will have to await future studies.
Availability of data and materials
All data generated and analyzed during this study are included in this published article.
Droplet digital PCR
Heavy strand promoter
Heavy strand promoter 1
Heavy strand promoter 2
Light strand promoter
Ojala D, Montoya J, Attardi G. tRNA punctuation model of RNA processing in human mitochondria. Nature. 1981;290:470–4.
Gelfand R, Attardi G. Synthesis and turnover of mitochondrial ribonucleic acid in HeLa cells: the mature ribosomal and messenger ribonucleic acid species are metabolically unstable. Mol Cell Biol. 1981;1(6):497–511. http://mcb.asm.org/content/1/6/497.abstract.
Mercer TR, Neph S, Dinger ME, Crawford J, Smith MA, Shearwood AMJ, et al. The human mitochondrial transcriptome. Cell. 2011;146(4):645–58. https://doi.org/10.1016/j.cell.2011.06.051.
Montoya J, Gaines GL, Attardi G. The pattern of transcription of the human mitochondrial rRNA genes reveals two overlapping transcription units. Cell. 1983;34(1):151–9.
Lodeiro MF, Uchida A, Bestwick M, Moustafa IM, Arnold JJ, Shadel GS, et al. Transcription from the second heavy-strand promoter of human mtDNA is repressed by transcription factor a in vitro. Proc Natl Acad Sci. 2012;109(17):6513–8. http://www.pnas.org/content/109/17/6513.abstract.
Zollo O, Tiranti V, Sondheimer N. Transcriptional requirements of the distal heavy-strand promoter of mtDNA. Proc Natl Acad Sci. 2012;109(17):6508–12. http://www.pnas.org/content/109/17/6508.abstract.
Kruse B, Narasimhan N, Attardi G. Termination of transcription in human mitochondria: Identification and purification of a DNA binding protein factor that promotes termination. Cell. 1989;58(2):391–7.
Hess JF, Parisi MA, Bennett JL, Clayton DA. Impairment of mitochondrial transcription termination by a point mutation associated with the MELAS subgroup of mitochondrial encephalomyopathies. Nature. 1991;351(6323):236–9. https://doi.org/10.1038/351236a0.
D’Souza AR, Minczuk M. Mitochondrial transcription and translation: overview. Garone C, Minczuk M, editors. Essays Biochem. 2018;62(3):309–20. https://doi.org/10.1042/EBC20170102.
Okimoto R, Macfarlane JL, Clary DO, Wolstenholme DR. The mitochondrial genomes of two nematodes, Caenorhabditis elegans and Ascaris suum. Genetics. 1992;130(3):471–98.
Elphick MR, Mirabeau O, Larhammar D. Evolution of neuropeptide signalling systems. J Exp Biol. 2018;221(3):jeb151092. http://jeb.biologists.org/content/221/3/jeb151092.abstract.
Elliott DJ, Jacobs HT. Mutually exclusive synthetic pathways for sea urchin mitochondrial rRNA and mRNA. Mol Cell Biol. 1989; 9(3):1069–82. http://mcb.asm.org/content/9/3/1069.abstract.
Blumberg A, Rice EJ, Kundaje A, Danko CG, Mishmar D. Initiation of mtDNA transcription is followed by pausing, and diverges across human cell types and during evolution. Genome Res. 2017;27(3):362–73.
We thank the members of the Patel Laboratory (Bryan L Gitschlag, Cassidy A Johnson, Cait S Kirby, Claudia V Pereira & Nikita Tsyba) for their valuable feedback on the manuscript. We thank Ann T Tate for help with statistics, and Rama D Gangula in the laboratory of Simon A Mallal for her assistance with ddPCR.
This work was supported by R01 GM123260 (M.R.P.) and by the support provided to J.P.H. by the Training Program in Environmental Toxicology (5T32ES007028-45).
Ethics approval and consent to participate
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.
About this article
Cite this article
Held, J.P., Patel, M.R. Functional conservation of mitochondrial RNA levels despite divergent mtDNA organization. BMC Res Notes 13, 334 (2020). https://doi.org/10.1186/s13104-020-05177-0
- Ribosomal RNA
- Droplet digital PCR
- Caenorhabditis elegans