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

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.

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.

Transcriptome profiling of Staphylococci-infected cow mammary gland parenchyma

BACKGROUND

Genome-wide gene expression profiling allows for identification of genes involved in the defense response of the host against pathogens. As presented here, transcriptomic analysis and bioinformatics tools were applied in order to identify genes expressed in the mammary gland parenchyma of cows naturally infected with coagulase-positive and coagulase-negative Staphylococci.

RESULTS

In cows infected with coagulase-positive Staphylococci, being in 1st or 2nd lactation, 1700 differentially expressed genes (DEGs) were identified. However, examination of the 3rd or 4th lactations revealed 2200 DEGs. Gene ontology functional classification showed the molecular functions of the DEGs overrepresented the activity of cytokines, chemokines, and their receptors. In cows infected with coagulase-negative Staphylococci, in the 1st or 2nd lactations 418 DEGs, while in the 3rd or 4th lactations, 1200 DEGs were identified that involved in molecular functions such as protein, calcium ion and lipid binding, chemokine activity, and protein homodimerization. Gene network analysis showed DEGs associated with inflammation, cell migration, and immune response to infection, development of cells and tissues, and humoral responses to infections caused by both types of Staphylococci.

CONCLUSION

A coagulase-positive Staphylococci infection caused a markedly stronger host response than that of coagulase-negative, resulting in vastly increased DEGs. A significant increase in the expression of the FOS, TNF, and genes encoding the major histocompatibility complex proteins (MHC) was observed. It suggests these genes play a key role in the synchronization of the immune response of the cow's parenchyma against mastitis-causing bacteria. Moreover, the following genes that belong to several physiological pathways (KEGG pathways) were selected for further studies as candidate genes of mammary gland immune response for use in Marker Assisted Selection (MAS): chemokine signaling pathway (CCL2, CXCL5, HCK, CCR1), cell adhesion molecules (CAMs) pathway (BOLA-DQA2, BOLA-DQA1, F11R, ITGAL, CD86), antigen processing and presentation pathway (CD8A, PDIA3, LGMN, IFI30, HSPA1A), and NOD-like receptor signaling pathway (TNF, IL8, IL18, NFKBIA).

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.

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).

Gene expression analysis of protein synthesis pathways in bovine mammary epithelial cells purified from milk during lactation and short-term restricted feeding

The objective of the study was to investigate selected key regulatory pathways of milk protein biosynthesis in primary bovine mammary epithelial cells (MECs) of dairy cows during the first 155 days of lactation. In addition, cows were exposed to feed restriction for a short period (FR) during different stages of lactation (week 4 and 21 pp) to study adjustment processes of molecular protein biosynthesis to metabolic challenge. Morning milk samples from twenty-four Holstein-Friesian cows were collected throughout the experimental period (n = 10 per animal). MEC from raw milk were purified using an immunomagnetic separation technique and used for real-time quantitative PCR analyses. As was seen in transcript abundances of all major milk proteins, mRNA levels of E74-like factor 5 (ELF5), an enhancer of signal transducer and activator of transcription (STAT) action, concomitantly decreased towards mid-lactation. Expression of ELF5 as well as of all milk protein genes showed a similar increase during FR in early lactation. Occasional changes in expression could be seen in other Janus kinase (JAK)/STAT factors and in mammalian target of rapamycin (mTOR) pathway elements. Amino acid transfer and glucose transporter and the β-casein expression were also partially affected. In conclusion, our findings suggest a pivotal role of the transcription factor ELF5 in milk protein mRNA expression with complementary JAK/STAT and mTOR signalling for the regulation of protein biosynthesis in the bovine mammary gland.

Validation of RNA isolated from milk fat globules to profile mammary epithelial cell expression during lactation and transcriptional response to a bacterial infection

Mastitis, an inflammation of the mammary gland, is the most costly infectious disease of dairy ruminants worldwide. Although it receives considerable attention, the early steps of the host response remain poorly defined. Here, we report a noninvasive method using milk fat globules (MFG) as a source of mammary RNA to follow the dynamics of the global transcriptional response of mammary epithelial cells (MEC) during the course of a bacterial infection. We first assessed that RNA isolated from MFG were representative of MEC RNA; we then evaluated whether MFG RNA could be used to monitor the MEC response to infection. Sufficiently high yields of good-quality RNA (RNA integrity numbers ranging between 6.7 and 8.7) were obtained from goat MFG for subsequent analyses. Contamination of MFG by macrophages and neutrophils, which can be trapped during creaming, was assessed and when using quantitative real-time PCR for cell-type specific markers, was shown to be weak enough (<8%) to affect MFG gene expression profiling. Using microarrays, we showed that RNA extracted from MFG and from mammary alveolar parenchyma shared approximately 90% of the highlighted probes corresponding in particular to genes encoding milk proteins (CSN, BLG, LALBA) and enzymes involved in milk fat synthesis and secretion (FASN, XDH, ADRP, SCD, and DGAT1). In addition, a gene involved in the acute-phase reaction, coding for the serum amyloid A3 (SAA3) protein, was found within the first 50 most highly expressed genes in a noninfectious context in both mammary alveolar parenchyma and MFG, strongly suggesting that SAA3 is expressed in MEC. We took advantage of this noninvasive RNA sampling to follow the early proinflammatory response of MEC during the course of a bacterial infection and showed that the levels of mRNA encoding SAA3 sharply increased at 24h postinfection. Taken together, our results demonstrate that MFG represent a unique source of MEC RNA to noninvasively sample sufficient amounts of high-quality RNA to assess the dynamics of MEC gene expression in vivo, especially during the first steps of infection, thereby paving the way for the discovery of early biomarkers for the control of intramammary infections. Furthermore, this noninvasive technique could be used to provide mammary transcriptomic data on a large scale, thus filling the gap between genomic and phenotypic data.