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

Transcriptome analysis of mammary epithelial cell between Sewa sheep and East FriEsian sheep from different localities

Mammary epithelial cells, the only milk-producing cell type in the mammary gland, undergo dynamic proliferation and differentiation during pregnancy, culminating in lactation postpartum. The East FriEsian sheep ranks among the world's most prolific dairy breeds, while the Sewa sheep, a unique dual-purpose breed autochthonous to the Qinghai-Tibet Plateau, exhibits significantly lower milk production. Employing tissue culture methods, we successfully established mammary epithelial cell lines from both breeds. Morphological assessment of mammary epithelial cells and immunofluorescence identification of Cytokeratin 7 and Cytokeratin 8 confirmed the epithelial identity of the isolated cells. Subsequent RNA-seq analysis of these in vitro epithelial cell lines revealed 1813 differentially expressed genes (DEGs). Among these, 1108 were significantly up-regulated and 705 were down-regulated in Sewa epithelial sheep cells compared to East FriEsian epithelial cells. KEGG enrichment analysis identified cellular processes, environmental information processing, human diseases, metabolism, and organismal systems as the primary functional categories associated with DEGs. Gene ontology (GO) terms annotation, categorized into molecular function, biological processes, and cellular component, yielded "binding and catalytic activity," "molecular function regulator activity," and "cellular process," "biological regulation," and "regulation of biological process" as the top three terms within each domain, respectively. Clusters of Orthologous Groups of proteins (KOG) classification further revealed that "signal transduction mechanisms" accounted for the largest proportion of DEGs among all KOG categories. Finally, based on these analyses, ATF3 and MPP7 were identified as promising candidate genes for regulating lactation.

Global changes of miRNA expression indicates an increased reprogramming efficiency of induced mammary epithelial cells by repression of miR-222-3p in fibroblasts

Background

Our previous studies have successfully reported the reprogramming of fibroblasts into induced mammary epithelial cells (iMECs). However, the regulatory relationships and functional roles of MicroRNAs (miRNAs) in the progression of fibroblasts achieving the cell fate of iMECs are insufficiently understood.

Methods

First, we performed pre-and post-induction miRNAs sequencing analysis by using high-throughput sequencing. Following that, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment studies were used to determine the primary roles of the significantly distinct miRNAs and targeted genes. Finally, the effect of miR-222-3p on iMECs fate reprogramming in vitro by transfecting.

Results

As a result goat ear fibroblasts (GEFs) reprogramming into iMECs activates a regulatory program, involving 79 differentially expressed miRNAs. Besides, the programming process involved changes in multiple signaling pathways such as adherens junction, TGF-β signaling pathway, GnRH secretion and the prolactin signaling pathway, etc. Furthermore, it was discovered that the expression of miR-222-3p downregulation by miR-222-3p inhibitor significantly increase the reprogramming efficiency and promoted lipid accumulation of iMECs.

Mammary fat globules as a source of mRNA to model alterations in the expression of some milk component genes during lactation in bovines

BACKGROUND

The milk's nutritional value is determined by its constituents, including fat, protein, carbohydrates, and minerals. The mammary gland's ability to produce milk is controlled by a complex network of genes. Thereby, the fat, protein, and lactose synthesis must be boost in milk to increase milk production efficiency. This can be accomplished by fusing genetic advancements with proper management practices. Therefore, this study aimed to investigate the association between the Lipoprotein lipase (LPL), kappa casein CSN3, and Glucose transporter 1 (GLUT1) genes expression levels and such milk components as fat, protein, and lactose in different dairy breeds during different stages of lactation.

METHODS

To achieve such a purpose, 94 milk samples were collected (72 samples from 36 multiparous black-white and red-white Holstein-Friesian (HF) cows and 22 milk samples from 11 Egyptian buffaloes) during the early and peak lactation stages. The milk samples were utilized for milk analysis and genes expressions analyses using non- invasive approach in obtaining milk fat globules (MFGs) as a source of Ribonucleic acid (RNA).

RESULTS

LPL and CSN3 genes expressions levels were found to be significantly higher in Egyptian buffalo than Holstein-Friesian (HF) cows as well as fat and protein percentages. On the other hand, GLUT1 gene expression level was shown to be significantly higher during peak lactation than early lactation. Moreover, lactose % showed a significant difference in peak lactation phase compared to early lactation phase. Also, fat and protein percentages were significantly higher in early lactation period than peak lactation period but lactose% showed the opposite pattern of Egyptian buffalo.

CONCLUSION

Total RNA can be successfully obtained from MFGs. The results suggest that these genes play a role in glucose absorption and lactose synthesis in bovine mammary epithelial cells during lactation. Also, these results provide light on the differential expression of these genes among distinct Holstein-Friesian cow breeds and Egyptian buffalo subspecies throughout various lactation phases.

TWAS revealed significant causal loci for milk production and its composition in Murrah buffaloes

