White, Emily

Abstract Background Physical activity (PA) is an established protective factor for colorectal cancer (CRC), but it is unclear if genetic variants modify this effect. To investigate this possibility, we conducted a genome-wide gene–PA interaction analysis. Methods Using logistic regression and two-step and joint tests, we analyzed interactions between common genetic variants across the genome and PA in relation to CRC risk. Self-reported PA levels were categorized as active (≥ 8.75 MET-h/wk) vs. inactive (< 8.75 MET-h/wk) and as study- and sex-specific quartiles of activity. Results PA had an overall protective effect on CRC (OR [active vs. inactive] = 0.85; 95%CI = 0.81–0.90). The two-step GxE method identified an interaction between rs4779584, an intergenic variant near the GREM1 and SCG5 genes, and PA for CRC risk (p-interaction = 2.6×10 − 8 ). Stratification by genotype at this locus showed a significant reduction in CRC risk by 20% in active vs. inactive participants with the CC genotype (OR = 0.80; 95%CI = 0.75–0.85), but no significant PA–CRC association among CT or TT carriers. When PA was modeled as quartiles, the 1-d.f. GxE test identified that rs56906466, an intergenic variant near the KCNG1 gene, modified the association between PA and CRC (p-interaction = 3.5×10 − 8 ). Stratification at this locus showed that increase in PA (highest vs. lowest quartile) was associated with a lower CRC risk solely among TT carriers (OR = 0.77; 95%CI = 0.72–0.82). Conclusions In summary, we identified two genetic variants that modified the association between PA and CRC risk. One of them, related to GREM1 and SCG5 , suggests that the bone morphogenetic protein (BMP)-related, inflammatory, and/or insulin signaling pathways may be associated with the protective influence of PA on colorectal carcinogenesis.

ABSTRACT Lung cancer is the leading cause of cancer-related deaths worldwide and, despite treatment advances, immune suppression remains an obstacle to effective therapy. Effector CD4 + T cells (CD4 + Teffs) are critical for antitumor immunity, but their function is often inhibited by regulatory T cells (Tregs), which accumulate in lung tumors and perform suppressive functions through multiple mechanisms. This suppression leads to tumor progression and poor patient outcomes. However, the mechanisms underlying Treg-mediated suppression are not fully understood. Here, we identify the extracellular ATP receptor P2RX7 as a key regulator of Treg function in lung tumors. Using a murine lung cancer model induced by Lewis lung carcinoma cells, we demonstrate that P2RX7 enhances the suppressive capacity of tumor-infiltrating Tregs, promoting tumor growth. In T cell-specific P2RX7-KO mice, reduced Treg infiltration was accompanied by increased CD4 + Teff accumulation and improved tumor control. Treg-specific P2RX7-KO mice exhibit reduced tumor growth, confirming a cell-intrinsic role of P2RX7 in Tregs. Suppression assays revealed that tumor-infiltrating WT Tregs have greater suppressive activity compared to P2RX7-KO Tregs, which failed to inhibit type 1 and Tfh-like responses. This was associated with increased tumor-specific IgG production by lung B cells in P2RX7-KO mice. We also observed that WT Tregs express higher levels of the immunosuppressive surface molecule CTLA-4 when compared to P2RX7-KO Tregs. In summary, we show that P2RX7 expression on Tregs is essential for their suppressive function in lung cancer, and targeting of P2RX7 may constitute a novel strategy to improve lung cancer treatment by alleviating Treg-mediated immune suppression.

