Concurrence of Iridovirus, Polyomavirus, and a Unique Member of a New Group of Fish Papillomaviruses in Lymphocystis Disease-Affected Gilthead Sea Bream

The Role of aDNA in Understanding the Coevolutionary Patterns of Human Sexually Transmitted Infections

Analysis of pathogen genome data sequenced from clinical and historical samples has made it possible to perform phylogenetic analyses of sexually transmitted infections on a global scale, and to estimate the diversity, distribution, and coevolutionary host relationships of these pathogens, providing insights into pathogen emergence and disease prevention. Deep-sequenced pathogen genomes from clinical studies and ancient samples yield estimates of within-host and between-host evolutionary rates and provide data on changes in pathogen genomic stability and evolutionary responses. Here we examine three groups of pathogens transmitted mainly through sexual contact between modern humans to provide insight into ancient human behavior and history with their pathogens. Exploring ancient pathogen genomic divergence and the ancient viral-host parallel evolutionary histories will help us to reconstruct the origin of present-day geographical distribution and diversity of clinical pathogen infections, and will hopefully allow us to foresee possible environmentally induced pathogen evolutionary responses. Lastly, we emphasize that ancient pathogen DNA research should be combined with modern clinical pathogen data, and be equitable and provide advantages for all researchers worldwide, e.g., through shared data.

Molecular identification of papillomavirus in ducks

Papillomaviruses infect many vertebrates, including birds. Persistent infections by some strains can cause malignant proliferation of cells (i.e. cancer), though more typically infections cause benign tumours, or may be completely subclinical. Sometimes extensive, persistent tumours are recorded-notably in chaffinches and humans. In 2016, a novel papillomavirus genotype was characterized from a duck faecal microbiome, in Bhopal, India; the sixth papillomavirus genotype from birds. Prompted by this finding, we screened 160 cloacal swabs and 968 faecal samples collected from 299 ducks sampled at Ottenby Bird Observatory, Sweden in 2015, using a newly designed real-time PCR. Twenty one samples (1.9%) from six individuals (2%) were positive. Eighteen sequences were identical to the published genotype, duck papillomavirus 1. One additional novel genotype was recovered from three samples. Both genotypes were recovered from a wild strain domestic mallard that was infected for more than 60 days with each genotype. All positive individuals were adult (P = 0.004). Significantly more positive samples were detected from swabs than faecal samples (P < 0.0001). Sample type data suggests transmission may be via direct contact, and only infrequently, via the oral-faecal route. Infection in only adult birds supports the hypothesis that this virus is sexually transmitted, though more work is required to verify this.

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.

Oncogenic human papillomaviruses

Human papillomaviruses (HPVs) are an ancient group of viruses with small, double-stranded DNA circular genomes. They are species-specific and have a strict tropism for mucosal and cutaneous stratified squamous epithelial surfaces of the host. A subset of these viruses has been demonstrated to be the causative agent of several human cancers. Here, we review the biology, natural history, evolution and cancer association of the oncogenic HPVs.This article is part of the themed issue 'Human oncogenic viruses'.

Invertebrate Iridoviruses: A Glance over the Last Decade

Members of the family Iridoviridae (iridovirids) are large dsDNA viruses that infect both invertebrate and vertebrate ectotherms and whose symptoms range in severity from minor reductions in host fitness to systemic disease and large-scale mortality. Several characteristics have been useful for classifying iridoviruses; however, novel strains are continuously being discovered and, in many cases, reliable classification has been challenging. Further impeding classification, invertebrate iridoviruses (IIVs) can occasionally infect vertebrates; thus, host range is often not a useful criterion for classification. In this review, we discuss the current classification of iridovirids, focusing on genomic and structural features that distinguish vertebrate and invertebrate iridovirids and viral factors linked to host interactions in IIV6 (Invertebrate iridescent virus 6). In addition, we show for the first time how complete genome sequences of viral isolates can be leveraged to improve classification of new iridovirid isolates and resolve ambiguous relations. Improved classification of the iridoviruses may facilitate the identification of genus-specific virulence factors linked with diverse host phenotypes and host interactions.

Fish polyomaviruses belong to two distinct evolutionary lineages

