The APOBEC3B cytidine deaminase is an adenovirus restriction factor.

Human adenoviruses (HAdVs) are a large family of DNA viruses counting more than a hundred strains divided into seven species (A to G). HAdVs induce respiratory tract infections, gastroenteritis and conjunctivitis. APOBEC3B is a cytidine deaminase that restricts several DNA viruses. APOBEC3B is also implicated in numerous cancers where it is responsible for the introduction of clustered mutations into the cellular genome. In this study, we demonstrate that APOBEC3B is an adenovirus restriction factor acting through a deaminase-dependent mechanism. APOBEC3B introduces C-to-T clustered mutations into the adenovirus genome. APOBEC3B reduces the propagation of adenoviruses by limiting viral genome replication, progression to late phase, and production of infectious virions. APOBEC3B restriction efficiency varies between adenoviral strains, the A12 strain being more sensitive to APOBEC3B than the B3 or C2 strains. In A12-infected cells, APOBEC3B clusters in the viral replication centers. Importantly, we show that adenovirus infection leads to a reduction of the quantity and/or enzymatic activity of the APOBEC3B protein depending on the strains. The A12 strain seems less able to resist APOBEC3B than the B3 or C2 strains, a characteristic which could explain the strong depletion of the APOBEC3-targeted motifs in the A12 genome. These findings suggest that adenoviruses evolved different mechanisms to antagonize APOBEC3B. Elucidating these mechanisms could benefit the design of cancer treatments. This study also identifies adenoviruses as triggers of the APOBEC3B-mediated innate response. The involvement of certain adenoviral strains in the genesis of the APOBEC3 mutational signature observed in tumors deserves further study.

Infection of Bronchial Epithelial Cells by the Human Adenoviruses A12, B3, and C2 Differently Regulates the Innate Antiviral Effector APOBEC3B.

Human adenoviruses (HAdVs) are a large family of DNA viruses that include more than 100 genotypes divided into seven species (A to G) and induce respiratory tract infections, gastroenteritis, and conjunctivitis. Genetically modified adenoviruses are also used as vaccines, gene therapies, and anticancer treatments. The APOBEC3s are a family of cytidine deaminases that restrict viruses by introducing mutations in their genomes. Viruses developed different strategies to cope with the APOBEC3 selection pressure, but nothing is known on the interplay between the APOBEC3s and the HAdVs. In this study, we focused on three HAdV strains: the B3 and C2 strains, as they are very frequent, and the A12 strain, which is less common but is oncogenic in animal models. We demonstrated that the three HAdV strains induce a similar APOBEC3B upregulation at the transcriptional level. At the protein level, however, APOBEC3B is abundantly expressed during HAdV-A12 and -C2 infection and shows a nuclear distribution. On the contrary, APOBEC3B is barely detectable in HAdV-B3-infected cells. APOBEC3B deaminase activity is detected in total protein extracts upon HAdV-A12 and -C2 infection. Bioinformatic analysis demonstrates that the HAdV-A12 genome bears a stronger APOBEC3 evolutionary footprint than that of the HAdV-C2 and HAdV-B3 genomes. Our results show that HAdV infection triggers the transcriptional upregulation of the antiviral innate effector APOBEC3B. The discrepancies between the APOBEC3B mRNA and protein levels might reflect the ability of some HAdV strains to antagonize the APOBEC3B protein. These findings point toward an involvement of APOBEC3B in HAdV restriction and evolution. IMPORTANCE The APOBEC3 family of cytosine deaminases has important roles in antiviral innate immunity and cancer. Notably, APOBEC3A and APOBEC3B are actively upregulated by several DNA tumor viruses and contribute to transformation by introducing mutations in the cellular genome. Human adenoviruses (HAdVs) are a large family of DNA viruses that cause generally asymptomatic infections in immunocompetent adults. HAdVs encode several oncogenes, and some HAdV strains, like HAdV-A12, induce tumors in hamsters and mice. Here, we show that HAdV infection specifically promotes the expression of the APOBEC3B gene. We report that infection with the A12 strain induces a strong expression of an enzymatically active APOBEC3B protein in bronchial epithelial cells. We provide bioinformatic evidence that HAdVs' genomes and notably the A12 genome are under APOBEC3 selection pressure. Thus, APOBEC3B might contribute to adenoviral restriction, diversification, and oncogenic potential of particular strains.

