Genome-wide interaction analysis identified low-frequency variants with sex disparity in lung cancer risk.

Differences by sex in lung cancer incidence and mortality have been reported which cannot be fully explained by sex differences in smoking behavior, implying existence of genetic and molecular basis for sex disparity in lung cancer development. However, the information about sex dimorphism in lung cancer risk is quite limited despite the great success in lung cancer association studies. By adopting a stringent two-stage analysis strategy, we performed a genome-wide gene-sex interaction analysis using genotypes from a lung cancer cohort including ~ 47 000 individuals with European ancestry. Three low-frequency variants (minor allele frequency < 0.05), rs17662871 [odds ratio (OR) = 0.71, P = 4.29×10-8); rs79942605 (OR = 2.17, P = 2.81×10-8) and rs208908 (OR = 0.70, P = 4.54×10-8) were identified with different risk effect of lung cancer between men and women. Further expression quantitative trait loci and functional annotation analysis suggested rs208908 affects lung cancer risk through differential regulation of Coxsackie virus and adenovirus receptor gene expression in lung tissues between men and women. Our study is one of the first studies to provide novel insights about the genetic and molecular basis for sex disparity in lung cancer development.

Identification of genetically predicted DNA methylation markers associated with non-small cell lung cancer risk among 34,964 cases and 448,579 controls.

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.

Genome-wide association meta-analysis identifies pleiotropic risk loci for aerodigestive squamous cell cancers.

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.

Iam hiQ-a novel pair of accuracy indices for imputed genotypes.

BACKGROUND: Imputation of untyped markers is a standard tool in genome-wide association studies to close the gap between directly genotyped and other known DNA variants. However, high accuracy with which genotypes are imputed is fundamental. Several accuracy measures have been proposed and some are implemented in imputation software, unfortunately diversely across platforms. In the present paper, we introduce Iam hiQ, an independent pair of accuracy measures that can be applied to dosage files, the output of all imputation software. Iam (imputation accuracy measure) quantifies the average amount of individual-specific versus population-specific genotype information in a linear manner. hiQ (heterogeneity in quantities of dosages) addresses the inter-individual heterogeneity between dosages of a marker across the sample at hand. RESULTS: Applying both measures to a large case-control sample of the International Lung Cancer Consortium (ILCCO), comprising 27,065 individuals, we found meaningful thresholds for Iam and hiQ suitable to classify markers of poor accuracy. We demonstrate how Manhattan-like plots and moving averages of Iam and hiQ can be useful to identify regions enriched with less accurate imputed markers, whereas these regions would by missed when applying the accuracy measure info (implemented in IMPUTE2). CONCLUSION: We recommend using Iam hiQ additional to other accuracy scores for variant filtering before stepping into the analysis of imputed GWAS data.

ImputAccur: fast and user-friendly calculation of genotype-imputation accuracy-measures.

BACKGROUND: ImputAccur is a software tool to measure genotype-imputation accuracy. Imputation of untyped markers is a standard approach in genome-wide association studies to close the gap between directly genotyped and other known DNA variants. However, high accuracy for imputed genotypes is fundamental. Several accuracy measures have been proposed, but unfortunately, they are implemented on different platforms, which is impractical. RESULTS: With ImputAccur, the accuracy measures info, Iam-hiQ and r2-based indices can be derived from standard output files of imputation software. Sample/probe and marker filtering is possible. This allows e.g. accurate marker filtering ahead of data analysis. CONCLUSIONS: The source code (Python version 3.9.4), a standalone executive file, and example data for ImputAccur are freely available at https://gitlab.gwdg.de/kolja.thormann1/imputationquality.git .

Integration of multiomic annotation data to prioritize and characterize inflammation and immune-related risk variants in squamous cell lung cancer.

