Next Article in Journal
Identification of the p34 Protein of African Swine Fever Virus as a Novel Viral Antigen with Protection Potential
Next Article in Special Issue
Canis Familiaris Papillomavirus Type 26: A Novel Papillomavirus of Dogs and the First Canine Papillomavirus within the Omegapapillomavirus Genus
Previous Article in Journal
Susceptibility and Resistance of SARS-CoV-2 Variants to LCB1 and Its Multivalent Derivatives
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Viruses
Volume 16
Issue 1
10.3390/v16010037
viruses-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Brief Report

Interspecies Papillomavirus Type Infection and a Novel Papillomavirus Type in Red Ruffed Lemurs (Varecia rubra)

by
Elise N. Paietta
1,*,
Simona Kraberger
2,
Melanie Regney
2,
Joy M. Custer
2,
Erin Ehmke
3,
Anne D. Yoder
1 and
Arvind Varsani
2,4,*
1
Department of Biology, Duke University, Durham, NC 27708, USA
2
The Biodesign Center for Fundamental and Applied Microbiomics, Center for Evolution and Medicine, School of Life Sciences, Arizona State University, Tempe, AZ 85287, USA
3
Duke Lemur Center, Durham, NC 27705, USA
4
Structural Biology Research Unit, Department of Integrative Biomedical Sciences, University of Cape Town, Cape Town 7925, South Africa
*
Authors to whom correspondence should be addressed.
Viruses 2024, 16(1), 37; https://doi.org/10.3390/v16010037
Submission received: 19 November 2023 / Revised: 16 December 2023 / Accepted: 20 December 2023 / Published: 25 December 2023
(This article belongs to the Special Issue Animal Papillomaviruses Research)
Browse Figures
Versions Notes

Abstract

:
The Papillomaviridae are a family of vertebrate-infecting viruses of oncogenic potential generally thought to be host species- and tissue-specific. Despite their phylogenetic relatedness to humans, there is a scarcity of data on papillomaviruses (PVs) in speciose non-human primate lineages, particularly the lemuriform primates. Varecia variegata (black-and-white ruffed lemurs) and Varecia rubra (red ruffed lemurs), two closely related species comprising the Varecia genus, are critically endangered with large global captive populations. Varecia variegata papillomavirus (VavPV) types −1 and −2, the first PVs in lemurs with a fully identified genome, were previously characterized from captive V. variegata saliva. To build upon this discovery, saliva samples were collected from captive V. rubra with the following aims: (1) to identify PVs shared between V. variegata and V. rubra and (2) to characterize novel PVs in V. rubra to better understand PV diversity in the lemuriform primates. Three complete PV genomes were determined from V. rubra samples. Two of these PV genomes share 98% L1 nucleotide identity with VavPV2, denoting interspecies infection of V. rubra by VavPV2. This work represents the first reported case of interspecies PV infection amongst the strepsirrhine primates. The third PV genome shares <68% L1 nucleotide identity with that of all PVs. Thus, it represents a new PV species and has been named Varecia rubra papillomavirus 1 (VarPV1). VavPV1, VavPV2, and VarPV1 form a new clade within the Papillomaviridae family, likely representing a novel genus. Future work diversifying sample collection (i.e., lemur host species from multiple genera, sample type, geographic location, and wild populations) is likely to uncover a world of diverse lemur PVs.

1. Introduction

Papillomaviruses (PVs) are double-stranded DNA viruses with icosahedral capsids that infect diverse vertebrates, including mammal, avian, fish, and reptile species [1,2]. Depending on the PV type, infected hosts can experience a variety of disease outcomes ranging from asymptomatic cases—the majority of PV infections—to benign papillomas and even progressing to invasive cancer [3,4,5,6]. PV genomes are ~6–8 kb in size and encode at least four proteins (i.e., two early proteins—E1, E2, and two late proteins—L1, L2) with high conservation of the L1 capsid protein across PVs [1,7]. PVs have co-evolved with their hosts and are generally recognized as host species- and tissue-specific [4,8]. However, the assumption that PVs are primarily host species-specific is being increasingly challenged by the detection of several PV types in species other than their host species [9]. While host–pathogen co-evolution, intra-host duplication, adaptive radiation, and recombination have been found to be drivers of PV evolution, mounting evidence of interspecies infection suggests that host-switching could also serve as a crucial evolutionary force in the history of PVs [9,10,11,12,13,14].
Cross-species PV infections leading to tumor disease are described for several PVs, especially those in the genus Deltapapillomavirus, including bovine papillomavirus (BPV) types −1, −2, or −13 in equine sarcoids [15,16,17], BPV-1 in captive tapir sarcoids [18], BPV-2 in bladder tumors of water buffalo [19], BPV-1 in cutaneous and perivulvar fibropapillomas in water buffalo [20], and BPV-14 in feline sarcoids in domestic cats and captive African lions [21,22]. Cross-species experimental inoculation with BPV1 has been shown to induce fibromas and fibrosarcomas in hamsters and mice [23,24,25] and oral papillomas in domestic dogs from coyote oral papillomavirus [26]. Overall, BPVs, in particular, have been repeatedly observed to have interspecies hosts with pathogenic effects on new host species. This may in part be due to the tissue tropism of PVs in the genus Deltapapillomavirus, as they can infect dermal fibroblasts [27,28].
Beyond BPVs, naturally occurring interspecies PV infection has been seen across numerous mammalian orders (Artiodactyla, Carnivora, Chiroptera, Lagomorpha, Primates, and Rodentia) and avian orders (Anseriformes, Charadriiformes, and Passeriformes) (Figure 1) [9,10,12,13,16,18,19,20,21,22,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46]. Interspecies PV infection in non-human primates has only been detected in very closely phylogenetically related species in the same genus. Four alphapapillomavirus types, Macaca fascicularis papillomavirus (MfPV) types −1, −8, −11, and Macaca mulatta papillomavirus (MmPV) −1, have been detected across rhesus macaques (Macaca mulatta) and cynomolgus macaques (Macaca fascicularis). Additionally, Pan paniscus papillomavirus 1 (PpPV1) was isolated from both chimpanzees (Pan troglodytes) and bonobos (Pan paniscus), although it is unclear which species served as the first host of PpPV1 [35].
Only 16 out of >390 non-human primate species have been screened for PVs (Table 1). Acknowledging ongoing taxonomic adjustments across all non-human primate lineages, the New World monkeys (>100 species), Old World monkeys (>130 species), and lemurs (>100 species) are the most speciose non-human primate lineages [47,48,49,50]. For understanding PV diversity in these vast lineages, just 8 complete PV genomes across 5 New World monkey species are available on NCBI, while 23 complete PV genomes across 6 Old World monkey species are currently available (Table 1). Further, despite lemuriforms comprising ~20% of primate species, only three complete genomes of PVs in lemurs from just one host species had been characterized prior to this study.
We previously identified two PV types, Varecia variegata papillomavirus (VavPV) −1 and −2, from captive black-and-white ruffed lemur (Varecia variegata variegata) saliva samples, providing the first complete genomes of PVs in the lemuriform primates [51]. VavPV1 and VavPV2 share <64% L1 identity with one another and <66% L1 identity with all other PV L1 sequences. VavPV1 and −2 formed a distinct clade within the Firstpapillomavirinae sub-family and likely represent a novel genus [51].
As interspecies PV type infection in primates has been detected in sister species, we obtained saliva samples from red ruffed lemurs (Varecia rubra), the only other species within the Varecia genus. Our primary goals were (1) to identify PVs shared between V. variegata and V. rubra and (2) to identify any new PVs in V. rubra to better understand PV diversity in the lemuriform primates using viral metagenomics. Occurring in the eastern rainforests of Madagascar, V. variegata and V. rubra are diurnal, frugivorous lemurs vital to their ecosystems as seed dispersers and pollinators [47]. Both are classified as critically endangered and, thus, have extensive global captive populations [52,53]. The work presented here is relevant for the health of captive Varecia populations and for forming a baseline of our knowledge of lemur PVs from which future comparisons between captive and wild Varecia PVs will be possible.

