ARID3a expression in human hematopoietic stem cells is associated with distinct gene patterns in aged individuals

Background Immunologic aging leads to immune dysfunction, significantly reducing the quality of life of the elderly. Aged-related defects in early hematopoiesis result in reduced lymphoid cell development, functionally defective mature immune cells, and poor protective responses to vaccines and pathogens. Despite considerable progress understanding the underlying causes of decreased immunity in the elderly, the mechanisms by which these occur are still poorly understood. The DNA-binding protein ARID3a is expressed in a subset of human hematopoietic progenitors. Inhibition of ARID3a in bulk human cord blood CD34+ hematopoietic progenitors led to developmental skewing toward myeloid lineage at the expense of lymphoid lineage cells in vitro. Effects of ARID3a expression in adult-derived hematopoietic stem cells (HSCs) have not been analyzed, nor has ARID3a expression been assessed in relationship to age. We hypothesized that decreases in ARID3a could explain some of the defects observed in aging. Results Our data reveal decreased frequencies of ARID3a-expressing peripheral blood HSCs from aged healthy individuals compared with young donor HSCs. Inhibition of ARID3a in young donor-derived HSCs limits B lineage potential, suggesting a role for ARID3a in B lymphopoiesis in bone marrow-derived HSCs. Increasing ARID3a levels of HSCs from aged donors in vitro alters B lineage development and maturation. Finally, single cell analyses of ARID3a-expressing HSCs from young versus aged donors identify a number of differentially expressed genes in aged ARID3A-expressing cells versus young ARID3A-expressing HSCs, as well as between ARID3A-expressing and non-expressing cells in both young and aged donor HSCs. Conclusions These data suggest that ARID3a-expressing HSCs from aged individuals differ at both molecular and functional levels compared to ARID3a-expressing HSCs from young individuals.


