Analysis of miRNAs in milk of four livestock species

BACKGROUND

Milk is essential for mammalian nutrition because it provides vital nutrients for growth and development. Milk composition, which is influenced by genetic and environmental factors, supports lactation, a complex process crucial for milk production and quality. Recent research has focused on noncoding RNAs, particularly microRNAs (miRNAs), which are present in body fluids and regulate gene expression post-transcriptionally. This study comprehensively characterizes miRNAs in milk of four livestock species, namely Bubalus bubalis, Capra hircus, Equus asinus, and Ovis aries and identifies potential target genes.

RESULTS

High-throughput sequencing of milk RNA resulted in distinct read counts across species: B. bubalis (8,790,441 reads), C. hircus (12,976,275 reads), E. asinus (9,385,067 reads), and O. aries (7,295,297 reads). E. asinus had the highest RNA mapping rate (94.6%) and O. aries the lowest (84.8%). A substantially greater proportion of miRNAs over other small RNAs was observed for the donkey milk sample (7.74%) compared to buffalo (0.87%), goat (1.57%), and sheep (1.12%). Shared miRNAs, which included miR-200a, miR-200b, miR-200c, and miR-23a among others, showed varying expression levels across species, confirmed by qPCR analysis. Functional annotation of predicted miRNA target genes highlighted diverse roles, with an enrichment in functions linked to metabolism and immunity. Pathway analysis identified immune response pathways as significant, with several miRNAs targeting specific genes across species, suggesting their regulatory function in milk.

CONCLUSIONS

Both conserved and species-specific miRNAs were detected in milk of the investigated species. The identified target genes of these miRNAs have important roles in neonatal development, adaptation, growth, and immune response. Furthermore, they influence milk and meat production traits in livestock.

Effect of milk stasis on mammary gland involution and the microRNA profile

The presence of an autocrine factor in milk that can trigger mammary gland involution was proposed more than 50 yr ago. To provide evidence for the existence of one or more autocrine factors, 10 multiparous cows in late lactation were quarter-milked for 7 d. Following this baseline period, the right front quarter of each cow was left unmilked, and the other quarters were milked for 7 d. Before the last milking of that period, milk (mammary secretions) was collected aseptically from both front quarters. After that milking, 250 mL of the collected samples were infused in the cows' respective rear quarters. No quarters were milked for the following 7 d (milk stasis period), and quarter milking was then resumed in all quarters for the last 7 d of the experiment (remilking period). Quarter milk samples were collected during the baseline period, before the milk stasis period, and during the remilking period. These samples were used for measuring milk components and the concentration of involution markers (SCC, BSA, and lactoferrin). Samples of mammary secretions were collected manually from the quarters during the milk stasis period for involution marker determination. We extracted RNA from samples collected from front quarters before the last milking before the milk stasis period for microRNA (miRNA) determination. As anticipated, the longer milk stasis period implemented for the right front quarter resulted in a more advanced involution than in the left front quarter, based on the concentration of involution markers in the mammary secretions, lower milk production recovery, and changes in milk composition during the remilking period. All 3 involution marker concentrations in the mammary secretions increased in both rear quarters, but were greater in the right quarter secretions than in the left quarter secretions. Resuming milking reinitiated milk production in all quarters, but milk production recovery in the right rear quarters was less robust than that in the left rear quarters (54.3 ± 1.4% vs. 61.6 ± 1.4%, respectively). Milk from the quarters infused with mammary secretions (right rear) had a lower lactose content, but a higher milk protein content and higher SCC than the quarters infused with milk. We detected a total of 359 miRNAs, 76 of which were differentially expressed in milk and mammary secretions. Expression of bta-miR-221 and bta-miR-223 was upregulated in mammary secretions 34- and 40-fold, respectively. The results of the present experiment support the contention that milk stasis leads to the accumulation of one or more factors that trigger involution. The results also indicate that milk stasis leads to changes in the miRNA profile of the milk, but whether such changes are a cause or a consequence of the involution process remains to be established.

Milk fat globule size: Unraveling the intricate relationship between metabolism, homeostasis, and stress signaling

