Seasonal effects on miRNA and transcriptomic profile of oocytes and follicular cells in buffalo (Bubalus bubalis)

Season clearly influences oocyte competence in buffalo (Bubalus bubalis); however, changes in the oocyte molecular status in relation to season are poorly understood. This study characterizes the microRNA (miRNA) and transcriptomic profiles of oocytes (OOs) and corresponding follicular cells (FCs) from buffalo ovaries collected in the breeding (BS) and non-breeding (NBS) seasons. In the BS, cleavage and blastocyst rates are significantly higher compared to NBS. Thirteen miRNAs and two mRNAs showed differential expression (DE) in FCs between BS and NBS. DE-miRNAs target gene analysis uncovered pathways associated with transforming growth factor β (TGFβ) and circadian clock photoperiod. Oocytes cluster in function of season for their miRNA content, showing 13 DE-miRNAs between BS and NBS. Between the two seasons, 22 differentially expressed genes were also observed. Gene Ontology (GO) analysis of miRNA target genes and differentially expressed genes (DEGs) in OOs highlights pathways related to triglyceride and sterol biosynthesis and storage. Co-expression analysis of miRNAs and mRNAs revealed a positive correlation between miR-296-3p and genes related to metabolism and hormone regulation. In conclusion, season significantly affects female fertility in buffalo and impacts on oocyte transcriptomic of genes related to folliculogenesis and acquisition of oocyte competence.