GWAS reveals determinants of mobilization rate and dynamics of an active endogenous retrovirus of cattle.

Five to ten percent of mammalian genomes is occupied by multiple clades of endogenous retroviruses (ERVs), that may count thousands of members. New ERV clades arise by retroviral infection of the germline followed by expansion by reinfection and/or retrotransposition. ERV mobilization is a source of deleterious variation, driving the emergence of ERV silencing mechanisms, leaving "DNA fossils". Here we show that the ERVK[2-1-LTR] clade is still active in the bovine and a source of disease-causing alleles. We develop a method to measure the rate of ERVK[2-1-LTR] mobilization, finding an average of 1 per ~150 sperm cells, with >10-fold difference between animals. We perform a genome-wide association study and identify eight loci affecting ERVK[2-1-LTR] mobilization. We provide evidence that polymorphic ERVK[2-1-LTR] elements in four of these loci cause the association. We generate a catalogue of full length ERVK[2-1-LTR] elements, and show that it comprises 15% of C-type autonomous elements, and 85% of D-type non-autonomous elements lacking functional genes. We show that >25% of the variance of mobilization rate is determined by the number of C-type elements, yet that de novo insertions are dominated by D-type elements. We propose that D-type elements act as parasite-of-parasite gene drives that may contribute to the observed demise of ERV elements.

Transmission of SARS-CoV-2 After COVID-19 Screening and Mitigation Measures for Primary School Children Attending School in Liège, Belgium.

Importance: Recent data suggest a relatively low incidence of COVID-19 among children. The possible role that children attending primary school may play in the transmission of SARS-CoV-2 remains poorly understood. Objective: To gain a better understanding of the possible role of children in the transmission of SARS-CoV-2. Design, Setting, and Participants: This prospective cohort study was conducted from September 21 to December 31, 2020, in a primary school in Liège, Belgium, among a volunteer sample of 181 children, parents, and school employees. Exposures: Participants were tested for SARS-CoV-2 infection once a week for 15 weeks through throat washing, performed with 5 mL of saline and collected in a sterile tube after approximately 30 seconds of gargling. Quantitative reverse transcription-polymerase chain reaction was performed to detect SARS-CoV-2 infection. Main Outcomes and Measures: In case of test positivity, participants were asked to complete a questionnaire aimed at determining the timing of symptom onset and symptom duration. SARS-CoV-2 genetic sequencing was also performed. Confirmed cases were linked based on available information on known contacts and viral sequences. Results: A total of 181 individuals participated in this study, including 63 children (34 girls [54.0%]; mean [SD] age, 8.6 [1.9] years [range, 5-13 years]) and 118 adults (75 women [63.6%]; mean [SD] age, 42.5 [5.7] years [range, 30-59 years]). Forty-five individuals (24.9%) tested positive: 13 children (20.6%; 95% CI, 10.6%-30.6%) and 32 adults (27.1%; 95% CI, 19.1%-35.7%) (P = .34). Children were more often asymptomatic compared with adults (6 [46.2%; 95% CI, 19.1%-73.3%] vs 4 of 31 [12.9%; 95% CI, 1.3%-24.5%]; P = .04). The median duration of symptoms was shorter in children than in adults (0.00 days [IQR, 0.00-1.00 days] vs 15.00 days [IQR, 7.00-22.00 days]). A reconstruction of the outbreak revealed that most transmission events occurred between teachers and between children within the school. Of the observed household transmission events, most seemed to have originated from a child or teacher who acquired the infection at school. Conclusions and Relevance: Despite the implementation of several mitigation measures, the incidence of COVID-19 among children attending primary school in this study was comparable to that observed among teachers and parents. Transmission tree reconstruction suggests that most transmission events originated from within the school. Additional measures should be considered to reduce the transmission of SARS-CoV-2 at school, including intensified testing.

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.

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.

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.

A Recurrent Mutation at Position 26340 of SARS-CoV-2 Is Associated with Failure of the E Gene Quantitative Reverse Transcription-PCR Utilized in a Commercial Dual-Target Diagnostic Assay.

