Large-scale whole exome sequencing studies identify two genes,CTSL and APOE, associated with lung cancer.

Common genetic variants associated with lung cancer have been well studied in the past decade. However, only 12.3% heritability has been explained by these variants. In this study, we investigate the contribution of rare variants (RVs) (minor allele frequency <0.01) to lung cancer through two large whole exome sequencing case-control studies. We first performed gene-based association tests using a novel Bayes Factor statistic in the International Lung Cancer Consortium, the discovery study (European, 1042 cases vs. 881 controls). The top genes identified are further assessed in the UK Biobank (European, 630 cases vs. 172 864 controls), the replication study. After controlling for the false discovery rate, we found two genes, CTSL and APOE, significantly associated with lung cancer in both studies. Single variant tests in UK Biobank identified 4 RVs (3 missense variants) in CTSL and 2 RVs (1 missense variant) in APOE stongly associated with lung cancer (OR between 2.0 and 139.0). The role of these genetic variants in the regulation of CTSL or APOE expression remains unclear. If such a role is established, this could have important therapeutic implications for lung cancer patients.

Radiomics analysis to predict pulmonary nodule malignancy using machine learning approaches.

BACKGROUND: Low-dose CT screening can reduce lung cancer-related mortality. However, most screen-detected pulmonary abnormalities do not develop into cancer and it often remains challenging to identify malignant nodules, particularly among indeterminate nodules. We aimed to develop and assess prediction models based on radiological features to discriminate between benign and malignant pulmonary lesions detected on a baseline screen. METHODS: Using four international lung cancer screening studies, we extracted 2060 radiomic features for each of 16 797 nodules (513 malignant) among 6865 participants. After filtering out low-quality radiomic features, 642 radiomic and 9 epidemiological features remained for model development. We used cross-validation and grid search to assess three machine learning (ML) models (eXtreme Gradient Boosted Trees, random forest, least absolute shrinkage and selection operator (LASSO)) for their ability to accurately predict risk of malignancy for pulmonary nodules. We report model performance based on the area under the curve (AUC) and calibration metrics in the held-out test set. RESULTS: The LASSO model yielded the best predictive performance in cross-validation and was fit in the full training set based on optimised hyperparameters. Our radiomics model had a test-set AUC of 0.93 (95% CI 0.90 to 0.96) and outperformed the established Pan-Canadian Early Detection of Lung Cancer model (AUC 0.87, 95% CI 0.85 to 0.89) for nodule assessment. Our model performed well among both solid (AUC 0.93, 95% CI 0.89 to 0.97) and subsolid nodules (AUC 0.91, 95% CI 0.85 to 0.95). CONCLUSIONS: We developed highly accurate ML models based on radiomic and epidemiological features from four international lung cancer screening studies that may be suitable for assessing indeterminate screen-detected pulmonary nodules for risk of malignancy.

A risk prediction model for head and neck cancers incorporating lifestyle factors, HPV serology and genetic markers.

Head and neck cancer is often diagnosed late and prognosis for most head and neck cancer patients remains poor. To aid early detection, we developed a risk prediction model based on demographic and lifestyle risk factors, human papillomavirus (HPV) serological markers and genetic markers. A total of 10 126 head and neck cancer cases and 5254 controls from five North American and European studies were included. HPV serostatus was determined by antibodies for HPV16 early oncoproteins (E6, E7) and regulatory early proteins (E1, E2, E4). The data were split into a training set (70%) for model development and a hold-out testing set (30%) for model performance evaluation, including discriminative ability and calibration. The risk models including demographic, lifestyle risk factors and polygenic risk score showed a reasonable predictive accuracy for head and neck cancer overall. A risk model that also included HPV serology showed substantially improved predictive accuracy for oropharyngeal cancer (AUC = 0.94, 95% CI = 0.92-0.95 in men and AUC = 0.92, 95% CI = 0.88-0.95 in women). The 5-year absolute risk estimates showed distinct trajectories by risk factor profiles. Based on the UK Biobank cohort, the risks of developing oropharyngeal cancer among 60 years old and HPV16 seropositive in the next 5 years ranged from 5.8% to 14.9% with an average of 8.1% for men, 1.3% to 4.4% with an average of 2.2% for women. Absolute risk was generally higher among individuals with heavy smoking, heavy drinking, HPV seropositivity and those with higher polygenic risk score. These risk models may be helpful for identifying people at high risk of developing head and neck cancer.

