Search

Background: Increased lung-cancer risks for low socioeconomic status (SES) groups are only partially attributable to smoking habits. Little effort has been made to investigate the persistent risks related to low SES by quantification of potential biases., Methods: Based on 12 case-control studies, including 18 centers of the international SYNERGY project (16,550 cases, 20,147 controls), we estimated controlled direct effects (CDE) of SES on lung cancer via multiple logistic regression, adjusted for age, study center, and smoking habits, and stratified by sex. We conducted mediation analysis by inverse odds ratio weighting to estimate natural direct effects (NDE) and natural indirect effects via smoking habits. We considered misclassification of smoking status, selection bias, and unmeasured mediator-outcome confounding by genetic risk, both separately as well as by multiple quantitative bias analysis, using bootstrap to create 95% simulation intervals (SI)., Results: Mediation analysis of lung-cancer risks for SES estimated mean proportions of 43% in men and 33% in women attributable to smoking. Bias analyses decreased direct effects of SES on lung cancer, with selection bias showing the strongest reduction in lung-cancer risk in the multiple bias analysis. Lung-cancer risks remained increased for lower SES groups, with higher risks in men [4th versus 1st (highest) SES quartile: CDE 1.50 (SI 1.32-1.69)] than women [CDE 1.20 (SI 1.01-1.45)]. NDE were similar to CDE, particularly in men., Conclusions: Bias adjustment lowered direct lung-cancer risk estimates of lower SES groups. However, risks for low SES remained elevated, likely attributable to occupational hazards or other environmental exposures., Competing Interests: Competing interests: The authors do not declare any conflict of interest., (Copyright © 2024 Wolters Kluwer Health, Inc. All rights reserved.)

Genome-wide association studies (GWAS) identified over fifty loci associated with lung cancer risk. However, underlying mechanisms and target genes are largely unknown, as most risk-associated variants might regulate gene expression in a context-specific manner. Here, we generate a barcode-shared transcriptome and chromatin accessibility map of 117,911 human lung cells from age/sex-matched ever- and never-smokers to profile context-specific gene regulation. Identified candidate cis-regulatory elements (cCREs) are largely cell type-specific, with 37% detected in one cell type. Colocalization of lung cancer candidate causal variants (CCVs) with these cCREs combined with transcription factor footprinting prioritize the variants for 68% of the GWAS loci. CCV-colocalization and trait relevance score indicate that epithelial and immune cell categories, including rare cell types, contribute to lung cancer susceptibility the most. A multi-level cCRE-gene linking system identifies candidate susceptibility genes from 57% of the loci, where most loci display cell-category-specific target genes, suggesting context-specific susceptibility gene function., (© 2024. This is a U.S. Government work and not under copyright protection in the US; foreign copyright protection may apply.)

Tumor mutational signatures are informative for cancer diagnosis and treatment. However, targeted sequencing, commonly used in clinical settings, lacks specialized analytical tools and a dedicated catalogue of mutational signatures. Here, we introduce SATS, a scalable mutational signature analyzer for targeted sequencing data. SATS leverages tumor mutational burdens to identify and quantify signatures in individual tumors, overcoming the challenges of sparse mutations and variable gene panels. Validations across simulated data, pseudo-targeted sequencing data, and matched whole-genome and targeted sequencing samples show that SATS can accurately detect common mutational signatures and estimate their burdens. Applying SATS to 111,711 tumors from the AACR Project GENIE, we created a pan-cancer mutational signature catalogue specific to targeted sequencing. We further validated signatures in lung, breast and colorectal cancers using an additional 16,774 independent samples. This signature catalogue is a valuable resource for estimating signature burdens in individual targeted sequenced tumors, facilitating the integration of mutational signatures with clinical data.

Background: Lung adenocarcinoma (LUAD) among never-smokers is a public health burden especially prevalent in East Asian (EAS) women. Polygenic risk scores (PRSs), which quanefy geneec suscepebility, are promising for straefying risk, yet have mainly been developed in European (EUR) populaeons. We developed and validated single-and mule-ancestry PRSs for LUAD in EAS never-smokers, using the largest available genome-wide associaeon study (GWAS) dataset., Methods: We used GWAS summary staesecs from both EAS (8,002 cases; 20,782 controls) and EUR (2,058 cases; 5,575 controls) populaeons, as well as independent EAS individual level data. We evaluated several PRSs approaches: a single-ancestry PRS using 25 variants that reached genome-wide significance (PRS-25), a genome-wide Bayesian based approach (LDpred2), and a mule-ancestry approach that models geneec correlaeons across ancestries (CT-SLEB). PRS performance was evaluated based on the associaeon with LUAD and AUC values. We then esemated the lifeeme absolute risk of LUAD (age 30-80) and projected the AUC at different sample sizes using EAS-derived effect-size distribueon and heritability esemates., Findings: The CT-SLEB PRS showed a strong associaeon with LUAD risk (odds raeo=1.71, 95% confidence interval (CI): 1.61, 1.82) with an AUC of 0.640 (95% CI: 0.629, 0.653). Individuals in the 95 th percenele of the PRS had an esemated 6.69% lifeeme absolute risk of LUAD. Comparison of LUAD risk between individuals in the highest and lowest 20% PRS quaneles revealed a 3.92-fold increase. Projeceon analyses indicated that achieving an AUC of 0.70, which approaches the maximized prediceon poteneal of the PRS given the esemated geneec variance, would require a future study encompassing 55,000 EAS LUAD cases with a 1:10 case-control raeo., Interpretations: Our study underscores the poteneal of mule-ancestry PRS approaches to enhance LUAD risk straeficaeon in never-smokers, parecularly in EAS populaeons, and highlights the necessary scale of future research to uncover the geneec underpinnings of LUAD.

Objectives: The quantitative job-exposure matrix SYN-JEM consists of various dimensions: job-specific estimates, region-specific estimates, and prior expert ratings of jobs by the semi-quantitative DOM-JEM. We analyzed the effect of different JEM dimensions on the exposure-response relationships between occupational silica exposure and lung cancer risk to investigate how these variations influence estimates of exposure by a quantitative JEM and associated health endpoints., Methods: Using SYN-JEM, and alternative SYN-JEM specifications with varying dimensions included, cumulative silica exposure estimates were assigned to 16 901 lung cancer cases and 20 965 controls pooled from 14 international community-based case-control studies. Exposure-response relationships based on SYN-JEM and alternative SYN-JEM specifications were analyzed using regression analyses (by quartiles and log-transformed continuous silica exposure) and generalized additive models (GAM), adjusted for age, sex, study, cigarette pack-years, time since quitting smoking, and ever employment in occupations with established lung cancer risk., Results: SYN-JEM and alternative specifications generated overall elevated and similar lung cancer odds ratios ranging from 1.13 (1st quartile) to 1.50 (4th quartile). In the categorical and log-linear analyses SYN-JEM with all dimensions included yielded the best model fit, and exclusion of job-specific estimates from SYN-JEM yielded the poorest model fit. Additionally, GAM showed the poorest model fit when excluding job-specific estimates., Conclusion: The established exposure-response relationship between occupational silica exposure and lung cancer was marginally influenced by varying the dimensions of SYN-JEM. Optimized modelling of exposure-response relationships will be obtained when incorporating all relevant dimensions, namely prior rating, job, time, and region. Quantitative job-specific estimates appeared to be the most prominent dimension for this general population JEM.

Background: Although the associations between genetic variations and lung cancer risk have been explored, the epigenetic consequences of DNA methylation in lung cancer development are largely unknown. Here, the genetically predicted DNA methylation markers associated with non-small cell lung cancer (NSCLC) risk by a two-stage case-control design were investigated., Methods: The genetic prediction models for methylation levels based on genetic and methylation data of 1595 subjects from the Framingham Heart Study were established. The prediction models were applied to a fixed-effect meta-analysis of screening data sets with 27,120 NSCLC cases and 27,355 controls to identify the methylation markers, which were then replicated in independent data sets with 7844 lung cancer cases and 421,224 controls. Also performed was a multi-omics functional annotation for the identified CpGs by integrating genomics, epigenomics, and transcriptomics and investigation of the potential regulation pathways., Results: Of the 29,894 CpG sites passing the quality control, 39 CpGs associated with NSCLC risk (Bonferroni-corrected p ≤ 1.67 × 10 -6 ) were originally identified. Of these, 16 CpGs remained significant in the validation stage (Bonferroni-corrected p ≤ 1.28 × 10 -3 ), including four novel CpGs. Multi-omics functional annotation showed nine of 16 CpGs were potentially functional biomarkers for NSCLC risk. Thirty-five genes within a 1-Mb window of 12 CpGs that might be involved in regulatory pathways of NSCLC risk were identified., Conclusions: Sixteen promising DNA methylation markers associated with NSCLC were identified. Changes of the methylation level at these CpGs might influence the development of NSCLC by regulating the expression of genes nearby., Plain Language Summary: The epigenetic consequences of DNA methylation in lung cancer development are still largely unknown. This study used summary data of large-scale genome-wide association studies to investigate the associations between genetically predicted levels of methylation biomarkers and non-small cell lung cancer risk at the first time. This study looked at how well larotrectinib worked in adult patients with sarcomas caused by TRK fusion proteins. These findings will provide a unique insight into the epigenetic susceptibility mechanisms of lung cancer., (© 2023 American Cancer Society.)

