McKay, James

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.

Telomeres cap chromosome ends, protecting them from degradation, double‐strand breaks, and end‐to‐end fusions. Telomeres are maintained by telomerase, a reverse transcriptase encoded by TERT, and an RNA template encoded by TERC. Loci in the TERT and adjoining CLPTM1L region are associated with risk of multiple cancers. We therefore investigated associations between variants in 22 telomere structure and maintenance gene regions and colorectal, breast, prostate, ovarian, and lung cancer risk. We performed subset‐based meta‐analyses of 204,993 directly‐measured and imputed SNPs among 61,851 cancer cases and 74,457 controls of European descent. Independent associations for SNP minor alleles were identified using sequential conditional analysis (with gene‐level p value cutoffs ≤3.08 × 10−5). Of the thirteen independent SNPs observed to be associated with cancer risk, novel findings were observed for seven loci. Across the DCLRE1B region, rs974494 and rs12144215 were inversely associated with prostate and lung cancers, and colorectal, breast, and prostate cancers, respectively. Across the TERC region, rs75316749 was positively associated with colorectal, breast, ovarian, and lung cancers. Across the DCLRE1B region, rs974404 and rs12144215 were inversely associated with prostate and lung cancers, and colorectal, breast, and prostate cancers, respectively. Near POT1, rs116895242 was inversely associated with colorectal, ovarian, and lung cancers, and RTEL1 rs34978822 was inversely associated with prostate and lung cancers. The complex association patterns in telomere‐related genes across cancer types may provide insight into mechanisms through which telomere dysfunction in different tissues influences cancer risk.

Abstract Telomeres are complex structures that cap chromosome ends, protecting them from degradation, double strand breaks and end-to-end fusions. Telomeres are maintained by the enzyme telomerase, which is made up of a reverse transcriptase encoded by TERT, and an RNA template encoded by TERC. The telomere structure itself is composed of proteins encoded by ACD, ACYP2, BICD1, DKC1, DCLRE1B, MPHOSPH6, NAF1, NOLA1, NOLA2, NOLA3, OBFC1, PIK3C3, POT1, RTEL1, TEP1, TERF1, TERF2, TERF21P, TINF2, TNKS/PINX1, TPP1 and ZNF208. Several single nucleotide polymorphisms (SNPs) in the TERT and adjoining CLPTM1L gene region are associated with multiple cancer types, and some are associated with both increased and decreased risks across different cancer types. We systematically characterized the patterns of association between variants in these 25 telomere structure and maintenance genes and risk across five cancer types in the Genetic Association and Mechanisms in Oncology (GAME-ON) consortium. We performed a subset-based meta-analysis (ASSET) of 209,367 directly measured and imputed SNPs, one megabase up- and downstream of these genes, across genome-wide association studies of colorectal (5,100 cases, 4,831 controls), lung (12,160 cases, 16,838 controls), breast (15,748 cases, 18,084 controls), ovarian (4,369 cases, 9,123 controls) and prostate (14,160 cases, 12,724 controls) cancers. Correlations (r2) between SNPs were examined in Haploview using the 1000 Genomes Project CEU population. A total of 87 TERT, 123 TERC and 26 DCLRE1B SNPs were associated with cancer risk at gene-level Bonferroni-corrected p-values of 4.2-7.8×10−6. Patterns of association were similar for prostate and colorectal cancers in DCLRE1B. Much stronger associations were observed in TERT and TERC, with 63 and 24 SNPs reaching genome-wide significance (p<5.0×10−8), respectively. Of these, 9 TERT and 9 TERC SNPs were correlated at r2<0.75. The most strongly associated SNPs in TERT (rs37004, p = 2.6×10−11) and TERC (rs76925190, p = 1.5×10−15) are in regions containing documented risk loci for lung and prostate cancers, respectively. TERT rs37004 (p = 1.2×10−13) and 6 of the other 9 SNPs in TERT were associated only with lung cancer risk. For the other 2 SNPs, there was a suggestion that patterns of risk were opposite for lung and prostate cancers. While TERC rs76925190 was strongly associated with prostate cancer risk (p = 5.4 x10−17), it was also suggestively associated with colorectal cancer risk (p<0.005). This pattern was similar for another 2 of the 9 SNPs, while for the rest, associations were limited to prostate cancer. The complex patterns of association in telomere structure and maintenance genes observed across cancer types may provide insight about the mechanisms through which telomere dysfunction in different tissues influences cancer risk. Citation Format: Sara Karami, Younghun Han, Fredrick R. Schumacher, Zsofia Kote-Jarai, Sara Lindstrom, John S. Witte, Iona Cheng, Shenying Fang, Jiali Han, Peter Kraft, Fengju Song, Rayjean J. Hung, James McKay, Stephen J. Chanock, Mala Pande, Angela Risch, Hongbing Shen, Christopher A. Haiman, Lisa Boardman, Cornelia M. Ulrich, Graham Casey, Ulrike Peters, Nilanjan Chatterjee, Brandon Pierce, Wei Zheng, Christopher I. Amos, Jennifer A. Doherty. Risk loci in telomere structure and maintenance genes across five cancer types: GAME-ON Consortium. [abstract]. In: Proceedings of the 106th Annual Meeting of the American Association for Cancer Research; 2015 Apr 18-22; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Res 2015;75(15 Suppl):Abstract nr 4609. doi:10.1158/1538-7445.AM2015-4609

