- Research article
- Open access
- Published:
Characterization and functional analyses of the human HTR1A gene: 5’ regulatory region modulates gene expression in vitro
BMC Genetics volume 19, Article number: 115 (2018)
Abstract
Background
The serotonin neurotransmitter (5-HT) and its receptors have important roles in neuropsychiatric disorders such as schizophrenia. The aim of this study was to investigate the functional sequences of the 5′ regulation region of the human HTR1A gene to explore the effects on the expression of the 5-HT1A receptor.
Methods
Fourteen recombinant pGL3-basic vectors containing deletion fragments of the HTR1A gene regulatory region were transfected with HEK-293 and SK-N-SH cells. The relative chemiluminescence intensities of different length fragments were analyzed. The JASPAR software was used for the prediction of transcription factors.
Results
In the HEK-293 cells, the relative chemiluminescence intensity of the − 1649 bp to − 1550 bp (ATG + 1) fragment was significantly different. Two inhibitory activity regions were found in the − 1409 bp to − 1381 bp and − 1196 bp to − 1124 bp fragments, which might be bound to the GATA or SOX10 transcription factors as predicted by the JASPAR software. In addition, the fragments located from − 1124 bp to − 1064 bp and from − 908 bp to − 722 bp up-regulated protein expression. Only the sequence from − 1550 bp to − 1409 bp demonstrated a difference in luciferase expression in the both cell lines. According to the results of the 5’-UTR truncated vectors, there was a repression region at the distal end of the 5’-UTR, an enhancer region might be present at the proximal end of the transcription start site.
Conclusions
Although the functional sequences of the HTR1A gene regulatory region were confirmed, the regulatory factors and functional components require further investigation.
Background
Many studies have shown that the serotonin system plays an important role in neuropsychiatric disorders such as epilepsy [1], depression [2, 3], and schizophrenia [4, 5]. Abnormalities of the serotonin neurotransmitter (5-HT) have also been observed in these disorders [6]. Interestingly, it is the serotonin receptors that determine the function of 5-HT. Among the subtypes of the serotonin receptors, the 5-HT1A receptor has attracted much attention because it has the widest distribution and highest concentration in the human brain [7, 8]. Previous studies have shown that certain brain regions, such as the prefrontal cortex and the cerebellar vermis, of patients with schizophrenia had an abnormal quantity of 5-HT1A receptors [9,10,11].
In the neural cells, the basal expression of the 5-HT1A receptor is controlled by a proximal promoter and an upstream inhibitory region in the HTR1A gene 5′ regulatory region [12]. Recent research identified a negative regulatory area containing a novel dual repressor element located in the upstream promoter region (− 1590 bp to − 1517 bp) of the HTR1A gene in the rat. This region, which appears to be essential for the regulation of the 5-HT1A receptor expression, is approximately 80% homologous to the human HTR1A gene [13, 14]. Although some other functional fragments were confirmed in the same study, the localization of the regulatory sequence was inaccurate and imprecise, primarily due to the large range of fragments. Therefore, the functional characterization of the 5′ regulatory region of the HTR1A gene, which encodes the 5-HT1A receptor, is poorly understood.
Additionally, most studies have focused on 2 polymorphisms--rs6295 (−1019G/C) in the promoter region and rs878567 (*287 T/C/G) located in the 3’-UTR of the human HTR1A gene. It has been reported that the G allele of rs6295 up-regulates the expression of the 5-HT1A receptor and reduces 5-HT levels, which may contribute to severe depression [15], schizophrenia [16], and other neuropsychiatric disorders. Our previous research involved in the haplotype study also demonstrated that the G allele of rs6295 increased the gene expression and might contribute to the susceptibility of schizophrenia [17]. Conversely, other researchers have posited that this single nucleotide polymorphism (SNP) exerts its effects by changing the binding capacity of the transcription factors instead of directly regulating gene expression [18]. Another study also suggested that rs6295 did not impact the expression of HTR1A gene [19]. Meanwhile, some reports have proved that rs878567 locus is significantly associated with schizophrenia [16]. Notably, our previous case-control experiment revealed that the human HTR1A gene polymorphism rs113195492 (−1068G/A) G allele might be a susceptibility gene in female paranoid schizophrenia [20]. However, the subsequent haplotype did not confirm that this polymorphism had a significant impact on protein expression [17].
At present, studies on the functional sequences and polymorphisms of the HTR1A gene 5′-regulatory region involved in schizophrenia were not comprehensive and the conclusions remained controversial. Therefore, we constructed recombinant vectors containing different length sequences in the HTR1A promoter and 5’-UTR regions: we used the dual luciferase reporter assay to explore specific functional areas and their effects on gene expression.
Methods
Construction of target-fragment pGL3 recombinant vectors
Primers (Table 1) and the introduction of BglII, HindIII, or NheI restriction endonuclease cleavage sites at the 5′ end of the corresponding primers were used to amplify the target fragments. The longest target fragment was located in the human HTR1A gene from − 1902 bp to + 99 bp (ATG + 1). The fragment was used as an amplification template for other truncated bodies to synthesize 12 recombinant vectors with 5′ deletion. Then, the fragment located from − 1124 bp to + 99 bp (ATG + 1) was used to amplify the 2 target sequences missing the 3′ end (Although the two target fragments were truncated from the 3′ end, the transcription start site remained [21]). The purified target genes were cloned into a pGM-T vector using the pGM-T Ligation® Kit (TIANGEN, Biotech, Beijing, China), and the recombinant vectors were screened by transformation and sequencing. Finally, pGM-T vectors were subcloned into pGL3 vectors (Promega, Madison, Wisconsin, USA; Catalog number: E1751), and the recombined pGL3 vectors were selected in the same manner to complete eukaryotic cells transfection [17].
Cell culture
A total of 14 recombinant vectors were transfected into human embryonic kidney cell HEK-293 and human neuroblastoma cell SK-N-SH lines (Human embryonic kidney cell HEK-293 and human neuroblastoma cell SK-N-SH were purchased from the Chinese Academy of Sciences cell bank and catalog numbers were GNHu43 and TCHu51, respectively.). HEK-293 cells were incubated in HyClone® DMEM high-glucose medium containing 10% fetal bovine serum (Thermo Fisher Scientific, Waltham, Massachusetts, USA); they were maintained in 5% CO2 + 95% mixed air at 37 °C. SK-N-SH cells were treated with KeyGEN BioTECH® DMEM high-glucose medium with 0.110 g/L sodium pyruvate containing 15% fetal bovine serum under the same conditions. Transfection experiments were performed when the cell density reached 90% or more [17].
Transient transfection and dual luciferase reporter analysis
Cells with densities of 90% were inoculated into 24-well plates at a concentration of 2 × 105 cells per well. After incubation for 36 to 48 h, 14 recombinant vectors were co-transfected with the internal control Renilla luciferase-expressing vector (The internal control vector pRL-TK (Promega) minimizes the experimental variability caused by the differences in cell type and quantity, transfection efficiency or cell lysis efficiency.) using the Lipofectamine® 2000 reagent (Invitrogen, Carlsbad, California, USA). The cells were incubated for 24 h, and then the cell lysates were collected and the expressions of the Firefly luciferase (LUC) and the Renilla luciferase (TK) were detected. Each sample had three replicate wells in each experiment. Need to repeat three complete experiments [17].
Statistical analysis
The LUC/TK value represents the relative fluorescence intensity. Nine LUC/TK values were obtained for each recombinant vector. The relative chemiluminescence intensity of each recombinant vector is expressed as a mean ± standard deviation. The one-way analysis of variance (ANOVA) and the least significant difference (LSD)-T test were performed using SPSS version 20.0 software. One-way ANOVA was used to compare the differences among all vectors. The LSD-T test was only used to compare adjacent fragments, which could reduce the occurrence of the I type of error. A p-value less than 0.05 represented a statistically significant difference in relative chemiluminescence intensities [22].
Transcription factor prediction
Bioinformatics analysis was performed on functional fragments that had a significant effect on protein expression. JASPAR software was used to predict the transcription factors (http://jaspar.genereg.net/cgi-bin/jaspar_db.pl) [23]. A transcription factor and DNA sequence matching degree of 85% was considered perfect.
Results
Construction of pGL3 recombinant vectors with different length target fragments in the 5′ regulatory region of the HTR1A gene
A total of 14 different length sequences of the HTR1A gene were cloned into pGL3-Basic vectors. Details of the inserted fragments are presented in Figs. 1 and 2.
Comparison and analysis of relative fluorescence intensities of 14 recombined vectors of the HTR1A gene regulatory region
A total of 14 recombinant plasmids were transfected into 2 cell lines to verify the roles of different sequences in the HTR1A gene regulatory region. In the HEK-293 cells, the fragment − 1124 bp to + 99 bp (ATG + 1) had the highest relative chemiluminescence intensity. The relative chemiluminescence intensities of all recombinant vectors containing the target fragments were statistically significant with the p = 0.000 (degree of freedom (11, 96); F = 19.294). A comparison of the relative chemiluminescence intensities of 3 pairs of 4 fragments (− 1649 bp to + 99 bp, − 1550 bp to + 99 bp, − 1409 bp to + 99 bp and − 1381 bp to + 99 bp) were significantly different: − 1649 bp to + 99 bp vs. -1550 bp to + 99 bp (p = 0.01); − 1550 bp to + 99 bp vs. -1409 bp to + 99 bp (p = 0.000); − 1409 bp to + 99 bp vs. -1381 bp to + 99 bp (p = 0.000). Comparisons of the relative chemiluminescence intensities among 3 other fragments (− 1196 bp to + 99 bp, − 1124 bp to + 99 bp, and − 1064 bp to + 99 bp) also revealed significant differences: − 1196 bp to + 99 bp vs. -1124 bp to + 99 bp (p = 0.02) and − 1124 bp to + 99 bp vs. -1064 bp to + 99 bp (p = 0.000). Additionally, there was a significant difference between the relative chemiluminescence intensities of the − 908 bp to + 99 bp and − 722 bp to + 99 bp fragments (p = 0.036) (Fig. 3). Finally, the relative chemiluminescence intensities of the 3′ deletion sequences (− 1124 bp to + 99 bp, − 1124 bp to − 119 bp and − 1124 bp to − 372 bp) were also significantly different: − 1124 bp to + 99 bp vs. -1124 bp to − 119 bp (p = 0.000) and − 1124 bp to − 119 bp vs. -1124 bp to − 372 bp (p = 0.000) (Fig. 4).
In the SK-N-SH cell line, the results of one-way analysis of variance showed that the regulatory activities of all target fragments showed significant differences and the p value was 0.000 (degree of freedom (11, 96); F = 7.370). However, the only statistically significant difference was observed between the 5′ deletion fragments − 1550 bp to + 99 bp and − 1409 bp to + 99 bp (p = 0.000) (Fig. 5). Moreover, for the target fragment of 3′ deletion, the relative chemiluminescence intensities of − 1124 bp to + 99 bp, − 1124 bp to − 119 bp and − 1124 bp to − 372 bp fragments were significantly different: − 1124 bp to + 99 bp vs. -1124 bp to − 119 bp (p = 0.000) and − 1124 bp to − 119 bp vs. -1124 bp to − 372 bp (p = 0.000). Additionally, the trend in luciferase expression of the 14 truncated vectors was similar to that of the HEK-293 cells, but the relative chemiluminescence intensities were generally lower in the SK-N-SH cells (Fig. 6).
Prediction and screening of transcription factors
The results of JASPAR software showed that the functional sequences − 1409 bp to − 1381 bp and − 1196 bp to − 1124 bp might be recognized by the transcription factors SOX10, GATA or RUNx family and TFAP2C, which might play a crucial role in the development of the nervous system (Fig. 7). For the sequence − 1124 bp to − 1064 bp including the polymorphism rs113195492, we predicted that when the G allele was converted to an A, the potential binding transcription factors changed from E2F1 and NFIX to PAX5, NFE2L2, and FOS:: JUN (Fig. 8).
Discussion
In order to study the functional sequence of the regulatory region of the HTR1A gene, 14 recombinant pGL3 vectors with different length segments were constructed, transfected into 2 cell lines and analyzed by the luciferase reporter assay. The results showed that, in the HEK-293 and SK-N-SH cell lines, the relative chemiluminescence intensity was significantly increased when the − 1550 bp to − 1409 bp fragment was removed. The luciferase expression of the truncated sequence located from − 1649 bp to − 1550 bp was decreased in both cell lines, but a larger decrease was noted in the HEK-293 cells. These results suggest that there might be an enhancing regulatory factor in the − 1649 bp to − 1550 bp region and that an inhibitory region might be located in the − 1550 bp to − 1409 bp region. As far as we know, there was an inhibitory region of 107 bases or 71 bases (− 1624 bp to − 1517 bp and − 1590 bp to − 1519 bp [ATG + 1]) of the human or rat HTR1A gene: this region including two copies of the dual repressor element (DRE) and a repressor element-1 (RE-1) binds 2 transcription factors (Freud-1 and REST) and plays a potentially important role in the silencing of the HTR1A gene [13, 14]. Therefore, we essentially verified that the sequence located at − 1624 bp to − 1517 bp did posses a suppressor region. Importantly, it appeared that this regulatory area could be reduced to a narrower range (the overlap of the sequence − 1624 bp to − 1517 bp and − 1550 bp to − 1409 bp). Hence, the 33 bases between − 1550 bp and − 1517 bp seem to be the core of the inhibitory region.
By comparing the relative chemiluminescence intensities in the HEK-293 cells, a functional region was identified that had inhibitory activity in the − 1409 bp to − 1124 bp sequence. We then truncated this region further and demonstrated that there were 2 functional inhibitory regions located between − 1409 bp and − 1381 bp and between − 1196 bp and − 1124 bp, which suggests that there might be some inhibitory transcription factors binding at these bases. The JASPAR analysis led to the theory that there may be some transcription factors, such as GATA, SOX10 and RUNX1/X2/X3, binding the sequence from − 1409 bp to − 1381 bp. Meanwhile, the transcription factors most likely to recognize the sequence − 1196 bp to − 1124 bp might be the TFAP2C and GATA families (Fig. 7). A search of the PUBMED database revealed that GATA and SOX10 have been used as common transcription factors to regulate the expression of multiple genes and TFAP2C might play an important role in the development of the nervous system. Related studies have demonstrated that the overexpression of the transcription factor GATA in the hippocampus leaded to a decrease in spinal synaptic density and synaptic-related gene expression, which ultimately resulting in the depressive behaviors [24]. Additionally, GATA also played an important role in the Parkinson’s disease related genes [25]. For the SOX10 and RUNx-related transcription factors, they had a regulatory effect on the peripheral and central nervous system [26, 27]. We hypothesized that the predicted transcription factors might be the essential reasons of regulatory functions of the − 1409 bp to − 1381 bp and − 1196 bp to − 1124 bp sequences. However, subsequent experiments are needed to further prove.
According to our previous haplotype study, the C allele of the rs6295 locus significantly down-regulated protein expression and the G allele of rs113195492 slightly up-regulated gene expression [17]. In this study, we further explored the effects of these 2 SNPs on gene expression by constructing pGL3 recombinant vectors that lacked the rs113195492 G allele or the rs6295 C allele at the 5′ end. In the HEK-293 cells, the relative chemiluminescence intensity of the − 1124 bp to + 99 bp fragment was significantly different than the − 1064 bp to + 99 bp fragment. The deleted region, containing the rs113195492 (−1068G/A) locus, significantly down-regulated protein expression. These results are consistent with our previous study. However, compared to the − 1064 bp to + 99 bp fragment, the deleted region of the rs6295 C allele at the − 908 bp to + 99 bp fragment did not show an effect on protein expression [12]. Some studies have shown that the effect of rs6295 on gene expression results mainly from the combination of transcription factors such as Deaf-1 and Hes5 [28], so we speculated that rs113195492 may also be combine with several transcription factors. Therefore, the possible bindings of the A and G alleles of rs113195492 were predicted separately, one by one. The results showed that when the G allele was converted to an A, the potential binding transcription factors changed from E2F1 and NFIX to PAX5, NFE2L2, and FOS::JUN (Fig. 8). Unfortunately, neither the role of the 2 polymorphic sites nor the role of specific fragment lengths could be proved precisely. Thus, more research is needed to clarify the exact mechanisms of these 2 SNPs on protein expression.
Considering that the 5’-UTR region plays important roles in the stability of mRNA, ribosome transport, and initiation regulation of translation [29], 2 recombinant vectors containing 3′ end deletion fragments (− 1124 bp to − 119 bp and − 1124 bp to − 372 bp) were constructed. Comparison with the − 1124 bp to + 99 bp fragment, which showed the strongest transcriptional activity, revealed that the relative chemiluminescence intensities of the 3 truncated vectors were significantly different between the 2 cell lines. Related studies reported that the − 79 bp, − 115 bp and − 350 bp sites of the HTR1A gene (ATG + 1) bind 2 regulatory factors---Pet-1 and NFκB, while both enhance regulatory effects on the HTR1A gene [6]. However, in this study, the truncated fragment including the 2 elements located in the − 119 bp to + 99 bp fragment showed an inhibitory effect on luciferase expression in the 2 cell lines. We speculated that, on one hand, this conflicting result might be due to differences between species and cell lines, but, on the other hand, some regulatory factors could also have regional and temporal specificity. Related study indicated that Pet-1 transcription factor was highly expressed in the adrenal gland and brain. However, the expression of Pet-1 was lower in the intestines and eyes. Particularly, Pet-1 activated the gene transcription exhibited the beta43’ enhancer and the cell-type dependence [30]. In addition, the expression of Pet-1 in the brain was restricted to the developing or adult 5-HT neurons [31, 32]. Therefore, its effects in other tissues or cells might be quite distinct. By comparing the relative chemiluminescence intensities of the − 1124 bp to − 119 bp fragment and the − 1124 bp to − 372 bp fragment, we confirmed the conclusions of other studies that reported that NFκB enhances gene expression [33, 34]. However, when comparing the effects of − 1124 bp to + 99 bp and − 1124 bp to − 372 bp fragments on luciferase expression, we found that the relative chemiluminescence intensity values were opposite in the 2 cell lines. We believe that the ambiguity of this result might be due to the fact that the regulatory factors in several cell lines might be slightly different.
In this study, we synthesized pGL3 recombinant vectors containing different sequences of the human HTR1A gene regulatory region and more accurately mapped the functional fragments. However, only the base positions of the enhancer regions and the inhibitory regions were determined. The mechanisms of the regulatory factors and functional elements in each region are still unclear. In some regions, transcription factors were also predicted, and subsequent experiments were conducted to demonstrate the binding of these transcription factors and their effects on protein expression. The experimental results were controversial, possibly due to the fact that the experiments were conducted in a single cell lines, which is a limitation to this study.
Conclusions
By constructing pGL3 recombinant vectors with different length sequences of the HTR1A gene regulatory region, the enhancer regulatory region and the inhibitory regions located between − 1649 bp and − 1550 bp and between − 1550 bp and − 1409 bp, respectively, were identified. Down-regulated elements were also present in the fragments from − 1409 bp to − 1381 bp and − 1196 bp to − 1124 bp. The effects of rs113195492 and rs6295 polymorphisms on protein expression were further confirmed by truncation experiments. The influence of the 5’-UTR region of the HTR1A gene on gene expression was also explored. However, follow-up experiments are needed to identify the transcription factors and other regulatory elements that bind to the HTR1A gene sequences [35].
Abbreviations
- ANOVA:
-
One-way analysis of variance
- LSD:
-
Least significant difference
- SNP:
-
Single nucleotide polymorphism
References
Cuellar-Herrera M, Velasco AL, Velasco F, Trejo D, Alonso-Vanegas M, Nuche-Bricaire A, et al. Alterations of 5-HT1A receptor-induced G-protein functional activation and relationship to memory deficits in patients with pharmacoresistant temporal lobe epilepsy. Epilepsy Res. 2014;108:1853–63.
Nugent AC, Bain EE, Carlson PJ, Neumeister A, Bonne O, Carson RE, et al. Reduced post-synaptic serotonin type 1A receptor binding in bipolar depression. Eur Neuropsychopharmacol. 2013;23:822–9.
Albert PR. Transcriptional regulation of the 5-HT1A receptor: implications for mental illness. Philos Trans R Soc Lond Ser B Biol Sci. 2012;367:2402–15.
Terzic T, Kastelic M, Dolzan V, Plesnicar BK. Influence of 5-HT1A and 5-HTTLPR genetic variants on the schizophrenia symptoms and occurrence of treatment-resistant schizophrenia. Neuropsychiatr Dis Treat. 2015;11:453–9.
Takekita Y, Fabbri C, Kato M, Koshikawa Y, Tajika A, Kinoshita T, et al. HTR1A polymorphisms and clinical efficacy of antipsychotic drug treatment in schizophrenia: a meta-analysis. Int J Neuropsychopharmacol. 2016;19(5):1–10.
Albert PR, Fiori LM. Transcriptional dys-regulation in anxiety and major depression: 5-HT1A gene promoter architecture as a therapeutic opportunity. Curr Pharm Des. 2014;20:3738–50.
van Wijngaarden I, Tulp MT, Soudijn W. The concept of selectivity in 5-HT receptor research. Eur J Pharmacol. 1990;188:301–12.
Muller CP, Carey RJ, Huston JP, De Souza Silva MA. Serotonin and psychostimulant addiction: focus on 5-HT1A-receptors. Prog Neurobiol. 2007;81:133–78.
Bantick RA, Deakin JF, Grasby PM. The 5-HT1A receptor in schizophrenia: a promising target for novel atypical neuroleptics? J Psychopharmacol. 2001;15:37–46.
Selvaraj S, Arnone D, Cappai A, Howes O. Alterations in the serotonin system in schizophrenia: a systematic review and meta-analysis of postmortem and molecular imaging studies. Neurosci Biobehav Rev. 2014;45:233–45.
Tauscher J, Kapur S, Verhoeff NP, Hussey DF, Daskalakis ZJ, Tauscher-Wisniewski S, et al. Brain serotonin 5-HT(1A) receptor binding in schizophrenia measured by positron emission tomography and [11C]WAY-100635. Arch Gen Psychiatry. 2002;59:514–20.
Czesak M, Lemonde S, Peterson EA, Rogaeva A, Albert PR. Cell-specific repressor or enhancer activities of Deaf-1 at a serotonin 1A receptor gene polymorphism. J Neurosci. 2006;26:1864–71.
Ou XM, Jafar-Nejad H, Storring JM, Meng JH, Lemonde S, Albert PR. Novel dual repressor elements for neuronal cell-specific transcription of the rat 5-HT1A receptor gene. J Biol Chem. 2000;275:8161–8.
Lemonde S, Rogaeva A, Albert PR. Cell type-dependent recruitment of trichostatin A-sensitive repression of the human 5-HT1A receptor gene. J Neurochem. 2004;88:857–68.
Haslacher H, Michlmayr M, Batmyagmar D, Perkmann T, Ponocny-Seliger E, Scheichenberger V, et al. rs6295 [C]-allele protects against depressive mood in elderly endurance athletes. J Sport Exerc Psychol. 2015;37:637–45.
Guan F, Lin H, Chen G, Li L, Chen T, Liu X, et al. Evaluation of association of common variants in HTR1A and HTR5A with schizophrenia and executive function. Sci Rep. 2016;6:38048.
Wu X, Yao J, Ding M, Shi ZS, Xu FL, Zhang JJ, et al. 5-HT1A receptor (HTR1A) 5′ region haplotypes significantly affect protein expression in vitro. Neurosci Lett. 2017;638:51–4.
Lemonde S, Turecki G, Bakish D, Du L, Hrdina PD, Bown CD, et al. Impaired repression at a 5-hydroxytryptamine 1A receptor gene polymorphism associated with major depression and suicide. J Neurosci. 2003;23:8788–99.
Kishi T, Okochi T, Tsunoka T, Okumura T, Kitajima T, Kawashima K, et al. Serotonin 1A receptor gene, schizophrenia and bipolar disorder: an association study and meta-analysis. Psychiatry Res. 2011;185:20–6.
Zhou X, Ding M, Ding C, Yao J, Pang H, Xing J, et al. Relationship between genetic polymorphisms in the HTR1A gene and paranoid schizophrenia in a northern Han Chinese population. J Mol Neurosci. 2013;49:625–31.
Parks CL, Shenk T. The serotonin 1a receptor gene contains a TATA-less promoter that responds to MAZ and Sp1. J Biol Chem. 1996;271:4417–30.
Yao J, Ding M, Pang H, Xing JX, Xuan JF, Wang BJ. Dopamine D1 receptor (DRD1) 5′ region haplotypes significantly affect transcriptional activity in vitro. Neurosci Lett. 2015;609:120–3.
Krasnov GS, Dmitriev AA, Melnikova NV, Zaretsky AR, Nasedkina TV, Zasedatelev AS, et al. CrossHub: a tool for multi-way analysis of the Cancer genome atlas (TCGA) in the context of gene expression regulation mechanisms. Nucleic Acids Res. 2016;44:e62.
Choi M, Wang SE, Ko SY, Kang HJ, Chae SY, Lee SH, et al. Overexpression of human GATA-1 and GATA-2 interferes with spine formation and produces depressive behavior in rats. PLoS One. 2014;9:e109253.
Sun Y, Ye L, Zheng Y, Yang Z. Identification of crucial genes associated with Parkinson's disease using microarray data. Mol Med Rep. 2018;17:3775–82.
Weider M, Wegner M. SoxE factors: transcriptional regulators of neural differentiation and nervous system development. Semin Cell Dev Biol. 2017;63:35–42.
Wang JW, Stifani S. Roles of Runx genes in nervous system development. Adv Exp Med Biol. 2017;962:103–16.
Savitz J, Lucki I, Drevets WC. 5-HT(1A) receptor function in major depressive disorder. Prog Neurobiol. 2009;88:17–31.
Koundrioukoff S, Polo S, Almouzni G. Interplay between chromatin and cell cycle checkpoints in the context of ATR/ATM-dependent checkpoints. DNA Repair (Amst). 2004;3:969–78.
Fyodorov D, Nelson T, Deneris E. Pet-1, a novel ETS domain factor that can activate neuronal nAchR gene transcription. J Neurobiol. 1998;34:151–63.
Pattyn A, Vallstedt A, Dias JM, Samad OA, Krumlauf R, Rijli FM, et al. Coordinated temporal and spatial control of motor neuron and serotonergic neuron generation from a common pool of CNS progenitors. Genes Dev. 2003;17:729–37.
Craven SE, Lim KC, Ye W, Engel JD, de Sauvage F, Rosenthal A. Gata2 specifies serotonergic neurons downstream of sonic hedgehog. Development. 2004;131:1165–73.
Pierce JW, Lenardo M, Baltimore D. Oligonucleotide that binds nuclear factor NF-kappa B acts as a lymphoid-specific and inducible enhancer element. Proc Natl Acad Sci U S A. 1988;85:1482–6.
Barboric M, Nissen RM, Kanazawa S, Jabrane-Ferrat N, Peterlin BM. NF-kappaB binds P-TEFb to stimulate transcriptional elongation by RNA polymerase II. Mol Cell. 2001;8:327–37.
Press MO, Carlson KD, Queitsch C. The overdue promise of short tandem repeat variation for heritability. Trends Genet. 2014;30:504–12.
Acknowledgements
This manuscript has been edited by native English-speaking experts of BioMed Proofreading.
Funding
This work was supported by the National Natural Science Foundation of China (No. 81601653). The funding is owned and provided by Jun Yao. The funding body involved in the design of the study and the writing of the manuscript.
Availability of data and materials
All data generated or analysed during this study are include in this published article.
Author information
Authors and Affiliations
Contributions
X.W conducted the experimental operation, data summary and the preparation of this article. X.W, F.L.X and J.J.Z researched the study program. B.J.W, J.Y and M.D. gave the experimental design ideas and experimental guidance. All authors read and approved the final manuscript.
Corresponding authors
Ethics declarations
Ethics approval and consent to participate
The experiments complied with the national guidelines and were conducted under the permissions.
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Wu, X., Xu, Fl., Ding, M. et al. Characterization and functional analyses of the human HTR1A gene: 5’ regulatory region modulates gene expression in vitro. BMC Genet 19, 115 (2018). https://doi.org/10.1186/s12863-018-0708-6
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12863-018-0708-6