Cigarette smoke, containing both nicotine and carcinogens, causes lung cancer. However, not all smokers develop lung cancer, highlighting the importance of the interaction between host susceptibility and environmental exposure in tumorigenesis. Here, we aimed to delineate the interaction between metabolizing ability of tobacco carcinogens and smoking intensity in mediating genetic susceptibility to smoking-related lung tumorigenesis. Single-variant and gene-based associations of 43 tobacco carcinogen-metabolizing genes with lung cancer were analyzed using summary statistics and individual-level genetic data, followed by causal inference of Mendelian randomization, mediation analysis, and structural equation modeling. Cigarette smoke-exposed cell models were used to detect gene expression patterns in relation to specific alleles. Data from the International Lung Cancer Consortium (29,266 cases and 56,450 controls) and UK Biobank (2,155 cases and 376,329 controls) indicated that the genetic variant rs56113850 C>T located in intron 4 of CYP2A6 was significantly associated with decreased lung cancer risk among smokers (OR = 0.88, 95% confidence interval = 0.85-0.91, P = 2.18 × 10-16), which might interact (Pinteraction = 0.028) with and partially be mediated (ORindirect = 0.987) by smoking status. Smoking intensity accounted for 82.3% of the effect of CYP2A6 activity on lung cancer risk but entirely mediated the genetic effect of rs56113850. Mechanistically, the rs56113850 T allele rescued the downregulation of CYP2A6 caused by cigarette smoke exposure, potentially through preferential recruitment of transcription factor helicase-like transcription factor. Together, this study provides additional insights into the interplay between host susceptibility and carcinogen exposure in smoking-related lung tumorigenesis., Significance: The causal pathway connecting CYP2A6 genetic variability and activity, cigarette consumption, and lung cancer susceptibility in smokers highlights the need for behavior modification interventions based on host susceptibility for cancer prevention., (©2023 American Association for Cancer Research.)

Background: Although polygenic risk score (PRS) has emerged as a promising tool for predicting cancer risk from genome-wide association studies (GWAS), the individual-level accuracy of lung cancer PRS and the extent to which its impact on subsequent clinical applications remains largely unexplored., Methods: Lung cancer PRSs and confidence/credible interval (CI) were constructed using two statistical approaches for each individual: (1) the weighted sum of 16 GWAS-derived significant SNP loci and the CI through the bootstrapping method (PRS-16-CV) and (2) LDpred2 and the CI through posteriors sampling (PRS-Bayes), among 17,166 lung cancer cases and 12,894 controls with European ancestry from the International Lung Cancer Consortium. Individuals were classified into different genetic risk subgroups based on the relationship between their own PRS mean/PRS CI and the population level threshold., Results: Considerable variances in PRS point estimates at the individual level were observed for both methods, with an average standard deviation (s.d.) of 0.12 for PRS-16-CV and a much larger s.d. of 0.88 for PRS-Bayes. Using PRS-16-CV, only 25.0% of individuals with PRS point estimates in the lowest decile of PRS and 16.8% in the highest decile have their entire 95% CI fully contained in the lowest and highest decile, respectively, while PRS-Bayes was unable to find any eligible individuals. Only 19% of the individuals were concordantly identified as having high genetic risk (> 90th percentile) using the two PRS estimators. An increased relative risk of lung cancer comparing the highest PRS percentile to the lowest was observed when taking the CI into account (OR = 2.73, 95% CI: 2.12-3.50, P-value = 4.13 × 10 -15 ) compared to using PRS-16-CV mean (OR = 2.23, 95% CI: 1.99-2.49, P-value = 5.70 × 10 -46 ). Improved risk prediction performance with higher AUC was consistently observed in individuals identified by PRS-16-CV CI, and the best performance was achieved by incorporating age, gender, and detailed smoking pack-years (AUC: 0.73, 95% CI = 0.72-0.74)., Conclusions: Lung cancer PRS estimates using different methods have modest correlations at the individual level, highlighting the importance of considering individual-level uncertainty when evaluating the practical utility of PRS., (© 2024. The Author(s).)

Rationale: Benzene has been classified as carcinogenic to humans, but there is limited evidence linking benzene exposure to lung cancer. Objectives: We aimed to examine the relationship between occupational benzene exposure and lung cancer. Methods: Subjects from 14 case-control studies across Europe and Canada were pooled. We used a quantitative job-exposure matrix to estimate benzene exposure. Logistic regression models assessed lung cancer risk across different exposure indices. We adjusted for smoking and five main occupational lung carcinogens and stratified analyses by smoking status and lung cancer subtypes. Measurements and Main Results: Analyses included 28,048 subjects (12,329 cases, 15,719 control subjects). Lung cancer odds ratios ranged from 1.12 (95% confidence interval, 1.03-1.22) to 1.32 (95% confidence interval, 1.18-1.48) ( P trend  = 0.002) for groups with the lowest and highest cumulative occupational exposures, respectively, compared with unexposed subjects. We observed an increasing trend of lung cancer with longer duration of exposure ( P trend  < 0.001) and a decreasing trend with longer time since last exposure ( P trend  = 0.02). These effects were seen for all lung cancer subtypes, regardless of smoking status, and were not influenced by specific occupational groups, exposures, or studies. Conclusions: We found consistent and robust associations between different dimensions of occupational benzene exposure and lung cancer after adjusting for smoking and main occupational lung carcinogens. These associations were observed across different subgroups, including nonsmokers. Our findings support the hypothesis that occupational benzene exposure increases the risk of developing lung cancer. Consequently, there is a need to revisit published epidemiological and molecular data on the pulmonary carcinogenicity of benzene.

Background: While much research has been done to identify individual workplace lung carcinogens, little is known about joint effects on risk when workers are exposed to multiple agents., Objectives: We investigated the pairwise joint effects of occupational exposures to asbestos, respirable crystalline silica, metals (i.e., nickel, chromium-VI), and polycyclic aromatic hydrocarbons (PAH) on lung cancer risk, overall and by major histologic subtype, while accounting for cigarette smoking., Methods: In the international 14-center SYNERGY project, occupational exposures were assigned to 16,901 lung cancer cases and 20,965 control subjects using a quantitative job-exposure matrix (SYN-JEM). Odds ratios (ORs) and 95% confidence intervals (CIs) were computed for ever vs. never exposure using logistic regression models stratified by sex and adjusted for study center, age, and smoking habits. Joint effects among pairs of agents were assessed on multiplicative and additive scales, the latter by calculating the relative excess risk due to interaction (RERI)., Results: All pairwise joint effects of lung carcinogens in men were associated with an increased risk of lung cancer. However, asbestos/metals and metals/PAH resulted in less than additive effects; while the chromium-VI/silica pair showed marginally synergistic effect in relation to adenocarcinoma (RERI: 0.24; CI: 0.02, 0.46; p = 0.05). In women, several pairwise joint effects were observed for small cell lung cancer including exposure to PAH/silica (OR = 5.12; CI: 1.77, 8.48), and to asbestos/silica (OR = 4.32; CI: 1.35, 7.29), where exposure to PAH/silica resulted in a synergistic effect (RERI: 3.45; CI: 0.10, 6.8)., Discussion: Small or no deviation from additive or multiplicative effects was observed, but co-exposure to the selected lung carcinogens resulted generally in higher risk than exposure to individual agents, highlighting the importance to reduce and control exposure to carcinogens in workplaces and the general environment. https://doi.org/10.1289/EHP13380.

