Gharahkhani, Puya

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.

Abstract Background Mendelian randomization (MR) studies are susceptible to metadata errors (e.g. incorrect specification of the effect allele column) and other analytical issues that can introduce substantial bias into analyses. We developed a quality control (QC) pipeline for the Fatty Acids in Cancer Mendelian Randomization Collaboration (FAMRC) that can be used to identify and correct for such errors. Methods We collated summary association statistics from fatty acid and cancer genome-wide association studies (GWAS) and subjected the collated data to a comprehensive QC pipeline. We identified metadata errors through comparison of study-specific statistics to external reference data sets (the National Human Genome Research Institute-European Bioinformatics Institute GWAS catalogue and 1000 genome super populations) and other analytical issues through comparison of reported to expected genetic effect sizes. Comparisons were based on three sets of genetic variants: (i) GWAS hits for fatty acids, (ii) GWAS hits for cancer and (iii) a 1000 genomes reference set. Results We collated summary data from 6 fatty acid and 54 cancer GWAS. Metadata errors and analytical issues with the potential to introduce substantial bias were identified in seven studies (11.6%). After resolving metadata errors and analytical issues, we created a data set of 219 842 genetic associations with 90 cancer types, generated in analyses of 566 665 cancer cases and 1 622 374 controls. Conclusions In this large MR collaboration, 11.6% of included studies were affected by a substantial metadata error or analytical issue. By increasing the integrity of collated summary data prior to their analysis, our protocol can be used to increase the reliability of downstream MR analyses. Our pipeline is available to other researchers via the CheckSumStats package (https://github.com/MRCIEU/CheckSumStats).

Abstract Background: Over 20 susceptibility single-nucleotide polymorphisms (SNP) have been identified for esophageal adenocarcinoma (EAC) and its precursor, Barrett esophagus (BE), explaining a small portion of heritability. Methods: Using genetic data from 4,323 BE and 4,116 EAC patients aggregated by international consortia including the Barrett's and Esophageal Adenocarcinoma Consortium (BEACON), we conducted a comprehensive transcriptome-wide association study (TWAS) for BE/EAC, leveraging Genotype Tissue Expression (GTEx) gene-expression data from six tissue types of plausible relevance to EAC etiology: mucosa and muscularis from the esophagus, gastroesophageal (GE) junction, stomach, whole blood, and visceral adipose. Two analytical approaches were taken: standard TWAS using the predicted gene expression from local expression quantitative trait loci (eQTL), and set-based SKAT association using selected eQTLs that predict the gene expression. Results: Although the standard approach did not identify significant signals, the eQTL set–based approach identified eight novel associations, three of which were validated in independent external data (eQTL SNP sets for EXOC3, ZNF641, and HSP90AA1). Conclusions: This study identified novel genetic susceptibility loci for EAC and BE using an eQTL set–based genetic association approach. Impact: This study expanded the pool of genetic susceptibility loci for EAC and BE, suggesting the potential of the eQTL set–based genetic association approach as an alternative method for TWAS analysis.

Abstract Background: Many loci have been found to be associated with risk of epithelial ovarian cancer (EOC). However, although there is considerable variation in progression-free survival (PFS), no loci have been found to be associated with outcome at genome-wide levels of significance. Methods: We carried out a genome-wide association study (GWAS) of PFS in 2,352 women with EOC who had undergone cytoreductive surgery and standard carboplatin/paclitaxel chemotherapy. Results: We found seven SNPs at 12q24.33 associated with PFS (P < 5 × 10–8), the top SNP being rs10794418 (HR = 1.24; 95% CI, 1.15–1.34; P = 1.47 × 10–8). High expression of a nearby gene, ULK1, is associated with shorter PFS in EOC, and with poor prognosis in other cancers. SNP rs10794418 is also associated with expression of ULK1 in ovarian tumors, with the allele associated with shorter PFS being associated with higher expression, and chromatin interactions were detected between the ULK1 promoter and associated SNPs in serous and endometrioid EOC cell lines. ULK1 knockout ovarian cancer cell lines showed significantly increased sensitivity to carboplatin in vitro. Conclusions: The locus at 12q24.33 represents one of the first genome-wide significant loci for survival for any cancer. ULK1 is a plausible candidate for the target of this association. Impact: This finding provides insight into genetic markers associated with EOC outcome and potential treatment options. See related commentary by Peres and Monteiro, p. 1604

Abstract Genome-wide association studies (GWAS) of esophageal adenocarcinoma (EAC) and its precursor, Barrett’s esophagus (BE), have uncovered significant genetic components of risk, but most heritability remains unexplained. Targeted assessment of genetic variation in biologically relevant pathways using novel analytical approaches may identify missed susceptibility signals. Central obesity, a key BE/EAC risk factor, is linked to systemic inflammation, altered hormonal signaling and insulin-like growth factor (IGF) axis dysfunction. Here, we assessed IGF-related genetic variation and risk of BE and EAC. Principal component analysis was employed to evaluate pathway-level and gene-level associations with BE/EAC, using genotypes for 270 single-nucleotide polymorphisms (SNPs) in or near 12 IGF-related genes, ascertained from 3295 BE cases, 2515 EAC cases and 3207 controls in the Barrett’s and Esophageal Adenocarcinoma Consortium (BEACON) GWAS. Gene-level signals were assessed using Multi-marker Analysis of GenoMic Annotation (MAGMA) and SNP summary statistics from BEACON and an expanded GWAS meta-analysis (6167 BE cases, 4112 EAC cases, 17 159 controls). Global variation in the IGF pathway was associated with risk of BE (P = 0.0015). Gene-level associations with BE were observed for GHR (growth hormone receptor; P = 0.00046, false discovery rate q = 0.0056) and IGF1R (IGF1 receptor; P = 0.0090, q = 0.0542). These gene-level signals remained significant at q < 0.1 when assessed using data from the largest available BE/EAC GWAS meta-analysis. No significant associations were observed for EAC. This study represents the most comprehensive evaluation to date of inherited genetic variation in the IGF pathway and BE/EAC risk, providing novel evidence that variation in two genes encoding cell-surface receptors, GHR and IGF1R, may influence risk of BE.

Objective Oesophageal adenocarcinoma (OA) incidence has risen sharply in Western countries over recent decades. Local and systemic inflammation is considered an important contributor to OA pathogenesis. Established risk factors for OA and its precursor, Barrett's oesophagus (BE), include symptomatic reflux, obesity and smoking. The role of inherited genetic susceptibility remains an area of active investigation. Here, we explore whether germline variation related to inflammatory processes influences susceptibility to BE/OA. Design We used data from a genomewide association study of 2515 OA cases, 3295 BE cases and 3207 controls. Our analysis included 7863 single-nucleotide polymorphisms (SNPs) in 449 genes assigned to five pathways: cyclooxygenase (COX), cytokine signalling, oxidative stress, human leucocyte antigen and nuclear factor-κB. A principal components-based analytic framework was employed to evaluate pathway-level and gene-level associations with disease risk. Results We identified a significant signal for the COX pathway in relation to BE risk (p=0.0059, false discovery rate q=0.03), and in gene-level analyses found an association with microsomal glutathione-S-transferase 1 (MGST1); (p=0.0005, q=0.005). Assessment of 36 MGST1 SNPs identified 14 variants associated with elevated BE risk (q<0.05). Four of these were subsequently confirmed (p<5.5×10−5) in a meta-analysis encompassing an independent set of 1851 BE cases and 3496 controls, and are known strong expression quantitative trait loci for MGST1. Three such variants were associated with similar elevations in OA risk. Conclusions This study provides the most comprehensive evaluation of inflammation-related germline variation in relation to risk of BE/OA and suggests that variants in MGST1 influence disease susceptibility.