Identification of microRNA transcriptome throughout the lifespan of yak (Bos grunniens) corpus luteum

Abstract

The corpus luteum (CL) is a temporary organ that plays a critical role for female fertility by maintaining the estrous cycle. MicroRNA (miRNA) is a class of non-coding RNAs involved in various biological processes. However, there exists limited knowledge of the role of miRNA in yak CL. In this study, we used high-throughput sequencing to study the transcriptome dynamics of miRNA in yak early (eCL), middle (mCL) and late-stage CL (lCL). A total of 6,730 miRNAs were identified, including 5,766 known and 964 novels miRNAs. Three miRNAs, including bta-miR-126-3p, bta-miR-143 and bta-miR-148a, exhibited the highest expressions in yak CLs of all the three stages. Most of the miRNAs were 20–24 nt in length and the peak was at 22 nt. Besides, most miRNAs with different lengths displayed significant uracil preference at the 5'-end. Furthermore, 1,067, 280 and 112 differentially expressed (DE) miRNAs were found in eCL vs. mCL, mCL vs. lCL, and eCL vs. lCL, respectively. Most of the DE miRNAs were down-regulated in the eCL vs. mCL and eCL vs. lCL groups, and up-regulated in the mCL vs. lCL group. A total of 18,904 target genes were identified, with 18,843 annotated. Pathway enrichment analysis of the DE miRNAs target genes illustrated that the most enriched cellular process in each group included pathways in cancer, PI3K-Akt pathway, endocytosis, and focal adhesion. A total of 20 putative target genes in 47 DE miRNAs were identified to be closely associated with the formation, function or regression of CL. Three DE miRNAs, including bta-miR-11972, novel-miR-619 and novel-miR-153, were proved to directly bind to the 3'-UTR of their predicated target mRNAs, including CDK4, HSD17B1 and MAP1LC3C, respectively. Both of these DE miRNAs and their target mRNAs exhibited dynamic expression profiles across the lifespan of yak CL. This study presents a general basis for understanding of the regulation of miRNA on yak CL and also provides a novel genetic resource for future analysis of the gene network during the estrous cycle in the yak.