Bovine Papillomavirus: Veterinary Reference

Overview and Taxonomy of Bovine Papillomavirus: Deltapapillomavirus Types 1, 2, and Beyond

The bovine papillomaviruses (BPVs) represent a remarkably heterogeneous assemblage of oncogenic, species-specific viruses within the family Papillomaviridae, a lineage of small, non-enveloped, double-stranded DNA viruses with a pronounced tropism for cutaneous and mucosal epithelia. The global burden of BPV-induced disease is substantial, manifesting as benign proliferative lesions, cutaneous papillomas (warts) and fibropapillomas, that, while often self-limiting, can cause significant economic losses due to reduced hide quality, diminished weight gain, and, in severe cases, astasia necessitating euthanasia [2, 13, 14]. More critically, specific BPV types, particularly those within the Deltapapillomavirus genus, are implicated in malignant transformation of the urinary bladder and upper gastrointestinal tract, as well as being the aetiological agents of equine sarcoids, the most common mesenchymal neoplasm of horses worldwide [1, 3, 12, 18]. The World Organisation for Animal Health (WOAH) recognizes the economic and zoonotic potential of papillomaviruses, underscoring the importance of accurate taxonomic classification for epidemiological surveillance and control.

Taxonomic Classification within Papillomaviridae

The taxonomic framework for BPVs is rooted in the phylogenetic analysis of the highly conserved L1 major capsid protein gene, a gold standard endorsed by the International Committee on Taxonomy of Viruses (ICTV). This approach has delineated distinct genera, which generally correspond to biological properties, host range, and tissue tropism. The Deltapapillomavirus genus is of paramount clinical and veterinary significance, currently encompassing the archetypal BPV-1 and BPV-2, alongside BPV-13 and several emerging types [6, 18]. The designation of a new papillomavirus type requires that its complete L1 gene sequence share less than 90% identity with any known type, a criterion that has driven the discovery of a vast and growing diversity beyond the initially characterized types [8, 15]. While early classification relied heavily on antigenic properties and host species, modern molecular phylogenetics has revealed that BPVs are not a monophyletic group. Instead, they are distributed across at least five genera, including Deltapapillomavirus, Xipapillomavirus (e.g., BPV-3, -4, -6, -9, -10, -11, -15), Epsilonpapillomavirus (BPV-5, -8), and Dyoxipapillomavirus [8, 15]. The Deltapapillomavirus genus is unique for its capacity to induce both epithelial and mesenchymal proliferation, leading to the characteristic fibropapilloma, a tumor composed of proliferating fibroblasts overlaid by hyperplastic epithelium. This dual tropism underpins their role in equine sarcoids, where BPV-1 and BPV-2 DNA is consistently detected and considered causal [3, 10, 12].

Genetic Architecture and Phylogenetic Basis of Typing

The BPV genome is a circular, double-stranded DNA molecule of approximately 7.2 to 8.0 kilobase pairs, organized into an early region (E) encoding regulatory and oncogenic proteins, a late region (L) encoding the capsid proteins L1 and L2, and a non-coding upstream regulatory region (URR) controlling transcription and replication [2, 8]. For Deltapapillomavirus types 1 and 2, the early region comprises open reading frames (ORFs) for E5, E6, E7, E1, E2, and E4. The E5 protein, in particular, is a potent oncogene for BPV-1 and BPV-2, functioning through the activation of the platelet-derived growth factor β receptor (PDGFβR) and other signalling cascades that drive uncontrolled fibroblast proliferation. This is a defining molecular hallmark of the Deltapapillomavirus genus [6, 18]. Fundamental to the taxonomy is the L1 gene, which forms the viral capsid and is the primary target for phylogenetic classification. The L1 ORF of Deltapapillomavirus 4 (BPV-2) consistently exhibits high sequence conservation across global isolates, with studies from Egypt, Brazil, and China demonstrating over 99% nucleotide identity, suggesting a stable, globally circulating lineage [1]. Similarly, full-genome analyses of BPV-1 from severe papillomatosis cases in Japan and from Chinese indigenous cattle reveal a clear phylogenetic clustering into subtypes, with Asian isolates often forming more ancestral clades compared to European strains, indicative of long-term co-evolution and geographic segregation [2, 11]. The L1 gene is also the target for serological classification and vaccine development, as due to its ability to form virus-like particles (VLPs) that elicit a robust neutralizing antibody response [4, 5, 9].

Biological Properties and Pathogenesis of Deltapapillomaviruses

BPV-1 and BPV-2 are the principal agents of classical cutaneous papillomatosis in cattle, typically causing fibropapillomas on the head, neck, and trunk of young animals. However, their pathogenic reach extends far beyond the skin. A well-documented and economically devastating manifestation is their association with equine sarcoids, a condition that affects horses, donkeys, and mules globally. Crucially, BPV-1 and BPV-2 are not competent to complete their full replication cycle in equine tissues; the infection is abortive, with the viral genome persisting as extrac chromosomal episomes in fibroblasts, driving neoplastic transformation without producing infectious virions [3, 10, 12, 17]. This cross-species transmission, while rare among papillomaviruses, highlights the unique biological plasticity of the Deltapapillomavirus genus. In cattle, a second, profoundly important pathological pathway involves BPV-2 and the closely related BPV-13. These types are known to cause latent, abortive infections in the urothelium of the urinary bladder. When co-factors such as chronic ingestion of bracken fern (Pteridium aquilinum), which contains immunosuppressive and mutagenic compounds, are present, these latent infections can progress to invasive urothelial carcinoma [6, 18]. This model of viral co-carcinogenesis is a cornerstone in comparative oncology, illustrating the interplay between a viral oncogene (E5) and environmental mutagens. Research has elucidated the molecular underpinnings of this process, demonstrating that BPV-2 and BPV-13 E5 oncoprotein expression in urothelial cells triggers parkin-mediated mitophagy, involving the inner mitochondrial membrane protein Prohibitin 2 (PHB2), which is significantly overexpressed in neoplastic tissues [6]. Furthermore, these tumors express the innate immune receptor Mincle, a finding that opens avenues for understanding immune evasion and potential immunotherapeutic strategies, such as BCG therapy, which is already used in human bladder cancer [12, 18].

Beyond BPV-1 and BPV-2: Emerging Diversity in Deltapapillomavirus

The taxonomy of the Deltapapillomavirus genus continues to expand, driven by molecular surveillance of cattle populations. BPV-13 is now firmly established as a Deltapapillomavirus species, often co-infecting with BPV-2 and sharing a similar oncogenic profile, including the ability to drive E5-mediated urothelial carcinogenesis [6, 18]. Further complexity is revealed by multiple infection studies. Investigations using the FAP primer pair, which targets a conserved region of the L1 gene, have identified not only BPV-1 and BPV-2 but also putative novel Deltapapillomavirus types, such as BPV/BR-UEL3, BPV/BR-UEL4, and BPV/BR-UEL5, in Brazilian cattle herds [15]. These findings underscore that mixed infections with two or three distinct BPV types within a single animal are common, a phenomenon mirroring the diversity seen in human papillomavirus (HPV) infections. The presence of such high viral diversity within a single herd and even within individual lesions has profound implications for vaccine design and epidemiological tracking, as a vaccine targeting only BPV-1 and BPV-2 may not confer protection against emerging types [15]. The delineation of novel types, such as BPV-15 identified in Xinjiang dairy cows, which clusters within the Xipapillomavirus genus, further refines our understanding of the evolutionary relationships among BPVs [8]. As comprehensive genomic sequencing becomes more accessible, the catalogue of Deltapapillomavirus types will undoubtedly grow, revealing a complex viral ecosystem within global bovine populations. From a One Health perspective, understanding this taxonomy is not merely an academic exercise; it is foundational for the development of targeted diagnostic tools, effective prophylactic vaccines using L1 VLPs, and the rational design of therapeutic interventions for BPV-associated diseases in both cattle and horses [4, 5, 7, 9]. The World Health Organization (WHO) and WOAH both emphasize that robust surveillance and molecular characterization are critical for managing the impact of oncogenic viruses in animal populations and mitigating their potential zoonotic and economic consequences [16].

Molecular Pathogenesis and Host-Virus Interactions in Bovine Papillomavirus Infection

The molecular pathogenesis of bovine papillomavirus (BPV) infection represents a paradigm of viral oncogenesis, wherein a small, double-stranded DNA virus subverts host cellular machinery through a sophisticated interplay of early gene products, immune evasion strategies, and dysregulated signaling cascades. BPV, a member of the Papillomaviridae family, encompasses multiple genera, including Deltapapillomavirus, Xipapillomavirus, and Epsilonpapillomavirus, each exhibiting distinct tropism and pathogenic potential. The Deltapapillomavirus genus, particularly BPV types 1, 2, and 13, is of paramount clinical significance due to its causal association with both cutaneous fibropapillomas in cattle and the cross-species transmission to equids, where it induces sarcoids, the most common mesenchymal neoplasm of horses worldwide [3, 10]. Understanding the molecular underpinnings of BPV infection is not merely an academic exercise; it is essential for developing targeted prophylactic and therapeutic interventions, as recognized by the World Organisation for Animal Health (WOAH) in its classification of BPV-associated diseases as economically impactful conditions requiring surveillance and control.

Viral Genome Organization and the Oncogenic Arsenal

The BPV genome, approximately 7.2–8.0 kilobase pairs in size, is a circular episome that replicates within the nucleus of infected cells. The genome is functionally divided into three major regions: the early (E) region, encoding non-structural regulatory and transforming proteins; the late (L) region, encoding the L1 major and L2 minor capsid proteins; and the upstream regulatory region (URR), which contains cis-acting elements essential for viral transcription and replication [2, 8]. The early region of Deltapapillomavirus species, including BPV-1 and BPV-2, encodes five canonical early proteins, E1, E2, E5, E6, and E7, with E8 also identified in certain genotypes such as BPV-15 [8]. Unlike high-risk human papillomaviruses (HPVs), where E6 and E7 are the primary oncogenes, BPV-induced transformation is predominantly driven by the E5 oncoprotein, a small, highly hydrophobic transmembrane protein that functions as a potent growth factor receptor activator.

The E5 oncoprotein of BPV-1 and BPV-2 is a 44-amino-acid polypeptide that localizes to the Golgi apparatus, endoplasmic reticulum, and plasma membrane. Its primary mechanism of action involves the ligand-independent activation of the platelet-derived growth factor β receptor (PDGFβR) and the epidermal growth factor receptor (EGFR). By forming a stable complex with the 16-kDa subunit of the vacuolar H⁺-ATPase, E5 prevents receptor degradation and promotes sustained mitogenic signaling through the Ras-MAPK and PI3K-Akt pathways. This constitutive activation drives cellular proliferation and is considered the cornerstone of BPV-induced fibropapilloma formation. Critically, the E5 protein of BPV-2 and BPV-13 has been detected in urothelial cells of cattle with naturally occurring urinary bladder tumors, providing direct evidence of its role in abortive infections that lead to neoplastic transformation at sites distant from the primary cutaneous lesion [6, 18]. The ability of BPV to establish such abortive infections, where viral DNA persists but does not complete the productive life cycle, is a hallmark of its pathogenic strategy and is particularly relevant to the development of urinary bladder cancer in cattle grazing on bracken fern-infested pastures, a classic example of viral–environmental co-carcinogenesis.