Observational studies have suggested that smoking may increase the risk of cutaneous squamous cell carcinoma (cSCC) while decreasing the risks of basal cell carcinoma (BCC), and melanoma. However, it remains possible that confounding by other factors may explain these associations. The aim of this investigation was to use Mendelian randomization (MR) to test whether smoking is associated with skin cancer, independently of other factors. Two-sample MR analyses were conducted to determine the causal effect of smoking measures on skin cancer risk using genome-wide association study (GWAS) summary statistics. We used the inverse-variance-weighted estimator to derive separate risk estimates across genetic instruments for all smoking measures. A genetic predisposition to smoking initiation was associated with lower risks of all skin cancer types, although none of the effect estimates reached statistical significance (OR 95% CI BCC 0.91, 0.82-1.01; cSCC 0.82, 0.66-1.01; melanoma 0.91, 0.82-1.01). Results for other measures were similar to smoking initiation with the exception of smoking intensity which was associated with a significantly reduced risk of melanoma (OR 0.67, 95% CI 0.51-0.89). Our findings support the findings of observational studies linking smoking to lower risks of melanoma and BCC. However, we found no evidence that smoking is associated with an elevated risk of cSCC; indeed, our results are most consistent with a decreased risk, similar to BCC and melanoma., (© 2023. Springer Nature Limited.)

Male-pattern baldness (MPB) is related to dysregulation of androgens such as testosterone. A previously observed relationship between MPB and skin cancer may be due to greater exposure to ultraviolet radiation or indicate a role for androgenic pathways in the pathogenesis of skin cancers. We dissected this relationship via Mendelian randomization (MR) analyses, using genetic data from recent male-only meta-analyses of cutaneous melanoma (12,232 cases; 20,566 controls) and keratinocyte cancers (KCs) (up to 17,512 cases; >100,000 controls), followed by stratified MR analysis by body-sites. We found strong associations between MPB and the risk of KC, but not with androgens, and multivariable models revealed that this relationship was heavily confounded by MPB single nucleotide polymorphisms involved in pigmentation pathways. Site-stratified MR analyses revealed strong associations between MPB with head and neck squamous cell carcinoma and melanoma, suggesting that sun exposure on the scalp, rather than androgens, is the main driver. Men with less hair covering likely explains, at least in part, the higher incidence of melanoma in men residing in countries with high ambient UV., (© 2023. Springer Nature Limited.)

Background: The association between duration of smoking abstinence before non-small-cell lung cancer (NSCLC) diagnosis and subsequent survival can influence public health messaging delivered in lung-cancer screening. We aimed to assess whether the duration of smoking abstinence before diagnosis of NSCLC is associated with improved survival., Methods: In this retrospective, pooled analysis of cohort studies, we used 26 cohorts participating in Clinical Outcomes Studies of the International Lung Cancer Consortium (COS-ILCCO) at 23 hospitals. 16 (62%) were from North America, six (23%) were from Europe, three (12%) were from Asia, and one (4%) was from South America. Patients enrolled were diagnosed between June 1, 1983, and Dec 31, 2019. Eligible patients had smoking data before NSCLC diagnosis, epidemiological data at diagnosis (obtained largely from patient questionnaires), and clinical information (retrieved from medical records). Kaplan-Meier curves and multivariable Cox models (ie, adjusted hazard ratios [aHRs]) were generated with individual, harmonised patient data from the consortium database. We estimated overall survival for all causes, measured in years from diagnosis date until the date of the last follow-up or death due to any cause and NSCLC-specific survival., Findings: Of 42 087 patients with NSCLC in the COS-ILCCO database, 21 893 (52·0%) of whom were male and 20 194 (48·0%) of whom were female, we excluded 4474 (10·6%) with missing data. Compared with current smokers (15 036 [40·0%] of 37 613), patients with 1-3 years of smoking abstinence before NSCLC diagnosis (2890 [7·7%]) had an overall survival aHR of 0·92 (95% CI 0·87-0·97), patients with 3-5 years of smoking abstinence (1114 [3·0%]) had an overall survival aHR of 0·90 (0·83-0·97), and patients with more than 5 years of smoking abstinence (10 841 [28·8%]) had an overall survival aHR of 0·90 (0·87-0·93). Improved NSCLC-specific survival was observed in 4301 (44%) of 9727 patients who had quit cigarette smoking and was significant at abstinence durations of more than 5 years (aHR 0·87, 95% CI 0·81-0·93). Results were consistent across age, sex, histology, and disease-stage distributions., Interpretation: In this large, pooled analysis of cohort studies across Asia, Europe, North America, and South America, overall survival was improved in patients with NSCLC whose duration of smoking abstinence before diagnosis was as short as 1 year. These findings suggest that quitting smoking can improve overall survival, even if NSCLC is diagnosed at a later lung-cancer screening visit. These findings also support the implementation of public health smoking cessation strategies at any time., Funding: The Alan B Brown Chair, The Posluns Family Fund, The Lusi Wong Fund, and the Princess Margaret Cancer Foundation., Competing Interests: Declaration of interests We declare no competing interests., (Copyright © 2023 The Author(s). Published by Elsevier Ltd. This is an Open Access article under the CC BY-NC-ND 4.0 license. Published by Elsevier Ltd.. All rights reserved.)

Objectives: Workers in the construction industry have been exposed to asbestos in various occupations. In Italy, a National Mesothelioma Registry has been implemented more than 20 years ago. Using cases selected from this registry and exploiting existing control data sets, we estimated relative risks for pleural mesothelioma (PM) among construction workers., Design: Case-control study., Setting: Cases from the National Mesothelioma Registry (2000-2018), controls from three previous case-control studies., Methods: We selected male PM incident cases diagnosed in 2000-2018. Population controls were taken from three studies performed in six Italian regions within two periods (2002-2004 and 2012-2016). Age-adjusted and period-adjusted unconditional logistic regression models were fitted to estimate odds ratios (OR) for occupations in the construction industry. We followed two approaches, one (primary) excluding and the other (secondary) including subjects employed in other non-construction blue collar occupations for >5 years. For both approaches, we performed an overall analysis including all cases and, given the incomplete temporal and geographic overlap of cases and controls, three time or/and space restricted sensitivity analyses., Results: The whole data set included 15 592 cases and 2210 controls. With the primary approach (4797 cases and 1085 controls), OR was 3.64 (2181 cases) for subjects ever employed in construction. We found elevated risks for blue-collar occupations (1993 cases, OR 4.52), including bricklayers (988 cases, OR 7.05), general construction workers (320 cases, OR 4.66), plumbers and pipe fitters (305 cases, OR 9.13), painters (104 cases, OR 2.17) and several others. Sensitivity analyses yielded very similar findings. Using the secondary approach, we observed similar patterns, but ORs were remarkably lower., Conclusions: We found markedly increased PM risks for most occupations in the construction industry. These findings are relevant for compensation of subjects affected with mesothelioma in the construction industry., Competing Interests: Competing interests: DCo, SMa and CM served as consultants in trials concerning asbestos-related diseases., (© Author(s) (or their employer(s)) 2023. Re-use permitted under CC BY-NC. No commercial re-use. See rights and permissions. Published by BMJ.)

Background: The shared inherited genetic contribution to risk of different cancers is not fully known. In this study, we leverage results from 12 cancer genome-wide association studies (GWAS) to quantify pairwise genome-wide genetic correlations across cancers and identify novel cancer susceptibility loci., Methods: We collected GWAS summary statistics for 12 solid cancers based on 376 759 participants with cancer and 532 864 participants without cancer of European ancestry. The included cancer types were breast, colorectal, endometrial, esophageal, glioma, head and neck, lung, melanoma, ovarian, pancreatic, prostate, and renal cancers. We conducted cross-cancer GWAS and transcriptome-wide association studies to discover novel cancer susceptibility loci. Finally, we assessed the extent of variant-specific pleiotropy among cancers at known and newly identified cancer susceptibility loci., Results: We observed widespread but modest genome-wide genetic correlations across cancers. In cross-cancer GWAS and transcriptome-wide association studies, we identified 15 novel cancer susceptibility loci. Additionally, we identified multiple variants at 77 distinct loci with strong evidence of being associated with at least 2 cancer types by testing for pleiotropy at known cancer susceptibility loci., Conclusions: Overall, these results suggest that some genetic risk variants are shared among cancers, though much of cancer heritability is cancer-specific and thus tissue-specific. The increase in statistical power associated with larger sample sizes in cross-disease analysis allows for the identification of novel susceptibility regions. Future studies incorporating data on multiple cancer types are likely to identify additional regions associated with the risk of multiple cancer types., (© The Author(s) 2023. Published by Oxford University Press. All rights reserved. For permissions, please email: journals.permissions@oup.com.)

