# RNA Profiling of the Human and Mouse Spinal Cord Stem Cell Niches Reveals an Embryonic-like Regionalization with MSX1+ Roof-Plate-Derived Cells

Hussein Ghazale, Chantal Ripoll, Nicolas Leventoux, Laurent Jacob, Safa Azar, Daria Mamaeva, Yael Glasson, Charles-Félix Calvo, Jean-Léon Thomas, Sarah Meneceur, et al.

► To cite this version:

Hussein Ghazale, Chantal Ripoll, Nicolas Leventoux, Laurent Jacob, Safa Azar, et al. RNA Profiling of the Human and Mouse Spinal Cord Stem Cell Niches Reveals an Embryonic-like Regionalization with MSX1+ Roof-Plate-Derived Cells. Stem Cell Reports, 2019, 12 (5), pp.1159-1177. 10.1016/j.stemcr.2019.04.001. hal-02122271

HAL Id: hal-02122271

https://hal.sorbonne-universite.fr/hal-02122271

Submitted on 7 May 2019

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

# RNA Profiling of the Human and Mouse Spinal Cord Stem Cell Niches Reveals an Embryonic-like Regionalization with MSX1+ Roof-Plate-Derived Cells

Hussein Ghazale,<fnref n="1" />,<fnref n="9" /> Chantal Ripoll,<fnref n="1" />,<fnref n="9" /> Nicolas Leventoux,<fnref n="1" /> Laurent Jacob,<fnref n="2" /> Safa Azar,<fnref n="1" /> Daria Mamaeva,<fnref n="1" /> Yael Glasson,<fnref n="1" /> Charles-Felix Calvo,<fnref n="3" /> Jean-Leon Thomas,<fnref n="2" />,<fnref n="4" /> Sarah Menecœur,<fnref n="5" /> Yvan Lallemand,<fnref n="5" /> Valérie Rigau,<fnref n="1" />,<fnref n="6" /> Florence E. Perrin,<fnref n="8" /> Harun N. Noristani,<fnref n="8" /> Brenda Rocamonde,<fnref n="2" /> Emmanuelle Huillard,<fnref n="2" /> Luc Bauchet,<fnref n="1" />,<fnref n="7" /> and Jean-Philippe Hugnot<fnref n="1" />,<fnref n="8" />,*

## SUMMARY

Anamniotes, rodents, and young humans maintain neural stem cells in the ependymal zone (EZ) around the central canal of the spinal cord, representing a possible endogenous source for repair in mammalian lesions. Cell diversity and genes specific for this region are ill defined. A cellular and molecular resource is provided here for the mouse and human EZ based on RNA profiling, immunostaining, and fluorescent transgenic mice. This uncovered the conserved expression of 1,200 genes including 120 transcription factors. Unexpectedly the EZ maintains an embryonic-like dorsal-ventral pattern of expression of spinal cord developmental transcription factors (ARX, FOXA2, MSX1, and PAX6). In mice, dorsal and ventral EZ cells express *Vegfr3* and are derived from the embryonic roof and floor plates. The dorsal EZ expresses a high level of *Bmp6* and *Gdf10* genes and harbors a subpopulation of radial quiescent cells expressing MSX1 and ID4 tran- scription factors.

## INTRODUCTION

The adult central nervous system maintains neural stem cells in specific areas called niches (Gage and Temple, 2013). The main stem cell pools are in the subventricular zone (SVZ) and in the hippocampus. A third stem cell niche is found in the adult spinal cord around the central canal in anamniotes, rodents, and humans (Becker et al., 2018; Marichal et al., 2017; Stenudd et al., 2015). This niche orig- inates from the embryonic neuroepithelium and forms the spinal cord ependymal zone (EZ) organized as a pseudo- epithelium. As in the brain, this stem cell niche is highly organized and contains stem and non-stem cells. Four different cell types, namely ependymocytes, cerebrospi- nal-fluid-contacting neurons (CSF-N), vessels, and long radial cells have been described in mice. In particular, the EZ presents a dorsal-ventral regionalization with long radial glial fibrillary acidic protein (GFAP)⁺ cells in the dor- sal part (Sabourin et al., 2009). In contrast to brain, spinal cord ependymal cells slowly proliferate to self-renew (Al- faro-Cervello et al., 2012; Pfenninger et al., 2011). In cul- tures, a fraction of these cells can also generate passageable neurospheres (i.e. clonal expansion of neural precursor

cells) which can generate astrocytes, oligodendrocytes, and neurons after differentiation (Weiss et al., 1996). The identity of these neurosphere-forming cells in the EZ is still not completely clear as both GFAP⁺ and GFAP⁻ ependymal cells can behave as neural stem cells *in vitro* (Barnabé- Heider et al., 2010; Fiorelli et al., 2013; Sabourin et al., 2009; Xu et al., 2017). Recent single cell analysis has iden- tified neurogenesis in the adult spinal cord (Habib et al., 2016); however, whether these new neurons are derived from the EZ is not yet established. It has been known since 1962 (Adrian and Walker, 1962) that the EZ can readily acti- vate and produce new cells upon injury (Becker et al., 2018). Depending on the lesion type and severity, EZ- derived cells can significantly contribute to the glial scar formation (Ren et al., 2017; Stenudd et al., 2015).

In comparison with the brain niches, less is known about the adult spinal cord EZ. Reminiscent of the mouse niche, in human, ependymal cells around the central canal display immature features such as expression of NES (nes- tin), VIM (vimentin), and SOX2 (Becker et al., 2018). How- ever, with aging the central lumen can disappear and the EZ is disorganized (Garcia-Ovejero et al., 2015). Multipo- tent neurospheres with a limited proliferation ability

<footnote n="1">¹INSERM U1051, INM, Hopital Saint Eloi, 80 Avenue Augustin Fliche, 34091 Montpellier, France</footnote>

<footnote n="2">²Inserm U 1127, CNRS UMR 7225, Sorbonne Universités, UPMC Univ Paris 06 UMR S 1127, Institut du Cerveau et de la Moelle épinière, ICM, 75013, Paris, France</footnote>

<footnote n="3">³Neuroglial Interactions in Cerebral Physiopathology. CIRB, CNRS UMR 7241/INSERM U1050 Collège de France 11, Place Marcelin Berthelot, 75005 Paris, France</footnote>

<footnote n="4">⁴Department of Neurology, Yale University School of Medicine, New Haven, CT 06511, USA</footnote>

<footnote n="5">⁵Institut Pasteur, Department of Developmental and Stem Cell Biology, UMR3738 CNRS, Paris 75015, France</footnote>

<footnote n="6">⁶CHU of Montpellier, Hopital Gui de Chaulliac, Pathology Department, 80 Avenue Augustin Fliche, 34091 Montpellier, France</footnote>

<footnote n="7">⁷CHU of Montpellier, Hopital Gui de Chaulliac, Neurosurgery Department, 80 Avenue Augustin Fliche, 34091 Montpellier, France</footnote>

<footnote n="8">⁸University of Montpellier, Faculté des Sciences, Place Eugène Bataillon, 34095 Montpellier Cedex 5, France</footnote>

<footnote n="9">⁹Co-first author</footnote>

<footnote>*Correspondence: jean-philippe.hugnot@umontpellier.fr</footnote>

<footnote>https://doi.org/10.1016/j.stemcr.2019.04.001</footnote>
<figure><img src="image_3.png" /></figure>

A

Hoechst

<figure><img src="image_4.png" /></figure>

mEZ

<figure><img src="image_5.png" /></figure>

hEZ1 (17y.)

<figure><img src="image_6.png" /></figure>

hEZ2 (46y.)

B

mouse

<figure><img src="image_7.png" /></figure>

human

<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /></figure>

C

mouse

<figure><img src="image_10.png" /></figure>

human

<figure><img src="image_11.png" /></figure>

D

<figure><img src="image_12.png" /></figure>

mouse upregulated EZ genes (fold change ≥ 2)

human upregulated EZ genes (fold change ≥ 2)

7510

1223

899

E

Transcription factors

(p< 2.5 10⁻⁷)

<table><tr><td>AHCTF1</td><td>FOXJ1</td><td>KLF5</td><td>NFATC2</td><td>PIAS1</td><td>SMAD5</td><td>ZFP3</td></tr><tr><td>ARX</td><td>GABPA</td><td>MDFIC</td><td>NFIA</td><td>POU2F1</td><td>SOX11</td><td>ZFPM2</td></tr><tr><td>ATF1</td><td>GLIS3</td><td>MED1</td><td>NFIB</td><td>POU3F2</td><td>SOX2</td><td>ZFX</td></tr><tr><td>BACH2</td><td>GPBP1</td><td>MED13L</td><td>NFIX</td><td>PRDM5</td><td>SOX4</td><td>ZHX2</td></tr><tr><td>BCL10</td><td>HDGF</td><td>MED14</td><td>NFXL1</td><td>PSIP1</td><td>SOX5</td><td>ZKSCAN7</td></tr><tr><td>BCOR</td><td>HES1</td><td>MED17</td><td>NFYB</td><td>PTPN14</td><td>SOX6</td><td>ZKSCAN8</td></tr><tr><td>BPTF</td><td>HMGB2</td><td>MED25</td><td>NOTCH1</td><td>PURB</td><td>SOX9</td><td>ZSCAN18</td></tr><tr><td>CENPF</td><td>HOXA2</td><td>MED30</td><td>NR1D2</td><td>REST</td><td>SP4</td><td></td></tr><tr><td>CITED2</td><td>HOXA3</td><td>MED31</td><td>NR2C1</td><td>RFX2</td><td>SPEN</td><td></td></tr><tr><td>CREBBP</td><td>HOXB3</td><td>MEIS2</td><td>NR2C2</td><td>RFX3</td><td>TAF1</td><td></td></tr><tr><td>CREBRF</td><td>ID4</td><td>MINA</td><td>NR2F1</td><td>RFX4</td><td>TAF1B</td><td></td></tr><tr><td>CREBZF</td><td>IKZF2</td><td>MLXIP</td><td>NR2F2</td><td>RORA</td><td>TBL1X</td><td></td></tr><tr><td>CRY1</td><td>IKZF4</td><td>MSX1</td><td>NR4A3</td><td>RREB1</td><td>TCF7L1</td><td></td></tr><tr><td>CTNNB1</td><td>IRF2</td><td>MTA1</td><td>NRF1</td><td>RUNX1</td><td>TCF7L2</td><td></td></tr><tr><td>DDX20</td><td>IRF6</td><td>MTA3</td><td>PAX6</td><td>SALL1</td><td>TEAD1</td><td></td></tr><tr><td>E2F5</td><td>KAT6A</td><td>MTF1</td><td>PAXBP1</td><td>SALL3</td><td>TRPS1</td><td></td></tr><tr><td>E2F7</td><td>KCTD1</td><td>MYB</td><td>PBX1</td><td>SATB1</td><td>TRRAP</td><td></td></tr><tr><td>EP300</td><td>KDM3A</td><td>NCOA3</td><td>PER2</td><td>SIN3A</td><td>UBN2</td><td></td></tr><tr><td>FOXA1</td><td>KDM5C</td><td>NFATC1</td><td>PHF21A</td><td>SMAD2</td><td>YAP1</td><td></td></tr></table>

Cilia

(p< 4.6 10⁻⁵²)

