Anat Melamed

The host genomic environment of the provirus determines the abundance of HTLV-1-infected T-cell clones

Human T-lymphotropic virus type 1 (HTLV-1) persists by driving clonal proliferation of infected T lymphocytes. A high proviral load predisposes to HTLV-1-associated diseases. Yet the reasons for the variation within and between persons in the abundance of HTLV-1-infected clones remain unknown. We devised a high-throughput protocol to map the genomic location and quantify the abundance of > 91,000 unique insertion sites of the provirus from 61 HTLV-1(+) persons and > 2100 sites from in vitro infection. We show that a typical HTLV-1-infected host carries between 500 and 5000 unique insertion sites. We demonstrate that negative selection dominates during chronic infection, favoring establishment of proviruses integrated in transcriptionally silenced DNA: this selection is significantly stronger in asymptomatic carriers. We define a parameter, the oligoclonality index, to quantify clonality. The high proviral load characteristic of HTLV-1-associated inflammatory disease results from a larger number of unique insertion sites than in asymptomatic carriers and not, as previously thought, from a difference in clonality. The abundance of established HTLV-1 clones is determined by genomic features of the host DNA flanking the provirus. HTLV-1 clonal expansion in vivo is favored by orientation of the provirus in the same sense as the nearest host gene.

Genome-wide determinants of proviral targeting, clonal abundance and expression in natural HTLV-1 infection.

The regulation of proviral latency is a central problem in retrovirology. We postulate that the genomic integration site of human T lymphotropic virus type 1 (HTLV-1) determines the pattern of expression of the provirus, which in turn determines the abundance and pathogenic potential of infected T cell clones in vivo. We recently developed a high-throughput method for the genome-wide amplification, identification and quantification of proviral integration sites. Here, we used this protocol to test two hypotheses. First, that binding sites for transcription factors and chromatin remodelling factors in the genome flanking the proviral integration site of HTLV-1 are associated with integration targeting, spontaneous proviral expression, and in vivo clonal abundance. Second, that the transcriptional orientation of the HTLV-1 provirus relative to that of the nearest host gene determines spontaneous proviral expression and in vivo clonal abundance. Integration targeting was strongly associated with the presence of a binding site for specific host transcription factors, especially STAT1 and p53. The presence of the chromatin remodelling factors BRG1 and INI1 and certain host transcription factors either upstream or downstream of the provirus was associated respectively with silencing or spontaneous expression of the provirus. Cells expressing HTLV-1 Tax protein were significantly more frequent in clones of low abundance in vivo. We conclude that transcriptional interference and chromatin remodelling are critical determinants of proviral latency in natural HTLV-1 infection.

HTLV-1-infected T cells contain a single integrated provirus in natural infection

Human T lymphotropic virus type 1 (HTLV-1) appears to persist in the chronic phase of infection by driving oligoclonal proliferation of infected T cells. Our recent high-throughput sequencing study revealed a large number (often > 10(4)) of distinct proviral integration sites of HTLV-1 in each host that is greatly in excess of previous estimates. Here we use the highly sensitive, quantitative high-throughput sequencing protocol to show that circulating HTLV-1(+) clones in natural infection each contain a single integrated proviral copy. We conclude that a typical host possesses a large number of distinct HTLV-1-infected T-cell clones.

Estimating abundances of retroviral insertion sites from DNA fragment length data

MOTIVATION The relative abundance of retroviral insertions in a host genome is important in understanding the persistence and pathogenesis of both natural retroviral infections and retroviral gene therapy vectors. It could be estimated from a sample of cells if only the host genomic sites of retroviral insertions could be directly counted. When host genomic DNA is randomly broken via sonication and then amplified, amplicons of varying lengths are produced. The number of unique lengths of amplicons of an insertion site tends to increase according to its abundance, providing a basis for estimating relative abundance. However, as abundance increases amplicons of the same length arise by chance leading to a non-linear relation between the number of unique lengths and relative abundance. The difficulty in calibrating this relation is compounded by sample-specific variations in the relative frequencies of clones of each length. RESULTS A likelihood function is proposed for the discrete lengths observed in each of a collection of insertion sites and is maximized with a hybrid expectation-maximization algorithm. Patient data illustrate the method and simulations show that relative abundance can be estimated with little bias, but that variation in highly abundant sites can be large. In replicated patient samples, variation exceeds what the model implies-requiring adjustment as in Efron (2004) or using jackknife standard errors. Consequently, it is advantageous to collect replicate samples to strengthen inferences about relative abundance.