Milk yield is the most complex trait in dairy animals, and mapping all causal variants even with smallest effect sizes has been difficult with the genome-wide association study (GWAS) sample sizes available in geographical regions with small livestock holdings such as Indian sub-continent. However, Transcriptome-wide association studies (TWAS) could serve as an alternate for fine mapping of expression quantitative trait loci (eQTLs). This is a maiden attempt to identify milk production and its composition related genes using TWAS in Murrah buffaloes (Bubalus bubalis). TWAS was conducted on a test (N = 136) set of Murrah buffaloes genotyped through ddRAD sequencing. Their gene expression level was predicted using reference (N = 8) animals having both genotype and mammary epithelial cell (MEC) transcriptome information. Gene expression prediction was performed using Elastic-Net and Dirichlet Process Regression (DPR) model with fivefold cross-validation and without any cross-validation. DPR model without cross-validation predicted 80.92% of the total genes in the test group of Murrah buffaloes which was highest compared to other methods. TWAS in test individuals based on predicted gene expression, identified a significant association of one unique gene for Fat%, and two for SNF% at Bonferroni corrected threshold. The false discovery rates (FDR) corrected P-values of the top ten SNPs identified through GWAS were comparatively higher than TWAS. Gene ontology of TWAS-identified genes was performed to understand the function of these genes, it was revealed that milk production and composition genes were mainly involved in Relaxin, AMPK, and JAK-STAT signaling pathway, along with CCRI, and several key metabolic processes. The present study indicates that TWAS offers a lower false discovery rate and higher significant hits than GWAS for milk production and its composition traits. Hence, it is concluded that TWAS can be effectively used to identify genes and cis-SNPs in a population, which can be used for fabricating a low-density genomic chip for predicting milk production in Murrah buffaloes.

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.

Successful editing and maintenance of lactogenic gene expression in primary bovine mammary epithelial cells

In vitro investigation of bovine lactation processes is limited by a lack of physiologically representative cell models. This deficiency is most evident through the minimal or absent expression of lactation-specific genes in cultured bovine mammary tissues. Primary bovine mammary epithelial cells (pbMECs) extracted from lactating mammary tissue and grown in culture initially express milk protein transcripts at relatively representative levels. However, expression drops dramatically after only three or four passages, which greatly reduces the utility of primary cells to model and further examine lactogenesis. To investigate the effects of alternate alleles in pbMECs including effects on transcription, we have developed methods to deliver CRISPR-Cas9 gene editing reagents to primary mammary cells, resulting in very high editing efficiencies. We have also found that culturing the cells on an imitation basement membrane composed of Matrigel, results in the restoration of a more representative lactogenic gene expression profile and the cells forming three-dimensional structures in vitro. Here, we present data from four pbMEC lines recovered from pregnant cows and detail the expression profile of five key milk synthesis genes in these MECs grown on Matrigel. Additionally, we describe an optimised method for preferentially selecting CRISPR-Cas9-edited cells conferring a knock-out of DGAT1, using fluorescence-activated cell sorting (FACS). The combination of these techniques facilitates the use of pbMECs as a model to investigate the effects of gene introgressions and genetic variation in lactating mammary tissue.

Human Milk Lipids and Small Metabolites: Maternal and Microbial Origins

Although there has been limited application in the field to date, human milk omics research continues to gain traction. Human milk lipidomics and metabolomics research is particularly important, given the significance of milk lipids and metabolites for infant health. For researchers conducting compositional milk analyses, it is important to consider the origins of these compounds. The current review aims to provide a summary of the existing evidence on the sources of human milk lipids and small metabolites. Here, we describe five major sources of milk lipids and metabolites: de novo synthesis from mammary cells, production by the milk microbiota, dietary consumption, release from non-mammary tissue, and production by the gut microbiota. We synthesize the literature to provide evidence and understanding of these pathways in the context of mammary gland biology. We recommend future research focus areas to elucidate milk lipid and small metabolite synthesis and transport pathways. Better understanding of the origins of human milk lipids and metabolites is important to improve translation of milk omics research, particularly regarding the modulation of these important milk components to improve infant health outcomes.

Novel organic selenium source hydroxy-selenomethionine counteracts the blood-milk barrier disruption and inflammatory response of mice under heat stress

Heat stress (HS) in summer has caused huge economic losses to animal husbandry production recently. When mammary gland is exposed to high temperatures, it will cause blood-milk barrier damage. Hydroxy-selenomethionine (HMSeBA) is a new selenium source with better guarantee of animals' production performance under stress, but whether it has protective effect on heat stress-induced blood-milk damage is still unclear. We established mammary epithelial cells and mice heat stress injury models to fill this research gap, and hope to provide theoretical basis for using HMSeBA to alleviate heat stress damage mammary gland. The results showed that (1) Heat stress significantly decreases in vitro transepithelial electrical resistance (TEER) and cell viability (P < 0.01), and significantly decreases clinical score, histological score, and total alveoli area of mice mammary gland tissue (P < 0.01). (2) HMSeBA significantly increases TEER and fluorescein sodium leakage of HS-induced monolayer BMECs (P < 0.01), significantly improves the milk production and total area of alveoli (P < 0.01), and reduces clinical score, histological score, mRNA expression of heat stress-related proteins, and inflammatory cytokines release of heat-stressed mice (P < 0.01). (3) HMSeBA significantly improves tight junction structure damage, and significantly up-regulated the expression of tight junction proteins (ZO-1, claudin 1, and occludin) as well as signal molecules PI3K, AKT, and mTOR (P < 0.01) in heat-stressed mammary tissue. (4) HMSeBA significantly increases glutathione peroxidase (GSH-Px), total antioxidant capacity (T-AOC), and superoxide dismutase release (SOD) (P < 0.01) and significantly reduce malondialdehyde (MDA) expression (P < 0.01) in heat-stressed mammary tissue. In conclusion, this study implemented heat-stressed cell and mice model and showed that HMSeBA significantly regulate antioxidant capacity, inhibited inflammation, and regulate tight junction proteins expression in blood-milk barrier via PI3K/AKT/mTOR signaling pathway, so as to alleviate mammary gland damage and ensure its structure and function integrity.

A flow cytometric method for measuring and isolating mammary epithelial cells from bovine milk

•

Flow cytometry using an antibody against butyrophilin allows quantification of mammary epithelial cells in milk.

•

Sorting butyrophilin-positive or CD45-negative cells isolates epithelial cells from milk.

•

Selection for cytokeratin was not effective in this flow cytometry application.

