Recessive mutations in POLR1C cause a leukodystrophy by impairing biogenesis of RNA polymerase III

Abstract

A small proportion of 4H (Hypomyelination, Hypodontia and Hypogonadotropic Hypogonadism) or RNA polymerase III (POLR3)-related leukodystrophy cases are negative for mutations in the previously identified causative genes POLR3A and POLR3B. Here we report eight of these cases carrying recessive mutations in POLR1C, a gene encoding a shared POLR1 and POLR3 subunit, also mutated in some Treacher Collins syndrome (TCS) cases. Using shotgun proteomics and ChIP sequencing, we demonstrate that leukodystrophy-causative mutations, but not TCS mutations, in POLR1C impair assembly and nuclear import of POLR3, but not POLR1, leading to decreased binding to POLR3 target genes. This study is the first to show that distinct mutations in a gene coding for a shared subunit of two RNA polymerases lead to selective modification of the enzymes' availability leading to two different clinical conditions and to shed some light on the pathophysiological mechanism of one of the most common hypomyelinating leukodystrophies, POLR3-related leukodystrophy.

ARTICLE

Received 17 Nov 2014 | Accepted 26 May 2015 | Published 7 Jul 2015

DOI: 10.1038/ncomms8623 OPEN

Recessive mutations in POLR1C cause a leukodystrophy by impairing biogenesis of RNA polymerase III

Isabelle Thiffault1,2,3,*, Nicole I. Wolf4,*, Diane Forget5,*, Kether Guerrero1, Luan T. Tran1, Karine Choquet6,

Mathieu Lavalle-Adam7, Christian Poitras5, Bernard Brais6, Grace Yoon8, Laszlo Sztriha9, Richard I. Webster10,11, Dagmar Timmann12, Bart P. van de Warrenburg13, Jrgen Seeger14, Alz Zimmermann9, Adrienn Mt15, Cyril Goizet16, Eva Fung17, Marjo S. van der Knaap4, Sbastien Fribourg18,19, Adeline Vanderver20,21,22, Cas Simons23, Ryan J. Taft22,23,24,25, John R. Yates III7, Benoit Coulombe5,26 & Genevive Bernard1

A small proportion of 4H (Hypomyelination, Hypodontia and Hypogonadotropic Hypogonadism) or RNA polymerase III (POLR3)-related leukodystrophy cases are negative for mutations in the previously identied causative genes POLR3A and POLR3B. Here we report eight of these cases carrying recessive mutations in POLR1C, a gene encoding a shared POLR1 and POLR3 subunit, also mutated in some Treacher Collins syndrome (TCS) cases. Using shotgun proteomics and ChIP sequencing, we demonstrate that leukodystrophy-causative mutations, but not TCS mutations, in POLR1C impair assembly and nuclear import of POLR3, but not POLR1, leading to decreased binding to POLR3 target genes. This study is the rst to show that distinct mutations in a gene coding for a shared subunit of two RNA polymerases lead to selective modication of the enzymes availability leading to two different clinical conditions and to shed some light on the pathophysiological mechanism of one of the most common hypomyelinating leukodystrophies, POLR3-related leukodystrophy.

1 Department of Neurology and Neurosurgery, McGill University, Department of Medical Genetics, Montreal Childrens Hospital, Research Institute of the McGill University Health Center, 1001 boul Dcarie, Montreal, Quebec H4A 3J1, Canada. 2 Service de Gntique, Centre Hospitalier Universitaire Sainte-Justine, 3175 Chemin de la Cte-Sainte-Catherine, Montreal, Quebec H3T1C5, Canada. 3 Center for Pediatric Genomic Medicine, Childrens Mercy Hospital, 2420 Pershing Road, Suite 421, Kansas City, Missouri 64108, USA.

4 Department of Child Neurology, VU University Medical Center, Neuroscience Campus Amsterdam, Amsterdam 1081 HZ, The Netherlands. 5 Translational Proteomics Laboratory, Institut de recherches cliniques de Montral (IRCM), 110 avenue des Pins ouest, Montral, Qubec H2W 1R7, Canada. 6 Neurogenetics of Motion Laboratory, Montreal Neurological Institute, 3801 University Street, McGill University, Montreal, Quebec H3A 2B4, Canada. 7 Department of Chemical Physiology, The Scripps Research Institute, 10550 North Torrey Pines Road SR302, La Jolla, California 92037, USA. 8 Division of Neurology and Clinical and Metabolic Genetics, the Hospital for Sick Children, University of Toronto, 555 University Avenue, Toronto, Ontario M5G 1X8, Canada. 9 Department of Paediatrics, Faculty of Medicine, University of Szeged, Temesvri krt. 35-37, Szeged H-6726, Hungary.

10 T.Y. Nelson Department of Neurology and Neurosurgery, The Childrens Hospital at Westmead, Locked Bag 4001, Westmead, New South Wales 2145, Australia. 11 Institute for Neuroscience and Muscle Research, The Childrens Hospital at Westmead, Locked Bag 4001, Westmead New South Wales 2145, Australia. 12 Department of Neurology, University Clinic Essen, University of Duisburg-Essen, Hufelandstrasse 55, 45147 Essen, Germany. 13 Department of Neurology, Donders Institute for Brain, Cognition, and Behaviour, Radboud University Medical Center, PO Box 9101, Nijmegen 6500 HB, The Netherlands. 14 Department of Pediatrics and Adolescent Medicine, Deutsche KlinikfrDiagnostik, Wiesbaden 65191, Germany. 15 Department of Neurosurgery, Faculty of Medicine, University of Szeged, 6 Semmelweis Street, Szeged H-6725, Hungary. 16 Service de Gntique, Hpital Pellegrin, CHU Bordeaux and University Bordeaux, Laboratoire MRGM (EA4576), Bordeaux 33076, France. 17 Department of Paediatrics, The Chinese University of Hong Kong, Prince of Wales Hospital, Shatin, Hong Kong, SAR China. 18 Universit de Bordeaux, Institut Europen de Chimie et Biologie, ARNA Laboratory, Pessac F-33607, France. 19 Institut National de la Sant Et de la Recherche Mdicale, INSERMU869, ARNA Laboratory, Bordeaux F-33000, France. 20 Center for Genetic Medicine Research, Childrens National, 111 Michigan Avenue Northwest, Washington, District of Columbia 20010, USA. 21 Department of Neurology, Childrens National, 111 Michigan Avenue Northwest, Washington, District of Columbia 20010, USA. 22 George Washington University, School of Medicine, Washington, District of Columbia 20052, USA. 23 Institute for Molecular Bioscience, University of Queensland, Brisbane, Queensland 4072, Australia. 24 Departments of Integrative Systems Biology and Pediatrics, School of Medicine and Health Sciences, The George Washington University, Washington, District of Columbia 20037, USA. 25 Illumina Inc., 5200 Illumina Way, San Diego, California 92122, USA. 26 Department of Biochemistry, Universit de Montral, Pavillon Roger-Gaudry, CP 6128, Succ Centre-Ville, Montreal, Qubec H3C 3J7, Canada. * These authors contributed equally to this work. Correspondence and requests for materials should be addressed to B.C. (email: mailto:[email protected]

Web End [email protected] ) or to G.B. (email: mailto:[email protected]

Web End [email protected] ).

