Transcriptome analysis of equine sarcoids

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.

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 Induced Expression of BPV E4 Gene in Equine Adult Dermal Fibroblast Cells as a Potential Model of Skin Sarcoid-like Neoplasia

The equine sarcoid is one of the most common neoplasias in the Equidae family. Despite the association of this tumor with the presence of bovine papillomavirus (BPV), the molecular mechanism of this lesion has not been fully understood. The transgenization of equine adult cutaneous fibroblast cells (ACFCs) was accomplished by nucleofection, followed by detection of molecular modifications using high-throughput NGS transcriptome sequencing. The results of the present study confirm that BPV-E4- and BPV-E1^E4-mediated nucleofection strategy significantly affected the transcriptomic alterations, leading to sarcoid-like neoplastic transformation of equine ACFCs. Furthermore, the results of the current investigation might contribute to the creation of in vitro biomedical models suitable for estimating the fates of molecular dedifferentiability and the epigenomic reprogrammability of BPV-E4 and BPV-E4^E1 transgenic equine ACFC-derived sarcoid-like cell nuclei in equine somatic cell-cloned embryos. Additionally, these in vitro models seem to be reliable for thoroughly recognizing molecular mechanisms that underlie not only oncogenic alterations in transcriptomic signatures, but also the etiopathogenesis of epidermal and dermal sarcoid-dependent neoplastic transformations in horses and other equids. For those reasons, the aforementioned transgenic models might be useful for devising clinical treatments in horses afflicted with sarcoid-related neoplasia of cutaneous and subcutaneous tissues.

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.

Mobility and Invasion Related Gene Expression Patterns in Equine Sarcoid

Simple Summary

The current studies profiled the expression of five equine sarcoid cell genes related to cell mobility and invasion (cell cycle control binding protein alpha, coronin 1b, metalloproteinase 2, tissue inhibitor of metalloproteinases 3 and vimentin) and compared the expression of these genes in healthy skin fibroblasts. Cells were collected from healthy and sarcoid-affected skin biopsies obtained by a qualified veterinarian. Gene expression patterns were investigated under two different conditions of cell culture, with high and low availability of nutritional components in the culture medium. The results showed significant differences in the expression of the two analyzed genes (coronin 1b and vimentin) depending on culture conditions. The obtained results emphasize the complexity of the genomic background of sarcoids and indicate the importance of further research on genes related to the physiological changes that occur in sarcoids.

Abstract

Sarcoids are the most common skin neoplasm in the Equidae family. Sarcoids are benign, but may cause severe damage in affected animals. Due to the high risk of post-treatment recurrence and the lack of an effective method of treatment, it is reasonable to perform studies on the molecular aspects of this neoplasm. Therefore, the present studies analyzed five genes (cell cycle control binding protein alpha, coronin 1b, metalloproteinase 2, tissue inhibitor of metalloproteinases 3 and vimentin) related to cell mobility and invasion traits. Primary healthy fibroblasts and sarcoid cells were obtained from skin biopsies. Cell lines were cultured in two different medium types with different concentrations of foetal bovine serum (10% and 0.5% FBS) to study its influence on the analyzed genes. Gene expression was measured using the real-time PCR method. The results showed significant differences in two genes (coronin and vimentin) depending on culture conditions. In conclusion, the results enabled finding two new genes, related to cell motility and invasion traits, in which gene expression is deregulated. Results of the study may put new knowledge into the complexity of the genetic background of this disease and show the importance of further analysis on this subject.

Bone formation transcripts dominate the differential gene expression profile in an equine osteoporotic condition associated with pulmonary silicosis

Osteoporosis has been associated with pulmonary silicosis in California horses exposed to soils rich in cytotoxic silica dioxide crystals, a syndrome termed silicate associated osteoporosis (SAO). The causal mechanism for the development of osteoporosis is unknown. Osteoporotic lesions are primarily located in bone marrow-rich sites such as ribs, scapula and pelvis. Gene transcription patterns within bone marrow and pulmonary lymph nodes of affected horses may offer clues to disease pathobiology. Bone marrow core and tracheobronchial lymph node tissue samples harvested postmortem from affected and unaffected horses were examined histologically and subjected to RNA sequencing (RNA-seq). Sequenced data were analyzed for differential gene expression and gene ontology. Metatranscriptomic and metagenomic assays evaluated samples for infectious agents. Thirteen of 17 differentially expressed transcripts in bone marrow were linked to bone and cartilage formation such as integrin binding bone sialoprotein (log2FC = 3.39, PFDR = 0.013) and chondroadherin (log2FC = 4.48, PFDR = 0.031). Equus caballus solute carrier family 9, subfamily A2 (log2FC = 3.77, PFDR = 0.0034) was one of the four differentially expressed transcripts linked to osteoclast activity. Osteoblasts were hyperplastic and hypertrophic in bone marrow from affected horses. Biological pathways associated with skeletal morphogenesis were significantly enriched in affected horses. The 30 differentially expressed genes in affected lymph nodes were associated with inflammatory responses. Evidence of infectious agents was not found. The SAO affected bone marrow molecular signature demonstrated increased transcription and heightened activation of osteoblasts. Increased osteoblastic activity could be part of the pathological mechanism for osteoporosis or a compensatory response to the accelerated osteolysis. Transcriptome data offer gene targets for inquiries into the role of osteocytes and osteoblasts in SAO pathogenesis. Viral or bacterial infectious etiology in SAO is less likely based on metatranscriptomic and metagenomic data but cannot be completely ruled out.