Tourette syndrome (TS) is a model neuropsychiatric disorder thought to arise from abnormal development and/or maintenance of cortico-striato-thalamo-cortical circuits. TS is highly heritable, but its underlying genetic causes are still elusive, and no genome-wide significant loci have been discovered to date. We analyzed a European ancestry sample of 2,434 TS cases and 4,093 ancestry-matched controls for rare (< 1% frequency) copy-number variants (CNVs) using SNP microarray data. We observed an enrichment of global CNV burden that was prominent for large (> 1 Mb), singleton events (OR = 2.28, 95% CI [1.39-3.79], p = 1.2 × 10^-3) and known, pathogenic CNVs (OR = 3.03 [1.85-5.07], p = 1.5 × 10^-5). We also identified two individual, genome-wide significant loci, each conferring a substantial increase in TS risk (NRXN1 deletions, OR = 20.3, 95% CI [2.6-156.2]; CNTN6 duplications, OR = 10.1, 95% CI [2.3-45.4]). Approximately 1% of TS cases carry one of these CNVs, indicating that rare structural variation contributes significantly to the genetic architecture of TS.

Tourette syndrome (TS) is a neuropsychiatric disorder of complex genetic architecture involving multiple interacting genes. Here, we sought to elucidate the pathways that underlie the neurobiology of the disorder through genome-wide analysis. We analyzed genome-wide genotypic data of 3581 individuals with TS and 7682 ancestry-matched controls and investigated associations of TS with sets of genes that are expressed in particular cell types and operate in specific neuronal and glial functions. We employed a self-contained, set-based association method (SBA) as well as a competitive gene set method (MAGMA) using individual-level genotype data to perform a comprehensive investigation of the biological background of TS. Our SBA analysis identified three significant gene sets after Bonferroni correction, implicating ligand-gated ion channel signaling, lymphocytic, and cell adhesion and transsynaptic signaling processes. MAGMA analysis further supported the involvement of the cell adhesion and trans-synaptic signaling gene set. The lymphocytic gene set was driven by variants in FLT3, raising an intriguing hypothesis for the involvement of a neuroinflammatory element in TS pathogenesis. The indications of involvement of ligand-gated ion channel signaling reinforce the role of GABA in TS, while the association of cell adhesion and trans-synaptic signaling gene set provides additional support for the role of adhesion molecules in neuropsychiatric disorders. This study reinforces previous findings but also provides new insights into the neurobiology of TS.

Objective

Tourette's syndrome is polygenic and highly heritable. Genome-wide association study (GWAS) approaches are useful for interrogating the genetic architecture and determinants of Tourette's syndrome and other tic disorders. The authors conducted a GWAS meta-analysis and probed aggregated Tourette's syndrome polygenic risk to test whether Tourette's and related tic disorders have an underlying shared genetic etiology and whether Tourette's polygenic risk scores correlate with worst-ever tic severity and may represent a potential predictor of disease severity.

Methods

GWAS meta-analysis, gene-based association, and genetic enrichment analyses were conducted in 4,819 Tourette's syndrome case subjects and 9,488 control subjects. Replication of top loci was conducted in an independent population-based sample (706 case subjects, 6,068 control subjects). Relationships between Tourette's polygenic risk scores (PRSs), other tic disorders, ascertainment, and tic severity were examined.

Results

GWAS and gene-based analyses identified one genome-wide significant locus within FLT3 on chromosome 13, rs2504235, although this association was not replicated in the population-based sample. Genetic variants spanning evolutionarily conserved regions significantly explained 92.4% of Tourette's syndrome heritability. Tourette's-associated genes were significantly preferentially expressed in dorsolateral prefrontal cortex. Tourette's PRS significantly predicted both Tourette's syndrome and tic spectrum disorders status in the population-based sample. Tourette's PRS also significantly correlated with worst-ever tic severity and was higher in case subjects with a family history of tics than in simplex case subjects.

Conclusions

Modulation of gene expression through noncoding variants, particularly within cortico-striatal circuits, is implicated as a fundamental mechanism in Tourette's syndrome pathogenesis. At a genetic level, tic disorders represent a continuous spectrum of disease, supporting the unification of Tourette's syndrome and other tic disorders in future diagnostic schemata. Tourette's PRSs derived from sufficiently large samples may be useful in the future for predicting conversion of transient tics to chronic tic disorders, as well as tic persistence and lifetime tic severity.