A novel virus detected in papillomas and carcinomas of the endangered western barred bandicoot (Perameles bougainville) exhibits genomic features of both the Papillomaviridae and Polyomaviridae

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.

Structural similarity of human papillomavirus E4 and polyomaviral VP4 exhibited by genomic analysis of the common kestrel (Falco tinnunculus) polyomavirus

Polyomaviruses are widely distributed viruses of birds that may induce developmental deformities and internal organ disorders primarily in nestlings. In this study, polyomavirus sequence was detected in kidney and liver samples of a common kestrel (Falco tinnunculus) that succumbed at a rescue station in Hungary. The amplified 5025 nucleotide (nt) long genome contained the early (large and small T antigen, LTA and STA) and late (viral proteins, VP1, VP2, VP3) open reading frames (ORFs) typical for polyomaviruses. One of the additional putative ORFs (named VP4) showed identical localization with the VP4 and ORF-X of gammapolyomaviruses, but putative splicing sites could not be found in its sequence. Interestingly, the predicted 123 amino acid (aa) long protein sequence showed the highest similarity with human papillomavirus E4 early proteins in respect of the aa distribution and motif arrangement implying similar functions. The LTA of the kestrel polyomavirus shared <59.2% nt and aa pairwise identity with the LTA sequence of other polyomaviruses and formed a separated branch in the phylogenetic tree among gammapolyomaviruses. Accordingly, the kestrel polyomavirus may be the first member of a novel species within the Gammapolyomavirus genus, tentatively named Gammapolyomavirus faltin.

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.

Return to 1616: Multispecies Fauna Reconstruction Requires Thinking Outside the Box

Conservation translocations have become increasingly popular for 'rewilding' areas that have lost their native fauna. These multispecies translocations are complex and need to consider the requirements of each individual species as well as the influence of likely interactions among them. The Dirk Hartog Island National Park Ecological Restoration Project, Return to 1616, aspires to restore ecological function to Western Australia's largest island. Since 2012, pest animals have been eradicated, and conservation translocations of seven fauna species have been undertaken, with a further six planned. Here, we present a synthesis of the innovative approaches undertaken in restoring the former faunal assemblage of Dirk Hartog Island and the key learnings gathered as the project has progressed.

Interpreting and de-noising genetically engineered barcodes in a DNA virus

The concept of a nucleic acid barcode applied to pathogen genomes is easy to grasp and the many possible uses are straightforward. But implementation may not be easy, especially when growing through multiple generations or assaying the pathogen long-term. The potential problems include: the barcode might alter fitness, the barcode may accumulate mutations, and construction of the marked pathogens may result in unintended barcodes that are not as designed. Here, we generate approximately 5,000 randomized barcodes in the genome of the prototypic small DNA virus murine polyomavirus. We describe the challenges faced with interpreting the barcode sequences obtained from the library. Our Illumina NextSeq sequencing recalled much greater variation in barcode sequencing reads than the expected 5,000 barcodes-necessarily stemming from the Illumina library processing and sequencing error. Using data from defined control virus genomes cloned into plasmid backbones we develop a vetted post-sequencing method to cluster the erroneous reads around the true virus genome barcodes. These findings may foreshadow problems with randomized barcodes in other microbial systems and provide a useful approach for future work utilizing nucleic acid barcoded pathogens.

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.

Next-Generation Sequencing for Confronting Virus Pandemics

Virus pandemics have happened, are happening and will happen again. In recent decades, the rate of zoonotic viral spillover into humans has accelerated, mirroring the expansion of our global footprint and travel network, including the expansion of viral vectors and the destruction of natural spaces, bringing humans closer to wild animals. Once viral cross-species transmission to humans occurs, transmission cannot be stopped by cement walls but by developing barriers based on knowledge that can prevent or reduce the effects of any pandemic. Controlling a local transmission affecting few individuals is more efficient that confronting a community outbreak in which infections cannot be traced. Genetic detection, identification, and characterization of infectious agents using next-generation sequencing (NGS) has been proven to be a powerful tool allowing for the development of fast PCR-based molecular assays, the rapid development of vaccines based on mRNA and DNA, the identification of outbreaks, transmission dynamics and spill-over events, the detection of new variants and treatment of vaccine resistance mutations, the development of direct-acting antiviral drugs, the discovery of relevant minority variants to improve knowledge of the viral life cycle, strengths and weaknesses, the potential for becoming dominant to take appropriate preventive measures, and the discovery of new routes of viral transmission.

Evolution of pathogen tolerance and emerging infections: A missing experimental paradigm

