Development and validation of a new high-throughput method to investigate the clonality of HTLV-1-infected cells based on provirus integration sites

A macrophage-cell model of HIV latency reveals the unusual importance of the bromodomain axis

BACKGROUND

Although macrophages are now recognized as an essential part of the HIV latent reservoir, whether and how viral latency is established and reactivated in these cell types is poorly understood. To understand the fundamental mechanisms of viral latency in macrophages, there is an urgent need to develop latency models amenable to genetic manipulations and screening for appropriate latency-reversing agents (LRAs). Given that differentiated THP-1 cells resemble monocyte-derived macrophages in HIV replication mechanisms, we set out to establish a macrophage cell model for HIV latency using THP-1 cells.

METHODS

We created single-cell clones of THP-1 cells infected with a single copy of the dual-labeled HIVGKO in which a codon switched eGFP (csGFP) is under the control of the HIV-1 5' LTR promoter, and a monomeric Kusabira orange 2 (mKO2) under the control of cellular elongation factor one alpha promoter (EF1α). Latently infected cells are csGFP-, mKO2+, while cells with actively replicating HIV (or reactivated virus) are csGFP+,mKO2+. After sorting for latently infected cells, each of the THP-1 clones with unique integration sites for HIV was differentiated into macrophage-like cells with phorbol 12-myristate 13-acetate (PMA) and treated with established LRAs to stimulate HIV reactivation. Monocyte-derived macrophages (MDMs) harboring single copies of HIVGKO were used to confirm our findings.

RESULTS

We obtained clones of THP-1 cells with latently infected HIV with unique integration sites. When the differentiated THP-1 or primary MDMs cells were treated with various LRAs, the bromodomain inhibitors JQ1 and I-BET151 were the most potent compounds. Knockdown of BRD4, the target of JQ1, resulted in increased reactivation, thus confirming the pharmacological effect. The DYRK1A inhibitor Harmine and lipopolysaccharide (LPS) also showed significant reactivation across all three MDM donors. Remarkably, LRAs like PMA/ionomycin, bryostatin-1, and histone deacetylase inhibitors known to potently reactivate latent HIV in CD4 + T cells showed little activity in macrophages.

CONCLUSIONS

Our results indicate that this model could be used to screen for appropriate LRAs for macrophages and show that HIV latency and reactivation mechanisms in macrophages may be distinct from those of CD4 + T cells.

Determination of molecular epidemiologic pattern of human T-lymphotropic virus type 1 (HTLV-1) in Alborz province, Iran

Human T-cell lymphotropic virus type 1 (HTLV-1) is linked to two debilitating diseases, adult T-cell leukemia/lymphoma (ATLL) and HTLV-1 associated myelopathy tropical spastic paraparesis (HAM/TSP), which are prevalent in various parts of the world, including the Alborz province in Iran. Understanding the prevalence and evolutionary relationships of HTLV-1 infections in these endemic areas is of utmost importance. In the realm of phylogenetic studies, long terminal repeat (LTR) region of HTLV-1 stands out as highly conserved, yet more variable compared to other gene segments. Consequently, it is the primary focus for phylogenetic analyses. Additionally, trans-activator of transcription (Tax), an oncoprotein, holds a pivotal role in the regulation of gene expression. This cross-sectional study delved into the phylogenetic analysis of HTLV-1 among individuals in Alborz province of Iran. To confirm infection, we amplified partial sequence LTR (PLTR) and HTLV-1 bZIP factor (PHBZ). For phylogenetic analysis, we sequenced the full sequence LTR (FLTR) and full Tax sequence (FTax). The FLTR and FTax sequences underwent analysis using BioEdit, and phylogenetic trees were constructed using MEGA-X software. Out of the roughly 15,000 annual blood donors in Alborz, 19 samples tested positive for HTLV-1, indicating a 0.13% HTLV-1 positivity rate among blood donors. Furthermore, the HTLV-1 virus prevalent in the Alborz province belongs to subtype A (cosmopolitan) subgroup A. The findings revealed that while mutations were observed in both the LTR and Tax genes, they were not significant enough to bring about fundamental alterations. Despite positive selection detected in three Alborz isolates, it has not led to mutations affecting Tax function and virulence.

Preclinical assessment of an anti-HTLV-1 heterologous DNA/MVA vaccine protocol expressing a multiepitope HBZ protein

BACKGROUND

Human T-lymphotropic virus 1 (HTLV-1) is associated with the development of several pathologies and chronic infection in humans. The inefficiency of the available treatments and the challenge in developing a protective vaccine highlight the need to produce effective immunotherapeutic tools. The HTLV-1 basic leucine zipper (bZIP) factor (HBZ) plays an important role in the HTLV-1 persistence, conferring a survival advantage to infected cells by reducing the HTLV-1 proteins expression, allowing infected cells to evade immune surveillance, and enhancing cell proliferation leading to increased proviral load.

METHODS

We have generated a recombinant Modified Virus Vaccinia Ankara (MVA-HBZ) and a plasmid DNA (pcDNA3.1(+)-HBZ) expressing a multiepitope protein based on peptides of HBZ to study the immunogenic potential of this viral-derived protein in BALB/c mice model. Mice were immunized in a prime-boost heterologous protocol and their splenocytes (T CD4+ and T CD8+) were immunophenotyped by flow cytometry and the humoral response was evaluated by ELISA using HBZ protein produced in prokaryotic vector as antigen.

RESULTS

T CD4+ and T CD8+ lymphocytes cells stimulated by HBZ-peptides (HBZ42-50 and HBZ157-176) showed polyfunctional double positive responses for TNF-α/IFN-γ, and TNF-α/IL-2. Moreover, T CD8+ cells presented a tendency in the activation of effector memory cells producing granzyme B (CD44+High/CD62L-Low), and the activation of Cytotoxic T Lymphocytes (CTLs) and cytotoxic responses in immunized mice were inferred through the production of granzyme B by effector memory T cells and the expression of CD107a by CD8+ T cells. The overall data is consistent with a directive and effector recall response, which may be able to operate actively in the elimination of HTLV-1-infected cells and, consequently, in the reduction of the proviral load. Sera from immunized mice, differently from those of control animals, showed IgG-anti-HBZ production by ELISA.

CONCLUSIONS

Our results highlight the potential of the HBZ multiepitope protein expressed from plasmid DNA and a poxviral vector as candidates for therapeutic vaccine.

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.

Hematopoietic reconstitution dynamics of mobilized- and bone marrow-derived human hematopoietic stem cells after gene therapy

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.

Meeting FDA Guidance recommendations for replication-competent virus and insertional oncogenesis testing

Integrating vectors are associated with alterations in cellular function related to disruption of normal gene function. This has been associated with clonal expansion of cells and, in some instances, cancer. These events have been associated with replication-defective vectors and suggest that the inadvertent exposure to a replication-competent virus arising during vector manufacture would significantly increase the risk of treatment-related adverse events. These risks have led regulatory agencies to require specific monitoring for replication-competent viruses, both prior to and after treatment of patients with gene therapy products. Monitoring the risk of cell expansion and malignancy is also required. In this review, we discuss the rational potential approaches and challenges to meeting the US FDA expectations listed in current guidance documents.

Longitudinal single-cell profiling of chemotherapy response in acute myeloid leukemia

