Integration of Gene Expression Profile Data of Human Epicardial Adipose Tissue from Coronary Artery Disease to Verification of Hub Genes and Pathways

Abstract

Background. This study aim to identify the core pathogenic genes and explore the potential molecular mechanisms of human coronary artery disease (CAD). Methodology. Two gene profiles of epicardial adipose tissue from CAD patients including GSE 18612 and GSE 64554 were downloaded and integrated by $R$ software packages. All the coexpression of deferentially expressed genes (DEGs) were picked out and analyzed by DAVID online bioinformatic tools. In addition, the DEGs were totally typed into protein-protein interaction (PPI) networks to get the interaction data among all coexpression genes. Pictures were drawn by cytoscape software with the PPI networks data. CytoHubba were used to predict the hub genes by degree analysis. Finally all the top 10 hub genes and prediction genes in Molecular complex detection were analyzed by Gene ontology and Kyoto encyclopedia of genes and genomes pathway analysis. qRT-PCR were used to identified all the 10 hub genes. Results. The top 10 hub genes calculated by the degree method were AKT1, MYC, EGFR, ACTB, CDC42, IGF1, FGF2, CXCR4, MMP2 and LYN, which relevant with the focal adhesion pathway. Module analysis revealed that the focal adhesion was also acted an important role in CAD, which was consistence with cytoHubba. All the top 10 hub genes were verified by qRT-PCR which presented that AKT1, EGFR, CDC42, FGF2, and MMP2 were significantly decreased in epicardial adipose tissue of CAD samples ($p<0.05$) and MYC, ACTB, IGF1, CXCR4, and LYN were significantly increased ($p<0.05$). Conclusions. These candidate genes could be used as potential diagnostic biomarkers and therapeutic targets of CAD.

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.

1. Introduction

Coronary artery disease (CAD) is a common cause of morbidity and mortality Worldwide [1]. Medical therapy and interventional or surgical techniques seem to save lots of patients with acute myocardial infarction in emergency phrase. However, long term outcomes remain unsatisfactory [2]. Thus further understanding of the etiology may provide potential diagnosis and therapeutic method for CAD. Recently, many reports present that epicardial adipose tissue (EAT) plays an important role in the progression of many disease by secretion various bioactive molecules [3]. In addition, finding also shows that the changing state of EAT will significantly affect the cardiac function and increase cardiovascular risk in human beings [4, 5].

CAD is a kind of coronary artery disease and most studies aim to research through the blood sample. However, EAT is found to have relevant with the heart because it is a type of visceral fat depot. The EAT widely distributes between the pericardium and myocardium which serves with metabolically activation [6]. Owing to its special location and close proximity to coronary vasculature, comparison of different express gene between EAT[7] and subcutaneous adipose tissue (SAT) may provide important information about the state of the coronary artery and maybe have potential for diagnosis and therapy for human CAD.

mRNAs have been reported to participate in the regulation of pathophysiological conditions of cardiovascular disease (CVD) [8]. However, the existing results of mRNAs profiles of CAD can not identify the core pathogenic genes and the potential molecular mechanisms of human CAD [9–11]. Therefore, we downloaded the two mRNAs profiles data and screened out the co-expression differentially expressed genes (DEGs) between GSE 18612 and GSE 64554. After analyzing by biomathematical online tools and software, 10 hub genes were finally identified which could serve as new biomarkers and therapeutic targets for human CAD.

2. Materials and Methods

2.1. Microarray Analysis

