OL Oncology Letters 1792-1074 1792-1082 D.A. Spandidos 10.3892/ol.2016.5122 OL-0-0-5122 Articles Microarray analysis of differentially-expressed genes and linker genes associated with the molecular mechanism of colorectal cancer ShenXingjie 1 * YueMeng 1 * MengFansheng 2 ZhuJingyu 1 ZhuXiaoyan 1 JiangYakun 1 Department of Gastroenterology, Jinan Central Hospital Affiliated to Shandong University, Jinan, Shandong 250013, P.R. China Department of Pharmacy, Jinan Central Hospital Affiliated to Shandong University, Jinan, Shandong 250013, P.R. China Correspondence to: Dr Meng Yue, Department of Gastroenterology, Jinan Central Hospital Affiliated to Shandong University, 105 Jiefang Road, Lixia, Jinan, Shandong 250013, P.R. China, E-mail: yuemme@163.com

Contributed equally

 11 2016 12 09 2016 12 5 3250 3258 08042015 16062016 Copyright: © Shen et al. 2016 This is an open access article distributed under the terms of the Creative Commons Attribution-NonCommercial-NoDerivs License, which permits use and distribution in any medium, provided the original work is properly cited, the use is non-commercial and no modifications or adaptations are made.

Colorectal cancer (CRC) is one of the most prevalent malignancies worldwide and remains the third leading cause of cancer-associated mortality. The present study aimed to fully elucidate the pathogenesis of CRC and identify associated genes in tumor development. Microarray GSE44076, GSE41328 and GSE44861 datasets were downloaded from the Gene Expression Omnibus database and integrated with meta-analysis. Differentially-expressed genes (DEGs) were identified from CRC samples compared with adjacent non-cancerous controls using the Limma package in R, followed by functional analysis using the Database for Annotation, Visualization, and Integrated Discovery online tool. A protein-protein interaction (PPI) network of DEGs and linker genes was constructed using NetBox software and modules were also mined. Functional annotation was performed for modules with a maximum number of nodes. Subsequent to meta-analysis to pool the data, one dataset that included 327 samples involved in 11,081 genes was obtained. A total of 697 DEGs were identified between CRC samples and adjacent non-cancerous controls. In the PPI network, modules 1 and 5 contained the maximum number of nodes. Collagen, type I, α1 (COL1A1), COL1A2 and matrix metallopeptidase 9 (MMP9) in module 1 and UDP-glucose 6-dehydrogenase (UGDH), aldehyde dehydrogenase 1 family, member A1 (ALDH1A1), fatty acid binding protein 4 (FABP4) and monoglyceride lipase (MGLL) in module 5 exhibited a high degree of connectivity. Functional analysis indicated that the genes in module 1 were involved in extracellular matrix (ECM)-associated functions and that the genes in module 5 were involved in metabolism-related functions. Overall, significant DEGs and linker genes, namely COL1A1, COL1A2, MMP9, UGDH, ALDH1A1, FABP4 and MGLL, play a crucial role in the development of CRC via regulating the ECM and cell metabolism.

 colorectal cancer linker gene protein-protein interaction module Introduction

Colorectal cancer (CRC), also known as colon cancer, is one of the most prevalent malignancies worldwide and remains the third leading cause of cancer-associated mortality (1). In 2014, ~65,000 women and 71,830 men were estimated to be diagnosed with CRC (2).

Similar to the majority of other complex tumors, CRC has been the subject of multiple studies with regards to its pathogenesis, diagnosis and therapy. Much has been elucidated about the molecular mechanism of CRC in recent years. It is widely recognized that chromosomal instability is the most common genetic abnormality to occur in CRC and has been found in almost 85% of all CRC cases (3). Key genes involved in this pathway include Kirsten rat sarcoma viral oncogene homolog (KRAS), deleted in colorectal carcinoma (DCC), SMAD family member 2 (SMAD2) and SMAD4. KRAS is a proto-oncogene that plays a critical role in the transduction of intracellular signals. The activation of KRAS by binding to guanosine triphosphate could regulate downstream mediator mitogen-activated protein kinase, which is involved in cell division (4). Additionally, SMAD2 and SMAD4 play a vital role in the transforming growth factor-β signaling pathway, which in involved in the regulation of cell proliferation, differentiation and apoptosis (5). DCC has been shown to correlate with metastasis and a poor prognosis in CRC (6). Moreover, chronic inflammation has been shown to increase the incidence of bowel cancer (7). In the process of inflammation, cyclooxygenase-2 (COX2) is a key molecule highlighted by experiments (8). Previous studies into the downstream effects of COX have illustrated that basic fibroblast growth factor and vascular endothelial growth factor are activated by COX2 via prostaglandin E2, all of which are involved in the regulation of cell proliferation and angiogenesis contributing to tumor development (9–11). However, the pathogenesis of CRC is complex and multifactorial. The complete elucidation of its etiology remains to be defined.

