Independent Validation of Early-Stage Non-Small Cell Lung Cancer Prognostic Scores Incorporating Epigenetic and Transcriptional Biomarkers With Gene-Gene Interactions and Main Effects

Immunogenic cell death-related genes predict prognosis and response to immunotherapy in lung squamous cell carcinoma

Lung squamous cell carcinoma (LUSC) is a malignancy with limited therapeutic options. Immunogenic cell death (ICD) has the potential to enhance the efficacy of cancer therapy by triggering immune responses. We aimed to explore the potential of ICD-based classification in predicting prognosis and response to immunotherapy for LUSC. RNA-seq information and clinical data of LUSC patients were obtained from The Cancer Genome Atlas (TCGA) dataset. ICD-related gene expressions in LUSC samples were analyzed by consensus clustering. Subsequently, differentially expressed genes (DEGs) between different ICD-related subsets were analyzed. Tumor mutation burden, immune cell infiltration, and survival analyses were conducted between different ICD subsets. Finally, an ICD-related risk signature was constructed and evaluated in LUSC patients, and the immunotherapy responses based on the gene expressions were also forecasted. ICD-high and ICD-low groups were defined, and 1466 DEGs were identified between the two subtypes. These DEGs were mainly enriched in collagen-containing extracellular matrix, cytokine-cytokine receptor interaction, the PI3K-Akt signaling pathway, and neuroactive ligand-receptor interaction. Furthermore, the ICD-low group exhibited a favorable prognosis, enhanced TTN and MUC16 mutation frequencies, increased infiltrating immune cells, and downregulated immune checkpoint expressions. Furthermore, we demonstrated that an ICD-related model (based on CD4, NLRP3, NT5E, and TLR4 genes) could forecast the prognosis of LUSC, and ICD risk scores were lower in the responder group. In summary, the predicted values of ICD-related genes (CD4, NLRP3, NT5E, and TLR4) for the prognosis and response to immunotherapy in LUSC were verified in the study, which benefits immunotherapy-based interventions for LUSC patients.

A two-phase epigenome-wide four-way gene-smoking interaction study of overall survival for early-stage non-small cell lung cancer

High-order interactions associated with non-small cell lung cancer (NSCLC) survival may elucidate underlying molecular mechanisms and identify potential therapeutic targets. Our previous work has identified a three-way interaction among pack-year of smoking (the number of packs of cigarettes smoked per day multiplied by the number of years the person has smoked) and two DNA methylation probes (cg05293407TRIM27 and cg00060500KIAA0226). However, whether a four-way interaction exists remains unclear. Therefore, we adopted a two-phase design to identify the four-way gene-smoking interactions by a hill-climbing strategy on the basis of the previously detected three-way interaction. One CpG probe, cg16658473SHISA9, was identified with FDR-q ≤ 0.05 in the discovery phase and P ≤ 0.05 in the validation phase. Meanwhile, the four-way interaction improved the discrimination ability for the prognostic prediction model, as indicated by the area under the receiver operating characteristic curve (AUC) for both 3- and 5-year survival. In summary, we identified a four-way interaction associated with NSCLC survival among pack-year of smoking, cg05293407TRIM27, cg00060500KIAA0226 and g16658473SHISA9, providing novel insights into the complex mechanisms underlying NSCLC progression.

Tumor Methylation Burden (TMeB) in Non-Small Cell Lung Cancer: A New Way of Thinking About Epigenetics

Lung cancer represents a substantial proportion of cancer-associated mortality worldwide, with non-small cell lung cancer (NSCLC) accounting for most of these cases [...].

ARTEMIS: An independently validated prognostic prediction model of breast cancer incorporating epigenetic biomarkers with main effects and gene-gene interactions

INTRODUCTION

Breast cancer, a heterogeneous disease, is influenced by multiple genetic and epigenetic factors. The majority of prognostic models for breast cancer focus merely on the main effects of predictors, disregarding the crucial impacts of gene-gene interactions on prognosis.

OBJECTIVES

Using DNA methylation data derived from nine independent breast cancer cohorts, we developed an independently validated prognostic prediction model of breast cancer incorporating epigenetic biomarkers with main effects and gene-gene interactions (ARTEMIS) with an innovative 3-D modeling strategy. ARTEMIS was evaluated for discrimination ability using area under the receiver operating characteristics curve (AUC), and calibration using expected and observed (E/O) ratio. Additionally, we conducted decision curve analysis to evaluate its clinical efficacy by net benefit (NB) and net reduction (NR). Furthermore, we conducted a systematic review to compare its performance with existing models.

RESULTS

ARTEMIS exhibited excellent risk stratification ability in identifying patients at high risk of mortality. Compared to those below the 25th percentile of ARTEMIS scores, patients with above the 90th percentile had significantly lower overall survival time (HR = 15.43, 95% CI: 9.57-24.88, P = 3.06 × 10-29). ARTEMIS demonstrated satisfactory discrimination ability across four independent populations, with pooled AUC3-year = 0.844 (95% CI: 0.805-0.883), AUC5-year = 0.816 (95% CI: 0.775-0.857), and C-index = 0.803 (95% CI: 0.776-0.830). Meanwhile, ARTEMIS had well calibration performance with pooled E/O ratio 1.060 (95% CI: 1.038-1.083) and 1.090 (95% CI: 1.057-1.122) for 3- and 5-year survival prediction, respectively. Additionally, ARTEMIS is a clinical instrument with acceptable cost-effectiveness for detecting breast cancer patients at high risk of mortality (Pt = 0.4: NB3-year = 19‰, NB5-year = 62‰; NR3-year = 69.21%, NR5-year = 56.01%). ARTEMIS has superior performance compared to existing models in terms of accuracy, extrapolation, and sample size, as indicated by the systematic review. ARTEMIS is implemented as an interactive online tool available at http://bigdata.njmu.edu.cn/ARTEMIS/.

