Conti, David V.

Abstract Background Observational studies have shown higher folate consumption to be associated with lower risk of colorectal cancer (CRC). Understanding whether and how genetic risk factors interact with folate could further elucidate the underlying mechanism. Aggregating functionally relevant genetic variants in set-based variant testing has higher power to detect gene-environment (GxE) interactions and may provide information on the underlying biological pathway. Objective We investigated interactions between folate consumption and predicted gene expression on colorectal cancer risk across the genome. Methods We used variant weights from the PrediXcan models of colon tissue-specific gene expression as a priori variant information for a set-based GxE approach. We harmonized total folate intake (mcg/day) based on dietary intake and supplemental use across cohort and case-control studies and calculated sex and study specific quantiles. Analyses were performed using a mixed effects score tests for interactions between folate and genetically predicted expression of 4,839 genes with available genetically predicted expression. We meta-analyzed results across 23 studies for a total of 13,498 cases with colorectal tumors and 13,918 controls of European ancestry. Results We found suggestive evidence of interaction with folate intake for genes including glutathione S-Transferase Alpha 1 (GSTA1; p=4.3E-4), Tonsuko Like, DNA Repair Protein (TONSL; p=4.3E-4), and Aspartylglucosaminidase (AGA: p=4.5E-4). Glutathione is an antioxidant, preventing cellular damage and is a downstream metabolite of homocysteine and metabolized by GSTA1. TONSL is part of a complex that functions in the recovery of double strand breaks and AGA plays a role in lysosomal breakdown of glycoprotein. Conclusion We identified three genes involved in preventing or repairing DNA damage that may interact with folate consumption to alter CRC risk.

Objective An understanding of the etiologic heterogeneity of colorectal cancer (CRC) is critical for improving precision prevention, including individualized screening recommendations and the discovery of novel drug targets and repurposable drug candidates for chemoprevention. Known differences in molecular characteristics and environmental risk factors among tumors arising in different locations of the colorectum suggest partly distinct mechanisms of carcinogenesis. The extent to which the contribution of inherited genetic risk factors for CRC differs by anatomical subsite of the primary tumor has not been examined. Design To identify new anatomical subsite-specific risk loci, we performed genome-wide association study (GWAS) meta-analyses including data of 48 214 CRC cases and 64 159 controls of European ancestry. We characterised effect heterogeneity at CRC risk loci using multinomial modelling. Results We identified 13 loci that reached genome-wide significance (p<5×10 −8 ) and that were not reported by previous GWASs for overall CRC risk. Multiple lines of evidence support candidate genes at several of these loci. We detected substantial heterogeneity between anatomical subsites. Just over half (61) of 109 known and new risk variants showed no evidence for heterogeneity. In contrast, 22 variants showed association with distal CRC (including rectal cancer), but no evidence for association or an attenuated association with proximal CRC. For two loci, there was strong evidence for effects confined to proximal colon cancer. Conclusion Genetic architectures of proximal and distal CRC are partly distinct. Studies of risk factors and mechanisms of carcinogenesis, and precision prevention strategies should take into consideration the anatomical subsite of the tumour.

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: Persons of African ancestry (AA) experience a 1.5-2-fold risk of multiple myeloma (MM) compared to persons of European ancestry (EA). We assembled a set of MM patients with self-reported AA in order to evaluate the contribution of genetics to etiology in this high-risk group. Methods: Here we present the results of a meta-analysis of two GWAS in 1,813 cases and 8,871 controls of AA. We also conducted an admixture mapping scan to identify risk alleles associated with local ancestry, fine-mapped the 23 known susceptibility loci to find markers that could better capture MM risk in individuals of AA, and constructed a polygenic risk score (PRS) to assess the aggregated effect of known MM risk alleles. Finally, we conducted an eQTL analysis measuring gene expression in those genes harboring a risk variant in malignant plasma cells from 292 of the patients from a single site. Results: In GWAS analysis, we identified two suggestive novel loci located at 9p24.3 and 9p13.1 at P&lt;1 × 10-6, but no genome-wide significant association was noted. In admixture mapping, we observed a genome-wide significant inverse association between local AA at 2p24.1-23.1 and MM risk in AA individuals. 20 of the 23 known EA risk variants showed directional consistency and 9 replicated at P&lt;0.05 in AA individuals. In eight regions, we identified markers that better capture MM risk in persons of AA. AA individuals with a PRS in the top 10% had a 1.82-fold (95%CI: 1.56, 2.11) increased MM risk compared to those with average risk (25-75%). The strongest functional association was between the risk allele for variant rs56219066 at 5q15 and lower ELL2 expression (P= 5.1 × 10–12). Conclusion: Our study shows that common genetic variation contributes to MM risk individuals of AA. This abstract is also being presented as Poster C040. Citation Format: Zhaohui Du, Niels Weinhold, Gregory Chi Song, Kristen A. Rand, David J. Van Den Berg, Amie E. Hwang, Xin Sheng, Victor Hom, Sikander Ailawadhi, Ajay K. Nooka, Seema Singhal, Karen Pawlish, Edward Peters, Cathryn Bock, Ann Mohrbacher, Alexander Stram, Sonja I. Berndt, William J. Blot, Graham Casey, Victoria L. Stevens, Rick Kittles, Phyllis J. Goodman, W. Ryan Diver, Anselm Hennis, Barbara Nemesure, Eric A. Klein, Benjamin A. Rybicki, Janet L. Stanford, John S. Witte, Lisa Signorello, Esther M. John, Leslie Bernstein, Antoinette Stroup, Owen W. Stephens, Maurizio Zangari, Frits Van Rhee, Andrew Olshan, Wei Zheng, Jennifer J. Hu, Regina Ziegler, Sarah J. Nyante, Sue Ann Ingles, Michael Press, John David Carpten, Stephen Chanock, Jayesh Mehta, Graham A Colditz, Jeffrey Wolf, Thomas G. Martin, Michael Tomasson, Mark A. Fiala, Howard Terebelo, Nalini Janakiraman, Laurence Kolonel, Kenneth C. Anderson, Loic Le Marchand, Daniel Auclair, Brian C.-H. Chiu, Elad Ziv, Daniel Stram, Ravi Vij, Leon Bernal-Mizrachi, Gareth J. Morgan, Jeffrey A. Zonder, Carol Ann Huff, Sagar Lonial, Robert Z. Orlowski, David V. Conti, Christopher A. Haiman, Wendy Cozen. A meta-analysis of genome-wide association study and eQTL analysis of multiple myeloma among African Americans [abstract]. In: Proceedings of the Twelfth AACR Conference on the Science of Cancer Health Disparities in Racial/Ethnic Minorities and the Medically Underserved; 2019 Sep 20-23; San Francisco, CA. Philadelphia (PA): AACR; Cancer Epidemiol Biomarkers Prev 2020;29(6 Suppl_2):Abstract nr PR05.