Clinical trial results have recently demonstrated that inhibiting inflammation by targeting the interleukin-1ß pathway can offer a significant reduction in lung cancer incidence and mortality, highlighting a pressing and unmet need to understand the benefits of inflammation-focused lung cancer therapies at the genetic level. While numerous genome-wide association studies (GWAS) have explored the genetic etiology of lung cancer, there remains a large gap between the type of information that may be gleaned from an association study and the depth of understanding necessary to explain and drive translational findings. Thus, in this study we jointly model and integrate extensive multiomics data sources, utilizing a total of 40 genome-wide functional annotations that augment previously published results from the International Lung Cancer Consortium (ILCCO) GWAS, to prioritize and characterize single nucleotide polymorphisms (SNPs) that increase risk of squamous cell lung cancer through the inflammatory and immune responses. Our work bridges the gap between correlative analysis and translational follow-up research, refining GWAS association measures in an interpretable and systematic manner. In particular, reanalysis of the ILCCO data highlights the impact of highly associated SNPs from nuclear factor-κB signaling pathway genes as well as major histocompatibility complex mediated variation in immune responses. One consequence of prioritizing likely functional SNPs is the pruning of variants that might be selected for follow-up work by over an order of magnitude, from potentially tens of thousands to hundreds. The strategies we introduce provide informative and interpretable approaches for incorporating extensive genome-wide annotation data in analysis of genetic association studies.

Causal relationships between body mass index, smoking and lung cancer: Univariable and multivariable Mendelian randomization.

At the time of cancer diagnosis, body mass index (BMI) is inversely correlated with lung cancer risk, which may reflect reverse causality and confounding due to smoking behavior. We used two-sample univariable and multivariable Mendelian randomization (MR) to estimate causal relationships of BMI and smoking behaviors on lung cancer and histological subtypes based on an aggregated genome-wide association studies (GWASs) analysis of lung cancer in 29 266 cases and 56 450 controls. We observed a positive causal effect for high BMI on occurrence of small-cell lung cancer (odds ratio (OR) = 1.60, 95% confidence interval (CI) = 1.24-2.06, P = 2.70 × 10-4 ). After adjustment of smoking behaviors using multivariable Mendelian randomization (MVMR), a direct causal effect on small cell lung cancer (ORMVMR = 1.28, 95% CI = 1.06-1.55, PMVMR = .011), and an inverse effect on lung adenocarcinoma (ORMVMR = 0.86, 95% CI = 0.77-0.96, PMVMR = .008) were observed. A weak increased risk of lung squamous cell carcinoma was observed for higher BMI in univariable Mendelian randomization (UVMR) analysis (ORUVMR = 1.19, 95% CI = 1.01-1.40, PUVMR = .036), but this effect disappeared after adjustment of smoking (ORMVMR = 1.02, 95% CI = 0.90-1.16, PMVMR = .746). These results highlight the histology-specific impact of BMI on lung carcinogenesis and imply mediator role of smoking behaviors in the association between BMI and lung cancer.

Genome-wide association study of INDELs identified four novel susceptibility loci associated with lung cancer risk.

Genome-wide association studies (GWAS) have identified 45 susceptibility loci associated with lung cancer. Only less than SNPs, small insertions and deletions (INDELs) are the second most abundant genetic polymorphisms in the human genome. INDELs are highly associated with multiple human diseases, including lung cancer. However, limited studies with large-scale samples have been available to systematically evaluate the effects of INDELs on lung cancer risk. Here, we performed a large-scale meta-analysis to evaluate INDELs and their risk for lung cancer in 23,202 cases and 19,048 controls. Functional annotations were performed to further explore the potential function of lung cancer risk INDELs. Conditional analysis was used to clarify the relationship between INDELs and SNPs. Four new risk loci were identified in genome-wide INDEL analysis (1p13.2: rs5777156, Insertion, OR = 0.92, p = 9.10 × 10-8 ; 4q28.2: rs58404727, Deletion, OR = 1.19, p = 5.25 × 10-7 ; 12p13.31: rs71450133, Deletion, OR = 1.09, p = 8.83 × 10-7 ; and 14q22.3: rs34057993, Deletion, OR = 0.90, p = 7.64 × 10-8 ). The eQTL analysis and functional annotation suggested that INDELs might affect lung cancer susceptibility by regulating the expression of target genes. After conducting conditional analysis on potential causal SNPs, the INDELs in the new loci were still nominally significant. Our findings indicate that INDELs could be potentially functional genetic variants for lung cancer risk. Further functional experiments are needed to better understand INDEL mechanisms in carcinogenesis.

