
==== Front
Cell Death Dis
Cell Death Dis
Cell Death & Disease
2041-4889
Nature Publishing Group UK London

39251572
7042
10.1038/s41419-024-07042-5
Article
Cilia defects upon loss of WDR4 are linked to proteasomal hyperactivity and ubiquitin shortage
Burkhalter Martin D. 1
Stiff Tom 2
Maerz Lars D. 3
Casar Tena Teresa 3
Wiese Heike 4
Gerhards Julian 1
http://orcid.org/0000-0002-3804-0986
Sailer Steffen A. 1
Vu Linh Anna Trúc 1
Duong Phu Max 1
Donow Cornelia 3
Alupei Marius 5
http://orcid.org/0000-0002-5207-2861
Iben Sebastian 5
http://orcid.org/0000-0002-9199-8990
Groth Marco 6
http://orcid.org/0000-0001-5697-2608
Wiese Sebastian 4
Church Joseph A. 7
http://orcid.org/0000-0002-8481-1462
Jeggo Penelope A. 2
http://orcid.org/0000-0003-2714-965X
Philipp Melanie melanie.philipp@uni-tuebingen.de

1
1 https://ror.org/03a1kwz48 grid.10392.39 0000 0001 2190 1447 Department of Experimental and Clinical Pharmacology and Pharmacogenomics, Section of Pharmacogenomics, Eberhard-Karls-University Tübingen, 72074 Tübingen, Germany
2 https://ror.org/00ayhx656 grid.12082.39 0000 0004 1936 7590 Genome Damage and Stability Centre, University of Sussex, Brighton, BN1 9RQ UK
3 https://ror.org/032000t02 grid.6582.9 0000 0004 1936 9748 Institute for Biochemistry and Molecular Biology, Ulm University, 89081 Ulm, Germany
4 https://ror.org/032000t02 grid.6582.9 0000 0004 1936 9748 Core Unit Mass Spectrometry and Proteomics, Medical Faculty, Ulm University, 89081 Ulm, Germany
5 https://ror.org/032000t02 grid.6582.9 0000 0004 1936 9748 Department of Dermatology, Ulm University, 89081 Ulm, Germany
6 https://ror.org/039a53269 grid.418245.e 0000 0000 9999 5706 Core Facility Next Generation Sequencing, Leibniz Institute on Aging–Fritz Lipmann Institute, 07745 Jena, Germany
7 https://ror.org/00412ts95 grid.239546.f 0000 0001 2153 6013 Department of Pediatrics, Keck School of Medicine of University of Southern California, Children’s Hospital Los Angeles, Los Angeles, CA 90033 USA
9 9 2024
9 9 2024
9 2024
15 9 6603 2 2024
7 8 2024
29 8 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by/4.0/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
The WD repeat-containing protein 4 (WDR4) has repeatedly been associated with primary microcephaly, a condition of impaired brain and skull growth. Often, faulty centrosomes cause microcephaly, yet aberrant cilia may also be involved. Here, we show using a combination of approaches in human fibroblasts, zebrafish embryos and patient-derived cells that WDR4 facilitates cilium formation. Molecularly, we associated WDR4 loss-of-function with increased protein synthesis and concomitant upregulation of proteasomal activity, while ubiquitin precursor pools are reduced. Inhibition of proteasomal activity as well as supplementation with free ubiquitin restored normal ciliogenesis. Proteasome inhibition ameliorated microcephaly phenotypes. Thus, we propose that WDR4 loss-of-function impairs head growth and neurogenesis via aberrant cilia formation, initially caused by disturbed protein and ubiquitin homeostasis.

Subject terms

Ciliogenesis
Disease model
https://doi.org/10.13039/501100001659 Deutsche Forschungsgemeinschaft (German Research Foundation) PH144/4-1 PH144/4-3 PH144/6-1 467868420 SFB 1074 SFB 1074 Wiese Heike Wiese Sebastian Philipp Melanie issue-copyright-statement© Associazione Differenziamento e Morte Cellulare ADMC 2024
==== Body
pmcIntroduction

Primary microcephaly (PM) is a failure in the growth of the brain and skull, which can be diagnosed from 32 weeks of gestation using ultrasound. Microcephaly is defined by reduced head circumference by at least 2 standard deviations below the median. Such reduction in head circumference directly correlates with a smaller size of an underdeveloped brain. As a consequence, infants with PM suffer from a number of symptoms including cognitive and motoric impairments, epilepsy, cerebral palsy, a delay in speaking and defects of the eyes [1]. Some PM children also show overall growth reduction such as in the case of Seckel or Meier-Gorlin Syndromes, which is then called microcephalic primordial dwarfism [2].

Despite the rare occurrence of the condition, several genetic loci have been linked to PM. We and others have shown that proteins required for DNA replication (i.e., ORC1), DNA damage response (i.e., ATR, ATRIP) and DNA repair (i.e., PRKDC, XRCC4, DNA ligase IV) allow for faithful brain growth [3–10]. Interestingly, in addition to their canonical function, namely maintenance of intact DNA and cell division, the formation of aberrant centrosomes can be observed upon loss of ORC1 or ATR [11, 12]. In addition, centrosomal aberrations can provoke PM, too, through reduced self-renewal and premature differentiation of neuronal progenitors, altered duration of the cell cycle, as well as a delay in mitosis [13]. Centrosomes, however, provide also the basal bodies for cilia, which are crucial cell organelles involved in cellular signal sensing, signal transduction and cell cycle progression [14]. Thus, microcephaly may arise also in the context of impaired cilia formation or function and therefore could be a disorder with ciliary contribution [11, 12, 15]. Additionally, genes with as yet unknown or less investigated functions have been identified as drivers of microcephaly when mutated. One such gene is WDR4, which encodes for WD repeat-containing protein 4 (also known as TRM82, tRNA methyltransferase 82 homolog or protein Wuho homolog) [16]. WDR4 is a WD40 domain containing protein, which is best known for its interaction with the methyltransferase METTL1 [17–19]. WDR4 promotes activity of METTL1 to generate N7-methylguanosine (m7G) tRNAs at position 46 to stabilize tRNAs, that has been shown to be of importance for maintenance and differentiation of embryonic stem cells [20]. The Drosophila homolog to WDR4, Wuho, has been linked to impaired fertility defects, perhaps due to impaired genomic stability [21]. Increased levels of WDR4 and its interaction partner METTL1 have been observed in various tumors suggesting that these two factors together increase cancer susceptibility [22–24]. Mechanistically, METTL1/WDR4 have been connected to different pathways to increase carcinogenesis, such as the WNT pathway as well as autophagy [25, 26]. However, Wang and colleagues elegantly showed that overexpression of WDR4 alone in A459 cells is sufficient to induce metastatic nodules in the lung upon transplantation of these cells [27]. In this case, WDR4 works as a substrate adapter of the Cullin 4 E3 ubiquitin ligase complex, which specifically targets the tumor suppressor promyelocytic leukemia (PML) for degradation [27].

Genetic variations of WDR4 have been repeatedly reported in the context of microcephalic disorders such as Galloway-Mowat syndrome 6 (OMIM 618347) and microcephaly, growth deficiency seizures and brain malformations (OMIM 618346). Several genetic variants have been identified in affected individuals all over the world suggesting that WDR4 is an essential protein for normal brain and skull development [28–31]. Indeed, the Chen lab recently showed in mice that WDR4 supports development of the cerebellum by activating the small GTPase RAC1 through degradation of Arhgap17, a negative regulator of RAC1 [32]. This prevents premature cell cycle exit and differentiation of progenitor cells needed for cerebellum development [32]. However, since PM induced by WDR4 dysfunction affects the whole brain, we wondered whether other mechanisms contribute to the phenotypes observed in patients [28–31]. We here show that WDR4 is needed to properly generate cilia, organelles important for cell signaling and brain development. Our data suggest that this function of WDR4 can be separated from its canonical interaction partner METTL1. Instead, dysregulated protein homeostasis involving changing ubiquitin distribution presents as the most likely driver for disturbed ciliogenesis in absence of WDR4.

Results

Despite WDR4 having been associated to PM, the mechanism how this factor affects embryonal development so far has not been completely understood. Because of its many advantages including rapid development of translucent embryos as well as high conservation to mammalian species we turned to zebrafish as model organism, in which a homolog of human WDR4 exists. In situ hybridization experiments showed that wdr4 is expressed throughout development with enrichment in neural structures such as the developing eyes, brain and neural tube (Fig. 1a). Consistent with reports in patients and mice [28–32] and the expression pattern in zebrafish, we found that depletion of Wdr4 using an antisense morpholino oligonucleotide (MO)-mediated approach resulted in zebrafish with smaller heads and eyes (Figs. 1b–d and S1a, b). Two different MOs were tested, one interfering with translation of wdr4 mRNA into protein and one, which prevents regular splicing (Fig. S2). Both MOs provoked similar defects in the formation of anterior structures. To further characterize this defect, we performed Alcian blue staining of 4 days post fertilization (4 dpf) embryos and found that formation of cartilaginous structures was still occurring, however to a lesser extent with structures appearing shortened compared to controls (Fig. 1e). This is consistent with typical PM phenotypes.Fig. 1 Wdr4 is expressed in and required for the developing brain.

a Whole mount in situ hybridization (WMISH) revealed expression of wdr4 in anterior structures of the zebrafish embryo (arrow heads) including the developing brain. Tb tailbud stage; 8 ss 8 somite stage; hpf hours post fertilization. Scale bars: 200 µm. b Live images of 48 hpf zebrafish embryos, which were either non-injected (NI) or injected with one of the following MOs: splCTRL, splMO or a MO blocking translation of Wdr4 (Wdr4 ATG MO). Scale bar: 250 µm. c Scatter plot showing the incidence of microcephaly in zebrafish upon Wdr4 KD. Co-injection of capped RNA encoding zebrafish Wdr4 rescues head size, while the patient-resembling Wdr4 R84H does not. n = 3–9 experiments with 93–244 embryos in total. d Loss of Wdr4 results in zebrafish embryos with smaller eyes, which cannot be fully rescued by wdr4 R84H injection. n = 3–9 experiments with 93–244 embryos in total. e Reduced size of cartilage structures in the developing head as shown by Alcian blue staining. Graph shows measurement of the palatoquadrate (red bracket). Lower bar shows control. n = 16/15; ****p < 0.0001, Mann–Whitney test. f Images of the otoliths in the otic vesicle upon injection of splCTRL or Wdr4 splMO. Scale bar: 10 µm. g Knockdown of Wdr4 increases the percentage of embryos with fused otoliths. n = 3–9 experiments with 93–244 embryos in total. ****p < 0.0001. One-way ANOVA with Sidak’s multiple comparison test. h Loss of Wdr4 randomizes heart looping in 48 hpf zebrafish embryos. D, correct D-loop; N, no loop; L, inverse loop. n = 4 experiments with 63–94 embryos. ****p < 0.0001. Fisher’s exact test. Scale bar: 200 µm. i Pancreas placement is randomized upon Wdr4 KD. R, pancreas right from the midline (dotted line); L, pancreas left from midline. Arrows indicate insulin n = 4 experiments with 63–94 embryos in total. ****p < 0.0001. Fisher’s exact test. Scale bar: 200 µm. Red lines: median.

In our quest for understanding the underpinnings of WDR4-associated PM better, we performed rescue experiments. To this end, we co-injected capped RNA encoding either wild-type Wdr4 or one encoding a mutant version, which mimics PM patients with a mutation in WDR4. Cells of these patients have been archived at the cell repository of the University of Sussex as the patients displayed a Seckel syndrome-like phenotype. The patients carry a variant of WDR4, which results in the exchange of arginine at position 84 with histidine (wdr4 R84H) leading to reduced expression of WDR4 (Fig. S3). This residue is conserved in zebrafish wdr4. Reconstitution of Wdr4 knockdown (KD) embryos with wild-type wdr4 at least partially restored growth of anterior structures (Fig. 1c, d) suggesting that the phenotypes after Wdr4 MO injection are specific. The same approach using wdr4 R84H, however, was unable to rescue head and eye size in zebrafish embryos (Fig. 1c, d), while injection of the RNAs alone had no effect (Fig. S1c, d). These results link the phenotypes in zebrafish embryos to PM in patients carrying genetic variants of WDR4.

