Skip to main content
  • Research article
  • Open access
  • Published:

Escherichia coli infection induces distinct local and systemic transcriptome responses in the mammary gland

Abstract

Background

Coliform bacteria are the most common etiologic agents in severe mastitis of cows. Escherichia coli infections are mostly restricted to a single udder quarter whereas neighboring quarters stay clinically inapparent, implicating the presence of a systemic defense reaction. To address its underlying mechanism, we performed a transcriptome study of mammary tissue from udder quarters inoculated with E. coli (6 h and 24 h post infection), from neighboring quarters of the same animals, and from untreated control animals.

Results

After 6 h 13 probe sets of differentially expressed genes (DEG) were detected in infected quarters versus control animals. Eighteen hours later 2154 and 476 DEG were found in infected and in neighboring quarters vs. control animals. Cluster analysis revealed DEG found only in infected quarters (local response) and DEG detected in both infected and neighboring quarters (systemic response). The first group includes genes mainly involved in immune response and inflammation, while the systemic reaction comprises antigen processing and presentation, cytokines, protein degradation and apoptosis. Enhanced expression of antimicrobial genes (S100A8, S100A9, S100A12, CXCL2, GNLY), acute phase genes (LBP, SAA3, CP, BF, C6, C4BPA, IF), and indicators of oxidative stress (GPX3, MT1A, MT2A, SOD2) point to an active defense reaction in infected and neighboring healthy quarters. Its early onset is indicated by increased transcription of NFIL3 at 6 h. NFIL3 is a predicted regulator of many genes of the systemic response at 24 h. The significance of our transcriptome study was evidenced by some recent findings with candidate gene based approaches.

Conclusions

The discovery and holistic analysis of an extensive systemic reaction in the mammary gland significantly expands the knowledge of host-pathogen interactions in mastitis which may be relevant for the development of novel therapies and for genetic selection towards mastitis resistance.

Background

Udder infections cause considerable economic losses to the dairy industry [1] and often lead to culling of affected animals [2]. Due to the economic importance of mastitis and the health risk for consumers, large efforts have been made to identify factors involved in the susceptibility of dairy cows to infections of the mammary gland [3, 4] which occur most frequently at parturition [5]. In that period, infections with E. coli often cause severe clinical symptoms [6, 7] accompanied by reduction in milk yield, altered milk composition and extensive damage of mammary tissue [8]. E. coli infection predominantly affects a single udder quarter [9]. The other quarters are clinically asymptomatic and are considered to be healthy. Nevertheless, there are indications that neighboring apparently healthy udder quarters are infected without showing clinical changes. Epidemiological studies revealed, that the identical environmental strain of E. coli can cause recurrent mastitis, but in different udder quarters relative to the one infected first [10]. This phenomenon was attributed to initial infection of multiple quarters, each at a different level and with variations in lesion development [10].

The symptoms of a defined experimental mastitis vary between different animals from mild to lethal. Interestingly, the severity of clinical response seems to be relatively characteristic and constant for a given animal, as shown by repeated infection studies [11]. The most prominent early response of the mammary gland is a dramatic increase in somatic cell numbers in milk occurring after infusion of E. coli[12] or its lipopolysaccharide [13] solely in the infected udder quarter. Moreover the infection provokes systemic effects as e.g. evidenced by fever, increased level of TNF-alpha in blood [14] or enhanced expression of several genes like haptoglobin in neighboring quarters [15]. In addition, a systemic reaction is evidenced by altered gene expression in distant organs as shown for liver after infusion of lipopolysaccharide into the mammary gland [16]. Surprisingly, experimental infection of a second healthy udder quarter twelve hours after the primary infection does not change the somatic cell count in this quarter, while the increase in somatic cell numbers in the first infected quarter continues [12, 17].

These observations favor the concept, that infection with E. coli provokes two distinct types of response to the pathogen. The first comprises locally restricted reactions with acute symptoms and inflammation of the infected udder quarter. The second occurs in infected as well as in neighboring healthy udder quarters and furthermore influences even distant organs as shown for liver [16]. This systemic response may impair progression of subsequent infections of neighboring mammary quarters [12, 17]. The interplay of both local and systemic reactions may finally determine the course of infection in an individual animal. The elucidation of the underlying mechanisms may help to improve mastitis resistance of livestock populations by selective breeding and reveal new ways to treat infections with E. coli. Here we present a transcriptome profiling approach to study the impact of an experimental E. coli infection on spatial gene expression profiles of the bovine mammary gland and to distinguish between local and systemic effects. The experimental system comprises dairy cows with defined lactation and immune status as hosts and E. coli 1303 as pathogen. Microarray analyses of mammary gland tissue were performed in the early (6 h) and the late (24 h) phase of experimental infection. Udder tissue from healthy, untreated cows served as control for the discrimination between local and systemic effects in infected animals. We found, that infection of a single udder quarter with E. coli 1303 induced substantial transcriptome changes both in the infected and in the neighboring udder quarters with more than two thousand differentially expressed genes. Cluster analyses of DEG revealed that two different reactions took place. The first was locally restricted to the infected quarter whereas the second occurred systemically in infected as well as in their neighboring noninfected quarters. Both reactions were characterized by gene ontology analyses and by searching of putative regulatory elements in the promoter regions of the DEG. Our results offer reasonable hypotheses why, after initial infection of one udder quarter, massive expansion of the bacteria in neighboring quarters is suppressed.

Results

Establishment of intramammary infection

Prior to infection, the udder secretions of all animals had repeatedly been tested negative for major and minor mastitis pathogens. Following administration of 500 cfu of E. coli 1303 to a single udder quarter and the same volume of physiological saline into a neighboring quarter of each of five cows (Figure 1), all inoculated animals suffered from acute mastitis and showed a strong cellular response. Major clinical symptoms, such as udder swelling and elevated counts of somatic cells in milk, were restricted to pathogen-inoculated udder quarters (Figure 2A). All milk samples from the inoculated quarters were positive for E. coli bacteria whereas neighboring quarters and all quarters of control animals remained bacteriologically negative during the observation period.

Figure 1
figure 1

Experimental setup of infection with E. coli 1303. During estrus, 10 synchronized dairy cows were infected on one quarter with E. coli 1303 (star, i) in saline solution (shaded). One of the neighboring quarters (n) was treated with sterile saline solution (shaded). Five synchronized healthy animals were slaughtered 24 h after onset of estrus and served as external control (c). Tissue samples were collected from marked areas of the bovine udder.

Figure 2
figure 2

Results of clinical and paraclinical findings. Time course of A: somatic cell count in the milk of infected and neighboring noninfected udder quarters; B: rectal body temperature; and C: leukocyte counts in blood of infected and healthy control animals. D: milk yield of infected and neighboring noninfected (untreated or saline treated) udder quarters. The data shown are mean values and their corresponding standard errors and do not include values from animals of the 6 h infection group.

Body temperature, milk yield, differential blood cell counts

All infected cows developed fever (>39.2°C) within the first 14 h after infection with maximal values (41.1 ± 0.2°C) between 15 h and 18 h after inoculation (Figure 2B). This was accompanied by significant leukopenia, reaching minimal blood leukocyte counts 12 h after pathogen inoculation (Figure 2C). Milk yield dropped progressively within the 24-h period after infection (Figure 2D), with a stronger decrease in infected quarters than in neighboring quarters. In infected quarters the number of somatic cells in milk increased starting between 6 and 12 h after infection (Figure 2A). In mock inoculated and untreated neighboring udder quarters of infected animals no significant increase in the number of somatic cells was observed.

Microarray analysis