Host Cell Entry, Replication, and the Role of the L1 Capsid Protein

The infectious cycle of BPV begins with viral attachment to host cells, a process mediated by the L1 major capsid protein. The L1 protein self-assembles into pentameric capsomeres that form the icosahedral capsid, and its interaction with cell surface heparan sulfate proteoglycans (HSPGs) constitutes the initial binding step. Following attachment, the virion undergoes conformational changes that facilitate entry via clathrin-mediated endocytosis or caveolar uptake, depending on the cell type. The L2 minor capsid protein then orchestrates the delivery of the viral genome to the nucleus, where it establishes as a low-copy-number episome. The L1 protein is not only critical for infectivity but also serves as the primary target for neutralizing antibodies, making it the cornerstone of vaccine development. Recombinant BPV-1 L1 protein expressed in baculovirus systems self-assembles into virus-like particles (VLPs) that are morphologically indistinguishable from native virions and elicit robust humoral immune responses in murine models [5]. Furthermore, the production of BPV-1 pseudovirions (PsVs) in plant-based expression systems, such as Nicotiana benthamiana, has demonstrated the feasibility of producing functional PsVs capable of infecting mammalian cells and expressing encapsidated reporter genes, opening avenues for plant-derived veterinary vaccines [9].

The genetic diversity of the L1 gene has profound implications for viral pathogenesis and immune evasion. Sequence analysis of BPV-1 isolates from Iraq revealed two non-synonymous mutations, SER31/ASN and ALA55/ASP, located within the surface-exposed loops of the L1 protein, with the ALA55/ASP substitution mapping to Loop 1, a region critical for antibody recognition [4]. These mutations, while not altering the overall capsid architecture, may influence antigenic variability and the ability of the virus to escape pre-existing immunity. In equine sarcoids, comparative sequencing of the L1 gene from BPV-1 isolates infecting horses, donkeys, and cattle demonstrated 98.5% nucleotide identity across species, with non-conservative mutations (T348N and K351T) identified in the equine-derived viruses [10]. This high degree of conservation underscores the cross-species transmission potential of BPV-1 and its remarkable ability to adapt to heterologous hosts, a feature that complicates disease control in multi-species environments.

The E5 Oncoprotein and Dysregulation of Cellular Signaling

The transformative capacity of BPV E5 extends beyond growth factor receptor activation. In urothelial cells of cattle infected with BPV-2 and BPV-13, the E5 oncoprotein has been shown to induce profound mitochondrial dysfunction, triggering a cascade of events that culminate in parkin-mediated mitophagy [6]. This process involves the inner mitochondrial membrane protein prohibitin 2 (PHB2), which acts as a receptor for the PINK1/parkin mitophagy pathway. In BPV-infected urothelial cells, PHB2 is significantly overexpressed in mitochondrial fractions and forms functional complexes with PHB1, phosphorylated PINK1, parkin, and the autophagosomal marker LC3-II. The interaction between PHB2 and transcription factor EB (TFEB), a master regulator of lysosomal biogenesis, suggests that BPV infection activates a coordinated transcriptional program that enhances autophagic flux. Furthermore, the interaction of PHB2 with embryonic stem cell-expressed Ras (ERAS), a constitutively active Ras family member, implicates ERAS in the maturation of mitophagosomes. This BPV-driven mitophagy is not merely a cellular stress response; it represents a viral strategy to maintain mitochondrial quality control while preventing the release of pro-apoptotic factors, thereby ensuring host cell survival and persistent infection.

The expression of the innate immune receptor Mincle (macrophage-inducible C-type lectin) in urothelial cancer cells of BPV-associated bovine urinary bladder tumors adds another layer of complexity to host–virus interactions [18]. Mincle, which recognizes the nuclear protein SAP130 released from dying cells, is expressed in a membranous pattern in papillary urothelial cancers and shows strong cytoplasmic immunoreactivity in invasive cancers. The functional significance of Mincle expression in neoplastic cells is twofold: it may enable cancer cells to act as antigen-presenting cells, potentially modulating anti-tumor immune responses, and it may influence the efficacy of Bacillus Calmette-Guérin (BCG) immunotherapy, which is the standard treatment for non-muscle invasive bladder cancer in humans. This finding has direct translational relevance, as BCG immunotherapy is also employed in equine sarcoid treatment, albeit with variable outcomes [12]. The presence of Mincle on BPV-transformed cells could theoretically enhance the immunostimulatory effects of BCG, providing a mechanistic rationale for its use in veterinary oncology.

Immune Evasion and the Tumor Microenvironment

BPV has evolved multiple strategies to subvert host immune surveillance, a prerequisite for establishing persistent infections that can progress to neoplasia. The viral life cycle is tightly linked to the differentiation state of the host keratinocyte, with productive replication occurring only in terminally differentiated cells of the stratum granulosum and stratum corneum. This differentiation-dependent replication minimizes exposure to immune recognition, as the virus remains largely intracellular and non-lytic during the early stages of infection. However, the abortive infections characteristic of BPV-2 and BPV-13 in urothelial cells present a different immunological challenge, as the virus persists in a non-productive state that may evade detection by pattern recognition receptors.

The E5 oncoprotein itself contributes to immune evasion by downregulating major histocompatibility complex class I (MHC-I) expression on the surface of infected cells, thereby impairing CD8⁺ T-cell-mediated cytotoxicity. Additionally, BPV infection modulates the expression of cytokines and chemokines within the tumor microenvironment. In equine sarcoids, the presence of BPV-1 and BPV-2 DNA is associated with a fibroproliferative response characterized by the activation of cancer-associated fibroblasts and the deposition of extracellular matrix components. The matrix metalloproteinases (MMPs), particularly MMP-1, MMP-2, and MMP-9, play a critical role in sarcoid invasiveness, and their expression is modulated by viral oncoproteins. Recent in vitro studies have demonstrated that cannabidiol (CBD) treatment of primary equine sarcoid cells significantly reduces MMP-1 and MMP-2 concentrations, with MMP-9 also showing decreased levels, suggesting that targeting the extracellular matrix remodeling machinery may represent a novel therapeutic avenue [17].

Co-infections and the Role of the Viral Microbiome

The concept of multiple BPV infections in a single host has gained traction as molecular detection methods have improved. In Brazilian cattle herds, analysis of skin warts from six animals using FAP primer-based PCR and sequencing revealed the presence of up to three different BPV types in individual animals, including BPV-1, -2, -6, and -8, as well as three putative new types [15]. This high prevalence of mixed infections suggests that co-infection may be the rule rather than the exception, and it raises important questions about viral interference, synergy, and the potential for recombination events that could generate novel pathogenic variants. The detection of BPV-1 in 84.6% of cutaneous wart samples from Egyptian cattle, with concurrent absence of BPV-2, -4, -5, and -10, indicates that type distribution varies geographically and may be influenced by host genetics, management practices, and environmental factors [7].

The role of the host genetic background in susceptibility to BPV infection and disease progression is an area of active investigation. Full genome sequencing of BPV-1 isolates from Chinese native cattle breeds, which are known for their higher disease resistance, revealed that these isolates belong to subtype A and are phylogenetically closer to the common ancestor than European isolates, which form a more recently diverged group (subtype B) [11]. This suggests that co-evolution between BPV and its bovine host has shaped viral diversity, with certain lineages being better adapted to specific host populations. The identification of a single nucleotide insertion in the upstream regulatory region of BPV-1 from a calf with severe cutaneous papillomatosis further highlights the potential for minor genetic changes to have major phenotypic consequences, as alterations in the URR can affect transcription factor binding and viral gene expression [2].

Cross-Species Transmission and the Equine Sarcoid Paradigm

The transmission of BPV-1 and BPV-2 from cattle to horses, resulting in equine sarcoids, represents one of the most compelling examples of cross-species viral oncogenesis in veterinary medicine. The molecular mechanisms underlying this host jump are not fully understood, but they likely involve the ability of the BPV E5 protein to activate equine PDGFβR and EGFR, despite species-specific differences in receptor structure. The viral load in equine sarcoids varies considerably between lesions, and accurate quantification is essential for both diagnostic and research purposes. Fine-needle aspirates (FNA) have been shown to provide a more accurate approximation of the actual viral load than superficial swabs, with a strong correlation (r = 0.59) between FNA results and postoperative tissue biopsies in occult sarcoids [3]. This methodological refinement is critical for studies investigating the relationship between viral load and clinical outcome, as well as for monitoring the response to therapy.

The phylogenetic analysis of BPV-1 isolates from equine sarcoids in Egypt revealed that viruses from horses, donkeys, and cattle cluster together without species-specific differentiation, indicating frequent cross-species transmission events [10]. The high degree of sequence identity across all genes except E5 and L2, where substitutions of up to 5.18% (nucleotide) and 6.81% (amino acid) were observed, suggests that these two genes may be under selective pressure to adapt to the equine host. The E5 gene, in particular, showed the highest variability, which may reflect the need for the oncoprotein to optimally interact with equine growth factor receptors. This genetic plasticity underscores the potential for BPV to continue evolving and potentially expanding its host range, a concern that is particularly relevant in regions where cattle and equids are managed in close proximity.

Viral Pathogenesis in the Context of the Bovine Papillomavirus Life Cycle

The complete genome sequencing of novel BPV types, such as BPV-15 from dairy cows in southern Xinjiang, China, has expanded our understanding of the genetic diversity within the Xipapillomavirus genus [8]. BPV-15, which encodes E8, E7, E1, E2, E4, L1, and L2, but lacks the E5 and E6 genes typical of Deltapapillomavirus, provides a contrasting model of viral pathogenesis. The E7 protein of BPV-15 contains a consensus CX₂CX₂₉CX₂C zinc-binding domain and an LxCxE motif, which is known to mediate binding to the retinoblastoma (pRb) tumor suppressor protein. This suggests that BPV-15, like high-risk HPVs, may transform cells through the inactivation of pRb, rather than through the E5-driven growth factor receptor activation that characterizes BPV-1 and BPV-2. The existence of such divergent pathogenic mechanisms within the same viral family highlights the importance of genotype-specific approaches to diagnosis, treatment, and prevention.

The detection of BPV DNA in the upper gastrointestinal tract of a calf with severe papillomatosis, where the lesions were diagnosed as fibropapillomas lacking BPV antigen expression, illustrates the phenomenon of abortive infection at mucosal sites [2]. The presence of identical BPV-1 genomes in both cutaneous papillomas and gastrointestinal fibropapillomas from the same animal suggests that the virus can disseminate hematogenously or via the lymphatic system, establishing secondary sites of infection where the viral life cycle is arrested at an early stage. This abortive infection is characterized by the expression of early genes, particularly E5, in the absence of late gene expression and virion production, and it is this non-productive state that is most strongly associated with neoplastic progression.