The Polyomaviridae is a diverse family of circular double-stranded DNA viruses. Polyomaviruses have been isolated from a wide array of animal hosts. An understanding of the evolutionary and ecological dynamics of these viruses is essential to understanding the pathogenicity of polyomaviruses. Using a high throughput sequencing approach, we identified a novel polyomavirus in an emerald notothen (Trematomus bernacchii) sampled in the Ross sea (Antarctica), expanding the known number of fish-associated polyomaviruses. Our analysis suggests that polyomaviruses belong to three main evolutionary clades; the first clade is made up of all recognized terrestrial polyomaviruses. The fish-associated polyomaviruses are not monophyletic, and belong to two divergent evolutionary lineages. The fish viruses provide evidence that the evolution of the key viral large T protein involves gain and loss of distinct domains.

Complete genome sequence of shrimp hemocyte iridescent virus (SHIV) isolated from white leg shrimp, Litopenaeus vannamei

Infection with shrimp hemocyte iridescent virus (SHIV), a new virus of the family Iridoviridae isolated in China, results in a high mortality rate in white leg shrimp (Litopenaeus vannamei). The complete genome sequence of SHIV was determined and analyzed in this study. The genomic DNA was 165,809 bp long with 34.6% G+C content and 170 open reading frames (ORFs). Dotplot analysis showed that the longest repetitive region was 320 bp in length, including 11 repetitions of an 18-bp sequence and 3.1 repetitions of a 39-bp sequence. Two phylogenetic trees were constructed based on 27 or 16 concatenated sequences of proteins encoded by genes that are conserved between SHIV homologous and other iridescent viruses. The results of this study, suggest that SHIV should be considered a member of the proposed new genus "Xiairidovirus".

Diverse papillomaviruses identified in Weddell seals

Papillomaviridae is a diverse family of circular, double-stranded DNA (dsDNA) viruses that infect a broad range of mammalian, avian and fish hosts. While papillomaviruses have been characterized most extensively in humans, the study of non-human papillomaviruses has contributed greatly to our understanding of their pathogenicity and evolution. Using high-throughput sequencing approaches, we identified 7 novel papillomaviruses from vaginal swabs collected from 81 adult female Weddell seals (Leptonychotes weddellii) in the Ross Sea of Antarctica between 2014-2017. These seven papillomavirus genomes were amplified from seven individual seals, and six of the seven genomes represented novel species with distinct evolutionary lineages. This highlights the diversity of papillomaviruses among the relatively small number of Weddell seal samples tested. Viruses associated with large vertebrates are poorly studied in Antarctica, and this study adds information about papillomaviruses associated with Weddell seals and contributes to our understanding of the evolutionary history of papillomaviruses.

Viral insulin-like peptides activate human insulin and IGF-1 receptor signaling: A paradigm shift for host-microbe interactions

Viruses are the most abundant biological entities and carry a wide variety of genetic material, including the ability to encode host-like proteins. Here we show that viruses carry sequences with significant homology to several human peptide hormones including insulin, insulin-like growth factors (IGF)-1 and -2, FGF-19 and -21, endothelin-1, inhibin, adiponectin, and resistin. Among the strongest homologies were those for four viral insulin/IGF-1-like peptides (VILPs), each encoded by a different member of the family Iridoviridae VILPs show up to 50% homology to human insulin/IGF-1, contain all critical cysteine residues, and are predicted to form similar 3D structures. Chemically synthesized VILPs can bind to human and murine IGF-1/insulin receptors and stimulate receptor autophosphorylation and downstream signaling. VILPs can also increase glucose uptake in adipocytes and stimulate the proliferation of fibroblasts, and injection of VILPs into mice significantly lowers blood glucose. Transfection of mouse hepatocytes with DNA encoding a VILP also stimulates insulin/IGF-1 signaling and DNA synthesis. Human microbiome studies reveal the presence of these Iridoviridae in blood and fecal samples. Thus, VILPs are members of the insulin/IGF superfamily with the ability to be active on human and rodent cells, raising the possibility for a potential role of VILPs in human disease. Furthermore, since only 2% of viruses have been sequenced, this study raises the potential for discovery of other viral hormones which, along with known virally encoded growth factors, may modify human health and disease.

Structural Insights in Multifunctional Papillomavirus Oncoproteins

