The Ancient Evolutionary History of Polyomaviruses

Progressive multifocal leukoencephalopathy and the spectrum of JC virus-related disease

Progressive multifocal leukoencephalopathy (PML) is a devastating CNS infection caused by JC virus (JCV), a polyomavirus that commonly establishes persistent, asymptomatic infection in the general population. Emerging evidence that PML can be ameliorated with novel immunotherapeutic approaches calls for reassessment of PML pathophysiology and clinical course. PML results from JCV reactivation in the setting of impaired cellular immunity, and no antiviral therapies are available, so survival depends on reversal of the underlying immunosuppression. Antiretroviral therapies greatly reduce the risk of HIV-related PML, but many modern treatments for cancers, organ transplantation and chronic inflammatory disease cause immunosuppression that can be difficult to reverse. These treatments — most notably natalizumab for multiple sclerosis — have led to a surge of iatrogenic PML. The spectrum of presentations of JCV-related disease has evolved over time and may challenge current diagnostic criteria. Immunotherapeutic interventions, such as use of checkpoint inhibitors and adoptive T cell transfer, have shown promise but caution is needed in the management of immune reconstitution inflammatory syndrome, an exuberant immune response that can contribute to morbidity and death. Many people who survive PML are left with neurological sequelae and some with persistent, low-level viral replication in the CNS. As the number of people who survive PML increases, this lack of viral clearance could create challenges in the subsequent management of some underlying diseases.

In this Review, Cortese et al. provide an overview of the pathobiology and evolving presentations of progressive multifocal leukoencephalopathy and other diseases caused by JC virus, and discuss emerging immunotherapeutic approaches that could increase survival.

Progressive multifocal leukoencephalopathy (PML) is a rare, debilitating and often fatal disease of the CNS caused by JC virus (JCV).

JCV establishes asymptomatic, lifelong persistent or latent infection in immune competent hosts, but impairment of cellular immunity can lead to reactivation of JCV and PML.

PML most commonly occurs in patients with HIV infection or lymphoproliferative disease and in patients who are receiving natalizumab for treatment of multiple sclerosis.

The clinical phenotype of PML varies and is shaped primarily by the host immune response; changes in the treatment of underlying diseases associated with PML have changed phenotypes over time.

Other clinical manifestations of JCV infection have been described, including granule cell neuronopathy.

Survival of PML depends on reversal of the underlying immunosuppression; emerging immunotherapeutic strategies include use of checkpoint inhibitors and adoptive T cell transfer.

Prevalence and viral loads of polyomaviruses BKPyV, JCPyV, MCPyV, TSPyV and NJPyV and hepatitis viruses HBV, HCV and HEV in HIV-infected patients in China

Human polyomaviruses (PyVs) and hepatitis viruses are often more prevalent or persistent in human immunodeficiency virus (HIV)-infected persons and the associated diseases are more abundant than in immunocompetent individuals. Here, we evaluated seroreactivities and viral loads of human PyVs and hepatitis viruses in HIV/AIDS patients and the general population in China in the combination antiretroviral therapy (cART) era. A total of 810 HIV-1-infected patients and age- and sex-matched HIV-negative individuals were enrolled to assess seroprevalence of PyVs BKPyV, JCPyV, MCPyV, TSPyV, and NJPyV and hepatitis viruses HBV, HCV, and HEV. 583 (72%) patients received cART, and among them, 31.2% had undetectable HIV RNA. While no significant difference was observed in prevalence of anti-PyV antibodies between HIV-positive and -negative groups, serum DNA positivity and DNA copy level of MCPyV were higher in the HIV-positive group. Among HIV-infected patients, BKPyV DNA positivity was significantly higher in patients with CD4 + cell counts < 200 cells/mm³ compared to those with CD4 + cell counts > 500 cells/mm³, suggesting possible reactivation caused by HIV-induced immune suppression. Higher HBV and HCV seropositivities but not HEV seropositivity were also observed in the HIV-positive group. Further correlation analyses demonstrated that HBV and HEV are potential risk factors for increased prevalence of PyV infection.

You Will Never Walk Alone: Codispersal of JC Polyomavirus with Human Populations

JC polyomavirus (JCPyV) is one of the most prevalent human viruses. Findings based on the geographic distribution of viral subtypes suggested that JCPyV codiverged with human populations. This view was however challenged by data reporting a much more recent origin and expansion of JCPyV. We collected information on ∼1,100 worldwide strains and we show that their geographic distribution roughly corresponds to major human migratory routes. Bayesian phylogeographic analysis inferred a Subsaharan origin for JCPyV, although with low posterior probability. High confidence inference at internal nodes provided strong support for a long-standing association between the virus and human populations. In line with these data, pairwise FST values for JCPyV and human mtDNA sampled from the same areas showed a positive and significant correlation. Likewise, very strong relationships were found when node ages in the JCPyV phylogeny were correlated with human population genetic distances (nuclear-marker based FST). Reconciliation analysis detected a significant cophylogenetic signal for the human population and JCPyV trees. Notably, JCPyV also traced some relatively recent migration events such as the expansion of people from the Philippines/Taiwan area into Remote Oceania, the gene flow between North-Eastern Siberian and Ainus, and the Koryak contribution to Circum-Arctic Americans. Finally, different molecular dating approaches dated the origin of JCPyV in a time frame that precedes human out-of-Africa migration. Thus, JCPyV infected early human populations and accompanied our species during worldwide dispersal. JCPyV typing can provide reliable geographic information and the virus most likely adapted to the genetic background of human populations.

JCPyV VP1 Mutations in Progressive MultifocalLeukoencephalopathy: Altering Tropismor Mediating Immune Evasion?

Polyomaviruses are ubiquitous human pathogens that cause lifelong, asymptomatic infections in healthy individuals. Although these viruses are restrained by an intact immune system, immunocompromised individuals are at risk for developing severe diseases driven by resurgent viral replication. In particular, loss of immune control over JC polyomavirus can lead to the development of the demyelinating brain disease progressive multifocal leukoencephalopathy (PML). Viral isolates from PML patients frequently carry point mutations in the major capsid protein, VP1, which mediates virion binding to cellular glycan receptors. Because polyomaviruses are non-enveloped, VP1 is also the target of the host's neutralizing antibody response. Thus, VP1 mutations could affect tropism and/or recognition by polyomavirus-specific antibodies. How these mutations predispose susceptible individuals to PML and other JCPyV-associated CNS diseases remains to be fully elucidated. Here, we review the current understanding of polyomavirus capsid mutations and their effects on viral tropism, immune evasion, and virulence.

Multiple DNA viruses identified in multimammate mouse (Mastomys natalensis) populations from across regions of sub-Saharan Africa

The multimammate mouse (Mastomys natalensis; M. natalensis) serves as the main reservoir for the zoonotic arenavirus Lassa virus (LASV), and this has led to considerable investigation into the distribution of LASV and other related arenaviruses in this host species. In contrast to the situation with arenaviruses, the presence of other viruses in M. natalensis remains largely unexplored. In this study, herpesviruses and polyomaviruses were identified and partially characterized by PCR methods, sequencing, and phylogenetic analysis. In tissues sampled from M. natalensis populations in Côte d'Ivoire and Mali, six new DNA viruses (four betaherpesviruses, one gammaherpesvirus and one polyomavirus) were identified. Phylogenetic analysis based on glycoprotein B amino acid sequences showed that the herpesviruses clustered with cytomegaloviruses and rhadinoviruses of multiple rodent species. The complete circular genome of the newly identified polyomavirus was amplified by PCR. Amino acid sequence analysis of the large T antigen or VP1 showed that this virus clustered with a known polyomavirus from a house mouse (species Mus musculus polyomavirus 1). These two polyomaviruses form a clade with other rodent polyomaviruses, and the newly identified virus represents the third known polyomavirus of M. natalensis. This study represents the first identification of herpesviruses and the discovery of a novel polyomavirus in M. natalensis. In contrast to arenaviruses, we anticipate that these newly identified viruses represent a low zoonotic risk due to the normally highly restricted specificity of members of these two DNA virus families to their individual mammalian host species.

Structural Basis and Evolution of Glycan Receptor Specificities within the Polyomavirus Family

Virus attachment to cell surface receptors is critical for productive infection. In this study, we have used a structure-based approach to investigate the cell surface recognition event for New Jersey polyomavirus (NJPyV) and human polyomavirus 12 (HPyV12). These viruses belong to the polyomavirus family, whose members target different tissues and hosts, including mammals, birds, fish, and invertebrates. Polyomaviruses are nonenveloped viruses, and the receptor-binding site is located in their capsid protein VP1. The NJPyV capsid features a novel sialic acid-binding site that is shifted in comparison to other structurally characterized polyomaviruses but shared with a closely related simian virus. In contrast, HPyV12 VP1 engages terminal sialic acids in a manner similar to the human Trichodysplasia spinulosa-associated polyomavirus. Our structure-based phylogenetic analysis highlights that even distantly related avian polyomaviruses possess the same exposed sialic acid-binding site. These findings complement phylogenetic models of host-virus codivergence and may also reflect past host-switching events.

