Identification and characterization of a novel enhancer in the HTLV-1 proviral genome

Matsuo, Misaki; Ueno, Takaharu; Monde, Kazuaki; Sugata, Kenji; Tan, Benjy Jek Yang; Rahman, Akhinur; Miyazato, Paola; Uchiyama, Kyosuke; Islam, Saiful; Katsuya, Hiroo; Nakajima, Shinsuke; Tokunaga, Masahito; Nosaka, Kisato; Hata, Hiroyuki; Utsunomiya, Atae; Fujisawa, Jun-ichi; Satou, Yorifumi

doi:10.1038/s41467-022-30029-9

Download PDF

Article
Open access
Published: 03 May 2022

Identification and characterization of a novel enhancer in the HTLV-1 proviral genome

Nature Communications volume 13, Article number: 2405 (2022) Cite this article

4826 Accesses
14 Citations
12 Altmetric
Metrics details

Subjects

Abstract

Human T-cell leukemia virus type 1 (HTLV-1) is a retrovirus that causes adult T-cell leukemia/lymphoma (ATL), a cancer of infected CD4⁺ T-cells. There is both sense and antisense transcription from the integrated provirus. Sense transcription tends to be suppressed, but antisense transcription is constitutively active. Various efforts have been made to elucidate the regulatory mechanism of HTLV-1 provirus for several decades; however, it remains unknown how HTLV-1 antisense transcription is maintained. Here, using proviral DNA-capture sequencing, we found a previously unidentified viral enhancer in the middle of the HTLV-1 provirus. The transcription factors, SRF and ELK-1, play a pivotal role in the activity of this enhancer. Aberrant transcription of genes in the proximity of integration sites was observed in freshly isolated ATL cells. This finding resolves certain long-standing questions concerning HTLV-1 persistence and pathogenesis. We anticipate that the DNA-capture-seq approach can be applied to analyze the regulatory mechanisms of other oncogenic viruses integrated into the host cellular genome.

Engineered CRISPR-Cas12a for higher-order combinatorial chromatin perturbations

Article Open access 17 May 2024

Bidirectional epigenetic editing reveals hierarchies in gene regulation

Article Open access 17 May 2024

Genome-wide Cas9-mediated screening of essential non-coding regulatory elements via libraries of paired single-guide RNAs

Article Open access 22 May 2024

Introduction

Human T-cell leukemia virus type 1 (HTLV-1) is an exogenous retrovirus endemic to some tropical regions^1,2,3. HTLV-1 reverse transcribes its viral RNA genome into double-stranded DNA which is then integrated into the host genomic DNA, forming a provirus, which serves as a template for generating new viral particles. A characteristic of HTLV-1 infection is that the virus maintains its copy number during chronic infection not via the production of free viral particles but via clonal expansion and persistence of infected T-cell clones^4,5. This is the reason the majority of infected individuals remain asymptomatic throughout their lifetime. However, up to 10% of those infected eventually develop adult T-cell leukemia/lymphoma (ATL). ATL can be categorized into two major subtypes: (I) aggressive-type such as acute ATL and lymphoma-type ATL which progress rapidly; and (II) indolent-type ATL such as smoldering ATL and chronic ATL which have a slow disease progression. The pathogenic mechanisms of ATL still remain elusive but is believed to be fueled by the viral genes tax and HBZ which play key roles in the persistence and expansion of infected cells. Tax, a viral protein encoded in the plus strand of HTLV-1, possesses oncogenic functions such as anti-apoptosis and cell proliferation^6,7 while HBZ, which is encoded in the minus strand and transcribed from the 3′LTR, plays a pivotal role in viral persistence and pathogenesis^4,8. Studies have shown that there is a constitutively active antisense transcription from the 3′LTR at the population level while sense transcription from the 5′LTR is frequently silenced in ATL cells in vivo^9,10. This phenomenon suggests that this proviral expression pattern is beneficial for the virus to persist in the host and predisposes infected cells to malignant cellular transformation⁵.

We previously reported the presence of an insulator region in the HTLV-1 provirus¹¹. While this insulator may explain the contrasting transcriptional pattern between the 5′LTR and 3′LTR of the HTLV-1 provirus, it cannot explain the large difference in their transcriptional activity. In general, exogenous and endogenous retroviruses as mobile DNA elements in the genome can be dangerous to host cells because they may act as genome mutagens and induce genomic instability. Therefore, the host cell has evolved defense mechanisms for transcriptional- and posttranscriptional-silencing of such mobile elements. For example, the KRAB ZnF-Trim28-Setdb1-ZFP809 complex induces transcriptional silencing of murine leukemia virus in embryonic stem cells^12,13. Constitutive activation of antisense transcription from integrated HTLV-1 provirus may be an exceptional case, indicating the possibility that there may exist a regulatory mechanism that actively maintains transcription from the 3′LTR.

In this study, we screened transcriptional regulatory regions within the HTLV-1 provirus to identify nucleosome-free regions (NFRs) using a highly sensitive micrococcal nuclease sequencing (MNase-seq) approach, followed by our recently developed HTLV-1 DNA-capture-seq protocol^14,15. The results revealed an internal viral enhancer in the HTLV-1 provirus. It is noteworthy that the region has been intensively analyzed as a coding region of the oncogenic viral gene tax without knowing the enhancer function.

Results

MNase-seq with HTLV-1 DNA-capture-seq identified a significant nucleosome-free region in the HTLV-1 provirus

Transcriptional regulatory regions in the genome, such as promoters and enhancers, are generally nucleosome-free because they need to be accessed by transcription factors, epigenetic modifiers, or chromatin remodelers to exert their regulatory functions. We utilized our recently developed HTLV-1 DNA-capture-seq approach, which allows us to increase the detection sensitivity for HTLV-1 sequences up to several thousandfold^14,15. We first analyzed two HTLV-1-infected T-cell lines, ED and TBX-4B. ED is an ATL cell line derived from an ATL patient, in which sense transcription of the provirus is silenced by DNA methylation and nonsense mutation of the tax gene, while the antisense transcription remains active^16,17 (Fig. 1a). On the other hand, TBX-4B, also a T-cell clone derived from an ATL patient; is a non-malignant clone (not derived from ATL cells)¹⁸. Possibly due to ex vivo cultivation, TBX-4B cells exhibit a higher abundance of HTLV-1 provirus sense transcription compared to antisense strand (Fig. 1a)¹⁹. To identify previously uncharacterized transcriptional regulatory regions, we screened for NFRs in the HTLV-1 provirus in an unbiased manner by performing MNase-seq analysis, where MNase preferentially digests genomic DNA lacking nucleosomes (Fig. 1b). MNase-seq demonstrated a sharp NFR signal between position 7000–7200 of HTLV-1 in ED cells, close to the insulator region we recently reported¹¹ (Fig. 1c). Because insulator regions are known to have a regulatory function, they generally possess open (nucleosome-free) chromatin, and thus can be frequently identified using MNase-seq. In ED cells, we observed that the most nucleosome-depleted region falls between the insulator region and the 3′LTR (Fig. 1c). This region is part of exon 3 of the tax gene (pX region); however, there have been no previous reports regarding its possible function as a regulatory DNA element. We further asked if the NFR is also observed in in vivo samples in addition to the in vitro cell lines by analyzing peripheral mononuclear cells (PBMCs) freshly isolated from ATL patients and an asymptomatic carrier (AC) (Supplementary Table. 1). Consistent with different DNA methylation status of HTLV-1 provirus between cell lines and in naturally virus-infected cells obtained from PBMCs of infected individuals¹⁷, cells from HTLV-1-infected individuals exhibited higher degrees of nucleosome freeness throughout the provirus in comparison to ED cells. Nevertheless, we discovered that the NFR was also present in the same region as in the HTLV-1 infected cell lines (Fig. 1d), indicating that the NFR is present in vivo in virus-infected cells in individuals as well as in vitro cell lines.

**Fig. 1: Transcriptome pattern and nucleosome positions in the HTLV-1 provirus.**

The nucleosome-free region harbors enhancer-related histone modifications and produces enhancer RNAs

To investigate the functional role of the NFR, we performed promoter assays using the promoter of the HBZ gene²⁰ (Fig. 2a). Promoter activity was enhanced by the insertion of the NFR either upstream or downstream of the promoter regardless of orientation, indicating that the NFR has an enhancer function (Fig. 2b). We also evaluated the effect of the NFR on the 5′LTR, which is the promoter for sense transcription in the HTLV-1 provirus (Fig. 2c). The promoter activity of the 5’LTR was also enhanced but at a much smaller effect than that observed for the 3′LTR (Fig. 2c). We observed a similar result when we used the U3 region of 5′LTR as a minimal promoter of sense transcription (Supplementary Fig. 1a) which indicated that the NFR exerts its function mainly on the 3′LTR. Since HTLV-1 Tax is known as a transactivator of 5′LTR, we further analyzed the effect of the NFR on the 5′LTR in the presence of Tax and found that Tax increased the enhancer effect just slightly but significantly (Supplementary Fig. 1b). We next analyzed the promoter/enhancer activity with or without the viral CTCF insulator¹¹. The presence of the insulator region increased the promoter/enhancer activity (Fig. 2c). We then stimulated the T-cells with TNF-α or Phorbol 12- Myristate 13-Acetate (PMA) but no enhancement of promoter activity (Fig. 2d) is observed, indicating that the NFR enhancer activity is not related to T-cell activation status.

