Characterisation and comparison of lactating mouse and bovine mammary gland miRNomes

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.

CeRNA Network Reveals the Circular RNA Characterization in Goat Ear Fibroblasts Reprogramming into Mammary Epithelial Cells

Circular RNAs (circRNAs) are a type of non-coding RNA that play a crucial role in the development and lactation of mammary glands in mammals. A total of 107 differentially expressed circRNAs (DE circRNAs) were found, of which 52 were up-regulated and 55 were down-regulated. We also found that DE circRNA host genes were mainly involved in GO terms related to the development process of mammary epithelial cells and KEGG pathways were mostly related to mammary epithelial cells, lactation, and gland development. Protein network analysis found that DE circRNAs can competitively bind to miRNAs as key circRNAs by constructing a circRNA-miRNA-mRNA network. CircRNAs competitively bind to miRNAs (miR-10b-3p, miR-671-5p, chi-miR-200c, chi-miR-378-3p, and chi-miR-30e-5p) involved in goat mammary gland development, mammary epithelial cells, and lactation, affecting the expression of core genes (CDH2, MAPK1, ITGB1, CAMSAP2, and MAPKAPK5). Here, we generated CiMECs and systematically explored the differences in the transcription profile for the first time using whole-transcriptome sequencing. We also analyzed the interaction among mRNA, miRNA, and cirRNA and predicted that circRNA plays an important role in the maintenance of mammary epithelial cells.

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.

Dynamic miRNA Landscape Links Mammary Gland Development to the Regulation of Milk Protein Expression in Mice

Simple Summary

Milk synthesis is vital for maintaining the normal growth of newborn animals. Abnormal mammary gland development leads to a decrease in female productivity and the overall productivity of animal husbandry. This study characterized the dynamic miRNA expression profile during the process of mammary gland development, and identified a novel miRNA regulating expression of β-casein—an important milk protein. The results are valuable for studying mammary gland development, increasing milk production, improving the survival rate of pups, and promoting the development of animal husbandry.

Abstract

Mammary gland morphology varies considerably between pregnancy and lactation status, e.g., virgin to pregnant and lactation to weaning. Throughout these critical developmental phases, the mammary glands undergo remodeling to accommodate changes in milk production capacity, which is positively correlated with milk protein expression. The purpose of this study was to investigate the microRNA (miRNA) expression profiles in female ICR mice’s mammary glands at the virgin stage (V), day 16 of pregnancy (P16d), day 12 of lactation (L12d), day 1 of forced weaning (FW1d), and day 3 of forced weaning (FW3d), and to identify the miRNAs regulating milk protein gene expression. During the five stages of testing, 852 known miRNAs and 179 novel miRNAs were identified in the mammary glands. Based on their expression patterns, the identified miRNAs were grouped into 12 clusters. The expression pattern of cluster 1 miRNAs was opposite to that of milk protein genes in mammary glands in all five different stages. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses revealed that the predicted target genes of cluster 1 miRNAs were related to murine mammary gland development and lactation. Furthermore, fluorescence in situ hybridization (FISH) analysis revealed that the novel-mmu-miR424-5p, which belongs to the cluster 1 miRNAs, was expressed in murine mammary epithelial cells. The dual-luciferase reporter assay revealed that an important milk protein gene—β-casein (CSN2)—was regarded as one of the likely targets for the novel-mmu-miR424-5p. This study analyzed the expression patterns of miRNAs in murine mammary glands throughout five critical developmental stages, and discovered a novel miRNA involved in regulating the expression of CSN2. These findings contribute to an enhanced understanding of the developmental biology of mammary glands, providing guidelines for increasing lactation efficiency and milk quality.

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

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

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

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

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.

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.

Bta-miR-34b controls milk fat biosynthesis via the Akt/mTOR signaling pathway by targeting RAI14 in bovine mammary epithelial cells

BACKGROUND

The biosynthesis of milk fat affects both the technological properties and organoleptic quality of milk and dairy products. MicroRNAs (miRNAs) are endogenous small non-coding RNAs that inhibit the expression of their mRNA targets and are involved in downstream signaling pathways that control several biological processes, including milk fat synthesis. miR-34b is a member of the miR-34 miRNA cluster, which is differentially expressed in the mammary gland tissue of dairy cows during lactation and dry periods. Previous studies have indicated miR-34b is a potential candidate gene that plays a decisive role in regulating milk fat synthesis; therefore, it is important to focus on miR-34b and investigate its regulatory effect on the biosynthesis of milk fat in bovine mammary epithelial cells (BMECs).

RESULTS

In this study, elevated miR-34b levels reduced milk fat synthesis, upregulated 1,999 genes, and downregulated 2,009 genes in BMECs. Moreover, Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis of differentially expressed genes suggested that miR-34b may play an inhibitory role in milk fat synthesis via the protein kinase B (Akt)/mammalian target of rapamycin (mTOR) signaling pathway by reducing phosphorylation levels. Notably, the mTOR activator MHY1485 rescued the inhibitory effect of miR-34b. Furthermore, we demonstrated that retinoic acid-induced protein 14 (RAI14) is a target of miR-34b via TargetScan and immunofluorescence assays. RAI14 mRNA and protein levels were significantly decreased by the miR-34b mimic and increased by the miR-34b inhibitor. Moreover, the reduction in RAI14 levels led to the inhibition of the Akt/mTOR signaling pathway.

