BioMed Research International
Volume 2020 (2020), Article ID 8062902, 8 pages
https://doi.org/10.1155/2020/8062902

Research Article

Comprehensive Study of Different Expressed Genes and Their Functional Modules in Anesthesia for Off-Pump Coronary Artery Bypass Grafting

Hui Zhou,¹ Wang Min,² and Zhihua Zhu¹

¹Department of Anesthesia, China-Japan Union Hospital of Jilin University, Changchun, Jilin Province, China
²Department of Pathology, Jilin Provincial Cancer Hospital, Changchun, Jilin Province, China

Correspondence should be addressed to Zhihua Zhu; zhuizi732025492@163.com

Received 9 May 2020; Accepted 1 June 2020; Published 6 July 2020

Guest Editor: Tao Huang

Copyright © 2020 Hui Zhou et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Purpose. The effect of preoperative anesthesia on coronary artery bypass grafting without extracorporeal circulation is not apparent. We want to investigate the effects and molecular mechanisms of two anesthesia methods on the treatment of coronary artery bypass grafting (OPCABG) under extracorporeal circulation. Patients and Methods. The data of inhaled anesthesia and intravenous anesthesia before coronary artery bypass grafting were downloaded from the GEO database, and the differences were analyzed with the control group. The combination of multiple analytical methods can decipher the mechanism of anesthesia on surgery, including protein interaction network analysis, enrichment analysis, and regulatory subprediction. Results. This study obtained 6699 differential genes under two kinds of anesthesia before OPCABG. By constructing a protein interaction network of differentially expressed genes, we obtained 14 functional module networks. By predicting regulators of functional module genes, we revealed a series of ncRNAs (miR-129-5p, miR-340-5p, and miR-410-3p) and transcription factors (VHL and YBX1). Conclusion. Based on functional module network analysis, we identified the effects of preoperative inhalation anesthesia and intravenous anesthesia on OPCABG, which provides a valuable theoretical reference for subsequent clinical studies.

1. Introduction

Off-pump coronary artery bypass grafting (OPCABG) is the latest surgical procedure in cardiac surgery [1]. OPCABG has proven to be an effective surgical revascularization program [2]. Surgery has determined the safety of OPCABG and its short-term efficacy, but long-term clinical outcomes are uncertain. OPCABG is still the treatment of choice in modern cardiac surgery, considering that surgical practice is continuously changing [3]. Compared with cardiopulmonary bypass grafting, OPCABG has particular advantages in reducing postoperative complications, including systemic inflammatory response, myocardial injury, renal injury, and brain injury [4]. In theory, OPCABG can improve long-term survival by reducing perioperative complications, such as stroke, myocardial damage, cardiac-related mortality, and neurocognitive impairment [5].

Most importantly, OPCABG reduces the risk of neurological complications compared with extracorporeal coronary artery bypass grafting [6]. OPCABG, a high-tech myocardial revascularization program, can reduce sympathetic stress and improve hemodynamic changes [7]. OPCABG has become the standard surgery in Japan [8]. The transplantation of OPCABG in high-risk patients remains controversial, but studies have shown its potential benefits. Therefore, in high-risk patients, the application of its technology is still valued and studied [9]. For patients undergoing OPCABG surgery, the measurement of anesthesia depth is clinically significant and can achieve the effects of avoiding intraoperative perception and cardiac suppression. The recent introduction of entropy as a monitor of anesthesia depth determines the amount of anesthetic used in patients with OPCABG surgery [10]. Given the high hemodynamic instability and the risk of organ damage in OPCABG surgery, xenon anesthesia has become an attractive aesthetic with good hemodynamic and organ protection properties [11]. Studies have shown that during OPCABG, compared with total intravenous anesthesia, volatile anesthetics can reduce myocardial damage, which is measured via cardiac troponin levels [12]. Compared with the sputum, the anesthetic function of sevoflurane in OPCABG reduces the risk of OPCABG surgery due to vasopressin [13]. Changes in magnesium levels can cause abnormal blood coagulation, leading to bleeding complications after OPCABG [14]. In OPCABG, ropivacaine and fentanyl can be used to temporarily reduce arterial pressure, optimize myocardial function, and affect perioperative fluid and vasoactive treatment [15]. Recent studies have shown that dexmedetomidine, as an effective adjuvant, can reduce hospital stays in patients receiving OPCABG [16]. Different expressed genes and their functional modules in anesthesia for OPCABG need to be further explored. Based on a functional modular network, we analyze the effects of preoperative inhalation anesthesia and intravenous anesthesia on OPCABG to explore the underlying molecular mechanisms. It provides not only valuable research directions but also abundant valuable theoretical references for further research on the clinical application of anesthesia.

