B Cells and T Follicular Helper Cells Mediate Response to Checkpoint Inhibitors in High Mutation Burden Mouse Models of Breast Cancer.

This study identifies mechanisms mediating responses to immune checkpoint inhibitors using mouse models of triple-negative breast cancer. By creating new mammary tumor models, we find that tumor mutation burden and specific immune cells are associated with response. Further, we developed a rich resource of single-cell RNA-seq and bulk mRNA-seq data of immunotherapy-treated and non-treated tumors from sensitive and resistant murine models. Using this, we uncover that immune checkpoint therapy induces T follicular helper cell activation of B cells to facilitate the anti-tumor response in these models. We also show that B cell activation of T cells and the generation of antibody are key to immunotherapy response and propose a new biomarker for immune checkpoint therapy. In total, this work presents resources of new preclinical models of breast cancer with large mRNA-seq and single-cell RNA-seq datasets annotated for sensitivity to therapy and uncovers new components of response to immune checkpoint inhibitors.

A Pan-Cancer Analysis of Enhancer Expression in Nearly 9000 Patient Samples.

The role of enhancers, a key class of non-coding regulatory DNA elements, in cancer development has increasingly been appreciated. Here, we present the detection and characterization of a large number of expressed enhancers in a genome-wide analysis of 8928 tumor samples across 33 cancer types using TCGA RNA-seq data. Compared with matched normal tissues, global enhancer activation was observed in most cancers. Across cancer types, global enhancer activity was positively associated with aneuploidy, but not mutation load, suggesting a hypothesis centered on "chromatin-state" to explain their interplay. Integrating eQTL, mRNA co-expression, and Hi-C data analysis, we developed a computational method to infer causal enhancer-gene interactions, revealing enhancers of clinically actionable genes. Having identified an enhancer ∼140 kb downstream of PD-L1, a major immunotherapy target, we validated it experimentally. This study provides a systematic view of enhancer activity in diverse tumor contexts and suggests the clinical implications of enhancers.

The ectonucleotidase CD39 identifies tumor-reactive CD8+ T cells predictive of immune checkpoint blockade efficacy in human lung cancer.

Improved identification of anti-tumor T cells is needed to advance cancer immunotherapies. CD39 expression is a promising surrogate of tumor-reactive CD8+ T cells. Here, we comprehensively profiled CD39 expression in human lung cancer. CD39 expression enriched for CD8+ T cells with features of exhaustion, tumor reactivity, and clonal expansion. Flow cytometry of 440 lung cancer biospecimens revealed weak association between CD39+ CD8+ T cells and tumoral features, such as programmed death-ligand 1 (PD-L1), tumor mutation burden, and driver mutations. Immune checkpoint blockade (ICB), but not cytotoxic chemotherapy, increased intratumoral CD39+ CD8+ T cells. Higher baseline frequency of CD39+ CD8+ T cells conferred improved clinical outcomes from ICB therapy. Furthermore, a gene signature of CD39+ CD8+ T cells predicted benefit from ICB, but not chemotherapy, in a phase III clinical trial of non-small cell lung cancer. These findings highlight CD39 as a proxy of tumor-reactive CD8+ T cells in human lung cancer.

Tumor and Microenvironment Evolution during Immunotherapy with Nivolumab.

The mechanisms by which immune checkpoint blockade modulates tumor evolution during therapy are unclear. We assessed genomic changes in tumors from 68 patients with advanced melanoma, who progressed on ipilimumab or were ipilimumab-naive, before and after nivolumab initiation (CA209-038 study). Tumors were analyzed by whole-exome, transcriptome, and/or T cell receptor (TCR) sequencing. In responding patients, mutation and neoantigen load were reduced from baseline, and analysis of intratumoral heterogeneity during therapy demonstrated differential clonal evolution within tumors and putative selection against neoantigenic mutations on-therapy. Transcriptome analyses before and during nivolumab therapy revealed increases in distinct immune cell subsets, activation of specific transcriptional networks, and upregulation of immune checkpoint genes that were more pronounced in patients with response. Temporal changes in intratumoral TCR repertoire revealed expansion of T cell clones in the setting of neoantigen loss. Comprehensive genomic profiling data in this study provide insight into nivolumab's mechanism of action.