CONCLUSION

ARTEMIS is an efficient and practical tool for breast cancer prognostic prediction.

Cross-attention enables deep learning on limited omics-imaging-clinical data of 130 lung cancer patients

Deep-learning tools that extract prognostic factors derived from multi-omics data have recently contributed to individualized predictions of survival outcomes. However, the limited size of integrated omics-imaging-clinical datasets poses challenges. Here, we propose two biologically interpretable and robust deep-learning architectures for survival prediction of non-small cell lung cancer (NSCLC) patients, learning simultaneously from computed tomography (CT) scan images, gene expression data, and clinical information. The proposed models integrate patient-specific clinical, transcriptomic, and imaging data and incorporate Kyoto Encyclopedia of Genes and Genomes (KEGG) and Reactome pathway information, adding biological knowledge within the learning process to extract prognostic gene biomarkers and molecular pathways. While both models accurately stratify patients in high- and low-risk groups when trained on a dataset of only 130 patients, introducing a cross-attention mechanism in a sparse autoencoder significantly improves the performance, highlighting tumor regions and NSCLC-related genes as potential biomarkers and thus offering a significant methodological advancement when learning from small imaging-omics-clinical samples.

Deep learning model for the prediction of all-cause mortality among long term care people in China: a prospective cohort study

This study aimed to develop a deep learning model to predict the risk stratification of all-cause death for older people with disability, providing guidance for long-term care plans. Based on the government-led long-term care insurance program in a pilot city of China from 2017 and followed up to 2021, the study included 42,353 disabled adults aged over 65, with 25,071 assigned to the training set and 17,282 to the validation set. The administrative data (including baseline characteristics, underlying medical conditions, and all-cause mortality) were collected to develop a deep learning model by least absolute shrinkage and selection operator. After a median follow-up time of 14 months, 17,565 (41.5%) deaths were recorded. Thirty predictors were identified and included in the final models for disability-related deaths. Physical disability (mobility, incontinence, feeding), adverse events (pressure ulcers and falls from bed), and cancer were related to poor prognosis. A total of 10,127, 25,140 and 7086 individuals were classified into low-, medium-, and high-risk groups, with actual risk probabilities of death of 9.5%, 45.8%, and 85.5%, respectively. This deep learning model could facilitate the prevention of risk factors and provide guidance for long-term care model planning based on risk stratification.

Biomarkers and experimental models for cancer immunology investigation

The rapid advancement of tumor immunotherapies poses challenges for the tools used in cancer immunology research, highlighting the need for highly effective biomarkers and reproducible experimental models. Current immunotherapy biomarkers encompass surface protein markers such as PD-L1, genetic features such as microsatellite instability, tumor-infiltrating lymphocytes, and biomarkers in liquid biopsy such as circulating tumor DNAs. Experimental models, ranging from 3D in vitro cultures (spheroids, submerged models, air-liquid interface models, organ-on-a-chips) to advanced 3D bioprinting techniques, have emerged as valuable platforms for cancer immunology investigations and immunotherapy biomarker research. By preserving native immune components or coculturing with exogenous immune cells, these models replicate the tumor microenvironment in vitro. Animal models like syngeneic models, genetically engineered models, and patient-derived xenografts provide opportunities to study in vivo tumor-immune interactions. Humanized animal models further enable the simulation of the human-specific tumor microenvironment. Here, we provide a comprehensive overview of the advantages, limitations, and prospects of different biomarkers and experimental models, specifically focusing on the role of biomarkers in predicting immunotherapy outcomes and the ability of experimental models to replicate the tumor microenvironment. By integrating cutting-edge biomarkers and experimental models, this review serves as a valuable resource for accessing the forefront of cancer immunology investigation.

Construction and validation of cuproptosis-related lncRNA prediction signature for bladder cancer and immune infiltration analysis

Bladder cancer (BC) is a common urologic tumor with a high recurrence rate. Cuproptosis and long noncoding RNAs (lncRNAs) have demonstrated essential roles in the tumorigenesis of many malignancies. Nevertheless, the prognostic value of cuproptosis-related lncRNA (CRLs) in BC is still unclear. The public data used for this study were acquired from the Cancer Genome Atlas database. A comprehensive exploration of the expression profile, mutation, co-expression, and enrichment analyses of cuproptosis-related genes was performed. A total of 466 CRLs were identified using Pearson's correlation analysis. 16 prognostic CRLs were then retained by univariate Cox regression. Unsupervised clustering divided the patients into two clusters with diverse survival outcomes. The signature consists of 7 CRLs was constructed using the least absolute shrinkage and selection operator (LASSO) Cox regression analyses. Survival curves and receiver operating characteristics showed the prognostic signature possessed good predictive value, which was validated in the testing and entire sets. The reliability and stability of our signature were further confirmed by stratified analysis. Additionally, the signature-based risk score was confirmed as an independent prognostic factor. Gene set enrichment analysis showed molecular alteration in the high-risk group was closely associated with cancer. We then developed the clinical nomogram using independent prognostic indicators. Notably, the infiltration of immune cells and expression of immune checkpoints were higher in the high-risk group, suggesting that they may benefit more from immunotherapy. In summary, the prognostic signature might effectively predict the prognosis and provide new insight into the clinical treatment of BC patients.

