Abstract: Sarcoids are the most common cutaneous tumors in horses, representing up to 90% (35%-90%) of skin neoplasms. Mostly caused by Bovine Papillomavirus (BPVs) infections, sarcoids are highly resistant to therapy and prone to recurring, posing a significant threat to equine health. The aim of this study is to explore molecular pathogenetic mechanisms underlying the development of equine sarcoids, by applying transcriptomic approach. After testing samples for viral DNA, both mRNA and small RNA expression was analyzed via high-throughput Illumina sequencing comparing 12 sarcoids and 12 healthy skin samples as controls. Differentially expressed genes (DEGs), DE miRNAs (sarcoids vs. controls) and miRNA-DEG couples with opposite expression trends, were retrieved and subjected to a functional analysis. Over 6K DEGs emerged, 3620 down-regulated and 2415 up-regulated along with 145 DE miRNAs, 56 downregulated and 89 upregulated. Among the enriched biological processes for DEGs, some were related to growth factors production and collagen binding, cell migration and proliferation, tissue morphogenesis and inflammatory response. Interestingly, "Pathways in cancer" and "Hippo signaling pathway" were enriched KEGG pathways for the miRNA-DEG couples. Our data identified a great transcription discrepancy between sarcoid lesions and healthy skin with an overall enrichment for processes related to cellular transformation. RNA-seq sequencing depth allowed the search for candidate chimeric transcripts associated with viral integration events. Chimeric RNAs can influence gene regulation and may contribute to tumor growth and immune modulation. Via computational analysis we identified six fusion loci in tumor samples and in two sarcoid margins, with the most frequent event involving WNT10B and FKBP11. This fusion, detected in 6/10 sarcoids, is of particular interest since WNT10B activates the WNT/β-catenin cascade, while FKBP11 has been implicated in osteosarcoma progression. Although functional validation is ongoing, this represents the first report of chimeric transcripts in equine sarcoids, opening new perspectives on BPV-driven oncogenesis.
The Equine Research Bank provides access to a large database of publicly available scientific literature. Inclusion in the Research Bank does not imply endorsement of study methods or findings by Mad Barn.
This research summary has been generated with artificial intelligence and may contain errors and omissions. Refer to the original study to confirm details provided. Submit correction.
Overview
This study investigates the molecular mechanisms behind equine sarcoids, common skin tumors in horses primarily caused by Bovine Papillomavirus type 1 (BPV1), using transcriptomic profiling to identify changes in gene and microRNA expression and potential viral integration events contributing to cancer development.
Background
Equine Sarcoids: The most frequent skin tumors in horses, accounting for up to 90% of skin neoplasms.
Cause: Primarily associated with infection by Bovine Papillomaviruses, especially BPV1.
Clinical Challenge: Sarcoids are difficult to treat, resistant to therapies, and prone to recurrence, posing a major health issue in equine populations.
Objective
To explore the molecular and pathogenetic mechanisms involved in the development of equine sarcoids, focusing on gene expression changes and the presence of chimeric (fusion) transcripts due to viral integration, by applying a transcriptomic approach.
Methods
Sample Collection: 12 equine sarcoid tumor samples and 12 healthy skin samples were collected.
Viral DNA Testing: Samples were tested for BPV DNA to confirm viral presence.
Sequencing: High-throughput Illumina sequencing was performed to analyze both messenger RNA (mRNA) and small RNA (including microRNAs) expression profiles.
Analysis: Differential expression analysis compared sarcoids versus healthy controls to identify differentially expressed genes (DEGs) and microRNAs (DE miRNAs).
Functional Enrichment: Biological processes and pathways enriched among DEGs and miRNA-DEG pairs (with opposite expression) were identified, focusing on cancer-related pathways.
Chimeric Transcript Search: Deep RNA sequencing data was computationally analyzed to detect fusion transcripts potentially arising from viral integration events.
Key Results
Differentially Expressed Genes: More than 6,000 DEGs were identified — 3,620 genes were downregulated and 2,415 were upregulated in sarcoids compared to controls.
Differentially Expressed microRNAs: 145 DE miRNAs were found, with 56 downregulated and 89 upregulated in tumor samples.
Biological Processes Enriched: DEGs were enriched in processes such as:
Growth factor production
Collagen binding
Cell migration and proliferation
Tissue morphogenesis
Inflammatory response
KEGG Pathways: Notably enriched pathways included “Pathways in cancer” and the “Hippo signaling pathway,” both involved in cell growth regulation and oncogenesis.
Significance of Transcriptional Findings
The study demonstrated significant transcriptional differences in sarcoid lesions versus normal skin, highlighting molecular changes consistent with cellular transformation and tumor progression.
Using deep RNA-seq, six fusion loci were identified in tumor samples, and two in sarcoid margins, representing potential sites of viral integration.
Most Frequent Fusion: A chimeric transcript involving genes WNT10B and FKBP11 was detected in 6 out of 10 tumors analyzed.
Biological Relevance:
WNT10B is important because it activates the WNT/β-catenin signaling pathway, which is critical for regulating cell proliferation and differentiation and frequently implicated in cancer.
FKBP11 has been linked to tumor progression in osteosarcoma, suggesting a role in oncogenic processes.
This discovery is the first report of chimeric transcripts in equine sarcoids, potentially implicating new mechanisms of BPV-driven tumorigenesis.
Impact and Future Directions
The identified gene expression changes and discovery of fusion transcripts open new avenues to understand how BPV infection promotes sarcoid tumor development.
Chimeric RNAs may alter gene regulation, contribute to tumor growth, and affect immune responses within the tumor microenvironment.
Ongoing functional validation is necessary to confirm the biological roles of these fusion transcripts.
This research enhances understanding of viral oncogenesis in horses and may guide future diagnostic or therapeutic strategies targeting molecular pathways and viral integration events.
Cite This Article
APA
Mecocci S, Capomaccio S, Porcellato I, Dell'Anno F, Ratto R, Mechelli L, De Paolis L, Fruscione F, Passeri B, Gialletti R, Pepe M, Ghelardi A, Razzuoli E, Cappelli K.
(2026).
Integrated transcriptomic profiling reveals oncogenic pathways and chimeric transcripts in equine sarcoid lesions with predominant BPV1 detection.
Front Mol Biosci, 13, 1818241.
https://doi.org/10.3389/fmolb.2026.1818241
Department of Veterinary Medicine, University of Perugia, Perugia, Italy.
Sports Horse Research Center (CRCS), University of Perugia, Perugia, Italy.
Capomaccio, Stefano
Department of Veterinary Medicine, University of Perugia, Perugia, Italy.
Sports Horse Research Center (CRCS), University of Perugia, Perugia, Italy.
Porcellato, Ilaria
Department of Veterinary Medicine, University of Perugia, Perugia, Italy.
Sports Horse Research Center (CRCS), University of Perugia, Perugia, Italy.
Dell'Anno, Filippo
Department of Veterinary Medicine, University of Perugia, Perugia, Italy.
National Reference Center of Veterinary and Comparative Oncology (CEROVEC), Genoa, Italy.
Ratto, Roberta
Department of Veterinary Medicine, University of Perugia, Perugia, Italy.
Sports Horse Research Center (CRCS), University of Perugia, Perugia, Italy.
Mechelli, Luca
Department of Veterinary Medicine, University of Perugia, Perugia, Italy.
Sports Horse Research Center (CRCS), University of Perugia, Perugia, Italy.
De Paolis, Livia
National Reference Center of Veterinary and Comparative Oncology (CEROVEC), Genoa, Italy.
Fruscione, Floriana
National Reference Center of Veterinary and Comparative Oncology (CEROVEC), Genoa, Italy.
Passeri, Benedetta
Department of Veterinary Medicine, University of Parma, Parma, Italy.
Gialletti, Rodolfo
Department of Veterinary Medicine, University of Parma, Parma, Italy.
Pepe, Marco
Department of Veterinary Medicine, University of Perugia, Perugia, Italy.
Sports Horse Research Center (CRCS), University of Perugia, Perugia, Italy.
Ghelardi, Alessandro
UOC Ostetricia e Ginecologia, Azienda Usl Toscana Nord-Ovest, Massa, Italy.
Razzuoli, Elisabetta
National Reference Center of Veterinary and Comparative Oncology (CEROVEC), Genoa, Italy.
Cappelli, Katia
Department of Veterinary Medicine, University of Perugia, Perugia, Italy.
Sports Horse Research Center (CRCS), University of Perugia, Perugia, Italy.
Conflict of Interest Statement
The author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
References
This article includes 64 references
Alshahrani S. H., Rakhimov N., Gupta J., Hassan Z. F., Alsalamy A., Saleh E. A. M.. The mechanisms, functions and clinical applications of miR-542–3p in human cancers. Pathol. - Res. Pract. 248, 154724.
Bortnik V., Wu M., Julcher B., Salinas A., Nikolic I., Simpson K. J.. Loss of HPV type 16 E7 restores cGAS-STING responses in human papilloma virus-positive oropharyngeal squamous cell carcinomas cells. J. Microbiol. Immunol. Infect. 54, 733–739.
Cassaro A., Grillo G., Notaro M., Gliozzo J., Esposito I., Reda G.. FZD6 triggers Wnt–signalling driven by WNT10B expression and highlights new targets in T‐cell acute lymphoblastic leukemia. Hematol. Oncol. 39, 364–379.
Cheng D. D., Yu T., Hu T., Yao M., Fan C. Y., Yang Q. C.. MiR-542-5p is a negative prognostic factor and promotes osteosarcoma tumorigenesis by targeting HUWE1. Oncotarget 6 (40), 42761–42772.
De Falco F., Cutarelli A., Pellicanò R., Brandt S., Roperto S.. Molecular detection and quantification of ovine papillomavirus DNA in equine sarcoid. Transbound. Emerg. Dis. 2024, 6453158.
Dey M., Skipar P., Bartnik E., Piątkowski J., Sulejczak D., Czarnecka A. M.. MicroRNA signatures in osteosarcoma: diagnostic insights and therapeutic prospects. Mol. Cell. Biochem. 480, 2065–2075.
Doncheva N. T., Morris J. H., Gorodkin J., Jensen L. J.. Cytoscape StringApp: network analysis and visualization of proteomics data. J. Proteome Res. 18, 623–632.
Gysens L., Vanmechelen B., Haspeslagh M., Maes P., Martens A.. New approach for genomic characterisation of equine sarcoid-derived BPV-1/-2 using nanopore-based sequencing. Virol. J. 19, 8.
Hainisch E. K., Brandt S., Shafti-Keramat S., Van Den Hoven R., Kirnbauer R.. Safety and immunogenicity of BPV-1 L1 virus-like particles in a dose-escalation vaccination trial in horses: VLP-based sarcoid prophylaxis. Equine Vet. J. 44, 107–111.
Hsing E.-W., Shiah S.-G., Hsiao J.-R., Lyu P.-C., Chang J.-Y.. The role of miR-450a in oral squamous cell carcinoma carcinogenesis.. FASEB J. 29, 711–722.
Hu X., Schwarz J. K., Lewis J. S., Huettner P. C., Rader J. S., Deasy J. O.. A MicroRNA expression signature for cervical cancer prognosis.. Cancer Res. 70, 1441–1448.
Jayachandran K., Pillai R., Inkman M., Ruiz F., Webster J., Schwarz J.. 13. HPV forms chimeric virus-human transcripts that affect host gene expression in cervical tumors.. Cancer Genet. 278–279.
Kamstock D. A., Ehrhart E. J., Getzy D. M., Bacon N. J., Rassnick K. M., Moroff S. D.. Recommended guidelines for submission, trimming, margin evaluation, and reporting of tumor biopsy specimens in veterinary surgical pathology.. Vet. Pathol. 48 (1), 19–31.
Kong Q., Han J., Deng H., Wu F., Guo S., Ye Z.. miR-431-5p alters the epithelial-to-mesenchymal transition markers by targeting UROC28 in hepatoma cells.. OncoTargets Ther. 11, 6489–6503.
Liao Y., Smyth G. K., Shi W.. featureCounts: an efficient general purpose program for assigning sequence reads to genomic features.. Bioinformatics 30, 923–930.
Liu M., Han Z., Zhi Y., Ruan Y., Cao G., Wang G.. Long-read sequencing reveals oncogenic mechanism of HPV-human fusion transcripts in cervical cancer.. Transl. Res. 253, 80–94.
Mosseri S, Hetzel U, Hahn S, Michaloupoulou E, Sallabank H C, Knottenbelt D C. Equine sarcoid: demonstration of matrix metalloproteinase expression. Vet. J. 202, 279–285.
Munday J S, Thomson N, Dunowska M, Knight C G, Laurie R E, Hills S. Genomic characterisation of the feline sarcoid-associated papillomavirus and proposed classification as Bos taurus papillomavirus type 14. Vet. Microbiol. 177, 289–295.
Nguyen N D, Deshpande V, Luebeck J, Mischel P S, Bafna V. ViFi: accurate detection of viral integration and mRNA fusion reveals indiscriminate and unregulated transcription in proximal genomic regions in cervical cancer. Nucleic Acids Res. 46, 3309–3325.
Patterson M R, Cogan J A, Cassidy R, Theobald D A, Wang M, Scarth J A. The hippo pathway transcription factors YAP and TAZ play HPV-type dependent roles in cervical cancer. Nat. Commun. 15, 5809.
Pan Y, Zhao Y, Lihui L, Xie Y, Zou Q. MiR-337-3p suppresses migration and invasion of breast cancer cells by downregulating ESRP1. Acta Histochem. 123 (7), 151777.
Peter Y Y, Balkhi M Y, Ladner K J, Alder H, Yu L, Mo X. A selective screening platform reveals unique global expression patterns of microRNAs in a cohort of human soft-tissue sarcomas. Lab Invest. 96 (4), 481–491.
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 13, 2498–2504.
Shinozaki A, Sakatani T, Ushiku T, Hino R, Isogai M, Ishikawa S. Downregulation of MicroRNA-200 in EBV-associated gastric carcinoma. Cancer Res. 70, 4719–4727.
Smith T, Heger A, Sudbery I. UMI-tools: modeling sequencing errors in unique molecular identifiers to improve quantification accuracy. Genome Res. 27, 491–499.
Uhrig S, Ellermann J, Walther T, Burkhardt P, Fröhlich M, Hutter B. Accurate and efficient detection of gene fusions from RNA sequencing data. Genome Res. 31 (3), 448–460.
Wu Y, Wang W, Yang A-G, Zhang R. The microRNA-424/503 cluster: a master regulator of tumorigenesis and tumor progression with paradoxical roles in cancer. Cancer Lett. 494, 58–72.
Zeng D, Li J, Yuan X, Cai F, Yu B, Liu L. FKBP11 improves the malignant property of osteosarcoma cells and acts as a prognostic factor of osteosarcoma. Aging 15, 2450–2459.
Zhang Z, Zhang L, Wang B, Wei R, Wang Y, Wan J. MiR-337–3p suppresses proliferation of epithelial ovarian cancer by targeting PIK3CA and PIK3CB. Cancer Lett 469, 54–67.