Concluding Perspectives on Molecular Pathogenesis

The molecular pathogenesis of BPV infection is a multifaceted process that integrates viral oncogene function, host cell signaling dysregulation, immune evasion, and environmental co-factors. The E5 oncoprotein of Deltapapillomavirus species stands as a master regulator of cellular transformation, driving proliferation through growth factor receptor activation and promoting cell survival through the induction of mitophagy. The abortive infections that characterize BPV-2 and BPV-13 in urothelial cells represent a unique pathogenic niche where the virus persists in a non-productive state, yet continues to exert oncogenic effects through the sustained expression of early proteins. The cross-species transmission of BPV to equids, resulting in sarcoids, underscores the zoonotic and transboundary nature of this pathogen and highlights the need for integrated surveillance and control strategies that align with WOAH guidelines. As our understanding of BPV pathogenesis deepens, the development of L1-based VLP vaccines, plant-derived pseudovirions, and targeted immunotherapies offers hope for reducing the economic burden of BPV-associated diseases in cattle and horses worldwide.

Epidemiological Patterns and Global Distribution of Bovine Papillomavirus Strains

Global Phylodynamic Diversity and Taxonomic Considerations

Bovine papillomaviruses (BPVs) constitute a remarkably heterogeneous assemblage of circular double-stranded DNA viruses belonging to the family Papillomaviridae, with at least 15 formally recognized types distributed across multiple genera, including Deltapapillomavirus, Xipapillomavirus, and Epsilonpapillomavirus. The epidemiological landscape of BPV infection is defined by a complex interplay of viral genotype, host species susceptibility, anatomical tropism, and geographic partitioning. Among the most clinically and economically significant are the members of the Deltapapillomavirus genus, BPV types 1, 2, and 13, which exhibit a unique capacity for cross-species transmission, most notably to equids, where they serve as the primary etiological agents of equine sarcoids [3, 10]. The global distribution of these viruses is not uniform; rather, it reflects both historical patterns of cattle domestication and trade, as well as contemporary husbandry practices that facilitate viral maintenance and dissemination.

The molecular characterization of BPV genomes has revealed substantial intratypic diversity that often corresponds to geographic provenance. Phylogenetic analyses of BPV-1 isolates from China, for instance, have demonstrated that strains recovered from indigenous Chinese native cattle breeds cluster within a distinct ancestral lineage, designated subtype A, which is phylogenetically basal to the European isolates that form the more recently diverged group B [11]. This pattern suggests that BPV-1 likely originated in Asia and subsequently spread to Europe through the translocation of infected animals, a hypothesis supported by the observation that the Chinese isolates share a closer genetic relationship with the common ancestor of all BPV-1 strains [11]. Similarly, the first identification of BPV-2 in Egypt revealed that the circulating strain shared 99.6% nucleotide identity with reference strains from Brazil and China, indicating a high degree of sequence conservation across continents and implying either recent introduction or strong purifying selection on the L1 gene [1]. Such findings underscore the importance of continuous molecular surveillance to differentiate between endemic circulation and novel incursions.

Regional Epidemiological Dissections

Africa and the Middle East: The epidemiological profile of BPV in North Africa and the Middle East is characterized by a predominance of BPV-1, with recent studies documenting its widespread presence across Egypt and Iraq. In a targeted molecular survey of 308 cattle from Egypt’s New Valley Province, BPV-1 DNA was detected in 84.6% of cutaneous wart biopsies, while BPV-2, -4, -5, and -10 were conspicuously absent [7]. This striking predominance of a single type suggests either a founder effect, where BPV-1 was introduced and subsequently amplified within the population, or a competitive advantage of BPV-1 over other types in the Egyptian ecological niche. Phylogenetic analyses further revealed that the Egyptian BPV-1 strains are closely related to those recovered from buffalo in India, pointing to a potential historical route of viral dissemination along ancient or modern trade corridors [7]. The recent identification of BPV-2 in Egypt, the first such report in the country, from bovines presenting with generalized papillomatosis in the Menya and Ismailia provinces indicates that the BPV landscape is dynamic, with new types continuously emerging or being introduced [1]. In Iraq, an exhaustive molecular survey of 50 tumor samples from Babylon, Wasit, and Al-Qadisiyah provinces reported an 84% BPV DNA positivity rate, with sequence analysis of the complete L1 gene confirming BPV-1 as the predominant circulating type. Notably, the Iraqi BPV-1 strains exhibited amino acid substitutions at positions SER31/ASN and Ala55/ASP within the L1 protein’s surface-exposed loops, mutations that may have implications for antigenicity and vaccine design [4]. The presence of these mutations in geographically distinct herds suggests they may represent stable polymorphisms characteristic of Middle Eastern BPV-1 lineages.

Asia and the Far East: China, as one of the world’s largest beef and dairy producers, has experienced significant BPV-associated economic losses. The complete genomic characterization of BPV-1 isolates from Chinese native cattle breeds (specifically the JX180408, LA150909, HX160815, and BS160810 strains) has provided critical insights into the evolutionary history of the virus in East Asia. These isolates were unambiguously classified within the Deltapapillomavirus genus and were found to belong to an ancestral subtype A, which is phylogenetically distinct from the European subtype B [11]. This deep phylogenetic divergence suggests that BPV-1 has been enzootic in Chinese cattle for a considerable period, evolving independently from European lineages following the historical separation of cattle populations. Furthermore, the complete genome of a novel putative type, BPV-15 (strain BPV Aks-02), was determined from a cutaneous neoplastic lesion on a dairy cow in Southern Xinjiang, China. This virus, clustering within the Xipapillomavirus genus, possessed a 7189-base-pair genome encoding five early genes (E8, E7, E1, E2, E4) and two late genes (L1, L2), with the L1 open reading frame sharing 99% nucleotide identity with the putative type BAPV-3 [8]. The detection of BPV-15 adds to the growing list of xipapillomaviruses associated with bovine cutaneous lesions and highlights the underappreciated diversity of BPVs in Central Asian cattle populations.

Europe and the Americas: In Brazil, a seminal study employing the FAP primer system for partial L1 gene amplification followed by cloning and sequencing provided the first robust evidence of multiple BPV infections within individual animals. Among six bovines from three herds in Paraná state, four characterized BPV types (BPV-1, -2, -6, and -8) and three previously described putative new types (BPV/BR-UEL3, BPV/BR-UEL4, and BPV/BR-UEL5) were identified, with double infections documented in four of the six animals [15]. This work unequivocally demonstrated that mixed infections with two or three different BPV types are as frequent in cattle as they are in humans, a finding that has profound implications for understanding viral interference, recombination potential, and the evolution of pathogenicity. The presence of high viral diversity within a single herd, six different types in one farm, suggests that co-infection is a common phenomenon that may facilitate the emergence of novel recombinant strains [15]. At the molecular level, the E5 oncoprotein, particularly from BPV-2 and BPV-13, has been shown to play a pivotal role in abortive infections in urothelial cells, triggering parkin-dependent mitophagy through interactions with prohibitin 2 (PHB2) and the PINK1/parkin pathway [6]. This mechanistic understanding of BPV pathogenesis is critical for interpreting epidemiological patterns, as the ability of certain types to establish persistent, non-productive infections in the urinary bladder may contribute to their maintenance in cattle populations and their association with urinary bladder neoplasia, which has been documented in regions with high exposure to bracken fern.

Cross-Species Transmission and Ecological Drivers

The epidemiological significance of BPV extends far beyond its bovine host. Deltapapillomavirus types 1 and 2 are the causative agents of equine sarcoids, the most common mesenchymal neoplasm of horses worldwide, and affect not only domestic horses but also donkeys, mules, and wild equids. In a molecular investigation of sarcoids in Egypt, BPV-1 was identified in 14 of 25 equids (eight horses and six donkeys) presented at the Zagazig University Veterinary Clinic. Remarkably, whole-genome sequencing via Illumina technology revealed that the papillomaviruses recovered from donkeys and horses were genetically identical, sharing 98.5% identity with the closest reference sequence (KX907623) and demonstrating that viral strains can circulate seamlessly between equine species without marked phylogenetic segregation [10]. Furthermore, comparison of the L1 gene across BPV-1 strains from cattle, horses, and donkeys identified only minor substitutions, including two non-conservative mutations (T348N and K351T) that did not produce species-specific clustering [10]. This lack of host-associated phylogenetic differentiation strongly supports the hypothesis that BPV-1 is a generalist pathogen capable of crossing species barriers with relative ease, likely facilitated by direct contact between cattle and equids or through arthropod vectors. The World Organisation for Animal Health (WOAH) recognizes equine sarcoids as a significant welfare and economic concern in many regions, particularly where working equids are integral to agricultural livelihoods.

The role of arthropod vectors in BPV transmission merits careful consideration from an epidemiological perspective. The mechanical transmission of BPV by flies, particularly Haematobia irritans (horn flies) and Stomoxys calcitrans (stable flies), has been implicated in the rapid spread of papillomatosis within herds and between species. Cutaneous lesions, which often occur on the head, neck, and trunk, are common sites for fly feeding, and the presence of ticks has also been noted in clinical cases of bovine papillomatosis [13, 14]. The spatial clustering of sarcoid cases in horses grazing alongside BPV-infected cattle further supports vector-borne transmission, and the WHO has highlighted the importance of vector control in managing viral diseases of livestock. The economic impact of BPV-associated disease is substantial: in cattle, papillomatosis can lead to reduced weight gain, decreased milk production, hide damage, and, in severe cases, astasia necessitating euthanasia [2]. In equids, sarcoids impair working ability, cause cosmetic disfigurement, and are notoriously difficult to treat, with high recurrence rates following conventional therapies [12].

Temporal Trends and Emerging Threats

Longitudinal data on BPV prevalence are scarce, but recent molecular surveys indicate that the landscape of BPV infection is changing, with the emergence of new types and the expansion of previously restricted genotypes into new geographic regions. The detection of BPV-2 for the first time in Egypt in 2021, despite the prior documentation of BPV-1 in the country, suggests that the viral population is not static and that international livestock trade may be introducing novel strains [1, 7]. Similarly, the identification of BPV-15 in China adds to the growing roster of xipapillomaviruses that may be associated with distinct clinical outcomes [8]. The increasing application of next-generation sequencing technologies in veterinary diagnostics is likely to reveal even greater diversity, as evidenced by the discovery of multiple putative new types in Brazilian cattle [15]. The Food and Agriculture Organization (FAO) has emphasized the need for enhanced surveillance of livestock viruses in the context of globalized trade, as the movement of live animals continues to erode geographic barriers to pathogen dissemination.

