Stanford, Janet L.

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<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<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.

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 < 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 < .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 Multiple myeloma (MM) is 2-3 times more common among African-Americans compared to non-Hispanic whites. The 2-3-fold increased risk among family members of cases suggests a genetic contribution to risk. Genome-wide association studies (GWAS) in populations of European ancestry have identified seven novel risk loci at 2p23.3 (rs6746082), 3q26.2 (rs10936599), 3p22.1 (rs1052501), 6p21.32 (rs2285803), 7p15.3 (rs4487645), 17p11.2 (rs4273077) and 22q13.1 (rs877529) (Broderick, et al. Nat Genet, 2011, Chubb, et al. Nat Genet, 2013), three of which were replicated in another European series (Martino et al., Br J Haematol, 2012). Here we examined the index signals and conducted fine-mapping for each locus in a case-control study of 1,049 multiple myeloma cases and 7,084 controls of African ancestry to identify better markers of risk and novel independent loci in seven previously reported regions in this high risk population. Incident cases were recruited from 10 clinical centers and SEER cancer registries from 2011 to 2013 and genotyped using the Illumina HumanCore GWAS array. Control data were obtained from previous genome-wide studies of breast and prostate cancer, genotyped using the Illumina 1M-Duo in 4425 male controls from the African Ancestry Prostate Cancer Consortium (consisting of 14 independent studies) and 2632 female controls from a breast cancer GWAS of African-American women (consisting of 9 independent studies). Imputation to 1000 Genomes (March 2012 release) was conducted for regions around six of the previously identified single nucleotide polymorphisms [SNPs] (the HLA region harboring rs2285803 is still being imputed, results will be presented). A case-control analysis of SNPs/indels >1% frequency within 250 kb of each index variant was conducted using unconditional multivariable logistic regression adjusting for age, sex and five leading principal components. Region-specific alpha levels were determined through permutation tests. The minimum alpha level across the six regions was α=0.002. All previously reported risk variants were common in African-Americans (minor allele frequency [MAF] >0.05). For five of the six SNPs, we had ≥94% power to detect the same effect observed in non-Hispanic whites, and 64% power for the less common variant rs10936599 (MAF=0.07). We observed directionally consistent effects (odds ratio [OR]>1) for the six risk variants tested, with three replicating at p≤0.05 (7p15.3, p=1.4x10-7; 17p11.2, p=0.05; 22q13.1, p=0.02). For three of the six regions, we observed better markers of risk in African-Americans that were correlated with the index SNP in Europeans (7p15.3, rs56333627, p=1.5x10-5, r2=0.89; 17p11.2, rs34562254, p=2.9x10-3, r2=0.90; 22q13.1, rs2092410, p=1.1x10-4 r2=.71). The missense variant identified in the 17p11.2 region (rs34562254, Pro251Leu) is located in TNFRSF13B, which encodes the protein TACI, a B cell surface receptor which plays a role in B cell maturation, apoptosis and antibody production by inducing activation of transcription factors including NFAT and NFκβ. In addition, there is evidence suggesting that TACI is involved in MM pathogenesis. Our results demonstrate that many of the risk loci for MM found in European ancestry populations are also risk loci in men and women of African ancestry and that by fine-mapping, we are able to identify variants that better capture risk in populations of African ancestry. Disclosures Terebelo: Celgene Corp: Membership on an entity's Board of Directors or advisory committees. Lonial:Millennium: The Takeda Oncology Company: Consultancy, Research Funding; Celgene: Consultancy, Research Funding; Novartis: Consultancy, Research Funding; Bristol-Myers Squibb: Consultancy, Research Funding; Onyx Pharmaceuticals: Consultancy, Research Funding.

Abstract Background: Whether associations between circulating metabolites and prostate cancer are causal is unknown. We report on the largest study of metabolites and prostate cancer (2,291 cases and 2,661 controls) and appraise causality for a subset of the prostate cancer–metabolite associations using two-sample Mendelian randomization (MR). Methods: The case–control portion of the study was conducted in nine UK centers with men ages 50–69 years who underwent prostate-specific antigen screening for prostate cancer within the Prostate Testing for Cancer and Treatment (ProtecT) trial. Two data sources were used to appraise causality: a genome-wide association study (GWAS) of metabolites in 24,925 participants and a GWAS of prostate cancer in 44,825 cases and 27,904 controls within the Association Group to Investigate Cancer Associated Alterations in the Genome (PRACTICAL) consortium. Results: Thirty-five metabolites were strongly associated with prostate cancer (P < 0.0014, multiple-testing threshold). These fell into four classes: (i) lipids and lipoprotein subclass characteristics (total cholesterol and ratios, cholesterol esters and ratios, free cholesterol and ratios, phospholipids and ratios, and triglyceride ratios); (ii) fatty acids and ratios; (iii) amino acids; (iv) and fluid balance. Fourteen top metabolites were proxied by genetic variables, but MR indicated these were not causal. Conclusions: We identified 35 circulating metabolites associated with prostate cancer presence, but found no evidence of causality for those 14 testable with MR. Thus, the 14 MR-tested metabolites are unlikely to be mechanistically important in prostate cancer risk. Impact: The metabolome provides a promising set of biomarkers that may aid prostate cancer classification.

