Search

Gene-corrected cells in Gene Therapy (GT) treated patients can be tracked in vivo by means of vector integration site (IS) analysis, since each engineered clone becomes univocally and stably marked by an individual IS. As the proper IS identification and quantification is crucial to accurately perform clonal tracking studies, we designed a customizable and tailored pipeline to analyze LAM-PCR amplicons sequenced by Illumina MiSeq/HiSeq technology. The sequencing data are initially processed through a series of quality filters and cleaned from vector and Linker Cassette (LC) sequences with customizable settings. Demultiplexing is then performed according to the recognition of specific barcodes combination used upon library preparation and the sequences are aligned to the reference genome. Importantly, the human genome assembly Hg19 is composed of 93 contigs, among which the mitochondrial genome, unlocalized and unplaced contigs and some alternative haplotypes of chr6. While previous approaches aligned IS sequences only to the standard 24 human chromosomes, using the whole assembled genome allowed improving alignment accuracy and concomitantly increased the amount of detectable ISs. To date, we have processed 28 independent human sample sets retrieving 260,994 ISs from 189,270,566 sequencing reads. Although, sequencing read counts at each IS have been widely used to estimate the relative IS abundance, this method carries inherent accuracy constraints due to the rounds of exponential amplification required by LAM-PCR that might generate unbalances on the original clonal representation. More recently, a method based on genomic sonication has been proposed exploiting shear site counts to tag the number of original fragments belonging to each IS before PCR amplification. However, the number of cells composing a given clone could far exceed the number of fragments of different lengths that can be generated upon fragmentation in proximity of that given IS. This would rapidly saturate the available diversity of shear sites and progressively generate more and more same-site shearing on independent genomes. In order to overcome the described biases and reliably quantify ISs, we designed and tested a new LC encoding random barcodes. The new LC is composed of a known sequence of 29nt used as binding site for the primers upon amplification steps, a 6nt-random barcode, a fixed-anchor sequence of 6nt, a second 6nt-random barcode and a final known sequence of 22nt containing sticky ends for the three main restriction enzymes in use (MluI, HpyCH4IV and AciI). This peculiar design allowed increasing the accuracy of clonal diversity estimation since the fixed-anchor sequence acts as a control for sequencing reliability in the barcode area. The theoretical number of different available barcodes per clone (412=16,777,216) far exceeds the requirements for not saturating the original diversity of the analyzed sample (on average composed by around 50.000 cells). We validated this novel approach by performing assays on serial dilutions of individual clones carrying known ISs. The precision rate obtained was averagely around 99.3%, while the worst error rate reaches at most the 1.86%, confirming the reliability of IS quantification. We successfully applied the barcoded-LC system to the analysis of clinical samples from a Wiskott Aldrich Syndrome GT patient, collecting to date 50,215 barcoded ISs from 94,052,785 sequencing reads.

Upon gene therapy (GT) for adenosine deaminase (ADA) deficient-SCID and Wiskott-Aldrich Syndrome (WAS), gene-corrected hematopoietic stem/progenitor cells (HSPC) generated a stable genetically engineered hematopoietic system where each vector-marked cell is univocally barcoded by a vector integration site (IS). To study the dynamics of human hematopoietic system, we collected by LAM-PCR+Illumina sequencing 28.539.414 sequence reads corresponding to 89.373 IS tagging clones belonging to 13 different cell types purified from the bone marrow and the peripheral blood of 4 WAS patients up to 48 months after GT. We firstly identified identical IS shared among CD34+ progenitors, and mature Myeloid/Lymphoid cells as markers of the real-time clonal output of vector-marked HSPC clones in vivo. We unraveled the nature and timing of short, intermediate and long term HSPC output showing that CD34+ clones active at 3-6 months after GT are not detectable at later follow up. Distinct waves of HSPC diversity were observed during the first 6-9 months after GT reaching a homeostatic equilibrium only by 12 months. By unsupervised clustering of IS similarities among lineages we unveiled diverse HSPC output towards lymphoid, myeloid and megakaryo-erythroid cells showing hints of a NK cells origin distinct from T and B cells. We exploited IS similarities to infere and test hematopoietic hierarchies by combining conditional probability distributions and static/dynamic graphical models of dependencies. We also estimated by mark-recapture approaches that just few thousands clones (1185-2884) are responsible for the long-term maintenance of the whole genetically engineered hematopoietic system. The ongoing analyses on IS collected from 7 distinct CD34+ subtypes isolated from GT patients will further increase HSPC tracking resolution. To evaluate the preservation of activity by transplanted HSPC we exploited data derived from the IS-based tracking of 4.845 clones in ADA-SCID patients performed for up to 6 years after GT. We showed that identical IS are consistently detected at multiple lineages level even several years after GT. Strikingly, by semi-quantitative PCRs on specific vector-genome junctions we tracked a fluctuating but consistent output of marked HSPC over a period of 5 years without the manifestation of clonal quiescence phases. Since gamma-retroviral vector used in this GT trial actively transduce only replicating cells, this provide the first evidence that HSPC, awaken from dormancy in vitro, can still retain in vivo long-term activity. Overall our work constitute the first molecular tracking of hematopoietic system in humans. The information gathered are crucial for the design of therapeutic approaches for a broad spectrum of hematological diseases and tumors.