NATURE COMMUNICATIONS | 6:7623 | DOI: 10.1038/ncomms8623 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 1

& 2015 Macmillan Publishers Limited. All rights reserved.

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms8623

Leukodystrophies are a heterogeneous group of genetically determined disorders characterized by abnormal white matter on brain imaging1,2. They are classied as

hypomyelinating and non-hypomyelinating leukodystrophies based on magnetic resonance imaging (MRI) characteristics2, depending on whether the principal problem appears to be a lack of myelin deposition during development or altered myelin homeostasis. RNA polymerase III (POLR3)-related leukodystrophy or 4H (Hypomyelination, Hypodontia and Hypogonadotropic Hypogonadism) leukodystrophy (MIM 607694, 614381)3 was found to be caused by recessive mutations in POLR3A (MIM 614258) or POLR3B (MIM 614366)49 and is characterized by an expanding spectrum of clinical3,10 and radiological features1012. POLR3A and POLR3B are, respectively, the largest and second largest of the 17 subunits that constitute POLR3. Together, they form the catalytic centre of the enzyme. POLR3 synthetizes small non-coding RNAs, including tRNAs, 5S RNA, 7SK RNA and U6 RNA, that are involved in the regulation of essential cellular processes, including transcription, RNA processing and translation13. A subset of patients (B5%) presenting with compatible clinical and/or radiological features of POLR3-related leukodystrophy have no detectible mutations in either POLR3A or POLR3B, suggesting that mutations in one or more additional genes may result in this presentation. We hypothesized that genes coding for other POLR3 subunits or for proteins interacting with POLR3 would be strong candidates in these cases.

In this study, we identied recessive mutations in POLR1C, a gene encoding for a subunit common to POLR1 and POLR3, which has thus far been known only to be associated with autosomal recessive Treacher Collins syndrome (TCS)14. We also demonstrated that the leukodystrophy-causing mutations affect POLR3, but not POLR1 assembly and nuclear import, leading to decrease binding to POLR3 target genes, whereas one TCS mutation leads to normal assembly of both polymerases, but rather affects POLR1 targeting to the nucleolus, the site of POLR1 gene transcription.

ResultsWhole-exome and Sanger sequencing for gene identication. To investigate the genetic aetiology of these unexplained cases, we performed exome sequencing in three cases with typical clinical and/or radiological features of POLR3-related leukodystrophy negative for POLR3A or POLR3B mutation. Genome coverage for these three cases exceeded 56 (Supplementary Table 1)

and yielded more than 4.6 million variants per genome. Review of the data using in-house lters revealed a homozygous mutation in POLR1C (NM_203290; GRCh37/hg19) in two patients (c.221A4G (p.Asn74Ser) and c.95A4T (p.Asn32Ile), respectively), and compound heterozygous mutations in one(c.436T4C (p.Cys146Arg) and c.883_885delAAG (p.Lys295-del)), all consistent with autosomal recessive inheritance (Fig. 1, Supplementary Table 1). These four variants were validated using Sanger sequencing. Co-segregation analysis conrmed that the mutations were inherited from heterozygous carrier parents. These variants were not observed in over 6,500 individuals of European and African American ancestry in the NHLBI Exome Sequencing Project database, in the ExAc data set nor in our internal variant databases. Furthermore, these variants were not present in more than 300 Centre dEtude du Polymorphisme Humain (CEPH) control chromosomes. In silico analysis predicted those to be disease-causing (Supplementary Table 2c). These variants were found to affect highly conserved amino-acid residues (Fig. 1). To identify additional mutations, we sequenced all exons, exonintron boundaries and 30 and 50 untranslated repeat of POLR1C (see Supplementary Table 3 for primers) in 16 other individuals selected on the basis of clinical (compatible neurological features with at least one non-neurological feature such as dental abnormalities or hypogonadotropic hypogonadism, together with hypomyelination on MRI) and/or radiological characteristics (typical MRI features of 4H leukodystrophy)3,10,11 but negative for POLR3A and POLR3B mutations and uncovered ve additional compound heterozygous or homozygous cases (Supplementary Table 2, Fig. 1). Sequencing of the mutated exons was performed in family members for whom DNA was available

a

POLR1C - Polymerase (RNA) I polypeptice C, 30 kDA GRCh37 (Chr 6)

Treacher collins syndrome

c.87deIT [p.Gly31ValfsX23]

c.835C>T [p.Arg279Trp]

c.836G>A [p.Arg279Gln]

c.923+3_c.922+6del

c.979A>T [p.Lys327X]

1 2 3 4 5 6 7 8 9

c.77C>T [p.Thr26IIe]c.95A>T [p.Asn32IIe]

c.193A>G [p.Met65Val] c.221A>G [p.Asn74Ser]

c.281T>C [p.Val94Ala] c.326G>A [p.Arg109Hls]

c.395G>A [p.Gly132Asp]

c.436T>C [p.Cys146Arg]

c.461_462delAA [p.Lys154Argfs*4]

c.572G>A [p.Arg191Gln]

c.785T>C [p.lle262Thr]

c. c.883_885delAAG [p.Lys295del]

c.970G>A [p.Glu324Lys]

POLR3-related leukodystrophy

Mutations

b

Figure 1 | POLR1C mutations in leukodystrophy and TCS cases. (a) Genomic organization of POLR1C in humans (UCSC Genome Browser hg19): mutations and their positions within the POLR1C gDNA; in light blue are mutations that cause TCS, mutations in black cause POLR3-related leukodystrophy.(b) POLR1C mutations in patients with leukodystrophy affect amino acids that are conserved through species.

2 NATURE COMMUNICATIONS | 6:7623 | DOI: 10.1038/ncomms8623 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

& 2015 Macmillan Publishers Limited. All rights reserved.

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms8623 ARTICLE

to conrm segregation. In total, 13 POLR1C mutations were detected in eight cases (Supplementary Table 2a, Fig. 1). Similar to what is observed in cases with the disease caused by mutations in POLR3A or POLR3B, clinical and radiological characteristics of these eight cases were compatible with POLR3-related leukodystrophy; however, patients did not necessarily have all clinical (that is, neurological, dental, ophthalmological and endocrine abnormalities) and MRI features1,3,12 (Table 1, Fig. 2, Supplementary Tables 4 and 5) of the disease.