BRN2 suppresses apoptosis, reprograms DNA damage repair, and is associated with a high somatic mutation burden in melanoma.

Whether cell types exposed to a high level of environmental insults possess cell type-specific prosurvival mechanisms or enhanced DNA damage repair capacity is not well understood. BRN2 is a tissue-restricted POU domain transcription factor implicated in neural development and several cancers. In melanoma, BRN2 plays a key role in promoting invasion and regulating proliferation. Here we found, surprisingly, that rather than interacting with transcription cofactors, BRN2 is instead associated with DNA damage response proteins and directly binds PARP1 and Ku70/Ku80. Rapid PARP1-dependent BRN2 association with sites of DNA damage facilitates recruitment of Ku80 and reprograms DNA damage repair by promoting Ku-dependent nonhomologous end-joining (NHEJ) at the expense of homologous recombination. BRN2 also suppresses an apoptosis-associated gene expression program to protect against UVB-, chemotherapy- and vemurafenib-induced apoptosis. Remarkably, BRN2 expression also correlates with a high single-nucleotide variation prevalence in human melanomas. By promoting error-prone DNA damage repair via NHEJ and suppressing apoptosis of damaged cells, our results suggest that BRN2 contributes to the generation of melanomas with a high mutation burden. Our findings highlight a novel role for a key transcription factor in reprogramming DNA damage repair and suggest that BRN2 may impact the response to DNA-damaging agents in BRN2-expressing cancers.

Construction of a prognostic model for breast cancer based on moonlighting genes.

Breast cancer (BRCA) is a highly heterogeneous disease, with significant differences in prognosis among patients. Existing biomarkers and prognostic models have limited ability to predict BRCA prognosis. Moonlighting genes regulate tumor progression and are associated with cancer prognosis. This study aimed to construct a moonlighting gene-based prognostic model for BRCA. We obtained differentially expressed genes (DEGs) in BRCA from The Cancer Genome Atlas and intersected them with moonlighting genes from MoonProt to acquire differential moonlighting genes. GO and KEGG results showed main enrichment of these genes in the response of BRCA cells to environmental stimuli and pentose phosphate pathway. Based on moonlighting genes, we conducted drug prediction and validated results through cellular experiments. After ABCB1 knockdown, viability and proliferation of BRCA cells were significantly enhanced. Based on differential moonlighting genes, BRCA was divided into three subgroups, among which cluster2 had the highest survival rate and immunophenoscore and relatively low tumor mutation burden. TP53 had the highest mutation frequency in cluster2 and cluster3, while PIK3CA had a higher mutation frequency in cluster1, with the majority being missense mutations. Subsequently, we established an 11-gene prognostic model in the training set based on DEGs among subgroups using univariate Cox regression, LASSO regression, and multivariable Cox regression analyses. Model prognostic performance was verified in GEO, METABRIC and ICGC validation sets. In summary, this study obtained three BRCA moonlighting gene-related subtypes and constructed an 11-gene prognostic model. The 11-gene BRCA prognostic model has good predictive performance, guiding BRCA prognosis for clinical doctors.

TMBstable: a variant caller controls performance variation across heterogeneous sequencing samples.

