TMBcat: A multi-endpoint p-value criterion on different discrepancy metrics for superiorly inferring tumor mutation burden thresholds

DeepLION2: deep multi-instance contrastive learning framework enhancing the prediction of cancer-associated T cell receptors by attention strategy on motifs

Introduction

T cell receptor (TCR) repertoires provide valuable insights into complex human diseases, including cancers. Recent advancements in immune sequencing technology have significantly improved our understanding of TCR repertoire. Some computational methods have been devised to identify cancer-associated TCRs and enable cancer detection using TCR sequencing data. However, the existing methods are often limited by their inadequate consideration of the correlations among TCRs within a repertoire, hindering the identification of crucial TCRs. Additionally, the sparsity of cancer-associated TCR distribution presents a challenge in accurate prediction.

Methods

To address these issues, we presented DeepLION2, an innovative deep multi-instance contrastive learning framework specifically designed to enhance cancer-associated TCR prediction. DeepLION2 leveraged content-based sparse self-attention, focusing on the top k related TCRs for each TCR, to effectively model inter-TCR correlations. Furthermore, it adopted a contrastive learning strategy for bootstrapping parameter updates of the attention matrix, preventing the model from fixating on non-cancer-associated TCRs.

Results

Extensive experimentation on diverse patient cohorts, encompassing over ten cancer types, demonstrated that DeepLION2 significantly outperformed current state-of-the-art methods in terms of accuracy, sensitivity, specificity, Matthews correlation coefficient, and area under the curve (AUC). Notably, DeepLION2 achieved impressive AUC values of 0.933, 0.880, and 0.763 on thyroid, lung, and gastrointestinal cancer cohorts, respectively. Furthermore, it effectively identified cancer-associated TCRs along with their key motifs, highlighting the amino acids that play a crucial role in TCR-peptide binding.

Conclusion

These compelling results underscore DeepLION2's potential for enhancing cancer detection and facilitating personalized cancer immunotherapy. DeepLION2 is publicly available on GitHub, at https://github.com/Bioinformatics7181/DeepLION2, for academic use only.

Molecular subtype identification and prognosis stratification based on golgi apparatus-related genes in head and neck squamous cell carcinoma

BACKGROUND

Abnormal dynamics of the Golgi apparatus reshape the tumor microenvironment and immune landscape, playing a crucial role in the prognosis and treatment response of cancer. This study aims to investigate the potential role of Golgi apparatus-related genes (GARGs) in the heterogeneity and prognosis of head and neck squamous cell carcinoma (HNSCC).

METHODS

Transcriptional data and corresponding clinical information of HNSCC were obtained from public databases for differential expression analysis, consensus clustering, survival analysis, immune infiltration analysis, immune therapy response assessment, gene set enrichment analysis, and drug sensitivity analysis. Multiple machine learning algorithms were employed to construct a prognostic model based on GARGs. A nomogram was used to integrate and visualize the multi-gene model with clinical pathological features.

RESULTS

A total of 321 GARGs that were differentially expressed were identified, out of which 69 were associated with the prognosis of HNSCC. Based on these prognostic genes, two molecular subtypes of HNSCC were identified, which showed significant differences in prognosis. Additionally, a risk signature consisting of 28 GARGs was constructed and demonstrated good performance for assessing the prognosis of HNSCC. This signature divided HNSCC into the high-risk and low-risk groups with significant differences in multiple clinicopathological characteristics, including survival outcome, grade, T stage, chemotherapy. Immune response-related pathways were significantly activated in the high-risk group with better prognosis. There were significant differences in chemotherapy drug sensitivity and immune therapy response between the high-risk and low-risk groups, with the low-risk group being more suitable for receiving immunotherapy. Riskscore, age, grade, and radiotherapy were independent prognostic factors for HNSCC and were used to construct a nomogram, which had good clinical applicability.

CONCLUSIONS

We successfully identified molecular subtypes and prognostic signature of HNSCC that are derived from GARGs, which can be used for the assessment of HNSCC prognosis and treatment responses.

