Skip to main content

Genes within the serotonergic system are differentially expressed in human brain

Abstract

Background

Serotonin is an important neurotransmitter with wide-ranging functions throughout the central nervous system. There is strong evidence to suggest that regulation of serotonergic gene expression might be related to genetic variability, and several studies have focused on understanding the functional effects of specific polymorphisms within these genes on expression levels. However, the combination of genotype together with gender and brain region could have an overall effect on gene expression. In this study, we report expression patterns of five serotonergic genes (TPH1, TPH2, 5-HT2A, 5-HT2C, 5-HTT) in seven different human post-mortem brain regions (superior frontal gyrus, superior temporal gyrus, striatum, cerebellum, hippocampus, midbrain and thalamus) using TaqManâ„¢ real-time quantitative PCR. In addition, the effect of genotype and gender on their expression levels was determined.

Results

The data revealed that mRNA from the five genes investigated was detected in all brain regions and showed an overall significant difference in expression levels. Furthermore, the expression of 5-HT2C, 5-HT2A and TPH2 was found to be significantly different between the various brain regions. However, neither gender nor genotype showed significant effects on the expression levels of any of the genes assayed. Interestingly, TPH1 and TPH2 were expressed in all brain regions similarly except for within the striatum and cerebellum, where TPH1 was expressed at a significantly higher level than TPH2.

Conclusion

The effect of brain region has a greater influence on serotonergic gene expression than either genotype or gender. These data add to the growing body of evidence that effects of functional polymorphisms on gene expression in vitro are not observed ex vivo, and provide information that will aid in the design of expression studies of the serotonergic gene system within human post-mortem brain.

Background

Serotonin (5-HT) is an important neurotransmitter found in neurons which project to various areas in the central nervous system. Consequently it is associated with several physiological mechanisms [1]. It mediates cellular effects through several proteins that are involved with its neurotransmission, synthesis, metabolism and membrane re-uptake [2, 3]. Disturbances in brain serotonergic systems result in a range of phenotypes such as depression, suicide and anxiety disorders [3–5].

Many functional imaging and binding studies have demonstrated that these disorders can, in part, be explained by differential binding and/or function of components of the serotonergic system such as receptors and enzymes. However, whether these observed differences correspond to differences in gene expression is largely unknown. Regulation of expression occurs via several different processes, but most studies have focused on genetic variation as a potential factor [6]. However, although many polymorphic variations have been demonstrated to elicit changes in expression in vitro (studies in cell cultures), little information on their role in ex vivo investigations using human post-mortem brain tissue exists. Furthermore, despite much attempt to determine a link between 5-HT pathway gene polymorphisms and disease phenotype, few consistent findings have been made including the promoter polymorphism in the 5-HTT gene and anxiety-related traits [7, 8]. Although there are many methodological and biological reasons for this, one explanation might be that differential regulation could be elicited via additional factors such as tissue type and gender as well as background genetic variation.

Some of the most widely studied serotonergic genes include tryptophan hydroxylase 1 and 2 (TPH1 and TPH2; serotonin synthesis), serotonin receptors 2A and 2C (5-HT2A and 5-HT2C; post-synaptic receptors) and serotonin transporter (5-HTT, serotonin re-uptake). TPH is the rate-limiting enzyme in the biosynthetic pathway of serotonin and thus has a major function in regulating the serotonergic system. Studies in mice have revealed the existence of two isoforms of the gene termed as TPH1 and TPH2, with rat and human homologues of the second isoform also cloned [9]. Although the discovery of TPH2 is relatively recent, evidence has emerged in support of its expression in human post-mortem brain. Firstly, TPH2 expression was observed in specific post-mortem brain regions including cortex, thalamus, hippocampus, amygdala and hypothalamus [10]. Secondly, differential expression of TPH1 and TPH2 genes in various post-mortem brain regions has been demonstrated, with TPH2 expression being significantly higher in raphe nuclei [10]. However, TPH1 expression was found to be significantly higher than that of TPH2 in hypothalamus and amygdala. The functional effects of genetic variation in TPH1 and TPH2 genes have had little investigation. However, there is evidence that in the TPH1 gene the TGT haplotype of 3 SNPs (T-1607C, G-1066A, T-346G) within the 5'regulatory region repressed transcriptional activity in human cell lines [11]. Similarly, [12] have demonstrated that a haplotype of 3 SNPs within the TPH2 promoter (G-703T, T-473A, and A90G) also influenced transcription in rat and human cell lines, the TTA haplotype resulting in reduced gene expression.

5-HT2A and 5-HT2C receptors are both expressed post-synaptically on serotonergic neurons. 5-HT2A receptors are found mainly in the cortex and hippocampus [13, 14] whilst 5-HT2C receptors are found in greatest density in the choroid plexus, although they are also located in other brain regions such as the prefrontal cortex [15, 16]. It has been suggested that there is some similarity between the neuronal distribution of 5-HT2A binding sites and 5-HT2A mRNA [14]. More recently, quantitative RT-PCR has revealed more abundant 5-HT2A mRNA in the cortex compared with the hippocampus in rat brain [17]. Interestingly, 5-HT2A receptor density has been shown to be under gender-specific influence, where men demonstrated significantly higher binding capacity of receptors than women [18].

There are a number of polymorphisms in the 5-HT2A gene, some of which have been investigated in terms of their effect on transcription, particularly the two highly-linked SNPs T102C and A-1436G. Both the A and G alleles of the latter have significant basal promoter activity in HeLa cells [19]. More recently, the transcriptional activity of the G allele was found to be significantly lower than the A allele when the more downstream of the two promoters was included in the constructs and in the presence of an enhancer region [20]. This effect was only observed in endogenous 5-HT2A expressing cell-lines, suggesting the presence of cell specific transcription factors is necessary to elicit this functional effect. However, no differences were found in allele-specific expression of the A-1436G SNP in human post-mortem brain tissue [21].