Mammary tissue samples were used to generate four sets of transcriptome profiles from: i) inoculated quarters 6 h post infection; ii) inoculated quarters 24 h post infection; iii) noninfected quarters neighboring the inoculated quarters 24 h post infection; and iv) quarters of untreated healthy control animals. Samples of sets iii and iv were generated by combining tissues of two udder quarters from each cow in order to equalize potential quarter specific effects. Based on somatic cell counts there was no indication for a different behavior of mock inoculated and untreated udder quarters of infected cows. Therefore tissues of two quarters with the lowest somatic cell counts were combined (see Additional File 1) provided their cell counts did not differ by more than 10.000 per ml. Otherwise tissue from the udder quarter with the lowest cell count was used. Thus four samples of set iii were pooled from one mock inoculated and one untreated udder quarter. The fifth sample (sample ID n24h-4) was derived from tissue of a single untreated quarter since the noninfected udder quarters of this cow differed in somatic cells counts.

The impact of the E. coli infection on gene expression was estimated by hierarchical clustering of the transcriptome profiles (Figure 3). The heatmap obtained with the transcriptome profiles 6 h after inoculation showed a uniform coloration with no separation of healthy controls and infected samples (Figure 3A). This indicated that the gene expression of the mammary tissue was only marginally affected in the early stage of infection. The same comparison 24 h after infection with E. coli resulted in a sharp separation of the infected and the control group (Figure 3B) and a distinct blue-red coloration indicated substantial transcriptome changes. At the same time point, the effect of E. coli on transcriptome changes in neighboring noninfected udder quarters was less pronounced (Figure 3C). Compared with Figure 3B, the length of the vertical branch between the two principal groups is considerable lower, indicating a less pronounced difference of the expression profiles.

Figure 3
figure 3

Hierarchical clustering of transcriptome profiles after experimental infection with E. coli. Tissues of udder quarters of five healthy untreated cows (c24h-1 to 5) were used to generate three replicates of expression profiles as controls for three separate comparisons. A: Replicate 1 (c124h-1 to 5) was compared to patterns of infected udder quarters 6 h after inoculation with E. coli 1303 (i6h-6 to 10). B: Replicate 2 (c224h-1 to 5) was compared to patterns of infected udder quarters 24 h after inoculation (i24h-1 to 5). C: Replicate 3 (c324h-1 to 5) was compared to patterns of noninfected quarters of infected cows 24 h after infection (n24h-1 to 5). Samples for n24h-1, 2, 3 and 5 were pooled from one mock inoculated and one untreated udder quarter, while n24h-4 was derived from tissue of a single untreated quarter. The degree of similarity of their expression patterns is shown as a heatmap of the distance matrix (red: more similar, blue: less similar). The trees depict the pair wise relations with branch lengths proportional to the degree of similarity.

Nevertheless two almost perfectly separated groups were obtained: three of four noninfected samples pooled from one mock inoculated and one untreated udder quarter (n24h-1, -2, -3) clustered together with the sample derived from a single untreated quarter (n24-4), while the fourth pooled sample (n24h-5) was placed into the well-separated control group.

DEG (see Additional File 2) were identified by comparison of transcriptome profiles of infected and neighboring quarters with the profiles of the control animals. The highest number of DEG was found 24 h post infection in infected quarters (2154 at 0.01 fdr), followed by 476 (0.01 fdr) in their neighboring noninfected quarters and only 13 (0.1 fdr) in infected quarters at 6 h after inoculation. Despite the noisy data at 6 h after infection 6 of the DEG were also detected as differentially expressed at 24 h in infected quarters and 2 in neighboring quarters. At 24 h post infection, 294 DEG of the infected udder quarters were found among the DEG of the neighboring noninfected quarters.

This overlap indicates either the existence of a systemic response to a local infection with E. coli or a local effect after application of 2 ml PBS in mock inoculated and infected udder quarters. The first reaction should take place in all udder quarters of infected cows, the second however should be absent in their untreated ones.

Local and systemic responses to E. coli revealed by cluster analysis

In order to look for local and systemic responses to E. coli infection, all DEG detected in infected quarters (2154) and noninfected neighboring quarters (476) were subjected to self organizing tree algorithm (SOTA) [18] clustering (Figure 4).

Figure 4
figure 4

Cluster analysis separates systemic and local response to infection. Expression levels of DEG 24 h after E. coli infection were converted to log2 fold changes relative to the mean of untreated control animals. The four samples labeled with i24h represent i24h-1 to 4 of infected udder quarters. c124h comprises the four samples c24h-1, 2, 4 and 5 derived from udder quarters of untreated controls. The four samples subscripted with n24h were isolated from udder quarters neighboring to the quarters i24h-1 to 4 and the four samples labeled with c224h are replicates of c124h. These technical replicates allow a common representation of the two separately performed analyses of i24h and n24h. Expression values were clustered using the SOTA algorithm with Pearson's correlation metric in the program MeV4.2. Individual probe sets are represented as lines, red coloration indicates up-regulation upon treatment, green down-regulation. Magenta represents the per sample average.

The calculated clusters revealed distinct reactions of the mammary gland in response to inoculation with E. coli. The first occurred exclusively in infected quarters (clusters 1 and 6) and comprised the local restricted response of the infected tissue. The second reaction took place both in infected and neighboring noninfected udder quarters (clusters 2 and 3) and pointed to a systemic response of the mammary gland to inoculation with E. coli in 2 ml PBS. Another type of reaction was recorded in the clusters 4 and 5. These 78 DGE were found induced specifically in the neighboring quarters, whereas in the infected quarters they were repressed (cluster 4) or unchanged (cluster 5). Their transcription was apparently regulated in a location specific manner. Nevertheless the different transcript levels in infected and neighboring noninfected udder quarters indicated that systemic signals are involved in their regulation.

The extent of apparently systemic gene regulation displayed in clusters 2 und 3 is much greater than the overlap of DEG from infected and neighboring quarters (294 genes). Roughly half of the 2335 genes fall into clusters that showed expression changes in infected as well as in neighboring quarters. This discrepancy is caused by the fact that DEG in infected quarters generally showed higher changes in expression than DEG in neighboring quarters. Due to the fold-change cut-off, genes with low fold-change are not assigned to the DEG but are nevertheless considered by SOTA clustering. Conversely, some genes differentially expressed in infected as well as in neighboring noninfected quarters were clustered by SOTA with the genes expressed locally in infected quarters (Figure 4, cluster 1), most likely because of their much higher induction in infected tissue than in tissue of neighboring quarters (e.g. S100A9, S100A12, SAA3, SDS, CP). For interpretation, those DEG, however, were treated as systemically reacting since the differential gene expression was regarded superior to clustering.

Based on annotations of the DEG more detailed analyses were performed. The genes of each cluster were analyzed in search of regulatory elements within their promoter region. Furthermore we identified biological processes involved in the systemic or local response of the mammary gland to the pathogen, which were represented by significantly more genes than expected by random sampling.

Both functional descriptions and transcription factor binding site (TFBS) enrichment analyses showed specific results for the members of the six clusters (Figure 4 and Table 1). Many results were specific for either local or systemic clusters and suggested that this dichotomy is not random but is biologically meaningful.

Table 1 Transcription factor binding site enrichment in clusters

Characterization of the local response 24 h after infection

Analysis of the enrichment of functional annotations revealed that the locally up-regulated genes of cluster 1 (Figure 4) have a distinct profile of functions that is not found among the members of the other clusters and that is dominated very prominently by ontology terms related to immune response and inflammation. In addition ontology terms like chemokine/cytokine signaling and diapedesis were among the most prominent. This correlated with our clinical findings of a sharp increase in leukocytes solely in the milk of infected udder quarters. Many genes of this cluster were associated with classical inflammatory events (Figure 5) and TFBS enrichment analysis (Table 1) predicted the classical inflammatory transcription factors NFκB and STAT1 as regulatory elements. Major genes related to fever induction, IL6 and IL1B, were still up-regulated (16-fold and 4-fold respectively) in infected udder quarters 24 h after infection although the inner body temperature was already declining and approached basal levels at this time point. Furthermore, the majority of observed acute phase genes (e.g. HP, PTX3, FTL) was assigned to cluster 1. As evidenced in Figure 5, almost all DEG of the ontology response to infection were found locally restricted to the infected quarters and only a few were observed also in neighboring quarters. The highly significant enrichment of biologically meaningful terms (e.g. adjusted p-value 1e-38 for 'immunity and defense') is also an important indication for the validity of the presented microarray data.