Development of a risk assessment model for cardiac injury in patients newly diagnosed with acute myeloid leukemia based on a multicenter, real-world analysis in China

BACKGROUND

Studies have revealed that acute myeloid leukemia (AML) patients are prone to combined cardiac injury. We aimed to identify hematological risk factors associated with cardiac injury in newly diagnosed AML patients before chemotherapy and develop a personalized predictive model.

METHODS

The population baseline, blood test, electrocardiogram, echocardiograph, and genetic and cytogenetic data were collected from newly diagnosed AML patients. The data were subdivided into training and validation cohorts. The independent risk factors were explored by univariate and multivariate logistic regression analysis respectively, and data dimension reduction and variable selection were performed using the least absolute shrinkage and selection operator (LASSO) regression models. The nomogram was generated and the reliability and generalizability were verified by receiver operating characteristic (ROC) curves, the area under the curve (AUC) and calibration curves in an external validation cohort.

RESULTS

Finally, 499 AML patients were included. After univariate logistic regression, LASSO regression and multivariate logistic regression analysis, abnormal NT-proBNP, NPM1 mutation, WBC, and RBC were independent risk factors for cardiac injury in AML patients (all P < 0.05). The nomogram was constructed based on the above four variables with high accuracy. The area under the curve was 0.742, 0.750, and 0.706 in the training, internal validation, and external validation cohort, respectively. The calibration curve indicated that the model has good testing capability. The Kaplan-Meier curve showed that the higher the risk of combined cardiac injury in AML patients, the lower their probability of survival.

CONCLUSIONS

This prediction nomogram identifies hematological risk factors associated with cardiac injury in newly diagnosed AML patients and can help hematologists identify the risk and provide precise treatment options.

What makes TMB an ambivalent biomarker for immunotherapy? A subtle mismatch between the sample-based design of variant callers and real clinical cohort

Tumor mutation burden (TMB) is a widely recognized biomarker for predicting the efficacy of immunotherapy. However, its use still remains highly controversial. In this study, we examine the underlying causes of this controversy based on clinical needs. By tracing the source of the TMB errors and analyzing the design philosophy behind variant callers, we identify the conflict between the incompleteness of biostatistics rules and the variety of clinical samples as the critical issue that renders TMB an ambivalent biomarker. A series of experiments were conducted to illustrate the challenges of mutation detection in clinical practice. Additionally, we also discuss potential strategies for overcoming these conflict issues to enable the application of TMB in guiding decision-making in real clinical settings.

TMBserval: a statistical explainable learning model reveals weighted tumor mutation burden better categorizing therapeutic benefits

A high tumor mutation burden (TMB) is known to drive the response to immune checkpoint inhibitors (ICI) and is associated with favorable prognoses. However, because it is a one-dimensional numerical representation of non-synonymous genetic alterations, TMB suffers from clinical challenges due to its equal quantification. Since not all mutations elicit the same antitumor rejection, the effect on immunity of neoantigens encoded by different types or locations of somatic mutations may vary. In addition, other typical genomic features, including complex structural variants, are not captured by the conventional TMB metric. Given the diversity of cancer subtypes and the complexity of treatment regimens, this paper proposes that tumor mutations capable of causing various degrees of immunogenicity should be calculated separately. TMB should therefore, be segmented into more exact, higher dimensional feature vectors to exhaustively measure the foreignness of tumors. We systematically reviewed patients' multifaceted efficacy based on a refined TMB metric, investigated the association between multidimensional mutations and integrative immunotherapy outcomes, and developed a convergent categorical decision-making framework, TMBserval (Statistical Explainable machine learning with Regression-based VALidation). TMBserval integrates a multiple-instance learning concept with statistics to create a statistically interpretable model that addresses the broad interdependencies between multidimensional mutation burdens and decision endpoints. TMBserval is a pan-cancer-oriented many-to-many nonlinear regression model with discrimination and calibration power. Simulations and experimental analyses using data from 137 actual patients both demonstrated that our method could discriminate between patient groups in a high-dimensional feature space, thereby rationally expanding the beneficiary population of immunotherapy.