2. Material and Methods

2.1. Data Resource

The STRING database [17] is a search tool specifically designed to retrieve protein-protein interactions, which provides comprehensive insight into the currently available PPIs and can, therefore, be used for a wide range of PPI analyses. In this study, the STRING database was utilized for the access of all human protein interaction data, and we set the , involving 405,916 interaction pairs of 10,514 proteins.

2.2. Differential Expression Analysis

We collected an expression microarray data set for inhaled anesthesia and intravenous anesthesia before coronary artery bypass grafting from the NCBI Gene Expression Omnibus database [18], numbered GSE4386. In this data matrix, patients scheduled for off-pump CABG were randomized into a group with the anesthetic gas sevoflurane () or the intravenous anesthetic propofol (). Atrial samples were collected at the beginning and end of bypass surgery to determine gene expression profiles. A two-difference analysis (anaerobic gas sevoflurane-intravenous anesthesia with propofol) was performed on the collected samples and calculated using “limma” package in R [19]. The 6699 differential genes were combined to construct an expression profile matrix for nonexternal coronary artery bypass grafting.

2.3. Functional and Pathway Enrichment Analysis

We explored the function of differential genes and their involved signaling pathways. We performed a GO function (with a cutoff value of 0.01 and a cutoff value of 0.01) and a KEGG pathway enrichment analysis (with a cutoff value of 0.05 and a cutoff value of 0.2) for the 14 modules of the gene using clusterProfiler [20]. We screened the process-related functions and pathways during anesthesia for OPCABG and mapped the bubbles.

2.4. Identification of Transcription Factors (TFs) and Modules of ncRNA Regulation

We defined the pivot regulator as a kind of modulator that significantly affects the functional modular network of off-peak coronary artery bypass grafting. Noncoding genes and TFs often drive the transcription, as well as the posttranscriptional regulation of the genes. We downloaded the relevant TF target data from the TRRUST database [21] and finally obtained 33 interaction pairs of 25 TFs. Afterwards, ncRNA-mRNA data was accessed from the RAID database [22], and 660 interaction pairs including 486 ncRNAs were also obtained. Pivot analysis was then performed according to these interactions to further identify the regulatory effects of these TFs and ncRNAs in the module. The saliency of the interactions among the drivers and the modules is calculated based on the hypergeometric test. TF and ncRNA are pivoted for the essential regulatory module, with a value < 0.01 as the screening standard. We performed a statistical analysis of the pivot, and the central function of the more dysfunctional module was identified as the core pivot.

3. Results

3.1. To Determine the Effect of Preoperative Anesthesia on Gene Expression during an OPCABG

We performed a differential analysis of the gene expression profiles of OPCABG under two kinds of anesthesia to further understand the effect of anesthesia on OPCABG. The 5701 differential genes in the CABG state of sevoflurane anesthesia were obtained, and 3210 differential genes were anesthetized by propofol in the CABG state. We combined the two sets of differences to obtain 6699 differential genes, which were considered to be critical genes for the effects of anesthesia on OPCABG.

3.2. Identify Functional Module Networks

Further, we explained the effects of two anesthesia methods on OPCABG from the gene network level. Based on PPI analysis, a protein interaction network was constructed using 6699 differential genes from OPCABG. As a result, a total of 14 interactive networks were obtained, called functional modules (Figure 1). According to functional modules, we identify the critical genes for each module. Analysis of its essential genes may be involved in different functions and pathways, causing different effects of OPCABG in both anesthesia situations.

Figure 1: Based on the human interaction protein network, we obtained a protein interaction network map of 6699 differential genes clustered into 14 functional module networks. This network may represent the mechanism of anesthesia affecting nonexternal coronary artery bypass grafting.

3.3. Functions and Pathways Involved in Specific Genes