Researchers worldwide are repeatedly warning us against future zoonotic diseases resulting from humankind's insurgence into natural ecosystems. The same zoonotic pathogens that cause severe infections in a human host frequently fail to produce any disease outcome in their natural hosts. What precise features of the immune system enable natural reservoirs to carry these pathogens so efficiently? To understand these effects, we highlight the importance of tracing the evolutionary basis of pathogen tolerance in reservoir hosts, while drawing implications from their diverse physiological and life-history traits, and ecological contexts of host-pathogen interactions. Long-term co-evolution might allow reservoir hosts to modulate immunity and evolve tolerance to zoonotic pathogens, increasing their circulation and infectious period. Such processes can also create a genetically diverse pathogen pool by allowing more mutations and genetic exchanges between circulating strains, thereby harboring rare alive-on-arrival variants with extended infectivity to new hosts (i.e., spillover). Finally, we end by underscoring the indispensability of a large multidisciplinary empirical framework to explore the proposed link between evolved tolerance, pathogen prevalence, and spillover in the wild.

A novel lineage of polyomaviruses identified in bark scorpions

Polyomaviruses are non-enveloped viruses with circular double-stranded DNA genomes (~4-7 kb). Initially identified in mammals, polyomaviruses have now been identified in birds and a few fish species. Although fragmentary polyomavirus-like sequences have been detected as apparent 'hitchhikers' in shotgun genomics datasets of various arthropods, the possible diversity of these viruses in invertebrates remains unclear. Scorpions are predatory arachnids that are among the oldest terrestrial animals. Using high-throughput sequencing and traditional molecular techniques we determine the genome sequences of eight novel polyomaviruses in scorpions (Centruroides sculpturatus) from the greater Phoenix area, Arizona, USA. Analysis of Centruroides transcriptomic datasets elucidated the splicing of the viral late gene array, which is more complex than that of vertebrate polyomaviruses. Phylogenetic analysis provides further evidence of co-divergence of polyomaviruses with their hosts, suggesting that at least one ancestral species of polyomaviruses was circulating amongst the primitive common ancestors of arthropods and chordates.

Adintoviruses: a proposed animal-tropic family of midsize eukaryotic linear dsDNA (MELD) viruses

Polintons (also known as Mavericks) were initially identified as a widespread class of eukaryotic transposons named for their hallmark type B DNA polymerase and retrovirus-like integrase genes. It has since been recognized that many polintons encode possible capsid proteins and viral genome-packaging ATPases similar to those of a diverse range of double-stranded DNA viruses. This supports the inference that at least some polintons are actually viruses capable of cell-to-cell spread. At present, there are no polinton-associated capsid protein genes annotated in public sequence databases. To rectify this deficiency, we used a data-mining approach to investigate the distribution and gene content of polinton-like elements and related DNA viruses in animal genomic and metagenomic sequence datasets. The results define a discrete family-like clade of viruses with two genus-level divisions. We propose the family name Adintoviridae, connoting similarities to adenovirus virion proteins and the presence of a retrovirus-like integrase gene. Although adintovirus-class PolB sequences were detected in datasets for fungi and various unicellular eukaryotes, sequences resembling adintovirus virion proteins and accessory genes appear to be restricted to animals. Degraded adintovirus sequences are endogenized into the germlines of a wide range of animals, including humans.

Meta-transcriptomic identification of Trypanosoma spp. in native wildlife species from Australia

Background

Wildlife species carry a remarkable diversity of trypanosomes. The detection of trypanosome infection in native Australian fauna is central to understanding their diversity and host-parasite associations. The implementation of total RNA sequencing (meta-transcriptomics) in trypanosome surveillance and diagnosis provides a powerful methodological approach to better understand the host species distribution of this important group of parasites.

Methods

We implemented a meta-transcriptomic approach to detect trypanosomes in a variety of tissues (brain, liver, lung, skin, gonads) sampled from native Australian wildlife, comprising four marsupials (koala, Phascolarctos cinereus; southern brown bandicoot, Isoodon obesulus; swamp wallaby, Wallabia bicolor; bare-nosed wombat, Vombatus ursinus), one bird (regent honeyeater, Anthochaera phrygia) and one amphibian (eastern dwarf tree frog, Litoria fallax). Samples corresponded to both clinically healthy and diseased individuals. Sequencing reads were de novo assembled into contigs and annotated. The evolutionary relationships among the trypanosomatid sequences identified were determined through phylogenetic analysis of 18S rRNA sequences.

Results

We detected trypanosome sequences in all six species of vertebrates sampled, with positive samples in multiple organs and tissues confirmed by PCR. Phylogenetic analysis indicated that the trypanosomes infecting marsupials were related to those previously detected in placental and marsupial mammals, while the trypanosome in the regent honeyeater grouped with avian trypanosomes. In contrast, we provide the first evidence for a trypanosome in the eastern dwarf tree frog that was phylogenetically distinct from those described in other amphibians.

Conclusions

To our knowledge, this is the first meta-transcriptomic analysis of trypanosomes in native Australian wildlife, expanding the known genetic diversity of these important parasites. We demonstrated that RNA sequencing is sufficiently sensitive to detect low numbers of Trypanosoma transcripts and from diverse hosts and tissues types, thereby representing an effective means to detect trypanosomes that are divergent in genome sequence.

