DYT6 dystonia: Review of the literature and creation of the UMD locusspecific database (LSDB) for mutations in the THAP1 gene HAL Id: hal-01670069 https://hal.archives-ouvertes.fr/hal-01670069 Submitted on 21 Dec 2017 HAL is a multi-disciplinary open access archive for the deposit and dissemination of sci- entific research documents, whether they are pub- lished or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers. L’archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d’enseignement et de recherche français ou étrangers, des laboratoires publics ou privés. DYT6 dystonia: Review of the literature and creation of the UMD locus-specific database (LSDB) for mutations in the THAP1 gene Arnaud Blanchard, Vuthy Ea, Agathe Roubertie, Mélanie Martin, Coline Coquart, Mireille Claustres, Christophe Béroud, Gwenaelle Collod-Beroud To cite this version: Arnaud Blanchard, Vuthy Ea, Agathe Roubertie, Mélanie Martin, Coline Coquart, et al.. DYT6 dystonia: Review of the literature and creation of the UMD locus-specific database (LSDB) for mutations in the THAP1 gene: review. Human Mutation, Wiley, 2011, 32 (11), pp.1213 - 1224. �10.1002/humu.21564�. �hal-01670069� https://hal.archives-ouvertes.fr/hal-01670069 https://hal.archives-ouvertes.fr DATABASES www.hgvs.org DYT6 Dystonia: Review of the Literature and Creation of the UMD Locus-Specific Database (LSDB) for Mutations in the THAP1 Gene Arnaud Blanchard,1,2 Vuthy Ea,1,2 Agathe Roubertie,1,2,3 Mélanie Martin,1,2 Coline Coquart,4 Mireille Claustres,1,2,4 Christophe Béroud,1,2,4 and Gwenaëlle Collod-Béroud1,2∗ 1 INSERM U827, Montpellier, F-34000, France; 2 Université Montpellier1, UFR de Médecine, Montpellier, F-34000, France; 3 CHU Montpellier, Hôpital Gui de Chauliac, Service de Neuropédiatrie, Montpellier, F-34000, France; 4 CHU Montpellier, Hôpital Arnaud de Villeneuve, Laboratoire de Génétique Moléculaire, Montpellier, F-34000, France Communicated by Alastair F. Broun Received 24 January 2011; accepted revised manuscript 20 June 2011. Published online 26 July 2011 in Wiley Online Library (www.wiley.com/humanmutation).DOI: 10.1002/humu.21564 ABSTRACT: By family-based screening, first Fuchs and then many other authors showed that mutations in THAP1 (THAP [thanatos-associated protein] domain- containing, apoptosis-associated protein 1) account for a substantial proportion of familial, early-onset, nonfo- cal, primary dystonia cases (DYT6 dystonia). THAP1 is the first transcriptional factor involved in primary dysto- nia and the hypothesis of a transcriptional deregulation, which was primarily proposed for the X-linked dystonia- parkinsonism (DYT3 dystonia), provided thus a new way to investigate the possible mechanism underlying the de- velopment of dystonic movements. Currently, 56 fami- lies present with a THAP1 mutation; however, no geno- type/phenotype relationship has been found. Therefore, we carried out a systematic review of the literature on the THAP1 gene to colligate all reported patients with a spe- cific THAP1 mutation and the associated clinical signs in order to describe the broad phenotypic continuum of this disorder. To facilitate the comparison of the identified mutations, we created a Locus-Specific Database (UMD- THAP1 LSDB) available at http://www.umd.be/THAP1/. Currently, the database lists 56 probands and 43 rela- tives with the associated clinical phenotype when avail- able. The identification of a larger number of THAP1 mutations and collection of high-quality clinical informa- tion for each described mutation through international collaborative effort will help investigating the structure– function and genotype–phenotype correlations in DYT6 dystonia. Hum Mutat 00:1–12, 2011. C© 2011 Wiley-Liss, Inc. ∗Correspondence to: Gwenaëlle Collod-Béroud, INSERM U827, Institut Universitaire de Recherche Clinique, 641 av du doyen Gaston Giraud, 34093 Montpellier cedex 5, France. E-mail: gwenaelle.collod-beroud@inserm.fr Contract grant sponsors: Université Montpellier1; INSERM; French Ministry of Health (National PHRC 2007-A00614-49); AMADYS-LFCD, Alliance France Dystonie; Lions Club La Grande Motte; the French Dystonia Network; the European Community’s Sev- enth Framework Programme (FP7/2007–2013 under grant agreement no. 200754—the GEN2PHEN project). KEY WORDS: dystonia; DYT6; THAP1, database; UMD Introduction Dystonia defines a heterogeneous group of movement disorders due to central nervous system (CNS) dysfunction leading to sus- tained involuntary muscle contractions that cause abnormal pos- tures and twisting movements. Since secondary dystonia is often linked to striatal injuries, it has been hypothesized that hereditary dystonia also could be due to a dysfunction of these structures and more broadly of the basal ganglia. Accordingly, many stud- ies have reported electrophysiological, metabolic, and structural abnormalities in these brain areas [Vidailhet et al., 2009]. How- ever, basal ganglia are not the only brain region that has been found to be altered in hereditary dystonia. Indeed, dysfunction of cerebral cortex, cerebellum, spinal cord, and other structures that are mainly involved in motor function has also been reported [Breakefield et al., 2008; Defazio et al., 2007]. Therefore, the con- cept that dystonia is a disorder linked to a specific CNS struc- ture appears to be too restrictive and it should rather be consid- ered as the result of a general disturbance of the motor circuit functions. Hereditary dystonias are clinically and genetically heterogeneous. To date, 20 different genetically forms are known (DYT1 to DYT21; with DYT14 = DYT5) and mutations in nine genes have been in- volved [Breakefield et al., 2008; Norgren et al., 2011]. All the inheri- tance modes are found (autosomic recessive, autosomic dominant, and X-linked) but no epistatic relation between these genes was discovered except for DYT1 and DYT6 forms [Gavarini et al., 2010; Kaiser et al., 2010]. Since their products have different cell functions, it is very difficult up to now to propose a scheme of cell dysfunction common to all the dystonia forms and particularly to propose a unitary pathophysiological mechanism even if links with dopamine deregulation are often found [Blanchard et al., 2010; Park et al., 2005]. The last identified gene, THAP1 for THAP [thanatos-associated protein] domain-containing, apoptosis-associated protein 1, was first found to be mutated in Amish–Mennonite families with pri- mary dystonia and its implication was then extended to other C© 2011 WILEY-LISS, INC. populations [Almasy et al., 1997; Fuchs et al., 2009]. THAP1 muta- tions cause DYT6 dystonia, an autosomal dominant primary form with about 60% penetrance. The main clinical characteristics of this form are early onset (mean = 17.8 years of age [calculated from 78 patients with available data]) and symptoms often localized to one upper limb at the beginning with tendency to extend to other body regions. Fifty-three different mutations in the THAP1 gene have been re- ported so far in 56 families (Table 1). These mutations are mainly private, essentially missense or out-of-frame deletions, usually non- recurrent, and widely distributed throughout the gene. To date, no genotype/phenotype relationship has been observed. From this per- spective, we have compiled a database that gather all the available information on this gene and on the patients with THAP1 mutations and created a software package to facilitate the mutational analysis of THAP1. This database is available at http://www.umd.be/THAP1/. THAP1 Expression Two spliced mRNA variants that produce functional proteins have been reported (THAP1a: CCDS6136 and THAP1b: CCDS6137) [Girard et al., 2010]. The first 2.2-kb isoform contains three exons (Fig. 1), whereas the second corresponds to an alternatively spliced isoform that lacks exon 2 (2-kb mRNA). This second isoform en- codes a truncated THAP1 protein without the C-terminus of the THAP domain. The two isoforms are expressed in many tissues, suggesting that THAP1 has a widespread (although not ubiquitous) distribution in humans. In mouse brain tissue, immunoblot analysis revealed a highest concentration in embryonic whole brain tissue (at E16), which declines after birth in the different tested brain regions (at P60) [Gavarini et al., 2010]. THAP1a, which contains the THAP domain with DNA-binding properties, is the most studied. THAP1 Is a Transcription Factor THAP1 was first identified by two-hybrid screen of a high en- dothelial venule cell cDNA library using a cytokine bait [Roussigne et al., 2003]. THAP1 is a 213 amino-acid-long protein characterized by an N-terminal THAP domain (amino acids 1 to 81) [Bessiere et al., 2008] with DNA-binding properties, followed by a proline- rich region (amino acids 90 to 110), and a nuclear localization signal (amino acids 146 to 162) (Fig. 1A). The protein is only found in the nucleus where it can display both a diffuse distribution and a discrete localization in nuclear dots. THAP1 colocalization with DAXX, a well-characterized pro- tein that is expressed in promyelocytic leukemia (PML) nuclear bodies (NBs), indicates that these subnuclear regions are PML NBs. PML NBs are spheres of 0.1–1.0 μm, heterogeneous in composition, mobility, and function, that are present in most mammalian cell nu- clei [Lallemand-Breitenbach and de The, 2010]. The PML protein is the key organizer that recruits an ever-growing number of proteins, residing constitutively or more often transiently. PML NBs facil- itate partner protein posttranslational modification (acetylation, sumoylation, or phosphorylation) resulting in partner sequestra- tion, activation, or degradation. They are proposed to fine-tune a wide variety of processes: induction of apoptosis and cellular senes- cence, inhibition of proliferation, maintenance of genomic stability, and antiviral responses [Bernardi and Pandolfi, 2007]. THAP1 interacts with the prostate apoptosis response-4 (PAR-4, HGNC-approved gene symbols Pawr), a pro-apoptotic protein that acts as a transcriptional regulator [Roussigne et al., 2003]. PAWR is present in low amounts in dendrites and synapses throughout the brain under normal conditions and its amount in- creases rapidly in neurons subjected to several triggers of apop- tosis including trophic factor deprivation, oxidative stress, and excitotoxins [Duan et al., 1999a; Guo et al., 1998]. Pawr induc- tion has been linked to neuronal death in various neurodegenera- tive diseases as Parkinson’s [Duan et al., 1999b], Alzheimer [Guo et al., 1998], or amyotrophic lateral sclerosis [Pedersen et al., 2000]. More interestingly, Pawr has been shown to serve other important function in nervous system. It interacts with D2 dopamine recep- tor (D2DR) suggesting its potential role in modulating calcium- mediated dopaminergic signaling [Park et al., 2005]. THAP1 and PAWR exhibit pro-apoptotic activities (they increase the apoptosis sensitivity of mouse 3T3 fibroblasts to TNF-α and serum withdrawal when overexpressed) and they are both indirectly recruited to PML NBs by PML [Roussigne et al., 2003]. Roussigne and coworkers proposed that PML NBs regulate the functions of THAP1 and PAWR. They may be the site of assembly and/or post- translational modifications of the THAP1/PAWR complex or they may serve as a site of storage or degradation, thus regulating the free nucleoplasmic pool of PAWR and THAP1. Extensive database search using the THAP domain identified about 100 THAP proteins among seven model animal organisms [Clouaire et al., 2005]. Interestingly, some of these proteins are tightly linked to the pRB (retinoblastoma protein)/E2F pathway, which regulates commitment of mammalian cells to DNA replica- tion. The THAP domain is present in the zebrafish ortholog of E2F6 and in five Caenorhabditis elegans proteins (LIN-36, LIN-15B, LIN- 15A, HIM-17, and GON-14), which are known to be genetically linked to LIN-35, the only Rb ortholog found in C. elegans. Based on these data, Cayrol et al. [2007] investigated whether THAP1 could be linked to the pRB/E2F pathway. They found that silenc- ing as well as overexpression of THAP1 led to cell cycle arrest at the G1/S transition. The similar effects following THAP1 silenc- ing/overexpression suggest that, like many other proteins, THAP1 needs an optimal range of concentration to fulfill its physiologi- cal roles. DNA microarray analysis showed in both cases (upreg- ulation or downregulation of THAP1) a reduction of the mRNA levels of cell cycle regulators and about 40 pRB/E2F-target genes [Cayrol et al., 2007], including RRM1 (ribonucleotide reductase M1), BIRC5 (survivin), CDC2 (cell division cycle 2), CCNB1 (cy- clin B1) [Lacroix, 2007]. RRM1 is activated at the G1/S transition and contains two predicted THAP1-binding sites (THABS) in its promoter. DNase I footprinting, electrophoretic mobility shift assay (EMSA) protein/DNA binding, and chromatin immunoprecipita- tion (ChIP) in proliferating endothelial cells demonstrated that en- dogenous THAP1 directly binds to the promoter of RRM1 [Cayrol et al., 2007]. Moreover, Mazars et al. [2010] showed that THAP1 recruits HCF-1 (Host cell factor 1) to the RRM1 promoter to allow its expression, thus adding further evidence for a role of THAP1 in cell cycle regulation. Although the list of cofactors that participate in THAP1 regulation of target genes is certainly not exhaustive, the available data clearly link THAP1 to apoptosis (through its interac- tion with PAWR) and to cell cycle (through HCF-1). THAP1 Structure The most extensive structural studies have focused on the THAP domain of THAP1 and implicitly serve as a model for the other hundred THAP proteins discovered in the animal kingdom. In many aspects, the THAP domain, which is about 80-amino-acid long, 2 HUMAN MUTATION, Vol. 00, No. 0, 1–12, 2011 Ta bl e 1. M ut at io ns R ep or te d in th e TH A P1 G en e in Ea ch Fa m il y/ P ro ba nd M u ta ti o n n am e N o m en cl at u re (p ro te in ) W il d T yp e co d o n M u ta n t co d o n M u ta ti o n al ev en t A m in o ac id p o si ti o n E xo n n o . St ru ct u re H C D R el at iv es in th e d at ab as e R ef er en ce n o . In sd el m u ta ti o n c. 1 3 5 _ 1 3 9 d el in sG G G T T T A ,a p .P h e4 5 L eu fs X 2 8 C C C In d el s In d el s st o p at 7 3 4 7 2 L 3 -T H A P d o m . 2 5 [F u ch s et al ., 2 0 0 9 ;B re ss m an et al ., 2 0 0 9 ] M is se n se m u ta ti o n s c. 1 A > G p .M et 1 ? A T G G T G A -> G 1 1 T H A P d o m . 1 [D e C ar va lh o A gu ia r et al ., 2 0 1 0 ] c. 1 6 T > C p .S er 6 P ro T C C C C C T -> C 6 1 L 1 -T H A P d o m . 1 [C lo t et al ., 2 0 1 0 ] c. 1 7 C > T p .S er 6 P h e T C C T T C C -> T 6 1 L 1 -T H A P d o m . 0 [H o u ld en et al ., 2 0 1 0 ] c. 2 3 A > G p .T yr 8 C ys A T G T G C A -> G 8 1 L 1 -T H A P d o m . 0 [H o u ld en et al ., 2 0 1 0 ] c. 2 5 G > T p .G ly 9 C ys G G C T G C G -> T 9 1 L 1 -T H A P d o m . 1 [X ia o et al ., 2 0 1 0 ] c. 3 6 C > A p .A sn 1 2 L ys A A C A A A C -> A 1 2 1 L 1 -T H A P d o m . 0 [B re ss m an et al ., 2 0 0 9 ] c. 3 8 G > A p .A rg 1 3 H is C G C C A C G -> A 1 3 1 L 1 -T H A P d o m . 1 [Z it te l et al ., 2 0 1 0 ] c. 5 0 A > G p .A sp 1 7 G ly G A C G G C A -> G 1 7 1 L 1 -T H A P d o m . 0 [X ia o et al ., 2 0 1 0 ] c. 6 1 T > A p .S er 2 1 T h r T C T A C T T -> A 2 1 1 L 1 -T H A P d o m . 0 [B re ss m an et al ., 2 0 0 9 ] c. 6 8 A > C p .H is 2 3 P ro C A C C C C A -> C 2 3 1 B S1 -T H A P d o m . 0 [K ai se r et al ., 2 0 1 0 ] c. 7 7 C > G p .P ro 2 6 A rg C C T C G T C -> G 2 6 2 L 2 -T H A P d o m . D N A b in d in g 0 [H o u ld en et al ., 2 0 1 0 ] c. 8 6 G > A p .A rg 2 9 G ln C G A C A A G -> A 2 9 2 L 2 -T H A P d o m . D N A b in d in g 0 [P ai sa n -R u iz et al ., 2 0 0 9 ] c. 8 6 G > A p .A rg 2 9 G ln C G A C A A G -> A 2 9 2 L 2 -T H A P d o m . D N A b in d in g 0 [P ai sa n -R u iz et al ., 2 0 0 9 ] c. 8 6 G > C p .A rg 2 9 P ro C G A C C A G -> C 2 9 2 L 2 -T H A P d o m . D N A b in d in g 0 [B re ss m an et al ., 2 0 0 9 ] c. 8 9 C > G p .P ro 3 0 A rg C C C C G C C -> G 3 0 2 L 2 -T H A P d o m . 1 [J ec h et al ., 2 0 1 1 ] c. 9 5 T > A p .L eu 3 2 H is C T T C A T T -> A 3 2 2 L 2 -T H A P d o m . A A in te ra ct io n s 2 [S ch n ei d er et al ., 2 0 1 1 ] c. 1 1 5 G > A p .A la 3 9 T h r G C T A C T G -> A 3 9 2 A H 1 -T H A P d o m . A A in te ra ct io n s 0 [B re ss m an et al ., 2 0 0 9 ] c. 1 6 1 G > A p .C ys 5 4 T yr T G T T A T G -> A 5 4 2 B S2 -T H A P d o m . C 2 C H m o ti f 0 [G av ar in i et al ., 2 0 1 0 ] c. 1 6 9 C > A p .H is 5 7 A sn C A C A A C C -> A 5 7 2 A H 2 -T H A P d o m . C 2 C H m o ti f 0 [S ö h n et al ., 2 0 1 0 ] c. 1 7 6 C > T p .T h r5 9 Il e A C T A T T C -> T 5 9 2 T H A P d o m . 0 [G ro en et al ., 2 0 1 0 ] c. 2 1 5 T > G p .L eu 7 2 A rg C T G C G G T -> G 7 2 2 L 4 -T H A P d o m . 0 [C lo t et al ., 2 0 1 0 ] c. 2 2 4 A > T p .A sn 7 5 Il e A A T A T T A -> T 7 5 2 L 4 -T H A P d o m . 1 [C h en g et al ., 2 0 1 0 ] c. 2 4 1 T > C p .P h e8 1 L eu T T T C T T T -> C 8 1 2 A H 4 -T H A P d o m . A A in te ra ct io n s 3 [F u ch s et al ., 2 0 0 9 ] c. 2 4 7 T > C p .C ys 8 3 A rg T G T C G T T -> C 8 3 2 0 [S ö h n et al ., 2 0 1 0 ] c. 2 6 6 A > G p .L ys 8 9 A rg A A G A G G A -> G 8 9 2 0 [B re ss m an et al ., 2 0 0 9 ] c. 3 9 5 T > C p .P h e1 3 2 Se r T T C T C C T -> C 1 3 2 3 0 [X ia o et al ., 2 0 1 0 ] c. 4 0 7 A > G p .A sn 1 3 6 Se r A A C A G C A -> G 1 3 6 3 H C F -1 b in d in g 0 [H o u ld en et al ., 2 0 1 0 ] c. 4 0 8 C > G p .A sn 1 3 6 L ys A A C A A G C -> G 1 3 6 3 H C F -1 b in d in g 1 [G ro en et al ., 2 0 1 0 ] c. 4 1 0 A > G p .T yr 1 3 7 C ys T A T T G T A -> G 1 3 7 3 H C F -1 b in d in g 0 [S ö h n et al ., 2 0 1 0 ] c. 4 2 7 A > G p .M et 1 4 3 V al A T G G T G A -> G 1 4 3 3 C o il ed -c o il d o m . 0 [S ö h n et al ., 2 0 1 0 ] c. 4 4 6 T > C b p .I le 1 4 9 T h r A T T A C T T -> C 1 4 9 3 N L S 1 [X ia o et al ., 2 0 1 0 ] c. 4 4 9 A > C p .H is 1 5 0 P ro C A T C C T A -> C 1 5 0 3 N L S 0 [C h en g et al ., 2 0 1 0 ] c. 4 9 6 G > A p .A la 1 6 6 T h r G C A A C A G -> A 1 6 6 3 C o il ed -c o il d o m . 0 [X ia o et al ., 2 0 1 0 ] c. 5 0 6 G > A p .A rg 1 6 9 G ln C G A C A A G -> A 1 6 9 3 C o il ed -c o il d o m . 0 [H o u ld en et al ., 2 0 1 0 ] c. 5 0 8 T > C p .C ys 1 7 0 A rg A T G C G C T -> C 1 7 0 3 C o il ed -c o il d o m . 0 [B o n et ti et al ., 2 0 0 9 ] c. 5 5 9 C > A p .G ln 1 8 7 L ys C A G A A G C -> A 1 8 7 3 C o il ed -c o il d o m . 0 [X ia o et al ., 2 0 1 0 ] c. 5 7 4 G > A p .A sp 1 9 2 A sn G A C A A C G -> A 1 9 2 3 0 [S ö h n et al ., 2 0 1 0 ] (C on ti n u ed ) HUMAN MUTATION, Vol. 00, No. 0, 1–12, 2011 3 Ta bl e 1. C on ti nu ed M u ta ti o n n am e N o m en cl at u re (p ro te in ) W il d T yp e co d o n M u ta n t co d o n M u ta ti o n al ev en t A m in o ac id p o si ti o n E xo n n o . St ru ct u re H C D R el at iv es in th e d at ab as e R ef er en ce n o . N o n se n se m u ta ti o n s c. 7 C > T p .G ln 3 X C A G T A G C -> T 3 1 L 1 -T H A P d o m . 0 [H o u ld en et al ., 2 0 1 0 ] c. 8 5 C > T p .A rg 2 9 X C G A T G A C -> T 2 9 2 L 2 -T H A P d o m . D N A b in d in g 0 [B re ss m an et al ., 2 0 0 9 ] c. 1 5 0 T > G p .T yr 5 0 X T A T T A G T -> G 5 0 2 L 3 -T H A P d o m . D N A b in d in g 0 [H o u ld en et al ., 2 0 1 0 ] c. 3 7 0 C > T p .G ln 1 2 4 X C A G T A G C -> T 1 2 4 3 0 [S ö h n et al ., 2 0 1 0 ] S m al lo u t- o f- fr am e d el et io n s c. 2 d el T p .M et 1 ? A T G d el 1 b St o p at 7 2 1 1 T H A P d o m . 0 [B re ss m an et al ., 2 0 0 9 ] c. 2 0 _ 3 3 d el p .A la 7 G lu fs X 2 3 G C C d el 1 4 b St o p at 2 9 7 1 L 1 -T H A P d o m . 0 [C lo t et al ., 2 0 1 0 ] c. 1 7 4 d el T p .P h e5 8 L eu fs X 1 5 T T T d el 1 c St o p at 7 2 5 8 2 T H A P d o m . D N A b in d in g 0 [H o u ld en et al ., 2 0 1 0 ] c. 1 9 7 _ 1 9 8 d el A G p .G lu 6 6 V al fs X 1 9 G A G d el 2 b St o p at 8 4 6 6 2 L 4 -T H A P d o m . 0 [S ö h n et al ., 2 0 1 0 ] c. 2 3 6 d el C p .T h r7 9 L ys fs X 4 1 A C A d el 1 b St o p at 1 1 9 7 9 2 A H 4 -T H A P d o m . A V P T IF m o ti f 0 [H o u ld en et al ., 2 0 1 0 ] c. 3 7 7 _ 3 7 8 d el C T p .P ro 1 2 6 A rg fs X 2 C C T d el 2 b St o p at 1 2 7 1 2 6 3 1 [B la n ch ar d et al ., 2 0 1 1 ] c. 3 8 8 _ 3 8 9 d el T C p .V al 1 3 1 P h ef sX 3 T C A d el 2 a St o p at 1 3 3 1 3 0 3 0 [S ö h n et al ., 2 0 1 0 ] c. 3 8 8 _ 3 8 9 d el T C b p .V al 1 3 1 P h ef sX 3 T C A d el 2 a St o p at 1 3 3 1 3 0 3 0 [D ja rm at i et al ., 2 0 1 0 ] c. 4 3 6 _ 4 4 3 d el p .A rg 1 4 6 A sp fs X 9 C G G d el 8 a St o p at 1 5 4 1 4 6 3 N L S 0 [C lo t et al ., 2 0 1 0 ] c. 4 6 0 d el C p .G ln 1 5 4 Se rf sX 2 7 C A G d el 1 a St o p at 1 8 0 1 5 4 3 N L S 0 [B re ss m an et al ., 2 0 0 9 ] c. 4 7 4 d el A b p .L ys 1 5 8 A sn fs X 2 3 A A A d el 1 c St o p at 1 8 0 1 5 8 3 N L S 2 [D ja rm at i et al ., 2 0 1 0 ] S m al li n -f ra m e d el et io n c. 2 0 7 _ 2 0 9 d el C A A p .A sn 6 9 _ A sn 6 9 d el A A C d el 3 c In -f ra m e d el 6 9 2 L 4 -T H A P d o m . 0 [G ro en et al ., 2 0 1 0 ] c. 2 0 7 _ 2 0 9 d el C A A p .A sn 6 9 _ A sn 6 9 d el A A C d el 3 c In -f ra m e d el 6 9 2 L 4 -T H A P d o m . 2 [C lo t et al ., 2 0 1 0 ] S m al lo u t- o f- fr am e in se rt io n c. 5 1 4 d u p p .A rg 1 7 2 L ys fs X 7 A G G in s1 b St o p at 1 7 8 1 7 2 3 C o il ed -c o il d o m . 0 [B la n ch ar d et al ., 2 0 1 1 ] F o r ea ch m u ta ti o n ,t h e n u cl eo ti d e an d p ro te in n am es ar e gi ve n ac co rd in g to th e n o m en cl at u re gu id el in es (h tt p :/ /w w w .h gv s. o rg /) . a N o te th at th e c. 1 3 4 _ 1 3 5 in sG G G T T ;1 3 7 _ 1 3 9 d el A A C A m is h m u ta ti o n h as b ee n re n am ed in c. 1 3 5 _ 1 3 9 d el in sG G G T T T A to fo ll o w th e gu id el in e fo r “i n sd el ” m u ta ti o n (h tt p :/ /w w w .h gv s. o rg /) . b Sa m e p at ie n t d es cr ib ed in th e tw o gi ve n re fe re n ce s; “T H A P d o m .” :T H A P d o m ai n ; “L 1 -T H A P d o m .” :T H A P -d o m ai n lo o p 1 ; “L 2 -T H A P d o m .” :T H A P -d o m ai n lo o p 2 ; “L 3 -T H A P d o m .” :T H A P -d o m ai n lo o p 3 ; “L 4 -T H A P d o m .” :T H A P -d o m ai n lo o p 4 ; “A H 1 -T H A P d o m .” :T H A P -d o m ai n α -h el ix 1 ; “H 2 -T H A P d o m .” :T H A P -d o m ai n 3 1 0 h el ix 2 ; “H 3 -T H A P d o m .” :T H A P -d o m ai n 3 1 0 h el ix 3 ; “H 4 -T H A P d o m .” :T H A P -d o m ai n 3 1 0 h el ix 4 ; “C o il ed -c o il d o m .” :c o il ed -c o il d o m ai n ; “N L S” :n u cl ea r lo ca li za ti o n se q u en ce ; “D N A b in d in g -” :a m in o ac id s p ro ve n to b e in vo lv ed in D N A -b in d in g ac ti vi ty ; “C 2 C H m o ti f” :a m in o ac id s in vo lv ed in th e zi n c- d ep en d en t D N A -b in d in g ac ti vi ty o f T H A P 1 th ro u gh th e C 2 C H m o ti f; “H C F -1 b in d in g” :a m in o ac id s in vo lv ed in th e in te ra ct io n w it h H C F -1 ; “A A in te ra ct io n s” :i n te ra ct io n s b et w ee n T H A P 1 am in o ac id s as d es cr ib ed in F ig u re 1 B .E xa m p le :T H A P 1 am in o ac id 3 9 is in vo lv ed in in te ra ct io n s w it h T H A P 1 am in o ac id s 4 0 ,4 5 ,a n d 8 1 . F o r ea ch p ro b an d ,t h e n u m b er o f re la ti ve s in cl u d ed in th e d at ab as e is gi ve n . 4 HUMAN MUTATION, Vol. 00, No. 0, 1–12, 2011 1 96 108 146 162 21381Amino acid: 1 72 267 268 64271 Exon 3 Nucleotide: Exon 1 Exon 2 139 190 THAP domain Proline-rich region Nuclear localization sequence Coiled-coil domain L1 L2 L3 L4B1 B2H1 H2 H3 H4 B1: MVQSCSAYGCKNRYDKDKPVSFHKFPL RPSLCKEWEAAVRRKNFKPTKYSSICSEHFTPDCFKRECNNKLLKENAVPTIFT C C K P R W R F T C HF P K L E V Y 1 5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 No binding activity: Binding activity diminished: Known amino acid interactions: 36 2526 A: THAP domain 1 5424 267 Amino acid: THAP1b THAP1a Amino acid: S S Y K PSLCKE R KN KP K SS ST No change in activity when mutated to Ala: Loop Beta strand Alpha helice THAP1 B2: 29 32 39 40 58 45 63 78 76 81 Figure 1. THAP1 and its THAP domain. A: Schematic representation of the THAP1 gene and its two isoforms. Isoform 1 (THAP1a) and isoform 2 (THAP1b) result from alternative splicing of exon 2; B1: Schematic representation of the THAP domain with amino acid positions for each structure. The upper figure represents the α-helix H1, the three 310 helices (H2, H3, and H4) (all in orange), the two β-sheets (B1 and B2) (in red), and the four loops (L1 to L4, in gray). “No binding activity”: when these amino acids of the THAP domain are mutated to alanine, THAP1 shows loss of zinc-dependent DNA-binding activity; amino acids that are part of the C2CH motif are in red; “Binding activity diminished”: when these amino acids are mutated to alanine, there is partial loss of DNA-binding activity; “No binding activity modification when mutated to Ala”: mutation of these amino acids to alanine does not change the DNA-binding activity of THAP1 [Bessiere et al., 2008; Clouaire et al., 2005]; “Known amino acid interactions”: structural representation of the known amino acid interactions [Bessiere et al., 2008]; B2: Topology diagram of the THAP domain of human THAP1 showing the secondary structure elements: the α-helix H1, the three 310 helices (H2, H3, and H4) (all in orange), the two β-sheets (B1 and B2) (in red), and the location of the four loops (L1 to L4). The zinc finger and the four ligands are shown in gray. Modified from Bessiere et al. [2008]. exhibits some specific features (Fig. 1B1). It is characterized by a C2CH signature (Cys-X2–4-Cys-X35–53-Cys-X2-His) associated with four invariant residues (Pro26, Trp36, Phe58, Pro78) in the THAP1 sequence. Direct mutagenesis of each of these eight amino acids showed their critical identity for the zinc-dependent specific binding of the THAP domain to a precise DNA sequence (TXXXGGCA: THABS consensus sequence) [Clouaire et al., 2005] (Fig. 1B2). An AVPTIF box (Ala76-Phe81, in the THAP1 sequence), also essential for DNA binding, is located at the C-terminus of the THAP domain. The NMR solution structure of the THAP zinc finger of THAP1 [Bessiere et al., 2008] indicates that this domain consists of a first loop (L1: Gln3-Ser21), followed by two short antiparallel β-sheets (B1: Phe22-Lys24 and B2: Ser52-Cys54) separated by a loop (L2: Phe25-Lys32)-helix (H1: Cys33-Val40)-loop (L3: Arg41-Ser51) struc- ture (Fig. 1B2). The four residues of the C2CH motif (Cys5-Cys10- Cys54-H57) are three-dimensionally brought together to form the zinc-dependent DNA-binding site. After the second β-strand, there are three 310helices: H2 (Ser55-His57), H3 (Arg60-Phe63), and H4 (Thr79-Phe81). H4 contains part of the AVPTIF motif and is pre- ceded by a last mobile loop (L4: Lys64-Pro78). The alpha-helix H1 can be considered as an important inter- action hub and it is essential for the maintenance of the three- dimensional (3-D) shape of the THAP domain (Fig. 1B1). Indeed, it makes contacts with many residues distributed throughout the do- main allowing a spatial connection between H1 (Trp36, Ala39, Val40), L2 (Phe25, Pro26, Arg29), L3 (Phe45), ∼H2 (Phe58), H3 (Phe63), L4 (Pro78), and H4 (Phe81). Residue Trp36, which is located in the core of H1, seems to be a key amino acid for this function. It makes hydrophobic contacts with four amino acids (Phe25, Arg29, Phe58, and Phe63), and NMR data strongly suggest the presence of interac- tion also with the two highly conserved prolines (Pro26 and Pro78). Four of these residues have been functionally tested (Pro26, Arg29, Phe58, Pro78) and their replacement by alanine strongly decreases the THAP domain ability to bind to DNA, suggesting that they have a critical role in the maintenance of THAP1 3-D shape [Clouaire et al., 2005]. The 3-D structure of the DNA–THAP domain complex has re- cently been solved [Campagne et al., 2010], showing that the THAP domain contacts DNA by insertion of its two antiparallel β-sheets into the major DNA groove. This ensures the connection of the second half of loop L3 (residues Lys46 to Ser51) with the core of the THABS sequence. Upon binding to DNA, the folding of THAP1 is changed, enabling loop L4 to contact the DNA minor groove. Based on their nonconservation in the THAP family, specific base-contact with the core of the THABS motif and chemical behavior with non- specific DNA, some THAP domain residues have been proposed to play a crucial role in the specificity of DNA recognition: Tyr50 and Ser51 of loop L3, Ser52 of B2, Gln3 and Ser4 of the N-terminal tail and, finally Arg65 of L4. Other residues are involved in the nonspecific interaction with DNA, but their role remains absolutely necessary for the targeting of THAP1 to the THABS sequence. For example, Thr48, which is located in the middle of the DNA-interaction area of loop L3, shows a chemical shift perturbation when the THAP domain is incubated with an unrelated DNA sequence, suggesting that Thr48 is involved in nonspecific interactions. However, muta- tion of this residue led to total loss of the DNA-binding activity of the THAP domain and it has been proposed that, with Lys46, it may be important for positioning the protein onto DNA. HUMAN MUTATION, Vol. 00, No. 0, 1–12, 2011 5 Table 2. Clinical Characteristics of Carriers of THAP1 Mutations Amish patients a N = 25 Non-Amish patients N = 81 All patients b N = 106 Sex: F 15 (60%) F 46c (57%) F 61c (58%) Age at onset (y): Median 14.5 (5–38) Median 13 (2–62) Median 13 (2–62) Age at last examination (y): Median 40 (10–66) Median 44 (7–86) Median 44 (7–86) Family history: 100% 58.5% 61.4% Site at onset: 25 patients 77 patientsd 102 patientsd Upper limb 11 (44%)e 36 (47%)e 48 (47%) Lower limb 1 (4%)e 14 (18%)e 17 (17%) Cervical 5 (20%)e 17 (22%)e 23 (23%) Cranial 8 (32%) 17 (22%) 25 (25%) Site at examination: 25 patients 81 patients 106 patients Upper limb 22 (88%) 60 (74%) 82 (77%) Lower limb 12 (48%) 39 (48%) 51 (48%) Cervical 14 (56%) 59 (73%) 73 (69%) Cranial 18 (72%)e 40 (49%)e 58 (55%) Speech 17 (68%)e 48 (59%)e,f 65 (61%)e,f Distribution: 25 patients 81 patients 106 patients FD: 3 (12%) FD: 16 (20%) FD: 19 (18%) SD: 10 (40%) SD: 20 (25%) SD: 30 (28%) MD: 4 (16%) MD: 5 (6%) MD: 9 (8%) GD: 8 (32%) GD: 40 (49%) GD: 48 (45%) Gender: F, female; M, male. Dystonia distribution: “GD”: generalized dystonia; “FD”: focal dystonia; “MD”: multifocal dystonia; “SD”: segmental dystonia. a Four families with common ancestry, 25 patients. Data from Fuchs et al. [Söhn et al., 2010]. b Thirty-five families, 22 sporadic, total 106 patients. c Sex is unknown for one patient from Söhn et al. [Söhn et al., 2010]. d For some patients, data were unavailable [Bressman et al., 2009; Gavarini et al., 2010; Kaiser et al., 2010]. e Data were estimated from Bressman et al. [2009] and Fuchs et al. [2009] when possible. f Of note, two patients from Groen et al. [2010] have oromandibular dystonia without knowing if there was speech disturbance or not. THAP1 Is Involved in DYT6 Dystonia DYT6 dystonia has been first described in the genetically iso- lated Amish–Mennonite population [Almasy et al., 1997] as an autosomal dominant trait with incomplete penetrance. In two Amish–Mennonite families, the site of onset of the disease was lo- calized with similar frequency to upper limb, cervical, or cranial muscles and symptoms began between 5 and 38 years of age (aver- age 18.1) with a tendency to progressive involvement of other body regions. Hence, DYT6 dystonia was named “idiopathic torsion dys- tonia of ‘mixed’ type.” These families were large enough to identify by linkage analysis the morbid locus, which maps in both families to chromosome 8p21-22. The locus was further narrowed by the discovery of other affected Amish–Mennonite families [Saunders-Pullman et al., 2007] and, in 2009, Fuchs et al. [2009] demonstrated the involvement of a THAP1 mutation in DYT6 dystonia: an insertion/deletion founder muta- tion found in four related Amish–Mennonite families. A point mu- tation was then identified in an unrelated non-Amish–Mennonite German family (family S, c.214T>C), indicating that DYT6 dysto- nia was present also in populations of different ancestry. Indeed, the complex insertion/deletion mutation was later identified also in a German family with no known Amish–Mennonites origin but sharing the same haplotype [Bressman et al., 2009]. Several teams have reported to date 53 different mutations in the coding sequence of THAP1 in 56 families (56 probands and 43 rela- tives) of American [Bressman et al., 2009; Xiao et al., 2010], Chinese [Cheng et al., 2010; Xiao et al., 2010], Czech [Jech et al., 2011], Dutch [Groen et al., 2010], English [Houlden et al., 2010], French [Blanchard et al., 2011; Clot et al., 2010; Houlden et al., 2010], German [Bressman et al., 2009; Djarmati et al., 2009; Fuchs et al., 2009; Houlden et al., 2010; Söhn et al., 2010; Zittel et al., 2010], Greek [Bonetti et al., 2009; Houlden et al., 2010], Iranian [Schneider et al., 2011], Jewish [Houlden et al., 2010], Mauritian/Indian [Houlden et al., 2010], Spanish [Paisan-Ruiz et al., 2009], and Brazilian [De Carvalho Aguiar et al., 2010] ancestry. Unfortunately, clinical data are not available for 11 probands implemented in the database [Bressman et al., 2009; Gavarini et al., 2010; Kaiser et al., 2010] and insufficient clinical and molecular data for 10 relatives do not allow to implement them in the database [Bressman et al., 2009]. One of the mutation has been described in an asymptomatic individual control (c.197_198delAG) [Söhn et al., 2010]. The enlargement of the group of patients with THAP1 muta- tions (108 up to now without the asymptomatic carrier) has pro- vided new data for delineating the phenotype of DYT6 dystonia, which nevertheless has hardly changed since the original description (Table 2). DYT6 dystonia has been reported slightly more in females than males. The age distribution at onset is broad, ranging from 2 to 62 years (mean 17.8 years [calculated from 78 patients with available data]). Onset before 10 years of age occurs in 35.8% of the patients, and onset after 28 years of age in one-sixth (16%) of the patients. In almost half (47%) of the patients, the first symptoms of dystonia concern an upper limb, while cranial or cervical onset is almost equally observed in the other half of the patients. Onset in a lower limb is less common (17/102 patients reported in the litera- ture). DYT6 dystonia tends to spread slowly to other body regions. Although the symptom distribution varies widely, the involvement of the cranial region is frequent, and disability mainly results from cranial dystonia with speech difficulties. At last follow-up, almost half of the reported patients with DYT6 dystonia suffered from gen- eralized dystonia, but most of them remained ambulatory. To date, 9 patients underwent deep brain stimulation; one good response is reported in one patient [Jech et al., 2011], although deep brain stimulation (DBS) resulted in only mild or moderate improvement 6 HUMAN MUTATION, Vol. 00, No. 0, 1–12, 2011 in the other patients, with poor effect on speech or dysarthria [Clot et al., 2010; Djarmati et al., 2009; Groen et al., 2010; Jech et al., 2011; Zittel et al., 2010]. Early onset in a limb and progressive generalization in few years are shared features between DYT6 and DYT1 dystonias [Valente and Albanese, 2010]. Nevertheless, several differences emerge. In DYT6 dystonia, cervical or cranial muscles are involved most of the time, whereas lower limbs are frequently spared. The DYT6 clinical spectrum is more variable, and some adult-onset cases with stable focal dystonia are reported as well. Onset (mean age) is somewhat later and with a wider age range in DYT6 than in DYT1 patients. The penetrance in DYT6 is difficult to evaluate; it was performed only in the Amish–Mennonite population (around 60% [Saunders- Pullman et al., 2007]). Recently, TOR1A has been demonstrated as a direct target of THAP1 by EMSA, ChIP, and luciferase reporter gene assays [Gavarini et al., 2010; Kaiser et al., 2010]. Specific modulation of torsinA expression was not reported in nonneuronal cells after THAP1 knockdown or overexpression, nor in fibroblasts or lym- phoblast cells from DYT6 patients [Gavarini et al., 2010; Kaiser et al., 2010]. The lack of effect could be explained by tissue and/or developmental stage specificity. However, THAP1 is predicted to po- tentially regulate numerous gene targets [data not published] and TOR1A deregulation in DYT6 patients should then be responsible of only one part of the phenotype. Nevertheless, human wild-type torsinA can lead to a cellular dysfunction of mice neurons when overexpressed at high levels (abundant inclusion-like structures, ab- normalities of the NE ultrastructure, behavioral abnormalities, sig- nificant reduction of dopamine levels, serotonin, and homoranillic acid) as shown in transgenic mice line hWT24 [Grundmann et al., 2007]. The UMD-THAP1 Database We have used the Universal Mutation Database (UMD) software [Beroud et al., 2000, 2005] to create a computerized locus-specific database that contains information about the published mutations of the THAP1 gene. Codons were numbered (e.g., +1 = A of ATG) based on the cDNA sequence of the THAP1 variant 1 (isoform with three exons) (accession number NM_018105.2, “Homo sapiens THAP domain containing, apoptosis associated protein 1 (THAP1), transcript variant 1, mRNA”) obtained from the GenBank database. Intron–exon boundaries as well as intronic sequences were defined by matching the cDNA sequence to the corresponding genome sequence (NT_167187.1) and to the module organization from SwissProt (accession number Q9NVV9). The database follows the guidelines on mutation databases of the Hugo Mutation Database Initiative including the latest nomenclature (http://www.hgvs.org/). The THAP1 mutation database can list point mutations, large and small deletions, insertions, indels, and mutations affecting splicing (intronic mutations) in the THAP1 gene. It cannot accommodate mutations from the UTR and promotor regions. In addition, two mutations that affect the same allele are entered as two different records linked by the same sample ID. It allows to link specific poly- morphisms identified in the course of sequencing to one mutation for a specific patient. Separate entries are made for each probands and for each relative (with different sample ID) in order to de- scribe specific phenotypes. If the same mutation in a single patient has been reported in different papers (example: the same patients 7021, L-2257, L-3736, and L-3637 were first reported in [Djarmati et al., 2009] and subsequently in [Zittel et al., 2010]), only one entry was made. If the same mutation has been reported in apparently unrelated patients (example c.207_209delCAA identified in French and Dutch patients [Clot et al., 2010; Groen et al., 2010]), separate entries were made for each patient as recurrent mutations, in the absence of haplotypes demonstrating common ancestor. For each mutation, different types of information are provided: genetic (exon and codon number, wild type and mutant codon, mutational event, mutation name), protein (wild type and mutant amino acids, af- fected domain, mutation name), clinical data (age of onset, site of onset, dystonia distribution, sites involved, genealogic tree when available), and experimental data. We have annotated the THAP1 sequence with indirect arguments in order to determine whether potential missense mutations are really causative using the UMD-Predictor R© tool [Frederic et al., 2009]. The UMD-Predictor R© tool is dedicated to the analysis of nucleotide substitutions in cDNA sequences. It provides a com- binatorial approach that includes localization of the mutation in the protein, conservation, biochemical properties of the mutant and wild type residue, and potential impact of the variation on mRNA (http://www.umd.be). The indirect arguments, listed in the “Highly conserved domain” (HCD) of UMD-Predictor R© are: (1) amino acids involved in the zinc-dependent DNA-binding activ- ity of THAP1 through the C2CH motif [Clouaire et al., 2005]; (2) amino acids involved in DNA binding of the THAP domain of hu- man THAP1 [Bessiere et al., 2008; Clouaire et al., 2005]; (3) amino acids involved in the interaction with HCF-1 (HBM consensus se- quence) [Mazars et al., 2010]; and (4) THAP1 amino acids involved in interactions with other THAP1 amino acids [Bessiere et al., 2008] (Fig. 1B1). Mutations that potentially affect splicing can be highlighted with dedicated algorithms that have been recently implemented in the Human Splicing Finder (HSF) tool [Desmet et al., 2009]. HSF is also available for any gene as a web resource (http://www.umd.be/HSF/) and is dedicated to the prediction of the impact of a mutation on splicing signals and/or to the identification of these motifs in any human sequence. It contains all available matrices for auxiliary sequences prediction, as well as new position weight matrices to evaluate the strength of 5′ and 3′ splice sites and of branch point sequences. The UMD-THAP1 Database can be used with different routines described at http://www.umd.be/THAP1/. It is possible to analyze the mutation distribution by exon or by mutation type and the mechanism for deletion/insertion. With the “HCD” annotation, it is possible to know whether an amino acid has a known reported function in the protein. The UMD software has already been used successfully since 1994 for the analysis of multiple genes implicated in genetic diseases, such as LDLR [Villeger et al., 2002], TGFBR2 [Frederic et al., 2008], DYS [Krahn et al., 2009], and in cancer, such as TP53 [Cariello et al., 1994a, b; Hamroun et al., 2006] or MEN1 [Wautot et al., 2002]. Twenty-nine UMDs are available through the web. More information concerning the UMD software is available at http://www.umd.be. The current database and the ensuing updated versions are avail- able online at http://www.umd.be. Notifications of omissions and errors in the current version as well as specific phenotypic data will be gratefully received by the corresponding author. The software will be expanded as the database grows and according to the users’ requirements, and new functions may be implemented. Mutation Analysis To date, a total of 53 different THAP1 mutations are known (Table 1) since two of the three patients and their specific mutation reported by Zittel et al. [2010] were already described by Djarmati HUMAN MUTATION, Vol. 00, No. 0, 1–12, 2011 7 et al. [2009], and the patient reported by Van Gerpen et al. [2010] is part of the cohort described by Xiao et al. [2010]. These mutations are mainly private, missense, and widely distributed throughout the gene. Only three apparently recurrent mutations have been de- scribed. For the c.86G>A, described in two patients by Paisan-Ruiz et al. [2009], no other polymorphism has been described ruling out common ancestor [Coro Paisan-Ruiz, personal communication]. The c.388_389delTC mutation [Djarmati et al., 2009; Söhn et al., 2010] and the c.207_209delCAA mutation [Clot et al., 2010; Groen et al., 2010] have been characterized by different teams, but until haplotype analyses become available, it is unclear whether these are truly recurrent mutations or whether they result from a founder effect. These mutations have then been considered as recurrent. It is noteworthy that two mutations, c.407A>G [Houlden et al., 2010] and c.95T>A [Schneider et al., 2011] have been found to be ho- mozygous. Insdel Mutations The heterozygous Amish mutation, described as a 5-base (GGGTT) insertion followed by a 3-base deletion (AAC) (c.134_135insGGGTT;137_139delAAC) [Fuchs et al., 2009], has been detected in four related Amish–Mennonite families (fami- lies M, C, R, and W). This mutation was later identified also in a German family without known Amish–Mennonite ancestry (fam- ily 2) but sharing the Amish haplotype [Bressman et al., 2009]. This patient has then been implemented as a relative. To follow the international guideline for “insdel” mutation nomenclature (http://www.hgvs.org/), the name of the mutation has been cor- rected to c.135_139delinsGGGTTTA, that is deletion of five bases (TAAAC) and insertion of seven bases (GGGTTTA), which results in the same predicted mutant protein. Nonsense and Missense Mutation Four reported mutations are nonsense mutation: c.7C>T [Houlden et al., 2010], c.85C>T [Bressman et al., 2009], c.150T>G [Houlden et al., 2010], and c.370C>T [Söhn et al., 2010]. Thirty- seven are missense mutations. Five mutations are predicted to im- pair THAP1 DNA-binding activity: c.77C>G [Houlden et al., 2010]; the two recurrent mutations c.86G>A [Paisan-Ruiz et al., 2009]; c.86G>C [Bressman et al., 2009]; and c.169C>A in the C2CH motif [Söhn et al., 2010]. Three mutations are located in the sequence involved in THAP1 binding to HCF-1: c.407A>G [Houlden et al., 2010], c.408C>G [Groen et al., 2010], and c.410A>G [Söhn et al., 2010]. One mutation (c.241T>C) is located in the AVPTIF motif [Fuchs et al., 2009]. The homozygous mutation c.407A>G [Houlden et al., 2010] creates a potential acceptor splice site in exon 3, but its strength and the presence of many upstream good candidate accep- tor splice sites make it unlikely that it will be used. The c.266A>G mutation [Bressman et al., 2009] is located in the penultimate nucleotide of exon 2. This mutation does not impair the donor acceptor splice site, neither auxiliary splicing sequences (exonic sequence enhancer [ESE] or silencer [ESS]). For the 24 remaining mutations, the mutated codon does not seem to be clearly involved in a particular function. Among all the missense mutations, only nine were not predicted by the UMD- Predictor R© tool to be pathogenic based on a combinatorial ap- proach [Frederic et al., 2009] (Table 3): c.38G>A [Zittel et al., 2010], c.61T>A [Bressman et al., 2009], c.169C>A [Söhn et al., 2010], c.241T>C [Fuchs et al., 2009], c.266A>G [Bressman et al., 2009], c.427A>G [Söhn et al., 2010], c.506G>A [Houlden et al., 2010], c.559C>A [Xiao et al., 2010], and c.574G>A [Söhn et al., 2010]. Nevertheless, mutation c.241T>C, reported in dbSNP without frequency data, shows a decreased DNA-binding activity in vitro [Fuchs et al., 2009]. EMSA experiments for c.61T>A show that the mutation abolishes binding of THAP1 to TOR1A [Bressman et al., 2009; Gavarini et al., 2010], while mutation c.38G>A is shown to significantly reduce THAP1-mediated repression of TOR1A (57% of residual activity)[Kaiser et al., 2010; Zittel et al., 2010]. UMD-Predictor R© tool has a very high positive predictive value and a high negative predictive value. Mutations experimentally demonstrated by gel-shift assays (Table 1) [Bessiere et al., 2008; Clouaire et al., 2005] to have (1) no binding activity (p.Cys5Ala, p.Cys10Ala, p.Lys24Ala, p.Pro26Ala, p.Arg29Ala, p.Trp36Ala, p.Arg42Ala, p.Phe45Ala, p.Thr48Ala, p.Cys54Ala, p.His57Ala, p.Phe58Ala, and p.Pro78Ala) are all predicted as pathogenic; (2) a diminished binding activity were also all predicted as pathogenic (p.Glu37Ala, p.Lys11Ala, p.Leu27Ala, and p.Tyr50Ala) or proba- bly pathogenic (p.Val40Ala); (3) no change in activity were ei- ther predicted: as pathogenic (p.Ser4Ala, p.Tyr8Ala, p.Lys16Ala, p.Leu32Ala, p.Arg41Ala, p.Asn44Ala, p.Lys46Ala, p.Lys49Ala); as probably pathogenic (p.Thr28Ala, p.Pro30Ala, p.Glu35Ala, p.Lys34Ala, p.Lys43Ala); as polymorphism (p.Ser55Ala); or as prob- able polymorphism (p.Ser6Ala, p.Pro47Ala, p.Ser31Ala, p.Cys33Ala, p.Ser51Ala, p.Ser52Ala). Nevertheless, these gel-shift assays demon- strate no modification of the DNA-binding activity “within the lim- its of detection of the present assay” and do not exclude “a role in the binding affinity and selectivity” [Bessiere et al., 2008]. Amino acids Ser4, Ser51, and Ser52 are described to be involved in the specificity of DNA recognition and Lys46 proposed, with Trp48, to be impor- tant for positioning the protein onto DNA [Campagne et al., 2010]. Therefore, the predictions are in full agreement with in vitro stud- ies for pathogenic mutations (18/18). For the alanine substitutions classified as “non modifying binding activities,” some discrepancies are observed but their interpretation is difficult as they could result from wrong predictions by UMD-predictor or wrong interpretation of experimental data. These experimental results show that in spite of the high efficiency of the UMD-predictor tool [Frederic et al., 2009] and other predic- tion tools [SIFT, PolyPhen, data not shown] for several genes (FBN1, FBN2, TGFBR1, TGFBR2 [Frederic et al., 2009], DYS [Krahn et al., 2009], CDKN2A [Kannengiesser et al. 2009], and ALAS2 [Ducamp et al. 2011]), more functional data will have to be collected to im- prove the prediction of missense variations in THAP1 gene. Insertion/Duplication Only one reported mutation is an insertion that corresponds to the duplication of a single base (A514). This insertion is predicted to create a premature termination codon (PTC), codon seven amino acids downstream of the mutation (position 178) [Blanchard et al., 2011]. Deletions Among the 13 small deletions, 11 create a PTC and two are in-frame deletions resulting in loss of the ASN69 residues (c.207_209delCAA). Two single-base deletions, c.174delT [Houlden et al., 2010] and c.474delA [Djarmati et al., 2009], might be the result of slipped mispairing. Five mutations are deletions of a re- peated sequence: c.197_198delAG [Söhn et al., 2010], the recurrent mutations c.388_389delTC [Djarmati et al., 2009; Söhn et al., 2010] and c.207_209delCAA [Clot et al., 2010; Groen et al., 2010]. One mutation, c.436_443del, is flanked by direct repeats [Clot et al., 2010]. For the five other mutations, the mechanisms are unknown 8 HUMAN MUTATION, Vol. 00, No. 0, 1–12, 2011 Ta bl e 3. A na ly si s of M is se ns e M ut at io ns w it h th e U M D -P re di ct or R© To ol N o m en cl at u re c. N o m en cl at u re p . St ru ct u re H C D C o n se rv at io n SI F T p ro b ab il it y B L O SU M 6 2 B io ch em ic al va lu es E SE m o d ifi ca ti o n Sp li ce si te P at h o ge n ic it y C o n cl u si o n c. 1 7 C > T p .S er 6 P h e L 1 -T H A P d o m . 0 .7 1 0 .0 6 -2 .0 0 0 .3 3 N o im p ac t 8 2 P at h o ge n ic c. 2 3 A > G p .T yr 8 C ys L 1 -T H A P d o m . 0 .7 1 0 .0 8 -2 .0 0 0 .3 8 PA SS C [7 3 .1 4 ] 7 6 P at h o ge n ic c. 2 5 G > T p .G ly 9 C ys L 1 -T H A P d o m . 0 .7 1 0 .0 2 -3 .0 0 0 .4 6 N o im p ac t 9 3 P at h o ge n ic c. 5 0 A > G p .A sp 1 7 G ly L 1 -T H A P d o m . 0 .7 1 0 .0 5 -1 .0 0 0 .2 9 N o im p ac t 7 9 P at h o ge n ic c. 6 8 A > C p .H is 2 3 P ro B S1 -T H A P d o m . 0 .5 0 .0 0 -2 .0 0 0 .2 5 N o im p ac t 8 8 P at h o ge n ic c. 7 7 C > G p .P ro 2 6 A rg L 2 -T H A P d o m . D N A b in d in g - 0 .7 9 0 .0 2 -2 .0 0 0 .1 7 -S C 3 5 [2 .7 3 ] N o im p ac t 1 0 0 P at h o ge n ic c. 8 6 G > C p .A rg 2 9 P ro L 2 -T H A P d o m . D N A b in d in g - 0 .7 9 0 .0 4 -2 .0 0 0 .1 7 PA SS C [7 5 .5 ] 1 0 0 P at h o ge n ic c. 8 9 C > G p .P ro 3 0 A rg L 2 -T H A P d o m 0 .7 9 0 .0 2 -2 .0 0 0 .1 7 N o im p ac t 8 7 P at h o ge n ic c. 9 5 T > A p .L eu 3 2 H is L 2 -T H A P d o m L 2 -T H A P d o m 0 .7 9 0 .0 1 -3 .0 0 0 .3 3 N o im p ac t 1 0 0 P at h o ge n ic c. 1 1 5 G > A p .A la 3 9 T h r A H 1 -T H A P d o m . L 2 -T H A P d o m 0 .7 9 0 .1 9 0 .0 0 0 .2 9 -S F 2 /A SF [2 .4 8 ] N o im p ac t 8 8 P at h o ge n ic c. 1 6 1 G > A a p .C ys 5 4 T yr B S2 -T H A P d o m C 2 C H m o ti f. 0 .7 9 0 .0 1 -2 .0 0 0 .3 8 N o im p ac t 9 9 P at h o ge n ic c. 2 1 5 T > G p .L eu 7 2 A rg L 4 -T H A P d o m . 0 .7 9 0 .1 3 -2 .0 0 0 .2 5 N o im p ac t 7 6 P at h o ge n ic c. 2 2 4 A > T p .A sn 7 5 Il e L 4 -T H A P d o m . 0 .7 1 0 .0 2 -3 .0 0 0 .3 3 N o im p ac t 9 3 P at h o ge n ic c. 2 4 7 T > C p .C ys 8 3 A rg 0 .8 6 0 .3 5 -3 .0 0 0 .2 5 N o im p ac t 8 2 P at h o ge n ic c. 3 9 5 T > C p .P h e1 3 2 Se r 1 0 .7 7 -2 .0 0 0 .3 3 N o im p ac t 7 6 P at h o ge n ic c. 4 0 8 C > G p .A sn 1 3 6 L ys H C F -1 b in d in g 1 0 .0 6 0 .0 0 0 .4 6 N o im p ac t 7 6 P at h o ge n ic c. 4 1 0 A > G p .T yr 1 3 7 C ys H C F -1 b in d in g 1 0 .1 0 -2 .0 0 0 .3 8 N o im p ac t 8 8 P at h o ge n ic c. 4 4 6 T > C p .I le 1 4 9 T h r N L S 0 .9 3 0 .0 2 -3 .0 0 0 .4 2 N o im p ac t 8 7 P at h o ge n ic c. 4 4 9 A > C p .H is 1 5 0 P ro N L S 0 .9 3 0 .2 4 -2 .0 0 0 .2 5 N o im p ac t 7 6 P at h o ge n ic c. 5 0 8 T > C p .C ys 1 7 0 A rg C o il ed -c o il d o m . 0 .8 6 0 .0 0 -3 .0 0 0 .2 5 N o im p ac t 9 4 P at h o ge n ic c. 1 A > G p .M et 1 ? T H A P d o m . 0 .7 1 0 .0 0 1 .0 0 0 .7 1 N o im p ac t 6 5 P ro b ab ly p at h o ge n ic c. 1 6 T > C p .S er 6 P ro L 1 -T H A P d o m . 0 .7 1 0 .2 0 -1 .0 0 0 .3 8 N o im p ac t 7 1 P ro b ab ly p at h o ge n ic c. 3 6 C > A p .A sn 1 2 L ys L 1 -T H A P d o m . 0 .7 1 0 .2 2 0 .0 0 0 .4 6 N o im p ac t 6 5 P ro b ab ly p at h o ge n ic c. 8 6 G > A p .A rg 2 9 G ln L 2 -T H A P d o m . D N A b in d in g - 0 .7 9 0 .0 6 1 .0 0 0 .5 0 PA SS C [7 1 .6 3 ] 6 5 P ro b ab ly p at h o ge n ic c. 1 7 6 C > T p .T h r5 9 Il e T H A P d o m . 0 .7 9 0 .0 7 -1 .0 0 0 .4 2 N o im p ac t 6 5 P ro b ab ly p at h o ge n ic c. 4 0 7 A > G p .A sn 1 3 6 Se r H C F -1 b in d in g 1 0 .1 3 1 .0 0 0 .6 3 P A S S C [8 5 .8 8 ] 7 1 P ro b ab ly p at h o ge n ic c. 4 9 6 G > A p .A la 1 6 6 T h r C o il ed -c o il d o m . 0 .9 3 0 .4 2 0 .0 0 0 .2 9 N o im p ac t 6 5 P ro b ab ly p at h o ge n ic c. 3 8 G > A a p .A rg 1 3 H is L 1 -T H A P d o m . 0 .7 1 0 .0 1 0 .0 0 0 .5 8 N o im p ac t 6 4 P ro b ab le p o ly m o rp h is m c. 1 6 9 C > A p .H is 5 7 A sn A H 2 -T H A P d o m . C 2 C H m o ti f 0 .7 9 0 .0 5 1 .0 0 0 .5 8 N o im p ac t 5 9 P ro b ab le p o ly m o rp h is m c. 2 4 1 T > C a p .P h e8 1 L eu A H 4 -T H A P d o m . A A in te ra ct io n s 0 .7 9 0 .4 5 0 .0 0 0 .7 1 N o im p ac t 5 9 P ro b ab le p o ly m o rp h is m c. 6 1 T > A a p .S er 2 1 T h r L 1 -T H A P d o m . 0 .7 1 0 .0 0 1 .0 0 0 .7 9 N o im p ac t 4 7 P o ly m o rp h is m c. 2 6 6 A > G p .L ys 8 9 A rg 0 .7 9 0 .3 4 2 .0 0 0 .9 6 P A S S C [7 2 .9 8 ] 3 5 P o ly m o rp h is m c. 4 2 7 A > G p .M et 1 4 3 V al C o il ed -c o il d o m . 0 .9 3 1 .0 0 1 .0 0 0 .7 1 N o im p ac t 3 5 P o ly m o rp h is m c. 5 0 6 G > A p .A rg 1 6 9 G ln C o il ed -c o il d o m . 0 .8 6 0 .1 8 1 .0 0 0 .5 0 N o im p ac t 4 7 P o ly m o rp h is m c. 5 5 9 C > A p .G ln 1 8 7 L ys C o il ed -c o il d o m . 0 .8 6 1 .0 0 1 .0 0 0 .5 4 PA SS C [7 8 .2 6 ] 4 1 P o ly m o rp h is m c. 5 7 4 G > A p .A sp 1 9 2 A sn 0 .9 3 0 .4 3 1 .0 0 0 .7 5 N o im p ac t 3 5 P o ly m o rp h is m T h e U M D -p re d ic to r R© to o l co m p u te s al l el em en ts (s tr u ct u re ,b io ch em is tr y, sp li ci n g, an d co n se rv at io n ) an d p ro vi d es fo r ea ch o f th em a sp ec ifi c st re n gt h b as ed o n it s re la ti ve im p ac t. “C o n se rv at io n ”: T h e sc o re ra n ge s fr o m 0 fo r n o n co n se rv ed to 1 fo r fu ll y co n se rv ed . “S IF T ”: (s o rt in g in to le ra n t fr o m to le ra n t) :u se s se q u en ce h o m o lo gy to p re d ic t w h et h er an am in o ac id su b st it u ti o n w il la ff ec t p ro te in fu n ct io n an d ,h en ce ,p o te n ti al ly co n fe r a p h en o ty p e [N g an d H en ik o ff ,2 0 0 3 ]. Sc o re s th at ar e ≤ 0 .0 5 ar e co n si d er ed to b e d el et er io u s (i n b o ld ). “B L O SU M 6 2 ” is an am in o ac id su b st it u ti o n m at ri x b as ed o n lo ca l m u lt ip le al ig n m en ts o f an u n se le ct ed p ro te in se t o f re la te d se q u en ce s [H en ik o ff an d H en ik o ff ,1 9 9 2 ]. P o si ti ve sc o re s ar e as so ci at ed w it h co n se rv at iv e ch an ge s an d n eg at iv e sc o re s (i n b o ld ) w it h le ss co n se rv at iv e ch an ge s. “B io ch em ic al va lu e” is an am in o ac id su b st it u ti o n m at ri x d ep en d in g o n 4 8 q u al it at iv e p h ys ic o ch em ic al p ro p er ti es d es cr ib in g si d e ch ai n st ru ct u re an d fu n ct io n al gr o u p s, o p ti ca l p ro p er ti es ,h yd ro p h o b ic it y/ ch ar ge /a ci d – b as e p ro p er ti es ,a n d si ze (v o lu m e an d si d e ch ai n le n gt h ) [Y u ,2 0 0 1 ]. A va lu e b el o w 0 .0 5 d efi n es an in va li d su b st it u ti o n co n ce rn in g th e p h ys ic o -c h em ic al p ro p er ti es (i n b o ld ). “E SE ” (e xo n ic sp li ci n g en h an ce r) an d “S p li ce si te ” m o d ifi ca ti o n s ar e ev al u at ed b y th e U M D to o l. “P A SS C ”: p o te n ti al ac ce p to r sp li ce si te cr ea te d ;p o te n ti al d el et er io u s m o d ifi ca ti o n s ar e sh o w n in b o ld . “P at h o ge n ic it y” :t h e U M D -p re d ic to r R© to o l co m p u te s al l th e el em en ts (s tr u ct u re ,b io ch em is tr y, sp li ci n g, an d co n se rv at io n ) p ro vi d in g fo r ea ch a sp ec ifi c st re n gt h b as ed o n th ei r re la ti ve im p ac t. P re d ic ti o n s ar e n o rm al iz ed o n a sc al e fr o m 0 to 1 0 0 . In “C o n cl u si o n ”: a va lu e o f le ss th an 5 0 is as so ci at ed w it h th e p re d ic ti o n o f a n o n p at h o ge n ic m u ta ti o n an n o ta te d as “p o ly m o rp h is m ”; a va lu e o f 5 0 – 6 4 is as so ci at ed w it h th e p re d ic ti o n o f a “p ro b ab le p o ly m o rp h is m ”; a va lu e o f 6 5 – 7 4 is as so ci at ed w it h th e p re d ic ti o n o f a “p ro b ab ly p at h o ge n ic ” m u ta ti o n ;w h il e a va lu e ab o ve 7 4 is as so ci at ed w it h th e p re d ic ti o n o f a “p at h o ge n ic ” m u ta ti o n .F o r “S tr u ct u re ” an d “H C D ” d es cr ip ti o n s se e th e le ge n d s to T ab le 1 .F o r m o re d et ai ls se e F re d er ic et al .[ 2 0 0 9 ]. a T h es e m u ta ti o n s h av e b ee n st u d ie d in vi tr o [F u ch s et al ., 2 0 0 9 ;G av ar in i et al ., 2 0 1 0 ;K ai se r et al ., 2 0 1 0 ]. HUMAN MUTATION, Vol. 00, No. 0, 1–12, 2011 9 and have to be determined by searching, among others, for the presence of quasi-palindromic sequences, inverted repeats, or sym- metric elements that facilitate the formation of secondary-structure intermediates [Krawczak and Cooper, 1991]. Splice-Site Mutations Only one splice site variation is described in the literature: c.72-4T>C (or c.IVS1-4T>C) in association with the missense mu- tation c.77C>G [Houlden et al., 2010]. The splice site variation has been analyzed with the UMD algorithm implemented in HSF tool [Desmet et al., 2009]. The software computed the variation of the consensus value (CV) of the mutant “splice site” (CV = 87.25) rel- ative to the wild-type site (CV = 86.59). As the variation between these values was not higher than 10%, it is highly probable that such CV variation does not disrupt the wild-type splice site and thus this variation corresponds to a polymorphism as suggested also by the authors [Houlden et al., 2010]. The pathogenic mutation is thus likely to be c.77C>G and consequently the c.72-4T>C (or c.IVS1- 4T>C) polymorphism was associated to the c.77C>G mutation in the database. Conclusion The UMD-THAP1 database contains 53 different mutations (Ta- ble 1) in 56 probands and 43 relatives. These mutations are mainly private, essentially missense, usually nonrecurrent, and widely dis- tributed in the THAP domain and the rest of the THAP1 gene. The currently reported THAP1 mutations are mainly missense mu- tations (64.9%, 37/57) and small out-of-frame deletions (19.3%, 11/57), with a minority of other mutation types (7% of nonsense [4/57], 3.5% of small in-frame deletions [2/57], 1.8% of small out- of-frame insertions [1/57], and 3.5% of complex mutations). To date, no clear genotype/phenotype relationship has been identified. This is not surprising when considering past experience with other disease genes, such as the Fibrillin-1 gene (FBN1), for which many mutations were accumulated before the genotype/phenotype rela- tionships emerged. Indeed, the UMD-FBN1 database was created in 1995 [Collod et al., 1996], and the first report about the correlations between genotype and phenotype was published in 2007 [Faivre et al., 2007] and more followed [Detaint et al., 2010; Faivre et al., 2007, 2008, 2009a, b, c; Stheneur et al., 2009]. From this perspective, the UMD-THAP1 locus-specific database will facilitate the muta- tional analysis of THAP1 gene as previously demonstrated for other genes. Acknowledgments A.B. and V.E. are supported by PhD studentships from MESR (Ministère de l’Enseignement Supérieur et de la Recherche). The research leading to these results has also received funding from the Université Montpellier1, INSERM, the French Ministry of Health (National PHRC 2007-A00614-49), AMADYS-LFCD, Alliance France Dystonie, Lions Club La Grande Motte, the French Dystonia Network, and from the European Community’s Sev- enth Framework Programme (FP7/2007-2013) under grant agreement no. 200754—the GEN2PHEN project. References Almasy L, Bressman S, Raymond D, Kramer P, Greene P, Heiman G, Ford B, Yount J, de Leon D, Chouinard S, Saunders-Pullman R, Brin M, Kapoor R, Jones A, Shen H, Fahn S, Risch N, Nygaard T. 1997. Idiopathic torsion dystonia linked to chromosome 8 in two Mennonite families. Ann Neurol 42:670–673. Bernardi R, Pandolfi PP. 2007. Structure, dynamics and functions of promyelocytic leukaemia nuclear bodies. Nat Rev Mol Cell Biol 8:1006–1016. Beroud C, Collod-Beroud G, Boileau C, Soussi T, Junien C. 2000. UMD (Universal mu- tation database): a generic software to build and analyze locus-specific databases. Hum Mutat 15:86–94. Beroud C, Hamroun D, Collod-Beroud G, Boileau C, Soussi T, Claustres M. 2005. UMD (Universal mutation database): 2005 update. Hum Mutat 26:184–191. Bessiere D, Lacroix C, Campagne S, Ecochard V, Guillet V, Mourey L, Lopez F, Czaplicki J, Demange P, Milon A, Girard JP, Gervais V. 2008. Structure-function analysis of the THAP zinc finger of THAP1, a large C2CH DNA-binding module linked to Rb/E2F pathways. J Biol Chem 283:4352–4363. Blanchard A, Roubertie A, Frederic MY, Claustres M, Collod-Beroud G. 2010. Mono- genetic dystonia: revisiting the dopaminergic hypothesis. Rev Neurol (Paris) 166:389–399. Blanchard A, Roubertie A, Simonetta-Moreau M, Ea V, Coquart C, Frederic M, Gallouedec G, Adenis J, Benatru I, Borg M, Burbaud P, Calvas P, Cif L, Damier P, Destee A, Faivre L, Guyant-Marechal L, Janik P, Janoura S, Kreisler A, Lusakowska A, Odent S, Potulska-Chromik A, Rudzińska M, Thobois S, Vuillaume I, Tran- chant C, Tuffery-Giraud S, Coubes P, Sablonnière B, Claustres M, Collod-Béroud G. 2011. Singular DYT6 phenotypes in association with new THAP1 frameshift mutations. Mov Disord, DOI: 10.1002/mds.23641 in press. Bonetti M, Barzaghi C, Brancati F, Ferraris A, Bellacchio E, Giovanetti A, Ialongo T, Zorzi G, Piano C, Petracca M, Albanese A, Nardocci N, Dallapiccola B, Bentivoglio AR, Garavaglia B, Valente EM. 2009. Mutation screening of the DYT6/THAP1 gene in Italy. Mov Disord 24:2424–2427. Breakefield XO, Blood AJ, Li Y, Hallett M, Hanson PI, Standaert DG. 2008. The patho- physiological basis of dystonias. Nat Rev Neurosci 9:222–234. Bressman SB, Raymond D, Fuchs T, Heiman GA, Ozelius LJ, Saunders-Pullman R. 2009. Mutations in THAP1 (DYT6) in early-onset dystonia: a genetic screening study. Lancet Neurol 8:441–446. Campagne S, Saurel O, Gervais V, Milon A. 2010. Structural determinants of specific DNA-recognition by the THAP zinc finger. Nucleic Acids Res 38:3466–3476. Cariello N, Cui L, Béroud C, Soussi T. 1994a. Database and software for the analysis of mutations in the human p53 gene. Cancer Res 54:4454–4460. Cariello NF, Beroud C, Soussi T. 1994b. Database and software for the analysis of mutations at the human p53 gene. Nucleic Acids Res 22:3549–3550. Cayrol C, Lacroix C, Mathe C, Ecochard V, Ceribelli M, Loreau E, Lazar V, Dessen P, Mantovani R, Aguilar L, Girard JP. 2007. The THAP-zinc finger protein THAP1 regulates endothelial cell proliferation through modulation of pRB/E2F cell-cycle target genes. Blood 109:584–594. Cheng FB, Wan XH, Feng JC, Wang L, Yang YM, Cui LY. 2010. Clinical and genetic evaluation of DYT1 and DYT6 primary dystonia in China. Eur J Neurol 18:497– 503. Clot F, Grabli D, Burbaud P, Aya M, Derkinderen P, Defebvre L, Damier P, Krys- tkowiak P, Pollak P, Leguern E, San C, Camuzat A, Roze E, Vidailhet M, Durr A, Brice A. 2010. Screening of the THAP1 gene in patients with early-onset dys- tonia: myoclonic jerks are part of the dystonia 6 phenotype. Neurogenetics 12: 87–89. Clouaire T, Roussigne M, Ecochard V, Mathe C, Amalric F, Girard JP. 2005. The THAP domain of THAP1 is a large C2CH module with zinc-dependent sequence-specific DNA-binding activity. Proc Natl Acad Sci U S A 102:6907–6912. Collod G, Beroud C, Soussi T, Junien C, Boileau C. 1996. Software and database for the analysis of mutations in the human FBN1 gene. Nucleic Acids Res 24:137–140. De Carvalho Aguiar P, Fuchs T, Borges V, Lamar KM, Silva SM, Ferraz HB, Ozelius L. 2010. Screening of Brazilian families with primary dystonia reveals a novel THAP1 mutation and a de novo TOR1A GAG deletion. Mov Disord 25:2854–2857. Defazio G, Berardelli A, Hallett M. 2007. Do primary adult-onset focal dystonias share aetiological factors? Brain 130:1183–1193. Desmet FO, Hamroun D, Lalande M, Collod-Beroud G, Claustres M, Beroud C. 2009. Human splicing finder: an online bioinformatics tool to predict splicing signals. Nucleic Acids Res 37:e67. Detaint D, Faivre L, Collod-Beroud G, Child AH, Loeys BL, Binquet C, Gautier E, Arbustini E, Mayer K, Arslan-Kirchner M, Stheneur C, Halliday D, Beroud C, Bonithon-Kopp C, Claustres M, Plauchu H, Robinson PN, Kiotsekoglou A, De Backer J, Ades L, Francke U, De Paepe A, Boileau C, Jondeau G. 2010. Cardiovas- cular manifestations in men and women carrying a FBN1 mutation. Eur Heart J 31:2223–2229. Djarmati A, Schneider SA, Lohmann K, Winkler S, Pawlack H, Hagenah J, Bruggemann N, Zittel S, Fuchs T, Rakovic A, Schmidt A, Jabusch HC, Wilcox R, Kostic VS, Siebner H, Altenmuller E, Munchau A, Ozelius LJ, Klein C. 2009. Mutations in THAP1 (DYT6) and generalised dystonia with prominent spasmodic dysphonia: a genetic screening study. Lancet Neurol 8:447–452. Duan W, Rangnekar VM, Mattson MP. 1999a. Prostate apoptosis response-4 produc- tion in synaptic compartments following apoptotic and excitotoxic insults: evi- dence for a pivotal role in mitochondrial dysfunction and neuronal degeneration. J Neurochem 72:2312–2322. 10 HUMAN MUTATION, Vol. 00, No. 0, 1–12, 2011 Duan W, Zhang Z, Gash DM, Mattson MP. 1999b. Participation of prostate apoptosis response-4 in degeneration of dopaminergic neurons in models of Parkinson’s disease. Ann Neurol 46:587–597. Ducamp S, Kannengiesser C, Touati M, Garcon L, Guerci-Bresler A, Guichard JF, Vermylen C, Dochir J, Poirel HA, Fouyssac F, Mansuy L, Leroux G, Tertian G, Girot R, Heimpel H, Matthes T, Talbi N, Deybach JC, Beaumont C, Puy H, Grandchamp B. 2011. Sideroblastic anemia: molecular analysis of the ALAS2 gene in a series of 29 probands and functional studies of 10 missense mutations. Hum Mutat 32:590–597. Faivre L, Collod-Beroud G, Callewaert B, Child A, Binquet C, Gautier E, Loeys BL, Arbustini E, Mayer K, Arslan-Kirchner M, Stheneur C, Kiotsekoglou A, Comeglio P, Marziliano N, Wolf JE, Bouchot O, Khau-Van-Kien P, Beroud C, Claustres M, Bonithon-Kopp C, Robinson PN, Ades L, De Backer J, Coucke P, Francke U, De Paepe A, Jondeau G, Boileau C. 2009a. Clinical and mutation-type analysis from an international series of 198 probands with a pathogenic FBN1 exons 24-32 mutation. Eur J Hum Genet 17:491–501. Faivre L, Collod-Beroud G, Callewaert B, Child A, Loeys BL, Binquet C, Gautier E, Arbustini E, Mayer K, Arslan-Kirchner M, Kiotsekoglou A, Comeglio P, Grasso M, Beroud C, Bonithon-Kopp C, Claustres M, Stheneur C, Bouchot O, Wolf JE, Robinson PN, Ades L, De Backer J, Coucke P, Francke U, De Paepe A, Boileau C, Jondeau G. 2009b. Pathogenic FBN1 mutations in 146 adults not meeting clinical diagnostic criteria for Marfan syndrome: further delineation of type 1 fibrillinopathies and focus on patients with an isolated major criterion. Am J Med Genet A 149A:854–860. Faivre L, Collod-Beroud G, Child A, Callewaert B, Loeys BL, Binquet C, Gautier E, Arbustini E, Mayer K, Arslan-Kirchner M, Stheneur C, Kiotsekoglou A, Comeglio P, Marziliano N, Halliday D, Beroud C, Bonithon-Kopp C, Claustres M, Plauchu H, Robinson PN, Ades L, De Backer J, Coucke P, Francke U, De Paepe A, Boileau C, Jondeau G. 2008. Contribution of molecular analyses in diagnosing Marfan syndrome and type I fibrillinopathies: an international study of 1009 probands. J Med Genet 45:384–390. Faivre L, Collod-Beroud G, Loeys BL, Child A, Binquet C, Gautier E, Callewaert B, Arbustini E, Mayer K, Arslan-Kirchner M, Kiotsekoglou A, Comeglio P, Marziliano N, Dietz HC, Halliday D, Beroud C, Bonithon-Kopp C, Claustres M, Muti C, Plauchu H, Robinson PN, Ades LC, Biggin A, Benetts B, Brett M, Holman KJ, De Backer J, Coucke P, Francke U, De Paepe A, Jondeau G, Boileau C. 2007. Effect of mutation type and location on clinical outcome in 1,013 probands with Marfan syndrome or related phenotypes and FBN1 mutations: an international study. Am J Hum Genet 81:454–466. Faivre L, Masurel-Paulet A, Collod-Beroud G, Callewaert BL, Child AH, Stheneur C, Binquet C, Gautier E, Chevallier B, Huet F, Loeys BL, Arbustini E, Mayer K, Arslan-Kirchner M, Kiotsekoglou A, Comeglio P, Grasso M, Halliday DJ, Beroud C, Bonithon-Kopp C, Claustres M, Robinson PN, Ades L, De Backer J, Coucke P, Francke U, De Paepe A, Boileau C, Jondeau G. 2009c. Clinical and molecular study of 320 children with Marfan syndrome and related type I fibrillinopathies in a series of 1009 probands with pathogenic FBN1 mutations. Pediatrics 123:391–398. Frederic MY, Hamroun D, Faivre L, Boileau C, Jondeau G, Claustres M, Beroud C, Collod-Beroud G. 2008. A new locus-specific database (LSDB) for mutations in the TGFBR2 gene: UMD-TGFBR2. Hum Mutat 29:33–38. Frederic MY, Lalande M, Boileau C, Hamroun D, Claustres M, Beroud C, Collod- Beroud G. 2009. UMD-predictor, a new prediction tool for nucleotide substitution pathogenicity—application to four genes: FBN1, FBN2, TGFBR1, and TGFBR2. Hum Mutat 30:952–959. Fuchs T, Gavarini S, Saunders-Pullman R, Raymond D, Ehrlich ME, Bressman SB, Ozelius LJ. 2009. Mutations in the THAP1 gene are responsible for DYT6 primary torsion dystonia. Nat Genet 41:286–288. Gavarini S, Cayrol C, Fuchs T, Lyons N, Ehrlich ME, Girard JP, Ozelius LJ. 2010. Direct interaction between causative genes of DYT1 and DYT6 primary dystonia. Ann Neurol 68:549–553. Girard JP, Amalric F, Roussigne M, Clouaire T. 2010. Novel death associated proteins, and THAP1 and PAR4 pathway in apoptosis control. United States, Patent version number 20100021482. Groen JL, Ritz K, Contarino MF, van de Warrenburg BP, Aramideh M, Foncke EM, van Hilten JJ, Schuurman PR, Speelman JD, Koelman JH, de Bie RM, Baas F, Tijssen MA. 2010. DYT6 dystonia: mutation screening, phenotype, and response to deep brain stimulation. Mov Disord 25:2420–2427. Grundmann K, Reischmann B, Vanhoutte G, Hubener J, Teismann P, Hauser TK, Bonin M, Wilbertz J, Horn S, Nguyen HP, Kuhn M, Chanarat S, Wolburg H, Van der Linden A, Riess O. 2007. Overexpression of human wildtype torsinA and human DeltaGAG torsinA in a transgenic mouse model causes phenotypic abnormalities. Neurobiol Dis 27:190–206. Guo Q, Fu W, Xie J, Luo H, Sells SF, Geddes JW, Bondada V, Rangnekar VM, Mattson MP. 1998. Par-4 is a mediator of neuronal degeneration associated with the pathogenesis of Alzheimer disease. Nat Med 4:957–962. Hamroun D, Kato S, Ishioka C, Claustres M, Beroud C, Soussi T. 2006. The UMD TP53 database and website: update and revisions. Hum Mutat 27:14–20. Henikoff S, Henikoff JG. 1992. Amino acid substitution matrices from protein blocks. Proc Natl Acad Sci USA 89:10915–10919. Houlden H, Schneider SA, Paudel R, Melchers A, Schwingenschuh P, Edwards M, Hardy J, Bhatia KP. 2010. THAP1 mutations (DYT6) are an additional cause of early-onset dystonia. Neurology 74:846–850. Jech R, Bares M, Krepelova A, Urgosik D, Havrankova P, Ruzicka E. 2011. DYT 6—a novel THAP1 mutation with excellent effect on pallidal DBS. Mov Disord 26:924– 925. Kaiser FJ, Osmanoric A, Rakovic A, Erogullari A, Uflacker N, Braunholz D, Lohnau T, Orolicki S, Albrecht M, Gillessen-Kaesbach G, Klein C, Lohmann K. 2010. The dystonia gene DYT1 is repressed by the transcription factor THAP1 (DYT6). Ann Neurol 68:554–559. Kannengiesser C, Brookes S, del Arroyo AG, Pham D, Bombled J, Barrois M, Mauffret O, Avril MF, Chompret A, Lenoir GM, Sarasin A, Peters G, Bressac-de Paillerets B. 2009. Functional, structural, and genetic evaluation of 20 CDKN2A germ line mutations identified in melanoma-prone families or patients. Hum Mutat 30:564– 574. Krahn M, Beroud C, Labelle V, Nguyen K, Bernard R, Bassez G, Figarella-Branger D, Fernandez C, Bouvenot J, Richard I, Ollagnon-Roman E, Bevilacqua JA, Salvo E, Attarian S, Chapon F, Pellissier JF, Pouget J, Hammouda el H, Laforet P, Urtizberea JA, Eymard B, Leturcq F, Levy N. 2009. Analysis of the DYSF mutational spectrum in a large cohort of patients. Hum Mutat 30:E345–E375. Krawczak M, Cooper D. 1991. Gene deletions causing human genetic disease: mecha- nisms of mutagenesis and the role of the local DNA sequence environment. Hum Genet 86:425–441. Lacroix C. 2007. THAP1, un régulateur clé de la prolifération des cellules endothéliales: relations structure/fonction et gènes cibles. Toulouse: Université Toulouse III-Paul Sabatier . 154 p. Lallemand-Breitenbach V, de The H. 2010. PML nuclear bodies. Cold Spring Harb Perspect Biol 2:a000661. Mazars R, Gonzalez-de-Peredo A, Cayrol C, Lavigne AC, Vogel JL, Ortega N, Lacroix C, Gautier V, Huet G, Ray A, Monsarrat B, Kristie TM, Girard JP. 2010. The THAP- zinc finger protein THAP1 associates with coactivator HCF-1 and O-GlcNAc transferase: a link between DYT6 and DYT3 dystonias. J Biol Chem 285:13364– 13371. Ng PC, Henikoff S. 2003. SIFT: predicting amino acid changes that affect protein function. Nucleic Acids Res 31:3812–3814. Norgren N, Mattson E, Forsgren L, Holmberg M. 2011. A high-penetrance form of late- onset torsion dystonia maps to a novel locus (DYT21) on chromosome 2q14.3- q21.3. Neurogenetics 12:137–143. Paisan-Ruiz C, Ruiz-Martinez J, Ruibal M, Mok KY, Indakoetxea B, Gorostidi A, Masso JF. 2009. Identification of a novel THAP1 mutation at R29 amino-acid residue in sporadic patients with early-onset dystonia. Mov Disord 24:2428– 2429. Park SK, Nguyen MD, Fischer A, Luke MP, Affar el B, Dieffenbach PB, Tseng HC, Shi Y, Tsai LH. 2005. Par-4 links dopamine signaling and depression. Cell 122:275– 287. Pedersen WA, Luo H, Kruman I, Kasarskis E, Mattson MP. 2000. The prostate apoptosis response-4 protein participates in motor neuron degeneration in amyotrophic lateral sclerosis. FASEB J 14:913–924. Roussigne M, Cayrol C, Clouaire T, Amalric F, Girard JP. 2003. THAP1 is a nuclear proapoptotic factor that links prostate-apoptosis-response-4 (Par-4) to PML nu- clear bodies. Oncogene 22:2432–2442. Saunders-Pullman R, Raymond D, Senthil G, Kramer P, Ohmann E, Deligtisch A, Shanker V, Greene P, Tabamo R, Huang N, Tagliati M, Kavanagh P, Soto-Valencia J, Aguiar Pde C, Risch N, Ozelius L, Bressman S. 2007. Narrowing the DYT6 dystonia region and evidence for locus heterogeneity in the Amish-Mennonites. Am J Med Genet A 143A:2098–2105. Schneider SA, Ramirez A, Shafiee K, Kaiser FJ, Erogullari A, Bruggemann N, Winkler S, Bahman I, Osmanovic A, Shafa MA, Bhatia KP, Najmabadi H, Klein C, Lohmann K. 2011. Homozygous THAP1 mutations as cause of early-onset generalized dystonia. Mov Disord 26:858–861. Söhn AS, Glockle N, Doetzer AD, Deuschl G, Felbor U, Topka HR, Schols L, Riess O, Bauer P, Muller U, Grundmann K. 2010. Prevalence of THAP1 sequence variants in German patients with primary dystonia. Mov Disord 25:1982–1986. Stheneur C, Collod-Beroud G, Faivre L, Buyck JF, Gouya L, Le Parc JM, Moura B, Muti C, Grandchamp B, Sultan G, Claustres M, Aegerter P, Chevallier B, Jondeau G, Boileau C. 2009. Identification of the minimal combination of clinical features in probands for efficient mutation detection in the FBN1 gene. Eur J Hum Genet 17:1121–1128. Valente EM, Albanese A. 2010. Advances in the genetics of primary torsion dystonia. F1000 Biol Rep 2:41. Van Gerpen JA, Ledoux MS, Wszolek ZK. 2010. Adult-onset leg dystonia due to a missense mutation in THAP1. Mov Disord 25:1306–1307. Vidailhet M, Grabli D, Roze E. 2009. Pathophysiology of dystonia. Curr Opin Neurol 22:406–413. HUMAN MUTATION, Vol. 00, No. 0, 1–12, 2011 11 Villeger L, Abifadel M, Allard D, Rabes JP, Thiart R, Kotze MJ, Beroud C, Junien C, Boileau C, Varret M. 2002. The UMD-LDLR database: additions to the software and 490 new entries to the database. Hum Mutat 20:81–87. Wautot V, Vercherat C, Lespinasse J, Chambe B, Lenoir GM, Zhang CX, Porchet N, Cordier M, Beroud C, Calender A. 2002. Germline mutation profile of MEN1 in multiple endocrine neoplasia type 1: search for correlation between phenotype and the functional domains of the MEN1 protein. Hum Mutat 20:35–47. Xiao J, Zhao Y, Bastian RW, Perlmutter JS, Racette BA, Tabbal SD, Karimi M, Paniello RC, Wszolek ZK, Uitti RJ, Van Gerpen JA, Simon DK, Tarsy D, Hedera P, Truong DD, Frei KP, Dev Batish S, Blitzer A, Pfeiffer RF, Gong S, LeDoux MS. 2010. Novel THAP1 sequence variants in primary dystonia. Neurology 74:229–238. Yu K. 2001. Theoretical determination of amino acid substitution groups based on qualitative physicochemical properties. http://cmgm.stanford.edu/biochem218/ Projects%202001/Yu.pdf. Zittel S, Moll CK, Bruggemann N, Tadic V, Hamel W, Kasten M, Lohmann K, Lohnau T, Winkler S, Gerloff C, Schonweiler R, Hagenah J, Klein C, Munchau A, Schneider SA. 2010. Clinical neuroimaging and electrophysiological assessment of three DYT6 dystonia families. Mov Disord 25:2405–2412. 12 HUMAN MUTATION, Vol. 00, No. 0, 1–12, 2011