Identification of potential immunotherapy biomarkers for breast cancer by bioinformatics analysis

Breast cancer is a serious malignancy with a high incidence worldwide and a tendency to relapse. We used integrated bioinformatics analysis to identify potential biomarkers in breast carcinoma in the present study. Microarray data, 127breast tumor samples and 23 non-tumor samples, received from the Gene Expression Omnibus (GEO) dataset; 121 differentially expressed genes (DEGs) were selected. Functional analysis using DAVID revealed that these DEGs were highly gathered in endodermal cell differentiation and proteinaceous extracellular matrix. Five bioactive compounds (prostaglandin J2, tanespimycin, semustine, 5182598, and flunarizine) were identified using Connectivity Map. We used Cytoscape software and STRING dataset to structure a protein–protein interaction (PPI) network. The expression of CD24, MMP1, SDC1, and SPP1 was much higher in breast carcinoma tissue than in Para cancerous tissues analyzed by Gene Expression Profiling Interactive Analysis (GEPIA) and ONCOMINE. Overexpression ofCD24, MMP1, SDC1, and SPP1 indicated the poor prognosis in breast carcinoma patients analyzed by Kaplan–Meier (KM) Plotter. Immunohistochemistry microarray was used to further confirm that protein expression of CD24, MMP1, SDC1, and SPP1 was much higher in tumor sections than in Para cancerous tissues. Hub genes expression at the protein level was correlated tothe breast cancer subtype and grade. Furthermore, immunity analysis showed that CD24, MMP1, SDC1, and SPP1 were potentially associated with five immune cell types infiltration (CD8+ T cells, CD4+ T cells, neutrophils, macrophages,and dendritic cells) by TIMER. Thus, this study indicates potential biomarkers that could have applications in the development of immune therapy for breast cancer. However, further studies are required for verifying these results in vivo and vitro.     

Identification of prognostic biomarkers for breast cancer brain metastases based on the bioinformatics analysis

PurposeThe prognosis of breast cancer (BC) patients who develop into brain metastases (BMs) is very poor. Thus, it is of great significance to explore the etiology of BMs in BC and identify the key genes involved in this process to improve the survival of BC patients with BMs.Patients and methodsThe gene expression data and the clinical information of BC patients were downloaded from TCGA and GEO database. Differentially expressed genes (DEGs) in TCGA-BRCA and GSE12276 were overlapped to find differentially expressed metastatic genes (DEMGs). The protein-protein interaction (PPI) network of DEMGs was constructed via STRING database. ClusterProfiler R package was applied to perform the gene ontology (GO) enrichment analysis of DEMGs. The univariate Cox regression analysis and the Kaplan-Meier (K-M) curves were plotted to screen DEMGs associated with the overall survival and the metastatic recurrence survival, which were identified as the key genes associated with the BMs in BC. The immune infiltration and the expressions of immune checkpoints for BC patients with brain relapses and BC patients with other relapses were analyzed respectively. The correlations among the expressions of key genes and the differently infiltrated immune cells or the differentially expressed immune checkpoints were calculated. The gene set enrichment analysis (GSEA) of each key gene was conducted to investigate the potential mechanisms of key genes involved in BC patients with BMs. Moreover, CTD database was used to predict the drug-gene interaction network of key genes.ResultsA total of 154 DEGs were identified in BC patients at M0 and M1 in TCGA database. A total of 667 DEGs were identified in BC patients with brain relapses and with other relapses. By overlapping these DEGs, 17 DEMGs were identified, which were enriched in the cell proliferation related biological processes and the immune related molecular functions. The univariate Cox regression analysis and the Kaplan-Meier curves revealed that CXCL9 and GPR171 were closely associated with the overall survival and the metastatic recurrence survival and were identified as key genes associated with BMs in BC. The analyses of immune infiltration and immune checkpoint expressions showed that there was a significant difference of the immune microenvironment between brain relapses and other relapses in BC. GSEA indicated that CXCL9 and GPR171 may regulate BMs in BC via the immune-related pathways.ConclusionOur study identified the key genes associated with BMs in BC patients and explore the underlying mechanisms involved in the etiology of BMs in BC. These findings may provide a promising approach for the treatments of BC patients with BMs.     

