Abstract In bovines, cryopreserved semen is used for artificial insemination; however, the fertility of cryopreserved semen is far lower than that of fresh semen. Although cryopreservation alters sperm phenotypic characteristics, its effect on sperm molecular health is not thoroughly understood. The present study applied next-generation sequencing to investigate the effect of cryopreservation on the sperm transcriptomic composition of bull spermatozoa. While freshly ejaculated bull spermatozoa showed 14,280 transcripts, cryopreserved spermatozoa showed only 12,375 transcripts. Comparative analysis revealed that 241 genes were upregulated, 662 genes were downregulated, and 215 genes showed neutral expression in cryopreserved spermatozoa compared to fresh spermatozoa. Gene ontology analysis indicated that the dysregulated transcripts were involved in nucleic acid binding, transcription-specific activity, and protein kinase binding involving protein autophosphorylation, ventricular septum morphogenesis, and organ development. Moreover, the dysregulated genes in cryopreserved spermatozoa were involved in pathways associated with glycogen metabolism, MAPK signalling, embryonic organ morphogenesis, ectodermal placode formation, and regulation of protein auto-phosphorylation. These findings suggest that the cryopreservation process induced alterations in the abundance of sperm transcripts related to potential fertility-associated functions and pathways, which might partly explain the reduced fertility observed with cryopreserved bull spermatozoa. Keywords: cryopreservation, spermatozoa, pathways, transcripts, fertility, bovine Introduction Artificial insemination using cryopreserved semen is widely and routinely practiced for the genetic improvement of livestock. Ultralow freezing allows the preservation of semen for several years, which can be used later for artificial breeding ([40]Barbas and Mascarenhas, 2009). Although the cryopreservation of sperm is an important process in assisted reproduction technologies, the fertility of cryopreserved spermatozoa is reportedly not as good as fresh spermatozoa ([41]Watson, 1995; [42]Bailey et al., 2000; [43]Kadirvel et al., 2009). Several reasons are attributed to the decreased fertility of cryopreserved semen, primarily damage that occurs in spermatozoa during cryopreservation (cryodamage), which alters its fertilising potential ([44]Cormier et al., 1997; [45]Watson, 2000; [46]Yeste, 2016). The process of cryopreservation results in the death of almost 50% of spermatozoa, while the remaining sperm population shows altered functional competencies ([47]Kumaresan et al., 2011; [48]Kumaresan et al., 2012; [49]Singh et al., 2016; [50]Kumaresan et al., 2017; [51]Saraf et al., 2019; [52]Rather et al., 2020; [53]Vignesh et al., 2020; [54]Nag et al., 2021), which might be linked to the reduced fertility of cryopreserved semen. However, the fertility of cryopreserved semen subjected to sperm selection methods (for the selection of viable, active, and phenotypically normal spermatozoa) was also not as good as that of fresh semen, although some improvement in fertility was observed ([55]Said and Land, 2011; [56]Marzano et al., 2020). Therefore, besides cryopreservation-induced sperm structural and functional alterations, other inherent factors in spermatozoa might also be altered during cryopreservation as cryopreserved spermatozoa with normal phenotypic characteristics also show altered fertilising potential compared to fresh spermatozoa ([57]Kadirvel et al., 2009; [58]Elango et al., 2021). Therefore, the assessment of the molecular alterations induced by the process of cryopreservation is essential for understanding the decreased fertility associated with cryopreserved semen. Sperm was previously believed to deliver only the paternal DNA to the oocyte after fertilisation. However, several recent studies have demonstrated the role of sperm transcripts in male fertility ([59]Paul et al., 2020; [60]Prakash et al., 2021; [61]Saraf et al., 2021) and subsequent embryonic development ([62]Zhang et al., 2018). After fertilisation, spermatozoal mRNAs are transferred to the oocyte and play important roles in embryonic development, morphogenesis, and implantation ([63]Ostermeier et al., 2004). Bull spermatozoa harbours a repertoire of transcripts ([64]Prakash et al., 2021) that vary with semen quality and fertility ([65]Karuthadurai et al., 2022). Earlier studies reported different transcriptomic profiles between epididymal and ejaculated spermatozoa, and also among different seasons ([66]AL-Sahaf, 2012), indicating the influence of the micro-environment on sperm transcripts. Any changes in the transcript content of the sperm affect sperm properties such as motility, DNA intactness, and acrosome integrity ([67]Bissonnette et al., 2009). Although structural damage to sperm during cryopreservation has been studied extensively, the molecular alterations induced by the process of cryopreservation are not well documented. We hypothesised that the transcriptomic composition of sperm could be influenced by the process of cryopreservation, which might be associated with the reduced fertility of cryopreserved semen. Therefore, in the present study, we assessed the global transcriptomic composition of spermatozoa immediately after ejaculation and after cryopreservation using next-generation sequencing. We compared the transcriptomic profile of fresh spermatozoa to that of cryopreserved spermatozoa and identified common and dysregulated transcripts. Using functional annotation of these genes, we report the alterations in sperm transcripts and important pathways associated with fertility, which were induced by the cryopreservation process. Materials and methods Ethical approval statement The current study was carried out at the Theriogenology Laboratory of the Southern Regional Station of ICAR-National Dairy Research Institute, Bengaluru, Karnataka. All the experiments and procedures performed in this study were approved by the Animal Ethical Committee of the institute (CPCSEA/IAEC/LA/SRS-ICAR-NDRI-2019/No.04). Experimental bulls and sample preparation This investigation was conducted on Holstein-Friesian crossbred bulls (n = 6; age 4–6 years). All experimental bulls had passed breeding soundness evaluations and were routinely used for artificial breeding. Ejaculates were collected using an artificial vagina; after preliminary evaluation, only ejaculates with minimum sperm concentrations of 600 million/ml, +3.0 mass activity (0–5 scale), ≥70% progressive individual motility, and <20% sperm abnormalities were utilised for further processing. The ejaculates were divided into two aliquots; one aliquot was used fresh while the other aliquot was subjected to cryopreservation as per the standard procedure. Briefly, the ejaculates were diluted using pre-warmed (34°C) Tris-egg yolk glycerol extender (20% egg yolk and 7% glycerol fractions) and then further processed for cryopreservation. The diluted semen was then filled and sealed in 0.25 ml mini straws (20 × 10^6 sperm per dose) using an automatic filling and sealing machine. The straws were then equilibrated in a cold handling cabinet (IMV Technologies, France) for 4 h at 4°C. Post-equilibration, the doses were loaded into a programmable Biofreezer (Digitcool, IMV Technologies, France) for cryopreservation as per standard protocol. After reaching −140°C, the straws were directly plunged into liquid nitrogen. After cryopreservation, the frozen semen was thawed at 37°C for 30 s and used for further processing. Ejaculates from three bulls were pooled to obtain one representative sample. Therefore, we obtained two representative samples from six bulls for each condition (fresh and cryopreserved), which were individually subjected to transcriptomic analysis. RNA extraction and synthesis of cDNA Discontinuous Percoll gradients (90–45%) were used to fractionate pure sperm by eliminating epithelial cells and seminal plasma. The method described by Parthipan et al. ([68]2015), with minor modifications was used for the extraction of total sperm RNA from fresh and cryopreserved samples using TRIzol (Ambion, Thermo Fisher Scientific, United States). The RNA was quantified using Nanodrop (ND-1000, Thermo-scientific, United States). For cDNA library preparation, RNA samples with 260/280 ratios of 1.85–2.0 were used. An initial concentration of 50–100 ng of RNA was used for cDNA synthesis using the RevertAid First Strand cDNA Synthesis Kit (Thermo Fisher Scientific, United States). The synthesised cDNA was stored at −20°C until further use. Transcriptomics library preparation Total RNA (1 μg) enriched for mRNA using the NEB Magnetic mRNA Isolation Kit (Illumina, United States). RNA library preparation was performed using the NEB Ultra II RNA library prep kit (Illumina, United States) and sequencing (paired-end technology) was performed on an Illumina NextSeq 500 instrument (Illumina, United States. The enriched mRNA was fragmented to 200 bp using fragmentation buffer. Complementary RNA sequence hybridisation was performed by adding random hexamer primers. Reverse transcriptase enzyme and dNTPs were used for synthesising the first strand of cDNA from fragments. DNA polymerase I and RNase H were used to convert the single-strand cDNA into double-stranded cDNA, which was purified using 1.8x AMPure beads. Adaptor-ligated cDNA was purified using AMPure beads and was enriched with specific primers for sequencing on the Illumina platforms. RNA sequencing and data analysis The Galaxy online server tool was used to analyse the sequences. The quality of the raw generated data was checked using the Fast QC program with a Phred quality score cut-off of Q30. The Cutadapt tool was used to remove adaptor sequences from the FASTQ files. HISAT2 (version2.1) was used to align the sequences to the reference genome. The samples were aligned to the reference genome of Bos taurus (version UCD 3.1.94). Cufflinks (version—2.2.1.2) was used to identify and estimate the abundance of the transcripts. After normalisation, the transcript expression levels were calculated as FPKM (fragments per kilobase of transcript per million mapped reads). Full-length transcript analysis was performed using the depth of coverage program version-0.0.2 in GATK. Functional annotation and Gene Ontology analysis Functional annotation and gene ontology (GO) analysis were performed using the Database for Annotation, Visualization, and Integrated Discovery (DAVID) Bioinformatics Resources (v6.8) ([69]https://david.ncifcrf.gov/). Molecular Function (MF), Biological Process (BP), cellular components (CC), and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were performed to identify the genes with variations. For network analysis, genes related to sperm quality were selected from MF, BP, CC, and KEGG. ClueGo (Version 2.5.4) plugins in the open-source Cytoscape software (version 3.7.1) ([70]Cluego.org) were used to identify the interactions of novel genes and their associated pathways. Statistical analysis DAVID was used to perform functional enrichment analysis of the Gene Ontology categories and KEGG pathway elements. The EASE score was calculated using the Fisher’s exact p-value based on the input list of genes and genes participating in a certain pathway, resulting in a strongly enriched word. A p-value threshold of 0.01 was used to determine the EASE score cut-off. After enrichment analysis, terms with p-values <0.01 were identified as enriched terms. Thus, the lower the p-value, the more enriched the phrase. At least two genes were required for the examination of connection. The corrected p-values obtained after adjustments using the Benjamini and Bonferroni tests were defined as the false discovery rate (FDR). The processed reads were mapped to the reference genome for analysis of the raw transcriptomic data. The number of reads mapped to the exonic region indicated the gene expression. The read counts were used by DESeq to determine the differential expression based on the number of mapped reads. DESeq was used to create size factors (using the size factors function) and fit the data with a negative binomial distribution (using the nbinom test function) to compare the control (fresh) and treated (cryopreserved) groups and return log2-fold changes and significant p-values. Results Differentially regulated transcripts The primary analysis of fresh and cryopreserved spermatozoa resulted in 14,280 and 12,375 transcripts, respectively ([71]Supplementary File S1 and [72]Supplementary Figure S1) indicating that the process of cryopreservation altered the sperm transcriptomic profile. After normalisation of the read counts, a total of 1118 genes were common between fresh and cryopreserved spermatozoa. In cryopreserved spermatozoa, 241 genes were upregulated, 662 genes were downregulated, and 215 genes were neutral (neither up- regulated nor downregulated) in expression compared to fresh spermatozoa ([73]Figure 1). The genes that were upregulated or downregulated based on log2 fold-change values are shown in a volcano plot ([74]Figure 2). A heat map of the top 20 upregulated and top 20 downregulated genes in the cryopreserved sample ([75]Figure 3) showed the genes with highly dysregulated expression. C11H9orf50, ARL2BP, and HSPB1 were the most downregulated, while RPRD2, ENSBTAG00000030892, and ENSBTAG00000030838 were the most upregulated in cryopreserved spermatozoa. FIGURE 1. FIGURE 1 [76]Open in a new tab Venn diagram of the transcriptomic compositions of fresh and cryopreserved bull spermatozoa. FIGURE 2. FIGURE 2 [77]Open in a new tab Volcano plot representing the gene log-fold changes and p-value. Green dots: significantly differentially expressed genes. FIGURE 3. [78]FIGURE 3 [79]Open in a new tab Heat map representing the top 20 upregulated (green) and top 20 downregulated (red and yellow) genes in cryopreserved spermatozoa compared to fresh spermatozoa. Gene Ontology The dysregulated genes were selected for Gene Ontology (GO) analysis based on a log2 fold-change cut-off of ±1. These dysregulated genes were then subjected to DAVID analysis to annotate GO terms including molecular function (MF), biological process (BP), and cellular components (CC). A total of 54 upregulated and 411 downregulated genes were annotated for GO analysis. From these, the top MF, BP, and CC for upregulated and downregulated genes ([80]Figures 4A,B) ([81]Supplementary File S2) were selected based on the counts of genes involved in each GO term. The important upregulated genes in the cryopreserved spermatozoa were involved in specific DNA binding, nucleic acid binding, transcription-specific activity, protein kinase binding, core promotor sequence-specific DNA binding, and inward rectifier potassium channel activity. These were mostly localised in the nucleus, cytosol, and potassium channel complex and were involved in biological processes such as the regulation of transcription from the RNA polymerase II promotor, the regulation of neuronal differentiation, the positive regulation of protein autophosphorylation, ventricular septum morphogenesis, organ development, Rap protein signal transduction, neuroblast proliferation, and negative regulation of myeloid cell differentiation. The downregulated genes were involved in poly(A) RNA binding, structural constituents of ribosomes, DNA binding, protein binding, RNA binding, ubiquitin protein ligase binding, and nucleotide binding. FIGURE 4. [82]FIGURE 4 [83]Open in a new tab Gene Ontology chart of differentially expressed genes between fresh and cryopreserved spermatozoa. All three components (molecular function, biological process, and cellular functions) are represented in the donut chart. (A) Upregulated genes. (B) Downregulated genes. KEGG pathway and network analysis KEGG pathway enrichment analysis of the downregulated genes ([84]Supplementary File S3) in cryopreserved spermatozoa revealed the involvement of genes in important pathways related to ribosomes (50 genes), Huntington’s disease (17 genes), thermogenesis (16 genes), spliceosomes (16 genes), and endocytosis (16 genes). KEGG pathway enrichment analysis of the upregulated genes ([85]Supplementary File S3) in cryopreserved spermatozoa revealed the involvement of genes in pathways related to the neurotrophin signalling pathway (4 genes), transcriptional misregulation in cancer (4 genes), and neuroactive ligand-receptor interaction (5 genes). A network analysis of GO terms and KEGG pathways for upregulated and downregulated genes related to fertility was also performed. The network analysis of downregulated ([86]Figure 5) genes in cryopreserved spermatozoa showed their involvement in pathways such as glycogen metabolism, the MAPK signalling pathway, fluid shear and stress atherosclerosis, cytoskeletal protein binding, the oestrogen signalling pathway, amphetamine addiction, transcription factor binding, Huntington’s disease, ubiquitin protein ligase binding, translation factor activity RNA binding, mRNA binding, transcription regulation region sequence-specific DNA binding, RNA polymerase II core promoter sequence-specific DNA binding, oxidoreductase activity, and glycolysis/gluconeogenesis. The network analysis of upregulated ([87]Figure 6) genes in cryopreserved spermatozoa showed their involvement in pathways such as embryonic organ morphogenesis, ectodermal placode formation, regulation of protein auto-phosphorylation, and death receptor binding ([88]Supplementary File S4). FIGURE 5. [89]FIGURE 5 [90]Open in a new tab Network analysis of all downregulated genes in the cryopreserved samples. MAPK signalling pathway, fluid shear and stress atherosclerosis, cytoskeletal protein binding, and the oestrogen signalling pathway are mapped. FIGURE 6. [91]FIGURE 6 [92]Open in a new tab Network analysis of all upregulated genes in the cryopreserved samples. Embryonic organ morphogenesis, ectodermal placode formation, regulation of protein auto-phosphorylation, and death receptor binding are mapped. Discussion The process of cryopreservation alters sperm phenome and functions, which affect the quality and fertility of spermatozoa. Recent studies have shown that sperm transcripts are associated with sperm quality and fertility, as well as embryonic development ([93]Zhang et al., 2018; [94]Paul et al., 2020; [95]Prakash et al., 2021). Since the process of cryopreservation alters sperm structure, this process may also affect the transcriptomic composition of sperm. The results of the present study demonstrated the significant alterations in the transcriptomic composition of bull sperm by the process of cryopreservation, which might explain the reduced fertility of cryopreserved semen. The mature spermatozoal RNA undergoes certain cellular modifications during spermiogenesis. The required mRNAs are then packed in spermatozoa before it reaches the transcriptionally nascent stage in the epididymis ([96]Grunewald et al., 2005; [97]Vijayalakshmy et al., 2018; [98]Paul et al., 2020) for delivery to the oocyte after fertilisation. We observed highly downregulated expression of C11H9orf50, ARL2BP, and HSPB1 and highly upregulated expression of RPRD2, ENSBTAG00000030892, and ENSBTAG00000030838 in cryopreserved spermatozoa. ARL2BP is required for ciliary microtubule structure. A mouse ARL2BP knockout showed structural impairments such as abnormal head and misassembled tail of sperm ([99]Moye et al., 2019). HSBP1 is responsible for organising the muscle cytoskeleton. HSBP1 (also known as HSP27) is important for muscle formation. A mouse HSBP1 knockout showed destructured myofibrils and increased gaps between myofibrils ([100]Kammoun et al., 2016). RPRD2 is an RNA polymerase II interacting protein that co-purifies with RPRD1A, RPRD1B, and RPRD2 and contains serine and proline-rich regions ([101]Ni et al., 2011). The proline-rich domains function as docking sites for signalling protein modules and play important roles in protein-protein interactions ([102]Chin et al., 1997; [103]Elias et al., 2020). The overexpression of this gene might alter transcription, thus hampering protein interactions. The results of the network analysis showed downregulation of cGMP/PKG signalling, which is responsible for sperm capacitation through Ca^2+ and tyrosine phosphorylation in the presence of C-type natriuretic peptide (CNP) ([104]Wu et al., 2019). CNP, which is localised in the acrosomal region of the sperm head and tail, has a dose-dependent role in acrosome reaction and motility ([105]Xia et al., 2016). CNP might be affected by the process of cryopreservation, as fast freezing increases the ice-liquid interface with protein molecules, thereby increasing protein damage ([106]Cao et al., 2003), which in turn affects sperm motility and acrosome reaction. The MAPK signalling pathway is involved in sperm development and function through spermatogenesis and fertilising potential. In mature spermatozoa, the MAPK signalling pathway plays an important role in acrosome reaction, hyperactivation, and motility ([107]Almog and Naor, 2010; [108]Ebenezer Samuel King et al., 2022). We also observed the downregulation of cytoskeletal protein binding pathways in the cryopreserved sample. These pathways are linked to tubulin binding. Actin, along with the cytoskeletal protein tubulin, is involved in the regulation of sperm motility, capacitation, and acrosome reaction ([109]Salvolini et al., 2013). In normal spermatogenesis, oestrogen signalling through ESR1 is essential. Human spermatozoa contain functional aromatase, which is expressed even after ejaculation. Aromatase, in the presence of oestrogen receptor, has functions including sperm mobility and fertilising ability ([110]Carreau et al., 2011; [111]Cooke and Walker, 2022). The ubiquitin ligase complex, which is required for caspase regulation, was also downregulated in cryopreserved spermatozoa. Caspase is an important protease required for apoptosis. The activation of caspase during sperm differentiation is regulated by the ubiquitin ligase complex ([112]Arama et al., 2007). mRNA binding proteins are needed for the differentiation of spermatids into spermatozoa by modulating the expression of specific mRNAs ([113]Wishart and Dixon, 2002). All the pathways important for spermatozoa function, especially motility, capacitation, and acrosome reaction, were downregulated in the cryopreserved spermatozoa in the present study. Genes involved in the embryo organ morphogenesis pathway were upregulated in cryopreserved spermatozoa. This pathway is required for tissue and organ development to perform special functions during the embryonic stage. Among these, SPRY2 is required for the normal development of genitalia; TSHZ1 is required for the development of the soft palate, middle ear, and axial skeleton; and SATB2 codes for the jaw development ([114]Britanova et al., 2006; [115]Coré et al., 2007; [116]Ching et al., 2014). The overexpression of these genes might lead to abnormal embryo development. PRICKLE 1 is involved in ectodermal placode formation and also plays a role in nervous system development and nerve cell movement. A polymorphism in PRICKLE 1 reportedly affected acrosome integrity and DNA quality in cryopreserved spermatozoa ([117]Mańkowska et al., 2020). The regulation of protein autophosphorylation, one of the upregulated pathways in cryopreserved spermatozoa, is involved in the positive regulation of protein autophosphorylation. RAP2B in protein autophosphorylation codes for miR-205; miR-205 overexpression inhibits the PI3K/AKT signalling pathway ([118]Cui et al., 2020). Earlier studies also showed that cryopreservation suppressed the expression of important pathways and genes in sperm, causing DNA fragmentation and morphological deformation ([119]Hossen et al., 2021). Moreover, sperm stored in liquid nitrogen for longer times showed significantly increased abnormalities ([120]Malik et al., 2015). A study in zebrafish showed that cryopreservation was the major cause of genetic and epigenetic changes in germ cells ([121]Riesco and Robles, 2013). Our findings demonstrated the alteration of the transcriptional abundances of genes involved in important pathways in sperm function and fertility by the process of cryopreservation. These findings open new avenues for targeted studies on individual pathways and the development of ameliorative measures to minimise the effect of cryopreservation on the molecular health of sperm. Conclusion In conclusion, the transcriptional abundances of many genes were altered in cryopreserved spermatozoa, leading to changes in pathways important for sperm function and fertility. Earlier studies also indicated that the cryopreservation process altered the sperm expression levels of genes crucial for fertilisation and early embryo development ([122]Valcarce et al., 2013), consistent with our observations. The cryopreservation process induced alterations in the abundance of sperm transcripts and potential fertility-associated pathways, including glycogen metabolism, cGMP/PKG signalling, MAPK signalling, and the ubiquitin ligase complex, which could be a possible explanation for the reduced fertility observed in cryopreserved bull spermatozoa. Acknowledgments