Control of the ongoing severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic requires accurate laboratory testing to identify infected individuals while also clearing essential staff to continue to work. At the current time, a number of quantitative real-time PCR (qRT-PCR) assays have been developed to identify SARS-CoV-2, targeting multiple positions in the viral genome. While the mutation rate of SARS-CoV-2 is moderate, given the large number of transmission chains, it is prudent to monitor circulating viruses for variants that might compromise these assays. Here, we report the identification of a C-to-U transition at position 26340 of the SARS-CoV-2 genome that is associated with failure of the cobas SARS-CoV-2 E gene qRT-PCR in eight patients. As the cobas SARS-CoV-2 assay targets two positions in the genome, the individuals carrying this variant were still called SARS-CoV-2 positive. Whole-genome sequencing of SARS-CoV-2 showed all to carry closely related viruses. Examination of viral genomes deposited on GISAID showed this mutation has arisen independently at least four times. This work highlights the necessity of monitoring SARS-CoV-2 for the emergence of single-nucleotide polymorphisms that might adversely affect RT-PCRs used in diagnostics. Additionally, it argues that two regions in SARS-CoV-2 should be targeted to avoid false negatives.

Essential Role of Human T Cell Leukemia Virus Type 1 orf-I in Lethal Proliferation of CD4+ Cells in Humanized Mice.

Human T cell leukemia virus type 1 (HTLV-1) is the ethological agent of adult T cell leukemia/lymphoma (ATLL) and a number of lymphocyte-mediated inflammatory conditions, including HTLV-1-associated myelopathy/tropical spastic paraparesis. HTLV-1 orf-I encodes two proteins, p8 and p12, whose functions in humans are to counteract innate and adaptive responses and to support viral transmission. However, the in vivo requirements for orf-I expression vary in different animal models. In macaques, the ablation of orf-I expression by mutation of its ATG initiation codon abolishes the infectivity of the molecular clone HTLV-1p12KO In rabbits, HTLV-1p12KO is infective and persists efficiently. We used humanized mouse models to assess the infectivity of both wild-type HTLV-1 (HTLV-1WT) and HTLV-1p12KO We found that NOD/SCID/γC-/- c-kit+ mice engrafted with human tissues 1 day after birth (designated NSG-1d mice) were highly susceptible to infection by HTLV-1WT, with a syndrome characterized by the rapid polyclonal proliferation and infiltration of CD4+ CD25+ T cells into vital organs, weight loss, and death. HTLV-1 clonality studies revealed the presence of multiple clones of low abundance, confirming the polyclonal expansion of HTLV-1-infected cells in vivo HTLV-1p12KO infection in a bone marrow-liver-thymus (BLT) mouse model prone to graft-versus-host disease occurred only following reversion of the orf-I initiation codon mutation within weeks after exposure and was associated with high levels of HTLV-1 DNA in blood and the expansion of CD4+ CD25+ T cells. Thus, the incomplete reconstitution of the human immune system in BLT mice may provide a window of opportunity for HTLV-1 replication and the selection of viral variants with greater fitness.IMPORTANCE Humanized mice constitute a useful model for studying the HTLV-1-associated polyclonal proliferation of CD4+ T cells and viral integration sites in the human genome. The rapid death of infected animals, however, appears to preclude the clonal selection typically observed in human ATLL, which normally develops in 2 to 5% of individuals infected with HTLV-1. Nevertheless, the expansion of multiple clones of low abundance in these humanized mice mirrors the early phase of HTLV-1 infection in humans, providing a useful model to investigate approaches to inhibit virus-induced CD4+ T cell proliferation.

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.

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.

Deep sequencing reveals abundant noncanonical retroviral microRNAs in B-cell leukemia/lymphoma.

Viral tumor models have significantly contributed to our understanding of oncogenic mechanisms. How transforming delta-retroviruses induce malignancy, however, remains poorly understood, especially as viral mRNA/protein are tightly silenced in tumors. Here, using deep sequencing of broad windows of small RNA sizes in the bovine leukemia virus ovine model of leukemia/lymphoma, we provide in vivo evidence of the production of noncanonical RNA polymerase III (Pol III)-transcribed viral microRNAs in leukemic B cells in the complete absence of Pol II 5'-LTR-driven transcriptional activity. Processed from a cluster of five independent self-sufficient transcriptional units located in a proviral region dispensable for in vivo infectivity, bovine leukemia virus microRNAs represent ∼40% of all microRNAs in both experimental and natural malignancy. They are subject to strong purifying selection and associate with Argonautes, consistent with a critical function in silencing of important cellular and/or viral targets. Bovine leukemia virus microRNAs are strongly expressed in preleukemic and malignant cells in which structural and regulatory gene expression is repressed, suggesting a key role in tumor onset and progression. Understanding how Pol III-dependent microRNAs subvert cellular and viral pathways will contribute to deciphering the intricate perturbations that underlie malignant transformation.