Distinct stage-specific transcriptional states of B cells derived from human tonsillar tissue.

B cells within secondary lymphoid tissues encompass a diversity of activation states and multiple maturation processes that reflect antigen recognition and transition through the germinal center (GC) reaction, in which mature B cells differentiate into memory and antibody-secreting cells (ASCs). Here, utilizing single-cell RNA-seq, we identify a range of distinct activation and maturation states of tonsil-derived B cells. In particular, we identify what we believe is a previously uncharacterized CCL4/CCL3 chemokine-expressing B cell population with an expression pattern consistent with B cell receptor/CD40 activation. Furthermore, we present a computational method that leverages regulatory network inference and pseudotemporal modeling to identify upstream transcription factor modulation along a GC-to-ASC axis of transcriptional maturation. Our data set provides valuable insight into diverse B cell functional profiles and will be a useful resource for further studies into the B cell immune compartment.

Comparison of busulfan and total body irradiation conditioning on hematopoietic clonal dynamics following lentiviral gene transfer in rhesus macaques.

The clonal dynamics following hematopoietic stem progenitor cell (HSPC) transplantation with busulfan conditioning are of great interest to the development of HSPC gene therapies. Compared with total body irradiation (TBI), busulfan is less toxic and more clinically relevant. We used a genetic barcoded HSPC autologous transplantation model to investigate the impact of busulfan conditioning on hematopoietic reconstitution in rhesus macaques. Two animals received lower busulfan dose and demonstrated lower vector marking levels compared with the third animal given a higher busulfan dose, despite similar busulfan pharmacokinetic analysis. We observed uni-lineage clonal engraftment at 1 month post-transplant, replaced by multilineage clones by 2 to 3 months in all animals. The initial multilineage clones in the first two animals were replaced by a second multilineage wave at 9 months; this clonal pattern disappeared at 13 months in the first animal, though was maintained in the second animal. The third animal maintained stable multilineage clones from 3 months to the most recent time point. In addition, busulfan animals exhibit more rapid HSPC clonal mixing across bone marrow sites and less CD16+ NK-biased clonal expansion compared with TBI animals. Therefore, busulfan conditioning regimens can variably impact the marrow niche, resulting in differences in clonal patterns with implications for HSPC gene therapies.

Immune aging in multiple sclerosis is characterized by abnormal CD4 T cell activation and increased frequencies of cytotoxic CD4 T cells with advancing age.

