- 1Shanxi Key Laboratory of Artificial Intelligence Assisted Diagnosis and Treatment for Mental Disorder, First Hospital of Shanxi Medical University, Taiyuan, China
- 2Department of Psychiatry, First Hospital/First Clinical Medical College of Shanxi Medical University, Taiyuan, China
- 3National Key Disciplines, Key Laboratory for Cellular Physiology of Ministry of Education, Department of Neurobiology, Shanxi Medical University, Taiyuan, China
- 4Department of Humanities and Social Science, Shanxi Medical University, Taiyuan, China
Schizophrenia (SZ) and depression (DEP) are two common major psychiatric disorders that are associated with high risk of suicide. These disorders affect not only physical and mental health, but they also affect the social function of the individual. However, diagnoses of SZ and DEP are mainly based on symptomatic changes and the clinical experience of psychiatrists. These rather subjective measures can induce misdiagnoses and missed diagnoses. Therefore, it is necessary to further explore objective indexes for improving the early diagnoses and prognoses of SZ and DEP. Current research indicates that non-coding RNA (ncRNA) may play a role in the occurrence and development of SZ and DEP. Circular RNA (circRNA), as an important component of ncRNA, is associated with many biological functions, especially post-transcriptional regulation. Since circRNA is easily detected in peripheral blood and has a high degree of spatiotemporal tissue specificity and stability, these attributes provide us with a new idea to further explore the potential value for the diagnosis and treatment of SZ and DEP. Here, we summarize the classification, characteristics, and biological functions of circRNA and the most significant results of experimental studies, aiming to highlight the involvement of circRNA in SZ and DEP.
Introduction
Among the great quantity molecular regulatory factors affecting gene expression, non-coding RNA (ncRNA) is known to be a class of important regulatory factors. NcRNA refers to a kind of RNAs that do not encode proteins, but can participate in genetic regulatory processes at multiple levels, including interactions with DNA, RNA, and proteins through a variety of mechanisms (1). Gene regulation of ncRNAs via epigenetics plays an important role in many major diseases, such as schizophrenia (SZ) and depression (DEP) (2–4). At present, SZ is considered to represent an umbrella of genetic diseases that are mainly characterized by a certain number of genetic-interaction networks and their differential clinical symptoms (5). For DEP, studies have proposed a pathogenic model that involves interactions among the environment, genetics, and epigenetics (6). Moreover, the high rate of comorbidity observed in psychiatric disorders indicates a large genetic correlation between each other, and the existence of shared molecular mechanisms involved in the occurrence and development of diseases (7). Sha Liu et al.’s genome-wide association study (GWAS) of seven psychiatric disorders (8), including SZ, DEP, bipolar disorder (BD), autism spectrum disorder (ASD), attention-deficit/hyperactivity disorder (ADHD), anxiety disorder (ANX), and neuroticism, also confirmed this. And they further found that chromosomes 5q14.3, 11q23.2, and 7p22.3 are the three genomic regions with the highest pleiotropic effects, suggesting that genetic factors have significant influence on SZ and DEP. Although few in-depth studies have been published on specific signaling pathways and networks regulated by ncRNAs in SZ and DEP, especially in terms of circular RNAs (circRNAs), it is possible that further elucidation of ncRNAs may be useful for better understanding the pathogenesis of SZ and DEP.
CircRNA, as an important type of ncRNA, forms an annular closed structure by undergoing reverse splicing, which was first discovered in plants in the 1970s (9). Since it was originally found only in a few transcriptional genes, circRNAs were considered to be a product of splicing errors from low-abundance transcriptional processes. However, this hypothesis has since been refuted. CircRNAs are widely present in human cells, some of which are expressed at as much as 10-fold higher levels compared to those of their linear isomer counterparts (10–12). CircRNA differs from linear RNA in that it is directly joined together to form a covalently closed structure without a 3′ tail and 5′ cap element (13–15). Due to the circular structure of circRNA, it is not easily degraded by RNA enzymes and also exhibits a high stability within cells (16, 17), enabling it to stably perform cellular functions over a longer period of time.
At present, the diagnostic criteria for SZ and DEP are still based solely on symptomology. Although diagnostic criteria—such as from the Diagnostic and Statistical Manual of Mental Disorders (Fifth Edition) (DSM-5) of American Psychiatric Association and the International Classification of Diseases (ICD-11)—are currently available, diagnoses and evaluations of treatment efficacies for SZ and DEP still primarily rely on clinical meetings with psychiatrists and lack objective physiological, biochemical, and pathological indicators. Hence, this more subjective diagnostic system can easily cause misdiagnoses and missed diagnoses. The above characteristics of circRNAs provide us with a novel approach to explore the biological basis for the occurrence and development of SZ and DEP. Therefore, the focus of this review is on circRNAs and their possible involvement in both SZ and DEP.
Methods
We searched PubMed and GeenMedical for publications over the period of 1976 to 2019 and included the following the key phrases during searching: “circular RNA and neuropsychiatric disorders” “circular RNA and schizophrenia” “circular RNA and depression” and “circular RNA and biological function”. Only papers in English were used in the preparation of this review. The references obtained were screened by the authors to determine which ones would be included for discussion in the present review.
Overview of circRNAs
Classifications and Looping Mechanisms of circRNAs
CircRNAs are divided into the following four categories: (1) whole-exonic-type circRNAs (EcircRNA); (2) circRNAs with introns and exons (EIcirRNA) (18); (3) lasso-type circRNAs composed of introns (ciRNA); and (4) circRNAs produced by cyclization of viral RNA genomes, ribosomal RNAs, small-nuclear RNAs, and transfer RNAs (tricRNA) (Figure 1). Based on their multiple biological mechanisms, circRNAs can retain the ability to either stay inside the nucleus or undergo nuclear export. This selective retention is important for determining the biological roles of different circRNAs and further provides us with research ideas for exploring the potential of circRNAs in the occurrence and development of SZ, DEP, and many other diseases. CircRNA is generally considered to be the result of reverse splicing. This indicates that the downstream 5′ donor site is covalently cyclized with the upstream 3′ receptor site (19). At present, there are three main looping mechanisms of circRNAs.
Figure 1 The overview diagram shows four biosynthetic pathways related to the source of circular RNAs: (1) whole-exonic-type circRNAs (EcircRNA); (2) circRNAs with introns and exons (EIcirRNA); (3) lasso-type circRNAs composed of introns (ciRNA); and (4) circRNAs produced by cyclization of viral RNA genomes, ribosomal RNAs, small-nuclear RNAs, and transfer RNAs (tricRNA).
Cable-Tail-Insertion Cyclization Dependent on a Shearing Body
In eukaryotic cells, the formation of circRNAs by exon cable-tail-insertion cyclization (20, 21) is the earliest and most common way of looping, and more than 80% of circRNAs contain exons. Exon cable-tail-insertion cyclization relies on a typical shearing-body mechanism. On the mRNA precursor, the downstream 5′ donor site of the exon is catalyzed to link to the upstream 3′ receptor site by successively assembling small-nuclear ribosomal proteins to form cable-tail-insertion cyclization. Eventually, RNA is formed by further shearing. However, the specific mechanism underlying this process remains unclear.
Promotion of circRNA Formation by Cyclic-Acting Elements
Partial introns on both sides of circRNA exons contain reverse complementary sequences that can form RNA double chains side by side at shear sites (22). Finally, two kinds of different circRNAs with or without introns are formed by variable shearing. In general, a nucleotide sequence of 30–40 nt in length can promote the cyclization of circRNAs. Furthermore, exons inside can also compete for RNA pairing and ultimately form different types of circRNAs via variable shearing.
Regulation of circRNA Formation by RNA-Binding Proteins (RBPs)
In Drosophila melanogaster, MBL protein (23) can promote the formation of circRNAs by binding to the introns flanking an exon. Some studies (24) have shown that hundreds of circRNAs are regulated during the process of the epithelial-mesenchymal transition (EMT), and the looping formations of more than one-third of circRNAs are dynamically regulated by the splicing factor, Quaking (QKI). The effect of QKI on circRNA abundance depends on the binding motif of the QKI intron. Critically, the addition of a QKI motif is sufficient to re-induce the formation of circRNA from a linear splicing transcript. In contrast to the effect of QKI protein, high expression of ADAR1 protein (22) can inhibit the formation of circRNAs by destroying the RNA pairing flanking the exon.
Some exon-derived circRNAs are localized in the cytoplasm and contain both ribosome-binding sites and evolutionarily conserved termination codons. Exon-derived circRNAs also share the start codon of the host mRNA, which indicates that hundreds of endogenous circRNAs have translation potential (25), and some of their translation products can carry out functions. Previous studies have confirmed that many (31 out of 132) proteins encoded in ribosomal-related circRNAs (ribo-circRNAs) contain at least one identifiable protein domain. Furthermore, circRNA, such as circMbl, can be translated in a cap-independent manner in vitro (26). The same phenomenon is found for circ-ZNF609, which is associated with multiple ribosomes. Circ-ZNF609 translates protein in a splicing-dependent and cap-independent manner. This gene controls the proliferation of primary myoblasts and is the key regulatory factor of their growth (27).
Molecular Biological Functions of circRNAs
Efficient microRNA Sponges
One of the most important molecular biological functions of circRNAs is that they indirectly affect the expression of target gene mRNAs by binding to microRNAs (miRNAs). As an important regulator at the post-transcriptional level, the function of miRNAs is mainly to downregulate or hinder the translation of target genes (mRNAs) by binding to the complementary 3′-end non-coding region located on the target mRNA (28, 29). As miRNA sponges, most circRNAs play an important role in transcriptional and translational regulation of mRNAs by competing for intracellular miRNAs (30, 31).
Among these miRNA sponges, Cdr1as circRNA is highly expressed and stable. It contains more than 70 conserved miRNA response elements (MREs) of miR-7 and one MRE of miR-671, which are then widely combined with AGO2, a key element of gene silencing. Bezzi and colleagues (32) provided evidence to support the functional association among Cdr1as, miR-7, and its target gene. A plasmid expressing Cdr1as was injected into zebrafish embryos to knockdown the expression of miR-7. This effect was partially rescued by injecting miR-7 precursor. At the same time, miR-671 indirectly influenced the expression level of the miR-7 by triggering the degradation of Cdr1as. The interaction between circRNAs and miRNAs plays an important role in tumors, such as breast cancer, liver cancer and cervical cancer (31–35). Han et al. found that circMTO1 was significantly downregulated in human hepatocellular carcinoma (HCC) (36). Silencing of circMTO1, a miRNA-9 sponge, in patients with HCC downregulated the expression of the miRNA9 target gene, p21, thereby promoting the proliferation and invasion of HCC cells. Furthermore, the survival time of HCC cells with low expression of circMTo1 was reduced. This study also found that miRNA9 inhibitors could block the silencing effect of circMTO1. These findings suggest that cycMTO1 may be used as an objective evaluation index of survival prognosis of HCC and may be a potential target for liver cancer treatment. In addition to circMTO1 acting as a sponge for miRNA-9 in HCC, hsa-circRNA-103809 (37) was also found to participate in the pathological mechanism of HCC. It could be used as a sponge for miR-377-3p to increase the expression of its target gene [fibroblast growth factor receptor 1 (FGFR1)] and promoted the proliferation and invasion of HCC cells. Moreover, hsa-circRNA-103809 short-hairpin RNA could act as a tumor inhibitor by downregulating FGFR1 in HCC.
Templates for Translating Proteins
N6-methyladenosine (M6A) is the most abundant base modification of RNA. A previous study found (38) that the m6A motif was enriched in circRNAs, and that a single m6A site was sufficient to drive the initiation of circRNA protein translation in human cells. However, this protein translation of circRNA driven by M6A did not function without the following: the activation factor, eIF4G2; the M6A reader, YTHDF 3; enhancement of the methyltransferase, METTL 3/14; inhibition and upregulation of demethylase during heat shock, FTO. Interestingly, further analysis indicated that m6A-driven circRNA translation is widespread.
Interactions With RBPs
Interactions between circRNAs and RBPs are not only related to their own biological occurrences, but they are also involved in the post-transcriptional regulation of RNAs (39). In general, the cell-cycle protein, cyclin-dependent kinase 2 [also known as cell-division protein kinase 2 (CDK2)], interacts with the cell-cycle proteins A and E to facilitate the cell cycle. However, cyclin-dependent kinase inhibitor 1 (p21) can inhibit these interactions and prevent cell-cycle progression. Du et al. (40) found that circ-Foxo3 was highly expressed in non-cancer cells, forming a ternary complex to inhibit the cell cycle by binding to CDK2 and p21. Additionally, silencing of endogenous circ-Foxo3 promoted cell proliferation. Meanwhile, Zhu et al. (41) showed that circZKSCAN1 exerted its inhibitory effect by competitively binding to the RNA-binding protein (RBP), FMRP. This mechanism blocked the binding of FMRP to the β-catenin-binding protein cell cycle and apoptosis regulator 1 (CCAR1) mRNA, and subsequently inhibited the transcriptional activity of the Wnt signal transduction pathway. In addition, downregulation of the RNA-splicing protein, Quaking5, in HCC tissue caused a decrease of circZKSCAN1. These findings suggest the formation of a Qki5-circZKSCAN1-FMRP-CCAR1-Wnt signaling pathway, which may be useful as a potential therapeutic target for HCC (Figure 2).
The above results indicate that circRNAs act not only as miRNA sponges by competing with intracellular miRNA to influence the expression of target genes, but they also act as translational templates for proteins. More importantly, circRNAs lead to altered gene expression outcomes by splicing their corresponding precursor transcripts, which may play an important role in the occurrence and development of various diseases (Table 1).
CircRNAs and SZ
SZ is one of the most common mental disorders, with a global prevalence of 1%, which brings a heavy burden to affected families and society (42). According to the latest research conducted by Huang et al. (43), the prevalence of SZ is 0.559% in China. Symptoms of SZ are divided into positive, negative, and cognitive categories (44, 45). The diagnosis of SZ is mainly based on symptomatic changes and the clinical experience of psychiatrists. However, these rather subjective measures can cause misdiagnoses and missed diagnoses. Therefore, it is necessary to further explore objective indexes to improve the early diagnosis and prognosis of SZ. Since circRNAs are easily detected in peripheral blood and have a high degree of spatiotemporal tissue specificity and stability (46, 47), these properties provide us with a novel approach to further explore the mechanisms underlying the occurrence and development of SZ.
Molecular Mechanisms of circRNAs in SZ
As representative antisense nucleic-acid sponge sequences, circRNAs can competitively bind to miRNAs via MREs, which negatively influence the regulatory effects of miRNAs at the post-transcriptional level. MiRNA-320a-3p, miRNA-320b, miRNA-181b-5p, 21-5p, 195-5p, 137, 346, 34a-5p, and hsa-miRNA-206 (48–50) have all been confirmed to play important roles in the occurrence and development of SZ. Wei et al. (51) found that eight miRNAs were upregulated in patients with SZ, as compared to corresponding levels in control participants. Through quantitative reverse-transcription polymerase chain reaction (qRT-PCR) assays and a follow-up study on 400 patients with SZ who received regular atypical antipsychotic treatment for 12 months, this study ultimately showed that miR-130b and miR-193a-3p were upregulated in SZ, and the possibility of time-dependent changes were ruled out. Hence, miR-130b and miR-193a-3p may be used as state-independent biomarkers for SZ. A recent study (52) has confirmed that the expression of the network axis composed of miR-30a-5p, its transcription factor, EGR1 and target gene, NEUROD1, changed with the disease state in the patients with SZ before and after treatment. Further receiver operating characteristic (ROC) curve analysis revealed that compared with the single miR-30a-5p, the EGR1-miR-30a-5p-NEUROD1 molecular regulatory axis had higher diagnostic value for predicting SZ. Hence, these findings suggest a possibility for the miRNA-TF-gene regulatory network/axis as a novel diagnostic marker for SZ.
Due to their circular structure, circRNAs are highly conserved across species. Some studies have confirmed that circRNAs are highly expressed in the brain and are differentially expressed in distinct tissues and disease states (53, 54). Moreover, circRNAs exhibit dynamic expression levels during neural development, including neural differentiation and maturation. For instance, circRNAs are upregulated during neuronal differentiation, and compared with their mRNA homologues, they generally exhibit differential expression (55–57). SZ often occurs in late adolescence and early adulthood. Previous studies have suggested that excessive pruning of synapses in the brain during growth and development may contribute to the pathogenesis of SZ (58).
Some studies have suggested synaptic localization of circRNAs. Wolf et al. (56) found that circRNAs in the nervous system were generally upregulated compared to those in the thyroid, liver, and muscle, and that expression varied among different brain regions. Especially in the cerebellum, the expression ratio of circular to linear RNAs was significantly increased, which correlated with a larger number of neurons in the cerebellum compared to that in other brain regions. Furthermore, by measuring high-purity synaptosomal components, this study found that brain-expressed circRNAs had a strong enriching effect at synapses. Another study found that circStau2a was mainly located at synapses, whereas its linear transcript, mRNAStau2, was almost completely localized to the cytoplasm (56). Similarly, circRNA derived from the known neuronal differentiation regulator, RMST, exhibits a high synaptic enrichment rate (59). It has also been found that during the development of the hippocampus, upregulated circRNAs are produced by the gene locus that simultaneously encodes proteins rich in synapse-related functions. Conversely, no enrichment of encoding of any other functional class of proteins was found in the gene locus that produced a downregulated dynamic expression pattern of circRNAs (29, 60, 61). In summary, the expression levels of circRNAs are regulated by neuronal development. Many circRNAs also change the structures of neurons, which are largely independent of the function of their linear transcripts.
Expression of circRNAs in Postmortem SZ Brains
Our literature review showed that the most extensive studies of circRNAs and SZ have involved the analysis of the expression profiles of circRNAs based on control studies of patients with SZ and the healthy controls matched to them. Mahmoudi et al. used circRNA enrichment sequencing to analyze the expression profiles of circRNAs in the cerebral cortex (BA46) of 35 postmortem patients with SZ as well as those of healthy controls (62). They found that more than 95,000 types of circRNAs in the human dorsolateral prefrontal cortex (DLPFC) showed significant diversity, and that half of them had not previously been reported, such as circHomer1, circKlhl2, circMpped1, and circNell2 (63, 64). A differential analysis of these circRNAs revealed that 390 circRNAs were downregulated and 240 were completely deleted in patients with SZ. In addition to the reduced overall level of gene expression in circRNAs, the total numbers of unique circRNAs found in SZ were also reduced compared to those in heathy controls, and there was a large overlap between them. For instance, circNEll2, which is upregulated after synaptic plasticity induction, was downregulated in patients with SZ. The researchers also explored the correlation between differentially expressed circRNAs and their linear isoforms. Surprisingly, the expression levels of more than half of the circRNAs showed an inverse relationship with those of their linear RNAs, suggesting that both forms of transcripts may sometimes antagonize each other’s biosynthesis. Zimmerman et al. (65) also discovered circHomer1a, which is highly expressed and neuron-enriched in the frontal cortex. Interestingly, in the prefrontal cortex (PFC) and induced pluripotent stem cell-derived neuronal cultures in patients with SZ and BD, circHomer1a was found to be significantly downregulated. Additionally, changes in circHomer1a expression in the DLPFC and orbitofrontal cortex (OFC) were positively correlated with the age of onset of SZ. In addition, circHomer1a was found to interact with RNA-binding protein, HuD, through animal-level verification, which was involved in regulating synaptic gene expression and cognitive flexibility. These findings are of great significance for exploring the molecular mechanisms underlying the pathogenesis of mental disorders.
The above results suggest that the overall expression levels of circRNAs are downregulated in patients with SZ, while there are still some unique and significantly differentially expressed circRNAs that differ from this overall trend. These circRNAs may play an important role in the clinical phenotypes of SZ by regulating corresponding cellular metabolic pathways underlying this mental disorder. Furthermore, there is a strong possibility of abnormal regulation of circRNAs in other mental illnesses and behavioral disorders.
Expression of circRNAs in Peripheral Blood of Patients With SZ
CircRNA plays a particularly important role in the regulation of gene expression in SZ (62, 66, 67). Previous studies (68) have identified that differentially expressed circRNAs from peripheral blood perhaps contribute on the development of the diagnosis and treatment of SZ.
A case-control study carried out by Yao and colleagues (69) screened nine kinds of significantly expressed circRNAs. Further verification by qRT-PCR revealed that three circRNAs were downregulated and two circRNAs were upregulated in the SZ group. The expression levels of circRNAs in patients with SZ after 4 and 8 weeks of conventional antipsychotic treatment were then re-quantified to see if there was any change. ROC curve analysis showed that three circRNAs (hsa circRNA 103704, hsa circRNA 101836, and hsa circRNA 104597) were of diagnostic significance, with a sensitivity of 87.25% and specificity of 85.44%. In particular, the sensitivity and specificity of hsa circRNA 104597 were 84.31% and 86.41%, respectively, indicating that it may be useful in the diagnostic and therapeutic evaluation of SZ.
The above results suggest that circRNAs may provide novel ideas for exploring objective biomarkers for SZ diagnosis. However, there has only been one study on this topic, which had some limitations. For instance, the sample size of screening differentially expressed circRNAs was relatively small. Since patients with SZ had been treated with standardized antipsychotic drugs for a relatively short period of time, changes in the expression level of circRNAs before and after treatment may also be due to their own time-dependent changes. Furthermore, only one type of circRNA was found to may represent a novel objective diagnostic value for SZ. Finally, there are still many significantly expressed circRNAs related to SZ that require further exploration (69–71).
CircRNAs and Depression
Depression (DEP) is one of the most common mental disorders in humans. The World Health Organization (WHO) has reported that more than 800 million people worldwide suffer from DEP (72, 73). Approximately 40–60% of patients with DEP show suicidal ideation or behavior, and the suicide rate has been reported to be as high as 15% (74, 75). The disease burden of DEP is increasing year by year. In 1990, DEP ranked fifth in terms of its global burden as a disease. However, DEP is expected to become the top disease burden in China by 2030 (76). The process of diagnosing this disease is similar to that for SZ and, as with SZ, it can easily lead to misdiagnoses and missed diagnoses. Therefore, it is necessary to further explore the pathological mechanisms and specific molecular markers related to DEP to improve the early diagnosis and prognosis of DEP.
Expression of circRNAs in Peripheral Blood in Depression
Cui and colleagues randomly selected peripheral blood mononuclear cells (PBMCs) from five patients with DEP and five healthy controls to analyze the expression profile of circRNAs. First, 15 differentially expressed circRNAs (including four upregulated and 11 downregulated circRNAs) were screened (77). Then, four significantly differentially expressed circRNAs were identified through verification in 100 DEP patients, 103 healthy controls, and 30 randomly selected DEP patients after 4 and 8 weeks of antidepressant therapy. However, only hsa-circRNA-103636 showed significant changes in DEP patients after 8 weeks of antidepressant treatment. Further ROC curve analysis showed that the sensitivity and specificity of this gene were 0.73 and 0.65, respectively. Therefore, hsa-circRNA-103636 may have potential diagnostic value of DEP.
Many studies have shown that type-2 diabetes (T2DM) is closely related to the occurrence and development of DEP (78, 79). Tian et al. collected venous plasma from five patients with T2DM and five cases of T2DM with DEP to determine the expression profiles of ncRNAs (including lncRNAs, circRNAs, and mRNAs) in order to discover specific molecular markers associated with DEP. The results of screening identified 28 lncRNAs, 107 circRNAs (including 75 upregulated and 32 downregulated circRNAs), and 89 mRNAs in the differential expression profiles of patients with DEP (80). Compared with that in the T2DM group, circRNA-TFRC had a significantly higher expression level in the DEP group. This circRNA is related to insulin resistance, which indicates that patients with DEP may have more severe insulin resistance symptoms. Some studies (81) have found that the expression level of TFRC is associated with mental disorders, suggesting that TFRC may be useful as a molecular target for studying the pathogenesis of DEP. Another upregulated circRNA-TNIK [Traf 2- and Nck-interacting kinase (TNIK)] is also associated with DEP. TNIK is a serine/threonine kinase that exhibits a high expression level in the brain, which has a good effect on the development and synaptic transmission of dendritic cells. A previous study (80) found that compared with a T2DM group, the expression level of TNIK was significantly increased in the DEP group. This suggests that the function of TNIK merits further investigation.
Jiang and colleagues (82) used microarrays to analyze plasma samples from seven patients with T2DM and seven patients with T2DM and DEP. Compared with the T2DM group, 183 circRNAs were upregulated and 64 circRNAs were downregulated in T2DM patients with DEP. Among them, the upregulated gene, DCP2 (hsa-circRNA-001520), may be used as a type of de-capping enzyme and plays an important role in neuronal development and mental retardation by degrading mRNAs. Another downregulated gene, CSGALNACT1 (hsa-circRNA-001781), is associated with an antidepressant response (83). In addition, Zhang et al. demonstrated that circDYM was also significantly decreased in MDD patients and two depressive-like mouse models [induced via chronic unpredictable pressure (CUS) and lipopolysaccharide (LPS)]. Recovery of circDYM levels significantly ameliorated CUS- or LPS-induced depressive-like behavior, which suggests that circDYM may be a novel therapeutic target for DEP (84).
According to the target genes of circRNAs related to DEP predicted by TargetScan and Miranda, 18 miRNAs and 529 mRNAs were found to interact with 4 circRNAs. Interestingly, hsa-miR-761 and hsa-miR-298 are common targets for hsa-circRNA-003251, hsa-circRNA-005019, and hsa-circRNA-015115. A previous study found that hsa-circRNA-003251 and hsa-circRNA-015115 may play important roles in the circRNA-miRNA-mRNA network associated with DEP by acting as an hsa-miR-761 sponge. Hsa-miR-761 (85, 86) has been shown to be involved in the regulation of mitochondrial networks and to promote the generation of learning and memory. Furthermore, overexpression of hsa-miR-761 inhibits the p38-MAPK signaling pathway (87). Therefore, hsa-miR-761 may participate in the pathological mechanisms of DEP. In addition, KEGG pathway analysis predicted that upregulation of hsa-circRNA-003251, hsa-circRNA-015115, hsa-circRNA-100918, and hsa-circRNA-001520 might be involved in the Wnt signaling pathway implicated in thyroid hormones. Another study showed that thyroid hormones affected mood and promoted signal pathways in the brain (88). These findings suggest that the occurrence and development of DEP might also be related to changes and metabolic disorders involving thyroid hormones (Table 2).
Based on the above findings, circRNA is a kind of biological signal that is measurable in physiological, pathological, and therapeutic process and perhaps contribute on predicting risk assessments, early diagnoses, and monitoring of responses in patients with SZ or DEP. At present, although some studies on the specific expression levels of circRNAs in various diseases have emerged, most of them have been at the level of expression profiles. Therefore, further studies are needed to better elucidate the precise associations between circRNAs and clinical symptoms in various diseases.
Limitations and Future Research
At present, an increasing number of studies on circRNAs related to SZ and DEP are emerging. However, there are still many questions to be answered, such as the following: whether changes in a single circRNA or among a group of circRNAs are specific to SZ and DEP; whether levels of circRNAs change with age; whether the expression of circRNAs are affected by smoking or drug abuse; and whether the circRNAs in cellular compartments have any differences in terms of their molecular bases of SZ and DEP (e.g., extra-vesicular circRNAs). These limitations indicate fruitful directions for future studies that may help elucidate the potential diagnostic and therapeutic value of circRNAs for SZ and DEP.
Author Contributions
YX designed and supervised the study. ZL and SL drafted the manuscript. XL, WZ, and JL collected some literature.
Funding
This work was supported by the National Natural Science Foundation of China, (81701326 and 81971601), the Support Program of the Youth Sanjin Scholars, and the 136 Medical Rejuvenation Project of Shanxi Province.
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.
Acknowledgments
The authors are grateful to the National Natural Science Foundation of China for research funding.
References
1. Stefanska B, MacEwan DJ. Epigenetics and pharmacology. Br J Pharmacol (2015) 172(11):2701–4. doi: 10.1111/bph.13136
2. Wang Q, Roy B, Dwivedi Y. Co-expression network modeling identifies key long non-coding RNA and mRNA modules in altering molecular phenotype to develop stress-induced depression in rats. Transl Psychiatry (2019) 9(1):125. doi: 10.1038/s41398-019-0448-z
3. Nestler EJ, Peña CJ, Kundakovic M, Mitchell A, Akbarian S. Epigenetic Basis of Mental Illness. Neuroscientist (2016) 22(5):447–63. doi: 10.1177/1073858415608147
4. Perkins DO, Jeffries C, Sullivan P. Expanding the ‘central dogma’: the regulatory role of nonprotein coding genes and implications for the genetic liability to schizophrenia. Mol Psychiatry (2005) 10(1):69–78. doi: 10.1038/sj.mp.4001577
5. Arnedo J, Svrakic DM, Del Val C, Romero-Zaliz R, Hernández-Cuervo H, Fanous AH, et al. Uncovering the hidden risk architecture of the schizophrenias: confirmation in three independent genome-wide association studies. Am J Psychiatry (2015) 172(2):139–53. doi: 10.1176/appi.ajp.2014.14040435
6. Mill J, Petronis A. Molecular studies of major depressive disorder: the epigenetic perspective. Mol Psychiatry (2007) 12(9):799–814. doi: 10.1038/sj.mp.4001992
7. Kessler RC, Chiu WT, Demler O, Merikangas KR, Walters EE. Prevalence, severity, and comorbidity of 12-month DSM-IV disorders in the National Comorbidity Survey Replication. Arch Gen Psychiatry (2005) 62(6):617–27. doi: 10.1001/archpsyc.62.6.617
8. Liu S, Rao S, Xu Y, Li J, Huang H, Zhang X, et al. Identifying common genome-wide risk genes for major psychiatric traits. Hum Genet (2020) 139(2):185–98. doi: 10.1007/s00439-019-02096-4
9. Sanger HL, Klotz G, Riesner D, Gross HJ, Kleinschmidt AK. Viroids are single-stranded covalently closed circular RNA molecules existing as highly base-paired rod-like structures. Proc Natl Acad Sci U S A (1976) 73(11):3852–6. doi: 10.1073/pnas.73.11.3852
10. Jeck WR, Sorrentino JA, Wang K, Slevin MK, Burd CE, Liu J, et al. Circular RNAs are abundant, conserved, and associated with ALU repeats. RNA (2013) 19(2):141–57. doi: 10.1261/rna.035667.112
11. Wang PL, Bao Y, Yee MC, Barrett SP, Hogan GJ, Olsen MN, et al. Circular RNA is expressed across the eukaryotic tree of life. PloS One (2014) 9(6):e90859. doi: 10.1371/journal.pone.0090859
12. Jeck WR, Sharpless NE. Detecting and characterizing circular RNAs. Nat Biotechnol (2014) 32(5):453–61. doi: 10.1038/nbt.2890
13. Capel B, Swain A, Nicolis S, Hacker A, Walter M, Koopman P, et al. Circular transcripts of the testis-determining gene Sry in adult mouse testis. Cell (1993) 73(5):1019–30. doi: 10.1016/0092-8674(93)90279-Y
14. Danan M, Schwartz S, Edelheit S, Sorek R. Transcriptome-wide discovery of circular RNAs in Archaea. Nucleic Acids Res (2012) 40(7):3131–42. doi: 10.1093/nar/gkr1009
15. Puttaraju M, Been MD. Group I permuted intron-exon (PIE) sequences self-splice to produce circular exons. Nucleic Acids Res (1992) 20(20):5357–64. doi: 10.1093/nar/20.20.5357
16. Salzman J, Gawad C, Wang PL, Lacayo N, Brown PO. Circular RNAs are the predominant transcript isoform from hundreds of human genes in diverse cell types. PloS One (2012) 7(2):e30733. doi: 10.1371/journal.pone.0030733
17. Vincent HA, Deutscher MP. Substrate recognition and catalysis by the exoribonuclease RNase R. J Biol Chem (2006) 281(40):29769–75. doi: 10.1074/jbc.M606744200
18. Li Z, Huang C, Bao C, et al. Exon-intron circular RNAs regulate transcription in the nucleus. Nat Struct Mol Biol (2015) 22(3):256–64. doi: 10.1038/nsmb.2959
19. Wilusz JE, Sharp PA. Molecular biology. A circuitous route to noncoding RNA. Science (2013) 340(6131):440–1. doi: 10.1126/science.1238522
20. Lu T, Cui L, Zhou Y, Zhu C, Fan D, Gong H, et al. Transcriptome-wide investigation of circular RNAs in rice. RNA (2015) 21(12):2076–87. doi: 10.1261/rna.052282.115
21. Barrett SP, Wang PL, Salzman J. Circular RNA biogenesis can proceed through an exon-containing lariat precursor. Elife (2015) 4:e07540. doi: 10.7554/eLife.07540
22. Ivanov A, Memczak S, Wyler E, Torti F, Porath HT, Orejuela MR, et al. Analysis of intron sequences reveals hallmarks of circular RNA biogenesis in animals. Cell Rep (2015) 10(2):170–7. doi: 10.1016/j.celrep.2014.12.019
23. Ashwal-Fluss R, Meyer M, Pamudurti NR, Ivanov A, Bartok O, Hanan M, et al. circRNA biogenesis competes with pre-mRNA splicing. Mol Cell (2014) 56(1):55–66. doi: 10.1016/j.molcel.2014.08.019
24. Conn SJ, Pillman KA, Toubia J, Conn VM, Salmanidis M, Phillips CA, et al. The RNA binding protein quaking regulates formation of circRNAs. Cell (2015) 160(6):1125–34. doi: 10.1016/j.cell.2015.02.014
25. Chen CY, Sarnow P. Initiation of protein synthesis by the eukaryotic translational apparatus on circular RNAs. Science (1995) 268(5209):415–7. doi: 10.1126/science.7536344
26. Pamudurti NR, Bartok O, Jens M, Ashwal-Fluss R, Stottmeister C, Ruhe L, et al. Translation of CircRNAs. Mol Cell (2017) 66(1):9–21.e7. doi: 10.1016/j.molcel.2017.02.021
27. Legnini I, Di Timoteo G, Rossi F, Morlando M, Briganti F, Sthandier O, et al. Circ-ZNF609 Is a Circular RNA that Can Be Translated and Functions in Myogenesis. Mol Cell (2017) 66(1):22–37.e9. doi: 10.1016/j.molcel.2017.02.017
28. Kim VN. MicroRNA biogenesis: coordinated cropping and dicing. Nat Rev Mol Cell Biol (2005) 6(5):376–85. doi: 10.1038/nrm1644
29. Bartel DP. MicroRNAs: target recognition and regulatory functions. Cell (2009) 136(2):215–33. doi: 10.1016/j.cell.2009.01.002
30. Salmena L, Poliseno L, Tay Y, Kats L, Pandolfi PP. A ceRNA hypothesis: the Rosetta Stone of a hidden RNA language. Cell (2011) 146(3):353–8. doi: 10.1016/j.cell.2011.07.014
31. Memczak S, Jens M, Elefsinioti A, Torti F, Krueger J, Rybak A, et al. Circular RNAs are a large class of animal RNAs with regulatory potency. Nature (2013) 495(7441):333–8. doi: 10.1038/nature11928
32. Bezzi M, Guarnerio J, Pandolfi PP. A circular twist on microRNA regulation. Cell Res (2017) 27(12):1401–2. doi: 10.1038/cr.2017.136
33. Hansen TB, Jensen TI, Clausen BH, Bramsen JB, Finsen B, Damgaard CK, et al. Natural RNA circles function as efficient microRNA sponges. Nature (2013) 495(7441):384–8. doi: 10.1038/nature11993
34. Tang W, Ji M, He G, Yang L, Niu Z, Jian M, et al. Silencing CDR1as inhibits colorectal cancer progression through regulating microRNA-7. Onco Targets Ther (2017) 10:2045–56. doi: 10.2147/OTT.S131597
35. Piwecka M, Glažar P, Hernandez-Miranda LR, Memczak S, Wolf SA, Rybak-Wolf A. Loss of a mammalian circular RNA locus causes miRNA deregulation and affects brain function. Science (2017) 357(6357). doi: 10.1126/science.aam8526
36. Han D, Li J, Wang H, Su X, Hou J, Gu Y, et al. Circular RNA circMTO1 acts as the sponge of microRNA-9 to suppress hepatocellular carcinoma progression. Hepatology (2017) 66(4):1151–64. doi: 10.1002/hep.29270
37. Zhan W, Liao X, Chen Z, et al. Circular RNA hsa_circRNA_103809 promoted hepatocellular carcinoma development by regulating miR-377-3p/FGFR1/ERK axis. J Cell Physiol (2019) 235(2):1733–45. doi: 10.1002/jcp.29092
38. Yang Y, Fan X, Mao M, Song X, Wu P, Zhang Y, et al. Extensive translation of circular RNAs driven by N6-methyladenosine. Cell Res (2017) 27(5):626–41. doi: 10.1038/cr.2017.31
39. Janga SC, Mittal N. Construction, structure and dynamics of post-transcriptional regulatory network directed by RNA-binding proteins. Adv Exp Med Biol (2011) 722:103–17. doi: 10.1007/978-1-4614-0332-6_7
40. Du WW, Yang W, Liu E, Yang Z, Dhaliwal P, Yang BB. Foxo3 circular RNA retards cell cycle progression via forming ternary complexes with p21 and CDK2. Nucleic Acids Res (2016) 44(6):2846–58. doi: 10.1093/nar/gkw027
41. Zhu YJ, Zheng B, Luo GJ, Ma XK, Lu XY, Lin XM, et al. Circular RNAs negatively regulate cancer stem cells by physically binding FMRP against CCAR1 complex in hepatocellular carcinoma. Theranostics (2019) 9(12):3526–40. doi: 10.7150/thno.32796
42. Cohen CI, Meesters PD, Zhao J. New perspectives on schizophrenia in later life: implications for treatment, policy, and research. Lancet Psychiatry (2015) 2(4):340–50. doi: 10.1016/S2215-0366(15)00003-6
43. Huang Y, Wang Y, Wang H, Liu Z, Yu X, Yan J, et al. Prevalence of mental disorders in China: a cross-sectional epidemiological study. Lancet Psychiatry (2019) 6(3):211–24. doi: 10.1016/S2215-0366(18)30511-X
44. Schultz SK, Andreasen NC. Schizophrenia. Lancet (1999) 353(9162):1425–30. doi: 10.1016/S0140-6736(98)07549-7
45. Bassett AS, Collins EJ, Nuttall SE, Honer WG. Positive and negative symptoms in families with schizophrenia. Schizophr Res (1993) 11(1):9–19. doi: 10.1016/0920-9964(93)90033-F
46. Xia S, Feng J, Lei L, Hu J, Xia L, Wang J, et al. Comprehensive characterization of tissue-specific circular RNAs in the human and mouse genomes. Brief Bioinform (2017) 18(6):984–92.
47. Maass PG, Glažar P, Memczak S, Dittmar G, Hollfinger I, Schreyer L, et al. A map of human circular RNAs in clinically relevant tissues. J Mol Med (Berl) (2017) 95(11):1179–89. doi: 10.1007/s00109-017-1582-9
48. Liu S, Zhang F, Wang X, Shugart YY, Zhao Y, Li X, et al. Diagnostic value of blood-derived microRNAs for schizophrenia: results of a meta-analysis and validation. Sci Rep (2017) 7(1):15328. doi: 10.1038/s41598-017-15751-5
49. Wang Y, Wang J, Guo T, Peng Y, Wang K, Bai K, et al. Screening of schizophrenia associated miRNAs and the regulation of miR-320a-3p on integrin β1. Med (Baltimore) (2019) 98(8):e14332. doi: 10.1097/MD.0000000000014332
50. Du Y, Yu Y, Hu Y, Li X-W, Wei Z-X, Pan R-Y, et al. Genome-Wide, Integrative Analysis Implicates Exosome-Derived MicroRNA Dysregulation in Schizophrenia. Schizophr Bull (2019) 45(6):1257–66. doi: 10.1093/schbul/sby191
51. Wei H, Yuan Y, Liu S, Wang C, Yang F, Lu Z, et al. Detection of circulating miRNA levels in schizophrenia. Am J Psychiatry (2015) 172(11):1141–7. doi: 10.1176/appi.ajp.2015.14030273
52. Liu S, Zhang F, Shugart YY, Yang L, Li X, Liu Z, et al. The early growth response protein 1-miR-30a-5p-neurogenic differentiation factor 1 axis as a novel biomarker for schizophrenia diagnosis and treatment monitoring. Transl Psychiatry (2017) 7(1):e998. doi: 10.1038/tp.2016.268
53. Ebbesen KK, Kjems J, Hansen TB. Circular RNAs: Identification, biogenesis and function. Biochim Biophys Acta (2016) 1859(1):163–8. doi: 10.1016/j.bbagrm.2015.07.007
54. Meng X, Li X, Zhang P, Wang J, Zhou Y, Chen M. Circular RNA: an emerging key player in RNA world. Brief Bioinform (2017) 18(4):547–57.
55. Venø MT, Hansen TB, Venø ST, Clausen BH, Grebing M, Finsen B, et al. Spatio-temporal regulation of circular RNA expression during porcine embryonic brain development. Genome Biol (2015) 16:245. doi: 10.1186/s13059-015-0801-3
56. Rybak-Wolf A, Stottmeister C, Glažar P, Jens M, Pino N, Giusti S, et al. Circular RNAs in the Mammalian Brain Are Highly Abundant, Conserved, and Dynamically Expressed. Mol Cell (2015) 58(5):870–85. doi: 10.1016/j.molcel.2015.03.027
57. Hanan M, Soreq H, Kadener S. CircRNAs in the brain. RNA Biol (2017) 14(8):1028–34. doi: 10.1080/15476286.2016.1255398
58. Kumra S, Charles Schulz S. Editorial: research progress in early-onset schizophrenia. Schizophr Bull (2008) 34(1):15–7. doi: 10.1093/schbul/sbm123
59. Ng SY, Bogu GK, Soh BS, Stanton LW. The long noncoding RNA RMST interacts with SOX2 to regulate neurogenesis. Mol Cell (2013) 51(3):349–59. doi: 10.1016/j.molcel.2013.07.017
60. You X, Vlatkovic I, Babic A, Will T, Epstein I, Tushev G, et al. Neural circular RNAs are derived from synaptic genes and regulated by development and plasticity. Nat Neurosci (2015) 18(4):603–10. doi: 10.1038/nn.3975
61. Ahmad R, Sportelli V, Ziller M, Spengler D, Hoffmann A. Tracing Early Neurodevelopment in Schizophrenia with Induced Pluripotent Stem Cells. Cells (2018), 7(9. doi: 10.3390/cells7090140
62. Mahmoudi E, Fitzsimmons C, Geaghan MP, Shannon Weickert C, Atkins JR, Wang X, et al. Circular RNA biogenesis is decreased in postmortem cortical gray matter in schizophrenia and may alter the bioavailability of associated miRNA. Neuropsychopharmacology (2019) 44(6):1043–54. doi: 10.1038/s41386-019-0348-1
63. Glažar P, Papavasileiou P, Rajewsky N. circBase: a database for circular RNAs. RNA (2014) 20(11):1666–70. doi: 10.1261/rna.043687.113
64. Liu Z, Ran Y, Tao C, Li S, Chen J, Yang E. Detection of circular RNA expression and related quantitative trait loci in the human dorsolateral prefrontal cortex. Genome Biol (2019) 20(1):99. doi: 10.1186/s13059-019-1701-8
65. Zimmerman AJ, Hafez AK, Amoah SK, Rodriguez BA, Dell'Orco M, Lozano E. A psychiatric disease-related circular RNA controls synaptic gene expression and cognition. Mol Psychiatry (2020). doi: 10.1038/s41380-020-0653-4
66. Yang Q, Wu J, Zhao J, Xu T, Zhao Z, Song X, et al. Circular RNA expression profiles during the differentiation of mouse neural stem cells. BMC Syst Biol (2018) 12(Suppl 8):128. doi: 10.1186/s12918-018-0651-1
67. Zaghlool A, Ameur A, Wu C, Westholm JO, Niazi A, Manivannan M, et al. Expression profiling and in situ screening of circular RNAs in human tissues. Sci Rep (2018) 8(1):16953. doi: 10.1038/s41598-018-35001-6
68. Yao G, Niu W, Zhu X, et al. hsa_circRNA_104597: a novel potential diagnostic and therapeutic biomarker for schizophrenia. Biomark Med (2019) 13(5):331–40. doi: 10.2217/bmm-2018-0447
69. Werfel S, Nothjunge S, Schwarzmayr T, Strom TM, Meitinger T, Engelhardt S. Characterization of circular RNAs in human, mouse and rat hearts. J Mol Cell Cardiol (2016) 98:103–7. doi: 10.1016/j.yjmcc.2016.07.007
70. Zhang S, Zhu D, Li H, Li H, Feng C, Zhang W. Characterization of circRNA-Associated-ceRNA Networks in a Senescence-Accelerated Mouse Prone 8 Brain. Mol Ther (2017) 25(9):2053–61. doi: 10.1016/j.ymthe.2017.06.009
71. Floris G, Zhang L, Follesa P, Sun T. Regulatory Role of Circular RNAs and Neurological Disorders. Mol Neurobiol (2017) 54(7):5156–65. doi: 10.1007/s12035-016-0055-4
72. Smith K. Mental health: a world of depression. Nature (2014) 515(7526):181. doi: 10.1038/515180a
73. Harris EC, Barraclough B. Excess mortality of mental disorder. Br J Psychiatry (1998) 173:11–53. doi: 10.1192/bjp.173.1.11
74. Courtet P, Lopez-Castroman J. Antidepressants and suicide risk in depression. World Psychiatry (2017) 16(3):317–8. doi: 10.1002/wps.20460
75. Ponsoni A, Branco LD, Cotrena C, Shansis FM, Grassi-Oliveira R, Fonseca RP. Self-reported inhibition predicts history of suicide attempts in bipolar disorder and major depression. Compr Psychiatry (2018) 82:89–94. doi: 10.1016/j.comppsych.2018.01.011
76. Holden C. Mental health. Global survey examines impact of depression. Science (2000) 288(5463):39–40. doi: 10.1126/science.288.5463.39
77. Cui X, Niu W, Kong L, He M, Jiang K, Chen S, et al. hsa_circRNA_103636: potential novel diagnostic and therapeutic biomarker in Major depressive disorder. Biomark Med (2016) 10(9):943–52. doi: 10.2217/bmm-2016-0130
78. Danna SM, Graham E, Burns RJ, Deschênes SS, Schmitz N. Association between Depressive Symptoms and Cognitive Function in Persons with Diabetes Mellitus: A Systematic Review. PloS One (2016) 11(8):e0160809. doi: 10.1371/journal.pone.0160809
79. Katon W. Depression and diabetes: unhealthy bedfellows. Depress Anxiety. (2010) 27: (4):323–6. doi: 10.1002/da.20683
80. An T, Zhang J, Ma Y, Lian J, Wu YX, Lv BH, et al. Relationships of Non-coding RNA with diabetes and depression. Sci Rep (2019) 9(1):10707. doi: 10.1038/s41598-019-47077-9
81. Abasolo N, Torrell H, Roig B, Moyano S, Vilella E, Martorell L. RT-qPCR study on post-mortem brain samples from patients with major psychiatric disorders: reference genes and specimen characteristics. J Psychiatr Res (2011) 45(11):1411–8. doi: 10.1016/j.jpsychires.2011.06.001
82. Jiang G, Ma Y, An T, Pan Y, Mo F, Zhao D, et al. Relationships of circular RNA with diabetes and depression. Sci Rep (2017) 7(1):7285. doi: 10.1038/s41598-017-07931-0
83. Hunter AM, Leuchter AF, Power RA, Muthén B, McGrath PJ, Lewis CM, et al. A genome-wide association study of a sustained pattern of antidepressant response. J Psychiatr Res (2013) 47(9):1157–65. doi: 10.1016/j.jpsychires.2013.05.002
84. Zhang Y, Du L, Bai Y, Han B, He C, Gong L. CircDYM ameliorates depressive-like behavior by targeting miR-9 to regulate microglial activation via HSP90 ubiquitination. Mol Psychiatry (2018). doi: 10.1038/s41380-018-0285-0
85. Cohen JE, Lee PR, Fields RD. Systematic identification of 3′-UTR regulatory elements in activity-dependent mRNA stability in hippocampal neurons. Philos Trans R Soc Lond B Biol Sci (2014) 369(1652). doi: 10.1098/rstb.2013.0509
86. Long B, Wang K, Li N, Murtaza I, Xiao JY, Fan YY, et al. miR-761 regulates the mitochondrial network by targeting mitochondrial fission factor. Free Radic Biol Med (2013) 65:371–9. doi: 10.1016/j.freeradbiomed.2013.07.009
87. Xu Y, Zhao C, Sun X, Liu Z, Zhang J. MicroRNA-761 regulates mitochondrial biogenesis in mouse skeletal muscle in response to exercise. Biochem Biophys Res Commun (2015) 467(1):103–8. doi: 10.1016/j.bbrc.2015.09.113
Keywords: schizophrenia (SZ), depression (DEP), circular RNA (circRNA), biological function, expression, epigenetic characteristics
Citation: Li Z, Liu S, Li X, Zhao W, Li J and Xu Y (2020) Circular RNA in Schizophrenia and Depression. Front. Psychiatry 11:392. doi: 10.3389/fpsyt.2020.00392
Received: 05 November 2019; Accepted: 17 April 2020;
Published: 07 May 2020.
Edited by:
Maria S. Garcia-Gutierrez, Universidad Miguel Hernández de Elche, SpainReviewed by:
Bhaskar Roy, University of Alabama at Birmingham, United StatesCarlos M. Opazo, University of Melbourne, Australia
Copyright © 2020 Li, Liu, Li, Zhao, Li and Xu. 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: Yong Xu, xuyongsmu@vip.163.com
†These authors have contributed equally to this work