Comparative expression profiling of E. coli and S. aureus inoculated primary mammary gland cells sampled from cows with different genetic predispositions for somatic cell score

Analysis of potential genes, immunological and antioxidant profiles associated with trypanosomiasis susceptibility in dromedary camels

Trypanosomiasis is associated with tissue damage and may trigger an immunological response. These tissue lesions are linked to metabolic issues and oxidative stress. The current study aimed to investigate the immunological, antioxidant, and metabolic changes that may be connected to camel trypanosomiasis. Blood samples were collected from 54 camels and allocated into two groups: The control group (35 camels) and the infected group (19 camels). The genes TLR2, TLR5, IL-17, MARCHF3, RASGRP1, EPS15L1, PPIE, ASB16, CMPK2, LPCAT1, FPGT, GPHN, TNNI3K, DIO3, keap1, and OXSR1 were significantly up-regulated in trypanosomiasis camels. However, down-regulation was observed for the genes Nrf2, PRDX6, and NDUFS5. PCR-DNA sequencing was used to identify nucleotide sequence polymorphisms in the immune (TLR2, TLR5, IL-17, MARCHF3, RASGRP1, and EPS15L1), metabolic (PPIE, ASB16, CMPK2, LPCAT1, FPGT, GPHN, TNNI3K, and DIO3), and antioxidant (Nrf2, Keap1, PRDX6, NDUFS5, and OXSR1) genes between healthy and trypanosomiasis-affected camels. Exploring the serum profile also showed a significant (P ˂ 0.05) increase in Hp, SAA, Cp, IL-1β, IL-6, IL 10, TNF-α, and MDA, with significant (P ˂ 0.05) reduction in the serum levels of CAT, SOD, GSH, T3, and T4 in diseased camels compared with healthy ones. Our findings confirm the significance of nucleotide variations, gene expression patterns, and the biochemical profile of the investigated markers as indicators for the susceptibility of trypanosomiasis in dromedary camels and may be utilized to create management strategies.

Genetic parameters and genome-wide association studies for mozzarella and milk production traits, lactation length, and lactation persistency in Murrah buffaloes

Genetic and genomic analyses of longitudinal traits related to milk production efficiency are paramount for optimizing water buffaloes breeding schemes. Therefore, this study aimed to (1) compare single-trait random regression models under a single-step genomic BLUP setting based on alternative covariance functions (i.e., Wood, Wilmink, and Ali and Schaeffer) to describe milk (MY), fat (FY), protein (PY), and mozzarella (MZY) yields, fat-to-protein ratio (FPR), somatic cell score (SCS), lactation length (LL), and lactation persistency (LP) in Murrah dairy buffaloes (Bubalus bubalis); (2) combine the best functions for each trait under a multiple-trait framework; (3) estimate time-dependent SNP effects for all the studied longitudinal traits; and (4) identify the most likely candidate genes associated with the traits. A total of 323,140 test-day records from the first lactation of 4,588 Murrah buffaloes were made available for the study. The model included the average curve of the population nested within herd-year-season of calving, systematic effects of number of milkings per day, and age at first calving as linear and quadratic covariates, and additive genetic, permanent environment, and residual as random effects. The Wood model had the best goodness of fit based on the deviance information criterion and posterior model probabilities for all traits. Moderate heritabilities were estimated over time for most traits (0.30 ± 0.02 for MY; 0.26 ± 0.03 for FY; 0.45 ± 0.04 for PY; 0.28 ± 0.05 for MZY; 0.13 ± 0.02 for FPR; and 0.15 ± 0.03 for SCS). The heritability estimates for LP ranged from 0.38 ± 0.02 to 0.65 ± 0.03 depending on the trait definition used. Similarly, heritabilities estimated for LL ranged from 0.10 ± 0.01 to 0.14 ± 0.03. The genetic correlation estimates across days in milk (DIM) for all traits ranged from -0.06 (186-215 DIM for MY-SCS) to 0.78 (66-95 DIM for PY-MZY). The SNP effects calculated for the random regression model coefficients were used to estimate the SNP effects throughout the lactation curve (from 5 to 305 d). Numerous relevant genomic regions and candidate genes were identified for all traits, confirming their polygenic nature. The candidate genes identified contribute to a better understanding of the genetic background of milk-related traits in Murrah buffaloes and reinforce the value of incorporating genomic information in their breeding programs.

Genomic Characterization and Initial Insight into Mastitis-Associated SNP Profiles of Local Latvian Bos taurus Breeds

Latvia has two local Bos taurus breeds-Latvian Brown (LBG) and Latvian Blue (LZG)-characterized by a good adaptation to the local climate, longevity, and high fat and protein contents in milk. Since these are desired traits in the dairy industry, this study investigated the genetic background of the LBG and LZG breeds and identified the genetic factors associated with mastitis. Blood and semen samples were acquired, and whole genome sequencing was then performed to acquire a genomic sequence with at least 35× or 10× coverage. The heterozygosity, nucleotide diversity, and LD analysis indicated that LBG and LZG cows have similar levels of genetic diversity compared to those of other breeds. An analysis of the population structure revealed that each breed clustered together, but the overall differentiation between the breeds was small. The highest genetic variance was observed in the LZG breed compared with the LBG breed. Our results show that SNP rs721295390 is associated with mastitis in the LBG breed, and SNPs rs383806754, chr29:43998719CG>C, and rs462030680 are associated with mastitis in the LZG breed. This study shows that local Latvian LBG and LZG breeds have a pronounced genetic differentiation, with each one suggesting its own mastitis-associated SNP profile.

Neutrophils expressing major histocompatibility complex class II molecules circulate in blood and milk during mastitis and show high microbicidal activity

Bovine mastitis is mainly caused by bacterial infection and is responsible for important economic losses as well as alterations of the health and welfare of animals. The increase in somatic cell count (SCC) in milk during mastitis is mainly due to the influx of neutrophils, which have a crucial role in the elimination of pathogens. For a long time, these first-line defenders have been viewed as microbe killers, with a limited role in the orchestration of the immune response. However, their role is more complex: we recently characterized a bovine neutrophil subset expressing major histocompatibility complex class II (MHC-II) molecules (MHC-II^pos), usually distributed on antigen-presenting cells, as having regulatory capacities in cattle. In this study, our objective was to evaluate the implication of different neutrophils subsets in the mammary gland immunity during clinical and subclinical mastitis. Using flow cytometry, we analyzed the presence of MHC-II^pos neutrophils in blood and in milk during clinical mastitis at different time points of inflammation (n = 10 infected quarters) and during subclinical mastitis, defined as the presence of bacteria and an SCC >150,000 cells/mL (n = 27 infected quarters). Our results show, for the first time, that in blood and milk, neutrophils are a heterogeneous population and encompass at least 2 subsets distinguishable by their expression of MHC-II. In milk without mastitis, we observed higher production of reactive oxygen species and higher phagocytosis capacity of MHC-II^pos neutrophils compared with their MHC-II^neg counterparts, indicating the high bactericidal capacities of MHC-II^pos neutrophils. MHC-II^pos neutrophils are enriched in milk compared with blood during subclinical mastitis but not during clinical mastitis. Moreover, we observed a positive and highly significant correlation between MHC-II^pos neutrophils and T lymphocytes present in milk during subclinical mastitis. Our experiments involved a total of 47 cows (40 Holstein and 7 Normande cows). To conclude, our study opens the way to the discovery of new biomarkers of mastitis inflammation.

Candidate genes for mastitis resistance in dairy cattle: a data integration approach

BACKGROUND