Abstract Observational studies have suggested a protective role for eosinophils in colorectal cancer (CRC) development and implicated neutrophils, but the causal relationships remain unclear. Here, we aimed to estimate the causal effect of circulating white blood cell (WBC) counts (N = ~550 000) for basophils, eosinophils, monocytes, lymphocytes and neutrophils on CRC risk (N = 52 775 cases and 45 940 controls) using Mendelian randomisation (MR). For comparison, we also examined this relationship using individual‐level data from UK Biobank (4043 incident CRC cases and 332 773 controls) in a longitudinal cohort analysis. The inverse‐variance weighted (IVW) MR analysis suggested a protective effect of increased basophil count and eosinophil count on CRC risk [OR per 1‐SD increase: 0.88, 95% CI: 0.78‐0.99, P = .04; OR: 0.93, 95% CI: 0.88‐0.98, P = .01]. The protective effect of eosinophils remained [OR per 1‐SD increase: 0.88, 95% CI: 0.80‐0.97, P = .01] following adjustments for all other WBC subtypes, to account for genetic correlation between the traits, using multivariable MR. A protective effect of increased lymphocyte count on CRC risk was also found [OR: 0.84, 95% CI: 0.76‐0.93, P = 6.70e‐4] following adjustment. Consistent with MR results, a protective effect for eosinophils in the cohort analysis in the fully adjusted model [RR per 1‐SD increase: 0.96, 95% CI: 0.93‐0.99, P = .02] and following adjustment for the other WBC subtypes [RR: 0.96, 95% CI: 0.93‐0.99, P = .001] was observed. Our study implicates peripheral blood immune cells, in particular eosinophils and lymphocytes, in CRC development, highlighting a need for mechanistic studies to interrogate these relationships.

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.

Salicylic acid (SA) has observationally been shown to decrease colorectal cancer (CRC) risk. Aspirin (acetylsalicylic acid, that rapidly deacetylates to SA) is an effective primary and secondary chemopreventive agent. Through a Mendelian randomization (MR) approach, we aimed to address whether levels of SA affected CRC risk, stratifying by aspirin use. A two-sample MR analysis was performed using GWAS summary statistics of SA (INTERVAL and EPIC-Norfolk, N = 14,149) and CRC (CCFR, CORECT, GECCO and UK Biobank, 55,168 cases and 65,160 controls). The DACHS study (4410 cases and 3441 controls) was used for replication and stratification of aspirin-use. SNPs proxying SA were selected via three methods: (1) functional SNPs that influence the activity of aspirin-metabolising enzymes; (2) pathway SNPs present in enzymes’ coding regions; and (3) genome-wide significant SNPs. We found no association between functional SNPs and SA levels. The pathway and genome-wide SNPs showed no association between SA and CRC risk (OR: 1.03, 95% CI: 0.84–1.27 and OR: 1.08, 95% CI: 0.86–1.34, respectively). Results remained unchanged upon aspirin use stratification. We found little evidence to suggest that an SD increase in genetically predicted SA protects against CRC risk in the general population and upon stratification by aspirin use.

Abstract Background: Evidence for aspirin's chemopreventative properties on colorectal cancer (CRC) is substantial, but its mechanism of action is not well-understood. We combined a proteomic approach with Mendelian randomization (MR) to identify possible new aspirin targets that decrease CRC risk. Methods: Human colorectal adenoma cells (RG/C2) were treated with aspirin (24 hours) and a stable isotope labeling with amino acids in cell culture (SILAC) based proteomics approach identified altered protein expression. Protein quantitative trait loci (pQTLs) from INTERVAL (N = 3,301) and expression QTLs (eQTLs) from the eQTLGen Consortium (N = 31,684) were used as genetic proxies for protein and mRNA expression levels. Two-sample MR of mRNA/protein expression on CRC risk was performed using eQTL/pQTL data combined with CRC genetic summary data from the Colon Cancer Family Registry (CCFR), Colorectal Transdisciplinary (CORECT), Genetics and Epidemiology of Colorectal Cancer (GECCO) consortia and UK Biobank (55,168 cases and 65,160 controls). Results: Altered expression was detected for 125/5886 proteins. Of these, aspirin decreased MCM6, RRM2, and ARFIP2 expression, and MR analysis showed that a standard deviation increase in mRNA/protein expression was associated with increased CRC risk (OR: 1.08, 95% CI, 1.03–1.13; OR: 3.33, 95% CI, 2.46–4.50; and OR: 1.15, 95% CI, 1.02–1.29, respectively). Conclusions: MCM6 and RRM2 are involved in DNA repair whereby reduced expression may lead to increased DNA aberrations and ultimately cancer cell death, whereas ARFIP2 is involved in actin cytoskeletal regulation, indicating a possible role in aspirin's reduction of metastasis. Impact: Our approach has shown how laboratory experiments and population-based approaches can combine to identify aspirin-targeted proteins possibly affecting CRC risk.

