Genome-wide association studies have identified hundreds of genetic variants associated with specific cancers. A few of these risk regions have been associated with more than one cancer site; however, a systematic evaluation of the associations between risk variants for other cancers and lung cancer risk has yet to be performed.

Lung cancer etiology is multifactorial, and growing evidence has indicated that long non-coding RNAs (lncRNAs) are important players in lung carcinogenesis. We performed a large-scale meta-analysis of 690,564 SNPs in 15,531 autosomal lncRNAs by using datasets from six previously published genome-wide association studies (GWASs) from the Transdisciplinary Research in Cancer of the Lung (TRICL) consortium in populations of European ancestry. Previously unreported significant SNPs (P value < 1 × 10-7) were further validated in two additional independent lung cancer GWAS datasets from Harvard University and deCODE. In the final meta-analysis of all eight GWAS datasets with 17,153 cases and 239,337 controls, a novel risk SNP rs114020893 in the lncRNA NEXN-AS1 region at 1p31.1 remained statistically significant (odds ratio = 1.17; 95% confidence interval = 1.11-1.24; P = 8.31 × 10-9). In further in silico analysis, rs114020893 was predicted to change the secondary structure of the lncRNA. Our finding indicates that SNP rs114020893 of NEXN-AS1 at 1p31.1 may contribute to lung cancer susceptibility.

Body mass index (BMI) is inversely associated with lung cancer risk in observational studies, even though it increases the risk of several other cancers, which could indicate confounding by tobacco smoking or reverse causality. We used the two-sample Mendelian randomization (MR) approach to circumvent these limitations of observational epidemiology by constructing a genetic instrument for BMI, based on results from the GIANT consortium, which was evaluated in relation to lung cancer risk using GWAS results on 16,572 lung cancer cases and 21,480 controls. Results were stratified by histological subtype, smoking status and sex. An increase of one standard deviation (SD) in BMI (4.65 Kg/m(2)) raised the risk for lung cancer overall (OR = 1.13; P = 0.10). This was driven by associations with squamous cell (SQ) carcinoma (OR = 1.45; P = 1.2 × 10(-3)) and small cell (SC) carcinoma (OR = 1.81; P = 0.01). An inverse trend was seen for adenocarcinoma (AD) (OR = 0.82; P = 0.06). In stratified analyses, a 1 SD increase in BMI was inversely associated with overall lung cancer in never smokers (OR = 0.50; P = 0.02). These results indicate that higher BMI may increase the risk of certain types of lung cancer, in particular SQ and SC carcinoma.

Recent genome-wide association studies (GWASs) have identified common genetic variants at 5p15.33, 6p21-6p22 and 15q25.1 associated with lung cancer risk. Several other genetic regions including variants of CHEK2 (22q12), TP53BP1 (15q15) and RAD52 (12p13) have been demonstrated to influence lung cancer risk in candidate- or pathway-based analyses. To identify novel risk variants for lung cancer, we performed a meta-analysis of 16 GWASs, totaling 14 900 cases and 29 485 controls of European descent. Our data provided increased support for previously identified risk loci at 5p15 (P = 7.2 × 10(-16)), 6p21 (P = 2.3 × 10(-14)) and 15q25 (P = 2.2 × 10(-63)). Furthermore, we demonstrated histology-specific effects for 5p15, 6p21 and 12p13 loci but not for the 15q25 region. Subgroup analysis also identified a novel disease locus for squamous cell carcinoma at 9p21 (CDKN2A/p16(INK4A)/p14(ARF)/CDKN2B/p15(INK4B)/ANRIL; rs1333040, P = 3.0 × 10(-7)) which was replicated in a series of 5415 Han Chinese (P = 0.03; combined analysis, P = 2.3 × 10(-8)). This large analysis provides additional evidence for the role of inherited genetic susceptibility to lung cancer and insight into biological differences in the development of the different histological types of lung cancer.

Genome-wide association studies (GWAS) identified the chromosome 15q25.1 locus as a leading susceptibility region for lung cancer. However, the pathogenic pathways, through which susceptibility SNPs within chromosome 15q25.1 affects lung cancer risk, have not been explored. We analyzed three cohorts with GWAS data consisting 42,901 individuals and lung expression quantitative trait loci (eQTL) data on 409 individuals to identify and validate the underlying pathways and to investigate the combined effect of genes from the identified susceptibility pathways. The KEGG neuroactive ligand receptor interaction pathway, two Reactome pathways, and 22 Gene Ontology terms were identified and replicated to be significantly associated with lung cancer risk, with P values less than 0.05 and FDR less than 0.1. Functional annotation of eQTL analysis results showed that the neuroactive ligand receptor interaction pathway and gated channel activity were involved in lung cancer risk. These pathways provide important insights for the etiology of lung cancer.

