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

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.

Sperm DNA Fragmentation: Unraveling Its Imperative Impact on Male Infertility Based on Recent Evidence

Male factors may be present in up to 50-70% of infertile couples and the prevalence of male infertility accounts for 20-30% of infertility cases. Understanding the mechanisms and causes behind male infertility remains a challenge, but new diagnostic tools such as DNA fragmentation might aid in cases where the routine semen analysis is insufficient. DNA fragmentation, which refers to damages or breaks of the genetic material of the spermatozoa, is considered one of the main causes of male infertility due to impaired functional capability of sperm. The aim of the present narrative review is to investigate and enlighten the potential correlation between DNA fragmentation and male infertility parameters such as the seminal profile and the reproductive outcomes. Comprehensive research in PubMed/Medline and Scopus databases was conducted and 28 studies were included in the present review. Fourteen studies provided data regarding the impact of DNA fragmentation and seminal parameters and showed a correlation of significantly lower sperm count, lower concentration, motility, and abnormal morphology with an increased DNA fragmentation index (DFI). Similarly, 15 studies provided data regarding the impact of DFI on reproductive outcomes. Two studies showed higher aneuploidy rates with higher DFI values, and seven studies showed significantly lower pregnancy rates and live birth rates with higher DFI values. Ultimately, the studies included in this review highlight, collectively, the importance of measuring sperm DFI in the assessment of male infertility. Further studies are needed to explore the effectiveness of interventions aiming to reduce DFI levels.

Sperm long non-coding RNAs as markers for ram fertility

It is critical in sheep farming to accurately estimate ram fertility for maintaining reproductive effectiveness and for production profitability. However, there is currently a lack of reliable biomarkers to estimate semen quality and ram fertility, which is hindering advances in animal science and technology. The objective of this study was to uncover long non-coding RNAs (lncRNAs) in sperm from rams with distinct fertility phenotypes. Mature rams were allocated into two groups: high and low fertility (HF; n = 31; 94.5 ± 2.8%, LF; n = 25; 83.1 ± 5.73%; P = 0.028) according to the pregnancy rates sired by the rams (average pregnancy rate; 89.4 ± 7.2%). Total RNAs were isolated from sperm of the highest- and lowest-fertility rams (n = 4, pregnancy rate; 99.2 ± 1.6%, and 73.6 ± 4.4%, respectively) followed by next-generation sequencing of the transcripts. We uncovered 11,209 lncRNAs from the sperm of rams with HF and LF. In comparison to each other, there were 93 differentially expressed (DE) lncRNAs in sperm from the two distinct fertility phenotypes. Of these, 141 mRNAs were upregulated and 134 were downregulated between HF and LF, respectively. Genes commonly enriched for 9 + 2 motile cilium and sperm flagellum were ABHD2, AK1, CABS1, ROPN1, SEPTIN2, SLIRP, and TEKT3. Moreover, CABS1, CCDC39, CFAP97D1, ROPN1, SLIRP, TEKT3, and TTC12 were commonly enriched in flagellated sperm motility and sperm motility. Differentially expressed mRNAs were enriched in the top 16 KEGG pathways. Targets of the differentially expressed lncRNAs elucidate functions in cis and trans manner using the genetic context of the lncRNA locus, and lncRNA sequences revealed 471 mRNAs targets of 10 lncRNAs. This study illustrates the existence of potential lncRNA biomarkers that can be implemented in analyzing the quality of ram sperm and determining the sperm fertility and is used in breeding soundness exams for precision livestock farming to ensure food security on a global scale.

A comprehensive analysis of spermatozoal RNA elements in idiopathic infertile males undergoing fertility treatment

Current approaches to diagnosing male infertility inadequately assess the complexity of the male gamete. Beyond the paternal haploid genome, spermatozoa also deliver coding and non-coding RNAs to the oocyte. While sperm-borne RNAs have demonstrated potential involvement in embryo development, the underlying mechanisms remain unclear. In this study, 47 sperm samples from normozoospermic males undergoing fertility treatment using donor oocytes were sequenced and analyzed to evaluate associations between sperm RNA elements (exon-sized sequences) and blastocyst progression. A total of 366 RNA elements (REs) were significantly associated with blastocyst rate (padj < 0.05), some of which were linked to genes related to critical developmental processes, including mitotic spindle formation and both ectoderm and mesoderm specification. Of note, 27 RE-associated RNAs are predicted targets of our previously reported list of developmentally significant miRNAs. Inverse RE-miRNA expression patterns were consistent with miRNA-mediated down-regulation. This study provides a comprehensive set of REs which differ by the patient's ability to produce blastocysts. This knowledge can be leveraged to improve clinical screening of male infertility and ultimately reduce time to pregnancy.

Genome-Wide Association Analysis of Semen Characteristics in Piétrain Boars

Since artificial insemination is common practice in pig breeding, the quality and persistence of the semen are decisive for the usability of individual boars. In the current study, genome-wide association analyses were performed to investigate the genetic variability underlying phenotypic variations in semen characteristics. These traits comprise sperm morphology and sperm motility under different temporal and thermal storage conditions, in addition to standard semen quality parameters. Two consecutive samples of the fourth and fifth ejaculates from the same boar were comprehensively analyzed in a genotyped Piétrain boar population. A total of 13 genomic regions on different chromosomes were identified that contain single-nucleotide polymorphisms significantly associated with these traits. Subsequent analysis of the genomic regions revealed candidate genes described to be involved in spermatogenesis, such as FOXL3, GPER1, PDGFA, PRKAR1B, SNRK, SUN1, and TSPO, and sperm motility, including ARRDC4, CEP78, DNAAF5, and GPER1. Some of these genes were also associated with male fertility or infertility in mammals (e.g., CEP78, GPER1). The analyses based on these laboriously determined and valuable phenotypes contribute to a better understanding of the genetic background of male fertility traits in pigs and could prospectively contribute to the improvement of sperm quality through breeding approaches.

An Assessment of the Genetic Parameters of Boars' Reproductive Traits

