Transcriptomic profiling of mare endometrium at different stages of endometrosis.

In the current study, transcriptome profiles of mare endometrium, classified into categories I, IIA, and IIB according to Kenney and Doig, were compared using RNA sequencing, analyzed, and functionally annotated using in silico analysis. In the mild stage (IIA) of endometrosis compared to category I endometrium, differentially expressed genes (DEGs) were annotated to inflammation, abnormal metabolism, wound healing, and quantity of connective tissue. In the moderate stage (IIB) of endometrosis compared to category I endometrium, DEGs were annotated to inflammation, fibrosis, cellular homeostasis, mitochondrial dysfunction, and pregnancy disorders. Ingenuity pathway analysis (IPA) identified cytokines such as transforming growth factor (TGF)-ß1, interleukin (IL)-4, IL-13, and IL-17 as upstream regulators of DEGs associated with cellular homeostasis, metabolism, and fibrosis signaling pathways. In vitro studies showed the effect of these cytokines on DEGs such as ADAMTS1, -4, -5, -9, and HK2 in endometrial fibroblasts at different stages of endometrosis. The effect of cytokines on ADAMTS members' gene transcription in fibroblasts differs according to the severity of endometrosis. The identified transcriptomic changes associated with endometrosis suggest that inflammation and metabolic changes are features of mild and moderate stages of endometrosis. The changes of ADAMTS-1, -4, -5, -9, in fibrotic endometrium as well as in endometrial fibroblast in response to TGF-ß1, IL-4, IL-13, and IL-17 suggest the important role of these factors in the development of endometrosis.

Analysis of morphological disorders and ploidy in domestic cat blastocysts.

This study describes, for the first time, the relationship between morphology and ploidy in domestic cat embryos. Blastocyst morphology and quality were assessed using time-lapse recordings, while ploidy was analyzed using fluorescence in situ hybridization. Out of 54 blastocysts, clear fluorescence signals for all the molecular probes used were observed in 24 (44.4%) blastocysts, while in another 14 (25.9%) blastocysts, fluorescence signals only allowed for sex assessment. No clear signals were observed in the remaining 16 blastocysts (29.7%). Of the 24 blastocysts with clear signals, normal ploidy was detected in 10 (41.4%), 7 (29.2%) were diagnosed as haploid, and the remaining 7 blastocysts (29.2%) were mosaics. Additionally, results showed the distribution of diploid, haploid, and mosaic blastocysts in relation to the occurrence of morphological disorders and to embryo quality. The presence of abnormal embryo morphology and karyotype disorders may affect further development and the pregnancy rate. Due to the comparable proportion of good and poor quality blastocysts with disturbed ploidy, it is important to implement new methods of embryo assessment, especially when techniques used in humans, such as pronuclear observation, cannot be used.

Screening and detection of chromosomal copy number alterations in the domestic horse using SNP-array genotyping data.

Chromosomal abnormalities are a common cause of infertility in horses. However, they are difficult to detect using automated methods. Here, we propose a simple methodology based on single nucleotide polymorphism (SNP)-array data that allows us to detect the main chromosomal abnormalities in horses in a single procedure. As proof of concept, we were able to detect chromosomal abnormalities in 33 out of 268 individuals, including monosomies, chimerisms, and male and female sex-reversions, by analyzing the raw signal intensity produced by an SNP array-based genotyping platform. We also demonstrated that the procedure is not affected by the SNP density of the array employed or by the inbreeding level of the individuals. Finally, the methodology proposed in this study could be performed in an open bioinformatic environment, thus permitting its integration as a flexible screening tool in diagnostic laboratories and genomic breeding programs.

Severe asynapsis in spermatocytes of interspecific hybrids of the silver fox (Vulpes vulpes) and the blue fox (Alopex lagopus) leads to pachytene I arrest as a result of sustained H2AXγ phosphorylation.

Infertility is frequently associated with meiotic anomalies which can result in the production of chromosomally abnormal gametes or be concomitant with meiotic arrest. We investigated whether spermatocytes of male interspecific hybrids of the red fox (Vulpes vulpes) and the arctic fox (Alopex lagopus) presented alterations in chromosomal synapses and meiotic checkpoint signalling. Using the immunofluorescence technique with SP1 and SP3 proteins, bivalent structures and their deviations (multivalents, univalents and not fully conjugated bivalents) were analyzed on meiotic preparations. This technique allowed the localization of frequent foci of phosphorylated histones H2AHγ (Ser 139) to the meiotic block in late pachytene. These results indicate a disruption of meiotic division in male fox hybrids, which leads to a high percentage of apoptotic cells in the gonads of these animals and, consequently, sterility.

