Multimodal Associations of FKBP5 Methylation with Emotion-Regulatory Brain Circuits

BACKGROUND

Understanding the biological processes underlying individual differences in emotion regulation and stress responsivity is a key challenge for translational neuroscience. The gene FKBP5 is a core regulator in molecular stress signaling that is implicated in the development of psychiatric disorders. Yet it remains unclear how FKBP5 DNA methylation (DNAm) in peripheral blood relates to individual differences in measures of neural structure and function, and their relevance to daily-life stress responsivity.

METHODS

Here, we characterize multimodal correlates of FKBP5 DNAm by combining epigenetic data with neuroimaging and Ambulatory Assessment in a sample of 395 healthy individuals.

RESULTS

First, we show that FKBP5 demethylation as a psychiatric risk factor relates to an anxiety-associated reduction of gray matter volume in the ventromedial prefrontal cortex (vmPFC), a brain area that is involved in emotion regulation and mental health risk and resilience. This effect of epigenetic upregulation of FKBP5 on neuronal structure is more pronounced where FKBP5 is epigenetically downregulated at baseline. Leveraging 208 functional MRI scans during a well-established emotion processing task we find that FKBP5 DNAm in peripheral blood is associated with functional difference of prefrontal-limbic circuits modulating affective responsivity to daily stressors, which we measured using ecological momentary assessment in daily life.

CONCLUSIONS

Overall, we demonstrate how FKBP5 contributes to interindividual differences in neural and real-life affect regulation via structural and functional changes in prefrontal-limbic brain circuits.

Genome-wide association analysis of chronic lymphocytic leukaemia, Hodgkin lymphoma and multiple myeloma identifies pleiotropic risk loci

B-cell malignancies (BCM) originate from the same cell of origin, but at different maturation stages and have distinct clinical phenotypes. Although genetic risk variants for individual BCMs have been identified, an agnostic, genome-wide search for shared genetic susceptibility has not been performed. We explored genome-wide association studies of chronic lymphocytic leukaemia (CLL, N = 1,842), Hodgkin lymphoma (HL, N = 1,465) and multiple myeloma (MM, N = 3,790). We identified a novel pleiotropic risk locus at 3q22.2 (NCK1, rs11715604, P = 1.60 × 10-9) with opposing effects between CLL (P = 1.97 × 10-8) and HL (P = 3.31 × 10-3). Eight established non-HLA risk loci showed pleiotropic associations. Within the HLA region, Ser37 + Phe37 in HLA-DRB1 (P = 1.84 × 10-12) was associated with increased CLL and HL risk (P = 4.68 × 10-12), and reduced MM risk (P = 1.12 × 10-2), and Gly70 in HLA-DQB1 (P = 3.15 × 10-10) showed opposing effects between CLL (P = 3.52 × 10-3) and HL (P = 3.41 × 10-9). By integrating eQTL, Hi-C and ChIP-seq data, we show that the pleiotropic risk loci are enriched for B-cell regulatory elements, as well as an over-representation of binding of key B-cell transcription factors. These data identify shared biological pathways influencing the development of CLL, HL and MM. The identification of these risk loci furthers our understanding of the aetiological basis of BCMs.

Genetic factors influencing the risk of multiple myeloma bone disease

A major complication of multiple myeloma (MM) is the development of osteolytic lesions, fractures and bone pain. To identify genetic variants influencing the development of MM bone disease (MBD), we analyzed MM patients of European ancestry (totaling 3774), which had been radiologically surveyed for MBD. Each patient had been genotyped for ~6 00 000 single-nucleotide polymorphisms with genotypes for six million common variants imputed using 1000 Genomes Project and UK10K as reference. We identified a locus at 8q24.12 for MBD (rs4407910, OPG/TNFRSF11B, odds ratio=1.38, P=4.09 × 10(-9)) and a promising association at 19q13.43 (rs74676832, odds ratio=1.97, P=9.33 × 10(-7)). Our findings demonstrate that germline variation influences MBD and highlights the importance of RANK/RANKL/OPG pathway in MBD development. These findings will contribute to the development of future strategies for prevention of MBD in the early precancerous phases of MM.

Multi-ancestry genome-wide association study of 21,000 cases and 95,000 controls identifies new risk loci for atopic dermatitis

Genetic association studies have identified 21 loci associated with atopic dermatitis risk predominantly in populations of European ancestry. To identify further susceptibility loci for this common complex skin disease, we performed a meta-analysis of >15 million genetic variants in 21,399 cases and 95,464 controls from populations of European, African, Japanese and Latino ancestry, followed by replication in 32,059 cases and 228,628 controls from 18 studies. We identified 10 novel risk loci, bringing the total number of known atopic dermatitis risk loci to 31 (with novel secondary signals at 4 of these). Notably, the new loci include candidate genes with roles in regulation of innate host defenses and T-cell function, underscoring the important contribution of (auto-)immune mechanisms to atopic dermatitis pathogenesis.

The association between lower educational attainment and depression owing to shared genetic effects? Results in ~25,000 subjects