To further understand the biological aspects of the module genes, GO function and KEGG analyses were performed on 14 functional module networks. We obtained GO terminology, including 1314 cell composition entries, 1640 molecular functional terms, and 13,160 biological processes (Figure 2(a)). Most of the modular genes are associated with the regulation of the cytokine-mediated signaling pathway and positive regulation of lymphocyte differentiation. The results of KEGG pathway enrichment reveal that the functional module genes are mainly responsible for platelet activation, cytokine-cytokine receptor interaction, and apoptosis-multiple species (Figure 2(b)). The above results of modular gene enrichment are firmly related to OPCABG; we therefore identified the 14 dysfunction modules.

(a)

(b)

Figure 2: Functional pathways included in modular gene identification of dysfunction modules for osteoarthritis. (a) GO analysis for module genes; (b) KEGG pathway enrichment analysis. The color depth represents the degree of enrichment, and the size of the circle represents the proportion of the genes in the module. FDR: false discovery rate.

3.4. Regulation of TF and ncRNA of Genes Associated with OPCABG under Anesthesia

To explore the regulatory factors, we applied regulator-predictive analysis to functional modules according to the regulatory relationships in transcriptional and posttranscriptional processes. Based on the number of regulatory modules and the significance of the values, we obtained 486 ncRNAs involving 660 ncRNA-module target pairs and 25 TFs involving 33 TF-module target pairs. Also, we analyzed the predicted results and found that miR-129-5p targets up to six dysfunction modules with a significant impact on OPCABG (Figure 3). The miR-129-5p and miR-410-3p also regulate five dysfunction modules, having remarkable regulatory effects on the module network. VHL and YBX1 have significant regulatory effects on three dysfunction modules and play a key role in the impact of anesthesia on OPCABG (Figure 4). Other TFs also have their specific regulatory effects on the functional module genes, which may have a significant impact on surgery.

Figure 3: Noncoding RNA (ncRNA) regulatory network of nonexternal coronary artery bypass grafting. The blue triangles represent the module while the pink ovals represent the ncRNA corresponding with the module.

Figure 4: Regulatory network map constructed by transcription factors (TFs) for nonexternal coronary artery bypass grafting. The cyan rectangles represent the module while the green ovals represent the TF corresponding with the module.

4. Discussion

