Marchand, Loïc Le

Abstract Background The associations between different types of diabetes, characterized by distinct pathophysiology and genetic architecture, and pancreatic ductal adenocarcinoma (PDAC) risk are not understood. Methods We investigated associations of genetic susceptibility to type 2 diabetes (T2D), 8 T2D mechanistic clusters, type 1 diabetes (T1D), and maturity-onset diabetes of the young (MODY) with PDAC risk. We used genome-wide association study (GWAS) summary-level statistics for T2D (242 283 cases, 1 569 734 controls), T1D (18 942 cases, 501 638 controls), and PDAC (10 244 cases and 360 535 controls) in individuals of European ancestry. Results Two-sample Mendelian randomization (MR) using the Robust Adjusted Profile Score (MR-RAPS) method indicated that genetically predicted T2D was associated with PDAC risk (OR = 1.10; 95% CI = 1.05 to 1.15), particularly the T2D obesity (OR = 1.28; 95% CI = 1.15 to 1.42) and lipodystrophy (OR = 1.25; 95% CI = 1.03 to 1.51) clusters. No association was observed for T1D with PDAC risk (OR = 1.01; 95% CI = 0.99 to 1.02). Pathway/gene-set analysis using the summary-based Adaptive Rank Truncated Product (sARTP) method revealed a significant association between the MODY gene-sets and PDAC risk (P = 1.5 × 10−8), which remained after excluding 20 known PDAC GWAS loci (P = 7.6 × 10−4). HNF1A, FOXA3, and HNF4A were the top contributing genes after excluding the previously identified GWAS loci regions. Conclusions Our results from this genetic association study support that T2D, particularly the obesity and lipodystrophy mechanistic clusters, and MODY genomic susceptibility regions play a role in the etiology of PDAC.

Abstract Menopausal users of hormone replacement therapy (HRT) are at increased breast cancer risk and decreased colorectal cancer (CRC) risk compared with individuals who have never used HRT, but these opposing associations may differ by familial risk of breast cancer and CRC. We harmonized data from 3 cohorts and generated separate breast cancer and CRC familial risk scores based on cancer family history. We defined moderate or strong family history as a risk score of 0.4 or higher, where 0.4 was equivalent to a 50-year-old woman with 1 parent diagnosed with either breast cancer or CRC at 55 years of age. Of 24 486 women assessed, 1243 and 405 were diagnosed with incident breast cancer and CRC, respectively. For breast cancer, menopausal HRT ever use versus never use hazard ratios were 1.27 (95% CI = 1.11 to 1.45) for a breast cancer familial risk score below 0.4 and 1.01 (95% CI = 0.82 to 1.25) for a breast cancer familial risk score of 0.4 or higher (Pdifference = .08). For CRC, menopausal HRT hazard ratios were 0.63 (95% CI = 0.50 to 0.78) for a CRC familial risk score below 0.4 and 1.21 (95% CI = 0.73 to 2.00) for a CRC familial risk score of 0.4 or higher (Pdifference = .03). Associations with menopausal HRT use that apply to the general population may not hold for women at moderate or strong familial risk of these cancers.

Abstract Background: Epidemiological and laboratory-based studies have provided conflicting evidence for a role of ghrelin in colorectal cancer development. We conducted two-sample Mendelian randomization (MR) analyses to evaluate evidence for an association of circulating ghrelin and colorectal cancer risk overall and by sex, cancer subsite, and age at diagnosis. Methods: Genetic instruments proxying plasma total ghrelin levels were obtained from a recent genome-wide association study of 54,219 participants. Summary data for colorectal cancer risk were obtained from a recent meta-analysis of several genetic consortia (up to 73,673 cases and 86,854 controls). A two-sample MR approach and several sensitivity analyses were applied. Results: We found no evidence for an association of genetically predicted plasma total ghrelin levels and colorectal cancer risk (0.95, 95% confidence interval, 0.81–1.12; R2 of ghrelin genetic instruments: 4.6%), with similarly null results observed when stratified by sex, anatomical subsite, and for early-onset colorectal cancer. Conclusions: Our study suggests that plasma ghrelin levels are unlikely to have a causal relationship with overall, early-onset, and sex- and cancer subsite–stratified colorectal cancer risk. Impact: This large-scale analysis adds to the growing body of evidence that plasma total ghrelin levels are not associated with colorectal cancer risk.

