A differential requirement for ciliary transition zone proteins in human and mouse neural progenitor fate specification

Wiegering, Antonia; Anselme, Isabelle; Brunetti, Ludovica; Metayer-Derout, Laura; Calderon, Damelys; Thomas, Sophie; Nedelec, Stéphane; Eschstruth, Alexis; Serpieri, Valentina; Catala, Martin; Antoniewski, Christophe; Schneider-Maunoury, Sylvie; Stedman, Aline

doi:10.1038/s41467-025-58554-3

Download PDF

Article
Open access
Published: 05 April 2025

A differential requirement for ciliary transition zone proteins in human and mouse neural progenitor fate specification

Nature Communications volume 16, Article number: 3258 (2025) Cite this article

2923 Accesses
2 Altmetric
Metrics details

Subjects

Abstract

Studying ciliary genes in the context of the human central nervous system is crucial for understanding the underlying causes of neurodevelopmental ciliopathies. Here, we use pluripotent stem cell-derived spinal organoids to reveal distinct functions of the ciliopathy gene RPGRIP1L in humans and mice, and uncover an unexplored role for cilia in human axial patterning. Previous research has emphasized Rpgrip1l critical functions in mouse brain and spinal cord development through the regulation of SHH/GLI pathway. Here, we show that RPGRIP1L is not required for SHH activation or motoneuron lineage commitment in human spinal progenitors and that this feature is shared by another ciliopathy gene, TMEM67. Furthermore, human RPGRIP1L-mutant motoneurons adopt hindbrain and cervical identities instead of caudal brachial identity. Temporal transcriptome analysis reveals that this antero-posterior patterning defect originates in early axial progenitors and correlates with cilia loss. These findings provide important insights into the role of cilia in human neural development.

Enhanced FGFR3 activity in postmitotic principal neurons during brain development results in cortical dysplasia and axonal tract abnormality

Article Open access 28 October 2020

Reconstituting human somitogenesis in vitro

Article 21 December 2022

AUTS2 disruption causes neuronal differentiation defects in human cerebral organoids through hyperactivation of the WNT/β-catenin pathway

Article Open access 22 August 2024

Introduction

Primary cilia are highly conserved signaling organelles that protrude from the surface of most vertebrate cells. Constructed from a microtubule scaffold, the ciliary axoneme is anchored in the plasma membrane via the basal body, a modified mother centriole. A specific structure above the basal body, the ciliary transition zone (TZ), regulates the entry and exit of proteins, thereby acting as a “ciliary gate” that regulates the composition of cilia and creates an independent cell compartment. By modulating several signaling pathways such as SHH, WNT and PDGFR, cilia play essential roles in embryonic development and adult tissue homeostasis. Consequently, defects in primary cilia can lead to a large group of rare genetic diseases called ciliopathies^1,2, which can affect the development of multiple organs, particularly the retina, kidneys, lungs, skeleton, heart and the central nervous system (CNS). Ciliopathies are characterized by considerable phenotypic and genotypic overlap and, rather than being considered separate clinical entities, they are now recognized as a spectrum of genetic diseases. Additionally, there is extensive phenotypic and genetic heterogeneity amongst ciliopathies, and mutations in the same gene can cause strikingly different phenotypes, without necessarily any genotype-phenotype correlation.

Although this is not completely deterministic, mutations in genes encoding ciliary TZ proteins often lead to neurodevelopmental ciliopathies such as Joubert Syndrome (JS) and the fetal lethal Meckel Gruber syndrome (MKS)^1,3. Most of the current knowledge about neurodevelopmental defects in ciliopathy patients and their underlying mechanisms was gained by in vivo studies on animal models and in vitro analyzes. Using these approaches, it was demonstrated that primary cilia are present throughout the neuronal lineage, from neural stem cells to progenitor cells, mature neurons and glia, where they play an essential role in CNS formation by modulating signaling pathways and cell cycle progression^{4,5,6,7,8,9,10}. We and others have demonstrated that the ciliopathy causal gene RPGRIP1L encodes a large scaffolding protein acting as a keystone of TZ construction^{11,12,13,14,15,16,17,18,19}. In addition, it regulates proteasomal activity at the ciliary base and autophagy in murine cells^20,21 and is involved in planar cell polarity in zebrafish and mice²². Furthermore, an essential function for Rpgrip1l in SHH signaling was demonstrated in mouse embryonic tissues and cultured cells^11,15,16,23. However, our understanding of the role of RPGRIP1L and other ciliary proteins in human neural fate is still very partial and has been hampered due to the paucity of human models, limiting our ability to understand the role of cilia in human development, and the origin of neural developmental defects in ciliopathies. Studies using increasingly diverse hiPSC-based organoid models in which the role of ciliary proteins can be investigated are therefore of growing interest to the field. These include studies that investigate ciliary proteins in retinal²⁴, cortical^25,26 or renal²⁷ organoids.

Moreover, evidence suggests that primary cilia function could be only partially overlapping between human and animal models, partly due to specificities in human CNS development. For example, while Rpgrip1l mutant mouse embryos, like other ciliary mutants, display a severe disorganization of the ventral diencephalic and hypothalamic regions and of the ventral spinal cord, no such phenotype was described in ciliopathy patients^6,7,28,29,30.

In this study, we took advantage of pluripotent stem cell-based organoid approaches to study the role of the TZ proteins RPGRIP1L and TMEM67 in human neural progenitor (NP) fate. Both proteins are involved in the gating function of the TZ and have been shown to regulate SHH-dependent spinal patterning in mouse models^8,31,32,33. However, they are structurally and functionally different: RPGRIP1L is a scaffolding protein and a cornerstone of TZ assembly^18,19, while TMEM67 is a Frizzled-like transmembrane receptor^8,31, and the two proteins belong to different TZ protein complexes³⁴. Since the role for ciliary proteins in SHH-dependent patterning of the ventral spinal cord has been largely studied³⁵, we used wild-type (WT) and Rpgrip1l-deficient mouse embryonic stem cell (mESC)- and human induced pluripotent stem cell (hiPSC)-derived spinal organoids to compare the functions of RPGRIP1L in ciliogenesis and spinal progenitors between these two models. We show that mouse Rpgrip1l^-/- organoids exhibit altered dorso-ventral patterning in response to SHH, similar to mutant mouse embryos, emphasizing the cell-intrinsic requirement for Rpgrip1l in spinal progenitors for cilia-dependent SHH pathway activation and commitment into ventral neuronal cell types. In contrast, RPGRIP1L-deficient human iPSCs submitted to a similar differentiation protocol are still able to form primary cilia and adopt SHH-dependent ventral spinal fates. Interestingly, this phenotype is not unique to RPGRIP1L deficiency, but is also observed in TMEM67-deficient human ventral organoids. Remarkably, in addition to the species-specific differences of TZ proteins in cilia formation and signaling between mouse and human, we discovered an unexpected role for the human RPGRIP1L protein in the precise assignment of antero-posterior (rostro-caudal) identities of human spinal progenitors. Indeed, mutant motoneurons (MNs) derived from early neuromesodermal progenitors (NMPs) shift their identity from brachial to more rostral/hindbrain identities. This is accompanied by a temporary loss of cilia in early RPGRIP1L-deficient axial progenitor cells. With this study, we bring a proof of concept that inter-species differences can be dissected by dynamically comparing mouse and human ex vivo organoids and we highlight human-specific functions for the ciliopathy gene RPGRIP1L in antero-posterior patterning of axial progenitors.

Results

Defective SHH-induced motoneuron specification in Rpgrip1l KO mouse spinal organoids

In vertebrates, graded SHH signaling along the dorso-ventral (DV) axis of the spinal cord is necessary for the specification of the floorplate (FP) and the four different ventral progenitor domains p3, pMN, p2 and p1^{36,37,38,39,40,41} (Fig. 1a). The study of mouse ciliary mutant embryos demonstrated that primary cilia are required for SHH transduction and the correct establishment of spinal DV patterning^15,35,42,43. In Rpgrip1l^-/- (KO) mouse embryos, brain and spinal cord progenitors lack a functional cilium, leading to altered SHH pathway activation and repression^6,15,16. As a consequence, in the spinal cord, the FP and p3 territories are not established, fewer pMNs are specified and the dorsal domains are shifted ventrally^6,15. In order to test whether this phenotype can be reproduced ex vivo, even under conditions where SHH activation is forced by the addition of a constant dose of SHH agonist, we established WT and Rpgrip1l KO mESC lines from mouse blastocysts that we differentiated into ventral spinal organoids^44,45,46. Briefly, mESC were differentiated into embryoid bodies (EBs) in N2B27 neuralizing medium for 2 days. EBs were then targeted to a neural spinal fate by combined exposure to the WNT agonist CHIR99021 and to retinoic acid (RA), and further ventralized by addition of SAG (Smoothened AGonist; SHH activation). By day 4, embryoid bodies had organized into spinal organoids that comprised spinal neuronal progenitors from different ventral identities, in particular from the MN lineage. Organoids were let to differentiate until day 7, when most progenitors underwent differentiation (Fig. 1b). WT and Rpgrip1l KO spinal organoids were collected at different time points of differentiation, sectioned and immunostained to assess the identity of the different progenitors obtained. We found that on average, 40% of spinal progenitors from control organoids committed to a pMN fate by day 4 of differentiation, as evidenced by Olig2 staining (Fig. 1c, d). pMNs progressively differentiated into post-mitotic MNs until day 6, a stage when only Islet1-positive MNs could be detected (Fig. 1c). Control organoids were also composed of progenitors from the more ventral domains p3 and FP, since we could detect the presence of Nkx6.1 + /Olig2-, Nkx2.2+ and Foxa2+ cells (Fig. 1d). Moreover, control organoids nicely recapitulated the sequential activation of ventral neural identities observed in vivo in the mouse spinal cord, as Olig2 expression preceded the initiation of Nkx2.2 and Foxa2.

**Fig. 1: Rpgrip1l-deficient mouse spinal progenitors fail to adopt SHH-dependent ventral fates.**

In contrast, Rpgrip1l-deficient progenitors failed to commit to ventral spinal identities. Indeed, the percentage of Olig2+ and Islet1+ cells of the MN lineage as that of Nkx6.1+ cells, was strongly reduced in mutant organoids (Fig. 1c–g), while Nkx2.2+ (p3) and Foxa2+ (FP) progenitor cells were nearly absent from mutant organoids (Fig. 1d). Instead, 80% of the progenitors were positive for Pax6, indicating a shift towards more dorsal identities (Fig. 1d, h). To gain a deeper understanding of the changes in cell fate occurring in Rpgrip1l-deficient organoids, we performed bulk RNAseq on WT and mutant organoids from day 2 to day 7, obtained from 3 independent differentiation experiments (Supplementary Fig. 1a). Our transcriptomic analysis confirmed the shift towards dorsal identities in Rpgrip1l mutant organoids, and showed that SHH-dependent ventral cell types were replaced by progenitors expressing markers characteristics of the p0-p2 domains such as Dbx1 (Fig. 1i and Supplementary Fig. 1b–c). qPCR performed on day 4 organoids confirmed these results and showed a more than 10-fold decrease of Olig2 expression while Dbx1 expression was increased 5 times in mutant compared with control organoids (Fig. 1k–l). Consistently, the molecular signature of post-mitotic neurons changed from V3-MNs in the control, to V1-V0 in Rpgrip1l KO organoids (Fig. 1j). The temporal analysis of the pan-neural markers Sox1 and Tuj1 expression in the time course of the differentiation (Supplementary Fig. 1e–f) showed that the reduction in MNs correlated with a transient decrease in neurogenesis between days 3 and 7 in Rpgrip1l KO compared with control organoids. Since MNs are the first neurons to differentiate in the spinal cord, this result is in line with the dorsalization occurring in the mutant organoids.

In Rpgrip1l mutant mice, the defective patterning of the ventral spinal cord has been attributed to an absence of primary cilia in spinal progenitors, and the consequent defect in SHH pathway activation^6,15. Ensemble Gene Set Enrichment Analysis (EGSEA) performed on differentially expressed genes in day 3 and day 4 organoids showed an enrichment in genes related to spinal patterning (Supplementary Fig. 1c) and SHH signaling (Supplementary Fig. 1d). Consistently, qPCR performed on control and mutant organoids showed that while expression of direct targets of the pathway, Gli1 and Patched1 was upregulated in controls from day 2 upon SAG addition into the EB culture medium, their expression barely increased in the mutants, showing impaired activation of SHH transduction in absence of Rpgrip1l (Fig. 2a, b and Supplementary Fig. 1g, h). To test whether this phenotype was correlated with defective ciliogenesis we analyzed the ciliary marker ADP Ribosylation Factor Like GTPase 13b (Arl13b) by immunostaining on organoid sections. We found that, similar to what was observed in vivo, Rpgrip1l deficiency led to a drastic reduction of Arl13b+ cilia (Fig. 2c, d). This phenotype might reflect a ciliary gating defect, or an absence of ciliary axoneme. To discriminate between these scenarios, we labeled the axonemal protein Ift88 and measured the length of Ift88 staining to evaluate the axonemal length. We found that the axoneme of mutant progenitors was abnormally stunted, with a mean length of 0.2 µm compared to 1 µm in the control (Fig. 2c, e), suggesting that only the TZ might be present above the basal body. This ciliary phenotype was correlated with the complete absence of Rpgrip1l at the ciliary TZ in mutant organoids, consistent with the absence of Rpgrip1l protein in the KO mouse ES cells as shown by western-blot analysis (Fig. 2c, f). Overall, our results demonstrate that, as in vivo, Rpgrip1l is required in vitro in mouse spinal progenitors for cilia structure, SHH pathway activation and the subsequent acquisition of spinal ventral identities. Previous studies showed that spinal organoids can be applied to dissect the activity of the Hedgehog pathway in NPs⁴⁶. Our results provide a proof of concept that they also provide a relevant model for the study of ciliary functions in neural cell fate.