Asymptomatic infections with polyomaviruses in humans are common, but these small viruses can cause severe diseases in immunocompromised hosts. New Jersey polyomavirus (NJPyV) was identified via a muscle biopsy in an organ transplant recipient with systemic vasculitis, myositis, and retinal blindness, and human polyomavirus 12 (HPyV12) was detected in human liver tissue. The evolutionary origins and potential diseases are not well understood for either virus. In order to define their receptor engagement strategies, we first used nuclear magnetic resonance (NMR) spectroscopy to establish that the major capsid proteins (VP1) of both viruses bind to sialic acid in solution. We then solved crystal structures of NJPyV and HPyV12 VP1 alone and in complex with sialylated glycans. NJPyV employs a novel binding site for a α2,3-linked sialic acid, whereas HPyV12 engages terminal α2,3- or α2,6-linked sialic acids in an exposed site similar to that found in Trichodysplasia spinulosa-associated polyomavirus (TSPyV). Gangliosides or glycoproteins, featuring in mammals usually terminal sialic acids, are therefore receptor candidates for both viruses. Structural analyses show that the sialic acid-binding site of NJPyV is conserved in chimpanzee polyomavirus (ChPyV) and that the sialic acid-binding site of HPyV12 is widely used across the entire polyomavirus family, including mammalian and avian polyomaviruses. A comparison with other polyomavirus-receptor complex structures shows that their capsids have evolved to generate several physically distinct virus-specific receptor-binding sites that can all specifically engage sialylated glycans through a limited number of contacts. Small changes in each site may have enabled host-switching events during the evolution of polyomaviruses.

Bat-borne polyomaviruses in Europe reveal an evolutionary history of intrahost divergence with horseshoe bats distributed across the African and Eurasian continents

Polyomaviruses (PyVs) are small, circular dsDNA viruses carried by diverse vertebrates, including bats. Although previous studies have reported several horseshoe bat PyVs collected in Zambia and China, it is still unclear how PyVs evolved in this group of widely dispersed mammals. Horseshoe bats (genus Rhinolophus) are distributed across the Old World and are natural reservoirs of numerous pathogenic viruses. Herein, non-invasive bat samples from European horseshoe bat species were collected in Hungary for PyV identification and novel PyVs with complete genomes were successfully recovered from two different European horseshoe bat species. Genomic and phylogenetic analysis of the Hungarian horseshoe bat PyVs supported their classification into the genera Alphapolyomavirus and Betapolyomavirus. Notably, despite the significant geographical distances between the corresponding sampling locations, Hungarian PyVs exhibited high genetic relatedness with previously described Zambian and Chinese horseshoe bat PyVs, and phylogenetically clustered with these viruses in each PyV genus. Correlation and virus-host relationship analysis suggested that these PyVs co-diverged with their European, African and Asian horseshoe bat hosts distributed on different continents during their evolutionary history. Additionally, assessment of selective pressures over the major capsid protein (VP1) of horseshoe bat PyVs showed sites under positive selection located in motifs exposed to the exterior of the capsid. In summary, our findings revealed a pattern of stable intrahost divergence of horseshoe bat PyVs with their mammalian hosts on the African and Eurasian continents over evolutionary time.

JC Polyomavirus, progressive multifocal leukoencephalopathy and immune reconstitution inflammatory syndrome: a review

The human neurotropic virus JC Polyomavirus, a member of the Polyomaviridae family, is the opportunistic infectious agent causing progressive multifocal leukoencephalopathy, typically in immunocompromised individuals. The spectrum of underlying reasons for the systemic immunosuppression that permits JCV infection in the central nervous system has evolved over the past 2 decades, and therapeutic immunosuppression arousing JCV infection in the brain has become increasingly prominent as a trigger for PML. Effective immune restoration subsequent to human immunodeficiency virus-related suppression is now recognized as a cause for unexpected deterioration of symptoms in patients with PML, secondary to a rebound inflammatory phenomenon called immune reconstitution inflammatory syndrome, resulting in significantly increased morbidity and mortality in a disease already infamous for its lethality. This review addresses current knowledge regarding JC Polyomavirus, progressive multifocal leukoencephalopathy, progressive multifocal leukoencephalopathy-related immune reconstitution inflammatory syndrome, and the immunocompromised states that incite JC Polyomavirus central nervous system infection, and discusses prospects for the future management of these conditions.

The characterization of two novel neotropical primate papillomaviruses supports the ancient within-species diversity model

Papillomaviruses (PVs) are non-enveloped icosahedral viruses with a circular double-stranded DNA genome of ∼8,000 base pairs (bp). More than 200 different PV types have been identified to date in humans, which are distributed in five genera, with several strains associated with cancer development. Although widely distributed in vertebrates, Neotropical Primates (NP) PV infection was described for the first time only in 2016. Currently, four complete genomes of NP PVs have been characterized, three from Saimiri sciureus (SscPV1 to SscPV3) and one from Alouatta guariba (AgPV1). In this work, we describe two novel PV strains infecting Callithrix penicillata (provisionally named CpenPV1 and CpenPV2), using anal swab samples from animals residing at the Brasilia Primatology Center and next generation sequencing. The genomes of CpenPV1 (7,288 bp; 41.5% guanine-cytosine content - GC) and CpenPV2 (7,250 bp; 40.7% GC) contain the characteristic open reading frames (ORFs) for the early (E6, E7, E1, E2, and E4) and late (L2 and L1) PV genes. The L1 ORFs, commonly used for phylogenetic identification, share 76 per cent similarity with each other and differ 32 per cent from any other known PV, indicating that these new strains meet the criteria for defining novel species. PV genes phylogenetic variance was analyzed and different degrees of saturation revealed similar levels of topological heterogeneity, ruling out saturation as primary etiological factor for this phenomenon. Interestingly, the two CpenPV strains form a monophyletic clade within the Gammapapillomavirus genus (provisionally named gammapapillomavirus 32). Unlike for other NP PV strains, which grouped into a new sister genus of Alphapapillomavirus, this is the first report of NP PV strains grouping into a genus previously considered to exclusively comprise Old World Primates (OWP) PVs, including human PVs. These findings confirm the existence of a common ancestor for Gammapapillomavirus already infecting primates before the split of OWP and NP at ∼40 million years ago. Finally, our findings are consistent with an ancient within-species diversity model and emphasize the importance of increasing sampling to help understanding the PV-primate codivergence dynamics and pathogenic potential.

Molecular Mechanisms of Merkel Cell Polyomavirus Transformation and Replication

Viral infection underlies a significant share of the global cancer burden. Merkel cell polyomavirus (MCPyV) is the newest member of the human oncogenic virus family. Its discovery over a decade ago marked the beginning of an exciting era in human tumor virology. Since then, significant evidence has emerged to support the etiologic role of MCPyV in Merkel cell carcinoma (MCC), an extremely lethal form of skin cancer. MCPyV infection is widespread in the general population. MCC diagnoses have tripled over the past 20 years, but effective treatments are currently lacking. In this review, we highlight recent discoveries that have shaped our understanding of MCPyV oncogenic mechanism and host cellular tropism, as well as the molecular events occurring in the viral infectious life cycle. These insights will guide future efforts in developing novel virus-targeted therapeutic strategies for treating the devastating human cancers associated with this new tumorigenic virus.

Extensive Genetic Diversity of Polyomaviruses in Sympatric Bat Communities: Host Switching versus Coevolution

Polyomaviruses (PyVs) are small DNA viruses carried by diverse vertebrates. The evolutionary relationships of viruses and hosts remain largely unclear due to very limited surveillance in sympatric communities. In order to investigate whether PyVs can transmit among different mammalian species and to identify host-switching events in the field, we conducted a systematic study of a large collection of bats (n = 1,083) from 29 sympatric communities across China which contained multiple species with frequent contact. PyVs were detected in 21 bat communities, with 192 PyVs identified in 186 bats from 15 species within 6 families representing at least 28 newly described PyVs. Surveillance results and phylogenetic analyses surprisingly revealed three interfamily PyV host-switching events in these sympatric bat communities: two distinct PyVs were identified in two bat species in restricted geographical locations, while another PyV clustered phylogenetically with PyVs carried by bats from a different host family. Virus-host relationships of all discovered PyVs were also evaluated, and no additional host-switching events were found. PyVs were identified in different horseshoe bat species in sympatric communities without observation of host-switching events, showed high genomic identities, and clustered with each other. This suggested that even for PyVs with high genomic identities in closely related host species, the potential for host switching is low. In summary, our findings revealed that PyV host switching in sympatric bat communities can occur but is limited and that host switching of bat-borne PyVs is relatively rare on the predominantly evolutionary background of codivergence with their hosts.IMPORTANCE Since the discovery of murine polyomavirus in the 1950s, polyomaviruses (PyVs) have been considered highly host restricted in mammals. Sympatric bat communities commonly contain several different bat species in an ecological niche facilitating viral transmission, and they therefore represent a model to identify host-switching events of PyVs. In this study, we screened PyVs in a large number of bats in sympatric communities from diverse habitats across China. We provide evidence that cross-species bat-borne PyV transmission exists, though is limited, and that host-switching events appear relatively rare during the evolutionary history of these viruses. PyVs with close genomic identities were also identified in different bat species without host-switching events. Based on these findings, we propose an evolutionary scheme for bat-borne PyVs in which limited host-switching events occur on the background of codivergence and lineage duplication, generating the viral genetic diversity in bats.

