Expression genome-wide association study identifies key regulatory variants enriched with metabolic and immune functions in four porcine tissues.

BACKGROUND: Integration of high throughput DNA genotyping and RNA-sequencing data enables the discovery of genomic regions that regulate gene expression, known as expression quantitative trait loci (eQTL). In pigs, efforts to date have been mainly focused on purebred lines for traits with commercial relevance as such growth and meat quality. However, little is known on genetic variants and mechanisms associated with the robustness of an animal, thus its overall health status. Here, the liver, lung, spleen, and muscle transcriptomes of 100 three-way crossbred female finishers were studied, with the aim of identifying novel eQTL regulatory regions and transcription factors (TFs) associated with regulation of porcine metabolism and health-related traits. RESULTS: An expression genome-wide association study with 535,896 genotypes and the expression of 12,680 genes in liver, 13,310 genes in lung, 12,650 genes in spleen, and 12,595 genes in muscle resulted in 4,293, 10,630, 4,533, and 6,871 eQTL regions for each of these tissues, respectively. Although only a small fraction of the eQTLs were annotated as cis-eQTLs, these presented a higher number of polymorphisms per region and significantly stronger associations with their target gene compared to trans-eQTLs. Between 20 and 115 eQTL hotspots were identified across the four tissues. Interestingly, these were all enriched for immune-related biological processes. In spleen, two TFs were identified: ERF and ZNF45, with key roles in regulation of gene expression. CONCLUSIONS: This study provides a comprehensive analysis with more than 26,000 eQTL regions identified that are now publicly available. The genomic regions and their variants were mostly associated with tissue-specific regulatory roles. However, some shared regions provide new insights into the complex regulation of genes and their interactions that are involved with important traits related to metabolism and immunity.

An ABCA4 loss-of-function mutation causes a canine form of Stargardt disease.

Autosomal recessive retinal degenerative diseases cause visual impairment and blindness in both humans and dogs. Currently, no standard treatment is available, but pioneering gene therapy-based canine models have been instrumental for clinical trials in humans. To study a novel form of retinal degeneration in Labrador retriever dogs with clinical signs indicating cone and rod degeneration, we used whole-genome sequencing of an affected sib-pair and their unaffected parents. A frameshift insertion in the ATP binding cassette subfamily A member 4 (ABCA4) gene (c.4176insC), leading to a premature stop codon in exon 28 (p.F1393Lfs*1395), was identified. In contrast to unaffected dogs, no full-length ABCA4 protein was detected in the retina of an affected dog. The ABCA4 gene encodes a membrane transporter protein localized in the outer segments of rod and cone photoreceptors. In humans, the ABCA4 gene is associated with Stargardt disease (STGD), an autosomal recessive retinal degeneration leading to central visual impairment. A hallmark of STGD is the accumulation of lipofuscin deposits in the retinal pigment epithelium (RPE). The discovery of a canine homozygous ABCA4 loss-of-function mutation may advance the development of dog as a large animal model for human STGD.

A systems biology framework integrating GWAS and RNA-seq to shed light on the molecular basis of sperm quality in swine.

