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

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.

Unraveling the Genomic Association for Milk Production Traits and Signatures of Selection of Cattle in a Harsh Tropical Environment

A study was designed to identify the genomic regions associated with milk production traits in a dairy cattle population reared by smallholder farmers in the harsh and challenging tropical savanna climate of Bengaluru, India. This study is a first-of-its-kind attempt to identify the selection sweeps for the dairy cattle breeds reared in such an environment. Two hundred forty lactating dairy cows reared by 68 farmers across the rural-urban transiting regions of Bengaluru were selected for this study. A genome-wide association study (GWAS) was performed to identify candidate genes for test-day milk yield, solids-not-fat (SNF), milk lactose, milk density and clinical mastitis. Furthermore, the cross-population extended haplotype homozygosity (XP-EHH) methodology was adopted to scan the dairy cattle breeds (Holstein Friesian, Jersey and Crossbred) in Bengaluru. Two SNPs, rs109340659 and rs41571523, were observed to be significantly associated with test-day milk yield. No significant SNPs were observed for the remaining production traits. The GWAS for milk lactose revealed one SNP (rs41634101) that was very close to the threshold limit, though not significant. The potential candidate genes fibrosin-like 1 (FBRSL) and calcium voltage-gated channel auxiliary subunit gamma 3 (CACN) were identified to be in close proximity to the SNP identified for test-day milk yield. These genes were observed to be associated with milk production traits based on previous reports. Furthermore, the selection signature analysis revealed a number of regions under selection for the breed-group comparisons (Crossbred-HF, Crossbred-J and HF-J). Functional analysis of these annotated genes under selection indicated pathways and mechanisms involving ubiquitination, cell signaling and immune response. These findings point towards the probable selection of dairy cows in Bengaluru for thermotolerance.

Identifying Suitable Reference Gene Candidates for Quantification of DNA Damage-Induced Cellular Responses in Human U2OS Cell Culture System

DNA repair pathways trigger robust downstream responses, making it challenging to select suitable reference genes for comparative studies. In this study, our goal was to identify the most suitable housekeeping genes to perform comparable molecular analyses for DNA damage-related studies. Choosing the most applicable reference genes is important in any kind of target gene expression-related quantitative study, since using the housekeeping genes improperly may result in false data interpretation and inaccurate conclusions. We evaluated the expressional changes of eight well-known housekeeping genes (i.e., 18S rRNA, B2M, eEF1α1, GAPDH, GUSB, HPRT1, PPIA, and TBP) following treatment with the DNA-damaging agents that are most frequently used: ultraviolet B (UVB) non-ionizing irradiation, neocarzinostatin (NCS), and actinomycin D (ActD). To reveal the significant changes in the expression of each gene and to determine which appear to be the most acceptable ones for normalization of real-time quantitative polymerase chain reaction (RT-qPCR) data, comparative and statistical algorithms (such as absolute quantification, Wilcoxon Rank Sum Test, and independent samples T-test) were conducted. Our findings clearly demonstrate that the genes commonly employed as reference candidates exhibit substantial expression variability, and therefore, careful consideration must be taken when designing the experimental setup for an accurate and reproducible normalization of RT-qPCR data. We used the U2OS cell line since it is generally accepted and used in the field of DNA repair to study DNA damage-induced cellular responses. Based on our current data in U2OS cells, we suggest using 18S rRNA, eEF1α1, GAPDH, GUSB, and HPRT1 genes for UVB-induced DNA damage-related studies. B2M, HPRT1, and TBP genes are recommended for NCS treatment, while 18S rRNA, B2M, and PPIA genes can be used as suitable internal controls in RT-qPCR experiments for ActD treatment. In summary, this is the first systematic study using a U2OS cell culture system that offers convincing evidence for housekeeping gene selection following treatment with various DNA-damaging agents. Here, we unravel an indispensable issue for performing and assessing trustworthy DNA damage-related differential gene expressional analyses, and we create a "zero set" of potential reference gene candidates.

Intracellular Staphylococcus aureus Infection Decreases Milk Protein Synthesis by Preventing Amino Acid Uptake in Bovine Mammary Epithelial Cells