Intra-patient viral evolution in polyomavirus-related diseases

Human polyomaviruses show relatively little genetic polymorphism between isolates, indicating that these viruses are genetically stable between hosts. However, it has become increasingly clear that intra-host molecular evolution is a feature of some polyomavirus (PyV) infections in humans. Mutations inducing premature stop codons in the early region of the integrated Merkel cell PyV genome lead to the expression of a truncated form of the large tumour (LT) antigen that is critical for the transformation of Merkel cell carcinoma (MCC) cells. Non-coding control region (NCCR) rearrangements and point mutations in virion protein (VP) 1 have been described in both JCPyV and BKPyV infections. In the context of JCPyV infection, molecular evolution at both these loci allows the virus to replicate effectively in the central nervous system, thereby leading to the development of progressive multifocal leukoencephalopathy (PML). In BKPyV infection, NCCR rearrangements have been linked to higher rates of virus replication in the kidney, and are proposed to play a direct causal role in the development of PyV-associated nephropathy. In all three of these infections, therefore, intra-host viral evolution appears to be an essential component of the disease process. This article is part of the theme issue 'Silent cancer agents: multi-disciplinary modelling of human DNA oncoviruses'.

Coinfections of Novel Polyomavirus, Anelloviruses and a Recombinant Strain of Myxoma Virus-MYXV-Tol Identified in Iberian Hares

Viruses are ubiquitous in nature; however, very few have been identified in the Leporid species. In the fall of 2018, an outbreak of myxomatosis in Iberian hares (Lepus granatensis) was reported in Spain and a novel recombinant myxoma virus strain (MYXV-Tol) was identified. To investigate variability within the recombinant region of the MYXV-Tol and identify any potential viral coinfections, samples (ear, eyelid or vaginal) of Iberian hares were collected from Spain and analyzed. The presence of the recombinant region of the MYXV-Tol was confirmed in six out of eleven samples analyzed. Additionally, a polyomavirus (family Polyomaviridae), representing a putative new species, and anelloviruses (family Anelloviridae) belonging to two putative species were identified, some as coinfection with the recombinant MYXV-Tol. The two polyomavirus genomes were identified in two hares and share >99% genome-wide identity. Based on the analysis of their large T-antigen, the new polyomavirus clusters in a distant clade from other mammals sharing <64% amino acid identity. A total of 14 anelloviruses were identified, which share 63-99% genome-wide identity. Overall, our results show a coinfection of different DNA viruses in the studied samples and raise awareness regarding the extensive unsampled diversity of viruses in hares.

Discovery of several thousand highly diverse circular DNA viruses

Although millions of distinct virus species likely exist, only approximately 9000 are catalogued in GenBank's RefSeq database. We selectively enriched for the genomes of circular DNA viruses in over 70 animal samples, ranging from nematodes to human tissue specimens. A bioinformatics pipeline, Cenote-Taker, was developed to automatically annotate over 2500 complete genomes in a GenBank-compliant format. The new genomes belong to dozens of established and emerging viral families. Some appear to be the result of previously undescribed recombination events between ssDNA and ssRNA viruses. In addition, hundreds of circular DNA elements that do not encode any discernable similarities to previously characterized sequences were identified. To characterize these ‘dark matter’ sequences, we used an artificial neural network to identify candidate viral capsid proteins, several of which formed virus-like particles when expressed in culture. These data further the understanding of viral sequence diversity and allow for high throughput documentation of the virosphere.

When scientists hunt for new DNA sequences, sometimes they get a lot more than they bargained for. Such is the case in metagenomic surveys, which analyze not just DNA of a particular organism, but all the DNA in an environment at large. A vexing problem with these surveys is the overwhelming number of DNA sequences detected that are so different from any known microbe that they cannot be classified using traditional approaches. However, some of these “known unknowns” are undoubtedly viral sequences, because only a fraction of the enormous diversity of viruses has been characterized.

This “viral dark matter” is a major obstacle for those studying viruses. This led Tisza et al. to attempt to classify some of the unknown viral sequences in their metagenomic surveys. The search, which specifically focused on viruses with circular DNA genomes, detected over 2,500 circular viral genomes. Intensive analysis revealed that many of these genomes had similar makeup to previously discovered viruses, but hundreds of them were totally different from any known virus, based on typical methods of comparison.

Computational analysis of genes that were conserved among some of these brand-new circular sequences often revealed virus-like features. Experiments on a few of these genes showed that they encoded proteins capable of forming particles reminiscent of characteristic viral shells, implying that these new sequences are indeed viruses.

Tisza et al. have added the 2,500 newly characterized viral sequences to the publicly accessible GenBank database, and the sequences are being considered for the more authoritative RefSeq database, which currently contains around 9,000 complete viral genomes. The expanded databases will hopefully now better equip scientists to explore the enormous diversity of viruses and help medics and veterinarians to detect disease-causing viruses in humans and other animals.

DETECTing Merkel Cell Polyomavirus in Merkel Tumors

Merkel cell carcinoma (MCC) is a rare, aggressive skin cancer caused either by Merkel cell polyomavirus (MCV) T antigen expression, post-integration (~80% cases), or by UV-mediated DNA damage. Interestingly, overall survival of MCV-positive Merkel cell carcinoma patients is better, making this differential information of significant diagnostic and prognostic value. Also, MCV provides a direct target for therapy in MCC patients. Currently, the methods used for diagnosis of MCV in tumors are often discordant and unreliable. Here we used a guided molecular scissors based-DNA Endonuclease Targeted CRISPR Trans Reporter (DETECTR) technique to develop an in vitro molecular diagnostic tool for MCV-positive MCC. DETECTR couples recombinase polymerase based amplification of target MCV DNA with Cas12a mediated detection. CRISPR diagnostics couple specific detection followed by cutting of the pathogenic DNA by the Cas enzyme-gRNA complex, with non-specific cutting of ssDNA that provides a measurable visual cue. To detect MCV DNA in MCC, we designed Cas12a gRNAs targeting the MCV DNA and tested their targeting efficiency, and sensitivity using a fluorophore quencher labeled reporter assay. We show that MCV DETECTR system can detect MCV integrated in Merkel tumor rapidly, specifically and with femto-molar sensitivity. Our study is a preliminary, proof-of-principle analysis showing the use of CRISPR for MCV diagnosis. Further validation in human tumor samples is needed for its clinical use in the near future. This new system is promising and we hope it can be coupled with immunohistochemical studies to diagnose the viral status of MCC in clinics soon.

Viral miRNAs as Active Players and Participants in Tumorigenesis

The theory that viruses play a role in human cancers is now supported by scientific evidence. In fact, around 12% of human cancers, a leading cause of morbidity and mortality in some regions, are attributed to viral infections. However, the molecular mechanism remains complex to decipher. In recent decades, the uncovering of cellular miRNAs, with their invaluable potential as diagnostic and prognostic biomarkers, has increased the number of studies being conducted regarding human cancer diagnosis. Viruses develop clever mechanisms to succeed in the maintenance of the viral life cycle, and some viruses, especially herpesviruses, encode for miRNA, v-miRNAs. Through this viral miRNA, the viruses are able to manipulate cellular and viral gene expression, driving carcinogenesis and escaping the host innate or adaptive immune system. In this review, we have discussed the main viral miRNAs and virally influenced cellular pathways, and their capability to drive carcinogenesis.

Human polyomavirus modulation of the host DNA damage response

The human DNA damage response (DDR) is a complex signaling network constituting many factors responsible for the preservation of genomic integrity. Human polyomaviruses (HPyVs) are able to harness the DDR machinery during their infectious cycle by expressing an array of tumor (T) antigens. These molecular interactions between human polyomavirus T antigens and the DDR create conditions that promote viral replication at the expense of host genomic stability to cause disease as well as carcinogenesis in the cases of the Merkel cell polyomavirus and BK polyomavirus. This review focuses on the six HPyVs with disease association, emphasizing strain-dependent differences in their selective manipulation of the DDR. Appreciation of the HPyV-DDR interface at a molecular scale is conducive to the development of novel therapeutic approaches.

Merkel cell polyomavirus and non-Merkel cell carcinomas: guilty or circumstantial evidence?

Merkel cell polyomavirus (MCPyV) is the major causative factor of the rare but aggressive cancer, Merkel cell carcinoma (MCC). Two characteristics of MCPyV-positive MCCs are integration of the viral genome and expression of a truncated version of one of its oncogenic proteins, namely large T antigen. The strong association of MCPyV with MCC development has incited researchers to further investigate a possible role of this virus in other cancers. However, many of the examples displaying the presence of the virus in the various non-MCC cancers are not able to clearly demonstrate a direct connection between cellular transformation and the presence of the virus. The prevalence of the virus is significantly lower in non-MCC cancers compared to MCCs, with a lower level of viral load and sparse viral protein expression. Moreover, the state of the viral genome, and whether a truncated large T antigen is expressed, has rarely been investigated. Nonetheless, considering the strong oncogenic potential of MCPyV proteins in MCC, the plausible contribution of MCPyV to transformation and cancer growth in non-MCC tumors cannot be ruled out. Furthermore, the absence of MCPyV in cancers does not exclude a hit-and-run mechanism, or the oncoproteins of MCPyV may potentiate the neoplastic process mediated by co-infecting oncoviruses such as high-risk human papillomaviruses and Epstein-Barr virus. The current review is focusing on the available data describing the presence of MCPyV in non-MCC tumors, with an aim to provide a comprehensive overview of the corresponding literature and to discuss the potential contribution of MCPyV to non-MCC cancer in light of this.

