In search of epigenetic hallmarks of different tissues: an integrative omics study of horse liver, lung, and heart.

DNA methylation and microRNA (miRNA) expression are epigenetic mechanisms essential for regulating tissue-specific gene expression and metabolic processes. However, high-resolution transcriptome, methylome, or miRNAome data is only available for a few model organisms and selected tissues. Up to date, only a few studies have reported on gene expression, DNA methylation, or miRNA expression in adult equine tissues at the genome-wide level. In the present study, we used RNA-Seq, miRNA-seq, and reduced representation bisulfite sequencing (RRBS) data from the heart, lung, and liver tissues of healthy cold-blooded horses to identify differentially expressed genes (DEGs), differentially expressed miRNA (DE miRNA) and differentially methylated sites (DMSs) between three types of horse tissues. Additionally, based on integrative omics analysis, we described the observed interactions of epigenetic mechanisms with tissue-specific gene expression alterations. The obtained data allowed identification from 4067 to 6143 DMSs, 9733 to 11,263 mRNAs, and 155 to 185 microRNAs, differentially expressed between various tissues. We pointed out specific genes whose expression level displayed a negative correlation with the level of CpG methylation and miRNA expression and revealed biological processes that they enrich. Furthermore, we confirmed and validated the accuracy of the Next-Generation Sequencing (NGS) results with bisulfite sequencing PCR (BSP) and quantitative PCR (qPCR). This comprehensive analysis forms a strong foundation for exploring the epigenetic mechanisms involved in tissue differentiation, especially the growth and development of the equine heart, lungs, and liver.

Polymorphism within IGFBP Genes Affects the Acidity, Colour, and Shear Force of Rabbit Meat.

Rabbits are important livestock animals, popular for their nutritional value. Nowadays, the molecular background of traits influencing the quality of meat and meat products is in high demand. Therefore, in the current study, we analyse the sequences of IGFBP1, IGFBP2, IGFBP4, IGFBP5, and IGFBP6 for possible polymorphisms. Based on a bioinformatics analysis in an association study on 466 animals of different breeds (New Zealand White × Flemish Giant crossbreed (9NZWxFG), Termond White (TW), Popielno White (PW), and Flemish Giant (FG)), we analyse the influence of five polymorphisms within the IGFBP genes. Statistically significant differences were found among the carcass and meat quality traits but not for all of the analysed rabbit breeds. The most promising polymorphism was g.158093018A>T within the IGFBP5 gene. The values of pH24 of m.longissimus lumborum (m.l.l.) and biceps femoris muscles (m.b.f.) were higher for the AT genotypes compared to the AA genotypes for the TW and NZWxFG crossbreeds. Also, for pH24, we found differences in ing.41594308T>C for NZWxFG, where the TT genotype values were higher than the TC values. We found differences in L*24 on m.l.l. for g.41592248A>C for NZWxFG. For m.b.f., significant differences were found in b*45 for g.3431insAC in the FG population and a*45 for g.41592248A>C and g.158093018A>T in the TW population. The shear force statistically differed for g.158093018A>T in TW rabbits and g.41592248A>C for NZWxFG. We conclude that this polymorphism may be promising for better quality rabbit meat and may be implemented in selection processes.

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.

Application of Nanopore Sequencing for High Throughput Genotyping in Horses.

Nanopore sequencing is a third-generation biopolymer sequencing technique that relies on monitoring the changes in an electrical current that occur as nucleic acids are passed through a protein nanopore. Increasing quality of reads generated by nanopore sequencing systems encourages their application in genome-wide polymorphism detection and genotyping. In this study, we employed nanopore sequencing to identify genome-wide polymorphisms in the horse genome. To reduce the size and complexity of genome fragments for sequencing in a simple and cost-efficient manner, we amplified random DNA fragments using a modified DOP-PCR and sequenced the resulting products using the MinION system. After initial filtering, this generated 28,426 polymorphisms, which were validated at a 3% error rate. Upon further filtering for polymorphism and reproducibility, we identified 9495 SNPs that reflected the horse population structure. To conclude, the use of nanopore sequencing, in conjunction with a genome enrichment step, is a promising tool that can be practical in a variety of applications, including genotyping, population genomics, association studies, linkage mapping, and potentially genomic selection.

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.