In cancer genomics, variant calling has advanced, but traditional mean accuracy evaluations are inadequate for biomarkers like tumor mutation burden, which vary significantly across samples, affecting immunotherapy patient selection and threshold settings. In this study, we introduce TMBstable, an innovative method that dynamically selects optimal variant calling strategies for specific genomic regions using a meta-learning framework, distinguishing it from traditional callers with uniform sample-wide strategies. The process begins with segmenting the sample into windows and extracting meta-features for clustering, followed by using a pre-trained meta-model to select suitable algorithms for each cluster, thereby addressing strategy-sample mismatches, reducing performance fluctuations and ensuring consistent performance across various samples. We evaluated TMBstable using both simulated and real non-small cell lung cancer and nasopharyngeal carcinoma samples, comparing it with advanced callers. The assessment, focusing on stability measures, such as the variance and coefficient of variation in false positive rate, false negative rate, precision and recall, involved 300 simulated and 106 real tumor samples. Benchmark results showed TMBstable's superior stability with the lowest variance and coefficient of variation across performance metrics, highlighting its effectiveness in analyzing the counting-based biomarker. The TMBstable algorithm can be accessed at https://github.com/hello-json/TMBstable for academic usage only.

Defining the single base importance of human mRNAs and lncRNAs.

As the fundamental unit of a gene and its transcripts, nucleotides have enormous impacts on the gene function and evolution, and thus on phenotypes and diseases. In order to identify the key nucleotides of one specific gene, it is quite crucial to quantitatively measure the importance of each base on the gene. However, there are still no sequence-based methods of doing that. Here, we proposed Base Importance Calculator (BIC), an algorithm to calculate the importance score of each single base based on sequence information of human mRNAs and long noncoding RNAs (lncRNAs). We then confirmed its power by applying BIC to three different tasks. Firstly, we revealed that BIC can effectively evaluate the pathogenicity of both genes and single bases through single nucleotide variations. Moreover, the BIC score in The Cancer Genome Atlas somatic mutations is able to predict the prognosis of some cancers. Finally, we show that BIC can also precisely predict the transmissibility of SARS-CoV-2. The above results indicate that BIC is a useful tool for evaluating the single base importance of human mRNAs and lncRNAs.

Comprehensive bioinformatics analysis of human cytomegalovirus pathway genes in pan-cancer.

BACKGROUND: Human cytomegalovirus (HCMV) is a herpesvirus that can infect various cell types and modulate host gene expression and immune response. It has been associated with the pathogenesis of various cancers, but its molecular mechanisms remain elusive. METHODS: We comprehensively analyzed the expression of HCMV pathway genes across 26 cancer types using the Cancer Genome Atlas (TCGA) and The Genotype-Tissue Expression (GTEx) databases. We also used bioinformatics tools to study immune invasion and tumor microenvironment in pan-cancer. Cox regression and machine learning were used to analyze prognostic genes and their relationship with drug sensitivity. RESULTS: We found that HCMV pathway genes are widely expressed in various cancers. Immune infiltration and the tumor microenvironment revealed that HCMV is involved in complex immune processes. We obtained prognostic genes for 25 cancers and significantly found 23 key genes in the HCMV pathway, which are significantly enriched in cellular chemotaxis and synaptic function and may be involved in disease progression. Notably, CaM family genes were up-regulated and AC family genes were down-regulated in most tumors. These hub genes correlate with sensitivity or resistance to various drugs, suggesting their potential as therapeutic targets. CONCLUSIONS: Our study has revealed the role of the HCMV pathway in various cancers and provided insights into its molecular mechanism and therapeutic significance. It is worth noting that the key genes of the HCMV pathway may open up new doors for cancer prevention and treatment.

Transcription-independent functions of p53 in DNA repair pathway selection.

Recently discovered transcription-independent features of p53 involve the choice of DNA damage repair pathway after PARylation, and p53's complex formation with phosphoinositide lipids, PI(4,5)P2 . PARylation-mediated rapid accumulation of p53 at DNA damage sites is linked to the recruitment of downstream repair factors and tumor suppression. This links p53's capability to sense damaged DNA in vitro and its relevant functions in cells. Further, PI(4,5)P2 rapidly accumulates at damage sites like p53 and complexes with p53, while it is required for ATR recruitment. These findings help explain how p53 and PI(4,5)P2 maintain genome stability by directing DNA repair pathway choice. Additionally, there is a strong correlation between p53 sequence homology, genome mutation rates as well as lifespans across various mammalian species. Further investigation is required to better understand the connections between genome stability, tumor suppression, longevity and the transcriptional-independent function of p53.

