Holcátová, Ivana

Abstract Understanding tumor cell dynamics can improve prognosis and treatment but remains limited for lung adenocarcinoma in people who have never smoked (NS-LUAD). With RNA sequencing data from 684 NS-LUAD cases and validation in an independent dataset, we identified three subtypes with distinct phenotypic traits and cell compositions. Additional genomic and histologic data further characterized the subtypes. “Steady,” marked by low proliferation, high alveolar cell fraction, moderate-to-well differentiation, and fewer driver gene alterations, is linked to prolonged survival and low immune evasion. “Proliferative” shows high proliferation markers, TP53 mutations, and gene fusions. “Chaotic,” with high epithelial-to-mesenchymal transition markers, has the worst prognosis, even within stage I tumors. Lacking known molecular or histologic characteristics, this aggressive subtype is solely identified by transcriptomic data. A 60-gene signature recapitulates the classification and predicts survival even within subgroups based on tumor stage or known genomic features, emphasizing its potential for improving early-stage NS-LUAD prognostication in clinical settings. Significance: The transcriptome of 684 NS-LUAD identifies three subtypes with different cellular dynamics and genomic and morphologic features. A 60-gene signature accurately stratifies subjects for mortality risk, even in stage I, offering a potential clinically applicable tool for treatment decision-making in patients with NS-LUAD. See related commentary by Azizi et al., p. 423

ABSTRACT Understanding lung cancer evolution can identify tools for intercepting its growth. In a landscape analysis of 1024 lung adenocarcinomas (LUAD) with deep whole-genome sequencing integrated with multiomic data, we identified 542 LUAD that displayed diverse clonal architecture. In this group, we observed an interplay between mobile elements, endogenous and exogenous mutational processes, distinct driver genes, and epidemiological features. Our results revealed divergent evolutionary trajectories based on tobacco smoking exposure, ancestry, and sex. LUAD from smokers showed an abundance of tobacco-related C:G>A:T driver mutations in KRAS plus short subclonal diversification. LUAD in never smokers showed early occurrence of copy number alterations and EGFR mutations associated with SBS5 and SBS40a mutational signatures. Tumors harboring EGFR mutations exhibited long latency, particularly in females of European-ancestry (EU_N). In EU_N, EGFR mutations preceded the occurrence of other driver genes, including TP53 and RBM10. Tumors from Asian never smokers showed a short clonal evolution and presented with heterogeneous repetitive patterns for the inferred mutational order. Importantly, we found that the mutational signature ID2 is a marker of a previously unrecognized mechanism for LUAD evolution. Tumors with ID2 showed short latency and high L1 retrotransposon activity linked to L1 promoter demethylation. These tumors exhibited an aggressive phenotype, characterized by increased genomic instability, elevated hypoxia scores, low burden of neoantigens, propensity to develop metastasis, and poor overall survival. Re-activated L1 retrotransposition-induced mutagenesis can contribute to the origin of the mutational signature ID2, including through the regulation of the transcriptional factor ZNF695 , a member of the KZFP family. The complex nature of LUAD evolution creates both challenges and opportunities for screening and treatment plans.

Abstract Background: Somatic EGFR mutations define a subset of non–small cell lung cancers (NSCLC) that have clinical impact on NSCLC risk and outcome. However, EGFR-mutation-status is often missing in epidemiologic datasets. We developed and tested pragmatic approaches to account for EGFR-mutation-status based on variables commonly included in epidemiologic datasets and evaluated the clinical utility of these approaches. Methods: Through analysis of the International Lung Cancer Consortium (ILCCO) epidemiologic datasets, we developed a regression model for EGFR-status; we then applied a clinical-restriction approach using the optimal cut-point, and a second epidemiologic, multiple imputation approach to ILCCO survival analyses that did and did not account for EGFR-status. Results: Of 35,356 ILCCO patients with NSCLC, EGFR-mutation-status was available in 4,231 patients. A model regressing known EGFR-mutation-status on clinical and demographic variables achieved a concordance index of 0.75 (95% CI, 0.74–0.77) in the training and 0.77 (95% CI, 0.74–0.79) in the testing dataset. At an optimal cut-point of probability-score = 0.335, sensitivity = 69% and specificity = 72.5% for determining EGFR-wildtype status. In both restriction-based and imputation-based regression analyses of the individual roles of BMI on overall survival of patients with NSCLC, similar results were observed between overall and EGFR-mutation-negative cohort analyses of patients of all ancestries. However, our approach identified some differences: EGFR-mutated Asian patients did not incur a survival benefit from being obese, as observed in EGFR-wildtype Asian patients. Conclusions: We introduce a pragmatic method to evaluate the potential impact of EGFR-status on epidemiological analyses of NSCLC. Impact: The proposed method is generalizable in the common occurrence in which EGFR-status data are missing.