Transcriptome-wide association study reveals candidate causal genes for lung cancer.

We have recently completed the largest GWAS on lung cancer including 29,266 cases and 56,450 controls of European descent. The goal of our study has been to integrate the complete GWAS results with a large-scale expression quantitative trait loci (eQTL) mapping study in human lung tissues (n = 1,038) to identify candidate causal genes for lung cancer. We performed transcriptome-wide association study (TWAS) for lung cancer overall, by histology (adenocarcinoma, squamous cell carcinoma and small cell lung cancer) and smoking subgroups (never- and ever-smokers). We performed replication analysis using lung data from the Genotype-Tissue Expression (GTEx) project. DNA damage assays were performed in human lung fibroblasts for selected TWAS genes. As expected, the main TWAS signal for all histological subtypes and ever-smokers was on chromosome 15q25. The gene most strongly associated with lung cancer at this locus using the TWAS approach was IREB2 (pTWAS = 1.09E-99), where lower predicted expression increased lung cancer risk. A new lung adenocarcinoma susceptibility locus was revealed on 9p13.3 and associated with higher predicted expression of AQP3 (pTWAS = 3.72E-6). Among the 45 previously described lung cancer GWAS loci, we mapped candidate target gene for 17 of them. The association AQP3-adenocarcinoma on 9p13.3 was replicated using GTEx (pTWAS = 6.55E-5). Consistent with the effect of risk alleles on gene expression levels, IREB2 knockdown and AQP3 overproduction promote endogenous DNA damage. These findings indicate genes whose expression in lung tissue directly influences lung cancer risk.

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.

Detecting responses to treatment with fenofibrate in pedigrees.

BACKGROUND: Fenofibrate (Fb) is a known treatment for elevated triglyceride (TG) levels. The Genetics of Lipid Lowering Drugs and Diet Network (GOLDN) study was designed to investigate potential contributors to the effects of Fb on TG levels. Here, we summarize the analyses of 8 papers whose authors had access to the GOLDN data and were grouped together because they pursued investigations into Fb treatment responses as part of GAW20. These papers report explorations of a variety of genetics, epigenetics, and study design questions. Data regarding treatment with 160 mg of micronized Fb per day for 3 weeks included pretreatment and posttreatment TG and methylation levels (ML) at approximately 450,000 epigenetic markers (cytosine-phosphate-guanine [CpG] sites). In addition, approximately 1 million single-nucleotide polymorphisms (SNPs) were genotyped or imputed in each of the study participants, drawn from 188 pedigrees. RESULTS: The analyses of a variety of subsets of the GOLDN data used a number of analytic approaches such as linear mixed models, a kernel score test, penalized regression, and artificial neural networks. CONCLUSIONS: Results indicate that (a) CpG ML are responsive to Fb; (b) CpG ML should be included in models predicting the TG level responses to Fb;

Genetic modifiers of radon-induced lung cancer risk: a genome-wide interaction study in former uranium miners.

