Genetic Basis of Follicle Development in Dazu Black Goat by Whole-Transcriptome Sequencing

SIMPLE SUMMARY: The follicle development (FD) of a goat is precisely regulated by various noncoding RNAs (ncRNAs), especially by the regulatory mechanism of competing endogenous RNAs (ceRNAs). This study aimed to determine the expression patterns of messenger RNA (mRNA), long noncoding RNA, microRNA, and circular RNA during the FD of Dazhu black goats by whole-transcriptomic sequencing and analyze the regulatory mechanism of the ncRNA and ceRNA regulatory network. The results may lay a foundation for further research on FD and improving the reproductive performance of goats. ABSTRACT: The follicle development (FD) is an important factor determining litter size in animals. Recent studies have found that noncoding RNAs (ncRNAs) play an important role in FD. In particular, the role of the regulatory mechanism of competing endogenous RNAs (ceRNAs) that drive FD has attracted increasing attention. Therefore, this study explored the genetic basis of goat FD by obtaining the complete follicular transcriptome of Dazu black goats at different developmental stages. Results revealed that 128 messenger RNAs (mRNAs), 4 long noncoding RNAs (lncRNAs), 49 microRNAs (miRNAs), and 290 circular RNAs (circRNAs) were significantly differentially expressed (DE) between large and small follicles. Moreover, DEmRNAs were enriched in many signaling pathways related to FD, as well as GO terms related to molecular binding and enzyme activity. Based on the analysis of the ceRNA network (CRN), 34 nodes (1 DElncRNAs, 10 DEcircRNAs, 14 DEmiRNAs, and 9 DEmRNAs) and 35 interactions (17 DEcircRNAs–DEmRNAs, 2 DElncRNAs–DEmiRNAs, and 16 DEmRNA–DEmiRNAs) implied that the CRN could be involved in the FD of goats. In conclusion, we described gene regulation by DERNAs and lncRNA/circRNA–miRNA–mRNA CRNs in the FD of goats. This study provided insights into the genetic basis of FD in precise transcriptional regulation.