Impact of POLR1C recessive mutations on POLR1 and POLR3. To investigate the potential pathogenic role of these mutations, we evaluated the impact of the two homozygous mutations (Table 1) on the function of nuclear POLR1 and POLR3. FLAG-tagged versions of the wild-type (WT) form of POLR1C and its variants having the p.Asn32Ile (N32I) or p.Asn74Ser (N74S) substitution were expressed in HeLa cells. Anti-FLAG afnity purication was performed on cell extracts and the puried proteins were analysed using shotgun proteomics. These experiments were performed in triplicate. The expression level of the various forms of POLR1C (that is, WT and mutants) were equivalent and comparable (see Supplementary Fig. 2), and protein expression levels were normalized by the expression level of the bait in each purication (Supplementary Table 6). The tagged WT POLR1C pulled down all subunits of both POLR1 and POLR3 (Fig. 3a), a nding that was expected since POLR1C is a shared subunit of both polymerases (see Fig. 3c for a schematic representation). Both tagged mutated POLR1C (N32I and N74S) pulled down amounts of POLR1-specic (POLR1A, 1B, 1E, CD3EAP, TWISTNB and ZNRD1) and POLR1/POLR3-shared subunits (POLR2E, 2F, 2H, 2K, 2L, 1C and 1D) that were not signicantly different from those pulled down by the WT (this is especially true for the POLR1-specic subunits). However, both mutated subunits pulled down lower amounts of POLR3 (most specic subunits) relative to WT. This nding suggests that the mutations lead to a selective defect in POLR3 assembly, and not in POLR1. Indeed, assembly of nuclear RNA polymerases has previously been shown to occur in the cells cytoplasm and defects in RNA polymerase assembly caused by functional disruption of the RNA Polymerase-Associated proteins was previously shown to lead to cytoplasmic accumulation of polymerase subunits1521. Notably, the position of mutated residues in the proposed structure of POLR3 is compatible with defects in enzyme

assembly and/or folding (see Supplementary Fig. 1). To further conrm that mutated POLR1C variants (N32I and N74S) are impaired in supporting enzyme assembly and nuclear import, immunouorescence studies were performed using anti-FLAG antibodies. The results reveal an accumulation of both mutated POLR1C subunits, but not the WT subunit, in the cytoplasm (Fig. 3d). We then performed chromatin immunoprecipitation (ChIP) of FLAG-tagged POLR1C followed by high-throughput sequencing (ChIP-Seq), as a proxy of gene transcription activity2224, to investigate the impact of POLR1C mutations on gene occupancy by POLR1 and POLR3. After alignment of the reads to the human reference genome (hg19), we compared occupancy of WT and mutated POLR1C variants over 659 POLR3-transcribed genes, including all transfer RNA (tRNA) and 5S ribosomal RNA genes (Supplementary Table 7). As expected, mutated POLR1C variants displayed reduced binding to POLR3-transcribed genes compared with WT POLR1C for all three classes of POLR3-transcribed genes (classied according to their regulatory elements; see legend to Fig. 4a). In contrast, there were no differences in WT and mutated POLR1C occupancy over the ribosomal RNA gene transcribed by POLR1 (Fig. 4c). Together, these results indicate that the N32I and N74S substitutions in the POLR1/POLR3-shared subunit POLR1C specically interfere with assembly, nuclear import and chromatin association of POLR3. To compare the roles of leukodystrophy versus TCS-causing mutations in the biogenesis of POLR1 and POLR3, we expressed FLAG-tagged POLR1C with the p.Arg279Gln (R279Q) mutation in HeLa cells, afnity-puried the tagged subunit and identied the puried interactors using mass spectrometry.Contrary to POLR1C (N32I) and POLR1C (N74S), none of the subunits of POLR1 and POLR3 were pulled down by tagged POLR1C (R279Q) in amounts that were statistically signicantly different from the WT (Fig. 5a,b, Supplementary Table 6), suggesting that this TCS-causing mutation does not affect the assembly of these polymerases. Notably, however, immunouorescence results indicate that POLR1C (R279Q) targeting to the nucleolus is impaired as compared with the WT subunit (Fig. 5c) and the N32I- and N74S-mutated subunits (see Fig. 3d).

DiscussionWith the advent of exome sequencing, it is becoming increasingly apparent that allelic heterogeneity in genes encoding essential

Table 1 | Major clinical and MRI ndings in index patients with mutations in POLR1C.

Patient numbers

Family nos

Gene Mutation 1 Mutation 2 Main clinical characteristics Main MRI characteristics

Neurological abnormalities

Myopia Dental abnormality

Hypogonadotropic hypogonadism

Hypomyelination

Thin corpus callosum

Cerebellar atrophy

1 I POLR1C c.95A4T;p.Asn32Ile

c.95A4T;p.Asn32Ile

Too young

2 II POLR1C c.221A4G;

p.Asn74Ser

c.221A4G;p.Asn74Ser

Too young

3 III POLR1C c.436T4C;

p.Cys146Arg

c.883_885delAAG;p.Lys295del

Too young

4 IV POLR1C c.77C4T;

p.Thr26Ile

c.326G4A;p.Arg109His

Too young

5 V POLR1C c.193A4G;

p.Met65Val

c.572G4A;p.Arg191Gln

Too young

c.970G4A;p.Glu324Lys 7 VII POLR1C c.395G4A;

p.Gly132Asp

c.461_462delAA;p.Lys154Argfs*4 8 VIII POLR1C c.281T4C;

p.Val94Ala

c.785T4C;p.Ile262Thr

6 VI POLR1C c.326G4A;

p.Arg109His

Too young

, absent; , present; MRI, magnetic resonance imaging.

NATURE COMMUNICATIONS | 6:7623 | DOI: 10.1038/ncomms8623 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 3

& 2015 Macmillan Publishers Limited. All rights reserved.

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms8623

a b c

d e f

g h i

Figure 2 | MRI characteristics of POLR3-related leukodystrophy caused by POLR1C mutations. Axial T2-weighted (a,b,d,e,g,h) and sagittal T1-weighted (c,f,i) images of case 1 aged 6 years (ac) and case 2 aged 4.5 years (df) compared with a healthy control aged 4 years (gi). Diffuse hyperintense signal of the supratentorial (red arrow, a,d) and cerebellar (blue arrow, b,e) white matter is visible on the T2-weighted images, indicating hypomyelination. There is no cerebellar atrophy. As typical for POLR3-related leukodystrophy, the ventrolateral thalamus (white arrow, a,d), the optic radiation (thick arrowhead blue, d) and the dentate nucleus (open red arrowhead, b) show a relative hypointense signal on the T2-weighted images resulting in an easily visible dentate nucleus (b) as compared with the control (h) as well as a small dot in the posterior limb of the internal capsule (red arrowhead, d). The corpus callosum is slightly thinned in case 1 and thinned in case 2 (open red arrowhead, c,f).