In the study of OPCABG, anesthesia can help to reduce intraoperative blood loss and intraoperative blood transfusion demands to varying degrees, shortening hospital stay. The most important is to reduce myocardial enzyme leakage, inflammation, and kidney and nerve damage [23]. However, the use of different anesthesia in OPCABG has different effects on postoperative outcomes. For example, the use of volatile anesthetics to maintain anesthesia in OPCABG surgery contributes to adequate anesthesia depth and reduces the need for analgesia [24]. Helium shows less hemodynamic instability in OPCABG than in conventionally used anesthetics, and it reduces the need for vasopressors [25]. For patients receiving OPCABG, it is important for them to reduce myocardial damage during surgery. The volatility induction and maintenance of sevoflurane has a certain degree of improvement in cardiac function in patients receiving OPCABG. Therefore, compared with intravenous propofol, sevoflurane has a more protective effect on heart damage [26]. The volatile drugs sevoflurane and desflurane can be considered to have a more protective effect on the myocardium of the whole vein [27]. These studies have validated the important role of inhaled anesthesia and intravenous anesthesia in the process of OPCABG, which prompted us to explore the effects of anesthesia in OPCABG. We constructed a PPI network according to differentially expressed RNA data sets for OPCABG after anesthesia and explored 14 functional modules. We found that these modules are mainly involved in positive regulation of lymphocyte differentiation and cytokine-cytokine receptor interaction. The Holler E. study found that after OPCABG, circulating endothelial cell number and apoptotic endothelial cell death are markers of endothelial cell activation and injury. Studies have shown that in OPCABG surgery, in order to protect circulating lymphocytes, the combination of intravenous propofol and gas anesthesia with sevoflurane is superior to the use of sevoflurane to maintain anesthesia [28]. One of the potential advantages of OPCABG is to attenuate systemic inflammatory responses caused by neutrophil activation [29]. Moreover, OPCABG is beneficial to the secretion of endogenous erythropoietin in patients, reducing the blood loss of surgery [30]. Studies have found that OPCABG is associated with decreased troponin I levels and activation, while elevated levels of troponin I and proinflammatory cytokines are present in most cardiac surgery patients [31]. To elucidate the transcriptional regulatory factors of OPCABG by anesthetic gases and intravenous anesthesia, we performed an analysis according to the regulatory relationships both in transcription and posttranscription. It was found that microRNAs (miR-129-5p, miR-340-5p, and miR-410-3p) and transcription factors (VHL and YBX1) have remarkable regulatory effects on functional modules. MicroRNAs play a key role in a variety of cellular processes [32]. miR-129-5p acts on neurons, which have an inhibitory effect on blocking the synaptic reduction in vitro and reducing the severity of seizures in vivo [33]. This also indicates that anesthesia has an important influence on postoperative cell regulation. Abnormal expression of long noncoding RNA plays a crucial part in regulating the progression and drug resistance of human tumors [34]. LGR5 expression was downregulated by the Wnt/β-catenin pathway. It leads to overexpression of miR-340-5p, which inhibits apoptosis, cell proliferation, and its resistance [35]. This validated the regulatory changes in genes affected by gene expression after OPCABG. Also, the transcription factor VHL plays a key role in the inactivation of the hypoxia-inducible factor (HIF-1) through regulating the PI3K/AKT/mTOR pathway activity [36]. YBX1 is a multifunctional protein [37]. It is involved in all kinds of DNA/RNA-dependent events, such as DNA repair, mRNA packaging, mRNA transcription and splicing, translational regulation, and mRNA stability. At the cellular level, various activities of YBX1 appear to be involved in stress response, cell differentiation and proliferation, and malignant cell transformation [38]. By regulating cell differentiation, cytokine synthesis, and chemokine synthesis, monocyte YBX1 has a prominent and unique role in cell feed-forward crosstalk and regression of inflammatory processes [39]. Based on the understanding of the anesthetic gas and intravenous anesthesia for the pivot regulator of OPCABG, we can identify the effects of these significant regulators on postoperative outcomes. Different functions are performed in the development of OPCABG after anesthesia.

5. Conclusion

A comprehensive functional module network was obtained. The modules, which offer some proven genes and transcription factors to be tested for OPCABG after anesthesia, provides a theoretical basis and reference for further research.

Data Availability

The data associated with our study was downloaded from GSE4386.

Conflicts of Interest

There are no conflicts of interest.

References

K. K. Saha, “Off-pump coronary artery bypass grafting in India,” Indian Heart Journal, vol. 66, no. 2, pp. 203–207, 2014.

B. Keeling, V. Thourani, G. Aliawadi et al., “Conversion from off-pump coronary artery bypass grafting to on-pump coronary artery bypass grafting,” The Annals of Thoracic Surgery, vol. 104, no. 4, pp. 1267–1274, 2017.

R. Cartier, “Current trends and technique in OPCAB surgery,” Journal of Cardiac Surgery, vol. 18, no. 1, pp. 32–46, 2003.

K. Nakazato and A. Sakamoto, “OPCAB,” Masui, vol. 63, no. 5, pp. 506–512, 2014.

S. Shaefi, A. Mittel, D. Loberman, and H. Ramakrishna, “Off-pump versus on-pump coronary artery bypass grafting-a systematic review and analysis of clinical outcomes,” Journal of Cardiothoracic and Vascular Anesthesia, vol. 33, no. 1, pp. 232–244, 2019.

W. Pawliszak, M. Kowalewski, G. M. Raffa et al., “Cerebrovascular events after no-touch off-pump coronary artery bypass grafting, conventional side-clamp off-pump coronary artery bypass, and proximal anastomotic devices: a meta-analysis,” Journal of the American Heart Association, vol. 5, no. 2, 2016.

C. Rajakaruna, C. Rogers, K. Pike et al., “Superior haemodynamic stability during off-pump coronary surgery with thoracic epidural anaesthesia: results from a prospective randomized controlled trial,” Interactive Cardiovascular and Thoracic Surgery, vol. 16, no. 5, pp. 602–607, 2013.

H. Yaku, K. Doi, and K. Okawa, “Off-pump coronary artery bypass grafting revisited: experience and evidence from Japan,” Annals of Thoracic and Cardiovascular Surgery, vol. 19, no. 2, pp. 83–94, 2013.

