Characterization of microRNA profile in mammary tissue of dairy and beef breed heifers

Milk-derived exosome nanovesicles: recent progress and daunting hurdles

Raw milk is the foundation of quality and safety in the dairy industry, and improving milk source management is the fundamental guarantee. Milk-derived exosomes (MDEs) are nanoscale information transfer molecules secreted by mammary cells with unique content and high stability, which can be used not only as potential markers to analyze key traits of lactation, reproduction, nutrition and health of animals, but also help farm managers to take timely interventions to improve animal welfare, milk quality, and functional traits. Our review first outlines the latest advances in MDEs isolation and purification, compositional analysis and characterization tools. We then provide a comprehensive summary of recent applications of MDEs liquid biopsy in breed selection, disease prevention and control, and feeding management. Finally, we evaluate the impact of processing on the stability of MDEs to offer guidance for dairy production and storage. The limitations and challenges in the development and use of MDEs markers are also discussed. As a noninvasive marker with high sensitivity and specificity, the MDEs-mediated assay technology is expected to be a powerful tool for measuring cow health and raw milk quality, enabling dynamic and precise regulation of dairy cows and full traceability of raw milk.

A Comprehensive Review of Bovine Colostrum Components and Selected Aspects Regarding Their Impact on Neonatal Calf Physiology

Colostrum contains macro- and micronutrients necessary to meet the nutritional and energy requirements of the neonatal calf, bioactive components that intervene in several physiological aspects, and cells and microorganisms that modulate the calf's immune system and gut microbiome. Colostrum is sometimes mistaken as transition milk, which, although more nutritive than whole milk, has a distinct biochemical composition. Furthermore, most research about colostrum quality and colostrum management focuses on the transfer of maternal IgG to the newborn calf. The remaining components of colostrum and transition milk have not received the same attention, despite their importance to the newborn animal. In this narrative review, a large body of literature on the components of bovine colostrum was reviewed. The variability of these components was summarized, emphasizing specific components that warrant deeper exploration. In addition, the effects of each component present in colostrum and transition milk on several key physiological aspects of the newborn calf are discussed.

Whole-transcriptome analysis of longissimus dorsi muscle in cattle-yaks reveals the regulatory functions of ADAMTS6 gene in myoblasts

Cattle-yak, which is the hybrid F1 generation of cattle and yak, demonstrates better production performance compared to yak. However, there is limited research on the molecular mechanisms responsible for the muscle development of cattle-yak. To address this knowledge gap, a comprehensive transcriptomic survey of the longissimus dorsi muscle in cattle-yak was conducted. Three transcript types, namely lncRNAs, miRNAs, and circRNAs, along with protein-coding genes were characterized at two developmental stages (6 m, 18 m) of cattle-yak. The results revealed significant enrichment of these transcripts into pathways related to myoblast differentiation and muscle development signaling. Additionally, the study identified the TCONS00024465/circHIPK3-bta-miR-499-ADAMTS6 regulatory network, which may play a crucial role in the muscle development of cattle-yak by combining the transcriptome data of yak and constructing the ceRNA co-expression network. HEK 293 T cells were used to validate that TCONS00024465 and circHIPK3 are located upstream of bta-miR-499, and can competitively bind to bta-miR-499 as ceRNA. The study also verified that ADAMTS6 regulates skeletal muscle development by inhibiting myoblast proliferation, promoting myoblast differentiation, and positively regulating the apoptosis of myoblasts. Taken together, this study provides new insights into the advantages of cattle-yak production performance and offers a molecular basis for further research on muscle development.

How Does Nutrition Affect the Epigenetic Changes in Dairy Cows?

Dairy cows require a balanced diet that provides enough nutrients to support milk production, growth, and reproduction. Inadequate nutrition can lead to metabolic disorders, impaired fertility, and reduced milk yield. Recent studies have shown that nutrition can affect epigenetic modifications in dairy cows, which can impact gene expression and affect the cows' health and productivity. One of the most important epigenetic modifications in dairy cows is DNA methylation, which involves the addition of a methyl group to the DNA molecule. Studies have shown that the methylation status of certain genes in dairy cows can be influenced by dietary factors such as the level of methionine, lysine, choline, and folate in the diet. Other important epigenetic modifications in dairy cows are histone modification and microRNAs as regulators of gene expression. Overall, these findings suggest that nutrition can have a significant impact on the epigenetic regulation of gene expression in dairy cows. By optimizing the diet of dairy cows, it may be possible to improve their health and productivity by promoting beneficial epigenetic modifications. This paper reviews the main nutrients that can cause epigenetic changes in dairy cattle by analyzing the effect of diet on milk production and its composition.

Identification of Key Genes and Potential Pathways Associated with Mastitis Induced by E. coli

Mastitis is one of the most important diseases of dairy cattle. It is an infectious disease leading to an inflammatory reaction in the cow's mammary gland. Escherichia coli is one of the common bacteria which induce mastitis in cows. The aim of this study was to identify key genes and potential pathways associated with mastitis induced by E. coli in dairy cattle using bioinformatics analysis. The gene expression profile of ten samples including five adjacent tissues from a quarter infected with Escherichia coli and five tissues from a healthy quarter of dairy cattle was assessed using GEO2R. Gene ontology and pathway analysis were performed using bioinformatics tools. A total of 156 differentially expressed genes were detected which 95 genes were upregulated and 61 genes were downregulated in adjacent tissue of quarter infected compared with healthy tissue. Cellular oxidant detoxification and oxidation-reduction process were the most significant biological process terms in gene ontology analysis. The most important pathways of DEGs were the biosynthesis of amino acids, p53 signaling pathway, and Metabolic pathways. Three important modules were identified and their path enrichment analysis was performed. There are 10 core genes, among which SOD2, COL1A2, COL3A1, POSTN, ALDH18A1, and CBS may be the main genes associated with mastitis, which can be considered as candidate genes in the prevention and carly diagnosis program of mastitis.