Lung adenocarcinoma is the most common type of lung cancer. Known risk variants explain only a small fraction of lung adenocarcinoma heritability. Here, we conducted a two-stage genome-wide association study of lung adenocarcinoma of East Asian ancestry (21,658 cases and 150,676 controls; 54.5% never-smokers) and identified 12 novel susceptibility variants, bringing the total number to 28 at 25 independent loci. Transcriptome-wide association analyses together with colocalization studies using a Taiwanese lung expression quantitative trait loci dataset (n = 115) identified novel candidate genes, including FADS1 at 11q12 and ELF5 at 11p13. In a multi-ancestry meta-analysis of East Asian and European studies, four loci were identified at 2p11, 4q32, 16q23, and 18q12. At the same time, most of our findings in East Asian populations showed no evidence of association in European populations. In our studies drawn from East Asian populations, a polygenic risk score based on the 25 loci had a stronger association in never-smokers vs. individuals with a history of smoking (P interaction  = 0.0058). These findings provide new insights into the etiology of lung adenocarcinoma in individuals from East Asian populations, which could be important in developing translational applications., (© 2023. This is a U.S. Government work and not under copyright protection in the US; foreign copyright protection may apply.)

There is limited evidence regarding the exposure-effect relationship between lung-cancer risk and hexavalent chromium (Cr(VI)) or nickel. We estimated lung-cancer risks in relation to quantitative indices of occupational exposure to Cr(VI) and nickel and their interaction with smoking habits. We pooled 14 case-control studies from Europe and Canada, including 16 901 lung-cancer cases and 20 965 control subjects. A measurement-based job-exposure-matrix estimated job-year-region specific exposure levels to Cr(VI) and nickel, which were linked to the subjects' occupational histories. Odds ratios (OR) and associated 95% confidence intervals (CI) were calculated by unconditional logistic regression, adjusting for study, age group, smoking habits and exposure to other occupational lung carcinogens. Due to their high correlation, we refrained from mutually adjusting for Cr(VI) and nickel independently. In men, ORs for the highest quartile of cumulative exposure to CR(VI) were 1.32 (95% CI 1.19-1.47) and 1.29 (95% CI 1.15-1.45) in relation to nickel. Analogous results among women were: 1.04 (95% CI 0.48-2.24) and 1.29 (95% CI 0.60-2.86), respectively. In men, excess lung-cancer risks due to occupational Cr(VI) and nickel exposure were also observed in each stratum of never, former and current smokers. Joint effects of Cr(VI) and nickel with smoking were in general greater than additive, but not different from multiplicative. In summary, relatively low cumulative levels of occupational exposure to Cr(VI) and nickel were associated with increased ORs for lung cancer, particularly in men. However, we cannot rule out a combined classical measurement and Berkson-type of error structure, which may cause differential bias of risk estimates., (© 2022 The Authors. International Journal of Cancer published by John Wiley & Sons Ltd on behalf of UICC.)

Background: Spitzoid morphology in familial melanoma has been associated with germline variants in POT1 , a telomere maintenance gene (TMG), suggesting a link between telomere biology and spitzoid differentiation., Objective: To assess if familial melanoma cases associated with germline variants in TMG ( POT1 , ACD , TERF2IP , and TERT ) commonly exhibit spitzoid morphology., Methods: In this case series, melanomas were classified as having spitzoid morphology if at least 3 of 4 dermatopathologists reported this finding in ≥25% of tumor cells. Logistic regression was used to calculate odds ratios (OR) of spitzoid morphology compared to familial melanomas from unmatched noncarriers that were previously reviewed by a National Cancer Institute dermatopathologist., Results: Spitzoid morphology was observed in 77% (23 of 30), 75% (3 of 4), 50% (2 of 4), and 50% (1 of 2) of melanomas from individuals with germline variants in POT1 , TERF2IP , ACD , and TERT , respectively. Compared to noncarriers ( n  = 139 melanomas), POT1 carriers (OR = 225.1, 95% confidence interval: 51.7-980.5; P  < .001) and individuals with TERF2IP, ACD, and TERT variants (OR = 82.4, 95% confidence interval: 21.3-494.6; P  < .001) had increased odds of spitzoid morphology., Limitations: Findings may not be generalizable to nonfamilial melanoma cases., Conclusion: Spitzoid morphology in familial melanoma could suggest germline alteration of TMG., Competing Interests: None disclosed.

Background: Data are scarce about tumor mutational burden (TMB) as a biomarker in never smokers with non-small cell lung cancer (NSCLC)., Methods: TMB was assessed by whole-genome sequencing (WGS) and compared with in silico reduced whole-exome sequencing (WES) and targeted commercial next-generation sequencing (NGS) gene panels in 92 paired tumor-normal samples from never smokers who underwent NSCLC resection with curative intent. Analyses were performed to test for association with survival after surgery and to identify the optimal prognostic TMB cutoff., Results: Tumors of never smokers with NSCLC had low TMB scores (median 1.57 mutations/Mb; range, 0.13-17.94). A TMB cutoff of 1.70 mutations/Mb was associated with a 5-year overall survival of 58% in the high-TMB (42% of cases) compared with 86% in low-TMB patients (Wald P = 0.0029). TMB scores from WGS and WES were highly correlated (Spearman ρ = 0.93, P < 2.2e-16). TMB scores from NGS panels demonstrated high intraindividual fluctuations and identified high-TMB patients with 65% concordance in average compared with WGS., Conclusions: In resected NSCLC of never smokers, high TMB was associated with worse prognosis. WES provided a good estimate of TMB while targeted NGS panels seem to lack adequate depth and resolution in the setting of low mutation burden., Impact: TMB is a prognostic indicator of survival in resected NSCLC from individuals who never smoked. In this setting of low mutation counts, TMB can be accurately measured by WGS or WES, but not NGS panels., (©2022 The Authors; Published by the American Association for Cancer Research.)

The most recent genome-wide association study (GWAS) of cutaneous melanoma identified 54 risk-associated loci, but functional variants and their target genes for most have not been established. Here, we performed massively parallel reporter assays (MPRAs) by using malignant melanoma and normal melanocyte cells and further integrated multi-layer annotation to systematically prioritize functional variants and susceptibility genes from these GWAS loci. Of 1,992 risk-associated variants tested in MPRAs, we identified 285 from 42 loci (78% of the known loci) displaying significant allelic transcriptional activities in either cell type (FDR < 1%). We further characterized MPRA-significant variants by motif prediction, epigenomic annotation, and statistical/functional fine-mapping to create integrative variant scores, which prioritized one to six plausible candidate variants per locus for the 42 loci and nominated a single variant for 43% of these loci. Overlaying the MPRA-significant variants with genome-wide significant expression or methylation quantitative trait loci (eQTLs or meQTLs, respectively) from melanocytes or melanomas identified candidate susceptibility genes for 60% of variants (172 of 285 variants). CRISPRi of top-scoring variants validated their cis-regulatory effect on the eQTL target genes, MAFF (22q13.1) and GPRC5A (12p13.1). Finally, we identified 36 melanoma-specific and 45 melanocyte-specific MPRA-significant variants, a subset of which are linked to cell-type-specific target genes. Analyses of transcription factor availability in MPRA datasets and variant-transcription-factor interaction in eQTL datasets highlighted the roles of transcription factors in cell-type-specific variant functionality. In conclusion, MPRAs along with variant scoring effectively prioritized plausible candidates for most melanoma GWAS loci and highlighted cellular contexts where the susceptibility variants are functional., Competing Interests: Declaration of interests The authors declare no competing interests., (Published by Elsevier Inc.)