G. A. Guida, P. Chivasso, D. Fudulu et al., “Off-pump coronary artery bypass grafting in high-risk patients: a review,” Journal of Thoracic Disease, vol. 8, Supplement 10, pp. S795–S798, 2016.

M. Jiahai, W. Xueyan, X. Yonggang et al., “Spectral entropy monitoring reduces anesthetic dosage for patients undergoing off-pump coronary artery bypass graft surgery,” Journal of Cardiothoracic and Vascular Anesthesia, vol. 26, no. 5, pp. 818–821, 2012.

L. Al tmimi, J. Van Hemelrijck, M. Van de Velde et al., “Xenon anaesthesia for patients undergoing off-pump coronary artery bypass graft surgery: a prospective randomized controlled pilot trial†,” British Journal of Anaesthesia, vol. 115, no. 4, pp. 550–559, 2015.

F. Guarracino, G. Landoni, L. Tritapepe et al., “Myocardial damage prevented by volatile anesthetics: a multicenter randomized controlled study,” Journal of Cardiothoracic and Vascular Anesthesia, vol. 20, no. 4, pp. 477–483, 2006.

T. M. Hemmerling, G. Romano, N. Terrasini, and N. Noiseux, “Anesthesia for off-pump coronary artery bypass surgery,” Annals of Cardiac Anaesthesia, vol. 16, no. 1, pp. 28–39, 2013.

F. Sabzi and R. Faraji, “The effect of magnesium sulfate on post off-pump coronary artery bypass grafting bleeding,” Indian Heart Journal, vol. 68, no. 3, pp. 349–354, 2016.

M. Y. Kirov, A. V. Eremeev, A. A. Smetkin, and L. J. Bjertnaes, “Epidural anesthesia and postoperative analgesia with ropivacaine and fentanyl in off-pump coronary artery bypass grafting: a randomized, controlled study,” BMC Anesthesiology, vol. 11, no. 1, p. 17, 2011.

M. A. Khalil and M. S. Abdel Azeem, “The impact of dexmedetomidine infusion in sparing morphine consumption in off-pump coronary artery bypass grafting,” Seminars in Cardiothoracic and Vascular Anesthesia, vol. 17, no. 1, pp. 66–71, 2013.

D. Szklarczyk, J. H. Morris, H. Cook et al., “The STRING database in 2017: quality-controlled protein-protein association networks, made broadly accessible,” Nucleic Acids Research, vol. 45, no. D1, pp. D362–D368, 2017.

T. Barrett, S. E. Wilhite, P. Ledoux et al., “NCBI GEO: archive for functional genomics data sets--update,” Nucleic Acids Research, vol. 41, no. Database issue, pp. D991–D995, 2013.

M. E. Ritchie, B. Phipson, D. Wu et al., “limma powers differential expression analyses for RNA-sequencing and microarray studies,” Nucleic Acids Research, vol. 43, no. 7, p. e47, 2015.

G. Yu, L. G. Wang, Y. Han, and Q. Y. He, “clusterProfiler: an R package for comparing biological themes among gene clusters,” OMICS, vol. 16, no. 5, pp. 284–287, 2012.

H. Han, J. W. Cho, S. Lee et al., “TRRUST v2: an expanded reference database of human and mouse transcriptional regulatory interactions,” Nucleic Acids Research, vol. 46, no. D1, pp. D380–D386, 2018.

Y. Yi, Y. Zhao, C. Li et al., “RAID v2.0: an updated resource of RNA-associated interactions across organisms,” Nucleic Acids Research, vol. 45, no. D1, pp. D115–D118, 2017.

J. M. Brown, R. S. Poston, J. S. Gammie et al., “Off-pump versus on-pump coronary artery bypass grafting in consecutive patients: decision-making algorithm and outcomes,” The Annals of Thoracic Surgery, vol. 81, no. 2, pp. 555–561, 2006, discussion 561.

B. G. Venkatesh, Y. Mehta, A. Kumar, and N. Trehan, “Comparison of sevoflurane and isoflurane in OPCAB surgery,” Annals of Cardiac Anaesthesia, vol. 10, no. 1, pp. 46–50, 2007.

