Figures and data
Dynamic regulation of genes during oligodendrocytes lineage progression.
A) Schematic diagram showing magnetic-activated cell sorting (MACS)-based isolation of oligodendrocyte progenitor cells (OPCs), pre-oligodendrocytes (PreOLs), and oligodendrocytes (OLs) from young rat brains (2–3 months). The schema illustrates the order in which these cell types were most effectively isolated: mature OLs were isolated first using the OL-specific marker MOG, followed by OPCs using A2B5, and finally PreOLs using O4. B) Principal component analysis (PCA) (left) and average-linkage hierarchical clustering (right) based on global gene expression show that each cell type forms a distinct cluster. OPCs, PreOLs, and OLs were isolated from four different rats. C) Using gSWITCH, we identified four distinct patterns of switch genes during OL lineage progression. (B-C) GEO accession number of RNA-seq data: GSE301578.
Down-regulation of switch genes due to ageing.
A) Principal component analysis (PCA) shows that aged and young OPCs are distinguishable based on global gene expression (meta-analysis of RNA-seq datasets, GEO accession number: GSE303317). B) Volcano plot showing differentially expressed genes in aged OPCs compared to young OPCs. n= 6 (aged rats), 7 (young rats), FC: fold change; Padj: adjusted p-value (Benjamini–Hochberg FDR after Wald test). C) Venn diagram (top) showing the subset of Type 3 and Type 4 switch genes affected by ageing. Their expression patterns are shown in the heat map (bottom), where we observed down-regulation with ageing. * P < 0.05, ** P < 0.001, Chi-square test (two sided) with Yates’ correction.
Transcription factor Bcl11a-mediated network potentially affected duuring ageing.
A) Transcription factor binding site enrichment analysis using Type 3 and Type 4 switch genes as input identified 27 transcription factors (TFs) detected in both ENCODE and ReMap2, which are ChIP-seq–based experimental databases. Among these 27 TFs, we found 5 TFs that are themselves present as Type 3 or Type 4 switches (TF switches extracted from gSWITCH). ** P < 0.001, hypergeometric test. B) Heat map (top) showing expression (RNA-seq, GEO accession number: GSE303317) of these 5 TFs. Their differential expression in aged OPCs relative to young OPCs plotted as log2FC (bottom). FC: Fold change, n= 6 (aged rats), 7 (young rats), mean+SEM, * adjusted p-value < 0.05, ** adjusted p-value < 1×10⁻⁶, Benjamini–Hochberg FDR after Wald test. Note: All five TFs are classified as Type 3 switches, and Bcl11a shows the strongest repression in aged OPCs, with the highest statistical significance. C) Bcl11a, as a Type 3 TF switch, shows significantly repressed expression as OL lineage progresses. DESeq2-normalised values relative to OPCs plotted. n= 4 different rats, mean+SEM, ** adjusted p-value < 1×10⁻⁶, Benjamini–Hochberg FDR after Wald test. GEO accession number of RNA-seq data: GSE301578.
Effect of Bcl11a inhibition on OPCs proliferation and differentiation.
A) Reverse transcription (RT)–qPCR analysis of Bcl11a expression in OPCs transfected with siRNA targeting Bcl11a. siBcl11a: siGENOME Rat Bcl11a; siControl: siGENOME non-targeting control siRNA (Dharmacon™). B) Effect on proliferation. Left: Immunofluorescence analysis using antibodies to OLIG2 and EdU detection (Click-iT protocol), followed by EdU proliferation assay. Right: Quantification of EdU⁺OLIG2⁺ cells plotted as a percentage of OLIG2⁺ cells. P: p-value, Student’s t test (unpaired and two tailed) C) Effect on differentiation. Left: Immunofluorescence analysis using antibody to MBP, OLIG2 upon inhibition of Bcl11a in OPCs. Right: quantification of MBP+ cells were plotted as a percentage of OLIG2+ cells. (A-C) n = 3 independent experiments (each time 4 P7 brains were pooled for OPC isolation). mean+SEM. ** p<0.01, Student’s t test (unpaired, two tailed) (B-C) scale bar: 22 µm.
Effect of Bcl11a transient over-expression on proliferation and differentiation of aged OPCs.
A) Effect on proliferation. Left: Immunofluorescence analysis using antibodies to OLIG2 and EdU detection (Click-iT protocol), followed by EdU proliferation assay. Right: Quantification of EdU⁺OLIG2⁺ cells plotted as a percentage of OLIG2⁺ cells. P: p-value, Student’s t test (unpaired, two tailed). B) Effect on differentiation. Left: Immunofluorescence analysis using antibodies to MBP and OLIG2 after transfection of Bcl11a mRNA (mBcl11a) or Gfp mRNA (mGfp; control) in aged OPCs. Right: Quantification of MBP⁺ cells plotted as a percentage of OLIG2⁺ cells. * P<0.05, Student’s t test (unpaired, two tailed). (A-B) n = 3 independent experiments (each using OPCs isolated from pooled brains of 2 aged rats). mean ± SEM. scale bar: 22 µm.
OL lineage specific Bcl11a expression in aged mice.
A) Left: SOX10-driven construct carrying either Gfp alone (control) or Bcl11a linked to Gfp via a self-cleaving peptide element (P2A). This construct ensures SOX10-driven expression of both Bcl11a and Gfp. P2A: Porcine teschovirus-1 2A peptide, bGH p(A): bovine growth hormone poly(A) signal. The AAV containing this construct was delivered to 18-month-old mice via tail-vein injection. Right: schematic showing the experimental timeline. B) RNA in situ analysis of Plp1. In the lesioned area, the number of Plp1⁺ cells per mm² was counted and plotted. n=5-6 mice, mean+SEM, ** P<0.01, Student’s t test (unpaired, two tailed). scale bar: 100 µm. The white dotted line demarcates the lesion area.
Test for evolutionary selection on Bcl11a.
Schematic representation of the logic of the gSWITCH web application
(for detail see ‘Materials and Methods’ section). The gSWITCH application is accessible at https://altoslabs.shinyapps.io/gSWITCH/. Red circles indicate the start and the end state of a series (time or lineage stages or any other variable dependent series). GLM fit and hypothesis testing followed by contrast analysis between start and the end state, generates a subset of genes on which we apply maximum (max) or minimum (min) criteria within a closed interval while taking account of monotonicity trend as well as further specific contrasts between groups. This allows identification of four types of switches. Detail criteria are shown inside the boxes and the trend pattern are represented by coloured dashed-straight lines. Arc or solid curve lines indicate statistical significance for specific contrast. Adj p value: Adjusted p value [Benjamini–Hochberg FDR after LRT or Benjamini–Hochberg FDR after Wald test for specific contrast analysis]. u.d.: user defined value. Note that to run gSWITCH, RNA-seq data are required (bulk or scRNA-seq flattened to pseudo bulk), where the number of series or time points or group ≥3 and the replicates per time point or group ≥2. gSWITCH also identifies transcription factors (TFs) within the switches and provides a summary of the total number of switches and TFs per switch type.
Fidelity of protein expression from in vitro–transcribed Bcl11a mRNA; mRNA transfection efficiency in aged OPCs; and fidelity of AAV construct expression in the OL lineage and infection efficiency in aged mouse brain.
A) Western blot analysis using antibodies to BCL11A and β-actin (ACTB), with ACTB serving as a loading control, in MCF7 cells transfected with Bcl11a mRNA or Gfp mRNA (control). B) Immunofluorescence analysis using antibodies to OLIG2 and GFP in aged OPCs transfected with Gfp mRNA. Quantification of GFP+OLIG2+ cells were plotted as percentage of OLIG2+ cells. n = 3 independent experiments (each using OPCs isolated from pooled brains of 2 aged rats). mean ± SEM. Scale bar: 22 µm. C) AAV carrying constructs are shown in Figure 6A. Left: Immunofluorescence analysis of OLIG2 and GFP immunostaining. Right top: Quantification of OLIG2+GFP+ cells were plotted as a percentage of GFP+DAPI+ cells. Right bottom: quantification of OLIG2+GFP+ cells were plotted as a percentage to OLIG2+ cells. n= 3 mice, mean+SEM. P: p-value, Student’s t test (unpaired, two tailed). Scale bar: 20 µm.
Immunofluorescence analysis using antibody to OLIG2 after OL lineage specific Bcl11a expression in aged mice.
The AAV containing the Sox10-driven construct carrying either Gfp alone (control) or Bcl11a linked to Gfp via a self-cleaving peptide element (P2A) was delivered by tail-vein injection (see Figure 6A). In the lesioned area, the number of OLIG2+ cells per mm2 was counted and plotted. In the same region, the number of Plp1⁻ OLIG2⁺ cells per mm² was also calculated and plotted (see Figure 6B for Plp1⁺ cells per mm²). n=5-6 mice, mean+SEM, P: p-value, ** P<0.01, Student’s t test (unpaired, two tailed). scale bar: 100 µm. The white dotted line demarcates the lesion area.
