Edlund, Christopher K.

Abstract Background: Inherited susceptibility is an important contributor to colorectal cancer risk, and rare variants in key genes or pathways could account in part for the missing proportion of colorectal cancer heritability. Methods: We conducted an exome-wide association study including 2,327 cases and 2,966 controls of European ancestry from three large epidemiologic studies. Single variant associations were tested using logistic regression models, adjusting for appropriate study-specific covariates. In addition, we examined the aggregate effects of rare coding variation at the gene and pathway levels using Bayesian model uncertainty techniques. Results: In an exome-wide gene-level analysis, we identified ST6GALNAC2 as the top associated gene based on the Bayesian risk index (BRI) method [summary Bayes factor (BF)BRI = 2604.23]. A rare coding variant in this gene, rs139401613, was the top associated variant (P = 1.01 × 10–6) in an exome-wide single variant analysis. Pathway-level association analyses based on the integrative BRI (iBRI) method found extreme evidence of association with the DNA repair pathway (BFiBRI = 17852.4), specifically with the nonhomologous end joining (BFiBRI = 437.95) and nucleotide excision repair (BFiBRI = 36.96) subpathways. The iBRI method also identified RPA2, PRKDC, ERCC5, and ERCC8 as the top associated DNA repair genes (summary BFiBRI ≥ 10), with rs28988897, rs8178232, rs141369732, and rs201642761 being the most likely associated variants in these genes, respectively. Conclusions: We identified novel variants and genes associated with colorectal cancer risk and provided additional evidence for a role of DNA repair in colorectal cancer tumorigenesis. Impact: This study provides new insights into the genetic predisposition to colorectal cancer, which has potential for translation into improved risk prediction.

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.

Abstract Introduction The 2-3 fold excess risk of multiple myeloma (MM) among family members of cases suggests a heritable contribution to risk. Recently, a genome-wide association study (GWAS) identified two genome-wide significant and one promising novel loci associated with multiple myeloma risk. To confirm these associations and identify additional novel risk loci, we performed a four-center, genome-wide association meta-analysis. Methods A fixed effects model was used for the meta-analysis which included a total of 1248 cases and 1485 controls, all of European descent, genotyped and analyzed at four separate centers with samples contributed by 10 studies. After quality control and imputation using the 1000 Genomes Project, the analysis included ∼9.5 million variants (λ=1.024). Associations between (single nucleotide polymorphisms) SNP genotypes and MM risk were evaluated under a log-additive model of inheritance, with each study adjusting for age, sex, and up to 10 principal components to control for population stratification. Promising results were replicated in an independent set of 1587 cases and 1770 controls using TaqMan, for a total of 2835 and 3255 cases and controls, respectively, in a combined meta-analysis. Results The discovery meta-analysis did not reveal any genome-wide significant associations (defined as p<5 x 10-8). We used a novel pruning algorithm to identify the top 35 most promising single nucleotide polymorphisms (SNPs) to advance to replication. We successfully genotyped 22 SNPs in the replication set. In the combined discovery and replication meta-analysis, rs1345359 at 12q23.1 was the most strongly associated SNP (P=9 x 10-8, Table 1). The variant allele C was associated with reduced risk (odds ratio discovery set [OR]= 0.69, OR replication set = 0.78, OR combined = 0.74). A second locus at 20q13.2 (rs150220835), was associated with a two-fold increased risk (P=1.22 x 10-6), a borderline increased risk (P=0.0900) and 45% increased risk (P=2.44 X 10-5) in the discovery, replication, and combined analysis sets respectively (Table 1). We also confirmed the association between MM risk and two previously published SNPs (rs4487645, p=0.0007and rs105251, p=0.0044) (Broderick et al., Nat. Genet., 2011). The third previously suggested SNP (rs6746082) was of nominal significance (p=0.0517) in the meta-analysis. Discussion We confirmed the association between MM risk and two previously published SNPs and identified a possible association with a novel SNP in chromosome 12q23.1 (rs1345359). This SNP is not located in a gene nor associated with biofeatures in ENCODE, thus further examination of correlated SNPs is necessary to identify a functional SNP linked to this locus. We also found suggestive evidence for a second locus at 20q13.2 requiring additional replication. Larger studies would improve risk variant discovery for this rare hematologic malignancy. Disclosures: Wolf: Celgene: Honoraria, Research Funding; Millenium: Honoraria; Onyx: Honoraria. Anderson:Celgene, Millennium, BMS, Onyx: Membership on an entity’s Board of Directors or advisory committees; Acetylon, Oncopep: Scientific Founder , Scientific Founder Other.

