The first complete papillomavirus genome characterized from a marsupial host: a novel isolate from Bettongia penicillata

Widespread Horizontal Gene Transfer Among Animal Viruses

The initial objective of this study was to shed light on the evolution of small DNA tumor viruses by analyzing de novo assemblies of publicly available deep sequencing datasets. The survey generated a searchable database of contig snapshots representing more than 100,000 Sequence Read Archive records. Using modern structure-aware search tools, we iteratively broadened the search to include an increasingly wide range of other virus families. The analysis revealed a surprisingly diverse range of chimeras involving different virus groups. In some instances, genes resembling known DNA-replication modules or known virion protein operons were paired with unrecognizable sequences that structural predictions suggest may represent previously unknown replicases and novel virion architectures. Discrete clades of an emerging group called adintoviruses were discovered in datasets representing humans and other primates. As a proof of concept, we show that the contig database is also useful for discovering RNA viruses and candidate archaeal phages. The ancillary searches revealed additional examples of chimerization between different virus groups. The observations support a gene-centric taxonomic framework that should be useful for future virus-hunting efforts.

Divergent hepaciviruses, delta-like viruses, and a chu-like virus in Australian marsupial carnivores (dasyurids)

Although Australian marsupials are characterised by unique biology and geographic isolation, little is known about the viruses present in these iconic wildlife species. The Dasyuromorphia are an order of marsupial carnivores found only in Australia that include both the extinct Tasmanian tiger (thylacine) and the highly threatened Tasmanian devil. Several other members of the order are similarly under threat of extinction due to habitat loss, hunting, disease, and competition and predation by introduced species such as feral cats. We utilised publicly available RNA-seq data from the National Center for Biotechnology Information (NCBI) Sequence Read Archive (SRA) database to document the viral diversity within four Dasyuromorph species. Accordingly, we identified fifteen novel virus sequences from five DNA virus families (Adenoviridae, Anelloviridae, Gammaherpesvirinae, Papillomaviridae, and Polyomaviridae) and three RNA virus taxa: the order Jingchuvirales, the genus Hepacivirus, and the delta-like virus group. Of particular note was the identification of a marsupial-specific clade of delta-like viruses that may indicate an association of deltaviruses with marsupial species. In addition, we identified a highly divergent hepacivirus in a numbat liver transcriptome that falls outside of the larger mammalian clade. We also detect what may be the first Jingchuvirales virus in a mammalian host-a chu-like virus in Tasmanian devils-thereby expanding the host range beyond invertebrates and ectothermic vertebrates. As many of these Dasyuromorphia species are currently being used in translocation efforts to reseed populations across Australia, understanding their virome is of key importance to prevent the spread of viruses to naive populations.

Identification of a novel polyomavirus from a marsupial host

We report the identification and analysis of a full sequence of a novel polyomavirus from a brushtail possum (Trichosurus vulpecula ) termed possum polyomavirus (PPyV). The sequence was obtained from the next-generation sequencing assembly during an investigation into the aetiological agent for a neurological disease of possums termed wobbly possum disease (WPD), but the virus was not aetiologically involved in WPD. The PPyV genome was 5,224 nt long with the organisation typical for polyomaviruses, including early (large and small T antigens) and late (Viral Protein 1 (VP1), VP2, and VP3) coding regions separated by the non-coding control region of 465 nt. PPyV clustered with betapolyomaviruses in the WUKI clade but showed less than 60 per cent identity to any of the members of this clade. We propose that PPyV is classified within a new species in the genus Betapolyomavirus . These data add to our limited knowledge of marsupial viruses and their evolution.

Ancient viral integrations in marsupials: a potential antiviral defence