Papillomaviruses infecting cetaceans exhibit signs of genome adaptation following a recombination event

Papillomaviruses (PVs) have evolved through a complex evolutionary scenario where virus-host co-evolution alone is not enough to explain the phenotypic and genotypic PV diversity observed today. Other evolutionary processes, such as host switch and recombination, also appear to play an important role in PV evolution. In this study, we have examined the genomic impact of a recombination event between distantly related PVs infecting Cetartiodactyla (even-toed ungulates and cetaceans). Our phylogenetic analyses suggest that one single recombination was responsible for the generation of extant 'chimeric' PV genomes infecting cetaceans. By correlating the phylogenetic relationships to the genomic content, we observed important differences between the recombinant and non-recombinant cetartiodactyle PV genomes. Notably, recombinant PVs contain a unique set of conserved motifs in the upstream regulatory region (URR). We interpret these regulatory changes as an adaptive response to drastic changes in the PV genome. In terms of codon usage preferences (CUPrefs), we did not detect any particular differences between orthologous open reading frames in recombinant and non-recombinant PVs. Instead, our results are in line with previous observations suggesting that CUPrefs in PVs are rather linked to gene expression patterns as well as to gene function. We show that the non-coding URR of PVs infecting cetaceans, the central regulatory element in these viruses, exhibits signs of adaptation following a recombination event. Our results suggest that also in PVs, the evolution of gene regulation can play an important role in speciation and adaptation to novel environments.

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.

Genome Plasticity in Papillomaviruses and De Novo Emergence of E5 Oncogenes

The clinical presentations of papillomavirus (PV) infections come in many different flavors. While most PVs are part of a healthy skin microbiota and are not associated to physical lesions, other PVs cause benign lesions, and only a handful of PVs are associated to malignant transformations linked to the specific activities of the E5, E6, and E7 oncogenes. The functions and origin of E5 remain to be elucidated. These E5 open reading frames (ORFs) are present in the genomes of a few polyphyletic PV lineages, located between the early and the late viral gene cassettes. We have computationally assessed whether these E5 ORFs have a common origin and whether they display the properties of a genuine gene. Our results suggest that during the evolution of Papillomaviridae, at least four events lead to the presence of a long noncoding DNA stretch between the E2 and the L2 genes. In three of these events, the novel regions evolved coding capacity, becoming the extant E5 ORFs. We then focused on the evolution of the E5 genes in AlphaPVs infecting primates. The sharp match between the type of E5 protein encoded in AlphaPVs and the infection phenotype (cutaneous warts, genital warts, or anogenital cancers) supports the role of E5 in the differential oncogenic potential of these PVs. In our analyses, the best-supported scenario is that the five types of extant E5 proteins within the AlphaPV genomes may not have a common ancestor. However, the chemical similarities between E5s regarding amino acid composition prevent us from confidently rejecting the model of a common origin. Our evolutionary interpretation is that an originally noncoding region entered the genome of the ancestral AlphaPVs. This genetic novelty allowed to explore novel transcription potential, triggering an adaptive radiation that yielded three main viral lineages encoding for different E5 proteins, displaying distinct infection phenotypes. Overall, our results provide an evolutionary scenario for the de novo emergence of viral genes and illustrate the impact of such genotypic novelty in the phenotypic diversity of the viral infections.

Cancer in wildlife: patterns of emergence

Cancer is ubiquitous in wildlife, affecting animals from bivalves to pachyderms and cetaceans. Reports of increasing frequency demonstrate that neoplasia is associated with substantial mortality in wildlife species. Anthropogenic activities and global weather changes are shaping new geographical limitations for many species, and alterations in living niches are associated with visible examples of genetic bottlenecks, toxin exposures, oncogenic pathogens, stress and immunosuppression, which can all contribute to cancers in wild species. Nations that devote resources to monitoring the health of wildlife often do so for human-centric reasons, including for the prediction of the potential for zoonotic disease, shared contaminants, chemicals and medications, and for observing the effect of exposure from crowding and loss of habitat. Given the increasing human footprint on land and in the sea, wildlife conservation should also become a more important motivating factor. Greater attention to the patterns of the emergence of wildlife cancer is imperative because growing numbers of species are existing at the interface between humans and the environment, making wildlife sentinels for both animal and human health. Therefore, monitoring wildlife cancers could offer interesting and novel insights into potentially unique non-age-related mechanisms of carcinogenesis across species.

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.

Transcriptome sequencing of the long-nosed bandicoot (Perameles nasuta) reveals conservation and innovation of immune genes in the marsupial order Peramelemorphia