PURPOSE: Radon is a risk factor for lung cancer and uranium miners are more exposed than the general population. A genome-wide interaction analysis was carried out to identify genomic loci, genes or gene sets that modify the susceptibility to lung cancer given occupational exposure to the radioactive gas radon. METHODS: Samples from 28 studies provided by the International Lung Cancer Consortium were pooled with samples of former uranium miners collected by the German Federal Office of Radiation Protection. In total, 15,077 cases and 13,522 controls, all of European ancestries, comprising 463 uranium miners were compared. The DNA of all participants was genotyped with the OncoArray. We fitted single-marker and in multi-marker models and performed an exploratory gene-set analysis to detect cumulative enrichment of significance in sets of genes. RESULTS: We discovered a genome-wide significant interaction of the marker rs12440014 within the gene CHRNB4 (OR = 0.26, 95% CI 0.11-0.60, p = 0.0386 corrected for multiple testing). At least suggestive significant interaction of linkage disequilibrium blocks was observed at the chromosomal regions 18q21.23 (p = 1.2 × 10-6), 5q23.2 (p = 2.5 × 10-6), 1q21.3 (p = 3.2 × 10-6), 10p13 (p = 1.3 × 10-5) and 12p12.1 (p = 7.1 × 10-5). Genes belonging to the Gene Ontology term "DNA dealkylation involved in DNA repair" (GO:0006307; p = 0.0139) or the gene family HGNC:476 "microRNAs" (p = 0.0159) were enriched with LD-blockwise significance. CONCLUSION: The well-established association of the genomic region 15q25 to lung cancer might be influenced by exposure to radon among uranium miners. Furthermore, lung cancer susceptibility is related to the functional capability of DNA damage signaling via ubiquitination processes and repair of radiation-induced double-strand breaks by the single-strand annealing mechanism.

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.

Inherited variation in circadian rhythm genes and risks of prostate cancer and three other cancer sites in combined cancer consortia.

Circadian disruption has been linked to carcinogenesis in animal models, but the evidence in humans is inconclusive. Genetic variation in circadian rhythm genes provides a tool to investigate such associations. We examined associations of genetic variation in nine core circadian rhythm genes and six melatonin pathway genes with risk of colorectal, lung, ovarian and prostate cancers using data from the Genetic Associations and Mechanisms in Oncology (GAME-ON) network. The major results for prostate cancer were replicated in the Prostate, Lung, Colorectal and Ovarian (PLCO) cancer screening trial, and for colorectal cancer in the Genetics and Epidemiology of Colorectal Cancer Consortium (GECCO). The total number of cancer cases and controls was 15,838/18,159 for colorectal, 14,818/14,227 for prostate, 12,537/17,285 for lung and 4,369/9,123 for ovary. For each cancer site, we conducted gene-based and pathway-based analyses by applying the summary-based Adaptive Rank Truncated Product method (sARTP) on the summary association statistics for each SNP within the candidate gene regions. Aggregate genetic variation in circadian rhythm and melatonin pathways were significantly associated with the risk of prostate cancer in data combining GAME-ON and PLCO, after Bonferroni correction (ppathway < 0.00625). The two most significant genes were NPAS2 (pgene = 0.0062) and AANAT (pgene = 0.00078); the latter being significant after Bonferroni correction. For colorectal cancer, we observed a suggestive association with the circadian rhythm pathway in GAME-ON (ppathway = 0.021); this association was not confirmed in GECCO (ppathway = 0.76) or the combined data (ppathway = 0.17). No significant association was observed for ovarian and lung cancer. These findings support a potential role for circadian rhythm and melatonin pathways in prostate carcinogenesis. Further functional studies are needed to better understand the underlying biologic mechanisms.

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.

Genetic determinants of telomere length and risk of common cancers: a Mendelian randomization study.