BACKGROUND: Genetic pressure in animal breeding is sparking the interest of breeders for selecting elite boars with higher sperm quality to optimize ejaculate doses and fertility rates. However, the molecular basis of sperm quality is not yet fully understood. Our aim was to identify candidate genes, pathways and DNA variants associated to sperm quality in swine by analysing 25 sperm-related phenotypes and integrating genome-wide association studies (GWAS) and RNA-seq under a systems biology framework. RESULTS: By GWAS, we identified 12 quantitative trait loci (QTL) associated to the percentage of head and neck abnormalities, abnormal acrosomes and motile spermatozoa. Candidate genes included CHD2, KATNAL2, SLC14A2 and ABCA1. By RNA-seq, we identified a wide repertoire of mRNAs (e.g. PRM1, OAZ3, DNAJB8, TPPP2 and TNP1) and miRNAs (e.g. ssc-miR-30d, ssc-miR-34c, ssc-miR-30c-5p, ssc-miR-191, members of the let-7 family and ssc-miR-425-5p) with functions related to sperm biology. We detected 6128 significant correlations (P-value ≤ 0.05) between sperm traits and mRNA abundances. By expression (e)GWAS, we identified three trans-expression QTL involving the genes IQCJ, ACTR2 and HARS. Using the GWAS and RNA-seq data, we built a gene interaction network. We considered that the genes and interactions that were present in both the GWAS and RNA-seq networks had a higher probability of being actually involved in sperm quality and used them to build a robust gene interaction network. In addition, in the final network we included genes with RNA abundances correlated with more than four semen traits and miRNAs interacting with the genes on the network. The final network was enriched for genes involved in gamete generation and development, meiotic cell cycle, DNA repair or embryo implantation. Finally, we designed a panel of 73 SNPs based on the GWAS, eGWAS and final network data, that explains between 5% (for sperm cell concentration) and 36% (for percentage of neck abnormalities) of the phenotypic variance of the sperm traits. CONCLUSIONS: By applying a systems biology approach, we identified genes that potentially affect sperm quality and constructed a SNP panel that explains a substantial part of the phenotypic variance for semen quality in our study and that should be tested in other swine populations to evaluate its relevance for the pig breeding sector.

A history of why fathers' RNA matters.

Having been debated for many years, the presence and role of spermatozoal RNAs is resolving, and their contribution to development is now appreciated. Data from different species continue show that sperm contain a complex suite of coding and noncoding RNAs that play a role in an individual's life course. Mature sperm RNAs provide a retrospective of spermatogenesis, with their presence and abundance reflecting sperm maturation, fertility potential, and the paternal contribution to the developmental path the offspring may follow.Sperm RNAs delivered upon fertilization provide some of the initial contacts with the oocyte, directly confront the maternal with the paternal contribution as a prelude to genome consolidation. Following syngamy, early embryo development may in part be modulated by paternal RNAs that can include epidydimal passengers. This provides a direct path to relay an experience and then initiate a paternal response to the environment to the oocyte and beyond. Their epigenetic impact is likely felt prior to embryonic genome activation when the population of sperm delivered transcripts markedly changes. Here, we review the insights gained from sperm RNAs over the years, the subtypes, and the caveats of the RNAs described. We discuss the role of sperm RNAs in fertilization and embryo development, and their possible mechanism(s) influencing offspring phenotype. Approaches to meet the future challenges as the study of sperm RNAs continues, include, elucidating the potential mechanisms underlying how paternal allostatic load, the constant adaptation of health to external conditions, may be relayed by sperm RNAs to affect future generations.

Genomic, transcriptomic and epigenomic analysis towards the understanding of porcine semen quality traits. Past, current and future trends.

The importance of boar reproductive traits, including semen quality, in the sustainability of pig production system is increasingly being acknowledged by academic and industrial sectors. Research is needed to understand the biology and genetic components underlying these traits so that they can be incorporated into selection schemes and managerial decisions. This article reviews our current understanding of genome biology and technologies for genome, transcriptome and epigenome analysis which now facilitate the identification of causal variants affecting phenotypes more than ever before. Genetic and transcriptomic analysis of candidate genes, Genome-Wide Association Studies, expression microarrays, RNA-Seq of coding and noncoding genes and epigenomic evaluations have been conducted to profile the molecular makeups of pig sperm. These studies have provided insightful information for a several semen-related parameters. Nonetheless, this research is still incipient. The spermatozoon harbors a reduced transcriptome and highly modified epigenome, and it is assumed to be transcriptionally silent for nuclear gene expression. For this reason, the extent to which the sperm's RNA and epigenome recapitulate sperm biology and function is unclear. Hence, we anticipate that single-cell level analyses of the testicle and other male reproductive organs, which can reveal active transcription and epigenomic profiles in cells influencing sperm quality, will gain popularity and markedly advance our understanding of sperm-related traits. Future research will delve deeper into sperm fertility, boar resilience to environmental changes or harsh conditions, especially in the context of global warming, and also in transgenerational inheritance and how the environment influences the sperm transcriptome and epigenome.

