Jun Wu1,2†
Hongyu Sun- 1Department of General Surgery and Pancreatic Injury and Repair Key Laboratory of Sichuan Province, The General Hospital of Western Theater Command, Chengdu, China
- 2College of Medicine, Southwest Jiaotong University, Chengdu, China
- 3Laboratory of Basic Medicine, The General Hospital of Western Theater Command, Chengdu, China
As the most abundant internal modification in eukaryotic cells, N6-methyladenosine (m6A) in mRNA has shown widespread regulatory roles in a variety of physiological processes and disease progressions. Circular RNAs (circRNAs) are a class of covalently closed circular RNA molecules and play an essential role in the pathogenesis of various diseases. Recently, accumulating evidence has shown that m6A modification is widely existed in circRNAs and found its key biological functions in regulating circRNA metabolism, including biogenesis, translation, degradation and cellular localization. Through regulating circRNAs, studies have shown the important roles of m6A modification in circRNAs during immunity and multiple diseases, which represents a new layer of control in physiological processes and disease progressions. In this review, we focused on the roles played by m6A in circRNA metabolism, summarized the regulatory mechanisms of m6A-modified circRNAs in immunity and diseases, and discussed the current challenges to study m6A modification in circRNAs and the possible future directions, providing a comprehensive insight into understanding m6A modification of circRNAs in RNA epigenetics.
Introduction
Epigenetic modifications in RNA have been found to playimportant regulatory roles in a variety of physiological processes and disease progressions. To date, over 160 types of chemical modifications have been identified in RNA molecules, such as N6-methyladenosine (m6A), 5-methylcytosine (m5C), N1-methyladenosine (m1A), 5-hydroxymethylcytosine (5hmC), N6, 2′-Odimethyladenosine (m6Am), 7-methylguanine (m7G), and so on (Roundtree et al., 2017a; Frye et al., 2018). Among them, m6A has been considered as the most abundant internal modification in RNA, which was discovered from methylated nucleosides in mRNA of Novikoff hepatoma cells in the early 1970s (Desrosiers et al., 1974; Roundtree et al., 2017a). With the development of next-generation sequencing and bioinformatics, a large number of m6A modifications in the transcriptome were discovered and annotated, including mRNAs and non-coding RNAs (Zhou et al., 2017; Frye et al., 2018; He et al., 2020; Sekar and Lakshmanan, 2020; Zhu et al., 2020). In mRNA, the m6A modification has been extensively studied and reviewed by several excellent reviews (Roundtree et al., 2017a; Zaccara et al., 2019; Zhu et al., 2020; Jiang et al., 2021).
Up to now, accumulating evidence has demonstrated that circular RNAs (circRNAs) play important roles in the occurrence, development and prognosis of various diseases, such as cardiovascular diseases (Aufiero et al., 2019; Zhang and Huang, 2020), immune diseases (Chen X. et al., 2019), tumors (Meng S. et al., 2017), skin diseases (Wu et al., 2020), and so on. However, the regulatory role of m6A in circRNAs remains unclear. It was not until 2017 that circRNAs were found to be modified by m6A (Zhou et al., 2017). Based on the important roles of m6A in RNA epigenetic modifications, predictably and rapidly, a large number of recent studies have demonstrated that m6A acts as a key regulator to affect circRNAs functions, thereby participating in development and disease progressions. Therefore, in this review, we focused on the roles of m6A modification in circRNAs during circRNA metabolism and discussed their biological consequences in human development and diseases based on current studies. Importantly, this review will provide a comprehensive understanding of m6A in circRNA biology.
Biogenesis, Characteristics and Functions of circRNAS
circRNAs are a class of covalently closed circular RNA molecules thatlack 5′ caps and 3′ tails, which were discovered in the 1970s (Hsu and Coca-Prados, 1979; Chen and Yang, 2015). In a long time afterward, they were thought to be byproducts of splicing (Sanger et al., 1976; Panda et al., 2017). In recent years, a large number of circRNAs have been discovered and annotated with the development of next-generation sequencing technologies and bioinformatic approaches (Xia et al., 2017). The circRNAs are mainly generated by the back-splicing of pre-mRNAs. Currently, three mechanisms were proposed regarding the biogenesis of circRNAs, namely lariat-driven circularization (Zaphiropoulos, 1996; Jeck et al., 2013), intron pairing-driven circularization (Jeck et al., 2013; Zhang et al., 2014; Ivanov et al., 2015) and RBP-driven circularization (Conn et al., 2015; Kramer et al., 2015). In addition, circRNAs are also generated by pre-transfer ribonucleic acid (tRNA), and pre-ribosomal ribonucleic acid (rRNA) molecules (Zhang et al., 2017). According to their origin, the formed circRNAs can be classified into four categories (Zhang et al., 2017): intronic circRNAs (CiRNAs), exon-intron circRNAs (EIciRNAs), exonic circRNAs (EcRNAs), and others (such as the circRNAs that are derived from tRNA and rRNA).
circRNAs were expressed widely in tissue, cell, and developmentalstage-specific manners (Memczak et al., 2013; Nicolet et al., 2018) and existedin most mammals (Memczak et al., 2013), plants (Wang et al., 2014), and virus (Nahand et al., 2020). For example, in human heart, about 9% of the expressed genes can produce circRNAs (Aufiero et al., 2018), while in brain, 20% of the genes can produce circRNAs (Rybak-Wolf et al., 2015). In addition, different circular junctions can be generated from the same gene (Nicolet et al., 2018). These discoveries indicated the important regulatory roles of circRNAs. Another characteristic is that circRNAs are more stable than the linear RNAs and have long half-lives than 48 h (Jeck et al., 2013). Due to the covalently closed circular structure without 5′ caps and 3′ tails, circRNAs can resist the foreign chemicals or exonucleases (Jeck et al., 2013; Aufiero et al., 2019). Therefore, despite the low expression level, circRNAs can be accumulated to the relatively high levels.
A large number of studies have revealed the functions of circRNAs. (1) microRNA (miRNA) sponges. Cytoplasmic circRNAs can bind specific miRNAs by miRNA response elements to prevent the interplays with target mRNAs (Thomson and Dinger, 2016; Cheng et al., 2019). (2) Interaction with proteins. Some circRNAs that contain protein binding sites have the protein binding abilities to regulate their activity and localization (Du W. W. et al., 2016; Barbagallo et al., 2019). (3) Templates for protein synthesis. Due to the lack of 5′-end, the translation of circRNAs only is initiated by the cap-independent mechanisms. This function of circRNAs will be discussed in detail in following section. (4) Regulation of gene transcription. Some circRNAs located in nucleus, such as circ-EIF3J and circ-PAIP2, can combine with the U1 snRNP to further interact with RNA Pol II and enhance the expression of their parental genes in HeLa and HEK293 cells (Li et al., 2015).
m6A Modification of circRNA
Widely Existed m6A Modification of circRNAs
Through depletion of rRNA, treatment with RNase R to digest linear RNAs and m6A-modified RNA immunoprecipitation sequencing (MeRIP-seq), m6A modification of circRNAs has been found to exist widely. In plant, Wang et al. (2020) found that about 10% of the EcRNAs contained m6A sites in moso bamboo. In human embryonic stem cells (hESCs) and HeLa cells, Zhou et al. (2017) identified 1,404 m6A circRNAs and 987 m6A circRNAs, respectively, and they found that m6A circRNAs were expressed in cell-type-specific methylation patterns. In a rat model of hypoxia mediated pulmonary hypertension (HPH), Su et al. (2020) identified 1130 m6A circRNAs in total. In the lens epithelium cells from age-related cataract (ARC), Li et al. (2020) identified 4876 m6A peaks within circRNAs. Besides, in HPH, 80% of m6A circRNAs were derived from protein-coding genes, while in ARC, over 70% of circRNAs were derived from sense overlapping. These studies suggested that m6A modification of circRNAs was widely existed and might play important roles in diseases. In addition, the m6A sites in circRNAs may be distinct from the corresponding mRNAs. In hESCs cells, Zhou et al. (2017) reported that 33% of m6A circRNAs were produced from genes that encode m6A mRNAs methylated on different exons, while 26% of m6A circRNAs were produced from genes that encode mRNAs without detectable m6A modification (Zhou et al., 2017). Another phenomenon is that m6A sites in mRNAs are most common in the last exon (Meyer et al., 2012), while circularization of the last exon of genes is uncommon (Zhang et al., 2014). These evidence suggested that circRNAs exhibit the patterns of m6A modification that are distinct from those of mRNAs.
m6A Regulators of circRNA
As the most abundant mRNA internal modification, the functions of m6A on mRNA are mediated through three homologous factors (Zaccara et al., 2019), namely so called “writers” (methyltransferase), “erasers” (demethylase), and “readers” (recognition). In circRNAs, m6A modification was found to be written, read and erased by the same regulators with mRNA (Zhou et al., 2017).
m6A Writers
The m6A is installed by a multi-subunit complex (Zaccara et al., 2019). This complex consists of two core components (methyltransferase-like 3 protein (METTL3) and METTL14) and other accessory regulatory subunits [Wilm’s tumor-1-associated protein (WTAP), VIRMA (Virilizer), E3 ubiquitin-protein ligase HAKAI, RNA recognition motif 15/15B (RBM15/15B), Zinc finger CCCH domain-containing protein 13 (ZC3H13)] (Zaccara et al., 2019; Gu et al., 2021; Jiang et al., 2021). METTL3 is the catalytic core of the multi-subunit complex (Bokar et al., 1997; Liu et al., 2014). METTL14 contains the RNA-binding site and is an allosteric activator of the enzymatic activity of METTL3 (Liu et al., 2014; Zaccara et al., 2019). WTAP is essential for m6A formation and it is responsible for localizing METTL3–METTL14 heterodimers to transcription sites (Ping et al., 2014). VIRMA may modulate region-selective methylation of sites in 3′ UTR and location near stop codons (Yue et al., 2018). RBM15/15B are the mediators of methylation specificity, which facilitate the recruitment of the writer complex to specific sites (Patil et al., 2016; Zaccara et al., 2019). ZC3H13 is a WTAP interactor, which acts as a bridge between RBM15 and WTAP (Knuckles et al., 2018). HAKAI functions as an E3 ubiquitin ligase, which was first identified in the WTAP interaction proteome (Horiuchi et al., 2013). Deletion of HAKAI results in a partial reduction in global m6A levels (Růžička et al., 2017). Notably, METTL16 is the new identified methyltransferase that can regulate S-adenosylmethionine homeostasis and can methylate long non-coding RNA and U6 small nuclear RNA (Brown et al., 2016; Pendleton et al., 2017; Warda et al., 2017). Currently, only few m6A residues in poly(A) RNA were found to be catalyzed by METTL16 (Brown et al., 2016; Pendleton et al., 2017; Warda et al., 2017).
m6A Erasers
m6A is a reversible process, which can be removed by demethylases. There are two demethylases, namely obesity-associated protein (FTO) and AlkB homolog 5 (ALKBH5). Both of them belong to the α-ketoglutarate-dependent dioxygenase AlkB-like superfamily and they remove m6A through Fe2+ and α-ketoglutarate-dependent manner (Jia et al., 2011; Zheng et al., 2013). Under physiological conditions, the demethylases appear to have a limited role, indicating that they function in specific tissues, or under specific stress and disease-relevant conditions (Zhao et al., 2017).
m6A Readers
m6A readers are a class of RNA binding proteins, which can recognize and bind to m6A sites, thereby resulting in different downstream effects, such as splicing, export and translation of m6A RNAs (Xiao et al., 2016; Roundtree et al., 2017b; Yang et al., 2017). In YT521-B homology (YTH) domain family, there are five m6A readers identified, namely YTHDC1, YTHDC2, YTHDF1, YTHDF2, and YTHDF3 (Liu et al., 2015). Among them, YTHDF1, YTHDF2, and YTHDF3 are mainly located in the cytoplasm. YTHDC1 is mainly found in the nucleus and YTHDC2 is located in both cytoplasm and nucleus. YTHDC1 affects their splicing (Kasowitz et al., 2018) and export (Roundtree et al., 2017b); YTHDC2 affects their degradation and translation (Hsu et al., 2017); YTHDF1 promotes translation (Wang et al., 2015); YTHDF2 promotes degradation (Du H. et al., 2016) and YTHDF3 promotes the translation and degradation (Shi et al., 2017). The other m6A readers include insulin-like growth factor 2 mRNA binding proteins 1-3 (IGF2BP1-3), heterogeneous nuclear ribonucleoprotein (HNRNP) family (HNRNPC, HNRNPA2B1, and HNRNPG), fragile X mental retardation protein (FMRP) and eukaryotic initiation factor 3 (eIF3). IGF2BPs can regulate gene expression by enhancing the stability of target RNAs (Huang H. et al., 2018). FMRP can promote nuclear export (Edens et al., 2019). eIF3 can facilitate cap-independent translation (Meyer et al., 2015).
Roles of m6A in circRNA Metabolism
Currently, several studies have found that m6A modification can affect the “fate” of m6A-modified circRNAs (Figure 1), including circRNAs biogenesis, translation, degradation and cellular localization, indicating that m6A plays important roles in circRNA metabolism. Through participating in circRNA metabolism, m6A can regulate many physiological and pathological processes. For example, during the nuclear phase, m6A can be read by specific nuclear readers and promote the nuclear export of circRNAs. Upon exporting to the cytoplasm, m6A can be read by specific cytosolic reader proteins, thereby affecting the function of circRNAs. Therefore, this section aimed to discuss the roles of m6A in circRNA metabolism comprehensively.
Figure 1. Roles of m6A in circRNA metabolism. It’s has confirmed that m6A can regulate circRNA metabolism, including circRNAs biogenesis, translation, degradation and cellular localization. On the one hand, m6A sites located around the start and stop codons in linear mRNAs can recruit spliceosome, leading to back splicing and circRNA production. On the other hand, m6A deposition and YTHDC1 binding to exons can regulate circularization. In the nucleus, the m6A can bind specific nuclear reader proteins, mainly YTHDC1, which can promote the export of circRNAs. Upon circRNAs export to the cytoplasm, m6A binds to specific reader proteins and other proteins to stabilize some mRNAs. The nuclear export of circRNAs also affects its miRNA sponges. The translation of circRNAs is only in cap-independent translation initiation mechanisms: IRES-dependent initiation of translation and m6A-dependent initiation of translation. m6A-driven translation of circRNA requires eIF4G2 and YTHDF3 and is enhanced by METTL3/14, inhibited by FTO. IRES-driven translation and m6A-driven translation may have interplays. Finally, m6A-modified circRNAs are endoribonuclease-cleaved via the YTHDF2-HRSP12-RNase P/MRP axis.
m6A and Biogenesis of circRNA
circRNA is conserved evolutionarily and is expressed in cell, tissue and developmental stage-specific expression patterns (Memczak et al., 2013; Hsiao et al., 2017; Nicolet et al., 2018), suggesting that circRNAs play an important regulatory role in different physiological and pathological processes. Therefore, it is essential to understand the regulation regarding the biogenesis of circRNAs. The biogenesis of circRNA is distinct from canonical splicing, and it is generated through a back-splicing orchestrated by the spliceosome machinery (Aufiero et al., 2019). Recently, Tang et al. (2020) reported that m6A could promote the biogenesis of circRNA in male germ cells. They found that, for open reading frames (ORFs)-containing circRNAs during murine spermatogenesis, the back splicing occurred mostly at m6A enriched sites, while these m6A sites were usually located around the start and stop codons in linear mRNAs. To further establish the cause–effect relationship between m6A and circRNAs, Tang et al. (2020) knocked out ALKBH5 and METTL3, respectively. After knocking out ALKBH5 in spermatogenic cells, the m6A level was significantly increased compared with wild-type controls. Consistently, the circRNAs abundance was increased in Alkbh5-null spermatogenic cells. After knocking out METTL3, they identified much fewer circRNAs. Similarly, Di Timoteo et al. (2020) also reported that the level of circRNA was decreased in HeLa cells after knocking out METTL3. Interestingly, they found the consistent decrease of circ-ZNF609 after knocking out YTHDC1 in HeLa, RD, RH4, and HEK293T cells. Further analysis found that m6A deposition and YTHDC1 binding to exons that undergo circularization were significantly correlative, and circRNAs containing total-only m6As were more represented among those decreasing upon depletion of both YTHDC1 and METTL3 compared with those affected by either METTL3 or YTHDC1 knockdown alone. These evidence indicated that m6A regulated the biogenesis of circRNAs in an METTL3/YTHDC1-dependent manner. In colorectal cancer, METTL3 could increase the expression of circ1662 (Chen et al., 2021). Therefore, based on the above evidence, we can confirm that m6A modification can increase the biogenesis of circRNAs. However, more detailed information is needed to understand how m6A regulates the biogenesis of circRNAs.
m6A and Translation of circRNA
Although circRNA is considered as non-coding RNA for a long time, accumulating evidence has shown that some circRNAs haveprotein-coding potential and can be translated (Legnini et al., 2017; Pamudurti et al., 2017; Yang et al., 2017, 2018; Zhang et al., 2018a, b; Liang et al., 2019). More recently, Huang et al. (2021) reported a database, TransCirc,1 to help investigate the circRNAs that have translation capacity. However, the mechanisms about the translation of circRNAs remain largely unknown. As we all know, the translation initiation of mRNA is dependent on the cap structure of 5′-end, which contains a 7-methyl guanosine (m7G) that can be recognized by the eukaryotic translation initiation factor 4E (eIF4E) (Gross et al., 2003). However, circRNA is a covalently closed RNA molecule without 5′ caps and 3′ tails (Meng X. et al., 2017), so that circRNA is translated in cap-independent translation initiation mechanisms. Currently, two mechanisms are reported in translation initiation of circRNA: internal ribosome entry site (IRES)-dependent initiation of translation and m6A-dependent initiation of translation. IRESs are sequences that form secondary structures on RNA and can initiate translation through recruiting ribosomes by IRES-transacting factors (ITAFs) in the absence of canonical translation initiation factors (Godet et al., 2019). At first, the fact that circRNAs can be translated is proved by IRES-driven pathway (Chen and Sarnow, 1995). Afterward, more and more circRNAs that can be translated by IRES-driven pathway were reported, such as circZNF609 (Legnini et al., 2017), circMbl (Pamudurti et al., 2017), circSHPRH (Zhang et al., 2018a), circFBXW7 (Yang et al., 2018), circLINC-PINT (Zhang et al., 2018b), circb-catenin (Liang et al., 2019).
Interestingly, the discovery of m6A that can initiate the translation of circRNA expands the coding landscape of human transcriptome (Yang et al., 2017). Yang et al. (2017) found that the m6A-driven translation of circRNAs was widespread, with hundreds of endogenous circRNAs having translation potential, by the analysis of polysome profiling, computational prediction and mass spectrometry. Mechanistically, they found that the m6A-driven translation of circRNA required initiation factor eIF4G2 and m6A reader YTHDF3, and was enhanced by methyltransferase METTL3/14, inhibited by demethylase FTO (Yang et al., 2017). When inserting a short fragment (19nt) containing different copies of consensus m6A motifs (RRACH) before the start codon of circRNA reporter in 293 cells, the protein could be detected. Moreover, single m6A site has similar translation efficiency compared to circRNA with two m6A sites (Yang et al., 2017). These evidence suggested that m6A motif was essential for initiating translation. In human papillomavirus (HPV), m6A-modified circE7 could be translated to produce E7 oncoprotein (Zhao et al., 2019).
In addition to initiating the translation of circRNAs, m6A may also regulate the translation of circRNAs (Yang et al., 2017). For example, circZNF609 could be translated by IRES-driven pathway (Legnini et al., 2017). A recent study reported that two mutant at m6A sites displayed approximately 50% reduction of circ-ZNF609 translation (Di Timoteo et al., 2020). Further analysis found that m6A regulated its translation through recognition by YTHDF3 and eIF4G2 (Di Timoteo et al., 2020). This study suggested that the two cap-independent translation of circRNA might have interplays. However, the relationships between the two cap-independent translations need further investigations.
In mRNAs, studies have shown that m6A regulates their translation under stress response (Zhou et al., 2015, 2018). Yang et al. (2017) found that m6A-mediated circRNA translation was increased under the heat shock condition and its mechanism might be the translocation of YTHDF2 from cytosol into nucleus upon heat shock to block the m6A “eraser” FTO (Zhou et al., 2015; Yang et al., 2017). These evidence suggested that m6A-mediated circRNA translation played important roles in cellular stress response. In the future, to broaden our understanding regarding the significance of circRNA translation in the specific conditions, it is essential to explore the conditions that affect the m6A-mediated circRNA translation.
m6A and Degradation of circRNA
Due to the covalently closed circular structure without 5′ caps and3′ tails, circRNAs can resist the foreign chemicals or exonucleases (Jeck et al., 2013; Aufiero et al., 2019). Therefore, circRNAs are more stable than the linear RNAs and are not degraded readily by RNase R, resulting in their long half-lives exceeding 48 h (Jeck et al., 2013; Suzuki and Tsukahara, 2014). Thus, there are few studies that reported the degradation of circRNA.
Hansen et al. (2011) reported that circRNAs with near perfect complementary miRNA target sites could be degraded in an Ago2-slicer-dependent manner. Ago2 is a member of the Argonaute protein family, which is important for RNA interference (Suzuki and Tsukahara, 2014) and can mediate the degradation of circRS-7 (Hansen et al., 2011). However, only a few circRNAs have miRNA sponge function or specific miRNA target sites. Jia et al. (2019) found that the depletion of GW182 (a key component of the P-body and RNA interference machine) could lead to the accumulation of endogenous circular transcripts. However, the depletion of other P-body components or RNAi complex factors did not have similar effects. So, there are more studies to explain this phenomenon.
A recent study reported that m6A-containing circRNAs could be degraded through YTHDF2–HRSP12–RNase P/MRP axis (Park et al., 2019). Mechanistically, human heat-responsive protein 12 (HRSP12) acted as an adaptor protein that linked YTHDC2 and RNase P/MRP to form YTHDF2-HRSP12-RNase P/MRP complex. Notably, after downregulating the expression YTHDF2, HRSP12, or POP1, Park et al. (2019) observed only three out of all tested 11 circRNAs commonly increased, indicating that this mechanism could only cleave a subset of m6A-containing circRNAs.
m6A and Cellular Localization of circRNA
As we all know, a large number of circRNAs function in cytoplasm, such as miRNA sponges (Thomson and Dinger, 2016). Thus, it is essential to understand how the circRNAs export from the nucleus to the cytoplasm. Huang C. et al. (2018) found that Drosophila Hel25E and its human homologs, UAP56/URH49, could regulate circRNAs localization and control the efficiency of nuclear export by measuring the lengths of mature circRNAs. Chen et al. (2021) reported that YTHDC1 could bind to the m6A motif GAACU of circNSUN2 and facilitate circNSUN2 export from the nucleus to the cytoplasm, further promote colorectal liver metastasis by stabilizing HMGA2 (Chen R. X. et al., 2019). In addition, the m6A-modified circ1662 could bind to YAPI and accelerate its nuclear transport (Chen et al., 2021). However, more studies are required to investigate the details that m6A regulates the cellular localization of circRNAs.
m6A circRNAS and Immunity
m6A in mRNA can control various aspects of immunity, including immune recognition, activation of innate and adaptive immune responses (Shulman and Stern-Ginossar, 2020). However, its roles in circRNA immunity are unknown. Several recent studies have discovered that m6A can also control circRNA immunity (Figure 2), including the roles in innate immune and anti-tumor immunity.
Figure 2. Roles of m6A-modified circRNAs in immunity. m6A can mark “self” circRNAs and these endogenous circRNAs inhibite innate immune responses, whereas exogenous circRNAs can activate innate immunity by activating RIG-I in the presence of K63-polyubiquitin. In tumor immunity, circRNAs may inhibit tumor progression by inducing anti-tumor immunity.
The innate immune system relies on pattern recognition receptors(PRRs) recognizing pathogen associated molecular patterns (PAMPs) to recognize “non-self” (Reikine et al., 2014; Wu and Chen, 2014). The PRRs, such as RIG-I and MDA5, can sense foreign nucleic acids. MDA5 can detect long double-stranded RNA (dsRNA) and RIG-I can recognize 5′ triphosphate on short dsRNAs (Schlee, 2013; Reikine et al., 2014; Wu and Chen, 2014). Reportedly, some exogenous circRNAs can activate antiviral and immune gene expression (Chen et al., 2017), while some endogenous circRNAs can inhibit protein kinase R and set the threshold for innate immunity upon virus infection (Liu et al., 2019). However, the mechanisms that the innate immune system recognizes self and non-self circRNAs remain unclear. Chen et al. (2021) found that m6A can mark circRNA as “self” and these endogenous circRNAs inhibited innate immune responses (Chen Y. G. et al., 2019). This immunosuppression was mediated by binding to YTHDF2 directly and inhibiting RIG-I activation. In addition, exogenous circRNAs activated innate immunity by directly binding RIG-I and K63-polyubiquitin chain and they together formed a three-component signaling-competent complex for immune signaling (Chen Y. G. et al., 2019). These evidence suggested that, for human circRNAs, m6A modification was essential for the recognition function of innate immunity.
Besides participating in the recognition of innate immunity, m6A-modified circRNAs can inhibit the tumor growth by inducing anti-tumor immunity. Chen et al. (2021) reported that circFOREIGN could induce CD8+ T cell responses, and after receiving circFOREIGN, the growth of tumor in mice model was lower than the negative control, indicating that circFOREIGN would can induce adaptive immunity against OVA-expressing tumors. Notably, Li et al. (2021) reported that the immune responses induced by m6A-modified circNDUFB2 could inhibit NSCLC progression. The overexpression of circNDUFB2 markedly inhibited the tumorigenicity of LLC1 cells (a murine lung carcinoma cell line) in vivo, and the level of IFN-β in serum was significantly increased in the mice with circNDUFB2-overexpressed LLC1 cells. Moreover, the infiltration of CD8+ T cells and DCs in the tumor tissues were significantly increased. These results suggested that m6A-modified circNDUFB2 inhibited tumor progression by activating anti-tumor immunity. However, more studies are needed to explore the roles of m6A modification in circRNAs in the field of tumor immunity.
These evidence demonstrated the roles of m6A modification in circRNAs during recognizing innate immunity and inducing anti-tumor immunity, but major knowledge gaps in field of immunity remain to be filled, such as its specific roles in different immune diseases and immune cells.
m6A circRNA and Diseases
Although it was not until 2017 that circRNAs were found to be modified by m6A, many studies have shown that m6A modification of circRNAs plays important roles in physiological processes and disease progressions (Figure 3), such as tumors (Chen R. X. et al., 2019; Chen et al., 2021; Li et al., 2021), acute coronary syndrome (Guo et al., 2020), pulmonary hypertension (Su et al., 2020), and age-related cataract (Li et al., 2020). In this section, we focused on the roles of m6A-modified circRNAs in physiological processes and disease progressions.
Figure 3. Roles of m6A-modified circRNAs in diseases. It has confirmed that m6A modification of circRNAs plays important roles in the occurrence and development of human diseases. In tumors, m6A-modified circRNAs are not only involved in their occurrence, development and metastasis, but also in the resistance of chemotherapy and radiotherapy by regulating the metabolism and functions of circRNAs. In acute coronary syndrome, IRF-1 facilitates macrophage pyroptosis and inflammation by inhibiting m6A-modified circ_0029589, thereby promoting the formation of necrotic core. In hypoxia mediated pulmonary hypertension, m6A may affect the interaction between circRNAs and miRNAs to participate in its development. In age-related cataract, many differentially expressed m6A-modified circRNAs are identified, but the specific roles remain unclear.
m6A circRNA and Tumor
m6A modification of circRNAs has been shown to be involved in theprogression of many tumors. In patients with colorectal carcinoma (CRC), metastasis is the overwhelming cause of death, and the liver is the most frequent distant metastatic site (Hofheinz and Stintzing, 2019). Recently, Chen R. X. et al. (2019) reported a novel metastatic mechanism to liver of CRC. First, YTHDC1 binds to the m6A motif of circNSUN2 and facilitates circNSUN2 export from the nucleus to the cytoplasm. Next, IGF2BP2 binds to the CAUCAU motif of circNSUN2 through the KH3-4 di-domain and the AAACA site inside of circNSUN2 can directly bind to the 3′UTR of HMGA2 with AU-Rich Elements, ultimately forming circNSUN2/IGF2BP2/HMGA2 RNA-protein ternary complex to stabilize HMGA2 mRNA, thereby promoting colorectal liver metastasis. This study provided the comprehensive evidence that m6A modification of circRNAs may contribute to cancer therapy (He et al., 2021). Another study found that m6A could increase the expression of circ1662 in CRC and circ1662 directly bound to yes-associated protein 1 (YAP1) and accelerated its nuclear accumulation to regulate the SMAD3 pathway, thereby promoting the invasion and migration of CRC cells (Chen et al., 2021).
In Hepatitis B virus (HBV)-associated hepatocellular carcinoma (HCC), circARL3 was upregulated in HBV positive HCC tissues and knockdown of circARL3 inhibited the proliferation and invasion of HBV + HCC cells (Rao et al., 2021). Rao et al. (2021) found that HBx protein could upregulate the expression of METTL3, increasing the m6A modification of circARL3. m6A-modified circARL3 can be read by YTHDC1 and regulate its biogenesis. Finally, circARL3 promoted HBV + HCC progression by sponging miR-1305.
In non-small cell lung cancer (NSCLC), circNDUFB2 was downregulated and elevated circNDUFB2 could inhibit growth and metastasis of NSCLC cells (Li et al., 2021). circNDUFB2 functions as a scaffold to enhance the interaction between TRIM25 and IGF2BPs and can from a TRIM25/circNDUFB2/IGF2BPs ternary complex. By the way, TRIM25 is an RNA-binding protein and belongs to the Tripartite Motif (TRIM) family of E3 ubiquitin ligases, which catalyzes the addition of polyubiquitin chains to its substrates for degradation (Lee et al., 2018). This ternary complex facilitates ubiquitination and degradation of IGF2BPs during NSCLC progression. m6A modification of circNDUFB2 can enhance this effect by regulating the interactions between circNDUFB2 and IGF2BPs. In addition, circNDUFB2 participates in the activation of anti-tumor immunity during NSCLC progression by activating the RIG-I–MAVS pathway (Li et al., 2021), indicating that circNDUFB2 inhibits NSCLC progression in a dual-role pattern.
In cervical cancer, Zhou et al. (2017) identified a large number of m6A-modified circRNAs in HeLa cells, indicating that m6A-modified circRNAs may play an important regulatory role in cervical cancer. As we all know, human papillomavirus (HPV) infection is the main cause that leads to the occurrence of cervical cancer (Chibwesha and Stringer, 2019). Zhao et al. (2019) reported that HPV could generate the circE7, which was an oncoprotein-encoding circRNA. circE7 is m6A-modified and can be translated to produce E7 oncoprotein in an m6A-independent manner. Specific disruption of circE7 can reduce E7 protein levels in CaSki cervical carcinoma cells and can inhibit cancer cell growth both in vitro and in tumor xenografts.
In poorly differentiated gastric adenocarcinoma (PDGA), Zhang et al. (2020) identified 20 differentiated expressed m6A-modified circRNAs by SRAMP analysis, m6A-immunoprecipitation and real-time quantitative PCR. They found that the trend of m6A modification alteration was mainly consistent with the expression level of circRNAs. However, the roles of these m6A-modified circRNAs in PDGA remain unclear and need to be explore extensively.
In human ameloblastoma, Niu et al. (2020) investigated the global m6A modification, including in mRNAs, lncRNAs and circRNAs. They found 3,673 differentially methylated m6A sites in mRNAs, 4,975 differentially methylated m6A sites in lncRNAs and 364 differentially methylated m6A sites within circRNAs. Although a global data of m6A in ameloblastoma was provided, their roles and the relationships among these m6A-modified mRNAs, lncRNAs and circRNAs need to be further explored.
Based on the above analysis, we can confirm that m6A modification of circRNAs plays critical roles in occurrence, development and metastasis of tumors. Of note, besides the above roles, m6A modification of circRNAs is also related to resistance of chemotherapy and radiotherapy. In sorafenib-resistant HCC, m6A could increase the expression of circRNA-SORE and increased circRNA-SORE sequestered miR-103a-2-5p and miR-660-3p by acting as a miRNA sponge, thereby competitively activating the Wnt/β-catenin pathway and inducing sorafenib resistance (Xu et al., 2020). In hypopharyngeal squamous cell carcinoma (HPSCC), m6A could stabilize the expression of circCUX1 and knocking down circCUX1 promoted the sensitivity of hypopharyngeal cancer cells to radiotherapy. Mechanistically, circCUX1 bound to caspase1 and inhibited its expression, leading to a decrease in the release of inflammatory factors, thereby developing tolerance to radiotherapy (Wu et al., 2021). As a key inflammation-related molecule, caspase1 can affect the occurrence, development, invasion, and metastasis of tumors by regulating the tumor inflammatory microenvironment (Hu et al., 2010), and caspase-1 actively induced tumor cell programmed death and anti-tumor immune surveillance (Zitvogel et al., 2012). Therefore, these evidence provided the new insights into understanding the resistant to radiotherapy.
m6A circRNA and Acute Coronary Syndrome
Reportedly, macrophage death plays a key role in promoting the formation of necrotic core and plaque disruption in late atherosclerotic lesions (Hizir et al., 2017; Xu et al., 2018; Guo et al., 2019). Pyroptosis is a newly identified inflammatory programmed cell death with the activation of caspase-1 and the consequent release of inflammatory cytokines (Liu et al., 2016), and has been found to have a strong correlation with ox-LDL-induced human macrophage death and aggravates the instability of atherosclerotic lesions (Yan et al., 2016; Martinet et al., 2019). Guo et al. (2019) found that IFN regulatory factor (IRF)-1 could activate ox-LDL-induced macrophage pyroptosis in atherosclerosis (AS) and acute coronary syndrome (ACS). However, the underlying mechanism is unknown. A recent study found that the expression of hsa_circ_0029589 was significantly decreased in ACS and the relative m6A level of circ_0029589 was obviously elevated (Guo et al., 2020). Inhibition of IRF-1 obviously induced the expression of circ_0029589 and reduced the m6A level of circ_0029589, and knocking down the expression of METTL3 markedly increased circ_0029589 expression and decreased the m6A level of circ_0029589 and caspase-1 activity in AS and ACS (Guo et al., 2020). These results suggested that IRF-1 facilitated macrophage pyroptosis and inflammation in ACS by inhibiting circ_0029589 through promoting its m6A modification.
m6A circRNA and Pulmonary Hypertension
Pulmonary hypertension is a lethal disease and hypoxia mediated pulmonary hypertension (HPH) belongs to group III pulmonary hypertension, which is induced by chronic hypoxia (Pugliese et al., 2015; Galiè et al., 2016). To investigate the significance of m6A circRNAs in HPH, Su et al. (2020) analyzed transcriptome-wide map of m6A circRNAs in a rat model by high-throughput m6A and circRNAs sequencing. They found that m6A abundance in circRNAs was significantly reduced under hypoxia in vitro and m6A-modified circXpo6 and circTmtc3 might be involved in HPH through regulating Wnt and FoxO signaling pathways by influencing the interactions between circRNAs and miRNAs. However, there are more evidence to explore the roles of m6A-modified circXpo6 and circTmtc3 in HPH and their bindings to miRNAs need to be validated by dual-luciferase reporter assay.
m6A circRNA and Age-Related Cataract
Age-related cataract (ARC) is one of the leading causes of vision impairment and accounts for the majority of senile blindness worldwide (Khairallah et al., 2015). However, the pathogenesis of ARC remains unclear. A recent study found that there are a lot of m6A circRNAs in the lens epithelium cells from the patients with ARC (Li et al., 2020). The level of m6A abundance in total circRNAs was decreased compared with the controls and ALKBH5 was significantly upregulated in ARCs. The authors predicted the potential functions of m6A modified circRNAs and found the relevant pathways that may be associated with m6A modified circRNAs. These results indicated that m6A modification of circRNAs may play an important role in ARC. However, a large number of studies are needed to investigate their roles and validate the predicted the potential functions of m6A modified circRNAs.
Based on the above analysis regarding the roles of m6A modified circRNAs in different diseases, we can deduce that m6A modification in circRNAs has following common mechanisms to participate in disease progression. (1) Under pathological conditions, the expression of methyltransferases and/or demethylases is upregulated or downregulated, leading to the changes in m6A levels of circRNAs. (2) Through recognizing m6A motif by different m6A readers, m6A modification can regulate the metabolism and functions of circRNAs, thereby participating in disease progression (Figures 2, 3). However, the relationships among m6A modification, circRNA functions and diseases remain largely unclear. It is confirmed that m6A modification is as the upstream regulatory mechanism to affect downstream molecules by regulating circRNAs directly. In the future, the conjoint analysis of MeRIP-seq and RNA-seq may be helpful for seeking disease-specific m6A modified circRNAs. After that, qRT-PCR, RIP, RNA pull down, and dual-luciferase reporter assay may be used for exploring the mechanism of selected m6A modified circRNAs.
Conclusion and Future Directions
No doubt, increasing number of studies have validated that circRNAs play key roles in human development and disease progressions; and epigenetic modification, especially m6A in mRNA, has emerged as widespread regulatory mechanisms that control gene expression in diverse physiological and pathological processes. However, the roles of epigenetic modifications in circRNAs remain largely unknown. Despite there are only 4 years regarding the researches of m6A in circRNAs, accumulating studies have identified numerous m6A-modified circRNAs and have found their important regulatory roles in development and diseases. In conclusion, current studies have found that these m6A-modified circRNAs are expressed in cell-type and disease-specific methylation patterns. In functions, m6A modification can regulate circRNAs metabolism, including circRNAs biogenesis, translation, degradation and cellular localization. Importantly, these m6A-modified circRNAs participate in diverse physiological and pathological processes, such as immunity, tumors, acute coronary syndrome, pulmonary hypertension and age-related cataract, by regulating circRNAs metabolism and functions. In general, the discovery of m6A modification in circRNAs broadened our horizons in RNA epigenetics and we provided a comprehensive understanding about the m6A modification of circRNAs and their roles based on current reports.
In addition, m6A-modified circRNAs may become the potential therapeutic target of diseases. For example, m6A-modified circNSUN2 has been found to play a key role in promoting colorectal liver metastasis (Chen R. X. et al., 2019; He et al., 2021). Targeted inhibition of circNSUN2 may block colorectal liver metastasis. Here, we listed the prospect targets that may be used for clinical treatment in the future in Table 1 based on current reports. However, the research regarding targeting m6A-modified circRNAs in diseases is still in infancy and many studies are needed to seek disease-specific targets and to develop effective methods for detection of m6A-modified circRNAs.
Table 1. m6A-modified circRNAs as the potential therapeutic target of different diseases based on current reports.
Currently, many questions are about when and how m6A is added on or removed from circRNAs, and about how m6A regulates circRNAs metabolism and disease progressions. For example, Tang et al. (2020) found that, during murine spermatogenesis, for a subset of circRNAs, the back splicing occurs mostly at m6A enriched sites, which are usually located around the start and stop codons in linear mRNAs. Zhou et al. (2017) found that numerous m6A circRNAs were generated from exons that didn’t contain m6A peaks in corresponding mRNAs. The question is whether m6A was added before or after the circRNAs were formed. In addition, the causes that lead to the change of m6A modification in circRNAs are unclear.
Another challenge is the technologies to study m6A modification. For example, MeRIP-Seq is the main method to identify m6A modifications using m6A-specific antibodies. A limitation of this method is that it cannot provide the precise location of m6A at single-nucleotide resolution. Ho-Xuan et al. (2020) found that endogenous proteins could not be detected when they constructed the overexpression of circZNF609. Comprehensive mutational analysis found that deletion constructs, which were deficient in producing circZNF609, still generated the observed protein products, suggesting that the apparent circZNF609 translation originated from transsplicing byproducts of the overexpression plasmids. Therefore, the overexpression construction of circRNAs need to be evaluated carefully. In the future, studies should overcome and solve these challenges and questions in order to better understand the roles of m6A in circRNAs biology further.
Author Contributions
HS and LT designed and executed the study. JW and XG collected and analyzed the data, and wrote the manuscript. YW and SH made the figures. XY was responsible for language quality. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by the National Natural Science Foundation of China (Grant No. 81772001) and the National Clinical Key Subject of China (Grant No. 41732113).
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.
Footnotes
References
Aufiero, S., Reckman, Y. J., Pinto, Y. M., and Creemers, E. E. (2019). Circular RNAs open a new chapter in cardiovascular biology. Nat. Rev. Cardiol. 16, 503–514. doi: 10.1038/s41569-019-0185-2
Aufiero, S., Van Den Hoogenhof, M. M. G., Reckman, Y. J., Beqqali, A., Van Der Made, I., Kluin, J., et al. (2018). Cardiac circRNAs arise mainly from constitutive exons rather than alternatively spliced exons. RNA 24, 815–827. doi: 10.1261/rna.064394.117
Barbagallo, D., Caponnetto, A., Brex, D., Mirabella, F., Barbagallo, C., Lauretta, G., et al. (2019). CircSMARCA5 regulates VEGFA mRNA splicing and angiogenesis in glioblastoma multiforme through the binding of SRSF1. Cancers (Basel) 11:194. doi: 10.3390/cancers11020194
Bokar, J. A., Shambaugh, M. E., Polayes, D., Matera, A. G., and Rottman, F. M. (1997). Purification and cDNA cloning of the AdoMet-binding subunit of the human mRNA (N6-adenosine)-methyltransferase. RNA 3, 1233–1247.
Brown, J. A., Kinzig, C. G., Degregorio, S. J., and Steitz, J. A. (2016). Methyltransferase-like protein 16 binds the 3′-terminal triple helix of MALAT1 long noncoding RNA. Proc. Natl. Acad. Sci. U.S.A. 113, 14013–14018. doi: 10.1073/pnas.1614759113
Chen, C., Yuan, W., Zhou, Q., Shao, B., Guo, Y., Wang, W., et al. (2021). N6-methyladenosine-induced circ1662 promotes metastasis of colorectal cancer by accelerating YAP1 nuclear localization. Theranostics 11, 4298–4315. doi: 10.7150/thno.51342
Chen, C. Y., and Sarnow, P. (1995). Initiation of protein synthesis by the eukaryotic translational apparatus on circular RNAs. Science 268, 415–417. doi: 10.1126/science.7536344
Chen, L. L., and Yang, L. (2015). Regulation of circRNA biogenesis. RNA Biol 12, 381–388. doi: 10.1080/15476286.2015.1020271
Chen, R. X., Chen, X., Xia, L. P., Zhang, J. X., Pan, Z. Z., Ma, X. D., et al. (2019). N(6)-methyladenosine modification of circNSUN2 facilitates cytoplasmic export and stabilizes HMGA2 to promote colorectal liver metastasis. Nat. Commun. 10:4695. doi: 10.1038/s41467-019-12651-2
Chen, X., Yang, T., Wang, W., Xi, W., Zhang, T., Li, Q., et al. (2019). Circular RNAs in immune responses and immune diseases. Theranostics 9, 588–607. doi: 10.7150/thno.29678
Chen, Y. G., Chen, R., Ahmad, S., Verma, R., Kasturi, S. P., Amaya, L., et al. (2019). N6-Methyladenosine modification controls circular RNA immunity. Mol. Cell 76, 96.e9–109.e9.
Chen, Y. G., Kim, M. V., Chen, X., Batista, P. J., Aoyama, S., Wilusz, J. E., et al. (2017). Sensing self and foreign circular RNAs by intron identity. Mol. Cell 67, 228.e5–238.e5.
Cheng, Z., Yu, C., Cui, S., Wang, H., Jin, H., Wang, C., et al. (2019). circTP63 functions as a ceRNA to promote lung squamous cell carcinoma progression by upregulating FOXM1. Nat. Commun. 10:3200.
Chibwesha, C. J., and Stringer, J. S. A. (2019). Cervical cancer as a global concern: contributions of the dual epidemics of HPV and HIV. JAMA 322, 1558–1560. doi: 10.1001/jama.2019.16176
Conn, S. J., Pillman, K. A., Toubia, J., Conn, V. M., Salmanidis, M., Phillips, C. A., et al. (2015). The RNA binding protein quaking regulates formation of circRNAs. Cell 160, 1125–1134. doi: 10.1016/j.cell.2015.02.014
Desrosiers, R., Friderici, K., and Rottman, F. (1974). Identification of methylated nucleosides in messenger RNA from Novikoff hepatoma cells. Proc. Natl. Acad. Sci. U.S.A. 71, 3971–3975. doi: 10.1073/pnas.71.10.3971
Di Timoteo, G., Dattilo, D., Centrón-Broco, A., Colantoni, A., Guarnacci, M., Rossi, F., et al. (2020). Modulation of circRNA metabolism by m(6)A modification. Cell Rep. 31:107641. doi: 10.1016/j.celrep.2020.107641
Du, H., Zhao, Y., He, J., Zhang, Y., Xi, H., Liu, M., et al. (2016). YTHDF2 destabilizes m(6)A-containing RNA through direct recruitment of the CCR4-NOT deadenylase complex. Nat. Commun. 7:12626.
Du, W. W., Yang, W., Liu, E., Yang, Z., Dhaliwal, P., and Yang, B. B. (2016). Foxo3 circular RNA retards cell cycle progression via forming ternary complexes with p21 and CDK2. Nucleic Acids Res. 44, 2846–2858. doi: 10.1093/nar/gkw027
Edens, B. M., Vissers, C., Su, J., Arumugam, S., Xu, Z., Shi, H., et al. (2019). FMRP modulates neural differentiation through m(6)A-Dependent mRNA nuclear export. Cell Rep. 28, 845.e5–854.e5.
Frye, M., Harada, B. T., Behm, M., and He, C. (2018). RNA modifications modulate gene expression during development. Science 361, 1346–1349. doi: 10.1126/science.aau1646
Galiè, N., Humbert, M., Vachiery, J. L., Gibbs, S., Lang, I., Torbicki, A., et al. (2016). 2015 ESC/ERS guidelines for the diagnosis and treatment of pulmonary hypertension. Rev. Esp. Cardiol. (Engl Ed) 69:177. doi: 10.1016/j.rec.2016.01.002
Godet, A. C., David, F., Hantelys, F., Tatin, F., Lacazette, E., Garmy-Susini, B., et al. (2019). IRES trans-acting factors, key actors of the stress response. Int. J. Mol. Sci. 20:924. doi: 10.3390/ijms20040924
Gross, J. D., Moerke, N. J., Von Der Haar, T., Lugovskoy, A. A., Sachs, A. B., Mccarthy, J. E., et al. (2003). Ribosome loading onto the mRNA cap is driven by conformational coupling between eIF4G and eIF4E. Cell 115, 739–750. doi: 10.1016/s0092-8674(03)00975-9
Gu, J., Zhan, Y., Zhuo, L., Zhang, Q., Li, G., Li, Q., et al. (2021). Biological functions of m(6)A methyltransferases. Cell Biosci. 11:15.
Guo, M., Yan, R., Ji, Q., Yao, H., Sun, M., Duan, L., et al. (2020). IFN regulatory Factor-1 induced macrophage pyroptosis by modulating m6A modification of circ_0029589 in patients with acute coronary syndrome. Int. Immunopharmacol. 86:106800. doi: 10.1016/j.intimp.2020.106800
Guo, M., Yan, R., Yao, H., Duan, L., Sun, M., Xue, Z., et al. (2019). IFN regulatory factor 1 mediates macrophage pyroptosis induced by oxidized low-density lipoprotein in patients with acute coronary syndrome. Mediators Inflamm. 2019:2917128.
Hansen, T. B., Wiklund, E. D., Bramsen, J. B., Villadsen, S. B., Statham, A. L., Clark, S. J., et al. (2011). miRNA-dependent gene silencing involving Ago2-mediated cleavage of a circular antisense RNA. EMBO J. 30, 4414–4422. doi: 10.1038/emboj.2011.359
He, R. Z., Jiang, J., and Luo, D. X. (2020). The functions of N6-methyladenosine modification in lncRNAs. Genes Dis. 7, 598–605. doi: 10.1016/j.gendis.2020.03.005
He, R. Z., Jiang, J., and Luo, D. X. (2021). M6A modification of circNSUN2 promotes colorectal liver metastasis. Genes Dis. 8, 6–7. doi: 10.1016/j.gendis.2019.12.002
Hizir, Z., Bottini, S., Grandjean, V., Trabucchi, M., and Repetto, E. (2017). RNY (YRNA)-derived small RNAs regulate cell death and inflammation in monocytes/macrophages. Cell Death Dis. 8:e2530. doi: 10.1038/cddis.2016.429
Hofheinz, R. D., and Stintzing, S. (2019). Study evidence confirms current clinical practice in refractory metastatic colorectal cancer: the ReDOS trial. Lancet Oncol. 20, 1036–1037. doi: 10.1016/s1470-2045(19)30390-0
Horiuchi, K., Kawamura, T., Iwanari, H., Ohashi, R., Naito, M., Kodama, T., et al. (2013). Identification of Wilms’ tumor 1-associating protein complex and its role in alternative splicing and the cell cycle. J. Biol. Chem. 288, 33292–33302. doi: 10.1074/jbc.m113.500397
Ho-Xuan, H., Glažar, P., Latini, C., Heizler, K., Haase, J., Hett, R., et al. (2020). Comprehensive analysis of translation from overexpressed circular RNAs reveals pervasive translation from linear transcripts. Nucleic Acids Res. 48, 10368–10382. doi: 10.1093/nar/gkaa704
Hsiao, K. Y., Sun, H. S., and Tsai, S. J. (2017). Circular RNA - New member of noncoding RNA with novel functions. Exp. Biol. Med. (Maywood) 242, 1136–1141. doi: 10.1177/1535370217708978
Hsu, M. T., and Coca-Prados, M. (1979). Electron microscopic evidence for the circular form of RNA in the cytoplasm of eukaryotic cells. Nature 280, 339–340. doi: 10.1038/280339a0
Hsu, P. J., Zhu, Y., Ma, H., Guo, Y., Shi, X., Liu, Y., et al. (2017). Ythdc2 is an N(6)-methyladenosine binding protein that regulates mammalian spermatogenesis. Cell Res. 27, 1115–1127. doi: 10.1038/cr.2017.99
Hu, B., Elinav, E., Huber, S., Booth, C. J., Strowig, T., Jin, C., et al. (2010). Inflammation-induced tumorigenesis in the colon is regulated by caspase-1 and NLRC4. Proc. Natl. Acad. Sci. U.S.A. 107, 21635–21640. doi: 10.1073/pnas.1016814108
Huang, C., Liang, D., Tatomer, D. C., and Wilusz, J. E. (2018). A length-dependent evolutionarily conserved pathway controls nuclear export of circular RNAs. Genes Dev. 32, 639–644. doi: 10.1101/gad.314856.118
Huang, H., Weng, H., Sun, W., Qin, X., Shi, H., Wu, H., et al. (2018). Recognition of RNA N(6)-methyladenosine by IGF2BP proteins enhances mRNA stability and translation. Nat. Cell. Biol. 20, 285–295. doi: 10.1038/s41556-018-0045-z
Huang, W., Ling, Y., Zhang, S., Xia, Q., Cao, R., Fan, X., et al. (2021). TransCirc: an interactive database for translatable circular RNAs based on multi-omics evidence. Nucleic Acids Res. 49, D236–D242.
Ivanov, A., Memczak, S., Wyler, E., Torti, F., Porath, H. T., Orejuela, M. R., et al. (2015). Analysis of intron sequences reveals hallmarks of circular RNA biogenesis in animals. Cell Rep. 10, 170–177. doi: 10.1016/j.celrep.2014.12.019
Jeck, W. R., Sorrentino, J. A., Wang, K., Slevin, M. K., Burd, C. E., Liu, J., et al. (2013). Circular RNAs are abundant, conserved, and associated with ALU repeats. RNA 19, 141–157. doi: 10.1261/rna.035667.112
Jia, G., Fu, Y., Zhao, X., Dai, Q., Zheng, G., Yang, Y., et al. (2011). N6-methyladenosine in nuclear RNA is a major substrate of the obesity-associated FTO. Nat. Chem. Biol. 7, 885–887. doi: 10.1038/nchembio.687
Jia, R., Xiao, M. S., Li, Z., Shan, G., and Huang, C. (2019). Defining an evolutionarily conserved role of GW182 in circular RNA degradation. Cell Discov. 5:45.
Jiang, X., Liu, B., Nie, Z., Duan, L., Xiong, Q., Jin, Z., et al. (2021). The role of m6A modification in the biological functions and diseases. Signal Transduct. Target Ther. 6:74.
Kasowitz, S. D., Ma, J., Anderson, S. J., Leu, N. A., Xu, Y., Gregory, B. D., et al. (2018). Nuclear m6A reader YTHDC1 regulates alternative polyadenylation and splicing during mouse oocyte development. PLoS Genet. 14:e1007412. doi: 10.1371/journal.pgen.1007412
Khairallah, M., Kahloun, R., Bourne, R., Limburg, H., Flaxman, S. R., Jonas, J. B., et al. (2015). Number of people blind or visually impaired by cataract worldwide and in world regions, 1990 to 2010. Invest Ophthalmol. Vis. Sci. 56, 6762–6769. doi: 10.1167/iovs.15-17201
Knuckles, P., Lence, T., Haussmann, I. U., Jacob, D., Kreim, N., Carl, S. H., et al. (2018). Zc3h13/Flacc is required for adenosine methylation by bridging the mRNA-binding factor Rbm15/Spenito to the m(6)A machinery component Wtap/Fl(2)d. Genes Dev. 32, 415–429. doi: 10.1101/gad.309146.117
Kramer, M. C., Liang, D., Tatomer, D. C., Gold, B., March, Z. M., Cherry, S., et al. (2015). Combinatorial control of Drosophila circular RNA expression by intronic repeats, hnRNPs, and SR proteins. Genes Dev. 29, 2168–2182. doi: 10.1101/gad.270421.115
Lee, J. M., Choi, S. S., Lee, Y. H., Khim, K. W., Yoon, S., Kim, B. G., et al. (2018). The E3 ubiquitin ligase TRIM25 regulates adipocyte differentiation via proteasome-mediated degradation of PPARγ. Exp. Mol. Med. 50, 1–11. doi: 10.1038/s12276-018-0162-6
Legnini, I., Di Timoteo, G., Rossi, F., Morlando, M., Briganti, F., Sthandier, O., et al. (2017). Circ-ZNF609 is a circular RNA that can be translated and functions in myogenesis. Mol. Cell 66, 22.e9–37.e9.
Li, B., Zhu, L., Lu, C., Wang, C., Wang, H., Jin, H., et al. (2021). circNDUFB2 inhibits non-small cell lung cancer progression via destabilizing IGF2BPs and activating anti-tumor immunity. Nat. Commun. 12:295.
Li, P., Yu, H., Zhang, G., Kang, L., Qin, B., Cao, Y., et al. (2020). Identification and characterization of N6-Methyladenosine CircRNAs and methyltransferases in the lens epithelium cells from age-related cataract. Invest. Ophthalmol. Vis. Sci. 61:13. doi: 10.1167/iovs.61.10.13
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
Liang, W. C., Wong, C. W., Liang, P. P., Shi, M., Cao, Y., Rao, S. T., et al. (2019). Translation of the circular RNA circβ-catenin promotes liver cancer cell growth through activation of the Wnt pathway. Genome Biol. 20:84.
Liu, C. X., Li, X., Nan, F., Jiang, S., Gao, X., Guo, S. K., et al. (2019). Structure and degradation of circular RNAs regulate PKR activation in innate immunity. Cell 177, 865.e1–880.e1. doi: 10.1016/j.cell.2019.03.046
Liu, J., Yue, Y., Han, D., Wang, X., Fu, Y., Zhang, L., et al. (2014). A METTL3-METTL14 complex mediates mammalian nuclear RNA N6-adenosine methylation. Nat. Chem. Biol. 10, 93–95. doi: 10.1038/nchembio.1432
Liu, N., Dai, Q., Zheng, G., He, C., Parisien, M., and Pan, T. (2015). N(6)-methyladenosine-dependent RNA structural switches regulate RNA-protein interactions. Nature 518, 560–564. doi: 10.1038/nature14234
Liu, X., Zhang, Z., Ruan, J., Pan, Y., Magupalli, V. G., Wu, H., et al. (2016). Inflammasome-activated gasdermin D causes pyroptosis by forming membrane pores. Nature 535, 153–158. doi: 10.1038/nature18629
Martinet, W., Coornaert, I., Puylaert, P., and De Meyer, G. R. Y. (2019). Macrophage death as a pharmacological target in atherosclerosis. Front. Pharmacol. 10:306. doi: 10.3389/fphar.2019.00306
Memczak, S., Jens, M., Elefsinioti, A., Torti, F., Krueger, J., Rybak, A., et al. (2013). Circular RNAs are a large class of animal RNAs with regulatory potency. Nature 495, 333–338. doi: 10.1038/nature11928
Meng, S., Zhou, H., Feng, Z., Xu, Z., Tang, Y., Li, P., et al. (2017). CircRNA: functions and properties of a novel potential biomarker for cancer. Mol. Cancer 16:94.
Meng, X., Li, X., Zhang, P., Wang, J., Zhou, Y., and Chen, M. (2017). Circular RNA: an emerging key player in RNA world. Brief Bioinform. 18, 547–557. doi: 10.1093/bib/bbw045
Meyer, K. D., Patil, D. P., Zhou, J., Zinoviev, A., Skabkin, M. A., Elemento, O., et al. (2015). 5’ UTR m(6)A promotes cap-independent translation. Cell 163, 999–1010. doi: 10.1016/j.cell.2015.10.012
Meyer, K. D., Saletore, Y., Zumbo, P., Elemento, O., Mason, C. E., and Jaffrey, S. R. (2012). Comprehensive analysis of mRNA methylation reveals enrichment in 3’ UTRs and near stop codons. Cell 149, 1635–1646. doi: 10.1016/j.cell.2012.05.003
Nahand, J. S., Jamshidi, S., Hamblin, M. R., Mahjoubin-Tehran, M., Vosough, M., Jamali, M., et al. (2020). Circular RNAs: new epigenetic signatures in viral infections. Front. Microbiol. 11:1853. doi: 10.3389/fmicb.2020.01853
Nicolet, B. P., Engels, S., Aglialoro, F., Van Den Akker, E., Von Lindern, M., and Wolkers, M. C. (2018). Circular RNA expression in human hematopoietic cells is widespread and cell-type specific. Nucleic Acids Res. 46, 8168–8180. doi: 10.1093/nar/gky721
Niu, X., Xu, J., Liu, J., Chen, L., Qiao, X., and Zhong, M. (2020). Landscape of N(6)-Methyladenosine modification patterns in human ameloblastoma. Front. Oncol. 10:556497. doi: 10.3389/fonc.2020.556497
Pamudurti, N. R., Bartok, O., Jens, M., Ashwal-Fluss, R., Stottmeister, C., Ruhe, L., et al. (2017). Translation of CircRNAs. Mol. Cell 66, 9.e7–21.e7.
Panda, A. C., Grammatikakis, I., Munk, R., Gorospe, M., and Abdelmohsen, K. (2017). Emerging roles and context of circular RNAs. Wiley Interdiscip. Rev. RNA 8:10.1002/wrna.1386. doi: 10.1002/wrna.1386
Park, O. H., Ha, H., Lee, Y., Boo, S. H., Kwon, D. H., Song, H. K., et al. (2019). Endoribonucleolytic cleavage of m(6)A-Containing RNAs by RNase P/MRP complex. Mol. Cell 74, 494.e8–507.e8. doi: 10.1016/j.molcel.2019.02.034
Patil, D. P., Chen, C. K., Pickering, B. F., Chow, A., Jackson, C., Guttman, M., et al. (2016). m(6)A RNA methylation promotes XIST-mediated transcriptional repression. Nature 537, 369–373. doi: 10.1038/nature19342
Pendleton, K. E., Chen, B., Liu, K., Hunter, O. V., Xie, Y., Tu, B. P., et al. (2017). The U6 snRNA m(6)A Methyltransferase METTL16 regulates SAM synthetase intron retention. Cell 169, 824.e4–835.e4.
Ping, X. L., Sun, B. F., Wang, L., Xiao, W., Yang, X., Wang, W. J., et al. (2014). Mammalian WTAP is a regulatory subunit of the RNA N6-methyladenosine methyltransferase. Cell Res. 24, 177–189. doi: 10.1038/cr.2014.3
Pugliese, S. C., Poth, J. M., Fini, M. A., Olschewski, A., El Kasmi, K. C., and Stenmark, K. R. (2015). The role of inflammation in hypoxic pulmonary hypertension: from cellular mechanisms to clinical phenotypes. Am. J. Physiol. Lung Cell Mol. Physiol. 308, L229–L252.
Rao, X., Lai, L., Li, X., Wang, L., Li, A., and Yang, Q. (2021). N(6) -methyladenosine modification of circular RNA circ-ARL3 facilitates Hepatitis B virus-associated hepatocellular carcinoma via sponging miR-1305. IUBMB Life 73, 408–417. doi: 10.1002/iub.2438
Reikine, S., Nguyen, J. B., and Modis, Y. (2014). Pattern recognition and signaling mechanisms of RIG-I and MDA5. Front. Immunol. 5:342. doi: 10.3389/fimmu.2014.00342
Roundtree, I. A., Evans, M. E., Pan, T., and He, C. (2017a). Dynamic RNA modifications in gene expression regulation. Cell 169, 1187–1200. doi: 10.1016/j.cell.2017.05.045
Roundtree, I. A., Luo, G. Z., Zhang, Z., Wang, X., Zhou, T., Cui, Y., et al. (2017b). YTHDC1 mediates nuclear export of N(6)-methyladenosine methylated mRNAs. Elife 6:e31311.
Růžička, K., Zhang, M., Campilho, A., Bodi, Z., Kashif, M., Saleh, M., et al. (2017). Identification of factors required for m(6) A mRNA methylation in Arabidopsis reveals a role for the conserved E3 ubiquitin ligase HAKAI. New Phytol. 215, 157–172. doi: 10.1111/nph.14586
Rybak-Wolf, A., Stottmeister, C., Glažar, P., Jens, M., Pino, N., Giusti, S., et al. (2015). Circular RNAs in the mammalian brain are highly abundant, conserved, and dynamically expressed. Mol. Cell 58, 870–885. doi: 10.1016/j.molcel.2015.03.027
Sanger, H. L., Klotz, G., Riesner, D., Gross, H. J., and Kleinschmidt, A. K. (1976). Viroids are single-stranded covalently closed circular RNA molecules existing as highly base-paired rod-like structures. Proc. Natl. Acad. Sci. U.S.A. 73, 3852–3856. doi: 10.1073/pnas.73.11.3852
Schlee, M. (2013). Master sensors of pathogenic RNA - RIG-I like receptors. Immunobiology 218, 1322–1335. doi: 10.1016/j.imbio.2013.06.007
Sekar, D., and Lakshmanan, G. (2020). Methylation of N6-adenosine (m6A) modification in miRNAs and its implications in immunity. Epigenomics 12, 1083–1085. doi: 10.2217/epi-2020-0131
Shi, H., Wang, X., Lu, Z., Zhao, B. S., Ma, H., Hsu, P. J., et al. (2017). YTHDF3 facilitates translation and decay of N(6)-methyladenosine-modified RNA. Cell Res. 27, 315–328. doi: 10.1038/cr.2017.15
Shulman, Z., and Stern-Ginossar, N. (2020). The RNA modification N(6)-methyladenosine as a novel regulator of the immune system. Nat. Immunol. 21, 501–512. doi: 10.1038/s41590-020-0650-4
Su, H., Wang, G., Wu, L., Ma, X., Ying, K., and Zhang, R. (2020). Transcriptome-wide map of m(6)A circRNAs identified in a rat model of hypoxia mediated pulmonary hypertension. BMC Genomics 21:39. doi: 10.1186/s12864-020-6462-y
Suzuki, H., and Tsukahara, T. (2014). A view of pre-mRNA splicing from RNase R resistant RNAs. Int. J. Mol. Sci. 15, 9331–9342. doi: 10.3390/ijms15069331
Tang, C., Xie, Y., Yu, T., Liu, N., Wang, Z., Woolsey, R. J., et al. (2020). m(6)A-dependent biogenesis of circular RNAs in male germ cells. Cell Res. 30, 211–228. doi: 10.1038/s41422-020-0279-8
Thomson, D. W., and Dinger, M. E. (2016). Endogenous microRNA sponges: evidence and controversy. Nat. Rev. Genet. 17, 272–283. doi: 10.1038/nrg.2016.20
Wang, P. L., Bao, Y., Yee, M. C., Barrett, S. P., Hogan, G. J., Olsen, M. N., et al. (2014). Circular RNA is expressed across the eukaryotic tree of life. PLoS One 9:e90859. doi: 10.1371/journal.pone.0090859
Wang, X., Zhao, B. S., Roundtree, I. A., Lu, Z., Han, D., Ma, H., et al. (2015). N(6)-methyladenosine modulates messenger RNA translation efficiency. Cell 161, 1388–1399. doi: 10.1016/j.cell.2015.05.014
Wang, Y., Wang, H., Xi, F., Wang, H., Han, X., Wei, W., et al. (2020). Profiling of circular RNA N(6) -methyladenosine in moso bamboo (Phyllostachys edulis) using nanopore-based direct RNA sequencing. J. Integr. Plant Biol. 62, 1823–1838. doi: 10.1111/jipb.13002
Warda, A. S., Kretschmer, J., Hackert, P., Lenz, C., Urlaub, H., Höbartner, C., et al. (2017). Human METTL16 is a N(6)-methyladenosine (m(6)A) methyltransferase that targets pre-mRNAs and various non-coding RNAs. EMBO Rep. 18, 2004–2014. doi: 10.15252/embr.201744940
Wu, J., and Chen, Z. J. (2014). Innate immune sensing and signaling of cytosolic nucleic acids. Annu. Rev. Immunol. 32, 461–488. doi: 10.1146/annurev-immunol-032713-120156
Wu, P., Fang, X., Liu, Y., Tang, Y., Wang, W., Li, X., et al. (2021). N6-methyladenosine modification of circCUX1 confers radioresistance of hypopharyngeal squamous cell carcinoma through caspase1 pathway. Cell Death Dis. 12:298.
Wu, X., Xiao, Y., Ma, J., and Wang, A. (2020). Circular RNA: a novel potential biomarker for skin diseases. Pharmacol Res 158, 104841. doi: 10.1016/j.phrs.2020.104841
Xia, S., Feng, J., Lei, L., Hu, J., Xia, L., Wang, J., et al. (2017). Comprehensive characterization of tissue-specific circular RNAs in the human and mouse genomes. Brief Bioinform. 18, 984–992. doi: 10.1093/bib/bbw081
Xiao, W., Adhikari, S., Dahal, U., Chen, Y. S., Hao, Y. J., Sun, B. F., et al. (2016). Nuclear m(6)A reader YTHDC1 regulates mRNA splicing. Mol. Cell 61, 507–519. doi: 10.1016/j.molcel.2016.01.012
Xu, J., Wan, Z., Tang, M., Lin, Z., Jiang, S., Ji, L., et al. (2020). N(6)-methyladenosine-modified CircRNA-SORE sustains sorafenib resistance in hepatocellular carcinoma by regulating β-catenin signaling. Mol. Cancer 19:163.
Xu, Y. J., Zheng, L., Hu, Y. W., and Wang, Q. (2018). Pyroptosis and its relationship to atherosclerosis. Clin. Chim. Acta 476, 28–37. doi: 10.1016/j.cca.2017.11.005
Yan, H., Li, Y., Peng, X., Huang, D., Gui, L., and Huang, B. (2016). Resistance of mitochondrial DNA-depleted cells against oxidized low-density lipoprotein-induced macrophage pyroptosis. Mol. Med. Rep. 13, 4393–4399. doi: 10.3892/mmr.2016.5077
Yang, Y., Fan, X., Mao, M., Song, X., Wu, P., Zhang, Y., et al. (2017). Extensive translation of circular RNAs driven by N(6)-methyladenosine. Cell Res. 27, 626–641. doi: 10.1038/cr.2017.31
Yang, Y., Gao, X., Zhang, M., Yan, S., Sun, C., Xiao, F., et al. (2018). Novel role of FBXW7 circular RNA in repressing glioma tumorigenesis. J. Natl. Cancer Inst. 110, 304–315. doi: 10.1093/jnci/djx166
Yue, Y., Liu, J., Cui, X., Cao, J., Luo, G., Zhang, Z., et al. (2018). VIRMA mediates preferential m(6)A mRNA methylation in 3’UTR and near stop codon and associates with alternative polyadenylation. Cell Discov. 4:10.
Zaccara, S., Ries, R. J., and Jaffrey, S. R. (2019). Reading, writing and erasing mRNA methylation. Nat. Rev. Mol. Cell Biol. 20, 608–624. doi: 10.1038/s41580-019-0168-5
Zaphiropoulos, P. G. (1996). Circular RNAs from transcripts of the rat cytochrome P450 2C24 gene: correlation with exon skipping. Proc. Natl. Acad. Sci. U.S.A. 93, 6536–6541. doi: 10.1073/pnas.93.13.6536
Zhang, C., Wang, J., Geng, X., Tu, J., Gao, H., Li, L., et al. (2020). Circular RNA expression profile and m6A modification analysis in poorly differentiated adenocarcinoma of the stomach. Epigenomics 12, 1027–1040. doi: 10.2217/epi-2019-0153
Zhang, M., Huang, N., Yang, X., Luo, J., Yan, S., Xiao, F., et al. (2018a). A novel protein encoded by the circular form of the SHPRH gene suppresses glioma tumorigenesis. Oncogene 37, 1805–1814. doi: 10.1038/s41388-017-0019-9
Zhang, M., Zhao, K., Xu, X., Yang, Y., Yan, S., Wei, P., et al. (2018b). A peptide encoded by circular form of LINC-PINT suppresses oncogenic transcriptional elongation in glioblastoma. Nat. Commun. 9:4475.
Zhang, T. R., and Huang, W. Q. (2020). Angiogenic circular RNAs: a new landscape in cardiovascular diseases. Microvasc. Res. 129:103983. doi: 10.1016/j.mvr.2020.103983
Zhang, X. O., Wang, H. B., Zhang, Y., Lu, X., Chen, L. L., and Yang, L. (2014). Complementary sequence-mediated exon circularization. Cell 159, 134–147. doi: 10.1016/j.cell.2014.09.001
Zhang, Y., Liang, W., Zhang, P., Chen, J., Qian, H., Zhang, X., et al. (2017). Circular RNAs: emerging cancer biomarkers and targets. J. Exp. Clin. Cancer Res. 36:152.
Zhao, B. S., Roundtree, I. A., and He, C. (2017). Post-transcriptional gene regulation by mRNA modifications. Nat. Rev. Mol. Cell Biol. 18, 31–42. doi: 10.1038/nrm.2016.132
Zhao, J., Lee, E. E., Kim, J., Yang, R., Chamseddin, B., Ni, C., et al. (2019). Transforming activity of an oncoprotein-encoding circular RNA from human papillomavirus. Nat. Commun. 10:2300.
Zheng, G., Dahl, J. A., Niu, Y., Fedorcsak, P., Huang, C. M., Li, C. J., et al. (2013). ALKBH5 is a mammalian RNA demethylase that impacts RNA metabolism and mouse fertility. Mol. Cell 49, 18–29. doi: 10.1016/j.molcel.2012.10.015
Zhou, C., Molinie, B., Daneshvar, K., Pondick, J. V., Wang, J., Van Wittenberghe, N., et al. (2017). Genome-wide maps of m6A circRNAs identify widespread and cell-type-specific methylation patterns that are distinct from mRNAs. Cell Rep. 20, 2262–2276. doi: 10.1016/j.celrep.2017.08.027
Zhou, J., Wan, J., Gao, X., Zhang, X., Jaffrey, S. R., and Qian, S. B. (2015). Dynamic m(6)A mRNA methylation directs translational control of heat shock response. Nature 526, 591–594. doi: 10.1038/nature15377
Zhou, J., Wan, J., Shu, X. E., Mao, Y., Liu, X. M., Yuan, X., et al. (2018). N(6)-Methyladenosine guides mRNA alternative translation during integrated stress response. Mol. Cell 69, 636.e7–647.e7.
Zhu, Z. M., Huo, F. C., and Pei, D. S. (2020). Function and evolution of RNA N6-methyladenosine modification. Int. J. Biol. Sci. 16, 1929–1940. doi: 10.7150/ijbs.45231
Keywords: N6-methyladenosine, circular RNAs, metabolism, immunity, diseases
Citation: Wu J, Guo X, Wen Y, Huang S, Yuan X, Tang L and Sun H (2021) N6-Methyladenosine Modification Opens a New Chapter in Circular RNA Biology. Front. Cell Dev. Biol. 9:709299. doi: 10.3389/fcell.2021.709299
Received: 13 May 2021; Accepted: 06 July 2021;
Published: 23 July 2021.
Edited by:
Jing Zhang, Shanghai Jiao Tong University, ChinaReviewed by:
Xiao-Yu Chen, Shanghai Sixth People’s Hospital, ChinaPeifeng Li, Qingdao University, China
Copyright © 2021 Wu, Guo, Wen, Huang, Yuan, Tang and Sun. 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: Lijun Tang, tanglj2016@163.com; Hongyu Sun, shongyu2008@163.com
†These authors have contributed equally to this work