Epidemiological studies have reported inconsistent associations between telomere length (TL) and risk for various cancers. These inconsistencies are likely attributable, in part, to biases that arise due to post-diagnostic and post-treatment TL measurement. To avoid such biases, we used a Mendelian randomization approach and estimated associations between nine TL-associated SNPs and risk for five common cancer types (breast, lung, colorectal, ovarian and prostate cancer, including subtypes) using data on 51 725 cases and 62 035 controls. We then used an inverse-variance weighted average of the SNP-specific associations to estimate the association between a genetic score representing long TL and cancer risk. The long TL genetic score was significantly associated with increased risk of lung adenocarcinoma (P = 6.3 × 10(-15)), even after exclusion of a SNP residing in a known lung cancer susceptibility region (TERT-CLPTM1L) P = 6.6 × 10(-6)). Under Mendelian randomization assumptions, the association estimate [odds ratio (OR) = 2.78] is interpreted as the OR for lung adenocarcinoma corresponding to a 1000 bp increase in TL. The weighted TL SNP score was not associated with other cancer types or subtypes. Our finding that genetic determinants of long TL increase lung adenocarcinoma risk avoids issues with reverse causality and residual confounding that arise in observational studies of TL and disease risk. Under Mendelian randomization assumptions, our finding suggests that longer TL increases lung adenocarcinoma risk. However, caution regarding this causal interpretation is warranted in light of the potential issue of pleiotropy, and a more general interpretation is that SNPs influencing telomere biology are also implicated in lung adenocarcinoma risk.

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.

Genetic variant in DNA repair gene GTF2H4 is associated with lung cancer risk: a large-scale analysis of six published GWAS datasets in the TRICL consortium.

DNA repair pathways maintain genomic integrity and stability, and dysfunction of DNA repair leads to cancer. We hypothesize that functional genetic variants in DNA repair genes are associated with risk of lung cancer. We performed a large-scale meta-analysis of 123,371 single nucleotide polymorphisms (SNPs) in 169 DNA repair genes obtained from six previously published genome-wide association studies (GWASs) of 12160 lung cancer cases and 16838 controls. We calculated odds ratios (ORs) with 95% confidence intervals (CIs) using the logistic regression model and used the false discovery rate (FDR) method for correction of multiple testing. As a result, 14 SNPs had a significant odds ratio (OR) for lung cancer risk with P FDR < 0.05, of which rs3115672 in MSH5 (OR = 1.20, 95% CI = 1.14-1.27) and rs114596632 in GTF2H4 (OR = 1.19, 95% CI = 1.12-1.25) at 6q21.33 were the most statistically significant (P combined = 3.99×10(-11) and P combined = 5.40×10(-10), respectively). The MSH5 rs3115672, but not GTF2H4 rs114596632, was strongly correlated with MSH5 rs3131379 in that region (r (2) = 1.000 and r (2) = 0.539, respectively) as reported in a previous GWAS. Importantly, however, the GTF2H4 rs114596632 T, but not MSH5 rs3115672 T, allele was significantly associated with both decreased DNA repair capacity phenotype and decreased mRNA expression levels. These provided evidence that functional genetic variants of GTF2H4 confer susceptibility to lung cancer.

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.

Informed genome-wide association analysis with family history as a secondary phenotype identifies novel loci of lung cancer.

Lung cancer is the leading cause of cancer death worldwide. Although several genetic variants associated with lung cancer have been identified in the past, stringent selection criteria of genome-wide association studies (GWAS) can lead to missed variants. The objective of this study was to uncover missed variants by using the known association between lung cancer and first-degree family history of lung cancer to enrich the variant prioritization for lung cancer susceptibility regions. In this two-stage GWAS study, we first selected a list of variants associated with both lung cancer and family history of lung cancer in four GWAS (3,953 cases, 4,730 controls), then replicated our findings for 30 variants in a meta-analysis of four additional studies (7,510 cases, 7,476 controls). The top ranked genetic variant rs12415204 in chr10q23.33 encoding FFAR4 in the Discovery set was validated in the Replication set with an overall OR of 1.09 (95% CI=1.04, 1.14, P=1.63×10(-4)). When combining the two stages of the study, the strongest association was found in rs1158970 at Ch4p15.2 encoding KCNIP4 with an OR of 0.89 (95% CI=0.85, 0.94, P=9.64×10(-6)). We performed a stratified analysis of rs12415204 and rs1158970 across all eight studies by age, gender, smoking status, and histology, and found consistent results across strata. Four of the 30 replicated variants act as expression quantitative trait loci (eQTL) sites in 1,111 nontumor lung tissues and meet the genome-wide 10% FDR threshold.