High-throughput, sensitive quantification of repopulating hematopoietic stem cell clones

Longitudinal clonal tracking in humanized mice reveals sustained polyclonal repopulation of gene-modified human-HSPC despite vector integration bias

BACKGROUND

Current understanding of hematopoiesis is largely derived from mouse models that are physiologically distant from humans. Humanized mice provide the most physiologically relevant small animal model to study human diseases, most notably preclinical gene therapy studies. However, the clonal repopulation dynamics of human hematopoietic stem and progenitor cells (HSPC) in these animal models is only partially understood. Using a new clonal tracking methodology designed for small sample volumes, we aim to reveal the underlying clonal dynamics of human cell repopulation in a mouse environment.

METHODS

Humanized bone marrow-liver-thymus (hu-BLT) mice were generated by transplanting lentiviral vector-transduced human fetal liver HSPC (FL-HSPC) in NOD.Cg-PrkdcscidIl2rgtm1Wjl/SzJ (NSG) mice implanted with a piece of human fetal thymus. We developed a methodology to track vector integration sites (VIS) in a mere 25 µl of mouse blood for longitudinal and quantitative clonal analysis of human HSPC repopulation in mouse environment. We explored transcriptional and epigenetic features of human HSPC for possible VIS bias.

RESULTS

A total of 897 HSPC clones were longitudinally tracked in hu-BLT mice-providing a first-ever demonstration of clonal dynamics and coordinated expansion of therapeutic and control vector-modified human cell populations simultaneously repopulating in the same humanized mice. The polyclonal repopulation stabilized at 19 weeks post-transplant and the contribution of the largest clone doubled within 4 weeks. Moreover, 550 (~ 60%) clones persisted over 6 weeks and were highly shared between different organs. The normal clonal profiles confirmed the safety of our gene therapy vectors. Multi-omics analysis of human FL-HSPC revealed that 54% of vector integrations in repopulating clones occurred within ± 1 kb of H3K36me3-enriched regions.

CONCLUSIONS

Human repopulation in mice is polyclonal and stabilizes more rapidly than that previously observed in humans. VIS preference for H3K36me3 has no apparent negative effects on HSPC repopulation. Our study provides a methodology to longitudinally track clonal repopulation in small animal models extensively used for stem cell and gene therapy research and with lentiviral vectors designed for clinical applications. Results of this study provide a framework for understanding the clonal behavior of human HPSC repopulating in a mouse environment, critical for translating results from humanized mice models to the human settings.

Modeling large fluctuations of thousands of clones during hematopoiesis: The role of stem cell self-renewal and bursty progenitor dynamics in rhesus macaque

In a recent clone-tracking experiment, millions of uniquely tagged hematopoietic stem cells (HSCs) and progenitor cells were autologously transplanted into rhesus macaques and peripheral blood containing thousands of tags were sampled and sequenced over 14 years to quantify the abundance of hundreds to thousands of tags or “clones.” Two major puzzles of the data have been observed: consistent differences and massive temporal fluctuations of clone populations. The large sample-to-sample variability can lead clones to occasionally go “extinct” but “resurrect” themselves in subsequent samples. Although heterogeneity in HSC differentiation rates, potentially due to tagging, and random sampling of the animals’ blood and cellular demographic stochasticity might be invoked to explain these features, we show that random sampling cannot explain the magnitude of the temporal fluctuations. Moreover, we show through simpler neutral mechanistic and statistical models of hematopoiesis of tagged cells that a broad distribution in clone sizes can arise from stochastic HSC self-renewal instead of tag-induced heterogeneity. The very large clone population fluctuations that often lead to extinctions and resurrections can be naturally explained by a generation-limited proliferation constraint on the progenitor cells. This constraint leads to bursty cell population dynamics underlying the large temporal fluctuations. We analyzed experimental clone abundance data using a new statistic that counts clonal disappearances and provided least-squares estimates of two key model parameters in our model, the total HSC differentiation rate and the maximum number of progenitor-cell divisions.

Hematopoiesis of virally tagged cells in rhesus macaques is analyzed in the context of a mechanistic and statistical model. We find that the clone size distribution and the temporal variability in the abundance of each clone (viral tag) in peripheral blood are consistent with (i) stochastic HSC self-renewal during bone marrow repair, (ii) clonal aging that restricts the number of generations of progenitor cells, and (iii) infrequent and small-size samples. By fitting data, we infer two key parameters that control the level of fluctuations of clone sizes in our model: the total HSC differentiation rate and the maximum proliferation capacity of progenitor cells. Our analysis provides insight into the mechanisms of hematopoiesis and a framework to guide future multiclone barcoding/lineage tracking measurements.

