Oligo-dT selected spermatozoal transcript profiles differ among higher and lower fertility dairy sires

RNA Sequencing of Sperm from Healthy Cattle and Horses Reveals the Presence of a Large Bacterial Population

RNA molecules within ejaculated sperm can be characterized through whole-transcriptome sequencing, enabling the identification of pivotal transcripts that may influence reproductive success. However, the profiling of sperm transcriptomes through next-generation sequencing has several limitations impairing the identification of functional transcripts. In this study, we explored the nature of the RNA sequences present in the sperm transcriptome of two livestock species, cattle and horses, using RNA sequencing (RNA-seq) technology. Through processing of transcriptomic data derived from bovine and equine sperm cell preparations, low mapping rates to the reference genomes were observed, mainly attributed to the presence of ribosomal RNA and bacteria in sperm samples, which led to a reduced sequencing depth of RNAs of interest. To explore the presence of bacteria, we aligned the unmapped reads to a complete database of bacterial genomes and identified bacteria-associated transcripts which were characterized. This analysis examines the limitations associated with sperm transcriptome profiling by reporting the nature of the RNA sequences among which bacterial RNA was found. These findings can aid researchers in understanding spermatozoal RNA-seq data and pave the way for the identification of molecular markers of sperm performance.

Bull Sperm SWATH-MS-Based Proteomics Reveals Link between High Fertility and Energy Production, Motility Structures, and Sperm-Oocyte Interaction

The prediction of male or semen fertility potential remains a persistent challenge that has yet to be fully resolved. This work analyzed several in vitro parameters and proteome of spermatozoa in bulls cataloged as high- (HF; n = 5) and low-field (LF; n = 5) fertility after more than a thousand artificial inseminations. Sperm motility was evaluated by computer-assisted sperm analysis. Sperm viability, mitochondrial membrane potential (MMP) and reactive oxygen species (mROS) of spermatozoa were assessed by flow cytometry. Proteome was evaluated by the SWATH-MS procedure. Spermatozoa of HF bulls showed significantly higher total motility than the LF group (41.4% vs 29.7%). Rates of healthy sperm (live, high MMP, and low mROS) for HF and LF bull groups were 49% and 43%, respectively (p > 0.05). Spermatozoa of HF bulls showed a higher presence of differentially abundant proteins (DAPs) related to both energy production (COX7C), mainly the OXPHOS pathway, and the development of structures linked with the motility process (TPPP2, SSMEM1, and SPAG16). Furthermore, we observed that equatorin (EQTN), together with other DAPs related to the interaction with the oocyte, was overrepresented in HF bull spermatozoa. The biological processes related to protein processing, catabolism, and protein folding were found to be overrepresented in LF bull sperm in which the HSP90AA1 chaperone was identified as the most DAP. Data are available via ProteomeXchange with identifier PXD042286.

Effect of arginine-induced motility and capacitation on RNA population in goat spermatozoa

INTRODUCTION

In vitro capacitation is essential in assisted reproductive technologies (ART) for embryo production. Recently, arginine has been proven to enhance capacitation in mammalian spermatozoa. However, the detailed mechanism of action of arginine remains elusive.

AIM

This study investigated the effect of arginine-induced capacitation and motility enhancement on the spermatozoal RNA (spRNA) population in goats.

MATERIAL AND METHODS

Goat spermatozoa were treated with arginine for up to six hours and compared with non-treated or PHE (penicillamine, hypotaurine, and epinephrine)-treated spermatozoa at different intervals (0, 1, 2, 4, and 6 hours). Sperm parameters, including viability, individual motility, capacitation, acrosome reaction, and ROS production, were evaluated. The spRNA population was analyzed by short-read RNA sequencing (RNA-seq).

RESULTS

The percentage of capacitated (73.21 ± 4.22%) and acrosome reacted (18.35 ± 0.56%) spermatozoa was highest in arginine treatment, while PHE treatment showed the highest percentage (79.82 ± 4.31%) of motile spermatozoa from 0 to 4 hours of incubation. RNA-seq analysis identified 1,321 differentially expressed genes (DEGs) in arginine-treated spermatozoa compared to the control. The PGK2, RNASE10, ODF1, and ROPN1L genes involved in sperm motility and ACR, DKKL1, KCNJ11, and PRND genes involved in the capacitation process were upregulated in arginine-treated spermatozoa. The DEGs regulate sperm capacitation-related cAMP-PKA, PI3-Akt, calcium, and MAPK signaling pathways.

CONCLUSION

The arginine-induced capacitation and enhanced sperm motility were associated with the upregulation of several genes involved in sperm motility and capacitation pathways. The comparative study also suggests that arginine may be used in lieu of PHE for motility enhancement and in vitro capacitation of goat spermatozoa.

Integration of omics studies indicates that species-dependent molecular mechanisms govern male fertility

BACKGROUND

Comparative and comprehensive omics studies have recently been conducted to provide a comprehensive understanding of the biological mechanisms underlying infertility. However, because these huge omics datasets often contain irrelevant information, editing strategies for summarizing and filtering the data are necessary prerequisite steps for identifying biomarkers of male fertility. Here, we attempted to integrate omics data from spermatozoa with normal and below-normal fertility from boars and bulls, including transcriptomic, proteomic, and metabolomic data. Pathway enrichment analysis was conducted and visualized using g:Profiler, Cytoscape, EnrichmentMap, and AutoAnnotation to determine fertility-related biological functions according to species.