2. Materials and Methods

To characterize additional novel lemur PVs and to identify potential interspecies infection by PV types, two saliva samples from captive red ruffed lemurs (Varecia rubra) were collected at the Duke Lemur Center in Durham (NC, USA). The sampled lemurs appeared healthy with no apparent symptoms and continued to be monitored by veterinary staff. V. rubra saliva samples were collected under IACUC #A161-21-08 in August–September 2022. Saliva was obtained by allowing the lemurs to chew on a SalivaBio Children’s Swab (Salivametrics, Carlsbad, CA, USA). Swabs saturated with saliva were placed within a SalivaBio Swab Storage Tube (Salivametrics, Carlsbad, CA, USA) and centrifuged to collect the saliva. Saliva samples were stored at −80 °C until viral DNA extraction.
Immediately prior to extraction, SM buffer (0.1 M NaCl, 50 mM Tris-HCl [pH 7.4]) was added to each saliva sample to obtain a final volume of 400 µL. Viral DNA was extracted from 200 µL of diluted sample using the High Pure Viral Nucleic Acid Kit (Roche Diagnostics, Indianapolis, IN, USA), and circular DNA in the extract was amplified using the Illustra TempliPhi Kit (GE Healthcare, Chicago, IL, USA). Illumina sequencing libraries were generated using the Illumina DNA Prep Kit (with Tagmentation) and sequenced on an Illumina HiSeq 2500 at Psomagen Inc. (Rockville, MD, USA).
Paired-end reads (2 × 150 bp) were trimmed using Trimmomatic-0.39 [54] and de novo assembled with MEGAHITv.1.2.9 [55]. After circular contigs were identified based on terminal redundancy, contigs >1000 nts were analyzed for viral-like sequences using Diamond [56] BLASTx against a local viral protein RefSeq database (release 220; downloaded September 2023). Potential PV-like contigs were confirmed using BLASTn [57]. Genomes were annotated using CenoteTaker2 [58] and refined with PaVE [59].
To determine the genera of the PV genomes characterized in this study, datasets of papillomavirus E1, E2, and L1 protein sequences were constructed using PaVE reference and non-reference sequences. The datasets were aligned with the Varecia PVs using MAFFT v7.113 [60] in AUTO mode. Alignments were trimmed with TrimAL [61] (0.2 gap threshold). ProtTest 3 [62] was used to determine the best-fit amino acid substitution models for each dataset. A partitioned maximum likelihood phylogenetic tree of concatenated E1 + E2 + L1 amino acid sequences was built using IQ-TREE 2 [63] with partition models LG + I + G for E1, LG + I + G + F for E2, and LG + I + G + F for L1. The tree was rooted with avian and reptilian PV sequences and edited in iTOL v6 [64].
Mitochondrial genomes were annotated using MITOS Web Server [65]. For mitochondrial genome comparisons, available mitochondrial genomes for lemur species within the Lemuridae family (and Indriidae as an outgroup) were aligned with mitochondrial genomes characterized in this study using MAFFT [60]. The mitochondrial genome maximum likelihood phylogenetic tree was built using IQ-TREE 2 with ModelFinder and ultrafast bootstrap (UFBoot) approximation (1000 bootstrap replicates) and edited in iTOL v6 [64].

3. Results

Our viral metagenomic workflows enabled the identification of three circular contigs, ranging in size from 7452 to 7770 nts in length, that represented complete PV genome sequences based on terminal redundancy. The mapped reads have been deposited at SRA under SRR26324874 and SRR26324875. The genome sequences are deposited in GenBank under accessions OR734654-OR734656. For OR734654 (7452 nts), the depth of coverage is 110× with 5463 reads, and for OR734655 (7770 nts) and OR734656 (7770 nts), the depth is 196× and 296× with 10134 and 15325 reads, respectively.
All three genomes contain open reading frames for L1, L2, E1, E2, E6, and E7 (Figure 2). As PV-type demarcation is determined by L1 sequence similarity <70%, we compared L1 nucleotide identity between the V. rubra-derived papillomavirus L1 sequences characterized in this study and the previously characterized VavPV1 and VavPV2 L1 sequences detected from V. variegata.
OR734655 and OR734656 share 100% L1 nucleotide identity and were isolated from twin V. rubra lemurs M1 and J1, female lemurs housed together at the time of sampling. Further, OR734655 and OR734656 share 98% L1 nucleotide identity with VavPV2. Therefore, these two PVs belong to the same PV type (i.e., VavPV2). As VavPV2 has been found in two V. rubra and two V. variegata female individuals’ saliva from the Duke Lemur Center, this serves as a case of interspecies PV-type infection between non-human primate sister species.
The third complete PV genome sequence identified in V. rubra, OR734654, shares <68% L1 nucleotide identity with VavPV1, VavPV2, and all other PVs. Based on the PV species demarcation threshold of 70%, this PV represents a novel type which we have named Varecia rubra papillomavirus (VarPV) type 1. VarPV1 is the third PV type and species to be characterized in the lemuriform primates. Furthermore, VarPV1 was identified in the individual M1 from which we also identified VavPV2, thus suggesting an oral co-infection.
In the E6 and E7 proteins of VarPV1 and VarPV2 from V. rubra, we identified the conserved zinc-binding domains (CxxC) (Figure 3). Unlike the pRB-binding motif (Lx[C/S]xE) identified in VavPV1 from V. variegata, we were unable to identify this motif in VarPV1 or VavPV2. We also identified the conserved regions 1 and 2 in VarPV1 and VavPV2 that are homologous to those in the E1A protein of human adenovirus 5 (family Adenoviridae) and the large tumor antigen of simian virus 40 (family Polyomaviridae) [66].
Based on their E1 + E2 + L1 protein sequence phylogeny, the two VavPVs and VarPV1 form a new cluster within a well-supported clade with sequences in the genera Dyoxipapillomavirus, Gammapapillomavirus, Pipapillomavirus, Taupapillomavirus, and Treisetapapillomavirus, in addition to six non-classified PV types, including BPV19, BPV21, and BPV27 (Figure 4). Overall, due to the divergence of the Varecia PV cluster from known PV genera, this new cluster likely represents at least one new genus.
The results presented in this study highlight the use of viral metagenomics for determining the complete genomes of viruses potentially relevant to endangered species’ health. In addition, metagenomic protocols employing rolling circle amplification are also ideal for the determination of mitochondrial genomes. Two complete mitochondrial genomes were characterized in V. rubra saliva from individuals M1 and J1, which are female twins. The raw reads are deposited under SRA accessions SRR26324872 and SRR26324873, and genomes are deposited in GenBank under accessions OR711366 and OR711367. Both mitochondrial genomes are 16972 nt in length and share 100% nucleotide identity. The cytochrome-b (cytb) genes, the gene primarily used for phylogenetic species identification in lemurs, in both mitochondrial genomes are 1140 nt in length; the same length has been previously found in V. rubra cytb (accession number AY441450). Although sequences of V. rubra mitochondrial genome cytb and D-loops [67] have been published, there were no complete V. rubra mitochondrial genomes available on NCBI prior to this study. The complete mitochondrial genomes of V. rubra share 96.8% nucleotide identity with V. variegata (accession numbers AB371089 and KJ944176), emphasizing the close evolutionary relationship between the two species (Figure 5).