Purging Exhausted Virus-Specific CD8 T Cell Phenotypes by Somatic Cell Reprogramming

Cytotoxic T cells are critical in controlling virus infections. However, continuous antigen stimulation and negative regulatory factors cause CD8 T cells to enter a dysfunctional state (T cell exhaustion), resulting in viral persistence. We hypothesized that the exhausted T cell state could be molecularly rejuvenated using a somatic cell reprogramming technology, which is technically able to convert any types of cells to induced pluripotent stem cells (iPSCs), to regenerate functional T cells capable of purging chronic infection. We generated a new mouse line (B6/129^OKSM) in which every somatic cell contains four doxycycline-inducible reprogramming genes (Oct4, Klf4, Sox2, and c-Myc: OKSM), and infected them with lymphocytic choriomeningitis virus (LCMV) clone 13 to establish chronic infection. Exhausted LCMV-specific T cells isolated by flow sorting were successfully reprogrammed ex vivo into iPSCs in the presence of doxycycline. Upon injection into blastocysts and subsequent transfer into foster females, the reprogrammed cells differentiated into functional naive T cells that maintained their original antigen specificity. These results provide proof of concept that somatic cell reprogramming of exhausted T cells into iPSCs can erase imprints of their previous exhausted state and in turn regenerate functional virus-specific T cells.

Bidirectional Retroviral Integration Site PCR Methodology and Quantitative Data Analysis Workflow

Integration Site (IS) assays are a critical component of the study of retroviral integration sites and their biological significance. In recent retroviral gene therapy studies, IS assays, in combination with next-generation sequencing, have been used as a cell-tracking tool to characterize clonal stem cell populations sharing the same IS. For the accurate comparison of repopulating stem cell clones within and across different samples, the detection sensitivity, data reproducibility, and high-throughput capacity of the assay are among the most important assay qualities. This work provides a detailed protocol and data analysis workflow for bidirectional IS analysis. The bidirectional assay can simultaneously sequence both upstream and downstream vector-host junctions. Compared to conventional unidirectional IS sequencing approaches, the bidirectional approach significantly improves IS detection rates and the characterization of integration events at both ends of the target DNA. The data analysis pipeline described here accurately identifies and enumerates identical IS sequences through multiple steps of comparison that map IS sequences onto the reference genome and determine sequencing errors. Using an optimized assay procedure, we have recently published the detailed repopulation patterns of thousands of Hematopoietic Stem Cell (HSC) clones following transplant in rhesus macaques, demonstrating for the first time the precise time point of HSC repopulation and the functional heterogeneity of HSCs in the primate system. The following protocol describes the step-by-step experimental procedure and data analysis workflow that accurately identifies and quantifies identical IS sequences.

Clonal Analysis of Cells with Cellular Barcoding: When Numbers and Sizes Matter

Cellular barcoding is a recently rediscovered tool to trace the clonal output of individual cells with genetically distinct and heritable DNA sequences. Each year a few dozens of papers are published using the cellular barcoding technique. Those publications largely focus on mutually related issues, namely: counting cells capable of clonal proliferation and expansion, monitoring clonal dynamics in time, tracing the origin of differentiated cells, characterizing the differentiation potential of stem cells and similar topics. Apart from their biological content, claims and conclusions, these studies show remarkable diversity in technical aspects of the barcoding method and sometimes in major conclusions. Although a diversity of approaches is quite usual in data analysis, deviant handling of barcode data might directly affect experimental results and their biological interpretation. Here, we will describe typical challenges and caveats in cellular barcoding publications available so far.

Mechanisms of blood homeostasis: lineage tracking and a neutral model of cell populations in rhesus macaques

Background

How a potentially diverse population of hematopoietic stem cells (HSCs) differentiates and proliferates to supply more than 10¹¹ mature blood cells every day in humans remains a key biological question. We investigated this process by quantitatively analyzing the clonal structure of peripheral blood that is generated by a population of transplanted lentivirus-marked HSCs in myeloablated rhesus macaques. Each transplanted HSC generates a clonal lineage of cells in the peripheral blood that is then detected and quantified through deep sequencing of the viral vector integration sites (VIS) common within each lineage. This approach allowed us to observe, over a period of 4-12 years, hundreds of distinct clonal lineages.

