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.

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.

Tracking the Molecular Scenarios for Tumorigenic Remodeling of Extracellular Matrix Based on Gene Expression Profiling in Equine Skin Neoplasia Models.

An important component of tissues is the extracellular matrix (ECM), which not only forms a tissue scaffold, but also provides the environment for numerous biochemical reactions. Its composition is strictly regulated, and any irregularities can result in the development of many diseases, including cancer. Sarcoid is the most common skin cancer in equids. Its formation results from the presence of the genetic material of the bovine papillomavirus (BPV). In addition, it is assumed that sarcoid-dependent oncogenic transformation arises from a disturbed wound healing process, which may be due to the incorrect functioning of the ECM. Moreover, sarcoid is characterized by a failure to metastasize. Therefore, in this study we decided to investigate the differences in the expression profiles of genes related not only to ECM remodeling, but also to the cell adhesion pathway, in order to estimate the influence of disturbances within the ECM on the sarcoid formation process. Furthermore, we conducted comparative research not only between equine sarcoid tissue bioptates and healthy skin-derived explants, but also between dermal fibroblast cell lines transfected and non-transfected with a construct encoding the E4 protein of the BP virus, in order to determine its effect on ECM disorders. The obtained results strongly support the hypothesis that ECM-related genes are correlated with sarcoid formation. The deregulated expression of selected genes was shown in both equine sarcoid tissue bioptates and adult cutaneous fibroblast cell (ACFC) lines neoplastically transformed by nucleofection with gene constructs encoding BPV1-E1^E4 protein. The identified genes (CD99, ITGB1, JAM3 and CADM1) were up- or down-regulated, which pinpointed the phenotypic differences from the backgrounds noticed for adequate expression profiles in other cancerous or noncancerous tumors as reported in the available literature data. Unravelling the molecular pathways of ECM remodeling and cell adhesion in the in vivo and ex vivo models of epidermal/dermal sarcoid-related cancerogenesis might provide powerful tools for further investigations of genetic and epigenetic biomarkers for both silencing and re-initiating the processes of sarcoid-dependent neoplasia. Recognizing those biomarkers might insightfully explain the relatively high capacity of sarcoid-descended cancerous cell derivatives to epigenomically reprogram their nonmalignant neoplastic status in domestic horse cloned embryos produced by somatic cell nuclear transfer (SCNT).

Assessment of BPV-1 Mediated Matrix Metalloproteinase Genes Deregulation in the In Vivo and In Vitro Models Designed to Explore Molecular Nature of Equine Sarcoids.

Matrix metalloproteinases (MMPs) represent a family of enzymes capable of biocatalytically breaking down the structural and functional proteins responsible for extracellular matrix (ECM) integrity. This capability is widely used in physiological processes; however, imbalanced MMP activity can trigger the onset and progression of various pathological changes, including the neoplasmic transformation of different cell types. We sought to uncover molecular mechanisms underlying alterations in transcriptional profiles of genes coding for MMPs, which were comprehensively identified in equine adult dermal tissue bioptates, sarcoid-derived explants, and ex vivo expanded adult cutaneous fibroblast cell (ACFC) lines subjected to inducible oncogenic transformation into sarcoid-like cells. The results strongly support the hypothesis that the transcriptional activity of MMP genes correlates with molecular modifications arising in equine dermal cells during their conversion into sarcoid cells. The alterations in MMP transcription signatures occurs in both sarcoid tissues and experimentally transformed equine ACFC lines expressing BPV1-E4^E1 transgene, which were characterized by gene up- and down-regulation patterns.

The G-Protein-Coupled Membrane Estrogen Receptor Is Present in Horse Cryptorchid Testes and Mediates Downstream Pathways.

Cryptorchidism in horses is a commonly occurring malformation. The molecular basis of this pathology is not fully known. In addition, the origins of high intratesticular estrogen levels in horses remain obscure. In order to investigate the role of the G-protein-coupled membrane estrogen receptor (GPER) and establish histological and biochemical cryptorchid testis status, healthy and cryptorchid horse testes were subjected to scanning electron microscopy analysis, histochemical staining for total protein (with naphthol blue black; NBB), acid content (with toluidine blue O; TBO), and polysaccharide content (with periodic acid-Schiff; PAS). The expression of GPER was analyzed by immunohistochemistry and Western blot. GPER-mediated intracellular cAMP and calcium (Ca2+) signaling were measured immunoenzymatically or colorimetrically. Our data revealed changes in the distribution of polysaccharide content but not the protein and acid content in the cryptorchid testis. Polysaccharides seemed to be partially translocated from the interstitial compartment to the seminiferous tubule compartment. Moreover, the markedly decreased expression of GPER and GPER downstream molecules, cAMP and Ca2+, suggests their potential role in testis pathology. Increased estrogen levels in cryptorchid conditions may be linked to disturbed GPER signaling. We postulate that GPER is a prominent key player in testis development and function and may be used as a new biomarker of horse testis in health and disease.

Comprehensive characteristics of microRNA expression profile of equine sarcoids.

Equine sarcoids are the most common neoplasms occurring in horses. Despite frequent occurrence, they are still not well described at the molecular level. Thus, in the present study, we performed a comprehensive comparative analysis of sarcoid miRNAome profile to identify aberrantly expressed microRNAs, along with their structural variants, potentially useful as biomarkers and, in a wider perspective, broaden the knowledge about this tumor and underlying mechanisms. To this end, we conducted next generation sequencing and as a result we identified both known and potentially novel miRNAs. Differential expression analysis revealed the existence of almost one hundred miRNAs being over- or underexpressed in sarcoids in comparison to healthy tissue (p-adj<0.05), of which many are known for their involvement in processes crucial for neoplastic transformation. Among upregulated miRNAs there were those associated with decreased cell adhesion abilities as well as engaged in global protein production, while downregulation of some miRs i.a. increased cell expansion abilities. Moreover, we identified altered expression levels of miRNA variants (isomiRs) between the investigated tissues. Further analysis revealed that 5' isomiRs comprise different seed sequences leading to target gene switching followed by activation of different biological pathways. Our results are the first which revealed the complexity of microRNA profiles in equine sarcoids and skin tissue, along with the dynamism of their growing in importance concomitants, namely isomiRs. They also showed miRNA molecules and biological pathways important from the sarcoid oncogenesis point of view.