Staphylococcus aureus (S. aureus) is one of the main pathogens in cow mastitis, colonizing mammary tissues and being internalized into mammary epithelial cells, causing intracellular infection in the udder. Milk that is produced by cows that suffer from mastitis due to S. aureus is associated with decreased production and changes in protein composition. However, there is limited information on how mastitis-inducing bacteria affect raw milk, particularly with regard to protein content and protein composition. The main purpose of this work was to examine how S. aureus infection affects milk protein synthesis in bovine mammary epithelial cells (BMECs). BMECs were infected with S. aureus, and milk protein and amino acid levels were determined by ELISA after S. aureus invasion. The activity of mTORC1 signaling and the transcription factors NF-κB and STAT5 and the expression of the amino acid transporters SLC1A3 and SLC7A5 were measured by western blot or immunofluorescence and RT-qPCR. S. aureus was internalized by BMECs in vitro, and the internalized bacteria underwent intracellular proliferation. Eight hours after S. aureus invasion, milk proteins were downregulated, and the level of BMECs that absorbed Glu, Asp, and Leu from the culture medium and the exogenous amino acids induced β-casein synthesis declined. Further, the activity of mTORC1 signaling, NF-κB, and STAT5 was impaired, and SLC1A3 and SLC7A5 were downregulated. Eight hours of treatment with 100 nM rapamycin inhibited NF-κB and STAT5 activity, SLC1A3 and SLC7A5 expression, and milk protein synthesis in BMECs. Thus mTORC1 regulates the expression of SLC1A3 and SLC7A5 through NF-κB and STAT5. These findings constitute a model by which S. aureus infection suppresses milk protein synthesis by decreasing amino acids uptake in BMECs.

Proteomics study on the protective mechanism of soybean isoflavone against inflammation injury of bovine mammary epithelial cells induced by Streptococcus agalactiae

This study aimed to verify the anti-inflammatory effect of soybean isoflavones (SI) on the inflammatory response induced by Streptococcus agalactiae (S. agalactiae) of bovine mammary epithelial cells (bMECs) and to elucidate its possible mechanism. BMECs were pretreated with SI of different concentrations (20, 40, 60, 80, 100 μg/mL) for 0.5, 3, 6, 9, 12, 15, 18, 24 h. And then, S. agalactiae was used to infect bMECs for 6 h (MOI = 50:1) to establish the inflammation model. Cell viability, growth curves of S. agalactiae, cytotoxicity, and S. agalactiae invasion rate were determined. A proteomics technique was used to further detect differential proteins and enrichment pathways. SI (40 μg/mL) improved the viability of bMECs at 12 h (p < 0.05) and 60 and 80 μg/mL of SI greater (p < 0.01). Moreover, 60 μg/mL of SI protects cells from bacterial damage (p < 0.05). SI could inhibit S. agalactiae growth and internalization into bMECs in a time- and dose-dependent manner. In addition, proteomics results showed that 133 proteins were up-regulated and 89 proteins were down-regulated significantly. The differentially significantly expressed proteins (DSEPs) were mainly related to cell proliferation, differentiation, apoptosis, and migration. GO annotation showed that 222 DSEPs were divided into 23 biological processes (BP) terms, 14 cell components (CC) terms, and 12 molecular functions (MF) terms. DSEPs were significantly enriched in 10 pathways, of which the immune pathway was the main enrichment pathway.

Feeding level regulates the expression of some genes involved with programed cell death and remodeling in goat and sheep mammary tissue

Mammary tissue (MT) turnover is characterized by programed cell death and remodeling which might be affected by both feeding level and animal species. Thus, twenty-four dairy goats and the same number of sheep were assigned to three homogenous sub-groups per animal species and fed the same diet in quantities which met 70% (FL70), 100% (FL100) and 130% (FL130) of their daily energy and crude protein requirements. Individual MT samples were taken by biopsy from the animals on the 30th and 60th experimental day. The results showed, in the first sampling time, a significant reduction in the mRNA abundance for selected genes involved in programed cell death in both FL 70 fed goats (STAT3 and BECN1) and sheep (CASPASE8 and BECN1) compared with the respective FL100 groups. The FL130, in comparison with the FL100, caused a significant increase in transcripts accumulation of STAT3 gene in both sampling times and CASPASE8 gene in the second sampling time in goat MT, while the opposite happened for the mRNA expression of CASPASE8 and BECN1 genes in sheep MT, but only in the first sampling time. Moreover, a significant up regulation in the mRNA levels of MMP2 gene in MT of FL130 fed sheep was observed. The FL130, in comparison with the FL70, caused an enhancement in the mRNA expression levels of BECN1, CASPASE8, BAX and STAT3 genes in goat MT only. It was also shown that apoptosis and autophagy can be affected simultaneously by the feeding level. Overfeeding affects MT programed cell death and remodeling by a completely different way in goats than sheep. In conclusion, feeding level and animal species have strong effects on both MT programed cell death (apoptosis and autophagy) and remodeling but the molecular mechanisms need further investigation.