Figure 5
figure 5

Gene Ontology and pathway enrichment of clusters. The genes constituting the clusters in Figure 4 were subjected to enrichment analysis using the Pantherdb gene expression tool and the CoPub keyword enrichment calculator. Biological processes, pathways and molecular functions that were significantly enriched (adjusted p-value < 0.05) in at least one cluster are shown as a heatmap. P-values are color coded ranging from white for non-significant results to blue colors saturating with increasing significance. Rows indicate the enriched terms and columns the p-values of the respective clusters. Clusters 4 and 5 did not yield any significant enrichment and were omitted.

Besides the response to infection TFBS enrichment analysis gave hints towards other processes. Among the genes of cluster 1, ELF5 was predicted as involved transcription factor. ELF5 is related to development of mammary epithelium [19] and milk production [20].

Characterization of the apparent systemic response 24 h after infection

The transcriptional response in the quarters neighboring the infected ones could in part be characterized as an immune response. Components of the antigen processing and presenting machinery as well as some cytokines are up-regulated. Of note are CXCL2 and CXCL14, whose products show direct antimicrobial properties against E. coli[21] and could probably function as a barrier against invasion of E. coli into the neighboring udder quarters. DEG of other important immune-related agents, as for example the CC-type of chemotactic cytokines, were restricted to the infected udder quarter. A number of genes, like GPX3, GLRX2, GSTM1, GSS, MT1A, MT2A and SOD2, indicate the presence of oxidative stress. Differential gene expression related to xenobiotic stress is evident by up-regulation of ATF5, which is responsible for the induction of the detoxifying cytochrome CYP2B6 [22], the latter being up-regulated only in the neighboring quarters (cluster 2).

Finally, there are signs of building up a first line of defense, like the up-regulation of genes with antimicrobial functions (S100A8/9/12, GNLY, CXCL2, CXCL14) or acute phase genes (SAA3) involved in recognition of pathogens (LBP) or in the complement cascade (C6, BF, C4BPA, IF).

NFIL3 as a candidate for signal transduction

Among the predicted transcription factors for the DEG, NFIL3 was the most interesting candidate. Its binding motif was enriched in the clusters of both systemically up- and down-regulated DEG (clusters 2 and 3), in line with its function as activator/repressor of inflammatory and apoptotic genes. Furthermore, NFIL3 was found in infected quarters as 3.5-fold up-regulated DEG 24 h after infection. Interestingly, up-regulation of NFIL3 expression was found as soon as 6 h after infection. Five of the DEG at 6 h belong to the apparently systemically regulated genes at 24 h and only 3 of them were found local in infected udder quarters. Interestingly at 6 h the DEG also comprised a couple of predicted NFIL3 target genes, and TFBS for NFIL3 were significantly enriched. The putative target genes regulated by NFIL3 suggest a systemic up-regulation of genes involved in cell proliferation (e.g. KLF6, FGF2) as well as apoptosis and a down-regulation of genes related to cell adhesion and extracellular matrix.

Analyses of mRNA expression by quantitative real-time RT-PCR

In order to verify the results of microarray analyses, we selected 14 genes as representatives of all SOTA-clusters (Figure 4) and used RT-qPCR to quantify their expression in infected udder quarters and in the neighboring saline treated or untreated quarters. The results were referenced against the expression in udder quarters of healthy animals. As summarized in Table 2 (Additional File 3), the values for induction and repression of transcription obtained by microarray analyses were clearly confirmed by quantitative PCR. In general the calculated fold-changes were higher for qPCR. For example, the expression of SAA3 was found to be 1825-fold induced in infected udder quarter 24 h after inoculation with E. coli 1303, whereas microarray hybridization detected a 52-fold induction.

Table 2 Validation of microarray results by quantitative real-time RT-PCR (qPCR)

In order to investigate the influence of the inoculation solution (2 ml of physiological saline), the expression of the selected 14 genes was determined in all sampled udder quarters of infected animals and the values of untreated udder quarters and mock inoculated udder quarters were compared. As shown in Additional File 4, the treatment of the neighboring udder quarters with saline did not provoke a significant difference in expression of these 14 genes.

The results of qPCR allowed an estimation of animal specific variation of gene expression. The ΔCP values summarized in Additional File 5 clearly demonstrate that the transcript levels of some genes like HP or MT2A differed up to one order of magnitude in infected udder quarters 24 h hours after inoculation with E. coli 1303, whereas other genes like PIR or LBP were found to be uniformly expressed among the different animals.

Discussion

In this work, the transcriptome response of the bovine mammary gland was analyzed after experimental infection with E. coli 1303. Previous observations in cattle [10, 12, 23] favor the hypothesis of a systemic defense reaction, which protects unaffected udder quarters. In order to study this phenomenon, we investigated the transcriptome of bovine mammary tissue after infection with E. coli 1303 in the early stage of infection (6 h) before onset of fever and in the late stage (24 h) after its decline. As described for E. coli 1303 [17] infection of a single udder quarter induced acute clinical mastitis in all inoculated animals with comparable clinical signs such as strong fever, leukopenia, decreased milk yield and influx of somatic cells into the infected quarter. The somatic cell count (SCC) in mock inoculated as well as untreated udder quarters of the infected animals remained unchanged at all time points.

At 6 h after infection only 13 DEG were captured solely in infected quarters, indicating subtle effects of either the pathogen or the inoculation solution itself on the mammary gland at this time point. Furthermore we did not observe enhanced expression of proinflammatory cytokines like TNF-alpha. The missing induction of proinflammatory molecules in the early phase of infection and the unchanged SCC in mock inoculated udder quarters up to 24 h proves that the inoculation procedure neither provoked inflammation nor an increase of somatic cells in milk.

At 24 h after infection we observed marked changes in gene expression in neighboring quarters. These could be caused by a systemic response to the pathogen or a direct effect of the mock inoculation with 2 ml of PBS. As described by Paape et al. [24] the inoculation of even 100 ml saline solution provoked only a mild increase of SCC. Our procedure based on the application of 2 ml saline did not change the SCC of mock inoculated udder quarters as shown in Figure 2A. Furthermore the analysis of transcript levels of 14 selected genes by qPCR did not show significant differences between untreated and mock inoculated udder quarters (Additional File 4) and proved the existence of a systemic reaction. A detailed look on the setup of the microarray hybridizations (Additional File 1) reveals, that four profiles (n24h-1, 2, 3, 5) are based on pooling of samples from a mock inoculated and an untreated udder quarter, while the profile n24h-4 was obtained from a single untreated quarter. Assuming that the transcriptome changes observed in the neighboring quarters are predominantly caused by the saline treatment, the profile n24h-4 should cluster together with the controls c24h1-5. However, the highly similar expression profiles of the untreated sample n24h-4 and other neighboring samples (Figure 3C) show that the systemic effect of the pathogen by far dominates a possible effect of the saline treatment. The notion that 2 ml PBS has little influence is further substantiated by the assignment of profile n24h-5 to the cluster of the controls. Hence the saline treatment did not markedly alter the transcriptome of the pooled sample n24h-5 compared to the control group. Taken together the above observations strongly suggest that the experimental infection of an udder quarter caused an effect on the neighboring udder quarters, which we believe to be mediated systemically.