A polyomavirus peptide binds to the capsid VP1 pore and has potent antiviral activity against BK and JC polyomaviruses

In pursuit of therapeutics for human polyomaviruses, we identified a peptide derived from the BK polyomavirus (BKV) minor structural proteins VP2/3 that is a potent inhibitor of BKV infection with no observable cellular toxicity. The thirteen-residue peptide binds to major structural protein VP1 with single-digit nanomolar affinity. Alanine-scanning of the peptide identified three key residues, substitution of each of which results in ~1000 fold loss of binding affinity with a concomitant reduction in antiviral activity. Structural studies demonstrate specific binding of the peptide to the pore of pentameric VP1. Cell-based assays demonstrate nanomolar inhibition (EC₅₀) of BKV infection and suggest that the peptide acts early in the viral entry pathway. Homologous peptide exhibits similar binding to JC polyomavirus VP1 and inhibits infection with similar potency to BKV in a model cell line. Lastly, these studies validate targeting the VP1 pore as a novel strategy for the development of anti-polyomavirus agents.

Prisoners of war - host adaptation and its constraints on virus evolution

Recent discoveries of contemporary genotypes of hepatitis B virus and parvovirus B19 in ancient human remains demonstrate that little genetic change has occurred in these viruses over 4,500-6,000 years. Endogenous viral elements in host genomes provide separate evidence that viruses similar to many major contemporary groups circulated 100 million years ago or earlier. In this Opinion article, we argue that the extraordinary conservation of virus genome sequences is best explained by a niche-filling model in which fitness optimization is rapidly achieved in their specific hosts. Whereas short-term substitution rates reflect the accumulation of tolerated sequence changes within adapted genomes, longer-term rates increasingly resemble those of their hosts as the evolving niche moulds and effectively imprisons the virus in co-adapted virus-host relationships. Contrastingly, viruses that jump hosts undergo strong and stringent adaptive selection as they maximize their fit to their new niche. This adaptive capability may paradoxically create evolutionary stasis in long-term host relationships. While viruses can evolve and adapt rapidly, their hosts may ultimately shape their longer-term evolution.

Evolution and molecular epidemiology of polyomaviruses

Polyomaviruses (PyVs) are small DNA viruses that infect several species, including mammals, birds and fishes. Their study gained momentum after the report of previously unidentified viral species in the past decade, and especially, since the description of the first polyomavirus clearly oncogenic for humans. The aim of this work was to review the most relevant aspects of the evolution and molecular epidemiology of polyomaviruses, allowing to reveal general evolutionary patterns and to identify some unaddressed issues and future challenges. The main points analysed included: 1) the species and genera assignation criteria; 2) the hypotheses, mechanisms and timescale of the ancient and recent evolutionary history of polyomaviruses; and 3) the molecular epidemiology of human viruses, with special attention to JC, BK and Merkel cell polyomaviruses.

The case for BK polyomavirus as a cause of bladder cancer

In 2014, the International Agency for Research on Cancer judged Merkel cell polyomavirus (MCPyV) to be a probable human carcinogen. BK polyomavirus (BKPyV, a distant cousin of MCPyV) was ruled a possible carcinogen. In this review, we argue that it has recently become reasonable to view both of these viruses as known human carcinogens. In particular, several complementary lines of evidence support a causal role for BKPyV in the development of bladder carcinomas affecting organ transplant patients. The expansion of inexpensive deep sequencing has opened new approaches to investigating the important question of whether BKPyV causes urinary tract cancers in the general population.

Codon usage patterns of LT-Ag genes in polyomaviruses from different host species

Background

Polyomaviruses (PyVs) have a wide range of hosts, from humans to fish, and their effects on hosts vary. The differences in the infection characteristics of PyV with respect to the host are assumed to be influenced by the biochemical function of the LT-Ag protein, which is related to the cytopathic effect and tumorigenesis mechanism via interaction with the host protein.

Methods

We carried out a comparative analysis of codon usage patterns of large T-antigens (LT-Ags) of PyVs isolated from various host species and their functional domains and sequence motifs. Parity rule 2 (PR2) and neutrality analysis were applied to evaluate the effects of mutation and selection pressure on codon usage bias. To investigate evolutionary relationships among PyVs, we carried out a phylogenetic analysis, and a correspondence analysis of relative synonymous codon usage (RSCU) values was performed.

Results

Nucleotide composition analysis using LT-Ag gene sequences showed that the GC and GC3 values of avian PyVs were higher than those of mammalian PyVs. The effective number of codon (ENC) analysis showed host-specific ENC distribution characteristics in both the LT-Ag gene and the coding sequences of its domain regions. In the avian and fish PyVs, the codon diversity was significant, whereas the mammalian PyVs tended to exhibit conservative and host-specific evolution of codon usage bias. The results of our PR2 and neutrality analysis revealed mutation bias or highly variable GC contents by showing a narrow GC12 distribution and wide GC3 distribution in all sequences. Furthermore, the calculated RSCU values revealed differences in the codon usage preference of the LT-AG gene according to the host group. A similar tendency was observed in the two functional domains used in the analysis.

Conclusions

Our study showed that specific domains or sequence motifs of various PyV LT-Ags have evolved so that each virus protein interacts with host cell targets. They have also adapted to thrive in specific host species and cell types. Functional domains of LT-Ag, which are known to interact with host proteins involved in cell proliferation and gene expression regulation, may provide important information, as they are significantly related to the host specificity of PyVs.

Mash Screen: high-throughput sequence containment estimation for genome discovery

The MinHash algorithm has proven effective for rapidly estimating the resemblance of two genomes or metagenomes. However, this method cannot reliably estimate the containment of a genome within a metagenome. Here, we describe an online algorithm capable of measuring the containment of genomes and proteomes within either assembled or unassembled sequencing read sets. We describe several use cases, including contamination screening and retrospective analysis of metagenomes for novel genome discovery. Using this tool, we provide containment estimates for every NCBI RefSeq genome within every SRA metagenome and demonstrate the identification of a novel polyomavirus species from a public metagenome.

Single-Stranded DNA Viruses in Antarctic Cryoconite Holes

Antarctic cryoconite holes, or small melt-holes in the surfaces of glaciers, create habitable oases for isolated microbial communities with tightly linked microbial population structures. Viruses may influence the dynamics of polar microbial communities, but the viromes of the Antarctic cryoconite holes have yet to be characterized. We characterize single-stranded DNA (ssDNA) viruses from three cryoconite holes in the Taylor Valley, Antarctica, using metagenomics. Half of the assembled metagenomes cluster with those in the viral family Microviridae (n = 7), and the rest with unclassified circular replication associated protein (Rep)-encoding single-stranded (CRESS) DNA viruses (n = 7). An additional 18 virus-like circular molecules encoding either a Rep, a capsid protein gene, or other unidentified but viral-like open reading frames were identified. The samples from which the genomes were identified show a strong gradient in microbial diversity and abundances, and the number of viral genomes detected in each sample mirror that gradient. Additionally, one of the CRESS genomes assembled here shares ~90% genome-wide pairwise identity with a virus identified from a freshwater pond on the McMurdo Ice Shelf (Antarctica). Otherwise, the similarity of these viruses to those previously identified is relatively low. Together, these patterns are consistent with the presence of a unique regional virome present in fresh water host populations of the McMurdo Dry Valley region.

Virus-Host Coevolution with a Focus on Animal and Human DNA Viruses

Viruses have been infecting their host cells since the dawn of life, and this extremely long-term coevolution gave rise to some surprising consequences for the entire tree of life. It is hypothesised that viruses might have contributed to the formation of the first cellular life form, or that even the eukaryotic cell nucleus originates from an infection by a coated virus. The continuous struggle between viruses and their hosts to maintain at least a constant fitness level led to the development of an unceasing arms race, where weapons are often shuttled between the participants. In this literature review we try to give a short insight into some general consequences or traits of virus–host coevolution, and after this we zoom in to the viral clades of adenoviruses, herpesviruses, nucleo-cytoplasmic large DNA viruses, polyomaviruses and, finally, circoviruses.

Novel Polyomaviruses in Mammals from Multiple Orders and Reassessment of Polyomavirus Evolution and Taxonomy