The aim of this study was to estimate genetic parameters for the reproductive traits of boars based on single-nucleotide polymorphism data. A total of 109,836 semen samples from 2249 boars were collected between 2010 and 2022. Five basic traits were assessed: sperm volume, sperm concentration, motility, number of abnormal sperm, and, for the first time for the local population, libido. In addition, two derived traits were assessed: total sperm count and number of functional sperm. Genetic parameters were estimated using the single-step genomic best linear unbiased prediction method (ssGBLUP). Dam and sire breeds were evaluated separately. The five basic traits were evaluated using five-trait models, while the two derived traits were evaluated using single-trait models. The heritability coefficients had lower values for all sperm quality traits with both methods. For the basic traits, the heritability ranged from 0.099 to 0.342. The greatest difference between dam and sire breeds was observed for the heritability of the sperm concentration trait (0.099 and 0.271, respectively). The heritability of the libido trait was twice as high for boars of sire breeds as it was for boars of dam breeds. The genetic parameters estimated with ssGBLUP can be used in routine genetic evaluations to improve the pig breeding process.

Correlation between Post-Thaw Spermatozoa Quality of the Endangered Javan Banteng with OPN Gene Expression

The role of ex situ conservation facilities or captivity through captive breeding programs is essential in the conservation of the endangered Javan banteng. The development of semen cryopreservation may assist on one side of the conservation plan. However, the male Javan banteng reproductive capability must be considered as it influences the targeted outputs. Studying the potential biomarker for fertility such as osteopontin gene expression is also expected to help predict male fertility. Therefore, this study aimed to analyze the quality of spermatozoa after thawing to help predict the male reproductive capability of Javan banteng. Furthermore, this study investigated the potential role of osteopontin gene expression in male Javan banteng fertility. A positive reinforcement approach was used to accustom the male and female animals as we focused on establishing a collection procedure using neither sedation nor anaesthesia. Semen samples were collected at Taman Safari Indonesia, Bogor, in accordance with the female banteng receptivity. Semen samples were then evaluated and then cryopreserved under field conditions. Our study showed the different predicted reproductive capability of the Javan banteng based on the post-thaw spermatozoa quality, which showed significant differences. The OPN gene showed positive correlations with the progressive motility (r = 0.711, p = 0.048), viability (r = 0.822, p = 0.012), and acrosomal integrity (r = 0.665, p = 0.072) of Javan banteng spermatozoa after thawing. Our study demonstrated the predicted Javan banteng reproductive capability based on various post-thaw spermatozoa variables. This finding is also the first report on the OPN gene potential to be developed as the assessment tool of post-thaw spermatozoa quality of the male Javan banteng. The findings in our study may help give recommendations for future breeding programs, especially in the ex situ conservation sites.

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.

Proteomic analysis reveals the potential positive effects of Mito-TEMPO on ram sperm motility and fertility during cryopreservation

The aim of this study was to investigate the effects of mitochondria-targeted antioxidant Mito-TEMPO on the protein profile of ram sperm during cryopreservation and evaluate the cryoprotective roles of Mito-TEMPO on ram sperm quality and fertilization capacity. Semen collected from 8 Dorper rams was cryopreserved in TCG-egg yolk extender supplemented with various concentrations of Mito-TEMPO (0, 20, 40 and 60 μM). After thawing, sperm characteristics, antioxidant status and the abundance of hexose transporters (GLUT 3 and 8) were analyzed. The cervical artificial insemination (AI) was performed to evaluate the fertilization ability of cryopreserved ram sperm. The alterations of sperm proteomic profile between the control and MT₄₀ groups were determined using iTRAQ-coupled LC-MS. Supplementation with 40 μM of Mito-TEMPO resulted in the highest post-thaw sperm motility and kinematics. Sperm quality, antioxidant capacity and glucose transporter abundance of frozen-thawed ram sperm were elevated in the MT₄₀ group. The inclusion of 40 μM Mito-TEMPO in freezing extender also resulted in the higher pregnancy rate of ewes. A total of 457 proteins including 179 upregulated proteins and 278 downregulated proteins were defied as differentially expressed proteins (DEPs) using fold change (FC) > 1.2 with P < 0.05. Sixty-one DEPs with (FC > 1.5) were dramatically regulated by Mito-TEMPO. These DEPs are mainly involved in sperm motility, energy metabolism and capacitation. Our data suggest that the beneficial effects of Mito-TEMPO on sperm motility and fertility potential of cryopreserved ram semen are achieved by regulating sperm antioxidant capacity and sperm proteins related to energy metabolism and fertility.

Preimplantation embryo gene expression: 56 years of discovery, and counting

The biology of preimplantation embryo gene expression began 56 years ago with studies of the effects of protein synthesis inhibition and discovery of changes in embryo metabolism and related enzyme activities. The field accelerated rapidly with the emergence of embryo culture systems and progressively evolving methodologies that have allowed early questions to be re-addressed in new ways and in greater detail, leading to deeper understanding and progressively more targeted studies to discover ever more fine details. The advent of technologies for assisted reproduction, preimplantation genetic testing, stem cell manipulations, artificial gametes, and genetic manipulation, particularly in experimental animal models and livestock species, has further elevated the desire to understand preimplantation development in greater detail. The questions that drove enquiry from the earliest years of the field remain drivers of enquiry today. Our understanding of the crucial roles of oocyte-expressed RNA and proteins in early embryos, temporal patterns of embryonic gene expression, and mechanisms controlling embryonic gene expression has increased exponentially over the past five and a half decades as new analytical methods emerged. This review combines early and recent discoveries on gene regulation and expression in mature oocytes and preimplantation stage embryos to provide a comprehensive understanding of preimplantation embryo biology and to anticipate exciting future advances that will build upon and extend what has been discovered so far.

Molecular insights to the sperm-cervix interaction and the consequences for cryopreserved sperm

Cryopreserved ram spermatozoa are limited in their capacity to traverse the ovine cervix and achieve fertilization. This altered interaction may be related to modified molecular communication between frozen-thawed ram spermatozoa, seminal plasma, and the female tract. As such, this review aims to identify the biological processes which underpin sperm maturation and transport throughout the female reproductive tract to elucidate factors which may alter this natural process in cryopreserved ram spermatozoa. We also assess critical barriers to ram spermatozoa specific to the ovine cervix and the role of seminal plasma in mitigating these barriers. Transcriptomics is explored as a new approach to understand the sperm-cervix interaction. Recent studies have demonstrated that both spermatozoa and seminal plasma contain a complex profile of coding and non-coding RNAs. These molecular species have clear links with functional fertility, and mounting evidence suggests they may be altered by cryopreservation. Emerging in vitro cell culture models are also investigated as a "next step" in studying this interaction, utilizing transcriptomics to identify subtle changes in female tract gene expression in response to spermatozoa. The application of such models is proposed as an exciting opportunity to investigate the unique challenges faced by cryopreserved spermatozoa traversing the ovine cervix prior to fertilization.

