Email updates

Keep up to date with the latest news and content from BMC Genomics and BioMed Central.

This article is part of the supplement: Proceedings of the Avian Genomics Conference and Gene Ontology Annotation Workshop

Open Access Proceedings

The transcriptional response of Pasteurella multocida to three classes of antibiotics

Bindu Nanduri12*, Leslie A Shack12, Shane C Burgess1234 and Mark L Lawrence12

Author affiliations

1 College of Veterinary Medicine, Mississippi State University, Mississippi State, MS 39762, USA

2 Institute for Digital Biology, Mississippi State University, Mississippi State, MS 39762, USA

3 Life Sciences and Biotechnology Institute, Mississippi State University, Mississippi State, MS 39762, USA

4 Mississippi Agriculture and Forestry Experiment Station, Mississippi State University, Mississippi State, MS 39762, USA

For all author emails, please log on.

Citation and License

BMC Genomics 2009, 10(Suppl 2):S4  doi:10.1186/1471-2164-10-S2-S4


The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2164/10/S2/S4


Published:14 July 2009

© 2009 Nanduri et al; licensee 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.

Abstract

Background

Pasteurella multocida is a gram-negative bacterial pathogen that has a broad host range. One of the diseases it causes is fowl cholera in poultry. The availability of the genome sequence of avian P. multocida isolate Pm70 enables the application of functional genomics for observing global gene expression in response to a given stimulus. We studied the effects of three classes of antibiotics on the P. multocida transcriptome using custom oligonucleotide microarrays from NimbleGen Systems. Hybridizations were conducted with RNA isolated from three independent cultures of Pm70 grown in the presence or absence of sub-minimum inhibitory concentration (sub-MIC) of antibiotics. Differentially expressed (DE) genes were identified by ANOVA and Dunnett's test. Biological modeling of the differentially expressed genes (DE) was conducted based on Clusters of Orthologous (COG) groups and network analysis in Pathway Studio.

Results

The three antibiotics used in this study, amoxicillin, chlortetracycline, and enrofloxacin, collectively influenced transcription of 25% of the P. multocida Pm70 genome. Some DE genes identified were common to more than one antibiotic. The overall transcription signatures of the three antibiotics differed at the COG level of the analysis. Network analysis identified differences in the SOS response of P. multocida in response to the antibiotics.

Conclusion

This is the first report of the transcriptional response of an avian strain of P. multocida to sub-lethal concentrations of three different classes of antibiotics. We identified common adaptive responses of P. multocida to antibiotic stress. The observed changes in gene expression of known and putative P. multocida virulence factors establish the molecular basis for the therapeutic efficacy of sub-MICs. Our network analysis demonstrates the feasibility and limitations of applying systems modeling to high throughput datasets in 'non-model' bacteria.

Background

Pasteurella multocida is a gram–negative bacterial pathogen that has a unique history of serving as a model species for new discoveries. In his seminal experiments in the late eighteenth century, Louis Pasteur demonstrated both attenuation as well as protective immune response utilizing P. multocida in birds [1,2]. P. multocida is also an excellent model species for studying the effects of antibiotics on gram-negative bacteria because its short lipopolysaccharide O side chains make it more permeable, allowing investigations on the effects of antibiotics on bacterial metabolism [3].

Of the different P. multocida subspecies, multocida has a broad host range and causes diseases in poultry, cattle, pigs, and rabbits [4]. Zoonosis with P. multocida is caused by cat and dog bites and scratches and respiratory infection [5]. Fowl cholera (avian pasteurellosis) caused by P. multocida can be either chronic or acute and is a septicemic disease that causes significant economic loss to the poultry industry. Due to the importance of this pathogen to the avian community, the genome of an isolate recovered from chicken fowl cholera (strain Pm70) was sequenced in 2001 [6]. The genome sequence has enabled subsequent functional genomics research with this species [7-9]. It has also enabled investigations on the proteomic and transcriptomic response of P. multocida Pm70 to sub-MICs of antibiotics [10-12].

Much attention has been focused on characterization of specific antibiotic resistance mechanisms; however, recent studies are shedding light on non-target effects of antibiotics. In the current study, we used P. multocida as a model species to investigate the effects of three antibiotics with disparate mechanisms of action (amoxicillin, chlortetracycline, and enrofloxacin) on bacterial metabolism. Because biological systems utilize highly complex, interrelated metabolic and signaling networks, a thorough understanding of bacterial physiology requires studying the global response to a given stimulus in the context of interacting gene networks. Computational systems biology facilitates this level of analysis. Systems level analysis can identify key regulatory elements of a molecular interaction network; dynamic changes in these elements govern a response [13].