Abstract Background Previous genome-wide association studies (GWAS) have identified 42 loci (P < 5 × 10−8) associated with risk of colorectal cancer (CRC). Expanded consortium efforts facilitating the discovery of additional susceptibility loci may capture unexplained familial risk. Methods We conducted a GWAS in European descent CRC cases and control subjects using a discovery–replication design, followed by examination of novel findings in a multiethnic sample (cumulative n = 163 315). In the discovery stage (36 948 case subjects/30 864 control subjects), we identified genetic variants with a minor allele frequency of 1% or greater associated with risk of CRC using logistic regression followed by a fixed-effects inverse variance weighted meta-analysis. All novel independent variants reaching genome-wide statistical significance (two-sided P < 5 × 10−8) were tested for replication in separate European ancestry samples (12 952 case subjects/48 383 control subjects). Next, we examined the generalizability of discovered variants in East Asians, African Americans, and Hispanics (12 085 case subjects/22 083 control subjects). Finally, we examined the contributions of novel risk variants to familial relative risk and examined the prediction capabilities of a polygenic risk score. All statistical tests were two-sided. Results The discovery GWAS identified 11 variants associated with CRC at P < 5 × 10−8, of which nine (at 4q22.2/5p15.33/5p13.1/6p21.31/6p12.1/10q11.23/12q24.21/16q24.1/20q13.13) independently replicated at a P value of less than .05. Multiethnic follow-up supported the generalizability of discovery findings. These results demonstrated a 14.7% increase in familial relative risk explained by common risk alleles from 10.3% (95% confidence interval [CI] = 7.9% to 13.7%; known variants) to 11.9% (95% CI = 9.2% to 15.5%; known and novel variants). A polygenic risk score identified 4.3% of the population at an odds ratio for developing CRC of at least 2.0. Conclusions This study provides insight into the architecture of common genetic variation contributing to CRC etiology and improves risk prediction for individualized screening.

Abstract Background: Common cancers develop through a multistep process often including inherited susceptibility. Collaboration among multiple institutions, and funding from multiple sources, has allowed the development of an inexpensive genotyping microarray, the OncoArray. The array includes a genome-wide backbone, comprising 230,000 SNPs tagging most common genetic variants, together with dense mapping of known susceptibility regions, rare variants from sequencing experiments, pharmacogenetic markers, and cancer-related traits. Methods: The OncoArray can be genotyped using a novel technology developed by Illumina to facilitate efficient genotyping. The consortium developed standard approaches for selecting SNPs for study, for quality control of markers, and for ancestry analysis. The array was genotyped at selected sites and with prespecified replicate samples to permit evaluation of genotyping accuracy among centers and by ethnic background. Results: The OncoArray consortium genotyped 447,705 samples. A total of 494,763 SNPs passed quality control steps with a sample success rate of 97% of the samples. Participating sites performed ancestry analysis using a common set of markers and a scoring algorithm based on principal components analysis. Conclusions: Results from these analyses will enable researchers to identify new susceptibility loci, perform fine-mapping of new or known loci associated with either single or multiple cancers, assess the degree of overlap in cancer causation and pleiotropic effects of loci that have been identified for disease-specific risk, and jointly model genetic, environmental, and lifestyle-related exposures. Impact: Ongoing analyses will shed light on etiology and risk assessment for many types of cancer. Cancer Epidemiol Biomarkers Prev; 26(1); 126–35. ©2016 AACR.

Abstract Background: Genome-wide association studies (GWAS) in European populations have identified genetic risk variants associated with multiple myeloma. Methods: We performed association testing of common variation in eight regions in 1,318 patients with multiple myeloma and 1,480 controls of European ancestry and 1,305 patients with multiple myeloma and 7,078 controls of African ancestry and conducted a meta-analysis to localize the signals, with epigenetic annotation used to predict functionality. Results: We found that variants in 7p15.3, 17p11.2, 22q13.1 were statistically significantly (P &lt; 0.05) associated with multiple myeloma risk in persons of African ancestry and persons of European ancestry, and the variant in 3p22.1 was associated in European ancestry only. In a combined African ancestry–European ancestry meta-analysis, variation in five regions (2p23.3, 3p22.1, 7p15.3, 17p11.2, 22q13.1) was statistically significantly associated with multiple myeloma risk. In 3p22.1, the correlated variants clustered within the gene body of ULK4. Correlated variants in 7p15.3 clustered around an enhancer at the 3′ end of the CDCA7L transcription termination site. A missense variant at 17p11.2 (rs34562254, Pro251Leu, OR, 1.32; P = 2.93 × 10−7) in TNFRSF13B encodes a lymphocyte-specific protein in the TNF receptor family that interacts with the NF-κB pathway. SNPs correlated with the index signal in 22q13.1 cluster around the promoter and enhancer regions of CBX7. Conclusions: We found that reported multiple myeloma susceptibility regions contain risk variants important across populations, supporting the use of multiple racial/ethnic groups with different underlying genetic architecture to enhance the localization and identification of putatively functional alleles. Impact: A subset of reported risk loci for multiple myeloma has consistent effects across populations and is likely to be functional. Cancer Epidemiol Biomarkers Prev; 25(12); 1609–18. ©2016 AACR.