Upregulation of POC1A in lung adenocarcinoma promotes tumour progression and predicts poor prognosis.

Lung adenocarcinoma (LUAD) is characterized by a high incidence rate and mortality. Recently, POC1 centriolar protein A (POC1A) has emerged as a potential biomarker for various cancers, contributing to cancer onset and development. However, the association between POC1A and LUAD remains unexplored. We extracted The Cancer Genome Atlas (TCGA) and the Gene Expression Omnibus (GEO) data sets to analyse the differential expression of POC1A and its relationship with clinical stage. Additionally, we performed diagnostic receiver operator characteristic (ROC) curve analysis and Kaplan-Meier (KM) survival analysis to assess the diagnostic and prognostic value of POC1A in LUAD. Furthermore, we investigated the correlation between POC1A expression and immune infiltration, tumour mutation burden (TMB), immune checkpoint expression and drug sensitivity. Finally, we verified POC1A expression using real-time quantitative polymerase chain reaction (RT-qPCR) and immunohistochemistry (IHC). Cell experiments were conducted to validate the effect of POC1A expression on the proliferation, migration and invasion of lung cancer cells. POC1A exhibited overexpression in most tumour tissues, and its overexpression in LUAD was significantly correlated with late-stage presentation and poor prognosis. The high POC1A expression group showed lower levels of immune infiltration but higher levels of immune checkpoint expression and TMB. Moreover, the high POC1A expression group demonstrated sensitivity to multiple drugs. In vitro experiments confirmed that POC1A knockdown led to decreased proliferation, migration, and invasion of lung cancer cells. Our findings suggest that POC1A may contribute to tumour development by modulating the cell cycle and immune cell infiltration. It also represents a potential therapeutic target and marker for the diagnosis and prognosis of LUAD.

Prognostic biomarkers associated with immune checkpoint inhibitors in hepatocellular carcinoma.

The treatment of hepatocellular carcinoma (HCC), particularly advanced HCC, has been a serious challenge. Immune checkpoint inhibitors (ICIs) are landmark drugs in the field of cancer therapy in recent years, which have changed the landscape of cancer treatment. In the field of HCC treatment, this class of drugs has shown good therapeutic prospects. For example, atezolizumab in combination with bevacizumab has been approved as first-line treatment for advanced HCC due to significant efficacy. However, sensitivity to ICI therapy varies widely among HCC patients. Therefore, there is an urgent need to search for determinants of resistance/sensitivity to ICIs and to screen biomarkers that can predict the efficacy of ICIs. This manuscript reviews the research progress of prognostic biomarkers associated with ICIs in HCC in order to provide a scientific basis for the development of clinically individualised precision medication regimens.

Deleterious and Adaptive Mutations in Plant Germplasm Conserved Ex Situ.

Conserving more than 7 million plant germplasm accessions in 1,750 genebanks worldwide raises the hope of securing the food supply for humanity for future generations. However, there is a genetic cost for such long-term germplasm conservation, which has been largely unaccounted for before. We investigated the extent and variation of deleterious and adaptive mutations in 490 individual plants representing barley, wheat, oat, soybean, maize, rapa, and sunflower collections in a seed genebank using RNA-Seq technology. These collections were found to have a range of deleterious mutations detected from 125 (maize) to 83,695 (oat) with a mean of 13,537 and of the averaged sample-wise mutation burden per deleterious locus from 0.069 to 0.357 with a mean of 0.200. Soybean and sunflower collections showed that accessions acquired earlier had increased mutation burdens. The germplasm with more years of storage in several collections carried more deleterious and fewer adaptive mutations. The samples with more cycles of germplasm regeneration revealed fewer deleterious and more adaptive mutations. These findings are significant for understanding mutational dynamics and genetic cost in conserved germplasm and have implications for long-term germplasm management and conservation.