Role of the JAK-STAT Pathway in Bovine Mastitis and Milk Production

Simple Summary

The cytokine-activated Janus kinase (JAK)—signal transducer and activator of transcription (STAT) pathway has an important role in the regulation of immunity and inflammation. In addition, the signaling of this pathway has been reported to be associated with mammary gland development and milk production. Because of such important functions, the JAK-STAT pathway has been widely targeted in both human and animal diseases as a therapeutic agent. Recently, the JAK2, STATs, and inhibitors of the JAK-STAT pathway, especially cytokine signaling suppressors (SOCSs), have been reported to be associated with milk production and mastitis-resistance phenotypic traits in dairy cattle. Thus, in the current review, we attempt to overview the development of the JAK-STAT pathway role in bovine mastitis and milk production.

Abstract

The cytokine-activated Janus kinase (JAK)—signal transducer and activator of transcription (STAT) pathway is a sequence of communications between proteins in a cell, and it is associated with various processes such as cell division, apoptosis, mammary gland development, lactation, anti-inflammation, and immunity. The pathway is involved in transferring information from receptors on the cell surface to the cell nucleus, resulting in the regulation of genes through transcription. The Janus kinase 2 (JAK2), signal transducer and activator of transcription A and B (STAT5 A & B), STAT1, and cytokine signaling suppressor 3 (SOCS3) are the key members of the JAK-STAT pathway. Interestingly, prolactin (Prl) also uses the JAK-STAT pathway to regulate milk production traits in dairy cattle. The activation of JAK2 and STATs genes has a critical role in milk production and mastitis resistance. The upregulation of SOCS3 in bovine mammary epithelial cells inhibits the activation of JAK2 and STATs genes, which promotes mastitis development and reduces the lactational performance of dairy cattle. In the current review, we highlight the recent development in the knowledge of JAK-STAT, which will enhance our ability to devise therapeutic strategies for bovine mastitis control. Furthermore, the review also explores the role of the JAK-STAT pathway in the regulation of milk production in dairy cattle.

Single-step genome-wide association studies (GWAS) and post-GWAS analyses to identify genomic regions and candidate genes for milk yield in Brazilian Girolando cattle

Milk production is economically important to the Brazilian agribusiness, and the majority of the country's milk production derives from Girolando (Gir × Holstein) cows. This study aimed to identify quantitative trait loci (QTL) and candidate genes associated with 305-d milk yield (305MY) in Girolando cattle. In addition, we investigated the SNP-specific variances for Holstein and Gir breeds of origin within the sequence of candidate genes. A single-step genomic BLUP procedure was used to identify QTL associated with 305MY, and the most likely candidate genes were identified through follow-up analyses. Genomic breeding values specific for Holstein and Gir were estimated in the Girolando animals using a model that uses breed-specific partial relationship matrices, which were converted to breed of origin SNP effects. Differences between breed of origin were evaluated by comparing estimated SNP variances between breeds. From 10 genome regions explaining most additive genetic variance for 305MY in Girolando cattle, 7 candidate genes were identified on chromosomes 1, 4, 6, and 26. Within the sequence of these 7 candidate genes, Gir breed of origin SNP alleles showed the highest genetic variance. These results indicated QTL regions that could be further explored in genomic selection panels and which may also help in understanding the gene mechanisms involved in milk production in the Girolando breed.

Polymorphisms in JAK2 Gene are Associated with Production Traits and Mastitis Resistance in Dairy Cattle

