Up-regulation of Hox genes in pons of Rnf220+/− mice. (A) The heatmap of RNA-seq data showing Hox genes expression in pons of WT or Rnf220+/− mice (n=2 mice per group). (B-C) UMAP diagram showing 15 identified cell clusters annotated by snRNA-seq analysis of pons. Each dot represents a single cell, and cells are laid out to show similarities (n=3 mice per group). (D) Heatmap of snRNA-seq data showing Hox expression changes in each cell cluster. (E) Quantitative real time polymerase chain reaction, qRT-PCR) analysis showing mRNA levels of indicated Hox genes in P19 cells when endogenous Rnf220 was knocked down by siRNAs in the presence of RA. WT, wild-type. HE, heterozygote. RA, retinoic acid. **p < 0.01, ***p < 0.001.

Differently expressed genes identified using microarray between WT and Rnf220+/− mice.

Differently expressed genes identified using microarray between WT and Rnf220−/− mice.

Differently expressed genes of pons identified using snRNA-seq between WT and Rnf220+/− mice.

Hox gene expression was dysregulated and motor cortex projections were disorganized in PN of Rnf220+/− mice. (A) qRT-PCR analysis of relative expression levels of Hox3, Hox4, and Hox5 in rostral, middle, and caudal sections of PN in WT and RNF220+/− mice. Expression level of each gene in rostral section of WT PN was set to 1. (n=4 mice per group). (B) Diagram of experimental stereotactic injections. (C, D) Fluorescent signals showing projection from motor cortex to PNs in WT and Rnf220+/− mice (n=2 mice per group). WT, wild-type; HE, heterozygote; **p < 0.01; ***p< 0.001.

RNF220 mediates WDR5 degradation. (A-D) Western blots analysis showing the protein level of WDR5 in the indicated brain tissues of mice with different genotypes at different ages. (E) Western blot analysis of protein levels of WDR5 in P19 cells with Rnf220 knockdown or not in the presence or absence of RA. IB, immunoblot. WT, wild-type. HE, heterozygote. KO, knockout. PN, pontine nuclei. NC, negative control. RA, retinoic acid.

RNF220 interacts with and targets WDR5 for K48-linked polyubiquitination. (A) Endogenous co-immunoprecipitation analysis showing the interaction between RNF220 and WDR5 in hindbrains of WT mice. (B-C) Co-immunoprecipitation (co-IP) analysis of interactions between RNF220 and WDR5 in HEK293 cells. HEK293 cells were transfected with indicated plasmids and harvested after 48 h. Cell lysates were immunoprecipitated with anti-FLAG beads. Whole-cell lysate and immunoprecipitates were subjected to western blot analysis using indicated antibodies. (D) Western blots analysis shows the protein level of WDR5 when co-expressed with wild type or mutated RNF220 in HEK293 cells. (E) Western blots analysis shows the protein level of WDR5 when co-expressed with RNF220 in HEK293 cells in the presence of MG132 or not. MG132 treatment reversed the RNF220-induced decrease of WDR5 protein. (F) In vivo ubiquitination assays showing the ubiquitination status of WDR5 in hindbrains of WT or Rnf220+/− mice. (G) In vivo ubiquitination assays showing the ubiquitination status of WDR5 when co-expressed with WT or mutated RNF220 in HEK293 cells. (H) In vivo ubiquitination assays showing RNF220-induced polyubiquitination of WDR5 when the indicated ubiquitin mutations were used in HEK293 cells. (I) In vivo ubiquitination assays showing the ubiquitination status of the indicated WDR5 mutants when co-expressed with WT or ligase dead RNF220 in HEK293 cells. WT, wild-type. HE, heterozygote. KO, knockout. IB, immunoblot. UB, ubiquitin. WCL, whole cell lysate. △Ring, RNF220 Ring domain deletion. W539R, RNF220 ligase dead mutation. K48, ubiquitin with all lysines except the K48 mutated to arginine. K48R, ubiquitin with the K48 was substituted by an arginine. 3M, substitution of lysines at the positions of 109, 112, and 120 in WDR5 with arginines simultaneously.

WDR5 recovered Rnf220 deficiency-induced up-regulation of Hox genes in P19 cell line. (A-B) qRT-PCR analysis showing the expression levels of Rnf220 (A) and Wdr5 (B) when transfected the indicated combinations of siRNAs against Rnf220 or Wdr5 in the presence or absence of RA. (C-D) qRT-PCR analysis showing the expression levels of HoxA1, HoxB1, HoxA9, HoxB9 when transfected siRnf220 or both siRnf220 and siWdr5 when treated with RA (D) or not (C). RA: retinoic acid. n.s., not significant. *p< 0.05, **p < 0.01, ***p < 0.001.

Genetic and pharmacological ablation of WDR5 recovered Rnf220 deficiency-induced up-regulation of Hox genes. (A) Diagram of experimental strategy for in utero local injection of WDR5 inhibitors. (B) qRT-PCR analysis of expression levels of Rnf220, HoxA9, and HoxB9 in hindbrains of WT and Rnf220+/− mouse embryos treated with WDR5 inhibitors or not at E18.5 (n=3 mice per group). (C) qRT-PCR analysis of expression levels of Rnf220, Wdr5, HoxA9, and HoxB9 in hindbrains of adult mice with indicated genotypes (n=2 mice per group). WT, wild-type. HE, heterozygote. ***p < 0.001.