Abstract Background: Diffuse large B-cell lymphoma (DLBCL) is the most common lymphoma subtype and is clinically aggressive. Relatives of DLBCL patients are at elevated risk for DLBCL. In small genome-wide association studies (GWAS) in which all NHL subtypes were combined, no conclusive loci for DLBCL have been previously identified in individuals of European background, whereas a recent study conducted in East Asia identified a locus at 3q27. Methods: To identify genetic susceptibility loci for DLBCL, we conducted a meta-analysis of three new GWAS and one prior scan, totaling 3,857 DLBCL cases and 7,666 controls of European ancestry. Imputation was used to combine all the data using the 1000 Genomes Project release v3 and IMPUTE2. Based on these findings, we selected 6 promising SNPs for genotyping in an additional 202 cases and 3,431 controls, which were then included in a second meta-analysis. SNPTEST was used to estimate odds ratios (ORs), using the additive model and adjusting for age, gender and significant eigenvectors. Results: In our multi-stage analysis, five independent SNPs in four loci achieved genome-wide significance marked by rs116446171 at 6p25.3 (EXOC2; OR=2.33, P=1.60x10-20), rs13255292 and rs4733601 at 8q24.21 (PVT1; OR=1.23, P=4.86x10-12 and OR=1.18, P=1.63x10-8, respectively), rs2523607 at 6p21.33 (HLA-B; OR=1.34, 3.35x10-9) and rs2681416 at 3q13.33 (CD86; OR=1.20, P=6.04x10-9). The two 8q24.21 SNPs displayed minimal linkage disequilibrium (r2=0.03) and conditional analysis supported their independence. We did not observe a notable signal for a locus previously reported for DLBCL on 3q27 in East Asia, rs6773854 (OR=1.06, P=0.81). We estimated that common SNPs, including but not limited to the loci discovered in this study, explain approximately 16% of the variance for DLBCL risk overall. Discussion: The susceptibility locus at 6p25.3 maps near a plausible DLBCL candidate gene, EXOC2 (exocyst complex component 2), which interacts with Ral proteins, and disruption of the Ral-exocyst regulatory node impacts cellular and developmental processes associated with malignant transformation and progression. The two 8q24.21 variants are approximately 1Mb telomeric to the 8q24 region linked with multiple cancers including CLL. Both variants are positioned in close proximity to PVT1, which is a non-coding RNA implicated in the MYC activation, and MYC is known to be deregulated in some DLBCLs. The locus at 3q13.33 maps to CD86 (B7-2), a B-7 family member that along with CD80 (B7-1) is essential in cognate interactions between T- and B-lymphocytes; both play key roles in immunodeficiency, autoimmune diseases, and in anti-tumor immunity. Finally, we had multiple signals in the HLA region, the strongest of which was at 6p21.33 near HLA-B. HLA-B encodes the HLA class I heavy chain paralogue, which heterodimerizes with a light chain (β2 microglobulin) and plays a central role in presenting intracellularly processed self or foreign antigens to CD8+ cytotoxic T lymphocytes. Class I molecules have been extensively linked to the risk of a variety of immune-mediated diseases and cancers including DLBCL. Conclusions: These data provide substantial new evidence for genetic susceptibility to this B-cell malignancy, and point towards pathways involved in immune recognition and immune function in the pathogenesis of DLBCL. Citation Format: James R. Cerhan, Sonja I. Berndt, Joseph Vijai, Hervé Ghesquières, James McKay, Sophia S. Wang, Zhaoming Wang, Meredith Yeager, Alexandra Nieters, David Cox, Alain Monnereau, Christopher Flowers, Anneclaire J. De Roos, Angela R. Brooks-Wilson, Qing Lan, Gianluca Severi, Mads Melbye, Rebecca D. Jackson, Lauren R. Teras, Mark Purdue, Claire Vajdic, Demetrius Albanes, Kimberly A. Bertrand, Anne Zeleniuch-Jacquotte, Simon Crouch, Yawei Zhang, Paolo Vineis, Susan L. Slager, Karin E. Smedby, Gilles Salles, Christine F. Skibola, Nathaniel Rothman, Stephen J. Chanock, on behalf of the NHL GWAS Project. Genome-wide association study identifies multiple susceptibility loci for diffuse large B-cell lymphoma. [abstract]. In: Proceedings of the 105th Annual Meeting of the American Association for Cancer Research; 2014 Apr 5-9; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2014;74(19 Suppl):Abstract nr LB-272. doi:10.1158/1538-7445.AM2014-LB-272

