Complete Genome Sequences of Six Novel Macaca mulatta Papillomavirus Types Isolated from Genital Sites of Rhesus Monkeys in Hong Kong SAR, China

Identification and evolutionary analysis of papillomavirus sequences in New World monkeys (genera Sapajus and Alouatta) from Argentina

OBJECTIVE

In this study, we investigated the occurrence of papillomavirus (PV) infection in non-human primates (NHPs) in northeastern Argentina. We also explored their evolutionary history and evaluated the co-speciation hypothesis in the context of primate evolution.

METHODS

We obtained DNA samples from 57 individuals belonging to wild and captive populations of Alouatta caraya, Sapajus nigritus, and Sapajus cay. We assessed PV infection by PCR amplification with the CUT primer system and sequencing of 337 bp (112 amino acids) of the L1 gene. The viral sequences were analyzed by phylogenetic and Bayesian coalescence methods to estimate the time to the most common recent ancestor (t_MRCA) using BEAST, v1.4.8 software. We evaluated viral/host tree congruence with TreeMap v3.0.

RESULTS

We identified two novel putative PV sequences of the genus Gammapapillomavirus in Sapajus spp. and Alouatta caraya (SPV1 and AcPV1, respectively). The t_MRCA of SPV1 was estimated to be 11,941,682 years before present (ybp), and that of AcPV1 was 46,638,071 ybp, both before the coalescence times of their hosts (6.4 million years ago [MYA] and 6.8 MYA, respectively). Based on the comparison of primate and viral phylogenies, we found that the PV tree was no more congruent with the host tree than a random tree would be (P > 0.05), thus allowing us to reject the model of virus-host coevolution.

CONCLUSION

This study presents the first evidence of PV infection in platyrrhine species from Argentina, expands the range of described hosts for these viruses, and suggests new scenarios for their origin and dispersal.

Molecular and Phylogenetic Characterization of Novel Papillomaviruses Isolated from Oral and Anogenital Neoplasms of Japanese Macaques (Macaca fuscata)

Papillomaviruses (PVs) are a diverse group of host species-specific DNA viruses, etiologically linked with various benign and malignant neoplasms of cutaneous and mucosal epithelia. Here, we describe the detection and characterization of the first two PVs naturally infecting Japanese macaques (Macaca fuscata), including the determination of their etiological association(s) with the development of original neoplasms. The molecular and phylogenetic analyses were performed on complete genome sequences of Macaca fuscata PV types 1 (MfuPV1) and 2 (MfuPV2), which were completely sequenced in samples of a malignant oral tumor and benign anogenital neoplasm of Japanese macaques, respectively. Subsequently, two type-specific quantitative real-time PCRs were developed to estimate viral loads of MfuPV1 and MfuPV2 and to evaluate their etiological roles. The in silico molecular analyses revealed that both viral genomes encode characteristic PV proteins with conserved functional domains and have a non-coding genomic region with regulatory sequences to regulate and complete the viral life cycle. However, additional experimental evidence is needed to finally confirm the presence and biological functionality of the molecular features of both novel PVs. While MfuPV1, together with PVs identified in other macaques, is classified into the Alphapapillomavirus (Alpha-PV) species 12, MfuPV2 is most likely a representative of the novel viral species within the Alpha-PV genus. Their relatively high viral loads suggest that both PVs are etiologically linked with the development of the original neoplasms.

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

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

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

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