By comparison with expression profiles of healthy animals, large numbers of DEG were identified at 24 h in infected and a lesser number in neighboring noninfected quarters. Cluster analysis identified at least two categories of DEG. The first was found in infected quarters and comprised the local response of the infected udder quarter. DEG of the second category were found in inoculated and noninfected udder quarters of infected animals. This presumptive systemic response was discovered by comparison of noninfected neighboring quarters and samples of healthy animals. In contrast, most experimental designs used in mastitis research lack external controls and are thus not able to differentiate between local and systemic responses to the pathogen [12, 17, 25, 26]. A substantial systemic transcriptome response in noninfected neighboring quarters was proposed [27], but the experimental outline precluded further characterization. Recently, systemic effects of microbial induced mastitis have been found in the bovine udder by analysis of composition and viability of somatic cells in milk [23] and in the liver transcriptome of intramammary LPS-challenged cows [16]. Comparison of the latter with our findings by gene set enrichment analyses resulted in a higher coincidence with the systemic response than with the local response of infected udder quarters (data not shown).

Cluster analyses revealed very similar expression profiles of all infected udder quarters and a pronounced response to infection, probably due to the pathogenicity of E. coli 1303. In contrast, the expression profiles of noninfected udder quarters were more variable and one was rather similar to untreated animals (n24h-5 in Figure 3C). This observation clearly demonstrates that animal specific factors affected the systemic reaction, although the animals included in our study were selected according to strict criteria regarding the udder health.

The DEG were analyzed in the context of biological pathways, gene ontologies and predicted transcription factor binding sites to facilitate a better understanding of the local and the systemic reaction. Generally we found a clear division of biological functions between both types of reactions (Figure 5). DEG of the local reaction are mainly involved in immune response, inflammation, acute phase response and chemokine/cytokine signaling. TFBS analyses based on human orthologs predicted the corresponding transcription factors like NFκB and STAT1. As shown by Zadissa et al. [28], the use of promoter sequences from human genes for the inference of bovine transcriptional regulators is appropriate. The predicted factor ELF5 probably refers to DEG related to the restoration of damaged epithelial tissue, which is also supported by the enrichment of ontology terms like blood coagulation, cell proliferation and differentiation in cluster 1 (Figure 5). The majority of the DEG involved in immune response was assigned to the local reaction and only a few were also found in neighboring noninfected quarters. In this context, a detailed analysis of expression of chemokine genes led to a remarkable result. Enhanced expression of CC-type chemokine genes was exclusively assigned to the local reaction, whereas many CXC chemokine genes were regulated systemically. In contrast to the CC-type, CXC chemokines are known to attract neutrophils. The enhanced systemic expression of genes involved in response to oxidative stress accounts for the influx of active neutrophils in all quarters. This corresponds with the finding, that acute microbial mastitis caused a significant increase of the percentage of neutrophils in neighboring unaffected quarters [23]. The systemic reaction comprises biological processes like antigen processing and presentation, protein degradation, apoptosis, autophagy and genes of cytokines. Compared with the local reaction it includes only parts of an immune reaction and could be considered as well balanced conditioning process facilitating the progression of the immune response. This is for example illustrated by the induction of genes related to apoptosis, which was described to counteract adherence of bacteria in epithelia by exfoliation of the outermost layers [29]. The early onset of this defense reaction is indicated by enhanced expression of NFIL3 already 6 h after infection. In the late stage of infection, NFIL3 was predicted as regulatory element of many DEG involved in the local as well as the systemic reaction. Thus, the NFIL3 signal might promote a remodeling of the epithelium in both, infected and adjacent udder quarters. Studies in mice suggested that Nfil3 is involved in the functional regulation of the mammary epithelium [30]. Altered epithelial cell functions might cause the decrease in milk production observed in our study in all quarters (Figure 2D). Furthermore the ontologies 'lactation' and 'lipid metabolism' were prominent among systemically down-regulated DEG (cluster 3). This correlates with the finding that infection with E. coli caused reduction in milk yield in infected as well as in noninfected neighboring udder quarters (Figure 2D). Interestingly, at 24 h after infection we observed enhanced systemic expression of antiapoptotic genes (e.g. CFLAR, MCL1). This may be part of a counterbalance reaction to normalize the function of the mammary epithelium after the observed decline of fever.

Our results favor the hypothesis that very early after contact with E. coli signals are transmitted to the surrounding tissue and to the neighboring quarters. They induce a systemic reaction in all udder quarters that apparently prevents multiplication of the pathogen in neighboring quarters and influences the progression in the infected quarter. It explains, why time-shifted inoculations of different udder quarters with the same (low) pathogen titer caused influx of leukocytes only in the udder quarter first infected [12, 17]. Our hypothesis is also in line with the phenomenon of recurrent mastitis in different udder quarters caused by the same environmental strain of E. coli[10]. The pathogen may be transferred from the first infected udder quarter to neighboring quarters, but the ongoing systemic reaction inhibits multiplication of the pathogen in the time delayed infected quarters.

By integration of our findings into a consensus scheme of the transcriptional response to infections [31] we could provide a detailed picture of the processes occurring in the bovine mammary gland 24 h after inoculation with E. coli 1303 (Figure 6) and add a new functional dimension by differentiating locally and systemically regulated genes. The latter group points to altered milk composition and synthesis of protection factors, such as antimicrobials and reactive oxygen species in the late stage of infection. The systemic reaction could effectively prevent growth of E. coli and formation of clinical signs in neighboring udder quarters and could thereby moderate neutrophil influx. The significance of our results is evidenced by some recent findings with candidate gene based approaches. We found that TLR2 was 6-fold up-regulated 24 h after infection solely in infected quarters and accordingly TLR4 2-fold. Enhanced expression of TLR2 and TLR4 was measured by RT-PCR in the pathogen inoculated udder quarters and furthermore expression of TLR2 was detected by immunohistochemistry in mammary epithelial cells [17]. Milk of infected animals was described to contain elevated levels of the acute phase proteins haptoglobin (HP), serum amyloid A (SAA), lipopolysaccharide binding protein (LBP) [32] and alpha-1-antitrypsin (SERPINA1) [33]. We found HP 200-fold up-regulated locally in infected quarters, while SAA3, LBP and SERPINA1 showed systemically increased expression in infected (SAA3: 60-fold; LBP: 20-fold; SERPINA1: 10-fold) as well as in neighboring uninfected quarters (SAA3: 3-fold; LBP: 4-fold; SERPINA1: 3-fold). S100 calcium-binding protein A12 (S100A12) which was shown to have a direct antimicrobial effect on E. coli[34] was also identified in milk of cows with E. coli induced mastitis [33]. In our study S100A12 was up-regulated in the infected (110-fold) as well as in the neighboring quarters (4-fold). Previously high levels of interleukin 8 were measured in milk of infected cows [35]. We found 120-fold enhanced expression of IL8 in infected udder quarters 24 h after infection and no significant change in neighboring quarters. IL8 was shown to activate and attract neutrophils [36]. The observed enhanced expression of IL8 in infected quarters might be involved in the dramatic influx of somatic cells observed in these quarters only (Figure 2A). The increase of SSC was described to correlate with enhanced expression of lingual antimicrobial peptide (LAP) and its expression in epithelium was confirmed by in situ hybridization solely in the infected udder quarter [37]. At 24 h after infection, we found a 20-fold induction of LAP solely in infected quarters. These examples demonstrate the correlation of our findings with results of previously published candidate gene based investigations of bovine mastitis.

Figure 6
figure 6

Local and systemic immune response of the bovine udder 24 h after infection. A graphical representation of a robust response to infection as detected by a series of microarray experiments was adapted from [31]. Genes that were differentially expressed in the bovine udder upon E. coli infection are listed, further differentiating between systemic reactions in blue and local reactions in red. Bold letters indicate accordance with [31]. Genes were up-regulated unless marked with an arrow, indicating that expression is repressed upon infection.

Conclusion

In summary, this first comprehensive description of complex systemic effects in the mammary gland after local E. coli infection introduces a new concept to mastitis research and may yield new targets for prevention and therapy of mastitis. Moreover, the observed variation in the systemic response between individual animals could be the basis for selective breeding towards resistance against E. coli induced mastitis.

