- Research article
- Open Access
Escherichia coli- and Staphylococcus aureus-induced mastitis differentially modulate transcriptional responses in neighbouring uninfected bovine mammary gland quarters
BMC Genomics volume 14, Article number: 36 (2013)
The most important disease of dairy cattle is mastitis, caused by the infection of the mammary gland by various micro-organisms. Although the transcriptional response of bovine mammary gland cells to in vitro infection has been studied, the interplay and consequences of these responses in the in vivo environment of the mammary gland are less clear. Previously mammary gland quarters were considered to be unaffected by events occurring in neighbouring quarters. More recently infection of individual quarters with mastitis causing pathogens, especially Escherichia coli, has been shown to influence the physiology of neighbouring uninfected quarters. Therefore, the transcriptional responses of uninfected mammary gland quarters adjacent to quarters infected with two major mastitis causing pathogens, E. coli and Staphylococcus aureus, were compared.
The bacteriologically sterile, within-animal control quarters exhibited a transcriptional response to the infection of neighbouring quarters. The greatest response was associated with E. coli infection, while a weaker, yet significant, response occurred during S. aureus infection. The transcriptional responses of these uninfected quarters included the enhanced expression of many genes previously associated with mammary gland infections. Comparison of the transcriptional response of uninfected quarters to S. aureus and E. coli infection identified 187 differentially expressed genes, which were particularly associated with cellular responses, e.g. response to stress. The most affected network identified by Ingenuity Pathway analysis has the immunosuppressor transforming growth factor beta 1 (TGFB1) at its hub and largely consists of genes more highly expressed in control quarters from S. aureus infected cows.
Uninfected mammary gland quarters reacted to the infection of neighbouring quarters and the responses were dependent on pathogen type. Therefore, bovine udder quarters exhibit interdependence and should not be considered as separate functional entities. This suggests that mastitis pathogens not only interact directly with host mammary cells, but also influence discrete sites some distance away, which will affect their response to the subsequent spread of the infection. Understanding the underlying mechanisms may provide further clues for ways to control mammary gland infections. These results also have implications for the design of experimental studies investigating immune regulatory mechanisms in the bovine mammary gland.
Mastitis, the inflammation of the mammary gland, is the most important disease of dairy cattle, with respect to the frequency of occurrence, animal welfare and economic cost, which is estimated to approach $2 billion annually in the US alone (reviewed by ). Mastitis is most frequently caused by infection by a variety of bacteria, of which Staphylococcus aureus and Escherichia coli are amongst the most important gram-positive and gram-negative bacteria respectively (reviewed by ). There is considerable variation in the mastitic disease caused by these two pathogens. E. coli intramammary infections often result in acute mastitis with severe clinical manifestations, which usually resolves within a few days. In contrast, symptoms induced by S. aureus infection are usually less severe, even asymptomatic, and the infection can persist for long periods of time.
Current control strategies largely involve hygiene practices on the farm and the direct administration of antibiotics into udder quarters exhibiting clinical signs of infection. However, many infections are asymptomatic and there is a desire to minimize the use of antibiotics. Therefore, there is a need to develop new and improved control strategies. To achieve this goal our understanding of events that occur during infection must be improved. The use of global transcriptomic studies, either involving the use of microarrays or more recently Next Generation Sequencing, is ideally suited to investigate the complex host-pathogen interplay occurring during mastitis, by simultaneously quantifying the expression of thousands of genes. Several studies have used microarrays to investigate the bovine mammary gland response to infection with particular bacteria, including E. coli[3, 4], Streptococcus uberis[5, 6] and S. aureus. A direct comparison of the induced transcriptional responses identified in these studies is difficult because of differences in the experimental design and the microarray platforms used. However, meta-analyses of such transcriptional studies have recently been reported [8, 9]. The meta-analysis reported by Genini et al.  identified a common transcriptional response in the mammary glands of several species to different pathogens and identified several new pathways not previously identified in the individual studies. However, a considerable amount of information was excluded from this meta-analysis because of the limited overlap of probes included on the microarrays used in the original studies. Therefore, we designed an infection study to investigate and compare the early transcriptional response of the bovine mammary gland to challenge with S. aureus and E. coli using a single microarray platform.
The experimental model used in this study involved the sequential inoculation of udder quarters with bacteria over the course of the study period, leaving one quarter uninfected as a within animal control . The use of within animal controls is a common and accepted practice because the udder quarters are generally considered to be separate, independent anatomical structures [3–7, 11–16] and reduces the between animal variation seen in outbred species such as cattle, as well as addressing the 3Rs, Reduction, Replacement and Refinement, the principles behind the ethical framework for conducting animal experiments. However, analysis of the microarray experiment reported here revealed that E. coli infection of mammary gland quarters can have a profound effect on gene expression in other, uninfected quarters. This conclusion is supported by another transcriptional study comparing uninfected quarters from healthy and E. coli infected cattle . Furthermore, analysis of the microarray data presented here suggests that a similar, but distinct, effect also occurs during S. aureus infection.
Microarray analysis of the transcriptional response of the bovine mammary gland to E. coli and S. aureusinfection
A strictly defined model for E. coli- and S. aureus-induced mastitis was developed involving the sequential infection of three udder quarters whilst the fourth quarter acted as a control . RNA samples collected from each quarter of the twelve animals included in this study were interrogated using the ARK-Genomics B. taurus 20 K cDNA microarray. Initially, the intramammary transcriptional response to E. coli and S. aureus infection was investigated by comparing the expression of genes in infected quarters with those expressed in the uninfected, within animal control quarters. Only five clones exhibited differential expression during S. aureus infection (false discovery rate (FDR) ≤ 0.05, fold difference ≥ 2) (Table 1). The clone AJ813772, which is currently unannotated and only matches EST sequences and an intergenic region of bovine chromosome 6, was significantly down-regulated 2.6 fold at 6 hours post infection. Similarly, two of the four clones that exhibit differential expression at 12 hours post infection are currently unannotated. These clones, AJ819694 and BF775518, were up-regulated 2.2 and 3.7 fold respectively. The other clones match lipopolysaccharide binding protein (LBP) and superoxide dismutase 2 (SOD2), which were up-regulated 4.8 and 5.9 fold respectively. No clones were identified as being differentially expressed at 24 or 72 hours post infection with S. aureus.
The E. coli infection induced clinical symptoms in cattle within 12 hours of infection, including; increased somatic cell count (SCC) in milk from infected quarters, decreased milk production, fever and decreased blood leukocyte counts . It was, therefore, surprising that large numbers of differentially expressed genes were only detected in quarters infected for 24 hours (Table 1). No clones exhibited differential expression at 6 hours post infection, whilst only six were differentially expressed at 12 hours post infection (FDR ≤ 0.05, fold difference ≥ 2). The currently unannotated clone C0882029 was down-regulated 2.1 fold. The remaining five clones were up-regulated and represent; chemokine (C-C motif) ligand 2 (3.9 fold), R3H domain containing 1 (3.3 fold), BTG family, member 2 (3.0 fold), Tribbles homolog 1 (2.4 fold) and chromosome 17 open reading frame 42 (2.3 fold). In contrast, 1048 clones were differentially expressed at 24 hours post infection. The two most up-regulated clones represent S100A8 and S100A12, which exhibited 21.9 and 19.8 fold up-regulation respectively. The two most down-regulated clones represent lipin 1 and WD repeat domain, phosphoinositide interacting 1 (WIPI1), which exhibited 5.4 and 5.2 fold down-regulation respectively.
Preliminary comparison of control quarters
The low numbers of differentially expressed genes (DEG) detected in infected udder quarters, especially those infected with E. coli, were unexpected. Recently over 70 genes were found to be differentially expressed in the bovine mammary gland 12 hours after E. coli infection . Furthermore, 183 genes were found to be differentially expressed 16 hours post S. aureus infection of the bovine mammary gland  using the smaller bovine innate immune microarray . The discrepancy between these results may relate to the strain of bacterium used or the tissue collection procedure. Another potential explanation is that the low number of DEG may be due to the base-line provided by the control quarters. To investigate this possibility the transcriptomes of the E. coli control quarters (EC24T0) and control quarters from both S. aureus infection experiments (SAT0) were compared. Preliminary analysis using GeneSpring (Agilent) highlighted differences between the transcription profiles. Log intensity plots (sample/reference) of EC24T0 and SAT0 had different shapes (Figure 1). Whilst both had normal type distributions centred around one, EC24T0 had a broader distribution than SAT0. This suggests that the transcription profile of SAT0 is more similar to the reference sample, a mixed pool of all the RNA samples, than EC24T0.
Genes exhibiting differential expression in control quarters from animals infected with E. coli and S. aureus
Due to the observed difference in the transcription profiles of RNA from the control quarters of E. coli and S. aureus infected animals, the transcriptomes of these samples were compared. Only the control quarters from the four animals infected with S. aureus for up to 24 hours (SA24T0) were included in this comparison. The analysis revealed that 255 clones exhibited significantly different expression levels in the control quarters (FDR ≤ 0.05, fold difference ≥ 2). The 255 clones, of which 212 are currently annotated, represent 187 genes and are listed in Additional file 1. Therefore, the uninfected quarters, which remained bacteriologically negative , exhibited significantly different transcriptional profiles, which were dependent on whether the animals had been infected in neighbouring quarters with E. coli or S. aureus. Ninety nine clones, representing 77 genes, were expressed at significantly higher levels in SA24T0 than in EC24T0. The clone AJ814393, which exhibited the highest expression level compared to EC24T0, a 4.6 fold difference, is currently unannotated, matching only EST sequences and an intergenic region of bovine chromosome 14. The remaining 156 clones, representing 110 genes, were expressed at significantly higher levels in EC24T0 than in SA24T0. The clone exhibiting the highest expression level compared to SA24T0, with a 6.6 fold difference, represents FK506 binding protein 5 (FKBP5).
The DEG list was examined using the web-based Database of Annotation, Visualization and Integrated Discovery (DAVID) [19, 20]. Human RefSeq accession numbers were used in the analysis due to the limitations in the current bovine annotation. Gene Ontology (GO) Biological Process classifications were assigned to 144 (77%) of the DEG. Gene enrichment analysis was carried out using the default human background and the probability of each functional category was calculated using the modified Fisher Extract test. The two GO terms exhibiting the greatest enrichment were; Cellular Compartment Organization and Biogenesis and Response to Stress (Table 2). Other enriched GO terms included Cell Death and Cell Activation (Table 2). These results indicate that a large cellular response to the infection of neighbouring quarters was occurring in EC24T0 and/or SA24T0.
Closer analysis of the gene list revealed that several genes previously associated with the mammary gland response to bacterial infection are differentially expressed between SA24T0 and EC24T0. These include signal transducer and activator of transcription 3 (STAT3), LBP, CXCL2, BCL2-associated X protein (BAX), B-cell translocation gene 1 (BTG1), metallothioneins and LCN2, which are all expressed at higher levels in EC24T0 samples. This suggests that the EC24T0 samples are responding to infection of neighbouring udder quarters and this accounts for the majority of the transcriptional differences observed between the control quarters.
The DEG were further analyzed using the Ingenuity Pathway Analysis (IPA) software (V7.0, Ingenuity Systems) as described previously . The most striking functional network found to be over-represented in the DEG list involved transforming growth factor beta 1 (TGFB1) and associated molecules (Figure 2). Of the 23 genes in this network that were in the DEG list, 17 were expressed at higher levels in SA24T0 samples than EC24T0 samples. It is unclear if this is due to up-regulation of these genes in SA24T0 or down-regulation in EC24T0, but the involvement of TGFB1 in immunosuppression suggests that regulation of these genes may have a profound effect on the ability of the mammary gland quarter to respond to infection.
Validation of microarray results comparing EC24T0 and SA24T0
Ten genes that exhibited differential expression between EC24T0 and SA24T0 were chosen for validation studies by quantitative RT-PCR (RT-qPCR). Analysis of the microarray results suggested that five of these genes; BTG1, CXCL2, FKBP5, LBP and LCN2, were expressed at higher levels in EC24T0 samples. The RT-qPCR results for all five genes were statistically significant and agreed with the microarray results (Figure 3). The other five genes; ATP-binding cassette, sub-family G (WHITE), member 2 (ABCG2), CXCL12, isocitrate dehydrogenase 1 (NADP+), soluble (IDH1), TGFB1 and xanthine dehydrogenase (XDH), were chosen from those that exhibited greater expression in SA24T0 than EC24T0 samples. The RT-qPCR analysis revealed that all five are expressed at statistically significantly higher levels in SA24T0 than EC24T0 samples (Figure 3), including ABCG2 which was expressed on average 20.9 fold higher levels in SA24T0 than EC24T0 samples.
Genes exhibiting differential expression in control quarters from animals infected with S. aureusfor 24 and 72 hours
To confirm that the variation between udder quarters was due only to events occurring in EC24T0, the transcriptome of SA24T0 was compared with that of the control quarters from the 72 hours S. aureus infection (SA72T0). No clones exhibited differential expression using a FDR less than 0.05. However, using less stringent statistical analyses revealed 115 clones that exhibited differential expression (P value ≤ 0.01, fold difference ≥ 2), which are listed in Additional file 2. Eighty one of these clones are annotated, representing 76 genes. Twenty two annotated genes (24 clones) were expressed at significantly higher levels in SA72T0. The remaining 54 annotated genes (57 clones), over two thirds of the differentially expressed genes, were expressed at significantly higher levels in SA24T0. The clone exhibiting the highest expression level compared to SA72T0 represents the cell surface marker CD4, with 4.5 fold higher expression. This suggests that the cellular composition of the control udder quarters differed depending on the time of infection of neighbouring quarters. SOD2, which was identified in the earlier analysis as an up-regulated gene during S. aureus infection, was expressed 2.3 fold higher in SA24T0 than SA72T0. Similarly, other genes associated with the mammary gland response to infection are expressed at higher levels in SA24T0, e.g. S100A8, S100A12, CXCL10, interleukin (IL) 1B and lactotransferrin (LTF).
Gene enrichment analysis using DAVID revealed that several GO Biological Process terms are over-represented in the DEG list (Table 3). The GO term exhibiting the greatest enrichment was Response to Stress. Other over-represented GO terms include; Cell Proliferation, Cell Development, Immune Response and Defence Response. The majority of genes associated with these GO terms are expressed at higher levels in SA24T0 samples. Furthermore, subsequent IPA analysis revealed that the DEG list contained 21 genes that are involved in the functional network with the pro-inflammatory cytokine IL1 at its hub and 17 of these were expressed at higher levels in SA24T0 than SA72T0 samples (Figure 4). These results suggest that uninfected quarters respond to the infection of neighbouring quarters with S. aureus for at least 24 hours, but this response had waned by 72 hours post infection. Therefore, of the three sets of control samples generated in this study the transcriptome of SA72T0 samples most closely resembles that of an uninfected quarter. This explains the low number of DEG identified in the initial analysis of the transcriptional response of infected quarters to infection (Table 1). The EC24T0 and SA24T0 samples were used as the baseline to identify DEG, but the response of these quarters to the infection of neighbouring quarters was sufficient to mask the transcriptional response of the infected quarters.
Validation of microarray results comparing SA24T0 and SA72T0
The validity of the DEG identified from the comparison of SA24T0 and SA72T0 using less stringent statistical analysis was investigated by RT-qPCR analysis of the expression levels of six genes in the RNA samples. Analysis of the microarray results suggested that CD4 was the most highly expressed gene in SA24T0 compared to SA72T0. However, the RT-qPCR revealed the opposite pattern of expression, with 2.1 fold higher levels in SA72T0 samples (P < 0.005) (Figure 5A). Furthermore, mannosyl (alpha-1,3-)-glycoprotein beta-1,2-N-acetylglucosaminyltransferase (MGAT1), which according to the microarray analysis was the most highly expressed gene in SA72T0 compared to SA24T0, exhibited no differential expression by RT-qPCR analysis (Figure 5A). However, frequently rearranged in advanced T-cell lymphomas (FRAT1) and LTF exhibited the same pattern of expression by RT-qPCR and microarray analyses. FRAT1 was expressed at on average 2.9 fold higher levels in SA72T0 samples (P < 0.05), whilst LTF was expressed at on average 2.4 fold higher levels in SA24T0 samples (P < 0.05) (Figure 5A). Furthermore, although S100A8 and S100A12 did not exhibit statistically significant differential expression by RT-qPCR analysis, the results followed the same trend as observed with the microarray results, with 2.7 fold and 3.6 fold higher expression in SA24T0 samples respectively (Figure 5A). The disparity between the microarray and RT-qPCR analyses may relate to the use of different RNA samples from the same udder quarters suggesting that there are localized differences in the response to infection of neighbouring quarters. Alternatively, the disparity may be due to the small number of out-bred animals used in this study.
In addition, the gene expression of the ten genes used to validate the microarray results for the comparison of EC24T0 and SA24T0 were also measured in the SA72T0 samples. No significant difference in expression was observed in the five genes more highly expressed in EC24T0 samples than SA24T0 samples (data not shown). However, four of the five genes more highly expressed in SA24T0 than EC24T0 samples were found to be more highly expressed in SA72T0 than SA24T0 samples and three showed statistical significance (Figure 5B). None of these four genes were identified from the microarray analysis comparison of SA24T0 and SA72T0. Levels of ABCG2, CXCL12 and IDH1 were all statistically significantly higher in SA72T0 than SA24T0 samples, with on average 1.9, 2.3 and 2.0 fold greater expression. Although not statistically significant XDH was on average 2.2 fold higher in SA72T0 compared to SA24T0 samples (P = 0.053). In contrast, TGFB1 levels were higher in SA24T0 samples compared to SA72T0 samples, albeit not to a statistically significant level. Therefore, although not in complete agreement with the microarray data, the RT-qPCR results confirm that the transcriptional profiles of uninfected mammary gland quarters differ during the course of S. aureus infection of neighbouring quarters.
Studies involving farm animals require the use of a large number of animals to obtain statistically meaningful data because these out-bred populations exhibit considerable genetic variation. Therefore, it is appealing to use intra-animal controls where possible, since this reduces genetic variation and fewer biological replicates are required. The quarters of the bovine mammary gland have been considered as separate physiological entities and therefore the use of one or more quarters as controls is a generally accepted experimental design [3–7, 11–16]. However, the results reported here demonstrate that the mammary gland quarters are interdependent during infection, with a profound change in the transcriptome of one quarter when neighbouring quarters are infected with S. aureus and especially E. coli. The control quarters remained bacteriologically negative throughout the study period  and therefore the observed differences are not due to contamination. The transcriptional response of the uninfected quarters was sufficiently high that it resulted in the masking of the response in infected quarters, when they were used as the base-line for DEG analysis (Table 1). In contrast, previous studies investigating the bovine mammary gland response to infection which have used within-animal controls have identified large numbers of differentially expressed genes [4–7]. However, in the earlier studies the animals were only challenged once, whilst in this study the animals were inoculated three times over 24 hours, which may have resulted in a greater signal triggering the response in the uninfected quarters.
Although a direct comparison between uninfected quarters from healthy and uninfected animals was not made in this study, the effect of infection in neighbouring quarters can be inferred from the differences observed between the various control samples. Recently, a microarray experiment comparing the transcriptome of uninfected quarters from healthy and E. coli infected animals was reported . A total of 476 DEG were identified, confirming that E. coli infection alters the transcriptional response of neighbouring mammary gland quarters. Over twenty eight percent (53) of the DEG identified in our comparison of E. coli and S. aureus control quarters were also found in this study, which utilized the bovine Affymetrix platform, including ABCG2 and FKBP5. The direction of the transcriptional change observed in these 53 genes was the same in both studies, i.e. those up-regulated during E. coli infection of neighbouring quarters  were found at higher levels in EC24T0 samples than SA24T0 samples, suggesting that the differential expression of at least these 53 genes is associated with E. coli infection and not S. aureus infection. The broader intensity plots of the EC24T0 samples illustrated that there was a greater range of transcriptional change in these samples compared to the reference than observed in the SAT0 samples. This may suggest that the vast majority of the differential expression is due to the E. coli infection of neighbouring quarters. Furthermore, the transcriptome of the uninfected quarters was indistinguishable from that of quarters infected with E. coli for six hours.
However, the design of the experiment reported here does not allow us to conclusively say that the transcriptional differences are all due to changes in expression in the E. coli control quarters. Indeed, the comparison of control quarters from animals infected with S. aureus for 24 and 72 hours has revealed that S. aureus also induces transcriptional changes in neighbouring quarters 24 hours post infection, which may account for some of the differences observed between SA24T0 and EC24T0. Furthermore, XDH, which is involved in the generation of the bactericidal agent peroxynitrate and is believed to play an important role in the mammary gland immune response (reviewed by ), is expressed at significantly higher levels in SA24T0 than EC24T0 and therefore may be induced during S. aureus infection of neighbouring quarters. Moreover, additional evidence that S. aureus infection induces a response in neighbouring quarters has recently been reported. The expression of the acute phase proteins serum amyloid A3 (SAA3) and haptoglobin (HP) as well as the β-defensin DEFB5 was elevated in uninfected control quarters from animals infected with a low virulence strain of S. aureus compared to tissue from healthy animals .
The initial statistical analysis of the microarray data comparing the S. aureus control quarters from the different time courses failed to detect any significantly differentially expressed genes. Although this could be a genuine result, S. aureus infection generally induces a milder immune response than E. coli, modulating gene expression to a lesser extent at the site of infection, which would impact on the signals being detected by the neighbouring quarters. This variation, compounded with the well-known fact that individual cattle differ greatly in their response to udder infections (reviewed by ), would result in considerable variability in the response of neighbouring quarters to S. aureus infection and could account for the lack of statistical significance, which would be further exacerbated by the low number of biological replicates included in this study. Indeed, lowering the stringency of the microarray analysis allowed the identification of over 70 genes, including several genes known to be associated with the mammary gland response to infection, e.g. SOD2, IL1B, LTF, S100A8 and S100A12. Nearly one sixth of the differentially expressed genes are associated with the response of cells to stress. Furthermore, IPA revealed that a network of IL1 and associated molecules were expressed at higher levels in SA24T0 samples than in the SA72T0 samples, which are a well-known trigger for inflammation. Paradoxically, the cattle infected with S. aureus did not show signs of inflammation, even in the infused udder quarters, 24 hours post infection , when the elevated levels of IL1 and related molecules were detected in this study. Conversely, at 72 hours, when the IL1 response was lower in uninfected quarters, the infected quarters showed up-regulation of inflammation-associated genes . Therefore, S. aureus infection may initially elicit a slight systemic inflammatory response that induces IL1 expression in uninfected quarters, which disappears during sustained infection. It is unclear if S. aureus actively contributes to the down-regulation of the systemic response, which may contribute to the typically sub-clinical disease pattern observed in S. aureus mastitis.
Some of the transcriptional changes may relate to changes in the cellular composition of the mammary gland. Bacterial infections induce the influx of immune cells, predominantly neutrophils, into the infected mammary gland quarter, which leads to increasing SCC in milk. The SCC measured in the control quarters from animals included in this study increased during the time course of both E. coli and S. aureus infection . This is in agreement with previous studies which found significantly higher SCC in milk from uninfected quarters of naturally infected animals compared to those from healthy animals [25–27]. Moreover, the percentage of neutrophils was significantly higher, whilst the percentage of macrophages was significantly lower, in milk from uninfected quarters from animals with mastitis compared to those from healthy animals . Differential expression of cell surface markers, e.g. CD74, Fc receptors and CD4, are indicative of differences in the cellular composition of the control quarters from E. coli and S. aureus infected animals. Levels of CD4+ T lymphocytes increase in milk during S. aureus infection [26, 28] and the results reported here suggest that some cells may originate from the uninfected quarters.
The influx of immune cells may not account for all the transcriptional changes of uninfected quarters. The response of the first quarter to infection with E. coli was significantly different from that of the second infected quarter after the same period of infection. The profound increase in SCC observed in milk from the first quarter 12 hours post infection was not observed in milk samples from the second quarter infected 12 hours later . The earlier infection appears to have primed the neighbouring mammary gland quarters to respond differently to infection. Similar results have been reported recently in an experimental infection study where one animal showed signs of infection in one udder quarter before the experiment started . This animal was still included in the study, using the other, uninfected udder quarters. The proinflammatory cytokine response, i.e. up-regulation of IL1B, IL6, IL8 and tumour necrosis factor (TNF), of this animal was more pronounced than that observed in other animals included in the study . The priming of the mammary gland may persist after infection. Experimental E. coli infection of individual quarters fourteen days after infection of neighbouring quarters resulted in significantly less bacteria being isolated from the milk than observed with the initial infection . Furthermore, it was recently reported that priming udders with a low dose of E. coli lipopolysaccharide (LPS) protects against experimental infection with E. coli pathogens . Moreover, LPS priming induces sustained expression of bactericidal factors (e.g. β-defensins) by bovine mammary epithelial cells (MEC), whilst reducing the over expression of pro-inflammatory cytokines which may be harmful to the animal . However, there is little overlap in DEG from primed MEC and those observed in the uninfected quarters included in this study, suggesting that different mechanisms may be involved in priming neighbouring quarters.
TGFB1 levels were observed to be much lower in the control quarters from E. coli infected animals than those from animals infected with S. aureus. Therefore, it is interesting to hypothesize that the priming of the udder quarters involves the suppression of levels of TGFB1 and associated molecules which was revealed by IPA. Moreover, the same pattern of TGFB3 expression has previously been reported using the same experimental samples . The TGFBs are pleiotropic cytokines which play an important role in multiple biological processes including regulating the adaptive and innate immune responses. TGFB1 inhibits inflammation by, amongst other mechanisms, suppressing the activation of several immune cells, e.g. NK cells and macrophages, and although chemotactic, can inhibit neutrophil transmigration (reviewed by ). Further analysis has shown that significantly lower levels of TGFB1 mRNA is present in the uninfected quarters of E. coli infected animals included in this study than in the quarters of healthy animals involved in another study  (Additional file 3). It appears that animals down regulate the expression of TGFB1 in the udder when a neighbouring quarter has been challenged with E. coli pathogens. This reduction in the level of the immunosuppressor TGFB1 would lower the threshold for triggering the induction of a variety of immunologically relevant genes, resulting in a faster and more robust response to infection. In contrast, the weaker stimulus elicited by invading S. aureus pathogens maintains, or possibly enhances (Additional file 3), the relatively high levels of TGFB1 present in uninfected quarters, thereby buffering immune stimuli in the healthy, uninfected udder quarter.
Response to Stress was one of the major GO Biological Processes categories identified in both the comparisons of EC24T0 and SA24T0 and that of S. aureus control quarters. Many of the genes within this category have previously been identified as being involved in the immune response to mastitis, including LBP, STAT3, S100A8 and SOD2. However, many of the most differentially regulated genes, e.g. ABCG2 and FKBP5, have not, to our knowledge, previously been associated with mastitis. Both ABCG2 and FKBP5 were also identified as DEG in the recent transcriptional comparison of uninfected quarters from healthy and mastitic animals , suggesting that their differential expression was due to the infection of neighbouring quarters with E. coli. Furthermore, the expression of both genes was modulated in E. coli infected mammary gland quarters .
ABCG2, also known as breast cancer resistance gene (BCRP), is a xenobiotic transporter that is expressed at high levels in lactating mammary glands of cows, humans and mice [34, 35]. ABCG2 plays an important role in the transport of nutrients into milk, especially riboflavin (vitamin B2)  and has been identified as the candidate gene underlying a quantitative trait locus on bovine chromosome 6 that affects milk yield and composition [37, 38]. Therefore, the down-regulation of ABCG2 in quarters neighbouring E. coli infected quarters may be associated with the decrease in milk yield and casein synthesis, which was observed in these animals [10, 39] and may imply a nutrient deficiency in milk from these animals.
The immunophilin FKBP5 has several functions, including mediating the action of the immunosuppressive drug FK506. The FKBP5-FK506 complex binds protein phosphatase 3, catalytic subunit, alpha isozyme (PPP3CA), also known as calcineurin, thus preventing it from phosphorylating and activating NFAT transcription factors, which in turn prevents T cell activation and proliferation . FKBP5 is highly expressed in murine T cells  and therefore the up-regulation observed during E. coli infection may be due to lymphocyte migration into mammary gland tissue. In addition, FKBP5 is induced by glucocorticoids and is believed to be involved in the negative feedback loop for glucocorticoid receptor signalling . FKBP5, along with heat shock protein 90 (HSP90), is a component of the glucocortcoid receptor complex. The HSP90 family member HSP90AB1 exhibits a similar expression profile as FKBP5 in both this study and Mitterheumer et al. . Glucocorticoids are able to modulate NF-κB activity by a FKBP5-dependent pathway . FKBP5 interacts with conserved helix-loop-helix ubiquitous kinase (IKKα), which phosphorylates nuclear factor of kappa light polypeptide gene enhancer in B-cells inhibitor, alpha (IκBα), the inhibitor of NF-κB, resulting in IκBα proteosome-mediated degradation and the nuclear translocation and activation of NF-κB .
The mechanism by which the infection of neighbouring quarters affects the transcriptome of uninfected quarters is unknown. Mastitis induces a systemic effect, inducing transcriptional changes in the liver , including the production of acute phase proteins which are released into the blood that circulates around all the udder quarters. This has been postulated to explain the transcriptional profile in healthy quarters of animals infected with E. coli. Furthermore, it has become clear that endogenous damage or danger signals  also lead to inflammation  and evidence now supports their role in the response to pathogens . Therefore, it is possible that the uninfected udder quarters are responding to such danger or damage signals released into the extracellular environment following infection in the neighbouring quarters. However, analysis of TLR and β-defensin expression in the local and peripheral lymph nodes of animals included in this study revealed that a response was only observed in the lymph nodes draining the quarters infected with E. coli for 12 and 24 hours , suggesting that the systemic effect is quite limited. Local cross-talk between the udder quarters must, therefore, be responsible for priming the neighbouring udder quarters. Indeed, lymphocytes have been shown to be capable of migrating between mammary gland quarters  and therefore the presence of a more localized interaction between quarters is probable and requires further investigation.
It has previously been shown that during E. coli infection the immune status of neighbouring, sterile quarters is modified . The results presented here provide further evidence of this and also illustrate that a similar phenomenon occurs during S. aureus infection. However, the details of the transcriptional response induced by these pathogens appear to be quite distinct, but are difficult to separate out fully in this study due to the experimental design used. However, the data suggest that the induced transcriptional changes prime the uninfected quarters to prevent or limit the spread of infection. The outcome of mastitis infections is influenced by multiple factors, including the virulence of the pathogen and host genetics. It is clear that the previous exposure of animals to pathogens is also important. Considering the frequency of mastitic infections on dairy farms the pathogen primed state may be the predominant status of uninfected quarters and therefore worthy of further investigation. Furthermore, the results described here illustrate that the four quarters of the bovine mammary gland cannot be considered as separate experimental units during infection studies. This has important ramifications for the design of experimental challenge studies, as the use of inappropriate controls can result in the masking of transcriptional changes during infection. Therefore caution should be taken when using within-animal controls.
Experimental design of the infection study
The design of the infection study has been described in detail previously in accordance with ARRIVE guidelines . This work was carried out under the approval of the ethics committee of the regional government in Hannover, Germany (No. 509.6-42502-03/678). Briefly, twelve Holstein cows in the middle of their first lactation were challenged over 24 or 72 hours with 500 cfu E. coli strain 1303 or 10,000 cfu S. aureus strain 1027. The animals were randomly divided into three groups. Four animals were inoculated intramammary with E. coli in the front right, hind right and hind left udder quarters at 0 (T24), 12 (T12) and 18 (T6) hours after the trial started respectively. The front left quarter received 2 ml sterile saline at the start of the experiment. The animals were euthanized after 24 hours, with a penetrating captive bolt gun followed by exsanguination, and tissue samples were collected asceptically from each udder quarter. The second group of four animals were inoculated with S. aureus over the same 24 hour period. The third group of four animals were inoculated with S. aureus in both right quarters (T72) and inoculated with 2 ml sterile saline in the front left quarter at 0 hours. The hind left quarter was inoculated at 60 hours (T12) after the start of the experiment. The cows were euthanized after 72 hours as described above and tissue samples were collected from each udder quarter. The control quarters from animals in these three groups are referred to as EC24T0, SA24T0 and SA72T0 respectively in this paper. All animals infected with E. coli developed acute mastitis within 12 hours, although the major clinical symptoms, e.g. udder swelling and elevated SCC, where only observed in the inoculated mammary gland quarters . In contrast, clinical signs of mastitis were only observed in animals infected with S. aureus for 72 hours .
A 5 cm3 tissue sample were removed asceptically from the inner area of each udder quarter, 10 cm dorsal of the milk cistern, within 5–10 minutes of culling. A small section of the centre of this sample (1.0 x 0.5 x 0.5 cm) was stored in RNAlater (Ambion) for RNA isolation. Total RNA was extracted as described previously  by caesium chloride gradient followed by phenol/chloroform extraction and sodium acetate precipitation. One hundred micrograms of each RNA sample was purified further using RNeasy mini columns (Qiagen) using the manufacturer’s protocols. The quality and quantity of the resulting RNA was determined by gel electrophoresis and Agilent 2100 Bioanalyzer.
The RNA samples from each mammary gland quarter were analyzed separately and hybridized in competition with a pooled reference sample, made up of all the RNA samples. A common reference design was used to allow the multifactorial analysis of the complete experiment. The reference sample was labelled with cyanine (Cy) 3 and the treatment sample with Cy5 on each microarray slide. The microarray experiment therefore consisted of a total of 48 microarray slides. The microarray data has been submitted to ArrayExpress and assigned the accession number E-TABM-704.
RNA labelling & hybridizations
Fluorescently labelled cDNA was generated from 20 μg total RNA using the FairPlay II microarray labelling kit (Stratagene), following the manufacturer’s protocol. The labelled cDNA was purified using a DyeEx spin column (Qiagen) and the labelling efficiency was determined by running 0.5 μl of each sample on a 1% agarose gel. Labelled cDNA was hybridized onto the ARK-Genomics Bos taurus 20 K v1.0 microarray (ArrayExpress accession no. A-MEXP-1402). Hybridizations were carried out in a GeneTac automated hybridization station (Genomic Solutions) to control for reactions conditions and washing of the arrays. The Cy3 and Cy5 labelled cDNA were mixed and added to 125 μl hybridization solution (ARK-Genomics protocols: http://www.ark-genomics.org) and hybridized onto the microarray for 12 hours. The microarray slides were then washed in wash buffers of increasing stringency (Genomic Solutions). After removal from the hybridization station the microarray slides were sequentially washed in post-wash buffer and isopropanol for one minute before being dried by centrifugation at 220 g for six minutes. The dried slides were scanned in a Scanarray 5000 XL scanner (GSI Lumonics) at constant laser power of 75% and 73% for Cy3 and Cy5 respectively.
Microarray spot intensity and quality data were extracted from the scanned images using the BlueFuse software version II (BlueGnome). Preliminary analysis was carried out using the GeneSpring software package (Agilent). Further analysis was carried out using modifications of the Limma package of Bioconductor , with additional plotting from the Marray package. The microarray data exhibited low signal to noise ratios, which resulted in excessive variability at low intensity values. To reduce this each channel was augmented by 256. The normalization and analysis steps were weighted using the BlueFuse quality measure, confidence values between 0 and 1 , and manually excluded spots were given the minimum confidence value. Each slide was normalized separately, using the log2-ratios of treatment to reference intensities for all the non-control spots. The data was normalized using a 2-step process of spatial followed by intensity dependent bias correction for each print-tip, as described previously . The spatial bias was corrected by subtracting corresponding row and column means from each spot . The intensity dependent bias was removed by block-lowess .
Genes that were differentially expressed between the infected and control samples were identified by analyzing the means of replicate spots. The Limma eBayes correction  was used to shrink the residual variances of each gene towards their median value. Comparisons between the infected and control samples were assessed by 1-way ANOVA followed by t-tests. FDR values for each comparison were calculated using the method of Benjamini and Hochberg . Genes with FDR values less than or equal to 0.05 were considered as significantly different. To refine the list of differentially expressed genes further to those of probable biological relevance, only those genes exhibiting a 2 fold or greater average difference were considered.
Validation of microarray experiments
RT-qPCR analysis was used to validate the microarray results for sixteen differentially expressed genes. Oligonucleotides were designed for each gene using Primer3  and Netprimer (Biosoft International) software (Table 4). Total RNA was isolated from a different sample of mammary gland tissue than that used for the microarray experiment. The tissue was powdered under liquid nitrogen and mixed with the chaotropic RLT lysis buffer (Qiagen). RNA was extracted with Phenol (Ambion) and purified using RNeasy spin columns (Qiagen) including a DNase digestion step (Qiagen) according to the manufacturer’s instructions. The RNA quality was assessed by agarose gel electrophoresis and NanoDrop ND-1000 spectrophotometer (Thermo Scientific).
The RT-qPCRs for twelve genes were carried out as described previously . The relative quantities of mRNA were calculated using the method described by Pfaffl . The RT-qPCR results for chromosome alignment maintaining phosphoprotein 1 (CHAMP1) were used to calculate differences in the template RNA levels and thereby standardize the results for the genes of interest. The RT-qPCR for the remaining four genes; chemokine (C-X-C motif) ligand (CXCL) 2, lipocalin 2 (LCN2), S100 calcium bind proteins (S100) A8 and S100A12, were carried out as described previously  and relative copy numbers were calculated. These results were converted to fold differences to standardize the results for all sixteen genes. All samples were measured twice from two independent cDNA preparations. The quality and authenticity of the resulting RT-qPCR products were assessed by agarose gel electrophoresis and sequencing. The relative quantity values were transformed on the log2 scale before statistical analyses to stabilize the variance and make them comparable to the log (intensity ratios) from the microarray analysis. The effect of infection in neighbouring quarters was examined by t-test analysis (Genstat 8.1, Lawes Agricultural Trust, Rothamsted).
Viguier C, Arora S, Gilmartin N, Welbeck K, O’Kennedy R: Mastitis detection: current trends and future perspectives. Trends Biotech. 2009, 27: 486-493.
Schukken YH, Günther J, Fitzpatrick J, Fontaine MC, Goetze L, Holst O, Leigh J, Petzl W, Schuberth HJ, Sipka A, Smith DGE, Quesnell R, Watts J, Yancey R, Zerbe H, Gurjar A, Zadoks RN, Seyfert HM: Host-response patterns of intramammary infections in dairy cows. Vet Immunol Immunopath. 2011, 144: 270-289.
Rinaldi M, Li RW, Bannerman DD, Daniels KM, Evock-Clover C, Silva MVB, Paape MJ, Van Ryssen B, Burvenich C, Capuco AV: A sentinel function for teat tissues in dairy cows: dominant innate immune response elements define early response to E. coli mastitis. Funct Integr Genomics. 2010, 10: 21-38.
Buitenhuis B, Røntved CM, Edwards SM, Ingvartsen KL, Sørensen P: In depth analysis of genes and pathways of the mammary gland involved in the pathogenesis of bovine Escherichia coli-mastitis. BMC Genomics. 2011, 12: 130-
Moyes KM, Drackley JK, Morin DE, Bionaz M, Rodriguez-Zas SL, Everts RE, Lewin HA, Loor JJ: Gene network and pathway analysis of bovine mammary tissue challenged with Streptococcus uberis reveals induction of cell proliferation and inhibition of PPARγ signalling as potential mechanism for the negative relationships between immune response and lipid metabolism. BMC Genomics. 2009, 10: 542-
Swanson KM, Stelwagen K, Dobson J, Henderson HV, Davis SR, Farr VC, Singh K: Transcriptome profiling of Streptococcus uberis-induced mastitis reveals fundamental differences between immune gene expression in the mammary gland and in a primary cell culture model. J Dairy Sci. 2009, 92: 117-129.
Lutzow YCS, Donaldson L, Gray CP, Vuocolo T, Pearson RD, Reverter A, Byrne KA, Sheehy PA, Windon R, Tellam RL: Identification of immune genes and proteins involved in the response of bovine mammary tissue to Staphylococcus aureus infection. BMC Vet Res. 2008, 4: 18-
Genini S, Badaoui B, Sclep G, Bishop SC, Waddington D, van der Pinard LMH, Klopp C, Cabau C, Seyfert HM, Petzl W, Jensen K, Glass EJ, de Greeff A, Smith HE, Smits MA, Olsaker I, Boman GM, Pisoni G, Moroni P, Castiglioni B, Cremonesi P, Del Corvo M, Foulon E, Foucras G, Rupp R, Giuffra E: Strengthening insights into host reponses to mastitis infection in ruminants by combining heterogeneous microarray data sources. BMC Genomics. 2011, 12: 225-
Loor JJ, Moyes KM, Bionaz M: Functional adaptations of the transcriptome to mastitis-causing pathogens: the mammary gland and beyond. J Mammary Gland Biol Neoplasia. 2011, 16: 305-322.
Petzl W, Zerbe H, Günther J, Yang W, Seyfert H-M, Nürnberg G, Schuberth HJ: Escherichia coli, but not Staphylococcus aureus, triggers an early increased expression of factors contributing to the innate immune defense in the udder of the cow. Vet Res. 2008, 39: 18-
Chang CC, Winter AJ, Norcross NL: Immune response in the bovine mammary gland after intestinal, local and systemic immunization. Infect Immun. 1981, 31: 650-659.
Mehrzad J, Dosogne H, Meyer E, Burvenich C: Local and systemic effects of endotoxin mastitis on the chemiluminescence of milk and blood neutrophils in dairy cows. Vet Res. 2001, 32: 131-144.
Bannerman DD, Paape MJ, Lee JW, Zhao X, Hope JC, Rainard P: Escherichia coli and Staphylococcus aureus elicit differential innate immune responses following intramammary infection. Clin Diagn Lab Immunol. 2004, 11: 463-472.
Schmitz S, Pfaffl MW, Meyer HHD, Bruckmaier RM: Short-term changes of mRNA expression of various inflammatory factors and milk proteins in mammary tissue during LPS-induced mastitis. Domest Anim Endocrin. 2004, 26: 111-126.
Silva LFP, Liesman JS, Etchbarne BE, Nielsen MSW, VandeHaar MJ: Intramammary infusion of IGF-1 increases bromodeoxyuridine labelling in mammary epithelial cells of prepubertal heifers. J Dairy Sci. 2005, 88: 2771-2773.
Grönlund U, Johannisson A, Waller KP: Changes in blood and milk lymphocyte sub-populations during acute and chronic phases of Staphylococcus aureus induced bovine mastitis. Res Vet Sci. 2006, 80: 147-154.
Mitterhuemer S, Petzl W, Krebs S, Mehne D, Klanner A, Wolf E, Zerbe H, Blum H: Escherichia coli infection induces distinct local and systemic transcriptome responses in the mammary gland. BMC Genomics. 2010, 11: 138-
Donaldson L, Vuocolo T, Gray C, Strandberg Y, Reverter A, McWilliam S, Wang Y, Byrne K, Tellam R: Construction and validation of a bovine innate immune microarray. BMC Genomics. 2005, 6: 135-
Dennis G, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, Lempicki RA: DAVID: database for annotation, visualization and integrated discovery. Genome Biol. 2003, 4: P3-
Huang DW, Sherman BT, Lempicki RA: Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009, 4: 44-57.
Günther J, Esch K, Poschadel N, Petzl W, Zerbe H, Mitterhuemer S, Blum H, Seyfert HM: Comparative kinetics of Escherichia coli- and Staphylococcus aureus-specific activation of key immune pathways in mammary epithelial cells demonstrates that S. aureus elicits a delayed response dominated by interleukin-6 (IL-6) but not IL-1A or tumor necrosis factor alpha. Infect Immun. 2011, 79: 695-707.
Rainard P, Riollet C: Innate immunity of the bovine mammary gland. Vet Res. 2006, 37: 369-400.
Whelehan CJ, Meade KG, Eckersall PD, Young FJ, O’Farrelly CO: Experimental Staphylococcus aureus infection of the mammary gland induces region-specific changes in innate immune gene expression. Vet Immunol Immunopathol. 2011, 140: 181-189.
Burvenich C, Van Merris V, Mehrzad J, Diez-Fraile A, Duchateau L: Severity of E. coli mastitis is mainly determined by cow factors. Vet Res. 2003, 34: 521-564.
Bansal BK, Hamann J, Grabowski NT, Singh KB: Variation in the composition of selected milk fraction samples from healthy and mastitic quarters and its significance for mastitis diagnosis. J Dairy Res. 2005, 72: 144-152.
Hamann J, Schröder A, Merle R: Differential cell count and interdependence of udder quarters. Mastitis in dairy production: current knowledge and future solutions. Edited by: Hogeveen H. 2005, The Netherlands: Wageningen Academic Publishers, 190-195.
Merle R, Schröder A, Hamann J: Cell function in the bovine mammary gland: a preliminary study on interdependence of healthy and infected udder quarters. J Dairy Res. 2007, 74: 174-179.
Riollet C, Rainard P, Poutrel B: Cell subpopulations and cytokine expression in cow milk in response to chronic Staphylococcus aureus infection. J Dairy Sci. 2001, 84: 1077-1084.
Suojala L, Orro T, Järvinen H, Saatsi J, Pyörälä S: Acute phase response in two consecutive experimentally induced E. coli intramammary infections in dairy cows. Acta Vet Scand. 2008, 50: 18-
Petzl W, Günther J, Pfister T, Sauter-Louis C, Goetze L, von Aulock S, Hafner-Marx A, Schuberth H-J, Seyfert H-M, Zerbe H: Lipopolysaccharide pretreatment of the udder protects against experimental Escherichia coli mastitis. Innate Immun. 2012, 18: 467-477.
Günther J, Petzl W, Zerbe H, Schuberth H-J, Koczan D, Goetze L, Seyfert H-M: Lipopolysaccharide priming enhances expression of effectors of immune defence while decreasing expression of pro-inflammatory cytokines in mammary epithelia cells from cows. BMC Genomics. 2012, 13: 17-
Yang W, Zerbe H, Petzl W, Brunner RM, Gunther J, Draing C, von Aulock S, Schuberth HJ, Seyfert HM: Bovine TLR2 and TLR4 properly transduce signals from Staphylococcus aureus and E. coli, but S. aureus fails to both activate NF-κB in mammary epithelial cells and to quickly induce TNFα and interleukin-8 (CXCL8) expression in the udder. Mol Immunol. 2008, 45: 1385-1397.
Li MO, Wan YY, Sanjabi S, Robertson AKL, Flavell RA: Transforming growth factor-β regulation of immune responses. Annu Rev Immunol. 2006, 24: 99-146.
Jonker JW, Merino G, Musters S, van Herwaarden AE, Bolscher E, Wagenaar E, Mesman E, Dale TC, Schinkel AH: The breast cancer resistance protein BCRP (ABCG2) concentrates drugs and carcinogenic xenotoxins into milk. Nat Med. 2005, 11: 127-129.
Mani O, Sorensen MT, Sejrsen K, Bruckmaier RM, Albrecht C: Differential expression and localization of lipid transporters in the bovine mammary gland during the pregnancy-lactation cycle. J Dairy Sci. 2009, 92: 3744-3756.
van Herwaarden AE, Wagenaar E, Merino G, Jonker JW, Rosing H, Beijnen JH, Schinkel AH: Multidrug transporter ABCG2/breast cancer resistance protein secretes riboflavin (vitamin B2) into milk. Mol Cell Biol. 2007, 27: 1247-1253.
Cohen-Zinder M, Seroussi E, Larkin DM, Loor JJ, Loor JJ, van der Everts WA, Lee J-H, Drackley JK, Band MR, Hernandez AG, Shani M, Lewin HA, Weller JI, Ron M: Identification of a missense mutation in the bovine ABCG2 gene with a major effect on the QTL on chromosome 6 affecting milk yield and composition in Holstein cattle. Genome Res. 2005, 15: 936-944.
Weikard R, Widmann P, Buitkamp J, Emmerling R, Kuehn C: Revisiting the quantitative trait loci for milk production traits on BTA6. Anim Genet. 2011, 43: 318-323.
Vanselow J, Yang W, Herrmann J, Zerbe H, Schuberth H-J, Petzl W, Tomek W, Seyfert H-M: DNA-remethylation around a STAT5-binding enhancer in the αS1-casein promoter is associated with abrupt shutdown of αS1-casein synthesis during acute mastitis. J Mol Endocrinol. 2006, 37: 463-477.
Baughman G, Wiederrecht GJ, Campbell NF, Martin MM, Bourgeois S: FKBP51, a novel T-cell-specific immunophilin capable of calcineurin inhibition. Mol Cell Biol. 1995, 15: 4395-4402.
Denny WB, Valentine DL, Reynolds PD, Smith DF, Scammell JG: Squirrel monkey immunophilin FKBP51 is a potent inhibitor of glucocorticoid receptor binding. Endocrinology. 2000, 141: 4107-4113.
Park J, Kim M, Na G, Jeon I, Kwon Y, Kim J, Youn H, Koo Y: Glucocorticoids modulate NF-κB-dependent gene expression by up-regulating FKBP51 expression in Newcastle disease virus-infected chickens. Mol Cell Endo. 2007, 278: 7-17.
Bouwmeester T, Bauch A, Ruffner H, Angrand P-O, Bergamini G, Croughton K, Cruciat C, Eberhard D, Gagneur J, Ghidelli S, Hopf C, Huhse B, Mangano R, Michon AM, Schirle M, Schleg J, Schwab M, Stein MA, Bauer A, Casari G, Drewes G, Gavin AC, Jackson DB, Joberty G, Naubauer G, Rick J, Kuster B, Superti-Furga G: A physical and functional map of the human TNF-α/NF-κB signal transduction pathway. Nat Cell Biol. 2004, 6: 97-105.
Jiang L, Sørensen P, Røntved C, Vels L, Ingvartsen KL: Gene expression profiling of liver from dairy cows treated intra-mammary with lipopolysaccharide. BMC Genomics. 2008, 9: 443-
Matzinger P: Tolerance, danger, and the extended family. Annu Rev Immunol. 1994, 12: 991-1045.
Manson J, Thiemermann C, Brohi K: Trauma alarmins as activators of damage-induced inflammation. Brit J Surg. 2012, 99 (Suppl 1): 12-20.
Fontana MF, Vance RE: Two signal models in innate immunity. Immunol Rev. 2011, 243: 26-39.
Kimura K, Harp JA, Goff JP, Olsen SC: Lymphocytes from one side of the bovine mammary gland migrate to the contra lateral gland and lymph node tissue. Vet Immunol Immunopathol. 2005, 108: 409-415.
Chirgwin JM, Przybyla AE, MacDonald RJ, Rutter WJ: Isolation of biologically active ribonucleic acid from sources enriched in ribonuclease. Biochemistry. 1979, 18: 5294-5299.
Smyth GK: Bioinformatics and Computational Biology Solutions using R and Bioconductor. Bioinformatics and Computational Biology Solutions using R and Bioconductor. Edited by: Gentleman R, Carey V, Dudoit S, Irizarry R, Huber W. 2005, New York: Springer, 397-420.
Jaffrézic F, de Koning DJ, Boettcher PJ, Bonnet A, Buitenhuis B, Closset R, Déjean A, Delmas C, Detilleux JC, Dovč P, Duval M, Foulley J-L, Hedegaard J, Hornshøj H, Hulsegge I, Janss L, Jensen K, Jiang L, Lavrič M, Lê Cao K-A, Lund MS, Malinverni R, Marot G, Nie H, Petzl W, Pool MH, Robert-Granié C, San Cristobal M, van Schothorst EM, Schuberth H-J, et al: Analysis of the real EADGENE dataset: comparison of methods and guidelines for data normalisation and selection of differentially expressed genes. Genet Sel Evol. 2007, 39: 633-650.
Jensen K, Talbot R, Paxton E, Waddington D, Glass EJ: Development and validation of a bovine macrophage specific cDNA microarray. BMC Genomics. 2006, 7: 224-
Baird D, Johnstone P, Wilson T: Normalization of microarray data using a spatial mixed model analysis which includes splines. Bioinformatics. 2004, 20: 3196-3205.
Smyth GK, Speed TP: Normalization of cDNA microarray data. Methods. 2003, 31: 265-273.
Smyth GK: Linear models and empirical Bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004, 3: 3-
Benjamini Y, Hochberg Y: Controlling the false discovery rate: a practical and powerful approach to multiple testing. J Roy Stat Soc B. 1995, 57: 289-300.
Rozen S, Skaletsky HJ: Primer3 on the WWW for general users and for biologist programmers. Bioinformatics Methods and Protocols: Methods in Molecular Biology. Edited by: Krawetz S, Misener S. 2000, Totowa: Humana Press, 365-386.
Jensen K, Paxton E, Waddington D, Talbot R, Darghouth M, Glass EJ: Differences in the transcriptional responses induced by Theileria annulata infection in bovine monocytes derived from resistant and susceptible cattle breeds. Int J Parasitol. 2008, 38: 313-325.
Pfaffl MW: A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. 2001, 29: e45-
Goldammer T, Zerbe H, Molenaar A, Schuberth HJ, Brunner RM, Kata SR, Seyfert HM: Mastitis increases mammary mRNA abundance of beta-defensin 5, toll-like-receptor 2 (TLR2), and TLR4 but not TLR9 in cattle. Clin Diagn Lab Immunol. 2004, 11: 174-185.
This work was supported by the European Animal Disease Genomics Network of Excellence for Animal Health and Food Safety (EADGENE) (EU Contract No. FOOD-CT-2004-506416). EJG was also supported by the Biotechnology and Biological Sciences Research Council (Institute Strategic Programme Grant). HMS was also supported by the Deutsche Forschungsgemeinschaft, through the Researcher Group FOR585 (grant Se 326/14-3). We would like to thank Dave Waddington from The Roslin Institute & Royal (Dick) School of Veterinary Studies, University of Edinburgh for the statistical analysis of the microarray data. We would also like to thank Mark Fell from the ARK-Genomics facility for helping with the submission of the experimental data to ArrayExpress.
The authors declare that they have no competing interests.
KJ designed the microarray experiment, analyzed the differentially expressed genes and drafted the manuscript. RT constructed the microarray and oversaw the microarray experiments. KJ and JG carried out the RT-qPCR analysis. JG analyzed the gene lists using Ingenuity. HMS, WP, HZ and HJS designed and carried out the infection study and collected the samples. HMS conceived the study and coordinated the use of the infection samples in the microarray experiment. EJG was involved in the microarray design and coordination of the study. All authors read and approved the final manuscript.