Systematic analyses of regulatory variants in DNase I hypersensitive sites identified two novel lung cancer susceptibility loci.

DNase I hypersensitive sites (DHS) are abundant in regulatory elements, such as promoter, enhancer and transcription factor binding sites. Many studies have revealed that disease-associated variants were concentrated in DHS-related regions. However, limited studies are available on the roles of DHS-related variants in lung cancer. In this study, we performed a large-scale case-control study with 20 871 lung cancer cases and 15 971 controls to evaluate the associations between regulatory genetic variants in DHS and lung cancer susceptibility. The expression quantitative trait loci (eQTL) analysis and pathway-enrichment analysis were performed to identify the possible target genes and pathways. In addition, we performed motif-based analysis to explore the lung-cancer-related motifs using sequence kernel association test. Two novel variants, rs186332 in 20q13.3 (C>T, odds ratio [OR] = 1.17, 95% confidence interval [95% CI]: 1.10-1.24, P = 8.45 × 10-7) and rs4839323 in 1p13.2 (T>C, OR = 0.92, 95% CI: 0.89-0.95, P = 1.02 × 10-6) showed significant association with lung cancer risk. The eQTL analysis suggested that these two SNPs might regulate the expression of MRGBP and SLC16A1, respectively. What's more, the expression of both MRGBP and SLC16A1 was aberrantly elevated in lung tumor tissues. The motif-based analysis identified 10 motifs related to the risk of lung cancer (P < 1.71 × 10-4). Our findings suggested that variants in DHS might modify lung cancer susceptibility through regulating the expression of surrounding genes. This study provided us a deeper insight into the roles of DHS-related genetic variants for lung cancer.

Genome-wide association study of familial lung cancer.

To identify genetic variation associated with lung cancer risk, we performed a genome-wide association analysis of 685 lung cancer cases that had a family history of two or more first or second degree relatives compared with 744 controls without lung cancer that were genotyped on an Illumina Human OmniExpressExome-8v1 array. To ensure robust results, we further evaluated these findings using data from six additional studies that were assembled through the Transdisciplinary Research on Cancer of the Lung Consortium comprising 1993 familial cases and 33 690 controls. We performed a meta-analysis after imputation of all variants using the 1000 Genomes Project Phase 1 (version 3 release date September 2013). Analyses were conducted for 9 327 222 SNPs integrating data from the two sources. A novel variant on chromosome 4p15.31 near the LCORL gene and an imputed rare variant intergenic between CDKN2A and IFNA8 on chromosome 9p21.3 were identified at a genome-wide level of significance for squamous cell carcinomas. Additionally, associations of CHRNA3 and CHRNA5 on chromosome 15q25.1 in sporadic lung cancer were confirmed at a genome-wide level of significance in familial lung cancer. Previously identified variants in or near CHRNA2, BRCA2, CYP2A6 for overall lung cancer, TERT, SECISPB2L and RTEL1 for adenocarcinoma and RAD52 and MHC for squamous carcinoma were significantly associated with lung cancer.

Two BRM promoter polymorphisms predict poor survival in patients with hepatocellular carcinoma.

Polymorphisms in the promoter of the BRM gene, a critical subunit of the chromatin remodeling SWI/SNF complex, have previously been implicated in risk and prognosis in Caucasian-predominant lung, head and neck, esophageal, and pancreatic cancers, and in hepatocellular cancers in Asians. We investigated the role of these polymorphisms in hepatocellular carcinoma (HCC) risk and prognosis. HCC cases were recruited in a comprehensive cancer center while the matched controls were recruited from family practice units from the same catchment area. For risk analyses, unconditional logistic regression analyses were performed in HCC patients and matched healthy controls. Overall survival analyses were performed using Cox proportional hazard models, Kaplan-Meier curves, and log-rank tests. In 266 HCC cases and 536 controls, no association between either BRM promoter polymorphism (BRM-741 or BRM-1321) and risk of HCC was identified (P > 0.10 for all comparisons). There was significant worsening of overall survival as the number of variant alleles increased: BRM-741 per variant allele adjusted hazards ratio (aHR) 5.77, 95% confidence interval (CI) 2.89-11.54 and BRM-1321 per variant allele aHR 4.09, 95%CI 2.22-7.51. The effects of these two polymorphisms were at least additive, where individuals who were double homozygotes for the variant alleles had a 45-fold increase in risk of death when compared to those who were double wild-type for the two polymorphisms. Two BRM promoter polymorphisms were strongly associated with HCC prognosis but were not associated with increased HCC susceptibility. The association was strongest in double homozygotes for the allele variants.