CONCLUSIONS

Overall, our results identified a miR-34b-RAI14-Akt/mTOR regulatory network, while also providing a theoretical basis for the molecular breeding of dairy cows.

Dynamics of cattle sperm sncRNAs during maturation, from testis to ejaculated sperm

Background

During epididymal transit, spermatozoa go through several functional maturation steps, resulting from interactions with epididymal secretomes specific to each region. In particular, the sperm membrane is under constant remodeling, with sequential attachment and shedding of various molecules provided by the epididymal lumen fluid and epididymosomes, which also deliver sncRNA cargo to sperm. As a result, the payload of sperm sncRNAs changes during the transit from the epididymis caput to the cauda. This work was designed to study the dynamics of cattle sperm sncRNAs from spermatogenesis to final maturation.

Results

Comprehensive catalogues of sperm sncRNAs were obtained from testicular parenchyma, epididymal caput, corpus and cauda, as well as ejaculated semen from three Holstein bulls. The primary cattle sncRNA sperm content is markedly remodeled as sperm mature along the epididymis. Expression of piRNAs, which are abundant in testis parenchyma, decreases dramatically at epididymis. Conversely, sperm progressively acquires miRNAs, rsRNAs, and tsRNAs along epididymis, with regional specificities. For instance, miRNAs and tsRNAs are enriched in epididymis cauda and ejaculated sperm, while rsRNA expression peaks at epididymis corpus. In addition, epididymis corpus contains mainly 20 nt long piRNAs, instead of 30 nt in all other locations. Beyond the bulk differences in abundance of sncRNAs classes, K-means clustering was performed to study their spatiotemporal expression profile, highlighting differences in specific sncRNAs and providing insights into their putative biological role at each maturation stage. For instance, Gene Ontology analyses using miRNA targets highlighted enriched processes such as cell cycle regulation, response to stress and ubiquitination processes in testicular parenchyma, protein metabolism in epididymal sperm, and embryonic morphogenesis in ejaculated sperm.

Conclusions

Our findings confirm that the sperm sncRNAome does not simply reflect a legacy of spermatogenesis. Instead, sperm sncRNA expression shows a remarkable level of plasticity resulting probably from the combination of multiple factors such as loss of the cytoplasmic droplet, interaction with epididymosomes, and more surprisingly, the putative in situ production and/or modification of sncRNAs by sperm. Given the suggested role of sncRNA in epigenetic trans-generational inheritance, our detailed spatiotemporal analysis may pave the way for a study of sperm sncRNAs role in embryo development.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13072-021-00397-5.

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

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

Nutritional Regulation of Mammary Gland Development and Milk Synthesis in Animal Models and Dairy Species

In mammals, milk is essential for the growth, development, and health. Milk quantity and quality are dependent on mammary development, strongly influenced by nutrition. This review provides an overview of the data on nutritional regulations of mammary development and gene expression involved in milk component synthesis. Mammary development is described related to rodents, rabbits, and pigs, common models in mammary biology. Molecular mechanisms of the nutritional regulation of milk synthesis are reported in ruminants regarding the importance of ruminant milk in human health. The effects of dietary quantitative and qualitative alterations are described considering the dietary composition and in regard to the periods of nutritional susceptibly. During lactation, the effects of lipid supplementation and feed restriction or deprivation are discussed regarding gene expression involved in milk biosynthesis, in ruminants. Moreover, nutrigenomic studies underline the role of the mammary structure and the potential influence of microRNAs. Knowledge from three lactating and three dairy livestock species contribute to understanding the variety of phenotypes reported in this review and highlight (1) the importance of critical physiological stages, such as puberty gestation and early lactation and (2) the relative importance of the various nutrients besides the total energetic value and their interaction.

Epigenetics: New Insights into Mammary Gland Biology

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

The micro-RNA content of unsorted cryopreserved bovine sperm and its relation to the fertility of sperm after sex-sorting

BACKGROUND

The use of sex-sorted sperm in cattle assisted reproduction is constantly increasing. However, sperm fertility can substantially differ between unsorted (conventional) and sex-sorted semen batches of the same sire. Sperm microRNAs (miRNA) have been suggested as promising biomarkers of bull fertility the last years. In this study, we hypothesized that the miRNA profile of cryopreserved conventional sperm is related to bull fertility after artificial insemination with X-bearing sperm. For this purpose, we analyzed the miRNA profile of 18 conventional sperm samples obtained from nine high- (HF) and nine low-fertility (LF) bulls that were contemporaneously used to produce conventional and sex-sorted semen batches. The annual 56-day non-return rate for each semen type (NRRconv and NRRss, respectively) was recorded for each bull.

RESULTS