**Fig. 2: Rpgrip1l-deficient mouse spinal organoids fail to activate the SHH pathway and lack functional cilia.**

Human RPGRIP1L-deficient neural progenitors differentiate into MNs in spinal organoids

To tackle the function of primary cilia and the TZ protein RPGRIP1L in human spinal development and SHH transduction, we generated human iPSC clones deficient for RPGRIP1L (RPGRIP1L^-/-) via CRISPR/Cas9-mediated genome editing that lead to premature stop codons in two independent hiPSC lines, PCIi033-A and UCSFi001-A (Supplementary Fig. 2a, b). After genome editing, two RPGRIP1L^+/+ (WT) and two RPGRIP1L^-/- (KO) clones of each hiPSC line (WT: n = 4; KO: n = 4) were validated (Supplementary Fig. 2b–e) and used for differentiation into ventral spinal organoids following a well-established 3D differentiation protocol^47,48 (Fig. 3a). Starting the differentiation of ventral spinal organoids from single cell hiPSCs, EBs were formed after one day in a floating culture in neural induction medium (dual SMAD inhibition) containing the GSK-3β antagonist CHIR99021 for WNT activation and axial progenitor specification^47,48 (Fig. 3a). From day 4 onwards, axial progenitors were treated with RA and high concentrations of SAG in order to induce a ventral neural fate (Fig. 3a). At day 6 of differentiation, the majority of WT and RPGRIP1L-deficient cells expressed SOX2 and PAX6, both markers of early NPs, indicating that control cells as well as RPGRIP1L-deficient cells adopted a neural fate (Supplementary Fig. 3a). In response to SAG stimulation, the pMN marker OLIG2 started to be expressed between day 4 and day 6 of differentiation with an expression peak reached at day 9 where 75% of the cells within organoids expressed OLIG2 (Fig. 3b and Supplementary Fig. 3b). Combined to that, cells expressed the broad ventral marker NKX6.1 at day 9, further evidence of ventral spinal fate specification (Supplementary Fig. 3b). Surprisingly, no difference in OLIG2 and NKX6.1 pattern was observed between control and RPGRIP1L-deficient organoids (Fig. 3b and Supplementary Fig. 3b). In line with this, control and RPGRIP1L-deficient cells started to express ISLET1 and ISLET2, markers for postmitotic MNs, from day 9 onwards (Supplementary Fig. 3c–e). A final treatment of organoids with the γ-Secretase Inhibitor DAPT was used to inhibit NOTCH signaling and accelerate the neural differentiation into postmitotic MNs from day 9 onwards (Fig. 3a). Consequently, ISLET RNA and protein levels were strongly upregulated in WT and RPGRIP1L-deficient organoids between day 11 and day 14 (Fig. 3b and Supplementary Fig. 3c–f). Taken together, these results show that RPGRIP1L-deficient hiPSCs are able to process differentiation from stem cells to NPs and further from pMNs to postmitotic MNs. Nevertheless, our data reveal 15% less ISLET+ MNs at the end of differentiation in RPGRIP1L-deficient hiPSC-derived spinal organoids compared to controls (Fig. 3b and Supplementary Fig. 3f).

**Fig. 3: RPGRIP1L-deficient hiPSC-derived spinal organoids adopt the MN fate.**

By analyzing the expression of various transcription factors that are distinctly expressed in specific DV regions of the developing spinal cord using bulk RNAseq data, we tested whether, despite the normal pMN production, human iPSC-derived organoids deficient for RPGRIP1L exhibited a shift in DV patterning towards the dorsal domains, as it was demonstrated for mESC-derived Rpgrip1l KO spinal organoids (Fig. 1). The unaltered expression of the broad ventral marker NKX6.1, the intermediate markers DBX1 and DBX2, and the dorsal spinal markers PAX3, PAX7, and ZIC1 (Fig. 3c, d) during differentiation, indicated that, contrary to mouse Rpgrip1l-deficient spinal organoids, human RPGRIP1L-deficient spinal organoids did not display a more dorsal progenitor type specification. We noted minor changes among the markers MSX1 and OLIG3 at day 2 and day 4 in RPGRIP1L KO organoids (Fig. 3c), likely reflecting subtle differences in axial progenitors, before DV neural specification.

To investigate whether the mildly reduced number of MNs in RPGRIP1L-deficient organoids could be caused by a delayed neurogenesis, we analyzed the expression of the pan-neuronal markers TUJ1 and MAP2. As expected, expression of both markers increased during differentiation with the highest expression of TUJ1 and MAP2 between day 11 and day 14 (Fig. 3e, f) and high protein levels of TUJ1 at day 14 (Fig. 3g). Comparing WT and RPGRIP1L-deficient organoids, no significant differences between MAP2 and TUJ1 expression on RNA or protein level was observed. In conclusion, neurogenesis appears unaltered and seems not to be responsible for the slightly reduced number of MNs observed in RPGRIP1L-deficient spinal organoids. Surprisingly, NKX2.2 RNA and protein expression was significantly increased at day 11 of MN differentiation in RPGRIP1L-deficient cells (Fig. 3d, h). Since NKX2.2 is expressed in various neural progenitors during differentiation, its increased expression could be caused by generation of different cell types in RPGRIP1L-deficient organoids. For instance, NKX2.2 is expressed in visceral pMNs of the hindbrain that do not express OLIG2^49,50 as well as in a specific subpopulation of ventral spinal pMNs in humans, where it is co-expressed with OLIG2^51,52. Furthermore, NKX2.2 is expressed in the most ventral NPs (p3) and neuronal subtype (V3) within the neural tube (Fig. 1a), where its expression plays a primary role in ventral neuronal patterning and requires high SHH signaling⁵³. To test different hypotheses, we first performed co-labeling of OLIG2 and NKX2.2 at day 9 and day 11 of differentiation. We did not detect any increase in OLIG2 + /NKX2.2+ cells in RPGRIP1L-deficient organoids compared to WT controls (Fig. 3h and Supplementary Fig. 3g). Consequently, the increased NKX2.2 amount in RPGRIP1L-deficient spinal organoids is unlikely to be caused by increased numbers of ventral and human-specific MN subtypes. However, this does not rule out the possibility of an increased production of hindbrain pMNs in RPGRIP1L-deficient spinal organoids. Next, to test whether increased levels of NXK2.2 could be related to the formation of more ventral NPs and spinal cord neurons, we determined SHH pathway activity in WT and RPGRIP1L-deficient spinal organoids by analyzing the expression profiles of SHH target genes GLI1 and PTCH1 along spinal differentiation.

RPGRIP1L deficient spinal neural progenitor possess cilia and transduce SHH signaling

In both WT and RPGRIP1L-deficient spinal organoids, an increase in GLI1 and PTCH1 expression in response to SAG stimulation was observed between day 4 and day 9 of differentiation, followed by an expected decrease in their expression in agreement with the removal of SAG from the differentiation medium (Fig. 4a and Supplementary Fig. 3h, i). No significant difference in SHH target gene expression was detected between WT and RPGRIP1L-deficient organoids. This result is in line with the functional differentiation of RPGRIP1L-deficient hiPSCs into OLIG2+ pMNs and ISLET+ postmitotic MNs that highly depends on SHH signaling transduction. In contrast, the earlier described increase in the number of NKX2.2+ cells in RPGRIP1L-deficient spinal organoids (Fig. 3d, h and Supplementary Fig. 3g) cannot be explained by an increased SHH pathway activity and a ventralization effect of NPs. The remaining hypothesis that could explain the increased amount of NKX2.2 in RPGRIP1L-deficient spinal organoids is therefore the production of visceral NKX2.2+ pMNs.

**Fig. 4: RPGRIP1L-deficient hiPSC-derived neural progenitors harbor cilia and transduce SHH signaling.**

Based on the results of functional SHH transduction and MN differentiation in human RPGRIP1L-deficient spinal organoids, we expected RPGRIP1L-deficient spinal organoids to possess cilia. We analyzed cilia formation at day 6 of MN differentiation, a time point where organoids were stimulated with SAG in order to activate SHH signaling and started to adopt the pMN fate. In WT organoids, short primary cilia were detectable on cells lining the inner lumens of the organoids, thereby pointing into internal cavities (Fig. 4b). These cilia were IFT81-positive and displayed a mean length of 1.3 µm (Fig. 4c, d). RPGRIP1L was present at the ciliary TZ in WT but not RPGRIP1L-deficient spinal organoids (Fig. 4f, g and Supplementary Fig. 2c). Despite the loss of RPGRIP1L, cilia were present lining apical surfaces (Fig. 4b) and displaying an unaltered length and density (Fig. 4d, f) in RPGRIP1L-deficient organoids. However, while WT cilia accumulated the ciliary membrane markers ARL13B, Inositol Polyphosphate-5-Phosphatase E (INPP5E) and Adenylate Cyclase 3 (AC3) (Fig. 4h, i), reduced amounts of ARL13B and AC3 (Fig. 4j, k) but not INPP5E (Fig. 4l) were observed in RPGRIP1L-deficient cilia. This indicates a disrupted ciliary gating in human spinal NPs upon RPGRIP1L deficiency. These results are in line with data from C.elegans, different mouse tissues and human cell lines, in which a disrupted TZ organization and in consequence an altered ciliary gating has been described under Rpgrip1l deficiency^{17,18,19,54,55,56}. Since ciliary gating can be involved in SHH signaling function^2,54,57 e.g., by regulating ciliary shuttling of pathway components such as SMO (Smoothened) thereby regulating ciliary GPR161 (G Protein-Coupled Receptor 161) amounts^58,59, we tested whether GPR161-trafficking was disrupted in RPGRIP1L mutant spinal organoids. GPR161 was enriched in the ciliary membrane in SHH pathway-inactive WT and RPGRIP1L-deficient cells (Supplementary Fig. 3j). In WT and RPGRIP1L mutant NPs stimulated with SAG, GPR161 did not accumulate within primary cilia and thus did not inhibit SHH pathway activation (Supplementary Fig. 3j). These results are in perfect agreement with the previously described functional SHH pathway activation in WT and RPGRIP1L-deficient spinal organoids and show that ciliary gating defects did not affect SHH signaling transduction in RPGRIP1L-deficient NPs.

Given the mild phenotypes observed so far during spinal differentiation of RPGRIP1L^-/- hiPSCs compared to the results obtained in the murine system, we wondered whether mutations could lead to residual RPGRIP1L activity. To validate the phenotype of our RPGRIP1L KO clones, we generated an additional RPGRIP1L-deficient hiPSC clone, this time carrying a full deletion of RPGRIP1L (Δexon3-exon27) (Supplementary Fig. 4a, b). The phenotype of this RPGRIP1L-deletion clone was consistent with our previous results obtained from clones carrying InDel mutations in RPGRIP1L. We observed the loss of RPGRIP1L at the TZ of cilia in NPs (Supplementary Fig. 4c, d) with cilia displaying an unaltered length (Supplementary Fig. 4e) and density (Supplementary Fig. 4f). Furthermore, SHH signaling was unaltered as SHH target gene expression of GLI1 and PTCH1 was induced as efficiently in the full deletion RPGRIP1L clone as in WT controls (Supplementary Fig. 4g, h). ISLET1 expression was highly upregulated in control and RPGRIP1L-deficient organoids at day 14 (Supplementary Fig. 4i) and quantification of ISLET+ cells at day 14 validated the efficient generation of MNs in control and RPGRIP1L-deletion clones (Supplementary Fig. 4j, k). Taken together, the consistency of the results demonstrates the authenticity of the phenotype along our RPGRIP1L mutant spinal organoids.