Automated Quantitative CD8+ Tumor-Infiltrating Lymphocytes and Tumor Mutation Burden as Independent Biomarkers in Melanoma Patients Receiving Front-Line Anti-PD-1 Immunotherapy.

BACKGROUND: CD8+ tumor-infiltrating lymphocyte (TIL) predicts response to anti-PD-(L)1 therapy. However, there remains no standardized method to assess CD8+ TIL in melanoma, and developing a specific, cost-effective, reproducible, and clinically actionable biomarker to anti-PD-(L)1 remains elusive. We report on the development of automatic CD8+ TIL density quantification via whole slide image (WSI) analysis in advanced melanoma patients treated with front-line anti-PD-1 blockade, and correlation immunotherapy response. METHODS: Seventy-eight patients treated with PD-1 inhibitors in the front-line setting between January 2015 and May 2023 at the University of Pittsburgh Cancer Institute were included. CD8+ TIL density was quantified using an image analysis algorithm on digitized WSI. Targeted next-generation sequencing (NGS) was performed to determine tumor mutation burden (TMB) in a subset of 62 patients. ROC curves were used to determine biomarker cutoffs and response to therapy. Correlation between CD8+ TIL density and TMB cutoffs and response to therapy was studied. RESULTS: Higher CD8+ TIL density was significantly associated with improved response to front-line anti-PD-1 across all time points measured. CD8+ TIL density ≥222.9 cells/mm2 reliably segregated responders and non-responders to front-line anti-PD-1 therapy regardless of when response was measured. In a multivariate analysis, patients with CD8+ TIL density exceeding cutoff had significantly improved PFS with a trend toward improved OS. Similarly, increasing TMB was associated with improved response to anti-PD-1, and a cutoff of 14.70 Mut/Mb was associated with improved odds of response. The correlation between TMB and CD8+ TIL density was low, suggesting that each represented independent predictive biomarkers of response. CONCLUSIONS: An automatic digital analysis algorithm provides a standardized method to quantify CD8+ TIL density, which predicts response to front-line anti-PD-1 therapy. CD8+ TIL density and TMB are independent predictors of response to anti-PD-1 blockade.

Constitutional Microsatellite Instability, Genotype, and Phenotype Correlations in Constitutional Mismatch Repair Deficiency.

BACKGROUND & AIMS: Constitutional mismatch repair deficiency (CMMRD) is a rare recessive childhood cancer predisposition syndrome caused by germline mismatch repair variants. Constitutional microsatellite instability (cMSI) is a CMMRD diagnostic hallmark and may associate with cancer risk. We quantified cMSI in a large CMMRD patient cohort to explore genotype-phenotype correlations using novel MSI markers selected for instability in blood. METHODS: Three CMMRD, 1 Lynch syndrome, and 2 control blood samples were genome sequenced to >120× depth. A pilot cohort of 8 CMMRD and 38 control blood samples and a blinded cohort of 56 CMMRD, 8 suspected CMMRD, 40 Lynch syndrome, and 43 control blood samples were amplicon sequenced to 5000× depth. Sample cMSI score was calculated using a published method comparing microsatellite reference allele frequencies with 80 controls. RESULTS: Thirty-two mononucleotide repeats were selected from blood genome and pilot amplicon sequencing data. cMSI scoring using these MSI markers achieved 100% sensitivity (95% CI, 93.6%-100.0%) and specificity (95% CI 97.9%-100.0%), was reproducible, and was superior to an established tumor MSI marker panel. Lower cMSI scores were found in patients with CMMRD with MSH6 deficiency and patients with at least 1 mismatch repair missense variant, and patients with biallelic truncating/copy number variants had higher scores. cMSI score did not correlate with age at first tumor. CONCLUSIONS: We present an inexpensive and scalable cMSI assay that enhances CMMRD detection relative to existing methods. cMSI score is associated with mismatch repair genotype but not phenotype, suggesting it is not a useful predictor of cancer risk.

Homologous recombination deficiency reflects the heterogeneity and monitoring treatment response for patients with breast cancer.