Sampling frequent time points of mammary signaling pathways is not possible with tissue biopsies. We have validated a flow cytometry and cell sorting procedure for isolating live bovine mammary epithelial cells from somatic cell populations in milk using butyrophilin 1A1 as a marker for mammary epithelial cells and CD45 as a marker for hematopoietic cells. Hoechst 33342 staining and propidium iodide exclusion were used to select for nucleated live cells. Positive selection of butyrophilin (BTN)-expressing cells was performed by fluorescence-activated cell sorting. Quantitative real-time PCR performed on mRNA isolated from these cells showed a 226-fold increase in κ-casein (CSN3) mRNA expression in BTN single-positive cells compared with unsorted cells, whereas CD45 single-positive cells showed a significant decrease. A negative selection strategy for cells not expressing the hematopoietic cell marker CD45 also resulted in a cell population with a 196-fold increase in CSN3 mRNA expression compared with unsorted cells. We found no enrichment of CSN3 mRNA expression after sorting cells using cytokeratin antibodies. The noninvasive assays described here allow for daily or more frequent sampling time points for measurement of mammary epithelial cells during the course of lactation.

SNPs in Mammary Gland Epithelial Cells Unraveling Potential Difference in Milk Production Between Jersey and Kashmiri Cattle Using RNA Sequencing

Deep RNA sequencing experiment was employed to detect putative single nucleotide polymorphisms (SNP) in mammary epithelial cells between two diverse cattle breeds (Jersey and Kashmiri) to understand the variations in the coding regions that reflect differences in milk production traits. The low milk-producing Kashmiri cattle are being replaced by crossbreeding practices with Jersey cattle with the aim of improving milk production. However, crossbred animals are prone to infections and various other diseases resulting in unsustainable milk production. In this study, we tend to identify high-impact SNPs from Jersey and Kashmiri cows (utilizing RNA-Seq data) to delineate key pathways mediating milk production traits in both breeds. A total of 607 (442 SNPs and 169 INDELs) and 684 (464 SNPs and 220 INDELs) high-impact variants were found specific to Jersey and Kashmir cattle, respectively. Based on our results, we conclude that in Jersey cattle, genes with high-impact SNPs were enriched in nucleotide excision repair pathway, ABC transporter, and metabolic pathways like glycerolipid metabolism, pyrimidine metabolism, and amino acid synthesis (glycine, serine, and threonine). Whereas, in Kashmiri cattle, the most enriched pathways include endocytosis pathway, innate immunity pathway, antigen processing pathway, insulin resistance pathway, and signaling pathways like TGF beta and AMPK which could be a possible defense mechanism against mammary gland infections. A varied set of SNPs in both breeds, suggests a clear differentiation at the genomic level; further analysis of high-impact SNPs are required to delineate their effect on these pathways.

Transdifferentiation of goat ear fibroblasts into lactating mammary epithelial cells induced by small molecule compounds

Mammary epithelial cells are the only cells in the mammary glands that are capable of lactation and they are ideal for studying cellular and molecular biology mechanisms during growth, development and lactation of the mammary glands. The limiting factors in most of the currently available mammary epithelial cells are low cell viability, transgenerational efficiency and lactation function that renders them unsuitable for subsequent studies on mammary gland's cellular and lactation mechanisms and utilizing them as bioreactors. Hence, new methods are required to obtain mammary epithelial cells with high transgenerational efficiency and lactation function. In this study, transdifferentiation of goat ear fibroblasts (GEFs) into goat mammary epithelial cells (CiMECs) was induced in only eight days by five small molecule compounds, including 500 μg/mL VPA, 10 μM Tranylcypromine, 10 μM Forskolin, 1 μM TTNPB, 10 μM RepSox. Morphological observation, marker genes comparison, specific antigen expression and comparison of gene expression levels by transcriptome sequencing between the two types of cells that led to the primary deduction that CiMECs have similar biological properties to goat mammary epithelial cells (GMECs) and comparatively more lactation capacity. Therefore, we establish a novel reprogramming route to convert fibroblasts into CiMECs under fully chemically conditions. This study is expected to provide an in vitro platform for understanding cellular mechanisms such as mammary epithelial cells' fate determination and developmental differentiation, and also to find a new way to obtain a large number of functional mammary epithelial cells in vitro.

Single-cell RNA sequencing of freshly isolated bovine milk cells and cultured primary mammary epithelial cells

Bovine mammary function at molecular level is often studied using mammary tissue or primary bovine mammary epithelial cells (pbMECs). However, bulk tissue and primary cells are heterogeneous with respect to cell populations, adding further transcriptional variation in addition to genetic background. Thus, understanding of the variation in gene expression profiles of cell populations and their effect on function are limited. To investigate the mononuclear cell composition in bovine milk, we analyzed a single-cell suspension from a milk sample. Additionally, we harvested cultured pbMECs to characterize gene expression in a homogeneous cell population. Using the Drop-seq technology, we generated single-cell RNA datasets of somatic milk cells and pbMECs. The final datasets after quality control filtering contained 7,119 and 10,549 cells, respectively. The pbMECs formed 14 indefinite clusters displaying intrapopulation heterogeneity, whereas the milk cells formed 14 more distinct clusters. Our datasets constitute a molecular cell atlas that provides a basis for future studies of milk cell composition and gene expression, and could serve as reference datasets for milk cell analysis.

Single Cell RNA Sequencing of Human Milk-Derived Cells Reveals Sub-Populations of Mammary Epithelial Cells with Molecular Signatures of Progenitor and Mature States: a Novel, Non-invasive Framework for Investigating Human Lactation Physiology

Cells in human milk are an untapped source, as potential "liquid breast biopsies", of material for investigating lactation physiology in a non-invasive manner. We used single cell RNA sequencing (scRNA-seq) to identify milk-derived mammary epithelial cells (MECs) and their transcriptional signatures in women with diet-controlled gestational diabetes (GDM) with normal lactation. Methodology is described for coordinating milk collections with single cell capture and library preparation via cryopreservation, in addition to scRNA-seq data processing and analyses of MEC transcriptional signatures. We comprehensively characterized 3740 cells from milk samples from two mothers at two weeks postpartum. Most cells (>90%) were luminal MECs (luMECs) expressing lactalbumin alpha and casein beta and positive for keratin 8 and keratin 18. Few cells were keratin 14⁺ basal MECs and a small immune cell population was present (<10%). Analysis of differential gene expression among clusters identified six potentially distinct luMEC subpopulation signatures, suggesting the potential for subtle functional differences among luMECs, and included one cluster that was positive for both progenitor markers and mature milk transcripts. No expression of pluripotency markers POU class 5 homeobox 1 (POU5F1, encoding OCT4) SRY-box transcription factor 2 (SOX2) or nanog homeobox (NANOG), was observed. These observations were supported by flow cytometric analysis of MECs from mature milk samples from three women with diet-controlled GDM (2-8 mo postpartum), indicating a negligible basal/stem cell population (epithelial cell adhesion molecule (EPCAM)^-/integrin subunit alpha 6 (CD49f)⁺, 0.07%) and a small progenitor population (EPCAM⁺/CD49f⁺, 1.1%). We provide a computational framework for others and future studies, as well as report the first milk-derived cells to be analyzed by scRNA-seq. We discuss the clinical potential and current limitations of using milk-derived cells as material for characterizing human mammary physiology.

Short communication: Molecular markers for epithelial cells across gastrointestinal tissues and fecal RNA in preweaning dairy calves

The objective of this study was to compare the transcription of gene markers for gastrointestinal (GI) epithelial cells, including fatty acid binding protein 2 (FABP2) and cytokeratin 8 (KRT8), and tight junction complex genes (TJP1, CLDN1, CLDN4) in fecal RNA against several GI tract tissue sections in dairy calves. Eight healthy Jersey calves were euthanized at 5 wk of age, and postmortem samples were collected from rumen, duodenum, jejunum, ileum, large intestine, cecum, and feces for total RNA isolation. Tissues and fecal samples were immediately frozen in liquid nitrogen until RNA isolation. A real-time quantitative PCR analysis was performed using a single standard curve composited of equal amounts of all samples, including cDNA from fecal and GI tract tissues. The mRNA expression of the tight junctions TJP1, CLDN1, and CLDN4 was greater in fecal RNA compared with lower GI tract tissues (i.e., duodenum, jejunum, ileum, large intestine, and cecum). Similar to fecal RNA, rumen tissue had greater expression of tight junctions CLDN1 and CLDN4 than lower GI tract tissues. Similarly, rumen tissue had greater expression of TPJ1 than all lower GI tract tissues except duodenum. The expression of TJP1 and CLDN4 was greater in fecal RNA than in rumen tissue; in contrast, CLDN1 mRNA expression was greater in rumen tissue than in the fecal RNA. The expression of FABP2 was greater in duodenum in comparison to all tissue except ileum. The mRNA expression of FABP2 in fecal samples was similar to jejunum and ileum. The expression of KRT8 in fecal samples was similar to duodenum, large intestine, and cecum. The fecal RNA had a greater expression of KRT8 in comparison to jejunum and ileum. The rumen tissue had the lowest mRNA expression of KRT8. The expression levels of FABP2, KRT8, and tight junction genes observed in fecal transcripts suggest that a considerable amount of RNA derived from GI tract epithelial cells can be detected in fecal RNA, which is in agreement with previous data in neonatal dairy calves and other biological models including humans, rodents, and primates. The greater expression of tight junctions in fecal RNA in comparison to sections of the low GI remains to be understood, and due to the importance of tight junctions in GI physiology, further clarification of this effect is warranted. The similarities in mRNA expression of FABP2 and KRT8 between fecal RNA and intestinal sections add up to the accumulating evidence that fecal RNA can be used to investigate molecular alterations in the GI tract of neonatal dairy calves. Further research in this area should include high-throughput transcriptomic analysis via RNA-seq to uncover novel molecular markers for specific sections of the GI tract of neonates.

Effects of micronutrient supplementation on performance and epigenetic status in dairy cows

The postpartum period is crucial in dairy cows and is marked by major physiological and metabolic changes that affect milk production, immune response and fertility. Nutrition remains the most important lever for limiting the negative energy balance and its consequences on general health status in highly selected dairy cows. In order to analyze the effect of a commercial micronutrient on intrinsic parameters, performances and the epigenome of dairy cows, 2 groups of 12 Holstein cows were used: 1 fed a standard diet (mainly composed of corn silage, soybean meal and non-mineral supplement) and the other 1 fed the same diet supplemented with the commercial micronutrient (µ-nutrient supplementation) for 4 weeks before calving and 8 weeks thereafter. Milk production and composition, BW, body condition score (BCS), DM intake (DMI) and health (calving score, metritis and mastitis) were recorded over the study period. Milk samples were collected on D15 and D60 post-calving for analyses of casein, Na+ and K+ contents and metalloprotease activity. Milk leukocytes and milk mammary epithelial cells (mMECs) were purified and counted. The viability of mMECs was assessed, together with their activity, through an analysis of gene expression. At the same time points, peripheral blood mononuclear cells (PBMCs) were purified and counted. Using genomic DNA extracted from PBMCs, mMECs and milk leukocytes, we assessed global DNA methylation (Me-CCGG) to evaluate the epigenetic imprinting associated with the µ-nutrient-supplemented diet. The µ-nutrient supplementation increased BCS and BW without modifying DMI or milk yield and composition. It also improved calving condition, reducing the time interval between calving and first service. Each easily collectable cell type displayed a specific pattern of Me-CCGG with only subtle changes associated with lactation stages in PBMCs. In conclusion, the response to the µ-nutrient supplementation improved the body condition without alteration of global epigenetic status in dairy cows.