Novel genetic variants in the P38MAPK pathway gene ZAK and susceptibility to lung cancer.

The P38MAPK pathway participates in regulating cell cycle, inflammation, development, cell death, cell differentiation, and tumorigenesis. Genetic variants of some genes in the P38MAPK pathway are reportedly associated with lung cancer risk. To substantiate this finding, we used six genome-wide association studies (GWASs) to comprehensively investigate the associations of 14 904 single nucleotide polymorphisms (SNPs) in 108 genes of this pathway with lung cancer risk. We identified six significant lung cancer risk-associated SNPs in two genes (CSNK2B and ZAK) after correction for multiple comparisons by a false discovery rate (FDR) <0.20. After removal of three CSNK2B SNPs that are located in the same locus previously reported by GWAS, we performed the LD analysis and found that rs3769201 and rs7604288 were in high LD. We then chose two independent representative SNPs of rs3769201 and rs722864 in ZAK for further analysis. We also expanded the analysis by including these two SNPs from additional GWAS datasets of Harvard University (984 cases and 970 controls) and deCODE (1319 cases and 26 380 controls). The overall effects of these two SNPs were assessed using all eight GWAS datasets (OR = 0.92, 95%CI = 0.89-0.95, and P = 1.03 × 10-5 for rs3769201; OR = 0.91, 95%CI = 0.88-0.95, and P = 2.03 × 10-6 for rs722864). Finally, we performed an expression quantitative trait loci (eQTL) analysis and found that these two SNPs were significantly associated with ZAK mRNA expression levels in lymphoblastoid cell lines. In conclusion, the ZAK rs3769201 and rs722864 may be functional susceptibility loci for lung cancer risk.

Functional variants in DCAF4 associated with lung cancer risk in European populations.

Cullin-RING ubiquitin ligases (CRLs) responsible for substrate specificity of ubiquitination play a key role in cell-cycle control and DNA damage response. In this study, we assessed associations between 16 599 SNPs in 115 CRL genes and lung cancer risk by using summary data of six published genome-wide association studies (GWASs) of 12 160 cases and 16 838 cases of European ancestry. As a result, we identified three independent SNPs in DCAF4 (rs117781739, rs12587742 and rs2240980) associated with lung cancer risk (odds ratio = 0.91, 1.09 and 1.09, respectively; 95% confidence interval = 0.88-0.95, 1.05-1.14 and 1.05-1.13, respectively; and P = 3.99 × 10-6, 4.97 × 10-5 and 1.44 × 10-5, respectively) after multiple comparison correction by a false discovery rate <0.05. Since SNP rs12587742 is located within the promoter region and one CpG island of DCAF4, we further performed in silico functional analyses and found that the rs12587742 variant A allele was associated with an increased mRNA expression (P = 2.20 × 10-16, 1.79 × 10-13 and 0.001 in blood cells, normal lung tissues and tumor tissues of lung squamous carcinoma, respectively) and a decreased methylation status (P = 2.48 × 10-9 and 0.032 in adipose and lung tumor tissues, respectively). Moreover, evidence from differential expression analyses further supported an oncogenic effect of DCAF4 on lung cancer, with higher mRNA levels in both lung squamous carcinoma and adenocarcinoma (P = 4.48 × 10-11 and 1.22 × 10-9, respectively) than in adjacent normal tissues. Taken together, our results suggest that rs12587742 is associated with an increased lung cancer risk, possibly by up-regulating mRNA expression and decreasing methylation status of DCAF4.

Identification of shared and unique susceptibility pathways among cancers of the lung, breast, and prostate from genome-wide association studies and tissue-specific protein interactions.