As the phylogenetic organization of mammalian polyomaviruses is complex and currently incompletely resolved, we aimed at a deeper insight into their evolution by identifying polyomaviruses in host orders and families that have either rarely or not been studied. Sixteen unknown and two known polyomaviruses were identified in animals that belong to 5 orders, 16 genera, and 16 species. From 11 novel polyomaviruses, full genomes could be determined. Splice sites were predicted for large and small T antigen (LTAg, STAg) coding sequences (CDS) and examined experimentally in transfected cell culture. In addition, splice sites of seven published polyomaviruses were analyzed. Based on these data, LTAg and STAg annotations were corrected for 10/86 and 74/86 published polyomaviruses, respectively. For 25 polyomaviruses, a spliced middle T CDS was observed or predicted. Splice sites that likely indicate expression of additional, alternative T antigens, were experimentally detected for six polyomaviruses. In contrast to all other mammalian polyomaviruses, three closely related cetartiodactyl polyomaviruses display two introns within their LTAg CDS. In addition, the VP2 of Glis glis (edible dormouse) polyomavirus 1 was observed to be encoded by a spliced transcript, a unique experimental finding within the Polyomaviridae family. Co-phylogenetic analyses based on LTAg CDS revealed a measurable signal of codivergence when considering all mammalian polyomaviruses, most likely driven by relatively recent codivergence events. Lineage duplication was the only other process whose influence on polyomavirus evolution was unambiguous. Finally, our analyses suggest that an update of the taxonomy of the family is required, including the creation of novel genera of mammalian and non-mammalian polyomaviruses.

Characterization of the self-assembly of New Jersey polyomavirus VP1 into virus-like particles and the virus seroprevalence in Japan

New Jersey polyomavirus (NJPyV) was discovered in 2014 in a pancreatic transplant recipient’s vascular endothelial cells. Here, in the recombinant baculovirus system, VP1 protein of NJPyV expressed in insect cells was processed. The protein self-assembled into virus-like particles (NJPyV-LPs) in a cell-type-dependent manner, and the particles were then released into the culture media. Spherical ~50-nm-dia. NJPyV-LPs of uniform size with morphology resembling that of the native particles of polyomaviruses were purified from the fraction at 1.33 g/cm³ in supernatants of VP1-expressing Sf9 cells. We investigated the antigenic properties of purified NJPyV-LPs and performed a VLP-based enzyme immunoassay to determine the age-specific prevalence of NJPyV infection in a general Japanese population aged 1–70 years. The overall seropositivity rate of anti-NJPyV antibodies was only 1.8%. This might be explained by the low circulation of NJPyV in Japan. This is the first report of a large-scale serological survey of NJPyV in Asia (n = 1,050).

Non-human Primate Papillomaviruses Share Similar Evolutionary Histories and Niche Adaptation as the Human Counterparts

Given high genetic diversity of papillomaviruses (PV) and complex scenario of virus-host interaction, the genetic basis underlying the mechanisms of HPV carcinogenicity is not well understood. In an effort to evaluate the origin and evolution of PV pathogenicity, we collected paired oral, perianal, and genital swabs from a wild macaque population. Of the 117 surveyed macaques, 88 (75.2%) were positive for PV DNA in one or more sites, mostly common from genital swabs, followed by oral and perianal sites. All putative macaque PV types phylogenetically clustered into the genera Alpha-, Beta-, and Gammapapillomavirus, with a strong phylogeny-tropism association as observed in HPVs. Using a Bayesian Markov Chain Monte Carlo framework, we demonstrated ancient intra-host divergence of primate PVs in which multiple ancestors had split and adapted to specific host ecosystems at least 41 million years ago, prior to the speciation events of primate host species. Following subsequent divergence and niche adaptation, distinct but phylogenetically related PV types were transmitted to similar host ecosystems by closely related host animals when host speciation occurred, which may explain in part the origin of carcinogenicity of HPV type 16 (HPV16) and Macaca fascicularis PV type 3 (MfPV3) that evolved from a most recent common ancestor containing the determinants for cervicovaginal colonization and cervical cancer. The findings identifying evolutionary and biological relatedness between human and non-human primate PVs lay a genetic foundation for research on parasite-host interactions and carcinogenic outcomes, which will prove useful in further study of viral pathogenesis and host specificity.

Effect of the Large and Small T-Antigens of Human Polyomaviruses on Signaling Pathways

Viruses are intracellular parasites that require a permissive host cell to express the viral genome and to produce new progeny virus particles. However, not all viral infections are productive and some viruses can induce carcinogenesis. Irrespective of the type of infection (productive or neoplastic), viruses hijack the host cell machinery to permit optimal viral replication or to transform the infected cell into a tumor cell. One mechanism viruses employ to reprogram the host cell is through interference with signaling pathways. Polyomaviruses are naked, double-stranded DNA viruses whose genome encodes the regulatory proteins large T-antigen and small t-antigen, and structural proteins that form the capsid. The large T-antigens and small t-antigens can interfere with several host signaling pathways. In this case, we review the interplay between the large T-antigens and small t-antigens with host signaling pathways and the biological consequences of these interactions.

Human Polyomavirus JCPyV and Its Role in Progressive Multifocal Leukoencephalopathy and Oncogenesis

The human neurotropic virus JCPyV, a member of the Polyomaviridiae family, is the opportunistic infectious agent of Progressive Multifocal Leukoencephalopathy (PML), a fatal disease seen in severe immunosuppressive conditions and, during the last decade, in patients undergoing immunotherapy. JCPyV is a ubiquitous pathogen with up to 85% of the adult population word-wide exhibiting antibodies against it. Early experiments demonstrated that direct inoculation of JCPyV into the brain of different species resulted in the development of brain tumors and other neuroectodermal-derived neoplasias. Later, several reports showed the detection of viral sequences in medulloblastomas and glial tumors, as well as expression of the viral protein T-Antigen. Few oncogenic viruses, however, have caused so much controversy regarding their role in the pathogenesis of brain tumors, but the discovery of new Polyomaviruses that cause Merkel cell carcinomas in humans and brain tumors in racoons, in addition to the role of JCPyV in colon cancer and multiple mechanistic studies have shed much needed light on the role of JCPyV in cancer. The pathways affected by the viral protein T-Antigen include cell cycle regulators, like p53 and pRb, and transcription factors that activate pro-proliferative genes, like c-Myc. In addition, infection with JCPyV causes chromosomal damage and T-Antigen inhibits homologous recombination, and activates anti-apoptotic proteins, such as Survivin. Here we review the different aspects of the biology and physiopathology of JCPyV.

Lyon-IARC Polyomavirus DNA in Feces of Diarrheic Cats

A viral metagenomic analysis of feces from an unexplained outbreak of feline diarrhea revealed the presence of Lyon-IARC polyomavirus (LIPyV) DNA. LIPyV, whose genome was originally sequenced from swabs of human skin, was fecally shed by three out of five diarrheic cats.

A viral metagenomic analysis of feces from an unexplained outbreak of feline diarrhea revealed the presence of Lyon-IARC polyomavirus (LIPyV) DNA. LIPyV, whose genome was originally sequenced from swabs of human skin, was fecally shed by three out of five diarrheic cats.

Battistuzzi F, Kumar S. A Machine Learning Method for Detecting Autocorrelation of Evolutionary Rates in Large Phylogenies

New species arise from pre-existing species and inherit similar genomes and environments. This predicts greater similarity of the tempo of molecular evolution between direct ancestors and descendants, resulting in autocorrelation of evolutionary rates in the tree of life. Surprisingly, molecular sequence data have not confirmed this expectation, possibly because available methods lack the power to detect autocorrelated rates. Here, we present a machine learning method, CorrTest, to detect the presence of rate autocorrelation in large phylogenies. CorrTest is computationally efficient and performs better than the available state-of-the-art method. Application of CorrTest reveals extensive rate autocorrelation in DNA and amino acid sequence evolution of mammals, birds, insects, metazoans, plants, fungi, parasitic protozoans, and prokaryotes. Therefore, rate autocorrelation is a common phenomenon throughout the tree of life. These findings suggest concordance between molecular and nonmolecular evolutionary patterns, and they will foster unbiased and precise dating of the tree of life.

Merkel Cell Carcinoma in the HIV-1/AIDS Patient

Merkel cell carcinoma (MCC) is a highly aggressive, primary neuroendocrine cancer of the skin. The majority of MCC cases are associated with the recently discovered Merkel cell polyomavirus (MCPyV), while the remaining are caused by ultraviolet (UV) light-induced mutations from excessive sunlight exposure. The risk of developing MCC is much higher in the white population relative to all other races. Approximately 10% of all patients with MCC have some form of immunosuppression including HIV-1/AIDS, chronic inflammatory conditions, solid organ transplantation, or hematological malignancies. The age of onset of MCC is lower and the mortality is higher in immunosuppressed individuals than in immune-competent patients. It is plausible that HIV-1/AIDS predisposes to virus-positive MCC, but it should be noted that HIV-1/AIDS increases the risk for developing of UV-induced skin cancers such as cutaneous squamous cell carcinoma and basal cell carcinoma and therefore may also increase the risk for virus-negative MCC. Surgical management is considered standard of care for localized Merkel cell carcinoma with current recommendations advising a wide local excision of the lesion. Most international guidelines support the use of local adjuvant radiotherapy coupled with tumor staging to improve the frequency of cure. For advanced, metastatic, and recurrent MCC, checkpoint blockade inhibitors targeting PD-1 and PD-L1 have shown remarkable activity including durable long-term. MCC in patients living with HIV-1/AIDS are treated with similar modalities as HIV-1 uninfected individuals with MCC.