From a clinical and pathological perspective, the severity of BPV-induced disease varies markedly with viral type and host immune status. Severe cutaneous multiple papillomatosis in calves, as documented in a case involving BPV-1, can result in the development of hundreds of papillomas on the skin and nodules in the upper gastrointestinal tract, leading to euthanasia [2]. The full genome characterization of this BPV-1 strain revealed a single nucleotide insertion in the upstream regulatory region relative to the reference sequence, a mutation that may alter viral gene expression and account for the unusually aggressive phenotype [2]. The interplay between viral genetics and host factors is further underscored by the observation that BPV-2 and BPV-13 have been implicated in the pathogenesis of urothelial tumors in cattle, where the E5 oncoprotein drives abortive infections characterized by mitochondrial dysfunction and dysregulated mitophagy [6]. The expression of the innate immune receptor Mincle (macrophage-inducible C-type lectin) in BPV-associated urothelial cancer cells, but not in normal urothelium, suggests that the virus may modulate the tumor microenvironment in ways that impact disease progression and host immune surveillance [18]. These molecular insights are essential for understanding why certain BPV types are associated with particular disease syndromes and why their global distribution patterns correlate with the prevalence of specific clinical outcomes.

Clinical Manifestations and Pathological Spectrum: From Cutaneous Papillomas to Gastrointestinal Fibropapillomas

Bovine papillomavirus (BPV) infection manifests across a remarkably broad clinicopathological continuum, ranging from benign, self-limiting cutaneous papillomas to more complex fibropapillomatous proliferations that can involve both the integumentary and alimentary systems. Understanding this spectrum is critical for veterinary practitioners, as the clinical presentation, histopathological features, and biological behavior of lesions vary considerably depending on viral genotype, host immune status, anatomical site, and the presence of co-infections. The World Organisation for Animal Health (WOAH) recognizes bovine papillomatosis as a condition of significant economic consequence, particularly in dairy and beef production systems where the disease can impair animal welfare, reduce productivity, and necessitate costly interventions [1, 7]. The Food and Agriculture Organization (FAO) has similarly emphasized the importance of understanding the epidemiological and pathological dimensions of BPV infections to inform control strategies in resource-limited settings.

Cutaneous Papillomas: The Classical Presentation

The most frequently encountered clinical manifestation of BPV infection is the cutaneous papilloma, commonly referred to as a wart or verruca. These lesions are typically exophytic, cauliflower-like proliferations that can range in size from a few millimeters to several centimeters in diameter [1, 13, 14, 20]. In affected cattle, papillomas may appear as solitary growths or, more commonly, as multiple, coalescing masses distributed across the head, neck, shoulders, and trunk [7, 11, 15]. The anatomical distribution is not random; trauma and abrasions, often facilitated by arthropod vectors such as ticks, provide portals of entry for the virus, explaining the predilection for areas prone to mechanical injury [13, 14]. Clinically, these lesions are firm, hyperkeratotic, and often pedunculated, with a rough, fissured surface that may accumulate dirt and debris. In advanced or heavily colonized cases, an offensive odor and secondary fly infestation can ensue, compounding animal distress [14].

The natural history of cutaneous papillomas follows a predictable pattern of growth, maturation, and, in immunocompetent hosts, spontaneous regression. Following an incubation period of several weeks to months, lesions initially appear as small, smooth papules that gradually enlarge and become hyperkeratotic. The proliferative phase is driven by viral replication within the keratinocytes of the stratum spinosum and granulosum, where the viral genome is maintained as an episome and early gene expression, particularly the E5, E6, and E7 oncoproteins, disrupts normal cell cycle regulation. Regression typically occurs within 6 to 12 months and is mediated by a cell-mediated immune response directed against viral antigens. This phenomenon has been well documented in both natural infections and experimental models, where biopsies or mechanical disruption of lesions can accelerate the regression process [13, 19]. The mechanisms underlying regression involve infiltration of CD8+ cytotoxic T lymphocytes and activation of innate immune pathways, though the precise molecular triggers remain an active area of investigation.

Histopathological Hallmarks of Cutaneous Papillomatosis

Histopathological examination of BPV-induced cutaneous papillomas reveals a constellation of characteristic features that reflect the virus's tropism for epithelial cells and its capacity to induce both hyperplastic and dysplastic changes. The most consistent findings include marked hyperkeratosis (both orthokeratotic and parakeratotic forms), acanthosis (thickening of the stratum spinosum), and papillomatosis (projections of fibrovascular cores covered by hyperplastic epithelium) [15, 20]. Within the granular layer, hypergranulosis is commonly observed, accompanied by the presence of koilocytes, enlarged keratinocytes with perinuclear halos and wrinkled, pyknotic nuclei that are pathognomonic for papillomavirus infection. These koilocytes are the cytopathic hallmark of productive viral replication and are most abundant in the superficial layers of the epidermis, where viral assembly and release occur.

Immunohistochemical (IHC) staining using anti-BPV antibodies consistently demonstrates nuclear positivity within the granular and, to a lesser extent, basal layers of the epidermis [20]. This pattern of antigen distribution correlates with the site of active viral transcription and capsid protein synthesis. Additionally, p53 tumor suppressor protein is frequently overexpressed in papilloma tissues, appearing as strong cytoplasmic and perinuclear staining predominantly in the basal and parabasal layers [20]. The dysregulation of p53, likely a consequence of E6-mediated degradation pathways being overwhelmed or functionally altered, underscores the early oncogenic potential of BPV even in benign lesions. It is important to note that while these histopathological features are highly suggestive of BPV etiology, definitive diagnosis requires molecular confirmation, as similar proliferative lesions can arise from other etiologies [15, 20].

Variability in Cutaneous Presentation: Generalized Papillomatosis and Mixed Infections

While solitary or few papillomas are common, a subset of animals, particularly young calves and immunocompromised individuals, may develop generalized papillomatosis, characterized by hundreds of lesions distributed across the body surface [1, 2]. This severe form can lead to significant morbidity, including astasia (inability to stand) that necessitates euthanasia in extreme cases [2]. The underlying virological basis for such widespread dissemination often involves a combination of high viral load, multiple BPV genotypes infecting the same host, and deficiencies in cell-mediated immunity. Source [15] provides compelling evidence that mixed infections with two or more BPV types (e.g., BPV-1, -2, -6, -8, and several putative novel types) are common in Brazilian cattle herds, mirroring the complexity seen in human papillomavirus coinfections. The clinical significance of these mixed infections is that they may broaden the spectrum of tissue tropism, alter the kinetics of lesion development, and potentially influence the likelihood of malignant transformation.

The macroscopic appearance of lesions can also vary depending on the specific BPV genotype involved. Xipapillomaviruses, such as BPV-15, are known to produce purely epithelial papillomas with minimal dermal involvement, whereas Deltapapillomaviruses (BPV-1, -2, -13) frequently induce fibropapillomas, mixed tumors featuring both epithelial hyperplasia and robust fibroblastic proliferation of the underlying dermis [2, 8, 11]. This distinction is of more than academic interest; fibropapillomas are more locally invasive, have a higher recurrence rate following surgical excision, and are associated with a more complex interplay between viral oncoproteins and host cellular signaling pathways.

Gastrointestinal and Internal Organ Involvement: The Fibropapilloma Spectrum

Beyond the cutaneous surface, BPV, particularly BPV-1 and BPV-2, can induce fibropapillomatous lesions in the upper gastrointestinal tract, including the oral cavity, esophagus, and forestomachs [2, 11]. This internal manifestation is frequently underdiagnosed because lesions are not externally visible and may only be detected at necropsy or during endoscopic examination. Source [2] describes a calf with severe cutaneous papillomatosis that also harbored multiple nodules in the upper gastrointestinal tract; histopathologically, these were diagnosed as fibropapillomas. Critically, immunohistochemistry failed to detect BPV antigen within the gastrointestinal lesions, despite the presence of identical viral genomes to those found in the skin papillomas. This suggests that in the gastrointestinal microenvironment, the viral life cycle may be restricted, with viral DNA persisting in a latent or abortive state without productive replication, yet the early oncoproteins (particularly E5) continue to drive mesenchymal proliferation.

The gastrointestinal fibropapillomas represent a distinct pathological entity from their cutaneous counterparts. Grossly, they appear as firm, sessile or pedunculated nodules that may protrude into the lumen of the esophagus or rumen. Histologically, they are characterized by a predominant fibroblastic component, with attenuated epithelial covering that may show only mild hyperplasia. The connective tissue proliferation is composed of spindle-shaped cells arranged in interlacing fascicles, often with varying degrees of collagen deposition. This fibroproliferative response is directly attributable to the paracrine effects of BPV E5 oncoprotein, which constitutively activates the platelet-derived growth factor beta receptor (PDGFβR) and other receptor tyrosine kinases, driving uncontrolled fibroblast proliferation [6, 18].

Of particular concern is the association between BPV-2 and BPV-13 infection and urothelial tumors of the urinary bladder, including transitional cell carcinoma [6, 11, 18]. While not strictly a fibropapilloma, these urothelial neoplasms share a common etiological thread with gastrointestinal fibropapillomas: both arise from abortive BPV infections in which the virus does not complete its replicative cycle but persists in a transcriptionally active state. Source [6] elegantly demonstrates that BPV E5 oncoprotein is expressed in urothelial cells of cattle with bladder cancer and is associated with profound mitochondrial dysfunction and activation of parkin-mediated mitophagy. This pathway, involving prohibitin 2 (PHB2) overexpression and interaction with PINK1/parkin and LC3-II, represents a novel mechanism by which BPV contributes to cellular transformation and survival under metabolic stress. Furthermore, source [18] reports the expression of the innate immune receptor Mincle in urothelial cancer cells from BPV-associated bladder tumors, suggesting that these neoplastic cells may acquire antigen-presenting capabilities that could influence both tumor progression and host immune surveillance.

The clinical implications of gastrointestinal and bladder involvement are substantial. Unlike cutaneous papillomas, which frequently regress spontaneously, fibropapillomas of the gastrointestinal tract and urothelial tumors tend to persist and may progress to malignancy. Esophageal and ruminal fibropapillomas can lead to dysphagia, ruminal tympany, and weight loss, while bladder tumors may cause hematuria, stranguria, and secondary urinary tract infections. The economic impact is further compounded by the fact that affected animals are often condemned at slaughter due to visceral lesions, representing a total loss of the animal's value [11, 16]. The WOAH has classified these BPV-associated neoplasms as emerging threats in certain production systems, particularly where bracken fern (Pteridium aquilinum) exposure acts as a co-carcinogen by immunosuppressing the host and activating latent BPV infections.

Diagnostic Approaches: Electron Microscopy, Molecular Characterization, Histopathology, and Immunohistochemistry