Germline copy number variants (CNVs) increase risk for many diseases, yet detection of CNVs and quantifying their contribution to disease risk in large-scale studies is challenging due to biological and technical sources of heterogeneity that vary across the genome within and between samples.

Impaired lung function is often caused by cigarette smoking, making it challenging to disentangle its role in lung cancer susceptibility. Investigation of the shared genetic basis of these phenotypes in the UK Biobank and International Lung Cancer Consortium (29,266 cases, 56,450 controls) shows that lung cancer is genetically correlated with reduced forced expiratory volume in one second (FEV1: rg = 0.098, p = 2.3 × 10-8) and the ratio of FEV1 to forced vital capacity (FEV1/FVC: rg = 0.137, p = 2.0 × 10-12). Mendelian randomization analyses demonstrate that reduced FEV1 increases squamous cell carcinoma risk (odds ratio (OR) = 1.51, 95% confidence intervals: 1.21-1.88), while reduced FEV1/FVC increases the risk of adenocarcinoma (OR = 1.17, 1.01-1.35) and lung cancer in never smokers (OR = 1.56, 1.05-2.30). These findings support a causal role of pulmonary impairment in lung cancer etiology. Integrative analyses reveal that pulmonary function instruments, including 73 novel variants, influence lung tissue gene expression and implicate immune-related pathways in mediating the observed effects on lung carcinogenesis.

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.

Germline variation and smoking are independently associated with pancreatic ductal adenocarcinoma (PDAC). We conducted genome-wide smoking interaction analysis of PDAC using genotype data from four previous genome-wide association studies in individuals of European ancestry (7,937 cases and 11,774 controls). Examination of expression quantitative trait loci data from the Genotype-Tissue Expression Project followed by colocalization analysis was conducted to determine whether there was support for common SNP(s) underlying the observed associations. Statistical tests were two sided and P < 5 × 10-8 was considered statistically significant. Genome-wide significant evidence of qualitative interaction was identified on chr2q21.3 in intron 5 of the transmembrane protein 163 (TMEM163) and upstream of the cyclin T2 (CCNT2). The most significant SNP using the Empirical Bayes method, in this region that included 45 significantly associated SNPs, was rs1818613 [per allele OR in never smokers 0.87, 95% confidence interval (CI), 0.82-0.93; former smokers 1.00, 95% CI, 0.91-1.07; current smokers 1.25, 95% CI 1.12-1.40, P interaction = 3.08 × 10-9). Examination of the Genotype-Tissue Expression Project data demonstrated an expression quantitative trait locus in this region for TMEM163 and CCNT2 in several tissue types. Colocalization analysis supported a shared SNP, rs842357, in high linkage disequilibrium with rs1818613 (r 2 = 0. 94) driving both the observed interaction and the expression quantitative trait loci signals. Future studies are needed to confirm and understand the differential biologic mechanisms by smoking status that contribute to our PDAC findings. SIGNIFICANCE: This large genome-wide interaction study identifies a susceptibility locus on 2q21.3 that significantly modified PDAC risk by smoking status, providing insight into smoking-associated PDAC, with implications for prevention.

Clonal hematopoiesis (CH) is a phenomenon caused by expansion of white blood cells descended from a single hematopoietic stem cell. While CH can be associated with leukemia and some solid tumors, the relationship between CH and lung cancer remains largely unknown. To help clarify this relationship, we analyzed whole-exome sequencing (WES) data from 1,958 lung cancer cases and controls. Potential CH mutations were identified by a set of hierarchical filtering criteria in different exonic regions, and the associations between the number of CH mutations and clinical traits were investigated. Family history of lung cancer (FHLC) may exert diverse influences on the accumulation of CH mutations in different age groups. In younger subjects, FHLC was the strongest risk factor for CH mutations. Association analysis of genome-wide genetic variants identified dozens of genetic loci associated with CH mutations, including a candidate SNP rs2298110, which may promote CH by increasing expression of a potential leukemia promoter gene OTUD3. Hundreds of potentially novel CH mutations were identified, and smoking was found to potentially shape the CH mutational signature. Genetic variants and lung cancer risk factors, especially FHLC, correlated with CH. These analyses improve our understanding of the relationship between lung cancer and CH, and future experimental studies will be necessary to corroborate the uncovered correlations. SIGNIFICANCE: Analysis of whole-exome sequencing data uncovers correlations between clonal hematopoiesis and lung cancer risk factors, identifies genetic variants correlated with clonal hematopoiesis, and highlights hundreds of potential novel clonal hematopoiesis mutations.