Since their discovery in the mid-eighties, the main papillomavirus oncoproteins E6 and E7 have been recalcitrant to high-resolution structure analysis. However, in the last decade a wealth of three-dimensional information has been gained on both proteins whether free or complexed to host target proteins. Here, we first summarize the diverse activities of these small multifunctional oncoproteins. Next, we review the available structural data and the new insights they provide about the evolution of E6 and E7, their multiple interactions and their functional variability across human papillomavirus (HPV) species.

Ancient DNA provides evidence of 27,000-year-old papillomavirus infection and long-term codivergence with rodents

The long-term evolutionary history of many viral lineages is poorly understood. Novel sources of ancient DNA combined with phylogenetic analyses can provide insight into the time scale of virus evolution. Here we report viral sequences from ancient North American packrat middens. We screened samples up to 27,000-years old and found evidence of papillomavirus (PV) infection in Neotoma cinerea (Bushy-tailed packrat). Phylogenetic analysis placed the PV sequences in a clade with other previously published PV sequences isolated from rodents. Concordance between the host and virus tree topologies along with a correlation in branch lengths suggests a shared evolutionary history between rodents and PVs. Based on host divergence times, PVs have likely been circulating in rodents for at least 17 million years. These results have implications for our understanding of PV evolution and for further research with ancient DNA from Neotoma middens.

Unique genome organization of non-mammalian papillomaviruses provides insights into the evolution of viral early proteins

The family Papillomaviridae contains more than 320 papillomavirus types, with most having been identified as infecting skin and mucosal epithelium in mammalian hosts. To date, only nine non-mammalian papillomaviruses have been described from birds (n = 5), a fish (n = 1), a snake (n = 1), and turtles (n = 2). The identification of papillomaviruses in sauropsids and a sparid fish suggests that early ancestors of papillomaviruses were already infecting the earliest Euteleostomi. The Euteleostomi clade includes more than 90 per cent of the living vertebrate species, and progeny virus could have been passed on to all members of this clade, inhabiting virtually every habitat on the planet. As part of this study, we isolated a novel papillomavirus from a 16-year-old female Adélie penguin (Pygoscelis adeliae) from Cape Crozier, Ross Island (Antarctica). The new papillomavirus shares ∼64 per cent genome-wide identity to a previously described Adélie penguin papillomavirus. Phylogenetic analyses show that the non-mammalian viruses (expect the python, Morelia spilota, associated papillomavirus) cluster near the base of the papillomavirus evolutionary tree. A papillomavirus isolated from an avian host (Northern fulmar; Fulmarus glacialis), like the two turtle papillomaviruses, lacks a putative E9 protein that is found in all other avian papillomaviruses. Furthermore, the Northern fulmar papillomavirus has an E7 more similar to the mammalian viruses than the other avian papillomaviruses. Typical E6 proteins of mammalian papillomaviruses have two Zinc finger motifs, whereas the sauropsid papillomaviruses only have one such motif. Furthermore, this motif is absent in the fish papillomavirus. Thus, it is highly likely that the most recent common ancestor of the mammalian and sauropsid papillomaviruses had a single motif E6. It appears that a motif duplication resulted in mammalian papillomaviruses having a double Zinc finger motif in E6. We estimated the divergence time between Northern fulmar-associated papillomavirus and the other Sauropsid papillomaviruses be to around 250 million years ago, during the Paleozoic-Mesozoic transition and our analysis dates the root of the papillomavirus tree between 400 and 600 million years ago. Our analysis shows evidence for niche adaptation and that these non-mammalian viruses have highly divergent E6 and E7 proteins, providing insights into the evolution of the early viral (onco-)proteins.

HPV8-E6 Interferes with Syntenin-2 Expression through Deregulation of Differentiation, Methylation and Phosphatidylinositide-Kinase Dependent Mechanisms