Bandicoots are omnivorous marsupials of the order Peramelemorphia. Conservation concerns and their unique biological characteristics suggest peramelomorphs are worthy research subjects, but knowledge of their genetics and immunology has lagged behind that of other high-profile marsupials. Here, we characterise the transcriptome of the long-nose bandicoot (Perameles nasuta), the first high-throughput data set from any peramelomorph. We investigate the immune gene repertoire of the bandicoot, with a focus on key immune gene families, and compare to previously characterised marsupial and mammalian species. We find that the immune gene complement in bandicoot is often conserved with respect to other marsupials; however, the diversity of expressed transcripts in several key families, such as major histocompatibility complex, T cell receptor μ and natural killer cell receptors, appears greater in the bandicoot than other Australian marsupials, including devil and koala. This transcriptome is an important first step for future studies of bandicoots and the bilby, allowing for population level analysis and construction of bandicoot-specific immunological reagents and assays. Such studies will be critical to understanding the immunology and physiology of Peramelemorphia and to inform the conservation of these unique marsupials.

Metastatic Squamous Cell Carcinoma in a Northern Brown Bandicoot (Isoodon macrourus)

Aside from a handful of notable exceptions, neoplasia is not reported as a major cause of mortality in wild animal populations and often goes undetected. For northern brown bandicoots specifically, there are few reported tumors in the literature and on file in the Australian Registry of Wildlife Health. This report describes a case of squamous cell carcinoma in a northern brown bandicoot (Isoodon macrourus), with metastases to the draining lymph nodes and lung. This neoplasm consisted predominantly of well-differentiated squamous cells and multifocal keratin pearls, with areas possibly consistent with epithelial to mesenchymal transition, as identified by positive immunohistochemical staining by both pancytokeratin (AE1/AE3) and vimentin. Additional investigations were negative for bandicoot papillomatosis carcinomatosis viruses.

A Bat-Derived Putative Cross-Family Recombinant Coronavirus with a Reovirus Gene

The emergence of severe acute respiratory syndrome coronavirus (SARS-CoV) in 2002 and Middle East respiratory syndrome coronavirus (MERS-CoV) in 2012 has generated enormous interest in the biodiversity, genomics and cross-species transmission potential of coronaviruses, especially those from bats, the second most speciose order of mammals. Herein, we identified a novel coronavirus, provisionally designated Rousettus bat coronavirus GCCDC1 (Ro-BatCoV GCCDC1), in the rectal swab samples of Rousettus leschenaulti bats by using pan-coronavirus RT-PCR and next-generation sequencing. Although the virus is similar to Rousettus bat coronavirus HKU9 (Ro-BatCoV HKU9) in genome characteristics, it is sufficiently distinct to be classified as a new species according to the criteria defined by the International Committee of Taxonomy of Viruses (ICTV). More striking was that Ro-BatCoV GCCDC1 contained a unique gene integrated into the 3'-end of the genome that has no homologs in any known coronavirus, but which sequence and phylogeny analyses indicated most likely originated from the p10 gene of a bat orthoreovirus. Subgenomic mRNA and cellular-level observations demonstrated that the p10 gene is functional and induces the formation of cell syncytia. Therefore, here we report a putative heterologous inter-family recombination event between a single-stranded, positive-sense RNA virus and a double-stranded segmented RNA virus, providing insights into the fundamental mechanisms of viral evolution.

Characterization of novel human papillomavirus types 157, 158 and 205 from healthy skin and recombination analysis in genus γ-Papillomavirus

Gammapapillomavirus (γ-PV) is a diverse and rapidly expanding genus, currently consisting of 79 fully characterized human PV (HPV) types. In this study, three novel types, HPV157, HPV158 and HPV205, obtained from healthy sun-exposed skin of two immunocompetent individuals, were amplified by the "Hanging droplet" long PCR technique, cloned, sequenced and characterized. HPV157, HPV158 and HPV205 genomes comprise 7154-bp, 7192-bp and 7298-bp, respectively, and contain four early (E1, E2, E6 and E7) and two late genes (L1 and L2). Phylogenetic analysis of the L1 ORF placed all novel types within the γ-PV genus: HPV157 was classified as a new member of species γ-12 while HPV158 and HPV205 belong to species γ-1. We then explored potential recombination events in genus γ-PV with the RDP4 program in a dataset of 74 viruses (71 HPV types with available full-length genomes and the 3 novel types). Two events, both located in the E1 ORF, met the inclusion criterion (p-values <0.05 with at least four methods) and persisted in different ORF combinations: an inter-species recombination in species γ-8 (major and minor parents: species γ-24 and γ-11, respectively), and an intra-species recombination in species γ-7 (recombinant strain: HPV170; major and minor parents: HPV-109 and HPV-149, respectively). These findings were confirmed by phylogenetic tree incongruence analysis. An additional incongruence was found in members of species γ-9 but it was not detected by the RDP4. This report expands our knowledge of the family Papillomaviridae and provides for the first time in silico evidence of recombination in genus γ-PV.

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 taxonomy update for the family Polyomaviridae