COVPRIG robustly predicts the overall survival of IDH wild-type glioblastoma and highlights METTL1+ neural-progenitor-like tumor cell in driving unfavorable outcome

BACKGROUND

Accurately predicting the outcome of isocitrate dehydrogenase (IDH) wild-type glioblastoma (GBM) remains hitherto challenging. This study aims to Construct and Validate a Robust Prognostic Model for IDH wild-type GBM (COVPRIG) for the prediction of overall survival using a novel metric, gene-gene (G × G) interaction, and explore molecular and cellular underpinnings.

METHODS

Univariate and multivariate Cox regression of four independent trans-ethnic cohorts containing a total of 800 samples. Prediction efficacy was comprehensively evaluated and compared with previous models by a systematic literature review. The molecular underpinnings of COVPRIG were elucidated by integrated analysis of bulk-tumor and single-cell based datasets.

RESULTS

Using a Cox-ph model-based method, six of the 93,961 G × G interactions were screened to form an optimal combination which, together with age, comprised the COVPRIG model. COVPRIG was designed for RNA-seq and microarray, respectively, and effectively identified patients at high risk of mortality. The predictive performance of COVPRIG was satisfactory, with area under the curve (AUC) ranging from 0.56 (CGGA693, RNA-seq, 6-month survival) to 0.79 (TCGA RNAseq, 18-month survival), which can be further validated by decision curves. Nomograms were constructed for individual risk prediction for RNA-seq and microarray-based cohorts, respectively. Besides, the prognostic significance of COVPRIG was also validated in GBM including the IDH mutant samples. Notably, COVPRIG was comprehensively evaluated and externally validated, and a systemic review disclosed that COVPRIG outperformed current validated models with an integrated discrimination improvement (IDI) of 6-16%. Moreover, integrative bioinformatics analysis predicted an essential role of METTL1+ neural-progenitor-like (NPC-like) malignant cell in driving unfavorable outcome.

CONCLUSION

This study provided a powerful tool for the outcome prediction for IDH wild-type GBM, and preliminary molecular underpinnings for future research.

ATHENA: an independently validated autophagy-related epigenetic prognostic prediction model of head and neck squamous cell carcinoma

