CircRNA screening and ceRNA network construction for milk fat metabolism in dairy cows

Abstract

Background: Milk fat is one of the main reference elements for evaluating milk quality and is a primary objective trait in dairy cattle breeding. In recent years, circular RNAs (circRNAs) have been found to play crucial roles in many biological processes. However, the function and expression profiles of circRNAs in milk fat synthesis in cows are not completely understood. We performed RNA sequencing to analyze the genome-wide expression of circRNA transcripts in bovine mammary epithelial cells (BMECs) from cows with extreme differences in milk fat percentage. We identified candidate differential circRNAs associated with milk fat metabolism using functional enrichment analysis and constructed a lipid metabolism-related competing endogenous RNA (ceRNA) interactive regulatory network.

Results: A total of 290 circRNAs were significantly differentially expressed (DE-circRNAs) in high milk fat percentage (HMF) cows compared to that in low milk fat percentage (LMF) cows. Of the 290 circRNAs, 142 were significantly upregulated and 148 were significantly downregulated. Enrichment analysis (Gene Ontology and Kyoto Encyclopedia of Genes and Genomes) identified four DE-circRNAs (circ_0001122, circ_0007367, circ_0018269, and circ_0015179) that potentially regulate milk fat metabolism. Among them, circ_0001122, circ_0007367, and circ_0015179 had relatively high expression levels in cow mammary gland tissue compared to other tissues (heart, liver, kidney, uterus, ovaries, and small intestine) of cows. The regulatory networks circ_0001122:miR-12043:LIPG, circ_0007367:miR-331-3p:CIDEA/PML, and circ_0018269:miR-11989:RORC/HPX are potential networks to explore the mechanism of milk fat regulation.

Conclusions: These results reveal the possible role of circRNAs in milk fat metabolism in dairy cows. Several important circRNAs and ceRNAs affecting milk fat synthesis were identified, providing insights into the complex biology of milk fat synthesis as well as a novel theoretical perspective for future research on lactation, milk quality, and breed improvement in dairy cows.

Keywords: RNA-seq; ceRNAs; circRNAs; dairy cows; milk fat percentage.