We next analyzed enhancer-related histone modifications within the HTLV-1 proviral region. Chromatin immunoprecipitation sequencing (ChIP-seq) signals of enhancer-related histone modifications²¹, including H3K27Ac, H3K4me1, and H3K4me2, showed peaks within a 3-kb distance from the NFR in ED cells (Fig. 2e, f), whereas the peaks were not observed in an infected clone, TBX-4B, without HTLV-1 integration in this locus (Supplementary Fig. 2a). Consistent with the high level of transcriptional activity from the 5′LTR in TBX-4B cells, in which both the 5′ and 3′LTRs are transcriptionally active (Fig. 1a), there was a wide distribution of enhancer-related histone modifications in this clone, not only in the provirus but also in the host genome nearby the viral integration site (IS) (Supplementary Fig. 2b, c). It has been reported that enhancer regions express enhancer RNAs (eRNAs)—non-coding RNAs with divergent orientation from the center of the enhancer²². Thus, we performed a native elongating transcript-cap analysis of gene expression (NET-CAGE) to detect eRNAs²³. NET-CAGE identifies the sequence of the 5′ region of mRNAs or non-coding RNA adjacent to the cap structure using nascent RNA, which is useful in identifying transcriptional start sites and eRNAs at high resolution. eRNAs from the intragenic HTLV-1 enhancer region were detected in ED cells (Fig. 2g) and TBX-4B cells (Supplementary Fig. 2d). These findings demonstrated that the NFR in the HTLV-1 pX region harbors several fundamental features of an enhancer region.

The host transcription factors SRF and ELK-1 bind to the intragenic HTLV-1 enhancer

The NFR region we identified in this study is 160 bp in length. We performed transcription factor binding prediction with the NFR sequence based on the consensus binding motif of various transcription factors and found several candidates (Fig. 3a). We performed a ChIP assay to examine whether these transcription factors localize to the NFR region. Initial evaluation was performed by ChIP-qPCR, and the results were confirmed by ChIP-seq. We found that ELK-1 (Ets like-l protein) and SRF (serum response factor) clearly bound to the NFR but other transcription factors did not (Fig. 3b and data not shown). Since SRF is also involved in the regulation of the 5′LTR²⁴, we also observed the SRF signal in the 5′LTR region in TBX-4B cells, in which HTLV-1 expression in sense orientation is active. Most importantly, the binding of SRF and ELK-1 to the NFR was observed in PBMCs freshly isolated from HTLV-1-infected individuals using highly sensitive ChIP-seq analysis with an HTLV-1 DNA-capture approach¹⁴ (Fig. 3b). These results indicated that this molecular mechanism is actually ongoing in vivo in infected individuals.

**Fig. 3: SRF and ELK-1 bind to the NFR in a DNA sequence-dependent manner.**

Next, we performed electrophoretic mobility shift assays (EMSA) to investigate whether SRF and ELK-1 binding to the NFR depends on DNA sequence. We generated oligonucleotide probes for the NFR using the wild-type (WT) sequence (NFR-wt) and negative control (NC) probes targeting viral regions other than the NFR (Fig. 3c). We observed a band shift when combining the NFR-wt probe and nuclear extract of 293 T cells transfected with SRF and ELK-1 expression vectors (Fig. 3c). Addition of either anti-SRF or anti-ELK-1 antibodies induced a band supershift, demonstrating the involvement of SRF and ELK-1 in the detected band (Fig. 3c). We further generated oligonucleotide probes with mutations in the SRF and/or ELK-1 consensus binding sequence. Mutant_1 (mut_1), mutant_2 (mut_2), and mutant_3 (mut_3) contain mutations in the SRF, ELK-1, or both SRF and ELK-1 binding sites, respectively (Fig. 3d and Supplementary Fig. 3a). To investigate whether the mutations alter transcription factor binding to the NFR, we performed competitive EMSA and found a marked reduction in the binding activity of mutant probes to SRF and ELK-1 compared with that of the wt probe (Fig. 3e). The result showed that there was clear competition observed with unlabeled wt probes but the competition activity was remarkably decreased when we used mut_1, mut_2, and mut_3 probes. As all mutants showed a marked decrease in the formation of the SRF/ELK-1 ternary complex on the NFR DNA, we used mut_3 for subsequent experiments and found a remarkable reduction in the enhancer activity of the NFR after introducing the mutation (Fig. 3f). These results demonstrated that SRF and ELK-1 binding to the NFR plays an indispensable role in enhancer activity.

The SRF and ELK-1 play a critical role in the HTLV-1 enhancer function

Next, we investigated the functional role of SRF/ELK-1 binding to the NFR in the context of the whole viral sequence. As the NFR is located in the coding region of the tax gene, we generated mutations in the SRF/ELK-1 binding site without altering the amino acid sequence of the Tax protein. To investigate the possibility that the nucleotide substitutions might alter the tax mRNA stability and/or translational efficiency resulting in changes in Tax protein levels, we performed western blotting and confirmed that Tax expression with mut_1, mut_2, or mut_3 sequence was equivalent to that with wt sequence (Supplementary Fig. 3b). We constructed HTLV-1 mutant molecular clones containing the same mutations as mut_3 (Fig. 3d) and then transfected HTLV-1-wt or mut_3 plasmids into 293 T cells. We collected cells one day after transfection and quantified viral gene expression in transfected cells and viral production in the culture supernatant (Fig. 4a). We observed a marginal decrease in p19 production in the supernatant of mut_3 plasmid-transfected cells; however, there was no statistically significant difference (Fig. 4b). Nevertheless, there was a significant reduction in tax and HBZ mRNA expression (Fig. 4c). Next, we generated T-cell lines infected with HTLV-1-wt or mut_3 by co-culturing the transfected 293 T cells with JET cells—Jurkat T cells stably transfected with a reporter plasmid to monitor Tax expression—as host cells (Fig. 4d). We sorted Tax-expressing cells 3 days after infection and then analyzed provirus sequences, proviral load, and the distribution of HTLV-1 ISs in the sorted cell populations in bulk. We performed DNA sequencing of the whole integrated provirus by DNA-capture-seq and confirmed that the proviral sequences of JET cells infected with HTLV-1-wt and mut_3 were the same as the plasmid sequences used for transfection (Fig. 4e). The proviral load of HTLV-1-mut_3-infected JET cells were lower than that of HTLV-1-wt-infected ones (Fig. 4f), although that may be at least partially induced by the sorting step. Since we sorted infected cells using a reporter fluorescent protein driven by Tax, the cell sorting efficiency in HTLV-1-mut_3 might be lower than that in HTLV-1-wt due to lower tax expression (Fig. 4c). We next analyzed whether mutations in the SRF/ELK-1 binding site actually reduced SRF/ELK-1 binding to the NFR in the infected cells in vivo. We performed ChIP-seq analysis for SRF and ELK-1 and observed SRF/ELK-1 binding in HTLV-1-wt-infected JET cells but not in mutant virus-infected JET cells (Fig. 4g). Viral IS analysis demonstrated that there were hundreds of different ISs in each condition, i.e., JET cells infected with HTLV-1-wt or mut_3 (Fig. 4h). The distribution of viral IS was not so different between the WT and mutant HTLV-1-infected JET cells in terms of the relationship with the host gene and epigenetic environment (Supplementary Fig. 4a, b and Supplementary Data 1). We then evaluated expression levels of tax and HBZ in JET cells infected with HTLV-1-wt or mut_3 and found that infected cells with HTLV-1-mut_3 showed a significant reduction in tax and HBZ expression (P < 0.05; Fig. 4i). Taking into consideration the similar distribution of ISs between HTLV-1-wt and mut_3-infected cells, the difference in proviral expression can be attributed to the mutation introduced in the NFR of the HTLV-1 provirus and not due to the different distribution of HTLV-1 ISs. HBZ was previously reported to confer an anti-apoptotic phenotype to Jurkat T cells²⁵; therefore, we analyzed the susceptibility to apoptosis induced by T-cell activation and found that JET cells infected with HTLV-1-mut_3 were more susceptible to activation-induced T-cell death than those infected with HTLV-1-wt (Fig. 4j). We further analyzed the effect of mutations in the NFR on chromatin status and found that the mutations induced a decrease in the chromatin openness of the NFR (Fig. 4k). These findings demonstrate that SRF and ELK-1 binding to the enhancer plays a critical role in the enhancer function.

