Detection of elusive DNA copy-number variations in hereditary disease and cancer through the use of noncoding and off-target sequencing reads.

Copy-number variants (CNVs) play a substantial role in the molecular pathogenesis of hereditary disease and cancer, as well as in normal human interindividual variation. However, they are still rather difficult to identify in mainstream sequencing projects, especially involving exome sequencing, because they often occur in DNA regions that are not targeted for analysis. To overcome this problem, we developed OFF-PEAK, a user-friendly CNV detection tool that builds on a denoising approach and the use of "off-target" DNA reads, which are usually discarded by sequencing pipelines. We benchmarked OFF-PEAK on data from targeted sequencing of 96 cancer samples, as well as 130 exomes of individuals with inherited retinal disease from three different populations. For both sets of data, OFF-PEAK demonstrated excellent performance (>95% sensitivity and >80% specificity vs. experimental validation) in detecting CNVs from in silico data alone, indicating its immediate applicability to molecular diagnosis and genetic research.

NK cells and monocytes modulate primary HTLV-1 infection.

We investigated the impact of monocytes, NK cells, and CD8+ T-cells in primary HTLV-1 infection by depleting cell subsets and exposing macaques to either HTLV-1 wild type (HTLV-1WT) or to the HTLV-1p12KO mutant unable to infect replete animals due to a single point mutation in orf-I that inhibits its expression. The orf-I encoded p8/p12 proteins counteract cytotoxic NK and CD8+ T-cells and favor viral DNA persistence in monocytes. Double NK and CD8+ T-cells or CD8 depletion alone accelerated seroconversion in all animals exposed to HTLV-1WT. In contrast, HTLV-1p12KO infectivity was fully restored only when NK cells were also depleted, demonstrating a critical role of NK cells in primary infection. Monocyte/macrophage depletion resulted in accelerated seroconversion in all animals exposed to HTLV-1WT, but antibody titers to the virus were low and not sustained. Seroconversion did not occur in most animals exposed to HTLV-1p12KO. In vitro experiments in human primary monocytes or THP-1 cells comparing HTLV-1WT and HTLV-1p12KO demonstrated that orf-I expression is associated with inhibition of inflammasome activation in primary cells, with increased CD47 "don't-eat-me" signal surface expression in virus infected cells and decreased monocyte engulfment of infected cells. Collectively, our data demonstrate a critical role for innate NK cells in primary infection and suggest a dual role of monocytes in primary infection. On one hand, orf-I expression increases the chances of viral transmission by sparing infected cells from efferocytosis, and on the other may protect the engulfed infected cells by modulating inflammasome activation. These data also suggest that, once infection is established, the stoichiometry of orf-I expression may contribute to the chronic inflammation observed in HTLV-1 infection by modulating monocyte efferocytosis.

PDZ domain-binding motif of Tax sustains T-cell proliferation in HTLV-1-infected humanized mice.

Human T-cell leukemia virus type 1 (HTLV-1) is the etiological agent of adult T-cell leukemia/lymphoma (ATLL), an aggressive malignant proliferation of activated CD4+ T lymphocytes. The viral Tax oncoprotein is critically involved in both HTLV-1-replication and T-cell proliferation, a prerequisite to the development of ATLL. In this study, we investigated the in vivo contribution of the Tax PDZ domain-binding motif (PBM) to the lymphoproliferative process. To that aim, we examined T-cell proliferation in humanized mice (hu-mice) carrying a human hemato-lymphoid system infected with either a wild type (WT) or a Tax PBM-deleted (ΔPBM) provirus. We observed that the frequency of CD4+ activated T-cells in the peripheral blood and in the spleen was significantly higher in WT than in ΔPBM hu-mice. Likewise, human T-cells collected from WT hu-mice and cultivated in vitro in presence of interleukin-2 were proliferating at a higher level than those from ΔPBM animals. We next examined the association of Tax with the Scribble PDZ protein, a prominent regulator of T-cell polarity, in human T-cells analyzed either after ex vivo isolation or after in vitro culture. We confirmed the interaction of Tax with Scribble only in T-cells from the WT hu-mice. This association correlated with the presence of both proteins in aggregates at the leading edge of the cells and with the formation of long actin filopods. Finally, data from a comparative genome-wide transcriptomic analysis suggested that the PBM-PDZ association is implicated in the expression of genes regulating proliferation, apoptosis and cytoskeletal organization. Collectively, our findings suggest that the Tax PBM is an auxiliary motif that contributes to the sustained growth of HTLV-1 infected T-cells in vivo and in vitro and is essential to T-cell immortalization.