Abstract Background: Accumulating evidence suggests a relationship between endometrial cancer and ovarian cancer. Independent genome-wide association studies (GWAS) for endometrial cancer and ovarian cancer have identified 16 and 27 risk regions, respectively, four of which overlap between the two cancers. We aimed to identify joint endometrial and ovarian cancer risk loci by performing a meta-analysis of GWAS summary statistics from these two cancers. Methods: Using LDScore regression, we explored the genetic correlation between endometrial cancer and ovarian cancer. To identify loci associated with the risk of both cancers, we implemented a pipeline of statistical genetic analyses (i.e., inverse-variance meta-analysis, colocalization, and M-values) and performed analyses stratified by subtype. Candidate target genes were then prioritized using functional genomic data. Results: Genetic correlation analysis revealed significant genetic correlation between the two cancers (rG = 0.43, P = 2.66 × 10−5). We found seven loci associated with risk for both cancers (PBonferroni &lt; 2.4 × 10−9). In addition, four novel subgenome-wide regions at 7p22.2, 7q22.1, 9p12, and 11q13.3 were identified (P &lt; 5 × 10−7). Promoter-associated HiChIP chromatin loops from immortalized endometrium and ovarian cell lines and expression quantitative trait loci data highlighted candidate target genes for further investigation. Conclusions: Using cross-cancer GWAS meta-analysis, we have identified several joint endometrial and ovarian cancer risk loci and candidate target genes for future functional analysis. Impact: Our research highlights the shared genetic relationship between endometrial cancer and ovarian cancer. Further studies in larger sample sets are required to confirm our findings.

Abstract Background: Whether circulating polyunsaturated fatty acid (PUFA) levels are associated with pancreatic cancer risk is uncertain. Mendelian randomization (MR) represents a study design using genetic instruments to better characterize the relationship between exposure and outcome. Methods: We utilized data from genome-wide association studies within the Pancreatic Cancer Cohort Consortium and Pancreatic Cancer Case–Control Consortium, involving approximately 9,269 cases and 12,530 controls of European descent, to evaluate associations between pancreatic cancer risk and genetically predicted plasma n-6 PUFA levels. Conventional MR analyses were performed using individual-level and summary-level data. Results: Using genetic instruments, we did not find evidence of associations between genetically predicted plasma n-6 PUFA levels and pancreatic cancer risk [estimates per one SD increase in each PUFA-specific weighted genetic score using summary statistics: linoleic acid odds ratio (OR) = 1.00, 95% confidence interval (CI) = 0.98–1.02; arachidonic acid OR = 1.00, 95% CI = 0.99–1.01; and dihomo-gamma-linolenic acid OR = 0.95, 95% CI = 0.87–1.02]. The OR estimates remained virtually unchanged after adjustment for covariates, using individual-level data or summary statistics, or stratification by age and sex. Conclusions: Our results suggest that variations of genetically determined plasma n-6 PUFA levels are not associated with pancreatic cancer risk. Impact: These results suggest that modifying n-6 PUFA levels through food sources or supplementation may not influence risk of pancreatic cancer.

