Stevens, Victoria L.

Abstract Background Lung adenocarcinoma in never-smokers is a major public health burden, especially among East Asian women. Polygenic risk scores (PRSs) are promising for risk stratification but are primarily developed in European-ancestry populations. We aimed to develop and validate single- and multi-ancestry PRSs for East Asian never-smokers to improve lung adenocarcinoma risk prediction. Methods PRSs were developed using genome-wide association study summary statistics from East Asian (8002 cases; 20 782 controls) and European (2058 cases; 5575 controls) populations. Single-ancestry models included a PRS with 25 single nucleotide variants (PRS-25), PRS clumping and thresholding (CT), and LDpred2; multi-ancestry models included LDpred2 plus European PRS-128, PRS-CSx, and clumping and thresholding, super learning, empirical Bayes (CT-SLEB). Performance was evaluated in independent East Asian data from the Female Lung Cancer Consortium and externally validated in the Nanjing Lung Cancer Cohort. We assessed predictive accuracy via area under the curve operating procedure (AUC), with 10-year and (ages 30-80 years) absolute risks estimates. Results The best multi-ancestry PRS, using East Asian and European data via CT-SLEB, outperformed the best East Asian-only (LDpred2; AUC = 0.629, 95% confidence interval [CI] = 0.618 to 0.641), achieving an AUC of 0.640 (95% CI = 0.629 to 0.653) and odds ratio of 1.71 (95% CI = 1.61 to 1.82) per standard deviation increase. Nanjing Lung Cancer Cohort validation confirmed robust performance (AUC = 0.649, 95% CI = 0.623 to 0.676). The top 20% PRS group had a 3.92-fold higher lung adenocarcinoma risk than the bottom 20%. Further, the top 5% PRS group reached a 6.69% lifetime absolute risk. Notably, this group reached the average population 10-year lung adenocarcinoma risk at age 50 years (0.42%) by age 41 years, 9 years earlier. Conclusions Multi-ancestry PRS approaches enhance lung adenocarcinoma risk stratification in East Asian never-smokers, with consistent external validation, suggesting future clinical utility.

Abstract BACKGROUND High-grade glioma (HGG) is the most common and aggressive form of primary brain cancer. Standard of care surgery employs Magnetic Resonance Imaging to identify the contrast enhancing (CE) core for resection, leaving a non-enhancing (NE) invasive rim of tumor which contributes to recurrence. The hypoxic tumor microenvironment is thought to drive invasion into healthy brain, and understanding cellular responses to hypoxia may provide insight into combatting recurrence. Here, we apply a graph network to integrate regionally heterogeneous data from a unique dataset of spatially-localized CE and NE biopsies from a large cohort of patients with HGG to identify gene networks regulating hypoxia response and invasion. METHODS RNA-sequencing of 1159 differentially expressed genes across 159 IDH-wt biopsies (74 patients) from NE/CE regions were represented as nodes in a Neo4j graph network. Candidate genes for network random walks occurred exclusively within NE communities detected by modularity optimization and were localized by cell-type using single-cell RNA-sequencing (scRNAseq). CSF2RA was selected as source (walk length = 5; 10,000 walks/source node). Frequently occurring genes (one-tailed test) in walks were ranked by correlation with CSF2RA expression in NE samples and analyzed by gene set enrichment analysis for biological processes (clusterProfiler). CSF2RA and hypoxia gene expression were visualized with DoHeatmap in composite NE/CE scRNAseq profiles. RESULTS CSF2RA, an immune receptor, was upregulated (logFC = 0.789) in NE samples, scRNAseq myeloid populations, and Glycolytic/Plurimetabolic glioma subtype. CSF2RA random walk genes (n=46) enriched metabolic/proliferative signatures. CSF2RA expression correlated negatively with hypoxia genes globally and regionally (NE/CE), while NE/CE scRNAseq revealed distinct cellular expression patterns of CSF2RA and hypoxia response genes. CONCLUSION CSF2RA expression may prime NE glioma cells for invasion by altering hypoxia response. The genetic network bridging CSF2RA and hypoxia response may be a targetable mechanism of tumor recurrence.