Prevalence of a vertically transmitted single-stranded DNA virus in spinybacked orbweavers (Gasteracantha cancriformis) from Florida, USA

Spiders (order Araneae, class Arachnida) are an important group of predatory arthropods in terrestrial ecosystems that have been recently identified as an untapped reservoir of single-stranded (ss)DNA viruses. Specifically, spiders harbour a diversity of ssDNA viruses encoding a replication-associated protein (Rep) within a circular genome. However, little is known about the ecology of novel circular Rep-encoding ssDNA (CRESS DNA) viruses. Here we investigated two CRESS DNA viruses recently identified in spinybacked orbweavers (Gasteracantha cancriformis), namely spinybacked orbweaver circular virus (SpOrbCV) 1 and 2. SpOrbCV-1 was detected in the majority (> 65 %) of spider specimens from all life stages, including eggs, spiderlings and adults, demonstrating that this virus is active within spinybacked orbweavers. In contrast, SpOrbCV-2 was only detected in adults at a lower (36 %) prevalence. Since we also detected SpOrbCV-2 in other spider species and this virus has been reported from a dragonfly, we suggest that SpOrbCV-2 is accumulated in these predators through common insect prey. The prevalence of SpOrbCV-1 in collected specimens allowed us to design assays to characterize this virus, which represents a new group of CRESS DNA viruses, the 'circularisviruses'. To our knowledge, SpOrbCV-1 is the first example of a vertically transmitted virus in spiders, which may explain its high prevalence in spinybacked orbweavers. Since vertically transmitted viruses infecting insects (class Insecta) can manipulate their host's behaviour and physiology, future studies should investigate the ecological role of vertically transmitted viruses in spiders.

Fecal Viral Diversity of Captive and Wild Tasmanian Devils Characterized Using Virion-Enriched Metagenomics and Metatranscriptomics

The Tasmanian devil is an endangered carnivorous marsupial threatened by devil facial tumor disease (DFTD). While research on DFTD has been extensive, little is known about viruses in devils and whether any are of potential conservation relevance for this endangered species. Using both metagenomics based on virion enrichment and sequence-independent amplification (virion-enriched metagenomics) and metatranscriptomics based on bulk RNA sequencing, we characterized and compared the fecal viromes of captive and wild devils. A total of 54 fecal samples collected from two captive and four wild populations were processed for virome characterization using both approaches. In total, 24 novel marsupial-related viruses, comprising a sapelovirus, astroviruses, rotaviruses, picobirnaviruses, parvoviruses, papillomaviruses, polyomaviruses, and a gammaherpesvirus, were identified, as well as known mammalian pathogens such as rabbit hemorrhagic disease virus 2. Captive devils showed significantly lower viral diversity than wild devils. Comparison of the two virus discovery approaches revealed substantial differences in the number and types of viruses detected, with metatranscriptomics better suited for RNA viruses and virion-enriched metagenomics largely identifying more DNA viruses. Thus, the viral communities revealed by virion-enriched metagenomics and metatranscriptomics were not interchangeable and neither approach was able to detect all viruses present. An integrated approach using both virion-enriched metagenomics and metatranscriptomics constitutes a powerful tool for obtaining a complete overview of both the taxonomic and functional profiles of viral communities within a sample.IMPORTANCE The Tasmanian devil is an iconic Australian marsupial that has suffered an 80% population decline due to a contagious cancer, devil facial tumor disease, along with other threats. Until now, viral discovery in this species has been confined to one gammaherpesvirus (dasyurid herpesvirus 2 [DaHV-2]), for which captivity was identified as a significant risk factor. Our discovery of 24 novel marsupial-associated RNA and DNA viruses, and that viral diversity is lower in captive than in wild devils, has greatly expanded our knowledge of gut-associated viruses in devils and provides important baseline information that will contribute to the conservation and captive management of this endangered species. Our results also revealed that a combination of virion-enriched metagenomics and metatranscriptomics may be a more comprehensive approach for virome characterization than either method alone. Our results thus provide a springboard for continuous improvements in the way we study complex viral communities.

Fecal Viral Diversity of Captive and Wild Tasmanian Devils Characterized Using Virion-Enriched Metagenomics and Metatranscriptomics

The Tasmanian devil is an iconic Australian marsupial that has suffered an 80% population decline due to a contagious cancer, devil facial tumor disease, along with other threats. Until now, viral discovery in this species has been confined to one gammaherpesvirus (dasyurid herpesvirus 2 [DaHV-2]), for which captivity was identified as a significant risk factor. Our discovery of 24 novel marsupial-associated RNA and DNA viruses, and that viral diversity is lower in captive than in wild devils, has greatly expanded our knowledge of gut-associated viruses in devils and provides important baseline information that will contribute to the conservation and captive management of this endangered species. Our results also revealed that a combination of virion-enriched metagenomics and metatranscriptomics may be a more comprehensive approach for virome characterization than either method alone. Our results thus provide a springboard for continuous improvements in the way we study complex viral communities.

The Tasmanian devil is an endangered carnivorous marsupial threatened by devil facial tumor disease (DFTD). While research on DFTD has been extensive, little is known about viruses in devils and whether any are of potential conservation relevance for this endangered species. Using both metagenomics based on virion enrichment and sequence-independent amplification (virion-enriched metagenomics) and metatranscriptomics based on bulk RNA sequencing, we characterized and compared the fecal viromes of captive and wild devils. A total of 54 fecal samples collected from two captive and four wild populations were processed for virome characterization using both approaches. In total, 24 novel marsupial-related viruses, comprising a sapelovirus, astroviruses, rotaviruses, picobirnaviruses, parvoviruses, papillomaviruses, polyomaviruses, and a gammaherpesvirus, were identified, as well as known mammalian pathogens such as rabbit hemorrhagic disease virus 2. Captive devils showed significantly lower viral diversity than wild devils. Comparison of the two virus discovery approaches revealed substantial differences in the number and types of viruses detected, with metatranscriptomics better suited for RNA viruses and virion-enriched metagenomics largely identifying more DNA viruses. Thus, the viral communities revealed by virion-enriched metagenomics and metatranscriptomics were not interchangeable and neither approach was able to detect all viruses present. An integrated approach using both virion-enriched metagenomics and metatranscriptomics constitutes a powerful tool for obtaining a complete overview of both the taxonomic and functional profiles of viral communities within a sample.

IMPORTANCE The Tasmanian devil is an iconic Australian marsupial that has suffered an 80% population decline due to a contagious cancer, devil facial tumor disease, along with other threats. Until now, viral discovery in this species has been confined to one gammaherpesvirus (dasyurid herpesvirus 2 [DaHV-2]), for which captivity was identified as a significant risk factor. Our discovery of 24 novel marsupial-associated RNA and DNA viruses, and that viral diversity is lower in captive than in wild devils, has greatly expanded our knowledge of gut-associated viruses in devils and provides important baseline information that will contribute to the conservation and captive management of this endangered species. Our results also revealed that a combination of virion-enriched metagenomics and metatranscriptomics may be a more comprehensive approach for virome characterization than either method alone. Our results thus provide a springboard for continuous improvements in the way we study complex viral communities.

Detection and genome characterization of two novel papillomaviruses and a novel polyomavirus in tree shrew (Tupaia belangeri chinensis) in China

Background

Papillomaviruses (PVs) and polyomaviruses (PyVs) infect diverse vertebrates including human and cause a broad spectrum of outcomes from asymptomatic infection to severe disease. There has been no PV and only one PyV detected in tree shrews, though the genomic properties of tree shrews are highly similar to those of the primates.

Methods

Swab and organ samples of tree shrews collected in the Yunnan Province of China, were tested by viral metagenomic analysis and random PCR to detect the presence of PVs and PyVs. By PCR amplification using specific primers, cloning, sequencing and assembling, genomes of two PVs and one PyV were identified in the samples.

Results

Two novel PVs and a novel PyV, named tree shrew papillomavirus 1 and 2 (TbelPV1 and TbelPV2) and polyomavirus 1 (TbelPyV1) were characterized in the Chinese tree shrew (Tupaia belangeri chinensis). The genomes of TbelPV1, TbelPV2, and TbelPyV1 are 7410 bp, 7526 bp, and 4982 bp in size, respectively. The TbelPV1 genome contains 7 putative open-reading frames (ORFs) coding for viral proteins E1, E2, E4, E6, E7, L1, and L2; the TbelPV2 genome contains 6 ORFs coding for viral proteins E1, E2, E6, E7, L1, and L2; and the TbelPyV1 genome codes for the typical small and large T antigens of PyV, as well as the VP1, VP2, and VP3 capsid proteins. Genomic comparison and phylogenetic analysis indicated that TbelPV1 and TbelPV2 represented 2 novel PV genera of Papillomaviridae, and TbelPyV1 represented a new species of genus Alphapolyomavirus. Our epidemiologic study indicated that TbelPV1 and TbelPV2 were both detected in oral swabs, while TbelPyV1 was detected in oral swabs and spleens.

Conclusion

Two novel PVs (TbelPV1 and TbelPV2) and a novel PyV (TbelPyV) were discovered in tree shrews and their genomes were characterized. TbelPV1, TbelPV2, and TbelPyV1 have the highest similarity to Human papillomavirus type 63, Ursus maritimus papillomavirus 1, and Human polyomavirus 9, respectively. TbelPV1 and TbelPV2 only showed oral tropism, while TbelPyV1 showed oral and spleen tropism.

