Whiteman, David C.

ABSTRACT Bipolar disorder’s (BD) clinical heterogeneity has an unresolved genetic basis. We meta-analyzed genome-wide association studies (GWAS) of 16 BD subphenotypes in 226,032 individuals from 57 cohorts (38,022 cases); 10 advanced to multivariate and multi-trait analyses. Four factors (compulsive, psychotic, dysregulated, internalizing) explained 82.8% of shared genetic variance. BD1 and BD2 loaded on distinct factors despite a high genetic correlation; 87.0% of common-factor loci were significant in neither subtype. Unipolar mania aligned with psychosis over internalizing, and was distinguishable from BD1, and rapid cycling showed heritable cross-domain liability. We identified 356 risk loci, 158 novel, including the first univariate-GWAS associations for psychosis, unipolar mania, rapid cycling and schizoaffective disorder—and 249 credible genes (89 high-confidence), 12 with approved-drug or clinical-phase annotations. Cell-type association showed a midbrain dopaminergic–GABAergic gradient along the psychotic factor. BD’s genetic architecture appears hierarchical—a general liability resolving into dimensions of course and comorbidity beyond diagnostic subtypes.

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 Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disease with a life-time risk of 1 in 350 people and an unmet need for disease-modifying therapies. We conducted a cross-ancestry GWAS in ALS including 29,612 ALS patients and 122,656 controls which identified 15 risk loci in ALS. When combined with 8,953 whole-genome sequenced individuals (6,538 ALS patients, 2,415 controls) and the largest cortex-derived eQTL dataset (MetaBrain), analyses revealed locus-specific genetic architectures in which we prioritized genes either through rare variants, repeat expansions or regulatory effects. ALS associated risk loci were shared with multiple traits within the neurodegenerative spectrum, but with distinct enrichment patterns across brain regions and cell-types. Of the environmental and life-style risk factors obtained from literature, Mendelian randomization analyses indicated a causal role for high cholesterol levels. All ALS associated signals combined reveal a role for perturbations in vesicle mediated transport and autophagy, and provide evidence for cell-autonomous disease initiation in glutamatergic neurons.

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.

Abstract Eating disorders and substance use disorders frequently co‐occur. Twin studies reveal shared genetic variance between liabilities to eating disorders and substance use, with the strongest associations between symptoms of bulimia nervosa and problem alcohol use (genetic correlation [ r g ], twin‐based = 0.23‐0.53). We estimated the genetic correlation between eating disorder and substance use and disorder phenotypes using data from genome‐wide association studies (GWAS). Four eating disorder phenotypes (anorexia nervosa [AN], AN with binge eating, AN without binge eating, and a bulimia nervosa factor score), and eight substance‐use‐related phenotypes (drinks per week, alcohol use disorder [AUD], smoking initiation, current smoking, cigarettes per day, nicotine dependence, cannabis initiation, and cannabis use disorder) from eight studies were included. Significant genetic correlations were adjusted for variants associated with major depressive disorder and schizophrenia. Total study sample sizes per phenotype ranged from ~2400 to ~537 000 individuals. We used linkage disequilibrium score regression to calculate single nucleotide polymorphism‐based genetic correlations between eating disorder‐ and substance‐use‐related phenotypes. Significant positive genetic associations emerged between AUD and AN ( r g = 0.18; false discovery rate q = 0.0006), cannabis initiation and AN ( r g = 0.23; q < 0.0001), and cannabis initiation and AN with binge eating ( r g = 0.27; q = 0.0016). Conversely, significant negative genetic correlations were observed between three nondiagnostic smoking phenotypes (smoking initiation, current smoking, and cigarettes per day) and AN without binge eating ( r gs = −0.19 to −0.23; qs < 0.04). The genetic correlation between AUD and AN was no longer significant after co‐varying for major depressive disorder loci. The patterns of association between eating disorder‐ and substance‐use‐related phenotypes highlights the potentially complex and substance‐specific relationships among these behaviors.

42 Background: Methods of stratifying esophageal adenocarcinoma patients into prognostic groups are needed, as are new insights into genetic determinants of disease behaviour. Prognosis is likely to have non-negligible genetic influences, as mediated by host responses to tumor, resistance to therapeutic side-effects, and/or an influence on tumor development. Prior studies have used candidate-gene approaches. We took an alternative approach, using an unbiased, genome-wide approach, and novel analytic methods that may be better able to detect multi-gene interactions, which may contribute the majority of genetic effects for many clinical phenotypes. Methods: Germline DNA from a Toronto-based cohort of EAC patients (n=270) was analyzed by Omni1 Quad microarray as part of the BEAGESS initiative. Quality control and analysis was performed using PLINK, R, and GenABEL software packages. A Cox proportional hazards (CPH) model for progression-free survival tested each polymorphism for independent effects at a genome-wide significance level of P < 1E-07, adjusting for population stratification. While classical analysis has limited ability to detect gene-gene interactions, a Random Survival Forest algorithm was used to detect effects based on the complex interactions among top 1,000 polymorphisms by p-value ranking. Results: After data cleaning and standard GWAS quality control procedures, there were 735,309 SNPs and 245 patients remaining for analysis. The CPH model, adjusted for population stratification, produced a satisfactory Q-Q plot, and showed one SNP (rs7844673, Chr 8) that was significant at p=7.8E-8. In addition, Random Forest based variable selection produced a set of 20 polymorphisms that (1) reproduced 86% of the predictive ability of the full 1000 variables, and (2) also included the #3 ranked polymorphism by CPH modeling (rs9290822, Chr 3) upstream of the IGF2BP2 gene. Conclusions: A genome-wide approach has discovered two previously undescribed SNPs with a potential influence on EAC prognosis via a combination of independent and interactive effects. Validation in an independent cohort is currently being pursued.

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.