Many distinct polyomaviruses infecting a variety of vertebrate hosts have recently been discovered, and their complete genome sequence could often be determined. To accommodate this fast-growing diversity, the International Committee on Taxonomy of Viruses (ICTV) Polyomaviridae Study Group designed a host- and sequence-based rationale for an updated taxonomy of the family Polyomaviridae. Applying this resulted in numerous recommendations of taxonomical revisions, which were accepted by the Executive Committee of the ICTV in December 2015. New criteria for definition and creation of polyomavirus species were established that were based on the observed distance between large T antigen coding sequences. Four genera (Alpha-, Beta, Gamma- and Deltapolyomavirus) were delineated that together include 73 species. Species naming was made as systematic as possible - most species names now consist of the binomial name of the host species followed by polyomavirus and a number reflecting the order of discovery. It is hoped that this important update of the family taxonomy will serve as a stable basis for future taxonomical developments.

Assessing parallel gene histories in viral genomes

BACKGROUND

The increasing abundance of sequence data has exacerbated a long known problem: gene trees and species trees for the same terminal taxa are often incongruent. Indeed, genes within a genome have not all followed the same evolutionary path due to events such as incomplete lineage sorting, horizontal gene transfer, gene duplication and deletion, or recombination. Considering conflicts between gene trees as an obstacle, numerous methods have been developed to deal with these incongruences and to reconstruct consensus evolutionary histories of species despite the heterogeneity in the history of their genes. However, inconsistencies can also be seen as a source of information about the specific evolutionary processes that have shaped genomes.

RESULTS

The goal of the approach here proposed is to exploit this conflicting information: we have compiled eleven variables describing phylogenetic relationships and evolutionary pressures and submitted them to dimensionality reduction techniques to identify genes with similar evolutionary histories. To illustrate the applicability of the method, we have chosen two viral datasets, namely papillomaviruses and Turnip mosaic virus (TuMV) isolates, largely dissimilar in genome, evolutionary distance and biology. Our method pinpoints viral genes with common evolutionary patterns. In the case of papillomaviruses, gene clusters match well our knowledge on viral biology and life cycle, illustrating the potential of our approach. For the less known TuMV, our results trigger new hypotheses about viral evolution and gene interaction.

CONCLUSIONS

The approach here presented allows turning phylogenetic inconsistencies into evolutionary information, detecting gene assemblies with similar histories, and could be a powerful tool for comparative pathogenomics.

A Comprehensive Analysis of Replicating Merkel Cell Polyomavirus Genomes Delineates the Viral Transcription Program and Suggests a Role for mcv-miR-M1 in Episomal Persistence

Merkel cell polyomavirus (MCPyV) is considered the etiological agent of Merkel cell carcinoma and persists asymptomatically in the majority of its healthy hosts. Largely due to the lack of appropriate model systems, the mechanisms of viral replication and MCPyV persistence remain poorly understood. Using a semi-permissive replication system, we here report a comprehensive analysis of the role of the MCPyV-encoded microRNA (miRNA) mcv-miR-M1 during short and long-term replication of authentic MCPyV episomes. We demonstrate that cells harboring intact episomes express high levels of the viral miRNA, and that expression of mcv-miR-M1 limits DNA replication. Furthermore, we present RACE, RNA-seq and ChIP-seq studies which allow insight in the viral transcription program and mechanisms of miRNA expression. While our data suggest that mcv-miR-M1 can be expressed from canonical late strand transcripts, we also present evidence for the existence of an independent miRNA promoter that is embedded within early strand coding sequences. We also report that MCPyV genomes can establish episomal persistence in a small number of cells for several months, a time period during which viral DNA as well as LT-Ag and viral miRNA expression can be detected via western blotting, FISH, qPCR and southern blot analyses. Strikingly, despite enhanced replication in short term DNA replication assays, a mutant unable to express the viral miRNA was severely limited in its ability to establish long-term persistence. Our data suggest that MCPyV may have evolved strategies to enter a non- or low level vegetative stage of infection which could aid the virus in establishing and maintaining a lifelong persistence.

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.

Identification of an avian polyomavirus associated with Adélie penguins (Pygoscelis adeliae)

Little is known about viruses associated with Antarctic animals, although they are probably widespread. We recovered a novel polyomavirus from Adélie penguin (Pygoscelis adeliae) faecal matter sampled in a subcolony at Cape Royds, Ross Island, Antarctica. The 4988 nt Adélie penguin polyomavirus (AdPyV) has a typical polyomavirus genome organization with three ORFs that encoded capsid proteins on the one strand and two non-structural protein-coding ORFs on the complementary strand. The genome of AdPyV shared ~60 % pairwise identity with all avipolyomaviruses. Maximum-likelihood phylogenetic analysis of the large T-antigen (T-Ag) amino acid sequences showed that the T-Ag of AdPyV clustered with those of avipolyomaviruses, sharing between 48 and 52 % identities. Only three viruses associated with Adélie penguins have been identified at a genomic level, avian influenza virus subtype H11N2 from the Antarctic Peninsula and, respectively, Pygoscelis adeliae papillomavirus and AdPyV from capes Crozier and Royds on Ross Island.