Closer examination of Wdr4 KD embryos revealed additional phenotypes such as defects in otolith seeding (Figs. 1f, g and S1e, f). Otolith formation and positioning is facilitated through concerted beating of cilia within the otic vesicle resulting in the accumulation of calcium carbonate crystals on particular matrix proteins [33]. In case of cilia defects, the number of otoliths may be changed or, as observed here, otoliths become symmetrical and centered in the lower half of the otic placode (Fig. 1f). Reconstitution with wild-type Wdr4 significantly rescued the otolith defect, while the patient variant failed to do so (Fig. 1g). To test whether loss of Wdr4 causes also other defects in cilium function, we analyzed the situs of the injected embryos The situs of inner organs is determined early during embryogenesis through the development of an internal left-right asymmetry along the anterior-posterior axis. This is mediated by cilia in a temporal organ of laterality, which is called Kupffer’s vesicle (KV) in zebrafish [14]. When cilia do not form properly or become dysfunctional in the KV, left-right asymmetry will not be established properly and organs may end up on the wrong side of the body. Experimentally, this can be scored by heart looping direction and placement of abdominal organs. Wdr4 KD indeed produced more often embryos with inversely looped or even unlooped hearts than controls (Figs. 1h and S1g). Similarly, Wdr4 KD also randomized pancreas placement as observed by the expression of insulin, which marks the primary endocrine islet (Figs. 1i and S1h). These results suggest an unanticipated function of WDR4 in cilium biology.

Our results prompted us to more closely investigate cilia after Wdr4 depletion. We started with the analysis of the ciliated organ of laterality, the KV in zebrafish (Fig. 2a–c). Neither area nor cell number of the KV was altered (Fig. S4a–c). However, while we found only a subtle reduction in cilium length in the KV (Fig. 2b), we counted substantially fewer cilia per KV (Figs. 2c and S3d) suggesting a defect in ciliogenesis. To corroborate this finding, we also assessed cilia in other ciliated organs of the zebrafish. In the otic placode as well as the neural tube we counted considerably fewer cilia after Wdr4 KD (Fig. 2d–h). Moreover, cilia in the pronephric duct showed length shortening upon Wdr4 KD (Fig. 2i, j). Increased frequency of renal cysts in embryos with reduced Wdr4 suggests physiological relevance of this ciliogenesis defect (Fig. 2k, l).Fig. 2 Lack of Wdr4 affects cilia in multiple organs in zebrafish.

a Representative images of cilia in the Kupffer’s vesicle (KV) of 8 ss zebrafish embryos. Cilia in magenta (acetylated tubulin), green counterstain (PKCz) to visualize the KV outline. Scale bar: 10 µm. b KV cilia are shorter in the absence of Wdr4. n = 2 experiments with 14–17 KVs and 283–379 cilia in total. **p = 0.0066. Mann–Whitney test. c Zebrafish embryos depleted of Wdr4 extend fewer cilia into the KV. n = 2 experiments with 14–17 KVs in total. **p = 0.0034. Welch’s t test. Mean ± SEM. d Loss of Wdr4 reduces cilia in the otic placode. Staining as in (a). Scale bar: 25 µm. e Quantification of otic placode cilia in 3 and 5 embryos injected with CTRL MO or Wdr4 ATG MO, respectively. ***p = 0.0003, Welch’s t test. Mean SEM. Scale bar: 10 µm. f Visualization of cilia in the neural tube by detection of acetylated tubulin (magenta). Shown are typical sections of CTRL MO or WDR4 ATG morphants. Scale bar: 10 µm. Wdr4 KD reduces primary (g) and motile (h) cilia. g n = 6/5, **p = 0.0043, Mann–Whitney test; h n = 6/5, ****p < 0.0001, Welch’s t test. i Wdr4 KD affects cilia in the pronephros. Staining as in (a). Scale bar: 10 µm. j Pronephric cilia are shortened upon loss of Wdr4. n = 6 with 170 and 193 cilia, ****p < 0.0001, Mann–Whitney test. k Induction of glomerular cysts upon loss of Wdr4. Pictures show cyst-free and cystic (arrows) glomeruli using the tg(wt1b:eGFP) reporter line. Scale bar: 100 µm. l Quantification of glomerular cyst formation after injection of two different MOs. n = 3 with 69–115 embryos; CTRL MO vs. Wdr4 ATG MO ****p < 0.0001; splCTRL MO vs. Wdr4 splMO ****p < 0.0001. One-way ANOVA with Sidak correction. Red line: median.

We turned then to human fibroblasts (Fig. 3), in which cilium formation can be induced through prolonged serum starvation. We transfected cells with a pool of four siRNAs against WDR4 or respective control siRNAs and then serum-starved for 3 days. qPCR analysis showed successful KD of WDR4 after RNA interference (Fig. 3a). Immunofluorescence analysis of these cells revealed shortened cilia (Fig. 3b, c) and a significant reduction of cells extending cilia in general (Fig. 3d). Overexpression of zebrafish Wdr4 fully rescued cilia shortening in siWDR4 fibroblasts and increased slightly, yet significantly the number of ciliated cells (Fig. S5). We further analyzed primary fibroblasts derived from the Seckel syndrome-like patients carrying the R84H variant, which display a markedly reduced expression of WDR4 (Fig. S3). Similar to our KD experiments, we observed both, shorter as well as fewer cilia compared to control fibroblasts (Fig. 3e–g). Nucleofection of patient fibroblasts with a plasmid encoding human WDR4 rescued the number of cells extending cilia (Fig. 3h) and fully restored cilium length (Fig. 3i). We further found that CP110 and its interaction partner CEP97 accumulated more frequently at the base of cilia in patient cells than in control cells, suggesting a cilium assembly defect (Fig. 3j, k). Several studies have demonstrated that both proteins need to be removed from the mother centriole to allow for cilium extension. When CP110 and CEP97 persist at the mother centriole, cilium formation is impeded [34, 35]. Could therefore be a centrosome defect the underlying cause for altered cilium formation? To test this, we treated fibroblasts and zebrafish embryos with ciliobrevin D (CBD), which has been shown to deplete cells of functional cilia [36]. CBD treatment reduced ciliation of fibroblasts also in our hands, yet did not affect centrioles as judged by γ-Tubulin staining (Fig. S6). Treatment of zebrafish embryos with CBD, induced highly similar phenotypes as WDR4 KD, e.g., small heads, small eyes, and otolith defects (Fig. S6). Yet, combination of WDR4 KD with CBD treatment did not significantly worsen this phenotype. This suggests that loss of WDR4 seems to trigger its effect mostly through cilia defects and less by causing centrosome defects. To further test a possible impact of WDR4 on centrioles, we analyzed spindle poles in mitotic HEK 293T cells after KD of WDR4 and did not detect changes as compared to controls (Fig. S7). This again suggests that centrioles remain largely unaffected by loss of WDR4.Fig. 3 Cilia defects in the absence of WDR4 in human fibroblasts.

a siRNA transfection of human fibroblasts results in reduced WDR4 expression. qPCR normalized to housekeeping gene (SDHA) and to siCTRL transfected cells. n = 3. ***p = 0.0003, paired t-test. b Immunofluorescence of cilia in fibroblasts transfected with control or WDR4 siRNA. Magenta: acetylated tubulin (cilia); green: pericentrin (PCNT, centrioles). Scale bar: 5 µm. c Cilia are shorter upon WDR4 KD. n = 3 experiments with 95–99 cilia in total. ****p = 0.0001. Welch’s t test. d WDR4 KD cells have fewer cilia. n = 3 experiments. *p = 0.0108. Unpaired, Welch’s t test. e Immunofluorescence of control and patient-derived primary fibroblasts. Staining as in (a). Scale bar: 5 µm. f Wdr4 patient fibroblasts have fewer cilia. n = 4 experiments. *p = 0.0255. Welch’s t test. g WDR4 patient fibroblasts have shorter cilia. n = 4 experiments with 108–174 cilia in total. ****p < 0.0001. Mann–Whitney test. h Nucleofection of WDR4 patient fibroblasts with a plasmid encoding human WDR4 rescues the percentage of ciliated cells. n = 3 experiments. **p = 0.01. Welch’s t test. i Reconstitution with WDR4 restores cilium length in patient fibroblasts. n = 3 experiments with 99 and 100 cilia. ****p < 0.0001. Welch’s t test. j Wdr4 patient cells show increased frequency of aberrant CP110 staining (green). n = 4 experiments with 407 and 411 cells in total. *p = 0.0455. Paired t-test. Mean + SEM. Scale bar: 3 µm. k Wdr4 patient cells show increased frequency of aberrant CEP97 staining (green). n = 6 experiments with 625 and 617 cells in total. *p = 0.0342. Paired t-test. Mean + SEM. Scale bar: 3 µm. Red line indicates median.

In order to rule out that the differences between control and patient cells are due to changes in the cell cycle or leaky arrest in G0, we analyzed cyclin A2 levels in arrested cells by immunofluorescence and Western blotting (Fig. S8). Both approaches did not reveal any differences in G0-arrested cells. However, patient cells showed a non-significant delay in cell cycle entry upon release, which did not translate into a slowed cell cycle, since control and patient cells grow with similar kinetics (Fig. S8). Moreover, cilium resorption occurred at the same pace as in control cells (Fig. S8). Hence, dysfunction or loss of WDR4 likely prevents normal cilium extension.

WDR4 also could affect cilia via disturbed actin organization, since it activates the small GTPase RAC1 through degradation of Arhgap17 [32]. To test this, we analyzed stress fibers in control and patient fibroblasts. As we did not detect a difference, we consider actin organization to be largely unaffected by WDR4 loss in our hands (Fig. S9).

Length alterations are known to affect ciliary function [37]. The occurrence of renal glomerular cysts indeed suggested dysfunctionality of cilia upon Wdr4 KD (Fig. 2k, l). To more directly assess the ability of cilia to properly serve as signaling platforms, we assessed the formation of U-shaped somites in Wdr4 morphant embryos. This defect arises, if Hedgehog (Hh) signaling is altered, a pathway that relies on functional primary cilia [38]. Loss of Wdr4 induced U-shaped somite formation, what could be partially rescued by co-injection of wdr4 RNA, but not wdr4 R84H RNA (Fig. 4a, b). Furthermore, expression of Hh target genes gli1 and nkx2.2 was significantly reduced after Wdr4 KD (Fig. 4c, d). Additionally, WDR4 patient fibroblasts showed reduced capability to activate AKT upon stimulation with PDGF AA, a process that also depends on the cilium [39] (Fig. 4e, f). These data suggest that WDR4 enables cilium function.Fig. 4 Loss of WDR4 dampens ciliary signaling capacity.

a Embryos were left non-injected or were injected with either splCTRL MO or Wdr4 splMO and scored for U-shaped somites at the age of 28 hpf as a read-out for Hedgehog signaling. Only embryos injected with Wdr4 splMO show this defect. Lower panels show same embryos with higher magnification. Scale bars: upper panels 500 µm, lower panels 100 µm. b Co-injection of wdr4 RNA partially rescues induction of U-shaped somites upon loss of Wdr4. n = 3–10 experiments with 103–189 embryos. splCTRL vs. splMO **p = 0.0033; splMO vs. splMO + wdr4 RNA *p = 0.015; Kruskal–Wallis test with Dunn’s post-test. c Expression of gli1 in injected zebrafish embryos. Normalized to gapdh and splCTRL. n = 6; **p = 0.0022. Mann–Whitney test. d Expression of nkx2.2 in injected zebrafish embryos. Normalized to gapdh and splCTRL. n = 6; **p = 0.0079. Mann–Whitney test. e WDR4 patient cells show lower capacity for AKT activation upon stimulation with PDGF AA. f Quantification of AKT activation upon stimulation with PDGF AA. Signals were normalized to unphosphorylated AKT and nonstimulated. n = 3; **p = 0.0012. Welch’s t test. Red line: median.