To test whether the results obtained in RPGRIP1L KO spinal organoids could be extended to other neurodevelopmental ciliopathy genes, we generated TMEM67-deficient hiPSCs via CRISPR/Cas9-mediated genome editing in the PCIi033-A hiPSC line. One control (TMEM67^+/-) and one KO (TMEM67^-/-) clone were used for further analyzes (Supplementary Fig. 5a-d). Cilia were present on control and TMEM67-deficient organoids at day 6 of differentiation where they are supposed to transduce SHH signaling (Supplementary Fig. 5c, e). No alterations of cilia length or quantity have been observed (Supplementary Fig. 5f, g), nonetheless gating defects might be present indicated by lower ciliary ARL13B intensities (Supplementary Fig. 5c). As observed in RPGRIP1L-deficient NPs, gating defects in TMEM67-deficient organoids did not affect SHH transduction demonstrated by unaltered SHH target gene expression of GLI1 and PTCH1 (Supplementary Fig. 5h, i). Furthermore, control and TMEM67 KO organoids started to express ISLET1 between day 9 and day 14 of differentiation (Supplementary Fig. 5j) without significant differences in the percentage of ISLET+ cells per organoid at day 14 (Supplementary Fig. 5k). Thus, human iPSCs mutant for two distinct ciliopathy genes encoding TZ proteins, RPGRIP1L and TMEM67, are able to form cilia, transduce the SHH pathway and differentiate into MNs, strongly suggesting that human and mouse spinal progenitors react differently to the loss of TZ proteins.

RPGRIP1L-deficient spinal organoids show antero-posterior patterning defects with MNs adopting hindbrain identities

To extend the analyzes of WT and RPGRIP1L KO organoids beyond consideration of MN differentiation, we analyzed the bulk RNAseq data in greater detail. Principal Component Analysis (PCA) demonstrated that WT and KO clones clustered closely together at each analyzed time point, indicating low variation between the genotypes (Supplementary Figs. 6a, b). Further testing for overrepresented gene categories via gene ontology (GOSEQ) did not reveal significant changes in biological processes, molecular functions or cellular components between WT and RPGRIP1L-deficient samples at different time points. Among the few genes that were differentially expressed between WT and RPGRIP1L-deficient organoids, expression of HOX genes such as HOXC6 and HOXA7 were significantly downregulated in RPGRIP1L-deficient organoids at day 11 and day 14 (Supplementary Fig. 6c, d, Supplementary Table 3). In addition, ensemble of gene set enrichment analyses (EGSEA) on GeneSetDB gene ontology (GO) at day 2 revealed significant downregulation of the biological process anterior/posterior pattern specification (GO:0009952) involving downregulation of HOX genes as well as the CDX genes CDX1, CDX2 and CDX4 (Supplementary Fig. 6e). In addition to the pairwise comparison of WT and RPGRIP1L KO gene expressions at each selected time point, we performed differential temporal profiles analyzes on our data. These analyzes revealed additional genes whose expression varied significantly over time (Supplementary Fig. 6f), and many of which were related to antero-posterior pattern specification, such as HOXC8 and HOXC9 (Supplementary Figs. 6f, g).

HOX genes are expressed in four clusters along the rostro-caudal axis of bilateral organisms: HOXA, HOXB, HOXC and HOXD. In humans, HOX proteins are encoded by 39 genes whose regionalized expression is realized by their sequential and collinear activation in axial progenitors, which gradually give rise to mesodermal and neuroectodermal structures with rostral and later caudal identity, each expressing a specific HOX code^{48,60,61,62,63,64}. Importantly, upon differentiation HOXB genes become restricted to dorsal spinal neurons and are not expressed in MNs, therefore they were not considered in our analyzes^62,65. By analyzing the dynamic expression pattern of the remaining 3 HOX gene clusters over time, we identified early differences in the expression of posterior HOX genes in RPGRIP1L-deficient organoids compared to controls. In WT, early axial progenitors at day 2 started to express anterior HOX genes such as HOXA1 followed by activation of later and more posterior HOX genes such as HOXA7 at day 4 (Fig. 5a), as it was previously reported⁴⁸. At later stages, WT organoids expressed anterior as well as more posterior HOX genes (Fig. 5a, d and Supplementary Fig. 7a) as HOXA genes from HOXA1 to HOXA7 (Fig. 5a–c and Supplementary Fig. 7b, c) and HOXC genes from HOXC4 to HOXC8 (Fig. 5d-f and Supplementary Fig. 7d, e). Consistently to the HOX gene expression data, MNs within WT organoids adopted a caudal brachial antero-posterior identity at day 14, that could be visualized by immunofluorescence labeling of HOXA7 and HOXC8, markers for brachial MNs⁶² (Fig. 5g, h), whereas more anterior HOX proteins such as HOXA5, a cervical spinal marker^48,62, were not expressed (Fig. 5i). During early stages of differentiation, RPGRIP1L-deficient progenitors expressed anterior HOX genes comparable to the WT situation while expression of the more posterior HOX genes HOXA7 and HOXC8 was already reduced at day 4 (Fig. 5a-f, Supplementary Table 3). This difference was even more striking at later stages, where the latest most caudal HOX genes expressed in RPGRIP1L-deficient MNs were HOXA5 instead of HOXA7 (Fig. 5a, b, c and Supplementary Fig. 7b, c) and HOXC4/HOXC6 instead of HOXC8 (Fig. 5d, e, f and Supplementary Fig. 7d, e). These results were confirmed by immunofluorescence analyzes, which showed that RPGRIP1L-deficient MNs do not express HOXA7 and HOXC8 (Fig. 5g, h), but significantly increased amounts of HOXA5 (Fig. 5i). Taken together, the differences in HOX gene expression indicate a shift in antero-posterior patterning in RPGRIP1L-deficient organoids in which MNs adopt more cervical identities than the caudal brachial identity observed in WT organoids (Fig. 5j) [57, 65, 67].

**Fig. 5: RPGRIP1L-deficient spinal organoids show antero-posterior patterning defects with MNs adopting hindbrain identity.**

We further confirmed the anterior shift in MN identities upon loss of RPGRIP1L by analyzing hindbrain specific pMN and MN markers in WT and RPGRIP1L KO organoids. Expression of PHOX2B, a marker of visceral pMNs and MNs of the hindbrain⁶⁶, and TBX20, a late hindbrain-specific MN marker in mice⁶⁷, was increased in RPGRIP1L KO compared to WT spinal organoids (Fig. 5k, l and Supplementary Fig. 7f). In line with this, the percentage of PHOX2B+ nuclei per organoid was significantly elevated in RPGRIP1L-deficient organoids (Fig. 5m). The results were confirmed in the full deletion RPGRIP1L KO clone (Supplementary Fig. 4l-p), where an increased expression of HOXA5 (Supplementary Fig. 4l) combined with decreased expressions of HOXA7 and HOXC8 was observed (Supplementary Fig. 4l-o). In line with this, the percentage of PHOX2B+ cells per organoid was significantly increased in the full deletion clone of RPGRIP1L (Supplementary Fig. 4p).

Interestingly, the results emphasize the previously discussed hypothesis that the increase of NKX2.2+ cells in RPGRIP1L-deficient organoids (Fig. 3d, h) could be explained by the formation of visceral PHOX2B+ MNs.

Analyzes of early hindbrain and spine-specific markers further confirmed the shift in antero-posterior patterning upon RPGRIP1L loss. The expressions of the hindbrain patterning markers EGR2/KROX20^68,69,70,71, FGF3^71,72,73, PHOX2A^66,68,74 and the visceral pMN marker NKX2.8/NKX2.9^50,75 were strikingly increased in RPGRIP1L-deficient organoids at day 2, day 4 and day 6 (Supplementary Fig. 7g–j). In contrast, expression of WNT5A, involved in axis elongation in mammals⁷⁶, was decreased in early RPGRIP1L KO progenitors (Supplementary Fig. 7k). In line with this, expression of the spine-specific markers CDX1 and CDX4 was strongly reduced in RPGRIP1L-deficient organoids (Supplementary Fig. 7l, m), further emphasizing the caudal to cranial shift in antero-posterior identities of RPGRIP1L KO progenitors^{77,78,79,80,81,82}.

RPGRIP1L controls axial progenitor fate specification and ciliogenesis in early human spinal organoids

The sequential and collinear expression of HOX genes that leads to specification of rostral and later caudal identities is described to be set-up within neuromesodermal progenitors (NMPs)^{60,61,62,63,83,84,85,86,87,88}. It has previously been shown that cells in the human spinal organoid model used in this study transit through an early NMP-like identity, characterized by the co-expression of CDX, BRACHYURY and SOX2^47,48.

To validate that WT and RPGRIP1L-deficient cells adopted the axial NMP-like fate during differentiation, we analyzed the co-expression of BRACHYURY and SOX2 in organoids at day 2 and day 4. SOX2 was equally expressed in WT and RPGRIP1L KO organoids (Supplementary Fig. 8a) with nearly 100% of SOX2+ cells per organoid at day 2 and day 4 (Supplementary Fig. 8c, d, g, h). In contrast, while BRACHYURY expression was equally induced in WT and RPGRIP1L-deficient organoids at day 2 (Supplementary Fig. 8b, c, e, f), its expression dropped as soon as day 4 in RPGRIP1L-deficient cells (Supplementary Fig. 8g, i, j) indicating an altered induction of the NMP-like fate upon RPGRIP1L loss. In addition to BRACHYURY and SOX2, axial NMP-like progenitors in WT organoids expressed the caudal marker CDX2 at day 2 and day 4 of differentiation (Fig. 6a-d). Analyzes of CDX2 gene expression profiles showed its strong upregulation between day 0 and day 2 of differentiation, followed by its downregulation from day 4 onwards (Fig. 6e). In RPGRIP1L-deficient organoids, CDX2 expression was reduced (Fig. 6e, Supplementary Table 3) and significantly less CDX2+ cells per organoid were detectable at day 2 and day 4 of differentiation (Fig. 6a–d), with clusters of CDX2-negative cells in RPGRIP1L-deficient but not WT organoids (Fig. 6c). Together, the results of reduced BRACHYURY and CDX2 expression in KO organoids point to a defective induction of the axial NMP-like cell fate upon RPGRIP1L loss.

**Fig. 6: Altered axial patterning and reduced ciliogenesis in early RPGRIP1L-deficient progenitors.**

We next analyzed the ciliary status of axial progenitors at day 2 and day 4 of differentiation. At day 2, cilia were present in WT organoids in an unpolarized fashion (Fig. 6f–h) whereas cilia density was significantly reduced in RPGRIP1L KO organoids (Fig. 6f, h). The same difference was observed at day 4 of differentiation. Here, cilia were present at the apical sites of WT organoids (Fig. 6i-l), but absent from RPGRIP1L KO cells (Fig. 6i-l).

Since the loss of cilia in early RPGRIP1L-deficient organoids correlated with the stage of NMP-like fate induction and antero-posterior specification of progenitors,we hypothesize that cilia loss on early RPGRIP1L-deficient progenitors is responsible for the reduced specification of the NMP-like axial fate that lead to the antero-posterior patterning defects observed in RPGRIP1L KO organoids. To further investigate this correlation, we analyzed TMEM67-deficient organoids. Cilia were present on control and KO progenitors at day 4 of differentiation (Supplementary Fig. 9a, b), with TMEM67 KO cilia showing the same length and density as controls (Supplementary Fig. 9c, d). Furthermore, the NMP-like fate was correctly acquired in control and TMEM67 KO cells, demonstrated by BRACHYURY and SOX2 expression at day 4 (Supplementary Fig. 9e). Consequently, no antero-posterior patterning defect was observed in TMEM67-deficient organoids at day 14. Expression of the HOX genes HOXA5, HOXA7, HOXC6 and HOXC8 was induced between day 4 and day 14 in control and TMEM67 KO organoids (Supplementary Fig. 9f) and the percentage of HOXA7+ and HOXC8+ cells per organoid was unaltered in TMEM67-deficient organoids compared to controls (Supplementary Fig. 9g). In both control and KO organoids a few HOXA5+ nuclei were present (Supplementary Fig. 9g). Furthermore, low expression of PHOX2B was detected in control and TMEM67 KO organoids at day 14 (Supplementary Fig. 9h), consistent with a few PHOX2B+ cells in control and KO organoids (Supplementary Fig. 9i). These results strengthen the correlation between cilia loss, altered NMP-like identity and perturbed antero-posterior patterning of spinal progenitors.

Interestingly, no spinal-to-cranial shift of progenitor identity was observed in Rpgrip1l^-/- mouse organoids, as evidenced by the unaltered expression of Cdx1-2-4, and Cyp26a1 in mutant progenitors compared to controls (Supplementary Fig. 10a-d). Furthermore, no difference in the efficiency of axial progenitors induction was observed in mutant organoids, where at day 3 similarly to control organoids, the majority of cells was positive for Cdx2 (Supplementary Fig. 10e-f). To gain more insight into the antero-posterior identity of control and Rpgrip1l^-/- mouse spinal progenitors present in our organoids, we also analyzed Hox gene expression by bulk RNAseq and qPCR analyzes. We found that in both genotypes, organoids expressed predominantly brachial and anterior thoracic Hox genes (Supplementary Fig. 10g-k). Moreover, although we noticed a slight decrease in Hox10 gene expression, no major alterations were observed in Hox gene expression that might indicate severe antero-posterior specification defects in the mutants (Supplementary Fig. 10g-k). Since the antero-posterior specification defect in RPGRIP1L KO human organoids was correlated with a reduction in NMP-like cells co-expressing BRACHYURY and SOX2, we analyzed this population in mouse organoids. At day 3, although the percentage of progenitors co-expressing Brachyury and Sox2 was only 20% in average, it did not change between control and mutant organoids (Supplementary Fig. 10l-m), nor did the levels of Brachyury expression (Supplementary Fig. 10n), suggesting that at least part of the progenitors migrate through a bi-potent NMP-like state in both genotypes. Finally, analysis of cilia in organoids at day 3 revealed that mouse axial progenitors deficient for Rpgrip1l presented abnormally shortened cilia (Supplementary Fig. 10o-p), similar to later stages.

