Prognostic value of a lactate metabolism gene signature in lung adenocarcinoma and its associations with immune checkpoint blockade therapy response.

Lung adenocarcinoma (LUAD) is a study that examines the prognostic value of lactate metabolism genes in tumor cells, which are associated with clinical prognosis. We analyzed the expression and clinical data for LUAD from The Cancer Genome Atlas database, using the GSE68465 dataset from the Gene Expression Omnibus and the MSigDB database. LASSO Cox regression and stepwise Cox regression were used to identify the optimal lactate metabolism gene signature. Differences in immune infiltration, tumor mutation burden (TMB), and response to immune checkpoint blockade (ICB) therapy were evaluated between groups. LASSO and Cox regression analyses showed an eight-lactate metabolism gene signature for model construction in both TCGA cohort and GSE68465 data, with higher survival outcomes in high-risk groups. The lactate metabolism risk score had an independent prognostic value (hazard ratio: 2.279 [1.652-3.146], P < .001). Immune cell infiltration differed between the risk groups, such as CD8+ T cells, macrophages, dendritic cells, mast cells, and neutrophils. The high-risk group had higher tumor purity, lower immune and stromal scores, and higher TMB. High-risk samples had high tumor immune dysfunction and exclusion (TIDE) scores and low cytolytic activity (CYT) scores, indicating a poor response to ICB therapy. Similarly, most immune checkpoint molecules, immune inhibitors/stimulators, and major histocompatibility complex (MHC) molecules were highly expressed in the high-risk group. The 8-lactate metabolism gene-based prognostic model predicts patient survival, immune infiltration, and ICB response in patients with LUAD, driving the development of therapeutic strategies to target lactate metabolism.

Huang T ，Lian D ，Chen M ，Liu Y ，Zhang M ，Zeng D ，Zhou SK ，Ying W ... -  《-》

PCDH11X mutation as a potential biomarker for immune checkpoint therapies in lung adenocarcinoma.

Immune checkpoint inhibitors (ICIs) have achieved impressive success in lung adenocarcinoma (LUAD). However, the response to ICIs varies among patients, and predictive biomarkers are urgently needed. PCDH11X is frequently mutated in LUAD, while its role in ICI treatment is unclear. In this study, we curated genomic and clinical data of 151 LUAD patients receiving ICIs from three independent cohorts. Relations between PCDH11X and treatment outcomes of ICIs were examined. A melanoma cohort collected from five published studies, a pan-cancer cohort, and non-ICI-treated TCGA-LUAD cohort were also examined to investigate whether PCDH11X mutation is a specific predictive biomarker for LUAD ICI treatment. Among the three ICI-treated LUAD cohorts, PCDH11X mutation (PCDH11X-MUT) was associated with better clinical response compared to wild-type PCDH11X (PCDH11X-WT). While in ICI-treated melanoma cohort, the pan-cancer cohort excluding LUAD, and the non-ICI-treated TCGA-LUAD cohort, no significant differences in overall survival (OS) were observed between the PCDH11X-MUT and PCDH11X-WT groups. PCDH11X mutation was associated with increased PD-L1 expression, tumor mutation burden (TMB), neoantigen load, DNA damage repair (DDR) mutations, and hot tumor microenvironment in TCGA-LUAD cohort. Our findings suggested that the PCDH11X mutation might serve as a specific biomarker to predict the efficacy of ICIs for LUAD patients. Considering the relatively small sample size of ICI-treated cohorts, future research with larger cohorts and prospective clinical trials will be essential for validating and further exploring the role of PCDH11X mutation in the context of immunotherapy outcomes in LUAD. KEY MESSAGES: PCDH11X mutation is associated with better clinical response compared to wild type PCDH11X in three ICIs-treated LUAD cohorts. In ICIs-treated melanoma cohort, the pan-cancer cohort excluding LUAD, and non-ICIs-treated TCGA-LUAD cohorts PCDH11X mutation is not associated with better clinical response, suggesting PCDH11X mutation might be a specific biomarker to predict the efficacy of ICIs treatment for LUAD patients. PCDH11X mutation is associated with increased PD-L1 expression, tumor mutation burden, and neoantigen load in TCGA-LUAD cohort. PCDH11X mutation is associated with hot tumor microenvironment in TCGA-LUAD cohort.

Liu M ，Yang M ，Zhang B ，Xia S ，Zhao J ，Yan L ，Ren Y ，Guo H ，Zhao J ... -  《-》