Mutational signature analysis is commonly performed in cancer genomic studies. Here, we present SigProfilerExtractor, an automated tool for de novo extraction of mutational signatures, and benchmark it against another 13 bioinformatics tools by using 34 scenarios encompassing 2,500 simulated signatures found in 60,000 synthetic genomes and 20,000 synthetic exomes. For simulations with 5% noise, reflecting high-quality datasets, SigProfilerExtractor outperforms other approaches by elucidating between 20% and 50% more true-positive signatures while yielding 5-fold less false-positive signatures. Applying SigProfilerExtractor to 4,643 whole-genome- and 19,184 whole-exome-sequenced cancers reveals four novel signatures. Two of the signatures are confirmed in independent cohorts, and one of these signatures is associated with tobacco smoking. In summary, this report provides a reference tool for analysis of mutational signatures, a comprehensive benchmarking of bioinformatics tools for extracting signatures, and several novel mutational signatures, including one putatively attributed to direct tobacco smoking mutagenesis in bladder tissues., Competing Interests: M.V. is an employee of NVIDIA corporation. L.B.A. is a compensated consultant and has equity interest in io9, LLC. His spouse is an employee of Biotheranostics, Inc. L.B.A. and B.S.A. are inventors of a US patent 10,776,718. E.N.B. and L.B.A. declare US provisional applications with serial numbers 63/289,601 and 63/269,033. L.B.A. and A.A. declare US provisional patent applications with serial numbers 63/366,392 and 63/367,846. All other authors declare no competing interests., (© 2022 The Author(s).)

At a time when complex diseases affect globally 280 million people and claim 14 million lives every year, there is an urgent need to rapidly increase our knowledge into their underlying etiologies. Though critical in identifying the people at risk, the causal environmental factors (microbiome and/or pollutants) and the affected pathophysiological mechanisms are not well understood. Herein, we consider the variations of autophagy-related ( ATG ) genes at the heart of mechanisms of increased susceptibility to environmental stress. A comprehensive autophagy genomic resource is presented with 263 single nucleotide polymorphisms (SNPs) for 69 autophagy-related genes associated with 117 autoimmune, inflammatory, infectious, cardiovascular, neurological, respiratory, and endocrine diseases. We thus propose the term 'autophagopathies' to group together a class of complex human diseases the etiology of which lies in a genetic defect of the autophagy machinery, whether directly related or not to an abnormal flux in autophagy, LC3-associated phagocytosis, or any associated trafficking. The future of precision medicine for common diseases will lie in our ability to exploit these ATG SNP x environment relationships to develop new polygenetic risk scores, new management guidelines, and optimal therapies for afflicted patients. Abbreviations: ATG, autophagy-related; ALS-FTD, amyotrophic lateral sclerosis-frontotemporal dementia; ccRCC, clear cell renal cell carcinoma; CD, Crohn disease; COPD, chronic obstructive pulmonary disease; eQTL, expression quantitative trait loci; HCC, hepatocellular carcinoma; HNSCC, head and neck squamous cell carcinoma; GTEx, genotype-tissue expression; GWAS, genome-wide association studies; LAP, LC3-associated phagocytosis; LC3-II, phosphatidylethanolamine conjugated form of LC3; LD, linkage disequilibrium; LUAD, lung adenocarcinoma; MAF, minor allele frequency; MAP1LC3/LC3: microtubule associated protein 1 light chain 3; NSCLC, non-small cell lung cancer; OS, overall survival; PtdIns3K CIII, class III phosphatidylinositol 3 kinase; PtdIns3P, phosphatidylinositol-3-phosphate; SLE, systemic lupus erythematosus; SNPs, single-nucleotide polymorphisms; mQTL, methylation quantitative trait loci; ULK, unc-51 like autophagy activating kinase; UTRs, untranslated regions; WHO, World Health Organization.

Purpose: Pancreatic ductal adenocarcinoma (PDAC) is a component of familial melanoma due to germline pathogenic variants (GPVs) in CDKN2A . However, it is unclear what role this gene or other genes play in its etiology., Materials and Methods: We analyzed 189 cancer predisposition genes using parametric rare-variant association (RVA) tests and nonparametric permutation tests to identify gene-level associations in PDAC for patients with ( CDKN2A + ) and without ( CDKN2A- ) GPV. Exome sequencing was performed on 84 patients with PDAC, 47 CDKN2A+ and 37 CDKN2A- . After variant filtering, various RVA tests and permutation tests were run separately by CDKN2A status. Genes with the strongest nominal associations were evaluated in patients with PDAC from The Cancer Genome Atlas and the UK Biobank (UKB). A secondary analysis including only GPV from UKB was also performed., Results: In RVA tests, ERCC4 and RET showed the most compelling evidence as plausible PDAC candidate genes for CDKN2A+ patients. In contrast, the findings in CDKN2A- patients provided evidence for HMBS , EPCAM , and MRE11 as potential new candidate genes and confirmed ATM, BRCA2 , and PALB2 as PDAC genes, consistent with findings in The Cancer Genome Atlas and the UKB. As expected, CDKN2A- patients were more likely to harbor GPVs from the 189 genes investigated. When including only GPVs from UKB, significant associations with PDAC were seen for ATM, BRCA2, and CDKN2A ., Conclusion: These results suggest that variants in other genes likely play a role in PDAC in all patients and that PDAC in CDKN2A+ patients has a distinct etiology from PDAC in CDKN2A- patients.

Although multiple common susceptibility loci for lung cancer (LC) have been identified by genome-wide association studies, they can explain only a small portion of heritability. The etiological contribution of rare deleterious variants (RDVs) to LC risk is not fully characterized and may account for part of the missing heritability. Here, we sequenced the whole exomes of 2777 participants from the Environment and Genetics in Lung cancer Etiology study, a homogenous population including 1461 LC cases and 1316 controls. In single-variant analyses, we identified a new RDV, rs77187983 [EHBP1, odds ratio (OR) = 3.13, 95% confidence interval (CI) = 1.34-7.30, P = 0.008] and replicated two previously reported RDVs, rs11571833 (BRCA2, OR = 2.18; 95% CI = 1.25-3.81, P = 0.006) and rs752672077 (MPZL2, OR = 3.70, 95% CI = 1.04-13.15, P = 0.044). In gene-based analyses, we confirmed BRCA2 (P = 0.007) and ATM (P = 0.014) associations with LC risk and identified TRIB3 (P = 0.009), involved in maintaining genome stability and DNA repair, as a new candidate susceptibility gene. Furthermore, cases were enriched with RDVs in homologous recombination repair [carrier frequency (CF) = 22.9% versus 19.5%, P = 0.017] and Fanconi anemia (CF = 12.5% versus 10.2%, P = 0.036) pathways. Our results were not significant after multiple testing corrections but were enriched in cases versus controls from large scale public biobank resources, including The Cancer Genome Atlas, FinnGen and UK Biobank. Our study identifies novel candidate genes and highlights the importance of RDVs in DNA repair-related genes for LC susceptibility. These findings improve our understanding of LC heritability and may contribute to the development of risk stratification and prevention strategies., (Published by Oxford University Press 2022.)

Background: The role of germline genetic factors in determining survival from cutaneous melanoma (CM) is not well understood., Objective: To perform a genome-wide association study (GWAS) meta-analysis of melanoma-specific survival (MSS), and test whether a CM-susceptibility polygenic risk score (PRS) is associated with MSS., Methods: We conducted two Cox proportional-hazard GWAS of MSS using data from the Melanoma Institute Australia, a high ultraviolet (UV) radiation setting (MIA; 5,762 patients with melanoma; 800 melanoma deaths) and UK Biobank (UKB: 5,220 patients with melanoma; 241 melanoma deaths), and combined them in a fixed-effects meta-analysis. Significant (P < 5 × 10-8) results were investigated in the Leeds Melanoma Cohort (LMC; 1,947 patients with melanoma; 370 melanoma deaths). We also developed a CM-susceptibility PRS using a large independent GWAS meta-analysis (23,913 cases, 342,870 controls). The PRS was tested for an association with MSS in the MIA and UKB cohorts., Results: Two loci were significantly associated with MSS in the meta-analysis of MIA and UKB with lead SNPs rs41309643 (G allele frequency 1.6%, HR = 2.09, 95%CI = 1.61-2.71, P = 2.08 × 10-8) on chromosome 1, and rs75682113 (C allele frequency 1.8%, HR = 2.38, 95%CI = 1.77-3.21, P = 1.07 × 10-8) on chromosome 7. While neither SNP replicated in the LMC, rs75682113 was significantly associated in the combined discovery and replication sets. After adjusting for age at diagnosis, sex and the first ten principal components, a one standard deviation increase in the CM-susceptibility PRS was associated with improved MSS in the discovery meta-analysis (HR = 0.88, 95% CI = 0.83-0.94, P = 6.93 × 10-5; I2 = 88%). However, this was only driven by the high UV setting cohort (MIA HR = 0.84, 95% CI = 0.78-0.90)., Conclusion: We found two loci potentially associated with MSS. Increased genetic susceptibility to develop CM is associated with improved MSS in a high UV setting., (© 2022. The Author(s).)

