|Home | About | Journals | Submit | Contact Us | Français|
Borderline personality disorder (BPD) is a severe disorder with high morbidity and mortality, but unknown etiology. Childhood abuse has been proposed as an etiological factor, but the mechanism by which an abuse history could influence risk for BPD has not been determined. The aim of this study was to determine whether the Tryptophan Hydroxylase 1 gene is related to BPD in a clinical sample, and whether TPH1 genotypes or haplotypes moderate the relationship between abuse history and BPD.
Three hundred ninety-eight patients diagnosed with mood disorders were genotyped for TPH1 G-6526A promoter polymorphism (rs4537731) and the A218C intron 7 polymorphism (rs1800532) and a set of ancestry informative markers, assessed for DSM IV diagnoses, and assessed for a history of physical and/or sexual abuse.
Patients with a diagnosis of BPD were more likely to be risk allele carriers (A alleles at both loci) than the non-BPD group. Logistic regression analysis predicting BPD diagnosis with both single SNPs and haplotypes showed significant interaction effects between genotype and abuse history. Poisson regression predicting the number of BPD diagnostic criteria met with the same predictor set also included a significant interaction term. Risk allele carriers with a history of abuse had an increased likelihood of a BPD diagnosis.
Variation in TPH1may increase risk for developing BPD as a result of childhood abuse. Elements of BPD pathology may be due in part to a genetically influenced serotonergic dysfunction, which in turn may lead to a differential response to environmental stressors.
Borderline Personality Disorder (BPD) is a serious, disabling condition with prevalence rates of approximately 0.7 – 1.8% of the population overall (Samuels et al. 2002; Torgersen et al. 2001). It is the most commonly diagnosed personality disorder in clinical samples, affecting an estimated 15 % of all patients in outpatient treatment (Widiger and Weissman 1991). It is also a condition with significant morbidity and mortality, with up to 80% of patients engaging in some form of self-injurious behavior (Zisook et al. 1994), and a completed suicide rate of approximately 10% (McGlashan 1986). Although little is known about the etiology of BPD, a substantial literature exists on the putative etiological correlates of the disorder. These can be broadly grouped into three main categories: family history of psychopathology, personality traits, and environmental stressors (Trull et al. 2000). Most research attention has been focused on the role of childhood trauma.
Since Herman and van der Kolk (Herman and van der Kolk 1987) first hypothesized that trauma might be an important etiological factor in BPD, many researchers have reported an association between childhood trauma and BPD, and several meta-analyses (Fossati et al. 1999) and literature reviews (Goodman and Yehuda 2002) have concluded that traumatic experiences in childhood are an important etiological factor in at least a subset of BPD cases. However, estimates of the size of this subset vary considerably. For example, in their 1999 meta-analysis, Fossati (Fossati et al. 1999) and colleagues reported only a modest overall effect size (r=.28), and concluded that their findings did not seem to support the notion that childhood sexual abuse was a major psychological risk factor or antecedent of BPD. In a later review, Goodman and Yehuda (Goodman and Yehuda 2002) concluded that BPD patients report a history of childhood sexual abuse with much greater frequency than non-BPD patients, and estimate that approximately 40 – 70 % of BPD patients have a history of childhood trauma. However, they also point out that these types of experiences are not uniquely associated with BPD.
The fact that a large percentage of individuals who report a history of childhood abuse will not develop BPD in later life (Collishaw et al. 2007) suggests that there are important moderators of the relationship between abuse and BPD that remain unelucidated. Recently, Goodman, New and Siever (Goodman et al. 2004) proposed an etiological model of BPD in which inherited susceptibility interacts with environmental stressors such as abuse and trauma to increase the risk for the development of BPD. In this study, we sought to test this model using single nucleotide polymorphisms in the Tryptophan Hydroxylase I gene (TPH1). Specifically, we tested the hypotheses that 1) previously identified TPH1 risk alleles will be associated with BPD, and that 2) a history of childhood abuse will be associated with a BPD diagnosis, but only in the presence of the TPH1 risk alleles.
We chose to examine TPH1 polymorphisms for several reasons. Tryptophan hydroxylase is the rate-limiting enzyme in the serotonin metabolic pathway, and is responsible for catalyzing the conversion of tryptophan to 5-hydroxytryptophan in the periphery, and to a lesser extent in the CNS (TPH2 is the predominant isoform in the CNS). Because of the importance of serotonin in the stress response (Holmes 2008), investigation of genes related to the serotonin system may provide important information concerning susceptibility to the effects of stress-inducing life events such as childhood trauma. Also, we have recently shown that the TPH1 A218C polymorphism could distinguish BPD patients from healthy controls (Wilson et al. 2009), suggesting that mutations in this gene may play a role in the etiology of BPD. In addition to its association with BPD, TPH1 polymorphisms have been associated with several traits that are characteristic of patients with the disorder, including hostility (Hennig et al. 2005; New et al. 1998; Reuter and Hennig 2005), anger (Baud et al. 2009), harm avoidance (Anghelescu et al. 2005), and impulsiveness (Staner et al. 2002). Finally, associations between TPH1 polymorphisms and suicidal (Bellivier et al. 2004; Liu et al. 2006a; Paik et al. 2000) and non-suicidal self-injurious behavior (Evans et al. 1327; Pooley et al. 2003) have been reported by several investigators. The fact that self-injurious behavior occurs with high frequency in patients with BPD (Simeon et al. 1992) also suggests an association with the disorder. While several studies have published results that conflict with some of these findings (Bennett et al. 2000; Kunugi et al. 1999; Zalsman et al. 2001), very few studies of these conflicting studies have accounted for BPD rates in their samples. As our initial findings with TPH1 and BPD suggest that the inconsistency could be due to varying rates of BPD in these samples, we also aimed to test this hypothesis with the current study. While both the G-6526A and the A218C polymorphisms have been extensively studied, neither has been examined in a sample including suicide attempters and nonattempters, with and without BPD. With the inclusiveness of our sample, our aim was to clarify some of the conflicting findings on TPH1, in the hope of better identifying the phenotype most strongly associated with this gene.
A total of 398 unrelated research participants were included in this study, selected from a larger sample admitted to a multisite project on mood disorders and suicidal behavior. All participants were diagnosed with an Axis I mood disorder, either major depressive disorder (n=312) or a bipolar disorder (n=86). Ninety-eight participants were diagnosed with borderline personality disorder, while the remaining 300 either had no Axis II diagnosis (76%, n=226) or one or more non-BPD diagnoses (24%, n=74). Of the 24% with another Axis II disorder, 3.3% (n=10) had a Cluster A diagnosis, 6% (n=18) had a cluster B diagnosis, and 15% (n=45) had a Cluster C diagnosis. The BPD sample has been previously studied and compared to a sample of healthy controls (Wilson et al. 2009). All participants were of Hispanic (20.9%) or non-Hispanic Caucasian (79.1%) descent. We limited the racial and ethnic composition of the sample in an effort to reduce the possibility of sample stratification. The decision to include Hispanic Caucasians was made based on recent research showing relatively minimal genetic distance between self-identified Hispanics and individuals of Caucasian descent (Tang et al. 2005; Wilson et al. 2009). As Hispanic individuals generally have a background that includes varying proportions of European, Native American, and/or African ancestry (Bertoni et al. 2003) and can be of any race (American Community Survey, US Census Bureau, http://www.census.gov/prod/2007pubs/acs-03.pdf), eliminating individuals with self-reported Black Hispanic heritage as we did in this study likely further minimized this distance. Exclusion criteria were the presence of persisting psychotic symptoms or psychotic disorder diagnosis, history of severe head trauma, organic mental syndromes, medical illnesses that could affect the brain or increase intracranial pressure, or the presence of mental retardation or cognitive impairment that might interfere with the completion of the assessments or the obtaining of informed consent. The study was approved by all applicable institutional review boards, and all participants provided written informed consent.
Diagnoses were determined using the Structured Clinical Interview for DSM-III-R/DSM-IV, Patient Edition (Spitzer et al. 1990) and the Structured Clinical Interview for DSM-IV Axis II Personality Disorders (First et al. 1995). The presence or absence of physical or sexual childhood abuse before the age of 15 was determined by responses to three direct questions from the Columbia Demographic and Treatment History Interview (Brodsky et al., 2001): “Any history of physical and/or sexual abuse over lifetime” (yes/no); “If yes, describe and code: physical/sexual/both/not applicable”; “If yes, did abuse take place before age 15 years?” (yes/no/not applicable). Physical abuse was defined as any physical disciplinary measure inflicted by a parental figure or caregiver that was severe enough to leave marks such as welts or bruises. Sexual abuse was defined as any sexually inappropriate contact occurring between the participant and an adult, or between the participant and another minor at least five years older.
Participants were grouped according to whether they reported no history of abuse, a history of either physical or sexual abuse, or a history of both physical and sexual abuse. This categorical grouping was recently used by Binder and colleagues (Binder et al. 2008) in a study investigating gene × environment interactions and the risk of post-traumatic stress disorder with good results. Because of cell size limitations of some of the analyses, a dichotomous Yes/No abuse variable was used in all analyses. Models with significant interaction terms were rerun with the categorical variable if appropriate, in order to determine whether severity of abuse determined the strength of the effect. All interviews were conducted by at least a master’s level clinicians with specific training in the use of the SCID instruments and clinical measures.
All participants were genotyped for the G-6526A promoter polymorphism (rs4537731) and the A218C intron 7 polymorphism (rs1800532) of the TPH1 gene. A218C genotypes were available for the full sample, while G-6526A genotypes were available for 300 out of 398 cases (75/98 BPD, 225/300 without BPD). Genomic DNA was extracted from lymphocytes as described by Huang et al. (2003) and from epithelial cells from cheek swabs according to the manufacturer’s procedure, and amplified via polymerase chain reaction (PCR). PCR amplification was performed in 20 µl of reaction mixture containing 1× PCR buffer, 40–100 ng DNA, 2 mM MgCl2, 4% of DMSO (dimethyl sulfoxide), 50 nM of each dNTP, 0.8 U RedTaq polymerase (Sigma, St Louis, MO), and 40 ng of each primer. Forward primer (TPH5F) 5’-TGGCATTGAAGTAAGAGCAC-3’ and reverse primer (TPH5R) 5’-GTTTCATGCAGGTATTAGTG-3’) were used for rs4537731. Forward primer (Mann-45), 5’-GCCAGGAATTCATCAATGG-3', and reverse primer (Mann-46), 5’-CCACCACATACACACCCAAA-3’ were used for rs1800532. The reaction conditions were as follows: 95°C for 4 min, 30 cycles at 95°C for 40 s, 40 s at 54°C, 40s at 72°C and a final extension step at 72°C for 3 min. PCR amplifications were performed in a DNA Robocycler (Stratagene, La Jolla, CA). PCR fragments for rs4537731 were digested with the Sau3AI (or MboI) restriction enzyme, while fragments for rs1800532 were digested with the NheI restriction enzyme. Electrophoresis was performed in 1.4% agarose and the PCR amplification products visualized with UV light to determine genotype.
In order to control for any hidden population stratification, we genotyped a set of 186 unlinked markers using a custom designed Illumina GoldenGate 96-well format Sentrix® arrays (Hodgkinson et al. 2008). A total of 500ng of sample DNA was used per assay. All pre-PCR processing was performed using a TECAN liquid handling robot running Illumina protocols. Arrays were imaged using an Illumina Beadstation GX500 and the data analyzed using GenCall v184.108.40.206 and GTS Reports software v220.127.116.11 (Illumina). Ethnicity factor scores for each individual were estimated by STRUCTURE using 1,051 CEPH subjects as reference. Proportion of membership of each subject in each of 7 clusters was estimated. These scores correspond to 7 geographic regions and sum to 1. Factors were labeled by geographic region of ancestry ("Africa", "Europe", "Middle East", "Central Asia", "Far East Asia", "Oceania", and "America"). To adjust for ethnic stratification in the SNP-wise analyses with a minimum degrees of freedom lost, it was necessary to determine which of these factors differed by TPH1 genotypes. Logistic regression was performed for each individual SNP, with the presence of the minor allele as a response variable and “Europe”, “Africa”, “America”, “Middle East”, and “Central Asia” as predictor variables. “Far East Asia” and “Oceania” factors had very low mean scores overall in the sample (means< 0.01) and were excluded from all analyses a priori.
As the TPH1 gene is covered by one large haplotype block, we created haplotypes based on the two polymorphisms. Haplotypes were imputed using the SIMHAP package (Carter et al. 2008) for the R statistical software program (http://cran.r-project.org), using genotyping information from the subset of 300 subjects with both genotypes. Haplotypes were imputed using genotyping information from the entire subsample. For individuals heterozygous at either of the two loci, chromosomal assignments were determined using the haplotype with the highest estimated posterior probability. These data were used to create codes indicating the number of haplotypes for each individual (0, 1, or 2), which were then used as predictors in our statistical models. Statistical tests of haplotype main effects and interactions were conducted using the −6526G and 218C alleles as the base haplotype. This two-SNP haplotype is the most common in the Hapmap CEPH sample, which is a population from Utah with primarily northern and western European ancestry (Thorisson et al. 2005), with a frequency of .442. Linkage disequilibrium coefficients were determined using the JLIN Java program(Carter et al. 2006).
Means and standard deviations were computed for all variables in the analyses, and all variables were checked for normality and outliers. In order to test for group differences on the demographic variables, chi-squares were computed for gender and race (Hispanic versus non-Hispanic), and a t-test for independent samples was computed for age. Chi-squares were computed for genotype, as well as for a dominant risk allele models (AA/GA versus GG, AA/AC versus CC). Although the G allele is the minor allele for the G-6526A polymorphism, previous research has identified the A wild-type allele as the risk allele at this locus (Liu et al. 2006a; Rotondo et al. 1999).
We used the Generalized Linear Model to run several types of analyses. Analyses were run with both individual genotypes and each of the two SNP haplotypes in the subset. We used a logistic binomial model to analyze the relationship between abuse history, genotype, and BPD diagnosis, controlling for age, gender, and attempter status, testing both additive and dominant genetic models. We also used a Poisson model to analyze the relationship between abuse history, genotype, and number of DSM-IV BPD criteria met (BPD criteria was available for a subset of the full sample, 373 of 398) with similar covariates and model assumptions. We used a dichotomous abuse variable (presence/absence) in all primary GLM models. Models with significant findings were reanalyzed using a categorical abuse history variable (none, physical or sexual, both physical and sexual) to assess whether severity of abuse history was a significant determinant of outcome. In order to adjust for multiple testing, we conducted a permutation-based analysis to empirically determine adjusted p-values for all covariates and interaction terms in the GLM models with significant findings based on the method of Wagner et al (2008). We conducted 1,000 permutations for each model, and present both original and adjusted p-values for the interaction terms. Differences in the levels of self-reported impulsiveness and hostility were analyzed with univariate ANOVA, comparing both the risk allele groupings and genotypes in the full sample as well as in the individual comparison groups (BPD versus no-BPD).
Demographic information comparing the two groups is presented in Table 1. The BPD group was younger than the control group and had more females. There was no difference in racial/ethnic composition of the groups (BPD 23.5% vs. No BPD 22.6% Hispanic Caucasian, χ2 = <1.0, df = 1, p = ns). Genotype frequencies for both SNPs were comparable in males and females in the full sample (G-6526A χ2 = 1.26, df = 2, ns, A218C χ2 = 1.25, df = 2, ns), the BPD group (G-6526A χ2 < 1.00, df = 2, ns, A218C χ2 < 1.00, df = 2, ns), and the non-BPD group (G-6526A χ2 = 1.05, df = 2, ns, A218C χ2 = 3.47, df = 2, p = .18). Consequently, males and females were combined for all analyses. There was no difference in Axis I mood disorder diagnoses between the groups, although there was more than twice the percentage of suicide attempters in the BPD group (72.4% vs. 31.0%, χ2 = 52.38, df = 1, p = .000). Because of this substantial disparity, in addition to including attempter status as a covariate in our analyses, we also confirmed our findings with statistical models including only suicide attempters.
Genotype and haplotype frequencies and statistics are presented in Tables 2, ,3,3, and and4.4. Genotypes were in Hardy Weinberg equilibrium in the full sample (G-6526A χ2 =< 1.0, df=2, ns; A218C χ2 =< 1.0, df=2, ns) as well as in BPD (G-6526A χ2 < 1.0, df=2, ns; A218C χ2 =2.65, df=2, p=.27) and non-BPD (G-6526A χ2 < 1.0, df=2, ns; A218C χ2 < 1.0, df=2, ns) cases independently. As an additional data check, we also examined HWE in males and females, as well as in Hispanics and non-Hispanics in the full sample and found no deviations. The SNPs were in strong linkage disequilibrium in our sample (D’=.79). There were no significant group differences in either the number of −6526A allele carriers (AA/GA genotypes, χ2<1.0, df=1, ns) or the G-6526A genotype frequencies (χ2<1.0, df=2, ns). The BPD group had a higher frequency of 218A allele carriers (AA/AC genotypes) than the control group (χ2=5.61, df=1, p=.02), as well as a significant difference in A218C genotype frequencies (χ2=5.86, df=2, p=.05). Although neither SNP was related to suicidal behavior in the BPD group, the A218C polymorphism was associated with suicide in the N-BPD group (χ2=9.06, df=2, p=.01). Interestingly, for both SNPs, only genotype frequencies in the N-BPD nonattempter group were similar to data published from the HapMap project on the CEPH population. All other subgroups deviated to some degree toward increased percentages of risk allele carriers.
Our first set of analyses modeled the relationship between BPD diagnosis (Y/N), abuse history (Y/N), and genotype for each of the SNPs independently, controlling for attempter status, gender, and age, under both additive and dominant risk allele models. In the first subset, predicting BPD diagnosis by G-6526A, the dominant risk allele model (AA/GA versus GG) provided the best fit to the data. The overall model was significant (χ2=65.50, df=6, p=.000), as was the interaction term (Wald=4.46, df=1, p=.035/adjp=.036). Neither the main effect for abuse history nor the main effect for genotype was significant once the interaction term was included in the model. Results were similar when the model included the categorical abuse variable (genotype × abuse interaction term Wald=4.45, df=1, p=.023/adjp=.025). In the second subset, predicting BPD diagnosis by A218C, the dominant risk allele model (AA/AC versus CC) also provided the best fit to the data. The overall model was significant (χ2=100.32, df=7, p=.000), as was the interaction term (Wald = 8.06, df=1, p = .005/adjp=.003). The main effect for abuse history was significant (Wald = 6.06, df=1, p = .014/adjp=.011), although the main effect for genotype was not. As with the −6256 analysis, results were similar when including the categorical abuse variable (genotype × abuse interaction term Wald=5.70, df=1, p=.017/adjp=.016). Risk of a BPD diagnosis for carriers of the A allele at both loci was increased in the presence of an abuse history. In addition, an increase in severity of abuse was associated with increased risk of a BPD diagnosis (see Figures 1 and and2).2). Neither the main effect for genotype nor the interaction term was significant in any of the additive risk models. Of note was the fact that in all of these analyses, the main effect for abuse history was not significant when the interaction term was included in the model.
The second set of analyses modeled the relationship between SCID-II BPD criteria, abuse history, and genotype for each of the SNPs independently as well as for the two-SNP haplotype, controlling for attempter status, gender, and age, again under both additive and dominant risk allele models. For the G-6526A polymorphism, none of the models resulted in a significant genetic effects or gene × abuse interaction terms. For the A218C polymorphism, the overall test of the dominant model was significant (χ2=68.95, df=7, p=.000), as was the interaction term (Wald=4.56, df=1, p=.033). The interaction between genotype and abuse history predicted the number of DSM-IV BPD criteria met across the entire sample, after controlling for age, sex, and attempter status. There were no significant effects in the model run with the categorical abuse variable.
Haplotypes were imputed for a subset of cases with genotypes for both SNPS (n=300). The AA risk haplotype was used in all analyses. As with the single SNP analyses, the dominant genetic risk model provided the best fit to the data. In the analysis predicting BPD diagnosis by the −6256/218 haplotypes the overall model was significant (χ2=71.87, df=7, p=.000), as was the interaction term (Wald=5.74, df=1, p=.017/adjp=.024). Individuals with one or two AA haplotypes were at increased risk for a BPD diagnosis when a history of abuse was present. As with the single SNP analyses, none of the genetic or interaction terms were significant under the additive model, and neither of the Poisson models included significant genetic or genotype × abuse interaction terms. As with the individual SNP models, neither the main effect for abuse history nor haplotype were significant once the interaction term was included in the model.
As several investigators (Liu et al. 2006a; Roy et al. 2001; Tsai et al. 1999), including members of our group (Galfalvy et al. 2009), have reported an association between TPH1 polymorphisms and suicidal behavior, in addition to controlling for suicide attempter status in our analyses, we also reran the logistic haplotype analysis in the subset of cases and controls with a history of suicidal behavior (n=164) under both additive and dominant risk allele models. Results confirmed the associations found in the original analyses, with the dominant model providing the best fit to the data, a significant main effect for the AA haplotype (Wald=4.96, df=1, p=.026) and a significant interaction term (Wald=9.17, df=1, p=.002).
Our results suggest that TPH1 moderates the association between a history of childhood abuse and the risk for developing BPD later in life, and that increasing severity of abuse may be associated with increasing BPD risk in those carrying the TPH1 risk alleles. In addition, the interaction between abuse history and individual genotypes, as well as the two SNP haplotype, not only increased risk for a BPD diagnosis, but also seemed to influence severity of psychopathology, as indicated by the number of DSM-IV diagnostic criteria met for BPD. Notably, these associations were found both in the full sample, as well as in a restricted sample comprised of only suicide attempters. This finding suggests that the observed interaction between TPH1 genotype and abuse history were not due to the association between TPH1 and suicidal behavior observed in the non-BPD group. Also, as there was no significant difference in Axis I mood disorder diagnoses between the BPD and Non-BPD groups, it is unlikely that our results reflect an association between TPH1 and Axis I pathology.
Although there is limited research investigating the interaction between TPH1 and a history of abuse, an interesting study by Keltikangas-Jarvinen and colleagues (Keltikangas-Jarvinen et al. 2007) reported results that lend support to our findings in a study exploring whether TPH1 polymorphisms interact with hostile early childhood environment to increase the personality trait of harm avoidance in healthy controls. Using a measure that assessed for the presence of emotional rejection, inconsistent discipline, and maternal neglect, the authors found that the A allele of the A218C interacted with a hostile childhood rearing environment to increase harm avoidance in otherwise psychiatrically healthy women. As high levels of harm avoidance have been reported in personality studies of BPD patients (Korner et al. 2007), as well as in the offspring of mothers with BPD (Barnow et al. 2006), the results of this research lend support to our findings.
Similar to our previous study comparing this sample to healthy controls, we also found that BPD patients had higher frequencies of both TPH1 polymorphisms and the two SNP haplotype in comparison to depressed patients without BPD. While this is not a true replication of the initial finding, since we are comparing the same BPD sample to a second control sample of patients with mood disorders but without BPD, the fact that comparison to a clinical comparison sample also yielded significant findings further supports the view that TPH1 is associated with BPD. In addition, we found that while neither of the polymorphisms distinguished suicide attempters from non-attempters in the BPD group, the intron 7 polymorphism was significantly associated with suicidal behavior in the non-BPD group. Although these associations appear to be somewhat conflicting, they provide support for our view that TPH1 may be associated with an endophenotype that is associated with BPD, but that can also be more frequently found in suicidal individuals without BPD. Indeed, visual inspection of the genotype frequencies for the polymorphisms, and for the G-6526A polymorphism in particular, strongly suggests that this is the case. Genotype frequencies for the BPD groups and the non-BPD attempters are similar, whereas the non-BPD nonattempters have frequencies that are similar to healthy controls.
Our findings may help to integrate some of the inconsistent findings in the literature regarding the association between TPH1 and suicidal behavior. As with single SNP association studies, several investigators have reported associations between TPH1 haplotypes and suicidal behavior in several populations. For example, Lui and colleges (Liu et al. 2006b) reported that a haplotype containing A alleles of both the G-6526A and the A218C SNPs was most strongly related to suicide in a Chinese psychiatric sample of predominantly patients with schizophrenia. Also, in a recent study by our group (Galfalvy et al. 2009) examining the relationship between TPH1 and suicide in 343 patients with mood disorders, we found the A/A haplotype was associated with history of suicidal behavior, and significantly predicted suicide risk prospectively. Although the research including the G-6526A promoter polymorphism is limited, many studies including the A218C polymorphism have also found no association between TPH1 and suicidal behavior (Kunugi et al. 1999; Souery et al. 2001) or suicidal ideation in depressed individuals(Du et al. 2001). In addition, a recent meta-analysis of studies investigating the association between TPH1 and suicide (Saetre et al. 2009) reported that results from association studies were highly heterogeneous, and concluded that TPH1 polymorphisms seem to be associated with increased risk for disorders that have high rates of suicide, rather than specifically and directly increasing risk for suicidal behavior. This heterogeneity could be explained in part by varying rates of BPD in the samples included in these studies.
Although little is known about the effects of mutations in the TPH1 gene, the G-6526A promoter polymorphism has been shown to be functional. The polymorphism resides within a c-ETS transcription factor binding domain (Rotondo et al. 1999), and the A allele has been associated with an approximately 3 fold increase in gene transcription (Sun et al. 2005). As sequence analysis shows that the presence of the mutant allele activates the binding site, the finding that the A allele is associated with increased TPH1 expression suggests that the mutation leads to increased transcription of TPH1. In contrast to the G-6526A polymorphism, no functional studies have been undertaken with the A218C polymorphism. Because it is located in an intronic region, it does not lead to any functional change in protein coding. However, the polymorphism is also in a potential transcription factor (GATA-1) binding site. The GATA-1 transcription factor has been shown to increase transcription rates of some genes up to 100 fold (Wong et al. 2004), but also can act as a repressor (Raich et al. 1995). In either case, a binding site created via the transmission of the A allele at this locus could potentially have a significant effect on TPH gene expression.
The mechanism by which differential expression levels of the TPH1 gene could exert an effect on CNS serotonergic function is not well understood. Tryptophan hydroxylase is a critical precursor enzyme in the serotonin metabolic pathway, and TPH1 mutations, particularly if they led to changes in gene expression or enzyme functionality, could influence serotonin levels. However, as research has shown that TPH2 is the predominant isoform in the adult CNS, while TPH1 is primarily active in the periphery (Sakowski et al. 2006), it is difficult to establish a direct link between TPH1 activity and CNS serotonin. Anecdotal evidence from animal studies suggests several possible indirect connections. First, animal studies have shown that chronic restraint stress can lead to significant increases in CNS TPH1 expression in adulthood (Abumaria et al. 2008; Abumaria et al. 2007), which suggests that differences in TPH1 functionality may be more evident in the context of stress induced gene transcription. Alternatively, animal studies have shown that TPH1 is expressed much more prominently in the brain during the late prenatal developmental period, and that TPH1 may play a role in the development of the serotonin system via an influence on neuron number and growth (Nakamura et al. 2006; Rind et al. 2000). Differences in TPH1 expression or functionality could have an effect on prenatal serotonin levels, which could in turn lead to subtle deficiencies in serotonergic functioning in later life. Research on the effects of antidepressant medication use during the prenatal period has demonstrated the deleterious effects of elevated serotonin exposure in utero, and has been associated with several adverse outcomes, including alterations in HPA reactivity (Oberlander et al. 2008).
The significant interaction terms in our analyses raise the question of how the putative genetic liability created by mutations in the TPH1 gene would impact the experience of, and the neurobiological consequences of childhood trauma. Many investigators have shown that childhood trauma is associated with significant neurobiological sequelae (for review, see Teicher et al. (Teicher et al. 2002). While the mechanism underlying these sequelae are not entirely clear, recent research has shown an association between post-trauma cortisol levels and hippocampal volume reductions over a 12–18 month period in children with PTSD (Carrion et al. 2007). This suggests that the prolonged activation of the HPA axis and increased exposure to stress hormones that are consequences of ongoing abuse can have long term neurobiological consequences on brain development and function. Because of serotonin’s role in the stress response system (Chaouloff 2000; Chaouloff et al. 1999), individual differences in serotonin production could alter the capacity to mitigate this effect. Animal studies have shown that among serotonin’s many functions, its release during stress serves to attenuate the damaging effects of corticosteroids by increasing mineralocorticoid receptor density in the hippocampus (Gesing et al. 2000; Mitchell et al. 1990) via binding at the 5-HT7 receptor (Laplante et al. 2002), and thus serves to increase the resiliency through increasing the modulatory capacity of the stress response system. It is possible that any reduction in the capacity for serotonin synthesis under basal and/or stressful conditions could lead to a reduced capacity to offset the damaging effects of stress, and thus to more severe neurobiological consequences. Again, more research is needed if we are to better understand what is likely a complex interplay between genetic liabilities and environmental factors.
This study has several limitations. First, the modest size of the sample included in this research means our findings must be replicated and validated. While association studies are often difficult to replicate, the fact that our findings are consistent with the large literature on TPH1 polymorphism associations increases the likelihood that our findings are valid. Also, the fact that we studied only two polymorphisms means that we did not capture all the known variation in the TPH1 gene. Several additional tagging SNPs are required in order to achieve this. Finally, TPH1 allele frequencies are known to vary considerably between ethnic groups, and even between northern and southern Europe (Nielsen et al. 1997). While we made an effort to control for hidden stratification using structured association methods, it is possible that there is additional unknown stratification in our sample that may have contributed to the results.
Despite these limitations, our results suggest that the TPH1 is associated with an underlying risk for both suicidal behavior and borderline personality disorder, and that TPH1 moderates the relationship between childhood trauma and risk for developing BPD later in life. Additional research on the neurobiological effects of genetic variation in the TPH1 gene will be important in understanding the manner in which TPH1 influences risk for both suicidal behavior and BPD, as well as in determining more definitively whether TPH1 mutations create a genetic vulnerability that leads to greater neurobiological sequelae as a result of childhood abuse or trauma, and a consequent increase in risk for BPD. Research focused on determining the effects of variation at the A218C locus is also of particular importance.
This work was funded in part by grants from the National Institute of Mental Health (MH 48514, MH62185, MH59710, MH61017, and MH62665) and the National Institute on Alcohol Abuse and Alcoholism (AA15630).
Publisher's Disclaimer: This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting proof before it is published in its final citable form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain.
The authors of this manuscript have no financial or other relationships that are relevant to the subject matter of this research or the results described herein.