Reduced adhesion of aged intestinal stem cells contributes to an accelerated clonal drift.

Upon aging, the function of the intestinal epithelium declines with a concomitant increase in aging-related diseases. ISCs play an important role in this process. It is known that ISC clonal dynamics follow a neutral drift model. However, it is not clear whether the drift model is still valid in aged ISCs. Tracking of clonal dynamics by clonal tracing revealed that aged crypts drift into monoclonality substantially faster than young ones. However, ISC tracing experiments, in vivo and ex vivo, implied a similar clonal expansion ability of both young and aged ISCs. Single-cell RNA sequencing for 1,920 high Lgr5 ISCs from young and aged mice revealed increased heterogeneity among subgroups of aged ISCs. Genes associated with cell adhesion were down-regulated in aged ISCs. ISCs of aged mice indeed show weaker adhesion to the matrix. Simulations applying a single cell-based model of the small intestinal crypt demonstrated an accelerated clonal drift at reduced adhesion strength, implying a central role for reduced adhesion for affecting clonal dynamics upon aging.

The Reconstitution Dynamics of Cultivated Hematopoietic Stem Cells and Progenitors Is Independent of Age.

Hematopoietic stem cell transplantation (HSCT) represents the only curative treatment option for numerous hematologic malignancies. While the influence of donor age and the composition of the graft have already been examined in clinical and preclinical studies, little information is available on the extent to which different hematological subpopulations contribute to the dynamics of the reconstitution process and on whether and how these contributions are altered with age. In a murine model of HSCT, we therefore simultaneously tracked different cultivated and transduced hematopoietic stem and progenitor cell (HSPC) populations using a multicolor-coded barcode system (BC32). We studied a series of age-matched and age-mismatched transplantations and compared the influence of age on the reconstitution dynamics. We show that reconstitution from these cultured and assembled grafts was substantially driven by hematopoietic stem cells (HSCs) and multipotent progenitors (MPPs) independent of age. The reconstitution patterns were polyclonal and stable in all age groups independently of the variability between individual animals, with higher output rates from MPPs than from HSCs. Our experiments suggest that the dynamics of reconstitution and the contribution of cultured and individually transduced HSPC subpopulations are largely independent of age. Our findings support ongoing efforts to expand the application of HSCT in older individuals as a promising strategy to combat hematological diseases, including gene therapy applications.

genBaRcode: a comprehensive R-package for genetic barcode analysis.

