Andrea Calabria

Hematopoietic Stem Cell (HSC) gene therapy (GT) may provide lifelong reconstitution of the hematopoietic system with gene-corrected cells¹. However, the effects of underlying genetic diseases, replication stress, and aging on hematopoietic reconstitution and lineage specification remain unclear. In this study, we analyzed hematopoietic reconstitution in 53 patients treated with lentiviral-HSC-GT for diverse conditions such as metachromatic leukodystrophy^2,3 (MLD), Wiskott-Aldrich syndrome^4,5 (WAS), and β-thalassemia⁶ (β-Thal) over a follow-up period of up to 8 years, using vector integration sites as markers of clonal identity. We found that long-term hematopoietic reconstitution was supported by 770 to 35,000 active HSCs. While 50% of transplanted clones demonstrated multilineage potential across all conditions, the remaining clones exhibited a disease-specific preferential lineage output and long-term commitment: myeloid for MLD, lymphoid for WAS, and erythroid for β-Thal, particularly in adult patients. Our results indicate that HSC clonogenic activity, lineage output, long-term lineage commitment, and rates of somatic mutations are influenced by the underlying disease, patient age at the time of therapy, the extent of genetic defect correction, and the hematopoietic stress imposed by the inherited disease. This suggests that HSCs adapt to the pathological condition during hematopoietic reconstitution.

In physiological conditions, few circulating hematopoietic stem/progenitor cells (cHSPCs) are present in the peripheral blood, but their contribution to human hematopoiesis remain unsolved. By integrating advanced immunophenotyping, single-cell transcriptional and functional profiling, and integration site (IS) clonal tracking, we unveiled the biological properties and the transcriptional features of human cHSPC subpopulations in relationship to their bone marrow (BM) counterpart. We found that cHSPCs reduced in cell count over aging and are enriched for primitive, lymphoid, and erythroid subpopulations, showing preactivated transcriptional and functional state. Moreover, cHSPCs have low expression of multiple BM-retention molecules but maintain their homing potential after xenotransplantation. By generating a comprehensive human organ-resident HSPC data set based on single-cell RNA sequencing data, we detected organ-specific seeding properties of the distinct trafficking HSPC subpopulations. Notably, circulating multi-lymphoid progenitors are primed for seeding the thymus and actively contribute to T-cell production. Human clonal tracking data from patients receiving gene therapy (GT) also showed that cHSPCs connect distant BM niches and participate in steady-state hematopoietic production, with primitive cHSPCs having the highest recirculation capability to travel in and out of the BM. Finally, in case of hematopoietic impairment, cHSPCs composition reflects the BM-HSPC content and might represent a biomarker of the BM state for clinical and research purposes. Overall, our comprehensive work unveiled fundamental insights into the in vivo dynamics of human HSPC trafficking and its role in sustaining hematopoietic homeostasis. GT patients’ clinical trials were registered at ClinicalTrials.gov (NCT01515462 and NCT03837483) and EudraCT (2009-017346-32 and 2018-003842-18).

Hematopoietic stem cell gene therapy (GT) using a γ-retroviral vector (γ-RV) is an effective treatment for Severe Combined Immunodeficiency due to Adenosine Deaminase deficiency. Here, we describe a case of GT-related T-cell acute lymphoblastic leukemia (T-ALL) that developed 4.7 years after treatment. The patient underwent chemotherapy and haploidentical transplantation and is currently in remission. Blast cells contain a single vector insertion activating the LIM-only protein 2 (LMO2) proto-oncogene, confirmed by physical interaction, and low Adenosine Deaminase (ADA) activity resulting from methylation of viral promoter. The insertion is detected years before T-ALL in multiple lineages, suggesting that further hits occurred in a thymic progenitor. Blast cells contain known and novel somatic mutations as well as germline mutations which may have contributed to transformation. Before T-ALL onset, the insertion profile is similar to those of other ADA-deficient patients. The limited incidence of vector-related adverse events in ADA-deficiency compared to other γ-RV GT trials could be explained by differences in transgenes, background disease and patient’s specific factors.

