Dbx2, an Aging-Related Homeobox Gene, Inhibits the Proliferation of Adult Neural Progenitors

In the adult mouse brain, the subventricular zone (SVZ) underlying the lateral ventricles harbours a population of quiescent neural stem cells, which can be activated (aNSCs) to initiate proliferation and generate a neurogenic lineage consisting of transit amplifying progenitors (TAPs), neuroblasts (NBs) and newborn neurons. This process is markedly reduced during aging. Recent studies suggest that the aged SVZ niche decreases the pool of proliferating neural/stem progenitor cells (NSPCs), and hence adult neurogenesis, by causing transcriptomic changes that promote NSC quiescence. The transcription factors that mediate these changes, however, remain unclear. We previously found that the homeobox gene Dbx2 is upregulated in NSPCs of the aged mouse SVZ and can inhibit the growth of NSPC cultures. Here, we further investigate its role as a candidate transcriptional regulator of neurogenic decline. We show that Dbx2 expression is downregulated by Epidermal Growth Factor receptor signaling, which promotes NSPC proliferation and decreases in the aged SVZ. By means of transgenic NSPC lines overexpressing Dbx2, we also show that this gene inhibits NSPC proliferation by hindering the G2/M transition. Furthermore, we exploit RNA sequencing of transgenic NSPCs to elucidate the transcriptomic networks modulated by Dbx2. Among the top hits, we report the downregulation of the molecular pathways implicated in cell cycle progression. Accordingly, we find that Dbx2 function is negatively correlated with the transcriptional signatures of proliferative NSPCs (aNSCs, TAPs and early NBs). These results point to Dbx2 as a transcription factor relaying the anti-neurogenic input of the aged niche to the NSPC transcriptome. Graphical Abstract Supplementary Information The online version contains supplementary material available at 10.1007/s12015-023-10600-7.