An association between lower educational attainment (EA) and an increased risk for depression has been confirmed in various western countries. This study examines whether pleiotropic genetic effects contribute to this association. Therefore, data were analyzed from a total of 9662 major depressive disorder (MDD) cases and 14,949 controls (with no lifetime MDD diagnosis) from the Psychiatric Genomics Consortium with additional Dutch and Estonian data. The association of EA and MDD was assessed with logistic regression in 15,138 individuals indicating a significantly negative association in our sample with an odds ratio for MDD 0.78 (0.75-0.82) per standard deviation increase in EA. With data of 884,105 autosomal common single-nucleotide polymorphisms (SNPs), three methods were applied to test for pleiotropy between MDD and EA: (i) genetic profile risk scores (GPRS) derived from training data for EA (independent meta-analysis on ~120,000 subjects) and MDD (using a 10-fold leave-one-out procedure in the current sample), (ii) bivariate genomic-relationship-matrix restricted maximum likelihood (GREML) and (iii) SNP effect concordance analysis (SECA). With these methods, we found (i) that the EA-GPRS did not predict MDD status, and MDD-GPRS did not predict EA, (ii) a weak negative genetic correlation with bivariate GREML analyses, but this correlation was not consistently significant, (iii) no evidence for concordance of MDD and EA SNP effects with SECA analysis. To conclude, our study confirms an association of lower EA and MDD risk, but this association was not because of measurable pleiotropic genetic effects, which suggests that environmental factors could be involved, for example, socioeconomic status.

Novel genetic matching methods for handling population stratification in genome-wide association studies

Background

A usually confronted problem in association studies is the occurrence of population stratification. In this work, we propose a novel framework to consider population matchings in the contexts of genome-wide and sequencing association studies. We employ pairwise and groupwise optimal case-control matchings and present an agglomerative hierarchical clustering, both based on a genetic similarity score matrix. In order to ensure that the resulting matches obtained from the matching algorithm capture correctly the population structure, we propose and discuss two stratum validation methods. We also invent a decisive extension to the Cochran-Armitage Trend test to explicitly take into account the particular population structure.

Results

We assess our framework by simulations of genotype data under the null hypothesis, to affirm that it correctly controls for the type-1 error rate. By a power study we evaluate that structured association testing using our framework displays reasonable power. We compare our result with those obtained from a logistic regression model with principal component covariates. Using the principal components approaches we also find a possible false-positive association to Alzheimer’s disease, which is neither supported by our new methods, nor by the results of a most recent large meta analysis or by a mixed model approach.

Conclusions

Matching methods provide an alternative handling of confounding due to population stratification for statistical tests for which covariates are hard to model. As a benchmark, we show that our matching framework performs equally well to state of the art models on common variants.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-015-0521-4) contains supplementary material, which is available to authorized users.

Rare variant testing of imputed data: an analysis pipeline typified

Important methodological advancements in rare variant association testing have been made recently, among them collapsing tests, kernel methods and the variable threshold (VT) technique. Typically, rare variants from a region of interest are tested for association as a group ('bin'). Rare variant studies are already routinely performed as whole-exome sequencing studies. As an alternative approach, we propose a pipeline for rare variant analysis of imputed data and develop respective quality control criteria. We provide suggestions for the choice and construction of analysis bins in whole-genome application and support the analysis with implementations of standard burden tests (COLL, CMAT) in our INTERSNP-RARE software. In addition, three rare variant regression tests (REG, FRACREG and COLLREG) are implemented. All tests are accompanied with the VT approach which optimizes the definition of 'rareness'. We integrate kernel tests as implemented in SKAT/SKAT-O into the suggested strategies. Then, we apply our analysis scheme to a genome-wide association study of Alzheimer's disease. Further, we show that our pipeline leads to valid significance testing procedures with controlled type I error rates. Strong association signals surrounding the known APOE locus demonstrate statistical power. In addition, we highlight several suggestive rare variant association findings for follow-up studies, including genomic regions overlapping MCPH1, MED18 and NOTCH3. In summary, we describe and support a straightforward and cost-efficient rare variant analysis pipeline for imputed data and demonstrate its feasibility and validity. The strategy can complement rare variant studies with next generation sequencing data.

Epistatic interaction of genetic depression risk variants in the human subgenual cingulate cortex during memory encoding

Recent genome-wide association studies have pointed to single-nucleotide polymorphisms (SNPs) in genes encoding the neuronal calcium channel CaV1.2 (CACNA1C; rs1006737) and the presynaptic active zone protein Piccolo (PCLO; rs2522833) as risk factors for affective disorders, particularly major depression. Previous neuroimaging studies of depression-related endophenotypes have highlighted the role of the subgenual cingulate cortex (CG25) in negative mood and depressive psychopathology. Here, we aimed to assess how recently associated PCLO and CACNA1C depression risk alleles jointly affect memory-related CG25 activity as an intermediate phenotype in clinically healthy humans. To investigate the combined effects of rs1006737 and rs2522833 on the CG25 response, we conducted three functional magnetic resonance imaging studies of episodic memory formation in three independent cohorts (N=79, 300, 113). An epistatic interaction of PCLO and CACNA1C risk alleles in CG25 during memory encoding was observed in all groups, with carriers of no risk allele and of both risk alleles showing higher CG25 activation during encoding when compared with carriers of only one risk allele. Moreover, PCLO risk allele carriers showed lower memory performance and reduced encoding-related hippocampal activation. In summary, our results point to region-specific epistatic effects of PCLO and CACNA1C risk variants in CG25, potentially related to episodic memory. Our data further suggest that genetic risk factors on the SNP level do not necessarily have additive effects but may show complex interactions. Such epistatic interactions might contribute to the 'missing heritability' of complex phenotypes.