Micrococcal nuclease sequencing of porcine sperm suggests enriched co-location between retained histones and genomic regions related to semen quality and early embryo development.

The mammalian spermatozoon has a unique chromatin structure in which the majority of histones are replaced by protamines during spermatogenesis and a small fraction of nucleosomes are retained at specific locations of the genome. The sperm's chromatin structure remains unresolved in most animal species, including the pig. However, mapping the genomic locations of retained nucleosomes in sperm could help understanding the molecular basis of both sperm development and function as well as embryo development. This information could then be useful to identify molecular markers for sperm quality and fertility traits. Here, micrococcal nuclease digestion coupled with high throughput sequencing was performed on pig sperm to map the genomic location of mono- and sub-nucleosomal chromatin fractions in relation to a set of diverse functional elements of the genome, some of which were related to semen quality and early embryogenesis. In particular, the investigated elements were promoters, the different sections of the gene body, coding and non-coding RNAs present in the pig sperm, potential transcription factor binding sites, genomic regions associated to semen quality traits and repeat elements. The analysis yielded 25,293 and 4,239 peaks in the mono- and sub-nucleosomal fractions, covering 0.3% and 0.02% of the porcine genome, respectively. A cross-species comparison revealed positional conservation of the nucleosome retention in sperm between the pig data and a human dataset that found nucleosome enrichment in genomic regions of importance in development. Both gene ontology analysis of the genes mapping nearby the mono-nucleosomal peaks and the identification of putative transcription factor binding motifs within the mono- and the sub- nucleosomal peaks showed enrichment for processes related to sperm function and embryo development. There was significant motif enrichment for Znf263, which in humans was suggested to be a key regulator of genes with paternal preferential expression during early embryogenesis. Moreover, enriched positional intersection was found in the genome between the mono-nucleosomal peaks and both the RNAs present in pig sperm and the RNAs related to sperm quality. There was no co-location between GWAS hits for semen quality in swine and the nucleosomal sites. Finally, the data evidenced depletion of mono-nucleosomes in long interspersed nuclear elements and enrichment of sub-nucleosomes in short interspersed repeat elements.These results suggest that retained nucleosomes in sperm could both mark regulatory elements or genes expressed during spermatogenesis linked to semen quality and fertility and act as transcriptional guides during early embryogenesis. The results of this study support the undertaking of ambitious research using a larger number of samples to robustly assess the positional relationship between histone retention in sperm and the reproductive ability of boars.

Urinary microRNAome in healthy cats and cats with pyelonephritis or other urological conditions.

MicroRNAs (miRNAs) are short non-coding RNAs that regulate gene expression at the post-transcriptional level. miRNAs have been found in urine and have shown diagnostic potential in human nephropathies. Here, we aimed to characterize, for the first time, the feline urinary miRNAome and explore the use of urinary miRNA profiles as non-invasive biomarkers for feline pyelonephritis (PN). Thirty-eight cats were included in a prospective case-control study and classified in five groups: healthy Control cats (n = 11), cats with PN (n = 10), cats with subclinical bacteriuria or cystitis (SB/C, n = 5), cats with ureteral obstruction (n = 7) and cats with chronic kidney disease (n = 5). By small RNA sequencing we identified 212 miRNAs in cat urine, including annotated (n = 137) and putative novel (n = 75) miRNAs. The 15 most highly abundant urinary miRNAs accounted for nearly 71% of all detected miRNAs, most of which were previously identified in feline kidney. Ninety-nine differentially abundant (DA) miRNAs were identified when comparing Control cats to cats with urological conditions and 102 DA miRNAs when comparing PN to other urological conditions. Tissue clustering analysis revealed that the majority of urine samples clustered close to kidney, which confirm the likely cellular origin of the secreted urinary miRNAs. Relevant DA miRNAs were verified by quantitative real-time PCR (qPCR). Eighteen miRNAs discriminated Control cats from cats with a urological condition. Of those, seven miRNAs were DA by both RNAseq and qPCR methods between Control and PN cats (miR-125b-5p, miR-27a-3p, miR-21-5p, miR-27b-3p, miR-125a-5p, miR-17-5p and miR-23a-3p) or DA between Control and SB/C cats (miR-125b-5p). Six additional miRNAs (miR-30b-5p, miR-30c, miR-30e-5p, miR-27a-3p, miR-27b-39 and miR-222) relevant for discriminating PN from other urological conditions were identified by qPCR alone (n = 4) or by both methods (n = 2) (P<0.05). This panel of 13 miRNAs has potential as non-invasive urinary biomarkers for diagnostic of PN and other urological conditions in cats.