Changes in DNA Methylation and mRNA Expression in Lung Tissue after Long-Term Supplementation with an Increased Dose of Cholecalciferol.

Maintaining an appropriate concentration of vitamin D is essential for the proper functioning of the body, regardless of age. Nowadays, there are more and more indications that vitamin D supplementation at higher than standard doses may show protective and therapeutic effects. Our study identified differences in the body's response to long-term supplementation with cholecalciferol at an increased dose. Two groups of pigs were used in the experiment. The first group received a standard dose of cholecalciferol (grower, 2000 IU/kg feed, and finisher, 1500 IU/kg feed), and the second group received an increased dose (grower, 3000 IU/kg feed, and finisher, 2500 IU/kg feed). After slaughter, lung samples were collected and used for RRBS and mRNA sequencing. Analysis of the methylation results showed that 2349 CpG sites had significantly altered methylation patterns and 1116 (47.51%) identified DMSs (Differentially Methylated Sites) were related to genes and their regulatory sites. The mRNA sequencing results showed a significant change in the expression of 195 genes. The integrated analysis identified eleven genes with DNA methylation and mRNA expression differences between the analyzed groups. The results of this study suggested that an increased vitamin D intake may be helpful for the prevention of lung cancer and pulmonary fibrosis. These actions may stem from the influence of vitamin D on the expression of genes associated with collagen production, such as SHMT1, UGT1A6, and ITIH2.The anti-cancer properties of vitamin D are also supported by changes in KLHL3 and TTPA gene expression.

Another lesson from unmapped reads: in-depth analysis of RNA-Seq reads from various horse tissues.

In recent years, a vast amount of sequencing data has been generated and large improvements have been made to reference genome sequences. Despite these advances, significant portions of reads still do not map to reference genomes and these reads have been considered as junk or artificial sequences. Recent studies have shown that these reads can be useful, e.g., for refining reference genomes or detecting contaminating microorganisms present in the analyzed biological samples. A special case of this is RNA sequencing (RNA-Seq) reads that come from tissue transcriptomes. Unmapped reads from RNA-Seq have received much less attention than those from whole-genome sequencing. In particular, in the horse, an analysis of unmapped RNA reads has not been performed yet. Thus, in this study, we analyzed the unmapped reads originating from the RNA-Seq performed through the Functional Annotation of Animal Genomes (FAANG) project in the horse, using eight different tissues from two mares. We demonstrated that unmapped reads from RNA-Seq could be easily assembled into transcripts relating to many important genes present in the sequences of other mammals. Large portions of these transcripts did not have coding potential and, thus, can be considered as non-coding RNA. Moreover, reads that were not mapped to the reference genome but aligned to the entries in NCBI database of horse proteins were enriched for biological processes that largely correspond to the functions of organ from which RNA was isolated and thus are presumably true transcripts of genes associated with cell metabolism in those tissues. In addition, a portion of reads aligned to the common pathogenic or neutral microbiota, of which the most common was Brucella spp. These data suggest that unmapped reads can be an important target for in-depth analysis that may substantially enrich results of initial RNA-Seq experiments for various tissues and organs.

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 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.

Application of the targeted sequencing approach reveals the single nucleotide polymorphism (SNP) repertoire in microRNA genes in the pig genome.

MicroRNAs (miRNAs) are recognized as gene expression regulators, indirectly orchestrating a plethora of biological processes. Single nucleotide polymorphism (SNP), one of the most common genetic variations in the genome, is established to affect miRNA functioning and influence complex traits and diseases. SNPs in miRNAs have also been associated with important production traits in livestock. Thus, the aim of our study was to reveal the SNP variability of miRNA genes in the genome of the pig, which is a significant farm animal and large-mammal human model. To this end, we applied the targeted sequencing approach, enabling deep sequencing of specified genomic regions. As a result, 73 SNPs localized in 50 distinct pre-miRNAs were identified. In silico analysis revealed that many of the identified SNPs influenced the structure and energy of the hairpin precursors. Moreover, SNPs localized in the seed regions were shown to alter targeted genes and, as a result, enrich different biological pathways. The obtained results corroborate a significant impact of SNPs on the miRNA processing and broaden the state of knowledge in the field of animal genomics. We also report the targeted sequencing approach to be a promising alternative for the whole genome sequencing in miRNA genes focused studies.

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.