Deucravacitinib, a novel, selective inhibitor of TYK2 is currently under review at the FDA and EMA for treatment of moderate-to-severe plaque psoriasis. It is unclear whether recent safety concerns (ie, elevated rates of lung cancer and lymphoma) related to similar medications (ie, other JAK inhibitors) are shared with this novel TYK2 inhibitor. We used a partial loss-of-function variant in TYK2 (rs34536443), previously shown to protect against psoriasis and other autoimmune diseases, to evaluate the potential effect of therapeutic TYK2 inhibition on risk of lung cancer and non-Hodgkin lymphoma. Summary genetic association data on lung cancer risk were obtained from a GWAS meta-analysis of 29 266 cases and 56 450 controls in the Integrative Analysis of Lung Cancer Risk and Aetiology (INTEGRAL) consortium. Summary genetic association data on non-Hodgkin lymphoma risk were obtained from a GWAS meta-analysis of 8489 cases and 374 506 controls in the UK Biobank and InterLymph consortium. In the primary analysis, each copy of the minor allele of rs34536443, representing partial TYK2 inhibition, was associated with an increased risk of lung cancer (OR 1.15, 95% CI 1.09-1.23, P = 2.29 × 10-6 ) and non-Hodgkin lymphoma (OR 1.18, 95% CI 1.05-1.33, P = 5.25 × 10-3 ). Our analyses using an established partial loss-of-function mutation to mimic TYK2 inhibition provide genetic evidence that therapeutic TYK2 inhibition may increase risk of lung cancer and non-Hodgkin lymphoma. These findings, consistent with recent reports from postmarketing trials of similar JAK inhibitors, could have important implications for future safety assessment of deucravacitinib and other TYK2 inhibitors in development.

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).

Tumour-promoting inflammation is a "hallmark" of cancer and conventional epidemiological studies have reported links between various inflammatory markers and cancer risk. The causal nature of these relationships and, thus, the suitability of these markers as intervention targets for cancer prevention is unclear.

We conducted a genome-wide association study of oral cavity and pharyngeal cancer in 6,034 cases and 6,585 controls from Europe, North America and South America. We detected eight significantly associated loci (P < 5 × 10-8), seven of which are new for these cancer sites. Oral and pharyngeal cancers combined were associated with loci at 6p21.32 (rs3828805, HLA-DQB1), 10q26.13 (rs201982221, LHPP) and 11p15.4 (rs1453414, OR52N2-TRIM5). Oral cancer was associated with two new regions, 2p23.3 (rs6547741, GPN1) and 9q34.12 (rs928674, LAMC3), and with known cancer-related loci-9p21.3 (rs8181047, CDKN2B-AS1) and 5p15.33 (rs10462706, CLPTM1L). Oropharyngeal cancer associations were limited to the human leukocyte antigen (HLA) region, and classical HLA allele imputation showed a protective association with the class II haplotype HLA-DRB1*1301-HLA-DQA1*0103-HLA-DQB1*0603 (odds ratio (OR) = 0.59, P = 2.7 × 10-9). Stratified analyses on a subgroup of oropharyngeal cases with information available on human papillomavirus (HPV) status indicated that this association was considerably stronger in HPV-positive (OR = 0.23, P = 1.6 × 10-6) than in HPV-negative (OR = 0.75, P = 0.16) cancers.