Results from genome-wide association studies (GWAS) have indicated that strong single-gene effects are the exception, not the rule, for most diseases. We assessed the joint effects of germline genetic variations through a pathway-based approach that considers the tissue-specific contexts of GWAS findings. From GWAS meta-analyses of lung cancer (12 160 cases/16 838 controls), breast cancer (15 748 cases/18 084 controls) and prostate cancer (14 160 cases/12 724 controls) in individuals of European ancestry, we determined the tissue-specific interaction networks of proteins expressed from genes that are likely to be affected by disease-associated variants. Reactome pathways exhibiting enrichment of proteins from each network were compared across the cancers. Our results show that pathways associated with all three cancers tend to be broad cellular processes required for growth and survival. Significant examples include the nerve growth factor (P = 7.86 × 10(-33)), epidermal growth factor (P = 1.18 × 10(-31)) and fibroblast growth factor (P = 2.47 × 10(-31)) signaling pathways. However, within these shared pathways, the genes that influence risk largely differ by cancer. Pathways found to be unique for a single cancer focus on more specific cellular functions, such as interleukin signaling in lung cancer (P = 1.69 × 10(-15)), apoptosis initiation by Bad in breast cancer (P = 3.14 × 10(-9)) and cellular responses to hypoxia in prostate cancer (P = 2.14 × 10(-9)). We present the largest comparative cross-cancer pathway analysis of GWAS to date. Our approach can also be applied to the study of inherited mechanisms underlying risk across multiple diseases in general.

Pathway-analysis of published genome-wide association studies of lung cancer: A potential role for the CYP4F3 locus.

The fatty acids (FAs) metabolism is suggested to play a pivotal role in the development of lung cancer, and we explored that by conducting a pathway-based analysis. We performed a meta-analysis of published datasets of six genome wide association studies (GWASs) from the Transdisciplinary Research in Cancer of the Lung (TRICL) consortium, which included 12 160 cases with lung cancer and 16 838 cancer-free controls. A total of 30 722 single-nucleotide polymorphisms (SNPs) from 317 genes relevant to FA metabolic pathways were identified. An additional dataset from the Harvard Lung Cancer Study with 984 cases and 970 healthy controls was also added to the final meta-analysis. In the initial meta-analysis, 26 of 28 SNPs that passed false discovery rate multiple tests were mapped to the CYP4F3 gene. Among the 26 top ranked hits was a proxy SNP, CYP4F3 rs4646904 (P = 8.65 × 10-6 , FDR = 0.018), which is suggested to change splicing pattern/efficiency and to be associated with gene expression levels. However, after adding data of rs4646904 from the Harvard GWAS, the significance in the combined analysis was reduced to P = 3.52 × 10-3 [odds ratio (OR) = 1.07, 95% confidence interval (95%CI) = 1.03-1.12]. Interestingly, the small Harvard dataset also pointed to the same direction of the association in subgroups of smokers (OR = 1.07) and contributed to a combined OR of 1.13 (95% CI = 1.06-1.20, P = 6.70 × 10-5 ). The results suggest that a potentially functional SNP in CYP4F3 (rs4646904) may contribute to the etiology of lung cancer, especially in smokers. Additional mechanistic studies are warranted to unravel the potential biological significance of the finding.

Susceptibility loci of CNOT6 in the general mRNA degradation pathway and lung cancer risk-A re-analysis of eight GWASs.

PURPOSE: mRNA degradation is an important regulatory step for controlling gene expression and cell functions. Genetic abnormalities involved in mRNA degradation genes were found to be associated with cancer risks. Therefore, we systematically investigated the roles of genetic variants in the general mRNA degradation pathway in lung cancer risk. EXPERIMENTAL DESIGN: Meta-analyses were conducted using summary data from six lung cancer genome-wide association studies (GWASs) from the Transdisciplinary Research in Cancer of the Lung and additional two GWASs from Harvard University and deCODE in the International Lung Cancer Consortium. Expression quantitative trait loci analysis (eQTL) was used for in silico functional validation of the identified significant susceptibility loci. RESULTS: This pathway-based analysis included 6816 single nucleotide polymorphisms (SNP) in 68 genes in 14 463 lung cancer cases and 44 188 controls. In the single-locus analysis, we found that 20 SNPs were associated with lung cancer risk with a false discovery rate threshold of <0.05. Among the 11 newly identified SNPs in CNOT6, which were in high linkage disequilibrium, the rs2453176 with a RegulomDB score "1f" was chosen as the tagSNP for further analysis. We found that the rs2453176 T allele was significantly associated with lung cancer risk (odds ratio = 1.11, 95% confidence interval = 1.04-1.18) in the eight GWASs. In the eQTL analysis, we found that levels of CNOT6 mRNA expression were significantly correlated with the rs2453176 T allele, which provided additional biological basis for the observed positive association. CONCLUSION: The CNOT6 rs2453176 SNP may be a new functional susceptible locus for lung cancer risk. © 2016 Wiley Periodicals, Inc.