The present study was designed to investigate the effects of single nucleotide polymorphisms (SNPs) in the JAK2 gene on the production and mastitis related traits in dairy cattle. Blood and milk samples were collected from 201 lactating dairy cattle of three breeds, i.e. Holstein Friesian (HF), Jersey (J) and Achai (A) and their crosses maintained at well-established dairy farms in Khyber Pakhtunkhwa, Pakistan. Generalized linear model was used to evaluate the association between genotypes and the studied traits. A DNA pool was made from randomly selected 30 samples which revealed three SNPs, i.e. SNP 1 in 5’ upstream region (G>A, rs379754157), SNP 2 in intron 15 (A>G, rs134192265), and SNP 3 in exon 20 (A>G, rs110298451) that were further validated in the population under study using SNaPshot technique. Of the three SNPs, SNP 1 did not obey Hardy-Weinberg equilibrium (P<0.05). SNP 2 and SNP 3 were found to be in strong linkage disequilibrium and allele G was highly prevalent compared to allele A in these SNPs. in SNP 1, the GG genotype was associated with significantly (P<0.01) higher SCC, whereas SNP 2 and SNP 3 were significantly (P<0.01) associated with higher lactose percentage compared to the other geno-types. The haplogroups association analysis revealed that H1H2 (GG GG AG) has significantly lower SCC than H2H2 (GG GG GG). The results infer that JAK2 could be an important candidate gene and the studied SNPs might be useful genetic markers for production and mastitis related traits.

Review: the cellular mechanisms underlying mammary tissue plasticity during lactation in ruminants

The mammary tissue is characterized by its capacity to adapt in response to a wide variety of changing conditions. This adaptation capacity is referred to as the plasticity of mammary tissue. In dairy ruminants, lactation is challenged by modifications that can either be induced on purpose, such as by modifying management practices, or occur involuntarily, when adverse environmental constraints arise. These modifications can elicit both immediate changes in milk yield and composition and carryover effects that persist after the end of the challenge. This review focuses on the current knowledge concerning the cellular mechanisms underlying mammary tissue plasticity. The main mechanisms contributing to this phenomenon are changes in the activity and number of mammary epithelial cells (MECs). Changes in the number of these cells result from variations in the rates of cell proliferation and death as well as changes in the rate MEC exfoliation. The number of MECs also depends on the number of resident adult mammary stem cells and their progenitors, which can regenerate the pools of the various mammary cells. Several challenges, including changes in milking frequency, changes in level of feed supply and hormonal manipulations, have been shown to modulate milk yield together with changes in mammary cell activity, turnover and exfoliation. Epigenetic changes may be an additional mechanism of adaptation. Indeed, changes in DNA methylation and reductions in milk yield have been observed during once-daily milking and during mastitis in dairy cows and may affect cell activity persistently. In contrast to what has been assumed for a long time, no carryover effect on milk yield were observed after feed supply challenges in dairy cows and modification of milking frequency in dairy goats, even though the number of mammary cells was affected. In addition, mammary tissue plasticity has been shown to be influenced by the stage of lactation, health status and genetic factors. In conclusion, the cellular mechanisms underlying mammary tissue plasticity are diverse, and the mammary tissue either does or does not show elastic properties (with no permanent deformation), in response to environmental changes.

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.

Short communication: Characterization of gene expression profiles related to yak milk protein synthesis during the lactation cycle

This research assessed the gene expression patterns related to the synthesis of milk in yak, which is characterized by high fat and protein content but low yield. The yak (Bos grunniens) is one of the most crucial domestic animals in Tibetan life; however, the genetic and molecular factors underlying yak milk protein synthesis remain understudied. Yak mammary biopsies harvested during late-pregnancy (d -15) through the end of subsequent lactation (d 1, 15, 30, 60, 180, and 240) were used to evaluate gene expression via real-time quantitative PCR. The expression pattern of 41 genes encompassing multiple pathways integral to milk protein synthesis including insulin, mammalian target of rapamycin (mTOR), 5' AMP-activated protein kinase, Jak2-Stat5 signaling, and the expression of glucose and AA transporters was evaluated. Our results confirmed that most upregulated genes increased from d -15 and peaked at d 30 or 60 and then remained relatively highly expressed. Specifically, there was an increased expression of mTOR-related amino acid transporters (SLC1A5, SLC7A5, and SLC36A1), glucose transporters (SLC2A1, SLC2A3, and SLC2A8), Jak2-Stat5 pathway (ELF5), and insulin signaling pathway components (IRS1, PDPK1, and AKT1). For activation of proteins synthesis, MTOR was significantly increased only at d 1. Among inhibitors of mTOR signaling, TSC1 and PRKAA2 were significantly upregulated during lactation. The RPL23 was downregulated among ribosomal components. In conclusion, a critical role for AA and glucose transporters and insulin signaling through mTOR for regulation of yak milk protein synthesis was revealed in this study of the yak mammary gland.