Bta-miR-106b Regulates Bovine Mammary Epithelial Cell Proliferation, Cell Cycle, and Milk Protein Synthesis by Targeting the CDKN1A Gene

Our previous studies found that bta-miR-106b and its corresponding target gene, CDKN1A, were differentially expressed between the mammary epithelium of lactating Holstein cows with extremely high and low milk protein and fat percentage, implying the potential role of bta-miR-106b in milk composition synthesis. In this study, with luciferase assay experiment, bta-miR-106b was validated to target the 3'-untranslated region (UTR) of bovine CDKN1A, thereby regulating its expression. Moreover, in bovine mammary epithelial cells (BMECs), over-expression of bta-miR-106b significantly down-regulated the CDKN1A expression at both mRNA and protein levels, and inhibitors of bta-miR-106b increased CDKN1A expression. Of note, we observed that bta-miR-106b accelerated cell proliferation and cell cycle, and changed the expressions of protein synthesis related pathways such as JAK-STAT and PI3K/AKT/mTOR through regulating CDKN1A expression. Our findings highlight the important regulatory role of bta-miR-106b in milk protein synthesis by targeting CDKN1A in dairy cattle.

Comparative Study of the Expression Profiles of miRNAs of Milk-Derived Exosomes of Yak and Jeryak

The Jeryak is the hybrid offspring of yaks and Jersey cattle and exhibit improved milk and meat yields. Biomolecules carried within milk exosomes are important for cell growth, development, immune regulation, and various pathophysiological processes. Previous studies showed that miRNAs regulate mammary gland development, lactation, and milk quality. This study explored the relationship between milk exosomes miRNAs and lactation performance. A comparison of the milk content showed that yak milk was of a better quality compared to Jeryak milk (casein, fat, TS, SNF, lactose). Milk collected in December was superior to that collected in June for both yak and Jeryak, except for lactose concentrations. Exosomes were extracted by density gradient centrifugation and miRNA expression profiles in milk exosomes from three yaks and three Jeryaks collected in June and December were detected by small RNA sequencing. In all, 22, 120, 78, and 62 differentially expressed miRNAs (DEMs) were identified in Jun_ JY vs. Jun_ Y (P1: Jeryak in June vs. Yak in June), Jun_ JY vs. Dec_ JY (P2: Jeryak in June vs. Jeryak in December), Dec_ JY vs. Dec_ Y (P3: Jeryak in December vs. Yak in December), and Jun_ Y vs. Dec_ Y (P4: Yak in June vs. Yak in December) groups. These DEMs were enriched in functions and signaling pathways related to lactation performance. In conclusion, these findings are a reference tool to study the molecular basis of lactation performance.

Determining extracellular vesicles properties and miRNA cargo variability in bovine milk from healthy cows and cows undergoing subclinical mastitis

Background

Subclinical mastitis, the inflammation of the mammary gland lacking clinical symptoms, is one of the most prevalent and costly diseases in dairy farming worldwide. Milk microRNAs (miRNAs) encapsulated in extracellular vesicles (EVs) have been proposed as potential biomarkers of different mammary gland conditions, including subclinical mastitis. However, little is known about the robustness of EVs analysis regarding sampling time-point and natural infections. To estimate the reliability of EVs measurements in raw bovine milk, we first evaluated changes in EVs size and concentration using Tunable Resistive Pulse Sensing (TRPS) during three consecutive days of sampling. Then, we analysed daily differences in miRNA cargo using small RNA-seq. Finally, we compared milk EVs differences from naturally infected udder quarters with their healthy adjacent quarters and quarters from uninfected udders, respectively.

Results

We found that the milk EV miRNA cargo was very stable over the course of three days regardless of the health status of the quarter, and that infected quarters did not induce relevant changes in milk EVs of adjacent healthy quarters. Chronic subclinical mastitis induced changes in milk EV miRNA cargo, but neither in EVs size nor concentration. We observed that the changes in immunoregulatory miRNAs in quarters with chronic subclinical mastitis were cow-individual, however, the most upregulated miRNA was bta-miR-223-3p across all individuals.

Conclusions

Our results showed that the miRNA profile and particle size characteristics remained constant throughout consecutive days, suggesting that miRNAs packed in EVs are physiological state-specific. In addition, infected quarters were solely affected while adjacent healthy quarters remained unaffected. Finally, the cow-individual miRNA changes pointed towards infection-specific alterations.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08377-z.

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.

Integrated Small RNA Sequencing, Transcriptome and GWAS Data Reveal microRNA Regulation in Response to Milk Protein Traits in Chinese Holstein Cattle

Milk protein is one of the most important economic traits in the dairy industry. Yet, the regulatory network of miRNAs for the synthesis of milk protein in mammary is poorly understood. Samples from 12 Chinese Holstein cows with three high ( ≥ 3.5%) and three low ( ≤ 3.0%) phenotypic values for milk protein percentage in lactation and non-lactation were examined through deep small RNA sequencing. We characterized 388 known and 212 novel miRNAs in the mammary gland. Differentially expressed analysis detected 28 miRNAs in lactation and 52 miRNAs in the non-lactating period with a highly significant correlation with milk protein concentration. Target prediction and correlation analysis identified some key miRNAs and their targets potentially involved in the synthesis of milk protein. We analyzed for enrichments of GWAS signals in miRNAs and their correlated targets. Our results demonstrated that genomic regions harboring DE miRNA genes in lactation were significantly enriched with GWAS signals for milk protein percentage traits and that enrichments within DE miRNA targets were significantly higher than in random gene sets for the majority of milk production traits. This integrated study on the transcriptome and posttranscriptional regulatory profiles between significantly differential phenotypes of milk protein concentration provides new insights into the mechanism of milk protein synthesis, which should reveal the regulatory mechanisms of milk secretion.

