Modeling the cell-type-specific mesoscale murine connectome with anterograde tracing experiments.

The Allen Mouse Brain Connectivity Atlas consists of anterograde tracing experiments targeting diverse structures and classes of projecting neurons. Beyond regional anterograde tracing done in C57BL/6 wild-type mice, a large fraction of experiments are performed using transgenic Cre-lines. This allows access to cell-class-specific whole-brain connectivity information, with class defined by the transgenic lines. However, even though the number of experiments is large, it does not come close to covering all existing cell classes in every area where they exist. Here, we study how much we can fill in these gaps and estimate the cell-class-specific connectivity function given the simplifying assumptions that nearby voxels have smoothly varying projections, but that these projection tensors can change sharply depending on the region and class of the projecting cells. This paper describes the conversion of Cre-line tracer experiments into class-specific connectivity matrices representing the connection strengths between source and target structures. We introduce and validate a novel statistical model for creation of connectivity matrices. We extend the Nadaraya-Watson kernel learning method that we previously used to fill in spatial gaps to also fill in gaps in cell-class connectivity information. To do this, we construct a "cell-class space" based on class-specific averaged regionalized projections and combine smoothing in 3D space as well as in this abstract space to share information between similar neuron classes. Using this method, we construct a set of connectivity matrices using multiple levels of resolution at which discontinuities in connectivity are assumed. We show that the connectivities obtained from this model display expected cell-type- and structure-specific connectivities. We also show that the wild-type connectivity matrix can be factored using a sparse set of factors, and analyze the informativeness of this latent variable model.

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.

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.

Clonal tracking of rhesus macaque hematopoiesis highlights a distinct lineage origin for natural killer cells.

Analysis of hematopoietic stem cell function in nonhuman primates provides insights that are relevant for human biology and therapeutic strategies. In this study, we applied quantitative genetic barcoding to track the clonal output of transplanted autologous rhesus macaque hematopoietic stem and progenitor cells over a time period of up to 9.5 months. We found that unilineage short-term progenitors reconstituted myeloid and lymphoid lineages at 1 month but were supplanted over time by multilineage clones, initially myeloid restricted, then myeloid-B clones, and then stable myeloid-B-T multilineage, long-term repopulating clones. Surprisingly, reconstitution of the natural killer (NK) cell lineage, and particularly the major CD16(+)/CD56(-) peripheral blood NK compartment, showed limited clonal overlap with T, B, or myeloid lineages, and therefore appears to be ontologically distinct. Thus, in addition to providing insights into clonal behavior over time, our analysis suggests an unexpected paradigm for the relationship between NK cells and other hematopoietic lineages in primates.