AbstractMotivationThe use of single nucleotide polymorphism (SNP) interactions to predict complex diseases is getting more attention during the past decade, but related statistical methods are still immature. We previously proposed the SNP Interaction Pattern Identifier (SIPI) approach to evaluate 45 SNP interaction patterns/patterns. SIPI is statistically powerful but suffers from a large computation burden. For large-scale studies, it is necessary to use a powerful and computation-efficient method. The objective of this study is to develop an evidence-based mini-version of SIPI as the screening tool or solitary use and to evaluate the impact of inheritance mode and model structure on detecting SNP–SNP interactions.ResultsWe tested two candidate approaches: the ‘Five-Full’ and ‘AA9int’ method. The Five-Full approach is composed of the five full interaction models considering three inheritance modes (additive, dominant and recessive). The AA9int approach is composed of nine interaction models by considering non-hierarchical model structure and the additive mode. Our simulation results show that AA9int has similar statistical power compared to SIPI and is superior to the Five-Full approach, and the impact of the non-hierarchical model structure is greater than that of the inheritance mode in detecting SNP–SNP interactions. In summary, it is recommended that AA9int is a powerful tool to be used either alone or as the screening stage of a two-stage approach (AA9int+SIPI) for detecting SNP–SNP interactions in large-scale studies.Availability and implementationThe ‘AA9int’ and ‘parAA9int’ functions (standard and parallel computing version) are added in the SIPI R package, which is freely available at https://linhuiyi.github.io/LinHY_Software/.Supplementary informationSupplementary data are available at Bioinformatics online.

ABSTRACTNext‐generation sequencing technologies have afforded unprecedented characterization of low‐frequency and rare genetic variation. Due to low power for single‐variant testing, aggregative methods are commonly used to combine observed rare variation within a single gene. Causal variation may also aggregate across multiple genes within relevant biomolecular pathways. Kernel‐machine regression and adaptive testing methods for aggregative rare‐variant association testing have been demonstrated to be powerful approaches for pathway‐level analysis, although these methods tend to be computationally intensive at high‐variant dimensionality and require access to complete data. An additional analytical issue in scans of large pathway definition sets is multiple testing correction. Gene set definitions may exhibit substantial genic overlap, and the impact of the resultant correlation in test statistics on Type I error rate control for large agnostic gene set scans has not been fully explored. Herein, we first outline a statistical strategy for aggregative rare‐variant analysis using component gene‐level linear kernel score test summary statistics as well as derive simple estimators of the effective number of tests for family‐wise error rate control. We then conduct extensive simulation studies to characterize the behavior of our approach relative to direct application of kernel and adaptive methods under a variety of conditions. We also apply our method to two case‐control studies, respectively, evaluating rare variation in hereditary prostate cancer and schizophrenia. Finally, we provide open‐source R code for public use to facilitate easy application of our methods to existing rare‐variant analysis results.

Abstract Motivation Testing SNP–SNP interactions is considered as a key for overcoming bottlenecks of genetic association studies. However, related statistical methods for testing SNP–SNP interactions are underdeveloped. Results We propose the SNP Interaction Pattern Identifier (SIPI), which tests 45 biologically meaningful interaction patterns for a binary outcome. SIPI takes non-hierarchical models, inheritance modes and mode coding direction into consideration. The simulation results show that SIPI has higher power than MDR (Multifactor Dimensionality Reduction), AA_Full, Geno_Full (full interaction model with additive or genotypic mode) and SNPassoc in detecting interactions. Applying SIPI to the prostate cancer PRACTICAL consortium data with approximately 21 000 patients, the four SNP pairs in EGFR-EGFR, EGFR-MMP16 and EGFR-CSF1 were found to be associated with prostate cancer aggressiveness with the exact or similar pattern in the discovery and validation sets. A similar match for external validation of SNP–SNP interaction studies is suggested. We demonstrated that SIPI not only searches for more meaningful interaction patterns but can also overcome the unstable nature of interaction patterns. Availability and Implementation The SIPI software is freely available at http://publichealth.lsuhsc.edu/LinSoftware/. Supplementary information Supplementary data are available at Bioinformatics online.