RESULTS

In particular, gamete production and protein biogenesis-associated pathways were enriched in bull spermatozoa with below-normal fertility, whereas mitochondrial-associated metabolic pathways were enriched in boar spermatozoa with normal fertility. These results indicate that below-normal fertility may be determined by aberrant regulation of protein synthesis during spermatogenesis, and the modulation of reactive oxygen species generation to maintain capacitation and the acrosome reaction governs boar sperm fertility.

CONCLUSION

Overall, this approach demonstrated that distinct molecular pathways drive sperm fertility in mammals in a species-dependent manner. Moreover, we anticipate that searching for species-specific signaling pathways may aid in the discovery of fertility-related biomarkers within large omics datasets.

Testicular hypoplasia in Nellore Cattle: Genetic analysis and functional analysis of genome-wide association study results

Characterized by the incomplete development of the germinal epithelium of the seminiferous tubules, Testicular hypoplasia (TH) leads to decreased sperm concentration, increased morphological changes in sperm and azoospermia. Economic losses resulting from the disposal of affected bulls reduce the efficiency of meat production systems. A genome-wide association study and functional analysis were performed to identify genomic windows and the underlying positional candidate genes associated with TH in Nellore cattle. Phenotypic and pedigree data from 207,195 animals and genotypes (461,057 single nucleotide polymorphism, SNP) from 17,326 sires were used in this study. TH was evaluated as a binary trait measured at 18 months of age. A possible correlated response on TH resulting from the selection for scrotal circumference was evaluated by using a two-trait analysis. Thus, estimated breeding values were calculated by fitting a linear-threshold animal model in a Bayesian approach. The SNP effects were estimated using the weighted single-step genomic BLUP method. Twelve non-overlapping windows of 20 adjacent SNP that explained more than 1% of the additive genetic variance were selected for candidate gene annotation. Functional and gene prioritization analysis of the candidate genes identified genes (KHDRBS3, GPX5, STAR, ERLIN2), which might play an important role in the expression of TH due to their known roles in the spermatogenesis process, synthesis of steroids and lipid metabolism.

Molecular Markers: A New Paradigm in the Prediction of Sperm Freezability

For decades now, sperm cryopreservation has been a pillar of assisted reproduction in animals as well as humans. Nevertheless, the success of cryopreservation varies across species, seasons, and latitudes and even within the same individual. With the dawn of progressive analytical techniques in the field of genomics, proteomics, and metabolomics, new options for a more accurate semen quality assessment have become available. This review summarizes currently available information on specific molecular characteristics of spermatozoa that could predict their cryotolerance before the freezing process. Understanding the changes in sperm biology as a result of their exposure to low temperatures may contribute to the development and implementation of appropriate measures to assure high post-thaw sperm quality. Furthermore, an early prediction of cryotolerance or cryosensitivity may lead to the establishment of customized protocols interconnecting adequate sperm processing procedures, freezing techniques, and cryosupplements that are most feasible for the individual needs of the ejaculate.

Sperm Transcriptome Analysis Accurately Reveals Male Fertility Potential in Livestock

Nowadays, selection of superior male candidates in livestock as a source of frozen semen based on sperm quality at the cellular level is not considered accurate enough for predicting the potential of male fertility. Sperm transcriptome analysis approaches, such as messenger RNA levels, have been shown to correlate with fertility rates. Using this technology in livestock growth has become the principal method, which can be widely applied to predict male fertility potential in the livestock industry through the analysis of the sperm transcriptome. It provides the gene expression to validate the function of sperm in spermatogenesis, fertilization, and embryo development, as the parameters of male fertility. This review proposes a transcriptomic analysis approach as a high-throughput method to predict the fertility potential of livestock more accurately in the future.

Identification of differentially expressed mRNAs and miRNAs in spermatozoa of bulls of varying fertility

Bulls used in artificial insemination, with apparently normal semen quality, can vary significantly in their field fertility. This study aimed to characterize the transcriptome of spermatozoa from high (HF) and low (LF) fertility bulls at the mRNA and miRNA level in order to identify potential novel markers of fertility. Holstein-Friesian bulls were assigned to either the HF or LF group (n = 10 per group) based on an adjusted national fertility index from a minimum of 500 inseminations. Total RNA was extracted from a pool of frozen-thawed spermatozoa from three different ejaculates per bull, following which mRNA-seq and miRNA-seq were performed. Six mRNAs and 13 miRNAs were found differentially expressed (P < 0.05, FC > 1.5) between HF and LF bulls. Of particular interest, the gene pathways targeted by the 13 differentially expressed miRNAs were related to embryonic development and gene expression regulation. Previous studies reported that disruptions to protamine 1 mRNA (PRM1) had deleterious consequences for sperm chromatin structure and fertilizing ability. Notably, PRM1 exhibited a higher expression in spermatozoa from LF than HF bulls. In contrast, Western Blot analysis revealed a decrease in PRM1 protein abundance for spermatozoa from LF bulls; this was not associated with increased protamine deficiency (measured by the degree of chromatin compaction) or DNA fragmentation, as assessed by flow cytometry analyses. However, protamine deficiency was positively and moderately correlated with the percentage of spermatozoa with DNA fragmentation, irrespective of fertility group. This study has identified potential biomarkers that could be used for improving semen quality assessments of bull fertility.

