From Genetics to Epigenetics: New Perspectives in Tourette Syndrome Research

Abstract

Gilles de la Tourette Syndrome (TS) is a neurodevelopmental disorder marked by the appearance of multiple involuntary motor and vocal tics. TS presents high comorbidity rates with other disorders such as attention deficit hyperactivity disorder (ADHD) and obsessive compulsive disorder (OCD). TS is highly heritable and has a complex polygenic background. However, environmental factors also play a role in the manifestation of symptoms. Different epigenetic mechanisms may represent the link between these two causalities. Epigenetic regulation has been shown to have an impact in the development of many neuropsychiatric disorders, however very little is known about its effects on Tourette Syndrome. This review provides a summary of the recent findings in the genetic background of TS, followed by an overview on different epigenetic mechanisms, such as DNA methylation, histone modifications and non-coding RNAs in the regulation of gene expression. Epigenetic studies in other neurological and psychiatric disorders are discussed along with the TS-related epigenetic findings available in the literature to date. Moreover, we are proposing that some general epigenetic mechanisms seen in other neuropsychiatric disorders may also play a role in the pathogenesis of TS.

Introduction

Gilles de la Tourette Syndrome (TS) is a neurodevelopmental disorder characterized by one vocal and multiple motor tics, lasting longer than a year (Robertson, 2000). The prevalence of TS is estimated to be ~1% and it occurs more in males than females, with a ratio of 4 to 1, without differences between social classes. Furthermore, in almost 90% of cases, TS arises along with comorbid neuropsychiatric disorders: in 45–60% with obsessive-compulsive disorder (OCD), while in 60% of cases with attention deficit hyperactivity disorder (ADHD). Anxiety, behavioral disorders, autism spectrum disorders, and learning disabilities are also quite common among individuals with TS (Baron-Cohen et al., 1999; Coffey et al., 2000; Kurlan et al., 2002; Burd et al., 2005; Robertson and Orth, 2006; Cavanna and Termine, 2012; Robertson, 2012).

Tics are defined as sudden movements or vocalizations that are recurrent, rapid, arrhythmic, and stereotyped. They decline in situations of distraction and relaxation, while they increase under stress and anxiety. It is generally preceded by premonitory urges or a sense of inner tension that is reduced or relieved by the performance of the tic (Kwak et al., 2003). Motor tics can be classified as (i) simple, involving single muscle or small group of muscles such as the case of blinking, eye rolling, nose twitching; or (ii) complex, requiring a coordinated pattern of movement or sound like touching, squatting, jumping. In the case of vocal tics, sniffing, throat clearing, snorting, gulping, and coughing are classified as simple; while barking, making of animal noises, and uttering strings of words are classified as complex (Robertson, 2012). The most common manifestation is the eye blinking, which is often the first to appear as well. The typical age of onset ranges from 2 to 21 years with a mean at 7 years. Tourette Syndrome is characterized by a waxing and waning course, peak symptom intensity is usually noted in late childhood. The decline of symptoms is usually in late adolescence or early adulthood (Fusco et al., 2006; Stillman et al., 2009).

The Genetic Background of Tourette Syndrome

Like many other neuropsychiatric disorders, TS also has a complex etiology. Several environmental risk factors have been identified, such as perinatal hypoxia, maternal smoking during pregnancy, exposure to androgens, heat, and fatigue (Swain et al., 2007). These environmental factors may interact with many underlying genetic risk factors. To assess the overall genetic contribution in developing a disorder or trait, heritability estimates can be determined, which usually vary between different studies. Based on a recent meta-analysis of all human twin studies to date, the heritability of tics/tic disorders is 0.45 (Polderman et al., 2015). The two most recent studies on the heritability of TS and tic disorders estimated them to be 0.58 and 0.77, respectively (Davis et al., 2013; Mataix-Cols et al., 2015). The risk for first-degree relatives was significantly higher than for second-degree relatives (odds ratios 18.7 and 4.6, respectively). Despite this relatively large overall genetic contribution, most individual gene variants have only very small effects. Hence, most of the earlier studies attempting to unravel the genetic architecture of TS have been hampered by low statistical power due to small sample sizes and clinical heterogeneity. Linkage and candidate gene association studies have identified a number of chromosomal regions and gene polymorphisms possibly implicated in Tourette Syndrome. The candidates for these small scale studies have traditionally been variations in genes involved in the dopaminergic and serotonergic neurotransmitter systems, generally due to their suspected contribution to the etiology of other psychiatric disorders often comorbid with TS, such as ADHD, OCD, autism, etc.

The studied genetic variations are classified in different categories based on the extent of DNA alteration, including single nucleotide polymorphisms (SNPs), as well as shorter or longer repeat variants, such as variable number of tandem repeats (VNTRs) and copy number variations (CNVs). SNPs are the most common form of genetic variation in humans. A SNP represents the change of a single base-pair (bp) in the individual's DNA sequence. The frequency of a particular SNP can vary from private mutations (i.e., only one individual possesses the variation due to a de novo mutation) to very common polymorphisms found in almost half of the population. Millions of SNPs exist in an individual's genome, however, most of them are presumed to be neutral, but some may alter gene expression or cause structural changes in the encoded proteins (and other transcripts). Due to this, many SNPs are known to be associated with various human traits and diseases. VNTRs are polymorphisms in the genome where a short sequence of nucleotides (usually up to 100 bp) is organized in multiple copies, which are clustered together and oriented in the same direction. The copy number can vary between individuals and this may influence gene expression or protein structure. CNVs represent much larger genetic rearrangements (usually thousands to millions of base-pairs in length) and therefore change the copy number of genes within the region due to the deletion or duplication of a chromosomal segment. Gene variants implicated in TS are summarized in Tables 1, 2.

TABLE 1

Table 1. Multiple (positive) findings.

TABLE 2

Table 2. Single (positive) findings.

Most genetic findings implicated in the pathogenesis of Tourette Syndrome to date are inconsistent and studies yielding positive results lack replication in other independent cohorts. On the other hand, some nominally significant candidate gene associations might suffer from lack of proper statistical correction for multiple testing (e.g., Bonferroni correction). However, there are some examples for positive replications as well, these are summarized in Table 1. During the recent years meta-analyses confirmed the role of SLITRK1, NTN4, DRD2, DRD4 and AADAC gene polymorphisms in TS (Liu et al., 2014; Paschou et al., 2014; Yuan et al., 2015; Bertelsen et al., 2016). Here, we will review most of the genes that have been consistently replicated in TS. Single positive findings are summarized in Table 2 and are not discussed in detail.

SLITRK1 is a member of the SLIT and TRK family, a type-I transmembrane protein with an extracellular leucine-rich repeat domain homologous to SLIT (Aruga and Mikoshiba, 2003). It is involved in the control of neurite outgrowth and it is expressed during both embryonic and postnatal development of the cortex, thalamus and the basal ganglia, the neuroanatomical structures believed to be affected in TS (Proenca et al., 2011). The gene coding for SLITRK1 is one of the most studied in relation with Tourette Syndrome. In 2005, a de novo inversion at chromosome 13q33.1 and two additional rare variants were identified in TS patients in the region including a single nucleotide deletion causing a frameshift and a truncated protein and a mutation in the 3′ untranslated region (3′UTR) at the putative binding site for microRNA hsa-miR-189 (Abelson et al., 2005). Since the original discovery several studies in various populations have attempted to replicate the association with controversial results. In a family study of Canadians, a common variation rs9593835 and two haplotypes were found to be associated with the disorder (Miranda et al., 2009). These results were also confirmed in a large sample of European trios with TS (Karagiannidis et al., 2012). A recent study on SLITRK1 found a significant difference in the distribution of haplotypes consisting of SNPs rs9546538, rs9531520, and rs9593835 between Japanese patients and controls (Inai et al., 2015).

Only one genome-wide association study (GWAS) has been published to date in TS by a large international consortium including a sample of 1285 cases and 4964 ancestry-matched controls (Scharf et al., 2013). After correction for multiple testing, none of the half a million studied SNPs reached genome wide significance (p < 5*10⁻⁸), however top hits were enriched for gene variants expressed in the brain and some of them coincide with previous candidate genes. The top hit was rs7868992 (p = 1.85*10⁻⁶), which is located in an intronic region of the COL27A1 gene (collagen, type XXVII, alpha 1). COL27A1 is a fibrillar collagen primarily expressed in cartilage, but it is also expressed in the cerebellum during several stages of human development (Pace et al., 2003), however, its function in the developing nervous system is unknown (Fox, 2008). Another study in 260 Chinese trios assessed the preferential transmission of the rs7868992 and two other COL27A1 gene variants (rs4979357 and rs7868992) by transmission disequilibrium test (TDT) and found the latter two variants nominally significant, however these results did not survive correction for multiple testing (Liu et al., 2015b).

A replication study with the top 42 SNPs of the original GWAS was performed with a sample of over 600 cases and 600 controls and the meta-analysis yielded a top signal at rs2060546 with p = 5.8*10⁻⁷ in proximity of the NTN4 gene on chromosome 12q22, which codes for netrin 4, an axon guidance protein expressed in the developing striatum. Many of the previous findings (26 out of 42) showed a similar trend underlining the reliability of the GWAS hits as true risk factors for TS (Paschou et al., 2014).

Dopamine receptors, especially DRD2 and DRD4 are two of the most widely studied candidate genes in the field of psychiatric genetics. The dopamine D2 receptor (DRD2) located on chromosome 11p23.2 was characterized previously with TaqIA, TaqIB and TaqID SNPs based on restriction digestion of these polymorphic sites with the TaqI restriction endonuclease enzyme (Vereczkei et al., 2013). Later it turned out that the TaqIA cleavage site is located ~10 kilobases downstream from the DRD2 gene, in exon 8 of the ANKK1 (ankyrin repeat and kinase domain containing 1) gene (Neville et al., 2004), which is a member of the serine/threonine kinase family. The TaqIA polymorphism (rs1800497) causes an amino acid change in ANKK1 (Glu713Lys), which seems to have a significant effect on the specificity of substrate binding. The protein product of the ANKK1 gene was considered as a negative regulator of the transcription factor NF-kB (Nuclear Factor-Kappa B) (Meylan and Tschopp, 2005). Moreover, the expression levels of NF-kB-regulated genes were shown to be altered by the TaqIA variant in an in vitro luciferase system (Huang et al., 2009). Since DRD2 is regulated by NF-kB (Fiorentini et al., 2002; Bontempi et al., 2007) it could be assumed that this ANKK1 variant can indirectly affect DRD2 receptor density. It is also possible, however, that the TaqIA SNP is only a marker of other functional DRD2 variants associated with a number of psychiatric disorders, such as the strongly linked TaqIB (Ponce et al., 2009).

A recent meta-analysis on the association of the TaqIA SNP rs1800497 and TS based on a number of previous case-control studies (Table 1) comprising a sample of over 500 cases and 500 controls, as well as TS trios concluded that this variant is indeed a risk factor for the development of the disorder (Yuan et al., 2015).

Another widely studied polymorphic dopamine receptor is thedopamine D4 receptor gene (DRD4) The DRD4 gene located on chromosome 11p15.5 is highly polymorphic, containing over 200 SNPs and several VNTRs. A 48 bp long VNTR ranging from 2 to 11 repeats in exon 3 of the gene changes the length of the third intracellular loop of the receptor (Vereczkei et al., 2013) with a possible effect on downstream signaling efficiency by inhibiting the adenylyl cyclase enzyme (Van Tol et al., 1992). The DRD4 7 repeat allele seems to show decreased sensitivity to dopamine compared to the 4 repeat allele (Asghari et al., 2002) and according to more recent neurobiological findings it does not form heteromers with D2 receptors in the striatum (Borroto-Escuela et al., 2011). The D4 receptor is mainly expressed in cortical and limbic regions in the CNS and carriers of the 7 repeat allele were shown to have higher susceptibility to various addictive behaviors, ADHD, as well as several other psychiatric traits.

A number of candidate gene studies have addressed the issue of possible relevance of the exon 3 VNTR in DRD4 in TS, and a couple of small family studies have found a positive association (Table 1). A combined family and case-control study in a Han Chinese TS population revealed significant transmission disequilibrium for the 2 repeat and the 7 repeat alleles. The results suggest that the 2 repeat allele might play a protective role, while the common 4 repeat may predispose to TS (Liu et al., 2014).

Arylacetamide deacetylase (AADAC) is an enzyme involved in neutral lipid lipolysis, detoxification, and drug metabolism and is mainly expressed in the liver (Quiroga and Lehner, 2011). However, its expression in different brain regions has also been confirmed with a so far unknown physiological relevance. A recent CNV-analysis in a smaller Danish TS cohort followed by a meta-analysis of a large European sample of over 1000 TS patients and 100,000 controls confirmed the role of AADAC deletions in TS pathogenesis (Bertelsen et al., 2016).

The IMMP2L gene located on chromosome 7q31 encodes a protein involved in processing the signal peptide sequences used to target mitochondrial proteins into the mitochondrium. Its association with TS was first discovered in a family with a balanced 7;18 translocation (Boghosian-Sell et al., 1996). A recent CNV study in a Danish cohort of 188 TS patients reported a 5′-end intragenic deletion of IMMP2L in seven of these patients (Bertelsen et al., 2014). Interestingly, in 4 of the 7 cases, the deletion was within intron 3. The frequency of this IMMP2L deletion was significantly higher than that of the control population. Notably, IMMP2L has been implicated in other neuropsychiatric disorders, such as autism and ADHD (International Molecular Genetic Study of Autism Consortium, 1998; Elia et al., 2010).

