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

Hormonal regulation of miRNA during mammary gland development

The mammary gland is a unique organ as most of its development occurs after birth through stages of proliferation, differentiation and apoptosis that are tightly regulated by circulating hormones and growth factors. Throughout development, hormonal cues induce the regulation of different pathways, ultimately leading to differential transcription and expression of genes involved in this process, but also in the activation or inhibition of post-transcriptional mechanisms of regulation. However, the role of microRNAs (miRNAs) in the different phases of mammary gland remodeling is still poorly understood. The objectives of this study were to analyze the expression of miRNA in key stages of mammary gland development in mice and to determine whether it could be associated with hormonal variation between stages. To do so, miRNAs were isolated from mouse mammary glands at stages of adulthood, pregnancy, lactation and involution, and sequenced. Results showed that 490, 473, 419, and 460 miRNAs are detected in adult, pregnant, lactating and involuting mice, respectively, most of them being common to all four groups, and 58 unique to one stage. Most genes could be divided into six clusters of expression, including two encompassing the highest number of miRNA (clusters 1 and 3) and showing opposite profiles of expression, reaching a peak at adulthood and valley at lactation, or showing the lowest expression at adulthood and peaking at lactation. GO and KEGG analyses suggest that the miRNAs differentially expressed between stages influence the expression of targets associated with mammary gland homeostasis and hormone regulation. To further understand the links between miRNA expression and hormones involved in mammary gland development, miRNAs were then sequenced in breast cells exposed to estradiol, progesterone, prolactin and oxytocin. Four, 38, 24 and 66 miRNAs were associated with progesterone, estradiol, prolactin, and oxytocin exposure, respectively. Finally, when looking at miRNAs modulated by the hormones, differentially expressed during mammary gland development, and having a pattern of expression that could be correlated with the relative levels of hormones known to be found in vivo, 16 miRNAs were identified as likely regulated by circulating hormones. Overall, our study brings a better understanding of the regulation of miRNAs throughout mammary gland development and suggests that there is a relationship between their expression and the main hormones involved in mammary gland development. Future studies will examine this role more in detail.

Identification of circulatory microRNA based biomarkers for early pregnancy diagnosis in buffalo

Introduction: The most crucial factor in improving animal reproduction efficiency is early pregnancy diagnosis. Early diagnosis not only reduces the time interval between two calvings but also aids farmers in identifying open animals, thereby preventing significant milk production losses. Therefore, the objective of this study was to discover circulatory miRNAs that would be useful for early pregnancy diagnosis in buffalo. Material and methods: Blood samples were taken on 0, 6th, 12th, and 18th day after artificial insemination from pregnant animals (n = 30) and non-pregnant animals (n = 20). During these stages of pregnancy, total RNA was extracted, and a small RNA library was subsequently generated and sequenced on the Illumina platform. Subsequently, Real-time PCR was used to validate the findings. Results and discussion: There were 4,022 miRNAs found during the pregnancy, with 15 of those lacking sequences and 4,007 having sequences already in the database. From the beginning of pregnancy until the 18th day, 25 of these miRNAs showed a substantial shift in expression levels in the maternal blood, with a change more than two logs. Furthermore, based on qPCR results, 19 miRNAs were found to be more abundant in pregnant animals than in non-pregnant animals. We used target prediction analysis to learn how maternally expressed miRNAs relate to fetal-maternal communication. In conclusion, miRNA based biomarkers that could be associated with the diagnosis of pregnancy were identified including miR-181a and miR-486 highly upregulated on the 18th day of pregnancy. This study also provides a comprehensive profile of the entire miRNA population in maternal buffalo blood during the early stages of pregnancy.

MicroRNA-2285f regulates milk fat metabolism by targeting MAP2K2 in bovine mammary epithelial cells