Relative gene expression analysis of β-casein gene and its transcription regulatory genes in primary buffalo mammary epithelial cells of Surti and Jaffarabadi buffaloes

The present study was designed to estimate the relative gene expression of major genes responsible for β-casein and its transcription regulation during early lactation cycle in Surti and Jaffarabadi buffaloes. Ten buffaloes of each breed maintained at Livestock Research Station, Navsari and Cattle Breeding Farm, Junagadh, respectively were selected for the study. The RNA was extracted from noninvasively isolated primary buffalo mammary epithelial cells (pBMEC) obtained using antibody mediated magnetic separation method from milk samples collected at day 15 and 60 postpartum (pp). Primers used for amplification of β-casein gene (CSN2) and its transcription regulatory genes were bovine specific. The mean relative expressions of C/EBPβ, RUNX2 and STAT5A genes were significantly higher in Jaffarabadi buffaloes as compared to Surti buffaloes at day 15 pp. However, between breed differences in the mean relative expressions of C/EBPβ, RUNX2 and STAT5A genes were not observed at 60 days pp. The relative expressions of CSN2 and YY1 genes did not differ significantly within or between breeds. The mean relative expression of KRT8 gene was almost comparable among all the groups under present study. Thus, only β- casein transcription regulatory genes C/EBPβ, RUNX2 and STAT5A amongst all studied genes showed significant between breed variability at day 15 pp, however this variability too got subsided as lactation stage advanced to day 60 pp.

Establishment of a 3D cell culture model of primary bovine mammary epithelial cells extracted from fresh milk

For the investigation of molecular processes underlying diseases of the bovine mammary gland, primary bovine mammary epithelial cells (pbMEC) are used. They are known to contribute to the innate immune system of the bovine mammary gland. The functionality of pbMEC depends on the maintenance of in vivo characteristics. So far, the optimization of pbMEC culture conditions was intended in a variety of experiments. For this purpose, most of the studies used stable cell lines or primary cells obtained from udder biopsies of slaughtered animals. By contrast, within our study, pbMEC of healthy and first lactating Brown Swiss cows were non-invasively isolated from fresh milk. The non-invasively isolated pbMEC were cultivated on the extracellular matrix-like scaffold Matrigel®_. Further, they were challenged with different compositions of proliferation media, containing lactogenic hormones and/or the essential amino acid L-lysine. Changes in expression levels of genes coding for milk proteins and for components of the janus kinase/signal transducers and activators of transcription (JAK-STAT) and mTOR pathways were analyzed by RT-qPCR. The secreted proteins were analyzed by LC-MS/MS measurements. We showed for the first time the establishment of a physiologically functional 3D cell culture model of pbMEC isolated from fresh milk. This represents a primary cell culture model system, based on non-invasive cell collection, that can be used to unravel physiological processes in an unbiased manner.

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

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

The effect of long term under- and overfeeding on the expression of genes related to glucose metabolism in the mammary tissue of goats

The glucose metabolism in the mammary tissue (MT) is controlled by several genes whose nutritional regulation is poorly documented. In this study, the effect of long-term under- and overfeeding on the expression of glucose transporter 1 (GLUT1), glucose transporter 3 (GLUT3), sodium glucose contransporter 1 (SGLT1), two isoforms of β- (1, 4) galactosyltransferase [β- (1, 4) GAT1 and β- (1, 4) GAT3], and α-lactalbumin (LALBA), related to glucose metabolism in goat MT, was examined by using real-time RT-PCR. Twenty- four lactating goats were divided into three homogenous subgroups and fed the same ration in quantities covering 70% (underfeeding), 100% (control) and 130% (overfeeding) of their energy and crude protein requirements, respectively. The results showed that the feeding level did not affect the GLUT1, GLUT3 and SGLT1 genes expression in goats MT. On the contrary, a significant reduction and a trend for reduction on mRNA of β- (1, 4) GAT1 and on β- (1, 4) GAT3, respectively, in the MT of underfed goats, compared with the overfed ones, was observed. Moreover, a significant decrease in the LALBA mRNA accumulation in the MT of underfed goats compared with the overfed was found. In conclusion, the MT of goats, unlike cows, adapt to changes in glucose or energy supply from different levels of feeding by changing the utilisation of glucose for the synthesis of lactose.

Differences during the first lactation between cows cloned by somatic cell nuclear transfer and noncloned cows

Lactation performance is dependent on both the genetic characteristics and the environmental conditions surrounding lactating cows. However, individual variations can still be observed within a given breed under similar environmental conditions. The role of the environment between birth and lactation could be better appreciated in cloned cows, which are presumed to be genetically identical, but differences in lactation performance between cloned and noncloned cows first need to be clearly evaluated. Conflicting results have been described in the literature, so our aim was to clarify this situation. Nine cloned Prim' Holstein cows were produced by the transfer of nuclei from a single fibroblast cell line after cell fusion with enucleated oocytes. The cloned cows and 9 noncloned counterparts were raised under similar conditions. Milk production and composition were recorded monthly from calving until 200d in milk. At 67d in milk, biopsies were sampled from the rear quarter of the udder, their mammary epithelial cell content was evaluated, and mammary cell renewal, RNA, and DNA were then analyzed in relevant samples. The results showed that milk production did not differ significantly between cloned and noncloned cows, but milk protein and fat contents were less variable in cloned cows. Furthermore, milk fat yield and contents were lower in cloned cows during early lactation. At around 67 DIM, milk fat and protein yields, as well as milk fat, protein, and lactose contents, were also lower in cloned cows. These lower yields could be linked to the higher apoptotic rate observed in cloned cows. Apoptosis is triggered by insulin-like factor growth binding protein 5 (IGFBP5) and plasminogen activator inhibitor (PAI), which both interact with CSN1S2. During our experiments, CSN1S2 transcript levels were lower in the mammary gland of cloned cows. The mammary cell apoptotic rate observed in cloned cows may have been related to the higher levels of DNA (cytosine-5-)-methyltransferase 1 (DNMT1) transcripts, coding for products that maintain the epigenetic status of cells. We conclude, therefore, that milk production in cloned cows differs slightly from that of noncloned cows. These differences may be due, in part, to a higher incidence of subclinical mastitis. They were associated with differences in cell apoptosis and linked to variations in DNMT1 mRNA. However, milk protein and fat contents were more similar among cloned cows than among noncloned cows.

The effect of long-term under- and overfeeding on the expression of six major milk proteins' genes in the mammary tissue of goats

Milk protein synthesis in the mammary gland involves expression of six major milk proteins' genes whose nutritional regulation remains poorly defined. In this study, the effect of long-term under- and overfeeding on the expression of as1-casein: CSN1S1, as2-casein: CSN1S2, β-casein: CSN2, κ-casein: CSN3, α-lactalbumin: LALBA and β-lactoglobulin: BLG gene in goat mammary tissue (MT) was examined. Twenty-four lactating dairy goat, at 90-98 days in milk, were divided into three homogenous subgroups and fed the same ration, for 60 days, in quantities which met 70% (underfeeding), 100% (control) and 130% (overfeeding) of their energy and crude protein requirements. The results showed a significant decrease in mRNA of CSN1S2, CSN2, CSN3 and LALBA genes in the MT of underfed goats compared with the overfed and on the CSN1S1 and BLG gene expressions in the MT of underfed goats compared with the respective control and overfed. CSN2 was the most abundant transcript in goat MT relative to the other milk proteins' genes. Significantly positive correlations were observed between the mRNA levels of caseins' and BLG genes with the milk yield. Moreover, a significant correlation was found between the mRNA levels of CSN1S2 with the milk protein, lactose content and lactose yield and also between the LALBA gene expression with the lactose content and lactose yield respectively. In conclusion, the feeding level and consequently the nutrients availability affected the milk lactose content, protein and lactose yield as well as the milk volume by altering the CSN1S1, CSN1S2, CSN2, CSN3, LALBA and BLG gene expression involved in their metabolic pathways.

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