Methylation Marks of Blood Leukocytes of Native Hucul Mares Differentiated in Age.

Horses are one of the longest-living species of farm animals. Advanced age is often associated with a decrease in body condition, dysfunction of immune system, and late-onset disorders. Due to this, the search for new solutions in the prevention and treatment of pathological conditions of the advanced age of horses is desirable. That is why the identification of aging-related changes in the horse genome is interesting in this respect. In the recent years, the research on aging includes studies of age-related epigenetic effects observed on the DNA methylation level. We applied reduced representation bisulfite sequencing (RRBS) to uncover a range of age DMR sites in genomes of blood leukocytes derived from juvenile and aged horses of native Hucul breed. Genes colocated with age-related differentially methylated regions (age DMRs) are the members of pathways involved in cellular signal transduction, immune response, neurogenesis, differentiation, development, and cancer progression. A positive correlation was found between methylation states and gene expression in particular loci from our data set. Some of described age DMR-linked genes were also reported elsewhere. Obtained results contribute to the knowledge about the molecular basis of aging of equine blood cells.

Chondrogenic expression and DNA methylation patterns in prolonged passages of chondrocyte cell lines of the horse.

We investigated the activity of chondrogenic markers and variation of methylation patterns in equine cartilaginous cells cultivated in monolayer. The transcriptional and epigenetic effect of the long-term culture of chondrocytes has been evaluated using several passages of chondrocyte cell-lines derived from equine articular cartilage. Using 3 genes as endogenous control we tested the expression of 7 genes important for different stages of chondrocyte differentiation and maturation. CpG islands in RUNX3 locus were inspected for the evaluation of differential methylation state of passaged cell-lines. The general decline of transcript abundance of marker loci was detected in passage 11 which is the sign of dedifferentiation of cultivated chondrocytes in prolonged monolayer culture. Passages 13 and 14 were characterized by the upregulation of a number of genes, possibly due to the heterogeneity of developed cell lines at this stage of the culture. Instead, gradual increase of methylation percent at particular CpG sites of RUNX3 locus was associated with the growing number of passage. This finding led us to the conclusion that epigenetic alterations better describe the stage of cultivated chondrocytes.

Characterisation of genome-wide structural aberrations in canine mammary tumours using single nucleotide polymorphism (SNP) genotyping assay.

Characterisation of copy number variation (CNV) and loss of heterozygosity (LOH) has pro- vided evidence for the relationship of this type of genetic variation with the occurrence of a broad spectrum of diseases, including cancer lesions. The role of CNVs and germinal or somatic LOHs in canine mammary tumours is still unknown. Therefore, the aim of this study was to identify CNVs and LOHs in canine mammary tumours. Forty-eight samples obtained from normal (n=24) and tumour (n=24) tissues of dogs were analysed. In the study, we used CanineHD BeadChip assay (Illumina) and OncoSNP software to identify copy number alternations in genomes of dif- ferent dog breeds and in different mammary cancer types occurring in this species. The analyses revealed that, in the case of CNV, the amplification-type variants were longer and more frequent than deletions. Based on the analysis of the frequency of different types of aberrations in the in- dividual parts of the genome, regions that are particularly susceptible to structural aberrations were indicated. The fraction of genes identified within these regions was associated with major processes of neoplastic transformation. Association analysis of such traits as tumour grading as well as the size and age of dogs demonstrated that structural aberrations were more frequent in dogs diagnosed with tumour malignancy grade II and III, in dogs with a larger body size, and in large dogs aged 7-8. The promising results of these pioneering investigations prompt continuation thereof to analyse other types of cancer.

Diversifying selection signatures among divergently selected subpopulations of Polish Red cattle.

Polish Red cattle is one of the few indigenous breeds of European red cattle which is characterized by several desired features, such as high disease resistance, good health, longevity, good fertility, and high nutritional value of milk. Currently, Polish Red cattle population is a subject of two independent breeding programs: (i) improvement program and (ii) genetic resources conservation program. The aim of the improvement program is the genetic progress in terms of milk production and body conformation traits, while the conservation program mainly focuses on protection of the genetic resources of Polish Red cattle and preservation of the existing, original gene pool. By the analysis of FST genetic distances across genome-wide SNP panel, we detected diversifying selection signatures among these two subpopulations and indicated (among others) the significance of DGAT1 and FGF2 genes for milk production traits in these cattle. We also found that among genes being presumably under selection in terms of milk production, there are genes responsible, for example, for mammary gland development (e.g., SOSTDC1, PYGO2, MED1, and CCND1) and immune system response (e.g., IL10RA, IL12B, and IL21). The most important finding of this study is that the most pronounced genetic differences between the analyzed populations were associated with ß-defensin genes (e.g., DEFB1, DEFB4A, DEFB5, DEFB7, DEFB10, DEFB13, EBD, BNBD-6, and LAP) located within so-called bovine cluster D on BTA27. The ß-defensins are expressed mainly in the mammary gland and are antimicrobial peptides against the Gram-negative and Gram-positive bacteria, viruses, and other unicellular parasites. This suggests that antimicrobial resistance of mammary gland is of high importance during selection towards increased milk production and that genes responsible for this process are selected together with increasing levels of productivity.

