Dataset Information

Coexpression analysis of large cancer datasets provides insight into the cellular phenotypes of the tumour microenvironment.

ABSTRACT:

Background

Biopsies taken from individual tumours exhibit extensive differences in their cellular composition due to the inherent heterogeneity of cancers and vagaries of sample collection. As a result genes expressed in specific cell types, or associated with certain biological processes are detected at widely variable levels across samples in transcriptomic analyses. This heterogeneity also means that the level of expression of genes expressed specifically in a given cell type or process, will vary in line with the number of those cells within samples or activity of the pathway, and will therefore be correlated in their expression.

Results

Using a novel 3D network-based approach we have analysed six large human cancer microarray datasets derived from more than 1,000 individuals. Based upon this analysis, and without needing to isolate the individual cells, we have defined a broad spectrum of cell-type and pathway-specific gene signatures present in cancer expression data which were also found to be largely conserved in a number of independent datasets.

Conclusions

The conserved signature of the tumour-associated macrophage is shown to be largely-independent of tumour cell type. All stromal cell signatures have some degree of correlation with each other, since they must all be inversely correlated with the tumour component. However, viewed in the context of established tumours, the interactions between stromal components appear to be multifactorial given the level of one component e.g. vasculature, does not correlate tightly with another, such as the macrophage.

SUBMITTER: Doig TN

PROVIDER: S-EPMC3721986 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:BackgroundThe TGF-β pathway plays critical roles in numerous malignancies. Nevertheless, its potential role in prognosis prediction and regulating tumour microenvironment (TME) characteristics require further elucidation in bladder cancer (BLCA).MethodsTGF-β-related genes were comprehensively summarized from several databases. The TCGA-BLCA cohort (training cohort) was downloaded from the Cancer Genome Atlas, and the independent validation cohorts were gathered from Xiangya Hospital (Xinagya cohort) and Gene Expression Omnibus. Initially, we identified differentially expressed TGF-β genes (DEGs) between cancer and normal tissues. Subsequently, univariate Cox analysis was applied to identify prognostic DEGs, which were further used to develop the TGF-β risk score by performing LASSO and multivariate Cox analyses. Then, we studied the role of the TGF-β risk score in predicting prognosis and the TME phenotypes. In addition, the role of the TGF-β risk score in guiding precision treatments for BLCA has also been assessed.ResultsWe successfully constructed a TGF-β risk score with an independent prognostic prediction value. A high TGF-β risk score indicated an inflamed TME, which was supported by the positive relationships between the risk score, enrichment scores of anticancer immunity steps, and the infiltration levels of tumour-infiltrating immune cells. In addition, the risk score positively correlated with the expression of several immune checkpoints and the T cell inflamed score. Consistently, the risk score was positively related to the enrichment scores of most immunotherapy-positive pathways. In addition, the sensitivities of six common chemotherapeutic drugs were positively associated with the risk score. Furthermore, higher risk score indicated higher sensitivity to radiotherapy and EGFR-targeted therapy. On the contrary, patients with low-risk scores were more sensitive to targeted therapies, including the blockade of FGFR3 and WNT-β-catenin networks.ConclusionsWe first constructed and validated a TGF-β signature that could predict the prognosis and TME phenotypes for BLCA. More importantly, the TGF-β risk score could aid in individual precision treatment for BLCA.

Project description:PurposeThe aim of this study was to systematically establish a comprehensive tumour microenvironment (TME)-relevant prognostic gene and target miRNA network for breast cancer patients.MethodsBased on a large-scale screening of TME-relevant prognostic genes (760 genes) for breast cancer patients, the prognostic model was established. The primary TME prognostic genes were selected from the constructing database and verified in the testing database. The internal relationships between the potential TME prognostic genes and the prognosis of breast cancer patients were explored in depth. The associated miRNAs for the TME prognostic genes were generated, and the functions of each primary TME member were investigated in the breast cancer cell line.ResultsCompared with sibling controls, breast cancer patients showed 55 differentially expressed TME prognostic genes, of which 31 were considered as protective genes, while the remaining 24 genes were considered as risk genes. According to the lambda values of the LASSO Cox analysis, the 15 potential TME prognostic genes were as follows: ENPEP, CCDC102B, FEZ1, NOS2, SCG2, RPLP2, RELB, RGS3, EMP1, PDLIM4, EPHA3, PCDH9, VIM, GFI1, and IRF1. Among these, there was a remarkable linear internal relationship for CCDC102B but non-linear relationships for others with breast cancer patient prognosis. Using the siRNA technique, we silenced the expression of each TME prognostic gene. Seven of the 15 TME prognostic genes (NOS2, SCG2, RGS3, EMP1, PDLIM4, PCDH9, and GFI1) were involved in enhancing cell proliferation, destroying cell apoptosis, promoting cell invasion, or migration in breast cancer. Six of them (CCDC102B, RPLP2, RELB, EPHA3, VIM, and IRF1) were favourable for maintaining cell invasion or migration. Only two of them (ENPEP and FEZ1) were favourable for the processes of cell proliferation and apoptosis.ConclusionsThis integrated study hypothesised an innovative TME-associated genetic functional network for breast cancer patients. The external relationships between these TME prognostic genes and the disease were measured. Meanwhile, the internal molecular mechanisms were also investigated.