Characterization of microRNA expression in B cells derived from Japanese black cattle naturally infected with bovine leukemia virus by deep sequencing

Bovine leukemia virus (BLV) is the causative agent of enzootic bovine leukosis (EBL), a malignant B cell lymphoma. However, the mechanisms of BLV-associated lymphomagenesis remain poorly understood. Here, after deep sequencing, we performed comparative analyses of B cell microRNAs (miRNAs) in cattle infected with BLV and those without BLV. In BLV-infected cattle, BLV-derived miRNAs (blv-miRNAs) accounted for 38% of all miRNAs in B cells. Four of these blv-miRNAs (blv-miR-B1-5p, blv-miR-B2-5p, blv-miR-B4-3p, and blv-miR-B5-5p) had highly significant positive correlations with BLV proviral load (PVL). The read counts of 90 host-derived miRNAs (bta-miRNAs) were significantly down-regulated in BLV-infected cattle compared to those in uninfected cattle. Only bta-miR-375 had a positive correlation with PVL in BLV-infected cattle and was highly expressed in the B cell lymphoma tissue of EBL cattle. There were a few bta-miRNAs that correlated with BLV tax/rex gene expression; however, BLV AS1 expression had a significant negative correlation with many of the down-regulated bta-miRNAs that are important for tumor development and/or tumor suppression. These results suggest that BLV promotes lymphomagenesis via AS1 and blv-miRNAs, rather than tax/rex, by down-regulating the expression of bta-miRNAs that have a tumor-suppressing function, and this downregulation is linked to increased PVL.

Nutrigenomic analyses reveal miRNAs and mRNAs affected by feed restriction in the mammary gland of midlactation dairy cows

The objective of this study was to investigate the effects of feed restriction on mammary miRNAs and coding gene expression in midlactation cows. Five Holstein cows and 6 Montbéliarde cows underwent 6 days of feed restriction, during which feed allowance was reduced to meet 50% of their net energy for lactation requirements. Mammary biopsies were performed before and at the end of the restriction period. Mammary miRNA and mRNA analyses were performed using high-throughput sequencing and microarray analyses, respectively. Feed restriction induced a negative energy balance and decreased milk production and fat and protein yields in both breeds. Feed restriction modified the expression of 27 miRNAs and 374 mRNAs in mammary glands from Holstein cows, whereas no significant miRNA change was observed in Montbéliarde cows. Among the 27 differentially expressed miRNAs, 8 miRNAs were associated with dairy QTL. Analysis of target genes indicate that the 8 most abundantly expressed miRNAs control transcripts related to lipid metabolism, mammary remodeling and stress response. A comparison between the mRNAs targeted by the 8 most strongly expressed miRNAs and 374 differentially expressed mRNAs identified 59 mRNAs in common. The bioinformatic analyses of these 59 mRNAs revealed their implication in lipid metabolism and endothelial cell proliferation. These effects of feed restriction on mammary miRNAs and mRNAs observed in Holstein cows suggest a potential role of miRNAs in mammary structure and lipid biosynthesis that could explain changes in milk production and composition.

Epigenetics: New Insights into Mammary Gland Biology

The mammary gland undergoes important anatomical and physiological changes from embryogenesis through puberty, pregnancy, lactation and involution. These steps are under the control of a complex network of molecular factors, in which epigenetic mechanisms play a role that is increasingly well described. Recently, studies investigating epigenetic modifications and their impacts on gene expression in the mammary gland have been performed at different physiological stages and in different mammary cell types. This has led to the establishment of a role for epigenetic marks in milk component biosynthesis. This review aims to summarize the available knowledge regarding the involvement of the four main molecular mechanisms in epigenetics: DNA methylation, histone modifications, polycomb protein activity and non-coding RNA functions.

Characterization of Breed Specific Differences in Spermatozoal Transcriptomes of Sheep in Australia

Reduced reproductive efficiency results in economic losses to the Australian sheep industry. Reproductive success, particularly after artificial insemination, is dependent on a number of contributing factors on both ewe and ram sides. Despite considerable emphasis placed on characterising ewe side contributions, little emphasis has been placed on characterising ram side contributions to conception success. Over 14,000 transcripts are in spermatozoa of other species, which are transferred to the ova on fertilisation. These transcripts conceivably influence early embryonic development and whether conception is successful. Semen was collected (n = 45) across three breeds; Merino, Dohne, and Poll Dorset. Following collection, each ejaculate was split in two; an aliquot was assessed utilising Computer Assisted Semen Analysis (CASA) and the remaining was utilised for RNA extraction and subsequent next-generation sequencing. Overall, 754 differentially expressed genes were identified in breed contrasts and contrast between ejaculates of different quality. Downstream analysis indicated that these genes could play significant roles in a broad range of physiological functions, including maintenance of spermatogenesis, fertilisation, conception, embryonic development, and offspring production performance. Overall results provide evidence that the spermatozoal transcriptome could be a crucial contributing factor in improving reproductive performance as well as in the overall productivity and profitability of sheep industries.

Deep small RNA-Seq reveals microRNAs expression profiles in lactating mammary gland of 2 sheep breeds with different milk performance