Expression of TXNRD1, HSPA4L and ATP1B1 Genes Associated with the Freezability of Boar Sperm

Cryopreservation is associated with increased oxidative stress, which is responsible for sperm damage. We analyzed the effect of cryopreservation on mRNA and protein expression of thioredoxin reductase 1 (TXNRD1), heat shock protein family A (HSP 70) member 4 like (HSPA4L) and sodium/potassium-transporting ATPase subunit beta-1 (ATP1B1) genes in boar sperm with different freezability. Boars were classified as having good and poor semen freezability (GSF and PSF, respectively), according to the assessment of post-thaw sperm motility. Total RNA was isolated from fresh pre-freeze (PF) and frozen-thawed (FT) sperm from five boars of the GSF and PSF groups, respectively. Quantification of TXNRD1, HSPA4L and ATP1B1 gene expression was performed by RT-qPCR analysis. Proteins extracted from sperm were subjected to Western blotting and SDS-PAGE analyses. Poor freezability ejaculates were characterized by significantly higher relative mRNA expression levels of TXNRD1 and HSPA4L in FT sperm compared with the fresh PF sperm. Furthermore, the relative mRNA expression level of ATP1B1 was significantly higher in the fresh PF sperm of the GSF group. Western blotting analysis revealed significantly higher relative expression of TXNRD1 protein in the fresh PF sperm of the GSF group, while HSPA4L protein expression was markedly increased in FT sperm of the PSF group. Electrophoretic and densitometric analyses revealed a higher number of proteins in the fresh PF and FT sperm of the PSF and GSF groups, respectively. The results of this study indicate that ATP1B1 mRNA expression in the fresh PF sperm is a promising cryotolerance marker, while the variations of TXNRD1 and HSPA4L protein expression in the fresh PF or FT sperm provide useful information that may help to elucidate their biological significance in cryo-damage.

Sperm Transcripts Associated With Odorant Binding and Olfactory Transduction Pathways Are Altered in Breeding Bulls Producing Poor-Quality Semen

Spermatozoa carries a reservoir of mRNAs regulating sperm functions and fertilizing potential. Although it is well recognized that a considerable proportion of high genetic merit breeding bulls produce poor-quality semen, the transcriptomic alterations in spermatozoa from such bulls are not understood. In the present study, comparative high-throughput transcriptomic profiling of spermatozoa from good and poor-quality semen-producing bulls was carried out to identify the transcripts associated with semen quality. Using next-generation sequencing (NGS), we identified 11,632 transcripts in Holstein Friesian bull spermatozoa; after total hit normalization, a total of 544 transcripts were detected, of which 185 transcripts were common to both good and poor-quality semen, while 181 sperm transcripts were unique to good quality semen, and 178 transcripts were unique to poor-quality semen. Among the co-expressed transcripts, 31 were upregulated, while 108 were downregulated, and 46 were neutrally expressed in poor-quality semen. Bioinformatics analysis revealed that the dysregulated transcripts were predominantly involved in molecular function, such as olfactory receptor activity and odor binding, and in biological process, such as detection of chemical stimulus involved in sensory perception, sensory perception of smell, signal transduction, and signal synaptic transmission. Since a majority of the dysregulated transcripts were involved in the olfactory pathway (85% of enriched dysregulated genes were involved in this pathway), the expression of selected five transcripts associated with this pathway (OR2T11, OR10S1, ORIL3, OR5M11, and PRRX1) were validated using real-time qPCR, and it was found that their transcriptional abundance followed the same trend as observed in NGS; the sperm transcriptional abundance of OR2T11 and OR10S1 differed significantly (p < 0.05) between good and poor-quality semen. It is concluded that poor-quality semen showed altered expression of transcripts associated with olfactory receptors and pathways indicating the relationship between olfactory pathway and semen quality in bulls.

Significance and Relevance of Spermatozoal RNAs to Male Fertility in Livestock

Livestock production contributes to a significant part of the economy in developing countries. Although artificial insemination techniques brought substantial improvements in reproductive efficiency, male infertility remains a leading challenge in livestock. Current strategies for the diagnosis of male infertility largely depend on the evaluation of semen parameters and fail to diagnose idiopathic infertility in most cases. Recent evidences show that spermatozoa contains a suit of RNA population whose profile differs between fertile and infertile males. Studies have also demonstrated the crucial roles of spermatozoal RNA (spRNA) in spermatogenesis, fertilization, and early embryonic development. Thus, the spRNA profile may serve as unique molecular signatures of fertile sperm and may play pivotal roles in the diagnosis and treatment of male fertility. This manuscript provides an update on various spRNA populations, including protein-coding and non-coding RNAs, in livestock species and their potential role in semen quality, particularly sperm motility, freezability, and fertility. The contribution of seminal plasma to the spRNA population is also discussed. Furthermore, we discussed the significance of rare non-coding RNAs (ncRNAs) such as long ncRNAs (lncRNAs) and circular RNAs (circRNAs) in spermatogenic events.

Cellular and Molecular Insights Into the Etiology of Subfertility/Infertility in Crossbred Bulls (Bos taurus × Bos indicus): A Review

Crossbreeding of indigenous cattle (Bos indicus) with improved (Bos taurus) breeds gained momentum and economic relevance in several countries to increase milk production. While production performance of the crossbred offspring is high due to hybrid vigor, they suffer from a high incidence of reproductive problems. Specifically, the crossbred males suffer from serious forms of subfertility/infertility, which can have a significant effect because semen from a single male is used to breed several thousand females. During the last two decades, attempts have been made to understand the probable reasons for infertility in crossbred bulls. Published evidence indicates that testicular cytology indices, hormonal concentrations, sperm phenotypic characteristics and seminal plasma composition were altered in crossbred compared to purebred males. A few recent studies compared crossbred bull semen with purebred bull semen using genomics, transcriptomics, proteomics and metabolomics; molecules potentially associated with subfertility/infertility in crossbred bulls were identified. Nevertheless, the precise reason behind the poor quality of semen and high incidence of sub-fertility/infertility in crossbred bulls are not yet well defined. To identify the underlying etiology for infertility in crossbred bulls, a thorough understanding of the magnitude of the problem and an overview of the prior art is needed; however, such systematically reviewed information is not available. Therefore, the primary focus of this review is to compile and analyze earlier findings on crossbred bull fertility/infertility. In addition, the differences between purebred and crossbred males in terms of testicular composition, sperm phenotypic characteristics, molecular composition, environmental influence and other details are described; future prospects for research on crossbred males are also outlined.

Comparative Transcriptomic Analysis of Spermatozoa From High- and Low-Fertile Crossbred Bulls: Implications for Fertility Prediction

Crossbred bulls produced by crossing Bos taurus and Bos indicus suffer with high incidence of infertility/subfertility problems; however, the etiology remains poorly understood. The uncertain predictability and the inability of semen evaluation techniques to maintain constant correlation with fertility demand for alternate methods for bull fertility prediction. Therefore, in this study, the global differential gene expression between high- and low-fertile crossbred bull sperm was assessed using a high-throughput RNA sequencing technique with the aim to identify transcripts associated with crossbred bull fertility. Crossbred bull sperm contained transcripts for 13,563 genes, in which 2,093 were unique to high-fertile and 5,454 were unique to low-fertile bulls. After normalization of data, a total of 776 transcripts were detected, in which 84 and 168 transcripts were unique to high-fertile and low-fertile bulls, respectively. A total of 176 transcripts were upregulated (fold change > 1) and 209 were downregulated (<1) in low-fertile bulls. Gene ontology analysis identified that the sperm transcripts involved in the oxidative phosphorylation pathway and biological process such as multicellular organism development, spermatogenesis, and in utero embryonic development were downregulated in low-fertile crossbred bull sperm. Sperm transcripts upregulated and unique to low-fertile bulls were majorly involved in translation (biological process) and ribosomal pathway. With the use of RT-qPCR, selected sperm transcripts (n = 12) were validated in crossbred bulls (n = 12) with different fertility ratings and found that the transcriptional abundance of ZNF706, CRISP2, TNP2, and TNP1 genes was significantly (p < 0.05) lower in low-fertile bulls than high-fertile bulls and was positively (p < 0.05) correlated with conception rate. It is inferred that impaired oxidative phosphorylation could be the predominant reason for low fertility in crossbred bulls and that transcriptional abundance of ZNF706, CRISP2, TNP2, and TNP1 genes could serve as potential biomarkers for fertility in crossbred bulls.

Gene promoter polymorphisms in boar spermatozoa differing in freezability

Single nucleotide polymorphisms (SNPs) in the 5'-flanking regulatory regions of genes could affect their expression levels. This is a follow-up study aimed to identify polymorphic variants in the 5'-flanking regulatory regions of genes expressed in boar spermatozoa, and to predict the interactions of such variants with transcription factors (TFs) on the gene promoter activity, using bioinformatics. Five and six boars were classified as having good and poor semen freezability (GSF and PSF, respectively) according to post-thaw (PT) assessment of sperm motility and membrane integrity characteristics. The 5'-flanking region sequences of the 14 genes (FOS, NFATC3, EAF2, FGF-14, BAMBI, RAB33B, CKS2, LARS2, SLC25A16, ACADM, CPT2, CCT3, DTD2 and CCDC85A) were PCR amplified and analyzed by Sanger sequencing method. A total of 32 polymorphic variants were identified in the 5'-flanking regions of the genes, including 4 insertion/deletion (indel) polymorphisms, and 8 unknown (novel) SNPs. Multiple sequence alignment analysis revealed a 26-bp indel variant in the 5'-flanking region of the LARS2 gene, which showed greater protein expression in spermatozoa from boars of the PSF group. It was found that 17 polymorphic variants, observed in the differentially expressed (DE) genes, showed significant allele frequency differences between the GSF and PSF groups. Polymorphic variants in the 5'-flanking regulatory regions of the genes contributed to the decrease or increase in the binding affinity for different testis-specific TFs, such as SMAD1, NF-1, FOXMI, RXRA, STAT4 and C/EBPβ. This study provides more insights into the mechanisms responsible for variations in transcriptional activity in promoters of genes expressed in boar spermatozoa. The allelic variants are promising genetic markers for predicting the freezability of boar spermatozoa.

Single-molecule long-read sequencing reveals a conserved intact long RNA profile in sperm

Sperm contributes diverse RNAs to the zygote. While sperm small RNAs have been shown to impact offspring phenotypes, our knowledge of the sperm transcriptome, especially the composition of long RNAs, has been limited by the lack of sensitive, high-throughput experimental techniques that can distinguish intact RNAs from fragmented RNAs, known to abound in sperm. Here, we integrate single-molecule long-read sequencing with short-read sequencing to detect sperm intact RNAs (spiRNAs). We identify 3440 spiRNA species in mice and 4100 in humans. The spiRNA profile consists of both mRNAs and long non-coding RNAs, is evolutionarily conserved between mice and humans, and displays an enrichment in mRNAs encoding for ribosome. In sum, we characterize the landscape of intact long RNAs in sperm, paving the way for future studies on their biogenesis and functions. Our experimental and bioinformatics approaches can be applied to other tissues and organisms to detect intact transcripts.

Impact of Cryopreservation on Spermatozoa Freeze-Thawed Traits and Relevance OMICS to Assess Sperm Cryo-Tolerance in Farm Animals

Sperm cryopreservation is a powerful tool for the livestock breeding program. Several technical attempts have been made to enhance the efficiency of spermatozoa cryopreservation in different farm animal species. However, it is well-recognized that mammalian spermatozoa are susceptible to cryo-injury caused by cryopreservation processes. Moreover, the factors leading to cryo-injuries are complicated, and the cryo-damage mechanism has not been methodically explained until now, which directly influences the quality of frozen–thawed spermatozoa. Currently, the various OMICS technologies in sperm cryo-biology have been conducted, particularly proteomics and transcriptomics studies. It has contributed while exploring the molecular alterations caused by cryopreservation, identification of various freezability markers and specific proteins that could be added to semen diluents before cryopreservation to improve sperm cryo-survival. Therefore, understanding the cryo-injury mechanism of spermatozoa is essential for the optimization of current cryopreservation processes. Recently, the application of newly-emerged proteomics and transcriptomics technologies to study the effects of cryopreservation on sperm is becoming a hotspot. This review detailed an updated overview of OMICS elements involved in sperm cryo-tolerance and freeze-thawed quality. While also detailed a mechanism of sperm cryo-injury and utilizing OMICS technology that assesses the sperm freezability potential biomarkers as well as the accurate classification between the excellent and poor freezer breeding candidate.

Transcriptomic Profiling of Buffalo Spermatozoa Reveals Dysregulation of Functionally Relevant mRNAs in Low-Fertile Bulls

Although, it is known that spermatozoa harbor a variety of RNAs that may influence embryonic development, little is understood about sperm transcriptomic differences in relation to fertility, especially in buffaloes. In the present study, we compared the differences in sperm functional attributes and transcriptomic profile between high- and low-fertile buffalo bulls. Sperm membrane and acrosomal integrity were lower (P < 0.05), while protamine deficiency and lipid peroxidation were higher (P < 0.05) in low- compared to high-fertile bulls. Transcriptomic analysis using mRNA microarray technology detected a total of 51,282 transcripts in buffalo spermatozoa, of which 4,050 transcripts were differentially expressed, and 709 transcripts were found to be significantly dysregulated (P < 0.05 and fold change >1) between high- and low-fertile bulls. Majority of the dysregulated transcripts were related to binding activity, transcription, translation, and metabolic processes with primary localization in the cell nucleus, nucleoplasm, and in cytosol. Pathways related to MAPK signaling, ribosome pathway, and oxidative phosphorylation were dysregulated in low-fertile bull spermatozoa. Using bioinformatics analysis, we observed that several genes related to sperm functional attributes were significantly downregulated in low-fertile bull spermatozoa. Validation of the results of microarray analysis was carried out using real-time qPCR expression analysis of selected genes (YBX1, ORAI3, and TFAP2C). The relative expression of these genes followed the same trend in both the techniques. Collectively, this is the first study to report the transcriptomic profile of buffalo spermatozoa and to demonstrate the dysregulation of functionally relevant transcripts in low-fertile bull spermatozoa. The results of the present study open up new avenues for understanding the etiology for poor fertility in buffalo bulls and to identify fertility biomarkers.

Comparison of spermatozoal RNA extraction methods in goats

RNA sequencing (RNAseq) has divulged newer role of spermatozoal RNA in male fertility. This study aimed to evaluate different sperm purification and RNA extraction methods for long-read RNA sequencing of poly(A) transcriptome in goat spermatozoa. Sperm samples were purified by swim-up separation using different purification medium. Spermatozoal RNA was extracted by seven different methods with additional supplementation of reducing agents in lysis buffer. poly(A) selected RNA was used for cDNA library preparation and long-read RNAseq in Nanopore sequencer. Sperm purification by 1 h swim-up resulted in higher recovery (89.20 ± 1.15%), motility (82.33 ± 1.53%), viability (88.10 ± 5.03%) and plasma membrane integrity (71.33 ± 4.51%) in sperm TALP (sp-TL) medium. A monophasic solution of GITC with phenol and DTT resulted in the highest yield of large sized RNAs (3.89 ± 0.46 ng/million cells) suitable for long-read RNAseq of poly(A) transcripts. RNAseq resulted in reads of length, ranging from 500bp to 2 Kb. A total of 123 transcripts were identified in goat spermatozoa by de novo assembly and included sperm-specific transcripts such as CATSPERG, PRM2, CYLC2, SPATA6, PLCZ1 etc. This study is the first report of long-read RNAseq of poly(A) transcriptome in goat spermatozoa.

RNA-Seq analysis reveals functionally relevant coding and non-coding RNAs in crossbred bull spermatozoa

•

RNA-Seq analysis was done to characterize the transcriptome of crossbred bull spermatozoa.

•

Among the 13,814 transcripts detected, 431 had FPKM > 1 and 13,673 had FPKM > 0 or < 1.

•

Coding and non-coding RNAs account for 13,145 (95.15%) and 152 (1.1%), respectively.

•

Sperm transcripts were mainly related to ribosome, oxidative phosphorylation and spliceosome pathways.

•

qPCR analysis showed individual variations in transcriptional abundance of selected genes.

Sperm, which are believed to be transcriptionally and translationally inactive, synthesize RNA and proteins before there is gradual disappearance of the ribosome during chromatin compaction. Sperm transfer several functionally relevant transcripts to the oocyte, controlling maternal-zygotic transition and embryonic development. The present study was undertaken to profile and analyze sperm transcripts comprehensively using Next Generation Ribonucleic acid sequencing technology in Holstein Friesian x Tharparkar crossbred bulls. The results from global transcriptomic profiling revealed transcripts for 13,814 genes; of which 431 transcripts were expressed with >1 FPKM and 13,383 transcripts were expressed with >0 or <1 FPKM. The abundant mRNA transcripts of crossbred bull sperm were PRM1 and HMGB4. Gene ontology of transcripts with>1 FPKM revealed there was a major involvement in the structural constituent of ribosomes and translation. Results from pathway enrichment indicated the connection between ribosome, oxidative phosphorylation and spliceosome pathways and the transcripts of crossbred bull spermatozoa. The transcriptional abundance of selected genes, validated using RT-qPCR, indicated significant variations between bulls. Collectively, it may be inferred that the transcripts in crossbred bull sperm were heavily implicated in functions such as the structural constituent of ribosomes and translation, and pathways such as ribosome, oxidative phosphorylation and spliceosome. Further studies using larger sample sizes are required to understand the possible implications of transcriptomic variations on semen quality and fertility.

Dynamics of paternal contributions to early embryo development in large animals

This review focuses on current knowledge of paternal contributions to preimplantation embryonic development with particular emphasis on large animals. Specifically, the included content aims to summarize genomic and epigenomic contributions of paternally expressed genes, their regulation, and chromatin structure that are indispensable for early embryo development. The accumulation of current knowledge will summarize conserved allelic function among species to include functional molecular and genomic studies across large domestic animals in context with reference to founding experimental models.

Cryopreservation and its effects on motility and gene expression patterns and fertilizing potential of bovine epididymal sperm

Despite encountering new challenges in using epididymal sperm recovered from cauda epididymides, this accessible and, in some species, worthwhile sample makes inevitable the further development of a suitable cryopreservation protocol. In this study, sperm was recovered from the epididymis of 4°C overnight stored slaughtered bulls' testes and the effects of cryopreservation on the bovine epididymal sperm motility (with CASA) and gene expression patterns (with quantitative Real time-PCR) were evaluated. Moreover the fertilizing potential of cryopreserved epididymal sperm was used in in vitro fertilization (IVF). After freezing and thawing of epididymal sperm, total and slow progressive sperm motility, VCL, VAP, MAD, ALH and BCF were significantly decreased (p < .05), while in the parameters of fast progressive motility, VSL, LIN, WOB and STR there were not any significant variations in the frozen sperm compared to fresh (non-frozen) counterpart. The assessment of abundance of transcripts encoding motility (TSSK6) and fertility (PRM1 and PRM2)-related genes in epididymal sperm, showed that these transcripts were affected by freezing especially in slow progressive motility status (p < .01). Furthermore, cleavage and blastocyst rate did not present any significant differences between bovine embryos produced in vitro by fresh or frozen-thawed epididymal sperm. It can be concluded that epididymal sperm has enough freezability after overnight testes storage, and cryopreservation could not affect the percentage of in vitro produced embryos in spite of the changes of relative abundance of some transcripts and direction progressive motility pattern of sperm.

Identification of markers associated with estimated breeding value and horn colour in Hungarian Grey cattle

OBJECTIVE

This study was conducted to estimate effect of single nucleotide polymorphisms (SNP) on the estimated breeding value of Hungarian Grey (HG) bulls and to find markers associated with horn colour.

METHODS

Genotypes 136 HG animals were determined on Geneseek high-density Bovine SNP 150K BeadChip. A multi-locus mixed-model was applied for statistical analyses.