Application of systems approaches to P. multocida is still in the initial stages, but we have started using these methods to characterize the P. multocida response to antibiotics at the proteome level [11]. Here we describe the use of systems analysis at the transcriptome level to investigate the P. multocida response to sub-MICs of antibiotics.

Results and discussion

Differentially expressed genes

Statistical analysis of microarray data revealed that 1/4 MIC of AMX, CTC, and ENR resulted in significant changes in expression of approximately 25% of the genome (525 genes, Additional file 1, supplemental table 1). The differences in gene expression that were determined to be statistically significant ranged from as small as a 5% decrease (relative to control) to as high as an 11.5-fold increase in expression (recN, Table 1). Earlier microarray studies with P. multocida arbitrarily chose a 1.5-fold significant change in gene expression as differentially regulated genes. Instead of choosing a pre-determined threshold for determining the biological relevance of changes in gene expression, we considered every significant change in expression to be a valid change based on our rigorous statistical testing. We validated changes as small as a 15% decrease or increase (for example, expression of murA in response to AMX and ENR, respectively) (Additional file 1, supplemental table 1) by qPCR. Figure 1 shows the correlation between the trends observed in microarray and qPCR for the subset of genes with qPCR validations.

Additional file 1. Significant changes identified in P. multocida gene expression with AMX, CTC and ENR.

Format: XLS Size: 189KB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

Table 1. Differential expression of P. multocida genes involved in SOS response

thumbnailFigure 1. Validation of microarray differential gene expression data by quantitative real time RT-PCR (qPCR). Ratio of treated vs. control were calculated for microarray data (ordinate) and qPCR data (abscissa). A total of six genes that were up or down regulated with amoxicillin, chlortetracycline or enrofloxacin were compared. The correlation coefficient r2 was 0.84.

Compared to the untreated control, expression of 413, 392, and 473 genes had significantly altered expression in response to AMX, CTC, and ENR, respectively, and the overlap between treatments is shown in the Venn diagram (Figure 2). The extensive overlap of regulated genes with the three different antibiotics is consistent with the recent transcriptome analysis of ampicillin and ofloxacin effects in E. coli, which showed that these two unrelated antibiotics had a significant overlap of regulated genes [14].

thumbnailFigure 2. Venn diagram of gene expression in response to sub-MIC antibiotics. The data intersections between significant changes in gene expression of P. multocida in response to sub-MIC amoxicillin (AMX), chlortetracycline (CTC) and enrofloxacin (ENR).

Functional analysis of DE genes: COGs

The three antibiotics used in this study have distinct modes of action; AMX is a cell wall biosynthesis inhibitor, CTC is a protein synthesis inhibitor, and ENR inhibits DNA gyrase and DNA topoisomerase IV activities [15]. At the COG category level, the three antibiotics had different effects on gene expression (Figure 3). The AMX and CTC treatments resulted in an overall decrease in gene expression in all COG categories. The overall suppression of gene expression could indicate that the antibiotics have a marked detrimental effect on the fitness of P. multocida at doses below MIC, or it is also possible that the overall transcriptional shutdown is a compensatory response by slowing metabolism. ENR had varying effects on gene expression in different COG categories. In particular, ENR caused a pronounced increase in expression of genes in categories I (lipid metabolism), J (translation, ribosomal structure and biogenesis), and L (DNA replication, recombination and repair). All three antibiotics had reduced expression of genes in COGs V (defense mechanisms), G (carbohydrate transport and metabolism), P (inorganic ion transport and metabolism) and Q (secondary metabolites biosynthesis, transport and catabolism).

thumbnailFigure 3. Summary of significant changes in P. multocida gene expression grouped by COGs. Numbers on the y-axis represent the percentage of genes in each COG category whose expression either significantly increased (black bars) or decreased (gray bars) in response to antibiotic treatment. COG category descriptions are: A, RNA processing and modification; K, Transcription; L, Replication, recombination and repair; D, Cell cycle control, mitosis and meiosis; V, Defense mechanisms; T, Signal transduction mechanisms; M, Cell wall/membrane biogenesis; N, Cell motility; W, Extracellular structures; U, Intracellular trafficking and secretion; O, Posttranslational modification, protein turnover, chaperones; C, Energy production and conversion; G, Carbohydrate transport and metabolism; E, Amino acid transport and metabolism; F, Nucleotide transport and metabolism; H, Coenzyme transport and metabolism; I, Lipid transport and metabolism P, Inorganic ion transport and metabolism; Q, Secondary metabolites biosynthesis, transport and catabolism R, General function prediction only; S, Function unknown; -, not in COGs