MOTIVATION: Genetic barcodes have been established as an efficient method to trace clonal progeny of uniquely labeled cells by introducing artificial genetic sequences into the corresponding genomes. The assessment of those sequences relies on next generation sequencing and the subsequent analysis aiming to identify sequences of interest and correctly quantifying their abundance. RESULTS: We developed the genBaRcode package as a toolbox combining the flexibility of digesting next generation sequencing reads with or without a sophisticated barcode structure, with a variety of error-correction approaches and the availability of several types of visualization routines. Furthermore, a graphical user interface was incorporated to allow also less experienced R users package-based analyses. Finally, the provided tool is intended to bridge the gap between generating and analyzing barcode data and thereby supporting the establishment of standardized and reproducible analysis strategies. AVAILABILITY AND IMPLEMENTATION: The genBaRcode package is available at CRAN (https://cran.r-project.org/package=genBaRcode).

Engrafted parenchymal brain macrophages differ from microglia in transcriptome, chromatin landscape and response to challenge.

Microglia are yolk sac-derived macrophages residing in the parenchyma of brain and spinal cord, where they interact with neurons and other glial. After different conditioning paradigms and bone marrow (BM) or hematopoietic stem cell (HSC) transplantation, graft-derived cells seed the brain and persistently contribute to the parenchymal brain macrophage compartment. Here we establish that graft-derived macrophages acquire, over time, microglia characteristics, including ramified morphology, longevity, radio-resistance and clonal expansion. However, even after prolonged CNS residence, transcriptomes and chromatin accessibility landscapes of engrafted, BM-derived macrophages remain distinct from yolk sac-derived host microglia. Furthermore, engrafted BM-derived cells display discrete responses to peripheral endotoxin challenge, as compared to host microglia. In human HSC transplant recipients, engrafted cells also remain distinct from host microglia, extending our finding to clinical settings. Collectively, our data emphasize the molecular and functional heterogeneity of parenchymal brain macrophages and highlight potential clinical implications for HSC gene therapies aimed to ameliorate lysosomal storage disorders, microgliopathies or general monogenic immuno-deficiencies.

Cellular Barcoding Identifies Clonal Substitution as a Hallmark of Local Recurrence in a Surgical Model of Head and Neck Squamous Cell Carcinoma.

Local recurrence after surgery for head and neck squamous cell carcinoma (HNSCC) remains a common event associated with a dismal prognosis. Improving this outcome requires a better understanding of cancer cell populations that expand from postsurgical minimal residual disease (MRD). Therefore, we assessed clonal dynamics in a surgical model of barcoded HNSCC growing in the submental region of immunodeficient mice. Clonal substitution and massive reduction of clonal heterogeneity emerged as hallmarks of local recurrence, as the clones dominating in less heterogeneous recurrences were scarce in their matched primary tumors. These lineages were selected by their ability to persist after surgery and competitively expand from MRD. Clones enriched in recurrences exhibited both private and shared genetic features and likely originated from ancestors shared with clones dominating in primary tumors. They demonstrated high invasiveness and epithelial-to-mesenchymal transition, eventually providing an attractive target for obtaining better local control for these tumors.

Hematopoietic stem cells can differentiate into restricted myeloid progenitors before cell division in mice.

Hematopoietic stem cells (HSCs) continuously replenish all blood cell types through a series of differentiation steps and repeated cell divisions that involve the generation of lineage-committed progenitors. However, whether cell division in HSCs precedes differentiation is unclear. To this end, we used an HSC cell-tracing approach and Ki67RFP knock-in mice, in a non-conditioned transplantation model, to assess divisional history, cell cycle progression, and differentiation of adult HSCs. Our results reveal that HSCs are able to differentiate into restricted progenitors, especially common myeloid, megakaryocyte-erythroid and pre-megakaryocyte progenitors, without undergoing cell division and even before entering the S phase of the cell cycle. Additionally, the phenotype of the undivided but differentiated progenitors correlated with the expression of lineage-specific genes and loss of multipotency. Thus HSC fate decisions can be uncoupled from physical cell division. These results facilitate a better understanding of the mechanisms that control fate decisions in hematopoietic cells.

Genetic Barcodes Facilitate Competitive Clonal Analyses In Vivo.

Monitoring the fate of individual cell clones is an important task to better understand normal tissue regeneration, for example after hematopoietic stem cell (HSC) transplantation, but also cancerogenesis. Based on their integration into the host cell's genome, retroviral vectors are commonly used to stably mark target cells and their progeny. The development of genetic barcoding techniques has opened new possibilities to determine clonal composition and dynamics in great detail. A modular genetic barcode was recently introduced consisting of 32 variable positions (BC32) with a customized backbone, and its advantages were demonstrated with regard to barcode calling and quantification. The study presented applied the BC32 system in a complex in vivo situation, namely to analyze clonal reconstitution dynamics for HSC grafts consisting of up to three cell populations with distinguishable barcodes using different alpha- and lentiviral vectors. In a competitive transplantation setup, it was possible to follow the differently marked cell populations within individual animals. This enabled the clonal contribution of the different BC32 constructs during reconstitution and long-term hematopoiesis in the peripheral blood and the spatial distribution in bone marrow and spleen to be identified. Thus, it was demonstrated that the system allows the output of individually marked cells to be tracked in vivo and their influence on clonal dynamics to be analyzed. Successful application of the BC32 system in a complex, competitive in vivo situation provided proof-of-principle that its high complexity and the large Hamming distance between individual barcodes, combined with the easy customization, facilitate efficient and precise quantification, even without prior knowledge of individual barcode sequences. Importantly, simultaneous high-sensitivity analyses of different cell populations in single animals may significantly reduce numbers of animals required to investigate specific scientific questions in accordance with RRR principles. It is concluded that this BC32 system will be excellently suited for various research applications in regenerative medicine and cancer biology.

Clonal competition in BcrAbl-driven leukemia: how transplantations can accelerate clonal conversion.

BACKGROUND: Clonal competition in cancer describes the process in which the progeny of a cell clone supersedes or succumbs to other competing clones due to differences in their functional characteristics, mostly based on subsequently acquired mutations. Even though the patterns of those mutations are well explored in many tumors, the dynamical process of clonal selection is underexposed. METHODS: We studied the dynamics of clonal competition in a BcrAbl-induced leukemia using a γ-retroviral vector library encoding the oncogene in conjunction with genetic barcodes. To this end, we studied the growth dynamics of transduced cells on the clonal level both in vitro and in vivo in transplanted mice. RESULTS: While we detected moderate changes in clonal abundancies in vitro, we observed monoclonal leukemias in 6/30 mice after transplantation, which intriguingly were caused by only two different BcrAbl clones. To analyze the success of these clones, we applied a mathematical model of hematopoietic tissue maintenance, which indicated that a differential engraftment capacity of these two dominant clones provides a possible explanation of our observations. These findings were further supported by additional transplantation experiments and increased BcrAbl transcript levels in both clones. CONCLUSION: Our findings show that clonal competition is not an absolute process based on mutations, but highly dependent on selection mechanisms in a given environmental context.

Limitations and challenges of genetic barcode quantification.

Genetic barcodes are increasingly used to track individual cells and to quantitatively assess their clonal contributions over time. Although barcode quantification relies entirely on counting sequencing reads, detailed studies about the method's accuracy are still limited. We report on a systematic investigation of the relation between barcode abundance and resulting read counts after amplification and sequencing using cell-mixtures that contain barcodes with known frequencies ("miniBulks"). We evaluated the influence of protocol modifications to identify potential sources of error and elucidate possible limitations of the quantification approach. Based on these findings we designed an advanced barcode construct (BC32) to improved barcode calling and quantification, and to ensure a sensitive detection of even highly diluted barcodes. Our results emphasize the importance of using curated barcode libraries to obtain interpretable quantitative data and underline the need for rigorous analyses of any utilized barcode library in terms of reliability and reproducibility.

Measles virus envelope pseudotyped lentiviral vectors transduce quiescent human HSCs at an efficiency without precedent.

Hematopoietic stem cell (HSC)-based gene therapy trials are now moving toward the use of lentiviral vectors (LVs) with success. However, one challenge in the field remains: efficient transduction of HSCs without compromising their stem cell potential. Here we showed that measles virus glycoprotein-displaying LVs (hemagglutinin and fusion protein LVs [H/F-LVs]) were capable of transducing 100% of early-acting cytokine-stimulated human CD34+ (hCD34+) progenitor cells upon a single application. Strikingly, these H/F-LVs also allowed transduction of up to 70% of nonstimulated quiescent hCD34+ cells, whereas conventional vesicular stomatitis virus G (VSV-G)-LVs reached 5% at the most with H/F-LV entry occurring exclusively through the CD46 complement receptor. Importantly, reconstitution of NOD/SCIDγc-/- (NSG) mice with H/F-LV transduced prestimulated or resting hCD34+ cells confirmed these high transduction levels in all myeloid and lymphoid lineages. Remarkably, for resting CD34+ cells, secondary recipients exhibited increasing transduction levels of up to 100%, emphasizing that H/F-LVs efficiently gene-marked HSCs in the resting state. Because H/F-LVs promoted ex vivo gene modification of minimally manipulated CD34+ progenitors that maintained stemness, we assessed their applicability in Fanconi anemia, a bone marrow (BM) failure with chromosomal fragility. Notably, only H/F-LVs efficiently gene-corrected minimally stimulated hCD34+ cells in unfractionated BM from these patients. These H/F-LVs improved HSC gene delivery in the absence of cytokine stimulation while maintaining their stem cell potential. Thus, H/F-LVs will facilitate future clinical applications requiring HSC gene modification, including BM failure syndromes, for which treatment has been very challenging up to now.

Model Based Analysis of Clonal Developments Allows for Early Detection of Monoclonal Conversion and Leukemia.

The availability of several methods to unambiguously mark individual cells has strongly fostered the understanding of clonal developments in hematopoiesis and other stem cell driven regenerative tissues. While cellular barcoding is the method of choice for experimental studies, patients that underwent gene therapy carry a unique insertional mark within the transplanted cells originating from the integration of the retroviral vector. Close monitoring of such patients allows accessing their clonal dynamics, however, the early detection of events that predict monoclonal conversion and potentially the onset of leukemia are beneficial for treatment. We developed a simple mathematical model of a self-stabilizing hematopoietic stem cell population to generate a wide range of possible clonal developments, reproducing typical, experimentally and clinically observed scenarios. We use the resulting model scenarios to suggest and test a set of statistical measures that should allow for an interpretation and classification of relevant clonal dynamics. Apart from the assessment of several established diversity indices we suggest a measure that quantifies the extension to which the increase in the size of one clone is attributed to the total loss in the size of all other clones. By evaluating the change in relative clone sizes between consecutive measurements, the suggested measure, referred to as maximum relative clonal expansion (mRCE), proves to be highly sensitive in the detection of rapidly expanding cell clones prior to their dominant manifestation. This predictive potential places the mRCE as a suitable means for the early recognition of leukemogenesis especially in gene therapy patients that are closely monitored. Our model based approach illustrates how simulation studies can actively support the design and evaluation of preclinical strategies for the analysis and risk evaluation of clonal developments.

Multiplexing clonality: combining RGB marking and genetic barcoding.

RGB marking and DNA barcoding are two cutting-edge technologies in the field of clonal cell marking. To combine the virtues of both approaches, we equipped LeGO vectors encoding red, green or blue fluorescent proteins with complex DNA barcodes carrying color-specific signatures. For these vectors, we generated highly complex plasmid libraries that were used for the production of barcoded lentiviral vector particles. In proof-of-principle experiments, we used barcoded vectors for RGB marking of cell lines and primary murine hepatocytes. We applied single-cell polymerase chain reaction to decipher barcode signatures of individual RGB-marked cells expressing defined color hues. This enabled us to prove clonal identity of cells with one and the same RGB color. Also, we made use of barcoded vectors to investigate clonal development of leukemia induced by ectopic oncogene expression in murine hematopoietic cells. In conclusion, by combining RGB marking and DNA barcoding, we have established a novel technique for the unambiguous genetic marking of individual cells in the context of normal regeneration as well as malignant outgrowth. Moreover, the introduction of color-specific signatures in barcodes will facilitate studies on the impact of different variables (e.g. vector type, transgenes, culture conditions) in the context of competitive repopulation studies.

Comparative clonal analysis of reconstitution kinetics after transplantation of hematopoietic stem cells gene marked with a lentiviral SIN or a γ-retroviral LTR vector.

Retroviral gene marking has been used successfully in preclinical and clinical transplantation settings. Highly sensitive techniques for vector insertion-site determination, such as linear amplification-mediated polymerase chain reaction (LAM-PCR) in conjunction with next-generation sequencing, have been introduced to assess the composition of gene-marked hematopoiesis at a single-cell level. Here we used these novel techniques for directly comparing clonal reconstitution kinetics in mice transplanted with bone-marrow-derived stem cells genetically marked with either a standard, spleen focus-forming virus long terminal repeat (LTR)-driven γ-retroviral, or a lentiviral self-inactivating vector containing an identical but internal spleen focus-forming virus-derived enhancer/promoter. We observed that the use of the lentiviral self-inactivating vector for gene marking was associated with a broader repertoire of differently marked hematopoietic clones. More importantly, we found a significantly higher probability of insertions in growth-promoting, clonal-dominance-associated genes in the spleen focus-forming virus LTR-driven γ-retroviral vector at later time points of analysis. Based on our data, we suggest that the combined use of LAM-PCR and next-generation sequencing represents a potent tool for the analysis of clonal reconstitution kinetics in the context of gene marking with integrated vectors. At the same time, our findings prove that the use of multiple restriction enzymes for LAM-PCR is indispensable to detect most or ideally all individual stem cell clones contributing to hematopoiesis. We have also found that techniques such as quantitative PCR can be helpful to retrospectively analyze reconstitution kinetics for individual hematopoietic stem cell clones. Finally, our results confirm the notion that marking with lentiviral self-inactivating vectors is associated with a lower risk of genotoxicity as compared with γ-retroviral LTR vectors.

Cellular aging leads to functional heterogeneity of hematopoietic stem cells: a modeling perspective.

Hematopoietic stem cells (HSCs) are the source for the life-long supply of functional cells in peripheral blood while they simultaneously maintain their own reserve pool. However, there is accumulating evidence that HSCs are themselves subject to quantitative and qualitative exhaustion. Although several processes linked to mitotic activity can potentially account for the observed aging phenomena (e.g., DNA damage, telomere shortening, epigenetic modification), a precise understanding of HSC exhaustion is still missing. It is particularly unclear how individual aging processes on the single-cell level translate on the phenotypic level of the overall tissue and whether there is a functional implication of an age-structured HSC population. We address these issues by applying a novel mathematical model of HSC organization in which division-specific, cumulative alterations of stem cell quality determine the phenotypic and functional appearance of the overall cell population. Adapting the model to a number of basic experimental findings, we quantify the level of additional heterogeneity that is introduced by a population of individually aging cells. Based on this model, we are able to conclude that division-dependent processes of cellular aging explain a wide range of phenomena on HSC exhaustion and that HSC aging needs to be considered as a highly heterogeneous process. We furthermore report that functional heterogeneity between young and old HSCs appears closely similar to the phenomena described for long- and short-term repopulating cells. We speculate whether differential, division-coupled stem cell aging introduces an intra-animal variability that also accounts for heterogeneity with respect to the repopulation ability of HSCs.

Stem cell proliferation and quiescence--two sides of the same coin.

The kinetics of label uptake and dilution in dividing stem cells, e.g., using Bromodeoxyuridine (BrdU) as a labeling substance, are a common way to assess the cellular turnover of all hematopoietic stem cells (HSCs) in vivo. The assumption that HSCs form a homogeneous population of cells which regularly undergo cell division has recently been challenged by new experimental results. For a consistent functional explanation of heterogeneity among HSCs, we propose a concept in which stem cells flexibly and reversibly adapt their cycling state according to systemic needs. Applying a mathematical model analysis, we demonstrate that different experimentally observed label dilution kinetics are consistently explained by the proposed model. The dynamically stabilized equilibrium between quiescent and activated cells leads to a biphasic label dilution kinetic in which an initial and pronounced decline of label retaining cells is attributed to faster turnover of activated cells, whereas a secondary, decelerated decline results from the slow turnover of quiescent cells. These results, which support our previous model prediction of a reversible activation/deactivation of HSCs, are also consistent with recent findings that use GFP-conjugated histones as a label instead of BrdU. Based on our findings we interpret HSC organization as an adaptive and regulated process in which the slow activation of quiescent cells and their possible return into quiescence after division are sufficient to explain the simultaneous occurrence of self-renewal and differentiation. Furthermore, we suggest an experimental strategy which is suited to demonstrate that the repopulation ability among the population of label retaining cells changes during the course of dilution.