The E6 oncoproteins of high-risk human papillomaviruses (HPV) of genus alpha contain a short peptide sequence at the carboxy-terminus, the PDZ binding domain, with which they interact with the corresponding PDZ domain of cellular proteins. Interestingly, E6 proteins from papillomaviruses of genus beta (betaPV) do not encode a comparable PDZ binding domain. Irrespective of this fact, we previously showed that the E6 protein of HPV8 (betaPV type) could circumvent this deficit by targeting the PDZ protein Syntenin-2 through transcriptional repression (Lazic et al., 2012). Despite its high binding affinity to phosphatidylinositol-4,5-bisphosphate (PI(4,5)P₂), very little is known about Syntenin-2. This study aimed to extend the knowledge on Syntenin-2 and how its expression is controlled. We now identified that Syntenin-2 is expressed at high levels in differentiating and in lower amounts in keratinocytes cultured in serum-free media containing low calcium concentration. HPV8-E6 led to a further reduction of Syntenin-2 expression only in cells cultured in low calcium. In the skin of patients suffering from Epidermodysplasia verruciformis, who are predisposed to betaPV infection, Syntenin-2 was expressed in differentiating keratinocytes of non-lesional skin, but was absent in virus positive squamous tumors. Using 5-Aza-2′-deoxycytidine, which causes DNA demethylation, Syntenin-2 transcription was profoundly activated and fully restored in the absence and presence of HPV8-E6, implicating that E6 mediated repression of Syntenin-2 transcription is due to promoter hypermethylation. Since Syntenin-2 binds to PI(4,5)P₂, we further tested whether the PI(4,5)P₂ metabolic pathway might govern Syntenin-2 expression. PI(4,5)P₂ is generated by the activity of phosphatidylinositol-4-phosphate-5-kinase type I (PIP5KI) or phosphatidylinositol-5-phosphate-4-kinase type II (PIP4KII) isoforms α, β and γ. Phosphatidylinositide kinases have recently been identified as regulators of gene transcription. Surprisingly, transfection of siRNAs directed against PIP5KI and PIP4KII resulted in higher Syntenin-2 expression with the highest effect mediated by siPIP5KIα. HPV8-E6 was able to counteract siPIP4KIIα, siPIP4KIIβ and siPIP5KIγ mediated Syntenin-2 re-expression but not siPIP5KIα. Finally, we identified Syntenin-2 as a key factor regulating PIP5KIα expression. Collectively, our data demonstrates that Syntenin-2 is regulated through multiple mechanisms and that downregulation of Syntenin-2 expression may contribute to E6 mediated dedifferentiation of infected skin cells.

The Enigmatic Origin of Papillomavirus Protein Domains

Almost a century has passed since the discovery of papillomaviruses. A few decades of research have given a wealth of information on the molecular biology of papillomaviruses. Several excellent studies have been performed looking at the long- and short-term evolution of these viruses. However, when and how papillomaviruses originate is still a mystery. In this study, we systematically searched the (sequenced) biosphere to find distant homologs of papillomaviral protein domains. Our data show that, even including structural information, which allows us to find deeper evolutionary relationships compared to sequence-only based methods, only half of the protein domains in papillomaviruses have relatives in the rest of the biosphere. We show that the major capsid protein L1 and the replication protein E1 have relatives in several viral families, sharing three protein domains with Polyomaviridae and Parvoviridae. However, only the E1 replication protein has connections with cellular organisms. Most likely, the papillomavirus ancestor is of marine origin, a biotope that is not very well sequenced at the present time. Nevertheless, there is no evidence as to how papillomaviruses originated and how they became vertebrate and epithelium specific.

Mechanisms and strategies of papillomavirus replication

Animal and human papillomaviruses (HPVs) replicate persistently in specific types of stratified epithelia of their host. After the initial infection, the viral genome replicates at low levels in the dividing cells of the epithelium, and these cells form a reservoir of infection that can last for decades. When the infected cells differentiate, viral genomes replicate to high levels to form progeny virus that is released from the surface of the epithelium. This complex life cycle requires several different modes of viral DNA replication, but papillomaviruses are masters at hijacking key cellular processes to facilitate their own reproduction.

Contribution of Next-Generation Sequencing to Aquatic and Fish Virology

The recent technological advances in nucleic acid sequencing, called next-generation sequencing (NGS), have revolutionized the field of genomics and have also influenced viral research. Aquatic viruses, and especially those infecting fish, have also greatly benefited from NGS technologies, which provide a huge amount of molecular information at a low cost in a relatively short period of time. Here, we review the use of the current high-throughput sequencing platforms with a special focus on the associated challenges (regarding sample preparation and bioinformatics) in their applications to the field of aquatic virology, especially for: (i) discovering novel viruses that may be associated with fish mortalities, (ii) elucidating the mechanisms of pathogenesis, and finally (iii) studying the molecular epidemiology of these pathogens.