Inflammation of the mammary tissue (mastitis) is one of the most detrimental health conditions in dairy ruminants and is considered the most economically important infectious disease of the dairy sector. Improving mastitis resistance is becoming an important goal in dairy ruminant breeding programmes. However, mastitis resistance is a complex trait and identification of mastitis-associated alleles in livestock is difficult. Currently, the only applicable approach to identify candidate loci for complex traits in large farm animals is to combine different information that supports the functionality of the identified genomic regions with respect to a complex trait.

METHODS

To identify the most promising candidate loci for mastitis resistance we integrated heterogeneous data from multiple sources and compiled the information into a comprehensive database of mastitis-associated candidate loci. Mastitis-associated candidate genes reported in association, expression, and mouse model studies were collected by searching the relevant literature and databases. The collected data were integrated into a single database, screened for overlaps, and used for gene set enrichment analysis.

RESULTS

The database contains candidate genes from association and expression studies and relevant transgenic mouse models. The 2448 collected candidate loci are evenly distributed across bovine chromosomes. Data integration and analysis revealed overlaps between different studies and/or with mastitis-associated QTL, revealing promising candidate genes for mastitis resistance.

CONCLUSION

Mastitis resistance is a complex trait influenced by numerous alleles. Based on the number of independent studies, we were able to prioritise candidate genes and propose a list of the 22 most promising. To our knowledge this is the most comprehensive database of mastitis associated candidate genes and could be helpful in selecting genes for functional validation studies.

New Insights on Nucleotide Sequence Variants and mRNA Levels of Candidate Genes Assessing Resistance/Susceptibility to Mastitis in Holstein and Montbéliarde Dairy Cows

A major factor in the propagation of an infectious disease is host genetics. In this study, 180 dairy cows (90 of each breed: Holstein and Montbéliarde) were used. Each breed's tested dairy cows were divided into two groups of comparable size (45 cows each), mastitis-free and mastitis-affected groups. Each cow's jugular vein was punctured to obtain blood samples for DNA and RNA extraction. In the examined Holstein and Montbéliarde dairy cows, single nucleotide polymorphisms (SNPs) related with mastitis resistance/susceptibility were found in the RASGRP1, NFkB, CHL1, MARCH3, PDGFD, MAST3, EPS15L1, C1QTNF3, CD46, COX18, NEURL1, PPIE, and PTX3 genes. Chi-square analysis of identified SNPs revealed a significant difference in gene frequency between mastitic and healthy cows. Except for CHL1, mastitic dairy cows of two breeds had considerably higher mRNA levels of the examined genes than did healthy ones. Marker-assisted selection and monitoring of dairy cows' susceptibility to mastitis may be accomplished through the use of discovered SNPs and changes in the gene expression profile of the studied genes. These findings also point to a possible method for reducing mastitis in dairy cows through selective breeding of animals using genetic markers linked to an animal's ability to resist infection.

Bta-miR-199a-3p Inhibits LPS-Induced Inflammation in Bovine Mammary Epithelial Cells via the PI3K/AKT/NF-κB Signaling Pathway

Mastitis is characterized by inflammatory damage to mammary gland tissue, which could decline milk production and quality and significantly affect the economic benefits of ranching. MicroRNAs (miRNAs), such as miR-199a-3p, are novel therapeutic targets in inflammation, and their regulation is an effective strategy for inflammation control. Despite its importance in humans and animals, the molecular mechanism of bovine miR-199a-3p (bta-miR-199a-3p) in dairy cow mastitis and bovine mammary epithelial cell (bMEC) inflammation is unclear. In our study, a bovine mammary epithelial cell line (MAC-T) induced by lipopolysaccharide (LPS) was used as an inflammatory cell model to investigate the molecular mechanism of bta-miR-199a-3p in the MAC-T inflammatory response. bta-miR-199a-3p was up-regulated in the LPS-induced MAC-T cells, while CD2-associated protein (CD2AP) was revealed as its target gene in a double luciferase reporter gene experiment. In addition, the overexpression of bta-miR-199a-3p negatively regulated the expression of CD2AP and the activation of the phosphatidylinositol 3-kinase (PI3K)/AKT/nuclear factor kappa-B (NF-κB) signaling pathway. These subsequently inhibited the secretion of related inflammatory factors (TNF-α, IL-1β, and IL-6) and the expression of apoptotic genes (CASP3 and CASP9), thereby alleviating the LPS-challenged inflammatory response in the MAC-T cells. Silencing of bta-miR-199a-3p, however, reversed the above effects. Thus, bta-miR-199a-3p inhibits LPS-induced inflammation in bMECs by directly targeting CD2AP and regulating the PI3K/AKT/NF-κB signaling pathway. This study reveals the potential regulatory mechanism of bta-miR-199a-3p in bMEC inflammatory immune response and may serve as a useful target for the treatment of mastitis.

Integrative Systems Biology Analysis Elucidates Mastitis Disease Underlying Functional Modules in Dairy Cattle

Background: Mastitis is the most prevalent disease in dairy cattle and one of the most significant bovine pathologies affecting milk production, animal health, and reproduction. In addition, mastitis is the most common, expensive, and contagious infection in the dairy industry. Methods: A meta-analysis of microarray and RNA-seq data was conducted to identify candidate genes and functional modules associated with mastitis disease. The results were then applied to systems biology analysis via weighted gene coexpression network analysis (WGCNA), Gene Ontology, enrichment analysis for the Kyoto Encyclopedia of Genes and Genomes (KEGG), and modeling using machine-learning algorithms. Results: Microarray and RNA-seq datasets were generated for 2,089 and 2,794 meta-genes, respectively. Between microarray and RNA-seq datasets, a total of 360 meta-genes were found that were significantly enriched as "peroxisome," "NOD-like receptor signaling pathway," "IL-17 signaling pathway," and "TNF signaling pathway" KEGG pathways. The turquoise module (n = 214 genes) and the brown module (n = 57 genes) were identified as critical functional modules associated with mastitis through WGCNA. PRDX5, RAB5C, ACTN4, SLC25A16, MAPK6, CD53, NCKAP1L, ARHGEF2, COL9A1, and PTPRC genes were detected as hub genes in identified functional modules. Finally, using attribute weighting and machine-learning methods, hub genes that are sufficiently informative in Escherichia coli mastitis were used to optimize predictive models. The constructed model proposed the optimal approach for the meta-genes and validated several high-ranked genes as biomarkers for E. coli mastitis using the decision tree (DT) method. Conclusion: The candidate genes and pathways proposed in this study may shed new light on the underlying molecular mechanisms of mastitis disease and suggest new approaches for diagnosing and treating E. coli mastitis in dairy cattle.

Comparison of the response of mammary gland tissue from two divergent lines of goat with high and low milk somatic cell scores to an experimental Staphylococcus aureus infection

Mastitis represents one of the major economic and health threats to the livestock sector associated with reduction in milk quality, loss of production and is a major reason for culling. Somatic cell score (SCS) is used as a criterion in breeding programmes to select cows genetically less susceptible to mastitis. The relevance of SCS as a predictor of udder health and susceptibility to mastitis is still untested in goats. In this study, two lines of French Alpine goats selected for extreme breeding values for somatic cell scores, one line with high SCS (HSCS) and the other with low SCS (LSCS), were used to test the hypothesis that the mammary response and function differed between the lines. The aim of the present study was to investigate differences in the early immune response in caprine mammary gland tissues challenged with Staphylococcus aureus, one of the main pathogens responsible for the intra-mammary infection in small ruminants, using transcriptomic and histopathology analyses. The comparison between HSCS and LSCS goat lines, showed differences in the response at the histological level for inflammation, presence of neutrophils and micro-abscess formation, and at the molecular level in the expression of CXCL8, IL-6, NFKBIZ and IL-1β. CXCL8 and CXCL2 genes, which showed a higher level of expression in the experimentally infected HSCS line. The molecular data and histopathology both suggested that following S. aureus infection, mobilization, recruitment, infiltration, and chemotaxis of neutrophil, leads to a more severe inflammation in the HSCS compared to LSCS animals. Our results represent an initial basis for further studies to unravel the genetic basis of early mastitis inflammatory responses and the selection of dairy animals more resistant to bacterial mastitis.

