DNA and histone modifications as potent diagnostic and therapeutic targets to advance non-small cell lung cancer management from the perspective of 3P medicine

Deciphering key roles of B cells in prognostication and tailored therapeutic strategies for lung adenocarcinoma: a multi-omics and machine learning approach towards predictive, preventive, and personalized treatment strategies

Background

Lung adenocarcinoma (LUAD) remains a significant global health challenge, with an urgent need for innovative predictive, preventive, and personalized medicine (PPPM) strategies to improve patient outcomes. This study leveraged multi-omics and machine learning approaches to uncover the prognostic roles of B cells in LUAD, thereby reinforcing the PPPM approach.

Methods

We integrated multi-omics data, including bulk RNA, ATAC-seq, single-cell RNA, and spatial transcriptomics sequencing, to characterize the B cell landscape in LUAD within the PPPM framework. Subsequently, we developed an integrative machine learning program that generated the Scissor+ related B cell score (SRBS). This score was validated in the training and validation sets, and its prognostic value was assessed along with clinical features to develop predictive nomograms. This study further assessed the role of SRBS and SRBS genes in response to immunotherapy and identified personalized drug targets for distinct risk subgroups, with gene expression verified experimentally to ensure tailored medical interventions.

Results

Our analysis identified 79 Scissor+ B cell genes linked to LUAD prognosis, supporting the predictive aspect of PPPM. The SRBS model, which utilizes multiple machine learning algorithms, performed excellently in predicting prognosis and clinical transformation, embodying the preventive and personalized aspects of PPPM. Multifactorial analysis confirmed that SRBS was an independent prognostic factor. We observed varying biological functions and immune cell infiltration in the tumor immune microenvironment (TIME) between the high- and low-SRBS groups, underscoring personalized treatment approaches. Notably, patients with elevated SRBS may exhibit resistance to immunotherapy but show increased sensitivity to chemotherapy and targeted therapies. Additionally, we found that LDHA, as an SRBS gene with significant clinical implications, may regulate the sensitivity of LUAD cells to cisplatin.

Conclusion

This study presents a B cell-associated gene signature that serves as a prognostic marker to facilitate personalized treatment for patients with LUAD, adhering to the principles of PPPM.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13167-024-00390-4.

Clinically relevant immune subtypes based on alternative splicing landscape of immune-related genes for lung cancer advanced PPPM approach

Background

Alternative splicing (AS) occurs in the process of gene post-transcriptional process, which is very important for the correct synthesis and function of protein. The change of AS pattern may lead to the change of expression level or function of lung cancer-related genes, and then affect the occurrence and development of lung cancers. The specific AS pattern might be used as a biomarker for early warning and prognostic assessment of a cancer in the framework of predictive, preventive, and personalized medicine (PPPM; 3PM). AS events of immune-related genes (IRGs) were closely associated with tumor progression and immunotherapy. We hypothesize that IRG-AS events are significantly different in lung adenocarcinomas (LUADs) vs. controls or in lung squamous cell carcinomas (LUSCs) vs. controls. IRG-AS alteration profiling was identified to construct IRG-differentially expressed AS (IRG-DEAS) signature models. Study on the selective AS events of specific IRGs in lung cancer patients might be of great significance for further exploring the pathogenesis of lung cancer, realizing early detection and effective monitoring of lung cancer, finding new therapeutic targets, overcoming drug resistance, and developing more effective therapeutic strategies, and better used for the prediction, diagnosis, prevention, and personalized medicine of lung cancer.

Methods

The transcriptomic, clinical, and AS data of LUADs and LUSCs were downloaded from TCGA and its SpliceSeq databases. IRG-DEAS events were identified in LUAD and LUSC, followed by their functional characteristics, and overall survival (OS) analyses. OS-related IRG-DEAS prognostic models were constructed for LUAD and LUSC with Lasso regression, which were used to classify LUADs and LUSCs into low- and high-risk score groups. Furthermore, the immune cell distribution, immune-related scores, drug sensitivity, mutation status, and GSEA/GSVA status were analyzed between low- and high-risk score groups. Also, low- and high-immunity clusters and AS factor (SF)-OS-related-AS co-expression network and verification of cell function of CELF6 were analyzed in LUAD and LUSC.