Exploring the ovine sperm transcriptome by RNAseq techniques. I Effect of seasonal conditions on transcripts abundance.

Understanding the cell molecular changes occurring as a results of climatic circumstances is crucial in the current days in which climate change and global warming are one of the most serious challenges that living organisms have to face. Sperm are one of the mammals' cells most sensitive to heat, therefore evaluating the impact of seasonal changes in terms of its transcriptional activity can contribute to elucidate how these cells cope with heat stress events. We sequenced the total sperm RNA from 64 ejaculates, 28 collected in summer and 36 collected in autumn, from 40 Manchega rams. A highly rich transcriptome (11,896 different transcripts) with 90 protein coding genes that exceed an average number of 5000 counts were found. Comparing transcriptome in the summer and autumn ejaculates, 236 significant differential abundance genes were assessed, most of them (228) downregulated. The main functions that these genes are related to sexual reproduction and negative regulation of protein metabolic processes and kinase activity. Sperm response to heat stress supposes a drastic decrease of the transcriptional activity, and the upregulation of only a few genes related with the basic functions to maintain the organisms' homeostasis and surviving. Rams' spermatozoids carry remnant mRNAs which are retrospectively indicators of events occurring along the spermatogenesis process, including abiotic factors such as environmental temperature.

Fine Mapping of a Major Backfat QTL Reveals a Causal Regulatory Variant Affecting the CCND2 Gene.

Backfat is an important trait in pork production, and it has been included in the breeding objectives of genetic companies for decades. Although adipose tissue is a good energy storage, excessive fat results in reduced efficiency and economical losses. A large QTL for backfat thickness on chromosome 5 is still segregating in different commercial pig breeds. We fine mapped this QTL region using a genome-wide association analysis (GWAS) with 133,358 genotyped animals from five commercial populations (Landrace, Pietrain, Large White, Synthetic, and Duroc) imputed to the porcine 660K SNP chip. The lead SNP was located at 5:66103958 (G/A) within the third intron of the CCND2 gene, with the G allele associated with more backfat, while the A allele is associated with less backfat. We further phased the QTL region to discover a core haplotype of five SNPs associated with low backfat across three breeds. Linkage disequilibrium analysis using whole-genome sequence data revealed three candidate causal variants within intronic regions and downstream of the CCND2 gene, including the lead SNP. We evaluated the association of the lead SNP with the expression of the genes in the QTL region (including CCND2) in a large cohort of 100 crossbred samples, sequenced in four different tissues (lung, spleen, liver, muscle). Results show that the A allele increases the expression of CCND2 in an additive way in three out of four tissues. Our findings indicate that the causal variant for this QTL region is a regulatory variant within the third intron of the CCND2 gene affecting the expression of CCND2.

Association of fecal and serum microRNA profiles with gastrointestinal cancer and chronic inflammatory enteropathy in dogs.

BACKGROUND: Reliable biomarkers to differentiate gastrointestinal cancer (GIC) from chronic inflammatory enteropathy (CIE) in dogs are needed. Fecal and serum microRNAs (miRNAs) have been proposed as diagnostic and prognostic markers of GI disease in humans and dogs. HYPOTHESIS/OBJECTIVES: Dogs with GIC have fecal and serum miRNA profiles that differ from those of dogs with CIE. AIMS: (a) identify miRNAs that differentiate GIC from CIE, (b) use high-throughput reverse transcription quantitative real-time PCR (RT-qPCR) to establish fecal and serum miRNA panels to distinguish GIC from CIE in dogs. ANIMALS: Twenty-four dogs with GIC, 10 dogs with CIE, and 10 healthy dogs, all client-owned. METHODS: An international multicenter observational prospective case-control study. Small RNA sequencing was used to identify fecal and serum miRNAs, and RT-qPCR was used to establish fecal and serum miRNA panels with the potential to distinguish GIC from CIE. RESULTS: The best diagnostic performance for distinguishing GIC from CIE was fecal miR-451 (AUC: 0.955, sensitivity: 86.4%, specificity: 100%), miR-223 (AUC: 0.918, sensitivity: 90.9%, specificity: 80%), and miR-27a (AUC: 0.868, sensitivity: 81.8%, specificity: 90%) and serum miR-20b (AUC: 0.905, sensitivity: 90.5%, specificity: 90%), miR-148a-3p (AUC: 0.924, sensitivity: 85.7%, specificity: 90%), and miR-652 (AUC: 0.943, sensitivity: 90.5%, specificity: 90%). Slightly improved diagnostic performance was achieved when combining fecal miR-451 and miR-223 (AUC: 0.973, sensitivity: 95.5%, specificity: 90%). CONCLUSIONS AND CLINICAL IMPORTANCE: When used as part of a diagnostic RT-qPCR panel, the abovementioned miRNAs have the potential to function as noninvasive biomarkers for the differentiation of GIC and CIE in dogs.

Characterization of the Impact of Density Gradient Centrifugation on the Profile of the Pig Sperm Transcriptome by RNA-Seq.

RNA-Seq data from human semen suggests that the study of the sperm transcriptome requires the previous elimination from the ejaculates of somatic cells carrying a larger load of RNA. Semen purification is also carried to study the sperm transcriptome in other species including swine and it is often done by density gradient centrifugation to obtain viable spermatozoa from fresh ejaculates or artificial insemination doses, thereby limiting the throughput and remoteness of the samples that can be processed in one study. The aim of this work was to evaluate the impact of purification with density gradient centrifugation by BoviPureTM on porcine sperm. Four boar ejaculates were purified with BoviPureTM and their transcriptome sequenced by RNA-Seq was compared with the RNA-Seq profiles of their paired non-purified sample. Seven thousand five hundred and nineteen protein coding genes were identified. Correlation, cluster, and principal component analysis indicated high-although not complete-similarity between the purified and the paired non-purified ejaculates. 372 genes displayed differentially abundant RNA levels between treatments. Most of these genes had lower abundances after purification and were mostly related to translation, transcription and metabolic processes. We detected a significant change in the proportion of genes of epididymal origin within the differentially abundant genes (1.3%) when compared with the catalog of unaltered genes (0.2%). In contrast, the proportion of testis-specific genes was higher in the group of unaltered genes (4%) when compared to the list of differentially abundant genes (0%). No proportion differences were identified for prostate, white blood, lymph node, tonsil, duodenum, skeletal muscle, liver, and mammary gland. Altogether, these results suggest that the purification impacts on the RNA levels of a small number of genes which are most likely caused by the removal of epididymal epithelial cells but also premature germinal cells, immature or abnormal spermatozoa or seminal exosomes with a distinct load of RNAs.

Whole genome sequencing identifies allelic ratio distortion in sperm involving genes related to spermatogenesis in a swine model.

Transmission Ratio Distortion (TRD), the uneven transmission of an allele from a parent to its offspring, can be caused by allelic differences affecting gametogenesis, fertilization or embryogenesis. However, TRD remains vaguely studied at a genomic scale. We sequenced the diploid and haploid genomes of three boars from leukocytes and spermatozoa at 50x to shed light into the genetic basis of spermatogenesis-caused Allelic Ratio Distortion (ARD). We first developed a Binomial model to identify ARD by simultaneously analysing all three males. This led to the identification of 55 ARD SNPs, most of which were animal-specific. We then evaluated ARD individually within each pig by a Fisher's exact test and identified two shared genes (TOP3A and UNC5B) and four shared genomic regions harbouring distinct ARD SNPs in the three boars. The shared genomic regions contained candidate genes with functions related to spermatogenesis including AK7, ARID4B, BDKRB2, GSK3B, NID1, NSMCE1, PALB2, VRK1 and ZC3H13. Using the Fisher's test, we also identified 378 genes containing variants with protein damaging potential in at least one boar, a high proportion of which, including FAM120B, TDRD15, JAM2 or AOX4 among others, are associated to spermatogenesis. Overall, our results show that sperm is subjected to ARD with variants associated to a wide variety of genes involved in different stages of spermatogenesis.

Identification of circular RNAs in porcine sperm and evaluation of their relation to sperm motility.

Circular RNAs (circRNAs) are emerging as a novel class of noncoding RNAs which potential role as gene regulators is quickly gaining interest. circRNAs have been studied in different tissues and cell types across several animal species. However, a thorough characterization of the circRNAome in ejaculated sperm remains unexplored. In this study, we profiled the sperm circRNA catalogue using 40 porcine ejaculates. A complex population of 1,598 circRNAs was shared in at least 30 of the 40 samples. Generally speaking, the predicted circRNAs presented low abundances and were tissue-specific. Around 80% of the circRNAs identified in the boar sperm were reported as novel. Results from abundance correlation between circRNAs and miRNAs together with the prediction of microRNA (miRNA) target sites in circRNAs suggested that circRNAs may act as miRNA sponges. Moreover, we found significant correlations between the abundance of 148 exonic circRNAs and sperm motility parameters. Two of these correlations, involving ssc_circ_1458 and ssc_circ_1321, were confirmed by RT-qPCR using 36 additional samples with extreme and opposite sperm motility values. Our study provides a thorough characterization of circRNAs in sperm and suggests that circRNAs hold potential as noninvasive biomarkers for sperm quality and male fertility.

A pilot RNA-seq study in 40 pietrain ejaculates to characterize the porcine sperm microbiome.

The microbiome plays a key role in homeostasis and health and it has been also linked to fertility and semen quality in several animal species including swine. Despite the more than likely importance of sperm bacteria on the boar's reproductive ability and the dissemination of pathogens and antimicrobial resistance genes, the high throughput characterization of the swine sperm microbiome remains scarce. We carried RNA-seq on 40 ejaculates each from a different Pietrain boar and found that a proportion of the sequencing reads did not map to the Sus scrofa genome. The current study aimed at using these reads not belonging to pig to carry a pilot study to profile the boar sperm bacterial population and its relation with 7 semen quality traits. We found that the boar sperm contains a broad population of bacteria. The most abundant phyla were Proteobacteria (39.1%), Firmicutes (27.5%), Actinobacteria (14.9%) and Bacteroidetes (5.7%). The predominant species contaminated sperm after ejaculation from soil, faeces and water sources (Bacillus megaterium, Brachybacterium faecium, Bacillus coagulans). Some potential pathogens were also found but at relatively low levels (Escherichia coli, Clostridioides difficile, Clostridium perfringens, Clostridium botulinum and Mycobacterium tuberculosis). We also identified 3 potential antibiotic resistant genes from E. coli against chloramphenicol, Neisseria meningitidis against spectinomycin and Staphylococcus aureus against linezolid. None of these genes were highly abundant. Finally, we classified the ejaculates into categories according to their bacterial features and semen quality parameters and identified two categories that significantly differed for 5 semen quality traits and 13 bacterial features including the genera Acinetobacter, Stenotrophomonas and Rhodobacter. Our results show that boar semen contains a bacterial community, including potential pathogens and putative antibiotic resistance genes, and that these bacteria may affect its reproductive performance.

A RNA-Seq Analysis to Describe the Boar Sperm Transcriptome and Its Seasonal Changes.

Understanding the molecular basis of cell function and ultimate phenotypes is crucial for the development of biological markers. With this aim, several RNA-seq studies have been devoted to the characterization of the transcriptome of ejaculated spermatozoa in relation to sperm quality and fertility. Semen quality follows a seasonal pattern and decays in the summer months in several animal species. The aim of this study was to deeply profile the transcriptome of the boar sperm and to evaluate its seasonal changes. We sequenced the total and the short fractions of the sperm RNA from 10 Pietrain boars, 5 collected in summer and 5 five sampled in winter, and identified a complex and rich transcriptome with 4,436 coding genes of moderate to high abundance. Transcript fragmentation was high but less obvious in genes related to spermatogenesis, chromatin compaction and fertility. Short non-coding RNAs mostly included piwi-interacting RNAs, transfer RNAs and microRNAs. We also compared the transcriptome of the summer and the winter ejaculates and identified 34 coding genes and 7 microRNAs with a significantly distinct distribution. These genes were mostly related to oxidative stress, DNA damage and autophagy. This is the deepest characterization of the boar sperm transcriptome and the first study linking the transcriptome and the seasonal variability of semen quality in animals. The annotation described here can be used as a reference for the identification of markers of sperm quality in pigs.

A technical assessment of the porcine ejaculated spermatozoa for a sperm-specific RNA-seq analysis.

The study of the boar sperm transcriptome by RNA-seq can provide relevant information on sperm quality and fertility and might contribute to animal breeding strategies. However, the analysis of the spermatozoa RNA is challenging as these cells harbor very low amounts of highly fragmented RNA, and the ejaculates also contain other cell types with larger amounts of non-fragmented RNA. Here, we describe a strategy for a successful boar sperm purification, RNA extraction and RNA-seq library preparation. Using these approaches our objectives were: (i) to evaluate the sperm recovery rate (SRR) after boar spermatozoa purification by density centrifugation using the non-porcine-specific commercial reagent BoviPureTM; (ii) to assess the correlation between SRR and sperm quality characteristics; (iii) to evaluate the relationship between sperm cell RNA load and sperm quality traits and (iv) to compare different library preparation kits for both total RNA-seq (SMARTer Universal Low Input RNA and TruSeq RNA Library Prep kit) and small RNA-seq (NEBNext Small RNA and TailorMix miRNA Sample Prep v2) for high-throughput sequencing. Our results show that pig SRR (~22%) is lower than in other mammalian species and that it is not significantly dependent of the sperm quality parameters analyzed in our study. Moreover, no relationship between the RNA yield per sperm cell and sperm phenotypes was found. We compared a RNA-seq library preparation kit optimized for low amounts of fragmented RNA with a standard kit designed for high amount and quality of input RNA and found that for sperm, a protocol designed to work on low-quality RNA is essential. We also compared two small RNA-seq kits and did not find substantial differences in their performance. We propose the methodological workflow described for the RNA-seq screening of the boar spermatozoa transcriptome. ABBREVIATIONS: FPKM: fragments per kilobase of transcript per million mapped reads; KRT1: keratin 1; miRNA: micro-RNA; miscRNA: miscellaneous RNA; Mt rRNA: mitochondrial ribosomal RNA; Mt tRNA: mitochondrial transference RNA; OAZ3: ornithine decarboxylase antizyme 3; ORT: osmotic resistance test; piRNA: Piwi-interacting RNA; PRM1: protamine 1; PTPRC: protein tyrosine phosphatase receptor type C; rRNA: ribosomal RNA; snoRNA: small nucleolar RNA; snRNA: small nuclear RNA; SRR: sperm recovery rate; tRNA: transfer RNA.