Four variants within the promoter of the 5-HT2C gene have been identified and their influence on gene expression investigated [22]. Two haplotypes each consisting of a particular allele of an upstream GTn microsatellite (Z-6) were found to increase expression of the gene relative to the wild-type. Furthermore, this effect was not cell-specific, suggesting that it would be present in any tissue expressing 5-HT2C.

5-HTT is responsible for re-uptake of serotonin at the pre-synaptic neuron, and is the site of action of several selective serotonin reuptake inhibitors and tricyclic antidepressants. The gene is expressed widely in the brain with highest expression detected in the dorsal raphe nuclei [23]. One of the most documented variants within 5-HTT is the serotonin transporter-linked polymorphic region (5-HTTLPR), a repeat region which exists as two common alleles, the short allele ('S') and long ('L'). There is evidence that this polymorphic region has functional significance, where in human carcinoma cells and lymphoblastoid cells, the L allele confers higher activity than the S allele [24]. More recently, it has been shown that the S and L alleles are themselves variable in that SNPs occur within their repeating units [25]. Another functional VNTR (variable number of tandem repeat) within intron 2 of the 5-HTT gene has revealed differential enhancer activities in the mouse embryo [26] and embryonic stem cells [27].

Here, we focused on first determining whether the mRNA expression levels of TPH1, TPH2, 5-HT2A, 5-HT2C and 5-HTT genes are affected by brain region, genotype or gender. Real-time quantitative PCR (qPCR) was utilised to determine the relative expression levels for these genes in seven different human post-mortem brain regions. The investigated polymorphisms are described above and are located either within, or close to, the postulated promoter region of the genes of interest suggesting a plausible functional effect. Given the extensive literature concerning the significance of TPH1 in disease aetiology, and the potential complication of the existence of its paralogue TPH2, the second aim was to compare the levels of expression of the two loci in the various brain regions.

Results

All genes assayed were expressed in all brain regions tested, and these data are represented in Figure 1. Preliminary analyses of the relative expression across all brain regions (grouping individual 2-ΔCt values) showed an overall significant difference in the mRNA levels of the five genes investigated (p < 0.001 F = 108.782). 5-HT2A has an overall highest expression level, whereas 5-HTT had the lowest expression level.

Figure 1
figure 1

Bar graph representing the expression (2-ΔCt) of five serotonergic genes before controlling for brain region, genotype or gender. All genes are expressed in the brain, 5-HTT showing the lowest expression and 5-HT2A the highest. Bars represent mean and error bars represent standard error of the mean (s.e.m).

Effect of gender on gene expression

Gender did not have any significant effect on the overall brain expression levels of the genes we assayed in this sample (data combined for all seven brain regions). However, expression of most genes was generally higher in males than in females (see Figure 2). Furthermore, there were no significant gender differences on gene expression levels within any of the brain regions we assayed (data not shown).

Figure 2
figure 2

Bar graph representing the effect of gender on expression (2-ΔCt) of the five genes in all brain regions. Males (grey bars) and females (white bars) represent the mean and error bars representing s.e.m. No significant differences in expression are observed.

Effect of brain region on gene expression

There were significant differences in expression levels of three genes within different brain regions. These genes included 5-HT2C (F = 19.297, p < 0.001), 5-HT2A (F = 17.952, p < 0.001) and TPH2 (F = 4.294, p = 0.001), and are represented in the bar charts of Figure 3, 4 and 5, respectively. The significant differences were consistently described by lower expression in the cerebellum compared to other brain regions as shown in Table 1.

Figure 3
figure 3

Bar graph describing the expression (2-ΔCt) of 5-HT2C in various brain regions. There are significant differences in expression (denoted by ** [p < 0.001]) across the brain regions, mainly due to lower expression of genes in the cerebellum. The symbol * represents a significant difference P = 0.003 between 5-HT2C expression in striatum and STG brain regions (see Table 1). Bars represent the mean and error bars represent s.e.m. The groups representing expression in the thalamus have N = 2.

Figure 4
figure 4

Bar graph describing the expression (2-ΔCt) of 5-HT2A in various brain regions. There are significant differences in expression (denoted by ** [p < 0.001] and * [p ≤ 0.04]) across the brain regions, mainly due to lower expression of genes in the cerebellum (see Table 1). Bars represent the mean and error bars represent s.e.m.

Figure 5
figure 5

Bar graph describing the expression (2-ΔCt) of TPH2 in various brain regions. There are significant differences in expression (denoted by ** [p = 0.001] and * [p = 0.011]) across the brain regions, mainly due to lower expression of genes in the cerebellum (see Table 1). Bars represent the mean and error bars represent s.e.m. The groups representing expression in the thalamus have N = 2.

Table 1 Results of comparisons of gene expression levels for those genes showing differences between brain regions

Effect of genotype on gene expression

Genotype frequencies and results of ANOVA of the effects of genotype on gene expression are shown in Table 2. As 5-HT2C is localised to the X-chromosome, we split the sample by gender when assessing genotype effects. TPH1 did not show any variation in genotype (all samples were "GG" after successful digestion, N = 11), and so no further analyses could be performed for this gene. In general, there were no significant effects of genotype on overall brain expression of any of the genes that were analysed. There appears to be a dose-dependant effect on expression of the 5-HT2A gene by A-1438G genotype (mean 2-ΔCt ± S.D.; AA genotype = 0.041 ± 0.048, AG genotype = 0.048 ± 0.052, GG genotype = 0.054 ± 0.065; ANOVA F = 0.02, p = 0.983). However, this failed to reach significance. Due to the small numbers in the groups stratified by genotype and allele frequency, further analysis of the effects of polymorphisms on within-brain region expression could not be reliably performed.