Cancer heterogeneities hold the key to a deeper understanding of cancer etiology and progression and the discovery of more precise cancer therapy. Modern pathological and molecular technologies offer a powerful set of tools to profile tumor heterogeneities at multiple levels in large patient populations, from DNA to RNA, protein and epigenetics, and from tumor tissues to tumor microenvironment and liquid biopsy. When coupled with well-validated epidemiologic methodology and well-characterized epidemiologic resources, the rich tumor pathological and molecular tumor information provide new research opportunities at an unprecedented breadth and depth. This is the research space where Molecular Pathological Epidemiology (MPE) emerged over a decade ago and has been thriving since then. As a truly multidisciplinary field, MPE embraces collaborations from diverse fields including epidemiology, pathology, immunology, genetics, biostatistics, bioinformatics, and data science. Since first convened in 2013, the International MPE Meeting series has grown into a dynamic and dedicated platform for experts from these disciplines to communicate novel findings, discuss new research opportunities and challenges, build professional networks, and educate the next-generation scientists. Herein, we share the proceedings of the Fifth International MPE meeting, held virtually online, on May 24 and 25, 2021. The meeting consisted of 21 presentations organized into the three main themes, which were recent integrative MPE studies, novel cancer profiling technologies, and new statistical and data science approaches. Looking forward to the near future, the meeting attendees anticipated continuous expansion and fruition of MPE research in many research fronts, particularly immune-epidemiology, mutational signatures, liquid biopsy, and health disparities., (© 2022. The Author(s), under exclusive licence to Springer Nature Switzerland AG.)

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., rs521828 C6orf10 and rs204999 PRRT1 , OR interaction  = 1.17, p = 6.57 × 10 -13 ; rs3135369 BTNL2 and rs2858859 HLA-DQA1 , OR interaction  = 1.17, p = 2.43 × 10 -13 ; rs2858859 HLA-DQA1 and rs9275572 HLA-DQA2 , OR interaction  = 1.15, p = 2.84 × 10 -13 ; rs2853668 TERT and rs62329694 CLPTM1L , OR interaction  = 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 (rs521828 C6orf10 and rs204999 PRRT1 , OR interaction  = 1.13, p = 0.008; rs3135369 BTNL2 and rs2858859 HLA-DQA1 , OR interaction  = 1.11, p = 5.23 × 10 -4 ; rs3135369 BTNL2 and rs9271300 HLA-DQA1 , OR interaction  = 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., (Copyright © 2022 International Association for the Study of Lung Cancer. Published by Elsevier Inc. All rights reserved.)

To identify new susceptibility loci to lung cancer among diverse populations, we performed cross-ancestry genome-wide association studies in European, East Asian and African populations and discovered five loci that have not been previously reported. We replicated 26 signals and identified 10 new lead associations from previously reported loci. Rare-variant associations tended to be specific to populations, but even common-variant associations influencing smoking behavior, such as those with CHRNA5 and CYP2A6, showed population specificity. Fine-mapping and expression quantitative trait locus colocalization nominated several candidate variants and susceptibility genes such as IRF4 and FUBP1. DNA damage assays of prioritized genes in lung fibroblasts indicated that a subset of these genes, including the pleiotropic gene IRF4, potentially exert effects by promoting endogenous DNA damage., (© 2022. The Author(s), under exclusive licence to Springer Nature America, Inc.)

The microbiome is the collection of all microbial genes and can be investigated by sequencing highly variable regions of 16S ribosomal RNA (rRNA) genes. Evidence suggests that environmental factors and host genetics may interact to impact human microbiome composition. Identifying host genetic variants associated with human microbiome composition not only provides clues for characterizing microbiome variation but also helps to elucidate biological mechanisms of genetic associations, prioritize genetic variants, and improve genetic risk prediction. Since a microbiota functions as a community, it is best characterized by β diversity; that is, a pairwise distance matrix. We develop a statistical framework and a computationally efficient software package, microbiomeGWAS, for identifying host genetic variants associated with microbiome β diversity with or without interacting with an environmental factor. We show that the score statistics have positive skewness and kurtosis due to the dependent nature of the pairwise data, which makes p -value approximations based on asymptotic distributions unacceptably liberal. By correcting for skewness and kurtosis, we develop accurate p -value approximations, whose accuracy was verified by extensive simulations. We exemplify our methods by analyzing a set of 147 genotyped subjects with 16S rRNA microbiome profiles from non-malignant lung tissues. Correcting for skewness and kurtosis eliminated the dramatic deviation in the quantile-quantile plots. We provided preliminary evidence that six established lung cancer risk SNPs were collectively associated with microbiome composition for both unweighted ( p = 0.0032) and weighted ( p = 0.011) UniFrac distance matrices. In summary, our methods will facilitate analyzing large-scale genome-wide association studies of the human microbiome.

Background: Exposure to polycyclic aromatic hydrocarbons (PAH) occurs widely in occupational settings. We investigated the association between occupational exposure to PAH and lung cancer risk and joint effects with smoking within the SYNERGY project., Methods: We pooled 14 case-control studies with information on lifetime occupational and smoking histories conducted between 1985 and 2010 in Europe and Canada. Exposure to benzo[a]pyrene (BaP) was used as a proxy of PAH and estimated from a quantitative general population job-exposure matrix. Multivariable unconditional logistic regression models, adjusted for smoking and exposure to other occupational lung carcinogens, estimated ORs, and 95% confidence intervals (CI)., Results: We included 16,901 lung cancer cases and 20,965 frequency-matched controls. Adjusted OR for PAH exposure (ever) was 1.08 (CI, 1.02-1.15) in men and 1.20 (CI, 1.04-1.38) in women. When stratified by smoking status and histologic subtype, the OR for cumulative exposure ≥0.24 BaP μg/m3-years in men was higher in never smokers overall [1.31 (CI, 0.98-1.75)], for small cell [2.53 (CI, 1.28-4.99)] and squamous cell cancers [1.33 (CI, 0.80-2.21)]. Joint effects between PAH and smoking were observed. Restricting analysis to the most recent studies showed no increased risk., Conclusions: Elevated lung cancer risk associated with PAH exposure was observed in both sexes, particularly for small cell and squamous cell cancers, after accounting for cigarette smoking and exposure to other occupational lung carcinogens., Impact: The lack of association between PAH and lung cancer in more recent studies merits further research under today's exposure conditions and worker protection measures., (©2022 The Authors; Published by the American Association for Cancer Research.)

While several high-penetrance melanoma risk genes are known, variation in these genes fail to explain melanoma susceptibility in a large proportion of high-risk families. As part of a melanoma family sequencing study, including 435 families from Mediterranean populations we identified a novel NRAS variant (c.170A > C, p.D57A) in an Italian melanoma-prone family. This variant is absent in exomes in gnomAD, ESP, UKBiobank, and the 1000 Genomes Project, as well as in 11,273 Mediterranean individuals and 109 melanoma-prone families from the US and Australia. This variant occurs in the GTP-binding pocket of NRAS. Differently from other RAS activating alterations, NRAS D57A expression is unable to activate MAPK-pathway both constitutively and after stimulation but enhances EGF-induced PI3K-pathway signaling in serum starved conditions in vitro. Consistent with in vitro data demonstrating that NRAS D57A does not enrich GTP binding, molecular modeling suggests that the D57A substitution would be expected to impair Mg2 + binding and decrease nucleotide-binding and GTPase activity of NRAS. While we cannot firmly establish NRAS c.170A > C (p.D57A) as a melanoma susceptibility variant, further investigation of NRAS as a familial melanoma gene is warranted., (© 2021. This is a U.S. government work and not under copyright protection in the U.S.; foreign copyright protection may apply.)

Limited efforts have been made in assessing the effect of genome-wide profiling of RNA splicing-related variation on lung cancer risk. In the present study, we first identified RNA splicing-related genetic variants linked to lung cancer in a genome-wide profiling analysis and then conducted a two-stage (discovery and replication) association study in populations of European ancestry. Discovery and validation were conducted sequentially with a total of 29,266 cases and 56,450 controls from both the Transdisciplinary Research in Cancer of the Lung and the International Lung Cancer Consortium as well as the OncoArray database. For those variants identified as significant in the two datasets, we further performed stratified analyses by smoking status and histological type and investigated their effects on gene expression and potential regulatory mechanisms. We identified three genetic variants significantly associated with lung cancer risk: rs329118 in JADE2 (P = 8.80E-09), rs2285521 in GGA2 (P = 4.43E-08), and rs198459 in MYRF (P = 1.60E-06). The combined effects of all three SNPs were more evident in lung squamous cell carcinomas (P = 1.81E-08, P = 6.21E-08, and P = 7.93E-04, respectively) than in lung adenocarcinomas and in ever smokers (P = 9.80E-05, P = 2.70E-04, and P = 2.90E-05, respectively) than in never smokers. Gene expression quantitative trait analysis suggested a role for the SNPs in regulating transcriptional expression of the corresponding target genes. In conclusion, we report that three RNA splicing-related genetic variants contribute to lung cancer susceptibility in European populations. However, additional validation is needed, and specific splicing mechanisms of the target genes underlying the observed associations also warrants further exploration., (© 2022. The Author(s).)

Genome-wide association studies (GWAS) have identified a number of risk loci for cutaneous melanoma. Cutaneous melanoma shares overlapping genetic risk (genetic correlation) with a number of other traits, including its risk factors such as sunburn propensity. This genetic correlation can be exploited to identify additional cutaneous melanoma risk loci by multitrait analysis of GWAS (MTAG). We used bivariate linkage disequilibrium-score regression score regression to identify traits that are genetically correlated with clinically confirmed cutaneous melanoma and then used publicly available GWAS for these traits in a multitrait analysis of GWAS. Multitrait analysis of GWAS allows GWAS to be combined while accounting for sample overlap and incomplete genetic correlation. We identified a total of 74 genome-wide independent loci, 19 of them were not previously reported in the input cutaneous melanoma GWAS meta-analysis. Of these loci, 55 were replicated (P < 0.05/74, Bonferroni-corrected P-value in two independent cutaneous melanoma replication cohorts from Melanoma Institute Australia and 23andMe, Inc. Among the, to our knowledge, previously unreported cutaneous melanoma loci are ones that have also been associated with autoimmune traits including rs715199 near LPP and rs10858023 near AP4B1. Our analysis indicates genetic correlation between traits can be leveraged to identify new risk genes for cutaneous melanoma., (Copyright © 2021 The Authors. Published by Elsevier Inc. All rights reserved.)

For de novo mutational signature analysis, the critical first step is to decide how many signatures should be expected in a cancer genomics study. An incorrect number could mislead downstream analyses. Here we present SUITOR (Selecting the nUmber of mutatIonal signaTures thrOugh cRoss-validation), an unsupervised cross-validation method that requires little assumptions and no numerical approximations to select the optimal number of signatures without overfitting the data. In vitro studies and in silico simulations demonstrated that SUITOR can correctly identify signatures, some of which were missed by other widely used methods. Applied to 2,540 whole-genome sequenced tumors across 22 cancer types, SUITOR selected signatures with the smallest prediction errors and almost all signatures of breast cancer selected by SUITOR were validated in an independent breast cancer study. SUITOR is a powerful tool to select the optimal number of mutational signatures, facilitating downstream analyses with etiological or therapeutic importance., Competing Interests: The authors have declared that no competing interests exist.

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., (©2022 American Association for Cancer Research.)

Lung cancer is the major leading cause of cancer-related mortality worldwide. Multiple epigenetic factors-in particular, DNA methylation-have been associated with the development of lung cancer. In this review, we summarize the current knowledge on DNA methylation alterations in lung tumorigenesis, as well as their associations with different histological subtypes, common cancer driver gene mutations (e.g., KRAS , EGFR , and TP53 ), and major epidemiological risk factors (e.g., sex, smoking status, race/ethnicity). Understanding the mechanisms of DNA methylation regulation and their associations with various risk factors can provide further insights into carcinogenesis, and create future avenues for prevention and personalized treatments. In addition, we also highlight outstanding questions regarding DNA methylation in lung cancer to be elucidated in future studies.

Background: Immune cell transcriptome signatures have been widely used to study the lung tumor microenvironment (TME). However, it is unclear to what extent the immune cell composition in the lung TME varies across histological and molecular subtypes (intertumor heterogeneity [inter-TH]) and within tumors (intratumor heterogeneity [ITH]) and whether ITH has any prognostic relevance., Methods: Using RNA sequencing in 269 tumor samples from 160 lung cancer patients, we quantified the inter-TH of immune gene expression and immune cell abundance and evaluated the association of median immune cell abundance with clinical and pathological features and overall survival. In 39 tumors with 132 multiregion samples, we also analyzed the ITH of immune cell abundance in relation to overall survival using a variance-weighted multivariate Cox model not biased by the number of samples per tumor., Results: Substantial inter-TH of 14 immune cell types was observed even within the same histological and molecular subtypes, but early stage tumors had higher lymphocyte infiltration across all tumor types. In multiregion samples, an unbiased estimate of low ITH of overall immune cell composition (hazard ratio [HR] = 0.40, 95% confidence interval [CI] = 0.21 to 0.78; P = .007), dendritic cells (HR = 0.24, 95% CI = 0.096 to 0.58; P = .002), and macrophages (HR = 0.50, 95% CI = 0.30 to 0.84; P = .009) was strongly associated with poor survival. The ITH of 3 markers, including CD163 and CD68 (macrophages) and CCL13 (dendritic cells), was enough to characterize the ITH of the corresponding immune cell abundances and its association with overall survival., Conclusion: This study suggests that lack of immune cell diversity may facilitate tumor evasion and progression. ITH inferred from CCL13, CD163, and CD68 could be used as a prognostic tool in clinical practice., (Published by Oxford University Press 2021. This work is written by US Government employees and is in the public domain in the US.)

Clustered somatic mutations are common in cancer genomes and previous analyses reveal several types of clustered single-base substitutions, which include doublet- and multi-base substitutions 1-5 , diffuse hypermutation termed omikli 6 , and longer strand-coordinated events termed kataegis 3,7-9 . Here we provide a comprehensive characterization of clustered substitutions and clustered small insertions and deletions (indels) across 2,583 whole-genome-sequenced cancers from 30 types of cancer 10 . Clustered mutations were highly enriched in driver genes and associated with differential gene expression and changes in overall survival. Several distinct mutational processes gave rise to clustered indels, including signatures that were enriched in tobacco smokers and homologous-recombination-deficient cancers. Doublet-base substitutions were caused by at least 12 mutational processes, whereas most multi-base substitutions were generated by either tobacco smoking or exposure to ultraviolet light. Omikli events, which have previously been attributed to APOBEC3 activity 6 , accounted for a large proportion of clustered substitutions; however, only 16.2% of omikli matched APOBEC3 patterns. Kataegis was generated by multiple mutational processes, and 76.1% of all kataegic events exhibited mutational patterns that are associated with the activation-induced deaminase (AID) and APOBEC3 family of deaminases. Co-occurrence of APOBEC3 kataegis and extrachromosomal DNA (ecDNA), termed kyklonas (Greek for cyclone), was found in 31% of samples with ecDNA. Multiple distinct kyklonic events were observed on most mutated ecDNA. ecDNA containing known cancer genes exhibited both positive selection and kyklonic hypermutation. Our results reveal the diversity of clustered mutational processes in human cancer and the role of APOBEC3 in recurrently mutating and fuelling the evolution of ecDNA., (© 2022. The Author(s).)

Background: Aberrant Wnt signalling, regulating cell development and stemness, influences the development of many cancer types. The Aryl hydrocarbon receptor (AhR) mediates tumorigenesis of environmental pollutants. Complex interaction patterns of genes assigned to AhR/Wnt-signalling were recently associated with lung cancer susceptibility., Aim: To assess the association and predictive ability of AhR/Wnt-genes with lung cancer in cases and controls of European descent., Methods: Odds ratios (OR) were estimated for genomic variants assigned to the Wnt agonist and the antagonistic genes DKK2, DKK3, DKK4, FRZB, SFRP4 and Axin2. Logistic regression models with variable selection were trained, validated and tested to predict lung cancer, at which other previously identified SNPs that have been robustly associated with lung cancer risk could also enter the model. Furthermore, decision trees were created to investigate variant × variant interaction. All analyses were performed for overall lung cancer and for subgroups., Results: No genome-wide significant association of AhR/Wnt-genes with overall lung cancer was observed, but within the subgroups of ever smokers (e.g., maker rs2722278 SFRP4; OR  = 1.20; 95% CI 1.13-1.27; p  = 5.6 × 10 -10 ) and never smokers (e.g., maker rs1133683 Axin2; OR  = 1.27; 95% CI 1.19-1.35; p  = 1.0 × 10 -12 ). Although predictability is poor, AhR/Wnt-variants are unexpectedly overrepresented in optimized prediction scores for overall lung cancer and for small cell lung cancer. Remarkably, the score for never-smokers contained solely two AhR/Wnt-variants. The optimal decision tree for never smokers consists of 7 AhR/Wnt-variants and only two lung cancer variants., Conclusions: The role of variants belonging to Wnt/AhR-pathways in lung cancer susceptibility may be underrated in main-effects association analysis. Complex interaction patterns in individuals of European descent have moderate predictive capacity for lung cancer or subgroups thereof, especially in never smokers., (© 2022. The Author(s).)

The tumor microenvironment (TME), including immune cells, cancer-associated fibroblasts, endothelial cells, adjacent normal cells, and others, plays a crucial role in influencing tumor behavior and progression. Here, we characterized the TME in 83 primary renal tumors and matched metastatic or recurrence tissue samples (n = 15) from papillary renal cell carcinoma (pRCC) types 1 (n = 20) and 2 (n = 49), collecting duct carcinomas (CDC; n = 14), and high-grade urothelial carcinomas (HGUC; n = 5). We investigated 10 different markers of immune infiltration, vasculature, cell proliferation, and epithelial-to-mesenchymal transition by using machine learning image analysis in conjunction with immunohistochemistry. Marker expression was compared by Mann-Whitney and Kruskal-Wallis tests and correlations across markers using Spearman's rank correlation coefficient. Multivariable Poisson regression analysis was used to compare marker expression between histological types, while accounting for variation in tissue size. Several immune markers showed different rates of expression across histological types of renal carcinoma. Using pRCC1 as reference, the incidence rate ratio (IRR) of CD3+ T cells (IRR [95% confidence interval, CI] = 2.48 [1.53-4.01]) and CD20+ B cells (IRR [95% CI] = 4.38 [1.22-5.58]) was statistically significantly higher in CDC. In contrast, CD68+ macrophages predominated in pRCC1 (IRR [95% CI] = 2.35 [1.42-3.9]). Spatial analysis revealed CD3+ T-cell and CD20+ B-cell expressions in CDC to be higher at the proximal (p < 0.0001) and distal (p < 0.0001) tumor periphery than within the central tumor core. In contrast, expression of CD68+ macrophages in pRCC2 was higher in the tumor center compared to the proximal (p = 0.0451) tumor periphery and pRCC1 showed a distance-dependent reduction, from the central tumor, in CD68+ macrophages with the lowest expression of CD68 marker at the distal tumor periphery (p = 0.004). This study provides novel insights into the TME of rare kidney cancer types, which are often understudied. Our findings of differences in marker expression and localization by histological subtype could have implications for tumor progression and response to immunotherapies or other targeted therapies., (© 2021 The Authors. The Journal of Pathology: Clinical Research published by The Pathological Society of Great Britain and Ireland and John Wiley & Sons Ltd. This article has been contributed to by US Government employees and their work is in the public domain in the USA.)

Purpose: Previous studies have shown an approximately two-fold elevation in the relative risk of urinary bladder cancer (UBC) among people with a family history that could not be entirely explained by shared environmental exposures, thus suggesting a genetic component in its predisposition. Multiple genome-wide association studies and recent gene panel sequencing studies identified several genetic loci that are associated with UBC risk; however, the list of UBC-associated variants and genes is incomplete., Materials and Methods: We exome sequenced eight patients from three multiplex UBC pedigrees and a group of 77 unrelated familial UBC cases matched to 241 cancer-free controls. In addition, we examined pathogenic germline variation in 444 candidate genes in 392 The Cancer Genome Atlas UBC cases., Results: In the pedigrees, segregating variants were family-specific although the identified genes clustered in common pathways, most notably DNA repair ( MLH1 and MSH2 ) and cellular metabolism ( IDH1 and ME1 ). In the familial UBC group, the proportion of pathogenic and likely pathogenic variants was significantly higher in cases compared with controls ( P = .003). Pathogenic and likely pathogenic variant load was also significantly increased in genes involved in cilia biogenesis ( P = .001). In addition, a pathogenic variant in CHEK2 (NM&#95;007194.4:c.1100del; p.T367Mfs*15) was over-represented in cases (variant frequency = 2.6%; 95% CI, 0.71 to 6.52) compared with controls (variant frequency = 0.21%; 95% CI, 0.01 to 1.15), but was not statistically significant., Conclusion: These results point to a complex polygenic predisposition to UBC. Despite heterogeneity, the genes cluster in several biologically relevant pathways and processes, for example, DNA repair, cilia biogenesis, and cellular metabolism. Larger studies are required to determine the importance of CHEK2 in UBC etiology., Competing Interests: Talia Wegman-Ostrosky Travel, Accommodations, Expenses: Pfizer Bin Zhu Employment: GeneDx BioReference Douglas R. Stewart Employment: Genome Medical No other potential conflicts of interest were reported. Talia Wegman-Ostrosky Travel, Accommodations, Expenses: Pfizer Bin Zhu Employment: GeneDx BioReference Douglas R. Stewart Employment: Genome Medical No other potential conflicts of interest were reported., (© 2021 by American Society of Clinical Oncology.)

Introduction: The optimal length for clinical follow-up of renal cell carcinoma (RCC) patients is unclear. We evaluated the impact of ISUP/WHO tumor grade and histological subtype on short- and long-term survival and risk of recurrence/metastasis in a large cohort of RCC patients., Patients and Methods: We studied 1679 RCC patients from a single referral center in Italy. Adjusted hazard ratios for overall survival were estimated using Cox regression models. Adjusted absolute risk of developing recurrence or metastasis was computed considering competing risks of mortality., Results: During up to 13 years of follow-up, 175 (10.4%) RCC patients died, of whom 92% beyond 5 years. Hazard ratio of grade IV clear cell carcinomas (ccRCC) was 3.82 compared to grade II. Notably, 33% of recurrences and 56% of distant metastases occurred beyond 5 years of follow-up. The estimated probabilities of recurrence/metastasis were 15% and 45% within and beyond 5 years of follow-up, respectively. After 5 years, the absolute risk of recurrences increased also for papillary renal cell carcinoma type I (35.2%) and grade I ccRCC (17%)., Conclusion: After 5 years of follow-up, both risk of mortality and recurrences or metastases were high and were modified by histological types and tumor grade. These data strongly support histology- and grade-tailored surveillance strategies and long-term follow-up for RCC patients., (Copyright © 2021. Published by Elsevier Inc.)

Genome-wide association studies (GWASs) have identified a melanoma-associated locus on chromosome band 7p21.1 with rs117132860 as the lead SNP and a secondary independent signal marked by rs73069846. rs117132860 is also associated with tanning ability and cutaneous squamous cell carcinoma (cSCC). Because ultraviolet radiation (UVR) is a key environmental exposure for all three traits, we investigated the mechanisms by which this locus contributes to melanoma risk, focusing on cellular response to UVR. Fine-mapping of melanoma GWASs identified four independent sets of candidate causal variants. A GWAS region-focused Capture-C study of primary melanocytes identified physical interactions between two causal sets and the promoter of the aryl hydrocarbon receptor (AHR). Subsequent chromatin state annotation, eQTL, and luciferase assays identified rs117132860 as a functional variant and reinforced AHR as a likely causal gene. Because AHR plays critical roles in cellular response to dioxin and UVR, we explored links between this SNP and AHR expression after both 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) and ultraviolet B (UVB) exposure. Allele-specific AHR binding to rs117132860-G was enhanced following both, consistent with predicted weakened AHR binding to the risk/poor-tanning rs117132860-A allele, and allele-preferential AHR expression driven from the protective rs117132860-G allele was observed following UVB exposure. Small deletions surrounding rs117132860 introduced via CRISPR abrogates AHR binding, reduces melanocyte cell growth, and prolongs growth arrest following UVB exposure. These data suggest AHR is a melanoma susceptibility gene at the 7p21.1 risk locus and rs117132860 is a functional variant within a UVB-responsive element, leading to allelic AHR expression and altering melanocyte growth phenotypes upon exposure., Competing Interests: Declaration of interests The authors declare no competing interests., (Published by Elsevier Inc.)