Anti-mastitis SNV identification of NFκB1 in Chinese Holstein cows and the possible anti-inflammation role of NFκB1/p105 in bovine MECs

NFκB1/p105 is the critical member of the NFκB family which can suppress inflammation, ageing, and cancer when p50/p50 homodimer is formed. Currently, the research about the role of NFκB1/p105 during cow mastitis is limited. Here, we analyzed the correlation of six single-nucleotide variants of the NFκB1 gene with somatic cell count, milk yield, milk fat content, and milk protein content in 547 Chinese Holstein cows, and explored the mRNA expression profiles of the NFκB family and ubiquitin ligases (βTrCP1, βTrCP2, KPC1, KPC2) in LPS-induced bovine mammary epithelial cells (MECs) by transcriptome-Seq. The association analysis showed that cows with SNV2-TT and SNV6-CC in the NFκB1 gene had significantly higher milk protein content (P < 0.05), while cows with SNV5-TT in the NFκB1 gene had significantly lower somatic cell score (SCS), but CC genotype at SNV5 locus was not detected in our Holstein cows. The transcriptome-Seq results demonstrated the mRNA expression of NFκB1 was increased and peaked at 4 h post-induction, while the mRNA expressions of both KPC1 and BCL3 that promote the anti-inflammation function of NFκB1/p105 were decreased in LPS-induced bovine MECs. TNFAIP3, an inhibitor of both degradation and processing of p105 precursor, was markedly increased by more than 3 folds. Furthermore, bta-miR-125b which targets at the 3'UTR of TNFAIP3 was reduced by 50%. These results indicated that SNV5-TT of the NFκB1 gene with lower SCS may be an anti-mastitis genotype that could cope with infection more efficiently in Chinese Holstein cows. In addition, the anti-inflammation role of NFκB1/p105 seemed to be inhibited in LPS-induced-bovine MECs because the formation of the p50/p50 homodimer was arrested. This study provides a new perspective to understand the inflammatory mechanism in dairy mastitis.

Identification of Gene Modules and Hub Genes Involved in Mastitis Development Using a Systems Biology Approach

Objective

Mastitis is defined as the inflammation of the mammary gland, which impact directly on the production performance and welfare of dairy cattle. Since, mastitis is a multifactorial complex disease and the molecular pathways underlying this disorder have not been clearly understood yet, a system biology approach was used in this study to a better understanding of the molecular mechanisms behind mastitis.

Methods

Publicly available RNA-Seq data containing samples from milk of five infected and five healthy Holstein cows at five time points were retrieved. Gene Co-expression network analysis (WGCNA) approach and functional enrichment analysis were then applied with the aim to find the non-preserved module of genes that their connectivity were altered under infected condition. Hub genes were identified in the non-preserved modules and were subjected to protein-protein interactions (PPI) network construction.

Results

Among the 25 modules identified, eight modules were non-preserved and were also biologically associated with inflammation, immune response and mastitis development. Interestingly most of the hub genes in the eight modules were also densely connected in the PPI network. Of the hub genes, 250 genes were hubs in both co-expression and PPI networks and most of them were reported to play important roles in immune response or inflammatory pathways. The blue module was highly enriched in inflammatory responses and STAT1 was suggested to play an important role in mastitis development by regulating the immune related genes in this module. Moreover, a set of highly connected genes were identified such as BIRC3, PSMA6, FYN, F11R, NFKBIZ, NFKBIA, GRO1, PHB, CD3E, IL16, GSN, SOCS2, HCK, VAV1 and TLR6, which have been established to be critical for mastitis pathogenesis.

Conclusion

This study improved the understanding of the mechanisms underlying bovine mastitis and suggested eight non-preserved modules along with several most important genes with promising potential in etiology of mastitis.

Transcriptome Analysis of The Inflammatory Responses of Bovine Mammary Epithelial Cells: Exploring Immunomodulatory Target Genes for Bovine Mastitis

Bovine mastitis is the inflammatory reaction of the mammary gland and is commonly caused by bacterial infections in high-yielding dairy cows. The detailed investigation of the immunotranscriptomic response of bovine mammary epithelial (BME) cells to pattern recognition receptors (PRRs) activation by microbial-associated molecular patterns (MAMPs) can be of great importance for understanding the innate immune defense mechanisms, and for exploring the immunomodulatory candidate genes. In this work, we investigated the transcriptome modifications of BME cells after the in vitro stimulation with Escherichia coli derived lipopolysaccharide (LPS) and heat-killed Staphylococcus aureus JE2 and S. aureus SA003. In addition, the effect of Pam3CSK4 (a synthetic triacylated lipopeptide that activates Toll-like receptor 2 (TLR2)), and the intracellular chemotactic protein cyclophilin A (CyPA), which is secreted by BME cells during mastitis, in the expression changes of selected cytokines and chemokines were evaluated by qPCR. Microarray analysis identified 447, 465 and 520 differentially expressed genes (DEGs) in the BME cells after LPS, S. aureus JE2 and S. aureus SA003 stimulation, respectively. A major differential response in the inflammatory gene expression was noticed between the stimulation of LPS and S. aureus strains. Unlike the S. aureus strains, LPS stimulation resulted in significant upregulation of CCL2, CXCL2, CXCL3, CXCL8,IL1α and IL1β, which were confirmed by qPCR analysis. Pam3CSK4 was not able to induce significant changes in the expression of cytokines and chemokines in challenged BME cells. The exogenous CyPA administration was able to upregulate CXCL2, CXCL3, CXCL8, IL1α and IL1β expression in BME cells indicating its ability to promote inflammation. The identification of transcriptional markers of mastitis specific for individual inflammatory factors such as LPS, Pam3CSK4 or CyPA, which can be evaluated in vitro in BME cells, may enable the development of novel diagnostics and/or immunomodulatory treatments, providing new tools for the effective management of mastitis in dairy cows. The results of this work are an advance in this regard.

In vivo model to study the impact of genetic variation on clinical outcome of mastitis in uniparous dairy cows

Background

In dairy herds, mastitis causes detrimental economic losses. Genetic selection offers a sustainable tool to select animals with reduced susceptibility towards postpartum diseases. Studying underlying mechanisms is important to assess the physiological processes that cause differences between selected haplotypes. Therefore, the objective of this study was to establish an in vivo infection model to study the impact of selecting for alternative paternal haplotypes in a particular genomic region on cattle chromosome 18 for mastitis susceptibility under defined conditions in uniparous dairy cows.

Results

At the start of pathogen challenge, no significant differences between the favorable (Q) and unfavorable (q) haplotypes were detected. Intramammary infection (IMI) with Staphylococcus aureus 1027 (S. aureus, n = 24, 96 h) or Escherichia coli 1303 (E. coli, n = 12, 24 h) was successfully induced in all uniparous cows. This finding was confirmed by clinical signs of mastitis and repeated recovery of the respective pathogen from milk samples of challenged quarters in each animal. After S. aureus challenge, Q-uniparous cows showed lower somatic cell counts 24 h and 36 h after challenge (P < 0.05), lower bacterial shedding in milk 12 h after challenge (P < 0.01) and a minor decrease in total milk yield 12 h and 24 h after challenge (P < 0.01) compared to q-uniparous cows.

Conclusion