Abstract Background: Pancreatic cancer is a leading cause of cancer-related death globally. Risk factors for pancreatic cancer include common genetic variants and potentially heavy alcohol consumption. We assessed if genetic variants modify the association between heavy alcohol consumption and pancreatic cancer risk. Methods: We conducted a genome-wide interaction analysis of single-nucleotide polymorphisms (SNP) by heavy alcohol consumption (more than three drinks per day) for pancreatic cancer in European ancestry populations from genome-wide association studies. Our analysis included 3,707 cases and 4,167 controls from case–control studies and 1,098 cases and 1,162 controls from cohort studies. Fixed-effect meta-analyses were conducted. Results: A potential novel region of association on 10p11.22, lead SNP rs7898449 (interaction P value (Pinteraction) = 5.1 × 10−8 in the meta-analysis; Pinteraction = 2.1 × 10−9 in the case–control studies; Pinteraction = 0.91 in the cohort studies), was identified. An SNP correlated with this lead SNP is an expression quantitative trait locus for the neuropilin 1 gene. Of the 17 genomic regions with genome-wide significant evidence of association with pancreatic cancer in prior studies, we observed suggestive evidence that heavy alcohol consumption modified the association for one SNP near LINC00673, rs11655237 on 17q25.1 (Pinteraction = 0.004). Conclusions: We identified a novel genomic region that may be associated with pancreatic cancer risk in conjunction with heavy alcohol consumption located near an expression quantitative trait locus for neuropilin 1, a protein that plays an important role in the development and progression of pancreatic cancer. Impact: This work can provide insights into the etiology of pancreatic cancer, particularly in heavy drinkers.

Abstract Deep learning may detect biologically important signals embedded in tumor morphologic features that confer distinct prognoses. Tumor morphologic features were quantified to enhance patient risk stratification within DNA mismatch repair (MMR) groups using deep learning. Using a quantitative segmentation algorithm (QuantCRC) that identifies 15 distinct morphologic features, we analyzed 402 resected stage III colon carcinomas [191 deficient (d)-MMR; 189 proficient (p)-MMR] from participants in a phase III trial of FOLFOX-based adjuvant chemotherapy. Results were validated in an independent cohort (176 d-MMR; 1,094 p-MMR). Association of morphologic features with clinicopathologic variables, MMR, KRAS, BRAFV600E, and time-to-recurrence (TTR) was determined. Multivariable Cox proportional hazards models were developed to predict TTR. Tumor morphologic features differed significantly by MMR status. Cancers with p-MMR had more immature desmoplastic stroma. Tumors with d-MMR had increased inflammatory stroma, epithelial tumor-infiltrating lymphocytes (TIL), high-grade histology, mucin, and signet ring cells. Stromal subtype did not differ by BRAFV600E or KRAS status. In p-MMR tumors, multivariable analysis identified tumor-stroma ratio (TSR) as the strongest feature associated with TTR [HRadj 2.02; 95% confidence interval (CI), 1.14–3.57; P = 0.018; 3-year recurrence: 40.2% vs. 20.4%; Q1 vs. Q2–4]. Among d-MMR tumors, extent of inflammatory stroma (continuous HRadj 0.98; 95% CI, 0.96–0.99; P = 0.028; 3-year recurrence: 13.3% vs. 33.4%, Q4 vs. Q1) and N stage were the most robust prognostically. Association of TSR with TTR was independently validated. In conclusion, QuantCRC can quantify morphologic differences within MMR groups in routine tumor sections to determine their relative contributions to patient prognosis, and may elucidate relevant pathophysiologic mechanisms driving prognosis. Significance: A deep learning algorithm can quantify tumor morphologic features that may reflect underlying mechanisms driving prognosis within MMR groups. TSR was the most robust morphologic feature associated with TTR in p-MMR colon cancers. Extent of inflammatory stroma and N stage were the strongest prognostic features in d-MMR tumors. TIL density was not independently prognostic in either MMR group.