Comparison of Two Modern Survival Prediction Tools, SORG-MLA and METSSS, in Patients With Symptomatic Long-bone Metastases Who Underwent Local Treatment With Surgery Followed by Radiotherapy and With Radiotherapy Alone.

Survival estimation for patients with symptomatic skeletal metastases ideally should be made before a type of local treatment has already been determined. Currently available survival prediction tools, however, were generated using data from patients treated either operatively or with local radiation alone, raising concerns about whether they would generalize well to all patients presenting for assessment. The Skeletal Oncology Research Group machine-learning algorithm (SORG-MLA), trained with institution-based data of surgically treated patients, and the Metastases location, Elderly, Tumor primary, Sex, Sickness/comorbidity, and Site of radiotherapy model (METSSS), trained with registry-based data of patients treated with radiotherapy alone, are two of the most recently developed survival prediction models, but they have not been tested on patients whose local treatment strategy is not yet decided. (1) Which of these two survival prediction models performed better in a mixed cohort made up both of patients who received local treatment with surgery followed by radiotherapy and who had radiation alone for symptomatic bone metastases? (2) Which model performed better among patients whose local treatment consisted of only palliative radiotherapy? (3) Are laboratory values used by SORG-MLA, which are not included in METSSS, independently associated with survival after controlling for predictions made by METSSS? Between 2010 and 2018, we provided local treatment for 2113 adult patients with skeletal metastases in the extremities at an urban tertiary referral academic medical center using one of two strategies: (1) surgery followed by postoperative radiotherapy or (2) palliative radiotherapy alone. Every patient's survivorship status was ascertained either by their medical records or the national death registry from the Taiwanese National Health Insurance Administration. After applying a priori designated exclusion criteria, 91% (1920) were analyzed here. Among them, 48% (920) of the patients were female, and the median (IQR) age was 62 years (53 to 70 years). Lung was the most common primary tumor site (41% [782]), and 59% (1128) of patients had other skeletal metastases in addition to the treated lesion(s). In general, the indications for surgery were the presence of a complete pathologic fracture or an impending pathologic fracture, defined as having a Mirels score of ≥ 9, in patients with an American Society of Anesthesiologists (ASA) classification of less than or equal to IV and who were considered fit for surgery. The indications for radiotherapy were relief of pain, local tumor control, prevention of skeletal-related events, and any combination of the above. In all, 84% (1610) of the patients received palliative radiotherapy alone as local treatment for the target lesion(s), and 16% (310) underwent surgery followed by postoperative radiotherapy. Neither METSSS nor SORG-MLA was used at the point of care to aid clinical decision-making during the treatment period. Survival was retrospectively estimated by these two models to test their potential for providing survival probabilities. We first compared SORG to METSSS in the entire population. Then, we repeated the comparison in patients who received local treatment with palliative radiation alone. We assessed model performance by area under the receiver operating characteristic curve (AUROC), calibration analysis, Brier score, and decision curve analysis (DCA). The AUROC measures discrimination, which is the ability to distinguish patients with the event of interest (such as death at a particular time point) from those without. AUROC typically ranges from 0.5 to 1.0, with 0.5 indicating random guessing and 1.0 a perfect prediction, and in general, an AUROC of ≥ 0.7 indicates adequate discrimination for clinical use. Calibration refers to the agreement between the predicted outcomes (in this case, survival probabilities) and the actual outcomes, with a perfect calibration curve having an intercept of 0 and a slope of 1. A positive intercept indicates that the actual survival is generally underestimated by the prediction model, and a negative intercept suggests the opposite (overestimation). When comparing models, an intercept closer to 0 typically indicates better calibration. Calibration can also be summarized as log(O:E), the logarithm scale of the ratio of observed (O) to expected (E) survivors. A log(O:E) > 0 signals an underestimation (the observed survival is greater than the predicted survival); and a log(O:E) < 0 indicates the opposite (the observed survival is lower than the predicted survival). A model with a log(O:E) closer to 0 is generally considered better calibrated. The Brier score is the mean squared difference between the model predictions and the observed outcomes, and it ranges from 0 (best prediction) to 1 (worst prediction). The Brier score captures both discrimination and calibration, and it is considered a measure of overall model performance. In Brier score analysis, the "null model" assigns a predicted probability equal to the prevalence of the outcome and represents a model that adds no new information. A prediction model should achieve a Brier score at least lower than the null-model Brier score to be considered as useful. The DCA was developed as a method to determine whether using a model to inform treatment decisions would do more good than harm. It plots the net benefit of making decisions based on the model's predictions across all possible risk thresholds (or cost-to-benefit ratios) in relation to the two default strategies of treating all or no patients. The care provider can decide on an acceptable risk threshold for the proposed treatment in an individual and assess the corresponding net benefit to determine whether consulting with the model is superior to adopting the default strategies. Finally, we examined whether laboratory data, which were not included in the METSSS model, would have been independently associated with survival after controlling for the METSSS model's predictions by using the multivariable logistic and Cox proportional hazards regression analyses. Between the two models, only SORG-MLA achieved adequate discrimination (an AUROC of > 0.7) in the entire cohort (of patients treated operatively or with radiation alone) and in the subgroup of patients treated with palliative radiotherapy alone. SORG-MLA outperformed METSSS by a wide margin on discrimination, calibration, and Brier score analyses in not only the entire cohort but also the subgroup of patients whose local treatment consisted of radiotherapy alone. In both the entire cohort and the subgroup, DCA demonstrated that SORG-MLA provided more net benefit compared with the two default strategies (of treating all or no patients) and compared with METSSS when risk thresholds ranged from 0.2 to 0.9 at both 90 days and 1 year, indicating that using SORG-MLA as a decision-making aid was beneficial when a patient's individualized risk threshold for opting for treatment was 0.2 to 0.9. Higher albumin, lower alkaline phosphatase, lower calcium, higher hemoglobin, lower international normalized ratio, higher lymphocytes, lower neutrophils, lower neutrophil-to-lymphocyte ratio, lower platelet-to-lymphocyte ratio, higher sodium, and lower white blood cells were independently associated with better 1-year and overall survival after adjusting for the predictions made by METSSS. Based on these discoveries, clinicians might choose to consult SORG-MLA instead of METSSS for survival estimation in patients with long-bone metastases presenting for evaluation of local treatment. Basing a treatment decision on the predictions of SORG-MLA could be beneficial when a patient's individualized risk threshold for opting to undergo a particular treatment strategy ranged from 0.2 to 0.9. Future studies might investigate relevant laboratory items when constructing or refining a survival estimation model because these data demonstrated prognostic value independent of the predictions of the METSSS model, and future studies might also seek to keep these models up to date using data from diverse, contemporary patients undergoing both modern operative and nonoperative treatments. Level III, diagnostic study.

