Oxytocin receptor gene, childhood maltreatment and borderline personality disorder features among male inmates in China
BMC Psychiatry volume 20, Article number: 332 (2020)
Borderline personality disorder (BPD) is caused by a variety of biological and environmental factors. Accumulating evidence suggests that childhood maltreatment is a risk environmental factor in the development of BPD, but research on the genetic pathology of BPD is still in its early stages, and very little is known about the oxytocin receptor (OXTR) gene. The purpose of this study is to further explore the interactive effects between OXTR gene polymorphisms and childhood maltreatment on BPD risk.
Among the 1804 Chinese Han male inmates, 765 inmates who had BPD or antisocial personality disorder (ASPD) or highly impulsive or violent crime were considered as high-risk inmates and included in this study. Childhood maltreatment, BPD, antisocial personality disorder (ASPD) and impulsivity were measured by self-reported questionnaires. Peripheral venous blood was collected for the genotype test.
Analyses revealed that the BP group (inmates with BPD features) had higher rs53576 AA genotype frequency and rs237987 AA genotype frequency than the non-BP group, while the statistical significances were lost after Bonferroni correction. Total childhood maltreatment score, emotional abuse and neglect could positively predict BPD risk. Among the high-risk samples, rs53576 GG genotype carriers had higher BPD scores at higher levels of physical abuse and sexual abuse and had lower BPD scores at lower levels of physical abuse and sexual abuse.
The findings suggest that the interaction between OXTR gene variations and childhood maltreatment is an important mechanism for the development of BPD. The moderating role of the OXTR gene provides evidence for gene plasticity.
Borderline personality disorder (BPD) is characterized by a pervasive pattern of emotional instability, interpersonal tension, high impulsiveness, identity disorder and cognitive impairment . BPD has been estimated to have 10% mortality caused by suicide . Therefore, BPD imposes enormous economic and psychological burdens, highlighting the necessity of elucidating the pathological mechanism of BPD. Decades of research demonstrate that BPD is caused by a variety of biological and environmental factors, and the interaction between these two factors has been emphasized recently .
Twin studies suggested that 35–67% of the variance in BPD risk is attributable to genetic factors [4, 5]. Great efforts have been made to elucidate the possible genetic background of BPD, and the current polymorphism studies mainly focused on the serotonin receptor (HTR), tryptophan hydroxylase (TPH), monoamine oxidase A (MAOA), catechol-O-methyltransferase (COMT) and dopamine transporter (DAT) genes [6,7,8,9]. However, very few studies have investigated the association between oxytocin receptor (OXTR) gene and BPD.
Oxytocin is a peptide hormone synthesized in the paraventricular and the supraoptic nucleus of the hypothalamus. Oxytocinergic neurotransmission is mediated by the oxytocin receptors. The OXTR gene located on chromosome 3p25 is associated with encoding and activity of oxytocin receptors . Several single nucleotide polymorphisms (SNPs) in the OXTR gene have been proven to be related to social cognition and behavior as well as greater susceptibility to psychiatric disorders. The OXTR SNPs rs53576 and rs2254298 were frequently studied in human [11, 12]. Specifically, previous studies found that OXTR rs53576 was associated with prosocial behavior including positive parenting  and better empathic ability , as well as psychiatry disorders such as schizophrenia  and depression . Similar to the SNP rs53576, numerous studies have demonstrated that the SNP rs2254298 was implicated in infant attachment security , higher social cognition , depression , autism spectrum disorder (ASD) . However, the associations between the two SNPs and BPD remain to be explored.
Although other OXTR SNPs were relatively less studied, evidence has accumulated suggesting that OXTR gene polymorphisms are associated with the susceptibility to psychopathology. For instance, a previous study found that the haplotype ATTGA (rs237897-rs13316193-rs237889-rs2254298-rs2268494) was associated with ASD . In addition, six SNPs in the OXTR gene including rs1042778, rs2268490, rs13316193, rs2254298, rs237889 and rs13311693 were linked to callous-unemotional (CU) traits . Malik et al. found that rs6770632 and rs1042778 might be associated with aggressive behaviors in females and males, respectively . In the present study, we focused on the eight SNPs in the OXTR gene (rs237897, rs13316193, rs237889, rs2254298, rs2268494, rs1042778, rs53576 and rs6770632), in order to explore the relationship between OXTR genetic variants and BPD.
As mentioned above, OXTR SNPs not only affect the prosocial behavior of individuals, but also affect the development of psychiatry disorders. The contradictory effect may be attributable to environmental factors, which means the interaction between the gene and the environment (G × E) may play a role in the development of BPD. The function of OXTR gene may vary according to social stimulus which conformed to the social sensitivity hypothesis [24, 25]. Specifically, OXTR gene is related to increased sensitivity to the surroundings, both in favorable or negative contexts. A recent study confirmed this hypothesis that individual who carried OXTR rs53576 A-allele had higher levels of BPD symptoms under negative family conditions and lower levels under positive conditions . Considering that the function of OXTR is environment-based, it is necessary to explore the environmental factors of BPD.
Regarding the influence of environmental factors on BPD, previous studies had investigated the associations between family environment , childhood maltreatment , early life stress  and BPD. In the present study, we sought to shed more light on childhood maltreatment and its interaction with OXTR gene. Several lines of evidence show that there are interactions between the OXTR gene polymorphisms and childhood maltreatment; for instance, OXTR rs53576 G-allele carriers who had experienced maltreatment in childhood were observed to increase conduct problems in adolescence , and maltreated individual who carried the rs237885 TT genotype in OXTR gene had an increased risk of aggression . In addition, the interaction between the OXTR gene and childhood maltreatment could affect the function or structure of the brain, especially the amygdala , which is responsible for emotional regulation what BPD patients lack. However, research on the interaction between the OXTR gene and childhood maltreatment for the development of BPD is very scarce [26, 33].
Reviewing researches on the interaction between OXTR gene and childhood maltreatment, we found that rs53576 has emerged as an important SNP in the understanding of the G × E interactions. Previous studies provide evidence that the GG genotype carriers are assumed to show higher trust-related behaviors and empathic concern in comparison to A allele carriers [34, 35]. However, in face of adversity such as childhood abuse, the GG genotype carriers may be more vulnerable to be affected by negative early-life experiences. For instance, Bradley et al. found that the GG genotype carriers exposed to childhood maltreatment exhibited emotional dysregulation and disorganized attachments styles . Hostinar et al. showed that GG carriers perceived significantly lower social support and reported higher levels of internalizing symptoms than A-carriers among maltreated adolescents .
Previous study suggested that most inmates had a personality disorder when they committed the crime. BPD and antisocial personality disorder (ASPD) are the most common personality disorders in prison setting . Several studies have provided strong evidence that ASPD and BPD are risk factors of violence [39,40,41], and the association between the personality disorders and violence may be a result of impulsiveness . Therefore, individuals with BPD and ASPD features, violent criminals and high impulsive prisoners are all high-risk groups in prison, which may lead to a great deal of hardship for prison management, especially among male inmates. Although the prevalence of BPD in women is higher than that in men , the crime rate of men is higher than that of women . Considering the absence of such study on this group in China, male high-risk inmates were chosen as our study subjects.
To explore the role of the OXTR gene in contributing risk to BPD features, and whether these polymorphisms which may be related to the variations of OXTR gene have different effects on BPD features depends on the levels of childhood maltreatment in male high-risk inmates, we sought to test three hypotheses: (1) participants with borderline personality features would show differences in OXTR genotype frequency compared with those without borderline personality features; (2) childhood maltreatment and its dimensions could positively predict BPD among male high-risk inmates in China; and (3) the OXTR gene polymorphisms would interact with childhood maltreatment or its dimensions to predict BPD risk.
A total of 1804 male inmates from an adult male prison in Jiangsu Province, China initially participated in the present study, and all participants were Chinese Han individuals. The study was approved by the Institutional Review Board of Nanjing Brain Hospital, and a consent form for this study was prepared. We carried out the study following the principle of the Helsinki Declaration. For the compensation for the participants, some lectures, psychoeducation about mental health issues for inmates and some didactic lessons and supervise program about mental illness for prison clinical faculty and staff were regularly arranged for their support of the research. After obtaining informed consent to conduct the self-report survey, all of them were investigated by a structure questionnaire asking about sociodemographic characteristics, personality disorders, history of childhood maltreatment and impulsivity. Considering that the present study aimed to investigate the high-risk population in prison, participants who met one of the following conditions would be included: (1) BPD or ASPD score in the Personality Diagnostic Questionnaire-4 (PDQ-4) ≥ 4 ; (2) Barratt Impulsiveness Scale-11 (BIS-11) score ≥ 62 ; (3) history of violent crime, such as intentional homicide, intentional injury, and robbery. Among the 1804 participants, 1050 (58.2%) inmates were considered as high-risk population. 765 (72.86%) valid blood samples were obtained from the 1050 inmates after obtaining informed consent for blood collection. Finally, 765 subjects with questionnaire data and blood samples were included in the study. Those who had previously been diagnosed with chronic heart, liver and kidney diseases; had a history of nervous system diseases or mental disorders; or had a long history of medication use were excluded.
Among the 765 inmates with high risk, 342 (44.71%) inmates met the criteria for borderline personality features (BP group), while 423 (55.29%) inmates did not meet the criteria (non-BP group), i.e., a score of 4 points on the BPD subscale of the PDQ-4. The average age of the high-risk subjects was 33.50 years old (SD = 9.01). Most subjects had a low level of education, of which 33.7% (n = 258) had attended school for less than or equal to 6 years, 62.2% (n = 476) for 7–12 years, and the remaining 4.1% (n = 31) for 13 years or more. Nearly half of the subjects had relatively stable partner relationships, as 46.4% (n = 355) of them were in married or cohabiting, 38.7% (n = 296) of them were unmarried and the remaining 14.9% (n = 114) were divorced or widowed. Among the high-risk sample, 57.6% (n = 441) had children.
Personality diagnostic Questionnaire-4 (PDQ-4)
The PDQ-4 was administered to assess BPD and ASPD which is a further revision to PDQ based on DSM-IV . This self-reported questionnaire consisted of 85 items, which are used to assess 10 types of personality disorders (paranoid, schizophrenic, split, hysteric, narcissistic, borderline, antisocial, avoidant, dependent and compulsive). In this study, the BPD subscale (9 items) of the PDQ was used to assess the borderline personality traits of the subjects instead of clinical diagnosis, with higher scores reflecting correspondingly higher levels of BPD. The ASPD subscale (8 items) was used to screen whether the inmates had the risk of ASPD. PDQ-4 demonstrated good psychometric properties when applied to clinical samples and general population [47, 48]. The internal consistency as reflected by Cronbach’s alpha for the 85 items in the present study was 0.919, and the internal consistency for BPD and ASPD subscales was 0.689 and 0.675, respectively.
Childhood trauma questionnaire short form (CTQ-SF)
The CTQ-SF, compiled by Bernstein DP and Flink L , was administered to assess five categories of childhood abuse and neglect including emotional abuse, physical abuse, sexual abuse, emotional neglect and physical neglect. This instrument is a 28-item retrospective self-reported questionnaire, and the items follow a Likert-type response style from 1 to 5, which is organized to reflect the frequency of abusive experiences (never true, rarely true, sometimes true, often true, very often true). Each dimension has 5 items, and the remaining 3 items are validity evaluation, so each dimension is scored between 5 ~ 25 and the total score is 25 ~ 125. The CTQ and its dimensions (other than physical neglect) had good reliability and validity in Chinese population . The internal consistency as reflected by Cronbach’s alpha for the 28 items in the present study was 0.897, and the internal consistency for each subscale was 0.801 (emotional abuse), 0.875 (physical abuse), 0.878 (sexual abuse), 0.843 (emotional neglect), and 0.628 (physical neglect), respectively.
Barratt impulsiveness Scale-11 (BIS-11)
The BIS-11 was administered to assess impulsivity which was originally developed by Barratt in 1959. This self-reported questionnaire consists of 26 items, which is answered on a 4-point Likert scale ranging from 1 to 4 (rarely, occasionally, often, almost always). Three types of impulsivity, including attentional impulsiveness, motor impulsiveness, and non-planning impulsiveness, are measured, with the higher scores of each dimension reflecting corresponding higher levels of such impulsiveness. In this study, the total score of BIS-11 was used to screen impulsivity of subjects in prison. According to the previous studies in Chinese population, the BIS-11 and its dimensions demonstrated good reliability and validity . The internal consistency as reflected by Cronbach’s alpha for the 26 items in the present study was 0.743.
SNPs selection, DNA extraction and genotyping
According to the previous studies [11,12,13,14,15,16,17,18,19,20,21,22,23], the SNPs which were found to be related to the development of psychiatric disorders were selected from the 1000 genome database of Han Chinese in Beijing, China (CHB) (http://www.internationalgenome.org/). We focused on eight SNPs in the OXTR gene, including rs1042778, rs13316193, rs2254298, rs2268494, rs237889, rs237897, rs53576, rs6770632.
Venous blood samples were collected from 765 subjects (342 subjects in BP group and 423 subjects in non-BP group). Genetic DNA was extracted using Sunshinebio TM blood genomic DNA Extraction Kit (Sunshinebio, Nanjing, China) following the manufacturer’s instructions. The extracted DNA was diluted to a concentration of 50 ng/μL. Genotyping of the OXTR SNPs applied multiple PCR technology and high-throughput sequencing technology. The primers showed in Additional file 1: Table S1. Multiple PCR amplification was performed using PERKIN ELMER Gene Amp PCR system 9600 (Applied Biosystems, Shanghai, China). Different samples were distinguished by different Barcode primers, and the high flux sequencing of the amplified products was performed using ABI PRISM 377 DNA Sequencer (Applied Biosystem, Shanghai, China). Several samples failed to be genotyped, but the call rates for all SNPs almost reached 95% (Table 2).
Statistical analyses were performed using Statistical Product and Service Solutions (SPSS) 22.0. Linkage disequilibrium (LD) and haplotype construction were performed using Haploview 4.1 software. First, descriptive statistics were used to present the demographic characteristics of the subjects and the distribution of alleles of all investigated SNPs. Second, the Hardy-Weinberg equilibrium (HWE) for genotypes was assessed by a chi-squared goodness-of-fit test to examine whether the sample was a genetic balance population. Third, the differences in demographic characteristics and genotypes between the BP group and the non-BP group were examined by performing a series of t-tests and chi-squared goodness-of-fit tests. Fourth, one-way analysis of variance (ANOVA) and chi-squared goodness-of-fit tests were used to test for differences in demographic characteristics stratified by genotypes. Fifth, association analysis between childhood maltreatment subtypes and BPD level was performed with a multiple linear regression. Lastly, the interactions between OXTR gene polymorphisms (rs53565: GG/AG + AA; rs237897: GG/AG + AA) and childhood maltreatment were analyzed using hierarchical linear regressions.
All reported p values were two-sided. The nominal significance threshold was set at 0.05 and Bonferroni corrected p-value threshold for multiple testing was set at 0.00625 (0.05/8 SNPs) or 0.0125 (0.05/4 haplotypes).
According to the BPD scores, the high-risk subjects were divided into a BP group and a non-BP group, 342 (44.71%) cases and 423 (55.29%) cases respectively. As presented in Table 1, groups were similar in age, education, marital status, and reproductive status (p > 0.05). Subjects in the BP group obtained higher scores on the overall CTQ as well as its five subtypes than those in the non-BP group (p < 0.001).
Association analysis between OXTR SNPs and BP
Table 2 presents the basic information about the eight OXTR SNPs which were genotyped in the present study. The minor allele frequencies (MAF) of these eight SNPs were similar to those in the 1000 genome database of CHB. The genotypes in the high-risk inmates conformed to Hardy-Weinberg equilibrium (HWE) (p > 0.05). HWE tests also conducted in the subgroups. The results indicated that, except OXTR rs1042778 in the non-BP group (χ2 = 10.396, p = 0.001), all other SNPs conformed to HWE in both groups (χ2 = 0.019–2.197, p > 0.05). Considering the SNP rs1042778 did not deviate from HWE in the whole sample of high-risk inmates (Table 2), such deviation in the non-BP group may be due to the reduction of sample size.
As showed in Table 3, the distributions of age showed significant differences in the genotype frequencies of rs2268494 (p < 0.007). Chi-squared tests showed that there were significant differences in the genotype frequencies of rs237889 among different levels of education (p < 0.05). In addition, significant differences in the genotype frequencies of rs53576 were observed among different marital status (p < 0.05). However, after Bonferroni correction, only the relationship between age and rs2268494 remained statistically significant. Apart from the results mentioned above, age (F = 0.047–1.134, p > 0.05), education (χ2 = 0.728–7.264, p > 0.05), marital status (χ2 = 3.605–8.685, p > 0.05), and reproductive status (χ2 = 0.081–5.495, p > 0.05) were not found to be statistically significant in the genotype frequencies of other SNPs.
As shown in Table 4, there were significant differences in the genotype frequencies of rs53576 and rs237897 between the two groups (p < 0.05), while the statistical significances were lost after Bonferroni correction. However, the other six markers, namely, rs13316193 (χ2 = 1.581, p = 0.454), rs2254298 (χ2 = 1.231, p = 0.540), rs2268494 (χ2 = 3.016, p = 0.221), rs237889 (χ2 = 3.812, p = 0.149), rs1042778 (χ2 = 5.063, p = 0.080), and rs6770632 (χ2 = 1.068, p = 0.477), were not found to be statistically significant between the two groups.
Figure 1 presented the pairwise linkage disequilibrium (LD) results of the eight SNPs, and a haplotype block consisting of three SNPs (rs2268494, rs2254298, rs237897) was identified. There were no significant differences in the frequencies of these haplotypes between BP group and non-BP group (Table 5).
Association analysis between childhood maltreatment and BPD
The impact of childhood maltreatment on BPD was examined by a linear regression analysis, in which BPD score was the dependent variable and total score of childhood maltreatment was the independent variable. The analysis revealed that childhood maltreatment was a significant positive predictor of BPD (R2 = 0.102, β = 0.322, S. E = 0.004 t = 9.384, p < 0.001). Because the total CTQ score had a significant independent effect on BPD, the five subtypes of childhood maltreatment were all entered into the whole model to determine whether any of the subtypes had an independent effect on BPD. The method of stepwise was used to control the multiple linear relationship between the independent variables. The results indicated that emotional abuse (β = 0.201, S. E = 0.021, t = 5.575, p < 0.001), and emotional neglect (β = 0.205, S. E = 0.014, t = 5.762, p < 0.001) could significantly predict BPD (R2 = 0.106, F = 46.444,p < 0.001), while physical abuse (t = 0.514, p = 0.607), sexual abuse (t = 0.542, p = 0.588) and physical neglect (t = 1.056, p = 0.291) were excluded from the multivariate linear regression model.
Gene-environment interactions with respect to BPD
According to the results in 3.2, the two SNPs (rs53576 and rs237897) showed significant differences in genotype frequencies between the BP and non-BP groups, while these differences were not significant after Bonferroni correction. This indicates that the two SNPs may be related to the OXTR genetic variants, which further lead to the development of BPD. Therefore, we focus on the moderation roles of the two SNPs (rs237897 and rs53576) on the relationship between childhood maltreatment and BPD.
Interaction between OXTR rs53576 and childhood maltreatment
To examine the interaction between OXTR rs53576 and childhood maltreatment on BPD features in the high-risk inmates, we constructed a hierarchical linear regression. Genotype and maltreatment were entered on the first step, and the genotype × maltreatment interaction term was entered on the second step. The moderating role of OXTR rs53576 genotype on the relationships between physical abuse (R2 = 0.066, t = 2.723, p = 0.007), sexual abuse (R2 = 0.053, t = 2.914, p = 0.004) and BPD was significant, while childhood maltreatment total score (R2 = 0.105, t = 1.575, p = 0.116), emotional abuse (R2 = 0.077, t = 1.931, p = 0.054), emotional neglect (R2 = 0.071, t = − 1.016, p = 0.310) and physical neglect (R2 = 0.076, t = 0.990, p = 0.323) all failed to interact with genotype to predict BPD scores. As shown in Fig. 2, when subjects had lower physical abuse scores, those who carried the GG genotype had a lower risk of BPD than those who carried the A allele. While, when subjects had higher physical abuse scores, those who carried GG genotype had a higher risk of BPD than those who carried A allele. The interaction between OXTR rs53576 genotype and sexual abuse was similar. As shown in Fig. 3, at lower levels of sexual abuse, BPD scores were lower among those with GG genotype, while at higher levels of sexual abuse, BPD scores were lower among those with A allele.
The interactions between OXTR rs53576 and childhood abuse were also examined in the BP and non-BP groups. The results showed that neither childhood maltreatment total score (R2 = 0.033, t = − 0.216, p = 0.829; R2 = 0.040, t = 0.935, p = 0.351, respectively), nor emotional abuse (R2 = 0.021, t = 0.591, p = 0.555; R2 = 0.057, t = 1.318, p = 0.188, respectively), physical abuse (R2 = 0.015, t = − 0.311, p = 0.756; R2 = 0.043, t = 1.608, p = 0.109, respectively), sexual abuse (R2 = 0.003, t = − 0.339, p = 0.735; R2 = 0.030, t = 1.683, p = 0.093, respectively), emotional neglect (R2 = 0.035, t = 0.568, p = 0.570; R2 = 0.003, t = − 0.289, p = 0.773, respectively), physical neglect (R2 = 0.040, t = − 0.914, p = 0.361; R2 = 0.013, t = 0.572, p = 0.568, respectively) interacted with genotype to predict BPD scores in both group.
Interaction between OXTR rs237897 and childhood maltreatment
A two-way interaction model was conducted to examine the moderating effect of OXTR rs237897 on the relationship between childhood maltreatment and BPD features in the whole high-risk sample. The analysis revealed that the interactive effect of OXTR gene by childhood maltreatment total score was non-significant for BPD in high-risk inmates (R2 = 0.097, t = 0.341, p = 0.734). Similar non-significant results were obtained for emotional abuse (R2 = 0.070, t = − 0.402, p = 0.688), physical abuse (R2 = 0.056, t = 0.475, p = 0.635), sexual abuse (R2 = 0.040, t = 0.765, p = 0.444), emotional neglect (R2 = 0.066, t = 0.658, p = 0.511) and physical neglect (R2 = 0.072, t = 1.525, p = 0.128).
To examine the moderating effect of OXTR rs237897 on the relationship between childhood maltreatment and BPD features in the BP and non-BP groups respectively, we also constructed two-way interaction models. There were no significant interactions between genotype and childhood maltreatment total score (R2 = 0.029, t = 1.244, p = 0.214; R2 = 0.047, t = 1.608, p = 0.109, respectively), emotional abuse (R2 = 0.016, t = 0.772, p = 0.441; R2 = 0.060, t = 1.271, p = 0.204, respectively), physical abuse (R2 = 0.011, t = 0.471, p = 0.638; R2 = 0.044, t = 1.255, p = 0.210, respectively), sexual abuse (R2 = − 0.004, t = 1.114, p = 0.266; R2 = 0.030, t = 1.331, p = 0.184, respectively), emotional neglect (R2 = 0.032, t = 1.457, p = 0.146; R2 = 0.005, t = 0.199, p = 0.842, respectively), physical neglect (R2 = 0.034, t = 1.772, p = 0.077; R2 = 0.015, t = 1.550, p = 0.122, respectively) in both groups.
This study investigated genetic (OXTR) and environmental (childhood maltreatment) factors of BPD and highlighted G × E interactions relevant to BPD risk. In terms of genetic test results, the BP group had an increased rs53576 AA genotype frequency compared with the non-BP group. Although the difference did not survive after Bonferroni correction, this finding indicates that the high AA genotype frequency may increase the risk of BPD. This result is inconsistent with previous finding that the genotype groups for OXTR did not demonstrate differences in BPD scores . However, these discrepancies may be due to the different analysis methods used to explore the relationship between rs53576 and BPD. Another possible explanation is ethnicity-related differences in allele frequency distribution. The frequency distribution of rs53576 alleles in this Asian population is opposite to that in Caucasian population, which indicates that the OXTR genetic variants in the two races are different. In the SNP rs237987, the higher AA genotype frequencies were found in BP group, while the difference was not significant after Bonferroni correction. The association between the specific oxytocin receptor marker rs237987 and BPD susceptibility has not previously been reported, so the results need to be verified in a larger and more representative sample.
With respect to the environmental aspect, the results showed that childhood maltreatment could positively predict BPD risk in the present study. This is not surprising given that evidence has accumulated suggesting that early childhood maltreatment is a risk factor for BPD. For instance, a retrospective study suggested that children who had been maltreated had an increased risk of presenting borderline personality features . Another study in adult reported the association between experiencing traumatic events in childhood and an increased clinical severity of BPD in adulthood .
However, the present study only found that CTQ scores, emotional abuse and neglect could positively predict BPD risk. This may be due to the fact that a certain type of traumatic experience may increase the risk of specific mental illness in the trajectory of susceptible individuals from childhood to adulthood. For instance, a clinical study has found that emotional and sexual abuse may lead to BPD or BPD and attention deficit hyperactivity disorder (ADHD) comorbidity, while the history of physical abuse was related to the persistence of ADHD in adulthood . From another perspective, the symptoms of BPD are heterogeneous, with multiple phenotypes, which may be caused by different factors in the development of individuals. A previous study suggested that affective instability was linked with childhood sexual abuse, while identity disturbance was associated with non-sexual developmental trauma among BPD patients .
The present study found OXTR gene×childhood maltreatment interactions only in the high-risk inmates, while no significant interaction was found in BP group and non-BP group. This indicates there are no significant gene×maltreatment×group three-way interactions, and the pattern of interactions in this study are only applicable to high-risk prison inmates. In the OXTR rs53576, the GG carriers were strongly affected whether the physical abuse or sexual abuse was low or high, which conforms to the hypothesis of gene plasticity . These results are consistent with the prior literature showing that GG carriers may be more sensitive to social cues and more prosocial than A-allele carriers, while this attunement for social cues may be a disadvantage under adverse circumstances such as early childhood maltreatment . In contrast, a recent study reported that A-allele carriers had higher levels of BPD symptoms under negative family conditions and lower levels under positive conditions, GG genotype carriers had average levels of BPD symptoms regardless of their family quality . Considering there are ethnic-related differences in the genotype frequency of rs53576, genetic susceptibility may also vary among different races. However, previous studies in Caucasian , African American  and multi-ethnic groups  found that carriers of at least one rs53576 G-allele were more susceptible to negative life experiences. In addition, a meta-analysis suggested that GG homozygotes had higher general sociality than A-allele carriers both in Caucasians and Asians . These results indicate that GG genotype carriers of OXTR rs53576 might be more sensitivity to environment, which was not affected by ethnicity. This may partially explain the findings of the present study.
Another possibility is that the contradictory results of genetic susceptibility might be attributable to gender differences. A study in low-income children of different ethnicities found that girls who carried the A allele were more susceptible than those with GG genotype, while boys with the GG genotype were more susceptible than those with A allele . A recent study confirmed that the A allele carriers were more susceptible to childhood maltreatment in females . Thus, it may mean there are different patterns of OXTR gene and childhood maltreatment interactions in different genders. Future research needs to clarify the exact mechanism.
It is noted that the present study found the total CTQ scores, emotional abuse, emotional neglect and physical neglect did not show significant interactions, which could significantly predict BPD scores. In contrast, physical and sexual abuse interacted with genotype to predict BPD scores, while these two types of abuse did not significantly predict BPD scores. These findings indicate that total CTQ scores, emotional abuse and emotional neglect have independent predictive effects on BPD, while the effects of physical and sexual abuse on BPD depend on the genetic characteristics of individuals. Another plausible explanation is that OXTR rs53576 interacts with physical and sexual abuse to influence OXTR transcription and expression. There is evidence that physical exposure, including physical abuse and sexual abuse, showed strongest association with methylated variation .
The explanation for the significant interactions between gene and environment should be cautious. The prior literature reported that adolescents whose mothers had BPD experienced more maltreatment compared with those whose mothers were normal , which revealed that the intergenerational transmission of BPD may be through childhood maltreatment. Another study showed that mothers’ maltreatment history was a more strongly prediction for mother-infant attachment disorganization score among mothers with more plasticity alleles of OXTR gene , which suggested that maternal maltreatment history would influence offspring depending on the maternal genetic characteristics. Hence, both genetic and environmental factors are contributed to the development of BPD. Whether adverse environment lead to changes in biomarkers or whether genetic basis increases susceptibility to adverse environment remains future research to explore.
There are several strengths in the study. First, eight OXTR SNPs were genotyped based on the previous literature, which is helpful to understand the effect of OXTR genetic variants on BPD. In addition, to our knowledge, this is the first study to demonstrate interactive effects between OXTR polymorphisms and childhood maltreatment subtypes on BPD in Chinese male high-risk inmates. Given the characteristics and high risk for BPD, this population is worth investigating. The findings that the patterns of G × E interactions in prison sample contribute both to prison management and to treatment or even prevention of BPD.
This study also has several limitations that need to be acknowledged. The current sample consists of male high-risk inmates, a population that is at risk for BPD. Thus, it is possible that the results would differ in the general population. Given that the different effects of oxytocin in different genders , the findings also cannot be generalized to females. Participants were only assessed with self-reported questionnaires. Although the BPD and ASPD subscales in PDQ-4 as well as CTQ-SF had a good reliability in this study, self-reporting may lead to potential bias due to their own memory of events. In addition, the subjects who had higher levels of BPD were likely to exaggerate the adverse events they experienced in early life, which would inflate the association between childhood maltreatment and BPD. The onset and duration of childhood maltreatment were not considered, so the interpretation and application of the results should be cautious.
In summary, the findings of the present study indicated that OXTR gene variations and childhood maltreatment are both risk factors for the development of BPD and that OXTR polymorphisms interact with childhood maltreatment subtypes to predict BPD scores. These interaction findings reinforced the view that OXTR is a plasticity gene that is associated with increased sensitivity to the environment, regardless of negative or positive [56, 62]. The results may contribute to understand the pathological mechanism of BPD, further guiding the early identification of personality disorders and the selection of interventions.
Availability of data and materials
The datasets used and analyzed during the current study are available from the corresponding author on reasonable request.
Borderline personality disorder
Antisocial personality disorder
Monoamine oxidase A
Single nucleotide polymorphism
Autism Spectrum Disorder
- G × E:
Interaction between the gene and the environment
Personality Diagnostic Questionnaire-4
Barratt Impulsiveness Scale-11
Childhood Trauma Questionnaire Short Form
Han Chinese in Beijing, China
Statistical Product and Service Solutions
Analysis of variance
Minor allele frequencies
attention deficit hyperactivity disorder
Association, AP. Diagnostic and statistical manual of mental disorders. 2013.
Oldham JM. Borderline Personality Disorder and Suicidality. Am J Psychiatry. 2006;163(1):20–6.
Amad A, Ramoz N, Thomas P, Jardri R, Gorwood P. Genetics of borderline personality disorder: systematic review and proposal of an integrative model. Neurosci Biobehav Rev. 2014;40(40):6–19.
Torgersen S, Czajkowski N, Jacobson K, Reichborn-Kjennerud T, Røysamb E, Neale MC, Kendler KS. Dimensional representations of DSM-IV cluster B personality disorders in a population-based sample of Norwegian twins: a multivariate study. Psychol Med. 2008;38(11):1617–25.
Torgersen S, Myers J, Reichborn-Kjennerud T, Røysamb E, Kubarych T, Kendler K. The heritability of cluster B personality disorders assessed both by personal interview and questionnaire. J Personal Disord. 2012;26(6):848–66.
Joyce P, Stephenson J, Kennedy M, Mulder R, McHugh P. The presence of both serotonin 1A receptor (HTR1A) and dopamine transporter (DAT1) gene variants increase the risk of borderline personality disorder. Front Genet. 2014;4:313.
Perez-Rodriguez M, Weinstein S, New A, Bevilacqua L, Yuan Q, Zhou Z, et al. Tryptophan-hydroxylase 2 haplotype association with borderline personality disorder and aggression in a sample of patients with personality disorders and healthy controls. J Psychiatr Res. 2010;44(15):1075–81.
Kolla N, Vinette S. Monoamine oxidase a in antisocial personality disorder and borderline personality disorder. Curr Behav Neurosci Rep. 2017;4(1):41–8.
Tadić A, Victor A, Başkaya O, von Cube R, Hoch J, Kouti I, et al. Interaction between gene variants of the serotonin transporter promoter region (5-HTTLPR) and catechol O-methyltransferase (COMT) in borderline personality disorder. Am. J. Med. Genet. B Neuropsychiatr. Genet. 2009;150B(4):487–95.
Tops S, Habel U, Radke S. Genetic and epigenetic regulatory mechanisms of the oxytocin receptor gene (OXTR) and the (clinical) implications for social behavior. Horm Behav. 2018;108:84–93.
Feldman R, Monakhov M, Pratt M, Ebstein RP. Oxytocin pathway genes: evolutionary ancient system impacting on human affiliation, sociality, and psychopathology. Biol Psychiatry. 2016;79(3):174–84.
Cataldo I, Azhari A, Lepri B, Esposito G. Oxytocin receptors (OXTR) and early parental care: an interaction that modulates psychiatric disorders. Res Dev Disabil. 2018;82:27–38.
Michalska KJ, Decety J, Liu C, Chen Q, Martz ME, Jacob S, et al. Genetic imaging of the association of oxytocin receptor gene (OXTR) polymorphisms with positive maternal parenting. Front Behav Neurosci. 2014;8:21.
Smith KE, Porges EC, Norman GJ, Connelly JJ, Decety J. Oxytocin receptor gene variation predicts empathic concern and autonomic arousal while perceiving harm to others. Soc Neurosci. 2014;9(1):1–9.
Montag C, Brockmann EM, Bayerl M, Rujescu D, Müller DJ, Gallinat J. Oxytocin and oxytocin receptor gene polymorphisms and risk for schizophrenia: a case-control study. World J Biol Psychiatry. 2013;14(7):500–8.
Thompson SM, Hammen C, Starr LR, Najman JM. Oxytocin receptor gene polymorphism (rs53576) moderates the intergenerational transmission of depression. Psychoneuroendocrinology. 2014;43:11–9.
Chen FS, Barth ME, Johnson SL, Gotlib IH, Johnson SC. Oxytocin receptor (OXTR) polymorphisms and attachment in human infants. Front Psychol. 2011;2:200.
Wade M, Hoffmann TJ, Wigg K, Jenkins JM. Association between the oxytocin receptor (OXTR) gene and children’s social cognition at 18 months. Genes Brain Behav. 2014;13(17):603–10.
Brüne M. Does the oxytocin receptor (OXTR) polymorphism (rs2254298) confer “vulnerability” for psychopathology or “differential susceptibility”? Insights from evolution. BMC Med. 2012;10:38.
LoParo D, Waldman ID. The oxytocin receptor gene (OXTR) is associated with autism spectrum disorder: a meta-analysis. Mol Psychiatry. 2014;20(5):640–6.
Lerer E, Levi S, Salomon S, Darvasi A, Yirmiya N, Ebstein RP. Association between the oxytocin receptor (OXTR) gene and autism: relationship to Vineland adaptive behavior scales and cognition. Mol Psychiatry. 2008;13(10):980–8.
Dadds MR, Moul C, Cauchi A, Dobson-Stone C, Hawes DJ, Brennan J, Urwin R, Ebstein RE. Polymorphisms in the oxytocin receptor gene are associated with the development of psychopathy. Dev Psychopathol. 2014;26(1):21–31.
Malik AI, Zai CC, Berall L, Abu Z, Din F, Nowrouzi B, et al. The role of genetic variants in genes regulating the oxytocin-vasopressin neurohumoral system in childhood-onset aggression. Psychiatr Genet. 2014;24(5):201–10.
Bartz JA, Zaki J, Bolger N, Ochsner KN. Social effects of oxytocin in humans: context and person matter. Trends Cogn Sci. 2011;15(7):301–9.
Kumsta R, Heinrichs M. Oxytocin, stress and social behavior: neurogenetics of the human oxytocin system. Curr Opin Neurobiol. 2013;23(1):11–6.
Hammen C, Bower JE, Cole SW. Oxytocin receptor gene variation and differential susceptibility to family environment in predicting youth borderline symptoms. J Personal Disord. 2015;29(2):177–92.
Mahan R, Kors S, Simmons M, Macfie J. Maternal psychological control, maternal borderline personality disorder, and adolescent borderline features. Personal Disord. 2018;9(4):297–304.
Elices M, Pascual JC, Carmona C, Martín-Blanco A, Feliu-Soler A, et al. Exploring the relation between childhood trauma, temperamental traits and mindfulness in borderline personality disorder. BMC Psychiatry. 2015;15(1):180.
Kim M, Kim J, Park H, Choi S, Oh W, Seok J. Early life stress, resilience and emotional dysregulation in major depressive disorder with comorbid borderline personality disorder. J Affect Disord. 2018;236:113–9.
Andreou D, Comasco E, Åslund C, Nilsson K, Hodgins S. Maltreatment, the oxytocin receptor gene, and conduct problems among male and female teenagers. Front Hum Neurosci. 2018;12:112.
Zhang Y, Wu C, Chang H, Yan Q, Wu L, Yuan S, et al. Genetic variants in oxytocin receptor gene (OXTR) and childhood physical abuse collaborate to modify the risk of aggression in chinese adolescents. J Affect Disord. 2018;229:105–10.
Dannlowski U, Kugel H, Grotegerd D, Redlich R, Opel N, Dohm K, et al. Disadvantage of social sensitivity: interaction of oxytocin receptor genotype and child maltreatment on brain structure. Biol Psychiatry. 2016;80(5):398–405.
Cicchetti D, Rogosch F, Hecht K, Crick N, Hetzel S. Moderation of maltreatment effects on childhood borderline personality symptoms by gender and oxytocin receptor and FK506 binding protein 5 genes. Dev Psychopathol. 2014;26(3):831–49.
Krueger F, Parasuraman R, Iyengar V, Thornburg M, Weel J, Lin M, et al. Oxytocin receptor genetic variation promotes human trust behavior. Front Hum Neurosci. 2012;6:4.
Huetter FK, Bachmann HS, Reinders A, Siffert D, Stelmach P, Knop D, et al. Association of a common oxytocin receptor gene polymorphism with selfreported ‘empathic concern’ in a large population of healthy volunteers. PLoS One. 2016;11:e0160059.
Bradley B, Western D, Mercer KB, Binder EB, Jovanovic T, Crain D, et al. Association between childhood maltreatment and adult emotional dysregulation in a low-income, urban, African American sample: moderation by oxytocin receptor gene. Dev Psychopathol. 2011;23:439–52.
Hostinar CE, Cicchetti D, Rogosch FA. Oxytocin receptor gene polymorphism, perceived social support, and psychological symptoms in maltreated adolescents. Dev Psychopathol. 2014;26:465–77.
Sepehrmanesh Z, Ahmadvand A, Akkasheh G, Saei R. Prevalence of psychiatric disorders and related factors in male prisoners. Iran Red Crescent Med J. 2014;16(1):e15205.
Falcus C, Johnson D. The violent accounts of men diagnosed with comorbid antisocial and borderline personality disorders. Int J Offender Ther Comp Criminol. 2018;62(9):2817–30.
Filov I. Antisocial personality traits as a risk factor of violence between individuals with mental disorders. Open Access Maced J Med Sci. 2019;7(4):657–62.
González RA, Igoumenou A, Kallis C, Coid JW. Borderline personality disorder and violence in the UK population: categorical and dimensional trait assessment. BMC Psychiatry. 2016;16:180.
Freestone M, Howard R, Coid JW, Ullrich S. Adult antisocial syndrome co-morbid with borderline personality disorder is associated with severe conduct disorder, substance dependence and violent antisociality. Personal Ment Health. 2013;7(1):11–21.
van den Brink C, Harte J, Denzel A. Men and women with borderline personality disorder resident in Dutch special psychiatric units in prisons: a descriptive and comparative study. Crim Behav Ment Health. 2018;28(4):324–34.
de Vries B, van Busschbach J, van der Stouwe E, Aleman A, van Dijk J, Lysaker P, et al. Prevalence rate and risk factors of victimization in adult patients with a psychotic disorder: a systematic review and meta-analysis. Schizophr Bull. 2018;45(1):114–26.
Edimansyah A, Kenneth GK, Mythily S, MengEn G, Timothy L, Colin T. Validity of the personality diagnostic questionnaire-4 (pdq-4+). J Pers Disord. 2012;25(6):834–41.
Zhou L, Xiao SY, He XY, Li J, Liu HM. Reliability and validity of Chinese version of Barratt impulsiveness Scale-11. J Clin Psychol. 2006;14(4):343–4.
Post R, Leverich G, McElroy S, Kupka R, Suppes T, Altshuler L, et al. Prevalence of axis II comorbidities in bipolar disorder: relationship to mood state. Bipolar Disord. 2018;20(4):303–12.
Taylor J, James LM, Bobadilla L, Reeves MD. Screening for disinhibited disorder cases in a college population: performance of the SMAST, DAST, SCID-II-Q, and PDQ-4. Psychol Assess. 2008;20(4):351–60.
Bernstein D, Fink L, Handelsman L, Foote J, Lovejoy M, Wenzel K, et al. Initial reliability and validity of a new retrospective measure of child abuse and neglect. Am J Psychiatry. 1994;151(8):1132–6.
Zhao XF, Zhang YL, Li LF. Evaluation on reliability and validity of Chinese version of childhood trauma questionnaire. Chin J Clin Rehabil. 2005;16:209–11.
Ibrahim J, Cosgrave N, Woolgar M. Childhood maltreatment and its link to borderline personality disorder features in children: a systematic review approach. Clin Child Psychol Psychiatry. 2018;23(1):57–76.
Ferrer M, Andión Ó, Calvo N, Ramos-Quiroga J, Prat M, Corrales M, Casas M. Differences in the association between childhood trauma history and borderline personality disorder or attention deficit/hyperactivity disorder diagnoses in adulthood. Eur Arch Psychiatry Clin Neurosci. 2017;267(6):541–9.
Parker G, McCraw S, Bayes A. Borderline personality disorder: does its clinical features show specificity to differing developmental risk factors? Australas Psychiatry. 2018;26(4):410–3.
Camelia EH, Dante C, Fred AR. Oxytocin receptor gene polymorphism, perceived social support, and psychological symptoms in maltreated adolescents. Dev Psychopathol. 2014;26(2):465–77.
Bradley B, Westen D, Mercer KB, Binder EB, Jovanovic T, Crain D, et al. Association between childhood maltreatment and adult emotional dysregulation in a low-income, urban, African American sample: moderation by oxytocin receptor gene. Dev Psychopathol. 2011;23(2):439–52.
Mcquaid RJ, Mcinnis OA, Stead JD, Kimberly M, Hymie A. A paradoxical association of an oxytocin receptor gene polymorphism: early-life adversity and vulnerability to depression. Front Neurosci. 2013;7(128):A128.
Li J, Zhao Y, Li R, Broster LS, Zhou C, Yang S. Association of oxytocin receptor gene (OXTR) rs53576 polymorphism with sociality: a meta-analysis. PLoS One. 2015;10:e0131820.
Flasbeck V, Moser D, Kumsta R, Brüne M. The OXTR single-nucleotide polymorphism rs53576 moderates the impact of childhood maltreatment on empathy for social pain in female participants: evidence for differential susceptibility. Front Psychiatry. 2018;9:A359.
Cecil CA, Smith RG, Walton E, Mill J, Mccrory EJ, Viding E. Epigenetic signatures of childhood abuse and neglect: implications for psychiatric vulnerability. J Psychiatr Res. 2016;83(Complete):184–94.
Kurdziel G, Kors S, Macfie J. Effect of maternal borderline personality disorder on adolescents' experience of maltreatment and adolescent borderline features. Personal Disord. 2018;9(4):385–9.
Ludmer J, Gonzalez A, Kennedy J, Masellis M, Meinz P, Atkinson L. Association between maternal childhood maltreatment and mother-infant attachment disorganization: moderation by maternal oxytocin receptor gene and cortisol secretion. Horm Behav. 2018;102:23–33.
Belsky J, Jonassaint C, Pluess M, Stanton M, Brummett B, Williams R. Vulnerability genes or plasticity genes ? Mol Psychiatry. 2009;14(8):746–54.
The authors would like to thank the subjects who participated in the study. We also thank the prison authorities and personnel for their assistant.
This work was supported by National Natural Science Foundation of China (81901390) and Postgraduate Research & Practice Innovation Program of Jiangsu Province (SJCX19–0332). The funding body had no involvement with study design; data collection, analysis, or interpretation; or preparing the manuscript.
Ethics approval and consent to participate
The study was approved by the Institutional Review Board of Nanjing Brain Hospital, affiliated to Nanjing Medical University (Number of ethical approval: 2019-KY078–01). We carried out the study following the principle of the Helsinki Declaration. The subjects signed the informed consent form after totally understanding the purpose of our study.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Zhang, M., Liu, N., Chen, H. et al. Oxytocin receptor gene, childhood maltreatment and borderline personality disorder features among male inmates in China. BMC Psychiatry 20, 332 (2020). https://doi.org/10.1186/s12888-020-02710-0