Background Gilles de la Tourette Syndrome is a neurodevelopmental disorder that is caused by the interaction of environment with a complex genetic background. The genetic etiology of the disorder remains, so far, elusive, although multiple promising leads have been recently reported. The recent implication of the histamine decarboxylase (HDC) gene, the key enzyme in histamine production, raises the intriguing hypothesis of a possible role of histaminergic dysfunction leading to TS onset.
Methods Following up on the finding of a nonsense mutation in a single family with TS, we investigated variation across the HDC gene for association with TS. As a result of a collaborative international effort, we studied a large sample of 520 nuclear families originating from seven European populations (Greek, Hungarian, Italian, Polish, German, Albanian, Spanish) as well as a sample collected in Canada.
Results and Conclusions Interrogating 12 tagging SNPs (tSNP) across the HDC region, we find strong over-transmission of alleles at two SNPs (rs854150 and rs1894236) in the complete sample, as well as a statistically significant associated haplotypes. Analysis of individual populations also reveals signals of association in the Canadian, German and Italian samples. Our results provide strong support for the histaminergic hypothesis in TS etiology and point to a possible role of histamine pathways in neuronal development.
- Other Psychiatry
- Complex traits
Statistics from Altmetric.com
Gilles de la Tourette Syndrome (TS) is a multifactorial neurodevelopmental disorder resulting from the complex interaction of environment with genetic backround. TS is characterised by the appearance of multiple motor and vocal tics, and high comorbidity rates with other neuropsychiatric disorders of childhood such as obsessive compulsive disorder (OCD), attention deficit hyperactivity disorder (ADHD) and autism, suggesting the possibility of a common aetiological background. Due to lack of education of health professionals and the general public, TS remains poorly understood and severely underdiagnosed in many countries despite the fact that recent investigations have reported a prevalence of 0.4–1% in populations of European origin.1 Results from large-scale studies and collaborative efforts have only very recently started to become available (eg, refs. 2–4) and the genetic basis of the disorder remains, so far, elusive, although multiple chromosomal regions have been implicated, providing several promising leads that warrant further investigation.5
Tics are thought to result from dysfunctions in cortical and subcortical regions that are involved in habit formation, including the basal ganglia and related cortical and thalamic structures.6 TS has traditionally been viewed as a disorder of dopaminergic neurotransmission; however, other transmitters and neuromodulators have also been implicated.6 In fact, the recent implication of the l-histidine decarboxylase (HDC) gene in TS aetiology has raised the intriguing hypothesis of the involvement of histaminergic neural pathways in the onset of the disorder.7 HDC is highly conserved throughout different species and catalyses the oxidative decarboxylation of histidine to histamine.8 Ercan–Sencicek et al7 studied a family with one father and eight of his children affected with TS, and found a linkage signal on chromosome 15. Sequencing of all genes within the linked interval led to the discovery of a single rare coding mutation, a premature termination codon (p. W317*, c.951G>A) in exon 9 of the HDC gene. The mutation could not be found in any of the 3000 control individuals of Western European origin who were screened. At the same time, resequencing of the coding region of HDC in 720 patients with TS and 360 controls revealed no additional nonsense variants, demonstrating the fact that the nonsense mutation identified in the index family is extremely rare.
Since the original study that implicated HDC in TS aetiology, Lei et al9 screened the HDC gene for exonic mutations in 100 Chinese Han patients with TS, and could only find three variants, not predicted to result in amino acid changes: a C>T intronic transition (IVS1 +52C>T), which did not affect the splicing site, a synonymous C>A transition (c.426C>A) in exon 4 and a synonymous G>A transition (c.1743G>A) in exon 12.9 The first genomewide association study for TS only produced a weak association signal close to the HDC gene (p value of 0.02 with rs7166052).3 However, a recent genomewide scan for de novo or transmitted rare CNVs in TS found enrichment of genes within histamine receptor signalling pathways.10 Furthermore, a genomewide study of 95 French Canadian trios with familial history of TS11 showed association to a chromosome 15 microsatellite marker (D15S1016) that lies within the same interval found to be linked with TS in the original study that implicated the HDC gene in TS aetiology.7
Following up on the novel hypothesis of the involvement of HDC and histaminergic neural pathways in TS aetiology, we investigated the possible association of variation across the HDC gene in a large sample of families originating from Canada as well as multiple countries in Europe. We find significant association to the HDC gene in our sample, providing strong support for the role of histamine in the genesis and mediation of tics, and pointing to new lines of research in this area.
Material and methods
Genetic association study samples
Five independently collected samples of TS nuclear families were analysed (total of 520 families). Details of age of onset per sample, sex distribution as well as comorbid OCD and ADHD are shown in online supplementary table S1. The Tourette Syndrome Genetics—Southern and Eastern Europe Initiative sample included European-descent families of Polish (36 trios), Italian (50 trios), Hungarian (84 trios), Greek (10 trios) and Albanian origin (6 trios). Assessment was performed by on-site clinicians using the tools provided by the TS Association International Consortium for Genetics.2 TS was ascertained according to DSM-IV-TR criteria for Italy, Hungary, Albania and Greece, and DSM-IV for Poland. Our study also included two independent samples of German descent. The first was the German Collaborative TS Research Group sample (96 trios), collected using DSM-III-R criteria.12 ,13 Our second sample of German ancestry was collected in Munich (69 trios), using DSM-IV-TR criteria. Our fourth European sample was of Spanish ancestry (19 families), collected using DSM-IV-TR criteria. Finally, we also analysed a large sample of TS trios collected in Canada (150 families), using DSM-III-R criteria. Differences between DSM-III-R, DSM-IV and DSM-IV-TR are minimal; the upper age limit of onset is 18 in DSM-IV (and DSM-IV-TR) and 21 in DSM-III-R, and the ‘marked distress’ criterion, possibly pointing to more severe cases, only appears in DSM-IV (applied only to the Polish families).14 For all samples, collection was approved by local Ethics Boards and informed consent was taken from all participating individuals or their parents.
Selection of single nucleotide polymorphisms, genotyping and statistical analysis
Using the HapMap CEPH European population as reference, we selected tagging single nucleotide polymorphisms (tSNPs) in order to capture variation for a 40 kb region including the HDC gene. A total of 12 tSNPs were selected (see online supplementary table S2) capturing variation at an additional nine SNPs with a mean r2 of 0.968. All selected SNPs are intronic. Genotyping was performed using the KASP genotyping method, a competitive allele-specific assay based on Fluorescent Resonance Energy Transfer from LGC Genomics and Kbioscience (http://www.lgcgenomics.com). The average genotyping success was greater than 99%. The linkage disequilibrium (LD) structure of each gene in the different studied populations was analysed using Haploview.15 In order to test for association of the studied SNPs with the TS phenotype, the transmission test for linkage disequilibrium (TDT) was performed, as also implemented in Haploview.15 We report here association results on single SNPs as well as SNP haplotypes of haplotype blocks defined with the Gabriel et al16 criteria (tests with one degree of freedom, as implemented by Haploview). In order to correct for multiple comparisons, 1000 permutations were performed (using Haploview) and the adjusted p values are also reported here (table 1).
Genotyping results for all SNPs conformed to Hardy–Weinberg equilibrium proportions and Mendelian transmission within families. Rare allele frequencies for each analysed sample are shown in online supplementary figure S1. Since the transmission disequilibrium test (TDT) is robust to population stratification, we were able to perform analysis for the complete sample that was available to us, as well as families within each individual population. Single SNP analysis for association with the TS phenotype in the complete sample (520 families) produced strong results of overtransmission of alleles for SNP rs854150 and rs1894236 (table 1). These results remained statistically significant after performing 1000 permutations of the data. Using the Gabriel et al16 criteria for the combined sample, four haplotype blocks were defined across the studied region. The first of the two TS-associated SNPs in our analysis (rs854150) resides within a two-SNP haplotype block, which is also significantly associated with the TS phenotype, both for the protective and the susceptibility alleles (table 1). Significant association is also found with the haplotype block that carries the second TS-associated SNP in our analysis (rs1894236). This is a five-SNP haplotype spanning 9.2 Kb (SNPs rs854151, rs2070596, rs2070595, rs2853766 and rs1894236). The rs1894236 allele that is overtransmitted to TS patients is found on four different haplotypes. On the other hand, the undertransmitted, and thus, protective allele, is only found on a single haplotype at a frequency of 0.21%, which is also significantly undertransmitted to TS patients (permutation p=0.04) (table 1).
Individual population analysis also revealed interesting results that further support a possible role for HDC in TS aetiology (see online supplementary tables S3 and S4). In our combined sample of German origin (our largest individual population sample with 165 trios), SNP rs854150 retained a trend of overtransmission to TS patients (uncorrected p value of 0.04), which, however, did not withstand permutation testing. Nevertheless, an additional association signal was found, including the haplotype block of the two SNPs most distal to our studied region (rs9920021 and rs1365503) with a permutation p value of 0.048. Interestingly, in our Canadian sample (150 families), these two same SNPs produced the strongest signals of association with permutation p values of 0.001 for rs1365503 and 0.021 for rs9920021. SNP rs1894236 was also significantly associated with TS in the Canadian sample. In the Italian population (50 trios), SNP rs7182203 was associated with the TS phenotype (permutation p value of 0.043), as well as a three-site haplotype including this SNP as well as rs854158 and rs854163. Individual analysis of the remaining samples did not produce any remarkable results, which could be due to small individual sample size. We should also note that analysis of the Canadian, German and Italian samples revealed an identical haplotype block structure, while a slightly different structure was observed for the remaining populations that we analysed.
Following up on the recent evidence about the involvement of the HDC gene and histaminergic neural pathways in the aetiology of TS, we studied variation across HDC for association with TS in a large sample of 520 nuclear families from seven European populations (Greece, Hungary, Poland, Italy, Albania, Germany and Spain) as well as a sample collected in Canada. We employed a tSNPs approach and identified strong overtransmission of alleles at two SNPs (rs854150 and rs1894236). Strong association was also found with a five-SNP haplotype carrying the protective rs1894236 allele, as well as two-SNP haplotypes carrying rs854150. Both SNPs that were found to be associated with TS in our study reside in intronic regions of HDC. However, it is worth noting that rs1894236 lies under one of the Encyclopedia of DNA Elements (ENCODE) project peaks for an enhancer-associated histone mark (UCSC genome browser NCBI36/hg18, http://genome.ucsc.edu). The ENCODE project has systematically mapped regions of transcription, transcription factor association, chromatin structure and histone modification, and histone marks are specific chemical modifications in histone sequence which may, for instance, enhance gene transcription.17 Thus, this particular region may be directly involved in regulation of transcription of the HDC gene. Individual population analysis produced additional signals of positive association. Thus, our results strongly support the histaminergic hypothesis for TS development.
Histamine is an evolutionary conserved signalling molecule that plays a central role in gastric acid secretion, innate and acquired immunity and immunomodulation, bronchoconstriction, vasodilation and neurotransmission. The neuronal histaminergic system is involved in a number of basic physiological functions, such as circadian rhythmicity, energy metabolism, neuroendocrine homeostasis, stress, sensory and motor functions, cognition, attention, learning and memory.18 Neuronal histamine is exclusively synthesised in the tuberomamillary nucleus of the hypothalamus (TMN), while histamine-positive fibres have been described projecting to cortex, hippocampus, basal ganglia, thalamus and virtually all other regions of the central nervous system.18 HDC is the key enzyme for histamine production from histidine. In the brain, its mRNA is expressed exclusively in the posterior hypothalamus. The specific regulation of HDC gene expression in the brain is not understood. However, neuroactive peptides and steroids are among the factors that control HDC transcription and degradation in various tissues and contexts (eg, oxidative stress).8 ,18 Recently, diurnal variation in the levels of HDC mRNA in the TMN has been observed in agreement with a role of neuronal histamine in regulating day–night patterns.19
Histamine binds to and acts through four known G protein–coupled receptors (H1 through H4) and a polyamine binding site on glutamatergic N-Methyl-d-Aspartate (NMDA) receptors.18 High densities of H2R and H3R are found in the basal ganglia,18 ,20 especially on the principal neurons of the striatum, the GABAergic medium spiny neurons (MSN).18 ,21 Furthermore, H3R mRNAs in the cortex and in the substantia nigra pars compacta indicate the presence of H3 heteroreceptors on the major inputs to the striatum.20 Thus, histaminergic innervations of anatomical structures that have been related to TS has been confirmed. Furthermore, H3R regulates a variety of neurotransmitters, including dopamine and serotonin. On the other hand, Hdc−/− deficient mice have several traits relevant to features of TS and have shown decreased brain histamine and increased sensitivity to stereotypic behaviours upon administration of dopamine agonists.22 Such stimulant-induced movements, including rearing, sniffing and biting have previously been proposed as a model of human tics.23
So far, the histaminergic system has not received as much attention as other monoaminergic systems of the brain. Classically established as a ‘peripherally’ important mediator of inflammation, the importance of histamine in neurotransmission and its role in neuropsychiatric disorders are only recently starting to be appreciated. For instance, it has been found severely affected in age-related neurodegenerative diseases such as Parkinson disease (PD) and Alzheimer disease (AD), while histamine receptors are becoming a prime target for pharmaceutical research and development and a number of agents have entered into clinical assessment as cognition enhancers for ADHD, schizophrenia, AD and PD.24 Adding to the multiple facets of neural histamine, our findings further support the implication of histaminergic dysfunction in TS aetiology, providing new insight. An interesting next step could be the functional characterisation of RNA transcripts in samples that carry the associated HDC genetic variants in order to help elucidate their role and clarify the possible involvement of histamine pathways in neurodevelopmental disorders.
This study was made possible thanks to the collaborative efforts of Tourette Syndrome Genetics—Southern and Eastern Europe Initiative (TSGeneSEE) and COST Action BM905: European Network for the Study of GTS (EUNETGTS). We would like to acknowledge the following members of these consortia who have also contributed in this work with sample collection (in alphabetical order by country): Mira Kapisyzi, Entela Basha (University Hospital Center ‘Mother Theresa’, Tirana, Albania); Inmaculada Bernal-Bernal, Marta Bonilla-Toribio, Fátima Carrillo, María T Cáceres-Redond, Marta Correa-Vela Ismael Huertas-Fernández, Laura Vargas-González (Unidad de Trastornos del Movimiento, Servicio de Neurología y Neurofisiología Clínica/Sección de Neuropediatría. Instituto de Biomedicina de Sevilla, Hospital Universitario Virgen del Rocío/CSIC/Universidad de Sevilla, Seville, Spain), Pilar Gómez-Garre (Unidad de Trastornos del Movimiento, Servicio de Neurología y Neurofisiología Clínica/Sección de Neuropediatría. Instituto de Biomedicina de Sevilla, Hospital Universitario Virgen del Rocío/CSIC/Universidad de Sevilla, Seville, Spain; Centro de Investigación Biomédica en Red sobre Enfermedades Neurodegenerativas, Spain).
This web only file has been produced by the BMJ Publishing Group from an electronic file supplied by the author(s) and has not been edited for content.
Files in this Data Supplement:
- Data supplement 1 - Online supplement
Contributors IK performed experiments and wrote the paper. SD contributed samples, participated in the design of experiments and analysis and in manuscript preparation. PS contributed samples, participated in the design of experiments and analysis and in manuscript preparation. ZT contributed samples, participated in the design of experiments and analysis and in manuscript preparation. RR contributed samples, participated in the design of experiments and analysis and in manuscript preparation. TW contributed samples, participated in the design of experiments and analysis and in manuscript preparation. MM contributed samples, participated in the design of experiments and analysis and in manuscript preparation. JH contributed samples, participated in the design of experiments and analysis and in manuscript preparation. MMN contributed samples, participated in the design of experiments and analysis and in manuscript preparation. GL contributed samples, participated in the design of experiments and analysis and in manuscript preparation. LF contributed samples, participated in the design of experiments and analysis and in manuscript preparation. PN contributed samples, participated in the design of experiments and analysis and in manuscript preparation. US contributed samples, participated in the design of experiments and analysis and in manuscript preparation. ZA performed experiments and analysis. Vassilios Stathias performed experiments and analysis. CA contributed samples, participated in the design of experiments and analysis and in manuscript preparation. VT contributed samples, participated in the design of experiments and analysis and in manuscript preparation. AK contributed samples, participated in the design of experiments and analysis and in manuscript preparation. CB contributed samples, participated in the design of experiments and analysis and in manuscript preparation. PZ participated in the design of experiments and analysis and in manuscript preparation. PM contributed samples, participated in the design of experiments and analysis and in manuscript preparation. NM contributed samples, participated in the design of experiments and analysis and in manuscript preparation. CB contributed samples, participated in the design of experiments and analysis and in manuscript preparation. PP conceived the study, designed experiments and analysis and prepared the manuscript.
Funding This work was supported by COST Action BM0905 (http://tourette-eu.org) and Tourette Syndrome Association, USA Research Grant Awards (2008–2010).
Competing interests None.
Ethics approval Research Ethics Committee, Democritus University of Thrace.
Provenance and peer review Not commissioned; externally peer reviewed.
Data sharing statement All data are available upon request to the authors and the TSGeneSEE consortium.