Footprint of the host restriction factors APOBEC3 on the genome of human viruses.

APOBEC3 enzymes are innate immune effectors that introduce mutations into viral genomes. These enzymes are cytidine deaminases which transform cytosine into uracil. They preferentially mutate cytidine preceded by thymidine making the 5'TC motif their favored target. Viruses have evolved different strategies to evade APOBEC3 restriction. Certain viruses actively encode viral proteins antagonizing the APOBEC3s, others passively face the APOBEC3 selection pressure thanks to a depleted genome for APOBEC3-targeted motifs. Hence, the APOBEC3s left on the genome of certain viruses an evolutionary footprint. The aim of our study is the identification of these viruses having a genome shaped by the APOBEC3s. We analyzed the genome of 33,400 human viruses for the depletion of APOBEC3-favored motifs. We demonstrate that the APOBEC3 selection pressure impacts at least 22% of all currently annotated human viral species. The papillomaviridae and polyomaviridae are the most intensively footprinted families; evidencing a selection pressure acting genome-wide and on both strands. Members of the parvoviridae family are differentially targeted in term of both magnitude and localization of the footprint. Interestingly, a massive APOBEC3 footprint is present on both strands of the B19 erythroparvovirus; making this viral genome one of the most cleaned sequences for APOBEC3-favored motifs. We also identified the endemic coronaviridae as significantly footprinted. Interestingly, no such footprint has been detected on the zoonotic MERS-CoV, SARS-CoV-1 and SARS-CoV-2 coronaviruses. In addition to viruses that are footprinted genome-wide, certain viruses are footprinted only on very short sections of their genome. That is the case for the gamma-herpesviridae and adenoviridae where the footprint is localized on the lytic origins of replication. A mild footprint can also be detected on the negative strand of the reverse transcribing HIV-1, HIV-2, HTLV-1 and HBV viruses. Together, our data illustrate the extent of the APOBEC3 selection pressure on the human viruses and identify new putatively APOBEC3-targeted viruses.

Evolution of retrovirus-infected premalignant T-cell clones prior to adult T-cell leukemia/lymphoma diagnosis.

Adult T-cell leukemia/lymphoma (ATL) is an aggressive hematological malignancy caused by human T-cell leukemia virus type-1 (HTLV-1). ATL is preceded by decades of chronic HTLV-1 infection, and the tumors carry both somatic mutations and proviral DNA integrated into the tumor genome. In order to gain insight into the oncogenic process, we used targeted sequencing to track the evolution of the malignant clone in 6 individuals, 2 to 10 years before the diagnosis of ATL. Clones of premalignant HTLV-1-infected cells bearing known driver mutations were detected in the blood up to 10 years before individuals developed acute and lymphoma subtype ATL. Six months before diagnosis, the total number and variant allele fraction of mutations increased in the blood. Peripheral blood mononuclear cells from premalignant cases (1 year prediagnosis) had significantly higher mutational burden in genes frequently mutated in ATL than did high-risk, age-matched HTLV-1 carriers who remained ATL-free after a median of 10 years of follow-up. These data show that HTLV-1-infected T-cell clones carrying key oncogenic driver mutations can be detected in cases of ATL years before the onset of symptoms. Early detection of such mutations may enable earlier and more effective intervention to prevent the development of ATL.

Long-term clinical remission maintained after cessation of zidovudine and interferon-α therapy in chronic adult T-cell leukemia/lymphoma.

Globally, > 5-10 million people are estimated to be infected with Human T-lymphotropic virus type 1 (HTLV-1), of whom ~ 5% develop adult T-cell leukemia/lymphoma (ATL). Despite advances in chemotherapy, overall survival (OS) has not improved in the 35 years since HTLV-1 was first described. In Europe/USA, combination treatment with zidovudine and interferon-α (ZDV/IFN-α) has substantially changed the management of patients with the leukemic subtypes of ATL (acute or unfavorable chronic ATL) and is under clinical trial evaluation in Japan. However, there is only a single published report of long-term clinical remission on discontinuing ZDV/IFN-α therapy and the optimal duration of treatment is unknown. Anecdotal cases where therapy is discontinued due to side effects or compliance have been associated with rapid disease relapse, and it has been widely accepted that the majority of patients will require life-long therapy. The development of molecular methods to quantify minimal residual disease is essential to potentially guide therapy for individual patients. Here, for the first time, we report molecular evidence that supports long-term clinical remission in a patient who was previously treated with ZDV/IFN-α for 5 years, and who has now been off all therapy for over 6 years.

