
==== Front
Dis Markers
Dis Markers
DM
Disease Markers
0278-0240
1875-8630
Hindawi

10.1155/2023/5493415
Research Article
Identification of Inflammatory Gene in the Congenital Heart Surgery Patients following Cardiopulmonary Bypass via the Way of WGCNA and Machine Learning Algorithms
Cai Liang
https://orcid.org/0000-0002-1817-4346
Zhang Bingdong zbd2019@126.com

Department of Anesthesiology in Cardiovascular Institute, The First Affiliated Hospital of Guangxi Medical University, Nanning, Guangxi, China
Academic Editor: Jian Song

2023
12 4 2023
2023 549341527 9 2022
7 11 2022
24 11 2022
Copyright © 2023 Liang Cai and Bingdong Zhang.
2023
https://creativecommons.org/licenses/by/4.0/ 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.
Performing cardiopulmonary bypass (CPB) to reduce ischemic injury during surgery is a common approach to cardiac surgery. However, this procedure can lead to systemic inflammation and multiorgan dysfunction. Therefore, elucidating the molecular mechanisms of CPB-induced inflammatory cytokine release is essential as a critical first step in identifying new targets for therapeutic intervention. The GSE143780 dataset which is mRNA sequencing from total circulating leukocytes of the neonatorum was downloaded from the Gene Expression Omnibus (GEO) database. A total of 21 key module genes were obtained by analyzing the intersection of differentially expressed gene (DEG) and gene coexpression network analysis (WGCNA), and then, 4 genes (TRAF3IP2-AS1, PPARGC1B, CD4, and PDLIM5) were further confirmed after the least absolute shrinkage and selection operator (LASSO) and support vector machine recursive feature elimination (SVM-RFE) screening and were used as hub genes for CPB-induced inflammatory cytokine release in patients with congenital heart defects. The enrichment analysis revealed 21 key module genes mainly related to the functions of developmental cell growth, regulation of monocyte differentiation, regulation of myeloid leukocyte differentiation, ERK1 and ERK2 cascade, volume-sensitive anion channel activity, and estrogen receptor binding. The result of gene set enrichment analysis (GSEA) showed that the hub genes were related to different physiological functions of cells. The ceRNA network established for hub genes includes 3 hub genes (PPARGC1B, CD4, and PDLIM5), 55 lncRNAs, and 34 miRNAs. In addition, 4 hub genes have 215 potential therapeutic agents. Finally, expression validation of the four hub genes revealed that they were all significantly low expressed in the surgical samples than before.
==== Body
pmc1. Introduction

Cardiopulmonary bypass (CPB) is a technique for maintaining survival, which is vital to the success of modern cardiac surgeries. The main reason for application is most sufferers undergoing cardiac surgery need CPB to reduce complications caused by the surgery itself. The blood will be exposed to high shear stress and plastic tubes during blood drawing by the artificial cardiopulmonary machine [1]. During the surgery, the temperature of patient is decreased to reduce metabolic needs, and the temperature rises rapidly after the procedure, which could induce a variety of complications. However, the current treatment methods, such as corticosteroids, CPB circuit coating, and modified ultrafiltration, fail to significantly reduce the complications secondary to CPB [2–5].

The way of CPB-induced inflammation involves a variety of mechanisms. Firstly, the coagulation cascade and complement system will be activated when blood flows through the foreign body surface of the extracorporeal circulation circuit, which causes the release of proinflammatory factors [6]. Secondly, aortic occlusion leads to ischemia-reperfusion injury in multiple organs, which makes these tissues and organs produce a series of inflammatory factors. Thirdly, CPB results in endotoxemia, because CPB can lead to intestinal injury and bacterial translocation. Fourthly, gaseous microemboli in CPB lead to local thrombosis and activation of complement systems, such as C3a and C5a. Finally, iNOS or ROS produced during CPB, calcium metabolism disorder, and endothelial injury will further expand the inflammatory response [7–9] .

Systemic inflammation and multiorgan dysfunction are the most common and severe CPB-induced complications, especially in neonatal and pediatric patients. Cell infiltration, vasoconstriction, and glomerular fibrin deposition caused by inflammatory cytokines lead to renal injury [8]. Sasser et al. also found that inflammatory cytokines can increase vascular permeability, which will lead to multiple organ edema. For example, pulmonary edema can affect gas exchange [10]. The mechanism of CPB-induced inflammation is extremely complex, and there is no consensus on how to effectively reduce the inflammation. Considering inhibition of inflammation is expected to reduce incidence rate and mortality after surgery, it is important to clarify the pathological mechanism of inflammation during CPB.