Characterization of novel Bovine Leukemia Virus (BLV) antisense transcripts by deep sequencing reveals constitutive expression in tumors and transcriptional interaction with viral microRNAs.

BACKGROUND: Bovine Leukemia Virus (BLV) is a deltaretrovirus closely related to the Human T cell leukemia virus-1 (HTLV-1). Cattle are the natural host of BLV where it integrates into B-cells, producing a lifelong infection. Most infected animals remain asymptomatic but following a protracted latency period about 5 % develop an aggressive leukemia/lymphoma, mirroring the disease trajectory of HTLV-1. The mechanisms by which these viruses provoke cellular transformation remain opaque. In both viruses little or no transcription is observed from the 5'LTR in tumors, however the proviruses are not transcriptionally silent. In the case of BLV a cluster of RNA polymerase III transcribed microRNAs are highly expressed, while the HTLV-1 antisense transcript HBZ is consistently found in all tumors examined. RESULTS: Here, using RNA-seq, we demonstrate that the BLV provirus also constitutively expresses antisense transcripts in all leukemic and asymptomatic samples examined. The first transcript (AS1) can be alternately polyadenylated, generating a transcript of ~600 bp (AS1-S) and a less abundant transcript of ~2200 bp (AS1-L). Alternative splicing creates a second transcript of ~400 bp (AS2). The coding potential of AS1-S/L is ambiguous, with a small open reading frame of 264 bp, however the transcripts are primarily retained in the nucleus, hinting at a lncRNA-like role. The AS1-L transcript overlaps the BLV microRNAs and using high throughput sequencing of RNA-ligase-mediated (RLM) 5'RACE, we show that the RNA-induced silencing complex (RISC) cleaves AS1-L. Furthermore, experiments using altered BLV proviruses with the microRNAs either deleted or inverted point to additional transcriptional interference between the two viral RNA species. CONCLUSIONS: The identification of novel viral antisense transcripts shows the BLV provirus to be far from silent in tumors. Furthermore, the consistent expression of these transcripts in both leukemic and nonmalignant clones points to a vital role in the life cycle of the virus and its tumorigenic potential. Additionally, the cleavage of the AS1-L transcript by the BLV encoded microRNAs and the transcriptional interference between the two viral RNA species suggest a shared role in the regulation of BLV.

Full-Length Single-Cell RNA-Sequencing with FLASH-seq.