An in vivo infection model to study the impact of genetic selection for mastitis susceptibility under defined conditions in uniparous dairy cows was successfully established and revealed significant differences between the two genetically selected haplotype groups. This result might explain their differences in susceptibility towards IMI. These clinical findings form the basis for further in-depth molecular analysis to clarify the underlying genetic mechanisms for mastitis resistance.

Transcriptomic analysis on the promoter regions discover gene networks involving mastitis in cattle

Mastitis is one of the costliest diseases in dairy farms caused by infection of different microorganisms such as Escherichia coli, Streptococcus uberis and Staphylococcus aureus. Promoters are significantly involved in regulating gene expression and shedding light on the mechanisms of transcriptional regulation in physiological and immunological processes of the infections. Exploiting regulatory elements such as transcription factor binding sites (TFBSs modules) on the promoter region could reveal co-regulated genes, which allow screating regulatory models and executing a cross-sectional analysis on several databases. In this study, the promoter regions of 11 genes associated with contagious mastitis including CCL4, CXCL8, STAT3, IKBKB, MAPK14, NFKBIA, NFKB1, TNF, IL18, IL6, and HCK were investigated to predict the activating regulatory modules on promoters and to discover the key related transcription factors. By exploring the promoter regions, 228 genes were discovered comprising the same transcription factors modules. Out of 228 genes, 36 were validated using five microarray datasets. The promoter research of these genes revealed that as many as 7 down-regulated and 12 up-regulated genes are predictable in the network. The genes whose functions were associated with the initial gene list (11 genes), were identified by DAVID queries with TFBSs models implying that the approach provides a clear image of the underlying regulatory mechanism of gene expression profile and offers a novel approach in designing gene networks in cattle.

Genetic selection for bovine chromosome 18 haplotypes associated with divergent somatic cell score affects postpartum reproductive and metabolic performance

The susceptibility of animals to periparturient diseases has a great effect on the economic efficiency of dairy industries, on the frequency of antibiotic treatment, and on animal welfare. The use of selection for breeding cows with reduced susceptibility to diseases offers a sustainable tool to improve dairy cattle farming. Several studies have focused on the association of distinct bovine chromosome 18 genotypes or haplotypes with performance traits. The aim of this study was to test whether selection of Holstein Friesian heifers via SNP genotyping for alternative paternal chromosome 18 haplotypes associated with favorable (Q) or unfavorable (q) somatic cell scores influences postpartum reproductive and metabolic diseases. Thirty-six heifers (18 Q and 18 q) were monitored from 3 wk before calving until necropsy on d 39 (± 4 d) after calving. Health status and rectal temperature were measured daily, and body condition score and body weight were assessed once per week. Blood samples were drawn twice weekly, and levels of insulin, nonesterified fatty acids, insulin-like growth factor-I, growth hormone, and β-hydroxybutyrate were measured. Comparisons between the groups were performed using Fisher's exact test, chi-squared test, and the GLIMMIX procedure in SAS. Results showed that Q-heifers had reduced incidence of metritis compared with q-heifers and were less likely to develop fever. Serum concentrations of β-hydroxybutyrate were lower and insulin-like growth factor-I plasma concentrations were higher in Q- compared with q-heifers. However, the body condition score and withers height were comparable between haplotypes, but weight loss tended to be lower in Q-heifers compared with q-heifers. No differences between the groups were detected concerning retained fetal membranes, uterine involution, or onset of cyclicity. In conclusion, selection of chromosome 18 haplotypes associated with a reduced somatic cell score resulted in a decreased incidence of postpartum reproductive and metabolic diseases in this study. The presented data add to the existing knowledge aimed at avoiding negative consequences of genetic selection strategies in dairy cattle farming. The underlying causal mechanisms modulated by haplotypes in the targeted genomic region and immune competence necessitate further investigation.

Characterization of functional traits with focus on udder health in heifers with divergent paternally inherited haplotypes on BTA18

Background

A major challenge in modern medicine and animal husbandry is the issue of antimicrobial resistance. One approach to solving this potential medical hazard is the selection of farm animals with less susceptibility to infectious diseases. Recent advances in functional genome analysis and quantitative genetics have opened the horizon to apply genetic marker information for efficiently identifying animals with preferential predisposition regarding health traits. The current study characterizes functional traits with a focus on udder health in dairy heifers. The animals were selected for having inherited alternative paternal haplotypes for a genomic region on Bos taurus chromosome (BTA) 18 genetically associated with divergent susceptibility to longevity and animal health, particularly mastitis.

Results

In the first weeks of lactation, the q heifers which had inherited the unfavorable (q) paternal haplotype displayed a significantly higher number of udder quarters with very low somatic cell count (< 10,000 cells / ml) compared to their paternal half-sib sisters with the favorable (Q) paternal haplotype. This might result in impaired mammary gland sentinel function towards invading pathogens. Furthermore, across the course of the first lactation, there was indication that q half-sib heifers showed higher somatic cell counts, a surrogate trait for udder health, in whole milkings compared to their paternal half-sib sisters with the favorable (Q) paternal haplotype. Moreover, heifers with the haplotype Q had a higher feed intake and higher milk yield compared to those with the q haplotype. Results of this study indicate that differences in milk production and calculated energy balance per se are not the main drivers of the genetically determined differences between the BTA18 Q and q groups of heifers.

Conclusions

The paternally inherited haplotype from a targeted BTA18 genomic region affect somatic cell count in udder quarters during the early postpartum period and might also contribute to further aspects of animal’s health and performance traits due to indirect effects on feed intake and metabolism.

Electronic supplementary material

The online version of this article (10.1186/s12917-019-1988-4) contains supplementary material, which is available to authorized users.

Genetic Selection for Mastitis Resistance

Mastitis is a prevalent and costly disease on dairy farms. Improved management and hygiene can reduce the risk of infection by contagious or environmental pathogens, and genetic selection can confer permanent improvement in mastitis resistance. National veterinary recording systems in the Nordic countries have allowed direct selection for sire families with low incidence of clinical mastitis for 3 decades, whereas other countries have practiced indirect selection for lower somatic cell count. Recently, pooling of producer-recorded data from on-farm herd management software programs has enabled selection for reduced incidence of clinical mastitis in the United States and other leading dairy countries.

Identification of several key genes by microarray data analysis of bovine mammary gland epithelial cells challenged with Escherichia coli and Staphylococcus aureus

PURPOSE

This study was aimed at exploring the mechanisms and identifying the key candidate genes associated with S. aureus and E. coli mastitis.

METHODS

A public microarray dataset GSE24560 was downloaded. Differentially expressed genes (DEGs) were screened in E. coli- and S. aureus-infected primary bovine mammary gland epithelial cell (pBMEC) samples, and compared with control samples at 1 h, 6 h, and 24 h. A functional enrichment analysis was performed, and construction of a gene co-expression network was performed based on genes that showed consistent changes over time, which were identified using time series expression analysis. Then, a miRNA/TF regulatory network was constructed based on the genes in the co-expression network. The genes in the miRNA/TF regulatory network were screened for involvement in related diseases. Furthermore, the expression of several selected DEGs was further validated using real-time RT-PCR.

RESULTS

In total, 92 and 81 DEGs showed continuous differential expression over time in the E. coli- and S. aureus-inoculated groups. DEGs in the E. coli-inoculated group were associated with the RIG-I-like receptor signaling pathway, and those in the S. aureus-inoculated group were associated with the lysosome pathway. Time series expression analysis identified two gene clusters. NFKBIZ and GRO1 in the gene co-expression network were associated with inflammatory and defense responses. Moreover, several genes such as CXADR, APP, and CXCL2 in the miRNA/TF regulatory network, were associated with infection, inflammation, or stress-related diseases.

CONCLUSION

