Xu Liu
Ge Shan- Hefei National Laboratory for Physical Sciences at Microscale, The CAS Key Laboratory of Innate Immunity and Chronic Disease, Division of Life Science and Medicine, Department of Clinical Laboratory, The First Affiliated Hospital of USTC, School of Basic Medical Sciences, University of Science and Technology of China, Hefei, China
Mitochondria are the powerhouses of mammalian cells, which participate in series of metabolic processes and cellular events. Mitochondria have their own genomes, and it is generally acknowledged that human mitochondrial genome encodes 13 proteins, 2 rRNAs and 22 tRNAs. However, the complexity of mitochondria derived transcripts is just starting to be envisaged. Currently, there are at least 8 lncRNAs, some dsRNAs, various small RNAs, and hundreds of circRNAs known to be generated from mitochondrial genome. These non-coding RNAs either translocate into cytosol/nucleus or reside in mitochondria to play various biological functions. Here we present an overview of regulatory non-coding RNAs encoded by the mammalian mitochondria genome. For overall understandings of non-coding RNAs in mitochondrial function, a brief summarization of nuclear-encoded non-coding RNAs in mitochondria is also included. We discuss about roles of these non-coding RNAs in cellular physiology and the communication between mitochondria and the nucleus.
Introduction
Mitochondria are vital to cells and are involved in multiple essential cellular events. On top of generating ATP, mitochondria also produce reactive oxygen species (ROS), redox molecules, and intermediates required for the synthesis of biomolecules. Furthermore, they are key hubs of intracellular signaling pathways and participate in response to external environmental changes (McBride et al., 2006).
Mitochondria have independent circular genomes (mtDNA) (Figure 1). Human mtDNA is mostly coding (more than 90%) and transcribed entirely in a bidirectional, polycistronic manner. Both strands of mtDNA harbor coding information. The heavy strand encodes 12 protein genes which are all subunits of the oxidative phosphorylation system (OXPHOS), two ribosomal RNAs (12S and 16S), and 14 tRNAs. The light strand encodes for a single mRNA (ND6) and 8 tRNAs (Scarpulla, 2008; Gustafsson et al., 2016; Barshad et al., 2018). Transcription is initiated from two heavy-strand promoters (HSP1, HSP2) and one light-strand promoter (LSP), located in the major non-coding region named “D-loop”. Two long polycistronic RNAs produced from HSP2 and LSP cover almost the full-length of each strand and another transcript transcribed from HSP1 consists of the two rRNA genes (Taanman, 1999; Asin-Cayuela and Gustafsson, 2007). The long mitochondrial precursor transcripts undergo the generally accepted “tRNA punctuations” processing by RNase P and RNase Z (ELAC2) cleavage to generate individual mRNAs, rRNAs, and tRNAs (Ojala et al., 1981; Holzmann et al., 2008; Lopez Sanchez et al., 2011; Barchiesi and Vascotto, 2019).
Figure 1. Map of the human mitochondrial genome with the following features: Genes corresponding to 13 mRNAs, 2 rRNAs, and 22 tRNAs within heavy (outer track) and light (inner track) strands; D-loop region and promoters (HSP1, HSP2, LSP) of both strands; eight mitochondria-encoded lncRNAs shown in the inner side of each ring (lncND5, lncCyt b, lncND6, MDL1S, MDL1AS, SmtncRNA, ASmtncRNA, and LIPCAR); four recently reported functional circRNAs (mecciND1, mecciND5, mc-COX2, circRNA SCAR).
Recent years have witnessed the fast growing in understanding the great complexity of transcripts from the nuclear genome (Djebali et al., 2012; Iyer et al., 2015), and growing evidence has also demonstrated the presence of mitochondria-encoded non-coding RNAs (ncRNAs) such as long non-coding RNA (lncRNA) (Jusic and Devaux, 2020), circular RNAs (circRNAs) (Liu et al., 2020; Wu et al., 2020; Zhao et al., 2020), small non-coding RNAs (sncRNAs) (Mercer et al., 2011; Jusic and Devaux, 2020), and double-stranded RNAs (dsRNAs) (Dhir et al., 2018) with potential regulatory functions, respectively (Figure 1).
In this review we will focus on known regulatory non-coding RNAs encoded by the mitochondrial genome, and we also present a summarization of nuclear-encoded non-coding RNAs found in mitochondria for overall understandings about ncRNAs in mitochondria-nucleus communication. We highlight and discuss about their functional relevance in biomedicine, and in the context of cross-talk between mitochondria and the nucleus.
Mitochondria-Encoded lncRNAs
LncRNA is a major class of non-coding RNAs that is widely expressed in cells. Nuclear-encoded lncRNAs play diverse regulatory roles in mammalian cells, including transcriptional regulation, translational regulation, protein scaffolding, chromosome remodeling etc. (Gomes et al., 2013; Rashid et al., 2016; Marchese et al., 2017; Chen and Shan, 2020). A series of lncRNAs transcribed from the mtDNA have been recently reported, and these lncRNAs may also participate in multiple biological processes (Figure 2 and Table 1).
Figure 2. Proposed functions and functional mechanisms of mitochondria-encoded noncoding RNAs.
Table 1. Mitochondria-encode non-coding RNAs.
Three lncRNAs (lncND5, lncCyt b, lncND6) generated from the mitochondrial genome were identified by Rackham et al. (2011) from deep sequencing data of human HeLa mitochondria (Rackham et al., 2011). The positions of these three lncRNAs corresponded to the regions complementary to the mitochondrial ND5, ND6, and Cyt b genes. There was no significant open reading frame present in them. It was shown that the three lncRNAs might play a functional role to stabilize ND5, ND6, and Cyt b mRNAs or to regulate the expression of the three genes in mitochondria through forming RNA-RNA duplexes (Rackham et al., 2011). The abundance of the three lncRNAs was cell and tissue-specific, indicating that they might be subjected to cell specific regulation and to play physiological roles.
Villegas et al. (2007) reported a 2,374 nt human mitochondrial transcript containing an inverted repeat (IR) of 820 nt linked to the 5′ end of the 16S mt-rRNA and named it SncmtRNA (sense non-coding mitochondrial RNA) (Villegas et al., 2007). In SncmtRNA, the IR generates a stem-loop structure with a double stranded region and a 40 nt loop. SncmtRNA was over-expressed in several tumor cell lines but not in resting cells (Villegas et al., 2007). The same researcher team later identified 2 transcripts containing IRs linked to the 5′ region of the antisense transcript of 16S mt-rRNA gene in the L-strand of the mtDNA, and designated them antisense ncmtRNA-1 (ASncmtRNA-1, 310 nt IR) and antisense ncmtRNA-2 (ASncmtRNA-2, 545 nt IR) (Burzio et al., 2009). AsncmtRNAs were expressed in normal proliferating cells but were down-regulated in different types of human tumor cells (Burzio et al., 2009). The subcellular location of SncmtRNA and ASncmtRNAs seemed to be nucleus (Landerer et al., 2011). Knockdown of ASncmtRNAs (ASK for short) by antisense oligonucleotide (ASO-1537S) targeting the single-stranded loop region resulted in Dicer-mediated release of hsa-miR-4485. hsa-miR-4485 in combination with nuclear miRNAs (mainly hsa-miR-5096 and hsa-miR-3609), which were induced by ASK, inhibited translation of 5 important cell cycle proteins (Cyclin B1, Cyclin D1, CDK1, CDK4, Survivin), and thereby induced growth inhibition of breast cancer cells (Fitzpatrick et al., 2019). The same researcher group later reported that exosomes derived from ASO-1537S-treated MDA-MB-231 breast cancer cells inhibited tumorigenesis of recipient cells (Lobos-González et al., 2020). And they also observed ASK induced bladder cancer cell death and inhibition of tumor growth (Borgna et al., 2020). The above results indicated that ASncmtRNAs could be potential therapeutic targets for breast cancer and bladder cancer.
Global transcriptomic profiling of human left ventricular (LV) from heart failing patients and controls revealed that in LV, relatively high abundance (71%) of lncRNAs were encoded by mitochondrial genome (mito-lncRNAs), and the abundance of mito-lncRNAs showed some but not statistically significant reduction in both non-ischemic and ischemic human failing hearts (Yang et al., 2014). A mitochondria-encoded chimeric transcript, long intergenic non-coding RNA predicting cardiac remodeling (LIPCAR) with the 5′-half (1–392 nt) mapped to lncCyt b and 3′-half (385–781 nt) mapped to antisense of mt-COX2 gene was identified in plasma from myocardial infarction patients (Kumarswamy et al., 2014). The expression level of LIPCAR was downregulated at the early stage of myocardial infarction but upregulated during later stages; in addition, high LIPCAR levels were associated with future cardiovascular death. The circulating LIPCAR has a potential to be a risk-free biomarker of LV remodeling and a predictor of patient survival of heart failure (Kumarswamy et al., 2014).
Using a PacBio full-length third-generation sequencing transcriptome dataset, Gao et al. (2018) discovered two polycistronic transcripts hsa-MDL1 (Mitochondrial D-loop 1) and hsa-MDL1AS (Mitochondrial D-loop 1 antisense) generated from the region covering the tRNAPro gene and full length of the human D-loop region (Gao et al., 2018). The presence of these transcripts was further confirmed by pan RNA-seq analysis (Xu et al., 2019). The mature hsa-MDL1 was generated from the H-strand polycistronic transcript, and hsa-MDLAS was antisense to hsa-MDL1 with much lower expression level. A great number of small RNAs reversely aligned to the mitochondrial D-loop region and a relatively lower abundant forwardly aligned small RNA (5% of reversely aligned ones) were obtained by searching public small RNA-sequencing data. hsa-MDL1 and hsa-MDL1AS were believed to be precursors of some of these small RNAs (Gao et al., 2018). Pan RNA-seq analysis also revealed the ubiquitous existence of 5′ and 3′ end small RNAs of MDL1 and MDL1AS (Xu et al., 2019).
Mitochondria-Encoded dsRNAs
Dhir et al. (2018) discovered highly unstable mitochondrial double-stranded RNA (mt-dsRNA) in HeLa cells, which was supposed to be the natural outcomes of the nearly complete transcription of both heavy and light strands of mtDNA. It has been known that almost the entire L-strand transcript undergoes rapid RNA turnover by the mitochondrial RNA degradosome (Borowski et al., 2013), and indeed the level of mt-dsRNAs was restricted by the RNA degradosome components SUV3 and PNPT1 (Dhir et al., 2018). Pathological PNPT1 mutations led to abnormal mt-dsRNA accumulation. The Bax-Bak dependent release of mt-dsRNAs into cytoplasm triggered the upregulation of interferon-stimulated genes and the activation of innate immune defenses through the MDA5–MAVS axis (Dhir et al., 2018). In Trp53 mutant mouse embryonic fibroblasts, mitochondrial dsRNAs cleaved by RNase L were immunogenic, and could activate the type I interferon (IFN) pathway via RIG-I-like receptors (Wiatrek et al., 2019).
Mitochondria-Encoded circRNAs
circRNAs encoded by mammalian mitochondrial genome had not been revealed until recent report about the identification of hundreds of mitochondria encoded circRNAs (mecciRNAs) in human and murine cells through second-generation sequencing of mitochondrial RNAs (Liu et al., 2020, 2021). In three cell lines and one tumor tissue (hepatocellular carcinoma) from human and two cell lines together with two tissues (skeletal muscle and heart) from mice, 248 human mecciRNAs and 268 murine mecciRNAs were identified, respectively (Liu et al., 2020). Actually, the presence of large amount of mecciRNAs was not mammalian specific, as over one hundred mecciRNAs were also identified in zebrafish (Liu et al., 2020). mecciRNAs were encoded by both the light and heavy strands of the mtDNA, although the heavy strand of mtDNA encoded the majority (Liu et al., 2020). The presence of mitochondria encoded circRNAs in human HEK293 cells was examined by RT-PCR assay to amplify over one hundred individual mecciRNAs (Mance et al., 2020). With third-generation sequencing, mecciRNAs from murine brain were also analyzed for their junctions and full length (Zhang et al., 2021).
It was found that some mecciRNAs distributed both inside the mitochondria and outside in the cytosol (Liu et al., 2020). Functional studies of two mecciRNAs, mecciND1 (encoded by the mitochondrial ND1 gene) and mecciND5 (encoded by the mitochondrial ND5 gene) demonstrated that these mecciRNAs conducted essential physiological functions. mecciND1 bound to RPA1 and RPA2 proteins, which were involved in mtDNA replication. The expression level of mecciND1 was positively related to mitochondrial RPA protein levels and mtDNA copy numbers. mtDNA copy numbers were increased after mecciND1 overexpression and decreased under mecciND1 knockdown. mecciND5 interacted with hnRNPA1, hnRNPA2B1 and hnRNPA3 and promoted their mitochondrial importation (Liu et al., 2020). Both in vivo and in vitro evidence showed that the two mecciRNAs, mecciND1 and mecciND5, could interact with TOM40 and PNPASE, and function as molecular chaperons to facilitate the mitochondrial entry of newly synthesized polypeptides encoded by the nuclear genome (Liu et al., 2020; Figure 2). mecciND1 and meccciND5 were upregulated in hepatocellular carcinoma (HCC) tissues, and moreover, mecciRNA levels were regulated in stress conditions, suggesting that they were critical for cells to cope with physiological and pathological changes (Liu et al., 2020).
When investigating the role of mitochondria-located circRNAs in metaflammation, Zhao et al. (2020) observed that 3 mitochondria-encoded circRNAs were downregulated in liver fibroblasts from patients with non-alcoholic steatohepatitis (NASH). One of the three mecciRNAs, termed SCAR (Steatohepatitis-associated circRNA ATP5B Regulator), was an antisense RNA from the locus COX2. SCAR bound directly to ATP5B, a mitochondrial permeability transition pore (mPTP) regulator. The interaction of ATP5B and SCAR shut down mPTP by blocking Cyclophilin D-mPTP interaction, and therefore, inhibited mitochondrial ROS (mROS) output (Zhao et al., 2020; Figure 2). mROS output was required for liver fibroblast activation, and this circRNA had a potential to serve as a therapeutic target for NASH (Zhao et al., 2020). Another highly expressed mecciRNA, mc-COX2, a sense RNA from the locus COX2, was found in the plasma exosomes of chronic lymphocytic leukemia (CLL) patients (Wu et al., 2020). mc-COX2 was closely correlated to prognosis of CLL, and moreover, it seemed that higher expression levels of mc-COX2 could promote cell proliferation and protect cells from apoptosis (Wu et al., 2020; Figure 2).
Mitochondria-Encoded sncRNAs
In mammalian cells, small non-coding RNAs (sncRNAs) are diverse and abundant, mainly including microRNAs (miRNAs), endogenous-short interfering RNA (siRNAs), PIWI-interacting RNAs (piRNAs), and other types of small RNAs derived from tRNAs, rRNAs, and snoRNAs (Gomes et al., 2013). Recent studies have identified several types of sncRNAs that are generated from mammalian mitochondrial genome (Table 1).
Mercer et al. (2011) provided a comprehensive description of the human mitochondrial transcriptome. They sequenced mitochondrial small RNAs (sRNAs) isolated from human bone osteosarcoma cell line 143B and revealed 31 sRNAs of two distinct classes of 21 and 26 nt from 17 loci of mtDNA. The majority (84%) of them were produced from tRNA genes. Expression of these mitochondria-encoded sncRNAs changed dynamically in different cell types, however, the expression of sncRNAs did not seem to correlate significantly with the expression of overlapping genes (Mercer et al., 2011). Ro et al. (2013) reported that thousands of small RNAs were encoded from murine and human mitochondrial genomes, and major size of these sRNAs ranged between 30 and 39 nt (mouse) and 20 and 29 nt (human). These RNAs were mainly transcribed from the sense direction of the mitochondrial genes (host genes), and only a small portion were from the antisense direction of their host genes. Those sRNAs appeared to target antisense transcripts and promoted the expression of their host genes in vitro (Ro et al., 2013). Interestingly, all these researches provided evidence that these mitochondria-encoded sncRNAs were generated in a Dicer-independent manner, and might be products of some unknown ribonucleases within the mitochondria (Mercer et al., 2011; Ro et al., 2013).
In human HEK293 and HeLa mitochondria, 4 known miRNAs (hsa-miR-4461, hsa-miR-4463, hsa-miR-4484, and hsa-miR-4485) and 24 putative novel miRNAs could be aligned to mitochondrial genome at the positions corresponding to 16S rRNA, tRNA, and mRNA (Sripada et al., 2012). These miRNAs might regulate site-specific turnover of target mRNAs. However, it was not conclusive that whether these microRNAs were actually transcribed from mitochondrial genome (Sripada et al., 2012). A series of small RNAs that aligned to the D-loop region of mtDNA were discovered by Gao et al. (2018). These small RNAs might be generated from mitochondrial lncRNAs hsa-MDL1AS and hsa-MDL1. The most abundant small RNA was hsa-tir-MDL1AS-18, derived from lncRNA hsa-MDL1AS. It belonged to transcription initial RNAs (tiRNAs), as 18 nucleotides from its 5′ end precisely overlapped with transcription initiation (IT) sites of the L-strand promoter. hsa-tir-MDL1AS-18 was down-regulated in hepatocellular carcinoma (HCC) tissues, suggesting that the balance of tiRNAs/lncRNAs regulation might be abnormal in cancer cells (Gao et al., 2018).
Through mapping known piRNA sequences to the human mtDNA, Kwon et al. (2014) identified 29 piRNAs, and 12 out of these 29 piRNAs matched to the stem-loop fragments of seven mitochondrial tRNAs with asymmetric tRNA fragment usage and cell type specific expression. They also reported the presence of PIWI through Western blots and three abundant mature piRNAs by RT-PCR in the mitochondria of HeLa-S3 cells (Kwon et al., 2014). Mitochondria-encoded piRNAs might also be present in mouse (Larriba et al., 2018). Among all the small RNA reads from oocytes and zygotes of mice, about 20% of the reads were mapped to mitochondrial DNA, and 80–90% of the reads were classified as piRNAs based on the size; these piRNAs were encoded by both strands of the mouse mitochondrial genome, and were highly expressed from the loci corresponding to mitochondrial tRNAs, 16S rRNA, and D-loop region (Larriba et al., 2018). Functional study of mitochondrial piRNAs is lacking, and these piRNAs are speculated to function in stress responses (Kwon et al., 2014), gamete differentiation, and fertilization (Larriba et al., 2018).
Nuclear-Encoded ncRNAs in Mitochondria
In mammalian mitochondria, besides some well-known nuclear encoded tRNAs (tRNALeuUAA, tRNAGlnUUG, and tRNAGlnCUG) (Rubio et al., 2008; Mercer et al., 2011; Gowher et al., 2013), 5S rRNA (Yoshionari et al., 1994; Magalhães et al., 1998; Entelis et al., 2001; Mercer et al., 2011; Zelenka et al., 2014; Autour et al., 2018), RMRP RNA (Chang and Clayton, 1987; Li et al., 1994; Wang et al., 2010; Noh et al., 2016), some nuclear-encoded ncRNAs including miRNAs and lncRNAs were also found to play critical roles in mitochondria (Kim et al., 2017; Jeandard et al., 2019; Table 2). miRNAs were described in purified mitochondria of rat and mouse through miRNA microarray analysis (Kren et al., 2009; Bian et al., 2010). Bandiera et al. (2011) identified 13 nuclear-encoded miRNAs enriched in mitochondria of HeLa cells, and defined these mitochondria-located miRNAs as mitomiRs. Meanwhile, Barrey et al. (2011) predicted 25 pre-miRNAs and 33 miRNAs in human mitochondrial in silico sequencing analysis, and a set of them were demonstrated to be present in the mitochondria of skeletal muscular cells by in situ hybridization and RT-PCR (Barrey et al., 2011). Moreover, both studies revealed that these miRNAs might post-transcriptionally regulate mt-RNAs in mitochondria through RNA interfering pathway (Bandiera et al., 2011; Barrey et al., 2011). Nuclear DNA-encoded miR-18c was found to be translocated into the mitochondria of rat cardiac myocytes (Das et al., 2012). miR-18c loaded in AGO2n bound to the 3′ UTR of mt-COX1, and inhibited mt-COX1 translation, ultimately leading to complex IV remodeling and mitochondrial dysfunction (Das et al., 2012). mitomiR-378 conducted a similar function that downregulated mt-ATP6 translation in mouse HL-1 cells (Jagannathan et al., 2015). Fan and co-workers identified mitomiR-2392 was upregulated in cisplatin-resistant tongue squamous cell carcinoma (TSCC) cells and TSCC tumors (Fan et al., 2019). Surprisingly, mitomiR-2392 was found to partially repress the transcription of mtDNA rather than the translation, through miRNA:mtDNA base pairing with an AGO2-dependent manner (Fan et al., 2019). The downstream genes regulated by mitomiR-2392 were different between two TSCC cell lines, suggesting a cell-specific regulation (Fan et al., 2019). In addition to inhibitory effects, some mitomiRs were found playing positive regulatory functions. For example, miR-1 was induced and imported into mitochondria during myogenesis to stimulate the translation of mt-ND1 and mt-COX1, through specific miRNA:mRNA base-pairing with the participation of AGO2 in mouse heart and C2C12 cells (Zhang et al., 2014). Most mitomiRs with demonstrated functions played AGO2-dependent roles, and an array of evidence showed that AGO2 was present in mitochondria (Bandiera et al., 2011, 2013; Zhang et al., 2014; Vendramin et al., 2017).
Table 2. Nuclear-encoded non-coding RNAs in mitochondria.
The malignant lncRNA SAMMSON (survival associated mitochondrial melanoma specific oncogenic non-coding RNA) was predominantly localized in the cytoplasm of human melanoblasts and melanoma cells, and interestingly, a large fraction of cytoplasmic SAMMSON was found to co-localize and co-purify with mitochondria (Leucci et al., 2016). SAMMSON played pro-oncogenic roles in mitochondrial homeostasis and metabolism by interacting with and stabilizing p32 protein, a regulator of the maturation of mitochondrial 16S rRNA (Leucci et al., 2016; Vendramin et al., 2018). Human nuclear RNase P RNA component (H1 RNA) was reported in mitochondria of HeLa cells (Bartkiewicz et al., 1989; Puranam and Attardi, 2001). However, as human mitochondrial RNase P was a protein-only complex and a trans-acting RNA component was not required for catalysis (Holzmann et al., 2008), the role of H1 RNA in mitochondria was still unclear. Cheng et al. (2018) found that human telomerase RNA component (hTERC) was imported into mitochondria, processed to a shorter form hTERC-53, and then exported back to the cytosol (Cheng et al., 2018). The underlaying mechanism and cellular significance of this cytosol and mitochondria import-export process of hTERC-hTERC-53 is elusive and needs further investigation.
Discussion
It is apparent that the known complexity of mitochondrial transcriptome of human and mammals has been greatly expanded. Mitochondrial genome encodes not only very limited numbers of mRNAs, rRNAs and tRNAs as previously understood, but also a variety of non-coding RNAs such as lncRNAs, circRNAs, sncRNAs, and dsRNAs with diverse regulatory functions. Mitochondria seem to have sacrificed the protein coding ability, although maintained the capacity of generating their “own” complex ncRNA profiles, with the small size but highly utilized genome.
Mitochondria-encoded lncRNAs seem function through RNA-RNA interaction or as precursors of small RNAs to regulate the stabilization or expression of their corresponding mt-RNAs in mitochondria (Rackham et al., 2011; Gao et al., 2018). Several mt-lncRNAs can be detected outside the mitochondria, such as LIPCAR in plasma, SncmtRNA and ASncmtRNAs in nucleus, but the mechanisms of their transport or function are still elusive (Landerer et al., 2011; Kumarswamy et al., 2014). circRNAs encoded by the mammalian nuclear genome are known to function as microRNA sponges, transcriptional regulators, protein binding partners, and templates for protein translation (Chen et al., 2015; Li et al., 2015; Kristensen et al., 2019; Chen and Shan, 2021). Mitochondria-encoded circRNAs may reside in or shuttle in and out of the mitochondria to play critical roles, and the functions and functional mechanisms of mecciRNAs are waiting for further in-depth investigation (Liu et al., 2020; Zhao et al., 2020). The sRNAs, miRNAs, piRNAs, and other snRNAs encoded by mitochondrial genome are generally classified by their size, and the functions of these non-coding RNAs need further elucidation. Some nuclear genome encoded miRNAs (mitomiRs) can also function in mitochondria. Studies of these mitomiRs may provide some implications for future investigation on mitochondria-encoded sncRNAs. The presence of AGO2 in mitochondria makes it reasonable to speculate that certain sncRNAs can also function through RNAi machinery in mitochondria (Bandiera et al., 2011; Bandiera et al., 2013; Zhang et al., 2014; Vendramin et al., 2017). The functional characterizations of mitochondria-encoded ncRNAs are just emerging, and the biogenesis and metabolism of mitochondria-encoded ncRNAs, including circRNAs, lncRNAs, and sncRNAs, remain largely obscure.
Due to the limited coding capacity of mitochondrial genome, it is known that a large number of proteins encoded by the nuclear genome are imported into mitochondria from cytosol to sustain the biogenesis and function of these organelles (Harbauer et al., 2014; Wiedemann and Pfanner, 2017). Mitochondria generate energy, ROS, and other metabolic molecules for multiple essential cellular events, and also participate in multiple pathways regulating nuclear gene expression, cell death, proliferation, differentiation, etc. (McBride et al., 2006). The growing evidence and examples reviewed in this article have shown that non-coding RNAs are engaged in the bidirectional communication between mitochondria and the nucleus (Vendramin et al., 2017; Gusic and Prokisch, 2020). Nuclear DNA-encoded miRNAs and lncRNAs can also be transported into mitochondria and regulate mitochondrial gene expression transcriptionally and post-transcriptionally. However, the import mechanisms are still poorly understood. Several lines of evidence show that tRNAs are imported into mitochondria by an ATP dependent mechanism that is distinct from protein import (Rubio et al., 2008; Gowher et al., 2013). Disruption of the mitochondrial membrane potential, which is crucial for protein import, has no effect on tRNA import (Rubio et al., 2008). There is only one protein PNPT1 (also named PNPASE), which locates in the intermembrane space (IMS) of mitochondria, known to facilitate the transport of RMRP, H1 RNA, and hTERC through binding to small stem-loops of these RNAs (Wang et al., 2010; Kim et al., 2017; Jeandard et al., 2019). Mitochondria-encoded ncRNAs contribute to the nucleus-mitochondria communication, and some of them not only function inside the mitochondria but also play roles outside of the mitochondria. For example, mecciRNAs shuttle between mitochondria and cytosol to facilitate the mitochondrial entry of newly synthesized polypeptides encoded by the nuclear genome (Liu et al., 2020). ASncmtRNAs are reported with nuclear and cytosolic localization, and regulate the translation of cell cycle proteins (Landerer et al., 2011; Fitzpatrick et al., 2019). But how these ncRNAs get transported through the bilayer membrane of mitochondria is still unclear, and further investigations about the mechanisms will contribute to a deeper understanding of mitochondria-nucleus communication.
Mitochondrial dysfunction relates to a series of diseases, such as cardiovascular diseases, cancers, and neurodegeneration (Herst et al., 2017; Jusic and Devaux, 2020). The cell/tissue specific expression of mt-ncRNAs also suggests that they may take part in mitochondria-related diseases. The investigations of mitochondrial ncRNAs regulatory network would contribute to better understanding of etiology, and lead to novel diagnostic and therapeutic approaches for mitochondrial dysfunction-related diseases. Mitochondria-encoded non-coding RNAs have no doubt pointed to a new direction for the study of mitochondria-nucleus communication, and study about these RNAs is an important field of biomedicine.
Author Contributions
GS and XL conceived the scope of the manuscript and wrote the manuscript. Both authors contributed to the article and approved the submitted version.
Funding
This work was supported by grants from the National Key R&D Program of China (2019YFA0802600 and 2018YFC1004500), the National Natural Science Foundation of China (32000438, 31725016, 31930019, and 91940303), and the Strategic Priority Research Program “Biological basis of aging and therapeutic strategies” of the Chinese Academy of Sciences (XDB39010400).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
Asin-Cayuela, J., and Gustafsson, C. M. (2007). Mitochondrial transcription and its regulation in mammalian cells. Trends Biochem. Sci. 32, 111–117. doi: 10.1016/j.tibs.2007.01.003
Autour, A., C Y Jeng, S., D Cawte, A., Abdolahzadeh, A., Galli, A., Panchapakesan, S. S. S., et al. (2018). Fluorogenic RNA Mango aptamers for imaging small non-coding RNAs in mammalian cells. Nat. Commun. 9:656. doi: 10.1038/s41467-018-02993-8
Bandiera, S., Matégot, R., Girard, M., Demongeot, J., and Henrion-Caude, A. (2013). MitomiRs delineating the intracellular localization of microRNAs at mitochondria. Free Radic. Biol. Med. 64, 12–19. doi: 10.1016/j.freeradbiomed.2013.06.013
Bandiera, S., Rüberg, S., Girard, M., Cagnard, N., Hanein, S., Chrétien, D., et al. (2011). Nuclear outsourcing of RNA interference components to human mitochondria. PLoS One 6:e20746. doi: 10.1371/journal.pone.0020746
Barchiesi, A., and Vascotto, C. (2019). Transcription, processing, and decay of mitochondrial RNA in health and disease. Int. J. Mol. Sci. 20:2221. doi: 10.3390/ijms20092221
Barrey, E., Saint-Auret, G., Bonnamy, B., Damas, D., Boyer, O., and Gidrol, X. (2011). Pre-microRNA and mature microRNA in human mitochondria. PLoS One 6:e20220. doi: 10.1371/journal.pone.0020220
Barshad, G., Marom, S., Cohen, T., and Mishmar, D. (2018). Mitochondrial DNA transcription and its regulation: an evolutionary perspective. Trends Genet. 34, 682–692. doi: 10.1016/j.tig.2018.05.009
Bartkiewicz, M., Gold, H., and Altman, S. (1989). Identification and characterization of an RNA molecule that copurifies with RNase P activity from HeLa cells. Genes Dev. 3, 488–499. doi: 10.1101/gad.3.4.488
Bian, Z., Li, L., Tang, R., Hou, D., Chen, X., Zhang, C., et al. (2010). Identification of mouse liver mitochondria-associated miRNAs and their potential biological functions. Cell Res. 20, 1076–1078. doi: 10.1038/cr.2010.119
Borgna, V., Lobos-González, L., Guevara, F., Landerer, E., Bendek, M., Ávila, R., et al. (2020). Targeting antisense mitochondrial noncoding RNAs induces bladder cancer cell death and inhibition of tumor growth through reduction of survival and invasion factors. J. Cancer 11, 1780–1791. doi: 10.7150/jca.38880
Borowski, L. S., Dziembowski, A., Hejnowicz, M. S., Stepien, P. P., and Szczesny, R. J. (2013). Human mitochondrial RNA decay mediated by PNPase–hSuv3 complex takes place in distinct foci. Nucleic Acids Res. 41, 1223–1240. doi: 10.1093/nar/gks1130
Burzio, V. A., Villota, C., Villegas, J., Landerer, E., Boccardo, E., Villa, L. L., et al. (2009). Expression of a family of noncoding mitochondrial RNAs distinguishes normal from cancer cells. Proc. Natl. Acad. Sci. U.S.A. 106, 9430–9434. doi: 10.1073/pnas.0903086106
Chang, D. D., and Clayton, D. A. (1987). A mammalian mitochondrial RNA processing activity contains nucleus-encoded RNA. Science 235, 1178–1184. doi: 10.1126/science.2434997
Chen, H., and Shan, G. (2020). The physiological function of long-noncoding RNAs. Noncoding RNA Res. 5, 178–184. doi: 10.1016/j.ncrna.2020.09.003
Chen, L., and Shan, G. (2021). CircRNA in cancer: fundamental mechanism and clinical potential. Cancer Lett. 505, 49–57. doi: 10.1016/j.canlet.2021.02.004
Chen, L., Huang, C., Wang, X. L., and Shan, G. (2015). Circular RNAs in eukaryotic cells. Curr. Genomics 16, 312–318. doi: 10.2174/1389202916666150707161554
Cheng, Y., Liu, P. P., Zheng, Q., Gao, G., Yuan, J. P., Wang, P. F., et al. (2018). Mitochondrial trafficking and processing of telomerase RNA TERC. Cell Rep. 24, 2589–2595. doi: 10.1016/j.celrep.2018.08.003
Das, S., Ferlito, M., Kent, O. A., Fox-Talbot, K., Wang, R., Liu, D., et al. (2012). Nuclear miRNA regulates the mitochondrial genome in the heart. Circ. Res. 110, 1596–1603. doi: 10.1161/circresaha.112.267732
Dhir, A., Dhir, S., Borowski, L. S., Jimenez, L., Teitell, M., Rötig, A., et al. (2018). Mitochondrial double-stranded RNA triggers antiviral signalling in humans. Nature 560, 238–242. doi: 10.1038/s41586-018-0363-0
Djebali, S., Davis, C. A., Merkel, A., Dobin, A., Lassmann, T., Mortazavi, A., et al. (2012). Landscape of transcription in human cells. Nature 489, 101–108. doi: 10.1038/nature11233
Entelis, N. S., Kolesnikova, O. A., Dogan, S., Martin, R. P., and Tarassov, I. A. (2001). 5 S rRNA and tRNA import into human mitochondria. comparison of in vitro requirements. J. Biol. Chem. 276, 45642–45653. doi: 10.1074/jbc.M103906200
Esteller, M. (2011). Non-coding RNAs in human disease. Nat. Rev. Genet. 12, 861–874. doi: 10.1038/nrg3074
Fan, S., Tian, T., Chen, W., Lv, X., Lei, X., Zhang, H., et al. (2019). Mitochondrial miRNA determines chemoresistance by reprogramming metabolism and regulating mitochondrial transcription. Cancer Res. 79, 1069–1084. doi: 10.1158/0008-5472.can-18-2505
Fitzpatrick, C., Bendek, M. F., Briones, M., Farfán, N., Silva, V. A., Nardocci, G., et al. (2019). 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. 10:423. doi: 10.1038/s41419-019-1649-3
Gao, S., Tian, X., Chang, H., Sun, Y., Wu, Z., Cheng, Z., et al. (2018). Two novel lncRNAs discovered in human mitochondrial DNA using PacBio full-length transcriptome data. Mitochondrion 38, 41–47. doi: 10.1016/j.mito.2017.08.002
Gomes, A. Q., Nolasco, S., and Soares, H. (2013). Non-coding RNAs: multi-tasking molecules in the cell. Int. J. Mol. Sci. 14, 16010–16039. doi: 10.3390/ijms140816010
Gowher, A., Smirnov, A., Tarassov, I., and Entelis, N. (2013). Induced tRNA import into human mitochondria: implication of a host aminoacyl-tRNA-synthetase. PLoS One 8:e66228. doi: 10.1371/journal.pone.0066228
Gusic, M., and Prokisch, H. (2020). ncRNAs: new players in mitochondrial health and disease? Front. Genet. 11:95. doi: 10.3389/fgene.2020.00095
Gustafsson, C. M., Falkenberg, M., and Larsson, N. G. (2016). Maintenance and expression of mammalian mitochondrial DNA. Annu. Rev. Biochem. 85, 133–160. doi: 10.1146/annurev-biochem-060815-014402
Harbauer, A. B., Zahedi, R. P., Sickmann, A., Pfanner, N., and Meisinger, C. (2014). The protein import machinery of mitochondria-a regulatory hub in metabolism, stress, and disease. Cell Metab. 19, 357–372. doi: 10.1016/j.cmet.2014.01.010
Herst, P. M., Rowe, M. R., Carson, G. M., and Berridge, M. V. (2017). Functional mitochondria in health and disease. Front. Endocrinol. 8:296. doi: 10.3389/fendo.2017.00296
Holzmann, J., Frank, P., Löffler, E., Bennett, K. L., Gerner, C., and Rossmanith, W. (2008). RNase P without RNA: identification and functional reconstitution of the human mitochondrial tRNA processing enzyme. Cell 135, 462–474. doi: 10.1016/j.cell.2008.09.013
Iyer, M. K., Niknafs, Y. S., Malik, R., Singhal, U., Sahu, A., Hosono, Y., et al. (2015). The landscape of long noncoding RNAs in the human transcriptome. Nat. Genet. 47, 199–208. doi: 10.1038/ng.3192
Jagannathan, R., Thapa, D., Nichols, C. E., Shepherd, D. L., Stricker, J. C., Croston, T. L., et al. (2015). Translational regulation of the mitochondrial genome following redistribution of mitochondrial MicroRNA in the diabetic heart. Circ. Cardiovasc. Genet. 8, 785–802. doi: 10.1161/circgenetics.115.001067
Jeandard, D., Smirnova, A., Tarassov, I., Barrey, E., Smirnov, A., and Entelis, N. (2019). Import of non-coding RNAs into human mitochondria: a critical review and emerging approaches. Cells 8:286. doi: 10.3390/cells8030286
Jusic, A., and Devaux, Y. (2020). Mitochondrial noncoding RNA-regulatory network in cardiovascular disease. Basic. Res. Cardiol. 115:23. doi: 10.1007/s00395-020-0783-5
Kim, K. M., Noh, J. H., Abdelmohsen, K., and Gorospe, M. (2017). Mitochondrial noncoding RNA transport. BMB Rep. 50, 164–174. doi: 10.5483/bmbrep.2017.50.4.013
Kren, B. T., Wong, P. Y. P., Sarver, A., Zhang, X., Zeng, Y., and Steer, C. J. (2009). MicroRNAs identified in highly purified liver-derived mitochondria may play a role in apoptosis. RNA Biol. 6, 65–72. doi: 10.4161/rna.6.1.7534
Kristensen, L. S., Andersen, M. S., Stagsted, L. V., Ebbesen, K. K., Hansen, T. B., and Kjems, J. (2019). The biogenesis, biology and characterization of circular RNAs. Nat. Rev. Genet. 20, 675–691. doi: 10.1038/s41576-019-0158-7
Kumarswamy, R., Bauters, C., Volkmann, I., Maury, F., Fetisch, J., Holzmann, A., et al. (2014). Circulating long noncoding RNA, LIPCAR, predicts survival in patients with heart failure. Circ. Res. 114, 1569–1575. doi: 10.1161/CIRCRESAHA.114.303915
Kwon, C., Tak, H., Rho, M., Chang, H. R., Kim, Y. H., Kim, K., et al. (2014). Detection of PIWI and piRNAs in the mitochondria of mammalian cancer cells. Biochem. Biophys. Res. Commun. 446, 218–223. doi: 10.1016/j.bbrc.2014.02.112
Landerer, E., Villegas, J., Burzio, V. A., Oliveira, L., Villota, C., Lopez, C., et al. (2011). Nuclear localization of the mitochondrial ncRNAs in normal and cancer cells. Cell. Oncol. 34, 297–305. doi: 10.1007/s13402-011-0018-8
Larriba, E., Rial, E., and Del Mazo, J. (2018). The landscape of mitochondrial small non-coding RNAs in the PGCs of male mice, spermatogonia, gametes and in zygotes. BMC Genomics 19:634. doi: 10.1186/s12864-018-5020-3
Leucci, E., Vendramin, R., Spinazzi, M., Laurette, P., Fiers, M., Wouters, J., et al. (2016). Melanoma addiction to the long non-coding RNA SAMMSON. Nature 531, 518–522. doi: 10.1038/nature17161
Li, K., Smagula, C. S., Parsons, W. J., Richardson, J. A., Gonzalez, M., Hagler, H. K., et al. (1994). Subcellular partitioning of MRP RNA assessed by ultrastructural and biochemical analysis. J. Cell Biol. 124, 871–882. doi: 10.1083/jcb.124.6.871
Li, Z., Huang, C., Bao, C., Chen, L., Lin, M., Wang, X., et al. (2015). Exon-intron circular RNAs regulate transcription in the nucleus. Nat. Struct. Mol. Biol. 22, 256–264. doi: 10.1038/nsmb.2959
Liu, X., Wang, X., Li, J., Hu, S., Deng, Y., Yin, H., et al. (2020). Identification of mecciRNAs and their roles in the mitochondrial entry of proteins. Sci. China Life Sci. 63, 1429–1449. doi: 10.1007/s11427-020-1631-9
Liu, X., Yang, Y., and Shan, G. (2021). Identification and detection of mecciRNAs. Methods 21, 43–48. doi: 10.1016/j.ymeth.2021.02.006
Lobos-González, L., Bustos, R., Campos, A., Silva, V., Silva, V., Jeldes, E., et al. (2020). Exosomes released upon mitochondrial ASncmtRNA knockdown reduce tumorigenic properties of malignant breast cancer cells. Sci. Rep. 10:343. doi: 10.1038/s41598-019-57018-1
Lopez Sanchez, M. I., Mercer, T. R., Davies, S. M., Shearwood, A.-M. J., Nygård, K. K., Richman, T. R., et al. (2011). RNA processing in human mitochondria. Cell Cycle 10, 2904–2916. doi: 10.4161/cc.10.17.17060
Magalhães, P. J., Andreu, A. L., and Schon, E. A. (1998). Evidence for the presence of 5S rRNA in mammalian mitochondria. Mol. Biol. Cell. 9, 2375–2382. doi: 10.1091/mbc.9.9.2375
Mance, L. G., Mawla, I., Shell, S. M., and Cahoon, A. B. (2020). Mitochondrial mRNA fragments are circularized in a human HEK cell line. Mitochondrion 51, 1–6. doi: 10.1016/j.mito.2019.11.002
Marchese, F. P., Raimondi, I., and Huarte, M. (2017). The multidimensional mechanisms of long noncoding RNA function. Genome Biol. 18:206. doi: 10.1186/s13059-017-1348-2
McBride, H. M., Neuspiel, M., and Wasiak, S. (2006). Mitochondria: more than just a powerhouse. Curr. Biol. 16, R551–R560. doi: 10.1016/j.cub.2006.06.054
Mercer, T. R., Neph, S., Dinger, M. E., Crawford, J., Smith, M. A., Shearwood, A.-M. J., et al. (2011). The human mitochondrial transcriptome. Cell 146, 645–658. doi: 10.1016/j.cell.2011.06.051
Noh, J. H., Kim, K. M., Abdelmohsen, K., Yoon, J. H., Panda, A. C., Munk, R., et al. (2016). HuR and GRSF1 modulate the nuclear export and mitochondrial localization of the lncRNA RMRP. Genes Dev. 30, 1224–1239. doi: 10.1101/gad.276022.115
Ojala, D., Montoya, J., and Attardi, G. (1981). tRNA punctuation model of RNA processing in human mitochondria. Nature 290, 470–474. doi: 10.1038/290470a0
Puranam, R. S., and Attardi, G. (2001). The RNase P associated with HeLa cell mitochondria contains an essential RNA component identical in sequence to that of the nuclear RNase P. Mol. Cell Biol. 21, 548–561. doi: 10.1128/mcb.21.2.548-561.2001
Rackham, O., Shearwood, A. M., Mercer, T. R., Davies, S. M., Mattick, J. S., and Filipovska, A. (2011). Long noncoding RNAs are generated from the mitochondrial genome and regulated by nuclear-encoded proteins. RNA 17, 2085–2093. doi: 10.1261/rna.029405.111
Rashid, F., Shah, A., and Shan, G. (2016). Long non-coding RNAs in the cytoplasm. Genomics Proteomics Bioinformatics 14, 73–80. doi: 10.1016/j.gpb.2016.03.005
Ro, S., Ma, H. Y., Park, C., Ortogero, N., Song, R., Hennig, G. W., et al. (2013). The mitochondrial genome encodes abundant small noncoding RNAs. Cell Res. 23, 759–774. doi: 10.1038/cr.2013.37
Rubio, M. A., Rinehart, J. J., Krett, B., Duvezin-Caubet, S., Reichert, A. S., Söll, D., et al. (2008). Mammalian mitochondria have the innate ability to import tRNAs by a mechanism distinct from protein import. Proc. Natl. Acad. Sci. U.S.A. 105, 9186–9191. doi: 10.1073/pnas.0804283105
Scarpulla, R. C. (2008). Transcriptional paradigms in mammalian mitochondrial biogenesis and function. Physiol. Rev. 88, 611–638. doi: 10.1152/physrev.00025.2007
Sripada, L., Tomar, D., Prajapati, P., Singh, R., Singh, A. K., and Singh, R. (2012). Systematic analysis of small RNAs associated with human mitochondria by deep sequencing: detailed analysis of mitochondrial associated miRNA. PLoS One 7:e44873. doi: 10.1371/journal.pone.0044873
Taanman, J.-W. (1999). The mitochondrial genome: structure, transcription, translation and replication. Biochim. Biophys. Acta. 1410, 103–123. doi: 10.1016/s0005-2728(98)00161-3
Vendramin, R., Marine, J. C., and Leucci, E. (2017). Non-coding RNAs: the dark side of nuclear-mitochondrial communication. EMBO J. 36, 1123–1133. doi: 10.15252/embj.201695546
Vendramin, R., Verheyden, Y., Ishikawa, H., Goedert, L., Nicolas, E., Saraf, K., et al. (2018). SAMMSON fosters cancer cell fitness by concertedly enhancing mitochondrial and cytosolic translation. Nat. Struct. Mol. Biol. 25, 1035–1046. doi: 10.1038/s41594-018-0143-4
Villegas, J., Burzio, V., Villota, C., Landerer, E., Martinez, R., Santander, M., et al. (2007). Expression of a novel non-coding mitochondrial RNA in human proliferating cells. Nucleic Acids Res. 35, 7336–7347. doi: 10.1093/nar/gkm863
Wang, G., Chen, H. W., Oktay, Y., Zhang, J., Allen, E. L., Smith, G. M., et al. (2010). PNPASE regulates RNA import into mitochondria. Cell 142, 456–467. doi: 10.1016/j.cell.2010.06.035
Wiatrek, D. M., Candela, M. E., Sedmík, J., Oppelt, J., Keegan, L. P., and O’Connell, M. A. (2019). Activation of innate immunity by mitochondrial dsRNA in mouse cells lacking p53 protein. RNA 25, 713–726. doi: 10.1261/rna.069625.118
Wiedemann, N., and Pfanner, N. (2017). Mitochondrial machineries for protein import and assembly. Annu. Rev. Biochem. 86, 685–714. doi: 10.1146/annurev-biochem-060815-014352
Wu, Z., Sun, H., Wang, C., Liu, W., Liu, M., Zhu, Y., et al. (2020). Mitochondrial genome-derived circRNA mc-COX2 functions as an oncogene in chronic lymphocytic leukemia. Mol. Ther. Nucleic Acids 20, 801–811. doi: 10.1016/j.omtn.2020.04.017
Xu, X., Ji, H., Jin, X., Cheng, Z., Yao, X., Liu, Y., et al. (2019). Using Pan RNA-Seq Analysis to Reveal the Ubiquitous Existence of 5′ and 3′ End Small RNAs. Front. Genet. 10:105. doi: 10.3389/fgene.2019.00105
Yang, K. C., Yamada, K. A., Patel, A. Y., Topkara, V. K., George, I., Cheema, F. H., et al. (2014). Deep RNA sequencing reveals dynamic regulation of myocardial noncoding RNAs in failing human heart and remodeling with mechanical circulatory support. Circulation 129, 1009–1021. doi: 10.1161/CIRCULATIONAHA.113.003863
Yoshionari, S., Koike, T., Yokogawa, T., Nishikawa, K., Ueda, T., Miura, K., et al. (1994). Existence of nuclear-encoded 5S-rRNA in bovine mitochondria. FEBS Lett. 338, 137–142. doi: 10.1016/0014-5793(94)80351-x
Zelenka, J., Alán, L., Jabůrek, M., and Ježek, P. (2014). Import of desired nucleic acid sequences using addressing motif of mitochondrial ribosomal 5S-rRNA for fluorescent in vivo hybridization of mitochondrial DNA and RNA. J. Bioenerg. Biomembr. 46, 147–156. doi: 10.1007/s10863-014-9543-2
Zhang, J., Hou, L., Zuo, Z., Ji, P., Zhang, X., Xue, Y., et al. (2021). Comprehensive profiling of circular RNAs with nanopore sequencing and CIRI-long. Nat. Biotechnol. 39, 836–845. doi: 10.1038/s41587-021-00842-6
Zhang, X., Zuo, X., Yang, B., Li, Z., Xue, Y., Zhou, Y., et al. (2014). MicroRNA directly enhances mitochondrial translation during muscle differentiation. Cell 158, 607–619. doi: 10.1016/j.cell.2014.05.047
Keywords: mitochondria, mitochondria-encoded non-coding RNA, lncRNA, dsRNA, small ncRNA, circRNA
Citation: Liu X and Shan G (2021) Mitochondria Encoded Non-coding RNAs in Cell Physiology. Front. Cell Dev. Biol. 9:713729. doi: 10.3389/fcell.2021.713729
Received: 24 May 2021; Accepted: 12 July 2021;
Published: 30 July 2021.
Edited by:
Veronica Andrea Burzio, Andres Bello University, ChileReviewed by:
Amaresh Chandra Panda, Institute of Life Sciences (ILS), IndiaChandrasekhar Kanduri, University of Gothenburg, Sweden
Piotr Stepien, University of Warsaw, Poland
Copyright © 2021 Liu and Shan. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Ge Shan, shange@ustc.edu.cn