WDR4 is generally known as a co-factor for the methyltransferase METTL1, which modifies several RNAs including tRNAs. tRNA modification in turn regulates the efficacy and precision of protein synthesis [20, 40]. We hence hypothesized that loss of WDR4 in patients or KD models could potentially cause insufficient and imprecise protein synthesis. Interestingly, however, measurements of de novo protein synthesis in control and patient derived fibroblasts showed a marked increase, which was even bigger in quiescent, ciliated cells (Fig. 5a). Yet, despite the increase in synthesis, overall protein content was not changed between control and patient-derived fibroblasts (Fig. 5b). Moreover, mass spectrometry analysis demonstrated that proteins were synthesized until the very C-terminal end suggesting that the translation machinery was still able to produce regular full-length proteins (Fig. 5c). A luciferase-based assay further demonstrated that WDR4 patient fibroblasts exhibit normal translation and folding fidelity (Fig. 5d). In line with this, we also did not find a higher frequency of amino acid exchanges using mass spectrometry (Fig. 5e) and concluded that despite increased translation, proteins were still synthesized largely without errors. We next applied an algorithm on our spectrometry data to predict post-translational protein modifications. This analysis revealed a substantial increase in protein ubiquitination in patient cells carrying the WDR4 R84H variant (Fig. 5f).Fig. 5 WDR4 patient fibroblasts are not deficient in protein synthesis.

a Nascent protein synthesis is increased in WDR4 patient fibroblasts under cycling and resting conditions. n = 4–5. p = 0.0038 (**) and 0.0002 (***). Welch’s t test. b Overall protein content is not changed in WDR4 patient fibroblasts. n = 5–6. Mann–Whitney test. c WDR4 patient cells do not have more premature abortions during protein synthesis. Bar graph shows the distribution of peptides detected by mass spectrometry along the whole length of the protein. n = 2. d Assessment of translation fidelity in control and WDR4 patient fibroblasts using a luciferase based assay. Cells were transfected with firefly luciferase reporter constructs including the indicated mutations in the active center of the luciferase. Elevated luciferase activity indicates erroneous translation and reduced translational fidelity. Firefly luciferase activity was normalized to co-transfected renilla luciferase activity. e Bar graph showing amino acid exchange rates in control and WDR4 patient fibroblasts. n = 2. f Bar graph depicting posttranslational modifications in control and WDR4 patient fibroblasts. Inset shows ubiquitination modification. n = 2.

Because of these unanticipated findings, we questioned the importance of the functional interaction of WDR4 with METTL1 regarding the ciliopathy-like phenotypes in fish and cells. We hence knocked down METTL1 in fibroblasts. A pool of siRNAs against METTL1 almost completely depleted all METTL1 RNA, while also producing a very moderate reduction of WDR4 (Fig. S10). Previous experiments showed that KD of WDR4 with such low efficacies does not impact cilia length or number (data not shown). Importantly, efficient METTL1 KD did neither affect ciliation nor cilia length (Fig. S10). These results suggest that WDR4 influences ciliogenesis likely independently of METTL1.

To better understand the underlying mechanism, we analyzed transcriptome changes induced by the WDR4-R84H variant via RNAseq. We identified 3431 differentially expressed genes (DEGs, Fig. S11a), which associated with transcriptional signatures of cell types present specifically in the central nervous system [41] (Fig. S11b). We also detected DEGs known to be connected to cilia [42] as well as DEGs involved in ubiquitin-associated processes [43], yet both sets of genes were not enriched in the group of DEGs (Fig. S11c, d). However, since differential expression of genes with ubiquitin-association seemed more pronounced (Fig. S11e), we focused our attention on ubiquitin-related processes.

Protein modification by covalent addition of single or multiple ubiquitin molecules serves many different purposes in the cell. The most commonly known ubiquitin modification is polyubiquitination of proteins, which targets them for proteasomal degradation. This marks dysfunctional proteins for removal, yet also to recycle the components for new synthesis including supplying the cell with free ubiquitin again [44]. Consistent with the increased appearance of ubiquitin modifications we detected an increased proteasomal activity in WDR4 patient fibroblasts in G0 (Fig. 6a). Addition of the proteasome inhibitor MG-132 completely abolished proteasomal activity in control and in patient cells (Fig. 6a). Western blot analysis of the same cells further demonstrated an increase in protein ubiquitination (Fig. 6b). We then tested, whether inhibition of the proteasome has an influence on cilia in patient cells and observed increased cilium length upon treatment with MG-132 (Fig. 6c). Furthermore, MG-132 treatment of zebrafish embryos was in itself sufficient to at least partially rescue the size of head and eyes in zebrafish (Fig. 6d, e). MG-132 also increased the number of neuronal cells in the neural tube, which were reduced upon Wdr4 KD (Fig. 6f). These data suggest that uncontrolled proteasomal degradation of proteins and subsequent cilium extension deficits due to WDR4 loss-of-function could potentially reduce the number of neuronal cells during embryogenesis and in turn predispose to microcephaly.Fig. 6 WDR4 patient fibroblasts display increased proteasomal activity.

a Proteasome activity is increased in WDR4 patient fibroblasts compared to control cells. Specificity of the assay is shown by complete signal abrogation by treatment with the proteasome inhibitor MG-132 (25 µM). n = 3 experiments. **p < 0.01. One-way ANOVA with Sidak’s post-test. b WDR4 patient cells display higher ubiquitination. Western blot for ubiquitin and GAPDH as loading control. **p = 0.0079. Mann–Whitney test. c Proteasome inhibition rescues cilia length in WDR4 patient fibroblasts. n = 3 experiments with 91–98 cilia in total. ****p < 0.0001. Mann–Whitney test. d Proteasome inhibition partially rescues the microcephaly in zebrafish embryos at 48 hpf. n = 4 experiments with 76–116 embryos in total. *p < 0.05, ****p < 0.0001. One-way ANOVA with Sidak’s post-test. e Treatment with a proteasome inhibitor also rescues the small eye phenotype in 48 hpf zebrafish. n = 4 experiments with 76–116 embryos in total. *p < 0.05, ****p < 0.0001. One-way ANOVA with Sidak’s post-test. f Wdr4 KD reduces the number of neuronal cells (Elavl3 + 4, green) in the neural tube of 24 hpf zebrafish, which can be rescued by proteasomal inhibition. Counterstaining for acetylated tubulin (magenta). n = 3 experiments with 45 embryos in each condition. ****p < 0.0001, Kruskal–Wallis test with Dunn’s post-test. Scale bar: 50 µm. g Simultaneous knockdown of HSPA1A partially rescues ciliation defect in human fibroblasts. n = 3, siCTRL1/siCTRL2 vs. siWDR4/siCTRL2 ***p = 0.0004; siWDR4/siCTRL2 vs. siHSPA1A/siCTRL2 ***p = 0.0007; siWDR4/siCTRL2 vs. siWDR4/siHSPA1A *p = 0.0255. One-way ANOVA with Sidak correction. h Simultaneous knockdown of HSPA1A rescues cilia length in human fibroblasts. n = 3, 96–104 cilia, siCTRL1/siCTRL2 vs. siWDR4/siCTRL2 ****p < 0.0001; siWDR4/siCTRL2 vs. siHSPA1A/siCTRL2 *p = 0.0119; siWDR4/siCTRL2 vs. siWDR4/siHSPA1A ****p < 0.0001; Kruskal–Wallis test with Dunn’s post-test. i The HSPA1A inhibitor pifithrin-µ restores cilium length in WDR4 patient cells. Fibroblasts were treated with 5 µM pifithrin-µ during the whole time of serum starvation. n = 3 with 103 (CTRL) and 112 (pifithrin-µ) cilia in total. ****p < 0.0001. Two-tailed Welch’s test. j qPCR analyses shows reduced expression of ubiquitin precursors in serum-starved patient fibroblasts. *p = 0.0156 (UBA52, UBB, UBC) and 0.0312 (UBA80). Wilcoxon tests. k Ciliation rate can be rescued in patient cells nucleofected with a ubiquitin plasmid. n = 3 experiments. Per experiment at least 100 cells were counted. *p = 0.0246. Paired t-test. l Cilium length in WDR4 patient fibroblasts can be rescued by nucleofection of a ubiquitin plasmid; n = 3 with 105 cilia each, ****p < 0.0001. Welch’s t test. m Cilium length in WDR4 patient fibroblasts can be rescued by inhibition of ubiquitin activating enzymes with 5 nM TAK-243; n = 3 with 96 and 99 cilia. **p < 0.0062. Mann–Whitney test. Red line: median.

In order to find an explanation for the increased proteasomal activity we turned to our mass spectrometry analysis. HSPA1A was the most abundant protein in patient cells and 5-fold upregulated compared to control cells (Table S1). Interestingly, high levels of HSPA1A have been reported to stimulate the proteasomal activity [45]. Therefore, we performed double KD of WDR4 and HSPA1A in human fibroblasts. KD efficiency could be verified and neither smartpool altered the expression of the other transcript (Fig. S12a, b). While depletion of HSPA1A alone affected neither ciliation nor cilium length, we observed a complete rescue of cilium formation and cilium length in WDR4 siRNA cells co-transfected with siRNAs against HSPA1A (Fig. 6g, h). In addition, pharmacological inhibition of HSPA1A using the small molecule pifithrin-µ [46] similarly restored cilium length in WDR4 patient cells (Fig. 6i). We further hypothesized that the production of ubiquitin precursors should be upregulated to render the increased proteasomal activity possible. Contrary to that we detected a reduction of all four ubiquitination precursors using qPCR on control and WDR4 patient fibroblasts after serum starvation (Fig. 6j). Yet, nucleofection of patient cells with a construct encoding for free ubiquitin, which cannot be conjugated [47] increased the percentage of cells extending cilia and these cilia were significantly longer (Fig. 6k, i). Additionally, we suppressed ubiquitination in WDR4 patient fibroblasts by treatment with TAK-243, an inhibitor of ubiquitin activating enzymes. This was sufficient to trigger a subtle, yet significant rescue of cilia length (Fig. 6m). Together, these data suggest that loss of WDR4 may lead to a shortage of free ubiquitin, which the cell attempts to counter-regulate by increased proteasomal activity.

Discussion

WDR4 has originally been identified as a co-factor of METTL1 in the modification of tRNAs at position 46 [40]. We initially hypothesized the WDR4-METTL1 complex to be essential for protein synthesis as its dysfunction in yeast led to the degradation of tRNAs [48]. One explanation for impaired brain and head growth in WDR4 patients and KD zebrafish would thus be a failure in overall or selected protein synthesis. This idea was further bolstered by a report about another methyltransferase (NSUN2), which is implicated in microcephaly and associated with decreased translation of mRNAs into protein [49]. Furthermore, deregulation of protein synthesis has been postulated to be causative for other neurodevelopmental disorders, e.g., Rett syndrome [50]. WDR4 patient fibroblasts, however, did not display any impairment in protein synthesis, instead de novo translation was elevated. Moreover, the increase in de novo protein synthesis did not yield a higher protein content, which may be the consequence of the observed higher proteasomal activity. Our subsequent hypothesis that loss of WDR4 may lead to errors in protein sequences and would hence activate the proteasome to clear out the faulty proteins could neither be confirmed. WDR4 patient cells possess a normal translation and folding fidelity. These observations suggest that the cilia defects upon WDR4 loss-of-function are not necessarily due to a general lack in proteins or accumulation of dysfunctional proteins. Intriguingly, though, inhibition of the proteasome by MG-132 rescued cilium length as well as partially the microcephaly and neurogenesis defect in zebrafish. The proteasome can impact on cilia in multiple ways, such as prevention as well as promotion of ciliogenesis, control of axoneme extension and cilium resorption (reviewed in [51]). Additionally, MG-132 has been used to prevent degradation of specific factors that affect ciliogenesis such as trichoplein [52] or applied to globally affect proteostasis to rescue cilium length in contexts such as autophagy-defective cells [53]. Interestingly, MG-132 enriches mostly newly synthesized proteins, which characteristically are short-lived such as signaling and membrane proteins on which the cilium depends greatly [54]. Hence, the cilium defect in WDR4 cells could possibly be caused by a shortage of certain cilia components. We do not have experimental evidence to confirm or reject this hypothesis, but similar observations have been made for cilium defects resulting from decreased mTOR signaling, which was suggested to provide a specific subset of proteins required for ciliogenesis [55].

