Brennan, Paul

Abstract There is growing evidence showing that sex disparity in lung cancer risk 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 genotype from a lung cancer cohort including ~ 47,000 individuals with European ancestry. Three low-frequency variants (minor allele frequency < 0.05), rs17662871 (OR=0.71, p=4.29x10-8); rs79942605 (OR=2.17, p=2.81x10-8); and rs208908 (OR=0.70, p=4.54x10-8), were identified with different risk effect of lung cancer between men and women. Further eQTL and functional annotation analysis suggested rs208908 affect lung cancer risk through differential regulation of CXADR (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.

Abstract 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 &lt; 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, Pinteraction = 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 (r2 = 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.

Abstract Background: Whether circulating polyunsaturated fatty acid (PUFA) levels are associated with pancreatic cancer risk is uncertain. Mendelian randomization (MR) represents a study design using genetic instruments to better characterize the relationship between exposure and outcome. Methods: We utilized data from genome-wide association studies within the Pancreatic Cancer Cohort Consortium and Pancreatic Cancer Case–Control Consortium, involving approximately 9,269 cases and 12,530 controls of European descent, to evaluate associations between pancreatic cancer risk and genetically predicted plasma n-6 PUFA levels. Conventional MR analyses were performed using individual-level and summary-level data. Results: Using genetic instruments, we did not find evidence of associations between genetically predicted plasma n-6 PUFA levels and pancreatic cancer risk [estimates per one SD increase in each PUFA-specific weighted genetic score using summary statistics: linoleic acid odds ratio (OR) = 1.00, 95% confidence interval (CI) = 0.98–1.02; arachidonic acid OR = 1.00, 95% CI = 0.99–1.01; and dihomo-gamma-linolenic acid OR = 0.95, 95% CI = 0.87–1.02]. The OR estimates remained virtually unchanged after adjustment for covariates, using individual-level data or summary statistics, or stratification by age and sex. Conclusions: Our results suggest that variations of genetically determined plasma n-6 PUFA levels are not associated with pancreatic cancer risk. Impact: These results suggest that modifying n-6 PUFA levels through food sources or supplementation may not influence risk of pancreatic cancer.

Abstract Background: Lung cancer remains the leading cause of cancer death worldwide, with 15% to 20% occurring in never smokers. To assess genetic determinants for prognosis among never smokers, we conducted a genome-wide investigation in the International Lung Cancer Consortium (ILCCO). Methods: Genomic and clinical data from 1,569 never-smoking patients with lung cancer of European ancestry from 10 ILCCO studies were included. HRs and 95% confidence intervals of overall survival were estimated. We assessed whether the associations were mediated through mRNA expression–based 1,553 normal lung tissues from the lung expression quantitative trait loci (eQTL) dataset and Genotype-Tissue Expression (GTEx). For cross-ethnicity generalization, we assessed the associations in a Japanese study (N = 887). Results: One locus at 13q22.2 was associated with lung adenocarcinoma survival at genome-wide level, with carriers of rs12875562-T allele exhibiting poor prognosis [HR = 1.71 (1.41–2.07), P = 3.60 × 10−8], and altered mRNA expression of LMO7DN in lung tissue (GTEx, P = 9.40 × 10−7; Lung eQTL dataset, P = 0.003). Furthermore, 2 of 11 independent loci that reached the suggestive significance level (P &lt; 10−6) were significant eQTL affecting mRNA expression of nearby genes in lung tissues, including CAPZB at 1p36.13 and UBAC1 at 9q34.3. One locus encoding NWD2/KIAA1239 at 4p14 showed associations in both European [HR = 0.50 (0.38–0.66), P = 6.92 × 10−7] and Japanese populations [HR = 0.79 (0.67–0.94), P = 0.007]. Conclusions: Based on the largest genomic investigation on the lung cancer prognosis of never smokers to date, we observed that lung cancer prognosis is affected by inherited genetic variants. Impact: We identified one locus near LMO7DN at genome-wide level and several potential prognostic genes with cis-effect on mRNA expression. Further functional genomics work is required to understand their role in tumor progression.

Abstract 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 &lt; 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.

Abstract Background: Obesity and diabetes are major modifiable risk factors for pancreatic cancer. Interactions between genetic variants and diabetes/obesity have not previously been comprehensively investigated in pancreatic cancer at the genome-wide level. Methods: We conducted a gene–environment interaction (GxE) analysis including 8,255 cases and 11,900 controls from four pancreatic cancer genome-wide association study (GWAS) datasets (Pancreatic Cancer Cohort Consortium I–III and Pancreatic Cancer Case Control Consortium). Obesity (body mass index ≥30 kg/m2) and diabetes (duration ≥3 years) were the environmental variables of interest. Approximately 870,000 SNPs (minor allele frequency ≥0.005, genotyped in at least one dataset) were analyzed. Case–control (CC), case-only (CO), and joint-effect test methods were used for SNP-level GxE analysis. As a complementary approach, gene-based GxE analysis was also performed. Age, sex, study site, and principal components accounting for population substructure were included as covariates. Meta-analysis was applied to combine individual GWAS summary statistics. Results: No genome-wide significant interactions (departures from a log-additive odds model) with diabetes or obesity were detected at the SNP level by the CC or CO approaches. The joint-effect test detected numerous genome-wide significant GxE signals in the GWAS main effects top hit regions, but the significance diminished after adjusting for the GWAS top hits. In the gene-based analysis, a significant interaction of diabetes with variants in the FAM63A (family with sequence similarity 63 member A) gene (significance threshold P &lt; 1.25 × 10−6) was observed in the meta-analysis (PGxE = 1.2 ×10−6, PJoint = 4.2 ×10−7). Conclusions: This analysis did not find significant GxE interactions at the SNP level but found one significant interaction with diabetes at the gene level. A larger sample size might unveil additional genetic factors via GxE scans. Impact: This study may contribute to discovering the mechanism of diabetes-associated pancreatic cancer.

1512 Background: Tobacco smoking profoundly impacts LC risk; however, data are limited as to what extent SC prior to diagnosis impacts LC overall survival (OS) and lung cancer specific survival (LCSS). LC screening offers a possible teachable moment, but there is uncertainty of SC benefits after a lifetime of smoking. We use the ILCCO database to answer if SC prior to LC dx is associated with better OS and LCSS, considering time since smoking cessation (TSSC). Methods: Using individual data, analysis was performed on 17 ILCCO studies with available TSSC to estimate survival using univariable analysis and models of stage-adjusted and cumulative smoking-adjusted multivariable analysis. Adjusted Hazard Ratios (aHR) from Cox models, cubic spline smooth curves and Kaplan-Meier curves were created. Sensitivity analysis was performed for TSSC and LCSS on 13 studies. Results: Of 34649 patients, 14322 (41%) were current smokers 14273 (41%) ex-smokers and 6054 (18%) never smokers at diagnosis. We confirmed that ex-smokers (aHR 0.88 CI 0.86-0.91) and never smokers (aHR 0.76 CI 0.73-0.8) improved OS compared to current smokers. Amongst ex-smokers, < 2y TSSC (aHR 0.88 CI 0.82-0.94), 2-5y TSSC (aHR 0.83 CI 0.77-0.90) and > 5y TSSC (aHR 0.8 CI 0.76-0.84) had improved OS compared to CS. Sensitivity analysis showed a trend towards improved LCSS survival for < 2y TSSC (aHR 0.95 CI 0.86-1.05) and 2-5y TSSC (aHR 0.93 CI 0.83-1.04), whereas > 5y TSSC significantly improved LCSS by 15% (aHR 0.85 CI 0.78-0.92). To mimic the LC screening participants, in analysis of > 30 pack-years individuals, associations were strikingly strong: < 2y TSSC had improved OS by 14% (aHR 0.86 CI 0.80-0.93); 2-5y TSSC by 17% (aHR 0.83 CI 0.76-0.90); and > 5 TSSC by 22% (aHR 0.78 CI 0.74-0.83), compared to current smokers; for < 30 packs-years, a trend towards better OS was observed for < 2y TSSC (aHR 0.95 CI 0.92-1.02) and 2-5y TSSC (aHR 0.86 CI 0.74-1.01), whereas > 5y TSSC improved OS by 23% (aHR 0.77 CI 0.72-0.82). Conclusions: Among ex-smokers, the risk of overall death was reduced by 12% on < 2y TSSC, 17% on 2-5y TSSC and 20% > 5y TSSC, whereas for LCSS, the benefit was significant only for > 5y TSCC, compared to current smokers at time of diagnosis. Here we demonstrate that convincing screening participants to quit smoking at any point of their trajectory, even just prior to dx such as < 2y TSSC, improved OS, and LCSS benefit was present beyond 5y of quitting. These relationships are independent of pack-years, age, across all stages and other prognostic variables.

Abstract Background Although 20 pancreatic cancer susceptibility loci have been identified through genome-wide association studies in individuals of European ancestry, much of its heritability remains unexplained and the genes responsible largely unknown. Methods To discover novel pancreatic cancer risk loci and possible causal genes, we performed a pancreatic cancer transcriptome-wide association study in Europeans using three approaches: FUSION, MetaXcan, and Summary-MulTiXcan. We integrated genome-wide association studies summary statistics from 9040 pancreatic cancer cases and 12 496 controls, with gene expression prediction models built using transcriptome data from histologically normal pancreatic tissue samples (NCI Laboratory of Translational Genomics [n = 95] and Genotype-Tissue Expression v7 [n = 174] datasets) and data from 48 different tissues (Genotype-Tissue Expression v7, n = 74–421 samples). Results We identified 25 genes whose genetically predicted expression was statistically significantly associated with pancreatic cancer risk (false discovery rate &lt; .05), including 14 candidate genes at 11 novel loci (1p36.12: CELA3B; 9q31.1: SMC2, SMC2-AS1; 10q23.31: RP11-80H5.9; 12q13.13: SMUG1; 14q32.33: BTBD6; 15q23: HEXA; 15q26.1: RCCD1; 17q12: PNMT, CDK12, PGAP3; 17q22: SUPT4H1; 18q11.22: RP11-888D10.3; and 19p13.11: PGPEP1) and 11 at six known risk loci (5p15.33: TERT, CLPTM1L, ZDHHC11B; 7p14.1: INHBA; 9q34.2: ABO; 13q12.2: PDX1; 13q22.1: KLF5; and 16q23.1: WDR59, CFDP1, BCAR1, TMEM170A). The association for 12 of these genes (CELA3B, SMC2, and PNMT at novel risk loci and TERT, CLPTM1L, INHBA, ABO, PDX1, KLF5, WDR59, CFDP1, and BCAR1 at known loci) remained statistically significant after Bonferroni correction. Conclusions By integrating gene expression and genotype data, we identified novel pancreatic cancer risk loci and candidate functional genes that warrant further investigation.

Abstract Background Recent genome wide association studies (GWAS) of Hodgkin lymphoma (HL) have identified several associations at both HLA and non-HLA loci. However, much of HL heritability remains unexplained. Methods To identify novel risk loci, we performed a meta-analysis of 3 HL GWAS including a total of 1,810 cases and 7,879 controls. Results were replicated in an independent set of 1,163 cases and 2,580 controls, for a total of 3,097 and 11,097 cases and controls combined, respectively. participants in discovery and replication stages were of European descent. quality control and imputation we conducted a meta-analysis addressing 1,004,829 variants (λ= 1.10, λ1000= 1.03). Associations between SNP genotypes and HL risk were evaluated under a log-additive model of inheritance adjusting for sex, study center and significant principal components to control for population stratification. We performed an analysis with all HL cases and then conducted stratified analyses by histological subtype (classical, nodular sclerosis and mixed cellularity), age at diagnosis (nodular sclerosis among those diagnosed at 15- 35 years in all studies, and those diagnosed at 35 and older years in the European Study only) and EBV tumor status (negative and positive). We then used a bioinformatic approach (FunciSNP) to identify potential functional variants associated with HD risk correlated with risk loci of interest. We extracted publically available ENCODE data on biofeatures to identify potential functional motifs associated with the index SNP or correlated SNPs. Finally, we measured expression levels of the two alternative mRNA transcripts in lymphoblastoid cell lines (LCLs) from 49 post-therapy HL patients and 25 unaffected controls. RT-PCR was carried out in triplicate. Relative expression levels were calculated relative to TBP as housekeeping gene. Linear models were used to assess correlation between genotype and TCF3expression levels. Results The meta-analysis identified a novel susceptibility variant at chromosome 19p13.3 (rs1860661) associated with HL risk (Odds Ratio [OR]= 0.78, P=2.0*10-8, I2=0%). variant is located in intron 2 of TCF3 (also known as E2A), a regulator of B- and T-cell lineage commitment. was also significantly associated with HL (OR= 0.85, P=0.002) in the replication series of 1,281 cases and 3,218 controls. the combined analysis consisting of the discovery and replication sets, rs1860661was strongly associated with HL (OR=0.81,=3.5*10-10), with no evidence of heterogeneity between contributing studies (Phom=0.41, I2=0%). The number of G alleles defined by rs1860661 was significantly associated with a reduced risk of each HL subtype except EBV positive HL. rs1860661 and two correlated SNPs, rs10413888 (r2=0.90) and rs8103453 (r2=0.89) identified by FunciSNP analysis map in or near marks of open chromatin and in DNAse hypersensitivity sites in TCF3 in CD20+ B cell lines., the protective minor alleles of these SNPs as defined by the G-G-G haplotype map to the binding sites for ZBTB7a (rs10413888 and rs1860661) and (rs8103453) transcription factors, likely improving the binding efficiency to the sites which may result in increased transcription rates of TCF3. TCF3 is encoded by two alternative transcripts (E12 and E47). Higher expression levels of TCF3-E47, whose transcription start site is located close to rs1860661, was associated with the rs1860661-G allele in controls (P=0.02), but not in HL patients (P=0.22). Conclusion/Discussion TCF3 is essential for the commitment of lymphoid progenitors to both B-cell and T-cell lineage development. A molecular and phenotypic hallmark of classical HL is the loss of the B-cell phenotype in HRS cells, including lack of demonstrable B-cell receptor and most B-cell specific markers such as CD19 or CD20. HRS cells have a low level of TCF3, particularly homodimers of the isoform E47, due to expression of the ABF-1 and ID2 inhibitors that bind to TCF3. Thus, higher TCF3 levels in HRS precursor cells may lead to enhanced retention of the B cell phenotype, thereby conferring a protective effect. These data suggest a link between the 19p13.3 locus including TCF3 and HL risk, indicating that TCF3 could be relevant to HL etiology and pathogenesis. Disclosures: Link: Genentech: Consultancy; Millenium: Consultancy; Pharmacyclics: Consultancy; Spectrum: Consultancy.

AbstractEpidemiologic studies show an increased risk of non‐Hodgkin lymphoma (NHL) in patients with autoimmune disease (AD), due to a combination of shared environmental factors and/or genetic factors, or a causative cascade: chronic inflammation/antigen‐stimulation in one disease leads to another. Here we assess shared genetic risk in genome‐wide‐association‐studies (GWAS).Secondary analysis of GWAS of NHL subtypes (chronic lymphocytic leukemia, diffuse large B‐cell lymphoma, follicular lymphoma, and marginal zone lymphoma) and ADs (rheumatoid arthritis, systemic lupus erythematosus, and multiple sclerosis). Shared genetic risk was assessed by (a) description of regional genetic of overlap, (b) polygenic risk score (PRS), (c)"diseasome", (d)meta‐analysis.Descriptive analysis revealed few shared genetic factors between each AD and each NHL subtype. The PRS of ADs were not increased in NHL patients (nor vice versa). In the diseasome, NHLs shared more genetic etiology with ADs than solid cancers (p = .0041). A meta‐analysis (combing AD with NHL) implicated genes of apoptosis and telomere length.This GWAS‐based analysis four NHL subtypes and three ADs revealed few weakly‐associated shared loci, explaining little total risk. This suggests common genetic variation, as assessed by GWAS in these sample sizes, may not be the primary explanation for the link between these ADs and NHLs.