Table 2 Results of expression analyses by genotype

TPH1 vs. TPH2 expression

Overall expression of TPH1 is not significantly different to that of TPH2 in this sample (t = 1.613, p = 0.101, mean 2-ΔCt ± S.D.; TPH1 = 0.001 ± 0.0008, N = 68; TPH2 = 0.001 ± 0.002, N = 63). However, when the data were further analysed for each brain region the expression of TPH1 was significantly higher than TPH2 in the striatum (t = 2.487, p = 0.022) and cerebellum (t = 4.297, p < 0.001) (See Figure 6). There appear to be differences in the expression of TPH1 and TPH2 in the hippocampus and thalamus but as only a small number of data points were obtained for these brain regions, no significant differences could be determined. In addition, no significant differences exist between the relative expression of TPH1 and TPH2 when stratified by gender (data not shown).

Figure 6
figure 6

Bar graph of mean expression (2-ΔCt) of TPH1 (grey bars) versus TPH2 (white bars) and error bars representing s.e.m. Both genes are expressed in all brain regions, and TPH1 expression is significantly higher than TPH2 in both the striatum (t = 2.487, p = 0.022, denoted by *) and cerebellum (t = 4.297, p < 0.001, denoted by **). The groups representing TPH2 expression in the hippocampus and thalamus have a much larger variance than that of TPH1, probably due to the smaller number of individuals in these groups (hippocampus N = 5, thalamus N = 2).

As the possible presence of Alzheimer's disease at death was recorded for 8 individuals, we tested for its possible impact on transcription. However, there was no significant influence of this phenotype on the difference in gene expression levels in all brain regions tested (TPH1, t = 1.33, d.f. = 13, p = 0.21, TPH2, t = 0.62, d.f. = 13, p = 0.55; 5-HT2A, t = 2.17, d.f. = 8, p = 0.06; 5-HT2C, t = 0.85, d.f. = 13, p = 0.41; 5-HTT, t = 2.04, d.f. = 13, p = 0.06).

Discussion

This is the first investigation, to our knowledge, that compares the expression patterns of different serotonergic genes in various human post-mortem brain regions. The results suggest that all five serotonergic genes are expressed in human post-mortem brain tissue and significant overall differences in expression of 5-HT2A, 5-HT2C and TPH2 were observed between the various brain regions analysed. This is in agreement with previous findings that show variation in receptor densities between brain regions [28–30], and suggests that these differences in protein levels might be determined by transcriptional regulation to some degree. Differences in abundance of these various serotonergic transcripts throughout the brain suggest differential gene regulation between specific brain regions. Therefore, it is important to take into consideration the brain region from which the mRNA is derived when performing comparative post-mortem brain expression studies.

Perhaps most surprisingly, there appeared to be little effect of genetic variation on the overall brain expression of the genes we analysed. Although it has been suggested that expression of all five genes is regulated by promoter polymorphisms in vitro [11, 12, 20, 22, 31], we did not observe the same effect ex vivo. Similar results have been observed in previous expression studies in human post-mortem brain. For example, differences in 5-HT2A expression were not observed between the 'A' and 'G' allele of the 5-HT2A A-1438G SNP [21], and allelic ratios of 5-HTT mRNA are not correlated to 5-HTTLPR genotype in human pons [32] or brainstem [33]. Whilst our data are in agreement with these findings, interpretation should be with some caution due the small sample size for each genotype group. This limits our ability to determine any brain region-specific effects of genetic variation on gene expression. Therefore these findings will need to be confirmed by using a larger sample-size in future studies. In addition, other potential confounders that were not measured could have an effect. For example, there has been a report of 5-HT2A imprinting in human fibroblasts and brain [34], although other independent studies found no further evidence for this phenomenon [21, 35]. Overall, the data described here suggest that the postulated functional polymorphisms do not affect gene expression and this may reflect the complexity of gene regulation ex vivo.

No effects of gender on expression levels were observed for any of the genes within this sample. Interestingly, little data exist to date regarding the effect of gender on the expression of human serotonergic genes. However, there is evidence to suggest that gender-specific hormones might exert some regulatory effect on expression in rodents and primates. In Rhesus macaques, administration of oestrogen and progesterone resulted in a decrease in expression of 5-HT2C in the hypothalamus [36], whilst treatment of rats with oestradiol lead to a decrease in 5-HT2C expression in the hippocampus [37]. Subsequently, the expression of both TPH1 and TPH2 was shown to increase in the dorsal raphe nuclei of both mice [38] and macaques [39] following oestradiol administration. In addition, both testosterone and oestrogen were capable of increasing the number of 5-HTT cells in castrated male rats [40]. A limitation of the observations of our studies regarding the situation in humans is that the tissue donors were elderly, and reports of the normal levels of sex hormones within this age group are conflicting [41, 42].

Finally, the results observed here reveal that TPH1 expression is detected in all investigated brain regions at a comparable or higher level than that of TPH2. These data confirm recent findings [10], where TPH1 expression was observed in human brain, and that expression of TPH1 was significantly higher than that of TPH2 in some of the brain regions we also examined (frontal cortex, thalamus, hippocampus, hypothalamus and amygdala). Our data show that TPH1 expression is higher than that of TPH2 in the striatum and cerebellum, but not the hippocampus. This inconsistency could be related to variables we could not measure, or factors such as the differing mean ages of samples analysed in this study and that by Zill and colleagues [10], since there is evidence that increasing age enlarges the heterogeneity of gene expression contributing to greater inter-individual variability in expression [43]. Additionally, the discrepancy in our results with that of Zill and colleagues could be due to the existence of transcriptional variants for TPH2 mRNA in humans. Recent evidence in rats has revealed that TPH2 has two transcriptional variants one of which is expressed at a very low level [44]. The potential of TPH2 mRNA splice variants has not been demonstrated in humans to date and warrants future investigation.