Marsupial viruses are understudied compared to their eutherian mammal counterparts, although they may pose severe threats to vulnerable marsupial populations. Genomic viral integrations, termed 'endogenous viral elements' (EVEs), could protect the host from infection. It is widely known past viral infections and EVEs play an active role in antiviral defence in invertebrates and plants. This study aimed to characterise actively transcribed EVEs in Australian marsupial species, because they may play an integral role in cellular defence against viruses. This study screened publicly available RNA sequencing data sets (n = 35) and characterised 200 viral transcripts from thirteen Australian marsupial species. Of the 200 transcripts, 188 originated from either Bornaviridae, Filoviridae, or Parvoviridae EVEs. The other twelve transcripts were from putative active infections from members of the Herpesviridae and Anelloviridae, and Hepadnaviridae. EVE transcripts (n = 188) were mapped to marsupial genomes (where available, n = 5/13) to identify the genomic insertion sites. Of the 188 transcripts, 117 mapped to 39 EVEs within the koala, bare-nosed wombat, tammar wallaby, brushtail possum, and Tasmanian devil genomes. The remaining eight animals had no available genome (transcripts n = 71). Every marsupial has Bornaviridae, Filoviridae, and Parvoviridae EVEs, a trend widely observed in eutherian mammals. Whilst eutherian bornavirus EVEs are predominantly nucleoprotein-derived, marsupial bornavirus EVEs demonstrate a surprising replicase gene bias. We predicted these widely distributed EVEs were conserved within marsupials from ancient germline integrations, as many were over 65 million years old. One bornavirus replicase EVE, present in six marsupial genomes, was estimated to be 160 million years old, predating the American-Australian marsupial split. We considered transcription of these EVEs through small non-coding RNA as an ancient viral defence. Consistent with this, in koala small RNA sequence data sets, we detected Bornaviridae replicase and Filoviridae nucleoprotein produced small RNA. These were enriched in testis tissue, suggesting they could protect marsupials from vertically transmitted viral integrations.

Origin and evolution of papillomavirus (onco)genes and genomes

Papillomaviruses (PVs) are ancient viruses infecting vertebrates, from fishes to mammals. Although the genomes of PVs are small and show conserved synteny, PVs display large genotypic diversity and ample variation in the phenotypic presentation of the infection. Most PV genomes contain two small early genes E6 and E7. In a bunch of closely related human papillomaviruses (HPVs), the E6 and E7 proteins provide the viruses with oncogenic potential. The recent discoveries of PVs without E6 and E7 in different fish species place a new root on the PV tree, and suggest that ancestral PVs consisted of the minimal PV backbone E1-E2-L2-L1. Bayesian phylogenetic analyses date the most recent common ancestor of the PV backbone to 424 million years ago (Ma). Common ancestry tests on extant E6 and E7 genes indicate that they share a common ancestor dating back to at least 184 Ma. In AlphaPVs infecting Old World monkeys and apes, the appearance of the E5 oncogene 53-58 Ma concurred with (i) a significant increase in substitution rate, (ii) a basal radiation and (iii) key gain of functions in E6 and E7. This series of events was instrumental to construct the extant phenotype of oncogenic HPVs. Our results assemble the current knowledge on PV diversity and present an ancient evolutionary timeline punctuated by evolutionary innovations in the history of this successful viral family. This article is part of the theme issue 'Silent cancer agents: multi-disciplinary modelling of human DNA oncoviruses'.

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.

Association of papillomavirus E6 proteins with either MAML1 or E6AP clusters E6 proteins by structure, function, and evolutionary relatedness

Papillomavirus E6 proteins bind to LXXLL peptide motifs displayed on targeted cellular proteins. Alpha genus HPV E6 proteins associate with the cellular ubiquitin ligase E6AP (UBE3A), by binding to an LXXLL peptide (ELTLQELLGEE) displayed by E6AP, thereby stimulating E6AP ubiquitin ligase activity. Beta, Gamma, and Delta genera E6 proteins bind a similar LXXLL peptide (WMSDLDDLLGS) on the cellular transcriptional co-activator MAML1 and thereby repress Notch signaling. We expressed 45 different animal and human E6 proteins from diverse papillomavirus genera to ascertain the overall preference of E6 proteins for E6AP or MAML1. E6 proteins from all HPV genera except Alpha preferentially interacted with MAML1 over E6AP. Among animal papillomaviruses, E6 proteins from certain ungulate (SsPV1 from pigs) and cetacean (porpoises and dolphins) hosts functionally resembled Alpha genus HPV by binding and targeting the degradation of E6AP. Beta genus HPV E6 proteins functionally clustered with Delta, Pi, Tau, Gamma, Chi, Mu, Lambda, Iota, Dyokappa, Rho, and Dyolambda E6 proteins to bind and repress MAML1. None of the tested E6 proteins physically and functionally interacted with both MAML1 and E6AP, indicating an evolutionary split. Further, interaction of an E6 protein was insufficient to activate degradation of E6AP, indicating that E6 proteins that target E6AP co-evolved to separately acquire both binding and triggering of ubiquitin ligase activation. E6 proteins with similar biological function clustered together in phylogenetic trees and shared structural features. This suggests that the divergence of E6 proteins from either MAML1 or E6AP binding preference is a major event in papillomavirus evolution.