RIG-I like receptor pathway and several DEGs such as NFKBIZ, GRO1, CXCL2, and CXADR may play critical roles in the response to infection in pBMECs.

Innate immune response of bovine mammary epithelial cells to Mycoplasma bovis

Mycoplasma spp. are contagious bacteria, and mycoplasmal mastitis is a serious productivity problem on dairy farms. Bovine mammary epithelial cells (bMECs) have an important role in the elimination of pathogens, but the effect of Mycoplasma bovis on bMECs has not been fully described. To elucidate the immune response against intramammary infection by M. bovis, we undertook microarray analysis to examine and profile mRNA expression in bMECs after stimulation with M. bovis. We also compared the effects of M. bovis, Staphylococcus aureus, and Escherichia coli on immune-related mRNA expression in bMECs. Transcriptome analysis indicated a significant decrease in the level of mRNA-encoding lysine-specific demethylase 4D, suggesting that the immune response is suppressed by a decrease in histone demethylase activity. Interleukin (IL)-1β, IL-6, tumor necrosis factor alpha, toll-like receptor (TLR) 2, and TLR4 mRNA expression levels were significantly increased in bMECs stimulated with heat-killed M. bovis, but the expression levels were lower than those following stimulation by heat-killed S. aureus or E. coli. Our results suggest that M. bovis weakly affects mRNA expression in bMECs compared to the effects of E. coli or S. aureus. Moreover, live M. bovis may induce suppression of the immune response in bMECs.

Invited review: Low milk somatic cell count and susceptibility to mastitis

An enduring controversy exists about low milk cell counts and susceptibility to mastitis. The concentration of milk leukocytes, or somatic cell count (SCC), is a well-established direct indicator of mammary gland inflammation that is highly correlated with the presence of a mammary infection. The SCC is also used as a trait for the selection of dairy ruminants less prone to mastitis. As selection programs favor animals with less SCC, and as milk cells contribute to the defense of the mammary gland, the idea that susceptibility to mastitis could possibly be increased in the long term has been put forward and is still widely debated. Epidemiological and experimental studies aimed at relating SCC to susceptibility to mastitis have yielded results that seem contradictory at first sight. Nevertheless, by taking into account the immunobiology of milk and mammary tissue cells and their role in the defense against infection, along with recent studies on SCC-based divergent selection of animals, the issue can be settled. Apparent SCC-linked susceptibility to mastitis is a phenotypic trait that may be linked to immunomodulation but not to selection.

Integration of machine learning and meta-analysis identifies the transcriptomic bio-signature of mastitis disease in cattle

Gram-negative bacteria such as Escherichia coli (E. coli) are assumed to be among the main agents that cause severe mastitis disease with clinical signs in dairy cattle. Rapid detection of this disease is so important in order to prevent transmission to other cows and helps to reduce inappropriate use of antibiotics. With the rapid progress in high-throughput technologies, and accumulation of various kinds of '-omics' data in public repositories, there is an opportunity to retrieve, integrate, and reanalyze these resources to improve the diagnosis and treatment of different diseases and to provide mechanistic insights into host resistance in an efficient way. Meta-analysis is a relatively inexpensive option with good potential to increase the statistical power and generalizability of single-study analysis. In the current meta-analysis research, six microarray-based studies that investigate the transcriptome profile of mammary gland tissue after induced mastitis by E. coli infection were used. This meta-analysis not only reinforced the findings in individual studies, but also several novel terms including responses to hypoxia, response to drug, anti-apoptosis and positive regulation of transcription from RNA polymerase II promoter enriched by up-regulated genes. Finally, in order to identify the small sets of genes that are sufficiently informative in E. coli mastitis, the differentially expressed gene introduced by meta-analysis were prioritized by using ten different attribute weighting algorithms. Twelve meta-genes were detected by the majority of attribute weighting algorithms (with weight above 0.7) as most informative genes including CXCL8 (IL8), NFKBIZ, HP, ZC3H12A, PDE4B, CASP4, CXCL2, CCL20, GRO1(CXCL1), CFB, S100A9, and S100A8. Interestingly, the results have been demonstrated that all of these genes are the key genes in the immune response, inflammation or mastitis. The Decision tree models efficiently discovered the best combination of the meta-genes as bio-signature and confirmed that some of the top-ranked genes -ZC3H12A, CXCL2, GRO, CFB- as biomarkers for E. coli mastitis (with the accuracy 83% in average). This research properly indicated that by combination of two novel data mining tools, meta-analysis and machine learning, increased power to detect most informative genes that can help to improve the diagnosis and treatment strategies for E. coli associated with mastitis in cattle.

Differentiating Staphylococcus aureus from Escherichia coli mastitis: S. aureus triggers unbalanced immune-dampening and host cell invasion immediately after udder infection

The etiology determines quality and extent of the immune response after udder infection (mastitis). Infections with Gram negative bacteria (e.g. Escherichia coli) will quickly elicit strong inflammation of the udder, fully activate its immune defence via pathogen receptor driven activation of IκB/NF-κB signaling. This often eradicates the pathogen. In contrast, Gram-positive bacteria (e.g. Staphylococcus aureus) will slowly elicit a much weaker inflammation and immune response, frequently resulting in chronic infections. However, it was unclear which immune regulatory pathways are specifically triggered by S. aureus causing this partial immune subversion. We therefore compared in first lactating cows the earliest (1–3 h) udder responses against infection with mastitis causing pathogens of either species. Global transcriptome profiling, bioinformatics analysis and experimental validation of key aspects revealed as S. aureus infection specific features the (i) failure to activating IκB/NF-κB signaling; (ii) activation of the wnt/β-catenin cascade resulting in active suppression of NF-κB signaling and (iii) rearrangement of the actin-cytoskeleton through modulating Rho GTPase regulated pathways. This facilitates invasion of pathogens into host cells. Hence, S. aureus mastitis is characterized by eliciting unbalanced immune suppression rather than inflammation and invasion of S. aureus into the epithelial cells of the host causing sustained infection.

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

Meta-Analysis of Transcriptional Responses to Mastitis-Causing Escherichia coli

Bovine mastitis is a widespread disease in dairy cows, and is often caused by bacterial mammary gland infection. Mastitis causes reduced milk production and leads to excessive use of antibiotics. We present meta-analysis of transcriptional profiles of bovine mastitis from 10 studies and 307 microarrays, allowing identification of much larger sets of affected genes than any individual study. Combining multiple studies provides insight into the molecular effects of Escherichia coli infection in vivo and uncovers differences between the consequences of E. coli vs. Staphylococcus aureus infection of primary mammary epithelial cells (PMECs). In udders, live E. coli elicits inflammatory and immune defenses through numerous cytokines and chemokines. Importantly, E. coli infection causes downregulation of genes encoding lipid biosynthesis enzymes that are involved in milk production. Additionally, host metabolism is generally suppressed. Finally, defensins and bacteria-recognition genes are upregulated, while the expression of the extracellular matrix protein transcripts is silenced. In PMECs, heat-inactivated E. coli elicits expression of ribosomal, cytoskeletal and angiogenic signaling genes, and causes suppression of the cell cycle and energy production genes. We hypothesize that heat-inactivated E. coli may have prophylactic effects against mastitis. Heat-inactivated S. aureus promotes stronger inflammatory and immune defenses than E. coli. Lipopolysaccharide by itself induces MHC antigen presentation components, an effect not seen in response to E. coli bacteria. These results provide the basis for strategies to prevent and treat mastitis and may lead to the reduction in the use of antibiotics.

Comparison of the pathogen species-specific immune response in udder derived cell types and their models