Effects of antibiotics on known target genes

We observed altered expression of some of the known targets (direct or inferred based on mechanism of action) with each antibiotic. CTC reduced the expression of genes encoding ribosomal proteins S4 and S7, which interact with tetracyclines [16]. One quarter-MIC of ENR resulted in a 1.5 fold increase in the expression of gyrB, which encodes a known ENR target, DNA gyrase B (Table 2). We had previously reported that the expression of RecA protein increased in response to sub-MIC of ENR [10], and in the current study we detected increased expression of recA. With AMX, except for the expression of ponC, there were no detectable differences in the expression of genes encoding penicillin binding proteins (pbp). There was a 50% decrease in ponC (pbp) gene expression. However, this was a not an AMX specific effect, as all three antibiotics caused a decrease in ponC gene expression.

Table 2. Changes in specific target gene expression in response to sub-MIC of antibiotics

While regulation of expression of known targets is an expected response to antibiotics, it is not a consistent finding; expression of the target is usually not affected. For example, microarray analysis showed that enrofloxacin, trimethoprim, brodimoprim, and cefquinome had no effects on expression of their targets in P. multocida [12]. However, even though expression of the known antibiotic target is often not affected, bacteria typically have a "signature" gene expression response to specific antibiotics. Indeed, bacterial transcription profiles are now sometimes used to obtain initial indications of the mechanism of action for new compounds [17,18].

Common trends in response to AMX, CTC and ENR

Despite differences in their mechanism of action, the majority of genes with significantly altered expression (282) were common to all three antibiotics (Figure 2). Beyond the effects on specific target gene expression, antibiotics are known to cause secondary effects on genes involved in general physiology as part of the adaptive response to antibiotic stress. A hallmark of recently described antibiotic mediated bacterial cell death common to bactericidal antibiotics in E. coli is the generation of hydroxyl radicals with subsequent induction of SOS response [19,20]. These effects were noted for MIC as well as sub-MIC concentrations of the antibiotics. There is increasing appreciation for the fact that fine tuning of these responses is unique to each organism and needs to be evaluated as such to facilitate identification of targets that potentiate the bactericidal activity of antibiotics. The P. multocida genome is known to contain genes necessary for antibiotic mediated SOS response [6] (for example, sulA, which is necessary for beta-lactam mediated SOS response).

Antibiotic stress in P. multocida resulted in increased ATP synthesis (Additional file 2, supplemental table 2). While this effect has been described for gyrase inhibitors [21], in P. multocida this response was common to all three antibiotics. In particular, all three antibiotics influenced gene expression involved in de novo nucleotide biosynthesis. The overall effects on expression of purM (increase), purK (decrease), and purF (decrease) could reduce the intracellular concentrations of IMP, while increased expression of pyrG and pyrF with ENR could result in increased concentrations of CTP and UDP.

Additional file 2. Supplemental table 2. Significant changes identified in P. multocida gene expression common to more than one antibiotic.

Format: DOC Size: 117KB Download file

This file can be viewed with: Microsoft Word ViewerOpen Data

Both AMX and CTC impaired the expression of hsp90, thus affecting proper protein folding under stress conditions in P. multocida. All three classes of antibiotics adversely affected putrescine (potE), ribose (rbscA and rbscC), and molybdate (modC) transport systems. Alterations to modC expression could have wide ranging effects on bacterial physiology due to effects on the synthesis of molybdoenzymes (Additional file 2, supplemental table 2). Reduced levels of polyamines like putrescine could have profound negative effects on protein translation and impair the ability to cope with oxidative stress [22].

Expression of 121 genes annotated as coding for hypothetical proteins were altered in response to the antibiotics (Additional file 1, supplemental table 1). Lack of functional information for these genes hampers our understanding of their role in adaptation of P. multocida to antibiotic stress and demonstrates the need for continuous functional annotation of genes beyond the initial annotation that is described with the genome sequence.

One of the demonstrated outcomes of antibiotic stress in E. coli is alterations in iron homeostasis; namely, expression of genes involved in siderophore mediated iron transport are affected (for example, fecC) [19]. In P. multocida, expression of genes encoding iron transporters fecC and fecE were reduced when it was grown in the presence of 1/4 MIC of AMX or ENR (Additional file 2, supplemental table 2).