Registry-based epidemiologic studies suggest associations between chronic inflammatory intestinal diseases and pancreatic ductal adenocarcinoma (PDAC). As genetic susceptibility contributes to a large proportion of chronic inflammatory intestinal diseases, we hypothesize that the genomic regions surrounding established genome-wide associated variants for these chronic inflammatory diseases are associated with PDAC. We examined the association between PDAC and genomic regions (±500 kb) surrounding established common susceptibility variants for ulcerative colitis, Crohn's disease, inflammatory bowel disease, celiac disease, chronic pancreatitis, and primary sclerosing cholangitis. We analyzed summary statistics from genome-wide association studies data for 8,384 cases and 11,955 controls of European descent from two large consortium studies using the summary data-based adaptive rank truncated product method to examine the overall association of combined genomic regions for each inflammatory disease group. Combined genomic susceptibility regions for ulcerative colitis, Crohn disease, inflammatory bowel disease, and chronic pancreatitis were associated with PDAC at P values < 0.05 (0.0040, 0.0057, 0.011, and 3.4 × 10-6, respectively). After excluding the 20 PDAC susceptibility regions (±500 kb) previously identified by GWAS, the genomic regions for ulcerative colitis, Crohn disease, and inflammatory bowel disease remained associated with PDAC (P = 0.0029, 0.0057, and 0.0098, respectively). Genomic regions for celiac disease (P = 0.22) and primary sclerosing cholangitis (P = 0.078) were not associated with PDAC. Our results support the hypothesis that genomic regions surrounding variants associated with inflammatory intestinal diseases, particularly, ulcerative colitis, Crohn disease, inflammatory bowel disease, and chronic pancreatitis are associated with PDAC. SIGNIFICANCE: The joint effects of common variants in genomic regions containing susceptibility loci for inflammatory bowel disease and chronic pancreatitis are associated with PDAC and may provide insights to understanding pancreatic cancer etiology.

Recent technological advancements have permitted high-throughput measurement of the human genome, epigenome, metabolome, transcriptome, and proteome at the population level. We hypothesized that subsets of genes identified from omic studies might have closely related biological functions and thus might interact directly at the network level. Therefore, we conducted an integrative analysis of multi-omic datasets of non-small cell lung cancer (NSCLC) to search for association patterns beyond the genome and transcriptome. A large, complex, and robust gene network containing well-known lung cancer-related genes, including EGFR and TERT, was identified from combined gene lists for lung adenocarcinoma. Members of the hypoxia-inducible factor (HIF) gene family were at the center of this network. Subsequent sequencing of network hub genes within a subset of samples from the Transdisciplinary Research in Cancer of the Lung-International Lung Cancer Consortium (TRICL-ILCCO) consortium revealed a SNP (rs12614710) in EPAS1 associated with NSCLC that reached genome-wide significance (OR = 1.50; 95% CI: 1.31-1.72; p = 7.75 × 10-9). Using imputed data, we found that this SNP remained significant in the entire TRICL-ILCCO consortium (p = .03). Additional functional studies are warranted to better understand interrelationships among genetic polymorphisms, DNA methylation status, and EPAS1 expression.

This study investigated nicotine dependence as an independent risk factor for upper aerodigestive tract (UADT) cancers, including lung and head and neck cancers (HNC). The study aimed to isolate the direct effect of nicotine dependence, independent of tobacco smoking.

Germline genetic variation contributes to lung cancer (LC) susceptibility. Previous genome-wide association studies (GWAS) have implicated susceptibility loci involved in smoking behaviors and DNA repair genes, but further work is required to identify susceptibility variants.

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 (Pinteraction = 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.

Squamous cell carcinomas (SqCC) of the aerodigestive tract have similar etiological risk factors. Although genetic risk variants for individual cancers have been identified, an agnostic, genome-wide search for shared genetic susceptibility has not been performed. To identify novel and pleotropic SqCC risk variants, we performed a meta-analysis of GWAS data on lung SqCC (LuSqCC), oro/pharyngeal SqCC (OSqCC), laryngeal SqCC (LaSqCC) and esophageal SqCC (ESqCC) cancers, totaling 13,887 cases and 61,961 controls of European ancestry. We identified one novel genome-wide significant (Pmeta<5x10-8) aerodigestive SqCC susceptibility loci in the 2q33.1 region (rs56321285, TMEM273). Additionally, three previously unknown loci reached suggestive significance (Pmeta<5x10-7): 1q32.1 (rs12133735, near MDM4), 5q31.2 (rs13181561, TMEM173) and 19p13.11 (rs61494113, ABHD8). Multiple previously identified loci for aerodigestive SqCC also showed evidence of pleiotropy in at least another SqCC site, these include: 4q23 (ADH1B), 6p21.33 (STK19), 6p21.32 (HLA-DQB1), 9p21.33 (CDKN2B-AS1) and 13q13.1(BRCA2). Gene-based association and gene set enrichment identified a set of 48 SqCC-related genes rel to DNA damage and epigenetic regulation pathways. Our study highlights the importance of cross-cancer analyses to identify pleiotropic risk loci of histology-related cancers arising at distinct anatomical sites.