The accurate and definitive diagnosis of bovine papillomavirus (BPV) infection is a multifaceted endeavor that relies on a synergistic integration of virological, molecular, and histopathological techniques. Given the diverse clinical presentations of BPV, ranging from benign cutaneous papillomas to malignant urinary bladder neoplasms and the cross-species pathogenesis observed in equine sarcoids, a comprehensive diagnostic algorithm is essential for both clinical management and epidemiological surveillance. The World Organisation for Animal Health (WOAH) recognizes the economic impact of BPV and the need for standardized diagnostic protocols. This section provides an exhaustive analysis of the core diagnostic pillars: electron microscopy, molecular characterization, histopathology, and immunohistochemistry, detailing their individual strengths, inherent limitations, and synergistic applications in the context of BPV infection.

Electron Microscopy: Direct Visualization of Viral Particles

Electron microscopy (EM) represents the classical gold standard for the direct visualization of viral particles, offering a morphology-first approach that requires no prior knowledge of the viral genome. The diagnostic utility of EM in BPV infection is predicated on the characteristic icosahedral structure of the papillomavirus capsid, approximately 55-60 nm in diameter, which is readily identifiable in negatively stained preparations of wart homogenates or tissue extracts [1]. This technique was instrumental in the initial discovery and classification of papillomaviruses and remains a powerful tool for confirming the presence of a papovavirus-like agent in clinical samples. In a study characterizing a novel BPV-2 strain in Egypt, transmission electron microscopy (TEM) of clarified wart homogenates definitively revealed the presence of non-enveloped, icosahedral viral particles, confirming the papovavirus morphology and providing the initial virological evidence before molecular characterization was undertaken [1].

Beyond simple detection, EM has advanced to provide ultrastructural insights into virus-host cell interactions and viral morphogenesis. The in planta production of BPV-1 pseudovirions (PsVs) and virus-like particles (VLPs) for vaccine development relies critically on TEM to confirm the structural integrity and correct assembly of the capsid proteins into the requisite icosahedral shells [5, 9]. Without this morphological validation, the functionality of these immunogens cannot be assured. Furthermore, advanced applications of TEM have elucidated the pathobiological mechanisms of BPV infection. For instance, in BPV-2/-13-associated urothelial carcinogenesis, TEM has revealed severe ultrastructural abnormalities of the inner mitochondrial membrane (IMM) in infected cells, providing direct visual evidence of mitochondrial dysfunction linked to the expression of the E5 oncoprotein [6]. These observations of disrupted cristae architecture are foundational to understanding the role of parkin-mediated mitophagy in BPV-driven oncogenesis, a finding that would be impossible to ascertain with molecular techniques alone [6]. While EM is not a practical tool for high-throughput screening or field diagnostics due to its cost, technical expertise, and relatively low sensitivity compared to molecular amplification, its unparalleled ability to provide direct visual proof of viral particles and detailed subcellular pathology ensures its continued relevance in confirmatory diagnostics and advanced pathogenesis research.

Molecular Characterization: The Cornerstone of Genotyping and Quantification

Molecular techniques, particularly nucleic acid amplification technologies (NAATs), have become the cornerstone of BPV diagnostics, offering unparalleled sensitivity, specificity, and the capacity for precise genotyping and viral load quantification. The World Organisation for Animal Health (WOAH) recommends PCR-based methods for the confirmatory diagnosis of BPV and for epidemiological studies. The fundamental approach involves the extraction of total DNA from lesional tissue, followed by PCR amplification of conserved genomic regions, most commonly the L1 gene, which is the most phylogenetically informative [1, 4, 15].

Conventional PCR and Genotyping: The use of consensus or degenerate primers, such as the FAP primer pair targeting the L1 gene, allows for the detection of a broad spectrum of known and novel BPV types. This approach has been instrumental in revealing the high diversity and frequency of multiple BPV infections within single animals and herds. For instance, in Brazilian cattle herds, PCR with FAP primers followed by cloning and sequencing of amplicons identified four characterized BPV types (BPV-1, -2, -6, -8) and three previously undescribed putative new types, demonstrating that mixed infections with two or three different BPV types are common [15]. This technique is critical for understanding the complex epidemiology of BPV, as the pathogenic potential and tissue tropism can vary significantly between genotypes. Species-specific or type-specific PCR assays, such as those designed for BPV-1, -2, -4, -5, and -10, provide a more targeted approach for clinical diagnostics in known endemic regions [7].

Quantitative Real-Time PCR (qPCR) and Viral Load: The advent of qPCR has revolutionized the study of BPV pathogenesis, particularly in the context of equine sarcoids, where the presence and quantity of BPV-1/-2 DNA are causally linked to disease. qPCR provides an accurate and sensitive measure of viral load (VL), which is essential for stratifying lesions, monitoring disease progression, and assessing therapeutic responses. Research has demonstrated that superficial swab sampling, while convenient, significantly underestimates the true VL in sarcoids compared to fine-needle aspirates (FNA) or tissue biopsies. A pivotal study by Gysens et al. established that qPCR of FNA samples showed a strong correlation with the reference standard of postoperative tissue biopsies, particularly for occult sarcoids (r=0.59), identifying FNA as the superior sampling method for accurate VL quantification [3]. This finding has profound implications for clinical trials and research into BPV pathogenesis, as relying on less sensitive samples could lead to erroneous conclusions about viral dynamics and treatment efficacy.

Digital Droplet PCR (ddPCR) and Enhanced Sensitivity: For applications requiring absolute quantification without the need for standard curves and with enhanced sensitivity for low-abundance targets, digital droplet PCR (ddPCR) has emerged as a superior platform. In a landmark comparative study of papillomavirus detection, ddPCR demonstrated a significantly higher detection rate (26.6%) for ovine papillomavirus (OaPV) DNA in equine vaginal swabs compared to conventional qPCR (11.7%) [21]. This statistically significant difference (p<0.0005) underscores the ability of ddPCR to detect viral genomes present at concentrations below the reliable detection threshold of qPCR, which is particularly relevant for identifying asymptomatic carriers or resolving latent BPV infections [21]. The adoption of ddPCR in BPV research promises to improve the sensitivity of epidemiological surveys and diagnostic surveillance.

Sequencing and Phylogenetic Analysis: The ultimate confirmation of BPV genotype and strain identity is achieved through nucleotide sequencing of PCR amplicons. Sanger sequencing of the L1 gene is the standard method for genotyping, while full-genome sequencing provides a comprehensive assessment of genetic variation, including the detection of mutations in oncogenes (E5, E6, E7) and regulatory regions (URR). Phylogenetic analysis of these sequences is essential for tracing the origins and spread of BPV strains, understanding host-specificity, and identifying novel types. This approach has been critical in characterizing emerging and geographically distinct strains. For example, the first report of BPV-2 in Egypt was confirmed by sequencing the L1 gene, which showed 99.6% identity with reference strains from Brazil and China, indicating a high degree of conservation despite geographical separation [1]. Similarly, full-genome characterization of BPV type 15 (Xipapillomavirus) from Chinese dairy cows and BPV-1 from native Chinese cattle has provided crucial insights into the molecular diversity and ancestral relationships of these viruses [8, 11]. Next-generation sequencing (NGS) platforms, such as Illumina, are now being applied to BPV diagnostics, enabling the simultaneous detection of multiple pathogens in a single sample and the discovery of novel viral sequences that would be missed by targeted PCR [10, 23]. This metagenomic approach is particularly powerful for investigating complex disease syndromes like bovine respiratory disease, where co-infections are common [22, 23].

Isothermal Amplification and Point-of-Need Diagnostics: To address the challenges of deploying molecular diagnostics in resource-limited settings, isothermal amplification techniques such as Recombinase Polymerase Amplification (RPA) have been developed for BPV detection. These methods operate at a constant low temperature, eliminating the need for expensive thermal cyclers. An RPA assay targeting the BPV-1 L1 gene, combined with a lateral flow strip for visual detection of the amplicon, has been validated against conventional PCR, demonstrating comparable diagnostic performance in under one hour [7]. This point-of-need molecular test holds immense promise for field-based screening and rapid outbreak response, making confirmatory BPV diagnosis accessible in areas lacking sophisticated laboratory infrastructure.

Histopathology: The Microscopic Hallmark of Papillomavirus Infection

Histopathological examination of biopsy specimens remains an indispensable diagnostic pillar, providing a static but highly informative snapshot of the tissue-level changes induced by BPV infection. The characteristic morphology of a cutaneous papilloma or fibropapilloma serves as the initial diagnostic clue, guiding the selection of confirmatory ancillary tests. The hallmark histopathological features of a BPV-induced papilloma include a marked, exophytic proliferation of well-differentiated squamous epithelium arranged in complex, finger-like projections over a central fibrovascular core [1, 11, 20]. These projections are covered by a variably thickened layer of keratin, manifesting as orthokeratotic and/or parakeratotic hyperkeratosis [20].

A critical cytopathic effect that is highly suggestive of papillomavirus infection is the presence of koilocytes. These are large, atypical squamous cells found primarily in the upper stratum spinosum and granular layer, characterized by a shrunken, pyknotic nucleus surrounded by a prominent, clear perinuclear halo [16, 20]. This vacuolization is a direct cytopathic effect of the viral replication cycle, which disrupts the cellular cytoskeleton and causes the cytoplasm to retract. Additional features include varying degrees of acanthosis (thickening of the stratum spinosum), hypergranulosis (thickening of the granular layer), and often, an intact basal layer, which is a key feature distinguishing a benign papilloma from an invasive squamous cell carcinoma [11, 20].

In the case of fibropapillomas (e.g., BPV-1 and -2), the histopathology is distinctly different. The epithelial component may be less pronounced, and the dominant feature is a massive proliferation of dermal fibroblasts and collagenous connective tissue, forming the bulk of the tumor mass [2]. The overlying epidermis may show varying degrees of acanthosis and hyperkeratosis, but the key diagnostic feature is the dense, haphazardly arranged fibroblastic network. In urinary bladder lesions associated with BPV-2, histopathology reveals areas of chronic inflammation, epithelial hyperplasia, and in more severe cases, papillary or invasive urothelial carcinoma [18]. The evaluation of these microscopic features, including the assessment of invasion, cellular atypia, and mitotic index, is crucial for grading lesions and predicting their biological behavior, providing information that molecular diagnostics alone cannot. The presence of intranuclear inclusion bodies, though sometimes present, is not a consistent feature of BPV lesions and is more characteristic of other viral dermatitides [19, 20].

Immunohistochemistry: Detecting Viral Antigen and Host Response

Immunohistochemistry (IHC) bridges the gap between morphology and molecular identity by allowing the in situ localization of specific viral or cellular antigens within the tissue architecture. This technique is invaluable for confirming the etiology of a lesion when histopathological features are suggestive but not pathognomonic, and for studying the expression patterns of viral oncoproteins and host cellular factors in formalin-fixed, paraffin-embedded (FFPE) tissues.