In this study, we searched for the gene and signaling pathway between CPB and inflammation through comprehensive bioinformatics analysis [11]. Initially, we got the GSE143780 dataset from the Gene Expression Omnibus (GEO) database and performed gene coexpression network analysis (WGCNA) and differential analysis. Then, the least absolute shrinkage and selection operator (LASSO) and support vector machine recursive feature elimination (SVM-RFE) screening methods were performed on key module genes, and the results obtained were used as hub genes for CPB-induced inflammatory cytokine release in patients with congenital heart defects. The possible biological functions and involved signaling pathways were screened via immune and enrichment analysis of hub genes. The potential regulatory relationships of hub genes were explored by constructing endogenous RNA (ceRNA) networks. Ultimately, the expression levels of hub genes in tissues were examined by applying quantitative reverse transcription polymerase chain reaction (qRT-PCR).

2. Materials and Methods

2.1. Data Source

The GSE143780 dataset was obtained from the GEO database, including 5 CPB samples of 0 hour as control and 10 CPB-related samples of 1 hour or end as case [12].

2.2. WGCNA

The coexpression network was established by WGCNA [13]. First of all, we cluster samples according to the presence of obvious outliers. Second, automatic network construction function was used to construct the coexpression network. pickSoftThreshold in R was applied to compute the soft thresholding power, to which the proximity of coexpression was raised to compute adjacency. Then, detect modules by hierarchical clustering and the dynamic tree cut function. Finally, modules will be related to clinical traits according to gene significance (GS) and module membership (MM). The information on associated module genes was used for further analysis.

2.3. Differential Expression Analysis

So as to find out the DEGs between control and case samples and to do more in-depth functional mining, the limma package was used to identify genes with P < 0.05 and |log2 FC| > 1 as differentially expressed genes [14]. Then, the DEGs were intersected with the key module genes obtained by WGCNA to gain the key genes.

2.4. Functional Enrichment Analysis

For exploring the signaling pathways and characteristic biological attributes involved in key genes, Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) were analyzed by the function of clusterProfiler package in R. Significance was set at P < 0.05. In addition, to further search the potential biological functions of the signaling pathways involved and the selected hub genes, we carried out GSEA [15].

2.5. Ingenuity Pathway Analysis (IPA)

To explore the interaction of key genes and the diseases or functional pathways involved, we performed IPA analysis on all genes [16]. Briefly, the expression matrix was entered into the software and screened for differential genes according to a threshold (|log2 (fold change)| > 1 and P < 0.05), after which the canonical pathways of the genes and the interaction network between genes and other substances, such as chemicals and drugs, could be obtained.

2.6. LASSO Regression and SVM-RFE Analyses

For further screening out the hub gene from the key genes, the R software glmnet package was used to set the parameter “family” to “binomial” to execute LASSO logistic regression [17]. The receiver operating characteristic curve (ROC) was used to validate the obtained LASSO model. Similarly, the SVM-RFE approach was applied to rank the features of the obtained key genes. Finally, take the intersection of the marker molecules obtained by the LASSO and SVM as the hub gene. To investigate the distinguishing ability of hub genes between control and case samples, we plotted the ROC curves of hub genes.

2.7. Immunoassay Based on Hub Genes

The ssGSEA algorithm was applied to compute the abundance of 28 kinds of immune cell infiltrates in all groups, and then, boxplots were drawn using R package ggplot2 and wilcox.test to screen immune cells with differences between case and control samples. Mcp_counter algorithm in immunedeconv package (version 2.0.4) was used to compute the proportion of 28 kinds of immune cells between the case and control groups, and the results were output as heat maps by pheatmap (version 1.0.12) [18, 19]. Then, the correlations between the proportions of 28 types of immune cells in the case and control groups were obtained through the Pearson correlation analysis and visualized. Finally, the Pearson analysis was further performed to analyze the intimacy between hub genes and the 28 kinds of immune cells.

2.8. Competing Endogenous RNA (ceRNA) Network Construction