Mobilized peripheral blood is increasingly used instead of bone marrow as a source of autologous hematopoietic stem/progenitor cells for ex vivo gene therapy. Here, we present an unplanned exploratory analysis evaluating the hematopoietic reconstitution kinetics, engraftment and clonality in 13 pediatric Wiskott-Aldrich syndrome patients treated with autologous lentiviral-vector transduced hematopoietic stem/progenitor cells derived from mobilized peripheral blood (n = 7), bone marrow (n = 5) or the combination of the two sources (n = 1). 8 out of 13 gene therapy patients were enrolled in an open-label, non-randomized, phase 1/2 clinical study (NCT01515462) and the remaining 5 patients were treated under expanded access programs. Although mobilized peripheral blood- and bone marrow- hematopoietic stem/progenitor cells display similar capability of being gene-corrected, maintaining the engineered grafts up to 3 years after gene therapy, mobilized peripheral blood-gene therapy group shows faster neutrophil and platelet recovery, higher number of engrafted clones and increased gene correction in the myeloid lineage which correlate with higher amount of primitive and myeloid progenitors contained in hematopoietic stem/progenitor cells derived from mobilized peripheral blood. In vitro differentiation and transplantation studies in mice confirm that primitive hematopoietic stem/progenitor cells from both sources have comparable engraftment and multilineage differentiation potential. Altogether, our analyses reveal that the differential behavior after gene therapy of hematopoietic stem/progenitor cells derived from either bone marrow or mobilized peripheral blood is mainly due to the distinct cell composition rather than functional differences of the infused cell products, providing new frames of references for clinical interpretation of hematopoietic stem/progenitor cell transplantation outcome.

Adeno-associated Viral vectors (AAVs) have been successful exploited in gene therapy applications for the treatment of several genetic disorders. AAV is considered an episomal vector but it has been shown to integrate within the host-cell genome following generation of double-strand DNA breaks or nicks. While AAV integration raises some safety concerns, it can also provide therapeutic benefit; the direct intrathymic injection of an AAV harboring a therapeutic transgene results in integration in T cell progenitors and long-term T cell immunity. To assess the mechanisms of AAV integration, we retrieved and analyzed hundreds of AAV integration sites from lymph node-derived mature T-cells as compared to liver and brain tissue from treated mice. Notably, we found that while AAV integrations in the liver and brain were distributed across the entire mouse genome, >90% of the integrations in T-cells were clustered within the T cell receptor α, β and γ genes. More precisely, the insertion mapped to DNA breaks created by the enzymatic activity of recombination activating genes (RAG) during V(D)J recombination. Our data indicate that RAG activity during T-cell receptor maturation induces a site-specific integration of AAV genomes and opens new therapeutic avenues for achieving long-term AAV-mediated gene transfer in dividing cells.