Abstract Background. Follicular lymphoma (FL) is an indolent B-cell malignancy with a variable clinical course that may transform to more aggressive forms of lymphoma. Previous modestly powered GWAS have reported multiple FL susceptibility loci in the HLA class I and II regions at 6p21.32-33, but novel non-HLA susceptibility loci and the genetic architecture of FL associations in the HLA region remain to be elucidated. Methods. We conducted the largest GWAS of FL to date consisting of 2,142 cases and 6,221 controls of European ancestry from 22 studies scanned on the Illumina OmniExpress. After imputation of common SNPs using IMPUTE2 and 1,000 Genomes Project v3 data, we conducted a meta-analysis of these data with two previous FL GWAS (n=586 cases, 1,537 controls) and replicated top hits in an additional 629 cases and 4,283 controls. Expression quantitative trait loci (eQTL) analyses were performed to evaluate the effects of associated variants on gene expression. Classical HLA allele imputations and stepwise regression analyses are underway to further characterize the HLA associations. Results. In the joint meta-analysis, the strongest association with FL risk was observed in the HLA region, where a large number of SNPs showed genome-wide significance (P<5x10-8 to P= 1.84x10-84), particularly in HLA class II at 6p21.32 where several HLA eQTLs were identified (P < 0.05). Outside HLA, we identified four novel loci associated with FL at 11q23.3 (CXCR5, rs4938573; P=3.14x10-15), 11q24.3 (ETS1, rs4937362; P=3.33x10-9), 3q13.33 (CD86, rs2681416; P=9.79x10-11) and 3q28 (LPP, rs6444305, P=2.55x10-8). These novel susceptibility loci were located in four biologically relevant genes involved in immune regulation: CXCR5 and its ligand, CXCL13, are involved in guiding B-cells into the B-cell zones of secondary lymphoid organs as well as T-cell migration. ETS1 is expressed in lymphoid cells and regulates immune cell function including differentiation, survival, and proliferation. CD86 and CD80 are classic members of the B7 costimulatory pathway that is important in maintaining immune function, suppression of autoimmunity and antitumor surveillance. LPP encodes a LIM domain-containing protein of the zyxin family and participates in cell adhesion, cell migration, proliferation and transcription dynamics. Conclusions. This large GWAS provides further support for the important role of common genetic variation in non-HLA genes and further evidence of the key role that HLA immune-regulatory genes play in the pathogenesis of FL. Citation Format: Christine F. Skibola, Sonja I. Berndt, James R. Cerhan, Zhaoming Wang, Joseph Vijai, Lucia Conde, Paul de Bakker, Sophia S. Wang, Claire M. Vajdic, Brenda M. Birmann, Susan L. Slager, James McKay, Paige M. Bracci, Alexandra Nieters, Qing Lan, Angela R. Brooks-Wilson, Martha S. Linet, Demetrius Albanes, John J. Spinelli, Roel C.H. Vermeulen, Mark P. Purdue, Meredith Yaeger, Lauren R. Teras, Silvia de Sanjose, Alain Monnereau, Simon Crouch, Jia Nee Foo, Henrik Hjalgrim, Gianluca Severi, Brian K. Link, Kimberly A. Bertrand, Yawei Zhang, Karin E. Smedby, Stephen J. Chanock, Nathaniel Rothman, NHL GWAS Consortium. Meta-analysis of genome-wide association studies identifies novel susceptibility loci for follicular lymphoma. [abstract]. In: Proceedings of the 105th Annual Meeting of the American Association for Cancer Research; 2014 Apr 5-9; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2014;74(19 Suppl):Abstract nr 5072. doi:10.1158/1538-7445.AM2014-5072