The outcome of an udder infection (mastitis) largely depends on the species of the invading pathogen. Gram-negative pathogens, such as Escherichia coli often elicit acute clinical mastitis while Gram-positive pathogens, such as Staphylococcus aureus tend to cause milder subclinical inflammations. It is unclear which type of the immune competent cells residing in the udder governs the pathogen species-specific physiology of mastitis and which established cell lines might provide suitable models. We therefore profiled the pathogen species-specific immune response of different cell types derived from udder and blood. Primary cultures of bovine mammary epithelial cells (pbMEC), mammary derived fibroblasts (pbMFC), and bovine monocyte-derived macrophages (boMdM) were challenged with heat-killed E. coli, S. aureus and S. uberis mastitis pathogens and their immune response was scaled against the response of established models for MEC (bovine MAC-T) and macrophages (murine RAW 264.7). Only E. coli provoked a full scale immune reaction in pbMEC, fibroblasts and MAC-T cells, as indicated by induced cytokine and chemokine expression and NF-κB activation. Weak reactions were induced by S. aureus and none by S. uberis challenges. In contrast, both models for macrophages (boMdM and RAW 264.7) reacted strongly against all the three pathogens accompanied by strong activation of NF-κB factors. Hence, the established cell models MAC-T and RAW 264.7 properly reflected key aspects of the pathogen species-specific immune response of the respective parental cell type. Our data imply that the pathogen species-specific physiology of mastitis likely relates to the respective response of MEC rather to that of professional immune cells.

Extracellular milieu grossly alters pathogen-specific immune response of mammary epithelial cells

BACKGROUND

Considerably divergent data have been published from attempts to model the E. coli vs. S. aureus specific immune reaction of the udder using primary cultures of bovine mammary epithelial cells from cows (pbMEC). Some groups reported a swift, strong and transient inflammatory response against challenges with E. coli and only a weak and retarded response against S. aureus, in agreement with the respective reaction of the udder. Others found almost the reverse. Presence or absence of fetal calf serum distinguished the experimental setting between both groups. We examined here if this causes the divergent reaction of the pbMEC towards both pathogen species. We challenged pbMEC with proteins from heat killed E. coli or S. aureus pathogens or purified TLR2 and TLR4 ligands. The stimuli were applied in normal growth medium with (SM10) or without (SM0) 10% fetal calf serum, or in the basal medium supplemented with 10 mg/ml milk proteins (SM Milk).

RESULTS

Withdrawal of FCS slowed down and decreased the extent by which E. coli or LPS enhanced the expression of cyto- and chemokine encoding genes through impaired TLR4 signalling but enforced their expression during stimulation with S. aureus. SM Milk strongly quenched the induction of those genes. S. aureus strain specific differences in the reaction of the pbMEC could only be recorded in SM0. NF-κB factors were activated by E. coli in all stimulation media, but only to a small extent by S. aureus, solely in SM0. Purified ligands for TLR2 stimulated expression of those genes and activated NF-κB equally well in SM10 and SM0. The mRNA destabilizing factor tristetraproline was only induced by E. coli in SM10 and by purified PAMPs.

CONCLUSIONS

Our data cross validate the correctness of previously published divergent data on the pathogen-specific induction of key immune genes in pbMEC. The differences are due to the presence of FCS, modulating signalling through TLR4 and TLR-unrelated pathogen receptors. S. aureus does not substantially activate any TLR signalling in MEC. Rather, receptors distinct from TLRs perceive the presence of S. aureus and control the immune response against this pathogen in MEC.

Inducible expression of enhanced green fluorescent protein by interleukin-1α, interleukin-1β and Toll-like receptor 2 promoters in goat mammary epithelial cells in response to bacterial challenges

The development of a bacteria-inducible expression system has several advantages compared with persistent expression of anti-bacterial proteins in milk to prevent and treat mastitis. The present study determined whether mastitis responsive promoters could regulate enhanced green fluorescent protein (EGFP) expression in goat mammary epithelial cells (GMECs) in response to challenges with Escherichia coli, Staphylococcus aureus or Streptococcus agalactiae. The level of expression of interleukin (IL)-1α was significantly increased in GMECs challenged with E. coli, S. aureus or S. agalactiae compared with untreated GMECs. IL-1β was induced by E. coli and S. aureus, while Toll-like receptor 2 (TLR2) was induced by E. coli only. GMECs were transfected with IL-1α, IL-1β and TLR2 promoter-EGFP reporter gene lentiviral expression vectors and the levels of expression of EGFP were measured by flow cytometry and Western blot analysis after bacterial challenge. EGFP expression driven by the IL-1α and IL-1β promoters was higher in GMECs challenged with E. coli, S. aureus or S. agalactiae than in untreated GMECs. There were no differences in EGFP expression driven by the TLR2 promoter between GMECs challenged with S. aureus or S. agalactiae and untreated GMECs, but EGFP expression was significantly increased in GMECs challenged with E. coli. Overall, these results indicate that the promoters of some bacteria-inducible genes can regulate EGFP expression in GMECs in response to bacterial challenges. This bacteria-inducible expression strategy could be used for production of mastitis resistant animals by regulating the expression of anti-bacterial proteins in the mammary gland.

Bovine mastitis: frontiers in immunogenetics

Mastitis is one of the most prevalent and costly diseases in the dairy industry with losses attributable to reduced milk production, discarded milk, early culling, veterinary services, and labor costs. Typically, mastitis is an inflammation of the mammary gland most often, but not limited to, bacterial infection, and is characterized by the movement of leukocytes and serum proteins from the blood to the site of infection. It contributes to compromised milk quality and the potential spread of antimicrobial resistance if antibiotic treatment is not astutely applied. Despite the implementation of management practises and genetic selection approaches, bovine mastitis control continues to be inadequate. However, some novel genetic strategies have recently been demonstrated to reduce mastitis incidence by taking advantage of a cow's natural ability to make appropriate immune responses against invading pathogens. Specifically, dairy cattle with enhanced and balanced immune responses have a lower occurrence of disease, including mastitis, and they can be identified and selected for using the high immune response (HIR) technology. Enhanced immune responsiveness is also associated with improved response to vaccination, increased milk, and colostrum quality. Since immunity is an important fitness trait, beneficial associations with longevity and reproduction are also often noted. This review highlights the genetic regulation of the bovine immune system and its vital contributions to disease resistance. Genetic selection approaches currently used in the dairy industry to reduce the incidence of disease are reviewed, including the HIR technology, genomics to improve disease resistance or immune response, as well as the Immunity(+)™ sire line. Improving the overall immune responsiveness of cattle is expected to provide superior disease resistance, increasing animal welfare and food quality while maintaining favorable production levels to feed a growing population.

Non-classical proIL-1beta activation during mammary gland infection is pathogen-dependent but caspase-1 independent

Infection of the mammary gland with live bacteria elicits a pathogen-specific host inflammatory response. To study these host-pathogen interactions wild type mice, NF-kappaB reporter mice as well as caspase-1 and IL-1beta knockout mice were intramammarily challenged with Escherichia coli (E. coli) and Staphylococcus aureus (S. aureus). The murine mastitis model allowed to compare the kinetics of the induced cytokine protein profiles and their underlying pathways. In vivo and ex vivo imaging showed that E. coli rapidly induced NF-kappaB inflammatory signaling concomitant with high mammary levels of TNF-alpha, IL-1 alpha and MCP-1 as determined by multiplex analysis. In contrast, an equal number of S. aureus bacteria induced a low NF-kappaB activity concomitant with high mammary levels of the classical IL-1beta fragment. These quantitative and qualitative differences in local inflammatory mediators resulted in an earlier neutrophil influx and in a more extensive alveolar damage post-infection with E. coli compared to S. aureus. Western blot analysis revealed that the inactive proIL-1beta precursor was processed into pathogen-specific IL-1beta fragmentation patterns as confirmed with IL-1beta knockout animals. Additionally, caspase-1 knockout animals allowed to investigate whether IL-1beta maturation depended on the conventional inflammasome pathway. The lack of caspase-1 did not prevent extensive proIL-1beta fragmentation by either of S. aureus or E. coli. These non-classical IL-1beta patterns were likely caused by different proteases and suggest a sentinel function of IL-1beta during mammary gland infection. Thus, a key signaling nodule can be defined in the differential host innate immune defense upon E. coli versus S. aureus mammary gland infection, which is independent of caspase-1.