Introduction
In adult mice, the subventricular zone (SVZ) of the lateral ventricles and the subgranular zone (SGZ) of the dentate gyrus harbour neurogenic neural stem/progenitor cells (NSPCs) within a specialized niche.The niche milieu includes other cell types, such as ependymal, vascular, glial and neuronal cells, as well as locally produced cues and systemic factors, acting together to support NSPC maintenance and neuronal differentiation [1].SVZ and SGZ NSPCs are organized in a hierarchical cell lineage encompassing the neurogenic process [2,3].Neural stem cells (NSCs) lie at the top of the lineage and are largely quiescent (qNSCs).Upon activation (aNSCs) by appropriate stimuli, they can enter the cell cycle, to self-renew or to generate daughter cells that give rise to transient amplifying progenitors (TAPs), neuroblasts (NBs) and newborn neurons [4,5].
Andrea Giuliani and Valerio Licursi contributed equally to this work.
An age-related reduction of adult neurogenesis might be a contributing factor to the cognitive decline of elderly humans [6,7].This is best documented in mice, wherein proliferating NSPCs and their neuronal output are markedly decreased in the aged SVZ and SGZ [8,9].Two explanations have been proposed for this reduction.Earlier models suggest that the age-dependent decrease in neurogenesis is driven by the progressive exhaustion of the qNSC pool due to their activation and differentiation [10,11].According to recent studies, however, a qNSC reservoir persists in the aged brain [12][13][14][15].Moreover, some NSCs may return to quiescence after activation and self-renewal [16][17][18].These observations lead to alternative models, whereby the neurogenic decline may be caused by the increased quiescence of the aged NSC pool.Considering the heterogeneous nature of NSC populations [12,14,15], and the time-dependent changes in neurogenesis [16,19], both models may be represented in certain NSC subpopulations and/or ages.
Transcriptomic comparison of qNSCs and aNSCs has revealed clear differences in the expression of genes related to cell cycle, protein synthesis and inflammatory response [16,18,22].Transcriptomic changes have also been detected between young adult and aged NSPCs [12,14,15,[28][29][30].This suggests that transcriptional regulation in NSPCs is key to control their proliferation and to the alterations of aged NSPCs.It is therefore important to identify the transcription factors that relay aged niche signals to the NSPC transcriptome and uncover how they function in this context.
The homeobox gene Dbx2 encodes a transcription factor involved in the development of spinal cord interneurons, and is also expressed in selected regions of the embryonic brain, limb buds and tooth germs [31,32].We previously reported that Dbx2 is expressed in the adult mouse SVZ, and established that Dbx2 is upregulated in aged NSPCs and inhibits the growth of NSPC cultures [30]; however, the extracellular signals controlling Dbx2 expression and the mechanisms mediating Dbx2 function remained unclear.Here, we show that Dbx2 transcription is negatively modulated by EGFR signaling.We provide evidence that Dbx2 can hinder the G2/M transition in adult NSPCs, thus reducing their proliferation.Finally, we describe the Dbx2-regulated transcriptome in adult NSPC cultures.Notably, this expression program includes cell cycle and inflammatory response gene networks, and is negatively correlated with the signatures of the proliferative neurogenic cell types (aNSCs, TAPs and early NBs).

Real-time qRT-PCR Assays
For the analysis of Dbx2 mRNA levels by real-time qRT-PCR, young adult SVZ NSPCs were seeded in T25 flasks or 6-well plates at a density of approximately 40000 cells/cm 2 in adherent proliferative culture conditions.On the next day, cultures were rinsed once with medium devoid of EGF and FGF2, then incubated in fresh medium with both EGF and FGF2, none of them, or only one of them.Treatments with human recombinant Transforming Growth Factor α (TGFα, Cell Guidance Systems GFH39-5) were performed by diluting a 50 µg/ml stock in PBS containing 0.1% Bovine Serum Albumin with culture medium to 100 ng/ml.Treatments with the EGFR inhibitor AG1478 (Calbiochem 658548) were performed by diluting a 10 mM stock in DMSO with culture medium to 2 µM; control cultures were treated with equal volumes of DMSO.At the end of the experimental treatment, cells cultured in T25 flasks were enzymatically detached with Accutase (Corning 25-058-CI), pelleted in 2 ml tubes and lysed for RNA purification, whereas cells in 6-well plates were directly lysed in the culture plates.Total RNA purification and quantification, reverse transcription and real-time qPCR were performed using Qiagen kits, a NanoDrop 2000 (Thermo Scientific) and a Rotor-Gene Q (Qiagen), as previously described [33].The following primers were used for Dbx2 and for the reference gene Rpl19: Dbx2_forward CCC GCC ATT CTA CTC TGC AT, Dbx2_ reverse GAG TCC TGG GTC AGC AAA GG, Rpl19_forward AGA CCA AGG AAG CAC GAA AG, Rpl19_reverse GCC GCT ATG TAC AGA CAC GA.

Cell Proliferation Analyses by Flow Cytometry and Immunocytochemistry
For cell cycle analysis, transgenic NSPCs were seeded in T25 flasks at a density of 5000 cells/cm 2 in non-adherent culture conditions supplemented with 600 ng/ml of puromycin.Neurospheres generated by control or Dbx2-overexpressing NSPCs were pelleted by centrifugation and dissociated with Accutase.The resulting cell suspensions were rinsed with PBS, fixed in cold PBS-methanol 1:1 and stored in this solution at 4 °C.Staining with propidium iodide (PI) and quantification of the fraction of cells in G0/G1, S and G2/M by flow cytometry were performed as previously described [35].
To quantify the fraction of cells expressing Ki67 and phosphorylated histone H3 (pH3), control and Dbx2-overexpressing neurospheres were collected by centrifugation and dissociated, followed by cell plating on poly-ornithine/ laminin coated glass coverslips in 24-well plates and culture for 2-4 hours (2-4 h) in media for adherent proliferative conditions to allow attachment to the glass.After attachment, cells were fixed for 20 min with methanol-free formaldehyde (Pierce 28908) diluted to 4% with PBS, then rinsed a few times with PBS and stored in PBS at 4 °C.Immunofluorescence analysis with an anti-Ki67 mouse monoclonal antibody (Leica Biosystems ACK02, 1:100) or with an anti-pH3 (Ser 10) rabbit polyclonal antibody (ThermoFisher PA5-17869, 1:100) was performed as previously described [30,35].In each experiment, 12 random fields per coverslip were selected based on Hoechst staining and photographed using a Nikon Eclipse TE300 microscope and a Nikon DS-U3 digital camera, followed by the quantification of the ratio between Ki67-positive (Ki67+) or pH3-positive (pH3+) cells and the total cell number in 1-3 coverslips for each experimental condition.Several hundred to a few thousand cells were counted for each condition in each experiment.Automated immunofluorescence analysis was performed with ImageJ software [36].Colocalization images were generated using the ImageJ plug-in "Colocalization Image Creator" [37].For Hoechst signal, the "Binary Element" option was used, and the automatic thresholding was set on "YEN"; for Ki67 and pH3 signals, the "Grayscale Element" option was chosen.The colocalization images generated were then used for automated counting of Ki67+ or pH3+ cells using the ImageJ plug-in "Colocalization Object Counter" [37].Only Ki67 and pH3 signals effectively colocalizing with Hoechst signal were counted.Automated analysis of signal intensity was performed using the ImageJ tool "Analyse Particles".A nuclear mask was generated using Hoechst signal to specifically select nuclear regions.This mask was then merged with the Ki67 or pH3 signals and the fluorescence intensity per pixel in each nucleus was measured with the ImageJ tool "Measure".

Statistical Analysis of Gene Expression and Cell Proliferation Assays
The experimental data obtained from qRT-PCR, flow cytometry and immunofluorescence assays were analysed and graphically represented using GraphPad Prism 9 software.The type of statistical test performed, the resulting p values and the number of independent experimental replicates performed with NSPCs from different culture passages (n) are indicated in the figure legends.

RNA-sequencing and Differential Gene Expression Analysis
For RNA-sequencing (RNA-seq), total RNA was purified from frozen cell pellets of control and Dbx2-overexpressing non-adherent cultures using Qiagen RNeasy kits.RNA-seq libraries from total RNA (100 ng) from each sample were prepared using QuantSeq 3' mRNA-Seq Library prep kit (Lexogen, Vienna, Austria), according to manufacturer's instructions, at Next Generation Diagnostics (Pozzuoli, Italy).The amplified fragmented cDNAs of 300 bp in size were sequenced in single-end mode using the NextSeq500 (Illumina) with a read length of 101 bp.Sequence read quality was evaluated using FastQC version 0.11.8 (Babraham Institute, Cambridge, UK), followed by trimming using bbduk software to remove adapter sequences, poly-A tails and low-quality end bases (Q < 20).Reads were then mapped to the mouse Ensembl GRCm38 (mm10) build reference genome with STAR version 2.5.0a[38]; gene annotations corresponding to Ensembl annotation release 96 were used to build a transcriptome index, which was provided to STAR during the alignment.
To identify differentially expressed genes (DEGs), data were filtered to remove from the analysis the genes having > 1 counts per million in less than 4 out of 16 total samples.Data normalization and differential gene expression analysis were performed using Bioconductor, R package edgeR version 3.36 [39], assigning the cell line batch as a covariate to fit the generalized log-linear model with the glmQLFit function of edgeR.DEGs were assessed by comparing Dbx2-NSPC and GFP-NSPC samples using a moderated t-test with a false discovery rate (FDR) < 0.05 threshold.Volcano plot was created using Bioconductor [40,41], R package EnhancedVolcano version 1.12.0.Principal component analysis (PCA) was performed using R base prcomp function, plotting the first two PCs using R package ggplot2 version 3.3.5.

Gene Ontology, Gene Set and Transcription Factor Enrichment Analyses
DEGs in Dbx2-NSPCs vs GFP-NSPCs were clustered by functional annotation in gene ontology (GO) and pathway enrichment analysis using Bioconductor R package cluster-Profiler version 4.2.0 [42] with annotation of GO database [43], and with annotation of REACTOME [44] and Kyoto Encyclopedia of Genes and Genomes (KEGG) [45] for pathways.R package pathview version 1.34.0 was used to integrate the RNA-seq data with KEGG pathway plots.
Broad Institute gene set enrichment analysis (GSEA) [46] was used to assess the enrichment of the Dbx2-associated signature list (12533 expressed genes in Dbx2-NSPCs and GFP-NSPCs without FDR threshold, ranked according to fold change (FC), versus the curated "Hallmark" and "C2" gene set collections from the Broad Molecular Signatures Database (MSigDB) version 7.4.1 (https:// www.gseamsigdb.org/ gsea/ msigdb/), or versus recently published mouse SVZ single-cell RNA-seq datasets [14,15,47].The mouse version of the gene set collections "Hallmark" and "C2" (REACTOME subcategory, selecting signatures related to cell cycle) were obtained from R package msigdbr version 7.4.1.GSEA enrichment score (ES) was calculated by walking down the ranked list of genes, increasing a running-sum statistic when a gene was in the gene set and decreasing it when it was not.A normalized ES (NES) was also calculated with GSEA, by considering differences in pathway size (i.e., gene set size) and allowing for comparisons between pathways within the analysis.
Transcription factor enrichment analysis was performed with Enrichr method [48], using as references TRANSFAC and JASPAR library [49], containing manually curated transcription factors binding profiles as position weight matrices, TRRUST library [50], a database of reference transcription factor-target regulatory interactions, and ENCODE and CHEA library, containing experimental transcription factor binding data from ENCODE Project Consortium [51] and CHEA database [52].

Dbx2 Expression is Upregulated in Aged or Quiescent NSPCs and is Repressed by EGFR Signaling
We previously reported that Dbx2 mRNA levels were increased both in NSPC cultures and in freshly isolated NSPCs obtained from the SVZ of aged (18 mo) mice as compared to adult (3-7 mo) samples [30].To confirm these results, we checked a recent transcriptomic dataset of NSPCs isolated from the SVZ of 2 mo, 7 mo and 19 mo mice, noticing an age-dependent upregulation of Dbx2 [14] (Fig. 1A).We also examined published transcriptomic datasets of qNSCs and aNSCs isolated from the adult mouse SVZ [15,18,22], consistently observing Dbx2 upregulation in qNSCs (Fig. 1B to D).
We then investigated whether EGFR signaling may be involved in Dbx2 expression changes upon NSC activation or NSPC aging.To this goal, we analysed Dbx2 mRNA levels in cultures of young adult (3 mo) mouse SVZ NSPCs that were maintained for 6 h in standard proliferation-supporting media, supplemented with EGF and FGF2, in media devoid of both growth factors, or in media containing only EGF or FGF2.As shown in Fig. 1E, Dbx2 was upregulated in NSPC cultures devoid of exogenous EGF, when compared with EGF-supplemented cultures.Furthermore, Dbx2 was downregulated when NSPCs were cultured without growth factors for 6 h and then treated with EGF for 2 h, as compared to growth factor withdrawal for 8 h (Fig. 1F).TGFα is the predominant EGFR ligand in the adult brain and its levels decrease in the aged SVZ [21].We observed a comparable decrease in Dbx2 mRNA levels in NSPC cultures treated for 6 h with either EGF or TGFα, in comparison with NSPCs cultured without exogenous growth factors; this decrease was prevented by including AG1478, an EGFR inhibitor, in the culture media (Fig. 1G).These results suggest that EGFR signaling represses Dbx2 expression in SVZ NSPCs; consequently, a reduction of EGFR activation during quiescence or aging may lead to Dbx2 upregulation in NSPCs (Fig. 1H).

Dbx2 Overexpression in Young Adult Mouse SVZ NSPCs Inhibits Cell Proliferation by Hindering the G2/M Transition
Dbx2 overexpression reduces the growth of NSPC cultures derived from the SVZ of young adult mice, without major effects on NSPC viability, suggesting that Dbx2 inhibits NSPC proliferation [30].We speculated that Dbx2 may promote cell cycle exit and entry in G0, or decrease the division rate of NSPCs remaining in the cell cycle.To test this hypotheses, we took advantage of previously described transgenic young adult NSPC lines with constitutive expression of either a mouse Dbx2 transgene (Dbx2-NSPCs) or a control GFP transgene (GFP-NSPCs), and cultured them for 4-6 days (4-6 d) in non-adherent conditions to generate neurospheres [30].We then used these cultures to quantify the percentage of Dbx2-NSPCs and GFP-NSPCs that were positive for Ki67, a nuclear marker expressed throughout the cell cycle, but not in cells that have exited the cell cycle.Moreover, we estimated the fraction of Dbx2-NSPCs and GFP-NSPCs in each phase of the cell cycle (G0/G1, S, G2/M) by flow cytometry quantification of DNA content in cells stained with PI.To make sure that the results were reproducible, we performed these analyses with two independent pairs of Dbx2-NSPC and GFP-NSPC lines [30].Cell cycle exit of NSPCs in response to Dbx2 overexpression would be expected to cause an increase in the fraction of cells in G0/G1 and a decrease of the fraction of Ki67+ cells; however, these effects were not detectable in the Dbx2-NSPC lines as compared to GFP-NSPC cultures (Fig. 2A, D, E, H; Fig. S1A, D, E, H; Fig. S2A and  D; representative images of Ki67-stained cultures are shown in Fig. 2I and J, and in Fig. S1I and J).Surprisingly, we detected a reproducible increase of the G2/M cell fraction upon Dbx2 overexpression (Fig. 2C, D, H; Fig. S1C, D, H).To distinguish between cells in G2 and in M phases, we quantified the fraction of cells positive for pH3, a marker of cells undergoing mitosis.Mitotic histone H3 phosphorylation starts in pericentric heterochromatin in late G2, then spreads throughout the condensing chromatin during prophase, reaching a peak at metaphase; dephosphorylation of histone H3 begins in anaphase and ends in telophase [53].In agreement with the dynamics of mitotic pH3 accumulation, we observed different patterns of pH3 immunostaining during NSPC progression through mitosis.Some of the pH3+ NSPCs, likely corresponding to late G2 cells, showed a weaker staining in the form of isolated nuclear spots; during prophase, these pH3+ foci increased in size and intensity, until the whole nucleus appeared to be stained (Fig. S3A to D).A peak in signal intensity was observed in NSPC nuclei between metaphase and early anaphase, followed by a progressive decrease between anaphase and telophase (Fig. S3I to L).We quantified  (I, M, K, O) and Dbx2-NSPC (J, N, L, P) cultures stained with anti-Ki67 (I, J) or anti-pH3 (K, L) antibodies.Hoechst nuclear staining is shown in (M to P). Scale bar, 40 µm the total percentage of pH3+ NSPCs, which included all the cells engaging in the G2/M transition, as well as the percentage of NSPCs showing pH3 staining patterns typical of cells that have progressed beyond early prophase (Fig. S3D and I to L), which we named pH3+ late cells.When compared with GFP-NSPCs, a decrease in the pH3+ total cell fraction was detected in one of the Dbx2-NSPC lines, but not the other; however, both lines showed a reduction in the pH3+ late cell fraction (Fig. 2F and G; Fig. S1F and G; Fig. S2B and E; representative images of pH3stained cultures are shown in Fig. 2K and L, and in Fig S1K and  L).The fraction of pH3+ cells within the top 40% of the fluorescence intensity range, which likely include pH3 + late cells, was also reduced in both Dbx2-NSPC lines (Fig. S1C and F).These results indicate that the accumulation in G2/M observed by flow cytometry is due to an altered mitotic progression in Dbx2-NSPC cultures.Therefore, they lead us to suggest that Dbx2 negatively regulates the G2/M transition in NSPCs, hindering the progression through mitosis when its expression levels increase.The observation that one of the Dbx2-NSPC lines showed a weaker effect on the G2/M transition than the other is consistent with the comparatively weaker growth phenotype of the neurosphere cultures generated from this cell line [30], and might be explained by differences in the genetic or epigenetic background of the two cell lines, as suggested by PCA results (see below).Nonetheless, the number of cells undergoing cell division is decreased in both lines, which provides an explanation for the reduced size of the neurospheres formed by Dbx2-NSPCs.

Dbx2 Regulates NSPC Transcriptional Programmes Associated with Cell Cycle, Inflammatory Response and Lineage Progression
To gain insight into the molecular mechanisms mediating the effects of Dbx2 overexpression in SVZ NSPCs, we performed a transcriptomic analysis with the transgenic NSPC lines used for cell cycle assays.To this aim, we performed RNA-seq and differential gene expression analysis with RNA samples obtained from 4-6 d non-adherent cultures of both lines of Dbx2-NSPCs and GFP NSPCs (4 independent experiments for each pair).PCA confirmed that the cell culture samples used for transcriptomic analysis were mainly clustered according to transgene expression (Dbx2 or GFP) (PC1, 63% of variance), although samples belonging to the same experimental condition, but to independent cell lines, could also be distinguished (PC2, 21% of variance) (Fig. 3C).By setting a threshold of FDR < 0.05, we identified 2854 DEGs between Dbx2-NSPCs and GFP-NSCPs (1368 upregulated and 1486 downregulated genes in Dbx2-NSPCs) (Fig. 3A and Table S1).
We then performed GO enrichment analysis of the DEGs using the GO Biological Process database.As shown in Fig. 3B, the genes downregulated upon Dbx2 overexpression were enriched for GO terms related to NSPC proliferation (e.g."neural precursor cell proliferation", "cell cycle phase transition"), lineage progression (e.g."cell fate commitment", "regulation of neuron differentiation"), and RNA processing (e.g."regulation of RNA splicing"); the genes upregulated in Dbx2-overexpressing NSPCs were enriched for GO terms related to inflammatory response (e.g."inflammatory response", "cytokine production").DEG analysis using REAC-TOME and KEGG databases also indicated that pathways related to cell cycle and RNA processing were enriched in the downregulated genes, and pathways related to inflammatory response and cell cycle inhibition ("p53 signaling pathway") were enriched in the upregulated genes (Fig. 3D and Fig. S4).
To characterise these transcriptional effects in more detail, we analysed the transcriptome of Dbx2-overexpressing NSPCs by GSEA, a computational method that investigates the biological effects associated with a ranked gene list by determining whether the members of a gene set linked to a specific biological process tend to occur at the top or the bottom of the ranked gene list [46].Since the size of the analysed gene list influences GSEA accuracy, we performed this analysis with a wider Dbx2-associated signature of 12533 genes, which included all the genes with detectable mRNA expression in Dbx2-NSPCs and GFP-NSPCs, irrespectively of FDR values.This gene list was ranked according to FC, followed by GSEA with "Hallmark" gene sets in MSigDB.Notably, several cell cycle-related gene sets ("E2F targets", "G2M checkpoint", "Myc targets") correlated with the Dbx2associated signature (Fig. 4A), with negative NES values resulting from the enrichment of these gene sets among the genes downregulated in Dbx2-overexpressing NSPCs (Fig. 4B, C and D; Table S2).Moreover, several gene sets linked to inflammatory response (e.g."Interferon gamma response", "IL6 JAK STAT3 signaling") or to cell cycle inhibition ("p53 pathway") correlated with the Dbx2-associated signature (Fig. 4A), with positive NES values resulting from the enrichment of these gene sets among the genes upregulated in Dbx2-overexpressing NSPCs (Fig. 4E, F and G, Table S2).Transcription factor motif enrichment analysis confirmed an enrichment of E2F and Myc motifs in the DEGs downregulated in Dbx2-overexpressing NSPCs, and an enrichment of p53 motifs in the upregulated DEGs (Fig. S5).
Given Dbx2 effects on NSPC proliferation, we performed a more detailed analysis of the transcriptional changes related to cell cycle progression in Dbx2-overexpressing NSPCs.To this aim, we carried out GSEA using the Dbx2associated signature and various gene sets related to different cell cycle phases as available in the "C2" curated gene sets collection in MSigDB.This analysis revealed a general negative correlation of the Dbx2-associated signature with cell cycle-related gene sets (Fig. S6).A characterization of the DEGs between Dbx2-NSPCs and GFP-NSCPs (FDR < 0.05) with KEGG cell cycle and p53 signaling pathway gene sets showed that several genes crucially implicated in the G1/S transition (Cdk2, Cdk4) or G2/M transition (CycA, Cdk1) pressing NSPCs is indicated in the X axis (down, FC < 0; up, FC > 0).The heat maps at the bottom of the plots highlight the proportion of Dbx2-downregulated genes (blue colour) and Dbx2-upregulated genes (red colour) in each gene set ◂ were downregulated in Dbx2-overexpressing NSPCs, and key cell cycle inhibitors acting in the p53 pathway (p21, Gadd45) were upregulated (Fig. S7).These results suggest that elevated Dbx2 expression broadly inhibits the transcriptional pathways involved in cell cycle progression.
The results of GSEA with cell cycle-related gene sets prompted us to use the same approach to investigate the effects of Dbx2 overexpression on the transcriptional programmes associated with NSPC lineage progression.To this aim, we took advantage of the recently described single-cell transcriptomic signatures of qNSCs, aNSCs, TAPs and NBs of the adult mouse SVZ [14,15,47]; we then used the gene sets related to each of these cell populations to perform GSEA with the Dbx2-associated signature.Remarkably, the gene sets related to aNSCs, TAPs and NBs negatively correlated with the Dbx2-associated signature (Fig. 5, Table S3).Although the signature related to the whole qNSC population was not correlated with the Dbx2-associated signature (Fig. 5A), a positive correlation (NES > 0) was observed when the analysis was carried out using the gene sets related to specific qNSC subpopulations (Fig. 5C and E).Since aNSCs, TAPs and NBs encompass the proliferating NSPC pool of the SVZ neurogenic lineage, we wondered whether the negative correlation that we observed between the gene sets related to these cell populations and the Dbx2-associated signature may be driven by cell cycle-related genes.To address this question, we quantified the percentage of REACTOME and KEGG cell cycle genes among the aNSC-related and TAP-related gene sets.This quantification suggested that cell cycle genes alone did not explain the correlation between the aNSC/TAP-related and the Dbx2-associated signatures; for example, cell cycle genes accounted for only 26-28% of the genes in the signatures related to proliferative TAPs (mTAPs and L1) (Fig. 5B  and D).We then filtered out cell cycle genes from aNSC/ TAP signatures, selected the Dbx2-downregulated DEGs (FDR < 0.05) among the remaining genes, and performed GO enrichment analysis on the resulting gene lists.Notably, GO terms related to RNA processing (e.g."RNA splicing") prominently featured in all the analysed gene lists (Fig. 6).These data suggest that Dbx2 may inhibit NSPC proliferation by regulating cell cycle-related genes, and NSPC lineage progression by modulating genes involved in RNA processing.

Discussion
In this study, we show that elevated Dbx2 expression in NSPCs impairs the G2/M transition, causing accumulation in G2/M by hindering the progression through mitosis.Although we could not detect a clear effect of Dbx2 overexpression on the G1/S transition, transcriptomic analyses revealed that Dbx2 broadly regulates the cell cycle molecular networks, repressing genes implicated in cell cycle progression and activating genes involved in cell cycle arrest.Of note, signatures related to E2F, Myc and p53 were among the top hits, which is consistent with the recently described association between these factors and the transcriptomic profiles of aged NSPCs [15].
Recent studies suggest that the increase in the time spent by NSCs in a quiescent state, which is generally thought to be in G0, is a main driver of the neurogenic decline of aged mice [12][13][14][15][16].Although Drosophila qNSCs can also be arrested in G2 [54], evidence that this may happen in mice is lacking.Transcriptomic comparison of qNSCs and aNSCs of the mouse SVZ, however, showed that both G1/S and G2/M signatures are downregulated in qNSCs [18]; thus, the inhibitory effects of Dbx2 on the transcriptional networks associated with both the G1/S and G2/M transitions is compatible with a role in the increased quiescence of aged NSCs.Notably, the gene set downregulated upon Dbx2 overexpression correlates with the transcriptomic signatures of the proliferative populations of the SVZ neurogenic lineage (aNSCs, TAPs and early NBs), supporting the implication of Dbx2 in NSC quiescence from a different angle.The limited percentage of cell cycle genes contributing to this correlation suggests that it is not driven predominantly by cell cycle inhibition; furthermore, GO analysis indicates an additional negative effect of Dbx2 on genes related to RNA processing.We speculate that increased Dbx2 expression during SVZ aging may facilitate NSC progressing into quiescence by inhibiting cell cycle genes along with RNA processing genes.
Several categories associated with the inflammatory response were top hits among the genes upregulated by Dbx2 overexpression.This is noteworthy, since inflammatory transcriptional signatures are associated with NSC quiescence, aged neurogenic niches and brain aging [8,9,18,55].We did not detect an increase in Dbx2 expression levels in young adult NSPC cultures treated with IFNαγ, IL2/6 or TNFα (our unpublished observations); however, gene sets related to immune response (especially IFN signaling) were upregulated in SVZ NSPC cultures treated with BMP4 and FGF2 to induce quiescence [25].Thus, these signatures may be activated in NSPCs independently of canonical cytokine-driven responses to regulate their proliferative state.
In conclusion, our study uncovers several links between increased Dbx2 function and key features of NSPC aging.These include: i) the downregulation of gene networks implicated in cell cycle progression, the upregulation of those involved in cell cycle arrest and the inhibition of NSPC proliferation; ii) the upregulation of inflammatory response pathways; iii) the downregulation of the transcriptomic signatures of the proliferative neurogenic cell populations (aNSCs, TAPs and early NBs), which are depleted in the aged SVZ; iv) the negative regulation of Dbx2 expression by EGFR signaling, which is reduced in the aged SVZ.This work has some limitations, since our functional experiments have been performed using in vitro NPSC cultures and gain-of-function approaches.Although some of the genes modulated upon Dbx2 overexpression showed opposite changes in aged NSPC cultures expressing Dbx2-targeting shRNAs, this approach yielded a modest Dbx2 knockdown in our hands [30].Of note, Dbx2 is a NSC-specific component of the recently reported molecular aging clocks in the mouse SVZ, and its age-associated modulation can be reversed by exercise [56], pointing to Dbx2 as a potential biomarker of neurogenic niche rejuvenation.showing that Dbx2 can promote age-associated changes in NSPC cultures, our in vitro analyses prompt the investigation of the functional role of Dbx2 in the aging of in vivo neurogenic niches.