Micro ribonucleic acid (miRNA) is a type of noncoding RNA, and it has been revealed to play important roles in the activity of the mammary gland (MG) in some species. However, the function of miRNAs in MG of sheep is poorly understood. In the study, Gansu Alpine Merino (GAM; n = 9) and Small-tailed Han sheep (STH; n = 9) with different milk production traits were investigated. Microstructures and the expression profile of miRNAs of MG tissues at peak lactation were studied. Mature alveolar lumens of MG in appearance were larger in STH than GAM. The expression levels of CSN2 and the content of rough endoplasmic reticulum were also higher in STH ewes than GAM ewes. A total of 124 mature miRNAs were expressed, and 18 of these were differentially expressed between the 2 breeds. The KEGG analysis results showed that the targeted genes of differentially expressed miRNAs were mainly involved in some metabolic pathways and signaling pathways related to MG development, milk protein, and fat synthesis. The findings in the study can improve our understanding of the roles of miRNAs in the development and lactation of MG in sheep.

Small RNA deep sequencing reveals the expressions of microRNAs in ovine mammary gland development at peak-lactation and during the non-lactating period

MicroRNAs (miRNAs) are small non-coding RNAs that are involved in mammary gland development and lactation in livestock. Little is known about the roles of miRNAs in ovine mammary gland development, hence in this study the expression profiles of miRNAs of the mammary gland tissues of ewes at peak-lactation and during the non-lactating period were investigated using RNA sequencing. A total of 147 mature miRNAs were expressed in the two periods. Compared with peak-lactation, eight miRNAs in the non-lactating ewe mammary gland were significantly up-regulated, whereas fifteen miRNAs were down-regulated. A KEGG analysis revealed that the target genes of the up-regulated miRNAs were significantly enriched in lysosome, Wnt and MAPK signaling pathways, while the target genes of down-regulated miRNAs were significantly enriched in the PI3K-Akt signaling pathway, protein processing in endoplasmic reticulum and axon guidance. These results suggest that further study of the differentially expressed miRNAs could provide a better understanding of the molecular mechanisms of mammary development and lactation in sheep.

RumimiR: a detailed microRNA database focused on ruminant species

The ever-increasing use of next-generation sequencing technologies to explore the genome has generated large quantities of data in recent years. Numerous publications have described several thousand sequences of microRNAs, all species included. A new database (RumimiR) has been created from the literature to provide a detailed description of microRNAs for three ruminant species: cattle, goats and sheep. To date, 2887, 2733 and 5095 unique microRNAs from bovine, caprine and ovine species, respectively, are included. In addition to the most recent reference genomic position and sequence of each microRNA, this database contains details about the animals, tissue origins and experimental conditions mentioned in the publications. Identity to human or mouse microRNA is also indicated. The RumimiR database allows data filtering by selecting microRNAs on the basis of defined criteria such as animal status or tissue origin. For ruminant studies, RumimiR supplements the widely used miRBase database, by using complementary criteria to allow browsing and filtering, and integrates all newly described published sequences. The principal goal of this database is to provide easy access to all the ruminant microRNAs described in the literature.

Characterization of Holstein and Normande whole milk miRNomes highlights breed specificities

The concept of milk as a healthy food has opened the way for studies on milk components, from nutrients to microRNAs, molecules with broad regulatory properties present in large quantities in milk. Characterization of these components has been performed in several species, such as humans and bovine, depending on the stages of lactation. Here, we have studied the variation in milk microRNA composition according to genetic background. Using high throughput sequencing, we have characterized and compared the milk miRNomes of Holstein and Normande cattle, dairy breeds with distinct milk production features, in order to highlight microRNAs that are essential for regulation of the lactation process. In Holstein and Normande milk, 2,038 and 2,030 microRNAs were identified, respectively, with 1,771 common microRNAs, of which 1,049 were annotated and 722 were predicted. The comparison of the milk miRNomes of two breeds allowed to highlight 182 microRNAs displaying significant differences in the abundance. They are involved in the regulation of lipid metabolism and mammary morphogenesis and development, which affects lactation. Our results provide new insights into the regulation of molecular mechanisms involved in milk production.

Complexity of the microRNA transcriptome of cow milk and milk-derived extracellular vesicles isolated via differential ultracentrifugation

MicroRNAs (miRNAs) are small gene-regulatory noncoding RNA that are highly enriched in cow milk. They are encapsulated in different extracellular vesicle (EV) subsets that protect them from the extracellular milieu and the harsh conditions of the gastrointestinal tract during digestion. Here, we isolated pellets enriched in 4 different EV subsets, via differential ultracentrifugation of commercial cow milk: 12,000 × g (P12K), 35,000 × g (P35K), 70,000 × g (P70K), and 100,000 × g (P100K). Small RNA sequencing (sRNA-Seq) analyses revealed an unprecedented level of diversity in the complete miRNA repertoire and features of unfractionated cow milk and derived EV subsets. Although 5 miRNA sequences represented more than 50% of all miRNAs, milk EV exhibited heterogeneous content of miRNAs and isomeric variants (termed isomiR): P100K EV were enriched in reference miRNA sequences, and P12K and P35K EV in related isomiR. Incubation of milk EV with human cultured HeLa cells led to cellular enrichment in miRNA miR-223, which was concomitant with decreased expression of a reporter gene placed under the control of miR-223, thereby demonstrating the functionality of miR-223. These results suggest that cow milk EV may transfer their miRNAs to human cells and regulate recipient cell gene expression programming in a manner as complex as that of their miRNA transcriptome. The biological activity and relevance of the different milk EV subsets and bioactive mediators, including small noncoding RNA, in health and disease, warrants further investigation.