Abstract Background: Clinical, molecular, and genetic epidemiology studies displayed remarkable differences between ever- and never-smoking lung cancer. Methods: We conducted a stratified multi-population (European, East Asian, and African descent) association study on 44,823 ever-smokers and 20,074 never-smokers to identify novel variants that were missed in the non-stratified analysis. Functional analysis including expression quantitative trait loci (eQTL) colocalization and DNA damage assays, and annotation studies were conducted to evaluate the functional roles of the variants. We further evaluated the impact of smoking quantity on lung cancer risk for the variants associated with ever-smoking lung cancer. Results: Five novel independent loci, GABRA4, intergenic region 12q24.33, LRRC4C, LINC01088, and LCNL1 were identified with the association at two or three populations (P < 5 × 10−8). Further functional analysis provided multiple lines of evidence suggesting the variants affect lung cancer risk through excessive DNA damage (GABRA4) or cis-regulation of gene expression (LCNL1). The risk of variants from 12 independent regions, including the well-known CHRNA5, associated with ever-smoking lung cancer was evaluated for never-smokers, light-smokers (packyear ≤ 20), and moderate-to-heavy-smokers (packyear > 20). Different risk patterns were observed for the variants among the different groups by smoking behavior. Conclusions: We identified novel variants associated with lung cancer in only ever- or never-smoking groups that were missed by prior main-effect association studies. Impact: Our study highlights the genetic heterogeneity between ever- and never-smoking lung cancer and provides etiologic insights into the complicated genetic architecture of this deadly cancer.

Abstract Background: DEPendency of association on the number of Top Hits (DEPTH) is an approach to identify candidate susceptibility regions by considering the risk signals from overlapping groups of sequential variants across the genome. Methods: We applied a DEPTH analysis using a sliding window of 200 SNPs to colorectal cancer data from the Colon Cancer Family Registry (CCFR; 5,735 cases and 3,688 controls), and Genetics and Epidemiology of Colorectal Cancer Consortium (GECCO; 8,865 cases and 10,285 controls) studies. A DEPTH score > 1 was used to identify candidate susceptibility regions common to both analyses. We compared DEPTH results against those from conventional genome-wide association study (GWAS) analyses of these two studies as well as against 132 published susceptibility regions. Results: Initial DEPTH analysis revealed 2,622 (CCFR) and 3,686 (GECCO) candidate susceptibility regions, of which 569 were common to both studies. Bootstrapping revealed 40 and 49 candidate susceptibility regions in the CCFR and GECCO data sets, respectively. Notably, DEPTH identified at least 82 regions that would not be detected using conventional GWAS methods, nor had they been identified by previous colorectal cancer GWASs. We found four reproducible candidate susceptibility regions (2q22.2, 2q33.1, 6p21.32, 13q14.3). The highest DEPTH scores were in the human leukocyte antigen locus at 6p21 where the strongest associated SNPs were rs762216297, rs149490268, rs114741460, and rs199707618 for the CCFR data, and rs9270761 for the GECCO data. Conclusions: DEPTH can identify candidate susceptibility regions for colorectal cancer not identified using conventional analyses of larger datasets. Impact: DEPTH has potential as a powerful complementary tool to conventional GWAS analyses for discovering susceptibility regions within the genome.