Investigation of the mechanisms leading to human sperm DNA damage based on transcriptome analysis by RNA-seq techniques

RESEARCH QUESTION

What are the molecular mechanisms leading to human sperm DNA damage?

DESIGN

Semen samples were collected and the sperm DNA fragmentation index (DFI) was assessed. Differentially expressed RNA in spermatozoa with a high (DFI ≥30%, experimental group) or normal (DFI <30%, control group) DFI were identified by RNA-sequencing (RNA-seq) technology, and Gene Ontology and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis was performed. Three differentially expressed RNA related to sperm DNA damage and repair, namely PMS1, TP53BP1 and TLK2, were validated using real-time quantitative (RT-qPCR).

RESULTS

A total of 19,970 expressed RNA were detected in the two groups. Compared with the control group, the expression levels of 189 RNA in the experimental group were significantly increased and those of 163 genes decreased. Gene Ontology enrichment analysis showed that these RNA were mainly concentrated in the ATPase-dependent transmembrane transport complex, extracellular exosome, somatic cell DNA recombination, protein binding, cytoplasm and regulation of localization. KEGG pathway analysis showed that these RNA were mainly related to the PI3K-Akt signalling pathway, endocytosis, p53 signalling pathway and cGMP-PKG signalling pathway. The RT-qPCR results showed that the expression levels of PMS1, TP53BP1 and TLK2 in the experimental group were significantly lower than in the control group (P = 0.01, 0.015 and 0.004, respectively), which was identical to the results of RNA sequencing.

CONCLUSIONS

Differentially expressed RNA related to sperm DNA damage and repair may be identified by RNA-seq technology, which provides new insights into the understanding of sperm DNA damage and repair, and will help to discover new biomarkers related to sperm DNA damage.

Unraveling the harmful effect of oxidative stress on male fertility: A mechanistic insight

Male infertility is a widely debated issue that affects males globally. There are several mechanisms involved. Oxidative stress is accepted to be the main contributing factor, with sperm quality and quantity affected by the overproduction of free radicals. Excess reactive oxygen species (ROS) cannot be controlled by the antioxidant system and, thus, potentially impact male fertility and hamper sperm quality parameters. Mitochondria are the driving force of sperm motility; irregularities in their function may lead to apoptosis, alterations to signaling pathway function, and, ultimately, compromised fertility. Moreover, it has been observed that the prevalence of inflammation may arrest sperm function and the production of cytokines triggered by the overproduction of ROS. Further, oxidative stress interacts with seminal plasma proteomes that influence male fertility. Enhanced ROS production disturbs the cellular constituents, particularly DNA, and sperms are unable to impregnate the ovum. Here, we review the latest information to better understand the relationship between oxidative stress and male infertility, the role of mitochondria, the cellular response, inflammation and fertility, and the interaction of seminal plasma proteomes with oxidative stress, as well as highlight the influence of oxidative stress on hormones; collectively, all of these factors are assumed to be important for the regulation of male infertility. This article may help improve our understanding of male infertility and the strategies to prevent it.

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.

The ribosome inhibitor chloramphenicol induces motility deficits in human spermatozoa: A proteomic approach identifies potentially involved proteins

Mature spermatozoa are almost completely devoid of cytoplasm; as such it has long been believed that they do not contain ribosomes and are therefore not capable of synthesising proteins. However, since the 1950s, various studies have shown translational activity within spermatozoa, particularly during their in vitro capacitation. But the type of ribosomes involved (cytoplasmic or mitochondrial) is still debated. Here, we investigate the presence and activity of the two types of ribosomes in mature human spermatozoa. By targeting ribosomal RNAs and proteins, we show that both types of ribosomes are localized in the midpiece as well as in the neck and the base of the head of the spermatozoa. We assessed the impact of cycloheximide (CHX) and chloramphenicol (CP), inhibitors of cytoplasmic and mitochondrial ribosomes, respectively, on different sperm parameters. Neither CHX, nor CP impacted sperm vitality, mitochondrial activity (measured through the ATP content), or capacitation (measured through the content in phosphotyrosines). However, increasing CP concentrations induced a decrease in total and progressive motilities as well as on some kinematic parameters while no effect was observed with CHX. A quantitative proteomic analysis was performed by mass spectrometry in SWATH mode to compare the proteomes of spermatozoa capacitated in the absence or presence of the two ribosome inhibitors. Among the ∼700 proteins identified in the different tested conditions, 3, 3 and 25 proteins presented a modified abundance in the presence of 1 and 2 mg/ml of CHX, and 1 mg/ml of CP, respectively. The observed abundance variations of some CP-down regulated proteins were validated using Multiple-Reaction Monitoring (MRM). Taken together, our results are in favor of an activity of mitochondrial ribosomes. Their inhibition by CP results in a decrease in the abundance of several proteins, at least FUNDC2 and QRICH2, and consequently induces sperm motility deficits.

Developmental RNA-Seq transcriptomics of haploid germ cells and spermatozoa uncovers novel pathways associated with teleost spermiogenesis

In non-mammalian vertebrates, the molecular mechanisms involved in the transformation of haploid germ cells (HGCs) into spermatozoa (spermiogenesis) are largely unknown. Here, we investigated this process in the marine teleost gilthead seabream (Sparus aurata) through the examination of the changes in the transcriptome between cell-sorted HGCs and ejaculated sperm (SPZEJ). Samples were collected under strict quality controls employing immunofluorescence microscopy as well as by determining the sperm motion kinematic parameters by computer-assisted sperm analysis. Deep sequencing by RNA-seq identified a total of 7286 differentially expressed genes (DEGs) (p-value < 0.01) between both cell types, of which nearly half were upregulated in SPZEJ compared to HCGs. In addition, approximately 9000 long non-coding RNAs (lncRNAs) were found, of which 56% were accumulated or emerged de novo in SPZEJ. The upregulated transcripts are involved in transcriptional and translational regulation, chromatin and cytoskeleton organization, metabolic processes such as glycolysis and oxidative phosphorylation, and also include a number of ion and water channels, exchangers, transporters and receptors. Pathway analysis conducted on DEGs identified 37 different signaling pathways enriched in SPZEJ, including 13 receptor pathways, from which the most predominant correspond to the chemokine and cytokine, gonadotropin-releasing hormone receptor and platelet derived growth factor signaling pathways. Our data provide new insight into the mRNA and lncRNA cargos of teleost spermatozoa and uncover the possible involvement of novel endocrine mechanisms during the differentiation and maturation of spermatozoa.