Results

While the distinct clone sizes varied by three orders of magnitude, we found that collectively, they form a steady-state clone size-distribution with a distinctive shape. Steady-state solutions of our model show that the predicted clone size-distribution is sensitive to only two combinations of parameters. By fitting the measured clone size-distributions to our mechanistic model, we estimate both the effective HSC differentiation rate and the number of active HSCs.

Conclusions

Our concise mathematical model shows how slow HSC differentiation followed by fast progenitor growth can be responsible for the observed broad clone size-distribution. Although all cells are assumed to be statistically identical, analogous to a neutral theory for the different clone lineages, our mathematical approach captures the intrinsic variability in the times to HSC differentiation after transplantation.

Electronic supplementary material

The online version of this article (doi:10.1186/s12915-015-0191-8) contains supplementary material, which is available to authorized users.

Quantitative shearing linear amplification polymerase chain reaction: an improved method for quantifying lentiviral vector insertion sites in transplanted hematopoietic cell systems

In gene therapy trials targeting blood disorders, it is important to detect dominance of transduced hematopoietic stem cell (HSC) clones arising from vector insertion site (VIS) effects. Current methods for VIS analysis often do not have defined levels of quantitative accuracy and therefore can fail to detect early clonal dominance. We have developed a rapid and inexpensive method for measuring clone size based on random shearing of genomic DNA, minimal exponential PCR amplification, and shear site counts as a quantitative endpoint. This quantitative shearing linear amplification PCR (qsLAM PCR) assay utilizes an internal control sample containing 19 lentiviral insertion sites per cell that is mixed with polyclonal samples derived from transduced human CD34+ cells. Samples were analyzed from transplanted pigtail macaques and from a participant in our X-linked severe combined immunodeficiency (XSCID) lentiviral vector trial and yielded controlled and quantitative results in all cases. One case of early clonal dominance was detected in a monkey transplanted with limiting numbers of transduced HSCs, while the clinical samples from the XSCID trial participant showed highly diverse clonal representation. These studies demonstrate that qsLAM PCR is a facile and quantitative assay for measuring clonal repertoires in subjects enrolled in human gene therapy trials using lentiviral-transduced HSCs.

Dynamics of HSPC repopulation in nonhuman primates revealed by a decade-long clonal-tracking study

In mice, clonal tracking of hematopoietic stem cells (HSCs) has revealed variations in repopulation characteristics. However, it is unclear whether similar properties apply in primates. Here, we examined this issue through tracking of thousands of hematopoietic stem and progenitor cells (HSPCs) in rhesus macaques for up to 12 years. Approximately half of the clones analyzed contributed to long-term repopulation (over 3-10 years), arising in sequential groups and likely representing self-renewing HSCs. The remainder contributed primarily for the first year. The long-lived clones could be further subdivided into functional groups contributing primarily to myeloid, lymphoid, or both myeloid and lymphoid lineages. Over time, the 4%-10% of clones with robust dual lineage contribution predominated in repopulation. HSPCs expressing a CCR5 shRNA transgene behaved similarly to controls. Our study therefore documents HSPC behavior in a clinically relevant model over a long time frame and provides a substantial system-level data set that is a reference point for future work.

Lentiviral and targeted cellular barcoding reveals ongoing clonal dynamics of cell lines in vitro and in vivo

Background

Cell lines are often regarded as clonal, even though this simplifies what is known about mutagenesis, transformation and other processes that destabilize them over time. Monitoring these clonal dynamics is important for multiple areas of biomedical research, including stem cell and cancer biology. Tracking the contributions of individual cells to large populations, however, has been constrained by limitations in sensitivity and complexity.

Results

We utilize cellular barcoding methods to simultaneously track the clonal contributions of tens of thousands of cells. We demonstrate that even with optimal culturing conditions, common cell lines including HeLa, K562 and HEK-293 T exhibit ongoing clonal dynamics. Starting a population with a single clone diminishes but does not eradicate this phenomenon. Next, we compare lentiviral and zinc-finger nuclease barcode insertion approaches, finding that the zinc-finger nuclease protocol surprisingly results in reduced clonal diversity. We also document the expected reduction in clonal complexity when cells are challenged with genotoxic stress. Finally, we demonstrate that xenografts maintain clonal diversity to a greater extent than in vitro culturing of the human non-small-cell lung cancer cell line HCC827.