High-Throughput Mapping and Clonal Quantification of Retroviral Integration Sites.

We describe here a method to identify the position of retroviral insertion sites and simultaneously to quantify the absolute abundance of each clone, i.e., the number of cells having the provirus inserted at a given place in the host genome. The method is based on random shearing of the host cell DNA, followed by a linker-mediated PCR to amplify the genomic regions flanking the proviruses, and high-throughput sequencing of the amplicons. The quantification of the abundance of each infected clone allowed us to develop two new metrics: i. the oligoclonality index, which quantifies the nonuniformity of the distribution of clone abundance, and ii. an estimator of the total number of clones in the body of the host. These new tools are valuable for the study of retroviral infections and can also be adapted for the tracking of gene-edited cells.

Whole genome sequencing of 51 breast cancers reveals that tumors are devoid of bovine leukemia virus DNA.

Controversy exists regarding the association of bovine leukemia virus (BLV) and breast cancer. PCR-based experimental evidence indicates that BLV DNA is present in breast tissue and that as many as 37% of cancer cases may be attributable to viral exposure. Since this association might have major consequences for human health, we evaluated 51 whole genomes of breast cancer samples for the presence of BLV DNA. Among 32 billion sequencing reads retrieved from the NCBI database of genotype and phenotype, none mapped on different strains of the BLV genome. Controls for sequence divergence and proviral loads further validated the approach. This unbiased analysis thus excludes a clonal insertion of BLV in breast tumor cells and strongly argues against an association between BLV and breast cancer.

Bovine Leukemia Virus Small Noncoding RNAs Are Functional Elements That Regulate Replication and Contribute to Oncogenesis In Vivo.

Retroviruses are not expected to encode miRNAs because of the potential problem of self-cleavage of their genomic RNAs. This assumption has recently been challenged by experiments showing that bovine leukemia virus (BLV) encodes miRNAs from intragenomic Pol III promoters. The BLV miRNAs are abundantly expressed in B-cell tumors in the absence of significant levels of genomic and subgenomic viral RNAs. Using deep RNA sequencing and functional reporter assays, we show that miRNAs mediate the expression of genes involved in cell signaling, cancer and immunity. We further demonstrate that BLV miRNAs are essential to induce B-cell tumors in an experimental model and to promote efficient viral replication in the natural host.

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.

Massive depletion of bovine leukemia virus proviral clones located in genomic transcriptionally active sites during primary infection.

Deltaretroviruses such as human T-lymphotropic virus type 1 (HTLV-1) and bovine leukemia virus (BLV) induce a persistent infection that remains generally asymptomatic but can also lead to leukemia or lymphoma. These viruses replicate by infecting new lymphocytes (i.e. the infectious cycle) or via clonal expansion of the infected cells (mitotic cycle). The relative importance of these two cycles in viral replication varies during infection. The majority of infected clones are created early before the onset of an efficient immune response. Later on, the main replication route is mitotic expansion of pre-existing infected clones. Due to the paucity of available samples and for ethical reasons, only scarce data is available on early infection by HTLV-1. Therefore, we addressed this question in a comparative BLV model. We used high-throughput sequencing to map and quantify the insertion sites of the provirus in order to monitor the clonality of the BLV-infected cells population (i.e. the number of distinct clones and abundance of each clone). We found that BLV propagation shifts from cell neoinfection to clonal proliferation in about 2 months from inoculation. Initially, BLV proviral integration significantly favors transcribed regions of the genome. Negative selection then eliminates 97% of the clones detected at seroconversion and disfavors BLV-infected cells carrying a provirus located close to a promoter or a gene. Nevertheless, among the surviving proviruses, clone abundance positively correlates with proximity of the provirus to a transcribed region. Two opposite forces thus operate during primary infection and dictate the fate of long term clonal composition: (1) initial integration inside genes or promoters and (2) host negative selection disfavoring proviruses located next to transcribed regions. The result of this initial response will contribute to the proviral load set point value as clonal abundance will benefit from carrying a provirus in transcribed regions.

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.