4. Discussion

Ongoing habitat loss has dramatically impacted Varecia populations in Madagascar. Additionally, given Varecia’s vital roles in primate evolutionary history and ecosystem services, large captive populations of around 800 V. variegata and 600 V. rubra are maintained globally [53,68]. To maintain healthy captive Varecia populations, understanding the pathogens that impact them is imperative for the survival of these species. The Duke Lemur Center in Durham (NC, USA) houses the largest and most diverse population of lemurs outside Madagascar and is currently home to 13 V. variegata and 12 V. rubra individuals. The Duke Lemur Center plays an active role in the Species Survival Plans and Population Analysis and Breeding and Transfer Plans for both V. variegata and V. rubra. Thus, the study of viruses circulating in captive populations is essential for learning about viral diversity within Varecia and implementing strategies to reduce the burden of pathogenic viruses that can be detrimental to their health and conservation.
Although the majority of PV infections are asymptomatic, information about PVs is relevant for lemur health as some types have been associated with invasive cancer in humans and non-human primates. In our previous work, we determined the first complete genomes of PVs in lemurs, VavPV1 and VavPV2, in V. variegata saliva. In this study, we present a case of interspecies infection by VavPV2 in V. rubra. In addition, a new PV type has been identified in V. rubra and termed VarPV1. In non-human primates, the two previously known instances of interspecies PV infection were between Macaca mulatta and M. fascicularis, and Pan paniscus and P. troglodytes (Figure 1). As non-human primate interspecies infection has only been detected in Old World monkeys and apes thus far, this study represents the first characterization of interspecies PV infection in strepsirrhines.
The Varecia PVs form a distinct cluster, likely representing at least one new genus, within a well-supported clade consisting of numerous established genera (Figure 4). In the genera Dyoxipapillomavirus and Pipapillomavirus, respectively, BPV7 and Phodopus sungorus papillomavirus (PsuPV)-1 types are members of this clade known to be capable of interspecies infection and tumor induction. BPV7 has been isolated from cutaneous papillomas in cattle [33], while PsuPV1 infection has been found to occasionally result in oral squamous cell carcinoma in hamsters [9]. Based on health exams completed by the Duke Lemur Center veterinary staff, the VavPVs and VarPV1 infections appear to be asymptomatic. Prior to this study, there were no instances of interspecies infection between non-human primates in the aforementioned clade, as MfPV3, MfPV8, MfPV11, MmPV1, and PpPV1 (Figure 1) all belong to the Alphapapillomavirus genus. This study expands our understanding of the diversity of PVs that can undergo interspecies infection in primate sister species.
Interspecies PV infection in non-human primates has thus far only been detected in evolutionarily closely related species within the same genus. V. variegata and V. rubra are closely related lemur species comprising the Varecia genus within the Lemuridae family. V. rubra was previously considered a subspecies of V. variegata until evidence supported its classification into a separate species, beginning in 2001 [69,70]. In Madagascar, the range of V. rubra is primarily restricted to the Masoala Peninsula in northeastern Madagascar, whereas that of V. variegata stretches through rainforest parcels from northeastern to southeastern Madagascar [52,53,71]. Providing evidence for their evolutionary relatedness, the geographic ranges of V. variegata and V. rubra have historically overlapped, resulting in reports of potential hybridization throughout time, based primarily on intermediate fur patterns, with hybridization likely to have been a rare occurrence [70,71]. V. variegata and V. rubra hybrids have also occurred in the past for the captive population studied, as Varecia were allowed to hybridize in the earlier history of the Duke Lemur Center, when it was known as the Duke University Primate Center, resulting in 53 hybrids (including 50/50, 7/8th, and 15/16th hybrids), 12 of which did not survive (stillborn or died shortly after birth). Although the M1 and J1 V. rubra individuals focused on for this study have no known level of hybridization in their pedigree, it is possible that relatives of M1 and J1 had interacted with hybrid individuals in the past. V. rubra and V. variegata’s close genetic relationship and exposure to one another through a captive environment are likely the drivers of the interspecies PV infection seen in this study.
Future work will delve more deeply into PV diversity across lemuriform primates through the sampling of additional host species, body sites (e.g., skin, anogenital region), and populations (i.e., wild versus captive). The same lemur individuals found to harbor PVs may be targeted for additional sampling across body regions to understand the cell tropism of VavPV1, VavPV2, and VarPV1.

5. Conclusions

In addition to expanding the known diversity of PVs, this study represents the first case of interspecies PV infection in strepsirrhines (VavPV2) and characterizes the third complete PV genome isolated from lemurs (VarPV1). The diversity of PVs characterized from two highly evolutionarily related lemur species, V. variegata and V. rubra, is likely just a glimpse into the undiscovered PV diversity circulating in over 100 species of lemurs. As V. variegata and V. rubra are critically endangered species with extensive global captive populations, understanding viral diversity in the Varecia genus is vital for the continued success of maintaining the health of captive populations and to provide a foundation for future comparisons to viruses found in wild Varecia. Lastly, the work here shows the value of viral metagenomics in recovering the complete genomes of viruses relevant to animal health, particularly for animals in which remarkably limited viral research has been conducted.