Our study is fortified by data obtained in mice, where WDR4 has previously been connected to very specific ubiquitination events. Here, WDR4 was instrumental to mark Pml [27] and Arhgap17 [32] for ubiquitin-mediated degradation, the latter being necessary for proper development of the murine cerebellum [32]. Using patient-derived cells, we identify a much further reaching impact of WDR4 on protein synthesis, ubiquitination and hence protein homeostasis. We propose that heatshock proteins as well as ubiquitin pools potentially drive this process. High levels of the heatshock protein HSPA1A, which is upregulated in WDR4 patient cells, provoke increased proteasomal activity [45]. Protein degradation by the proteasome, on the other hand, has been shown to generate free ubiquitin [56]. These data together with our observation that ubiquitination is elevated in WDR4 patient cells renders the question whether WDR4 loss-of-function cells lack sufficient amounts of free ubiquitin. Reduced expression of ubiquitin precursors also points in this direction. Sufficiently high ubiquitin levels are important during development. Lack of UBC results in embryonic lethality at midgestation [57] and deletion of UBB impairs neural differentiation [58]. Moreover, free ubiquitin levels appear to be particularly high in the brain [59]. Indeed, supplementation with free ubiquitin did rescue the cilium defect in cells (Fig. 6h, i).

Cilia are at the heart of development and have previously also been connected to PM [60]. WDR4 dysfunction does not only impair cilium formation in many cells and organs, but impairs also their functionality as we detected dampened activity of different signaling pathways that rely on the cilium. The Hh pathway may be the most important in this context, since it is crucial for brain development. Yet, previous studies failed to detect changes in Hh pathway activity after loss of WDR4 [32]. This may be explained by the lack of agonist stimulation as this is necessary to induce measurable Hh signaling changes in cell culture. The symptoms observed in patients carrying WDR4 variants, however, most likely are due to an interplay of more than one mechanism. The above-mentioned impact of WDR4 on Arhgap17 degradation is an important example of a contributing factor [32] and it would be interesting to assess cilia in the presence of accumulated Arhgap17. Additionally, METTL1 as well as WDR4 propagates the expansion of neuronal cells and the differentiation of neuronal cells from iPSC cells [20] suggesting an additional impact on normal brain development by the canonical METTL1-WDR4 complex. There may also be a centrosomal contribution as in other PM cases. This, however, will be subject of subsequent studies.

Taken together, we here report that WDR4 is indispensable for proper cilium generation. This provides a new axis along which WDR4 contributes to faithful brain and skull development and helps to better understand how WDR4 variants contribute to PM.

Materials and methods

Materials

The following, already established cell lines were obtained from the University of Sussex’s cell line collection at the Genome Damage and Stability Centre: h-Tert immortalized 1BR3 forearm fibroblasts, h-Tert immortalized WDR4 patient-derived fibroblasts, non-immortalized WDR4 patient-derived fibroblasts and lymphoblastoid cell lines from WDR4 patients (c.251G > A; p.R84H). HEK 293T cells were from the American Type Culture Collection. MG-132 and TAK-243 were purchased from Selleckchem, pifithrin-µ from Sigma-Aldrich and ciliobrevin D from MedChemExpress. Antibody sources are indicated in the respective methods sections.

Zebrafish maintenance and manipulation

Zebrafish were maintained in a ZebTec rack system with automated water recycling and adjustment of water parameters based on automated pH, conductivity and temperature measurements (Tecniplast). Fish were maintained under a 14 h light and 10 h dark cycle and fed three to four times a day with living artemia in the morning (Aqua Schwarz) and pelleted food at least twice thereafter (Zebrafeed, Sparos). Used lines in here were AB and EK wild-type lines as well as a reporter line expressing GFP in the glomerulus and proximal pronephros. WDR4 loss-of-function was achieved by microinjection of antisense morpholino oligonucleotides (MO) into the yolk of 1–2 cell stage embryos for ubiquitous targeting. MOs were designed and synthesized based on submitted sequences by Gene Tools Inc (Oregon, USA) and had the following sequences: Wdr4 splMO: 5’-CCCACCTGCAAGAGGACAATCATAC, Wdr4 splCTRL: 5’-CACACCTCCAACAGGAAAATAATAC, Wdr4 ATG MO: 5’-CACTACAGCCATGTTGAGTGTTGTT. In addition, a standard control MO provided by GeneTools was used. For reconstitution experiments, capped RNA encoding different variants of zebrafish Wdr4 (which resembled human wild-type and patient variants) was co-injected with MO. In all experiments, uninjected wild-type controls as well as embryos injected with a non-targeting control MO (Wdr4 splCTRL or the standard control MO provided by GeneTools) of the same clutch of eggs were included to assess clutch quality and to identify injection-based off-target effects. Eggs of the same clutch were distributed randomly into the different treatment groups. Embryos were raised at 28.5 °C and incubated to the desired stages. No blinding occurred as the phenotypes were too obvious to ensure blinding during analyses. Splice blocking was verified by RT-PCR as follows: Total RNA was isolated from 24 hpf embryos using the Quick-RNA Miniprep kit (Zymo Research), equal amounts were transcribed into cDNA using Superscript III reverse transcriptase (Thermo Fisher) and oligodTTPs (IDT). PCR was performed using Tag polymerase (NEB) and the following primers: Forward: 5’- GGA GGG ATC AGA AGA GAA GGA, Reverse: 5’- GCC TGA ATG TTG TAG GGT GAA. In case the MO inhibits splicing, exon 4 should be excised and a smaller PCR band occurs. Treatment with 2.5 µM ciliobrevin D was from tailbud stage on.

Analysis of cartilage formation

Zebrafish embryos were fixed at 4 days post fertilization (dpf) using 4% buffered paraformaldehyde (PFA). After fixation, embryos were shortly rinsed with PBS and transferred first to 50% EtOH in PBS for 10 min and then to staining solution (0.02% Alcian blue, 70% EtOH, 50 mM MgCl2). After overnight incubation at room temperature zebrafish were once rinsed with H2O and then bleached for 20 min (bleaching solution: 1.5% H2O2 in 1% KOH) and gradually cleared (30 min 20% Glycerol/0.25% KOH, 2 h 50% Glycerol/0.1% KOH). Stained embryos were stored at 4 °C in 50% Glycerol/0.1% KOH.

In situ hybridization

Whole mount in situ hybridization (ISH) was performed according to standard protocols using DIG-labeled antisense probes, an alkaline phosphatase-coupled anti-DIG antibody (Roche) and NBT/BCIP (Roche) for color development. A probe was generated to detect wdr4 by NotI linearization of the plasmid described below and in vitro transcription using SP6 RNA polymerase (NEB) along with Roche’s DIG labeling system. All other probes used have been described previously [61].

Cloning

For generation of an antisense probe against zebrafish wdr4 a 986 bp fragment of NM_001126435 was amplified using the Expand long template PCR System (Roche) from cDNA of 24 hours post fertilization (hpf) zebrafish embryos. The PCR was cleaned up with the Qiagen nucleotide removal system and cloned into pCRII by TOPO TA cloning (Life technologies). For generation of capped RNA zebrafish Wdr4 (NM_001126435.2) was cloned via BamHI and EcoRI into pCS2-GFP resulting in a C-terminal fusion of GFP to Wdr4. Site-directed mutagenesis was used to generate the patient-resembling pCS2-GFP-Wdr4-R84H. Capped RNA was synthesized using Ambion’s SP6 mMessage mMachine Kit after linearization of the template plasmid using NotI. A plasmid containing human WDR4 was obtained from Origene. Site-directed mutagenesis was used to insert the R84H mutation.

Cell culture and transfection

hTert-immortalized fibroblasts were cultured in MEM alpha containing essential amino acids, penicillin/streptomycin and 10% fetal calf serum (FCS) (all Thermo Fisher) at 37 °C and 5% CO2. Primary cells were cultured just the same except for 15% FCS in the medium. Lymphoblastoid cells were cultured in RPMI 1640 medium (Thermo Fisher) supplemented with 15% FCS, penicillin, and streptomycin. siRNA-mediated transfections were carried out using pools of siRNAs (either Dharmacon’s SmartPools or siRNA pools from siTools) and Lipofectamine RNAiMAX (Thermo Fisher) according to the manufacturer’s protocol. Nucleofections were carried out with the AMAXA nucleofector IIb and the AMAXA Nucleofector Kit R (both Lonza). Cilium formation was induced by serum deprivation (0.1% FCS) for 72 h. Cells were regularly monitored for mycoplasma contamination and authenticated using highly polymorphic short tandem repeat loci sequencing. A plasmid carrying GFP-tagged zebrafish Wdr4 was used for rescue experiments.

qPCR

RNA was isolated using Zymo’s Quick-RNA MiniPrep system, which includes genomic DNA removal and an on-column DNase digest. cDNA was transcribed with equal amounts of total RNA using either Thermo Fisher’s Super Script II or NEB’s Protoscript II kit, in both cases using oligodTTP for priming and following the manufacturer’s protocol. qPCR was performed with Roche’s Universal Probe System, Luna Universal Probe qPCR Master Mix (NEB) and a LightCycler 480 (Roche) or Quantstudio 5 (Thermo Fisher Scientific). qPCR assays were designed with the help of the Universal Probe Assay Design Center (Roche). Primers and corresponding probes are given in Supplementary Table 1.

Protein synthesis assay

Newly synthesized protein amounts were measured with the help of the Click-iT® Plus OPP Protein Synthesis Assay Kits (cat. no. C10456, Molecular Probes) according to the manufacturer’s instructions.

Proteasome activity assay

Proteasome activity in human wild-type and patient-derived fibroblasts was analyzed with the Proteasome 20S Activity Assay Kit (cat. no. MAK172, Sigma).

Translation fidelity

Translation fidelity was measured essentially as described before ([62] and references therein). 106 cells were nucleofected using the Amaxa system (Lonza), using 5 µg firefly control/mutant reporter plasmid and 0.1 µg renilla luciferase plasmid. Cells were plated in a 96-well plate (5 × 104 cells/well in 75 µl) and allowed to grow for 24 h. Luciferase activity was measured using Dual-Glo assay (Promega). The ratio of firefly to renilla luciferase activity was calculated and used as an indicator of translational fidelity.

Measurement of protein content

1BR3 fibroblasts and primary WDR4-derived fibroblasts were cultured in a 10 cm dish. When cells reached confluency, cells were trypsinized, counted using a Countess II (Thermo Fisher) and 350,000 cells each pelleted in an Eppendorf tube. Cells were washed once with pre-warmed PBS, pelleted again and all liquid was carefully aspirated. Cells were lysed then in SDS lysis buffer (see also Western blotting) and protein content relative to a bovine serum albumin calibration curve was measured with the help of a bicinchoninic acid (BCA) kit (Sigma-Aldrich).

Sample preparation and mass spectrometry (MS) analysis

Cell pellets for whole proteome analysis were lysed using 7 M Urea, 2 M Thiourea, 30 mM Tris-HCl, pH 8.5. After 1:6 dilution, proteins were reduced, alkylated and subsequently digested overnight using a trypsin to protein ratio of 1:50. The resulting peptides were lyophilized, reconstituted in 15 µl 5% TFA and analyzed via HPLC-MS. Both, chromatography and MS analysis was performed as described previously [63], using an LTQ Orbitrap Velos Pro system (Thermo Fisher Scientific). Data have been deposited at MassIVE (Id: MSV000092847).

MS data analysis and statistics

For peptide distribution analysis, database search was performed using MaxQuant version 1.6.3.4 (www.maxquant.org) [64]. Employing the build-in Andromeda search engine [65], MS/MS spectra were correlated with the UniProt human reference proteome set (www.uniprot.org) for peptide identification. Carbamidomethylated cysteine was considered as a fixed modification along with oxidation (M), and acetylated protein N-termini as variable modifications. False Discovery rates were set on both, peptide and protein level, to 0.01.