Project description:BackgroundAlterations in brain-derived neurotrophic factor (BDNF) gene expression contribute to serious pathologies such as depression, epilepsy, cancer, Alzheimer's, Huntington and Parkinson's disease. Therefore, exploring the mechanisms of BDNF regulation represents a great clinical importance. Studying BDNF expression remains difficult due to its multiple neural activity-dependent and tissue-specific promoters. Thus, microarray data could provide insight into the regulation of this complex gene. Conventional microarray co-expression analysis is usually carried out by merging the datasets or by confirming the re-occurrence of significant correlations across datasets. However, co-expression patterns can be different under various conditions that are represented by subsets in a dataset. Therefore, assessing co-expression by measuring correlation coefficient across merged samples of a dataset or by merging datasets might not capture all correlation patterns.ResultsIn our study, we performed meta-coexpression analysis of publicly available microarray data using BDNF as a "guide-gene" introducing a "subset" approach. The key steps of the analysis included: dividing datasets into subsets with biologically meaningful sample content (e.g. tissue, gender or disease state subsets); analyzing co-expression with the BDNF gene in each subset separately; and confirming co- expression links across subsets. Finally, we analyzed conservation in co-expression with BDNF between human, mouse and rat, and sought for conserved over-represented TFBSs in BDNF and BDNF-correlated genes. Correlated genes discovered in this study regulate nervous system development, and are associated with various types of cancer and neurological disorders. Also, several transcription factor identified here have been reported to regulate BDNF expression in vitro and in vivo.ConclusionThe study demonstrates the potential of the "subset" approach in co-expression conservation analysis for studying the regulation of single genes and proposes novel regulators of BDNF gene expression.

Project description:BackgroundThe cellular dynamics in the tumour microenvironment (TME) along with non-small cell lung cancer (NSCLC) progression remain unclear.MethodsMultiplex immunofluorescence test detecting 10 immune-related markers on 553 primary tumour (PT) samples of NSCLC was conducted and spatial information in TME was assessed by the StarDist depth learning model. The single-cell transcriptomic atlas of PT (n = 4) and paired tumour-draining lymph nodes (TDLNs) (n = 5 for tumour-invaded, n = 3 for tumour-free) microenvironment was profiled. Various bioinformatics analyses based on Gene Expression Omnibus, TCGA and Array-Express databases were also used to validate the discoveries.ResultsSpatial distances of CD4+ T cells-CD38+ T cells, CD4+ T cells-neutrophils and CD38+ T cells-neutrophils prolonged and they were replaced by CD163+ macrophages in PT along with tumour progression. Neutrophils showed unique stage and location-dependent prognostic effects. A high abundance of stromal neutrophils improved disease-free survival in the early-stage, whereas high intratumoural neutrophil infiltrates predicted poor prognosis in the mid-to-late-stage. Significant molecular and functional reprogramming in PT and TDLN microenvironments was observed. Diverse interaction networks mediated by neutrophils were found between positive and negative TDLNs. Five phenotypically and functionally heterogeneous subtypes of tumour-associated neutrophil (TAN) were further identified by pseudotime analysis, including TAN-0 with antigen-presenting function, TAN-1 with strong expression of interferon (IFN)-stimulated genes, the pro-tumour TAN-2 subcluster, the classical subset (TAN-3) and the pro-inflammatory subtype (TAN-4). Loss of IFN-stimulated signature and growing angiogenesis activity were discovered along the transitional trajectory. Eventually, a robust six neutrophil differentiation relevant genes-based model was established, showing that low-risk patients had longer overall survival time and may respond better to immunotherapy.ConclusionsThe cellular composition, spatial location, molecular and functional changes in PT and TDLN microenvironments along with NSCLC progression were deciphered, highlighting the immunoregulatory roles and evolutionary heterogeneity of TANs.