- Open Access
Construction of microRNA and transcription factor regulatory network based on gene expression data in cardiomyopathy
European Journal of Medical Research volume 19, Article number: 57 (2014)
Cardiomyopathy is a progressive myocardial disorder. Here, we attempted to reveal the possible mechanism of cardiomyopathy at the transcription level with the roles of microRNAs (miRNAs) and transcription factors (TFs) taken into account.
We firstly identified differentially expressed genes (DEGs) between cardiomyopathy patients and controls with data from the gene expression omnibus (GEO) database. DEGs were associated with the canonical pathways, molecular and cellular functions, physiological system development and function in the Ingenuity Knowledge Base by using the Ingenuity Pathway Analysis (IPA) software. TFs and miRNAs that DEGs significantly enriched were identified and a double-factor regulatory network was constructed.
A total of 1,680 DEGs were identified. The DEGs were enriched for various pathways, with glucocorticoid receptor signaling as the most significant. A double-factor regulatory network was constructed, including seven TFs and two miRNAs. A subnetwork under the regulation of MEF2C and SRF was also constructed to illustrate their regulatory effects on cardiac functions.
Our results may provide new understanding of cardiomyopathy and may facilitate further therapeutic studies.
Cardiomyopathy is a progressive myocardial disorder, usually leading to cardiovascular death or heart failure-related disability . It develops at any age, in either sex, and in any population ,. The etiology of cardiomyopathy is highly complex and improving its treatments has become a research hotspot.
Using bioinformatics combined with gene expression data to identify potential therapeutic targets has shown great application prospects. The majority of the previous studies mainly focused on the analyses of differentially expressed genes (DEGs), without considering microRNAs (miRNAs) and transcription factors (TFs) that regulate the expression of DEGs. miRNAs are small non-coding RNAs that control various biological processes through affecting the stability and translation of target mRNAs. Previous studies have proposed several miRNAs as being involved in the pathogenesis of cardiomyopathy, such as miR-1  and miR-21 . TFs can regulate gene expression through binding to the cis-elements in target genes' promoter regions. TFs, such as MEF2, have been reported to be associated with cardiomyopathy . Considering the important regulatory roles of miRNAs and TFs in the pathogenesis of cardiomyopathy, identification of miRNAs and TFs that enriched with target DEGs and construction of a double-factor regulatory network may provide new understanding of the molecular mechanism of cardiomyopathy.
In the current study, based on gene expression data from the Gene Expression Omnibus (GEO) database, we acquired DEGs, miRNAs and TFs that enriched with target DEGs and constructed a double-factor regulatory network. Our results may reveal the possible mechanism of cardiomyopathy at the transcription level.
This study was approved by the institutional review board of the Xi'an Children's Hospital (20140518). Written informed consent was obtained from all patients for the publication of this report and any accompanying images.
The gene expression profile GSE5406 from the GEO database was used. This dataset includes transcription profiles of 210 left ventricular myocardial tissue samples, 194 of which were from patients with advanced cardiomyopathy and 16 of which were from healthy donors. All patients had New York Heart Association class 3 to 4 symptoms and left ventricular systolic dysfunction, with ejection fraction of 14 ± 8% (mean ± SD). Patients suffered from heart failure due to ischemic (n =86) or idiopathic dilated (n =108) cardiomyopathy. Control samples had normal left ventricular function with ejection fraction of 56 ± 7% (P =0.0001 versus patients). None of the subjects received mechanical support with left ventricular assist devices. Myocardial tissue samples were obtained from patients undergoing heart transplantation and from controls deemed unsuitable for transplantation. Whole hearts were removed after preservation in cold cardioplegia at the time of transplantation or donor harvest. Then, segments of noninfarcted left ventricular free wall were snap-frozen in liquid nitrogen. For each sample, RNA was isolated by using Trizol reagent (Invitrogen, Carlsbad, CA, USA). The dataset was generated by using the (HG-U133A) Affymetrix Human Genome U133A Array (Affymetrix, Santa Clara, CA, USA).
Identification of differentially expressed genes (DEGs)
Raw data from all arrays were normalized using Robust Multi-array Analysis (RMA)  in the R software (version 3.0.0). The resulting expression values were used to identify DEGs with the limma package (3.12.1) in R. DEGs were detected by using t-tests and multiple test corrections were carried out with the Benjamini-Hochberg method . The threshold for significance was set as P <0.01.
To explore the functions and pathways of DEGs, DEGs were associated with the canonical pathways, molecular and cellular functions, physiological system development and function in the Ingenuity Knowledge Base by using the Ingenuity Pathway Analysis (IPA) software (Ingenuity® Systems, http://www.ingenuity.com).
Construction of miRNA-TF regulatory network
We also acquired miRNAs and TFs which were overrepresented with target DEGs. For miRNA analysis, we used Targetscan , miRanda  and Pita  for miRNA target DEGs prediction. To avoid false positive results, miRNA-DEG prediction results supported by all three prediction methods were considered to be confidential. For TF analysis, we obtained TF binding sites and the coordinate's position information of human being from the University of California Santa Cruz (UCSC) database. TFs and miRNAs that DEGs significantly enriched were identified by using the hypergeometric distribution test. P <0.01 was set as the threshold. Based on the regulatory relationship between miRNAs or TFs and DEGs, a network was constructed by using the IPA software.
A total of 1,680 DEGs were identified, including 963 down-regulated genes and 717 up-regulated ones in advanced cardiomyopathy patients.
As shown in Figure 1, the IPA results of canonical pathways showed that the DEGs were enriched for various pathways, among which glucocorticoid receptor signaling was the most significant. Molecular and cellular functions analysis revealed that `Cell Death and Survival' (P =2.53 × 10–24 to 1.36 × 10–3) was the top molecular function affected by DEGs followed by `Cellular Growth and Proliferation' (P =3.37 × 10–24 to 1.36 × 10–3) (Table 1). `Cardiovascular System Development and Function' (P =4.17 × 10–12 to 1.20 × 10–3) was the top physiological function mediated by DEGs (Table 1).
The results of miRNA and TF enrichment analysis are listed in Table 2. A network was constructed to illustrate the regulatory relationships (Figure 2). Considering the important role of SRF and MEF2C in the progression of cardiomyopathy, a subnetwork under their regulation was constructed and the biological functions they may affect were also indicated (Figure 3).
In the current study, based on GSE5406 from the GEO database, 1,680 DEGs were identified between advanced cardiomyopathy patients and healthy controls. DEGs were mapped into the Ingenuity Knowledge Base and several canonical pathways were screened out, with glucocorticoid receptor signaling as the most significant. Glucocorticoid receptors are involved in cardiovascular homeostasis  and glucocorticoid could protect rodent hearts from ischemia/reperfusion injury . Our results confirmed the essential roles of this pathway in the pathogenesis of cardiomyopathy. Molecular and cellular functions analysis revealed that `Cell Death and Survival' (P =2.53 × 10–24 to 1.36 × 10–3) was the top molecular function affected by DEGs, suggesting that dysregulation of the cell function may contribute to the progressive heart failure induced by cardiomyopathy.
TFs or miRNAs that DEGs significantly enriched were identified and a double-factor regulatory network was constructed, including seven TFs and two miRNAs (Figure 2). All these miRNAs and TFs have previously been implicated in cardiomyogenesis or cardiac function. Considering their important regulatory roles in controlling various biological processes, these miRNAs and TFs may provide new avenues for the therapeutic strategies of cardiomyopathy. In addition, since all these miRNAs and TFs were detected in patients suffered from New York Heart Association class 3 to 4 symptoms, they may also be considered as prognosis markers in clinical practice. A subnetwork under the regulation of SRF and MEF2C was also constructed and the functions they may affect were also indicated (Figure 3) considering the considerable evidence in support of their involvement in cardiomyopathy -. As shown in Figure 3, both SRF and MEF2C may inhibit the expression of MYH6 and CNN1. MYH6 is a famous cardiac muscle myosin gene and mutations of this gene are implicated in cardiomyopathy in various studies . Inhibition of CNN1 expression may disrupt its suppression of cardiomyopathy through the εPKC pathway . In addition, both SRF and MEF2C may regulate the expression of NPPA, which was overexpressed in patients in our results. This gene is involved in many cardiac dysfunctions (Figure 3), including myocardial infarction, heart failure, and high ventricular pressure ,.
In summary, with DEGs identified by gene chip data generated by high-throughput technologies from the GEO database, we constructed a cardiomyopathy double-factor regulatory network and revealed the possible mechanism on regulation level. A subnetwork under the regulation of SRF and MEF2C was also constructed to illustrate their possible regulatory roles. Our results may facilitate further therapeutic studies of cardiomyopathy.
LW conceived of the study and drafted the manuscript. JH, HX and MS performed the statistical analysis. JW, QJ and HY revised the manuscript. All authors read and approved the final manuscript.
differentially expressed genes
Gene Expression Omnibus
Ingenuity Pathway Analysis
Robust Multi-array Analysis
University of California Santa Cruz
Maron BJ, Towbin JA, Thiene G, Antzelevitch C, Corrado D, Arnett D, Moss AJ, Seidman CE, Young JB: Contemporary definitions and classification of the cardiomyopathies: an American Heart Association scientific statement from the council on clinical cardiology, heart failure and transplantation committee; quality of care and outcomes research and functional genomics and translational biology interdisciplinary working groups; and council on epidemiology and prevention. Circulation 2006, 113: 1807–1816. 10.1161/CIRCULATIONAHA.106.174287
Towbin JA, Lowe AM, Colan SD, Sleeper LA, Orav EJ, Clunie S, Messere J, Cox GF, Lurie PR, Hsu D, Canter C, Wilkinson JD, Lipshultz SE: Incidence, causes, and outcomes of dilated cardiomyopathy in children. JAMA 2006, 296: 1867–1876. 10.1001/jama.296.15.1867
Taylor MR, Carniel E, Mestroni L: Cardiomyopathy, familial dilated. Orphanet J Rxare Dis 2006, 1: 27. 10.1186/1750-1172-1-27
Ferreira LR, Frade AF, Santos RH, Teixeira PC, Baron MA, Navarro IC, Benvenuti LA, Fiorelli AI, Bocchi EA, Stolf NA, Chevillard C, Kalil J, Cunha-Neto E: MicroRNAs miR-1, miR-133a, miR-133b, miR-208a and miR-208b are dysregulated in chronic Chagas disease cardiomyopathy. Int J Cardiol 2014, 175: 409–417. 10.1016/j.ijcard.2014.05.019
Xu HF, Ding YJ, Zhang ZX, Wang ZF, Luo CL, Li BX, Shen YW, Tao LY, Zhao ZQ: MicroRNA21 regulation of the progression of viral myocarditis to dilated cardiomyopathy. Mol Med Reports 2014, 10: 161–168.
Hannenhalli S, Putt ME, Gilmore JM, Wang J, Parmacek MS, Epstein JA, Morrisey EE, Margulies KB, Cappola TP: Transcriptional genomics associates FOX transcription factors with human heart failure. Circulation 2006, 114: 1269–1276. 10.1161/CIRCULATIONAHA.106.632430
Irizarry RA, Hobbs B, Collin F, Beazer-Barclay YD, Antonellis KJ, Scherf U, Speed TP: Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostatistics 2003, 4: 249–264. 10.1093/biostatistics/4.2.249
Benjamini Y, Hochberg Y: Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser B Methodol 1995, 57: 289–300.
Grimson A, Farh KK, Johnston WK, Garrett-Engele P, Lim LP, Bartel DP: MicroRNA targeting specificity in mammals: determinants beyond seed pairing. Mol Cell 2007, 27: 91–105. 10.1016/j.molcel.2007.06.017
John B, Enright AJ, Aravin A, Tuschl T, Sander C, Marks DS: Human MicroRNA targets. PLoS Biol 2004, 2: e363. 10.1371/journal.pbio.0020363
Kertesz M, Iovino N, Unnerstall U, Gaul U, Segal E: The role of site accessibility in microRNA target recognition. Nat Genet 2007, 39: 1278–1284. 10.1038/ng2135
Rickard AJ, Funder JW, Morgan J, Fuller PJ, Young MJ: Does glucocorticoid receptor blockade exacerbate tissue damage after mineralocorticoid/salt administration? Endocrinology 2007, 148: 4829–4835. 10.1210/en.2007-0209
Tokudome S, Sano M, Shinmura K, Matsuhashi T, Morizane S, Moriyama H, Tamaki K, Hayashida K, Nakanishi H, Yoshikawa N, Shimizu N, Endo J, Katayama T, Murata M, Yuasa S, Kaneda R, Tomita K, Eguchi N, Urade Y, Asano K, Utsunomiya Y, Suzuki T, Taguchi R, Tanaka H, Fukuda K: Glucocorticoid protects rodent hearts from ischemia/reperfusion injury by activating lipocalin-type prostaglandin D synthase-derived PGD2 biosynthesis. J Clin Invest 2009, 119: 1477–1488. 10.1172/JCI37413
Alonso-Montes C, Naves-Diaz M, Fernandez-Martin JL, Rodriguez-Reguero J, Moris C, Coto E, Cannata-Andia JB, Rodriguez I: New polymorphisms in human MEF2C gene as potential modifier of hypertrophic cardiomyopathy. Mol Biol Rep 2012, 39: 8777–8785. 10.1007/s11033-012-1740-7
Xu J, Gong NL, Bodi I, Aronow BJ, Backx PH, Molkentin JD: Myocyte enhancer factors 2A and 2C induce dilated cardiomyopathy in transgenic mice. J Biol Chem 2006, 281: 9152–9162. 10.1074/jbc.M510217200
Dasgupta T, Stillwagon SJ, Ladd AN: Gene expression analyses implicate an alternative splicing program in regulating contractile gene expression and serum response factor activity in mice. PLoS One 2013, 8: e56590. 10.1371/journal.pone.0056590
Tritsch E, Mallat Y, Lefebvre F, Diguet N, Escoubet B, Blanc J, De Windt LJ, Catalucci D, Vandecasteele G, Li Z, Mericskay M: An SRF/miR-1 axis regulates NCX1 and annexin A5 protein levels in the normal and failing heart. Cardiovasc Res 2013, 98: 372–380. 10.1093/cvr/cvt042
Jiang J, Wakimoto H, Seidman JG, Seidman CE: Allele-specific silencing of mutant Myh6 transcripts in mice suppresses hypertrophic cardiomyopathy. Science 2013, 342: 111–114. 10.1126/science.1236921
Lu D, Zhang L, Bao D, Lu Y, Zhang X, Liu N, Ge W, Gao X, Li H, Zhang L: Calponin1 inhibits dilated cardiomyopathy development in mice through the epsilonPKC pathway. Int J Cardiol 2014, 173: 146–153. 10.1016/j.ijcard.2014.02.032
Asanuma H, Kitakaze M: Carperitide and adiponectin - how are they connected each other to benefit acute decompensated heart failure? Circ J 2009, 73: 2206–2207. http://www.ncbi.nlm.nih.gov/pubmed/19940369 10.1253/circj.CJ-09-0785
Tarazon E, Rosello-Lleti E, Rivera M, Ortega A, Molina-Navarro MM, Trivino JC, Lago F, Gonzalez-Juanatey JR, Orosa P, Montero JA, Salvador A, Portolés M: RNA sequencing analysis and atrial natriuretic peptide production in patients with dilated and ischemic cardiomyopathy. PLoS One 2014, 9: e90157. 10.1371/journal.pone.0090157
This study is supported by Provincial Natural Science Basic Research Foundation of Shaanxi (2014JM4152) and the Science and Technology Project of Xi'an Municipal Health Bureau (Number 2013028).
The authors declare that they have no competing interests.
About this article
Cite this article
Wang, L., Hu, J., Xing, H. et al. Construction of microRNA and transcription factor regulatory network based on gene expression data in cardiomyopathy. Eur J Med Res 19, 57 (2014). https://doi.org/10.1186/s40001-014-0057-5
- Transcription factors
- Gene expression