Comparative analysis of DNA methylation patterns of equine sarcoid and healthy skin samples

Methylome and transcriptome data integration reveals aberrantly regulated genes in equine sarcoids

DNA methylation is a key mechanism in transcription regulation, and aberrant methylation is a common and important mechanism in tumor initiation, maintenance, and progression. To find genes that are aberrantly regulated by altered methylation in horse sarcoids, we used reduced representation bisulfite sequencing (RRBS) accompanied by RNA sequencing (RNA-Seq) for methylome (whole genome DNA methylation sequencing) and transcriptome profiling, respectively. We found that the DNA methylation level was generally lower in lesion samples than in controls. In the analyzed samples, a total of 14,692 differentially methylated sites (DMSs) in the context of CpG (where cytosine and guanine are separated by a phosphate), and 11,712 differentially expressed genes (DEGs) were identified. The integration of the methylome and transcriptome data suggests that aberrant DNA methylation may be involved in the deregulation of expression of the 493 genes in equine sarcoid. Furthermore, enrichment analysis of the genes demonstrated the activation of multiple molecular pathways related to extracellular matrix (ECM), oxidative phosphorylation (OXPHOS), immune response, and disease processes that can be related to tumor progression. The results provide further insight into the epigenetic alterations in equine sarcoids and provide a valuable resource for follow-up studies to identify biomarkers for predicting susceptibility to this common condition in horses.

Transcriptome and methylome sequencing reveals altered long non-coding RNA genes expression and their aberrant DNA methylation in equine sarcoids

Recent publications confirmed that long non-coding RNAs (lncRNAs) perform an essential function in gene-specific transcription regulation. Nevertheless, despite its important role, lncRNA has not yet been described in equine sarcoids, the skin neoplasia of horses. Therefore, the aim of this study is to deepen the knowledge about lncRNA expression in the pathogenesis of equine sarcoids and provide new insight into the regulatory function of lncRNA in the bovine papillomavirus-dependent neoplasia of horse dermal tissues. RNA sequencing (RNA-seq) data from 12 equine sarcoid samples and the corresponding controls were reanalyzed in this study. A total of 3396 differentially expressed (DE) lncRNAs and 128 DElncRNA-DE genes (DEGs) pairs were identified. Differentially expressed lncRNAs predicted target genes were enriched in pathways associated with inter alia the extracellular matrix disassembly and cancer pathways. Furthermore, methylation data from the same samples were integrated into the analysis, and 12 DElncRNAs were described as potentially disturbed by aberrant methylation. In conclusion, this study presents novel data about lncRNA's role in the pathogenesis of equine sarcoids.

Methylation Status of Gene Bodies of Selected microRNA Genes Associated with Neoplastic Transformation in Equine Sarcoids

Horses are of great importance in recreation, livestock production, as working animals in poorly developed countries, and for equine-assisted therapy. Equine sarcoids belong to the most commonly diagnosed tumors in this species. They may cause discomfort, pain, and can lead to the permanent impairment of motor function. The molecular bases of their formation are still under investigation. Our previous studies revealed altered microRNA (miRNA) expression and DNA methylation levels in sarcoid tumors. Abnormal patterns of methylation may be responsible for changes in gene expression levels, including microRNAs. Recently, the DNA methylation of gene bodies has also been shown to have an impact on gene expression. Thus, the aim of the study was to investigate the methylation pattern of gene bodies of chosen miRNAs identified in sarcoid tissue (miR-101, miR-10b, miR-200a, and miR-338-3p), which have also been established to play roles in neoplastic transformation. To this end, we applied qRT-PCR, Bisulfite Sequencing PCR (BSP), and Mquant methods. As a result, we identified the statistically significant downregulation of pri-mir-101-1, pri-mir-10b, and pri-mir-200a in the sarcoid samples in comparison to the control. The DNA methylation analysis revealed their hypermethylation. This suggests that DNA methylation may be one mechanism responsible for the downregulation of theses miRNAs. However, the identified differences in the methylation levels are not very high, which implies that other mechanisms may also underlie the downregulation of the expression of these miRNAs in equine sarcoids. For the first time, the results obtained shed light on microRNA expression regulation by gene body methylation in equine sarcoids and provide bases for further deeper studies on other mechanisms influencing the miRNA repertoire.

