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

Profiling several key milk miRNAs and analysing their signalling pathways in dairy sheep breeds during peak and late lactation

BACKGROUND

The comprehensive understanding of microRNAs (miRNAs) in sheep milk during various lactation periods and their impact on milk yield and composition remains limited.

OBJECTIVES

This study aimed to investigate the expression patterns of four highly expressed miRNAs in sheep milk and their association with milk composition and yield parameters during peak and late lactation stages.

METHODS

A total of 40 healthy 4-year-old Akkaraman (n = 20) and Awassi (n = 20) ewes registered with the Ministry of Agriculture and Forestry of the Republic of Türkiye were used in the present study. For miRNA isolation from milk, the Qiagen miRNeasy Serum/Plasma Advanced Kit was utilised following the manufacturer's instructions. The expression levels of miRNAs were assessed using Qiagen miRNA PCR Assays.

RESULTS

The significant fold changes in the expression levels of oar-miR-30a-5p, oar-miR-148a and oar-miR-181a were observed between peak and late lactation periods in the Awassi sheep breed. Conversely, only oar-miR-30a-5p and oar-miR-148a exhibited statistically significant changes in the Akkaraman sheep breed during the same lactation periods. Furthermore, oar-miR-21-5p demonstrated a significant fold change exclusively in peak lactation compared to Akkaraman and Awassi ewes.

CONCLUSIONS

The findings suggest that the expression of the analysed miRNAs is influenced by both the lactation stage and different sheep breeds. This study offers valuable insights into the relationship between key miRNA expressions in sheep milk and milk composition and yield parameters during peak and late lactation, contributing to the existing knowledge in this field.

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.

Comparison of miR-106b, miR-191, and miR-30d expression dynamics in milk with regard to its composition in Holstein and Ayrshire cows

OBJECTIVE

Milk composition varies considerably and depends on paratypical, genetic, and epigenetic factors. MiRNAs belong to the class of small non-coding RNAs; they are one of the key tools of epigenetic control because of their ability to regulate gene expression at the post-transcriptional level. We compared the relative expression levels of miR-106b, miR-191, and miR-30d in milk to demonstrate the relationship between the content of these miRNAs with protein and fat components of milk in Holstein and Ayrshire cattle.

METHODS

Milk fat, protein, and casein contents were determined in the obtained samples, as well as the content of the main fatty acids (g/100 g milk), including: saturated acids, such as myristic (C14:0), palmitic (C16:0), and stearic (C18:0) acids; monounsaturated acids, including oleic (C18:1) acid; as well as long-, medium- and short-chain, polyunsaturated, and trans fatty acids. Real-time stem-loop one-tube reverse transcription polymerase chain reaction with TaqMan probes was used to measure the miRNA expression levels.

RESULTS

The miRNA expression levels in milk samples were found to be decreased in the first two months in Holstein breed, and in the first four months in Ayrshire breed. Correlation analysis did not reveal any dependence between changes in the expression level of miRNA and milk fat content, but showed a multidirectional relationship with individual milk fatty acids. Positive associations between the expression levels of miR-106b and miR-30d and protein and casein content were found in the Ayrshire breed. Receiver operating characteristic curve analysis showed that miR-106b and miR-30d expression levels can cause changes in fatty acid and protein composition of milk in Ayrshire cows, whereas miR-106b expression level determines the fatty acid composition in Holsteins.

CONCLUSION

The data obtained in this study showed that miR-106b, miR-191, and miR-30d expression levels in milk samples have peculiarities associated with breed affiliation and the lactation period.

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.

Gestational age at birth influences protein and RNA content in human milk extracellular vesicles

Human milk extracellular vesicles (HM EVs) are proposed to protect against disease development in infants. This protection could in part be facilitated by the bioactive EV cargo of proteins and RNA. Notably, mothers birth infants of different gestational ages with unique needs, wherein the EV cargo of HM may diverge. We collected HM from lactating mothers within two weeks of a term or preterm birth. Following purification of EVs, proteins and mRNA were extracted for proteomics and sequencing analyses, respectively. Over 2000 protein groups were identified, and over 8000 genes were quantified. The total number of proteins and mRNA did not differ significantly between the two conditions, while functional bioinformatics of differentially expressed cargo indicated enrichment in immunoregulatory cargo for preterm HM EVs. In term HM EVs, significantly upregulated cargo was enriched in metabolism-related functions. Based on gene expression signatures from HM-contained single cell sequencing data, we proposed that a larger portion of preterm HM EVs are secreted by immune cells, whereas term HM EVs contain more signatures of lactocyte epithelial cells. Proposed differences in EV cargo could indicate variation in mother's milk based on infants' gestational age and provide basis for further functional characterisation.

MiRNome variations in milk fractions during feed restrictions of different intensities in dairy cows

BACKGROUND

In dairy cows, diet is one factor that can affect their milk production and composition. However, the effect of feed restriction on milk miRNome has not yet been described. Indeed, milk is the body fluid with the highest RNA concentration, which includes numerous microRNA. Its presence in the four different milk fractions, whole milk, fat globules, mammary epithelial cells and extracellular vesicles, is still poorly documented. This study aimed to describe the effects of different feed restrictions on the miRNome composition of different milk fractions.

RESULTS

Two feed restrictions were applied to lactating dairy cows, one of high intensity and one of moderate intensity. 2,896 mature microRNA were identified in the different milk fractions studied, including 1,493 that were already known in the bovine species. Among the 1,096 microRNA that were sufficiently abundant to be informative, the abundance of 1,027 of them varied between fractions: 36 of those were exclusive to one milk fraction. Feed restriction affected the abundance of 155 microRNA, with whole milk and milk extracellular vesicles being the most affected, whereas milk fat globules and exfoliated mammary epithelial cells were little or not affected at all. The high intensity feed restriction led to more microRNA variations in milk than moderate restriction. The target prediction of known microRNA that varied under feed restriction suggested the modification of some key pathways for lactation related to milk fat and protein metabolisms, cell cycle, and stress responses.

CONCLUSIONS

This study highlighted that the miRNome of each milk fraction is specific, with mostly the same microRNA composition but with variations in abundance between fractions. These specific miRNomes were affected differently by feed restrictions, the intensity of which appeared to be a major factor modulating milk miRNomes. These findings offer opportunities for future research on the use of milk miRNA as biomarkers of energy status in dairy cows, which is affected by feed restrictions.

Transcriptome analysis of the Bactrian camel (Camelus bactrianus) reveals candidate genes affecting milk production traits

BACKGROUND

Milk production traits are complex traits with vital economic importance in the camel industry. However, the genetic mechanisms regulating milk production traits in camels remain poorly understood. Therefore, we aimed to identify candidate genes and metabolic pathways that affect milk production traits in Bactrian camels.

METHODS

We classified camels (fourth parity) as low- or high-yield, examined pregnant camels using B-mode ultrasonography, observed the microscopic changes in the mammary gland using hematoxylin and eosin (HE) staining, and used RNA sequencing to identify differentially expressed genes (DEGs) and pathways.

RESULTS

The average standard milk yield over the 300 days during parity was recorded as 470.18 ± 9.75 and 978.34 ± 3.80 kg in low- and high-performance camels, respectively. Nine female Junggar Bactrian camels were subjected to transcriptome sequencing, and 609 and 393 DEGs were identified in the low-yield vs. high-yield (WDL vs. WGH) and pregnancy versus colostrum period (RSQ vs. CRQ) comparison groups, respectively. The DEGs were compared with genes associated with milk production traits in the Animal Quantitative Trait Loci database and in Alashan Bactrian camels, and 65 and 46 overlapping candidate genes were obtained, respectively. Functional enrichment and protein-protein interaction network analyses of the DEGs and candidate genes were conducted. After comparing our results with those of other livestock studies, we identified 16 signaling pathways and 27 core candidate genes associated with maternal parturition, estrogen regulation, initiation of lactation, and milk production traits. The pathways suggest that emerged milk production involves the regulation of multiple complex metabolic and cellular developmental processes in camels. Finally, the RNA sequencing results were validated using quantitative real-time PCR; the 15 selected genes exhibited consistent expression changes.