Cortisol levels in skimmed milk during the first 22 weeks of lactation and response to short-term metabolic stress and lameness in dairy cows

Background

Cortisol is secreted into blood in reaction to acute stress, but also in phases of diminished feed intake and changed animal behavior. As cows do not always show clear signs of discomfort, reliable diagnostic markers could be used to provide information regarding individual cows’ distress. The objective of this study was to establish an ether free immunoassay for the detection of cortisol and to determine values during the first 22 weeks of lactation. Furthermore, the response in milk cortisol levels was assessed during times of metabolic stress and pain associated symptoms of lameness.

Methods

Milk yield and composition, blood serum glucose, NEFA and BHBA as well as milk cortisol were determined in 24 multiparous Holstein-Friesian cows over the course of the first 22 weeks of lactation. Animals were further checked for signs of clinical diseases on a daily basis. Two feed restrictions over three days (FR; 70 % of precious ad libitum intake) were performed during the 4^th wk and the 21^st wk, respectively. An ELISA for cortisol measurement in easily accessible bovine skimmed milk was established and applied.

Results

On the last day of FR in early lactation, a reduction in milk yield and changes in serum metabolites compared to respective previous values were detected. The FR in mid-lactation resulted in no changes in milk production and serum metabolites. Milk cortisol was highest during first wk of lactation and remained on comparable levels thereafter. Milk yield and composition were not influenced by FR. Lameness resulted in enhanced milk cortisol levels.

Conclusion

Milk cortisol could be used as an indicator of painful symptoms such as lameness. Higher values of milk cortisol levels during first wk of lactation should be taken into account for interpretation.

The effect of long term under- and over-feeding on the expression of six major milk protein genes in the mammary tissue of sheep

Milk protein synthesis in the mammary gland involves expression of six major milk protein genes whose nutritional regulation remains poorly defined. In this study, the effect of long term under- and over-feeding on the expression of αs1-casein: CSN1S1, αs2-casein: CSN1S2, β-casein: CSN2, κ-casein: CSN3, α-lactalbumin: LALBA and β-lactoglobulin: BLG gene in sheep mammary tissue (MT) was examined. Twenty-four lactating dairy sheep, at 90-98 d in milk, were divided into three groups and fed the same ration, for 60 d, in quantities which met 70% (underfeeding), 100% (control) and 130% (overfeeding) of their energy and crude protein requirements. The results showed a significant reduction on mRNA of CSN1S1, CSN1S2, CSN2 and BLG gene in the MT of underfed sheep compared with the overfed ones and a significant reduction in CSN3 and LALBA gene expression compared with the respective control animals. Significant positive correlations were observed between the mRNA levels of milk proteins' genes with the milk protein yield and milk yield respectively. In conclusion, the feeding level and consequently the nutrients availability, affected the milk protein yield and milk volume by altering the CSN1S1, CSN1S2, CSN2, CSN3, LALBA and BLG gene expression involved in their metabolic pathways.

The effect of long term under- and over-feeding on the expression of genes related to glucose metabolism in mammary tissue of sheep

Glucose utilisation for lactose synthesis in the mammary gland involves expression of a large number of genes whose nutritional regulation remains poorly defined. In this study, the effect of long term under- and over-feeding on the expression of genes [glucose transporter 1: GLUT1, glucose transporter 3: GLUT3, Sodium glucose contransporter 1: SGLT1, two isoforms of β- (1,4) galactosyltransferase: β- (1,4) GAT1, β- (1,4) GAT3 and α-lactalbumin: LALBA] related to glucose metabolism in sheep mammary tissue (MT) was examined. Twenty-four lactating dairy sheep were divided into three homogenous sub-groups and fed the same ration in quantities which met 70% (underfeeding), 100% (control) and 130% (overfeeding) of their energy and crude protein requirements. The results showed a significant reduction on mRNA of GLUT1 and LALBA gene in the MT of underfed sheep, compared with the respective controls and overfed and a significant reduction on mRNA level of SGLT1 and β- (1,4) GAT1 in the MT of underfed sheep, compared with the overfed ones. A significant increase in the GLUT3 mRNA accumulation in the MT of both under- and over- fed sheep was found. Additionally, a trend of reduction on β- (1,4) GAT3 mRNA level in the MT of the underfed sheep, compared with the overfed, was observed. A close positive relationship was obtained between the mRNA transcripts accumulation of GLUT1, SGLT1, β- (1,4) GAT1 and LALBA gene with the milk lactose content and milk lactose yield respectively. In conclusion, feeding level and consequently nutrient availability, may affect glucose uptake and utilisation in sheep MT by altering the GLUT1, GLUT3, SGLT1, β- (1,4) GAT1 and LALBA gene expression involved in their metabolic pathways.