Polymorphisms of the centrosomal gene (FGFR1OP) and lung cancer risk: a meta-analysis of 14,463 cases and 44,188 controls.

Centrosome abnormalities are often observed in premalignant lesions and in situ tumors and have been associated with aneuploidy and tumor development. We investigated the associations of 9354 single-nucleotide polymorphisms (SNPs) in 106 centrosomal genes with lung cancer risk by first using the summary data from six published genome-wide association studies (GWASs) of the Transdisciplinary Research in Cancer of the Lung (TRICL) (12,160 cases and 16 838 controls) and then conducted in silico replication in two additional independent lung cancer GWASs of Harvard University (984 cases and 970 controls) and deCODE (1319 cases and 26,380 controls). A total of 44 significant SNPs with false discovery rate (FDR) ≤ 0.05 were mapped to one novel gene FGFR1OP and two previously reported genes (TUBB and BRCA2). After combined the results from TRICL with those from Harvard and deCODE, the most significant association (P combined = 8.032 × 10(-6)) was with rs151606 within FGFR1OP. The rs151606 T>G was associated with an increased risk of lung cancer [odds ratio (OR) = 1.10, 95% confidence interval (95% CI) = 1.05-1.14]. Another significant tagSNP rs12212247 T>C (P combined = 9.589 × 10(-6)) was associated with a decreased risk of lung cancer (OR = 0.93, 95% CI = 0.90-0.96). Further in silico functional analyzes revealed that rs151606 might affect transcriptional regulation and result in decreased FGFR1OP expression (P trend = 0.022). The findings shed some new light on the role of centrosome abnormalities in the susceptibility to lung carcinogenesis.

Identification of lung cancer histology-specific variants applying Bayesian framework variant prioritization approaches within the TRICL and ILCCO consortia.

Large-scale genome-wide association studies (GWAS) have likely uncovered all common variants at the GWAS significance level. Additional variants within the suggestive range (0.0001> P > 5×10(-8)) are, however, still of interest for identifying causal associations. This analysis aimed to apply novel variant prioritization approaches to identify additional lung cancer variants that may not reach the GWAS level. Effects were combined across studies with a total of 33456 controls and 6756 adenocarcinoma (AC; 13 studies), 5061 squamous cell carcinoma (SCC; 12 studies) and 2216 small cell lung cancer cases (9 studies). Based on prior information such as variant physical properties and functional significance, we applied stratified false discovery rates, hierarchical modeling and Bayesian false discovery probabilities for variant prioritization. We conducted a fine mapping analysis as validation of our methods by examining top-ranking novel variants in six independent populations with a total of 3128 cases and 2966 controls. Three novel loci in the suggestive range were identified based on our Bayesian framework analyses: KCNIP4 at 4p15.2 (rs6448050, P = 4.6×10(-7)) and MTMR2 at 11q21 (rs10501831, P = 3.1×10(-6)) with SCC, as well as GAREM at 18q12.1 (rs11662168, P = 3.4×10(-7)) with AC. Use of our prioritization methods validated two of the top three loci associated with SCC (P = 1.05×10(-4) for KCNIP4, represented by rs9799795) and AC (P = 2.16×10(-4) for GAREM, represented by rs3786309) in the independent fine mapping populations. This study highlights the utility of using prior functional data for sequence variants in prioritization analyses to search for robust signals in the suggestive range.

Cannabis smoking and lung cancer risk: Pooled analysis in the International Lung Cancer Consortium.