Acute myeloid leukemia may be characterized by a fraction of leukemia stem cells (LSCs) that sustain disease propagation eventually leading to relapse. Yet, the contribution of LSCs to early therapy resistance and AML regeneration remains controversial. We prospectively identify LSCs in AML patients and xenografts by single-cell RNA sequencing coupled with functional validation by a microRNA-126 reporter enriching for LSCs. Through nucleophosmin 1 (NPM1) mutation calling or chromosomal monosomy detection in single-cell transcriptomes, we discriminate LSCs from regenerating hematopoiesis, and assess their longitudinal response to chemotherapy. Chemotherapy induced a generalized inflammatory and senescence-associated response. Moreover, we observe heterogeneity within progenitor AML cells, some of which proliferate and differentiate with expression of oxidative-phosphorylation (OxPhos) signatures, while others are OxPhos (low) miR-126 (high) and display enforced stemness and quiescence features. miR-126 (high) LSCs are enriched at diagnosis in chemotherapy-refractory AML and at relapse, and their transcriptional signature robustly stratifies patients for survival in large AML cohorts.

Diagnosis and Early Prediction of Lymphoma Using High-Throughput Clonality Analysis of Bovine Leukemia Virus-Infected Cells

Bovine leukemia virus (BLV), a retrovirus, infects B cells of ruminants and is integrated into the host genome as a provirus for lifelong infection. After a long latent period, 1% to 5% of BLV-infected cattle develop aggressive lymphoma, enzootic bovine leukosis (EBL). Since the clonal expansion of BLV-infected cells is essential for the development of EBL, the clonality of proviral integration sites could be a molecular marker for diagnosis and early prediction of EBL. Recently, we developed Rapid Amplification of the Integration Site without Interference by Genomic DNA Contamination (RAISING) and an analysis software of clonality value (CLOVA) to analyze the clonality of transgene-integrated cells. RAISING-CLOVA is capable of assessing the risk of adult T-cell leukemia/lymphoma development in human T-cell leukemia virus-I-infected individuals through the clonality analysis of proviral integration sites. Thus, we herein examined the performance of RAISING-CLOVA for the clonality analysis of BLV-infected cells and conducted a comprehensive clonality analysis by RAISING-CLOVA in EBL and non-EBL cattle. RAISING-CLOVA targeting BLV was a highly accurate and reproducible method for measuring the clonality value. The comprehensive clonality analysis successfully distinguished EBL from non-EBL specimens with high sensitivity and specificity. A longitudinal clonality analysis in BLV-infected sheep, an experimental model of lymphoma, also confirmed the effectiveness of RAISING-CLOVA for early detection of EBL development. Therefore, our study emphasizes the usefulness of RAISING-CLOVA as a routine clinical test for monitoring virus-related cancers. IMPORTANCE Bovine leukemia virus (BLV) infection causes aggressive B-cell lymphoma in cattle and sheep. The virus has spread to farms around the world, causing significant economic damage to the livestock industry. Thus, the identification of high-risk asymptomatic cattle before they develop lymphoma can be effective in reducing the economic damage. Clonal expansion of BLV-infected cells is a promising marker for the development of lymphoma. Recently, we have developed a high-throughput method to amplify random integration sites of transgenes in host genomes and analyze their clonality, named as RAISING-CLOVA. As a new application of our technology, in this study, we demonstrate the value of the RAISING-CLOVA method for the diagnosis and early prediction of lymphoma development by BLV infection in cattle. RAISING-CLOVA is a reliable technology for monitoring the clonality of BLV-infected cells and would contribute to reduce the economic losses by EBL development.

ISAnalytics enables longitudinal and high-throughput clonal tracking studies in hematopoietic stem cell gene therapy applications

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.

Clonal reconstruction from co-occurrence of vector integration sites accurately quantifies expanding clones in vivo

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.

High transduction rates of viral vectors ensure good gene delivery; however multiple integration events can occur in the same cell. Here the authors use correlations between repeated measurements of integration site abundances to estimate their mutual similarity and identify clusters of co-occurring sites.

RAISING is a high-performance method for identifying random transgene integration sites

Both natural viral infections and therapeutic interventions using viral vectors pose significant risks of malignant transformation. Monitoring for clonal expansion of infected cells is important for detecting cancer. Here we developed a novel method of tracking clonality via the detection of transgene integration sites. RAISING (Rapid Amplification of Integration Sites without Interference by Genomic DNA contamination) is a sensitive, inexpensive alternative to established methods. Its compatibility with Sanger sequencing combined with our CLOVA (Clonality Value) software is critical for those without access to expensive high throughput sequencing. We analyzed samples from 688 individuals infected with the retrovirus HTLV-1, which causes adult T-cell leukemia/lymphoma (ATL) to model our method. We defined a clonality value identifying ATL patients with 100% sensitivity and 94.8% specificity, and our longitudinal analysis also demonstrates the usefulness of ATL risk assessment. Future studies will confirm the broad applicability of our technology, especially in the emerging gene therapy sector.

Elucidation of the Mechanism of Host NMD Suppression by HTLV-1 Rex: Dissection of Rex to Identify the NMD Inhibitory Domain

The human retrovirus human T-cell leukemia virus type I (HTLV-1) infects human T cells by vertical transmission from mother to child through breast milk or horizontal transmission through blood transfusion or sexual contact. Approximately 5% of infected individuals develop adult T-cell leukemia/lymphoma (ATL) with a poor prognosis, while 95% of infected individuals remain asymptomatic for the rest of their lives, during which time the infected cells maintain a stable immortalized latent state in the body. It is not known why such a long latent state is maintained. We hypothesize that the role of functional proteins of HTLV-1 during early infection influences the phenotype of infected cells in latency. In eukaryotic cells, a mRNA quality control mechanism called nonsense-mediated mRNA decay (NMD) functions not only to eliminate abnormal mRNAs with nonsense codons but also to target virus-derived RNAs. We have reported that HTLV-1 genomic RNA is a potential target of NMD, and that Rex suppresses NMD and stabilizes viral RNA against it. In this study, we aimed to elucidate the molecular mechanism of NMD suppression by Rex using various Rex mutant proteins. We found that region X (aa20-57) of Rex, the function of which has not been clarified, is required for NMD repression. We showed that Rex binds to Upf1, which is the host key regulator to detect abnormal mRNA and initiate NMD, through this region. Rex also interacts with SMG5 and SMG7, which play essential roles for the completion of the NMD pathway. Moreover, Rex selectively binds to Upf3B, which is involved in the normal NMD complex, and replaces it with a less active form, Upf3A, to reduce NMD activity. These results revealed that Rex invades the NMD cascade from its initiation to completion and suppresses host NMD activity to protect the viral genomic mRNA.

Efficient and safe correction of hemophilia A by lentiviral vector-transduced BOECs in an implantable device

Hemophilia A (HA) is a rare bleeding disorder caused by deficiency/dysfunction of the FVIII protein. As current therapies based on frequent FVIII infusions are not a definitive cure, long-term expression of FVIII in endothelial cells through lentiviral vector (LV)-mediated gene transfer holds the promise of a one-time treatment. Thus, here we sought to determine whether LV-corrected blood outgrowth endothelial cells (BOECs) implanted through a prevascularized medical device (Cell Pouch) would rescue the bleeding phenotype of HA mice. To this end, BOECs from HA patients and healthy donors were isolated, expanded, and transduced with an LV carrying FVIII driven by an endothelial-specific promoter employing GMP-like procedures. FVIII-corrected HA BOECs were either directly transplanted into the peritoneal cavity or injected into a Cell Pouch implanted subcutaneously in NSG-HA mice. In both cases, FVIII secretion was sufficient to improve the mouse bleeding phenotype. Indeed, FVIII-corrected HA BOECs reached a relatively short-term clinically relevant engraftment being detected up to 16 weeks after transplantation, and their genomic integration profile did not show enrichment for oncogenes, confirming the process safety. Overall, this is the first preclinical study showing the safety and feasibility of transplantation of GMP-like produced LV-corrected BOECs within an implantable device for the long-term treatment of HA.

Retrieval of vector integration sites from cell-free DNA