Purpose: Multiple myeloma (MM) treatment has changed tremendously, with significant improvement in patient out-comes. One group with a suboptimal benefit is patients with high-risk cytogenetics, as tested by conventional karyotyping or fluorescence in situ hybridization (FISH). Methodology for these tests has been published, but not necessarily standardized. Methods: We address variability in the testing and reporting methodology for MM cytogenetics in the United States using the ongoing African American Multiple Myeloma Study (AAMMS). We evaluated clinical and cytogenetic data from 1,221 patients (1,161 with conventional karyotyping and 976 with FISH) tested between 1998 and 2016 across 58 laboratories nationwide. Results: Interlab and intralab variability was noted for the number of cells analyzed for karyotyping, with a significantly higher number of cells analyzed in patients in whom cytogenetics were normal (P 5.0025). For FISH testing, CD138-positive cell enrichment was used in 29.7% of patients and no enrichment in 50% of patients, whereas the remainder had unknown status. A significantly smaller number of cells was analyzed for patients in which CD138 cell enrichment was used compared with those without such enrichment (median, 50 v 200; P, .0001). A median of 7 loci probes (range, 1-16) were used for FISH testing across all laboratories, with variability in the loci probed even within a given laboratory. Chromosome 13–related abnormalities were the most frequently tested abnormality (n5956; 97.9%), and t(14;16) was the least frequently tested abnormality (n 5 119; 12.2%). Conclusions: We report significant variability in cytogenetic testing across the United States for MM, potentially leading to variability in risk stratification, with possible clinical implications and personalized treatment approaches.

AbstractBackgroundBody mass index (BMI) and diabetes are established risk factors for colorectal cancer (CRC), likely through perturbations in metabolic traits (e.g. insulin resistance and glucose homeostasis). Identification of interactions between variation in genes and these metabolic risk factors may identify novel biologic insights into CRC etiology.MethodsTo improve statistical power and interpretation for gene‐environment interaction (G × E) testing, we tested genetic variants that regulate expression of a gene together for interaction with BMI (kg/m2) and diabetes on CRC risk among 26 017 cases and 20 692 controls. Each variant was weighted based on PrediXcan analysis of gene expression data from colon tissue generated in the Genotype‐Tissue Expression Project for all genes with heritability ≥1%. We used a mixed‐effects model to jointly measure the G × E interaction in a gene by partitioning the interactions into the predicted gene expression levels (fixed effects), and residual G × E effects (random effects). G × BMI analyses were stratified by sex as BMI‐CRC associations differ by sex. We used false discovery rates to account for multiple comparisons and reported all results with FDR <0.2.ResultsAmong 4839 genes tested, genetically predicted expressions of FOXA1 (P = 3.15 × 10−5), PSMC5 (P = 4.51 × 10−4) and CD33 (P = 2.71 × 10−4) modified the association of BMI on CRC risk for men; KIAA0753 (P = 2.29 × 10−5) and SCN1B (P = 2.76 × 10−4) modified the association of BMI on CRC risk for women; and PTPN2 modified the association between diabetes and CRC risk in both sexes (P = 2.31 × 10−5).ConclusionsAggregating G × E interactions and incorporating functional information, we discovered novel genes that may interact with BMI and diabetes on CRC risk.