Invited review: MicroRNAs in bovine colostrum-Focus on their origin and potential health benefits for the calf

Colostrum is the first milk produced by a cow after she gives birth. Compared with mature milk, it has a high concentration of immunoglobulin G. Calves are born without circulating antibodies, thus ingestion of colostrum is necessary to protect the calf against pathogens in the first challenging weeks of life. In addition to the life-saving supply of antibodies, colostrum contains minerals, vitamins, growth factors, and immune cells. Recently, microRNAs (miRNAs) were added to that list. MicroRNAs are short, non-coding RNA molecules that can regulate gene expression at the post-transcriptional level. They are thought to act as key regulators of diverse biological and developmental processes. Colostrum contains higher amounts of miRNAs than mature milk; immune- and development-related miRNAs are prominent. Their expression pattern in milk is likely to be influenced by maternal nutrition and environment. The fat content of the maternal diet appears to have a major effect on expression of miRNAs in milk and in the neonate. The immunological state of the mammary gland seems to affect miRNA expression as well. In cows diagnosed with subclinical mastitis, alterations in the expression of miRNAs in milk have been observed. It is believed that miRNAs in colostrum and milk are signaling molecules passed from mother to newborn. They are packaged in extracellular vesicles, which makes them resistant to the harsh conditions in the gastrointestinal tract. Therefore, they can reach the small intestine, where they are absorbed and transferred into the bloodstream. MicroRNAs are important for the development of the intestines. For example, miRNAs stimulate cell viability, proliferation, and stem cell activity of the intestinal epithelium. Furthermore, miRNAs seem to act as key players in the development of the complete immune system. They can, among other things, regulate B- and T-cell differentiation and affect interleukin production of macrophages. The abundance of miRNAs in colostrum and milk and the possibility for their absorption in the intestines of the neonate supports the hypothesis that these tiny molecules are important for the development of the newborn. The probable relation of diet to the expression of miRNAs by the mother creates a possible avenue to optimize expression of miRNAs and improve neonatal maturation.

Deep RNA-Seq reveals miRNome differences in mammary tissue of lactating Holstein and Montbéliarde cows

Background

Genetic polymorphisms are known to influence milk production and composition. However, the genomic mechanisms involved in the genetic regulation of milk component synthesis are not completely understood. MicroRNAs (miRNAs) regulate gene expression. Previous research suggests that the high developmental potential of the mammary gland may depend in part on a specific miRNA expression pattern. The objective of the present study was to compare the mammary gland miRNomes of two dairy cow breeds, Holstein and Montbéliarde, which have different mammogenic potentials that are related to differences in dairy performance.

Results

Milk, fat, protein, and lactose yields were lower in Montbéliarde cows than in Holstein cows. We detected 754 distinct miRNAs in the mammary glands of Holstein (n = 5) and Montbéliarde (n = 6) midlactating cows using RNA-Seq technology, among which 738 were known and 16 were predicted miRNAs. The 25 most abundant miRNAs accounted for 90.6% of the total reads. The comparison of their abundances in the mammary glands of Holstein versus Montbéliarde cows identified 22 differentially expressed miRNAs (P_adj ≤ 0.05). Among them, 11 presented a fold change ≥2, and 2 (miR-100 and miR-146b) were highly expressed. Among the most abundant miRNAs, miR-186 is known to inhibit cell proliferation and epithelial-to-mesenchymal transition.

Data mining showed that 17 differentially expressed miRNAs with more than 20 reads were involved in the regulation of mammary gland plasticity. Several of them may potentially target mRNAs involved in signaling pathways (such as mTOR) and lipid metabolism, thereby indicating that they could influence milk composition.

Conclusion

We found differences in the mammary gland miRNomes of two dairy cattle breeds. These differences suggest a potential role for miRNAs in mammary gland plasticity and milk component synthesis, both of which are related to milk production and composition. Further research is warranted on the genetic regulation of miRNAs and their role in milk synthesis.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5987-4) contains supplementary material, which is available to authorized users.

Co-Expression Networks Reveal Potential Regulatory Roles of miRNAs in Fatty Acid Composition of Nelore Cattle

Fatty acid (FA) content affects the sensorial and nutritional value of meat and plays a significant role in biological processes such as adipogenesis and immune response. It is well known that, in beef, the main FAs associated with these biological processes are oleic acid (C18:1 cis9, OA) and conjugated linoleic acid (CLA-c9t11), which may have beneficial effects on metabolic diseases such as type 2 diabetes and obesity. Here, we performed differential expression and co-expression analyses, weighted gene co-expression network analysis (WGCNA) and partial correlation with information theory (PCIT), to uncover the complex interactions between miRNAs and mRNAs expressed in skeletal muscle associated with FA content. miRNA and mRNA expression data were obtained from skeletal muscle of Nelore cattle that had extreme genomic breeding values for OA and CLA. Insulin and MAPK signaling pathways were identified by WGCNA as central pathways associated with both of these fatty acids. Co-expression network analysis identified bta-miR-33a/b, bta-miR-100, bta-miR-204, bta-miR-365-5p, bta-miR-660, bta-miR-411a, bta-miR-136, bta-miR-30-5p, bta-miR-146b, bta-let-7a-5p, bta-let-7f, bta-let-7, bta-miR 339, bta-miR-10b, bta-miR 486, and the genes ACTA1 and ALDOA as potential regulators of fatty acid synthesis. This study provides evidence and insights into the molecular mechanisms and potential target genes involved in fatty acid content differences in Nelore beef cattle, revealing new candidate pathways of phenotype modulation that could positively benefit beef production and human consumption.