Temporal Patterns in the Abundance of a Critically Endangered Marsupial Relates to Disturbance by Roads and Agriculture

The aim of this study was to investigate how landscape disturbance associated with roads, agriculture and forestry influenced temporal patterns in woylie (Bettongia penicillata) abundance before, during and after periods of rapid population change. Data were collected from an area of approximately 140,000 ha of forest within the Upper Warren region in south-western Australia. Woylie abundance was measured using cage trapping at 22 grid and five transect locations with varying degrees of landscape disturbance between 1994 and 2012. We found evidence that the distribution and abundance of woylies over time appears to be related to the degree of fragmentation by roads and proximity to agriculture. Sites furthest from agriculture supported a greater abundance of woylies and had slower rates of population decline. Sites with fewer roads had a greater abundance of woylies generally and a greater rate of increase in abundance after the implementation of invasive predator control. The results of this study suggest that landscape disturbance is less important at peak population densities, but during times of environmental and population change, sites less dissected by roads and agriculture better support woylie populations. This may be due to the role these factors play in increasing the vulnerability of woylies to introduced predators, population fragmentation, weed species invasion, mortality from road collisions or a reduction in available habitat. Strategies that reduce the impact of disturbance on woylie populations could include the rationalisation of forest tracks and consolidation of contiguous habitat through the acquisition of private property. Reducing the impact of disturbance in the Upper Warren region could improve the resilience of this critically important woylie population during future environmental change.

The Ancient Evolutionary History of Polyomaviruses

Polyomaviruses are a family of DNA tumor viruses that are known to infect mammals and birds. To investigate the deeper evolutionary history of the family, we used a combination of viral metagenomics, bioinformatics, and structural modeling approaches to identify and characterize polyomavirus sequences associated with fish and arthropods. Analyses drawing upon the divergent new sequences indicate that polyomaviruses have been gradually co-evolving with their animal hosts for at least half a billion years. Phylogenetic analyses of individual polyomavirus genes suggest that some modern polyomavirus species arose after ancient recombination events involving distantly related polyomavirus lineages. The improved evolutionary model provides a useful platform for developing a more accurate taxonomic classification system for the viral family Polyomaviridae.

Polyomaviruses are a family of DNA-based viruses that are known to infect various terrestrial vertebrates, including humans. In this report, we describe our discovery of highly divergent polyomaviruses associated with various marine fish. Searches of public deep sequencing databases unexpectedly revealed the existence of polyomavirus-like sequences in scorpion and spider datasets. Our analysis of these new sequences suggests that polyomaviruses have slowly co-evolved with individual host animal lineages through an established mechanism known as intrahost divergence. The proposed model is similar to the mechanisms through with other DNA viruses, such as papillomaviruses, are thought to have evolved. Our analysis also suggests that distantly related polyomaviruses sometimes recombine to produce new chimeric lineages. We propose a possible taxonomic scheme that can account for these inferred ancient recombination events.

A Century of Shope Papillomavirus in Museum Rabbit Specimens