First, we used the miRWalk website (http://mirwalk.umm.uni-heidelberg.de/) to predict miRNAs for hub genes, and then, the starbase website (http://starbase.sysu.edu.cn/) was used to predict lncRNAs based on the predicted miRNAs [20]. The lncRNA-miRNA screening criteria were clipExpNum > 10. The ceRNA network maps were visualized by Cytoscape software [21].

2.9. Drug Prediction Analysis

To focus on which drugs affect the hub gene, we predicted drugs for the hub gene using the Comparative Toxicogenomics Database (CTD, http://ctdbase.org/) [22].

2.10. Verification of Gene Expression

The mRNA expression levels of hub genes were detected in 10 CPB-related blood samples including 5 cases at 1 hour or end and 5 control samples of 0 hour from the First Affiliated Hospital of Guangxi Medical University. This study was allowed by the Ethics Committee of the First Affiliated Hospital of Guangxi Medical University. All patients had approved for the use of clinical tissue for research purposes. Total RNA was isolated using RNAiso Blood (Takara No. 9113). PrimeScript RT reagent Kit with gDNA Eraser (Perfect Real Time) (Takara No. 047A) was used for first-strand cDNA synthesis. For the analysis of the target genes' mRNA levels, qPCR was performed using TB Green Premix Ex Taq II (Tli RNaseH Plus) (No. 820A) according to the manufacturer's instructions (Takara). The thermocycling protocol is predenaturation at 95°C for 30 seconds and 40 cycles of 5 seconds at 95°C and 34 seconds at 60°C, and melting were done at 95°C. The relative expression of mRNA was calculated by 2−ΔΔCt method with the normalization to GAPDH [23]. The primers were synthesized by Tsingke Biotechnology, and the primer sequences are given in Table 1.

3. Results

3.1. Construction of Gene Coexpression Modules

To explore the overall correlation of all samples in the dataset, we first clustered the samples and removed outliers to ensure the accuracy of the analysis. Next, sample clustering and clinical trait heat maps were constructed (Figure 1(a)). A soft threshold was determined for the data to ensure that the gene interactions maximally conform to the scale-free distribution, and an optimal power value of 25 was determined (Figure 1(b)). Then, a total of 14 modules were aggregated on the basis of the criteria of dynamic tree cutting algorithm and MEDissThres set to 0.2 (Figure 1(c)). The correlation analysis between modules and traits found that the MEblue module had the highest correlation (Figure 1(d)), which contained 1765 genes, after excluding the gray modules that could not be classified. After screening according to the criteria of |MM| > 0.8 and |GS| > 0.8, a total of 25 key modular genes were obtained (Figure 1(e)).

3.2. Obtained Differentially Expressed Key Modular Genes and Enrichment Analysis

A total of 1450 DEGs were distinguished by comparing the control (0 h) with CPB (1 h or end) samples, of which 1390 were downregulated while the other 60 were upregulated (Figures 2(a) and 2(b)). Then, the overlap analysis between 1450 DEGs and the 25 key module genes from WGCNA analysis found 21 key genes (Figure 2(c)). Enrichment analysis identified 21 genes enriched to 6 GO entries, including developmental cell growth, regulation of monocyte differentiation, regulation of myeloid leukocyte differentiation, ERK1 and ERK2 cascade, volume-sensitive anion channel activity, and estrogen receptor binding (Figure 2(d)). Nevertheless, no KEGG pathway was found. Therefore, we performed ingenuity pathway analysis to explore potential signaling pathways, which resulted in a total of 514 significant pathways (Figure 2(e)). Analysis of the interaction network between genes and other substances such as chemicals and drugs revealed that nine genes appeared in the network (Supplementary Figure 1A-G).

3.3. Acquisition of Four Hub Genes

The 21 obtained genes were deconstructed for LASSO regression, and four signature genes (TRAF3IP2-AS1, PPARGC1B, CD4, and PDLIM5) were screened out (Figure 3(a)). The ROC curve analysis of the LASSO model showed that the AUC value was 1 (Figure 3(b)), indicating that the model performed well. Meanwhile, 17 signature genes (TRAF3IP2-AS1, CD4, POLR1B, MYC, OXNAD1, GVINP1, PDLIM5, DPYSL2, FBXO32, SLC35B4, CA5B, PPARGC1B, LRRC8B, TIAM1, MINA, ARRB1, and NEK11) were obtained using the SVM-RFE method, and the model has withstood the test (Figure 3(c)). Finally, combining the two results, this study obtained four genes (TRAF3IP2-AS1, PPARGC1B, CD4, and PDLIM5) as the most final hub genes (Figure 3(d)). ROC curves were plotted for hub genes and found that all the AUC values were 1, indicating that hub genes can accurately differentiate between normal and patient samples (Figure 3(e)).

3.4. Hub Gene Functional Enrichment Results

By GSEA analysis of four hub genes, we discovered that CD4 is mainly related to chromatin organization, mRNA processing, noncoding RNA conversion process, noncoding RNA processing, biological functions of the ribosome, and other functions (Figure 4(a)). It can be found that CD4 was enriched in antigen processing and presentation, apoptosis, hematopoietic cell lineage, neurotrophin signaling pathway, NOD-like receptor signaling pathway, etc. (Figure 4(b)). Hub gene PDLIM5A is involved in biological functions such as chromatin organization, mRNA processing, and ribosome biogenesis (Figure 4(c)) and participates in adherens junction, antigen recognition, apoptosis, and Fc γ R-mediated phagocytosis (Figure 4(d)). The PPARGC1B gene was related to noncoding RNA conversion process, noncoding RNA processing, biological functions of the ribosome, etc. (Figure 4(e)), participating in β-alanine metabolism, Fc γ R-mediated phagocytosis, hematopoietic cell lineage, N glycan biosynthesis, and other signaling pathways (Figure 4(f)). TRAF3IP2-AS1 is associated with functions such as cytoplasmic translation, mRNA processing, and noncoding RNA metabolic process (Figure 4(g)). The signaling pathways involved include antigen processing and presentation, citrate cycle, TCA cycle, hematopoietic cell lineage, spliceosome, and ribosome (Figure 4(h)). Based on the functional analysis of the above hub genes, we found that the four key genes possessed many of the same biological functions and participate in the same signaling pathways.

3.5. Association between Hub Genes and Immune

It can be found from the boxplot of ssGSEA that the proportions of activated CD4 T cell, MDSC, T follicular helper cell, and type 1 T helper cell were significantly different between the two groups. In more detail, the proportions of MDSC, T follicular helper cell, and type 1 T helper cell were higher in the control group, but the proportion of activated CD4 T cell was higher in the case group (Figure 5(a)). Moreover, the immunocyte proportional heat map also showed the similar result (Figure 5(b)). In addition, most of the immune cells were positively correlated between the case and control groups, but some of them were negatively correlated such as immature B cell-immature dendritic cell, activated B cell-immature dendritic cell, and immature B cell-activated dendritic cell (Figure 5(c)). Finally, only activated CD4 T cell, central memory CD8 T cell, MDSC, plasmacytoid dendritic cell, T follicular helper cell, and type 1 T helper cell were significant with certain hub genes. For instance, type 1 T helper cell was positively interrelated to all hub genes, and MDSC and T follicular helper cell were positively bound up with CD4, PDLIM5, and PPARGC1B, but activated CD4 T cell was negatively bound up with CD4, PDLIM5, and TRAF3IP2-AS1 (Figure 5(d)).

3.6. Hub Gene ceRNA Regulatory Network

By predicting the miRNAs and lncRNAs associated with hub genes, a total of 92 nodes (including 3 hub genes PPARGC1B, CD4, PDLIM5, 55 lncRNAs, and 34 miRNAs) as well as 268 edges were found in the ceRNA network (Figure 6).

3.7. Hub Gene Targeted Drugs

A total number of 215 drugs that matched the relationships between increase-expression and decrease-expression were predicted by the 4 hub genes in the CTD database. Additionally, the drug network contains 219 nodes and 278 edges (Figure 7).

3.8. Hub Gene Expression Validation

Analysis of the amplification levels of the four hub genes in the GSE143780 dataset disclosed that the mRNA of these genes was significantly downregulated in CPB-related congenital heart defect patients compared with controls (Figure 8(a)). Similarly, assays in clinical blood samples revealed that the expression of TRAF3IP2-AS1, PPARGC1B, CD4, and PDLIM5 was lower in CPB-related samples than in control samples (Figure 8(b)).

4. Discussion

CPB can effectively reduce the mortality of patients undergoing cardiac surgery, especially children. However, despite the continuous improvement in CPB, the synthetic circuit in CPB is still considered a foreign body in the blood circulation, which eventually leads to inflammation [24]. However, the specific mechanism by which CPB induces inflammatory responses remains incompletely understood. To explore these molecular mechanisms, in this study, we conducted a differential analysis of the GSE143780 dataset and found a total of 1,450 significantly differentially expressed genes in CPB-related samples. Furthermore, 21 key genes were obtained by WGCNA.

GO analysis of 21 key genes identified six functional items, namely, developmental cell growth, regulation of monocyte differentiation, regulation of myeloid leucocyte differentiation, ERK1 and EEK2 cascade, volume-sensitive anion channel activity, and estrogen receptor binding. Furthermore, 514 significant pathways were found in ingenuity pathway analysis. Among these pathways, volume-sensitive anion channel activity affects neuroinflammation by regulating glial cells [25]. Estrogen receptor binding can participate in inflammatory responses by regulating nuclear factor-κB, interferon, and toll-like receptor immune signaling pathways as well as immune cell differentiation; moreover, transcription factors such as SP1 and AP-1 also have mutual regulatory effects on estrogen receptors [26]. Additionally, the identification of regulation of monocyte differentiation, regulation of myeloid leukocyte differentiation, macropinocytosis signaling, and Th1 and Th2 pathways suggested that the 21 key genes directly affect the differentiation and secretion of immune cells [27, 28]. Moreover, it has been shown that the ERK1/2 signaling pathway and ephrin receptor signaling are involved in inflammatory responses in multiple organs and participate in processes such as development, activation, and migration of immune cells [28–33].

Furthermore, the 21 key module genes were obtained by WGCNA and differential expression analysis, and four genes (TRAF3IP2-AS1, PPARGC1B, CD4, and PDLIM5) were identified as the hub genes for CPB-induced inflammatory cytokine release in patients with congenital heart defects. The analysis of GSEA reveals that hub genes were involved in mRNA processing, noncoding RNA processing, apoptosis, and immune reaction. Among these, TRAF3IP2-AS1 mainly affects inflammatory responses by negatively regulating IL-17 signaling via the SRSF10-IRF1-Act1 axis. Several researches have found that IL-17 is an important factor in the pathogenesis of inflammatory and autoimmune diseases, for example, multiple sclerosis and rheumatoid arthritis [34].

PPPARGC1B and PPARGC1A are coactivators of immunomodulator PPARG. In addition, the product of the PPARGC1B gene, PGC1b, not only inhibits NLRP3 transcription but also indirectly or directly regulates the function of NLRP3. PPARGC1B suppresses IL-1β-induced inflammation by inhibiting the formation of NLRP3 inflammatory bodies [35, 36]. Moreover, CD4 is the key gene for the normal development of CD4 T lymphocytes. Some CD4 cells are regulatory T cells, which can produce TGF-β and IL-10 to inhibit inflammatory responses [37]. In the GSE143780 dataset and clinical samples, the expressions of TRAF3IP2-AS1, PPARGC1B, CD4, and PDLIM5 genes were lower in CPB-related samples, thereby supporting the possibility that the abovementioned four genes play important roles in CPB-induced inflammation.

By predicting the miRNAs and lncRNAs associated with hub genes, a total of 92 nodes (including 3 hub genes PPARGC1B, CD4, PDLIM5, 55 lncRNAs, and 34 miRNAs) as well as 268 edges were found in the ceRNA network. This could lay the foundation for the study of regulatory mechanisms of four key genes in diseases. In addition, Xu et al. found that PPARGC1B can alleviate IL-1β-mediated osteoarthritis [35]. If PDLIM5 is downregulated, some cancer progression rates will be suppressed [38]. More ceRNA-regulated mechanisms need to be further explored and validated. By prediction, the 4 hub genes were screened in the CTD database, resulting in a total of 215 drugs that matched the relationships between increase-expression and decrease-expression. The drug network contains 219 nodes and 278 edges. These analyses provide the basis for targeted therapy of key genes.

On the other hand, some shortcomings of the present study still exist. First, we only compared the gene expression between before and after CPB, but did not intensely investigate how these genes affect inflammation or the degree of the effect on inflammation. After that, the available dataset is lacking, and subsequent analysis in additional samples is required.

In conclusion, using CPB-associated cohort profile datasets and integrated bioinformatics analysis, 4 CPB-induced inflammatory cytokine release-associated hub genes, TRAF3IP2-AS1, PPARGC1B, CD4, and PDLIM5, were identified. These new potential targets may provide novel methods for preventing CPB-induced inflammation. Our future studies will aim at promulgating the potential diagnostic and prognostic value of these hub genes, which may ultimately support the translation of these targets into clinical practice.

5. Conclusion

We got four potential target genes (TRAF3IP2-AS1, PPARGC1B, CD4, and PDLIM5) about CPB-induced inflammation by analyzing the GSE143780 dataset and found that the expression of TRAF3IP2-AS1, PPARGC1B, CD4, and PDLIM5 was lower in CPB-related samples. Our research provides precious information on pathological mechanisms of CPB-induced inflammation.

Abbreviations

CPB: Cardiopulmonary bypass

GEO: Gene Expression Omnibus

lncRNA: Long noncoding RNA

miRNA: MicroRNA

iNOS: Inducible nitric oxide synthase

ROS: Reactive oxygen species

qRT-PCR: Quantitative real-time polymerase chain reaction

KEGG: Kyoto Encyclopedia of Genes and Genomes

GO: Gene Ontology

GSEA: Gene set enrichment analysis

IPA: Ingenuity pathway analysis

MDSCs: Myeloid-derived suppressor cells

NLRP3: NOD-like receptor thermal protein domain-associated protein 3

LASSO: Least absolute shrinkage and selection operator

SVM-RFE: Support vector machine recursive feature elimination

WGCNA: Gene coexpression network analysis

DEGs: Differentially expressed genes

ROC: Receiver operating characteristic curve

CTD: Comparative Toxicogenomics Database.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Ethical Approval

The experiment was approved by the Medical Ethics Committee of First Affiliated Hospital of Guangxi Medical University.

Conflicts of Interest

The authors have no conflicts of interest to declare.

Authors' Contributions

Liang Cai and Bingdong Zhang performed the experiments and wrote the paper. Cai Liang performed data analysis. All the authors approved the manuscript for publication.

Supplementary Materials

Supplementary Materials Supplementary Figure 1: (A)–(G) show the representative networks of chemicals or drugs associated with LCRC8, LCRC8B, CA5B, FBXO32, OXNAD1, PDLIM5, SLC25A43, and SLC35B4, respectively. The right side of the image shows the representation of the geometric diagram of the network diagram.

Figure 1 Screening of key modules and genes by WGCNA. (a) Clustering dendrogram of samples with trait heat map. (b) Analysis of network topology for various soft-thresholding powers. The left panel shows the scale-free fit index (y-axis) as a function of the soft-thresholding power (x-axis). The right panel displays the mean connectivity (degree, y-axis) as a function of the soft-thresholding power (x-axis). (c) Clustering dendrogram of genes, with dissimilarity based on topological overlap, together with assigned module colors. (d) Module-trait associations: each row corresponds to a module eigengene and each column to a trait. Each cell contains the corresponding correlation and P value. (e) MM and GS scatter plots for the key module blue. The vertical line is |MM| = 0.8, and the horizontal line is |GS| = 0.8. The key genes of the module are in the box in the upper right corner of the figure.

Figure 2 21 differentially expressed key module gene identification and enrichment analysis in congenital heart defects following cardiopulmonary bypass. (a) Volcano plot of the DEGs between CPB samples and control samples. Downregulated genes are green and upregulated genes come as red. (b) Heat map of the DEGs. Red is high expression and blue is low expression. (c) Venn diagram of DEGs and key module genes. (d, e) Bar chart of GO enrichment results (TOP10) and canonical pathway analysis for 21 key genes.

Figure 3 Identification of hub genes for CPB. (a) Log λ value of the differentially expressed CPB-related genes in the LASSO model (left). The most proper log λ value in the LASSO model (right). (b) The area under the ROC curve of LASSO model is 1.0. (c) Support vector machine-recursive feature elimination algorithm to screen 17 characteristic genes. (d) Venn diagram showed the 4 overlapped genes obtained by the LASSO algorithm and SVM-RFE algorithm. (e) The area under the ROC curve of all 4 hub genes is 1.0.

Figure 4 GSEA results for hub genes. (a) GO biological functions of CD4. (b) KEGG signaling pathways of CD4. (c) GO biological functions of PDLIM5A. (d) KEGG signaling pathways of PDLIM5A. (e) GO biological functions of PPARGC1B. (f) KEGG signaling pathways of PPARGC1B. (g) GO biological functions of TRAF3IP2-AS1. (h) KEGG signaling pathways of TRAF3IP2-AS1.

Figure 5 Evaluation and visualization of immune cell infiltration. (a) Boxplot of 28 types of immune cells proportions in the case and control groups. (b) Immunocyte proportional heat map of 28 immune cells between the case and control groups. (c) Correlations between immune cells in the case and control groups; the red represents the positive correlation and blue represents negative correlation. (d) Correlations between immune cells and hub genes. ∗P < 0.05 and ∗∗P < 0.01; ns represents no significant difference.

Figure 6 The ceRNA network based on 4 hub genes. Red ellipse indicate hub gene, purple quads indicate lncRNA, and green rhombus indicates miRNA.

Figure 7 Drug hub gene network. The orange-red elliptic nodes represent hub genes. The green diamond nodes represent drugs. The orange-red line shows that the drug upregulates the expression of the target gene. The blue line shows that the drug downregulates the expression of the target gene.

Figure 8 The expression levels of hub genes. (a) The expression levels of hub genes between CPB samples of 0 hour and CPB samples of 1 hour and end in the GSE143780 dataset. (b) The expression levels of hub genes in the between CPB samples of 0 hour and CPB samples of 1 hour and end in clinical samples. ∗∗∗P < 0.001 and ∗∗∗∗P < 0.0001.

Table 1 Specific primer series for hub genes.

Genes	Forward primer	Reverse primer	
CD4	5′-AACTGGAGAACAAGGAGGCAAAG-3′	5′-AAGAAGATGCCTAGCCCAATGAA-3′	
PDLIM5	5′-GAGGATCCCAGGGTGACAGTAAA-3′	5′-CGCCAGTCTTCAGTATCCTCAAT-3′	
TRAF3IP2-AS1	5′-TTCCTTCACGCCATCACAAG-3′	5′-CCTTTGTTGGTGCTCTGTCA-3′	
PPARGC1B	5′-GAGCAGACCTTGACAGTGGA-3′	5′-GCTGAGAGCTATTTCTTTGCCT-3′	
GAPDH	5′-TGTCAAGCTCATTTCCTGGTATG-3′	5′-TCTCTCTTCCTCTTGTGCTCTTG-3′
==== Refs
1 Whiting D. Yuki K. DiNardo J. A. Cardiopulmonary bypass in the pediatric population Best Practice & Research. Clinical Anaesthesiology 2015 29 2 241 256 10.1016/j.bpa.2015.03.006 2-s2.0-84947026680 26060033
2 Otto K. A. Therapeutic hypothermia applicable to cardiac surgery Veterinary Anaesthesia and Analgesia 2015 42 6 559 569 10.1111/vaa.12299 2-s2.0-84945464740 26361886
3 Graham E. M. Atz A. M. McHugh K. E. Preoperative steroid treatment does not improve markers of inflammation after cardiac surgery in neonates: results from a randomized trial The Journal of Thoracic and Cardiovascular Surgery 2014 147 3 902 908 10.1016/j.jtcvs.2013.06.010 2-s2.0-84894042066 23870160
4 Horton S. B. Butt W. W. Mullaly R. J. IL-6 and IL-8 levels after cardiopulmonary bypass are not affected by surface coating The Annals of Thoracic Surgery 1999 68 5 1751 1755 10.1016/S0003-4975(99)01066-8 2-s2.0-0032759656 10585054
5 Chew M. S. Brix-Christensen V. Ravn H. B. Effect of modified ultrafiltration on the inflammatory response in paediatric open-heart surgery: a prospective, randomized study Perfusion 2002 17 5 327 333 10.1191/0267659102pf595oa 2-s2.0-18544375284 12243435
6 Levi M. van der Poll T. Inflammation and coagulation Critical Care Medicine 2010 38 2 Supplement S26 S34 10.1097/CCM.0b013e3181c98d21 2-s2.0-74949117097 20083910
7 Aljure O. D. Fabbro M. N. Cardiopulmonary bypass and inflammation: the hidden enemy Journal of Cardiothoracic and Vascular Anesthesia 2019 33 2 346 347 10.1053/j.jvca.2018.05.030 2-s2.0-85050400102 30057250
8 Paparella D. Yau T. M. Young E. Cardiopulmonary bypass induced inflammation: pathophysiology and treatment. An update European Journal of Cardio-Thoracic Surgery 2002 21 2 232 244 10.1016/S1010-7940(01)01099-5 2-s2.0-0036168128 11825729
9 Giacinto O. Satriano U. Nenna A. Inflammatory response and endothelial dysfunction following cardiopulmonary bypass: pathophysiology and pharmacological targets Recent Patents on Inflammation & Allergy Drug Discovery 2019 13 2 158 173 10.2174/1872213X13666190724112644 31339081
10 Sasser W. C. Dabal R. J. Askenazi D. J. Prophylactic peritoneal dialysis following cardiopulmonary bypass in children is associated with decreased inflammation and improved clinical outcomes Congenital Heart Disease 2014 9 2 106 115 10.1111/chd.12072 2-s2.0-84897575352 23647999
11 Xie C. Yin Z. Liu Y. Analysis of characteristic genes and ceRNA regulation mechanism of endometriosis based on full transcriptional sequencing Frontiers in Genetics 2022 13, article 902329 10.3389/fgene.2022.902329 35938015
12 Tu L. N. Hsieh L. Kajimoto M. Shear stress associated with cardiopulmonary bypass induces expression of inflammatory cytokines and necroptosis in monocytes JCI Insight 2021 6 1 10.1172/jci.insight.141341
13 Langfelder P. Horvath S. WGCNA: an R package for weighted correlation network analysis BMC Bioinformatics 2008 9 1 p. 559 10.1186/1471-2105-9-559 2-s2.0-60549111634 19114008
14 Liu G. J. Cimmino L. Jude J. G. Pax5 loss imposes a reversible differentiation block in B progenitor acute lymphoblastic leukemia Genes & Development 2014 28 12 1337 1350 10.1101/gad.240416.114 2-s2.0-84902504214 24939936
15 Subramanian A. Tamayo P. Mootha V. K. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles Proceedings of the National Academy of Sciences of the United States of America 2005 102 43 15545 15550 10.1073/pnas.0506580102 2-s2.0-27344435774 16199517
16 Krämer A. Green J. Pollard J. Jr. Tugendreich S. Causal analysis approaches in ingenuity pathway analysis Bioinformatics 2014 30 4 523 530 10.1093/bioinformatics/btt703 2-s2.0-84894571860 24336805
17 Engebretsen S. Bohlin J. Statistical predictions with glmnet Clin Epigenetics 2019 11 1 p. 123
18 Sturm G. Finotello F. List M. Boegel S. Immunedeconv: an R package for unified access to computational methods for estimating immune cell fractions from bulk RNA-sequencing data Bioinformatics for Cancer Immunotherapy 2020 2120 New York, NY Humana 223 232 Methods in Molecular Biology 10.1007/978-1-0716-0327-7_16
19 Kolde R. Kolde M. R. Package ‘pheatmap’ 2018 1 R Package
20 Sticht C. De La Torre C. Parveen A. Gretz N. miRWalk: an online resource for prediction of microRNA binding sites PLoS One 2018 13 10, article e0206239 10.1371/journal.pone.0206239 2-s2.0-85055080177 30335862
21 Shannon P. Markiel A. Ozier O. Cytoscape: a software environment for integrated models of biomolecular interaction networks Genome Research 2003 13 11 2498 2504 10.1101/gr.1239303 2-s2.0-0242490780 14597658
22 Davis A. P. Grondin C. J. Johnson R. J. Comparative toxicogenomics database (CTD): update 2021 Nucleic Acids Research 2021 49 D1 D1138 D1143 10.1093/nar/gkaa891 33068428
23 Schmittgen T. D. Zakrajsek B. A. Mills A. G. Gorn V. Singer M. J. Reed M. W. Quantitative reverse transcription-polymerase chain reaction to study mRNA decay: comparison of endpoint and real-time methods Analytical Biochemistry 2000 285 2 194 204 10.1006/abio.2000.4753 2-s2.0-0034668118 11017702
24 Hirata Y. Cardiopulmonary bypass for pediatric cardiac surgery General Thoracic and Cardiovascular Surgery 2018 66 2 65 70 10.1007/s11748-017-0870-1 2-s2.0-85035124083 29185163
25 Kittl M. Helm K. Beyreis M. Acid- and volume-sensitive chloride currents in microglial cells International Journal of Molecular Sciences 2019 20 14 p. 3475 10.3390/ijms20143475 2-s2.0-85069990408 31311135
26 Kovats S. Estrogen receptors regulate innate immune cells and signaling pathways Cellular Immunology 2015 294 2 63 69 10.1016/j.cellimm.2015.01.018 2-s2.0-84933053329 25682174
27 Ji X. Han M. Yun Y. Li G. Sang N. Acute nitrogen dioxide (NO2) exposure enhances airway inflammation via modulating Th1/Th2 differentiation and activating JAK-STAT pathway Chemosphere 2015 120 722 728 10.1016/j.chemosphere.2014.10.039 2-s2.0-84912116064 25462318
28 Romagnani S. T-cell subsets (Th1 versus Th2) Annals of Allergy, Asthma & Immunology 2000 85 1 9 18 10.1016/s1081-1206(10)62426-x 2-s2.0-0033861859
29 Darling T. K. Lamb T. J. Emerging roles for Eph receptors and ephrin ligands in immunity Frontiers in Immunology 2019 10 p. 1473 10.3389/fimmu.2019.01473 2-s2.0-85069488879 31333644
30 Theofanous S. A. Florens M. V. Appeltans I. Ephrin-B2 signaling in the spinal cord as a player in post-inflammatory and stress-induced visceral hypersensitivity Neurogastroenterology and Motility 2020 32 4, article e13782 10.1111/nmo.13782 32004400
31 Coulthard M. G. Morgan M. Woodruff T. M. Eph/ephrin signaling in injury and inflammation The American Journal of Pathology 2012 181 5 1493 1503 10.1016/j.ajpath.2012.06.043 2-s2.0-84868106303 23021982
32 Abdel-dayem S. I. A. Khalil M. N. A. Abdelrahman E. H. El-Gohary H. M. Kamel A. S. Sesquiterpene lactones; damsin and neoambrosin suppress cytokine-mediated inflammation in complete Freund's adjuvant rat model via shutting Akt/ERK1/2/STAT3 signaling Journal of Ethnopharmacology 2021 266, article 113407 10.1016/j.jep.2020.113407
33 Zhang T. Ouyang H. Mei X. Erianin alleviates diabetic retinopathy by reducing retinal inflammation initiated by microglial cells via inhibiting hyperglycemia-mediated ERK1/2–NF-κB signaling pathway The FASEB Journal 2019 33 11 11776 11790 10.1096/fj.201802614RRR 31365278
34 He R. Wu S. Gao R. Identification of a long noncoding RNA TRAF3IP2-AS1 as key regulator of IL-17 signaling through the SRSF10-IRF1-Act1 axis in autoimmune diseases Journal of Immunology 2021 206 10 2353 2365 10.4049/jimmunol.2001223 33941656
35 Xu J. Pei Y. Lu J. lncRNA SNHG7 alleviates IL-1β-induced osteoarthritis by inhibiting miR-214-5p-mediated PPARGC1B signaling pathways International Immunopharmacology 2021 90, article 107150 10.1016/j.intimp.2020.107150 33296783
36 Merriman T. Terkeltaub R. PPARGC1B: insight into the expression of the gouty inflammation phenotype: PPARGC1B and gouty inflammation Rheumatology 2017 56 3 323 325 10.1093/rheumatology/kew453 2-s2.0-85019986555 28003496
37 Agita A. Alsagaff M. T. Inflammation, immunity, and hypertension Acta Medica Indonesiana 2017 49 2 158 165 28790231
38 Xie P. Liu M. Chen F. Long non-coding RNA AGAP2-AS1 silencing inhibits PDLIM5 expression impeding prostate cancer progression via up-regulation of microRNA-195-5p Frontiers in Genetics 2020 11 p. 1030 10.3389/fgene.2020.01030 33101368
