Previous Article | Next Article ![]()
Infection and Immunity, December 2002, p. 6871-6879, Vol. 70, No. 12
0019-9567/02/$04.00+0 DOI: 10.1128/IAI.70.12.6871-6879.2002
Copyright © 2002, American Society for Microbiology. All Rights Reserved.
Bacterial Pathogenesis Research Group, Department of Microbiology, School of Biomedical Sciences, Monash University, Victoria 3800,1 Veterinary Pathology and Anatomy, School of Veterinary Science and Animal Production, The University of Queensland, Queensland 4072, Australia,2 Department of Veterinary Pathobiology and Biomedical Genomics Center, University of Minnesota, St. Paul, Minnesota 551083
Received 24 June 2002/ Returned for modification 23 August 2002/ Accepted 13 September 2002
|
|
|---|
|
|
|---|
Several methods have been devised to identify genes important or essential for growth within a host. Two of the most powerful methods are in vivo expression technology (IVET) (26) and signature-tagged mutagenesis (STM) (14), both of which have been used to identify Pasteurella multocida virulence genes (9, 17). Recently, the development of DNA microarray methods for studying gene expression on a whole-genome scale has provided the opportunity to analyze gene expression changes directly in response to growth within a host. There is a large amount of published data on genes that are differentially regulated in vitro under conditions that mimic those within host organisms. Such conditions include low iron (27, 31) and growth at different temperatures (38, 41), at different pHs (1, 38), and on different media (32, 43). However, at present, the only large-scale transcription-profiling data available on bacterial growth within a natural host organism are those on the growth of Borrelia burgdorferi within dialysis membranes implanted into rat peritoneal cavities (38). Furthermore, there are only a very limited number of data available on how the transcription of even a small number of pathogen genes responds during growth within host organisms (47, 48).
P. multocida is a gram-negative bacterial pathogen that causes a range of diseases in mammals and birds. It is the etiological agent of several economically important diseases, including fowl cholera, atrophic rhinitis in pigs, snuffles in rabbits, and hemorrhagic septicemia in cattle (24). Despite considerable research, the molecular mechanisms by which P. multocida can survive and multiply within a host are poorly understood. Indeed, only a small number of true bacterial virulence or virulence-associated genes (46) have been definitively identified. The identified true virulence genes of P. multocida include those involved in the production of toxins (restricted to a few toxin-producing strains) (22), capsules (2), and hemagglutinins and hemolysins (9), while the virulence-associated genes include those involved in amino acid, nucleotide, and iron transport and metabolism (8, 9, 15). Thus, it is likely that many important P. multocida virulence genes remain unidentified and uncharacterized.
We are interested in identifying P. multocida virulence genes in order to understand pathogenesis at the molecular level and also because of the potential for vaccine development. To this end, we have undertaken a genomic-scale comparison of P. multocida gene expression during growth in rich medium and growth within the chicken host. We suggest that genes expressed at higher levels during growth within a natural host are likely virulence genes. Therefore, we believe that the genes identified by this analysis will form the basis for future directed vaccine approaches as targets both for attenuating mutations and for producing recombinant antigens. Furthermore, the information will substantially increase our understanding of the fundamental interaction between this bacterial pathogen and its host and provide a broad framework for addressing how bacterial transcription is regulated during the course of infection.
|
|
|---|
RNA isolation. Bacteria were harvested from duplicate BHI cultures at late log phase (5 x 109 CFU/ml), added to 0.1 volume of ice-cold killing buffer (0.05 M Tris-HCl [pH 7.5], 15 mg of sodium azide/ml, 0.6 mg of chloramphenicol/ml), and pelleted by centrifugation. RNA was isolated from bacteria by using Trizol reagent (Gibco/BRL) as described by the manufacturer. Purified RNA was treated with DNase (15 U for 10 min at 37°C), and the RNA was further purified on RNeasy minicolumns (Qiagen). For the isolation of in vivo-grown bacteria, chickens (outbred Leghorn cross commercial layers) were infected with 5 x 104 CFU of P. multocida strain X-73 by injection into the breast muscle. Blood was recovered from infected chickens during the final stages of disease (between 17 and 22 h after infection), when the level of bacteremia was observed to be between 109 and 1010 CFU/ml. Clinical signs at the time of blood recovery indicated that each infection was in the terminal phase. Blood (30 to 40 ml) was recovered by terminal heart puncture and added to 0.1 volume of ice-cold killing buffer, and the bacteria were separated from the red blood cells by sucrose density centrifugation for 1 h at 3,000 x g in the presence of 12.5% (wt/vol) sucrose. Bacteria were removed from above the sucrose cushion and concentrated by centrifugation for 5 min at 12,000 x g. All manipulations were carried out on ice or at 4°C and took less than 1.5 h. RNA was purified as described above for in vitro-grown bacteria.
Microarray hybridizations. P. multocida strain PM70 arrays were constructed as described previously (32). These arrays were generated by PCR amplification of each of the 2,015 annotated P. multocida genes and contained triplicate spots for each gene. Fluorescently labeled cDNA was prepared by the incorporation of amino allyl dUTP (Sigma) during reverse transcription (RT) and the subsequent coupling of monofunctional cy3 and cy5 dyes (Amersham) to the amino allyl-labeled cDNA. RT reactions contained 20 µg of total RNA; 30 µg of random hexamers; 10 U of Superscript II (Gibco/BRL) reverse transcriptase; 500 µM concentrations each of dATP, dCTP, and dGTP; 200 µM dTTP; and 300 µM amino allyl dUTP. The amino allyl-labeled cDNA was coupled to either the cy3 or cy5 dyes as described previously (41), and experimental and control probes were hybridized to the same microarray. Typically, 10 µg of each labeled cDNA probe (cy3- and cy5-labeled cDNA) was used in a single hybridization. Microarray hybridizations were carried out at 42°C for 16 h in a solution of 25% formamide, 5x SSC (1x SSC is 0.15 M NaCl plus 0.015 M sodium citrate), 0.1% sodium dodecyl sulfate (SDS), 0.5 mg of yeast tRNA/ml, and 1 mg of sheared salmon sperm DNA/ml under a coverslip in a humidified chamber. The slides were given one 5-min wash in 2x SSC-0.1% SDS at 42°C, one 10-min wash in 0.1x SSC-0.1% SDS at 20°C, and four 1-min washes in 0.1% SDS at 20°C and then rinsed in distilled H2O and ethanol. Slides were dried by centrifugation at 1,000 x g for 3 min and scanned immediately.
Real-time RT-PCR. Primers for real-time RT-PCR were designed with Primer Express software (ABI) (Table 1). RT reactions were carried out as described for microarray hybridizations (see above). The synthesized cDNA samples were diluted 40- to 80-fold prior to real-time RT-PCR, which was carried out on an ABI PRISM model 7700 sequence detector with product accumulation quantified by incorporation of the fluorescent dye SYBR Green. Triplicate real-time RT-PCRs were performed on 2.4 µl of cDNA with SYBR Green PCR master mix (ABI) and 50 nM concentrations of each gene-specific primer. The concentration of template in each reaction was determined by comparison with a gene-specific standard curve constructed from known concentrations of P. multocida strain X-73 genomic DNA. gyrB was used as the normalizer for all reactions. All RT-PCRs amplified a single product as determined by melting curve analysis.
|
View this table: [in a new window] |
TABLE 1. Sequences of primers used for real-time RT-PCR
|
|
|
|---|
Microarrays were hybridized with control and experimental RNA labeled with either cy3 or cy5. Two hybridizations were carried out with each RNA sample; one hybridization had the experimental sample labeled with cy3, and the other hybridization had the experimental sample labeled with cy5. Hybridizations derived from the different chicken infections were analyzed separately. We designated all genes with average expression ratios of 1.5-fold (up or down) differentially expressed. We chose the 1.5-fold value, as other studies have indicated that a 1.5-fold change in gene expression may be biologically relevant (16, 41). Furthermore, we used the Genesight (Biodiscovery) confidence analyzer to identify genes that had a 95% likelihood of being differentially expressed at above or below the 1.5-fold threshold.
Real-time RT-PCR was used to confirm a number of the expression ratios identified by microarray analysis (Fig. 1). Thirteen genes (including up- and down-regulated and constitutively expressed genes) were analyzed with bacterial RNA from each of the three infected chickens. For RNA derived from all three infections, there was a strong positive correlation between the two methods (r2 = 0.83). The slope of the regression line (m = 1.25) indicated that generally slightly higher expression changes were measured by real-time RT-PCR rather than by microarray hybridization.
![]() View larger version (13K): [in a new window] |
FIG. 1. Correlation between microarray and real-time RT-PCR expression values. Log2-transformed expression values for 13 genes from bacteria isolated from each of the three infected chickens (c1, c2, and c3) were compared. The best-fit linear regression line is shown together with the r2 value and the calculated equation.
|
![]() View larger version (15K): [in a new window] |
FIG. 2. Venn diagrams showing the number of P. multocida genes up-regulated (A) or down-regulated (B) in bacterial populations isolated from three infected chickens (c1, c2, and c3).
|
|
View this table: [in a new window] |
TABLE 2. P. multocida genes up-regulated during growth within all three infected chickens
|
|
View this table: [in a new window] |
TABLE 3. P. multocida genes down-regulated during growth within all three infected chickens
|
The glutamate dehydrogenase gene (gdhA) was very strongly up-regulated in bacteria isolated from all three infected chickens (8-fold, 3-fold, and 10-fold [Table 2]). GdhA converts ammonia from the environment to glutamate (NH3 +
-ketoglutarate + NADPH
glutamate + NADP). In Escherichia coli, this reaction is one of the primary pathways for the assimilation of nitrogen and the glutamate derived from this process is required directly or indirectly for the synthesis of all amino acids; half the nitrogen atoms for the pyrimidine, purine, and imidazole ring; and the amino group of adenine (36). Although glutamate can be synthesized by two different pathways (13), the GdhA pathway is preferred under energy-limited conditions (12, 13). Therefore, in an environment containing ammonia but poor in amino acids and energy limited, this reaction is of critical importance. Glutamate is also a prime cellular defense against osmotic stress and can serve as a carbon source (28). The gene asnA was also strongly up-regulated during growth in all three chickens; this gene encodes an ammonia-dependent asparagine synthetase that produces asparagine from aspartate and ammonia. It has been shown previously in Klebsiella aerogenes that levels of this enzyme are high when cells are grown in an amino acid-limited but ammonia-rich environment (37). The gene aspC, which was also strongly up-regulated in all chickens, is predicted to encode aspartate transaminase, which catalyses the conversion of glutamate to aspartate. In E. coli, this reaction is the prime mechanism of aspartate formation and aspartate is the substrate for the reaction catalyzed by AsnA (see above). It should be noted that glutamate can also function as a major carbon source in E. coli and that the conversion of glutamate to aspartate by AspC is part of the glutamate energy utilization pathway. Furthermore, AspC can function to synthesize phenylalanine and tyrosine (36). Thus, it is likely that the products of gdhA, asnA, and aspC function to convert host-derived ammonia to amino acids. The dppA gene, which was significantly up-regulated in all three chickens, encodes part of a periplasmic dipeptide transport system, while another significantly up-regulated gene, ilvH, encodes an enzyme involved in amino acid biosynthesis. Taken together, these data clearly show that, during growth in the host, P. multocida must modify biosynthetic pathways to cope with an environment that is comparatively amino acid poor but nitrogen (ammonia) rich. To this end, most common up-regulated genes are involved in critical (probably rate-limiting) steps in the transport and biosynthesis of amino acids.
Many of the common up-regulated genes were predicted to be involved in energy production and conversion (Table 2). The nitrate reductase (nap) operon in P. multocida has an arrangement (napFDAGHBC) identical to that found in other gram-negative pathogens (34). The napF, napA, and napC genes were highly up-regulated in the bacteria from all three chickens, while napB was significantly up-regulated in the c2 and c3 chickens and fell just below the level of significance in c1. The napD gene was observed to be significantly up-regulated in the c1 infection but only slightly up-regulated in the other two infections. Taken together, these data indicate that the nitrate reductase operon is significantly up-regulated during growth in the chicken host. The function of the periplasmic nitrate reductase is to reduce nitrate to nitrite, with nitrate being the preferred terminal electron acceptor in the absence of oxygen. Therefore, nitrate reductases are generally assumed to be necessary for anaerobic growth (6). It has also been suggested that an important function of these proteins in bacterial pathogens is to scavenge nitrogen from bodily fluids (34). Given that the nitrite produced by the periplasmic nitrate reductases is rapidly converted to ammonia by nitrite reductases, we hypothesize that this pathway may aid the biosynthesis of glutamate (see above) by increasing the ammonia pool through nitrate scavenging.
Two other genes involved in energy production were up-regulated during growth within all of the three infected chickens (Table 2). Citrate synthase (gltA) has been shown to be the rate-limiting step for the tricarboxylic acid (TCA) cycle when E. coli cells are grown on acetate (45). It is likely that up-regulation of this gene in P. multocida leads to a concomitant increase in the rate of energy production (or biosynthetic functions branching from the TCA cycle) under the carbon-limited conditions encountered during growth within the host. The ppc gene is predicted to encode phosphoenolpyruvate (PEP) carboxylase, which functions to convert PEP to oxaloacetate. Oxaloacetate is the precursor of several intermediates of both the TCA cycle and other biosynthetic pathways.
Seven other genes were identified as significantly up-regulated during growth in all three chickens (Table 2). The purH gene is predicted to be involved in purine biosynthesis. Genes involved in purine or pyrimidine biosynthesis have been identified as virulence-associated genes by STM techniques in a number of bacterial species, including P. multocida (9, 33). The ptfA gene is the type IV fimbrial subunit gene (7, 40). Fimbriae are well-characterized virulence factors (10, 19), and there is evidence that some types of fimbriae may be expressed in blood (35). The pm1460 and pm0092 genes are both likely to be involved in gene regulation. The pm0092 gene is predicted to encode a sigma-54 modulation protein, while pm1460 is predicted to encode a phosphoglycerate transport protein that contains a sigma-54 interaction domain. Genes controlled by sigma-54 are generally involved in nitrogen fixation and utilization, although pilin synthesis is also often controlled by sigma-54 (29). The other three up-regulated genes have either been poorly characterized (pm0601) or have no known function (dcaA and pm0287), although dcaA is likely to be transcribed with purH. Taken together, these data indicate that a clear common theme in the genes up-regulated in bacteria during growth in all three chickens is their central importance in the metabolism of nitrogen-containing compounds.
Genes down-regulated in all infections. In contrast to the genes up-regulated during growth in all three chickens, most (15 of 23) of the genes down-regulated in all three infections (Table 3) were of unknown or poorly characterized function. Two genes, rpL34 and rpL35, encode ribosomal proteins; their lower level of transcription during growth within the host may be an indicator of a reduced bacterial growth rate. The reduced expression of murZ, which encodes UDP-GlcNAc enolpyruvyl transferase, the enzyme which catalyses the first committed step in peptidoglycan synthesis (3), may also indicate reduced growth rate in the chicken host. Two genes, glpB and glpC, encode subunits of the anaerobically expressed glycerol-3-phosphate dehydrogenase, while the gene glpA, which encodes the third subunit, was also significantly down-regulated in two chickens.
Analysis of all differentially expressed genes by functional groupings. Separate analysis of bacteria derived from each of the three chickens indicated that 522 P. multocida genes (26%) were differentially expressed in bacteria isolated from at least one chicken. The expression ratios of these genes (discussed in the following three sections) in each of the three infections are available online (http://www.med.monash.edu.au/microbiology/staff/adler/PastyArray.html). Of these, 270 were up-regulated and 268 were down-regulated during growth in the chicken compared with those under in vitro growth conditions. Sixteen genes were present in both the up-regulated and down-regulated gene sets. The 522 differentially regulated genes were categorized into their clusters-of-orthologous-groups (COGs) functional categories (44) (Fig. 3). COGs categories highly represented in the up-regulated genes included those involved in amino acid, carbohydrate, nucleotide, and inorganic ion transport and metabolism. COGs categories highly represented in the down-regulated genes included those involved in protein synthesis (translation, ribosomal structure, and biogenesis), coenzyme and lipid metabolism, and cell envelope and outer membrane biogenesis. Notably, there were also 169 genes of unknown function identified as differentially regulated during bacterial growth in at least one chicken; 70 of these were up-regulated.
![]() View larger version (22K): [in a new window] |
FIG. 3. Number of differentially expressed genes by COGs functional categories, including all genes identified as differentially expressed in any of the three animals. The bars for category S have been truncated, and the numbers above each bar indicate the number of genes in each category. The COGs functional categories are as follows: J, translation, ribosomal structure, and biogenesis; K, transcription; L, DNA replication, recombination, and repair; D, cell division and chromosome partitioning; O, posttranslational modification, protein turnover, and chaperone functions; M, cell envelope and outer membrane biogenesis; N, cell motility and secretion; P, inorganic ion transport and metabolism; T, signal transduction mechanisms; C, energy production and conversion; G, carbohydrate transport and metabolism; E, amino acid transport and metabolism; F, nucleotide transport and metabolism; H, coenzyme metabolism; I, lipid metabolism; Q, secondary metabolite biosynthesis, transport, and catabolism; S, unknown function (which includes categories R, general function prediction only, and category -, not in COGs).
|
Signal transduction and transcription genes. Many of the genes involved in signal transduction and transcription were differentially regulated in bacteria from at least one infected chicken (Fig. 3). Most of these differentially regulated genes were predicted to be involved in the sensing of the redox state and the presence and/or absence of oxygen and included fnr, oxyR, arcA, and arcB. These findings are consistent with a difference in redox potential between growth in vivo and in aerated broth culture. Surprisingly, the differential regulation of these genes was not consistent across the three infections. Indeed, fnr was significantly up-regulated in vivo in the c2 and c3 infections but significantly down-regulated in the c1 infection. The expression of this gene in the three infections was confirmed by real-time PCR, which gave identical results (data not shown).
Many transcription factors were down-regulated in at least one infection (Fig. 3), including rpoE, rpoH, PM0209, rho, nusA, and nusB. The sigma factors RpoH and RpoE are predicted to be involved in the heat and extracytoplasmic shock responses. Interestingly, rseB, the gene encoding a regulator of RpoE, was also significantly down-regulated during growth in chickens. Paradoxically, the heat shock-induced chaperones HtpG, DnaK, DnaJ, GroES, and PrlC, which are positively regulated by RpoH, were generally up-regulated in vivo. The up-regulation of heat shock proteins during infection has been reported previously (18, 21, 23), and P. multocida has been shown to express heat shock proteins during infection (25). These results may indicate that, in P. multocida, as in E. coli (4), induction of RpoH synthesis is primarily regulated posttranscriptionally. Furthermore, dnaK, dnaJ, and groES transcription is known to be induced by a range of stresses other than heat shock, including oxidative stress, osmotic stress, and amino acid starvation (5, 42).
Differences between infections. Each of the three infections was initiated with the injection of 5 x 104 CFU of P. multocida strain X-73 into the breast muscle. At the time of bacterial harvest, the attained cell density in the blood varied only slightly (between approximately 2 x 109 and 8 x 109 CFU/ml) in birds from which bacterial RNA was taken for study. Furthermore, each of the three chickens displayed identical clinical signs and was in the terminal phase of disease. It is possible that blood homeostasis may have broken down in some of the animals at this time, but we regard this as a natural consequence of infection. Thus, although it is impossible to know if bacteria were harvested at identical stages during the infection, every effort was taken to isolate the bacteria at the same time with regard to disease progression.
Overall, the differentially expressed gene sets identified for bacteria harvested from each chicken were all highly similar (chi-square test; P < 0.001), although there were also clear differences at the level of functional gene subsets. Cluster analysis indicated clearly that the variation within replicate hybridizations of RNA from the same infection was significantly lower than the variation between hybridizations of RNA samples from different infections (data not shown). The differentially regulated genes from each infection were categorized into their COGs functional categories (44) (Fig. 4).
![]() View larger version (26K): [in a new window] |
FIG. 4. The proportion of differentially regulated genes in each COGs category expressed as a percentage of the total number of differentially expressed genes for each infection. The percentages of genes from each infection that were up-regulated (c1 up, c2 up, and c3 up) or down-regulated (c1 down, c2 down, and c3 down) are shown for each COGs category. The COGs functional categories are as described for Fig. 3.
|
The response of P. multocida to low iron concentration in vitro has been studied previously (27, 31), and a number of genes have been identified as being differentially regulated in response to limiting free iron. We compared the set of genes identified as up-regulated in vitro under iron limitation (31), with the set of genes up-regulated in bacteria grown within each of the three infected chickens (Table 4). Only 3 of the 71 genes identified as up-regulated in the c1 infection have previously been shown to respond to low iron in vitro, indicating that these gene sets are not significantly similar (Fisher's exact test; P > 0.9). However, 13 of the 93 genes identified as up-regulated in the c2 infection and 25 of the 180 genes identified as up-regulated in the c3 infection have previously been shown to respond to low iron in vitro, indicating that each of these gene sets is significantly similar (Fisher's exact test; P < 0.001). These data indicate that only the c2 and c3 infections displayed a gene expression profile similar to that seen during in vitro iron starvation. It is often assumed that gene expression in vivo would be similar to gene expression under in vitro conditions in the absence of iron. The gene expression levels in the c1 infection indicate that this may not always be the case, although it is possible that some genes involved in growth during iron-starved conditions are expressed only transiently.
|
View this table: [in a new window] |
TABLE 4. P. multocida genes previously shown to respond to low iron levels in vitro (31) that were up-regulated during growth within infected chickens
|
A number of large-scale screening methods have been used previously to identify potential virulence genes in P. multocida, including STM (9) and IVET (17). In this study, 7 of the 23 P. multocida genes identified by STM were also identified as differentially regulated while 3 of the 14 genes identified by IVET were identified as up-regulated during growth in chickens. Clearly many genes were identified in our study but not by the STM and IVET methods, indicating that neither of these methods produces a saturation coverage of the genome. One of the clear benefits of whole-genome expression profiling with microarrays is the ability to analyze the gene expression profiles of every gene under the conditions being tested. However, some genes were identified by STM and IVET methods that were not identified in this work. This is not surprising, as STM and IVET analyses were carried out on different P. multocida strains and both used mouse infection models. Thus, we believe our infection model is significantly more relevant to natural fowl cholera. Furthermore, STM identifies genes that are necessary for in vivo survival and it is clear that such genes need not be differentially expressed under in vivo and in vitro conditions. Another clear benefit of the microarray approach, compared to that of IVET, is the ability to apply statistical measures of reproducibility to datasets derived from multiple infections.
A key finding of our work was that the expression patterns of bacteria may vary considerably between infected individuals. The reasons for this are not clear but could include genetic host variability or differences in the breakdown of blood homeostasis in animals. However, we suggest that experiments in which in vivo-derived bacteria are pooled from more than one animal may lead to overly simplified conclusions regarding gene expression in the host. An important finding resulting from this work is the identification of 70 genes of unknown function that were significantly up-regulated in at least one infected individual. Future detailed analyses of these genes and their products will lead to an improved understanding of the pathogenesis of fowl cholera at the molecular level.
This work was funded in part by a project grant from the Australian Research Council. M.P. is supported by an NIH National Institute of General Medical Sciences Training for Future Biotechnology Development grant (T32 GM08347). Funding for the preparation of the microarrays was provided by research grants from the Minnesota Turkey Growers Association, the Minnesota Agricultural Experimentation Station, the University of Minnesota Academic Health Center, and the U.S. Department of Agriculture National Research Initiative (to V.K.).
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»