Detection of BPV Structural Antigens: The most common IHC approach for diagnosing BPV involves the use of a pan-reactive antibody directed against the papillomavirus group-specific antigen (GSA), which is a conserved epitope on the major capsid protein L1. This antibody, raised against BPV, cross-reacts with PVs of many species, making it a versatile tool in veterinary diagnostics. In classic skin papillomas, IHC with an anti-BPV antibody reveals a specific, intense nuclear staining concentrated in the cells of the upper stratum spinosum and granular layer, precisely where productive viral replication and capsid assembly occur [1, 19, 20]. The staining is often absent in the basal and parabasal layers, where the virus is maintained in a latent, non-productive state. This topographical distribution of viral antigen is a powerful diagnostic correlate, confirming the active presence of the virus in the lesional tissue [1]. In a study of canine papillomatosis, the same anti-BPV antibody demonstrated immunostaining in the nuclei of keratinocytes within the finger-like projections of the lesion, confirming the diagnosis [19].

Expression of Host Cell Proteins: IHC also allows for the investigation of host cellular proteins that are dysregulated by BPV infection, providing insights into the oncogenic process. The tumor suppressor protein p53, a critical guardian of the genome, is frequently a target of viral oncoproteins. In cutaneous BPV lesions, a significant proportion of cells have been shown to exhibit strong, aberrant cytoplasmic, and perinuclear staining for p53, particularly in the basal and parabasal layers [20]. This pattern of expression, which is distinct from the normal nuclear distribution, suggests that BPV infection may sequester and inactivate p53, contributing to genomic instability and the inhibition of apoptosis. Similarly, IHC has been used to demonstrate the expression of the innate immune receptor Mincle on the cell surface of urothelial carcinoma cells in BPV-2/-13-associated bladder tumors [18]. The finding that Mincle, a C-type lectin typically associated with professional antigen-presenting cells, is expressed by cancer cells has profound implications for our understanding of the tumor microenvironment and the potential for immunotherapy, such as BCG therapy, in this model [12, 18]. These IHC data provide crucial spatial context that cannot be gleaned from bulk molecular analyses, pinpointing the exact cell types involved in the host’s response to the virus.

In summary, a robust diagnostic algorithm for BPV infection relies on the sequential and complementary use of these four core techniques. Electron microscopy provides the definitive viral morphology, molecular characterization offers unparalleled sensitivity, high-throughput genotyping, and viral load quantification, histopathology defines the tissue-level pathology and tumor grade, and immunohistochemistry localizes the antigenic footprint of the virus and its effect on the host cell within the diseased tissue. The integration of these modalities ensures a comprehensive and accurate diagnosis, which is the cornerstone of effective clinical management, epidemiological surveillance, and the development of targeted control strategies.

Genomic Diversity and Evolutionary Dynamics: Full Genome Analysis and Strain Variability

Bovine papillomaviruses (BPVs) constitute a genetically heterogeneous assemblage of circular double-stranded DNA viruses belonging to the family Papillomaviridae. Their genomes, typically 7.2–8.0 kb in size, encode five to eight early genes (E1–E8) and two late genes (L1, L2), flanked by a non-coding upstream regulatory region (URR). The past decade has witnessed an exponential increase in full-genome sequencing of BPV isolates from diverse geographic regions and host species, revealing substantial intra-type and inter-type variability that underpins their tropism, pathogenicity, and evolutionary trajectories. This section synthesizes the current state of knowledge regarding BPV genomic architecture, strain diversity, and evolutionary dynamics, drawing upon full-genome analyses from cattle, equids, and wildlife reservoirs.

Full-Genome Architecture and Conserved Functional Domains

The complete nucleotide sequences of several BPV types have been determined, providing a foundation for understanding essential genetic elements. The prototype BPV-1 genome (reference strain X02346.1) comprises 7945 base pairs, with a G+C content of approximately 53%. Within the Deltapapillomavirus genus, BPV-1 and BPV-2 share a canonical organization: early genes E1 (helicase), E2 (transcriptional regulator), E5 (major oncoprotein), and E6/E7 (transforming proteins), together with the late capsid genes L1 and L2 [2, 11]. The URR contains binding sites for cellular transcription factors and the viral E2 protein, and minor alterations here can profoundly affect replication efficiency. Shimakura et al. [2] identified a single nucleotide insertion in the URR of a BPV-1 isolate from a severely affected calf, suggesting that even subtle URR polymorphisms may modulate viral gene expression and lesion severity.

For the Xipapillomavirus genus, full-genome analysis of BPV-15 (strain Aks-02) from a dairy cow in southern Xinjiang, China, disclosed a 7189 bp genome with a G+C content of 42.5% [8]. This genome encodes five early genes (E8, E7, E1, E2, E4) and two late genes, but notably lacks E5 and E6, a hallmark of xipapillomaviruses. The E7 protein contains a consensus CX₂CX₂₉CX₂C zinc-binding domain and an LxCxE motif, confirming its capacity to bind retinoblastoma protein and drive cell cycle progression [8]. The L1 ORF of BPV-15 is 99% identical to that of the putative BAPV-3 reference, illustrating the close relatedness among xipapillomaviruses.

Intra-Type Sequence Variability and Geographic Clades

Phylogenetic analyses consistently demonstrate that BPV-1 and BPV-2 isolates cluster into geographically distinct clades. Peng et al. [11] performed full-genome sequencing of four BPV-1 strains from Chinese native cattle (JX180408, LA150909, HX160815, BS160810) and revealed that all belong to subtype A, which appears more ancestral relative to European isolates (subtype B). The European clade exhibited more recent divergence, consistent with a model of co-evolution with cattle populations following domestication and transcontinental dispersal. Within the Chinese subtype A, nucleotide identity among isolates exceeded 98.5%, while comparison with European references showed approximately 96–97% identity, underscoring the presence of distinct phylogeographic lineages [11].

Similarly, in Egypt, BPV-1 sequences from both cattle and equine sarcoids fall into two clades without host-specific differentiation [10]. Sequencing of the L1 gene (amino acids 319–454) in 14 isolates from horses and donkeys revealed three conservative mutations (D346N, Q398E, F441Y) and two non-conservative substitutions (T348N, K351T) compared to reference strain KX907623 [10]. Illumina whole-genome sequencing confirmed that isolates from donkeys and horses were identical, with 98.5% identity to the closest reference. Notably, the highest variability was observed in the E5 and L2 ORFs: substitution rates reached 5.18% (nucleotide) and 6.81% (amino acid) in E5, while E4 remained relatively conserved (0.5% nt, 0.89% aa). This hypervariability in E5, the primary oncoprotein driving mesenchymal proliferation, may reflect immune selection pressure or adaptation to different host cellular environments [10].

BPV-2 also displays marked intra-type diversity. The first molecular characterization of BPV-2 in Egypt identified a strain with 99.6% L1 nucleotide identity to reference delta papillomavirus-4 sequences from Brazil and China [1]. Partial L1 sequences (accessions MW289843 and MW289844) clustered tightly with global BPV-2 isolates, suggesting that BPV-2 has circulated broadly with minimal divergence, yet even 0.4% variation may correspond to amino acid substitutions affecting neutralization epitopes. In urothelial tumors of cattle, BPV-2 and BPV-13 E5 oncogenes are frequently co-detected, and comparative analyses of their E5 sequences reveal conserved transmembrane domains essential for activation of growth factor receptors [6, 18].

Novel Types and Putative New Species

The use of degenerate PCR primers (e.g., FAP primers targeting the L1 gene) has uncovered a previously unrecognized diversity of BPV types, particularly in cattle with multiple papillomatosis. Claus et al. [15] examined fifteen skin warts from six Brazilian cattle and identified four established BPV types (BPV-1, -2, -6, and -8) plus three novel putative types designated BPV/BR-UEL3, BPV/BR-UEL4, and BPV/BR-UEL5. Double infections were found in four animals, and a single herd harbored six distinct viral types. This level of co-infection and diversity rivals that seen in human papillomavirus infections and suggests that the BPV virome is far richer than previously appreciated [15]. The presence of multiple types within individual animals complicates epidemiological studies but highlights the need for metagenomic sequencing to capture the full spectrum of circulating viruses.

Full genome characterization of BPV-15 [8] further expands the known diversity within the Xipapillomavirus genus. The close relationship (99% L1 identity) to BAPV-3 indicates that several putative types await formal classification. Ongoing surveillance in China, Brazil, and the Middle East continues to yield isolates that do not cluster with any described type, implying that the current taxonomy underestimates true BPV diversity [4, 11, 15].

Evolutionary Dynamics and Selection Pressures

The evolutionary forces shaping BPV genomes can be inferred from dN/dS ratios across ORFs. In BPV-1, the L1 capsid gene is generally under purifying selection, as expected for a structural protein essential for host entry and immune recognition. However, Abdul-Zahra et al. [4] identified two non-synonymous mutations in the L1 gene of Iraqi BPV-1 isolates: SER31ASN and ALA55ASP. The ALA55ASP substitution lies in Loop 1 of the L1 hexameric capsid, a region critical for antibody neutralization. While the overall dN/dS ratio remains low (<0.5), these localized changes may reflect ongoing adaptation to host antibody pressure, with potential implications for vaccine design [4]. In contrast, the E5 oncogene exhibits elevated nucleotide and amino acid diversity, with dN/dS approaching or exceeding 1 in some lineages, indicative of diversifying selection [10]. This pattern is consistent with E5’s role as a transmembrane adaptor protein that interacts with host signaling complexes, where changes in amino acid sequence can alter binding partners or evade immune detection.

Phylogenetic analyses of L1 and E5 sequences from diverse geographic origins support a model of periodic introduction and local expansion. The Egyptian BPV-2 strain shows 99.6% identity to Brazilian and Chinese isolates, suggesting recent transcontinental movement, likely through international livestock trade [1]. Similarly, BPV-1 detected in Egyptian sarcoids clusters with European and Indian lineages, pointing to multiple independent introductions [7, 10]. The lack of strict host-specific clustering (cattle vs. horses) in BPV-1 indicates that cross-species transmission is frequent and that the virus can circulate without species barriers, a finding with profound implications for equine sarcoid epidemiology and control [10].

Genomic Variability and Pathogenicity

The relationship between genomic variation and clinical phenotype remains an active area of investigation. The single nucleotide insertion in the URR of a BPV-1 strain from a calf with severe generalized papillomatosis [2] provides a potential link between regulatory region changes and disease severity. Insertions in the URR can alter transcription factor binding sites, potentially upregulating E5 or E6/E7 expression and enhancing tumorigenicity. In BPV-15, the absence of E5 and E6 suggests a distinct mechanism of transformation, possibly relying solely on E7-mediated Rb inactivation [8]. Such genomic diversity likely explains the spectrum of clinical presentations, from benign cutaneous warts (xipapillomaviruses) to invasive fibropapillomas and urinary bladder neoplasia (deltapapillomaviruses) [6, 8, 15].