RESULTS

Six SNPs were identified to be associated (-log10P>10) with green and white horn. These loci are located on chromosome 1, 3, 9, 18, and 25. Seven loci (on chromosome 1, 3, 6, 9, 10, 28) showed considerable association (-log10P>10) with the estimated breeding value.

CONCLUSION

Analysis provides markers for further research of horn colour and supplies markers to achieve more effective selection work regarding estimated breeding value of HG.

Integration of whole-genome DNA methylation data with RNA sequencing data to identify markers for bull fertility

Predicting bull fertility prior to breeding is a current challenge for the dairy industry. The use of molecular biomarkers has been previously assessed. However, the integration of this information has not been performed to extract biologically relevant markers. The goal of this study was to integrate DNA methylation data with previously published RNA-sequencing results in order to identify candidate markers for sire fertility. A total of 1765 differentially methylated cytosines were found between high- and low-fertility sires. Ten genes associated with 11 differentially methylated cytosines were found in a previous study of gene expression between high- and low-fertility sires. Additionally, two of these genes code for proteins found exclusively in bull seminal plasma. Collectively, our results reveal 10 genes that could be used in the future as a panel for predicting bull fertility.

Characterization and functional roles of paternal RNAs in 2-4 cell bovine embryos

Embryos utilize oocyte-donated RNAs until they become capable of producing RNAs through embryonic genome activation (EGA). The sperm's influence over pre-EGA RNA content of embryos remains unknown. Recent studies have revealed that sperm donate non-genomic components upon fertilization. Thus, sperm may also contribute to RNA presence in pre-EGA embryos. The first objective of this study was to investigate whether male fertility status is associated with the RNAs present in the bovine embryo prior to EGA. A total of 65 RNAs were found to be differentially expressed between 2-4 cell bovine embryos derived from high and low fertility sires. Expression patterns were confirmed for protein phosphatase 1 regulatory subunit 36 (PPP1R36) and ataxin 2 like (ATXN2L) in three new biological replicates. The knockdown of ATXN2L led to a 22.9% increase in blastocyst development. The second objective of this study was to characterize the parental origin of RNAs present in pre-EGA embryos. Results revealed 472 sperm-derived RNAs, 2575 oocyte-derived RNAs, 2675 RNAs derived from both sperm and oocytes, and 663 embryo-exclusive RNAs. This study uncovers an association of male fertility with developmentally impactful RNAs in 2-4 cell embryos. This study also provides an initial characterization of paternally-contributed RNAs to pre-EGA embryos. Furthermore, a subset of 2-4 cell embryo-specific RNAs was identified.

Transcriptome analysis of boar spermatozoa with different freezability using RNA-Seq

Semen freezability is associated with genetic markers, and there is a diverse set of sperm transcripts that have been attributed to various cellular functions. RNA-Seq was performed to compare the transcript profiles of spermatozoa from boars with different semen freezability. We examined ejaculates from the Polish large white (PLW) boars that were classified as having good and poor semen freezability (GSF and PSF, respectively; n = 3 boars per group) by assessing post-thaw motility characteristics, mitochondrial membrane potential, plasma membrane and acrosome integrity. Total RNA was isolated from fresh spermatozoa from boars of the GSF and PSF groups and subjected to RNA-Seq (Illumina NextSeq 500 platform). Transcript abundance was assessed with the DESeq2, DESeq, and EdgeR Bioconductor R packages, and varying numbers of differentially expressed gene (DEG) transcripts were detected in the spermatozoa of each boar. Using RNA-Seq, we identified several genes associated with inflammation and apoptosis (FOS, NFATC3, ITGAL, EAF2 and ZDHHC14), spermatogenesis (FGF-14 and BAMBI), autophagy (RAB33B), protein phosphorylation (PTPRU and PTPN2) and energy metabolism (ND6 and ACADM) that were predominantly up-regulated in poor freezability ejaculates. Quantitative reverse transcription polymerase chain reaction (RT-qPCR) validated the transcript expression levels detected by RNA-Seq and thus confirmed the reliability of this technique. Subsequent validation with western blotting showed that the expression of three proteins was in accordance with the transcript abundance. Overall, we demonstrated that the up-regulation of the DEG transcripts in spermatozoa was associated with poor semen freezability. We suggest that spermatozoa transcriptome profiling provides a foundation to further elucidate the relevance of sperm-related transcripts on cryo-survival. The sperm-related transcripts, namely FOS, NFATC3, EAF2, BAMBI, PTPRU, PTPN2, ND6 and ACADM, are potential markers for predicting the freezability of boar semen.

Integrated analysis of mRNAs and long noncoding RNAs in the semen from Holstein bulls with high and low sperm motility