Gene therapy (GT) has rapidly attracted renewed interest as a treatment for otherwise incurable diseases, with several GT products already on the market and many more entering clinical testing for selected indications. Clonal tracking techniques based on vector integration enable monitoring of the fate of engineered cells in the blood of patients receiving GT and allow assessment of the safety and efficacy of these procedures. However, owing to the limited number of cells that can be tested and the impracticality of studying cells residing in peripheral organs without performing invasive biopsies, this approach provides only a partial snapshot of the clonal repertoire and dynamics of genetically modified cells and reduces the predictive power as a safety readout. In this study, we developed liquid biopsy integration site sequencing, or LiBIS-seq, a polymerase chain reaction technique optimized to quantitatively retrieve vector integration sites from cell-free DNA released into the bloodstream by dying cells residing in several tissues. This approach enabled longitudinal monitoring of in vivo liver-directed GT and clonal tracking in patients receiving hematopoietic stem cell GT, improving our understanding of the clonal composition and turnover of genetically modified cells in solid tissues and, in contrast to conventional analyses based only on circulating blood cells, enabling earlier detection of vector-marked clones that are aberrantly expanding in peripheral tissues.

Epigenomic regulation of human T-cell leukemia virus by chromatin-insulator CTCF

Human T-cell leukemia virus type 1 (HTLV-1) is a retrovirus that causes an aggressive T-cell malignancy and a variety of inflammatory conditions. The integrated provirus includes a single binding site for the epigenomic insulator, CCCTC-binding protein (CTCF), but its function remains unclear. In the current study, a mutant virus was examined that eliminates the CTCF-binding site. The mutation did not disrupt the kinetics and levels of virus gene expression, or establishment of or reactivation from latency. However, the mutation disrupted the epigenetic barrier function, resulting in enhanced DNA CpG methylation downstream of the CTCF binding site on both strands of the integrated provirus and H3K4Me3, H3K36Me3, and H3K27Me3 chromatin modifications both up- and downstream of the site. A majority of clonal cell lines infected with wild type HTLV-1 exhibited increased plus strand gene expression with CTCF knockdown, while expression in mutant HTLV-1 clonal lines was unaffected. These findings indicate that CTCF binding regulates HTLV-1 gene expression, DNA and histone methylation in an integration site dependent fashion.

Human T-cell leukemia virus type 1 (HTLV-1) is a cause of leukemia and lymphoma as well as several inflammatory medical disorders. The virus integrates in the host cell DNA, and it has a single binding site for a protein designated CTCF. This protein is important in the regulation of many DNA viruses, as well as many properties of normal and malignant cells. In order to define the role of CTCF binding to HTLV, we analyzed a mutant virus lacking the binding site. We found that this mutation variably affected gene expression, DNA and histone modification, suggesting a key role in regulation of virus replication in infected cells.

Quantification of T cell clonality in human T cell leukaemia virus type-1 carriers can detect the development of adult T cell leukaemia early

Adult T cell leukaemia/lymphoma (ATL) arises from clonally expanded T cells that are infected with human T cell leukaemia virus type-1 (HTLV-1). Here, we show that ATL can be detected early in HTLV-1-carriers through quantification of T-cell receptor (TCR)Vβ subunit diversity on T-cells infected with HTLV-1 (CD3+ CCR4+ CD26- T-cells) using an 'oligoclonality index' (OCI-flow). We established a reference range for OCI-flow by analysing peripheral blood mononuclear cells (PBMCs) from HTLV-1-carriers who had not developed ATL in a median of 10.5 years follow up (n = 38) and patients with ATL (n = 30). In the third cohort of HTLV-1-carriers with no history or clinical evidence of ATL (n = 106), 19% of high proviral load (PVL, ≥4 copies of HTLV-1/100 PBMCs) carriers had an OCI-flow in the ATL range, >0.770. Carriers with an OCI-flow >0.770 (n = 14) had higher lymphocyte counts and PVLs and were more likely to have a family history of ATL than carriers with OCI-flow ≤0.770. ATL subsequently developed in two of these 14 carriers but no carriers with OCI-flow ≤0.770 (p = 0.03, cumulative follow-up 129 person-years). This method can be used to identify a subset of high-PVL HTLV-1-carriers at increased risk of developing ATL who may benefit from intervention therapy, prior to the detection of disease.

Overexpression of aberrant Wnt5a and its effect on acquisition of malignant phenotypes in adult T-cell leukemia/lymphoma (ATL) cells

Wnt5a is a ligand of the non-canonical Wnt signaling pathway involved in cell differentiation, motility, and inflammatory response. Adult T-cell leukemia/lymphoma (ATL) is one of the most aggressive T-cell malignancies caused by infection of human T-cell leukemia virus type1 (HTLV-1). Among subtypes of ATL, acute-type ATL cells are particularly resistant to current multidrug chemotherapies and show remarkably high cell-proliferative and invasive phenotypes. Here we show a dramatic increase of WNT5A gene expression in acute-type ATL cells compared with those of indolent-type ATL cells. Treatment with IWP-2 or Wnt5a-specific knockdown significantly suppressed cell growth of ATL-derived T-cell lines. We demonstrated that the overexpression of c-Myb and FoxM1 was responsible for the synergistic activation of the WNT5A promoter. Also, a WNT5A transcript variant without the exon4 (the ΔE4-WNT5A mRNA), encoding ΔC-Wnt5 (1-136aa of 380aa), is overexpressed in acute-type ATL cells. The ΔC-Wnt5a is secreted extracellularly and enhances cellular migration/invasion to a greater extent compared with wildtype (WT)-Wnt5a. Moreover, the ΔC-Wnt5a secretion was not suppressed by IWP-2, indicating that this mutant Wnt5a is secreted via a different pathway from the WT-Wnt5a. Taken together, synergistic overexpression of the ΔC-Wnt5a by c-Myb and FoxM1 may be responsible for the malignant phenotype of acute-type ATL cells.

A Methodology for Assessing Tumor Clonality of Adult T Cell Leukemia/Lymphoma

While clonal heterogeneity has been demonstrated in most cancers, quantitative assessment of individual tumor clones has not been translated to inform clinical practice. A few methods have been developed to investigate the tumor clonality of adult T cell leukemia/lymphoma (ATLL), but currently there is no clinically translatable method available for quantifying individual tumor clones in ATLL patients. Here, we present a methodology to assess the tumor clonality of ATLL and quantify patient-specific tumor clones in a clinical setting. The methodology consists of three steps: (1) selective amplification of restriction fragments containing a human T cell leukemia virus type 1 (HTLV-1) integration site, (2) amplicon deep sequencing to estimate the clonal structure and identify HTLV-1 integration sites of dominant clones, and (3) digital PCR targeting the HTLV-1 integration sites of the dominant clones to quantify specific tumor clones. We successfully tracked individual tumor clones using this approach and demonstrated that each clone had a distinct response to therapies. The procedure is straightforward and clinically feasible, which should facilitate the proper assessment and management of ATLL.

Sleeping Beauty-engineered CAR T cells achieve antileukemic activity without severe toxicities

BACKGROUNDChimeric antigen receptor (CAR) T cell immunotherapy has resulted in complete remission (CR) and durable response in highly refractory patients. However, logistical complexity and high costs of manufacturing autologous viral products limit CAR T cell availability.METHODSWe report the early results of a phase I/II trial in B cell acute lymphoblastic leukemia (B-ALL) patients relapsed after allogeneic hematopoietic stem cell transplantation (HSCT) using donor-derived CD19 CAR T cells generated with the Sleeping Beauty (SB) transposon and differentiated into cytokine-induced killer (CIK) cells.RESULTSThe cellular product was produced successfully for all patients from the donor peripheral blood (PB) and consisted mostly of CD3+ lymphocytes with 43% CAR expression. Four pediatric and 9 adult patients were infused with a single dose of CAR T cells. Toxicities reported were 2 grade I and 1 grade II cytokine-release syndrome (CRS) cases at the highest dose in the absence of graft-versus-host disease (GVHD), neurotoxicity, or dose-limiting toxicities. Six out of 7 patients receiving the highest doses achieved CR and CR with incomplete blood count recovery (CRi) at day 28. Five out of 6 patients in CR were also minimal residual disease negative (MRD-). Robust expansion was achieved in the majority of the patients. CAR T cells were measurable by transgene copy PCR up to 10 months. Integration site analysis showed a positive safety profile and highly polyclonal repertoire in vitro and at early time points after infusion.CONCLUSIONSB-engineered CAR T cells expand and persist in pediatric and adult B-ALL patients relapsed after HSCT. Antileukemic activity was achieved without severe toxicities.TRIAL REGISTRATIONClinicalTrials.gov NCT03389035.FUNDINGThis study was supported by grants from the Fondazione AIRC per la Ricerca sul Cancro (AIRC); Cancer Research UK (CRUK); the Fundación Científica de la Asociación Española Contra el Cáncer (FC AECC); Ministero Della Salute; Fondazione Regionale per la Ricerca Biomedica (FRRB).

An Improved Sequencing-Based Bioinformatics Pipeline to Track the Distribution and Clonal Architecture of Proviral Integration Sites

The combined application of linear amplification-mediated PCR (LAM-PCR) protocols with next-generation sequencing (NGS) has had a large impact on our understanding of retroviral pathogenesis. Previously, considerable effort has been expended to optimize NGS methods to explore the genome-wide distribution of proviral integration sites and the clonal architecture of clinically important retroviruses like human T-cell leukemia virus type-1 (HTLV-1). Once sequencing data are generated, the application of rigorous bioinformatics analysis is central to the biological interpretation of the data. To better exploit the potential information available through these methods, we developed an optimized bioinformatics pipeline to analyze NGS clonality datasets. We found that short-read aligners, specifically designed to manage NGS datasets, provide increased speed, significantly reducing processing time and decreasing the computational burden. This is achieved while also accounting for sequencing base quality. We demonstrate the utility of an additional trimming step in the workflow, which adjusts for the number of reads supporting each insertion site. In addition, we developed a recall procedure to reduce bias associated with proviral integration within low complexity regions of the genome, providing a more accurate estimation of clone abundance. Finally, we recommend the application of a “clean-and-recover” step to clonality datasets generated from large cohorts and longitudinal studies. In summary, we report an optimized bioinformatics workflow for NGS clonality analysis and describe a new set of steps to guide the computational process. We demonstrate that the application of this protocol to the analysis of HTLV-1 and bovine leukemia virus (BLV) clonality datasets improves the quality of data processing and provides a more accurate definition of the clonal landscape in infected individuals. The optimized workflow and analysis recommendations can be implemented in the majority of bioinformatics pipelines developed to analyze LAM-PCR-based NGS clonality datasets.

Hematopoietic Tumors in a Mouse Model of X-linked Chronic Granulomatous Disease after Lentiviral Vector-Mediated Gene Therapy

Chronic granulomatous disease (CGD) is a rare inherited disorder due to loss-of-function mutations in genes encoding the NADPH oxidase subunits. Hematopoietic stem and progenitor cell (HSPC) gene therapy (GT) using regulated lentiviral vectors (LVs) has emerged as a promising therapeutic option for CGD patients. We performed non-clinical Good Laboratory Practice (GLP) and laboratory-grade studies to assess the safety and genotoxicity of LV targeting myeloid-specific Gp91^phox expression in X-linked chronic granulomatous disease (XCGD) mice. We found persistence of gene-corrected cells for up to 1 year, restoration of Gp91^phox expression and NADPH oxidase activity in XCGD phagocytes, and reduced tissue inflammation after LV-mediated HSPC GT. Although most of the mice showed no hematological or biochemical toxicity, a small subset of XCGD GT mice developed T cell lymphoblastic lymphoma (2.94%) and myeloid leukemia (5.88%). No hematological malignancies were identified in C57BL/6 mice transplanted with transduced XCGD HSPCs. Integration pattern analysis revealed an oligoclonal composition with rare dominant clones harboring vector insertions near oncogenes in mice with tumors. Collectively, our data support the long-term efficacy of LV-mediated HSPC GT in XCGD mice and provide a safety warning because the chronic inflammatory XCGD background may contribute to oncogenesis.

The Nature of the HTLV-1 Provirus in Naturally Infected Individuals Analyzed by the Viral DNA-Capture-Seq Approach

The retrovirus human T-cell leukemia virus type 1 (HTLV-1) integrates into the host DNA, achieves persistent infection, and induces human diseases. Here, we demonstrate that viral DNA-capture sequencing (DNA-capture-seq) is useful to characterize HTLV-1 proviruses in naturally virus-infected individuals, providing comprehensive information about the proviral structure and the viral integration site. We analyzed peripheral blood from 98 naturally HTLV-1-infected individuals and found that defective proviruses were present not only in patients with leukemia, but also in those with other clinical entities. We further demonstrated that clones with defective-type proviruses exhibited a higher degree of clonal abundance than those with full-length proviruses. The frequency of defective-type proviruses in HTLV-1-infected humanized mice was lower than that in infected individuals, indicating that defective proviruses were rare at the initial phase of infection but preferentially selected during persistent infection. These results demonstrate the robustness of viral DNA-capture-seq for HTLV-1 infection and suggest potential applications for other virus-associated cancers in humans.

A high-throughput detection method for the clonality of Human T-cell leukemia virus type-1-infected cells in vivo

Approximately 10-20 million of Human T-cell leukemia virus type-1 (HTLV-1)-infected carriers have been previously reported, and approximately 5% of these carriers develop adult T-cell leukemia/lymphoma (ATL) with a characteristic poor prognosis. In Japan, Southern blotting has long been routinely performed for detection of clonally expanded ATL cells in vivo, and as a confirmatory diagnostic test for ATL. However, alternative methods to Southern blotting, such as sensitive, quantitative, and rapid analytical methods, are currently required in clinical practice. In this study, we developed a high-throughput method called rapid amplification of integration site (RAIS) that could amplify HTLV-1-integrated fragments within 4 h and detect the integration sites in > 0.16% of infected cells. Furthermore, we established a novel quantification method for HTLV-1 clonality using Sanger sequencing with RAIS products, and the validity of the quantification method was confirmed by comparing it with next-generation sequencing in terms of the clonality. Thus, we believe that RAIS has a high potential for use as an alternative routine molecular confirmatory test for the clonality analysis of HTLV-1-infected cells.

Phagocytosis-shielded lentiviral vectors improve liver gene therapy in nonhuman primates

Liver-directed gene therapy for the coagulation disorder hemophilia showed safe and effective results in clinical trials using adeno-associated viral vectors to replace a functional coagulation factor, although some unmet needs remain. Lentiviral vectors (LVs) may address some of these hurdles because of their potential for stable expression and the low prevalence of preexisting viral immunity in humans. However, systemic LV administration to hemophilic dogs was associated to mild acute toxicity and low efficacy at the administered doses. Here, exploiting intravital microscopy and LV surface engineering, we report a major role of the human phagocytosis inhibitor CD47, incorporated into LV cell membrane, in protecting LVs from uptake by professional phagocytes and innate immune sensing, thus favoring biodistribution to hepatocytes after systemic administration. By enforcing high CD47 surface content, we generated phagocytosis-shielded LVs which, upon intravenous administration to nonhuman primates, showed selective liver and spleen targeting and enhanced hepatocyte gene transfer compared to parental LV, reaching supraphysiological activity of human coagulation factor IX, the protein encoded by the transgene, without signs of toxicity or clonal expansion of transduced cells.