In total, 85 miRNAs were detected. MiR-34b-3p and miR-100-5p were the two most highly expressed miRNAs with their relative abundance reaching 30% in total. MiR-10a-5p and miR-9-5p were differentially expressed in LF and HF samples (false discovery rate < 10%). The expression levels of miR-9-5p, miR-34c, miR-423-5p, miR-449a, miR-5193-5p, miR-1246, miR-2483-5p, miR-92a, miR-21-5p were significantly correlated to NRRss but not to NRRconv. Based on robust regression analysis, miR-34c, miR-7859 and miR-342 showed the highest contribution to the prediction of NRRss.

CONCLUSIONS

A set of miRNAs detected in conventionally produced semen batches were linked to the fertilizing potential of bovine sperm after sex-sorting. These miRNAs should be further evaluated as potential biomarkers of a sire's suitability for the production of sex-sorted sperm.

Different miRNA contents between mammary epithelial cells and milk fat globules: a random or a targeted process?

MicroRNAs (miRNAs) are small noncoding RNAs present in milk-derived extracellular vesicles and milk fat globules (MFG). Nucleic acid content between the lactating mammary tissue (MT) and MFG are quite similar but discrepancies exist in their miRNA content. Our objective was to identify the origin of these discrepancies, and to evaluate the existence of a possible mechanism sorting miRNAs that will or will not be exported from the mammary epithelial cells (MECs) in bovine MFG. miR-125b-5p, miR-126-3p, miR-141-3p, and miR-204-5p, chosen on the basis of their abundance in the MT, were quantified using RT-qPCR in lactating cow MT, MFG, and laser capture-microdissected MECs. Two miRNAs (miR-125b-5p and miR-141-3p) were detected in the MT as well as in MFG and MECs. miR-204-5p was detected only in the MT, suggesting that it is very likely expressed in a cell type other than MECs. miR-126-3p was detected both in the MT and in MECs but not in MFG, suggesting a targeting mechanism for miRNAs in MECs. These results highlights differences in miRNA content between MECs and MFG may be due to a possibly not random mechanism for loading MFG with miRNA cargos that could involve a variable distribution in MECs or a sorting mechanism.

Bta-miR-223 Targeting CBLB Contributes to Resistance to Staphylococcus aureus Mastitis Through the PI3K/AKT/NF-κB Pathway

Bovine mastitis is an inflammatory condition of the mammary gland often caused by (Staphylococcus aureus) S. aureus infection. The aim of this study was to identify mastitis-related miRNAs and their downstream target genes, and therefore elucidate the regulatory mechanisms involved in disease progression and resistance. Three healthy and three mastitic cows were identified on the basis of the somatic cell count and bacterial culture of their milk, and the histological examination of udder tissues. High-throughput RNA sequencing and bioinformatic analyses revealed that 48 differentially expressed miRNAs (DEMs) in the mastitic udder tissues relative to the healthy tissues. Among 48 DEMs, the expression level of bta-miR-223 was the most up-regulated. Overexpression of the bta-miR-223 in Mac-T cells mitigated the inflammatory pathways induced by S. aureus-derived lipoteichoic acid (LTA). The Cbl proto-oncogene B (CBLB) was identified as the target gene of bta-miR-223, and the direct binding of the miRNA to the CBLB promoter was confirmed by dual luciferase reporter assay using wild-type and mutant 3'-UTR constructs. Furthermore, overexpression of CBLB in the LTA-stimulated Mac-T cells significantly upregulated PI3K, AKT, and phosphorylated NF-κB p65, whereas CBLB knockdown had the opposite effect. Consistent with the in vitro findings, the mammary glands of mice infected with 10⁸CFU/100 μL S. aureus showed high levels of CBLB, PI3K, AKT, and p-NF-κB p65 48 h after infection. Taken together, bta-miR-223 is a predominant miRNA involved in mastitis, and bta-miR-223 likely mitigates the inflammatory progression by targeting CBLB and inhibiting the downstream PI3K/AKT/NF-κB pathway.

A comprehensive overview of bull sperm-borne small non-coding RNAs and their diversity across breeds

Background

Mature sperm carry thousands of RNAs, including mRNAs, lncRNAs, tRNAs, rRNAs and sncRNAs, though their functional significance is still a matter of debate. Growing evidence suggests that sperm RNAs, especially sncRNAs, are selectively retained during spermiogenesis or specifically transferred during epididymis maturation, and are thus delivered to the oocyte at fertilization, providing resources for embryo development. However , a deep characterization of the sncRNA content of bull sperm and its expression profile across breeds is currently lacking. To fill this gap, we optimized a guanidinium–Trizol total RNA extraction protocol to prepare high-quality RNA from frozen bull sperm collected from 40 representative bulls from six breeds. Deep sequencing was performed (40 M single 50-bp reads per sample) to establish a comprehensive repertoire of cattle sperm sncRNA.

Results

Our study showed that it comprises mostly piRNAs (26%), rRNA fragments (25%), miRNAs (20%) and tRNA fragments (tsRNA, 14%). We identified 5p-halves as the predominant tsRNA subgroup in bull sperm, originating mostly from Gly and Glu isoacceptors. Our study also increased by ~ 50% the sperm repertoire of known miRNAs and identified 2022 predicted miRNAs. About 20% of sperm miRNAs were located within genomic clusters, expanding the list of known polycistronic pri-miRNA clusters and defining several networks of co-expressed miRNAs. Strikingly, our study highlighted the great diversity of isomiRs, resulting mainly from deletions and non-templated additions (A and U) at the 3p end. Substitutions within miRNA sequence accounted for 40% of isomiRs, with G>A, U>C and C>U substitutions being the most frequent variations. In addition, many sncRNAs were found to be differentially expressed across breeds.