Results

Comprehensive analysis of transcriptomic, clinical, and AS data of LUADs and LUSCs identified IRG-AS events in LUAD (n = 1607) and LUSC (n = 1656), including OS-related IRG-AS events in LUAD (n = 127) and LUSC (n = 105). A total of 66 IRG-DEAS events in LUAD and 89 IRG-DEAS events in LUSC were identified compared to controls. The overlapping analysis between IRG-DEASs and OS-related IRG-AS events revealed 14 OS-related IRG-DEAS events for LUAD and 16 OS-related IRG-DEAS events for LUSC, which were used to identify and optimize a 12-OS-related-IRG-DEAS signature prognostic model for LUAD and an 11-OS-related-IRG-DEAS signature prognostic model for LUSC. These two prognostic models effectively divided LUAD or LUSC samples into low- and high-risk score groups that were closely associated with OS, clinical characteristics, and tumor immune microenvironment, with significant gene sets and pathways enriched in the two groups. Moreover, weighted gene co-expression network (WGCNA) and nonnegative matrix factorization method (NMF) analyses identified four OS-relevant subtypes of LUAD and six OS-relevant subtypes of LUSC, and ssGSEA identified five immunity-relevant subtypes of LUAD and five immunity-relevant subtypes of LUSC. Interestingly, splicing factors-OS-related-AS network revealed hub molecule CELF6 was significantly related to the malignant phenotype in lung cancer cells.

Conclusions

This study established two reliable IRG-DEAS signature prognostic models and constructed interesting splicing factor-splicing event networks in LUAD and LUSC, which can be used to construct clinically relevant immune subtypes, patient stratification, prognostic prediction, and personalized medical services in the PPPM practice.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13167-024-00366-4.

Discovering single cannabidiol or synergistic antitumor effects of cannabidiol and cytokine-induced killer cells on non-small cell lung cancer cells

Introduction

A multitude of findings from cell cultures and animal studies are available to support the anti-cancer properties of cannabidiol (CBD). Since CBD acts on multiple molecular targets, its clinical adaptation, especially in combination with cancer immunotherapy regimen remains a serious concern.

Methods

Considering this, we extensively studied the effect of CBD on the cytokine-induced killer (CIK) cell immunotherapy approach using multiple non-small cell lung cancer (NSCLC) cells harboring diverse genotypes.

Results

Our analysis showed that, a) The Transient Receptor Potential Cation Channel Subfamily V Member 2 (TRPV2) channel was intracellularly expressed both in NSCLC cells and CIK cells. b) A synergistic effect of CIK combined with CBD, resulted in a significant increase in tumor lysis and Interferon gamma (IFN-g) production. c) CBD had a preference to elevate the CD25+CD69+ population and the CD62L_CD45RA+terminal effector memory (EMRA) population in NKT-CIK cells, suggesting early-stage activation and effector memory differentiation in CD3+CD56+ CIK cells. Of interest, we observed that CBD enhanced the calcium influx, which was mediated by the TRPV2 channel and elevated phosphor-Extracellular signal-Regulated Kinase (p-ERK) expression directly in CIK cells, whereas ERK selective inhibitor FR180204 inhibited the increasing cytotoxic CIK ability induced by CBD. Further examinations revealed that CBD induced DNA double-strand breaks via upregulation of histone H2AX phosphorylation in NSCLC cells and the migration and invasion ability of NSCLC cells suppressed by CBD were rescued using the TRPV2 antagonist (Tranilast) in the absence of CIK cells. We further investigated the epigenetic effects of this synergy and found that adding CBD to CIK cells decreased the Long Interspersed Nuclear Element-1 (LINE-1) mRNA expression and the global DNA methylation level in NSCLC cells carrying KRAS mutation. We further investigated the epigenetic effects of this synergy and found that adding CBD to CIK cells decreased the Long Interspersed Nuclear Element-1 (LINE-1) mRNA expression and the global DNA methylation level in NSCLC cells carrying KRAS mutation.