The transcript integrity index (TII) provides a standard measure of sperm RNA

Standardizing RNA quality is key to interpreting RNA-seq data as a compromised sample can mask the underlying biology. The challenge remains when evaluating RNA quality in samples with high RNA fragmentation. For example, programmed fragmentation and cytoplasmic expulsion, integral to sperm maturation, is a prime example of the complexities of interpreting RNA-seq data, given that fragmentation can be random and\or targeted. To meet this challenge, we developed an algorithm that accurately measures RNA quality in samples with high fragmentation, such as spermatozoa. The integrity of 1,000 previously identified abundant sperm transcripts were independently visualized and evaluated using the Transcript Integrity Index (TII) algorithm to identify intact transcripts. Full-length transcripts from visual and the TII algorithm were evaluated for testis preference in humans using the GTEx tissues database. Samples were then filtered by the Interquartile Range (IQR), identifying those in which the greatest number of transcripts failed to pass the visual or TII thresholds. Transcript lists were overlapped, forming the set of intact transcripts used as TII standards. Each sample was re-evaluated as a function of this TII set of intact transcripts, with poor quality samples identified as those failing in the largest number of transcripts. While ontologically enriched in roles related to spermatogenesis and/or fertilization, samples did not segregate based on birth outcome. The TII algorithm proved an effective means to identify samples of similar quality from sperm, a cell type enriched in biologically fragmented RNAs. The algorithm should facilitate other studies using samples compromised by high levels of RNA fragmentation, such as Formalin-Fixed Paraffin-Embedded samples. Requisite to assessing male health, TII provides a solution to the long-sought-after standard that identifies samples of similar quality.

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.

Heat shock protein family D member 1 in boar spermatozoa is strongly related to the litter size of inseminated sows

Background

Sperm quality evaluation is the logical first step in increasing field fertility. Spermatozoa contain cytoplasmic organelles and biomolecules known as sperm-intrinsic factors, which play key roles in sperm maturation, sperm-oocyte fusion, and embryo development. In particular, sperm membrane proteins [e.g., arginine vasopressin receptor 2, beta-actin, prohibitin, and heat shock protein family D member 1 (HSPD1)] and RNA could be used as functional indicators of male fertility. We sought to clarify the effects of differential mRNA expression of selected genes on several fertilisation parameters, including sperm motility, motion kinematics, capacitation, and litter size, in a porcine model.

Results

Our results demonstrated that HSPD1 expression was significantly correlated with male fertility, as measured by the litter size of inseminated sows. The expression of HSPD1 mRNA was linked to sperm motility and other motion kinematic characteristics. Furthermore, HSPD1 had a 66.7% overall accuracy in detecting male fertility, and the high-litter size group which was selected with the HSPD1 marker had a 1.34 greater litter size than the low-litter size group.

Conclusions

Our findings indicate that HSPD1 might be a helpful biomarker for superior boar selection for artificial insemination, which could boost field fertility.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40104-022-00689-0.

Heritability and genetic correlation estimates of semen production traits with litter traits and pork production traits in purebred Duroc pigs

We estimated heritabilities of semen production traits and their genetic correlations with litter traits and pork production traits in purebred Duroc pigs. Semen production traits were semen volume, sperm concentration, proportion of morphologically normal sperms, total number of sperm, and total number of morphologically normal sperm. Litter traits at farrowing were total number born, number born alive, number stillborn, total litter weight at birth, mean litter weight at birth, and piglet survival rate at birth. Litter traits at weaning were litter size at weaning, total litter weight at weaning, mean litter weight at weaning, and piglet survival rate from birth to weaning. Pork production traits were average daily gain, backfat thickness, and loin muscle area. We analyzed 45,913 semen collection records of 896 boars, 6,950 farrowing performance records of 1,400 sows, 2,237 weaning performance records of 586 sows, and individual growth performance records of 9,550 animals measured at approximately 5 mo of age. Heritabilities were estimated using a single-trait animal model. Genetic correlations were estimated using a 2-trait animal model. Estimated heritabilities of semen production traits ranged from 0.20 for sperm concentration to 0.29 for semen volume and were equal to or higher than those of litter traits, ranging from 0.06 for number stillborn and piglet survival rate at birth to 0.25 for mean litter weight at birth, but lower than those of pork production traits, ranging from 0.50 for average daily gain to 0.63 for backfat thickness. In many cases, the absolute values of estimated genetic correlations between semen production traits and other traits were smaller than 0.3. These estimated genetic parameters provide useful information for establishing a comprehensive pig breeding scheme.

Genetic parameters of 5 semen production traits, 10 litter traits, and 3 pork production traits in purebred Duroc pigs was estimated. Heritabilities of semen production traits ranged from 0.20 for sperm concentration to 0.29 for semen volume and were equal to or higher than those of litter traits, ranging from 0.06 for number stillborn and piglet survival rate at birth to 0.25 for mean litter weight at birth, but lower than those of pork production traits, ranging from 0.50 for average daily gain to 0.63 for backfat thickness. In many cases, the absolute values of genetic correlations between semen production traits and other traits were smaller than 0.3. These estimated genetic parameters provide useful information for establishing a comprehensive pig breeding scheme.

Weighted single-step GWAS for body mass index and scans for recent signatures of selection in Yorkshire pigs