Abstract BACKGROUND High-grade glioma (HGG) is the most common and aggressive form of primary brain cancer with poor prognosis. Standard of care surgery employs Magnetic Resonance Imaging (MRI) to identify the contrast enhancing (CE) core for resection, leaving behind an unresected non-enhancing (NE) invasive rim of tumor which borders healthy brain tissue and invariably contributes to recurrence. Recent multi-omic analyses identified predominant neuronal and immune signatures within NE tumor; understanding these signatures could provide key insights to combatting recurrence and improving patient outcomes. To this end, we utilized a novel graph database method to integrate multi-omic data from a unique dataset of multi-regional spatially-localized biopsies from the CE core and invasive NE rim from a large cohort of patients with HGG. MATERIAL AND METHODS 159 IDH-wildtype tumor biopsy samples from the NE and CE regions of 74 glioma patients were profiled by bulk RNA sequencing. Differentially expressed genes (DEG) (n=1159) across NE and CE samples were selected using EdgeR. Samples were hierarchically clustered (method = median) for each DEG and represented as nodes in a Neo4j graph network. Edges were generated between nodes sharing ≥ 1 sample. Modularity optimization was used to identify communities of densely connected nodes. Genes occurring exclusively within NE communities were localized by cell type using external single-cell RNA sequencing datasets and selected as candidate source nodes for NE network random walks. RBFOX1 was selected as a source (walk length = 5; 10,000 walks/source node). Random walks traversing samples with significantly elevated RBFOX1 expression (one-tailed test) were kept for analysis. Genes occurring at significantly high frequency (one-tailed test) within random walks were ranked by correlation with RBFOX1 expression in NE samples and analyzed by gene set enrichment analysis (GSEA) for biological processes using clusterProfiler. RESULTS RBFOX1, a neurodevelopmental regulator, was upregulated (logFC = 1.85) in NE samples and uniquely occurred in NE communities. We show RBFOX1 is differential expressed in the aggressive Neuronal glioma subtype and enriched in neuroepithelial/astrocyte subpopulations across single cell HGG datasets. NE network random walks identified 201 DEGs associated with RBFOX1. GSEA revealed upregulation of neurodevelopmental signatures and downregulation of proliferative pathways amongst RBFOX1’s network. CONCLUSION RBFOX1’s genetic network within NE tumor represents neurodevelopmental processes potentially underlying glioma invasiveness within the NE region. As such, understanding the mechanistic role of key genes within the RBFOX1 network may provide insight into novel strategies for targeting these invasive cells within an aggressive subtype of glioma in the NE rim.

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

AbstractDeficient intake of micronutrients involved in one‐carbon metabolism (eg, choline, methionine, vitamin B12 and folic acid) leads to hepatocellular carcinoma (HCC) development in rodents, but is under‐investigated in humans. We investigated the association between one‐carbon metabolism‐related micronutrient intake and HCC risk in a prospective cohort of 494 860 participants with 16 years of follow‐up in the NIH‐AARP study. Dietary intakes and supplement use were ascertained at baseline using a food‐frequency questionnaire. Total intake (diet plus supplements) of the following one‐carbon metabolism‐related micronutrients were calculated: folate, methionine and vitamins B2 (riboflavin), B3 (niacin), B6 and B12. These micronutrients were examined both individually and simultaneously, with adjustment for covariates. Cox proportional hazard models were used to calculate hazard ratios (HRs) and 95% confidence intervals (CIs). Over the 16‐year follow‐up period, 647 incident HCC cases were diagnosed. When examined individually, higher total vitamin B3 intake was associated with a lower HCC risk (HRQ5 vs Q1 = 0.60; 95% CI = 0.42‐0.85; Ptrend = .008), and the association remained significant when all six micronutrients were examined simultaneously (HRQ5 vs Q1 = 0.32; 95% CI = 0.18‐0.55; Ptrend < .0001). Among participants with >3 years of follow‐up, higher total vitamin B3 intake was again associated with lower risk (HRQ5 vs Q1 = 0.37; 95% CI = 0.20‐0.68; Ptrend = .001), whereas higher total vitamin B6 intake was associated with higher risk (HRQ5 vs Q1 = 2.04; 95% CI = 1.02‐4.07; Ptrend = .04). Restricted cubic spline analyses showed a dose‐response inverse association between total vitamin B3 intake and HCC risk, and dose‐response positive association between total vitamin B6 intake and HCC risk. The study suggests that higher vitamin B3 intake is associated with lower HCC risk, whereas higher vitamin B6 intake is associated with increased risk.

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