High-throughput sequencing technologies have begun to reveal the extent of subgenomic variability. Illumina sequencing of equine sarcoid-associated BPV-1 genomes detected variations at frequencies as low as 1–5% within a single lesion, consistent with a quasispecies model [10]. These minor variants may constitute a reservoir for rapid adaptation under selective pressure from host immunity or antiviral therapy. The presence of both wild-type and variant L1 sequences within the same tumor underscores the dynamic evolutionary landscape of BPV infections and the need for deep sequencing to capture the full mutational spectrum.

Global Perspectives and Economic Significance

The World Organisation for Animal Health (WOAH) lists bovine papillomatosis as a notifiable disease in many countries due to its economic impact on the cattle industry and its association with equine sarcoids, which reduce working ability of draft animals and cause cosmetic defects in sport horses [10, 16]. Genomic diversity directly influences diagnostic sensitivity and vaccine efficacy. For instance, the high sequence conservation of the L1 gene across BPV-1 isolates has enabled the development of broad-spectrum molecular diagnostics, including recombinase polymerase amplification assays that detect multiple genotypes [7]. However, the discovery of novel types such as BPV-15 and the putative Brazilian types necessitates periodic reassessment of primer target regions to avoid false negatives [8, 15].

The evolutionary dynamics of BPVs also intersect with oncogenesis in ways that mirror human papillomavirus biology. The World Health Organization (WHO) recognizes that understanding PV genomic diversity is essential for developing prophylactic L1 virus-like particle vaccines. In the bovine system, L1-based vaccines have shown promise in experimental settings [5], but vaccine coverage may be compromised if circulating strains diverge from the vaccine type at critical epitopes. Continued genomic surveillance, supported by open-access databases, is imperative to monitor the emergence of antigenic variants and to inform control strategies.

In summary, full-genome analyses have transformed our understanding of BPV diversity and evolution. From the URR insertion that escalates pathogenicity to the hypervariable E5 that drives tumor progression, each genomic region contributes to a complex interplay between virus and host. The identification of geographically structured clades, novel types, and cross-species transmission events underscores the need for a global, One Health approach to BPV research and management. Further metagenomic studies in under-sampled regions will undoubtedly unveil additional diversity and refine our phylogenetic frameworks, ultimately enhancing diagnostics, vaccines, and therapeutic interventions.

Transmission Pathways, Risk Factors, and Economic Impact in Cattle Populations

Transmission Pathways

The epidemiology of bovine papillomavirus (BPV) in cattle populations is governed by a complex interplay of direct, indirect, and potentially vertical transmission routes, each shaped by the virus's remarkable stability in the environment and its strict epithelial tropism. BPV is an epitheliotropic, non-enveloped DNA virus that relies on micro-abrasions or breaches in the skin or mucosal epithelium to initiate infection, as it cannot penetrate intact keratinized surfaces. This fundamental biological requirement dictates the primary transmission pathways observed in both endemic and outbreak settings.

Direct contact remains the most well-documented and significant route of horizontal transmission. The virus is shed in high concentrations from exophytic papillomas, which are fragile, cauliflower-like growths that readily desquamate viral particles into the environment. Calves that are co-mingled with infected animals, particularly in confinement operations or during transport, are at elevated risk. The presence of generalized papillomatosis, where multiple warts cover the head, neck, and trunk, creates a high viral shedding burden within the herd. Importantly, BPV transmission is not limited to overtly diseased animals; subclinical shedders may perpetuate the infection cycle, although the degree to which this occurs remains poorly quantified in field studies. The "nose-to-flank" contact common during grooming, suckling, and group housing facilitates inoculation of the virus into minor skin wounds, a mechanism corroborated by the frequent localization of warts on the head, neck, and teats [13, 14].

Indirect transmission via fomites and arthropod vectors constitutes a second major pathway, particularly in pasture-based systems. BPV is remarkably stable outside the host and can persist in the environment on contaminated bedding, feed troughs, water sources, and grooming equipment. The virus's resistance to desiccation and many common disinfectants prolongs its infectivity in barn environments, enabling transmission even in the absence of direct animal-to-animal contact. Evidence from clinical case reports in Nigeria and other tropical regions implicates arthropod vectors, including ticks and biting flies, in mechanical transmission [13, 14]. Ticks have been observed on the ears and perineum of affected animals, and their hematophagous behavior can create the requisite micro-wounds while simultaneously carrying viral particles from an infected host to a susceptible one. This vector-mediated route is particularly insidious because it circumvents standard biosecurity measures aimed at controlling direct contact.

Vertical transmission, while less well-documented than horizontal pathways, is suggested by observations of papillomatosis in very young calves, sometimes at or shortly after birth. The occurrence of severe, generalized papillomatosis in a seven-month-old calf described in one study raises the possibility of perinatal infection, though the authors emphasize that the gastrointestinal lesions in that case were fibropapillomas that lacked detectable BPV antigen, suggesting potential hematogenous dissemination of viral DNA rather than true transplacental infection [2]. Further research is needed to clarify whether true in utero transmission occurs or if neonatal infections arise from contaminated environments or colostrum. Notably, a study of the vaginal virobiota of healthy mares using ultra-sensitive digital droplet PCR demonstrated the presence of ovine papillomavirus DNA in 26.6% of samples, raising the analogous question of whether reproductive tract carriage occurs in cattle and could facilitate transmission at parturition [21]. Such a hypothesis deserves investigation, particularly given the capacity of BPV to infect urothelial cells and cause abortive infections in the bladder [6, 18].

Risk Factors for Infection and Disease Progression

The transition from BPV exposure to clinical papillomatosis is governed by a constellation of host, pathogen, and environmental risk factors. Understanding these factors is critical for designing effective control programs, as many are modifiable through management interventions.

Host-Related Risk Factors: Age and Immunocompetence

Age is one of the most consistently identified risk factors. Young cattle, particularly those between six months and two years of age, exhibit the highest incidence of cutaneous papillomatosis. This age-related susceptibility is likely multifactorial, reflecting the immaturity of the adaptive immune system, the waning of maternally derived antibodies, and the increased likelihood of skin trauma during weaning and social integration. In the severe case of a seven-month-old Holstein calf that developed numerous skin papillomas and gastrointestinal nodules, the authors noted that the animal's age aligned with the typical peak of BPV-1 susceptibility [2]. Conversely, adult cattle that have experienced prior exposure often develop immunity, although the duration of protection is variable and may be type-specific, explaining why recurrent infections with different BPV genotypes are possible.

The immune status of the host is a critical determinant of both susceptibility to infection and the likelihood of disease regression. Spontaneous regression of papillomas, observed in both cattle and other species, is mediated by cell-mediated immunity, particularly CD8+ T-cell responses directed against viral early proteins. In dogs with canine papillomavirus type 2 footpad papillomatosis, regression occurred after a biopsy, suggesting that even minor immune stimulation can tip the balance toward clearance in immunocompetent individuals [19]. In cattle, this phenomenon is exploited therapeutically by autogenous vaccines. Conversely, immunosuppression, whether induced by concurrent viral infections (e.g., bovine viral diarrhea virus, bovine leukemia virus), malnutrition, stress from transport or parturition, or genetic immunodeficiency, markedly increases the risk of persistent, generalized, and severe papillomatosis. The economic consequences of such severe cases are profound, as illustrated by astasia (inability to stand) necessitating euthanasia [2].

Pathogen-Related Risk Factors: Viral Genotype and Coinfections

BPV encompasses at least 15 genotypes classified into five genera, and genotype exerts a significant influence on lesion phenotype, tissue tropism, and clinical outcome. The Xipapillomavirus genus (e.g., BPV-3, -4, -6, -9, -10, -11, -15) produces purely epithelial papillomas that are often self-limiting, whereas the Deltapapillomavirus genus (BPV-1, -2, -13) induces fibroblastic proliferation, yielding fibropapillomas that are more locally invasive and prone to persistence. This distinction has important implications for transmission, as fibroblastic lesions may be more resistant to immune-mediated regression and may shed virus for longer periods. The detection of BPV-2 for the first time in Egypt highlights the global dispersal of these genotypes and the potential for emerging strains to alter local epidemiological patterns [1].

Multiple infections with two or more BPV genotypes within a single animal are surprisingly common and constitute a major risk factor for severe, recalcitrant disease. A seminal study of Brazilian cattle herds revealed double infections in four of six animals examined, with a single herd harboring six distinct viral types, including BPV-1, -2, -6, -8, and three putative novel types [15]. This high viral diversity suggests that cattle are regularly exposed to multiple genotypes simultaneously, and coinfection may overwhelm the host's immune capacity, delay lesion regression, and increase the total viral load shed into the environment. The phenomenon of multiple BPV infections mirrors that seen in human papillomavirus infections and underscores the need for polyvalent prophylactic vaccines.

Environmental and Management Risk Factors

Husbandry practices profoundly shape BPV transmission dynamics. Intensification of production, with high stocking densities, confinement housing, and frequent mixing of animals, creates ideal conditions for virus spread. The risk is amplified by poor biosecurity, including the sharing of contaminated equipment (e.g., ear tag applicators, dehorning tools, tattooing instruments) and the failure to isolate newly introduced or quarantined animals. In a study of Egyptian cattle, the majority of cutaneous warts were observed on the head and neck, a distribution consistent with transmission during feeding from shared troughs and during social grooming [7].

Arthropod infestation is a modifiable risk factor of particular importance in tropical and subtropical regions. The case series from Nigeria noted the concomitant presence of ticks on affected animals and recommended ivermectin as both an acaricidal and adjunctive therapeutic agent [13, 14]. While ivermectin lacks direct antiviral activity, its ability to reduce tick burden may break the vector-borne transmission cycle. Other management factors include the provision of rough surfaces or poorly maintained fencing that predisposes to skin abrasions, and the practice of surgical excision of lesions without appropriate biosecurity, which can result in iatrogenic inoculation if contaminated instruments are reused.

Economic Impact in Cattle Populations

The economic burden of bovine papillomavirus infection is substantial, multifactorial, and often underestimated due to its insidious nature. Unlike acute epidemic diseases that cause overt mortality, BPV imposes a chronic drag on productivity, and the cumulative losses across the production cycle can be devastating for individual herds and national economies alike.

Direct Costs: Morbidity, Mortality, and Reduced Productivity

The most visible direct cost is the loss of milk production in dairy cattle. Papillomas on the teats and udder interfere with milking, cause pain and discomfort, and predispose the gland to secondary bacterial mastitis. Although precise estimates of milk loss attributable to BPV are scarce, the economic impact can be gauged by the fact that mastitis itself, to which BPV lesions predispose, is the costliest disease of dairy cattle globally, with losses estimated by the Food and Agriculture Organization (FAO) at billions of dollars annually. In beef operations, severe cutaneous papillomatosis reduces weight gain, delays feedlot finishing, and can lead to condemnation of hides at slaughter. The value of a bull hide afflicted with extensive warts may be reduced by 50% or more, and the cosmetic defects render the leather unsuitable for high-quality products.

