The genetic architecture of non-small cell lung cancer (NSCLC) is relevant to smoking status. However, the genetic contribution of long-term smoking cessation to the prognosis of NSCLC patients remains largely unknown. We conducted a genome-wide association study primarily on the prognosis of 1299 NSCLC patients of long-term former smokers from independent discovery (n = 566) and validation (n = 733) sets, and used in-silico function prediction and multi-omics analysis to identify single nucleotide polymorphisms (SNPs) on prognostics with NSCLC. We further detected SNPs with at least moderate association strength on survival within each group of never, short-term former, long-term former, and current smokers, and compared their genetic similarity at the SNP, gene, expression quantitative trait loci (eQTL), enhancer, and pathway levels. We identified two SNPs, rs34211819TNS3 at 7p12.3 (P = 3.90 × 10-9) and rs1143149SEPT7 at 7p14.2 (P = 9.75 × 10-9), were significantly associated with survival of NSCLC patients who were long-term former smokers. Both SNPs had significant interaction effects with years of smoking cessation (rs34211819TNS3: Pinteraction = 8.0 × 10-4; rs1143149SEPT7: Pinteraction = 0.003). In addition, in silico function prediction and multi-omics analysis provided evidence that these QTLs were associated with survival. Moreover, comparison analysis found higher genetic similarity between long-term former smokers and never-smokers, compared to short-term former smokers or current smokers. Pathway enrichment analysis indicated a unique pattern among long-term former smokers that was related to immune pathways. This study provides important insights into the genetic architecture associated with long-term former smoking NSCLC.

Recent studies suggest that rare variants exhibit stronger effect sizes and might play a crucial role in the etiology of lung cancers (LC). Whole exome plus targeted sequencing of germline DNA was performed on 1045 LC cases and 885 controls in the discovery set. To unveil the inherited causal variants, we focused on rare and predicted deleterious variants and small indels enriched in cases or controls. Promising candidates were further validated in a series of 26,803 LCs and 555,107 controls. During discovery, we identified 25 rare deleterious variants associated with LC susceptibility, including 13 reported in ClinVar. Of the five validated candidates, we discovered two pathogenic variants in known LC susceptibility loci, ATM p.V2716A (Odds Ratio [OR] 19.55, 95%CI 5.04-75.6) and MPZL2 p.I24M frameshift deletion (OR 3.88, 95%CI 1.71-8.8); and three in novel LC susceptibility genes, POMC c.*28delT at 3' UTR (OR 4.33, 95%CI 2.03-9.24), STAU2 p.N364M frameshift deletion (OR 4.48, 95%CI 1.73-11.55), and MLNR p.Q334V frameshift deletion (OR 2.69, 95%CI 1.33-5.43). The potential cancer-promoting role of selected candidate genes and variants was further supported by endogenous DNA damage assays. Our analyses led to the identification of new rare deleterious variants with LC susceptibility. However, in-depth mechanistic studies are still needed to evaluate the pathogenic effects of these specific alleles.

Malignant pleural mesothelioma (MPM) is a rare but lethal pleural cancer with high intratumor heterogeneity (ITH). A recent study in lung adenocarcinoma has developed a clonal gene signature (ORACLE) from multiregional transcriptomic data and demonstrated high prognostic values and reproducibility. However, such a strategy has not been tested in other types of cancer with high ITH. We aimed to identify biomarkers from multi-regional data to prognostically stratify MPM patients. We generated a multiregional RNA-seq dataset for 78 tumor samples obtained from 26 MPM patients, each with one sample collected from a superior, lateral, and inferior region of the tumor. By integrating this dataset with the Cancer Genome Atlas MPM RNA-seq data, we selected 29 prognostic genes displaying high variability across different tumors but low ITH, which named PRACME (Prognostic Risk Associated Clonal Mesothelioma Expression). We evaluated PRACME in two independent MPM datasets and demonstrated its prognostic values. Patients with high signature scores are associated with poor prognosis after adjusting established clinical factors. Interestingly, the PRACME and the ORACLE signatures defined respectively from MPM and lung adenocarcinoma cross-predict prognosis between the two cancer types. Further investigation indicated that the cross-prediction ability might be explained by the high similarity between the two cancer types in their genomic regions with copy number variation, which host many clonal genes. Overall, our clonal signature PRACME provided prognostic stratification in MPM and this study emphasized the importance of multi-regional transcriptomic data for prognostic stratification based on clonal genes.