Genome-wide Burden of Rare Short Deletions is Enriched in Major Depressive Disorder in Four Cohorts

Major depressive disorder (MDD) is moderately heritable, with a high prevalence and a presumed high heterogeneity. Copy number variants (CNVs) could contribute to the heritable component of risk, but the two previous genome-wide association studies of rare CNVs did not report significant findings. In this meta-analysis of four cohorts (5780 patients and 6626 control subjects), we analyzed the association of MDD to 1) genome-wide burden of rare deletions and duplications, partitioned by length (100 kb) and other characteristics, and 2) individual rare exonic CNVs and CNV regions. Patients with MDD carried significantly more short deletions than control subjects (p = .0059) but not long deletions or short or long duplications. The confidence interval for long deletions overlapped with that for short deletions, but long deletions were 70% less frequent genome-wide, reducing the power to detect increased burden. The increased burden of short deletions was primarily in intergenic regions. Short deletions in cases were also modestly enriched for high-confidence enhancer regions. No individual CNV achieved thresholds for suggestive or significant association after genome-wide correction. p values < .01 were observed for 15q11.2 duplications (TUBGCP5, CYFIP1, NIPA1, and NIPA2), deletions in or near PRKN or MSR1, and exonic duplications of ATG5. The increased burden of short deletions in patients with MDD suggests that rare CNVs increase the risk of MDD by disrupting regulatory regions. Results for longer deletions were less clear, but no large effects were observed for long multigenic CNVs (as seen in schizophrenia and autism). Further studies with larger sample sizes are warranted. [Abstract copyright: Copyright © 2019 Society of Biological Psychiatry. All rights reserved.

Classical Human Leukocyte Antigen Alleles and C4 Haplotypes Are Not Significantly Associated With Depression

Background: The prevalence of depression is higher in individuals with autoimmune diseases, but the mechanisms underlying the observed comorbidities are unknown. Shared genetic etiology is a plausible explanation for the overlap, and in this study we tested whether genetic variation in the major histocompatibility complex (MHC), which is associated with risk for autoimmune diseases, is also associated with risk for depression. Methods: We fine-mapped the classical MHC (chr6: 29.6–33.1 Mb), imputing 216 human leukocyte antigen (HLA) alleles and 4 complement component 4 (C4) haplotypes in studies from the Psychiatric Genomics Consortium Major Depressive Disorder Working Group and the UK Biobank. The total sample size was 45,149 depression cases and 86,698 controls. We tested for association between depression status and imputed MHC variants, applying both a region-wide significance threshold (3.9 × 10−6) and a candidate threshold (1.6 × 10−4). Results: No HLA alleles or C4 haplotypes were associated with depression at the region-wide threshold. HLA-B*08:01 was associated with modest protection for depression at the candidate threshold for testing in HLA genes in the meta-analysis (odds ratio = 0.98, 95% confidence interval = 0.97–0.99). Conclusions: We found no evidence that an increased risk for depression was conferred by HLA alleles, which play a major role in the genetic susceptibility to autoimmune diseases, or C4 haplotypes, which are strongly associated with schizophrenia. These results suggest that any HLA or C4 variants associated with depression either are rare or have very modest effect sizes

Genome-wide by Environment Interaction Studies of Depressive Symptoms and Psychosocial Stress in UK Biobank and Generation Scotland

Stress is associated with poorer physical and mental health. To improve our understanding of this link, we performed genome-wide association studies (GWAS) of depressive symptoms and genome-wide by environment interaction studies (GWEIS) of depressive symptoms and stressful life events (SLE) in two UK population-based cohorts (Generation Scotland and UK Biobank). No SNP was individually significant in either GWAS, but gene-based tests identified six genes associated with depressive symptoms in UK Biobank (DCC, ACSS3, DRD2, STAG1, FOXP2 and KYNU; p < 2.77 x 10(-6)). Two SNPs with genome-wide significant GxE effects were identified by GWEIS in Generation Scotland: rs12789145 (53-kb downstream PIWIL4; p = 4.95 x 10(-9); total SLE) and rs17070072 (intronic to ZCCHC2; p = 1.46 x 10(-8); dependent SLE). A third locus upstream CYLC2 (rs12000047 and rs12005200, p < 2.00 x 10(-8); dependent SLE) when the joint effect of the SNP main and GxE effects was considered. GWEIS gene-based tests identified: MTNR1B with GxE effect with dependent SLE in Generation Scotland; and PHF2 with the joint effect in UK Biobank (p < 2.77 x 10(-6)). Polygenic risk scores (PRSs) analyses incorporating GxE effects improved the prediction of depressive symptom scores, when using weights derived from either the UK Biobank GWAS of depressive symptoms (p = 0.01) or the PGC GWAS of major depressive disorder (p = 5.91 x 10(-3)). Using an independent sample, PRS derived using GWEIS GxE effects provided evidence of shared aetiologies between depressive symptoms and schizotypal personality, heart disease and COPD. Further such studies are required and may result in improved treatments for depression and other stress-related conditions

Integrated analysis of environmental and genetic influences on cord blood DNA methylation in new-borns

Epigenetic processes, including DNA methylation (DNAm), are among the mechanisms allowing integration of genetic and environmental factors to shape cellular function. While many studies have investigated either environmental or genetic contributions to DNAm, few have assessed their integrated effects. Here we examine the relative contributions of prenatal environmental factors and genotype on DNA methylation in neonatal blood at variably methylated regions (VMRs) in 4 independent cohorts (overall n = 2365). We use Akaike’s information criterion to test which factors best explain variability of methylation in the cohort-specific VMRs: several prenatal environmental factors (E), genotypes in cis (G), or their additive (G + E) or interaction (GxE) effects. Genetic and environmental factors in combination best explain DNAm at the majority of VMRs. The CpGs best explained by either G, G + E or GxE are functionally distinct. The enrichment of genetic variants from GxE models in GWAS for complex disorders supports their importance for disease risk

Candidate genes in panic disorder: Meta-analyses of 23 common variants in major anxiogenic pathways

The utilization of molecular genetics approaches in examination of panic disorder (PD) has implicated several variants as potential susceptibility factors for panicogenesis. However, the identification of robust PD susceptibility genes has been complicated by phenotypic diversity, underpowered association studies and ancestry-specific effects. In the present study, we performed a succinct review of case–control association studies published prior to April 2015. Meta-analyses were performed for candidate gene variants examined in at least three studies using the Cochrane Mantel–Haenszel fixed-effect model. Secondary analyses were also performed to assess the influences of sex, agoraphobia co-morbidity and ancestry-specific effects on panicogenesis. Meta-analyses were performed on 23 variants in 20 PD candidate genes. Significant associations after correction for multiple testing were observed for three variants, TMEM132D rs7370927 (T allele: odds ratio (OR)=1.27, 95% confidence interval (CI): 1.15–1.40, P=2.49 × 10−6), rs11060369 (CC genotype: OR=0.65, 95% CI: 0.53–0.79, P=1.81 × 10−5) and COMT rs4680 (Val (G) allele: OR=1.27, 95% CI: 1.14–1.42, P=2.49 × 10−5) in studies with samples of European ancestry. Nominal associations that did not survive correction for multiple testing were observed for NPSR1 rs324891 (T allele: OR=1.22, 95% CI: 1.07–1.38, P=0.002), TPH1 rs1800532 (AA genotype: OR=1.46, 95% CI: 1.14–1.89, P=0.003) and HTR2A rs6313 (T allele: OR=1.19, 95% CI: 1.07–1.33, P=0.002) in studies with samples of European ancestry and for MAOA-uVNTR in female PD (low-active alleles: OR=1.21, 95% CI: 1.07–1.38, P=0.004). No significant associations were observed in the secondary analyses considering sex, agoraphobia co-morbidity and studies with samples of Asian ancestry. Although these findings highlight a few associations, PD likely involves genetic variation in a multitude of biological pathways that is diverse among populations. Future studies must incorporate larger sample sizes and genome-wide approaches to further quantify the observed genetic variation among populations and subphenotypes of PD