To investigate the association between cannabis smoking and lung cancer risk, data on 2,159 lung cancer cases and 2,985 controls were pooled from 6 case-control studies in the US, Canada, UK, and New Zealand within the International Lung Cancer Consortium. Study-specific associations between cannabis smoking and lung cancer were estimated using unconditional logistic regression adjusting for sociodemographic factors, tobacco smoking status and pack-years; odds-ratio estimates were pooled using random effects models. Subgroup analyses were done for sex, histology and tobacco smoking status. The shapes of dose-response associations were examined using restricted cubic spline regression. The overall pooled OR for habitual versus nonhabitual or never users was 0.96 (95% CI: 0.66-1.38). Compared to nonhabitual or never users, the summary OR was 0.88 (95%CI: 0.63-1.24) for individuals who smoked 1 or more joint-equivalents of cannabis per day and 0.94 (95%CI: 0.67-1.32) for those consumed at least 10 joint-years. For adenocarcinoma cases the ORs were 1.73 (95%CI: 0.75-4.00) and 1.74 (95%CI: 0.85-3.55), respectively. However, no association was found for the squamous cell carcinoma based on small numbers. Weak associations between cannabis smoking and lung cancer were observed in never tobacco smokers. Spline modeling indicated a weak positive monotonic association between cumulative cannabis use and lung cancer, but precision was low at high exposure levels. Results from our pooled analyses provide little evidence for an increased risk of lung cancer among habitual or long-term cannabis smokers, although the possibility of potential adverse effect for heavy consumption cannot be excluded.

Cancer patients' acceptability of incorporating an epidemiology questionnaire within a clinical trial.

BACKGROUND/AIMS: Understanding the influence or impact of epidemiological factors on cancer outcomes in clinical trials can broaden our knowledge of disease, trial populations and therapeutic effects thus leading to improved patient care. However, there is a lack of data on cancer patients' compliance with an epidemiology questionnaire in the context of a clinical trial. PATIENTS AND METHODS: Cancer patients were provided with a hypothetical scenario and surveyed regarding their willingness and preferences to complete an epidemiology questionnaire if incorporated into a cancer therapy trial. Patient compliance with completing a voluntary epidemiology questionnaire and trial coordinators perceptions therein were separately determined in the NCIC Clinical Trials Group HN.6 clinical trial, an ongoing randomized phase III trial comparing two first-line treatment regimens in patients with locoregionally advanced head and neck cancer. RESULTS: Of 617 cancer patients from community, academic and tertiary cancer centres, the majority were willing to complete an epidemiology questionnaire either unconditionally (45%), or provided it did not inconvenience them (31%); 4% would refuse. Patients preferred shorter questionnaires of 30-50 questions requiring 10-20 min to complete, administered over 1-3 sessions. Patients were less willing, but still compliant, to answer questions relating to sexual history (71%) and annual household income (66%) relative to other questions (>90%). Eighteen percent thought that the questionnaire should be mandatory, with 31% believing that they may benefit personally from such research. In the HN.6 trial, compliance averaged 94.8% per question. CONCLUSIONS: Cancer patients are very willing to complete epidemiology questions in clinical trials.

Circulating proteome for pulmonary nodule malignancy.

BACKGROUND: Although lung cancer screening with low-dose computed tomography is rolling out in many areas of the world, differentiating indeterminate pulmonary nodules remains a major challenge. We conducted one of the first systematic investigations of circulating protein markers to differentiate malignant from benign screen-detected pulmonary nodules. METHODS: Based on 4 international low-dose computed tomography screening studies, we assayed 1078 protein markers using prediagnostic blood samples from 1253 participants based on a nested case-control design. Protein markers were measured using proximity extension assays, and data were analyzed using multivariable logistic regression, random forest, and penalized regressions. Protein burden scores (PBSs) for overall nodule malignancy and imminent tumors were estimated. RESULTS: We identified 36 potentially informative circulating protein markers differentiating malignant from benign nodules, representing a tightly connected biological network. Ten markers were found to be particularly relevant for imminent lung cancer diagnoses within 1 year. Increases in PBSs for overall nodule malignancy and imminent tumors by 1 standard deviation were associated with odds ratios of 2.29 (95% confidence interval: 1.95 to 2.72) and 2.81 (95% confidence interval: 2.27 to 3.54) for nodule malignancy overall and within 1 year of diagnosis, respectively. Both PBSs for overall nodule malignancy and for imminent tumors were substantially higher for those with malignant nodules than for those with benign nodules, even when limited to Lung Computed Tomography Screening Reporting and Data System (LungRADS) category 4 (P < .001). CONCLUSIONS: Circulating protein markers can help differentiate malignant from benign pulmonary nodules. Validation with an independent computed tomographic screening study will be required before clinical implementation.

Design and methodological considerations for biomarker discovery and validation in the Integrative Analysis of Lung Cancer Etiology and Risk (INTEGRAL) Program.