Abstract Background: There are conflicting data on whether nonalcoholic fatty liver disease (NAFLD) is associated with susceptibility to pancreatic cancer. Using Mendelian randomization (MR), we investigated the relationship between genetic predisposition to NAFLD and risk for pancreatic cancer. Methods: Data from genome-wide association studies (GWAS) within the Pancreatic Cancer Cohort Consortium (PanScan; cases n = 5,090, controls n = 8,733) and the Pancreatic Cancer Case Control Consortium (PanC4; cases n = 4,163, controls n = 3,792) were analyzed. We used data on 68 genetic variants with four different MR methods [inverse variance weighting (IVW), MR-Egger, simple median, and penalized weighted median] separately to predict genetic heritability of NAFLD. We then assessed the relationship between each of the four MR methods and pancreatic cancer risk, using logistic regression to calculate ORs and 95% confidence intervals (CI), adjusting for PC risk factors, including obesity and diabetes. Results: No association was found between genetically predicted NAFLD and pancreatic cancer risk in the PanScan or PanC4 samples [e.g., PanScan, IVW OR, 1.04; 95% confidence interval (CI), 0.88–1.22; MR-Egger OR, 0.89; 95% CI, 0.65–1.21; PanC4, IVW OR, 1.07; 95% CI, 0.90–1.27; MR-Egger OR, 0.93; 95% CI, 0.67–1.28]. None of the four MR methods indicated an association between genetically predicted NAFLD and pancreatic cancer risk in either sample. Conclusions: Genetic predisposition to NAFLD is not associated with pancreatic cancer risk. Impact: Given the close relationship between NAFLD and metabolic conditions, it is plausible that any association between NAFLD and pancreatic cancer might reflect host metabolic perturbations (e.g., obesity, diabetes, or metabolic syndrome) and does not necessarily reflect a causal relationship between NAFLD and pancreatic cancer.

Abstract Colorectal cancer risk can be impacted by genetic, environmental, and lifestyle factors, including diet and obesity. Gene-environment interactions (G × E) can provide biological insights into the effects of obesity on colorectal cancer risk. Here, we assessed potential genome-wide G × E interactions between body mass index (BMI) and common SNPs for colorectal cancer risk using data from 36,415 colorectal cancer cases and 48,451 controls from three international colorectal cancer consortia (CCFR, CORECT, and GECCO). The G × E tests included the conventional logistic regression using multiplicative terms (one degree of freedom, 1DF test), the two-step EDGE method, and the joint 3DF test, each of which is powerful for detecting G × E interactions under specific conditions. BMI was associated with higher colorectal cancer risk. The two-step approach revealed a statistically significant G×BMI interaction located within the Formin 1/Gremlin 1 (FMN1/GREM1) gene region (rs58349661). This SNP was also identified by the 3DF test, with a suggestive statistical significance in the 1DF test. Among participants with the CC genotype of rs58349661, overweight and obesity categories were associated with higher colorectal cancer risk, whereas null associations were observed across BMI categories in those with the TT genotype. Using data from three large international consortia, this study discovered a locus in the FMN1/GREM1 gene region that interacts with BMI on the association with colorectal cancer risk. Further studies should examine the potential mechanisms through which this locus modifies the etiologic link between obesity and colorectal cancer. Significance: This gene-environment interaction analysis revealed a genetic locus in FMN1/GREM1 that interacts with body mass index in colorectal cancer risk, suggesting potential implications for precision prevention strategies.

Abstract Background The shared inherited genetic contribution to risk of different cancers is not fully known. In this study, we leverage results from 12 cancer genome-wide association studies (GWAS) to quantify pairwise genome-wide genetic correlations across cancers and identify novel cancer susceptibility loci. Methods We collected GWAS summary statistics for 12 solid cancers based on 376 759 participants with cancer and 532 864 participants without cancer of European ancestry. The included cancer types were breast, colorectal, endometrial, esophageal, glioma, head and neck, lung, melanoma, ovarian, pancreatic, prostate, and renal cancers. We conducted cross-cancer GWAS and transcriptome-wide association studies to discover novel cancer susceptibility loci. Finally, we assessed the extent of variant-specific pleiotropy among cancers at known and newly identified cancer susceptibility loci. Results We observed widespread but modest genome-wide genetic correlations across cancers. In cross-cancer GWAS and transcriptome-wide association studies, we identified 15 novel cancer susceptibility loci. Additionally, we identified multiple variants at 77 distinct loci with strong evidence of being associated with at least 2 cancer types by testing for pleiotropy at known cancer susceptibility loci. Conclusions Overall, these results suggest that some genetic risk variants are shared among cancers, though much of cancer heritability is cancer-specific and thus tissue-specific. The increase in statistical power associated with larger sample sizes in cross-disease analysis allows for the identification of novel susceptibility regions. Future studies incorporating data on multiple cancer types are likely to identify additional regions associated with the risk of multiple cancer types.