MiR-24-3p regulates cell proliferation and milk protein synthesis of mammary epithelial cells through menin in dairy cows

MiR-24-3p, a broadly conserved, small, noncoding RNA, is abundantly expressed in mammary tissue. However, its regulatory role in this tissue remains poorly understood. It was predicted that miR-24-3p targets the 3' untranslated region (3'-UTR) of multiple endocrine neoplasia type 1 (MEN1), an important regulatory factor in mammary tissue. The objective of this study was to investigate the function of miR-24-3p in mammary cells. Using a luciferase assay in mammary epithelial cells (MAC-T), miR-24-3p was confirmed to target the 3'-UTR of MEN1. Furthermore, miR-24-3p negatively regulated the expression of the MEN1 gene and its encoded protein, menin. miR-24-3p enhanced proliferation of MAC-T by promoting G1/S phase progression. MiR-24-3p also regulated the expression of key factors involved in phosphatidylinositol-3-kinase/protein kinase B/mammalian target of rapamycin and Janus kinase/signal transducer and activators of transcription signaling pathways, therefore controlling milk protein synthesis in epithelial cells. Thus, miR-24-3p appears to act on MAC-T by targeting MEN1. The expression of miR-24-3p was controlled by MEN1/menin, indicating a negative feedback loop between miR-24-3p and MEN1/menin. The negatively inhibited expression pattern of miR-24-3p and MEN1 was active in mammary tissues at different lactation stages. The feedback mechanism is a new concept to further understand the lactation cycle of mammary glands and can possibly to be manipulated to improve milk yield and quality.

Characterization of miRNA profiles in the mammary tissue of dairy cattle in response to heat stress

Background

MicroRNAs (miRNAs) are a class of small noncoding RNAs that play important roles in the regulation of gene expression. However, the role of miRNAs in bovine mammary gland responses to heat stress is not well understood.

Results

In the present study, we performed a deep RNA sequencing analysis to identify miRNAs associated with the heat stress potential of the bovine mammary gland. We identified 27 miRNAs that were differentially expressed significantly between the mammary tissue of Holstein cattle heat stress and normal conditions. Twenty miRNAs had higher expression in the mammary tissue of heat-stressed Holstein cattle. The seven highest differentially expressed candidate miRNAs (bta-miR-21-5p, bta-miR-99a-5p, bta-miR-146b, bta-miR-145, bta-miR-2285 t, bta-miR-133a, and bta-miR-29c) identified by deep RNA sequencing were additionally evaluated by stem-loop qPCR. Enrichment analyses for targeted genes revealed that the major differences between miRNAs expression in the mammary gland of heat-stressed versus control were associated with the regulation of Wnt, TGF-β, MAPK, Notch, and JAK-STAT.

Conclusions

These data indicated that the differentially expressed miRNAs identified in this study may act as dominant regulators during heat stress. We might reduce heat stress damage of Holstein cows by up-regulating or down-regulating these differentially expressed miRNAs.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5298-1) contains supplementary material, which is available to authorized users.

Transcriptome Analysis of Long Non-Coding RNA in the Bovine Mammary Gland Following Dietary Supplementation with Linseed Oil and Safflower Oil

This study aimed to characterize the long non-coding RNA (lncRNA) expression in the bovine mammary gland and to infer their functions in dietary response to 5% linseed oil (LSO) or 5% safflower oil (SFO). Twelve cows (six per treatment) in mid lactation were fed a control diet for 28 days followed by a treatment period (control diet supplemented with 5% LSO or 5% SFO) of 28 days. Mammary gland biopsies were collected from each animal on day-14 (D-14, control period), D+7 (early treatment period) and D+28 (late treatment period) and were subjected to RNA-Sequencing and subsequent bioinformatics analyses. Functional enrichment of lncRNA was performed via potential cis regulated target genes located within 50 kb flanking regions of lncRNAs and having expression correlation of >0.7 with mRNAs. A total of 4955 lncRNAs (325 known and 4630 novel) were identified which potentially cis targeted 59 and 494 genes in LSO and SFO treatments, respectively. Enrichments of cis target genes of lncRNAs indicated potential roles of lncRNAs in immune function, nucleic acid metabolism and cell membrane organization processes as well as involvement in Notch, cAMP and TGF-β signaling pathways. Thirty-two and 21 lncRNAs were differentially expressed (DE) in LSO and SFO treatments, respectively. Six genes (KCNF1, STARD13, BCL6, NXPE2, HHIPL2 and MMD) were identified as potential cis target genes of six DE lncRNAs. In conclusion, this study has identified lncRNAs with potential roles in mammary gland functions and potential candidate genes and pathways via which lncRNAs might function in response to LSO and SFA.

Co-Expression Network Analysis Identifies miRNA⁻mRNA Networks Potentially Regulating Milk Traits and Blood Metabolites

MicroRNAs (miRNA) regulate mRNA networks to coordinate cellular functions. In this study, we constructed gene co-expression networks to detect miRNA modules (clusters of miRNAs with similar expression patterns) and miRNA–mRNA pairs associated with blood (triacylglyceride and nonesterified fatty acids) and milk (milk yield, fat, protein, and lactose) components and milk fatty acid traits following dietary supplementation of cows’ diets with 5% linseed oil (LSO) (n = 6 cows) or 5% safflower oil (SFO) (n = 6 cows) for 28 days. Using miRNA transcriptome data from mammary tissues of cows for co-expression network analysis, we identified three consensus modules: blue, brown, and turquoise, composed of 70, 34, and 86 miRNA members, respectively. The hub miRNAs (miRNAs with the most connections with other miRNAs) were miR-30d, miR-484 and miR-16b for blue, brown, and turquoise modules, respectively. Cell cycle arrest, and p53 signaling and transforming growth factor–beta (TGF-β) signaling pathways were the common gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways enriched for target genes of the three modules. Protein percent (p = 0.03) correlated with the turquoise module in LSO treatment while protein yield (p = 0.003) and milk yield (p = 7 × 10⁻⁰⁴) correlated with the turquoise model, protein and milk yields and lactose percent (p < 0.05) correlated with the blue module and fat percent (p = 0.04) correlated with the brown module in SFO treatment. Several fatty acids correlated (p < 0.05) with the blue (CLA:9,11) and brown (C4:0, C12:0, C22:0, C18:1n9c and CLA:10,12) modules in LSO treatment and with the turquoise (C14:0, C18:3n3 and CLA:9,11), blue (C14:0 and C23:0) and brown (C6:0, C16:0, C22:0, C22:6n3 and CLA:10,12) modules in SFO treatment. Correlation of miRNA and mRNA data from the same animals identified the following miRNA–mRNA pairs: miR-183/RHBDD2 (p = 0.003), miR-484/EIF1AD (p = 0.011) and miR-130a/SBSPON (p = 0.004) with lowest p-values for the blue, brown, and turquoise modules, respectively. Milk yield, protein yield, and protein percentage correlated (p < 0.05) with 28, 31 and 5 miRNA–mRNA pairs, respectively. Our results suggest that, the blue, brown, and turquoise modules miRNAs, hub miRNAs, miRNA–mRNA networks, cell cycle arrest GO term, p53 signaling and TGF-β signaling pathways have considerable influence on milk and blood phenotypes following dietary supplementation of dairy cows’ diets with 5% LSO or 5% SFO.

Screening of miRNA profiles and construction of regulation networks in early and late lactation of dairy goat mammary glands

In recent years, studies related to the expression profiles of miRNAs in the dairy goat mammary gland were performed, but regulatory mechanisms in the physiological environment and the dynamic homeostasis of mammary gland development and lactation are not clear. In the present study, sequencing data analysis of early and late lactation uncovered a total of 1,487 unique miRNAs, including 45 novel miRNA candidates and 1,442 known and conserved miRNAs, of which 758 miRNAs were co-expressed and 378 differentially expressed with P < 0.05. Moreover, 76 non-redundant target genes were annotated in 347 GO consortiums, with 3,143 candidate target genes grouped into 33 pathways. Additionally, 18 predicted target genes of 214 miRNAs were directly annotated in mammary gland development and used to construct regulatory networks based on GO annotation and the KEGG pathway. The expression levels of seven known miRNAs and three novel miRNAs were examined using quantitative real-time PCR. The results showed that miRNAs might play important roles in early and late lactation during dairy goat mammary gland development, which will be helpful to obtain a better understanding of the genetic control of mammary gland lactation and development.

MicroRNA roles in signalling during lactation: an insight from differential expression, time course and pathway analyses of deep sequence data

The study examined microRNA (miRNA) expression and regulatory patterns during an entire bovine lactation cycle. Total RNA from milk fat samples collected at the lactogenesis (LAC, day1 [D1] and D7), galactopoiesis (GAL, D30, D70, D130, D170 and D230) and involution (INV, D290 and when milk production dropped to 5 kg/day) stages from 9 cows was used for miRNA sequencing. A total of 475 known and 238 novel miRNAs were identified. Fifteen abundantly expressed miRNAs across lactation stages play regulatory roles in basic metabolic, cellular and immunological functions. About 344, 366 and 209 miRNAs were significantly differentially expressed (DE) between GAL and LAC, INV and GAL, and INV and LAC stages, respectively. MiR-29b/miR-363 and miR-874/miR-6254 are important mediators for transition signals from LAC to GAL and from GAL to INV, respectively. Moreover, 58 miRNAs were dynamically DE in all lactation stages and 19 miRNAs were significantly time-dependently DE throughout lactation. Relevant signalling pathways for transition between lactation stages are involved in apoptosis (PTEN and SAPK/JNK), intracellular signalling (protein kinase A, TGF-β and ERK5), cell cycle regulation (STAT3), cytokines, hormones and growth factors (prolactin, growth hormone and glucocorticoid receptor). Overall, our data suggest diverse, temporal and physiological signal-dependent regulatory and mediator functions for miRNAs during lactation.

Comparative Analysis of the miRNome of Bovine Milk Fat, Whey and Cells

Abundant miRNAs have been identified in milk and mammary gland tissues of different species. Typically, RNA in milk can be extracted from different fractions including fat, whey and cells and the mRNA transcriptome of milk could serve as an indicator of the transcriptome of mammary gland tissue. However, it has not been adequately validated if the miRNA transcriptome of any milk fraction could be representative of that of mammary gland tissue. The objectives of this study were to (1) characterize the miRNA expression spectra from three milk fractions- fat, whey and cells; (2) compare miRNome profiles of milk fractions (fat, whey and cells) with mammary gland tissue miRNome, and (3) determine which milk fraction miRNome profile could be a better representative of the miRNome profile of mammary gland tissue. Milk from four healthy Canadian Holstein cows in mid lactation was collected and fractionated. Total RNA extracted from each fraction was used for library preparation followed by small RNA sequencing. In addition, miRNA transcripts of mammary gland tissues from twelve Holstein cows in our previous study were used to compare our data. We identified 210, 200 and 249 known miRNAs from milk fat, whey and cells, respectively, with 188 universally expressed in the three fractions. In addition, 33, 31 and 36 novel miRNAs from milk fat, whey and cells were identified, with 28 common in the three fractions. Among 20 most highly expressed miRNAs in each fraction, 14 were expressed in common and 11 were further shared with mammary gland tissue. The three milk fractions demonstrated a clear separation from each other using a hierarchical cluster analysis with milk fat and whey being most closely related. The miRNome correlation between milk fat and mammary gland tissue (rmean = 0.866) was significantly higher than the other two pairs (p < 0.01), whey/mammary gland tissue (rmean = 0.755) and milk cell/mammary gland tissue (rmean = 0.75), suggesting that milk fat could be an alternative non-invasive source of RNA in assessing miRNA activities in bovine mammary gland. Predicted target genes (1802) of 14 highly expressed miRNAs in milk fractions were enriched in fundamental cellular functions, infection, organ and tissue development. Furthermore, some miRNAs were highly enriched (FDR <0.05) in milk whey (3), cells (11) and mammary gland tissue (14) suggesting specific regulatory functions in the various fractions. In conclusion, we have obtained a comprehensive miRNA profile of the different milk fractions using high throughput sequencing. Our comparative analysis showed that miRNAs from milk fat accurately portrayed the miRNome of mammary gland tissue. Functional annotation of the top expressed miRNAs in milk confirmed their critical regulatory roles in mammary gland functions and potentially to milk recipients.

Milk: an epigenetic amplifier of FTO-mediated transcription? Implications for Western diseases

Single-nucleotide polymorphisms within intron 1 of the FTO (fat mass and obesity-associated) gene are associated with enhanced FTO expression, increased body weight, obesity and type 2 diabetes mellitus (T2DM). The N⁶-methyladenosine (m⁶A) demethylase FTO plays a pivotal regulatory role for postnatal growth and energy expenditure. The purpose of this review is to provide translational evidence that links milk signaling with FTO-activated transcription of the milk recipient. FTO-dependent demethylation of m⁶A regulates mRNA splicing required for adipogenesis, increases the stability of mRNAs, and affects microRNA (miRNA) expression and miRNA biosynthesis. FTO senses branched-chain amino acids (BCAAs) and activates the nutrient sensitive kinase mechanistic target of rapamycin complex 1 (mTORC1), which plays a key role in translation. Milk provides abundant BCAAs and glutamine, critical components increasing FTO expression. CpG hypomethylation in the first intron of FTO has recently been associated with T2DM. CpG methylation is generally associated with gene silencing. In contrast, CpG demethylation generally increases transcription. DNA de novo methylation of CpG sites is facilitated by DNA methyltransferases (DNMT) 3A and 3B, whereas DNA maintenance methylation is controlled by DNMT1. MiRNA-29s target all DNMTs and thus reduce DNA CpG methylation. Cow´s milk provides substantial amounts of exosomal miRNA-29s that reach the systemic circulation and target mRNAs of the milk recipient. Via DNMT suppression, milk exosomal miRNA-29s may reduce the magnitude of FTO methylation, thereby epigenetically increasing FTO expression in the milk consumer. High lactation performance with increased milk yield has recently been associated with excessive miRNA-29 expression of dairy cow mammary epithelial cells (DCMECs). Notably, the galactopoietic hormone prolactin upregulates the transcription factor STAT3, which induces miRNA-29 expression. In a retrovirus-like manner milk exosomes may transfer DCMEC-derived miRNA-29s and bovine FTO mRNA to the milk consumer amplifying FTO expression. There is compelling evidence that obesity, T2DM, prostate and breast cancer, and neurodegenerative diseases are all associated with increased FTO expression. Maximization of lactation performance by veterinary medicine with enhanced miRNA-29s and FTO expression associated with increased exosomal miRNA-29 and FTO mRNA transfer to the milk consumer may represent key epigenetic mechanisms promoting FTO/mTORC1-mediated diseases of civilization.

Mammary epithelial cells isolated from milk are a valuable, non-invasive source of mammary transcripts

Milk is produced in the udder by mammary epithelial cells (MEC). Milk contains MEC, which are gradually exfoliated from the epithelium during lactation. Isolation of MEC from milk using immunomagnetic separation may be a useful non-invasive method to investigate transcriptional regulations in ruminants' udder. This review aims to describe the process of isolating MEC from milk, to provide an overview on the studies that use this method to analyze gene expression by qRT PCR and to evaluate the validity of this method by analyzing and comparing the results between studies. In several goat and cow studies, consistent reductions in alpha-lactalbumin mRNA levels during once-daily milking (ODM) and in SLC2A1 mRNA level during feed restriction are observed. The effect of ODM on alpha-lactalbumin mRNA level was similarly observed in milk isolated MEC and mammary biopsy. Moreover, we and others showed decreasing alpha-lactalbumin and increasing BAX mRNA levels with advanced stages of lactation in dairy cows and buffalo. The relevance of using the milk-isolated MEC method to analyze mammary gene expression is proven, as the transcript variations were also consistent with milk yield and composition variations under the effect of different factors such as prolactin inhibition or photoperiod. However, the RNA from milk-isolated MEC is particularly sensitive to degradation. This could explain the differences obtained between milk-isolated MEC and mammary biopsy in two studies where gene expression was compared using qRT-PCR or RNA Sequencing analyses. As a conclusion, when the RNA quality is conserved, MEC isolated from milk are a valuable, non-invasive source of mammary mRNA to study various factors that impact milk yield and composition (ODM, feeding level, endocrine status, photoperiod modulation, and stage of lactation).