Equine Sarcoids-Causes, Molecular Changes, and Clinicopathologic Features: A Review

Equine sarcoid is the most common skin tumor of horses. Clinically, it occurs as a locally invasive, fibroblastic, wart-like lesion of equine skin, which has 6 clinical classes: occult, verrucose, nodular, fibroblastic, mixed, and malignant. Sarcoids may be single but multiple lesions are more frequent. The typical histological feature is increased density of dermal fibroblasts which form interlacing bundles and whorls within the dermis. Lesions are mostly persistent, resist therapy, and tend to recur following treatment. In general, sarcoids are not fatal but their location, size, and progression to the more aggressive form may lead to the withdrawal of a horse from use and serious infringement of their welfare leading to the loss of valuable animals. Bovine papillomavirus (BPV) type 1 and less commonly type 2 contribute to the development of equine sarcoid. The viral genome and proteins are detected in a high percentage of cases. Furthermore, viral oncoprotein activity leads to changes in the fibroblastic tissue similar to changes seen in other types of tumors. Equine sarcoids are characterized by a loss of tumor suppressor activity and changes allowing abnormal formation of the affected tissue, as well as y immune defense abnormalities that weaken the host's immune response. This impaired immune response to BPV infection appears to be crucial for the development of lesions that do not spontaneously regress, as occurs in BPV-infected cows.

Molecular approaches to equine sarcoids

Sarcoids are the most commonly diagnosed skin tumours in equines. Bovine papillomaviruses (BPVs) are the primary causative agent of sarcoids. There has been intensive research to discover the molecular mechanisms that may contribute to the aetiopathogenesis of this disease and tumour suppressors and proto-oncogenes known to play a role in human neoplastic conditions have been investigated in equine sarcoids. Current approaches include the identification of gene expression profiles, characterising sarcoid and normal skin tissues, and an assessment of epigenetic alterations such as microRNA differential expression and DNA methylation status. This review focuses on selected groups of genes that contribute to the molecular mechanisms of sarcoid formation. These genes have the potential to complement current clinical examinations of equine sarcoid disease in diagnosis, prognosis, therapeutic response and screening.

Differences in miRNA differential expression in whole blood between horses with sarcoid regression and progression

Background

Currently no methods are available to predict the clinical outcome of individual horses with equine sarcoid (ES) disease.

Objective

To investigate if whole blood microRNA (miRNA) profiles can predict the long‐term development of ES tumors.

Animals

Five horses with regression and 5 with progression of ES lesions monitored over 5‐7 years and 5 control horses free of ES for at least 5 years.

Methods

For this cohort study, RNA extracted from whole blood samples from the regression, progression, and control groups was used for high throughput sequencing. Known and novel miRNAs were identified using miRDeep2 and differential expression analysis was carried out by the DESeq2 algorithm. Target gene and pathway prediction as well as enrichment and network analyses were conducted using TarBase, mirPath, and metaCore from GeneGo.

Results

Fourteen miRNAs were differentially expressed between regression and progression groups after accounting for the control condition: 4 miRNAs (28.6%) were upregulated and 10 miRNAs (71.4%) were downregulated with >2‐fold change. Seven of the 10 downregulated miRNAs are encoded in an miRNA cluster on equine chromosome 24, homologous to the well‐known 14q32 cluster in humans. Their target genes show enrichment for pathways involved in viral carcinogenesis.

Conclusions and Clinical Importance

Whole blood miRNA expression profiles are associated with long‐term ES growth in horses and warrant further validation as prognostic biomarkers in a larger study cohort. Deregulation of miRNAs on equine chromosome 24 might represent a trigger for ES development.