Novel papillomaviruses in free-ranging Iberian bats: no virus-host co-evolution, no strict host specificity, and hints for recombination

Papillomaviruses (PVs) are widespread pathogens. However, the extent of PV infections in bats remains largely unknown. This work represents the first comprehensive study of PVs in Iberian bats. We identified four novel PVs in the mucosa of free-ranging Eptesicus serotinus (EserPV1, EserPV2, and EserPV3) and Rhinolophus ferrumequinum (RferPV1) individuals and analyzed their phylogenetic relationships within the viral family. We further assessed their prevalence in different populations of E. serotinus and its close relative E. isabellinus. Although it is frequent to read that PVs co-evolve with their host, that PVs are highly species-specific, and that PVs do not usually recombine, our results suggest otherwise. First, strict virus–host co-evolution is rejected by the existence of five, distantly related bat PV lineages and by the lack of congruence between bats and bat PVs phylogenies. Second, the ability of EserPV2 and EserPV3 to infect two different bat species (E. serotinus and E. isabellinus) argues against strict host specificity. Finally, the description of a second noncoding region in the RferPV1 genome reinforces the view of an increased susceptibility to recombination in the E2-L2 genomic region. These findings prompt the question of whether the prevailing paradigms regarding PVs evolution should be reconsidered.

A novel papillomavirus in Adélie penguin (Pygoscelis adeliae) faeces sampled at the Cape Crozier colony, Antarctica

Papillomaviruses are epitheliotropic viruses that have circular dsDNA genomes encapsidated in non-enveloped virions. They have been found to infect a variety of mammals, reptiles and birds, but so far they have not been found in amphibians. Using a next-generation sequencing de novo assembly contig-informed recovery, we cloned and Sanger sequenced the complete genome of a novel papillomavirus from the faecal matter of Adélie penguins (Pygoscelis adeliae) nesting on Ross Island, Antarctica. The genome had all the usual features of a papillomavirus and an E9 ORF encoding a protein of unknown function that is found in all avian papillomaviruses to date. This novel papillomavirus genome shared ~60 % pairwise identity with the genomes of the other three known avian papillomaviruses: Fringilla coelebs papillomavirus 1 (FcPV1), Francolinus leucoscepus papillomavirus 1 (FlPV1) and Psittacus erithacus papillomavirus 1. Pairwise identity analysis and phylogenetic analysis of the major capsid protein gene clearly indicated that it represents a novel species, which we named Pygoscelis adeliae papillomavirus 1 (PaCV1). No evidence of recombination was detected in the genome of PaCV1, but we did detect a recombinant region (119 nt) in the E6 gene of FlPV1 with the recombinant region being derived from ancestral FcPV1-like sequences. Previously only paramyxoviruses, orthomyxoviruses and avian pox viruses have been genetically identified in penguins; however, the majority of penguin viral identifications have been based on serology or histology. This is the first report, to our knowledge, of a papillomavirus associated with a penguin species.

The miRNA world of polyomaviruses

Polyomaviruses are a family of non-enveloped DNA viruses infecting several species, including humans, primates, birds, rodents, bats, horse, cattle, raccoon and sea lion. They typically cause asymptomatic infection and establish latency but can be reactivated under certain conditions causing severe diseases. MicroRNAs (miRNAs) are small non-coding RNAs that play important roles in several cellular processes by binding to and inhibiting the translation of specific mRNA transcripts. In this review, we summarize the current knowledge of microRNAs involved in polyomavirus infection. We review in detail the different viral miRNAs that have been discovered and the role they play in controlling both host and viral protein expression. We also give an overview of the current understanding on how host miRNAs may function in controlling polyomavirus replication, immune evasion and pathogenesis.

Evolution of the papillomaviridae

Viruses belonging to the Papillomaviridae family have been isolated from a variety of mammals, birds and non-avian reptiles. It is likely that most, if not all, amniotes carry a broad array of viral types. To date, the complete genomic sequence of more than 240 distinct viral types has been characterized at the nucleotide level. The analysis of this sequence information has begun to shed light on the evolutionary history of this important virus family. The available data suggests that many different evolutionary mechanisms have influenced the papillomavirus phylogenetic tree. Increasing evidence supports that the ancestral papillomavirus initially specialized to infect different ecological niches on the host. This episode of niche sorting was followed by extensive episodes of co-speciation with the host. This review attempts to summarize our current understanding of the papillomavirus evolution.

A cornucopia of human polyomaviruses

During the past 6 years, focused virus hunting has led to the discovery of nine new human polyomaviruses, including Merkel cell polyomavirus, which has been linked to Merkel cell carcinoma, a lethal skin cell cancer. The discovery of so many new and highly divergent human polyomaviruses raises key questions regarding their evolution, tropism, latency, reactivation, immune evasion and contribution to disease. This Review describes the similarities and differences among the new human polyomaviruses and discusses how these viruses might interact with their human host.

