Breast cancer represents the number one cause of cancer-associated mortality globally. The most aggressive molecular subtype is triple negative breast cancer (TNBC), of which limited therapeutic options are available. It is well known that breast cancer prognosis and tumor sensitivity toward immunotherapy are dictated by the tumor microenvironment. Breast cancer gene expression profiles were extracted from the METABRIC dataset and two TNBC clusters displaying unique immune features were identified. Activated immune cells formed a large proportion of cells in the high infiltration cluster, which correlated to a good prognosis. Differentially expressed genes (DEGs) extracted between two heterogeneous subtypes were used to further explore the underlying immune mechanism and to identify prognostic biomarkers. Functional enrichment analysis revealed that the DEGs were predominately related to some processes involved in activation and regulation of innate immune signaling. Using network analysis, we identified two modules in which genes were selected for further prognostic investigation. Validation by independent datasets revealed that CXCL9 and CXCL13 were good prognostic biomarkers for TNBC. We also performed comparisons between the above two genes and immune markers (CYT, APM, TILs, and TIS), as well as cell checkpoint marker expressions, and found a statistically significant correlation between them in both METABRIC and TCGA datasets. The potential of CXCL9 and CXCL13 to predict chemotherapy sensitivity was also evaluated. We found that the CXCL9 and CXCL13 were good predictors for chemotherapy and their expressions were higher in chemotherapy-responsive patients in contrast to those who were not responsive. In brief, immune infiltrate characterization on TNBC revealed heterogeneous subtypes with unique immune features allowed for the identification of informative and reliable characteristics representative of the local immune tumor microenvironment and were potential candidates to guide the management of TNBC patients.
Background: Atherosclerotic cardiovascular diseases accounted for a quarter of global deaths. Most of these fatal diseases like coronary atherosclerotic disease (CAD) and stroke occur in the advanced stage of atherosclerosis, during which candidate therapeutic targets have not been fully established. This study aims to identify hub genes and possible regulatory targets involved in treatment of advanced atherosclerotic plaques.
Material/Methods: Microarray dataset GSE43292 and GSE28829, both containing advanced atherosclerotic plaques group and early lesions group, were obtained from the Gene Expression Omnibus database. Weighted gene co-expression network analysis (WGCNA) was conducted to identify advanced plaque-related modules. Module conservation analysis was applied to assess the similarity of advanced plaque-related modules between GSE43292 and GSE28829. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of these modules were performed by Metascape. Differentially expressed genes (DEGs) were mapped into advanced plaque-related modules and module membership values of DEGs in each module were calculated to identify hub genes. Hub genes were further validated for expression in atherosclerotic samples, for distinguishing capacity of CAD and for potential functions in advanced atherosclerosis.
Results: The lightgreen module (MElightgreen) in GSE43292 and the brown module (MEbrown) in GSE28829 were identified as advanced plaque-related modules. Conservation analysis of these two modules showed high similarity. GO and KEGG enrichment analysis revealed that genes in both MElightgreen and MEbrown were enriched in immune cell activation, secretory granules, cytokine activity, and immunoinflammatory signaling. RBM47, HCK, CD53, TYROBP, and HAVCR2 were identified as common hub genes, which were validated to be upregulated in advanced atherosclerotic plaques, to well distinguish CAD patients from non-CAD people and to regulate immune cell function-related mechanisms in advanced atherosclerosis.
Conclusions: We have identified RBM47, HCK, CD53, TYROBP, and HAVCR2 as immune-responsive hub genes related to advanced plaques, which may provide potential intervention targets to treat advanced atherosclerotic plaques.
Aims: In the cancer-related research field, there is currently a major need for a greater number of valuable biomarkers to predict the prognosis of hepatocellular carcinoma (HCC). In this study, we aimed to screen hub genes related to immune cell infiltration and explore their prognostic value for HCC.
Methods: We analyzed five datasets (GSE46408, GSE57957, GSE74656, GSE76427, and GSE87630) from the Gene Expression Omnibus database to screen the differentially expressed genes (DEGs). A protein–protein interaction network of the DEGs was constructed using the Search Tool for the Retrieval of Interacting Genes; then, the hub genes were identified. Functional enrichment of the genes was performed on the Metascape website. Next, the expression of these hub genes was validated in several databases, including Oncomine, Gene Expression Profiling Interactive Analysis 2 (GEPIA2), and Human Protein Atlas. We explored the correlations between the hub genes and infiltrated immune cells in the TIMER2.0 database. The survival curves were generated in GEPIA2, and the univariate and multivariate Cox regression analyses were performed using TIMER2.0.
Results: The top ten hub genes [DNA topoisomerase II alpha (TOP2A), cyclin B2 (CCNB2), protein regulator of cytokinesis 1 (PRC1), Rac GTPase-activating protein 1 (RACGAP1), aurora kinase A (AURKA), cyclin-dependent kinase inhibitor 3 (CDKN3), nucleolar and spindle-associated protein 1 (NUSAP1), cell division cycle-associated 5 (CDCA5), abnormal spindle microtubule assembly (ASPM), and non-SMC condensin I complex subunit G (NCAPG)] were identified in subsequent analysis. These genes are most markedly enriched in cell division, suggesting their close association with tumorigenesis. Multi-database analyses validated that the hub genes were upregulated in HCC tissues. All hub genes positively correlated with several types of immune infiltration, including B cells, CD4+ T cells, macrophages, and dendritic cells. Furthermore, these hub genes served as independent prognostic factors, and the expression of these hub genes combing with the macrophage levels could help predict an unfavorable prognosis of HCC.
Conclusion: In sum, these hub genes (TOP2A, CCNB2, PRC1, RACGAP1, AURKA, CDKN3, NUSAP1, CDCA5, ASPM, and NCAPG) may be pivotal markers for prognostic prediction as well as potentially work as targets for immune-based intervention strategies in HCC.
As a data-driven dimensionality reduction and visualization tool, t-distributed stochastic neighborhood embedding (t-SNE) has been successfully applied to a variety of fields. In recent years, it has also received increasing attention for classification and regression analysis. This study presented a t-SNE based classification approach for compositional microbiome data, which enabled us to build classifiers and classify new samples in the reduced dimensional space produced by t-SNE. The Aitchison distance was employed to modify the conditional probabilities in t-SNE to account for the compositionality of microbiome data. To classify a new sample, its low-dimensional features were obtained as the weighted mean vector of its nearest neighbors in the training set. Using the low-dimensional features as input, three commonly used machine learning algorithms, logistic regression (LR), support vector machine (SVM), and decision tree (DT) were considered for classification tasks in this study. The proposed approach was applied to two disease-associated microbiome datasets, achieving better classification performance compared with the classifiers built in the original high-dimensional space. The analytic results also showed that t-SNE with Aitchison distance led to improvement of classification accuracy in both datasets. In conclusion, we have developed a t-SNE based classification approach that is suitable for compositional microbiome data and may also serve as a baseline for more complex classification models.
Cancer stem cells (CSCs), characterized by infinite proliferation and self-renewal, greatly challenge tumor therapy. Research into their plasticity, dynamic instability, and immune microenvironment interactions may help overcome this obstacle. Data on the stemness indices (mRNAsi), gene mutations, copy number variations (CNV), tumor mutation burden (TMB), and corresponding clinical characteristics were obtained from The Cancer Genome Atlas (TCGA) and UCSC Xena Browser. The infiltrating immune cells in stomach adenocarcinoma (STAD) tissues were predicted using the CIBERSORT method. Differentially expressed genes (DEGs) between the normal and tumor tissues were used to construct prognostic models with weighted gene co-expression network analysis (WGCNA) and Lasso regression. The association between cancer stemness, gene mutations, and immune responses was evaluated in STAD. A total of 6,739 DEGs were identified between the normal and tumor tissues. DEGs in the brown (containing 19 genes) and blue (containing 209 genes) co-expression modules were used to perform survival analysis based on Cox regression. A nine-gene signature prognostic model (ARHGEF38-IT1, CCDC15, CPZ, DNASE1L2, NUDT10, PASK, PLCL1, PRR5-ARHGAP8, and SYCE2) was constructed from 178 survival-related DEGs that were significantly related to overall survival, clinical characteristics, tumor microenvironment immune cells, TMB, and cancer-related pathways in STAD. Gene correlation was significant across the prognostic model, CNVs, and drug sensitivity. Our findings provide a prognostic model and highlight potential mechanisms and associated factors (immune microenvironment and mutation status) useful for targeting CSCs.