Fat is an important component of milk which delivers energy, nutrients, and bioactive molecules from the lactating mother to the suckling neonate. Milk fat consists of a complex mixture of different types of lipids; hundreds of fatty acids, triglycerides, phospholipids, sphingolipids, cholesterol and cholesteryl ester, and glycoconjugates, secreted by the mammary gland epithelial cells (MEC) in the form of a lipid-protein assembly termed the milk fat globule (MFG). The mammary gland in general, and specifically that of modern dairy cows, faces metabolic stress once lactation commences, which changes the lipogenic capacity of MECs directly by reducing available energy and reducing factors required for both lipid synthesis and secretion or indirectly by activating a proinflammatory response. Both processes have the capacity to change the morphometric features (e.g., number and size) of the secreted MFG and its precursor-the intracellular lipid droplet (LD). The MFG size is tightly associated with its lipidome and proteome and also affects the bioavailability of milk fat and protein. Thus, MFG size has the potential to regulate the bioactivity of milk and dairy products. MFG size also plays a central role in the functional properties of milk and dairy products such as texture and stability. To understand how stress affects the structure-function of the MFG, we cover: (i) The mechanism of production and secretion of the MFG and the implications of MFG size, (ii) How the response mechanisms to stress can change the morphometric features of MFGs, and (iii) The possible consequences of such modifications.

Organelles coordinate milk production and secretion during lactation: Insights into mammary pathologies

The mammary gland undergoes a spectacular series of changes during its development and maintains a remarkable capacity to remodel and regenerate during progression through the lactation cycle. This flexibility of the mammary gland requires coordination of multiple processes including cell proliferation, differentiation, regeneration, stress response, immune activity, and metabolic changes under the control of diverse cellular and hormonal signaling pathways. The lactating mammary epithelium orchestrates synthesis and apical secretion of macromolecules including milk lipids, milk proteins, and lactose as well as other minor nutrients that constitute milk. Knowledge about the subcellular compartmentalization of these metabolic and signaling events, as they relate to milk production and secretion during lactation, is expanding. Here we review how major organelles (endoplasmic reticulum, Golgi apparatus, mitochondrion, lysosome, and exosome) within mammary epithelial cells collaborate to initiate, mediate, and maintain lactation, and how study of these organelles provides insight into options to maintain mammary/breast health.

The Role of microRNAs in the Mammary Gland Development, Health, and Function of Cattle, Goats, and Sheep

Milk is an integral and therefore complex structural element of mammalian nutrition. Therefore, it is simple to conclude that lactation, the process of producing milk, is as complex as the mammary gland, the organ responsible for this biochemical activity. Nutrition, genetics, epigenetics, disease pathogens, climatic conditions, and other environmental variables all impact breast productivity. In the last decade, the number of studies devoted to epigenetics has increased dramatically. Reports are increasingly describing the direct participation of microRNAs (miRNAs), small noncoding RNAs that regulate gene expression post-transcriptionally, in the regulation of mammary gland development and function. This paper presents a summary of the current state of knowledge about the roles of miRNAs in mammary gland development, health, and functions, particularly during lactation. The significance of miRNAs in signaling pathways, cellular proliferation, and the lipid metabolism in agricultural ruminants, which are crucial in light of their role in the nutrition of humans as consumers of dairy products, is discussed.

MicroRNA-193a-5p Regulates the Synthesis of Polyunsaturated Fatty Acids by Targeting Fatty Acid Desaturase 1 (FADS1) in Bovine Mammary Epithelial Cells

Cardiovascular diseases (CVDs) are seriously threatening to human life and health. Polyunsaturated fatty acids (PUFAs) are known for their role in preventing CVDs. It is beneficial to population health to promote the content of PUFAs in bovine milk. In recent years, limited research based on molecular mechanisms has focused on this field. The biological roles of numerous microRNAs (miRNAs) remain unknown. In this study, a promising and negatively correlated pair of the miRNA (miRNA-193a-5p) and a fatty acid desaturase 1 (FADS1) gene are identified and screened to explore whether they are potential factors of PUFAs’ synthesis in bovine milk. The targeted relationship between miRNA-193a-5p and FADS1 in bovine mammary epithelial cells (BMECs) is demonstrated by dual luciferase reporter assays. qRT-PCR and western blot assays indicate that both the expression of mRNA and the protein FADS1 show a negative correlation with miRNA-193a-5p expression in BMECs. Also, miR-193a-5p expression is positively correlated with the expression of genes associated with milk fatty acid metabolism, including ELOVL fatty acid elongase 6 (ELOVL6) and diacylglycerol O-acyltransferase 2 (DGAT2). The expression of fatty acid desaturase 2 (FADS2) is negatively correlated with miR-193a-5p expression in BMECs. The contents of triglycerides (TAG), eicosapentaenoic acid (EPA), docosapentaenoic acid (DPA) and docosahexaenoic acid (DHA) have a significant positive correlation with the expression of FADS1 and a significant negative correlation with the expression of miR-193a-5p in BMECs. For the first time, this study confirms that miRNA-193a-5p regulates PUFAs metabolism in BMECs by targeting FADS1, indicating that miRNA-193a-5p and FADS1 are underlying factors that improve PUFAs content in bovine milk.