Intrathymic adeno-associated virus gene transfer rapidly restores thymic function and long-term persistence of gene-corrected T cells

BACKGROUND

Patients with T-cell immunodeficiencies are generally treated with allogeneic hematopoietic stem cell transplantation, but alternatives are needed for patients without matched donors. An innovative intrathymic gene therapy approach that directly targets the thymus might improve outcomes.

OBJECTIVE

We sought to determine the efficacy of intrathymic adeno-associated virus (AAV) serotypes to transduce thymocyte subsets and correct the T-cell immunodeficiency in a zeta-associated protein of 70 kDa (ZAP-70)-deficient murine model.

METHODS

AAV serotypes were injected intrathymically into wild-type mice, and gene transfer efficiency was monitored. ZAP-70^-/- mice were intrathymically injected with an AAV8 vector harboring the ZAP70 gene. Thymus structure, immunophenotyping, T-cell receptor clonotypes, T-cell function, immune responses to transgenes and autoantibodies, vector copy number, and integration were evaluated.

RESULTS

AAV8, AAV9, and AAV10 serotypes all transduced thymocyte subsets after in situ gene transfer, with transduction of up to 5% of cells. Intrathymic injection of an AAV8-ZAP-70 vector into ZAP-70^-/- mice resulted in a rapid thymocyte differentiation associated with the development of a thymic medulla. Strikingly, medullary thymic epithelial cells expressing the autoimmune regulator were detected within 10 days of gene transfer, correlating with the presence of functional effector and regulatory T-cell subsets with diverse T-cell receptor clonotypes in the periphery. Although thymocyte reconstitution was transient, gene-corrected peripheral T cells harboring approximately 1 AAV genome per cell persisted for more than 40 weeks, and AAV vector integration was detected.

CONCLUSIONS

Intrathymic AAV-transduced progenitors promote a rapid restoration of the thymic architecture, with a single wave of thymopoiesis generating long-term peripheral T-cell function.

Assessing the Impact of Cyclosporin A on Lentiviral Transduction and Preservation of Human Hematopoietic Stem Cells in Clinically Relevant Ex Vivo Gene Therapy Settings

Improving hematopoietic stem and progenitor cell (HSPC) permissiveness to lentiviral vector (LV) transduction without compromising their biological properties remains critical for broad-range implementation of gene therapy as a treatment option for several inherited diseases. This study demonstrates that the use of one-hit ex vivo LV transduction protocols based on either cyclosporin A (CsA) or rapamycin enable as efficient gene transfer as the current two-hit clinical standard into bone marrow-derived CD34⁺ cells while better preserving their engraftment capacity in vivo. CsA was additive with another enhancer of transduction, prostaglandin E2, suggesting that tailored enhancer combinations may be applied to overcome multiple blocks to transduction simultaneously in HSPC. Interestingly, besides enhancing LV transduction, CsA also significantly reduced HSPC proliferation, preserving the quiescent G₀ fraction and the more primitive multipotent progenitors, thereby yielding the highest engraftment levels in vivo. Importantly, no alterations in the vector integration profiles could be detected between CsA and control transduced HSPC. Overall, the present findings contribute to the development of more efficient and sustainable LV gene therapy protocols, underscoring the benefits of scaling down required vector doses, as well as shortening the HSPC ex vivo culture time.

RNA sequencing identifies clonal structure of T-cell repertoires in patients with adult T-cell leukemia/lymphoma

The diversity of T-cell receptor (TCR) repertoires, as generated by somatic DNA rearrangements, is central to immune system function. High-throughput sequencing technologies now allow examination of antigen receptor repertoires at single-nucleotide and, more recently, single-cell resolution. The TCR repertoire can be altered in the context of infections, malignancies or immunological disorders. Here we examined the diversity of TCR clonality and its association with pathogenesis and prognosis in adult T-cell leukemia/lymphoma (ATL), a malignancy caused by infection with human T-cell leukemia virus type-1 (HTLV-1). We analyzed 62 sets of high-throughput RNA sequencing data from 59 samples of HTLV-1−infected individuals—asymptomatic carriers (ACs), smoldering, chronic, acute and lymphoma ATL subtypes—and three uninfected controls to evaluate TCR distribution. Based on these TCR profiles, CD4-positive cells and ACs showed polyclonal patterns, whereas ATL patients showed oligo- or monoclonal patterns (with 446 average clonotypes across samples). Expression of TCRα and TCRβ genes in the dominant clone differed among the samples. ACs, CD4-positive samples and smoldering patients showed significantly higher TCR diversity compared with chronic, acute and lymphoma subtypes. CDR3 sequence length distribution, amino acid conservation and gene usage variability for ATL patients resembled those of peripheral blood cells from ACs and healthy donors. Thus, determining monoclonal architecture and clonal diversity by RNA sequencing might be useful for prognostic purposes and for personalizing ATL diagnosis and assessment of treatments.

Mutational Intratumor Heterogeneity is a Complex and Early Event in the Development of Adult T-cell Leukemia/Lymphoma

The clonal architecture of tumors plays a vital role in their pathogenesis and invasiveness; however, it is not yet clear how this clonality contributes to different malignancies. In this study we sought to address mutational intratumor heterogeneity (ITH) in adult T-cell leukemia/lymphoma (ATL). ATL is a malignancy with an incompletely understood molecular pathogenesis caused by infection with human T-cell leukemia virus type-1 (HTLV-1). To determine the clonal structure through tumor genetic diversity profiles, we investigated 142 whole-exome sequencing data of tumor and matched normal samples from 71 ATL patients. Based on SciClone analysis, the ATL samples showed a wide spectrum of modes over clonal/subclonal frequencies ranging from one to nine clusters. The average number of clusters was six across samples, but the number of clusters differed among different samples. Of these ATL samples, 94% had more than two clusters. Aggressive ATL cases had slightly more clonal clusters than indolent types, indicating the presence of ITH during earlier stages of disease. The known significantly mutated genes in ATL were frequently clustered together and possibly coexisted in the same clone. IRF4, CCR4, TP53, and PLCG1 mutations were almost clustered in subclones with a moderate variant allele frequency (VAF), whereas HLA-B, CARD11, and NOTCH1 mutations were clustered in subclones with lower VAFs. Taken together, these results show that ATL displays a high degree of ITH and a complex subclonal structure. Our findings suggest that clonal/subclonal architecture might be a useful measure for prognostic purposes and personalized assessment of the therapeutic response.

VISPA2: a scalable pipeline for high-throughput identification and annotation of vector integration sites

Background

Bioinformatics tools designed to identify lentiviral or retroviral vector insertion sites in the genome of host cells are used to address the safety and long-term efficacy of hematopoietic stem cell gene therapy applications and to study the clonal dynamics of hematopoietic reconstitution. The increasing number of gene therapy clinical trials combined with the increasing amount of Next Generation Sequencing data, aimed at identifying integration sites, require both highly accurate and efficient computational software able to correctly process “big data” in a reasonable computational time.

Results

Here we present VISPA2 (Vector Integration Site Parallel Analysis, version 2), the latest optimized computational pipeline for integration site identification and analysis with the following features: (1) the sequence analysis for the integration site processing is fully compliant with paired-end reads and includes a sequence quality filter before and after the alignment on the target genome; (2) an heuristic algorithm to reduce false positive integration sites at nucleotide level to reduce the impact of Polymerase Chain Reaction or trimming/alignment artifacts; (3) a classification and annotation module for integration sites; (4) a user friendly web interface as researcher front-end to perform integration site analyses without computational skills; (5) the time speedup of all steps through parallelization (Hadoop free).

Conclusions

