Gene Expression Profile in Similar Tissues Using Transcriptome Sequencing Data of Whole-Body Horse Skeletal Muscle

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.

Identification of Differentially Expressed Genes after Endurance Runs in Karbadian Horses to Determine Candidates for Stress Indicators and Performance Capability

RNA sequencing makes it possible to uncover genetic mechanisms that underlie certain performance traits. In order to gain a deeper insight into the genetic background and biological processes involved in endurance performance in horses, the changes in the gene expression profiles induced by endurance runs over long (70 km) and short (15 km) distances in the blood of Kabardian horses (Equus caballus) were analyzed. For the long-distance runs, we identified 1484 up- and 691 downregulated genes, while after short-distance runs, only 13 up- and 8 downregulated genes (FC > |1.5|; p < 0.05) were found. These differentially expressed genes (DEGs) are involved in processes and pathways that are primarily related to stress response (interleukin production, activation of inflammatory system) but also to metabolism (carbohydrate catabolic process, lipid biosynthesis, NADP metabolic process). The most important genes involved in these processes therefore represent good candidates for the monitoring and evaluation of the performance of horses in order to avoid excessive demands when endurance performance is required, like ACOD1, CCL5, CD40LG, FOS, IL1R2, IL20RA, and IL22RA2, on the one hand, and, on the other hand, for assessing the suitability of a horse for endurance races, like GATA2, GYG1, HIF1A, MOGAT1, PFKFB3, PLIN5, SIK1, and STBD1.

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.