The majority of these existing prognostic models of head and neck squamous cell carcinoma (HNSCC) have unsatisfactory prediction accuracy since they solely utilize demographic and clinical information. Leveraged by autophagy-related epigenetic biomarkers, we aim to develop a better prognostic prediction model of HNSCC incorporating CpG probes with either main effects or gene-gene interactions. Based on DNA methylation data from three independent cohorts, we applied a 3-D analysis strategy to develop An independently validated auTophagy-related epigenetic prognostic prediction model of HEad and Neck squamous cell carcinomA (ATHENA). Compared to prediction models with only demographic and clinical information, ATHENA has substantially improved discriminative ability, prediction accuracy and more clinical net benefits, and shows robustness in different subpopulations, as well as external populations. Besides, epigenetic score of ATHENA is significantly associated with tumor immune microenvironment, tumor-infiltrating immune cell abundances, immune checkpoints, somatic mutation and immunity-related drugs. Taken together these results, ATHENA has the demonstrated feasibility and utility of predicting HNSCC survival ( http://bigdata.njmu.edu.cn/ATHENA/ ).

Mutational and Transcriptional Characterization Establishes Prognostic Models for Resectable Lung Squamous Cell Carcinoma

Background

The prognosis of non-small cell lung cancer (NSCLC) patients has been comprehensively studied. However, the prognosis of resectable (stage I-IIIA) lung squamous cell carcinoma (LUSC) has not been thoroughly investigated at genomic and transcriptional levels.

Methods

Data of genomic alterations and transcriptional-level changes of 355 stage I-IIIA LUSC patients were downloaded from The Cancer Genome Atlas (TCGA) database, together with the clinicopathological information (training cohort). A validation cohort of 91 patients was retrospectively recruited. Data were analyzed and figures were plotted using the R software.

Results

Training cohort was established with 355 patients. TP53 (78%), TTN (68%), CSMD3 (39%), MUT16 (36%) and RYR2 (36%) were genes with the highest mutational frequency. BRINP3, COL11A1, GRIN2B, MUC5B, NLRP3 and TENM3 exhibited significant higher mutational frequency in stage III (P < 0.05). Patients with stage III also exhibited significantly higher tumor mutational burden (TMB) than those with stage I (P < 0.01). The mutational status of 10 genes were found to have significant stratification on patient prognosis. TMB at threshold of 25 percentile (TMB = 2.39 muts/Mb) also significantly stratified the patient prognosis (P = 0.0003). Univariate and multivariate analyses revealed TTN, ADGRB3, MYH7 and MYH15 mutational status and TMB as independent risk factors. Further analysis of transcriptional profile revealed many significantly up- and down-regulated genes, and multivariate analysis found the transcriptional levels of seven genes as independent risk factors. Significant factors from the multivariate analyses were used to establish a Nomogram model to quantify the risk in prognosis of individual LUSC patients. The model was validated with a cohort containing 91 patients, which showed good predicting efficacy and consistency.

Conclusion

The influencing factors of prognosis of stage I-III LUSC patients have been revealed. Risk factors including gender, T stage, cancer location, and the mutational and transcriptional status of several genes were used to establish a Nomogram model to assess the patient prognosis. Subsequent validation proved its effectiveness.

A trans-omics assessment of gene-gene interaction in early-stage NSCLC

Epigenome-wide gene-gene (G × G) interactions associated with non-small-cell lung cancer (NSCLC) survival may provide insights into molecular mechanisms and therapeutic targets. Hence, we proposed a three-step analytic strategy to identify significant and robust G × G interactions that are relevant to NSCLC survival. In the first step, among 49 billion pairs of DNA methylation probes, we identified 175 775 G × G interactions with PBonferroni  ≤ 0.05 in the discovery phase of epigenomic analysis; among them, 15 534 were confirmed with P ≤ 0.05 in the validation phase. In the second step, we further performed a functional validation for these G × G interactions at the gene expression level by way of a two-phase (discovery and validation) transcriptomic analysis, and confirmed 25 significant G × G interactions enriched in the 6p21.33 and 6p22.1 regions. In the third step, we identified two G × G interactions using the trans-omics analysis, which had significant (P ≤ 0.05) epigenetic cis-regulation of transcription and robust G × G interactions at both the epigenetic and transcriptional levels. These interactions were cg14391855 × cg23937960 (βinteraction  = 0.018, P = 1.87 × 10-12 ), which mapped to RELA × HLA-G (βinteraction  = 0.218, P = 8.82 × 10-11 ) and cg08872738 × cg27077312 (βinteraction  = -0.010, P = 1.16 × 10-11 ), which mapped to TUBA1B × TOMM40 (βinteraction =-0.250, P = 3.83 × 10-10 ). A trans-omics mediation analysis revealed that 20.3% of epigenetic effects on NSCLC survival were significantly (P = 0.034) mediated through transcriptional expression. These statistically significant trans-omics G × G interactions can also discriminate patients with high risk of mortality. In summary, we identified two G × G interactions at both the epigenetic and transcriptional levels, and our findings may provide potential clues for precision treatment of NSCLC.

Gene Screening in High-Throughput Right-Censored Lung Cancer Data

Background

Advances in sequencing technologies have allowed collection of massive genome-wide information that substantially advances lung cancer diagnosis and prognosis. Identifying influential markers for clinical endpoints of interest has been an indispensable and critical component of the statistical analysis pipeline. However, classical variable selection methods are not feasible or reliable for high-throughput genetic data. Our objective is to propose a model-free gene screening procedure for high-throughput right-censored data, and to develop a predictive gene signature for lung squamous cell carcinoma (LUSC) with the proposed procedure.

Methods

A gene screening procedure was developed based on a recently proposed independence measure. The Cancer Genome Atlas (TCGA) data on LUSC was then studied. The screening procedure was conducted to narrow down the set of influential genes to 378 candidates. A penalized Cox model was then fitted to the reduced set, which further identified a 6-gene signature for LUSC prognosis. The 6-gene signature was validated on datasets from the Gene Expression Omnibus.

Results

Both model-fitting and validation results reveal that our method selected influential genes that lead to biologically sensible findings as well as better predictive performance, compared to existing alternatives. According to our multivariable Cox regression analysis, the 6-gene signature was indeed a significant prognostic factor (p-value < 0.001) while controlling for clinical covariates.

Conclusions

Gene screening as a fast dimension reduction technique plays an important role in analyzing high-throughput data. The main contribution of this paper is to introduce a fundamental yet pragmatic model-free gene screening approach that aids statistical analysis of right-censored cancer data, and provide a lateral comparison with other available methods in the context of LUSC.

A Large-Scale Genome-Wide Gene-Gene Interaction Study of Lung Cancer Susceptibility in Europeans With a Trans-Ethnic Validation in Asians

INTRODUCTION

Although genome-wide association studies have been conducted to investigate genetic variation of lung tumorigenesis, little is known about gene-gene (G × G) interactions that may influence the risk of non-small cell lung cancer (NSCLC).

METHODS

Leveraging a total of 445,221 European-descent participants from the International Lung Cancer Consortium OncoArray project, Transdisciplinary Research in Cancer of the Lung and UK Biobank, we performed a large-scale genome-wide G × G interaction study on European NSCLC risk by a series of analyses. First, we used BiForce to evaluate and rank more than 58 billion G × G interactions from 340,958 single-nucleotide polymorphisms (SNPs). Then, the top interactions were further tested by demographically adjusted logistic regression models. Finally, we used the selected interactions to build lung cancer screening models of NSCLC, separately, for never and ever smokers.

RESULTS

With the Bonferroni correction, we identified eight statistically significant pairs of SNPs, which predominantly appeared in the 6p21.32 and 5p15.33 regions (e.g., rs521828C6orf10 and rs204999PRRT1, ORinteraction = 1.17, p = 6.57 × 10-13; rs3135369BTNL2 and rs2858859HLA-DQA1, ORinteraction = 1.17, p = 2.43 × 10-13; rs2858859HLA-DQA1 and rs9275572HLA-DQA2, ORinteraction = 1.15, p = 2.84 × 10-13; rs2853668TERT and rs62329694CLPTM1L, ORinteraction = 0.73, p = 2.70 × 10-13). Notably, even with much genetic heterogeneity across ethnicities, three pairs of SNPs in the 6p21.32 region identified from the European-ancestry population remained significant among an Asian population from the Nanjing Medical University Global Screening Array project (rs521828C6orf10 and rs204999PRRT1, ORinteraction = 1.13, p = 0.008; rs3135369BTNL2 and rs2858859HLA-DQA1, ORinteraction = 1.11, p = 5.23 × 10-4; rs3135369BTNL2 and rs9271300HLA-DQA1, ORinteraction = 0.89, p = 0.006). The interaction-empowered polygenetic risk score that integrated classical polygenetic risk score and G × G information score was remarkable in lung cancer risk stratification.

CONCLUSIONS

Important G × G interactions were identified and enriched in the 5p15.33 and 6p21.32 regions, which may enhance lung cancer screening models.

Epigenome-wide gene–age interaction study reveals reversed effects of MORN1 DNA methylation on survival between young and elderly oral squamous cell carcinoma patients

DNA methylation serves as a reversible and prognostic biomarker for oral squamous cell carcinoma (OSCC) patients. It is unclear whether the effect of DNA methylation on OSCC overall survival varies with age. As a result, we performed a two-phase gene–age interaction study of OSCC prognosis on an epigenome-wide scale using the Cox proportional hazards model. We identified one CpG probe, cg11676291_MORN1, whose effect was significantly modified by age (HR_discovery = 1.018, p = 4.07 × 10⁻⁰⁷, FDR-q = 3.67 × 10⁻⁰²; HR_validation = 1.058, p = 8.09 × 10⁻⁰³; HR_combined = 1.019, p = 7.36 × 10⁻¹⁰). Moreover, there was an antagonistic interaction between hypomethylation of cg11676291_MORN1 and age (HR_interaction = 0.284; 95% CI, 0.135–0.597; p = 9.04 × 10⁻⁰⁴). The prognosis of OSCC patients was well discriminated by the prognostic score incorporating cg11676291_MORN1–age interaction (HR_{high vs. low} = 3.66, 95% CI: 2.40–5.60, p = 1.93 × 10⁻⁰⁹). By adding 24 significant gene–age interactions using a looser criterion, we significantly improved the area under the receiver operating characteristic curve (AUC) of the model at 3- and 5-year prognostic prediction (AUC_3-year = 0.80, AUC_5-year = 0.79, C-index = 0.75). Our study identified a significant interaction between cg11676291_MORN1 and age on OSCC survival, providing a potential therapeutic target for OSCC patients.

A Novel Cuproptosis-Related Prognostic Gene Signature and Validation of Differential Expression in Clear Cell Renal Cell Carcinoma

Clear cell renal cell carcinoma (ccRCC) is the most prevalent subtype of renal cell carcinoma, which is characterized by metabolic reprogramming. Cuproptosis, a novel form of cell death, is highly linked to mitochondrial metabolism and mediated by protein lipoylation. However, the clinical impacts of cuproptosis-related genes (CRGs) in ccRCC largely remain unclear. In the current study, we systematically evaluated the genetic alterations of cuproptosis-related genes in ccRCC. Our results revealed that CDKN2A, DLAT, DLD, FDX1, GLS, PDHA1 and PDHB exhibited differential expression between ccRCC and normal tissues (|log2(fold change)| > 2/3 and p < 0.05). Utilizing an iterative sure independence screening (SIS) method, we separately constructed the prognostic signature of CRGs for predicting the overall survival (OS) and progression-free survival (PFS) in ccRCC patients. The prognostic score of CRGs yielded an area under the curve (AUC) of 0.658 and 0.682 for the prediction of 5-year OS and PFS, respectively. In the Kaplan−Meier survival analysis of OS, a higher risk score of cuproptosis-related gene signature was significantly correlated with worse overall survival (HR = 2.72 (2.01−3.68), log-rank p = 1.76 × 10−7). Patients with a higher risk had a significantly shorter PFS (HR = 2.83 (2.08−3.85), log-rank p = 3.66 × 10−7). Two independent validation datasets (GSE40435 (N = 101), GSE53757 (N = 72)) were collected for meta-analysis, suggesting that CDKN2A (log2(fold change) = 1.46, 95%CI: 1.75−2.35) showed significantly higher expression in ccRCC tissues while DLAT (log2(fold change) = −0.54, 95%CI: −0.93−−0.15) and FDX1 (log2(fold change) = −1.01, 95%CI: −1.61−−0.42) were lowly expressed. The expression of CDKN2A and FDX1 in ccRCC was also significantly associated with immune infiltration levels and programmed cell death protein 1 (PD-1) expression (CDKN2A: r = 0.24, p = 2.14 × 10−8; FDX1: r = −0.17, p = 1.37 × 10−4). In conclusion, the cuproptosis-related gene signature could serve as a potential prognostic predictor for ccRCC patients and may offer novel insights into the cancer treatment.

APOLLO: An accurate and independently validated prediction model of lower-grade gliomas overall survival and a comparative study of model performance

BACKGROUND

Virtually few accurate and robust prediction models of lower-grade gliomas (LGG) survival exist that may aid physicians in making clinical decisions. We aimed to develop a prognostic prediction model of LGG by incorporating demographic, clinical and transcriptional biomarkers with either main effects or gene-gene interactions.

METHODS

Based on gene expression profiles of 1,420 LGG patients from six independent cohorts comprising both European and Asian populations, we proposed a 3-D analysis strategy to develop and validate an Accurate Prediction mOdel of Lower-grade gLiomas Overall survival (APOLLO). We further conducted decision curve analysis to assess the net benefit (NB) of identifying true positives and the net reduction (NR) of unnecessary interventions. Finally, we compared the performance of APOLLO and the existing prediction models by the first systematic review.

FINDINGS

APOLLO possessed an excellent discriminative ability to identify patients at high mortality risk. Compared to those with less than the 20th percentile of APOLLO risk score, patients with more than the 90th percentile of APOLLO risk score had significantly worse overall survival (HR=54·18, 95% CI: 34·73-84·52, P=2·66 × 10-69). Further, APOLLO can accurately predict both 36- and 60-month survival in six independent cohorts with a pooled AUC36-month=0·901 (95% CI: 0·879-0·923), AUC60-month=0·843 (95% CI: 0·815-0·871) and C-index=0·818 (95% CI: 0·800-0·835). Moreover, APOLLO offered an effective screening strategy for detecting LGG patients susceptible to death (NB36-month=0·166, NR36-month=40·1% and NB60-month=0·258, NR60-month=19·2%). The systematic comparisons revealed APOLLO outperformed the existing models in accuracy and robustness.

INTERPRETATION

APOLLO has the demonstrated feasibility and utility of predicting LGG survival (http://bigdata.njmu.edu.cn/APOLLO).

FUNDING

National Key Research and Development Program of China (2016YFE0204900); Natural Science Foundation of Jiangsu Province (BK20191354); National Natural Science Foundation of China (81973142 and 82103946); China Postdoctoral Science Foundation (2020M681671); National Institutes of Health (CA209414, CA249096, CA092824 and ES000002).

MiR-30a-3p Suppresses the Growth and Development of Lung Adenocarcinoma Cells Through Modulating GOLM1/JAK-STAT Signaling

A considerable amount of people succumbs to lung adenocarcinoma (LUAD) due to its high incidence and mortality. This study attempted to reveal the impacts of GOLM1 on LUAD. This work analyzed GOLM1 expression in LUAD and normal tissue and studied its prognostic value utilizing data from The Cancer Genome Atlas. RNA and protein levels were, respectively, determined utilizing qRT-PCR and western blot. Cell-aggressive behaviors were assessed employing Cell Counting Kit-8, scratch healing, and Transwell assays. The targetting relationship between GOLM1 and miR-30a-3p was assayed by dual-luciferase method. GOLM1 up-regulation in LUAD was found in TCGA and it was also a negative factor for survival in patients. GOLM1 overexpression promoted cell progression in LUAD. Down-regulated miR-30a-3p in LUAD was an upstream regulatory miRNA of GOLM1 in terms of molecular mechanism. Further, rescue assays illustrated that miR-30a-3p overexpression attenuated the GOLM1 facilitating impacts on LUAD progression. Finally, we proved that miR-30a-3p/GOLM1 regulated progression of LUAD cells via JAK-STAT pathway. Collectively, the inhibitory impacts of miR-30a-3p on LUAD growth may be mediated by GOLM1/JAK-STAT, which may contribute to the diagnosis of LUAD therapy and the development of therapeutic tools.

The effect and apoptosis mechanism of 6-methoxyflavone in HeLa cells

INTRODUCTION

Tumor cell apoptosis is a crucial indicator for judging the antiproliferative effects of anti-cancer drugs. The detection of optical and macromolecular biomarkers is the most common method for assessing the level of apoptosis. We aimed to explore the anti-tumor mechanisms of 6-methoxyflavone.

MATERIAL AND METHODS

Three optical methods, including the percentage of apoptotic cells, cell morphology, and subcellular ultrastructure changes, were obtained using flow cytometry, inverted fluorescence microscopy, and transmission electron microscope imaging. The mRNA or protein expression of macromolecular biomarkers related to common apoptotic pathways was determined via polymerase chain reactions or western blot assays. The functional role of the core gene biomarker was investigated through overexpression, knockdown, and phosphorylation inhibitor (GSK2656157).

RESULTS

Transcriptome sequencing and the optical biomarkers assays demonstrated that 6-methoxyflavone could induce apoptosis in HeLa cells. The expression of macromolecular biomarkers indicated that 6-methoxyflavone induced apoptosis through the PERK/EIF2α/ATF4/CHOP pathway. Phosphorylated PERK was identified as the core biomarker of this pathway. Both overexpression and GSK2656157 significantly altered the expression level of phosphorylated PERK in 6-methoxyflavone-treated HeLa cells.

DISCUSSION AND CONCLUSION

Macromolecular biomarkers such as phosphorylated PERK and phosphorylated EIF2α are of great significance for assessing the therapeutic effects of 6-methoxyflavone.

Combined Consideration of Tumor-Associated Immune Cell Density and Immune Checkpoint Expression in the Peritumoral Microenvironment for Prognostic Stratification of Non-Small-Cell Lung Cancer Patients

Given the complexity and highly heterogeneous nature of the microenvironment and its effects on antitumor immunity and cancer immune evasion, the prognostic value of a single immune marker is limited. Here, we show how the integration of immune checkpoint molecule expression and tumor-associated immune cell distribution patterns can influence prognosis prediction in non-small-cell lung cancer (NSCLC) patients. We analyzed tissue microarray (TMA) data derived from multiplex immunohistochemistry results and measured the densities of tumor-infiltrating CD8+ and FOXP3+ immune cells and tumor cells (PanCK+), as well as the densities of programmed cell death 1 (PD-1)+ and programmed cell death ligand 1 (PD-L1)+ cells in the peritumor and intratumor subregions. We found a higher density of infiltrating CD8+ and FOXP3+ immune cells in the peritumoral compartment than in the intratumoral compartment. In addition, unsupervised hierarchical clustering analysis of these markers revealed that the combination of high CD8/FOXP3 expression, low PD-1 and PD-L1 immune checkpoint expression, and lack of epidermal growth factor receptor (EGFR) mutation could be a favorable predictive marker. On the other hand, based on the clustering analysis, low CD8/FOXP3 and immune checkpoint (PD-1 and PD-L1) expression might be a marker for patients who are likely to respond to strategies targeting regulatory T (Treg) cells. Furthermore, an immune risk score model was established based on multivariate Cox regression, and the risk score was determined to be an independent prognostic factor for NSCLC patients. These results indicate that the immune context is heterogeneous because of the complex interactions of different components and that using multiple factors in combination might be promising for predicting the prognosis of and stratifying NSCLC patients.

A two-phase comprehensive NSCLC prognostic study identifies lncRNAs with significant main effect and interaction

Long noncoding RNA (lncRNA) are involved in regulating physiological behaviors for various malignant tumors, including non-small-cell lung cancer (NSCLC). However, few studies comprehensively evaluated both lncRNA-lncRNA interaction effects and main effects of lncRNA on overall survival of NSCLC. Hence, we performed a two-phase designed study of lncRNA expression in tumor tissues using 604 NSCLC patients from The Cancer Genome Atlas as the discovery phase and 839 patients from Gene Expression Omnibus as the validation phase. In the discovery phase, we adopted a two-step strategy, Screening before Testing, for dimension reduction and signal detection. These candidate lncRNAs first screened out by the weighted random forest (Ranger), were then tested through the Cox proportional hazards model adjusted for covariates. Significant lncRNAs with either type of effects aforementioned were carried forward into the validation phase to confirm their significances again. As a result, in the discovery phase, 19 lncRNAs were identified by Ranger, among which five lncRNAs and one pair of lncRNA-lncRNA interaction exhibited significant effects (FDR-q ≤ 0.05) main and interaction effects on NSCLC survival, respectively, through Cox model. After the independent validation, we finally observed that one lncRNA (ENSG00000227403.1) with main effect was robustly associated with NSCLC prognosis (HR_discovery = 0.90, P = 1.20 × 10^-3; HR_validation = 0.94, P = 4.11 × 10^-3) and one pair of lncRNAs (ENSG00000267121.4 and ENSG00000272369.1) had significant interaction effect on NSCLC survival (HR_discovery = 1.12, P = 3.07 × 10^-4; HR_validation = 1.11, P = 0.0397). Our comprehensive NSCLC prognostic study of lncRNA provided population-level evidence for further functional study.

Epigenome-wide three-way interaction study identifies a complex pattern between TRIM27, KIAA0226, and smoking associated with overall survival of early-stage NSCLC

The interaction between DNA methylation of tripartite motif containing 27 (cg05293407TRIM27 ) and smoking has previously been identified to reveal histologically heterogeneous effects of TRIM27 DNA methylation on early-stage non-small-cell lung cancer (NSCLC) survival. However, to understand the complex mechanisms underlying NSCLC progression, we searched three-way interactions. A two-phase study was adopted to identify three-way interactions in the form of pack-year of smoking (number of cigarettes smoked per day × number of years smoked) × cg05293407TRIM27  × epigenome-wide DNA methylation CpG probe. Two CpG probes were identified with FDR-q ≤ 0.05 in the discovery phase and P ≤ 0.05 in the validation phase: cg00060500KIAA0226 and cg17479956EXT2 . Compared to a prediction model with only clinical information, the model added 42 significant three-way interactions using a looser criterion (discovery: FDR-q ≤ 0.10, validation: P ≤ 0.05) had substantially improved the area under the receiver operating characteristic curve (AUC) of the prognostic prediction model for both 3-year and 5-year survival. Our research identified the complex interaction effects among multiple environment and epigenetic factors, and provided therapeutic target for NSCLC patients.

A narrative review of prognosis prediction models for non-small cell lung cancer: what kind of predictors should be selected and how to improve models?

Objective

To discover potential predictors and explore how to build better models by summarizing the existing prognostic prediction models of non-small cell lung cancer (NSCLC).

Background

Research on clinical prediction models of NSCLC has experienced explosive growth in recent years. As more predictors of prognosis are discovered, the choice of predictors to build models is particularly important, and in the background of more applications of next-generation sequencing technology, gene-related predictors are widely used. As it is more convenient to obtain samples and follow-up data, the prognostic model is preferred by researchers.

Methods

PubMed and the Cochrane Library were searched using the items “NSCLC”, “prognostic model”, “prognosis prediction”, and “survival prediction” from 1 January 1980 to 5 May 2021. Reference lists from articles were reviewed and relevant articles were identified.

Conclusions

The performance of gene-related models has not obviously improved. Relative to the innovation and diversity of predictors, it is more important to establish a highly stable model that is convenient for clinical application. Most of the prevalent models are highly biased and referring to PROBAST at the beginning of the study may be able to significantly control the bias. Existing models should be validated in a large external dataset to make a meaningful comparison.

Performance of a Machine Learning Algorithm Using Electronic Health Record Data to Identify and Estimate Survival in a Longitudinal Cohort of Patients With Lung Cancer

IMPORTANCE

Electronic health records (EHRs) provide a low-cost means of accessing detailed longitudinal clinical data for large populations. A lung cancer cohort assembled from EHR data would be a powerful platform for clinical outcome studies.

OBJECTIVE

To investigate whether a clinical cohort assembled from EHRs could be used in a lung cancer prognosis study.

DESIGN, SETTING, AND PARTICIPANTS

In this cohort study, patients with lung cancer were identified among 76 643 patients with at least 1 lung cancer diagnostic code deposited in an EHR in Mass General Brigham health care system from July 1988 to October 2018. Patients were identified via a semisupervised machine learning algorithm, for which clinical information was extracted from structured and unstructured data via natural language processing tools. Data completeness and accuracy were assessed by comparing with the Boston Lung Cancer Study and against criterion standard EHR review results. A prognostic model for non-small cell lung cancer (NSCLC) overall survival was further developed for clinical application. Data were analyzed from March 2019 through July 2020.

EXPOSURES

Clinical data deposited in EHRs for cohort construction and variables of interest for the prognostic model were collected.

MAIN OUTCOMES AND MEASURES

The primary outcomes were the performance of the lung cancer classification model and the quality of the extracted variables; the secondary outcome was the performance of the prognostic model.

RESULTS

Among 76 643 patients with at least 1 lung cancer diagnostic code, 42 069 patients were identified as having lung cancer, with a positive predictive value of 94.4%. The study cohort consisted of 35 375 patients (16 613 men [47.0%] and 18 756 women [53.0%]; 30 140 White individuals [85.2%], 1040 Black individuals [2.9%], and 857 Asian individuals [2.4%]) after excluding patients with lung cancer history and less than 14 days of follow-up after initial diagnosis. The median (interquartile range) age at diagnosis was 66.7 (58.4-74.1) years. The area under the receiver operating characteristic curves of the prognostic model for overall survival with NSCLC were 0.828 (95% CI, 0.815-0.842) for 1-year prediction, 0.825 (95% CI, 0.812-0.836) for 2-year prediction, 0.814 (95% CI, 0.800-0.826) for 3-year prediction, 0.814 (95% CI, 0.799-0.828) for 4-year prediction, and 0.812 (95% CI, 0.798-0.825) for 5-year prediction.

CONCLUSIONS AND RELEVANCE

These findings suggest the feasibility of assembling a large-scale EHR-based lung cancer cohort with detailed longitudinal clinical measurements and that EHR data may be applied in cancer progression with a set of generalizable approaches.

Prognostic impact of lung adenocarcinoma second predominant pattern from a large European database

BACKGROUND AND OBJECTIVES

Adenocarcinoma patterns could be grouped based on clinical behaviors: low- (lepidic), intermediate- (papillary or acinar), and high-grade (micropapillary and solid). We analyzed the impact of the second predominant pattern (SPP) on disease-free survival (DFS).

METHODS

We retrospectively collected data of surgically resected stage I and II adenocarcinoma.

SELECTION CRITERIA

anatomical resection with lymphadenectomy and pathological N0. Pure adenocarcinomas and mucinous subtypes were excluded. Recurrence rate and factors affecting DFS were analyzed according to the SPP focusing on intermediate-grade predominant pattern adenocarcinomas.

RESULTS

Among 270 patients, 55% were male. The mean age was 68.3 years. SPP pattern appeared as follows: lepidic 43.0%, papillary 23.0%, solid 14.4%, acinar 11.9%, and micropapillary 7.8%. The recurrence rate was 21.5% and 5-year DFS was 71.1%. No difference in DFS was found according to SPP (p = .522). In patients with high-grade SPP, the percentage of SPP, age, and tumor size significantly influenced DFS (p = .016). In patients with lepidic SPP, size, male gender, and lymph-node sampling (p = .005; p = .014; p = .038, respectively) significantly influenced DFS.

CONCLUSIONS

The impact of SPP on DFS is not homogeneous in a subset of patients with the intermediate-grade predominant patterns. The influence of high-grade SPP on DFS is related to its proportion in the tumor.

Epigenetic-smoking interaction reveals histologically heterogeneous effects of TRIM27 DNA methylation on overall survival among early-stage NSCLC patients

Tripartite motif containing 27 (TRIM27) is highly expressed in lung cancer, including non-small-cell lung cancer (NSCLC). Here, we profiled DNA methylation of lung adenocarcinoma (LUAD) and lung squamous cell carcinoma (LUSC) tumours from 613 early-stage NSCLC patients and evaluated associations between CpG methylation of TRIM27 and overall survival. Significant CpG probes were confirmed in 617 samples from The Cancer Genome Atlas. The methylation of the CpG probe cg05293407_TRIM27 was significantly associated with overall survival in patients with LUSC (HR = 1.65, 95% CI: 1.30-2.09, P = 4.52 × 10^-5), but not in patients with LUAD (HR = 1.08, 95% CI: 0.87-1.33, P = 0.493). As incidence of LUSC is associated with higher smoking intensity compared to LUAD, we investigated whether smoking intensity impacted on the prognostic effect of cg05293407_TRIM27 methylation in NSCLC. LUSC patients had a higher average pack-year of smoking (37.49_LUAD vs 54.79_LUSC, P = 1.03 × 10^-19) and included a higher proportion of current smokers than LUAD patients (28.24%_LUAD vs 34.09%_LUSC, P = 0.037). cg05293407_TRIM27 was significantly associated with overall survival only in NSCLC patients with medium-high pack-year of smoking (HR = 1.58, 95% CI: 1.26-1.96, P = 5.25 × 10^-5). We conclude that cg05293407_TRIM27 methylation is a potential predictor of LUSC prognosis, and smoking intensity may impact on its prognostic value across the various types of NSCLC.