Escherichia coli infection induces distinct local and systemic transcriptome responses in the mammary gland
© Mitterhuemer et al; licensee BioMed Central Ltd. 2010
Received: 30 June 2009
Accepted: 25 February 2010
Published: 25 February 2010
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.
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.
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.
Udder infections cause considerable economic losses to the dairy industry  and often lead to culling of affected animals . 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 . 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 . E. coli infection predominantly affects a single udder quarter . 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 . This phenomenon was attributed to initial infection of multiple quarters, each at a different level and with variations in lesion development .
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 . 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 or its lipopolysaccharide  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  or enhanced expression of several genes like haptoglobin in neighboring quarters . 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 . 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 . 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.
Establishment of intramammary infection
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.
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.
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
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.
Transcription factor binding site enrichment in clusters
Promoter range [kb]a
1 (local, 936 DEG)
2 (systemic, 703 DEG)
3 (systemic, 367 DEG)
4 (systemic, 41 DEG)
5 (systemic, 36 DEG)
6 (local, 252 DEG)
Characterization of the local response 24 h after infection
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  and milk production .
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 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 , 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
Validation of microarray results by quantitative real-time RT-PCR (qPCR)
Mean fold change qPCR
Mean fold change Affymetrix
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.
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  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.  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 , 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  and in the liver transcriptome of intramammary LPS-challenged cows . 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. , 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 . 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 . 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 . 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. 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.
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.
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.
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.
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  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  and filtered for a detection criterion (MAS5 'present'-calls determined by Bioconductor's affy-package in R). DEG were identified by the packages LIMMA  and SAM . 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  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  and CoPub  for GO terms and pathways and oPOSSUM  for TFBS. These analyses were based on human orthologs of the DEG as described by Hintermair  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.
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 . 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 . 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 . 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' . The normal distribution was tested by the Kolmogorow-Smirnov method, followed by a Student's t-test to confirm significant differences between groups.
udder quarter of healthy cows
colony forming units
complementary ribonucleic acid
probe sets of differentially expressed genes
- E. coli :
false discovery rate
udder quarter inoculated with E. coli
linear models for microarray analysis
udder quarter neighboring to E. coli- infected udder quarter
quantitative polymerase chain reaction
robust multi-array analysis
significance analysis of microarray data
somatic cell count
self organizing tree algorithm
transcription factor binding site.
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).
- 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.PubMedView Article
- 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.PubMedView Article
- 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.PubMedView Article
- 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.PubMedView Article
- Pyorala S: Mastitis in post-partum dairy cows. Reprod Domest Anim. 2008, 43 (Suppl 2): 252-259.PubMed
- 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.PubMedView Article
- 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.PubMedView Article
- Hill AW: Escherichia coli mastitis. 1994, CAB international, Wallingford
- 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.PubMedView Article
- 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.PubMedView Article
- 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.PubMed
- 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.PubMedView Article
- 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.PubMedView Article
- 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.PubMedView Article
- 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.PubMedView Article
- 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.PubMed CentralPubMedView Article
- 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.PubMedView Article
- 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.PubMedView Article
- 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.PubMed CentralPubMedView Article
- 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.PubMedView Article
- 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.PubMedView Article
- 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.PubMedView Article
- 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.PubMedView Article
- 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.PubMedView Article
- 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.PubMedView Article
- 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.PubMed CentralPubMedView Article
- 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.PubMed CentralPubMedView Article
- 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.PubMed CentralPubMedView Article
- 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.PubMedView Article
- 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.PubMedView Article
- Jenner RG, Young RA: Insights into host responses against pathogens from transcriptional profiling. Nat Rev Microbiol. 2005, 3 (4): 281-294. 10.1038/nrmicro1126.PubMedView Article
- 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.PubMed CentralPubMedView Article
- 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.PubMedView Article
- 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.PubMed CentralPubMedView Article
- 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.PubMedView Article
- 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.PubMedView Article
- 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.PubMed CentralPubMedView Article
- 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.PubMed CentralPubMedView Article
- 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.PubMed CentralPubMedView Article
- Smyth GK: Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004, 3 (12): Article3
- 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.PubMed CentralPubMedView Article
- Multiexperiment Viewer MeV V4.2.02. [http://www.tm4.org/mev.html]
- Pantherdb. [http://www.pantherdb.org/]
- 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
- 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
- Hintermair V: Comparative Analysis of Bos taurus and Homo sapiens DNA Microarrays. Master's Thesis, Technical University, Munich;. 2007
- 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.PubMedView Article
- 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.PubMed CentralPubMedView Article
- 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.PubMedView Article
- 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.PubMedView Article
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.