proteins, such as those involved in transcription, results in highly variable phenotypes. A previous report of mutations in POLR1C (ref. 14) highlighted the discovery of the rst cases of TCS with an autosomal recessive mode of inheritance. TCS (MIM 154500, 613717, 248390) is characterized by an abnormal craniofacial development and is caused, from most to least frequent, by mutations in TCOF1 (dominant), POLR1D (dominant) or POLR1C (recessive)14. TCS caused by mutations in POLR1D or POLR1C has been proposed to arise as a consequence of a decreased quantity of functional ribosomes in the neuroepithelium and the neuronal crest cells during critical points of embryogenesis14,25. We assessed the role of a TCS-causing mutation (R279Q) in the biogenesis of POLR1 and POLR3. Our results (see Fig. 5, Supplementary Table 6) indicate that this mutation does not impair polymerase assembly, as opposed to leukodystrophy-causing mutations (see Fig. 3), but affect targeting to the nucleolus, the site for Pol I transcription.

Our ndings suggest that improper assembly and nuclear import of POLR3 resulting from leukodystrophy-causative mutations lead to decreased availability of the complex at the chromatin. As POLR3 binding is well correlated to tRNA

expression2224, this decreased POLR3 occupancy is likely to cause reduced transcription of tRNAs and other essential small non-coding RNAs. One hypothesis is that mutations in POLR3A, POLR3B or POLR1C lead to decreased levels of certain tRNAs crucial for the synthesis of proteins essential for central nervous system myelin development. tRNA function has also been suggested to be impaired in other white matter disorders2630 caused by mutations in tRNA-aminoacyl synthetases, including hypomyelinating leukodystrophies such as RARS-associated hypomyelination30 and Hypomyelination with Brain Stem and spinal cord involvement and Leg spasticity29. Of note, tRNA synthetases have not been found in our POLR3 purications, suggesting that tRNA aminoacylation is not coupled with POLR3 transcription. An alternative hypothesis would involve changes in the expression of other essential small non-coding RNAs synthesized by POLR3.

In conclusion, our sequencing study of 18 cases with compatible clinical and/or radiological features with 4H or POLR3-related leukodystrophy identied 13 different POLR1C mutations in eight cases. POLR1C joins an emerging group of genes with dual roles in pathogenesis of human diseases31,

4 NATURE COMMUNICATIONS | 6:7623 | DOI: 10.1038/ncomms8623 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

& 2015 Macmillan Publishers Limited. All rights reserved.

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms8623 ARTICLE

a b

WT N321 N74S

Ratio<0.67

Ratio>1.5

POLR1A

POLR1B POLR1E

CD3EAP

ZNRD1

POLR3A

POLR3B POLR3C

POLR3D

POLR3E

POLR3F

POLR3G

POLR3H

CRCP *POLR1C

POLR1D

POLR2E

POLR2F

POLR2H POLR2K

POLR2L

2.5

2

1.5

1

0.5

0 3.5 3 2.5 2 1.5 1 0.5 0 0.5 1 1.5 2 2.5 3 3.5

Specific Pol l subunits

Specific Pol lll subunits

Shared subunits (Pol l and

Pol lll)

POLR3E

POLR3B

POLR1D

Log 10(t-test adjusted Pvalue)

POLR3G

POLR3F

POLR3H POLR3F POLR3E

POLR3A

POLR3A

Log2(average ratio N32I / WT)

Ratio<0.67 Ratio>1.5

Log 10(t-test adjusted p-value)

2.5

2

1.5

1

0.5

0 4.5 4 3.5 3 2.5 2 1.5 1 0.5 0 0.5 1 1.5 2 2.5 3

CRCP

MAF1

POLR2B

GPN1 POLR3D

POLR3G

TCOF1

POLR3C

RUVBL1

4 3.5 3 2.5 2 1.5 1 0.5 0 0.5 1

Log2(average ratio N74S / WT)

Pvalue <0.05Pvalue < 0.05

c d

FLAG DNA Merge

CD3EAP

ZNRD1

1E

2L

1B

2K

1D

2H

POLR1C (WT)

POLR1C (N32I)

POLR1C (N74S)

2E

1A

2F

TWISTNB Pol l

3D 3E

3B

2L

1C

2K

2E

3K 3C 3F 3G

3A

3H

1C

1D

2H

2F

CRCP Pol lll

20 m

Figure 3 | Impact of POLR1C mutations on polymerase assembly and nuclear import. (a) FLAG-tagged POLR1C variants, either the wild-type (1C) polypeptide or mutated versions having a N32I or a N74S substitution, were expressed in HeLa cells and puried using anti-FLAG afnity chromatography. The co-puried proteins were identied using LC-MS/MS mass spectrometry. The heatmap contains the log2-transformed average spectral count ratios

N32I or N74S/WTacross all three replicates. Spectral counts were computed with Mascot (see Supplementary Table 6 for the complete data set). Specic and shared POLR1 (Pol I) and POLR3 (Pol III) subunits are identied on the left. POLR1C (the bait) is identied by an asterisk. (b) Volcano plots of the log2-

transformed average spectral count ratios N32I or N74S/WT (x axis) and the log10-transformed P values (adjusted with the BenjaminiHochberg procedure) resulting from the two-tailed one-sample t-tests of the high-condence interactors of POLR1C. Red proteins show a level of differential interaction with POLR1C that is statistically signicant, while blue proteins do not. (c) Schematic representation of the subunit composition of POLR1 (Pol I) and POLR3 (Pol III; see refs 53,54 for details). Shared subunits are in grey and POLR1C in black. (d) Immunouorescence experiments showing the cellular localization of tagged POLR1C variants. Nuclei are stained using TO-PRO-3 iodine. Scale bar, 20 mm.

expanding the clinical phenotype associated with constitutional mutations in this gene and opening new aspects in the annotation and assessment of pathogenicity of sequence variants.

Furthermore, our functional analyses on the mutational impact of POLR1C bring the rst insights into the pathophysiology of POLR3-related disorders.

NATURE COMMUNICATIONS | 6:7623 | DOI: 10.1038/ncomms8623 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 5

& 2015 Macmillan Publishers Limited. All rights reserved.

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms8623

a

200 Type 1

100 Type 2

200 Type 3

Aggregate

signal

POLR1C (WT)

POLR1C (N32I)

POLR1C (N74S)

01,000 0 1,000

01,000 0 1,000

01,000 0 1,000

Distance to TSS (bp)

b c

26,286,000 26,286,500 26,287,000 26,287,500 10 kb 20 kb 30 kb 40 kb

5.8s

100 bp 2 kb

POLR1C (WT)

POLR1C (N32I)