The Integrative Analysis of Lung Cancer Etiology and Risk (INTEGRAL) program is an NCI-funded initiative with an objective to develop tools to optimize low-dose CT (LDCT) lung cancer screening. Here, we describe the rationale and design for the Risk Biomarker and Nodule Malignancy projects within INTEGRAL. The overarching goal of these projects is to systematically investigate circulating protein markers to include on a panel for use (i) pre-LDCT, to identify people likely to benefit from screening, and (ii) post-LDCT, to differentiate benign versus malignant nodules. To identify informative proteins, the Risk Biomarker project measured 1161 proteins in a nested-case control study within 2 prospective cohorts (n = 252 lung cancer cases and 252 controls) and replicated associations for a subset of proteins in 4 cohorts (n = 479 cases and 479 controls). Eligible participants had a current or former history of smoking and cases were diagnosed up to 3 years following blood draw. The Nodule Malignancy project measured 1078 proteins among participants with a heavy smoking history within four LDCT screening studies (n = 425 cases diagnosed up to 5 years following blood draw, 430 benign-nodule controls, and 398 nodule-free controls). The INTEGRAL panel will enable absolute quantification of 21 proteins. We will evaluate its performance in the Risk Biomarker project using a case-cohort study including 14 cohorts (n = 1696 cases and 2926 subcohort representatives), and in the Nodule Malignancy project within five LDCT screening studies (n = 675 cases, 680 benign-nodule controls, and 648 nodule-free controls). Future progress to advance lung cancer early detection biomarkers will require carefully designed validation, translational, and comparative studies.

Deciphering associations between three RNA splicing-related genetic variants and lung cancer risk.

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.

Assessing Lung Cancer Absolute Risk Trajectory Based on a Polygenic Risk Model.

Lung cancer is the leading cause of cancer-related death globally. An improved risk stratification strategy can increase efficiency of low-dose CT (LDCT) screening. Here we assessed whether individual's genetic background has clinical utility for risk stratification in the context of LDCT screening. On the basis of 13,119 patients with lung cancer and 10,008 controls with European ancestry in the International Lung Cancer Consortium, we constructed a polygenic risk score (PRS) via 10-fold cross-validation with regularized penalized regression. The performance of risk model integrating PRS, including calibration and ability to discriminate, was assessed using UK Biobank data (N = 335,931). Absolute risk was estimated on the basis of age-specific lung cancer incidence and all-cause mortality as competing risk. To evaluate its potential clinical utility, the PRS distribution was simulated in the National Lung Screening Trial (N = 50,772 participants). The lung cancer ORs for individuals at the top decile of the PRS distribution versus those at bottom 10% was 2.39 [95% confidence interval (CI) = 1.92-3.00; P = 1.80 × 10-14] in the validation set (P trend = 5.26 × 10-20). The OR per SD of PRS increase was 1.26 (95% CI = 1.20-1.32; P = 9.69 × 10-23) for overall lung cancer risk in the validation set. When considering absolute risks, individuals at different PRS deciles showed differential trajectories of 5-year and cumulative absolute risk. The age reaching the LDCT screening recommendation threshold can vary by 4 to 8 years, depending on the individual's genetic background, smoking status, and family history. Collectively, these results suggest that individual's genetic background may inform the optimal lung cancer LDCT screening strategy. SIGNIFICANCE: Three large-scale datasets reveal that, after accounting for risk factors, an individual's genetics can affect their lung cancer risk trajectory, thus may inform the optimal timing for LDCT screening.

Germline determinants of humoral immune response to HPV-16 protect against oropharyngeal cancer.

Although several oropharyngeal cancer (OPC) susceptibility loci have been identified, most previous studies lacked detailed information on human papillomavirus (HPV) status. We conduct a genome-wide analysis by HPV16 serology status in 4,002 oral cancer cases (OPC and oral cavity cancer (OCC)) and 5,256 controls. We detect four susceptibility loci pointing to a distinct genetic predisposition by HPV status. Our most notable finding in the HLA region, that is now confirmed to be specific of HPV(+)OPC risk, reveal two independent loci with strong protective effects, one refining the previously reported HLA class II haplotype association. Antibody levels against HPV16 viral proteins strongly implicate the protective HLA variants as major determinants of humoral response against L1 capsid protein or E6 oncoprotein suggesting a natural immune response against HPV(+)OPC promoted by HLA variants. This indicates that therapeutic vaccines that target E6 and attenuate viral response after established HPV infections might protect against HPV(+)OPC.