Discussion

Previous analysis of Rpgrip1l and Tmem67 mutant mice demonstrated that deficiency for these TZ proteins led to DV mis-patterning of the neural tube marked by the strong decrease of ventral neuronal subtypes like V3 and MNs^6,15,32. These phenotypes were correlated to the absence of primary cilia and functional SHH signaling, and were also described for Ift mutant mice^15,43,89. Since these mice are lacking or displaying a reduced FP, a question remained as to whether altered patterning of the neural tube resulted only from a reduced SHH secretion ventrally, or from impairment of SHH signal transduction within ventral spinal progenitors themselves. Here, by exposing control and Rpgrip1l mutant spinal organoids to high dose of SHH agonist, we show that primary cilia are intrinsically required within mouse spinal progenitors to transduce SHH. Upon inactive SHH signaling, mutant progenitors switch on intermediate-dorsal identities instead of ventral ones.

Strikingly, we find that contrary to the mouse, human RPGRIP1L and TMEM67 are not required for the establishment of ventral spinal neuronal cell types in organoids. Despite a clear ciliary gating defect, with decreased ARL13B or AC3 ciliary localization, human mutant spinal NPs were still able to grow primary cilia and transduce SHH signaling similarly to controls, adopting MN identities. Ciliogenesis in the absence of Rpgrip1l and Tmem67 has also been observed in other cellular contexts, like in mouse limb bud progenitors, MEFs, HEK cells^18,19,20,32 or human fibroblasts derived from RPGRIP1L and TMEM67 mutated JS patients^17,90. Whether the difference in ciliary phenotype observed between mouse and human spinal progenitors stems from differences in ciliary composition or stability, and accounts for the difference in SHH activation in these two systems remains to be addressed. In addition to the species-specific differences in the ciliary phenotype, our results highlight a context-specific role for ciliary TZ proteins in human neural progenitors. In RPGRIP1L-deficient human organoids, NPs at day 6 display cilia of normal length and a moderately altered ciliary content, while at day 2 and day 4, progenitors display a strong reduction of the cilia density. Differences in ciliary structure and composition between tissues, already suggested by the high heterogeneity of phenotypes in human ciliopathies, have been described in other models^18,91,92,93. However, here we uncover a change in cilia sensitivity to the loss of a TZ protein along the process of neural progenitor specification in humans, at a key step of cell fate transition.

Whether this change parallels differences in ciliary structure, and which signaling and molecular cascades lead to these changes in cilia stability during the differentiation process are important questions raised by our study. In this regard, recent data show that cell fate can regulate ciliogenesis by different mechanisms such as transcriptional regulation⁹⁴, niche condition-specific regulation^95,96, and several non-canonical signaling pathway activities such as SHH⁹⁷ and WNT⁹⁸. It will be a key future task to examine cell type specific cilia features and explore whether different modes of ciliogenesis regulation play a role in diverse neural progenitors to explain at least in part the regional specificity of defects found in neurodevelopmental ciliopathies.

Quite unexpectedly, we found that although they do not display SHH-dependent dorsoventral patterning defects, human RPGRIP1L mutant spinal organoids shift to more anterior (hindbrain/cervical spinal) identities as assessed by altered HOX expression patterns and activation of genes specific for hindbrain neurons such as PHOX2B. This switch occurs at early stages of HOX activation, at the axial progenitor stage (day 2 in our spinal organoids) and is paralleled by the above discussed transient loss of cilia and a reduction in the expression of CDX genes. Later (at day 4), reduced numbers of BRACHYURY+/SOX2+ NMP-like cells are found in RPGRIP1L mutants. In TMEM67 mutants, which do not lose cilia at early stages, no perturbation of AP patterning nor NMP-like fate establishment is observed; Thus, although this is not a proof of causal relationship, we find a correlation between cilia loss at axial progenitor and NMP-like stage and the activation of a more anterior progenitor program in spinal organoids. Why TMEM67 depletion does not lead to cilia loss at day 4, while RPGRIP1L depletion does, remains to be understood. As mentioned earlier, these two TZ proteins have different structures and functions and participate in distinct complexes at the TZ^8,18,31,34. Thus, it is possible that they contribute differently to TZ stability and ciliary gating. The analysis of other ciliary gene mutants will be essential to understand the contribution of different TZ proteins in cilia stability.

Another remaining question is the role of Rpgrip1l and more generally cilia, in early antero-posterior specification in the mouse embryo. Our study suggests that, unlike in humans, Rpgrip1l is not required for the specification of the antero-posterior identity of murine axial progenitors. This could reveal an intrinsic difference in the role of Rpgrip1l in ciliary construction and function between human and mouse at the level of axial progenitors. In support of this hypothesis, the ciliary phenotypes observed in axial progenitors appear to differ between the two species as in human organoids, a decrease in cilia density is observed, while in murine organoids cilia are present but truncated. Furthermore, we cannot exclude that variations in the protocols used to generate human and murine organoids may also be at the root of the antero-posterior phenotype differences. In mouse organoids, NMP-like progenitors were observed at a much lower frequency compared to humans, probably because of early addition of retinoic acid at day 2. Such differences in NMP-like cell kinetics might mask the phenotype in mouse organoids. Re-analysis of the phenotype of Rpgrip1l KO mice and other ciliary mutants in the light of these results will be essential to tackle the role of cilia in antero-posterior regionalization of the embryonic axis in an in vivo context. The altered antero-posterior patterning that was observed in human RPGRIP1L KO organoids could arise from several mechanisms. WNT, RA, Fibroblast Growth Factor (FGF), and more recently NOTCH signaling, have been involved in the determination of axial progenitors and thereby regulation of gradual HOX gene expression^{47,64,77,99,100,101,102,103}. From our transcriptome data, we could not observe any indication of perturbations in FGF or NOTCH pathways upon RPGRIP1L deficiency. Analysis of Wnt3a- and Cdx-deficient mouse embryos^78,104,105 or experimental manipulation of CDX or BRACHYURY activity^{106,107,108,109} lead to a working model in which early BRACHYURY and CDX expression activates and fine-tunes HOX gene expression profiles along the antero-posterior axis downstream of WNT/β-Catenin signaling in axial progenitors^{77,79,83,85,99,100,104,106,109,110,111,112}. Interestingly, it was demonstrated recently that WNT signaling during early differentiation of hiPSCs into spinal MNs allows progressive, collinear HOX activation and that a switch to RA treatment at any time point during this process stops collinear HOX activation and favors transition into neuroectodermal progenitors with a discrete HOX gene profile^48,100. Thus, we hypothesize that the shift in antero-posterior patterning that we describe in RPGRIP1L-deficient human spinal organoids results from reduced WNT transduction during early stages of axial progenitor differentiation, leading to a reduction in BRACHYURY, CDX2, CDX1 and CDX4 expression and the consequent decreased expression of posterior HOX genes, so that MNs adopt hindbrain or cervical rather than caudal brachial identities. Whether RPGRIP1L or cilia act directly on maintaining WNT signaling in axial progenitors remains to be demonstrated. While cilia and TZ proteins, and specifically RPGRIP1L and TMEM67, have been involved in modulating WNT pathways, their effect is complex and depends on the tissue and cell type^{22,113,114,115,116,117}. Further analyzes will show if and how cilia regulate WNT signaling and consequently BRACHYURY, CDX and HOX gene expression during axial progenitor fate specification, or whether other pathways such as FGF and RA are implicated.

Another important question raised by our study is whether the defect in antero-posterior patterning upon loss of RPGRIP1L could be related to phenotypic features of neurodevelopmental ciliopathies. Since the antero-posterior patterning defect of MNs described here most likely lead to a mild shift of boundaries between the posterior hindbrain region and the anterior spinal region with an increase in anterior identities, a corresponding phenotype could be present in ciliopathy patients. Unfortunately, analyzes of such patterning defects in patients are not available. A first important step would be therefore to test the differentiation capacity of patient-derived hiPSCs and analyze their phenotype related to antero-posterior patterning. Apart from that, our study provides interesting entry points for further investigations where it will be important to analyze antero-posterior patterning more carefully and in additional CNS regions such as the hindbrain, as it is highly sensitive to antero-posterior patterning defects and one of the major targets of neurodevelopmental ciliopathies. Here, the increased expression of PHOX2B appears to be particularly interesting, since its deregulation was shown to cause breathing pattern regulation defects^{118,119,120,121,122}, a feature that can also be observed in JS patients¹²³.

Methods

Ethics

The work performed in this manuscript complies with all relevant ethical regulations. Mouse breeding, mating and embryo collection were performed under ethical agreement number APAFIS-1382 from the French Ministry of Higher Education and Research. The obtaining and manipulation of human iPSC lines were performed under ethical agreements from the French Ministry of Higher Education and Research DC-2021-4734 (for the UCSFi001-A line), DC-2015-2595 and DC-2016-A00773-480 (for the PCIi033-A line) and GMO authorization number DUO-9509.

Generation of mESCs lines

Rpgrip1l^-/- and Rpgrip1l^+/+ mouse ES cell lines were derived from blastocysts according to standard procedures using 2-4 months-aged male and female Ftm-deficient mice produced and genotyped as described previously¹⁵ and maintained in the C57BL/6 J background. One ES clone from each genotype was used in this study. Both Rpgrip1l^-/- and Rpgrip1l^+/+ clones are males, which avoids sex-dependent bias in RNAseq data.

mESCs culture and differentiation into ventral spinal organoids

mESCs were cultured in DMEM-GlutaMAX (Gibco) supplemented with 15% FCS, 100 µM β-mercaptoethanol, 100 U/ml penicillin, 0.1 mg/ml streptomycin, and 10³ U/ml LIF (Millipore). To generate mouse spinal organoids, low-passage (inferior to 10) mESCs were dissociated and differentiated into embryoid bodies by seeding them in ultra-low attachment petri dishes (Corning) at a density of 9 × 10⁴ cells ml⁻¹ in N2B27 differentiation medium containing Advanced DMEMF12 and Neurobasal media (1:1, Life Technologies), B27 w/o vitamin A (Life Technologies), N2 (Life Technologies), 2 mM L-glutamine (Life Technologies), 0.1 mM β-mercaptoethanol, penicillin and streptomycin (Life Technologies). On day 2, the medium was changed and supplemented with 3 µM CHIR, 10 nM RA and 500 nM SAG to induce ventral spinal fate. The medium was changed daily, with addition of 10 nM RA only.

hiPSCs culture