This study was conducted in accordance with the Declaration of Helsinki and was approved by the Ethics Committee of the Second Hospital of Jilin University. The two mRNA profiles were searched by CAD and EAT as well as human in gene expression omnibus (GEO) database (http://www.ncbi.nlm.nih.gov/geo). All the data were provided from GPL96 platform. GSE 18612 contained 7 EAT from CAD patients and 6 SAT from noncoronary artery disease (NC) and there were 13 EAT from CAD patients and 10 SAT from NC patients in GSE 64554. All the data were analyzed with GEO2R online tool with two classifications including EAT and SAT.

2.2. Data Processing and Identifcation of DEGs

Downregulated and upregulated genes were selected in each profile with $p<0.05$ and fold change >1.2 or < −1.2. $R$ sofware including affy package [12]and limma package [13] were used to identify the final co-expression DEGs.

2.3. Gene Ontology and Pathway Enrichment Analyses

All DEGs and hub genes as well as predict genes in Molecular Complex Detection (MCODE) were analysis by DAVID (the Database for Annotation, Visualization, and Integrated Discovery) online bio-informatics database [14]. Gene ontology (GO) enrichment including cellular component (CC), biological processes (BP) and molecular function (MF) were acquired from DAVID. Pathway analysis of all DEGs was analysised by KEGG pathways tool. All the results were selected with gene count >2 and $p<0.05$.

2.4. Integration of Protein-Protein Interaction (PPI) Network Analysis

The DEGs were totally typed into STRING (https://string-db.org/cgi/input.pl) database [15] to get PPI networks and the interaction data among all co-expression genes were picked with a significant confidence score >0.9. Then the interaction data were typed into the Cytoscape software [16] to structure a PPI network. Based on the above data, we used MCODE [17], a built-in APP in Cytoscape software, to analyze the interactive relationship of the DEGs encoding proteins and screen hub gene. The parameters of network scoring and cluster finding were set as follows: find clusters = in whole network, degree cutoff = 2, cluster finding = haircut, node score cutoff = 0.2, k-core = 2, and max depth = 100. The parameters of cytoHubba were set as follows: Hubba nodes = top 10 nodes ranked by degree, display options = check the first-stage nodes, display the shortest path and display the expanded sub network.

2.5. Quantitative Reverse Transcription-PCR (qRT-PCR) Validation and Statistical Analysis

qRT-PCR was used to verify the top 10 hug genes. PrimeScript RT reagent Kit with gDNA Eraser (TaKaRa, Japan) was used to provide the process of RNA reverse-transcribed to cDNA. Primers were designed from Primer-BLAST (https://www.ncbi.nlm.nih.gov/tools/primer-blast) and listed in Table 1. QuantStudio 7 Flex real-time PCR system (Applied Biosystems, Carlsbad, CA, USA) was used [18, 19] and samples were normalized to GAPDH. All the operations were according to the manufacturer’s instructions. The relative expression levels of each gene were calculated using 2^−ΔΔCt methods.

Table 1

The primers of the top 10 hub genes.

Gene name	Forward primer	Reverse primer
AKT1	AGCGACGTGGCTATTGTGAAG	GCCATCATTCTTGAGGAGGAAGT
MYC	GGCTCCTGGCAAAAGGTCA	CTGCGTAGTTGTGCTGATGT
EGFR	AGGCACGAGTAACAAGCTCAC	ATGAGGACATAACCAGCCACC
ACTB	CATGTACGTTGCTATCCAGGC	CTCCTTAATGTCACGCACGAT
CDC42	CCATCGGAATATGTACCGACTG	CTCAGCGGTCGTAATCTGTCA
IGF1	GCTCTTCAGTTCGTGTGTGGA	GCCTCCTTAGATCACAGCTCC
FGF2	AGAAGAGCGACCCTCACATCA	CGGTTAGCACACACTCCTTTG
CXCR4	ACTACACCGAGGAAATGGGCT	CCCACAATGCCAGTTAAGAAGA
MMP2	TACAGGATCATTGGCTACACACC	GGTCACATCGCTCCAGACT
LYN	GCTTTTGGCACCAGGAAATAGC	TCATGTCGCTGATACAGGGAA
GAPDH	CGGACCAATACGACCAAATCCG	AGCCACATCGCTCAGACACC

3. Results

3.1. Identification of DEGs in the CAD and NC Groups

All of 20 EAT from CAD patients and 16 SAT from NC patients were analyzed with $R$ software with a threshold setting at $p<0.05$ and fold change >1.2. There were 2617 and 2192 DEGs in GSE 18612 and GSE 64554 mRNAs profiles, respectively. Finally, all 179 co-expression DEGs were selected through integrating analysis, which included 105 downregulated and 74 upregulated DEGs in EAT from CAD samples.

3.2. GO Functional Enrichment Analysis

All DEGs including downregulated and upregulated co-expression genes were typed into DAVID online tool respectively. The functions of GO analysis were divided into three classifications: BP, CC, and MF (Figure 1). As shown in Figure 1 and Table 2, in the biological processes group, the upregulated DEGs were mainly enriched in B cell receptor signaling pathway, cellular response to drug, adaptive immune response, inflammatory response, and negative regulation of protein phosphorylation, and the downregulated DEGs were mainly enriched in collagen catabolic process, cellular response to amino acid stimulus, endodermal cell differentiation and cell proliferation. In the cellular component group, the upregulated DEGs were mainly enriched in plasma membrane, cytosol, extracellular space, anchored component of membrane and extracellular region, and the downregulated DEGs were mainly enriched in extracellular matrix, extracellular exosome, protein complex, cytoplasm and cell-cell adherens junction. In the molecular function group, the upregulated DEGs were mainly enriched in identical protein binding, copper ion binding, protein complex binding, integrin binding and G-protein alpha-subunit binding, and the downregulated DEGs were mainly enriched in Rho GDP-dissociation inhibitor binding, ubiquitin protein ligase binding, cadherin binding involved in cell-cell adhesion, protein domain specific binding and lipoprotein lipase activity.

[figure omitted; refer to PDF]

Table 2

The significantly enriched analysis of differentially expressed genes in coronary artery disease.

Expression	Category	Term	Description	Gene count	$P$ -value
UP-DEGs	BP	B cell receptor signaling pathway	GO:0050853	6	2.74E-06
	BP	Adaptive immune response	GO:0002250	6	3.55E-04
	BP	GO:inflammatory response	GO:0006954	8	9.15E-04
	BP	Negative regulation of protein phosphorylation	GO:0001933	4	2.02E-03
	BP	Lipopolysaccharide-mediated signaling pathway	GO:0031663	3	7.62E-03
	CC	Plasma membrane	GO:0005886	35	5.28E-06
	CC	Cytosol	GO:0005829	26	7.67E-04
	CC	Extracellular space	GO:0005615	14	2.47E-03
	CC	Anchored component of membrane	GO:0031225	4	1.05E-02
	CC	Extracellular region	GO:0005576	13	2.56E-02
	MF	Identical protein binding	GO:0042802	10	2.44E-03
	MF	Copper ion binding	GO:0005507	3	2.02E-02
	MF	Protein complex binding	GO:0032403	4	4.68E-02
	MF	Integrin binding	GO:0005178	3	6.36E-02
	MF	G-protein alpha-subunit binding	GO:0001965	2	7.18E-02
DOWN-DEGs	BP	Collagen catabolic process	GO:0030574	6	3.55E-05
	BP	GO:cellular response to amino acid stimulus	GO:0071230	5	1.60E-04
	BP	Endodermal cell differentiation	GO:0035987	4	5.09E-04
	BP	Cell proliferation	GO:0008283	9	1.35E-03
	BP	Receptor-mediated endocytosis	GO:0006898	6	4.68E-03
	CC	Extracellular matrix	GO:0031012	10	3.94E-05
	CC	Extracellular exosome	GO:0070062	29	1.09E-03
	CC	Protein complex	GO:0043234	7	2.71E-02
	CC	Cytoplasm	GO:0005737	39	3.18E-02
	CC	Cell-cell adherens junction	GO:0005913	6	3.39E-02
	MF	Rho GDP-dissociation inhibitor binding	GO:0051022	2	2.21E-02
	MF	Ubiquitin protein ligase binding	GO:0031625	6	2.21E-02
	MF	Cadherin binding involved in cell-cell adhesion	GO:0098641	6	2.30E-02
	MF	Protein domain specific binding	GO:0019904	5	2.89E-02
	MF	Lipoprotein lipase activity	GO:0004465	2	4.37E-02

3.3. Signaling Pathway Analysis

KEGG pathway analysis on DAVID online tool indicated that upregulated genes were mainly enriched in HIF-1 signaling pathway, Rap1 signaling pathway, Ras signaling pathway, Salmonella infection and NF-kappa B signaling pathway. The downregulated genes were mainly enriched in Focal adhesion, Proteoglycans in cancer, PI3K-Akt signaling pathway, Ras signaling pathway, and Neurotrophin signaling pathway (Figure 2).

[figure omitted; refer to PDF]

3.4. PPI Network and Modular Analysis

STRING database presented a result with 178 nodes and 363 edges (Figure 3) after analysis of all the DEGs. The nodes and edges were put into the cytoscape software in order to acquire a PPI network. CytoHubba analysis section showed 10 hub genes, including AKT1, MYC, EGFR, ACTB, CDC42, IGF1, FGF2, CXCR4, MMP2, and LYN after degree calculation. In these 10 hub genes, AKT1 presented with the highest degree (degree = 38). The cytoscape plugin MCODE showed the top three modules with scores of 11.429, 3.333, and 3.000 (Figure 4). MCODE 1 contained 15 nodes including IL1B, AKT1, ANGPT2, MYC, LYN, FGF2, EGFR, TFRC, AR, ACTB, MMP2, IGF1, CDC42, CXCR4, and ANGPT1 with 80 edges. Then, all the 15 prediction genes in module 1 were analyzed by GO and KEGG analysis, which mainly aimed to Focal adhesion, MAPK signaling pathway, Salmonella infection, Adherens junction and Regulation of actin cytoskeleton. Module 2 was mainly associated with the Ubiquitin mediated proteolysis and Protein processing in endoplasmic reticulum. Module 3 was mainly associated with Protein digestion and absorption, ECM-receptor interaction, Focal adhesion and PI3K-Akt signaling pathway.

[figure omitted; refer to PDF]

[figures omitted; refer to PDF]

All the top 10 hub gene were verified by qRT-PCR which presented that AKT1, EGFR, CDC42, FGF2 and MMP2 were significantly decreased in EAT of CAD samples ($p<0.05$) and MYC, ACTB, IGF1, CXCR4, and LYN were significantly increased ($p<0.05$) (Figure 5). All validation data were consistent with the microarray data and analytical results in this study.

[figure omitted; refer to PDF]

4. Discussion

Adipose tissue has been thought to act as an endocrine organ to participate and regulate the inflammatory process of CVD through paracrine or endocrine pathway. Recent research points that EAT, an atypical fat depot surrounding the heart, plays essential roles in CAD not only highly metabolic paracrine and endocrine functions, but also by its blood supply from coronary circulation [20]. Therefore some changes of EAT may directly indicate the different state of the heart as they supplies from the same coronary blood. In the meantime, studies have points that metabolic function of EAT may participate the initiation and progression of atherosclerosis [21]. So researching on the EAT can be able to provide a promising therapeutic target for the treatment of CVD including atherosclerosis.

With the fast development of sequencing technology, finding the changing genes between disease and normal tissue has become convenient. Up to now, there are many mRNAs profiles have been carried out and a lot of different genes have been detected between EAT and SAT from CAD patients [9–11]. mRNAs act as protein-coding genes, are believed to play an important role on modulating the metabolic function of adipose tissue and CAD. However, the initially change and the core pathogenic gene of CAD are not identified and the key pathogenesis of CAD has not been confirmed until now.

As different expression profiles of CAD provided various results with each other, we integrated the mRNAs profiles and analysised the co-expression DEGs. In this study, 105 downregulated DEGs and 74 upregulated DEGs were identified after integrating profile datasets. Similar to other studies [22, 23], these DEGs were dealt with GO and pathway analysis and a lot of KEGG pathways were found. However, the hub gene in CAD and the most relevant pathway and mechanism of CAD could not be identified. Therefore, we used more bioinformation tools such as PPI Network online tool, cytoscape analysis with cytoHubba and MCODE analysis for further analyses. The cytoHubba analysis showed that the top ten hub genes were AKT1, MYC, EGFR, ACTB, CDC42, IGF1, FGF2, CXCR4, MMP2, and LYN. DAVID on line analysis tool predicted that these genes were highly relevant with the focal adhesion pathway in CAD. In addition, the MCODE provided another analysis of these DEGs. The top 1 MCODE contained 15 nodes including IL1B, AKT1, ANGPT2, MYC, LYN, FGF2, EGFR, TFRC, AR, ACTB, MMP2, IGF1, CDC42, CXCR4, ANGPT1 and 80 edges. These 15 nodes also indicated that Focal adhesion might act as the most relevant pathway after analysis by DAVID. Even better, the predicted genes in top 1 MCODE also contained all the top 10 hug genes. Therefore, we aimed to analyze the hub gene AKT1 and the Focal adhesion in CAD.

Akt kinase has been identified as one of the member of AGC kinases (AMP/GMP kinase and PKC subfamily of proteins), which is reported to play an important role in cell proliferation, migration and growth [24]. Although Akt1, Akt2, and Akt3 are the 3 isoforms of Akt, they play different functions in various cell environment. Recent study shows that Akt acts critical roles in the CVD including atherosclerosis [25]. The reduction of endothelial cell migration was found when Akt phosphorylation influenced by oxidized low density lipoprotein (LDL) [26]. In addition, inhibition of Akt in macrophages was also presented with reducing atherosclerosis [27]through decreasing in proinflammation and immune cell migration. Recent study also presented that inhibition of vascular endothelial growth factor (VEGF) through inducing Akt phosphorylation pathway might protect arteries from atherosclerosis through changing the peroxisome proliferator-activated receptors (PPARs). Therefore, the toppest hub gene Akt was believed to serve as a hub gene in the initial and progressive pathological in atherosclerosis.

Akt has been reported to activate or over expression with focal adhesion kinase (FAK) pathway signal. Reports show that interaction between FAK and Akt1 is found in many cancers. This interaction is always needed for subsequent FAK autophosphorylation, activation, and translocation to the focal adhesion complex (FAC) [28]. FAK was reported to participate in regulating cell adhesion, mobility and proliferation during many cell process. The focal adhesion signaling had provided insight into the role of FAK in cardiomyocyte signal transduction and had been reported to contribute to myocardial remodeling and the progression to heart failure [29]. In addition, FAK was also reported taking part in regulating the function of endothelial barrier [30]. Low expression of FAK in endothelial cells would result in formation of a tighter endothelial monolayer through raising cell attachment and connection. In addition, disruption of FAC would also lead to Ado/HC-induced endothelial apoptosis [31] and lipid would deposition following endothelial cell damage. Study also pointed that FAK played a role in reduction of lipopolysaccharide-induced inflammation injury through inactivation of the Wnt and NF-kappaB pathways [32] and inflammation injury was also an important factor during CAD. So focal adhesion signaling might act an important role in vascular endothelial dysfunction and pro-inflammatory reactions, and might be a core change during CAD.

In summary, by integration of different high-throughput sequencing profiles and data processing as well as hub genes, these genes may have the potential to be used as drug targets and diagnostic markers of CAD. The focal adhesion maybe the key pathway in CAD. However, there are still some limitations: further experimental studies with larger sample sizes are needed to confirm the role of the focal adhesion pathway in CAD.

Additional Points

The elements in the manuscript have not been published or are under consideration for publication elsewhere.

Conflicts of Interest

The authors declare no competing financial and/or non financial interests in relation to the work described.

Authors’ Contributions

W.T.W, Q.L, B.L, H.L.P, Q.L, and Y.W wrote the paper. Z.C.Z, D.L, T.C.W, R.H.X, and K.X.L checked the References. All authors reviewed the final manuscript.

Acknowledgments

This work was supported by the project supported by the Bethune Medical Department Doctoral Postgraduate Excellent Talents Training Program Project of Jilin Province, China (Grant no. 20181201), the Project of Direct Health Project of Jilin Provincial, China (Grant no. 20170101), and Project of Direct Health Project of Jilin Provincial, China (Grant no. 20180101). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Glossary

Abbreviations

CAD:Coronary artery disease

EAT:Epicardial adipose tissue

SAT:Subcutaneous adipose tissue

DEGs:The differentially expressed genes

GEO:The gene expression omnibus

NC:Noncoronary artery disease

GO:Gene ontology;

KEGG:Kyoto encyclopedia of genes and genomes

CC:Cellular component

BP:Biological processes

MF:Molecular function

PPI:Protein-protein interaction

MCODE:Molecular complex detection

qRT-PCR:Quantitative reverse transcription-PCR

CVD:Cardiovascular disease

AGC kinases:AMP/GMP kinase and PKC subfamily of proteins

LDL:Low density lipoprotein

VEGF:Vascular endothelial growth factor

PPARs:Peroxisome proliferator-activated receptors

FAK:Focal adhesion kinase

FAC:Focal adhesion complexes

SDF1:Stromal derived factor 1

EPCs:Endothelial progenitor cells

GHSR:Growth hormone secretagogue receptor.

References

[1] M. Coutinho, H. C. Gerstein, Y. Wang, S. Yusuf, "The relationship between glucose and incident cardiovascular events. A metaregression analysis of published data from 20 studies of 95,783 individuals followed for 12.4 years," Diabetes Care, vol. 22 no. 2, pp. 233-240, DOI: 10.2337/diacare.22.2.233, 1999.

[2] M. Chang, C. W. Lee, J.-M. Ahn, R. Cavalcante, Y. Sotomi, Y. Onuma, M. Han, S.-W. Park, P. W. Serruys, S.-J. Park, "Predictors of long-term outcomes after bypass grafting versus drug-eluting stent implantation for left main or multivessel coronary artery disease," Catheterization and Cardiovascular Interventions, vol. 90 no. 2, pp. 177-185, DOI: 10.1002/ccd.26927, 2017.

[3] D. K. Song, Y. S. Hong, H. Lee, J.-Y. Oh, Y.-A. Sung, Y. Kim, "Increased epicardial adipose tissue thickness in type 2 diabetes mellitus and obesity," Diabetes & Metabolism Journal, vol. 39 no. 5, pp. 405-413, 2015.

[4] Y. Liu, W. Fu, M. Lu, S. Huai, Y. Song, Y. Wei, "Role of miRNAs in epicardial adipose tissue in CAD patients with T2DM," Biomed Research International, vol. 2016,DOI: 10.1155/2016/1629236, 2016.

[5] Q. Ji, J. Zhang, Y. Du, E. Zhu, Z. Wang, B. Que, H. Miao, S. Shi, X. Qin, Y. Zhao, Y. Zhou, F. Huang, S. Nie, "Human epicardial adipose tissue-derived and circulating secreted frizzled-related protein 4 (SFRP4) levels are increased in patients with coronary artery disease," Cardiovascular Diabetology, vol. 16 no. 1,DOI: 10.1186/s12933-017-0612-9, 2017.

[6] A. M. Noyes, "Cardiac adipose tissue and its relationship to diabetes mellitus and cardiovascular disease," World Journal of Diabetes, vol. 6 no. 5, pp. 868-876, DOI: 10.4239/wjd.v5.i6.868, 2014.

[7] G. Iacobellis, A. C. Bianco, "Epicardial adipose tissue: emerging physiological, pathophysiological and clinical features," Trends in Endocrinology & Metabolism, vol. 22 no. 11, pp. 450-457, DOI: 10.1016/j.tem.2011.07.003, 2011.

[8] A. B. Herman, M. V. Autieri, "Inflammation-regulated mRNA stability and the progression of vascular inflammatory diseases," Clinical Sciences, vol. 131 no. 22, pp. 2687-2699, DOI: 10.1042/cs20171373, 2017.

[9] J. Ribalta, A. Alipour, F. Sánchez-Cabo, J.-C. Vallvé, T. Njo, R. Álvarez, H. Janssen, A. Liem, A. Dopazo, M. Castro-Cabezas, "Differential leucocyte RNA expression in the coronary arteries compared to systemic circulation discriminates between patients with and those without coronary artery disease," Clínica e Investigación en Arteriosclerosis, vol. 29 no. 2, pp. 60-66, DOI: 10.1016/j.arteri.2016.11.003, 2017.

[10] R. Gong, M. Chen, C. Zhang, M. Chen, H. Li, "A comparison of gene expression profiles in patients with coronary artery disease, type 2 diabetes, and their coexisting conditions," Diagnostic Pathology, vol. 12 no. 1,DOI: 10.1186/s13000-017-0630-7, 2017.

[11] J. Li, C. Zhou, J. Li, Y. Wan, T. Li, P. Ma, Y. Wang, H. Sang, "Gene expression profiling of CD133-positive cells in coronary artery disease," Molecular Medicine Reports, vol. 12 no. 5, pp. 7512-7516, 2015.

[12] L. Gautier, L. Cope, B. M. Bolstad, R. A. Irizarry, "Affy—analysis of Affymetrix GeneChip data at the probe level," Bioinformatics, vol. 20 no. 3, pp. 307-315, 2004.

[13] M. E. Ritchie, B. Phipson, D. Wu, Y. Hu, C. W. Law, W. Shi, G. K. Smyth, "Limma powers differential expression analyses for RNA-sequencing and microarray studies," Nucleic Acids Research, vol. 43 no. 7, pp. e47-e47, DOI: 10.1093/nar/gkv007, 2015.

[14] D. W. Huang, B. T. Sherman, Q. Tan, J. R. Collins, W. G. Alvord, J. Roayaei, R. Stephens, M. W. Baseler, H. C. Lane, R. A. Lempicki, "The DAVID gene functional classifcation tool: a novel biological module-centric algorithm to functionally analyze large gene lists," Genome Biology, vol. 8 no. 9,DOI: 10.1186/gb-2007-8-9-r183, 2007.

[15] D. Szklarczyk, A. Franceschini, M. Kuhn, M. Simonovic, A. Roth, P. Minguez, T. Doerks, M. Stark, J. Muller, P. Bork, L. J. Jensen, C. V. Mering, "The STRING database in 2011: functional interaction networks of proteins, globally integrated and scored," Nucleic Acids Research, vol. 39 no. Database, pp. D561-D568, DOI: 10.1093/nar/gkq973, 2011.

[16] P. Shannon, "Cytoscape: a sofware Environment for integrated models of biomolecular interaction networks," Genome Research, vol. 13 no. 11, pp. 2498-2504, DOI: 10.1101/gr.1239303, 2003.

[17] G. D. Bader, C. W. Hogue, "An automated method for fnding molecular complexes in large protein interaction networks," BMC Bioinformatics, vol. 13,DOI: 10.1186/1471-2105-4-2, 2003.

[18] Z. Li, Q. Wang, G. Chen, X. Li, Q. Yang, Z. Du, M. Ren, Y. Song, G. Zhang, "Integration of gene expression profile data to screen and verify hub genes involved in osteoarthritis," BioMed Research International, vol. 2018,DOI: 10.1155/2018/9482726, 2018.

[19] Z. Li, L. Zhong, Z. Du, G. Chen, J. Shang, Q. Yang, Q. Wang, Y. Song, G. Zhang, "Network analyses of differentially expressed genes in osteoarthritis to identify hub genes," BioMed Research International, vol. 2019,DOI: 10.1155/2019/8340573, 2019.

[20] M. L. McKenney-Drake, S. D. Rodenbeck, R. S. Bruning, A. Kole, K. W. Yancey, M. Alloosh, H. S. Sacks, M. Sturek, "Epicardial adipose tissue removal potentiates outward remodeling and arrests coronary atherogenesis," The Annals of Thoracic Surgery, vol. 103 no. 5, pp. 1622-1630, DOI: 10.1016/j.athoracsur.2016.11.034, 2017.

[21] M. Vacca, M. Di Eusanio, M. Cariello, G. Graziano, S. D'Amore, F. D. Petridis, A. D'orazio, L. Salvatore, A. Tamburro, G. Folesani, D. Rutigliano, F. Pellegrini, C. Sabbà, G. Palasciano, R. Di Bartolomeo, A. Moschetta, "Integrative miRNA and whole-genome analyses of epicardial adipose tissue in patients with coronary atherosclerosis," Cardiovascular Research, vol. 109 no. 2, pp. 228-239, DOI: 10.1093/cvr/cvv266, 2016.

[22] R. Fadaei, E. Parvaz, S. Emamgholipour, N. Moradi, A. Vatannejad, M. Najafi, M. Doosti, "The mRNA expression and circulating levels of visfatin and their correlation with coronary artery disease severity and 25-hydroxyvitamin D," Hormone and Metabolic Research, vol. 48 no. 4, pp. 269-274, DOI: 10.1055/s-0035-1564133, 2016.

[23] L. Shao, P. Zhang, Y. Zhang, Q. Lu, A. Ma, "TLR3 and TLR4 as potential clinically biomarkers of cardiovascular risk in coronary artery disease (CAD) patients," Heart and Vessels, vol. 29 no. 5, pp. 690-698, DOI: 10.1007/s00380-013-0421-3, 2014.

[24] G. Manning, D. B. Whyte, R. Martinez, T. Hunter, S. Sudarsanam, "The protein kinase complement of the human genome," Science, vol. 298 no. 5600, pp. 1912-1934, DOI: 10.1126/science.1075762, 2002.

[25] Prasanna Abeyrathna, Su Yunchao, "The critical role of akt in cardiovascular function," Vascular Pharmacology, vol. 74, pp. 38-48, DOI: 10.1016/j.vph.2015.05.008, 2015.

[26] E. Chavakis, E. Dernbach, C. Hermann, U. F. Mondorf, A. M. Zeiher, S. Dimmeler, "Oxidized LDL inhibits vascular endothelial growth factor-induced endothelial cell migration by an inhibitory effect on the Akt/endothelial nitric oxide synthase pathway," Circulation, vol. 103 no. 16, pp. 2102-2107, DOI: 10.1161/01.CIR.103.16.2102, 2001.

[27] V. R. Babaev, K. E. Hebron, C. B. Wiese, C. L. Toth, L. Ding, Y. Zhang, J. M. May, S. Fazio, K. C. Vickers, M. F. Linton, "Macrophage deficiency of Akt2 reduces atherosclerosis in Ldlrnull mice," Journal of Lipid Research, vol. 55 no. 11, pp. 2296-2308, DOI: 10.1194/jlr.m050633, 2014.

[28] M. Higuchi, R. Kihara, T. Okazaki, I. Aoki, S. Suetsugu, "Akt1 promotes focal adhesion disassembly and cell motility through phosphorylation of FAK in growth factor-stimulated cells," Journal of Cell Science, vol. 126 no. 3, pp. 745-55, 2013.

[29] A. M. Samarel, "Focal adhesion signaling in heart failure," Pflugers Archiv, vol. 466 no. 6, pp. 1101-1111, DOI: 10.1007/s00424-014-1456-8, 2014.

[30] K. M. Arnold, Z. M. Goeckeler, R. B. Wysolmerski, "Loss of focal adhesion kinase enhances endothelial barrier function and increases focaladhesions," Microcirculation, vol. 20 no. 7, pp. 637-649, DOI: 10.1111/micc.12063, 2013.

[31] Q. Lu, S. Rounds, "Focal adhesion kinase and endothelial cell apoptosis," Microvascular Research, vol. 83 no. 1, pp. 56-63, DOI: 10.1016/j.mvr.2011.05.003, 2012.

[32] D. Bai, S. Cong, L. P. Zhu, "Attenuation of focal adhesion kinase reduces lipopolysaccharide-induced inflammation injury through inactivation of the wnt and NF-kB pathways in A549 Cells," Biochemistry (Moscow), vol. 82 no. 4, pp. 446-453, DOI: 10.1134/s0006297917040058, 2017.