We tested VISPA2 performances using simulated and real datasets of lentiviral vector integration sites, previously obtained from patients enrolled in a hematopoietic stem cell gene therapy clinical trial and compared the results with other preexisting tools for integration site analysis. On the computational side, VISPA2 showed a > 6-fold speedup and improved precision and recall metrics (1 and 0.97 respectively) compared to previously developed computational pipelines. These performances indicate that VISPA2 is a fast, reliable and user-friendly tool for integration site analysis, which allows gene therapy integration data to be handled in a cost and time effective fashion. Moreover, the web access of VISPA2 (http://openserver.itb.cnr.it/vispa/) ensures accessibility and ease of usage to researches of a complex analytical tool. We released the source code of VISPA2 in a public repository (https://bitbucket.org/andreacalabria/vispa2).

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-017-1937-9) contains supplementary material, which is available to authorized users.

Human T Cell Leukemia Virus Type 1: Persistence and Pathogenesis

Human T cell leukemia virus type 1 (HTLV-1), also known as human T lymphotropic virus type 1, was the first exogenous human retrovirus discovered. Unlike the distantly related lentivirus HIV-1, HTLV-1 causes disease in only 5-10% of infected people, depending on their ethnic origin. But whereas HIV-1 infection and the consequent diseases can be efficiently contained in most cases by antiretroviral drug treatment, there is no satisfactory treatment for the malignant or inflammatory diseases caused by HTLV-1. The purpose of the present article is to review recent advances in the understanding of the mechanisms by which the virus persists in vivo and causes disabling or fatal diseases.

Analyzing the Genotoxicity of Retroviral Vectors in Hematopoietic Cell Gene Therapy

Retroviral vectors, including those derived from gammaretroviruses and lentiviruses, have found their way into the clinical arena and demonstrated remarkable efficacy for the treatment of immunodeficiencies, leukodystrophies, and globinopathies. Despite these successes, gene therapy unfortunately also has had to face severe adverse events in the form of leukemias and myelodysplastic syndromes, related to the semi-random vector integration into the host cell genome that caused deregulation of neighboring proto-oncogenes. Although improvements in vector design clearly lowered the risk of this insertional mutagenesis, analysis of potential genotoxicity and the consequences of vector integration remain important parameters for basic and translational research and most importantly for the clinic. Here, we review current assays to analyze biodistribution and genotoxicity in the pre-clinical setting and describe tools to monitor vector integration sites in vector-treated patients as a biosafety readout.

Inferring clonal structure in HTLV-1-infected individuals: towards bridging the gap between analysis and visualization

Background

Human T cell leukemia virus type 1 (HTLV-1) causes adult T cell leukemia (ATL) in a proportion of infected individuals after a long latency period. Development of ATL is a multistep clonal process that can be investigated by monitoring the clonal expansion of HTLV-1-infected cells by isolation of provirus integration sites. The clonal composition (size, number, and combinations of clones) during the latency period in a given infected individual has not been clearly elucidated.

Methods

We used high-throughput sequencing technology coupled with a tag system for isolating integration sites and measuring clone sizes from 60 clinical samples. We assessed the role of clonality and clone size dynamics in ATL onset by modeling data from high-throughput monitoring of HTLV-1 integration sites using single- and multiple-time-point samples.

Results

From four size categories analyzed, we found that big clones (B; 513–2048 infected cells) and very big clones (VB; >2048 infected cells) had prognostic value. No sample harbored two or more VB clones or three or more B clones. We examined the role of clone size, clone combination, and the number of integration sites in the prognosis of infected individuals. We found a moderate reverse correlation between the total number of clones and the size of the largest clone. We devised a data-driven model that allows intuitive representation of clonal composition.

Conclusions

This integration site-based clonality tree model represents the complexity of clonality and provides a global view of clonality data that facilitates the analysis, interpretation, understanding, and visualization of the behavior of clones on inter- and intra-individual scales. It is fully data-driven, intuitively depicts the clonality patterns of HTLV-1-infected individuals and can assist in early risk assessment of ATL onset by reflecting the prognosis of infected individuals. This model should assist in assimilating information on clonal composition and understanding clonal expansion in HTLV-1-infected individuals.

Electronic supplementary material

The online version of this article (doi:10.1186/s40246-017-0112-8) contains supplementary material, which is available to authorized users.

Clonality of HTLV-1-infected T cells as a risk indicator for development and progression of adult T-cell leukemia

Adult T-cell leukemia (ATL) is an aggressive T-cell malignancy caused by human T-cell leukemia virus type 1 (HTLV-1) that develops along a carcinogenic process involving 5 or more genetic events in infected cells. The lifetime incidence of ATL among HTLV-1-infected individuals is approximately 5%. Although epidemiologic studies have revealed risk factors for ATL, the molecular mechanisms that determine the fates of carriers remain unclear. A better understanding of clonal composition and related longitudinal dynamics would clarify the process of ATL leukemogenesis and provide insights into the mechanisms underlying the proliferation of a malignant clone. Genomic DNA samples and clinical information were obtained from individuals enrolled in the Joint Study for Predisposing Factors for ATL Development, a Japanese prospective cohort study. Forty-seven longitudinal samples from 20 individuals (9 asymptomatic carriers and 11 patients with ATL at enrollment) were subjected to a clonality analysis. A method based on next-generation sequencing was used to characterize clones on the basis of integration sites. Relationships were analyzed among clonal patterns, clone sizes, and clinical status, including ATL onset and progression. Among carriers, those exhibiting an oligoclonal or monoclonal pattern with largely expanded clones subsequently progressed to ATL. All indolent patients who progressed to acute-type ATL exhibited monoclonal expansion. In both situations, the major expanded clone after progression was derived from the largest pre-existing clone. This study has provided the first detailed information regarding the dynamics of HTLV-1-infected T-cell clones and collectively suggests that the clonality of HTLV-1-infected cells could be a useful predictive marker of ATL onset and progression.

Adult T-cell leukemia: molecular basis for clonal expansion and transformation of HTLV-1-infected T cells

Adult T-cell leukemia (ATL) is an aggressive T-cell malignancy caused by human T-cell leukemia virus type 1 (HTLV-1) that develops through a multistep carcinogenesis process involving 5 or more genetic events. We provide a comprehensive overview of recently uncovered information on the molecular basis of leukemogenesis in ATL. Broadly, the landscape of genetic abnormalities in ATL that include alterations highly enriched in genes for T-cell receptor-NF-κB signaling such as PLCG1, PRKCB, and CARD11 and gain-of function mutations in CCR4 and CCR7 Conversely, the epigenetic landscape of ATL can be summarized as polycomb repressive complex 2 hyperactivation with genome-wide H3K27 me3 accumulation as the basis of the unique transcriptome of ATL cells. Expression of H3K27 methyltransferase enhancer of zeste 2 was shown to be induced by HTLV-1 Tax and NF-κB. Furthermore, provirus integration site analysis with high-throughput sequencing enabled the analysis of clonal composition and cell number of each clone in vivo, whereas multicolor flow cytometric analysis with CD7 and cell adhesion molecule 1 enabled the identification of HTLV-1-infected CD4+ T cells in vivo. Sorted immortalized but untransformed cells displayed epigenetic changes closely overlapping those observed in terminally transformed ATL cells, suggesting that epigenetic abnormalities are likely earlier events in leukemogenesis. These new findings broaden the scope of conceptualization of the molecular mechanisms of leukemogenesis, dissecting them into immortalization and clonal progression. These recent findings also open a new direction of drug development for ATL prevention and treatment because epigenetic marks can be reprogrammed. Mechanisms underlying initial immortalization and progressive accumulation of these abnormalities remain to be elucidated.

Multidisciplinary insight into clonal expansion of HTLV-1-infected cells in adult T-cell leukemia via modeling by deterministic finite automata coupled with high-throughput sequencing

