Abstract
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
| PDCM group (n = 25) | Control group (n = 25) | P-value | |
|---|---|---|---|
| Sex | |||
| Male | 10 | 11 | |
| Female | 15 | 14 | |
| Age (years) (median, P25, P75) | 1.42 (0.71, 7.42) | 1.5 (0.63, 7.67) | |
| Echo | |||
| LVEDD Z-score (mean ± SD) | 6.74 ± 2.69 | −0.50 ± 0.84 | *** |
| LVEF (%) (mean ± SD) | 29.4 ± 7.88 | 63.92 ± 0.64 | *** |
| NT-ProBNP (pg/ml) (mean ± SD) | 12,003 ± 11,373 | 50.7 ± 37.02 | *** |
Clinical characteristics of the pediatric dilated cardiomyopathy (PDCM) and control groups.
Echo, echocardiographic; LVEDD, left ventricular end-diastolic diameter; Z-score, after body surface area correction, the distance from the average (normal range, 0 ± 2); LVEF, left ventricular ejection fraction (normal value >60%); NT-ProBNP, brain natriuretic peptide, an index of heart failure (normal range, 0–450 pg/ml).
p < 0.001 vs. control.
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.
Table 2
| Primer name | Forward primer sequence (5′-3′) | Reverse primer sequence (5′-3′) |
|---|---|---|
| homoβ-actin | AGTTGCGTTACACCCTTTCTTG | CACCTTCACCGTTCCAGTTTT |
| hsa_circ_0067735 | CAACGTGGCTGATCAAGATG | GCTGGCTGATTATTGAAGCCT |
| hsa_circ_0070186 | GGGCATTTTGAAGACTTACTG | ATGGATTTCTTTAGCTGCTTT |
| hsa_circ_0069972 | TCTAAACTTCATGGAAACAAGGAAA | TCTCTGCTGAAGACTGGGAA |
| hsa_circ_0005495 | TGTGGCCAGAAATTTGCCAG | AATGGGGTCCGGAACAAAGC |
Specific primers for qRT-PCR.
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
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
| CircRNA | p-values | FDR | Fold change | Regulation | Chromosome | Sequence length | Host gene |
|---|---|---|---|---|---|---|---|
| hsa_circ_0033063 | 0.01686 | 0.44991 | 38.38073 | Up | chr14 | 587 | IFI27 |
| hsa_circ_0088174 | 0.01637 | 0.45149 | 15.19901 | Up | chr9 | 471 | ORM1 |
| hsa_circ_0088177 | 0.01418 | 0.44940 | 12.08522 | Up | chr9 | 292 | ORM2 |
| hsa_circ_0070185 | 0.02702 | 0.43012 | 11.05281 | Up | chr4 | 531 | ANXA3 |
| hsa_circ_0070186 | 0.01789 | 0.44561 | 9.46793 | Up | chr4 | 591 | ANXA3 |
| hsa_circ_0070187 | 0.01300 | 0.45041 | 8.94111 | Up | chr4 | 396 | ANXA3 |
| hsa_circ_0070170 | 0.01325 | 0.44894 | 6.76701 | Up | chr4 | 1,364 | FRAS1 |
| hsa_circ_0025518 | 0.01989 | 0.43710 | 4.148068 | Up | chr12 | 927 | PLBD1 |
| hsa_circ_0074239 | 0.02866 | 0.43501 | 4.10679 | Up | chr5 | 692 | C5orf32 |
| hsa_circ_0078682 | 0.02615 | 0.43219 | 4.08332 | Up | chr6 | 641 | MLLT4 |
| hsa_circ_0058021 | 0.00107 | 0.46791 | 0.06619 | Down | chr2 | 348 | CPS1 |
| hsa_circ_0024749 | 0.00065 | 0.53912 | 0.10711 | Down | chr11 | 433 | FEZ1 |
| hsa_circ_0069972 | 0.04774 | 0.44405 | 0.11357 | Down | chr4 | 2,014 | CXCL5 |
| hsa_circ_0052078 | 0.03754 | 0.43460 | 0.12151 | Down | chr19 | 2,045 | PPP2R1A |
| hsa_circ_0027275 | 0.00724 | 0.45933 | 0.12608 | Down | chr12 | 1,520 | MBD6 |
| hsa_circ_0046702 | 0.00409 | 0.44719 | 0.14144 | Down | chr18 | 279 | YES1 |
| hsa_circ_0058810 | 0.03100 | 0.43439 | 0.14659 | Down | chr2 | 169 | AGAP1 |
| hsa_circ_0067735 | 0.04257 | 0.43888 | 0.15006 | Down | chr3 | 457 | MED12L |
| hsa_circ_0002515 | 0.02177 | 0.43474 | 0.18675 | Down | chr4 | 336 | AFAP1 |
| hsa_circ_0030584 | 0.00476 | 0.44922 | 0.19012 | Down | chr13 | 2,004 | ABCC4 |
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
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
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
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.
Statements
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 1Differentially expressed circRNAs in pediatric dilated cardiomyopathy.
Supplementary File 2Differentially expressed mRNAs in pediatric dilated cardiomyopathy.
Supplementary File 3miRNA binding sites for circRNAs.
Supplementary File 4Gene Ontology analysis.
Supplementary File 5Kyoto Encyclopedia of Genes and Genomes pathway.
References
1
AufieroS.ReckmanY. J.PintoY. M.CreemersE. E. (2019). Circular RNAs open a new chapter in cardiovascular biology. Nat. Rev. Cardiol.16, 503–514. 10.1038/s41569-019-0185-2
2
Beltran-GarciaJ.Osca-VerdegalR.Nacher-SendraE.PallardoF. V.Garcia-GimenezJ. L. (2020). Circular RNAs in sepsis: biogenesis, function, and clinical significance. Cells9:1544. 10.3390/cells9061544
3
BergmannM. W. (2010). WNT signaling in adult cardiac hypertrophy and remodeling: lessons learned from cardiac development. Circ. Res.107, 1198–1208. 10.1161/CIRCRESAHA.110.223768
4
ChenL. (2016). The biogenesis and emerging roles of circular RNAs. Nat. Rev. Mol. Cell Biol.17, 205–211. 10.1038/nrm.2015.32
5
DevauxY.CreemersE. E.BoonR. A.WerfelS.ThumT.EngelhardtS.et al. (2017). Circular RNAs in heart failure. Eur. J. Heart Fail.19, 701–709. 10.1002/ejhf.801
6
DuW. W.YangW.ChenY.WuZ. K.FosterF. S.YangZ.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. 10.1093/eurheartj/ehw001
7
EverittM. D.SleeperL. A.LuM.CanterC. E.PahlE.WilkinsonJ. D.et al. (2014). Recovery of echocardiographic function in children with idiopathic dilated cardiomyopathy. J. Am. Coll. Cardiol.63, 1405–1413. 10.1016/j.jacc.2013.11.059
8
EverlyM. J. (2008). Cardiac transplantation in the United States: an analysis of the UNOS registry. Clin. Transpl. 2008, 35–43.
9
FangY. (2018). Circular RNAs as novel biomarkers with regulatory potency in human diseases. Fut. Sci. OA4:O314. 10.4155/fsoa-2018-0036
10
GriffinM. L.HernandezA.MartinT. C.GoldringD.BolmanR. M.SprayT. L.et al. (1988). Dilated cardiomyopathy in infants and children. J. Am. Coll. Cardiol.11, 139–144. 10.1016/0735-1097(88)90179-9
11
GuoJ. U.AgarwalV.GuoH.BartelD. P. (2014). Expanded identification and characterization of mammalian circular RNAs. Genome Biol.15:409. 10.1186/s13059-014-0409-z
12
HansenT. B.JensenT. I.ClausenB. H.BramsenJ. B.FinsenB.DamgaardC. K.et al. (2013). Natural RNA circles function as efficient microRNA sponges. Nature495, 384–388. 10.1038/nature11993
13
HershbergerR. E.HedgesD. J.MoralesA. (2013). Dilated cardiomyopathy: the complexity of a diverse genetic architecture. Nat. Rev. Cardiol.10, 531–547. 10.1038/nrcardio.2013.105
14
HertzM. I.AuroraP.ChristieJ. D.DobbelsF.EdwardsL. B.KirkR.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. 10.1016/j.healun.2009.08.005
15
JappA. G.GulatiA.CookS. A.CowieM. R.PrasadS. K. (2016). the diagnosis and evaluation of dilated cardiomyopathy. J. Am. Coll. Cardiol.67, 2996–3010. 10.1016/j.jacc.2016.03.590
16
LiX.YangL.ChenL. (2018). The biogenesis, functions, and challenges of circular RNAs. Mol. Cell71, 428–442. 10.1016/j.molcel.2018.06.034
17
LuD.BaoD.DongW.LiuN.ZhangX.GaoS.et al. (2016). Dkk3 prevents familial dilated cardiomyopathy development through wnt pathway. Lab. Invest.96, 239–248. 10.1038/labinvest.2015.145
18
MalekarP.HagenmuellerM.AnyanwuA.BussS.StreitM. R.WeissC. S.et al. (2010). Wnt signaling is critical for maladaptive cardiac hypertrophy and accelerates myocardial remodeling. Hypertension55, 939–945. 10.1161/HYPERTENSIONAHA.109.141127
19
MemczakS.JensM.ElefsiniotiA.TortiF.KruegerJ.RybakA.et al. (2013). Circular RNAs are a large class of animal RNAs with regulatory potency. Nature495, 333–338. 10.1038/nature11928
20
NoutsiasM.PauschingerM.SchultheissH. P.KuhlU. (2003). Cytotoxic perforin+ and TIA-1+ infiltrates are associated with cell adhesion molecule expression in dilated cardiomyopathy. Eur. J. Heart Fail.5, 469–479. 10.1016/S1388-9842(03)00037-0
21
NoutsiasM.SeebergB.SchultheissH. P.KuhlU. (1999). Expression of cell adhesion molecules in dilated cardiomyopathy: evidence for endothelial activation in inflammatory cardiomyopathy. Circulation99, 2124–2131. 10.1161/01.CIR.99.16.2124
22
PamudurtiN. R.BartokO.JensM.Ashwal-FlussR.StottmeisterC.RuheL.et al. (2017). Translation of CircRNAs. Mol. Cell66, 9–21. 10.1016/j.molcel.2017.02.021
23
PanX.XiongK.AnthonC.HyttelP.FreudeK. K.JensenL. J.et al. (2018). WebCircRNA: classifying the circular RNA potential of coding and noncoding RNA. Genes9:536. 10.3390/genes9110536
24
PollerW.KuhlU.TschoepeC.PauschingerM.FechnerH.SchultheissH. P. (2005). Genome-environment interactions in the molecular pathogenesis of dilated cardiomyopathy. J. Mol. Med.83, 579–586. 10.1007/s00109-005-0664-2
25
RichardsonP.McKennaW.BristowM.MaischB.MautnerB.O'ConnellJ.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. Circulation93, 841–842. 10.1161/01.CIR.93.5.841
26
TowbinJ. A. (1999). Pediatric myocardial disease. Pediatr. Clin. North Am.46, 289–312. 10.1016/S0031-3955(05)70119-X
27
TowbinJ. A.LoweA. M.ColanS. D.SleeperL. A.OravE. J.ClunieS.et al. (2006). Incidence, causes, and outcomes of dilated cardiomyopathy in children. JAMA296, 1867–1876. 10.1001/jama.296.15.1867
28
WangK.LongB.LiuF.WangJ. X.LiuC. Y.ZhaoB.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. 10.1093/eurheartj/ehv713
29
WangZ.LiaoY. H.YuanJ.ZhangJ. H.LiuZ. P.DongJ. 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. Autoimmunity39, 455–460. 10.1080/08916930600845915
30
WeintraubR. G.SemsarianC.MacdonaldP. (2017). Dilated cardiomyopathy. Lancet390, 400–414. 10.1016/S0140-6736(16)31713-5
31
WiluszJ. E.SharpP. A. (2013). A circuitous route to noncoding RNA. Science340, 440–441. 10.1126/science.1238522
32
WuB.LiJ.NiH.ZhuangX.QiZ.ChenQ.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. 10.1155/2018/3181278
33
WuH.LeeJ.VincentL. G.WangQ.GuM.LanF.et al. (2015). Epigenetic regulation of phosphodiesterases 2A and 3A underlies compromised beta-adrenergic signaling in an iPSC model of dilated cardiomyopathy. Cell Stem Cell17, 89–100. 10.1016/j.stem.2015.04.020
34
YouH.JiangW.JiaoM.WangX.JiaL.YouS.et al. (2019). Association of soluble ST2 serum levels with outcomes in pediatric dilated cardiomyopathy. Can. J. Cardiol.35, 727–735. 10.1016/j.cjca.2019.02.016
35
ZangJ.LuD.XuA. (2020). The interaction of circRNAs and RNA binding proteins: an important part of circRNA maintenance and function. J. Neurosci. Res.98, 87–97. 10.1002/jnr.24356
36
ZhangL.HanB.WangJ.LiuQ.KongY.JiangD.et al. (2019). Differential expression profiles and functional analysis of circular RNAs in children with fulminant myocarditis. Epigenomics11, 1129–1141. 10.2217/epi-2019-0101
37
ZhengQ.BaoC.GuoW.LiS.ChenJ.ChenB.et al. (2016). Circular RNA profiling reveals an abundant circHIPK3 that regulates cell growth by sponging multiple miRNAs. Nat. Commun.7:11215. 10.1038/ncomms11215
Summary
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
Volume
7 - 2020
Edited by
Jose Luis García-Giménez, Center for Biomedical Research in Rare Diseases Network (CIBERER), Spain
Reviewed by
Jesus Beltran, University of Valencia, Spain; Nadiah Abu, National University of Malaysia, Malaysia
Updates
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 hanbo35@163.com
This article was submitted to Molecular Diagnostics and Therapeutics, a section of the journal Frontiers in Molecular Biosciences
Disclaimer
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.