GeIST: a pipeline for mapping integrated DNA elements

IS-Seq: a bioinformatics pipeline for integration sites analysis with comprehensive abundance quantification methods

BACKGROUND

Integration site (IS) analysis is a fundamental analytical platform for evaluating the safety and efficacy of viral vector based preclinical and clinical Gene Therapy (GT). A handful of groups have developed standardized bioinformatics pipelines to process IS sequencing data, to generate reports, and/or to perform comparative studies across different GT trials. Keeping up with the technological advances in the field of IS analysis, different computational pipelines have been published over the past decade. These pipelines focus on identifying IS from single-read sequencing or paired-end sequencing data either using read-based or using sonication fragment-based methods, but there is a lack of a bioinformatics tool that automatically includes unique molecular identifiers (UMI) for IS abundance estimations and allows comparing multiple quantification methods in one integrated pipeline.

RESULTS

Here we present IS-Seq a bioinformatics pipeline that can process data from paired-end sequencing of both old restriction sites-based IS collection methods and new sonication-based IS retrieval systems while allowing the selection of different abundance estimation methods, including read-based, Fragment-based and UMI-based systems.

CONCLUSIONS

We validated the performance of IS-Seq by testing it against the most popular  analytical workflow available in the literature (INSPIIRED) and using different scenarios. Lastly, by performing extensive simulation studies and a comprehensive wet-lab assessment of our IS-Seq pipeline we could show that in clinically relevant scenarios, UMI quantification provides better accuracy than the currently most widely used sonication fragment counts as a method for IS abundance estimation.

Low incidence of hepatocellular carcinoma in mice and cats treated with systemic adeno-associated viral vectors

Adeno-associated viral (AAV) vectors have emerged as the preferred platform for in vivo gene transfer because of their combined efficacy and safety. However, insertional mutagenesis with the subsequent development of hepatocellular carcinomas (HCCs) has been recurrently noted in newborn mice treated with high doses of AAV, and more recently, the association of wild-type AAV integrations in a subset of human HCCs has been documented. Here, we address, in a comprehensive, prospective study, the long-term risk of tumorigenicity in young adult mice following delivery of single-stranded AAVs targeting liver. HCC incidence in mice treated with therapeutic and reporter AAVs was low, in contrast to what has been previously documented in mice treated as newborns with higher doses of AAV. Specifically, HCCs developed in 6 out 76 of AAV-treated mice, and a pathogenic integration of AAV was found in only one tumor. Also, no evidence of liver tumorigenesis was found in juvenile AAV-treated mucopolysaccharidosis type VI (MPS VI) cats followed as long as 8 years after vector administration. Together, our results support the low risk of tumorigenesis associated with AAV-mediated gene transfer targeting juvenile/young adult livers, although constant monitoring of subjects enrolled in AAV clinical trial is advisable.

INSPIIRED: Quantification and Visualization Tools for Analyzing Integration Site Distributions

Analysis of sites of newly integrated DNA in cellular genomes is important to several fields, but methods for analyzing and visualizing these datasets are still under development. Here, we describe tools for data analysis and visualization that take as input integration site data from our INSPIIRED pipeline. Paired-end sequencing allows inference of the numbers of transduced cells as well as the distributions of integration sites in target genomes. We present interactive heatmaps that allow comparison of distributions of integration sites to genomic features and that support numerous user-defined statistical tests. To summarize integration site data from human gene therapy samples, we developed a reproducible report format that catalogs sample population structure, longitudinal dynamics, and integration frequency near cancer-associated genes. We also introduce a novel summary statistic, the UC50 (unique cell progenitors contributing the most expanded 50% of progeny cell clones), which provides a single number summarizing possible clonal expansion. Using these tools, we characterize ongoing longitudinal characterization of a patient from the first trial to treat severe combined immunodeficiency-X1 (SCID-X1), showing successful reconstitution for 15 years accompanied by persistence of a cell clone with an integration site near the cancer-associated gene CCND2. Software is available at https://github.com/BushmanLab/INSPIIRED.

INSPIIRED: A Pipeline for Quantitative Analysis of Sites of New DNA Integration in Cellular Genomes

Integration of new DNA into cellular genomes mediates replication of retroviruses and transposons; integration reactions have also been adapted for use in human gene therapy. Tracking the distributions of integration sites is important to characterize populations of transduced cells and to monitor potential outgrow of pathogenic cell clones. Here, we describe a pipeline for quantitative analysis of integration site distributions named INSPIIRED (integration site pipeline for paired-end reads). We describe optimized biochemical steps for site isolation using Illumina paired-end sequencing, including new technology for suppressing recovery of unwanted contaminants, then software for alignment, quality control, and management of integration site sequences. During library preparation, DNAs are broken by sonication, so that after ligation-mediated PCR the number of ligation junction sites can be used to infer abundance of gene-modified cells. We generated integration sites of known positions in silico, and we describe optimization of sample processing parameters refined by comparison to truth. We also present a novel graph-theory-based method for quantifying integration sites in repeated sequences, and we characterize the consequences using synthetic and experimental data. In an accompanying paper, we describe an additional set of statistical tools for data analysis and visualization. Software is available at https://github.com/BushmanLab/INSPIIRED.

Amplification, Next-generation Sequencing, and Genomic DNA Mapping of Retroviral Integration Sites

Retroviruses exhibit signature integration preferences on both the local and global scales. Here, we present a detailed protocol for (1) generation of diverse libraries of retroviral integration sites using ligation-mediated PCR (LM-PCR) amplification and next-generation sequencing (NGS), (2) mapping the genomic location of each virus-host junction using BEDTools, and (3) analyzing the data for statistical relevance. Genomic DNA extracted from infected cells is fragmented by digestion with restriction enzymes or by sonication. After suitable DNA end-repair, double-stranded linkers are ligated onto the DNA ends, and semi-nested PCR is conducted using primers complementary to both the long terminal repeat (LTR) end of the virus and the ligated linker DNA. The PCR primers carry sequences required for DNA clustering during NGS, negating the requirement for separate adapter ligation. Quality control (QC) is conducted to assess DNA fragment size distribution and adapter DNA incorporation prior to NGS. Sequence output files are filtered for LTR-containing reads, and the sequences defining the LTR and the linker are cropped away. Trimmed host cell sequences are mapped to a reference genome using BLAT and are filtered for minimally 97% identity to a unique point in the reference genome. Unique integration sites are scrutinized for adjacent nucleotide (nt) sequence and distribution relative to various genomic features. Using this protocol, integration site libraries of high complexity can be constructed from genomic DNA in three days. The entire protocol that encompasses exogenous viral infection of susceptible tissue culture cells to integration site analysis can therefore be conducted in approximately one to two weeks. Recent applications of this technology pertain to longitudinal analysis of integration sites from HIV-infected patients.

A 3D Searchable Database of Transgenic Zebrafish Gal4 and Cre Lines for Functional Neuroanatomy Studies

Transgenic methods enable the selective manipulation of neurons for functional mapping of neuronal circuits. Using confocal microscopy, we have imaged the cellular-level expression of 109 transgenic lines in live 6 day post fertilization larvae, including 80 Gal4 enhancer trap lines, 9 Cre enhancer trap lines and 20 transgenic lines that express fluorescent proteins in defined gene-specific patterns. Image stacks were acquired at single micron resolution, together with a broadly expressed neural marker, which we used to align enhancer trap reporter patterns into a common 3-dimensional reference space. To facilitate use of this resource, we have written software that enables searching for transgenic lines that label cells within a selectable 3-dimensional region of interest (ROI) or neuroanatomical area. This software also enables the intersectional expression of transgenes to be predicted, a feature which we validated by detecting cells with co-expression of Cre and Gal4. Many of the imaged enhancer trap lines show intrinsic brain-specific expression. However, to increase the utility of lines that also drive expression in non-neuronal tissue we have designed a novel UAS reporter, that suppresses expression in heart, muscle, and skin through the incorporation of microRNA binding sites in a synthetic 3′ untranslated region. Finally, we mapped the site of transgene integration, thus providing molecular identification of the expression pattern for most lines. Cumulatively, this library of enhancer trap lines provides genetic access to 70% of the larval brain and is therefore a powerful and broadly accessible tool for the dissection of neural circuits in larval zebrafish.