Abstract Pancreatic ductal adenocarcinoma (PDAC) is one of most lethal malignancies with few known risk factors and biomarkers. Identification of disease biomarkers is critical for understanding the pathogenesis of this cancer and identifying high risk individuals for close surveillance. Several blood protein biomarkers have been linked to PDAC in previous studies, but these studies have assessed only a limited number of biomarkers usually in small samples. To identify novel circulating protein biomarkers of PDAC, we studied 8,280 cases and 6,728 controls of European descent from the Pancreatic Cancer Cohort Consortium and the Pancreatic Cancer Case-Control Consortium, by using genetic instruments. Protein quantitative trait loci (pQTLs) for 1,226 plasma proteins identified in a large INTERVAL study of 3,301 healthy European descendants were used as instruments to evaluate associations between genetically predicted protein levels and PDAC. For proteins showing a significant association, we further conducted conditional analysis with adjustments for previously identified risk variants to assess whether the observed associations between genetically predicted protein concentrations and PDAC risk were independent of the risk variants identified in genome-wide association studies (GWAS). Furthermore, for the proteins that were associated with PDAC risk, we performed an enrichment analysis of the genes encoding these proteins to examine whether they are enriched in specific pathways, functions or networks. We observed associations between predicted concentrations of 38 proteins and PDAC risk at a false discovery rate of &lt; 0.05, including those of 23 proteins that showed a significant association even after Bonferroni correction (4.08 × 10−5). These include Histo-blood group ABO system transferase encoded by ABO, which has been previously implicated as a potential target gene of PDAC risk variant identified in GWAS. Eight of the identified proteins (Beta-crystallin B2, Dedicator of cytokinesis protein 9, VIP36-like protein, Erythrocyte band 7 integral membrane protein, Tensin-2, Transmembrane protease serine 11D, Alcohol dehydrogenase 1B, and C-X-C motif chemokine 10) were associated with PDAC risk after conditioning on previously reported pancreatic cancer risk variants (odds ratios ranged from 0.79 to 1.52, P-values from 1.28 × 10−3 to 6.47 × 10−4). Pathway enrichment analysis showed that the encoding genes for the implicated proteins were significantly enriched in cancer-related pathways, such as STAT3 and IL-15 production. In conclusion, we identified 38 protein biomarker candidates for PDAC risk, which if validated by additional studies, may contribute to the etiological understanding of PDAC tumor development. Citation Format: Jingjing Zhu, Xiang Shu, Xingyi Guo, Duo Liu, Jiandong bao, Roger Milne, Graham G Giles, Chong Wu, Mengmeng Du, Emily White, Harvey A Risch, Nuria Malats, Eric J. Duell, Phyllis J. Goodman, Donghui Li, Paige Bracci, Verena Katzke, Rachel E Neale, Steven Gallinger, Stephen Van Den Eeden, Alan Arslan, Federico Canzian, Charles Kooperberg, Brian Wolpin, Laura Beane-Freeman, Ghislaine Scelo, Kala Visvanatha, Christopher A. Haiman, Loïc Le Marchand, Herbert Yu, Gloria M Petersen, Rachael Stolzenberg-Solomon, Alison P Klein, Laufey T Amundadottir, Qiuyin Cai, Jirong Long, Xiao-Ou Shu, Wei Zheng, Lang Wu. Associations between genetically predicted blood protein biomarkers and pancreatic ductal adenocarcinoma risk [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr 1200.

Abstract Background: Regular use of nonsteroidal anti-inflammatory drugs (NSAID) is associated with lower risk of colorectal cancer. Genome-wide interaction analysis on single variants (G × E) has identified several SNPs that may interact with NSAIDs to confer colorectal cancer risk, but variations in gene expression levels may also modify the effect of NSAID use. Therefore, we tested interactions between NSAID use and predicted gene expression levels in relation to colorectal cancer risk. Methods: Genetically predicted gene expressions were tested for interaction with NSAID use on colorectal cancer risk among 19,258 colorectal cancer cases and 18,597 controls from 21 observational studies. A Mixed Score Test for Interactions (MiSTi) approach was used to jointly assess G × E effects which are modeled via fixed interaction effects of the weighted burden within each gene set (burden) and residual G × E effects (variance). A false discovery rate (FDR) at 0.2 was applied to correct for multiple testing. Results: Among the 4,840 genes tested, genetically predicted expression levels of four genes modified the effect of any NSAID use on colorectal cancer risk, including DPP10 (PG×E = 1.96 × 10−4), KRT16 (PG×E = 2.3 × 10−4), CD14 (PG×E = 9.38 × 10−4), and CYP27A1 (PG×E = 1.44 × 10−3). There was a significant interaction between expression level of RP11-89N17 and regular use of aspirin only on colorectal cancer risk (PG×E = 3.23 × 10−5). No interactions were observed between predicted gene expression and nonaspirin NSAID use at FDR &lt; 0.2. Conclusions: By incorporating functional information, we discovered several novel genes that interacted with NSAID use. Impact: These findings provide preliminary support that could help understand the chemopreventive mechanisms of NSAIDs on colorectal cancer.