Lee CC ，Chen CW ，Yen HK ，Lin YP ，Lai CY ，Wang JL ，Groot OQ ，Janssen SJ ，Schwab JH ，Hsu FM ，Lin WH ... -  《-》

Machine learning-based prognostic model of lactylation-related genes for predicting prognosis and immune infiltration in patients with lung adenocarcinoma.

Histone lactylation is a novel epigenetic modification that is involved in a variety of critical biological regulations. However, the role of lactylation-related genes in lung adenocarcinoma has yet to be investigated. RNA-seq data and clinical information of LUAD were downloaded from TCGA and GEO datasets. Unsupervised consistent cluster analysis was performed to identify differentially expressed genes (DEGs) between the two clusters, and risk prediction models were constructed by Cox regression analysis and LASSO analysis. Kaplan-Meier (KM) survival analysis, ROC curves and nomograms were used to validate the accuracy of the models. We also explored the differences in risk scores in terms of immune cell infiltration, immune cell function, TMB, TIDE, and anticancer drug sensitivity. In addition, single-cell clustering and trajectory analysis were performed to further understand the significance of lactylation-related genes. We further analyzed lactate content and glucose uptake in lung adenocarcinoma cells and tissues. Changes in LUAD cell function after knockdown of lactate dehydrogenase (LDHA) by CCK-8, colony formation and transwell assays. Finally, we analyzed the expression of KRT81 in LUAD tissues and cell lines using qRT-PCR, WB, and IHC. Changes in KRT81 function in LUAD cells were detected by CCK-8, colony formation, wound healing, transwell, and flow cytometry. A nude mouse xenograft model and a KrasLSL-G12D in situ lung adenocarcinoma mouse model were used to elucidate the role of KRT81 in LUAD. After identifying 26 lactylation-associated DEGs, we constructed 10 lactylation-associated lung adenocarcinoma prognostic models with prognostic value for LUAD patients. A high score indicates a poor prognosis. There were significant differences between the high-risk and low-risk groups in the phenotypes of immune cell infiltration rate, immune cell function, gene mutation frequency, and anticancer drug sensitivity. TMB and TIDE scores were higher in high-risk score patients than in low-risk score patients. MS4A1 was predominantly expressed in B-cell clusters and was identified to play a key role in B-cell differentiation. We further found that lactate content was abnormally elevated in lung adenocarcinoma cells and cancer tissues, and glucose uptake by lung adenocarcinoma cells was significantly increased. Down-regulation of LDHA inhibits tumor cell proliferation, migration and invasion. Finally, we verified that the model gene KRT81 is highly expressed in LUAD tissues and cell lines. Knockdown of KRT81 inhibited cell proliferation, migration, and invasion, leading to cell cycle arrest in the G0/G1 phase and increased apoptosis. KRT81 may play a tumorigenic role in LUAD through the EMT and PI3K/AKT pathways. In vivo, KRT81 knockdown inhibited tumor growth. We successfully constructed a new prognostic model for lactylation-related genes. Lactate content and glucose uptake are significantly higher in lung adenocarcinoma cells and cancer tissues. In addition, KRT81 was validated at cellular and animal levels as a possible new target for the treatment of LUAD, and this study provides a new perspective for the individualized treatment of LUAD.