Conclusions

Taken together, CBD holds a great potential for treating NSCLC with CIK cell immunotherapy. In addition, we utilized NSCLC with different driver mutations to investigate the efficacy of CBD. Our findings might provide evidence for CBD-personized treatment with NSCLC patients.

Clinically relevant stratification of lung squamous carcinoma patients based on ubiquitinated proteasome genes for 3P medical approach

Relevance

The proteasome is a crucial mechanism that regulates protein fate and eliminates misfolded proteins, playing a significant role in cellular processes. In the context of lung cancer, the proteasome's regulatory function is closely associated with the disease's pathophysiology, revealing multiple connections within the cell. Therefore, studying proteasome inhibitors as a means to identify potential pathways in carcinogenesis and metastatic progression is crucial in in-depth insight into its molecular mechanism and discovery of new therapeutic target to improve its therapy, and establishing effective biomarkers for patient stratification, predictive diagnosis, prognostic assessment, and personalized treatment for lung squamous carcinoma in the framework of predictive, preventive, and personalized medicine (PPPM; 3P medicine).

Methods

This study identified differentially expressed proteasome genes (DEPGs) in lung squamous carcinoma (LUSC) and developed a gene signature validated through Kaplan-Meier analysis and ROC curves. The study used WGCNA analysis to identify proteasome co-expression gene modules and their interactions with the immune system. NMF analysis delineated distinct LUSC subtypes based on proteasome gene expression patterns, while ssGSEA analysis quantified immune gene-set abundance and classified immune subtypes within LUSC samples. Furthermore, the study examined correlations between clinicopathological attributes, immune checkpoints, immune scores, immune cell composition, and mutation status across different risk score groups, NMF clusters, and immunity clusters.

Results

This study utilized DEPGs to develop an eleven-proteasome gene-signature prognostic model for LUSC, which divided samples into high-risk and low-risk groups with significant overall survival differences. NMF analysis identified six distinct LUSC clusters associated with overall survival. Additionally, ssGSEA analysis classified LUSC samples into four immune subtypes based on the abundance of immune cell infiltration with clinical relevance. A total of 145 DEGs were identified between high-risk and low-risk score groups, which had significant biological effects. Moreover, PSMD11 was found to promote LUSC progression by depending on the ubiquitin-proteasome system for degradation.

Conclusions

Ubiquitinated proteasome genes were effective in developing a prognostic model for LUSC patients. The study emphasized the critical role of proteasomes in LUSC processes, such as drug sensitivity, immune microenvironment, and mutation status. These data will contribute to the clinically relevant stratification of LUSC patients for personalized 3P medical approach. Further, we also recommend the application of the ubiquitinated proteasome system in multi-level diagnostics including multi-omics, liquid biopsy, prediction and targeted prevention of chronic inflammation and metastatic disease, and mitochondrial health-related biomarkers, for LUSC 3PM practice.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13167-024-00352-w.

Identifying oxidative stress-related biomarkers in idiopathic pulmonary fibrosis in the context of predictive, preventive, and personalized medicine using integrative omics approaches and machine-learning strategies

Background