Sperm motility is the main index used to assess the quality of bull semen. It may also be used to evaluate the fertility potential of bulls. Protein-coding mRNA and long noncoding RNA (lncRNA) participate in the regulation of spermatogenesis. Here, we employed strand-specific RNA sequencing to profile the semen transcriptome (mRNA and lncRNA) of six paired full-sibling Holstein bulls with divergent sperm motility and to determine the functions of mRNA and lncRNA in sperm motility. Among 20,875 protein-encoding genes detected in semen, 19 were differentially expressed between the high sperm motility group (H: H1, H2, and H3) and low sperm motility group (L: L1, L2, and L3). Of the 11,561 lncRNAs identified in sperm, 2,517 were differentially expressed between the H and L groups. We found that TCONS_00041733 lncRNA targets the node gene EFNA1 (ephrin A1), involved in male reproductive physiology. Our study provides a global mRNA and lncRNA transcriptome of bull semen, as well as novel insights into the regulation of neighboring protein coding by lncRNAs and the influence of mRNAs on sperm motility.

Database on spermatozoa transcriptogram of catagorised Frieswal crossbred (Holstein Friesian X Sahiwal) bulls

Bull spermatozoa contain different functional genes and many of them plays important roles in different stages of spermatogenesis, spermatozoa kinetics, fertilization as well as embryonic development. RNA deep sequencing is one of the preferred tools for absolute quantification of messenger RNA. The intention of the current study was to investigate the abundance of spermatozoal transcripts in categorized Frieswal (Holstein-Friesian X Sahiwal) crossbred bull semen through RNA deep sequencing. A total 1546561 and 1019308 numbers of reads were identified among good and poor quality bull spermatozoa based on their conception rate. Post mapping with Bos taurus reference genome identified 1,321,236 and 842,022 number of transcripts among good and poor quality RNA libraries, respectively. However, a total number of 3510 and 6759 functional transcripts were identified among good and poor quality bull spermatozoa, respectively. Most of the identified transcripts were related to spermatozoa functions, embryonic development and other functional aspects of fertilization. Wet laboratory validation of the top five selected transcripts (AKAP4, PRM1, ATP2B4, TRIM71 and SLC9B2) illustrated the significant (p < 0.01) level of expression in the good quality crossbred bull semen than the poor quality counterparts. The present study with comprehensive profiling of spermatozoal transcripts provides a useful non-invasive tool to understand the causes of as well as an effective way to predict male infertility in crossbred bulls.

Exploration of miRNA and mRNA Profiles in Fresh and Frozen-Thawed Boar Sperm by Transcriptome and Small RNA Sequencing

Due to lower farrowing rate and reduced litter size with frozen-thawed semen, over 90% of artificial insemination (AI) is conducted using liquid stored boar semen. Although substantial progress has been made towards optimizing the cryopreservation protocols for boar sperm, the influencing factors and underlying mechanisms related to cryoinjury and freeze tolerance of boar sperm remain largely unknown. In this study, we report the differential expression of mRNAs and miRNAs between fresh and frozen-thawed boar sperm using high-throughput RNA sequencing. Our results showed that 567 mRNAs and 135 miRNAs were differentially expressed (DE) in fresh and frozen-thawed boar sperm. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses revealed that the majority of DE mRNAs were enriched in environmental information processing such as cytokine-cytokine receptor interactions, PI3K-Akt signaling, cell adhesion, MAPK, and calcium signaling pathways. Moreover, the targets of DE miRNAs were enriched in significant GO terms such as cell process, protein binding, and response to stimuli. In conclusion, we speculate that DE mRNAs and miRNAs are heavily involved in boar sperm response to environment stimuli, apoptosis, and metabolic activities. The differences in expression also reflect the various structural and functional changes in sperm during cryopreservation.

Review: Genomics of bull fertility

Fertility is one of the most economically important traits in both beef and dairy cattle production; however, only female fertility is typically subjected to selection. Male and female fertility have only a small positive genetic correlation which is likely due to the existence of a relatively small number of genetic variants within each breed that cause embryonic and developmental losses. Genomic tools have been developed that allow the identification of lethal recessive loci based upon marker haplotypes. Selection against haplotypes harbouring lethal alleles in conjunction with selection to improve female fertility will result in an improvement in male fertility. Genomic selection has resulted in a two to fourfold increase in the rate of genetic improvement of most dairy traits in US Holstein cattle, including female fertility. Considering the rapidly increasing rate of adoption of high-throughput single nucleotide polymorphism genotyping in both the US dairy and beef industries, genomic selection should be the most effective of all currently available approaches to improve male fertility. However, male fertility phenotypes are not routinely recorded in natural service mating systems and when artificial insemination is used, semen doses may be titrated to lower post-thaw progressively motile sperm numbers for high-merit and high-demand bulls. Standardization of sperm dosages across bull studs for semen distributed from young bulls would allow the capture of sire conception rate phenotypes for young bulls that could be used to generate predictions of genetic merit for male fertility in both males and females. These data would allow genomic selection to be implemented for male fertility in addition to female fertility within the US dairy industry. While the rate of use of artificial insemination is much lower within the US beef industry, the adoption of sexed semen in the dairy industry has allowed dairy herds to select cows from which heifer replacements are produced and cows that are used to produce terminal crossbred bull calves sired by beef breed bulls. Capture of sire conception rate phenotypes in dairy herds utilizing sexed semen will contribute data enabling genomic selection for male fertility in beef cattle breeds. As the commercial sector of the beef industry increasingly adopts fixed-time artificial insemination, sire conception rate phenotypes can be captured to facilitate the development of estimates of genetic merit for male fertility within US beef breeds.