POLR1C (N74S)

chr6.tRNA2-MetCAT

IGS

rRNA gene

Figure 4 | Impact of POLR1C mutations on polymerase association with chromatin. (a) ChIP-Seq experiments of FLAG-tagged POLR1C variants (wild type, N32I or N74S). Aggregate prole produced with the annotation mode of the Versatile Aggregate Proler shows ChIP-Seq data sets over the three classes of POLR3-transcribed genes, as dened by the promoter structure. Type 1 genes have an internal promoter composed of A and C boxes (5S rRNA genes). Type 2 genes have an internal promoter composed of A and B boxes (tRNA genes, for example), while the promoter of type 3 genes (U6, 7SK, RNase P and others) is located upstream of the transcription start site (TSS)55. The TSS was used as the reference point. (b) IGV view of a tRNA-Met gene transcribed by POLR3. POLR1C binding is decreased in mutated variants compared with wild type. IGS, intergenic spacer. (c) IGV view of one ribosomal DNA (rDNA) repeat. The rDNA gene encodes a 45S pre-rRNA precursor that will generate the 5.8S, 18S and 28S rRNAs. There are B400 copies of the rDNA gene arranged in tandem repeats in the human genome. rDNA repeats are not present in the reference genome assemblies; therefore, unique reads were aligned directly to the human rDNA reference sequence (NCBI accession number: HSU13369)56. No differences were observed in POLR1C occupancy between wild-type and mutant variants. A schematic of a rDNA repeat is included below the graph.

Methods