Electronic supplementary material

The online version of this article (10.1186/s12985-019-1141-9) contains supplementary material, which is available to authorized users.

Polyomavirus T Antigen Induces APOBEC3B Expression Using an LXCXE-Dependent and TP53-Independent Mechanism

APOBEC3B is a single-stranded DNA cytosine deaminase with beneficial innate antiviral functions. However, misregulated APOBEC3B can also be detrimental by inflicting APOBEC signature C-to-T and C-to-G mutations in genomic DNA of multiple cancer types. Polyomavirus and papillomavirus oncoproteins induce APOBEC3B overexpression, perhaps to their own benefit, but little is known about the cellular mechanisms hijacked by these viruses to do so. Here we investigate the molecular mechanism of APOBEC3B upregulation by the polyomavirus large T antigen. First, we demonstrate that the upregulated APOBEC3B enzyme is strongly nuclear and partially localized to virus replication centers. Second, truncated T antigen (truncT) is sufficient for APOBEC3B upregulation, and the RB-interacting motif (LXCXE), but not the p53-binding domain, is required. Third, genetic knockdown of RB1 alone or in combination with RBL1 and/or RBL2 is insufficient to suppress truncT-mediated induction of APOBEC3B Fourth, CDK4/6 inhibition by palbociclib is also insufficient to suppress truncT-mediated induction of APOBEC3B Last, global gene expression analyses in a wide range of human cancers show significant associations between expression of APOBEC3B and other genes known to be regulated by the RB/E2F axis. These experiments combine to implicate the RB/E2F axis in promoting APOBEC3B transcription, yet they also suggest that the polyomavirus RB-binding motif has at least one additional function in addition to RB inactivation for triggering APOBEC3B upregulation in virus-infected cells.IMPORTANCE The APOBEC3B DNA cytosine deaminase is overexpressed in many different cancers and correlates with elevated frequencies of C-to-T and C-to-G mutations in 5'-TC motifs, oncogene activation, acquired drug resistance, and poor clinical outcomes. The mechanisms responsible for APOBEC3B overexpression are not fully understood. Here, we show that the polyomavirus truncated T antigen (truncT) triggers APOBEC3B overexpression through its RB-interacting motif, LXCXE, which in turn likely modulates the binding of E2F family transcription factors to promote APOBEC3B expression. This work strengthens the mechanistic linkage between active cell cycling, APOBEC3B overexpression, and cancer mutagenesis. Although this mutational mechanism damages cellular genomes, viruses may leverage it to promote evolution, immune escape, and pathogenesis. The cellular portion of the mechanism may also be relevant to nonviral cancers, where genetic mechanisms often activate the RB/E2F axis and APOBEC3B mutagenesis contributes to tumor evolution.

Metagenomic Discovery of 83 New Human Papillomavirus Types in Patients with Immunodeficiency

Although some members of the viral family Papillomaviridae cause benign skin warts (papillomas), many human papillomavirus (HPV) infections are not associated with visible symptoms. For example, most healthy adults chronically shed Gammapapillomavirus (Gamma) virions from apparently healthy skin surfaces. To further explore the diversity of papillomaviruses, we performed viromic surveys on immunodeficient individuals suffering from florid skin warts. Our results nearly double the number of known Gamma HPV types and suggest that WHIM syndrome patients are uniquely susceptible to Gamma HPV-associated skin warts. Preliminary results suggest that treatment with the drug plerixafor may promote resolution of the unusual Gamma HPV skin warts observed in WHIM patients.

Several immunodeficiencies are associated with high susceptibility to persistent and progressive human papillomavirus (HPV) infection leading to a wide range of cutaneous and mucosal lesions. However, the HPV types most commonly associated with such clinical manifestations in these patients have not been systematically defined. Here, we used virion enrichment, rolling circle amplification, and deep sequencing to identify circular DNA viruses present in skin swabs and/or wart biopsy samples from 48 patients with rare genetic immunodeficiencies, including patients with warts, hypogammaglobulinemia, infections, myelokathexis (WHIM) syndrome, or epidermodysplasia verruciformis (EV). Their profiles were compared with the profiles of swabs from 14 healthy adults and warts from 6 immunologically normal children. Individual patients were typically infected with multiple HPV types; up to 26 different types were isolated from a single patient (multiple anatomical sites, one time point). Among these, we identified the complete genomes of 83 previously unknown HPV types and 35 incomplete genomes representing possible additional new types. HPV types in the genus Gammapapillomavirus were common in WHIM patients, whereas EV patients mainly shed HPVs from the genus Betapapillomavirus. Preliminary evidence based on three WHIM patients treated with plerixafor, a leukocyte mobilizing agent, suggest that longer-term therapy may correlate with decreased HPV diversity and increased predominance of HPV types associated with childhood skin warts.

IMPORTANCE Although some members of the viral family Papillomaviridae cause benign skin warts (papillomas), many human papillomavirus (HPV) infections are not associated with visible symptoms. For example, most healthy adults chronically shed Gammapapillomavirus (Gamma) virions from apparently healthy skin surfaces. To further explore the diversity of papillomaviruses, we performed viromic surveys on immunodeficient individuals suffering from florid skin warts. Our results nearly double the number of known Gamma HPV types and suggest that WHIM syndrome patients are uniquely susceptible to Gamma HPV-associated skin warts. Preliminary results suggest that treatment with the drug plerixafor may promote resolution of the unusual Gamma HPV skin warts observed in WHIM patients.

Genetic variability and oncogenic risk association of human papillomavirus type 58 E6 and E7 genes in Taizhou area, China

It is well recognized that the association of human papillomavirus (HPV) and cervical carcinogenesis is based on the presence of HPV DNA sequence. The E6 and E7 oncoproteins encoded by high-risk HPV types play a key role in carcinogenesis. HPV58 type accounts for a larger share of cervical disease in China, whereas data on HPV58 genetic variability in China is limited. We aimed to evaluate the diversity of HPV58 genetic variants by sequencing the entire E6 and E7 genes. Phylogenetic trees were constructed by Maximum likelihood method by MEGA 5.05 software. In this study, the overall HPV infection rate was 22.6% (2891/12780) in Southeast China and the prevalence of HPV58 infection rate was 2.6% (335/12780). 26 nucleotides substitutions were observed in E6 and E7 genes with 10 novel substitutions and 17 non-synonymous substitutions. We obtained 25 distinct variation patterns which the accession GenBank numbers as MH348918-MH348942. All of HPV58 variants belong to lineage A, while no lineage B, C and D were detected in Taizhou area, Southeast China. The sublineage A1, A2, and A3 variants were found in 136 (68.3%), 39 (19.6%), and 24 (12.1%) of HPV58 isolates, respectively. The sublineage A3 variants with T20I/G63S substitutions at E7 oncoprotein carried a significantly higher risk for high-grade cervical intraepithelial neoplasia (CIN2 or worse, CIN2+) when compared with other HPV58 variants (odds ratio = 4.41, P < 0.05). Nevertheless, there was no association between HPV58 (sub) lineages and cervical lesions. These data provide the critical characteristics of HPV58 variants to assist further investigation of carcinogenic association and the development of next generation vaccines and diagnostic assays in China.

Sialic Acids in Nonenveloped Virus Infections

Sialic acid-based glycoconjugates cover the surfaces of many different cell types, defining key properties of the cell surface such as overall charge or likely interaction partners. Because of this prominence, sialic acids play prominent roles in mediating attachment and entry to viruses belonging to many different families. In this review, we first describe how interactions between viruses and sialic acid-based glycan structures can be identified and characterized using a range of techniques. We then highlight interactions between sialic acids and virus capsid proteins in four different viruses, and discuss what these interactions have taught us about sialic acid engagement and opportunities to interfere with binding.

Niche adaptation and viral transmission of human papillomaviruses from archaic hominins to modern humans

Recent discoveries on the origins of modern humans from multiple archaic hominin populations and the diversity of human papillomaviruses (HPVs) suggest a complex scenario of virus-host evolution. To evaluate the origin of HPV pathogenesis, we estimated the phylogeny, timing, and dispersal of HPV16 variants using a Bayesian Markov Chain Monte Carlo framework. To increase precision, we identified and characterized non-human primate papillomaviruses from New and Old World monkeys to set molecular clock models. We demonstrate specific host niche adaptation of primate papillomaviruses with subsequent coevolution with their primate hosts for at least 40 million years. Analyses of 212 HPV16 complete genomes and 3582 partial sequences estimated ancient divergence of HPV16 variants (between A and BCD lineages) from their most recent common ancestors around half a million years ago, roughly coinciding with the timing of the split between archaic Neanderthals and modern Homo sapiens, and nearly three times longer than divergence times of modern Homo sapiens. HPV16 A lineage variants were significantly underrepresented in present African populations, whereas the A sublineages were highly prevalent in European (A1-3) and Asian (A4) populations, indicative of viral sexual transmission from Neanderthals to modern non-African humans through multiple interbreeding events in the past 80 thousand years. Remarkably, the human leukocyte antigen B*07:02 and C*07:02 alleles associated with increased risk in cervix cancer represent introgressed regions from Neanderthals in present-day Eurasians. The archaic hominin-host-switch model was also supported by other HPV variants. Niche adaptation and virus-host codivergence appear to influence the pathogenesis of papillomaviruses.