The original discovery of the novel TPH2 isoform led to the hypothesis that it was TPH2 not TPH1 that was the relevant 5-HT catabolic gene in brain tissue based on its much higher expression level (150 times more) in the brain stem of mice [9]. This region contains raphe nuclei which are rich in serotonergic neurons and our results differ from these previous studies as we focused on several other brain areas which did not include the brain stem. Prior to discovery of TPH2, brain regional differences in TPH densities were traditionally explained by differences in translational efficiencies of splice variants [45]. However, this might now be explained by region specific expression of the two TPH isoforms. Our data support this hypothesis as TPH1 is expressed in all brain regions, and that there are brain region specific differences in expression of TPH1 and TPH2. Although emerging data suggest that TPH1 is not sufficient to drive serotonergic metabolism in the brain, it is expressed, albeit at differing levels, and as such does not preclude its role in the regulation of serotonergic functioning. Furthermore, TPH1 in the pineal gland is essential in the synthesis of melatonin, a mechanism by which it can, in turn, influence serotonergic functioning and behaviour [46].

Although we attempted to control for factors that could have an effect on results described here, there are other confounding issues to consider. Firstly, the sample size is small, which might lead to bias in group comparisons. Furthermore, specific conditions at the time of death and post-mortem (e.g. Post mortem interval, PMI) significantly affect RNA integrity and abundance [47–50]. However, these studies show that post-mortem intervals of 48 hours do not affect expression levels. The PMI of the post-mortem samples utilised in this study were 38.67 ± 16.62 hours (mean ± S.D.) and within the limit so that gene expression levels would not be affected.

Real-time qPCR ('TaqMan') was utilised to quantify gene expression since it is a highly reproducible and sensitive technique [51, 52]. The relative quantification method of analysis was employed which compares the expression level of a target gene to the expression of an endogenous control gene. This high-throughput method has been successfully used in previous gene expression studies with post-mortem human brain and biopsies [53, 54]. However, there are some caveats to be made regarding choice of suitable endogenous control gene, especially in regards to expression stability. In this study, SDHA, UBC and GAPDH housekeeping genes were used as endogenous controls based on previous data which showed these genes to be most stably expressed in human neuroblastoma cells [55]. The results did not reveal any significant differences between relative expression values when normalised to each of the endogenous control genes separately (data not shown), thus differential expression of control gene is unlikely to be a confounding factor.

Conclusion

This investigation has demonstrated the differential expression of 5-HT genes in several human post-mortem brain tissue regions using real-time qPCR. These data provide additional information that will be relevant for appropriate sample selection when performing gene expression studies involving human brain. Both gender and transcriptionally-relevant polymorphisms appear to have little influence on ex vivo expression within this sample suggesting that complex interactions present in human brain tissue might eclipse the genetic effects that are consistently observed in vitro. These results, taken together with those from alternative approaches such as animal models, in vitro and neuroimaging studies, will contribute to the understanding of these interactions and inform researchers of the mechanisms that regulate serotonergic gene expression in human brain.

Methods

The sample consisted of brain tissue collected from 15 individuals post-mortem (9 Male, 6 female) provided by the Maudsley Brain Bank (Department of Neuropathology, Institute of Psychiatry, London, UK) and stored at -70°C prior to use. These individuals were normal (N = 7) adults and Alzheimer's patients (N = 8) with agonal states of either cardiac failure, bronchopneumonia, myocardial infarction, hypertension, coronary occlusion, carcinoma of left kidney, ischaemic heart disease, pulmonary embolus. Mean age at death for the sample was 78.0 ± 7.4 years (mean ± s.d.), and mean post-mortem delay was 38.7 ± 16.6 (mean ± s.d.) hours. No significant differences were observed between males and females with regard to either age at death (males: 77.6 ± 7.0 years (mean ± s.d.); females: 79.2 ± 8.6 years (mean ± s.d.) t = -0.381, p = 0.712) or post-mortem delay (males 38.2 ± 14.1 hrs (mean ± s.d.); females 39.3 ± 21.3 hrs (mean ± s.d.); t = -0.112, p = 0.913). From each individual, samples from up to 7 different brain regions were collected and included the following: superior frontal gyrus (SFG; N = 15), superior temporal gyrus (STG; N = 15), hippocampus (N = 7), cerebellum (N = 15), midbrain (N = 9), thalamus (N = 3) and striatum (N = 14). The Ethical Committee (Research) at the Institute of Psychiatry considered and approved this experimental project and the reference study number is 167/00.

Nucleic acid preparation

Brain tissue was stored at -80°C prior to total RNA and DNA extraction using TRI reagent (Sigma, UK) following the manufacturer's instructions. Total RNA was DNase I treated using Qiagen RNeasy kit (Qiagen, UK) and DNase I (Qiagen, UK) following manufacturer's recommendations. cDNA was generated from total RNA using the TaqMan RT-PCR kit (Applied Biosystems, UK). Briefly, 1 μg of total RNA was reverse transcribed with 5.5 mM of MgCl, 500 μM of each dNTP, 2.5 μM random hexamers and 0.4 U/μL RNase inhibitor in a total volume of 20 μl. Samples were incubated at 25°C for 10 minutes, followed by 48°C for 30 minutes and then denatured at 95°C for 5 minutes in a thermocycler.

Relative quantitative real-time PCR using 'TaqMan'

Primers and probes for each gene were obtained via 'Assays-on-demand' (Applied Biosystems, UK). All probes were tested for equal reaction efficiencies prior to use by making a serial dilution of a single cDNA sample (range 1–103 dilution). All probes showed an amplification efficiencies in the accepted range (0.98 – 1.02). Assay ID numbers are listed in Table 3. All probes are designed to span adjacent exons to minimise effects of any residual DNA contamination.