DNA methylation patterns of the S100A14, POU2F3 and SFN genes in equine sarcoid tissues.

Genetic and epigenetic alterations in the equine sarcoid, a locally invasive skin tumour of equids, are still poorly characterized. Numerous studies have provided reliable evidence for the relationship between the development of cancer and the loss of function of a number of tumour suppressor genes. In the present study, we assessed methylation levels in the promoter region of SFN, S100A14 and POU2F3 genes in sarcoid samples to clarify whether DNA methylation may be associated with previously identified changes in the expression level of these genes during the course of tumour progression. Using bisulfite sequencing and clone sequencing, we detected that lesional samples had a significantly higher rate of DNA methylation in the analyzed S100A14A region than the corresponding normal skin tissue. A frequent methylation of the SFN and POU2F3 promoter sequences were observed in both the tumour samples and the control skin tissues. Further studies are needed to evaluate the role of aberrant methylation in sarcoid progression and to understand the mechanisms involved in reduced expression of SFN, S100A14 and POU2F3 genes in the lesional tissues.

Molecular characterization of the apoptosis-related SH3RF1 and SH3RF2 genes and their association with exercise performance in Arabian horses.

BACKGROUND: Apoptosis plays an important role in the regulation of healthy tissue growth and development as well as in controlling the maintenance of homeostasis in exercising muscles. During an intensive physical effort, the regulation of cell death by apoptosis results in the replacement of unaccustomed muscle cells by new cells that are better suited to exercise. The aim of this study was to determine the expression of two genes (SH3FR1 and SH3RF2) that control apoptosis in muscle tissues during training periods characterized by different intensities. The gene expression levels were estimated using real-time PCR method in skeletal muscle biopsies collected from 15 Arabian horses (untrained, after an intense gallop phase, and at the end of the racing season). An association study was performed on 250 Arabian horses to assess the effect of the SH3RF2:c.796 T > C (p.Ser266Pro) variant on race performance traits in flat gallop-racing. RESULTS: A gene expression analysis confirmed a significant decrease (p < 0.01) in the anti-apoptotic SH3RF2 (POSHER) gene during training periods that differed in intensity. The highest SH3RF2 expression level was detected in the muscles of untrained horses, whereas the lowest expression was identified at the end of the racing season in horses that were fully adapted to the exercise. A non-significant decrease in SH3RF1 gene expression following the training periods was observed. Moreover, a serine substitution by proline at amino acid position 266 (CC genotype) was negatively associated with the probability of winning races, the number of races in which a horse occurred and the financial value of the prizes. Horses with the TT genotype achieved the highest financial benefits, both for total winnings and for winnings per race in which the horses participated. CONCLUSIONS: The present study showed the supposed regulation mechanism of exercise-induced apoptosis in horses at the molecular level. The identified SH3RF2: c.796 T > C missense variant was associated with selected racing performance traits, which is important information during the evaluation of horses' exercise predisposition. The association results and frequencies of the CT and TT genotypes suggest the possibility of using SH3RF2 variant in selection to improve the racing performance of Arabian horses.

Transcriptomic hallmarks of bone remodelling revealed by RNA-Seq profiling in blood of Arabian horses during racing training regime.

The impact of exercises on young developing organisms is still of interest to researchers. Similarly like Thoroughbreds, Arabian horses competing at the race track. The high percent of lameness and loss of days in training are often the result of weakness in the condition of the musculoskeletal system. The objective of the presented study was to identify by RNA-Seq method, the possible skeletal system originating transcriptomic profile in peripheral blood of Arabian horses undergoing race training. Obtained results showed that one of the most significantly deregulated pathway involved in bone homeostasis was those involved in osteoclast differentiation. Among the significantly expressed molecules, we recognized twelve genes potentially involved in the metabolism of the skeletal system: BGLAP, CTSK, TYROBP, PDLIM7, SLC9B2, TWSG1, NOTCH2, IL6ST, VAV3, NFATc1, CLEC5A, TXLNG. The panel of identified genes should be evaluated as candidate biomarkers for bone homeostasis indicators of Arabians performing on race tracks to assess bone remodelling states during training for race track competitions.

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

