PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of oncotargetLink to Publisher's site
 
Oncotarget. 2017 September 1; 8(36): 60280–60290.
Published online 2017 August 4. doi:  10.18632/oncotarget.19941
PMCID: PMC5601138

Circular RNAs promote TRPM3 expression by inhibiting hsa-miR-130a-3p in coronary artery disease patients

Abstract

We investigated the differential expression of circular RNAs (circRNAs) in plasma samples from three coronary artery disease (CAD) patients to identify putative therapeutic targets. We identified 24 differentially expressed circRNAs (18 up-regulated and 6 down-regulated) and 7 differentially expressed mRNAs (6 up-regulated and 1 down-regulated) in CAD patients based on competing endogenous RNA (ceRNA) microarray analysis. MiR-221(p = 0.001), miR-155(p = 0.049), and miR-130a (p = 0.001) were downregulated in CAD patients based on qRT-PCR analysis of another independent population of 932 study subjects (648 CAD subjects and 284 controls). We constructed a hsa-miR-130a-3p-mediated circRNA-mRNA ceRNA network using the miRanda database. This included 9 circRNAs (hsa_circ_0089378, hsa_circ_0083357, hsa_circ_0082824, hsa_circ_0068942, hsa_circ_0057576, hsa_circ_0054537, hsa_circ_0051172, hsa_circ_0032970, and hsa_circ_0006323) and 1 mRNA (transient receptor potential cation channel subfamily M member 3 [TRPM3]). We have shown that 9 circRNAs promote TRPM3 expression by inhibiting hsa-miR-130a-3p in CAD patients.

Keywords: microRNA, coronary heart disease, circRNA, ceRNA

INTRODUCTION

Cardiovascular disease (CVD) is the leading cause of human morbidity and mortality worldwide. According to World Health Organization (WHO), approximately 17.5 million people die annually from CVD [1]. This represents 31% of all global deaths. Stroke and coronary artery disease (CAD) are the major causes of CVD related deaths [1]. In China, CVD accounts for 300 deaths out of every 100,000 individuals [2], while, nearly a quarter of the Western population suffers from CAD and stroke [34]. Therefore, it is essential to identify risk factors associated with CVD that could be used for diagnosis and treatment at an early stage. The role of non-coding RNAs (ncRNAs) has been recognized in various human pathologies and they represent great potential as CAD biomarkers.

Human genome sequencing shows that although only 3% of the human genome codes for proteins, 80% of the human genome are transcribed [5]. The numbers of non-coding transcripts greatly exceed protein-coding mRNAs and represent species complexity [67]. Long noncoding RNAs (lncRNAs) are noncoding RNAs that are longer than 200 nucleotides that regulate diverse cellular functions. Recently, lncRNAs have been detected in human plasma, which can be used as disease biomarkers [8]. Aberrant levels of lncRNAs have also been reported in plasma from CAD patients [9]. Circular RNAs (circRNAs) are another class of non-coding RNAs that regulate gene expression in eukaryotes [10]. The circRNAs are competing endogenous RNAs (ceRNAs) that act as a sponge for microRNAs (miRNAs) by complementary base paring and therefore regulate gene transcription [1112]. However, the expression and biological functions of ceRNAs in CAD are unknown. Therefore, in the present study, we generated plasma ceRNA expression profiles of three pairs of CAD and control samples to investigate their potential role as diagnostic markers in CAD.

RESULTS

Differentially expressed plasma circRNAs and mRNAs in CAD patients