Changes in the ratio of NAD+/NADH are known to influence the levels of toxic superoxide formation, which ultimately leads to hydroxyl radical formation. NADH I levels are directly linked to the activities of TCA cycle enzymes that generate NADH from NAD+ [19]. ENR treatment caused decreased expression of sucB, and all three antibiotics caused decreased lpdA expression, which could reduce hydroxyl radical formation. Expression of the gene encoding anaerobic respiration electron transfer enzyme molybdoenzyme dimethylsulfoxide reductase was reduced in response to AMX and ENR. AMX treatment reduced the expression of acetyl-CoA carboxylase enzyme (accB, accD), which catalyzes the first step of type II fatty acid biosynthesis as well as the step catalyzed by acyl carrier protein S-malonyltransferase (fabD); by contrast, ENR increased the expression of these genes, which could result in increased fatty acid biosynthesis.

DE genes: network analysis

COG analysis of DE genes was useful for identifying overall trends at a broad functional category level. Gene ontology (GO) analysis of P. multocida, which is available through Uniprot, was useful to some extent in adding specific functions to the genes within each of the broad functional categories represented by COGs. However, neither of these two approaches could actually show exactly how the genes within a functional category were interacting with each other to bring about a specific function. To enable network analysis, we used Pathway Studio to build interaction networks for all DE P. multocida genes with AMX, CTC and ENR. Careful analysis of these networks revealed expression cascades for some of the DE genes common to CTC and ENR.

Expression of fnr increased in response to sub-MIC of CTC and ENR (Figure 4). Fnr regulates the transcription of many genes that are involved in either aerobic or anaerobic respiration, including Nap nitrate reductase, which facilitates nitrate reduction in the periplasmic space. Expression of napD and napF increased in response to ENR and decreased with CTC (Additional file 3, supplemental figure 1). NapD is a chaperone protein that is specific for NapA and is required for the stability of NapA protein. NapF stimulates the Nap reductase. NapF promoters regulate the transcription of cmABCDEFGH operon in E. coli, which constitutes the cytochrome biogenesis system. Expression of ccmH increased with CTC and decreased with ENR (Table 3). The identification of this cascade of gene regulation was possible by combining the network analysis with the functional information available at EcoCyc. The implications of up regulation of enzymes that facilitate anaerobic respiration in response to antibiotics remain unclear at present.

Additional file 3. Supplemental figures. Interaction networks. The network diagrams show changes in gene expression in response to sub-MIC. Red and green indicate significant increase and decrease in expression respectively. Genes whose expression did not change (pink) and genes in the network from other bacteria (gray) are also included. Supplemental figure 1. fnr regulated genes with CTC administration. Supplemental figure 2. Interaction network with genes involved in SOS response with AMX administration. Supplemental figure 3. Interaction network with genes involved in SOS response with CTC administration.

Format: PPT Size: 2MB Download file

This file can be viewed with: Microsoft PowerPoint ViewerOpen Data

Table 3. Fnr regulated genes and their response to sub-MIC of chlortetracycline and enrofloxacin

thumbnailFigure 4. Fumarate/nitrate response regulator network in P. multocida in response to enrofloxacin visualized in Pathway Studio. P. multocida sub-MIC ENR response was marked by significant changes in fnr response regulated genes. Red nodes are genes with increased expression and green nodes are genes with decreased expression. The rest of the interacting nodes are shown in gray. These nodes had no significant changes in expression or were orthologs from additional gram-negative species in the molecular interaction database in Pathway Studio.

Network analysis also facilitated the identification of differences in the SOS response with different antibiotics in P. multocida (Figure 5 (ENR) and Additional file 3, supplemental figures 2 (AMX) and 3 (CTC)). In particular, the response to ENR showed the signature gene expression pattern typical of SOS response (Figure 5). Increased RecA protein expression results in the cleavage of LexA repressor, thus removing the LexA mediated repression of DNA repair genes like uvrD. Expression of both recA and lexA genes increased in response to ENR (Table 1). Expression of recN was increased, which also relieves LexA repression [23-25]. Genes involved in fine-tuning the actions of RecA, namely recX and uvrD, also had increased expression. Expression of primosome genes (dnaJ and dnaK) that encode enzymes required for re-starting the stalled replication fork and other replicative enzymes like DNA polymerase III increased with ENR.