Deep sequencing-based transcriptional analysis of bovine mammary epithelial cells gene expression in response to in vitro infection with Staphylococcus aureus stains

Staphylococcus aureus (S. aureus) is an important etiological organism in chronic and subclinical mastitis in lactating cows. Given the fundamental role the primary bovine mammary epithelial cells (pBMECs) play as a major first line of defense against invading pathogens, their interactions with S. aureus was hypothesized to be crucial to the establishment of the latter’s infection process. This hypothesis was tested by investigating the global transcriptional responses of pBMECs to three S. aureus strains (S56,S178 and S36) with different virulent factors, using a tag-based high-throughput transcriptome sequencing technique. Approximately 4.9 million total sequence tags were obtained from each of the three S. aureus-infected libraries and the control library. Referenced to the control, 1720, 219, and 427 differentially expressed unique genes were identified in the pBMECs infected with S56, S178 and S36 S. aureus strains respectively. Gene ontology (GO) and pathway analysis of the S56-infected pBMECs referenced to those of the control revealed that the differentially expressed genes in S56-infected pBMECs were significantly involved in inflammatory response, cell signalling pathways and apoptosis. In the same vein, the clustered GO terms of the differentially expressed genes of the S178-infected pBMECs were found to comprise immune responses, metabolism transformation, and apoptosis, while those of the S36-infected pBMECs were primarily involved in cell cycle progression and immune responses. Furthermore, fundamental differences were observed in the levels of expression of immune-related genes in response to treatments with the three S. aureus strains. These differences were especially noted for the expression of important pro-inflammatory molecules, including IL-1α, TNF, EFNB1, IL-8, and EGR1. The transcriptional changes associated with cellular signaling and the inflammatory response in this study may reflect different immunomodulatory mechanisms that underlie the interaction between pBMECs and S. aureus strains during infection by the latter.

Sequencing the transcriptome of milk production: milk trumps mammary tissue

BACKGROUND

Studies of normal human mammary gland development and function have mostly relied on cell culture, limited surgical specimens, and rodent models. Although RNA extracted from human milk has been used to assay the mammary transcriptome non-invasively, this assay has not been adequately validated in primates. Thus, the objectives of the current study were to assess the suitability of lactating rhesus macaques as a model for lactating humans and to determine whether RNA extracted from milk fractions is representative of RNA extracted from mammary tissue for the purpose of studying the transcriptome of milk-producing cells.

RESULTS

We confirmed that macaque milk contains cytoplasmic crescents and that ample high-quality RNA can be obtained for sequencing. Using RNA sequencing, RNA extracted from macaque milk fat and milk cell fractions more accurately represented RNA from mammary epithelial cells (cells that produce milk) than did RNA from whole mammary tissue. Mammary epithelium-specific transcripts were more abundant in macaque milk fat, whereas adipose or stroma-specific transcripts were more abundant in mammary tissue. Functional analyses confirmed the validity of milk as a source of RNA from milk-producing mammary epithelial cells.

CONCLUSIONS

RNA extracted from the milk fat during lactation accurately portrayed the RNA profile of milk-producing mammary epithelial cells in a non-human primate. However, this sample type clearly requires protocols that minimize RNA degradation. Overall, we validated the use of RNA extracted from human and macaque milk and provided evidence to support the use of lactating macaques as a model for human lactation.

The mammary gland in domestic ruminants: a systems biology perspective

Milk and dairy products are central elements in the human diet. It is estimated that 108kg of milk per year are consumed per person worldwide. Therefore, dairy production represents a relevant fraction of the economies of many countries, being cattle, sheep, goat, water buffalo, and other ruminants the main species used worldwide. An adequate management of dairy farming cannot be achieved without the knowledge on the biological mechanisms behind lactation in ruminants. Thus, understanding the morphology, development and regulation of the mammary gland in health, disease and production is crucial. Presently, innovative and high-throughput technologies such as genomics, transcriptomics, proteomics and metabolomics allow a much broader and detailed knowledge on such issues. Additionally, the application of a systems biology approach to animal science is vastly growing, as new advances in one field of specialization or animal species lead to new lines of research in other areas or/and are expanded to other species. This article addresses how modern research approaches may help us understand long-known issues in mammary development, lactation biology and dairy production.

BIOLOGICAL SIGNIFICANCE

Dairy production depends upon the knowledge of the morphology and regulation of the mammary gland and lactation. High-throughput technologies allow a much broader and detailed knowledge on the biology of the mammary gland. This paper reviews the major contributions that genomics, transcriptomics, metabolomics and proteomics approaches have provided to understand the regulation of the mammary gland in health, disease and production. In the context of mammary gland "omics"-based research, the integration of results using a Systems Biology Approach is of key importance.

Recombinant bovine S100A8 and A9 enhance IL-1β secretion of interferon-gamma primed monocytes

Calgranulin A (S100A8) and B (S100A9) are found at high levels in inflamed tissue and have been associated with acute and chronic inflammatory disorders. Calgranulins are discussed as damage-associated molecular patterns (DAMPs). To analyze the role of calgranulins for inflammatory responses, bovine S100A8 and S100A9 were cloned, successfully expressed and FPLC-purified. Both molecules did not induce NF-κB activation in boTLR4-transfected HEK293 cells and stimulation of bovine monocytes with both proteins did not result in interleukin 1β (IL-1β) secretion or an upregulated mRNA expression of selected genes (IL1B, TNF, CXCL8, IL10, IL12). However, Interferon γ (IFN-γ) primed bovine monocytes released significantly higher amounts of IL-1β after stimulation with S100A8, S100A9, and co-stimulation with adenosine triphosphate (ATP). In IL-4/IL-13-primed monocytes, the IL-1β release was completely abrogated. The results imply that TLR4/MyD88/NF-κB-independent S100A8/A9-mediated activation of the inflammasome in cattle is favored in a Th1 environment and that S100A8 and S100A9 act as a DAMP in cattle.

Differential response of bovine mammary epithelial cells to Staphylococcus aureus or Escherichia coli agonists of the innate immune system

Mastitis caused by Escherichia coli and Staphylococcus aureus is a major pathology of dairy cows. To better understand the differential response of the mammary gland to these two pathogens, we stimulated bovine mammary epithelial cells (bMEC) with either E. coli crude lipopolysaccharide (LPS) or with S. aureus culture supernatant (SaS) to compare the transcriptomic profiles of the initial bMEC response. By using HEK 293 reporter cells for pattern recognition receptors, the LPS preparation was found to stimulate TLR2 and TLR4 but not TLR5, Nod1 or Nod2, whereas SaS stimulated TLR2. Biochemical analysis revealed that lipoteichoic acid, protein A and α-hemolysin were all present in SaS, and bMEC were found to be responsive to each of these molecules. Transcriptome profiling revealed a core innate immune response partly shared by LPS and SaS. However, LPS induced expression of a significant higher number of genes and the fold changes were of greater magnitude than those induced by SaS. Microarray data analysis suggests that the activation pathways and the early chemokine and cytokine production preceded the defense and stress responses. A major differential response was the activation of the type I IFN pathway by LPS but not by SaS. The higher upregulation of chemokines (Cxcl10, Ccl2, Ccl5 and Ccl20) that target mononuclear leucocytes by LPS than by SaS is likely to be related to the differential activation of the type I IFN pathway, and could induce a different profile of the initial recruitment of leucocytes. The MEC responses to the two stimuli were different, as LPS was associated with NF-κB and Fas signaling pathways, whereas SaS was associated with AP-1 and IL-17A signaling pathways. It is noteworthy that at the protein level secretion of TNF-α and IL-1β was not induced by either stimulus. These results suggest that the response of MEC to diffusible stimuli from E. coli and S. aureus contributes to the onset of the response with differential leucocyte recruitment and distinct inflammatory and innate immune reactions of the mammary gland to infection.