Table 3 Details of qPCR and genotyping assay primers and conditions

TaqMan reactions were performed using 0.5 μl Assays-on-demand primer/probe, 5 μl 2× TaqMan mastermix (Applied Biosystems, UK), 0.5 μl cDNA and appropriate volume of primer and probe in a total volume of 10 μl. Concurrently, three housekeeping endogenous controls were assayed for each individual sample for normalisation purposes (Succinate Deydrogenase (SDHA), Ubiquitin C (UBC) and Glyceraldehyde-3-phosphate dehydrogenase (GAPDH) [55]). Each experiment was performed in triplicate. The TaqMan reaction was carried out using an ABI7900HT sequence detection system set up in absolute quantification mode. Cycling conditions were one cycle of 95°C at 10 minutes, followed by 40 cycles of 95°C at 15 sec and 60°C for 1 min. On completion of PCR, Ct values were generated using SDS 2.1. software.

Genotyping

The following polymorphisms were genotyped: 5-HTTLPR and intron 2 VNTR (5-HTT), T-346G (TPH1), G-703T (TPH2), A-1438G (5-HT2A), GTn microsatellite (5-HT2C). Details of these markers, PCR reaction conditions and detection methods are given in Table 3. PCR reactions consisted of 2.5 μl × 10 PCR reaction buffer IV (Abgene, UK), 200 mM dNTPS (Abgene, UK), 1 unit Taq polymerase (Promega, UK), 10 pmol each primer, appropriate concentration of MgCl2 and 25 ng template DNA up to 25 μl. Cycling conditions consisted of initial denaturation at 94°C for 4 minutes followed by 30 cycles of 95°C for 1 minute, appropriate annealing temperature for 1 minute and 72°C for 1 minute, followed by a final elongation step at 72°C for 4 minutes. PCR reactions were carried out using an MJ research tetrad thermal cycler (MJ research, UK).

Enzyme digests were carried out using Msp I (Promega, UK) or Msl I (NEB, USA) following manufacturer's instructions. Digested samples were subjected to agarose gel electrophoresis using a 2% agarose gel (Abgene, UK) and visualised using UV transillumination (UVP, UK). Those assays requiring capillary gel electrophoresis were analysed using an ABI3100 (Applied Biosystems) set up in genotyping mode. 1 μl PCR product was prepared in 12 μl HI-DI formamide and 0.4 μl GSRox size standard (Applied Biosystems, UK), denatured at 95°C for 10 minutes and snap cooled on ice. Data were analysed using Genotyper V3.6 software (Applied Biosystems, UK).

Statistical analyses

Mean triplicate Ct values for each sample were calculated, and outliers within the triplicates removed using Grubb's method. Briefly, if the obtained Z value (the difference between the investigated Ct and the mean Ct divided by the SD) exceeded the critical value 1.15, the sample was considered an outlier and removed from further analysis. The mean Ct values were used to calculate relative expression using the ΔCt method [56] where normalisation using multiple reference genes was employed. The normalisation index was calculated as the arithmetic mean Ct of the three housekeeping genes (GAPDH, UBC, SDHA). The ΔCt was then calculated by subtracting the normalisation index from the mean target Ct, and then the relative expression is calculated using the formula 2-ΔCt. Relative expression values (2-ΔCt) were log-transformed prior to analysis since these values are not normally distributed [57]. Student's t test or ANOVA with post-hoc Bonferroni correction and Games-Howell test were performed.

The statistical package for social scientists (SPSS, USA) was employed in all analyses, and grouping variables were gender, genotype and brain region.