BACKGROUND: In breast cancer (BC), homologous recombination defect (HRD) is a common carcinogenic mechanism. It is meaningful to classify BC according to HRD biomarkers and to develop a platform for identifying BC molecular features, pathological features and therapeutic responses. METHODS: In total, 109 HRD genes were collected and screened by univariate Cox regression analysis to determine the prognostic genes, which were used to construct a consensus matrix to identify BC subtype. Differentially expressed genes (DEGs) were filtered by the Limma package and screened by random forest analysis to build a model to analyze the immunotherapy response and sensitivity and prognosis of patients suffering from BC to different drugs. RESULTS: Thirteen out of 109 HRD genes were prognostic genes of BC, and BC was classified into two subgroups based on their expression. Cluster 1 had a significantly backward survival outcome and a significantly higher adaptive immunity score relative to cluster 2. Six genes were identified by random forest analysis as factors for developing the model. The model provided a prediction called risk score, which showed a significant stratification effect on BC prognosis, immunotherapy response and IC50 values of 62 drugs. CONCLUSIONS: In the present study, two HRD subtypes of BC were successfully identified, for which mutation and immunological features were determined. A model based on differential genes of HRD subtypes was established, which was a potential predictor of prognosis, immunotherapy response and drug sensitivity of BC.

Genomic characterization and immunotherapy for microsatellite instability-high in cholangiocarcinoma.

BACKGROUND: Microsatellite instability-high (MSI-H) is a unique genomic status in many cancers. However, its role in the genomic features and immunotherapy in cholangiocarcinoma (CCA) is unclear. This study aimed to systematically investigate the genomic characterization and immunotherapy efficacy of MSI-H patients with CCA. METHODS: We enrolled 887 patients with CCA in this study. Tumor samples were collected for next-generation sequencing. Differences in genomic alterations between the MSI-H and microsatellite stability (MSS) groups were analyzed. We also investigated the survival of PD-1 inhibitor-based immunotherapy between two groups of 139 patients with advanced CCA. RESULTS: Differential genetic alterations between the MSI-H and MSS groups included mutations in ARID1A, ACVR2A, TGFBR2, KMT2D, RNF43, and PBRM1 which were enriched in MSI-H groups. Patients with an MSI-H status have a significantly higher tumor mutation burden (TMB) (median 41.7 vs. 3.1 muts/Mb, P < 0.001) and more positive programmed death ligand 1 (PD-L1) expression (37.5% vs. 11.9%, P < 0.001) than those with an MSS status. Among patients receiving PD-1 inhibitor-based therapy, those with MSI-H had a longer median overall survival (OS, hazard ratio (HR) = 0.17, P = 0.001) and progression-free survival (PFS, HR = 0.14, P < 0.001) than patients with MSS. Integrating MSI-H and PD-L1 expression status (combined positive score ≥ 5) could distinguish the efficacy of immunotherapy. CONCLUSIONS: MSI-H status was associated with a higher TMB value and more positive PD-L1 expression in CCA tumors. Moreover, in patients with advanced CCA who received PD-1 inhibitor-based immunotherapy, MSI-H and positive PD-L1 expression were associated with improved both OS and PFS. TRIAL REGISTRATION: This study was registered on ClinicalTrials.gov on 07/01/2017 (NCT03892577).

Concomitant genomic features stratify prognosis to patients with advanced EGFR mutant lung cancer.

This study aimed to explore the clinical significance of genomics features including tumor mutation burden (TMB) and copy number alteration (CNA) for advanced EGFR mutant lung cancer. We retrospectively identified 1378 patients with advanced EGFR mutant lung cancer and next-generation sequencing tests from three cohorts. Multiple co-occurring genomics alternations occurred in a large proportion (97%) of patients with advanced EGFR mutant lung cancers. Both TMB and CNA were predictive biomarkers for these patients. A joint analysis of TMB and CNA found that patients with high TMB and high CNA showed worse responses to EGFR-TKIs and predicted worse outcomes. TMBhighCNAhigh, as a high-risk genomic feature, showed predictive ability in most of the subgroups based on clinical characteristics. These patients had larger numbers of metastatic sites, and higher rates of EGFR copy number amplification, TP53 mutations, and cell-cycle gene alterations, which showed more potential survival gain from combination treatment. Furthermore, a nomogram based on genomic features and clinical features was developed to distinguish prognosis. Genomic features could stratify prognosis and guide clinical treatment for patients with advanced EGFR mutant lung cancer.

