Abstract
The stress-related gene FKBP5 has been related to dysregulated glucocorticoid receptor (GR) signaling, showing increased GR sensitivity in trauma-exposed subjects with post-traumatic stress disorder (PTSD) but not in those without PTSD. However, the neural mechanism underlying the effects of FKBP5 remains poorly understood. Two hundred and thirty-seven Han Chinese adults who had lost their only child were included. Four FKBP5 single nucleotide polymorphisms (rs3800373, rs9296158, rs1360780, and rs9470080) were genotyped. All 179 participants were successfully divided into three FKBP5 diplotype subgroups according to two major FKBP5 H1 and H2 yin yang haplotypes. Brain average spectral power was compared using a two-way (PTSD diagnosis and FKBP5 diplotypes) analysis of covariance within four separate frequency bands (slow-5, slow-4, slow-3, and slow-2). Adults with PTSD showed lower spectral power in bilateral parietal lobules in slow-4 and in left inferior frontal gyrus (IFG) in slow-5. There was significant FKBP5 diplotype main effect in anterior cingulate cortex (ACC) in slow-4 (H1/H1 higher than other two subgroups), and in precentral/postcentral gyri and middle cingulate cortex (MCC) in slow-3 (H2/H2 higher than other two subgroups). Also, there was a significant diagnosis × FKBP5 diplotype interaction effect in right parietal lobule in slow-3. These findings suggest that adults with PTSD have lower low-frequency power in executive control network regions. Lower power in ACC and greater power in the motor/sensory areas in FKBP5 high-risk diplotype group suggest a disturbance of emotional processing and hypervigilance/sensitization to threatening stimuli. The interaction effect of diagnosis × FKBP5 in parietal lobule may contribute to PTSD development.
Similar content being viewed by others
Introduction
Post-traumatic stress disorder (PTSD) is characterized by dysregulation of the hypothalamus-pituitary-adrenal (HPA) axis after exposure to traumatic events, displayed as enhanced negative-feedback inhibition of the HPA axis1,2. One typical clinical feature of this abnormality in PTSD is an exaggerated suppression of cortisol in response to dexamethasone administration and the enhanced sensitivity of glucocorticoid receptors (GR)3,4. Even so, only a minority of trauma-exposed individuals will eventually develop PTSD5, and the increased GR sensitivity following stress happens only in PTSD but not in non-PTSD subjects3,4. The underlying mechanisms of the association between PTSD development and GR sensitivity remain poorly understood.
Genetic factors that regulate GR signaling influence the individual differences in the HPA axis response to stress6,7. A critical regulatory gene is the FK506-binding protein 5 (FKBP5) gene which encodes the FKBP5 protein8. In the absence of cortisol, FKBP5 binds to the GR complex, resulting in decreased affinity for cortisol and less efficient nuclear translocation of GR. In the presence of cortisol, FKBP5 is exchanged with other co-chaperones, causing more efficient translocation of GR to the nucleus9. Healthy carriers of specific FKBP5 single nucleotide polymorphisms (SNPs) were tended to have higher FKBP5 protein expression and thus decreased cortisol affinity and nuclear translocation of GR, resulting in GR resistance6, suggesting the presence of protective features after response to stress. However, this functional association is switched in PTSD6, in which the same FKBP5 SNPs were associated with increased GR sensitivity, conferring “risk” SNPs for PTSD. So far, little is known about the exact mechanisms of different modulatory effects of FKBP5 risk SNPs in people who experienced traumatic events with and without PTSD.
Imaging genomics is significantly improving our understanding of human brain function, and offers a means to study molecular and cellular abnormalities that are implicated in many psychiatric disorders in which genes play an important role10,11,12,13. To the best of our knowledge, only one imaging genomics study with a small sample size (a total of 54 subjects) has investigated the effect of FKBP5 risk SNP polymorphisms on brain function of PTSD14. In that study, Fani et al. firstly reported that PTSD patients who carried two FKBP5 rs1360780 risk SNPs had lower fractional anisotropy in the cingulum, compared to other diagnostic and genotype groups. Also, trauma-exposed subjects who carried two FKBP5 rs1360780 risk SNPs had lower functional connectivity of hippocampus-anterior cingulate cortex than non-risk SNP carriers, regardless of PTSD diagnosis14. The FKBP5 rs1360780 has been demonstrated to function together with other FKBP5 SNPs which are in strong linkage disequilibrium, resulting in gene dose-dependent regulation on the HPA axis6. Several studies have focused on the two major functional FKBP5 haplotypes driven from four SNPs (rs3800373, rs9296158, rs1360780, and rs9470080)—the H1 (carrying zero risk SNPs: AGCC) and H2 (carrying four risk SNPs: CATT) yin yang haplotypes15,16. However, no imaging genomics study to date has used this desirable haplotype to investigate the modulatory effects of FKBP5 on brain function in PTSD development. In a recent fMRI study examining effects of FKBP5 risk SNPs, Bryant et al. divided the healthy participants into high-risk and low-risk genogroups based on the number of risk alleles carried (using the above-mentioned four FKBP5 SNPs) and then investigated the differences in resting-state brain function between these two genogroups17. They found that the high-risk genogroup showed less power in low frequency range but more power in higher frequency range in the frontotemporoparietal network, compared to the low-risk genogroup.
In this prospective study, we aimed to examine the modulatory effects of FKBP5 H1 and H2 yin yang haplotypes on resting state brain function using resting-state fMRI in trauma-exposed subjects with and without PTSD. A wide variation in the types of trauma and differences in race/ethnicity of the populations studied might potentially influence the findings of imaging genomics analyses18. Therefore, we conducted a PTSD survey in a specific group of Han Chinese adults who had lost their only child. By doing so, we restricted our study to participants of the same ethnicity (Han Chinese) and homogenous traumatic event (loss of only child). Power spectrum analysis19,20—an important fMRI approach that measures the strength of intrinsic brain activity—was used in this study for three reasons. Firstly, the power spectrum, especially within the low-frequency band, is a physiologically meaningful and powerful biomarker in exploring human spontaneous brain activity20,21,22, and has been successfully applied in prior genetic association studies in several mental and neuropsychiatric disorders such as mild cognitive impairment23,24 and major depressive disorder25. Secondly, in a prior imaging genomics study using four FKBP5 risk SNPs, only spectral power demonstrated significant differences between the high-risk and low-risk genogroups17. Thirdly, the power spectrum analysis is a whole-brain algorithm that needs no region-specific hypothesis, which might be more suitable for a preliminary study. We hypothesized that the FKBP5 haplotypes would exhibit different modulatory effects on brain spectral power in Han Chinese adults who had lost their only child with and without PTSD.
Methods
Subjects
The present study was approved by the Medical Research Ethics Committee of Jiangsu University. Each participant provided written informed consent. Between September 2016 and March 2017, we conducted a survey in Jiangsu Province, China on a population of the Han Chinese adults who had lost their only child during the period of “One-Child Policy”26,27. All 237 Han adults who had lost their only child—without other major traumatic exposures based on the clinician-administered PTSD scale (CAPS) life events checklist—were successfully interviewed and screened by the clinician-administered PTSD scale (CAPS). They were also further screened with the Chinese version of the structured clinical interview for DSM-IV (SCID)28 which was revised by Prof. Lipeng Fei from the Beijing Hui Long Guan Hospital. After doing this, 170 out of the total of 237 adults did not meet any diagnostic criteria for mental illness (including current or lifetime PTSD) or substance use disorders; 57 adults were diagnosed with PTSD; the remaining 10 adults were diagnosed with other psychiatric disorders (5 with major depressive disorder, 4 with generalized anxiety disorder, and 1 with both depression and anxiety diagnosis) and they were not included in the current study.
Exclusion criteria for the following MRI study were as follows: any current or history of brain injury or other major medical or neurological conditions (4 adults without PTSD were excluded for cerebral infarction or ischemia, and 1 adult without PTSD was ruled out for a history of major depressive disorder and antidepressant drug therapy); any MRI contraindication (none); and left-handedness (none).
MR scanning
MR imaging was performed with a 3-T scanner (Achieva 3.0 TTX; Philips, Amsterdam, the Netherlands). Each participant was instructed to stay still in the scanner during image acquisition, keep eyes closed, remain awake and move as little as possible. T1-weighted anatomical images were acquired using a three-dimensional turbo fast echo sequence (repetition time ms/echo time msec [TR/TE], 9.7/4.6; flip angle, 9°; field of view (FOV), 256 × 256 mm2; matrix size, 256 × 256; slice thickness, 1 mm; 160 sagittal slices). For resting-state functional imaging, we performed a single-shot, gradient-recalled echo-planar imaging sequence (TR/TE, 2000/30; flip angle, 90°; FOV = 192 × 192 mm2; matrix, 64 × 64; voxel size, 3 × 3 × 4 mm3; volume number, 230; 35 axial slices). The total scan time for resting-state functional imaging was 460 s.
Measures
All participants were assessed with neuropsychological tests, which included: the Hamilton Depression (HAMD)29 and Hamilton Anxiety (HAMA)30 rating scales, Mini-Mental State Examination (MMSE)31, Chinese Social Support Rating Scale (SSRS) with sections about subjective support, objective support and the utility of support32, and individual Simple Coping Style Questionnaire (SCSQ) with active and negative coping score, and the coping tendency score33. A detailed description is available in the online Supplementary Note 1.
Data preprocessing
MRI data was preprocessed with the Data Processing Assistant for Resting-State fMRI (DPARSF, http://rfmri.org/DPARSF) which is based on Statistical Parametric Mapping (http://www.fil.ion.ucl.ac.uk/spm). First, the initial ten volumes were excluded. Then, the slice timing and head motion correction were conducted on all the remaining volumes. Individual T1-weighted images were co-registered to the functional images and then segmented into gray matter, white matter, and cerebrospinal fluid, and transformed into the standard Montreal Neurological Institute (MNI) space using the Diffeomorphic Anatomical Registration Through Exponentiated Lie algebra (DARTEL) method34. Finally, the functional images were transformed into the MNI stereotaxic space (3 × 3 × 3 mm3), using the parameters of the T1-weighted image normalization, and then smoothed with an 8 mm full width at half maximum (FWHM) isotropic Gaussian kernel.
Quality control and nuisance regression
To minimize head motion confounds, we used the Friston 24-parameter model35 to regress out head motion effects. Individual head translations, rotations, and framewise displacement (using the Jenkinson formula) were calculated. Seven adults were excluded (5 without PTSD, and 2 with PTSD) for head translations > 1.5 mm or rotations > 1.5°, and 2 adults (without PTSD) were excluded for mean framewise displacement > 2.5 standard deviations. The mean framewise displacement of each participant was also included as a nuisance covariate in the statistical analysis of the fMRI data36,37. Mean signals from cerebrospinal fluid and white matter were also regressed out as spurious variance to restrict the analysis to gray matter.
After performing the quality control, a total of nine participants were excluded, and 55 adults with PTSD and 158 adults without PTSD remained.
Power spectrum analysis
Firstly, for a given voxel, the time course was converted to the frequency domain using a fast Fourier transform20,22. Then, the square root of the power spectrum was computed and averaged across a predefined frequency interval. Based on the most current knowledge, some brain disorders38,39, as well as some genes17, may selectively affect the brain power spectrum within a certain frequency band. Here, we have divided the full frequency range (0–0.25 Hz) into four narrowly-defined bands according to prior studies21,22: slow-5 (0.01–0.027 Hz), slow-4 (0.027–0.073 Hz), slow-3 (0.073–0.198 Hz), and slow-2 (0.198–0.25 Hz); and we calculated the average power spectrum in individual frequency band. For standardization purposes, the average power spectrum of each voxel was divided by the global mean values.
DNA genotyping
With the exception of three adults without PTSD who refused the blood collection procedure, we successfully collected DNA data from all other participants from peripheral blood samples. Four FKBP5 SNPs (rs3800373, rs9296158, rs1360780, and rs9470080) were genotyped using the Improved Multiple Ligase Detection Reaction (iMLDR) technique developed by Genesky Biotechnologies, Inc. (Shanghai, China)40. Detailed primers for these four SNPs of FKBP5 are listed in Supplementary Table 1. About 5% of the samples were randomly selected for confirmation, and the results were 100% concordant. For all 210 study participants with available gene data—and who passed quality control for their fMRI data—the distributions of rs3800373, rs9296158, rs1360780, and rs9470080 did not differ from Hardy-Weinberg equilibrium (P = 0.74, 0.45, 0.62, and 0.30, respectively) calculated using R, version 3.5.3 (https://www.r-project.org).
FKBP5 haplotype estimation
A haplotype is a combination of alleles that are located close together on the same chromosome and inherited together41. Here, we used Haploview version 4.242 and genotype data to compute the linkage disequilibrium blocks for these four FKBP5 SNPs, and found they were in strong linkage disequilibrium (Fig. 1) and in approximate allelic identity (Supplementary Table 2). There were 9 FKBP5 haplotypes (frequency ≥ 0.01) identified in this study, with two major functional haplotypes—the H1 and H2 yin yang haplotypes (Fig. 1). The H1 (yin) haplotype carries zero risk SNPs (AGCC), while the H2 (yang) haplotype carries all four risk SNPs (CATT)6,15,16. These two major functional haplotypes accounted for 92.1% of haplotype diversity in this study, which is similar to findings in prior studies15,16. Then, PHASE version 2.143,44,45 was performed to determine the most probable haplotype assignments for each individual by assessing the probability of each possible haplotype and determining a confidence score. A diplotype represents a pair of haplotypes on homologous chromosomes and could provide more complete genetic information46,47. For this reason, the FKBP5 diplotype was further estimated for 210 study participants; 179 of them carried homozygote and heterozygote combinations of the H1 and H2 yin yang haplotypes and were included in the final diplotype based neuroimaging analysis. From these 179 adults, 87 were H1/H1 (64 without PTSD, 23 with PTSD), 72 were H1/H2 (51 without PTSD, 21 with PTSD), and 20 were H2/H2 (15 without PTSD, 5 with PTSD) (Supplementary Table 3).
Statistical analysis
SPSS version 25 (IBM Corp, Armonk, New York, USA) was used to analyze the demographic and neuropsychological data. A two-way (diagnosis of PTSD, FKBP5 diplotypes) analysis of variance was used to evaluate the main effects of diagnosis, FKBP5 diplotypes and their interaction effects on demographic and neuropsychological data. A voxel-wise two-way (diagnosis, FKBP5 diplotypes) full factorial analysis of covariance (ANCOVA) was performed using SPM12 to assess the main effect of PTSD diagnosis, FKBP5 diplotypes and their interaction effects on brain power spectrum maps within the four different frequency bands, adjusting for effects of age, sex, educational level, duration since child-loss trauma, and framewise displacement due to head motion. The frequency band was not considered as a repeated factor in this paper as we wanted to primarily focus on the effects that diagnosis and genetic variation may exert on the power spectrum within different bands, but not the differences among different bands themselves. Results were corrected for multiple comparisons (corrected P < 0.05) using random-field theory (RFT)48 with initial voxel level P < 0.001 and cluster level P < 0.05.
A partial correlation analysis was performed to examine the relationship between regions with significantly different power spectra and CAPS, SSRS, SCSQ, HAMA, and HAMD, with the inclusion of age, sex, educational level, duration since child-loss trauma, and head motion as covariates. Correlation results were corrected for multiple comparisons using the Bonferroni correction for the number of regions where altered power spectra were detected from the two-way ANCOVA (cut-off P values of 0.05/8 = 0.006, corresponding to all eight regions showing differences in this study).
Results
Clinical and neuropsychological data
A flowchart describing the study population is detailed in Supplementary Fig. 1. All 179 adults (49 with PTSD and 130 without PTSD) who carried homozygote and heterozygote combinations of the H1 (yin) and H2 (yang) haplotypes were included in the final FKBP5 diplotypes based neuroimaging analyses (Table 1 and Supplementary Table 3).
There were no significant differences between PTSD and non-PTSD groups in age, educational level, SSRS, or SCSQ (all P > 0.05), but the PTSD group showed higher CAPS, HAMA and HAMD scores, higher male/female ratio, and shorter duration since losing the child (Table 1) than the non-PTSD groups. There was no significant diagnosis or FKBP5 diplotypes main effect, and no interaction effect on other clinical and neuropsychological data.
The main effect of PTSD diagnosis
Significant PTSD diagnosis effects on spectral power were found in bilateral parietal lobules in slow-4 band and in left inferior frontal gyrus (IFG) in slow-5 band. Post-hoc analysis showed that PTSD adults had lower spectral power in these regions, relative to trauma-exposed adults without PTSD (Fig. 2 and Supplementary Table 4).
The main effect of FKBP5 diplotypes
There was significant FKBP5 diplotypes main effect on spectral power in the anterior cingulate cortex (ACC) in slow-4 band, and in the bilateral precentral/postcentral gyri and middle cingulate cortex (MCC) in slow-3 band in all participants, irrespective of PTSD diagnosis (Fig. 3 and Supplementary Table 4). At ACC, the H1/H1 diplotype subgroup had higher spectral power than other two genogroups. For the precentral/postcentral gyri and MCC, the H2/H2 diplotype group had higher spectral power than other two genogroups. The FKBP5 diplotype modulation accounted for 14% of the variance in spectral power in left precentral/postcentral gyri, 10% in right precentral/postcentral gyri, 10% in the right MCC and 9% in the ACC (partial eta squared [η2] = 14, 10, 10, and 9%, respectively).
Interaction of FKBP5 diplotypes and PTSD diagnosis
There was a significant diagnosis × FKBP5 diplotype interaction effect on spectral power in right parietal lobule in slow-3 band (Fig. 4 and Supplementary Table 4). Post-hoc analysis showed that within non-PTSD adults, the H2/H2 diplotype subgroup had higher spectral power than other two genogroups; whereas within PTSD adults, the H2/H2 diplotype subgroup had lower spectral power than the other two genogroups (Fig. 4). This diagnosis × FKBP5 diplotype interaction accounted for 10% of the variance in spectral power in right parietal lobule (η2 = 10%).
Partial correlation results
A marginally negative partial correlation was found between the power spectrum in right parietal lobule in slow-4 band and CAPS scores, albeit only in the PTSD group (r = −0.27; P = 0.06, Supplementary Fig. 2). Power spectra in brain regions affected by FKBP5 diplotypes and diagnosis × FKBP5 interaction were not correlated with any clinical or neuropsychological indices.
Additional analyses of structural data
To evaluate whether the functional results in this study might be confounded by brain structural changes, we further performed voxel-based morphometry (VBM) to examine possible structural deficits, using the CAT12 Toolbox (http://dbm.neuro.uni-jena.de/cat12/). We used the default settings detailed in the manual for CAT12 (http://dbm.neuro.uni-jena.de/cat12/CAT12-Manual.pdf) except for applying the affine regularization using the International Consortium for Brain Mapping template for East Asian brains. The individual T1-weighted images were finally segmented into gray matter, white matter and cerebrospinal fluid. The segmented gray matter images were smoothed with an FWHM of 8 mm. Regional gray matter volume of each region showing significant differences in FKBP5 diplotypes based neuroimaging analyses was extracted and then compared with a similar two-way ANCOVA. No significant diagnosis or FKBP5 main effect or their interaction effect on gray matter volume was found in these regions.
Discussion
In this study, we investigated the effects of PTSD diagnosis and stress-related gene FKBP5 on spontaneous brain activity in Han Chinese adults who had lost their only child. We found a main effect of PTSD diagnosis on low-frequency power (slow-4 and slow-5 bands) in the parietal lobules and inferior frontal gyrus, and a main effect of FKBP5 diplotype on low-frequency power (slow-4 band) in ACC and middle-frequency power (slow-3 band) in the motor/sensory areas irrespective of PTSD diagnosis. Also, there was a significant diagnosis × FKBP5 interaction effect on middle-frequency power (slow-3 band) in parietal lobule.
The prominent neuroanatomical theory of PTSD points to a decreased prefrontal inhibitory control over the amygdala49, supported by quantitative neuroimaging studies50,51 and meta-analytic reviews49,52. Other evidence also supports the cognitive-affective imbalance theory in the pathology of PTSD53,54, as evidenced by underactivated regions within the brain executive system and overactivated regions within the emotional processing system. The parietal lobule is a core component of the brain executive control network55,56. Decreased parietal lobule activity has been demonstrated in adults with a history of early life stress exposure57,58 and patients with PTSD54,59. Its activity has also been found to be negatively correlated with PTSD symptoms54. Both right60 and left61 inferior frontal gyri play an important role in attentional monitoring and inhibiting inappropriate responses. Reduced inferior frontal gyrus activity during a proactive inhibition task was reported in veterans with PTSD, relative to control veterans without PTSD62. Thus, in this study, the lower low-frequency power in parietal lobule and inferior frontal gyrus aligns with the findings in prior studies about PTSD. The marginally negative correlation between parietal lobule and CAPS in PTSD adults in this study also supports our findings.
Twin studies posit that PTSD is moderately heritable, with heritability estimates in the range of 28–46%63,64,65. Increased GR sensitivity is only found in trauma-exposed subjects with PTSD rather than those without PTSD6. As a critical modulator of GR sensitivity, FKBP5 is thought to be an interesting candidate gene for PTSD development66. Individuals with PTSD who carried two FKBP5 rs1360780 risk SNPs (TT) exhibited the lowest cingulum fractional anisotropy, compared to the same risk SNPs carriers without PTSD, and other genotypes both with and without PTSD14. The FKBP5 diplotypes with homozygote and heterozygote combinations of the H1 and H2 yin yang haplotypes provide more complete genetic information46,47. A recent study divided a group of healthy participants into high and low-risk genogroups according to the number of risk alleles of four FKBP5 SNPs, and found that the FKBP5 high-risk allele group demonstrated lower low-frequency power spectrum but greater high-frequency power spectrum in brain frontotemporoparietal network, compared to the low-risk allele group17. To the best of our knowledge, our current study is the first to further describe the effect of FKBP5 diplotypes on resting-state brain activity in trauma-exposed subjects both with and without PTSD diagnosis. In this study, at ACC, the zero-risk diplotype group (H1/H1) displayed higher low-frequency power spectrum than the middle (H1/H2)- and high (H2/H2)- risk diplotype groups; while at precentral/postcentral gyri and MCC, the H2/H2 diplotype group had higher middle spectral power than other two genogroups. Although using different study samples and different FKBP5 gene grouping methods, both our current results and findings from prior study17 suggest that the FKBP5 risk alleles may have different modulatory effects on brain activity within different frequency bands, where high-risk alleles may associate with lower low-frequency power but greater middle or high-frequency power in frontotemporoparietal regions. In the FKBP5 high-risk genogroup, the lower power spectrum in ACC may indicate impaired emotional processing, while the higher power spectrum in precentral/postcentral gyri may suggest hypervigilance/sensitization to threatening stimuli17. This interpretation is partially supported by a prior study on FKBP5 rs1360780 in subjects recruited from the general medical clinics, where risk allele (TC/TT) carriers showed attention bias toward the threat, compared to non-risk allele carriers (CC)67.
In the present study, another important finding is the interaction modulatory effect of FKBP5 × diagnosis in the parietal lobule, where in PTSD adults, the high-risk diplotype group was associated with lower spectral power than the other two genogroups, whereas, in non-PTSD adults the opposite occurred. One possible interpretation for the different modulatory effects of FKBP5 risk diplotype in PTSD and non-PTSD subjects is that some brain regions—such as the parietal lobule—may have a resilient or compensatory role in trauma-exposed subjects without PTSD. Further research is needed to gain a precise understanding of the role of FKBP5 risk diplotype on PTSD development. For example, more quantitative neuroimaging-genetics studies are required to provide reliable and repeatable intermediate phenotypes that would improve our understanding of GR sensitivity differences in subjects with PTSD and without PTSD.
This study has several limitations. First, our study only focused on the influence of losing an only child in China, and so we urge caution when applying these results to other traumatic experiences and other populations of different race or ethnicity. Second, due to the cross-sectional nature of the current study, it was unclear if the brain differences were present before the traumatic experience or if they occurred after the traumatic event. Third, an increasing number of studies have demonstrated that the FKBP5 risk SNPs often interact with individual childhood trauma, resulting in an increased risk of PTSD development and greater PTSD symptoms66,68. This evidence suggests the importance of gene × childhood environment interactions for PTSD. However, childhood trauma data was not originally collected in this study and thus needs to be taken into account in further studies. Finally, given the significant history of failure to replicate candidate-gene studies69,70, the findings here should be considered as preliminary results and need to be validated by studies involving large replication samples or using data from genome-wide association studies (GWAS) of PTSD.
Conclusion
In conclusion, our findings suggest that PTSD may impair lower low-frequency power in executive control network regions in Han Chinese adults who had lost their only child. The FKBP5 high-risk diplotype group displayed lower spectral power in ACC and greater power in the motor/sensory areas, suggesting a disturbance of emotional processing and hypervigilance/sensitization to threatening stimuli. The interaction effect of diagnosis × FKBP5 in parietal lobule may contribute to PTSD development in adults who experienced the loss of an only child.
References
Yehuda, R., Giller, E. L., Southwick, S. M., Lowy, M. T. & Mason, J. W. Hypothalamic-pituitary-adrenal dysfunction in posttraumatic stress disorder. Biol. Psychiatry 30, 1031–1048 (1991).
Yehuda, R. Post-traumatic stress disorder. N. Engl. J. Med. 346, 108–114 (2002).
Yehuda, R., Golier, J. A., Halligan, S. L., Meaney, M. & Bierer, L. M. The ACTH response to dexamethasone in PTSD. Am. J. Psychiatry 161, 1397–1403 (2004).
Yehuda, R., Boisoneau, D., Lowy, M. T. & Giller, E. L. Dose-response changes in plasma cortisol and lymphocyte glucocorticoid receptors following dexamethasone administration in combat veterans with and without posttraumatic stress disorder. Arch. Gen. Psychiatry 52, 583–593 (1995).
Kessler, R. C., Sonnega, A., Bromet, E., Hughes, M. & Nelson, C. B. Posttraumatic stress disorder in the National Comorbidity Survey. Arch. Gen. Psychiatry 52, 1048–1060 (1995).
Binder, E. B. et al. Association of FKBP5 polymorphisms and childhood abuse with risk of posttraumatic stress disorder symptoms in adults. JAMA 299, 1291–1305 (2008).
Binder, E. B. The role of FKBP5, a co-chaperone of the glucocorticoid receptor in the pathogenesis and therapy of affective and anxiety disorders. Psychoneuroendocrinology 34, S186–S195 (2009).
Pape, J. C. & Binder, E. B. The role of genetics and epigenetics in the pathogenesis of posttraumatic stress disorder. Psychiatr. Ann. 46, 510–518 (2016).
Wochnik, G. M. et al. FK506-binding proteins 51 and 52 differentially regulate dynein interaction and nuclear translocation of the glucocorticoid receptor in mammalian cells. J. Biol. Chem. 280, 4609–4616 (2005).
Hariri, A. R. & Weinberger, D. R. Imaging genomics. Br. Med. Bull. 65, 259–270 (2003).
Huang, H. et al. Imaging genomics. Pac. Symp. Biocomput. 23, 304–306 (2018).
Thompson, P. M., Martin, N. G. & Wright, M. J. Imaging genomics. Curr. Opin. Neurol. 23, 368–373. (2010).
Glahn, D. C., Paus, T. & Thompson, P. M. Imaging genomics: mapping the influence of genetics on brain structure and function. Hum. Brain Mapp. 28, 461–463 (2007).
Fani, N. et al. Structural and functional connectivity in posttraumatic stress disorder: associations with Fkbp5. Depress Anxiety 33, 300–307 (2016).
Roy, A., Gorodetsky, E., Yuan, Q., Goldman, D. & Enoch, M. A. Interaction of FKBP5, a stress-related gene, with childhood trauma increases the risk for attempting suicide. Neuropsychopharmacology 35, 1674–1683 (2010).
Bevilacqua, L. et al. Interaction between FKBP5 and childhood trauma and risk of aggressive behavior. Arch. Gen. Psychiatry 69, 62–70 (2012).
Bryant, R. A., Felmingham, K. L., Liddell, B., Das, P. & Malhi, G. S. Association of FKBP5 polymorphisms and resting-state activity in a frontotemporal-parietal network. Transl. Psychiatry 6, e925 (2016).
Morey, R. A. et al. Genome-wide association study of subcortical brain volume in PTSD cases and trauma-exposed controls. Transl. Psychiatry 7, 1265 (2017).
Malinen, S. et al. Aberrant temporal and spatial brain activity during rest in patients with chronic pain. Proc. Natl Acad. Sci. USA 107, 6493–6497 (2010).
Zang, Y. F. et al. Altered baseline brain activity in children with ADHD revealed by resting-state functional MRI. Brain Dev. 29, 83–91 (2007).
Han, Y. et al. Frequency-dependent changes in the amplitude of low-frequency fluctuations in amnestic mild cognitive impairment: a resting-state fMRI study. Neuroimage 55, 287–295 (2011).
Zuo, X. N. et al. The oscillating brain: complex and reliable. Neuroimage 49, 1432–1445 (2010).
Bai, F. et al. Association of a GSK-3β polymorphism with brain resting-state function in amnestic-type mild cognitive impairment. J. Alzheimers Dis. 32, 387–396 (2012).
Liu, X. et al. The association between TOMM40 gene polymorphism and spontaneous brain activity in amnestic mild cognitive impairment. J. Neurol. 261, 1499–1507 (2014).
Cui, L. et al. Relationship between the LHPP gene polymorphism and resting-state brain activity in major depressive disorder. Neural Plast. 2016, 9162590 (2016).
Hesketh, T., Lu, L. & Xing, Z. W. The effect of China’s one-child family policy after 25 years. N. Engl. J. Med. 353, 1171–1176 (2005).
Basten, S. & Jiang, Q. China’s family planning policies: recent reforms and future prospects. Stud. Fam. Plann. 45, 493–509 (2014).
First, M. B., Spitzer, R. L., Gibbon, M. & Williams, J. B. Structured Clinical Interview for DSM-IV-TR Axis I Disorders, Research Version, Patient Edition. (SCID-I/P). (Biometrics Research, New York State Psychiatric Institute, New York, 2002).
Hamilton, M. A rating scale for depression. J. Neurol. Neurosurg. Psychiatry 23, 56 (1960).
Hamilton, M. The assessment of anxiety states by rating. Br. J. Med. Psychol. 32, 50–55 (1959).
Folstein, M. F., Robins, L. N. & Helzer, J. E. The mini-mental state examination. Arch. Gen. Psychiatry 40, 812 (1983).
Cheng, Y. et al. Social support plays a role in depression in Parkinson’s disease: a cross-section study in a Chinese cohort. Parkinsonism Relat. Disord. 14, 43–45 (2008).
Jiang, X.-R., Du, J.-J. & Dong, R.-Y. Coping style, job burnout and mental health of university teachers of the millennial generation. J. Math. Sci. Technol. Educ. 13, 3379–3392 (2017).
Ashburner, J. A fast diffeomorphic image registration algorithm. Neuroimage 38, 95–113 (2007).
Friston, K. J., Williams, S., Howard, R., Frackowiak, R. S. & Turner, R. Movement‐related effects in fMRI time‐series. Magn. Reson. Med. 35, 346–355 (1996).
Yan, C. G. et al. A comprehensive assessment of regional variation in the impact of head micromovements on functional connectomics. Neuroimage 76, 183–201 (2013).
Yan, C. G. et al. Reduced default mode network functional connectivity in patients with recurrent major depressive disorder. Proc. Natl Acad. Sci. USA 116, 9078–9083 (2019).
Wang, L. et al. Frequency-dependent changes in amplitude of low-frequency oscillations in depression: a resting-state fMRI study. Neurosci. Lett. 614, 105–111 (2016).
Wang, Z. et al. Frequency-dependent amplitude alterations of resting-state spontaneous fluctuations in idiopathic generalized epilepsy. Epilepsy Res. 108, 853–860 (2014).
Liu, Y. et al. A rapid improved multiplex ligation detection reaction method for the identification of gene mutations in hereditary hearing loss. PLoS ONE 14, e0215212 (2019).
Altshuler, D., Donnelly, P. & Consortium, I. H. A haplotype map of the human genome. Nature 437, 1299 (2005).
Barrett, J. C., Fry, B., Maller, J. & Daly, M. J. Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics 21, 263–265 (2005).
Stephens, M., Smith, N. J. & Donnelly, P. A new statistical method for haplotype reconstruction from population data. Am. J. Hum. Genet. 68, 978–989 (2001).
Stephens, M. & Donnelly, P. A comparison of bayesian methods for haplotype reconstruction from population genotype data. Am. J. Hum. Genet. 73, 1162–1169 (2003).
Stephens, M. & Scheet, P. Accounting for decay of linkage disequilibrium in haplotype inference and missing-data imputation. Am. J. Hum. Genet. 76, 449–462 (2005).
Hoehe, M. R. et al. Multiple haplotype-resolved genomes reveal population patterns of gene and protein diplotypes. Nat. Commun. 5, 5569 (2014).
Lingjun, Z., Kesheng, W. & Xingguang, L. Use of diplotypes–matched haplotype pairs from homologous chromosomes-in gene-disease association studies. Shanghai Arch. Psychiatry 26, 165 (2014).
Worsley, K. J., Taylor, J. E., Tomaiuolo, F. & Lerch, J. Unified univariate and multivariate random field theory. Neuroimage 23, S189–S195 (2004).
Etkin, A. & Wager, T. D. Functional neuroimaging of anxiety: a meta-analysis of emotional processing in PTSD, social anxiety disorder, and specific phobia. Am. J. Psychiatry 164, 1476–1488 (2007).
Shin, L. M. et al. A functional magnetic resonance imaging study of amygdala and medial prefrontal cortex responses to overtly presented fearful faces in posttraumatic stress disorder. Arch. Gen. Psychiatry 62, 273–281 (2005).
Williams, L. M. et al. Trauma modulates amygdala and medial prefrontal responses to consciously attended fear. Neuroimage 29, 347–357 (2006).
Boccia, M. et al. Different neural modifications underpin PTSD after different traumatic events: an fMRI meta-analytic study. Brain Imaging Behav. 10, 226–237 (2016).
Hayes, J. P., VanElzakker, M. B. & Shin, L. M. Emotion and cognition interactions in PTSD: a review of neurocognitive and neuroimaging studies. Front. Integr. Neurosci. 6, 89 (2012).
Morey, R. A., Petty, C. M., Cooper, D. A., Labar, K. S. & McCarthy, G. Neural systems for executive and emotional processing are modulated by symptoms of posttraumatic stress disorder in Iraq War veterans. Psychiatry Res. 162, 59–72 (2008).
Seeley, W. W. et al. Dissociable intrinsic connectivity networks for salience processing and executive control. J. Neurosci. 27, 2349–2356 (2007).
Cabeza, R., Ciaramelli, E. & Moscovitch, M. Cognitive contributions of the ventral parietal cortex: an integrative theoretical account. Trends Cogn. Sci. 16, 338–352 (2012).
Philip, N. S. et al. Regional homogeneity and resting state functional connectivity: associations with exposure to early life stress. Psychiatry Res. 214, 247–253 (2013).
Pechtel, P. & Pizzagalli, D. A. Effects of early life stress on cognitive and affective function: an integrated review of human literature. Psychopharmacology 214, 55–70 (2011).
Aupperle, R. L. et al. Dorsolateral prefrontal cortex activation during emotional anticipation and neuropsychological performance in posttraumatic stress disorder. Arch. Gen. Psychiatry 69, 360–371 (2012).
Hampshire, A., Chamberlain, S. R., Monti, M. M., Duncan, J. & Owen, A. M. The role of the right inferior frontal gyrus: inhibition and attentional control. Neuroimage 50, 1313–1319 (2010).
Swick, D., Ashley, V. & Turken, U. Left inferior frontal gyrus is critical for response inhibition. BMC Neurosci. 9, 102 (2008).
van Rooij, S. J. et al. Impaired right inferior frontal gyrus response to contextual cues in male veterans with PTSD during response inhibition. J. Psychiatry Neurosci. 39, 330–338 (2014).
Stein, M. B., Jang, K. L., Taylor, S., Vernon, P. A. & Livesley, W. J. Genetic and environmental influences on trauma exposure and posttraumatic stress disorder symptoms: a twin study. Am. J. Psychiatry 159, 1675–1681 (2002).
Sartor, C. E. et al. Common heritable contributions to low-risk trauma, high-risk trauma, posttraumatic stress disorder, and major depression. Arch. Gen. Psychiatry 69, 293–299 (2012).
Afifi, T. O., Asmundson, G. J., Taylor, S. & Jang, K. L. The role of genes and environment on trauma exposure and posttraumatic stress disorder symptoms: a review of twin studies. Clin. Psychol. Rev. 30, 101–112 (2010).
Hawn, S. E. et al. GxE effects of FKBP5 and traumatic life events on PTSD: a meta-analysis. J. Affect. Disord. 243, 455–462 (2019).
Fani, N. et al. FKBP5 and attention bias for threat: associations with hippocampal function and shape. JAMA Psychiatry 70, 392–400 (2013).
Tamman, A. J. F. et al. Attachment style moderates effects of FKBP5 polymorphisms and childhood abuse on post-traumatic stress symptoms: Results from the National Health and Resilience in Veterans Study. World J. Biol. Psychiatry 20, 289–300 (2019).
Border, R. et al. No support for historical candidate gene or candidate gene-by-interaction hypotheses for major depression across multiple large samples. Am. J. Psychiatry 176, 376–387 (2019).
Logue, M. W. et al. The Psychiatric Genomics Consortium Posttraumatic Stress Disorder Workgroup: posttraumatic stress disorder enters the age of large-scale genomic collaboration. Neuropsychopharmacology 40, 2287 (2015).
Acknowledgements
This work was funded by the grants from the National Nature Science Foundation of China [Nos. 81671672, 81301209 to R.Q., and 81801678 to Y.L.]; the Jiangsu Provincial Medical Youth Talent [Nos. QNRC2016888 to R.Q., and QNRC2016207 to Y.L.]; the Chinese Key Grant [No. BWS11J063 to G.M.L.]; the Youth Nature Science Foundation of Jiangsu Province [No. BK20170223 to Y.L.]; and NIH [No. U54 EB020403 to P.M.T. and N.J.].
Author information
Authors and Affiliations
Corresponding authors
Ethics declarations
Conflict of interest
All authors declare no competing interests. P.M.T. and N.J. received a research grant from Biogen, Inc., (Boston, USA), for research unrelated to this manuscript.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Qi, R., Luo, Y., Zhang, L. et al. FKBP5 haplotypes and PTSD modulate the resting-state brain activity in Han Chinese adults who lost their only child. Transl Psychiatry 10, 91 (2020). https://doi.org/10.1038/s41398-020-0770-5
Received:
Revised:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41398-020-0770-5
This article is cited by
-
FKBP51 modulates hippocampal size and function in post-translational regulation of Parkin
Cellular and Molecular Life Sciences (2022)