Fig. 1
Fig. 1 Activation and EGFR signaling downregulate Dbx2 expression in adult NSPCs.(A) Dbx2 mRNA levels in freshly isolated NSPCs from the SVZ of 2 mo, 7 mo and 19 mo mice as reported in a published RNA-seq dataset [14].Black lines show mean transcript levels ± standard deviation; blue and red dots show individual experimental replicates for young adult and aged mice, respectively (n = 2); **, adjusted p value < 0.01, 2 mo vs 19 mo NSPCs.(B) Dbx2 mRNA levels in freshly isolated qNSCs (red dots) and aNSCs (blue dots) from the SVZ of adult mice as reported in a published microarray dataset [22] (n = 3); *, adjusted p value < 0.05.(C) Box-and-whisker plots of Dbx2 mRNA levels in freshly isolated qNSCs and aNSCs from the SVZ of adult mice as reported in a published RNA-seq dataset [18].The lower and higher whiskers indicate the minimum and maximum values, respectively.The bottom and top of the box represent the first and third quartiles, respectively, and the band inside the box indicates the second quartile (the median).The + symbol indicates mean transcript levels.Red and blue triangles show individual experimental replicates for qNSCs and aNSCs, respectively (n = 4); ***, p < 0.001, Student's t-test.(D) Box-and-whisker plots of Dbx2 mRNA levels in freshly isolated GFP high (GFP-H, red triangles) and GFP low (GFP-L, blue triangles) NSCs (largely corresponding to qNSCs and aNSCs, respectively) from the SVZ of adult mice as