The circRNA expression profiles in human CAD and control plasma were compared by the scatter plot (Figure (Figure1)1) and volcano plot filtering (Figure (Figure2)2) to identify differentially expressed circRNAs. We identified 24 differentially expressed circRNAs (fold change ≥ 1.5 and P < 0.05) between CAD and control plasma (Table (Table1).1). Figure Figure33 shows the heat map of these 24 differentially expressed circRNAs. Among these, 18 circRNAs were up-regulated and 6 circRNAs were down-regulated in CAD plasma. Among the up-regulated circRNAs, 17 were exonic, and 1 intragenic. Among the down-regulated circRNAs, 6 were exonic and 1 intronic. The circRNA, hsa_circ_0051686 was both intronic and exonic. SBC Human (4*180 K) ceRNA microarray profiling showed 6 up-regulated and 1 down-regulated mRNAs (fold change ≥ 1.5 and P < 0.05) in CAD patients (Table (Table22).

Figure 1
Scatter plot of differentially expressed plasma circRNAs in CAD and control subjects
Figure 2
Volcano plot of the differentially expressed plasma circRNAs in CAD and control subjects
Table 1
List of deregulated cirRNAs in 3 CAD patients
Figure 3
Plasma circRNA profile of CAD and control subjects
Table 2
List of deregulated mRNAs in 3 CAD patients

Differentially expressed miRNAs in CAD

To establish a circRNA-miRNA-mRNA ceRNA network, we compared plasma miRNA expression profiles of CAD and control subjects in another independent population that was previously reported [13]. We identified 9 CAD-related miRNAs including miR-122, miR-133b, miR-214, miR-21, miR-106a, miR-130a, miR-155, miR-221, and miR-125b. Among these, miR-221 (p = 0.001), miR-155 (p = 0.049), and miR-130a (p = 0.001) were downregulated in CAD subjects than in non-CAD subjects (Table (Table33).

Table 3
Differentially expressed miRNAs in the 2nd set of CAD patients

MicroRNA-mRNA interaction

The hsa-miR-221-3p, hsa-miR-155-5p, and hsa-miR-130a-3p were selected for integrated analysis of miRNA and mRNA profiling data. The miRNA target predictions were based on miRanda (release August 2010; http://www.microrna.org/). We performed integrated analysis of the inverse relations of expressed miRNAs and mRNAs in conjunction with their predicted targets and identified transient receptor potential cation channel subfamily M member 3 (TRPM3) as a target gene for hsa-miR-130a-3p.

CircRNA-microRNA interactions

Recently, circRNAs have been identified as miRNA sponges that regulate gene expression. Therefore, we used miRanda database to investigate potential miRNAs that bind to circRNAs in CAD patients. We observed that hsa-miR-221-3p, hsa-miR-155-5p, and hsa-miR-130a-3p were bound by 10, 12, and 9 circRNAs, respectively (Table (Table44).

Table 4
MicroRNA-circRNA interactions in CAD

Construction of ceRNA network

Since circRNAs interact with miRNAs through miRNA response elements (MREs), we searched for putative hsa-miR-130a-3p MREs in the circRNAs using miRanda. We identified 9 circRNAs that had hsa-miR-130a-3p binding sites and negatively associated with hsa-miR-130a-3p. These included hsa_circ_0089378, hsa_circ_0083357, hsa_circ_0082824, hsa_circ_0068942, hsa_circ_0057576, hsa_circ_0054537, hsa_circ_0051172, hsa_circ_0032970, and hsa_circ_0006323. Based on these data, we constructed a hsa-miR-130a-3p-mediated circRNA-mRNA ceRNA network with 9 circRNAs and 1 mRNA (Figure (Figure44 and Table Table55).

Figure 4
Schematic representation of hsa-miR-130a-3p-mediated circRNA-mRNA ceRNA network
Table 5
Hsa-miR-130a-3p-mediated circRNA-mRNA ceRNA network

DISCUSSION

The ceRNAs play a key role in post-transcriptional regulation and have been implicated in cardiovascular disease [12, 1415]. In our previous study, we reported several miRNAs associated with CAD, but their actions were unknown [13]. In order to deciphere the ceRNA mechanisms related to CAD, we constructed a global ceRNA-miRNA-mRNA triple network using the miRanda database. We identified 9 circRNAs and 1 mRNA that formed a network with hsa-miR-130a-3p. This study sheds new insights to exploring the complex post-transcriptional regulatory networks via ceRNA interactions and identifying pathways that are altered in CAD. Subsequently, specific ceRNAs can be therapeutically used to modulate specific pathways that are involved in CAD pathology.

The miR-130 precursor is a small non-coding RNA that has been identified in mice (MI0000156, MI0000408), humans (MI0000448, MI0000748) and a range of vertebrate species (MIPF0000034). Mature miR-130 is generated by the dicer enzyme upon excision from the 3′ arm of the hairpin. MiR-130a-3p is inversely associated with coronary atherosclerosis with its down-regulation contributing to endothelial progenitor cell dysfunction in subjects suffering from coronary artery disease [1618]. Therefore, miR-130a-3p has therapeutic potential for the prevention and treatment of CAD. Its regulation by ceRNAs as shown in this study represents one possible mode of regulation that is clinically applicable.

TRPM3 belongs to the family of transient receptor potential (TRP) channels that regulate cellular calcium homeostasis. TRPM3 mediates calcium entry potentiated by calcium store depletion. Alternatively spliced transcript variants encoding different isoforms of TRPM3 have also been identified. TRPM3 regulates proliferation and contractility of vascular smooth muscle cells in co-ordination with cholesterol and is potentially involved in therapeutic vascular modulation [19]. In the present study, we identified TRPM3 mRNA as a hsa-miR-130a-3p target and is upregulated in CAD subjects.

Circular RNAs (circRNAs) are new members of ceRNAs that are involved in regulating gene expression [20]. However, their role in CAD pathogenesis has not been reported. In the present study, we identified 24 aberrantly expressed circRNAs (18 up-regulated and 6 down-regulated) in CAD patients. These included 9 circRNAs,namely, hsa_circ_0089378, hsa_circ_0083357, hsa_circ_0082824, hsa_circ_0068942, hsa_circ_0057576, hsa_circ_0054537, hsa_circ_0051172, hsa_circ_0032970, and hsa_circ_0006323, which sponge hsa-miR-130a-3p that regulates TRPM3. This suggests that the cohort of circRNAs negatively regulate miR-130A-3p, thereby resulting in upregulation of TRPM3.

This study has several limitations that need to be addressed while interpreting our results. First, the ceRNA microarray was based on a small group of 3 CAD and 3 control subjects. Since the expression of circRNA and ceRNA can vary in individuals due to a number of factors, the circRNA and mRNA profiles were verified in another independent cohort of 932 subjects by qRT-PCR. The in vivo relevance of our findings requires further comprehensive investigation. Second, the mechanism of circRNA regulation of hsa-miR-130a-3p and TRPM3 is based on the bioinformatics analysis and needs to be tested in in vitro and in vivo models.

In conclusion, we identified a network of 9 circRNAs that regulate TRPM3 expression by inhibiting hsa-miR-130a-3p in CAD patients.

MATERIALS AND METHODS

Study subjects

In the first group, we enrolled 3 CAD and 3 control (4 males and 2 females) subjects at the Yancheng TCM Hospital Affiliated to Nanjing University of Chinese Medicine in China. The study was performed as approved by the ethics committee of the Yancheng TCM Hospital Affiliated to Nanjing University of Chinese Medicine and the First Affiliated Hospital of Nanjing Medical University. All subjects provided written informed consent. The characteristics of the study subjects are shown in Table Table66.

Table 6
Characteristics of the CAD study population

In the second group, 932 consecutive adult subjects (681 males and 251 females, 648 CAD subjects and 284 controls) aged 32–84 years that underwent coronary angiography for suspected or known coronary atherosclerosis was analyzed. They were part of a previous study on CAD. Subjects with spastic angina pectoris, infectious processes within 2 weeks, heart failure, adrenal dysfunction, and thyroid dysfunction were excluded from this study. The plasma miRNA levels were confirmed by qRT-PCR analysis. Circulating levels of miRNAs were quantified using the 2−Δct method.

Microarray analysis

Total RNA was extracted from the plasma of the subjects and purified using the mirVanaTM PARISTM kit (Cat#AM1556, Ambion, Austin, TX, USA) according to the manufacturer's instructions. The RNA integration number (RIN) was determined by an Agilent Bioanalyzer 2100 (Agilent technologies, Santa Clara, CA, USA). Total RNA (1μg) was amplified and labeled by Low Input Quick Amp WT Labeling Kit (Cat. # 5190-2943, Agilent technologies, Santa Clara, CA, USA) according to the manufacturer's instructions. Labeled cRNA were purified by RNeasy mini kit (Cat.# 74106, QIAGEN, GmBH, Germany).

The cDNA was labeled and hybridized to the human SBC-ceRNA (4×180k) Array, which can detect 88,371 circRNAs, and 18,853 coding transcripts (Circbase (88371), GENCODE v21 /Ensembl (18,100), LNCipedia v3.1 (40,621), Lncrnadb (28), Noncode v4 (2,608), UCSC (25,919) databases). We hybridized 1μg Cy3-labeled cRNA onto each slide using Gene Expression Hybridization Kit (Cat.# 5188-5242, Agilent technologies, Santa Clara, CA, US) in an hybridization oven (Cat.# G2545A, Agilent technologies, Santa Clara, CA, US) according to the manufacturer's instructions for 17 h. Then, slides were washed in staining dishes (Cat.# 121, Thermo Shandon, Waltham, MA, USA) with Gene Expression Wash Buffer Kit (Cat.# 5188-5327, Agilent technologies, Santa Clara, CA, USA) according to the manufacturer's instructions.

The slides were scanned by Agilent Microarray Scanner (Cat#G2565CA, Agilent technologies, Santa Clara, CA, USA) with default settings (Dye channel: Green; Scan resolution = 3 μm; PMT 100%; 20 bit). Raw data was extracted with Feature Extraction software 10.7 (Agilent technologies, Santa Clara, CA, USA) and normalized by Quantile algorithm and LIMMA packages in R. The microarray work was performed by Shanghai Biotechnology Cooperation, Shanghai, P.R. China.

Construction of ceRNA network

Figure Figure55 shows the methodology used to identify the ceRNA interacting genes. The miRanda method, which is based on dynamic programming (SW algorithm) and computing free energy was used to identify the target genes of miRNAs and ceRNAs [2125]. Based on this analysis, we built a microRNA-cirRNA-mRNA interaction network. The relationship between the target and the microRNA was based on the adjacency matrix of microRNA and target A = [ai,j], where ai,j represents the weight of the relationship between the target (i) and its microRNA (j). In the microRNA-cirRNA-mRNA network, the circle represents one edge, whereas the center of the network represents a degree. The degree denotes the contribution of a microRNA to the target gene around it or the contribution of a target to the microRNAs around it. The key microRNAs and their targets in the network always had the largest degrees.

Figure 5
Flow chart of CAD related competing endogenous RNA network study

ACKNOWLEDGMENTS AND FUNDING

This study was supported by grants from the National Natural Science Foundations of China (grant numbers 81170180, 30400173 and 30971257) and the Priority Academic Program Development of Jiangsu Higher Education Institutions.

Footnotes

CONFLICTS OF INTEREST

The authors declare that no conflicts of interest exist.

REFERENCES

1. WHO Global Status Report on noncommunicable diseases 2014. Available at: http://www.thehealthwell.info/node/831799 Accessed December 18 2015.
2. WHO Cardiovascular diseases mortality: age-standardized death rate per 100,000 population, 2000–2012. Available at http://gamapserver.who.int/gho/interactive_charts/ncd/mortality/cvd/atlas.html Accessed: December 18, 2015.
3. Mozaffarian D, Benjamin EJ, Go AS, Arnett DK, Blaha MJ, Cushman M, de Ferranti S, Després JP, Fullerton HJ, Howard VJ, Huffman MD, Judd SE, Kissela BM, et al. Heart disease and stroke statistics—2015 update: a report from the American Heart Association. Circulation. 2015;131:e29–322. [PubMed]
4. Leal J, Luengo-Fernandez R, Gray A, Rayner M. Economic burden of cardiovascular diseases in the enlarged European Union. Eur Heart J. 2006;27:1610–9. [PubMed]
5. Genomics Pennisi E. ENCODE project writes eulogy for junk DNA. Science. 2012;337:1159–61. [PubMed]
6. Liu G, Mattick JS, Taft RJ. A meta-analysis of the genomic and transcriptomic composition of complex life. Cell Cycle. 2013;12:2061–72. [PMC free article] [PubMed]
7. Arner E, Daub CO, Vitting-Seerup K, Andersson R, Lilje B, Drabløs F, Lennartsson A, Rönnerblad M, Hrydziuszko O, Vitezic M, Freeman TC, Alhendi AM, Arner P, et al. Transcribed enhancers lead waves of coordinated transcription in transitioning mammalian cells. Science. 2015;347:1010–4. [PMC free article] [PubMed]
8. Schlosser K, Hanson J, Villeneuve PJ, Dimitroulakos J, McIntyre L, Pilote L, Stewart DJ. Assessment of Circulating LncRNAs Under Physiologic and Pathologic Conditions in Humans Reveals Potential Limitations as Biomarkers. Sci Rep. 2016;6:36596. https://doi.org/10.1038/srep36596 [PMC free article] [PubMed]
9. Yang Y, Cai Y, Wu G, Chen X, Liu Y, Wang X, Yu J, Li C, Chen X, Jose PA, Zhou L, Zeng C. Plasma long non-coding RNA, CoroMarker, a novel biomarker for diagnosis of coronary artery disease. Clin Sci (Lond) 2015;129:675–85. [PubMed]
10. Chen J, Li Y, Zheng Q, Bao C, He J, Chen B, Lyu D, Zheng B, Xu Y, Long Z, Zhou Y, Zhu H, Wang Y, et al. Circular RNA profile identifies circPVT1 as a proliferative factor and prognostic marker in gastric cancer. Cancer Lett. 2016 [PubMed]
11. Wang K, Long B, Liu F, Wang JX, Liu CY, Zhao B, Zhou LY, Sun T, Wang M, Yu T, Gong Y, Liu J, Dong YH, et al. A circular RNA protects the heart from pathological hypertrophy and heart failure by targeting miR-223. Eur Heart J. 2016;37:2602–11. https://doi.org/10.1093/eurheartj/ehv713 [PubMed]
12. Song C, Zhang J, Liu Y, Pan H, Qi HP, Cao YG, Zhao JM, Li S, Guo J, Sun HL, Li CQ. Construction and analysis of cardiac hypertrophy-associated lncRNA-mRNA network based on competitive endogenous RNA reveal functional lncRNAs in cardiac hypertrophy. Oncotarget. 2016;7:10827–40. https://doi.org/10.18632.oncotarget.7312 [PMC free article] [PubMed]
13. Jia QW, Chen ZH, Ding XQ, Liu JY, Ge PC, An FH, Li LH, Wang LS, Ma WZ, Yang ZJ, Jia EZ. Predictive Effects of Circulating miR-221, miR-130a and miR-155 for Coronary Heart Disease: A Multi-Ethnic Study in China. Cell Physiol Biochem. 2017;42:808–823. [PubMed]
14. Pang L, Hu J, Zhang G, Li X, Zhang X, Yu F, Lan Y, Xu J, Pang B, Han D, Xiao Y, Li X. Dysregulated long intergenic non-coding RNA modules contribute to heart failure. Oncotarget. 2016;7:59676–59690. https://doi.org/10.18632/oncotarget.10834 [PMC free article] [PubMed]
15. Ling S, Birnbaum Y, Nanhwan MK, Thomas B, Bajaj M, Li Y, Li Y, Ye Y. Dickkopf-1 (DKK1) phosphatase and tensin homolog on chromosome 10 (PTEN) crosstalk via microRNA interference in the diabetic heart. Basic Res Cardiol. 2013;108:352. https://doi.org/10.1007/s00395-013-0352-2 [PubMed]
16. de Gonzalo-Calvo D, Cenarro A, Garlaschelli K, Pellegatta F, Vilades D, Nasarre L, Camino-Lopez S, Crespo J, Carreras F, Leta R, Catapano AL, Norata GD, Civeira F, et al. Translating the microRNA signature of microvesicles derived from human coronary artery smooth muscle cells in patients with familial hypercholesterolemia and coronary artery disease. J Mol Cell Cardiol. 2017;106:55–67. https://doi.org/10.1016/j.yjmcc.2017.03.005 [PubMed]
17. Staszel T, Zapała B, Polus A, Sadakierska-Chudy A, Kieć-Wilk B, Stępień E, Wybrańska I, Chojnacka M, Dembińska-Kieć A. Role of microRNAs in endothelial cell pathophysiology. Pol Arch Med Wewn. 2011;121:361–6. [PubMed]
18. Zhang Q, Kandic I, Kutryk MJ. Dysregulation of angiogenesis-related microRNAs in endothelial progenitor cells from patients with coronary artery disease. Biochem Biophys Res Commun. 2011;405:42–6. https://doi.org/10.1016/j.bbrc.2010.12.119 [PubMed]
19. Naylor J, Li J, Milligan CJ, Zeng F, Sukumar P, Hou B, Sedo A, Yuldasheva N, Majeed Y, Beri D, Jiang S, Seymour VA, McKeown L, et al. Pregnenolone sulphate- and cholesterol-regulated TRPM3 channels coupled to vascular smooth muscle secretion and contraction. Circ Res. 2010;106:1507–15. [PMC free article] [PubMed]
20. Zhou B, Xu HY. A novel identified circular RNA, circRNA_010567, promotes myocardial fibrosis via suppressing miR-141 by targeting TGF-β1. Biochem Biophys Res Commun. 2017 https://doi.org/10.1016/j.bbrc.2017.04.044 [PubMed]
21. Lewis BP, Burge CB, Bartel DP. Conserved Seed Pairing, Often Flanked by Adenosines, Indicates that Thousands of Human Genes are MicroRNA Targets. Cell. 2005;120:15–20. [PubMed]
22. Friedman RC, Farh KK, Burge CB, Bartel DP. Most Mammalian mRNAs Are Conserved Targets of MicroRNAs. Genome Research. 2009;19:92–105. [PubMed]
23. Grimson A, Farh KK, Johnston WK, Garrett-Engele P, Lim LP, Bartel DP. MicroRNA Targeting Specificity in Mammals: Determinants beyond Seed Pairing. Molecular Cell. 2007;27:91–105. [PMC free article] [PubMed]
24. Garcia DM, Baek D, Shin C, Bell GW, Grimson A, Bartel DP. Weak Seed-Pairing Stability and High Target-Site Abundance Decrease the Proficiency of lsy-6 and Other miRNAs. Nat Struct Mol Biol. 2011;18:1139–1146. [PMC free article] [PubMed]
25. Myers EW, Miller W. Approximate matching of regular expressions. Bull Math Biol. 1989;51:5–37. [PubMed]

Articles from Oncotarget are provided here courtesy of Impact Journals, LLC