thumbnailFigure 5. SOS response gene network in P. multocida in response to enrofloxacin. P. multocida sub-MIC ENR response was marked by significant changes in SOS response genes. Red nodes are genes with increased expression and green nodes are genes with decreased expression. Pink nodes represent no significant changes in gene expression while gray nodes are orthologs from gram-negative species in the molecular interaction database in Pathway Studio.

It is reported that beta lactam mediated cell killing involves SOS response through SulA and DpiBA two component system [26]. Based on the genome annotation of Pm70 [6], all the components of this regulon are not present. There was no induction of the expression of recA or lexA with AMX (Additional file 3, supplemental figure 2). The expression profile of SOS response genes remained either unchanged or showed decreased expression, which indicates that sub-MIC AMX does not activate SOS response in P. multocida. Excluding the increased expression of lexA and gyrB, CTC (Additional file 3, supplemental figure 3) also did not activate SOS response.

DE virulence genes

Productive infection by bacterial pathogens relies on the expression of virulence factors that have wide ranging functions like competence, adherence, capsule synthesis and export, evading host immune responses etc. Transcription profiling of the response of a bovine P. multocida isolate (L386) to MIC of eight different antibiotics identified mostly reduced virulence gene expression [12]. In the current study, transcription profiling of the P. multocida response to AMX, CTC and ENR identified significant changes in the expression of known and putative virulence genes. AMX and CTC reduced the expression of ompA, a known virulence factor of P. multocida involved in binding to host cells [27] (Table 4). AMX and CTC decreased expression of the gene encoding detoxifying enzyme superoxide dismutase, which is one of the virulence genes identified in the P. multocida L386 study [12]. Additional virulence factors (based on L386 study) which had reduced expression in Pm70 are: clpB protease with AMX, capsule transport protein hexA with all three antibiotics, PM1714 regulator with AMX and CTC, capsule biosynthesis genes phyA and hyaE with AMX and CTC, phyB with CTC, and tight adherence gene tadB with AMX. Expression of dnaK, which is required for virulence in a number of bacterial pathogens [26], was decreased by CTC.

Table 4. Effects of sub-MIC of antibiotics on P. multocida putative virulence factors

Conclusion

Our global transcriptome analysis of P. multocida response to three antibiotics with differing modes of action identified gene expression changes in a little over a quarter of the annotated open reading frames in the genome. The antibiotics had varying effects on various cellular and metabolic functions. Amoxicillin reduced the overall transcription rate as reflected by the majority of the identified significant changes in gene expression showing a downward trend in expression. An interesting aspect of the AMX response was the marked lack of expression of SOS response genes. Enrofloxacin, a DNA gyrase inhibitor, resulted in the significant overexpression of its target gene and induced a typical SOS response. Analysis of DE genes in the context of interacting protein networks facilitated the identification of coordinated regulation of gene expression across different COG functional categories. Sub-MIC of antibiotics influenced the expression of virulence factors in P. multocida. Our results form the framework for understanding the global effects of antibiotics on P. multocida, which will aid rational drug design for containing as well as treating infections caused by this pathogen across multiple species.

Methods

Bacterial culture and RNA isolation

P. multocida Pm70 was cultivated in brain heart infusion broth (BHI) at 37°C with rotary aeration. Minimum inhibitory concentrations (MICs) of AMX, CTC and ENR for Pm70 are 0.5 μg/ml, 4 μg/ml, and 0.031 μg/ml, respectively [10]. Growth kinetics of Pm70 in the presence of 1/4 MIC of the three antibiotics were previously described [10]. Stationary phase cultures of Pm70 were used to inoculate 50 ml BHI to an initial A600 of 0.05; antibiotic treated cultures contained 1/4 MIC of AMX, CTC, or ENR, and control cultures were grown without antibiotics. For each treatment, cultures were grown in triplicate to mid-log phase (A600 of 0.8). RNAprotect (Qiagen) was added to culture samples, and bacteria were harvested by centrifugation (10,000 × g, 10 min, 4°C). Pellets were stored at -80°C.

Total RNA from three biological replicate cultures (for untreated control culture as well as antibiotic treated cultures) was isolated using Qiagen RNeasy kit using the manufacturer's protocol. The quality and concentration of RNA was determined by Agilent Bioanalyzer.

Microarrays and data analysis