Fig. 3
Fig. 3 Transcriptomic profiling of Dbx2-overexpressing NSPCs.(A) Volcano plot showing the distribution of -log 10 (FDR) values (Y axis) relative to log 2 (FC) values (X axis) resulting from the comparison of gene expression levels in Dbx2-NSPCs vs GFP-NSPCs.Red dots indicate DEGs associated with FDR < 0.05 and FC < -2 or FC > 2; blue dots indicate genes with FDR < 0.05 and -2 < FC < 2; green dots indicate genes with FC < -2 or FC > 2, and FDR > 0.05.Representative genes are highlighted in the plot.(B) Dot plot showing the GO terms enriched in the upregulated (left) and the downregulated (right) DEGs.The size of the dots is based on the count of the DEGs cor-

Fig. 4
Fig. 4 Dbx2 overexpression affects the expression of gene sets related to cell cycle and inflammatory response.(A) Bar plot showing the top "Hallmark" gene sets correlated with the Dbx2-associated transcriptional signature (12533 genes), according to GSEA.The size and the colour of the bars indicate the NES value and the qvalue associated with each gene set, respectively.(B to G) Enrichment plots showing the top 3 NES < 0 gene sets (B to D), or the top 3 NES > 0 gene sets (E to G), in the Dbx2-associated signature.Vertical black lines indicate individual members of each gene set and their position within the ranked Dbx2-associated signature; the FC direction in Dbx2-overexpressing NSPCs is indicated in the X axis (down, FC < 0; up, FC > 0).The heat maps at the bottom of the plots highlight the proportion of Dbx2-downregulated genes (blue colour) and Dbx2-upregulated genes (red colour) in each gene set

Fig. 5 Fig. 6
Fig.5Dbx2 function is negatively correlated with neurogenic lineage progression.(A, C and E) Bar plots showing the correlation between the Dbx2associated signature and the qNSC-specific, aNSC-specific, TAP-specific and NB-specific gene sets, as defined in recent studies [[47] (A),[15] (C),[14] (E)], according to GSEA.In (C), H0-H3 correspond to qNSC subpopulations, L0 to aNSCs, L1 to TAPs and L2 to NBs, as previously described[15].(B, D and F) Enrichment plots showing selected cell type-specific gene sets in the Dbx2-associated signature.The percentages indicate the fraction of cell cycle genes in each gene set