Background
The US Census Bureau estimates that nearly a quarter of the US population will be over the age of 65 by the year 2060 [1]. A major consequence of aging is a decline in immune function. Both murine and human studies revealed age-related defects in early hematopoietic development, and functional defects in mature immune cell populations, that result in decreased potentials to mount protective immune responses in aged individuals (reviewed in [2]), as exemplified by increased susceptibility to influenza and pneumonia in the elderly. Human hematopoiesis is a dynamic process requiring complex regulation of multiple gene expression pathways for lineage commitment and resulting in development of diverse numbers of blood cell types [3]. Several studies indicate that aged hematopoietic stem cell (HSC) frequencies are increased in both human and mouse [4,5]. Bulk CD34-expressing hematopoietic stem and progenitor cells (HSPCs) from aged donors exhibit epigenetic and transcriptional changes that promote self-renewal over differentiation [6,7]. In old age, cells in the hematopoietic progenitor pool accumulate decreased telomere lengths and DNA damage markers, and their developmental potential becomes increasingly skewed toward myeloid versus lymphoid lineage development [6][7][8][9]. Identification of the changes in old age that alter the development of mature immune cells, and possibly contribute to their dysfunction, will require mechanistic studies that better define potential differences in gene regulatory mechanisms critical for lineage choices.
The transcription factor ARID3a is an understudied member of a large family of proteins with epigenetic functions [10][11][12]. Previous studies indicated that ARID3a can contribute to repression and enhancement of transcription in a cell type-specific fashion [13][14][15]. Regulation of transcription through ARID3a may be associated with epigenetic functions that affect large subsets of genes and lineage decisions ( [13,[16][17][18], our unpublished data). ARID3a deficient mice die in utero between days 12 and 14 of gestation, but exhibited a 90% reduction in HSCs numbers in the fetal liver associated with defective erythropoiesis and B lymphopoiesis [19]. Functional loss of ARID3a in B lineage cells, either through ARID3a dominant negative transgenic mice or rare adult ARID3a −/− mice, revealed important roles for ARID3a in B1 B cell lineage development and function [19,20]. These data were recently confirmed using conditional knockout mice, showing definitively that ARID3a is required for B1 B lineage development in mice [21][22][23]. Loss of functional ARID3a in B lineage cells in mice directly impaired normal protective immune responses to infection with S. pneumoniae [20], an organism associated with pneumonia in aged individuals [24,25]. Forced expression of ARID3a in mouse B lineage cells resulted in enhanced development of B1 and MZ B cells versus conventional follicular B cells [26], suggesting ARID3a levels can modulate B lineage responses in mice. Mechanisms responsible for generating B1 lineage B cells in man remain controversial [27,28]. Together, these data identify ARID3a as an important regulator of B lymphopoiesis.
Roles for ARID3a in human hematopoiesis are less clear. We found that ARID3a is variably expressed in healthy human HSPCs, including total CD34 + HSPCs, HSCs, multipotent progenitor (MPP), multi-lymphoid progenitors (MLP), and multi-myeloid progenitors (MMP) derived from adult peripheral blood [29], but the functional significance of expression in those progenitors is not clear. In functional studies with human cord blood HSPCs, where ARID3a expression dominates the majority of those cells, manipulation of ARID3a resulted in skewing of lineage development with promotion of myeloid over lymphoid lineage differentiation upon loss of ARID3a expression and increased B lymphopoiesis upon over-expression of ARID3a [30]. ARID3a expression in circulating peripheral blood HSPCs from lupus erythematosus patients is upregulated compared to similar cells from healthy individuals, although the role of ARID3a in those cells is unknown [29]. These data suggest the need for further experiments to determine how ARID3a levels affect adult human hematopoiesis.
We hypothesized that one explanation for reduced B lymphopoiesis and increased numbers of myeloid cells in aged versus young individuals is that ARID3a expression is reduced in HSCs from healthy aged individuals compared to healthy young individuals, or that its function in those cells is impaired. Our results indicate that peripheral blood HSCs from aged donors exhibit reduced frequencies of ARID3a-expressing cells compared with young donors. Furthermore, modulation of ARID3a levels in both aged and young donor-derived HSCs altered B lymphopoiesis in vitro. Finally, single cell RNA-seq analyses revealed unexpected differences in gene expression patterns in ARID3A-expressing progenitors from aged versus young individuals.

Results
We assessed numbers of ARID3a-expressing HSCs (lin − CD34 + CD38 − CD45RA − CD49f + ) from the peripheral blood of 54 healthy adults of varying ages by flow cytometry (Fig. 1a), as previously described [29,30]. Total HSC numbers increased with age (p < 0.0001) (Fig. 1b), consistent with reports by others that aged individuals (red dots) have increased numbers of HSCs [4-6, 31, 32]. Further analyses of those cells for ARID3a expression, revealed decreases in frequencies of ARID3a + HSCs with increasing donor age (≤10% ARID3a + HSCs), in donors over the age of 50 (Fig. 1c). Based on these data, for our purposes, we defined "young" as under the age of 50 and "aged" as 50 years of age and over but show donors over 65 as red symbols. The frequency of ARID3a + HSCs in aged donors was reduced by an average of 80% compared to young donors (Fig. 1d), but total ARID3a + HSC numbers were not statistically different between young and aged donors (Fig. 1e). Similarly, the amounts of ARID3a within HSCs varied among individuals, but were not statistically different between young and aged individuals (Fig. 1f). Data from individuals over 65 (red symbols) did not differ from individuals over 50 in any of these characteristics. These data show average Fig. 1 PBMCs from aged individuals show reduced frequencies of ARID3a + HSCs. Human HSCs expressing ARID3a were identified by flow cytometry from total PBMCs of 54 healthy adults (ages 18-70). a Representative gating from young (left) and aged (center) donors, with an isotype control for ARID3a (right) are shown. b Numbers of HSCs/ml of blood are plotted against age with a line of best fit. Red samples indicate samples from donors ≥65 years of age. c Frequencies of ARID3a + HSCs are plotted against age. Frequencies (d) and numbers (e) of ARID3a + HSCs for young (triangles) and aged (diamonds) donors are shown. f Mean Fluorescence Intensities (MFI) of ARID3a + HSCs from young and aged donors are presented. Averages and standard error bars are shown. Statistical significance was determined by Mann Whitney U test and Pearson's Correlation absolute numbers of ARID3a-expressing HSCs are maintained with age, but that relative frequencies of ARID3a + HSCs in the peripheral blood are decreased with increased chronologic age.
Although numbers of HSCs increase with chronologic age (Fig. 1b), there was considerable variability in numbers of HSCs in aged individuals and in some young individuals. We considered that this variability might reflect variations in biological age versus chronological age. Telomere length has been used as a biologic measure of cellular age [33,34], and telomeres in hematopoietic cells from aged persons have been reported to be shortened compared to those from younger individuals [8]. Therefore, we chose individual HSC samples of both young and aged donor HSCs in Fig. 1a for assessment of telomere lengths (indicated in Fig. 2a). Aged samples with reduced numbers of HSCs (closed triangles) had longer telomeres than those samples with higher numbers of HSCs (closed boxes) (Fig. 2b) (Fig. 2b). Others reported an average of 18 Kb for telomere lengths in hematopoietic progenitors from young individuals [8]. Measurements of telomere lengths from other individuals of various ages (circles, Fig. 2c), showed good correlations between telomere length and chronologic age. Neither the frequency of ARID3a + cells (Fig. 2d), nor the total numbers of ARID3a + cells per sample (Fig. 2e) in the HSCs for which we performed telomere analyses were associated with telomere length. Therefore, no correlation was observed between ARID3a status and telomere length.
Another molecular marker associated with aging is accumulation of DNA damage as indicated by accumulation of H2A.xp139 (H2Ax) [9,35]. This assay allows direct comparison of ARID3a levels with the aging marker via flow cytometric analyses (Fig. 2f). As expected, aged donor HSCs displayed an increased frequency of H2Ax as compared to HSCs from young donors (Fig. 2g). There was also a significant positive correlation between age and H2Ax expression (Fig. 2h), and a negative correlation between frequencies of ARID3a and H2Ax (Fig. 2i). Total frequencies of ARID3a-expressing cells with H2Ax marks associated with DNA damage were not different in HSCs derived from young versus aged individuals (Fig. 2j). However, increased numbers of H2Ax-marked HSCs without ARID3a occurred in aged individuals compared to young individuals (Fig. 2k), suggesting the possibility that increases in total HSCs with age may occur as an expansion of cells that do not express ARID3a. These data suggest that ARID3a + HSCs from aged donors do not reflect cells that have accumulated increases in DNA damage to a greater degree than in young individuals.
We hypothesized that ARID3a-expressing cells from aged individuals could still be functionally deficient in generating B lymphoid responses that are dependent on ARID3a [20]. Therefore, we explored the developmental potential of aged versus young donor derived HSCs using in vitro cultures that allowed development of B lymphoid, natural killer (NK) and myeloid lineage cells [36]. Because ARID3a is an intracellular protein, it is not possible to isolate those cells via flow cytometry for use in vitro. After 3 weeks in culture, both young and aged donor HSCs gave rise to the same total numbers of cells (Fig. 3a). In these in vitro cultures, total percentages of CD33-expressing myeloid cells were equivalent in cultures initiated with HSCs from young or aged donors (Fig. 3b). HSCs from aged donors generated reduced frequencies of B lineage cells compared to those from young donors (Fig. 3c). Consistent with previous studies [37], aged HSCs generated increased percentages of NK cells compared to young HSCs (Fig. 3d). All HSCs originally express CD34 and lose that surface marker as they differentiate into various hematopoietic lineages [38][39][40] . Therefore, as a means of investigating maturation of cells in these cultures, we assessed numbers of cells of each lineage that retained CD34 expression (Fig. 3e). Total frequencies of CD34 + cells were greatly increased in cultures derived from aged versus young donor HSCs (Fig.  3f). In the aged donor HSCs, the majority of both CD33 and CD19-expressing cells also maintained CD34 expression, compared to lower percentages of monocyte and B lymphoid cells derived from young donor HSCs (Fig. 3f). Conversely, frequencies of CD34-expressing NK lineage cells were reduced in aged HSC donor derived cultures compared to those obtained from young HSC donors (Fig. 3f). These data confirm previous results demonstrating reduced B lineage development from aged donor HSCs [6], and extend those studies to suggest that the B lymphocyte and myeloid lineage cells derived from aged donor HSCs may be less mature than those derived from young donors as reflected by maintenance of CD34 expression.
To determine if modulation of ARID3a levels affected the developmental potential of HSCs, we first forced expression of ARID3a in aged donor HSCs. ARID3a expression was inferred from co-expression of RFP in the transfecting vector, which revealed > 80% RFPexpressing cells (Fig. 4a, b). ARID3a-transfected cultures had increased total cell numbers compared to both untransfected, media only and RFP-only vector controls (Fig. 4c). CD33-expressing myeloid lineage cells were decreased in frequency after ARID3a expression, but total cell numbers remained equivalent (Fig. 4d). As predicted, ARID3a over-expression in aged donor HSCs led Telomere lengths of HSCs from samples above and below (b), or on/near the trend line (c) are plotted against age. Telomere lengths are plotted against frequency (d), and total numbers of ARID3a + HSCs (e). f A representative flow cytometric plot of HSCs analyzed for ARID3a and H2A.xp139 is shown. g Frequencies of H2A.x + HSCs from young (n = 5) and aged (n = 7) donors are shown. h Frequencies of H2A.x + HSCs are plotted against age. i Frequencies of H2A.x + HSCs are plotted against frequencies of ARID3a + HSCs. Frequencies of H2A + ARID3a + HSCs (j), and H2A + ARID3a − HSCs (k) are plotted for young and old donors. Averages and standard error bars are indicated. Statistical significance was determined by Mann Whitney U test and Pearson's Correlation. **, p = 0.0051 to large increases in both B lineage cell frequencies and cell numbers (Fig. 4e), indicating that increased ARID3a expression in adult-derived HSCs from aged donors drove B lymphocyte development in vitro. NK frequencies did not differ among the three culture conditions, but total numbers of NK cells were increased on average in the ARID3a-expressing cultures (Fig. 4f). Examination of CD34 expression in these cultures revealed only slight differences in CD34 expression between ARID3a transfected cultures and the controls except in the CD19 + B lineage cells where CD34 expressing cells were significantly reduced in frequency compared to the controls ( Fig. 4g and h). These data reveal that increased ARID3a expression in aged donor Similarly, inhibition of ARID3a expression in young donor HSCs using shRNA vectors that co-expressed GFP revealed > 80% GFP-expressing cells, suggesting that ARID3a was reduced in those cells as well (Fig. 5a,  b). No differences were observed in total cell numbers between controls and ARID3a-inhibited cultures were observed (Fig. 5c), nor were differences observed in CD33expressing monocyte or NK lineage cells between vector control and ARID3a-inhibited cultures ( Fig. 5d and f). However, frequencies and total cell numbers of B lineage cells without ARID3a were dramatically reduced compared with media and control shRNA cultures (Fig. 5e). While total frequencies of CD34expressing cells did not differ appreciably among the three culture conditions (Fig. 5g), the majority of CD19 + B lineage cells also co-expressed CD34 (Fig. 5h). These data demonstrate that ARID3a is required for B lineage maturation of young donor-derived HSCs in vitro. Myeloid versus lymphoid ratios in both the aged and young cultures with ARID3a modulation (Table 1) suggest that ARID3a is important for B lineage maturation in both aged and young HSCs. However, these data also suggest that aged HSCs may exhibit defects in function despite ARID3a expression.
To identify genes associated with ARID3a expression in HSCs from aged versus young donors, we performed single cell RNA-seq analyses and analyzed each cell for the presence and level of ARID3A transcription as shown by the scatter plot (Fig. 6a). Analyses of ARID3A transcript by qPCR from bulk HSCs of known ARID3a protein expression suggest that transcript and protein expression in bulk HSCs correlate (data not shown). There were 153 ARID3a + and 148 ARID3a − cells from aged donors and 172 ARID3a + and 92 ARID3a − cells from the young donors. Three-dimensional t distributed stochastic neighbor embedding plots (tSNE) of 301 aged and 264 young HSCs from 8 donors revealed considerable spread in dimensionality in the aged (circles) versus young (squares), shown as overlays (Fig. 6b and c). This suggests that isolation of HSCs using standard surface markers (Fig. 1a) results in cells that are heterogeneous with respect to their transcriptomes in both aged and young donors. Identification of ARID3A-expressing cells, as indicated by blue symbols, reveals widespread ARID3A expression in both aged and young HSCs, with increased clustering of ARID3A-expressing cells toward the righthand side in the aged donor cells (Fig. 6b). Ingenuity Pathway Analyses (IPA) analyses of ARID3A-associated genes and non-ARID3A associated genes from aged donors revealed enrichment in pathways associated with cell cycle, regulation of B cell apoptosis, negative regulation of B cell activation, and positive regulation of histone methylation in the ARID3A cells (Fig. 6d). Similar analyses of young donor cells indicated enrichment of pathways associated with lymphocyte homeostasis, JAK-STAT signaling and nucleic acid binding in the ARID3A cells (Fig. 6e).
Direct comparison of ARID3A-expressing cells from aged versus young donors revealed 253 genes downregulated > 2-fold and 195 up-regulated genes. IPA of the ARID3A-associated genes from the young donors and the ARID3A-associated genes in the aged donors revealed enrichment in pathways associated with B1 B cell differentiation lymphocyte differentiation, type 1 interferon signaling, and regulation of gene expression in the aged donor cells (Fig. 6f). The most highly upregulated gene in the aged ARID3A + HSCs versus young donor ARID3A + HSCS is AIF1, a gene associated with macrophage activation (Fig. 6g), while the most highly down-regulated gene encodes a chemokine, PPBP. In addition, at least four transcription factors HES1, MAFB, ATF4 and JAK1 are differentially regulated in aged versus young ARID3A-expressing HSCs. Thus, alterations in gene expression profiles in young versus aged ARID3A-expressing cells may be the result of changes in these regulatory proteins. Together, ARID3A-expressing cells in aged donors differ in gene expression patterns from ARID3A-expressing young donor HSCs.

Discussion
The data presented here suggest that HSCs from both aged and young donors contain equivalent numbers of ARID3a-expressing cells, although total frequencies of ARID3a + cells were reduced in aged individuals. Analyses of hematopoietic lineage potential revealed that aged donor HSCs yielded reduced numbers of B cells compared to those derived from young donor HSCs, despite the presence of equivalent total numbers of ARID3a-expressing progenitors. These data suggest ARID3a expression alone is insufficient to cause ageassociated shifts in B cell generation. Surprisingly, the aged donor-derived B cells retained CD34 expression to a greater degree than B lymphoid lineage cells derived from young donors, implying the aged donor B lineage cells could be less mature than young donor derived B cells. Artificial over-expression of ARID3a in aged donor-derived HSCs suggested that increased ARID3a expression resulted in better B lineage development without retention of surface CD34. ARID3a expression was critical for B lineage development from young donor peripheral blood-derived HSCs, consistent with our previous findings using HSCs from cord blood [30]. Importantly, single cell RNA-seq of both aged and young donor HSCs expressing ARID3A revealed a number of transcriptome differences between these cells, suggesting that the ARID3a-expressing cells in aged individuals differ at the molecular level from ARID3a-expressing cells from young individuals. These differences may contribute to the functional differences observed in hematopoiesis in vitro. Together, these data suggest that ARID3a-expressing progenitor cells in aged individuals differ at the transcription and functional level from those present in young individuals. Previous studies on HSC-enriched HSPC populations from human report developmental skewing of aged cells toward myeloid lineage at the expense of lymphoid lineage cells [5,7,41]. While skewing toward the myeloid lineage in these in vitro cultures was minimal (Table 1), our studies differed from those previously published in several ways. We used isolated HSCs rather than the more heterogeneous HSPCs used by others. In addition, most other studies used myeloid colony forming assays in combination with liquid assays to enrich for myeloid development. Only one study used a liquid culture system to assess B cell and myeloid lineage development [7], though that study used a different feeder cell line and cytokine cocktail. Another possibility is that our aged precursors may have contained more pre-committed lineage progenitors for myeloid cells, as the frequencies of the ARID3a-expressing B cells were reduced in the aged HSCs. Thus, there may have been competition for B cell development due to pre-committed myeloid lineage cells. However, B lineage cytokines were present at equivalent levels in cultures of HSPCs from young donors that gave rise to more B lymphocytes. All of these factors are likely to contribute to the reduced myeloid/lymphoid ratio skewing we observed compared to previous reports.
We also observed increases in NK cells in the aged donor-derived cultures compared to the young donorderived cells. Intriguingly, NK cells develop from precursors shared by B lymphocytes [42], and our unpublished results suggest that NK cells transcribe ARID3a. Others reported increased numbers of NK cells in the peripheral blood of aged individuals [43]. Further experiments would be merited to determine the potential relationship between B lineage and NK cell development in aged HSCs. Together, these data demonstrate that there are clear differences in the developmental potential of young and aged HSCs.
The aged donor B cells that developed in vitro retained CD34 expression, suggesting they may be less mature than the B lineage cells from young donorderived cultures and highlighting at least one difference at the protein level between ARID3a-expressing cells from aged versus young donors. CD34 expression is reduced at the pre-B cell stage during normal B lymphopoiesis [44,45]. Our previous data suggest that ARID3a is important for B lineage development at all early and immature stages of differentiation [30]. In this study, we could not meaningfully assess for which types of B lineage progenitors were present due to the limited numbers of cells generated in these assays. Irrespective of the developmental stage of B lineage cells derived from unmanipulated aged donor HSCs, when ARID3a expression was dramatically increased in those HSCs, they generated increased numbers of B lineage cells that no longer expressed CD34 and were presumably more mature. However, we cannot rule out the possibility that the cells that over-expressed ARID3a in those cultures were derived from some different subset of cells that did not have ARID3a originally. For example, myeloid progenitors, which are increased in old age, may have been diverted to the B lineage by over-expression of ARID3a (reviewed in [46]). Together, these data suggest that increased ARID3a expression in old age could improve B lineage development and perhaps allow increases in humoral immune responses in aged individuals.
In support of our surface expression data suggesting that the ARID3a-expressing cells from aged individuals differ from ARID3a + HSCs in the young, single cell transcriptome analyses revealed major differences in gene expression as well. GO analyses of young ARID3Aexpressing HSCs and aged ARID3A-expressing HSCs reveal an enrichment in the pathways associated with B1 B cells differentiation in the aged cells (Fig. 6d). Several elegant studies from the Hardy group in the last few years have elucidated the role of ARID3a in B1 B cell development in fetal/neonatal murine development [21][22][23]. While human B1 cell identification remains controversial, others reported that cord blood HSPCs and mouse fetal liver HSPCs are very similar phenotypically and functionally [47,48]. We previously found that~75% of human umbilical cord blood HSCs express ARID3a where it is important for B lineage development [30]. These data suggest the interesting possibility that the ARID3a-expressing HSCs from aged donors could be enriched for B1 lineage-like development and might have fewer precursors that would ultimately result in conventional B lymphocytes. Murine bone marrow derived B1 B cells are suggested to be continually generated into old age [49], a population that has been reported to have separate progenitors than conventional B cells [41]. Alternatively, age-associated B cells (ABCs) are a unique heterogeneous subset of B cells that have been the focus of study for several laboratories in the last decade. First described separately by both the Cancro and Marrack research groups, these cells expand in old age and have been attributed to origins from multiple types of B lineage cells, including B1 B cells [50][51][52]. However, no one has demonstrated the origins of those cells from early hematopoietic progenitors. Our data suggest that alterations in humoral immunity attributed to alterations in B cell maturation in aged individuals are already evident at the gene level in early hematopoietic stem cells. The differences in gene expression that we observed between young and aged donor HSCs are generally consistent with published reports from others using young and aged HSPCs [35,53]. However, as shown in Fig. 6, considerable heterogeneity remains at the transcriptome level even in purified HSCs. The ability to directly compare populations with the B cell-associated protein ARID3a underscores additional differences in aged versus young donor-derived HSCs, and present a new caveat that even within the ARID3a-expressing cells, some progenitors may be predisposed to develop into B1 versus conventional B cells as suggested by the GO analyses. More studies are clearly needed to better define the discrete progenitors within HSCs and the associated transcriptome differences between young and aged individuals.

Conclusions
The data presented here demonstrate that ARID3aexpressing HSCs total numbers are equivalent between old and young donors, yet the aged donor HSCs produced fewer B cells in vitro. Functional defects in aged donor ARID3a-expressing HSCs were reflected both in expression of surface protein markers and at the transcriptome level, suggesting either that the ARID3a-expressing HSCs from aged donors are functionally defective in maturation responses, or that the HSCs expressing ARID3a in aged donors reflect skewing of a precursor subset of B lineage cells with different capacities to expand in these in vitro cultures. In addition, these data highlight the expansion of HSCs in aged individuals that do not express ARID3a, raising the possibility that these cells contribute to ageassociated defects in hematopoiesis. Together, these data suggest that age-related defects in B cell function and humoral antibody responses may reflect changes that occur as early as the hematopoietic stem cell stage.

Cells
Peripheral blood mononuclear cells (PBMCs) from 54 healthy adults (ranging in age from 18 to 70) were isolated from fresh heparinized peripheral blood, or from leukocyte reduction chambers, with Ficoll-Pague Plus (GE Healthcare) or Lymphoprep (StemCell Technologies). Leukocyte reduction chambers typically contained 1-2 × 10 9 PBMCs. Four million total PBMCs were stained for flow cytometric analyses and the remaining cells were either immediately enriched for CD34 + hematopoietic progenitors using the EasySep Human CD34 Positive Selection Kit (StemCell Technologies) or were cryopreserved for later enrichment. Thawed CD34 + enriched HSPCs were cultured overnight in Stem Pro media (Invitrogen) supplemented with 100 U penicillin and 100 μg streptomycin and a cytokine cocktail including stem cell factor (SCF) (100 ng/ml), FLT3 ligand (50 ng/ml), and thrombopoietin (TPO) (10 ng/ml) (R&D Systems) to acclimate cells prior to sorting. CD34enriched cells were then sorted using a FACS Jazz (BD Biosciences) for CD38 − CD45RA − CD49f + HSCs for in vitro experiments [54]. Post-sort analyses were typically > 98% HSCs.

Flow Cytometry
Cell surface markers were identified with the following fluorochromes: human hematopoietic lineage markers (CD2, CD3, CD14, CD16, CD19, CD56, CD235a) conjugated to APC (eBioscience), CD34 PE, CD38 Alexa Fluor 700, CD49f BV421, and CD45RA Brilliant Violet 570 (all from BioLegend). Appropriate isotype controls (BioLegend) were used for gating HSCs as described [29,30]. In some cases, following surface marker staining, cells were fixed with fixation buffer (BD Biosciences), permeabilized with Transcription Factor Fixation/Permeabilization Buffer kit (eBioscience) and stained with goat anti-human ARID3a, followed by rabbit anti-goat FITC secondary antibody (Invitrogen) to identify ARID3a-expressing cells and/ or with H2A.X phopspho (Ser139) APC-Fire750 (Biologend). Human anti-ARID3a peptide specific antibodies were generated and isolated over a peptide-specific column as described [55]. Cells were pretreated with anti-human CD32 antibody to block Fc receptor binding prior to surface staining to exclude non-specific binding. Doublet exclusion was used to ensure analyses of single cells prior to forward/side scatter gating. Data were collected using an LSRII (BD Biogenics) and FACSDiva (BD Biosciences) software version 4.1 or Stratedigm S1200Ex and CellCapTure acquisition software and were analyzed using FlowJo (Tree Star) software version 10.

Telomere length analyses
Telomere lengths of sorted HSCs (5 young, 7 aged) were measured using the Relative Human Telomere Length Quantification qPCR assay kit (ScienCell Research Laboratories) according to the manufacturer's protocol. Samples were chosen based on the trend line from estimated numbers of HSCs per mL of peripheral blood. Briefly, analyses were performed in duplicate for each sample, with control reference samples with predetermined telomere lengths provided by the manufacturer. Average telomere lengths were determined and the reference samples with a known telomere length was used to calculate the absolute telomere length of each individual HSC population.

ARID3a modulation
FACS-purified HSCs (approximately 2000 cells per samples total divided between 3 treatments) were cultured in X-VIVO 10 medium (BioWhittaker) supplemented with 1% BSA, L-glutamine (Gibco), 100 U/mL penicillinstreptomycin (Invitrogen) and a cocktail of cytokines to allow hematopoietic lineage differentiation [36]: 100 ng/ ml SCF (BioLegend), 100 ng/ml FLT3L, 15 ng/ml TPO, 10 ng/ml G-CSF (all R&D Systems) and 10 ng/ml IL-6 (Invitrogen) for transduction for 24 h in the presence of 4 μg/ml polybrene, as previously described [30]. Lentivirus containing shRNA specific for ARID3a or an unrelated scramble control shRNA with a co-expressed GFP reporter were prepared by Genecopoeia, Inc., as described previously [30]. Over-expression of full length human ARID3a was achieved using Lentiviral vectors or empty control vectors co-expressing an RFP reporter and prepared by ABM Inc. HSCs from young donors were treated with shRNA lentivirus using a multiplicity of infection (MOI) of 1. HSCs from aged donors were treated similarly with lentivirus expressing control and full-length ARID3a with RFP similarly, at an MOI of 1. Following transduction, HSCs were harvested, washed, counted and seeded at 100 cells per well in triplicate for each treatment (polybrene only, control vectors and experimental vectors) onto 1 × 10 4 MS-5 murine stromal cells per well. Seeding numbers were chosen to ensure sufficient differentiation down multiple lineages and based on published data [56]. Stromal cells were plated 24 h previously in 96 well plates in H5100 medium (StemCell Technologies) supplemented with P/S, Lglutamine, and the following cytokines [36]: 100 ng/mL SCF, 10 ng/mL IL-2, 20 ng/mL IL-7 (BioLegend), 20 ng/ mL IL-6 (Invitrogen), 10 ng/mL Flt3-L, 50 ng/mL TPO, 20 ng/mL G-CSF, and 20 ng/mL GM-CSF (R&D Systems) [36]. Cultures were fed by 50% volume change once a week. Cells were harvested after 3 weeks and assessed by flow cytometric analyses for myeloid lineage (CD33-PE-Cy5), NK lineage (CD56-BV421), and B lineage (CD19-BV650) cells. MS-5 cells were excluded from analyses using expression of mouse CD106-APC, as previously described [30]. Due to the small numbers of cells initially transfected, pre-culture assessments of transfection efficiencies could not be made at the initiation of culture, but were assessed at the conclusion of the cultures via flow cytometry for the presence of either GFP or RFP reporter expression, as these markers were co-expressed with the ARID3a shRNA or full length ARID3a, respectively.

Single-cell RNA-seq
HSCs from 4 young (ages 19-40) and 4 aged (ages 61-70) individuals (2 males and 2 females each) were FACSsorted and immediately used for single cell analyses. Smart-seq/C1 libraries were prepared on the OUHSC Consolidated Core Laboratory Fluidigm C1 system using the SMARTer Ultra Low RNA Kit (Clontech) according to the manufacturer's protocol. Cells were loaded on a 5-10 μm RNA-seq microfluidic IFC at a concentration of 200,000/ml. Capture site occupancy was surveyed using a standard light microscope and recorded to verify cell capture. Sequencing library amplification was performed using Nextera XT Index primers (Illumina) according to manufacturer's protocol. Barcoded library concentration and fragment size distribution was determined using Agilent High Sensitivity D1000 kit on an Agilent 2200 TapeStation (Agilent Technologies) at the OMRF Genomics Core Facility. Paired-end (2 × 50 bp) sequencing was performed on a NovaSeq 6000 platform by the OMRF Genomics Core Facility.
Paired-end reads were aligned to the hg38 genome assembly with STAR using default parameters. All sequencing was performed on the same day to minimize batch effects from individual donors and they were further reduced using ComBat (Johnson et al., Biostatistic vol 8 issue 1, pgs 118-127, 2007) within the Partek Genomics suite on the gene expression matrices for young and old donors. The Partek Genomics Suite was used to process the aligned reads to estimate gene expression levels. Cells that had greater than 25% of counts mapping to mitochondrial genes indicate stressed or dying cells and were excluded from analyses. The number of detected genes was 44,025 across 301 cells for aged samples and 39,417 genes across 264 cells for young samples. Log normalized counts per millions (CPM) were log 2 -transformed and used for tSNE visualization and differential expression analyses. Data from aged and young samples were interrogated separately based on ARID3a expression, with > 0.5 CPM denoting ARID3A + cells and < 0.5 CPM denoting ARID3A − cells for further analysis. ANOVA analysis was performed within the Partek Flow software to identify DEGs with 2-fold or greater changes in expression and a False Discovery Rate (FDR) < 0.05. IPA analyses were performed on significant DEGs. RNAseq data are publicly available through the GEO NCBI database under the accession number GSE138544.

Data analyses
Data were statistically evaluated using the non-parametric Mann Whitney U test or the non-parametric Wilcoxon paired test to compare distribution of variables between groups. Correlations were evaluated using Pearson's Correlations. Statistical analysis was performed with Prism (Graphpad) software version 8.2. Differential gene expression was analyzed using ANOVA within the Partek Genomics Suite. Gene ontology analyses on differentially expressed genes were analyzed within the Partek Flow Genomics Suite using the gene set enrichment tool. P values of less than 0.05 were considered significant.