Differential abundance of microRNAs in seminal plasma extracellular vesicles (EVs) in Sahiwal cattle bull related to male fertility

Sahiwal cattle, known for their high milk yield, are propagated through artificial insemination (AI) using male germplasm, largely contingent on semen quality. Spermatozoa, produced in the testes, carry genetic information and molecular signals essential for successful fertilization. Seminal plasma, in addition to sperm, contains nano-sized lipid-bound extracellular vesicles (SP-EVs) that carry key biomolecules, including fertility-related miRNAs, which are essential for bull fertility. The current study focused on miRNA profiling of SP-EVs from high-fertile (HF) and low-fertile (LF) Sahiwal bulls. SP-EVs were isolated using size exclusion chromatography (SEC) and characterized by dynamic light scattering (DLS) and nanoparticle tracking analysis (NTA). Western blotting detected the EV-specific protein markers TSG101 and CD63. The DLS analysis showed SP-EV sizes of 170–180 nm in HF and 130–140 nm in LF samples. The NTA revealed particle concentrations of 5.76 × 10¹⁰ to 5.86 × 10¹¹ particles/mL in HF and 5.31 × 10¹⁰ to 2.70 × 10¹¹ particles/mL in LF groups, with no significant differences in size and concentration between HF and LF. High-throughput miRNA sequencing identified 310 miRNAs in SP-EVs from both groups, with 61 upregulated and 119 downregulated in HF bull. Further analysis identified 41 miRNAs with significant fold changes and p-values, including bta-miR-1246, bta-miR-195, bta-miR-339b, and bta-miR-199b, which were analyzed for target gene prediction. Gene Ontology (GO) and KEGG pathway analyses indicated that these miRNAs target genes involved in transcription regulation, ubiquitin-dependent endoplasmic reticulum-associated degradation (ERAD) pathways, and signalling pathways. Functional exploration revealed that these genes play roles in spermatogenesis, motility, acrosome reactions, and inflammatory responses. qPCR analysis showed that bta-miR-195 had 80% higher expression in HF spermatozoa compared to LF, suggesting its association with fertility status (p < 0.05). In conclusion, this study elucidates the miRNA cargoes in SP-EVs as indicators of Sahiwal bull fertility, highlighting bta-miR-195 as a potential fertility factor among the various miRNAs identified.

Introduction

High genetic merit bulls have the capacity to propagate their progenies at a faster rate, as a single bull can be used to breed thousands of cows through artificial insemination (AI). The availability of good breeding bulls is a major bottleneck in the production of the requisite quantity of frozen straws. Limitation of producing sufficient number of frozen straws remains a concern in Sahiwal cattle farming. While both males and females contribute to the success or failure of pregnancy, the selection of highly fertile bulls holds greater significance in artificial breeding. Although bulls are selected based on breeding soundness evaluation (BSE) criteria, bull fertility remains diverse in bull farms (Tanga et al., 2021). A single ejaculate from a bull contains only 5% sperm and 95% seminal plasma. Seminal plasma is a blend of fluids produced by different glands and sexual organs, serving as a protective environment for sperm following ejaculation (Juyena and Stelletta, 2012; Druart and de Graaf, 2018). Seminal plasma plays a crucial role in sperm metabolism, function, survival, and transport. This indicates significant interactions and material exchanges occur between seminal plasma and sperm (Björndahl and Kvist, 2009; Bromfield, 2014). Additionally, seminal plasma is rich in various forms of RNA, including mRNA, microRNA, and tsRNAs (Li et al., 2012; Kuster and Althouse, 2016). The seminal plasma of bulls contains a large number of extracellular vesicles (SP-EVs), similar to that found in other animals (Talluri et al., 2017). The seminal plasma EVs are small, nano-sized (50–500 nm), lipid membrane-bound structures that are released mainly by the epididymis and the prostate (Saez et al., 2003; Aalberts et al., 2014). EVs transport different types of cargo like DNA, coding and non-coding RNA (miRNA), proteins, lipids, and carbohydrates (Takahashi et al., 2017) to the spermatozoa. EVs help in cell-cell communication and regulate sperm functions (Machtinger et al., 2016).

SP-EVs play a pivotal role in different sperm functions, such as sperm capacitation (Pons-Rejraji et al., 2011), sperm motility (Fabiani et al., 1995; Murdica et al., 2019), and the acrosome reaction (Piehl et al., 2013). Non-coding RNAs (nc-RNAs) in seminal plasma extracellular vesicles have been linked to the protection of sperm within the female genital tract (Aalberts et al., 2014) and also help in sperm-egg fusion (Franchi et al., 2020). The repertoire of non-coding RNAs (nc-RNAs), particularly microRNAs (miRNAs), and proteins in human semen extracellular vesicles (EVs) was analyzed using RNA and proteome sequencing. miRNAs are single-stranded non-coding RNA molecules, typically 18–25 base pairs in length, that regulate gene expression at the post-transcriptional level (Iacomino, 2023). miRNAs play a role in regulating numerous genes associated with reproductive processes, including the development and function of the reproductive tract, germ cell development and maturation, fertilization, and early embryonic development (Vojtech et al., 2014). The abundance of miRNAs in semen-derived extracellular vesicles (EVs) significantly increased from 7.5% in the caput epididymis to 31% in the cauda epididymis, suggesting that miRNAs in seminal plasma EVs played a crucial biological role (Chen et al., 2020). Using high-throughput sequencing, a study of the expression of miRNAs in the testes of chickens with low and high sperm motility revealed that 182 miRNAs were upregulated and 120 miRNAs were downregulated in the testes of chickens with high motility sperm (Xing et al., 2020). It also showed that has-miR-629-3p is linked to sperm motility in human seminal samples (Salas-Huetos et al., 2015), and that miRNA-122-5p plays a key role in sperm motility and the ability of sperm to fertilize an oocyte in bull sperm (da Silva et al., 2021). Furthermore, the expression of let-7a, let-7d, let-7e, and miR-22 was increased in porcine sperm with abnormal motility and morphology, while the expression of miR-15b was decreased (Curry et al., 2011). However, the mechanism by which miRNAs affect sperm motility and morphology has not yet been explored. SP-EVs are rich in bioactive molecules, including proteins and small non-coding RNAs (sncRNAs) such as miRNAs, Piwi-interacting RNAs (piRNAs), ribosomal RNAs (rRNAs), and transfer RNAs (tRNAs) (Lange-Consiglio et al., 2022). miR-31-5p has been identified as a biomarker for azoospermia (Lange-Consiglio et al., 2022), while two exosomal miRNAs, miR-196-5p and miR-501-3p, are known to be downregulated in prostate cancer.

The miRNA from SP-EVs of the various animals, including rats (Fornés et al., 1991), rabbits (Davis, 1978), bulls (Breitbart and Rubinstein, 1982) stallions (Arienti et al., 1998), and boars (Piehl et al., 2013) have been reported. However, the full miRNA profile of SP-EVs in Sahiwal cattle bulls has not yet been investigated, and it remains unclear how the dynamics of miRNA in SP-EVs compared to spermatozoa affect the fertilizing potential of the sperm. Given this background, in the current study we investigated the differentially expressed short non-coding miRNA profiles of SP-EVs in Sahiwal bulls of distinct fertility status. We report the detection of these specific miRNAs in both SP-EVs and spermatozoa, providing a comprehensive understanding of how the transfer of miRNAs from SP-EVs to spermatozoa affects fertility in Sahiwal bulls. The findings from this study could offer valuable insights into the role of miRNAs cargo in SP-EVs in bull fertility.

Materials and methods

Selection of distinct fertility bulls and collection of fresh semen

The standing bulls included in this study were sourced from the Artificial Breeding Research Centre (ABRC) at the National Dairy Research Institute (NDRI), India, based on their fertility status. Out of the 32 bulls assessed, six Sahiwal bulls with distinct fertility, each having over 50 records of artificial inseminations, were selected for further investigation. These bulls underwent evaluation for breeding soundness and semen quality parameters, including semen volume, sperm count, viability, and progressive motility. The chosen bulls were managed in accordance with the routine feeding and management plan of ABRC, NDRI.

The conception rates (CR) of the 32 bulls were examined for normal distribution. The CR demonstrated a normal distribution with a mean (M) of 44.485% and a standard deviation (S.D.) of 7.88%. Bulls with CRs below M-1 S.D. (36.6%) were classified as low fertile (LF), while those with CRs above M+1 S.D. (52.36%) were classified as high fertile (HF) (Supplementary Table S1).

For this study, six Sahiwal bulls (n = 6) were selected, with three falling into each fertility category (HF, n = 3, CR between 55.81% and 59.59%, and LF, n = 3, CR between 29.15% and 30.43%) shown in Supplementary Figure S1 (Karanwal et al., 2023; Batra et al., 2020; Prakash et al., 2021; Somashekar et al., 2017).

Isolation and processing of seminal plasma

The semen samples from HF (n = 3) and LF (n = 3) bulls were collected and transported to the lab at 37°C, ensuring no delay (within 30 min). They were then centrifuged at 1,520 g for 15 min. The resulting supernatant was transferred to a fresh tube and centrifuged again at 850 g for 5 min at 37°C. The final supernatant (seminal plasma) was either utilized for EV isolation or stored at −80°C (Badrhan et al., 2024).

Isolation of SP-EVs through size exclusion chromatography