Discovery of STL polyomavirus, a polyomavirus of ancestral recombinant origin that encodes a unique T antigen by alternative splicing

The family Polyomaviridae is comprised of circular double-stranded DNA viruses, several of which are associated with diseases, including cancer, in immunocompromised patients. Here we describe a novel polyomavirus recovered from the fecal microbiota of a child in Malawi, provisionally named STL polyomavirus (STLPyV). We detected STLPyV in clinical stool specimens from USA and The Gambia at up to 1% frequency. Complete genome comparisons of two STLPyV strains demonstrated 5.2% nucleotide divergence. Alternative splicing of the STLPyV early region yielded a unique form of T antigen, which we named 229T, in addition to the expected large and small T antigens. STLPyV has a mosaic genome and shares an ancestral recombinant origin with MWPyV. The discovery of STLPyV highlights a novel alternative splicing strategy and advances our understanding of the complex evolutionary history of polyomaviruses.

Evidence for an endogenous papillomavirus-like element in the platypus genome

Papillomaviruses (PVs) infect a wide range of vertebrates and have diversified into multiple genetic types, some of which have serious consequences for human health. Although PVs have to date only been characterized as exogenous viral forms, here we report the observation of an endogenous viral element (EPVLoa) in the genome of the platypus (Ornithorhynchus anatinus) that is related to PVs. Further data mining for endogenous PV-like elements is therefore warranted.

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.

Evolution and taxonomic classification of human papillomavirus 16 (HPV16)-related variant genomes: HPV31, HPV33, HPV35, HPV52, HPV58 and HPV67

Background

Human papillomavirus 16 (HPV16) species group (alpha-9) of the Alphapapillomavirus genus contains HPV16, HPV31, HPV33, HPV35, HPV52, HPV58 and HPV67. These HPVs account for 75% of invasive cervical cancers worldwide. Viral variants of these HPVs differ in evolutionary history and pathogenicity. Moreover, a comprehensive nomenclature system for HPV variants is lacking, limiting comparisons between studies.

Methods

DNA from cervical samples previously characterized for HPV type were obtained from multiple geographic regions to screen for novel variants. The complete 8 kb genomes of 120 variants representing the major and minor lineages of the HPV16-related alpha-9 HPV types were sequenced to capture maximum viral heterogeneity. Viral evolution was characterized by constructing phylogenic trees based on complete genomes using multiple algorithms. Maximal and viral region specific divergence was calculated by global and pairwise alignments. Variant lineages were classified and named using an alphanumeric system; the prototype genome was assigned to the A lineage for all types.

Results

The range of genome-genome sequence heterogeneity varied from 0.6% for HPV35 to 2.2% for HPV52 and included 1.4% for HPV31, 1.1% for HPV33, 1.7% for HPV58 and 1.1% for HPV67. 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. Each gene/region differs in sequence diversity, from most variable to least variable: noncoding region 1 (NCR1) /noncoding region 2 (NCR2) >upstream regulatory region (URR)> E6/E7 > E2/L2 > E1/L1.

Conclusions

These data define maximum viral genomic heterogeneity of HPV16-related alpha-9 HPV variants. The proposed nomenclature system facilitates the comparison of variants across epidemiological studies. Sequence diversity and phylogenies of this clinically important group of HPVs provides the basis for further studies of discrete viral evolution, epidemiology, pathogenesis and preventative/therapeutic interventions.

Taxonomical developments in the family Polyomaviridae

The Polyomaviridae Study Group of the International Committee on Taxonomy of Viruses (ICTV) has recommended several taxonomical revisions, as follows: The family Polyomaviridae, which is currently constituted as a single genus (Polyomavirus), will be comprised of three genera: two containing mammalian viruses and one containing avian viruses. The two mammalian genera will be designated Orthopolyomavirus and Wukipolyomavirus, and the avian genus will be named Avipolyomavirus. These genera will be created by the redistribution of species from the current single genus (Polyomavirus) and by the inclusion of several new species. In addition, the names of several species will be changed to reflect current usage.

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.

Merkel cell polyomavirus and two previously unknown polyomaviruses are chronically shed from human skin

Mounting evidence indicates that Merkel cell polyomavirus (MCV), a circular double-stranded DNA virus, is a causal factor underlying a highly lethal form of skin cancer known as Merkel cell carcinoma. To explore the possibility that MCV and other polyomaviruses commonly inhabit healthy human skin, we developed an improved rolling circle amplification (RCA) technique to isolate circular DNA viral genomes from human skin swabs. Complete MCV genomes were recovered from 40% of healthy adult volunteers tested, providing full-length, apparently wild-type cloned MCV genomes. RCA analysis also identified two previously unknown polyomavirus species that we name human polyomavirus-6 (HPyV6) and HPyV7. Biochemical experiments show that polyomavirus DNA is shed from the skin in the form of assembled virions. A pilot serological study indicates that infection or coinfection with these three skin-tropic polyomaviruses is very common. Thus, at least three polyomavirus species are constituents of the human skin microbiome.