Abstract Background: Results from epidemiologic studies examining polyunsaturated fatty acids (PUFA) and colorectal cancer risk are inconsistent. Mendelian randomization may strengthen causal inference from observational studies. Given their shared metabolic pathway, examining the combined effects of aspirin/NSAID use with PUFAs could help elucidate an association between PUFAs and colorectal cancer risk. Methods: Information was leveraged from genome-wide association studies (GWAS) regarding PUFA-associated SNPs to create weighted genetic scores (wGS) representing genetically predicted circulating blood PUFAs for 11,016 non-Hispanic white colorectal cancer cases and 13,732 controls in the Genetics and Epidemiology of Colorectal Cancer Consortium (GECCO). Associations per SD increase in the wGS were estimated using unconditional logistic regression. Interactions between PUFA wGSs and aspirin/NSAID use on colorectal cancer risk were also examined. Results: Modest colorectal cancer risk reductions were observed per SD increase in circulating linoleic acid [ORLA = 0.96; 95% confidence interval (CI) = 0.93–0.98; P = 5.2 × 10−4] and α-linolenic acid (ORALA = 0.95; 95% CI = 0.92–0.97; P = 5.4 × 10−5), whereas modest increased risks were observed for arachidonic (ORAA = 1.06; 95% CI = 1.03–1.08; P = 3.3 × 10−5), eicosapentaenoic (OREPA = 1.04; 95% CI = 1.01–1.07; P = 2.5 × 10−3), and docosapentaenoic acids (ORDPA = 1.03; 95% CI = 1.01–1.06; P = 1.2 × 10−2). Each of these effects was stronger among aspirin/NSAID nonusers in the stratified analyses. Conclusions: Our study suggests that higher circulating shorter-chain PUFAs (i.e., LA and ALA) were associated with reduced colorectal cancer risk, whereas longer-chain PUFAs (i.e., AA, EPA, and DPA) were associated with an increased colorectal cancer risk. Impact: The interaction of PUFAs with aspirin/NSAID use indicates a shared colorectal cancer inflammatory pathway. Future research should continue to improve PUFA genetic instruments to elucidate the independent effects of PUFAs on colorectal cancer.

PURPOSE In studies of men of European ancestry, rare pathogenic variants in DNA repair pathway genes have been shown to be associated with risk of aggressive prostate cancer. The contribution of rare coding variation to prostate cancer risk in men of African ancestry has not been established. METHODS We sequenced a panel of 19 DNA repair and cancer predisposition genes in 2,453 African American and 1,151 Ugandan cases and controls with prostate cancer. Rare variants were classified as pathogenic or putatively functionally disruptive and examined in association with prostate cancer risk and disease aggressiveness in gene and pathway-level association analyses. RESULTS Pathogenic variants were found in 75 of 2,098 cases (3.6%) and 31 of 1,481 controls (2.1%; odds ratio [OR], 1.82; 95% CI, 1.19 to 2.79; P = .0044), with the association being stronger for more aggressive disease phenotypes (OR, 3.10; 95% CI, 1.54 to 6.23; P = .0022). The highest risks for aggressive disease were observed with pathogenic variants in the ATM, BRCA2, PALB2, and NBN genes, with ORs ranging from approximately 4 to 15 in the combined study sample of African American and Ugandan men. Rare, nonpathogenic, nonsynonymous variants did not have a major impact on risk of overall prostate cancer or disease aggressiveness. CONCLUSION Rare pathogenic variants in DNA repair genes have appreciable effects on risk of aggressive prostate cancer in men of African ancestry. These findings have potential implications for panel testing and risk stratification in this high-risk population.

Abstract Persons of African ancestry (AA) have a twofold higher risk for multiple myeloma (MM) compared with persons of European ancestry (EA). Genome-wide association studies (GWASs) support a genetic contribution to MM etiology in individuals of EA. Little is known about genetic risk factors for MM in individuals of AA. We performed a meta-analysis of 2 GWASs of MM in 1813 cases and 8871 controls and conducted an admixture mapping scan to identify risk alleles. We fine-mapped the 23 known susceptibility loci to find markers that could better capture MM risk in individuals of AA and constructed a polygenic risk score (PRS) to assess the aggregated effect of known MM risk alleles. In GWAS meta-analysis, we identified 2 suggestive novel loci located at 9p24.3 and 9p13.1 at P &lt; 1 × 10−6; however, no genome-wide significant association was noted. In admixture mapping, we observed a genome-wide significant inverse association between local AA at 2p24.1-23.1 and MM risk in AA individuals. Of the 23 known EA risk variants, 20 showed directional consistency, and 9 replicated at P &lt; .05 in AA individuals. In 8 regions, we identified markers that better capture MM risk in persons with AA. AA individuals with a PRS in the top 10% had a 1.82-fold (95% confidence interval, 1.56-2.11) increased MM risk compared with those with average risk (25%-75%). The strongest functional association was between the risk allele for variant rs56219066 at 5q15 and lower ELL2 expression (P = 5.1 × 10−12). Our study shows that common genetic variation contributes to MM risk in individuals with AA.

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.