Abstract
Small non-coding RNAs have been linked to different phenotypes in bovine sperm, however attempts to identify sperm-borne molecular biomarkers of male fertility have thus far failed to identify a robust profile of expressed miRNAs related to fertility. We hypothesized that some differences in bull fertility may be reflected in the levels of different miRNAs in sperm. To explore such differences in fertility that are not due to differences in visible metrics of sperm quality, we employed Next Generation Sequencing to compare the miRNA populations in Bos taurus sperm from bulls with comparable motility and morphology but varying Sire Conception Rates. We identified the most abundant miRNAs in both populations (miRs -34b-3p; -100-5p; -191-5p; -30d-4p; -21-5p) and evaluated differences in the overall levels and specific patterns of isomiR expression. We also explored correlations between specific pairs of miRNAs in each population and identified 10 distinct pairs of miRNAs that were positively correlated in bulls with higher fertility and negatively correlated in comparatively less fertile individuals. Furthermore, 8 additional miRNA pairs demonstrated the opposite trend; negatively correlated in high fertility animals and positively correlated in less fertile bulls. Finally, we performed pathway analysis to identify potential roles of miRNAs present in bull sperm in the regulation of specific genes that impact spermatogenesis and embryo development. Together, these results present a comprehensive picture of the bovine sperm miRNAome that suggests multiple potential roles in fertility.