Methods

Animals

This study included 15 healthy German Holstein Frisian heifers in mid lactation (3 to 6 months post partum). The trials were conducted at the Clinic for Ruminants, LMU Munich (Oberschleißheim, Germany) with the approval of the ethics committee of the regional government of Upper Bavaria, Germany (No. 55.2-1-54-2531-108-05). Only animals without previous diagnosis of clinical or subclinical mastitis and a reported somatic cell count <50,000/ml were included in the study. Quarter milk samples were collected and tested weekly before the trial to ensure that they contained <50,000 somatic cells/ml and were free of mastitis pathogens. Two different infection models were used in which the animals were inoculated in one quarter with E. coli and killed after 6 h (n = 5) or 24 h (n = 5), respectively. Five heifers served as controls, received no treatment and were killed after 24 h.

Somatic cell count (SCC) and bacteriological examination

Three weeks before the trial the following parameters were documented twice daily for each animal: milk yield, somatic cell count (SCC, California Mastitis Test, CMT) and rectal body temperature. Quarter milk samples were collected weekly for determination of the SCC and bacteriological examination. The SCC was determined with a Fossomatic 5000® (FOSS Electric, Hillerod, Denmark) cell counter. To check for the absence of bacteria, 15 μl quarter milk samples were plated on Columbia sheep blood agar, Gassner agar and Edwards agar (Oxoid, Wesel, Germany). Samples were incubated for 48 h at 37°C.

Inoculum dose

Bovine mastitis isolate E. coli 1303 belongs to the major E. coli phylogenetic group A (E. coli collection of reference strains, ECOR-A). According to a multiplex PCR-based screening for virulence-associated genes of pathogenic E. coli, this strain does not represent an extraintestinal or intestinal pathogenic E. coli isolate. Only the genes coding for type 1 fimbriae, F17 fimbriae, antigen 43, the ferric citrate siderophore system and the EAST1 toxin could be detected. Bacteria were kept cryo-conserved (Mikrobank-System Cryobank™, Mast Diagnostika, Reinfeld, Germany) for subsequent infections. Bacteria were plated on Columbia sheep blood agar and incubated (37°C) for 24 h. A few colonies were transferred to a tube of brain-heart infusion broth (Oxoid, Wesel, Germany) and incubated for 6 h (37°C); then a 100-μl sample was transferred to a tube of trypticase soy broth (9.9 ml, Oxoid, Wesel, Germany). Serial dilutions were made after 18 h to prepare the desired inoculum dose of 500 cfu/2 ml 0.9% sterile, pyrogen-free saline. The inoculum dose was plated for control and ranged from 421-613 cfu/2 ml.

Experimental pathogen inoculation, sampling of milk, blood, and tissue

To achieve comparable hormonal conditions, cows were synchronized by two injections of (+)-cloprostenol (Dalmazin®, Selectavet, Weyarn, Germany) at 12-day intervals. Animals were challenged three days after the second injection during estrus.20 IU oxytocin (Veyx-Pharma, Schwarzenborn, Germany) were applied intravenously and the udder was entirely milked out. Teats were cleaned and disinfected with 70% ethanol and 500 cfu E. coli strain 1303 were administered intracisternally in one quarter through the teat canal. Two milliliters of 0.9% sterile pyrogen-free saline without bacteria were inoculated into another quarter as placebo. The control group was treated the same way but without inoculation of an udder quarter.

The udder secretions were sampled in 12 h intervals before regular milking. Bacteriological examination was performed from foremilk samples, the SCC was determined in 10 ml whole milk after each udder quarter had been milked separately with a quarter milker (WestfaliaSurge, Bönen, Germany).

Blood samples (10 ml) were taken from the jugular vein aseptically using EDTA-vacutainers (Becton Dickinson, Heidelberg, Germany) at 0, 3, 6, 12 and 24 h after start of the trial. Finally, cows were killed 6 h or 24 h after trial start with a penetrating captive bolt gun followed by exsanguination. Leukocyte counts in blood samples were carried out with a Sysmex pocH100i (Sysmex, Norderstedt, Germany). Tissue samples were collected aseptically from slaughtered cows within 10 min after killing. A piece of tissue (5 × 5 × 5 cm) was removed from a deeper location of the udder quarter, 7 cm dorsal of the milk cistern, of which a smaller tissue piece (1 × 1 × 0.5 cm) was transferred to a tube containing 5 ml RNAlater® (Applied Biosystems Incorporation, Foster City, USA) and served as original material for all further analyses.

Microarray hybridization

Total RNA was isolated by disrupting small pieces of approx. 500 mg of bovine mammary gland in 5 ml of TRIzol reagent (Invitrogen, Carlsbad, USA) with a tissue homogenizer (Heidolph, Schwabach, Germany). Homogenate was processed following the TRIzol manufacturer's protocol to yield highly pure total RNA. Purity (A260/A280 >1.9, A260/A230 >2.2) and integrity of RNA (28S:18S rRNA ratio ~1.5:1) was assessed by spectrophotometry (Nanodrop ND-100, NanoDrop Technologies LLC, Wilmington, DE, USA) and agarose gel electrophoresis.