Detection of runs of homozygosity in conserved and commercial pig breeds in Poland.

Runs of homozygosity (ROH) are continuous segments of the genome that arose as a result of inbreeding, resulting in the inheritance of identical haplotypes from both parents who shared a common ancestor. In the present study, we performed a detailed characterization and comparison of ROH in four pig breeds, including intensively selected Polish Landrace as well as native unselected animals of Pulawska and two Zlotnicka breeds (White and Spotted). We used a medium-density PorcineSNP60 BeadChip assay (Illumina) and cgaTOH software to detect ROH covering a minimum of 30 adjacent SNPs and maintaining a size over 1 Mb. By analysing ROH distribution and frequency across the genome, we also identified genomic regions with high ROH frequency (so-called "ROH hotspots"). The obtained results showed that the analysed conserved breeds were characterized by a higher ROH span and higher ROH-based inbreeding coefficients (FROH ), which likely result from past population bottlenecks, increasing the overall inbreeding level within these populations. The analysis of ROH distribution across the genomes revealed the presence of both shared and breed-specific ROH hotspots. These hotspots, presumably representing genome regions under selection, overlapped with a variety of genes associated with processes connected with immune system functioning, reproduction, glucose homeostasis and metabolism. The genome regions with ROH hotspots overlapping in all analysed populations, located on SSC4 (51.9-55.9 Mb) and 13 (92.6-97.8 Mb), covered thirty-one different genes, including MMP16, SLC7A13, ATP6V0D2, CNGB3, WWiP1, RiMDN1 and CPNE3. These genes are primarily associated with biological regulation and metabolism, processes that could be responsible for the variety of the selected production and functional features.

Genetic Differentiation of the Two Types of Polish Cold-blooded Horses Included in the National Conservation Program.

The current role of the horse as a companion animal resulted in a decrease of interest in breeding and usage of draft horses. This meant that the population of cold-blooded horses in Poland has been dramatically reduced during the last decades. To avoid impoverishment of the gene pool of the local horse population, a conservation program was established which involves draft horses and other local horse breeds. The draft horses bred in Poland can be subdivided in a few horse types of which the most widespread and consolidated are Sztumski and Sokólski horses. These two subpopulations are phenotypically diversified, however, the overall level of their genetic differentiation seems to be relatively low and not precisely determined, especially with the use of molecular markers. In reference to this, in this study we used Illumina genotyping arrays to describe in detail the genetic differentiation of these two cold-blooded horse populations. We describe the genetic distance between them, as well as within-population variation, admixture patterns, and level of relatedness within populations. We also made an attempt to detect genome regions divergently selected between those horses by the detection of diversifying selection signals. The results of this study provide initial evidence supporting breeding decisions that were made during conservation breeding program design and answer questions raised by the breeders of Sokólski and Sztumski horses concerning the level of their genetic variation and differentiation.

A genome-wide scan for diversifying selection signatures in selected horse breeds.

The genetic differentiation of the current horse population was evolutionarily created by natural or artificial selection which shaped the genomes of individual breeds in several unique ways. The availability of high throughput genotyping methods created the opportunity to study this genetic variation on a genome-wide level allowing detection of genome regions divergently selected between separate breeds as well as among different horse types sharing similar phenotypic features. In this study, we used the population differentiation index (FST) that is generally used for measuring locus-specific allele frequencies variation between populations, to detect selection signatures among six horse breeds maintained in Poland. These breeds can be classified into three major categories, including light, draft and primitive horses, selected mainly in terms of type (utility), exterior, performance, size, coat color and appearance. The analysis of the most pronounced selection signals found in this study allowed us to detect several genomic regions and genes connected with processes potentially important for breed phenotypic differentiation and associated with energy homeostasis during physical effort, heart functioning, fertility, disease resistance and motor coordination. Our results also confirmed previously described association of loci on ECA3 (spanning LCORL and NCAPG genes) and ECA11 (spanning LASP1 gene) with the regulation of body size in our draft and primitive (small size) horses. The efficiency of the applied FST-based approach was also confirmed by the identification of a robust selection signal in the blue dun colored Polish Konik horses at the locus of TBX3 gene, which was previously shown to be responsible for dun coat color dilution in other horse breeds. FST-based method showed to be efficient in detection of diversifying selection signatures in the analyzed horse breeds. Especially pronounced signals were observed at the loci responsible for fixed breed-specific features. Several candidate genes under selection were proposed in this study for traits selected in separate breeds and horse types, however, further functional and comparative studies are needed to confirm and explain their effect on the observed genetic diversity of the horse breeds.