Conclusions

We demonstrate the feasibility of tracking and quantifying the clonal dynamics of entire cell populations within multiple cultured cell lines. Our results suggest that cell heterogeneity should be considered in the design and interpretation of in vitro culture experiments. Aside from clonal cell lines, we propose that cellular barcoding could prove valuable in modeling the clonal behavior of heterogeneous cell populations over time, including tumor populations treated with chemotherapeutic agents.

No impact of lentiviral transduction on hematopoietic stem/progenitor cell telomere length or gene expression in the rhesus macaque model

The occurrence of clonal perturbations and leukemia in patients transplanted with gamma-retroviral (RV) vector-transduced autologous hematopoietic stem and progenitor cells (HSPCs) has stimulated extensive investigation, demonstrating that proviral insertions may perturb adjacent proto-oncogene expression. Although enhancer-deleted lentiviruses are less likely to result in insertional oncogenesis, there is evidence that they may perturb transcript splicing, and one patient with a benign clonal expansion of lentivirally transduced HPSC has been reported. The rhesus macaque model provides an opportunity for informative long-term analysis to ask whether transduction impacts on long-term HSPC properties. We used two techniques to examine whether lentivirally transduced HSPCs from eight rhesus macaques transplanted 1-13.5 years previously are perturbed at a population level, comparing telomere length as a measure of replicative history and gene expression profile of vector positive versus vector negative cells. There were no differences in telomere lengths between sorted GFP+ and GFP- blood cells, suggesting that lentiviral (LV) transduction did not globally disrupt replicative patterns. Bone marrow GFP+ and GF- CD34+ cells showed no differences in gene expression using unsupervised and principal component analysis. These studies did not uncover any global long-term perturbation of proliferation, differentiation, or other important functional parameters of transduced HSPCs in the rhesus macaque model.

Evaluating a ligation-mediated PCR and pyrosequencing method for the detection of clonal contribution in polyclonal retrovirally transduced samples

Retroviral gene transfer has proven therapeutic potential in clinical gene therapy trials but may also cause abnormal cell growth via perturbation of gene expression in the locus surrounding the insertion site. By establishing clonal marks, retroviral insertions are also used to describe the regenerative potential of individual cells. Deep sequencing approaches have become the method of choice to study insertion profiles in preclinical models and clinical trials. We used a protocol combining ligation-mediated polymerase chain reaction (LM-PCR) and pyrosequencing for insertion profiling and quantification in cells of various tissues transduced with various retroviral vectors. The presented method allows simultaneous analysis of a multitude of DNA-barcoded samples per pyrosequencing run, thereby allowing cost-effective insertion screening in studies with multiple samples. In addition, we investigated whether the number of pyrosequencing reads can be used to quantify clonal abundance. By comparing pyrosequencing reads against site-specific quantitative PCR and by performing spike-in experiments, we show that considerable variation exists in the quantification of insertion sites even when present in the same clone. Our results suggest that the protocol used here and similar approaches might misinterpret abundance clones defined by insertion sites, unless careful calibration measures are taken. The crucial variables causing this variation need to be defined and methodological improvements are required to establish pyrosequencing reads as a quantification measure in polyclonal situations.

Tracking single hematopoietic stem cells in vivo using high-throughput sequencing in conjunction with viral genetic barcoding

Disentangling cellular heterogeneity is a challenge in many fields, particularly in the stem cell and cancer biology fields. Here, we demonstrate how to combine viral genetic barcoding with high-throughput sequencing to track single cells in a heterogeneous population. We use this technique to track the in vivo differentiation of unitary hematopoietic stem cells (HSCs). The results are consistent with single cell transplantation studies, but require two orders of magnitude fewer mice. In addition to its high throughput, the high sensitivity of the technique allows for a direct examination of the clonality of sparse cell populations such as HSCs. We show how these capabilities offer a clonal perspective of the HSC differentiation process. In particular, our data suggests that HSCs do not equally contribute to blood cells after irradiation-mediated transplantation, and that two distinct HSC differentiation patterns co-exist in the same recipient mouse post irradiation. This technique can be applied to any viral accessible cell type for both in vitro and in vivo processes.