Gao M ，Wang M ，Zhou S ，Hou J ，He W ，Shu Y ，Wang X ... -  《Cancer Cell International》

Biomarkers of success of anti-PD-(L)1 immunotherapy for non-small cell lung cancer derived from RNA- and whole-exome sequencing: results of a prospective observational study on a cohort of 85 patients.

Immune checkpoint inhibitors (ICIs) treatment have shown high efficacy for about 15 cancer types. However, this therapy is only effective in 20-30% of cancer patients. Thus, the precise biomarkers of ICI response are an urgent need. We conducted a prospective observational study of the prognostic potential ofseveral existing and putative biomarkers of response to immunotherapy in acohort of 85 patients with lung cancer (LC) receiving PD-1 or PD-L1 targeted ICIs. Tumor biosamples were obtained prior to ICI treatment and profiled by whole exome and RNA sequencing. The entire 403 putative biomarkers were screened, including tumor mutation burden (TMB) and number of cancer neoantigens, 131 specific HLA alleles, homozygous state of 11 HLA alleles and their superfamilies; four gene mutation biomarkers, expression of 45 immune checkpoint genes and closely related genes, and three previously published diagnostic gene signatures; for the first time, activation levels of 188 molecular pathways containing immune checkpoint genes and activation levels of 19 pathways algorithmically generated using a human interactome model centered around immune checkpoint genes. Treatment outcomes and/or progression-free survival (PFS) times were available for 61 of 85 patients with LC, including 24 patients with adenocarcinoma and 27 patients with squamous cell LC, whose samples were further analyzed. For the rest 24 patients, both treatment outcomes and PFS data could not be collected. Of these, 54 patients were treated with PD1-specific and 7 patients with PD-L1-specific ICIs. We evaluated the potential of biomarkers based on PFS and RECIST treatment response data. In our sample, 45 biomarkers were statistically significantly associated with PFS and 44 with response to treatment, of which eight were shared. Five of these (CD3G and NCAM1 gene expression levels, and levels of activation of Adrenergic signaling in cardiomyocytes, Growth hormone signaling, and Endothelin molecular pathways) were used in our signature that showed an AUC of 0.73 and HR of 0.27 (p=0.00034) on the experimental dataset. This signature was also reliable (AUC 0.76, 0.87) for the independent publicly available LC datasets GSE207422, GSE126044 annotated with ICI response data and demonstrated same survival trends on independent dataset GSE135222 annotated with PFS data. In both experimental and one independent datasets annotated with samples' histotypes, the signature worked better for lung adenocarcinoma than for squamous cell LC. The high reliability of our signature to predict response and PFS after ICI treatment was demonstrated using experimental and 3 independent datasets. Additionally, annotated molecular profiles obtained in this study were made publicly accessible.

Poddubskaya E ，Suntsova M ，Lyadova M ，Luppov D ，Guryanova A ，Lyadov V ，Garazha A ，Sorokin M ，Semenova A ，Shatalov V ，Biakhova M ，Simonov A ，Moisseev A ，Buzdin A ... -  《Frontiers in Immunology》