Hox genes exhibited de-repression in hindbrain of Rnf220−/−and Rnf220+/− mouse embryos at late developmental stages. Heatmap of Hox expression in hindbrain of Rnf220−/− mouse embryos at E16.5 and E15.5. Expression of each Hox gene was analyzed by qRT-PCR against GAPDH. Each Hox gene in wild-type controls was set to 1 (n=2 mice per group). Blank color represents the expression of that Hox is low and exceed detection range. KO, knockout.

Hox genes were up-regulated in brainstem of Rnf220+/− mice. (A) Diagram of central neural system (CNS) in mice. (B) qRT-PCR analysis of Rnf220 expression in CNS of Rnf220+/− and control mice. (C-D) qRT-PCR analysis of HoxA9 (C) and HoxB9 (D) expression in each CNS section of Rnf220+/− and control mouse. (E) Diagram of mouse brain stem. (F) qRT-PCR analysis of mRNA levels of HoxA9 and HoxB9 in indicated sections of WT and Rnf220+/− mouse brainstems (n=2 mice per group). WT, wild-type. HE, heterozygote. SC, Spinal cord. n.s., not significant. *p < 0.05, **p < 0.01, ***p < 0.001.

Expression levels of Hox genes were not affected by Rnf220 knockdown in P19 cell line without RA induction. (A) qRT-PCR analysis of Rnf220 and indicated Hox genes expression after RA treatment in P19 cells. (B) qRT-PCR analysis of expression levels of HoxA1, HoxB1, HoxA9, and HoxB9 in Rnf220 knockdown P19 cells without RA induction. RA: retinoic acid. n.s., not significant. *p < 0.05, **p < 0.01, ***p < 0.001.

Hox genes were up-regulated in PN of Rnf220+/− mice. (A) Heatmap of snRNA-seq data showing expression levels of PN markers (Pax6, Barhl1, Unc5b, and Nfib) among 15 identified cell clusters. (B, C) qRT-PCR analysis of expression levels of indicated PN markers (B) and Hox genes (C) in PN of Rnf220+/− and WT mice (n= 3 mice per group). WT, wild-type. HE, heterozygote. n.s., not significant. *p < 0.05, ***p < 0.001.

Repressive epigenetic modification was down-regulated while activated epigenetic modification was up-regulated in promoter regions of indicated Hox genes in hindbrains of Rnf220+/− mice. (A-B) ChIP-qRT-PCR analysis of repressive epigenetic modification (H3K27me3) (A) and activated epigenetic modification (H3K4me3) (B) levels in promoter regions of indicated Hox genes in hindbrains of Rnf220+/− and WT mice (n=2 mice per group). WT, wild-type. HE, heterozygote. n.s., not significant. *p< 0.05, **p < 0.01, ***p < 0.001.

Protein levels of the indicated core components of PRC1 and PRC2 complex in indicated mouse brain tissues of different genotypes. (A) Western blot analysis of protein levels of core components of PRC1 and PRC2 complexes in hindbrain of WT, Rnf220+/−, and Rnf220−/− mouse embryos at E18.5. (B) Western blot analysis of protein levels of core components of PRC1 and PRC2 complexes in pons of adult Rnf220+/− and WT mice. (C, D) Western blot analysis of protein levels of WDR5 in cerebellum and cortex of adult Rnf220 +/−and WT mice. WT, wild-type; HE, heterozygote; KO, knockout; IB, immunoblot.

RNF220 interacted with and targeted WDR5 for polyubiquitination at multiple lysine sites. (A) Western blot analysis of levels of three WDR5 truncated proteins in HEK293 cells when co-transfected with RNF220 or not. (B) In vivo ubiquitination analysis of ubiquitination status of indicated WDR5 KR mutants when co-expressed with RNF220 or not in HEK293 cells. IB, immunoblot. IP, immunoprecipitation. WCL, whole cell lysate. K31R, K52R, K109R, K112R, K120R, K123R, or K126R, substitution of lysine with arginine in WDR5 at indicated positions.

WDR5 knockdown had no effect on Hox genes in P19 cell line. qRT-PCR analysis of mRNA levels of WDR5, HoxA1, HoxB1, HoxA9, and HoxB9 when WDR5 was knocked down by siRNA transfection in P19 cells with or without RA treatment. RA: retinoic acid. n.s., not significant. *p< 0.05.

Genetic ablation of Wdr5 reversed Hox gene de-repression caused by Rnf220 deficiency. qRT-PCR analysis of mRNA levels of HoxB8, HoxC4, HoxC5, HoxC6, HoxD8, HoxC9, and HoxD9 in hindbrains of WT, Wdr5fl/wt;Nestin-Cre; Rnf220fl/wt;Nestin-Cre, and Rnf220fl/wt;Wdr5fl/wt;Nestin-Cre mice (n=2 mice per group). WT, wild-type. n.s., not significant. ***p < 0.001.