Virus discovery in all three major lineages of terrestrial arthropods highlights the diversity of single-stranded DNA viruses associated with invertebrates

Viruses encoding a replication-associated protein (Rep) within a covalently closed, single-stranded (ss)DNA genome are among the smallest viruses known to infect eukaryotic organisms, including economically valuable agricultural crops and livestock. Although circular Rep-encoding ssDNA (CRESS DNA) viruses are a widespread group for which our knowledge is rapidly expanding, biased sampling toward vertebrates and land plants has limited our understanding of their diversity and evolution. Here, we screened terrestrial arthropods for CRESS DNA viruses and report the identification of 44 viral genomes and replicons associated with specimens representing all three major terrestrial arthropod lineages, namely Euchelicerata (spiders), Hexapoda (insects), and Myriapoda (millipedes). We identified virus genomes belonging to three established CRESS DNA viral families (Circoviridae, Genomoviridae, and Smacoviridae); however, over half of the arthropod-associated viral genomes are only distantly related to currently classified CRESS DNA viral sequences. Although members of viral and satellite families known to infect plants (Geminiviridae, Nanoviridae, Alphasatellitidae) were not identified in this study, these plant-infecting CRESS DNA viruses and replicons are transmitted by hemipterans. Therefore, members from six out of the seven established CRESS DNA viral families circulate among arthropods. Furthermore, a phylogenetic analysis of Reps, including endogenous viral sequences, reported to date from a wide array of organisms revealed that most of the known CRESS DNA viral diversity circulates among invertebrates. Our results highlight the vast and unexplored diversity of CRESS DNA viruses among invertebrates and parallel findings from RNA viral discovery efforts in undersampled taxa.

Rat polyomavirus 2 infection in a colony of X-linked severe combined immunodeficiency rats in Japan

Polyomaviruses (PyVs) infect a wide range of animals and provoke wasting diseases, particularly in immunosuppressed hosts. Recently, a novel Rattus norvegicus polyomavirus 2 (RatPyV2) has been identified in a colony of X-linked severe combined immunodeficiency (X-SCID) rats in the United States. Here, we describe the first report of the RatPyV2 infection in an X-SCID rat colony in Japan. The affected rats exhibited adult-onset wasting. Histologically, we observed large basophilic intranuclear inclusion bodies within the hyperplastic or dysplastic epithelial cells in the salivary glands, Harderian glands, extraorbital lacrimal glands, and in respiratory and reproductive tissues. Among these organs, the parotid salivary, Harderian, and extraorbital lacrimal glands were most obviously affected. In particular, the parotid salivary glands were the most severely and diffusely affected and atrophic lesions were prominent even at 1 month of age, which suggested that the parotid salivary glands would be highly susceptible to RatPyV2 in X-SCID rats. RatPyV2 inclusion bodies were also detected in the tail of the epididymis and deferent duct. Such reproductive lesions developed significantly in the later stage of breeding age, and therefore may be associated with the reduced fecundity observed in the infected X-SCID rats. We also established a simple, rapid, and non-invasive diagnostic method based on the Amp-FTA method, using buccal swabs for the detection of RatPyV2 in immunodeficient rats. Our findings contribute to the early detection and diagnosis of RatPyV2 infections.

Microscopic and Molecular Evidence of the First Elasmobranch Adomavirus, the Cause of Skin Disease in a Giant Guitarfish, Rhynchobatus djiddensis

Only eight families of double-stranded DNA (dsDNA) viruses are known to infect vertebrate animals. During an investigation of papillomatous skin disease in an elasmobranch species, the giant guitarfish (Rhynchobatus djiddensis), a novel virus, distinct from all known viral families in regard to particle size, morphology, genome organization, and helicase phylogeny was discovered. Large inclusion bodies containing 75-nm icosahedral viral particles were present within epithelial cell nuclei in the proliferative skin lesions. Deep metagenomic sequencing revealed a 22-kb circular dsDNA viral genome, tentatively named guitarfish "adomavirus" (GAdoV), with only distant homology to two other fish viruses, Japanese eel endothelial cell-infecting virus (JEECV) and a recently reported marbled eel virus. Phylogenetic analysis of the helicase domain places the guitarfish virus in a novel clade that is equidistant between members of the Papillomaviridae and Polyomaviridae families. Specific PCR, quantitative PCR, and in situ hybridization were used to detect, quantify, and confirm that GAdoV DNA was localized to affected epithelial cell nuclei. Changes in the viral titer, as well as the presence of a hybridization signal, coincided with the progression and then final resolution of gross and microscopic lesions. The results indicate that GAdoV is the causative agent of the proliferative skin lesions.IMPORTANCE Cartilaginous fish, including the sharks and rays, evolved from ancestral fish species at least 400 million years ago. Even though they are the descendants of one of the most ancient vertebrate lineages, reports of viral diseases in these species are rare and poorly documented. Deep sequencing revealed a highly divergent virus, tentatively named guitarfish adomavirus, that is distantly related to known papillomaviruses and polyomaviruses. Out of the eight predicted viral genes, only the helicase could be identified as viral by sequence homology searches (BLAST), exemplifying the difficulties of discovering novel viruses within seas of unidentifiable "dark matter" associated with deep sequencing data. The novel adomavirus represents the first viral genome shown to cause clinical disease in a cartilaginous fish species, the giant guitarfish. Our findings demonstrate that emerging fish viruses are fertile ground to expand our understanding of viral evolution in vertebrates.

Characterization of BK Polyomaviruses from Kidney Transplant Recipients Suggests a Role for APOBEC3 in Driving In-Host Virus Evolution

BK polyomavirus (BKV) frequently causes nephropathy (BKVN) in kidney transplant recipients (KTRs). BKV has also been implicated in the etiology of bladder and kidney cancers. We characterized BKV variants from two KTRs who developed BKVN followed by renal carcinoma. Both patients showed a swarm of BKV sequence variants encoding non-silent mutations in surface loops of the viral major capsid protein. The temporal appearance and disappearance of these mutations highlights the intra-patient evolution of BKV. Some of the observed mutations conferred resistance to antibody-mediated neutralization. The mutations also modified the spectrum of receptor glycans engaged by BKV during host cell entry. Intriguingly, all observed mutations were consistent with DNA damage caused by antiviral APOBEC3 cytosine deaminases. Moreover, APOBEC3 expression was evident upon immunohistochemical analysis of renal biopsies from KTRs. These results provide a snapshot of in-host BKV evolution and suggest that APOBEC3 may drive BKV mutagenesis in vivo.

Phylodynamics of Merkel-cell polyomavirus and human polyomavirus 6: A long-term history with humans

New human polyomaviruses have been described in the last years, including the Merkel-cell polyomavirus (MCPyV; Human polyomavirus 5) and the Human polyomavirus 6 (HPyV6). Although their infection is usually asymptomatic, in immunocompromised host can cause life-threatening pathologies, such as the Merkel cell carcinoma, an aggressive skin neoplasia associated to the MCPyV. Despite being prevalent viruses in population, epidemiological data from South America are scarce, as well as the characterization of the viral types circulating and their origin. The aims of this work were to describe MCPyV and HPyV6 from environmental samples with different geographical origin and to analyze their phylogenetic and evolutionary histories, particularly for MCPyV. Partial and complete genome sequences were obtained from sewage samples from Argentina, Uruguay and Spain. A total number of 87 sequences were obtained for MCPyV and 33 for HPyV6. Phylogenetic analysis showed that MCPyV sequences distributed according to their geographic origin in Europe/North America, Africa, Asia, South America and Oceania groups, suggesting that viral diversification might have followed human migrations across the globe. In particular, viruses from Argentina associated with Europe/North America and South America genotypes, whereas those from Uruguay and Spain also grouped with Africa genotype, reflecting the origin of the current population in each country, which could arrive not only during ancient human migration but also during recent migratory events. In addition, the South American group presented a high level of clusterization, showing internal clusters that could be related to specific locations, such as French Guiana and Brazil or the Southern region into South America, such as Argentina and Uruguay, suggesting a long term evolutionary process in the region. Additionally, in this work, we carried out the first analysis about the evolutionary history of MCPyV trough the integration of phylogenetic, epidemiological and historical data. Since a strong association is observed between the phylogenetic relationships and the origin of the sampled population, this analysis was based on the hypothesis of co-divergence between the virus and human populations. This analysis resulted in a substitution rate of 5.1 × 10^-8 s/s/y (∼5.1% of divergence per million years) for the complete genome of MCPyV, which is in the range of those estimated for other double-stranded DNA viruses. Regarding HPyV6, a South American group with clusterization was observed (sequences from Uruguay). Meanwhile, sequences from Argentina grouped with European ones (France and Spain) and remained separated from those isolated in China, USA or Australia. The analysis of viruses from the environment allowed us to deep characterize prevalent infections in different geographic regions, reveling that viruses circulating in each population reflected its origin and that there are specific lineages associated with South America.