Integrated bioinformatics analysis to identify 15 hub genes in breast cancer

Hub genes associated with immune cell infiltration in breast cancer, identified through bioinformatic analyses of multiple datasets

OBJECTIVE

The aim of this study was to identify hub genes associated with immune cell infiltration in breast cancer through bioinformatic analyses of multiple datasets.

METHODS

Nonparametric (NOISeq) and robust rank aggregation-ranked parametric (EdgeR) methods were used to assess robust differentially expressed genes across multiple datasets. Protein-protein interaction network, GO, KEGG enrichment, and sub-network analyses were performed to identify immune-associated hub genes in breast cancer. Immune cell infiltration was evaluated with the CIBERSORT, XCELL, and TIMER methods. The association between the hub gene-based risk signature and survival was determined through Kaplan-Meier survival analysis, multivariate Cox analysis, and a nomogram with external verification.

RESULTS

We identified 163 robust differentially expressed genes in breast cancer through applying both nonparametric and parametric methods to multiple GEO (n = 2,212) and TCGA (n = 1,045) datasets. Integrated bioinformatic analyses further identified 10 hub genes: CXCL10, CXCL9, CXCL11, SPP1, POSTN, MMP9, DPT, COL1A1, ADAMDEC1, and RGS1. The 10 hub-gene-based risk signature significantly correlated with the prognosis of patients with breast cancer. Moreover, these hub genes were strongly associated with the extent of infiltration of CD4+ T cells, CD8+ T cells, neutrophils, macrophages, and myeloid dendritic cells into breast tumors.

CONCLUSIONS

Integrated analyses of multiple databases led to the discovery of 10 robust hub genes that together may serve as a risk factor characteristic of the immune microenvironment in breast cancer.

Integrated bioinformatics and statistical approaches to explore molecular biomarkers for breast cancer diagnosis, prognosis and therapies

Integrated bioinformatics and statistical approaches are now playing the vital role in identifying potential molecular biomarkers more accurately in presence of huge number of alternatives for disease diagnosis, prognosis and therapies by reducing time and cost compared to the wet-lab based experimental procedures. Breast cancer (BC) is one of the leading causes of cancer related deaths for women worldwide. Several dry-lab and wet-lab based studies have identified different sets of molecular biomarkers for BC. But they did not compare their results to each other so much either computationally or experimentally. In this study, an attempt was made to propose a set of molecular biomarkers that might be more effective for BC diagnosis, prognosis and therapies, by using the integrated bioinformatics and statistical approaches. At first, we identified 190 differentially expressed genes (DEGs) between BC and control samples by using the statistical LIMMA approach. Then we identified 13 DEGs (AKR1C1, IRF9, OAS1, OAS3, SLCO2A1, NT5E, NQO1, ANGPT1, FN1, ATF6B, HPGD, BCL11A, and TP53INP1) as the key genes (KGs) by protein-protein interaction (PPI) network analysis. Then we investigated the pathogenetic processes of DEGs highlighting KGs by GO terms and KEGG pathway enrichment analysis. Moreover, we disclosed the transcriptional and post-transcriptional regulatory factors of KGs by their interaction network analysis with the transcription factors (TFs) and micro-RNAs. Both supervised and unsupervised learning's including multivariate survival analysis results confirmed the strong prognostic power of the proposed KGs. Finally, we suggested KGs-guided computationally more effective seven candidate drugs (NVP-BHG712, Nilotinib, GSK2126458, YM201636, TG-02, CX-5461, AP-24534) compared to other published drugs by cross-validation with the state-of-the-art alternatives top-ranked independent receptor proteins. Thus, our findings might be played a vital role in breast cancer diagnosis, prognosis and therapies.

SOCS3 is Related to Cell Proliferation in Neuronal Tissue: An Integrated Analysis of Bioinformatics and Experiments