Abstract Identifying genetic variants with pleiotropic associations can uncover common pathways influencing multiple cancers. We took a two-stage approach to conduct genome-wide association studies for lung, ovary, breast, prostate, and colorectal cancer from the GAME-ON/GECCO Network (61,851 cases, 61,820 controls) to identify pleiotropic loci. Findings were replicated in independent association studies (55,789 cases, 330,490 controls). We identified a novel pleiotropic association at 1q22 involving breast and lung squamous cell carcinoma, with eQTL analysis showing an association with ADAM15/THBS3 gene expression in lung. We also identified a known breast cancer locus CASP8/ALS2CR12 associated with prostate cancer, a known cancer locus at CDKN2B-AS1 with different variants associated with lung adenocarcinoma and prostate cancer, and confirmed the associations of a breast BRCA2 locus with lung and serous ovarian cancer. This is the largest study to date examining pleiotropy across multiple cancer-associated loci, identifying common mechanisms of cancer development and progression. Cancer Res; 76(17); 5103–14. ©2016 AACR.

Nodular sclerosing Hodgkin lymphoma (NSHL) is a distinct, highly heritable Hodgkin lymphoma subtype. We undertook a genome-wide meta-analysis of 393 European-origin adolescent/young adult NSHL patients and 3315 controls using the Illumina Human610-Quad Beadchip and Affymetrix Genome-Wide Human SNP Array 6.0. We identified 3 single nucleotide polymorphisms (SNPs) on chromosome 6p21.32 that were significantly associated with NSHL risk: rs9268542 (P = 5.35 × 10−10), rs204999 (P = 1.44 × 10−9), and rs2858870 (P = 1.69 × 10−8). We also confirmed a previously reported association in the same region, rs6903608 (P = 3.52 × 10−10). rs204999 and rs2858870 were weakly correlated (r2 = 0.257), and the remaining pairs of SNPs were not correlated (r2 < 0.1). In an independent set of 113 NSHL cases and 214 controls, 2 SNPs were significantly associated with NSHL and a third showed a comparable odds ratio (OR). These SNPs are found on 2 haplotypes associated with NSHL risk (rs204999-rs9268528-rs9268542-rs6903608-rs2858870; AGGCT, OR = 1.7, P = 1.71 × 10−6; GAATC, OR = 0.4, P = 1.16 × 10−4). All individuals with the GAATC haplotype also carried the HLA class II DRB1*0701 allele. In a separate analysis, the DRB1*0701 allele was associated with a decreased risk of NSHL (OR = 0.5, 95% confidence interval = 0.4, 0.7). These data support the importance of the HLA class II region in NSHL etiology.

AbstractOver 22,000 cases of ovarian cancer were diagnosed in 2007 in the United States, but only a fraction of them can be attributed to mutations in highly penetrant genes such as BRCA1. To determine whether low-penetrance genetic variants contribute to ovarian cancer risk, we genotyped 1,536 single nucleotide polymorphisms (SNP) in several candidate gene pathways in 848 epithelial ovarian cancer cases and 798 controls in the North Carolina Ovarian Cancer Study (NCO) using a customized Illumina array. The inflammation gene interleukin-18 (IL18) showed the strongest evidence for association with epithelial ovarian cancer in a gene-by-gene analysis (P = 0.002) with a &lt;25% chance of being a false-positive finding (q value = 0.240). Using a multivariate model search algorithm over 11 IL18 tagging SNPs, we found that the association was best modeled by rs1834481. Further, this SNP uniquely tagged a significantly associated IL18 haplotype and there was an increased risk of epithelial ovarian cancer per rs1834481 allele (odds ratio, 1.24; 95% confidence interval, 1.06-1.45). In a replication stage, 12 independent studies from the Ovarian Cancer Association Consortium (OCAC) genotyped rs1834481 in an additional 5,877 cases and 7,791 controls. The fixed effects estimate per rs1834481 allele was null (odds ratio, 0.99; 95% confidence interval, 0.94-1.05) when data from the 12 OCAC studies were combined. The effect estimate remained unchanged with the addition of the initial North Carolina Ovarian Cancer Study data. This analysis shows the importance of consortia, like the OCAC, in either confirming or refuting the validity of putative findings in studies with smaller sample sizes. (Cancer Epidemiol Biomarkers Prev 2008;17(12):3567–72)