- Department of Pediatric Cardiology, Shandong Provincial Hospital, Cheeloo College of Medicine, Shandong University, Jinan, China
Circular RNAs (circRNAs) have emerged as essential regulators and biomarkers in various diseases. To assess the different expression levels of circRNAs in pediatric dilated cardiomyopathy (PDCM) and explore their biological and mechanistic significance, we used RNA microarrays to identify differentially expressed circRNAs between three children diagnosed with PDCM and three healthy age-matched volunteers. The biological function of circRNAs was assessed with a circRNA–microRNA (miRNA)–mRNA interaction network constructed from Gene Ontology and the Kyoto Encyclopedia of Genes and Genomes. Differentially expressed circRNAs were validated by quantitative real-time polymerase chain reaction (qRT-PCR) in 25 children with PDCM and 25 healthy volunteers. We identified 257 up-regulated (fold change ≤ 0.5, P < 0.05) and 899 down-regulated (fold change ≥2, P < 0.05) circRNAs in PDCM patients when compared to healthy volunteers. The qRT-PCR experiments confirmed has_circ_0067735 down-regulation (0.45-fold, P < 0.001), has_circ_0070186 up-regulation (2.82-fold, P < 0.001), and has_circ_0069972 down-regulation (0.50-fold, P < 0.05). A functional analysis of these differentially expressed circRNAs suggests that they are associated with hypertrophy, remodeling, fibrosis, and autoimmunity. CircRNAs have been implicated in PDCM through largely unknown mechanisms. Here we report differentially expressed circRNAs in PDCM patients that may illuminate the mechanistic roles in the etiology of PDCM that could serve as non-invasive diagnostic biomarkers.
Introduction
Pediatric dilated cardiomyopathy (PDCM) is characterized by left ventricular dilation and systolic dysfunction and commonly results in progressive congestive heart failure, arrhythmia, and sudden cardiac death (Towbin et al., 2006; Weintraub et al., 2017; You et al., 2019). PDCM often has a poor prognosis and is the leading cause of cardiac transplantation worldwide (Everly, 2008; Hertz et al., 2009). It is believed that PDCM pathogenesis is caused by a combination of genetic susceptibility and environmental insults (Japp et al., 2016). However, the specific pathogenic mechanisms of the disease remain unclear. PDCM in children appears to have a wider spectrum of causes than dilated cardiomyopathy (DCM) in adults (Griffin et al., 1988; Towbin, 1999), including idiopathic, genetic mutations, myocarditis, neuromuscular disorders, and inborn metabolic dysfunction. Recently, gene modification and epigenetic regulation have become major sources of investigation to identify the mechanism of pathogenesis in PDCM (Hershberger et al., 2013; Wu et al., 2015).
Circular RNAs (circRNAs) are a class of endogenous coding and non-coding RNA created from precursor mRNA back-splicing in eukaryotes (Chen, 2016). Structurally, circRNAs form unique covalent rings without 5′ caps or 3′ polyadenylated tails (Memczak et al., 2013) and thus exhibit greater stability and sequence conservation than normal linear RNA molecules (Guo et al., 2014). Moreover, circRNAs are generally cell and tissue specific (Wilusz and Sharp, 2013). CircRNAs are involved in numerous regulatory processes, including transcriptional modulation, splicing interference, miRNA sequestration, and translation (Fang, 2018; Li et al., 2018). In recent years, a growing number of studies demonstrate that some circRNAs may act as regulatory “miRNA sponges” that naturally sequester and competitively suppress miRNA activity, suggesting that circRNAs might play important roles in post-transcriptional regulation (Memczak et al., 2013). For example, Zheng et al. reported that circHIPK3 directly binds to miR-124 and inhibits miR-124 activity, modulating cell growth (Zheng et al., 2016). Additionally, ciRS-7 was reported by Hansen TB et al. to strongly suppress miR-7 activity by acting as a sponge, resulting in increased levels of miR-7 targets (Hansen et al., 2013). CircRNAs have also been reported to regulate gene expression by interacting with RNA binding proteins and translational regulators or by binding directly to mRNAs (Beltran-Garcia et al., 2020; Zang et al., 2020). Moreover, circRNAs can be translated in vitro and in vivo (Pamudurti et al., 2017).
Recent research on circRNAs has advanced our understanding of the mechanistic roles they play in cardiovascular diseases (Devaux et al., 2017; Aufiero et al., 2019; Zhang et al., 2019). For example, the circRNA Foxo3 was found to effectively reduce doxorubicin-induced cardiomyopathies, and it plays an important role in the senescence of mouse embryonic fibroblasts (Du et al., 2017). Additionally, Wang K recently verified that heart-related circRNA could protect the heart from pathological hypertrophy and heart failure by inhibiting miR-223 activity (Wang et al., 2016). However, in-depth investigations into the role circRNAs play in PDCM pathogenesis are needed to develop early diagnostic techniques and advance new therapeutic targets.
Here we assess the expression patterns and functions of circRNAs in PDCM patients and compare them to healthy control participants using circRNA microarray analysis. We constructed a circRNA–microRNA (miRNA)–mRNA interaction network with Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG). Finally, we verified four differentially expressed DCM-associated circRNAs in blood samples from 25 PDCM patients and 25 healthy volunteers by quantitative real-time polymerase chain reaction (qRT-PCR). Our findings provide a framework for the role circRNAs may play in PDCM pathogenesis. Furthermore, we identified three potential biomarkers of pediatric DCM in the peripheral blood that may serve future diagnostic significance.
Materials and Methods
Patients and Peripheral Blood Samples
Twenty-five peripheral blood samples were obtained from children with PDCM and 25 from healthy children between March 2019 and June 2020. All PDCM cases were clinically diagnosed in strict accordance with the World Health Organization guidelines (Richardson et al., 1996). The inclusion criteria included the following: (i) younger than 18 years, (ii) left ventricular end-diastolic diameter z-score >+2, after body surface area correction (Everitt et al., 2014), and (iii) left ventricular ejection fraction ≤ 45%. Patients with hypertension, congenital heart disease, ischemic heart disease, and malformations were excluded. The healthy volunteers were age- and sex-matched with the PDCM cases. The clinical characteristics of the 25 patients and the 25 volunteers are summarized in Table 1. Six samples were subjected to microarray analysis, and all 50 samples were used for subsequent RT-qPCR validation. This study was approved by the Institutional Ethics Committee (NSFC: NO. 2018-115), and the participants provided informed consent or assent (parental informed consent for minors).
Table 1. Clinical characteristics of the pediatric dilated cardiomyopathy (PDCM) and control groups.
RNA Extraction and Quality Control
Leukocytes were isolated from whole peripheral blood via centrifugation (1,500 g for 15 min at 4°C) after lysing the red blood cells with Red Blood Lysis Buffer (Solarbio, China). RNA from the three paired samples was extracted with RNeasy Mini Kit (Qiagen, Germany) per the manufacturer's instructions. RNA from the other 44 samples was extracted using Sparkzol Reagent (SparkJade Science Co., Ltd., China), chloroform, and isopropanol precipitation. The RNA Integrity Number (RIN) of the three paired samples was determined using an Agilent 2100 Bioanalyzer (Agilent Technologies, CA, USA); all samples had RIN ≥7.5. The integrity of all 50 RNA samples was determined by agarose gel electrophoresis; RNA concentration and purity were quantified with a Nanodrop 2000 spectrophotometer (Thermo Scientific, Waltham, MA). The RNA samples were stored at −80°C until further use.
CircRNA Microarray Analysis
CircRNA microarray analysis was performed at Shanghai Sinomics Corporation (Shanghai, China), using Sino human ceRNA array V3.0 (Sinomics Corporation, China). cRNA synthesis and labeling, chip hybridization, washing, and image scanning were performed per manufacturer's instructions. Microarray data were extracted and visualized using the Feature Extraction software 10.7(Agilent Technologies), and the resulting raw data were subjected to quantile normalization using the limma package in R. Data analysis was performed according to Agilent Technologies at Sinotech Genomics Corporation. Differentially expressed circRNAs were identified using a fold change cutoff of 2 and a P-value of 0.05.
Functional Analysis of Differentially Expressed circRNAs
To determine the functional roles of differentially expressed circRNAs in PDCM, we predicted their respective miRNA response elements (MREs) and built a circRNA–miRNA–mRNA network with Arraystar's proprietary miRNA target prediction software. GO and KEGG pathway enrichment analyses were performed using the R package clusterProfiler.
Quantitative Real-Time Polymerase Chain Reaction
qRT-PCR was performed to quantify the circRNA expression levels on a LightCycler480 system (Roche Diagnostics, Switzerland) using SYBR Green Pro Taq HS Premix (AG11701, Accurate Biotechnology, Hunan, China). Briefly, 1,000 ng of total RNA was reverse-transcribed into cDNA with random primers. All reactions were performed in triplicate containing 2× SYBR Green Premix, 20 ng of template cDNA, and 8,000 nM primers in a final volume of 20 μl. Reaction mixes were analyzed in a 96-well optical reaction plate. Melting curves were analyzed, and PCR products were validated by Sanger sequencing. All primers for qRT-PCR were synthesized by BioSune Biotechnology Co., Ltd. (China) and are listed in Table 2. In addition, circRNA expression levels were normalized to the housekeeping gene β-actin and determined by the 2−ΔΔCT method.
Statistical Analysis
All statistical analyses were conducted using SPSS software (version 24.0, SPSS, IL, USA) and GraphPad Prism version 8.0 (GraphPad Software, CA, USA). Normal distribution of data was presented as mean ± standard deviation (mean ± SD) and compared by Student's t-test (two-tailed, unpaired, equal variance). P < 0.05 was considered statistically significant.
Results
Expression Profiles of circRNAs in Pediatric Dilated Cardiomyopathy
We identified 53,635 circRNAs in three children with PDCM and three healthy volunteers with circRNA microarray analysis. We observed differential circRNA expression levels with various p-values and fold changes between the PDCM case and control groups (Figure 1A), in addition to differential mRNA expression (Figure 1B). Hierarchical clustering analysis indicated a distinct circRNA (Figure 1C) and mRNA (Figure 1D) expression profile in PDCM patients when compared to healthy volunteers.
Figure 1. Circulating expression levels of circRNAs and mRNA in the pediatric dilated cardiomyopathy group and control group. (A) A volcano plot of differential circRNA expression with various p-values and fold changes. (B) A volcano plot of differential mRNA expression with various p-values and fold changes. (C) Hierarchical clustering analysis of circRNAs with altered expression (p < 0.05, fold change >2) between the two groups. Red strip, high relative expression; green strip, low relative expression; black strip, no change. Color brightness reflects the degree of expression increase or decrease. (D) Hierarchical clustering analysis of mRNAs with altered expression (p < 0.05, fold change >2) between the two groups. Red strip, high-relative expression; green strip, low-relative expression; black strip, no-change. Color brightness reflects the degree of expression increase or decrease. (E) The dysregulated circRNA distribution rule in human chromosomes. (F) Length distribution of the dysregulated circRNAs.
A total of 1,156 differentially expressed circRNAs were detected in PDCM cases compared with control patients (fold change >2, p < 0.05), including 257 up-regulated and 899 down-regulated circRNAs (Supplementary File 1). The 10 most up-regulated and 10 most down-regulated circRNAs according to fold change are presented in Table 3. Additionally, 482 mRNAs were up-regulated and 159 were down-regulated in the PDCM group compared to the control group (Supplementary File 2). To explore the molecular characteristics of circRNAs, we assessed the length and the chromosomal distribution of the differentially expressed circRNAs. We found that the differentially expressed circRNAs in PDCM patients were derived mainly from chromosome 2 (10.03%; 116/1,156), chromosome 1 (8.91%; 103/1,156), and chromosome 5 (6.49%; 75/1,156) in descending order (Figure 1E). Additionally, the lengths of these circRNAs were often <2,000 nucleotides (81.06%; 937/1,156), as shown in Figure 1F.
Table 3. Serial numbers, p-values, fold changes, chromosomes of origin, sequence lengths, and host genes of the top 10 up-regulated and down-regulated circular RNAs associated with pediatric dilated cardiomyopathy.
Construction of a Functional circRNA–miRNA–mRNA Interaction Network
According to previous studies (Hansen et al., 2013), circRNAs can bind to the MREs on miRNAs in a competitive manner to terminate their regulatory effects on target genes. We assessed the likely interactions of differentially expressed PDCM circRNAs with complementary miRNAs using Cytoscape 3.5 and present the five most likely miRNA binding sites for each circRNA in Supplementary File 3. We then constructed our functional network of the top 10 up-regulated and top 10 down-regulated circRNAs with their respective miRNA targets (Figure 2A). We validated our circRNA–miRNA–mRNA interaction network of three differentially expressed circRNAs in 50 samples by qRT-PCR (Figure 2B). Interestingly, the differentially expressed mRNAs that we identified, including CACNA2D2, IGF1, PRKCA, PIK3CA, VAV3, PRKCQ, TLR4, IL1B, TLR8, and CTNNBIP1, are involved in pathways relevant to DCM pathogenesis such as “dilated cardiomyopathy,” “leukocyte transendothelial migration,” “T cell receptor signaling pathways,” “Toll-like receptor signaling pathways,” and “WNT signaling pathways.”
Figure 2. CircRNA–miRNA–mRNA interaction network. (A) CircRNA–miRNA interaction network of the top 10 up-regulated and down-regulated circRNAs. (B) CircRNA–miRNA–mRNA interaction network of three selected circRNAs. Inverted triangles represent circRNAs; rectangles represent the predicted miRNAs; ellipse nodes represent mRNAs (red, up-regulated; green, down-regulated).
Validation of circRNA Expression
We designed primer pairs for the top 10 differentially up-regulated and down-regulated circRNAs in PDCM. Three pairs successfully amplified the target circRNA sequences, which were validated by RT-qPCR. Two circRNAs of interest, has_circ_0067735 (Figure 3A) and has_circ_0069972 (Figure 3B), were down-regulated 0.45-fold (p < 0.001) and 0.50-fold (p < 0.05) in the PDCM group, respectively, when compared to healthy participants. Another circRNA, named has_circ_0070186 (Figure 3C), was up-regulated 2.82-fold (p < 0.001) in PDCM patient samples. Meanwhile, has_circ_0005495 (Figure 3D) was confirmed to be down-regulated by 0.69-fold (p < 0.05), which was opposite to the microarray data.
Figure 3. Circulating expression levels of the four circRNAs assessed by qRT-PCR. (A) Circulating expression levels of hsa_circ_0067735 in control subjects and the children with pediatric dilated cardiomyopathy (PDCM). (B) Circulating expression levels of hsa_circ_0069972 in control subjects and the children with PDCM. (C) Circulating expression levels of hsa_circ_0070186 in control subjects and the children with PDCM. (D) Circulating expression levels of hsa_circ_0005495 in control subjects and the children with PDCM. Control group, n = 25; PDCM group, n = 25. *P < 0.05 vs. control; ***p < 0.001 vs. control.
Functional Prediction of Differentially Expressed circRNAs
To predict the underlying mechanisms of differentially expressed circRNAs in PDCM, we performed a functional annotation analysis of their target mRNAs with GO and the KEGG pathway enrichment tools (Supplementary Files 4, 5). These mRNA targets are likely modulated by circRNAs via competitive endogenous RNA (ceRNA) regulation. There are three GO categories, including biological process (BP), cellular component (CC), and molecular function (MF). The top 10 enriched GO terms (in BP, CC, and MF) and the top 30 KEGG pathways of the dysregulated circRNAs in PDCM are displayed in Figure 4. The three most enriched GO terms were “response to wounding,” “inflammatory response,” and “cytokine secretion” in BP; the top three terms in CC were “endomembrane system,” “ruffle membrane,” and “endosome membrane.” In MF, “enzyme binding,” “SMAD binding,” and “peroxidase activity” were the three most enriched GO terms. The KEGG pathways in differentially expressed circRNAs were enriched in “Toll-like receptor signaling” and “Fc gamma R-mediated phagocytosis,” processes commonly associated with PDCM.
Figure 4. Gene Ontology and the Kyoto Encyclopedia of Genes and Genomes pathway analyses of dysregulated circRNAs in pediatric dilated cardiomyopathy (PDCM). (A) Ten most enriched gene ontology terms in biological process, cellular component, and molecular function for differentially expressed circRNAs. (B) Top 30 Kyoto Encyclopedia of Gene and Genome pathways of dysregulated circRNAs. EGFR, epidermal growth factor receptor.
Discussion
It is widely believed that PDCM etiology is driven by both genetic and environmental factors (Poller et al., 2005), with the specific molecular mechanisms largely unknown. Circular RNAs (circRNAs) were recently discovered to be widely expressed across species and have been implicated in several diseases (Pan et al., 2018). However, the expression profiles and functions of circRNAs in PDCM remain elusive. Here we investigated the different expression patterns of circRNAs in peripheral blood leukocytes isolated from PDCM patients when compared to healthy age-matched individuals. A total of 1,156 circRNAs had differential expression in PDCM (fold change >2, p < 0.05), including 257 up-regulated and 899 down-regulated transcripts. This is the first study to assess the expression patterns of circRNAs in PDCM and may provide new insights into their mechanistic role in this disease. Four differentially expressed circRNAs from three patients were verified by qRT-PCR in 25 patient samples. Two circRNAs, has_circ_0067735 and has_circ_0069972, were significantly downregulated in DCM, while another, has_circ_0070186, was significantly upregulated. These three circRNAs may serve as potential biomarkers of PDCM in the future. We discovered that these three circRNAs possessed miRNA binding sites, suggesting that they carried a potential for post-transcriptional regulation. has_circ_0067735 likely regulates the expression of CACNA2D2 by binding hsa-miR-4262, while has_circ_0070186 likely regulates IGF1 expression by binding hsa-miR-4448. CACNA2D2 and IGF1 are associated with dilated cardiomyopathy pathway according to KEGG pathway enrichment analyses. However, these associations require further investigation.
The GO biological process and the KEGG pathway enrichment analyses were carried out to explore the potential mechanisms of dysregulated circRNAs in PDCM. PDCM-associated circRNAs were enriched in BP, such as “response to wounding,” “inflammatory response,” and “cytokine secretion.” Meanwhile, the KEGG analysis identified that DCM-associated circRNAs were enriched in pathways such as Toll-like receptor signaling, leukocyte transendothelial migration, T cell receptor signaling, and WNT signaling. These pathways have all previously been implicated in DCM pathogenesis. For example, TLR4 activation causes experimental autoimmune myocarditis progress to DCM in mice (Wu et al., 2018), which was closely related to PDCM. Additionally, WNT signaling is a critical pathway for cardiac hypertrophy and remodeling (Bergmann, 2010; Malekar et al., 2010; Lu et al., 2016). CD4+ T-cells may play a critical role in ADP/ATP carrier-caused mouse DCM (Wang et al., 2006). Inflammatory endothelial activation and migration of immunocompetent cells have been observed in 67% of DCM patients, a process mediated by cell adhesion (Noutsias et al., 1999, 2003). PDCM pathogenesis is usually associated with cardiac hypertrophy, remodeling, fibrosis, and autoimmunity. Therefore, given that the relevant circRNAs we identify relate to these cellular processes, it is likely that their differential expression contributes to PDCM pathogenesis via these pathways.
While our findings provide a framework for understanding the role circRNAs play in PDCM etiology, the differential expression of the circRNAs identified here should be verified first in myocardial tissue. Second, these findings should be verified across broader sociodemographic characteristics, including ethnicity. Our participants lacked sociodemographic breadth and were mostly of Asian ethnicity. The expression patterns of the circRNAs we identified here should also be analyzed in other cardiac pathologies to ensure that they are specific to PDCM. Finally, follow-up experiments are needed at the cellular and organismal level to determine the mechanistic functions of these circRNAs in PDCM.
Conclusion
This study provides the first profile of differentially expressed circRNAs in PDCM. We used GO and KEGG pathway enrichment analyses to construct a circRNA–miRNA–mRNA interaction network to preliminarily assess the roles and potential mechanisms of dysregulated circRNAs in the development of PDCM. We identified three relevant circRNAs in pediatric DCM: has_circ_0067735 and has_circ_0069972 were markedly down-regulated, while has_circ_0070186 was up-regulated, suggesting that they may constitute candidate biomarkers of PDCM. This work provides a foundation for further research on the mechanistic role circRNAs play in PDCM.
No effective circRNA has been utilized for early PDCM diagnosis. However, the circRNAs that we identified have potential as novel non-invasive biomarkers in PDCM.
Future Perspective
Our study provides a new foundation for understanding the roles circRNAs play in PDCM. We hope to explore their exact mechanism of action in a follow-up series of experiments. We believe that has_circ_0067735, has_circ_0070186, and has_circ_0069972 could serve as novel biomarkers and therapeutic targets in PDCM. Ultimately, understanding the role of differently expressed circRNAs has the potential to improve detection, prevention, and treatment of PDCM early when the maximum clinical benefit can be achieved.
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics Statement
The studies involving human participants were reviewed and approved by the Institutional Ethics Committee, Provincial Hospital affiliated to Shandong First Medical University. Written informed consent to participate in this study was provided by the participants' legal guardian/next of kin.
Author Contributions
BH and WS conceived and supervised the study. WS and JW designed the experiments. WS and DC performed the experiments. DJ and HJ analyzed and interpreted the data. WS reviewed and edited the manuscript for important intellectual content. BH and JW provided substantive revisions to the manuscript. All the authors approved the final version of the manuscript.
Funding
This study was supported by the National Natural Science Foundation of China (grant number 8187020860), Jinan Science and Technology Development Plan (grant number 201805020), and Special Expert of Taishan Scholars (grant number ts201511099).
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.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmolb.2020.600170/full#supplementary-material
Supplementary File 1. Differentially expressed circRNAs in pediatric dilated cardiomyopathy.
Supplementary File 2. Differentially expressed mRNAs in pediatric dilated cardiomyopathy.
Supplementary File 3. miRNA binding sites for circRNAs.
Supplementary File 4. Gene Ontology analysis.
Supplementary File 5. Kyoto Encyclopedia of Genes and Genomes pathway.
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
Beltran-Garcia, J., Osca-Verdegal, R., Nacher-Sendra, E., Pallardo, F. V., and Garcia-Gimenez, J. L. (2020). Circular RNAs in sepsis: biogenesis, function, and clinical significance. Cells 9:1544. doi: 10.3390/cells9061544
Bergmann, M. W. (2010). WNT signaling in adult cardiac hypertrophy and remodeling: lessons learned from cardiac development. Circ. Res. 107, 1198–1208. doi: 10.1161/CIRCRESAHA.110.223768
Chen, L. (2016). The biogenesis and emerging roles of circular RNAs. Nat. Rev. Mol. Cell Biol. 17, 205–211. doi: 10.1038/nrm.2015.32
Devaux, Y., Creemers, E. E., Boon, R. A., Werfel, S., Thum, T., Engelhardt, S., et al. (2017). Circular RNAs in heart failure. Eur. J. Heart Fail. 19, 701–709. doi: 10.1002/ejhf.801
Du, W. W., Yang, W., Chen, Y., Wu, Z. K., Foster, F. S., Yang, Z., et al. (2017). Foxo3 circular RNA promotes cardiac senescence by modulating multiple factors associated with stress and senescence responses. Eur. Heart J. 38, 1402–1412. doi: 10.1093/eurheartj/ehw001
Everitt, M. D., Sleeper, L. A., Lu, M., Canter, C. E., Pahl, E., Wilkinson, J. D., et al. (2014). Recovery of echocardiographic function in children with idiopathic dilated cardiomyopathy. J. Am. Coll. Cardiol. 63, 1405–1413. doi: 10.1016/j.jacc.2013.11.059
Everly, M. J. (2008). Cardiac transplantation in the United States: an analysis of the UNOS registry. Clin. Transpl. 2008, 35–43.
Fang, Y. (2018). Circular RNAs as novel biomarkers with regulatory potency in human diseases. Fut. Sci. OA 4:O314. doi: 10.4155/fsoa-2018-0036
Griffin, M. L., Hernandez, A., Martin, T. C., Goldring, D., Bolman, R. M., Spray, T. L., et al. (1988). Dilated cardiomyopathy in infants and children. J. Am. Coll. Cardiol. 11, 139–144. doi: 10.1016/0735-1097(88)90179-9
Guo, J. U., Agarwal, V., Guo, H., and Bartel, D. P. (2014). Expanded identification and characterization of mammalian circular RNAs. Genome Biol. 15:409. doi: 10.1186/s13059-014-0409-z
Hansen, T. B., Jensen, T. I., Clausen, B. H., Bramsen, J. B., Finsen, B., Damgaard, C. K., et al. (2013). Natural RNA circles function as efficient microRNA sponges. Nature 495, 384–388. doi: 10.1038/nature11993
Hershberger, R. E., Hedges, D. J., and Morales, A. (2013). Dilated cardiomyopathy: the complexity of a diverse genetic architecture. Nat. Rev. Cardiol. 10, 531–547. doi: 10.1038/nrcardio.2013.105
Hertz, M. I., Aurora, P., Christie, J. D., Dobbels, F., Edwards, L. B., Kirk, R., et al. (2009). Scientific registry of the international society for heart and lung transplantation: introduction to the 2009 annual reports. J. Heart Lung Transpl. 28, 989–992. doi: 10.1016/j.healun.2009.08.005
Japp, A. G., Gulati, A., Cook, S. A., Cowie, M. R., and Prasad, S. K. (2016). the diagnosis and evaluation of dilated cardiomyopathy. J. Am. Coll. Cardiol. 67, 2996–3010. doi: 10.1016/j.jacc.2016.03.590
Li, X., Yang, L., and Chen, L. (2018). The biogenesis, functions, and challenges of circular RNAs. Mol. Cell 71, 428–442. doi: 10.1016/j.molcel.2018.06.034
Lu, D., Bao, D., Dong, W., Liu, N., Zhang, X., Gao, S., et al. (2016). Dkk3 prevents familial dilated cardiomyopathy development through wnt pathway. Lab. Invest. 96, 239–248. doi: 10.1038/labinvest.2015.145
Malekar, P., Hagenmueller, M., Anyanwu, A., Buss, S., Streit, M. R., Weiss, C. S., et al. (2010). Wnt signaling is critical for maladaptive cardiac hypertrophy and accelerates myocardial remodeling. Hypertension 55, 939–945. doi: 10.1161/HYPERTENSIONAHA.109.141127
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
Noutsias, M., Pauschinger, M., Schultheiss, H. P., and Kuhl, U. (2003). Cytotoxic perforin+ and TIA-1+ infiltrates are associated with cell adhesion molecule expression in dilated cardiomyopathy. Eur. J. Heart Fail. 5, 469–479. doi: 10.1016/S1388-9842(03)00037-0
Noutsias, M., Seeberg, B., Schultheiss, H. P., and Kuhl, U. (1999). Expression of cell adhesion molecules in dilated cardiomyopathy: evidence for endothelial activation in inflammatory cardiomyopathy. Circulation 99, 2124–2131. doi: 10.1161/01.CIR.99.16.2124
Pamudurti, N. R., Bartok, O., Jens, M., Ashwal-Fluss, R., Stottmeister, C., Ruhe, L., et al. (2017). Translation of CircRNAs. Mol. Cell 66, 9–21. doi: 10.1016/j.molcel.2017.02.021
Pan, X., Xiong, K., Anthon, C., Hyttel, P., Freude, K. K., Jensen, L. J., et al. (2018). WebCircRNA: classifying the circular RNA potential of coding and noncoding RNA. Genes 9:536. doi: 10.3390/genes9110536
Poller, W., Kuhl, U., Tschoepe, C., Pauschinger, M., Fechner, H., and Schultheiss, H. P. (2005). Genome-environment interactions in the molecular pathogenesis of dilated cardiomyopathy. J. Mol. Med. 83, 579–586. doi: 10.1007/s00109-005-0664-2
Richardson, P., McKenna, W., Bristow, M., Maisch, B., Mautner, B., O'Connell, J., et al. (1996). Report of the 1995 world health organization/international society and federation of cardiology task force on the definition and classification of cardiomyopathies. Circulation 93, 841–842. doi: 10.1161/01.CIR.93.5.841
Towbin, J. A. (1999). Pediatric myocardial disease. Pediatr. Clin. North Am. 46, 289–312. doi: 10.1016/S0031-3955(05)70119-X
Towbin, J. A., Lowe, A. M., Colan, S. D., Sleeper, L. A., Orav, E. J., Clunie, S., et al. (2006). Incidence, causes, and outcomes of dilated cardiomyopathy in children. JAMA 296, 1867–1876. doi: 10.1001/jama.296.15.1867
Wang, K., Long, B., Liu, F., Wang, J. X., Liu, C. Y., Zhao, B., et al. (2016). A circular RNA protects the heart from pathological hypertrophy and heart failure by targeting miR-223. Eur. Heart J. 37, 2602–2611. doi: 10.1093/eurheartj/ehv713
Wang, Z., Liao, Y. H., Yuan, J., Zhang, J. H., Liu, Z. P., and Dong, J. H. (2006). Analysis of IgG subclass antibodies and expression of T-cell receptor signaling molecules in anti-CD4 monoclonal antibody treated mice with autoimmune cardiomyopathy. Autoimmunity 39, 455–460. doi: 10.1080/08916930600845915
Weintraub, R. G., Semsarian, C., and Macdonald, P. (2017). Dilated cardiomyopathy. Lancet 390, 400–414. doi: 10.1016/S0140-6736(16)31713-5
Wilusz, J. E., and Sharp, P. A. (2013). A circuitous route to noncoding RNA. Science 340, 440–441. doi: 10.1126/science.1238522
Wu, B., Li, J., Ni, H., Zhuang, X., Qi, Z., Chen, Q., et al. (2018). TLR4 activation promotes the progression of experimental autoimmune myocarditis to dilated cardiomyopathy by inducing mitochondrial dynamic imbalance. Oxid. Med. Cell. Longev. 2018, 1–15. doi: 10.1155/2018/3181278
Wu, H., Lee, J., Vincent, L. G., Wang, Q., Gu, M., Lan, F., et al. (2015). Epigenetic regulation of phosphodiesterases 2A and 3A underlies compromised beta-adrenergic signaling in an iPSC model of dilated cardiomyopathy. Cell Stem Cell 17, 89–100. doi: 10.1016/j.stem.2015.04.020
You, H., Jiang, W., Jiao, M., Wang, X., Jia, L., You, S., et al. (2019). Association of soluble ST2 serum levels with outcomes in pediatric dilated cardiomyopathy. Can. J. Cardiol. 35, 727–735. doi: 10.1016/j.cjca.2019.02.016
Zang, J., Lu, D., and Xu, A. (2020). The interaction of circRNAs and RNA binding proteins: an important part of circRNA maintenance and function. J. Neurosci. Res. 98, 87–97. doi: 10.1002/jnr.24356
Zhang, L., Han, B., Wang, J., Liu, Q., Kong, Y., Jiang, D., et al. (2019). Differential expression profiles and functional analysis of circular RNAs in children with fulminant myocarditis. Epigenomics 11, 1129–1141. doi: 10.2217/epi-2019-0101
Keywords: biomarkers, microarray, pediatric dilated cardiomyopathy, circular RNAs (circRNAs), gene expression profile (GEP)
Citation: Sun W, Han B, Cai D, Wang J, Jiang D and Jia H (2020) Differential Expression Profiles and Functional Prediction of Circular RNAs in Pediatric Dilated Cardiomyopathy. Front. Mol. Biosci. 7:600170. doi: 10.3389/fmolb.2020.600170
Received: 03 September 2020; Accepted: 03 November 2020;
Published: 18 December 2020.
Edited by:
Jose Luis García-Giménez, Center for Biomedical Research in Rare Diseases Network (CIBERER), SpainReviewed by:
Jesus Beltran, University of Valencia, SpainNadiah Abu, National University of Malaysia, Malaysia
Copyright © 2020 Sun, Han, Cai, Wang, Jiang and Jia. 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: Bo Han, aGFuYm8zNSYjeDAwMDQwOzE2My5jb20=