Coffee consumption has been shown in some studies to be associated with lower risk of prostate cancer. However, it is unclear if this association is causal or due to confounding or reverse causality. We conducted a Mendelian randomisation analysis to investigate the causal effects of coffee consumption on prostate cancer risk and progression. We used two genetic variants robustly associated with caffeine intake (rs4410790 and rs2472297) as proxies for coffee consumption in a sample of 46,687 men of European ancestry from 25 studies in the PRACTICAL consortium. Associations between genetic variants and prostate cancer case status, stage and grade were assessed by logistic regression and with all‐cause and prostate cancer‐specific mortality using Cox proportional hazards regression. There was no clear evidence that a genetic risk score combining rs4410790 and rs2472297 was associated with prostate cancer risk (OR per additional coffee increasing allele: 1.01, 95% CI: 0.98,1.03) or having high‐grade compared to low‐grade disease (OR: 1.01, 95% CI: 0.97,1.04). There was some evidence that the genetic risk score was associated with higher odds of having nonlocalised compared to localised stage disease (OR: 1.03, 95% CI: 1.01, 1.06). Amongst men with prostate cancer, there was no clear association between the genetic risk score and all‐cause mortality (HR: 1.00, 95% CI: 0.97,1.04) or prostate cancer‐specific mortality (HR: 1.03, 95% CI: 0.98,1.08). These results, which should have less bias from confounding than observational estimates, are not consistent with a substantial effect of coffee consumption on reducing prostate cancer incidence or progression.

Prostate cancer is the most common cancer in men in developed countries, and is a target for risk reduction strategies. The effects of alcohol consumption on prostate cancer incidence and survival remain unclear, potentially due to methodological limitations of observational studies. In this study, we investigated the associations of genetic variants in alcohol‐metabolising genes with prostate cancer incidence and survival. We analysed data from 23,868 men with prostate cancer and 23,051 controls from 25 studies within the international PRACTICAL Consortium. Study‐specific associations of 68 single nucleotide polymorphisms (SNPs) in 8 alcohol‐metabolising genes (Alcohol Dehydrogenases (ADHs) and Aldehyde Dehydrogenases (ALDHs)) with prostate cancer diagnosis and prostate cancer‐specific mortality, by grade, were assessed using logistic and Cox regression models, respectively. The data across the 25 studies were meta‐analysed using fixed‐effect and random‐effects models. We found little evidence that variants in alcohol metabolising genes were associated with prostate cancer diagnosis. Four variants in two genes exceeded the multiple testing threshold for associations with prostate cancer mortality in fixed‐effect meta‐analyses. SNPs within ALDH1A2 associated with prostate cancer mortality were rs1441817 (fixed effects hazard ratio, HRfixed = 0.78; 95% confidence interval (95%CI):0.66,0.91; p values = 0.002); rs12910509, HRfixed = 0.76; 95%CI:0.64,0.91; p values = 0.003); and rs8041922 (HRfixed = 0.76; 95%CI:0.64,0.91; p values = 0.002). These SNPs were in linkage disequilibrium with each other. In ALDH1B1, rs10973794 (HRfixed = 1.43; 95%CI:1.14,1.79; p values = 0.002) was associated with prostate cancer mortality in men with low‐grade prostate cancer. These results suggest that alcohol consumption is unlikely to affect prostate cancer incidence, but it may influence disease progression.

Background The rarity of mutations in PALB2, CHEK2 and ATM make it difficult to estimate precisely associated cancer risks. Population-based family studies have provided evidence that at least some of these mutations are associated with breast cancer risk as high as those associated with rare BRCA2 mutations. We aimed to estimate the relative risks associated with specific rare variants in PALB2, CHEK2 and ATM via a multicentre case-control study. Methods We genotyped 10 rare mutations using the custom iCOGS array: PALB2 c.1592delT, c.2816T>G and c.3113G>A, CHEK2 c.349A>G, c.538C>T, c.715G>A, c.1036C>T, c.1312G>T, and c.1343T>G and ATM c.7271T>G. We assessed associations with breast cancer risk (42 671 cases and 42 164 controls), as well as prostate (22 301 cases and 22 320 controls) and ovarian (14 542 cases and 23 491 controls) cancer risk, for each variant. Results For European women, strong evidence of association with breast cancer risk was observed for PALB2 c.1592delT OR 3.44 (95% CI 1.39 to 8.52, p=7.1×10−5), PALB2 c.3113G>A OR 4.21 (95% CI 1.84 to 9.60, p=6.9×10−8) and ATM c.7271T>G OR 11.0 (95% CI 1.42 to 85.7, p=0.0012). We also found evidence of association with breast cancer risk for three variants in CHEK2, c.349A>G OR 2.26 (95% CI 1.29 to 3.95), c.1036C>T OR 5.06 (95% CI 1.09 to 23.5) and c.538C>T OR 1.33 (95% CI 1.05 to 1.67) (p≤0.017). Evidence for prostate cancer risk was observed for CHEK2 c.1343T>G OR 3.03 (95% CI 1.53 to 6.03, p=0.0006) for African men and CHEK2 c.1312G>T OR 2.21 (95% CI 1.06 to 4.63, p=0.030) for European men. No evidence of association with ovarian cancer was found for any of these variants. Conclusions This report adds to accumulating evidence that at least some variants in these genes are associated with an increased risk of breast cancer that is clinically important.