Background

Clonal expansion of leukemic cells leads to onset of adult T-cell leukemia (ATL), an aggressive lymphoid malignancy with a very poor prognosis. Infection with human T-cell leukemia virus type-1 (HTLV-1) is the direct cause of ATL onset, and integration of HTLV-1 into the human genome is essential for clonal expansion of leukemic cells. Therefore, monitoring clonal expansion of HTLV-1–infected cells via isolation of integration sites assists in analyzing infected individuals from early infection to the final stage of ATL development. However, because of the complex nature of clonal expansion, the underlying mechanisms have yet to be clarified. Combining computational/mathematical modeling with experimental and clinical data of integration site–based clonality analysis derived from next generation sequencing technologies provides an appropriate strategy to achieve a better understanding of ATL development.

Methods

As a comprehensively interdisciplinary project, this study combined three main aspects: wet laboratory experiments, in silico analysis and empirical modeling.

Results

We analyzed clinical samples from HTLV-1–infected individuals with a broad range of proviral loads using a high-throughput methodology that enables isolation of HTLV-1 integration sites and accurate measurement of the size of infected clones. We categorized clones into four size groups, “very small”, “small”, “big”, and “very big”, based on the patterns of clonal growth and observed clone sizes. We propose an empirical formal model based on deterministic finite state automata (DFA) analysis of real clinical samples to illustrate patterns of clonal expansion.

Conclusions

Through the developed model, we have translated biological data of clonal expansion into the formal language of mathematics and represented the observed clonality data with DFA. Our data suggest that combining experimental data (absolute size of clones) with DFA can describe the clonality status of patients. This kind of modeling provides a basic understanding as well as a unique perspective for clarifying the mechanisms of clonal expansion in ATL.

Electronic supplementary material

The online version of this article (doi:10.1186/s12920-016-0241-2) contains supplementary material, which is available to authorized users.

Serological and Molecular Methods to Study Epidemiological Aspects of Human T-Cell Lymphotropic Virus Type 1 Infection

We estimated that at least 5-10 million individuals are infected with HTLV-1. Importantly, this number is based on the study of nearly 1.5 billion people living in known human T-cell lymphotropic virus type 1 (HTLV-1) endemic areas, for which reliable epidemiological data are available. However, for some highly populated regions including India, the Maghreb, East Africa, and some regions of China, no consistent data are yet available which prevents a more accurate estimation. Thus, the number of HTLV-1 infected people in the world is probably much higher. The prevalence of HTLV-1 prevalence varies depending on age, sex, and economic level in most HTLV-1 endemic areas. HTLV-1 seroprevalence gradually increases with age, especially in women. HTLV-1 has a simian origin and was originally acquired by humans through interspecies transmission from STLV-1 infected monkeys in the Old World. Three main modes of HTLV-1 transmission have been described; (1) from mother-to-child after prolonged breast-feeding lasting more than six months, (2) through sexual intercourse, which mainly, but not exclusively, occurs from male to female and lastly, (3) from contaminated blood products, which contain HTLV-1 infected lymphocytes. In specific areas, such as Central Africa, zoonotic transmission from STLV-1 infected monkeys to humans is still ongoing.The diagnostic methods used to study the epidemiological aspects of HTLV-1 infection mainly consist of serological assays for the detection of antibodies specifically directed against different HTLV-1 antigens. Screening tests are usually based on enzyme-linked immunoabsorbent assay (ELISA), chemiluminescence enzyme-linked immunoassay (CLEIA) or particle agglutination (PA). Confirmatory tests include mostly Western blots (WB)s or innogenetics line immunoassay (INNO-LIA™) and to a lesser extent immunofluorescence assay (IFA). The search for integrated provirus in the DNA from peripheral blood cells can be performed by qualitative and/or quantitative polymerase chain reaction (qPCR). qPCR is widely used in most diagnostic laboratories and quantification of proviral DNA is useful for the diagnosis and follow-up of HTLV-1 associated diseases such as adult T-cell leukemia (ATL) and tropical spastic paraparesis/HTLV-1 associated myelopathy (TSP/HAM). PCR also provides amplicons for further sequence analysis to determine the HTLV-1 genotype present in the infected person. The use of new generation sequencing methodologies to molecularly characterize full and/or partial HTLV-1 genomic regions is increasing. HTLV-1 genotyping generates valuable molecular epidemiological data to better understand the evolutionary history of this virus.

High-Throughput Mapping and Clonal Quantification of Retroviral Integration Sites

We describe here a method to identify the position of retroviral insertion sites and simultaneously to quantify the absolute abundance of each clone, i.e., the number of cells having the provirus inserted at a given place in the host genome. The method is based on random shearing of the host cell DNA, followed by a linker-mediated PCR to amplify the genomic regions flanking the proviruses, and high-throughput sequencing of the amplicons. The quantification of the abundance of each infected clone allowed us to develop two new metrics: i. the oligoclonality index, which quantifies the nonuniformity of the distribution of clone abundance, and ii. an estimator of the total number of clones in the body of the host. These new tools are valuable for the study of retroviral infections and can also be adapted for the tracking of gene-edited cells.

Transition of adult T-cell leukemia/lymphoma clones during clinical progression

Adult T-cell leukemia/lymphoma (ATLL) is a peripheral T-cell neoplasm caused by the transformation of HTLV-1-infected T cells. ATLL, especially its aggressive form, is known for its poor prognosis, even with intensive chemotherapy. ATLL cells are considered to be monoclonal; however, multiclonal proliferation or emergence of a new clone over time has been reported based on Southern blot analysis, although direct molecular evidence remains elusive. Furthermore, it is thought that clonal change may be a cause of early drug resistance in ATLL. To directly analyze potential clonal changes in ATLL during its clinical course, we used inverse PCR to detect integration sites in combination with a newly developed method using next-generation sequencing, and compared ATLL cell clonality at different time points. The results of inverse PCR indicated that the major clone was altered in three of 19 patients. Together with results from five patients, using this new method, we found direct evidence of clonal change occurring during the clinical course or in response to chemotherapy in ATLL. These results also highlight the importance of clonality analysis for understanding the mechanisms of ATLL development and drug resistance.

Transcriptional and Epigenetic Regulatory Mechanisms Affecting HTLV-1 Provirus

Human T-cell leukemia virus type 1 (HTLV-1) is a retrovirus associated with human diseases, such as adult T-cell leukemia (ATL) and HTLV-1-associated myelopathy/Tropic spastic paraparesis (HAM/TSP). As a retrovirus, its life cycle includes a step where HTLV-1 is integrated into the host genomic DNA and forms proviral DNA. In the chronic phase of the infection, HTLV‑1 is known to proliferate as a provirus via the mitotic division of the infected host cells. There are generally tens of thousands of infected clones within an infected individual. They exist not only in peripheral blood, but also in various lymphoid organs. Viral proteins encoded in HTLV-1 genome play a role in the proliferation and survival of the infected cells. As is the case with other chronic viral infections, HTLV-1 gene expression induces the activation of the host immunity against the virus. Thus, the transcription from HTLV-1 provirus needs to be controlled in order to evade the host immune surveillance. There should be a dynamic and complex regulation in vivo, where an equilibrium between viral antigen expression and host immune surveillance is achieved. The mechanisms regulating viral gene expression from the provirus are a key to understanding the persistent/latent infection with HTLV-1 and its pathogenesis. In this article, we would like to review our current understanding on this topic.

Dysregulation of c-Myb Pathway by Aberrant Expression of Proto-oncogene MYB Provides the Basis for Malignancy in Adult T-cell Leukemia/lymphoma Cells

PURPOSE

Adult T-cell leukemia/lymphoma (ATLL) is an aggressive human T-cell malignancy induced by human T-lymphotrophic virus-1 (HTLV-1) infection. The genetic alterations in infected cells that lead to transformation have not been completely elucidated, thus hindering the identification of effective therapeutic targets for ATL. Here, we present the first assessment of MYB proto-oncogene dysregulation in ATL and an exploration of its role in the onset of ATL.

EXPERIMENTAL DESIGN

We investigated the expression patterns of MYB splicing variants in ATL. The molecular characteristics of the c-Myb-9A isoform, which was overexpressed in ATL cells, were examined using chromatin immunoprecipitation and promoter assays. We further examined the biologic impacts of abnormal c-Myb overexpression in ATL using overall c-Myb knockdown with shRNA or c-Myb-9A knockdown with morpholino oligomers.

RESULTS

Both total c-Myb and c-Myb-9A, which exhibited strong transforming activity, were overexpressed in ATL cells in a leukemogenesis- and progression-dependent manner. Knockdown of either total c-Myb or c-Myb-9A induced ATL cell death. c-Myb transactivates nine genes that encode essential regulators of cell proliferation and NF-κB signaling. c-Myb-9A induced significantly stronger transactivation of all tested genes and stronger NF-κB activation compared with wild-type c-Myb.

CONCLUSIONS

Our data demonstrate that c-Myb pathway overactivation caused by unbalanced c-Myb-9A overexpression is associated with disorders in cellular homeostasis and consequently, accelerated transformation, cell proliferation, and malignancy in ATL cells. These data support the notion of the c-Myb pathway as a promising new therapeutic target for ATL. Clin Cancer Res; 22(23); 5915-28. ©2016 AACR.

Oncogenic viruses: Lessons learned using next-generation sequencing technologies

Fifteen percent of cancers are driven by oncogenic human viruses. Four of those viruses, hepatitis B virus, human papillomavirus, Merkel cell polyomavirus, and human T-cell lymphotropic virus, integrate the host genome. Viral oncogenesis is the result of epigenetic and genetic alterations that happen during viral integration. So far, little data have been available regarding integration mechanisms and modifications in the host genome. However, the emergence of high-throughput sequencing and bioinformatic tools enables researchers to establish the landscape of genomic alterations and predict the events that follow viral integration. Cooperative working groups are currently investigating these factors in large data sets. Herein, we provide novel insights into the initiating events of cancer onset during infection with integrative viruses. Although much remains to be discovered, many improvements are expected from the clinical point of view, from better prognosis classifications to better therapeutic strategies.

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.

Advanced human T-cell leukemia virus type 1 carriers and early-stage indolent adult T-cell leukemia-lymphoma are indistinguishable based on CADM1 positivity in flow cytometry

We previously reported that the cell adhesion molecule 1 (CADM1) versus CD7 plot in flow cytometry reflects disease progression in human T-cell leukemia virus type 1 (HTLV-1) infection. In CD4⁺ cells from peripheral blood, CADM1⁻CD7⁺ (P), CADM1⁺CD7^dim (D) and CADM1⁺CD7⁻ (N) subpopulations are observed. The D and N subpopulations increase as asymptomatic HTLV-1 carriers (AC) progress to indolent adult T-cell leukemia-lymphoma (ATL) and the N subpopulation then expands in aggressive ATL. In the present study we examined whether the analysis can estimate the risk of developing ATL in advanced AC. Peripheral blood samples from AC (N = 41) and indolent ATL patients (N = 19) were analyzed by flow cytometry using the CADM1 versus CD7 plot for CD4⁺ cells and inverse long PCR (clonality analysis) of FACS-sorted subpopulations. Almost all AC with a high HTLV-1 proviral load (>4 copies/100 cells) had a CADM1⁺ (D + N) frequency of >10%. AC with 25% < CADM1⁺ ≤ 50% contained expanded clones similar to smoldering-type ATL. In many patients in the 25% < CADM1⁺ ≤ 50% group, the proportion of abnormal lymphocytes was distributed around the 5% line, which divides AC and smoldering-type ATL in Shimoyama's classification. In conclusion, the CADM1 versus CD7 plot is useful for selection of putative high-risk AC. The characteristics of some AC and smoldering ATL are said to be similar; however, long-term follow up is required and the clinical outcome (e.g. rate of transformation) of these cases should be used to determine whether to include them in the same clinical category.

Bioinformatics and HIV latency

Despite effective treatment, HIV is not completely eliminated from the infected organism because of the existence of viral reservoirs. A major reservoir consists of infected resting CD4+ T cells, mostly of memory type, that persist over time due to the stable proviral insertion and a long cellular lifespan. Resting cells do not produce viral particles and are protected from viral-induced cytotoxicity or immune killing. However, these latently infected cells can be reactivated by stochastic events or by external stimuli. The present review focuses on novel genome-wide technologies applied to the study of integration, transcriptome, and proteome characteristics and their recent contribution to the understanding of HIV latency.

Quantitative shearing linear amplification polymerase chain reaction: an improved method for quantifying lentiviral vector insertion sites in transplanted hematopoietic cell systems

In gene therapy trials targeting blood disorders, it is important to detect dominance of transduced hematopoietic stem cell (HSC) clones arising from vector insertion site (VIS) effects. Current methods for VIS analysis often do not have defined levels of quantitative accuracy and therefore can fail to detect early clonal dominance. We have developed a rapid and inexpensive method for measuring clone size based on random shearing of genomic DNA, minimal exponential PCR amplification, and shear site counts as a quantitative endpoint. This quantitative shearing linear amplification PCR (qsLAM PCR) assay utilizes an internal control sample containing 19 lentiviral insertion sites per cell that is mixed with polyclonal samples derived from transduced human CD34+ cells. Samples were analyzed from transplanted pigtail macaques and from a participant in our X-linked severe combined immunodeficiency (XSCID) lentiviral vector trial and yielded controlled and quantitative results in all cases. One case of early clonal dominance was detected in a monkey transplanted with limiting numbers of transduced HSCs, while the clinical samples from the XSCID trial participant showed highly diverse clonal representation. These studies demonstrate that qsLAM PCR is a facile and quantitative assay for measuring clonal repertoires in subjects enrolled in human gene therapy trials using lentiviral-transduced HSCs.

Identification of TL-Om1, an adult T-cell leukemia (ATL) cell line, as reference material for quantitative PCR for human T-lymphotropic virus 1

Quantitative PCR (qPCR) for human T-lymphotropic virus 1 (HTLV-1) is useful for measuring the amount of integrated HTLV-1 proviral DNA in peripheral blood mononuclear cells. Many laboratories in Japan have developed different HTLV-1 qPCR methods. However, when six independent laboratories analyzed the proviral load of the same samples, there was a 5-fold difference in their results. To standardize HTLV-1 qPCR, preparation of a well-defined reference material is needed. We analyzed the integrated HTLV-1 genome and the internal control (IC) genes of TL-Om1, a cell line derived from adult T-cell leukemia, to confirm its suitability as a reference material for HTLV-1 qPCR. Fluorescent in situ hybridization (FISH) showed that HTLV-1 provirus was monoclonally integrated in chromosome 1 at the site of 1p13 in the TL-Om1 genome. HTLV-1 proviral genome was not transferred from TL-Om1 to an uninfected T-cell line, suggesting that the HTLV-1 proviral copy number in TL-Om1 cells is stable. To determine the copy number of HTLV-1 provirus and IC genes in TL-Om1 cells, we used FISH, digital PCR, and qPCR. HTLV-1 copy numbers obtained by these three methods were similar, suggesting that their results were accurate. Also, the ratio of the copy number of HTLV-1 provirus to one of the IC genes, RNase P, was consistent for all three methods. These findings indicate that TL-Om1 cells are an appropriate reference material for HTLV-1 qPCR.