Data analysis was repeated using PEAKS Studio (BSI Inc., Version 8.5) employing a database search with similar parameters and conducting a subsequent SPIDER search to allow for AA exchange and PTM discovery. For normalization, peptide counts were divided by the total peptide count in the respective sample.

RNA sequencing and bioinformatical analysis

Total RNA was isolated with the Quick-RNA MicroPrep system (Zymo) after cilia induction through serum starvation. Sequencing of RNA samples was performed using Illumina’s next-generation sequencing methodology [66]. In detail, total RNA was quantified and quality checked using Tapestation 4200 in combination with RNA ScreenTape (both Agilent Technologies). RIN was >9.5 for all samples (average 9.9). Libraries were prepared from 200 ng of input material (total RNA) using NEBNext Ultra II Directional RNA Library Preparation Kit in combination with NEBNext Poly(A) mRNA Magnetic Isolation Module and NEBNext Multiplex Oligos for Illumina (Unique Dual Index UMI Adaptors RNA) following the manufacturer’s instructions (New England Biolabs). Quantification and quality checked of libraries was done using an Tapestation 4200 and a D5000 assay (Agilent Technologies). Libraries were pooled and sequenced on a NovaSeq 6000. System runs in SP v1.5/101 cycle/single-end/standard loading workflow mode. Using 19 cycles in index 1 read, the index as well as the UMI sequences were obtained. Sequence information was converted to FASTQ format using bcl2fastq v2.20.0.422. Sequencing approach resulted in 58–75 m total reads per sample (av. 66 m).

Mapping of reads to the reference was done using TopHat v2.1.0 [67] using parameters --no-convert-bam --no-coverage-search -x 1 -g 1. As reference the human Ensemble genome version GRCh38 was used taking the gene annotation Release 92 into account [68]. Depending on the sample 77–83% of the reads could be mapped uniquely to the reference. PCR duplicons occurred in library preparation were removed after mapping of reads based on UMIs (unique molecular identifier) using umi_tools v1.1.1 [69]. Around 25–60% of reads per sample passed de-duplication process. Remaining reads per gene were counted for each sample using featureCounts v2.0.3 [70] using parameter -s2. Around 90% of reads per sample could be assigned for counting. Gene counts were further processed using the programming language R. In order to find differentially expressed genes (DEG) the statistical packageDESeq2 [71] was applied. p values were adjusted for multiple testing according to Benjamini–Hochberg correction (FDR). Genes with adjusted p value < 0.05 were regarded to be differentially expressed. For further analysis, changes of 20% in expression were regarded as relevant. Metascape (v3.5.20230501) [41], Syscilia [42], and Biogrid [43] platforms were used to analyze data and identify genes associated with cilia and ubiquitin-related processes, respectively.

Western blotting and PDGF AA assay

Fibroblasts were seeded at a density of 60,000 cells per well in a 12-well plate. On the next day, full medium was removed and cells were exposed to starvation medium for the following 72 h. Cells were then washed once in pre-warmed PBS and lysed in 60 µl SDS lysis buffer (2% SDS and 50 mM Tris pH 6.8) containing protease and phosphatase inhibitors (Roche, Germany). In case of PDGF AA assays, cells were treated with DMSO as control or PDGF AA (50 ng/ml, BioLegend) during 3 min before lysis. Lysates were cleared by nuclease treatment (Pierce universal nuclease, Thermo Fisher) and centrifugation. Protein content was assessed by BCA (Sigma-Aldrich) and equal amounts of protein were separated on Bolt Bis-Tris 4–12% gradient gels using MES running buffer (both Thermo Fisher) before blotting onto nitrocellulose membranes (Bio-Rad). Blots were blocked in 3% milk in Tris buffered saline containing 0.1% NP-40 (Applichem) and subsequently incubated in primary antibody over night at 4 °C. After three washes with Tris buffered saline containing 0.1% NP-40, secondary antibodies were applied for 2 h at room temperature or at 4 °C overnight. Signals were detected using a Li-COR Odyssey SA system and the preinstalled Li-COR software or alternatively with an iBright FL1500 and iBright Analysis Software (version 5.2.0, both Thermo Fisher). Chemiluminescence was used to detect signals in Westerns shown in Fig. S3a. Primary antibodies used were: mouse anti-GAPDH (Acris antibodies, clone 6C5, 1:1000 or Millipore cat. no. MAB6374, 1:500), mouse anti-phospho-AKT S473 (Rockland, cat. no. 200-301-268, 1:1000), rabbit anti-AKT (Cell Signaling, clone C67E7, cat. no. 4691, 1:1000), mouse anti-ubiquitinylated proteins (Merck Millipore, clone FK2, cat. No. 04-263, 1:1000), rabbit anti-WDR4 (Santa Cruz, cat. no. sc-83348, 1:500), rabbit anti cyclin A2 (Proteintech, cat. no. 18202-1-AP, 1:500), and mouse-anti-γTub (Sigma,1:1000) Secondary antibodies for Li-COR were IR-dye labeled and used at 1:10,000 dilution (Li-COR) or for iBright FL1500 Alexa Fluor plus 800 or 647 labeled (Thermo Fisher). Uncropped images of membranes are shown in Fig. S13.

Immunofluorescence and cilium analysis

Fibroblasts grown on cover slips were fixed either with 4% PFA at room temperature for 10 min or using ice-cold methanol for a duration of 4 min. Before antibody incubation, cells underwent a 10 min permeabilization step using 0.1% Triton X-100 in PBS and 1 h of blocking in PBS containing 10% (FCS) or 5% bovine serum albumin Fraction V. Primary antibodies were diluted in blocking buffer and applied to cells over night at 4 °C. After three washes, secondary antibodies diluted in blocking buffer were applied for 2 h at room temperature. After three more washes with PBS, coverslips were mounted onto glass slides using Vectashield mounting medium containing Dapi (Vector labs) and edges sealed with nail polish.

Zebrafish embryos were fixed over night at 4 °C in 4% PFA and stored in absolute methanol at −20 °C. Cilia were stained according to standard protocols [72]. In order to visualize cilia in the Kupffer’s vesicle the tail portion of the embryos was separated from the rest of the embryo and the majority of yolk was gently removed using forceps. This deyolked portion of the embryo was then flat-mounted in a small drop of Vectashield containing Dapi between two cover slips sealed with vacuum grease.

Antibodies used were: rabbit anti-CEP97 (Sigma-Aldrich, catalog no. HPA002980, 1:200), rabbit anti-CP110 (Sigma-Aldrich, cat. no. HPA039402, 1:100), rabbit anti-Elav3 + 4 (GeneTex, cat. no. GTX128365, 1:1000), rabbit-anti-pericentrin (Sigma-Aldrich, cat. no. HPA019887, 1:1000), rabbit anti-PKCζ (Santa Cruz Biotechnology, cat. no. sc-216, 1:500), rabbit anti cyclin A2 (Proteintech, cat. no. 18202-1-AP, 1:500), rabbit anti-centrin 1 (Proteintech, cat.no. 12794-1-AP, 1:300), mouse anti-acetylated tubulin (Sigma-Aldrich, clone 6-11B-1, 1:1000), mouse anti-γTubulin (Sigma-Aldrich, clone GTU-88, 1:5000), rabbit anti-γTubulin (Sigma-Aldrich, T6793, 1:5000), mouse anti-alpha1-tubulin-FITC (Sigma, F2168, 1:500). Unlabeled primary antibodies were detected with Alexa-coupled secondary antibodies (Thermo Fisher). CoraLite594-Phalloidin (Proteintech) was used to detect stress fibers.

Imaging

Live zebrafish embryos and such processed by whole mount ISH were imaged on a Leica M125 with an IC80HD or a MC190HD camera (Leica, Wetzlar, Germany). For the analysis of smaller heads, eyes and shorter bodies, injected embryos were compared to average non-injected and control-injected controls. KV cilia in flat-mounted tails of zebrafish embryos as well as primary cilia in human fibroblasts were imaged using a Leica TCS SP5II confocal microscope or Leica Stellaris 5. Here, z-stacks were captured every 0.3–4 µm and processed into a 3D stack with the Leica acquisition software LAS AF and LAS X, respectively. All image processing was done using Adobe Photoshop CC2020 (Adobe Systems Software Ireland Limited, Dublin, Ireland). Adjustments (brightness, contrast or color balance) were applied equally to the whole image and equally to same treatment groups (i.e., CTRL MO vs. Wdr4 MO). Final figures were assembled in Adobe Illustrator CC2021.

Cilium and neuronal cell measurements

To measure cilia length, cilia were traced in z-stacks and the length measured in FIJI (NIH). Neural tubes of 24 hpf zebrafish embryos were imaged on a M205FCA stereomicroscope equipped with a DFC9000 GT sCMOS camera and the LAS X software (all Leica). Elav3 + 4-positive cell bodies were counted along a 300 µm stretch of the neural tube.

Statistical analysis

Data were analyzed using Prism 8, 9, and 10 respectively (GraphPad). Data were first assessed for normal distribution with a Shapiro–Wilks test, before the respective two-tailed parametric or non-parametric test was applied. Whenever more than two conditions were compared, an ANOVA-based test was applied. Transformation was applied to make data samples obtained as percentage amenable for analysis by ANOVA [73].

Supplementary information

Supplemental Material

Supplementary information

The online version contains supplementary material available at 10.1038/s41419-024-07042-5.

Acknowledgements

We thank Lisa Woodbine and the Cell Culture Collection at the University of Sussex for providing the immortalized 1BR3 control fibroblast line and the WDR4 patient cell lines, Lennert Janssen for the ubiquitin plasmid and the animal care takers for excellent fish care. The Core Facility DNA sequencing (Ivonne Goerlich) of the FLI is gratefully acknowledged for their technological support in library preparation and sequencing. The Philipp lab has been supported by grants from Deutsche Forschungsgemeinschaft (DFG, German Research Foundation, grant numbers PH144/4-1, PH144/4-3, PH144/6-1, and 467868420), the medical school at University of Tuebingen (Fortüne grant no. F1534033, IZKF grant no. 2020-2-09), the Fundació La Marató de TV3 (grant number 202019-32) and the Million Dollar Bike Ride Program (MDBR-22-136-Cohen). LM and TCT were fellows of the International Graduate School in Molecular Medicine at Ulm University (funded by DFG). HW and parts of proteomics method development were funded by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation)—SFB 1074.

Author contributions

MDB, TS, LDM, TCT, JG, SAS, MDP, LATV, CD, and MP performed and analyzed experiments; HW and SW performed and analyzed mass spectrometry experiments; MA and SI performed and analyzed translation fidelity assay; MG performed and analyzed RNA-seq; MDB and MP performed functional analysis of transcriptome data; JAC collected patient cells; MDB and MP wrote the manuscript; PAJ and MP conceived the study.

Funding

Open Access funding enabled and organized by Projekt DEAL.

Data availability