Correlation of zinc finger protein 2, a prognostic biomarker,with immune infiltrates in liver cancer

Purpose: The expression and clinical value of zinc finger protein 2 gene (ZIC2) in hepatocellular carcinoma (HCC) were analyzed by mining gene information from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases.Methods: Gene chip data sets were retrieved from GEO and TCGA and screened for differentially expressed genes in HCC. Gene expression profile interaction analysis (GEPIA) and Kaplan–Meier curves were used to analyze the relationship between differentially expressed genes (DEGs) and survival and prognosis in patients with HCC. Moreover, the Genecards database was used to extract ZIC2-related proteins and to analyze the physiological process of protein enrichment. Furthermore, the relationships between ZIC2 gene and tumor cell immune invasion and that between immune cell infiltration and the 5-year survival rate were studied using the tumor immune evaluation resource (TIMER) database.Results: Datasets from GEO and TCGA revealed that ZIC2 was differentially expressed in HCC tissues and normal tissues (P<0.05). High ZIC2 expression was associated with overall survival (OS) and progress-free survival in HCC patients. Overall, 25 ZIC2 related proteins, including Gli3, PRKDC, and rnf180 were identified and protein enrichment analysis indicated these were associated with four types of cell components, six types of cell functions, and eight types of biological processes. ZIC2 was positively correlated with immune infiltration cells in patients with HCC, and higher expression of ZIC2 mRNA CD4+T cells is associated with a better 5-year survival.Conclusion: ZIC2 gene may be used as an immune response marker in liver cancer to predict the prognosis of HCC.     

Identification of hub genes and pathways in adrenocortical carcinoma by integrated bioinformatic analysis

Adrenocortical carcinoma (ACC), a rare malignant neoplasm originating from adrenal cortical cells, has high malignancy and few treatments. Therefore, it is necessary to explore the molecular mechanism of tumorigenesis, screen and verify potential biomarkers, which will provide new clues for the treatment and diagnosis of ACC. In this paper, three gene expression profiles (GSE10927, GSE12368 and GSE90713) were downloaded from the Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) were obtained using the Limma package. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were enriched by DAVID. Protein‐protein interaction (PPI) network was evaluated by STRING database, and PPI network was constructed by Cytoscape. Finally, GEPIA was used to validate hub genes’ expression. Compared with normal adrenal tissues, 74 up‐regulated DEGs and 126 down‐regulated DEGs were found in ACC samples; GO analysis showed that up‐regulated DEGs were enriched in organelle fission, nuclear division, spindle, et al, while down‐regulated DEGs were enriched in angiogenesis, proteinaceous extracellular matrix and growth factor activity; KEGG pathway analysis showed that up‐regulated DEGs were significantly enriched in cell cycle, cellular senescence and progesterone‐mediated oocyte maturation; Nine hub genes (CCNB1, CDK1, TOP2A, CCNA2, CDKN3, MAD2L1, RACGAP1, BUB1 and CCNB2) were identified by PPI network; ACC patients with high expression of 9 hub genes were all associated with worse overall survival (OS). These hub genes and pathways might be involved in the tumorigenesis, which will offer the opportunities to develop the new therapeutic targets of ACC.     

RNA sequence analysis identified bone morphogenetic protein-2 (BMP2) as a biomarker underlying form deprivation myopia