A multi-dimensional approach to unravel the intricacies of lactylation related signature for prognostic and therapeutic insight in colorectal cancer.

BACKGROUND: Lactylation, a novel contributor to post-translational protein modifications, exhibits dysregulation across various tumors. Nevertheless, its intricate involvement in colorectal carcinoma, particularly for non-histone lactylation and its intersection with metabolism and immune evasion, remains enigmatic. METHODS: Employing immunohistochemistry on tissue microarray with clinical information and immunofluorescence on colorectal cell lines, we investigated the presence of global lactylation and its association with development and progression in colorectal cancer as well as its functional location. Leveraging the AUCell algorithm alongside correlation analysis in single-cell RNA sequencing data, as well as cox-regression and lasso-regression analysis in TCGA dataset and confirmed in GEO dataset, we identified a 23-gene signature predicting colorectal cancer prognosis. Subsequently, we analyzed the associations between the lactylation related gene risk and clinical characteristics, mutation landscapes, biological functions, immune cell infiltration, immunotherapy responses, and drug sensitivity. Core genes were further explored for deep biological insights through bioinformatics and in vitro experiments. RESULTS: Our study innovatively reveals a significant elevation of global lactylation in colorectal cancer, particularly in malignant tumors, confirming it as an independent prognostic factor for CRC. Through a comprehensive analysis integrating tumor tissue arrays, TCGA dataset, GEO dataset, combining in silico investigations and in vitro experiments, we identified a 23-gene Lactylation-Related Gene risk model capable of predicting the prognosis of colorectal cancer patients. Noteworthy variations were observed in clinical characteristics, biological functions, immune cell infiltration, immune checkpoint expression, immunotherapy responses and drug sensitivity among distinct risk groups. CONCLUSIONS: The Lactylation-Related Gene risk model exhibits significant potential for improving the management of colorectal cancer patients and enhancing therapeutic outcomes, particularly at the intersection of metabolism and immune evasion. This finding underscores the clinical relevance of global lactylation in CRC and lays the groundwork for mechanism investigation and targeted therapeutic strategies given the high lactate concentration in CRC.

The Landscape and Prognosis of Microsatellite Stable (MSS) Esophageal, Gastro-Esophageal Junction and Gastric Adenocarcinomas with High Tumor Mutation Burden (TMB).

BACKGROUND: A minority of patients with MSS tumors present a high tumor mutation burden (TMB) without underlying MMR defects. METHODS: Publicly available genomic series were assessed for identification of patients with MSS gastric gastroesophageal junction, and esophageal adenocarcinomas and a high TMB, defined as more than 10 mutations per Mb. These were compared with MSS cancers and a low TMB for genetic alterations and for survival outcomes. RESULTS: Patients with MSS cancers with high TMB in the MSK series were older but did not differ in other clinicopathologic parameters compared with MSS patients with low TMB. Mutations in tumor suppressors TP53 and APC and oncogenes KRAS and ERBB4 as well as amplifications of ERBB2 were more prevalent in the high TMB group of MSS cancers. Mutations in DDR associated genes, in epigenetic modifiers and in genes associated with immune response were more prevalent in the hIgh TMB group patients. However, high TMB was not associated with an improved survival in MSS gastric/gastroesophageal junction/esophageal adenocarcinomas (Log Rank p = 0.5). CONCLUSION: MSS Gastric/gastroesophageal junction/esophageal adenocarcinomas with TMB above 10 mutations per Mb possess a genomic landscape with increased alteration frequencies in common gastroesophageal cancer genes and pathways.