Expression profile analysis based on DNA microarray for patients undergoing off‑pump coronary artery bypass surgery
- Authors:
- Published online on: January 18, 2016 https://doi.org/10.3892/etm.2016.3003
- Pages: 864-872
-
Copyright: © Sun et al. This is an open access article distributed under the terms of Creative Commons Attribution License.
Abstract
Introduction
Coronary artery disease is a defect with a narrowing of the coronary arteries that supply blood and oxygen to the heart (1). Off-pump coronary artery bypass (OPCAB) surgery, in which veins or arteries are grafted to the coronary arteries to bypass atherosclerotic narrowing and improve the blood supply to the heart muscle, reduces the risk of mortality from coronary artery disease (2). However, coronary patients who have undergone OPCAB grafting often suffer from disorders of coagulation and hemostasis (3). Therefore, it is necessary to research the molecular mechanisms underlying the effect of OPCAB.
Expression profile analyses have been used to examine the myocardial stress response to cardiac surgery (4–6). In addition, molecular biology studies have confirmed that several genes, including interleukin 6 (IL6), IL8 and activating transcription factor 3 (ATF3), and certain biological processes are associated with the effect of OPCAB grafting. IL6 and IL8 proteins, released by ischemic cardiac myocytes, are induced following acute ischemia (5,7). Moreover, Tomic et al found that the expression of IL6 tended to increase after OPCAB grafting (8), and Ghorbel et al showed that ventricular heart cells may undergo alterations in chemokine and plasma cytokine levels following OPCAB surgery (9). Furthermore, several important miRNA regulatory relationships in OPCAB samples have been investigated. For example, microRNA (MIR)-494 targets the 3′ untranslated region of ATF3 reporter and decreases ATF3 mRNA expression (10,11).
The identification of sensitive genes and miRNAs has become the focus for establishing suitable therapeutic strategies, and gene expression profiling techniques have been widely used to analyze the effects of OPCAB. A number of studies have reported some specific gene expression changes in samples following OPCAB surgery (9,12,13), while a variety of genes have been studied in the fundamental research of OPCAB. Studies have screened for genes associated with the levels of chemokines and cytokine following OPCAB surgery, but little attention has been paid to the corresponding miRNA transcripts (8). In addition, the molecular mechanisms of complications, including hemostasis and coagulation disorders, caused by OPCAB grafting are unclear (14).
In the present study, the data in the gene expression series record GSE12486 were downloaded, to screen for differentially expressed genes (DEGs) between samples obtained pre- and post-OPCAB surgery, followed by function and pathway enrichment analysis, as well as protein-protein interaction (PPI) network construction. Based on expression profiles obtained by microarray analysis, Gene Set Enrichment Analysis (GSEA) was used to enrich the miRNAs predicted to have a correlation with OPCAB grafting. Finally, miRNA regulatory networks were constructed to investigate the association between miRNA and DEGs.
Materials and methods
Microarray data
The gene expression profile GSE12486 was obtained from the Gene Expression Omnibus database (GEO; http://www.ncbi.nlm.nih.gov/geo/) (9). The platform was GPL570 Affymetrix Human Genome U133 Plus 2.0 Array (Affymetrix, Inc., Santa Clara, CA, USA). A total of 10 myocardial samples obtained prior to and after grafting were available; they were collected from 5 patients undergoing OPCAB grafting with cardiopulmonary bypass and cardiac arrest. The pre-surgery samples contained GSM313629, GSM313631, GSM313633, GSM313635 and GSM313637, while post-surgery samples included GSM313638, GSM313639, GSM313640, GSM313641 and GSM313642.
Data preprocessing and DEG screening
The original data were read by the Affymetrix package with R-based software and processed by the Robust Multi-array Analysis (RMA) method (15). The processed data then underwent normalization with a universal background in order to evaluate the expression values. The limma multiple linear regression package (16) was applied for the calculation and analysis of DEGs, with further rectification by the Bayes method (17). The genes with P<0.01 and |log (fold change)|>1 were selected as DEGs.
Principal component analysis (PCA)
The PCA was conducted on the basis of the expression profile of DEGs. PCA is a multivariate statistical analysis that uses linear transformation to convert multiple correlated variables into a set of important linearly uncorrelated variables (18).
Pathway enrichment analysis
Database for Annotation, Visualization and Integrated Discovery (DAVID) is a web-based tool for extracting the biological meaning of large numbers of genes (19). It was used to carry out Gene Ontology (GO) function and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis of upregulated DEGs. Based on the principle of hypergeometric distribution, GO and KEGG terms were enriched with the threshold of P<0.01.
Protein-protein interaction (PPI) network construction
A PPI network of DEGs was constructed using the Search Tool for Retrieval of Interacting Genes (STRING) database, which provides integrated knowledge of the known and predicted associations for protein networks (20). The PPI network comprised ‘nodes’ and ‘edges’; each protein was represented by a node, while the interaction of pairwise proteins was shown by an edge. The 5 nodes with the greatest number of connections (high degree) were further screened as the ‘hub’ genes.
miRNA prediction analysis
GSEA is a method for the enrichment analysis of a gene set based on whole genome expression profiling (21). On the basis of the expression profiles obtained by microarray analysis, GSEA was used to enrich the disease-related miRNAs. The related miRNAs were screened with a threshold of false discovery rate (FDR) of <0.05.
miRNA regulatory networks
Cytoscape is an open software platform used for visualizing complex networks and integrating these networks with any type of attribute data (22). It was applied for building the regulatory networks of DEGs that were regulated by miRNA.
Results
Data standardization and DEG screening
After processing by the RMA method, the expression data presented good standardization in the samples taken pre- and post-surgery (Fig. 1). On the basis of the differential expression analysis for GSE12486, a total of 64 DEGs with |log2FC| >1 and P<0.01 were obtained, including 63 upregulated genes and 1 downregulated gene. Perilipin 1 (PLIN1) was the only downregulated gene. The top 10 DEGs are presented in Table I.
PCA
The DEGs were processed by the PCA method. GSM313629 and GSM313640 samples were the first principal components pre- and post grafting, respectively. As shown in Fig. 2, the first two principal components (score, 84.88%) were well able to distinguish between the samples taken pre- and post-surgery.
GO and KEGG pathway enrichment analysis
GO enrichment analysis of upregulated and downregulated DEGs was performed using the DAVID tool. In total, 20 GO terms were found to be enriched by the upregulated DEGs (Table II), of which several important functions were associated with immune function, such as inflammatory response, defense response, response to wounding and chemokine receptor binding. IL6, IL8, chemokine (C-C motif) ligand 4 (CCL4), CCL3, CCL2, chemokine (C-X-C motif) ligand 2 (CXCL2) and nuclear factor of κ light polypeptide gene enhancer in B-cells inhibitor, ζ (NFKBIZ) were among the DEGs enriched in these significant functions.
The DAVID tool was also used for KEGG pathway analysis of the screened upregulated DEGs. Only 5 pathways were enriched (Table III), including the chemokine signaling pathway (P=2.23E-05), cytokine-cytokine receptor interaction (P=2.24E-05), Toll-like receptor signaling pathway (P=1.08E-04), NOD-like receptor signaling pathway (P=2.09E-04) and mitogen-activated protein kinase (MAPK) signaling pathway (P=0.008586). The most significant pathway was the chemokine signaling pathway, which was enriched by the following upregulated DEGs: CCL3, CCL2, IL8, CXCR4, CCL3L1, CXCL2, CCL3L3, CCL8 and CCL4.
PPI network construction
The PPI network, which contained 55 nodes and 234 edges, was constructed on the basis of the STRING database analysis (Fig. 3). Based on the degree (connectivity) of the nodes, the top 10 upregulated hub nodes with higher degrees were screened. These hub genes included FBJ murine osteosarcoma viral oncogene homolog (FOS), jun proto-oncogene (JUN), ATF3, IL6, FBJ murine osteosarcoma viral oncogene homolog B (FOSB), jun B proto-oncogene (JUNB), prostaglandin-endoperoxide synthase 2 (PTGS2), early growth response 1 (EGR1), early growth response 2 (EGR2) and nuclear receptor subfamily 4, group A, member 2 (NR4A2). Among these, FOS showed the highest node degree, which was 33 (Table IV).
miRNA prediction analysis
In total, 36 miRNAs (FDR<0.05) were screened by GSEA enrichment analysis (Table V). MIR-224, which was the most significant miRNA with an FDR of 0.013101, regulated 152 target genes and had a GTGACTT target site. MIR-7 was also screened as a significant miRNA (FDR, 0.13673); it regulated 163 target genes and had a GTCTTCC target site. Moreover, CAG GTCC was the target site of MIR-492 (FDR, 0.15391), which regulated 59 target genes. There were also several miRNAs that showed combined actions on one target site, such as MIR-34A, MIR-34C and MIR-449 (FDR, 0.017363), all of which acted on CACTGCC target sites and regulated a total of 273 target genes.
miRNA regulatory network
An miRNA regulatory network was constructed from the screened miRNAs and the DEGs that they regulated. There were 176 edges and 97 nodes in the network, and the nodes consisted of 29 DEGs and 68 miRNAs (Fig. 4). All of the DEGs were regulated by one or several miRNAs. For instance, EGR2 was regulated by 8 miRNAs, namely MIR-150, MIR-142-3P, MIR-367, MIR-363, MIR-32, MIR-92, MIR-25 and MIR-224, whereas ATF3 was only regulated by MIR-494. Similarly, each miRNA regulated one or more DEGs. MIR-142-5P regulated 6 DEGs simultaneously, namely CD69 molecule (CD69), ZFP36 ring finger protein (ZFP36), ADAM metallopeptidase with thrombospondin type 1 motif 1 (ADAMTS1), FOS, stanniocalcin 1 and early growth response 3 (EGR3), whereas MIR-194 only had a regulatory effect on solute carrier family 2 (facilitated glucose transporter), member 3 (SLC2A3).
Discussion
OPCAB is a type of surgery that improves blood flow to heart. In this study, GSE12486 was downloaded from the GEO to investigate the effect of this surgery on myocardial tissue on the basis of molecular mechanisms. A total of 64 DEGs, including 63 upregulated and 1 downregulated DEGs were screened. PLIN1 was the unique downregulated DEG in this study, which has been confirmed to be associated with anterolateral myocardial infarction through various metabolic pathways and axon guidance (23,24).
In addition, the present study also identified that the DEGs were mainly enriched in GO terms associated with chemokine activity, transcription factor activity and transcription regulator activity, and several KEGG pathways, such as chemokine signaling pathway, cytokine-cytokine receptor interaction, and Toll-like receptor signaling pathway. The FOS gene was enriched in various GO terms such as inflammatory response, defense response and transcription regulator activity. Wu et al have shown that FOS is an important target of immune inflammation responses in dendritic cells (25). The OPCAB surgery has been confirmed to activate an inflammatory response, specific hemostatic responses and immune mechanisms (26). Thus, the FOS gene may be involved in the changes that occur following OPCAB grafting via an immune mechanism. In addition, FOS, IL6 and JUN were the enriched genes in the Toll-like receptor and MAPK signaling pathways. The MAPK signaling pathway can be regulated by the activation of extracellular signal-regulated kinase (ERK)1 and 2, which can translocate to the nucleus and activate transcription factors, such as c-FOS (27). Inhibition of the activation of MAPK may attenuate warm ischemia-reperfusion injury (28), and MAPK inhibition is widely used to improve cardiac function after rewarming, cold storage and reperfusion (29). IL6 is an independent predictor for left ventricular systolic function in patients following OPCAB surgery (30). Furthermore, the ratio IL10/IL6 can be used as an inflammatory biomarker in OPCAB patients (31). In the present study, IL6 and IL8 were enriched in inflammation-related terms, such as chemotaxis, inflammatory response, defense response and response to wounding.
Moreover, a PPI network was constructed, and 10 hub nodes including FOS, JUN, ATF3, IL6, FOSB and JUNB were screened. Among them, FOS and JUN were the most important nodes with degrees of 33 and 25, respectively. A previous study has shown a direct association of signal transducer and activator of transcription 3 with c-FOS and c-JUN in response to IL6 (32). The Fos family, including c-FOS, FOSB and FOSB2, are able to form dimers with JUN proteins to regulate target genes (29). Furthermore, ATF3 has been confirmed to inhibit IL6 transcription by altering the structure of chromatin, and then restrict access to transcription factors (33). In the present study, FOS, JUN, IL6 and ATF3 were also screened and enriched in transcription factor activity. A correlation between FOS and EGR2 was also shown. EGR2, which is a human zinc finger-encoding gene, is induced with c-Fos-like kinetics by various mitogens in cells (34).
A total of 36 miRNAs were screened including MIR-494, MIR-501, MIR-524 and MIR-142-3P. Furthermore, an miRNA regulatory network was constructed based on the DEGs and screened miRNAs. It was found that ATF3 was only regulated by MIR-494. Similarly, Lan et al found that following reperfusion, the expression of ATF3 was inhibited by upregulated miRNA-494, and also associated with inflammation and adhesion (35). In addition, there was mutual regulated relationship between EGR2 and MIR-142-3P: ERG2 was encoded by the genes repressed by MIR-142-3P, while ERG2 associated with nerve growth factor-induced gene-A (NGFI-A) binding protein 2 binds to the pre-MIR-142-3p promoter to regulate its expression negatively (36).
In conclusion, the identified DEGs, such as FOS, IL6, EGR2 and JUN, might participate in biological processes including immune inflammation responses, the Toll-like receptor signaling pathway and the defense response of patients following OPCAB surgery. Moreover, the associations between DEGs and miRNA were screened. ATF3 was regulated by MIR-494, and JUN was regulated by MIR-501 and MIR494. However, these results require confirmation by experimental study.
Acknowledgements
The authors express warm thanks to Fenghe (Shanghai) Information Technology Co., Ltd., whose ideas and help provided a valuable added dimension to the study.
References
Worthley SG, Osende JI, Helft G, Badimon JJ and Fuster V: Coronary artery disease: Pathogenesis and acute coronary syndromes. Mt Sinai J Med. 68:167–181. 2001.PubMed/NCBI | |
Nalysnyk L, Fahrbach K, Reynolds MW, Zhao SZ and Ross S: Adverse events in coronary artery bypass graft (CABG) trials: A systematic review and analysis. Heart. 89:767–772. 2003. View Article : Google Scholar : PubMed/NCBI | |
Kim NY, Shim JK, Bang SO, Sim JS, Song JW and Kwak YL: Effects of ulinastatin on coagulation in high-risk patients undergoing off-pump coronary artery bypass graft surgery. Korean J Anesthesiol. 64:105–111. 2013. View Article : Google Scholar : PubMed/NCBI | |
Arab S, Konstantinov IE, Boscarino C, Cukerman E, Mori A, Li J, Liu PP, Redington AN and Coles JG: Early gene expression profiles during intraoperative myocardial ischemia-reperfusion in cardiac surgery. J Thorac Cardiovasc Surg. 134:74–81.e2. 2007. View Article : Google Scholar : PubMed/NCBI | |
Gabrielsen A, Lawler PR, Yongzhong W, Steinbrüchel D, Blagoja D, Paulsson-Berne G, Kastrup J and Hansson GK: Gene expression signals involved in ischemic injury, extracellular matrix composition and fibrosis defined by global mRNA profiling of the human left ventricular myocardium. J Mol Cell Cardiol. 42:870–883. 2007. View Article : Google Scholar : PubMed/NCBI | |
Konstantinov IE, Coles JG, Boscarino C, Takahashi M, Goncalves J, Ritter J and Van Arsdell GS: Gene expression profiles in children undergoing cardiac surgery for right heart obstructive lesions. J Thorac Cardiovasc Surg. 127:746–754. 2004. View Article : Google Scholar : PubMed/NCBI | |
Chandrasekar B, Mitchell DH, Colston JT and Freeman GL: Regulation of CCAAT/Enhancer binding protein, interleukin-6, interleukin-6 receptor, and gp130 expression during myocardial ischemia/reperfusion. Circulation. 99:427–433. 1999. View Article : Google Scholar : PubMed/NCBI | |
Tomic V, Russwurm S, Möller E, Claus RA, Blaess M, Brunkhorst F, Bruegel M, Bode K, Bloos F, Wippermann J, et al: Transcriptomic and proteomic patterns of systemic inflammation in on-pump and off-pump coronary artery bypass grafting. Circulation. 112:2912–2920. 2005.PubMed/NCBI | |
Ghorbel MT, Cherif M, Mokhtari A, Bruno VD, Caputo M and Angelini GD: Off-pump coronary artery bypass surgery is associated with fewer gene expression changes in the human myocardium in comparison with on-pump surgery. Physiol Genomics. 42:67–75. 2010. View Article : Google Scholar : PubMed/NCBI | |
Lan YF, Chen HH, Lai PF, Cheng CF, Huang YT, Lee YC, Chen TW and Lin H: MicroRNA-494 reduces ATF3 expression and promotes AKI. J Am Soc Nephrol. 23:2012–2023. 2012. View Article : Google Scholar : PubMed/NCBI | |
Wang X, Zhang X, Ren XP, Chen J, Liu H, Yang J, Medvedovic M, Hu Z and Fan GC: MicroRNA-494 targeting both proapoptotic and antiapoptotic proteins protects against ischemia/reperfusion-induced cardiac injury. Circulation. 122:1308–1318. 2010. View Article : Google Scholar : PubMed/NCBI | |
Wang Z, Shao J, Zhou Q, Liu J, Zhu Y, Yang J and Wei M: The −251A>T polymorphism of interleukin-8 is associated with longer mechanical ventilation and hospital staying after coronary surgery. Cytokine. 50:268–272. 2010. View Article : Google Scholar : PubMed/NCBI | |
Nathan N, Preux PM, Feiss P and Denizot Y: Plasma interleukin-4, interleukin-10, and interleukin-13 concentrations and complications after coronary artery bypass graft surgery. J Cardiothorac Vasc Anesth. 14:156–160. 2000. View Article : Google Scholar : PubMed/NCBI | |
Fromes Y, Daghildjian K, Caumartin L, Fischer M, Rouquette I, Deleuze P and Bical OM: A comparison of low vs conventional-dose heparin for minimal cardiopulmonary bypass in coronary artery bypass grafting surgery. Anaesthesia. 66:488–492. 2011. View Article : Google Scholar : PubMed/NCBI | |
Gautier L, Cope L, Bolstad BM and Irizarry RA: Affy-analysis of Affymetrix GeneChip data at the probe level. Bioinformatics. 20:307–315. 2004. View Article : Google Scholar : PubMed/NCBI | |
Diboun I, Wernisch L, Orengo CA and Koltzenburg M: Microarray analysis after RNA amplification can detect pronounced differences in gene expression using limma. BMC Genomics. 7:2522006. View Article : Google Scholar : PubMed/NCBI | |
Sartor MA, Tomlinson CR, Wesselkamper SC, Sivaganesan S, Leikauf GD and Medvedovic M: Intensity-based hierarchical Bayes method improves testing for differentially expressed genes in microarray experiments. BMC Bioinformatics. 7:5382006. View Article : Google Scholar : PubMed/NCBI | |
Jolliffe I: Principal component analysis. Encyclopedia of Statistics in Behavioral Science. Everitt BS and Howell D: (New York, NY). John Wiley & Sons. 2005. View Article : Google Scholar | |
da Huang W, Sherman BT and Lempicki RA: Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 4:44–57. 2009. View Article : Google Scholar : PubMed/NCBI | |
Franceschini A, Szklarczyk D, Frankild S, Kuhn M, Simonovic M, Roth A, Lin J, Minguez P, Bork P, von Mering C and Jensen LJ: STRING v9.1: Protein-protein interaction networks, with increased coverage and integration. Nucleic Acids Res. 41:D808–D815. 2013. View Article : Google Scholar : PubMed/NCBI | |
Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES and Mesirov JP: Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA. 102:15545–15550. 2005. View Article : Google Scholar : PubMed/NCBI | |
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B and Ideker T: Cytoscape: A software environment for integrated models of biomolecular interaction networks. Genome Res. 13:2498–2504. 2003. View Article : Google Scholar : PubMed/NCBI | |
Jaffer I, Riederer M, Shah P, Peters P, Quehenberger F, Wood A, Scharnagl H, März W, Kostner KM and Kostner GM: Expression of fat mobilizing genes in human epicardial adipose tissue. Atherosclerosis. 220:122–127. 2012. View Article : Google Scholar : PubMed/NCBI | |
Lei P, Baysa A, Nebb HI, Valen G, Skomedal T, Osnes JB, Yang Z and Haugen F: Activation of Liver X receptors in the heart leads to accumulation of intracellular lipids and attenuation of ischemia-reperfusion injury. Basic Res Cardiol. 108:3232013. View Article : Google Scholar : PubMed/NCBI | |
Wu C, Gong Y, Yuan J, Zhang W, Zhao G, Li H, Sun A, Hu Kai Zou Y and Ge J: microRNA-181a represses ox-LDL-stimulated inflammatory response in dendritic cell by targeting c-Fos. J Lipid Res. 53:2355–2363. 2012. View Article : Google Scholar : PubMed/NCBI | |
Levy JH and Tanaka KA: Inflammatory response to cardiopulmonary bypass. Ann Thorac Surg. 75(Suppl): S715–S720. 2003. View Article : Google Scholar : PubMed/NCBI | |
Junttila MR, Li SP and Westermarck J: Phosphatase-mediated crosstalk between MAPK signaling pathways in the regulation of cell survival. FASEB J. 22:954–965. 2008. View Article : Google Scholar : PubMed/NCBI | |
Kobayashi M, Takeyoshi I, Yoshinari D, Matsumoto K and Morishita Y: P38 mitogen-activated protein kinase inhibition attenuates ischemia-reperfusion injury of the rat liver. Surgery. 131:344–349. 2002. View Article : Google Scholar : PubMed/NCBI | |
Vassalli G, Milano G and Moccetti T: Role of mitogen-activated protein kinases in myocardial ischemia-reperfusion injury during heart transplantation. J Transplant. 2012:9289542012. View Article : Google Scholar : PubMed/NCBI | |
Mathew GM, Mathew DC, Lo SC, Alexios GM, Yang JC, Sashikumar JM, Shaikh TM and Huang CC: Synergistic collaboration of gut symbionts in Odontotermes formosanus for lignocellulosic degradation and bio-hydrogen production. Bioresour Technol. 145:337–344. 2013. View Article : Google Scholar : PubMed/NCBI | |
Mazandarani M, Mojtahedzadeh M, Yousefshahi F and Hamishehkar H: Effect of pre-treatment with hypertonic saline on IL-10/IL-6 ratio in CABG patient. Res Pharm Sci. 7:S8862012. | |
Schuringa JJ, Timmer H, Luttickhuizen D, Vellenga E and Kruijer W: c-Jun and c-Fos cooperate with STAT3 in IL-6-induced transactivation of the IL-6 response element (IRE). Cytokine. 14:78–87. 2001. View Article : Google Scholar : PubMed/NCBI | |
Gilchrist M, Thorsson V, Li B, Rust AG, Korb M, Roach JC, Kennedy K, Hai T, Bolouri H and Aderem A: Systems biology approaches identify ATF3 as a negative regulator of Toll-like receptor 4. Nature. 441:173–178. 2006. View Article : Google Scholar : PubMed/NCBI | |
Rangnekar VM, Aplin AC and Sukhatme VP: The serum and TPA responsive promoter and intron-exon structure of EGR2, a human early growth response gene encoding a zinc finger protein. Nucleic Acids Res. 18:2749–2757. 1990. View Article : Google Scholar : PubMed/NCBI | |
Lan HY and Chung AC: TGF-β/Smad signaling in kidney disease. Semin Nephrol. 32:236–243. 2012. View Article : Google Scholar : PubMed/NCBI | |
Lagrange B, Martin RZ, Droin N, Aucagne R, Paggetti J, Largeot A, Itzykson R, Solary E, Delva L and Bastie JN: A role for miR-142-3p in colony-stimulating factor 1-induced monocyte differentiation into macrophages. Biochim Biophys Acta. 1833:1936–1946. 2013. View Article : Google Scholar : PubMed/NCBI |