PurposeForm deprivation myopia (FDM) is an urgent public issue characterized by pathological changes, but the underlying mechanism remained unclear. The aim was to investigate bone morphogenetic proteins (BMPs) utilizing the pathogenesis of FDM.Material and methodsGene expression omnibus (GEO) database was used to analyze one mRNA profile (GSE89325) of FDM. Sixteen retina samples (8 FDM and 8 controls) were randomly divided into seven groups for differential gene expression analysis in R. software. The gene pathway and protein-protein interaction (PPI) analysis were performed by the DAVID and STRING databases. Cytoscape was used to draw the PPI network. The gene ontology (GO) enrichment and Kyoto encyclopedia of genes and Genomes (KEGG) analysis were determined to achieve gene annotation and visualization.ResultsA total of 18420 differentially expressed genes (DEGs) were identified associated with FDM. The only non-significant gene (BEND6) was separately analyzed between two groups. Thirteen hub genes were discovered, ACVR1, ACVR2A, ACVR2B, RGMB, BMPR2, BMPR1A, BMP2, BMPR1B, CHRD, PTH, PTH1R, PTHLH, and WNT9A. The expression alteration in FDM were mainly enriched in cytokine-cytokine, and neuroactive ligand receptor interaction pathways. BMP2 was the key gene in myopia progression.ConclusionsOf clinical perspective, our findings reveal that expression of BMP2 as an underlying mechanism of FDM, providing an insight for therapeutic interventions.     

Fifteen hub genes associated with progression and prognosis of clear cell renal cell carcinoma identified by coexpression analysis

Renal cell carcinoma (RCC) is the most common type of renal tumor, and the clear cell renal cell carcinoma (ccRCC) is the most frequent subtype. In this study, our aim is to identify potential biomarkers that could effectively predict the prognosis and progression of ccRCC. First, we used The Cancer Genome Atlas (TCGA) RNA-sequencing (RNA-seq) data of ccRCC to identify 2370 differentially expressed genes (DEGs). Second, the DEGs were used to construct a coexpression network by weighted gene coexpression network analysis (WGCNA). Moreover, we identified the yellow module, which was strongly related to the histologic grade and pathological stage of ccRCC. Then, the functional annotation of the yellow module and single-samples gene-set enrichment analysis of DEGs were performed and mainly enriched in cell cycle. Subsequently, 18 candidate hub genes were screened through WGCNA and protein–protein interaction (PPI) network analysis. After verification of TCGA’s ccRCC data set, Gene Expression Omnibus (GEO) data set (GSE73731) and tissue validation, we finally identified 15 hub genes that can actually predict the progression of ccRCC. In addition, by using survival analysis, we found that patients of ccRCC with high expression of each hub gene were more likely to have poor prognosis than those with low expression. The receiver operating characteristic curve showed that each hub gene could effectively distinguish between localized and advanced ccRCC. In summary, our study indicates that 15 hub genes have great predictive value for the prognosis and progression of ccRCC, and may contribute to the exploration of the pathogenesis of ccRCC.     

Key genes and pathways in tumor-educated dendritic cells by bioinformatical analysis

Specific tumor microenvironment signaling might prevent the maturation of dendritic cells (DCs) with tolerogenic and immunosuppressive potential accounting for antigen-specific unresponsiveness in the lymphoid organs and in the periphery. In the present study, dendritic cells treated with LLC lung cancer cell or 4T1 breast cancer cell culture supernatants significantly down-regulated the expression of co-stimulatory molecules MHC-II, CD40, CD80, but up-regulated the inhibitory molecule PD-L1/L2, VISTA, and increased the messengerRNA levels of interleukin (IL)-6, arginase I, and IL-10, but decreased tumor necrosis factor-α and IL-12a. RNA was isolated from the dendritic cells with or without tumor supernatant stimulation and RNA sequencing was done. Then the differential expression genes were sorted, the candidate genes were analyzed and pathway enrichment analysis was done, and the associated protein–protein interaction network (PPI) was established. After integrated bioinformatical analysis, 405 (279 up-regulated and 126 down-regulated) consistently differential expression genes were identified. Using gene ontology and pathway analysis, it was found that differential expression genes were mainly enriched in the immune response, cell–cell interaction, hemostasis, and cell surface interactions with the vascular wall. The PPI data demonstrated that 236 nodes were classified with 1072 edges, and the most remarkable three modules involved 53 central node genes associated with cell survival, cell-substrate adhesion, chemotaxis, migration, immune response, and complement receptor mediated signaling pathway. These findings revealed the immune status of dendritic cells in the tumor environment.     