The present study analyzed three microarray datasets, comparing between colon tumor samples and adjacent normal mucosa tissue samples. Differentially-expressed genes (DEGs) were identified and functional annotation was performed for significant genes, followed by protein-protein interaction (PPI) network construction. The study aimed to detect the molecular mechanisms and associated genes in the development of CRC.

 Materials and methods <sec> <title>Microarray data

A total of 3 microarray datasets were downloaded from the Gene Expression Omnibus database (http://www.ncbi.nlm.nih.gov/geo/), including GSE44076 (12), GSE41328 (13) and GSE44861 (14). Expression data from GSE44076, which included colon tumor samples from 98 patients and adjacent paired normal mucosa tissues from 50 healthy donors, were obtained using platform GPL13667 (Affymetrix Human Genome U219 Arrays). Microarray data from GSE41328, which included 5 colorectal adenocarcinomas samples and 5 matched normal colon tissue, were generated with the [HG-U133_Plus_2] Affymetrix Human Genome U133 Plus 2.0 Array platform. Expression data from GSE44861, which included 56 tumor samples and 55 adjacent non-cancerous tissue samples, were obtained through the [HT_HG-U133A] Affymetrix HT Human Genome U133A Array platform.

 Data preprocessing

Prior to analysis, probe identifications in each dataset were converted into standard gene symbols. For genes with more than one probe set in the array, the average value for the probes was obtained as the expression value of the gene. By contrast, the probe set was deleted when mapped to more than one gene. As the genes were different in the 3 datasets, meta-analysis was performed of these studies, pooling the microarray data across different platforms. In the combined process, batch effects are inevitable. To adjust the data for these batch effects, the surrogate variable analysis package (15) was applied, and normalization was performed using the preprocessCore package (16) in R.

 Identification of DEGs in CRC

To identify significant DEGs in colon tumor samples compared with adjacent non-cancerous controls, preprocessed data were exported to Limma package in R language (17). An adjusted P-value was estimated using the Benjamini & Hochberg (BH) method (18). Significant DEGs were identified as those with |log 2 FC (fold-change)|>1 and an adjusted P-value of <0.05.

 Functional annotation of DEGs in CRC

Testing for functional enrichment of DEGs in CRC was performed using the Database for Annotation, Visualization, and Integrated Discovery (DAVID) online tool (19). Categories analyzed included Gene Ontology (GO) terms (20) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways (21). Data from the GO annotations was used to construct a functional enrichment network, which was visualized by the enrichment map plugin in Cytoscape (22). The BH correction for multiple testing was performed with a cutoff for an adjusted P-value of <0.05.

 PPI network construction

NetBox software, which is written in the Java language, is used to store and establishment the Human Interaction Network based on public databases consisting of Reactome (23,24), the Human Protein Reference Database (25), Memorial Sloan-Kettering Cancer Center Cancer Cell Map (26) and the National Cancer Institute-Nature Pathway Interaction Database (27). Linker genes with statistical significance, which are not differentially-expressed in colon tumors, but interact with DEGs, were obtained through mapping DEGs onto the network. Cytoscape software (28) was used to visualize the molecular interaction. Besides the PPI network under the criteria, NetBox also divided the network into modules. The modules with the maximum number of nodes in the PPI network were subjected to GO terms and Swiss-Prot and Protein Information Resource Keywords enrichment analysis with the DAVID online tool.

 Results <sec> <title>Preprocessed results and DEGs in CRC

Following meta-analysis of these 3 studies to pool microarray data across the different platforms, one dataset was obtained that included 327 samples and 11,081 genes. The dataset was preprocessed and normalized, followed by further analysis. The normalized results are shown in Fig. 1. A total of 697 genes were selected as DEGs, including 286 upregulated and 411 downregulated genes, between CRC samples and adjacent non-cancerous control.

 Significant functions and pathways of DEGs

To annotate these DEGs in the tumor samples, DAVID was used for GO function and KEGG pathway analysis, with the threshold of the adjusted P-value at <0.05. Functional enrichment networks of upregulated and downregulated DEGs are shown in Figs. 2 and 3. The results showed that upregulated DEGs were significantly enriched in cell cycle-related functions, including the cell cycle process, the regulation of the mitotic cell cycle and the regulation of cell proliferation. Downregulated DEGs were mainly enriched in homeostasis-related functions, including chemical homeostasis and cellular ion homeostasis.

 PPI network analysis

Significant DEGs and linker genes were used to construct the PPI network (Figs. 4 and 5). In the PPI network for upregulated DEGs (Fig. 4), there were 2,508 edges and 296 genes, including 140 DEGs and 156 linker genes. The network was divided into 9 modules by NetBox, in which module 1 contained the maximum number of nodes. Additionally, in the PPI network for downregulated DEGs, there were 301 edges and 165 genes, including 113 DEGs and 42 linker genes. The network was divided into 18 modules, in which module 5 contained the maximum number of nodes. In the PPI network, the hub genes were mined with the top-five degrees of connectivity in the different modules (Table I). The upregulated minichromosome maintenance complex component 7 gene in module 0, linker genes collagen, type I, α1 (COL1A1) and COL1A2, and differentially-expressed matrix metallopeptidase 9 (MMP9) in module 1, and linker genes polo-like kinase 1 and exportin 1 in module 2 exhibited a connectivity degree of >20. Downregulated genes UDP-glucose 6-dehydrogenase (UGDH), aldehyde dehydrogenase 1 family, member A1 (ALDH1A1), fatty acid binding protein 4, adipocyte (FABP4) and monoglyceride lipase (MGLL) in module 5 exhibited a connectivity degree of >20.

The functional annotation results showed that the DEGs in module 1 were mainly enriched in extracellular region-related functions and extracellular matrix (ECM)-associated functions (Table II). Downregulated DEGs in module 5 were significantly enriched in metabolic process and biosynthetic process-related functions (Table III).

 Discussion

Using a meta-analysis approach to group 3 microarray datasets, including GSE44076, GSE41328 and GSE44861, DEGs were identified in CRC mucosa compared with adjacent normal mucosa samples. The results suggested that there were 697 DEGs, including 286 upregulated genes. Functional annotation results showed that the upregulated DEGs were involved in cell cycle-related functions, in comparison with the downregulated DEGs, which were enriched in homeostasis-associated functions. In the PPI network, the linker genes COL1A1 and COL1A2, and the DEGs MMP9, UGDH, ALDH1A1, FABP4 and MGLL, which exhibited a connectivity degree of >20, participated in the development of CRC.

After the upregulated and downregulated networks were divided into multiple modules, modules 1 and 5 with the maximum number of nodes were subjected to functional annotation. In module 1, the linker genes COL1A1 and COL1A2, and the DEG MMP9 exhibited the highest degree of connectivity. COL1A1 and COL1A2, two type I collagen members, are major components of the ECM. Growing evidence has shown that the ECM plays a critical role in promoting epithelial-to-mesenchymal transition (EMT), which is associated with tumor invasion and metastasis (29). Additionally, EMT is indicated to confer tumor cell resistance to apoptosis and to promote the escape of tumor cells from the senescence process (30,31). Moreover, a pioneer study uncovered the fact that EMT has the capacity of endowing tumor cells with cancer stem cell-like characteristics, which could promote tumor development and chemoresistance (32). MMP9 (also known as gelatinase B), a member of the MMP family, has been proven to degrade various components of the ECM, including type I collagen (33). Notably, an elevated level of MMP9 has been found in CRC (34), which is consistent with the present analysis. Numerous studies have shown that MMP9 plays crucial roles in invasion, metastasis, cell proliferation and angiogenesis (35,36). Angiogenesis and cell proliferation are critically important for tumor development and metastatic spreading (37). From the results of the functional annotation in the present study, the three important genes in module 1 were mainly enriched in ECM-related functions and the ECM-receptor interaction pathway. Accordingly, COL1A1, COL1A2 and MMP9 are involved in CRC tumorigenesis and metastasis via regulation of ECM-associated functions.

In module 5, UGDH, ALDH1A1, FABP4 and MGLL were downregulated in colorectal tumor samples and were significantly involved in metabolism-related functions. UGDH is the four-electron transfer enzyme and is associated with the biosynthesis of hyaluronan (HA), which participants in tissue organization, development and cell proliferation (38). A previous study showed that elevated levels of HA are directly involved in the progression of various cancers, and UGDH has been proposed as a biomarker for prostate cancer (39). In parallel, ALDH1A1, which belongs to a superfamily of enzymes, has been identified as a crystalline in the lens and cornea (40). Notably, ALDH1A1 also plays a critical role in regulating lipid metabolism and gluconeogenesis (41). In addition, FABP4, known as a new adipokine, is involved in fatty acid trafficking from the cytoplasm to the nucleus and in lipid metabolism (42). Furthermore, FABP4 is considered as a candidate biomarker of lipodystrophy and metabolic syndrome (43). It is well known that MGLL is a member of the serine hydrolase superfamily, which hydrolyze intracellular triglyceride and cholesteryl ester into free fatty acid as an important fuel in mammals (44). More recently, MGLL was found to be abnormally expressed in aggressive human cancer, and to promote cell proliferation and tumor growth (45). It is now clear that the conversion of cells from a normal to cancerous state requires metabolic alterations, including changes in lipid metabolism and gluconeogenesis, in order to support tumor growth and survival. As a result, UGDH, ALDH1A1, FABP4 and MGLL play a key role in metabolism-related functions and regulate the tumorigenesis of CRC.

Taken together, the present results suggest that COL1A1, COL1A2 and MMP9 in module 1, and UGDH, ALDH1A1, FABP4 and MGLL in module 5 serve as key hub genes in CRC development, where the genes regulate ECM and cell metabolism-associated functions that are important for tumor growth. However, additional experiments will be required to confirm the bioinformatic results.

 References KempZThirlwellCSieberOSilverATomlinsonIAn update on the genetics of colorectal cancerHum Mol Genet13R177R185200410.1093/hmg/ddh24715358723 SiegelRDeSantisCJemalAColorectal cancer statistics, 2014CA Cancer J Clin64104117201410.3322/caac.2122024639052 GradyWMCarethersJMGenomic and epigenetic instability in colorectal cancer pathogenesisGastroenterology13510791099200810.1053/j.gastro.2008.07.07618773902 HatzivassiliouGSongKYenIBrandhuberBJAndersonDJAlvaradoRLudlamMJStokoeDGloorSLVigersGRAF inhibitors prime wild-type RAF to activate the MAPK pathway and enhance growthNature464431435201010.1038/nature0883320130576 BellamNPascheBTgf-beta signaling alterations and colon cancerCancer Treat Res15585103201010.1007/978-1-4419-6033-7_520517689 ChangSCLinJKYangSHWangHSLiAFChiCWRelationship between genetic alterations and prognosis in sporadic colorectal cancerInt J Cancer11817211727200610.1002/ijc.2156316231316 HarpazNPolydoridesADColorectal dysplasia in chronic inflammatory bowel disease: Pathology, clinical implications and pathogenesisArch Pathol Lab Med134876895201020524866 ZismanTLRubinDTColorectal cancer and dysplasia in inflammatory bowel diseaseWorld J Gastroenterol1426622669200810.3748/wjg.14.266218461651 WangDDuboisRNProstaglandins and cancerGut55115122200610.1136/gut.2004.04710016118353 EisingerALPrescottSMJonesDAStafforiniDMThe role of cyclooxygenase-2 and prostaglandins in colon cancerProstaglandins Other Lipid Mediat82147154200710.1016/j.prostaglandins.2006.05.02617164142 DohertyGAByrneSMMolloyESMalhotraVAustinSCKayEWMurrayFEFitzgeraldDJProneoplastic effects of PGE2 mediated by EP4 receptor in colorectal cancerBMC Cancer9207200910.1186/1471-2407-9-20719558693 Sanz-PamplonaRBerenguerACorderoDMollevíDGCrous-BouMSoleXParé-BrunetLGuinoESalazarRSantosCAberrant gene expression in mucosa adjacent to tumor reveals a molecular crosstalk in colon cancerMol Cancer1346201410.1186/1476-4598-13-4624597571 LinGHeXJiHShiLDavisRWZhongSReproducibility probability score-incorporating measurement variability across laboratories for gene selectionNat Biotechnol2414761477200610.1038/nbt1206-147617160039 RyanBMZanettiKARoblesAISchetterAJGoodmanJHayesRBHuangWYGunterMJYeagerMBurdetteLGermline variation in NCF4, an innate immunity gene, is associated with an increased risk of colorectal cancerInt J Cancer13413991407201410.1002/ijc.2845723982929 LeekJTJohnsonWEParkerHSFertigEJJaffeAEStoreyJDPackage ‘SVA’: Surrogate Variable AnalysisR package version 32013https://www.bioconductor.org/packages/devel/bioc/manuals/sva/man/sva.pdf BolstadBMPackage ‘preprocessCore’: A collection of pre-processing functionsR package version 12013https://www.bioconductor.org/packages/devel/bioc/manuals/preprocessCore/man/preprocessCore.pdf SmythGKLimma: Linear models for microarray dataBioinformatics and computational biology solutions using R and BioconductorSpringerNY397420200510.1007/0-387-29362-0_23 FerreiraJAThe Benjamini-Hochberg method in the case of discrete test statisticsInt J Biostat311200710.2202/1557-4679.1065 DennisGJrShermanBTHosackDAYangJGaoWLaneHCLempickiRADAVID: Database for annotation, visualization, and integrated discoveryGenome Biol4P3200310.1186/gb-2003-4-5-p312734009 AshburnerMBallCABlakeJABotsteinDButlerHCherryJMDavisAPDolinskiKDwightSSEppigJTGene ontology: Tool for the unification of biology. The Gene Ontology ConsortiumNat Genet252529200010.1038/7555610802651 KanehisaMGotoSKEGG: Kyoto encyclopedia of genes and genomesNucleic Acids Res282730200010.1093/nar/28.1.2710592173 MericoDIsserlinRBaderGDVisualizing gene-set enrichment results using the Cytoscape plug-in enrichment mapMethods Mol Biol781257277201110.1007/978-1-61779-276-2_1221877285 Joshi-TopeGGillespieMVastrikID'EustachioPSchmidtEde BonoBJassalBGopinathGRWuGRMatthewsLReactome: A knowledgebase of biological pathwaysNucleic Acid Res33D428D432200510.1093/nar/gki07215608231 MatthewsLGopinathGGillespieMCaudyMCroftDde BonoBGarapatiPHemishJHermjakobHJassalBReactome knowledgebase of human biological pathways and processesNucleic Acid Res37D619D622200910.1093/nar/gkn86318981052 PrasadTKGoelRKandasamyKKeerthikumarSKumarSMathivananSTelikicherlaDRajuRShafreenBVenugopalAHuman protein reference database-2009 updateNucleic Acid Res37D767D772200910.1093/nar/gkn89218988627 SomwarRErdjument-BromageHLarssonEShumDLockwoodWWYangGSanderCOuerfelliOTempstPJDjaballahHVarmusHESuperoxide dismutase 1 (SOD1) is a target for a small molecule identified in a screen for inhibitors of the growth of lung adenocarcinoma cell linesProc Natl Acad Sci USA1081637516380201110.1073/pnas.111355410821930909 SchaeferCFAnthonyKKrupaSBuchoffJDayMHannayTBuetowKHPID: The pathway interaction databaseNucleic Acid Res37D674D679200910.1093/nar/gkn65318832364 KohlMWieseSWarscheidBCytoscape: software for visualization and analysis of biological networksMethods Mol Biol696291303201110.1007/978-1-60761-987-1_1821063955 De CraeneBBerxGRegulatory networks defining EMT during cancer initiation and progressionNat Rev Cancer1397110201310.1038/nrc344723344542 ValdesFAlvarezAMLocascioAVegaSHerreraBFernándezMBenitoMNietoMAFabregatIThe epithelial mesenchymal transition confers resistance to the apoptotic effects of transforming growth factor Beta in fetal rat hepatocytesMol Cancer Res16878200212496370 AnsieauSBastidJDoreauAMorelAPBouchetBPThomasCFauvetFPuisieuxIDoglioniCPiccininSInduction of EMT by twist proteins as a collateral effect of tumor-promoting inactivation of premature senescenceCancer Cell147989200810.1016/j.ccr.2008.06.00518598946 ManiSAGuoWLiaoMJEatonENAyyananAZhouAYBrooksMReinhardFZhangCCShipitsinMThe epithelial-mesenchymal transition generates cells with properties of stem cellsCell133704715200810.1016/j.cell.2008.03.02718485877 RoyRYangJMosesMAMatrix metalloproteinases as novel biomarkers and potential therapeutic targets in human cancerJ Clin Oncol2752875297200910.1200/JCO.2009.23.555619738110 Turpeenniemi-HujanenTGelatinases (MMP-2 and −9) and their natural inhibitors as prognostic indicators in solid cancersBiochimie87287297200510.1016/j.biochi.2005.01.01415781315 DeryuginaEIQuigleyJPPleiotropic roles of matrix metalloproteinases in tumor angiogenesis: Contrasting, overlapping and compensatory functionsBiochim Biophys Acta1803103120201010.1016/j.bbamcr.2009.09.01719800930 GialeliCTheocharisADKaramanosNKRoles of matrix metalloproteinases in cancer progression and their pharmacological targetingFEBS J2781627201110.1111/j.1742-4658.2010.07919.x21087457 CarmelietPJainRKMolecular mechanisms and clinical applications of angiogenesisNature473298307201110.1038/nature1014421593862 ViolaMVigettiDGenasettiARizziMKarousouEMorettoPClericiMBartoliniBPallottiFDe LucaGPassiAMolecular control of the hyaluronan biosynthesisConnect Tissue Res49111114200810.1080/0300820080214840518661323 HuangDCasaleGPTianJLeleSMPisarevVMSimpsonMAHemstreetGPIIIUdp-glucose dehydrogenase as a novel field-specific candidate biomarker of prostate cancerInt J Cancer126315327201010.1002/ijc.2482019676054 ChenYKoppakaVThompsonDCVasiliouVFocus on molecules: ALDH1A1: From lens and corneal crystallin to stem cell markerExp Eye Res102105106201210.1016/j.exer.2011.04.00821536030 KieferFWOrasanuGNallamshettySBrownJDWangHLugerPQiNRBurantCFDuesterGPlutzkyJRetinaldehyde dehydrogenase 1 coordinates hepatic gluconeogenesis and lipid metabolismEndocrinology15330893099201210.1210/en.2011-210422555438 WootanMGBernlohrDAStorchJMechanism of fluorescent fatty acid transfer from adipocyte fatty acid binding protein to membranesBiochemistry3286228627199310.1021/bi00084a0338357805 KarakasSEAlmarioRUKimKSerum fatty acid binding protein 4, free fatty acids, and metabolic risk markersMetabolism5810021007200910.1016/j.metabol.2009.02.02419394980 LongJZCravattBFThe metabolic serine hydrolases and their functions in mammalian physiology and diseaseChem Rev11160226063201110.1021/cr200075y21696217 NomuraDKLongJZNiessenSHooverHSNgSWCravattBFMonoacylglycerol lipase regulates a fatty acid network that promotes cancer pathogenesisCell1404961201010.1016/j.cell.2009.11.02720079333

Results of normalization. The upper image represents data prior to normalization, the middle image represents the data adjusted for the batch effect and the image below represents data after normalization.

Functional enrichment map of upregulated differentially-expressed genes (DEGs). The nodes represent the enriched functions of the upregulated genes. Edge thickness is proportional to the number of overlapped genes between the different functions. The size of the nodes is proportional to the number of DEGs.

Functional enrichment map of downregulated differentially-expressed genes (DEGs). The nodes represent the enriched functions of the downregulated genes. Edge thickness is proportional to the number of overlapped genes between the different functions. The size of the nodes is proportional to the number of DEGs.

Protein-protein interaction network of up-regulated genes and linker genes. The coloration of the nodes is representative of the genes in different modules. Circular nodes represent differentially-expressed genes and rhombic nodes represent linker genes.

Protein-protein interaction network of downregulated genes and linker genes. The coloration of the nodes is representative of the genes in different modules. Circular nodes represent differentially-expressed genes and rhombic nodes represent linker genes.

Connectivity degree of hub genes in the top-five modules.

Module no. Hub gene Degree
Upregulated gene modules
  0 MCM7 22
  0 TP53 21
  0 ORC1L 21
  0 ORC4L 20
  0 CDC45L 20
  1 COL1A1 25
  1 COL1A2 20
  1 MMP9 20
  1 FN1 19
  1 ITGB1 18
  2 PLK1 71
  2 XPO1 68
  2 CDC20 66
  2 BIRC5 64
  2 PAFAH1B1 63
  4 RRM2 11
  4 SLC27A5 9
  4 CYP39A1 8
  4 SQLE 7
  4 LIPE 6
  5 BMP7 7
  5 INHBA 7
  5 BMP4 6
  5 INHBB 5
  5 BAMBI 5
Downregulated gene modules
  5 UGDH 26
  5 ALDH1A1 22
  5 FABP4 21
  5 MGLL 21
  5 PPAP2A 21
  6 PLG 7
  6 MEP1A 6
  6 MEP1B 6
  6 C3 5
  6 NPY 5
  8 SHBG 8
  8 SPINK7 3
  8 MT1G 2
  8 MT2A 2
  8 MT1F 1
  9 CCRL1 5
  9 VCAN 5
  9 CCL5 5
  9 CCL21 5
  9 CCL19 4
  12 CALM1 11
  12 CAV1 5
  12 SCP2 3
  12 EDNRB 3
  12 EDN3 2

Functional annotation of genes in module 1.

Category Term Count Bonferroni
Annotation cluster 1 Enrichment score: 42.348195504843254
  GOTERM_CC_FAT GO:0044421~extracellular region part 63 2.07×10−45
  SP_PIR_KEYWORDS Secreted 65 8.47×10−42
  GOTERM_CC_FAT GO:0005576~extracellular region 71 2.53×10−35
Annotation Cluster 2 Enrichment score: 41.212870504723476
  SP_PIR_KEYWORDS Signal 80 6.46×10−42
  UP_SEQ_FEATURE Signal peptide 80 2.31×10−41
  GOTERM_CC_FAT GO:0005576~extracellular region 71 2.53×10−35
Annotation cluster 3 Enrichment score: 26.15393329757373
  SP_PIR_KEYWORDS Extracellular matrix 32 3.45×10−33
  GOTERM_CC_FAT GO:0031012~extracellular matrix 35 8.66×10−28
  GOTERM_CC_FAT GO:0005578~proteinaceous extracellular matrix 33 3.76×10−26
  GOTERM_CC_FAT GO:0044420~extracellular matrix part 15 1.57×10−11

GO, Gene Ontology; SP_PIR_KEYWORDS, Swiss-Prot and Protein Information Resource Keywords; CC, cellular component.

Functional annotation of genes in module 5.

Category Term Count Benjamini
Annotation cluster 1 Enrichment score: 28.81018243601853
  SP_PIR_KEYWORDS Oxidoreductase 28 1.58×10−31
  GOTERM_BP_FAT GO:0055114~oxidation reduction 28 2.57×10−25
  GOTERM_BP_FAT GO:0008202~steroid metabolic process 21 5.77×10−25
Annotation cluster 2 Enrichment score: 17.61044037865203
  GOTERM_BP_FAT GO:0008202~steroid metabolic process 21 5.77×10−25
  SP_PIR_KEYWORDS Nadp 14 3.94×10−17
  GOTERM_BP_FAT GO:0016125~sterol metabolic process 12 2.20×10−13
  GOTERM_BP_FAT GO:0008203~cholesterol metabolic process 10 1.82×10−10
Annotation cluster 3 Enrichment score: 10.153502326316204
  GOTERM_BP_FAT GO:0006694~steroid biosynthetic process 15 5.19×10−20
  GOTERM_BP_FAT GO:0016125~sterol metabolic process 12 2.20×10−13
  GOTERM_BP_FAT GO:0008610~lipid biosynthetic process 15 6.44×10−12
  SP_PIR_KEYWORDS Steroid biosynthesis 8 1.29×10−11
  GOTERM_BP_FAT GO:0016126~sterol biosynthetic process 7 2.68×10−8
  KEGG_PATHWAY hsa00100:steroid biosynthesis 7 1.05×10−8
  SP_PIR_KEYWORDS Sterol biosynthesis 6 1.17×10−8
  SP_PIR_KEYWORDS Lipid synthesis 7 2.69×10−7
  SP_PIR_KEYWORDS Cholesterol biosynthesis 3 2.98×10−3
  GOTERM_BP_FAT GO:0006695~cholesterol biosynthetic process 3 4.89×10−2

GO, Gene Ontology; SP_PIR_KEYWORDS, Swiss-Prot and Protein Information Resource Keywords; BP, biological process; KEGG, Kyoto Encyclopedia of Genes and Genomes.