Strongyloidiasis and infective dermatitis alter human T lymphotropic virus-1 clonality in vivo.

Human T-lymphotropic Virus-1 (HTLV-1) is a retrovirus that persists lifelong by driving clonal proliferation of infected T-cells. HTLV-1 causes a neuroinflammatory disease and adult T-cell leukemia/lymphoma. Strongyloidiasis, a gastrointestinal infection by the helminth Strongyloides stercoralis, and Infective Dermatitis associated with HTLV-1 (IDH), appear to be risk factors for the development of HTLV-1 related diseases. We used high-throughput sequencing to map and quantify the insertion sites of the provirus in order to monitor the clonality of the HTLV-1-infected T-cell population (i.e. the number of distinct clones and abundance of each clone). A newly developed biodiversity estimator called "DivE" was used to estimate the total number of clones in the blood. We found that the major determinant of proviral load in all subjects without leukemia/lymphoma was the total number of HTLV-1-infected clones. Nevertheless, the significantly higher proviral load in patients with strongyloidiasis or IDH was due to an increase in the mean clone abundance, not to an increase in the number of infected clones. These patients appear to be less capable of restricting clone abundance than those with HTLV-1 alone. In patients co-infected with Strongyloides there was an increased degree of oligoclonal expansion and a higher rate of turnover (i.e. appearance and disappearance) of HTLV-1-infected clones. In Strongyloides co-infected patients and those with IDH, proliferation of the most abundant HTLV-1⁺ T-cell clones is independent of the genomic environment of the provirus, in sharp contrast to patients with HTLV-1 infection alone. This implies that new selection forces are driving oligoclonal proliferation in Strongyloides co-infection and IDH. We conclude that strongyloidiasis and IDH increase the risk of development of HTLV-1-associated diseases by increasing the rate of infection of new clones and the abundance of existing HTLV-1⁺ 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.

HTLV-1 propels thymic human T cell development in "human immune system" Rag2⁻/⁻ gamma c⁻/⁻ mice.

Alteration of early haematopoietic development is thought to be responsible for the onset of immature leukemias and lymphomas. We have previously demonstrated that Tax(HTLV-1) interferes with ß-selection, an important checkpoint of early thymopoiesis, indicating that human T-cell leukemia virus type 1 (HTLV-1) infection has the potential to perturb thymic human αß T-cell development. To verify that inference and to clarify the impact of HTLV-1 infection on human T-cell development, we investigated the in vivo effects of HTLV-1 infection in a "Human Immune System" (HIS) Rag2⁻/⁻γ(c)⁻/⁻ mouse model. These mice were infected with HTLV-1, at a time when the three main subpopulations of human thymocytes have been detected. In all but two inoculated mice, the HTLV-1 provirus was found integrated in thymocytes; the proviral load increased with the length of the infection period. In the HTLV-1-infected mice we observed alterations in human T-cell development, the extent of which correlated with the proviral load. Thus, in the thymus of HTLV-1-infected HIS Rag2⁻/⁻γc⁻/⁻ mice, mature single-positive (SP) CD4⁺ and CD8⁺ cells were most numerous, at the expense of immature and double-positive (DP) thymocytes. These SP cells also accumulated in the spleen. Human lymphocytes from thymus and spleen were activated, as shown by the expression of CD25: this activation was correlated with the presence of tax mRNA and with increased expression of NF-kB dependent genes such as bfl-1, an anti-apoptotic gene, in thymocytes. Finally, hepato-splenomegaly, lymphadenopathy and lymphoma/thymoma, in which Tax was detected, were observed in HTLV-1-infected mice, several months after HTLV-1 infection. These results demonstrate the potential of the HIS Rag2⁻/⁻γ(c)⁻/⁻ animal model to elucidate the initial steps of the leukemogenic process induced by HTLV-1.

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.