Figures
Abstract
Small non-coding RNAs (ncRNAs), particularly miRNAs, play key roles in a plethora of biological processes both in health and disease. Although largely operative in the cytoplasm, emerging data indicate their shuttling in different subcellular compartments. Given the central role of mitochondria in cellular homeostasis, here we systematically profiled their small ncRNAs content across mouse tissues that largely rely on mitochondria functioning. The ubiquitous presence of piRNAs in mitochondria (mitopiRNA) of somatic tissues is reported for the first time, supporting the idea of a strong and general connection between mitochondria biology and piRNA pathways. Then, we found groups of tissue-shared and tissue-specific mitochondrial miRNAs (mitomiRs), potentially related to the “basic” or “cell context dependent” biology of mitochondria. Overall, this large data platform will be useful to deepen the knowledge about small ncRNAs processing and their governed regulatory networks contributing to mitochondria functions.
Citation: Siniscalchi C, Di Palo A, Petito G, Senese R, Manfrevola F, Leo ID, et al. (2024) A landscape of mouse mitochondrial small non-coding RNAs. PLoS ONE 19(1): e0293644. https://doi.org/10.1371/journal.pone.0293644
Editor: Eduardo Andrés-León, Institute of Parasitology and Biomedicine, SPAIN
Received: January 31, 2023; Accepted: October 17, 2023; Published: January 2, 2024
Copyright: © 2024 Siniscalchi et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: All Sequencing files are available from ArrayExpress database (https://www.ebi.ac.uk/biostudies/arrayexpress/studies/E-MTAB-12911?query=E-MTAB-12911). All relevant data are within the paper and its Supporting Information files.
Funding: The research was funded by the University of Campania “Luigi Vanvitelli” within VALERE 2019 program giving support to the ICOM project. The funder had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
The recent advances in high-throughput sequencing technologies shed light on a plethora of RNA species with limited coding potential, collectively called non-coding RNAs (ncRNA) [1]. Besides ribosomal RNA, ncRNA molecules are generally classified according to their size threshold in long ncRNA (lncRNAs), longer than 200 nt up to several kilobases, and small or short ncRNAs, from few to 200 nt. lncRNA are involved in the regulation of gene expression at transcriptional and post-transcriptional levels by interacting with DNA, proteins and other RNA molecules. Short ncRNAs include the well known tRNA (transfer RNA), involved in the translation of mRNA; snRNA (small nuclear RNA), involved in splicing; snoRNA (small nucleolar RNA), implicated in ribosomal RNA modification; Piwi-interacting RNA (piRNA) and microRNAs (miRNA) [2].
piRNAs, 24–35 nucleotide long, 2’-O-methyl modified at 3’-end, associate with PIWI proteins; they were discovered in germlines in 2001, representing the most abundant small RNAs in mammalian testes and are essential for fertility; they are classically associated with “genome defense” against transposable elements and viruses [3–5]; more recently, they have been involved in epigenetic modification and gene expression regulation, probably at lower extent also in somatic tissues [6]. piRNA biogenesis and functioning are still not completely known, probably due to the occurrence of different production pathways, the wide diversity of functions (in the nucleus and cytosol) and poorly conserved pathways specifically associated to features of germ cell development [6, 7].
miRNAs, the most studied group of small ncRNAs, with a size of about 20 nt, negatively regulate gene expression by binding to target transcripts, thus affecting their translation and/or stability [8]. In mammals, up to 50% of all protein-coding transcripts are predicted to be targeted by miRNAs; in addition, one miRNA can bind different mRNAs and one mRNA can targeted by different miRNAs, giving rise to complex regulatory networks that play key roles in almost all physiological pathways, but also in the pathogenesis of several diseases, especially cancer [9–11]. miRNAs biogenesis is compartmentalized between nucleus and cytoplasm, but they are largely operative in the cytoplasm. However, emergent data indicate that mature miRNAs can localize in the nucleus and in subcellular compartments in the cytoplasm, unveiling possible new and specialized functions related to specific cellular compartments [12]. In this regard, some studies have detected microRNAs in mitochondria isolated from rat, mouse and human cells, collectively called mitomiRs (mitochondria-localized microRNA) [13, 14]. In fact, in 2006 a few miRNAs were found in rat’s liver mitochondria; although initially thought as cytosolic contamination, few years later other studies performed by microarray profiling and more recently by RNAseq confirmed the localization of miRNAs in mitochondria [15–22].
Mitochondria are the powerhouses in the cell and essential organelles for the integration of several key metabolic processes such as respiration, fatty and amino acid metabolism [23]. They can be defined as central regulators of cellular homeostasis, and dysregulation of any activity has pathological consequences, like cancer, inflammation, neurodegeneration [24–27]. Mitochondria contain their own circular genome, about 16 kbp in length, with 37 genes producing 2 rRNAs, 22 tRNAs and 13 proteins of the oxidative phosphorylation (OXPHOS) system [28]. Mitochondrial homeostasis is dependent on the expression of its own genes, quite limited in number, but largely on the import of nuclear-encoded proteins from the cytoplasm [13]. It is becoming increasingly clear that nuclear genome contributes to mitochondria activities not only through proteins (about 1500 nuclear-encoded proteins), but also through a number of miRNAs that can affect transcripts encoding mitochondrial-localized and mitochondrial-related proteins. This crosstalk between nucleus and mitochondria is essential for regulating mitochondria dynamics (biogenesis, fusion, fission, mitophagy), metabolic pathways and mitochondria-mediated apoptosis, whose dysregulations have pathological consequences. Indeed, some miRNAs can be translocated to mammalian mitochondria. The proposed import mechanisms include a detachment of AGO2 and miRNA from RISC or just GW182 due to AGO2 phosphorylation or some other signal activation promoting their translocation together or separately; this process could be stimulated by P-bodies; translocation across mitochondrial membranes is unknown but suggested to be promoted by PNPase and occurs within TOM/TIM complexes, or relies on VDAC at the outer mitochondrial membrane [13, 29, 30].
Different gaps in the knowledge emerge from literature related to ncRNA contribution to mitochondria biology. As an example, mitomiR function still remains controversial, especially due to debated presence of RISC components required for their canonical mechanism of silencing activity [12, 31]. Moreover, a deep sequencing analyses of small RNA content in mitochondria, probably a foundation for a systematic understanding, is available only for HEK293 and HeLa human cells, and human myotubes [20, 22], and regarding to mouse model only for primordial male germ cells, spermatogonia, spermatozoa, oocytes and zygotes [21]. To gain insights into possible new and specialized function of small RNAs and in particular miRNAs related to mitochondrial biology, this work explores the content of small RNAs population by RNAseq of mitochondria isolated from different mouse tissues, also in order to identify a possible shared pool and distinctive mitomiR signatures potentially involved in basic or cell-context dependent mitochondria functions.
Materials and methods
Animals and ethics statement
Three wild-type male CD1 mice, ∼3 months old, were used in this study. They were maintained in a temperature-controlled room at 22 ± 2°C under a 12 h dark/light cycle with a standard pellet diet and unlimited access to water.
Mice were sacrificed under anesthesia by cervical dislocation. In detail, animals were placed in a plexiglas chamber with 4% isoflurane (Iso-Vet, Piramal Healthcare, United Kingdom) for 5 min and were sacrificed when fully sedated, as measured by a lack of heartbeat and active paw reflex. Tissues were rapidly removed, rinsed in PBS (pH 7.6) to remove blood contaminants and used for experimental procedures. Every effort was made to minimize animal pain and suffering.
Experiments were approved by the Italian Ministry of Education and the Italian Ministry of Health (authorization n°48/2022-PR issued on 28.01.2022). Procedure involving animal care were carried out in accordance with the National Research Council’s publication Guide for Care and Use of Laboratory Animals (National of Institutes of Health Guide).
Mitochondria purification
Mitochondria were purified from freshly excised tissues (liver, gastrocnemius skeletal muscle, testis, white adipose tissue) using Qproteome Mitochondrial Isolation kit (Qiagen) according to manufacturer’s instructions with minor modifications and performing all the steps at 4°C. Briefly, tissues were homogenized in Lysis Buffer supplemented with Protease Inhibitor solution by Potter homogenizer, incubated on an end-over-end shaker for 10 min and centrifuged at 1000 x g for 10 min. The pellet was resuspended in ice-cold Disruption Buffer supplemented with Protease Inhibitor, subjected to Potter homogenizer for cell disruption, followed by centrifugation at 1000 x g for 10 min. The supernatant was centrifuged at 6000 x g for 10 min to pellet down mitochondria. The mitochondria were further purified by resuspending the pellet in Purification buffer and carefully pipetting onto layers of Purification Buffer and Disruption Buffer. The solution was centrifuged at 14000 x g for 15 min. The mitochondrial ring at the interface of Purification Buffer/Disruption Buffer was collected and washed in Mitochondria Storage Buffer; the final mitochondrial pellet was resuspended in 20 μl Storage buffer.
In order to remove RNA of nuclear origin present in the cytosol and eventually adsorbed on the outer mitochondrial membrane, mitochondria were treated with RNase A solution, as already described [20, 22]. In brief, the 20 μl purified mitochondria from the above described procedure were mixed with 300 μl PBS containing RNase A (100 μg/ml) and incubated at 37°C for 1 h. Rnase A treated mithocondria were pelleted down at 8000 rpm for 10 min, washed twice in PBS and finally resuspended in Qiazol buffer (Qiagen) for RNA extraction.
RNA extraction, control of mitochondria purity and real-time PCR analyses
Mitochondria and cytosolic fractions obtained by the procedure described above according to Qproteome Mitochondrial Isolation kit (Qiagen) protocol were used for total RNA purification by miRNeasy Mini Kit (Qiagen) according to the manufacturer’s protocol. The yield and quality of RNA were evaluated by NanoDrop One spectrophotometer (NanoDrop), Qubit 4 Fluorometer (Invitrogen) and by TapeStation 4200 (Agilent Technologies).
The mitochondrial RNA enrichment and eventual cytosolic RNA contamination were checked by RT-qPCR. Mitochondria and cytosolic RNAs were retrotranscribed by SensiFAST cDNA Synthesis kit (Bioline). Then standard SYBR Green Real-time qPCR assays were performed for amplifying β-actin transcript as nucleus-derived RNA control and Cytochrome c as mitochondrion-derived control with the following primers: β-actin, 5’-CAACGGCTCCGGCAT -3’ and 5’- CTCTTGCTCTGGGCC -3’; Cytochrome b, 5′-TACCTGCCCCATCCAACATT-3′ and 5′-TAAGCCTCGTCCGACATGAA-3′. For the further experimental steps, we selected the RNA mitochondrial preparations giving the following qPCR parameters on triplicate analyses: Ct>35 for β-actin, which was fully detectable in the parallel cytosolic sample, indicating the absence of cytosolic contaminant RNA; Ct<28 for cytochrome b, indicating a good enrichment of mitochondrial RNA.
miR-122-5p, miR-133a-3p, miR-29a-3p, miR-15b-5p were quantified along with U6 (reference transcript) by RT-qPCR with TaqMan® miRNA assays from Applied Biosystems according to the manufacturer’s protocol. The expression levels of miRNAs were normalized to the reference transcript by using the 2-ΔCt method.
RNA sequencing and data analyses
Indexed libraries were prepared from 50 ng of each purified RNA sample with NEXSmallRNA Seq v3.0 (PerkinElmer) according to the manufacturer’s instructions. Libraries were quantified using the TapeStation 4200 (Agilent Technologies) and Qubit 4 Fluorometer (Invitrogen). Libraries were pooled in equimolar amounts with final concentration of 2 nM. The pooled samples were subject to cluster generation and sequencing using an Illumina NextSeq550Dx System (Illumina Inc.) in a 1x75 single-end format with an average depth of the sequences of 20 Million reads/sample. The raw sequence files generated (.fastq files) underwent quality control analysis using FastQC (http://www.bioinformatics.babraham.ac.uk/projects/fastqc/). sRNAbench by sRNAtoolBox was used to identify the non-coding RNAs signature in all sequenced samples, thus removing adapter sequences, based on reference kit instructions, and the low quality reads [32]. In particular, the following steps were performed: 4 nt barcode from 5’ end of reads were removed; the adapter sequence "TGGAATTCTCGGGTGCCAAGGG" were detected; minimum length of adapter sequence aligned considered as 10 and only 1 mismatch allowed between adapter and reference sequences; 4 nt random adapter from 3’ end of adapter trimmed reads were removed.
The samples were mapped respect to miRBase database (version22 –GRCh38 ‐ GCA_000001405.15) to Ensembl and to RNACentral annotatio, to obtain several classes of ncRNAs allowing only 1 mismatch. piRNA-related reads were also mapped by Bowtie2 on the IPpiRNAdb (https://github.com/OdBT/piRNA-IPdb_v2) [33] obtaining an overall alignment rate equal to 99.05% for gastrocnemius skeletal muscle samples, 98.81% for liver samples, 99.32% for testis samples, 98.85% for WAT samples.
Finally, reads were also aligned to mouse mitochondrial genome (NC_005089.1 from NCBI) by blastn command line applications with defaults parameters, finding reads with at least one reported alignment in the following percentages: 4.87% for gastrocnemius muscle, 6.17% for liver, 4.00% for testis and 3,47% for WAT. e-value≤0.01 was considered statistically significant. Finally, all piRNA sequences of piRbase (http://bigdata.ibp.ac.cn/piRBase/) and all miRNA sequences included in miRbase (https://www.mirbase.org/) were aligned to mitochondrial genome with parameters described above.
Sequencing data are available in the ArrayExpress database repository (https://www.ebi.ac.uk/arrayexpress/) with the following accession number: E-MTAB-12911.
mitomiR targetome analyses
The list of 126 tissue-shared mitomiRs and the lists of tissue-specific mitomiRs were used to launch a search on miRWalk platform (http://mirwalk.umm.uni-heidelberg.de/) to retrieve all their target transcripts. In detail, by selecting “miRTarbase”, or “TargetScan” and “miRDB”, and “5′UTR”, or “3′ UTR” or “coding region” and selecting 0.5 as score for miRNA–mRNA pairings, 6 different target lists for tissue-shared and the 4 tissue-specific mitomiRs subset were retrieved containing all experimentally validated targets (by miRTarbase tool), or those predicted by both TargetScan and miRDB tools. By exporting the lists in Excel, a list of non-redundant transcripts representing the whole targetome of tissue-shared and 4 lists of tissue-specific mitomiRs were produced. On those lists, enrichment of biological pathways supplied by Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway annotation was performed by Database for Annotation, Visualization and Integrated Discovery (DAVID, https://david.ncifcrf.gov) which facilitates the transition from data collection to biological analysis by systematically extracting biological meaning and gene functional groups from large gene lists. Considering this type of analysis could increase the false-positive rate, to control the false-positive rate in the results, a multiple test correction of enrichment p values was performed on the functional annotation categories by selecting the Benjamini–Hochberg procedure. Biological pathways with p values < 0.05 were considered statistically significant [34].
Targetomes of tissue-shared or–specific mitomiRs were also analyzed by MitoCarta3.0 to extract enriched biological pathways specifically related to mithocondria biology. In brief, targetome lists were matched with 1140 mouse genes encoding proteins with strong support of mitochondrial localization; the merging target transcript lists, plus those mitochondrial transcripts reported in PubMed as experimentally validated target of our detected mitomiRs, were analyzed by MitoCarta3.0 for the assignment to 149 hierarchical ’MitoPathways’ spanning seven broad functional categories relevant to mitochondria (http://www.broadinstitute.org/mitocarta) [35].
Results
Rationale of the study
In order to have a large platform of data regarding the population of small RNAs in mitochondria, arguing their potential contribution to mitochondrial activities, novel mechanisms of regulation and crosstalk with other cellular components, we decided to select (Fig 1):
- liver, since it can be seen as the metabolic hub of the body, in analogy to mitochondria being the metabolic hub of the cell; in fact, the liver is extremely rich in mitochondria, with each hepatocyte containing 1000–2000 mitochondria [36];
- gastrocnemius skeletal muscle, whose health is dependent on the optimal function of its mitochondria [37]; due to the high metabolic requirements for locomotion, mitochondrial ATP production is very important, with typically 3–8% of the skeletal muscle volume being mitochondria [38];
- testis, since spermatogenesis relies on mitochondrial dynamics and functions, whose disruption causes male infertility [39];
- white adipose tissue (WAT), where mitochondrial activity is essential for the choice of carbon source (use of carbohydrate through glycolysis versus use of lipids for oxidative phosphorylation) and impacts on tissue and systemic energy homeostasis; the important role of mitochondria in the regulation of white adipose tissue remodeling and energy balance is increasingly appreciated [40].
Illustration created with Servier Medical Art, provided by Servier, licensed under a Creative Commons Attribution 3.0 Unported license.
Mitochondria were purified as described in Material and Methods section, and smallRNA sequencing experiments were performed on extracted RNAs (Fig 1). The percentage of detected mitochondrial smallRNA populations on the annotated ones for each tissue are reported in Fig 2 and detailed below.
The light-blue indicates RNA biotypes percentages detected in indicated mitochondria tissue on the annotated ones in the different databases, i.e. RNA-Central, Ensembl (GRCm38_p5) and miRBase (v22). Misc_RNA: miscellaneous RNA; otherRNA: snRNA, snoRNA and RNA fragments are included.
mitopiRNAs
piRNAs are largely present in the mitochondria of testis, representing the most predominant population of small RNAs (Fig 2); their finding into mitochondria of testis is consistent with bioinformatic analyses on smallRNAseq libraries from primordial germ cells and spermatozoa of male mice [21, 41]. Surprisingly, here piRNAs were also detected in the mitochondria of the other tissues, representing a novelty in the current knowledge (S1 Table). In fact, although their presence has been already described for different somatic tissues [42] and related to epigenetic reprogramming and gene expression regulation [3, 43], piRNAs specifically located in mitochondria are only reported in mammalian cancer cells [22, 44] and inferred by bioinformatic analyses on RNA extracted from the whole cell lysates derived from murine primordial germ cells, spermatozoa oocytes, zygotes and embryonic gonadal cells [21, 41]. Two studies have demonstrated a role of the mitochondrial protein MitoPLD/Zucchini (mammalian/Drosophila homologues) in piRNA biogenesis and piRNA mediated silencing in mouse and fly germlines [45, 46]; finally, PIWIL-1, a human homologue of the mouse MIWI protein, has been detected in the mitochondria of human cells [44] pointing for a mitochondrial role in the biogenesis of piRNAs. Our discovery of piRNA populations in mitochondria of somatic tissues strongly supports the idea that the organelles may be an important and general player in piRNA pathways, beyond germ line tissues. Intriguingly, their expression pattern appears different for analyzed tissues, showing their massive presence in testis and some overlapping among analyzed tissues (Fig 3). Analogously to the so-called mitomiR, those piRNAs can be indicated as mitopiRNAs.
(A) Normalized expression unsupervised heatmap relative to piRNAs in the mitochondria of the indicated tissues, obtained by using the software ComplexHeatmap, a package of R. (B) The Venn diagram displaying the numbers of overlapping mitopiRNAs (Raw Count Mean ≥1) among the different tissues.
mitomiRNAs
The RNAseq analyses revealed similar percentages of miRNAs on the annotated ones in the mitochondria of gastrocnemius muscle, liver and WAT and a higher content in testis (Fig 2).
The overall expression pattern of mitomiRs is distinctive for the different tissues and ranges from few read counts to more than 1000, whose frequency counts is smaller compared to others (Fig 4A and 4B). We hypothesized that some miRNAs found in the mitochondria could be related to the basic functions of the organelle, whereas others could be mainly associated to specialized functions depending on the tissue. To explore that hypothesis, all the miRNAs found in the mitochondria of each tissue, from very low (> 1) to higher read counts were analyzed, resulting in a core of 126 shared mitomiRs, some others overlapping among the tissues and a number of tissue-specific mitomiRs (Fig 4C and S2 Table). Different papers and related datasets consistently report their expression also in tissues [47–49], and for some of them targets related to mitochondria functioning have been reported (Tables 1–5), thus further connecting their relevance for mitochondria biology.
(A) Normalized expression unsupervised heatmaps relative to mitomiRs of the indicated tissues, obtained by using the software ComplexHeatmap, a package of R. (B) Distribution of mitomiR levels with respect to frequency; numbers of sequence reads (showed on the bars) were taken as miRNA levels and the values are represented in the form of range of values in mitochondrial sRNA libraries of each tissue. (C) The Venn diagram displaying the numbers of overlapping mitomiRs among the different tissues.
It would be interesting to address the question whether miRNAs found in mitochondria reflect the general content of the cells or are differentially distributed between cytosol and mitochondria. To gain some insights into the issue, we quantify the level of specific miRNAs in cytosol and mitochondria fractions; in particular, we analyzed the level of the 1st and the 10th of tissues-shared top-10 mitomiRs according to their reads count mean, i.e. miR-122-5p and miR-133a-3p (Table 1), and two other miRNAs, miR-29a and miR-15b, respectively 12- and 17-fold less expressed in comparison to miR-133a (S2 Table). We found that in most cases microRNAs are differentially distributed between cytosol and mitochondria, appearing significantly enriched in mitochondria, (e.g., miR-122 in gastrocnemius, miR-29a in liver, miR-122 and miR-15b in testis, miR-133a in WAT), or in lower amount in comparison to the cytosol (e.g., miR-133a, miR-29a and miR-15b in gastrocnemius, miR-122 in liver, miR-29a in testis and WAT, miR-15b in WAT), pointing to a potential relevance of selective translocation into the organelles (Fig 5).
Expression analyses of the indicated miRNAs were performed by using Taqman assays on the RNA extracted from cytosol and mitochondria fractions of the indicated tissues. All analyses were performed at least in triplicate on each sample and mean + sd is reported. p-values at Student’s t-test were *p < 0.05, **p < 0.01, or ***p < 0.001.
mitomiR targetomes and pathways
miRNAs exert their function by down-regulating multiple target transcripts, since their pairing can tolerate mismatches, beyond the perfect complementarity required for miRNA seed sequence (nucleotides 2–8). By screening the published papers on PubMed (searching keywords: “miRNA name AND mitochondria”), we found that various mitomiRs detected in our analyses have experimentally validated targets related to mitochondria function (Tables 1–5); however, gaps in the knowledge emerge for many others. To gain insights into biological functions of identified mitomiRs, the whole targetomes of all tissue-shared and all tissue-specific mitomiRs and related biological pathways were analyzed. In brief, we obtained a list of 4840 non-redundant transcripts experimentally validated or predicted by both TargetScan and miRDB tools as targets of the 126 tissue-shared mitomiRs; a list of 705 non redundant transcripts for the 21 mitomiR specifically found in gastrocnemius muscle; a list of 632 for the 30 mitomiRs for liver; a list of 4897 for the 174 testis-specific mitomiRs; a list of 520 for 18 WAT-specific mitomiRs.
Enrichment of biological pathways supplied by KEGG on the targetome of tissue-shared mitomiRs shows an overrepresentation of different signaling pathways among the top p-values scored, beyond “Pathways in cancer”, the most studied field in relation to microRNAs (Fig 6 and S3 Table). That observation is also applicable on KEGG analysis performed on targetome of testis specific mitomiRs, although including some different signaling pathways. The enrichment of biological pathways implied in different signaling routes was probably expected for miRNAs associated to mitochondria, since the central role of the organelle in metabolism and homeostasis and the results here reported highlight the potential strong miRNA contribution. Behind those pathways enriched for mitomiRs arising from all analyzed tissues, others are specifically found for tissue-specific mitomiRs; e.g, “insuline resistance” and “Type II diabetes mellitus” pathways were found for liver-specific mitomiRs; “Wnt signling pathway”, well-known to be related to gluconeogenesis, was found for WAT; “Regulation of actin cytoskeleton” pathway emerges significantly enriched for mitomiRs in testis, where actin remodeling contributes to proper spermatogenesis in terms of morphology and then motility [74]. Overall, those results support the idea that regulation via miRNAs could be functionally relevant for mitochondria biology, since a pool of miRNAs could be related to the basic function of mitochondria shared by different cell types, and others may be related to the specialized functions of each tissue.
The whole targetomes (predicted and validated) of tissue-shared or tissue-specific mitomiRs were analyzed by the program DAVID. Significantly enriched pathways according to KEGG are showed for tissue-shared (A; top-20 p values), liver-specific (B), testis-specific (top-20 p values) and WAT-specific mitomiRs; “PI3K-Akt signaling pathway” was the only one significantly enriched for gastrocnemius-specific mitomiRs targetome; full lists are given as Supplementary material. Biological pathways were considered statistically significant if p value was less than 0.05 (Benjamini–Hochberg procedure for multiple correction).
KEGG annotations was helpful to put mitochondrial pathways potentially governed by miRNAs in a larger context, but we were also interested in exploring mitochondria-centric pathway annotations complementary to the broader pathway database. To that aim, the targetome lists were matched with 1140 transcripts encoding proteins with mitochondrial localization in the database of MitoCarta3.0, manually curated by adding experimentally validated transcripts targets as resulting from PubMed for our detected mitomiRs (Fig 7 and S4 Table). The assignment to the “MitoPathways” spanning the 7 categories related to mitochondria showed that they are all represented for tissue-shared mitomiRs targetome, with an overrepresentation of “Metabolism”, thus indicating a miRNA-mediated regulation contributing to the key roles of mitochondria. “Metabolism” is also prevailing on the 7 categories emerging from testis-specific mitomiRs targetome and on the 5 categories represented for gastrocnemius muscle. “Mitochondrial dynamic and surveillance”, “Metabolism”, “Oxphos” and “Mitochondrial central dogma” have prevailing and similar enrichment percentage for targetome of liver-specific mitomiRs, but also all the other categories are represented. With regard to WAT-specific mitomiRs targetome, 4 categories appeared enriched, with “OXPHOS” and “Mitochondrial central dogma” as the prevailing ones. These analyses point to a role of mitomiRs in regulating the different mitochondria-related pathways, but also show gaps in the knowledge (e.g, 16% of unknown pathways for liver) that should be filled with future experimentation.
The targetomes of tissue-shared (A) and tissue-specific mitomiRs (B, gastrocnemius muscle; C, liver; D, testis; E, WAT) were further processed as described in Materials and Methods and then analyzed by MitoCarta3.0 for assignment to the seven broad functional categories relevant to mitochondria (“MitoPathways”).
Overall, targetome analyses of detected mitomiRs indicated that they can strongly support mitochondria biology and their regulatory power should be further investigated in a mitochondria-centric view, given the key role of the organelle in physiology and pathology.
Discussion
Mitochondria are “multitasked” organelles, tailoring even their morphology, number and distribution to execute specialized functions in different cell types and/or different environments, behind the oxidative phosphorylation. Proteins and protein-driven pathways involved in mitochondria homeostasis are well known, whereas non-coding RNA contribution is still overlooked. In order to gain insights into possible new and specialized function of small RNAs related to mitochondria biology, here we performed a survey of small RNA content of mitochondria purified from different tissues that particularly relies on mitochondria functioning (Fig 1).
It is described for the first time the presence of mitopiRNAs in somatic tissues, behind testis tissues, thus supporting the idea of a strong and ubiquitous connection between mitochondria biology and piRNA pathways (Figs 2 and 3). piRNAs represents an emerging and fascinating class of short RNA, longer than miRNAs, usually 24–35 nt, whose biogenesis and functioning still require further investigation along with validation of those piRNAs identified only on their size, rather than on their PIWI-interaction or functional roles. It is well established that piRNAs are enriched in germ cells, where they modulate repression of transposable elements, and probably also gene expression regulation [5–7]. With regard to biogenesis, piRNAs are generally transcribed as longer precursors from specific genomic regions termed “piRNA clusters”; then, they are processed into mature piRNAs trough two interconnected pathways, ping-pong amplification mediated by PIWI proteins in ribonucleoprotein granules, and phasing, where PIWI-loaded piRNA precursors are translocated to the mitochondrial outer membrane by RNA helicases and further cleaved by endonuclease Zucchini/PLD [75]. Also PIWIL-1 was observed to localize to mitochondria, thus raising the possibility of a mitochondrial role in the biogenesis of piRNAs [44]. piRNAs in mitochondria have been already reported, but only for mammalian cancer cells [22, 44] and inferred by bioinformatic analyses on RNA purified from the whole cell lysates deriving from murine primordial germ cells, spermatozoa, oocytes, zygotes and embryonic gonadal cells [21, 41]. Those studies also suggest the possibility that some piRNAs may be generated also from mitochondrial genome. In the mitochondria purified from the tissues, we found different populations of piRNAs, shared or tissue-specific, but none of them seems to be produced by mitochondrial genome. To assess that, all piRNA sequences contained in piRBAse were aligned to mouse mitochondrial genome, finding 7.29% of perfectly matching piRNA full sequences (S5 Table), and thus possibly deriving from mitochondrial genome, consistently with data already reported [41]; however, no reads related to those perfectly matching piRNAs were found in our samples. Different reads related to the other piRNAs, but not completely covering their full sequences, cross-map to mitochondrial genome (S6 Table and S1 Fig), probably due to the presence of NUMTs, nuclear- mitochondrial identical DNA sequences, found in many metazoans [14, 76]. Overall, our data indicate that mitopiRNAs populations found here and commonly or specifically expressed in the different tissues are generated by nuclear genome, translocated to mitochondria and involved in a piRNA-mediated communication between mitochondria and nucleus, whose mechanisms and biological meaning deserves further investigation.
With regard to miRNAs, it is known that they can regulate different pathways related to mitochondria, e.g. by targeting transcripts encoding proteins involved in TCA cycle, OXPHOS, and mitochondrial dynamics [13, 75]. Although miRNAs are largely operative in the cytoplasm, some other papers report their localization in the nucleus and mitochondria [13, 15–22], whose significance and mechanism of action are still largely unknown, since the debated presence of RISC components inside the organelle. Here we describe a mitomiRs landscape, finding a group of mitomiRs shared by all tissues analyzed, and thus potentially endowed with a regulatory role required for general mitochondria functioning, and tissue-specific mitomiRs, potentially related to specialized functions dependent on cell context (Fig 4). Among the top-10 expressed miRNAs, different ones have experimentally validated targets related to mitochondria functioning, as reported in Tables 1–5; however, for many of them studies demonstrating possible functional links with mitochondria are still missing and could be inspired by their finding inside the organelle. To gain insights into biological functions of detected mitomiRs, we extended our analyses to the whole targetomes of all tissue-shared and all tissue-specific mitomiRs, considering both experimentally validated and predicted targets and their enrichment in biological pathways. KEGG analyses highlighted that the tissue-shared mitomiRs greatly contribute to the central role of mitochondria in different signaling routes, overrepresented among the top p-values scored pathways; mitomiRs detected in specific tissues potentially contribute to regulate biological pathways typical of those tissues, as detailed in Results, overall supporting the idea of a shared population of mitomiRs for basic functioning and others for specialized functions related to the different tissues. In a more mitochondria-centric view, exploitation of MitoCarta3.0 database further confirm the presence of targets encoding proteins related to mitochondria functioning in the whole targetomes and their involvement in pathways specific of the organelle such as “Metabolism” and “OXPHOS” (Fig 7 and S4 Table).
With regard to mapping of mitomiRs to nuclear genome, intriguingly, 8 out of 10 most expressed mitomiRs specifically found in testis is produced by X chromosome (Table 5), consistently with data already reported and arguing a special role of miRNAs from X chromosome in male fertility [77, 78].
Finally, it is debated if mitomiRs could also derive from mitochondrial genome, although the mainstream view is that they originate from the nucleus and then translocated to mitochondria as illlustrated in Introduction. In order to give our contribution to the debate, we followed the same workflow described above for piRNAs to assess the miRNAs origin in our sample: all murine miRNAs deposited in miRbase were aligned to mitochondrial genome, finding only mmu-miR-6481 and mmu-miR-12206-5p sequences perfectly matching as full sequences to mitochondrial genome. However, those 2 miRNAs are not represented in our RNAseq data, instead finding reads related to other miRNAs, although not covering their full sequences (S6 Table and S1 Fig). Also taking into consideration the absence of Dicer in mitochondria [18, 47, 79], the data suggest that in our samples the reads related to miRNAs and mapping to mitochondrial genome can be attributed to miRNAs produced by their well-known chromosomal location, and could cross-map to mitochondrial genome for the presence of NUMTs, as also discussed above for mitopiRNAs [14, 76]. Besides miRNAs and piRNAs, other small non-coding RNA biotypes (generally called mitosRNA) could be produced by mitochondrial genome, not affected by Dicer absence, but probably products of currently unidentified mitochondrial ribonucleases, and whose biogenesis and role still remain unknown requiring further investigation [47].
Overall, our analyses represent a large data platform obtained by RNAseq on highly purified mitochondria from different tissues, suggesting the idea of a “basic” small RNA population residing into mitochondria and shared by different tissues, and others restricted to specific tissues. The significantly different profiles between mitochondria and cytosol fractions for some analyzed mitomiRs (Fig 5) support the relevance of miRNAs for mitochondria functioning and a potential selective translocation. The next challenge will be the planning of functional studies to deepen our understanding of small RNA processing, the shuttling among different subcellular compartments and their governed regulatory networks contributing to mitochondria biology.
Supporting information
S1 Table. mitopiRNAs detected in the current study, as annotated in RNA central (4 excel sheets).
https://doi.org/10.1371/journal.pone.0293644.s001
(XLSX)
S2 Table. Full lists of tissues-shared and tissue-specific mitomiRs detected in the current study (5 excel sheets).
https://doi.org/10.1371/journal.pone.0293644.s002
(XLSX)
S3 Table. Significantly enriched pathways according to KEGG of tissues-shared and testis-specific mitomiR targetomes (2 excel sheets).
https://doi.org/10.1371/journal.pone.0293644.s003
(XLSX)
S4 Table. mitomiR targetome enrichments for the 149 MitoPathaways, as analyzed by MitoCarta3.0 (5 excel sheets).
https://doi.org/10.1371/journal.pone.0293644.s004
(XLSX)
S5 Table. piRNAs from piRBase aligned to mitochondrial genome.
https://doi.org/10.1371/journal.pone.0293644.s005
(XLSX)
S6 Table. piRNA and miRNA related reads aligned to mitochondrial genome, found in the different tissues and graphically reported in S1 Fig (8 excel sheets).
https://doi.org/10.1371/journal.pone.0293644.s006
(XLSX)
S1 Fig. Mapping of piRNA and miRNA related reads to murine mitochondrial genome.
Starting from the inner space, the four green tracks report reads (represented by radial bars) related to piRNAs found in gastrocnemius muscle, liver, testis and WAT tissues; then, the four blue tracks report reads related to miRNAs found in the same tissues with the same order as above. The red circle represents the mitochondrial genome (NC_005089) 16299 bp long and its genes with position indicated by numbers on external circle. Full information is included in S6 Table.
https://doi.org/10.1371/journal.pone.0293644.s007
(TIFF)
References
- 1. ENCODE Project Consortium. An integrated encyclopedia of DNA elements in the human genome. Nature. 2012;489: 57–74. pmid:22955616
- 2. Beermann J, Piccoli M-T, Viereck J, Thum T. Non-coding RNAs in Development and Disease: Background, Mechanisms, and Therapeutic Approaches. Physiol Rev. 2016;96: 1297–1325. pmid:27535639
- 3. Weick E-M, Miska EA. piRNAs: from biogenesis to function. Development. 2014;141: 3458–3471. pmid:25183868
- 4. Ku H-Y, Lin H. PIWI proteins and their interactors in piRNA biogenesis, germline development and gene expression. Natl Sci Rev. 2014;1: 205–218. pmid:25512877
- 5. Aravin AA, Naumova NM, Tulin AV, Vagin VV, Rozovsky YM, Gvozdev VA. Double-stranded RNA-mediated silencing of genomic tandem repeats and transposable elements in the D. melanogaster germline. Curr Biol. 2001;11: 1017–1027. pmid:11470406
- 6. Sun YH, Lee B, Li XZ. The birth of piRNAs: how mammalian piRNAs are produced, originated, and evolved. Mamm Genome. 2022;33: 293–311. pmid:34724117
- 7. Zhang J, Chen S, Liu K. Structural insights into piRNA biogenesis. Biochim Biophys Acta Gene Regul Mech. 2022;1865: 194799. pmid:35182819
- 8. Bartel DP. MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004;116: 281–297. pmid:14744438
- 9. Dragomir MP, Knutsen E, Calin GA. Classical and noncanonical functions of miRNAs in cancers. Trends Genet. 2022;38: 379–394. pmid:34728089
- 10. Chan JJ, Tay Y. Noncoding RNA:RNA Regulatory Networks in Cancer. Int J Mol Sci. 2018;19: E1310. pmid:29702599
- 11. Di Palo A, Siniscalchi C, Mosca N, Russo A, Potenza N. A Novel ceRNA Regulatory Network Involving the Long Non-Coding Antisense RNA SPACA6P-AS, miR-125a and its mRNA Targets in Hepatocarcinoma Cells. Int J Mol Sci. 2020;21: E5068. pmid:32709089
- 12. Leung AKL. The Whereabouts of microRNA Actions: Cytoplasm and Beyond. Trends Cell Biol. 2015;25: 601–610. pmid:26410406
- 13. Gusic M, Prokisch H. ncRNAs: New Players in Mitochondrial Health and Disease? Front Genet. 2020;11: 95. pmid:32180794
- 14. Ren B, Guan M-X, Zhou T, Cai X, Shan G. Emerging functions of mitochondria-encoded noncoding RNAs. Trends Genet. 2023;39: 125–139. pmid:36137834
- 15. Lung B, Zemann A, Madej MJ, Schuelke M, Techritz S, Ruf S, et al. Identification of small non-coding RNAs from mitochondria and chloroplasts. Nucleic Acids Res. 2006;34: 3842–3852. pmid:16899451
- 16. Kren BT, Wong PY-P, Sarver A, Zhang X, Zeng Y, Steer CJ. MicroRNAs identified in highly purified liver-derived mitochondria may play a role in apoptosis. RNA Biol. 2009;6: 65–72. pmid:19106625
- 17. Bian Z, Li L-M, Tang R, Hou D-X, Chen X, Zhang C-Y, et al. Identification of mouse liver mitochondria-associated miRNAs and their potential biological functions. Cell Res. 2010;20: 1076–1078. pmid:20733615
- 18. Bandiera S, Matégot R, Girard M, Demongeot J, Henrion-Caude A. MitomiRs delineating the intracellular localization of microRNAs at mitochondria. Free Radic Biol Med. 2013;64: 12–19. pmid:23792138
- 19. Bandiera S, Rüberg S, Girard M, Cagnard N, Hanein S, Chrétien D, et al. Nuclear outsourcing of RNA interference components to human mitochondria. PLoS One. 2011;6: e20746. pmid:21695135
- 20. Barrey E, Saint-Auret G, Bonnamy B, Damas D, Boyer O, Gidrol X. Pre-microRNA and mature microRNA in human mitochondria. PLoS One. 2011;6: e20220. pmid:21637849
- 21. Larriba E, Rial E, Del Mazo J. The landscape of mitochondrial small non-coding RNAs in the PGCs of male mice, spermatogonia, gametes and in zygotes. BMC Genomics. 2018;19: 634. pmid:30153810
- 22. Sripada L, Tomar D, Prajapati P, Singh R, Singh AK, Singh R. Systematic analysis of small RNAs associated with human mitochondria by deep sequencing: detailed analysis of mitochondrial associated miRNA. PLoS One. 2012;7: e44873. pmid:22984580
- 23. Spinelli JB, Haigis MC. The multifaceted contributions of mitochondria to cellular metabolism. Nat Cell Biol. 2018;20: 745–754. pmid:29950572
- 24. Tuppen HAL, Blakely EL, Turnbull DM, Taylor RW. Mitochondrial DNA mutations and human disease. Biochim Biophys Acta. 2010;1797: 113–128. pmid:19761752
- 25. Wang C, Liu X, Wei B. Mitochondrion: an emerging platform critical for host antiviral signaling. Expert Opin Ther Targets. 2011;15: 647–665. pmid:21476879
- 26. Kamp DW, Shacter E, Weitzman SA. Chronic inflammation and cancer: the role of the mitochondria. Oncology (Williston Park). 2011;25: 400–410, 413. pmid:21710835
- 27. Schon EA, Przedborski S. Mitochondria: the next (neurode)generation. Neuron. 2011;70: 1033–1053. pmid:21689593
- 28. Boore JL. Animal mitochondrial genomes. Nucleic Acids Res. 1999;27: 1767–1780. pmid:10101183
- 29. Salinas T, Duchêne A-M, Maréchal-Drouard L. Recent advances in tRNA mitochondrial import. Trends Biochem Sci. 2008;33: 320–329. pmid:18513973
- 30. Giordani C, Silvestrini A, Giuliani A, Olivieri F, Rippo MR. MicroRNAs as Factors in Bidirectional Crosstalk Between Mitochondria and the Nucleus During Cellular Senescence. Front Physiol. 2021;12: 734976. pmid:34566699
- 31. Geiger J, Dalgaard LT. Interplay of mitochondrial metabolism and microRNAs. Cell Mol Life Sci. 2017;74: 631–646. pmid:27563705
- 32. Aparicio-Puerta E, Lebrón R, Rueda A, Gómez-Martín C, Giannoukakos S, Jaspez D, et al. sRNAbench and sRNAtoolbox 2019: intuitive fast small RNA profiling and differential expression. Nucleic Acids Res. 2019;47: W530–W535. pmid:31114926
- 33. Barreñada O, Larriba E, Brieño-Enriquez MA, Mazo JD. piRNA-IPdb: a PIWI-bound piRNAs database to mining NGS sncRNA data and beyond. BMC Genomics. 2021;22: 765. pmid:34702185
- 34. Huang DW, Sherman BT, Lempicki RA. Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2009;37: 1–13. pmid:19033363
- 35. Rath S, Sharma R, Gupta R, Ast T, Chan C, Durham TJ, et al. MitoCarta3.0: an updated mitochondrial proteome now with sub-organelle localization and pathway annotations. Nucleic Acids Res. 2021;49: D1541–D1547. pmid:33174596
- 36. An P, Wei L-L, Zhao S, Sverdlov DY, Vaid KA, Miyamoto M, et al. Hepatocyte mitochondria-derived danger signals directly activate hepatic stellate cells and drive progression of liver fibrosis. Nat Commun. 2020;11: 2362. pmid:32398673
- 37. Carter HN, Chen CCW, Hood DA. Mitochondria, muscle health, and exercise with advancing age. Physiology (Bethesda). 2015;30: 208–223. pmid:25933821
- 38. Park S-Y, Gifford JR, Andtbacka RHI, Trinity JD, Hyngstrom JR, Garten RS, et al. Cardiac, skeletal, and smooth muscle mitochondrial respiration: are all mitochondria created equal? Am J Physiol Heart Circ Physiol. 2014;307: H346–352. pmid:24906913
- 39. Wang X, Yin L, Wen Y, Yuan S. Mitochondrial regulation during male germ cell development. Cell Mol Life Sci. 2022;79: 91. pmid:35072818
- 40. Zhu Q, An YA, Scherer PE. Mitochondrial regulation and white adipose tissue homeostasis. Trends Cell Biol. 2022;32: 351–364. pmid:34810062
- 41. Barreñada O, Larriba E, Fernández-Pérez D, Brieño-Enríquez MÁ, Del Mazo Martínez J. Unraveling mitochondrial piRNAs in mouse embryonic gonadal cells. Sci Rep. 2022;12: 10730. pmid:35750721
- 42. Yan Z, Hu HY, Jiang X, Maierhofer V, Neb E, He L, et al. Widespread expression of piRNA-like molecules in somatic tissues. Nucleic Acids Res. 2011;39: 6596–6607. pmid:21546553
- 43. Ross RJ, Weiner MM, Lin H. PIWI proteins and PIWI-interacting RNAs in the soma. Nature. 2014;505: 353–359. pmid:24429634
- 44. Kwon C, Tak H, Rho M, Chang HR, Kim YH, Kim KT, et al. Detection of PIWI and piRNAs in the mitochondria of mammalian cancer cells. Biochem Biophys Res Commun. 2014;446: 218–223. pmid:24602614
- 45. Watanabe T, Chuma S, Yamamoto Y, Kuramochi-Miyagawa S, Totoki Y, Toyoda A, et al. MITOPLD is a mitochondrial protein essential for nuage formation and piRNA biogenesis in the mouse germline. Dev Cell. 2011;20: 364–375. pmid:21397847
- 46. Huang H, Gao Q, Peng X, Choi S-Y, Sarma K, Ren H, et al. piRNA-associated germline nuage formation and spermatogenesis require MitoPLD profusogenic mitochondrial-surface lipid signaling. Dev Cell. 2011;20: 376–387. pmid:21397848
- 47. Ro S, Ma H-Y, Park C, Ortogero N, Song R, Hennig GW, et al. The mitochondrial genome encodes abundant small noncoding RNAs. Cell Res. 2013;23: 759–774. pmid:23478297
- 48. Landgraf P, Rusu M, Sheridan R, Sewer A, Iovino N, Aravin A, et al. A mammalian microRNA expression atlas based on small RNA library sequencing. Cell. 2007;129: 1401–1414. pmid:17604727
- 49. Isakova A, Fehlmann T, Keller A, Quake SR. A mouse tissue atlas of small noncoding RNA. Proc Natl Acad Sci U S A. 2020;117: 25634–25645. pmid:32978296
- 50. Shi Y, Zhang Z, Yin Q, Fu C, Barszczyk A, Zhang X, et al. Cardiac-specific overexpression of miR-122 induces mitochondria-dependent cardiomyocyte apoptosis and promotes heart failure by inhibiting Hand2. J Cell Mol Med. 2021;25: 5326–5334. pmid:33942477
- 51. Wang Y, Jin P, Liu J, Xie X. Exosomal microRNA-122 mediates obesity-related cardiomyopathy through suppressing mitochondrial ADP-ribosylation factor-like 2. Clin Sci (Lond). 2019;133: 1871–1881. pmid:31434696
- 52. He R, Ding C, Yin P, He L, Xu Q, Wu Z, et al. MiR-1a-3p mitigates isoproterenol-induced heart failure by enhancing the expression of mitochondrial ND1 and COX1. Exp Cell Res. 2019;378: 87–97. pmid:30853447
- 53. Fuseini H, Cephus J-Y, Wu P, Davis JB, Contreras DC, Gandhi VD, et al. ERα Signaling Increased IL-17A Production in Th17 Cells by Upregulating IL-23R Expression, Mitochondrial Respiration, and Proliferation. Front Immunol. 2019;10: 2740. pmid:31849948
- 54. Xu G, Li JY. ATP5A1 and ATP5B are highly expressed in glioblastoma tumor cells and endothelial cells of microvascular proliferation. J Neurooncol. 2016;126: 405–413. pmid:26526033
- 55. Park S, Lee M, Chun C-H, Jin E-J. The lncRNA, Nespas, Is Associated with Osteoarthritis Progression and Serves as a Potential New Prognostic Biomarker. Cartilage. 2019;10: 148–156. pmid:28805067
- 56. Sharma P, Sharma V, Ahluwalia TS, Dogra N, Kumar S, Singh S. Let-7a induces metabolic reprogramming in breast cancer cells via targeting mitochondrial encoded ND4. Cancer Cell Int. 2021;21: 629. pmid:34838007
- 57. Lu C-H, Chen D-X, Dong K, Wu Y-J, Na N, Wen H, et al. Inhibition of miR-143-3p alleviates myocardial ischemia reperfusion injury via limiting mitochondria-mediated apoptosis. Biol Chem. 2023. pmid:36780323
- 58. Huang S, Wu K, Li B, Liu Y. lncRNA UCA1 inhibits mitochondrial dysfunction of skeletal muscle in type 2 diabetes mellitus by sequestering miR-143-3p to release FGF21. Cell Tissue Res. 2023;391: 561–575. pmid:36602629
- 59. Jung SE, Kim SW, Jeong S, Moon H, Choi WS, Lim S, et al. MicroRNA-26a/b-5p promotes myocardial infarction-induced cell death by downregulating cytochrome c oxidase 5a. Exp Mol Med. 2021;53: 1332–1343. pmid:34518647
- 60. Cheng Y-C, Su L-Y, Chen L-H, Lu T-P, Chuang EY, Tsai M-H, et al. Regulatory Mechanisms and Functional Roles of Hypoxia-Induced Long Non-Coding RNA MTORT1 in Breast Cancer Cells. Front Oncol. 2021;11: 663114. pmid:34141617
- 61. Araujo HN, Lima TI, Guimarães DSPSF, Oliveira AG, Favero-Santos BC, Branco RCS, et al. Regulation of Lin28a-miRNA let-7b-5p pathway in skeletal muscle cells by peroxisome proliferator-activated receptor delta. Am J Physiol Cell Physiol. 2020;319: C541–C551. pmid:32697599
- 62. Li H, Dai B, Fan J, Chen C, Nie X, Yin Z, et al. The Different Roles of miRNA-92a-2-5p and let-7b-5p in Mitochondrial Translation in db/db Mice. Mol Ther Nucleic Acids. 2019;17: 424–435. pmid:31319246
- 63. Cagnin S, Brugnaro M, Millino C, Pacchioni B, Troiano C, Di Sante M, et al. Monoamine Oxidase-Dependent Pro-Survival Signaling in Diabetic Hearts Is Mediated by miRNAs. Cells. 2022;11: 2697. pmid:36078109
- 64. Canale P, Nicolini G, Pitto L, Kusmic C, Rizzo M, Balzan S, et al. Role of miR-133/Dio3 Axis in the T3-Dependent Modulation of Cardiac mitoK-ATP Expression. Int J Mol Sci. 2022;23: 6549. pmid:35742991
- 65. Kim S, Park JW, Lee MG, Nam KH, Park JH, Oh H, et al. Reversine promotes browning of white adipocytes by suppressing miR-133a. J Cell Physiol. 2019;234: 3800–3813. pmid:30132867
- 66. Willers IM, Martínez-Reyes I, Martínez-Diez M, Cuezva JM. miR-127-5p targets the 3’UTR of human β-F1-ATPase mRNA and inhibits its translation. Biochim Biophys Acta. 2012;1817: 838–848. pmid:22433606
- 67. Zhao Y, Xie Y, Yao W-Y, Wang Y-Y, Song N. Long non-coding RNA Opa interacting protein 5-antisense RNA 1 promotes mitochondrial autophagy and protects SH-SY5Y cells from 1-methyl-4-phenylpyridine-induced damage by binding to microRNA-137 and upregulating NIX. Kaohsiung J Med Sci. 2022;38: 207–217. pmid:35049152
- 68. Channakkar AS, Singh T, Pattnaik B, Gupta K, Seth P, Adlakha YK. MiRNA-137-mediated modulation of mitochondrial dynamics regulates human neural stem cell fate. Stem Cells. 2020;38: 683–697. pmid:32012382
- 69. Li W, Zhang X, Zhuang H, Chen H-G, Chen Y, Tian W, et al. MicroRNA-137 is a novel hypoxia-responsive microRNA that inhibits mitophagy via regulation of two mitophagy receptors FUNDC1 and NIX. J Biol Chem. 2014;289: 10691–10701. pmid:24573672
- 70. Zhang C, Ma C, Zhang L, Zhang L, Zhang F, Ma M, et al. MiR-449a-5p mediates mitochondrial dysfunction and phenotypic transition by targeting Myc in pulmonary arterial smooth muscle cells. J Mol Med (Berl). 2019;97: 409–422. pmid:30715622
- 71. Zou T, Zhang J, Liu Y, Zhang Y, Sugimoto K, Mei C. Antidepressant-Like Effect of Geniposide in Mice Exposed to a Chronic Mild Stress Involves the microRNA-298-5p-Mediated Nox1. Front Mol Neurosci. 2020;13: 131. pmid:33613190
- 72. Satoh T, Wang L, Espinosa-Diez C, Wang B, Hahn SA, Noda K, et al. Metabolic Syndrome Mediates ROS-miR-193b-NFYA-Dependent Downregulation of Soluble Guanylate Cyclase and Contributes to Exercise-Induced Pulmonary Hypertension in Heart Failure With Preserved Ejection Fraction. Circulation. 2021;144: 615–637. pmid:34157861
- 73. Yang F, Li T, Dong Z, Mi R. MicroRNA-410 is involved in mitophagy after cardiac ischemia/reperfusion injury by targeting high-mobility group box 1 protein. J Cell Biochem. 2018;119: 2427–2439. pmid:28914970
- 74. Manfrevola F, Potenza N, Chioccarelli T, Di Palo A, Siniscalchi C, Porreca V, et al. Actin remodeling driven by circLIMA1: sperm cell as an intriguing cellular model. Int J Biol Sci. 2022;18: 5136–5153. pmid:35982890
- 75. Sun W, Lu Y, Zhang H, Zhang J, Fang X, Wang J, et al. Mitochondrial Non-Coding RNAs Are Potential Mediators of Mitochondrial Homeostasis. Biomolecules. 2022;12: 1863. pmid:36551291
- 76. Pozzi A, Dowling DK. The Genomic Origins of Small Mitochondrial RNAs: Are They Transcribed by the Mitochondrial DNA or by Mitochondrial Pseudogenes within the Nucleus (NUMTs)? Genome Biol Evol. 2019;11: 1883–1896. pmid:31218347
- 77. Di Palo A, Siniscalchi C, Salerno M, Russo A, Gravholt CH, Potenza N. What microRNAs could tell us about the human X chromosome. Cell Mol Life Sci. 2020;77: 4069–4080. pmid:32356180
- 78. Guo X, Su B, Zhou Z, Sha J. Rapid evolution of mammalian X-linked testis microRNAs. BMC Genomics. 2009;10: 97. pmid:19257908
- 79. Das S, Ferlito M, Kent OA, Fox-Talbot K, Wang R, Liu D, et al. Nuclear miRNA regulates the mitochondrial genome in the heart. Circ Res. 2012;110: 1596–1603. pmid:22518031