环氧化物水解酶2在肝细胞癌中的表达水平与功能作用

本研究通过公共数据和实验数据,全面分析环氧化物水解酶2(epoxide hydrolase 2, EPHX2)在肝细胞癌中的表达情况、功能作用以及预后意义。利用GEO和MitoCarta数据集,筛选肝细胞癌中呈差异表达的线粒体相关基因;利用TCGA数据库分析EPHX2及其相关基因在肝细胞癌中的表达水平;运行R包绘制Kaplan-Meier生存曲线和功能富集分析;基于STRING和GSEA构建蛋白质互作网络和基因集富集分析;荧光定量PCR和GEO数据集验证EPHX2在肝细胞癌中的表达水平。本研究共筛选得到15个在肝细胞癌中呈差异表达的线粒体相关基因。EPHX2在肝细胞癌组织中的表达水平显著降低(P<0.01)。EPHX2表达水平与肝癌患者性别、分期和级别有关,而与年龄、T分期等因素无关。与EPHX2低表达组肝癌患者相比,EPHX2高表达组肝癌患者预后较好。功能富集结果显示,EPHX2与补体途径、脂肪酸降解等信号通路有关。蛋白质互作网络结果显示,EPHX2与HAO1、AGXT、ACOX1、GSTκ1、SCP-2、CAT、CYP2C8,CYP2C9,CYP2B6,和CYP2J2等密切相关。GSEA结果显示,EPHX2低表达组与肝癌细胞增殖、肝癌复发等基因集正相关。荧光定量PCR和GEO数据集验证结果显示,EPHX2在肝细胞癌组织和肝癌细胞株中呈显著低表达。EPHX2在肝细胞癌中呈显著低表达,提示其可能在肝细胞癌发生发展过程中发挥抑癌基因作用,但具体作用机制还需进一步验证。     

Development and validation of stable ferroptosis- and pyroptosis-related signatures in predicting prognosis and immune status in breast cancer

To develop and validate the predictive effects of stable ferroptosis- and pyroptosis-related features on the prognosis and immune status of breast cancer (BC). We retrieved as well as downloaded ferroptosis- and pyroptosis-related genes from the FerrDb and GeneCards databases. The minimum absolute contraction and selection operator (LASSO) algorithm in The Cancer Genome Atlas (TCGA) was used to construct a prognostic classifier combining the above two types of prognostic genes with differential expression, and the Integrated Gene Expression (GEO) dataset was used for validation. Seventeen genes presented a close association with BC prognosis. Thirteen key prognostic genes with prognostic value were considered to construct a new expression signature for classifying patients with BC into high- and low-risk groups. Kaplan–Meier analysis revealed a worse prognosis in the high-risk group. The receiver operating characteristic (ROC) curve and multivariate Cox regression analysis identified its predictive and independent features. Immune profile analysis showed that immunosuppressive cells were upregulated in the high-risk group, and this risk model was related to immunosuppressive molecules. We successfully constructed combined features of ferroptosis and pyroptosis in BC that are closely related to prognosis, clinicopathological and immune features, chemotherapy efficacy and immunosuppressive molecules. However, further experimental studies are required to verify these findings.     

Integrated bioinformatics analysis revealed the regulation of angiogenesis by tumor cells in hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is one of the leading causes of cancer mortality, metastasis accounts for most of the cases. Angiogenesis plays an important role in cancer metastasis, but how tumor cells affect the function of endothelial cells by dictating their microRNA (miRNA) expression remains largely unknown.Differentially expressed miRNAs (DEMs) were identified through dataset downloaded from the Gene Expression Omnibus (GEO) database and analyzed by GEO2R. We then used online tools to obtain potential targets of candidate miRNAs and functional enrichment analysis, as well as the protein-protein interaction (PPI). Finally, the function of miR-302c-3p was validated through in vitro assay.In the current study, we found that HCC cells altered miRNA expression profiles of human umbilical vein endothelial cells (HUVECs) and miR-302c-3p was the most down-regulated miRNA in HUVECs when they were co-cultured with HCC-LM3 cells. Functional enrichment analysis of the candidate targets revealed that these genes were involved in epigenetic regulation of gene expression, in particular, cytosine methylation. In addition, PPI network demonstrated distinct roles of genes targeted by miR-302c-3p. Importantly, inhibition of angiogenesis, migration and permeability by the most down-regulated miR-302c-3p in HUVECs was confirmed in vitro. These findings brought us novel insight into the regulation of angiogenesis by HCC cells and provided potential targets for the development of therapeutic strategies.     