T helper 17-associated cytokines are produced during antigen-specific inflammation in the mammary gland

Infectious mastitis cuts down milk production profitability and is a major animal welfare problem. Bacteria-induced inflammation in the mammary gland (MG) is driven by innate immunity, but adaptive immunity can modulate the innate response. Several studies have shown that it is possible to elicit inflammation in the MG by sensitization to an antigen subsequently infused into the lumen of the gland. The objective of our study was to characterize the inflammation triggered in the MG of cows sensitized to ovalbumin, by identifying the cytokines and chemokines likely to play a part in the reaction. Among immunized cows, responders mobilized locally high numbers of leukocytes. An overexpression of the genes encoding IL-17a, IL-17F, IL-21, IL-22 and INF-γ was found in milk cell RNA extracts in the early phase of the inflammatory response. At the protein level, IL-17A was detected in milk as soon as the first sampling time (8 h post-challenge), and both IL-17A and IFN-γ concentrations peaked at 12 to 24 h post-challenge. In mammary tissue from challenged quarters, overexpression of the genes encoding IL-17A, IL-17F, IL-21, IL-22, IL-26 and IFN-γ was observed. Neutrophil-attracting chemokines (CXCL3 and CXCL8) were found in milk, and overexpressed transcripts of chemokines attracting lymphocytes and other mononuclear leukocytes (CXCL10, CCL2, CCL5, CCL20) were detected in mammary tissue. Expression of IL-17A, as revealed by immunohistochemistry, was located in epithelial cells, in leukocytes in the connective tissue and in association with the epithelium, and in migrated alveolar leukocytes of challenged quarters. Altogether, these results show that antigen-specific inflammation in the MG was characterized by the production of IL-17 and IFN-γ. The orientation of the inflammatory response induced by the antigen-specific response has the potential to strongly impact the outcome of bacterial infections of the MG.

Effects of Bos taurus autosome 9-located quantitative trait loci haplotypes on the disease phenotypes of dairy cows with experimentally induced Escherichia coli mastitis

Several quantitative trait loci (QTL) affecting mastitis incidence and mastitis-related traits such as somatic cell score exist in dairy cows. Previously, QTL haplotypes associated with susceptibility to Escherichia coli mastitis in Nordic Holstein-Friesian (HF) cows were identified on Bos taurus autosome 9. In the present study, we induced experimental E. coli mastitis in Danish HF cows to investigate the effect of 2 E. coli mastitis-associated QTL haplotypes on the cows' disease phenotypes and recovery in early lactation. Thirty-two cows were divided in 2 groups bearing haplotypes with either low (HL) or high (HH) susceptibility to E. coli. In addition, biopsies (liver and udder) were collected from half of the cows (n=16), resulting in a 2 × 2 factorial design, with haplotype being one factor (HL vs. HH) and biopsy being the other factor (biopsies vs. no biopsies). Each cow was inoculated with a low E. coli dose (20 to 40 cfu) in one front quarter at time 0 h. Liver biopsies were collected at -144, 12, 24, and 192 h; udder biopsies were collected at 24h and 192 h post-E. coli inoculation. The clinical parameters: feed intake, milk yield, body temperature, heart rate, respiration rate, rumen motility; and the paraclinical parameters: bacterial counts, somatic cell count (SCC), and milk amyloid A levels in milk; and white blood cell count, polymorphonuclear neutrophilic leukocyte (PMNL) count, and serum amyloid A levels in blood were recorded at different time points post-E. coli inoculation. Escherichia coli inoculation changed the clinical and paraclinical parameters in all cows except one that was not infected. Clinically, the HH group tended to have higher body temperature and heart rate than the HL group did. Paraclinically, the HL group had faster PMNL recruitment and SCC recovery than the HH group did. However, we also found interactions between the effects of haplotype and biopsy for body temperature, heart rate, and PMNL. In conclusion, when challenged with E. coli mastitis, HF cows with the specific Bos taurus autosome 9-located QTL haplotypes were associated with differences in leukocyte kinetics, with low-susceptibility cows having faster blood PMNL recruitment and SCC recovery and a tendency for a milder clinical response than the high-susceptibility cows did.

Genetic susceptibility to S. aureus mastitis in sheep: differential expression of mammary epithelial cells in response to live bacteria or supernatant

Staphylococcus aureus is a prevalent pathogen for mastitis in dairy ruminants and is responsible for both clinical and subclinical mastitis. Mammary epithelial cells (MEC) represent not only a physical barrier against bacterial invasion but are also active players of the innate immune response permitting infection clearance. To decipher their functions in general and in animals showing different levels of genetic predisposition to Staphylococcus in particular, MEC from ewes undergoing a divergent selection on milk somatic cell count were stimulated by S. aureus. MEC response was also studied according to the stimulation condition with live bacteria or culture supernatant. The early MEC response was studied during a 5 h time course by microarray to identify differentially expressed genes with regard to the host genetic background and as a function of the conditions of stimulation. In both conditions of stimulation, metabolic processes were altered, the apoptosis-associated pathways were considerably modified, and inflammatory and immune responses were enhanced with the upregulation of il1a, il1b, and tnfa and several chemokines known to enhance neutrophil (cxcl8) or mononuclear leukocyte (ccl20) recruitment. Genes associated with oxidative stress were increased after live bacteria stimulation, whereas immune response-related genes were higher after supernatant stimulation in the early phase. Only 20 genes were differentially expressed between Staphylococcus spp-mastitis resistant and susceptible animals without any clearly defined role on the control of infection. To conclude, this suggests that MEC may not represent the cell type at the origin of the difference of mastitis susceptibility, at least as demonstrated in our genetic model. Supernatant or heat-killed S. aureus produce biological effects that are essentially different from those induced by live bacteria.

Functional adaptations of the transcriptome to mastitis-causing pathogens: the mammary gland and beyond

Application of microarrays to the study of intramammary infections in recent years has provided a wealth of fundamental information on the transcriptomics adaptation of tissue/cells to the disease. Due to its heavy toll on productivity and health of the animal, in vivo and in vitro transcriptomics works involving different mastitis-causing pathogens have been conducted on the mammary gland, primarily on livestock species such as cow and sheep, with few studies in non-ruminants. However, the response to an infectious challenge originating in the mammary gland elicits systemic responses in the animal and encompasses tissues such as liver and immune cells in the circulation, with also potential effects on other tissues such as adipose. The susceptibility of the animal to develop mastitis likely is affected by factors beyond the mammary gland, e.g. negative energy balance as it occurs around parturition. Objectives of this review are to discuss the use of systems biology concepts for the holistic study of animal responses to intramammary infection; providing an update of recent work using transcriptomics to study mammary and peripheral tissue (i.e. liver) as well as neutrophils and macrophage responses to mastitis-causing pathogens; discuss the effect of negative energy balance on mastitis predisposition; and analyze the bovine and murine mammary innate-immune responses during lactation and involution using a novel functional analysis approach to uncover potential predisposing factors to mastitis throughout an animal's productive life.