L. Al tmimi, S. Devroe, G. Dewinter et al., “Xenon as an adjuvant to propofol anesthesia in patients undergoing off-pump coronary artery bypass graft Surgery,” Anesthesia and Analgesia, vol. 125, no. 4, pp. 1118–1128, 2017.

X. Liu, X. Liu, R. Wang et al., “Circulating microRNAs indicate cardioprotection by sevoflurane inhalation in patients undergoing off-pump coronary artery bypass surgery,” Experimental and Therapeutic Medicine, vol. 11, no. 6, pp. 2270–2276, 2016.

S. Suryaprakash, M. Chakravarthy, G. Muniraju et al., “Myocardial protection during off pump coronary artery bypass surgery: a comparison of inhalational anesthesia with sevoflurane or desflurane and total intravenous anesthesia,” Annals of Cardiac Anaesthesia, vol. 16, no. 1, pp. 4–8, 2013.

L. Jia, R. Dong, F. Zhang et al., “Propofol provides more effective protection for circulating lymphocytes than sevoflurane in patients undergoing off-pump coronary artery bypass graft surgery,” Journal of Cardiothoracic and Vascular Anesthesia, vol. 29, no. 5, pp. 1172–1179, 2015.

S. Al-Ruzzeh, G. Hoare, N. Marczin et al., “Off-pump coronary artery bypass surgery is associated with reduced neutrophil activation as measured by the expression of CD11b: a prospective randomized study,” The Heart Surgery Forum, vol. 6, no. 2, pp. 89–93, 2003.

A. Drohomirecka, M. Kuśmierczyk, S. Kocańda et al., “Endogenic erythropoietin secretion in patients undergoing off-pump coronary artery bypass grafting,” Kardiologia Polska, vol. 75, no. 5, pp. 470–475, 2016.

N. Nesher, I. Frolkis, M. Vardi et al., “Higher levels of serum cytokines and myocardial tissue markers during on-pump versus off-pump coronary artery bypass surgery,” Journal of Cardiac Surgery, vol. 21, no. 4, pp. 395–402, 2006.

F. Tang, M. Kaneda, D. O'Carroll et al., “Maternal microRNAs are essential for mouse zygotic development,” Genes & Development, vol. 21, no. 6, pp. 644–648, 2007.

M. Rajman, F. Metge, R. Fiore et al., “A microRNA-129-5p/Rbfox crosstalk coordinates homeostatic downscaling of excitatory synapses,” The EMBO Journal, vol. 36, no. 12, pp. 1770–1787, 2017.

L. Song, Z. Zhou, Y. Gan et al., “Long noncoding RNA OIP5-AS1 causes cisplatin resistance in osteosarcoma through inducing the LPAATβ/PI3K/AKT/mTOR signaling pathway by sponging the miR-340-5p,” Journal of Cellular Biochemistry, vol. 120, no. 6, pp. 9656–9666, 2018.

S. Shi, X. Chen, H. Liu et al., “LGR5 acts as a target of miR-340-5p in the suppression of cell progression and drug resistance in breast cancer via Wnt/β-catenin pathway,” Gene, vol. 683, pp. 47–53, 2019.

M. Cuperlovic-Culf, K. Cormier, M. Touaibia et al., “(1)H NMR metabolomics analysis of renal cell carcinoma cells: effect of VHL inactivation on metabolism,” International Journal of Cancer, vol. 138, no. 10, pp. 2439–2449, 2016.

L. L. Zhou, J. Ni, W. T. Feng et al., “High YBX1 expression indicates poor prognosis and promotes cell migration and invasion in nasopharyngeal carcinoma,” Experimental Cell Research, vol. 361, no. 1, pp. 126–134, 2017.

D. N. Lyabin, I. A. Eliseeva, and L. P. Ovchinnikov, “YB-1 protein: functions and regulation,” Wiley Interdisciplinary Reviews: RNA, vol. 5, no. 1, pp. 95–110, 2014.

A. Bernhardt, A. Fehr, S. Brandt et al., “Inflammatory cell infiltration and resolution of kidney inflammation is orchestrated by the cold-shock protein Y-box binding protein-1,” Kidney International, vol. 92, no. 5, pp. 1157–1177, 2017.