Labeled cRNA probes for array hybridization of the 24 h samples (infected udder quarters, neighboring quarters and controls) were prepared with an Ambion Kit (MessageAmp™ II-Biotin Enhanced Single Round aRNA Amplification Kit). The probes for analyses of the 6 h time point (infected udder quarters and controls) were prepared with an Affymetrix kit (GeneChip® Expression 3' Amplification One-Cycle Target Labeling Kit). RNA samples from noninfected udder quarters of the same animal were pooled when they had comparable somatic cells counts (detailed description of samples: see Additional File 1). The controls were performed as replicates to ensure comparability of the sets.

10 μg of the resulting fragmented cRNA were hybridized overnight to GeneChip® Bovine Genome Arrays (Affymetrix, Santa Clara, CA, USA) processed in an Affymetrix Fluidics Station 450 and scanned with the Affymetrix 3000 7G scanner. Hybridizations were performed in three sets of 10 arrays each: (set 1) 5 infected quarters 24 h after infection and 5 controls; (set 2) 5 samples from neighboring quarters 24 h after infection (see Additional File 1) and 5 controls; (set 3) 5 challenged quarters 6 h after infection and 5 samples from untreated controls. The statistic software R and packages from the bioconductor [38] microarray suite (affy, affyPLM) were used for processing the CEL files and to identify possible outliers due to technical or handling artifacts. The microarray data has been deposited in the Gene Expression Omnibus database (accession GSE15025).

Statistical analysis and bioinformatics

Microarray raw data were summarized and normalized with the RMA method [39] and filtered for a detection criterion (MAS5 'present'-calls determined by Bioconductor's affy-package in R). DEG were identified by the packages LIMMA [40] and SAM [41]. LIMMA was performed for the analysis of the sample sets 24 h after infection, excluding profiles of animal 5 (samples i24h-5 and n24h-5) as a biological outlier and control animal 3 (sample c124h-3 and c224h-3) as a technical outlier. Due to its higher sensitivity in case of noisy data SAM was used for analysis of the 6 h time point. Cut-off criteria were a higher than 2-fold change and a false discovery rate of <1% for LIMMA and <10% for SAM analysis. The expression values of set 1 and set 2 for 2335 DEG were subjected to cluster analysis using the SOTA module of MeV4.2 [42] with the Pearson correlation as clustering metric and default values except for Winning Cell Migration Weight (0.006). The number of clusters was set to six, based on figure-of-merit calculations. The members of the resulting clusters were subjected to enrichment analyses with different online tools e.g. pantherdb [43] and CoPub [44] for GO terms and pathways and oPOSSUM [45] for TFBS. These analyses were based on human orthologs of the DEG as described by Hintermair [46] and the human genome or the orthologs present on the microarray as the background. Enriched functional annotations and pathways are represented as a heatmap of the Bonferroni-adjusted p-values with a cut-off of 0.05. The human ortholog gene symbols of differentially expressed bovine genes are listed in Additional File 2.

Real-Time RT-PCR

Quantitative real-time PCR reactions were performed with the same RNA samples as used for microarray hybridization using the LightCycler DNA Master SYBR Green I protocol [47]. The following primers were designed to amplify specific fragments referring to selected differentially expressed genes of each cluster (length of the PCR products in square brackets): cluster 1: HP (haptoglobin; for 5'-ACCTGGTATGCGGCCGGGA; rev 5'-TCCGAACCCAGTCCAGAATGGAGG [100 bp]), SAA3 (serum amyloid A3; for 5'-ATGACGCTGCCCGAAGGGGA; rev 5'-TGTCAGGCAGGCCAGCAGGT [210 bp]), PTX3 (pentraxin 3; for 5'-AAAGGGAGACTGGAGAAGGC; rev 5'-TGGCCAAAATGAAATTAACCA [173 bp]), cluster 2: S100A8 (S100 calcium binding protein A8; for 5'-AAAAAGGGAATTACCACGCC; rev 5'-ATCACCAGCACGAGGAACTC [163 bp]), MT2A (metallothionein 2A; for 5'-TCTCCGGACCCCAGCCTCCA; rev 5'-GCAGGAGGCGCACTTGCAATC [119 bp]), LBP (lipopolysaccharide binding protein; for 5'-CCTTGCCCTCAAACTCTCAG; rev 5'-CTTCCCCCTCCCTCTGTTAC [114 bp]), cluster 3: AQP3 (aquaporin 3; for 5'-TCGGTGGAGTTGGGTGGGGG; rev 5'-AGCCCCCTGAATAGAAGAAAGGGC [158 bp]), ALOX15 (arachidonate 15-lipoxygenase; for 5'-GGGTGGGATTCACCACGTGTCC; rev 5'-ACTGAGGCCAGATACCTCCAAC [170 bp]), cluster 4: ACAS2 (acyl-CoA synthetase short-chain family member 2; for 5'-GGCATGCACTTGCCCCGAGA; rev 5'-ACAGGCACTGCCATCCGGT [122 bp]), ADHFE1 (alcohol dehydrogenase iron containing 1; for 5'-AGGGCAGCCTGGACAAACGC; rev 5'-GATCCATGGAGCAAGGGCAGTC [102 bp]), cluster 5: PIR (pirin; for 5'-TTGGACCTGATGATGCACAGCAAA; rev 5'-ACGTGGACACTGTCACCTTCTCCC [76 bp]), THRSP (thyroid hormone responsive; for 5'-AAGTCAAAAGCATCTGGCATGT; rev 5'-ACTCAGAGTTGAGGACTCGGCTT [134 bp]), cluster 6: HPGD (hydroxyprostaglandin dehydrogenase; for 5'-TCTGTCTTCCACTGTAATGCTCAAAGC; rev 5'-AGTGAAGGCCACCAGTACAAAGACT [80 bp]), LPL (lipoprotein lipase; for 5'-TGACTTGTTGTTGGCATCCCCC; rev 5'-AAGTCAGAGTTCCCAGGGCCA [129 bp]) and SF3A1 (splicing factor 3a, subunit 1; for 5'-ACAAGGGTCCAGTGTCCATC; rev 5'-AGACCAGCACCTGTCCATTC [84 bp])) as housekeeping gene [48]. All amplified PCR fragments were sequenced (3100-Avant Genetic Analyzer; Applied Biosystems) to verify the resulting PCR product. In addition, the specific melting point (MP) of the amplified product carried out within the LightCycler standard PCR protocol confirmed product identity [49]. In each PCR reaction cDNA according to 50 ng total RNA was amplified in a total reaction volume of 20 μl (5 μM primer forward and reverse each, 1x LightCycler DNA Master SYBR Green I; Roche) using the LightCycler480 II instrument (Roche). The annealing temperature was 60°C for all PCRs.

Data analysis of real-time RT-PCR

The cycle number required to achieve a definite SYBR Green fluorescence signal (= crossing point; CP) was calculated by the second derivative maximum method (LightCycler software, version 1.5.0). The CP is correlated inversely with the logarithm of the initial template concentration. The CP of the housekeeping gene SF3A1 showed no significant statistical difference in all analyzed samples. Therefore, it was used to normalize the CP for the target genes (ΔCP). Differences between values obtained for udder quarters of the control group and infected or neighboring quarters were stated by the ΔΔCP as well as the 'mean fold changes qPCR' [50]. The normal distribution was tested by the Kolmogorow-Smirnov method, followed by a Student's t-test to confirm significant differences between groups.

Abbreviations

c:

udder quarter of healthy cows

cfu:

colony forming units

CP:

crossing point

cRNA:

complementary ribonucleic acid

DEG:

probe sets of differentially expressed genes

E. coli :

Escherichia coli

fdr:

false discovery rate

GO:

gene ontology

i:

udder quarter inoculated with E. coli

LIMMA:

linear models for microarray analysis

LPS:

lipopolysaccharide

MAS:

microarray suite

n:

udder quarter neighboring to E. coli- infected udder quarter

qPCR:

quantitative polymerase chain reaction

RMA:

robust multi-array analysis

SAM:

significance analysis of microarray data

SCC:

somatic cell count

SOTA:

self organizing tree algorithm

TFBS:

transcription factor binding site.

References

  1. Bar D, Tauer LW, Bennett G, Gonzalez RN, Hertl JA, Schukken YH, Schulte HF, Welcome FL, Grohn YT: The cost of generic clinical mastitis in dairy cows as estimated by using dynamic programming. J Dairy Sci. 2008, 91 (6): 2205-2214. 10.3168/jds.2007-0573.

    Article  CAS  PubMed  Google Scholar 

  2. Seegers H, Fourichon C, Beaudeau F: Production effects related to mastitis and mastitis economics in dairy cattle herds. Vet Res. 2003, 34 (5): 475-491. 10.1051/vetres:2003027.

    Article  PubMed  Google Scholar 

  3. Compton CW, Cursons RT, Barnett CM, McDougall S: Expression of innate resistance factors in mammary secretion from periparturient dairy heifers and their association with subsequent infection status. Vet Immunol Immunopathol. 2009, 127 (3-4): 357-364. 10.1016/j.vetimm.2008.10.331.

    Article  CAS  PubMed  Google Scholar 

  4. Takeshima S, Matsumoto Y, Chen J, Yoshida T, Mukoyama H, Aida Y: Evidence for cattle major histocompatibility complex (BoLA) class II DQA1 gene heterozygote advantage against clinical mastitis caused by Streptococci and Escherichia species. Tissue Antigens. 2008, 72 (6): 525-531. 10.1111/j.1399-0039.2008.01140.x.

    Article  CAS  PubMed  Google Scholar 

  5. Pyorala S: Mastitis in post-partum dairy cows. Reprod Domest Anim. 2008, 43 (Suppl 2): 252-259.

    PubMed  Google Scholar 

  6. 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 (5): 521-564. 10.1051/vetres:2003023.

    Article  PubMed  Google Scholar 

  7. Vandeputte-Van Messom G, Burvenich C, Roets E, Massart-Leen AM, Heyneman R, Kremer WD, Brand A: Classification of newly calved cows into moderate and severe responders to experimentally induced Escherichia coli mastitis. J Dairy Res. 1993, 60 (1): 19-29. 10.1017/S002202990002731X.

    Article  CAS  PubMed  Google Scholar 

  8. Hill AW: Escherichia coli mastitis. 1994, CAB international, Wallingford

    Google Scholar 

  9. Adkinson RW, Ingawa KH, Blouin DC, Nickerson SC: Distribution of clinical mastitis among quarters of the bovine udder. J Dairy Sci. 1993, 76 (11): 3453-3459.

    Article  CAS  PubMed  Google Scholar 

  10. Dopfer D, Barkema HW, Lam TJ, Schukken YH, Gaastra W: Recurrent clinical mastitis caused by Escherichia coli in dairy cows. J Dairy Sci. 1999, 82 (1): 80-85.

    Article  CAS  PubMed  Google Scholar 

  11. Hirvonen J, Eklund K, Teppo AM, Huszenicza G, Kulcsar M, Saloniemi H, Pyorala S: Acute phase response in dairy cows with experimentally induced Escherichia coli mastitis. Acta Vet Scand. 1999, 40 (1): 35-46.

    CAS  PubMed  Google Scholar 

  12. 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-kappaB in mammary epithelial cells and to quickly induce TNFalpha and interleukin-8 (CXCL8) expression in the udder. Mol Immunol. 2008, 45 (5): 1385-1397. 10.1016/j.molimm.2007.09.004.

    Article  CAS  PubMed  Google Scholar 

  13. 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 (2): 131-144. 10.1051/vetres:2001100.

    Article  CAS  PubMed  Google Scholar 

  14. Blum JW, Dosogne H, Hoeben D, Vangroenweghe F, Hammon HM, Bruckmaier RM, Burvenich C: Tumor necrosis factor-alpha and nitrite/nitrate responses during acute mastitis induced by Escherichia coli infection and endotoxin in dairy cows. Domest Anim Endocrinol. 2000, 19 (4): 223-235. 10.1016/S0739-7240(00)00079-5.

    Article  CAS  PubMed  Google Scholar 

  15. Hiss S, Mielenz M, Bruckmaier RM, Sauerwein H: Haptoglobin concentrations in blood and milk after endotoxin challenge and quantification of mammary Hp mRNA expression. J Dairy Sci. 2004, 87 (11): 3778-3784.

    Article  CAS  PubMed  Google Scholar 

  16. Jiang L, Sorensen P, Rontved C, Vels L, Ingvartsen KL: Gene expression profiling of liver from dairy cows treated intra-mammary with lipopolysaccharide. BMC Genomics. 2008, 9: 443-10.1186/1471-2164-9-443.

    Article  PubMed Central  PubMed  Google Scholar 

  17. Petzl W, Zerbe H, Gunther J, Yang W, Seyfert HM, Nurnberg 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 (2): 18-10.1051/vetres:2007057.

    Article  PubMed  Google Scholar 

  18. Herrero J, Valencia A, Dopazo J: A hierarchical unsupervised growing neural network for clustering gene expression patterns. Bioinformatics. 2001, 17 (2): 126-136. 10.1093/bioinformatics/17.2.126.

    Article  CAS  PubMed  Google Scholar 

  19. Zhou J, Chehab R, Tkalcevic J, Naylor MJ, Harris J, Wilson TJ, Tsao S, Tellis I, Zavarsek S, Xu D, Lapinskas EJ, Visvader J, Lindeman GJ, Thomas R, Ormandy CJ, Hertzog PJ, Kola I, Pritchard MA: Elf5 is essential for early embryogenesis and mammary gland development during pregnancy and lactation. Embo J. 2005, 24 (3): 635-644. 10.1038/sj.emboj.7600538.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  20. Harris J, Stanford PM, Sutherland K, Oakes SR, Naylor MJ, Robertson FG, Blazek KD, Kazlauskas M, Hilton HN, Wittlin S, Alexander WS, Lindeman GJ, Visvader JE, Ormandy CJ: Socs2 and elf5 mediate prolactin-induced mammary gland development. Mol Endocrinol. 2006, 20 (5): 1177-1187. 10.1210/me.2005-0473.

    Article  CAS  PubMed  Google Scholar 

  21. Yang D, Chen Q, Hoover DM, Staley P, Tucker KD, Lubkowski J, Oppenheim JJ: Many chemokines including CCL20/MIP-3alpha display antimicrobial activity. J Leukoc Biol. 2003, 74 (3): 448-455. 10.1189/jlb.0103024.

    Article  CAS  PubMed  Google Scholar 

  22. Pascual M, Gomez-Lechon MJ, Castell JV, Jover R: ATF5 is a highly abundant liver-enriched transcription factor that cooperates with constitutive androstane receptor in the transactivation of CYP2B6: implications in hepatic stress responses. Drug Metab Dispos. 2008, 36 (6): 1063-1072. 10.1124/dmd.107.019380.

    Article  CAS  PubMed  Google Scholar 

  23. Merle R, Schroder 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 (2): 174-179. 10.1017/S002202990600238X.

    Article  CAS  PubMed  Google Scholar 

  24. Paape MJ, Pearson RE, Wergin WP, Guidry AJ: Enhancement of chemotactic response of polymorphonuclear leukocytes into the mammary gland and isolation from milk. J Dairy Sci. 1977, 60 (1): 53-62.

    Article  CAS  PubMed  Google Scholar 

  25. 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 (1): 117-129. 10.3168/jds.2008-1382.

    Article  CAS  PubMed  Google Scholar 

  26. Zheng J, Watson AD, Kerr DE: Genome-wide expression analysis of lipopolysaccharide-induced mastitis in a mouse model. Infect Immun. 2006, 74 (3): 1907-1915. 10.1128/IAI.74.3.1907-1915.2006.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  27. Jaffrezic F, de Koning DJ, Boettcher PJ, Bonnet A, Buitenhuis B, Closset R, Dejean S, Delmas C, Detilleux JC, Dovc P, Duval M, Foulley JL, Hedegaard J, Hornshoj H, Hulsegge I, Janss L, Jensen K, Jiang L, Lavric M, Le Cao KA, Lund MS, Malinverni R, Marot G, Nie H, Petzl W, Pool MH, Robert-Granie C, San Cristobal M, van Schothorst EM, Schuberth HJ: Analysis of the real EADGENE data set: comparison of methods and guidelines for data normalisation and selection of differentially expressed genes (open access publication). Genet Sel Evol. 2007, 39 (6): 633-650. 10.1186/1297-9686-39-6-633.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  28. Zadissa A, McEwan JC, Brown CM: Inference of transcriptional regulation using gene expression data from the bovine and human genomes. BMC Genomics. 2007, 8: 265-10.1186/1471-2164-8-265.

    Article  PubMed Central  PubMed  Google Scholar 

  29. Mysorekar IU, Mulvey MA, Hultgren SJ, Gordon JI: Molecular regulation of urothelial renewal and host defenses during infection with uropathogenic Escherichia coli. J Biol Chem. 2002, 277 (9): 7412-7419. 10.1074/jbc.M110560200.

    Article  CAS  PubMed  Google Scholar 

  30. Rudolph MC, McManaman JL, Phang T, Russell T, Kominsky DJ, Serkova NJ, Stein T, Anderson SM, Neville MC: Metabolic regulation in the lactating mammary gland: a lipid synthesizing machine. Physiol Genomics. 2007, 28 (3): 323-336.

    Article  CAS  PubMed  Google Scholar 

  31. Jenner RG, Young RA: Insights into host responses against pathogens from transcriptional profiling. Nat Rev Microbiol. 2005, 3 (4): 281-294. 10.1038/nrmicro1126.

    Article  CAS  PubMed  Google Scholar 

  32. Suojala L, Orro T, Jarvinen H, Saatsi J, Pyorala S: Acute phase response in two consecutive experimentally induced E. coli intramammary infections in dairy cows. Acta Vet Scand. 2008, 50: 18-10.1186/1751-0147-50-18.

    Article  PubMed Central  PubMed  Google Scholar 

  33. Boehmer JL, Bannerman DD, Shefcheck K, Ward JL: Proteomic analysis of differentially expressed proteins in bovine milk during experimentally induced Escherichia coli mastitis. J Dairy Sci. 2008, 91 (11): 4206-4218. 10.3168/jds.2008-1297.

    Article  CAS  PubMed  Google Scholar 

  34. Lutzow YC, 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-10.1186/1746-6148-4-18.

    Article  PubMed Central  PubMed  Google Scholar 

  35. Bannerman DD, Kauf AC, Paape MJ, Springer HR, Goff JP: Comparison of Holstein and Jersey innate immune responses to Escherichia coli intramammary infection. J Dairy Sci. 2008, 91 (6): 2225-2235. 10.3168/jds.2008-1013.

    Article  CAS  PubMed  Google Scholar 

  36. Caswell JL, Middleton DM, Sorden SD, Gordon JR: Expression of the neutrophil chemoattractant interleukin-8 in the lesions of bovine pneumonic pasteurellosis. Vet Pathol. 1998, 35 (2): 124-131. 10.1177/030098589803500206.

    Article  CAS  PubMed  Google Scholar 

  37. Swanson K, Gorodetsky S, Good L, Davis S, Musgrave D, Stelwagen K, Farr V, Molenaar A: Expression of a beta-defensin mRNA, lingual antimicrobial peptide, in bovine mammary epithelial tissue is induced by mastitis. Infect Immun. 2004, 72 (12): 7311-7314. 10.1128/IAI.72.12.7311-7314.2004.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  38. Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, Ellis B, Gautier L, Ge Y, Gentry J, Hornik K, Hothorn T, Huber W, Iacus S, Irizarry R, Leisch F, Li C, Maechler M, Rossini AJ, Sawitzki G, Smith C, Smyth G, Tierney L, Yang JY, Zhang J: Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004, 5 (10): R80-10.1186/gb-2004-5-10-r80.

    Article  PubMed Central  PubMed  Google Scholar 

  39. Irizarry RA, Bolstad BM, Collin F, Cope LM, Hobbs B, Speed TP: Summaries of Affymetrix GeneChip probe level data. Nucleic Acids Res. 2003, 31 (4): e15-10.1093/nar/gng015.

    Article  PubMed Central  PubMed  Google Scholar 

  40. Smyth GK: Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004, 3 (12): Article3

    Google Scholar 

  41. Tusher VG, Tibshirani R, Chu G: Significance analysis of microarrays applied to the ionizing radiation response. Proc Natl Acad Sci USA. 2001, 98 (9): 5116-5121. 10.1073/pnas.091062498.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  42. Multiexperiment Viewer MeV V4.2.02. [http://www.tm4.org/mev.html]

  43. Pantherdb. [http://www.pantherdb.org/]

  44. Frijters R, Heupers B, van Beek P, Bouwhuis M, van Schaik R, de Vlieg J, Polman J, Alkema W: CoPub: a literature-based keyword enrichment tool for microarray data analysis. Nucleic Acids Res. 2008, W406-410. 10.1093/nar/gkn215. 36 Web Server

  45. Ho Sui SJ, Fulton DL, Arenillas DJ, Kwon AT, Wasserman WW: oPOSSUM: integrated tools for analysis of regulatory motif over-representation. Nucleic Acids Res. 2007, W245-252. 10.1093/nar/gkm427. 35 Web Server

  46. Hintermair V: Comparative Analysis of Bos taurus and Homo sapiens DNA Microarrays. Master's Thesis, Technical University, Munich;. 2007

    Google Scholar 

  47. Klein C, Bauersachs S, Ulbrich SE, Einspanier R, Meyer HH, Schmidt SE, Reichenbach HD, Vermehren M, Sinowatz F, Blum H, Wolf E: Monozygotic twin model reveals novel embryo-induced transcriptome changes of bovine endometrium in the preattachment period. Biol Reprod. 2006, 74 (2): 253-264. 10.1095/biolreprod.105.046748.

    Article  CAS  PubMed  Google Scholar 

  48. Perez R, Tupac-Yupanqui I, Dunner S: Evaluation of suitable reference genes for gene expression studies in bovine muscular tissue. BMC Mol Biol. 2008, 9: 79-10.1186/1471-2199-9-79.

    Article  PubMed Central  PubMed  Google Scholar 

  49. Bauersachs S, Rehfeld S, Ulbrich SE, Mallok S, Prelle K, Wenigerkind H, Einspanier R, Blum H, Wolf E: Monitoring gene expression changes in bovine oviduct epithelial cells during the oestrous cycle. J Mol Endocrinol. 2004, 32 (2): 449-466. 10.1677/jme.0.0320449.

    Article  CAS  PubMed  Google Scholar 

  50. Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001, 25 (4): 402-408. 10.1006/meth.2001.1262.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

The authors thank Sebastian Waszak for technical help in R programming, Angela Hafner-Marx for help in tissue sampling and Michael Salomon, Sirion Biotech GmbH for help in performing qPCR. This work was funded by the Deutsche Forschungsgemeinschaft (FOR 585).

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Helmut Blum.

Additional information

Authors' contributions

SM carried out the sample preparation, performed microarray analysis and qPCR, participated in bioinformatic analysis and contributed to writing of the manuscript. WP conducted the animal preparations and experimental infections, performed the evaluation of the clinical data and participated in writing the manuscript. SK performed bioinformatic and statistical analysis and contributed to drafting the manuscript. DM assisted in animal preparation and animal experiments. AK performed hybridization of microarrays. EW participated in the design of the study and performed writing and critical reading of the paper. HZ contributed to project conception and coordinated the animal experiments. HB conceived of the study, participated in interpretation of the results and drafted the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material

Additional file 1: Setup of microarray hybridizations and qPCR.(XLS 20 KB)

12864_2009_2732_MOESM2_ESM.XLS

Additional file 2: List of differentially expressed genes. List of differentially expressed genes Description Human ortholog gene symbols of differentially expressed bovine genes were used in the Additional File 2 and are based on unpublished data from Hintermair [46]. (XLS 440 KB)

12864_2009_2732_MOESM3_ESM.XLS

Additional file 3: Detailed data of validation of microarray results by quantitative real-time RT-PCR (qPCR). CP, crossing point; ΔCP, CP-CPSF3A1; ΔΔCP, ΔCPcontrol - ΔCPtreatment; Øn24h, average of untreated (n24h_u) and saline treated neighboring quarters (n24h_s) of 24 h infected quarters. (XLS 31 KB)

12864_2009_2732_MOESM4_ESM.PDF

Additional file 4: Expression profiles of 14 selected genes measured by real time RT-PCR. The mRNA expression of HP, SAA3, PTX3, S100A8, MT2A, LBP, AQP3, ALOX15, ACAS2, ADHFE1, PIR, THRSP, HPGD and LPL were detected using real time RT-PCR (qPCR). The height of the bars indicates the average fold changes for the i24h and n24h group - the latter separated in untreated (n24h_u) and saline treated (n24h_s) -relative to the mean expression of the control. The error bars indicate the standard error of mean. (PDF 18 KB)

12864_2009_2732_MOESM5_ESM.XLS

Additional file 5: Expression values of 14 selected genes in individual udder quarters measured by real time RT-PCR. The mRNA expression of HP, SAA3, PTX3, S100A8, MT2A, LBP, AQP3, ALOX15, ACAS2, ADHFE1, PIR, THRSP, HPGD and LPL were detected using real time RT-PCR (qPCR). The samples i24h-1 to 5 were derived from infected udder quarters of five animals 24 h after infection with E. coli 1303. The five samples subscripted with n24h1-5 were isolated from udder quarters neighboring to the infected quarters and were untreated (u) or mock inoculated with 2 ml PBS. The samples labeled with c24h were taken from two udder quarters of 5 untreated control animals. The listed CP values represent the expression of each mRNA in relation to SF3A1. (XLS 34 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Mitterhuemer, S., Petzl, W., Krebs, S. et al. Escherichia coli infection induces distinct local and systemic transcriptome responses in the mammary gland. BMC Genomics 11, 138 (2010). https://doi.org/10.1186/1471-2164-11-138

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2164-11-138

Keywords