Recombination events and variability among full-length genomes of co-circulating molluscum contagiosum virus subtypes 1 and 2

Molluscum contagiosum virus (MCV) is the sole member of the Molluscipoxvirus genus and causes a highly prevalent human disease of the skin characterized by the formation of a variable number of lesions that can persist for prolonged periods of time. Two major genotypes, subtype 1 and subtype 2, are recognized, although currently only a single complete genomic sequence corresponding to MCV subtype 1 is available. Using next-generation sequencing techniques, we report the complete genomic sequence of four new MCV isolates, including the first one derived from a subtype 2. Comparisons suggest a relatively distant evolutionary split between both MCV subtypes. Further, our data illustrate concurrent circulation of distinct viruses within a population and reveal the existence of recombination events among them. These results help identify a set of MCV genes with potentially relevant roles in molluscum contagiosum epidemiology and pathogenesis.

Identification of a polyomavirus in Weddell seal (Leptonychotes weddellii) from the Ross Sea (Antarctica)

Viruses are ubiquitous in nature, however, very few have been identified that are associated with Antarctic animals. Here we report the identification of a polyomavirus in the kidney tissue of a deceased Weddell seal from the Ross Sea, Antarctica. The circular genome (5186 nt) has typical features of polyomaviruses with a small and larger T-antigen open reading frames (ORFs) and three ORFs encoding VP1, VP2 and VP3 capsid proteins. The genome of the Weddell seal polyomavirus (WsPyV) shares 85.4% genome-wide pairwise identity with a polyomavirus identified in a California sea lion. To our knowledge WsPyV is the first viral genome identified in Antarctic pinnipeds and the third polyomavirus to be identified from an Antarctic animal, the other two being from Adélie penguin (Pygoscelis adeliae) and a sharp-spined notothen (Trematomus pennellii), both sampled in the Ross sea. The GenBank accession number: KX533457.

The Papillomavirus Episteme: a major update to the papillomavirus sequence database

The Papillomavirus Episteme (PaVE) is a database of curated papillomavirus genomic sequences, accompanied by web-based sequence analysis tools. This update describes the addition of major new features. The papillomavirus genomes within PaVE have been further annotated, and now includes the major spliced mRNA transcripts. Viral genes and transcripts can be visualized on both linear and circular genome browsers. Evolutionary relationships among PaVE reference protein sequences can be analysed using multiple sequence alignments and phylogenetic trees. To assist in viral discovery, PaVE offers a typing tool; a simplified algorithm to determine whether a newly sequenced virus is novel. PaVE also now contains an image library containing gross clinical and histopathological images of papillomavirus infected lesions. Database URL: https://pave.niaid.nih.gov/.

Transmission between Archaic and Modern Human Ancestors during the Evolution of the Oncogenic Human Papillomavirus 16

Every human suffers through life a number of papillomaviruses (PVs) infections, most of them asymptomatic. A notable exception are persistent infections by Human papillomavirus 16 (HPV16), the most oncogenic infectious agent for humans and responsible for most infection-driven anogenital cancers. Oncogenic potential is not homogeneous among HPV16 lineages, and genetic variation within HPV16 exhibits some geographic structure. However, an in-depth analysis of the HPV16 evolutionary history was still wanting. We have analyzed extant HPV16 diversity and compared the evolutionary and phylogeographical patterns of humans and of HPV16. We show that codivergence with modern humans explains at most 30% of the present viral geographical distribution. The most explanatory scenario suggests that ancestral HPV16 already infected ancestral human populations and that viral lineages co-diverged with the hosts in parallel with the split between archaic Neanderthal-Denisovans and ancestral modern human populations, generating the ancestral HPV16A and HPV16BCD viral lineages, respectively. We propose that after out-of-Africa migration of modern human ancestors, sexual transmission between human populations introduced HPV16A into modern human ancestor populations. We hypothesize that differential coevolution of HPV16 lineages with different but closely related ancestral human populations and subsequent host-switch events in parallel with introgression of archaic alleles into the genomes of modern human ancestors may be largely responsible for the present-day differential prevalence and association with cancers for HPV16 variants.