<table><tr><td>ARL13B</td><td>CEP164</td><td>IFT52</td><td>TTC8</td></tr><tr><td>B9D1</td><td>CEP290</td><td>IFT57</td><td>TMEM107</td></tr><tr><td>BBIP1</td><td>CEP41</td><td>IFT74</td><td>TMEM138</td></tr><tr><td>BBS1</td><td>CEP83</td><td>IFT80</td><td>TMEM17</td></tr><tr><td>BBS2</td><td>CFAP221</td><td>IFT81</td><td>TMEM216</td></tr><tr><td>BBS4</td><td>CFAP54</td><td>IFT88</td><td>TMEM231</td></tr><tr><td>BBS5</td><td>CC2D2A</td><td>INTU</td><td>TMEM237</td></tr><tr><td>DZIPI1</td><td>CCDC113</td><td>NOTCH1</td><td>UNC119B</td></tr><tr><td>DZIPI</td><td>DNAAF1</td><td>PCM1</td><td></td></tr><tr><td>FBF1</td><td>DNAH5</td><td>RFX2</td><td></td></tr><tr><td>IQUB</td><td>DYNC2H1</td><td>RFX3</td><td></td></tr><tr><td>NEK1</td><td>FAM161A</td><td>RPGR</td><td></td></tr><tr><td>OFD1</td><td>FOXJ1</td><td>SCLT1</td><td></td></tr><tr><td>RILPL2</td><td>HSPB11</td><td>SPAG16</td><td></td></tr><tr><td>TBC1D32</td><td>IFT122</td><td>CTTN1</td><td></td></tr><tr><td>WDR19</td><td>IFT172</td><td>CTTN3</td><td></td></tr><tr><td>WDR35</td><td>IFT20</td><td>TEKT3</td><td></td></tr><tr><td>CEP126</td><td>IFT43</td><td>TTC21A</td><td></td></tr><tr><td>CEP162</td><td>IFT46</td><td>TTC26</td><td></td></tr></table>

Division

(p< 1.4 10⁻³)

<table><tr><td>CDC14A</td><td>NEK2</td></tr><tr><td>CDK13</td><td>NEK4</td></tr><tr><td>CDK6</td><td>NUDC</td></tr><tr><td>CDK7</td><td>NUMA1</td></tr><tr><td>CENPF</td><td>NUP37</td></tr><tr><td>CEP164</td><td>PARD3B</td></tr><tr><td>CETN2</td><td>PARD6B</td></tr><tr><td>CETN3</td><td>PARD6G</td></tr><tr><td>CHEK2</td><td>PSRC1</td></tr><tr><td>CIB1</td><td>RUVBL1</td></tr><tr><td>CKS1B</td><td>SMC1A</td></tr><tr><td>CNTRL</td><td>SMC3</td></tr><tr><td>EVI5</td><td>SPDL1</td></tr><tr><td>KATNA1</td><td>SPICE1</td></tr><tr><td>KATNB1</td><td>STAG2</td></tr><tr><td>LATS1</td><td>STOX1</td></tr><tr><td>MAD2L1</td><td>TTC28</td></tr><tr><td>NEDD1</td><td></td></tr><tr><td>NEK1</td><td></td></tr></table>

Smo.

(p<5.5 10⁻¹⁰)

<table><tr><td>ARL13B</td><td>B9D1</td></tr><tr><td>DZIPI</td><td>IQUB</td></tr><tr><td>CLUAP1</td><td>CC2D2A</td></tr><tr><td>DYRK2</td><td>HSPB11</td></tr><tr><td>HE51</td><td>HIPK1</td></tr><tr><td>IFT172</td><td>IFT20</td></tr><tr><td>IFT27</td><td>IFT46</td></tr><tr><td>IFT52</td><td>IFT57</td></tr><tr><td>IFT80</td><td>PAX6</td></tr><tr><td>TCTN3</td><td></td></tr></table>

F

Dorsal

<figure><img src="image_13.png" /></figure>

Lateral

Ventral

G

VENTRAL EZ ENRICHED GENES

fold change ≥ 4.0

<table><thead><tr><th>Gene Symbol</th><th>Gene description</th><th>Fold change Ventral vs Lateral</th><th>pval</th></tr></thead><tbody><tr><td>Sulf1</td><td>sulfatase 1</td><td>14,4</td><td>0,03</td></tr><tr><td>Ucma</td><td>associated</td><td>11,7</td><td>0,03</td></tr><tr><td>Arx</td><td>aristaless related homeobox</td><td>11,5</td><td>0,03</td></tr><tr><td>Dcn</td><td>decorin</td><td>9,2</td><td>0,03</td></tr><tr><td>Rbp4</td><td>retinol binding protein 4, plasma</td><td>8,0</td><td>0,03</td></tr><tr><td>Sulf1</td><td>sulfatase 1</td><td>6,1</td><td>0,03</td></tr><tr><td>Foxa1</td><td>forkhead box A1</td><td>4,7</td><td>0,03</td></tr><tr><td>Prss23</td><td>serine protease 23</td><td>4,7</td><td>0,03</td></tr><tr><td>Fgl2</td><td>fibrinogen-like protein 2</td><td>4,3</td><td>0,03</td></tr></tbody></table>

DORSAL EZ ENRICHED GENES

fold change ≥ 2.0

<table><thead><tr><th>Gene Symbol</th><th>Gene description</th><th>Fold change Dorsal vs Lateral</th><th>pval</th></tr></thead><tbody><tr><td>A2m</td><td>alpha-2-macroglobulin</td><td>5,8</td><td>0,03</td></tr><tr><td>Gdf10</td><td>growth differentiation factor 10</td><td>5,3</td><td>0,03</td></tr><tr><td>Gcnt1</td><td>glucosaminyl (N-acetyl) transferase 1</td><td>3,3</td><td>0,03</td></tr><tr><td>Bmp6</td><td>bone morphogenetic protein 6</td><td>3,0</td><td>0,03</td></tr><tr><td>Efnb3</td><td>epsilon B3</td><td>2,6</td><td>0,03</td></tr><tr><td>Prokr2</td><td>prokineticin receptor 2</td><td>2,5</td><td>0,03</td></tr><tr><td>Dclk1</td><td>doublecortin-like kinase 1</td><td>2,2</td><td>0,03</td></tr><tr><td>Thbs2</td><td>thrombospondin 2</td><td>2,0</td><td>0,03</td></tr><tr><td>Hexb</td><td>hexosaminidase B</td><td>2,0</td><td>0,03</td></tr></tbody></table>

(legend on next page)
<figure><img src="image_14.png" /></figure>

have been derived from the human spinal cord (Dromard et al., 2008) and using alternative culture conditions, Mothe et al. (2011) were able to maintain a sustained proliferation of multipotent human-derived neural stem cells.

A detailed transcriptomic profiling of the human and mouse EZ is currently lacking. This would help us to understand the specificity and diversity of these cells as well as identify gene expressions and molecular pathways conserved between primates and rodents. It would also provide important insights into why, in contrast to anamniotes, mammalian ependymal cells cannot regenerate neurons after spinal cord injury (Becker et al., 2018).

Here we provide a cellular and molecular resource for the mouse and human EZ based on RNA profiling, immunostaining, and fluorescent transgenic mice. This uncovered the conserved expression of 1,200 genes specifically expressed in the EZ, including 120 transcription factors (TFs). Unexpectedly, the EZ maintains an embryonic-like dorsal-ventral pattern of expression of spinal cord developmental TFs. New subpopulations of cells expressing specific genes were identified in the dorsal and ventral part of the EZ. In mice, dorsal ependymal cells were found to be derived from the embryonic spinal cord roof plate.

## RESULTS

### Genes Enriched in the Adult Human and Mouse EZ

An epithelial organization of the EZ is observed both in human and mouse as evidenced by CTNNB1 (β-catenin) and CD24 stainings (Figures 2 and 3). To identify gene expression enriched in the EZ, we microdissected this region and adjacent tissue in two human samples and four mice (Figure 1A). For human, we selected two samples with a lumen from patients aged 17 and 46 years. Microarrays were used for RNA profiling and heatmaps indicated adequate clustering of EZ samples (Figure 1B). Volcano plots showed 8,733 and 2,122 genes enriched (fold change ≥2) in the mouse and human EZ, respectively, and 1,223 genes commonly enriched in both species (Figures 1C and 1D). Table 1 shows the top 15 genes commonly