Abstract Background: While numerous susceptibility loci for epithelial ovarian cancer (EOC) have been identified, few associations have been reported with overall survival. In the absence of common prognostic genetic markers, we hypothesize that rare coding variants may be associated with overall EOC survival and assessed their contribution in two exome-based genotyping projects of the Ovarian Cancer Association Consortium (OCAC). Methods: The primary patient set (Set 1) included 14 independent EOC studies (4,293 patients) and 227,892 variants, and a secondary patient set (Set 2) included six additional EOC studies (1,744 patients) and 114,620 variants. Because power to detect rare variants individually is reduced, gene-level tests were conducted. Sets were analyzed separately at individual variants and by gene, and then combined with meta-analyses (73,203 variants and 13,163 genes overlapped). Results: No individual variant reached genome-wide statistical significance. A SNP previously implicated to be associated with EOC risk and, to a lesser extent, survival, rs8170, showed the strongest evidence of association with survival and similar effect size estimates across sets (Pmeta = 1.1E−6, HRSet1 = 1.17, HRSet2 = 1.14). Rare variants in ATG2B, an autophagy gene important for apoptosis, were significantly associated with survival after multiple testing correction (Pmeta = 1.1E−6; Pcorrected = 0.01). Conclusions: Common variant rs8170 and rare variants in ATG2B may be associated with EOC overall survival, although further study is needed. Impact: This study represents the first exome-wide association study of EOC survival to include rare variant analyses, and suggests that complementary single variant and gene-level analyses in large studies are needed to identify rare variants that warrant follow-up study. Cancer Epidemiol Biomarkers Prev; 25(3); 446–54. ©2016 AACR.

The strong male predominance in Barrett's esophagus (BE) and esophageal adenocarcinoma (EAC) remains inadequately explained, but sex hormones might be involved. We hypothesized that single nucleotide polymorphisms (SNPs) in the androgen pathway influence risk of developing BE and EAC. This genetic‐epidemiological analysis included 14 studies from Australia, Europe and North America. Polymorphisms in 16 genes coding for the androgen pathway were analyzed using a gene‐based approach: versatile gene‐based test association study. This method evaluates associations between a trait and all SNPs within a specific gene rather than each SNP marker individually as in a conventional GWAS. The data were stratified for sex, body‐mass index, waist‐to‐hip ratio, tobacco smoking and gastroesophageal reflux status. Included were data from 1,508 EAC patients, 2,383 BE patients and 2,170 control participants. SNPs within the gene CYP17A1 were associated with risk of BE in the sexes combined (p = 0.002) and in males (p = 0.003), but not in females separately (p = 0.3). This association was found in tobacco smokers (p = 0.003) and in BE patients without reflux (p = 0.004), but not in nonsmokers (p = 0.2) or those with reflux (p = 0.036). SNPs within JMJD1C were associated with risk of EAC in females (p = 0.001). However, none of these associations replicated in a subsequent sample. Fourteen other genes studied did not reach statistically significant levels of association with BE, EAC or the combination of BE and EAC, after correcting for the number of genes included in the analysis. In conclusion, genetic variants in the androgen‐related genes CYP17A1 and JMJD1C might be associated with risk of BE and EAC, respectively, but replication data with larger sample sizes are needed.

Abstract Background: A recent association study identified a common variant (rs9790517) at 4q24 to be associated with breast cancer risk. Independent association signals and potential functional variants in this locus have not been explored. Methods: We conducted a fine-mapping analysis in 55,540 breast cancer cases and 51,168 controls from the Breast Cancer Association Consortium. Results: Conditional analyses identified two independent association signals among women of European ancestry, represented by rs9790517 [conditional P = 2.51 × 10−4; OR, 1.04; 95% confidence interval (CI), 1.02–1.07] and rs77928427 (P = 1.86 × 10−4; OR, 1.04; 95% CI, 1.02–1.07). Functional annotation using data from the Encyclopedia of DNA Elements (ENCODE) project revealed two putative functional variants, rs62331150 and rs73838678 in linkage disequilibrium (LD) with rs9790517 (r2 ≥ 0.90) residing in the active promoter or enhancer, respectively, of the nearest gene, TET2. Both variants are located in DNase I hypersensitivity and transcription factor–binding sites. Using data from both The Cancer Genome Atlas (TCGA) and Molecular Taxonomy of Breast Cancer International Consortium (METABRIC), we showed that rs62331150 was associated with level of expression of TET2 in breast normal and tumor tissue. Conclusion: Our study identified two independent association signals at 4q24 in relation to breast cancer risk and suggested that observed association in this locus may be mediated through the regulation of TET2. Impact: Fine-mapping study with large sample size warranted for identification of independent loci for breast cancer risk. Cancer Epidemiol Biomarkers Prev; 24(11); 1680–91. ©2015 AACR.