Background The association among gallbladder disease, cholecystectomy, and pancreatic cancer is unclear. Moreover, time interval between gallbladder disease or cholecystectomy and pancreatic cancer diagnosis is not considered in most previous studies. Aim To quantify the association among gallbladder disease, cholecystectomy, and pancreatic cancer, considering time since first diagnosis of gallbladder disease or cholecystectomy. Methods We used data from nine case-control studies within the Pancreatic Cancer Case-Control Consortium, including 5760 cases of adenocarcinoma of the exocrine pancreas and 8437 controls. We estimated pooled odds ratios and the corresponding 95% confidence intervals by estimating study-specific odds ratios through multivariable unconditional logistic regression models, and then pooling the obtained estimates using fixed-effects models. Results Compared with patients with no history of gallbladder disease, the pooled odds ratio of pancreatic cancer was 1.69 (95% confidence interval, 1.51–1.88) for patients reporting a history of gallbladder disease. The odds ratio was 4.90 (95% confidence interval, 3.45–6.97) for gallbladder disease diagnosed <2 years before pancreatic cancer diagnosis and 1.11 (95% confidence interval, 0.96–1.29) when ≥2 years elapsed. The pooled odds ratio was 1.64 (95% confidence interval, 1.43–1.89) for patients who underwent cholecystectomy, as compared to those without cholecystectomy. The odds ratio was 7.00 (95% confidence interval, 4.13–11.86) for a surgery <2 years before pancreatic cancer diagnosis and 1.28 (95% confidence interval, 1.08–1.53) for a surgery ≥2 years before. Conclusions There appears to be no long-term effect of gallbladder disease on pancreatic cancer risk, and at most a modest one for cholecystectomy. The strong short-term association can be explained by diagnostic bias and reverse causation.

A small number of common susceptibility loci have been identified for pancreatic cancer, one of which is marked by rs401681 in the TERT–CLPTM1L gene region on chromosome 5p15.33. Because this region is characterized by low linkage disequilibrium, we sought to identify whether additional single nucleotide polymorphisms (SNPs) could be related to pancreatic cancer risk, independently of rs401681. We performed an in‐depth analysis of genetic variability of the telomerase reverse transcriptase (TERT) and the telomerase RNA component (TERC) genes, in 5,550 subjects with pancreatic cancer and 7,585 controls from the PANcreatic Disease ReseArch (PANDoRA) and the PanScan consortia. We identified a significant association between a variant in TERT and pancreatic cancer risk (rs2853677, odds ratio = 0.85; 95% confidence interval = 0.80–0.90, p = 8.3 × 10−8). Additional analysis adjusting rs2853677 for rs401681 indicated that the two SNPs are independently associated with pancreatic cancer risk, as suggested by the low linkage disequilibrium between them (r2 = 0.07, D′ = 0.28). Three additional SNPs in TERT reached statistical significance after correction for multiple testing: rs2736100 (p = 3.0 × 10−5), rs4583925 (p = 4.0 × 10−5) and rs2735948 (p = 5.0 × 10−5). In conclusion, we confirmed that the TERT locus is associated with pancreatic cancer risk, possibly through several independent variants.