enriched in the human and mouse EZ and genes that are more specifically enriched in human or mice. Figure S1A shows examples of identified genes whose specific expression in the mouse EZ is confirmed in the Allen brain atlas (Lein et al., 2007). In mouse and human, GO (Gene Ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway analyses revealed an enrichment for genes involved in cilia formation, smoothened and hippo pathways, cell division, and transcription (Figure 1E and Tables S1, S2, and S3). In addition, consistent with their position at the interface between CSF and the nervous parenchyma, ependymal cells also expressed 129 and 34 members of the solute carrier family in mice and human, respectively (Tables S1, S2, and S3). Some of them are very specifically (fold change >30) expressed in the mouse EZ such as *Slc26A3*, *Slc14a1*, and *Slc16a12*, which transport chloride, urea, and monocarboxylic acid, respectively. Unexpectedly, *Cftr*, a chloride transporter responsible for cystic fibrosis, was found to be specifically expressed in the EZ in both species (Table 1). Using a single cell transcriptome approach in the mouse spinal cord, two recent studies (Rosenberg et al., 2018; Zeisel et al., 2018) provided a limited number of genes (<50) enriched in spinal cord ependymal cells and CSF-N, and most of them were identified in our study (Table S4). Interestingly, only *Tmem212*, a gene coding for a transmembrane protein with few annotations, was identified in the three studies (including ours) as being enriched in the ependymal cells. With regard to CSF-N, *Espn* (*Espin*) and *Pkd2l1* were identified in the three studies (Table S4).

We and others have reported that the dorsal and ventral parts of the mouse EZ have distinctive features such as the presence of radial cells expressing NES protein (Alfaro-Cervello et al., 2012; Becker et al., 2018; Hamilton et al., 2009) and the preferential expression of the TF ZEB1 in the dorsal part (Sabourin et al., 2009). Further exploration of regional gene expression was done by microdissecting the dorsal, lateral, and ventral parts of the mouse EZ (Figure 1F, n = 4 mice). Only few regionally-expressed genes were identified (Figure 1G and Table S5). In the ventral part, we found a strong expression of three genes

### Figure 1. RNA Profiles of the Mouse and Human EZ

(A) Aspect of EZ in the mouse (mEZ, thoraco level) and human samples (hEZ1 and hEZ2, thoracolumbar level). A lumen was present in mice and in the two human (aged 17 and 46 years) samples. Microdissected EZ and peri EZ regions are delimited with dotted circles. Scale bars, 100 µm.

(B) Heatmap of hierarchical clustering of genes expressed in EZ and peri EZ (pEZ) regions in the four mouse and two human samples.

(C) Volcano plots of genes whose expression is enriched in the mouse and human EZ (fold change ≥2).

(D) Venn diagram of genes enriched in the mouse and human EZ.

(E) Genes enriched for ciliogenesis, smoothened pathway (Smo.), division, and transcription factors in the mouse and human EZ. Expression of genes in red was subsequently confirmed at the protein level (Figures 2 and 3).

(F) Microdissected subregions of the mouse EZ.

(G) Genes enriched (top 9) in the ventral and dorsal EZ (full lists are in Table S5).
<figure><img src="image_15.png" /></figure>

Please cite this article in press as: Ghazale et al., RNA Profiling of the Human and Mouse Spinal Cord Stem Cell Niches Reveals an Embryonic-like Regionalization with MSX1+ Roof-Plate-Derived Cells, Stem Cell Reports (2019), https://doi.org/10.1016/j.stemcr.2019.04.001

<table><caption>Table 1. Three Lists of Genes (Top 15) Enriched in Mouse, Human, and Mouse and Human EZ</caption><thead><tr><th colspan="6">Top 15 Human EZ-Specific Genes</th><th colspan="6">Top 15 Mouse EZ-Specific Genes</th></tr><tr><th>Gene</th><th>Description</th><th>hEZ (log<sub>2</sub>)</th><th>Peri hEZ (log<sub>2</sub>)</th><th>Fold Change</th><th>p Value</th><th>Gene</th><th>Description</th><th>mEZ (log<sub>2</sub>)</th><th>Peri mEZ (log<sub>2</sub>)</th><th>Fold Change</th><th>p Value</th></tr></thead><tbody><tr><td>C7orf57</td><td>chromosome 7 open reading frame 57</td><td>10.9</td><td>4.9</td><td>61.8</td><td>0.0013</td><td>Rsph4a</td><td>radial spoke head 4 homolog A</td><td>10.9</td><td>3.8</td><td>133.1</td><td>2.32 x 10<sup>-13</sup></td></tr><tr><td>VWA3B</td><td>von Willebrand factor A domain-containing 3B</td><td>11.4</td><td>5.7</td><td>52.8</td><td>7.43 x 10<sup>-5</sup></td><td>1500015010Rik</td><td>RIKEN cDNA 1500015010 gene</td><td>11.1</td><td>4.0</td><td>132.7</td><td>1.38 x 10<sup>-12</sup></td></tr><tr><td>NEK5</td><td>NIMA-related kinase 5</td><td>11.1</td><td>5.5</td><td>47.9</td><td>6.43 x 10<sup>-5</sup></td><td>Stoml3</td><td>stomatin (Epb7.2)-like 3</td><td>10.4</td><td>3.5</td><td>121.7</td><td>8.47 x 10<sup>-13</sup></td></tr><tr><td>CD36</td><td>CD36 molecule</td><td>10.5</td><td>5.0</td><td>45.8</td><td>3.50 x 10<sup>-5</sup></td><td>1700007K13Rik</td><td>RIKEN cDNA 1700007K13 gene</td><td>11.1</td><td>4.2</td><td>120.3</td><td>2.66 x 10<sup>-11</sup></td></tr><tr><td>EFHB</td><td>EF-hand domain family, member B</td><td>11.2</td><td>5.7</td><td>44.0</td><td>0.0004</td><td>Slc26a3</td><td>solute carrier family 26, member 3</td><td>11.2</td><td>4.4</td><td>110.5</td><td>2.05 x 10<sup>-12</sup></td></tr><tr><td>MYLK3</td><td>myosin light-chain kinase 3</td><td>11.4</td><td>6.0</td><td>42.6</td><td>0.0044</td><td>Tnnc2</td><td>troponin C2, fast</td><td>10.1</td><td>3.4</td><td>105.2</td><td>8.39 x 10<sup>-5</sup></td></tr><tr><td>CCDC39</td><td>coiled-coil domain-containing 39</td><td>10.0</td><td>4.6</td><td>42.4</td><td>0.0001</td><td>2810047C21Rik1</td><td>RIKEN cDNA 2810047C21</td><td>9.8</td><td>3.2</td><td>103.2</td><td>4.85 x 10<sup>-11</sup></td></tr><tr><td>ODF3B</td><td>outer dense fiber of sperm tails 3B</td><td>9.9</td><td>4.5</td><td>41.5</td><td>0.0008</td><td>Chil3; Chil4</td><td>chitinase-like 3; chitinase-like 4</td><td>11.5</td><td>4.8</td><td>102.2</td><td>4.87 x 10<sup>-5</sup></td></tr><tr><td>ZBBX</td><td>zinc finger, B-box domain-containing</td><td>10.6</td><td>5.3</td><td>41.3</td><td>2.38 x 10<sup>-5</sup></td><td>C1qtnf3</td><td>C1q and tumor necrosis factor related protein 3</td><td>10.7</td><td>4.1</td><td>97.9</td><td>1.63 x 10<sup>-11</sup></td></tr><tr><td>CCDC114</td><td>coiled-coil domain-containing 114</td><td>11.3</td><td>6.0</td><td>40.3</td><td>0.0005</td><td>Cfap161</td><td>cilia- and flagella-associated protein 161</td><td>10.5</td><td>3.9</td><td>95.3</td><td>5.65 x 10<sup>-13</sup></td></tr><tr><td>FAM216B</td><td>family with sequence similarity 216, member B</td><td>9.9</td><td>4.5</td><td>39.8</td><td>0.0012</td><td>Fam183b</td><td>family with sequence similarity 183, member B</td><td>11.5</td><td>4.9</td><td>91.8</td><td>1.06 x 10<sup>-11</sup></td></tr><tr><td>CFAP43</td><td>cilia- and flagella-associated protein 43</td><td>11.3</td><td>6.0</td><td>39.1</td><td>0.0003</td><td>Gm11992</td><td>predicted gene 11.992</td><td>9.8</td><td>3.3</td><td>91.6</td><td>7.02 x 10<sup>-13</sup></td></tr><tr><td>FM03</td><td>flavin-containing monooxygenase 3</td><td>9.9</td><td>4.6</td><td>38.9</td><td>0.0014</td><td>Iqca</td><td>IQ motif containing with AAA domain</td><td>11.0</td><td>4.5</td><td>90.4</td><td>1.12 x 10<sup>-10</sup></td></tr></tbody></table>

(Continued on next page)
<table><caption>**Table 1. Continued**</caption><caption>Top 15 Human EZ-Specific Genes</caption><thead><tr><th>Gene</th><th>Description</th><th>hEZ (log<sub>2</sub>)</th><th>Peri hEZ (log<sub>2</sub>)</th><th>Fold Change</th><th>p Value</th></tr></thead><tbody><tr><td>CFAP70</td><td>cilia- and flagella-associated protein 70</td><td>11.3</td><td>6.0</td><td>38.5</td><td>0.0023</td></tr><tr><td>SPATA17</td><td>spermatogenesis-associated 17</td><td>10.3</td><td>5.1</td><td>37.1</td><td>0.0003</td></tr></tbody></table><table><caption>Top 15 Mouse EZ-Specific Genes</caption><thead><tr><th>Gene</th><th>Description</th><th>mEZ (log<sub>2</sub>)</th><th>Peri mEZ (log<sub>2</sub>)</th><th>Fold Change</th><th>p Value</th></tr></thead><tbody><tr><td>Capsl</td><td>calcyphosin-like</td><td>9.6</td><td>3.1</td><td>87.3</td><td>5.95 x 10<sup>-12</sup></td></tr><tr><td>Acta1</td><td>actin, alpha 1, skeletal muscle</td><td>10.7</td><td>4.2</td><td>85.5</td><td>0.0002</td></tr></tbody></table>

<table><caption>**TOP15 Human and mouse EZ-specific genes**</caption><thead><tr><th>Gene</th><th>Description</th><th>hEZ (log<sub>2</sub>)</th><th>peri hEZ (log<sub>2</sub>)</th><th>Fold Change</th><th>p Value</th><th>mEZ (log<sub>2</sub>)</th><th>peri mEZ (log<sub>2</sub>)</th><th>Fold Change</th><th>p Value</th></tr></thead><tbody><tr><td>CCDC39</td><td>coiled-coil domain containing 39</td><td>10.0</td><td>4.6</td><td>42.4</td><td>0.0001</td><td>9.6</td><td>3.5</td><td>68.5</td><td>9.13 x 10<sup>-12</sup></td></tr><tr><td>ODF3B</td><td>outer dense fiber of sperm tails 3B</td><td>9.9</td><td>4.5</td><td>41.5</td><td>0.0008</td><td>10.5</td><td>4.4</td><td>68.7</td><td>1.01 x 10<sup>-10</sup></td></tr><tr><td>RSPH4A</td><td>radial spoke head 4 homolog A (Chlamydomonas)</td><td>9.8</td><td>4.9</td><td>30.0</td><td>0.0004</td><td>10.9</td><td>3.8</td><td>133.1</td><td>2.32 x 10<sup>-13</sup></td></tr><tr><td>ARMC4</td><td>armadillo repeat containing 4</td><td>10.3</td><td>5.1</td><td>36.7</td><td>3.27 x 10<sup>-5</sup></td><td>9.9</td><td>4.1</td><td>58.9</td><td>9.43 x 10<sup>-12</sup></td></tr><tr><td>CAPSL</td><td>calcyphosine-like</td><td>11.3</td><td>6.4</td><td>31.0</td><td>0.0001</td><td>9.6</td><td>3.1</td><td>87.3</td><td>5.95 x 10<sup>-12</sup></td></tr><tr><td>MNS1</td><td>meiosis-specific nuclear structural 1</td><td>11.0</td><td>6.1</td><td>30.7</td><td>0.0014</td><td>10.4</td><td>4.2</td><td>72.3</td><td>3.96 x 10<sup>-13</sup></td></tr><tr><td>DYNLRB2</td><td>dynein, light chain, roadblock-type 2</td><td>11.2</td><td>6.2</td><td>30.6</td><td>0.0007</td><td>11.3</td><td>5.1</td><td>72.4</td><td>1.63 x 10<sup>-12</sup></td></tr><tr><td>TTC29</td><td>tetratricopeptide repeat domain 29</td><td>9.7</td><td>4.9</td><td>28.4</td><td>6.06 x 10<sup>-5</sup></td><td>9.5</td><td>3.4</td><td>67.9</td><td>8.98 x 10<sup>-13</sup></td></tr><tr><td>TEKT1</td><td>tektin 1</td><td>9.8</td><td>5.0</td><td>27.3</td><td>0.0004</td><td>10.1</td><td>4.1</td><td>64.7</td><td>2.59 x 10<sup>-11</sup></td></tr><tr><td>WDR63</td><td>WD repeat domain 63</td><td>10.0</td><td>5.3</td><td>26.3</td><td>0.0005</td><td>9.6</td><td>3.5</td><td>66.4</td><td>2.21 x 10<sup>-12</sup></td></tr><tr><td>ARMC3</td><td>armadillo repeat containing 3</td><td>9.7</td><td>5.0</td><td>26.0</td><td>0.0002</td><td>9.1</td><td>2.9</td><td>73.2</td><td>4.72 x 10<sup>-11</sup></td></tr><tr><td>ZMYND10</td><td>zinc finger, MYND-type containing 10</td><td>9.0</td><td>4.5</td><td>23.5</td><td>0.0002</td><td>9.8</td><td>3.8</td><td>62.8</td><td>5.05 x 10<sup>-12</sup></td></tr><tr><td>STOML3</td><td>stomatin (EPB72)-like 3</td><td>9.3</td><td>4.8</td><td>23.4</td><td>7.92 x 10<sup>-5</sup></td><td>10.4</td><td>3.5</td><td>121.7</td><td>8.47 x 10<sup>-13</sup></td></tr><tr><td>SPEF2</td><td>sperm flagellar 2</td><td>9.3</td><td>4.8</td><td>22.0</td><td>1.46 x 10<sup>-5</sup></td><td>7.2</td><td>3.2</td><td>15.6</td><td>2.52 x 10<sup>-9</sup></td></tr><tr><td>CFAP52</td><td>cilia and flagella associated protein 52</td><td>12.0</td><td>7.6</td><td>21.3</td><td>0.0013</td><td>10.6</td><td>4.3</td><td>81.3</td><td>1.07 x 10<sup>-11</sup></td></tr></tbody></table>

<table><caption>**Top 15 Human and Mouse EZ Transcription Factors**</caption><thead><tr><th>Gene</th><th>Description</th><th>hEZ (log<sub>2</sub>)</th><th>Peri hEZ (log<sub>2</sub>)</th><th>Fold Change</th><th>p Value</th><th>mEZ (log<sub>2</sub>)</th><th>Peri mEZ (log<sub>2</sub>)</th><th>Fold Change</th><th>p Value</th></tr></thead><tbody><tr><td>ARX</td><td>aristaless-related homeobox</td><td>8.6</td><td>4.6</td><td>16.4</td><td>0.0012</td><td>6.7</td><td>3.3</td><td>10.6</td><td>2.15 x 10<sup>-9</sup></td></tr><tr><td>RFX2</td><td>regulatory factor X, 2</td><td>9.6</td><td>5.6</td><td>16.2</td><td>0.0018</td><td>7.6</td><td>3.5</td><td>17.3</td><td>5.10 x 10<sup>-8</sup></td></tr></tbody></table>

(Continued on next page)
<figure><img src="image_16.png" /></figure>

Please cite this article in press as: Ghazale et al., RNA Profiling of the Human and Mouse Spinal Cord Stem Cell Niches Reveals an Embryonic-like Regionalization with MSX1⁺ Roof-Plate-Derived Cells, Stem Cell Reports (2019), https://doi.org/10.1016/j.stemcr.2019.04.001

<table><caption>**Table 1. Continued**</caption><caption>Top 15 Human and Mouse EZ Transcription Factors</caption><thead><tr><th>Gene</th><th>Description</th><th>hEZ (log<sub>2</sub>)</th><th>Peri hEZ (log<sub>2</sub>)</th><th>Fold Change</th><th>p Value</th><th>mEZ (log<sub>2</sub>)</th><th>Peri mEZ (log<sub>2</sub>)</th><th>Fold Change</th><th>p Value</th></tr></thead><tbody><tr><td>SOX6</td><td>SRY box 6</td><td>10.4</td><td>8.3</td><td>4.5</td><td>0.0016</td><td>9</td><td>5.3</td><td>13.1</td><td>1.84 × 10<sup>-8</sup></td></tr><tr><td>PAX6</td><td>paired box 6</td><td>9.7</td><td>7.5</td><td>4.5</td><td>0.0385</td><td>7.6</td><td>3.8</td><td>13.6</td><td>4.18 × 10<sup>-10</sup></td></tr><tr><td>NFIB</td><td>nuclear factor I/B</td><td>9.4</td><td>7.2</td><td>4.4</td><td>0.0582</td><td>9.3</td><td>5.7</td><td>12.8</td><td>1.18 × 10<sup>-8</sup></td></tr><tr><td>ID4</td><td>inhibitor of DNA binding 4</td><td>11.6</td><td>9.5</td><td>4.4</td><td>0.0724</td><td>11.8</td><td>8.5</td><td>10.1</td><td>2.79 × 10<sup>-8</sup></td></tr><tr><td>MYB</td><td>v-myb avian myeloblastosis viral oncogene homolog</td><td>6.4</td><td>4.3</td><td>4.3</td><td>0.0001</td><td>10.2</td><td>4.4</td><td>59</td><td>1.48 × 10<sup>-11</sup></td></tr><tr><td>NR4A3</td><td>nuclear receptor subfamily 4, member 3</td><td>6.3</td><td>4.3</td><td>4</td><td>0.0037</td><td>11.2</td><td>5.7</td><td>45.6</td><td>1.02 × 10<sup>-11</sup></td></tr><tr><td>REST</td><td>RE1-silencing transcription factor</td><td>9.1</td><td>7.2</td><td>3.8</td><td>0.0032</td><td>7.1</td><td>4.1</td><td>8</td><td>2.53 × 10<sup>-9</sup></td></tr><tr><td>KDM3A</td><td>lysine (K)-specific demethylase 3A</td><td>8.3</td><td>6.5</td><td>3.6</td><td>0.0262</td><td>8.9</td><td>5.8</td><td>8.9</td><td>1.86 × 10<sup>-9</sup></td></tr><tr><td>SOX2</td><td>SRY box 2</td><td>12.5</td><td>10.7</td><td>3.5</td><td>0.0301</td><td>12.1</td><td>9.2</td><td>7.3</td><td>3.92 × 10<sup>-8</sup></td></tr><tr><td>SOX9</td><td>SRY box 9</td><td>11.5</td><td>9.8</td><td>3.4</td><td>0.0255</td><td>9</td><td>4.4</td><td>24.2</td><td>2.16 × 10<sup>-10</sup></td></tr><tr><td>NR2F1</td><td>nuclear receptor subfamily 2, member 1</td><td>11.4</td><td>9.7</td><td>3.3</td><td>0.0165</td><td>11</td><td>7.8</td><td>8.9</td><td>1.46 × 10<sup>-9</sup></td></tr><tr><td>SALL1</td><td>spalt-like transcription factor 1</td><td>11.5</td><td>9.9</td><td>3.1</td><td>0.0278</td><td>9.9</td><td>4.6</td><td>38</td><td>2.76 × 10<sup>-10</sup></td></tr><tr><td>RFX3</td><td>regulatory factor X, 3</td><td>5.4</td><td>3.8</td><td>2.9</td><td>0.0031</td><td>9.2</td><td>5.2</td><td>16.5</td><td>1.01 × 10<sup>-9</sup></td></tr></tbody></table>

<table><caption>**Top 15 Human and Mouse EZ Transport-related genes**</caption><thead><tr><th>Gene</th><th>Description</th><th>hEZ (log<sub>2</sub>)</th><th>Peri hEZ (log<sub>2</sub>)</th><th>Fold Change</th><th>p Value</th><th>mEZ (log<sub>2</sub>)</th><th>Peri mEZ (log<sub>2</sub>)</th><th>Fold Change</th><th>p Value</th></tr></thead><tbody><tr><td>CFTR</td><td>cystic fibrosis transmembrane conductance regulator</td><td>9.5</td><td>5.1</td><td>21.2</td><td>3.86 × 10<sup>-6</sup></td><td>8.5</td><td>3.4</td><td>34.8</td><td>6.82 × 10<sup>-11</sup></td></tr><tr><td>SLC44A1</td><td>solute carrier family 44, member 1</td><td>9.5</td><td>7.3</td><td>4.5</td><td>0.0082</td><td>8.9</td><td>6.5</td><td>5.3</td><td>2.58 × 10<sup>-8</sup></td></tr><tr><td>OCA2</td><td>oculocutaneous albinism II</td><td>8.2</td><td>6.2</td><td>4.1</td><td>0.0026</td><td>6.2</td><td>4.0</td><td>4.6</td><td>6.24 × 10<sup>-6</sup></td></tr><tr><td>SLC15A2</td><td>solute carrier family 15, member 2</td><td>8.2</td><td>6.3</td><td>3.8</td><td>0.0069</td><td>7.5</td><td>4.7</td><td>7.1</td><td>0.0003</td></tr><tr><td>SLC38A6</td><td>solute carrier family 38, member 6</td><td>7.5</td><td>5.6</td><td>3.6</td><td>0.0842</td><td>7.0</td><td>5.5</td><td>3.0</td><td>5.09 × 10<sup>-5</sup></td></tr><tr><td>SLC40A1</td><td>solute carrier family 40, member 1</td><td>10.2</td><td>8.6</td><td>3.1</td><td>0.0473</td><td>7.2</td><td>5.8</td><td>2.7</td><td>0.0083</td></tr><tr><td>LRP6</td><td>low-density lipoprotein receptor-related protein 6</td><td>10</td><td>8.4</td><td>3</td><td>0.0099</td><td>7.7</td><td>6.4</td><td>2.4</td><td>0.0016</td></tr><tr><td>FOLR1</td><td>folate receptor 1 (adult)</td><td>6.6</td><td>5.1</td><td>2.9</td><td>0.0164</td><td>9.0</td><td>4.3</td><td>25.0</td><td>4.04 × 10<sup>-11</sup></td></tr></tbody></table>

(Continued on next page)
<figure><img src="image_17.png" /><figcaption>Table 1. Continued</figcaption></figure>

<table><caption>Top 15 Human and Mouse EZ Transport-related genes</caption><thead><tr><th>Gene</th><th>Description</th><th>hEZ (log<sub>2</sub>)</th><th>Peri hEZ (log<sub>2</sub>)</th><th>Fold Change</th><th>p Value</th><th>mEZ (log<sub>2</sub>)</th><th>Peri mEZ (log<sub>2</sub>)</th><th>Fold Change</th><th>p Value</th></tr></thead><tbody><tr><td>SLC16A9</td><td>solute carrier family 16, member 9</td><td>8.7</td><td>7.2</td><td>2.9</td><td>0.006</td><td>8.1</td><td>5.6</td><td>5.4</td><td>8.88 × 10<sup>-7</sup></td></tr><tr><td>SLC22A4</td><td>solute carrier family 22, member 4</td><td>9.3</td><td>7.8</td><td>2.9</td><td>0.0151</td><td>7.0</td><td>4.0</td><td>7.7</td><td>3.30 × 10<sup>-7</sup></td></tr><tr><td>SLC35D1</td><td>solute carrier family 35 D1</td><td>8.5</td><td>7.1</td><td>2.6</td><td>0.0042</td><td>7.7</td><td>4.6</td><td>8.4</td><td>6.28 × 10<sup>-9</sup></td></tr><tr><td>SLC12A7</td><td>solute carrier family 12, member 7</td><td>8.6</td><td>7.4</td><td>2.4</td><td>0.0143</td><td>8.1</td><td>5.3</td><td>7.0</td><td>1.03 × 10<sup>-7</sup></td></tr><tr><td>SLC2A12</td><td>solute carrier family 2, member 12</td><td>6.8</td><td>5.6</td><td>2.4</td><td>0.0335</td><td>8.4</td><td>4.8</td><td>12.0</td><td>1.18 × 10<sup>-9</sup></td></tr><tr><td>SLC30A5</td><td>solute carrier family 30, member 5</td><td>5.8</td><td>4.6</td><td>2.3</td><td>0.0341</td><td>8.5</td><td>7.0</td><td>2.7</td><td>0.0018</td></tr><tr><td>ABCC1</td><td>ATP-binding cassette, member 1</td><td>9.3</td><td>8.2</td><td>2.2</td><td>0.0194</td><td>8.2</td><td>5.2</td><td>7.9</td><td>5.16 × 10<sup>-9</sup></td></tr></tbody><tfoot><tr><td colspan="10">The table also shows lists of transporters and transcription factors (top 15) enriched in the mouse and human EZ, mEZ, mouse EZ; hEZ, human EZ.</td></tr></tfoot></table>

(Arx, Foxa1, Sulf1) involved in SHH signaling, a pathway known to be crucial for ventral specification of spinal cord during development. The proteoglycan Decorin (Dcn) was also expressed in the ventral part, a gene positively regulated by SHH signaling (Ingram et al., 2008). To identify the cellular origin of Dcn, we performed double immunofluorescence (IF) for DCN and PKD2L1, a specific marker for CSF-N (Becker et al., 2018) (Figure S2). Indeed, co-localization indicated that these neurons express a high level of DCN, a result supported by recent single cell RNA-sequencing (RNA-seq) databases (Zeisel et al., 2018) (Figure S5A). Contrasting with ventral SHH signaling, a robust expression of two bone morphogenetic protein (BMP) morphogens (namely Bmp6 and Gdf10/Bmp3b) involved in the dorsal patterning of the developing spinal cord (Wilson and Maden, 2005) were detected in the dorsal EZ (Figure 1G and Table S5). Two genes, Prokr2 and Thbs2, modulating neural stem cell fate in the SVZ (Benner et al., 2013; Prosser et al., 2007), were also enriched in the dorsal part. The dorsal or ventral expression of Arx, Bmp6, Dcn, Foxa1, Gdf10, Prokr2, Sulf1, and Thbs2 in the EZ was confirmed in the Allen brain atlas (Figure S3A).

## The Mouse and Human EZ Show a Conserved Dorsal-Ventral Regionalization of Transcription Factor Expression

One hundred twenty-one TFs were found more specifically-expressed (fold change EZ/periEZ ≥ 2) in the human and mouse EZ (Figure 1E). Examples of expression in the Allen brain atlas are illustrated in Figure S1B. IF confirmed the presence of corresponding proteins for some of them (ARX, FOXJ1, ID4, MEIS2, MSX1, NFIA, PAX6, PBX1, SOX4, SOX9, and SOX11) both in human (17-year-old) and mice (Figures 2 and 3). The histological quality of the second human spinal cord (46-year-old) was reduced compared with that of the first patient, although we could confirm protein expression for FOXJ1, ID4, MEIS2, NFIA, and PAX6 in this sample (Figure S4A). Unexpectedly, some of these TFs showed a dorsal-ventral asymmetric expression. Both in human (17-year-old) and in mouse, PAX6 protein stained dorsal and lateral cells while ventral cells were negative (Figures 2 and 3), which was also observed in the GENSAT gene expression atlas (Figure S1B). Another clear regionalized-expression was found for MSX1 protein in mouse and in human (17-year-old). In the human sample, MSX1 was weakly expressed by ependymal cells but cells in the dorsal part had a much stronger staining (Figure 2). In mice, MSX1 protein was confined to a few cells constituting the roof of the EZ (Figure 3), which was confirmed in gene expression atlases (Figure S3B). These MSX1+ cells typically exhibit a higher level of ID4 staining in mice (Figure 3). Contrasting with
<figure><img src="image_18.png" /></figure>

<figure><img src="image_19.png" /></figure>

(legend on next page)
<figure><img src="image_20.png" /></figure>

the dorsally-expressed proteins, we detected a ventral expression of ARX, a TF expressed by the floor-plate during spinal cord development. ARX protein was confined to a group of ventral cells both in human and mouse EZ (Figures 2 and 3). This led us to investigate by IF the expression of FOXA2, another TF involved in floor plate formation (Cho et al., 2014). Foxa2 specificity for the EZ was under our selected threshold in mice and in human (fold change EZ/periEZ = 1.7 and 1.3, respectively) however the FOXA2 protein was detected in a group of ventral and ventral/ventro-lateral cells in human and mice respectively (Figures 2 and 3). In human, the same population of ventral cells expressed ARX and FOXA2 (Figure S4B) whereas in mice, FOXA2 was expressed by ventral CSF-N expressing DCN and PKD2L1 but did not express ARX (Figure S2). In addition to these evolutionarily-conserved TFs, we also detected genes whose expression appears to be specifically expressed in the mouse or the human EZ (Table S3). This was checked at the protein level for *Tal1* (also known as *Scl*), a well-known TF involved in hematopoiesis but also expressed during spinal cord development (Smith et al., 2002). IF for TAL1 revealed the presence of positive subpopulations in the mouse EZ (Figure S2) which co-labeled with PDK2L1, thus identifying these cells as CSF-N. This was confirmed by single cell CNS RNA-seq database (Zeisel et al., 2018) (Figure S5A). In contrast, no convincing staining for TAL1 and PDK2L1 proteins was detected in the two human spinal cord EZs (not shown), and DCN labeling was restricted to vessels in the parenchyma (Figure S5B) as expected from the literature (Jarvelainen et al., 2015). A recent single cell CNS RNA-seq database (Zeisel et al., 2018) identified six markers for mouse CSF-N (*Crct1*, *Dcn*, *Espn*, *Pkd1l2*, *Pkd2l1*, and *Pdzk1ip1*), which we found enriched at the transcriptional level in the mouse EZ but not in the human EZ (Figure S5C).

Collectively, these results indicate a conserved dorsal-ventral regionalization of TF expression in the mouse and human EZ; however, no evidence could be found for the presence of CSF-N in human.

## Vegfr3-YFP Transgenic Mice Reveal Distinct MSX1⁺ and ARX⁺ Cells in the Roof and Floor of EZ

The regionalized expression of TFs in the EZ was indicative of specific cells located in the roof and floor of the EZ. This was also suggested by the specific expression of

the immature neural marker NES in these regions (Becker et al., 2018). The Vegfr3-YFP transgenic mouse has been used to detect subpopulations of neural precursor cells in the neurogenic SVZ (Calvo et al., 2011), which prompted us to use it in the spinal cord EZ context. Indeed, subpopulations of EZ cells expressing YFP were observed in all sections and were specifically enriched at the roof and floor portions (Figures 4A and 4B). These cells were positive for TFs strongly expressed by ependymal cells such as FOXJ1 and SOX2 while only the dorsal YFP⁺ cells express PAX6 (Figure 4C). In addition, in the roof part, 67% (41 among 61 cells, two mice) of YFP⁺ cells expressed the MSX1 TF (Figure 4C), which was further confirmed by generating a double transgenic mouse (Vegfr3-YFP × Msx1-Tomato, see below, Figure 6F) in which approximately 77% of the dorsal YFP⁺ cells were also Tomato⁺ (37 among 48 cells, two mice). It was previously reported that the c-RET receptor is present in a subpopulation of radial dorsal cells in the EZ (Pfenninger et al., 2011), and we observed that 71% of YFP⁺ (20 among 28 cells, two mice) were positive for this receptor (Figure 4C). In the EZ floor, 90% of YFP⁺ cells (60 among 66 cells, three mice) expressed the ARX TF (Figure 4C). Finally, some dorsal Vegfr3-YFP⁺ cells were clearly positive for GFAP and NES intermediate filaments (Figure S6A), whereas we could not establish this with confidence for ventral cells. These results support the existence of distinct EZ roof and floor cells expressing *Vegfr3* and specific TFs.

## Roof MSX1⁺ and Floor ARX⁺ Ependymal Cells are Born Early during Development

We next questioned the developmental origin of MSX1⁺ and ARX⁺ cells by performing IF at multiple developmental stages (embryonic day 13 [E13], E18, postnatal day 1 [P1], P6, P22, and P38) (Figure 5A). At E13, a group of MSX1⁺ cells was detected at the dorsalmost part of the developing spinal cord while a group of ARX⁺ cells were present in the ventralmost region. These cells, situated at positions corresponding respectively to the embryonic floor and roof plates, had elongated nuclei oriented along the dorsal-ventral axis, evocating migration. At E18, ARX⁺ and MSX1⁺ cells appeared to have migrated centrally to delimit a small group of cells presumptive of the adult EZ. At P1 and P6, the number of ARX⁺ and MSX1⁺ cells was reduced and at P22 and P38 only a few cells remained in the roof and floor of the EZ (Figure 5A). To analyze

## Figure 2. Human EZ Characterization

IF of the indicated proteins in the human EZ (patient aged 17 years). White arrowheads show negative cells indicative of the staining specificity. Yellow arrows show ventral ARX⁺ cells, dorsal MSX1⁺ cells, and ventral PAX6⁻ cells. Note that SOX4 is mostly cytoplasmic. Images are oriented with ventral part at bottom. Labelings were performed at the thoracic (control, ARX, CD24, FOXA2, MSX1, SOX4, SOX9) or lumbar level (CTNNB1, FOXJ1, ID4, MEIS2, NFIA, PAX6, PBX1, SOX11). These images are representative of at least eight sections. Scale bar, 50 µm (applies to all images).
<figure><img src="image_21.png" /></figure>

<figure><img src="image_22.png" /></figure>

(legend on next page)
<figure><img src="image_23.png" /></figure>

further the developmental origin of adult EZ MSX1⁺ cells, we used a genetic tracing approach based on a Msx1-CreERT2/Rosa-Loxed Tomato transgenic line (hereafter referred to as Msx1-Tomato [Lallemand et al., 2013]). Tamoxifen injection in adult mice revealed the presence of Tomato⁺ cells residing in the EZ roof (Figures 6A and S6B), which validated this tool. Tomato⁺ cells were very rarely observed in the ventral part of the EZ or in the parenchyma (not shown). Tamoxifen was then injected in pregnant transgenic mice at E11.5 to permanently label embryonic MSX1⁺ cells and explore their fate (Figure 5B). Analysis of sections taken from cervical and lumbar spinal cord levels of P30 mice (n = 2 mice) revealed the constant presence of Tomato⁺ cells situated in the EZ roof (Figure 5B). Altogether, these data demonstrate an early embryonic origin of ARX⁺ and MSX1⁺ spinal cord EZ cells.

## Roof MSX1⁺ Cells Are Radial Quiescent Cells

We previously reported an enrichment of neurosphere-forming cells in the dorsal half of the EZ (Sabourin et al., 2009). This was associated with a high content of GFAP⁺ radial glial-like cells. As MSX1⁺ cells were found almost exclusively in the roof of the EZ, further characterization of these cells was done. This was also motivated by the reported expression of *Msx1/2* in some stem cells during eye development (Bélanger et al., 2017) and its role in promoting regeneration of amputated tail including spinal cord (Beck et al., 2003) in *Xenopus*. The intense fluorescence observed in Msx1-Tomato mice enabled morphological characterization of these cells, notably using clarification of the whole spinal cord and 3D reconstruction. Tomato⁺ cells have a long radial morphology and send their process toward the pial surface (Figures 6A, 6B, and S6B). Their soma make contact with the lumen, but radial Tomato⁺ cells can also be observed more dorsally and at distance from the EZ (Figure 6A, yellow arrow). Horizontal sections also revealed the existence of Tomato⁺ processes running longitudinally along the EZ roof (Figure 6C).

By performing IF for c-RET, FOXJ1, GFAP, NES, and SOX2, we observed that most Tomato⁺ cells also express these markers (Figures 6D and S6C). We next crossed Msx1-Tomato mice with hGFAP-GFP transgenic mice, which have been widely used to purify adult neural stem cells (Nolte et al., 2001). Contrasting with GFAP immunostaining (Figure S6C), only a small fraction of Tomato⁺ cells

(15%, 8 double-positive cells among 54 Tomato⁺ cells; n = 2 mice) were GFP⁺ (Figure 6E). Often GFP⁺ and Tomato⁺ cells were found to be associated, suggesting close interactions between them. We then explored the proliferation rate of Tomato⁺ cells by injecting tamoxifen for 5 days in Msx1-Tomato mice to label MSX1⁺ cells, then 5-ethynyl-2'-deoxyuridine (EdU) injections were performed for 5 additional days (twice a day) before sacrificing on day 11 (Figure 6G). As previously reported, few EdU⁺ cells were found in the EZ (approximately one cell per section), however no Tomato⁺ EdU⁺ could be observed (Figure 6G). This indicates that MSX1⁺ cells are quiescent or proliferate at a much slower rate compared with the other ependymal cells.

## DISCUSSION

In this article, RNA profiling was used to generate a cellular and molecular resource for the adult human and mouse spinal cord EZ. Our findings reveal important and new characteristics of this poorly defined region.

In accordance with the presence of neural stem cells in the EZ, this region shows enrichment for smoothened/SHH and Hippo/YAP signaling genes, which are involved in stem cell maintenance (Alvarez-Buylla and Ihrie, 2014; Mo et al., 2014). Echoing their localization at the interface between CSF and the spinal cord parenchyma, ependymal cells also highly express genes involved in transport and ciliogenesis including three TFs of the RFX family and FOXJ1 (Figure 1E). Compared with brain, spinal cord ependymal cells proliferate (Alfaro-Cervello et al., 2012; Pfeninginger et al., 2011), which is also reflected at the RNA level by enrichment for genes involved in cell division (Figure 1E; Tables S1 and S2). A set of 120 TFs conserved between mice and humans was identified and the expression at the protein level was confirmed for 12 of them. MEIS2 and PBX1 expressions are particularly interesting, as these TFs dimerize and are essential regulators of adult SVZ neurogenesis (Grebbin et al., 2016). Ependymal cells also express high levels of NFIA and SOX9, two TFs involved in gliogenesis during development (Kang et al., 2012). Post lesion, these cells mostly generate astrocytes (Barnabé-Heider et al., 2010) and the expression of NFIA and SOX9 may be responsible for their glial fate restriction.

One unexpected observation was the conserved and regional expression of four homeodomain-containing TFs

### Figure 3. Mouse EZ Characterization

IF for the indicated proteins in the adult mouse EZ (lumbar level). Images are oriented with ventral part at bottom. The yellow arrows on PAX6 staining show negative cells in the ventral part. Note the higher expression of ID4 in cells localized in the dorsal part, which also express MSX1. FOXA2⁺ and TAL1⁺ cells are localized in a subependymal position and express PKD2L1 a marker specific for CSF-N (Figure S2). Note that SOX4 is mostly cytoplasmic. These images are representative of ten sections per animal, n = 4 mice analyzed. Scale bar, 20 µm (applies to all images).
<figure><img src="image_24.png" /><figcaption>A
Vegfr3-YFP/Hoechst</figcaption></figure>

<figure><img src="image_25.png" /><figcaption>B
Position of YFP+ cells</figcaption></figure>

<figure><img src="image_26.png" /><figcaption>C
ARX/YFP/Hoe.
c-RET/YFP/Hoe.
FOXJ1/YFP
MSX1/YFP/Hoe.
PAX6/YFP
SOX2/YFP</figcaption></figure>

<figure><img src="image_27.png" /></figure>

(legend on next page)
<figure><img src="image_28.png" /></figure>

in the niche, namely ARX, FOXA2, MSX1, and PAX6. This
situation is reminiscent of the developing spinal cord neu-
roepithelium where cells along the dorsal-ventral axis ex-
press different combinations of homeogenes (Wilson and
Maden, 2005). Both in human and mouse, ARX and
FOXA2 are expressed by ventralmost cells that do not co-
express PAX6. Of note, whereas FOXA2 was expressed by
ventral CSF-N in mice (Figure S2 and Petracca et al.,
2016), no such cells have been identified in human so far.
The two other markers we identified in mice for these cells,
DCN, an extracellular matrix protein (Jarvelainen et al.,
2015), and TAL1, a hematopoietic transcription factor,
were not observed in the human EZ. This indicates a diver-
gent organization of the ventral EZ between rodents and
humans. In the dorsal portion of the EZ, high MSX1 expres-
sion was restricted to few cells located in the EZ roof both in
mouse and human. By tracing the origin of ARX⁺ and
MSX1⁺ cells during spinal cord development in mice, we
observed that they are already present at E13.5 in the spinal
cord roof and floor plates. Their number is then reduced
and they appear to migrate centrally to generate the roof
and floor of the adult EZ. The early developmental origin
(at least E11.5) of MSX1⁺ cells was further demonstrated us-
ing genetic tracing. Thus, in contradiction to what was pre-
viously established (Fu et al., 2003), ependymal cells are
not entirely derived from the ventral neuroepithelium
but also incorporate cells from the dorsal part of the devel-
oping spinal cord. Another distinguishing feature of the
dorsal and ventral EZ cells is the expression of *Vegfr3* (also
known as *Flt4*, the receptor for VEGFC), which was
observed using *Vegfr3*-YFP transgenic mice. Part of these
*Vegfr3*-YFP⁺ cells express MSX1 and ARX. During brain
development *Vegfr3* is expressed in the ventricular zone
and radial glial cells (Ward and Cunningham, 2015). It is
also expressed in brain adult neural stem cells where it reg-
ulates their activation and proliferation (Calvo et al., 2011).
A similar role for *Vegfr3* may apply for adult spinal cord
ependymal cells.

These results provide evidence for the persistence of em-
bryonic floor and roof plate cells in the adult spinal cord
niche. Recent studies, based on a Wnt-reporter mouse la-
beling dorsal neural tube cells (Xing et al., 2018) and on a

Nato3-reporter mouse labeling floor plate cells (Khazanov et al., 2017), reached similar conclusions. During development, these floor and roof plate cells secrete morphogens such as SHH and BMP6, acting as growth factors and patterning signals (Wilson and Maden, 2005). Our results (Figure 1G) and expression atlas (Figure S3A) indicate that dorsal EZ cells express morphogen genes (Bmp6 and Gdf10). This suggests a specific role for these cells in the spinal cord niche, which warrants further exploration.

As the presence of MSX1⁺ cells has not been reported in
the spinal cord EZ, we characterized these cells further
and observed that they have distinct features compared
with the other cells of the niche. They have a long radial
morphology and express FOXJ1, GFAP, NES, and SOX2 pro-
teins. They also express c-RET, the receptor for GDNF and
NTN growth factors, which is important for hematopoietic
stem cells and neural crest cells (Kubota et al., 2004). By
crossing Msx1-CreERT2/Tomato and hGFAP-GFP trans-
genic mice, only a minority of Tomato⁺ cells were double-
positive, thus revealing the presence of a cellular heteroge-
neity in the dorsal part of the niche and of MSX1⁺ cells.

The existence of subpopulations of quiescent neural
stem cells has been reported in the hippocampus and the
SVZ niches (Codega et al., 2014; Llorens-Bobadilla and
Martin-Villalba, 2017; Llorens-Bobadilla et al., 2015).
Importantly, data mining of single cell RNA-seq analysis
performed in the SVZ revealed that *Msx1* and *Id4* are highly
enriched in quiescent neural stem cells (Llorens-Bobadilla
et al., 2015) (Figure S5D). These two TFs are regulated by
the BMP signaling (Ramos and Robert, 2005; Ruzinova
and Benezra, 2003) and, indeed, expression of *Bmp6* and
*Bmp3b/Gdf10* is also restricted to SVZ quiescent stem cells
(Figure S5D). A similar situation may also be present in
the dorsal part of the spinal cord EZ. Indeed, we found
that *Bmp6*, *Bmp3b/Gdf10*, and *Msx1* are highly expressed
in the dorsal part of the spinal cord EZ and that MSX1⁺ cells
express a higher level of ID4 transcription factor (Figure 3).
EdU incorporation also revealed that MSX1⁺ cells are less
proliferative than the other cells of the niche (Figure 6G).
This suggests that dorsal radial MSX1⁺ cells may behave
as quiescent neural stem cells in the adult spinal cord.
Further work is needed to support this hypothesis.

**Figure 4. EZ in Vegfr3-YFP Mice**

(A) Representative images of the spinal cord EZ (lumbar level) in *Vegfr3*-YFP mice (IF anti-GFP). Images are oriented with ventral part at bottom. YFP⁺ cells are mainly present in the dorsal and ventral regions (yellow arrows). Boxed areas show dorsal cells sending a process toward the lumen. These images are representative of at least 20 sections per animal (n = 3 mice analyzed). Ventral YFP⁺ cells are present in >95% of sections whereas dorsal YFP⁺ cells are present in approximately 25% of examined sections.

(B) Quantification of YFP⁺ cells in the lumbar EZ (187 YFP⁺ cells counted) indicated a preferential ventral and dorsal localization. One-way ANOVA + Tukey's post test (n = 3 mice).

(C) Phenotypic characterization of YFP⁺ cells with indicated protein (images are representative of 20 sections, n = 3 mice). Yellow arrows show double-positive cells. Images at the bottom are high magnification of arrow-pointed areas. Hoe., Hoechst. Scale bar, 20 µm (applies to all images).
<figure><img src="image_29.png" /></figure>

<figure><img src="image_30.png" /><figcaption>**Figure 5.** MSX1 and ARX Expression in the Developing Spinal Cord</figcaption></figure>

(A) IF for MSX1 and ARX during spinal cord embryonic development (E13 and E18) and postnatal stages (P1, P6, P22, P38). All images are oriented with ventral part at bottom. Images are representative of ten sections (n = 2 embryos and pups analyzed per stage).

(B) Representative images (n = 20 sections each levels, two mice) of the EZ of a P30 Msx1-CreERT2/Rosa-Tomato mouse derived from an embryo subjected to tamoxifen at E11.5.

In summary, our results uncovered that the adult spinal cord EZ region is conserved, regionalized, and composed of a mosaic of cells with different embryonic origin and ex-

pressing different types of TFs. This corpus of knowledge on the organization and genes expressed in the EZ will help to explore this adult stem cell niche further and will also be
<figure><img src="image_31.png" /></figure>

<figure><img src="image_32.png" /></figure>

<figure><img src="image_33.png" /></figure>

(legend on next page)
<figure><img src="image_34.png" /></figure>

useful to shed light on ependymoma, a rare type of tumor that can arise in the human spinal cord EZ.

# EXPERIMENTAL PROCEDURES

## Human Samples

Human spinal cords were collected at the Montpellier Hospital from two organ-donor patients (17 [male] and 46 [female] years old, accidental death) in strict agreement with the French bioethics laws (articles L1232-1 and -6) and after approval by the French institution for organ transplantation. An informed consent from the families was obtained by the organ procurement organization for this study. Surgery was performed as described previously (Bauchet et al., 2013), and the thoracolumbar segments were immediately placed in liquid nitrogen before processing for microdissection and immunofluorescence.

## Animals

Mice were handled following the guidelines of the Animal Care and Use Committee of the National Institute of Health and Medical Research (INSERM) who approved this study in accordance with the European Council directive (2010/63/UE) for the protection and use of vertebrate animals. Adult CD1 mice (3 months, Charles River, France) were used for microdissection, RNA profiling, and histology. Msx1-CreERT2/Rosa-Loxed Tomato transgenic line (Lallemand et al., 2013) was obtained from Y. Lallemand (Pasteur Institute, Paris). To induce recombination in Msx1-CreERT2 animals, we injected 100 µL of tamoxifen (Sigma, T5648, 20 mg/mL dissolved in corn oil) intraperitoneally for 4-5 days. hGFAP-GFP (Nolte et al., 2001) and Vegfr3-YFP (Calvo et al., 2011) transgenic mice were obtained from Prof. H. Kettenmann (MDC, Berlin) and Dr. J.L. Thomas (ICM, Paris), respectively.

## Tissue Microdissection

For human and mouse, after collection the spinal cords were flash frozen in N₂ without chemical fixation. Frozen sections (30 µm thick, T9-T10 thoracic part for mouse and thoracolumbar for human) were obtained at -23°C to prevent RNA degradation using a CM3050S microtome (Leica Microsystems, Wetzlar, Germany) and were mounted on PEN-membrane 1-mm glass slides (P.A.L.M. Microlaser Technologies, Bernried, Germany) that had been pretreated to inactivate RNase. Sections were then fixed in a series of pre-cooled ethanol baths (40 s in 95%, 75% and 30 s in 50%), stained with cresyl violet 1% for 30 s, and dehydrated in a

series of pre-cooled ethanol baths (30 s in 50%, 75%, and 40 s in 95% and 100%). Immediately after dehydration laser microdissection was performed using a PALM MicroBeam microdissection system version 4.6 equipped with PALM RoboSoftware (P.A.L.M. Microlaser Technologies). Laser power and duration were adjusted to optimize capture efficiency, and microdissection was performed at 63× magnification. Samples were collected in adhesive caps (P.A.L.M. Microlaser Technologies). To limit RNA degradation, we collected samples for up to 15 min per slide and lysed microdissected tissue with 250 µL of lysis buffer (Promega, Madison, WI, USA). The samples were stored at -80°C until extraction was performed using the ReliaPrep RNA cell Miniprep System (Promega) according to the manufacturer's protocol and eluted with 14 µL of RNase-free water. The concentration of RNA was determined using Nanodrop 1000 and the integrity of RNA was determined using the RNA 6000 Pico Kit and Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA). The RNA integrity number was above 7/8.

## RNA Profiling and Bioinformatics Analysis

RNA profiling was performed using Affymetrix microarray technology. Hybridization targets were obtained following a double-amplification procedure according to the protocol developed by Affymetrix (GeneChipTwo-Cycle Eukaryotic Target Labeling Assay; Affymetrix, USA). A hybridization mixture containing 10 µg of biotinylated cRNA was generated. The biotinylated cRNA was hybridized to HT_HG-U133_Plus_PM (human) and HT_MG-430_PM (mouse) Affymetrix microarrays. Four mouse and two human samples (EZ and peri EZ regions each, Figure 1A) were analyzed in total. The microarrays were scanned using the Affymetrix Gene Atlas scanner. The data files were generated with Affymetrix Expression Console v1.2.1 and gene expression data were normalized with the GC-RMA algorithm. Gene expression profiles were analyzed using the Affymetrix TAC 4.0 software (Transcriptome Analysis Console). The filter criteria were set to a linear fold change ≥2 between EZ and peri EZ regions. Gene lists were analyzed with DAVID Bioinformatics Resources 6.83 for gene enrichment analysis (Huang et al., 2008).

## Human and Mouse Histology

Mice were anesthetized by intraperitoneal injection of sodium pentobarbital (50 mg/kg) and perfused intracardially with 10 mL of PBS followed by 50 mL of 4% formaldehyde-PBS solution (pH 7.0). After dissection, spinal cords were post-fixed in the same solution for 1 h at 4°C and cryopreserved by successive immersion in 10%, 20%, and 30% sucrose solutions in PBS for at least 6 h. Cervical, thoracic and lumbar parts of the spinal cord were cut,

<caption>**Figure 6. Characterization of Dorsal MSX1⁺ Cells**</caption>

(A) Radial morphology of MSX1⁺ cells observed in Msx1-CreERT2/Rosa-Tomato mice. Yellow arrow shows rare radial Tomato⁺ cells outside EZ. Images are oriented with ventral part at bottom.
(B) 3D reconstruction of Tomato⁺ cells (coronal, intermediate, and lateral views).
(C) Dorsal view of the dorsal EZ region showing rostral-caudal oriented Tomato⁺ processes.
(D) Expression of c-RET receptor in Tomato⁺ cells (n = 10 sections).
(E and F) Representative images (n = 20 sections, two mice) of the EZ in double transgenic hGFAP-GFP/Msx1-CreERT2/Rosa-Tomato (E) and Vegfr3-YFP/Msx1-CreERT2/Rosa-Tomato (F) mice. White arrows indicate double-positive cells and yellow arrows Tomato⁺-only cells.
(G) EdU incorporation (5 days) in the EZ. Images show EdU⁺ cells (yellow arrows) not positive for Tomato (white arrows) in Msx1-CreERT2/Rosa-Tomato mouse sections. Quantifications are provided on right-hand graph (112 spinal cord sections, 5,120 cells examined, unpaired t test, n = 4 mice).
<figure><img src="image_35.png" /></figure>

embedded in OCT medium, rapidly frozen in liquid N₂-cooled isopentane, and cryosectioned (14 µm) (Leica apparatus). For developmental studies using time-mated embryos (Figure 5A), the day the plug was found was considered as E0.5. The embryos and pups (entire animal for P1, P6, and dissected spinal cord for P22) were fixed by direct immersion in 4% formaldehyde-PBS solution for 1 h then processed as for the adult spinal cord. P38 animals were perfused intracardially and post-fixed as for adults. To access the proliferation rate of Tomato+ cells (Figure 6G), we injected adult mice with EdU twice a day for 5 days (50 mg/kg), perfused them intracardially, and processed them for EdU staining (Baseclick kit).

For human, unfixed spinal cords at thoraco or thoracolumbar levels were directly cryosectioned (20 µm) and sections were immediately fixed by direct immersion in ice-cooled 4% formaldehyde-PBS solution for 20 min, followed by three washes with PBS.

Immunofluorescences were performed with primary antibodies (listed in Table S6) on sections permeabilized for 1 h with 0.1% Triton X-100 and 5% donkey serum. Secondary antibodies (Alexa 488- or Alexa 594-conjugated species-specific anti-mouse, -rabbit, or -goat) were purchased from The Jackson Laboratories. Incubations without primary antibody or with antibody recognizing antigens not present in the sections (monoclonal anti-DYKDDDDK tag or polyclonal antibodies against GFP) were used as negative controls. Nuclei (blue in all images) were stained with 1 µg/mL Hoechst for 10 min. The quality of staining was evaluated by two independent investigators (J.-P.H. and H.G. or C.R.). Images were taken using optical sectioning with structured illumination (Zeiss apotome microscope) or using a multi-photonic microscope (Zeiss LSM 7MP OPO) (Figure 6B). All presented images for mice are representative images, and the number of examined sections and animals are indicated in the figure legends. For human labelings, images are representative of at least eight sections for the first (17-year-old) patient and four sections for the second (46-year-old) patient. Clarification of mouse spinal cord was performed with protocol described in Tomer et al. (2014), and Imaris software was used for image treatments and 3D reconstruction.

## Statistical Analysis

All experiments and stainings were performed at least twice, most of them three times. Data are represented as means ± standard error of mean. Statistical differences in experiments were analyzed with tests indicated in the figure legends (GraphPad Prism software). Significance is denoted by \*\*\*p < 0.001, \*\*p < 0.01, and \*p ≤ 0.05.

## ACCESSION NUMBERS

Data are available at the functional genomics data Gene Expression Omnibus (GEO: GSE118445).

## SUPPLEMENTAL INFORMATION

Supplemental Information can be found online at https://doi.org/10.1016/j.stemcr.2019.04.001.

## AUTHOR CONTRIBUTIONS

C.R., H.G., N.L., D.M., and J.-P.H. performed and analyzed most of the experiments. J.-P.H. wrote the article. F.E.P., H.N.N., S.A., S.M.,

L.J., and B.R. contributed to spinal cord histology. Y.G. performed clarification and 3D reconstruction of spinal cord. L.B. performed surgery for human spinal cord and V.R. provided human spinal cord tissues. C.-F.C., J.-L.T., Y.L., and E.H. provided reagents and mouse models as well as intellectual input.

## ACKNOWLEDGMENTS

This work was supported by grants from IRP (Switzerland), IRME (France), AFM (France), ANR EU ERANET Neuroniche (J.-P.H.), and ANR Brainwash (J.-L.T., L.J.). H.G. was supported by an AFM PhD fellowship. We thank all Montpellier biocampus facilities for help (RHEM, MRI RIO, RAM) and excellent technical work. We are very grateful to Dr. H Boukhaddaoui (imaging), Dr. C Duperray (cytometry), M. Maistre (laser microdissection facilities, Bordeaux), V. Pantesco (Affymetrix facilities), and M. Goussard/P. Guigue (animals) for providing technical expertise in this work. We warmly thank Prof. Morohashi (Kyushu University, Japan) for anti-ArX antibody. The authors declare that there is no conflict of interest regarding the publication of this article.

Received: September 4, 2018

Revised: March 29, 2019

Accepted: April 1, 2019

Published: April 25, 2019

## REFERENCES

Adrian, E.K., Jr., and Walker, B.E. (1962). Incorporation of thymidine-H3 by cells in normal and injured mouse spinal cord. J. Neuropathol. Exp. Neurol. 21, 597-609.

Alfaro-Cervello, C., Soriano-Navarro, M., Mirzadeh, Z., Alvarez-Buylla, A., and Garcia-Verdugo, J.M. (2012). Biciliated ependymal cell proliferation contributes to spinal cord growth. J. Comp. Neurol. 520, 3528-3552.

Alvarez-Buylla, A., and Ihrie, R.A. (2014). Sonic hedgehog signaling in the postnatal brain. Semin. Cell Dev. Biol. 33, 105-111.

Barnabé-Heider, F., Göritz, C., Sabelström, H., Takebayashi, H., Pfrieger, F.W., Meletis, K., and Frisén, J. (2010). Origin of new glial cells in intact and injured adult spinal cord. Cell Stem Cell 7, 470-482.

Bauchet, L., Lonjon, N., Vachieri-Lahaye, F., Boularan, A., Privat, A., and Hugnot, J.-P. (2013). Isolation and culture of precursor cells from the adult human spinal cord. Methods Mol. Biol. 1059, 87-93.

Beck, C.W., Christen, B., and Slack, J.M. (2003). Molecular pathways needed for regeneration of spinal cord and muscle in a vertebrate. Dev. Cell 5, 429-439.

Becker, C.G., Becker, T., and Hugnot, J.-P. (2018). The spinal ependymal zone as a source of endogenous repair cells across vertebrates. Prog. Neurobiol. 170, 67-80.

Bélanger, M.-C., Robert, B., and Cayouette, M. (2017). Msx1-positive progenitors in the retinal ciliary margin give rise to both neural and non-neural progenies in mammals. Dev. Cell 40, 137-150.

Benner, E.J., Luciano, D., Jo, R., Abdi, K., Paez-Gonzalez, P., Sheng, H., Warner, D.S., Liu, C., Eroglu, C., and Kuo, C.T. (2013).
<figure><img src="image_36.png" /></figure>

Protective astrogenesis from the SVZ niche after injury is controlled by Notch modulator Thbs4. *Nature* **497**, 369.

Calvo, C.-F., Fontaine, R.H., Soueid, J., Tammela, T., Makinen, T., Alfaro-Cervello, C., Bonnaud, F., Miguez, A., Benhaim, L., and Xu, Y. (2011). Vascular endothelial growth factor receptor 3 directly regulates murine neurogenesis. *Genes Dev.* **25**, 831–844.

Cho, G., Lim, Y., Cho, I.T., Simonet, J.C., and Golden, J.A. (2014). Arx together with FoxA2, regulates Shh floor plate expression. *Dev. Biol.* **393**, 137–148.

Codega, P., Silva-Vargas, V., Paul, A., Maldonado-Soto, A.R., DeLeo, A.M., Pastrana, E., and Doetsch, F. (2014). Prospective identification and purification of quiescent adult neural stem cells from their in vivo niche. *Neuron* **82**, 545–559.

Dromard, C., Guillon, H., Rigau, V., Ripoll, C., Sabourin, J.-C., Perrin, F., Scamps, F., Bozza, S., Sabatier, P., and Lonjon, N. (2008). Adult human spinal cord harbors neural precursor cells that generate neurons and glial cells in vitro. *J. Neurosci.* Res. **86**, 1916–1926.

Fiorelli, R., Cebrian-Silla, A., Garcia-Verdugo, J.M., and Raineteau, O. (2013). The adult spinal cord harbors a population of GFAP-positive progenitors with limited self-renewal potential. *Glia* **61**, 2100–2113.

Fu, H., Qi, Y., Tan, M., Cai, J., Hu, X., Liu, Z., Jensen, J., and Qiu, M. (2003). Molecular mapping of the origin of postnatal spinal cord ependymal cells: evidence that adult ependymal cells are derived from Nkx6. 1+ ventral neural progenitor cells. *J. Comp. Neurol.* **456**, 237–244.

Gage, F.H., and Temple, S. (2013). Neural stem cells: generating and regenerating the brain. *Neuron* **80**, 588–601.

Garcia-Ovejero, D., Arevalo-Martin, A., Paniagua-Torija, B., Florensa-Vila, J., Ferrer, I., Grassner, L., and Molina-Holgado, E. (2015). The ependymal region of the adult human spinal cord differs from other species and shows ependymoma-like features. *Brain* **138**, 1583–1597.

Grebbin, B.M., Hau, A.-C., Groß, A., Anders-Maurer, M., Schramm, J., Koss, M., Wille, C., Mittelbronn, M., Selleri, L., and Schulte, D. (2016). Pbx1 is required for adult SVZ neurogenesis. *Development* **143**, 2281–2291.

Habib, N., Li, Y., Heidenreich, M., Swiech, L., Avraham-Davidi, I., Trombetta, J.J., Hession, C., Zhang, F., and Regev, A. (2016). DivSeq: single-nucleus RNA-Seq reveals dynamics of rare adult newborn neurons. *Science* **353**, 925–928.

Hamilton, L., Truong, M., Bednarczyk, M., Aumont, A., and Fernandes, K. (2009). Cellular organization of the central canal ependymal zone, a niche of latent neural stem cells in the adult mammalian spinal cord. *Neuroscience* **164**, 1044–1056.

Huang, D.W., Sherman, B.T., and Lempicki, R.A. (2008). Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. *Nat. Protoc.* **4**, 44–57.

Ingram, W., McCue, K., Tran, T., Hallahan, A., and Wainwright, B. (2008). Sonic Hedgehog regulates Hes1 through a novel mechanism that is independent of canonical Notch pathway signalling. *Oncogene* **27**, 1489–1500.

Jarvelainen, H., Sainio, A., and Wight, T.N. (2015). Pivotal role for decorin in angiogenesis. *Matrix Biol.* **43**, 15–26.

Kang, P., Lee, H.K., Glasgow, S.M., Finley, M., Donti, T., Gaber, Z.B., Graham, B.H., Foster, A.E., Novitch, B.G., and Gronostajski, R.M. (2012). Sox9 and NFIA coordinate a transcriptional regulatory cascade during the initiation of gliogenesis. *Neuron* **74**, 79–94.

Khazanov, S., Paz, Y., Hefetz, A., Gonzales, B.J., Netser, Y., Mansour, A.A., and Ben-Arie, N. (2017). Floor plate descendants in the ependyma of the adult mouse central nervous system. *Int. J. Dev. Biol.* **61**, 257–265.

Kubota, H., Avarbock, M.R., and Brinster, R.L. (2004). Growth factors essential for self-renewal and expansion of mouse spermatogonial stem cells. *Proc. Natl. Acad. Sci. U S A* **101**, 16489–16494.

Lallemand, Y., Moreau, J., Saint Cloment, C., Vives, F.L., and Robert, B. (2013). Generation and characterization of a tamoxifen inducible Msx1CreERT2 knock-in allele. *Genesis* **51**, 110–119.

Lein, E.S., Hawrylycz, M.J., Ao, N., Ayres, M., Bensinger, A., Bernard, A., Boe, A.F., Boguski, M.S., Brockway, K.S., Byrnes, E.J., et al. (2007). Genome-wide atlas of gene expression in the adult mouse brain. *Nature* **445**, 168–176.

Llorens-Bobadilla, E., and Martin-Villalba, A. (2017). Adult NSC diversity and plasticity: the role of the niche. *Curr. Opin. Neurobiol.* **42**, 68–74.

Llorens-Bobadilla, E., Zhao, S., Baser, A., Saiz-Castro, G., Zwadlo, K., and Martin-Villalba, A. (2015). Single-cell transcriptomics reveals a population of dormant neural stem cells that become activated upon brain injury. *Cell Stem Cell* **17**, 329–340.

Marichal, N., Reali, C., Trujillo-Cenóz, O., and Russo, R.E. (2017). Spinal cord stem cells in their microenvironment: the ependyma as a stem cell niche. *Adv. Med. Exp. Biol.* **1041**, 55–79.

Mo, J.S., Park, H.W., and Guan, K.L. (2014). The Hippo signaling pathway in stem cell biology and cancer. *EMBO Rep.* **15**, 642–656.

Mothe, A.J., Zahir, T., Santaguida, C., Cook, D., and Tator, C.H. (2011). Neural stem/progenitor cells from the adult human spinal cord are multipotent and self-renewing and differentiate after transplantation. *PLoS One* **6**, e27079.

Nolte, C., Matyash, M., Pivneva, T., Schipke, C.G., Ohlemeyer, C., Hanisch, U.K., Kirchhoff, F., and Kettenmann, H. (2001). GFAP promoter-controlled EGFP-expressing transgenic mice: a tool to visualize astrocytes and astrogliosis in living brain tissue. *Glia* **33**, 72–86.

Petracca, Y.L., Sartoretti, M.M., Di Bella, D.J., Marin-Burgin, A., Cagno, A.L., Schinder, A.F., and Lanuza, G.M. (2016). The late and dual origin of cerebrospinal fluid-contacting neurons in the mouse spinal cord. *Development* **143**, 880–891.

Pfenninger, C.V., Steinhoff, C., Hertwig, F., and Nuber, U.A. (2011). Prospectively isolated CD133/CD24-positive ependymal cells from the adult spinal cord and lateral ventricle wall differ in their long-term in vitro self-renewal and in vivo gene expression. *Glia* **59**, 68–81.

Prosser, H.M., Bradley, A., and Caldwell, M.A. (2007). Olfactory bulb hypoplasia in Prokr2 null mice stems from defective neuronal progenitor migration and differentiation. *Eur. J. Neurosci.* **26**, 3339–3344.

Ramos, C., and Robert, B. (2005). msh/Msx gene family in neural development. *Trends Genet.* **21**, 624–632.
<figure><img src="image_37.png" /></figure>

Ren, Y., Ao, Y., O'Shea, T.M., Burda, J.E., Bernstein, A.M., Brumm, A.J., Muthusamy, N., Ghashghaei, H.T., Carmichael, S.T., and Cheng, L. (2017). Ependymal cell contribution to scar formation after spinal cord injury is minimal, local and dependent on direct ependymal injury. *Sci. Rep.* 7, 41122.

Rosenberg, A.B., Roco, C.M., Muscat, R.A., Kuchina, A., Sample, P., Yao, Z., Graybuck, L.T., Peeler, D.J., Mukherjee, S., and Chen, W. (2018). Single-cell profiling of the developing mouse brain and spinal cord with split-pool barcoding. *Science* **360**, 176–182.

Ruzinova, M.B., and Benezra, R. (2003). Id proteins in development, cell cycle and cancer. *Trends Cell Biol.* **13**, 410–418.

Sabourin, J.C., Ackema, K.B., Ohayon, D., Guichet, P.O., Perrin, F.E., Garces, A., Ripoll, C., Charité, J., Simonneau, L., and Kettenmann, H. (2009). A mesenchymal-like ZEB1+ niche harbors dorsal radial glial fibrillary acidic protein-positive stem cells in the spinal cord. *Stem Cells* **27**, 2722–2733.

Smith, E., Hargrave, M., Yamada, T., Begley, C.G., and Little, M.H. (2002). Coexpression of SCL and GATA3 in the V2 interneurons of the developing mouse spinal cord. *Dev. Dyn.* **224**, 231–237.

Stenudd, M., Sabelström, H., and Frisén, J. (2015). Role of endogeanal neural stem cells in spinal cord injury and repair. *JAMA Neurol.* **72**, 235–237.

Tomer, R., Ye, L., Hsueh, B., and Deisseroth, K. (2014). Advanced CLARITY for rapid and high-resolution imaging of intact tissues. *Nat. Protoc.* **9**, 1682.

Ward, M., and Cunningham, A. (2015). Developmental expression of vascular endothelial growth factor receptor 3 and vascular endothelial growth factor C in forebrain. *Neuroscience* **303**, 544–557.

Weiss, S., Dunne, C., Hewson, J., Wohl, C., Wheatley, M., Peterson, A.C., and Reynolds, B.A. (1996). Multipotent CNS stem cells are present in the adult mammalian spinal cord and ventricular neuroaxis. *J. Neurosci.* **16**, 7599–7609.

Wilson, L., and Maden, M. (2005). The mechanisms of dorsoventral patterning in the vertebrate neural tube. *Dev. Biol.* **282**, 1–13.

Xing, L., Anbarchian, T., Tsai, J.M., Plant, G.W., and Nusse, R. (2018). Wnt/β-catenin signaling regulates ependymal cell development and adult homeostasis. *Proc. Natl. Acad. Sci. U S A* **115**, E5954–E5962.

Xu, W., Sachewsky, N., Azimi, A., Hung, M., Gappasov, A., and Morshead, C.M. (2017). Myelin basic protein regulates primitive and definitive neural stem cell proliferation from the adult spinal cord. *Stem Cells* **35**, 485–496.

Zeisel, A., Hochgerner, H., Lonnerberg, P., Johnsson, A., Memic, F., van der Zwan, J., Haring, M., Braun, E., Borm, L.E., La Manno, G., et al. (2018). Molecular architecture of the mouse nervous system. *Cell* **174**, 999–1014.e22.