Context
The tryptophan hydroxylase 2 (TPH2) gene encodes the first (also the rate-limiting) enzyme in the serotonin biosynthetic pathway. Despite reports of possible associations between polymorphisms in human TPH2 and many psychiatric disorders, including bipolar disorder (BPD), the functional effect and susceptibility loci of such polymorphisms for BPD have not yet been identified.
Objectives
To examine the association of TPH2 with BPD and to identify the functional variants that may be involved in the pathophysiological development of BPD.
Design, Setting, and Patients
We systematically screened all exons and promoters of the TPH2 gene in Han Chinese subjects to identify sequence variants. Association tests were conducted in 105 cases and 106 control subjects using single-locus, linkage disequilibrium, and haplotype analyses. Two promoter and one exon 2 single-nucleotide polymorphisms were examined for their functional role using a reporter gene system and enzyme activity assay, respectively. Additional statistical analysis was performed to study the interaction between the 2 TPH genes in 205 study participants with TPH1 and TPH2 genotype data.
Results
Significant haplotype association of TPH2 polymorphisms and BPD was identified (P < .001). In addition, allelic alteration of polymorphisms in the promoter region and exon 2 of TPH2 caused noteworthy functional losses in promoter and enzyme activities, respectively, indicating the potential susceptibility loci for BPD. We found that the odds ratio changed from 3.73 of the TAG haplotype to 4.81 or 1.68, depending on the combined effect of both TPH genotypes. These data suggested an interaction between the 2 TPH genes to confer a risk for BPD.
Conclusions
This study supports the involvement of TPH2 in the etiology of BPD, and the functional single-nucleotide polymorphisms identified herein might be the susceptibility loci for BPD. Although the interaction between the 2 TPH genes merits further investigation, our findings suggest that the interactive effect should be considered in future studies of serotonin-related disorders.
Serotonin (5-hydroxytryptamine) is a neurotransmitter synthesized in the raphe nuclei of the brainstem that controls a wide range of physiological events in the central nervous system.1 Accordingly, disruption of serotonergic function has been implicated in the pathogenesis of many psychiatric disorders, including bipolar disorder (BPD).2-4 The biosynthesis of serotonin is initiated by tryptophan hydroxylase (TPH), using tryptophan as a substrate to generate 5-hydroxytryptophan, followed by decarboxylation using aromatic amino acid decarboxylase to produce serotonin.5 Because TPH is the first as well as the rate-limiting enzyme in the synthetic pathway for serotonin, the TPH gene has been considered a candidate for many behavioral and psychiatric traits.6 Two TPH isoforms encoded by different genes have been identified in humans. The first, TPH1, was discovered earlier and was associated with various neuropsychiatric disorders such as manic-depressive illness,7 depression,8,9 suicidality, and alcoholism.10,11 Despite many reports of positive associations, reports finding no association between TPH1 polymorphisms and many neuropsychiatric disorders have also been published12-14; thus, the role of TPH1 in the development of different neural disorders is still inconclusive.
The gene for the second isoform, TPH2, was identified recently because TPH1-knockout mice expressed normal amounts of serotonin in the brain but not in the periphery,5,15 which resulted in a cardiac dysfunction phenotype.16 The human TPH2 protein shares 71% of its amino acid sequence identity with TPH1; all residues that have important functions for the structural or enzyme activity of TPH1 are also conserved.5 Using both in vitro and in vivo systems, studies confirmed that TPH2 controls brain serotonin synthesis.17 Later studies further demonstrated that the TPH2 gene is predominantly expressed in the brainstem, especially in the raphe nuclei neurons, in a rodent model18 and in postmortem brain tissue,19 whereas TPH1 is expressed in the pineal gland and peripheral tissue.5,18 The expression patterns of the 2 proteins were confirmed recently by using monospecific polyclonal antibodies against TPH1 and TPH2.20 These studies suggested a duality hypothesis for TPH1 and TPH2, namely, that they are responsible for the peripheral and major central nervous system serotonergic effects, respectively.5,18
The TPH2 polymorphisms were associated with many psychiatric disorders such as early-onset obsessive-compulsive disorder,21 suicide,22,23 attention-deficit/hyperactivity disorder,24,25 and repetitive behaviors in autism.26 In addition, disease-associated haplotypes have been identified for major depression27 and BPD.28 Furthermore, higher levels of TPH2 expression have been found in the raphe nuclei of suicidal patients with major depression29 and in the dorsolateral prefrontal cortex of patients with BPD30 than in control subjects without psychiatric disorders. Results from these studies strongly support the involvement of the TPH2 gene in the development of affective disorders. Because of the lack of potentially functional single-nucleotide polymorphisms (SNPs) in these studies, the polymorphisms that represent susceptibility loci of the TPH2 have not yet been identified.28
The discovery of TPH2 has changed ideas about linking the TPH1 polymorphism with a variety of psychiatric diseases15; however, the question of whether the 2 proteins are regulated independently—and thus do not interact—remains unanswered. Because the previous study results concerning TPH1 association with BPD were inconsistent,7,12 we propose that the TPH1 gene alone has little or no effect on the etiology of BPD, but it may interact with the TPH2 gene to influence BPD development. To test this hypothesis, the present study aimed to examine the effect of the TPH2 gene on the etiology of BPD in 105 patients with bipolar I disorder and 106 matched controls who had been previously involved in our TPH1 association study,12 and to combine both data sets for further statistical analysis. Results from our study confirm the important role of the TPH2 gene in the etiology of BPD. Moreover, we have identified common variants of the TPH2 gene that exhibit significant loss in TPH functions. Finally, our data provide statistical evidence, for the first time, that the interaction between these 2 TPH genes might confer the risk for BPD.
The patients were all Han Chinese subjects recruited for a series of studies that began in 1998 and sought to examine the association of genes involved in the serotonin system and the etiology of BPD.12,31,32 A detailed description of the participants is available in the supplementary “Methods” section on our Web site (http://140.116.6.160/SUNLab/data/TPH2/).
Sequences in the coding and up to 1 kilobase (kb) of the 5′ regulatory regions of the human TPH2 gene were screened using the direct sequencing method. Polymerase chain reaction (PCR) primers were designed to generate fragments to cover 11 exons from the main transcript (GenBank accession number, NM_173353) and 2 alternatively spliced exons from the transcript (GenBank accession number, AK094614) plus 1 kb upstream from the transcription start site. The primers used for polymorphism screening are available in supplementary eTable 1 on our Web site.
Marker selection and genotyping assays
Markers selected for further investigation were based on criteria provided in the supplementary “Methods” section available on our Web site. The TaqMan Custom SNP Genotyping Assay by Design system (Applied Biosystems, Foster City, California) was applied for genotyping in all participants in this study. Primers and probes designed for this study and genotyping data for the 6 markers investigated in all participants are given in the supplementary eTable 2 and eTable 3, respectively, available on our Web site.
Quantification of brain tph
Cell Culture and Human Brain Complementary DNA Panel
Two human neuroblastoma cell lines, IMR-32 and SH-SY5Y (American Type Culture Collection, Manassas, Virginia), were grown in appropriate media. The endogenous serotonin amounts in IMR-32 and SH-SY5Y cells were measured as 2.9 and 6.2 ng/mL, respectively. In addition, the expressions of TPH2 and Pou domain class 3 transcription factor 2 (POU3F2) in these cell lines were confirmed by quantitative reverse transcription PCR (Q–RT-PCR) and Western blot analyses (supplementary eFigure 1, available on our Web site), respectively.
Using Q–RT-PCR, 8 complementary DNA libraries made from various sections of the human brain (BioChain Institute, Hayward, California) were used to assay the relative messenger RNA (mRNA) expression level of the TPH1 and TPH2 genes.
RNA Isolation and Real-Time Q–RT-PCR
Total RNA from cultured cells was isolated using a reagent (REzol C&T; PROtech Technology, Taipei, Taiwan) according to the manufacturer's protocol. The RNA samples were used for complementary DNA synthesis. We performed Q–RT-PCR of the TPH1, TPH2, and 18S rRNA genes in a sequence detector (ABI PRISM 7900; Applied Biosystems) using the reagents (TaqMan Assay on Demand probe; Applied Biosystems) according to the manufacturer's protocol, and primer sets (HS00188220-m1, HS00998775-m1, and HS99999901-s1 for TPH1, TPH2, and 18S rRNA, respectively; Applied Biosystems). The relative levels of expression were presented as ΔCt for normalized cycle numbers or as ΔΔCt to further standardize ΔCt in each tissue with ΔCt in the brain.
Functional Assay of TPH2 5′-Promoter Polymorphisms
Details on the preparation of various reporter constructs and procedures of luciferase assay, electrophoretic mobility shift assay, and pull-down assay are provided in the supplementary “Methods” section available on our Web site.
Functional Assay of TPH2 Exon 2 C2755A Polymorphism
Detailed procedures for plasmid construction, transfection, and the enzyme activity assay are provided in the supplementary “Methods” section available on our Web site.
The significance level for all statistical tests was .05. We applied Bonferroni corrections for all multiple tests in the association study. Detailed procedures for the association study, linkage disequilibrium (LD) mapping, haplotype construction, interaction plot, and statistics for functional assays are given in the supplementary “Methods” section available on our Web site.
IDENTIFICATION OF TPH2 SEQUENCE VARIANTS
We systematically sequenced all functional regions of the human TPH2 gene in a panel of 54 samples and identified 10 single-nucleotide variants. Their positions in relation to the transcription start site of the TPH2 gene and the distribution of alternative alleles in the screening population are summarized in Table 1. The variant G33775A was observed only once in a heterozygote individual (minor allele frequency, 0.8%), and thus it could be a mutation. Of the rest, 2 of each polymorphism were in the promoter region (T − 703G and T − 473A), the 5′- and 3′-untranslated regions (A90G and G93329A), and the intron regions (C10662T and G26270A); and 3 coding SNPs contained 2 synonymous changes (G40237A and A83610T) and 1 nonsynonymous change (C2755A). To our knowledge, the SNPs C2755A and G26270A have not been reported in any public databases or literature, thus they are new and probably population-specific polymorphisms in Han Chinese.
Single-locus association analysis
Because this study aimed to identify the association between potentially functional SNPs and BPD, 6 SNPs—T − 703G (rs4570625), T − 473A (rs11178997), A90G (rs11178998), C2755A, C10662T (rs11179003), and G93329A (rs17110747)—were chosen for further genotyping. All markers were in Hardy-Weinberg equilibrium, except for SNP C2755A in patients who showed significant deviation from Hardy-Weinberg equilibrium (P < .001; available in supplementary eTable 4 on our Web site).
The distributions of alleles for each SNP and the Akaike information criterion values calculated from independent vs dependent models are given in Table 2. A significant difference in allele distribution between patients and controls was obtained for SNP C2755A (P = .03; not significant after Bonferroni correction). Because the Akaike information criterion values of the allele (3.333) and risk estimation (odds ratio [OR], 7.43; 95% confidence interval [CI], 0.91-60.95) were also positive, results from single-locus analysis suggest that BPD disease status might be dependent on the effect of the TPH2 C2755A polymorphism. All other SNPs were not associated with the disease, although few positive Akaike information criterion values were obtained from regression analysis for the few polymorphisms in allelic (Table 2) or genotypic (data not shown) levels.
The data from pairwise LD analyses between 6 SNPs in patient and control groups, as represented by the LD coefficient (D′) values,34 are given in Figure 1 and supplementary eTable 5 (available on our Web site). As demonstrated in Figure 1, patients and controls had distinctly different patterns of LD. Despite the SNPs from C2755A to G93329A being in complete LD in both groups (D′, 1), strong LD for SNPs from T − 703G to C2755A was seen only in patients with BPD (D′, 0.87-1). In contrast, complete LD in controls was present only in a shorter region between SNPs T − 703G (rs4570625) and T − 473A (rs11178997). In addition, the D′ values dropped for long-distance SNP pairs (ie, T − 703G and G93329A), and the pattern was similar in both groups. Previously Zhou et al22 used 15 SNPs covering 106 kb of the TPH2 gene to define 2, 3, and 1 haplotype blocks in US white, African American, and southwestern Native American subjects, respectively. In the present study, we defined 2 LD blocks in Han Chinese, and the boundary defined for block 1 is similar to that for the African American population in Zhou et al.22 The distinct LD block patterns across populations may suggest different evolutionary paths for these human lineages.
Haplotype association analysis
The strong LD between T − 703G (rs4570625) to C2755A was present only in patients and not in controls, suggesting a possible disease-related selection, and was further investigated in haplotype association analysis. The SNP C2755A was not used for haplotype construction because of its low population frequency. The test statistics for 3-locus haplotype analysis are given in Table 3. Three major haplotypes accounted for 76% and 96% of the total haplotypes in the controls and patients, respectively, and the data confirmed the strong LD within this region. Except for the major TTA haplotype, which was equally distributed between patients and controls, significant differences in distribution between groups were detected for haplotypes TAG, TAA, and GTG, even after a Bonferroni correction (P < .001). The TAG haplotype was associated with noteworthy risk (OR, 3.73; 95% CI, 1.89-7.36), and the TAA haplotype with significant protective effects (OR, 0.19; 95% CI, 0.08-0.46). Furthermore, P values from an overall comparison of all haplotypes and from permutation tests for disease and gene association were highly significant (P < .001). The statistical analyses for haplotypes constructed with 3 SNPs in the 3′ region or 6 SNPs (supplementary eTable 6 and eTable 7, available on our Web site) showed results similar to those in Table 3; therefore, these data strongly imply the involvement of the TPH2 gene in the development of BPD. Previous studies that used 10 SNPs, covering exons 5 to 7,27 and 5 SNPs covering exons 7 to 928 of the TPH2 gene have detected an association between a particular haplotype and major depression and between a particular haplotype and BPD, respectively. Although their results were significant, the SNPs used in the studies were either intronic or synonymous. This led Zill and colleagues27 to suggest that the true susceptibility locus may lie in other parts of the TPH2 gene, but that it was in LD with the markers these authors used. Our data independently replicate the results of Harvey et al.28 Because the potential functional variants present within the haplotype were located from the promoter region to exon 2, our data suggest that the susceptibility locus for BPD may be located toward the 5′ region of the TPH2 gene. The previous study by Zhou et al22 suggested that a functional locus for major depression in US white and African American populations might be located between introns 5 and 8 of the TPH2 gene. Together, these data suggested that the TPH2 gene is involved in the development of many psychiatric disorders and that the susceptibility locus responsible for various disorders may be different.
Promoter activity in t - 703g and t - 473a reporter constructs
To show the possible effect of promoter SNPs, we used the Transcription Element Search System (http://www.cbil.upenn.edu/cgi-bin/tess/tess) to predict whether the genomic positions of these SNPs are associated with any known transcription factors. The Transcription Element Search System prediction showed that SNPs T − 703G (rs4570625) and T − 473A (rs11178997) were both in the POU3F2 (also known as N-Oct-3) binding sites, and that the base changes at both positions had altered the consensus sequence; therefore, we predicted that it would abolish the binding of POU3F2. We used reporter gene systems with different haplotypes for 2 SNPs to assay promoter activity (Figure 2A). We found that luciferase activities obtained from 3 haplotype constructs containing at least 1 T allele were not significantly different (P = .85 in SH-SY5Y; P = .79 in IMR32) and were consistently more active than the plasmid constructed with the GA haplotype (P = .002 in SH-SY5Y; P = .001 in IMR32). These results suggest that both polymorphisms in the promoter region affect in vitro gene expression.
Determination of pou3f2 binding by electrophoretic mobility shift assay
To investigate whether the difference in the transcriptional activity observed in reporter gene assays is due to alteration in DNA-protein interaction, we used the electrophoretic mobility shift assay to determine whether the POU3F2 protein can bind to the predicted element of the TPH2 promoter. As shown in Figure 2B, 3 major DNA-protein complexes (C1, C2, and C3) were detected with biotin-labeled probes (lanes 2-9 and 11-18) irrespective of polymorphisms. The presence of the 3 complexes were in agreement with previous reports describing 2 additional proteins that can also bind to the POU3F2 sequence motif.35,36 The addition of anti-POU3F2 antibody but not mouse IgG inhibited POU3F2 binding to the probe (especially the C1 complex), indicating that the DNA-protein complex indeed contains POU3F2. Attenuation of the C1 complex by increasing amounts of unlabelled probes (lanes 5-6, 8-9, 14-15, and 17-18) further confirmed that the binding is specific. More important, polymorphisms did affect the binding of POU3F2 to its cis element when lanes 2 and 7 and lanes 11 and 16 were compared. When the C1 complex was contrasted, signals were 30% and 40% lower in lanes 7 and 16, respectively. Results of an additional pull-down assay demonstrated both − 473T and − 703T probes form complexes with POU3F2 in IMR32 cells (supplementary eFigure 2, available on our Web site). These data indicated that the 2 polymorphisms on the TPH2 promoter region can be bound by POU3F2 and that the − 473 T-to-A and − 703 T-to-G alterations indeed reduce POU3F2 protein binding affinity.
Tph enzyme activity in c2755 vs a2755 constructs
The novel coding SNP C2755A identified in this study changed the amino acid sequence from serine to tyrosine at peptide position 41. Although this S41Y polymorphism is on the nonpreserved N-terminal region between the human TPH1 and TPH2 proteins, the sequence alignment of TPH2 proteins in 9 species showed that S41 is highly conserved across species from zebrafish to humans (Figure 3A), but not in the mouse and rat. To analyze the substitution effect on TPH2 function, we transfected plasmids containing full-length TPH2 protein–encoding sequences with 2 alternative alleles into SH-SY5Y cells and measured serotonin production in the supernatants of transfected cells using an enzyme-linked immunosorbent assay for serotonin (Figure 3B). Our data demonstrated, with roughly equal recombinant TPH2 proteins in each assay, that the amount of serotonin in SH-SY5Y cells expressing the TPH2_41Y allele was about 36% (range, 30%-42%) lower than cells expressing the TPH2_41S allele (P = .002). Although we did not measure TPH2 enzyme activity, previous studies in humans38 and mice17 indicated a good correlation between serotonin levels measured in the supernatants of transfected cells and the rate of serotonin synthesis in the cells. The levels of serotonin in SH-SY5Y cells can, therefore, serve as indirect data for TPH2 enzyme activities. These data suggest that the TPH2_41Y allele loses the ability to synthesize serotonin.
Interaction between tph2 and tph1
Expression of mRNAs for TPH1 and TPH2 in the Human Brain
The expression of TPH2 and TPH1 was found to be nonoverlapped in the rat brain18 but coexpressed with different levels using tissue from postmortem human brain.19 To confirm the expression pattern of these 2 proteins in the human brain, we applied nonnormalized complementary DNA panels from human tissue and measured, using Q–RT-PCR, the relative mRNA amount. The measured relative mRNA level was justified with 18S rRNA in each tissue panel and present in a dual y-axis figure (Figure 4). Except in the total brain and cerebellum, where the relative levels of mRNA were similar for TPH1 and TPH2, expression of TPH2 was significantly more abundant in all investigated brain regions, including the amygdala, cerebral cortex, frontal lobes, medulla oblongata (caudal raphe), and pons (rostral raphe). The highest TPH2 mRNA expression was detected in the pons at about 45 times more abundantly than TPH1. These data are consistent with previous reports that TPH2 was the predominant form in the raphe nuclei in the brains of humans19 and rats.18 Results from our study support the overlapped pattern of the 2 TPHs in the human brain and suggest that the 2 genes may be coexpressed in many serotonergic brain regions.
Combined Genotyping Analysis of TPH1 and TPH2
The coexpression of TPH1 and TPH2 proteins in many brain regions but leading to different biological paths suggests a requirement for different regulatory controls. We hypothesized that the 2 proteins share common resources under normal physiological conditions and that they may compete for common factors when the regulatory network becomes unbalanced. To test this hypothesis, we estimated ORs from 205 study participants with TPH112 and TPH2 genotype data (Table 4 and Figure 5). Of those with the TPH2 TAG haplotype, the estimated risk for BPD was 3.73 times greater than that for individuals without the TAG haplotype (95% CI, 1.89-7.36; P < .001). Within the same group, the risk was 4.81 times greater for individuals who also had the TPH1 − 347G allele (95% CI, 2.62-8.83; P < .001). The risk for BPD was only 1.68 times greater (95% CI, 0.72-3.89) if the individual carrying the TPH2 risk haplotype had a TPH1 − 347T allele, and the difference between the 2 subgroups was not significant (P = .22). These data suggested that the effect of TPH2 on BPD etiology can be influenced by the presence of TPH1.
Although TPH2 polymorphisms were associated with many psychiatric disorders, including major depression and BPD,21-28 the functional polymorphisms that represent the susceptibility loci of the TPH2 gene have not yet been identified. Our study, which used single-locus analysis, LD mapping, and haplotype analysis, provided additional support for the association of TPH2 with BPD. In addition, we have demonstrated the function of 3 TPH2 polymorphisms, 2 common variants in the promoter region and 1 novel variant in the exon 2, the combination of which may underlie the risk for BPD.
We have undertaken a systematic approach to screen polymorphisms of the human TPH2 gene in the Han Chinese population. Of the 9 polymorphisms identified in this study, 2 SNPs have not been reported previously, to our knowledge, and thus may represent population-specific polymorphisms in Taiwan. The C2755A polymorphism is the only variant that causes the S41Y substitution. These data are consistent with other reports22,39 in which most of the variants identified from the TPH2 gene were in noncoding regions. Results from these studies point out the conserved nature of the human TPH2 gene. In addition, 3 nonsynonymous variants of L36P,39 P206S,22 and R441H38 are all present at very low frequencies and are not found in our population. Because the R441H variant identified in the US population38 results in approximately 80% loss of function in serotonin production, these data suggest that the S41Y polymorphism identified in our study may also be functional and thus worth further characterization.
The 2 promoter SNPs, T − 703G (rs4570625) and T − 473A (rs11178997), are commonly present in many populations,22 and it was predicted that base changes at both positions can affect POU3F2 binding. Using an in vitro luciferase reporter assay, we demonstrated that significant loss in promoter activity is indeed associated with the GA haplotype. Gel shift assays indicate a reduction in POU3F2 binding affinity of probes containing the − 703G and − 473A alleles when compared with the − 703T- and − 473T-containing probes, respectively. These data provide evidence that the 2 promoter SNPs can affect POU3F2 binding and modulate TPH2 gene expression. A previous study by De Luca et al40 showed abnormal expression of TPH2 mRNA in the brain in patients with BPD. A recent publication36 reported that a polymorphism on a GABRB3 promoter reduces POU3F2 binding affinity and therefore could be a cause for childhood absence epilepsy. Taken together, our data offer additional evidence and suggest that the 2 promoter polymorphisms may represent the common susceptibility loci of the TPH2 gene.
Although TPH1 and TPH2 are homologous proteins and share an overall 71% of amino acid identity,5 the N-terminal end and the regulatory region are quite divergent between these 2 proteins. For example, both TPH1 and TPH2 can be phosphorylated by protein kinase A but TPH2 is phosphorylated at Ser19, a site not even present in TPH1.41 Furthermore, TPH2 is more soluble than TPH1, with a higher molecular weight and different kinetic properties.41 These data strongly suggest that the diverse 5′ region of the TPH2 protein may be important for regulating serotonin production in the central nervous system. The C2755A (S41Y) polymorphism is located on the diverse N-terminal region of human TPH2 proteins, and sequence alignment indicated that the position is highly conserved in the TPH2 proteins across species. Further data using recombinant constructs demonstrated that the TPH2-41Y fusion protein exhibits lower enzyme activity than the TPH2-41S fusion protein. We assume that the reduced enzyme activity is due to a change in the regulatory level, perhaps through the alteration of the 3-dimensional structure in the N-terminal of the TPH2 protein, thus affecting its association with proteins like protein kinase A or 14-3-3. Alternatively, because the S41 residue is predicted on the protein kinase C–binding motif, the polymorphism itself may be an important phosphorylation site for activating TPH2. Thus, a change in this residue might significantly alter the amount of active TPH2 and result in a loss of final serotonin production.
Although SNPs on the promoter region and exon 2 show significant functional alterations in gene expression and TPH2 activity, respectively, their associations with BPD at single-locus levels are not significant. These results indicate that an association of partial functional loss in each SNP with BPD may be weak, but that the combination of several common variants may have a dramatic effect on the etiology of BPD. Our results support the common disease/common variants hypothesis.42 More important, our data confirm the results of Harvey et al28 and suggest that the functional polymorphisms detected in the promoter region and exon 2 might be the susceptibility loci for BPD.
The current conventional wisdom suggests that genes in the same pathway or involved in the same regulatory network tend to have risk interaction.43 Growing evidence also supports the existence of higher-order risk interactions between genes to cause common and complex diseases, for example, the interactions between 5,10-methylene tetrahydrofolate reductase and transcobalamin for spontaneous abortion,44 between the serotonin transporter and dopamine transporter for harm avoidance and reward dependence traits,45 and between prothrombin and factor XIII-A for myocardial infarction.46 The TPH1 and TPH2 proteins are paralogous and have distinct structural and kinetic properties41 but share common resources for serotonin production.5 The dual functions of serotonin correlated well with the expression levels of these 2 enzymes, in which serotonin production as a precursor for melatonin synthesis (the pineal gland projection area) or for central nervous system neurotransmission (raphe nuclei projection area) is synthesized by TPH1 or TPH2, respectively.5 In the present study, we applied Q–RT-PCR to demonstrate that both TPH1 and TPH2 mRNA are expressed in all brain regions under investigation. Expression of TPH2 is about 20 to 45 times greater than that of TPH1 in the raphe nuclei and about 10 to 12 times greater in the raphe nuclei projecting area. The detection of TPH1 in a few brain regions involved in serotonergic neurotransmission (eg, the pons and medulla oblongata) may come from pineal gland projection.47,48 Based on this observation, we propose a competing hypothesis for the 2 TPH proteins. Under normal conditions, serotonergic brain regions predominantly express TPH2 and lead to the generation of the neurotransmitter serotonin in the central nervous system. The existence of a minimal amount of TPH1 that competes with limited cofactors or cosubtracts to generate serotonin for melatonin synthesis can be tolerated. However, when the function of TPH2 in an individual is impaired, ie, gene expression is decreased or enzyme activity is reduced, the competition between the 2 TPHs becomes harmful and significant. The situation turns worse when TPH1 is present with increased amounts or activity. This hypothesis can be partially supported by the observation in our study that the OR increased when at-risk haplotypes of TPH2 combined with the TPH1 − 346G allele, which exhibits stronger promoter activity.10 To our knowledge, this is the first article to propose the interaction between the 2 TPH genes, and the theory can explain the lack of association in our previous study between TPH1 and BPD12: because the association of TPH1 with BPD is influenced by the genotype of TPH2, examination of TPH1 alone missed the gene effect. Although the model seems reasonable, our data cannot exclude the possibility that the observation may simply indicate that critical behavioral outputs have complex contributions from the brain and periphery; as such, they depend on both genes' actions but not direct molecular interaction. Further experiments are required to confirm the proposed model in in vivo biological systems.
Although the results of our genetic analysis were highly significant, the possibility of false-positive findings cannot be totally excluded because of the sampling limitations (ie, small sample size, confounders with other medical conditions, and subject assessment lacked support by structured interviews). However, we have provided additional statistics and functional analyses to support our conclusion. The discovery of TPH2 has changed scientists' views of TPH1 as the candidate gene for a variety of psychiatric diseases. We believe that the present study redefines the role of TPH1 in the development of BPD. Although its effect is indirect, our data suggest that TPH1 can interact with TPH2 to influence the risk for BPD. We propose that the interactive effect between the 2 TPH genes be considered in future studies of serotonin-related disorders.
Correspondence: H. Sunny Sun, PhD, Institute of Molecular Medicine, National Cheng Kung University Medical College, 1 University Rd, Tainan 70101, Taiwan (hssun@mail.ncku.edu.tw).
Submitted for Publication: July 13, 2006; final revision received October 11, 2006; accepted December 24, 2006.
Author Contributions: Drs Lin and Chao contributed equally to this study.
Financial Disclosure: None reported.
Funding/Support: This study was supported by grants 94-2320-B-006-032 and 94-3112-B-006-003-Y from the National Science Council of Taiwan.
Additional Information: The supplementary “Methods” section, eTables, and eFigures are available at http://140.116.6.160/SUNLab/data/TPH2/.
1.Hoyer
DLubbert
HBruns
C Molecular pharmacology of somatostatin receptors.
Naunyn Schmiedebergs Arch Pharmacol 1994;350
(5)
441- 453
PubMedGoogle ScholarCrossref 2.Van Praag
HM 5-HT–related, anxiety- and/or aggression-driven depression.
Int Clin Psychopharmacol 1994;9
((suppl 1))
5- 6
PubMedGoogle ScholarCrossref 3.Bellivier
FHenry
CSzoke
ASchurhoff
FNosten-Bertrand
MFeingold
JLaunay
J-MLeboyer
MLapianche
J-L Serotonin transporter gene polymorphisms in patients with unipolar or bipolar depression.
Neurosci Lett 1998;255
(3)
143- 146
PubMedGoogle ScholarCrossref 4.Brunner
HGNelen
MBreakefield
XORopers
HHvan Oost
BA Abnormal behavior associated with a point mutation in the structural gene for monoamine oxidase A.
Science 1993;262
(5133)
578- 580
PubMedGoogle ScholarCrossref 7.Bellivier
FLeboyer
MCourtet
PBuresi
CBeaufils
BSamolyk
DAllilaire
JFFeingold
JMallet
JMalafosse
A Association between the tryptophan hydroxylase gene and manic-depressive illness.
Arch Gen Psychiatry 1998;55
(1)
33- 37
PubMedGoogle ScholarCrossref 8.Neiswanger
KZubenko
GSGiles
DEFrank
EKupfer
DJKaplan
BB Linkage and association analysis of chromosomal regions containing genes related to neuroendocrine or serotonin function in families with early-onset, recurrent major depression.
Am J Med Genet 1998;81
(5)
443- 449
PubMedGoogle ScholarCrossref 9.Sun
HSTsai
HWKo
HCChang
FMYeh
TL Association of tryptophan hydroxylase gene polymorphism with depression, anxiety and comorbid depression and anxiety in a population-based sample of postpartum Taiwanese women.
Genes Brain Behav 2004;3
(6)
328- 336
PubMedGoogle ScholarCrossref 10.Sun
HSFann
CSLane
HYChang
YTChang
CJLiu
YLCheng
AT A functional polymorphism in the promoter region of the tryptophan hydroxylase gene is associated with alcohol dependence in one aboriginal group in Taiwan.
Alcohol Clin Exp Res 2005;29
(1)
1- 7
PubMedGoogle ScholarCrossref 11.Nielsen
DAVirkkunen
MLappalainen
JEggert
MBrown
GLLong
JCGoldman
DLinnoila
M A tryptophan hydroxylase gene marker for suicidality and alcoholism.
Arch Gen Psychiatry 1998;55
(7)
593- 602
PubMedGoogle ScholarCrossref 12.Lai
TJWu
CYTsai
HWLin
YMSun
HS Polymorphism screening and haplotype analysis of the tryptophan hydroxylase gene (TPH1) and association with bipolar affective disorder in Taiwan.
BMC Med Genet March 31, 2005;614
PubMedGoogle ScholarCrossref 13.Rujescu
DGiegling
ISato
THartmann
AMMoller
HJ Genetic variations in tryptophan hydroxylase in suicidal behavior: analysis and meta-analysis.
Biol Psychiatry 2003;54
(4)
465- 473
PubMedGoogle ScholarCrossref 14.Tang
GRen
DXin
RQian
YWang
DJiang
S Lack of association between the tryptophan hydroxylase gene A218C polymorphism and attention-deficit hyperactivity disorder in Chinese Han population.
Am J Med Genet 2001;105
(6)
485- 488
PubMedGoogle ScholarCrossref 15.Walther
DJPeter
JUBashammakh
SHortnagl
HVoits
MFink
HBader
M Synthesis of serotonin by a second tryptophan hydroxylase isoform.
Science 2003;299
(5603)
76
PubMedGoogle ScholarCrossref 16.Côté
FThevenot
EFligny
CFromes
YDarmon
MRipoche
MABayard
EHanoun
NSaurini
FLechat
PDandolo
LHamon
MMallet
JVodjdani
G Disruption of the nonneuronal tph1 gene demonstrates the importance of peripheral serotonin in cardiac function.
Proc Natl Acad Sci U S A 2003;100
(23)
13525- 13530
PubMedGoogle ScholarCrossref 17.Zhang
XBeaulieu
JMSotnikova
TDGainetdinov
RRCaron
MG Tryptophan hydroxylase-2 controls brain serotonin synthesis.
Science 2004;305
(5681)
217
PubMedGoogle ScholarCrossref 18.Patel
PDPontrello
CBurke
S Robust and tissue-specific expression of TPH2 versus TPH1 in rat raphe and pineal gland.
Biol Psychiatry 2004;55
(4)
428- 433
PubMedGoogle ScholarCrossref 19.Zill
PBüttner
AEisenmenger
WMöller
HJAckenheil
MBondy
B Analysis of tryptophan hydroxylase I and II mRNA expression in the human brain: a post-mortem study.
J Psychiatr Res 2007;41
(1-2)
168- 173
PubMedGoogle ScholarCrossref 20.Sakowski
SAGeddes
TJThomas
DMLevi
EHatfield
JSKuhn
DM Differential tissue distribution of tryptophan hydroxylase isoforms 1 and 2 as revealed with monospecific antibodies.
Brain Res 2006;1085
(1)
11- 18
PubMedGoogle ScholarCrossref 21.Mössner
RWalitza
SGeller
FScherag
AGutknecht
LJacob
CBogusch
LRemschmidt
HSimons
MHerpertz-Dahlmann
BFleischhaker
CSchulz
EWarnke
AHinney
AWewetzer
CLesch
KP Transmission disequilibrium of polymorphic variants in the tryptophan hydroxylase-2 gene in children and adolescents with obsessive-compulsive disorder.
Int J Neuropsychopharmacol 2006;9
(4)
437- 442
PubMedGoogle ScholarCrossref 22.Zhou
ZRoy
ALipsky
RKuchipudi
KZhu
GTaubman
JEnoch
MAVirkkunen
MGoldman
D Haplotype-based linkage of tryptophan hydroxylase 2 to suicide attempt, major depression, and cerebrospinal fluid 5-hydroxyindoleacetic acid in 4 populations.
Arch Gen Psychiatry 2005;62
(10)
1109- 1118
PubMedGoogle ScholarCrossref 23.Zill
PButtner
AEisenmenger
WMoller
HJBondy
BAckenheil
M Single nucleotide polymorphism and haplotype analysis of a novel tryptophan hydroxylase isoform (TPH2) gene in suicide victims.
Biol Psychiatry 2004;56
(8)
581- 586
PubMedGoogle ScholarCrossref 24.Walitza
SRenner
TJDempfle
AKonrad
KWewetzer
CHalbach
AHerpertz-Dahlmann
BRemschmidt
HSmidt
JLinder
MFlierl
LKnolker
UFriedel
SSchafer
HGross
CHebebrand
JWarnke
ALesch
KP Transmission disequilibrium of polymorphic variants in the tryptophan hydroxylase-2 gene in attention-deficit/hyperactivity disorder.
Mol Psychiatry 2005;10
(12)
1126- 1132
PubMedGoogle ScholarCrossref 25.Sheehan
KLowe
NKirley
AMullins
CFitzgerald
MGill
MHawi
Z Tryptophan hydroxylase 2 (TPH2) gene variants associated with ADHD.
Mol Psychiatry 2005;10
(10)
944- 949
PubMedGoogle ScholarCrossref 26.Coon
HDunn
DLainhart
JMiller
JHamil
CBattaglia
ATancredi
RLeppert
MFWeiss
RMcMahon
W Possible association between autism and variants in the brain-expressed tryptophan hydroxylase gene (TPH2).
Am J Med Genet B Neuropsychiatr Genet 2005;135
(1)
42- 46
PubMedGoogle ScholarCrossref 27.Zill
PBaghai
TCZwanzger
PSchule
CEser
DRupprecht
RMoller
HJBondy
BAckenheil
M SNP and haplotype analysis of a novel tryptophan hydroxylase isoform (TPH2) gene provide evidence for association with major depression.
Mol Psychiatry 2004;9
(11)
1030- 1036
PubMedGoogle ScholarCrossref 28.Harvey
MShink
ETremblay
MGagné
BRaymond
CLabbé
MWalther
DJBader
MBarden
N Support for the involvement of TPH2 gene in affective disorders.
Mol Psychiatry 2004;9
(11)
980- 981
PubMedGoogle ScholarCrossref 29.Bach-Mizrachi
HUnderwood
MDKassir
SABakalian
MJSibille
ETamir
HMann
JJArango
V Neuronal tryptophan hydroxylase mRNA expression in the human dorsal and median raphe nuclei: major depression and suicide.
Neuropsychopharmacology 2006;31
(4)
814- 824
PubMedGoogle ScholarCrossref 30.De Luca
VLikhodi
OVan Tol
HHKennedy
JLWong
AH Tryptophan hydroxylase 2 gene expression and promoter polymorphisms in bipolar disorder and schizophrenia.
Psychopharmacology (Berl) 2005;183
(3)
378- 382
PubMedGoogle ScholarCrossref 31.Sun
HSWang
HCLai
TJWang
TJLi
CM Sequence variants and haplotype analysis of serotonin transporter gene and association with bipolar affective disorder in Taiwan.
Pharmacogenetics 2004;14
(3)
173- 179
PubMedGoogle ScholarCrossref 32.Lin
YMJYang
HCLai
TJFann
CSSun
HS Receptor mediated effect of serotonergic transmission in patients with bipolar affective disorder.
J Med Genet 2003;40
(10)
781- 786
PubMedGoogle ScholarCrossref 35.Schreiber
ETobler
AMalipiero
USchaffner
WFontana
A cDNA cloning of human N-Oct3, a nervous-system specific POU domain transcription factor binding to the octamer DNA motif.
Nucleic Acids Res 1993;21
(2)
253- 258
PubMedGoogle ScholarCrossref 36.Urak
LFeucht
MFathi
NHornik
KFuchs
KA GABRB3 promoter haplotype associated with childhood absence epilepsy impairs transcriptional activity.
Hum Mol Genet 2006;15
(16)
2533- 2541[published correction appears in
Hum Mol Genet. 2006;15(21):3272].
PubMedGoogle ScholarCrossref 38.Zhang
XGainetdinov
RRBeaulieu
JMSotnikova
TDBurch
LHWilliams
RBSchwartz
DAKrishnan
KRCaron
MG Loss-of-function mutation in tryptophan hydroxylase-2 identified in unipolar major depression.
Neuron 2005;45
(1)
11- 16
PubMedGoogle ScholarCrossref 39.Breidenthal
SEWhite
DJGlatt
CE Identification of genetic variants in the neuronal form of tryptophan hydroxylase (TPH2).
Psychiatr Genet 2004;14
(2)
69- 72
PubMedGoogle ScholarCrossref 40.De Luca
VVoineskos
DWong
GWShinkai
TRothe
CStrauss
JKennedy
JL Promoter polymorphism of second tryptophan hydroxylase isoform (TPH2) in schizophrenia and suicidality.
Psychiatry Res 2005;134
(2)
195- 198
PubMedGoogle ScholarCrossref 41.McKinney
JKnappskog
PMHaavik
J Different properties of the central and peripheral forms of human tryptophan hydroxylase.
J Neurochem 2005;92
(2)
311- 320
PubMedGoogle ScholarCrossref 42.Collins
FSGuyer
MSCharkravarti
A Variations on a theme: cataloging human DNA sequence variation.
Science 1997;278
(5343)
1580- 1581
PubMedGoogle ScholarCrossref 43.Carlson
CSEberle
MAKruglyak
LNickerson
DA Mapping complex disease loci in whole-genome association studies.
Nature 2004;429
(6990)
446- 452
PubMedGoogle ScholarCrossref 44.Zetterberg
HZafiropoulos
ASpandidos
DARymo
LBlennow
K Gene-gene interaction between fetal MTHFR 677C>T and transcobalamin 776C>G polymorphisms in human spontaneous abortion.
Hum Reprod 2003;18
(9)
1948- 1950
PubMedGoogle ScholarCrossref 45.Kim
SJKim
YSLee
HSKim
SYKim
CH An interaction between the serotonin transporter promoter region and dopamine transporter polymorphisms contributes to harm avoidance and reward dependence traits in normal healthy subjects.
J Neural Transm 2006;113
(7)
877- 886
PubMedGoogle ScholarCrossref 46.Butt
CZheng
HRandell
ERobb
DParfrey
PXie
YG Combined carrier status of prothrombin 20210A and factor XIII-A Leu34 alleles as a strong risk factor for myocardial infarction: evidence of a gene-gene interaction.
Blood 2003;101
(8)
3037- 3041
PubMedGoogle ScholarCrossref 47.Míguez
JMMartín
FJAldegunde
M Melatonin effects on serotonin synthesis and metabolism in the striatum, nucleus accumbens, and dorsal and median raphe nuclei of rats.
Neurochem Res 1997;22
(1)
87- 92
PubMedGoogle ScholarCrossref 48.Sallanon
MClaustrat
BTouret
M Presence of melatonin in various cat brainstem nuclei determined by radioimmunoassay.
Acta Endocrinol (Copenh) 1982;101
(2)
161- 165
PubMedGoogle Scholar