Quantification of HTLV-1 Clonality and TCR Diversity.

Estimation of immunological and microbiological diversity is vital to our understanding of infection and the immune response. For instance, what is the diversity of the T cell repertoire? These questions are partially addressed by high-throughput sequencing techniques that enable identification of immunological and microbiological “species” in a sample. Estimators of the number of unseen species are needed to estimate population diversity from sample diversity. Here we test five widely used non-parametric estimators, and develop and validate a novel method, DivE, to estimate species richness and distribution. We used three independent datasets: (i) viral populations from subjects infected with human T-lymphotropic virus type 1; (ii) T cell antigen receptor clonotype repertoires; and (iii) microbial data from infant faecal samples. When applied to datasets with rarefaction curves that did not plateau, existing estimators systematically increased with sample size. In contrast, DivE consistently and accurately estimated diversity for all datasets. We identify conditions that limit the application of DivE. We also show that DivE can be used to accurately estimate the underlying population frequency distribution. We have developed a novel method that is significantly more accurate than commonly used biodiversity estimators in microbiological and immunological populations.

HTLV-1 clonality in adult T-cell leukaemia and non-malignant HTLV-1 infection

Human T lymphotropic virus type 1 (HTLV-1) causes a range of chronic inflammatory diseases and an aggressive malignancy of T lymphocytes known as adult T-cell leukaemia/lymphoma (ATLL). A cardinal feature of HTLV-1 infection is the presence of expanded clones of HTLV-1-infected T cells, which may persist for decades. A high viral burden (proviral load) is associated with both the inflammatory and malignant diseases caused by HTLV-1, and it has been believed that the oligoclonal expansion of infected cells predisposes to these diseases. However, it is not understood what regulates the clonality of HTLV-1 in vivo, that is, the number and abundance of HTLV-1-infected T cell clones. We review recent advances in the understanding of HTLV-1 infection and disease that have come from high-throughput quantification and analysis of HTLV-1 clonality in natural infection.

The retrovirus HTLV-1 inserts an ectopic CTCF-binding site into the human genome

Significance The retrovirus human T-lymphotropic virus type 1 (HTLV-1) causes inflammatory and malignant diseases in humans. To maintain latency and avoid immune detection in vivo, HTLV-1 minimizes expression of genes on the plus-strand of the integrated provirus but allows constitutive expression of the minus-strand gene, which maintains clonal persistence. It is not understood how this gene expression is regulated. We show that CTCF, a master regulator of chromatin structure and gene expression, binds to HTLV-1, forms loops between the provirus and host genome, and alters expression of proviral and host genes. Because a typical HTLV-1–infected host carries >104 infected T-cell clones, each containing a provirus integrated in a different genomic site, CTCF binding gives HTLV-1 the potential to cause widespread abnormalities in the human genome. Human T-lymphotropic virus type 1 (HTLV-1) is a retrovirus that causes malignant and inflammatory diseases in ∼10% of infected people. A typical host has between 104 and 105 clones of HTLV-1–infected T lymphocytes, each clone distinguished by the genomic integration site of the single-copy HTLV-1 provirus. The HTLV-1 bZIP (HBZ) factor gene is constitutively expressed from the minus strand of the provirus, whereas plus-strand expression, required for viral propagation to uninfected cells, is suppressed or intermittent in vivo, allowing escape from host immune surveillance. It remains unknown what regulates this pattern of proviral transcription and latency. Here, we show that CTCF, a key regulator of chromatin structure and function, binds to the provirus at a sharp border in epigenetic modifications in the pX region of the HTLV-1 provirus in T cells naturally infected with HTLV-1. CTCF is a zinc-finger protein that binds to an insulator region in genomic DNA and plays a fundamental role in controlling higher order chromatin structure and gene expression in vertebrate cells. We show that CTCF bound to HTLV-1 acts as an enhancer blocker, regulates HTLV-1 mRNA splicing, and forms long-distance interactions with flanking host chromatin. CTCF-binding sites (CTCF-BSs) have been propagated throughout the genome by transposons in certain primate lineages, but CTCF binding has not previously been described in present-day exogenous retroviruses. The presence of an ectopic CTCF-BS introduced by the retrovirus in tens of thousands of genomic locations has the potential to cause widespread abnormalities in host cell chromatin structure and gene expression.

Clonality of HTLV-2 in Natural Infection