Investigation on the suitability of milk-derived primary bovine mammary epithelial cells grown on permeable membrane supports as an in vitro model for lactation

This study aimed to establish an in vitro model for lipid synthesis in primary bovine mammary epithelial cells (pbMECs) extracted from milk and cultured on Transwell permeable supports (TW culture). The suitability of these cells as a functional model for lactation was assessed by measuring κ-casein (CSN3) and diacylglycerol acyl transferase 1 (DGAT1) gene expression, the presence of intracellular lipid droplets, and the concentration of triacylglycerol in the cell lysates. The functionality of the milk-derived pbMECs cultured under lactogenic conditions, with and without oleic acid supplementation, was evaluated by comparing the cells grown on Transwell supports to cells grown on an extracellular matrix (ECM) gel (3D culture) or a plastic surface (2D culture). Furthermore, the functionality of milk-derived cells was compared to pbMECs obtained from bovine mammary tissue. Here, we show that in both tissue and milk-derived pbMECs, 3D culture offered the most suitable in vitro environment and led to increased levels of CSN3 and DGAT1 gene expression along with increased intracellular triacylglycerol content. The TW culture conditions also resulted in increased DGAT1 gene expression compared to the 2D conditions and milk-derived pbMECs cultured on TW inserts showed the highest viability compared to cells grown under 2D or 3D treatments. However, this was not observed for tissue-derived pbMECs, suggesting that TW culture may offer a beneficial environment specifically for milk-derived cells. We suggest that with further optimization of the culture conditions, TW culture may present a suitable model for the study of milk lipid synthesis in pbMECs.

Heat stress induces apoptosis through disruption of dynamic mitochondrial networks in dairy cow mammary epithelial cells

Heat stress-induced reductions in milk yield and the dysfunction of mammary glands are economically important challenges that face the dairy industry, especially during summer. The aim of the present study is to investigate the effects of heat stress on mitochondrial function by using dairy cow mammary epithelial cells (DCMECs) as an in vitro model. Live cell imaging shows that the mitochondria continually change shape through fission and fusion. However, heat stress induces the fragmentation of mitochondria, as well as the decreased of ATP level, membrane potential, and anti-oxidant enzyme activity and the increased of respiratory chain complex I activity. In addition, the cytosolic Ca²⁺ concentration and cytochrome c expression (Cyto-c) were increased after heat stress treatment. Both qRT-PCR and western blot analysis indicate that mitofusin1/2 (Mfn1/2) and optic atrophy protein-1 (Opa-1) are downregulated after heat stress, whereas dynamin-related protein 1 (Drp1) and fission 1 (Fis-1) are upregulated, which explains the observed defect of mitochondrial network dynamics. Accordingly, the present study indicated that heat stress induced the dysfunction of DCMEC through disruption of the normal balance of mitochondrial fission and fusion.

Identification of the complete coding cDNAs and expression analysis of B4GALT1, LALBA, ST3GAL5, ST6GAL1 in the colostrum and milk of the Garganica and Maltese goat breeds to reveal possible implications for oligosaccharide biosynthesis

Background

Milk sialylated oligosaccharides (SOS) play crucial roles in many biological processes. The most abundant free SOS in goat’s milk are 3’sialyllactose (3′-SL), 6’sialyllactose (6′-SL) and disialyllactose (DSL). The production of these molecules is determined genetically by the expression of glycosyltransferases and by the availability of nucleotide sugar substrates, but the precise mechanisms regulating the differential patterns of milk oligosaccharides are not known. We aimed to identify the complete cDNAs of candidate genes implicated in SOS biosynthesis (B4GALT1, LALBA, ST3GAL5, ST6GAL1) and to analyse their expression during lactation in the Garganica and Maltese goat breeds. Moreover, we analysed the colostrum and milk contents of 3′-SL, 6′-SL and disialyllactose (DSL) and the possible correlations between expressed genes and SOS.

Results

We identified the complete coding cDNAs of B4GALT1 (HQ700335.1), ST3GAL5 (KF055858.2), and ST6GAL1 (HQ709167.1), the single nucleotide polymorphism (SNPs) of these genes and 2 splicing variants of the ST6GAL1 cDNA. RT-qPCR analysis showed that LALBA and ST6GAL1 were the genes with the highest and lowest expression in both breeds, respectively. The interaction effects of the breeds and sampling times were associated with higher levels of B4GALT1 and ST3GAL5 gene expression in Garganica than in Maltese goats at kidding. B4GALT1, LALBA, and ST3GAL5 gene expression changed from kidding to 60 and 120 days in Maltese goats, while in Garganica goats, a difference was observed only for the LALBA gene. Breed and lactation effects were also found for SOS contents. Positive correlations of B4GALT1, LALBA, ST3GAL5, and ST6GAL1 with 3′-SL/6′SL and DSL were found.

Conclusions

The genetic effect on the oligosaccharide content of milk was previously highlighted in bovines, and this study is the first to investigate this effect in two goat breeds (Garganica and Maltese) during lactation. The genetic variability of candidate genes involved in SOS biosynthesis highlights their potential role in affecting gene expression and ultimately biological function. The investigation of gene regulatory regions as well as the examination of other sialyltransferase genes will be needed to identify the genetic pattern leading to a higher SOS content in the autochtonous Garganica breed and to protect it using a focused breeding strategy.

Short communication: Comparative gene expression analysis on the enrichment of polymorphonuclear leukocytes and gastrointestinal epithelial cells in fecal RNA from nondiarrheic neonatal dairy calves

Increased understanding of the biology of the gastrointestinal tract (GIT) in neonatal dairy calves during their adaptation to an extrauterine environment will decrease health problems such as diarrhea while increasing feed efficiency and average daily gain in preweaned dairy calves. Within this context, a noninvasive method, based on fecal RNA, to study the GIT in neonatal dairy calves through the isolation of RNA from fecal samples for quantitative reverse-transcription PCR analysis can provide valuable information on GIT biological adaptations during the preweaning period. We aimed to evaluate the potential enrichment of RNA from immune cells or GIT epithelial cells during fecal RNA isolation. Eight neonatal Holstein calves less than 3 wk old (14.9 ± 5.5 d of age at sampling ± standard deviation) and a fecal score of 2.0 ± 0.7 (mean ± standard deviation) were used. During a single sampling, fecal and blood samples were taken simultaneously from each calf before the morning feeding. Fecal samples were immediately frozen in liquid nitrogen until RNA isolation, whereas polymorphonuclear leukocytes (PMN) were isolated from blood samples before RNA isolation. An quantitative reverse-transcription PCR analysis was performed using a single standard curve composited of equal amounts of all samples including cDNA from fecal and PMN. The genes myeloperoxidase (MPO) and L-selectin (SELL) were selected for their specific known function in PMN, whereas keratin 8 (KRT8) and aquaporin 3 (AQP3) have been associated with epithelial enterocytes. Our results showed a contrasting gene expression profile between PMN and fecal RNA; whereas greater mRNA expression of SELL was observed in PMN, a greater KRT8 expression was observed in fecal RNA. The mRNA expression of AQP3 tended to be greater in PMN than fecal RNA. Additionally, MPO was not amplified in fecal RNA. Our findings suggest that under nondiarrheic conditions RNA isolated from stool samples of neonatal dairy calves will have a considerable number of GIT epithelial cells, which confirms the reliability of this method under these conditions. However, further research needs to be done to determine if the same effects are observed during diarrhea or throughout the preweaning period of dairy calves.

NDRG1 negatively regulates proliferation and Milk bio-synthesis of bovine epithelial cells via the mTOR signaling pathway

The expression of N-myc downstream-regulated gene 1 (NDRG1) was significantly correlated with diverse processes such as cell growth and differentiation, lipid synthesis, stress and immune responses. Here we explored the role of NDRG1 expression in bovine mammary tissue and epithelial cells under an inflammatory condition. Results showed that NDRG1 expression was elevated in bovine mammary tissue with mastitis and mammary epithelial cells treated by heat inactivated Escherichia coli and Staphylococcus aureus compared to normal tissue and untreated cells. Overexpression of NDRG1 significantly inhibited cell proliferation and migration, β-casein secretion, gene expressions of inflammatory cytokines tumor necrosis factor-α, interleukin (IL)-6 and IL-8, and activation of mTOR signal pathway of mammary epithelial cells, and vice versa by NDRG1 knockdown. These findings suggest that NDRG1 has immense potential in the regulation of properties in bovine mammary epithelial cells under an inflammatory condition.

The impact of organic vs. inorganic selenium on dairy goat productivity and expression of selected genes in milk somatic cells

The aim of this study was to determine the effect of diet supplemented with selenized yeast (Se-yeast) on milk yield and milk composition of goats and expression of casein and mammary-gland-immune system genes in milk somatic cells (MSC). Twenty-four dairy goats in their second to fourth lactations were divided into control and experimental groups, balanced according to lactation number and breed (Polish White or Fawn Improved). Morning milk and blood samples were collected four times during lactation (on the 21st, 70th, 120th, 180th day after kidding). The control and experimental groups were fed diets with 0.7 mg inorganic Se/goat/day (sodium selenite) or 0.6 mg organic Se/goat/day (selenized yeast), respectively. Milk, fat and protein yields during lactation as well as average somatic cell count, fat, protein and lactose contents in milk were evaluated. Microelements in milk and blood serum and biochemical parameters in blood serum were determined at the beginning and the end of the experiment. The expression levels of the genes encoding αS1-casein (CSN1S1), αS2-casein (CSN1S2), κ-casein (CSN3), interleukin 8 (IL-8), serum amyloid A3 (SAA3), interleukin 1β (IL-1β), bactenecin 7.5 (BAC7.5), bactenecin 5 (BAC5), β2-defensin (GBD2), hepcidin (HAMP), chemokine 4 (CCL4), tumour necrosis factor α (TNFα), toll-like receptor 2 (TLR2), cathelicidin-7 (MAP34) and cathelicidin-6 (MAP28) were determined in MSC. Milk, fat, and protein yields were higher and somatic cell count (SCC expressed as natural logarithm) was lower in the milk of goats fed organic Se. The Se concentration in milk was twice as high in the organic vs. inorganic treatment groups at the end of the experiment, while there were no differences in studied biochemical parameters between groups. The transcript levels of CSN1S2 and BAC7.5 were higher and IL-8 was lower in MSC of Se-yeast treated groups. Such results may indicate better health status of mammary glands of goats treated with organic Se as well as positive impact of selenized yeast on the goat's milk composition. Differences in the IL-1β and IL-8 transcript levels were also noted between the stages of lactation, with the highest expression at the peak of lactation (day 70), highlighting the metabolic burden at this time. We concluded that the Se-yeast supplementation improved the productivity and health status of goats and could have significant economic impact on farmer's income.

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.

Milk yield loss in response to feed restriction is associated with mammary epithelial cell exfoliation in dairy cows

In dairy cows, feed restriction is known to decrease milk yield by reducing the number of mammary epithelial cells (MEC) in the udder through a shift in the MEC proliferation-apoptosis balance, by reducing the metabolic activity of MEC, or both. The exfoliation of MEC from the mammary epithelium into milk is another process that may participate in regulating the number of MEC during feed restriction. The aim of the present study was to clarify the mechanisms that underlie the milk yield loss induced by feed restriction. Nineteen Holstein dairy cows producing 40.0 ± 0.7 kg/d at 77 ± 5 d in milk were divided into a control group (n = 9) and a feed-restricted group (n = 10). Ad libitum dry matter intake (DMI) was recorded during a pre-experimental period of 2 wk. For 29 d (period 1), cows were fed either 100 (control) or 80% (feed-restricted) of their ad libitum DMI measured during the pre-experimental period. Then, all cows were fed ad libitum for 35 d (period 2). Milk production and DMI were recorded daily. Blood and milk samples were collected once during the pre-experimental period; on d 5, 9, and 27 of period 1; and on d 5, 9, and 30 of period 2. Mammary epithelial cells were purified from milk using an immunomagnetic method to determine the rate of MEC exfoliation. Mammary tissue samples were collected by biopsy at the end of each period to analyze the rates of cell proliferation and apoptosis and the expression of genes involved in synthesizing constituents of milk. Feed restriction decreased milk yield by 3 kg/d but had no effect on rates of proliferation and apoptosis in the mammary tissue or on the expression of genes involved in milk synthesis. The daily MEC exfoliation rate was 65% greater in feed-restricted cows than in control cows. These effects in feed-restricted cows were associated with reduced insulin-like growth factor-1 and cortisol plasma concentrations. When all cows returned to ad libitum feeding, no significant difference on milk yield or MEC exfoliation rate was observed between feed-restricted and control cows, but refeeding increased prolactin release during milking. These results show that the exfoliation process may play a role in regulating the number of MEC in the udders of dairy cows during feed restriction without any carryover effect on their milk production.

Small ruminant lentivirus infection influences expression of acute phase proteins and cathelicidin genes in milk somatic cells and peripheral blood leukocytes of dairy goats

The aim of the study was to analyze acute phase protein and cathelicidin gene responses to small ruminant lentivirus (SRLV) infection in goats. In uninfected goats, we found higher Cp and lower Fbγ mRNA levels in blood leucocytes (BL) than in milk somatic cells (MSC), as well as lower SAA, Hp, and CRP and higher Cp and AGP concentrations in blood serum than in milk. In SRLV-infected goats, we found higher Fbγ and MAP28 and lower Cp expression in MSC than in BL, and higher SAA, Hp, Fb, and MAP28 and lower AGP concentrations in milk than in blood serum. Higher SAA and Hp expressions in BL and Hp expression in MSC were found in SRLV-infected goats. In SRLV-infected goats, we observed a higher concentration of SAA in blood serum, while in milk, lower SAA, Cp, and MAP28 and higher MAP34 concentrations were observed. The expression profiles of the studied genes differed between BL/serum and MSC/milk. The elevated SAA concentration in blood serum was accompanied by a decreased concentration of SAA and Cp in the milk of infected goats. No differences in the expression of the other studied genes may mean that the SRLV has the ability to evade the immune system, continuing to replicate. The elevated concentration of SAA in blood serum may promote viral multiplication. This higher concentration of SAA in blood serum and simultaneous reduced concentration of SAA and Cp in milk may be additive indicators of this infection.

An Exploratory Search for Potential Molecular Targets Responsive to the Probiotic Lactobacillus salivarius PS2 in Women With Mastitis: Gene Expression Profiling vs. Interindividual Variability

Probiotics constitute an attractive alternative in the battle against microbial infections. Oral administration of certain strains of lactobacilli isolated from human milk has resulted in an effective reduction of the bacterial load as well as an improvement of the mastitis-associated symptoms. Nevertheless, little is yet known about the potential molecular mechanisms and specific targets implicated in these effects. Transcriptomic profiling has been used to search for disease-associated and therapy-responsive molecules in different disorders and experimental models. We have applied for the first time a gene expression-based molecular approach to explore for potential targets responsive to intervention with a probiotic in: (i) breast milk somatic cells (n = 17) and (ii) blood leukocytes (n = 19). Women with mastitis ingested a new strain of lactobacilli, Lactobacillus salivarius PS2 (3 × capsules per day, each capsule contained ~9.5 log10 CFU) for 21 days. We applied Affymetrix microarrays and Taqman one-step quantitative reverse transcription PCR (RT-qPCR) to analyze and compare gene expression changes between samples pre- and post-treatment. Our results substantiate the involvement of inflammatory and cell-growth related pathways and genes in the breast milk somatic cells following the intake of L. salivarius PS2. Individual analyses of selected genes: (1) supported the upregulation of STC1 and IL19 and the downregulation of PLAUR and IFNGR1 in the somatic cells of the patients as potential targets responsive to the probiotic, (2) detected a lack of a relationship between the gene expression responses in the two types of cells, and (3) evidenced a substantial interindividual variability in the gene expression changes in both types of cells. Our study provides an insight into the essentiality of incorporating the study of tissue-specific interindividual molecular responsivity into future clinical intervention trials to further understand the complexity of human gene expression responses to therapy and the potentiality of selecting appropriate responsive targets.

A two-step immunomagnetic separation of somatic cell subpopulations for a gene expression profile study in bovine milk

The objective of this study was to demonstrate the usefulness of an immunomagnetic method to purify subpopulations of milk somatic cells. The experiment was conducted on milk samples collected from healthy cows (n = 17) and from cows with clinical mastitis (n = 24) due to a Staphylococcus aureus natural infection. A two-step immunomagnetic purification was applied to simultaneously separate three somatic cell subpopulations from the same milk sample. Total RNA was extracted and qPCR was performed to determinate mRNA levels of innate immunity target genes in purified somatic cell subpopulations. Good quality and quantity of RNA allowed the reference gene analysis in each cell subpopulation. An up-regulation of the main genes involved in innate immune defence was detected in separated polymorphonuclear neutrophilic leucocytes-monocytes and lymphocytes of mastitic milk. These results and flow cytometric analysis suggest that the immunomagnetic purification is an efficient method for the isolation of the three populations from milk, allowing the cells to be studied separately.

Milk somatic cells, factors influencing their release, future prospects, and practical utility in dairy animals: An overview

Milk somatic cells (SCs) are a mixture of milk-producing cells and immune cells. These cells are secreted in milk during the normal course of milking and are used as an index for estimating mammary health and milk quality of dairy animals worldwide. Milk SC is influenced by cow productivity, health, parity, lactation stage, and breed of an animal. Any change in environmental conditions, poor management practices, and also stressful conditions significantly increases the amount of SC coming in milk. Better hygiene and proper nutrition help in reducing milk SC. Milk with low SC means better milk products with a longer shelf life. The present review describes the role of SCs (both secretory and immune) in milk, their role in maintaining the integrity of the mammary gland, and factors affecting their release in milk. This information may help to reduce milk somatic cell counts (SCCs) and to establish differential SCC standards.

Relative gene expression of fatty acid synthesis genes at 60 days postpartum in bovine mammary epithelial cells of Surti and Jafarabadi buffaloes

Aim:

Aim of the study was to study the relative gene expression of genes associated with fatty acid synthesis at 60 days postpartum (pp) in bovine mammary epithelial cells (MECs) of Surti and Jafarabadi buffaloes.

Materials and Methods:

A total of 10 healthy Surti and Jafarabadi buffaloes of each breed were selected at random from Livestock Research Station, Navsari and Cattle Breeding Farm, Junagadh, Gujarat, respectively, for this study. Milk sample was collected from each selected buffalo at day 60 pp from these two breeds to study relative gene expression of major milk fat genes using non-invasive approach of obtaining primary bovine MECs (pBMEC) from milk samples.

Results:

In this study overall, the relative expression of the six major milk lipogenic genes butyrophilin subfamily 1 member A1 (BTN1A1), stearoyl-CoA desaturase (SCD), lipoprotein lipase (LPL), glycerol-3-phosphate acyltransferase mitochondrial (GPAM), acetyl-coenzyme A carboxylase alpha (ACACA), and lipin (LPIN) did not show changes in expression patterns at 60^th day of lactation in both Surti and Jafarabadi buffaloes.

Conclusion:

The pBMEC can be successfully recovered from 1500 ml of milk of Surti and Jafarabadi buffaloes using antibody-mediated magnetic bead separation and can be further used for recovering RNA for down step quantification of major milk lipogenic gene expression. The relative expression of the six major milk lipogenic genes BTN1A1, SCD, LPL, GPAM, ACACA, and LPIN did not show changes in expression patterns in both Surti and Jafarabadi buffaloes, suggesting expression levels of lipogenic genes are maintained almost uniform till peak lactation without any significant difference.

RNA-Sequencing for profiling goat milk transcriptome in colostrum and mature milk

BACKGROUND

In this work we aimed at sequencing and assembling the goat milk transcriptome corresponding at colostrum and 120 days of lactation. To reconstruct transcripts we used both the genome as reference, and a de novo assembly approach. Additionally, we aimed at identifying the differentially expressed genes (DEGs) between the two lactation stages and at analyzing the expression of genes involved in oligosaccharides metabolism.

RESULTS

A total of 44,635 different transcripts, organized in 33,757 tentative genes, were obtained using the goat genome as reference. A significant sequence similarity match was found for 40,353 transcripts (90%) against the NCBI NT and for 35,701 (80%) against the NR databases. 68% and 69% of the de novo assembled transcripts, in colostrum and 120 days of lactation samples respectively, have a significant match with the merged transcriptome obtained using Cufflinks/Cuffmerge. CSN2, PAEP, CSN1S2, CSN3, LALBA, TPT1, FTH1, M-SAA3, SPP1, GLYCAM1, EEF1A1, CTSD, FASN, RPS29, CSN1S1, KRT19 and CHEK1 were found between the top fifteen highly expressed genes. 418 loci were differentially expressed between lactation stages, among which 207 and 122 were significantly up- and down-regulated in colostrum, respectively. Functional annotation and pathway enrichment analysis showed that in goat colostrum somatic cells predominate biological processes involved in glycolysis, carbohydrate metabolism, defense response, cytokine activity, regulation of cell proliferation and cell death, vasculature development, while in mature milk, biological process associated with positive regulation of lymphocyte activation and anatomical structure morphogenesis are enriched. The analysis of 144 different oligosaccharide metabolism-related genes showed that most of these (64%) were more expressed in colostrum than in mature milk, with eight expressed at very high levels (SLCA3, GMSD, NME2, SLC2A1, B4GALT1, B3GNT2, NANS, HEXB).

CONCLUSIONS

To our knowledge, this is the first study comparing goat transcriptome of two lactation stages: colostrum and 120 days. Our findings suggest putative differences of expression between stages and can be envisioned as a base for further research in the topic. Moreover because a higher expression of genes involved in immune defense response, carbohydrate metabolism and related to oligosaccharide metabolism was identified in colostrum we here corroborate the potential of goat milk as a natural source of lactose-derived oligosaccharides and for the development of functional foods.