Sylvilagus floridanus Papillomavirus (SfPV) causes growth of large horn-like tumors on rabbits. SfPV was described in cottontail rabbits (probably Sylvilagus floridanus) from Kansas and Iowa by Richard Shope in 1933, and detected in S. audubonii in 2011. It is known almost exclusively from the US Midwest. We explored the University of Kansas Natural History Museum for historical museum specimens infected with SfPV, using molecular techniques, to assess if additional wild species host SfPV, and whether SfPV occurs throughout the host range, or just in the Midwest. Secondary aims were to detect distinct strains, and evidence for strain spatio-temporal specificity. We found 20 of 1395 rabbits in the KU collection SfPV symptomatic. Three of 17 lagomorph species (S. nuttallii, and the two known hosts) were symptomatic, while Brachylagus, Lepus and eight additional Sylvilagus species were not. 13 symptomatic individuals were positive by molecular testing, including the first S. nuttallii detection. Prevalence of symptomatic individuals was significantly higher in Sylvilagus (1.8%) than Lepus. Half of these specimens came from Kansas, though new molecular detections were obtained from Jalisco—Mexico’s first—and Nebraska, Nevada, New Mexico, and Texas, USA. We document the oldest lab-confirmed case (Kansas, 1915), pre-dating Shope’s first case. SfPV amplification was possible from 63.2% of symptomatic museum specimens. Using multiple methodologies, rolling circle amplification and, multiple isothermal displacement amplification in addition to PCR, greatly improved detection rates. Short sequences were obtained from six individuals for two genes. L1 gene sequences were identical to all previously detected sequences; E7 gene sequences, were more variable, yielding five distinct SfPV1 strains that differing by less than 2% from strains circulating in the Midwest and Mexico, between 1915 and 2005. Our results do not clarify whether strains are host species specific, though they are consistent with SfPV specificity to genus Sylvilagus.

Papillomaviruses: Viral evolution, cancer and evolutionary medicine

Papillomaviruses (PVs) are a numerous family of small dsDNA viruses infecting virtually all mammals. PVs cause infections without triggering a strong immune response, and natural infection provides only limited protection against reinfection. Most PVs are part and parcel of the skin microbiota. In some cases, infections by certain PVs take diverse clinical presentations from highly productive self-limited warts to invasive cancers. We propose PVs as an excellent model system to study the evolutionary interactions between the immune system and pathogens causing chronic infections: genotypically, PVs are very diverse, with hundreds of different genotypes infecting skin and mucosa; phenotypically, they display extremely broad gradients and trade-offs between key phenotypic traits, namely productivity, immunogenicity, prevalence, oncogenicity and clinical presentation. Public health interventions have been launched to decrease the burden of PV-associated cancers, including massive vaccination against the most oncogenic human PVs, as well as systematic screening for PV chronic anogenital infections. Anti-PVs vaccines elicit protection against infection, induce cross-protection against closely related viruses and result in herd immunity. However, our knowledge on the ecological and intrapatient dynamics of PV infections remains fragmentary. We still need to understand how the novel anthropogenic selection pressures posed by vaccination and screening will affect viral circulation and epidemiology. We present here an overview of PV evolution and the connection between PV genotypes and the phenotypic, clinical manifestations of the diseases they cause. This differential link between viral evolution and the gradient cancer-warts-asymptomatic infections makes PVs a privileged playground for evolutionary medicine research.

Sudden and rapid decline of the abundant marsupial Bettongia penicillata in Australia

The woylie Bettongia penicillata is categorized as Critically Endangered, having declined by c. 90% between 1999 and 2006. The decline continues and the cause is not fully understood. Within a decline diagnosis framework we characterized the nature of the decline and identified potential causes, with a focus on the species’ largest populations, located in south-west Western Australia. We described the spatio-temporal pattern of the decline, and several attributes that are common across sites. We categorized the potential causes of the decline as resources, predators, disease and direct human interference. Based on the available evidence the leading hypothesis is that disease may be making woylies more vulnerable to predation but this remains to be tested. No substantial recoveries have been sustained to date, and one of the three remaining indigenous populations now appears to be extinct. Therefore, verifying the factors causing the decline and those limiting recovery is becoming increasingly urgent. Active adaptive management can be used to test putative agents, such as introduced predators. Insurance populations and ecological monitoring should also be included in an integrated conservation and management strategy for the species.

Human papillomavirus genome variants

Amongst the human papillomaviruses (HPVs), the genus Alphapapillomavirus contains HPV types that are uniquely pathogenic. They can be classified into species and types based on genetic distances between viral genomes. Current circulating infectious HPVs constitute a set of viral genomes that have evolved with the rapid expansion of the human population. Viral variants were initially identified through restriction enzyme polymorphisms and more recently through sequence determination of viral fragments. Using partial sequence information, the history of variants, and the association of HPV variants with disease will be discussed with the main focus on the recent utilization of full genome sequence information for variant analyses. The use of multiple sequence alignments of complete viral genomes and phylogenetic analyses have begun to define variant lineages and sublineages using empirically defined differences of 1.0-10.0% and 0.5-1.0%, respectively. These studies provide the basis to define the genetics of HPV pathogenesis.

Animal papillomaviruses

We provide an overview of the host range, taxonomic classification and genomic diversity of animal papillomaviruses. The complete genomes of 112 non-human papillomavirus types, recovered from 54 different host species, are currently available in GenBank. The recent characterizations of reptilian papillomaviruses extend the host range of the Papillomaviridae to include all amniotes. Although the genetically diverse papillomaviruses have a highly conserved genomic lay-out, deviations from this prototypic genome organization are observed in several animal papillomaviruses, and only the core ORFs E1, E2, L2 and L1 are present in all characterized papillomavirus genomes. The discovery of papilloma-polyoma hybrids BPCV1 and BPCV2, containing a papillomaviral late region but an early region encoding typical polyomaviral nonstructural proteins, and the detection of recombination breakpoints between the early and late coding regions of cetacean papillomaviruses, could indicate that early and late gene cassettes of papillomaviruses are relatively independent entities that can be interchanged by recombination.

Taxonomy and phylogeny of papillomaviruses: an overview and recent developments

For more than 30 years, papillomaviruses are standing in the center of medical and molecular interest as they cause several important cancers in humans. Research of the sheer unlimited number of different papillomavirus genomes, their host specificity and slow mutation rate is an important a branch of these efforts and has led to fascinating insight into the phylogeny of a virus family that can be traced back for several 100 million years.

Characterization of a novel papillomavirus species (ZcPV1) from two California sea lions (Zalophus californianus)

A seven-year old California sea lion (Zalophus californianus) presented with focally extensive, bilaterally symmetric, proliferative axillary skin lesions and preputial lesions. A second California sea lion in the same population presented with similar proliferative lesions on the underside of the tail. Histopathology revealed epidermal hyperplasia with severe hyperkeratosis, with proliferating keratinocytes forming broad, branching pegs that extended into the dermis. Pan-papillomaviral consensus PCR was used to obtain initial E1 sequence template and the complete genome was determined using a combination of rolling circle amplification and specific-primer PCR. Analysis revealed a novel papillomavirus, Zalophus californianus papillomavirus 1 (ZcPV1), with seven open reading frames encoding five early proteins (E6, E7, E1, E2 and E4) and two late proteins (L1 and L2). Phylogenetic analysis revealed that (ZcPV1) is most closely related to Equine papillomavirus 1 (EcPV1) in the genus Zetapapillomavirus, and Canine papillomaviruses 3 and 4 (CPV3, CPV4) in the genus Chipapillomavirus. The lesions regressed without intervention over a period of several months.

Classification and nomenclature system for human Alphapapillomavirus variants: general features, nucleotide landmarks and assignment of HPV6 and HPV11 isolates to variant lineages

BACKGROUND

Papillomaviruses constitute a family of viruses that can be classified into genera, species and types based on their viral genome heterogeneity. Currently circulating infectious human Alphapapillomaviruses (alpha-PVs) constitute a set of viral genomes that have evolved from archaic times and display features of host co-speciation. Viral variants are more recently evolved genomes that require a standardized classification and nomenclature.

OBJECTIVES

To describe a system for the classification and nomenclature of HPV viral variants and provide landmarks for the numbering of nucleotide positions.

METHODS

The complete 8 kb genomes of the alpha-9 species group and HPV6 and 11 types, collected from isolates throughout the world were obtained from published reports and GenBank. Complete genomes for each HPV type were aligned using the E1 start codon and sequence divergence was calculated by global and pairwise alignments using the MUSCLE program. Phylogenetic trees were constructed from the aligned sequences using a maximum likelihood method (RAxML).

RESULTS