Human T-lymphotropic virus type 1 (HTLV-1) and type 2 (HTLV-2) both cause lifelong persistent infections, but differ in their clinical outcomes. HTLV-1 infection causes a chronic or acute T-lymphocytic malignancy in up to 5% of infected individuals whereas HTLV-2 has not been unequivocally linked to a T-cell malignancy. Virus-driven clonal proliferation of infected cells both in vitro and in vivo has been demonstrated in HTLV-1 infection. However, T-cell clonality in HTLV-2 infection has not been rigorously characterized. In this study we used a high-throughput approach in conjunction with flow cytometric sorting to identify and quantify HTLV-2-infected T-cell clones in 28 individuals with natural infection. We show that while genome-wide integration site preferences in vivo were similar to those found in HTLV-1 infection, expansion of HTLV-2-infected clones did not demonstrate the same significant association with the genomic environment of the integrated provirus. The proviral load in HTLV-2 is almost confined to CD8+ T-cells and is composed of a small number of often highly expanded clones. The HTLV-2 load correlated significantly with the degree of dispersion of the clone frequency distribution, which was highly stable over ∼8 years. These results suggest that there are significant differences in the selection forces that control the clonal expansion of virus-infected cells in HTLV-1 and HTLV-2 infection. In addition, our data demonstrate that strong virus-driven proliferation per se does not predispose to malignant transformation in oncoretroviral infections.

Retroviruses integrate into a shared, non-palindromic DNA motif.

Many DNA-binding factors, such as transcription factors, form oligomeric complexes with structural symmetry that bind to palindromic DNA sequences1. Palindromic consensus nucleotide sequences are also found at the genomic integration sites of retroviruses2–6 and other transposable elements7–9, and it has been suggested that this palindromic consensus arises as a consequence of the structural symmetry in the integrase complex2,3. However, we show here that the palindromic consensus sequence is not present in individual integration sites of human T-cell lymphotropic virus type 1 (HTLV-1) and human immunodeficiency virus type 1 (HIV-1), but arises in the population average as a consequence of the existence of a non-palindromic nucleotide motif that occurs in approximately equal proportions on the plus strand and the minus strand of the host genome. We develop a generally applicable algorithm to sort the individual integration site sequences into plus-strand and minus-strand subpopulations, and use this to identify the integration site nucleotide motifs of five retroviruses of different genera: HTLV-1, HIV-1, murine leukaemia virus (MLV), avian sarcoma leucosis virus (ASLV) and prototype foamy virus (PFV). The results reveal a non-palindromic motif that is shared between these retroviruses.

HTLV-1 proviral integration sites differ between asymptomatic carriers and patients with HAM/TSP

BackgroundHTLV-1 causes proliferation of clonal populations of infected T cells in vivo, each clone defined by a unique proviral integration site in the host genome. The proviral load is strongly correlated with odds of the inflammatory disease HTLV-1-associated myelopathy/tropical spastic paraparesis (HAM/TSP). There is evidence that asymptomatic HTLV-1 carriers (ACs) have a more effective CD8 + T cell response, including a higher frequency of HLA class I alleles able to present peptides from a regulatory protein of HTLV-1, HBZ. We have previously shown that specific features of the host genome flanking the proviral integration site favour clone survival and spontaneous expression of the viral transactivator protein Tax in naturally infected PBMCs ex vivo. However, the previous studies were not designed or powered to detect differences in integration site characteristics between ACs and HAM/TSP patients. Here, we tested the hypothesis that the genomic environment of the provirus differs systematically between ACs and HAM/TSP patients, and between individuals with strong or weak HBZ presentation.MethodsWe used our recently described high-throughput protocol to map and quantify integration sites in 95 HAM/TSP patients and 68 ACs from Kagoshima, Japan, and 75 ACs from Kumamoto, Japan. Individuals with 2 or more HLA class I alleles predicted to bind HBZ peptides were classified ‘strong’ HBZ binders; the remainder were classified ‘weak binders’.ResultsThe abundance of HTLV-1-infected T cell clones in vivo was correlated with proviral integration in genes and in areas with epigenetic marks associated with active regulatory elements. In clones of equivalent abundance, integration sites in genes and active regions were significantly more frequent in ACs than patients with HAM/TSP, irrespective of HBZ binding and proviral load. Integration sites in genes were also more frequent in strong HBZ binders than weak HBZ binders.ConclusionClonal abundance is correlated with integration in a transcriptionally active genomic region, and these regions may promote cell proliferation. A clone that reaches a given abundance in vivo is more likely to be integrated in a transcriptionally active region in individuals with a more effective anti-HTLV-1 immune response, such those who can present HBZ peptides or those who remain asymptomatic.