CONCLUSIONS

This study identified DEGs and metabolic pathways affecting maternal parturition and milk production traits. The results provides a theoretical foundation for further research on the molecular mechanism of genes related to milk production traits in camels. Furthermore, these findings will help improve breeding strategies to achieve the desired milk yield in camels.

CircRNA-02191 regulating unsaturated fatty acid synthesis by adsorbing miR-145 to enhance CD36 expression in bovine mammary gland

CD36 functions as a receptor for long-chain fatty acids, promoting the absorption and transport of long-chain unsaturated fatty acids. However, the regulatory influence of upstream circRNAs or miRNAs on its expression in cow mammary gland remains unclear. Herein, we performed high-throughput sequencing to screen for differentially expressed miRNAs and mRNAs in bovine mammary tissue during the late-lactation and the dry period to screen and conducted bioinformatics analysis to identify 420 miRNA/mRNA pairs, including miR-145/CD36. Experimental results indicate that miR-145 can directly target CD36 and inhibit its expression. Additionally, the circRNA-02191 sequence is predicted to contain a miR-145 binding site. As shown by dual luciferase reporter system detection, circRNA-02191 bound to miR-145 and its overexpression significantly reduced the expression of miR-145. Furthermore, the overexpression of miR-145 inhibited triglyceride accumulation, while circRNA-02191 enhanced the expression of the miR-145 target gene CD36. The above results indicate that circRNA-02191 can regulate triglyceride and fatty acid components by binding miR-145 and subsequently alleviating the inhibitory effect of miR-145 on the expression of CD36. Taken together, these findings present a novel approach to improve milk quality by analyzing the regulatory effect and mechanism regulating the circ02191/miR-145/CD36 pathway on fatty acid synthesis in the mammary gland of dairy cows.

Analysis of immune-related microRNAs in cows and newborn calves

Bovine colostrum contains a high concentration of immune-related microRNAs (miRNAs) that are packaged in exosomes and are very stable. In this study, 5 immune-related miRNAs (miR-142-5p, miR-150, miR-155, miR-181a, and miR-223) were quantified in dam blood, colostrum, and calf blood using reverse transcription quantitative PCR. Their levels in calf blood after colostrum ingestion were investigated to assess whether miRNAs are transferred from the dam to newborn calves. Three groups of Holstein-Friesian bull calves were bottle-fed 2 L of colostrum or milk from different sources twice per day. The group A calves received colostrum from their own dam and the group B calves were fed foster dam colostrum. Each pair of group A and group B calves were fed identical colostrum from the same milking of the corresponding group A dam for 3 d and then bulk tank milk for 7 d after birth. Group C calves were fed only 2L of "pooled colostrum" from multiple dams d 0 to 4 postpartum, and then fed bulk tank milk thereafter for 7 d after birth. The groups were fed colostrum from different sources and different amounts to assess possible miRNA absorption from the colostrum. All miRNAs were at the highest level in colostrum at d 0 and then decreased rapidly after d 1. The level of miR-150 had the largest decrease from 489 × 106 copies/µL (d 0) to 78 × 106 copies/µL (d 1). MicroRNA-223 and miR-155 were the most abundant in both colostrum and milk. Dam colostrum had significantly higher levels of miR-142-5p, miR-155, and miR-181a than the bulk tank milk. However, only the miR-155 concentration was significantly higher in the dam colostrum than in the pooled colostrum. The concentrations of miRNAs in the colostrum were less than in the cow blood (100- to 1,000-fold less). There was no significant correlation between the level of miRNAs in the dam blood and their colostrum, suggesting that miRNA is synthesized locally by the mammary gland rather than being transferred from the blood. MicroRNA-223 had the highest level in both calf and cow blood compared with the other 4 immune-related miRNAs. Calves were born with high levels of immune-related miRNAs in their blood, and there were no significant differences in miRNA levels between the 3 calf groups at birth or after they were fed different colostrum. This suggests that these miRNAs were not transferred from the colostrum to the newborn calves.

Oar-miR-16b and oar-miR-27a: negatively correlated with milk yield and milk protein in sheep

The characterization of miRNAs from sheep milk and their effect on milk yield and composition in sheep are remain unclear. Therefore, the aim of this study was to determine the expression pattern of several important miRNAs, which are associated with lactation in the sheep milk between high and low lactating-yield ewe groups. In addition to experimentally obtained miRNA expression results, the miRNA target genes were determined by bioinformatics analysis to identify biological pathways involved. miRNAs found to differ significantly in the expression level between the groups were oar-miR-181a, oar-miR-23a, oar-miR-27a, oar-miR-16b and oar-miR-374. Also, oar-miR-27a was shown negative correlation with milk protein and lactose contents while oar-miR-16b was shown negative correlation with milk yield in the high milk yield group. The highest connected hub genes for miR-27a target genes were determined as MAPK14 and PPARG. Also, six genes (HSPA4L, DNAJA2, ATP6V1B2, PPP2R1A, PPP2R1B, and PRKAR2A) were detected as hub genes for miR-16b. In this study, the relationship between expression profiles of several important miRNAs in sheep milk and milk yield and milk composition were investigated for the first time in high and low lactating yield groups.

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.

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.

Multi-omic data integration for the study of production, carcass, and meat quality traits in Nellore cattle

Data integration using hierarchical analysis based on the central dogma or common pathway enrichment analysis may not reveal non-obvious relationships among omic data. Here, we applied factor analysis (FA) and Bayesian network (BN) modeling to integrate different omic data and complex traits by latent variables (production, carcass, and meat quality traits). A total of 14 latent variables were identified: five for phenotype, three for miRNA, four for protein, and two for mRNA data. Pearson correlation coefficients showed negative correlations between latent variables miRNA 1 (mirna1) and miRNA 2 (mirna2) (−0.47), ribeye area (REA) and protein 4 (prot4) (−0.33), REA and protein 2 (prot2) (−0.3), carcass and prot4 (−0.31), carcass and prot2 (−0.28), and backfat thickness (BFT) and miRNA 3 (mirna3) (−0.25). Positive correlations were observed among the four protein factors (0.45–0.83): between meat quality and fat content (0.71), fat content and carcass (0.74), fat content and REA (0.76), and REA and carcass (0.99). BN presented arcs from the carcass, meat quality, prot2, and prot4 latent variables to REA; from meat quality, REA, mirna2, and gene expression mRNA1 to fat content; from protein 1 (prot1) and mirna2 to protein 5 (prot5); and from prot5 and carcass to prot2. The relations of protein latent variables suggest new hypotheses about the impact of these proteins on REA. The network also showed relationships among miRNAs and nebulin proteins. REA seems to be the central node in the network, influencing carcass, prot2, prot4, mRNA1, and meat quality, suggesting that REA is a good indicator of meat quality. The connection among miRNA latent variables, BFT, and fat content relates to the influence of miRNAs on lipid metabolism. The relationship between mirna1 and prot5 composed of isoforms of nebulin needs further investigation. The FA identified latent variables, decreasing the dimensionality and complexity of the data. The BN was capable of generating interrelationships among latent variables from different types of data, allowing the integration of omics and complex traits and identifying conditional independencies. Our framework based on FA and BN is capable of generating new hypotheses for molecular research, by integrating different types of data and exploring non-obvious relationships.

Deciphering the Key Regulatory Roles of KLF6 and Bta-miR-148a on Milk Fat Metabolism in Bovine Mammary Epithelial Cells

MicroRNAs (miRNAs) are non-coding RNAs that regulate the expression of their target genes involved in many cellular functions at the post-transcriptional level. Previously, bta-miR-148a showed significantly high expression in bovine mammary epithelial cells (BMECs) of Chinese Holstein cows producing high milk fat compared to those with low milk fat content. Here, we investigated the role of bta-miR-148a through targeting Krüppel-like factor 6 (KLF6) and further analyzed the role of KLF6 in regulating fat metabolism through targeting PPARA, AMPK/mTOR/PPARG, and other fat marker genes in BMECs of Chinese Holstein. The bioinformatics analysis showed that the 3’ UTR of KLF6 mRNA possesses the binding sites for bta-miR-148a, which was further verified through dual-luciferase reporter assay. The BMECs were transfected with bta-miR-148a-mimic, inhibitor, and shNC, and the expression of KLF6 was found to be negatively regulated by bta-miR-148a. Moreover, the contents of triglyceride (TG), and cholesterol (CHO) in BMECs transfected with bta-miR-148a-mimic were significantly lower than the contents in BMECs transfected with bta-miR-148a-shNC. Meanwhile, the TG and CHO contents were significantly increased in BMECs transfected with bta-miR-148a-inhibitor than in BMECs transfected with bta-miR-148a-shNC. In addition, the TG and CHO contents were significantly decreased in BMECs upon the down-regulation of KLF6 through transfection with pb7sk-KLF6-siRNA1 compared to the control group. Contrarily, when KLF6 was overexpressed in BMECs through transfection with pBI-CMV3-KLF6, the TG and CHO contents were significantly increased compared to the control group. Whereas, the qPCR and Western blot evaluation of PPARA, AMPK/mTOR/PPARG, and other fat marker genes revealed that all of the genes were considerably down-regulated in the KLF6-KO-BMECs compared to the normal BMECs. Taking advantage of deploying new molecular markers and regulators for increasing the production of better-quality milk with tailored fat contents would be the hallmark in dairy sector. Hence, bta-miR-148a and KLF6 are potential candidates for increased milk synthesis and the production of valuable milk components in dairy cattle through marker-assisted selection in molecular breeding. Furthermore, this study hints at the extrapolation of a myriad of functions of other KLF family members in milk fat synthesis.

Emerging Roles of Noncoding RNAs in Bovine Mastitis Diseases

Non-coding RNAs (ncRNAs) are an abundant class of RNA with varying nucleotide lengths. They have been shown to have great potential in eutherians/human disease diagnosis and treatments and are now gaining more importance for the improvement of diseases in livestock. To date, thousands of ncRNAs have been discovered in the bovine genome and the continuous advancement in deep sequencing technologies and various bioinformatics tools has enabled the elucidation of their roles in bovine health. Among farm animals' diseases, mastitis, a common inflammatory disease in cattle, has caused devastating economic losses to dairy farmers over the last few decades. Here, we summarize the biology of bovine mastitis and comprehensively discuss the roles of ncRNAs in different types of mastitis infection. Based on our findings and relevant literature, we highlighted various evidence of ncRNA roles in mastitis. Different approaches (in vivo versus in vitro) for exploring ncRNA roles in mastitis are emphasized. More particularly, the potential applications of emerging genome editing technologies, as well as integrated omics platforms for ncRNA studies and implications for mastitis are presented.

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.

Emerging Roles of Non-Coding RNAs in the Feed Efficiency of Livestock Species

A global population of already more than seven billion people has led to an increased demand for food and water, and especially the demand for meat. Moreover, the cost of feed used in animal production has also increased dramatically, which requires animal breeders to find alternatives to reduce feed consumption. Understanding the biology underlying feed efficiency (FE) allows for a better selection of feed-efficient animals. Non-coding RNAs (ncRNAs), especially micro RNAs (miRNAs) and long non-coding RNAs (lncRNAs), play important roles in the regulation of bio-logical processes and disease development. The functions of ncRNAs in the biology of FE have emerged as they participate in the regulation of many genes and pathways related to the major FE indicators, such as residual feed intake and feed conversion ratio. This review provides the state of the art studies related to the ncRNAs associated with FE in livestock species. The contribution of ncRNAs to FE in the liver, muscle, and adipose tissues were summarized. The research gap of the function of ncRNAs in key processes for improved FE, such as the nutrition, heat stress, and gut–brain axis, was examined. Finally, the potential uses of ncRNAs for the improvement of FE were discussed.

Early downregulation of hsa-miR-144-3p in serum from drug-naïve Parkinson's disease patients

Advanced age represents one of the major risk factors for Parkinson's Disease. Recent biomedical studies posit a role for microRNAs, also known to be remodelled during ageing. However, the relationship between microRNA remodelling and ageing in Parkinson's Disease, has not been fully elucidated. Therefore, the aim of the present study is to unravel the relevance of microRNAs as biomarkers of Parkinson's Disease within the ageing framework. We employed Next Generation Sequencing to profile serum microRNAs from samples informative for Parkinson's Disease (recently diagnosed, drug-naïve) and healthy ageing (centenarians) plus healthy controls, age-matched with Parkinson's Disease patients. Potential microRNA candidates markers, emerging from the combination of differential expression and network analyses, were further validated in an independent cohort including both drug-naïve and advanced Parkinson's Disease patients, and healthy siblings of Parkinson's Disease patients at higher genetic risk for developing the disease. While we did not find evidences of microRNAs co-regulated in Parkinson's Disease and ageing, we report that hsa-miR-144-3p is consistently down-regulated in early Parkinson's Disease patients. Moreover, interestingly, functional analysis revealed that hsa-miR-144-3p is involved in the regulation of coagulation, a process known to be altered in Parkinson's Disease. Our results consistently show the down-regulation of hsa-mir144-3p in early Parkinson's Disease, robustly confirmed across a variety of analytical and experimental analyses. These promising results ask for further research to unveil the functional details of the involvement of hsa-mir144-3p in Parkinson's Disease.

Characterization of RNA Editome in the Mammary Gland of Yaks during the Lactation and Dry Periods

Simple Summary

In order to study the influence of RNA editing sites on lactation and mammary gland development process in yaks, we comprehensively characterized the RNA editome of the yak mammary gland during the lactation period and dry period by using the transcriptome and genome sequencing data. The results revealed 82,872 nonredundant RNA editing sites, 14,159 of which were differentially edited between the lactation period and dry period. Enrichment analysis showed that the genes harboring differential editing sites were mainly associated with mammary gland development-related pathways, such as MAPK pathway, PI3K-Akt pathway, FoxO signaling pathway, GnRH signaling pathway, and focal adhesion pathway. Our findings offer some novel insights into the RNA editing function in the mammary gland of yaks.

Abstract

The mammary gland is a complicated organ comprising several types of cells, and it undergoes extensive morphogenetic and metabolic changes during the female reproductive cycle. RNA editing is a posttranscriptional modification event occurring at the RNA nucleotide level, and it drives transcriptomic and proteomic diversities, with potential functional consequences. RNA editing in the mammary gland of yaks, however, remains poorly understood. Here, we used REDItools to identify RNA editing sites in mammary gland tissues in yaks during the lactation period (LP, n = 2) and dry period (DP, n = 3). Totally, 82,872 unique RNA editing sites were identified, most of which were detected in the noncoding regions with a low editing degree. In the coding regions (CDS), we detected 5235 editing sites, among which 1884 caused nonsynonymous amino acid changes. Of these RNA editing sites, 486 were found to generate novel possible miRNA target sites or interfere with the initial miRNA binding sites, indicating that RNA editing was related to gene regulation mediated by miRNA. A total of 14,159 RNA editing sites (involving 3238 common genes) showed a significant differential editing level in the LP when compared with that in the DP through Tukey’s Honest Significant Difference method (p < 0.05). According to the Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, genes that showed different RNA editing levels mainly participated in pathways highly related to mammary gland development, including MAPK, PI3K-Akt, FoxO, and GnRH signaling pathways. Collectively, this work demonstrated for the first time the dynamic RNA editome profiles in the mammary gland of yaks and shed more light on the mechanism that regulates lactation together with mammary gland development.

Longitudinal Human Milk miRNA Composition over the First 3 mo of Lactation in a Cohort of Healthy Mothers Delivering Term Infants

BACKGROUND

MicroRNAs (miRNAs) are small noncoding RNAs involved in posttranscriptional regulation. miRNAs can be secreted and found in many body fluids, and although they are particularly abundant in breastmilk, their functions remain elusive. Human milk (HM) miRNAs start to raise considerable interest, but a comprehensive understanding of the repertoire and expression profiles along lactation has not been well characterized.

OBJECTIVES

This study aimed to characterize the longitudinal profile of HM miRNA between the second week and third month postpartum.

METHODS

We used a new sensitive technology to measure HM miRNAs in a cohort of 44 French mothers [mean ± SD age: 31 ± 3.5; BMI (in kg/m2) 21.8 ± 2.3] who delivered at term and provided HM samples at 3 time points (17 ± 3 d, 60 ± 3 d, and 90 ± 3 d) during follow-up visits.

RESULTS

We detected 685 miRNAs, of which 35 showed a high and stable expression along the lactation period analyzed. We also described for the first time a set of 11 miRNAs with a dynamic expression profile. To gain insight into the potential functional relevance of this set of miRNAs, we selected miR-3126 and miR-3184 to treat undifferentiated Caco-2 human intestinal cells and then assessed differentially expressed genes and modulation of related biological pathways.

CONCLUSIONS

Overall, our study provides new insights into HM miRNA composition and, to our knowledge, the first description of its longitudinal dynamics in mothers who delivered at term. Our in vitro results obtained in undifferentiated Caco-2 human intestinal cells transfected with HM miRNAs also provide further support to the hypothesized mother-to-neonate signaling role of HM miRNAs. This trial was registered at clinicaltrials.gov as NCT01894893.

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.

Hematopoietic PBX-interacting protein is a novel regulator of mammary epithelial cell differentiation

Hematopoietic PBX-interacting protein (HPIP, also known as PBXIP1) is an estrogen receptor (ER) interacting protein that regulates estrogen-mediated breast cancer cell proliferation and tumorigenesis. However, its functional significance in the context of mammary gland development is unexplored. Here, we report that HPIP is required for prolactin (PRL)-induced lactogenic differentiation in vitro. Molecular analysis of HPIP expression in mice revealed its induced expression at pregnancy and lactation stages of mammary gland. Moreover, PRL is a lactogenic hormone that controls pregnancy as well as lactation and induces Hpip/Pbxip1 expression in a signal transducer and activator of transcription 5a-dependent manner. Using mammary epithelial and lactogenic-competent cell lines, we further show that HPIP plays a regulatory role in PRL-mediated mammary epithelial cell differentiation, which is measured by acini formation, β-casein synthesis, and lipid droplet formation. Further mechanistic studies using pharmacological inhibitors revealed that HPIP modulates PRL-induced β-casein synthesis via phosphoinositide 3-kinase (PI3K)/protein kinase B (AKT) activation. This study also identified HPIP as a critical regulator of autocrine PRL signaling as treatment with the PRL receptor antagonist Δ1-9-G129R-hPRL restrained HPIP-mediated PRL synthesis, AKT activation, and β-casein synthesis in cultured HC11 cells. Interestingly, we also uncovered that microRNA-148a (miR-148a) antagonizes HPIP-mediated mammary epithelial cell differentiation. Together, our study identified HPIP as a critical regulator of PRL signaling and revealed a novel molecular circuitry involving PRL, HPIP, PI3K/AKT, and miR-148a that controls mammary epithelial cell differentiation in vitro.

In silico identification of variations in microRNAs with a potential impact on dairy traits using whole ruminant genome SNP datasets

MicroRNAs are small noncoding RNAs that have important roles in the lactation process and milk biosynthesis. Some polymorphisms have been studied in various livestock species from the perspective of pathology or production traits. To target variants that could be the causal variants of dairy traits, genetic variants of microRNAs expressed in the mammary gland or present in milk and localized in dairy quantitative trait loci (QTLs) were investigated in bovine, caprine, and ovine species. In this study, a total of 59,124 (out of 28 millions), 13,427 (out of 87 millions), and 4761 (out of 38 millions) genetic variants in microRNAs expressed in the mammary gland or present in milk were identified in bovine, caprine, and ovine species, respectively. A total of 4679 of these detected bovine genetic variants are located in dairy QTLs. In caprine species, 127 genetic variants are localized in dairy QTLs. In ovine species, no genetic variant was identified in dairy QTLs. This study leads to the detection of microRNA genetic variants of interest in the context of dairy production, taking advantage of whole genome data to identify microRNA genetic variants expressed in the mammary gland and localized in dairy QTLs.

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.

The miR-200 family in normal mammary gland development

The miR-200 family of microRNAs plays a significant role in inhibiting mammary tumor growth and progression, and its members are being investigated as therapeutic targets. Additionally, if future studies can prove that miR-200s prevent mammary tumor initiation, the microRNA family could also offer a preventative strategy. Before utilizing miR-200s in a therapeutic setting, understanding how they regulate normal mammary development is necessary. No studies investigating the role of miR-200s in embryonic ductal development could be found, and only two studies examined the impact of miR-200s on pubertal ductal morphogenesis. These studies showed that miR-200s are expressed at low levels in virgin mammary glands, and elevated expression of miR-200s have the potential to impair ductal morphogenesis. In contrast to virgin mammary glands, miR-200s are expressed at high levels in mammary glands during late pregnancy and lactation. miR-200s are also found in the milk of several mammalian species, including humans. However, the relevance of miR-200s in milk remains unclear. The increase in miR-200 expression in late pregnancy and lactation suggests a role for miR-200s in the development of alveoli and/or regulating milk production. Therefore, studies investigating the consequence of miR-200 overexpression or knockdown are needed to identify the function of miR-200s in alveolar development and lactation.

Lifetime Impact of Cow's Milk on Overactivation of mTORC1: From Fetal to Childhood Overgrowth, Acne, Diabetes, Cancers, and Neurodegeneration

The consumption of cow's milk is a part of the basic nutritional habits of Western industrialized countries. Recent epidemiological studies associate the intake of cow's milk with an increased risk of diseases, which are associated with overactivated mechanistic target of rapamycin complex 1 (mTORC1) signaling. This review presents current epidemiological and translational evidence linking milk consumption to the regulation of mTORC1, the master-switch for eukaryotic cell growth. Epidemiological studies confirm a correlation between cow's milk consumption and birthweight, body mass index, onset of menarche, linear growth during childhood, acne vulgaris, type 2 diabetes mellitus, prostate cancer, breast cancer, hepatocellular carcinoma, diffuse large B-cell lymphoma, neurodegenerative diseases, and all-cause mortality. Thus, long-term persistent consumption of cow's milk increases the risk of mTORC1-driven diseases of civilization. Milk is a highly conserved, lactation genome-controlled signaling system that functions as a maternal-neonatal relay for optimized species-specific activation of mTORC1, the nexus for regulation of eukaryotic cell growth, and control of autophagy. A deeper understanding of milk´s impact on mTORC1 signaling is of critical importance for the prevention of common diseases of civilization.

Impacts of Epigenetic Processes on the Health and Productivity of Livestock

The dynamic changes in the epigenome resulting from the intricate interactions of genetic and environmental factors play crucial roles in individual growth and development. Numerous studies in plants, rodents, and humans have provided evidence of the regulatory roles of epigenetic processes in health and disease. There is increasing pressure to increase livestock production in light of increasing food needs of an expanding human population and environment challenges, but there is limited related epigenetic data on livestock to complement genomic information and support advances in improvement breeding and health management. This review examines the recent discoveries on epigenetic processes due to DNA methylation, histone modification, and chromatin remodeling and their impacts on health and production traits in farm animals, including bovine, swine, sheep, goat, and poultry species. Most of the reports focused on epigenome profiling at the genome-wide or specific genic regions in response to developmental processes, environmental stressors, nutrition, and disease pathogens. The bulk of available data mainly characterized the epigenetic markers in tissues/organs or in relation to traits and detection of epigenetic regulatory mechanisms underlying livestock phenotype diversity. However, available data is inadequate to support gainful exploitation of epigenetic processes for improved animal health and productivity management. Increased research effort, which is vital to elucidate how epigenetic mechanisms affect the health and productivity of livestock, is currently limited due to several factors including lack of adequate analytical tools. In this review, we (1) summarize available evidence of the impacts of epigenetic processes on livestock production and health traits, (2) discuss the application of epigenetics data in livestock production, and (3) present gaps in livestock epigenetics research. Knowledge of the epigenetic factors influencing livestock health and productivity is vital for the management and improvement of livestock productivity.

Exploration of microRNA profiles in human colostrum

Background

Colostrum is well known to have excellent nutritional value for newborns. The aim of this study was to investigate the dynamic expression pattern of microRNA in human colostrum and mature milk. Furthermore, we identified the specific microRNA in human colostrum and analyzed the regulatory function of human colostrum.

Methods

We collected breast milk samples from 18 lactating volunteers. The expression of microRNA in breast milk was detected by microarray analysis. The expression differences were characterized by log2FC (|log2fold change| >1.58) and associated P values (P<0.05). Furthermore, the prediction of microRNA targets, bioinformatics analysis and network generation were carried out using network database.

Results

Our results showed that during the human lactation process, the composition of microRNAs in human milk changes dynamically. Compared to the microRNA expression profile in human mature milk, the expression levels of 49 microRNAs were significantly different and 67 microRNAs were specifically expressed in human colostrum. Based on the results of Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis, the predicted target mRNAs of the identified colostrum-specific microRNAs were involved in the regulation of distinct biological processes, such as signal transduction, positive regulation of GTPase activity, and protein phosphorylation. Moreover, the predicted mRNA targets were from large spectrums of signaling pathways, such as the MAPK, Ras, Hippo, Wnt, and mTOR signaling pathways, as well as the longevity regulating pathway.

Conclusions

Our study illuminates the landscape of microRNA expressions in human colostrum and mature milk, and emphasizes the value of microRNAs as nutritional additives in milk-related commercial products.

Mammary tissue transcriptomic analysis for construction of integrated regulatory networks involved in lactogenesis of Ovis aries

The mammary gland experiences vast changes between the onset of lactation and pregnancy. This remodeling involves different functions such as lactation that is controlled by innumerable regulators and various gene networks which are still not completely understood. MicroRNAs (miRNAs) are one of the important non-coding gene regulators which control an extensive range of biological processes. Thus, exploring miRNAs functions is important for solving gene regulation complexity. The main purpose in the present study is to identify the various gene regulative integrated networks involved in lactation progress in mammary gland. We analyzed ovine mammary tissue data sets which included expression profiles of mRNA (genes) and miRNAs related to six ewes in different days of lactation and nutritional treatments. We combined two different types of information: the network that is module inference by mRNAs (RNA-seq data), miRNAs and transcription factors (TFs) expression matrix and prediction of targets via computational methods. To discover the miRNAs regulatory function, 134 modules were predicted by using gene expression data and 14 TFs and 20 miRNAs were allocated to these predicted modules. By applying this integrated computation-based method, 38 miRNA-modules and 35 TF-module interactions were identified from ovine mammary tissue data during lactogenesis. A lot of these modules were involved in lipid and protein metabolism, as well as steroids and vitamin biosynthesis, which would play key roles in mammary tissue and lactation development. These results present new information about the regulatory procedures at the miRNAs and TF levels throughout lactation.

MicroRNA Milk Exosomes: From Cellular Regulator to Genomic Marker

Recent advances in ruminants' milk-derived exosomes (EXO) have indicated a role of microRNAs (miRNAs) in cell-to-cell communication in dairy ruminants. The miRNAs EXO retain peculiar mechanisms of uptake from recipient cells, which enables the selective delivery of cargos, with a specific regulation of target genes. Although many studies have been published on the miRNAs contained in milk, less information is available on the role of miRNAs EXO, which are considered stable over time and resistant to digestion and milk processing. Several miRNAs EXO have been implicated in the cellular signaling pathway, as in the regulation of immune response. Moreover, they exert epigenetic control, as extenuating the expression of DNA methyltransferase 1. However, the study of miRNAs EXO is still challenging due to the difficulty of isolating EXO. In fact, there are not agreed protocols, and different methods, often time-consuming, are used, making it difficult to routinely process a large number of samples. The regulation of cell functions in mammary glands by miRNAs EXO, and their applications as genomic markers in livestock, is presented.

Leukocytosis and Expression of Bovine Leukemia Virus microRNAs in Cattle

Bovine Leukemia Virus (BLV) is an established model for studying retroviral infections, in particular the infection by the human T-cell leukemia type 1 (HTLV-1) virus. Here, we quantified gene expression of several BLV-related genes: effector protein of T and NK-killer cells NK-lysin (Nklys), reverse BLV transcriptase pol, BLV receptor (blvr), and also key enzymes of the microRNA maturation, Dicer (dc1) and Argonaut (ago2). The differences in the expression of the above genes were compared between five groups: (1) BLV infected cows with high and (2) low lymphocyte count, (3) with and (4) without BLV microRNA expressions, and (5) cows without BLV infections (control group). As compared to control, infected cows with high lymphocyte count and BLV microRNA expression had significantly decreased Nklys gene expression and increased dc1 and ago2 gene expressions. Few infected animals without pol gene expression nevertheless transcribed BLV microRNA, while others with pol gene expression didn't transcribe BLV microRNA. Notably, Pol expression significantly (P < 0.05) correlated with dc1 expression. For infected animals, there were no direct correlations between the number of leukocytes and pol, Nklys, and BLV microRNA gene expressions. Blvr gene expression is typical for juvenile lymphocytes and decreases during terminal differentiation. Our data suggest that BLV infects primarily juvenile lymphocytes, which further divide into two groups. One group expresses BLV DNA and another one expressed BLV microRNA that decreases host immune response against cells, expressing BLV proteins. It is suspected that regulatory microRNAs play a significant role in the bovine leukemia infections, yet the precise mechanisms and targets of the microRNAs remain poorly defined. Vaccines that are currently in use have a low response rate. Understanding of microRNA regulatory mechanisms and targets would allow to develop more effective vaccines for retroviral infections.

Identification and comparison of exosomal microRNAs in the milk and colostrum of two different cow breeds

Bovine milk and colostrum provide essential nutrients and immunologically active factors that are beneficial to a newborn calf. Milk-and-colostrum-derived exosomes are known as the most important for cellular communication. Exosomes also contain non-coding RNA, such as microRNA. However, there is limited information about exosomal miRNA derived from the milk and colostrum of Holstein and DAK cattle. This study aimed to identify and characterize the exosomal microRNA in the milk and colostrum of Holstein and Doğu Anadolu Kirmizisi (DAK) cows. For this purpose, total RNA isolation was carried out on the milk and colostrum samples that were collected from the Holstein and DAK cattle breeds. The RNA samples were subjected to RNA sequencing and the microRNAs were determined. Lastly, gene ontology analysis was performed for target genes. A total of 795 miRNAs that were expressed differently were identified. A total of 545 of these were known miRNAs and 260 were found to be novel miRNAs. In the functional enrichment analysis, the miRNAs expressed in Holstein milk were mostly associated with milk synthesis, and those in colostrum were mostly involved in the immunity pathways. It was also observed that the miRNAs expressed in DAK milk regulated milk fat and protein metabolism, and there were miRNAs that regulated immune pathways in the colostrum. In addition to this, many novel miRNAs were defined in DAK colostrum. When the target genes of exosomal miRNA in Holstein and DAK milk and colostrum were compared, it was suggested that the DAK breed had genes that were mostly associated with the immune system. As a result, the data obtained from this study will provide beneficial contributions to potential miRNA biomarker studies for milk yield and mastitis.

Prolactin-Responsive Circular RNA circHIPK3 Promotes Proliferation of Mammary Epithelial Cells from Dairy Cow

The highly expressed circHIPK3 is a circular RNA that has been previously reported to regulate the growth of human cells. In this study, we found an increased expression of circHIPK3 in bovine mammary epithelial cells treated with prolactin (PRL) in high-throughput sequencing data. Thus, we further investigated the effect of circHIPK3 on the proliferation and differentiation of mammary epithelial cells. We used qRT-PCR/Cell Counting Kit-8 (CCK-8) and a Western blotting analysis to evaluate the effects on cell proliferation. We found that circHIPK3 promotes the proliferation of mammary epithelial cells. The STAT5 signaling pathway was previously associated with the prolactin response and when the STAT5 was suppressed, the expression of circHIPK3 decreased. The results suggest that the response to prolactin and the associated STAT5 signaling pathway affect the expression of circHIPK3, which subsequently affects the proliferation of mammary epithelial cells in dairy cows.

Assessing selection signatures within and between selected lines of dual-purpose black and white and German Holstein cattle

The aim of this study was to detect selection signatures considering cows from the German Holstein (GH) and the local dual-purpose black and white (DSN) population, as well as from generated sub-populations. The 4654 GH and 261 DSN cows were genotyped with the BovineSNP50 Genotyping BeadChip. The geographical herd location was used as an environmental descriptor to create the East-DSN and West-DSN sub-populations. In addition, two further sub-populations of GH cows were generated, using the extreme values for solutions of residual effects of cows for the claw disorder dermatitis digitalis. These groups represented the most susceptible and most resistant cows. We used cross-population extended haplotype homozygosity methodology (XP-EHH) to identify the most recent selection signatures. Furthermore, we calculated Wright's fixation index (F_ST ). Chromosomal segments for the top 0.1 percentile of negative or positive XP-EHH scores were studied in detail. For gene annotations, we used the Ensembl database and we considered a window of 250 kbp downstream and upstream of each core SNP corresponding to peaks of XP-EHH. In addition, functional interactions among potential candidate genes were inferred via gene network analyses. The most outstanding XP-EHH score was on chromosome 12 (at 77.34 Mb) for DSN and on chromosome 20 (at 36.29-38.42 Mb) for GH. Selection signature locations harbored QTL for several economically important milk and meat quality traits, reflecting the different breeding goals for GH and DSN. The average F_ST value between GH and DSN was quite low (0.068), indicating shared founders. For group stratifications according to cow health, several identified potential candidate genes influence disease resistance, especially to dermatitis digitalis.

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.

MicroRNA profiles of dry secretions through the first three weeks of the dry period from Holstein cows

In dairy cows, the period from the end of lactation through the dry period and into the transition period, requires vast physiological and immunological changes critical to mammary health. The dry period is important to the success of the next lactation and intramammary infections during the dry period will adversely alter mammary function, health and milk production for the subsequent lactation. MicroRNAs (miRNAs) are small non-coding RNAs that can post transcriptionally regulate gene expression. We sought to characterize the miRNA profile in dry secretions from the last day of lactation to 3, 10, and 21 days post dry-off. We identified 816 known and 80 novel miRNAs. We found 46 miRNAs whose expression significantly changed (q-value < 0.05) over the first three weeks of dry-off. Additionally, we examined the slopes of random regression models of log transformed normalized counts and cross analyzed the 46 significantly upregulated and downregulated miRNAs. These miRNAs were found to be associated with important components of pregnancy, lactation, as well as inflammation and disease. Detailing the miRNA profile of dry secretions through the dry-off period provides insight into the biology at work, possible means of regulation, components of resistance and/or susceptibility, and outlets for targeted therapy development.

Milk exosomal miRNAs: potential drivers of AMPK-to-mTORC1 switching in β-cell de-differentiation of type 2 diabetes mellitus

Type 2 diabetes mellitus (T2DM) steadily increases in prevalence since the 1950's, the period of widespread distribution of refrigerated pasteurized cow's milk. Whereas breastfeeding protects against the development of T2DM in later life, accumulating epidemiological evidence underlines the role of cow's milk consumption in T2DM. Recent studies in rodent models demonstrate that during the breastfeeding period pancreatic β-cells are metabolically immature and preferentially proliferate by activation of mechanistic target of rapamycin complex 1 (mTORC1) and suppression of AMP-activated protein kinase (AMPK). Weaning determines a metabolic switch of β-cells from a proliferating, immature phenotype with low insulin secretion to a differentiated mature phenotype with glucose-stimulated insulin secretion, less proliferation, reduced mTORC1- but increased AMPK activity. Translational evidence presented in this perspective implies for the first time that termination of milk miRNA transfer is the driver of this metabolic switch. miRNA-148a is a key inhibitor of AMPK and phosphatase and tensin homolog, crucial suppressors of mTORC1. β-Cells of diabetic patients return to the postnatal phenotype with high mTORC1 and low AMPK activity, explained by continuous transfer of bovine milk miRNAs to the human milk consumer. Bovine milk miRNA-148a apparently promotes β-cell de-differentiation to the immature mTORC1-high/AMPK-low phenotype with functional impairments in insulin secretion, increased mTORC1-driven endoplasmic reticulum stress, reduced autophagy and early β-cell apoptosis. In contrast to pasteurized cow's milk, milk's miRNAs are inactivated by bacterial fermentation, boiling and ultra-heat treatment and are missing in current infant formula. Persistent milk miRNA signaling adds a new perspective to the pathogenesis of T2DM and explains the protective role of breastfeeding but the diabetogenic effect of continued milk miRNA signaling by persistent consumption of pasteurized cow's milk.

Bovine milk transcriptome analysis reveals microRNAs and RNU2 involved in mastitis

Mastitis is a common inflammatory infectious disease in dairy cows. To understand the microRNA (miRNA) expression profile changes during bovine mastitis, we undertook a genome-wide miRNA study of normal milk and milk that tested positive on the California mastitis test for bovine mastitis (CMT+). Twenty-five miRNAs were differentially expressed (23 miRNAs upregulated and two downregulated) during bovine mastitis relative to their expression in normal milk. Upregulated mature miR-1246 probably derived from a U2 small nuclear RNA rather than an miR-1246 precursor. The significantly upregulated miRNA precursors and RNU2 were significantly enriched on bovine chromosome 19, which is homologous to human chromosome 17. A gene ontology analysis of the putative mRNA targets of the significantly upregulated miRNAs showed that these miRNAs were involved in binding target mRNA transcripts and regulating target gene expression, and a Kyoto Encyclopedia of Genes and Genomes pathway analysis showed that the upregulated miRNAs were predominantly related to cancer and immune system pathways. Three novel miRNAs were associated with bovine mastitis and were relatively highly expressed in milk. We confirmed that one of the novel mastitis-related miRNAs was significantly upregulated using a digital PCR system. The differentially expressed miRNAs were involved in human cancers, infections, and immune-related diseases. The genome-wide analysis of miRNA profiles in this study provides insight into bovine mastitis and inflammatory diseases. DATABASES: The miRNAseq generated for this study can be found in the Sequence Read Archive (SRA) under BioProject Number PRJNA421075 and SRA Study Number SRP126134 (https://www.ncbi.nlm.nih.gov/bioproject/PRJNA421075).

Breast Milk Supply of MicroRNA Associated with Leptin and Adiponectin Is Affected by Maternal Overweight/Obesity and Influences Infancy BMI

Breast milk constitutes a dietary source of leptin, adiponectin and microRNAs (miRNAs) for newborns. Expression of miRNAs previously associated with maternal obesity, leptin or adiponectin function were assessed and their impact on infant weight analyzed. Milk samples were collected (at month 1, 2, and 3) from a cohort of 59 healthy lactating mothers (38 normal-weight and 21 overweight/obese (BMI ≥ 25)), and infant growth was followed up to 2 years of age. Thirteen miRNAs, leptin and adiponectin were determined in milk. Leptin, adiponectin and miRNA showed a decrease over time of lactation in normal-weight mothers that was altered in overweight/obesity. Furthermore, negative correlations were observed in normal-weight mothers between the expression of miRNAs in milk and the concentration of leptin or adiponectin, but were absent in overweight/obesity. Moreover, miRNAs negatively correlated with infant BMI only in normal-weight mothers (miR-103, miR-17, miR-181a, miR-222, miR-let7c and miR-146b). Interestingly, target genes of milk miRNAs differently regulated in overweight/obesity could be related to neurodevelopmental processes. In conclusion, a set of miRNAs present in breast milk, in close conjunction with leptin and adiponectin, are natural bioactive compounds with the potential to modulate infant growth and brain development, an interplay that is disturbed in the case of maternal overweight/obesity.

Associations between body condition score at parturition and microRNA profile in colostrum of dairy cows as evaluated by paired mapping programs

MicroRNA (miRNA) are abundant in milk, and likely have regulatory activity involving lactation and immunity. The objective of this study was to determine the miRNA profile in colostrum of overconditioned cows compared with cows of more moderate body condition score (BCS) at calving. Multiparous cows with either high (≥4.0 on a scale of 1 to 5; n = 7) or moderate BCS (2.75 to 3.50; n = 9) in the week before parturition were selected from a commercial dairy herd. Blood and colostrum were sampled within 24 h after calving. Blood serum was analyzed for free fatty acid (FFA) concentration. MicroRNA was isolated from colostrum samples after removing milk fat and cells. MicroRNA were sequenced, and reads were mapped to the bovine genome and to the existing database of miRNA at miRBase.org. Two programs, Oasis 2.0 and miRDeep2, were employed in parallel for read alignment, and analysis of miRNA count data was performed using DESeq2. Identification of differentially expressed miRNA from DESeq2 was not affected by the differences in miRNA detected by the 2 mapping programs. Most abundant miRNA included miR-30a, miR-148a, miR-181a, let-7f, miR-26a, miR-21, miR-22, and miR-92a. Large-scale shifts in miRNA profile were not observed; however, colostrum of cows with high BCS contained less miR-486, which has been linked with altered glucose metabolism. Colostrum from cows with elevated serum FFA contained less miR-885, which may be connected to hepatic function during the transition period. Potential functions of abundant miRNA suggest involvement in development and maintenance of cellular function in the mammary gland, with the additional possibility of influencing neonatal tissue and immune system development.

Decreased expression of microRNA-145 promotes the biological functions of fibroblasts in hypertrophic scar tissues by upregulating the expression of transcription factor SOX-9

The present study aimed to determine the expression of microRNA (miRNA or miR)-145 in hypertrophic scars at the tissue and cellular levels, and to investigate its biological functions and mechanism of action. A total of 36 patients who were diagnosed with hypertrophic scar were included in the present study. Reverse transcription-quantitative polymerase chain reaction was used to determine the expression of miR-145 in tissues and fibroblasts. Primary fibroblasts were transfected with negative control miRNA, miR-145 mimics or inhibitor. A Cell Counting Kit-8 assay was performed to determine the level of proliferation of fibroblasts. Flow cytometry was employed for cell cycles determination and apoptosis in fibroblasts. A Matrigel assay was used to evaluate the invasion ability of fibroblasts. Western blotting was used to determine the expression of the transcription factor SOX-9 (SOX-9) protein in fibroblasts. Rescue experiments were performed to examine the effect of SOX-9 on the regulation of fibroblasts by miR-145. The dual luciferase reporter assay was performed to identify the direct interaction between SOX-9 and miR-145. The expression of miR-145 was reduced in hypertrophic tissues and fibroblasts. Overexpression of miR-145 inhibited the proliferation, G1/S phase transition and invasion of fibroblasts, and promoted the apoptosis of fibroblasts. In addition, overexpression of miR-145 inhibited SOX-9 protein expression. By contrast, the expression of SOX-9 reversed the effects of miR-145 on the proliferation, cell cycle, apoptosis and invasion of fibroblasts. The miR-145 seed region was able to bind with the 3′-untranslated region of the SOX-9 mRNA to regulate its expression. The present study demonstrated that miR-145 expression is reduced in hypertrophic scar tissues and negatively associated with SOX-9 expression. In addition, miR-145 inhibits the proliferation, cell cycle and invasion, and promotes the apoptosis of fibroblasts by down-regulating the expression of SOX-9. The current study provides a potential target for the clinical diagnosis and treatment of hypertrophic scars.

Invited review: Advances and challenges in application of feedomics to improve dairy cow production and health

Dairy cattle science has evolved greatly over the past century, contributing significantly to the improvement in milk production achieved today. However, a new approach is needed to meet the increasing demand for milk production and address the increased concerns about animal health and welfare. It is now easy to collect and access large and complex data sets consisting of molecular, physiological, and metabolic data as well as animal-level data (such as behavior). This provides new opportunities to better understand the mechanisms regulating cow performance. The recently proposed concept of feedomics could help achieve this goal by increasing our understanding of interactions between the different components or levels and their impact on animal production. Feedomics is an emerging field that integrates a range of omics technologies (e.g., genomics, epigenomics, transcriptomics, proteomics, metabolomics, metagenomics, and metatranscriptomics) to provide these insights. In this way, we can identify the best strategies to improve overall animal productivity, product quality, welfare, and health. This approach can help research communities elucidate the complex interactions among nutrition, environment, management, animal genetics, metabolism, physiology, and the symbiotic microbiota. In this review, we summarize the outcomes of the most recent research on omics in dairy cows and highlight how an integrated feedomics approach could be applied in the future to improve dairy cow production and health. Specifically, we focus on 2 topics: (1) improving milk yield and milk quality, and (2) understanding metabolic physiology in transition dairy cows, which are 2 important challenges faced by the dairy industry worldwide.

Multi-population GWAS and enrichment analyses reveal novel genomic regions and promising candidate genes underlying bovine milk fatty acid composition

Background

The power of genome-wide association studies (GWAS) is often limited by the sample size available for the analysis. Milk fatty acid (FA) traits are scarcely recorded due to expensive and time-consuming analytical techniques. Combining multi-population datasets can enhance the power of GWAS enabling detection of genomic region explaining medium to low proportions of the genetic variation. GWAS often detect broader genomic regions containing several positional candidate genes making it difficult to untangle the causative candidates. Post-GWAS analyses with data on pathways, ontology and tissue-specific gene expression status might allow prioritization among positional candidate genes.

Results

Multi-population GWAS for 16 FA traits quantified using gas chromatography (GC) in sample populations of the Chinese, Danish and Dutch Holstein with high-density (HD) genotypes detects 56 genomic regions significantly associated to at least one of the studied FAs; some of which have not been previously reported. Pathways and gene ontology (GO) analyses suggest promising candidate genes on the novel regions including OSBPL6 and AGPS on Bos taurus autosome (BTA) 2, PRLH on BTA 3, SLC51B on BTA 10, ABCG5/8 on BTA 11 and ALG5 on BTA 12. Novel genes in previously known regions, such as FABP4 on BTA 14, APOA1/5/7 on BTA 15 and MGST2 on BTA 17, are also linked to important FA metabolic processes.

Conclusion

Integration of multi-population GWAS and enrichment analyses enabled detection of several novel genomic regions, explaining relatively smaller fractions of the genetic variation, and revealed highly likely candidate genes underlying the effects. Detection of such regions and candidate genes will be crucial in understanding the complex genetic control of FA metabolism. The findings can also be used to augment genomic prediction models with regions collectively capturing most of the genetic variation in the milk FA traits.

Electronic supplementary material

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

MicroRNA expression profiles across blood and different tissues in cattle

MicroRNAs (miRNAs) play essential roles in regulating gene expression involved in various biological functions. The knowledge of miRNA expression profiles across different tissues in cattle is still limited. Using the miRNAs data generated from 158 samples in three studies, we characterized the miRNA expression profiles of bovine sera, exosomes and 11 different tissues. Totally 639 miRNAs were identified and 159 miRNAs were expressed in all samples. After relative log expression normalization, four miRNA expression clusters were generated: 1) sera and exosomes; 2) liver; 3) mammary gland; 4) rumen and gut tissues. The top 10 most abundant miRNAs accounted for >55% of total miRNA expression in each tissue. In addition, this study described a detailed pipeline for identification of both tissue and circulating miRNAs, and the shareable datasets can be re-used by researchers to investigate miRNA-related biological questions in cattle. In addition, a web-based repository was developed, which enables researchers to access the distribution range and raw counts number of the miRNA expression data (https://www.cattleomics.com/micrornaome).

Comparative transcriptome analysis of mammary epithelial cells at different stages of lactation reveals wide differences in gene expression and pathways regulating milk synthesis between Jersey and Kashmiri cattle

Jersey and Kashmiri cattle are important dairy breeds that contribute significantly to the total milk production of the Indian northern state of Jammu and Kashmir. The Kashmiri cattle germplasm has been extensively diluted through crossbreeding with Jersey cattle with the goal of enhancing its milk production ability. However, crossbred animals are prone to diseases resulting to unsustainable milk production. This study aimed to provide a comprehensive transcriptome profile of mammary gland epithelial cells at different stages of lactation and to find key differences in genes and pathways regulating milk traits between Jersey and Kashmiri cattle. Mammary epithelial cells (MEC) isolated from milk obtained from six lactating cows (three Jersey and three Kashmiri cattle) on day 15 (D15), D90 and D250 in milk, representing early, mid and late lactation, respectively were used. RNA isolated from MEC was subjected to next-generation RNA sequencing and bioinformatics processing. Casein and whey protein genes were found to be highly expressed throughout the lactation stages in both breeds. Largest differences in differentially expressed genes (DEG) were between D15 vs D90 (1,805 genes) in Kashmiri cattle and, D15 vs D250 (3,392 genes) in Jersey cattle. A total of 1,103, 1,356 and 1,397 genes were differentially expressed between Kashmiri and Jersey cattle on D15, D90 and D250, respectively. Antioxidant genes like RPLPO and RPS28 were highly expressed in Kashmiri cattle. Differentially expressed genes in both Kashmiri and Jersey were enriched for multicellular organismal process, receptor activity, catalytic activity, signal transducer activity, macromolecular complex and developmental process gene ontology terms. Whereas, biological regulation, endopeptidase activity and response to stimulus were enriched in Kashmiri cattle and, reproduction and immune system process were enriched in Jersey cattle. Most of the pathways responsible for regulation of milk production like JAK-STAT, p38 MAPK pathway, PI3 kinase pathway were enriched by DEG in Jersey cattle only. Although Kashmiri has poor milk production efficiency, the present study suggests possible physicochemical and antioxidant properties of Kashmiri cattle milk that needs to be further explored.

Enhanced supply of methionine or arginine alters mechanistic target of rapamycin signaling proteins, messenger RNA, and microRNA abundance in heat-stressed bovine mammary epithelial cells in vitro

Heat stress (HS) causes reductions in milk production, but it is unclear whether this effect is due to reduced number or functional capacity (or both) of mammary cells. Methionine supplementation improves milk protein, whereas Arg is taken up in excess by mammary cells to produce energy and nonessential AA that can be incorporated into milk protein. To evaluate molecular mechanisms by which mammary functional capacity is affected by HS and Met or Arg, mammary alveolar (MAC-T) cells were incubated at thermal-neutral (37°C) or HS (42°C) temperatures. Treatments were optimal AA profiles (control; Lys:Met = 2.9:1.0; Lys:Arg = 2.1:1.0), control plus Met (Lys:Met = 2.5:1.0), or control plus Arg (Lys:Arg = 1.0:1.0). After incubation for 6 h, cells were harvested and RNA and protein were extracted for quantitative real-time PCR and Western blotting. Protein abundance of mechanistic target of rapamycin (MTOR), eukaryotic initiation factor 2a, serine-threonine protein kinase (AKT), 4E binding protein 1 (EIF4EBP1), and phosphorylated EIF4EBP1 was lower during HS. The lower phosphorylated EIF4EBP1 with HS would diminish translation initiation and reduce protein synthesis. Both Met and Arg had no effect on MTOR proteins, but the phosphorylated EIF4EBP1 decreased by AA, especially Arg. Additionally, Met but not Arg decreased the abundance of phosphorylated eukaryotic elongation factor 2, which could be positive for protein synthesis. Although HS upregulated the heat shock protein HSPA1A, the apoptotic gene BAX, and the translation inhibitor EIF4EBP1, the mRNA abundance of PPARG, FASN, ACACA (lipogenesis), and BCL2L1 (antiapoptotic) decreased. Greater supply of Met or Arg reversed most of the effects of HS occurring at the mRNA level and upregulated the abundance of HSPA1A. In addition, compared with the control, supply of Met or Arg upregulated genes related to transcription and translation (MAPK1, MTOR, SREBF1, RPS6KB1, JAK2), insulin signaling (AKT2, IRS1), AA transport (SLC1A5, SLC7A1), and cell proliferation (MKI67). Upregulation of microRNA related to cell growth arrest and apoptosis (miR-34a, miR-92a, miR-99, and miR-184) and oxidative stress (miR-141 and miR-200a) coupled with downregulation of fat synthesis-related microRNA (miR-27ab and miR-221) were detected with HS. Results suggest that HS has a direct negative effect on synthesis of protein and fat, mediated in part by coordinated changes in mRNA, microRNA, and protein abundance of key networks. The positive responses with Met and Arg raise the possibility that supplementation with these AA during HS might have a positive effect on mammary metabolism.

Exosomes of pasteurized milk: potential pathogens of Western diseases

Milk consumption is a hallmark of western diet. According to common believes, milk consumption has beneficial effects for human health. Pasteurization of cow's milk protects thermolabile vitamins and other organic compounds including bioactive and bioavailable exosomes and extracellular vesicles in the range of 40-120 nm, which are pivotal mediators of cell communication via systemic transfer of specific micro-ribonucleic acids, mRNAs and regulatory proteins such as transforming growth factor-β. There is compelling evidence that human and bovine milk exosomes play a crucial role for adequate metabolic and immunological programming of the newborn infant at the beginning of extrauterine life. Milk exosomes assist in executing an anabolic, growth-promoting and immunological program confined to the postnatal period in all mammals. However, epidemiological and translational evidence presented in this review indicates that continuous exposure of humans to exosomes of pasteurized milk may confer a substantial risk for the development of chronic diseases of civilization including obesity, type 2 diabetes mellitus, osteoporosis, common cancers (prostate, breast, liver, B-cells) as well as Parkinson's disease. Exosomes of pasteurized milk may represent new pathogens that should not reach the human food chain.

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.