Histidine Decarboxylase (HDC), located on chromosome 15q21.2, encodes a member of the group II decarboxylase family that converts the amino acid L-histidine into histamine, a biogenic amine that can act as a local mediator released from mast cells during an immune reaction, as well as a monoaminergic neurotransmitter in the CNS. A premature termination codon (W317X) in the HDC gene was discovered in a large multigenerational family, where the father and all eight of his children were affected with TS had the non-sense variant (Ercan-Sencicek et al., 2010). On the other hand, a study of 100 Han Chinese TS patients failed to confirm the association of this non-sense mutation in HDC with the disorder (Lei et al., 2012). However, a subsequent large family study of 520 European trios with Tourette Syndrome transmission disequilibrium test (TDT) found SNPs rs854150 and rs1894236 to be over-transmitted in patients, confirming the role of HDC in the development of TS (Karagiannidis et al., 2013). Interestingly, HDC knockout mice have also been proposed as a genetic animal model of TS. The KO animals exhibited potentiated tic-like stereotypies, recapitulating a core phenomenology of the disorder (Castellan Baldan et al., 2014; Xu et al., 2015).

The NRXN1 gene located on chromosome 2p16.3 codes for an important mediator of cell-cell interactions in the central nervous system and it has been implicated in neuropsychiatric disorders, such as autism and schizophrenia (Vrijenhoek et al., 2008; Glessner et al., 2009). A rearrangement of ~400 Kb within exons 1-3 of NRXN1 was recently found in TS patients (Nag et al., 2013). Notably, this result is also consistent with a previous analysis reporting the presence of a CNV comprising NRXN1 in a Danish TS cohort (Sundaram et al., 2010).

For a more detailed review on the genetic background of TS and the description of candidate genes where no positive replications were published in the last few years, such as BTBD9, SLC6A3 (DAT1) and MAO-A, see recent review papers by (Paschou, 2013; Pauls et al., 2014).

From Genetics to Epigenetics

As described in the previous section, etiology of Tourette Syndrome has a considerable genetic component. However, it is evident that environment also plays an important role in TS, since discordance between monozygotic twins is not rare, while both prenatal (smoke, alcohol abuse, low birth weight, etc), and perinatal (complicated birth) risk factors were reported. Furthermore, it has also been hypothesized that TS might arise as a consequence of autoimmune mechanisms following Group A β-hemolytic streptococcal infection (GABHS) (Swedo et al., 1998; Snider and Swedo, 2003), a condition labeled as pediatric autoimmune neuropsychiatric disorder associated with streptococcal infection (PANDAS). Finally, the course of the disease is also influenced by environmental factors.

As seen earlier, based on family and twin studies the heritability of TS lies approximately between 50 and 80%, which is a figure often seen in case of other psychiatric disorders. However, gene polymorphisms reported by association studies of these disorders rarely account for more than just a fraction of the overall estimated genetic variance. This phenomenon of “missing heritability” may, in part, be explained by various epigenetic mechanisms arising from the dynamic interaction between the environment and an individual's genome.

Cells are reacting to acute and chronic environmental changes by altering their gene expression state, which can be considered as their adaptive reaction. The first step in this adaptive reaction is the modification of the chromatin structure. Chromatin is the macromolecular complex containing DNA and nuclear proteins/histones. The nucleosome is the building unit of chromatin. It consists of two of each of the four core histones (H3, H4, H2A, and H2B) forming an octamer wrapped around twice by DNA (Figure 1). The positively charged N-terminal histone tails are prone to undergo several modifications (acetylation, methylation, phosphorylation, etc.). Similarly, DNA can be methylated. These covalent modifications influence the chromatin structure and lead to changes of the transcriptional state of genes. Even though these changes are dynamic, they can be maintained through cell divisions and from one generation to the next. This mechanism, which does not involve the DNA sequence, but allows the transmission of acquired traits through mitosis and sometimes through a few generations, is called epigenetic inheritance. Another main but chromatin-independent epigenetic mechanism is via non-coding RNAs, which also play a crucial role in gene expression regulation often together with chromatin-related mechanisms.

FIGURE 1

Figure 1. Overview of common covalent epigenetic modifications. A schematic nucleosome and examples of potential epigenetic modifications are shown. The histone octamer is represented in a cylindrical form with one pair of histones H3/H4 indicated. The protruding H3 histone tail and DNA are indicated in orange and purple, respectively. Oppositely, histone and DNA modifications are shown in purple and orange. The functional roles of histone modifications are indicated in colored boxes. The enzyme families catalyzing the modifications are listed in boxes below. The enzymatic links between the different cytosine modifications are shown in the upper right corner of the figure.

Epigenetic Modifications

DNA Methylation

The most common epigenetic modification of DNA is the covalent attachment of a methyl group to the C5 position of cytosines. In mammals methylation occurs almost exclusively in CpG dinucleotides (Ziller et al., 2011). Most of the CpG dinucleotides of the genome are present in regions with low GC abundance (Bird, 1993). These CG dinucleotides are generally methylated in all cell types. Methylation leads to high mutation rate, since the frequent oxidative deamination of 5-methyl-cytosine (5 mC) results in thymine (Antonarakis et al., 2000; Baba et al., 2011). This mutation is inefficiently repaired and thus its transmission rate is high. Therefore, low CpG frequency (10% of the expected) characterizes the non-coding part of the genome and CpG cytosines progressively disappear from these genomic regions. However, short CpG-rich sequences (CpG islands – CGI) are frequently found in gene regulatory regions (Bird et al., 1985; Deaton and Bird, 2011). CGIs often show tissue-specific methylation pattern and are frequently unmethylated (only) in germ cells (Smallwood et al., 2011; Zeng et al., 2014). Thus, the important gene regulatory sequences are preserved from a high mutation rate maintaining their biological function.

DNA methylation is tightly linked to gene expression. The more a regulatory region is methylated the less the gene is expressed. Several gene expression regulating DNA binding proteins (e.g., transcription factors) are sensitive to the methylation of their target sequence. Some of them cannot bind if the sequence is methylated (e.g., CTCF, E2F family, Myc, CREB; Hark et al., 2000; Blattler and Farnham, 2013), while others require methylated DNA for binding (e.g., MeCP2 methyl-CpG-binding protein 2 or Kaiso, a zinc finger protein; Lewis et al., 1992; Prokhortchouk et al., 2001; Arányi et al., 2005; Smith and Meissner, 2013). By the stabilization of different chromatin states, DNA methylation contributes to cell differentiation, cellular memory, X chromosome inactivation, and other nuclear processes. Although initially it was considered to be a static epigenetic mark, DNA methylation dynamically changes (Kangaspeska et al., 2008; Métivier et al., 2008; Yamagata et al., 2012b) and different enzymes ensure the equilibrium state.

DNA methylation is catalyzed by DNA methyl-transferase (DNMT) enzymes, which either establish or maintain the methylation pattern. DNMT1 is a maintenance methyl-transferase and it preferentially methylates DNA methylated only on one strand in order to preserve and reestablish the pattern of methylation after DNA replication (Pradhan et al., 1999; Mohan and Chaillet, 2013). By doing so, during DNA replication DNMT1 is enriched at the replication fork and reproduces the methylation pattern based on the original template strand (Schermelleh et al., 2007). DNMT3A and DNMT3B are de novo methyl-transferases capable of establishing new patterns of methylation. DNA demethylation is mainly performed by members of the ten-eleven-transferase (TET) enzyme family through the hydroxylation of the methyl group. This leads to the formation of 5-hydroxymethyl-cytosine (5 hmC), which is lost after further oxidation by the same enzyme (Hashimoto et al., 2015). While 5 mC nucleotides represent a few percent (typically between 3 and 7%) of all genomic cytosines in cells and cell lines, 5 hmC is much less abundant constituting only 0.01–1% of all cytosines (Globisch et al., 2010). Interestingly, 5 hmC is rather frequent in primary cells and particularly in the brain (Globisch et al., 2010). Different data also indicate that 5 hmC does not have a general transcriptional repression effect as 5 mC does (Wu et al., 2011; Kang et al., 2015; Li et al., 2016).

During development, genome-wide methylation changes occur before the differentiating cells acquire the adult-type methylation profiles. Similarly, induction of pluripotent stem cells or differentiation of stem cells in their physiological niches is accompanied by general DNA methylation changes. It is a dynamic and lifelong feature of DNA methylation. Local methylation changes happen in response to environmental factors, such as hormonal and metabolic effects or even early childhood stress (see later). The affected genes are silenced for long periods potentially through several generations. Toxic molecules, infections and hereditary or acquired diseases can have similar effects (Yamagata et al., 2012a). However, these alterations can undergo rapid reversion under appropriate environmental conditions. Due to the different environmental factors acting on each individual separately, DNA methylation changes throughout aging generate increasing methylation pattern differences between monozygotic twins, which leads to progressively appearing phenotypical variability (Fraga et al., 2005).

Histone Modifications

Histones are small, globular proteins, which are highly conserved in all eukaryotes. As mentioned earlier, the core histones building up the nucleosomes have N-terminal protruding tails, which are particularly prone to undergo posttranslational modifications (Allfrey et al., 1964). These epigenetic modifications include acetylation, methylation and phosphorylation. They play an important role in different nuclear processes, such as replication, DNA repair, transcription, and chromatin structure stabilization (Kouzarides, 2007; Bannister and Kouzarides, 2011). Although it was reported several decades ago that histones might undergo covalent modifications, their intensive investigation started only in the late 1990s.

The initial studies identified the lysine residues as targets of acetylation and revealed that they are essentially located on the H3 and H4 histone N-terminal tails. These modifications have rapid turnover (Zee et al., 2010). The reactions are catalysed by a high number of histone acetyl-transferases (HAT) and histone deacetylases (HDAC and Sirt) (Kuo and Allis, 1998; Legube and Trouche, 2003). Some complexes with HAT activity (e.g., p300 and CBP–CREB-binding protein) recruit also transcription factors and RNA PolII (Sakamoto et al., 2004). Thus, not surprisingly, lysine acetylation marks transcriptionally active euchromatic gene regulatory regions. For example, H3K27 (lysine 27 of histone H3) acetylation identifies active regulatory elements and separates active and inactive enhancers (ENCODE Project Consortium, 2012). H4 acetylation also indicates active chromatin territories; the acetylation of the two histones often occurs simultaneously. Acetylation profiles are inherited during DNA replication. The molecular mechanisms are still unclear and there might be several. According to one of them, histones with their epigenetic modifications are randomly distributed between the two new DNA molecules while the new chromatin is forming. Than newly synthesized histones are deposited and they rapidly acquire similar modifications to the old ones (Budhavarapu et al., 2013).

After understanding the role of histone acetylation, studies on histone methylation begun. Histone methylation has a much more complex pattern than acetylation since both arginines and lysines can be modified. Furthermore, arginines can be mono- or di-methylated, while lysines can be mono-, di-, or tri-methylated (Bannister and Kouzarides, 2011; Jahan and Davie, 2015; Greenblatt et al., 2016). Dozens of enzymes catalyse these post-translational modifications and their removal. The different enzymes are very selective and are capable to catalyse only one or two specific reactions (such as mono- and di-methylation of a specific lysine but not tri-methylation). Different histone methylations are specifically associated with various gene or chromatin regions. For example, H3K4me3 (trimethylation of lysine 4 of histone H3) marks transcription start sites (TSS) and it is characteristic of active promoters in the euchromatin (Santos-Rosa and Caldas, 2005). In addition, H3K36me3 is associated with actively transcribed gene bodies (Edmunds et al., 2008), while H3K27me3 typically occurs in transcriptionally repressed, heterochromatic regions (Bracken et al., 2006).

Networks of Epigenetic Modifications

In the different chromatin regions, complex patterns of histone and DNA modifications co-occur, which led to the formulation of the “histone code” hypothesis (Jenuwein and Allis, 2001). According to this hypothesis distinct patterns of chromatin modifications at any given genomic region would have the same meaning. These patterns are read by specific proteins, which then execute their local roles accordingly. However, it turned out that the histone code is most probably highly redundant.