The datasets generated during the current study are available from the corresponding author upon reasonable request. Datasets of transcriptome analysis have been deposited in NCBI’s Gene Expression Omnibus [74] and are accessible through GEO Series accession number GSE242969. Datasets of proteome analysis of the study have been deposited in the MassIVE repository with the identifier MSV000092847 (https://massive.ucsd.edu/ProteoSAFe/private-dataset.jsp?task=9c0bb17c77249dfb4cb2bfc5afd2c7a).

Competing interests

The authors declare no competing interests.

Ethics approval and consent to participate

All methods were performed in accordance with relevant guidelines and regulations. The study was approved by the University of Sussex School of Life Sciences Research Governance Committee. Ethical consent has been obtained from the individuals at the time of collection and also for the purpose of experimentation. Both, zebrafish husbandry as well as all experiments described in here have been approved by local authorities (Veterinary Care Units at Ulm University (registry no. 0140) and University of Tübingen and the animal welfare commissioner of the regional board for scientific animal experiments in Tübingen, Germany (35/9185.46/Uni Tü)). Experiments were performed according to the European Union Directive 2010/63/EU.

Edited by Anastasis Stephanou

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. Harris SR Measuring head circumference: update on infant microcephaly Can Fam Physician 2015 61 680 4 26505062
Harris SR. Measuring head circumference: update on infant microcephaly. Can Fam Physician. 2015;61:680–4.26505062
2. Alcantara D O’Driscoll M Congenital microcephaly Am J Med Genet C Semin Med Genet 2014 166C 124 39 10.1002/ajmg.c.31397 24816482
Alcantara D, O’Driscoll M. Congenital microcephaly. Am J Med Genet C Semin Med Genet. 2014;166C:124–39.24816482 10.1002/ajmg.c.31397
3. Bicknell LS Walker S Klingseisen A Stiff T Leitch A Kerzendorfer C Mutations in ORC1, encoding the largest subunit of the origin recognition complex, cause microcephalic primordial dwarfism resembling Meier-Gorlin syndrome Nat Genet 2011 43 350 5 10.1038/ng.776 21358633
Bicknell LS, Walker S, Klingseisen A, Stiff T, Leitch A, Kerzendorfer C, et al. Mutations in ORC1, encoding the largest subunit of the origin recognition complex, cause microcephalic primordial dwarfism resembling Meier-Gorlin syndrome. Nat Genet. 2011;43:350–5.21358633 10.1038/ng.776
4. Griffith E Walker S Martin CA Vagnarelli P Stiff T Vernay B Mutations in pericentrin cause Seckel syndrome with defective ATR-dependent DNA damage signaling Nat Genet 2008 40 232 6 10.1038/ng.2007.80 18157127
Griffith E, Walker S, Martin CA, Vagnarelli P, Stiff T, Vernay B, et al. Mutations in pericentrin cause Seckel syndrome with defective ATR-dependent DNA damage signaling. Nat Genet. 2008;40:232–6.18157127 10.1038/ng.2007.80
5. Harley ME Murina O Leitch A Higgs MR Bicknell LS Yigit G TRAIP promotes DNA damage response during genome replication and is mutated in primordial dwarfism Nat Genet 2016 48 36 43 10.1038/ng.3451 26595769
Harley ME, Murina O, Leitch A, Higgs MR, Bicknell LS, Yigit G, et al. TRAIP promotes DNA damage response during genome replication and is mutated in primordial dwarfism. Nat Genet. 2016;48:36–43.26595769 10.1038/ng.3451
6. Murray JE van der Burg M IJspeert H Carroll P Wu Q Ochi T Mutations in the NHEJ component XRCC4 cause primordial dwarfism Am J Hum Genet 2015 96 412 24 10.1016/j.ajhg.2015.01.013 25728776
Murray JE, van der Burg M, IJspeert H, Carroll P, Wu Q, Ochi T, et al. Mutations in the NHEJ component XRCC4 cause primordial dwarfism. Am J Hum Genet. 2015;96:412–24.25728776 10.1016/j.ajhg.2015.01.013
7. O’Driscoll M Ruiz-Perez VL Woods CG Jeggo PA Goodship JA A splicing mutation affecting expression of ataxia-telangiectasia and Rad3-related protein (ATR) results in Seckel syndrome Nat Genet 2003 33 497 501 10.1038/ng1129 12640452
O’Driscoll M, Ruiz-Perez VL, Woods CG, Jeggo PA, Goodship JA. A splicing mutation affecting expression of ataxia-telangiectasia and Rad3-related protein (ATR) results in Seckel syndrome. Nat Genet. 2003;33:497–501.12640452 10.1038/ng1129
8. Ogi T Walker S Stiff T Hobson E Limsirichaikul S Carpenter G Identification of the first ATRIP-deficient patient and novel mutations in ATR define a clinical spectrum for ATR-ATRIP Seckel Syndrome PLoS Genet 2012 8 e1002945 10.1371/journal.pgen.1002945 23144622
Ogi T, Walker S, Stiff T, Hobson E, Limsirichaikul S, Carpenter G, et al. Identification of the first ATRIP-deficient patient and novel mutations in ATR define a clinical spectrum for ATR-ATRIP Seckel Syndrome. PLoS Genet. 2012;8:e1002945.23144622 10.1371/journal.pgen.1002945
9. Woodbine L Neal JA Sasi NK Shimada M Deem K Coleman H PRKDC mutations in a SCID patient with profound neurological abnormalities J Clin Invest 2013 123 2969 80 10.1172/JCI67349 23722905
Woodbine L, Neal JA, Sasi NK, Shimada M, Deem K, Coleman H, et al. PRKDC mutations in a SCID patient with profound neurological abnormalities. J Clin Invest. 2013;123:2969–80.23722905 10.1172/JCI67349
10. O’Driscoll M Cerosaletti KM Girard PM Dai Y Stumm M Kysela B DNA ligase IV mutations identified in patients exhibiting developmental delay and immunodeficiency Mol Cell 2001 8 1175 85 10.1016/S1097-2765(01)00408-7 11779494
O’Driscoll M, Cerosaletti KM, Girard PM, Dai Y, Stumm M, Kysela B, et al. DNA ligase IV mutations identified in patients exhibiting developmental delay and immunodeficiency. Mol Cell. 2001;8:1175–85.11779494 10.1016/S1097-2765(01)00408-7
11. Stiff T Alagoz M Alcantara D Outwin E Brunner HG Bongers EM Deficiency in origin licensing proteins impairs cilia formation: implications for the aetiology of Meier-Gorlin syndrome PLoS Genet 2013 9 e1003360 10.1371/journal.pgen.1003360 23516378
Stiff T, Alagoz M, Alcantara D, Outwin E, Brunner HG, Bongers EM, et al. Deficiency in origin licensing proteins impairs cilia formation: implications for the aetiology of Meier-Gorlin syndrome. PLoS Genet. 2013;9:e1003360.23516378 10.1371/journal.pgen.1003360
12. Stiff T Casar Tena T O'Driscoll M Jeggo PA Philipp M ATR promotes cilia signalling: links to developmental impacts Hum Mol Genet 2016 25 1574 87 10.1093/hmg/ddw034 26908596
Stiff T, Casar Tena T, O'Driscoll M, Jeggo PA, Philipp M. ATR promotes cilia signalling: links to developmental impacts. Hum Mol Genet. 2016;25:1574–87.26908596 10.1093/hmg/ddw034
13. Phan TP Holland AJ Time is of the essence: the molecular mechanisms of primary microcephaly Genes Dev 2021 35 1551 78 10.1101/gad.348866.121 34862179
Phan TP, Holland AJ. Time is of the essence: the molecular mechanisms of primary microcephaly. Genes Dev. 2021;35:1551–78.34862179 10.1101/gad.348866.121
14. Casar Tena T Burkhalter MD Philipp M Left-right asymmetry in the light of TOR: an update on what we know so far Biol Cell 2015 107 306 18 10.1111/boc.201400094 25943139
Casar Tena T, Burkhalter MD, Philipp M. Left-right asymmetry in the light of TOR: an update on what we know so far. Biol Cell. 2015;107:306–18.25943139 10.1111/boc.201400094
15. Lovera M Luders J The ciliary impact of nonciliary gene mutations Trends Cell Biol 2021 31 876 87 10.1016/j.tcb.2021.06.001 34183231
Lovera M, Luders J. The ciliary impact of nonciliary gene mutations. Trends Cell Biol. 2021;31:876–87.34183231 10.1016/j.tcb.2021.06.001
16. Michaud J Kudoh J Berry A Bonne-Tamir B Lalioti MD Rossier C Isolation and characterization of a human chromosome 21q22.3 gene (WDR4) and its mouse homologue that code for a WD-repeat protein Genomics 2000 68 71 9 10.1006/geno.2000.6258 10950928
Michaud J, Kudoh J, Berry A, Bonne-Tamir B, Lalioti MD, Rossier C, et al. Isolation and characterization of a human chromosome 21q22.3 gene (WDR4) and its mouse homologue that code for a WD-repeat protein. Genomics. 2000;68:71–9.10950928 10.1006/geno.2000.6258
17. Jin X Guan Z Hu N He C Yin P Gong Z Structural insight into how WDR4 promotes the tRNA N7-methylguanosine methyltransferase activity of METTL1 Cell Discov 2023 9 65 10.1038/s41421-023-00562-y 37369656
Jin X, Guan Z, Hu N, He C, Yin P, Gong Z, et al. Structural insight into how WDR4 promotes the tRNA N7-methylguanosine methyltransferase activity of METTL1. Cell Discov. 2023;9:65.37369656 10.1038/s41421-023-00562-y
18. Li J Wang L Hahn Q Nowak RP Viennet T Orellana EA Structural basis of regulated m(7)G tRNA modification by METTL1-WDR4 Nature 2023 613 391 7 10.1038/s41586-022-05566-4 36599985
Li J, Wang L, Hahn Q, Nowak RP, Viennet T, Orellana EA, et al. Structural basis of regulated m(7)G tRNA modification by METTL1-WDR4. Nature. 2023;613:391–7.36599985 10.1038/s41586-022-05566-4
19. Ruiz-Arroyo VM Raj R Babu K Onolbaatar O Roberts PH Nam Y Structures and mechanisms of tRNA methylation by METTL1-WDR4 Nature 2023 613 383 90 10.1038/s41586-022-05565-5 36599982
Ruiz-Arroyo VM, Raj R, Babu K, Onolbaatar O, Roberts PH, Nam Y. Structures and mechanisms of tRNA methylation by METTL1-WDR4. Nature. 2023;613:383–90.36599982 10.1038/s41586-022-05565-5
20. Lin S Liu Q Lelyveld VS Choe J Szostak JW Gregory RI Mettl1/Wdr4-mediated m(7)G tRNA methylome is required for normal mRNA translation and embryonic stem cell self-renewal and differentiation Mol Cell 2018 71 244 55.e5 10.1016/j.molcel.2018.06.001 29983320
Lin S, Liu Q, Lelyveld VS, Choe J, Szostak JW, Gregory RI. Mettl1/Wdr4-mediated m(7)G tRNA methylome is required for normal mRNA translation and embryonic stem cell self-renewal and differentiation. Mol Cell. 2018;71:244–55.e5.29983320 10.1016/j.molcel.2018.06.001
21. Cheng IC Chen BC Shuai HH Chien FC Chen P Hsieh TS Wuho is a new member in maintaining genome stability through its interaction with flap endonuclease 1 PLoS Biol 2016 14 e1002349 10.1371/journal.pbio.1002349 26751069
Cheng IC, Chen BC, Shuai HH, Chien FC, Chen P, Hsieh TS. Wuho is a new member in maintaining genome stability through its interaction with flap endonuclease 1. PLoS Biol. 2016;14:e1002349.26751069 10.1371/journal.pbio.1002349
22. Cheng W Gao A Lin H Zhang W Novel roles of METTL1/WDR4 in tumor via m(7)G methylation Mol Ther Oncolytics 2022 26 27 34 10.1016/j.omto.2022.05.009 35784404
Cheng W, Gao A, Lin H, Zhang W. Novel roles of METTL1/WDR4 in tumor via m(7)G methylation. Mol Ther Oncolytics. 2022;26:27–34.35784404 10.1016/j.omto.2022.05.009
23. Deng L Hua RX Deng C Zhu J Zhang Z Cheng J WDR4 gene polymorphisms and Wilms tumor susceptibility in Chinese children: a five-center case-control study J Cancer 2023 14 1293 300 10.7150/jca.83747 37283791
Deng L, Hua RX, Deng C, Zhu J, Zhang Z, Cheng J, et al. WDR4 gene polymorphisms and Wilms tumor susceptibility in Chinese children: a five-center case-control study. J Cancer. 2023;14:1293–300.37283791 10.7150/jca.83747
24. Wang Z Yu P Zou Y Ma J Han H Wei W METTL1/WDR4-mediated tRNA m(7)G modification and mRNA translation control promote oncogenesis and doxorubicin resistance Oncogene 2023 42 1900 12 10.1038/s41388-023-02695-6 37185458
Wang Z, Yu P, Zou Y, Ma J, Han H, Wei W, et al. METTL1/WDR4-mediated tRNA m(7)G modification and mRNA translation control promote oncogenesis and doxorubicin resistance. Oncogene. 2023;42:1900–12.37185458 10.1038/s41388-023-02695-6
25. Chen B Jiang W Huang Y Zhang J Yu P Wu L N(7)-methylguanosine tRNA modification promotes tumorigenesis and chemoresistance through WNT/beta-catenin pathway in nasopharyngeal carcinoma Oncogene 2022 41 2239 53 10.1038/s41388-022-02250-9 35217794
Chen B, Jiang W, Huang Y, Zhang J, Yu P, Wu L, et al. N(7)-methylguanosine tRNA modification promotes tumorigenesis and chemoresistance through WNT/beta-catenin pathway in nasopharyngeal carcinoma. Oncogene. 2022;41:2239–53.35217794 10.1038/s41388-022-02250-9
26. Han H Yang C Ma J Zhang S Zheng S Ling R N(7)-methylguanosine tRNA modification promotes esophageal squamous cell carcinoma tumorigenesis via the RPTOR/ULK1/autophagy axis Nat Commun 2022 13 1478 10.1038/s41467-022-29125-7 35304469
Han H, Yang C, Ma J, Zhang S, Zheng S, Ling R, et al. N(7)-methylguanosine tRNA modification promotes esophageal squamous cell carcinoma tumorigenesis via the RPTOR/ULK1/autophagy axis. Nat Commun. 2022;13:1478.35304469 10.1038/s41467-022-29125-7
27. Wang YT Chen J Chang CW Jen J Huang TY Chen CM Ubiquitination of tumor suppressor PML regulates prometastatic and immunosuppressive tumor microenvironment J Clin Invest 2017 127 2982 97 10.1172/JCI89957 28691927
Wang YT, Chen J, Chang CW, Jen J, Huang TY, Chen CM, et al. Ubiquitination of tumor suppressor PML regulates prometastatic and immunosuppressive tumor microenvironment. J Clin Invest. 2017;127:2982–97.28691927 10.1172/JCI89957
28. Braun DA Shril S Sinha A Schneider R Tan W Ashraf S Mutations in WDR4 as a new cause of Galloway-Mowat syndrome Am J Med Genet A 2018 176 2460 5 10.1002/ajmg.a.40489 30079490
Braun DA, Shril S, Sinha A, Schneider R, Tan W, Ashraf S, et al. Mutations in WDR4 as a new cause of Galloway-Mowat syndrome. Am J Med Genet A. 2018;176:2460–5.30079490 10.1002/ajmg.a.40489
29. Chen X Gao Y Yang L Wu B Dong X Liu B Speech and language delay in a patient with WDR4 mutations Eur J Med Genet 2018 61 468 72 10.1016/j.ejmg.2018.03.007 29597095
Chen X, Gao Y, Yang L, Wu B, Dong X, Liu B, et al. Speech and language delay in a patient with WDR4 mutations. Eur J Med Genet. 2018;61:468–72.29597095 10.1016/j.ejmg.2018.03.007
30. Shaheen R Abdel-Salam GM Guy MP Alomar R Abdel-Hamid MS Afifi HH Mutation in WDR4 impairs tRNA m(7)G46 methylation and causes a distinct form of microcephalic primordial dwarfism Genome Biol 2015 16 210 10.1186/s13059-015-0779-x 26416026
Shaheen R, Abdel-Salam GM, Guy MP, Alomar R, Abdel-Hamid MS, Afifi HH, et al. Mutation in WDR4 impairs tRNA m(7)G46 methylation and causes a distinct form of microcephalic primordial dwarfism. Genome Biol. 2015;16:210.26416026 10.1186/s13059-015-0779-x
31. Trimouille A Lasseaux E Barat P Deiller C Drunat S Rooryck C Further delineation of the phenotype caused by biallelic variants in the WDR4 gene Clin Genet 2018 93 374 7 10.1111/cge.13074 28617965
Trimouille A, Lasseaux E, Barat P, Deiller C, Drunat S, Rooryck C, et al. Further delineation of the phenotype caused by biallelic variants in the WDR4 gene. Clin Genet. 2018;93:374–7.28617965 10.1111/cge.13074
32. Wu PR Chiang SY Midence R Kao WC Lai CL Cheng IC Wdr4 promotes cerebellar development and locomotion through Arhgap17-mediated Rac1 activation Cell Death Dis 2023 14 52 10.1038/s41419-022-05442-z 36681682
Wu PR, Chiang SY, Midence R, Kao WC, Lai CL, Cheng IC, et al. Wdr4 promotes cerebellar development and locomotion through Arhgap17-mediated Rac1 activation. Cell Death Dis. 2023;14:52.36681682 10.1038/s41419-022-05442-z
33. Colantonio JR Vermot J Wu D Langenbacher AD Fraser S Chen JN The dynein regulatory complex is required for ciliary motility and otolith biogenesis in the inner ear Nature 2009 457 205 9 10.1038/nature07520 19043402
Colantonio JR, Vermot J, Wu D, Langenbacher AD, Fraser S, Chen JN, et al. The dynein regulatory complex is required for ciliary motility and otolith biogenesis in the inner ear. Nature. 2009;457:205–9.19043402 10.1038/nature07520
34. Huang N Zhang D Li F Chai P Wang S Teng J M-Phase Phosphoprotein 9 regulates ciliogenesis by modulating CP110-CEP97 complex localization at the mother centriole Nat Commun 2018 9 4511 10.1038/s41467-018-06990-9 30375385
Huang N, Zhang D, Li F, Chai P, Wang S, Teng J, et al. M-Phase Phosphoprotein 9 regulates ciliogenesis by modulating CP110-CEP97 complex localization at the mother centriole. Nat Commun. 2018;9:4511.30375385 10.1038/s41467-018-06990-9
35. Spektor A Tsang WY Khoo D Dynlacht BD Cep97 and CP110 suppress a cilia assembly program Cell 2007 130 678 90 10.1016/j.cell.2007.06.027 17719545
Spektor A, Tsang WY, Khoo D, Dynlacht BD. Cep97 and CP110 suppress a cilia assembly program. Cell. 2007;130:678–90.17719545 10.1016/j.cell.2007.06.027
36. Hoing S Yeh TY Baumann M Martinez NE Habenberger P Kremer L Dynarrestin, a novel inhibitor of cytoplasmic dynein Cell Chem Biol 2018 25 357 69 e6 10.1016/j.chembiol.2017.12.014 29396292
Hoing S, Yeh TY, Baumann M, Martinez NE, Habenberger P, Kremer L, et al. Dynarrestin, a novel inhibitor of cytoplasmic dynein. Cell Chem Biol. 2018;25:357–69 e6.29396292 10.1016/j.chembiol.2017.12.014
37. Macarelli V Leventea E Merkle FT Regulation of the length of neuronal primary cilia and its potential effects on signalling Trends Cell Biol 2023 33 979 90 10.1016/j.tcb.2023.05.005 37302961
Macarelli V, Leventea E, Merkle FT. Regulation of the length of neuronal primary cilia and its potential effects on signalling. Trends Cell Biol. 2023;33:979–90.37302961 10.1016/j.tcb.2023.05.005
38. Huangfu D Liu A Rakeman AS Murcia NS Niswander L Anderson KV Hedgehog signalling in the mouse requires intraflagellar transport proteins Nature 2003 426 83 7 10.1038/nature02061 14603322
Huangfu D, Liu A, Rakeman AS, Murcia NS, Niswander L, Anderson KV. Hedgehog signalling in the mouse requires intraflagellar transport proteins. Nature. 2003;426:83–7.14603322 10.1038/nature02061
39. Schneider L Clement CA Teilmann SC Pazour GJ Hoffmann EK Satir P PDGFRalphaalpha signaling is regulated through the primary cilium in fibroblasts Curr Biol 2005 15 1861 6 10.1016/j.cub.2005.09.012 16243034
Schneider L, Clement CA, Teilmann SC, Pazour GJ, Hoffmann EK, Satir P, et al. PDGFRalphaalpha signaling is regulated through the primary cilium in fibroblasts. Curr Biol. 2005;15:1861–6.16243034 10.1016/j.cub.2005.09.012
40. Alexandrov A Martzen MR Phizicky EM Two proteins that form a complex are required for 7-methylguanosine modification of yeast tRNA RNA 2002 8 1253 66 10.1017/S1355838202024019 12403464
Alexandrov A, Martzen MR, Phizicky EM. Two proteins that form a complex are required for 7-methylguanosine modification of yeast tRNA. RNA. 2002;8:1253–66.12403464 10.1017/S1355838202024019
41. Zhou Y Zhou B Pache L Chang M Khodabakhshi AH Tanaseichuk O Metascape provides a biologist-oriented resource for the analysis of systems-level datasets Nat Commun 2019 10 1523 10.1038/s41467-019-09234-6 30944313
Zhou Y, Zhou B, Pache L, Chang M, Khodabakhshi AH, Tanaseichuk O, et al. Metascape provides a biologist-oriented resource for the analysis of systems-level datasets. Nat Commun. 2019;10:1523.30944313 10.1038/s41467-019-09234-6
42. Vasquez SSV van Dam J Wheway G An updated SYSCILIA gold standard (SCGSv2) of known ciliary genes, revealing the vast progress that has been made in the cilia research field Mol Biol Cell 2021 32 br13 10.1091/mbc.E21-05-0226 34613793
Vasquez SSV, van Dam J, Wheway G. An updated SYSCILIA gold standard (SCGSv2) of known ciliary genes, revealing the vast progress that has been made in the cilia research field. Mol Biol Cell. 2021;32:br13.34613793 10.1091/mbc.E21-05-0226
43. Stark C Breitkreutz BJ Reguly T Boucher L Breitkreutz A Tyers M BioGRID: a general repository for interaction datasets Nucleic Acids Res 2006 34 D535 9 10.1093/nar/gkj109 16381927
Stark C, Breitkreutz BJ, Reguly T, Boucher L, Breitkreutz A, Tyers M. BioGRID: a general repository for interaction datasets. Nucleic Acids Res. 2006;34:D535–9.16381927 10.1093/nar/gkj109
44. Park CW Ryu KY Cellular ubiquitin pool dynamics and homeostasis BMB Rep 2014 47 475 82 10.5483/BMBRep.2014.47.9.128 24924398
Park CW, Ryu KY. Cellular ubiquitin pool dynamics and homeostasis. BMB Rep. 2014;47:475–82.24924398 10.5483/BMBRep.2014.47.9.128
45. Morozov AV Astakhova TM Garbuz DG Krasnov GS Bobkova NV Zatsepina OG Interplay between recombinant Hsp70 and proteasomes: proteasome activity modulation and ubiquitin-independent cleavage of Hsp70 Cell Stress Chaperones 2017 22 687 97 10.1007/s12192-017-0792-y 28447215
Morozov AV, Astakhova TM, Garbuz DG, Krasnov GS, Bobkova NV, Zatsepina OG, et al. Interplay between recombinant Hsp70 and proteasomes: proteasome activity modulation and ubiquitin-independent cleavage of Hsp70. Cell Stress Chaperones. 2017;22:687–97.28447215 10.1007/s12192-017-0792-y
46. Leu JI Pimkina J Frank A Murphy ME George DL A small molecule inhibitor of inducible heat shock protein 70 Mol Cell 2009 36 15 27 10.1016/j.molcel.2009.09.023 19818706
Leu JI, Pimkina J, Frank A, Murphy ME, George DL. A small molecule inhibitor of inducible heat shock protein 70. Mol Cell. 2009;36:15–27.19818706 10.1016/j.molcel.2009.09.023
47. Dantuma NP Groothuis TA Salomons FA Neefjes J A dynamic ubiquitin equilibrium couples proteasomal activity to chromatin remodeling J Cell Biol 2006 173 19 26 10.1083/jcb.200510071 16606690
Dantuma NP, Groothuis TA, Salomons FA, Neefjes J. A dynamic ubiquitin equilibrium couples proteasomal activity to chromatin remodeling. J Cell Biol. 2006;173:19–26.16606690 10.1083/jcb.200510071
48. Alexandrov A Chernyakov I Gu W Hiley SL Hughes TR Grayhack EJ Rapid tRNA decay can result from lack of nonessential modifications Mol Cell 2006 21 87 96 10.1016/j.molcel.2005.10.036 16387656
Alexandrov A, Chernyakov I, Gu W, Hiley SL, Hughes TR, Grayhack EJ, et al. Rapid tRNA decay can result from lack of nonessential modifications. Mol Cell. 2006;21:87–96.16387656 10.1016/j.molcel.2005.10.036
49. Blanco S Dietmann S Flores JV Hussain S Kutter C Humphreys P Aberrant methylation of tRNAs links cellular stress to neuro-developmental disorders EMBO J 2014 33 2020 39 10.15252/embj.201489282 25063673
Blanco S, Dietmann S, Flores JV, Hussain S, Kutter C, Humphreys P, et al. Aberrant methylation of tRNAs links cellular stress to neuro-developmental disorders. EMBO J. 2014;33:2020–39.25063673 10.15252/embj.201489282
50. Rodrigues DC Mufteev M Weatheritt RJ Djuric U Ha KCH Ross PJ Shifts in ribosome engagement impact key gene sets in neurodevelopment and ubiquitination in Rett syndrome Cell Rep 2020 30 4179 96.e11 10.1016/j.celrep.2020.02.107 32209477
Rodrigues DC, Mufteev M, Weatheritt RJ, Djuric U, Ha KCH, Ross PJ, et al. Shifts in ribosome engagement impact key gene sets in neurodevelopment and ubiquitination in Rett syndrome. Cell Rep. 2020;30:4179–96.e11.32209477 10.1016/j.celrep.2020.02.107
51. Senatore E Iannucci R Chiuso F Delle Donne R Rinaldi L Feliciello A Pathophysiology of primary cilia: signaling and proteostasis regulation Front Cell Dev Biol 2022 10 833086 10.3389/fcell.2022.833086 35646931
Senatore E, Iannucci R, Chiuso F, Delle Donne R, Rinaldi L, Feliciello A. Pathophysiology of primary cilia: signaling and proteostasis regulation. Front Cell Dev Biol. 2022;10:833086.35646931 10.3389/fcell.2022.833086
52. Kasahara K Kawakami Y Kiyono T Yonemura S Kawamura Y Era S Ubiquitin-proteasome system controls ciliogenesis at the initial step of axoneme extension Nat Commun 2014 5 5081 10.1038/ncomms6081 25270598
Kasahara K, Kawakami Y, Kiyono T, Yonemura S, Kawamura Y, Era S, et al. Ubiquitin-proteasome system controls ciliogenesis at the initial step of axoneme extension. Nat Commun. 2014;5:5081.25270598 10.1038/ncomms6081
53. Wang S Livingston MJ Su Y Dong Z Reciprocal regulation of cilia and autophagy via the MTOR and proteasome pathways Autophagy 2015 11 607 16 10.1080/15548627.2015.1023983 25906314
Wang S, Livingston MJ, Su Y, Dong Z. Reciprocal regulation of cilia and autophagy via the MTOR and proteasome pathways. Autophagy. 2015;11:607–16.25906314 10.1080/15548627.2015.1023983
54. Yen HC Xu Q Chou DM Zhao Z Elledge SJ Global protein stability profiling in mammalian cells Science 2008 322 918 23 10.1126/science.1160489 18988847
Yen HC, Xu Q, Chou DM, Zhao Z, Elledge SJ. Global protein stability profiling in mammalian cells. Science. 2008;322:918–23.18988847 10.1126/science.1160489
55. Yuan S Li J Diener DR Choma MA Rosenbaum JL Sun Z Target-of-rapamycin complex 1 (Torc1) signaling modulates cilia size and function through protein synthesis regulation Proc Natl Acad Sci USA 2012 109 2021 6 10.1073/pnas.1112834109 22308353
Yuan S, Li J, Diener DR, Choma MA, Rosenbaum JL, Sun Z. Target-of-rapamycin complex 1 (Torc1) signaling modulates cilia size and function through protein synthesis regulation. Proc Natl Acad Sci USA. 2012;109:2021–6.22308353 10.1073/pnas.1112834109
56. Yao T Cohen RE A cryptic protease couples deubiquitination and degradation by the proteasome Nature 2002 419 403 7 10.1038/nature01071 12353037
Yao T, Cohen RE. A cryptic protease couples deubiquitination and degradation by the proteasome. Nature. 2002;419:403–7.12353037 10.1038/nature01071
57. Ryu KY Maehr R Gilchrist CA Long MA Bouley DM Mueller B The mouse polyubiquitin gene UbC is essential for fetal liver development, cell-cycle progression and stress tolerance EMBO J 2007 26 2693 706 10.1038/sj.emboj.7601722 17491588
Ryu KY, Maehr R, Gilchrist CA, Long MA, Bouley DM, Mueller B, et al. The mouse polyubiquitin gene UbC is essential for fetal liver development, cell-cycle progression and stress tolerance. EMBO J. 2007;26:2693–706.17491588 10.1038/sj.emboj.7601722
58. Ryu HW Park CW Ryu KY Disruption of polyubiquitin gene Ubb causes dysregulation of neural stem cell differentiation with premature gliogenesis Sci Rep 2014 4 7026 10.1038/srep07026 25391618
Ryu HW, Park CW, Ryu KY. Disruption of polyubiquitin gene Ubb causes dysregulation of neural stem cell differentiation with premature gliogenesis. Sci Rep. 2014;4:7026.25391618 10.1038/srep07026
59. Kaiser SE Riley BE Shaler TA Trevino RS Becker CH Schulman H Protein standard absolute quantification (PSAQ) method for the measurement of cellular ubiquitin pools Nat Methods 2011 8 691 6 10.1038/nmeth.1649 21743460
Kaiser SE, Riley BE, Shaler TA, Trevino RS, Becker CH, Schulman H, et al. Protein standard absolute quantification (PSAQ) method for the measurement of cellular ubiquitin pools. Nat Methods. 2011;8:691–6.21743460 10.1038/nmeth.1649
60. Thomas S Boutaud L Reilly ML Benmerah A Cilia in hereditary cerebral anomalies Biol Cell 2019 111 217 31 10.1111/boc.201900012 31177551
Thomas S, Boutaud L, Reilly ML, Benmerah A. Cilia in hereditary cerebral anomalies. Biol Cell. 2019;111:217–31.31177551 10.1111/boc.201900012
61. Burkhalter MD Fralish GB Premont RT Caron MG Philipp M Grk5l controls heart development by limiting mTOR signaling during symmetry breaking Cell Rep 2013 4 625 32 10.1016/j.celrep.2013.07.036 23972986
Burkhalter MD, Fralish GB, Premont RT, Caron MG, Philipp M. Grk5l controls heart development by limiting mTOR signaling during symmetry breaking. Cell Rep. 2013;4:625–32.23972986 10.1016/j.celrep.2013.07.036
62. Alupei MC Maity P Esser PR Krikki I Tuorto F Parlato R Loss of proteostasis is a pathomechanism in cockayne syndrome Cell Rep 2018 23 1612 9 10.1016/j.celrep.2018.04.041 29742419
Alupei MC, Maity P, Esser PR, Krikki I, Tuorto F, Parlato R, et al. Loss of proteostasis is a pathomechanism in cockayne syndrome. Cell Rep. 2018;23:1612–9.29742419 10.1016/j.celrep.2018.04.041
63. Frappart PO Walter K Gout J Beutel AK Morawe M Arnold F Pancreatic cancer-derived organoids—a disease modeling tool to predict drug response United Eur Gastroenterol J 2020 8 594 606 10.1177/2050640620905183
Frappart PO, Walter K, Gout J, Beutel AK, Morawe M, Arnold F, et al. Pancreatic cancer-derived organoids—a disease modeling tool to predict drug response. United Eur Gastroenterol J. 2020;8:594–606.10.1177/2050640620905183
64. Cox J Mann M MaxQuant enables high peptide identification rates, individualized p.p.b.-range mass accuracies and proteome-wide protein quantification Nat Biotechnol 2008 26 1367 72 10.1038/nbt.1511 19029910
Cox J, Mann M. MaxQuant enables high peptide identification rates, individualized p.p.b.-range mass accuracies and proteome-wide protein quantification. Nat Biotechnol. 2008;26:1367–72.19029910 10.1038/nbt.1511
65. Cox J Neuhauser N Michalski A Scheltema RA Olsen JV Mann M Andromeda: a peptide search engine integrated into the MaxQuant environment J Proteome Res 2011 10 1794 805 10.1021/pr101065j 21254760
Cox J, Neuhauser N, Michalski A, Scheltema RA, Olsen JV, Mann M. Andromeda: a peptide search engine integrated into the MaxQuant environment. J Proteome Res. 2011;10:1794–805.21254760 10.1021/pr101065j
66. Bentley DR Balasubramanian S Swerdlow HP Smith GP Milton J Brown CG Accurate whole human genome sequencing using reversible terminator chemistry Nature 2008 456 53 9 10.1038/nature07517 18987734
Bentley DR, Balasubramanian S, Swerdlow HP, Smith GP, Milton J, Brown CG, et al. Accurate whole human genome sequencing using reversible terminator chemistry. Nature. 2008;456:53–9.18987734 10.1038/nature07517
67. Kim D Pertea G Trapnell C Pimentel H Kelley R Salzberg SL TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions Genome Biol 2013 14 R36 10.1186/gb-2013-14-4-r36 23618408
Kim D, Pertea G, Trapnell C, Pimentel H, Kelley R, Salzberg SL. TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 2013;14:R36.23618408 10.1186/gb-2013-14-4-r36
68. Zerbino DR Achuthan P Akanni W Amode MR Barrell D Bhai J Ensembl 2018 Nucleic Acids Res 2018 46 D754 61 10.1093/nar/gkx1098 29155950
Zerbino DR, Achuthan P, Akanni W, Amode MR, Barrell D, Bhai J, et al. Ensembl 2018. Nucleic Acids Res. 2018;46:D754–61.29155950 10.1093/nar/gkx1098
69. Smith T Heger A Sudbery I UMI-tools: modeling sequencing errors in Unique Molecular Identifiers to improve quantification accuracy Genome Res 2017 27 491 9 10.1101/gr.209601.116 28100584
Smith T, Heger A, Sudbery I. UMI-tools: modeling sequencing errors in Unique Molecular Identifiers to improve quantification accuracy. Genome Res. 2017;27:491–9.28100584 10.1101/gr.209601.116
70. Liao Y Smyth GK Shi W featureCounts: an efficient general purpose program for assigning sequence reads to genomic features Bioinformatics 2014 30 923 30 10.1093/bioinformatics/btt656 24227677
Liao Y, Smyth GK, Shi W. featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics. 2014;30:923–30.24227677 10.1093/bioinformatics/btt656
71. Love MI Huber W Anders S Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2 Genome Biol 2014 15 550 10.1186/s13059-014-0550-8 25516281
Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15:550.25516281 10.1186/s13059-014-0550-8
72. Jaffe KM Thiberge SY Bisher ME Burdine RD Imaging cilia in zebrafish Methods Cell Biol 2010 97 415 35 10.1016/S0091-679X(10)97022-2 20719283
Jaffe KM, Thiberge SY, Bisher ME, Burdine RD. Imaging cilia in zebrafish. Methods Cell Biol. 2010;97:415–35.20719283 10.1016/S0091-679X(10)97022-2
73. Sokal RR, Rohlf FJ. Biometry: principals and practice of statistics in biological research. San Francisco, California: W. H. Freeman and Co Ltd; 1981.
74. Edgar R Domrachev M Lash AE Gene expression omnibus: NCBI gene expression and hybridization array data repository Nucleic Acids Res 2002 30 207 10 10.1093/nar/30.1.207 11752295
Edgar R, Domrachev M, Lash AE. Gene expression omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002;30:207–10.11752295 10.1093/nar/30.1.207