Controlling extra fat deposition is economically favorable in modern swine industry. Understanding the genetic architecture of fat deposition traits such as body mass index (BMI) can help in improving genomic selection for such traits. We utilized a weighted single-step genome-wide association study (WssGWAS) to detect genetic regions and candidate genes associated with BMI in a Yorkshire pig population. Three extended haplotype homozygosity (EHH)-related statistics were also incorporated within a de-correlated composite of multiple signals (DCMS) framework to detect recent selection signatures signals. Overall, the full pedigree consisted of 7,016 pigs, of which 5,561 had BMI records and 598 pigs were genotyped with an 80 K single nucleotide polymorphism (SNP) array. Results showed that the most significant windows (top 15) explained 9.35% of BMI genetic variance. Several genes were detected in regions previously associated with pig fat deposition traits and treated as potential candidate genes for BMI in Yorkshire pigs: FTMT, SRFBP1, KHDRBS3, FOXG1, SOD3, LRRC32, TSKU, ACER3, B3GNT6, CCDC201, ADCY1, RAMP3, TBRG4, CCM2. Signature of selection analysis revealed multiple candidate genes previously associated with various economic traits. However, BMI genetic variance explained by regions under selection pressure was minimal (1.31%). In conclusion, candidate genes associated with Yorkshire pigs’ BMI trait were identified using WssGWAS. Gene enrichment analysis indicated that the identified candidate genes were enriched in the insulin secretion pathway. We anticipate that these results further advance our understanding of the genetic architecture of BMI in Yorkshire pigs and provide information for genomic selection for fat deposition in this breed.

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.

In Silico Analysis of Sperm From Ejaculates with Different Semen Characteristics

Background:

Male infertility is associated with altered characteristics of the sperm within the ejaculate. It is possible to find molecular explanations for the observed phenotypes and their consequences. This study aimed to analyze, using a specialized software, a gene set of transcriptomic data from different types of ejaculates.

Methods:

Data from ejaculate samples categorized as normal, oligospermia, and teratozoospermia were obtained from Gene Expression Omnibus (GEO). After normalization, the data average for each sample category was calculated and analyzed independently using Ingenuity Pathway Analysis (IPA).

Results:

Five important canonical pathways are involved in normal and altered semen samples (Oligospermia and teratozoospermia) except sirtuin signaling and mitochondrial dysfunction pathways. The five most important biological processes are identified in all semen phenotypes, but the only difference is the genes connected with initiation of RNA transcription in oligospermic and asthenospermic samples.

Conclusion:

Surprisingly, different types of ejaculates share many pathways and biological processes; sperm proteomics as a new global approach gives clues for the development of strategies to explain the reason for observed phenotypes of ejaculated spermatozoa, their possible effect on fertility, and for implementing research strategies in the context of infertility diagnosis and treatment.

Sperm-borne sncRNAs: potential biomarkers for semen fertility?

Semen infertility or sub-fertility, whether in humans or livestock species, remains a major concern for clinicians and technicians involved in reproduction. Indeed, they can cause tragedies in human relationships or have a dramatic overall negative impact on the sustainability of livestock breeding. Understanding and predicting semen fertility issues is therefore crucial and quality control procedures as well as biomarkers have been proposed to ensure sperm fertility. However, their predictive values appeared to be too limited and additional relevant biomarkers are still required to diagnose sub-fertility efficiently. During the last decade, the study of molecular mechanisms involved in spermatogenesis and sperm maturation highlighted the regulatory role of a variety of small non-coding RNAs (sncRNAs) and led to the discovery that sperm sncRNAs comprise both remnants from spermatogenesis and post-testicular sncRNAs acquired through interactions with extracellular vesicles along epididymis. This has led to the hypothesis that sncRNAs may be a source of relevant biomarkers, associated either with sperm functionality or embryo development. This review aims at providing a synthetic overview of the current state of knowledge regarding implication of sncRNA in spermatogenesis defects and their putative roles in sperm maturation and embryo development, as well as exploring their use as fertility biomarkers.

MicroRNA-222 Transferred From Semen Extracellular Vesicles Inhibits Sperm Apoptosis by Targeting BCL2L11

Seminal plasma contains a large number of extracellular vesicles (EVs). However, the roles of these EVs and their interactions with sperm are not clear. To identify the important molecules affecting sperm motility in EVs, we analyzed RNA from seminal plasma EVs of boars with different sperm motility using whole-transcriptome sequencing and proteomic analysis. In total, 7 miRNAs, 67 lncRNAs, 126 mRNAs and 76 proteins were differentially expressed between the two groups. We observed that EV-miR-222 can obviously improve sperm motility. In addition, the results suggested that miR-222 was transferred into sperm by the EVs and that miR-222 affected sperm apoptosis by inhibiting the expression of EGFR, BCL2L11, BAX, CYCs, CASP9 and CASP3. The results of electron microscopy also showed that overexpression of miR-222 in EVs could reduce sperm apoptosis. The study of the whole transcriptomes and proteomes of EVs in boar semen revealed some miRNAs may play an important role in these EVs interactions with Duroc sperm, and the findings suggest that the release of miR-222 by semen EVs is an important mechanism by which sperm viability is maintained and sperm apoptosis is reduced. Our studies provide a new insight of miR-222 in EVs regulation for sperm motility and sperm apoptosis.

GH Overexpression Alters Spermatic Cells MicroRNAome Profile in Transgenic Zebrafish

Overexpression of growth hormone (GH) in gh-transgenic zebrafish of a highly studied lineage F0104 has earlier been reported to cause increased muscle growth. In addition to this, GH affects a broad range of cellular processes in transgenic fish, such as morphology, physiology, and behavior. Reports show changes such as decreased sperm quality and reduced reproductive performance in transgenic males. It is hypothesized that microRNAs are directly involved in the regulation of fertility potential during spermatogenesis. The primary aim of our study was to verify whether gh overexpression disturbs the sperm miRNA profile and influences the sperm quality in transgenic zebrafish. We report a significant increase in body weight of gh-transgenic males along with associated reduced sperm motility and other kinetic parameters in comparison to the non-transgenic group. MicroRNA transcriptome sequencing of gh-transgenic zebrafish sperms revealed expressions of 186 miRNAs, among which six miRNA were up-regulated (miR-146b, miR-200a-5p, miR-146a, miR-726, miR-184, and miR-738) and sixteen were down-regulated (miR-19d-3p, miR-126a-5p, miR-126b-5p, miR-22a-5p, miR-16c-5p, miR-20a-5p, miR-126b-3p, miR-107a-3p, miR-93, miR-2189, miR-202–5p, miR-221–3p, miR-125a, miR-125b-5p, miR-126a-3p, and miR-30c-5p) in comparison to non-transgenic zebrafish. Some of the dysregulated miRNAs were previously reported to be related to abnormalities in sperm quality and reduced reproduction ability in other species. In this study, an average of 134 differentially expressed miRNAs-targeted genes were predicted using the in silico approach. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis demonstrated that the genes of affected pathways were primarily related to spermatogenesis, sperm motility, and cell apoptosis. Our results suggested that excess GH caused a detrimental effect on sperm microRNAome, consequently reducing the sperm quality and reproductive potential of zebrafish males.

Single paternal dexamethasone challenge programs offspring metabolism and reveals multiple candidates in RNA-mediated inheritance

Single traumatic events that elicit an exaggerated stress response can lead to the development of neuropsychiatric conditions. Rodent studies suggested germline RNA as a mediator of effects of chronic environmental exposures to the progeny. The effects of an acute paternal stress exposure on the germline and their potential consequences on offspring remain to be seen. We find that acute administration of an agonist for the stress-sensitive Glucocorticoid receptor, using the common corticosteroid dexamethasone, affects the RNA payload of mature sperm as soon as 3 hr after exposure. It further impacts early embryonic transcriptional trajectories, as determined by single-embryo sequencing, and metabolism in the offspring. We show persistent regulation of tRNA fragments in sperm and descendant 2-cell embryos, suggesting transmission from sperm to embryo. Lastly, we unravel environmentally induced alterations in sperm circRNAs and their targets in the early embryo, highlighting this class as an additional candidate in RNA-mediated inheritance of disease risk.

Advances in sperm analysis: techniques, discoveries and applications

Infertility affects one in six couples worldwide, and fertility continues to deteriorate globally, partly owing to a decline in semen quality. Sperm analysis has a central role in diagnosing and treating male factor infertility. Many emerging techniques, such as digital holography, super-resolution microscopy and next-generation sequencing, have been developed that enable improved analysis of sperm motility, morphology and genetics to help overcome limitations in accuracy and consistency, and improve sperm selection for infertility treatment. These techniques have also improved our understanding of fundamental sperm physiology by enabling discoveries in sperm behaviour and molecular structures. Further progress in sperm analysis and integrating these techniques into laboratories and clinics requires multidisciplinary collaboration, which will increase discovery and improve clinical outcomes.

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 BoviPure^TM on porcine sperm. Four boar ejaculates were purified with BoviPure^TM 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.

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.

Seminal Plasma: Relevant for Fertility?

Seminal plasma (SP), the non-cellular component of semen, is a heterogeneous composite fluid built by secretions of the testis, the epididymis and the accessory sexual glands. Its composition, despite species-specific anatomical peculiarities, consistently contains inorganic ions, specific hormones, proteins and peptides, including cytokines and enzymes, cholesterol, DNA and RNA-the latter often protected within epididymis- or prostate-derived extracellular vesicles. It is beyond question that the SP participates in diverse aspects of sperm function pre-fertilization events. The SP also interacts with the various compartments of the tubular genital tract, triggering changes in gene function that prepares for an eventual successful pregnancy; thus, it ultimately modulates fertility. Despite these concepts, it is imperative to remember that SP-free spermatozoa (epididymal or washed ejaculated) are still fertile, so this review shall focus on the differences between the in vivo roles of the SP following semen deposition in the female and those regarding additions of SP on spermatozoa handled for artificial reproduction, including cryopreservation, from artificial insemination to in vitro fertilization. This review attempts, including our own results on model animal species, to critically summarize the current knowledge of the reproductive roles played by SP components, particularly in our own species, which is increasingly affected by infertility. The ultimate goal is to reconcile the delicate balance between the SP molecular concentration and their concerted effects after temporal exposure in vivo. We aim to appraise the functions of the SP components, their relevance as diagnostic biomarkers and their value as eventual additives to refine reproductive strategies, including biotechnologies, in livestock models and humans.

Sperm Methylome Profiling Can Discern Fertility Levels in the Porcine Biomedical Model

A combined Genotyping By Sequencing (GBS) and methylated DNA immunoprecipitation (MeDIP) protocol was used to identify—in parallel—genetic variation (Genomic-Wide Association Studies (GWAS) and epigenetic differences of Differentially Methylated Regions (DMR) in the genome of spermatozoa from the porcine animal model. Breeding boars with good semen quality (n = 11) and specific and well-documented differences in fertility (farrowing rate, FR) and prolificacy (litter size, LS) (n = 7) in artificial insemination programs, using combined FR and LS, were categorized as High Fertile (HF, n = 4) or Low Fertile (LF, n = 3), and boars with Unknown Fertility (UF, n = 4) were tested for eventual epigenetical similarity with those fertility-proven. We identified 165,944 Single Nucleotide Polymorphisms (SNPs) that explained 14–15% of variance among selection lines. Between HF and LF individuals (n = 7, 4 HF and 3 LF), we identified 169 SNPs with p ≤ 0.00015, which explained 58% of the variance. For the epigenetic analyses, we considered fertility and period of ejaculate collection (late-summer and mid-autumn). Approximately three times more DMRs were observed in HF than in LF boars across these periods. Interestingly, UF boars were clearly clustered with one of the other HF or LF groups. The highest differences in DMRs between HF and LF experimental groups across the pig genome were located in the chr 3, 9, 13, and 16, with most DMRs being hypermethylated in LF boars. In both HF and LF boars, DMRs were mostly hypermethylated in late-summer compared to mid-autumn. Three overlaps were detected between SNPs (p ≤ 0.0005, n = 1318) and CpG sites within DMRs. In conclusion, fertility levels in breeding males including FR and LS can be discerned using methylome analyses. The findings in this biomedical animal model ought to be applied besides sire selection for andrological diagnosis of idiopathic sub/infertility.

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.

Seasonal Changes of Nuclear DNA Fragmentation in Boar Spermatozoa in Spain

Simple Summary

Artificial insemination is widely used in pig production and currently a boar performs several thousand matings per year. Traditionally sperm quality is focused on the number of spermatozoa, their motility and morphology. However, the quality of sperm DNA, which contains genetic information, is also related to fertility problems. The aim of this research was to study the effect of natural light hours and age of the boar on the status of the sperm DNA. After a powerful statistical analysis, it was found that the percentage of spermatozoa with fragmented DNA decreases within the observed age range as the boar gets older. On the other hand, the amount of spermatozoa with fragmented DNA was the lowest in autumn while it was the highest in summer. This study demonstrates the remaining seasonality of boars in Spain and highlights the importance of controlling the environmental conditions in the farms. Sperm DNA testing provides a basis for improving the selection of AI boars by excluding males with higher DNA fragmentation due to their very young reproductive age that may pose a potential subfertility.

Abstract

There are numerous cases when conventional spermiogram parameters are all within an acceptable range but boar subfertility persists. The total sperm nuclear DNA fragmentation index (tDFI) is a trait related to fertility and prolificacy problems that is not routinely evaluated in commercial AI boars. The aim of this research was to study the effect of the photoperiod, season and reproductive age of the boar on tDFI (measured by SCSA) of 1279 ejaculates from 372 different boars belonging to 6 different breeds located in 6 AI studs in Spain. tDFI data ranged from 0.018% to 20.1%. Although there was a significant single boar effect in the tDFI occurrence, a negative correlation between the tDFI and the age of the boar was found (p < 0.001). tDFI would decrease due to aging of the boar 0.66% each year old within the observed age range. After including age as a covariate in the ANCOVA, no differences were found in tDFI between photoperiods when the sperm collection date was evaluated. However, when the date of the production of semen in the testis was evaluated, the total percentage of spermatozoa with fragmented nuclear DNA was 1.46% higher in the increasing photoperiod in comparison to the decreasing photoperiod (p < 0.0001). On the other hand, for both dates, the lowest tDFI values corresponded to minimum day length for decreasing photoperiod phase (autumn), while the highest tDFI values were found in summer (maximum day length for decreasing photoperiod phase).

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.

Characterisation of sperm piRNAs and their correlation with semen quality traits in swine

Piwi-interacting RNAs (piRNAs) are a class of non-coding RNAs that are essential in the transcriptional silencing of transposable elements and warrant genome stability in the mammalian germline. In this study, we have identified piRNAs in porcine sperm using male germline and zygote datasets from human, mice, cow and pig, and evaluated the relation between their abundances and sperm quality traits. In our analysis, we identified 283 382 piRNAs, 1355 of which correlated with P ≤ 0.01 to at least one semen quality trait. Fifty-seven percent of the correlated piRNAs mapped less than 50 kb apart from any other piRNA in the pig genome. Furthermore, piRNA location was significantly enriched near long interspersed nuclear elements. Moreover, some of the significant piRNAs mapped within or close to genes relevant for fertility or spermatogenesis such as CSNK1G2 and PSMF1.

Regulatory Potential of Long Non-Coding RNAs (lncRNAs) in Boar Spermatozoa with Good and Poor Freezability

Long non-coding RNAs (lncRNAs) are suggested to play an important role in the sperm biological processes. We performed de novo transcriptome assembly to characterize lncRNAs in spermatozoa, and to investigate the role of the potential target genes of the differentially expressed lncRNAs (DElncRNAs) in sperm freezability. We detected approximately 4007 DElncRNAs, which were differentially expressed in spermatozoa from boars classified as having good and poor semen freezability (GSF and PSF, respectively). Most of the DElncRNAs were upregulated in boars of the PSF group and appeared to significantly affect the sperm's response to the cryopreservation conditions. Furthermore, we predicted that the potential target genes were regulated by DElncRNAs in cis or trans. It was found that DElncRNAs of both freezability groups had potential cis- and trans-regulatory effects on different protein-coding genes, such as COX7A2L, TXNDC8 and SOX-7. Gene Ontology (GO) enrichment revealed that the DElncRNA target genes are associated with numerous biological processes, including signal transduction, response to stress, cell death (apoptosis), motility and embryo development. Significant differences in the de novo assembled transcriptome expression profiles of the DElncRNAs between the freezability groups were confirmed by quantitative real-time PCR analysis. This study reveals the potential effects of protein-coding genes of DElncRNAs on sperm functions, which could contribute to further research on their relevance in semen freezability.

Expression Pattern of Seminal Plasma Extracellular Vesicle Small RNAs in Boar Semen

Extracellular vesicles (EVs) regulate multiple physiological processes. Seminal plasma contains numerous EVs that may deliver functional molecules such as small RNAs (sRNAs) to the sperm. However, the RNA profiles in the boar seminal plasma extracellular vesicles (SP-EVs) and its function have not been characterized. The aim of this study was to characterize the functions and sRNA profiles in the boar SP-EVs using deep sequencing technology. Briefly, boar SP-EVs were isolated by differential ultracentrifugation and confirmed with a transmission electron microscope (TEM), nanoparticle tracking analysis (NTA), and Western blot. The isolated boar SP-EVs contained numerous and diverse sRNA families, including microRNAs (miRNAs, 9.45% of the total reads), PIWI-interacting RNAs (piRNAs, 15.25% of the total reads), messenger RNA fragments (mRNA, 25.30% of the total reads), and tRNA-derived small RNAs (tsRNA, 0.01% of the total reads). A total of 288 known miRNAs, 37 novel miRNA, and 19,749 piRNAs were identified in boar SP-EVs. The identified ssc-miR-21-5p may confer negative effects on sperm fertility based on a dual-luciferase reporter experiment. This study therefore provides an effective method to isolate SP-EVs and characterizes the sRNA profile.

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.

Period of Boar Ejaculate Collection Contributes to the Yearly Intra-Male Variability of Seminal Plasma Cytokines

The concentrations of cytokines in seminal plasma (SP) fluctuate over time in healthy males, weakening their practical usefulness as diagnostic tools. This study evaluated the relevance of intra-male variability in SP cytokines and to what extent the period of the year when ejaculate is collected contributes to such variability. Thirteen cytokines (GM-CSF, IFNγ, IL-1α, IL-1β, IL-1ra, IL-2, IL-4, IL-6, IL-8, IL-10, IL-12, IL-18, and TNFα) were measured using the Luminex xMAP^® technology for 180 SP samples of ejaculate collected over a year from nine healthy and fertile boars. The SP samples were grouped into two annual periods according to decreasing or increasing daylight and ambient temperature. Intra-male variability was higher than inter-male variability for all cytokines. All SP cytokines showed concentration differences between the two periods of the year, showing the highest concentration during the increasing daylength/temperature period, irrespective of the male. Similarly, some cytokines showed differences between daylength/temperature periods when focusing on their total amount in the ejaculate. No strong relationship (explaining more than 50% of the total variance) was found between annual fluctuations in SP-cytokine levels and semen parameters. In conclusion, the period of the year during which ejaculates were collected helps explain the intra-male variability of SP-cytokine levels in breeding boars.

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.

Progress in Research on Sperm DNA Fragmentation

With the increasing incidence of male infertility, routine detection of semen is insufficient to accurately assess male fertility. Infertile men, who have lower odds of conceiving naturally, exhibit high levels of sperm DNA fragmentation (SDF). The mechanisms driving SDF include abnormal spermatogenesis, oxidative stress damage, and abnormal sperm apoptosis. As these factors can induce SDF and subsequent radical changes leading to male infertility, detection of the extent of SDF has become an efficient routine method for semen analysis. Although it is still debated, SDF detection has become a research hotspot in the field of reproductive medicine as a more accurate indicator for assessing sperm quality and male fertility. SDF may be involved in male infertility, reproductive assisted outcomes, and growth and development of offspring. The effective detection methods of SDF are sperm chromatin structure analysis (SCSA), terminal transferase-mediated dUTP end labeling (TUNEL) assay, single-cell gel electrophoresis (SCGE) assay, and sperm chromatin dispersion (SCD) test, and all of these methods are valuable for assisted reproductive techniques. Currently, the preferred method for detecting sperm DNA integrity is SCSA. However, the regulation network of SDF is very complex because the sperm DNA differs from the somatic cell DNA with its unique structure. A multitude of molecular factors, including coding genes, non-coding genes, or methylated DNA, participate in the complex physiological regulation activities associated with SDF. Studying SDF occurrence and the underlying mechanisms may effectively improve its clinical treatments. This review aimed to outline the research status of SDF mechanism and detection technology-related issues, as well as the effect of increased SDF rate, aiming to provide a basis for clinical male infertility diagnosis and treatment.

Identification and quantification of coding and long non-coding RNAs in stallion spermatozoa separated by density

BACKGROUND

It is not unusual for stallions to have fertility problems. For many, artificial insemination with more dense spermatozoa (isolated by density gradient centrifugation) results in greater pregnancy rates compared with the rates when using unfractionated spermatozoa. RNAs in spermatozoa delivered to the oocyte at conception are required for embryo development. Novel molecular assays of spermatozoa that reflect function are needed to predict the fertility of stallions.

OBJECTIVES

To describe and compare the RNA populations in more dense and less dense spermatozoa from stallions.

MATERIALS AND METHODS

Spermatozoa from five stallions were separated into more dense and less dense populations by density gradient centrifugation. Complementary DNA libraries were made from each of the ten total RNA samples after ribosomal RNA removal. Next-generation sequencing characterized the RNA populations in more and less dense spermatozoa. Quantitative reverse transcription-PCR was used to confirm differential expression of selected RNAs.

RESULTS

Stallion spermatozoa contain 11 215 RNAs, with the most prevalent RNA being a 1492 base long non-coding RNA. The levels of 159 RNAs were greater in more dense spermatozoa, while levels of seven other RNAs were greater in less dense spermatozoa. Quantitative reverse transcription-PCR confirmed the threefold greater levels of solute carrier family 26 member 8 (SLC26A8) mRNA in less dense spermatozoa, and sixfold and threefold greater expression levels of the SCP2 sterol binding domain containing 1 (SCP2D1) and spermatogenesis-associated protein 31D1 (SPATA31D1) mRNAs in more dense spermatozoa, respectively.

DISCUSSION AND CONCLUSION

We identified 11 215 RNAs in stallion spermatozoa and 166 with differential expression between more dense and less dense fractions. Many prevalent RNAs were also found in bull, boar, and human spermatozoa. Many differentially expressed RNAs are known to be testis- or spermatozoa-specific. Our results may lead to identification of an RNA population in spermatozoa that is optimal for establishing successful pregnancies.

The Transcriptome of Pig Spermatozoa, and Its Role in Fertility

In the study presented here we identified transcriptomic markers for fertility in the cargo of pig ejaculated spermatozoa using porcine-specific micro-arrays (GeneChip^® miRNA 4.0 and GeneChip^® Porcine Gene 1.0 ST). We report (i) the relative abundance of the ssc-miR-1285, miR-16, miR-4332, miR-92a, miR-671-5p, miR-4334-5p, miR-425-5p, miR-191, miR-92b-5p and miR-15b miRNAs, and (ii) the presence of 347 up-regulated and 174 down-regulated RNA transcripts in high-fertility breeding boars, based on differences of farrowing rate (FS) and litter size (LS), relative to low-fertility boars in the (Artificial Insemination) AI program. An overrepresentation analysis of the protein class (PANTHER) identified significant fold-increases for C-C chemokine binding (GO:0019957): CCR7, which activates B- and T-lymphocytes, 8-fold increase), XCR1 and CXCR4 (with ubiquitin as a natural ligand, 1.24-fold increase), cytokine receptor activity (GO:0005126): IL23R receptor of the IL23 protein, associated to JAK2 and STAT3, 3.4-fold increase), the TGF-receptor (PC00035) genes ACVR1C and ACVR2B (12-fold increase). Moreover, two micro-RNAs (miR-221 and mir-621) were down- and up-regulated, respectively, in high-fertility males. In conclusion, boars with different fertility performance possess a wide variety of differentially expressed RNA present in spermatozoa that would be attractive targets as non-invasive molecular markers for predicting fertility.

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.

Sperm miR-15a and miR-29b are associated with bull fertility

MicroRNAs modulate male fertility by regulating gene expression. In this study, dynamics of sperm miR-15a, miR-29b and miR-34a from high fertility (HF) and low fertility (LF) bulls using RT-qPCR were evaluated. Bioinformatic tools were employed to ascertain genes of interest of the sperm miRNAs. The expression levels of p53, BCL2, BAX and DNMT1 in bull spermatozoa were determined by immunoblotting. MicroRNA levels of miR-15a and miR-29 were higher in LF sires when compared with those present in HF bulls. Expression levels of miR-34a did not differ between the two groups. We found an inverse correlation between miR-15a and bull fertility. MiR29-b was also negatively associated with fertility scores. BCL2 and DNMT1 were higher in HF bulls while BAX was higher in the LF group. Our data showed a positive correlation between BCL2 and bull fertility. In addition, DNMT1 was positively associated with bull fertility. Furthermore, levels of BAX were negatively linked with bull fertility scores. Identification of miRNAs found in the spermatozoa of sires with different in vivo fertility helps understand the alterations in the fertilising capacity from cattle and other mammals. These potential biomarkers can be used in reproductive biotechnology as fertility markers to assess semen quality and predict male fertility.