Conclusions

Our study provides a comprehensive overview of cattle sperm sncRNA, and these findings will pave the way for future work on the role of sncRNAs in embryo development and their relevance as biomarkers of semen fertility.

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.

Bta-miR-34b regulates milk fat biosynthesis by targeting mRNA decapping enzyme 1A (DCP1A) in cultured bovine mammary epithelial cells1

Milk fat is a main nutritional component of milk, and it has become one of the important traits of dairy cow breeding. Recently, there is increasing evidence that microRNAs (miRNA) play significant roles in the process of milk fat synthesis in the mammary gland. Primary bovine mammary epithelial cells (BMEC) were harvested from midlactation cows and cultured in DMEM/F-12 medium with 10% fetal bovine serum, 100 units/mL penicillin, 100 µg/mL streptomycin, 5 µg/mL bovine insulin, 1 µg/mL hydrocortisone, and 2 µg/mL bovine prolactin. We found that miR-34b mimic transfection in BMEC reduced the content of intracellular triacylglycerol (TAG) and lipid droplet accumulation via triacylglycerol assay and Oil Red O staining; meanwhile, overexpression of miR-34b inhibited mRNA expression of lipid metabolism-related genes such as peroxisome proliferator-activated receptor gamma (PPARγ), fatty acid synthase (FASN), fatty acid binding protein 4 (FABP4), and CCAAT enhancer binding protein alpha (C/EBPα). Whereas miR-34b inhibitor resulted in completely opposite results. Furthermore, q-PCR and western blot analysis revealed the mRNA and protein expression levels of DCP1A were downregulated in miR-34b mimic transfection group and upregulated in miR-34b inhibitor group. Moreover, luciferase reporter assays verified that DCP1A was the direct target of miR-34b and DCP1A gene silencing in BMEC-inhibited TAG accumulation and suppressed lipid droplet formation. In conclusion, these findings revealed a novel miR-34b-DCP1A axis that has a significant role in regulating milk fat synthesis and suggested that miR-34b may be used to improve the beneficial ingredients in milk.

A Comprehensive Analysis of Runs of Homozygosity of Eleven Cattle Breeds Representing Different Production Types

Simple Summary

Runs of homozygosity (ROH) regions are known to be common in the genomes of cattle and have become a subject of interest of various research in recent years. ROH can be used as a valuable tool to estimate inbreeding, which needs to be controlled in livestock populations. Moreover, analysis of ROH is considered to be an effective method of identifying genome regions that are a subject of selection pressure, which may help in understanding the genetic aspects of production traits under selection. In this study, we analyzed ROH characteristics of 11 cattle breeds, both commercial and native, maintained in Poland. We presented distinct differences in the length, quantity and frequency of ROH between the analyzed breeds as well as in the levels of genomic inbreeding. Higher levels of inbreeding were characteristic for commercial breeds, especially beef breeds. In addition, within ROH islands, we observed a number of genes with a confirmed influence on the level of production traits. The presented results and identified genes can be a basis for further research focused on the identification of genes and markers essential in the determination of the most important production traits in cattle.

Abstract

In the presented research, BovineSNP50 microarrays (Illumina) were applied to determine runs of homozygosity in the genomes of 11 cattle breeds maintained in Poland. These cattle breeds represent three basic utility types: milk, meat and dual purpose. Analysis of runs of homozygosity allowed the evaluation of the level of autozygosity within each breed in order to calculate the genomic inbreeding coefficient (F_ROH), as well as to identify regions of the genome with a high frequency of ROH occurrence, which may reflect traces of directional selectin left in their genomes. Visible differences in the length and distribution of runs of homozygosity in the genomes of the analyzed cattle breeds have been observed. The highest mean number and mean sums of lengths of runs of homozygosity were characteristic for Hereford cattle and intermediate for the Holstein-Friesian Black-and-White variety, Holstein-Friesian Red-and-White variety, Simmental, Limousin, Montbeliarde and Charolais breeds. However, lower values were observed for cattle of conserved breeds. Moreover, the selected livestock differed in the level of inbreeding estimated using the F_ROH coefficient. In regions of the genome with a high frequency of ROH occurrence, which may reflect the impact of directional selection, a number of genes were observed that can be potentially related to the production traits which are under selection pressure for specific production types. The most important detected genes were GHR, MSTN, DGAT1, FABP4, and TRH, with a known influence on the milk and meat traits of the studied cattle breeds.

MiR-16a Regulates Milk Fat Metabolism by Targeting Large Tumor Suppressor Kinase 1 (LATS1) in Bovine Mammary Epithelial Cells

Milk contains a number of beneficial fatty acids including short and medium chain and unsaturated conjugated and nonconjugated fatty acids. In this study, microRNA sequencing of mammary tissue collected in early-, peak-, mid-, and late-lactation periods was performed to determine the miRNA expression profiles. miR-16a was one of the differentially expressed miRNA and was selected for in-depth functional studies pertaining to fatty acid metabolism. The mimic of miR-16a impaired fat metabolism [triacylglycerol (TAG) and cholesterol] while knock-down of miR-16a promoted fat metabolism in vitro in bovine mammary epithelial cells (BMECs). In addition, the in vitro work with BMECs also revealed that miR-16a had a negative effect on the cellular concentration of cis 9-C18:1, total C18:1, C20:1, and C22:1 and long-chain polyunsaturated fatty acids. Therefore, these data suggesting a negative effect on fatty acid metabolism extend the discovery of the key role of miR-16a in mediating adipocyte differentiation. Through a combination of bioinformatics analysis, target gene 3' UTR luciferase reporter assays, and western blotting, we identified large tumor suppressor kinase 1 (LATS1) as a target of miR-16a. Transfection of siRNA-LATS1 into BMECs led to increases in TAG, cholesterol, and cellular fatty acid concentrations, suggesting a positive role of LATS1 in mammary cell fatty acid metabolism. In summary, data suggest that miR-16a regulates biological processes associated with intracellular TAG, cholesterol, and unsaturated fatty acid synthesis through LATS1. These data provide a theoretical and experimental framework for further clarifying the regulation of lipid metabolism in mammary cells of dairy cows.

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

Background

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

Results

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

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

Conclusion

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

Electronic supplementary material

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

Sequence-based GWAS, network and pathway analyses reveal genes co-associated with milk cheese-making properties and milk composition in Montbéliarde cows

BACKGROUND

Milk quality in dairy cattle is routinely assessed via analysis of mid-infrared (MIR) spectra; this approach can also be used to predict the milk's cheese-making properties (CMP) and composition. When this method of high-throughput phenotyping is combined with efficient imputations of whole-genome sequence data from cows' genotyping data, it provides a unique and powerful framework with which to carry out genomic analyses. The goal of this study was to use this approach to identify genes and gene networks associated with milk CMP and composition in the Montbéliarde breed.

RESULTS

Milk cheese yields, coagulation traits, milk pH and contents of proteins, fatty acids, minerals, citrate, and lactose were predicted from MIR spectra. Thirty-six phenotypes from primiparous Montbéliarde cows (1,442,371 test-day records from 189,817 cows) were adjusted for non-genetic effects and averaged per cow. 50 K genotypes, which were available for a subset of 19,586 cows, were imputed at the sequence level using Run6 of the 1000 Bull Genomes Project (comprising 2333 animals). The individual effects of 8.5 million variants were evaluated in a genome-wide association study (GWAS) which led to the detection of 59 QTL regions, most of which had highly significant effects on CMP and milk composition. The results of the GWAS were further subjected to an association weight matrix and the partial correlation and information theory approach and we identified a set of 736 co-associated genes. Among these, the well-known caseins, PAEP and DGAT1, together with dozens of other genes such as SLC37A1, ALPL, MGST1, SEL1L3, GPT, BRI3BP, SCD, GPAT4, FASN, and ANKH, explained from 12 to 30% of the phenotypic variance of CMP traits. We were further able to identify metabolic pathways (e.g., phosphate and phospholipid metabolism and inorganic anion transport) and key regulator genes, such as PPARA, ASXL3, and bta-mir-200c that are functionally linked to milk composition.

CONCLUSIONS

By using an approach that integrated GWAS with network and pathway analyses at the whole-genome sequence level, we propose candidate variants that explain a substantial proportion of the phenotypic variance of CMP traits and could thus be included in genomic evaluation models to improve milk CMP in Montbéliarde cows.

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.

Differential expression of miRNAs in milk exosomes of cows subjected to group relocation

This study examined microRNA (miRNAs) expression and regulatory patterns from milk exosomes of cows experienced group relocation, a common husbandry practice during the lactation period and used as a spontaneous model of stress. Total RNA from milk exosome samples was collected from 3 cows that showed an increased milk cortisol (HRC) after relocation (T2 vs T1) and 3 cows that did not show cortisol increase (LRC). A total of 69 known miRNAs were identified. Thirteen miRNAs were consistently down-regulated at T2 in comparison to T1, such as miR-2904-1, miR-142, miR-2284x and miR-30b-3p. Only two miRNAs, miR-2284z and miR-146a were significantly different between LRC and HRC group. Functional enrichment analyses highlighted that glucocorticoid receptor signaling and neurotrophic factor mediated TRK receptor signaling were among the biological pathways affected by differentially expressed miRNA target genes. Mir-135a-5p and miR-320a were involved in both biological pathways. miR-142-5p, miR-142-3p, miR-30b-5p and miR-320a shared the same target genes belonging to the RAS superfamily of small GTP-binding proteins (RAC1, RAP1A and RASA1), involved in neurotrophin-mediated cell survival. MiR-142, miR-135 and miR-320a in milk exosomes were the most responsive to group relocation of cows.

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

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

A High-Fat Diet Rich in Saturated and Mono-Unsaturated Fatty Acids Induces Disturbance of Thyroid Lipid Profile and Hypothyroxinemia in Male Rats

SCOPE

Increasing evidence has shown that the disturbance of lipid metabolism might make a possible contribution to the pathogenesis of organ dysfunction, including thyroid, yet it is unknown whether excess intake of dietary fat interferes in thyroid lipid profile. We investigate the effects of dietary fat toward the thyroid lipid profile and thyroid function.

METHODS AND RESULTS

Male Sprague-Dawley rats are fed with high-fat diet (HFD) rich in saturated and mono-unsaturated fatty acids or chow diet for 18 weeks. LC-MS analysis of thyroid shows that total free fatty acids (FFAs) content is significantly higher in HFD rats. The concentration of highly saturated triglycerides significantly increases in HFD rats, whereas the polyunsaturated triglyceride significantly decreases, indicating the decrease in unsaturation in the HFD group. Significant increase of lysophosphatidylcholines (LPC) is observed in HFD rats. Thyroid function tests show hypothyroxinemia (total thyroxine [TT₄ ] and free thyroxine [FT₄ ]) in HFD rats, and elevated thyrotropin (TSH) concentration. The HFD rats also show decreased thyroid uptake of iodine.

CONCLUSION

Excess intake of dietary fat induces disturbance of thyroid lipid profile and hypothyroxinemia, indicating thyroid dysfunction. We speculate that it may provide a new prospect in understanding the pathogenesis of hypothyroidism.

Sunflower oil supplementation affects the expression of miR-20a-5p and miR-142-5p in the lactating bovine mammary gland

Oil supplementation in dairy cattle diets is used to modulate milk fat composition, as well as the expression of mammary lipogenic genes, whose regulation remains unclear. MiRNAs are small non-coding RNA considered as crucial regulators of gene expression, offering clues to explain the mechanism underlying gene nutriregulation. The present study was designed to identify miRNAs whose expression in the cow mammary gland is modulated by sunflower oil supplementation. MiRNomes were obtained using RNAseq technology from the mammary gland of lactating cows receiving a low forage diet, supplemented or not with 4% sunflower oil. Among the 272 miRNAs characterized, eight were selected for RT-qPCR validations, showing the significant down-regulation of miR-142-5p and miR-20a-5p by sunflower supplementation. These two miRNAs are predicted to target genes whose expression was reported as differentially expressed by sunflower supplementation. Among their putative targets, ELOVL6 gene involved in lipid metabolism has been studied. However, a first analysis did not show its significant down-regulation, in response to the over-expression of miR-142-5p, of miR-20a-5p, or both, in a bovine mammary epithelial cell line. However, a clearer understanding of the miRNA expression by lipid supplementation would help to decipher the regulation of lactating cow mammary gland in response to nutrition.

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

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

The Use of "Omics" in Lactation Research in Dairy Cows

“Omics” is the application of genomics, transcriptomics, proteomics, and metabolomics in biological research. Over the years, tremendous amounts of biological information has been gathered regarding the changes in gene, mRNA and protein expressions as well as metabolites in different physiological conditions and regulations, which has greatly advanced our understanding of the regulation of many physiological and pathophysiological processes. The aim of this review is to comprehensively describe the advances in our knowledge regarding lactation mainly in dairy cows that were obtained from the “omics” studies. The “omics” technologies have continuously been preferred as the technical tools in lactation research aiming to develop new nutritional, genetic, and management strategies to improve milk production and milk quality in dairy cows.

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

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

Can milk cell or skim milk miRNAs be used as biomarkers for early pregnancy detection in cattle?

The most critical phase of pregnancy is the first three weeks following insemination. During this period about 50% of high yielding lactating cows suffer embryonic loss prior to implantation, which poses a high economic burden on dairy farmers. Early diagnosis of pregnancy in cattle is therefore essential for monitoring breeding outcome and efficient production intervals. Regulated microRNAs (miRNAs) that reach easily accessible body fluids via a ‘liquid biopsy’ could be a new class of pregnancy predicting biomarkers. As milk is obtained regularly twice daily and non-invasively from the animal, it represents an ideal sample material. Our aim was to establish a pregnancy test system based on the discovery of small RNA biomarkers derived from the bovine milk cellular fraction and skim milk of cows. Milk samples were taken on days 4, 12 and 18 of cyclic cows and after artificial insemination, respectively, of the same animals (n = 6). miRNAs were analysed using small RNA sequencing (small RNA Seq). The miRNA profiles of milk cells and skim milk displayed similar profiles despite the presence of immune cell related miRNAs in milk cells. Trends in regulation of miRNAs between the oestrous cycle and pregnancy were found in miR-cluster 25~106b and its paralog cluster 17~92, miR-125 family, miR-200 family, miR-29 family, miR-15a, miR-21, miR-26b, miR-100, miR-140, 193a-5p, miR-221, miR-223, miR-320a, miR-652, miR-2898 and let-7i. A separation of cyclic and pregnant animals was achieved in a principal component analysis. Bta-miRs-29b, -221, -125b and -200b were successfully technically validated using quantitative real-time PCR, however biological validation failed. Therefore we cannot recommend the diagnostic use of these miRNAs in milk as biomarkers for detection of bovine pregnancy for now.

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

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

Lactogenic hormones alter cellular and extracellular microRNA expression in bovine mammary epithelial cell culture

Background

Bovine milk contains not only a variety of nutritional ingredients but also microRNAs (miRNAs) that are thought to be secreted by the bovine mammary epithelial cells (BMECs). The objective of this study was to elucidate the production of milk-related miRNAs in BMECs under the influence of lactogenic hormones.

Results

According to a microarray result of milk exosomal miRNAs prior to cellular analyses, a total of 257 miRNAs were detected in a Holstein cow milk. Of these, 18 major miRNAs of interest in the milk were selected for an expression analysis in BMEC culture that was treated with or without dexamethasone, insulin, and prolactin (DIP) to induce a lactogenic differentiation. Quantitative polymerase chain reaction (qPCR) results showed that the expressions of miR-21–5p (P = 0.005), miR-26a (P = 0.016), and miR-320a (P = 0.011) were lower in the DIP-treated cells than in the untreated cells. In contrast, the expression of miR-339a (P = 0.017) in the cell culture medium were lower in the DIP-treated culture than in the untreated culture. Intriguingly, the miR-148a expression in cell culture medium was elevated by DIP treatment of BMEC culture (P = 0.018). The medium-to-cell expression ratios of miR-103 (P = 0.025), miR-148a (P < 0.001), and miR-223 (P = 0.013) were elevated in the DIP-treated BMECs, suggesting that the lactogenic differentiation-induced secretion of these three miRNAs in BMECs. A bioinformatic analysis showed that the miRNAs down-regulated in the BMECs were associated with the suppression of genes related to transcriptional regulation, protein phosphorylation, and tube development.

Conclusion

The results suggest that the miRNAs changed by lactogenic hormones are associated with milk protein synthesis, and mammary gland development and maturation. The elevated miR-148a level in DIP-treated BMECs may be associated with its increase in milk during the lactation period of cows.

Deep sequencing shows microRNA involvement in bovine mammary gland adaptation to diets supplemented with linseed oil or safflower oil

BACKGROUND

Bovine milk fat composition is responsive to dietary manipulation providing an avenue to modify the content of fatty acids and especially some specific unsaturated fatty acid (USFA) isomers of benefit to human health. MicroRNAs (miRNAs) regulate gene expression but their specific roles in bovine mammary gland lipogenesis are unclear. The objective of this study was to determine the expression pattern of miRNAs following mammary gland adaptation to dietary supplementation with 5 % linseed or safflower oil using next generation RNA-sequencing.

METHODS

Twenty-four Canadian Holstein dairy cows (twelve per treatment) in mid lactation were fed a control diet (total mixed ration of corn:grass silages) for 28 days followed by a treatment period (control diet supplemented with 5 % linseed or safflower oil) of 28 days. Milk samples were collected weekly for fat and individual fatty acid determination. RNA from mammary gland biopsies harvested on day-14 (control period) and on days +7 and +28 (treatment period) from six randomly selected cows per treatment was subjected to small RNA sequencing.

RESULTS

Milk fat percentage decreased significantly (P < 0.001) during treatment with the two diets as compared to the control period. The individual saturated fatty acids C4:0, C6:0, C8:0, C14:0 and C16:0 decreased significantly (P < 0.05) while five USFAs (C14:1, C18:1n11t, C20:3n3, C20:5n3 and CLA:t10c12) increased remarkably (P < 0.05) in response to both treatments. Analysis of 361 million sequence reads generated 321 known bovine miRNAs and 176 novel miRNAs. The expression of fourteen and twenty-two miRNAs was affected (P < 0.05) by linseed and safflower oil treatments, respectively. Seven miRNAs including six up-regulated (bta-miR-199c, miR-199a-3p, miR-98, miR-378, miR-148b and miR-21-5p) and one down-regulated (bta-miR-200a) were found to be regulated (P < 0.05) by both treatments, and thus considered core differentially expressed (DE) miRNAs. The gene targets of core DE miRNAs have functions related to gene expression and general cellular metabolism (P < 0.05) and are enriched in four pathways of lipid metabolism (3-phosphoinositide biosynthesis, 3-phosphoinositide degradation, D-myo-inisitol-5-phosphate metabolism and the superpathway of inositol phosphate compounds).

CONCLUSION

Our results suggest that DE miRNAs in this study might be important regulators of bovine mammary lipogenesis and metabolism. The novel miRNAs identified in this study will further enrich the bovine miRNome repertoire and contribute to understanding mammary gland biology.

Food Deprivation Affects the miRNome in the Lactating Goat Mammary Gland

Background

Nutrition affects milk composition thus influencing its nutritional properties. Nutrition also modifies the expression of mammary genes, whose regulation is not fully understood. MicroRNAs (miRNA) are small non coding RNA which are important post-transcriptional regulators of gene expression by targeting messenger RNAs. Our goal was to characterize miRNA whose expression is regulated by nutrition in the lactating goat mammary gland, which may provide clues to deciphering regulations of the biosynthesis and secretion of milk components.

Methodology/principal findings

Using high-throughput sequencing technology, miRNomes of the lactating mammary gland were established from lactating goats fed ad libitum or deprived of food for 48h affecting milk production and composition. High throughput miRNA sequencing revealed 30 miRNA with an expression potentially modulated by food deprivation; 16 were down-regulated and 14 were up-regulated. Diana-microT predictive tools suggested a potential role for several nutriregulated miRNA in lipid metabolism. Among the putative targets, 19 were previously identified as differently expressed genes (DEG). The functions of these 19 DEG revealed, notably, their involvement in tissue remodelling.

Conclusion/significance

In conclusion, this study offers the first evidence of nutriregulated miRNA in the ruminant mammary gland. Characterization of these 30 miRNA could contribute to a clearer understanding of gene regulation in the mammary gland in response to nutrition.

Epigenetic marks: regulators of livestock phenotypes and conceivable sources of missing variation in livestock improvement programs

Improvement in animal productivity has been achieved over the years through careful breeding and selection programs. Today, variations in the genome are gaining increasing importance in livestock improvement strategies. Genomic information alone, however, explains only a part of the phenotypic variance in traits. It is likely that a portion of the unaccounted variance is embedded in the epigenome. The epigenome encompasses epigenetic marks such as DNA methylation, histone tail modifications, chromatin remodeling, and other molecules that can transmit epigenetic information such as non-coding RNA species. Epigenetic factors respond to external or internal environmental cues such as nutrition, pathogens, and climate, and have the ability to change gene expression leading to emergence of specific phenotypes. Accumulating evidence shows that epigenetic marks influence gene expression and phenotypic outcome in livestock species. This review examines available evidence of the influence of epigenetic marks on livestock (cattle, sheep, goat, and pig) traits and discusses the potential for consideration of epigenetic markers in livestock improvement programs. However, epigenetic research activities on farm animal species are currently limited partly due to lack of recognition, funding and a global network of researchers. Therefore, considerable less attention has been given to epigenetic research in livestock species in comparison to extensive work in humans and model organisms. Elucidating therefore the epigenetic determinants of animal diseases and complex traits may represent one of the principal challenges to use epigenetic markers for further improvement of animal productivity.

Annotation of the goat genome using next generation sequencing of microRNA expressed by the lactating mammary gland: comparison of three approaches

Background

MicroRNAs (miRNA) are small endogenous non-coding RNA involved in the post-transcriptional regulation of specific mRNA targets. The first whole goat genome sequence became available in 2013, with few annotations. Our goal was to establish a list of the miRNA expressed in the mammary gland of lactating goats, thus enabling implementation of the goat miRNA repertoire and considerably enriching annotation of the goat genome.

Results

Here, we performed high throughput RNA sequencing on 10 lactating goat mammary glands. The bioinformatic detection of miRNA was carried out using miRDeep2 software. Three different methods were used to predict, quantify and annotate the sequenced reads. The first was a de novo approach based on the prediction of miRNA from the goat genome only. The second approach used bovine miRNA as an external reference whereas the last one used recently available goat miRNA. The three methods enabled the prediction and annotation of hundreds of miRNA, more than 95% were commonly identified. Using bovine miRNA, 1,178 distinct miRNA were detected, together with the annotation of 88 miRNA for which corresponding precursors could not be retrieved in the goat genome, and which were not detected using the de novo approach or with the use of goat miRNA. Each chromosomal coordinate of the precursors determined here were generated and depicted on a reference localisation map. Forty six goat miRNA clusters were also reported. The study revealed 263 precursors located in goat protein-coding genes, amongst which the location of 43 precursors was conserved between human, mouse and bovine, revealing potential new gene regulations in the goat mammary gland. Using the publicly available cattle QTL database, and cow precursors conserved in the goat and expressed in lactating mammary gland, 114 precursors were located within known QTL regions for milk production and composition.

Conclusions

The results reported here represent the first major identification study on miRNA expressed in the goat mammary gland at peak lactation. The elements generated by this study will now be used as references to decipher the regulation of miRNA expression in the goat mammary gland and to clarify their involvement in the lactation process.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1471-y) contains supplementary material, which is available to authorized users.

Functional Role of the microRNA-200 Family in Breast Morphogenesis and Neoplasia

Branching epithelial morphogenesis is closely linked to epithelial-to-mesenchymal transition (EMT), a process important in normal development and cancer progression. The miR-200 family regulates epithelial morphogenesis and EMT through a negative feedback loop with the ZEB1 and ZEB2 transcription factors. miR-200 inhibits expression of ZEB1/2 mRNA, which in turn can down-regulate the miR-200 family that further results in down-regulation of E-cadherin and induction of a mesenchymal phenotype. Recent studies show that the expression of miR-200 genes is high during late pregnancy and lactation, thereby indicating that these miRs are important for breast epithelial morphogenesis and differentiation. miR-200 genes have been studied intensively in relation to breast cancer progression and metastasis, where it has been shown that miR-200 members are down-regulated in basal-like breast cancer where the EMT phenotype is prominent. There is growing evidence that the miR-200 family is up-regulated in distal breast metastasis indicating that these miRs are important for colonization of metastatic breast cancer cells through induction of mesenchymal to epithelial transition. The dual role of miR-200 in primary and metastatic breast cancer is of interest for future therapeutic interventions, making it important to understand its role and interacting partners in more detail.