Mortality is uncommon but devastating when it occurs. The calf with severe papillomatosis that developed astasia and was euthanized exemplifies the extreme end of the economic spectrum: the loss of a replacement heifer represents not only the immediate cost of the animal but also the genetic investment and rearing costs [2]. Urinary bladder tumors associated with BPV-2 and BPV-13 infection, while less common than cutaneous papillomas, are often malignant and carry a poor prognosis, further contributing to mortality losses [6, 18].

Treatment Costs and Indirect Economic Losses

Treatment of cutaneous papillomatosis is often protracted and includes surgical excision, cryotherapy, immunotherapy (e.g., autogenous vaccines), and adjunctive therapies such as ivermectin [13, 14]. Each of these interventions carries labor, material, and professional veterinary costs. In resource-limited settings, where access to sophisticated therapies is restricted, farmers may resort to repeated surgical excision, which is traumatic, time-consuming, and associated with high recurrence rates if viral clearance is incomplete. The global burden of treatment costs is amplified by the high prevalence of infection; in one study from Iraq, 84% of tumor samples tested positive for BPV-1 DNA, indicating that infection is near-ubiquitous in affected herds [4].

Indirect Costs: Market Access, Trade, and Genetic Improvement

Although BPV is not listed as a notifiable disease by the World Organisation for Animal Health (WOAH), the presence of papillomatosis can have implications for international trade in live animals, semen, and embryos. Importing countries may impose health certification requirements that exclude animals with visible warts or require proof of herd freedom, which is technically difficult to achieve given the high prevalence of subclinical infection. Furthermore, the occurrence of BPV-associated tumors in valuable breeding stock, such as the elite bulls used in artificial insemination, represents a loss of genetic potential and a setback to national breeding programs. The economic impact is particularly acute in developing countries, where the cost of control measures may exceed the financial capacity of smallholder farmers, perpetuating a cycle of infection and lost productivity that undermines efforts to improve livestock genetics and food security [7, 16]. Comparative modeling studies for other bovine viral diseases, such as bovine viral diarrhea virus, have demonstrated that eradication programs, while costly upfront, yield net economic benefits within a decade through reduced mortality, improved reproductive performance, and enhanced market access [24]. Similar cost-benefit analyses are urgently needed for BPV control interventions.

Regional and Global Burden

The economic impact of BPV is distributed unevenly across the globe, with tropical and subtropical regions bearing the heaviest burden due to favorable environmental conditions for vector transmission, higher background prevalence of immunosuppressive diseases, and limitations in veterinary infrastructure. In Egypt, molecular surveys have documented BPV-1 as the dominant genotype responsible for cutaneous papillomatosis, and the virus has been detected in multiple provinces including New Valley, Menya, and Ismailia, suggesting widespread distribution [1, 7]. In Brazil, the identification of multiple novel BPV genotypes in a single herd hints at a viral diversity that may outpace diagnostic and control capabilities [15]. In China, where BPV-1 has been characterized in native cattle breeds, the virus is recognized as a cause of "economic losses in dairy and beef production industries," though systematic surveillance data remain sparse [11]. The International Livestock Research Institute (ILRI) and FAO have highlighted the importance of such endemic diseases in the context of sustainable livestock development, yet BPV remains a neglected pathogen that receives far less research and control investment than its economic impact warrants.

References

[1] Hassanien RT, Hamdy M, Elnomrosy SM, Hussein H, Afify A, Darwish FM, et al.. Molecular characterization and pathological identification of a novel strain of delta papillomavirus-4 (bovine papillomavirus-2) in Egypt. Veterinary World. 2021. DOI: https://doi.org/10.14202/vetworld.2021.2296-2305

[2] Shimakura H, Dong J, Zhu W, Chambers J, Uchida K, Kiriki K, et al.. Full genome analysis of bovine papillomavirus type 1 derived from a calf with severe cutaneous multiple papillomatosis. Journal of Veterinary Medical Science. 2018. DOI: https://doi.org/10.1292/jvms.18-0289

[3] Gysens L, Martens A, Haspeslagh M. Performance of fine-needle aspirate testing compared with superficial swab testing for quantification of BPV-1/-2 viral load in equine sarcoids.. Research in Veterinary Science. 2023. DOI: https://doi.org/10.1016/j.rvsc.2023.04.014

[4] Abdul-Zahra HH, Khudhair YI, Al-Hraishawi H. Molecular Analysis and Bioinformatics Assessment of Full-Length L1 Gene of Bovine Papillomavirus Type-1 as a Potential DNA Vaccine Study. Veterinary Medicine International. 2025. DOI: https://doi.org/10.1155/vmi/6785087

[5] Vráblíková A, Brezani V, Pšikal I, Fraiberk M, Šebela M, Fojtíková M, et al.. Development of modern immunization agent against bovine papillomavirus type 1 infection based on BPV1 L1 recombinant protein. Frontiers in Veterinary Science. 2023. DOI: https://doi.org/10.3389/fvets.2023.1116661

[6] Falco FD, Gentile I, Cerino P, Cutarelli A, Catoi C, Roperto S. Prohibitin 2 is Involved in Parkin-Mediated Mitophagy in Urothelial Cells of Cattle Infected with Bovine Papillomavirus. Pathogens. 2020. DOI: https://doi.org/10.3390/pathogens9080621

[7] El-Tholoth M, Mauk M, Elnaker YF, Mosad SM, Tahoun A, El-Sherif M, et al.. Molecular Characterization and Developing a Point-of-Need Molecular Test for Diagnosis of Bovine Papillomavirus (BPV) Type 1 in Cattle from Egypt. Animals. 2020. DOI: https://doi.org/10.3390/ani10101929

[8] Hu J, Zhang W, Chauhan S, Shi C, Song Y, Qiu Z, et al.. Complete genome and phylogenetic analysis of bovine papillomavirus type 15 in Southern Xinjiang dairy cow. bioRxiv. 2020. DOI: https://doi.org/10.4142/jvs.2020.21.e73

[9] Pietersen I, Zyl ARv, Rybicki E, Hitzeroth I. Novel Production of Bovine Papillomavirus Pseudovirions in Tobacco Plants. Pathogens. 2020. DOI: https://doi.org/10.3390/pathogens9120996

[10] Sobhy N, Refaai W, Kumar R, Youssef CRB, Goyal S. Molecular Characterization of Bovine Deltapapillomavirus in Equine Sarcoids in Egypt. Veterinary Medicine International. 2025. DOI: https://doi.org/10.1155/vmi/9773642

[11] Peng H, Wu C, Li J, Li C, Chen Z, Pei Z, et al.. Detection and Genomic characterization of Bovine papillomavirus isolated from Chinese Native Cattle.. Transboundary and Emerging Diseases. 2019. DOI: https://doi.org/10.1111/tbed.13285

[12] Monteiro MM, Castro ELAd, Pereira AJM, Thiesen R, Thiesen RMC, Salvarani FM. BCG Immunotherapy in Equine Sarcoid Treatment: Mechanisms, Clinical Efficacy, and Challenges in Veterinary Oncology. Viruses. 2025. DOI: https://doi.org/10.3390/v17101322

[13] . Management of cutaneous bovine papillomatosis in cross bred Holstein Friesian cow – Case Report. Journal of Sustainable Veterinary and Allied Sciences. 2024. DOI: https://doi.org/10.54328/covm.josvas.2024.177

[14] Suka I, Abdulrasheed A, Abdullahi N, Musa A, Balarabe A, Umar M, et al.. Management of cutaneous bovine papillomatosis in a year-old Sokoto Gudali heifer through marginal resection. Sokoto Journal of Veterinary Sciences. 2024. DOI: https://doi.org/10.4314/sokjvs.v22i1.9

[15] Claus MP, Lunardi M, Alfieri A, Otonel RA, Sartori D, Fungaro M, et al.. Multiple bovine papillomavirus infections associated with cutaneous papillomatosis in Brazilian cattle herds.. Brazilian Archives of Biology and Technology. 2009. DOI: https://doi.org/10.1590/S1516-89132009000700013

[16] Ibrahim AH, El-Habbaa A, El-Nahas E. Virological perspectives of the current situation of oncogenic viruses in Egypt: a review. Infectious Agents and Cancer. 2025. DOI: https://doi.org/10.1186/s13027-025-00706-7

[17] Semik-Gurgul E, Ocłoń E, Zubel-Łojek J, Pędziwiatr R, Pawlina-Tyszko K. Cannabidiol-induced cellular and matrix-associated responses in primary equine sarcoid cells. Journal of Veterinary Internal Medicine. 2026. DOI: https://doi.org/10.1093/jvimsj/aalaf015

[18] Roperto S, Russo V, Esposito I, Ceccarelli D, Paciello O, Avallone L, et al.. Mincle, an Innate Immune Receptor, Is Expressed in Urothelial Cancer Cells of Papillomavirus-Associated Urothelial Tumors of Cattle. PLoS ONE. 2015. DOI: https://doi.org/10.1371/journal.pone.0141624

[19] Iyori K, Inai K, Shimakura H, Haga T, Shimoura H, Imanishi I, et al.. Spontaneous regression of canine papillomavirus type 2-related papillomatosis on footpads in a dog. Journal of Veterinary Medical Science. 2019. DOI: https://doi.org/10.1292/jvms.19-0136

[20] Al-Dabhawi AH, Al-Salihi K, Erzuki IA, Ali TAH. Clinico-Histopathological and Immunohistochemistry Study of Ruminant’s Cutaneous Papillomavirus in Al-Muthanna Veterinary Hospital/ Iraq. . 2019. DOI: https://doi.org/10.20944/PREPRINTS201902.0177.V1

[21] Cutarelli A, Passantino G, Razzuoli E, Serpe F, Leonardi L, Zizzo N, et al.. Digital droplet PCR-based detection and quantification of ovine papillomavirus DNA from the vaginal virobiota of healthy mares. Scientific Reports. 2025. DOI: https://doi.org/10.1038/s41598-025-94279-5

[22] Johnston D, Earley B, Cormican P, Murray G, Kenny D, Waters S, et al.. Illumina MiSeq 16S amplicon sequence analysis of bovine respiratory disease associated bacteria in lung and mediastinal lymph node tissue. BMC Veterinary Research. 2017. DOI: https://doi.org/10.1186/s12917-017-1035-2

[23] Anis E, Hawkins I, Ilha M, Woldemeskel M, Saliki J, Wilkes R. Evaluation of Targeted Next-Generation Sequencing for Detection of Bovine Pathogens in Clinical Samples. Journal of Clinical Microbiology. 2018. DOI: https://doi.org/10.1128/JCM.00399-18

[24] Graham D, More S, O'Sullivan P, Lane E, Barrett D, Lozano J, et al.. The Irish Programme to Eradicate Bovine Viral Diarrhoea Virus, Organization, Challenges, and Progress. Frontiers in Veterinary Science. 2021. DOI: https://doi.org/10.3389/fvets.2021.674557