Glioma is the most common primary malignant tumor that occurs in the central nervous system. Gliomas are subdivided according to a combination of microscopic morphological, molecular, and genetic factors. Glioblastoma (GBM) is the most aggressive malignant tumor; however, efficient therapies or specific target molecules for GBM have not been developed. We accessed RNA-seq and clinical data from The Cancer Genome Atlas, the Chinese Glioma Genome Atlas, and the GSE16011 dataset, and identified differentially expressed genes (DEGs) that were common to both GBM and lower-grade glioma (LGG) in three independent cohorts. The biological functions of common DEGs were examined using NetworkAnalyst. To evaluate the prognostic performance of common DEGs, we performed Kaplan-Meier and Cox regression analyses. We investigated the function of SOCS3 in the central nervous system using three GBM cell lines as well as zebrafish embryos. There were 168 upregulated genes and 50 downregulated genes that were commom to both GBM and LGG. Through survival analyses, we found that SOCS3 was the only prognostic gene in all cohorts. Inhibition of SOCS3 using siRNA decreased the proliferation of GBM cell lines. We also found that the zebrafish ortholog, socs3b, was associated with brain development through the regulation of cell proliferation in neuronal tissue. While additional mechanistic studies are necessary, our results suggest that SOCS3 is an important biomarker for glioma and that SOCS3 is related to the proliferation of neuronal tissue.

A Multimodal Affinity Fusion Network for Predicting the Survival of Breast Cancer Patients

Accurate survival prediction of breast cancer holds significant meaning for improving patient care. Approaches using multiple heterogeneous modalities such as gene expression, copy number alteration, and clinical data have showed significant advantages over those with only one modality for patient survival prediction. However, existing survival prediction methods tend to ignore the structured information between patients and multimodal data. We propose a multimodal data fusion model based on a novel multimodal affinity fusion network (MAFN) for survival prediction of breast cancer by integrating gene expression, copy number alteration, and clinical data. First, a stack-based shallow self-attention network is utilized to guide the amplification of tiny lesion regions on the original data, which locates and enhances the survival-related features. Then, an affinity fusion module is proposed to map the structured information between patients and multimodal data. The module endows the network with a stronger fusion feature representation and discrimination capability. Finally, the fusion feature embedding and a specific feature embedding from a triple modal network are fused to make the classification of long-term survival or short-term survival for each patient. As expected, the evaluation results on comprehensive performance indicate that MAFN achieves better predictive performance than existing methods. Additionally, our method can be extended to the survival prediction of other cancer diseases, providing a new strategy for other diseases prognosis.

Identification of hub genes in colorectal cancer based on weighted gene co-expression network analysis and clinical data from The Cancer Genome Atlas

Colorectal cancer (CRC) is one of the most common tumors worldwide and is associated with high mortality. Here we performed bioinformatics analysis, which we validated using immunohistochemistry in order to search for hub genes that might serve as biomarkers or therapeutic targets in CRC. Based on data from The Cancer Genome Atlas (TCGA), we identified 4832 genes differentially expressed between CRC and normal samples (1562 up-regulated and 3270 down-regulated in CRC). Gene ontology (GO) analysis showed that up-regulated genes were enriched mainly in organelle fission, cell cycle regulation, and DNA replication; down-regulated genes were enriched primarily in the regulation of ion transmembrane transport and ion homeostasis. Weighted gene co-expression network analysis (WGCNA) identified eight gene modules that were associated with clinical characteristics of CRC patients, including brown and blue modules that were associated with cancer onset. Analysis of the latter two hub modules revealed the following six hub genes: adhesion G protein-coupled receptor B3 (BAI3, also known as ADGRB3), cyclin F (CCNF), cytoskeleton-associated protein 2 like (CKAP2L), diaphanous-related formin 3 (DIAPH3), oxysterol binding protein-like 3 (OSBPL3), and RERG-like protein (RERGL). Expression levels of these hub genes were associated with prognosis, based on Kaplan–Meier survival analysis of data from the Gene Expression Profiling Interactive Analysis database. Immunohistochemistry of CRC tumor tissues confirmed that OSBPL3 is up-regulated in CRC. Our findings suggest that CCNF, DIAPH3, OSBPL3, and RERGL may be useful as therapeutic targets against CRC. BAI3 and CKAP2L may be novel biomarkers of the disease.

DNAJB4 identified as a potential breast cancer marker: evidence from bioinformatics analysis and basic experiments

Background