OBJECTIVE: In this study, for the first time we report the genome-wide DNA methylation profile of skin tumour in horses and describe differentially methylated genomic regions (DMRs) with respect to healthy skin. MATERIALS & METHODS: The comparative analysis of DNA methylation patterns detected using Reduced Representation Bisulfite Sequencing (RRBS) technique, allowed identification of 136 regions showing differential methylation between sarcoid and normal skin tissue. RESULTS: Most of the identified DMRs were short fragments, less than 1 kb in size, located in the intergenic regions. Among identified DMRs there were also regions located within genes directly or indirectly related with oncogenesis. We additionally validated 9 CpG sites showing hypomethylation and 9 CpG sites that were hypermethylated in lesional sample, confirming the identified changes in the DNA methylation. CONCLUSION: Knowledge on the changes taking place in the process of DNA methylation may provide a basis for the development of new alternative diagnostic or therapeutic approaches to equine sarcoids.

First case of sterility associated with sex chromosomal abnormalities in a jenny.

Chromosomal abnormalities are one of the main causes of genetic infertility in horses. Currently, their detection rate is rising due to the use of new diagnostic tools employing molecular markers linked to the sex chromosome pair. Despite genetic similarities, there are no previous reports of sterility associated with chromosomal abnormalities in the domestic donkey (Equus asinus). Hereby, we determined the presence of a chromosomal mosaicism in a female donkey with reproductive problems using molecular methodologies developed for horses. A two-and-a-half-year-old jenny characterized by morphological abnormalities of the reproductive tract was cytogenetically analysed using conventional and fluorescent techniques and a group of microsatellite markers (short tandem repeat, STR). At the same time, five ultrasound measures of the reproductive tract were taken and compared with eight contemporary jennies of the same breed. After slaughter, morphological examinations showed that the case study had a blind vaginal vestibule defining an empty pouch that covered the entrance of the cervical os. Histopathological studies demonstrated that this abnormal structure was compatible with a remnant hymen. Molecular markers, STR and fluorescent in situ hybridization determinations revealed that the animal was a 62, XX/61,X mosaic and, therefore, the first case of chromosomal abnormalities in the sex pair reported in donkeys.

Transcriptome analysis of equine sarcoids.

Equine sarcoids are the most commonly detected skin tumours in Equidae. In the present research, a comparative transcriptomic analysis was performed which aimed at looking inside a tumour biology and identification of the expression profile as a potential source of cancer specific genes useful as biomarkers. We have used Horse Gene Expression Microarray data from matched equine sarcoids and tumour-distant skin samples. In total, 901 significantly differentially expressed genes (DEGs) between lesional and healthy skin samples have been identified (fold change ≥ 2; P < 0.05). The large subset of DEGs, with decreased expression, was associated with a suppression of malignant transformation, whereas several overexpressed genes were involved in the processes associated with growth and progression of a tumour or immune system activity. Our results, as a first to date, showed comprehensive transcriptome analysis of skin tumour in horses and pinpointed significant pathways and genes related with oncogenesis processes.

Relevance of Molecular Changes in the ND4 Gene in German Shepherd Dog Tumours.

The aim of the study was to identify polymorphisms and mutations in the mitochondrial ND4 gene and to analyse the associations between the occurrence of molecular changes in mtDNA and phenotypic traits in tumours in German Shepherd dogs. Fifty samples obtained from blood and tumour tissues of German Shepherd dogs with diagnosed tumours were analysed. DNA extraction, amplification, and sequencing of the mtDNA ND4 gene, and bioinformatics, statistical, and in silico protein coding SNP analyses were performed. ND4 mutations and/or polymorphisms were noted in eleven nucleotide positions in nearly half of the examined dogs. All the changes were substitution mutations. A majority of the changes identified were homoplasmic. In one dog with osteosarcoma, blood heteroplasmy was detected. In two positions of the ND4 gene, presence of non-synonymous mutations leading to amino acid changes in the ND4 protein was reported. Analyses carried out to determine the deleterious effect of mutations indicated an almost 97 and 62% probability that a single amino acid substitution (p.G239V and p.I401T, respectively) in the protein has a negative impact on its function. The results of statistical analyses indicate a significant association between the occurrence of mutations in three loci of the ND4 gene and the location of tumours. The mutations identified may be a result of cell adaptation to the changes in the environment occurring during carcinogenesis. The high frequency of mutations in the tumours may indicate genetic instability of mtDNA, which may also play a role in carcinogenesis.