BACKGROUND: Immunosenescence (ISC) describes age-related changes in immune-system composition and function. Multiple sclerosis (MS) is a lifelong inflammatory condition involving effector and regulatory T-cell imbalance, yet little is known about T-cell ISC in MS. We examined age-associated changes in circulating T cells in MS compared to normal controls (NC). METHODS: Forty untreated MS (Mean Age 43·3, Range 18-72) and 49 NC (Mean Age 48·6, Range 20-84) without inflammatory conditions were included in cross-sectional design. T-cell subsets were phenotypically and functionally characterized using validated multiparametric flow cytometry. Their aging trajectories, and differences between MS and NC, were determined using linear mixed-effects models. FINDINGS: MS patients demonstrated early and persistent redistribution of naïve and memory CD4 T-cell compartments. While most CD4 and CD8 T-cell aging trajectories were similar between groups, MS patients exhibited abnormal age-associated increases of activated (HLA-DR+CD38+; (P = 0·013) and cytotoxic CD4 T cells, particularly in patients >60 (EOMES: P < 0·001). Aging MS patients also failed to upregulate CTLA-4 expression on both CD4 (P = 0·014) and CD8 (P = 0·009) T cells, coupled with abnormal age-associated increases in frequencies of B cells expressing costimulatory molecules. INTERPRETATION: While many aspects of T-cell aging in MS are conserved, the older MS patients harbour abnormally increased frequencies of CD4 T cells with activated and cytotoxic effector profiles. Age-related decreased expression of T-cell co-inhibitory receptor CTLA-4, and increased B-cell costimulatory molecule expression, may provide a mechanism that drives aberrant activation of effector CD4 T cells that have been implicated in progressive disease. FUNDING: Stated in Acknowledgements section of manuscript.

Cellular and humoral immune responses following SARS-CoV-2 mRNA vaccination in patients with multiple sclerosis on anti-CD20 therapy.

SARS-CoV-2 messenger RNA vaccination in healthy individuals generates immune protection against COVID-19. However, little is known about SARS-CoV-2 mRNA vaccine-induced responses in immunosuppressed patients. We investigated induction of antigen-specific antibody, B cell and T cell responses longitudinally in patients with multiple sclerosis (MS) on anti-CD20 antibody monotherapy (n = 20) compared with healthy controls (n = 10) after BNT162b2 or mRNA-1273 mRNA vaccination. Treatment with anti-CD20 monoclonal antibody (aCD20) significantly reduced spike-specific and receptor-binding domain (RBD)-specific antibody and memory B cell responses in most patients, an effect ameliorated with longer duration from last aCD20 treatment and extent of B cell reconstitution. By contrast, all patients with MS treated with aCD20 generated antigen-specific CD4 and CD8 T cell responses after vaccination. Treatment with aCD20 skewed responses, compromising circulating follicular helper T (TFH) cell responses and augmenting CD8 T cell induction, while preserving type 1 helper T (TH1) cell priming. Patients with MS treated with aCD20 lacking anti-RBD IgG had the most severe defect in circulating TFH responses and more robust CD8 T cell responses. These data define the nature of the SARS-CoV-2 vaccine-induced immune landscape in aCD20-treated patients and provide insights into coordinated mRNA vaccine-induced immune responses in humans. Our findings have implications for clinical decision-making and public health policy for immunosuppressed patients including those treated with aCD20.

Comparative engraftment and clonality of macaque HSPCs expanded on human umbilical vein endothelial cells versus non-expanded cells.

Ex vivo hematopoietic stem and progenitor cell (HSPC) expansion platforms are under active development, designed to increase HSPC numbers and thus engraftment ability of allogeneic cord blood grafts or autologous HSPCs for gene therapies. Murine and in vitro models have not correlated well with clinical outcomes of HSPC expansion, emphasizing the need for relevant pre-clinical models. Our rhesus macaque HSPC competitive autologous transplantation model utilizing genetically barcoded HSPC allows direct analysis of the relative short and long-term engraftment ability of lentivirally transduced HSPCs, along with additional critical characteristics such as HSPC clonal diversity and lineage bias. We investigated the impact of ex vivo expansion of macaque HSPCs on the engineered endothelial cell line (E-HUVECs) platform regarding safety, engraftment of transduced and E-HUVEC-expanded HSPC over time compared to non-expanded HSPC for up to 51 months post-transplantation, and both clonal diversity and lineage distribution of output from each engrafted cell source. Short and long-term engraftment were comparable for E-HUVEC expanded and the non-expanded HSPCs in both animals, despite extensive proliferation of CD34+ cells during 8 days of ex vivo culture for the E-HUVEC HSPCs, and optimization of harvesting and infusion of HSPCs co-cultured on E-HUVEC in the second animal. Long-term hematopoietic output from both E-HUVEC expanded and unexpanded HSPCs was highly polyclonal and multilineage. Overall, the comparable HSPC kinetics of macaques to humans, the ability to study post-transplant clonal patterns, and simultaneous multi-arm comparisons of grafts without the complication of interpreting allogeneic effects makes our model ideal to test ex vivo HSPC expansion platforms, particularly for gene therapy applications.

Interrogation of clonal tracking data using barcodetrackR.

Clonal tracking methods provide quantitative insights into the cellular output of genetically labelled progenitor cells across time and cellular compartments. In the context of gene and cell therapies, clonal tracking methods have enabled the tracking of progenitor cell output both in humans receiving therapies and in corresponding animal models, providing valuable insight into lineage reconstitution, clonal dynamics, and vector genotoxicity. However, the absence of a toolbox for analysis of clonal tracking data has precluded the development of standardized analytical frameworks within the field. Thus, we developed barcodetrackR, an R package and accompanying Shiny app containing diverse tools for the analysis and visualization of clonal tracking data. We demonstrate the utility of barcodetrackR in exploring longitudinal clonal patterns and lineage relationships in a number of clonal tracking studies of hematopoietic stem and progenitor cells (HSPCs) in humans receiving HSPC gene therapy and in animals receiving lentivirally transduced HSPC transplants or tumor cells.

Clonal tracking of erythropoiesis in rhesus macaques.

The classical model of hematopoietic hierarchies is being reconsidered on the basis of data from in vitro assays and single cell expression profiling. Recent experiments suggested that the erythroid lineage might differentiate directly from multipotent hematopoietic stem cells / progenitors or from a highly biased subpopulation of stem cells, rather than transiting through common myeloid progenitors or megakaryocyte-erythrocyte progenitors. We genetically barcoded autologous rhesus macaque stem and progenitor cells, allowing quantitative tracking of the in vivo clonal output of thousands of individual cells over time following transplantation. CD34+ cells were lentiviral-transduced with a high diversity barcode library, with the barcode in an expressed region of the provirus, allowing barcode retrieval from DNA or RNA, with each barcode representing an individual stem or progenitor cell clone. Barcode profiles from bone marrow CD45-CD71+ maturing nucleated red blood cells were compared with other lineages purified from the same bone marrow sample. There was very high correlation of barcode contributions between marrow nucleated red blood cells and other lineages, with the highest correlation between nucleated red blood cells and myeloid lineages, whether at earlier or later time points post transplantation, without obvious clonal contributions from highly erythroid-biased or restricted clones. A similar profile occurred even under stressors such as aging or erythropoietin stimulation. RNA barcode analysis on circulating mature red blood cells followed over long time periods demonstrated stable erythroid clonal contributions. Overall, in this nonhuman primate model with great relevance to human hematopoiesis, we documented continuous production of erythroid cells from multipotent, non-biased hematopoietic stem cell clones at steady-state or under stress.

Impact of CMV Infection on Natural Killer Cell Clonal Repertoire in CMV-Naïve Rhesus Macaques.

Recent functional, gene expression, and epigenetic studies have suggested the presence of a subset of mature natural killer (NK) cells responsible for maintaining NK cell memory. The lack of endogenous clonal markers in NK cells impedes understanding the genesis of these cell populations. In humans, primates, and mice, this phenotype and memory or adaptive functions have been strongly linked to cytomegalovirus or related herpes virus infections. We have used transplantation of lentivirally-barcoded autologous hematopoietic stem and progenitor cells (HSPC) to track clonal hematopoiesis in rhesus macaques and previously reported striking oligoclonal expansions of NK-biased barcoded clones within the CD56-CD16+ NK cell subpopulation, clonally distinct from ongoing output of myeloid, B cell, T cell, and CD56+16- NK cells from HSPC. These CD56-CD16+ NK cell clones segregate by expression of specific KIR surface receptors, suggesting clonal expansion in reaction to specific environmental stimuli. We have now used this model to investigate the impact of rhesus CMV(RhCMV) infection on NK clonal dynamics. Following transplantation, RhCMVneg rhesus macaques display less dominant and oligoclonal CD16+ NK cells biased clones compared to RhCMVpos animals, however these populations of cells are still clearly present. Upon RhCMV infection, CD16+ NK cells proliferate, followed by appearance of new groups of expanded NK clones and disappearance of clones present prior to RhCMV infection. A second superinfection with RhCMV resulted in rapid viral clearance without major change in the mature NK cell clonal landscape. Our findings suggest that RhCMV is not the sole driver of clonal expansion and peripheral maintenance of mature NK cells; however, infection of macaques with this herpesvirus does result in selective expansion and persistence of specific NK cell clones, providing further information relevant to adaptive NK cells and the development of NK cell therapies.

Aberrant Clonal Hematopoiesis following Lentiviral Vector Transduction of HSPCs in a Rhesus Macaque.

Lentiviral vectors (LVs) are used for delivery of genes into hematopoietic stem and progenitor cells (HSPCs) in clinical trials worldwide. LVs, in contrast to retroviral vectors, are not associated with insertion site-associated malignant clonal expansions and, thus, are considered safer. Here, however, we present a case of markedly abnormal dysplastic clonal hematopoiesis affecting the erythroid, myeloid, and megakaryocytic lineages in a rhesus macaque transplanted with HSPCs that were transduced with a LV containing a strong retroviral murine stem cell virus (MSCV) constitutive promoter-enhancer in the LTR. Nine insertions were mapped in the abnormal clone, resulting in overexpression and aberrant splicing of several genes of interest, including the cytokine stem cell factor and the transcription factor PLAG1. This case represents the first clear link between lentiviral insertion-induced clonal expansion and a clinically abnormal transformed phenotype following transduction of normal primate or human HSPCs, which is concerning, and suggests that strong constitutive promoters should not be included in LVs.

Barcoding of Macaque Hematopoietic Stem and Progenitor Cells: A Robust Platform to Assess Vector Genotoxicity.

Gene therapies using integrating retrovirus vectors to modify hematopoietic stem and progenitor cells have shown great promise for the treatment of immune system and hematologic diseases. However, activation of proto-oncogenes via insertional mutagenesis has resulted in the development of leukemia. We have utilized cellular bar coding to investigate the impact of different vector designs on the clonal behavior of hematopoietic stem and progenitor cells (HSPCs) during in vivo expansion, as a quantitative surrogate assay for genotoxicity in a non-human primate model with high relevance for human biology. We transplanted two rhesus macaques with autologous CD34+ HSPCs transduced with three lentiviral vectors containing different promoters and/or enhancers of a predicted range of genotoxicities, each containing a high-diversity barcode library that uniquely tags each individual transduced HSPC. Analysis of clonal output from thousands of individual HSPCs transduced with these barcoded vectors revealed sustained clonal diversity, with no progressive dominance of clones containing any of the three vectors for up to almost 3 years post-transplantation. Our data support a low genotoxic risk for lentivirus vectors in HSPCs, even those containing strong promoters and/or enhancers. Additionally, this flexible system can be used for the testing of future vector designs.

Clonal expansion and compartmentalized maintenance of rhesus macaque NK cell subsets.

Natural killer (NK) cells recognize and eliminate infected and malignant cells. Their life histories are poorly understood, particularly in humans, due to lack of informative models and endogenous clonal markers. Here, we apply transplantation of barcoded rhesus macaque hematopoietic cells to interrogate the landscape of NK cell production, expansion, and life histories at a clonal level long term and after proliferative challenge. We identify oligoclonal populations of rhesus CD56-CD16+ NK cells that are characterized by marked expansions and contractions over time yet remained long-term clonally uncoupled from other hematopoietic lineages, including CD56+CD16- NK cells. Individual or groups of CD56-CD16+ expanded clones segregated with surface expression of specific killer immunoglobulin-like receptors. These clonally distinct NK cell subpopulation patterns persisted for more than 4 years, including after transient in vivo anti-CD16-mediated depletion and subsequent regeneration. Profound and sustained interleukin-15-mediated depletion was required to generate new oligoclonal CD56-CD16+ NK cells. Together, our results indicate that linear NK cell production from multipotent hematopoietic progenitors or less mature CD56+CD16- cells is negligible during homeostasis and moderate proliferative stress. In such settings, peripheral compartmentalized self-renewal can maintain the composition of distinct, differentiated NK cell subpopulations.

The impact of aging on primate hematopoiesis as interrogated by clonal tracking.

Age-associated changes in hematopoietic stem and progenitor cells (HSPCs) have been carefully documented in mouse models but poorly characterized in primates and humans. To investigate clinically relevant aspects of hematopoietic aging, we compared the clonal output of thousands of genetically barcoded HSPCs in aged vs young macaques after autologous transplantation. Aged macaques showed delayed emergence of output from multipotent (MP) clones, with persistence of lineage-biased clones for many months after engraftment. In contrast to murine aging models reporting persistence of myeloid-biased HSPCs, aged macaques demonstrated persistent output from both B-cell and myeloid-biased clones. Clonal expansions of MP, myeloid-biased, and B-biased clones occurred in aged macaques, providing a potential model for human clonal hematopoiesis of indeterminate prognosis. These results suggest that long-term MP HSPC output is impaired in aged macaques, resulting in differences in the kinetics and lineage reconstitution patterns between young and aged primates in an autologous transplantation setting.

Geographic clonal tracking in macaques provides insights into HSPC migration and differentiation.

The geographic distribution of hematopoiesis at a clonal level is of interest in understanding how hematopoietic stem and progenitor cells (HSPCs) and their progeny interact with bone marrow (BM) niches during regeneration. We tagged rhesus macaque autologous HSPCs with genetic barcodes, allowing clonal tracking over time and space after transplantation. We found marked geographic segregation of CD34+ HSPCs for at least 6 mo posttransplantation, followed by very gradual clonal mixing at different BM sites over subsequent months to years. Clonal mapping was used to document local production of granulocytes, monocytes, B cells, and CD56+ natural killer (NK) cells. In contrast, CD16+CD56- NK cells were not produced in the BM, and in fact were clonally distinct from multipotent progenitors producing all other lineages. Most surprisingly, we documented local BM production of CD3+ T cells early after transplantation, using both clonal mapping and intravascular versus tissue-resident T cell staining, suggesting a thymus-independent T cell developmental pathway operating during BM regeneration, perhaps before thymic recovery.

Quantitative stability of hematopoietic stem and progenitor cell clonal output in rhesus macaques receiving transplants.

Autologous transplantation of hematopoietic stem and progenitor cells lentivirally labeled with unique oligonucleotide barcodes flanked by sequencing primer targets enables quantitative assessment of the self-renewal and differentiation patterns of these cells in a myeloablative rhesus macaque model. Compared with other approaches to clonal tracking, this approach is highly quantitative and reproducible. We documented stable multipotent long-term hematopoietic clonal output of monocytes, granulocytes, B cells, and T cells from a polyclonal pool of hematopoietic stem and progenitor cells in 4 macaques observed for up to 49 months posttransplantation. A broad range of clonal behaviors characterized by contribution level and biases toward certain cell types were extremely stable over time. Correlations between granulocyte and monocyte clonalities were greatest, followed by correlations between these cell types and B cells. We also detected quantitative expansion of T cell-biased clones consistent with an adaptive immune response. In contrast to recent data from a nonquantitative murine model, there was little evidence for clonal succession after initial hematopoietic reconstitution. These findings have important implications for human hematopoiesis, given the similarities between macaque and human physiologies.

Acquired somatic mutations in PNH reveal long-term maintenance of adaptive NK cells independent of HSPCs.

Natural killer (NK) cells have long been considered short-lived effectors of innate immunity. However, recent animal models and human studies suggest that subsets of NK cells have adaptive features. We investigate clonal relationships of various NK-cell subsets, including the adaptive population, by taking advantage of naturally occurring X-linked somatic PIGA mutations in hematopoietic stem and progenitor cells (HSPCs) from patients with paroxysmal nocturnal hemoglobinuria (PNH). The affected HSPCs and their progeny lack expression of glycosylphosphatidylinositol (GPI) anchors on their cell surface, allowing quantification of PIGA-mutant (GPI-negative) HSPC-derived peripheral blood cell populations. The fraction of GPI-negative cells within the CD56dim NK cells was markedly lower than that of neutrophils and the CD56bright NK-cell compartments. This discrepancy was most prominent within the adaptive CD56dim NK-cell population lacking PLZF expression. The functional properties of these adaptive NK cells were similar in PNH patients and healthy individuals. Our findings support the existence of a long-lived, adaptive NK-cell population maintained independently from GPIposCD56dim.

SAMHD1 transcript upregulation during SIV infection of the central nervous system does not associate with reduced viral load.

Restriction of HIV-1 in myeloid-lineage cells is attributed in part to the nucleotidase activity of the SAM-domain and HD-domain containing protein (SAMHD1), which depletes free nucleotides, blocking reverse transcription. In the same cells, the Vpx protein of HIV-2 and most SIVs counteracts SAMHD1. Both Type I and II interferons may stimulate SAMHD1 transcription. The contributions of SAMHD1 to retroviral restriction in the central nervous system (CNS) have been the subject of limited study. We hypothesized that SAMHD1 would respond to interferon in the SIV-infected CNS but would not control virus due to SIV Vpx. Accordingly, we investigated SAMHD1 transcript abundance and association with the Type I interferon response in an SIV model. SAMHD1 transcript levels were IFN responsive, increasing during acute phase infection and decreasing during a more quiescent phase, but generally remaining elevated at all post-infection time points. In vitro, SAMHD1 transcript was abundant in macaque astrocytes and further induced by Type I interferon, while IFN produced a weaker response in the more permissive environment of the macrophage. We cannot rule out a contribution of SAMHD1 to retroviral restriction in relatively non-permissive CNS cell types. We encourage additional research in this area, particularly in the context of HIV-1 infection.