Genotyping-by-sequencing performance in selected livestock species.

Application of next generation sequencing for large scale genotyping in livestock is limited by high costs and challenging data analysis process. However, available restriction enzyme-based enrichment techniques like e.g. genotyping-by-sequencing (GBS) are promising tools allowing reduction of financial outlies by a high sample multiplexing and narrowing down the sequenced genome areas to the randomly distributed read tags. In this study, we tested the performance of standard, PstI endonuclease-adapted GBS protocol for population genetics in cattle, horse and sheep with application of different, including low-depth sequencing setups. It was found that the detected SNPs display desirable polymorphism parameters and are evenly scattered across the whole genome including gene coding regions. It was also shown that the SNPs can be successfully applied in population genetics, revealing the genetic differentiation of the studied breeds. The GBS approach represents a cost-effective alternative to existing genotyping methods which may find adoption in various research applications.

A genome-wide detection of selection signatures in conserved and commercial pig breeds maintained in Poland.

BACKGROUND: Identification of selection signatures can provide a direct insight into the mechanism of artificial selection and allow further disclosure of the candidate genes related to the animals' phenotypic variation. Domestication and subsequent long-time selection have resulted in extensive phenotypic changes in domestic pigs, involving a number of traits, like behavior, body composition, disease resistance, reproduction and coat color. In this study, based on genotypes obtained from PorcineSNP60 Illumina assay we attempt to detect both diversifying and within-breed selection signatures in 530 pigs belonging to four breeds: Polish Landrace, Pulawska, Zlotnicka White and Zlotnicka Spotted, of which the last three are a subject of conservative breeding and substantially represent the native populations. RESULTS: A two largely complementary statistical methods were used for signatures detection, including: pairwise FST and relative extended haplotype homozygosity (REHH) test. Breed-specific diversifying selection signals included several genes involved in processes connected with fertility, growth and metabolism which are potentially responsible for different phenotypes of the studied breeds. The diversifying selection signals also comprised PPARD gene that was previously found to have a large effect on the shape of the external ear in pigs or two genes encoding neuropeptide Y receptors (Y2 and Y5) involved in fat deposition and stress response which are important features differentiating the studied breeds. REHH statistics allowed detecting several within-breed selection signatures overlapping with genes connected with a range of functions including, among others: metabolic pathways, immune system response or implantation and development of the embryo. CONCLUSIONS: The study provides many potential candidate genes with implication for traits selected in the individual breeds and gives strong basis for further studies aiming at identification of sources of variation among the studied pig breeds.

Genomic landscape of copy number variation and copy neutral loss of heterozygosity events in equine sarcoids reveals increased instability of the sarcoid genome.

Although they are the most common neoplasms in equids, sarcoids are not fully characterized at the molecular level. Therefore, the objective of this study was to characterize the landscape of structural rearrangements, such as copy number variation (CNV) and copy neutral loss of heterozygosity (cnLOH), in the genomes of sarcoid tumor cells. This information will not only broaden our understanding of the characteristics of this genome but will also improve the general knowledge of this tumor and the mechanisms involved in its generation. To this end, Equine SNP64K Illumina microarrays were applied along with bioinformatics tools dedicated for signal intensity analysis. The analysis revealed increased instability of the genome of sarcoid cells compared with unaltered skin tissue samples, which was manifested by the prevalence of CNV and cnLOH events. Many of the identified CNVs overlapped with the other research results, but the simultaneously observed variability in the number and sizes of detected aberrations indicated a need for further studies and the development of more reliable bioinformatics algorithms. The functional analysis of genes co-localized with the identified aberrations revealed that these genes are engaged in vital cellular processes. In addition, a number of these genes directly contribute to neoplastic transformation. Furthermore, large numbers of cnLOH events identified in the sarcoids suggested that they may play no less significant roles than CNVs in the carcinogenesis of this tumor. Thus, our results indicate the importance of cnLOH and CNV in equine sarcoid oncogenesis and present a direction of future research.