Pairwise comparisons of nucleotide differences between complete genomes of each type from alpha-9 HPV isolates (HPV16, 31, 33, 35, 52, 58 and 67) revealed a trimodal distribution. Maximum heterogeneity for variants within a type varied from 0.6%-2.3%. Nucleotide differences of approximately 1.0%-10.0% and 0.5%-1.0% of the complete genomes were used to define variant lineages and sublineages, respectively. Analysis of 43 HPV6 complete genomes indicated the presence of 2 variant lineages, whereas 32 HPV11 isolates were highly similar and clustered into 2 sublineages. A table was constructed of the human alpha-PV landmark nucleotide sequences for future reference and alignments.

CONCLUSIONS

A proposed nomenclature system for viral variants and coordination of nucleotide positions will facilitate the comparison of variants across geographic regions and amongst different populations. In addition, this system will facilitate study of pathogenic, tissue tropism and functional differences amongst variant lineages of and polymorphisms within HPV variants.

Insights into Polyomaviridae microRNA function derived from study of the bandicoot papillomatosis carcinomatosis viruses

Several different members of the Polyomaviridae, including some human pathogens, encode microRNAs (miRNAs) that lie antisense with respect to the early gene products, the tumor (T) antigens. These miRNAs negatively regulate T antigen expression by directing small interfering RNA (siRNA)-like cleavage of the early transcripts. miRNA mutant viruses of some members of the Polyomaviridae express increased levels of early proteins during lytic infection. However, the importance of miRNA-mediated negative regulation of the T antigens remains uncertain. Bandicoot papillomatosis carcinomatosis virus type 1 (BPCV1) is associated with papillomas and carcinomas in the endangered marsupial the western barred bandicoot (Perameles bougainville). BPCV1 is the founding member of a new group of viruses that remarkably share distinct properties in common with both the polyomavirus and papillomavirus families. Here, we show that BPCV1 encodes, in the same orientation as the papillomavirus-like transcripts, a miRNA located within a long noncoding region (NCR) of the genome. Furthermore, this NCR serves the function of both promoter and template for the primary transcript that gives rise to the miRNA. Unlike the polyomavirus miRNAs, the BPCV1 miRNA is not encoded antisense to the T antigen transcripts but rather lies in a separate, proximal region of the genome. We have mapped the 3' untranslated region (UTR) of the BPCV1 large T antigen early transcript and identified a functional miRNA target site that is imperfectly complementary to the BPCV1 miRNA. Chimeric reporters containing the entire BPCV1 T antigen 3' UTR undergo negative regulation when coexpressed with the BPCV1 miRNA. Notably, the degree of negative regulation observed is equivalent to that of an identical reporter that is engineered to bind to the BPCV1 miRNA with perfect complementarity. We also show that this miRNA and this novel mode of early gene regulation are conserved with the related BPCV2. Finally, papillomatous lesions from a western barred bandicoot express readily detectable levels of this miRNA, stressing its likely importance in vivo. Combined, the alternative mechanisms of negative regulation of T antigen expression between the BPCVs and the polyomaviruses support the importance of miRNA-mediated autoregulation in the life cycles of some divergent polyomaviruses and polyomavirus-like viruses.

Novel DNA virus isolated from samples showing endothelial cell necrosis in the Japanese eel, Anguilla japonica

Economic loss due to viral endothelial cell necrosis of eel (VECNE) of Anguilla japonica is a serious problem for the cultured Japanese eel market. However, the viral genome responsible for VECNE is unknown. We recently developed a rapid determination system for viral nucleic acid sequences (RDV) to determine viral genome sequences. In this study, viral DNA fragments were obtained using RDV, and approximately 15-kbp circular full genome sequences were determined using a next-generation sequencing system, overlapping PCR, and Southern blot analysis. One open reading frame (ORF) was homologous to the large T-antigen of polyomavirus; other ORFs have no homology with any nucleic or amino acid sequences of polyomavirus. Therefore, as this DNA virus might comprise a novel virus family, we provisionally named it Japanese eel endothelial cells-infecting virus (JEECV). JEECV was detected in both naturally and experimentally infected eels, suggesting that JEECV potentially causes VECNE.