SEC columns were assembled using Sepharose CL-2B (100 mL, Sigma Aldrich; New Delhi, India) and rinsed with 1X PBS containing 0.32% trisodium citrate (pH 7.4, 0.22 μm filtered). The column was packed with 10 mL of Sepharose 2B beads. Subsequently, 1 mL of seminal plasma was applied to the column, and elution was carried out using PBS + 0.32% trisodium citrate. The eluted material was collected in 26 distinct fractions of 0.5 mL each (Badrhan et al., 2024; Böing et al., 2014). Each fraction was then subjected to dynamic light scattering and nanoparticle tracking analysis to determine size and concentration, following the methods described by Böing et al. (2014). The composite of various SEC fractions of extracellular vesicles (EVs) was employed for subsequent experiments (Böing et al., 2014; Pal et al., 2023).

Characterization of SP-EVs

EVs were isolated from seminal plasma sample of three HF and LF bulls as previously described and pooled together. The isolated EVs were characterized through several approaches following MISEV 2018 guidelines.

Differential light scattering

Using a Zetasizer Nano ZS ZEN3600 (Malvern Instruments, Malvern, United Kingdom), the size distribution of EVs was determined. The scattered light’s intensity was measured at 173°. The particle size measurement was made in triplicate at 25°C, followed by data analysis and processing, which were carried out with the help of the Zetasizer software, version 7.03. After EV isolation, 10 μL of purified EVs were diluted in 990 μL of filtered 1X PBS and sonicated in an ultrasonic water bath for 1 min at 50 Hz. The sample was immediately put in a disposable cuvette for size measurements to avoid aggregation of the EVs. Three independent measurements were recorded for each sample (Pal et al., 2023; Szatanek et al., 2017).

Nanoparticle tracking assay

All particle tracking analyses used a Malvern NS300 device with a 488 nm laser and a 500 nm long-pass filter for fluorescence detection. To perform the NTA count, all samples were diluted in a ratio of 1:100 with 1X PBS (10010023, Gibco). The sample infusion pump was set to a constant flow rate of 5 μL/min. The camera level was set at 14, as all particles were visible at this level without signal saturation, and the detection threshold was set at 5. The ideal measurement concentration was found by pre-testing the ideal particles per frame value (40–120 particles/frame). The sample was measured five times, and 30- to 60-sec videos were collected with a minimum of 200 valid tracks recorded per video. Ten of the particles tracked were calculated and graphed using GraphPad Prism. All studies were performed on Nanosight 3.0 software using the default settings (Doyle and Wang, 2019; Pal et al., 2023).

Transmission electron microscopy

Isolated EVs were diluted 1:100 in 1X PBS and fixed with an equal volume of 2% paraformaldehyde prepared in phosphate buffer for 1 h at 4°C. The carbon Formvar film-coated 300-mesh transmission electron microscopy grids were glow-discharged. The fixed EV samples were applied to the TEM grid and incubated for 15 min at room temperature. After washing with PBS, the samples were fixed with 1% glutaraldehyde for 5 min. Upon washing with distilled water, the grids were stained with 1% phosphotungstic acid for 45 s, wicked off with Whatman filter paper, and allowed to dry before viewing. TEM examination was performed using the JEM1400 FLASH transmission electron microscope (JEOL USA Inc., Peabody, MA, United States) at 120 kV and viewed under ×250,000 magnification with a highly sensitive Olympus sCMOS camera at the TEM facility of the Advanced Technology Platform Centre (ATPC), Regional Centre for Biotechnology, Faridabad (Pal et al., 2023).

Western blotting

Protein was isolated by incubating the isolated EVs in RIPA buffer (Sigma-Aldrich) for 15 min at room temperature, followed by sonication. Approximately 20 µg of protein was loaded per well into SDS-PAGE gel, and proteins were separated by SDS polyacrylamide gel electrophoresis using a mini gel tank electrophoresis system (Invitrogen Life Technologies) and then transferred to Immobilon-FL polyvinylidene difluoride membranes (Millipore, Billerica, MA, United States). The membrane containing the transferred protein was probed with primary mouse polyclonal anti-TSG-101 (1:1000 SC-7964, Santa Cruz Biotechnology, United States), primary goat anti-CD63 (1:3000, STJ140029, St. Johns Laboratory, London, United Kingdom), primary CD9 monoclonal antibody (IVA50) (1:6000, MA1-19301, Invitrogen), and primary anti-Calnexin (CNX) monoclonal antibody (CAA280Hu22, Cloud Clone Corp, United States). Membranes were washed in Tris-buffered saline (TBS, pH 7.6) and incubated for 2 h in TBST (TBS with Tween-20) containing 5% BSA along with horseradish peroxidase-conjugated anti-mouse secondary antibody (1:10000, SC-516102, Santa Cruz Biotechnology), anti-goat secondary antibody (1:20000, STJ99512, St. Johns Laboratory), and anti-mouse secondary antibody (for CD9-1:30000, for Calnexin- 1:10000, A9044, Sigma), respectively. The Substrate Pierce ECL Western Blotting kit (32106) was utilized in the next step for chemiluminescence, and subsequently, the membranes were exposed to X-ray film for 1–5 min before visualization (Karanwal et al., 2023).

RNA extraction from SP-EVs

Isolated SP-EV fractions (7–16) were pooled and suspended in TRIzol (Ambion, Thermo Fisher Scientific, United States). Then the phase separation was performed by adding 200 µL of chloroform, followed by an invert mix and centrifugation at 13,200 rpm for 15 min at 4°C. The aqueous phase was transferred to fresh microcentrifuge tubes, and RNA was precipitated overnight at −80°C with an equal volume of isopropanol and 1 µL (20 μg/μL) glycogen. The pellet recovered by centrifugation was washed with 70% ethanol. The air-dried pellet was re-suspended in nuclease-free water. The concentration and purity of RNA were quantified using the Nanodrop Spectrophotometer (Thermo Scientific, 2000).

Library preparation