Identification and validation of key genes associated with non-small-cell lung cancer

Non-small-cell lung cancer (NSCLC) is one of the main causes of death induced by cancer globally. However, the molecular aberrations in NSCLC patients remain unclearly. In the present study, four messenger RNA microarray datasets (GSE18842, GSE40275, GSE43458, and GSE102287) were downloaded from the Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) between NSCLC tissues and adjacent lung tissues were obtained from GEO2R and the overlapping DEGs were identified. Moreover, functional and pathway enrichment were performed by Funrich, while the protein–protein interaction (PPI) network construction were obtained from STRING and hub genes were visualized and identified by Cytoscape software. Furthermore, validation, overall survival (OS) and tumor staging analysis of selected hub genes were performed by GEPIA. A total of 367 DEGs (95 upregulated and 272 downregulated) were obtained through gene integration analysis. The PPI network consisted of 94 nodes and 1036 edges in the upregulated DEGs and 272 nodes and 464 edges in the downregulated DEGs, respectively. The PPI network identified 46 upregulated and 27 downregulated hub genes among the DEGs, and six (such as CENPE, NCAPH, MYH11, LRRK2, HSD17B6, and A2M) of that have not been identified to be associated with NSCLC so far. Moreover, the expression differences of the mentioned hub genes were consistent with that in lung adenocarcinoma and lung squamous cell carcinoma in the TCGA database. Further analysis showed that all the six hub genes were associated with tumor staging except MYH11, while only the upregulated DEG CENPE was associated with the worse OS of patients with NSCLC. In conclusion, the current study showed that CENPE, NCAPH, MYH11, LRRK2, HSD17B6, and A2M might be the key genes contributed to tumorigenesis or tumor progression in NSCLC, further functional study is needed to explore the involved mechanisms.     

GNG2 acts as a tumor suppressor in breast cancer through stimulating MRAS signaling

G-protein gamma subunit 2 (GNG2) is involved in several cell signaling pathways, and is essential for cell proliferation and angiogenesis. However, the role of GNG2 in tumorigenesis and development remains unclear. In this study, 1321 differentially expressed genes (DEGs) in breast cancer (BC) tissues were screened using the GEO and TCGA databases. KEGG enrichment analysis showed that most of the enriched genes were part of the PI3K-Akt signaling pathway. We identified GNG2 from the first five DEGs, its expression was markedly reduced in all BC subtype tissues. Cox regression analysis showed that GNG2 was independently associated with overall survival in patients with luminal A and triple-negative breast cancers (TNBC). GNG2 over-expression could significantly block the cell cycle, inhibit proliferation, and promote apoptosis in BC cells in vitro. In animal studies, GNG2 over-expression inhibited the growth of BC cells. Further, we found that GNG2 significantly inhibited the activity of ERK and Akt in an MRAS-dependent manner. Importantly, GNG2 and muscle RAS oncogene homolog (MRAS) were co-localized in the cell membrane, and the fluorescence resonance energy transfer (FRET) experiment revealed that they had direct interaction. In conclusion, the interaction between GNG2 and MRAS likely inhibits Akt and ERK activity, promoting apoptosis and suppressing proliferation in BC cells. Increasing GNG2 expression or disrupting the GNG2–MRAS interaction in vivo could therefore be a potential therapeutic strategy to treat BC.Subject terms: Breast cancer, Breast cancer     

Analyzing time-series microarray data reveals key genes in spinal cord injury