**Fig. 4: Generation and characterization of the HTLV-1 infectious clones with NFR-wt or mut_3 in the SRF/ELK-1 binding sites.**

The intragenic viral enhancer induces upregulation of host genome transcription near the viral IS

The presence of an intragenic viral enhancer in the HTLV-1 provirus raises the possibility that it acts as an ectopic enhancer to activate transcription of host cellular genomic DNA, resulting in changes in host gene expression near the viral IS. To investigate the effect of HTLV-1 integration on host gene expression near the ISs, we cloned JET cells infected with HTLV-1-wt or mut_3 by limiting dilution from bulk cell populations (Fig. 4d) and established five clones infected with HTLV-1-wt with one to four proviruses per clone (Fig. 5a). We also established four clones infected with HTLV-1-mut_3 containing one or two proviruses per clone. The characteristics of each individual clone are shown in Supplementary Table 2. We then performed RNA-seq analysis using these clones and found read-through transcripts around the IS of the JET HTLV-1-wt-infected clone (Fig. 5b) but not in the mutant infected clones (Fig. 5c). We further investigated whether the insertion of an ectopic enhancer by HTLV-1 would alter host gene expression. First, we performed principal component analysis (PCA) to investigate global transcriptional differences among JET clones used in the analysis. The result suggested that there was much less global transcriptional difference among each JET clone when compared with TBX-4B and ED cells (Fig. 5d). We analyzed the expression of host genes within the proximity of HTLV-1 provirus, which we defined as genes found within 100 kb up/downstream from the viral IS, and found that the proportion of upregulated genes in JET clones infected with HTLV-1-wt was significantly higher than those in mutant HTLV-1 clones (P < 0.01; Fig. 5e). It has been reported that the viral CTCF plays a role in chromatin looping with the host CTCF-binding site and induces changes in host gene transcription²⁶. To investigate whether CTCF plays a role in the transactivation of the host gene in cis by HTLV-1 integration, we analyzed CTCF binding to the host genes near ISs. We found a high frequency of CTCF-binding sites in the upregulated host genes, indicating that the transactivation of the host gene by HTLV-1 provirus could be partially mediated by the combination of CTCF and the HTLV-1 enhancer (Fig. 5e). We investigated the possibility that the upregulation of genes near the viral IS was induced via viral gene expression in trans. First, we analyzed the expression level of tax and HBZ genes in each JET clone (Supplementary Fig. 5a). Second, we analyzed publicly available RNA-seq data regarding Jurkat cells with an inducible tax or HBZ gene to investigate whether the genes we show in Fig. 5e would be induced by tax or HBZ expression (Supplementary Fig. 5b)²⁷. These data indicated that upregulation of the host genes near viral IS was not via viral gene expression but at least partially mediated by an ectopic presence of the enhancer inserted by HTLV-1. Further experiments are needed to understand the effect of HTLV-1 integration on the host genome.

**Fig. 5: Establishment and characterization of JET cell clones infected with HTLV-1-wt and mut_3.**