5 μL of the total RNA was taken for fragmentation and priming. Fragmented and primed RNA was further subjected to first-strand synthesis, followed by second-strand synthesis. The double-stranded cDNA was purified using NEBNext purification beads (NEBNext, Cat #E7767S). Purified cDNA was end-repaired, adenylated, and ligated to Illumina multiplex barcode adapters as per the NEBNext^® UltraTM II Directional RNA Library Prep protocol, followed by second-strand excision using the USER enzyme at 37°C for 15 min.

The Illumina Universal Adapters used in the study were:

5′AAT​GAT​ACG​GCG​ACC​ACC​GAG​ATC​TAC​ACT​CTT​TCC​CTA​CAC​GAC​GCT​CTT​CCG​ATC​T-3′ and Index Adapter 5′-GAT​CGG​AAG​AGC​ACA​CGT​CTG​AAC​TCC​AGT​CAC [INDEX] ATC​TCG​TAT​GCC​GTC​TTC​TGC​TTG-3' [INDEX] Unique sequence to identify sample-specific sequencing data. Adapter-ligated cDNA was purified using NEBNext purification beads and subjected to 12 cycles of indexing [98°C for 30 s, cycling (98°C for 10 s, 65°C for 75 s, and 65°C for 5 min)] to enrich the adapter-ligated fragments. The final PCR product (sequencing library) was purified with JetSeq beads, followed by a library quality control check (Prakash et al., 2021).

RNA sequencing and data analysis

The libraries were paired-end sequenced on an Illumina Novaseq 6000 sequencer for 150 cycles (Illumina, San Diego, United States) following the manufacturer’s instructions. Transcriptome analysis was done by processing the raw data for the removal of low-quality reads and adapter sequences. Expression analysis was performed using the high-quality reads after alignment with the reference genome [Bos indicus] using a splice-aware aligner. The raw reads were processed using FastQC1 for quality assessment and pre-processing, which includes removing adapter sequences and low-quality bases (< q30) using TrimGalore2. The pre-processed, high-quality data was aligned to the Bos indicus reference genome downloaded from the NCBI database using Hisat2 with default parameters. Reads were classified as aligned (which align to the reference genome) and unaligned reads. The feature count tool was used to estimate and calculate transcript abundance. Absolute counts for transcripts were estimated and used in differential expression analysis.

DESeq25 was used to calculate the differentially expressed transcripts. It performs an internal normalization where the geometric mean is calculated for each transcript across all samples. It fits negative binomial generalized linear models for each transcript and uses the Wald test for significance testing. Further, the expressed transcripts were categorized into up, down, and neutrally regulated based on the log2fold change cutoff of ±1 value.

miRNA target prediction

We employed pre-computed predictions from TargetScan (TS), miRmap (MM), microTCDS (mt), miRdb, and miRwalk, all of which are publicly accessible online. Target Scan, miRmap, takes into account the seed-based interactions in the target miRNA’s 3′UTR. Based on the miRTarget (MT) score determined by examining hundreds of miRNA target interactions from high-throughput sequencing research, targets in the miRDB are predicted. MicroTCDS recognizes miRNA targets in both the coding and 3′UTR regions (Van Rooij, 2011). These tools were chosen because their miRNA prediction capabilities were the most complete and they have an update strategy. In order to increase the sensitivity of target prediction for bta-miR-1246, bta-miR-195, and bta-miR-199b, we take into account the union of the minimum tools in the venn diagram of four tools for prediction. With the exception of bta-miR-339b, only three programs were used because this miRNA was not available in the database (Pal et al., 2023).

Gene ontology and pathway enrichment analysis

Gene ontology analyses were conducted on the chosen miRNA target genes and classified as biological processes (BP), molecular functions (MF), and cellular components (CC) through the Database for Annotation, Visualization, and Integrated Discovery (DAVID) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis. Significance was determined for pathways and GO terms with the highest number of genes, employing a p-value threshold of 0.05 (Nguyen et al., 2019; Fan et al., 2022).

Isolation of miRNA from sperm and cDNA synthesis

The miRNA was isolated from the sperm of HF (n = 3) and LF (n = 3) using the Qiagen Advanced miRNeasy Serum/plasma kit (cat. 217204, Qiagen, Germany) according to the manufacturer’s instructions, with the addition of bacteriophage MS2 RNA (10165948001, Roche) as a carrier RNA (to improve the RNA yield) in a 20 μL volume. Approximately 120 ng of the isolated RNA was converted into cDNA for each sample using the Reverse Transcriptase Kit II miScript (Cat# 218161 Qiagen, Valencia, CA) using HiSpec buffer (5X), miScript nucleic acid mix, and miScript enzyme (RT) in a final reaction volume of 20 µL for mature miRNA profiling. miRNA, cel-miR-39-3p, was used as an exogenous spike-in control. The RNA concentration for each sample was measured using a NanoDrop 1000 spectrophotometer (NanoDrop Technologies), and a 10–15 ng/μL concentration range was obtained from each sample. Approximately 120 ng of the isolated RNA was converted into cDNA for each sample using the Reverse Transcriptase Kit II miScript (Cat # 218161 Qiagen, Valencia, CA) using HiSpec buffer (5X), miScript nucleic acid mix, and miScript enzyme (RT) in a final reaction volume of 20 µL for mature miRNA profiling. The samples were incubated at 37°C for 1 h and then at 95°C for 5 min to inactivate the miScript enzyme (RT). Subsequently, the cDNA was diluted in nuclease-free water (2.5 ng/μL for qPCR) and stored at −20°C (Dubey et al., 2023; Pal et al., 2023).

Primer design for miRNA

Specific miRNAs were selected based on the NGS data that indicated their highest abundance in high-fertile cattle. A total of four miRNAs, i.e., bta-miR-195, bta-miR-1246, miR-195-5p, and miR-200a-3p, were selected, and miR-26a, let-7a, and miR-16a were initially chosen as reference controls. Primers were designed using the miR primer tool (Busk, 2014). The forward primers were made against the seed sequence of the miRNAs from the database available at http://www.mirbase.org/. The reference controls miR-23a, miR-26a, and U6 were later excluded due to their lesser stability in different biological samples as calculated by the online tool RefFinder (http://www.heartcure.com.au/refnder/). The mature miRNAs and the primer sequences utilized for the q-PCR assay have been listed in Supplementary Table S14. Primer oligos were procured from Sigma-Aldrich (Pal et al., 2023).

miRNA quantification using quantitative real-time PCR (RT-qPCR)

Relative quantification of selected miRNAs was performed using the miScript SYBR Green PCR Kit (QIAGEN, Valencia, CA, United States). Approximately 5 ng of cDNA was used as a template in each reaction. A no-template control (NTC) was also run on each plate. The reactions were performed in duplicate using the Bio-Rad c1000 thermal cycler detection system (United States) with the following conditions: 95°C for 5 min, 45 amplification cycles consisting of incubations at 94°C for 15 s, 55°C for 45 s, and 70°C for 30 s, followed by melt curve analysis from 65°C to 95°C. The relative expression of the target miRNAs was calculated using the 2−ΔΔCT method by considering let-7a and miR-16a as reference controls (Livak and Schmittgen, 2001). The differential expression levels of miRNAs between high-fertility and low-fertility bulls were analyzed statistically. The qPCR data was analyzed using GraphPad Prism 8.0. The datasets were analyzed by an unpaired t-test, and the difference was considered to be significant at P < 0.05 (Pal et al., 2023).

Results

Isolation and characterization of EVs from the seminal plasma of Sahiwal cattle bulls

SP-EVs were isolated from high fertile (HF) and low fertile (LF) bull seminal plasma using size exclusion chromatography (SEC). The isolated SP-EVs fractions were characterized using differential light scattering (DLS), nanoparticle tracking assay (NTA), transmission electron microscopy (TEM), and Western blot. The size measurements through DLS (Figure 1) revealed that the SP-EVs of HF Sahiwal bulls ranged from 179.3 to 202 nm in fractions 7-8, 177–190 nm in fractions 9–10, 160–174 nm in fractions 11–12, and 178–190 nm in fractions 13–14. For LF bulls, the size range was found to be from 160 nm to 186 nm in fractions 7–8, 150 nm–180 nm in fractions 9–10, 142 nm–160 nm in fractions 11–12, and 140 nm–155 nm in fractions 13–14, in accordance with Böing et al. (2014) (Supplementary Tables S2, S3). The NTA results revealed that the SP-EV population size ranged, with a mean average particle size of 152.2 ± 1.8 nm in 10–12 fractions of HF bulls and 146.4 ± 1.1 nm in LF bulls. The concentration of SP-EVs through NTA was found to be 7.19 × 10¹⁰ to 7.23 × 10¹⁰ particles/mL in HF bulls and 4.23 × 10¹⁰ to 6.85 × 10¹⁰ particles/mL in LF bulls (Supplementary Tables S4, S5). The characterization using DLS and NTA revealed that SP-EVs were predominantly enriched in fractions 7-14, leading to discard the other fractions. Based on NTA analysis, neither the size nor the concentration were found to be significantly different between the HF and LF bull SP-EVs (Figure 2). The general morphology and ultrastructure of seminal plasma-derived EVs were examined using transmission electron microscopy (TEM), which revealed round, cup-shaped vesicles with lipid bilayer structures, ranging in diameter from 120 to 200 nm (Figures 3A, B). The presence of small EVs was further confirmed through Western blot analysis using EV protein markers, including CD63 and TSG101 (Figures 3C, D). Full blot films are shown in Supplementary Figure S2.

Figure 1

Figure 1. Characterization of seminal EVs. Intensity-based size distribution of SP-EVs analyzed by Zetasizer nano zs particle sizer. Each graph shows mean ± SD (n = 3), PDI, and particle size distribution for (A) Fractions 9-10 (B) fractions 11-12 of SP- EVs of HF bulls, and (C) fractions 9-10 (D) fractions 11-12 of seminal EVs from LF bulls.

Figure 2

Figure 2. Characterization of seminal EVs Through Nanoparticle Tracking Analysis (NTA). The EVs were diluted 100X for NTA. The concentration vs. size graphs were plotted for (A) fractions 7-9, (B) 10-12 of HF bulls, and for fractions (C) 7-9, and (D) 10-12 of LF bulls.

Figure 3

Figure 3. TEM image of EVs derived from seminal plasma (A) HF bulls (B) LF bulls SP-EVs. EVs were negatively stained with 1% phosphotungstic acid after removing the extra moisture (Magnification-250,000×, Scale bar-50 nm, 120 kV). Identification of the CD63 (C) and TSG101 (D) EV-specific protein markers by the western blot analysis of isolated pooled EVs samples from the seminal plasma of Sahiwal bull.

miRNA profiling of SP-EVs of HF and LF bulls

DESeq normalized expression values were used to calculate the fold change for a given transcript based on the negative binomial distribution. In the case of HF samples, the transcripts that exhibited a log2-fold change less than −1 were represented as downregulated, while values greater than 1 were represented as upregulated, and values in between −1 and +1 were represented as neutrally regulated. A total of 310 miRNAs were fetched, out of which 61 were upregulated in HF, 119 were downregulated in HF, and 130 were neutrally regulated, i.e., there was no significant difference between the HF and LF groups (Figure 4A).

Figure 4

Figure 4. (A) Pie chart representation of the known miRNAs in SP-EVs. Total miRNAs identified in both HF and LF SP-EVs were 310 out of which 61 were upregulated in HF and 119 were downregulated in HF. (B) Bar graph showing significant upregulated and downregulated miRNAs in HF. The y-axis represents log fold change and the x-axis represents the shortlisted miRNAs in the decreasing order of their log fold changes.

Out of the entire pool of 310 miRNAs, a more refined selection was made by considering both the p-value and the log-fold change. On the basis of a p-value threshold of less than 0.05 and a log fold change exceeding 1, a subset of 41 miRNAs were identified, of which 20 were upregulated in HF and 21 were downregulated in HF. The miRNAs exhibited upregulation in HF bulls, as indicated in Table 1 and Figure 4B.

Table 1

Table 1. List of miRNAs that were upregulated in HF Sahiwal bulls.

In-silico target prediction of miRNA candidates

The target prediction of all the shortlisted 20 miRNAs that were upregulated in HF bulls was analyzed using online software such as Target Scan, miRwalk, miRmap, and Micro-T-CDS. Four miRNA candidates were selected, namely, bta-miR-339b, bta-miR-195, bta-miR1246, and bta-miR-199b. The clustering analysis of all predicted miRNA target genes is shown in Supplementary Tables S6–S9. Target genes identified by at least four databases were clustered using the Venn diagram (https://bioinformatics.psb.ugent.be/webtools/Venn/).

For bta-miR-1246, when considering all four tools, the number of targets predicted by the intersection of TargetScan, miRwalk, and miRmap was 14; the intersection of TargetScan, miRwalk, and micro-T-CDS yielded five genes; the intersection of TargetScan, miRmap, and micro-T-CDS exhibited an overlap of 26 genes; and the intersection of miRmap, miRwalk, and micro-T-CDS showed five genes. Additionally, when considering pairs of target prediction tools, TargetScan and miRwalk had 22 common genes, whereas 119 genes were common between TargetScan and miRmap. The overlap between TargetScan and micro-T-CDS was even more significant, with 143 common genes. Moreover, miRmap and miRwalk displayed an overlap of 33 genes, while miRwalk and micro-T-CDS had 18 genes in common. Finally, the tools micro-T-CDS and miRmap shared 47 genes (Figure 5A).

Figure 5

Figure 5. miRNA target prediction: (A) Venn diagram of the number of bta-miR-1246 targets predicted by each tool. (B) Venn diagram of the number of bta-miR-195 targets predicted by each tool. (C) Venn diagram of the number of bta-miR-199b targets predicted by each tool. (D) Venn diagram of the number of bta-miR-339b targets predicted by each tool.

In the case of bta-miR-195, one gene was common among all four tools. While considering the overlap of three softwares, four genes were found to be common in TargetScan, miRwalk, and micro-T-CDS, whereas 26 genes were found common in TargetScan, miRmap, and micro-T-CDS. An overlap of only one gene was observed between miRmap, miRwalk, and micro-T-CDS, and eight genes were identified as common in both TargetScan and miRwalk. A total of 63 genes were found common between TargetScan and miRmap; 198 genes overlapped in TargetScan and micro-T-CDS; and 28 genes were identified as common in miRmap and miRwalk. Among miRwalk and micro-T-CDS, 14 genes showed overlap, and lastly, 22 genes were found common between miRmap and micro-T-CDS (Figure 5B).

In the case of bta-miR-199b, none of the genes were found to be common among all four softwares. Two genes were common in TargetScan, miRmap, and miRwalk, while one gene was common in TargetScan, miRwalk, and micro-T-CDS. Fourteen genes were found to be common in TargetScan, miRmap, and micro-T-CDS. When considering two tools, our genes overlapped between TargetScan and miRwalk, while 38 genes were shared between TargetScan and miRmap. A total of 184 genes were found to be common in TargetScan and micro-T-CDS, and seven genes were found to overlap between miRmap and miRwalk. Among miRwalk and micro-T-CDS, twelve genes were found to be shared, while two genes were found to be common between miRmap and micro-T-CDS (Figure 5C).

In the case of bta-miR-339b, target prediction was done through miRmap, miRwalk, and micro-T-CDS (Figure 4). Due to the small number of gene candidates in the software Target Scan, it was not employed for this miRNA. Three genes were found to overlap among them, while 59 genes were common among miRmap and miRwalk. A total of 26 genes were found to be shared between miRwalk and micro-T-CDS (Figure 5D).

Gene ontology and functional pathway enrichment analysis for bta-miR-1246 targeted genes

For bta-miR-1246 target genes, the Gene Ontology (GO) term biological process revealed a significant involvement of genes in protein transport (Figure 6A). Notably, some of the identified genes, such as CFTR, ARL3b, and AXDND1 (p = 1.5E-4), have been involved in spermiogenesis and in the differentiation of spermatids into mature sperm Regarding the GO term analysis of molecular function for bta-miR-1246, a notable enrichment was observed for metal ion binding such as SMAD4, (Figure 6C). The GO term cellular components revealed a significant enrichment of genes in the nucleus and the cytoplasm such as GMCL1 (p = 6.92E-06), and PICK1 (p = 2.2E-6), has been found to be associated with spermatogenesis (Figure 6B).

Figure 6

Figure 6. Pie charts representing top Gene Ontology (GO) terms and pathway analysis for bta-miR-1246 with respect to (A) Biological Process (BP), (B) Cellular Component (CC) (C) Molecular Functions (MF) and (D) KEGG Pathway.

The KEGG pathway of bta-miR-1246 target genes such as FGF7 and GHR (p = 5.7E-4), revealed significant associations with the PI3-Akt pathway and the MAPK signaling pathway (Figure 6D; Supplementary Table S10).

GO and KEGG pathway enrichment analysis for bta-miR-195 targeted genes

The GO term biological processes in the case of bta-miR-195 target genes revealed a significant enrichment of these genes (PIM1, RAF, HMGA1, p = 2.3E-5) in positive regulation of transcription, DNA templated (p = 1.89E-04) (Figure 7A).

Figure 7

Figure 7. Pie charts representing top 10 Gene Ontology (GO) terms and pathway analysis for bta-miR-195 concerning (A) Biological Process (BP), (B) Cellular Component (CC) (C) Molecular Functions (MF) and (D) KEGG Pathway.

The GO term cellular components for bta-miR-195 target genes revealed that the maximum number of genes were present in the cytoplasm (ARL2, ARL3, BART, p = 1.4E-13) (Figure 7B). The GO term molecular functions revealed that the maximum number of genes were associated with several important functions, including (NDP, FBXL20 p = 1.7E-4), the sequence-specific SMAD binding, and protein kinase activity (Figure 7C).

The KEGG pathway enrichment analysis of bta-miR-195 target genes revealed that the maximum number of genes (LCRMP1, AKT3, p = 9.0E-6) were associated with the PI3K-Akt-signaling pathway (Figure 7D). Some of the identified genes of these pathways, such as. AKT3, play a crucial role in the proliferation of precursor cells in spermatogenesis, which needs tight regulation in order to maintain a normal sperm count without any abnormal spermatozoa (Wang et al., 2021). LCRMP 1 (p = 9.66E-06) is known to play crucial roles in the PI3-Akt and Ras signaling pathways. This gene is abundantly present in the testis, and upon its knockout in mice, it exhibits aberrant spermiation with apoptotic spermatids, suggesting its important role in spermatogenesis (Chang et al., 2023) (Supplementarty Table S11).

GO and KEGG pathway enrichment analysis for bta-miR-399b-targeted genes

GO-term biological processes for bta-miR-339b target genes revealed significant enrichment of these genes associated with the ubiquitin-dependent ERAD pathway (Figure 8A). The gene TSPAN1 (p 0.012549) and p53 (p = 0.03655) identified in GO biological process, which is known to be an EV biomarker, therefore suggesting a link between bta-miR-339b and EV production and The gene p53 (p = 0.03655) is involved in the apoptosis of abnormal sperm.

Figure 8

Figure 8. Pie charts representing top 10 Gene Ontology (GO) terms and pathway analysis for bta-miR-339b with respect to (A) Biological Process (BP), (B) Cellular Component (CC) (C) Molecular Functions (MF) and (D) KEGG Pathway.

The GO term molecular process for bta-miR-339b target genes revealed that a maximum number of genes were associated with protein homodimerization activity and indicated that bta-miR-339b is associated with specific genes, including CCNB2 (p = 0.008349) (Figure 8C).

Gene Ontology (GO) term cellular components showed enrichment of bta-miR-339b target genes within the nucleoplasm such as ALDP (p = 1.09E-04) and ZDHHC19 (p = 1.09E-04) (Figure 8B).

The functional analysis of bta-miR-339b target genes using KEGG pathway enrichment led to the identification of genes, i.e (p = 0.049008358), associated with the Ras signaling pathway (Figure 8D). The identification of the RABL2 target of bta-miR-339b highlights the potential regulatory role of this miRNA in signal transduction and during spermatogenesis (Supplementary Table S12).

Gene ontology and KEGG pathway enrichment analysis for bta-miR-199b-targeted genes

GO analysis of bta-miR-199b target genes highlights a significant enrichment of genes associated with crucial biological processes (Figure 9A), viz., regulation of transcription and signal transduction. Moreover, the functional analysis identified the genes SMAD2 and DAZ (p = 0.010716) as significant targets of bta-miR-199b. The GO term molecular function for bta-miR-199b revealed that the maximum number of genes were associated with RNA Polymerase II Core Promoter Proximal Region Sequence-Specific DNA Binding, GTPase Activity, and RNA Binding (Figure 9C), which play roles in transcriptional regulation as well as cell signaling in the cells.

Figure 9

Figure 9. Pie charts representing top 10 Gene Ontology (GO) terms and pathway analysis for bta-miR-199b concerning (A) Biological Process (BP), (B) Cellular Component (CC) (C) Molecular Functions (MF) and (D) KEGG Pathway.

The GO term cellular components for bta-miR-199b targets revealed that these are predominantly located in the nucleus Among the genes located in the nucleus, FXR1 was identified to be of particular significance (p = 6.78E-04) (Figure 9B).

The functional analysis of bta-miR-199b using the KEGG pathway analysis revealed that the Ras signaling pathway is enriched with the maximum number of genes associated with this miRNA. Within the Ras Signaling Pathway (Figure 9D), specific genes, such as cKIT and RB1 (p = 0.035763), were identified as significant targets of bta-miR-199b.

Relative expression analysis of selected miRNAs between HF and LF bull spermatozoa

Semen samples were collected from HF and LF bulls, and miRNA was isolated from the sperm cells using the miRneasy kit. RT-qPCR was performed to measure the relative expression levels of specific miRNAs in HF and LF bulls. The expression of bta-miR-195 was found to be significantly 80% higher (p-value <0.05) in HF bulls compared to LF bulls (Figure 10A). This indicates that bta-miR-195 is strongly associated with the uptake of bta-miR-195 on sperm, may play a role in regulating processes related to fertility, and could be associated with improved reproductive outcomes in HF bulls. However, no significant differences in expression were observed for bta-miR-1246 and bta-miR-339b between HF and LF bulls (Figures 10B, C). Though both miRNAs were differentially different in both HF and LF bulls. This suggests that these miRNAs may not be directly linked to fertility-related processes in the context of Sahiwal bulls. Unfortunately, due to a primer error, the amplification of bta-miR-199b could not be achieved and, hence, was excluded from this study.

Figure 10

Figure 10. Relative expression of (A) bta-miR-195, (B) bta-miR-1246, and (C) bta-miR-339b between the spermatozoa of HF and LF bulls. Expression values were normalized to let-7a miRNA. The -axis represents the bull fertility status and the y-axis represents the fold change. Vertical bars represent means of fold change and error bars represent the standard error of mean at p = 0.05.

Mounting evidence from this study suggests that the high expression of bta-miR-195 appears to have a potential association with the fertility status of Sahiwal bulls. It can be supported that the gene ontology and KEGG pathway analysis previously discussed have shown that overexpression of certain bta-miR-195 target genes, such as PIM1, AKT3, and FBXL20, could have detrimental effects on spermatozoa and the spermatogenesis process. As a result, it is hypothesized that bta-miR-195 might positively impact fertility in HF bulls by down-regulating the above genes.

Discussion

The unpredictability and limitations in current semen evaluation techniques in the assessment of bull fertility more precisely led us to explore the SP-EVs transcript-based alternative methods for predicting Sahiwal bull fertility. In this study, we hypothesize that the presence of miRNAs in SP-EVs plays a crucial role in regulating the fertility of Sahiwal bulls. In addition, the differential expression of miRNAs in SP-EVs and their role in fertility remain unexplained in indigenous cattle. So, this investigation was specifically designed to determine how the variations in the abundance of specific miRNAs in SP-EVs are related with Sahiwal bull fertility.

Herein, SP-EVs were isolated using SEC and characterized by DLS, NTA, TEM, and Western blotting. Fractions 7–14 contained EVs with a mean size of 200 nm in both fertility groups, consistent with findings from previous studies (Barranco et al., 2023). The particles in the sample were homogeneous, as indicated by a polydispersity index (PDI) ranging from 0.2 to 0.35 in the DLS data from EVs of both fertility groups. Additionally, the NTA assay not only confirmed the size of the seminal EVs but also verified their abundance in cattle seminal plasma. TEM also revealed the characteristic round, cup-shaped morphology of seminal EVs in both HF and LF bulls, within the desired size range (Xu et al., 2020). miRNA profiling of SP-EVs from high-fertility (HF) and low-fertility (LF) Sahiwal bulls was performed using high-throughput RNA sequencing. The present study identified a pool of 310 miRNAs. By refining the selection based on a p-value (< 0.05) and a log-fold change greater than one, a subset of 41 miRNAs were identified. Of these, 20 miRNAs were upregulated in high-fertility (HF) bulls and 21 were downregulated in HF bulls (Table 1), respectively. Our results are in concordance with the earlier studies by Dlamini et al. (2023).

The microRNAs have the ability to control target genes at the post-transcriptional level by either repressing their translation or causing the degradation of nucleic acids, thereby inhibiting gene expression (O’Brien et al., 2018). In accordance with previous literature, certain X-linked miRNAs are believed to play a crucial role in spermatogenesis and male reproductive processes (Li et al., 2014; Qing et al., 2017). Among the most abundant miRNAs, several were common across datasets, such as bta-miR-195 (Do et al., 2019), bta-miR-1246 (Salas-Huetos and Aston, 2021), bta-miR-486 (Donnellan et al., 2022), and bta-miR-199b (Kasimanickam et al., 2022). These miRNAs are anticipated to play crucial regulatory roles in diverse physiological processes within sperm. Additionally, gene targets for the identified miRNAs were predicted using various online tools (Pal et al., 2023). We selected four miRNAs—bta-miR-1246, bta-miR-195, bta-miR-339b, and bta-miR-199b—because they had the highest number of predicted target genes. As it is well known that, extracellular vesicles (EVs) play a crucial role in transmitting signals involved in cellular functions, including the maturation of germ cells, and EVs can transfer a variety of RNA molecules to sperm (Ding et al., 2021), Next we envisaged to uncover critical insights into the influence of miRNAs on sperm function and fertility by conducting gene ontology and KEGG pathway analyses.

The GO term biological process for bta-miR-1246 indicated a significant involvement of genes such as ARL3b and AXDND1 both of which are recognized for their critical roles in spermiogenesis and the differentiation of spermatids into mature sperm. These genes are particularly integral to protein transport pathways, suggesting that bta-miR-1246 may play a regulatory role in these essential reproductive processes (Ma et al., 2021; Qi et al., 2013). GO term analysis of bta-miR-1246 also revealed significant enrichment for metal ion binding. SMAD4, a metal ion-binding protein, mediates TGF-β signaling and is crucial for testicular development. Overexpression of SMAD4 has been linked to spermatogenic arrest and seminiferous tubule degradation (Narula et al., 2002). Our findings suggest a positive correlation, as bta-miR-1246 is upregulated in high fertility (HF), potentially regulating SMAD4 expression to facilitate normal spermatogenesis. GO-term analysis of cellular components revealed significant enrichment of genes in the nucleus and cytoplasm, including GMCL1, which is essential for nuclear envelope formation and spermatogenesis. GMCL1 knockout leads to abnormal or absent sperm production (Azhar et al., 2021). Another gene, PICK1 (p = 2.2E-6), is associated with spermatogenesis, and its downregulation has been linked to disruptions in the acrosome reaction (Du et al., 2023).

The pathway analysis of target genes regulated by bta-miR-1246 revealed significant associations with the PI3-Akt pathway and the MAPK signaling pathway. Genes involved, FGF and GHR, play roles in testicular maturation and contribute to germ cell production by releasing fibroblast growth factors and testosterone, crucial for testicle development. These genes also support sperm motility improvement (Garbarino Azúa et al., 2017; Jiang et al., 2013). Therefore, bta-miR-1246 miRNA is also implicated in the regulation of sperm motility Biological process of bta-miR-195 target genes showed significant enrichment in the positive regulation of transcription, DNA-templated processes. Genes such as HMGA1 and PIM1 are involved in spermatocyte-to-spermatid differentiation, with PIM1 upregulation potentially causing harmful inflammation to spermatozoa (Jiménez-Garciá et al., 2016). In cellular components for bta-miR-195 target genes, the majority of genes were found in the cytoplasm. ARL2 and BART are typically cytosolic but can translocate to mitochondria, where they bind to adenine nucleotide transporters, potentially influencing energy metabolism during spermatogenesis (Sharer et al., 2002). In the molecular function, a significant number of genes were associated with various important functions, including RNA polymerase activity, sequence-specific DNA binding, and protein homodimerization. Among these genes, NDPencodes nucleoside diphosphate (NDP) kinases that are crucial during spermiogenesis. Another identified gene, FBXL20, plays a role in inflammatory responses which can be detrimental to spermatozoa and ultimately lead to infertility (Yang et al., 2019). The pathway enrichment analysis of bta-miR-195 target genes revealed that the maximum number of genes were associated with the Ras signaling pathway.

The gene TSPAN1, identified in the GO biological process, encodes a tetraspanin protein known to serve as an extracellular vesicle (EV) biomarker, implying a potential association between bta-miR-339b and EV production (Jankovičová et al., 2020). The enrichment of TSPAN1 and p53 indicates that bta-miR-339b likely plays a critical role in regulating this pathway, which is essential for maintaining protein homeostasis and cellular integrity. In the molecular processes for bta-miR-339b target genes, a significant number of genes were associated with protein homodimerization activity. Specifically, bta-miR-339b was found to regulate genes such as CCNB2 (p = 8.5E-3), a member of the cyclin family critical for G2/M cell cycle progression through activation of CDK1. CCNB2 is expressed in testicular germ cells and is essential for spermatogenesis (Shen et al., 2022).

The cellular components revealed the enrichment of bta-miR-339b target genes within the nucleoplasm. Of particular interest, specific genes such as ALDP (p = 1.09E-04) are involved in the degradation process of very long-chain fatty acids during spermatogenesis (Linn et al., 2021). ZDHHC19 (p = 1.09E-04) is another gene identified as crucial for normal sperm function; its absence in spermatozoa leads to multiple defects, including abnormal sperm tail morphology, reduced motility, and disrupted acrosome reactions (Wang et al., 2021). The pathway functional analysis of bta-miR-339b target genes led to the identification of genes associated with the Ras signaling pathway. In the Ras signaling pathway, RABL2 was identified as significantly associated, playing a crucial role in signal transduction as a key mediator during spermatogenesis (Lo et al., 2012).

Furthermore, In the biological processes of bta-miR-199b target genes (SMAD2 and DAZ) are involved in transcription regulation and signal transduction. SMAD2 is present in testis, sertoli cells, and preleptotene to pachytene spermatocytes but not in post-meiotic germ cells, suggesting its involvement in spermatocyte differentiation (Wang and Zhao, 1999). DAZ, a highly conserved RNA-binding protein, is essential for gametogenesis in metazoans. It plays a role in regulating the translation of specific mRNAs, and its absence halts the process of spermatogenesis (Reynolds and Cooke, 2005). The GO analysis of cellular components for bta-miR-199b targets showed predominant localization in the nucleus. Notably, FXR1 emerged as significant among these nuclear genes. Downregulation or knockout of FXR1 has been associated with disrupted spermiogenesis, which is implicated in infertility (Kang et al., 2022).

In addition, the pathway functional analysis revealed that the Ras signaling pathway is enriched with the maximum number of genes associated with this miRNA. Specific genes, such as cKIT and RB1, were identified as significant targets of bta-miR-199b. The cKIT regulates primordial germ cell migration, proliferation, and apoptosis during fetal gonadal development (Mauduit et al., 1999). RB1 (Retinoblastoma 1) is a tumor suppressor gene known to regulate cell cycle progression. In the context of spermatogenesis, RB1 plays a critical role in controlling meiotic cell division, ensuring the proper development of haploid spermatocytes and oocytes, and maintaining the spermatogonial stem cell pool (Yang et al., 2019). The gene ontology and pathway analysis of the target genes for the four selected miRNAs, which are abundant in SP-EVs, indicate that many of these genes are associated with functions that either positively or negatively impact spermatogenesis and reproductive functions.

Following the functional analysis of various miRNAs associated with fertility, we proceeded to compare the relative expression of selected miRNAs between the spermatozoa of high-fertility (HF) and low-fertility (LF) bulls. The significant increase in the detection of bta-miR-195 in HF spermatozoa indicates that bta-miR-195 is closely associated with sperm uptake and may play a regulatory role in fertility-related processes, potentially contributing to improved reproductive outcomes in HF bulls. However, no significant differences in expression were observed for bta-miR-1246 and bta-miR-339b between HF and LF bulls, despite their differential expression patterns in both groups, highlighting that in HF spermatozoa these miRNAs abundance did not perfectly relate with bull fertility.

The findings of this study suggest a potential link between the high expression of bta-miR-195 indicating its capacity to regulate the target genes, such as PIM1, AKT3, and FBXL20 may negatively impact spermatozoa and the spermatogenesis process in highly fertile bulls with superior conception rates.

Additionally, due to logistical and financial limitations, we were able to collect semen samples from the bulls only twice a week, as the study relied entirely on standing bulls. This limited the availability of samples. For the purposes of characterization in the western blot experiment (for protein isolation) and sequencing (RNA isolation), we pooled fractions of SP-EVs. Pooling samples reduces the sample size while maintaining a high degree of confidence in the data (Somashekar et al., 2017). It increases the likelihood of detecting the maximum number of miRNAs with minimal inter-individual variation within a specific treatment group. Further experiments involving data analysis with a higher number of bulls within a specific (HF and LF) group will certainly help in gaining deeper and more robust insights into the findings of the current study.

Conclusion

In this study, we developed the miRNA profile of SP-EVs in contrasting fertility Sahiwal bulls using high-throughput RNA sequencing. Our findings suggest that the fertility status of Sahiwal bulls is not influenced by the size and concentration of SP-EVs. However, the miRNA profile of SP-EVs reveals numerous differentially abundant miRNAs. Based on target prediction of these miRNAs, followed by gene ontology and pathway analysis, indicate that the target of bta-miR-195 is associated with various sperm functions, such as sperm motility and acrosome reaction. Additionally, while examining the expression of these miRNAs in spermatozoa, bta-miR-195 exhibited approximately 80% higher expression in high-fertility bulls compared to the low-fertility group (p < 0.05), Highlighting a potential association between bta-miR-195 and the fertility status of Sahiwal bulls. Taken together leads clearly suggest that bta-miR-195 promises as a potential fertility indicator in Sahiwal bulls. Further validation of differentially abundant miRNAs as well as assessment of their transcripts in distinct fertility bulls could significantly enhance our understanding of miRNA in male fertility in Sahiwal cattle bulls.

Data availability statement

The datasets presented in this study are included within the manuscript and its Supplementary Material. Raw data supporting the findings of this study are openly available in an online repository. The name(s) of the repository/repositories and the corresponding details can be found below: https://www.ncbi.nlm.nih.gov/sra/PRJNA1147473, accession number PRJNA1147473.

Ethics statement

The Institutional Animal Ethics Committee (IAEC) at the National Dairy Research Institute in Karnal reviewed and approved the animal study. The study was conducted in accordance with the local legislation and institutional requirements.

Author contributions

VC: Conceptualization, Formal Analysis, Investigation, Methodology, Validation, Visualization, Writing–original draft, Writing–review and editing. PK: Conceptualization, Investigation, Methodology, Visualization, Writing–original draft, Writing–review and editing. JC: Data curation, Formal Analysis, Visualization, Writing–review and editing, Writing–original draft. AnP: Formal Analysis, Methodology, Validation, Writing–review and editing, Writing–original draft. AdP: Methodology, Validation, Writing–review and editing, Writing–original draft. SK: Methodology, Visualization, Writing–review and editing, Writing–original draft. SB: Methodology, Visualization, Writing–review and editing, Writing–original draft. FJ: Methodology, Writing–review and editing, Writing–original draft. SS: Conceptualization, Methodology, Writing–review and editing, Writing–original draft. MB: Resources, Supervision, Writing–review and editing, Writing–original draft. TD: Conceptualization, Funding acquisition, Methodology, Project administration, Writing–review and editing, Writing–original draft. RK: Conceptualization, Data curation, Funding acquisition, Methodology, Project administration, Resources, Supervision, Visualization, Writing–original draft, Writing–review and editing, Formal Analysis, Investigation, Validation.

Funding

The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. The experiments in this research were funded by the Department of Biotechnology, New Delhi, India (Grant No. BT/PR34220/AAQ/1/771/2020) and the Bill & Melinda Gates Foundation (Grant No. INV-008501).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcell.2024.1473825/full#supplementary-material

References

Aalberts, M., Stout, T. A. E., and Stoorvogel, W. (2014). Prostasomes: extracellular vesicles from the prostate. Reproduction 147, 1–14. doi:10.1530/REP-13-0358

PubMed Abstract | CrossRef Full Text | Google Scholar

Arienti, G., Carlini, E., De Cosmo, A. M., Di Profio, P., and Palmerini, C. A. (1998). Prostasome-like particles in stallion semen. Biol. Reprod. 59, 309–313. doi:10.1095/biolreprod59.2.309

PubMed Abstract | CrossRef Full Text | Google Scholar

Azhar, M., Altaf, S., Uddin, I., Cheng, J., Wu, L., Tong, X., et al. (2021). Towards post-meiotic sperm production: genetic insight into human infertility from mouse models. Int. J. Biol. Sci. 17, 2487–2503. doi:10.7150/ijbs.60384

PubMed Abstract | CrossRef Full Text | Google Scholar

Badrhan, S., Karanwal, S., Pal, A., Chera, J. S., Chauhan, V., Patel, A., et al. (2024). Differential protein repertoires related to sperm function identified in extracellular vesicles (EVs) in seminal plasma of distinct fertility buffalo (Bubalus bubalis) bulls. Front. Cell Dev. Biol. 12 (12), 1400323. doi:10.3389/fcell.2024.1400323

PubMed Abstract | CrossRef Full Text | Google Scholar

Barranco, I., Sanchez-López, C. M., Bucci, D., Alvarez-Barrientos, A., Rodriguez-Martinez, H., Marcilla, A., et al. (2023). The proteome of large or small extracellular vesicles in pig seminal plasma differs, defining sources and biological functions. Mol. Cell Proteomics 22, 100514. doi:10.1016/j.mcpro.2023.100514

PubMed Abstract | CrossRef Full Text | Google Scholar

Batra, V., Dagar, K., Nayak, S., Kumaresan, A., Kumar, R., and Datta, T. K. (2020). A higher abundance of O-linked glycans confers a selective advantage to high fertile Buffalo spermatozoa for immune-evasion from neutrophils. Front. Immunol. 11, 1928. Published 2020 Aug 28. doi:10.3389/fimmu.2020.01928

PubMed Abstract | CrossRef Full Text | Google Scholar

Björndahl, L., and Kvist, U. (2009). Human sperm chromatin stabilization: a proposed model including zinc bridges. Mol. Hum. Reprod. 16, 23–29. doi:10.1093/molehr/gap099

PubMed Abstract | CrossRef Full Text | Google Scholar

Böing, A. N., van der Pol, E., Grootemaat, A. E., Coumans, F. A. W., Sturk, A., and Nieuwland, R. (2014). Single-step isolation of extracellular vesicles by size-exclusion chromatography. J. Extracell. Vesicles 3, 23430. doi:10.3402/jev.v3.23430

PubMed Abstract | CrossRef Full Text | Google Scholar

Breitbart, H., and Rubinstein, S. (1982). Characterization of Mg2+- and Ca2+-ATPase activity in membrane vesicles from ejaculated ram seminal plasma. Arch. Androl. 9, 147–157. doi:10.3109/01485018208990233

PubMed Abstract | CrossRef Full Text | Google Scholar

Bromfield, J. J. (2014). Seminal fluid and reproduction: much more than previously thought. J. Assist. Reprod. Genet. 31, 31627–31636. doi:10.1007/s10815-014-0243-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Busk, P. K. (2014). A tool for design of primers for microRNA-specific quantitative RT-qPCR. BMC bioinform. 15, 1–9. doi:10.1186/1471-2105-15-29

PubMed Abstract | CrossRef Full Text | Google Scholar

Chang, J. H., Chou, C. H., Wu, J. C., Liao, K. M., Luo, W. J., Hsu, W. L., et al. (2023). LCRMP-1 is required for spermatogenesis and stabilises spermatid F-actin organization via the PI3K-Akt pathway. Commun. Biol. 6, 389. doi:10.1038/s42003-023-04778-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, X., Zheng, Y., Lei, A., Zhang, H., Niu, H., Li, X., et al. (2020). Early cleavage of preimplantation embryos is regulated by tRNAGln-TTG-derived small RNAs present in mature spermatozoa. J. Biol. Chem. 295 (32), 10885–10900. doi:10.1074/jbc.RA120.013003

PubMed Abstract | CrossRef Full Text | Google Scholar

Curry, E., Safranski, T. J., and Pratt, S. L. (2011). Differential expression of porcine sperm microRNAs and their association with sperm morphology and motility. Theriogenology 76 (8), 1532–1539. doi:10.1016/j.theriogenology.2011.06.025

PubMed Abstract | CrossRef Full Text | Google Scholar

Da Silva, L. D. S., Domingues, W. B., Barreto, B. F., Da Silveira, M., Dellagostin, E. N., Komninou, E. R., et al. (2021). Capillary electroporation affects the expression of miRNA-122-5p from bull sperm cells. Gene 768, 145286. doi:10.1016/j.gene.2020.145286

PubMed Abstract | CrossRef Full Text | Google Scholar

Davis, B. K. (1978). Uterine fluid from progesterone treated rabbits contains subcellular membranes. Experientia 34, 350–351. doi:10.1007/BF01923028

PubMed Abstract | CrossRef Full Text | Google Scholar

Ding, Y., Ding, N., Zhang, Y., Xie, S., Huang, M., Ding, X., et al. (2021). MicroRNA-222 transferred from semen extracellular vesicles inhibits sperm apoptosis by targeting BCL2L11. Front. Cell Dev. Biol. 9, 736864. doi:10.3389/fcell.2021.736864

PubMed Abstract | CrossRef Full Text | Google Scholar

Dlamini, N. H., Nguyen, T., Gad, A., Tesfaye, D., Liao, S. F., Willard, S. T., et al. (2023). Characterization of extracellular vesicle-coupled miRNA profiles in seminal plasma of boars with divergent semen quality status. Int. J. Mol. Sci. 24, 3194. doi:10.3390/ijms24043194

PubMed Abstract | CrossRef Full Text | Google Scholar

Do, D. N., Dudemaine, P. L., Fomenky, B. E., and Ibeagha-Awemu, E. M. (2019). Integration of miRNA weighted gene co-expression network and miRNA-mRNA co-expression analyses reveals potential regulatory functions of miRNAs in calf rumen development. Genomics 111, 849–859. doi:10.1016/j.ygeno.2018.05.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Donnellan, E. M., Perrier, J. P., Keogh, K., Štiavnická, M., Collins, C. M., Dunleavy, E. M., et al. (2022). Identification of differentially expressed mRNAs and miRNAs in spermatozoa of bulls of varying fertility. Front. Vet. Sci. 9, 993561. doi:10.3389/fvets.2022.993561

PubMed Abstract | CrossRef Full Text | Google Scholar

Doyle, L. M., and Wang, M. Z. (2019). Overview of extracellular vesicles, their origin, composition, purpose, and methods for exosome isolation and analysis. Cells 8, 727. doi:10.3390/cells8070727

PubMed Abstract | CrossRef Full Text | Google Scholar

Druart, X., and de Graaf, S. (2018). Seminal plasma proteomes and sperm fertility. Anim. Reprod. Sci. 194, 33–40. doi:10.1016/j.anireprosci.2018.04.061

PubMed Abstract | CrossRef Full Text | Google Scholar

Du, Y. Q., Shu, C. Y., Zheng, M., Xu, W., Sun, Y., Shen, L., et al. (2023). Truncating PICK1 variant identified in azoospermia affected mitochondrial dysfunction in knockout mice. Curr. Med. Sci. 43, 313–323. doi:10.1007/s11596-023-2704-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Dubey, P., Batra, V., Sarwalia, P., Nayak, S., Baithalu, R., Kumar, R., et al. (2023). miR-1246 is implicated as a possible candidate for endometrium remodelling facilitating implantation in buffalo (Bubalus bubalis). Vet. Med. Sci. 9, 443–456. doi:10.1002/vms3.968

PubMed Abstract | CrossRef Full Text | Google Scholar

Fabiani, R., Johansson, L., Lundkvist, Ö., and Ronquist, G. (1995). Prolongation and improvement of prostasome promotive effect on sperm forward motility. Eur. J. Obstet. Gynecol. Reprod. Biol. 58, 191–198. doi:10.1016/0028-2243(94)01997-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Fan, N., Yuan, S., Hai, Y., Du, P., Li, J., Kong, X., et al. (2022). Identifying the potential role of IL-1β in the molecular mechanisms of disc degeneration using gene expression profiling and bioinformatics analysis. J. Orthop. Surg. 30, 23094990211068203. doi:10.1177/23094990211068203

CrossRef Full Text | Google Scholar

Fornés, M. W., Barbieri, A., Sosa, M. A., and Bertini, F. (1991). First observations on enzymatic activity and protein content of vesicles separated from rat epididymal fluid. Andrologia 23, 347–351. doi:10.1111/j.1439-0272.1991.tb02578.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Franchi, A., Moreno-Irusta, A., Domínguez, E. M., Adre, A. J., and Giojalas, L. C. (2020). Extracellular vesicles from oviductal isthmus and ampulla stimulate the induced acrosome reaction and signaling events associated with capacitation in bovine spermatozoa. J. Cell Biochem. 121, 2877–2888. doi:10.1002/jcb.29522

PubMed Abstract | CrossRef Full Text | Google Scholar

Garbarino Azúa, D. J., Saucedo, L., Giordana, S., Magri, M. L., Buffone, M. G., Neuspiller, F., et al. (2017). Fibroblast growth factor 2 (FGF2) is present in human spermatozoa and is related with sperm motility. The use of recombinant FGF2 to improve motile sperm recovery. Andrology 5, 990–998. doi:10.1111/andr.12398

PubMed Abstract | CrossRef Full Text | Google Scholar

Iacomino, G. (2023). miRNAs: the road from bench to bedside. Genes (Basel) 14, 314. doi:10.3390/genes14020314

PubMed Abstract | CrossRef Full Text | Google Scholar

Jankovičová, J., Sečová, P., Michalková, K., and Antalíková, J. (2020). Tetraspanins, more than markers of extracellular vesicles in reproduction. Int. J. Mol. Sci. 21, 7568. doi:10.3390/ijms21207568

PubMed Abstract | CrossRef Full Text | Google Scholar

Jiang, X., Skibba, M., Zhang, C., Tan, Y., Xin, Y., and Qu, Y. (2013). The roles of fibroblast growth factors in the testicular development and tumor. J. Diabetes Res. 2013, 489095. doi:10.1155/2013/489095

PubMed Abstract | CrossRef Full Text | Google Scholar

Jiménez-Garciá, M. P., Lucena-Cacace, A., Robles-Friás, M. J., Narlik-Grassow, M., Blanco-Aparicio, C., and Carnero, A. (2016). The role of PIM1/PIM2 kinases in tumors of the male reproductive system. Sci. Rep. 6, 38079. doi:10.1038/srep38079

PubMed Abstract | CrossRef Full Text | Google Scholar

Juyena, N. S., and Stelletta, C. (2012). Seminal plasma: an essential attribute to spermatozoa. J. Androl. 33, 536–551. doi:10.2164/jandrol.110.012583

PubMed Abstract | CrossRef Full Text | Google Scholar

Kang, J. Y., Wen, Z., Pan, D., Zhang, Y., Li, Q., Zhong, A., et al. (2022). LLPS of FXR1 drives spermiogenesis by activating translation of stored mRNAs. Science 377, eabj6647. doi:10.1126/science.abj6647

PubMed Abstract | CrossRef Full Text | Google Scholar

Karanwal, S., Pal, A., Chera, J. S., Batra, V., Kumaresan, A., Datta, T. K., et al. (2023). Identification of protein candidates in spermatozoa of water buffalo (Bubalus bubalis) bulls helps in predicting their fertility status. Front. Cell Dev. Biol. 11, 111920. doi:10.3389/fcell.2023.1119220

CrossRef Full Text | Google Scholar

Kasimanickam, V., Kumar, N., and Kasimanickam, R. (2022). Investigation of sperm and seminal plasma candidate MicroRNAs of bulls with differing fertility and in silico prediction of miRNA-mRNA interaction network of reproductive function. Animals 12, 2360. doi:10.3390/ani12182360

PubMed Abstract | CrossRef Full Text | Google Scholar

Kuster, C. E., and Althouse, G. C. (2016). The impact of bacteriospermia on boar sperm storage and reproductive performance. Theriogenology 85, 21–26. doi:10.1016/j.theriogenology.2015.09.049

PubMed Abstract | CrossRef Full Text | Google Scholar

Lange-Consiglio, A., Capra, E., Monferini, N., Canesi, S., Bosi, G., Cretich, M., et al. (2022). Extracellular vesicles from seminal plasma to improve fertilizing capacity of bulls. Reprod. Fertil. 3, 313–327. doi:10.1530/raf-22-0037

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, H., Huang, S., Guo, C., Guan, H., and Xiong, C. (2012). Cell-free seminal mRNA and microRNA exist in different forms. PLoS One 7, 34566. doi:10.1371/journal.pone.0034566

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, J., Ma, W., Zeng, P., Wang, J., Geng, B., Yang, J., et al. (2014). LncTar: a tool for predicting the RNA targets of long noncoding RNAs. Brief. Bioinform 16, 806–812. doi:10.1093/bib/bbu048

PubMed Abstract | CrossRef Full Text | Google Scholar

Linn, E., Ghanem, L., Bhakta, H., Greer, C., and Avella, M. (2021). Genes regulating spermatogenesis and sperm function associated with rare disorders. Front. Cell Dev. Biol. 9, 6434536. doi:10.3389/fcell.2021.634536

PubMed Abstract | CrossRef Full Text | Google Scholar

Livak, K. J., and Schmittgen, T. D. (2001). Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods 25, 402–408. doi:10.1006/meth.2001.1262

PubMed Abstract | CrossRef Full Text | Google Scholar

Lo, J. C., Jamsai, D., O'Connor, A. E., Borg, C., Clark, B. J., Whisstock, J. C., et al. (2012). RAB-like 2 has an essential role in male fertility, sperm intra-flagellar transport, and tail assembly. PLoS Genet. 8, e1002969. doi:10.1371/journal.pgen.1002969

PubMed Abstract | CrossRef Full Text | Google Scholar

Ma, Q., Cao, C., Zhuang, C., Luo, X., Li, X., Wan, H., et al. (2021). AXDND1, a novel testis-enriched gene, is required for spermiogenesis and male fertility. Cell Death Discov. 7, 348. doi:10.1038/s41420-021-00738-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Machtinger, R., Laurent, L. C., and Baccarelli, A. A. (2016). Extracellular vesicles: roles in gamete maturation, fertilization and embryo implantation. Hum. Reprod. Update 22, 182–193. doi:10.1093/humupd/dmv055

PubMed Abstract | CrossRef Full Text | Google Scholar

Mauduit, C., Hamamah, S., and Benahmed, M. (1999). Stem cell factor/c-kit system in spermatogenesis. Hum. Reprod. Update 5, 535–545. doi:10.1093/humupd/5.5.535

PubMed Abstract | CrossRef Full Text | Google Scholar

Murdica, V., Giacomini, E., Alteri, A., Bartolacci, A., Cermisoni, G. C., Zarovni, N., et al. (2019). Re: seminal plasma of men with severe asthenozoospermia contain exosomes that affect spermatozoa motility and capacitation. J. Urol. 202, 1083. doi:10.1097/01.JU.0000585176.75690.34

PubMed Abstract | CrossRef Full Text | Google Scholar

Narula, A., Kilen, S., Ma, E., Kroeger, J., Goldberg, E., and Woodruff, T. K. (2002). Smad4 overexpression causes germ cell ablation and Leydig cell hyperplasia in transgenic mice. Am. J. Pathol. 161, 1723–1734. doi:10.1016/S0002-9440(10)64449-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Nguyen, T. M., Shafi, A., Nguyen, T., and Draghici, S. (2019). Identifying significantly impacted pathways: a comprehensive review and assessment. Genome Biol. 20, 203–215. doi:10.1186/s13059-019-1790-4

PubMed Abstract | CrossRef Full Text | Google Scholar

O’Brien, J., Hayder, H., Zayed, Y., and Peng, C. (2018). Overview of microRNA biogenesis, mechanisms of actions, and circulation. Front. Endocrinol. (Lausanne) 9, 402. doi:10.3389/fendo.2018.00402

PubMed Abstract | CrossRef Full Text | Google Scholar

Pal, A., Karanwal, S., Chera, J. S., Batra, V., Kumaresan, A., Sarwalia, P., et al. (2023). Circulatory extracellular vesicle derived miR-195-5p promotes cellular apoptosis and suppresses cell proliferation in the buffalo endometrial primary cell culture. Sci. Rep. 13, 16703. doi:10.1038/s41598-023-43530-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Piehl, L. L., Fischman, M. L., Hellman, U., Cisale, H., and Miranda, P. V. (2013). Boar seminal plasma exosomes: effect on sperm function and protein identification by sequencing. Theriogenology 79, 1071–1082. doi:10.1016/j.theriogenology.2013.01.028

PubMed Abstract | CrossRef Full Text | Google Scholar

Pons-Rejraji, H., Artonne, C., Sion, B., Brugnon, F., Canis, M., Janny, L., et al. (2011). Prostasomes: inhibitors of capacitation and modulators of cellular signalling in human sperm. Int. J. Androl. 34, 568–580. doi:10.1111/j.1365-2605.2010.01116.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Prakash, M. A., Kumaresan, A., Ebenezer Samuel King, J. P., Nag, P., Sharma, A., Sinha, M. K., et al. (2021). Comparative transcriptomic analysis of spermatozoa from high- and low-fertile crossbred bulls: implications for fertility prediction. Front. Cell Dev. Biol. 9, 647717. doi:10.3389/fcell.2021.647717

PubMed Abstract | CrossRef Full Text | Google Scholar

Qi, Y., Jiang, M., Yuan, Y., Bi, Y., Zheng, B., Guo, X., et al. (2013). Adp-ribosylation factor-like 3, a manchette-associated protein, is essential for mouse spermiogenesis. Mol. Hum. Reprod. 19, 327–335. doi:10.1093/molehr/gat001

PubMed Abstract | CrossRef Full Text | Google Scholar

Qing, X., Shi, J., Dong, T., Wu, C., Hu, L., and Li, H. (2017). Dysregulation of an X-linked primate-specific epididymal microRNA cluster in unexplained asthenozoospermia. Oncotarget 8, 56839–56849. doi:10.18632/oncotarget.18076

PubMed Abstract | CrossRef Full Text | Google Scholar

Reynolds, N., and Cooke, H. J. (2005). Role of the DAZ genes in male fertility. Reprod. Biomed. Online 10, 72–80. doi:10.1016/S1472-6483(10)60806-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Saez, F., Frenette, G., and Sullivan, R. (2003). Epididymosomes and prostasomes: their roles in post testicular maturation of the sperm cells. J. Androl. 24, 149–154. doi:10.1002/j.1939-4640.2003.tb02653.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Salas-Huetos, A., and Aston, K. I. (2021). Defining new genetic etiologies of male infertility: progress and future prospects. Transl. Androl. Urol. 10, 1486–1498. doi:10.21037/tau.2020.03.43

PubMed Abstract | CrossRef Full Text | Google Scholar

Salas-Huetos, A., Blanco, J., Vidal, F., Godo, A., Grossmann, M., Pons, M. C., et al. (2015). Spermatozoa from patients with seminal alterations exhibit a differential micro-ribonucleic acid profile. Fertil. Steril. 104 (3), 591–601. doi:10.1016/j.fertnstert.2015.06.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Sharer, J. D., Shern, J. F., Van Valkenburgh, H., Wallace, D. C., and Kahn, R. A. (2002). ARL2 and BART enter mitochondria and bind the adenine nucleotide transporter. Mol. Biol. Cell 13, 71–83. doi:10.1091/mbc.01-05-0245

PubMed Abstract | CrossRef Full Text | Google Scholar

Shen, Y., Wu, X., Li, Q., Huang, X., Wang, J., Zhao, L., et al. (2022). Identification and potential value of candidate genes in patients with non-obstructive azoospermia. Urology 164, 133–139. doi:10.1016/j.urology.2022.02.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Somashekar, L., Selvaraju, S., Parthipan, S., Patil, S. K., Binsila, B. K., Venkataswamy, M. M., et al. (2017). Comparative sperm protein profiling in bulls differing in fertility and identification of phosphatidylethanolamine-binding protein 4, a potential fertility marker. Andrology 5 (5), 1032–1051. doi:10.1111/andr.12404

PubMed Abstract | CrossRef Full Text | Google Scholar

Szatanek, R., Baj-Krzyworzeka, M., Zimoch, J., Lekka, M., Siedlar, M., and Baran, J. (2017). The methods of choice for extracellular vesicles (EVs) characterization. Int. J. Mol. Sci. 18, 1153. doi:10.3390/ijms18061153

PubMed Abstract | CrossRef Full Text | Google Scholar

Takahashi, A., Okada, R., Nagao, K., Kawamata, Y., Hanyu, A., Yoshimoto, S., et al. (2017). Exosomes maintain cellular homeostasis by excreting harmful DNA from cells. Nat. Commun. 8, 15287. doi:10.1038/ncomms15287

PubMed Abstract | CrossRef Full Text | Google Scholar

Talluri, T. R., Mal, G., and Ravi, S. K. (2017). Biochemical components of seminal plasma and their correlation to the fresh seminal characteristics in Marwari stallions and Poitou jacks. Vet. World 10, 214–220. doi:10.14202/vetworld.2017.214-220

PubMed Abstract | CrossRef Full Text | Google Scholar

Tanga, B. M., Qamar, A. Y., Raza, S., Bang, S., Fang, X., Yoon, K., et al. (2021). Semen evaluation: methodological advancements in sperm quality-specific fertility assessment - a review. Annu. Rev. Anim. Biosci. 34, 1253–1270. doi:10.5713/ab.21.0072

CrossRef Full Text | Google Scholar

Van Rooij, E. (2011). The art of microRNA research. Circ. Res. 108 (2), 219–234. doi:10.1161/circresaha.110.227496

PubMed Abstract | CrossRef Full Text | Google Scholar

Vojtech, L., Woo, S., Hughes, S., Levy, C., Ballweber, L., Sauteraud, R. P., et al. (2014). Exosomes in human semen carry a distinctive repertoire of small non-coding RNAs with potential regulatory functions. Nucleic Acids Res. 42, 7290–7304. doi:10.1093/nar/gku347

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, R.-A., and Zhao, G.-Q. (1999). Transforming growth factor signal transducer Smad2 is expressed in mouse meiotic germ cells, sertoli cells, and leydig cells during spermatogenesis. Biol. Reprod. 61, 999–1004. doi:10.1095/biolreprod61.4.999

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, S., Qiao, H., Wang, P., Wang, Y., and Qin, D. (2021). ZDHHC19 is dispensable for spermatogenesis, but is essential for sperm functions in mice. Int. J. Mol. Sci. 22, 8894. doi:10.3390/ijms22168894

PubMed Abstract | CrossRef Full Text | Google Scholar

Xing, K., Gao, M., Li, X., Feng, Y., Ge, Y., Qi, X., et al. (2020). An integrated analysis of testis miRNA and mRNA transcriptome reveals important functional miRNA-targets in reproduction traits of roosters. Reprod. Biol. 20 (3), 433–440. doi:10.1016/j.repbio.2020.03.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, Z., Xie, Y., Zhou, C., Hu, Q., Gu, T., Yang, J., et al. (2020). Expression pattern of seminal plasma extracellular vesicle small RNAs in boar semen. Front. Vet. Sci. 7, 585276. doi:10.3389/fvets.2020.585276

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, Q., Qiao, X., Li, D., Chen, B., Zhang, L., Yuan, C., et al. (2019). Expression and association of IL-21, FBXL20 and tumour suppressor gene PTEN in laryngeal cancer. Saudi J. Biol. Sci. 26, 2048–2051. doi:10.1016/j.sjbs.2019.08.013

PubMed Abstract | CrossRef Full Text | Google Scholar