Molecular correlates of immune cytolytic subgroups in colorectal cancer by integrated genomics analysis.
Ontology highlight
ABSTRACT:
SUBMITTER: Roufas C
PROVIDER: S-EPMC8210146 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:BackgroundMicrosatellite unstable colorectal cancers (MSI+ CRCs) expressing PD-L1, respond to anti-PD-1 or anti-PD-L1 checkpoint blockade, whereas microsatellite-stable tumors do not respond the same. Our aim was to examine how the immune landscape relates to different aspects of the CRC's biology, including neoepitope burden.MethodsWe used TCGA data to stratify patients based on a cytolytic T-cell activity expression index and correlated immune cytolytic activity (CYT) with mutational, structural, and neoepitope features of each tumor sample. The expression of several immune checkpoints was verified in an independent cohort of 72 CRC patients, relative to their MSI status, using immunohistochemistry and RT-qPCR.ResultsCRC exhibits a range of intertumoral cytolytic T-cell activity, with lower cytolytic levels in the tumor, compared to the normal tissue. We separated CRC patients into CYT-high and CYT-low subgroups. High cytolytic activity correlated with increased mutational load in colon tumors, the count of MHC-I/-II classically defined and alternatively defined neoepitopes, high microsatellite instability and deregulated expression of several inhibitory immune checkpoints (VISTA, TIGIT, PD-1, IDO1, CTLA-4, and PD-L1, among others). Many immune checkpoint molecules (IDO1, LAG3, TIGIT, VISTA, PD-1, PD-L1 and CTLA-4) expressed significantly higher in MSI+ CRCs compared to MSS tumors. The expression of Treg markers was also significantly higher in CYT-high tumors. Both individual and simultaneous high levels of CTLA-4 and PD-L1 had a positive effect on the patients' overall survival. On the reverse, simultaneous low expression of both genes led to a significant shift towards negative effect. Assessed globally, CYT-low CRCs contained more recurrent somatic copy number alterations. PD-L1 protein was absent in most samples in the independent cohort and stained lowly in 33% of MSI CRCs. PD-L1+ CRCs stained moderately for CD8 and weakly for FOXP3. CYT-high colon tumors had higher TIL load, whereas CYT-high rectum tumors had higher TAN load compared to their CYT-low counterparts.ConclusionsOverall, we highlight the link between different genetic events and the immune microenvironment in CRC, taking into consideration the status of microsatellite instability. Our data provide further evidence that MSI+ and CYT-high tumors are better candidates for combinatorial checkpoint inhibition.
| S-EPMC6701076 | biostudies-literature
Project description:IntroductionColorectal adenoma can develop into colorectal cancer. Determining the risk of tumorigenesis in colorectal adenoma would be critical for avoiding the development of colorectal cancer; however, genomic features that could help predict the risk of tumorigenesis remain uncertain.MethodsIn this work, DNA and RNA parallel capture sequencing data covering 519 genes from colorectal adenoma and colorectal cancer samples were collected. The somatic mutation profiles were obtained from DNA sequencing data, and the expression profiles were obtained from RNA sequencing data.ResultsDespite some similarities between the adenoma samples and the cancer samples, different mutation frequencies, co-occurrences, and mutually exclusive patterns were detected in the mutation profiles of patients with colorectal adenoma and colorectal cancer. Differentially expressed genes were also detected between the two patient groups using RNA sequencing. Finally, two random forest classification models were built, one based on mutation profiles and one based on expression profiles. The models distinguished adenoma and cancer samples with accuracy levels of 81.48% and 100.00%, respectively, showing the potential of the 519-gene panel for monitoring adenoma patients in clinical practice.ConclusionThis study revealed molecular characteristics and correlations between colorectal adenoma and colorectal cancer, and it demonstrated that the 519-gene panel may be used for early monitoring of the progression of colorectal adenoma to cancer.
| S-EPMC10401844 | biostudies-literature
Project description:Recurrent medulloblastoma is a therapeutic challenge because it is almost always fatal. Studies have confirmed that medulloblastoma consists of at least four distinct subgroups. We sought to delineate subgroup-specific differences in medulloblastoma recurrence patterns.We retrospectively identified a discovery cohort of all recurrent medulloblastomas at the Hospital for Sick Children (Toronto, ON, Canada) from 1994 to 2012 (cohort 1), and established molecular subgroups using a nanoString-based assay on formalin-fixed paraffin-embedded tissues or frozen tissue. The anatomical site of recurrence (local tumour bed or leptomeningeal metastasis), time to recurrence, and survival after recurrence were assessed in a subgroup-specific manner. Two independent, non-overlapping cohorts (cohort 2: samples from patients with recurrent medulloblastomas from 13 centres worldwide, obtained between 1991 and 2012; cohort 3: samples from patients with recurrent medulloblastoma obtained at the NN Burdenko Neurosurgical Institute [Moscow, Russia] between 1994 and 2011) were analysed to confirm and validate observations. When possible, molecular subgrouping was done on tissue obtained from both the initial surgery and at recurrence.Cohort 1 consisted of 30 patients with recurrent medulloblastomas; nine with local recurrences, and 21 with metastatic recurrences. Cohort 2 consisted of 77 patients and cohort 3 of 96 patients with recurrent medulloblastoma. Subgroup affiliation remained stable at recurrence in all 34 cases with available matched primary and recurrent pairs (five pairs from cohort 1 and 29 pairs from cohort 2 [15 SHH, five group 3, 14 group 4]). This finding was validated in 17 pairs from cohort 3. When analysed in a subgroup-specific manner, local recurrences in cohort 1 were more frequent in SHH tumours (eight of nine [89%]) and metastatic recurrences were more common in group 3 and group 4 tumours (17 of 20 [85%] with one WNT, p=0·0014, local vs metastatic recurrence, SHH vs group 3 vs group 4). The subgroup-specific location of recurrence was confirmed in cohort 2 (p=0·0013 for local vs metastatic recurrence, SHH vs group 3 vs group 4,), and cohort 3 (p<0·0001). Treatment with craniospinal irradiation at diagnosis was not significantly associated with the anatomical pattern of recurrence. Survival after recurrence was significantly longer in patients with group 4 tumours in cohort 1 (p=0·013) than with other subgroups, which was confirmed in cohort 2 (p=0·0075), but not cohort 3 (p=0·70).Medulloblastoma does not change subgroup at the time of recurrence, reinforcing the stability of the four main medulloblastoma subgroups. Significant differences in the location and timing of recurrence across medulloblastoma subgroups have potential treatment ramifications. Specifically, intensified local (posterior fossa) therapy should be tested in the initial treatment of patients with SHH tumours. Refinement of therapy for patients with group 3 or group 4 tumours should focus on metastases.
| S-EPMC3953419 | biostudies-literature
Project description:Rhabdomyosarcoma (RMS) is the most common soft-tissue sarcoma in childhood. Here we studied 60 RMSs using whole-exome/-transcriptome sequencing, copy number (CN) and DNA methylome analyses to unravel the genetic/epigenetic basis of RMS. On the basis of methylation patterns, RMS is clustered into four distinct subtypes, which exhibits remarkable correlation with mutation/CN profiles, histological phenotypes and clinical behaviours. A1 and A2 subtypes, especially A1, largely correspond to alveolar histology with frequent PAX3/7 fusions and alterations in cell cycle regulators. In contrast, mostly showing embryonal histology, both E1 and E2 subtypes are characterized by high frequency of CN alterations and/or allelic imbalances, FGFR4/RAS/AKT pathway mutations and PTEN mutations/methylation and in E2, also by p53 inactivation. Despite the better prognosis of embryonal RMS, patients in the E2 are likely to have a poor prognosis. Our results highlight the close relationships of the methylation status and gene mutations with the biological behaviour in RMS.
| S-EPMC4506514 | biostudies-literature
Project description:BACKGROUND:Gene expression profiling (GEP), next-generation sequencing (NGS) and copy number variation (CNV) analysis have led to an increasingly detailed characterization of the genomic profiles of DLBCL. The aim of this study was to perform a fully integrated analysis of mutational, genomic, and expression profiles to refine DLBCL subtypes. A comparison of our model with two recently published integrative DLBCL classifiers was carried out, in order to best reflect the current state of genomic subtypes. METHODS:223 patients with de novo DLBCL from the prospective, multicenter and randomized LNH-03B LYSA clinical trials were included. GEP data was obtained using Affymetrix GeneChip arrays, mutational profiles were established by Lymphopanel NGS targeting 34 key genes, CNV analysis was obtained by array CGH, and FISH and IHC were performed. Unsupervised independent component analysis (ICA) was applied to GEP data and integrated analysis of multi-level molecular data associated with each component (gene signature) was performed. FINDINGS:ICA identified 38 components reflecting transcriptomic variability across our DLBCL cohort. Many of the components were closely related to well-known DLBCL features such as cell-of-origin, stromal and MYC signatures. A component linked to gain of 19q13 locus, among other genomic alterations, was significantly correlated with poor OS and PFS. Through this integrated analysis, a high degree of heterogeneity was highlighted among previously described DLBCL subtypes. INTERPRETATION:The results of this integrated analysis enable a global and multi-level view of DLBCL, as well as improve our understanding of DLBCL subgroups.
| S-EPMC6838437 | biostudies-literature
Project description:BACKGROUND:Gene expression profiling (GEP), next-generation sequencing (NGS) and copy number variation (CNV) analysis have led to an increasingly detailed characterization of the genomic profiles of DLBCL. The aim of this study was to perform a fully integrated analysis of mutational, genomic, and expression profiles to refine DLBCL subtypes. A comparison of our model with two recently published integrative DLBCL classifiers was carried out, in order to best reflect the current state of genomic subtypes. METHODS:223 patients with de novo DLBCL from the prospective, multicenter and randomized LNH-03B LYSA clinical trials were included. GEP data was obtained using Affymetrix GeneChip arrays, mutational profiles were established by Lymphopanel NGS targeting 34 key genes, CNV analysis was obtained by array CGH, and FISH and IHC were performed. Unsupervised independent component analysis (ICA) was applied to GEP data and integrated analysis of multi-level molecular data associated with each component (gene signature) was performed. FINDINGS:ICA identified 38 components reflecting transcriptomic variability across our DLBCL cohort. Many of the components were closely related to well-known DLBCL features such as cell-of-origin, stromal and MYC signatures. A component linked to gain of 19q13 locus, among other genomic alterations, was significantly correlated with poor OS and PFS. Through this integrated analysis, a high degree of heterogeneity was highlighted among previously described DLBCL subtypes. INTERPRETATION:The results of this integrated analysis enable a global and multi-level view of DLBCL, as well as improve our understanding of DLBCL subgroups.
| S-EPMC6838437 | biostudies-literature
Project description:BACKGROUND:Gene expression profiling (GEP), next-generation sequencing (NGS) and copy number variation (CNV) analysis have led to an increasingly detailed characterization of the genomic profiles of DLBCL. The aim of this study was to perform a fully integrated analysis of mutational, genomic, and expression profiles to refine DLBCL subtypes. A comparison of our model with two recently published integrative DLBCL classifiers was carried out, in order to best reflect the current state of genomic subtypes. METHODS:223 patients with de novo DLBCL from the prospective, multicenter and randomized LNH-03B LYSA clinical trials were included. GEP data was obtained using Affymetrix GeneChip arrays, mutational profiles were established by Lymphopanel NGS targeting 34 key genes, CNV analysis was obtained by array CGH, and FISH and IHC were performed. Unsupervised independent component analysis (ICA) was applied to GEP data and integrated analysis of multi-level molecular data associated with each component (gene signature) was performed. FINDINGS:ICA identified 38 components reflecting transcriptomic variability across our DLBCL cohort. Many of the components were closely related to well-known DLBCL features such as cell-of-origin, stromal and MYC signatures. A component linked to gain of 19q13 locus, among other genomic alterations, was significantly correlated with poor OS and PFS. Through this integrated analysis, a high degree of heterogeneity was highlighted among previously described DLBCL subtypes. INTERPRETATION:The results of this integrated analysis enable a global and multi-level view of DLBCL, as well as improve our understanding of DLBCL subgroups.
| S-EPMC6838437 | biostudies-literature
Project description:Cancer is a complex and heterogeneous disease involving multiple somatic mutations that accumulate during its progression. In the past years, the wide availability of genomic data from patients' samples opened new perspectives in the analysis of gene mutations and alterations. Hence, visualizing and further identifying genes mutated in massive sets of patients are nowadays a critical task that sheds light on more personalized intervention approaches. Here, we extensively review existing tools for visualization and analysis of alteration data. We compare different approaches to study mutual exclusivity and sample coverage in large-scale omics data. We complement our review with the standalone software AVAtar ('analysis and visualization of alteration data') that integrates diverse aspects known from different tools into a comprehensive platform. AVAtar supplements customizable alteration plots by a multi-objective evolutionary algorithm for subset identification and provides an innovative and user-friendly interface for the evaluation of concurrent solutions. A use case from personalized medicine demonstrates its unique features showing an application on vaccination target selection. AVAtar is available at: https://github.com/sysbio-bioinf/avatar. hans.kestler@uni-ulm.de, phone: +49 (0) 731 500 24 500, fax: +49 (0) 731 500 24 502.
| S-EPMC8138884 | biostudies-literature
Project description:BackgroundSkin melanoma is a highly immunogenic cancer. The intratumoral immune cytolytic activity (CYT) reflects the ability of cytotoxic T and NK cells to eliminate cancer cells, and is associated with improved patient survival. Despite the enthusiastic clinical results seen in advanced-stage metastatic melanoma patients treated with immune checkpoint inhibitors, a subgroup of them will later relapse and develop acquired resistance. We questioned whether CYT associates with different genomic profiles and thus, patient outcome, in skin melanoma.MethodsWe explored the TCGA-SKCM dataset and stratified patients to distinct subgroups of cytolytic activity. The tumor immune contexture, somatic mutations and recurrent copy number aberrations were calculated using quanTIseq, MutSigCV and GISTIC2. Chromothriptic events were explored using CTLPScanner and cancer neoepitopes were predicted with antigen garnish. Each tumor's immunophenoscore was calculated using Immunophenogram. Mutational signatures and kataegis were explored using SigProfiler and compared to the known single or doublet base substitution signatures from COSMIC.ResultsMetastatic skin melanomas had significantly higher CYT levels compared to primary tumors. We assessed enrichment for immune-related gene sets within CYT-high tumors, whereas, CYT-low tumors were enriched for non-immune related gene sets. In addition, distinct mutational and neoantigen loads, primarily composed of C > T transitions, along with specific types of copy number aberrations, characterized each cytolytic subgroup. We found a broader pattern of chromothripsis across CYT-low tumors, where chromosomal regions harboring chromothriptic events, contained a higher number of cancer genes. SBS7a/b, SBS5 and SBS1 were the most prevalent mutational signatures across both cytolytic subgroups, but SBS1 differed significantly between them. SBS7a/b was mutually exclusive with SBS5 and SBS1 in both CYT subgroups. CYT-high patients had markedly higher immunophenoscore, suggesting that they should display a clinical benefit upon treatment with immune checkpoint inhibition therapy, compared to CYT-low patients.ConclusionsOverall, our data highlight the existence of distinct genomic features across cytolytic subgroups in skin melanoma, which might affect the patients' relapse rate or their acquisition of resistance to immune checkpoint inhibition therapies.
| S-EPMC8505325 | biostudies-literature
Project description:Breast cancer (BC) is the most common malignancy with high morbidity and mortality in females worldwide. Emerging evidence indicates that transferrin receptor 1 (TfR1) plays vital roles in regulating cellular iron import. However, the distinct role of TfR1 in BC remains elusive. TfR1 expression was investigated using the TCGA, GEO, TIMER, UALCAN and Oncomine databases. The prognostic potential of TfR1 was evaluated by Kaplan-Meier (KM) plotter and univariate and multivariate Cox regression analyses. Moreover, Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) and gene set enrichment analysis (GSEA) were used to explore the molecular mechanism of TfR1. The potential link between TfR1 expression and infiltrating abundances of immune cells was examined through the TIMER and CIBERSORT algorithm. The expression of TfR1 was dramatically upregulated in BC tissues. Increased TfR1 expression and decreased methylation levels of TfR1 were strongly correlated with multiple clinicopathological parameters. Elevated TfR1 expression was associated with a poor survival rate in BC patients. The nomogram model further confirmed that TfR1 could act as an independent prognostic biomarker in BC. The results of GO, KEGG and GSEA revealed that TfR1 was closely correlated with multiple signaling pathways and immune responses. Additionally, TfR1 was positively associated with the infiltration abundances of six major immune cells, including CD4+ T cells, CD8+ T cells, B cells, neutrophils, macrophages, and dendritic cells in BC. Interestingly, TfR1 influenced prognosis partially through immune infiltration. These comprehensive bioinformatics analyses suggest that TfR1 is a new independent prognostic biomarker and a potential target for immunotherapy in BC.
| S-EPMC8457555 | biostudies-literature