Abstract
Losing one’s only child is a major traumatic life event that may lead to posttraumatic stress disorder (PTSD); however, not all parents who experience this trauma develop PTSD. Genetic variants are associated with the risk of developing PTSD. Catechol-O-methyltransferase (COMT) rs4680 and brain-derived neurotrophic factor (BDNF) rs6265 are two most well-described single-nucleotide polymorphisms that relate to stress response; however, the neural mechanism underlying their effects on adults who lost an only child remains poorly understood. Two hundred and ten Han Chinese adults who had lost their only child (55 with PTSD and 155 without PTSD) were included in this imaging genetics study. Participants were divided into subgroups according to their COMT rs4680 and BDNF rs6265 genotypes. Degree Centrality (DC)—a resting-state fMRI index reflecting the brain network communication—was compared with a three-way (PTSD diagnosis, COMT, and BDNF polymorphisms) analysis of covariance. Diagnosis state had a significant effect on DC in bilateral inferior parietal lobules and right middle frontal gyrus (MFG), where PTSD adults showed weaker DC. BDNF × diagnosis interaction effect was found in the right MFG and hippocampus, and these two regions were reversely modulated. Also, there was a significant COMT × BDNF interaction effect in left cuneus, middle temporal gyrus, right inferior occipital gyrus, and bilateral putamen, independent of PTSD diagnosis. These findings suggest that the modulatory effect of BDNF polymorphism on the MFG and hippocampus may contribute to PTSD development in bereaved adults. Interactions of COMT × BDNF polymorphisms modulate some cortices and basal ganglia, irrespective of PTSD development.
Similar content being viewed by others
Introduction
The “One-Child Policy”—which permits each family to have only one child—was implemented in mainland China for more than 30 years1,2. Although this policy succeeded in slowing the rapid population growth rate in China, its associated problems are also becoming apparent. For example, it has brought about a new phenomenon of “childless” older Chinese adults, these parents lost their only child at a time when they were unable to bear another. This causes them to experience long-term grief, loneliness, anxiety, depression, and posttraumatic stress disorder (PTSD)3,4,5. Although the One-Child Policy ended in 2015, the number of families that might lose their only child will possibly continue to increase for a considerable period of time6, which could lead to significant and widespread public health problems.
PTSD is characterized by dysregulated fear conditioning, extinction, and stress responses7,8,9. Functional neuroimaging studies have identified the dysfunction of several brain regions involved in the regulation of fear and stress responses in PTSD, which mainly includes the amygdala, hippocampus, and prefrontal cortex9,10,11. These findings strongly support the prefrontal-limbic imbalance theory in the pathology of PTSD8,11,12,13. It should be noted that not all adults who experience the trauma of losing their only child will go on to develop PTSD14. The underlying mechanisms of PTSD development remain inadequately investigated and poorly understood. Twin and heritability studies have consistently suggested that at least one-third of the variance in PTSD risk is determined by genetic components15,16,17. A variety of candidate genes in the dopaminergic system, hypothalamic-pituitary-adrenal axis, serotonergic system, and neuroinflammation have been conducted in PTSD studies18,19,20. Also, several unbiased genome-wide association studies (GWAS) based on the clinical diagnosis of PTSD have been performed21,22,23,24. These studies have helped us to understand the genetic basis of PTSD development and its genetic overlap with other mental disorders such as schizophrenia24. However, prior candidate-gene studies and even GWAS have largely been unsuccessful at determining an overlapped or replicable risk gene for PTSD19. The failure to replicate these results could possibly be due to the wide variations in the types of trauma, the severity of PTSD symptoms, and race/ethnicity populations involved in these genetic studies of PTSD25.
One way to alleviate the clinical heterogeneity of PTSD is to take advantage of intermediate phenotypes, especially the quantitative neuroimaging phenotypes26,27,28. A growing body of imaging genetics research has proven the association between neuroimaging measures and candidate genes in PTSD29,30,31. Catechol-O-methyltransferase (COMT) Val158Met (rs4680) and brain-derived neurotrophic factor (BDNF) Val66Met (rs6265) polymorphisms are two typical and well-described genetic variants that have been proven to be associated with stress response and resilience32,33,34. These two genetic variants have been extensively investigated in populations consisting of healthy subjects33,35 and those with anxiety disorders36,37, including PTSD29,31. Meanwhile, both COMT rs4680 and BDNF rs6265 regulate the brain dopamine system38—an important system in the pathogenesis of PTSD—by their important role in regulating fear memory emotion and behaviors. Dopaminergic dysfunction has also been reported to play a role in the arousal symptoms39 and fear conditioning to an aversive stimulus40, which are seen in PTSD patients. COMT catalyzes the degradation of catecholamines, particularly dopamine, and the Met158 and Val158 alleles of COMT rs4680 polymorphism have low and high enzyme activity, respectively41,42. BDNF is a key neurotrophin that regulates neuronal plasticity, neuronal survival, and neurogenesis, and the Met66 allele of BDNF rs6265 polymorphism is associated with lower BDNF release compared with the Val66 allele, resulting in reduced activity-dependent dopamine release43. Both COMT rs4680 and BDNF rs6265 polymorphisms are suggested to moderate several neuroimaging structural phenotypes observed in PTSD29,31,44. In addition, recently, the interaction effect between COMT rs4680 and BDNF rs6265 on brain function has been reported in healthy subjects45,46. However, to the best of our knowledge, no study to date has investigated the effects of COMT rs4680, BDNF rs6265, or COMT by BDNF interaction on brain function in PTSD patients. So far, little is known about the effects of COMT rs4680 and BDNF rs6265 genotypes on brain function in Chinese adults who lost their only child.
In this study, we aimed to examine the main effects of COMT rs4680, BDNF rs6265, and their interaction effect on brain function as seen using resting-state functional magnetic resonance imaging (fMRI) in Chinese adults who had lost their only child. Degree Centrality (DC)47,48—an important fMRI-derived metric that measures the strength of intrinsic connectivity between a given voxel and others in the brain—was used in this study for three reasons. First, DC is a physiologically meaningful and powerful biomarker in exploring the communication potential/strength of a certain brain region47,48,49,50. Second, the DC algorithm has been successfully applied in studying many mental disorders, such as obsessive-compulsive disorder51,52, bipolar disorder53,54, and one published fMRI study involving Chinese adults who lost their only child (22 trauma-exposed adults without PTSD diagnosis were recruited in that study)55. Third, recent studies using DC have demonstrated modulatory effects of COMT56,57, BDNF45, and even COMT by BDNF interaction45 on the connectivity strength of the brain’s functional hubs. In this preliminary study, we hypothesized that COMT rs4680 and BDNF rs6265 would exhibit different modulatory effects on the DC in prefrontal-limbic regions in Han Chinese adults with loss of an only child and diagnosed with PTSD, when compared with those without PTSD diagnosis.
Subjects and methods
Subjects
We performed a PTSD survey in Han Chinese adults who had lost their only child from Jiangsu Province, China, between September 2016 and March 2017. This study was approved by the Medical Research Ethics Committee of Jiangsu University and all participants provided written informed consents. All 237 Han adults (ages 40–67 years, mean age of 58.54 years) who had lost their only child—without other major traumatic exposures detected on the clinician-administered PTSD scale (CAPS) life events checklist—were successfully interviewed and screened by the CAPS. The Chinese version of the structured clinical interview for DSM-IV58 (revised by Professor Lipeng Fei from Beijing Hui Long Guan Hospital) was used to screen all these participants. After this procedure, 57 trauma-exposed adults were diagnosed with PTSD (19 out of these 57 PTSD adults had comorbid major depressive disorder (MDD), 3 had comorbid generalized anxiety disorder (GAD), and one had comorbidities of both MDD and GAD). One hundred and seventy trauma-exposed adults did not meet any diagnostic criteria of mental illness or substance-use disorders. Ten trauma-exposed adults were diagnosed with other psychiatric disorders (five diagnosed with MDD, four with GAD, and one with mixed generalized anxiety and MDD) and they were excluded in this study.
The exclusion criteria for the subsequent fMRI study were as follows: any current or history of brain injury or other major medical or neurological conditions (five trauma-exposed adults without PTSD were ruled out: four had cerebral infarction or ischemia and one had a history of MDD and corresponding antidepressant drug therapy), any MRI contraindication (none), and left-handedness (none).
Measures
Each participant was assessed using a set of neuropsychological tests, which included the Hamilton Depression (HAMD)59 and Hamilton Anxiety (HAMA)60 rating scales, the Mini-Mental State Examination61, Chinese Social Support Rating Scale (SSRS)62, and individual Simple Coping Style Questionnaire (SCSQ)63. The SSRS contains 3 subscales of social support: subjective support (which reflects the perceived interpersonal network that an individual can count on, 4 items with scores ranging from 8 to 32); objective support (actual support an individual received, 3 items with scores ranging from 1 to 22); and the utility of support (the pattern of behavior that an individual exhibits when seeking social support, 3 items with scores ranging from 3 to 12). Higher scores for the SSRS indicate stronger social support and the total support score (ranging from 12 to 66) is the sum of all the sub-items. The SCSQ contains assessments of both active and negative coping, containing 12 and 8 items, respectively. The scale of each SCSQ item uses 4-level Likert score standards, in which “3” stands for regular use, whereas “0” stands for no use. Next, the scores for active and negative coping are calculated independently, and a higher score indicates the inclination to adopt the corresponding coping style, whereas the coping tendency score is defined as the active coping score minus the negative coping score.
MRI data acquisition
MRI data were acquired using a 3 Tesla MR scanner (Achieva 3.0 TTX; Philips, Amsterdam, The Netherlands). Foam pads were applied to minimize head motion. The participants were instructed to stay still during scanning and keep their eyes closed but not fall asleep. First, high-resolution T1-weighted structural images were acquired with the three-dimensional turbo fast echo sequence (repetition time/echo time [(TR/TE) = 9.7 ms/4.6 ms, field of view (FOV) = 256 × 256 mm2, flip angle = 9°, matrix size = 256 × 256, 160 sagittal slices with thickness of 1 mm). Second, rs-fMRI data were acquired with a single-shot, gradient-recalled echo-planar imaging sequence (TR/TE = 2000 ms/30 ms, FOV = 192 × 192 mm2, flip angle = 90°, matrix = 64 × 64, voxel size = 3 × 3 × 4 mm3). For each subject, 230 brain volumes with 35 axial slices of functional data were collected in 460 s.
Data preprocessing
MRI data were preprocessed using Data Processing Assistant for Resting-State fMRI (DPARSF; http://rfmri.org/DPARSF)64, which is essentially based on SPM12 software (http://www.fil.ion.ucl.ac.uk/spm). The first ten volumes of each subject were excluded for steady-state longitudinal magnetization; the remaining 220 volumes were corrected for temporal differences and head motion. Each subject’s T1-weighted image was co-registered to the functional images and then segmented into gray matter, white matter, and cerebrospinal fluid, and then transformed into the standard Montreal Neurological Institute (MNI) space using the Diffeomorphic Anatomical Registration Through Exponentiated Lie algebra. The fMRI data were then transformed into the MNI stereotaxic space of 3 × 3 × 3 mm3, using the parameters of the T1 image normalization and smoothed with an 8 mm full width at half maximum (FWHM) isotropic Gaussian kernel.
Quality control and nuisance regression
To address concerns about head motion, according to recent recommendations in minimizing head motion confounds65,66, we used the Friston 24-parameter model67 to regress out head motion effects. Head translations, rotations, and the framewise displacement (using the Jenkinson formula) from each participant were calculated. We excluded seven participants (two PTSD adults and five trauma-exposed controls) for head translations >1.5 mm or rotations >1.5°, and two participants (trauma-exposed controls) for mean framewise displacement >2.5 SD. We also included the mean framewise displacement as a nuisance covariate in the following imaging genetic analyses. Other sources of spurious variance (mean signals from cerebrospinal fluid and white matter) were also regressed out. The imaging data were also temporally filtered (bandpass: 0.01–0.1 Hz).
After the quality control, 9 subjects were excluded, and 55 PTSD adults and 158 trauma-exposed controls remained.
Degree centrality analysis
To measure the DC of each voxel throughout the brain, voxel-wise whole-brain correlation analysis with a threshold of r = 0.25 was performed51,54. We also analyzed the DC with other r thresholds (0.2, 0.3, and 0.4, respectively) (please refer to “Validation Analysis” in the Results section for details). The DC of a voxel was calculated as the sum of connectivity between a given voxel and other voxels, which exceeded the r threshold47,48. For standardization purposes, the DC of each voxel was divided by the global mean DC value.
DNA genotyping
Of all the subjects enrolled in this study, three trauma-exposed adults without PTSD refused blood collection. DNA for other participants was obtained from peripheral blood samples. Single-nucleotide polymorphism (SNP) genotyping of COMT rs4680 and BDNF rs6265 (detailed primers were listed in Supplementary Table S1) were conducted using the Improved Multiple Ligase Detection Reaction technique developed by Genesky Biotechnologies, Inc. (Shanghai, China)68. About 5% of the samples were randomly selected for confirmation and the results were 100% concordant.
Statistical analysis
Hardy–Weinberg equilibrium test of COMT rs4680 and BDNF rs6265 was performed using R version 3.5.3 (https://www.r-project.org). SPSS version 25 (IBM, Corp, Armonk, New York, USA) was used to analyze the clinical and psychological data. A three-way (diagnosis of PTSD, COMT rs4680, and BDNF rs6265 genotypes) analysis of variance was used to evaluate the main effects of PTSD diagnosis, COMT rs4680, BDNF rs6265 polymorphisms, and their interaction effects for all the demographic and psychological data. A voxel-wise three-way (diagnosis, COMT, and BDNF) analysis of covariance (ANCOVA) was performed with SPM12 to assess the individual effect of PTSD diagnosis, COMT rs4680, BDNF rs6265, and their interaction effects on brain DC map, adjusting for the effects of age, sex, educational level, duration since child-loss trauma, and head motion of framewise displacement. Significant clusters were determined by the Alphasim program, which was implemented in REST1.8 (http://www.restfmri.net/forum/index.php). The threshold was set at combined P < 0.005 for each voxel and a cluster size larger than 48 voxels (parameters: FWHM = 8 mm; the number of Monte Carlo simulations = 1000; within the gray matter mask consisting of 67,541 voxels), which corresponds to a corrected P < 0.05 (https://afni.nimh.nih.gov/pub/dist/doc/manual/AlphaSim.pdf).
A partial correlation analysis was performed to examine the relationship between DC values of those regions showing significant effects from above ANCOVA and CAPS, SSRS, SCSQ, HAMA, and HAMD, adjusting for age, sex, educational level, duration since child-loss trauma, and head motion. Correlation results were corrected for multiple comparisons using the Bonferroni correction for the number of regions where altered DC was detected from the ANCOVA (cutoff P-values of 0.05/10 = 0.005 in this study, corresponding to all ten regions showing significance).
Results
Clinical and psychological data
The flowchart of the study population is shown in Supplementary Fig. S1. All 55 PTSD adults and 155 trauma-exposed controls without PTSD were included in the final neuroimaging genetic analyses (Table 1). The distributions of both COMT rs4680 (PTSD group: 31 Val/Val, 22 Met/Val, and 2 Met/Met; trauma-exposed control group: 88 Val/Val, 56 Met/Val, and 11 Met/Met) and BDNF rs6265 (PTSD group: 17 Val/Val, 27 Met/Val, and 11 Met/Met; control group: 38 Val/Val, 83 Met/Val, and 34 Met/Met) were in Hardy–Weinberg equilibrium (all P > 0.05). For each polymorphism, participants were divided into two subgroups (Val-homozygous and Met-allele carriers) based on their genotypes, following methods in previous studies46,69 and then grouped based on both COMT rs4680 and BDNF rs6265 polymorphisms (as detailed in Supplementary Table S2). Among the PTSD samples, there were 8 COMT Met-BDNF Val/Val, 16 COMT Met-BDNF Met, 9 COMT Val/Val-BDNF Val/Val, and 22 COMT Val/Val-BDNF Met; within trauma-exposed controls, there were 16 COMT Met-BDNF Val/Val, 51 COMT Met-BDNF Met, 22 COMT Val/Val-BDNF Val/Val, and 66 COMT Val/Val-BDNF Met.
There were no significant differences between PTSD adults and trauma-exposed controls in age, educational level, SSRS, or SCSQ (P > 0.05), but PTSD adults had higher CAPS, HAMA, and HAMD scores, a lower male-to-female ratio, and shorter duration since losing the child (Table 1). There was a significant COMT × BDNF interaction effect on scores of active coping style (F = 4.43, P = 0.04) in all participants, irrespective of PTSD diagnosis, which was driven by the lower score in the COMT Met-BDNF Val/Val subgroup than the COMT Met-BDNF Met subgroup (P = 0.05), whereas no significant difference was found between COMT Val/Val-BDNF Val/Val and COMT Val/Val-BDNF Met subgroups (P = 0.20). There were no significant COMT or BDNF main effects, or two-way or three-way interactions for all other clinical and psychological data.
The main effect of PTSD diagnosis on DC
Significant PTSD diagnosis effects on DC were found in bilateral inferior parietal lobules (IPL) and right middle frontal gyrus (MFG). Post-hoc analysis showed that PTSD adults had weaker DC in these three regions, relative to trauma-exposed adults without PTSD (Fig. 1 and Supplementary Table S3).
The individual effect of COMT and BDNF genotypes on DC
No significant main effect of COMT or BDNF genotypes was found in this study. However, there was a significant BDNF genotype × diagnosis interaction effect in the right MFG (Fig. 2 and Supplementary Table S3), which accounted for 9% of the interindividual variance in DC in right MFG (partial η2 = 9%, estimated using SPSS 25). At this locus, in PTSD patients, the BDNF Val/Val genotype was associated with weaker DC compared with its Met carrier counterpart (P < 0.001), whereas in trauma-exposed controls, there was no significant difference between these two genogroups (P = 0.17).
In this study, the right hippocampus—a core region affected by the pathology of PTSD—only passed the P < 0.005 threshold but not the cluster size threshold (Supplementary Table S3); we then used a lenient statistical threshold to investigate its potential change56. A 20 mm-radius box centered at the peak location of the hippocampus was placed based on the whole brain results of BDNF × diagnosis interaction effect (P < 0.005), then a family-wise error small-volume correction (SVC), which implemented in SPM12, was conducted to correct for the multiple comparisons within the box56. After this procedure, a marginal significance of the right hippocampus was found (P = 0.058, SVC corrected). The BDNF genotype × diagnosis interaction accounted for 9% variance in DC in this region (η2 = 9%). Post-hoc analysis showed that at the hippocampus, in PTSD adults, the BDNF Val/Val genotype was associated with higher DC than its Met carrier counterpart (P < 0.001), whereas in trauma-exposed controls the opposite occurred (P = 0.046) (Fig. 2).
Interaction of COMT and BDNF genotypes on DC
There were significant interactions between the effects of the COMT and BDNF genotypes in the regions of left cuneus (η2 = 10%), left middle temporal gyrus (η2 = 11%), right inferior occipital gyrus (η2 = 10%), and bilateral putamen (right: η2 = 13%; left: η2 = 8%) in all participants, irrespective of PTSD diagnosis (Fig. 3 and Supplementary Table S3). Interestingly, after conducting the Curve Estimation procedure with SPSS 25, the quadratic regression was significant in all these regions by sorting according to COMT and BDNF genotypes (P < 0.05). The distribution of the DC was a U-shaped curve in the cortical regions, whereas for the putamen it showed an inverted U-shaped curve (Fig. 3), according to the presumed dopamine signaling from high to low (COMT Met-BDNF Val/Val genogroup > COMT Met-BDNF Met and COMT Val/Val-BDNF Val/Val > COMT Val/Val-BDNF Met)45. Specifically, in the regions of the left cuneus, left middle temporal gyrus, and right inferior occipital gyrus, the DC was higher in the COMT Met-BDNF Val/Val subgroup than in the COMT Met-BDNF Met subgroup, but weaker in COMT Val/Val-BDNF Val/Val than in COMT Val/Val-BDNF Met (except for the region of the left cuneus, where no significant difference was found between COMT Val/Val-BDNF Val/Val and COMT Val/Val-BDNF Met subgroups, with P = 0.33). In contrast, for the region of bilateral putamen, the DC was weaker in COMT Met-BDNF Val/Val than in COMT Met-BDNF Met subgroups, but higher in COMT Val/Val-BDNF Val/Val than in COMT Val/Val-BDNF Met subgroups (in the left putamen, there was no significant difference but a slight trend of COMT Met-BDNF Val/Val < COMT Met-BDNF Met, with P = 0.08). No significant effect of diagnosis × COMT × BDNF interaction was found in this study.
Correlation analysis
For regions showing significant effect of diagnosis × BDNF interaction, a negative partial correlation was found between the right hippocampus and right MFG DC values, only in the PTSD group (r = −0.34; P = 0.01; Fig. 4). However, after multiple testing correction, it was no longer statistically significant. No significant correlation was found between regions showing significant PTSD main effect, COMT × BDNF interaction, and clinical or psychological indices.
Validation analysis
We assessed the validity of our findings through several strategies: first, we used different correlation coefficient thresholds for the DC measurement, with r-values = 0.2, 0.3, and 0.4, respectively. Similar results were found, except for the loss of significance of COMT × BDNF interaction effect in the region of inferior occipital gyrus at r = 0.3 and 0.4, after the same multiple comparisons correction. Second, as PTSD adults showed higher HAMA and HAMD scores than those without PTSD, we also performed the statistical analysis by additionally including anxiety and depression (along with gender and age, education, duration, and head motion) as covariates to evaluate the effect of anxiety and depression on the results. We found that some regions’ effects became insignificant after the multiple comparisons (Supplementary Table S3). This finding may suggest that part of the changes could be accounted for by a higher level of anxiety and depression, which was in line with the frequently observed accompanying depression and anxiety symptoms in PTSD subjects in actual clinical practice70,71. Third, to assess the possible effect of volume differences in detected regions, we extracted regional gray matter volume of each region from our findings and compared each of them with a similar three-way ANCOVA, and found no significant effects on gray matter volume of these regions (a detailed description was provided in the Supplementary Note 1). Fourth, for head motion control, despite concerns that scrubbing could compromise graph construction for DC analysis by increasing the likelihood of extreme correlation values, we conducted the scrubbing (removing time points with framewise displacement >0.2 mm)—along with the same strategy at the analyses to verify the results. All results remained the same after using this aggressive head motion control strategy. Furthermore, to ensure that our results are specific to the effects of BDNF and COMT rather than other dopamine-related genes, we further included five other dopamine-related SNPs—two DRD2 SNPs (rs2075652 and rs2134655)72 and three DRD3 SNPs (rs4646996, rs7131056, and rs9868039)73—in a multiple liner regression model to investigate their effects on the brain regions from our current results (SNP main effect and SNP × diagnosis interaction effect). After adopting the Bonferroni correction threshold (P = 0.005), we did not find any statistically significant effects of these five dopamine SNPs (a detailed description was provided in the Supplementary Note 2).
Discussion
In this study, we investigated the effects of PTSD diagnosis, COMT rs4680, and BDNF rs6265 polymorphisms on brain functional connectivity strength in Han Chinese adults who had lost their only child. We found the main effect of PTSD diagnosis on DC in the regions of IPL and MFG, and the interaction effect of BDNF × diagnosis in MFG and hippocampus. In addition, there was significant COMT × BDNF interaction effect in regions of the left cuneus, left middle temporal gyrus, right inferior occipital gyrus, and bilateral putamen in all participants, irrespective of PTSD diagnosis.
Converging evidence has suggested prefrontal-limbic imbalance along with a hypoactive (prefrontal)–hyperactive (limbic) gradient in PTSD12. For example, an increasing number of neuroimaging studies74,75 and meta-analytic reviews12,76 consistently demonstrated decreased activity in prefrontal regions and increased activity in limbic regions in PTSD patients. In supplement to prefrontal-limbic imbalance, the cognitive-affective imbalance theory8,77 was also propounded, where brain executive system (under-activated) and emotional processing system (over-activated) are differentially affected by PTSD. The IPL and dorsolateral prefrontal cortex (DLPFC, lies within the MFG) are two core parts of the brain executive-control network78 and both of them have been demonstrated to be susceptible to the effects of stress, partially due to the fact that these two regions developmentally mature later than other regions79,80. Decreased IPL or DLPFC activity was often reported in adults with a history of early life stress exposure79,81 and patients with PTSD77,82. Besides, decreased connectivity strength in IPL and DLPFC were detected by Liu et al.55 in a small group of Chinese adults who lost their only child (compared with healthy controls without trauma exposure). Thus, the weaker DC in IPL and MFG observed in the current study aligns with the findings of prior studies and provides further insight into understanding the changes of their communication potential/strength when PTSD occurs in bereaved adults.
Besides, we found an interaction effect of BDNF × diagnosis in the present study. As for MFG, the BDNF Val/Val genotype was associated with weaker DC than its Met carrier counterpart in PTSD adults, but for trauma-exposed controls, there was an opposing trend. A nonlinear, inverted U-shaped relationship between dopamine levels with spatial working memory and related cortical (prefrontal) function has been robustly demonstrated by a series of studies, both in non-human primates83 and humans84. This inverted U-shaped curve predicts that too much or too little dopamine can impair cortical function, and this has been applied to investigate the disruption of brain function in several mental disorders85,86,87. For example, Prata et al.85 reported that the cortical efficiency in patients with schizophrenia may shift to the left limb of the inverted U-shaped curve. In Parkinson’s disease, the attention function and prefrontal function tended to shift to the right limb of the curve (based on putative prefrontal DA levels)87. In healthy controls, the BDNF Met polymorphism was thought to be associated with higher cognitive deficit risk and lower frontal activity, especially within the context of working memory88. Dopamine plays a crucial role in neurological processes such as reward, motivation, and stress89. Increasing evidence has demonstrated dopaminergic hyperactivity as an underlying mechanism in PTSD-related pathophysiology90. Single-photon emission computed tomography study also found increased dopamine transporter density in bilateral striatum in PTSD patients91, which was interpreted as a contributing factor to the perpetuation and potentiation of exaggerated fear responses to a particular event associated with the traumatic experience. In the current study, as the increased putative dopamine signaling in BDNF Val/Val genotype was associated with decreased MFG DC in the PTSD group, we speculate that the inverted U-shaped relationship between dopamine signaling and frontal function may shift to the right limb in these bereaved PTSD adults (Supplementary Fig. S2). Although a few studies reported that the BDNF Met allele may be related to increased PTSD susceptibility92, other studies93 and recent reviews94,95 did not support this relationship. A precise understanding of the role of BDNF rs6265 polymorphism on PTSD development is required in further studies.
Another interesting finding in this study was that the interaction effect of BDNF × diagnosis in the hippocampus was almost opposite to the MFG, where the BDNF Val/Val genotype was associated with higher hippocampal DC compared with its Met carrier counterpart in PTSD adults but reversed in trauma-exposed controls. Besides, the DC values of MFG and hippocampus were anticorrelated in the PTSD group. As the hippocampus here was just marginally significant after multiple testing corrections, we explain the underlying biological mechanisms with caution. However, several explanations could be proposed based on current knowledge. The hippocampus is essential for memory functions, especially to memorize facts and events, and memory consolidation96. A series of studies have demonstrated greater hippocampal engagement in response to negative stimuli and trauma-specific cues, but less engagement during exposure to positive images97. Prefrontal-limbic imbalance in PTSD is thought to result in the failure to inhibit the negative memory8,12. Here we speculate a potential U-shaped relationship between the presumed dopamine levels and hippocampus function (Supplementary Fig. S2), in which both lower dopamine in trauma-exposed controls and higher dopamine in PTSD adults would be associated with lower hippocampal efficiency. This is shown as increased hippocampal DC, which may be related to overgeneralized relevance to negative stimuli and trauma-specific cues (may indicate inefficiency/disturbance of hippocampal connectivity)12. This speculation was supported by a study in a large number of healthy subjects, in which the BDNF Val/Val genotype subgroup had higher hippocampal activity than the Met carrier genotype subgroup during a working memory task, but accomplished with poorer episodic memory43. However, this speculation needs to be validated with a larger sample and more studies in the future.
In the present study, COMT × BDNF interaction effects were found in several temporal and occipital regions and putamen, irrespective of PTSD diagnosis. Besides, the distribution of the DC is more likely to be a U-shaped curve in these cortical regions and an inverted U-shaped curve in the putamen, according to the presumed dopamine signaling from high to low. The U-shaped and inverted U-shaped curves modulated by dopamine signaling have recently been reported in several studies45,98,99. For example, Zhao et al.98 reported network-dependent modulations of dopamine signaling on brain function, in which the intra-network functional connectivities of the brain “processing system” and “control system” were differently modulated by the COMT × DRD2 interaction effects. Based on current substantiated knowledge that the prefrontal and striatal dopamine systems in the brain have an antagonistic relationship, particularly in pathological circumstances100,101, we propose that there might be an anti-relationship of brain dopamine modulation in cortical regions and basal ganglia regions. A future study that includes both trauma-exposed subjects and non-traumatized healthy controls would be suitable to validate whether this phenomenon is unique in trauma-exposed subjects.
This study has several limitations. First, our study only focused on the influence of losing an only child, a unique phenomenon accompanying One-Child Policy in China, and so we urge caution when generalizing these resulting interpretations to other traumatic events and other race/ethnicity populations. Second, the cross-sectional nature of our measurement did not allow us to ascertain whether the abnormalities in PTSD and the effect of gene modulation were present before the traumatic experience, or if they occurred after the traumatic event. Third, we note that several other genes—such as dopamine transporter and receptor genes—involved in dopamine functioning102 should be addressed in future studies to clarify the modulatory effects of dopamine pathway genes on brain function103. Finally, given the history of inconsistent replication of candidate-gene studies104,105, the findings in this study should be considered preliminary and need to be validated by studies using large replication samples or using data from GWAS of PTSD19,106.
In conclusion, our findings provide evidence to suggest that the modulatory effect of BDNF polymorphism on the MFG and hippocampus may contribute to PTSD development in bereaved adults. Interactions of COMT × BDNF polymorphisms modulate some cortices and basal ganglia, irrespective of PTSD development. In future, tests of interventions and treatment effects in patients with PTSD may benefit from stratifying patients into groups by genotype, or at least from modeling genotype effects on brain metrics that may differ by genotype. Eventually, there may be sufficient data to test for treatment × genotype interaction effects on brain phenotypes relevant to PTSD, a key goal of precision medicine.
References
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. Plan. 45, 493–509 (2014).
Zheng, Y., Lawson, T. R. & Anderson Head, B. “Our only child has died”—a study of bereaved older Chinese parents. OMEGA J. Death Dying 74, 410–425 (2017).
Wei, Y., Jiang, Q. & Gietel-Basten, S. The well-being of bereaved parents in an only-child society. Death Stud. 40, 22–31 (2016).
Song, Y. Losing an only child: the one-child policy and elderly care in China. Reprod. Health Matters 22, 113–124 (2014).
Zhang, Y. & Jia, X. A qualitative study on the grief of people who lose their only child: from the perspective of familism culture. Front. Psychol. 9, 869 (2018).
Norrholm, S. D. & Jovanovic, T. Fear processing, psychophysiology, and PTSD. Harv. Rev. Psychiatry 26, 129–141 (2018).
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).
Shin, L. M. & Liberzon, I. The neurocircuitry of fear, stress, and anxiety disorders. Neuropsychopharmacology 35, 169–191 (2010).
Bremner, J. D. Neuroimaging in posttraumatic stress disorder and other stress-related disorders. Neuroimaging Clin. N. Am. 17, 523–538, ix (2007).
Hughes, K. C. & Shin, L. M. Functional neuroimaging studies of post-traumatic stress disorder. Expert Rev. Neurother. 11, 275–285 (2011).
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).
Pitman, R. K. et al. Biological studies of post-traumatic stress disorder. Nat. Rev. Neurosci. 13, 769–787 (2012).
Wang, Q., Xu, W., Ren, L., Wang, W. & Wang, Y. The relationship between hope and post-traumatic stress disorder in Chinese shidu parents: the mediating role of perceived stress. J. Affect. Disord. 251, 23–30 (2019).
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 genetic and environmental contributions to post-traumatic stress disorder and alcohol dependence in young women. Psychol. Med. 41, 1497–1505 (2011).
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).
Sharma, S. & Ressler, K. J. Genomic updates in understanding PTSD. Prog. Neuropsychopharmacol. Biol. Psychiatry 90, 197–203 (2019).
Nievergelt, C. M. et al. Genomic approaches to posttraumatic stress disorder: the Psychiatric Genomic Consortium Initiative. Biol. Psychiatry 83, 831–839 (2018).
Zhang, K. et al. An overview of posttraumatic stress disorder genetic studies by analyzing and integrating genetic data into genetic database PTSDgene. Neurosci. Biobehav. Rev. 83, 647–656 (2017).
Logue, M. W. et al. A genome-wide association study of post-traumatic stress disorder identifies the retinoid-related orphan receptor alpha (RORA) gene as a significant risk locus. Mol. Psychiatry 18, 937–942 (2013).
Xie, P. et al. Genome-wide association study identifies new susceptibility loci for posttraumatic stress disorder. Biol. Psychiatry 74, 656–663 (2013).
van der Merwe, C. et al. Concordance of genetic variation that increases risk for anxiety disorders and posttraumatic stress disorders and that influences their underlying neurocircuitry. J. Affect. Disord. 245, 885–896 (2019).
Duncan, L. E. et al. Largest GWAS of PTSD (N = 20 070) yields genetic overlap with schizophrenia and sex differences in heritability. Mol. Psychiatry 23, 666 (2018).
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).
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).
Hariri, A. R. & Weinberger, D. R. Imaging genomics. Br. Med. Bull. 65, 259–270 (2003).
Thompson, P. M., Martin, N. G. & Wright, M. J. Imaging genomics. Curr. Opin. Neurol. 23, 368–373 (2010).
Jin, M. J., Jeon, H., Hyun, M. H. & Lee, S. H. Influence of childhood trauma and brain-derived neurotrophic factor Val66Met polymorphism on posttraumatic stress symptoms and cortical thickness. Sci. Rep. 9, 6028 (2019).
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).
Hayes, J. P. et al. COMT Val158Met polymorphism moderates the association between PTSD symptom severity and hippocampal volume. J. Psychiatry Neurosci. 42, 95–102 (2017).
Feder, A., Nestler, E. J. & Charney, D. S. Psychobiology and molecular genetics of resilience. Nat. Rev. Neurosci. 10, 446–457 (2009).
Kang, J. I., Kim, S. J., Song, Y. Y., Namkoong, K. & An, S. K. Genetic influence of COMT and BDNF gene polymorphisms on resilience in healthy college students. Neuropsychobiology 68, 174–180 (2013).
Witte, A. V. et al. Interaction of BDNF and COMT polymorphisms on paired-associative stimulation-induced cortical plasticity. J. Neurosci. 32, 4553–4561 (2012).
Jung, Y. H. et al. Influence of brain-derived neurotrophic factor and catechol O-methyl transferase polymorphisms on effects of meditation on plasma catecholamines and stress. Stress 15, 97–104 (2012).
Hovatta, I. & Barlow, C. Molecular genetics of anxiety in mice and men. Ann. Med. 40, 92–109 (2008).
Smoller, J. W., Gardner-Schuster, E. & Misiaszek, M. Genetics of anxiety: would the genome recognize the DSM? Depression Anxiety 25, 368–377 (2008).
Savitz, J., Solms, M. & Ramesar, R. The molecular genetics of cognition: dopamine, COMT and BDNF. Genes Brain Behav. 5, 311–328 (2006).
Charney, D. S., Deutch, A. Y., Krystal, J. H., Southwick, S. M. & Davis, M. Psychobiologic mechanisms of posttraumatic stress disorder. Arch. Gen. Psychiatry 50, 294–305 (1993).
Pezze, M. A. & Feldon, J. Mesolimbic dopaminergic pathways in fear conditioning. Prog. Neurobiol. 74, 301–320 (2004).
Lotta, T. et al. Kinetics of human soluble and membrane-bound catechol O-methyltransferase: a revised mechanism and description of the thermolabile variant of the enzyme. Biochemistry 34, 4202–4210 (1995).
Tunbridge, E. M., Harrison, P. J. & Weinberger, D. R. Catechol-o-methyltransferase, cognition, and psychosis: Val158Met and beyond. Biol. Psychiatry 60, 141–151 (2006).
Egan, M. F. et al. The BDNF val66met polymorphism affects activity-dependent secretion of BDNF and human memory and hippocampal function. Cell 112, 257–269 (2003).
Schulz-Heik, R. J. et al. Catechol-O-methyltransferase Val158Met polymorphism moderates anterior cingulate volume in posttraumatic stress disorder. Biol. Psychiatry 70, 1091–1096 (2011).
Li, W., Liu, B., Xu, J., Jiang, T. & Yu, C. Interaction of COMT rs4680 and BDNF rs6265 polymorphisms on functional connectivity density of the left frontal eye field in healthy young adults. Hum. Brain Mapp. 37, 2468–2478 (2016).
Chen, W. et al. Interaction effects of BDNF and COMT genes on resting-state brain activity and working memory. Front. Hum. Neurosci. 10, 540 (2016).
Tomasi, D. & Volkow, N. D. Functional connectivity density mapping. Proc. Natl Acad. Sci. USA 107, 9885–9890 (2010).
Zuo, X. N. et al. Network centrality in the human functional connectome. Cereb. Cortex 22, 1862–1875 (2012).
Tomasi, D., Shokri-Kojori, E. & Volkow, N. D. High-resolution functional connectivity density: hub locations, sensitivity, specificity, reproducibility, and reliability. Cereb. Cortex 26, 3249–3259 (2016).
Tomasi, D., Wang, R., Wang, G.-J. & Volkow, N. D. Functional connectivity and brain activation: a synergistic approach. Cereb. Cortex 24, 2619–2629 (2013).
Beucke, J. C. et al. Abnormally high degree connectivity of the orbitofrontal cortex in obsessive-compulsive disorder. JAMA Psychiatry 70, 619–629 (2013).
Li, P. et al. The effects of cognitive behavioral therapy on resting‐state functional brain network in drug‐naive patients with obsessive–compulsive disorder. Brain Behav. 8, e00963 (2018).
Collin, G. et al. Brain network analysis reveals affected connectome structure in bipolar I disorder. Hum. brain Mapp. 37, 122–134 (2016).
Wang, Y. et al. Disrupted resting-state functional connectivity in nonmedicated bipolar disorder. Radiology 280, 529–536 (2016).
Liu, W. et al. Abnormal degree centrality of functional hubs associated with negative coping in older Chinese adults who lost their only child. Biol. Psychol. 112, 46–55 (2015).
Tian, T., Qin, W., Liu, B., Jiang, T. & Yu, C. Functional connectivity in healthy subjects is nonlinearly modulated by the COMT and DRD2 polymorphisms in a functional system-dependent manner. J. Neurosci. 33, 17519–17526 (2013).
Zhang, X. et al. The catechol-o-methyltransferase Val 158 Met polymorphism modulates the intrinsic functional network centrality of the parahippocampal cortex in healthy subjects. Sci. Rep. 5, 10105 (2015).
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. Eurasia J. Math. Sci. Technol. Educ. 13, 3379–3392 (2017).
Yan, C. & Zang, Y. DPARSF: a MATLAB toolbox for” pipeline” data analysis of resting-state fMRI. Front. Syst. Neurosci. 4, 13 (2010).
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., Craddock, R. C., Zuo, X.-N., Zang, Y.-F. & Milham, M. P. Standardizing the intrinsic brain: towards robust measurement of inter-individual variation in 1000 functional connectomes. Neuroimage 80, 246–262 (2013).
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).
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).
Hong, S. B. et al. COMT genotype affects brain white matter pathways in attention-deficit/hyperactivity disorder. Hum. Brain Mapp. 36, 367–377 (2015).
Pietrzak, R. H., Goldstein, R. B., Southwick, S. M. & Grant, B. F. Prevalence and Axis I comorbidity of full and partial posttraumatic stress disorder in the United States: results from Wave 2 of the National Epidemiologic Survey on Alcohol and Related Conditions. J. Anxiety Disord. 25, 456–465 (2011).
Boscarino, J. A. Posttraumatic stress disorder and physical illness: results from clinical and epidemiologic studies. Ann. N. Y. Acad. Sci. 1032, 141–153 (2004).
Duan, Z. et al. Assessment of functional tag single nucleotide polymorphisms within the DRD2 gene as risk factors for post-traumatic stress disorder in the Han Chinese population. J. Affect Disord. 188, 210–217 (2015).
Wolf, E. J. et al. The dopamine D3 receptor gene and posttraumatic stress disorder. J. Trauma Stress 27, 379–387 (2014).
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).
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. Neuroimaging 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).
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).
Gogtay, N. et al. Dynamic mapping of human cortical development during childhood through early adulthood. Proc. Natl Acad. Sci. USA 101, 8174–8179 (2004).
Pechtel, P. & Pizzagalli, D. A. Effects of early life stress on cognitive and affective function: an integrated review of human literature. Psychopharmacology (Berl.) 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).
Williams, G. V. & Goldman-Rakic, P. S. Modulation of memory fields by dopamine D1 receptors in prefrontal cortex. Nature 376, 572–575 (1995).
Vijayraghavan, S., Wang, M., Birnbaum, S. G., Williams, G. V. & Arnsten, A. F. Inverted-U dopamine D1 receptor actions on prefrontal neurons engaged in working memory. Nat. Neurosci. 10, 376–384 (2007).
Prata, D. P. et al. Opposite effects of catechol-O-methyltransferase Val158Met on cortical function in healthy subjects and patients with schizophrenia. Biol. Psychiatry 65, 473–480 (2009).
Di Giorgio, A. et al. Catechol-O-methyltransferase Val 158 Met association with parahippocampal physiology during memory encoding in schizophrenia. Psychol. Med. 41, 1721–1731 (2011).
Fallon, S., Williams-Gray, C., Barker, R., Owen, A. & Hampshire, A. Prefrontal dopamine levels determine the balance between cognitive stability and flexibility. Cereb. Cortex 23, 361–369 (2012).
Chen, C. C. et al. BDNF Val66Met polymorphism on functional MRI during n-back working memory tasks. Medicine (Baltim.) 94, e1586 (2015).
Nieoullon, A. & Coquerel, A. Dopamine: a key regulator to adapt action, emotion, motivation and cognition. Curr. Opin. Neurol. 16(Suppl 2), S3–S9 (2003).
Bandelow, B. et al. Biological markers for anxiety disorders, OCD and PTSD: a consensus statement. Part II: neurochemistry, neurophysiology and neurocognition. World J. Biol. Psychiatry 18, 162–214 (2017).
Hoexter, M. Q. et al. Higher striatal dopamine transporter density in PTSD: an in vivo SPECT study with [99m Tc] TRODAT-1. Psychopharmacology 224, 337–345 (2012).
Zhang, L. et al. PTSD risk is associated with BDNF Val66Met and BDNF overexpression. Mol. psychiatry 19, 8 (2014).
Lee, H. J. et al. No association between the brain‐derived neurotrophic factor gene Val66Met polymorphism and post‐traumatic stress disorder. Stress Health. J. Int. Soc. Investig. Stress 22, 115–119 (2006).
Frielingsdorf, H. et al. Variant burain-derived neurotrophic factor Val66Met endophenotypes: implications for posttraumatic stress disorder. Ann. N. Y. Acad. Sci. 1208, 150–157 (2010).
Rakofsky, J. J., Ressler, K. J. & Dunlop, B. W. BDNF function as a potential mediator of bipolar disorder and post-traumatic stress disorder comorbidity. Mol. Psychiatry 17, 22–35 (2012).
Braun, K. The prefrontal-limbic system: development, neuroanatomy, function, and implications for socioemotional development. Clin. Perinatol. 38, 685–702 (2011).
Fitzgerald, J. M., DiGangi, J. A. & Phan, K. L. Functional neuroanatomy of emotion and its regulation in PTSD. Harv. Rev. Psychiatry 26, 116–128 (2018).
Zhao, F. et al. Network-dependent modulation of COMT and DRD2 polymorphisms in healthy young adults. Sci. Rep. 5, 17996 (2015).
Cui, L. et al. Spontaneous regional brain activity in healthy individuals is nonlinearly modulated by the interaction of ZNF804A rs1344706 and COMT rs4680 polymorphisms. Neurosci. Bull. 35, 735–742 (2019).
Gleich, T. et al. Prefrontal and striatal glutamate differently relate to striatal dopamine: potential regulatory mechanisms of striatal presynaptic dopamine function? J. Neurosci. 35, 9615–9621 (2015).
Goto, Y. & Grace, A. A. Dopamine-dependent interactions between limbic and prefrontal cortical plasticity in the nucleus accumbens: disruption by cocaine sensitization. Neuron 47, 255–266 (2005).
Set, E. et al. Dissociable contribution of prefrontal and striatal dopaminergic genes to learning in economic games. Proc. Natl Acad. Sci. USA 111, 9615–9620 (2014).
Wang, C. et al. Multilocus genetic profile in dopaminergic pathway modulates the striatum and working memory. Sci. Rep. 8, 5372 (2018).
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).
Nievergelt, C. M. et al. Largest genome-wide association study for PTSD identifies genetic risk loci in European and African ancestries and implicates novel biological pathways. bioRxiv 458562 (2018).
Acknowledgements
We acknowledge the support of grants from the National Nature Science Foundation of China [numbers 81671672 and 81301209 to R.Q. and 81801678 to Y.L.]; the Jiangsu Provincial Medical Youth Talent [numbers QNRC2016888 to R.Q. and QNRC2016207 to Y.L.]; the Chinese Key Grant [number BWS11J063 to G.M.L.]; and the Youth Nature Science Foundation of Jiangsu Province [number BK20170223 to Y.L.]. We thank Professor Paul Thompson of University of Southern California, for his helpful advice and editing contributions for this paper.
Author information
Authors and Affiliations
Corresponding authors
Ethics declarations
Conflict of interest
The authors declare that they have no conflict of interest.
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. Effects of COMT rs4680 and BDNF rs6265 polymorphisms on brain degree centrality in Han Chinese adults who lost their only child. Transl Psychiatry 10, 46 (2020). https://doi.org/10.1038/s41398-020-0728-7
Received:
Revised:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41398-020-0728-7