Abstract Background: The single-nucleotide polymorphism (SNP) 5p12-rs10941679 has been found to be associated with risk of breast cancer, particularly estrogen receptor (ER)-positive disease. We aimed to further explore this association overall, and by tumor histopathology, in the Breast Cancer Association Consortium. Methods: Data were combined from 37 studies, including 40,972 invasive cases, 1,398 cases of ductal carcinoma in situ (DCIS), and 46,334 controls, all of white European ancestry, as well as 3,007 invasive cases and 2,337 controls of Asian ancestry. Associations overall and by tumor invasiveness and histopathology were assessed using logistic regression. Results: For white Europeans, the per-allele OR associated with 5p12-rs10941679 was 1.11 (95% CI = 1.08–1.14, P = 7 × 10−18) for invasive breast cancer and 1.10 (95% CI = 1.01–1.21, P = 0.03) for DCIS. For Asian women, the estimated OR for invasive disease was similar (OR = 1.07, 95%CI = 0.99–1.15, P = 0.09). Further analyses suggested that the association in white Europeans was largely limited to progesterone receptor (PR)-positive disease (per-allele OR = 1.16, 95% CI = 1.12–1.20, P = 1 × 10−18 vs. OR = 1.03, 95% CI = 0.99–1.07, P = 0.2 for PR-negative disease; Pheterogeneity = 2 × 10−7); heterogeneity by ER status was not observed (P = 0.2) once PR status was accounted for. The association was also stronger for lower grade tumors [per-allele OR (95% CI) = 1.20 (1.14–1.25), 1.13 (1.09–1.16), and 1.04 (0.99–1.08) for grade 1, 2, and 3/4, respectively; Ptrend = 5 × 10−7]. Conclusion: 5p12 is a breast cancer susceptibility locus for PR-positive, lower grade breast cancer. Impact: Multicenter fine-mapping studies of this region are needed as a first step to identifying the causal variant or variants. Cancer Epidemiol Biomarkers Prev; 20(10); 2222–31. ©2011 AACR.