A genome-wide association study for clinical mastitis in first parity US Holstein cows using single-step approach and genomic matrix re-weighting procedure

Clinical mastitis (CM) is one of the health disorders with large impacts on dairy farming profitability and animal welfare. The objective of this study was to perform a genome-wide association study (GWAS) for CM in first-lactation Holstein. Producer-recorded mastitis event information for 103,585 first-lactation cows were used, together with genotype information on 1,361 bulls from the Illumina BovineSNP50 BeadChip. Single-step genomic-BLUP methodology was used to incorporate genomic data into a threshold-liability model. Association analysis confirmed that CM follows a highly polygenic mode of inheritance. However, 10-adjacent-SNP windows showed that regions on chromosomes 2, 14 and 20 have impacts on genetic variation for CM. Some of the genes located on chromosome 14 (LY6K, LY6D, LYNX1, LYPD2, SLURP1, PSCA) are part of the lymphocyte-antigen-6 complex (LY6) known for its neutrophil regulation function linked to the major histocompatibility complex. Other genes on chromosome 2 were also involved in regulating immune response (IFIH1, LY75, and DPP4), or are themselves regulated in the presence of specific pathogens (ITGB6, NR4A2). Other genes annotated on chromosome 20 are involved in mammary gland metabolism (GHR, OXCT1), antibody production and phagocytosis of bacterial cells (C6, C7, C9, C1QTNF3), tumor suppression (DAB2), involution of mammary epithelium (OSMR) and cytokine regulation (PRLR). DAVID enrichment analysis revealed 5 KEGG pathways. The JAK-STAT signaling pathway (cell proliferation and apoptosis) and the 'Cytokine-cytokine receptor interaction' (cytokine and interleukines response to infectious agents) are co-regulated and linked to the 'ABC transporters' pathway also found here. Gene network analysis performed using GeneMania revealed a co-expression network where 665 interactions existed among 145 of the genes reported above. Clinical mastitis is a complex trait and the different genes regulating immune response are known to be pathogen-specific. Despite the lack of information in this study, candidate QTL for CM were identified in the US Holstein population.

DNA methylation and transcription in a distal region upstream from the bovine AlphaS1 casein gene after once or twice daily milking

Once daily milking (ODM) induces a reduction in milk production when compared to twice daily milking (TDM). Unilateral ODM of one udder half and TDM of the other half, enables the study of underlying mechanisms independently of inter-individual variability (same genetic background) and of environmental factors. Our results show that in first-calf heifers three CpG, located 10 kb upstream from the CSN1S1 gene were methylated to 33, 34 and 28%, respectively, after TDM but these levels were higher after ODM, 38, 38 and 33%, respectively. These methylation levels were much lower than those observed in the mammary gland during pregnancy (57, 59 and 50%, respectively) or in the liver (74, 78 and 61%, respectively). The methylation level of a fourth CpG (CpG4), located close by (29% during TDM) was not altered after ODM. CpG4 methylation reached 39.7% and 59.5%, during pregnancy or in the liver, respectively. CpG4 is located within a weak STAT5 binding element, arranged in tandem with a second high affinity STAT5 element. STAT5 binding is only marginally modulated by CpG4 methylation, but it may be altered by the methylation levels of the three other CpG nearby. Our results therefore shed light on mechanisms that help to explain how milk production is almost, but not fully, restored when TDM is resumed (15.1±0.2 kg/day instead of 16.2±0.2 kg/day, p<0.01). The STAT5 elements are 100 bp away from a region transcribed in the antisense orientation, in the mammary gland during lactation, but not during pregnancy or in other reproductive organs (ovary or testes). We now need to clarify whether the transcription of this novel RNA is a consequence of STAT5 interacting with the CSN1S1 distal region, or whether it plays a role in the chromatin structure of this region.