hiPSCs were thawed in presence of Rock-inhibitor Y-27632 (5 µM, Stemgent-Ozyme #04-0012) and cultured under standard conditions at 37 °C in mTeSR+ medium (Stem Cell Technologies #100-0276) on Matrigel (Corning, VWR #354277) coated plates upon confluency of 70-80%. Passaging was performed using ReLeSR (Stem Cell Technologies #05872) and testing for potential mycoplasma contamination was performed regularly by Eurofins genomic (Mycoplasmacheck). Accutase (Stem Cell Technologies #07920) was used for the dissociation of hiPSC colonies into single cells.

Generation of mutant hiPSC lines

RPGRIP1L mutant and wildtype hiPSC clones were produced via CRISPR/Cas9-mediated genome editing in two different hiPSC lines: PCIi033-A (PHENOCELL (PLI), product reference PCi_CAU2) and UCSFi001-A (deposited at the Coriell Institute for Medical Research under the identifier GM25256). Both cell lines are reported to be males. The targeted exon 3 of RPGRIP1L and flanking sequences were sequenced prior to gene editing in both cell lines (see on- and off-target analyzes section). The gRNA targeting exon 3 of RPGRIP1L was designed using the online design tool CRISPOR¹²⁴. Annealing of crRNA (Alt-R® CRISPR-Cas9 crRNA, 400 nM, IDT) and tracrRNA (Alt-R® CRISPR-Cas9 tracrRNA, 400 nM, IDT #1072533) was performed and 400.000 hiPSCs of each line were transfected with 1 µl Cas9 protein (30 µM; kindly provided by TACGENE, Paris, France) and 1 µl gRNA (200 µM) using nucleofection. The transfected cells were cultured in mTeSR+ medium (Stem Cell Technologies #100-0276) containing Rock-inhibitor Y-27632 (10 µM, Stemgent-Ozyme #04-0012) for one day, replaced by mTeSR+ medium without Rock-inhibitor. Medium was changed every second day and serial dilutions were performed in mTeSR+ medium containing Rock-inhibitor Y-27632 and CloneR (1:1000, Stem Cell Technologies #05889). Rock-inhibitor was removed after two days and the medium containing CloneR was changed every second day until colonies were ready to be picked. Colonies were isolated in 24-well plates in mTesR+ medium containing CloneR. DNA was isolated for genotype analyzes (see on- and off-target analyzes section). Four WT (two clones of each hiPSC line) and four homozygous mutant hiPSC clones (two of each hiPSC line) carrying InDel mutations that lead to premature STOP were used for further analyzes (Supplementary Fig. 2).

A large deletion of RPGRIP1L was induced in the PCli033-A line by combining the gRNA targeting exon 3 of RPGRIP1L with a gRNA targeting exon 27 of RPGRIP1L (Δexon3-exon27). This approach was performed at the ICV-iPS platform of the Brain Institute in Paris (ICM, Paris). One RPGRIP1L KO clone and one heterozygous control clone were used for further analyzes (Supplementary Fig. 3 and Supplementary Fig. 11). No WT clone was recovered in this experiment. For the generation of TMEM67-deficient hiPSC clones from the PCli033-A cell line, the following strategy was used: The gRNA targeting the 5’UTR of TMEM67 was combined with a gRNA targeting intron 27 of TMEM67 to generate a large deletion of TMEM67 (Δexon1-exon27). One TMEM67 KO clone and one heterozygous control clone were used for further analyzes (Supplementary Fig. 5 and Supplementary Fig. 12). No WT clone was recovered in this experiment. For all clones, genomic stability was assessed by detection of recurrent genetic abnormalities using the iCS-digitalTM PSC test, provided as a service by Stem Genomics (https://www.stemgenomics.com), as described previously¹²⁵.

On- and off-target analyzes

Genomic DNA was extracted from both hiPSC lines (PCli033-A and UCSFi001-A) prior to CRISPR/Cas9 gene editing. Genotyping of on-targets was performed by PCR, and PCR products were sent for sequencing. After CRISPR/Cas9 gene editing, genomic DNA was isolated from individual hiPSC clones and the on-target sites were amplified by PCR. The PCR products were sent for sequencing to determine the specific mutations that were induced. Off-targets were identified using the CRISPOR online design tool¹²⁶ and ranked according to their CFD off-target score¹²⁷. Potential off-targets within exons and with CFD scores > 0.02 were selected for PCR and subsequent sequencing analyzes¹²⁸. The primers used for on- and off-target analyzes are listed in Supplementary Table 1.

Differentiation of hiPSCs into spinal organoids

Differentiation of wild-type and mutant hiPSCs was performed as previously described. After amplification, hiPSC lines were dissociated into single cells using Accutase (Stem Cell Technologies #07920) and resuspended in differentiation medium N2B27 [vol:vol; Advanced DMEM/F-12 (Gibco) and Neurobasal Medium (Gibco)] supplemented with N2 (Thermo Fisher #17502048), B27 without Vitamin A (Thermo Fisher #12587010), penicillin/streptomycin 1% (Thermo Fisher #15140122), β-mercaptoethanol 0.1% (Thermo Fisher #31350010). Cells were seeded in ultra-low attachment dishes (Corning #3261) to allow EB formation. The N2B27 differentiation medium was used throughout the whole differentiation process (Fig. 3a). Rock-inhibitor Y-27632 (5 μM; Stemgent-Ozyme #04-0012) was added from day 0 to day 2, CHIR-99021 (3 µM; Stemgent-Ozyme #04-0004) from day 0 to day 4, SB431542 (20 μM; Stemgent-Ozyme #04-0010) from day 0 to day 3 and LDN 193189 (0.1 μM; Stemgent-Ozyme #04-0074) from day 0 to day 4. The differentiation was preceded by adding Smoothened AGonist (SAG 500 nM; Merck #566660) and Retinoic Acid (100 nM RA; Sigma #R2625) to the N2B27 differentiation medium from day 4 to day 9. γ-Secretase inhibitor DAPT (10 µM; Tocris Bioscience #2634) was added from day 9 to the end of differentiation. Medium was changed every other day and additionally at day 3 of differentiation.

EB embedding and Cryosectioning

Murine and human spinal organoids were collected at different time points during differentiation, rinsed with PBS and fixed in cold PFA (4%) for 7–12 min at 4 °C. EBs were rinsed in PBS and incubated in 30% sucrose in PBS until completely saturated. Cryoprotected EBs were embedded in OCT embedding matrix (Cell Path #KMA-0100-00A) and stored at -80 °C. 12 µm cryostat sections were prepared and stored at -80 °C.

Western Blot

Mouse embryonic stem cells were collected during amplification and lysis was performed using RIPA buffer (150 mM NaCl, 25 mM Tris pH 7.5, 1% NP40, 0.1% SDS, Sodium Deoxycholate 1%) containing proteases inhibitors. Protein samples in 4x sample buffer were loaded on precast gels (Mini-PROTEAN TGX Gels, BioRad #4561093) and separated by SDS-PAGE in the Mini-PROTEAN Tetra Cell system (BioRad). The Spectra Multicolor High Range Protein Ladder (Thermo Fisher #26625) was used as a molecular weight marker. Blotting was performed on PVDF membranes (BioRad #162-0176) using the Mini Trans-Blot Cell tank blotting protocol (BioRad). Primary antibodies against Rpgrip1l and Actin, and secondary antibodies conjugated to horseradish peroxidase were used (Supplementary Table 2). Visualization of protein bands was realized using the ECL detection kit (SuperSignal WestPico PLUS Chemiluminescent substrate, ThermoScientific #34577) and the ChemiDoc MP Imaging System (BioRad).

Immunofluorescence

Cryostat sections of murine and human organoids were washed in PBS/0.1% Triton X-100. Optional steps for permeabilization with PBS/0.5 %Triton X-100 for 10 min, post-fixation with MeOH (100%) for 10 min or antigen-retrieval with citrate-buffer were performed depending on the antibody and specimen. After washing, cryosections were incubated with a blocking solution containing 10% NGS in PBS/0.1% Triton-X100 for at least 2 hours. The sections were incubated with primary antibodies (Supplementary Table 2) diluted in PBS/0.1% Triton-X100 and 1% NGS overnight at 4 °C. Next day, the slides were washed 3 × 5 min with PBS/0.1% Triton-X100 and incubated with the secondary antibodies (Supplementary Table 2) diluted in PBS/0.1% Triton X-100 and 1% NGS for 2 hours. After final washing steps with PBS/0.1%Triton X-100 for at least 1 hour, sections were mounted with Vecta-Shield (Vector #H-1000) or Mowiol (Roth #0713.2) optionally containing DAPI (Merck #1.24653).

Image acquisition

Fluorescence image acquisition of patterning and differentiation markers was performed at room temperature using a Zeiss Observer Z1 equipped with an Apotome module, an Axiocam 506 monochrome cooled CCD camera and a 40x oil objective with a NA 1.3. All images were processed with Zen software (blue edition). Image analyzes and figure preparation were performed using Fiji (ImageJ; National Institutes of Health). For cilia analyzes, fluorescence images were acquired at room temperature using an inverted confocal microscope (LSM 710, Carl Zeiss AG or TCS SP5, Leica), a 63x oil objective with a NA 1.4, and a monochrome charge-coupled device camera. Z-stacks with 0.2 µm steps were acquired. Images were processed with Zen software (black edition). Image analyzes and figure preparation were performed using Fiji (ImageJ; National Institutes of Health).

STED images were acquired on the STEDYCON confocal and STED module (Abberior Instruments GmbH, Germany) connected to the lateral port of an IX83 inverted microscope (Evident Scientific). The STEDYCON SmartControl software was used for data acquisition. Imaging was performed using a UPLXAPO 100x oil immersion objective (Evident Scientific). For STED, a pulsed 775 nm laser was used. STAR RED was excited with a pulsed 640 nm laser, STAR ORANGE was excited with a pulsed 561 nm laser, and DAPI was excited using a continuous wave 405 nm source. Pixel size was 20 nm in xy direction, with a 250 nm step size in the z direction.

Quantification of fluorescent staining and statistical analyzes

The quantitative analysis of nuclear stainings was done either by scoring the percentage of positive nuclei for the selected markers, or by measuring the stained area. The percentage of positive nuclei for selected stainings was determined by first applying automated nuclear segmentation using the cellpose python environment (https://www.cellpose.org/), followed by manual analysis and scoring of the stained nuclei. The positive stained areas were defined and measured via setting common thresholds in Fiji. At least 10 organoids per clone (n = number of clones) and experiment (N = number of experiments) were analyzed for quantifications. Intensity measurements of ciliary proteins were performed on unprocessed images (raw data) as described before^{18,19,20,21,54}. Analyzes of ciliary density were performed by manually counting cilia axonemes over the number of basal bodies (γ-Tubulin). Whenever possible, at least 100 cilia per clone were used for intensity, length and density measurements. After quantifications were performed, representative images were processed by means of background subtraction and contrast settings via Adobe Photoshop CS2.

Data are presented as mean ± SEM or as median with quartiles. Two-tailed t test with Welch’s correction was performed for all data in which two normally distributed datasets were compared. Mann-Whitney tests were performed for data sets that were not normally distributed. Two-tailed t test Analysis of variance (ANOVA) and Tukey honest significance difference (HSD) tests were used for all data in which more than two normally distributed datasets were compared. Kruskal-Wallis with Dunn’s multiple comparison test was used for comparison of more than two datasets with small or not normally distributed data. Two-sided tests with confidence intervals of 5% were used for all statistical analyzes. Following statistical significances were considered: *p < 0.05, **p < 0.01, ***p < 0.001 and ****p < 0.0001. All statistical data analysis and graph illustrations were performed using GraphPad Prism (GraphPad Software).

RNA isolation and qPCR analyzes

mESC- and hiPSC-derived organoids were collected at different stages during differentiation and washed with PBS. RNA was isolated using the RNeasy Kit (Qiagen #74104) and RNAse-free DNase Set (Qiagen #79254). Isolated RNA was transcribed into cDNA using the GoScript Reverse Transcriptase Kit (Promega #A5001). For quantitative real-time PCR, 50 ng of cDNA of each sample was used in a Maxima SYBR Green/ROX qPCR Master Mix 2x (Thermo Scientific #K0222). Reactions were run in a Step One Real-Time PCR System Thermal Cycling Block (Applied Biosystems #4376357). Primer pairs are listed in Supplementary Table 1. The analysis of real-time data was performed using the included StepOne Software version 2.0 and graphs were generated using GraphPad Prism (GraphPad Software).

Bulk RNAseq analysis on mouse and human spinal organoids

mESC- and hiPSC-derived organoids were collected at different stages during differentiation and washed with PBS. For mouse organoids, pooled organoids from Wt and Rpgrip1l KO were collected from 3 independent experiments (N = 3). For human organoids, material from two control clones of the PCIi033-A hiPSC line and the UCSFi001-A hiPSC line each (WT: n = 4) and from one RPGRIP1L-deficient clone of the PCIi033-A hiPSC line and the UCSFi001-A hiPSC line each (KO: n = 2) was collected. RNA was isolated using the RNeasy Kit (Qiagen #74104) and RNAse-free DNase Set (Qiagen #79254). RNA quality was determined by measuring the RNA integrity number (RIN) via the High Sensitive RNA Screen Tape Analysis kit (Agilent Technologies #5067) on the TapeStation system (Agilent Technologies). A RIN above 9.5 was considered as good quality RNA and 250 ng RNA in a total volume of 25 µl was prepared per sample for further procedure. Bulk RNAseq was realized by the Genotyping and Sequencing Core Facility of the Paris Brain Institute (iGenSeq, ICM Paris). RNAseq data processing was performed in Galaxy under supervision of the ARTBio platform (IBPS, Paris).

Paired-end RNA reads were aligned against the M. musculus mm39 and the homo sapiens hg38 genome, respectively, by using the STAR aligner (v2.7.10, Galaxy wrapper v4)¹²⁹ and the gene annotations gtf files GRCm39.109 and GRCh38.109, respectively. Quality of sequencing was controlled with fastQC (Galaxy wrapper v0.74+galaxy0) and MultiQC (Galaxy wrapper v1.11+galaxy1)¹³⁰. Gene expression was assessed by counting sequencing reads aligned to gene exons with featureCounts (Galaxy wrapper v2.0.3+galaxy2)¹³¹. Raw counts were further converted to normalized gene expression values using the log2(tpm+1) transformation where tpm is the count of transcript aligned reads per length of transcript (kb) per million of total mapped reads (Galaxy tool “cpm_tpm_rpk” v0.5.2). Principal Component Analyzes (PCA) were performed using r-factominer (v.2.9, Galaxy tool id “high_dimensions_visualisation” v4.3+galaxy0)¹³² and heatmap visualizations were produced with the Galaxy tool “high_dim_heatmap” (v3.1.3+galaxy0) using the normalized gene expression values. Differentially expressed genes (DEGs) were selected from the gene raw read counts using DESeq2 (Galaxy wrapper v2.11.40.8+galaxy0)¹³³ and the Benjamini-Hochberg p-adjusted cutoff 0.01. DESeq2 statistical tables were used for generation of Volcano Plots (Galaxy tool id “volcanoplot” v0.0.5)¹³⁴, for goseq tests of overrepresented gene ontology categories (Galaxy wrapper v1.50.0+galaxy0)¹²⁴, as well as for ensemble of gene set enrichment analyzes (EGSEA) (Galaxy wrapper v1.20.0)¹³⁵.

Analysis of differential temporal profiles on human spinal organoids

We obtained 62710 gene temporal profiles of expression levels [expressed in log2(tpm+1)], each composed of 7 time points in the WT and RPGRIP1L KO mutant contexts. We selected genes whose variance of expression across the 14 WT and KO time points is higher than 0.1, and we scaled each gene time point (WT and KO) using the minmax function on the TSPred R package, following which the scaled expression in each series is E’ = (E - E_min) / (E_max - E_min), the highest expression of series is 1 and the lowest expression of series is 0. Note that since we scaled wild type and mutant expressions together, their relationship is maintained during the minmax scaling.

We next split the data into gene time series of 7 points in WT or KO context, and for each gene, we computed the euclidean distance between the wild type and mutant series. The density function of the euclidean distance was empirically derived from the observed distances. Using this function, we selected genes with the 1% highest distances between time profiles in wild type or mutant context. Time profiles of selected genes were finally plotted using their unscaled expression values. Since the euclidean distance is not able to capture significant differences between all time series, we repeated the above analysis with 11 other distance metrics, including short time series (sts), dynamic time warping (dtw), time alignment measurement (tam), autocorrelation-based dissimilarity (acf), fourier discrete transform (fourier), compression-based dissimilarity (cdm), complexity-invariant (cid), Pearson’s correlation dissimilarities (cor), integrated periodogram-based dissimilarity (int.per), periodogram-based dissimilarity (per) and Frechet (frechet) distances.

Gene time profiles selected using the various distances metrics are compiled in table [distance mode] which indicates for each gene the metrics that returned a significant distance. In total, we selected 542 genes with significant distance between mutant and wild-type time profiles. These profiles are sorted according to the number of metrics returning significant distance (p < 0.01). R code used for this analysis is available here for download.

Reporting summary

Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.

Data availability

The RNA sequencing data for this study have been deposited in the European Nucleotide Archive (ENA) at EMBL-EBI under accession number PRJEB86489. Source data are provided with this paper.

References

Bachmann-Gagescu, R. et al. Healthcare recommendations for Joubert syndrome. Am. J. Med. Genet. A 182, 229–249 (2020).
Article PubMed Google Scholar
Reiter, J. & Leroux, M. Genes and molecular pathways underpinning ciliopathies. Nat. Rev. Mol. Cell Biol. 18, 533–547 (2017).
Article CAS PubMed PubMed Central MATH Google Scholar
Brancati, F., Dallapiccola, B. & Valente, E. Joubert Syndrome and related disorders. Orphanet J. Rare Dis. 5, 20 (2010).
Article PubMed PubMed Central Google Scholar
Goetz, S. & Anderson, K. The primary cilium: a signalling centre during vertebrate development. Nat. Rev. Genet. 11, 331–344 (2010).
Article CAS PubMed PubMed Central MATH Google Scholar
Spassky, N. et al. Primary cilia are required for cerebellar development and Shh-dependent expansion of progenitor pool. Dev. Biol. 317, 246–259 (2008).
Article CAS PubMed PubMed Central Google Scholar
Andreu-Cervera, A., Catala, M. & Schneider-Maunoury, S. Cilia, ciliopathies and hedgehog-related forebrain developmental disorders. Neurobiol. Dis. 150, 105236 (2021).
Article CAS PubMed Google Scholar
Andreu-Cervera, A. et al. The ciliopathy gene ftm/rpgrip1l controls mouse forebrain patterning via region-specific modulation of hedgehog/gli signaling. J. Neurosci. 39, 2398–2415 (2019).
PubMed PubMed Central MATH Google Scholar
Abdelhamed, Z. et al. The ciliary Frizzled-like receptor Tmem67 regulates canonical Wnt/β-catenin signalling in the developing cerebellum via Hoxb5. Sci. Rep. 9, 5446 (2019).
Article ADS PubMed PubMed Central MATH Google Scholar
Lancaster, M. et al. Defective Wnt-dependent cerebellar midline fusion in a mouse model of Joubert syndrome. Nat. Med. 17, 726–731 (2011).
Article CAS PubMed PubMed Central MATH Google Scholar
Basten, S. & Giles, R. Functional aspects of primary cilia in signaling, cell cycle and tumorigenesis. Cilia 2, 6 (2013).
Article CAS PubMed PubMed Central MATH Google Scholar
Delous, M. et al. The ciliary gene RPGRIP1L is mutated in cerebello-oculo-renal syndrome (Joubert syndrome type B) and Meckel syndrome. Nat. Genet. 39, 875–881 (2007).
Article CAS PubMed Google Scholar
Arts, H. et al. Mutations in the gene encoding the basal body protein RPGRIP1L, a nephrocystin-4 interactor, cause Joubert syndrome. Nat. Genet. 39, 882–888 (2007).
Article CAS PubMed MATH Google Scholar
Doherty, D. et al. Mutations in 3 genes (MKS3, CC2D2A and RPGRIP1L) cause COACH syndrome (Joubert syndrome with congenital hepatic fibrosis). J. Med. Genet. 47, 8–21 (2010).
Article CAS PubMed Google Scholar
Wolf, M. et al. Mutational analysis of the RPGRIP1L gene in patients with Joubert syndrome and nephronophthisis. Kidney Int 72, 1520–1526 (2007).
Article CAS PubMed MATH Google Scholar
Vierkotten, J., Dildrop, R., Peters, T., Wang, B. & Rüther, U. Ftm is a novel basal body protein of cilia involved in Shh signalling. Development 134, 2569–2577 (2007).
Article CAS PubMed Google Scholar
Besse, L. et al. Primary cilia control telencephalic patterning and morphogenesis via Gli3 proteolytic processing. Development 138, 2079–2088 (2011).
Article CAS PubMed MATH Google Scholar
Shi, X. et al. Super-resolution microscopy reveals that disruption of ciliary transition-zone architecture causes Joubert syndrome. Nat. Cell Biol. 19, 1178–1188 (2017).
Article CAS PubMed PubMed Central MATH Google Scholar
Wiegering, A. et al. Cell type-specific regulation of ciliary transition zone assembly in vertebrates. EMBO J. 37, e97791 (2018).
Article PubMed PubMed Central Google Scholar
Wiegering, A. et al. Rpgrip1l controls ciliary gating by ensuring the proper amount of Cep290 at the vertebrate transition zone. Mol. Biol. Cell 32, 675–689 (2021).
Article CAS PubMed PubMed Central Google Scholar
Gerhardt, C. et al. The transition zone protein Rpgrip1l regulates proteasomal activity at the primary cilium. J. Cell Biol. 210, 115–133 (2015).
Article CAS PubMed MATH Google Scholar
Struchtrup, A., Wiegering, A., Stork, B., Rüther, U. & Gerhardt, C. The ciliary protein RPGRIP1L governs autophagy independently of its proteasome-regulating function at the ciliary base in mouse embryonic fibroblasts. Autophagy 14, 567–583 (2018).
Article CAS PubMed PubMed Central Google Scholar
Mahuzier, A. et al. Dishevelled stabilization by the ciliopathy protein Rpgrip1l is essential for planar cell polarity. J. Cell Biol. 198, 927–940 (2012).
Article CAS PubMed PubMed Central MATH Google Scholar
Gerhardt, C., Wiegering, A., Leu, T. & Rüther, U. Control of Hedgehog signalling by the cilia-regulated proteasome. J. Dev. Biol. 4, 27 (2016).
Article PubMed PubMed Central Google Scholar
McDonald, A. & Wijnholds, J. Retinal ciliopathies and potential gene therapies: a focus on human ipsc-derived organoid models. Int. J. Mol. Sci. 25, 2887 (2024).
Article CAS PubMed PubMed Central MATH Google Scholar
Schembs, L. et al. The ciliary gene INPP5E confers dorsal telencephalic identity to human cortical organoids by negatively regulating Sonic hedgehog signaling. Cell Rep. 39, 110811 (2022).
Article CAS PubMed PubMed Central Google Scholar
Boutaud, L. et al. 2D and 3D human induced pluripotent stem cell-based models to dissect primary cilium involvement during neocortical development. J. Vis. Exp. 181, https://doi.org/10.3791/62667 (2022).
Cruz, N. et al. Modelling ciliopathy phenotypes in human tissues derived from pluripotent stem cells with genetically ablated cilia. Nat. Biomed. Eng. 6, 463–475 (2022).
Article CAS PubMed PubMed Central MATH Google Scholar
Willaredt, M. et al. A crucial role for primary cilia in cortical morphogenesis. J. Neurosci. 28, 12887–12900 (2008).
Article CAS PubMed PubMed Central MATH Google Scholar
Juric-Sekhar, G., Adkins, J., Doherty, D. & Hevner, R. Joubert syndrome: brain and spinal cord malformations in genotyped cases and implications for neurodevelopmental functions of primary cilia. Acta Neuropathol. 123, 695–709 (2012).
Article PubMed Google Scholar
ten Donkelaar, H., Hoevenaars, F. & Wesseling, P. A case of Joubert’s syndrome with extensive cerebral malformations. Clin. Neuropathol. 19, 85–93 (2000).
PubMed MATH Google Scholar
Abdelhamed, Z. et al. The Meckel-Gruber syndrome protein TMEM67 controls basal body positioning and epithelial branching morphogenesis in mice via the non-canonical Wnt pathway. Dis. Model Mech. 8, 527–541 (2015).
Article CAS PubMed PubMed Central Google Scholar
Abdelhamed, Z. et al. Variable expressivity of ciliopathy neurological phenotypes that encompass Meckel-Gruber syndrome and Joubert syndrome is caused by complex de-regulated ciliogenesis, Shh and Wnt signalling defects. Hum. Mol. Genet. 22, 1358–1372 (2013).
Article CAS PubMed PubMed Central Google Scholar
Adams, M. et al. A meckelin–filamin A interaction mediates ciliogenesis. Hum. Mol. Genet. 21, 1272–1286 (2012).
Article CAS PubMed MATH Google Scholar
Czarnecki, P. G. & Shah, J. V. The ciliary transition zone: from morphology and molecules to medicine. Trends Cell Biol. 22, 201–210 (2012).
Article CAS PubMed PubMed Central MATH Google Scholar
Bangs, F. & Anderson, K. Primary cilia and mammalian hedgehog signaling. Cold Spring Harb. Perspect. Biol. 9, a028175 (2017).
Article PubMed PubMed Central Google Scholar
Briscoe, J. & Ericson, J. Specification of neuronal fates in the ventral neural tube. J. Curr. Opin. Neurobiol. 11, 43–49 (2001).
Article CAS MATH Google Scholar
Balaskas, N. et al. Gene regulatory logic for reading the Sonic Hedgehog signaling gradient in the vertebrate neural tube. Cell 148, 273–284 (2012).
Article CAS PubMed PubMed Central Google Scholar
Ribes, V. et al. Distinct Sonic Hedgehog signaling dynamics specify floor plate and ventral neuronal progenitors in the vertebrate neural tube. Genes Dev. 24, 1186–1200 (2010).
Article CAS PubMed PubMed Central MATH Google Scholar
Dessaud, E. et al. Dynamic assignment and maintenance of positional identity in the ventral neural tube by the morphogen sonic hedgehog. PLoS Biol. 8, e1000382 (2010).
Article PubMed PubMed Central Google Scholar
Dessaud, E. et al. Interpretation of the sonic hedgehog morphogen gradient by a temporal adaptation mechanism. Nature 450, 717–720 (2007).
Article ADS CAS PubMed MATH Google Scholar
Lek, M. et al. A homeodomain feedback circuit underlies step-function interpretation of a Shh morphogen gradient during ventral neural patterning. Development 137, 4051–4060 (2010).
Article CAS PubMed MATH Google Scholar
Huangfu, D. et al. Hedgehog signalling in the mouse requires intraflagellar transport proteins. Nature 426, 83–87 (2003).
Article ADS CAS PubMed MATH Google Scholar
Haycraft, C. et al. Gli2 and Gli3 localize to cilia and require the intraflagellar transport protein polaris for processing and function. PLoS Genet 1, e53 (2005).
Article PubMed PubMed Central MATH Google Scholar
Wichterle, H., Lieberam, I., Porter, J. & Jessell, T. Directed differentiation of embryonic stem cells into motor neurons. Cell 110, 385–397 (2002).
Article CAS PubMed Google Scholar
Duval, N. et al. BMP4 patterns Smad activity and generates stereotyped cell fate organization in spinal organoids. Development 146, dev175430 (2019).
Article CAS PubMed Google Scholar
Holzner, M., Wutz, A. & Di Minin, G. Applying Spinal Cord Organoids as a quantitative approach to study the mammalian Hedgehog pathway. PLoS One 19, e0301670 (2024).
Article CAS PubMed PubMed Central Google Scholar
Maury, Y. et al. Combinatorial analysis of developmental cues efficiently converts human pluripotent stem cells into multiple neuronal subtypes. Nat. Biotechnol. 33, 89–96 (2015).
Article CAS PubMed MATH Google Scholar
Mouilleau, V. et al. Dynamic extrinsic pacing of the HOX clock in human axial progenitors controls motor neuron subtype specification. Development 146, dev194514 (2021).
Article Google Scholar
Ericson, J., Briscoe, J., Rashbass, P., van Heyningen, V. & Jessell, T. Graded sonic hedgehog signaling and the specification of cell fate in the ventral neural tube. Cold Spring Harb. Symp. Quant. Biol. 62, 451–466 (1997).
Article CAS PubMed Google Scholar
Pattyn, A. et al. Coordinated temporal and spatial control of motor neuron and serotonergic neuron generation from a common pool of CNS progenitors. Genes Dev. 17, 729–737 (2003).
Article CAS PubMed PubMed Central MATH Google Scholar
Jang, S., Gumnit, E. & Wichterle, H. A human-specific progenitor sub-domain extends neurogenesis and increases motor neuron production. Nat Neurosci. 27, 1945–1953 (2024).
Rayon, T., Maizels, R., Barrington, C. & Briscoe, J. Single-cell transcriptome profiling of the human developing spinal cord reveals a conserved genetic programme with human-specific features. Development 148, dev199711 (2021).
Article CAS PubMed PubMed Central Google Scholar
Briscoe, J. et al. Homeobox gene Nkx2.2 and specification of neuronal identity by graded Sonic hedgehog signalling. Nature 398, 622–627 (1999).
Article ADS CAS PubMed MATH Google Scholar
Garcia-Gonzalo, F. et al. A transition zone complex regulates mammalian ciliogenesis and ciliary membrane composition. Nat. Genet. 43, 776–784 (2011).
Article CAS PubMed PubMed Central MATH Google Scholar
Williams, C. et al. MKS and NPHP modules cooperate to establish basal body/transition zone membrane associations and ciliary gate function during ciliogenesis. J. Cell Biol. 192, 1023–1041 (2011).
Article CAS PubMed PubMed Central MATH Google Scholar
Jensen, V. et al. Formation of the transition zone by Mks5/Rpgrip1L establishes a ciliary zone of exclusion (CIZE) that compartmentalises ciliary signalling proteins and controls PIP2 ciliary abundance. EMBO J. 34, 2537–2556 (2015).
Article CAS PubMed PubMed Central MATH Google Scholar
Garcia-Gonzalo, F. et al. Phosphoinositides regulate ciliary protein trafficking to modulate hedgehog signaling. Dev. Cell 34, 400–409 (2015).
Article CAS PubMed PubMed Central Google Scholar
Mukhopadhyay, S. et al. The ciliary G-protein-coupled receptor Gpr161 negatively regulates the Sonic hedgehog pathway via cAMP signaling. Cell 152, 210–223 (2013).
Article CAS PubMed MATH Google Scholar
Pal, K. et al. Smoothened determines β-arrestin-mediated removal of the G protein-coupled receptor Gpr161 from the primary cilium. J. Cell Biol. 212, 861–875 (2016).
Article CAS PubMed PubMed Central Google Scholar
Cambray, N. & Wilson, V. Two distinct sources for a population of maturing axial progenitors. development. Development 134, 2829–2840 (2007).
Article CAS PubMed MATH Google Scholar
Cambray, N. & Wilson, V. Axial progenitors with extensive potency are localised to the mouse chordoneural hinge. Development 129, 4855–4866 (2002).
Article CAS PubMed MATH Google Scholar
Dasen, J., Liu, J. & Jessell, T. Motor neuron columnar fate imposed by sequential phases of Hox-c activity. Nature 425, 926–933 (2003).
Article ADS CAS PubMed Google Scholar
Deschamps, J. & Duboule, D. Embryonic timing, axial stem cells, chromatin dynamics, and the Hox clock. Genes Dev. 31, 1406–1416 (2017).
Article CAS PubMed PubMed Central MATH Google Scholar
Miller, A. & Dasen, J. Establishing and maintaining Hox profiles during spinal cord development. Semin Cell Dev. Biol. 152-153, 44–57 (2024).
Article CAS PubMed Google Scholar
Graham, A., Maden, M. & Krumlauf, R. The murine Hox-2 genes display dynamic dorsoventral patterns of expression during central nervous system development. Development 112, 255–264 (1991).
Article CAS PubMed Google Scholar
Pattyn, A., Hirsch, M., Goridis, C. & Brunet, J. Control of hindbrain motor neuron differentiation by the homeobox gene Phox2b. Development 127, 1349–1358 (2000).
Article CAS PubMed Google Scholar
Ahn, D., Ruvinsky, I., Oates, A., Silver, L. & Ho, R. tbx20, a new vertebrate T-box gene expressed in the cranial motor neurons and developing cardiovascular structures in zebrafish. Mech. Dev. 95, 253–258 (2000).
Article CAS PubMed Google Scholar
Guthrie, S. Patterning and axon guidance of cranial motor neurons. Nat. Rev. Neurosci. 8, 859–871 (2007).
Article CAS PubMed MATH Google Scholar
Seitanidou, T., Schneider-Maunoury, S., Desmarquet, C., Wilkinson, D. & Charnay, P. Krox-20 is a key regulator of rhombomere-specific gene expression in the developing hindbrain. Mechanisms Dev. 65, 31–42 (1997).
Article CAS PubMed Google Scholar
Schneider-Maunoury, S., Gilardi-Hebenstreit, P. & Charnay, P. How to build a vertebrate hindbrain. lessons from genetics. C. R. Acad. Sci. III 321, 819–834 (1998).
Article CAS PubMed Google Scholar
Frank, D. & Sela-Donenfeld, D. Hindbrain induction and patterning during early vertebrate development. Cell Mol. Life Sci. 76, 941–960 (2019).
Article CAS PubMed Google Scholar
Aragon, F. & Pujades, C. FGF signaling controls caudal hindbrain specification through Ras-ERK1/2 pathway. BMC Dev. Biol. 9, 61 (2009).
Walshe, J., Maroon, H., McGonnell, I., Dickson, C. & Mason, I. Establishment of hindbrain segmental identity requires signaling by FGF3 and FGF8. Curr. Biol. 12, 1117–1123 (2002).
Article CAS PubMed Google Scholar
Hirsch, M., Glover, J., Dufour, H., Brunet, J. & Goridis, C. Forced expression of Phox2 homeodomain transcription factors induces a branchio-visceromotor axonal phenotype. Dev. Biol. 303, 687–702 (2007).
Article CAS PubMed Google Scholar
Jarrar, W., Dias, J., Ericson, J., Arnold, H. & Holz, A. Nkx2.2 and Nkx2.9 are the key regulators to determine cell fate of branchial and visceral motor neurons in caudal hindbrain. PloS One 10, e0124408 (2015).
Article PubMed PubMed Central Google Scholar
Andre, P., Song, H., Kim, W., Kispert, A. & Yang, Y. Wnt5a and Wnt11 regulate mammalian anterior-posterior axis elongation. Development 142, 1516–1527 (2015).
CAS PubMed PubMed Central Google Scholar
Nordström, U., Maier, E., Jessell, T. & Edlund, T. An early role for WNT signaling in specifying neural patterns of Cdx and Hox gene expression and motor neuron subtype identity. PLoS Biol. 4, e252 (2006).
Article PubMed PubMed Central Google Scholar
van den Akker, E. et al. Cdx1 and Cdx2 have overlapping functions in anteroposterior patterning and posterior axis elongation. Development 129, 2181–2193 (2002).
Article PubMed MATH Google Scholar
Metzis, V. et al. Nervous system regionalization entails axial allocation before neural differentiation. Cell 175, 1105–1118.e1117 (2018).
Article CAS PubMed PubMed Central MATH Google Scholar
Chang, J., Skromne, I. & Ho, R. CDX4 and retinoic acid interact to position the hindbrain-spinal cord transition. Dev. Biol. 410, 178–189 (2016).
Article CAS PubMed PubMed Central Google Scholar
Skromne, I., Thorsen, D., Hale, M., Prince, V. & Ho, R. Repression of the hindbrain developmental program by Cdx factors is required for the specification of the vertebrate spinal cord. Development 134, 2147–2158 (2007).
Article CAS PubMed Google Scholar
Joshi, P., Darr, A. & Skromne, I. CDX4 regulates the progression of neural maturation in the spinal cord. Dev. Biol. 449, 132–142 (2019).
Article CAS PubMed Google Scholar
Wymeersch, F., Wilson, V. & Tsakiridis, A. Understanding axial progenitor biology in vivo and in vitro. Development 148, dev180612 (2021).
Article CAS PubMed Google Scholar
Forlani, S., Lawson, K. & Deschamps, J. Acquisition of Hox codes during gastrulation and axial elongation in the mouse embryo. Development 130, 3807–3819 (2003).
Article CAS PubMed MATH Google Scholar
Gouti, M. et al. A gene regulatory network balances neural and mesoderm specification during vertebrate trunk development. Dev. Cell 41, 243–261.e247 (2017).
Article CAS PubMed PubMed Central Google Scholar
Henrique, D., Abranches, E., Verrier, L. & Storey, K. Neuromesodermal progenitors and the making of the spinal cord. Development 142, 2864–2875 (2015).
Article CAS PubMed Google Scholar
Kondoh, H., Takada, S. & Takemoto, T. Axial level-dependent molecular and cellular mechanisms underlying the genesis of the embryonic neural plate. Dev. Growth Differ. 58, 427–436 (2016).
Article CAS PubMed MATH Google Scholar
Wymeersch, F. et al. Transcriptionally dynamic progenitor populations organised around a stable niche drive axial patterning. Development 146, dev168161 (2019).
PubMed PubMed Central Google Scholar
Huangfu, D. & Anderson, K. Cilia and Hedgehog responsiveness in the mouse. Proc. Natl Acad. Sci. USA 102, 11325–11330 (2005).
Article ADS CAS PubMed PubMed Central MATH Google Scholar
De Mori, R. et al. Joubert syndrome-derived induced pluripotent stem cells show altered neuronal differentiation in vitro. Cell Tissue Res 396, 255–267 (2024).
Article PubMed PubMed Central MATH Google Scholar
Bangs, F., Schrode, N., Hadjantonakis, A.-K. & Anderson, K. Lineage specificity of primary cilia in the mouse embryo. Nat. Cell Biol. 17, 113–122 (2015).
Article CAS PubMed PubMed Central Google Scholar
Soares, H., Carmona, B., Nolasco, S., Viseu Melo, L. & Gonçalves, J. Cilia distal domain: diversity in evolutionarily conserved structures. Cells 8, 160 (2019).
Article CAS PubMed PubMed Central MATH Google Scholar
Wang, J. et al. Variable phenotypes and penetrance between and within different zebrafish ciliary transition zone mutants. Dis. Model Mech. 15, dmm049568 (2022).
Article CAS PubMed PubMed Central Google Scholar
Chang, C. et al. Atoh1 controls primary cilia formation to allow for SHH-triggered granule neuron progenitor proliferation. Developmental Cell 48, 184–199 (2019).
Article CAS PubMed Google Scholar
Ong, T., Trivedi, N., Wakefield, R., Frase, S. & Solecki, D. Siah2 integrates mitogenic and extracellular matrix signals linking neuronal progenitor ciliogenesis with germinal zone occupancy. Nat. Commun. 11, 5312 (2020).
Article ADS CAS PubMed PubMed Central Google Scholar
Blaess, S. et al. Beta1-integrins are critical for cerebellar granule cell precursor proliferation. J. Neurosci. 24, 3402–3412 (2004).
Article CAS PubMed PubMed Central MATH Google Scholar
Akhshi, T. & Trimble, W. A non-canonical Hedgehog pathway initiates ciliogenesis and autophagy. J. Cell Biol. 220, e202004179 (2021).
Article CAS PubMed Google Scholar
Zhang, K. et al. Primary cilia are WNT-transducing organelles whose biogenesis is controlled by a WNT-PP1 axis. 58, 139–154(2023).
Gouti, M. et al. In vitro generation of neuromesodermal progenitors reveals distinct roles for wnt signalling in the specification of spinal cord and paraxial mesoderm identity. PLoS Biol. 12, e1001937 (2014).
Article PubMed PubMed Central Google Scholar
Lippmann, E. et al. Deterministic HOX patterning in human pluripotent stem cell-derived neuroectoderm. Stem Cell Rep. 14, 632–644 (2015).
Article MATH Google Scholar
Afzal, Z. & Krumlauf, R. Transcriptional regulation and implications for controlling Hox gene expression. J. Dev. Biol. 10, 4 (2022).
Article CAS PubMed PubMed Central MATH Google Scholar
Wind, M. et al. Defining the signalling determinants of a posterior ventral spinal cord identity in human neuromesodermal progenitor derivatives. Development 148, dev194415 (2021).
Article CAS PubMed PubMed Central MATH Google Scholar
Cooper, F. et al. Notch signalling influences cell fate decisions and HOX gene induction in axial progenitors. Development 151, dev202098 (2024).
Article CAS PubMed PubMed Central Google Scholar
Young, T. & Deschamps, J. Hox, Cdx, and anteroposterior patterning in the mouse embryo. J. Curr. Top. Dev. Biol. 88, 235–255 (2009).
Article CAS Google Scholar
Takada, S. et al. Wnt-3a regulates somite and tailbud formation in the mouse embryo. Genes Dev. 8, 174–189 (1994).
Article CAS PubMed MATH Google Scholar
Epstein, M., Pillemer, G., Yelin, R., Yisraeli, J. & Fainsod, A. Patterning of the embryo along the anterior-posterior axis: the role of the caudal genes. Development 124, 3805–3814 (1997).
Article CAS PubMed Google Scholar
Davidson, A. & Zon, L. The caudal-related homeobox genes cdx1a and cdx4 act redundantly to regulate hox gene expression and the formation of putative hematopoietic stem cells during zebrafish embryogenesis. Dev. Biol. 292, 506–518 (2006).
Article CAS PubMed Google Scholar
Faas, L. & Isaacs, H. Overlapping functions of Cdx1, Cdx2, and Cdx4 in the development of the amphibian Xenopus tropicalis. Dev. Dyn. 238, 835–852 (2009).
Article CAS PubMed PubMed Central MATH Google Scholar
Amin, S. et al. Cdx and T brachyury co-activate growth signaling in the embryonic axial progenitor niche. Cell Rep. 17, 3165–3177 (2016).
Article CAS PubMed MATH Google Scholar
Neijts, R., Amin, S., van Rooijen, C. & Deschamps, J. Cdx is crucial for the timing mechanism driving colinear Hox activation and defines a trunk segment in the Hox cluster topology. Dev. Biol. 422, 146–154 (2017).
Article CAS PubMed Google Scholar
Mazzoni, E. et al. Saltatory remodeling of Hox chromatin in response to rostrocaudal patterning signals. Nat. Neurosci. 16, 1191–1198 (2013).
Article CAS PubMed PubMed Central MATH Google Scholar
Schyr, R., Shabtai, Y., Shashikant, C. & Fainsod, A. Cdx1 is essential for the initiation of HoxC8 expression during early embryogenesis. FASEB J. 26, 2674–2684 (2012).
Article PubMed Google Scholar
Simons, M. et al. Inversin, the gene product mutated in nephronophthisis type II, functions as a molecular switch between Wnt signaling pathways. Nat. Genet 37, 537–543 (2005).
Article CAS PubMed PubMed Central MATH Google Scholar
Veland, I. et al. Inversin/Nephrocystin-2 is required for fibroblast polarity and directional cell migration. PLoS One 8, e60193 (2013).
Article ADS CAS PubMed PubMed Central MATH Google Scholar
Lancaster, M., Schroth, J. & Gleeson, J. Subcellular spatial regulation of canonical Wnt signalling at the primary cilium. Nat. Cell Biol. 13, 700–707 (2011).
Article PubMed PubMed Central Google Scholar
Ocbina, P., Tuson, M. & Anderson, K. Primary cilia are not required for normal canonical Wnt signaling in the mouse embryo. PLoS One 4, e6839 (2009).
Article ADS PubMed PubMed Central Google Scholar
Vuong, L. & Mlodzik, M. The complex relationship of Wnt-signaling pathways and cilia. Curr. Top. Dev. Biol. 155, 95–125 (2023).
Article PubMed PubMed Central MATH Google Scholar
Bishara, J., Keens, T. & Perez, I. The genetics of congenital central hypoventilation syndrome: clinical implications. Appl Clin. Genet 11, 135–144 (2018).
Article CAS PubMed PubMed Central Google Scholar
Dubreuil, V. et al. A human mutation in Phox2b causes lack of CO₂ chemosensitivity, fatal central apnea, and specific loss of parafacial neurons. Proc. Natl Acad. Sci. USA 105, 1067–1072 (2008).
Article ADS CAS PubMed PubMed Central MATH Google Scholar
Dubreuil, V. et al. Defective respiratory rhythmogenesis and loss of central chemosensitivity in Phox2b mutants targeting retrotrapezoid nucleus neurons. J. Neurosci. 29, 14836–14846 (2009).
Article CAS PubMed PubMed Central MATH Google Scholar
Li, A. & Nattie, E. CO2 dialysis in one chemoreceptor site, the RTN: stimulus intensity and sensitivity in the awake rat. Respir. Physiol. Neurobiol. 133, 11–22 (2002).
Article PubMed Google Scholar
Kanbar, R., Stornetta, R., Cash, D., Lewis, S. & Guyenet, P. Photostimulation of Phox2b medullary neurons activates cardiorespiratory function in conscious rats. Am. J. Respir. Crit. Care Med 182, 1184–1194 (2010).
Article CAS PubMed PubMed Central Google Scholar
Joubert, M., Eisenring, J., Robb, J. & Andermann, F. Familial agenesis of the cerebellar vermis. A syndrome of episodic hyperpnea, abnormal eye movements, ataxia, and retardation. Neurology 19, 813–825 (1969).
Article CAS PubMed Google Scholar
Young, M., Wakefield, M., Smyth, G. & Oshlack, A. Gene ontology analysis for RNA-seq: accounting for selection bias. Genome Biol. 11, R14 (2010).
Article PubMed PubMed Central Google Scholar
Assou, S. et al. Recurrent genetic abnormalities in human pluripotent stem cells: definition and routine detection in culture supernatant by targeted droplet digital PCR. Stem Cell Rep. 14, 1–8 (2020).
Article CAS Google Scholar
Concordet, J. & Haeussler, M. CRISPOR: intuitive guide selection for CRISPR/Cas9 genome editing experiments and screens. Nucleic Acids Res 26, W242–W245 (2018).
Article MATH Google Scholar
Doench, J. et al. Optimized sgRNA design to maximize activity and minimize off-target effects of CRISPR-Cas9. Nat. Biotechnol. 34, 184–191 (2016).
Article CAS PubMed PubMed Central Google Scholar
Haeussler, M. et al. Evaluation of off-target and on-target scoring algorithms and integration into the guide RNA selection tool CRISPOR. Genome Biol. 17, 148 (2016).
Article PubMed PubMed Central MATH Google Scholar
Dobin, A. et al. STAR: ultrafast universal RNA-seq aligner. Bioinformatics 29, 15–21 (2013).
Article CAS PubMed Google Scholar
Ewels, P., Magnusson, M., Lundin, S. & Käller, M. MultiQC: summarize analysis results for multiple tools and samples in a single report. Bioinformatics 32, 3047–3048 (2016).
Article CAS PubMed PubMed Central Google Scholar
Liao, Y., Smyth, G. & Shi, W. featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics 30, 923–930 (2014).
Article CAS PubMed MATH Google Scholar
Lê, S., Josse, J. & Husson, F. FactoMineR: an R package for multivariate analysis. J. Stat. Softw. 25, 1–18 (2008).
Article MATH Google Scholar
Love, M., Huber, W. & Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15, 550 (2014).
Article PubMed PubMed Central MATH Google Scholar
Zink, R., Wolfinger, R. & Mann, G. Summarizing the incidence of adverse events using volcano plots and time intervals. Clin. Trials 10, 398–406 (2013).
Article PubMed MATH Google Scholar
Alhamdoosh, M. et al. Combining multiple tools outperforms individual methods in gene set enrichment analyses. Bioinformatics 33, 414–424 (2017).
Article CAS PubMed Google Scholar

Download references

Acknowledgements

Image acquisition was carried out at the IBPS Imaging Facility. We particularly thank Chloé Chaumeton from the IBPS imaging platform, and Frédéric Eghiaian from Abberior Instruments GmbH, for setting up the STED imaging system. We further thank Jean-François Gilles for his help in setting up the python cellpose environment. We thank the ICV-iPS platform at the Paris Brain Institute (ICM) for the generation and characterization of RPGRIP1L-deficient hiPSC lines, and the iGenSeq sequencing platform at the ICM for reliable sequencing and fast processing of our samples. We are very grateful to Vanessa Ribes and Pascal Gilardi Hebenstreit for providing us with protocols and help in mESCs culture and differentiation approaches. We also thank Jean-François Brunet for kindly providing the Phox2b antibody and the TACGENE team of the Natural History Museum in Paris for providing the Cas9 protein. We thank Estelle Balissat for her help in the preparation of TMEM67 spinal cord organoids and cryosections. We also acknowledge the worldwide contributions of users and developers to the Galaxy project (https://galaxyproject.org/) and all the upstream authors and contributors of the software ecosystem we use and rely on. This work was supported by funding to SSM from the Fondation pour la Recherche Médicale (Equipe FRM EQU201903007943), from the Fondation pour la Recherche sur le Cancer (ARC; PJA 20171206591), from the Agence Nationale de la Recherche under the ERA-NET Neuron 2021 “NDCil” project (ANR-21-NEU2-0009-03), from the Fondation Maladies rares and the Association Mieux vivre avec le syndrome de Joubert under the reference “ASSOJoubert_2024_Schneider”. AW received funding from the German Research Foundation (DFG; WI 5451/1-1). ST and DC received funding from the Agence Nationale de la Recherche (ANR; ANR-17-CE16-0003-01) and from Investissements d’Avenir programs (ANR-10-IAHU-01 and “RHU-C’IL-LICO”/ANR-17-RHUS-0002) operated by the Agence Nationale de la Recherche. VS received funding from the ERA-NET Neuron 2021 “NDCil”. SN received funding from the Agence Nationale de la Recherche (ANR; ANR-19-CE16-0002 and ANR-23-CE16-0026).

Author information

These authors jointly supervised this work: Sylvie Schneider-Maunoury, Aline Stedman.

Authors and Affiliations

Sorbonne Université, CNRS, Inserm, Development, Adaptation and Aging, Dev2A, Paris, France
Antonia Wiegering, Isabelle Anselme, Ludovica Brunetti, Laura Metayer-Derout, Alexis Eschstruth, Martin Catala, Sylvie Schneider-Maunoury & Aline Stedman
Sorbonne Université, CNRS, Inserm, Institut de Biologie Paris Seine, Paris, France
Antonia Wiegering, Isabelle Anselme, Ludovica Brunetti, Laura Metayer-Derout, Alexis Eschstruth, Martin Catala, Christophe Antoniewski, Sylvie Schneider-Maunoury & Aline Stedman
INSERM UMR 1163, Institut Imagine, Université Paris Cité, Paris, France
Damelys Calderon & Sophie Thomas
Sorbonne Université, Inserm, Institut du Fer à Moulin, UMR-S 1270, Paris, France
Stéphane Nedelec
Université Paris Cité, CNRS, Inserm U1340, Institut Jacques Monod, Paris, France
Stéphane Nedelec
Department of Molecular Medicine, University of Pavia, Pavia, Italy
Valentina Serpieri

Authors

Antonia Wiegering
View author publications
You can also search for this author inPubMed Google Scholar
Isabelle Anselme
View author publications
You can also search for this author inPubMed Google Scholar
Ludovica Brunetti
View author publications
You can also search for this author inPubMed Google Scholar
Laura Metayer-Derout
View author publications
You can also search for this author inPubMed Google Scholar
Damelys Calderon
View author publications
You can also search for this author inPubMed Google Scholar
Sophie Thomas
View author publications
You can also search for this author inPubMed Google Scholar
Stéphane Nedelec
View author publications
You can also search for this author inPubMed Google Scholar
Alexis Eschstruth
View author publications
You can also search for this author inPubMed Google Scholar
Valentina Serpieri
View author publications
You can also search for this author inPubMed Google Scholar
Martin Catala
View author publications
You can also search for this author inPubMed Google Scholar
Christophe Antoniewski
View author publications
You can also search for this author inPubMed Google Scholar
Sylvie Schneider-Maunoury
View author publications
You can also search for this author inPubMed Google Scholar
Aline Stedman
View author publications
You can also search for this author inPubMed Google Scholar

Contributions

A.S., S.S.M. and A.W. designed the study, interpreted the data and wrote the manuscript. A.W. & A.S. supervised and conducted the experiments. I.A., L. M.-D. & L.B. conducted experiments for the production and analysis of mouse and human spinal organoids. D.C., S.T., A.E. & V.S. designed and conducted experiments to produce iPSC mutant models. C.A. designed and conducted RNAseq analysis. S.N. & M.C. participated in designing the study, interpreting the data and writing the manuscript.

Corresponding authors

Correspondence to Antonia Wiegering, Sylvie Schneider-Maunoury or Aline Stedman.

Ethics declarations

Competing interests

The authors declare no competing interests.

Peer review

Peer review information

Nature Communications thanks Natalia Rodriguez-Muela, and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. A peer review file is available.

Additional information

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

Supplementary information

Supplementary Information

Reporting Summary

Transparent Peer Review file

Source data

Source Data

Rights and permissions

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/.

Reprints and permissions

About this article

Cite this article

Wiegering, A., Anselme, I., Brunetti, L. et al. A differential requirement for ciliary transition zone proteins in human and mouse neural progenitor fate specification. Nat Commun 16, 3258 (2025). https://doi.org/10.1038/s41467-025-58554-3

Download citation

Received: 24 July 2024
Accepted: 24 March 2025
Published: 05 April 2025
DOI: https://doi.org/10.1038/s41467-025-58554-3