High-throughput clonal tracking in patients under hematopoietic stem cell gene therapy with integrating vector is instrumental in assessing bio-safety and efficacy. Monitoring the fate of millions of transplanted clones and their progeny across differentiation and proliferation over time leverages the identification of the vector integration sites, used as surrogates of clonal identity. Although γ-tracking retroviral insertion sites (γ-TRIS) is the state-of-the-art algorithm for clonal identification, the computational drawbacks in the tracking algorithm, based on a combinatorial all-versus-all strategy, limit its use in clinical studies with several thousands of samples per patient. We developed the first clonal tracking graph database, InCliniGene (https://github.com/calabrialab/InCliniGene), that imports the output files of γ-TRIS and generates the graph of clones (nodes) connected by arches if two nodes share common genomic features as defined by the γ-TRIS rules. Embedding both clonal data and their connections in the graph, InCliniGene can track all clones longitudinally over samples through data queries that fully explore the graph. This approach resulted in being highly accurate and scalable. We validated InCliniGene using an in vitro dataset, specifically designed to mimic clinical cases, and tested the accuracy and precision. InCliniGene allows extensive use of γ-TRIS in large gene therapy clinical applications and naturally realizes the full data integration of molecular and genomics data, clinical and treatment measurements and genomic annotations. Further extensions of InCliniGene with data federation and with application programming interface will support data mining toward precision, personalized and predictive medicine in gene therapy. Database URL: https://github.com/calabrialab/InCliniGene

Longitudinal clonal tracking studies based on high-throughput sequencing technologies supported safety and long-term efficacy and unraveled hematopoietic reconstitution in many gene therapy applications with unprecedented resolution. However, monitoring patients over a decade-long follow-up entails a constant increase of large data volume with the emergence of critical computational challenges, unfortunately not addressed by currently available tools. Here we present ISAnalytics, a new R package for comprehensive and high-throughput clonal tracking studies using vector integration sites as markers of cellular identity. Once identified the clones externally from ISAnalytics and imported in the package, a wide range of implemented functionalities are available to users for assessing the safety and long-term efficacy of the treatment, here described in a clinical trial use case for Hurler disease, and for supporting hematopoietic stem cell biology in vivo with longitudinal analysis of clones over time, proliferation and differentiation. ISAnalytics is conceived to be metadata-driven, enabling users to focus on biological questions and hypotheses rather than on computational aspects. ISAnalytics can be fully integrated within laboratory workflows and standard procedures. Moreover, ISAnalytics is designed with efficient and scalable data structures, benchmarked with previous methods, and grants reproducibility and full analytical control through interactive web-reports and a module with Shiny interface. The implemented functionalities are flexible for all viral vector-based clonal tracking applications as well as genetic barcoding or cancer immunotherapies.

While the role of common genetic variants in multiple sclerosis (MS) has been elucidated in large genome-wide association studies, the contribution of rare variants to the disease remains unclear. Herein, a whole-genome sequencing study in four affected and four healthy relatives of a consanguineous Italian family identified a novel missense c.1801T > C (p.S601P) variant in the GRAMD1B gene that is shared within MS cases and resides under a linkage peak (LOD: 2.194). Sequencing GRAMD1B in 91 familial MS cases revealed two additional rare missense and two splice-site variants, two of which (rs755488531 and rs769527838) were not found in 1000 Italian healthy controls. Functional studies demonstrated that GRAMD1B, a gene with unknown function in the central nervous system (CNS), is expressed by several cell types, including astrocytes, microglia and neurons as well as by peripheral monocytes and macrophages. Notably, GRAMD1B was downregulated in vessel-associated astrocytes of active MS lesions in autopsied brains and by inflammatory stimuli in peripheral monocytes, suggesting a possible role in the modulation of inflammatory response and disease pathophysiology.

High transduction rates of viral vectors in gene therapies (GT) and experimental hematopoiesis ensure a high frequency of gene delivery, although multiple integration events can occur in the same cell. Therefore, tracing of integration sites (IS) leads to mis-quantification of the true clonal spectrum and limits safety considerations in GT. Hence, we use correlations between repeated measurements of IS abundances to estimate their mutual similarity and identify clusters of co-occurring IS, for which we assume a clonal origin. We evaluate the performance, robustness and specificity of our methodology using clonal simulations. The reconstruction methods, implemented and provided as an R-package, are further applied to experimental clonal mixes and preclinical models of hematopoietic GT. Our results demonstrate that clonal reconstruction from IS data allows to overcome systematic biases in the clonal quantification as an essential prerequisite for the assessment of safety and long-term efficacy of GT involving integrative vectors.

Viral vectors are used to insert genetic material into semirandom genomic positions of hematopoietic stem cells which, after reinfusion into patients, regenerate the entire hematopoietic system. Hematopoietic cells originating from genetically modified stem cells will harbor insertions in specific genomic positions called integration sites, which represent unique genetic marks of clonal identity. Therefore, the analysis of vector integration sites present in the genomic DNA of circulating cells allows to determine the number of clones in the blood ecosystem. Shannon diversity index is adopted to evaluate the heterogeneity of the transduced population of gene corrected cells. However, this measure can be affected by several technical variables such as the DNA amount used and the sequencing depth of the library analyzed and therefore the comparison across samples may be affected by these confounding factors. We developed an advanced spline-regression approach that leverages on confounding effects to provide a normalized entropy index. Our proposed method was first validated and compared with two state of the art approaches in a specifically designed in vitro assay. Subsequently our approach allowed to observe the expected impact of vector genotoxicity on entropy level decay in an in vivo model of hematopoietic stem cell gene therapy based on tumor prone mice.