miRNA profiling from peripheral white blood cells in fertile and subfertile beef heifers at weaning
The miRNA profiles were measured using small-RNA sequencing in beef heifers sampled at weaning that was retrospectively classified as fertile or subfertile following the breeding protocol. To accomplish this, the miRNA profiles were generated from the blood samples (10 mL) collected from crossbred heifers (Angus-Simmental) at the time of weaning (~238 days after birth). Peripheral white blood cells (PWBC) were extracted from the blood samples and stored at -80°C until further processing. During the breeding season, all the heifers followed the same breeding protocol, estrus synchronization, and fixed-time artificial insemination (AI). Fourteen days following the fixed-time AI, the non-pregnant heifers were exposed to fertile bulls for 60 days. Depending on the presence or absence of conceptus at 75 days following AI, heifers were classified as fertile for those who were pregnant through artificial insemination, pregnant to natural breeding (P-NB), or subfertile for those who were not pregnant. Heifers from fertile (n = 7) and subfertile (n = 7) groups were considered for the study. Total RNA was extracted from the PWBC of 14 samples and was subjected to small RNA library preparation and sequencing. After quality control, adapter trimming, and alignment, mature miRNAs were used for differential expression analysis. The read counts were transformed to counts per million (CPM), and raw counts with CPM < 1 in 50% of the samples were filtered out. The filtered raw counts were analyzed using DESeq2 v 1.26.0 to identify differentially expressed miRNAs (DEMIs). The DEMIs identified with a p-value < 0.05 and absolute (log2 fold change) > 0.5 were considered significant. With the subfertile heifers as the reference group, we identified 16 DEMIs between fertile and subfertile groups. To determine the genes targeting the DEMIs, we downloaded the target genes for each DEMI and retained only those genes expressed in the PWBCs. For the miRNA-gene correlation, we used the partial correlation and information theory (PCIT) approach to identify the significant gene-miRNA correlated pairs. The significant genes correlated with the miRNAs identified pathways including MAPK, ErbB, HIF-1, FoxO, p53, mTOR, T-cell receptor, insulin and GnRH signaling, apoptosis, and pathways regulating pluripotency of stem cells in the fertile group while cell cycle, p53 signaling pathway and apoptosis pathways in the subfertile group. Overall design: Peripheral white blood cell miRNA profiles from 14 Angus-Simmental crossbred heifers (fertile = 7; subfertile = 7) were generated with the Illumina Next-Seq 500 platform using the single-end 50 bp chemistry.