The HTLV-1 ISs are different among JET clones infected with HTLV-1-wt or mut_3, therefore, we cannot exclude the possibility that different IS may generate the different transcriptome in the clones we analyzed in Fig. 5. To solve this issue, we introduced mutations that abrogate the SRF/ELK-1 binding to the enhancer region (Figs. 3d, e and 4e) in a clone infected with HTLV-1-wt (wt_#5, Supplementary Table 2) by using CRISPR/Cas9 technique. SRF/ELK-1 ChIP-seq peaks in HTLV-1-wt-infected cells were abolished in the CRISPR-mutated cells (C-mut), thereby reducing proviral transcription in both sense and antisense strands (Fig. 6a–c) as well as chromatin openness in the enhancer region (Fig. 6d). We further analyzed another clone with two copies of HTLV-1-wt provirus (wt_#1, Supplementary Table 2). We searched the CRISPR-mutated clones with the mut_3 sequence in the enhancer region of two proviruses by doing the whole proviral sequence by HTLV-1 DNA-capture-seq and identified the clone with mut_3 sequence in the NFR in two proviruses while the other regions were identical to the wt_#1 cells. We performed RNA-seq analysis and found that proviral transcription was remarkably decreased (Fig. 6e–g). Importantly, there were a clear decrease in the host transcriptome and splice junction near the viral integration site in C-mut clone (Fig. 6h, i). These data collectively provided the evidence to support the idea that SRF- and ELK-1-binding to the NFR play an important role in the enhancer function.

**Fig. 6: Establishment and characterization of HTLV-1-mut_3 clones by CRISPR/Cas9.**

SRF and ELK-1 localization to the enhancer and aberrant host genome transcription near the proviral integration site in fresh PBMCs

We further investigated the effect of HTLV-1 integration on viral and host genomes by performing mRNA-seq analysis using freshly isolated PBMCs from five ATL cases. All five cases had a high proviral load (Supplementary Table 1) and had a clonally expanded ATL clone (Supplementary Fig. 6). Consistent with the previous reports^9,28, proviral expression in the sense orientation was lower than that in the antisense orientation (Fig. 7a). There was read-through proviral transcription in the sample with HTLV-1 ISs in the host genomic region but not in other samples (Fig. 7b, c), as previously reported²⁹. Interestingly, an ATL case with a defective provirus lacking the 5′LTR also exhibited read-through transcription from the virus to the flanking host genome (Fig. 7d). More importantly, there were clear peaks of SRF and ELK-1 ChIP-seq signals in the integrated proviruses, indicating SRF and ELK-1 play a role in the transcriptional regulation (Fig. 7b–d). PBMCs contain not only ATL cells but also non-ATL infected T cells, uninfected T cells, and various non-T cells; thus, the mRNA-seq data shown in Fig. 7b–d represents the average expression of PBMC subsets. To observe the effect of HTLV-1 ISs on the host genome with high accuracy at single-cell resolution, we performed single-cell RNA-seq analysis using PBMCs from five ATL cases including the same ATL case as in Fig. 7b, c, and in other three ATL cases containing defective proviruses. Based on the T-cell receptor (TCR) clonotype and transcriptome data, we performed clustering analysis and found that the ATL clones, which were identified by the T-cell receptor (TCR) clonotype, clustered differently from the other CD4⁺ T-cell clones (Fig. 7e, f). We then compared the transcriptome near viral IS of CD4⁺ T cells among five ATL cases. There was remarkable upregulation of the local transcriptome only in the sample with viral integration (Fig. 7g, h and Supplementary Fig. 7a–c, left panels), which is consistent with previous reports showing read-through transcript from defective provirus^28,29. Furthermore, there was a significant increase in the local transcriptome in the ATL clone but not in non-ATL CD4⁺ T-cell clones (Fig. 7g, h and Supplementary Fig. 7a–c, right panels). These data were consistent with the idea that the intragenic viral enhancer we identified in this study plays a role in the persistent proviral expression and aberrant transcription of the integrated host genome by recruiting SRF and ELK-1.

**Fig. 7: Transcriptional characterization of the provirus and the flanking host genomes in freshly isolated PBMCs from infected individuals.**

Discussion

The HTLV-1 genome is just over 9000 bp in size but by alternative splicing, it encodes several viral genes which play a role to help the virus achieve persistent infection in the host. Additionally, the provirus is transcribed from both the 3′LTR and the 5′LTR^20,30,31. It has been reported that antisense transcription is frequently expressed in vivo, whereas sense transcription is typically silenced or expressed only intermittently^9,17,32,33. It has not been understood how HTLV-1 antisense transcription remains selectively active. In the present study, we demonstrated the presence of a previously uncharacterized viral enhancer in the HTLV-1 pX region by exploiting the high efficiency and resolution of the viral DNA-capture-seq approach. The enhancer we identified here is located at the 3′ side of the insulator region in the provirus (Fig. 2a). This enhancer is not a typical retroviral enhancer as retroviral enhancers are generally located in the LTR region^30,34. Thus, we propose that this internal enhancer region near the 3′LTR may have two distinct functions: first, to drive the frequent antisense transcription from the 3′LTR, and second, to co-operate with the viral insulator to inhibit the spread of heterochromatin from the 5′LTR towards the 3′LTR. The antisense transcript HBZ plays an indispensable role in viral persistence^9,35 and therefore the intragenic viral enhancer would also contribute to viral persistence via HBZ upregulation. Consistent with this notion, the intragenic viral enhancer and insulator are maintained even in defective type proviruses that are observed in 20–30% of ATL cells^15,36,37. Further experiments are required to understand how these viral regulatory elements, including 5′LTR, viral insulator, enhancer, and 3′LTR, in the small viral genome cooperatively regulate viral and host genome transcriptome.

There are several thousands of different HTLV-1-infected T-cell clones in an infected individual³⁸. After long-term clinical latency, a specific clone may undergo malignant transformation, causing the syndrome of ATL. A key question that remains is how a certain clone is selected as an ATL clone from various infected clones. Previous reports demonstrated that the HTLV-1 infected clones harboring proviruses near cancer-related genes are preferentially selected in ATL cells^29,39, indicating that aberrant host genome transcription caused by viral integration may contribute to the multistep oncogenic process induced by HTLV-1 infection. HTLV-1 contains a CTCF-binding site and therefore viral integration generates an ectopic CTCF-binding site in the host genome¹¹, which induces deregulation of host gene transcription via chromatin looping as Melamed A et al reported²⁶. We demonstrate here that HTLV-1 generates an ectopic enhancer region in addition to the viral CTCF insulator region. These findings indicate that the HTLV-1 enhancer can induce a distinct alteration of the host transcriptome via chromatin looping²⁶, and thereby upregulates cancer-related genes near ISs which might contribute to the preferential selection of a specific infected cell for clonal expansion during the early phase of leukemogenesis. This mechanism is similar to how endogenous retroviruses might contribute to the development of acute myeloid leukemia⁴⁰.

Mobile DNA elements, including endogenous retroviruses or foreign DNA elements introduced by exogenous retroviruses, can be dangerous for the host cell because they disturb cellular genomic homeostasis. Mammalian cells have an evolutionally acquired host defense system that silences such elements in the genomic DNA. For example, the murine leukemia virus (MLV) is silenced by Trim28—a well-characterized transcriptional co-repressor⁴¹—and ZFP809 to prevent further viral spread in embryonic stem cells¹³. Although little is known regarding the precise molecular mechanisms behind the silencing of HTLV-1 provirus in the host genome, the HTLV-1 5’LTR is frequently silenced by DNA methylation, histone modifications⁴², or only transcribed intermittently^32,33. This suggests that a host defense mechanism is playing a role in selecting infected clones with silenced HTLV-1 proviral DNA. As a result, there is no detectable viremia in the serum of HTLV-1-infected individuals. However, the virus maintains the ability to reactivate viral transcription when the virus needs to induce de novo infection from an infected host to an uninfected host. We showed here that HTLV-1 recruits the host transcription factors SRF and ELK-1 to an NFR in proviral DNA to sustain chromatin openness and proviral transcription in host cells. This molecular mechanism possibly enables the virus to be latent but at the same time allows it to maintain the ability to reactivate viral expression when infected cells need to induce de novo infection from the infected to the uninfected host.

HTLV-1 has co-existed with humans for the past 20,000–30,000 years⁴³. The virus may have evolved this strategy—the presence of an internal insulator and enhancer region in the provirus—to achieve persistent infection under pressure from the host system to silence foreign DNA elements as well as from the host immune response. Usage of lentiviral/retroviral vectors for gene therapy or for the generation of induced pluripotent stem (iPS) cells has been under intense research and development⁴⁴. Lentiviral and retroviral vectors integrate into host genomic DNA and form a provirus in the target cells; however, the provirus tends to be silenced by host defense mechanisms as described above. Various efforts have been made to optimize the lentiviral and retroviral vectors to prevent the silencing of the integrated provirus, such as the introduction of insulator or enhancer elements^44,45. Retrovirus vector insertion can trigger deregulated cell proliferation, most likely driven by the activity of retrovirus enhancers on cancer-related genes⁴⁶. It is surprising that an exogenous virus HTLV-1 has by itself evolved a similar system, obtaining an insulator, an enhancer, and a chromatin-opening element in the retroviral genome. This experiment of nature may provide insights into how an exogenous retrovirus achieves persistent infection in humans and also how to tackle the silencing of foreign DNA elements to maintain chromatin openness and transgene transcription without causing the transformation of host cells.

In conclusion, we have analyzed the HTLV-1 provirus integrated into the host genome with high resolution and efficiency using the HTLV-1-DNA-capture sequencing approach and discovered an internal viral enhancer in the HTLV-1 genome. This finding provides clues to help solve several long-lasting questions related to HTLV-1 persistence and pathogenesis. Viral DNA-capture-seq approaches can be applied to studies aiming to understand the transcriptional regulatory mechanism of other oncogenic viruses integrated into the host cellular genomic DNA.

Methods

Ethics statement and patient blood samples

All protocols involving human subjects were reviewed and approved by the Kumamoto University Institutional Review Board (approval number 263). The study was carried out in accordance with the guidelines proposed in the Declaration of Helsinki. Informed written consent was obtained from all subjects in this study. Peripheral blood mononuclear cells (PBMCs) were isolated from whole blood within 24 h of sample collection using Ficoll-Paque (GE Healthcare Life Sciences, Marlborough, MA) according to the manufacturer’s instructions. Characteristics of clinical samples were summarized in Supplementary Table 1.

Cell culture

JET cells⁴⁷ are Jurkat cells expressing tdTomato under the control of five times tandem repeat of Tax-responsive element (TRE). ED¹⁶, 293 T, Jurkat, and JET cells infected with WT or mutant HTLV-1 molecular clones were cultured in RPMI-1640 medium (Thermo Fisher Scientific, Waltham, MA) supplemented with 10% fetal bovine serum (FBS), 100 U/mL penicillin, and 100 μg/mL streptomycin. TBX-4B cells¹⁸ were cultured in RPMI-1640 supplemented with 20% FBS, interleukin-2 (200 U/mL; PeproTech, Cranbury, NJ), 100 U/mL penicillin, and 100 μg/mL streptomycin.

Generation of reporter constructs

The HBZ promoter, 3′LTR300²⁰, and 5′LTR were amplified from ED cells. The NFR was amplified from ED cells and the NFR mutant was generated by gBlocks^® Gene Fragments (Integrated DNA Technologies, Coralville, IA). Using XhoI and HindIII restriction sites, each promoter construct was inserted into pGL4-basic (Promega, Madison, WI) which includes the luciferase reporter gene. The NFR was inserted into pGL4-3′LTR300 or pGL4-5′LTR using BamHI or KpnI restriction sites while the NFR mutant was inserted into pGL4-3′LTR300 using the BamHI restriction site. NFR-CTCF fragments were cloned into pGL4-3′LTR300 or pGL4-5′LTR by NEBuilder HiFi DNA Assembly Master Mix (New England Biolabs, Ipswich, MA). Primers associated with each construct and the NFR mutant are listed in Supplementary Table 3.

mRNA-seq and qRT-PCR

RNA was extracted using the RNeasy Mini Kit (Qiagen, Hilden, Germany) according to the manufacturer’s instructions and treated with DNase. For mRNA-seq, mRNA libraries were prepared using NEBNext^® Ultra^TM II Directional RNA Library Prep Kit for Illumina^® Multiplex Oligos for Illumina (New England Biolabs) according to the manufacturer’s instructions. Libraries were run as 75-cycle-single end reads on a NextSeq 550 (Illumina, San Diego, CA) using a high-output flow cell. cDNA was synthesized using ReverTra Ace^® qPCR RT Master Mix (Toyobo, Osaka, Japan) according to the manufacturer’s instructions. qPCR was performed using Thunderbird SYBR qPCR mix (Toyobo) and run on an Applied Biosystems^® StepOnePlus^TM Real-Time PCR System (Thermo Fisher Scientific); primers used are listed in Supplementary Table 4.

Preparation and culture of HTLV-1-infected cells in vitro

293 T cells were transfected with a wt or enhancer-mutated HTLV-1 molecular clone⁴⁸ by polyethylenimine (PEI) and then irradiated with 30 Gy. The irradiated 293 T cells were co-cultured with JET cells for 3 days⁴⁷, after which tdTomato-positive cells were sorted by FACS Aria^TM (Becton, Dickinson and Company, Franklin Lakes, NJ), and cultured in RPMI supplemented with 10% FBS, 100 U/mL penicillin, and 100 μg/mL streptomycin for 2 weeks.

Proviral load (PVL) measurement

We estimated the number of infected cells by quantifying the copy number of the tax gene normalized to the copy number of the ALB gene by using digital droplet PCR as previously described but with minor modifications¹⁵. PVL was calculated as follows, PVL (%) = [(copy number of tax)/(copy number of albumin)/2] × 100. Primer sequences are listed in Supplementary Table 4.

HTLV-1 DNA-capture-seq

HTLV-1 DNA-capture-seq was performed as previously described¹⁵ with minor modifications. Briefly, 1 µg genomic DNA was fragmented by sonication using a Picoruptor (Diagenode s.a., Liège, Belgium) to produce 300–500-bp fragments. The DNA library was generated using a NEBNext Ultra II DNA Library Prep Kit for Illumina and Multiplex Oligos for Illumina (New England Biolabs). DNA-seq libraries were used for HTLV-1 sequence enrichment with HTLV-1 specific probes, after which enriched libraries were amplified by additional PCR. Enriched libraries were quantified using P5 and P7 primers and then sequenced via Illumina MiSeq or NextSeq.

MNase assay and MNase-seq

Cells (1.0 × 10⁶ for cell lines or 2.0 × 10⁶ for patient PBMCs) were lysed using cell lysis buffer (0.05% Triton X-100, 2 mM PMSF, 5 mM sodium butyrate, 100× protease inhibitor cocktail) or PBMC lysis buffer (10 mM Tris-HCl pH 7.4, 10 mM NaCl, 3 mM MgCl₂, 0.1% Nonidet-P40). Extracted nuclei were digested by MNase (TaKaRa Bio, Kusatsu, Japan) for 5–20 minutes at 37 °C after which the reaction was stopped by the addition of 20 mM ethylenediaminetetraacetic acid (EDTA). After deproteination with proteinase K solution (Nacalai Tesque, Kyoto, Japan), MNase digestion samples were purified using a PCR Purification Kit (Qiagen). We confirmed constant MNase treatment among different samples by gel electrophoresis after MNase digestion. MNase-digested DNA and input DNA were measured either by QX200 droplet digital PCR system (BIO-RAD, Hercules, CA) or by next-generation sequencing (NGS). Primer sequences for ddPCR are shown in Supplementary Table 6. For NGS, MNase-seq libraries were prepared by the NEBNext Ultra II DNA Library Prep Kit for Illumina and Multiplex Oligos for Illumina (New England Biolabs), after which the efficiency was quantified using P5 and P7 primers and then sequenced using Illumina MiSeq. We prepared fragmented DNAs by sonication using a Picoruptor (Diagenode s.a., Liège, Belgium) and use them as input DNAs. The degree of nucleosome freeness which was calculated as the MNase-seq value normalized to the input DNA value.

ChIP-seq

ChIP assays were performed using the SimpleChIP^® Enzymatic Chromatin IP Kit (Cell Signaling Technology, Danvers, MA) according to the manufacturer’s instructions. Briefly, cells (4 × 10⁶) were fixed in 1% formaldehyde for 10 min at room temperature, quenched in glycine solution, and washed in ice-cold PBS. Nuclei were extracted by lysis buffer (buffer A) and then samples were digested by MNase for 20 min at 37 °C and sonicated for 30 s on and 30 s off for 5–8 min using Bioruptor UCD-300 (Tosyodenki, Kanagawa, Japan) to break the nuclear membrane. Extracted chromatin was immunoprecipitated using anti-H3K27Ac (#07-360; Millipore, Burlington, MA), H3K4me1 (#ab8895; Abcam, UK, England), H3K4me2 (#ab7766, Abcam), SRF (#5147; Cell Signaling Technology), and ELK-1 (#ab32106; Abcam) antibodies. All antibodies were used at a 1:250 dilution. ChIP sample libraries were prepared by NEBNext Ultra II DNA Library Prep Kit for Illumina and Multiplex Oligos for Illumina (New England Biolabs), after which the efficiency was quantified using P5 and P7 primers and sequenced using Illumina MiSeq or NextSeq. To increase the detection sensitivity, we analyzed with the HTLV-1 DNA-capture method as reported previously¹⁵. Briefly, after library preparation, we mixed libraries and biotinylated DNA probes for the whole HTLV-1 provirus, and then HTLV-1 sequences were enriched by streptavidin. The enriched libraries were analyzed by Illumina MiSeq or NextSeq.

Luciferase reporter assays

Jurkat cells (2 × 10⁵) were harvested 24 h after transfection with 1 μg of each reporter construct, using 2 μl of Turbofect Transfection Reagent (Thermo Fisher Scientific). Luciferase assays were then performed using the Dual-Glo Luciferase Assay System (Promega) according to the manufacturer’s instructions, and luminescence was detected using GloMax^® 20/20 Luminometer (Promega).

NET-CAGE

Nascent RNAs were extracted from the nuclei of ED cells and TBX-4B cells following the previously described²³. NET-CAGE libraries were generated using the CAGE library preparation kit (K.K. DNAFORM, Yokohama, Japan) according to the manufacturer’s instructions. Briefly, cDNA was synthesized from 5 μg nascent RNAs. The 5′cap-structures of nascent RNAs were labeled by 4 μl of 10 mM biotin hydrazide for the cap-trapping step. After removing the Remaining RNA fragments without 5’cap structure by RNaseONE enzyme, enriched cDNA by cap-trapping was used for linker ligation and library generation. NET-CAGE Libraries were quantified by qPCR and sequenced using Illumina NextSeq.

EMSA (electrophoretic mobility shift assay)

293 T cells (2 × 10⁶) were harvested 24 h after transfection with 2 μg of pcDNA3-myc-SRF⁴⁹ and 2 μg of pCGN-ELK-1 (Addgene, Watertown, MA) using 16 μl of HilyMax (Dojindo Laboratories, Kumamoto, Japan). Cell lysates were extracted in 500 μl of cell lysis buffer (10 mM Tris-HCl pH 8.0, 60 mM KCl, 1 mM EDTA, 1 mM DTT, 100 µM PMSF, 0.1% NP-40) with 5 min incubation on ice. After cell lysis, nuclear lysates were extracted in 100 μl of nuclear extraction buffer (20 mM Tris-HCl pH 8.0, 420 mM NaCl, 1.5 mM MgCl₂, 0.2 mM EDTA, 25 % glycerol) with 10 min incubation on ice and then samples were sonicated for 20 on and 30 s off for 17 min using Bioruptor UCD-300 (Tosyodenki) to break overfilled DNA. EMSA was then performed with the 1 μl of extracted nuclear lysates, biotin-labeled NFR-wt probe, and NFR-wt or mut unlabeled probes using Perfect NT Gel which is a 3–12% gradient polyacrylamide gel (#NTH-5X5HP; DRC, Tokyo, Japan) and the LightShift Chemiluminescent EMSA Kit (#20148; Thermo Fisher Scientific) according to the manufacturer’s instructions. Nuclear lysates were mixed with 50 fmol biotin-labeled probes and 1 µg each of the anti-SRF (#5147; Cell Signaling Technology) and anti-ELK-1 (#ab32106; Abcam) antibodies. For the competition assay, NFR-wt or -mut unlabeled competitor probes (5, 10, and 15 pmol) were added in the mixture of nuclear lysates and biotin-labeled NFR-wt probe. Probe sequences are listed in Supplementary Table 5.

p19 ELISA

293 T cells (2 × 10⁵) were transfected with HTLV-1-wt or mut molecular clone (0.25, 0.5, and 1 μg) using 3 μl of HilyMax (Dojindo Laboratories). After 24 h, the supernatants were collected and measured p19 presence by RETROtek HTLV p19 Antigen ELISA (ZeptoMetrix Corporation, Buffalo, NY) following the manufacturer’s instruction.

Apoptosis analysis

JET cells infected with HTLV-1-wt or mut molecular clone were stimulated with 100 ng/ml PMA and 2 µM Ionomycin and incubated for 24 h. After incubation, apoptotic cells were stained with annexin V by MEBCYTO^® Apoptosis Kit (MBL, Nagoya, Japan) and detected by flow cytometry using BD FACSVerse^TM (Becton, Dickinson and Company). Flow cytometry data were analyzed using FlowJo ^TM (Becton, Dickinson and Company). Gating strategies for annexin V⁺ cells were shown in Supplementary Fig. 8

CRISPR/Cas9 mutagenesis

Guide sequences were designed with both edges of NFR in target and cloned into the pX330-U6-Chimeric BB-CBh-hSpCas9 plasmid (pX330; Addgene, 42230) as previously described⁵⁰. The oligonucleotides for constructing guide sequence were listed in Supplementary Table 7. HTLV-1-wt-infected JET clone (2 × 10⁶) was co-transfected with each 3 μg of two pX330 plasmids for each NFR edge, 1.5 μg of an expression vector with puromycin resistance gene and 3 μg of NFR_mut_3 cassette plasmid for HDR by electroporation using NEPA21 (NEPAGENE, Ichikawa, Japan) and 2 mm gap cuvette (EC-002S, NEPAGENE). The electroporation program was following; 275 V, 1 ms, six times, and a 50 ms interval for poring pulse, and 20 V, 50 ms, three times, and a 50 ms interval for transfer pulse. Transfected cells were selected by 1 μg/mL puromycin with 24 h incubation. After transfected cells were recovered cell damages for 3 days, limiting dilution was performed to get a single clone. CRISPR/Cas9 mediated mutant clone was confirmed the sequence which converted wt to mut by Sanger sequencing.

Single-cell RNA-seq analysis

Single-cell data acquisition was performed in a previous study⁵¹ and the sequencing data were obtained from the European Nucleotide Archive (ENA) (https://www.ebi.ac.uk/ena/browser/home) with the following accession numbers for each sample: ATL_1, ERX6294562; ATL_2, ERX6294563; ATL_8, ERX6294567; ATL_9, ERX6294566; and ATL_10, ERX6294565. Cell Ranger Single-Cell Software Suite (v3.1.0, 10x Genomics, Pleasanton, CA) was used to perform sequence alignment against a modified hg38 human reference genome which contains the HTLV-1 genome (Genbank accession no. AB513134) as a separate chromosome. The resulting gene-cell barcode matrix was imported into R (v4.0.3) and analyzed using Seurat (v4.0) according to the vignette on Seurat’s webpage available here (https://satijalab.org/seurat/articles/pbmc3k_tutorial.html). Clusters were annotated based on examination of known marker genes for each PBMC subsets. In this case, CD4 T-cells, CD8/NK cells, B cells and myeloid cells are CD3D⁺CD4⁺, CD3D⁺CD8A⁺NKG7⁺, CD79A⁺CD19⁺, and CD14⁺FCGR3A⁺ respectively⁵². T-cell clones are identified based on TCR information with the ATL clone defined as the most expanded CD4⁺ T-cell clone in the sample.

Western blot

We generated wt- or mutated-tax (tagged c-Myc) pcDNA3.1(-)-c-Myc vectors based on pcDNA3.1(-) (Invitrogen, Waltham, MA). 293 T cells (2 × 10⁶) were harvested 24 h after transfection with 6 μg of wt- or mutated-tax (tagged c-Myc) pcDNA3.1(-)-c-Myc vectors using 16 μl of HilyMax (Dojindo Laboratories). After cell lysis in 200 μl of RIPA buffer (50 mM Tris-HCl pH 7.4, 150 mM NaCl, 0.1% SDS, 1% Triton X-100, 1% sodium deoxycholate, 1 mM EDTA) with protease inhibitor cocktail and phosphatase inhibitor cocktail, 30 μg of cell proteins were separated on the precast gel for SDS-PAGE (197-15011; FUJIFILM, Osaka, Japan). Separated proteins were blotted onto polyvinylidene difluoride membrane (5400412; ATTO, Tokyo, Japan). About 0.1 μg/ml anti-Myc antibody (M192-3; MBL, Tokyo, Japan) and 0.1 μg/ml anti-Actin (C-2) antibody (sc-8432; Santa Cruz Biotechnology, Dallas, TX) were used as primary antibodies staining. secondary reactions with HRP were performed using Pierce^TM Fast Western Blot Kit, ECL Substrate (35055; Thermo Fisher Scientific) according to the manufacturer’s instruction. Chemiluminescent detection was performed on a ChemiDoc^TM Touch Imaging system (BIO-RAD).

Bioinformatic analysis

Prediction of transcription factors which bind to the NFR sequence was performed using TFBIND (https://tfbind.hgc.jp). Analysis of next-generation sequencing data was performed as follows. First, the quality of raw FASTQ files were checked with FastQC (v 0.10.0) followed by adapter trimming and removal of poor-quality reads using cutadapt (v 1.18) and PRINSEQ (v 0.20.4) respectively. For ChIP-seq analysis, peak calling was performed using MACS (v 1.4.2) as described previously¹¹. For the HTLV-1-DNA-seq data, alignment to reference genome was performed using BWA (v 0.7.12) followed by viral integration site and clonal abundance analysis using samtools (v 1.11), picard (v 2.0.1), and in-house perl scripts as we previously reported¹⁵. RefSeq gene data were obtained from UCSC tables (https://genome.ucsc.edu/). The relationship between viral integration site and host genes or epigenetic microenvironment were analyzed using the R package hiAnnotator (http://github.com/malnirav/hiAnnotator) as described previously⁵³. Gene expression for bulk RNA-seq data of cell lines (ED and TBX-4B), parent (JET), wild type, and mutated clones (see Supplementary Table 2) was quantified using kallisto⁵⁴ and exported into R for differential expression analysis using the R package DESeq2⁵⁵ (https://github.com/mikelove/DESeq2). Data were filtered to remove genes with low counts (<10) followed by log-transformation and visualization on a PCA plot. The same RNA-seq data were also aligned to a reference genome using STAR (v 2.7.3) andfor visualization on IGV (v 2.8.0). Data for doxycycline-induced Tax- or HBZ-expressing Jurkat cells were obtained from SRA1049749²⁷. Fold change of gene expression by tax or HBZ was calculated with respect to the non-induced cells (i.e. Doxycycline-induced HBZ-Jurkat against non-induced HBZ-Jurkat).

Statistical analysis

Data were analyzed using a chi-squared test with GraphPad Prism 7 software (GraphPad Software Inc., La Jolla, CA) unless otherwise stated. Statistical significance was defined as P < 0.05.

Reporting summary

Further information on research design is available in the Nature Research Reporting Summary linked to this article.

Data availability

All NGS sequencing data generated in this study have been deposited in the DDBJ Sequence Read Archive (DRA), which is associated with DNA DataBank of Japan (DDBJ) under accession code DRA013478 (https://ddbj.nig.ac.jp/resource/sra-submission/DRA013478). Processed data have been deposited in the DRA under the accession number DRA013591 (https://ddbj.nig.ac.jp/resource/sra-submission/DRA013591) and DRA013588 (https://ddbj.nig.ac.jp/resource/sra-submission/DRA013588). Peak call data of ChIP-seq have been deposited in the Genomic Expression Archive (GEA), which is associated with DDBJ under an accession number E-GEAD-481. Raw experimental data such as luciferase assay, EMSA, ELISA, qRT-PCR, PVL measurement, cell apoptosis assay, MNase assay, PCA plots, TPM of RNA-seq, cell clustering assay of scRNA-seq, integration site distribution analysis, western blots are available from Source Data file.

Code availability

The source code to reproduce our analysis has been uploaded to GitHub and linked to Zenodo: https://doi.org/10.5281/zenodo.6361824.

References

Uchiyama, T., Yodoi, J., Sagawa, K., Takatsuki, K. & Uchino, H. Adult T-cell leukemia: clinical and hematologic features of 16 cases. Blood 50, 481–492 (1977).
Article CAS PubMed Google Scholar
Poiesz, B. J. et al. Detection and isolation of type C retrovirus particles from fresh and cultured lymphocytes of a patient with cutaneous T-cell lymphoma. Proc. Natl Acad. Sci. USA 77, 7415–7419 (1980).
Article ADS CAS PubMed PubMed Central Google Scholar
Hinuma, Y. et al. Adult T-cell leukemia: antigen in an ATL cell line and detection of antibodies to the antigen in human sera. Proc. Natl Acad. Sci. USA 78, 6476–6480 (1981).
Article ADS CAS PubMed PubMed Central Google Scholar
Matsuoka, M. & Jeang, K. T. Human T-cell leukaemia virus type 1 (HTLV-1) infectivity and cellular transformation. Nat. Rev. Cancer 7, 270–280 (2007).
Article CAS PubMed Google Scholar
Bangham, C. R. M. Human T cell leukemia virus type 1: persistence and pathogenesis. Annu. Rev. Immunol. 36, 43–71 (2018).
Article CAS PubMed Google Scholar
Yoshida, M. Multiple viral strategies of HTLV-1 for dysregulation of cell growth control. Annu. Rev. Immunol. 19, 475–496 (2001).
Article CAS PubMed Google Scholar
Giam, C. Z. & Semmes, O. J. HTLV-1 infection and adult T-cell leukemia/lymphoma-A tale of two proteins: Tax and HBZ. Viruses https://doi.org/10.3390/v8060161 (2016).
Satou, Y. et al. HTLV-1 bZIP factor induces T-cell lymphoma and systemic inflammation in vivo. PLoS Pathog. 7, e1001274 (2011).
Article CAS PubMed PubMed Central Google Scholar
Satou, Y., Yasunaga, J., Yoshida, M. & Matsuoka, M. HTLV-I basic leucine zipper factor gene mRNA supports proliferation of adult T cell leukemia cells. Proc. Natl Acad. Sci. USA 103, 720–725 (2006).
Article ADS CAS PubMed PubMed Central Google Scholar
Usui, T. et al. Characteristic expression of HTLV-1 basic zipper factor (HBZ) transcripts in HTLV-1 provirus-positive cells. Retrovirology 5, 34 (2008).
Article PubMed PubMed Central CAS Google Scholar
Satou, Y. et al. The retrovirus HTLV-1 inserts an ectopic CTCF-binding site into the human genome. Proc. Natl Acad. Sci. USA 113, 3054–3059 (2016).
Article ADS CAS PubMed PubMed Central Google Scholar
Wolf, D. & Goff, S. P. TRIM28 mediates primer binding site-targeted silencing of murine leukemia virus in embryonic cells. Cell 131, 46–57 (2007).
Article CAS PubMed Google Scholar
Wolf, D. & Goff, S. P. Embryonic stem cells use ZFP809 to silence retroviral DNAs. Nature 458, 1201–1204 (2009).
Article ADS CAS PubMed PubMed Central Google Scholar
Miyazato, P. et al. Application of targeted enrichment to next-generation sequencing of retroviruses integrated into the host human genome. Sci. Rep. 6, 28324 (2016).
Article ADS CAS PubMed PubMed Central Google Scholar
Katsuya, H. et al. The nature of the HTLV-1 provirus in naturally infected individuals analyzed by the viral DNA-capture-seq approach. Cell Rep. 29, 724–735 (2019).
Article CAS PubMed Google Scholar
Maeda, M. et al. Origin of human T-lymphotrophic virus I-positive T cell lines in adult T cell leukemia. Analysis of T cell receptor gene rearrangement. J. Exp. Med. 162, 2169–2174 (1985).
Article CAS PubMed Google Scholar
Takeda, S. et al. Genetic and epigenetic inactivation of tax gene in adult T-cell leukemia cells. Int. J. Cancer 109, 559–567 (2004).
Article CAS PubMed Google Scholar
Cook, L. B., Rowan, A. G., Melamed, A., Taylor, G. P. & Bangham, C. R. HTLV-1-infected T cells contain a single integrated provirus in natural infection. Blood 120, 3488–3490 (2012).
Article CAS PubMed PubMed Central Google Scholar
Hanon, E. et al. Abundant tax protein expression in CD4+ T cells infected with human T-cell lymphotropic virus type I (HTLV-I) is prevented by cytotoxic T lymphocytes. Blood 95, 1386–1392 (2000).
Article CAS PubMed Google Scholar
Yoshida, M., Satou, Y., Yasunaga, J., Fujisawa, J. & Matsuoka, M. Transcriptional control of spliced and unspliced human T-cell leukemia virus type 1 bZIP factor (HBZ) gene. J. Virol. 82, 9359–9368 (2008).
Article CAS PubMed PubMed Central Google Scholar
Zhou, V. W., Goren, A. & Bernstein, B. E. Charting histone modifications and the functional organization of mammalian genomes. Nat. Rev. Genet. 12, 7–18 (2011).
Article PubMed CAS Google Scholar
Kim, T. K. et al. Widespread transcription at neuronal activity-regulated enhancers. Nature 465, 182–187 (2010).
Article ADS CAS PubMed PubMed Central Google Scholar
Hirabayashi, S. et al. NET-CAGE characterizes the dynamics and topology of human transcribed cis-regulatory elements. Nat. Genet. 51, 1369–1379 (2019).
Article CAS PubMed Google Scholar
Suzuki, T., Hirai, H., Fujisawa, J., Fujita, T. & Yoshida, M. A trans-activator Tax of human T-cell leukemia virus type 1 binds to NF-kappa B p50 and serum response factor (SRF) and associates with enhancer DNAs of the NF-kappa B site and CArG box. Oncogene 8, 2391–2397 (1993).
CAS PubMed Google Scholar
Tanaka-Nakanishi, A., Yasunaga, J., Takai, K. & Matsuoka, M. HTLV-1 bZIP factor suppresses apoptosis by attenuating the function of FoxO3a and altering its localization. Cancer Res. 74, 188–200 (2014).
Article CAS PubMed Google Scholar
Melamed, A. et al. The human leukemia virus HTLV-1 alters the structure and transcription of host chromatin in cis. Elife https://doi.org/10.7554/eLife.36245 (2018).
Vandermeulen, C. et al. The HTLV-1 viral oncoproteins Tax and HBZ reprogram the cellular mRNA splicing landscape. PLoS Pathog. 17, e1009919 (2021).
Article CAS PubMed PubMed Central Google Scholar
Kataoka, K. et al. Integrated molecular analysis of adult T cell leukemia/lymphoma. Nat. Genet. 47, 1304–1315 (2015).
Article CAS PubMed Google Scholar
Rosewick, N. et al. Cis-perturbation of cancer drivers by the HTLV-1/BLV proviruses is an early determinant of leukemogenesis. Nat. Commun. 8, 15264 (2017).
Article ADS CAS PubMed PubMed Central Google Scholar
Fujisawa, J., Seiki, M., Kiyokawa, T. & Yoshida, M. Functional activation of the long terminal repeat of human T-cell leukemia virus type I by a trans-acting factor. Proc. Natl Acad. Sci. USA 82, 2277–2281 (1985).
Article ADS CAS PubMed PubMed Central Google Scholar
Fujisawa, J., Seiki, M., Sato, M. & Yoshida, M. A transcriptional enhancer sequence of HTLV-I is responsible for trans-activation mediated by p40 chi HTLV-I. EMBO J. 5, 713–718 (1986).
Article CAS PubMed PubMed Central Google Scholar
Billman, M. R., Rueda, D. & Bangham, C. R. M. Single-cell heterogeneity and cell-cycle-related viral gene bursts in the human leukaemia virus HTLV-1. Wellcome Open Res. 2, 87 (2017).
Article PubMed PubMed Central CAS Google Scholar
Mahgoub, M. et al. Sporadic on/off switching of HTLV-1 Tax expression is crucial to maintain the whole population of virus-induced leukemic cells. Proc. Natl Acad. Sci. USA 115, E1269–E1278 (2018).
Article CAS PubMed PubMed Central Google Scholar
Sodroski, J., Patarca, R., Rosen, C., Wong-Staal, F. & Haseltine, W. Location of the trans-activating region on the genome of human T-cell lymphotropic virus type III. Science 229, 74–77 (1985).
Article ADS CAS PubMed Google Scholar
Arnold, J. et al. Enhancement of infectivity and persistence in vivo by HBZ, a natural antisense coded protein of HTLV-1. Blood 107, 3976–3982 (2006).
Article CAS PubMed PubMed Central Google Scholar
Tamiya, S. et al. Two types of defective human T-lymphotropic virus type I provirus in adult T-cell leukemia. Blood 88, 3065–3073 (1996).
Article CAS PubMed Google Scholar
Miyazaki, M. et al. Preferential selection of human T-cell leukemia virus type 1 provirus lacking the 5’ long terminal repeat during oncogenesis. J. Virol. 81, 5714–5723 (2007).
Article CAS PubMed PubMed Central Google Scholar
Gillet, N. A. et al. The host genomic environment of the provirus determines the abundance of HTLV-1-infected T-cell clones. Blood 117, 3113–3122 (2011).
Article CAS PubMed PubMed Central Google Scholar
Cook, L. B. et al. The role of HTLV-1 clonality, proviral structure, and genomic integration site in adult T-cell leukemia/lymphoma. Blood 123, 3925–3931 (2014).
Article CAS PubMed PubMed Central Google Scholar
Deniz, O. et al. Endogenous retroviruses are a source of enhancers with oncogenic potential in acute myeloid leukaemia. Nat. Commun. 11, 3506 (2020).
Article ADS CAS PubMed PubMed Central Google Scholar
O’Geen, H. et al. Genome-wide analysis of KAP1 binding suggests autoregulation of KRAB-ZNFs. PLoS Genet. 3, e89 (2007).
Article PubMed PubMed Central CAS Google Scholar
Taniguchi, Y. et al. Silencing of human T-cell leukemia virus type I gene transcription by epigenetic mechanisms. Retrovirology 2, 64 (2005).
Article PubMed PubMed Central CAS Google Scholar
Verdonck, K. et al. Human T-lymphotropic virus 1: recent knowledge about an ancient infection. Lancet Infect. Dis. 7, 266–281 (2007).
Article CAS PubMed Google Scholar
David, R. M. & Doherty, A. T. Viral vectors: the road to reducing genotoxicity. Toxicol. Sci. 155, 315–325 (2017).
Article CAS PubMed Google Scholar
Liu, M. et al. Genomic discovery of potent chromatin insulators for human gene therapy. Nat. Biotechnol. 33, 198–203 (2015).
Article PubMed CAS Google Scholar
Hacein-Bey-Abina, S. et al. LMO2-associated clonal T cell proliferation in two patients after gene therapy for SCID-X1. Science 302, 415–419 (2003).
Article ADS CAS PubMed Google Scholar
Furuta, R. et al. Human T-cell leukemia virus type 1 infects multiple lineage hematopoietic cells in vivo. PLoS Pathog. 13, e1006722 (2017).
Article PubMed PubMed Central CAS Google Scholar
Mitchell, M. S. et al. Phenotypic and genotypic comparisons of human T-cell leukemia virus type 1 reverse transcriptases from infected T-cell lines and patient samples. J. Virol. 81, 4422–4428 (2007).
Article CAS PubMed PubMed Central Google Scholar
Matsuzaki, K. et al. PML-nuclear bodies are involved in cellular serum response. Genes Cells 8, 275–286 (2003).
Article CAS PubMed Google Scholar
Ran, F. A. et al. Genome engineering using the CRISPR-Cas9 system. Nat. Protoc. 8, 2281–2308 (2013).
Article CAS PubMed PubMed Central Google Scholar
Tan, B. J. et al. HTLV-1 infection promotes excessive T cell activation and transformation into adult T cell leukemia/lymphoma. J. Clin. Invest. https://doi.org/10.1172/JCI150472 (2021).
Zheng, G. X. et al. Massively parallel digital transcriptional profiling of single cells. Nat. Commun. 8, 14049 (2017).
Article ADS CAS PubMed PubMed Central Google Scholar
Satou, Y. et al. Dynamics and mechanisms of clonal expansion of HIV-1-infected cells in a humanized mouse model. Sci. Rep. 7, 6913 (2017).
Article ADS PubMed PubMed Central CAS Google Scholar
Bray, N. L., Pimentel, H., Melsted, P. & Pachter, L. Near-optimal probabilistic RNA-seq quantification. Nat. Biotechnol. 34, 525–527 (2016).
Article CAS PubMed Google Scholar
Love, M. I., Huber, W. & Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15, 550 (2014).
Article PubMed PubMed Central CAS Google Scholar
Tsunoda, T. & Takagi, T. Estimating transcription factor bindability on DNA. Bioinformatics 15, 622–630 (1999).
Article CAS PubMed Google Scholar

Download references

Acknowledgements

We would like to thank M. Miura for providing the R program to perform quality checks of the index reads using the R program, M. Nakao for providing the SRF expression vector, N. Misawa, S. Nagaoka, and K. Sato for technical support and valuable discussion. We are also grateful to CRM. Bangham, S. Hino, and M. Ono for their critical reading of the manuscript. This study was supported by grants from the Japan Society for the Promotion of Science (JSPS) KAKENHI (JP20H03724, and JP18KK0230 to Y.S., 16KK0206 and JP18K16122 to H.K., JP18K08437 and JP18KK0452 to P.M., JP20K22783 and JP21K08494 to K.S.; JP21K15454 to M.M.) and Japan Agency for Medical Research and Development (AMED) (JP21jm0210074, JP21wm0325015, JP19fm0208012, JP21fk0410023, and JP21wm03250152to Y.S.) the Grant for Joint Research Project of the Institute of Medical Science, the University of Tokyo to Y.S., the grant from Kumamoto University Excellent Research Projects to Y.S., JST MIRAI (18077147) to Y.S., the program of the Joint Usage/Research Center for Developmental Medicine, Inter-University Research Network for Trans-Omics Medicine, Institute of Molecular Embryology and Genetics, Kumamoto University to Y.S. and Kumamoto University Fellowship for Excellent Graduate Students to M.M. The funders had no role in study design, data collection, data interpretation, or the discussion regarding submission for publication.

Author information

Saiful Islam
Present address: Viral Recombination Section, HIV Dynamics and Replication Program, National Cancer Institute, Frederick, MD, 21702, US
These authors contributed equally: Takaharu Ueno, Kazuaki Monde, Kenji Sugata.

Authors and Affiliations

Division of Genomics and Transcriptomics, Joint Research Center for Human Retrovirus Infection, Kumamoto University, Kumamoto, 860-8556, Japan
Misaki Matsuo, Kenji Sugata, Benjy Jek Yang Tan, Akhinur Rahman, Paola Miyazato, Kyosuke Uchiyama, Saiful Islam, Hiroo Katsuya & Yorifumi Satou
International Research Center for Medical Sciences (IRCMS), Kumamoto University, Kumamoto, 860-0811, Japan
Misaki Matsuo, Benjy Jek Yang Tan, Paola Miyazato, Kyosuke Uchiyama, Saiful Islam & Yorifumi Satou
Department of Microbiology, Kansai Medical University, Osaka, 573-1010, Japan
Takaharu Ueno, Shinsuke Nakajima & Jun-ichi Fujisawa
Department of Microbiology, Faculty of Life Sciences, Kumamoto University, Kumamoto, 860-8556, Japan
Kazuaki Monde
Division of Hematology, Respiratory Medicine and Oncology, Saga University, Saga, 849-8501, Japan
Hiroo Katsuya
Department of Hematology, Imamura General Hospital, Kagoshima, 890-0064, Japan
Masahito Tokunaga & Atae Utsunomiya
Department of Hematology, Rheumatology and Infectious Disease, Kumamoto University Hospital, Kumamoto, 860-8556, Japan
Kisato Nosaka
Cancer Center, Kumamoto University Hospital, Kumamoto, 860-8556, Japan
Kisato Nosaka
Division of Informative Clinical Sciences, Faculty of Life Sciences, Kumamoto University, Kumamoto, 862-0972, Japan
Hiroyuki Hata
Graduate School of Medical and Dental Sciences, Kagoshima University, Kagoshima, 890-8544, Japan
Atae Utsunomiya

Authors

Misaki Matsuo
View author publications
You can also search for this author in PubMed Google Scholar
Takaharu Ueno
View author publications
You can also search for this author in PubMed Google Scholar
Kazuaki Monde
View author publications
You can also search for this author in PubMed Google Scholar
Kenji Sugata
View author publications
You can also search for this author in PubMed Google Scholar
Benjy Jek Yang Tan
View author publications
You can also search for this author in PubMed Google Scholar
Akhinur Rahman
View author publications
You can also search for this author in PubMed Google Scholar
Paola Miyazato
View author publications
You can also search for this author in PubMed Google Scholar
Kyosuke Uchiyama
View author publications
You can also search for this author in PubMed Google Scholar
Saiful Islam
View author publications
You can also search for this author in PubMed Google Scholar
Hiroo Katsuya
View author publications
You can also search for this author in PubMed Google Scholar
Shinsuke Nakajima
View author publications
You can also search for this author in PubMed Google Scholar
Masahito Tokunaga
View author publications
You can also search for this author in PubMed Google Scholar
Kisato Nosaka
View author publications
You can also search for this author in PubMed Google Scholar
Hiroyuki Hata
View author publications
You can also search for this author in PubMed Google Scholar
Atae Utsunomiya
View author publications
You can also search for this author in PubMed Google Scholar
Jun-ichi Fujisawa
View author publications
You can also search for this author in PubMed Google Scholar
Yorifumi Satou
View author publications
You can also search for this author in PubMed Google Scholar

Contributions

M.M acquired funding for the project, designed and performed almost experiments, including epigenetic profiles, gene expression profiles, wt or mutant HTLV-1 profiles, functional analysis of the target proviral region, flow cytometry, genome editing, bioinformatic analysis, data curation, and wrote the paper. T.U. and K.M. established wt or mutant HTLV-1 infected clones for analysis of provirus profiles. K.S. acquired funding for the project, performed a functional analysis of the target proviral region, and wrote the paper. B.J.Y.T. performed a bioinformatic analysis of scRNA-seq and wrote the paper. A.R., K.U., and S.I. generated DNA libraries, performed DNA-captured-seq, and proviral load measurements. P.M. acquired funding for the project and performed NET-CAGE. H.K. acquired funding for the project and performed HTLV-1 integration-site (IS) analysis and IS characterization. M.T., K.N., and A.U. provided clinical samples. S.N. produced data for characterization of wt or mutant HTLV-1 infected cells. H.H. and J.F. supervised the project. Y.S. conceived and supervised the project, acquired funding for the project, performed data curation and bioinformatic analysis, and wrote the paper. All authors discussed the results and commented on the manuscript.

Corresponding author

Correspondence to Yorifumi Satou.

Ethics declarations

Competing interests

The authors declare no competing interests.

Peer review

Peer review information

Nature Communications thanks Angela Ciuffi and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.

Additional information

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Supplementary Information

Peer Review File

Reporting Summary

Source data

Source Data

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and permissions

About this article

Cite this article

Matsuo, M., Ueno, T., Monde, K. et al. Identification and characterization of a novel enhancer in the HTLV-1 proviral genome. Nat Commun 13, 2405 (2022). https://doi.org/10.1038/s41467-022-30029-9

Download citation

Received: 03 February 2021
Accepted: 07 April 2022
Published: 03 May 2022
DOI: https://doi.org/10.1038/s41467-022-30029-9

This article is cited by

Analysis of Programmed Cell Death-1 (PD-1) Gene Variations (re11568821 and rs41386349) in HTLV-1 Infection Using One Primer Pair and Proviral Load
- Yalda Amiri Hezave
- Zohreh Sharifi
- Fahime Ranjbar Kermani
Journal of Molecular Evolution (2023)
Viral, genetic, and immune factors in the oncogenesis of adult T-cell leukemia/lymphoma
- Jun-ichirou Yasunaga
International Journal of Hematology (2023)

Comments

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.