The single-cell RNA-sequencing (scRNA-seq) field has evolved tremendously since the first paper was published back in 2009 (Tang et al. Nat Methods 6:377-382, 2009). While the first methods analyzed just a handful of cells, the throughput and performance rapidly increased over a very short time span. However, it was not until the introduction of emulsion droplets methods, such as the well-known kits commercialized by 10x Genomics, that the robust and reproducible analysis of thousands of cells became feasible (Zheng et al Massively parallel digital transcriptional profiling of single cells. Nat Commun 8:14049, 2017). Despite generating data at a speed and a cost per cell that remains unmatched for full-length protocols like Smart-seq (Hagemann-Jensen et al Single-cell RNA counting at allele and isoform resolution using Smart-seq3. Nat Biotechnol 38:708-714, 2020; Picelli et al Smart-seq2 for sensitive full-length transcriptome profiling in single cells. Nat Methods 10:1096-1098, 2013), scRNA-seq in droplets still comes with the drawback of addressing only the terminal portion of the transcripts, thus lacking the required sensitivity for comprehensively analyzing the entire transcriptome.Building upon the existing Smart-seq2/3 workflows (Hagemann-Jensen et al Single-cell RNA counting at allele and isoform resolution using Smart-seq3. Nat Biotechnol 38:708-714, 2020; Picelli et al Smart-seq2 for sensitive full-length transcriptome profiling in single cells. Nat Methods 10:1096-1098, 2013), we developed FLASH-seq (FS), a new full-length scRNA-seq method capable of detecting a significantly higher number of genes than previous versions, requiring limited hands-on time and with a great potential for customization (Hahaut et al. Lightning Fast and Highly Sensitive Full-Length Single-cell sequencing using FLASH-Seq. http://biorxiv.org/lookup/doi/10.1101/2021.07.14.452217. https://doi.org/10.1101/2021.07.14.452217, 2021). Here, we present three variants of the FS protocol.Standard FLASH-seq (FS), which builds upon Smart-seq2 developed in the past, is non-stranded and does not use unique molecular identifiers (UMIs) but still remains the easiest method to measure gene expression in a cell population.FLASH-seq low-amplification (FS-LA) represents the fastest method, which generates sequencing-ready libraries in 4.5 h, without sacrificing performance.FLASH-seq with UMIs (FS-UMI) builds upon the same principle as Smart-seq3 and introduces UMIs for molecule counting and isoform reconstruction. The newly designed template-switching oligonucleotide (TSO) contains a 5-bp spacer, which allows the generation of high-quality data while minimizing the amount of strand-invasion artifacts.

Fast and highly sensitive full-length single-cell RNA sequencing using FLASH-seq.

We present FLASH-seq (FS), a full-length single-cell RNA sequencing (scRNA-seq) method with increased sensitivity and reduced hands-on time compared to Smart-seq3. The entire FS protocol can be performed in ~4.5 hours, is simple to automate and can be easily miniaturized to decrease resource consumption. The FS protocol can also use unique molecular identifiers (UMIs) for molecule counting while displaying reduced strand-invasion artifacts. FS will be especially useful for characterizing gene expression at high resolution across multiple samples.

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.

HTLV-1 infection promotes excessive T cell activation and transformation into adult T cell leukemia/lymphoma.

Human T cell leukemia virus type 1 (HTLV-1) mainly infects CD4+ T cells and induces chronic, persistent infection in infected individuals, with some developing adult T cell leukemia/lymphoma (ATL). HTLV-1 alters cellular differentiation, activation, and survival; however, it is unknown whether and how these changes contribute to the malignant transformation of infected cells. In this study, we used single-cell RNA-sequencing and T cell receptor-sequencing to investigate the differentiation and HTLV-1-mediated transformation of T cells. We analyzed 87,742 PBMCs from 12 infected and 3 uninfected individuals. Using multiple independent bioinformatics methods, we demonstrated the seamless transition of naive T cells into activated T cells, whereby HTLV-1-infected cells in an activated state further transformed into ATL cells, which are characterized as clonally expanded, highly activated T cells. Notably, the greater the activation state of ATL cells, the more they acquire Treg signatures. Intriguingly, the expression of HLA class II genes in HTLV-1-infected cells was uniquely induced by the viral protein Tax and further upregulated in ATL cells. Functional assays revealed that HTLV-1-infected cells upregulated HLA class II molecules and acted as tolerogenic antigen-presenting cells to induce anergy of antigen-specific T cells. In conclusion, our study revealed the in vivo mechanisms of HTLV-1-mediated transformation and immune escape at the single-cell level.

PCIP-seq: simultaneous sequencing of integrated viral genomes and their insertion sites with long reads.

The integration of a viral genome into the host genome has a major impact on the trajectory of the infected cell. Integration location and variation within the associated viral genome can influence both clonal expansion and persistence of infected cells. Methods based on short-read sequencing can identify viral insertion sites, but the sequence of the viral genomes within remains unobserved. We develop PCIP-seq, a method that leverages long reads to identify insertion sites and sequence their associated viral genome. We apply the technique to exogenous retroviruses HTLV-1, BLV, and HIV-1, endogenous retroviruses, and human papillomavirus.

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

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

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.