In this study, Holstein dairy cows raised in Ningxia were selected as the research object. Mammary epithelial cells (BMECs) were extracted from the milk of eight Holstein cows with significantly different milk fat expression rates and transcribed for sequencing. Bioinformatics analysis was used to analyse the correlation of fat milk percentage, and the critical miR-2285f regulating milk fat was screened out. The target gene binding sites were predicted, and 293T cells and mammary epithelial cells were used as miRNA and target gene models for functional verification in vitro. The tissue difference of miR-2285f Holstein cows was quantitatively analysed by transfecting miR-2285f mimic and inhibitor. Assay (dual luciferase reporter gene assay) and quantitative real-time PCR (quantitative real-time PCR, qRT-PCR), triglyceride (TAG) detection, oil red O detection of lipid droplets, Western Blot assay, Edu and Flow cytometry, The molecular regulatory effects of miR-2285f and target gene MAP2K2 on milk fat metabolism of Holstein dairy cows were studied. The wild-type vector and mutant vector of map2k2-3'utr were constructed, and double luciferase reporting experiments were conducted to verify that MAP2K2 was one of the target genes of miR-2285f. According to qRT-PCR and Western Blot analysis, miR-2285f mainly regulates the expression of MAP2K2 protein in BMECs at the translation level. Bta-miR-2285f can promote cell proliferation and slow cell apoptosis by regulating MAP2K2. Bta-miR-2285f can promote triglyceride (TAG) and lipid droplet accumulation in mammary epithelial cells by targeting MAP2K2. Bta-miR-2285f can regulate protein levels of fat milk marker gene PPARG by targeting MAP2K2. In conclusion, miR-2285f can target the expression of the MAP2K2 gene, promote the proliferation of dairy mammary epithelial cells, inhibit cell apoptosis and regulate the milk fat metabolism in dairy mammary epithelial cells. The results of this study revealed the function of miR-2285f in regulating the differential expression of fat milk in Holstein dairy cows at the cellular level. They provided a theoretical and experimental basis for analysing the regulation network of milk fat synthesis of Holstein dairy cows and the molecular breeding of dairy cows.

Milk exosomal microRNA profiling identified miR-375 and miR-199-5p for regulation of immune response during subclinical mastitis of crossbred cattle

BACKGROUND

The dairy industry has experienced significant economic losses as a result of mastitis, an inflammatory disease of cows, including both subclinical and clinical cases. Milk exosome microRNAs have gained attention due to their stable and selective wrapping nature, offering potential for the prognosis and diagnosis of bovine mastitis, the most common pathological condition of the mammary gland.

METHODS AND RESULTS

 In the present investigation, the microRNA profile of milk exosomes was explored using high-throughput small RNA sequencing data in sub-clinical mastitic and healthy crossbred Vrindavani cattle. In both groups, 349 microRNAs were identified, with 238 (68.19%) microRNAs co-expressed; however, 35 and 76 distinct microRNAs were found in subclinical mastitic and healthy cattle, respectively. Differential expression analysis revealed 11 microRNAs upregulated, and 18 microRNAs were downregulated in sub-clinical mastitic cattle. The functional annotation of the target genes of differentially expressed known and novel microRNAs including bta-miR-375, bta-miR-199-5p and bta-miR-12030 reveals their involvement in the regulation of immune response and inflammatory mechanisms and could be involved in development of mastitis.

CONCLUSIONS

The analysis of milk exosomal miRNAs cargos hold great promise as an approach to study the underlying molecular mechanisms associated with mastitis in high milk producing dairy cattle. Concurrently, the significantly downregulated miR-375 may upregulate key target genes, including CTLA4, IHH, IRF1, and IL7R. These genes are negative regulators of immune response pathways, which could be associated with impaired inflammatory mechanisms in mammary cells. According to the findings, bta-miR-375 could be a promising biomarker for the development of mastitis in dairy cattle.

Identification of Hub Genes and Target miRNAs Crucial for Milk Production in Holstein Friesian Dairy Cattle

Dairy milk production is a quantitative trait that is controlled by many biological and environmental factors. This study employs a network-driven systems approach and clustering algorithm to uncover deeper insights into its genetic associations. We analyzed the GSE33680 dataset from the GEO database to understand the biological importance of milk production through gene expression and modules. In this study, we employed CytoNCA and ClusterONE plugins within Cytoscape for network analysis. Moreover, miRWalk software was utilized to detect miRNAs, and DAVID was employed to identify gene ontology and pathways. The results revealed 140 up-regulated genes and 312 down-regulated genes. In addition, we have identified 91 influential genes and 47 miRNAs that are closely associated with milk production. Through our examination of the network connecting these genes, we have found significant involvement in important biological processes such as calcium ion transit across cell membranes, the BMP signaling pathway, and the regulation of MAPK cascade. The conclusive network analysis further reveals that GAPDH, KDR, CSF1, PYGM, RET, PPP2CA, GUSB, and PRKCA are closely linked to key pathways essential for governing milk production. Various mechanisms can control these genes, making them valuable for breeding programs aiming to enhance selection indexes.

Pigeon during the Breeding Cycle: Behaviors, Composition and Formation of Crop Milk, and Physiological Adaptation

Pigeon is an important economic poultry species in many countries. As an altricial bird, its growth and development are largely reliant on pigeon milk produced by the crop tissue in the first week. During the breeding cycle, pigeons undergo a series of behavioral changes. Pigeon milk is generally characterized by having high concentrations of proteins and lipids, and a complicated regulatory network is involved in the milk formation. Hormones, especially prolactin, could promote the proliferation of crop epidermal cells and nutrient accumulation. The expression of target genes associated with these important biological processes in the crop epidermis is affected by non-coding RNAs. Meanwhile, signaling pathways, such as target of rapamycin (TOR), Janus kinase/signal transducer and activator of transcription proteins (JAK/STAT), protein kinase B (Akt), etc., influence the production of crop milk by either enhancing protein synthesis in crop cells or inducing apoptosis of crop epidermal cells. In order to adapt to the different breeding periods, pigeons are physiologically changed in their intestinal morphology and function and liver metabolism. This paper reviews the behaviors and physiological adaptations of pigeon during the breeding cycle, the composition of pigeon crop milk, and the mechanism of its formation, which is important for a better understanding of the physiology of altricial birds and the development of artificial crop milk.

Identification and Characterization of circRNAs in Non-Lactating Dairy Goat Mammary Glands Reveal Their Regulatory Role in Mammary Cell Involution and Remodeling

This study conducted transcriptome sequencing of goat-mammary-gland tissue at the late lactation (LL), dry period (DP), and late gestation (LG) stages to reveal the expression characteristics and molecular functions of circRNAs during mammary involution. A total of 11,756 circRNAs were identified in this study, of which 2528 circRNAs were expressed in all three stages. The number of exonic circRNAs was the largest, and the least identified circRNAs were antisense circRNAs. circRNA source gene analysis found that 9282 circRNAs were derived from 3889 genes, and 127 circRNAs' source genes were unknown. Gene Ontology (GO) terms, such as histone modification, regulation of GTPase activity, and establishment or maintenance of cell polarity, were significantly enriched (FDR < 0.05), which indicates the functional diversity of circRNAs' source genes. A total of 218 differentially expressed circRNAs were identified during the non-lactation period. The number of specifically expressed circRNAs was the highest in the DP and the lowest in LL stages. These indicated temporal specificity of circRNA expression in mammary gland tissues at different developmental stages. In addition, this study also constructed circRNA-miRNA-mRNA competitive endogenous RNA (ceRNA) regulatory networks related to mammary development, immunity, substance metabolism, and apoptosis. These findings help understand the regulatory role of circRNAs in mammary cell involution and remodeling.

MicroRNAs in Ruminants and Their Potential Role in Nutrition and Physiology

The knowledge of how diet choices, dietary supplements, and feed intake influence molecular mechanisms in ruminant nutrition and physiology to maintain ruminant health, is essential to attain. In the present review, we focus on the role of microRNAs in ruminant health and disease; additionally, we discuss the potential of circulating microRNAs as biomarkers of disease in ruminants and the state of technology for their detection, also considering the major difficulties in the transition of biomarker development from bench to clinical practice. MicroRNAs are an inexhaustible class of endogenous non-protein coding small RNAs of 18 to 25 nucleotides that target either the 3' untranslated (UTR) or coding region of genes, ensuring a tight post-transcriptionally controlled regulation of gene expression. The development of new "omics" technologies facilitated a fresh perspective on the nutrition-to-gene relationship, incorporating more extensive data from molecular genetics, animal nutrition, and veterinary sciences. MicroRNAs might serve as important regulators of metabolic processes and may present the inter-phase between nutrition and gene regulation, controlled by the diet. The development of biomarkers holds the potential to revolutionize veterinary practice through faster disease detection, more accurate ruminant health monitoring, enhanced welfare, and increased productivity. Finally, we summarize the latest findings on how microRNAs function as biomarkers, how technological paradigms are reshaping this field of research, and how platforms are being used to identify novel biomarkers. Numerous studies have demonstrated a connection between circulating microRNAs and ruminant diseases such as mastitis, tuberculosis, foot-and-mouth disease, fasciolosis, and metabolic disorders. Therefore, the identification and analysis of a small number of microRNAs can provide crucial information about the stage of a disease, etiology, and prognosis.

Comprehensive Profiling of ceRNA (circRNA-miRNA-mRNA) Networks in Hypothalamic-Pituitary-Mammary Gland Axis of Dairy Cows under Heat Stress

Heat stress (HS) is directly correlated with mammary gland dysfunction and the hypothalamic-pituitary-mammary gland (HPM) axis is involved in regulating stress responses and lactation in dairy cows. Circular RNAs (circRNAs) play major roles in regulating transcription and post-transcription but their expression in the HPM axis of dairy cows under HS is still unclear. In the present study, we performed RNA sequencing to identify diferentially expressed (DE) circRNAs, DE microRNAs(miRNAs) and DEmRNAs, and performed bioinformatics analysis on those in HPM axis-related tissues of heat-stressed and normal cows. A total of 1680, 1112 and 521 DEcircRNAs, 120, 493 and 108 DEmiRNAs, 274, 6475 and 3134 DEmRNAs were identified in the hypothalamic, pituitary, and mammary gland tissues, respectively. Gene ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG) enrichment analyses indicated that the MAPK signaling pathway is potentially a key pathway. Competitive endogenous RNA (ceRNA) networks related to HS response and lactation regulation were established in three tissues. In conclusion, our results indicate that HS induces differential circRNA expression profiles in HPM axis-related tissues, and the predicted ceRNA network provides a molecular basis for regulating the stress response and lactation regulation in heat-stressed dairy cows.

miR-193-5p negatively regulates PIK3CD to promote crop fibrocyte proliferation in pigeon (Columba livia)

The crop of pigeon has specific characteristics as producing crop milk in the lactating period. However, the exact mechanisms underlying the regulation of crop lactation remain unclear. miRNAs, the essential regulators of gene expression, are implicated in various physiological and biological activities. In this study, we discovered a new miRNA that regulated phosphatidylinositol-4,5-bisphosphate 3-kinase catalytic subunit delta (PIK3CD) and crop fibrocyte proliferation. Results of the luciferase reporter assay suggested that miR-193-5p suppressed PIK3CD expression by targeting a conserved binding site in the 3'-untranslated region (UTR) of PIK3CD mRNA. MiR-193-5p promoted crop fibrocyte proliferation and migration, whereas PIK3CD inhibited these effects. These findings suggested an important regulatory role of miR-193-5p in crop fibrocyte proliferation, suggesting that miR-193-5p and PIK3CD might be important regulators of crop milk production.

Non-coding transcriptomic profiles in the sheep mammary gland during different lactation periods

Sheep milk production is a dynamic and multifactorial trait regulated by diverse biological mechanisms. To improve the quality and production of sheep milk, it is necessary to understand the underlying non-coding transcriptomic mechanisms. In this study, ribonucleic acid-sequencing (RNA-seq) was used to profile the expression of microRNAs (miRNAs) and circular RNAs (circRNAs) in the sheep mammary gland at three key lactation time points (perinatal period, PP; early lactation, EL; and peak lactation, PL). A total of 2,369 novel circRNAs and 272 miRNAs were profiled, of which 348, 373, and 36 differentially expressed (DE) circRNAs and 30, 34, and 7 DE miRNAs were detected in the comparison of EL vs. PP, PL vs. PP, and PL vs. EL, respectively. A series of bioinformatics analyses including functional enrichment, machine learning prediction, and competing endogenous RNA (ceRNA) network analyses were conducted to identify subsets of the potential candidate miRNAs (e.g., oar_miR_148a, oar_miR_362, and oar_miR_432) and circRNAs (e.g., novel_circ_0011066, novel_circ_0010460, and novel_circ_0006589) involved in sheep mammary gland development. Taken together, this study offers a window into the dynamics of non-coding transcriptomes that occur during sheep lactation and may provide further insights into miRNA and circRNA that influence sheep mammary gland development.

miRNA expression patterns in blood leukocytes and milk somatic cells of goats infected with small ruminant lentivirus (SRLV)

The study aims to determine the selected miRNAs expression in milk somatic cells (MSC) and blood leukocytes (BL) of SRLV-seronegative (SRLV-SN) and SRLV-seropositive (SRLV-SP) goats. A functional in silico analysis of their target genes was also conducted. MiR-93-5p and miR-30e-5p were expressed only in BL, while miR-144 was expressed only in MSC, regardless of SRLV infection. In the SRLV-SP goats, higher miR-214-3p and miR-221-5p levels were found in the MSC than in the BL. Only miR-30e-5p was influenced by the lactation stage in BL in both groups, while only miR-93-5p was altered in BL of SRLV-SN goats. The target gene protein products exhibited contradictory functions, protecting the host from virus on the one hand and assisting viruses in their life cycle on the other. The differential expression of the miRNAs observed between the MSC and BL of SRLV-SP goats may suggest that the local immune response to the infection in the udder differs from the systemic response, and acts independently. Some miRNAs demonstrated different expression between lactation stages. It may be influenced by the metabolic burden occurring in early lactation and its peak. Some of the studied miRNAs may influence viral infection by regulating the expression of their target genes.

Effect of Diet Enriched With Hemp Seeds on Goat Milk Fatty Acids, Transcriptome, and miRNAs

In dairy ruminants, a diet supplemented with feed rich in unsaturated fatty acids can be an effective medium to increase the health-promoting properties of milk, although their effect on the pathways/genes involved in these processes has not been properly and completely defined to date. To improve our knowledge of the cell’s activity in specific conditions, next-generation RNA-sequencing technology was used to allow whole transcriptome characterization under given conditions. In addition to this, microRNAs (miRNAs) have recently been known as post-transcriptional regulators in fatty acid and cholesterol metabolism by targeting lipid metabolism genes. In this study, to analyze the transcriptome and miRNAs in goat milk after a supplemental diet enriched with linoleic acid (hemp seeds), next-generation RNA-sequencing was used in order to point out the general biological mechanisms underlying the effects related to milk fat metabolism. Ten pluriparous Alpine goats were fed with the same pretreatment diet for 40 days; then, they were arranged to two dietary treatments consisting of control (C) and hemp seed (H)-supplemented diets. Milk samples were collected at 40 (time point = T0) and 140 days of lactation (time point = T1). Milk fatty acid (FA) profiles revealed a significant effect of hemp seeds that determined a strong increment in the preformed FA, causing a reduction in the concentration of de-novo FA. Monounsaturated and polyunsaturated n−3 FAs were increased by hemp treatment, determining a reduction in the n−6/n−3 ratio. After removing milk fats and proteins, RNA was extracted from the milk cells and transcriptomic analysis was conducted using Illumina RNA-sequencing. A total of 3,835 genes were highly differentially expressed (p-value &lt; 0.05, fold change &gt; 1.5, and FDR &lt; 0.05) in the H group. Functional analyses evidenced changes in metabolism, immune, and inflammatory responses. Furthermore, modifications in feeding strategies affected also key transcription factors regulating the expression of several genes involved in milk fat metabolism, such as peroxisome proliferator-activated receptors (PPARs). Moreover, 38 (15 known and 23 novel) differentially expressed miRNAs were uncovered in the H group and their potential functions were also predicted. This study gives the possibility to improve our knowledge of the molecular changes occurring after a hemp seed supplementation in the goat diet and increase our understanding of the relationship between nutrient variation and phenotypic effects.