Still, the epigenetic modifications are recognized by chromatin binding proteins. These proteins then often antagonize or promote the removal or catalysis of other covalent marks leading to the formation of complex patterns (Bannister and Kouzarides, 2011). Other proteins sensitive to the epigenetic pattern play important role in nuclear processes (e.g., γH2AX, a phosphorylated histone, participating in DNA repair). The proteins sensing and modifying the epigenetic marks are called readers, writers or erasers. Certainly, these proteins bind the chromatin with different affinity and therefore, they reside there for a shorter or longer time period depending on the local environmental context. The networks of modifications are therefore stochastically self-assembling and disassembling with various probability depending on the environmental conditions (Jeltsch and Jurkowska, 2014). This rapidly changing feature of the network makes the system particularly efficient in reacting to environmental stress (e.g., starvation, oxidative stress, or viral infection (Yamagata et al., 2012a) and modulating gene expression states in order to maintain the cellular homeostasis.

Therefore, it is not surprising that enzymes catalysing the addition or removal of covalent post-translational modifications of histones and DNA are closely linked to intermediary metabolism (Figure 2). To acetylate lysine residues, HAT enzymes use acetyl-CoA, a key molecule in carbohydrate and fat metabolism. Class III histone deacetylases (Sirt) need NAD+ for their activity (Vaquero et al., 2007). High level of energy intake leads to hyperacetylated histones, while low energy intake favors histone hypoacetylation. Methylation of histones and DNA needs S-adenosyl-methionin (SAM) as a cofactor, the methyl donor in biochemical reactions. The reaction also needs folate to regenerate SAM. The TeT (Ten-eleven Translocase) enzymes are oxygenases, which catalyse the demethylation reaction of DNA through the formation of 5-hydroxy-methyl cytosine. Jumonji family of histone demethylases have a similar reaction mechanism (Chen et al., 2006). Both TeT and Jumonji enzymes use α-ketoglutarate as a co-factor, which is a key metabolite of the citric-acid cycle. Enzymes catalysing the formation of α-ketoglutarate are different isoforms of isocitrate dehydrogenase (IDH). Some of the IDHs are mitochondrial, others are cytosolic, and they depend on NADP and NAD, respectively. SNPs of these enzymes are associated with TS (see later), while their mutations show frequent occurrence in tumors (gliomas, AML) (Dang et al., 2010). The gain-of-function mutant enzymes catalyse the formation of 2-hydroxyglutarate (2-HG), a potent inhibitor of demethylases. Altogether these data suggest a tight link between environmental factors and epigenetic modifications.

FIGURE 2

Figure 2. The relationship between epigenetic modifications and intermediary metabolism. Glycolysis, lipid metabolism, citric acid cycle, amino acid metabolism, and folate/SAM cycles are tightly linked to epigenetic modifications (shown in the middle), since their products and cofactors (shown in red) are substrates of enzymes catalyzing the epigenetic modifications. Acetyl-coenzyme A and NAD contribute to histone acetylation and deacetylation, respectively. Methyl groups and alfa-ketoglutarate participate in the methylation and demethylation of both histones and DNA. NAD: nicotinamide adenine dinucleotide, THF: tetrahydrofolate, SAH: S-adenosylhomocysteine, SAM: S-adenosylmethionine.

Non-coding RNA(ncRNA)

Evidence from the Encyclopedia of DNA Elements (ENCODE) suggests that at least 80% of our genome is transcribed. The human genome encodes for less than 3% of protein-coding transcripts and consists primarily of the non-coding RNAs (ncRNAs), which was previously regarded as “junk” DNA (https://www.encodeproject.org/). Non-coding RNAs play a role in gene expression regulation (see below). For example they are implicated in the regulation of genes coding for enzymes catalyzing epigenetic modifications. Furthermore, non-coding RNAs are also involved in the regulation of chromosomal domains in tight interaction with epigenetic covalent modifications (e.g., X-chromosome inactivation).

An arbitrary threshold of 200 nucleotides of transcript length was drawn to classify two groups of ncRNAs into small or long ncRNAs. Small ncRNAs include microRNAs (miRNA), transfer RNAs (tRNA) and small nucleolar RNAs (snRNA). A microRNA (miRNA) is a small non-coding RNA molecule found in plants, animals and some viruses, which functions in transcriptional and post-transcriptional regulation of gene expression.

The biogenesis of microRNA involves two different cleavage steps by protein complexes taking place in the nucleus and in the cytoplasm leading to the development of an ~22 nucleotide long mature single-stranded miRNA (Figure 3; Filipowicz et al., 2008). First, the gene coding for miRNAs is transcribed by RNA Polimerase II resulting in a long precursor pri-miRNA characterized by a hairpin or fold-back structure with an imperfectly base-paired stem and a terminal loop (Miyoshi et al., 2010). Then, the hairpin structure of the pri-miRNA is cleaved and the ~70 nucleotide long precursor called pre-miRNA is released (Miyoshi et al., 2010). Subsequently, the pre-miRNA is exported into the cytoplasm (Kim et al., 2009). Once in the cytoplasm, the pre-miRNA is processed by RNase III Dicer generating a ~22 nucleotide miRNA-miRNA*duplex. One strand of the duplex is loaded into a large multi-protein miRNA ribonucleoprotein complex (RISC complex), while the other strand, “passenger,” is degraded. Once incorporated into RISC, the miRNA guides the complex to its messenger RNA targets by base-pairing interactions. The binding to the miRNA target relies on the seed sequence, a 6–8 nucleotide domain located at the 5′ end of the miRNA. Based on the complementarity of the seed sequence and the target mRNA sequence located in the 3′UTRof the transcript, miRNAs down-regulate gene expression either through translational repression or mRNA degradation (Hutvágner and Zamore, 2002).

FIGURE 3

Figure 3. miRNA biogenesis. MicroRNA (miRNA) genes are transcribed as primary miRNAs (pri-miRNAs) by RNA polymerase II (Pol II) in the nucleus. The long pri-miRNAs are cleaved by Microprocessor, which includes DROSHA and DiGeorge syndrome critical region 8 (DGCR8), to produce precursor miRNAs (pre-miRNAs), which are then exported to the cytoplasm by Exportin 5 and further processed by the DICER/TRPB complex to produce an miRNA duplex. One strand of the mature miRNA (the guide strand) is loaded into the miRNA-induced silencing complex (RISC) mediating gene suppression by targeted mRNA degradation or translational repression.

Longer RNAs include ribosomal RNAs (rRNA), natural antisense transcripts and other long non-coding RNAs (lncRNAs). Currently, about 2500 miRNAs and 50,000 lncRNAs have been annotated in the human genome, besides the ~19,000 protein coding genes. As opposed to miRNAs, only a few lncRNAs show evolutionary conservation of the primary sequence, but most of them show tissue and cell type-specific expression, indicating that their expression must be tightly controlled. LncRNAs exhibit a diversity of molecular functions: they can act as transcriptional activators or repressors, as scaffolds for protein-protein interactions or as molecular decoys.

Genetic polymorphisms can involve both ncRNA sequences throughout the genome, as well as sequences in their target genes, which can affect ncRNA-mediated regulation. Polymorphisms in ncRNA genes can influence both their level of expression and the ncRNA function, thus resulting in differential regulation of their target genes.

Genetic variation may affect miRNA-mediated regulation in different ways. Altered miRNA transcription can be a result of polymorphisms in (a) miRNA promoters, (b) splice sites of the host gene where miRNAs reside, since many miRNAs are intronic, or (c) of polycistronic miRNA clusters (Calore et al., 2015). Mutations in the sequence of the transcribed miRNA may change its binding affinity to the biogenesis enzyme complexes, which changes their processing and may lead to misregulation of target genes. IsomiRs are variants within the processed miRNA sequence, which can affect specificity to their target, on the other hand target messenger RNA 3′UTR variants may either destroy existing miRNA seed regions or create new recognition sites. The presence of a SNP in the 3′UTR can theoretically result in three alterations in miRNA related regulation: (a) it can partially or completely disrupt the miRNA binding site, thus resulting in higher expression of the target gene, or (b) more rarely a SNP could either enhance the binding of a miRNA to the 3′UTR region through improvement of the original recognition site or (c) it may create a novel binding site for another miRNA. The latter will only affect the expression of the target gene if it coincides with the expression of the new miRNA spatially and temporally.

Although lncRNAs possess a much lower degree of conservation than miRNAs, their genetic polymorphisms may still be functional, i.e., SNPs within the lncRNA loci can change their expression or influence their downstream target genes. Altering lncRNA architecture may influence its ability to interact with proteins or other RNAs. In recent years, genetic variations in lncRNAs were implicated in several human diseases.

Epigenetics of Neuropsychiatric Disorders

Epimutations are epigenetic alterations, which have been linked to several diseases. These alterations can be classified as primary or secondary based on their origin (Oey and Whitelaw, 2014). Primary epimutations seem to be due to only environmental stress factors. It is often difficult to understand the pathomechanism of the resulting diseases or traits, which are generally less severe than in the case of secondary epimutations. Secondary epimutations are due to an initial genomic mutation. Most of these mutations target readers, writers and erasers of the epigenetic system and can lead to important changes of the global epigenetic profile (Lopez-Atalaya et al., 2014). In the following section we will describe some examples of epimutations leading to neuropsychiatric disorders (Plazas-Mayorca and Vrana, 2011).

Rubinstein-Taybi Syndrome

Rubinstein-Taybi syndrome has various clinical signs including moderate to severe learning difficulties (for a recent review see Lopez-Atalaya et al., 2014). The hereditary disease is autosomal dominant although very few documented cases of transmission exist. Most of the patients have de novo mutations. The development of the syndrome can be attributed to mutations in the CREB-binding protein (CBP) or more rarely in the highly homologous p300 protein encoding genes. Both proteins are transcriptional co-activators, have HAT activity and they bind the acetylated histones via their bromodomain. Thus, they provide platforms for other proteins (transcription factors and RNA PolII) necessary for transcription initiation.

Gliomas

Gliomas (Vigneswaran et al., 2015) are tumors arising from glial cells in the central nervous system (CNS). During tumor progression patients experience psychiatric, cognitive, and neurologic symptoms (Boele et al., 2015). Primary low-grade gliomas are typically diagnosed in the 40 s, and after treatment, these slow-growing tumors have a tendency to reappear and progress in grade to become grade III gliomas or glioblastomas (grade IV). In contrast to the secondary high-grade glioblastomas, the primary high-grade glioblastomas are diagnosed later and have very poor prognosis. Although histologically identical, the primary and secondary glioblastomas have different molecular characteristics.

Both primary low-grade and secondary high-grade gliomas are characterized by IDH1/2 mutation (see above). Approximately 90% of the mutations occur in the IDH1 gene and almost all of them are the R132H variant (Vigneswaran et al., 2015). As mentioned earlier, this enzyme variant catalyzes the formation of an oncometabolite (2-HG), which inhibits DNA and histone demethylation leading to general (secondary) alteration of the epigenetic profile (Cohen et al., 2013). Nevertheless, gliomas with mutated IDH1/2 gene have better prognosis than the others (Andronesi et al., 2013).

Rett Syndrome

Rett syndrome is a disease affecting only girls (1:10,000–15,000) (Chahrour and Zoghbi, 2007; Katz et al., 2012) and it is characterized by early onset (18 months) and variable neurological symptoms. Severe mental retardation and motor impairment such as ataxia, apraxia, and tremor (Chahrour and Zoghbi, 2007) accompanied by seizures and gastrointestinal symptoms are frequently present (Katz et al., 2012). Rett syndrome is an example of a severe disorder due to the mutation of the MeCP2 gene encoding an epigenetic “reader” i.e., a regulatory protein recognizing an epigenetic mark (Amir et al., 1999). Therefore, although Rett syndrome is considered to be a typical epigenetic disease, no major epigenetic alterations can be observed in the patients. Loss-of-function mutations of the gene coding for the transcriptional repressor MeCP2 are responsible for the development of the disease (Amir et al., 1999). MeCP2 is a member of the methyl-CpG binding domain (MBD) protein family. Upon binding of methylated DNA, MeCP2 recruits transcriptional repressor complexes and HDACs. Interestingly, MeCP2 has recently been found to bind and repress long genes implicated in neuronal differentiation and modulation of neuronal functions (Gabel et al., 2015).

Autism Spectrum Disorder (ASD)

Rett syndrome is also considered to be a rare form of autism spectrum disorder (ASD) (Mbadiwe and Millis, 2013). ASD is clinically characterized by social communication deficits and repetitive behavior, which appears as early as 2 years of age and causes clinically significant impairment. ASD has high heritability rates suggesting a substantial genetic background. It has a polygenic origin with hundreds of susceptibility genes. Most of them are common variants with small effects, while some rare de novo variants with large effects also exist (Loke et al., 2015). Apart from MeCP2, FMR1, and OXTR also have profound effects and they were repeatedly reported in relation with ASD (Mbadiwe and Millis, 2013). Both of them are also linked to epigenetic alterations, which is clearly secondary in case of FMR1.

Fragile X syndrome

Approximately half of the patients with Fragile X syndrome meet the criteria of autism and it is a relatively frequent cause of ASD (5% of all monogenic cases). The FMR1 gene encodes FMRP, a polyribosome associated protein playing an important role in protein translation (Penagarikano et al., 2007). The absence of the protein leads to perturbed neuronal development and intellectual disability (Contractor et al., 2015). Fragile X syndrome is caused by a CGG trinucleotide expansion in the regulatory region of the FMR1 gene located on the chromosome X. This repeat expansion (>200) leads to the attraction of DNA methylation and the loss of expression of the gene (Oberlé et al., 1991; Penagarikano et al., 2007).

The oxytocin receptor (OXTR) gene is also a candidate gene for ASD (Loke et al., 2015). The neurotransmitter and hormone oxytocin was found to play a role in anxiety, aggressive behavior, and other neural functions. Several observations indicate that OXTR plays a role in the development of ASD. For instance, four SNPs in the gene were suggested to be associated with ASD. Furthermore, several studies reported higher DNA methylation level in patients than in controls in the promoter region of the gene (Gregory et al., 2009; Jack et al., 2012; Ziegler et al., 2015). Very important methylation increase (20–40%) was observed both in temporal cortex and peripheral blood. This temporal hypermethylation was also correlated with lower OXTR mRNA levels in autists (Gregory et al., 2009). Thus, it is not surprising that OXTR methylation has been associated with anxiety disorder and other traits characterizing ASD.

A recent genome-wide analysis of DNA methylation studied a sample of 50 monozygotic (MZ) twin pairs including twins discordant, as well as concordant for ASD. Within-twin and between-group analyses identified a number of differentially methylated regions associated with ASD. In addition, the authors reported significant correlations between DNA methylation and quantitatively measured autistic trait scores across the cohort implicating a role for altered DNA methylation in autism (Wong et al., 2014).

Finally, strong evidence shows that ASD has primary epigenetic origin, as well (Tordjman et al., 2014). Children with in utero exposure to the HDAC inhibitor valproic acid (an anticonvulsive and mood stabilizer drug) were found in several studies (Moore et al., 2000; Bromley et al., 2013; Christensen et al., 2013) to have a significantly increased risk to develop autism relative to those who were not treated. Other environmental factors during pregnancy can be considered as risk factors for ASD, such as viral infection (e.g., rubella) (Ornoy et al., 2015) and the dietary folic acid supplementation, which is regarded as controversial (Yang et al., 1989; Mbadiwe and Millis, 2013). Finally, several studies indicate that prenatal maternal stress is also a risk factor for developing ASD (Kinney et al., 2008 and references therein).

A considerable number of studies indicate that early life adversities (ELA) (e.g., childhood abuse or even prenatal and/or maternal stress) are severe risk factors for the development of psychiatric disorders such as major depression, suicidal behavior, etc. (Hoffmann and Spengler, 2014; Palma-Gudiel et al., 2015; Cattaneo and Riva, 2016). This seems to be due to a difficulty in coping with stress in these patients. During stress reactions the hypothalamus-pituitary-adrenal axis (HPA) is activated and glucocorticoids (cortisol in humans and corticosterone) are released, which in turn activate the pathways regulated by glucocorticoid and mineralocorticoid receptors. The axis is inhibited by the feedback activation of the glucocorticoid receptor (GR) in the hippocampus (Palma-Gudiel et al., 2015; Cattaneo and Riva, 2016).

The molecular mechanisms of the development of ELA-related alterations have been deciphered in a rat model system (Weaver et al., 2004). Weaver and colleagues have compared pups from “good” nursing and “bad” nursing females, two maternal behaviors generally occurring in rats (Liu et al., 1997; Caldji et al., 1998). Offsprings of “good moms” were less fearful and had a lower stress response in their adulthood than those of “bad moms.” Weaver and colleagues observed higher DNA methylation levels in the hippocampus from early childhood (after postnatal day 1) until at least 3 months of age in the GR promoter in the offsprings of “bad moms” relative to those of “good moms.” This difference concentrated at a certain region of the promoter and more precisely at a single CpG, located at the binding site of transcription factor NGFI-A regulating GR expression. They also observed the decreased binding of NGFI-A and hypoacetylation of histones in the methylated region of the promoter. These findings were accompanied by lower GR expression. Interestingly, but not surprisingly for an epigenetic mark, both the molecular and the behavioral phenotypes were reversible. Treatment with Trichostatin A, an HDAC inhibitor, reversed histone hypoacetylation, increased NGFI-A binding, GR expression, and decreased DNA methylation to some extent. Similarly, changing the environment had a similar effect as shown by cross-fostering, demonstrating that this phenotype is not determined by the genetic background and should be considered as having a primary epigenetic origin.

Based on these observations, several studies in animal models confirmed these findings (McGowan et al., 2011). In human cohorts a very small, but systematic methylation increase was reported from the same region of the GR promoter in individuals undergoing stressful events (Palma-Gudiel et al., 2015).

Non-coding RNAs are also associated with a wide range of neurodevelopmental, neurodegenerative, and psychiatric diseases both in humans and in animal models (Lin et al., 2011; Johnson, 2012; Talkowski et al., 2012; Ziats and Rennert, 2012; Nishimoto et al., 2013; Petazzi et al., 2013; Barry et al., 2014). As a recent example, genetic variants of the long non-coding RNA MIAT were found to contribute to the risk of paranoid schizophrenia in a Han Chinese population (Rao et al., 2015). The authors performed a two-stage association analysis on 8 tagging SNPs covering the whole MIAT locus in two independent Han Chinese schizophrenia case–control cohorts. The discovery sample with over 1000 cases and 1000 controls yielded a significant increase of the minor T-allele of rs1894720 in patients and this association was confirmed in the replication cohort of a similar size.

MicroRNAs are also implicated in several neuropsychiatric disorders, such as schizophrenia and autism (Beveridge and Cairns, 2012; Mellios and Sur, 2012), neurodegenerative disorders like Alzheimer's and Parkinson's disease (Salta and De Strooper, 2012; Abe and Bonini, 2013; Tan et al., 2013), but also in other neurodevelopmental disorders such as Fragile X syndrome and Rett syndrome (Urdinguio et al., 2010; Wu et al., 2010; Im and Kenny, 2012; Sellier et al., 2013).

Epigenetics of Tourette Syndrome

As described earlier, the studies investigating TS mainly focused on the genetic background of the disease. Only few studies have been performed to date to investigate the role of epigenetic factors and non-coding RNA in the development of TS. One of these identified a nucleotide variant (var321) in the 3′ UTR of the SLITRK1 gene leading to its stronger repression by miR-189. This variant has been investigated in several studies in Tourette patients and reported to be rare (Abelson et al., 2005). However, the role of this variant in TS pathogenesis is questionable due to its very low frequency (Keen-Kim et al., 2006). The unique study reported to date on the role of microRNAs in Tourette Syndrome profiled the expression of 754 miRNAs in the sera of six TS patients and three unaffected controls (Rizzo et al., 2015). The study found that miR-429, which is involved in midbrain and hindbrain differentiation and synaptic transmission was significantly underexpressed in TS patients. Measurement of circulating miR-429 may in the future be useful as a molecular biomarker to aid TS diagnosis.

Two association studies on DNA methylation related to TS have been conducted so far. The first showed no methylation alteration in TS patients relative to controls in a region on chromosome 8 in KCNK9 and TRAPPC9 genes (Sánchez Delgado et al., 2014). The regions were identified recently by genome-wide screens and by mapping mutations in single families. The other study was the first Epigenome-Wide Association Study (EWAS) investigating DNA methylation differences between hundreds of controls and patients from the Netherlands Twin Registry with tic phenotype (Zilhão et al., 2015). Very small methylation differences were observed, however, an enrichment of differentially methylated neural genes previously linked to neuropsychiatric disorders or with brain specific function was found among the top hits.

Finally, a recent promising study investigated GWAS results on gene sets. The association of TS and a gene set related to carbohydrate metabolism and more particularly a group of 33 genes involved in “astrocyte-neuron metabolic coupling,” glycolysis and Krebs-cycle was demonstrated (de Leeuw et al., 2015). This is interesting because the genes identified include IDH2 (see above) and Malic enzyme 1, which is regulated by glucose level. Since these TS associated metabolic genes are known to play a role in epigenetic modifications, this suggeststhat the disorder is potentially characterized by altered neural epigenetic patterns.

Outlook

In the present review, we have shown that TS is a neuropsychiatric disorder with significant heritability. However, while very few rare genetic variants with large effects were described, it is plausible to assume that hundreds or even more frequent variants with small effects underlie the genetic susceptibility of the disorder. Furthermore, a considerable number of observations indicate that environmental factors also play a crucial role in the development of TS. As introduced above, environmental factors act via epigenetic modifications, including heritable covalent modifications of the chromatin and regulatory non-coding RNAs. In order to better understand the pathomechanism of TS, we propose here that more studies should be performed focusing on the role of epigenetics.

What questions should be asked? We think that since environmental (risk) factors are implicated in TS, these should be studied with particular interest. For instance, discordant monozygotic twins are very good candidates for finding epigenetic modifications implicated in the development of the disease, since they are genetically identical. Similarly, patients with known prenatal or perinatal antecedents probably also have important epigenetic grounding in the development of TS. We also consider that patients who developed TS due to streptococcus infection should also show epigenetic alteration relative to controls. Finally, since the disease is characterized by waxing and waning, kinetic analysis of epigenetic changes might also reflect important aspects of TS.

Animal models of TS or tic phenotype should also be studied for epigenetic alterations. While studies on human cohorts might be more descriptive, analyses of animal models could explain more directly the underlying molecular mechanisms.

How should these questions be addressed? Currently several techniques are available to study epigenetics. Genome-wide analyses are much more informative than investigations of single targets. These genome-wide approaches have recently become much less expensive and therefore affordable for most of the laboratories or consortia. Expression level and GWAS analyses can be performed on ncRNA, while chromatin immunoprecipitation (ChIP) followed by next generation sequencing (NGS) is useful for the investigation of histone modifications (Furey, 2012). Along with the ongoing technological advances in NGS, identifying genetic variation affecting ncRNA function associated with neuropsychiatric disease is likely to grow rapidly.

Finally, several genome-wide techniques exist for the study of DNA hydroxymethylation and DNA methylation. We propose to study DNA methylation rather than histone modifications, because it is more stable, than the latter one, although not as static as initially thought (Yamagata et al., 2012b). Furthermore, ChIP-like antibody-based techniques are much less quantitative than the chemical transformation (bisulfite conversion-based) techniques developed for DNA methylation. Therefore either array-based hybridization approaches after bisulfite conversion [Illumina 450 k (Bibikova et al., 2011) or the recent 850 k bead chips] or bisulfite conversion coupled to next generation sequencing could yield useful information (Kulis et al., 2012).

What cautions should be taken? First, we consider that the most important issues are sample size and tissue of origin. TS is a psychiatric disease mainly due to alterations in neuronal and potentially in glial cells. These cells in humans are rarely available for research, while blood cells and buccal or nasopharingeal epithelial cells have somewhat different epigenetic signatures. Although some of these cells can have similar epigenetic profiles in some genomic regions to the most important brain regions in TS (e.g., striatum), the results should always be interpreted cautiously (Hannon et al., 2015). In order to avoid such problems, brain samples should be investigated when possible, however that has obvious limitations in humans. Alternatively, as mentioned before discordant monozygotic twins can also be studied, finally, patients with TS presumably due to PANDAS origin are good candidates for blood sample analysis. It is hard to determine the ideal sample size. However, it is clear that already small sample sizes can be informative if repeated experiments give similar results and the samples are well-selected.

Second, new approaches available and proposed lead to the generation of “big data.” Their correct analysis is necessary and requires the intensive collaboration of the biomedical and bioinformatitian scientists with profound knowledge of statistics.

Finally, when interpreting data, attention and caution should be exercised. Although biologically meaningful cutoffs for methylation differences between patients and controls is hard to determine, the value of very small but statistically significant changes is questionable. Furthermore, statistically significant hits showing association with the phenotype does not neccessarily mean causality. From a single descriptive experiment causality cannot be concluded, but further experiments should be performed.

In conclusion, we consider that introducing epigenetic studies in TS research has great potential. These investigations based on the previous results, the animal models and the twin and patients registries already existing will certainly identify new molecular mechanisms and hits playing an important role in the development of the disease. The discovery of the molecular details of ncRNA and epigenetic modification mediated regulation of gene expression, proteins and pathways is likely to provide novel insights into the pathogenesis of neuropsychiatric disease including Tourette Syndrome. This will also open new avenues for genetic diagnosis, as well as targeted and personalized therapeutic approaches. These studies will also strengthen the importance of some already known and suspected hits and altogether this will lead to the better understanding of TS and the opening of new avenues for the development of treatments.

Author Contributions

LP, BV, TA, and CB were all involved in building up the concept of the paper, literature research, and writing of the manuscript.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

This project has been financed by FP7-People-2012-ITN, project: TS-Eurotrain, grant number 316978. CB was supported by the Merit-prize fellowship of Semmelweis University and the Bolyai János research fellowship of the Hungarian Academy of Sciences BO/00987/16/5.

References

Abe, M., and Bonini, N. M. (2013). MicroRNAs and neurodegeneration: role and impact. Trends Cell Biol. 23, 30–36. doi: 10.1016/j.tcb.2012.08.013

Abelson, J. F., Kwan, K. Y., O'Roak, B. J., Baek, D. Y., Stillman, A. A., Morgan, T. M., et al. (2005). Sequence variants in SLITRK1 are associated with Tourette's syndrome. Science 310, 317–320. doi: 10.1126/science.1116502

Allfrey, V. G., Faulkner, R., and Mirsky, A. E. (1964). Acetylation and methylation of histones and their possible role in the regulation of RNA synthesis. Proc. Natl. Acad. Sci. U.S.A. 51, 786–794.

Amir, R. E., Van den Veyver, I. B., Wan, M., Tran, C. Q., Francke, U., and Zoghbi, H. Y. (1999). Rett syndrome is caused by mutations in X-linked MECP2, encoding methyl-CpG-binding protein 2. Nat. Genet. 23, 185–188. doi: 10.1038/13810

Andronesi, O. C., Rapalino, O., Gerstner, E., Chi, A., Batchelor, T. T., Cahill, D. P., et al. (2013). Detection of oncogenic IDH1 mutations using magnetic resonance spectroscopy of 2-hydroxyglutarate. J. Clin. Invest. 123, 3659–3663. doi: 10.1172/JCI67229

Antonarakis, S. E., Krawczak, M., and Cooper, D. N. (2000). Disease-causing mutations in the human genome. Eur. J. Pediatr. 159(Suppl. S173–S178). doi: 10.1007/PL00014395

Arányi, T., Faucheux, B. A., Khalfallah, O., Vodjdani, G., Biguet, N. F., Mallet, J., et al. (2005). The tissue-specific methylation of the human tyrosine hydroxylase gene reveals new regulatory elements in the first exon. J. Neurochem. 94, 129–139. doi: 10.1111/j.1471-4159.2005.03173.x

Aruga, J., and Mikoshiba, K. (2003). Identification and characterization of Slitrk, a novel neuronal transmembrane protein family controlling neurite outgrowth. Mol. Cell. Neurosci. 24, 117–129. doi: 10.1016/S1044-7431(03)00129-5

Asghari, V., Sanyal, S., Buchwaldt, S., Paterson, A., Jovanovic, V., and Van Tol, H. H. M. (2002). Modulation of intracellular cyclic AMP levels by different human dopamine D4 receptor variants. J. Neurochem. 65, 1157–1165. doi: 10.1046/j.1471-4159.1995.65031157.x

Baba, Y., Iida, A., and Watanabe, S. (2011). Sall3 plays essential roles in horizontal cell maturation through regulation of neurofilament expression levels. Biochimie 93, 1037–1046. doi: 10.1016/j.biochi.2011.02.016

Bannister, A. J., and Kouzarides, T. (2011). Regulation of chromatin by histone modifications. Cell Res. 21, 381–395. doi: 10.1038/cr.2011.22

Baron-Cohen, S., Mortimore, C., Moriarty, J., Izaguirre, J., and Robertson, M. (1999). The prevalence of Gilles de la Tourette's syndrome in children and adolescents with autism. J. Child Psychol. Psychiatry 40, 213–218.

Barry, G., Briggs, J. A., Vanichkina, D. P., Poth, E. M., Beveridge, N. J., Ratnu, V. S., et al. (2014). The long non-coding RNA Gomafu is acutely regulated in response to neuronal activation and involved in schizophrenia-associated alternative splicing. Mol. Psychiatry 19, 486–494. doi: 10.1038/mp.2013.45

Bertelsen, B., Melchior, L., Jensen, L. R., Groth, C., Glenthøj, B., Rizzo, R., et al. (2014). Intragenic deletions affecting two alternative transcripts of the IMMP2L gene in patients with Tourette syndrome. Eur. J. Hum. Genet. 22, 1283–1289. doi: 10.1038/ejhg.2014.24

Bertelsen, B., Melchior, L., Jensen, L. R., Groth, C., Nazaryan, L., Debes, N. M., et al. (2015). A t(3;9)(q25.1;q34.3) translocation leading to OLFM1 fusion transcripts in Gilles de la Tourette syndrome, OCD and ADHD. Psychiatry Res. 225, 268–275. doi: 10.1016/j.psychres.2014.12.028

Bertelsen, B., Stefánsson, H., Riff Jensen, L., Melchior, L., Mol Debes, N., Groth, C., et al. (2016). Association of AADAC Deletion and Gilles de la Tourette Syndrome in a Large European Cohort. Biol. Psychiatry 79, 383–391. doi: 10.1016/j.biopsych.2015.08.027

Beveridge, N. J., and Cairns, M. J. (2012). MicroRNA dysregulation in schizophrenia. Neurobiol. Dis. 46, 263–271. doi: 10.1016/j.nbd.2011.12.029

Bibikova, M., Barnes, B., Tsan, C., Ho, V., Klotzle, B., Le, J. M., et al. (2011). High density DNA methylation array with single CpG site resolution. Genomics 98, 288–295. doi: 10.1016/j.ygeno.2011.07.007

Bird, A. P. (1993). Functions for DNA methylation in vertebrates. Cold Spring Harb. Symp. Quant. Biol. 58, 281–285.

Bird, A., Taggart, M., Frommer, M., Miller, O. J., and Macleod, D. (1985). A fraction of the mouse genome that is derived from islands of nonmethylated, CpG-rich DNA. Cell 40, 91–99.

Blattler, A., and Farnham, P. J. (2013). Cross-talk between site-specific transcription factors and DNA methylation states. J. Biol. Chem. 288, 34287–34294. doi: 10.1074/jbc.R113.512517

Boele, F. W., Rooney, A. G., Grant, R., and Klein, M. (2015). Psychiatric symptoms in glioma patients: from diagnosis to management. Neuropsychiatr. Dis. Treat. 11, 1413–1420. doi: 10.2147/NDT.S65874

Boghosian-Sell, L., Comings, D. E., and Overhauser, J. (1996). Tourette syndrome in a pedigree with a 7;18 translocation: identification of a YAC spanning the translocation breakpoint at 18q22.3. Am. J. Hum. Genet. 59, 999–1005.

Bontempi, S., Fiorentini, C., Busi, C., Guerra, N., Spano, P., and Missale, C. (2007). Identification and characterization of two nuclear factor-kappaB sites in the regulatory region of the dopamine D2 receptor. Endocrinology 148, 2563–2570. doi: 10.1210/en.2006-1618

Borroto-Escuela, D. O., Van Craenenbroeck, K., Romero-Fernandez, W., Guidolin, D., Woods, A. S., Rivera, A., et al. (2011). Dopamine D2 and D4 receptor heteromerization and its allosteric receptor-receptor interactions. Biochem. Biophys. Res. Commun. 404, 928–934. doi: 10.1016/j.bbrc.2010.12.083

Bracken, A. P., Dietrich, N., Pasini, D., Hansen, K. H., and Helin, K. (2006). Genome-wide mapping of Polycomb target genes unravels their roles in cell fate transitions. Genes Dev. 20, 1123–1136. doi: 10.1101/gad.381706

Bromley, R. L., Mawer, G. E., Briggs, M., Cheyne, C., Clayton-Smith, J., García-Fiñana, M., et al. (2013). The prevalence of neurodevelopmental disorders in children prenatally exposed to antiepileptic drugs. J. Neurol. Neurosurg. Psychiatry 84, 637–643. doi: 10.1136/jnnp-2012-304270

Budhavarapu, V. N., Chavez, M., and Tyler, J. K. (2013). How is epigenetic information maintained through DNA replication? Epigenetics Chromatin 6:32. doi: 10.1186/1756-8935-6-32

Burd, L., Freeman, R. D., Klug, M. G., and Kerbeshian, J. (2005). Tourette Syndrome and learning disabilities. BMC Pediatr. 5:34. doi: 10.1186/1471-2431-5-34

Caldji, C., Tannenbaum, B., Sharma, S., Francis, D., Plotsky, P. M., and Meaney, M. J. (1998). Maternal care during infancy regulates the development of neural systems mediating the expression of fearfulness in the rat. Proc. Natl. Acad. Sci. U.S.A. 95, 5335–5340.

Calore, M., De Windt, L. J., and Rampazzo, A. (2015). Genetics meets epigenetics: genetic variants that modulate noncoding RNA in cardiovascular diseases. J. Mol. Cell. Cardiol. 89, 27–34. doi: 10.1016/j.yjmcc.2015.10.028

Castellan Baldan, L., Williams, K. A., Gallezot, J.-D., Pogorelov, V., Rapanelli, M., Crowley, M., et al. (2014). Histidine decarboxylase deficiency causes tourette syndrome: parallel findings in humans and mice. Neuron 81, 77–90. doi: 10.1016/j.neuron.2013.10.052

Cattaneo, A., and Riva, M. A. (2016). Stress-induced mechanisms in mental illness: a role for glucocorticoid signalling. J. Steroid Biochem. Mol. Biol. 160, 169–174. doi: 10.1016/j.jsbmb.2015.07.021

Cavanna, A. E., and Termine, C. (2012). Tourette syndrome. Adv. Exp. Med. Biol. 724, 375–383. doi: 10.1007/978-1-4614-0653-2_28

Chahrour, M., and Zoghbi, H. Y. (2007). The story of Rett syndrome: from clinic to neurobiology. Neuron 56, 422–437. doi: 10.1016/j.neuron.2007.10.001

Che, F., Zhang, Y., Wang, G., Heng, X., Liu, S., and Du, Y. (2015). The role of GRIN2B in Tourette syndrome: results from a transmission disequilibrium study. J. Affect. Disord. 187, 62–65. doi: 10.1016/j.jad.2015.07.036

Chen, Z., Zang, J., Whetstine, J., Hong, X., Davrazou, F., Kutateladze, T. G., et al. (2006). Structural insights into histone demethylation by JMJD2 family members. Cell 125, 691–702. doi: 10.1016/j.cell.2006.04.024

Chou, I.-C., Lin, H.-C., Wang, C.-H., Lin, W.-D., Lee, C.-C., Tsai, C.-H., et al. (2010). Polymorphisms of interleukin 1 gene IL1RN are associated with Tourette syndrome. Pediatr. Neurol. 42, 320–324. doi: 10.1016/j.pediatrneurol.2010.01.006

Christensen, J., Grønborg, T. K., Sørensen, M. J., Schendel, D., Parner, E. T., Pedersen, L. H., et al. (2013). Prenatal valproate exposure and risk of autism spectrum disorders and childhood autism. JAMA 309, 1696–1703. doi: 10.1001/jama.2013.2270

Clarke, R. A., Fang, Z. M., Diwan, A. D., and Gilbert, D. L. (2009). Tourette syndrome and klippel-feil anomaly in a child with chromosome 22q11 duplication. Case Rep. Med. 2009:361518. doi: 10.1155/2009/361518

Coffey, B. J., Biederman, J., Smoller, J. W., Geller, D. A., Sarin, P., Schwartz, S., et al. (2000). Anxiety disorders and tic severity in juveniles with Tourette's disorder. J. Am. Acad. Child Adolesc. Psychiatry 39, 562–568. doi: 10.1097/00004583-200005000-00009

Cohen, A. L., Holmen, S. L., and Colman, H. (2013). IDH1 and IDH2 mutations in gliomas. Curr. Neurol. Neurosci. Rep. 13:345. doi: 10.1007/s11910-013-0345-4

Comings, D. E., Comings, B. G., Muhleman, D., Dietz, G., Shahbahrami, B., Tast, D., et al. (1991). The dopamine D2 receptor locus as a modifying gene in neuropsychiatric disorders. JAMA 266, 1793–1800.

Comings, D. E., Gade, R., Muhleman, D., Chiu, C., Wu, S., To, M., et al. (1996a). Exon and intron variants in the human tryptophan 2,3-dioxygenase gene: potential association with Tourette syndrome, substance abuse and other disorders. Pharmacogenetics 6, 307–318.

Comings, D. E., Muhleman, D., Dietz, G., Dino, M., LeGro, R., and Gade, R. (1993). Association between Tourette's syndrome and homozygosity at the dopamine D3 receptor gene. Lancet 341, 906.

Comings, D. E., Wu, S., Chiu, C., Ring, R. H., Gade, R., Ahn, C., et al. (1996b). Polygenic inheritance of Tourette syndrome, stuttering, attention deficit hyperactivity, conduct, and oppositional defiant disorder: the additive and subtractive effect of the three dopaminergic genes–DRD2, D beta H, and DAT1. Am. J. Med. Genet. 67, 264–288. doi: 10.1002/(SICI)1096-8628(19960531)67:3<264::AID-AJMG4>3.0.CO;2-N

Contractor, A., Klyachko, V. A., and Portera-Cailliau, C. (2015). Altered neuronal and circuit excitability in Fragile X syndrome. Neuron 87, 699–715. doi: 10.1016/j.neuron.2015.06.017

Consortium, T. E. P. (2012). An integrated encyclopedia of DNA elements in the human genome. Nature 489, 57. doi: 10.1038/nature11247

Crane, J., Fagerness, J., Osiecki, L., Gunnell, B., Stewart, S. E., Pauls, D. L., et al. (2011). Family-based genetic association study of DLGAP3 in Tourette Syndrome. Am. J. Med. Genet. B. Neuropsychiatr. Genet. 156B, 108–114. doi: 10.1002/ajmg.b.31134

Dang, L., Jin, S., and Su, S. M. (2010). IDH mutations in glioma and acute myeloid leukemia. Trends Mol. Med. 16, 387–397. doi: 10.1016/j.molmed.2010.07.002

Davis, L. K., Yu, D., Keenan, C. L., Gamazon, E. R., Konkashbaev, A. I., Derks, E. M., et al. (2013). Partitioning the heritability of Tourette syndrome and obsessive compulsive disorder reveals differences in genetic architecture. PLoS Genet. 9:e1003864. doi: 10.1371/journal.pgen.1003864

de Leeuw, C., Goudriaan, A., Smit, A. B., Yu, D., Mathews, C. A., Scharf, J. M., et al. (2015). Involvement of astrocyte metabolic coupling in Tourette syndrome pathogenesis. Eur. J. Hum. Genet. 23, 1519–1522. doi: 10.1038/ejhg.2015.22

Deaton, A. M., and Bird, A. (2011). CpG islands and the regulation of transcription. Genes Dev. 25, 1010–1022. doi: 10.1101/gad.2037511

Dehning, S., Müller, N., Matz, J., Bender, A., Kerle, I., Benninghoff, J., et al. (2010). A genetic variant of HTR2C may play a role in the manifestation of Tourette syndrome. Psychiatr. Genet. 20, 35–38. doi: 10.1097/YPG.0b013e32833511ce

Díaz-Anzaldúa, A., Joober, R., Rivière, J.-B., Dion, Y., Lespérance, P., Chouinard, S., et al. (2004). Association between 7q31 markers and Tourette syndrome. Am. J. Med. Genet. A 127A, 17–20. doi: 10.1002/ajmg.a.20631

Edmunds, J. W., Mahadevan, L. C., and Clayton, A. L. (2008). Dynamic histone H3 methylation during gene induction: HYPB/Setd2 mediates all H3K36 trimethylation. EMBO J. 27, 406–420. doi: 10.1038/sj.emboj.7601967

Elia, J., Gai, X., Xie, H. M., Perin, J. C., Geiger, E., Glessner, J. T., et al. (2010). Rare structural variants found in attention-deficit hyperactivity disorder are preferentially associated with neurodevelopmental genes. Mol. Psychiatry 15, 637–646. doi: 10.1038/mp.2009.57

Ercan-Sencicek, A. G., Stillman, A. A., Ghosh, A. K., Bilguvar, K., O'Roak, B. J., Mason, C. E., et al. (2010). L-histidine decarboxylase and Tourette's syndrome. N. Engl. J. Med. 362, 1901–1908. doi: 10.1056/NEJMoa0907006

Filipowicz, W., Bhattacharyya, S. N., and Sonenberg, N. (2008). Mechanisms of post-transcriptional regulation by microRNAs: are the answers in sight? Nat. Rev. Genet. 9, 102–114. doi: 10.1038/nrg2290

Fiorentini, C., Guerra, N., Facchetti, M., Finardi, A., Tiberio, L., Schiaffonati, L., et al. (2002). Nerve growth factor regulates dopamine D(2) receptor expression in prolactinoma cell lines via p75(NGFR)-mediated activation of nuclear factor-kappaB. Mol. Endocrinol. 16, 353–366. doi: 10.1210/mend.16.2.0773

Fox, M. A. (2008). Novel roles for collagens in wiring the vertebrate nervous system. Curr. Opin. Cell Biol. 20, 508–513. doi: 10.1016/j.ceb.2008.05.003

Fraga, M. F., Ballestar, E., Paz, M. F., Ropero, S., Setien, F., Ballestar, M. L., et al. (2005). Epigenetic differences arise during the lifetime of monozygotic twins. Proc. Natl. Acad. Sci. U.S.A. 102, 10604–10609. doi: 10.1073/pnas.0500398102

Furey, T. S. (2012). ChIP-seq and beyond: new and improved methodologies to detect and characterize protein-DNA interactions. Nat. Rev. Genet. 13, 840–852. doi: 10.1038/nrg3306

Fusco, C., Bertani, G., Caricati, G., and Della Giustina, E. (2006). Stress fracture of the peroneal bone secondary to a complex tic. Brain Dev. 28, 52–54. doi: 10.1016/j.braindev.2005.03.009

Gabel, H. W., Kinde, B., Stroud, H., Gilbert, C. S., Harmin, D. A., Kastan, N. R., et al. (2015). Disruption of DNA-methylation-dependent long gene repression in Rett syndrome. Nature 522, 89–93. doi: 10.1038/nature14319

Gade, R., Muhleman, D., Blake, H., MacMurray, J., Johnson, P., Verde, R., et al. (1998). Correlation of length of VNTR alleles at the X-linked MAOA gene and phenotypic effect in Tourette syndrome and drug abuse. Mol. Psychiatry 3, 50–60.

Glessner, J. T., Wang, K., Cai, G., Korvatska, O., Kim, C. E., Wood, S., et al. (2009). Autism genome-wide copy number variation reveals ubiquitin and neuronal genes. Nature 459, 569–573. doi: 10.1038/nature07953

Globisch, D., Münzel, M., Müller, M., Michalakis, S., Wagner, M., Koch, S., et al. (2010). Tissue distribution of 5-hydroxymethylcytosine and search for active demethylation intermediates. PLoS ONE 5:e15367. doi: 10.1371/journal.pone.0015367

Greenblatt, S. M., Liu, F., and Nimer, S. D. (2016). Arginine methyltransferases in normal and malignant hematopoiesis. Exp. Hematol. 44, 435–441. doi: 10.1016/j.exphem.2016.03.009

Gregory, S. G., Connelly, J. J., Towers, A. J., Johnson, J., Biscocho, D., Markunas, C. A., et al. (2009). Genomic and epigenetic evidence for oxytocin receptor deficiency in autism. BMC Med. 7:62. doi: 10.1186/1741-7015-7-62

Grice, D. E., Leckman, J. F., Pauls, D. L., Kurlan, R., Kidd, K. K., Pakstis, A. J., et al. (1996). Linkage disequilibrium between an allele at the dopamine D4 receptor locus and Tourette syndrome, by the transmission-disequilibrium test. Am. J. Hum. Genet. 59, 644–652.

Guo, Y., Su, L., Zhang, J., Lei, J., Deng, X., Xu, H., et al. (2012). Analysis of the BTBD9 and HTR2C variants in Chinese Han patients with Tourette syndrome. Psychiatr. Genet. 22, 300–303. doi: 10.1097/YPG.0b013e32835862b1

Hannon, E., Lunnon, K., Schalkwyk, L., and Mill, J. (2015). Interindividual methylomic variation across blood, cortex, and cerebellum: implications for epigenetic studies of neurological and neuropsychiatric phenotypes. Epigenetics 10, 1024–1032. doi: 10.1080/15592294.2015.1100786

Hark, A. T., Schoenherr, C. J., Katz, D. J., Ingram, R. S., Levorse, J. M., and Tilghman, S. M. (2000). CTCF mediates methylation-sensitive enhancer-blocking activity at the H19/Igf2 locus. Nature 405, 486–489. doi: 10.1038/35013106

Hashimoto, H., Zhang, X., Vertino, P. M., and Cheng, X. (2015). The mechanisms of generation, recognition, and Erasure of DNA 5-Methylcytosine and Thymine oxidations. J. Biol. Chem. 290, 20723–20733. doi: 10.1074/jbc.R115.656884

Herzberg, I., Valencia-Duarte, A. V., Kay, V. A., White, D. J., Müller, H., Rivas, I. C., et al. (2010). Association of DRD2 variants and Gilles de la Tourette syndrome in a family-based sample from a South American population isolate. Psychiatr. Genet. 20, 179–183. doi: 10.1097/YPG.0b013e32833a215a

Hoffmann, A., and Spengler, D. (2014). DNA memories of early social life. Neuroscience 264, 64–75. doi: 10.1016/j.neuroscience.2012.04.003

Huang, W., Payne, T. J., Ma, J. Z., Beuten, J., Dupont, R. T., Inohara, N., et al. (2009). Significant association of ANKK1 and detection of a functional polymorphism with nicotine dependence in an African-American sample. Neuropsychopharmacology 34, 319–330. doi: 10.1038/npp.2008.37

Huertas-Fernández, I., Gómez-Garre, P., Madruga-Garrido, M., Bernal-Bernal, I., Bonilla-Toribio, M., Martín-Rodríguez, J. F., et al. (2015). GDNF gene is associated with tourette syndrome in a family study. Mov. Disord. 30, 1115–1120. doi: 10.1002/mds.26279

Hutvágner, G., and Zamore, P. D. (2002). A microRNA in a multiple-turnover RNAi enzyme complex. Science 297, 2056–2060. doi: 10.1126/science.1073827

Im, H.-I., and Kenny, P. J. (2012). MicroRNAs in neuronal function and dysfunction. Trends Neurosci. 35, 325–334. doi: 10.1016/j.tins.2012.01.004

Inai, A., Tochigi, M., Kuwabara, H., Nishimura, F., Kato, K., Eriguchi, Y., et al. (2015). Analysis of SLITRK1 in Japanese patients with Tourette syndrome using a next-generation sequencer. Psychiatr. Genet. 25, 256–258. doi: 10.1097/YPG.0000000000000104

International Molecular Genetic Study of Autism Consortium (1998). A full genome screen for autism with evidence for linkage to a region on chromosome 7q. Hum. Mol. Genet. 7, 571–578.

Jack, A., Connelly, J. J., and Morris, J. P. (2012). DNA methylation of the oxytocin receptor gene predicts neural response to ambiguous social stimuli. Front. Hum. Neurosci. 6:280. doi: 10.3389/fnhum.2012.00280

Jahan, S., and Davie, J. R. (2015). Protein arginine methyltransferases (PRMTs): role in chromatin organization. Adv. Biol. Regul. 57, 173–184. doi: 10.1016/j.jbior.2014.09.003

Janik, P., Berdyński, M., Safranow, K., and Żekanowski, C. (2015). Association of ADORA1 rs2228079 and ADORA2A rs5751876 Polymorphisms with Gilles de la Tourette Syndrome in the Polish Population. PLoS ONE 10:e0136754. doi: 10.1371/journal.pone.0136754

Jeltsch, A., and Jurkowska, R. Z. (2014). New concepts in DNA methylation. Trends Biochem. Sci. 39, 310–318. doi: 10.1016/j.tibs.2014.05.002

Jenuwein, T., and Allis, C. D. (2001). Translating the histone code. Science 293, 1074–1080. doi: 10.1126/science.1063127

Johnson, R. (2012). Long non-coding RNAs in Huntington's disease neurodegeneration. Neurobiol. Dis. 46, 245–254. doi: 10.1016/j.nbd.2011.12.006

Kang, J., Lienhard, M., Pastor, W. A., Chawla, A., Novotny, M., Tsagaratou, A., et al. (2015). Simultaneous deletion of the methylcytosine oxidases Tet1 and Tet3 increases transcriptome variability in early embryogenesis. Proc. Natl. Acad. Sci. U.S.A. 112, E4236–E4245. doi: 10.1073/pnas.1510510112

Kangaspeska, S., Stride, B., Métivier, R., Polycarpou-Schwarz, M., Ibberson, D., Carmouche, R. P., et al. (2008). Transient cyclical methylation of promoter DNA. Nature 452, 112–115. doi: 10.1038/nature06640

Karagiannidis, I., Dehning, S., Sandor, P., Tarnok, Z., Rizzo, R., Wolanczyk, T., et al. (2013). Support of the histaminergic hypothesis in Tourette syndrome: association of the histamine decarboxylase gene in a large sample of families. J. Med. Genet. 50, 760–764. doi: 10.1136/jmedgenet-2013-101637

Karagiannidis, I., Rizzo, R., Tarnok, Z., Wolanczyk, T., Hebebrand, J., Nöthen, M. M., et al. (2012). Replication of association between a SLITRK1 haplotype and Tourette Syndrome in a large sample of families. Mol. Psychiatry 17, 665–668. doi: 10.1038/mp.2011.151

Katz, D. M., Berger-Sweeney, J. E., Eubanks, J. H., Justice, M. J., Neul, J. L., Pozzo-Miller, L., et al. (2012). Preclinical research in Rett syndrome: setting the foundation for translational success. Dis. Model. Mech. 5, 733–745. doi: 10.1242/dmm.011007

Keen-Kim, D., Mathews, C. A., Reus, V. I., Lowe, T. L., Herrera, L. D., Budman, C. L., et al. (2006). Overrepresentation of rare variants in a specific ethnic group may confuse interpretation of association analyses. Hum. Mol. Genet. 15, 3324–3328. doi: 10.1093/hmg/ddl408

Keszler, G., Kruk, E., Kenezloi, E., Tarnok, Z., Sasvari-Szekely, M., and Nemoda, Z. (2014). Association of the tumor necrosis factor -308 A/G promoter polymorphism with Tourette syndrome. Int. J. Immunogenet. 41, 493–498. doi: 10.1111/iji.12147

Kim, V. N., Han, J., and Siomi, M. C. (2009). Biogenesis of small RNAs in animals. Nat. Rev. Mol. Cell Biol. 10, 126–139. doi: 10.1038/nrm2632

Kinney, D. K., Munir, K. M., Crowley, D. J., and Miller, A. M. (2008). Prenatal stress and risk for autism. Neurosci. Biobehav. Rev. 32, 1519–1532. doi: 10.1016/j.neubiorev.2008.06.004

Kouzarides, T. (2007). Chromatin modifications and their function. Cell 128, 693–705. doi: 10.1016/j.cell.2007.02.005

Kroisel, P. M., Petek, E., Emberger, W., Windpassinger, C., Wladika, W., and Wagner, K. (2001). Candidate region for Gilles de la Tourette syndrome at 7q31. Am. J. Med. Genet. 101, 259–261. doi: 10.1002/1096-8628(20010701)101:3<259::AID-AJMG1374>3.0.CO;2-#

Kulis, M., Heath, S., Bibikova, M., Queirós, A. C., Navarro, A., Clot, G., et al. (2012). Epigenomic analysis detects widespread gene-body DNA hypomethylation in chronic lymphocytic leukemia. Nat. Genet. 44, 1236–1242. doi: 10.1038/ng.2443

Kuo, M. H., and Allis, C. D. (1998). Roles of histone acetyltransferases and deacetylases in gene regulation. Bioessays 20, 615–626. doi: 10.1002/(SICI)1521-1878(199808)20:8<615::AID-BIES4>3.0.CO;2-H

Kurlan, R., Como, P. G., Miller, B., Palumbo, D., Deeley, C., Andresen, E. M., et al. (2002). The behavioral spectrum of tic disorders: a community-based study. Neurology 59, 414–420. doi: 10.1212/WNL.59.3.414

Kwak, C., Dat Vuong, K., and Jankovic, J. (2003). Premonitory sensory phenomenon in Tourette's syndrome. Mov. Disord. 18, 1530–1533. doi: 10.1002/mds.10618

Lawson-Yuen, A., Saldivar, J.-S., Sommer, S., and Picker, J. (2008). Familial deletion within NLGN4 associated with autism and Tourette syndrome. Eur. J. Hum. Genet. 16, 614–618. doi: 10.1038/sj.ejhg.5202006

Lee, C.-C., Chou, I.-C., Tsai, C.-H., Wang, T.-R., Li, T.-C., and Tsai, F.-J. (2005). Dopamine receptor D2 gene polymorphisms are associated in Taiwanese children with Tourette syndrome. Pediatr. Neurol. 33, 272–276. doi: 10.1016/j.pediatrneurol.2005.05.005

Legube, G., and Trouche, D. (2003). Regulating histone acetyltransferases and deacetylases. EMBO Rep. 4, 944–947. doi: 10.1038/sj.embor.embor941

Lei, J., Deng, X., Zhang, J., Su, L., Xu, H., Liang, H., et al. (2012). Mutation screening of the HDC gene in Chinese Han patients with Tourette syndrome. Am. J. Med. Genet. B Neuropsychiatr. Genet. 159B, 72–76. doi: 10.1002/ajmg.b.32003

Lewis, J. D., Meehan, R. R., Henzel, W. J., Maurer-Fogy, I., Jeppesen, P., Klein, F., et al. (1992). Purification, sequence, and cellular localization of a novel chromosomal protein that binds to methylated DNA. Cell 69, 905–914.

Li, S., Papale, L. A., Zhang, Q., Madrid, A., Chen, L., Chopra, P., et al. (2016). Genome-wide alterations in hippocampal 5-hydroxymethylcytosine links plasticity genes to acute stress. Neurobiol. Dis. 86, 99–108. doi: 10.1016/j.nbd.2015.11.010

Lin, M., Pedrosa, E., Shah, A., Hrabovsky, A., Maqbool, S., Zheng, D., et al. (2011). RNA-Seq of human neurons derived from iPS cells reveals candidate long non-coding RNAs involved in neurogenesis and neuropsychiatric disorders. PLoS ONE 6:e23356. doi: 10.1371/journal.pone.0023356

Lin, W.-Y., Lee, C.-C., Liu, H.-P., Chou, I.-C., Sheu, J. J.-C., Wan, L., et al. (2012). Association of genetic variations in X-ray repair cross-complementing group 1 and Tourette syndrome. J. Clin. Lab. Anal. 26, 321–324. doi: 10.1002/jcla.21525

Liu, D., Diorio, J., Tannenbaum, B., Caldji, C., Francis, D., Freedman, A., et al. (1997). Maternal care, hippocampal glucocorticoid receptors, and hypothalamic-pituitary-adrenal responses to stress. Science 277, 1659–1662.

Liu, S., Cui, J., Niu, Z., Yi, M., Zhang, X., Che, F., et al. (2015a). Do obsessive-compulsive disorder and Tourette syndrome share a common susceptibility gene? An association study of the BDNF Val66Met polymorphism in the Chinese Han population. World J. Biol. Psychiatry 16, 602–609. doi: 10.3109/15622975.2015.1012226

Liu, S., Cui, J., Zhang, X., Wu, W., Niu, H., Ma, X., et al. (2014). Variable number tandem repeats in dopamine receptor D4 in Tourette's syndrome. Mov. Disord. 29, 1687–1691. doi: 10.1002/mds.26027

Liu, S., Yu, X., Xu, Q., Cui, J., Yi, M., Zhang, X., et al. (2015b). Support of positive association in family-based genetic analysis between COL27A1 and Tourette syndrome. Sci. Rep. 5:12687. doi: 10.1038/srep12687

Loke, Y. J., Hannan, A. J., and Craig, J. M. (2015). The role of epigenetic change in Autism spectrum disorders. Front. Neurol. 6:107. doi: 10.3389/fneur.2015.00107

Lopez-Atalaya, J. P., Valor, L. M., and Barco, A. (2014). Epigenetic factors in intellectual disability: the Rubinstein-Taybi syndrome as a paradigm of neurodevelopmental disorder with epigenetic origin. Prog. Mol. Biol. Transl. Sci. 128, 139–176. doi: 10.1016/B978-0-12-800977-2.00006-1

Mataix-Cols, D., Isomura, K., Pérez-Vigil, A., Chang, Z., Rück, C., Larsson, K. J., et al. (2015). Familial risks of Tourette Syndrome and Chronic Tic disorders. A population-based cohort study. JAMA Psychiatry 72, 787–793. doi: 10.1001/jamapsychiatry.2015.0627

Matsumoto, N., David, D. E., Johnson, E. W., Konecki, D., Burmester, J. K., Ledbetter, D. H., et al. (2000). Breakpoint sequences of an 1;8 translocation in a family with Gilles de la Tourette syndrome. Eur. J. Hum. Genet. 8, 875–883. doi: 10.1038/sj.ejhg.5200549

Mbadiwe, T., and Millis, R. M. (2013). Epigenetics and Autism. Autism Res. Treat. 2013, 1–9. doi: 10.1155/2013/826156

McGowan, P. O., Suderman, M., Sasaki, A., Huang, T. C. T., Hallett, M., Meaney, M. J., et al. (2011). Broad epigenetic signature of maternal care in the brain of adult rats. PLoS ONE 6:e14739. doi: 10.1371/journal.pone.0014739

Melchior, L., Bertelsen, B., Debes, N. M., Groth, C., Skov, L., Mikkelsen, J. D., et al. (2013). Microduplication of 15q13.3 and Xq21.31 in a family with Tourette syndrome and comorbidities. Am. J. Med. Genet. B Neuropsychiatr. Genet. 162B, 825–831. doi: 10.1002/ajmg.b.32186

Mellios, N., and Sur, M. (2012). The Emerging Role of microRNAs in Schizophrenia and Autism Spectrum Disorders. Front. Psychiatry 3:39. doi: 10.3389/fpsyt.2012.00039

Métivier, R., Gallais, R., Tiffoche, C., Le Péron, C., Jurkowska, R. Z., Carmouche, R. P., et al. (2008). Cyclical DNA methylation of a transcriptionally active promoter. Nature 452, 45–50. doi: 10.1038/nature06544

Meylan, E., and Tschopp, J. (2005). The RIP kinases: crucial integrators of cellular stress. Trends Biochem. Sci. 30, 151–159. doi: 10.1016/j.tibs.2005.01.003

Miranda, D. M., Wigg, K., Kabia, E. M., Feng, Y., Sandor, P., and Barr, C. L. (2009). Association of SLITRK1 to Gilles de la Tourette Syndrome. Am. J. Med. Genet. B Neuropsychiatr. Genet. 150B, 483–486. doi: 10.1002/ajmg.b.30840

Miyoshi, K., Miyoshi, T., and Siomi, H. (2010). Many ways to generate microRNA-like small RNAs: non-canonical pathways for microRNA production. Mol. Genet. Genomics 284, 95–103. doi: 10.1007/s00438-010-0556-1

Mohan, K. N., and Chaillet, J. R. (2013). Cell and molecular biology of DNA methyltransferase 1. Int. Rev. Cell Mol. Biol. 306, 1–42. doi: 10.1016/B978-0-12-407694-5.00001-8

Moore, S. J., Turnpenny, P., Quinn, A., Glover, S., Lloyd, D. J., Montgomery, T., et al. (2000). A clinical study of 57 children with fetal anticonvulsant syndromes. J. Med. Genet. 37, 489–97. doi: 10.1136/jmg.37.7.489

Mössner, R., Müller-Vahl, K. R., Döring, N., and Stuhrmann, M. (2007). Role of the novel tryptophan hydroxylase-2 gene in Tourette syndrome. Mol. Psychiatry 12, 617–619. doi: 10.1038/sj.mp.4002004

Moya, P. R., Wendland, J. R., Rubenstein, L. M., Timpano, K. R., Heiman, G. A., Tischfield, J. A., et al. (2013). Common and rare alleles of the serotonin transporter gene, SLC6A4, associated with Tourette's disorder. Mov. Disord. 28, 1263–1270. doi: 10.1002/mds.25460

Nag, A., Bochukova, E. G., Kremeyer, B., Campbell, D. D., Muller, H., Valencia-Duarte, A. V., et al. (2013). CNV analysis in Tourette syndrome implicates large genomic rearrangements in COL8A1 and NRXN1. PLoS ONE 8:e59061. doi: 10.1371/journal.pone.0059061

Neville, M. J., Johnstone, E. C., and Walton, R. T. (2004). Identification and characterization of ANKK1: a novel kinase gene closely linked to DRD2 on chromosome band 11q23.1. Hum. Mutat. 23, 540–545. doi: 10.1002/humu.20039

Nishimoto, Y., Nakagawa, S., Hirose, T., Okano, H. J., Takao, M., Shibata, S., et al. (2013). The long non-coding RNA nuclear-enriched abundant transcript 1_2 induces paraspeckle formation in the motor neuron during the early phase of amyotrophic lateral sclerosis. Mol. Brain 6:31. doi: 10.1186/1756-6606-6-31

O'Roak, B. J., Morgan, T. M., Fishman, D. O., Saus, E., Alonso, P., Gratacòs, M., et al. (2010). Additional support for the association of SLITRK1 var321 and Tourette syndrome. Mol. Psychiatry 15, 447–450. doi: 10.1038/mp.2009.105

Oberlé, I., Rousseau, F., Heitz, D., Kretz, C., Devys, D., Hanauer, A., et al. (1991). Instability of a 550-base pair DNA segment and abnormal methylation in fragile X syndrome. Science 252, 1097–1102.

Oey, H., and Whitelaw, E. (2014). On the meaning of the word “epimutation”. Trends Genet. 30, 519–520. doi: 10.1016/j.tig.2014.08.005

Ornoy, A., Weinstein-Fudim, L., and Ergaz, Z. (2015). Prenatal factors associated with autism spectrum disorder (ASD). Reprod. Toxicol. 56, 155–169. doi: 10.1016/j.reprotox.2015.05.007

Pace, J. M., Corrado, M., Missero, C., and Byers, P. H. (2003). Identification, characterization and expression analysis of a new fibrillar collagen gene, COL27A1. Matrix Biol. 22, 3–14. doi: 10.1016/S0945-053X(03)00007-6

Palma-Gudiel, H., Córdova-Palomera, A., Leza, J. C., and Fañanás, L. (2015). Glucocorticoid receptor gene (NR3C1) methylation processes as mediators of early adversity in stress-related disorders causality: a critical review. Neurosci. Biobehav. Rev. 55, 520–535. doi: 10.1016/j.neubiorev.2015.05.016

Paschou, P. (2013). The genetic basis of Gilles de la Tourette Syndrome. Neurosci. Biobehav. Rev. 37, 1026–1039. doi: 10.1016/j.neubiorev.2013.01.016

Paschou, P., Feng, Y., Pakstis, A. J., Speed, W. C., DeMille, M. M., Kidd, J. R., et al. (2004). Indications of linkage and association of Gilles de la Tourette syndrome in two independent family samples: 17q25 is a putative susceptibility region. Am. J. Hum. Genet. 75, 545–560. doi: 10.1086/424389

Paschou, P., Stylianopoulou, E., Karagiannidis, I., Rizzo, R., Tarnok, Z., Wolanczyk, T., et al. (2012). Evaluation of the LIM homeobox genes LHX6 and LHX8 as candidates for Tourette syndrome. Genes. Brain. Behav. 11, 444–451. doi: 10.1111/j.1601-183X.2012.00778.x

Paschou, P., Yu, D., Gerber, G., Evans, P., Tsetsos, F., Davis, L. K., et al. (2014). Genetic association signal near NTN4 in Tourette syndrome. Ann. Neurol. 76, 310–315. doi: 10.1002/ana.24215

Patel, C., Cooper-Charles, L., McMullan, D. J., Walker, J. M., Davison, V., and Morton, J. (2011). Translocation breakpoint at 7q31 associated with tics: further evidence for IMMP2L as a candidate gene for Tourette syndrome. Eur. J. Hum. Genet. 19, 634–639. doi: 10.1038/ejhg.2010.238

Pauls, D. L., Fernandez, T. V., Mathews, C. A., State, M. W., and Scharf, J. M. (2014). The inheritance of Tourette Disorder: a review. J. Obsessive. Compuls. Relat. Disord. 3, 380–385. doi: 10.1016/j.jocrd.2014.06.003

Penagarikano, O., Mulle, J. G., and Warren, S. T. (2007). The pathophysiology of fragile x syndrome. Annu. Rev. Genomics Hum. Genet. 8, 109–129. doi: 10.1146/annurev.genom.8.080706.092249

Petazzi, P., Sandoval, J., Szczesna, K., Jorge, O. C., Roa, L., Sayols, S., et al. (2013). Dysregulation of the long non-coding RNA transcriptome in a Rett syndrome mouse model. RNA Biol. 10, 1197–1203. doi: 10.4161/rna.24286

Petek, E., Windpassinger, C., Vincent, J. B., Cheung, J., Boright, A. P., Scherer, S. W., et al. (2001). Disruption of a novel gene (IMMP2L) by a breakpoint in 7q31 associated with Tourette syndrome. Am. J. Hum. Genet. 68, 848–858. doi: 10.1086/319523

Plazas-Mayorca, M. D., and Vrana, K. E. (2011). Proteomic investigation of epigenetics in neuropsychiatric disorders: a missing link between genetics and behavior? J. Proteome Res. 10, 58–65. doi: 10.1021/pr100463y

Polderman, T. J. C., Benyamin, B., de Leeuw, C. A., Sullivan, P. F., van Bochoven, A., Visscher, P. M., et al. (2015). Meta-analysis of the heritability of human traits based on fifty years of twin studies. Nat. Genet. 47, 702–709. doi: 10.1038/ng.3285

Ponce, G., Pérez-González, R., Aragüés, M., Palomo, T., Rodríguez-Jiménez, R., Jiménez-Arriero, M. A., et al. (2009). The ANKK1 kinase gene and psychiatric disorders. Neurotox. Res. 16, 50–59. doi: 10.1007/s12640-009-9046-9

Pradhan, S., Bacolla, A., Wells, R. D., and Roberts, R. J. (1999). Recombinant human DNA (cytosine-5) methyltransferase. I. Expression, purification, and comparison of de novo and maintenance methylation. J. Biol. Chem. 274, 33002–330010.

Proenca, C. C., Gao, K. P., Shmelkov, S. V., Rafii, S., and Lee, F. S. (2011). Slitrks as emerging candidate genes involved in neuropsychiatric disorders. Trends Neurosci. 34, 143–153. doi: 10.1016/j.tins.2011.01.001

Prokhortchouk, A., Hendrich, B., Jørgensen, H., Ruzov, A., Wilm, M., Georgiev, G., et al. (2001). The p120 catenin partner Kaiso is a DNA methylation-dependent transcriptional repressor. Genes Dev. 15, 1613–1618. doi: 10.1101/gad.198501

Prontera, P., Napolioni, V., Ottaviani, V., Rogaia, D., Fusco, C., Augello, B., et al. (2014). DPP6 gene disruption in a family with Gilles de la Tourette syndrome. Neurogenetics 15, 237–242. doi: 10.1007/s10048-014-0418-9

Quiroga, A. D., and Lehner, R. (2011). Role of endoplasmic reticulum neutral lipid hydrolases. Trends Endocrinol. Metab. 22, 218–225. doi: 10.1016/j.tem.2011.03.003

Rao, S.-Q., Hu, H.-L., Ye, N., Shen, Y., and Xu, Q. (2015). Genetic variants in long non-coding RNA MIAT contribute to risk of paranoid schizophrenia in a Chinese Han population. Schizophr. Res. 166, 125–130. doi: 10.1016/j.schres.2015.04.032

Rivière, J.-B., Xiong, L., Levchenko, A., St-Onge, J., Gaspar, C., Dion, Y., et al. (2009). Association of intronic variants of the BTBD9 gene with Tourette syndrome. Arch. Neurol. 66, 1267–1272. doi: 10.1001/archneurol.2009.213

Rizzo, R., Ragusa, M., Barbagallo, C., Sammito, M., Gulisano, M., Calì, P. V., et al. (2015). Circulating miRNAs profiles in tourette syndrome: molecular data and clinical implications. Mol. Brain 8, 44. doi: 10.1186/s13041-015-0133-y

Robertson, M. M. (2000). Tourette syndrome, associated conditions and the complexities of treatment. Brain 123, 425–462. doi: 10.1093/brain/123.3.425

Robertson, M. M. (2012). The Gilles de la Tourette syndrome: the current status. Arch. Dis. Child. Educ. Pract. Ed. 97, 166–175. doi: 10.1136/archdischild-2011-300585

Robertson, M. M., and Orth, M. (2006). Behavioral and affective disorders in Tourette syndrome. Adv. Neurol. 99, 39–60.

Sakamoto, S., Potla, R., and Larner, A. C. (2004). Histone deacetylase activity is required to recruit RNA polymerase II to the promoters of selected interferon-stimulated early response genes. J. Biol. Chem. 279, 40362–40367. doi: 10.1074/jbc.M406400200

Salta, E., and De Strooper, B. (2012). Non-coding RNAs with essential roles in neurodegenerative disorders. Lancet. Neurol. 11, 189–200. doi: 10.1016/S1474-4422(11)70286-1

Sánchez Delgado, M., Camprubí, C., Tümer, Z., Martínez, F., Milà, M., and Monk, D. (2014). Screening individuals with intellectual disability, autism and Tourette's syndrome for KCNK9 mutations and aberrant DNA methylation within the 8q24 imprinted cluster. Am. J. Med. Genet. B Neuropsychiatr. Genet. 165B, 472–478. doi: 10.1002/ajmg.b.32250

Santos-Rosa, H., and Caldas, C. (2005). Chromatin modifier enzymes, the histone code and cancer. Eur. J. Cancer 41, 2381–2402. doi: 10.1016/j.ejca.2005.08.010

Scharf, J. M., Yu, D., Mathews, C. A., Neale, B. M., Stewart, S. E., Fagerness, J. A., et al. (2013). Genome-wide association study of Tourette's syndrome. Mol. Psychiatry 18, 721–728. doi: 10.1038/mp.2012.69

Schermelleh, L., Haemmer, A., Spada, F., Rösing, N., Meilinger, D., Rothbauer, U., et al. (2007). Dynamics of Dnmt1 interaction with the replication machinery and its role in postreplicative maintenance of DNA methylation. Nucleic Acids Res. 35, 4301–4312. doi: 10.1093/nar/gkm432

Sellier, C., Freyermuth, F., Tabet, R., Tran, T., He, F., Ruffenach, F., et al. (2013). Sequestration of DROSHA and DGCR8 by expanded CGG RNA repeats alters microRNA processing in fragile X-associated tremor/ataxia syndrome. Cell Rep. 3, 869–880. doi: 10.1016/j.celrep.2013.02.004

Shen, C.-P., Chou, I.-C., Liu, H.-P., Lee, C.-C., Tsai, Y., Wu, B.-T., et al. (2014). Association of glutathione S-transferase P1 (GSTP1) polymorphism with Tourette syndrome in Taiwanese patients. Genet. Test. Mol. Biomarkers 18, 41–44. doi: 10.1089/gtmb.2013.0245

Smallwood, S. A., Tomizawa, S.-I., Krueger, F., Ruf, N., Carli, N., Segonds-Pichon, A., et al. (2011). Dynamic CpG island methylation landscape in oocytes and preimplantation embryos. Nat. Genet. 43, 811–814. doi: 10.1038/ng.864

Smith, Z. D., and Meissner, A. (2013). DNA methylation: roles in mammalian development. Nat. Rev. Genet. 14, 204–220. doi: 10.1038/nrg3354

Snider, L. A., and Swedo, S. E. (2003). Post-streptococcal autoimmune disorders of the central nervous system. Curr. Opin. Neurol. 16, 359–365. doi: 10.1097/01.wco.0000073938.19076.31

Stillman, A. A., Krsnik, Z., Sun, J., Rasin, M.-R., State, M. W., Sestan, N., et al. (2009). Developmentally regulated and evolutionarily conserved expression of SLITRK1 in brain circuits implicated in Tourette syndrome. J. Comp. Neurol. 513, 21–37. doi: 10.1002/cne.21919

Sundaram, S. K., Huq, A. M., Wilson, B. J., and Chugani, H. T. (2010). Tourette syndrome is associated with recurrent exonic copy number variants. Neurology 74, 1583–1590. doi: 10.1212/WNL.0b013e3181e0f147

Swain, J. E., Scahill, L., Lombroso, P. J., King, R. A., and Leckman, J. F. (2007). Tourette syndrome and tic disorders: a decade of progress. J. Am. Acad. Child Adolesc. Psychiatry 46, 947–968. doi: 10.1097/chi.0b013e318068fbcc

Swedo, S. E., Leonard, H. L., Garvey, M., Mittleman, B., Allen, A. J., Perlmutter, S., et al. (1998). Pediatric autoimmune neuropsychiatric disorders associated with streptococcal infections: clinical description of the first 50 cases. Am. J. Psychiatry 155, 264–271. doi: 10.1176/ajp.155.2.264

Talkowski, M. E., Maussion, G., Crapper, L., Rosenfeld, J. A., Blumenthal, I., Hanscom, C., et al. (2012). Disruption of a large intergenic noncoding RNA in subjects with neurodevelopmental disabilities. Am. J. Hum. Genet. 91, 1128–1134. doi: 10.1016/j.ajhg.2012.10.016

Tan, L., Yu, J.-T., Hu, N., and Tan, L. (2013). Non-coding RNAs in Alzheimer's disease. Mol. Neurobiol. 47, 382–393. doi: 10.1007/s12035-012-8359-5

Tarnok, Z., Ronai, Z., Gervai, J., Kereszturi, E., Gadoros, J., Sasvari-Szekely, M., et al. (2007). Dopaminergic candidate genes in Tourette syndrome: association between tic severity and 3′ UTR polymorphism of the dopamine transporter gene. Am. J. Med. Genet. B Neuropsychiatr. Genet. 144B, 900–905. doi: 10.1002/ajmg.b.30517

Tordjman, S., Somogyi, E., Coulon, N., Kermarrec, S., Cohen, D., Bronsard, G., et al. (2014). Gene × Environment interactions in autism spectrum disorders: role of epigenetic mechanisms. Front. Psychiatry 5:53. doi: 10.3389/fpsyt.2014.00053

Urdinguio, R. G., Fernandez, A. F., Lopez-Nieva, P., Rossi, S., Huertas, D., Kulis, M., et al. (2010). Disrupted microRNA expression caused by Mecp2 loss in a mouse model of Rett syndrome. Epigenetics 5, 656–663. doi: 10.4161/epi.5.7.13055

Van Tol, H. H., Wu, C. M., Guan, H. C., Ohara, K., Bunzow, J. R., Civelli, O., et al. (1992). Multiple dopamine D4 receptor variants in the human population. Nature 358, 149–152. doi: 10.1038/358149a0

Vaquero, A., Sternglanz, R., and Reinberg, D. (2007). NAD+-dependent deacetylation of H4 lysine 16 by class III HDACs. Oncogene 26, 5505–5520. doi: 10.1038/sj.onc.1210617

Vereczkei, A., Demetrovics, Z., Szekely, A., Sarkozy, P., Antal, P., Szilagyi, A., et al. (2013). Multivariate analysis of dopaminergic gene variants as risk factors of heroin dependence. PLoS ONE 8:e66592. doi: 10.1371/journal.pone.0066592

Verkerk, A. J. M. H., Mathews, C. A., Joosse, M., Eussen, B. H. J., Heutink, P., and Oostra, B. A. (2003). CNTNAP2 is disrupted in a family with Gilles de la Tourette syndrome and obsessive compulsive disorder. Genomics 82, 1–9. doi: 10.1016/S0888-7543(03)00097-1

Vigneswaran, K., Neill, S., and Hadjipanayis, C. G. (2015). Beyond the World Health Organization grading of infiltrating gliomas: advances in the molecular genetics of glioma classification. Ann. Transl. Med. 3, 95. doi: 10.3978/j.issn.2305-5839.2015.03.57

Vrijenhoek, T., Buizer-Voskamp, J. E., van der Stelt, I., Strengman, E., Sabatti, C., Geurts van Kessel, A., et al. (2008). Recurrent CNVs disrupt three candidate genes in schizophrenia patients. Am. J. Hum. Genet. 83, 504–510. doi: 10.1016/j.ajhg.2008.09.011

Weaver, I. C. G., Cervoni, N., Champagne, F. A., D'Alessio, A. C., Sharma, S., Seckl, J. R., et al. (2004). Epigenetic programming by maternal behavior. Nat. Neurosci. 7, 847–854. doi: 10.1038/nn1276

Wong, C. C. Y., Meaburn, E. L., Ronald, A., Price, T. S., Jeffries, A. R., Schalkwyk, L. C., et al. (2014). Methylomic analysis of monozygotic twins discordant for autism spectrum disorder and related behavioural traits. Mol. Psychiatry 19, 495–503. doi: 10.1038/mp.2013.41

Wu, B.-T., Lin, W.-Y., Chou, I.-C., Liu, H.-P., Lee, C.-C., Tsai, Y., et al. (2013a). Association of poly(ADP-ribose) polymerase-1 polymorphism with Tourette syndrome. Neurol. Sci. 34, 1911–1916. doi: 10.1007/s10072-013-1405-x

Wu, B.-T., Lin, W.-Y., Chou, I.-C., Liu, H.-P., Lee, C.-C., Tsai, Y., et al. (2013b). Association of tyrosyl-DNA phosphodiesterase 1 polymorphism with Tourette syndrome in Taiwanese patients. J. Clin. Lab. Anal. 27, 323–327. doi: 10.1002/jcla.21606

Wu, H., D'Alessio, A. C., Ito, S., Wang, Z., Cui, K., Zhao, K., et al. (2011). Genome-wide analysis of 5-hydroxymethylcytosine distribution reveals its dual function in transcriptional regulation in mouse embryonic stem cells. Genes Dev. 25, 679–684. doi: 10.1101/gad.2036011

Wu, H., Tao, J., Chen, P. J., Shahab, A., Ge, W., Hart, R. P., et al. (2010). Genome-wide analysis reveals methyl-CpG-binding protein 2-dependent regulation of microRNAs in a mouse model of Rett syndrome. Proc. Natl. Acad. Sci. U.S.A. 107, 18161–18166. doi: 10.1073/pnas.1005595107

Xu, M., Li, L., Ohtsu, H., and Pittenger, C. (2015). Histidine decarboxylase knockout mice, a genetic model of Tourette syndrome, show repetitive grooming after induced fear. Neurosci. Lett. 595, 50–53. doi: 10.1016/j.neulet.2015.03.067

Yamagata, Y., Parietti, V., Stockholm, D., Corre, G., Poinsignon, C., Touleimat, N., et al. (2012a). Lentiviral transduction of CD34(+) cells induces genome-wide epigenetic modifications. PLoS ONE 7:e48943. doi: 10.1371/journal.pone.0048943

Yamagata, Y., Szabó, P., Szüts, D., Bacquet, C., Arànyi, T., and Páldi, A. (2012b). Rapid turnover of DNA methylation in human cells. Epigenetics 7, 141–145. doi: 10.4161/epi.7.2.18906

Yang, Y. C., Ricciardi, S., Ciarletta, A., Calvetti, J., Kelleher, K., and Clark, S. C. (1989). Expression cloning of cDNA encoding a novel human hematopoietic growth factor: human homologue of murine T-cell growth factor P40. Blood 74, 1880–1884.

Yoon, D. Y., Rippel, C. A., Kobets, A. J., Morris, C. M., Lee, J. E., Williams, P. N., et al. (2007). Dopaminergic polymorphisms in Tourette syndrome: association with the DAT gene (SLC6A3). Am. J. Med. Genet. B Neuropsychiatr. Genet. 144B, 605–610. doi: 10.1002/ajmg.b.30466

Yuan, A., Su, L., Yu, S., Li, C., Yu, T., and Sun, J. (2015). Association between DRD2/ANKK1 TaqIA Polymorphism and Susceptibility with Tourette Syndrome: a meta-analysis. PLoS ONE 10:e0131060. doi: 10.1371/journal.pone.0131060

Zee, B. M., Levin, R. S., Dimaggio, P. A., and Garcia, B. A. (2010). Global turnover of histone post-translational modifications and variants in human cells. Epigenetics Chromatin 3:22. doi: 10.1186/1756-8935-3-22

Zeng, J., Nagrajan, H. K., and Yi, S. V. (2014). Fundamental diversity of human CpG islands at multiple biological levels. Epigenetics 9, 483–491. doi: 10.4161/epi.27654

Ziats, M. N., and Rennert, O. M. (2012). Aberrant Expression of Long Noncoding RNAs in Autistic Brain. J. Mol. Neurosci. 49, 589–593. doi: 10.1007/s12031-012-9880-8

Ziegler, C., Dannlowski, U., Bräuer, D., Stevens, S., Laeger, I., Wittmann, H., et al. (2015). Oxytocin receptor gene methylation: converging multilevel evidence for a role in social anxiety. Neuropsychopharmacology 40, 1528–1538. doi: 10.1038/npp.2015.2

Zilhão, N. R., Padmanabhuni, S. S., Pagliaroli, L., Barta, C., Smit, D. J. A., Cath, D., et al. (2015). Epigenome-wide association study of Tic disorders. Twin Res. Hum. Genet. 18, 1–11. doi: 10.1017/thg.2015.72

Ziller, M. J., Müller, F., Liao, J., Zhang, Y., Gu, H., Bock, C., et al. (2011). Genomic distribution and inter-sample variation of non-CpG methylation across human cell types. PLoS Genet. 7:e1002389. doi: 10.1371/journal.pgen.1002389