The use of runs of homozygosity for estimation of recent inbreeding in Holstein cattle.

Controlling inbreeding in livestock populations is of great importance because excess relatedness among animals leads to a rapid loss of genetic variation and to adverse phenotypical effects associated with an inbreeding depression. Recent advances in genotyping technology have made it possible to study inbreeding at a molecular level by the analysis of genome-wide single nucleotide polymorphism panels. In this study, we used BovineSNP50 assay (Illumina) to estimate genomic inbreeding coefficient in 298 Holstein cattle by the analysis of the genome portion in runs of homozygosity (FROH) or using genomic relationship matrix (FGRM), and compared this data with conventional pedigree-based inbreeding coefficients (FPED). Weak or moderate Spearman's rank correlations were observed between FROH and FPED which depended on the ROH length categories used for calculations and inclusion of animals with different number of complete generations registered in pedigrees. The highest correlations were observed when using ROH with lengths over 8 Mb (0.334). The correlations tended to increase as pedigree depth increased, and were the highest for animals with seven complete generations of pedigree data. FGRM correlated poorly with pedigree-based estimates, which suggests that ROH-based inbreeding coefficients better reflect recent relatedness among animals.

Age-related methylation profiles of equine blood leukocytes in the RNASEL locus.

Methylation profiles across three CpG islands of the RNASEL gene were determined in blood leukocyte samples of Anglo-Arabian and Hucul horses. Bisulfite sequencing revealed hypomethylated state of the RNASEL promoter coinciding with methylated CpG island placed inside the gene. Several CpG sites were identified for which the methylation state was influenced by DNA polymorphism. Two of them showed monoallelic methylation. One of the CpG sites revealed functional polymorphism. A number of partially methylated CpG sites have been observed in the promoter area of RNASEL, which were used for the comparison of breed- and age-related effects. Clone bisulfite sequencing of blood leukocyte samples collected at different ages from particular individuals of AA and HC breeds and, also, BSPCR sequencing of 50 samples of juvenile and old AA and HC horses revealed increased methylation in particular CpG sites during aging. The age-related heterogeneity of white blood cells was hypothesized as being one of the potential causes of observed variability of methylation profiles in the RNASEL promoter.

Identification of genome-wide selection signatures in the Limousin beef cattle breed.

The study is aimed at identifying selection footprints within the genome of Limousin cattle. With the use of Extended Haplotype Homozygosity test, supplemented with correction for variation in recombination rates across the genome, we created map of selection footprints and detected 173 significant (p < 0.01) core haplotypes being potentially under positive selection. Within these regions, a number of candidate genes associated inter alia with skeletal muscle growth (GDF15, BMP7, BMP4 and TGFB3) or postmortem proteolysis and meat maturation (CAPN1 and CAPN5) were annotated. Noticeable clusters of selection footprints were detected on chromosomes 1, 4, 8 and 14, which are known to carry several quantitative trait loci for growth traits and meat quality. The study provides information about the genes and metabolic pathways potentially modified under the influence of directional selection, aimed at improving beef production characteristics in Limousin cattle.

Genome-wide characteristics of copy number variation in Polish Holstein and Polish Red cattle using SNP genotyping assay.

Copy number variation (CNV), which results from deletions or amplifications of large fragments of genomic DNA, is widespread in mammalian genomes and apart from its potential pathogenic effect it is considered as a source of natural genetic diversity. In cattle populations, this kind of genetic variability remains still insufficiently elucidated and studies focusing on the detection of new structural genomic variants in different cattle populations may contribute to a better understanding of cattle breeds' diversity and genetic basis of production traits. In this study, by using BovineSNP50 assay and cnvPartition algorithm we identified CNVs in two different cattle breeds: Holstein (859 animals) and Polish Red (301). In Holstein cattle we found 648 CNVs which could be reduced to 91 non-redundant variable genomic regions (CNVRs) covering in total 168.6 Mb of the genomic sequence. In Polish Red cattle we detected 62 CNVs, localized in 37 variable regions encompassing 22.3 Mb of the sequence, corresponding to 0.89 % of the autosomal genome. Within the regions we identified 1,192 unique RefSeq genes which are engaged in a variety of biological processes. High concordance of the regions' distribution was found between the studied breeds, however copy number variants seemed to be more common in Holstein cattle. About 26 % of the regions described in this study could be classified as newly identified. The results of this study will broaden the knowledge of CNVs in genomes of cattle of different breeds and will provide foundations for further research aiming to identify a relationship between this type of genetic variation and phenotypic traits.