Evaluation of milk somatic cells as a source of mRNA for study of lipogenesis in the mammary gland of lactating beef cows supplemented with dietary high-linoleate safflower seeds

Effect of Dietary Eicosapentaenoic and Docosahexaenoic Fatty Acid Supplementation during the Last Month of Gestation on Fatty Acid Metabolism and Oxidative Status in Charolais Cows and Calves

Fatty acids (FAs) are of utmost importance in the peripartal period for the development of the central nervous and immune systems of the newborn. The transport of polyunsaturated fatty acids (PUFAs) through the placenta is considered to be minimal in ruminants. Nevertheless, the cow's FAs are the main source of FAs for the calf during gestation. This research aimed to investigate the influence of low-dose eicosapentaenoic acid (EPA) and docosahexaenoic acid (DHA) supplementation during late gestation on the FA metabolism of cows and their calves. A total of 20 Charolais cows during the last month of their gestation were included in the feeding trial and were divided into a control group (CON) and an experimental group (EPA + DHA). The latter received a supplement in the amount of 100 g/day (9.1 and 7.8 g/cow/day of EPA and DHA, respectively). Supplementation of low-dose EPA and DHA alters colostrum and milk fatty acid composition through the elevation of n-3 long-chain polyunsaturated fatty acids (LC-PUFAs) without affecting milk fat and protein concentrations and oxidative status. Plasma composition in cows was significantly altered, while the same effect was not detected in calf plasma. No significant change in mRNA expression was detected for the genes fatty acid synthase (FASN) and acetyl-CoA carboxylase alpha (ACACA).

Epigenetic regulation of functional candidate genes for milk production traits in dairy sheep subjected to protein restriction in the prepubertal stage

BACKGROUND

As the prepubertal stage is a crucial point for the proper development of the mammary gland and milk production, this study aims to evaluate how protein restriction at this stage can affect methylation marks in milk somatic cells. Here, 28 Assaf ewes were subjected to 42.3% nutritional protein restriction (14 animals, NPR) or fed standard diets (14 animals, C) during the prepubertal stage. During the second lactation, the milk somatic cells of these ewes were sampled, and the extracted DNA was subjected to whole-genome bisulfite sequencing.

RESULTS

A total of 1154 differentially methylated regions (DMRs) were identified between the NPR and C groups. Indeed, the results of functional enrichment analyses of the genes harboring these DMRs suggested their relevant effects on the development of the mammary gland and lipid metabolism in sheep. The additional analysis of the correlations of the mean methylation levels within these DMRs with fat, protein, and dry extract percentages in the milk and milk somatic cell counts suggested associations between several DMRs and milk production traits. However, there were no phenotypic differences in these traits between the NPR and C groups.

CONCLUSION

In light of the above, the results obtained in the current study might suggest potential candidate genes for the regulation of milk production traits in the sheep mammary gland. Further studies focusing on elucidating the genetic mechanisms affected by the identified DMRs may help to better understand the biological mechanisms modified in the mammary gland of dairy sheep as a response to nutritional challenges and their potential effects on milk production.

Effects of PUFA-Rich Dietary Strategies on Ruminants' Mammary Gland Gene Network: A Nutrigenomics Review

Although the inclusion of polyunsaturated fatty acids (PUFAs) in ruminants' diets appears to be a well-documented strategy to enrich milk with PUFAs, several gene networks that regulate milk synthesis and mammary gland homeostasis could be impaired. The objective of this literature review is to assess the effects of nutritional strategies focused on enriching milk with PUFAs on gene networks regulating mammary gland function and lipogenesis, as well as the impact of feed additives and bioactive compounds with prominent antioxidant potential on immune-oxidative transcriptional profiling, as a part of mammary gland homeostasis and health. The findings support the conclusion that PUFAs' inclusion in ruminants' diets more strongly downregulate the stearoyl-CoA desaturase (SCD) gene compared to other key genes involved in de novo fatty acid synthesis in the mammary gland. Additionally, it was revealed that seed oils rich in linoleic and linolenic acids have no such strong impact on networks that regulate lipogenic homeostasis compared to marine oils rich in eicosapentaenoic and docosahexaenoic acids. Furthermore, ample evidence supports that cows and sheep are more prone to the suppression of lipogenesis pathways compared to goats under the impact of dietary marine PUFAs. On the other hand, the inclusion of feed additives and bioactive compounds with prominent antioxidant potential in ruminants' diets can strengthen mammary gland immune-oxidative status. Considering that PUFA's high propensity to oxidation can induce a cascade of pro-oxidant incidences, the simultaneous supplementation of antioxidant compounds and especially polyphenols may alleviate any side effects caused by PUFA overload in the mammary gland. In conclusion, future studies should deeply investigate the effects of PUFAs on mammary gland gene networks in an effort to holistically understand their impact on both milk fat depression syndrome and homeostatic disturbance.

Single nucleotide polymorphisms in the bovine SLC2A12 and SLC5A1 glucose transporter genes - the effect on gene expression and milk traits of Holstein Friesian cows

In this study, novel single nucleotide polymorphisms (SNPs) were found in the 5'-regulatory regions (promoters) of the bovine glucose transporter (GT) genes SLC2A12 and SLC5A1. These polymorphisms were shown to associate with certain milk production traits in HF cows, including milk yield, milk composition, and somatic cell count. It was shown that the SNP g.-671C > G (NC_037336.1: g.72224078C > G) in the SLC2A12 gene could be an effective marker of cattle production traits and that genotypes CC and CG are associated with the best productivity. The polymorphisms found in the SLC5A1 gene promoter also influenced milk production traits in HF cows, albeit to a lesser extent, and we propose that these polymorphisms could be useful as genetic markers for milk production traits in marker-assisted selection (MAS); however, this must be confirmed on larger populations of cattle. In addition, the presence of polymorphisms within promoter regions appears to affect the expression of GT genes in the cow mammary gland and modify transcription factor (TF) binding capacity.

The effect of single-nucleotide polymorphism in the promoter region of bovine alpha-lactalbumin (LALBA) gene on LALBA expression in milk cells and milk traits of cows

Polymorphisms of milk protein genes have been proposed as candidate markers for dairy production traits in cattle. In the present study, a polymorphism was detected in the 5'-flanking (promoter) region of the bovine alpha-lactalbumin (LALBA) gene, a T/C transition located at nucleotide -1,001 relative to the transcription start site g.-1001T > C (NC_037332.1:g.31183170T > C), which is recognizable with PstI restriction endonuclease. In silico analyses showed that this mutation created novel retinoid X receptor alpha and vitamin D receptor transcription factor binding sites. Real-time PCR found that cows with different genetic variants of the promoter demonstrated different levels of expression of LALBA mRNA in milk somatic cells (MSCs). The TT genotype cows demonstrated low expression, whereas those with CT demonstrated much higher expression (P < 0.05). ELISA analysis found milk LALBA protein levels also differed between the TT and CT cows (P < 0.05) and that these levels were not correlated with the mRNA abundance in MSC. Association analysis found that the g.-1001T > C polymorphism in the promoter region of the LALBA gene influenced milk production traits in Polish Holstein-Friesian cows. High daily milk yield and dry matter yield, and high lactose yield and concentration were associated with the TT genotype. The TT genotype cows also had a lower number of somatic cells in the milk, considered as an indicator of udder health status. Therefore, the TT genotype could be more desirable from the breeder's perspective.

Intramammary inoculation with lactic acid bacteria at dry-off triggers an immunomodulatory response in dairy cows

The use of antibiotics to prevent bovine mastitis is responsible for the emergence and selection of resistant strains. Lactic acid bacteria (LAB) could be introduced into animal feed as an alternative prevention method that would bypass the risk of resistance development. In previous research, we demonstrated that two probiotic LAB strains isolated from bovine milk were capable of stimulating the production of antibodies and the host's immune cellular response in the udder. The present study aimed to elucidate whether the antibodies of animals inoculated with these strains were able to increase phagocytosis by neutrophils and inhibit the growth of different mastitis-causing pathogens. Moreover, the effect of LAB on the expression of pro-inflammatory cytokines was assessed. Ten animals were inoculated intramammarily with 10⁶ cells of the two strains at dry-off. The blood serum was tested for its ability to opsonize bovine mastitis pathogens, the in vitro bactericidal activity of bovine blood and milk against these pathogens was determined, and cytokine mRNA expression was quantified in milk somatic cells. The inoculated animals did not show abnormal signs of sensitivity to the LAB. Their blood serum significantly enhanced the phagocytosis of Staphylococcus spp. and the LAB. Escherichia coli and Streptococcus uberis were inhibited by the milk serum but not the blood serum, whereas Staphylococcus aureus and Staphylococcus haemolyticus were inhibited by both. In regard to cytokine expression, interleukin (IL)-1β increased markedly for up to 4 h post-inoculation, and an increase in IL-8 was observed 4, 12 and 24 h after inoculation. Tumour necrosis factor-α mRNA increased 1 and 2 h after inoculation and a significant difference was registered at 6 h for interferon-γ. This rapid immunomodulatory response shows that inoculating animals with LAB at dry-off, when they are especially susceptible, could be a useful strategy for the prevention of bovine mastitis.

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.

Transcript profiling in the milk of dairy ewes fed extruded linseed

The identification of genes regulating milk secretion by mammary glands is a key-step for exploiting changes in milk composition induced by different diet regimens. Linseed supplementation is a reliable feeding strategy to enhance polyunsaturated fatty acid content in milk fat from sheep, increasing milk quality and nutraceutical value. To investigate the molecular bases of diet induced differences in milk composition, we collected milk from dairy ewes at 90 day of lactation and after 3 weeks of diet supplementation with extruded linseed. The milk of dairy ewes contains milk somatic cells mostly derived by mammary glands. After isolating milk somatic cells by centrifugation, RNAs were purified from these cells, and Illumina RNA sequencing was performed to analyze RNA synthesis. Our data provide a resource (available at Gene Expression Omnibus database under GSE89163) to be employed for comparative analyses of gene expression in milk somatic cells in different breeds and different diets, with the long-term aim of developing strategies to improve sheep milk quality.

Dietary linseed oil with or without malate increases conjugated linoleic acid and oleic acid in milk fat and and gene expression in mammary gland and milk somatic cells of lactating goats

Supplementary dietary plant oils have the potential to alter milk fatty acid composition in ruminants as a result of changes in the amount and kind of fatty acid precursors. We hypothesized that linseed oil in combination with malate (a key propionate precursor in the rumen) would increase ∆9 unsaturated fatty acids and specific gene expression in somatic cells and mammary glands of lactating goats. Twelve lactating goats were used in a 3 × 3 Latin square design. Treatments included the basal diet (CON), the CON plus 4% linseed oil (LO), and the CON plus 4% linseed oil and 2% -malate (LOM). Relative to CON, the LO and LOM supplements increased the daily intake of palmitic (16:0), stearic (18:0), oleic (18:1-9), linoleic (18:2-6), α-linolenic (18:3-3), and γ-linolenic acids (18:2-6); α-linolenic acid intake was increased over 9-fold, from 6.77 to over 51 g/d ( < 0.02). The LO and LOM supplements increased daily milk yield, milk fat yield, and milk fat percentage ( < 0.05). The LOM supplement also increased milk lactose percentage and daily yield ( = 0.03). Both the LO and LOM supplements increased plasma glucose and total cholesterol and decreased plasma β-hydroxbutyrate concentrations ( = 0.03). The LO and LOM supplements increased concentrations of stearic acid; -vaccenic acid (TVA; 18:1-11); -9, -11 CLA; -10 -12 CLA; and α-linolenic acid in rumen fluid and increased the concentrations of oleic acid; TVA; -9, -11 CLA; -10, -12 CLA; and α-linolenic acid in plasma lipids and milk fat ( < 0.05). Conversely, the LO and LOM supplements decreased short- and medium-chain SFA, including lauric (12:0), myristic (14:0), and palmitic acids, in plasma and milk fat ( < 0.05). Relative mRNA levels for and () gene expression were increased in somatic cells and mammary gland tissue by LO and LOM ( < 0.05). We conclude that the higher intake and ruminal production of stearic acid promoted SCD gene expression in somatic cells and mammary tissue. Furthermore, milk somatic cells are a suitable substitute for documenting treatment effects of dietary oils on gene expression in goat mammary tissue.

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

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

Milk fat globule is an alternative to mammary epithelial cells for gene expression analysis in buffalo

Owing to the difficulty in obtaining mammary gland tissue from lactating animals, it is difficult to test the expression levels of genes in mammary gland. The aim of the current study was to identify if milk fat globule (MFG) in buffalo milk was an alternative to mammary gland (MG) and milk somatic cell (MSC) for gene expression analysis. Six buffalos in late lactation were selected to collect MFG and MSC, and then MG was obtained by surgery. MFG was stained with acridine orange to successfully visualise RNA and several cytoplasmic crescents in MFG. The total RNA in MFG was successfully isolated and the integrity was assessed by agarose gel electrophoresis. We analysed the cellular components in MFG, MG and MSC through testing the expression of cell-specific genes by qRT-PCR. The results showed that adipocyte-specific gene (AdipoQ) and leucocyte-specific genes (CD43, CSF1 and IL1α) in MFG were not detected, whereas epithelial cell marker genes (Keratin 8 and Keratin 18) in MFG were higher than in MSC and lower than in MG, fibroblast marker gene (vimentin) in MFG was significantly lower than in MG and MSC, milk protein genes (LALBA, BLG and CSN2) and milk fat synthesis-related genes (ACC, BTN1A1, FABP3 and FAS) in MFG were higher than in MG and MSC. In conclusion, the total RNA in MFG mainly derives from mammary epithelial cells and can be used to study the functional gene expression of mammary epithelial cells.

Advances in research on cis-9, trans-11 conjugated linoleic acid: a major functional conjugated linoleic acid isomer

Conjugated linoleic acid (CLA) consists of a group of positional and geometric conjugated isomers of linoleic acid. Since the identification of CLA as a factor that can inhibit mutagenesis and carcinogenesis, thousands of studies have been conducted in the last several decades. Among the many isomers discovered, cis-9, trans-11 CLA is the most intensively studied because of its multiple, isomer-specific effects in humans and animals. This paper provides an overview of the available data on cis-9, trans-11 CLA, including its isomer-specific effects, biosynthesis, in vivo/in vitro research models, quantification, and the factors influencing its content in ruminant products.

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

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

Screening for the Most Suitable Reference Genes for Gene Expression Studies in Equine Milk Somatic Cells

Apart from the well-known role of somatic cell count as a parameter reflecting the inflammatory status of the mammary gland, the composition of cells isolated from milk is considered as a valuable material for gene expression studies in mammals. Due to its unique composition, in recent years an increasing interest in mare's milk consumption has been observed. Thus, investigating the genetic background of horse's milk variability presents and interesting study model. Relying on 39 milk samples collected from mares representing three breeds (Polish Primitive Horse, Polish Cold-blooded Horse, Polish Warmblood Horse) we aimed to investigate the utility of equine milk somatic cells as a source of mRNA and to screen the best reference genes for RT-qPCR using geNorm and NormFinder algorithms. The results showed that despite relatively low somatic cell counts in mare's milk, the amount and the quality of the extracted RNA are sufficient for gene expression studies. The analysis of the utility of 7 potential reference genes for RT-qPCR experiments for the normalization of equine milk somatic cells revealed some differences between the outcomes of the applied algorithms, although in both cases the KRT8 and TOP2B genes were pointed as the most stable. Analysis by geNorm showed that the combination of 4 reference genes (ACTB, GAPDH, TOP2B and KRT8) is required for apropriate RT-qPCR experiments normalization, whereas NormFinder algorithm pointed the combination of KRT8 and RPS9 genes as the most suitable. The trial study of the relative transcript abundance of the beta-casein gene with the use of various types and numbers of internal control genes confirmed once again that the selection of proper reference gene combinations is crucial for the final results of each real-time PCR experiment.

Dietary sunflower oil modulates milk fatty acid composition without major changes in adipose and mammary tissue fatty acid profile or related gene mRNA abundance in sheep

There are very few studies in ruminants characterizing mammary and adipose tissue (AT) expression of genes and gene networks for diets causing variations in milk fatty acid (FA) composition without altering milk fat secretion, and even less complementing this information with data on tissue FA profiles. This work was conducted in sheep in order to investigate the response of the mammary gland and the subcutaneous and perirenal AT, in terms of FA profile and mRNA abundance of genes involved in lipid metabolism, to a diet known to modify milk FA composition. Ten lactating Assaf ewes were randomly assigned to two treatments consisting of a total mixed ration based on alfalfa hay and a concentrate (60 : 40) supplemented with 0 (control diet) or 25 (SO diet) g of sunflower oil/kg of diet dry matter for 7 weeks. Milk composition, including FA profile, was analysed after 48 days on treatments. On day 49, the animals were euthanized and tissue samples were collected to analyse FA and mRNA abundance of 16 candidate genes. Feeding SO did not affect animal performance but modified milk FA composition. Major changes included decreases in the concentration of FA derived from de novo synthesis (e.g. 12:0, 14:0 and 16:0) and increases in that of long-chain FA (e.g. 18:0, c9-18:1, trans-18:1 isomers and c9,t11-CLA); however, they were not accompanied by significant variations in the mRNA abundance of the studied lipogenic genes (i.e. ACACA, FASN, LPL, CD36, FABP3, SCD1 and SCD5) and transcription factors (SREBF1 and PPARG), or in the constituent FA of mammary tissue. Regarding the FA composition of AT, the little influence of SO did not appear to be linked to changes in gene mRNA abundance (decreases of GPAM and SREBF1 in both tissues, and of PPARG in the subcutaneous depot). Similarly, the great variation between AT (higher contents of saturated FA and trans-18:1 isomers in the perirenal, and of cis-18:1, c9,t11-CLA and n-3 PUFA in the subcutaneous AT) could not be related to differences in gene mRNA abundance due to tissue site (higher LPL and CD36, and lower SREBF1 in perirenal than in subcutaneous AT). Overall, these results suggest a marginal contribution of gene expression to the nutritional regulation of lipid metabolism in these tissues, at least with the examined diets and after 7 weeks on treatments. It cannot be ruled out, however, that the response to SO is mediated by other genes or post-transcriptional mechanisms.

Expression profiling of glucose transporter 1 (GLUT1) and apoptotic genes (BAX and BCL2) in milk enriched mammary epithelial cells (MEC) in riverine buffalo during lactation

Lactation is an important physiological process in dairy animals. During lactation, up to 85% of the body glucose is directed toward the mammary glands for milk synthesis. Studies related to lactation physiology are generally carried out on mammary biopsies, which may adversely affect animal health. In the present study, milk enriched MEC were used to study the expression pattern of GLUT1 and apoptotic genes (BAX and BCL2) across different stages of lactation in riverine buffalo in relation to milk yield. MEC were enriched from milk using cytokeratin-8 antibodies coated magnetic beads. Total RNA isolated from enriched MEC showed significant correlation (r(2) = 0.92 ± 0.02) with the milk yield at different stages of lactation. GLUT1 expression pattern correlated with the milk yield as highest GLUT1 expression (4.68 ± 0.79) was observed during peak-lactation (90 days post-parturition), whereas low GLUT1 expression (1.01 ± 0.1, 15 d; 0.71 ± 0.03, 30 d) was observed during early lactation. The BAX/BCL2 ratio was high (1.02 ± 0.2, 15 d; 0.94 ± 0.06, 30 d) during the early phase of lactation, indicating high rate of apoptosis, whereas low BAX/BCL2 ratio (0.25 ± 0.03, 60 d; 0.42 ± 0.04, 90 d) was observed during mid-lactation coinciding with the increase in RNA concentration and milk yield. Highest BAX/BCL2 ratio (1.41 ± 0.3, 120 d; 4.02 ± 0.6, 240 d) was observed during late lactation i.e., 240 days, which was also reflected as decline in milk yield and RNA concentration. Also, BAX/BCL2 ratio in milk enriched MEC was in accordance with RNA concentration in MEC and milk yield at different phases of lactation. Our study showed that expression pattern of genes under study (GLUT1, BAX, and BCL2) in milk enriched MEC correlated well with important physiological properties such as milk yield in buffalo.

Technical note: Selection of suitable reference genes for studying gene expression in milk somatic cell of yak (Bos grunniens) during the lactation cycle

Quantitative real-time PCR is the most sensitive technique for gene expression analysis. Data normalization is essential to correct for potential errors incurred in all steps from RNA isolation to PCR amplification. The commonly accepted approach for normalization is the use of reference gene. Until now, no suitable reference genes have been available for data normalization of gene expression in milk somatic cells of lactating yaks across lactation. In the present study, we evaluated the transcriptional stability of 10 candidate reference genes in milk somatic cells of lactating yak, including ACTB, B2M, GAPDH, GTP, MRPL39, PPP1R11, RPS9, RPS15, UXT, and RN18S1. Four genes, RPS9, PPP1R11, UXT, and MRPL39, were identified as being the most stable genes in milk somatic cells of lactating yak. Using the combination of RPS9, PPP1R11, UXT, and MRPL39 as reference genes, we further assessed the relative expression of 4 genes of interest in milk somatic cells of yak across lactation, including ELF5, ABCG2, SREBF2, and DGAT1. Compared with expression in colostrum, the overall transcription levels of ELF5, ABCG2, and SREBF2 in milk were found to be significantly upregulated in early, peak, and late lactation, and significantly downregulated thereafter, before the dry period. A similar pattern was observed in the relative expression of DGAT1, but no significant difference was revealed in its expression in milk from late lactation compared with colostrum. Based on these results, we suggest that the geometric mean of RPS9, PPP1R11, UXT, and MRPL39 can be used for normalization of real-time PCR data in milk somatic cells of lactating yak, if similar experiments are performed.

Unilateral once daily milking locally induces differential gene expression in both mammary tissue and milk epithelial cells revealing mammary remodeling

Once daily milking reduces milk yield, but the underlying mechanisms are not yet fully understood. Local regulation due to milk stasis in the tissue may contribute to this effect, but such mechanisms have not yet been fully described. To challenge this hypothesis, one udder half of six Holstein dairy cows was milked once a day (ODM), and the other twice a day (TDM). On the 8th day of unilateral ODM, mammary epithelial cells (MEC) were purified from the milk using immunomagnetic separation. Mammary biopsies were harvested from both udder halves. The differences in transcript profiles between biopsies from ODM and TDM udder halves were analyzed by a 22k bovine oligonucleotide array, revealing 490 transcripts that were differentially expressed. The principal category of upregulated transcripts concerned mechanisms involved in cell proliferation and death. We further confirmed remodeling of the mammary tissue by immunohistochemistry, which showed less cell proliferation and more apoptosis in ODM udder halves. Gene expression analyzed by RT-qPCR in MEC purified from milk and mammary biopsies showed a common downregulation of six transcripts (ABCG2, FABP3, NUCB2, RNASE1 and 5, and SLC34A2) but also some discrepancies. First, none of the upregulated transcripts in biopsies varied in milk-purified MEC. Second, only milk-purified MEC showed significant LALBA downregulation, which suggests therefore that they correspond to a mammary epithelial cell subpopulation. Our results, obtained after unilateral milking, suggest that cell remodeling during ODM is due to a local effect, which may be triggered by milk accumulation.

Development of a RNA extraction method from milk for gene expression study in the mammary gland of sheep

The aim of the study was to develop a reliable method for the RNA extraction from milk of Sarda sheep breed and to highlight if the extracted RNA can be used for expression study on mammary genes involved in milk fat synthesis using RT-qPCR. The main result is that a sample of 150 ml of milk provides an optimal amount of RNA (73.5 μg/ml). The highest RNA concentration has been found in the samples analysed within 4 h after collection. The RNA extracted was positively correlated to the number of somatic cells (P < 0.001). The efficiency of the extraction method was confirmed by the results obtained from qPCR which showed a Ct value, for SREBPF1 gene of 26.8 ± 0.15. This research demonstrated that the high-quality of the RNA obtained is suited to use for studies of mammary genes expression in sheep, avoiding any damage caused by mammary gland biopsy.

Selection of suitable reference genes for quantitative gene expression studies in milk somatic cells of lactating cows (Bos indicus)

We assessed the suitability of 9 internal control genes (ICG) in milk somatic cells of lactating cows to find suitable reference genes for use in quantitative PCR (qPCR). Eighteen multiparous lactating Sahiwal cows were used, 6 in each of 3 lactation stages: early (25 ± 5 d in milk), mid (160 ± 15 d in milk), and late (275 ± 25 d in milk) lactation. Nine candidate reference genes [glyceraldehyde 3-phosphate dehydrogenase (GAPDH), protein phosphatase 1 regulatory subunit 11 (PPP1R11), β-actin (ACTB), β-2 microglobulin (B2M), 40S ribosomal protein S15a (RPS15A), ubiquitously expressed transcript (UXT), mitochondrial GTPase 1 (MTG1), 18S rRNA (RN18S1), and ubiquitin (UBC)] were evaluated. Three genes, β-casein (CSN2), lactoferrin (LTF), and cathelicidin (CAMP) were chosen as target genes. Very high amplification was observed in 7 ICG and very low level amplification was observed in 2 ICG (UXT and MTG1). Thus, UXT and MTG1 were excluded from further analysis. The qPCR data were analyzed by 2 software packages, geNorm and NormFinder, to determine suitable reference genes, based on their stability and expression. Overall, PPP1R11, ACTB, UBC, and GAPDH were stably expressed among all candidate reference genes. Therefore, these genes could be used as ICG for normalization of qPCR data in milk somatic cells through lactation.

Use of milk epithelial cells to study regulation of cell activity and apoptosis during once-daily milking in goats

Generally, once-daily milking (ODM) decreases milk yield. This effect may be the consequence of a decrease in mammary epithelial cell (MEC) activity or a reduction in their number. The aim of this study was to determine the effect of ODM on the synthetic activity and rate of apoptosis of MEC using a non-invasive method. Eight Alpine goats were subjected to ODM or twice-daily milking for two 5-week periods. MECs were purified by centrifugation and immunocytochemical binding in milk after 1 and 5 weeks of each period. mRNA levels of some proteins involved in lactose and milk protein synthesis and in apoptosis were evaluated using real-time PCR. Isolation of MEC from milk was a useful method to investigate transcriptional regulation in a timeline study. ODM induced greater decreases in milk, lactose and protein yields after 1 week than after 5 weeks. This suggests an adaptation of the mammary gland to ODM, which reduces the inhibitory effect of this practice. Reductions in milk component yields were associated with lower α-lactalbumin transcripts, suggesting a transcriptional decrease of lactose synthesis during ODM. Glucose transporter GLUT1 transcripts were downregulated under ODM, suggesting that lactose precursor uptake by MEC might be involved in the regulation of lactose synthesis. κ-Casein mRNA levels tended to be lower during ODM. ODM increased levels of the pro-apoptotic transcript Bax after both 1 and 5 weeks, but no variation was observed in the Bax/Bcl-2 ratio. ODM affected cell synthetic activity through transcriptional regulation and may have induced apoptosis. The reduction of the negative effect of ODM on milk yield suggests that Alpine goats are able to adapt to ODM. Further studies are needed to investigate the effect of ODM on MEC turnover.

Alteration of gene expression in mammary gland tissue of dairy cows in response to dietary unsaturated fatty acids

The aim of this study was to determine the effects of supplementing unprotected dietary unsaturated fatty acids (UFAs) from different plant oils on gene expression in the mammary gland of grazing dairy cows. A total of 28 Holstein-Friesian dairy cows in mid-lactation were blocked according to parity, days in milk, milk yield and fat percentage. The cows were then randomly assigned to four UFA sources based on rapeseed, soybean, linseed or a mixture of the three oils for 23 days, after which, all 28 cows were switched to a control diet for an additional 28 days. On the last day of both periods, mammary gland biopsies were taken to study genome-wide differences in gene expression on Affymetrix GeneChip® Bovine Genome Arrays (no. 900493) by ServiceXS (Leiden, The Netherlands). Supplementation with UFAs resulted in increased milk yield but decreased milk fat and protein percentages. Furthermore, the proportion of de novo fatty acids (FAs) in the milk was reduced, whereas that of long-chain FAs increased. Applying a statistical cut-off of false discovery rate of q-values <0.05 together with an absolute fold change of 1.3, a total of 972 genes were found to be significantly affected through UFA supplementation, indicating that large transcriptional adaptations occurred in the mammary gland when grazing dairy cows were supplemented with unprotected dietary UFA. Gene sets related to cell development and remodeling, apoptosis, nutrient metabolic process, as well as immune system response were predominantly downregulated during UFA supplementation. Such molecular knowledge on the physiology of the mammary gland might provide the basis for further functional research on dairy cows.

Evaluation of a replacement method for mammary gland biopsies by comparing gene expression in udder tissue and mammary epithelial cells isolated from milk

Somatic cells isolated from milk offer an attractive non-invasive replacement of invasive udder biopsies for monitoring bovine mammary gland metabolism. However, for metabolic gene expression studies the mammary gland epithelial cells (MEC) isolated from milk have to be purified from the non-epithelial leukocyte fraction in milk samples. In our study, enrichment of MEC by using anti-cytokeratin peptide 18 (KRT18) antibody coated magnetic beads was evaluated. MEC showed a substantially increased expression of the epithelial-cell-specific KRT18 gene compared to udder tissue. The expression levels of genes specific for mammary gland epithelial cells (CSN3 and LALBA) showed a significant positive correlation in MEC and also in udder tissue. However, no significant correlation of the expression of a specific gene was found between udder and MEC samples. Therefore, MEC isolated from total milk samples via KRT18 antibodies probably do not reflect the true metabolic situation of the bovine udder. Thus, quantitative gene expression profiling of MEC isolated via KRT18 antibodies has to be interpreted carefully with respect to the situation in the udder.

Maintaining RNA integrity in a homogeneous population of mammary epithelial cells isolated by Laser Capture Microdissection

BACKGROUND

Laser-capture microdissection (LCM) that enables the isolation of specific cell populations from complex tissues under morphological control is increasingly used for subsequent gene expression studies in cell biology by methods such as real-time quantitative PCR (qPCR), microarrays and most recently by RNA-sequencing. Challenges are i) to select precisely and efficiently cells of interest and ii) to maintain RNA integrity. The mammary gland which is a complex and heterogeneous tissue, consists of multiple cell types, changing in relative proportion during its development and thus hampering gene expression profiling comparison on whole tissue between physiological stages. During lactation, mammary epithelial cells (MEC) are predominant. However several other cell types, including myoepithelial (MMC) and immune cells are present, making it difficult to precisely determine the specificity of gene expression to the cell type of origin. In this work, an optimized reliable procedure for producing RNA from alveolar epithelial cells isolated from frozen histological sections of lactating goat, sheep and cow mammary glands using an infrared-laser based Arcturus Veritas LCM (Applied Biosystems®) system has been developed. The following steps of the microdissection workflow: cryosectioning, staining, dehydration and harvesting of microdissected cells have been carefully considered and designed to ensure cell capture efficiency without compromising RNA integrity.

RESULTS

The best results were obtained when staining 8 μm-thick sections with Cresyl violet® (Ambion, Applied Biosystems®) and capturing microdissected cells during less than 2 hours before RNA extraction. In addition, particular attention was paid to animal preparation before biopsies or slaughtering (milking) and freezing of tissue blocks which were embedded in a cryoprotective compound before being immersed in isopentane. The amount of RNA thus obtained from ca.150 to 250 acini (300,000 to 600,000 μm2) ranges between 5 to 10 ng. RNA integrity number (RIN) was ca. 8.0 and selectivity of this LCM protocol was demonstrated through qPCR analyses for several alveolar cell specific genes, including LALBA (α-lactalbumin) and CSN1S2 (αs2-casein), as well as Krt14 (cytokeratin 14), CD3e and CD68 which are specific markers of MMC, lymphocytes and macrophages, respectively.

CONCLUSIONS

RNAs isolated from MEC in this manner were of very good quality for subsequent linear amplification, thus making it possible to establish a referential gene expression profile of the healthy MEC, a useful platform for tumor biomarker discovery.

Milking and feed restriction regulate transcripts of mammary epithelial cells purified from milk

Feed restriction and once-daily milking (ODM) reduce milk yield in dairy cows and the amount of glucose taken up by the mammary gland. The modulation of mammary glucose uptake may be the consequence of modifications to glucose transport, capacity for lactose synthesis, and cell death in mammary epithelial cells (MEC). The aim was to demonstrate the usefulness of a new method to purify MEC from milk somatic cells and to examine the effects of feed restriction and ODM on mammary transcripts. Five Holstein cows were subjected to a 2 x 2 factorial arrangement of 2 milking frequencies and 2 feeding levels, during which the cows were milked once or twice daily while fed a diet providing either 98 or 70% of requirements. The cows were equipped to study net mammary balance of glucose. On d 7 of each experimental week, milk and lactose yields and mammary glucose uptake were measured. Cells were isolated from fresh milk by centrifugation to generate total milk cell samples. Mammary epithelial cells were separated from total milk cells by using magnetic beads associated with anticytokeratin 8 antibodies. Total RNA was extracted from both total milk cells and purified MEC samples. Real-time reverse transcription PCR was performed to determine mRNA levels in purified MEC under feed restriction and under ODM. Purified MEC samples revealed higher total RNA quality (RNA integrity number = 8) and were better suited to the measurement of mammary transcripts than total milk cell samples (RNA integrity number = 4). Significant correlations were obtained between mRNA levels and net glucose balance data (0.465 < r < 0.680), demonstrating the validity of results obtained by using purified MEC. Feed restriction induced a significant reduction (by half) in type 1 glucose transporter mRNA levels without any effect on alpha-lactalbumin (alpha-LA), galactosyltransferase, kappa-casein, bcl2, or bax mRNA levels. When compared with twice daily milking, ODM reduced kappa-casein (-86%) and alphaLA (-73%) mRNA levels and up-regulated bax and bcl2 mRNA levels (7- and 9-fold). The results suggest that the regulation of glucose uptake and milk yield is dependent on the transcription of glucose transporters under feed restriction and on the transcription of alphaLA under ODM. Further studies are required to con-firm the suggested onset of cell death after 7 d of ODM.

Expression and nutritional regulation of lipogenic genes in the ruminant lactating mammary gland

The effect of nutrition on milk fat yield and composition has largely been investigated in cows and goats, with some differences for fatty acid (FA) composition responses and marked species differences in milk fat yield response. Recently, the characterization of lipogenic genes in ruminant species allowed in vivo studies focused on the effect of nutrition on mammary expression of these genes, in cows (mainly fed milk fat-depressing diets) and goats (fed lipid-supplemented diets). These few studies demonstrated some similarities in the regulation of gene expression between the two species, although the responses were not always in agreement with milk FA secretion responses. A central role for trans-10 C18:1 and trans-10, cis-12 CLA as regulators of milk fat synthesis has been proposed. However, trans-10 C18:1 does not directly control milk fat synthesis in cows, despite the fact that it largely responds to dietary factors, with its concentration being negatively correlated with milk fat yield response in cows and, to a lesser extent, in goats. Milk trans-10, cis-12 CLA is often correlated with milk fat depression in cows but not in goats and, when postruminally infused, acts as an inhibitor of the expression of key lipogenic genes in cows. Recent evidence has also proven the inhibitory effect of the trans-9, cis-11 CLA isomer. The molecular mechanisms by which nutrients regulate lipogenic gene expression have yet to be well identified, but a central role for SREBP-1 has been outlined as mediator of FA effects, whereas the roles of PPARs and STAT5 need to be determined. It is expected that the development of in vitro functional systems for lipid synthesis and secretion will allow future progress toward (1) the identification of the inhibitors and activators of fat synthesis, (2) the knowledge of cellular mechanisms, and (3) the understanding of differences between ruminant species.

Extraction and Quantitative Analysis of Stearoyl-Coenzyme A Desaturase mRNA from Dairy Cow Milk Somatic Cells

Study of stearoyl-coenzyme A desaturase (SCD) gene expression in the bovine mammary gland is limited by restricted availability of mammary tissue samples from biopsy or postmortem sampling of cows during temporal experiments. A technique was developed to isolate total RNA from somatic cells in bovine milk and to analyze SCD mRNA expression by quantitative reverse-transcription PCR. Total RNA yield was lower than in a previous goat study and was related to numbers of viable somatic cells. To obtain sufficient total RNA, 1-L milk samples were taken and stored for up to 24 h at 4 degrees C. Complementary DNA prepared from somatic cells showed a 99% match with the published sequence for SCD mRNA in bovine adipose tissue. Stearoyl-CoA desaturase mRNA abundance relative to beta-actin mRNA for 12 cows sampled across 4 time points varied (mean +/- SE) from 0.88 +/- 0.17 to 4.40 +/- 0.50. Fifty-five percent of variation was due to individual cows and 42% was due to daily variation within cows. Relative abundance of SCD mRNA was not related to the number of viable somatic cells or total RNA extracted from samples, but it was related to mammary desaturase activity, as indicated by changes in milk C14 fatty acid concentrations. We concluded that somatic cells provide a noninvasive and repeatable alternative to mammary tissue samples obtained by biopsy or postmortem.

Effects of postpartum dietary fat and body condition score at parturition on plasma, adipose tissue, and milk fatty acid composition of lactating beef cows

Two experiments were conducted with lactating Angus x Gelbvieh beef cows to determine the effects of postpartum lipid supplementation, BCS at parturition, and day of lactation on fatty acid profiles in plasma, adipose tissue, and milk. In Exp. 1, 36 pri-miparous cows (488 +/- 10 kg of initial BW; 5.5 +/- 0.02 initial BCS) were given ad libitum access to hay and assigned randomly to a low-fat (control) supplement or supplements with cracked, high-linoleate safflower seeds (linoleate) or cracked, high-oleate safflower seeds (oleate) from d 3 to 90 of lactation. Diets were formulated to be isonitrogenous and isocaloric; safflower seed diets provided 5% of DMI as fat. Plasma and milk samples were collected on d 30, 60, and 90 of lactation. Adipose tissue biopsies were collected near the tail-head region of cows on d 45 and 90 of lactation. In Exp. 2, 3-yr-old cows achieving a BCS of 4 +/- 0.07 (479 +/- 36 kg of BW) or 6 +/- 0.07 (580 +/- 53 kg of BW) at parturition were used in a 2-yr experiment (n = 36/yr). Beginning 3 d postpartum through d 61 of lactation, cows were fed diets similar to those of Exp. 1. Adipose tissue and milk samples were collected on d 30 and 60, and plasma was collected on d 31 and 61 of lactation. Responses to postpartum dietary treatment were comparable in both experiments. Cows fed linoleate and oleate had greater (P < 0.001) total fatty acid concentrations in plasma than cows fed control. Except for 15:1, milk fatty acids with <18 carbons were greatest (P < or = 0.01) for cows fed control, whereas milk from cows fed linoleate had the greatest (P < or = 0.02) 18:1trans-11, 18:2n-6, and cis-9, trans-11 CLA. Milk from cows fed oleate had the greatest (P < 0.001) 18:1cis-9. In Exp. 1, total fatty acid concentrations in adipose tissue samples decreased at d 90 compared with d 45 of lactation, but the fatty acid profile of cow adipose tissue was not affected (P = 0.14 to 0.80) by dietary treatment. In Exp. 2, the percentage of cis-9, trans-11 CLA in adipose tissue of cows with a BCS of 6 decreased (P = 0.001) from d 30 to 60 of lactation. Plasma and milk fatty acid composition reflected alterations in postpartum diet. Less medium-chain fatty acids and more 18-carbon fatty acids in milk were indicative of reduced de novo fatty acid synthesis in the mammary gland of beef cows fed lipid supplements; however, the metabolic demands of lactation prevented the deposition of exogenously derived fatty acids in adipose tissue through d 90 of lactation.