Confirming genes influencing risk to cleft lip with/without cleft palate in a case-parent trio study

A collection of 1,108 case-parent trios ascertained through an isolated, nonsyndromic cleft lip with or without cleft palate (CL/P) was used to replicate the findings from a genome-wide association study (GWAS) conducted by Beaty et al. (Nat Genet 42:525-529, 2010), where four different genes/regions were identified as influencing risk to CL/P. Tagging SNPs for 33 different genes were genotyped (1,269 SNPs). All four of the genes originally identified as showing genome-wide significance (IRF6, ABCA4 and MAF, plus the 8q24 region) were confirmed in this independent sample of trios (who were primarily of European and Southeast Asian ancestry). In addition, eight genes classified as 'second tier' hits in the original study (PAX7, THADA, COL8A1/FILIP1L, DCAF4L2, GADD45G, NTN1, RBFOX3 and FOXE1) showed evidence of linkage and association in this replication sample. Meta-analysis between the original GWAS trios and these replication trios showed PAX7, COL8A1/FILIP1L and NTN1 achieved genome-wide significance. Tests for gene-environment interaction between these 33 genes and maternal smoking found evidence for interaction with two additional genes: GRID2 and ELAVL2 among European mothers (who had a higher rate of smoking than Asian mothers). Formal tests for gene-gene interaction (epistasis) failed to show evidence of statistical interaction in any simple fashion. This study confirms that many different genes influence risk to CL/P.

Increased levels of glucocorticoid receptors and enhanced glucocorticoid receptor auto-regulation after hydrocortisone challenge in B-lymphoblastoids from patients with affective disorders

The stress response is mediated by a negative feedback effect of glucocorticoids on corticosteroid receptors. Here, we examine the potential contribution of these receptors and their response to a glucocorticoid challenge to dysfunctions of the hypothalamic-pituitary-adrenal axis reported for patients with affective disorders. In a pilot-study, we established B-lymphoblastoid cell lines from patients suffering from affective disorders and healthy subjects and measured the quantity of glucocorticoid receptors at steady state conditions after 12-weeks cell culture. After short-term incubation with 0.1 microM hydrocortisone for 48 h, the decrease of glucocorticoid receptors was also investigated. After 12-weeks cell culture, we found a significantly higher number of cytosolic glucocorticoid receptors in B-lymphoblastoids from patients (B(max)=804.9+/-342.5 fmol/mg protein) compared to those from healthy subjects (B(max)=576.9+/-190.3 fmol/mg protein: p=0.045; t-test). The increase of the glucocorticoid receptor level in the group of patients could be attributed largely to the higher number of these receptors measured in B-lymphoblastoids of patients suffering from major depressive disorder. The in vitro regulation of glucocorticoid receptors in response to 0.1 microM hydrocortisone for 48 h resulted in a significantly larger decrease in cultures of B-lymphoblastoids derived from patients (to 32.9+/-7.5%) than in those from healthy subjects (to 45.8+/-8.2%). The stronger decrease of glucocorticoid receptors in the group of patients (p=0.0001; t-test) was independent of the duration of illness and medication, suggesting a trait-like characteristic of the response.

Excess of allele1 for alpha3 subunit GABA receptor gene (GABRA3) in bipolar patients: a multicentric association study

The available data from preclinical and pharmacological studies on the role of gamma amino butyric acid (GABA) support the hypothesis that a dysfunction in brain GABAergic system activity contributes to the vulnerability to bipolar affective disorders (BPAD). Moreover, the localization of the alpha3 subunit GABA receptor GABRA3 gene on the Xq28, a region of interest in certain forms of bipolar illness, suggests that GABRA3 may be a candidate gene in BPAD. In the present study, we tested the genetic contribution of the GABRA3 dinucleotide polymorphism in a European multicentric case-control sample, matched for sex and ethnogeographical origin. Allele and genotype (in females) frequencies were compared in 185 BPAD patients and 370 controls. A significant increase of genotype 1-1 was observed in BPAD females compared to controls (P=0.0004). Furthermore, when considering recessivity of allele 1 (females with genotype 1-1 and males carrying allele 1), results were even more significant (P= 0.00002). Our findings suggest that the GABRA3 polymorphism may confer susceptibility to or may be in linkage disequilibrium with another gene involved in the genetic etiology of BPAD.