Search

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.

DNA methylation patterns and their relation with genetic polymorphisms were determined in the equine OAS1 locus. Genetic variants of OAS1 were previously found to be associated with susceptibility to West Nile virus infections in horses. The subject of the study were white blood cells of 13 juvenile and 13 old horses from AA and HC breed and a set of solid tissues from a single adult horse. The aim was to determine the degree of variation of CpG methylation profiles with concern for tissue type, horse breed and age. Results of direct BSPCR and cloned BSPCR sequencing revealed that all of determined CpG islands (CGIs) were hypermethylated in exception to CGI covering OAS1 promoter and exon 1. One of intragenic CGIs displayed variability of methylation patterns across eight tissue types. The variability of particular sub-types of white blood cells between AA and HC horses were considered as the possible cause of interbreed differences of methylation levels. Comparison of sequence variability between converted and unconverted DNAs of both horse breeds showed polymorphisms of CpG sites to be the source of monoallelic methylation in exception to the polymorphic CpGs located in the OAS1 promoter. Two of them are new polymorphic variants in the OAS1 promoter region. Application of methylation data in conjunction with genetic variation detected at the OAS1 locus might be useful to deepen the knowledge about mechanisms underlying immunity to viral infections in the horse.

In recent years, characterization of a copy number variation (CNV) of the genomic DNA has provided evidence for the relationship of this type of genetic variation with the occurrence of a broad spectrum of diseases, including cancer lesions. Copy number variants (CNVs) also occur in the genomes of healthy individuals as a result of abnormal recombination processes in germ cells and have a hereditary character contributing to the natural genetic diversity. Recent image analysis methods and advanced computational techniques allow for identification of CNVs using SNPs genotyping microarrays based on the analysis of signal intensity observed for markers located in the specific genomic regions. In this study we used CanineHD BeadChip assay (Illumina) to identify both natural and cancer-induced CNVs in the genomes of different dog breeds and in different cancer types occurring in this species. The obtained results showed that structural aberrations are a common phenomenon arising during a tumor progression and are more complex and widespread in tumors of mesenchymal tissue origin than in epithelial tissue originating tumors. The tumor derived CNVs, in comparison to healthy samples, were characterized by larger sizes of regions, higher number of amplifications, and in some cases encompassed genes with potential effect on tumor progression.

Animal genomics is currently undergoing dynamic development, which is driven by the flourishing of high-throughput genome analysis methods. Recently, a large number of animals has been genotyped with the use of whole-genome genotyping assays in the course of genomic selection programmes. The results of such genotyping can also be used for studies on different aspects of livestock genome functioning and diversity. In this article, we review the recent literature concentrating on various aspects of animal genomics, including studies on linkage disequilibrium, runs of homozygosity, selection signatures, copy number variation and genetic differentiation of animal populations. Our work is aimed at providing insight into certain achievements of animal genomics and to arouse interest in basic research on the complexity and structure of the genomes of livestock.

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.

MicroRNAs (miRNAs) are a class of small, noncoding RNAs, which play a vital role in the regulation of gene expression by binding to the 3' untranslated region (3'UTR) of a target mRNA. Despite a significant improvement in the identification of miRNAs in a variety of species, the coverage of the porcine miRNAome is still scarce. To identify porcine miRNAs potentially regulating processes taking place in the liver, we applied next generation sequencing. As a result, we detected 206 distinct miRNAs, of which 68 represented potential novel miRNAs. Among these new miRNAs, there were miRNAs deriving from the opposite arm of a hairpin precursor of already known miRNAs. Moreover, we observed 3' and 5' length and sequence variants, probably constituting so called isomiRs, as well as differentially mapped precursor loci, alternative precursor sequences and clustering of miRNA encoding genes. On the basis of expression levels, reflected by the number of sequence reads, we identified the most abundant miRNAs followed by gene target prediction and pathway analysis. The enriched pathways were connected with cellular and metabolic processes, growth factors as well as enzymatic activity. The obtained results are the first ones to concern the porcine liver miRNAome. Consequently, they will increase the number of known porcine miRNAs and facilitate further research on gene regulation mechanisms as well as biological processes associated with the liver functioning in pigs.