A custom oligonucleotide microarray consisting of all the 2015 annotated open reading frames (ORFs) of Pasteurella multocida strain Pm70 was constructed by NimbleGen Systems. The array design included 23 perfect match (PM) and mismatch (MM) probes per target gene, and each probe was spotted twice. Reverse transcription, cDNA labeling, and hybridization were conducted by NimbleGen Sytems according to their established protocols. Briefly, three independent RNA samples (10 μg each) from untreated and antibiotic treated Pm70 were reverse transcribed using Random Hexamer (Gibco) and SuperScript II Reverse Transcriptase (Invitrogen). Biotin end-labeling of cDNA utilized Terminal Deoxynucleotidyl Transferase (Promega) with Biotin-N6-ddATP (Perkin Elmer); hybridization and washes were conducted using NimbleGen Systems protocol. Detection of hybridized biotin-labeled probe used Fluorolink Cy3 Labeled Streptavidin (Amersham Pharmacia).

Raw intensity data (single channel) were log transformed, and normalization was done by quantile normalization [28] and robust multichip analysis (RMA) algorithm [29]. The experimental design and all microarray data have been deposited in the NCBI Gene Expression Omnibus (GEO, http://www.ncbi.nlm.nih.gov/geo/ webcite, accession number GSE12779). ANOVA followed by step-down Bonferroni correction for multiple testing was conducted to identify genes with significant changes in expression relative to untreated control (p ≤ 0.05) in any group. Significant changes between control vs specific antibiotic treatments were identified by performing Dunnett's test [30]. The differences in gene expression in response to antibiotics were calculated as the ratio of treatment vs control intensities.

Real time RT-PCR

The majority of significant changes in gene expression identified in this study were less than one-fold. Therefore, to validate changes in gene expression, we performed duplex real time quantitative RT-PCR (qPCR) of six genes: gyrb, impA, rpsg, hpkR, recQ, and murA using gapdh as the internal standard. The RNA template used for qPCR was the same RNA that was used for microarray hybridizations. The probes and primers for duplex RT-PCR were designed using Beacon Designer software (Additional file 4, supplemental table 3). Reactions were performed using the Platinum® SYBR® Green One-Step qRT-PCR Kit (Invitrogen Corporation, Carlsbad, CA). Amplification and detection of specific products were done using the iCycler iQ Real-Time PCR Detection System (Bio-Rad Laboratories, Inc., Hercules, CA). Regression analysis of the Ct values of the test RT-PCRs was used to generate standard curves. The mean Ct value for the GAPDH mRNA-specific reactions was used to normalize the test values. Regression analysis of the gene expression trends determined by qPCR and microarray analysis was performed in Microsoft Excel.

Additional file 4. Supplemental table 3. Quantitative RT-PCR probes and primers. The sequences of primers used for qPCR validations of selected gene expression changes in microarray.

Format: DOC Size: 25KB Download file

This file can be viewed with: Microsoft Word ViewerOpen Data

Functional analysis of differentially expressed genes

Initial functional analysis of genes with significant changes in expression was done using clusters of orthologous groups (COG) classification [31]. For each COG category in the genome, the percent increase or decrease in expression of genes belonging to that category for each antibiotic treatment was calculated. The COG functional categories described for P. multocida genes [6] along with the description for proteins encoded by the genes at Uniprot [32] were used to identify overall themes in antibiotic response. Specific functions of P. multocida proteins were deduced from information available for E. coli orthologs at the EcoCyc database [33]. Network analysis of DE genes was done in Pathway Studio (Ariadne Genomics) as described earlier [11]. Briefly, we built interaction networks in Pathway Studio with proteins of interest including the upstream regulators and downstream targets. In the interaction networks, different colors were used for nodes to indicate significant increase (red), decrease (green), or no change (pink) in gene expression in response to sub-MIC of antibiotic. Entities in the interaction map that were not present in the Pm70 genome were shown in gray color.

List of abbreviations used

MIC: minimum inhibitory concentration; CTC: chlortetracycline; AMX: amoxicillin; ENR: enrofloxacin; BHI: brain heart infusion; COG: Clusters of Orthologous Groups; DE: differentially expressed; recA: recombinase A; gyrB: DNA gyrase B.

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

BN cultured Pm70, did RNA isolations, data analysis and network modeling and wrote the draft of the manuscript. LAS carried out the qPCR validations and contributed to the methods and results of this aspect in the manuscript. ML and SCB conceived and coordinated this study and helped draft the manuscript. All authors read and approved the final manuscript.

Acknowledgements

This project was partially supported by a grant from the National Science Foundation (EPS-0556308-06040293) and a competitive grant from Mississippi State University Life Sciences and Biotechnology Institute (LSBI). We acknowledge Dr. Vivek Kapur for P. multocida Pm70 strain and Dr. Carolyn Boyle for statistical analysis of microarray datasets.

This article has been published as part of BMC Genomics Volume 10 Supplement 2, 2009: Proceedings of the Avian Genomics Conference and Gene Ontology Annotation Workshop. The full contents of the supplement are available online at http://www.biomedcentral.com/1471-2164/10?issue=S2

References

  1. Pasteur L: De l'att'enuation du virus du chol'era des poules.

    C R Acad Sci 1880, 91:673-680. OpenURL

  2. Pasteur L: Sur les virus-vaccins du chol'era des poules et du charbon.

    C R Trav Congr Int Dir Agronom session de Versailles 1881, 161-162. OpenURL

  3. Ellison ML, Champlin FR: Outer membrane permeability for nonpolar antimicrobial agents underlies extreme susceptibility of Pasteurella multocida to the hydrophobic biocide triclosan.

    Vet Microbiol 2007, 124(3–4):310-318. PubMed Abstract | Publisher Full Text OpenURL

  4. Harper M, Boyce JD, Adler B: Pasteurella multocida pathogenesis: 125 years after Pasteur.

    FEMS Microbiol Lett 2006, 265(1):1-10. PubMed Abstract | Publisher Full Text OpenURL

  5. Talan DA, Citron DM, Abrahamian FM, Moran GJ, Goldstein EJ: Bacteriologic analysis of infected dog and cat bites. Emergency Medicine Animal Bite Infection Study Group.

    N Engl J Med 1999, 340(2):85-92. PubMed Abstract | Publisher Full Text OpenURL

  6. May BJ, Zhang Q, Li LL, Paustian ML, Whittam TS, Kapur V: Complete genomic sequence of Pasteurella multocida, Pm70.

    Proc Natl Acad Sci USA 2001, 98(6):3460-3465. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  7. Boyce JD, Wilkie I, Harper M, Paustian ML, Kapur V, Adler B: Genomic scale analysis of Pasteurella multocida gene expression during growth within the natural chicken host.

    Infect Immun 2002, 70(12):6871-6879. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  8. Paustian ML, May BJ, Cao D, Boley D, Kapur V: Transcriptional response of Pasteurella multocida to defined iron sources.

    J Bacteriol 2002, 184(23):6714-6720. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  9. Paustian ML, May BJ, Kapur V: Transcriptional response of Pasteurella multocida to nutrient limitation.

    J Bacteriol 2002, 184(13):3734-3739. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  10. Nanduri B, Lawrence ML, Boyle CR, Ramkumar M, Burgess SC: Effects of subminimum inhibitory concentrations of antibiotics on the Pasteurella multocida proteome.

    J Proteome Res 2006, 5(3):572-580. PubMed Abstract | Publisher Full Text OpenURL

  11. Nanduri B, Lawrence ML, Peddinti DS, Burgess SC: Effects of Subminimum Inhibitory Concentrations of Antibiotics on the Pasteurella multocida Proteome: A Systems Approach.

    Comp Funct Genomics 2008, 254836. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  12. Melnikow E, Schoenfeld C, Spehr V, Warrass R, Gunkel N, Duszenko M, Selzer PM, Ullrich HJ: A compendium of antibiotic-induced transcription profiles reveals broad regulation of Pasteurella multocida virulence genes.

    Vet Microbiol 2008. PubMed Abstract | Publisher Full Text OpenURL

  13. Weston AD, Hood L: Systems biology, proteomics, and the future of health care: toward predictive, preventative, and personalized medicine.

    J Proteome Res 2004, 3(2):179-196. PubMed Abstract | Publisher Full Text OpenURL

  14. Kaldalu N, Mei R, Lewis K: Killing by ampicillin and ofloxacin induces overlapping changes in Escherichia coli transcription profile.

    Antimicrobial agents and chemotherapy 2004, 48(3):890-896. Publisher Full Text OpenURL

  15. Walsh C: Antibiotics: actions, origins, resistance. ASM press; 2003.

  16. Yonath A, Bashan A: Ribosomal crystallography: initiation, peptide bond formation, and amino acid polymerization are hampered by antibiotics.

    Annu Rev Microbiol 2004, 58:233-251. PubMed Abstract | Publisher Full Text OpenURL

  17. Freiberg C, Fischer HP, Brunner NA: Discovering the mechanism of action of novel antibacterial agents through transcriptional profiling of conditional mutants.

    Antimicrobial agents and chemotherapy 2005, 49(2):749-759. Publisher Full Text OpenURL

  18. Hutter B, Schaab C, Albrecht S, Borgmann M, Brunner NA, Freiberg C, Ziegelbauer K, Rock CO, Ivanov I, Loferer H: Prediction of mechanisms of action of antibacterial compounds by gene expression profiling.

    Antimicrobial agents and chemotherapy 2004, 48(8):2838-2844. Publisher Full Text OpenURL

  19. Dwyer DJ, Kohanski MA, Hayete B, Collins JJ: Gyrase inhibitors induce an oxidative damage cellular death pathway in Escherichia coli.

    Mol Syst Biol 2007, 3:91. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  20. Kohanski MA, Dwyer DJ, Hayete B, Lawrence CA, Collins JJ: A common mechanism of cellular death induced by bactericidal antibiotics.

    Cell 2007, 130(5):797-810. PubMed Abstract | Publisher Full Text OpenURL

  21. Dahan-Grobgeld E, Livneh Z, Maretzek AF, Polak-Charcon S, Eichenbaum Z, Degani H: Reversible induction of ATP synthesis by DNA damage and repair in Escherichia coli. In vivo NMR studies.

    J Biol Chem 1998, 273(46):30232-30238. PubMed Abstract | Publisher Full Text OpenURL

  22. Shah P, Swiatlo E: A multifaceted role for polyamines in bacterial pathogens.

    Mol Microbiol 2008, 68(1):4-16. PubMed Abstract | Publisher Full Text OpenURL

  23. Erill I, Campoy S, Barbe J: Aeons of distress: an evolutionary perspective on the bacterial SOS response.

    FEMS Microbiol Rev 2007, 31(6):637-656. PubMed Abstract | Publisher Full Text OpenURL

  24. Kelley WL: Lex marks the spot: the virulent side of SOS and a closer look at the LexA regulon.

    Mol Microbiol 2006, 62(5):1228-1238. PubMed Abstract | Publisher Full Text OpenURL

  25. Cox MM: Regulation of bacterial RecA protein function.

    Crit Rev Biochem Mol Biol 2007, 42(1):41-63. PubMed Abstract | Publisher Full Text OpenURL

  26. Miller C, Thomsen LE, Gaggero C, Mosseri R, Ingmer H, Cohen SN: SOS response induction by beta-lactams and bacterial defense against antibiotic lethality.

    Science 2004, 305(5690):1629-1631. PubMed Abstract | Publisher Full Text OpenURL

  27. Dabo SM, Confer AW, Quijano-Blas RA: Molecular and immunological characterization of Pasteurella multocida serotype A:3 OmpA: evidence of its role in P. multocida interaction with extracellular matrix molecules.

    Microb Pathog 2003, 35(4):147-157. PubMed Abstract | Publisher Full Text OpenURL

  28. Bolstad BM, Irizarry RA, Astrand M, Speed TP: A comparison of normalization methods for high density oligonucleotide array data based on variance and bias.

    Bioinformatics 2003, 19(2):185-193. PubMed Abstract | Publisher Full Text OpenURL

  29. Irizarry RA, Hobbs B, Collin F, Beazer-Barclay YD, Antonellis KJ, Scherf U, Speed TP: Exploration, normalization, and summaries of high density oligonucleotide array probe level data.

    Biostatistics 2003, 4(2):249-264. PubMed Abstract | Publisher Full Text OpenURL

  30. Shaffer JP: Multiple comparisons emphasizing selected contrasts: an extension and generalization of Dunnett's procedure.

    Biometrics 1977, 33(2):293-303. PubMed Abstract | Publisher Full Text OpenURL

  31. Tatusov RL, Koonin EV, Lipman DJ: A genomic perspective on protein families.

    Science 1997, 278(5338):631-637. PubMed Abstract | Publisher Full Text OpenURL

  32. The universal protein resource (UniProt)

    Nucleic Acids Res 2008, (36 Database):D190-195. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  33. Keseler IM, Collado-Vides J, Gama-Castro S, Ingraham J, Paley S, Paulsen IT, Peralta-Gil M, Karp PD: EcoCyc: a comprehensive database resource for Escherichia coli.

    Nucleic Acids Res 2005, (33 Database):D334-337. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL