Methylome-wide analysis of milk somatic cells upon subclinical mastitis in dairy cattle

Better understanding of the molecular mechanisms behind bovine mastitis is fundamental for improving the management of this disease, which continues to be of major concern for the dairy industry, especially in its subclinical form. Disease severity and progression depend on numerous aspects, such as livestock genetics, and the interaction between the causative agent, the host, and the environment. In this context, epigenetic mechanisms have proven to have a role in controlling the response of the animal to inflammation. Therefore, in this study we aimed to explore genome-wide DNA methylation of milk somatic cells (SC) in healthy cows (n = 15) and cows affected by naturally occurring subclinical mastitis by Streptococcus agalactiae (n = 12) and Prototheca spp. (n = 11), to better understand the role of SC methylome in the host response to disease. Differentially methylated regions (DMR) were evaluated comparing: (1) Strep. agalactiae-infected versus healthy; (2) Prototheca-infected versus healthy, and (3) mastitis versus healthy and (4) Strep. agalactiae-infected versus Prototheca-infected. The functional analysis was performed at 2 levels. To begin with, we extracted differentially methylated genes (DMG) from promoter DMR, which were analyzed using the Cytoscape ClueGO plug-in. Coupled with this DMG-driven approach, all the genes associated with promoter-methylated regions were fed to the Pathifier algorithm. From the DMR analysis, we identified 1,081 hypermethylated and 361 hypomethylated promoter regions in Strep. agalactiae-infected animals, while 1,514 hypermethylated and 358 hypomethylated promoter regions were identified in Prototheca-infected animals, when compared with the healthy controls. When considering infected animals as a whole group (regardless of the pathogen), we found 1,576 hypermethylated and 460 hypomethylated promoter regions. Both pathogens were associated with methylation differences in genes involved in pathways related to meiosis, reproduction and tissue remodeling. Exploring the whole methylome, in subclinically infected cows we observed a strong deregulation of immune-related pathways, such as nuclear factor kB and toll-like receptors signaling pathways, and of energy-related pathways such as the tricarboxylic acid cycle and unsaturated fatty acid biosynthesis. In conclusion, no evident pathogen-specific SC methylome signature was detected in the present study. Overall, we observed a clear regulation of host immune response driven by DNA methylation upon subclinical mastitis. Further studies on a larger cohort of animals are needed to validate our results and to possibly identify a unique SC methylome that signifies pathogen-specific alterations.

Expression profile of different classes of proteases in milk derived somatic cells across different lactation stages of indigenous cows (Bos indicus) and riverine buffaloes (Bubalus bubalis)

Proteases play a significant role in milk and its products by affecting flavor, texture and longevity. The expression of endogenous proteases varies across different stages of lactation. The study was conducted to understand the transcriptional pattern of different classes of protease-pathways associated genes (CTSB, CTSD, CTSH, CTSL, CTSK, CTSS, CTSZ, PLAU, PLAT) and potential protease inhibitors (SERPIN E2 and SERPIN F2) in 40 milk somatic cells (MSC) samples isolated during early, peak, mid and late lactation stages of Sahiwal cows and Murrah buffaloes - the two most important dairy breeds of India. In Sahiwal cows, except CTSK and PLAU, the expression of other proteases class was not affected significantly (p > 0.05) across lactation stages. However, in Murrah buffaloes, the expression of different proteases increased as the lactation progressed. Most of the proteases showed lower expression during early and peak lactation stages while their expression tends to increase during mid to late lactation stages. The overall trend was somewhat similar in both the dairy species albeit the level of expression was higher in buffalo MSC as compared to cow MSC. The study has provided valuable information on expression kinetics of different proteases in milk somatic cells of two major dairy breeds of India.

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.

Comparative Transcriptome Analysis of Milk Somatic Cells During Lactation Between Two Intensively Reared Dairy Sheep Breeds

In dairy sheep industry, milk production dictates the value of a ewe. Milk production is directly related to the morphology and physiology of the mammary gland; both being designated targets of breeding strategies. Although within a flock breeding parameters are mutual, large differences in milk production among individual ewes are usually observed. In this work, we tested two of the most productive dairy sheep breeds reared intensively in Greece, one local the Chios breed and one foreign the Lacaune breed. We used transcriptome sequencing to reveal molecular mechanisms that render the mammary gland highly productive or not. While highly expressed genes (caseins and major whey protein genes) were common among breeds, differences were observed in differentially expressed genes. ENSOARG00000008077, as a member of ribosomal protein 14 family, together with LPCAT2, CCR3, GPSM2, ZNF131, and ASIP were among the genes significantly differentiating mammary gland's productivity in high yielding ewes. Gene ontology terms were mainly linked to the inherent transcriptional activity of the mammary gland (GO:0005524, GO:0030552, GO:0016740, GO:0004842), lipid transfer activity (GO:0005319) and innate immunity (GO:0002376, GO:0075528, GO:0002520). In addition, clusters of genes affecting zinc and iron trafficking into mitochondria were highlighted for high yielding ewes (GO:0071294, GO:0010043). Our analyses provide insights into the molecular pathways involved in lactation between ewes of different performances. Results revealed management issues that should be addressed by breeders in order to move toward increased milk yields through selection of the desired phenotypes. Our results will also contribute toward the selection of the most resilient and productive ewes, thus, will strengthen the existing breeding systems against a spectrum of environmental threats.

Flow Cytometry-Detected Immunological Markers and on Farm Recorded Parameters in Composite Cow Milk as Related to Udder Health Status

Flow cytometry is a powerful technology used in many fields of cell biology. It is also used as a routine method to count somatic cells in milk and to characterize bovine milk leukocytes. In this study, we used flow cytometry to simultaneously assess the viability, the percentage of the single subsets of leukocytes and to quantify the expression of CD11b, an immunological marker of cell activation status. Immunological markers were then related with on farm recorded parameters as milk electrical conductivity (MEC) and average milk flow rate (MFR). Composite milk samples were collected from 43 cows, nine of which had naturally infected udders and 34 of which had no infected udders. First, the milk samples were classified according to bacteriological test in positive and negative. The results showed that the negative samples to bacteriological test had: (i) significantly higher percentages of live lymphocytes; (ii) significantly lower percentages of CD11b⁺ leukocytes; (iii) significantly lower MEC and higher MFR values. Then, samples were classified in three groups according to somatic cell count (SCC): Group A (n = 15), samples with SCC ≤ 100,000 cells/mL, all negative to bacteriological analysis; Group B (n = 11), samples with 100,000 < SCC < 300,000 cells/m, with four samples positive to bacteriological analysis; Group C (n = 17), samples with SCC ≥ 300,000 cell/mL with five samples positive to bacteriological analysis. Multivariate discriminant analysis was used to verify which flow cytometry immunological markers and on farm recorded parameters could better discriminate among the different groups of SCCs. Linear discriminant analysis (LDA) indicated that 5 of the 10 parameters could best be used to reveal the differences between positive and negative samples among the considered groups of SCCs. Furthermore, the Canonical discriminant analysis (CDA) indicated that composite milk samples with different SCC and infection status were clearly separated from each other in a two-dimensional space. In conclusion, the study highlighted that: (1) the conventional flow cytometry analysis applied on milk samples is a useful tool to investigate immunological parameters as potential indicators of udder health; (2) the combined evaluation of live milk leukocytes and recorded farm parameters could be useful to assess udder health status in dairy cows. The results obtained from this pilot study on few data require new and larger trials to be confirmed.

Effects of Dietary Vegetable Oils on Mammary Lipid-Related Genes in Holstein Dairy Cows

Simple Summary

This study analyzed effects of vegetable oils fed to dairy cows on abundance of genes related to lipid metabolism in milk somatic cells (MSC). During 63 days (9 weeks), 15 cows were allocated to 3 treatments: a control diet with no added lipid and the same diet supplemented with olive oil (OO, 30 g/kg DM) or hydrogenated vegetable oil (HVO, 30 g/kg DM). Dietary oil supplementation (3% DM) had a modest nutrigenomic effect on biological functions such as acetate and FA activation and intra-cellular transport, lipid droplet formation, and transcription regulation in MSC. Results suggest that long-term dietary monounsaturated and saturated lipids could alter mRNA abundance in MSC from mid-lactating cows.

Abstract

This study analyzed effects of vegetable oils fed to dairy cows on abundance of genes related to lipid metabolism in milk somatic cells (MSC). During 63 days, 15 cows were allocated to 3 treatments: a control diet with no added lipid the same diet supplemented with olive oil (OO, 30 g/kg DM) or hydrogenated vegetable oil (HVO, 30 g/kg DM). On days 21, 42 and 63, MSC were obtained from all cows. Relative abundance of genes involved in lipid metabolism in MSC from cows fed control on days 42 and 63 was compared with relative abundance at day 21 to evaluate fold-changes. Those genes without changes over the time were selected to analyze effects of OO and HVO. Compared with control, on day 42, PLIN2 and THRSP were upregulated by OO. Compared with control, on day 21, HVO up regulated ACACA, down regulated FABP3, and on day 63 THRSP and FABP4 were down regulated. Dietary oil supplementation (3% DM) had a modest nutrigenomic effect on different biological functions such as acetate and FA activation and intra-cellular transport, lipid droplet formation, and transcription regulation in MSC.

Gelsolin expression in sheep milk somatic cells during lactation

The identification of genes involved in phenotypes related to milk quality is important for both economic and health aspects in livestock production. The aim of this study was to assess the level of gelsolin gene expression in two breeds of dairy sheep - Sarda and Gentile - with pronounced differences in quantitative and qualitative milk traits. Gelsolin, a type of actin-modulating proteins is involved in the processes of actin remodeling during cell growth and apoptosis; therefore a role of this protein in mammary changes during lactation was here hypothesized. Individual milk samples were collected three times during lactation from 26 ewes of the two breeds. The differential gene expression of gelsolin in the two breeds and the three lactation times was estimated by quantitative PCR on RNA extracted from milk somatic cells. Correlations of gelsolin gene expression with milk yield and quality and days of lactation were also estimated. The results showed that gelsolin gene expression was significantly higher in the Sarda compared to the Gentile at each lactation stage, in agreement with the longer lactation duration and the higher daily milk yield of the first breed. Significant correlations of gelsolin gene expression were found with milk fat content in Sarda breed (-0.46, P<0.05). Gelsolin expression analysis confirmed the link between gelsolin gene function and milk fat content of sheep.

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.

New Rapid Method of DNA Isolation from Milk Somatic Cells

Isolation of genomic DNA is one of the basic steps in many different molecular analyses. There are a few reports on methods of DNA isolation from milk, but many of them are time consuming and expensive, and require relatively large volumes of raw milk. In this study a rapid, sensitive, and efficient method of DNA extraction from milk somatic cells of various mammals (cattle, sheep, goats, horses) is presented. It was found that milk is a good source of genomic DNA, and to obtain a sufficient amount and quality of DNA, suitable for molecular analysis such as PCR, 10 mL of raw milk is sufficient. Thanks to this method, stress in animals can be reduced during collection of researched material. Therefore, this method could be widely used in molecular analyses.

Evaluation of internal reference genes for quantitative expression analysis by real-time reverse transcription-PCR in somatic cells from goat milk

Reverse transcription (RT) quantitative real-time PCR (qPCR) is the most accurate and easy-to-perform technique to measure the expression level of a selected gene of interest by quantifying mRNA transcripts. The use of reference genes is commonly accepted as the most reliable approach to normalize RT-qPCR data and reduce possible errors generated in the quantification of gene expression. The optimal number and choice of reference genes are experimentally validated for specific tissues or cell types and experimental designs. To date, data on qPCR normalization in goats are scarce and the most suitable reference genes in this species have been identified for only a limited number of tissues. The aim of this study was to determine an optimal combination of stably expressed reference genes in caprine milk somatic cells (MSC) from healthy and infected mammary glands. For the purpose, we performed RT-qPCR for 10 commonly used reference genes from various functional classes and then determined their expression level in MSC from goats intramammary challenged with Staphylococcus aureus and in MSC from healthy controls, with a view to select genes whose stability would be unaffected under infection conditions. The geNorm and NormFinder algorithms were used for validating the reference genes. Furthermore, to demonstrate the importance of normalization of gene expression with appropriate reference genes, we tested the effect of using a combination of the least stable genes for expression analysis evaluation. On the basis of our evaluation, we recommend the use of a panel of reference genes that should include G6PD, YWHAZ, and ACTB for caprine MSC gene expression profiling. The expression of the 2 genes of interest, pentraxin-related protein (PTX3) and secreted phosphoprotein 1 (SPP1), was evaluated by RT-qPCR in all samples collected pre- and postinfection, and the recommended reference genes were used to normalize the data. Our study provides a validated panel of optimal reference genes for the identification of genes differentially expressed by qRT-PCR in caprine MSC. Moreover, we provided a set of intron-spanning primer sequences that could be suitable for gene expression experiments using SYBR Green chemistry on other caprine tissues and cells.

Caloric stress alters fat characteristics and Hsp70 expression in milk somatic cells of lactating beef cows

Selection for higher production rate in cattle inhabiting challenging habitats may be considered disadvantageous because of possible deleterious effects on immunity and reproduction and, consequently, on calf crop percentage. In Israel, free-grazing high productive beef cows experience reduction in nutritional quality of forage during up to 8 months of the year. As milk production by dams dictates calf performance, dam's nutritional needs and rebreeding rates, the aim of the present study was to test how lactating beef cows deal with combined caloric and protein stress both at the productive and self protective levels. For this purpose, we studied the effect of long-term caloric stress on milk characteristics and gene expression of stress and milk components producing proteins. Lactating dams responded to caloric stress by decreased body weight, milk, and milk protein production. To compensate for total energy loses in milk, they produced milk of higher fat concentration and shifted the proportions of its fatty acids towards long and unsaturated ones. This was reflected by increased mRNA transcription of the fatty acid binding protein. Prolonged low-energy diet promoted cell-specific heat shock protein (Hsp) response; whereas significant increase of Hsp90 but unchanged levels of Hsp70 proteins were observed in white blood cells, the expression of Hsp70 in milk somatic cells was markedly attenuated, in parallel with a marked increase of alpha(s1)-casein expression. At the mammary gland level, these results may indicate a decrease in turnover of proteins and a shift to an exclusive expression of milk components producing factors. Similar responses to caloric stress were revealed also in ketotic dairy cows. Ketosis promoted a shift towards long and unsaturated fatty acids and an increased expression of alpha(s1)-casein in milk somatic cells. These findings may reflect an evolutionary-preserved mechanism in lactating cows for coping with caloric restriction. Overall, our results provide an index to test suitability of beef cattle breeds to inadequate caloric demands.