Breast cancer (BC) is the leading cause of tumor-related death in women worldwide, but its pathogenesis is not clear. The efficient screening of new therapeutic targets for BC through bioinformatics and biological experimental techniques has become a hot topic in BC research.

Methods

The bioinformatics method was used to analyze the gene chips and obtain the hub genes, playing an important role in the development of BC. The biological processes (BP) involved in the hub genes were analyzed by Bingo, and the impact of each hub gene on disease-free survival (DFS) and overall survival (OS) in BC patients was evaluated in the Kaplan-Meier Plotter database. The expression of DNAJB4, the hub gene with the greatest degree and having an effect on the prognosis of BC patients, was detected in BC cell lines and clinicopathological specimens. And DNAJB4 was selected for further biological experiments and clinical prognosis verification.

Results

Ten hub genes including DNAJB4, the greatest degree genes, were found by bioinformatics analysis of BC gene chips. DNAJB4 expressions in both BC cell lines and clinicopathological specimens were detected and the results showed that DNAJB4 was significantly down-regulated in BC cell lines and tissues. After interfering with the expression of DNAJB4, it was found that the invasion and migration ability of MDA-MB-231 cell line was significantly enhanced in vitro. The clinical survival data of BC patients showed that patients with high DNAJB4 expression had longer DFS.

Conclusions

DNAJB4 may be a tumor suppressor gene in BC as it could regulate invasion and migration of BC cells and its expression level is related to the prognosis of BC patients. Nevertheless, further researches are still necessary to verify its role in BC so as to provide evidences for clinical guidance regarding diagnosis and treatment.

Improving existing analysis pipeline to identify and analyze cancer driver genes using multi-omics data

The cumulative of genes carrying mutations is vital for the establishment and development of cancer. However, this driver gene exploring research line has selected and used types of tools and models of analysis unsystematically and discretely. Also, the previous studies may have neglected low-frequency drivers and seldom predicted subgroup specificities of identified driver genes. In this study, we presented an improved driver gene identification and analysis pipeline that comprises the four most widely focused analyses for driver genes: enrichment analysis, clinical feature association with expression profiles of identified driver genes as well as with their functional modules, and patient stratification by existing advanced computational tools integrating multi-omics data. The improved pipeline's general usability was demonstrated straightforwardly for breast cancer, validated by some independent databases. Accordingly, 31 validated driver genes, including four novel ones, were discovered. Subsequently, we detected cancer-related significantly enriched gene ontology terms and pathways, probable drug targets, two co-expressed modules associated significantly with several clinical features, such as number of positive lymph nodes, Nottingham prognostic index, and tumor stage, and two biologically distinct groups of BRCA patients. Data and source code of the case study can be downloaded at https://github.com/hauldhut/drivergene.

Multi-Omics Analysis Detects Novel Prognostic Subgroups of Breast Cancer

The unprecedented proliferation of recent large-scale and multi-omics databases of cancers has given us many new insights into genomic and epigenomic deregulation in cancer discovery in general. However, we wonder whether or not there exists a systematic connection between copy number aberrations (CNA) and methylation (MET)? If so, what is the role of this connection in breast cancer (BRCA) tumorigenesis and progression? At the same time, the PAM50 intrinsic subtypes of BRCA have gained the most attention from BRCA experts. However, this classification system manifests its weaknesses including low accuracy as well as a possible lack of association with biological phenotypes, and even further investigations on their clinical utility were still needed. In this study, we performed an integrative analysis of three-omics profiles, CNA, MET, and mRNA expression, in two BRCA patient cohorts (one for discovery and another for validation) – to elucidate those complicated relationships. To this purpose, we first established a set of CNAcor and METcor genes, which had CNA and MET levels significantly correlated (and anti-correlated) with their corresponding expression levels, respectively. Next, to revisit the current classification of BRCA, we performed single and integrated clustering analyses using our clustering method PINSPlus. We then discovered two biologically distinct subgroups that could be an improved and refined classification system for breast cancer patients, which can be validated by a third-party data. Further studies were then performed and realized each-subgroup-specific genes and different interactions between each of the two identified subgroups with the age factor. These findings can show promise as diagnostic and prognostic values in BRCA, and a potential alternative to the PAM50 intrinsic subtypes in the future.