Although many scholars have utilized high-throughput microarrays to delineate gene expression patterns after spinal cord injury (SCI), no study has evaluated gene changes in raphe magnus (RM) and somatomotor cortex (SMTC), two areas in brain primarily affected by SCI. In present study, we aimed to analyze the differentially expressed genes (DEGs) of RM and SMTC between SCI model and sham injured control at 4, 24 h, 7, 14, 28 days, and 3 months using microarray dataset GSE2270 downloaded from gene expression omnibus and unpaired significance analysis of microarray method. Protein–protein interaction (PPI) network was constructed for DEGs at crucial time points and significant biological functions were enriched using DAVID. The results indicated that more DEGs were identified at 14 days in RM and at 4 h/3 months in SMTC after SCI. In the PPI network for DEGs at 14 days in RM, interleukin 6, glyceraldehyde-3-phosphate dehydrogenase (GAPDH), FBJ murine osteosarcoma viral oncogene homolog (FOS), tumor necrosis factor, and nuclear receptor subfamily 3, group C, member 1 (glucocorticoid receptor) were the top 5 hub genes; In the PPI network for DEGs at 3 months in SMTC, the top 5 hub genes were ubiquitin B, Ras‐related C3 botulinum toxin substrate 1 (rho family, small GTP binding protein Rac1), FOS, Janus kinase 2 and vascular endothelial growth factor A. Hedgehog and Wnt signaling pathways were the top 2 significant pathways in RM. These hub DEGs and pathways may be underlying therapeutic targets for SCI.     

Discovery of core gene families associated with liver metastasis in colorectal cancer and regulatory roles in tumor cell immune infiltration

In this study, we aimed to uncover genes that drive the pathogenesis of liver metastasis in colorectal cancer (CRC), and identify effective genes that could serve as potential therapeutic targets for treating with colorectal liver metastasis patients based on two GEO datasets. Several bioinformatics approaches were implemented. First, differential expression analysis screened out key differentially expressed genes (DEGs) across the two GEO datasets. Based on gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses, we identified the enrichment functions and pathways of the DEGs that were associated with liver metastasis in CRC. Second, immune infiltration analysis identified key immune signature gene sets associated with CRC liver metastasis, among which two key immune gene families (CD and CCL) identified as key DEGs were filtered by protein-protein interaction (PPI) network. Some of the members in these gene families were associated with disease free survival (DFS) or overall survival (OS) in two subtypes of CRC, namely COAD and READ. Finally, functional enrichment analysis of the two gene families and their neighboring genes revealed that they were closely associated with cytokine, leukocyte proliferation and chemotaxis. These results are valuable in comprehending the pathogenesis of liver metastasis in CRC, and are of seminal importance in understanding the role of immune tumor infiltration in CRC. Our study also identified potentially effective therapeutic targets for liver metastasis in CRC including CCL20, CCL24 and CD70.     

Identification of key microRNAs and hub genes in non-small-cell lung cancer using integrative bioinformatics and functional analyses

Non-small-cell lung cancer (NSCLC) is an extremely debilitating respiratory malignancy. However, the pathogenesis of NSCLC has not been fully clarified. The main objective of our study was to identify potential microRNAs (miRNAs) and their regulatory mechanism in NSCLC. Using a systematic review, two NSCLC-associated miRNA data sets (GSE102286 and GSE56036) were obtained from Gene Expression Omnibus, and the differentially expressed miRNAs (DE-miRNAs) were accessed by GEO2R. Survival analysis of candidate DE-miRNAs was conducted using the Kaplan-Meier plotter database. To further illustrate the roles of DE-miRNAs in NSCLC, their potential target genes were predicted by miRNet and were annotated by the Database for Annotation, Visualization and Integrated Discovery (DAVID) program. Moreover, the protein-protein interaction (PPI) and miRNA-hub gene regulatory network were established using the STRING database and Cytoscape software. The function of DE-miRNAs in NSCLC cells was evaluated by transwell assay. Compared with normal tissues, a total of eight DE-miRNAs was commonly changed in two data sets. The survival analysis showed that six miRNAs (miR-21-5p, miR-31-5p, miR-708-5p, miR-30a-5p, miR-451a, and miR-126-3p) were significantly correlated with overall survival. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis indicated that target genes of upregulated miRNAs were enriched in pathways in cancer, microRNAs in cancer and proteoglycans in cancer, while the target genes of downregulated miRNAs were mainly associated with pathways in cancer, the PI3K-Akt signaling pathway and HTLV-I infection. Based on the miRNA-hub gene network and expression analysis, PTEN, EGFR, STAT3, RHOA, VEGFA, TP53, CTNNB1, and KRAS were identified as potential target genes. Furthermore, all six miRNAs exhibited significant effects on NSCLC cell invasion. These findings indicate that six DE-miRNAs and their target genes may play important roles in the pathogenesis of NSCLC, which will provide novel information for NSCLC treatments.     

Restoration of miR-29b exerts anti-cancer effects on glioblastoma

Background

Methylation plays an important role in the etiology and pathogenesis of colorectal cancer (CRC). This study aimed to identify aberrantly methylated-differentially expressed genes (DEGs) and pathways in CRC by comprehensive bioinformatics analysis.

Methods

Data of gene expression microarrays (GSE68468, GSE44076) and gene methylation microarrays (GSE29490, GSE17648) were downloaded from GEO database. Aberrantly methylated-DEGs were obtained by GEO2R. Functional and enrichment analyses of selected genes were performed using DAVID database. Protein–protein interaction (PPI) network was constructed by STRING and visualized in Cytoscape. MCODE was used for module analysis of the PPI network.

Results

Totally 411 hypomethylation-high expression genes were identified, which were enriched in biological processes of response to wounding or inflammation, cell proliferation and adhesion. Pathway enrichment showed cytokine–cytokine receptor interaction, p53 signaling and cell cycle. The top 5 hub genes of PPI network were CAD, CCND1, ATM, RB1 and MET. Additionally, 239 hypermethylation-low expression genes were identified, which demonstrated enrichment in biological processes including cell–cell signaling, nerve impulse transmission, etc. Pathway analysis indicated enrichment in calcium signaling, maturity onset diabetes of the young, cell adhesion molecules, etc. The top 5 hub genes of PPI network were EGFR, ACTA1, SST, ESR1 and DNM2. After validation in TCGA database, most hub genes still remained significant.

Conclusion

In summary, our study indicated possible aberrantly methylated-differentially expressed genes and pathways in CRC by bioinformatics analysis, which may provide novel insights for unraveling pathogenesis of CRC. Hub genes including CAD, CCND1, ATM, RB1, MET, EGFR, ACTA1, SST, ESR1 and DNM2 might serve as aberrantly methylation-based biomarkers for precise diagnosis and treatment of CRC in the future.

     

Novel key genes in triple-negative breast cancer identified by weighted gene co-expression network analysis

Triple-negative breast cancer (TNBC) is a special subtype of breast cancer (BC) with poor prognosis. Although some molecular mechanisms of TNBC have been elucidated, the efficacy of current treatments is limited. Therefore, it is urgently demanded to screen for novel biomarkers and drug targets for TNBC. In this study, we obtained four independent data sets (GSE76250, GSE31448, GSE43358, and METABRIC) from the Gene Expression Omnibus (GEO) database and the cBioPortal website. In the GSE76250 data set, 890 differentially expressed genes were identified and weighted gene co-expression network analysis was performed based on them. Then, two preserved modules associated with the KI67 score were detected. Gene ontology and pathway enrichment analyses showed genes in the modules participated in some cancer-related biological processes or pathways. Non-SMC condensin I complex subunit G (NCAPG) and ATP-binding cassette subfamily A member 9 (ABCA9) were identified as hub genes of the modules, and the significance of hub genes was validated in the GSE43358 data set. Finally, their prognostic value was assessed by survival analysis. These findings suggested that NCAPG and ABCA9 may be the key genes of TNBC. Moreover, ABCA9 was first reported in TNBC. They deserved further studies.