Idiopathic pulmonary fibrosis (IPF) is a rare interstitial lung disease with a poor prognosis that currently lacks effective treatment methods. Preventing the acute exacerbation of IPF, identifying the molecular subtypes of patients, providing personalized treatment, and developing individualized drugs are guidelines for predictive, preventive, and personalized medicine (PPPM / 3PM) to promote the development of IPF. Oxidative stress (OS) is an important pathological process of IPF. However, the relationship between the expression levels of oxidative stress-related genes (OSRGs) and clinical indices in patients with IPF is unclear; therefore, it is still a challenge to identify potential beneficiaries of antioxidant therapy. Because PPPM aims to recognize and manage diseases by integrating multiple methods, patient stratification and analysis based on OSRGs and identifying biomarkers can help achieve the above goals.

Methods

Transcriptome data from 250 IPF patients were divided into training and validation sets. Core OSRGs were identified in the training set and subsequently clustered to identify oxidative stress-related subtypes. The oxidative stress scores, clinical characteristics, and expression levels of senescence-associated secretory phenotypes (SASPs) of different subtypes were compared to identify patients who were sensitive to antioxidant therapy to conduct differential gene functional enrichment analysis and predict potential therapeutic drugs. Diagnostic markers between subtypes were obtained by integrating multiple machine learning methods, their expression levels were tested in rat models with different degrees of pulmonary fibrosis and validation sets, and nomogram models were constructed. CIBERSORT, single-cell RNA sequencing, and immunofluorescence staining were used to explore the effects of OSRGs on the immune microenvironment.

Results

Core OSRGs classified IPF into two subtypes. Patients classified into subtypes with low oxidative stress levels had better clinical scores, less severe fibrosis, and lower expression of SASP-related molecules. A reliable nomogram model based on five diagnostic markers was constructed, and these markers' expression stability was verified in animal experiments. The number of neutrophils in the immune microenvironment was significantly different between the two subtypes and was closely related to the degree of fibrosis.

Conclusion

Within the framework of PPPM, this work comprehensively explored the role of OSRGs and their mediated cellular senescence and immune processes in the progress of IPF and assessed their capabilities aspredictors of high oxidative stress and disease progression,targets of the vicious loop between regulated pulmonary fibrosis and OS for targeted secondary and tertiary prevention, andreferences for personalized antioxidant and antifibrotic therapies.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13167-023-00334-4.

Underexplored reciprocity between genome-wide methylation status and long non-coding RNA expression reflected in breast cancer research: potential impacts for the disease management in the framework of 3P medicine

Breast cancer (BC) is the most common female malignancy reaching a pandemic scale worldwide. A comprehensive interplay between genetic alterations and shifted epigenetic regions synergistically leads to disease development and progression into metastatic BC. DNA and histones methylations, as the most studied epigenetic modifications, represent frequent and early events in the process of carcinogenesis. To this end, long non-coding RNAs (lncRNAs) are recognized as potent epigenetic modulators in pathomechanisms of BC by contributing to the regulation of DNA, RNA, and histones' methylation. In turn, the methylation status of DNA, RNA, and histones can affect the level of lncRNAs expression demonstrating the reciprocity of mechanisms involved. Furthermore, lncRNAs might undergo methylation in response to actual medical conditions such as tumor development and treated malignancies. The reciprocity between genome-wide methylation status and long non-coding RNA expression levels in BC remains largely unexplored. Since the bio/medical research in the area is, per evidence, strongly fragmented, the relevance of this reciprocity for BC development and progression has not yet been systematically analyzed. Contextually, the article aims at:consolidating the accumulated knowledge on both-the genome-wide methylation status and corresponding lncRNA expression patterns in BC andhighlighting the potential benefits of this consolidated multi-professional approach for advanced BC management. Based on a big data analysis and machine learning for individualized data interpretation, the proposed approach demonstrates a great potential to promote predictive diagnostics and targeted prevention in the cost-effective primary healthcare (sub-optimal health conditions and protection against the health-to-disease transition) as well as advanced treatment algorithms tailored to the individualized patient profiles in secondary BC care (effective protection against metastatic disease). Clinically relevant examples are provided, including mitochondrial health control and epigenetic regulatory mechanisms involved.