Patients and exome sequencing. Informed consent was obtained from all participants. The project was approved by the research ethics committee of the Montreal Children Hospital (11-105-PED), the institutional review board ofthe VU University Medical Center, Neuroscience Campus, Amsterdam, the Netherlands, the Childrens National Health System as part of the Myelin Disorders Bioregistry Project in Washington, DC and the University of Queensland, Australia. Genomic DNA was extracted from peripheral blood leukocytes of patients and family members using the Qiagen Gentra Puregene Blood Kit (Qiagen, Hilden, Germany) according to the manufacturers instructions. Exome sequencing was performed in two cases by PerkinElmer (Branford, Connecticut, USA) and in a third case by the Institute for Molecular Biology at the University of Queensland. Cases one and two were sequenced using PerkinElmers sequencing service using the Agilent Sure Select Human All Exon Capture V4 Kit and exome sequencing for these two cases was performed (two paired-end 100-bp reads) with the Illumina HiSeq 2000 system. Reads were aligned to the reference human genome (UCSC Genome Browser hg19) with the Genome Analysis Toolkit (GATK)32,33, SAMtools32,34, Picard (see web resource) and CASAVA v1.8 (ref. 35) and annotated using the snpEff software tool (http://snpeff.sourceforge.net/

Web End =http://snpeff.sourceforge.net/), as well as visualization tools from PerkinElmer. Subsequent analyses were performed using the Ingenuity software package (Qiagen, Redwood City, USA). Exome enrichment for case three was performed using the Illumina Nextera Rapid Capture kit and sequenced on an Illumina HiSeq 2000 (2 100 bp paired-end).

Reads were aligned to the reference human genome (GRCh37) and pedigree-informed variant calling was performed using the Real Time Genomics (RTG) integrated analysis tool rtg Family v3.2 (ref. 36). Variants were annotated using SnpEff v3.4 (ref. 37). In all cases, ltering queries were created as specic presets that allowed in silico reduction of variant lists down to candidates with correlation to phenotype, transmission mode of inheritance and alteration classication of pathogenicity. Sanger sequencing and co-segregation analysis were performed on genomic DNA using primer pairs designed with the primer3 software package and the genomic sequence of POLR1C (NM_203290; GRCh37/hg19). PCR products were forward- and reverse-sequenced at the McGill University and Genome Quebec Innovation Centre using an ABI 3730xl DNA Analyzer (ABI; Applied Biosystems, Foster City, CA, USA). Sequences were analysed using SeqMan 4.03 (DNAStar, Wisconsin, USA) and Chromas 1.62 (Technelysium Pty, Ltd, Australia).

Immunouorescence and western blot analyses. Transfection experiments for generating stable HeLa cell lines expressing FLAG-tagged versions of POLR1C variants used lipofectamine, as described by the supplier (Invitrogen, Carlsbad, CA, USA)18. The antibodies used in this study were obtained from the following sources: anti-FLAG monoclonal primary antibody (Sigma-Aldrich, St Louis,

Missouri, USA), Alexa Fluor 488 uorescence secondary antibody (Invitrogen) and anti-GAPDH (FL-335; Santa Cruz Biotechnology, Santa Cruz, CA, USA). Immunouorescence studies used an anti-FLAG (dilution 1/300) antibody and the secondary antibody Alexa Fluor 488 (1/200) to localize exogenously expressed FLAG-POLR1C variants in HeLa cells. For western blot analysis, anti-FLAG (dilution 1/3,000) was used to detect the FLAG-POLR1C variants and anti-GAPDH (dilution 1/2,000;) for loading control18.

Protein afnity purication coupled to mass spectrometry. Generation of cell lines stably expressing FLAG-tagged POLR1C subunits (WT and mutated) and afnity purication from the soluble fraction were performed using standard procedures38,39. The eluates were digested with trypsin and the resulting tryptic peptides were puried and identied with tandem mass spectrometry (LC-MS/MS) using a microcapillary reversed-phase high-pressure liquid chromatography-coupled LTQ-Orbitrap (ThermoElectron) quadrupole ion trap mass spectrometer with a nanospray interface, as recently described40. Protein database searching and protein spectral count quantication were performed with Mascot (version2.3.02)41. The NCBI_Human protein sequence database was downloaded on 20 February 2014. Known protein contaminants such as keratins, which are not expressed in HeLa cells, were excluded from the data set. Undistinguishable protein isoforms were considered as a single protein. For each LC-MS/MS analysis, protein spectral counts were normalized by the spectral count of the FLAG-POLR1C in order to allow the comparison of different purications. To simulate the background noise of the LC-MS/MS analysis of a given sample, spectral counts reported as 0 by Mascot were replaced by randomly generated spectral count values that are normally distributed with a mean and s.d. equal to those of the lowest 20% spectral count values from the LC-MS/MS analysis. Each replicate LC-MS/MS analysis of the afnity purications of the FLAG-POLR1C mutants (MUT; N32I, N74S and R279Q) was paired with a LC-MS/MS analysis of the WT FLAG-POLR1C that was performed at the same time. The set of high-condence interactors of POLR1C for a given mutant analysis was identied by comparing the spectral counts of the interactors obtained from the purications of the paired WT POLR1C and the MUT POLR1C to those of the proteins puried with an empty vector (EV) of the FLAG tag (nonspecic interaction). A protein is labelled as a high-condence interactor if it was identied and quantied in all three replicates of WT POLR1C or MUT POLR1C and that the ratio of the average spectral counts across the three replicates (WT/EV or MUT/EV) was greater than 5. These stringent criteria allow us to eliminate the vast majority of nonspecic interactors of POLR1C for the analysis of the interactions of each mutant. For each high-condence interacting protein, a two-tailed one-sample t-test was performed on the spectral count ratios (MUT/WT). The resulting P values were adjusted for multiple hypothesis testing using the BenjaminiHochberg procedure42. To maximize the

6 NATURE COMMUNICATIONS | 6:7623 | DOI: 10.1038/ncomms8623 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

& 2015 Macmillan Publishers Limited. All rights reserved.

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms8623 ARTICLE

a

WT R279Q

b

Ratio<0.67 Ratio>1.5

TUBA4A

POLR1A

POLR1B

POLR1E

CD3EAP

ZNRD1

POLR3A

POLR3B

POLR3C

POLR3D

POLR3E

POLR3F

POLR3G

POLR3H

CRCP *POLR1C

POLR1D

POLR2E

POLR2F

POLR2H

POLR2K

POLR2L

Log 10(t-test adjusted Pvalue)

3.5

3

2.5

2

1.5

1

0.5

0 4.5 3 1.5 0 1.5 3 4.5 Log2(average ratio R279Q / WT)

FLAG DNA Merge

POLR1C (WT)

POLR1C (R279Q)

Specific Pol l subunits

Specific Pol lll subunits

Shared subunits (Pol l and

pol lll)

Pvalue < 0.05

ROCK2

c

4

3.5

3

2.5

2

1.5

1

0.5

0

0.5

1

20 m

Figure 5 | Impact of a TCS-causative mutation in POLR1C on polymerase assembly and cellular localization. FLAG-tagged POLR1C variants, either the wild type (1C) or a mutated version with the R279Q substitution, were expressed, afnity-puried and used in anti-FLAG immunouorescence experiments as in Fig. 3. (a) Afnity purication coupled to mass spectrometry data is represented in the form of a heatmap that contains the log2-transformed average spectral count ratios R279Q/WTacross all three replicates. Spectral counts were computed with Mascot (see legend to Fig. 3 for details). (b) Volcano plot of the log2-transformed average spectral count ratios N279Q/WT (x axis) and the log10-transformed P-values (adjusted with the BenjaminiHochberg procedure) resulting from the two-tailed one-sample t-tests of the high-condence interactors of POLR1C. Red proteinsshow a level of differential interaction with POLR1C that is statistically signicant, while blue proteins do not. Proteins with a log2-transformed average spectral count ratio o 4.5 were capped to 4.5 for display purposes. (c) Immunouorescence data showing the cellular localization of tagged POLR1C

variants are shown. Scale bar, 20 mm.

specicity of our approach, a protein is deemed to show a level of differential interaction with POLR1C that is statistically signicant when its adjusted P value is o0.05 and that its average spectral count fold-change (MUT/WT) is o0. 6or 41.5.

ChIP-sequencing and data analysis. Stable HeLa cell lines expressing FLAG-tagged POLR1C (WT or mutated) were cultured to 80% conuence and cross-linked with 1% formaldehyde directly in the cell medium for 5 min followed by a 5-min quenching in 125 mM glycine. For ChIP experiments43, nuclei from 3 106

cells were lysed and re-suspended in sonication buffer (10 mM Tris-HCl pH 8, 140 mM NaCl, 1 mM EDTA, 0.5 mM EGTA, 0.5% Triton, 0.5% SDS and protease inhibitors). Chromatin was prepared by sonicating the nuclei in a Covaris E-Series E22 for 6 min at Duty 2. This generated chromatin fragments of 500 bp on average. Sonicated chromatin was immunoprecipitated using 25 ml of anti-FLAG M2

Magnetic beads (Sigma) for 4 h. The beads were then washed and eluted, followed by phenol:chloroform-isoamyl (25:24:1 pH 8, Invitrogen) extraction and ethanol precipitation. Sequencing libraries were prepared from input and ChIP eluates (WT and mutated) using the TruSeq DNA library preparation kit (Illumina). Libraries were sequenced on an Illumina HiSeq 2000 (2 50 cycles paired-end).

Quality control of the sequencing data was performed with FastQC and low-quality bases and adapter sequences were trimmed with Trimmomatic44. Unique reads were aligned to the reference human genome (build hg19) or to the human rDNA reference sequence (rDNA; NCBI accession number: HSU13369) with Bowtie version 2.1.0 (ref. 45). Input DNA signal was subtracted from ChIP signal. POLR1C peaks were called using MACS version 2.0.10 (ref. 46). Peaks were annotated with HOMER version 4.5.0 (ref. 47) using GENCODE Genes V19 annotations48 and the tRNA Genes track of the UCSC Genome Browser49,50. Data visualization was performed with the Integrative Genomics Viewer51. To compare chromatin binding between WT and mutated POLR1C variants over all POLR3-transcribed genes divided into three classes according to their regulatory elements

(see legend to Fig. 4), we used the annotation mode of the Versatile Aggregate Proler with the transcription start site as the only reference point and 100 windows of 10 base pairs each on both sides of the reference point52.

References

1. Schiffmann, R. & van der Knaap, M. S. The latest on leukodystrophies. Curr. Opin. Neurol. 17, 187192 (2004).

2. Schiffmann, R. & van der Knaap, M. S. Invited article: an MRI-based approach to the diagnosis of white matter disorders. Neurology 72, 750759 (2009).3. Bernard, G. & Vanderver, A. in GeneReviewss (eds Pagon, R. A., Adam, M. P., Ardinger, H. H., et al.) (University of Washington, Seattle, 19932015), http://www.ncbi.nlm.nih.gov/books/NBK99167/

Web End =http:// http://www.ncbi.nlm.nih.gov/books/NBK99167/

Web End =www.ncbi.nlm.nih.gov/books/NBK99167/ .

4. Bernard, G. et al. Mutations of POLR3A encoding a catalytic subunit of RNA polymerase Pol III cause a recessive hypomyelinating leukodystrophy. Am. J. Hum. Genet. 89, 415423 (2011).

5. Tetreault, M. et al. Recessive mutations in POLR3B, encoding the second largest subunit of Pol III, cause a rare hypomyelinating leukodystrophy. Am. J. Hum. Genet. 89, 652655 (2011).

6. Saitsu, H. et al. Mutations in POLR3A and POLR3B encoding RNA polymerase III subunits cause an autosomal-recessive hypomyelinating leukoencephalopathy. Am. J. Hum. Genet. 89, 644651 (2011).

7. Terao, Y. et al. Diffuse central hypomyelination presenting as 4H syndrome caused by compound heterozygous mutations in POLR3A encoding the catalytic subunit of polymerase III. J. Neurol. Sci. 320, 102105 (2012).

8. Potic, A., Brais, B., Choquet, K., Schiffmann, R. & Bernard, G. 4H syndrome with late-onset growth hormone deciency caused by POLR3A mutations. Arch. Neurol. 69, 920923 (2012).

NATURE COMMUNICATIONS | 6:7623 | DOI: 10.1038/ncomms8623 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 7

& 2015 Macmillan Publishers Limited. All rights reserved.

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms8623

9. Daoud, H. et al. Mutations in POLR3A and POLR3B are a major cause of hypomyelinating leukodystrophies with or without dental abnormalities and/or hypogonadotropic hypogonadism. J. Med. Genet. 50, 194197 (2013).

10. Wolf, N. I. et al. Clinical spectrum of 4H leukodystrophy caused by POLR3A and POLR3B mutations. Neurology 83, 18981905 (2014).

11. Steenweg, M. E. et al. Magnetic resonance imaging pattern recognition in hypomyelinating disorders. Brain 133, 29712982 (2010).

12. La Piana, R. et al. Brain magnetic resonance imaging (MRI) pattern recognition in Pol III-related leukodystrophies. J. Child Neurol. 29, 214220 (2014).

13. Dumay-Odelot, H., Durrieu-Gaillard, S., Da, S. D., Roeder, R. G. & Teichmann,M. Cell growth- and differentiation-dependent regulation of RNA polymerase III transcription. Cell Cycle 9, 36873699 (2010).14. Dauwerse, J. G. et al. Mutations in genes encoding subunits of RNA polymerases I and III cause Treacher Collins syndrome. Nat. Genet. 43, 2022 (2011).

15. Miron-Garcia, M. C. et al. The prefoldin bud27 mediates the assembly of the eukaryotic RNA polymerases in an rpb5-dependent manner. PLoS Genet. 9, e1003297 (2013).

16. Jeronimo, C. et al. Systematic analysis of the protein interaction network for the human transcription machinery reveals the identity of the 7SK capping enzyme. Mol. Cell 27, 262274 (2007).

17. Forget, D. et al. Nuclear import of RNA polymerase II is coupled with nucleocytoplasmic shuttling of the RNA polymerase II-associated protein 2. Nucleic Acids Res. 41, 68816891 (2013).

18. Forget, D. et al. The protein interaction network of the human transcription machinery reveals a role for the conserved GTPase RPAP4/GPN1 and microtubule assembly in nuclear import and biogenesis of RNA polymerase II. Mol. Cell. Proteom. 9, 28272839 (2010).

19. Forget, D., Cloutier, P., Domecq, C. & Coulombe, B. in Systems Analysis of Chromatin-Related Protein Complexes in Cancer. (eds Emili, Andrew, Greenblatt, Jack & Wodak, ShoshanaCh. 12, 227238 (Springer, 2014).

20. Cloutier, P. et al. High-resolution mapping of the protein interaction network for the human transcription machinery and afnity purication of RNA polymerase II-associated complexes. Methods 48, 381386 (2009).

21. Boulon, S. et al. HSP90 and its R2TP/Prefoldin-like cochaperone are involved in the cytoplasmic assembly of RNA polymerase II. Mol. Cell 39, 912924 (2010).

22. Canella, D., Praz, V., Reina, J. H., Cousin, P. & Hernandez, N. Dening the RNA polymerase III transcriptome: genome-wide localization of the RNA polymerase III transcription machinery in human cells. Genome Res. 20, 710721 (2010).

23. Kutter, C. et al. Pol III binding in six mammals shows conservation among amino acid isotypes despite divergence among tRNA genes. Nat. Genet. 43, 948955 (2011).

24. Raha, D. et al. Close association of RNA polymerase II and many transcription factors with Pol III genes. Proc. Natl Acad. Sci. USA 107, 36393644 (2010).

25. The Treacher Collins Syndrome Collaborative Group et al. Positional cloning of a gene involved in the pathogenesis of Treacher Collins syndrome. Nat. Genet. 12, 130136 (1996).

26. Scheper, G. C. et al. Mitochondrial aspartyl-tRNA synthetase deciency causes leukoencephalopathy with brain stem and spinal cord involvement and lactate elevation. Nat. Genet. 39, 534539 (2007).

27. Steenweg, M. E. et al. Leukoencephalopathy with thalamus and brainstem involvement and high lactate LTBL caused by EARS2 mutations. Brain 135, 13871394 (2012).

28. Gotz, A. et al. Exome sequencing identies mitochondrial alanyl-tRNA synthetase mutations in infantile mitochondrial cardiomyopathy. Am. J. Hum. Genet. 88, 635642 (2011).

29. Taft, R. J. et al. Mutations in DARS cause hypomyelination with brain stem and spinal cord involvement and leg spasticity. Am. J. Hum. Genet. 92, 774780 (2013).

30. Wolf, N. I. et al. Mutations in RARS cause hypomyelination. Ann. Neurol. 76, 134139 (2014).

31. Bilguvar, K. et al. Whole-exome sequencing identies recessive WDR62 mutations in severe brain malformations. Nature 467, 207210 (2010).

32. DePristo, M. A. et al. A framework for variation discovery and genotyping using next-generation DNA sequencing data. Nat. Genet. 43, 491498 (2011).

33. McKenna, A. et al. The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 20, 12971303 (2010).

34. Li, H. et al. The sequence alignment/map format and SAMtools. Bioinformatics 25, 20782079 (2009).

35. Cheng, A. Y., Teo, Y. Y. & Ong, R. T. Assessing single nucleotide variant detection and genotype calling on whole-genome sequenced individuals. Bioinformatics 30, 17071713 (2014).

36. Cleary, J. G. et al. Joint variant and de novo mutation identication on pedigrees from high-throughput sequencing data. J. Comput. Biol. 21, 405419 (2014).

37. Cingolani, P. et al. A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: SNPs in the genome of Drosophila melanogaster strain w1118; iso-2; iso-3. Fly 6, 8092 (2012).

38. Mellacheruvu, D. et al. The CRAPome: a contaminant repository for afnity purication-mass spectrometry data. Nat. Methods 10, 730736 (2013).

39. Chen, G. I. & Gingras, A. C. Afnity-purication mass spectrometry (AP-MS) of serine/threonine phosphatases. Methods 42, 298305 (2007).

40. Lavallee-Adam, M. et al. Discovery of cell compartment specic protein-protein interactions using afnity purication combined with tandem mass spectrometry. J. Proteome Res. 12, 272281 (2013).

41. Perkins, D. N., Pappin, D. J., Creasy, D. M. & Cottrell, J. S. Probability-based protein identication by searching sequence databases using mass spectrometry data. Electrophoresis 20, 35513567 (1999).

42. Benjamini, Y. & Hochberg, Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J. R. Stat. SocB Methodol.. 57, 289300 (1995).

43. Langlais, D., Couture, C., Sylvain-Drolet, G. & Drouin, J. A pituitary-specic enhancer of the POMC gene with preferential activity in corticotrope cells. Mol. Endocrinol. 25, 348359 (2011).

44. Lohse, M. et al. RobiNA: a user-friendly, integrated software solution for RNA-Seq-based transcriptomics. Nucleic Acids Res. 40, W622W627 (2012).

45. Langmead, B. & Salzberg, S. L. Fast gapped-read alignment with Bowtie 2. Nat. Methods 9, 357359 (2012).

46. Zhang, Y. et al. Model-based analysis of ChIP-Seq (MACS). Genome Biol. 9, R137 (2008).

47. Heinz, S. et al. Simple combinations of lineage-determining transcription factors prime cis-regulatory elements required for macrophage and B cell identities. Mol. cell 38, 576589 (2010).

48. Harrow, J. et al. GENCODE: the reference human genome annotation for The ENCODE Project. Genome Res. 22, 17601774 (2012).

49. Lowe, T. M. & Eddy, S. R. tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 25, 955964 (1997).

50. Chan, P. P. & Lowe, T. M. GtRNAdb: a database of transfer RNA genes detected in genomic sequence. Nucleic Acids Res. 37, D93D97 (2009).51. Thorvaldsdottir, H., Robinson, J. T. & Mesirov, J. P. Integrative Genomics Viewer (IGV): high-performance genomics data visualization and exploration. Brief. Bioinformatics 14, 178192 (2013).

52. Coulombe, C. et al. VAP: a versatile aggregate proler for efcient genome-wide data representation and discovery. Nucleic Acids Res. 42, W485W493 (2014).

53. Fernandez-Tornero, C. et al. Crystal structure of the 14-subunit RNA polymerase I. Nature 502, 644649 (2013).

54. Wild, T. & Cramer, P. Biogenesis of multisubunit RNA polymerases. Trends Biochem. Sci. 37, 99105 (2012).

55. Alla, R. K. & Cairns, B. R. RNA polymerase III transcriptomes in human embryonic stem cells and induced pluripotent stem cells, and relationships with pluripotency transcription factors. PLoS ONE 9, e85648 (2014).

56. Cong, R. et al. Interaction of nucleolin with ribosomal RNA genes and its role in RNA polymerase I transcription. Nucleic Acids Res. 40, 94419454 (2012).

Acknowledgements

We are grateful for the support and cooperation of the subjects, as well as their families. G.B. has received a Research Scholar Junior 1 award from the Fonds de Recherche du Qubec en Sant (FRQS). I.T. has received a scholarship from the Rseau de mdecine gntique applique du FRQS. N.I.W. and M.S.v.d.K. are supported by ZonMw (TOP grant 91211005). B.C. holds the Bell-Bombardier Research Chair of Excellence (IRCM). K.C. received a scholarship from the FRQS. D.T. received support from the German Research Council and European Union (Marie Curie ITN). R.J.T., C.S. and A.V. are supported by a National Health and Medical Research Council, Australia project grant (APP1068278). The Myelin Disorders Bioregistry Project was supported by the NIH National Center for Advancing Translational Sciences (award number UL1TR000075). S.F. was supported by the Ligue Rgionale contre le Cancer (comit Dordogne), INSERM and the Conseil Rgional dAquitaine. This work was supported by Canadian Institute for Health Research grant MOP-G-287547 to G.B., B.B. and B.C. We also thank the McGill University and Genome Quebec Innovation Center, the IRCM Molecular Biology and Functional Genomics Platform and Denis Faubert, Jose Champagne, Sylvain Tessier and Marguerite Boulos of the IRCM Proteomics Discovery Platform for their services and Dr Franois Robert and Alexis Blanchet-Cohen for bioinformatics support. We thank Dr Claudia Kleinman for her advice regarding ChIP-Seq data analysis. J.R.Y. has been funded by National Institutes of Health (NIH) grants P41 GM103533, R01 MH067880, R01 MH100175 and UCLA/NHLBI Proteomics Centers (HHSN268201000035C). M.L.-A. holds a postdoctoral fellowship from Fonds de Recherche du QubecNature et Technologies.

Author contributions

I.T. and N.I.W. collected and analysed the sequencing and clinical/radiological data and wrote the rst draft of the manuscript. D.F. performed the experiments, collected and

8 NATURE COMMUNICATIONS | 6:7623 | DOI: 10.1038/ncomms8623 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

& 2015 Macmillan Publishers Limited. All rights reserved.

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms8623 ARTICLE

analysed the data and wrote the rst draft of the manuscript. K.G. performed the experiments, collected and analysed the data and reviewed the manuscript. L.T.T. and K.C. collected and analysed the data and reviewed the manuscript. N.I.W., G.B., B.B., G.Y., L.S., R.I.W., D.T., B.P.v.d.W., J.S., A.Z., A.M., C.G., A.V., E.F. and M.S.v.d.K. recruited the patients, contributed to the phenotyping of the patients and reviewed the manuscript. A.V., R.J.T. and C.S. provided exome data on one of the patients and critically reviewed the manuscript. S.F. worked on the three-dimensional representation of the different POLR1C mutations and the predictions of their impact on the RNA polymerases I and III, and reviewed the manuscript. M.L.-A. and J.R.Y. performed the statistical analyses of LC-MS/MS data, in collaboration with C.P., and reviewed the manuscript. B.C. designed and supervised the functional studies and reviewed the manuscript. G.B. designed and supervised the study, and reviewed the manuscript.

Additional information

Accession codes: The high-throughput sequencing data generated in this study have been deposited into the NCBI Sequence Read Archive under the accession code SRP057978.

Supplementary Information accompanies this paper at http://www.nature.com/naturecommunications

Web End =http://www.nature.com/ http://www.nature.com/naturecommunications

Web End =naturecommunications

Competing nancial interests: The authors declare no competing nancial interests.

Reprints and permission information is available online at http://npg.nature.com/reprintsandpermissions/

Web End =http://npg.nature.com/ http://npg.nature.com/reprintsandpermissions/

Web End =reprintsandpermissions/

How to cite this article: Thiffault, I. et al. Recessive mutations in POLR1C cause a leukodystrophy by impairing biogenesis of RNA polymerase III. Nat. Commun. 6:7623 doi: 10.1038/ncomms8623 (2015).

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the articles Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/

Web End =http://creativecommons.org/licenses/by/4.0/

NATURE COMMUNICATIONS | 6:7623 | DOI: 10.1038/ncomms8623 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 9

& 2015 Macmillan Publishers Limited. All rights reserved.