Author Contributions

Conceptualization, E.N.P., A.D.Y. and A.V.; methodology, E.N.P., S.K., J.M.C., M.R., E.E., A.D.Y. and A.V.; software, E.N.P. and A.V.; validation, E.N.P., A.D.Y. and A.V.; formal analysis, E.N.P., A.D.Y. and A.V.; investigation, E.N.P., S.K., J.M.C., M.R., E.E., A.D.Y. and A.V.; resources, E.N.P., E.E., A.D.Y. and A.V.; data curation, E.N.P. and A.V.; writing—original draft preparation, E.N.P., A.D.Y. and A.V.; writing—review and editing, E.N.P., S.K., J.M.C., M.R., E.E., A.D.Y. and A.V.; visualization, E.N.P. and A.V.; supervision, E.E., A.D.Y. and A.V.; project administration, E.N.P. and A.D.Y.; funding acquisition, E.N.P. and A.D.Y. All authors have read and agreed to the published version of the manuscript.

Funding

The work described here was supported by TriCEM (Triangle Center for Evolutionary Medicine), Duke Lemur Center, Duke Biology, and Sigma Xi grants awarded to ENP.

Institutional Review Board Statement

Animal handling and sample collection protocols were approved by the Office of Animal Welfare Assurance at Duke University, USA (IACUC #A161-21-08).

Informed Consent Statement

Not applicable.

Data Availability Statement

The data for this study are available under BioProject PRJNA874427. The raw reads for this study are deposited under SRR26324872- SRR26324875. The papillomavirus genomes described in this study have been deposited in GenBank under accessions OR734654–OR734656. The mitochondrial genomes described in this study have been deposited in GenBank under accessions OR711366–OR711367.

Acknowledgments

We particularly thank the Duke Lemur Center research and veterinary teams who aided in lemur saliva sample collection and lemur health exams. This is Duke Lemur Center publication #1581.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Van Doorslaer, K.; Chen, Z.; Bernard, H.U.; Chan, P.K.; DeSalle, R.; Dillner, J.; Forslund, O.; Haga, T.; McBride, A.A.; Villa, L.L. ICTV Virus Taxonomy Profile: Papillomaviridae. J. Gen. Virol. 2018, 99, 989–990. [Google Scholar] [CrossRef] [PubMed]
  2. Rector, A.; Van Ranst, M. Animal papillomaviruses. Virology 2013, 445, 213–223. [Google Scholar] [CrossRef] [PubMed]
  3. Muñoz, N.; Bosch, F.X.; Jensen, O.M.; International Agency for Research on Cancer. Cancerregisteret (Denmark). In Human Papillomavirus and Cervical Cancer; International Agency for Research on Cancer; Oxford University Press: Lyon, NY, USA, 1989; Volume xii, 155p. [Google Scholar]
  4. Joh, J.; Hopper, K.; Van Doorslaer, K.; Sundberg, J.P.; Jenson, A.B.; Ghim, S.J. Macaca fascicularis papillomavirus type 1: A non-human primate betapapillomavirus causing rapidly progressive hand and foot papillomatosis. J. Gen. Virol. 2009, 90, 987–994. [Google Scholar] [CrossRef] [PubMed]
  5. Bergin, I.L.; Bell, J.D.; Chen, Z.; Zochowski, M.K.; Chai, D.; Schmidt, K.; Culmer, D.L.; Aronoff, D.M.; Patton, D.L.; Mwenda, J.M.; et al. Novel Genital Alphapapillomaviruses in Baboons (Papio hamadryas Anubis) With Cervical Dysplasia. Vet. Pathol. 2013, 50, 200–208. [Google Scholar] [CrossRef] [PubMed]
  6. Denny, L.A.; Franceschi, S.; de Sanjosé, S.; Heard, I.; Moscicki, A.B.; Palefsky, J. Human Papillomavirus, Human Immunodeficiency Virus and Immunosuppression. Vaccine 2012, 30, F168–F174. [Google Scholar] [CrossRef]
  7. Wang, J.; Guo, Y.; Wang, H.; Li, Y.; Zhang, L.; Wang, Z.; Song, L.; Liu, H. Genetic diversity of E6, E7 and the long control region in human papillomavirus type 16 variants in Beijing, China. Biochem. Biophys. Rep. 2022, 31, 101286. [Google Scholar] [CrossRef] [PubMed]
  8. Long, T.; Burk, R.D.; Chan, P.K.S.; Chen, Z. Non-human primate papillomavirus E6-mediated p53 degradation reveals ancient evolutionary adaptation of carcinogenic phenotype to host niche. PLoS Pathog. 2022, 18, e1010444. [Google Scholar] [CrossRef]
  9. Gimpelj Domjanič, G.; Hošnjak, L.; Lunar, M.M.; Skubic, L.; Zorec, T.M.; Račnik, J.; Cigler, B.; Poljak, M. First Report of Phodopus sungorus Papillomavirus Type 1 Infection in Roborovski Hamsters (Phodopus roborovskii). Viruses 2021, 13, 739. [Google Scholar] [CrossRef]
  10. Canuti, M.; Munro, H.J.; Robertson, G.J.; Kroyer, A.N.; Roul, S.; Ojkic, D.; Whitney, H.G.; Lang, A.S. New insight into avian papillomavirus ecology and evolution from characterization of novel wild bird papillomaviruses. Front. Microbiol. 2019, 10, 701. [Google Scholar] [CrossRef]
  11. Gottschling, M.; Stamatakis, A.; Nindl, I.; Stockfleth, E.; Alonso, Á.; Bravo, I.G. Multiple Evolutionary Mechanisms Drive Papillomavirus Diversification. Mol. Biol. Evol. 2007, 24, 1242–1258. [Google Scholar] [CrossRef]
  12. Trewby, H.; Ayele, G.; Borzacchiello, G.; Brandt, S.; Campo, M.S.; Del Fava, C.; Marais, J.; Leonardi, L.; Vanselow, B.; Biek, R.; et al. Analysis of the long control region of bovine papillomavirus type 1 associated with sarcoids in equine hosts indicates multiple cross-species transmission events and phylogeographical structure. J. Gen. Virol. 2014, 95, 2748–2756. [Google Scholar] [CrossRef] [PubMed]
  13. García-Pérez, R.; Ibáñez, C.; Godínez, J.M.; Aréchiga, N.; Garin, I.; Pérez-Suárez, G.; de Paz, O.; Juste, J.; Echevarría, J.E.; Bravo, I.G. Novel papillomaviruses in free-ranging Iberian bats: No virus-host co-evolution, no strict host specificity, and hints for recombination. Genome. Biol. Evol. 2014, 6, 94–104. [Google Scholar] [CrossRef] [PubMed]
  14. Varsani, A.; van der Walt, E.; Heath, L.; Rybicki, E.P.; Williamson, A.L.; Martin, D.P. Evidence of ancient papillomavirus recombination. J. Gen. Virol. 2006, 87, 2527–2531. [Google Scholar] [CrossRef] [PubMed]
  15. Otten, N.; von Tscharner, C.; Lazary, S.; Antczak, D.F.; Gerber, H. DNA of bovine papillomavirus type 1 and 2 in equine sarcoids: PCR detection and direct sequencing. Arch. Virol. 1993, 132, 121–131. [Google Scholar] [CrossRef] [PubMed]
  16. Lunardi, M.; Alcântara, B.K.d.; Otonel, R.A.A.; Rodrigues, W.B.; Alfieri, A.F.; Alfieri, A.A. Bovine Papillomavirus Type 13 DNA in Equine Sarcoids. J. Clin. Microbiol. 2013, 51, 2167–2171. [Google Scholar] [CrossRef] [PubMed]
  17. Martens, A.; De Moor, A.; Demeulemeester, J.; Peelman, L. Polymerase chain reaction analysis of the surgical margins of equine sarcoids for bovine papilloma virus DNA. Vet. Surg. 2001, 30, 460–467. [Google Scholar] [CrossRef] [PubMed]
  18. Kidney, B.A.; Berrocal, A. Sarcoids in two captive tapirs (Tapirus bairdii): Clinical, pathological and molecular study. Vet. Dermatol. 2008, 19, 380–384. [Google Scholar] [CrossRef] [PubMed]
  19. Roperto, S.; Russo, V.; Ozkul, A.; Corteggio, A.; Sepici-Dincel, A.; Catoi, C.; Esposito, I.; Riccardi, M.G.; Urraro, C.; Luca, R. Productive infection of bovine papillomavirus type 2 in the urothelial cells of naturally occurring urinary bladder tumors in cattle and water buffaloes. PLoS ONE 2013, 8, e62227. [Google Scholar] [CrossRef]
  20. Silvestre, O.; Borzacchiello, G.; Nava, D.; Iovane, G.; Russo, V.; Vecchio, D.; D’ausilio, F.; Gault, E.; Campo, M.; Paciello, O. Bovine papillomavirus type 1 DNA and E5 oncoprotein expression in water buffalo fibropapillomas. Vet. Pathol. 2009, 46, 636–641. [Google Scholar] [CrossRef]
  21. Munday, J.S.; Thomson, N.; Dunowska, M.; Knight, C.G.; Laurie, R.E.; Hills, S. Genomic characterisation of the feline sarcoid-associated papillomavirus and proposed classification as Bos taurus papillomavirus type 14. Vet. Microbiol. 2015, 177, 289–295. [Google Scholar] [CrossRef]
  22. Orbell, G.; Young, S.; Munday, J. Cutaneous sarcoids in captive African lions associated with feline sarcoid-associated papillomavirus infection. Vet. Pathol. 2011, 48, 1176–1179. [Google Scholar] [CrossRef] [PubMed]
  23. Pfister, H.; Fink, B.; Thomas, C. Extrachromosomal bovine papillomavirus type 1 DNA in hamster fibromas and fibrosarcomas. Virology 1981, 115, 414–418. [Google Scholar] [CrossRef] [PubMed]
  24. Robl, M.G.; Olson, C. Oncogenic Action of Bovine Papilloma Virus in Hamsters1. Cancer Res. 1968, 28, 1596–1604. [Google Scholar] [PubMed]
  25. Boiron, M.; Levy, J.P.; Thomas, M.; Friedmann, J.C.; Bernard, J. Some Properties of Bovine Papilloma Virus. Nature 1964, 201, 423–424. [Google Scholar] [CrossRef] [PubMed]
  26. Sundberg, J.; Reszka, A.; Williams, E.; Reichmann, M. An oral papillomavirus that infected one coyote and three dogs. Vet. Pathol. 1991, 28, 87–88. [Google Scholar] [CrossRef] [PubMed]
  27. De Falco, F.; Cuccaro, B.; De Tullio, R.; Alberti, A.; Cutarelli, A.; De Carlo, E.; Roperto, S. Possible etiological association of ovine papillomaviruses with bladder tumors in cattle. Virus. Res. 2023, 328, 199084. [Google Scholar] [CrossRef] [PubMed]
  28. Savini, F.; Gallina, L.; Prosperi, A.; Puleio, R.; Lavazza, A.; Di Marco, P.; Tumino, S.; Moreno, A.; Lelli, D.; Guercio, A.; et al. Bovine Papillomavirus 1 Gets Out of the Flock: Detection in an Ovine Wart in Sicily. Pathogens 2020, 9, 429. [Google Scholar] [CrossRef]
  29. Russo, V.; Roperto, F.; De Biase, D.; Cerino, P.; Urraro, C.; Munday, J.S.; Roperto, S. Bovine Papillomavirus Type 2 Infection Associated with Papillomatosis of the Amniotic Membrane in Water Buffaloes (Bubalus bubalis). Pathogens 2020, 9, 262. [Google Scholar] [CrossRef]
  30. Reid, S.; Smith, K.; Jarrett, W. Detection, cloning and characterisation of papillomaviral DNA present in sarcoid tumours of Equus asinus. Vet. Rec. 1994, 135, 430–432. [Google Scholar] [CrossRef]
  31. Savini, F.; Dal Molin, E.; Gallina, L.; Casà, G.; Scagliarini, A. Papillomavirus in healthy skin and mucosa of wild ruminants in the Italian Alps. J. Wildl. Dis. 2016, 52, 82–87. [Google Scholar] [CrossRef]
  32. Chen, L.; Liu, B.; Yang, J.; Jin, Q. DBatVir: The database of bat-associated viruses. Database 2014, 2014, bau021. [Google Scholar] [CrossRef] [PubMed]
  33. Ogawa, T.; Tomita, Y.; Okada, M.; Shirasawa, H. Complete genome and phylogenetic position of bovine papillomavirus type 7. J. Gen. Virol. 2007, 88, 1934–1938. [Google Scholar] [CrossRef] [PubMed]
  34. Lawson, B.; Robinson, R.A.; Fernandez, J.R.; John, S.K.; Benitez, L.; Tolf, C.; Risely, K.; Toms, M.P.; Cunningham, A.A.; Williams, R.A.J. Spatio-temporal dynamics and aetiology of proliferative leg skin lesions in wild British finches. Sci. Rep. 2018, 8, 14670. [Google Scholar] [CrossRef] [PubMed]
  35. Dunay, E.; Rukundo, J.; Atencia, R.; Cole, M.F.; Cantwell, A.; Emery Thompson, M.; Rosati, A.G.; Goldberg, T.L. Viruses in saliva from sanctuary chimpanzees (Pan troglodytes) in Republic of Congo and Uganda. PLoS ONE 2023, 18, e0288007. [Google Scholar] [CrossRef] [PubMed]
  36. Van Ranst, M.; Fuse, A.; Fiten, P.; Beuken, E.; Pfister, H.; Burk, R.D.; Opdenakker, G. Human papillomavirus type 13 and pygmy chimpanzee papillomavirus type 1: Comparison of the genome organizations. Virology 1992, 190, 587–596. [Google Scholar] [CrossRef] [PubMed]
  37. Root, J.J.; Hopken, M.W.; Gidlewski, T.; Piaggio, A.J. Cottontail Rabbit Papillomavirus Infection in a Desert cottontail (Sylvilagus audubonii) from Colorado, USA. J. Wildl. Dis. 2013, 49, 1060–1062. [Google Scholar] [CrossRef] [PubMed]
  38. Chen, Z.; van Doorslaer, K.; DeSalle, R.; Wood, C.E.; Kaplan, J.R.; Wagner, J.D.; Burk, R.D. Genomic diversity and interspecies host infection of α12 Macaca fascicularis papillomaviruses (MfPVs). Virology 2009, 393, 304–310. [Google Scholar] [CrossRef] [PubMed]
  39. Literák, I.; Tomita, Y.; Ogawa, T.; Shirasawa, H.; Šmid, B.; Novotný, L.; Adamec, M. Papillomatosis in a European bison. J. Wildl. Dis. 2006, 42, 149–153. [Google Scholar] [CrossRef]
  40. Munday, J.S.; Fairley, R.; Lowery, I. Detection of Ovis aries papillomavirus type 2 DNA sequences in a sarcoid-like mass in the mouth of a pig. Vet. Microbiol. 2020, 248, 108801. [Google Scholar] [CrossRef]
  41. Gysens, L.; Vanmechelen, B.; Maes, P.; Martens, A.; Haspeslagh, M. Complete genomic characterization of bovine papillomavirus type 1 and 2 strains infers ongoing cross-species transmission between cattle and horses. Vet. J. 2023, 298–299, 106011. [Google Scholar] [CrossRef]
  42. Bogaert, L.; Martens, A.; Van Poucke, M.; Ducatelle, R.; De Cock, H.; Dewulf, J.; De Baere, C.; Peelman, L.; Gasthuys, F. High prevalence of bovine papillomaviral DNA in the normal skin of equine sarcoid-affected and healthy horses. Vet. Microbiol. 2008, 129, 58–68. [Google Scholar] [CrossRef]
  43. Bengis, R.; Van Heerden, J.; Venter, E.; Bosman, A.; Van Dyk, E.; Williams, J.; Van Wilpe, E. Detection and characterisation of papillomavirus in skin lesions of giraffe and sable antelope in South Africa. J. South Afr. Vet. Assoc. 2011, 82, 80–85. [Google Scholar]
  44. Löhr, C.V.; Juan-Sallés, C.; Rosas-Rosas, A.; García, A.P.; Garner, M.M.; Teifke, J.P. Sarcoids in captive zebras (Equus burchellii): Association with bovine papillomavirus type 1 infection. J. Zoo Wildl. Med. 2005, 36, 74–81. [Google Scholar] [CrossRef] [PubMed]
  45. Gallina, L.; Savini, F.; Casà, G.; Bertoletti, I.; Bianchi, A.; Gibelli, L.R.; Lelli, D.; Lavazza, A.; Scagliarini, A. Epitheliotropic infections in wildlife ruminants from the Central Alps and Stelvio National Park. Front. Vet. Sci. 2020, 7, 229. [Google Scholar] [CrossRef] [PubMed]
  46. De Falco, F.; Cutarelli, A.; Cuccaro, B.; Catoi, C.; De Carlo, E.; Roperto, S. Evidence of a novel cross-species transmission by ovine papillomaviruses. Transbound Emerg. Dis. 2022, 69, 3850–3857. [Google Scholar] [CrossRef]
  47. Davies, N.; Johnson, S.; Louis, E.E.; Mittermeier, R.A.; Nash, S.D.; Rajaobelina, S.; Ratsimbazafy, J.; Razafindramanana, J.; Schwitzer, C. Lemurs of Madagascar: A Strategy for Their Conservation 2013–2016; (IUCN) International Union for Conservation of Nature; Bristol Zoo Gardens; Conservation International; IUCN Species Survival Commission (SSC); Primate Specialist Group: Bristol, UK, 2013. [Google Scholar]
  48. Mittermeier, R.A.; Wilson, D.E. Handbook of the Mammals of the World—Volume 3: Primates; LYNX Nature Book: Spain, Barcelona, 2013. [Google Scholar]
  49. Rylands, A.B.; Mittermeier, R.A. Primate taxonomy: Species and conservation. Evol. Anthropol. Issues News Rev. 2014, 23, 8–10. [Google Scholar] [CrossRef]
  50. Perez, S.I.; Tejedor, M.F.; Novo, N.M.; Aristide, L. Divergence Times and the Evolutionary Radiation of New World Monkeys (Platyrrhini, Primates): An Analysis of Fossil and Molecular Data. PLoS ONE 2013, 8, e68029. [Google Scholar] [CrossRef]
  51. Paietta, E.N.; Kraberger, S.; Custer, J.M.; Vargas, K.L.; Van Doorslaer, K.; Yoder, A.D.; Varsani, A. Identification of diverse papillomaviruses in captive black-and-white ruffed lemurs (Varecia variegata). Arch. Virol. 2022, 168, 13. [Google Scholar] [CrossRef]
  52. Borgerson, C.; Eppley, T.M.; Patel, E.; Johnson, S.; Louis, E.E.; Razafindramanana, J. Varecia rubra. The IUCN Red List of Threatened Species. 2020. Available online: https://www.iucnredlist.org/species/22920/115574598 (accessed on 18 November 2023).
  53. Louis, E.E.; Sefczek, T.M.; Raharivololona, B.; King, T.; Morelli, T.L.; Baden, A. Varecia variegata. The IUCN Red List of Threatened Species. 2020. Available online: https://www.iucnredlist.org/species/22918/115574178 (accessed on 18 November 2023).
  54. Bolger, A.M.; Lohse, M.; Usadel, B. Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinformatics 2014, 30, 2114–2120. [Google Scholar] [CrossRef]
  55. Li, D.; Liu, C.-M.; Luo, R.; Sadakane, K.; Lam, T.-W. MEGAHIT: An ultra-fast single-node solution for large and complex metagenomics assembly via succinct de Bruijn graph. Bioinformatics 2015, 31, 1674–1676. [Google Scholar] [CrossRef]
  56. Buchfink, B.; Xie, C.; Huson, D.H. Fast and sensitive protein alignment using DIAMOND. Nat. Methods 2015, 12, 59–60. [Google Scholar] [CrossRef] [PubMed]
  57. Altschul, S.F.; Gish, W.; Miller, W.; Myers, E.W.; Lipman, D.J. Basic local alignment search tool. J. Mol. Biol. 1990, 215, 403–410. [Google Scholar] [CrossRef] [PubMed]
  58. Tisza, M.J.; Belford, A.K.; Domínguez-Huerta, G.; Bolduc, B.; Buck, C.B. Cenote-Taker 2 democratizes virus discovery and sequence annotation. Virus Evol. 2021, 7, veaa100. [Google Scholar] [CrossRef] [PubMed]
  59. Van Doorslaer, K.; Li, Z.; Xirasagar, S.; Maes, P.; Kaminsky, D.; Liou, D.; Sun, Q.; Kaur, R.; Huyen, Y.; McBride, A.A. The Papillomavirus Episteme: A major update to the papillomavirus sequence database. Nucleic Acids Res. 2017, 45, D499–D506. [Google Scholar] [CrossRef] [PubMed]
  60. Katoh, K.; Standley, D.M. MAFFT multiple sequence alignment software version 7: Improvements in performance and usability. Mol. Biol. Evol. 2013, 30, 772–780. [Google Scholar] [CrossRef] [PubMed]
  61. Capella-Gutiérrez, S.; Silla-Martínez, J.M.; Gabaldón, T. trimAl: A tool for automated alignment trimming in large-scale phylogenetic analyses. Bioinformatics 2009, 25, 1972–1973. [Google Scholar] [CrossRef] [PubMed]
  62. Darriba, D.; Taboada, G.L.; Doallo, R.; Posada, D. ProtTest 3: Fast selection of best-fit models of protein evolution. Bioinformatics 2011, 27, 1164–1165. [Google Scholar] [CrossRef]
  63. Minh, B.Q.; Schmidt, H.A.; Chernomor, O.; Schrempf, D.; Woodhams, M.D.; von Haeseler, A.; Lanfear, R. IQ-TREE 2: New Models and Efficient Methods for Phylogenetic Inference in the Genomic Era. Mol. Biol. Evol. 2020, 37, 1530–1534. [Google Scholar] [CrossRef]
  64. Letunic, I.; Bork, P. Interactive Tree Of Life (iTOL) v5: An online tool for phylogenetic tree display and annotation. Nucleic Acids Res. 2021, 49, W293–W296. [Google Scholar] [CrossRef]
  65. Bernt, M.; Donath, A.; Jühling, F.; Externbrink, F.; Florentz, C.; Fritzsch, G.; Pütz, J.; Middendorf, M.; Stadler, P.F. MITOS: Improved de novo metazoan mitochondrial genome annotation. Mol. Phylogenetics Evol. 2013, 69, 313–319. [Google Scholar] [CrossRef]
  66. Roman, A.; Munger, K. The papillomavirus E7 proteins. Virology 2013, 445, 138–168. [Google Scholar] [CrossRef] [PubMed]
  67. Vega, R.; Hopper, J.; Kitchener, A.C.; Catinaud, J.; Roullet, D.; Robsomanitrandrasana, E.; Hollister, J.D.; Roos, C.; King, T. The mitochondrial DNA diversity of captive ruffed lemurs (Varecia spp.): Implications for conservation. Oryx 2023, 57, 649–658. [Google Scholar] [CrossRef]
  68. Schwitzer, C.; King, T.; Robsomanitrandrasana, E.; Chamberlan, C.; Rasolofoharivelo, T. Integrating ex situ and in situ conservation of lemurs. In Lemurs of Madagascar: A Strategy for Their Conservation; Research Gate: Berlin, Germany, 2013; Volume 2016, pp. 146–152. [Google Scholar]
  69. Groves, C. Primate Taxonomy; Smithsonian Institution Press: Washington, DC, USA, 2001. [Google Scholar]
  70. Vasey, N.; Tattersall, I. Do ruffed lemurs form a hybrid zone? Distribution and discovery of Varecia, with systematic and conservation implications. Am. Mus. Novit. 2002, 2002, 1–26. [Google Scholar] [CrossRef]
  71. Hekkala, E.R.; Rakotondratsima, M.; Vasey, N. Habitat and distribution of the ruffed lemur, Varecia, north of the Bay of Antongil in northeastern Madagascar. Primate Conserv. 2007, 22, 89–95. [Google Scholar] [CrossRef]
Viruses 16 00037 g001
Figure 1. Summary of evidence available throughout the literature and NCBI Virus (https://www.ncbi.nlm.nih.gov/labs/virus/vssi/ (accessed on 1 September 2023)) of naturally occurring interspecies infection in mammalian and avian orders. PVs are connected via solid vertical lines to the animal species they have been found to infect. Abbreviations and references for PV types are as follows: Anas platyrhynchos papillomavirus 1 (AplaPV1) [10], bovine papillomavirus types −1, −2, −7, −8, −13, −14 (BPV1, BPV2, BPV7, BPV8, BPV13, BPV14) [12,16,18,19,20,21,22,28,29,30,31,32,33,39,41,42,43,44,45], Eptesicus serotinus papillomavirus types −2, −3 (EsPV2, EsPV3) [13], Fringilla coelebs papillomavirus 1 (FcPV1) [34], Larus smithsonianus papillomavirus 1 (LsmiPV1) [10], Macaca fascicularis papillomavirus types −3, −8, −11 (MfPV3, MfPV8, MfPV11) [38], Macaca mulatta papillomavirus 1 (MmPV1) [38], Ovis aries papillomavirus type −1, −2, −4 (OaPV1, OaPV2, OaPV4) [27,40,46], Pan paniscus papillomavirus 1 (PpPV1) [35,36], Phodopus sungorus papillomavirus type 1 (PsuPV1) [9], Sylvilagus floridanus papillomavirus 1 (SfPV1) [37], Varecia variegata papillomavirus 2 (VavPV2).
Figure 1. Summary of evidence available throughout the literature and NCBI Virus (https://www.ncbi.nlm.nih.gov/labs/virus/vssi/ (accessed on 1 September 2023)) of naturally occurring interspecies infection in mammalian and avian orders. PVs are connected via solid vertical lines to the animal species they have been found to infect. Abbreviations and references for PV types are as follows: Anas platyrhynchos papillomavirus 1 (AplaPV1) [10], bovine papillomavirus types −1, −2, −7, −8, −13, −14 (BPV1, BPV2, BPV7, BPV8, BPV13, BPV14) [12,16,18,19,20,21,22,28,29,30,31,32,33,39,41,42,43,44,45], Eptesicus serotinus papillomavirus types −2, −3 (EsPV2, EsPV3) [13], Fringilla coelebs papillomavirus 1 (FcPV1) [34], Larus smithsonianus papillomavirus 1 (LsmiPV1) [10], Macaca fascicularis papillomavirus types −3, −8, −11 (MfPV3, MfPV8, MfPV11) [38], Macaca mulatta papillomavirus 1 (MmPV1) [38], Ovis aries papillomavirus type −1, −2, −4 (OaPV1, OaPV2, OaPV4) [27,40,46], Pan paniscus papillomavirus 1 (PpPV1) [35,36], Phodopus sungorus papillomavirus type 1 (PsuPV1) [9], Sylvilagus floridanus papillomavirus 1 (SfPV1) [37], Varecia variegata papillomavirus 2 (VavPV2).
Viruses 16 00037 g001
Viruses 16 00037 g002
Figure 2. Annotations of complete PV genomes characterized from V. rubra saliva. Two of the complete genomes can be classified as VavPV2, a PV previously determined from V. variegata saliva. One of the complete genomes belongs to a new type and species and has been named Varecia rubra papillomavirus 1 (VarPV1).
Figure 2. Annotations of complete PV genomes characterized from V. rubra saliva. Two of the complete genomes can be classified as VavPV2, a PV previously determined from V. variegata saliva. One of the complete genomes belongs to a new type and species and has been named Varecia rubra papillomavirus 1 (VarPV1).
Viruses 16 00037 g002
Viruses 16 00037 g003
Figure 3. Summary of the conserved zinc-binding motifs (CxxC) in the E6 and E7 proteins, and conserved regions 1 and 2 and pRB-binding motif (Lx[C/S]xE) in the E7 protein of all lemur papillomaviruses.
Figure 3. Summary of the conserved zinc-binding motifs (CxxC) in the E6 and E7 proteins, and conserved regions 1 and 2 and pRB-binding motif (Lx[C/S]xE) in the E7 protein of all lemur papillomaviruses.
Viruses 16 00037 g003
Viruses 16 00037 g004
Figure 4. Partitioned maximum likelihood phylogenetic tree of concatenated amino acid sequences of E1, E2, and L1 showing the relationship of the lemur papillomaviruses with their nearest neighbors.
Figure 4. Partitioned maximum likelihood phylogenetic tree of concatenated amino acid sequences of E1, E2, and L1 showing the relationship of the lemur papillomaviruses with their nearest neighbors.
Viruses 16 00037 g004
Viruses 16 00037 g005
Figure 5. Maximum likelihood phylogenetic tree of host mitochondrial genomes. Only a subset of species within the Lemuridae family (and outgroup Indriidae) for which mitochondrial genomes were available are depicted. V. rubra mitochondrial genomes characterized in this study and publicly available V. variegata mitochondrial genomes are depicted in blue to allow for comparison. Branch support was determined with ultrafast bootstrap (UFBoot) approximation using IQ-TREE 2.
Figure 5. Maximum likelihood phylogenetic tree of host mitochondrial genomes. Only a subset of species within the Lemuridae family (and outgroup Indriidae) for which mitochondrial genomes were available are depicted. V. rubra mitochondrial genomes characterized in this study and publicly available V. variegata mitochondrial genomes are depicted in blue to allow for comparison. Branch support was determined with ultrafast bootstrap (UFBoot) approximation using IQ-TREE 2.
Viruses 16 00037 g005
Table 1. For each non-human primate (NHP) superfamily, the approximate number of extant NHP species is compared to the NHP species with PV sequences, including both partial sequences and complete genomes, available in NCBI. Additionally, the number of complete PV genomes available in NCBI for each primate superfamily is displayed. The data, obtained from NCBI Virus (https://www.ncbi.nlm.nih.gov/labs/virus/vssi/ (accessed on 18 November 2023)), emphasizes the scarcity of PV data across speciose NHP lineages. The six PV genomes available for lemurs include the three genomes characterized previously and the three genomes described in this study.
Table 1. For each non-human primate (NHP) superfamily, the approximate number of extant NHP species is compared to the NHP species with PV sequences, including both partial sequences and complete genomes, available in NCBI. Additionally, the number of complete PV genomes available in NCBI for each primate superfamily is displayed. The data, obtained from NCBI Virus (https://www.ncbi.nlm.nih.gov/labs/virus/vssi/ (accessed on 18 November 2023)), emphasizes the scarcity of PV data across speciose NHP lineages. The six PV genomes available for lemurs include the three genomes characterized previously and the three genomes described in this study.
NHP SuperfamilyApprox. Number of Extant SpeciesSpecies with PV
Sequences (Partial and Complete) in NCBI		Number of
Complete PV
Genomes
Available in NCBI
Ceboidea
(New World Monkeys)		>100Alouatta caraya
Alouatta guariba
Callithrix penicillata
Saimiri sciureus
Sapajus sp.		8
Cercopithecoidea
(Old World Monkeys)		>130Colobus guereza
Macaca fascicularis
Macaca fuscata
Macaca mulatta
Papio hamadryas
Piliocolobus tephrosceles		23
Hominoidea
(Apes, excludes humans)		~25Gorilla gorilla
Pan paniscus
Pan troglodytes		4
Lemuroidea
(Lemurs)		>100Varecia variegata
Varecia rubra		6
Lorisoidea
(Lorisids & Galagos)		>25-0
Tarsioidea
(Tarsiers)		>10-0
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Paietta, E.N.; Kraberger, S.; Regney, M.; Custer, J.M.; Ehmke, E.; Yoder, A.D.; Varsani, A. Interspecies Papillomavirus Type Infection and a Novel Papillomavirus Type in Red Ruffed Lemurs (Varecia rubra). Viruses 2024, 16, 37. https://doi.org/10.3390/v16010037

AMA Style

Paietta EN, Kraberger S, Regney M, Custer JM, Ehmke E, Yoder AD, Varsani A. Interspecies Papillomavirus Type Infection and a Novel Papillomavirus Type in Red Ruffed Lemurs (Varecia rubra). Viruses. 2024; 16(1):37. https://doi.org/10.3390/v16010037

Chicago/Turabian Style

Paietta, Elise N., Simona Kraberger, Melanie Regney, Joy M. Custer, Erin Ehmke, Anne D. Yoder, and Arvind Varsani. 2024. "Interspecies Papillomavirus Type Infection and a Novel Papillomavirus Type in Red Ruffed Lemurs (Varecia rubra)" Viruses 16, no. 1: 37. https://doi.org/10.3390/v16010037

APA Style

Paietta, E. N., Kraberger, S., Regney, M., Custer, J. M., Ehmke, E., Yoder, A. D., & Varsani, A. (2024). Interspecies Papillomavirus Type Infection and a Novel Papillomavirus Type in Red Ruffed Lemurs (Varecia rubra). Viruses, 16(1), 37. https://doi.org/10.3390/v16010037

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop