Retrovirus Integration Database (RID): a public database for retroviral insertion sites into host genomes

Divergent populations of HIV-infected naive and memory CD4+ T cell clones in children on antiretroviral therapy

BACKGROUNDNaive cells comprise 90% of the CD4+ T cell population in neonates and exhibit distinct age-specific capacities for proliferation and activation. We hypothesized that HIV-infected naive CD4+ T cell populations in children on long-term antiretroviral therapy (ART) would thus be distinct from infected memory cells.METHODSPeripheral blood naive and memory CD4+ T cells from 8 children with perinatal HIV on ART initiated at age 1.7-17 months were isolated by FACS. DNA was extracted from sorted cells, and HIV proviruses were counted, evaluated for intactness, and subjected to integration site analysis (ISA).RESULTSNaive CD4+ T cells containing HIV proviruses were detected in children with 95% statistical confidence. A median 4.7% of long terminal repeat-containing naive CD4+ T cells also contained HIV genetic elements consistent with intactness. Full-length proviral sequencing confirmed intactness of 1 provirus. In the participant with the greatest degree of naive cell infection, ISA revealed infected expanded cell clones in both naive and memory T cells, with no common HIV integration sites detected between subsets. Divergent integration site profiles reflected differential gene expression patterns of naive and memory T cells.CONCLUSIONThese results demonstrate that HIV persisted in both naive and memory CD4+ T cells that underwent clonal expansion and harbored intact proviruses, and suggest that infected memory T cell clones do not frequently arise from naive cell differentiation in children with perinatal HIV on long-term ART.FUNDINGCenter for Cancer Research, NCI; Office of AIDS Research; NCI FLEX; Children's and Emory Junior Faculty Focused Award.

Longitudinal analysis of lentiviral and retroviral chimeric antigen receptors' integration sites reveals distinct clonal evolutionary patterns

T-cell malignancies following chimeric antigen receptor (CAR) therapies are partly related to insertional mutagenesis, but the longitudinal evolution of CAR-integration sites (IS) remains understudied. We performed an IS analysis in blood from three tisagenlecleucel (lentiviral), one axicabtagene-ciloleucel and one brexucabtagene-autoleucel (gammaretrovirals) patient at peak expansion and 1-year follow-up. All were complete responders. Peak expansion IS patterns were vector dependent: lentiviral CAR integrated mostly in introns and gammaretrovirals in intergenic regions, closer to transcription start sites. At 1-year post-infusion, lentiviral CAR showed no major clonal proliferation. Gammaretroviral CARs had divergent outcomes: no detectable CAR (axicabtagene-ciloleucel) or low-level oligoclonal persistence (brexucabtagene-autoleucel). Whether this latter evolution is at risk of further CAR malignancies needs further investigations.

SIV infection and ARV treatment reshape the transcriptional and epigenetic profile of naïve and memory T cells in vivo

Human and simian immunodeficiency viruses (HIV and SIV) are lentiviruses that reverse transcribe their RNA genome with subsequent integration into the genome of the target cell. How progressive infection and administration of antiretrovirals (ARVs) longitudinally influence the transcriptomic and epigenetic landscape of particular T cell subsets, and how these may influence the genetic location of integration are unclear. Here, we use RNAseq and ATACseq to study the transcriptomics and epigenetic landscape of longitudinally sampled naïve and memory CD4+ and CD8+ T cells in two species of non-human primates prior to SIV infection, during chronic SIV infection, and after administration of ARVs. We find that SIV infection leads to significant alteration to the transcriptomic profile of all T cell subsets that are only partially reversed by administration of ARVs. Epigenetic changes were more apparent in animals with longer periods of untreated SIV infection and correlated well with changes in corresponding gene expression. Known SIV integration sites did not vary due to SIV status but did contain more open chromatin in rhesus macaque memory T cells, and the expression of proteasome-related genes at the pre-SIV timepoint correlated with subsequent viremia.IMPORTANCEChronic inflammation during progressive human and simian immunodeficiency virus (HIV and SIV) infections leads to significant co-morbidities in infected individuals with significant consequences. Antiretroviral (ARV)-treated individuals also manifest increased levels of inflammation which are associated with increased mortalities. These data will help guide rational development of modalities to reduce inflammation observed in people living with HIV and suggest mechanisms underlying lentiviral integration site preferences.

InCliniGene enables high-throughput and comprehensive in vivo clonal tracking toward clinical genomics data integration

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

2-Hydroxypropyl-β-Cyclodextrin Treatment Induces Modest Immune Activation in Healthy Rhesus Macaques

Experimental simian immunodeficiency virus (SIV) infection of Asian macaques is an excellent model for HIV disease progression and therapeutic development. Recent coformulations of nucleoside analogs and an integrase inhibitor have been used for parenteral antiretroviral (ARV) administration in SIV-infected macaques, successfully resulting in undetectable plasma SIV RNA. In a cohort of SIVmac239-infected macaques, we recently observed that administration of coformulated ARVs resulted in an unexpected increase in plasma levels of soluble CD14 (sCD14), associated with stimulation of myeloid cells. We hypothesized that the coformulation solubilizing agent Kleptose (2-hydroxypropyl-β-cyclodextrin [HPβCD]) may induce inflammation with myeloid cell activation and the release of sCD14. Herein, we stimulated peripheral blood mononuclear cells (PBMCs) from healthy macaques with HPβCD from different commercial sources and evaluated inflammatory cytokine production in vitro. Treatment of PBMCs resulted in increased sCD14 release and myeloid cell interleukin-1β (IL-1β) production-with stimulation varying significantly by HPβCD source-and destabilized lymphocyte CCR5 surface expression. We further treated healthy macaques with Kleptose alone. In vivo, we observed modestly increased myeloid cell activation in response to Kleptose treatment without significant perturbation of the immunological transcriptome or epigenome. Our results demonstrate a need for vehicle-only controls and highlight immunological perturbations that can occur when using HPβCD in pharmaceutical coformulations. IMPORTANCE SIV infection of nonhuman primates is the principal model system for assessing HIV disease progression and therapeutic development. HPβCD has recently been incorporated as a solubilizing agent in coformulations of ARVs in SIV-infected nonhuman primates. Although HPβCD has historically been considered inert, recent findings suggest that HPβCD may contribute to inflammation. Herein, we investigate the contribution of HPβCD to healthy macaque inflammation in vitro and in vivo. We observe that HPβCD causes an induction of sCD14 and IL-1β from myeloid cells in vitro and demonstrate that HPβCD stimulatory capacity varies by commercial source. In vivo, we observe modest myeloid cell activation in blood and bronchoalveolar lavage specimens absent systemic immune activation. From our findings, it is unclear whether HPβCD stimulation may improve or diminish immune reconstitution in ARV-treated lentiviral infections. Our results demonstrate a need for vehicle-only controls and highlight immunological perturbations that can occur when using HPβCD in pharmaceutical coformulations.

Deep learning for detecting and elucidating human T-cell leukemia virus type 1 integration in the human genome

Human T-cell leukemia virus type 1 (HTLV-1), a retrovirus, is the causative agent for adult T cell leukemia/lymphoma and many other human diseases. Accurate and high throughput detection of HTLV-1 virus integration sites (VISs) across the host genomes plays a crucial role in the prevention and treatment of HTLV-1-associated diseases. Here, we developed DeepHTLV, the first deep learning framework for VIS prediction de novo from genome sequence, motif discovery, and cis-regulatory factor identification. We demonstrated the high accuracy of DeepHTLV with more efficient and interpretive feature representations. Decoding the informative features captured by DeepHTLV resulted in eight representative clusters with consensus motifs for potential HTLV-1 integration. Furthermore, DeepHTLV revealed interesting cis-regulatory elements in regulation of VISs that have significant association with the detected motifs. Literature evidence demonstrated nearly half (34) of the predicted transcription factors enriched with VISs were involved in HTLV-1-associated diseases. DeepHTLV is freely available at https://github.com/bsml320/DeepHTLV.

Linking Human Betaretrovirus with Autoimmunity and Liver Disease in Patients with Primary Biliary Cholangitis

Primary biliary cholangitis (PBC) is an autoimmune liver disease characterized by the production of diagnostic antimitochondrial antibodies (AMA) reactive to the pyruvate dehydrogenase complex. A human betaretrovirus (HBRV) resembling mouse mammary tumor virus has been characterized in patients with PBC. However, linking the viral infection with the disease is not a straight-forward process because PBC is a complex multifactorial disease influenced by genetic, hormonal, autoimmune, environmental, and other factors. Currently, PBC is assumed to have an autoimmune etiology, but the evidence is lacking to support this conjecture. In this review, we describe different approaches connecting HBRV with PBC. Initially, we used co-cultivation of HBRV with biliary epithelial cells to trigger the PBC-specific phenotype with cell surface expression of cryptic mitochondrial autoantigens linked with antimitochondrial antibody expression. Subsequently, we have derived layers of proof to support the role of betaretrovirus infection in mouse models of autoimmune biliary disease with spontaneous AMA production and in patients with PBC. Using Hill’s criteria, we provide an overview of how betaretrovirus infection may trigger autoimmunity and propagate biliary disease. Ultimately, the demonstration that disease can be cured with antiviral therapy may sway the argument toward an infectious disease etiology in an analogous fashion that was used to link H. pylori with peptic ulcer disease.

Isolation of a Human Betaretrovirus from Patients with Primary Biliary Cholangitis

A human betaretrovirus (HBRV) has been linked with the autoimmune liver disease, primary biliary cholangitis (PBC), and various cancers, including breast cancer and lymphoma. HBRV is closely related to the mouse mammary tumor virus, and represents the only exogenous betaretrovirus characterized in humans to date. Evidence of infection in patients with PBC has been demonstrated through the identification of proviral integration sites in lymphoid tissue, the major reservoir of infection, as well as biliary epithelium, which is the site of the disease process. Accordingly, we tested the hypothesis that patients with PBC harbor a transmissible betaretrovirus by co-cultivation of PBC patients’ lymph node homogenates with the HS578T breast cancer line. Because of the low level of HBRV replication, betaretrovirus producing cells were subcloned to optimize viral isolation and production. Evidence of infection was provided by electron microscopy, RT-PCR, in situ hybridization, cloning of the HBRV proviral genome and demonstration of more than 3400 integration sites. Further evidence of viral transmissibility was demonstrated by infection of biliary epithelial cells. While HBRV did not show a preference for integration proximal to specific genomic features, analyses of common insertion sites revealed evidence of integration proximal to cancer associated genes. These studies demonstrate the isolation of HBRV with features similar to mouse mammary tumor virus and confirm that patients with PBC display evidence of a transmissible viral infection.

Collaborative Mining of Whole Genome Sequences for Intelligent HIV-1 Sub-Strain(s) Discovery

BACKGROUND

Effective global antiretroviral vaccines and therapeutic strategies depend on the diversity, evolution, and epidemiology of their various strains as well as their transmission and pathogenesis. Most viral disease-causing particles are clustered into a taxonomy of subtypes to suggest pointers toward nucleotide-specific vaccines or therapeutic applications of clinical significance sufficient for sequence-specific diagnosis and homologous viral studies. These are very useful to formulate predictors to induce cross-resistance to some retroviral control drugs being used across study areas.

OBJECTIVE

This research proposed a collaborative framework of hybridized (Machine Learning and Natural Language Processing) techniques to discover hidden genome patterns and feature predictors, for HIV-1 genome sequences mining.

METHOD

630 human HIV-1 genome sequences above 8500 bps were excavated from the National Center for Biotechnology Information (NCBI) database (https://www.ncbi.nlm.nih.gov) for 21 countries across different continents, Antarctica exempt. These sequences were transformed and learned using a self-organizing map (SOM). To discriminate emerging/new sub-strain(s), the HIV-1 reference genome was included as part of the input isolates/samples during the training. After training the SOM, component planes defining pattern clusters of the input datasets were generated, for cognitive knowledge mining and subsequent labelling of the datasets. Additional genome features including dinucleotide transmission recurrences, codon recurrences, and mutation recurrences, were finally extracted from the raw genomes to construct output classification targets for supervised learning.

RESULTS

SOM training explains the inherent pattern diversity of HIV-1 genomes as well as inter- and intra-country transmissions in which mobility might play an active role, as corroborated by literature. Nine sub-strains were discovered after disassembling the SOM correlation hunting matrix space attributed to disparate clusters. Cognitive knowledge mining separated similar pattern clusters bounded by a certain degree of correlation range, discovered by the SOM. A Kruskal-Wallis rank-sum test and Wilcoxon rank-sum test showed statistically significant variations in dinucleotide, codon, and mutation patterns.

CONCLUSION

Results of the discovered sub-strains and response clusters visualizations corroborate existing literature, with significant haplotype variations. The proposed framework would assist in the development of decision support systems for easy contact tracing, infectious disease surveillance, and studying the progressive evolution of the reference HIV-1 genome.

Acute Myeloid Leukemia Case after Gene Therapy for Sickle Cell Disease

Gene therapy with LentiGlobin for sickle cell disease (bb1111, lovotibeglogene autotemcel) consists of autologous transplantation of a patient's hematopoietic stem cells transduced with the BB305 lentiviral vector that encodes the β^A-T87Q-globin gene. Acute myeloid leukemia developed in a woman approximately 5.5 years after she had received LentiGlobin for sickle cell disease as part of the initial cohort (Group A) of the HGB-206 study. An analysis of peripheral-blood samples revealed that blast cells contained a BB305 lentiviral vector insertion site. The results of an investigation of causality indicated that the leukemia was unlikely to be related to vector insertion, given the location of the insertion site, the very low transgene expression in blast cells, and the lack of an effect on expression of surrounding genes. Several somatic mutations predisposing to acute myeloid leukemia were present after diagnosis, which suggests that patients with sickle cell disease are at increased risk for hematologic malignant conditions after transplantation, most likely because of a combination of risks associated with underlying sickle cell disease, transplantation procedure, and inadequate disease control after treatment. (Funded by Bluebird Bio.).

ViMIC: a database of human disease-related virus mutations, integration sites and cis-effects

Molecular mechanisms of virus-related diseases involve multiple factors, including viral mutation accumulation and integration of a viral genome into the host DNA. With increasing attention being paid to virus-mediated pathogenesis and the development of many useful technologies to identify virus mutations (VMs) and viral integration sites (VISs), much research on these topics is available in PubMed. However, knowledge of VMs and VISs is widely scattered in numerous published papers which lack standardization, integration and curation. To address these challenges, we built a pilot database of human disease-related Virus Mutations, Integration sites and Cis-effects (ViMIC), which specializes in three features: virus mutation sites, viral integration sites and target genes. In total, the ViMIC provides information on 31 712 VMs entries, 105 624 VISs, 16 310 viral target genes and 1 110 015 virus sequences of eight viruses in 77 human diseases obtained from the public domain. Furthermore, in ViMIC users are allowed to explore the cis-effects of virus-host interactions by surveying 78 histone modifications, binding of 1358 transcription regulators and chromatin accessibility on these VISs. We believe ViMIC will become a valuable resource for the virus research community. The database is available at http://bmtongji.cn/ViMIC/index.php.

New Approaches to Multi-Parametric HIV-1 Genetics Using Multiple Displacement Amplification: Determining the What, How, and Where of the HIV-1 Reservoir

Development of potential HIV-1 curative interventions requires accurate characterization of the proviral reservoir, defined as host-integrated viral DNA genomes that drive rebound of viremia upon halting ART (antiretroviral therapy). Evaluation of such interventions necessitates methods capable of pinpointing the rare, genetically intact, replication-competent proviruses within a background of defective proviruses. This evaluation can be achieved by identifying the distinct integration sites of intact proviruses within host genomes and monitoring the dynamics of these proviruses and host cell lineages over longitudinal sampling. Until recently, molecular genetic approaches at the single proviral level have been generally limited to one of a few metrics, such as proviral genome sequence/intactness, host-proviral integration site, or replication competency. New approaches, taking advantage of MDA (multiple displacement amplification) for WGA (whole genome amplification), have enabled multiparametric proviral characterization at the single-genome level, including proviral genome sequence, host-proviral integration site, and phenotypic characterization of the host cell lineage, such as CD4 memory subset and antigen specificity. In this review, we will examine the workflow of MDA-augmented molecular genetic approaches to study the HIV-1 reservoir, highlighting technical advantages and flexibility. We focus on a collection of recent studies in which investigators have used these approaches to comprehensively characterize intact and defective proviruses from donors on ART, investigate mechanisms of elite control, and define cell lineage identity and antigen specificity of infected CD4+ T cell clones. The highlighted studies exemplify how these approaches and their future iterations will be key in defining the targets and evaluating the impacts of HIV curative interventions.

Antigen-driven clonal selection shapes the persistence of HIV-1-infected CD4+ T cells in vivo

Clonal expansion of infected CD4+ T cells is a major mechanism of HIV-1 persistence and a barrier to achieving a cure. Potential causes are homeostatic proliferation, effects of HIV-1 integration, and interaction with antigens. Here, we show that it is possible to link antigen responsiveness, the full proviral sequence, the integration site, and the T cell receptor β-chain (TCRβ) sequence to examine the role of recurrent antigenic exposure in maintaining the HIV-1 reservoir. We isolated CMV- and Gag-responding CD4+ T cells from 10 treated individuals. Proviral populations in CMV-responding cells were dominated by large clones, including clones harboring replication-competent proviruses. TCRβ repertoires showed high clonality driven by converging adaptive responses. Although some proviruses were in genes linked to HIV-1 persistence (BACH2, STAT5B, MKL1), the proliferation of infected cells under antigenic stimulation occurred regardless of the site of integration. Paired TCRβ and integration site analysis showed that infection could occur early or late in the course of a clone's response to antigen and could generate infected cell populations too large to be explained solely by homeostatic proliferation. Together, these findings implicate antigen-driven clonal selection as a major factor in HIV-1 persistence, a finding that will be a difficult challenge to eradication efforts.

rigrag: high-resolution mapping of genic targeting preferences during HIV-1 integration in vitro and in vivo

HIV-1 integration favors recurrent integration gene (RIG) targets and genic proviruses can confer cell survival in vivo. However, the relationship between initial RIG integrants and how these evolve in patients over time are unknown. To address these shortcomings, we built phenomenological models of random integration in silico, which were used to identify 3718 RIGs as well as 2150 recurrent avoided genes from 1.7 million integration sites across 10 in vitro datasets. Despite RIGs comprising only 13% of human genes, they harbored 70% of genic HIV-1 integrations across in vitro and patient-derived datasets. Although previously reported to associate with super-enhancers, RIGs tracked more strongly with speckle-associated domains. While depletion of the integrase cofactor LEDGF/p75 significantly reduced recurrent HIV-1 integration in vitro, LEDGF/p75 primarily occupied non-speckle-associated regions of chromatin, suggesting a previously unappreciated dynamic aspect of LEDGF/p75 functionality in HIV-1 integration targeting. Finally, we identified only six genes from patient samples-BACH2, STAT5B, MKL1, MKL2, IL2RB and MDC1-that displayed enriched integration targeting frequencies and harbored proviruses that likely contributed to cell survival. Thus, despite the known preference of HIV-1 to target cancer-related genes for integration, we conclude that genic proviruses play a limited role to directly affect cell proliferation in vivo.

Early Emergence and Long-Term Persistence of HIV-Infected T-Cell Clones in Children

HIV-1 integrates its genome into the DNA of host cells. Consequently, HIV-1 genomes are copied with the host cell DNA during cellular division.

Little is known about the emergence and persistence of human immunodeficiency virus (HIV)-infected T-cell clones in perinatally infected children. We analyzed peripheral blood mononuclear cells (PBMCs) for clonal expansion in 11 children who initiated antiretroviral therapy (ART) between 1.8 and 17.4 months of age and with viremia suppressed for 6 to 9 years. We obtained 8,662 HIV type 1 (HIV-1) integration sites from pre-ART samples and 1,861 sites from on-ART samples. Expanded clones of infected cells were detected pre-ART in 10/11 children. In 8 children, infected cell clones detected pre-ART persisted for 6 to 9 years on ART. A comparison of integration sites in the samples obtained on ART with healthy donor PBMCs infected ex vivo showed selection for cells with proviruses integrated in BACH2 and STAT5B. Our analyses indicate that, despite marked differences in T-cell composition and dynamics between children and adults, HIV-infected cell clones are established early in children, persist for up to 9 years on ART, and can be driven by proviral integration in proto-oncogenes.

Integration in oncogenes plays only a minor role in determining the in vivo distribution of HIV integration sites before or during suppressive antiretroviral therapy

HIV persists during antiretroviral therapy (ART) as integrated proviruses in cells descended from a small fraction of the CD4+ T cells infected prior to the initiation of ART. To better understand what controls HIV persistence and the distribution of integration sites (IS), we compared about 15,000 and 54,000 IS from individuals pre-ART and on ART, respectively, with approximately 395,000 IS from PBMC infected in vitro. The distribution of IS in vivo is quite similar to the distribution in PBMC, but modified by selection against proviruses in expressed genes, by selection for proviruses integrated into one of 7 specific genes, and by clonal expansion. Clones in which a provirus integrated in an oncogene contributed to cell survival comprised only a small fraction of the clones persisting in on ART. Mechanisms that do not involve the provirus, or its location in the host genome, are more important in determining which clones expand and persist.

In HIV-infected individuals, a small fraction of the infected cells persist and divide. This reservoir persists during fully suppressive ART and can rekindle the infection if ART is discontinued. Because the number of possible sites of HIV DNA integration is very large, each infected cell, and all of its descendants, can be identified by the site where the provirus is integrated (IS). To understand the selective forces that determine the fates of infected cells in vivo, we compared the distribution of HIV IS in freshly-infected cells to cells from HIV-infected donors sampled both before and during ART. We found that, as previously reported, integration favors highly-expressed genes. However, over time, the fraction of cells with proviruses integrated in highly-expressed genes decreases, implying that they grow less well. There are exceptions to this broad negative selection. When a provirus is integrated in a specific region in one of seven genes, the proviruses affect the expression of the target gene, promoting growth and/or survival of the cell. Although this effect is striking, it is only a minor component of the forces that promote the growth and survival of the population of infected cells during ART.

Transcriptional behavior of the HIV-1 promoter in context of the BACH2 prominent proviral integration gene

Chronic infection with human immunodeficiency virus (HIV)-1 is characterized by accumulation of proviral sequences in the genome of target cells. Integration of viral DNA in patients on long-term antiretroviral therapy selectively persists at preferential loci, suggesting site-specific crosstalk of viral sequences and human genes. This crosstalk likely contributes to chronic HIV disease through modulation of host immune pathways and emergence of clonal infected cell populations. To systematically interrogate such effects, we undertook genome engineering to generate Jurkat cell models that replicate integration of HIV-1 long terminal repeat (LTR) sequences at the BTB and CNC Homolog 2 (BACH2) integration locus. This locus is a prominent HIV-1 integration gene in chronic infection, found in 30 % of long-term treated patients with mapped proviral integrations. Using five clonal models carrying an LTR-driven reporter at different BACH2 intergenic regions, we here show that LTR transcriptional activity is repressed in BACH2 regions associated with proviral-DNA integrations in vivo but not in a control region. Our data indicates that this repression is in part epigenetically regulated, particularly through DNA methylation. Importantly, we demonstrate that transcriptional activity of the LTR is independent of BACH2 gene transcription and vice versa in our models. This suggests no transcriptional interference of endogenous and HIV-1 promoters. Taken together, our study provides first insights into how activity of HIV-1 LTR sequences is regulated at the BACH2 locus as prominent example for a recurrently-detected integration gene in chronic infection. Given the importance of integration-site dependent virus/host crosstalk for chronic HIV disease, our findings for the BACH2 locus have potential implications for future therapeutic strategies.

TADs enriched in histone H1.2 strongly overlap with the B compartment, inaccessible chromatin, and AT-rich Giemsa bands

Giemsa staining of metaphase chromosomes results in a characteristic banding useful for identification of chromosomes and its alterations. We have investigated in silico whether Giemsa bands (G bands) correlate with epigenetic and topological features of the interphase genome. Staining of G-positive bands decreases with GC content; nonetheless, G-negative bands are GC heterogeneous. High GC bands are enriched in active histone marks, RNA polymerase II, and SINEs and associate with gene richness, gene expression, and early replication. Low GC bands are enriched in repressive marks, lamina-associated domains, and LINEs. Histone H1 variants distribute heterogeneously among G bands: H1X is enriched at high GC bands and H1.2 is abundant at low GC, compacted bands. According to epigenetic features and H1 content, G bands can be organized in clusters useful to compartmentalize the genome. Indeed, we have obtained Hi-C chromosome interaction maps and compared topologically associating domains (TADs) and A/B compartments to G banding. TADs with high H1.2/H1X ratio strongly overlap with B compartment, late replicating, and inaccessible chromatin and low GC bands. We propose that GC content is a strong driver of chromatin compaction and 3D genome organization, that Giemsa staining recapitulates this organization denoted by high-throughput techniques, and that H1 variants distribute at distinct chromatin domains. DATABASES: Hi-C data on T47D breast cancer cells have been deposited in NCBI's Gene Expression Omnibus and are accessible through GEO Series accession number GSE147627.

Targeted deep sequencing reveals clonal and subclonal mutational signatures in Adult T-cell leukemia/lymphoma and defines an unfavorable indolent subtype

Adult T-cell leukemia/lymphoma (ATL) carries a poor prognosis even in indolent subtypes. We performed targeted deep sequencing combined with mapping of HTLV-1 proviral integration sites of 61 ATL patients of African and Caribbean origin. This revealed mutations mainly affecting TCR/NF-kB (74%), T-cell trafficking (46%), immune escape (29%), and cell cycle (26%) related pathways, consistent with the genomic landscape previously reported in a large Japanese cohort. To examine the evolution of mutational signatures upon disease progression while tracking the viral integration architecture of the malignant clone, we carried out a longitudinal study of patients who either relapsed or progressed from an indolent to an aggressive subtype. Serial analysis of relapsing patients identified several patterns of clonal evolution. In progressing patients, the longitudinal study revealed NF-kB/NFAT mutations at progression that were present at a subclonal level at diagnosis (allelic frequency < 5%). Moreover, the presence in indolent subtypes of mutations affecting the TCR/NF-kB pathway, whether clonal or subclonal, was associated with significantly shorter time to progression and overall survival. Our observations reveal the clonal dynamics of ATL mutational signatures at relapse and during progression. Our study defines a new subgroup of indolent ATLs characterized by a mutational signature at high risk of transformation.

VISDB: a manually curated database of viral integration sites in the human genome

Virus integration into the human genome occurs frequently and represents a key driving event in human disease. Many studies have reported viral integration sites (VISs) proximal to structural or functional regions of the human genome. Here, we systematically collected and manually curated all VISs reported in the literature and publicly available data resources to construct the Viral Integration Site DataBase (VISDB, https://bioinfo.uth.edu/VISDB). Genomic information including target genes, nearby genes, nearest transcription start site, chromosome fragile sites, CpG islands, viral sequences and target sequences were integrated to annotate VISs. We further curated VIS-involved oncogenes and tumor suppressor genes, virus–host interactions involved in non-coding RNA (ncRNA), target gene and microRNA expression in five cancers, among others. Moreover, we developed tools to visualize single integration events, VIS clusters, DNA elements proximal to VISs and virus–host interactions involved in ncRNA. The current version of VISDB contains a total of 77 632 integration sites of five DNA viruses and four RNA retroviruses. VISDB is currently the only active comprehensive VIS database, which provides broad usability for the study of disease, virus related pathophysiology, virus biology, host–pathogen interactions, sequence motif discovery and pattern recognition, molecular evolution and adaption, among others.

DeepHINT: understanding HIV-1 integration via deep learning with attention

MOTIVATION

Human immunodeficiency virus type 1 (HIV-1) genome integration is closely related to clinical latency and viral rebound. In addition to human DNA sequences that directly interact with the integration machinery, the selection of HIV integration sites has also been shown to depend on the heterogeneous genomic context around a large region, which greatly hinders the prediction and mechanistic studies of HIV integration.

RESULTS

We have developed an attention-based deep learning framework, named DeepHINT, to simultaneously provide accurate prediction of HIV integration sites and mechanistic explanations of the detected sites. Extensive tests on a high-density HIV integration site dataset showed that DeepHINT can outperform conventional modeling strategies by automatically learning the genomic context of HIV integration from primary DNA sequence alone or together with epigenetic information. Systematic analyses on diverse known factors of HIV integration further validated the biological relevance of the prediction results. More importantly, in-depth analyses of the attention values output by DeepHINT revealed intriguing mechanistic implications in the selection of HIV integration sites, including potential roles of several DNA-binding proteins. These results established DeepHINT as an effective and explainable deep learning framework for the prediction and mechanistic study of HIV integration.

AVAILABILITY AND IMPLEMENTATION

DeepHINT is available as an open-source software and can be downloaded from https://github.com/nonnerdling/DeepHINT.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Joint profiling of chromatin accessibility and CAR-T integration site analysis at population and single-cell levels

Chimeric antigen receptor (CAR)-T immunotherapy has yielded impressive results in several B cell malignancies, establishing itself as a powerful means to redirect the natural properties of T lymphocytes. In this strategy, the T cell genome is modified by the integration of lentiviral vectors encoding CAR that direct tumor cell killing. However, this therapeutic approach is often limited by the extent of CAR-T cell expansion in vivo. A major outstanding question is whether or not CAR-T integration itself enhances the proliferative competence of individual T cells by rewiring their regulatory landscape. To address this question, it is critical to define the identity of an individual CAR-T cell and simultaneously chart where the CAR-T vector integrates into the genome. Here, we report the development of a method called EpiVIA (https://github.com/VahediLab/epiVIA) for the joint profiling of the chromatin accessibility and lentiviral integration site analysis at the population and single-cell levels. We validate our technique in clonal cells with previously defined integration sites and further demonstrate the ability to measure lentiviral integration sites and chromatin accessibility of host and viral genomes at the single-cell resolution in CAR-T cells. We anticipate that EpiVIA will enable the single-cell deconstruction of gene regulation during CAR-T therapy, leading to the discovery of cellular factors associated with durable treatment.

An analytical pipeline for identifying and mapping the integration sites of HIV and other retroviruses

BACKGROUND

All retroviruses, including human immunodeficiency virus (HIV), must integrate a DNA copy of their genomes into the genome of the infected host cell to replicate. Although integrated retroviral DNA, known as a provirus, can be found at many sites in the host genome, integration is not random. The adaption of linker-mediated PCR (LM-PCR) protocols for high-throughput integration site mapping, using randomly-sheared genomic DNA and Illumina paired-end sequencing, has dramatically increased the number of mapped integration sites. Analysis of samples from human donors has shown that there is clonal expansion of HIV infected cells and that clonal expansion makes an important contribution to HIV persistence. However, analysis of HIV integration sites in samples taken from patients requires extensive PCR amplification and high-throughput sequencing, which makes the methodology prone to certain specific artifacts.

RESULTS

To address the problems with artifacts, we use a comprehensive approach involving experimental procedures linked to a bioinformatics analysis pipeline. Using this combined approach, we are able to reduce the number of PCR/sequencing artifacts that arise and identify the ones that remain. Our streamlined workflow combines random cleavage of the DNA in the samples, end repair, and linker ligation in a single step. We provide guidance on primer and linker design that reduces some of the common artifacts. We also discuss how to identify and remove some of the common artifacts, including the products of PCR mispriming and PCR recombination, that have appeared in some published studies. Our improved bioinformatics pipeline rapidly parses the sequencing data and identifies bona fide integration sites in clonally expanded cells, producing an Excel-formatted report that can be used for additional data processing.

CONCLUSIONS

We provide a detailed protocol that reduces the prevalence of artifacts that arise in the analysis of retroviral integration site data generated from in vivo samples and a bioinformatics pipeline that is able to remove the artifacts that remain.

Dynamic Shifts in the HIV Proviral Landscape During Long Term Combination Antiretroviral Therapy: Implications for Persistence and Control of HIV Infections

Combination antiretroviral therapy (cART) controls but does not eradicate HIV infection; HIV persistence is the principal obstacle to curing infections. The proportion of defective proviruses increases during cART, but the dynamics of this process are not well understood, and a quantitative analysis of how the proviral landscape is reshaped after cART is initiated is critical to understanding how HIV persists. Here, we studied longitudinal samples from HIV infected individuals undergoing long term cART using multiplexed Droplet Digital PCR (ddPCR) approaches to quantify the proportion of deleted proviruses in lymphocytes. In most individuals undergoing cART, HIV proviruses that contain gag are lost more quickly than those that lack gag. Increases in the fraction of gag-deleted proviruses occurred only after 1–2 years of therapy, suggesting that the immune system, and/or toxicity of viral re-activation helps to gradually shape the proviral landscape. After 10–15 years on therapy, there were as many as 3.5–5 times more proviruses in which gag was deleted or highly defective than those containing intact gag. We developed a provirus-specific ddPCR approach to quantify individual clones. Investigation of a clone of cells containing a deleted HIV provirus integrated in the HORMAD2 gene revealed that the cells underwent a massive expansion shortly after cART was initiated until the clone, which was primarily in effector memory cells, dominated the population of proviruses for over 6 years. The expansion of this HIV-infected clone had substantial effects on the overall proviral population.

The Impact of Cellular Proliferation on the HIV-1 Reservoir

Human immunodeficiency virus (HIV) is a chronic infection that destroys the immune system in infected individuals. Although antiretroviral therapy is effective at preventing infection of new cells, it is not curative. The inability to clear infection is due to the presence of a rare, but long-lasting latent cellular reservoir. These cells harboring silent integrated proviral genomes have the potential to become activated at any moment, making therapy necessary for life. Latently-infected cells can also proliferate and expand the viral reservoir through several methods including homeostatic proliferation and differentiation. The chromosomal location of HIV proviruses within cells influences the survival and proliferative potential of host cells. Proliferating, latently-infected cells can harbor proviruses that are both replication-competent and defective. Replication-competent proviral genomes contribute to viral rebound in an infected individual. The majority of available techniques can only assess the integration site or the proviral genome, but not both, preventing reliable evaluation of HIV reservoirs.

ISDB: a database toolkit for storing and analyzing viral integration site data

SUMMARY

We introduce ISDB, a set of software tools for the creation and administration of relational databases of viral integration site (IS) data. Using ISDB, investigators can curate a private database from any heterogeneous set of data sources, including previously-published datasets and internal, work-in-progress data. To make data visible and accessible to collaborators with varying degrees of computational expertise, ISDB automatically generates web sites describing database contents and data exports in several common formats. Compared to a public depository database, the ability to build local, private databases makes ISDB suitable for use in testing hypotheses and developing analyses in the long pre-publication phase of most research.

AVAILABILITY AND IMPLEMENTATION

Installation and usage documentation for ISDB are provided on our website https://mullinslab.microbiol.washington.edu/isdb/. Source code is available under the open source MIT license from https://github.com/MullinsLab/ISDB.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Combined HIV-1 sequence and integration site analysis informs viral dynamics and allows reconstruction of replicating viral ancestors

Understanding HIV-1 persistence despite antiretroviral therapy (ART) is of paramount importance. Both single-genome sequencing (SGS) and integration site analysis (ISA) provide useful information regarding the structure of persistent HIV DNA populations; however, until recently, there was no way to link integration sites to their cognate proviral sequences. Here, we used multiple-displacement amplification (MDA) of cellular DNA diluted to a proviral endpoint to obtain full-length proviral sequences and their corresponding sites of integration. We applied this method to lymph node and peripheral blood mononuclear cells from 5 ART-treated donors to determine whether groups of identical subgenomic sequences in the 2 compartments are the result of clonal expansion of infected cells or a viral genetic bottleneck. We found that identical proviral sequences can result from both cellular expansion and viral genetic bottlenecks occurring prior to ART initiation and following ART failure. We identified an expanded T cell clone carrying an intact provirus that matched a variant previously detected by viral outgrowth assays and expanded clones with wild-type and drug-resistant defective proviruses. We also found 2 clones from 1 donor that carried identical proviruses except for nonoverlapping deletions, from which we could infer the sequence of the intact parental virus. Thus, MDA-SGS can be used for "viral reconstruction" to better understand intrapatient HIV-1 evolution and to determine the clonality and structure of proviruses within expanded clones, including those with drug-resistant mutations. Importantly, we demonstrate that identical sequences observed by standard SGS are not always sufficient to establish proviral clonality.

Spatially clustered loci with multiple enhancers are frequent targets of HIV-1 integration

HIV-1 recurrently targets active genes and integrates in the proximity of the nuclear pore compartment in CD4⁺ T cells. However, the genomic features of these genes and the relevance of their transcriptional activity for HIV-1 integration have so far remained unclear. Here we show that recurrently targeted genes are proximal to super-enhancer genomic elements and that they cluster in specific spatial compartments of the T cell nucleus. We further show that these gene clusters acquire their location during the activation of T cells. The clustering of these genes along with their transcriptional activity are the major determinants of HIV-1 integration in T cells. Our results provide evidence of the relevance of the spatial compartmentalization of the genome for HIV-1 integration, thus further strengthening the role of nuclear architecture in viral infection.

Enhort: a platform for deep analysis of genomic positions

The rise of high-throughput methods in genomic research greatly expanded our knowledge about the functionality of the genome. At the same time, the amount of available genomic position data increased massively, e.g., through genome-wide profiling of protein binding, virus integration or DNA methylation. However, there is no specialized software to investigate integration site profiles of virus integration or transcription factor binding sites by correlating the sites with the diversity of available genomic annotations. Here we present Enhort, a user-friendly software tool for relating large sets of genomic positions to a variety of annotations. It functions as a statistics based genome browser, not focused on a single locus but analyzing many genomic positions simultaneously. Enhort provides comprehensive yet easy-to-use methods for statistical analysis, visualization, and the adjustment of background models according to experimental conditions and scientific questions. Enhort is publicly available online at enhort.mni.thm.de and published under GNU General Public License.

A majority of HIV persistence during antiretroviral therapy is due to infected cell proliferation

Antiretroviral therapy (ART) suppresses viral replication in people living with HIV. Yet, infected cells persist for decades on ART and viremia returns if ART is stopped. Persistence has been attributed to viral replication in an ART sanctuary and long-lived and/or proliferating latently infected cells. Using ecological methods and existing data, we infer that >99% of infected cells are members of clonal populations after one year of ART. We reconcile our results with observations from the first months of ART, demonstrating mathematically how a fossil record of historic HIV replication permits observed viral evolution even while most new infected cells arise from proliferation. Together, our results imply cellular proliferation generates a majority of infected cells during ART. Therefore, reducing proliferation could decrease the size of the HIV reservoir and help achieve a functional cure.

HIV infected cells persist for decades in patients under ART, but the mechanisms responsible remain unclear. Here, Reeves et al. use modeling approaches adapted from ecology to show that cellular proliferation, rather than viral replication, generates a majority of infected cells during ART.

The role of integration and clonal expansion in HIV infection: live long and prosper

Integration of viral DNA into the host genome is a central event in the replication cycle and the pathogenesis of retroviruses, including HIV. Although most cells infected with HIV are rapidly eliminated in vivo, HIV also infects long-lived cells that persist during combination antiretroviral therapy (cART). Cells with replication competent HIV proviruses form a reservoir that persists despite cART and such reservoirs are at the center of efforts to eradicate or control infection without cART. The mechanisms of persistence of these chronically infected long-lived cells is uncertain, but recent research has demonstrated that the presence of the HIV provirus has enduring effects on infected cells. Cells with integrated proviruses may persist for many years, undergo clonal expansion, and produce replication competent HIV. Even proviruses with defective genomes can produce HIV RNA and may contribute to ongoing HIV pathogenesis. New analyses of HIV infected cells suggest that over time on cART, there is a shift in the composition of the population of HIV infected cells, with the infected cells that persist over prolonged periods having proviruses integrated in genes associated with regulation of cell growth. In several cases, strong evidence indicates the presence of the provirus in specific genes may determine persistence, proliferation, or both. These data have raised the intriguing possibility that after cART is introduced, a selection process enriches for cells with proviruses integrated in genes associated with cell growth regulation. The dynamic nature of populations of cells infected with HIV during cART is not well understood, but is likely to have a profound influence on the composition of the HIV reservoir with critical consequences for HIV eradication and control strategies. As such, integration studies will shed light on understanding viral persistence and inform eradication and control strategies. Here we review the process of HIV integration, the role that integration plays in persistence, clonal expansion of the HIV reservoir, and highlight current challenges and outstanding questions for future research.

Selection for avian leukosis virus integration sites determines the clonal progression of B-cell lymphomas

Avian leukosis virus (ALV) is a simple retrovirus that causes a wide range of tumors in chickens, the most common of which are B-cell lymphomas. The viral genome integrates into the host genome and uses its strong promoter and enhancer sequences to alter the expression of nearby genes, frequently inducing tumors. In this study, we compare the preferences for ALV integration sites in cultured cells and in tumors, by analysis of over 87,000 unique integration sites. In tissue culture we observed integration was relatively random with slight preferences for genes, transcription start sites and CpG islands. We also observed a preference for integrations in or near expressed and spliced genes. The integration pattern in cultured cells changed over the course of selection for oncogenic characteristics in tumors. In comparison to tissue culture, ALV integrations are more highly selected for proximity to transcription start sites in tumors. There is also a significant selection of ALV integrations away from CpG islands in the highly clonally expanded cells in tumors. Additionally, we utilized a high throughput method to quantify the magnitude of clonality in different stages of tumorigenesis. An ALV-induced tumor carries between 700 and 3000 unique integrations, with an average of 2.3 to 4 copies of proviral DNA per infected cell. We observed increasing tumor clonality during progression of B-cell lymphomas and identified gene players (especially TERT and MYB) and biological processes involved in tumor progression.

The Avian Leukosis Virus (ALV) is a simple retrovirus that causes cancer in chickens. The virus integrates its genome into the host genome and induces changes in expression of nearby genes. Here, we determine the sites of viral integrations and their role in the progression of tumors. We report pathways and novel gene players that might cooperate and play a role in the progression of B-cell lymphomas. Our study provides new insights into the changes during lymphoma initiation, progression, and metastasis, as a result of selection for specific ALV integration sites.

Cis-perturbation of cancer drivers by the HTLV-1/BLV proviruses is an early determinant of leukemogenesis

Human T-cell leukaemia virus type-1 (HTLV-1) and bovine leukaemia virus (BLV) infect T- and B-lymphocytes, respectively, provoking a polyclonal expansion that will evolve into an aggressive monoclonal leukaemia in ∼5% of individuals following a protracted latency period. It is generally assumed that early oncogenic changes are largely dependent on virus-encoded products, especially TAX and HBZ, while progression to acute leukaemia/lymphoma involves somatic mutations, yet that both are independent of proviral integration site that has been found to be very variable between tumours. Here, we show that HTLV-1/BLV proviruses are integrated near cancer drivers which they affect either by provirus-dependent transcription termination or as a result of viral antisense RNA-dependent cis-perturbation. The same pattern is observed at polyclonal non-malignant stages, indicating that provirus-dependent host gene perturbation contributes to the initial selection of the multiple clones characterizing the asymptomatic stage, requiring additional alterations in the clone that will evolve into full-blown leukaemia/lymphoma.

Human T-cell leukaemia virus type-1 and bovine leukaemia virus infect T and B lymphocytes and lead to aggressive leukaemia. Here, the authors show these proviruses integrate near cancer drivers perturbing transcription termination or antisense RNA-dependent interaction, suggesting post-transcriptional mechanisms in some cases.

HIV-1 Integrates Widely throughout the Genome of the Human Blood Fluke Schistosoma mansoni

Schistosomiasis is the most important helminthic disease of humanity in terms of morbidity and mortality. Facile manipulation of schistosomes using lentiviruses would enable advances in functional genomics in these and related neglected tropical diseases pathogens including tapeworms, and including their non-dividing cells. Such approaches have hitherto been unavailable. Blood stream forms of the human blood fluke, Schistosoma mansoni, the causative agent of the hepatointestinal schistosomiasis, were infected with the human HIV-1 isolate NL4-3 pseudotyped with vesicular stomatitis virus glycoprotein. The appearance of strong stop and positive strand cDNAs indicated that virions fused to schistosome cells, the nucleocapsid internalized and the RNA genome reverse transcribed. Anchored PCR analysis, sequencing HIV-1-specific anchored Illumina libraries and Whole Genome Sequencing (WGS) of schistosomes confirmed chromosomal integration; >8,000 integrations were mapped, distributed throughout the eight pairs of chromosomes including the sex chromosomes. The rate of integrations in the genome exceeded five per 1,000 kb and HIV-1 integrated into protein-encoding loci and elsewhere with integration bias dissimilar to that of human T cells. We estimated ~ 2,100 integrations per schistosomulum based on WGS, i.e. about two or three events per cell, comparable to integration rates in human cells. Accomplishment in schistosomes of post-entry processes essential for HIV-1replication, including integrase-catalyzed integration, was remarkable given the phylogenetic distance between schistosomes and primates, the natural hosts of the genus Lentivirus. These enigmatic findings revealed that HIV-1 was active within cells of S. mansoni, and provided the first demonstration that HIV-1 can integrate into the genome of an invertebrate.