The infusion of donor lymphocytes transduced with a retroviral vector expressing the HSV-TK suicide gene in patients undergoing hematopoietic stem cell transplantation for leukemia/lymphoma promotes immune reconstitution and prevents infections and graft-versus-host disease. Analysis of the clonal dynamics of genetically modified lymphocytes in vivo is of crucial importance to understand the potential genotoxic risk of this therapeutic approach. We used linear amplification-mediated PCR and pyrosequencing to build a genome-wide, high-definition map of retroviral integration sites in the genome of peripheral blood T cells from two different donors and used gene expression profiling and bioinformatics to associate integration clusters to transcriptional activity and to genetic and epigenetic features of the T cell genome. Comparison with matched random controls and with integrations obtained from CD34(+) hematopoietic stem/progenitor cells showed that integration clusters occur within chromatin regions bearing epigenetic marks associated with active promoters and regulatory elements in a cell-specific fashion. Analysis of integration sites in T cells obtained ex vivo two months after infusion showed no evidence of integration-related clonal expansion or dominance, but rather loss of cells harboring integration events interfering with RNA post-transcriptional processing. The study shows that high-definition maps of retroviral integration sites are a powerful tool to analyze the fate of genetically modified T cells in patients and the biological consequences of retroviral transduction.

Clinical trials have demonstrated the potential of ex vivo hematopoietic stem cell gene therapy to treat X-linked severe combined immunodeficiency (SCID-X1) using γ-retroviral vectors, leading to immune system functionality in the majority of treated patients without pretransplant conditioning. The success was tempered by insertional oncogenesis in a proportion of the patients. To reduce the genotoxicity risk, a self-inactivating (SIN) lentiviral vector (LV) with improved expression of a codon optimized human interleukin-2 receptor γ gene (IL2RG) cDNA (coγc), regulated by its 1.1 kb promoter region (γcPr), was compared in efficacy to the viral spleen focus forming virus (SF) and the cellular phosphoglycerate kinase (PGK) promoters. Pretransplant conditioning of Il2rg(-/-) mice resulted in long-term reconstitution of T and B lymphocytes, normalized natural antibody titers, humoral immune responses, ConA/IL-2 stimulated spleen cell proliferation, and polyclonal T-cell receptor gene rearrangements with a clear integration preference of the SF vector for proto-oncogenes, contrary to the PGK and γcPr vectors. We conclude that SIN lentiviral gene therapy using coγc driven by the γcPr or PGK promoter corrects the SCID phenotype, potentially with an improved safety profile, and that low-dose conditioning proved essential for immune competence, allowing for a reduced threshold of cell numbers required.

Retroviral vectors are widely used in gene therapy to introduce therapeutic genes into patients' cells, since, once delivered to the nucleus, the genes of interest are stably inserted (integrated) into the target cell genome. There is now compelling evidence that integration of retroviral vectors follows non-random patterns in mammalian genome, with a preference for active genes and regulatory regions. In particular, Moloney Leukemia Virus (MLV)-derived vectors show a tendency to integrate in the proximity of the transcription start site (TSS) of genes, occasionally resulting in the deregulation of gene expression and, where proto-oncogenes are targeted, in tumor initiation. This has drawn the attention of the scientific community to the molecular determinants of the retroviral integration process as well as to statistical methods to evaluate the genome-wide distribution of integration sites. In recent approaches, the observed distribution of MLV integration distances (IDs) from the TSS of the nearest gene is assumed to be non-random by empirical comparison with a random distribution generated by computational simulation procedures. To provide a statistical procedure to test the randomness of the retroviral insertion pattern, we propose a probability model (Beta distribution) based on IDs between two consecutive genes. We apply the procedure to a set of 595 unique MLV insertion sites retrieved from human hematopoietic stem/progenitor cells. The statistical goodness of fit test shows the suitability of this distribution to the observed data. Our statistical analysis confirms the preference of MLV-based vectors to integrate in promoter-proximal regions.