Multiply primed rolling-circle amplification method for the amplification of circular DNA viruses

The use of whole genome amplification and analysis of viruses is of increasing importance, as data generated using these methods are currently used for clinical diagnostics, epidemiological studies, phylogenetic analyses, and studies of genome organization and evolution. The best known amplification method for DNA is the polymerase chain reaction (PCR). This technique, however, has drawbacks: PCR produces relatively small amplicons and also requires prior knowledge of sequence data for the construction of the required consensus or degenerate primers. For circular DNA templates, it is possible to overcome these drawbacks by using the multiply primed rolling-circle amplification (RCA) technique, which mimics the rolling-circle mechanism that occurs in nature for replication of circular DNA molecules, e.g., plasmids. The RCA protocol described here is optimized for the amplification of circular DNA virus genomes.

Classification of papillomaviruses (PVs) based on 189 PV types and proposal of taxonomic amendments

We present an expansion of the classification of the family Papillomaviridae, which now contains 29 genera formed by 189 papillomavirus (PV) types isolated from humans (120 types), non-human mammals, birds and reptiles (64, 3 and 2 types, respectively). To accommodate the number of PV genera exceeding the Greek alphabet, the prefix "dyo" is used, continuing after the Omega-PVs with Dyodelta-PVs. The current set of human PVs is contained within five genera, whereas mammalian, avian and reptile PVs are contained within 20, 3 and 1 genera, respectively. We propose standardizations to the names of a number of animal PVs. As prerequisite for a coherent nomenclature of animal PVs, we propose founding a reference center for animal PVs. We discuss that based on emerging species concepts derived from genome sequences, PV types could be promoted to the taxonomic level of species, but we do not recommend implementing this change at the current time.

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

The first fully sequenced papillomavirus (PV) of marsupials, tentatively named Bettongia penicillata papillomavirus type 1 (BpPV1), was detected in papillomas from a woylie (Bettongia penicillata ogilbyi). The circular, double-stranded DNA genome contains 7,737 bp and encodes 7 open reading frames (ORFs), E6, E7, E1, E2, E4, L2, and L1, in typical PV conformation. BpPV1 is a close-to-root PV with L1 and L2 ORFs most similar to European hedgehog PV and bandicoot papillomatosis carcinomatosis virus types 1 and 2 (BPCV1 and -2). It appears that the BPCVs arose by recombination between an ancient PV and an ancient polyomavirus more than 10 million years ago.

Papillomavirus-Associated Cutaneous Neoplasia in Mammals

For humans, strong evidence indicates that some mucosal papillomavirus (PV) types cause genital and oral neoplasia, and weaker evidence suggests that some cutaneous PVs may cause cutaneous squamous cell carcinomas (SCC). For nonhuman species, strong evidence supports a causal role of PVs in the development of feline and equine sarcoids. Likewise, PVs are believed to cause cutaneous SCCs in rabbits, western barred bandicoots, and some rodents. Furthermore, some evidence suggests that PVs may influence the development of both feline and canine cutaneous SCCs. This review discusses the evidence that PVs cause human cutaneous SCCs and the proposed mechanisms for this action. It then reviews preneoplastic and neoplastic skin diseases that are associated with PV infection in nonhuman mammals.

Prevalence, emergence, and factors associated with a viral papillomatosis and carcinomatosis syndrome in wild, reintroduced, and captive western barred bandicoots (Perameles bougainville)

Once widespread across western and southern Australia, wild populations of the western barred bandicoot (WBB) are now only found on Bernier and Dorre Islands, Western Australia. Conservation efforts to prevent the extinction of the WBB are presently hampered by a papillomatosis and carcinomatosis syndrome identified in captive and wild bandicoots, associated with infection with the bandicoot papillomatosis carcinomatosis virus type 1 (BPCV1). This study examined the prevalence and distribution of BPCV1 and the associated syndrome in two island and four mainland (reintroduced and captive) WBB populations in Western Australia, and factors that may be associated with susceptibility to this syndrome. BPCV1 and the syndrome were found in the wild WBB population at Red Cliff on Bernier Island, and in mainland populations established from all or a proportion of founder WBBs from Red Cliff. BPCV1 and the syndrome were not found in the wild population on Dorre Island or in the mainland population founded by animals exclusively from Dorre Island. Findings suggested that BPCV1 and the syndrome were disseminated into mainland WBB populations through the introduction of affected WBBs from Red Cliff. No difference in susceptibility to the syndrome was found between Dorre Island, Bernier Island, and island-cross individuals. Severity of lesions and the number of affected animals observed in captivity was greater than that observed in wild populations. This study provided epidemiological evidence to support the pathological and molecular association between BPCV1 infection and the papillomatosis and carcinomatosis syndrome and revealed increasing age as an additional risk factor for this disease.