References

  1. Lucki I: The spectrum of behaviors influenced by serotonin. Biological Psychiatry. 1998, 44: 151-162. 10.1016/S0006-3223(98)00139-5.

    Article  CAS  PubMed  Google Scholar 

  2. Veenstra-VanderWeele J, Anderson GM, Cook EH: Pharmacogenetics and the serotonin system: initial studies and future directions. Eur J Pharmacol. 2000, 410: 165-181. 10.1016/S0014-2999(00)00814-1.

    Article  CAS  PubMed  Google Scholar 

  3. Cravchik A, Goldman D: Neurochemical individuality – Genetic diversity among human dopamine and serotonin receptors and transporters. Archives of General Psychiatry. 2000, 57: 1105-1114. 10.1001/archpsyc.57.12.1105.

    Article  CAS  PubMed  Google Scholar 

  4. Arango V, Huang YY, Underwood MD, Mann JJ: Genetics of the serotonergic system in suicidal behavior. Journal of Psychiatric Research. 2003, 37: 375-386. 10.1016/S0022-3956(03)00048-7.

    Article  PubMed  Google Scholar 

  5. Gingrich JA, Hen R: Dissecting the role of the serotonin system in neuropsychiatric disorders using knockout mice. Psychopharmacology. 2001, 155: 1-10. 10.1007/s002130000573.

    Article  CAS  PubMed  Google Scholar 

  6. D'Souza UM, Craig IW: Functional polymorphisms in dopamine and serotonin pathway genes. Human Mutation. 2006, 27: 1-13. 10.1002/humu.20278.

    Article  PubMed  Google Scholar 

  7. Stoltenberg SF, Burmeister M: Recent progress in psychiatric genetics – some hope but no hype. Human Molecular Genetics. 2000, 9: 927-935. 10.1093/hmg/9.6.927.

    Article  CAS  PubMed  Google Scholar 

  8. Lesch KP, Bengel D, Heils A, Sabol SZ, Greenberg BD, Petri S, et al.: Association of anxiety-related traits with a polymorphism in the serotonin transporter gene regulatory region. Science. 1996, 274: 1527-1531. 10.1126/science.274.5292.1527.

    Article  CAS  PubMed  Google Scholar 

  9. Walther DJ, Bader M: A unique central tryptophan hydroxylase isoform. Biochemical Pharmacology. 2003, 66: 1673-1680. 10.1016/S0006-2952(03)00556-2.

    Article  CAS  PubMed  Google Scholar 

  10. Zill P, Buttner A, Eisenmenger W, Moller HJ, Ackenheil M, Bondy B: Analysis of tryptophan hydroxylase I and II mRNA expression in the human brain: A post-mortem study. Journal of Psychiatric Research. 2007, 41: 168-173. 10.1016/j.jpsychires.2005.05.004.

    Article  PubMed  Google Scholar 

  11. Sun HS, Fann CSJ, Lane HY, Chang YT, Chang CJ, Liu YL, et al.: A functional polymorphism in the promoter region of the tryptophan hydroxylase gene is associated with alcohol dependence in one aboriginal group in Taiwan. Alcoholism-Clinical and Experimental Research. 2005, 29: 1-7. 10.1097/01.ALC.0000150635.51934.6D.

    Article  Google Scholar 

  12. Chen GL, Vallender EJ, Miller GM: Functional characterization of the human TPH2 5 ' regulatory region: untranslated region and polymorphisms modulate gene expression in vitro. Human Genetics. 2008, 122: 645-657. 10.1007/s00439-007-0443-y.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  13. Frazer A, Maayani S, Wolfe BB: Subtypes of Receptors for Serotonin. Annu Rev Pharmacol Toxicol. 1990, 30: 307-348. 10.1146/annurev.pa.30.040190.001515.

    Article  CAS  PubMed  Google Scholar 

  14. Barnes NM, Sharp T: A review of central 5-HT receptors and their function. Neuropharmacology. 1999, 38: 1083-1152. 10.1016/S0028-3908(99)00010-6.

    Article  CAS  PubMed  Google Scholar 

  15. Pompeiano M, Palacios JM, Mengod G: Distribution of the Serotonin 5-Ht2 Receptor Family Messenger-Rnas – Comparison Between 5-Ht(2A) and 5-Ht(2C) Receptors. Molecular Brain Research. 1994, 23: 163-178. 10.1016/0169-328X(94)90223-2.

    Article  CAS  PubMed  Google Scholar 

  16. Abramowski D, Rigo M, Duc D, Hoyer D, Staufenbiel M: Localization of the 5-hydroxytryptamine2C receptor protein in human and rat brain using specific antisera. Neuropharmacology. 1995, 34: 1635-1645. 10.1016/0028-3908(95)00138-7.

    Article  CAS  PubMed  Google Scholar 

  17. Kulikov AV, Naumenko VS, Voronova IP, Tikhonova MA, Popova NK: Quantitative RT-PCR assay of 5-HT1A and 5-HT2A serotonin receptor mRNAs using genomic DNA as an external standard. J Neurosci Methods. 2005, 141: 97-101. 10.1016/j.jneumeth.2004.06.005.

    Article  CAS  PubMed  Google Scholar 

  18. Biver F, Lotstra F, Monclus M, Wikler D, Damhaut P, Mendlewicz J, et al.: Sex difference in 5HT(2) receptor in the living human brain. Neuroscience Letters. 1996, 204: 25-28. 10.1016/0304-3940(96)12307-7.

    Article  CAS  PubMed  Google Scholar 

  19. Spurlock G, Heils A, Holmans P, Williams J, D'Souza UM, Cardno A, et al.: A family based association study of T102C polymorphism in 5HT2A and schizophrenia plus identification of new polymorphisms in the promoter. Molecular Psychiatry. 1998, 3: 42-49. 10.1038/sj.mp.4000342.

    Article  CAS  PubMed  Google Scholar 

  20. Parsons MJ, D'Souza UM, Arranz MJ, Kerwin RW, Makoff AJ: The -1438A/G polymorphism in the 5-hydroxytryptamine type 2A receptor gene affects promoter activity. Biol Psychiatry. 2004, 56: 406-410. 10.1016/j.biopsych.2004.06.020.

    Article  CAS  PubMed  Google Scholar 

  21. Bray NJ, Buckland PR, Hall H, Owen MJ, O'Donovan MC: The serotonin-2A receptor gene locus does not contain common polymorphism affecting mRNA levels in adult brain. Mol Psychiatry. 2004, 9: 109-114. 10.1038/sj.mp.4001366.

    Article  CAS  PubMed  Google Scholar 

  22. Yuan X, Yamada K, Ishiyama-Shigemoto S, Koyama W, Nonaka K: Identification of polymorphic loci in the promoter region of the serotonin 5-HT2C receptor gene and their association with obesity and type II diabetes. Diabetologia. 2000, 43: 373-376. 10.1007/s001250050056.

    Article  CAS  PubMed  Google Scholar 

  23. Ozsarac N, Santha E, Hoffman BJ: Alternative non-coding exons support serotonin transporter mRNA expression in the brain and gut. J Neurochem. 2002, 82: 336-344. 10.1046/j.1471-4159.2002.00964.x.

    Article  CAS  PubMed  Google Scholar 

  24. Heils A, Teufel A, Petri S, Seemann M, Bengel D, Balling U, et al.: Functional promoter and polyadenylation site mapping of the human serotonin (5-HT) transporter gene. J Neural Transm Gen Sect. 1995, 102: 247-254. 10.1007/BF01281159.

    Article  CAS  PubMed  Google Scholar 

  25. Nakamura M, Ueno S, Sano A, Tanabe H: The human serotonin transporter gene linked polymorphism (5-HTTLPR) shows ten novel allelic variants. Mol Psychiatry. 2000, 5: 32-38. 10.1038/sj.mp.4000698.

    Article  CAS  PubMed  Google Scholar 

  26. MacKenzie A, Quinn J: A serotonin transporter gene intron 2 polymorphic region, correlated with affective disorders, has allele-dependent differential enhancer-like properties in the mouse embryo. Proc Natl Acad Sci USA. 1999, 96: 15251-15255. 10.1073/pnas.96.26.15251.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  27. Fiskerstrand CE, Lovejoy EA, Quinn JP: An intronic polymorphic domain often associated with susceptibility to affective disorders has allele dependent differential enhancer activity in embryonic stem cells. FEBS Lett. 1999, 458: 171-174. 10.1016/S0014-5793(99)01150-3.

    Article  CAS  PubMed  Google Scholar 

  28. Heisler LK, Cowley MA, Tecott LH, Fan W, Low MJ, Smart JL, et al.: Activation of central melanocortin pathways by fenfluramine. Science. 2002, 297: 609-611. 10.1126/science.1072327.

    Article  CAS  PubMed  Google Scholar 

  29. Forutan F, Estalji S, Beu M, Nikolaus S, Hamacher K, Coenen HH, et al.: Distribution of 5HT2A receptors in the human brain: comparison of data in vivo and post mortem. Nuklearmedizin. 2002, 41: 197-201.

    CAS  PubMed  Google Scholar 

  30. Marazziti D, Rossi A, Giannaccini G, Zavaglia KM, Dell'Osso L, Lucacchini A, et al.: Distribution and characterization of [3H]mesulergine binding in human brain postmortem. Eur Neuropsychopharmacol. 1999, 10: 21-26. 10.1016/S0924-977X(99)00045-0.

    Article  CAS  PubMed  Google Scholar 

  31. Lesch KP, Bengel D, Heils A, Sabol SZ, Greenberg BD, Petri S, et al.: Association of anxiety-related traits with a polymorphism in the serotonin transporter gene regulatory region. Science. 1996, 274: 1527-1531. 10.1126/science.274.5292.1527.

    Article  CAS  PubMed  Google Scholar 

  32. Lim JE, Papp A, Pinsonneault J, Sadee W, Saffen D: Allelic expression of serotonin transporter (SERT) mRNA in human pons: lack of correlation with the polymorphism SERTLPR. Mol Psychiatry. 2006, 11: 649-662. 10.1038/sj.mp.4001797.

    Article  CAS  PubMed  Google Scholar 

  33. Johnson AD, Zhang Y, Papp AC, Pinsonneault JK, Lim JE, Saffen D, et al.: Polymorphisms affecting gene transcription and mRNA processing in pharmacogenetic candidate genes: detection through allelic expression imbalance in human target tissues. Pharmacogenet Genomics. 2008, 18: 781-791. 10.1097/FPC.0b013e3283050107.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  34. Bunzel R, Blumcke I, Cichon S, Normann S, Schramm J, Propping P, et al.: Polymorphic imprinting of the serotonin-2A (5-HT2A) receptor gene in human adult brain. Brain Res Mol Brain Res. 1998, 59: 90-92. 10.1016/S0169-328X(98)00146-6.

    Article  CAS  PubMed  Google Scholar 

  35. Polesskaya OO, Sokolov BP: Differential expression of the "C" and "T" alleles of the 5-HT2A receptor gene in the temporal cortex of normal individuals and schizophrenics. J Neurosci Res. 2002, 67: 812-822. 10.1002/jnr.10173.

    Article  CAS  PubMed  Google Scholar 

  36. Gundlah C, Pecins-Thompson M, Schutzer WE, Bethea CL: Ovarian steroid effects on serotonin 1A, 2A and 2C receptor mRNA in macaque hypothalamus. Brain Res Mol Brain Res. 1999, 63: 325-339. 10.1016/S0169-328X(98)00295-2.

    Article  CAS  PubMed  Google Scholar 

  37. Birzniece V, Johansson IM, Wang MD, Backstrom T, Olsson T: Ovarian hormone effects on 5-hydroxytryptamine(2A) and 5-hydroxytryptamine(2C) receptor mRNA expression in the ventral hippocampus and frontal cortex of female rats. Neurosci Lett. 2002, 319: 157-161. 10.1016/S0304-3940(01)02570-8.

    Article  CAS  PubMed  Google Scholar 

  38. Gundlah C, Alves SE, Clark JA, Pai LY, Schaeffer JM, Rohrer SP: Estrogen receptor-beta regulates tryptophan hydroxylase-1 expression in the murine midbrain raphe. Biol Psychiatry. 2005, 57: 938-942. 10.1016/j.biopsych.2005.01.014.

    Article  CAS  PubMed  Google Scholar 

  39. Sanchez RL, Reddy AP, Centeno ML, Henderson JA, Bethea CL: A second tryptophan hydroxylase isoform, TPH-2 mRNA, is increased by ovarian steroids in the raphe region of macaques. Brain Res Mol Brain Res. 2005, 135: 194-203. 10.1016/j.molbrainres.2004.12.011.

    Article  CAS  PubMed  Google Scholar 

  40. McQueen JK, Wilson H, Sumner BE, Fink G: Serotonin transporter (SERT) mRNA and binding site densities in male rat brain affected by sex steroids. Brain Res Mol Brain Res. 1999, 63: 241-247. 10.1016/S0169-328X(98)00281-2.

    Article  CAS  PubMed  Google Scholar 

  41. Bjornerem A, Straume B, Midtby M, Fonnebo V, Sundsfjord J, Svartberg J, et al.: Endogenous sex hormones in relation to age, sex, lifestyle factors, and chronic diseases in a general population: the Tromso Study. J Clin Endocrinol Metab. 2004, 89: 6039-6047. 10.1210/jc.2004-0735.

    Article  CAS  PubMed  Google Scholar 

  42. Leifke E, Gorenoi V, Wichers C, Von Zur MA, Von Buren E, Brabant G: Age-related changes of serum sex hormones, insulin-like growth factor-1 and sex-hormone binding globulin levels in men: cross-sectional data from a healthy male cohort. Clin Endocrinol (Oxf). 2000, 53: 689-695. 10.1046/j.1365-2265.2000.01159.x.

    Article  CAS  Google Scholar 

  43. Somel M, Khaitovich P, Bahn S, Paabo S, Lachmann M: Gene expression becomes heterogeneous with age. Curr Biol. 2006, 16: R359-R360. 10.1016/j.cub.2006.04.024.

    Article  CAS  PubMed  Google Scholar 

  44. Abumaria N, Ribic A, Anacker C, Fuchs E, Flugge G: Stress upregulates TPH1 but not TPH2 mRNA in the rat dorsal raphe nucleus: identification of two TPH2 mRNA splice variants. Cell Mol Neurobiol. 2008, 28: 331-342. 10.1007/s10571-007-9259-5.

    Article  CAS  PubMed  Google Scholar 

  45. Boularand S, Darmon MC, Mallet J: The human tryptophan hydroxylase gene. An unusual splicing complexity in the 5'-untranslated region. J Biol Chem. 1995, 270: 3748-3756. 10.1074/jbc.270.8.3748.

    Article  CAS  PubMed  Google Scholar 

  46. Malhotra S, Sawhney G, Pandhi P: The therapeutic potential of melatonin: a review of the science. MedGenMed. 2004, 6 (2): 46.

    PubMed Central  PubMed  Google Scholar 

  47. Barton AJ, Pearson RC, Najlerahim A, Harrison PJ: Pre- and postmortem influences on brain RNA. J Neurochem. 1993, 61: 1-11. 10.1111/j.1471-4159.1993.tb03532.x.

    Article  CAS  PubMed  Google Scholar 

  48. Preece P, Cairns NJ: Quantifying mRNA in postmortem human brain: influence of gender, age at death, postmortem interval, brain pH, agonal state and inter-lobe mRNA variance. Brain Res Mol Brain Res. 2003, 118: 60-71. 10.1016/S0169-328X(03)00337-1.

    Article  CAS  PubMed  Google Scholar 

  49. Bahn S, Augood SJ, Ryan M, Standaert DG, Starkey M, Emson PC: Gene expression profiling in the post-mortem human brain–no cause for dismay. J Chem Neuroanat. 2001, 22: 79-94. 10.1016/S0891-0618(01)00099-0.

    Article  CAS  PubMed  Google Scholar 

  50. Tomita H, Vawter MP, Walsh DM, Evans SJ, Choudary PV, Li J, et al.: Effect of agonal and postmortem factors on gene expression profile: quality control in microarray analyses of postmortem human brain. Biol Psychiatry. 2004, 55: 346-352. 10.1016/j.biopsych.2003.10.013.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  51. Bustin SA: Quantification of mRNA using real-time reverse transcription PCR (RT-PCR): trends and problems. J Mol Endocrinol. 2002, 29: 23-39. 10.1677/jme.0.0290023.

    Article  CAS  PubMed  Google Scholar 

  52. Bustin SA, Benes V, Nolan T, Pfaffl MW: Quantitative real-time RT-PCR–a perspective. J Mol Endocrinol. 2005, 34: 597-601. 10.1677/jme.1.01755.

    Article  CAS  PubMed  Google Scholar 

  53. Medhurst AD, Harrison DC, Read SJ, Campbell CA, Robbins MJ, Pangalos MN: The use of TaqMan RT-PCR assays for semiquantitative analysis of gene expression in CNS tissues and disease models. J Neurosci Methods. 2000, 98: 9-20. 10.1016/S0165-0270(00)00178-3.

    Article  CAS  PubMed  Google Scholar 

  54. Castensson A, Emilsson L, Sundberg R, Jazin E: Decrease of serotonin receptor 2C in schizophrenia brains identified by high-resolution mRNA expression analysis. Biol Psychiatry. 2003, 54: 1212-1221. 10.1016/S0006-3223(03)00526-2.

    Article  CAS  PubMed  Google Scholar 

  55. Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, De Paepe A, et al.: Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. 2002, 3: RESEARCH0034-10.1186/gb-2002-3-7-research0034.

    Article  PubMed Central  PubMed  Google Scholar 

  56. Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001, 25: 402-408. 10.1006/meth.2001.1262.

    Article  CAS  PubMed  Google Scholar 

  57. Bengtsson M, Stahlberg A, Rorsman P, Kubista M: Gene expression profiling in single cells from the pancreatic islets of Langerhans reveals lognormal distribution of mRNA levels. Genome Res. 2005, 15: 1388-1392. 10.1101/gr.3820805.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

The research study was supported by the Medical Research Council (MRC) with an MRC Component Grant (#G0000197). We thank Dr. Michael Parsons for his kind help with the 5-HT2A genotyping. We also acknowledge Clare Wright and Mark Frost for technical assistance with sample preparation.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Ursula M D'Souza.

Additional information

Authors' contributions

KS participated in design of the study, carried out sample preparation, qPCR, statistical analyses and drafted the manuscript. AT participated in statistical analyses. NK performed sample collection. IWC participated in design of the study. UMD conceived of the study, participated in its design and coordination, was involved in statistical analyses and helped draft the manuscript. All authors read and approved the final manuscript.

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Sugden, K., Tichopad, A., Khan, N. et al. Genes within the serotonergic system are differentially expressed in human brain. BMC Neurosci 10, 50 (2009). https://0-doi-org.brum.beds.ac.uk/10.1186/1471-2202-10-50

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/1471-2202-10-50

Keywords