Previous Article | Next Article ![]()
Journal of Bacteriology, April 2004, p. 1959-1971, Vol. 186, No. 7
0021-9193/04/$08.00+0 DOI: 10.1128/JB.186.7.1959-1971.2004
Copyright © 2004, American Society for Microbiology. All Rights Reserved.
Department of Chemical Engineering, Northwestern University, Evanston, Illinois 60208
Received 16 October 2003/ Accepted 23 December 2003
|
|
|---|
|
|
|---|
DNA microarray analysis has become an important tool in elucidating the transcriptional programs of mutant strains. Fawcett et al. examined B. subtilis strains deficient in either spo0A or sigF to identify 520 genes whose expression levels differed at least threefold and were dependent on Spo0A but not
F and 66 genes whose expression depended on both (11). Characterization of SKO1 by DNA microarray analysis reveals downregulation of solventogenic, sporulation, and carbohydrate metabolism genes and upregulation of flagellar and chemotaxis genes (30). The advantage of performing transcriptional analysis of 824(pMSPOA) is that accelerated differentiation may reveal more genes that are differentially expressed in stationary phase. Also, stress response is associated with solvent production; this response cannot be observed in strain SKO1, since it does not produce significant amounts of solvents.
Butanol stress and tolerance have recently been examined in a groESL-overexpressing strain (31, 32). Overexpression of groESL imparts enhanced solvent production (32) and tolerance (31). Butanol toxicity and tolerance have typically been associated with effects on the cell membrane's biophysical properties (31). Since cell differentiation and sporulation involve drastic changes in membrane properties and spore coat formation, we hypothesized that changes in the cell's differentiation program by spo0A overexpression might alter the response of the cells to butanol stress. We examined this hypothesis by subjecting strain 824(pMSPOA) and its plasmid control strain [824(pIMP1)] to butanol challenges and examining the transcriptional differences of these challenges.
|
|
|---|
Analytical methods. Cell growth was determined by measuring absorbance at 600 nm (A600) with a Thermo Spectronic (Rochester, N.Y.) BioMate3 spectrophotometer. Absorbance readings were kept below 0.40 by dilution with the appropriate media. Culture supernatants from bioreactor samples were analyzed for glucose, acetate, butyrate, acetoin, ethanol, acetone, and butanol by using a Waters (Milford, Mass.) high-performance liquid chromatography system (5, 33).
Growth conditions and maintenance. Strains were stored at -85°C in clostridial growth medium (CGM) (35) with 15% glycerol and were grown at 37°C in an anaerobic chamber (Forma Scientific, Marietta, Ohio) in tubes containing 10 ml of CGM with 80 g of glucose/liter or on agar-solidified 2xYTG plates supplemented with erythromycin (100 µg/ml in liquid; 40 µg/ml on plates; Sigma, St. Louis, Mo.). Individual colonies were transferred to 10-ml CGM tubes and heat shocked at 70 to 80°C for 10 min to induce sporulation.
Bioreactor experiments. Preinoculum flasks were grown to an A600 of 0.2, and the cultures were transferred to either a BioStat M (working volume, 1.6 liters; B. Braun, Allentown, Pa.) or a BioFlo II (working volume, 3.6 liters; New Brunswick Scientific, Edison, N.J.) reactor in a 10-fold dilution of CGM supplemented with 75 µg of clarithromycin (Abbott Laboratories, Abbott Park, Ill.)/ml. Fermentations were pH controlled at a low end of 5.0 with 6 N NH4OH. Anaerobic conditions were maintained by using nitrogen at flow rates of 55 and 125 ml/min, respectively, for an A600 below 1.0 and 11 and 25 ml/min for an A600 above 1.0.
Butanol challenge experiments. Static flask cultures (800 ml) supplemented with 75 µg of clarithromycin/ml were grown to an A600 of 0.6 to 0.7 and then were divided equally (working volume, 200 ml) into four 250-ml bottles. Three bottles were challenged with butanol (0.2%, 0.4%, and 0.6% [vol/vol]; Sigma). RNA samples were taken from challenged and unchallenged cultures at 0.17, 0.5, 1, 2, 6, 12, and 24 h poststress, where the indicated time is when samples were first placed in the centrifuge to be spun down.
RNA purification and cDNA labeling. Isolation of RNA and reverse transcription labeling of cDNA probes were performed as previously described (30).
cDNA arrays and hybridizations. The targeted 1,019-gene microarrays used in these experiments and the hybridization protocol have been described previously (30). To minimize labeling biases, dyes were swapped when hybridizations were repeated for the same time point. Samples from the 824(pMSPOA) bioreactor experiments were hybridized against oppositely labeled (Cy3 or Cy5) 824(pIMP1) samples based on similarity in growth stages. Samples from butanol-challenged cultures were hybridized against oppositely labeled unchallenged samples from the same strain.
Microarray analysis. Microarray data were normalized by using a segmental nearest-neighbor approach that has been shown to be superior to other normalization methods (36). The microarray analysis methodology has been published previously (30) and was used here with the following modification. The constant ß used in the filtering criterion was set to 1.282 so that the channel intensity with background and nonspecific binding subtracted was greater than the noise of the background to a 90% confidence interval. All previous results were renormalized by using the same criterion, and it was verified that this reanalysis had little to no effect on the general expression patterns of the previously published data (data not shown).
Differentially expressed genes were grouped on the basis of similarity in expression profile by using self-organizing maps (SOMs) (GeneCluster, version 2.0 [28]) or by average linkage hierarchical clustering (Cluster [10]) and were visualized using TreeView (10).
Metabolic flux analysis. Metabolic flux analysis uses measurable substrate and product data to solve a system of linear equations based on metabolic reaction stoichiometry to determine in vivo fluxes (22). We used a program specifically designed for C. acetobutylicum (9) to examine pathway fluxes (expressed in millimolar concentrations per hour per gram of cell weight) obtained in this study.
|
|
|---|
![]() View larger version (19K): [in a new window] |
FIG. 1. Growth profiles for 824(pMSPOA) (solid squares) and 824(pIMP1) (open squares). Capital letters indicate the times at which microarray analyses were performed (this study), and the arrows show the times at which Northern blot analyses were performed previously (13).
|
![]() View larger version (26K): [in a new window] |
FIG. 2. Comparison of ratios between 824(pMSPOA) and 824(pIMP1) spo0A (a) and aad-ctfAB (b) transcripts from microarrays (open circles) (this study) and Northern blot analysis (solid diamonds) (13). Asterisks indicate differential expression with a confidence level of at least 95% at the time point for the microarray data. The aad-ctfAB Northern blots are shown below the plots. Lanes indicate time points: 1 and 2, exponential growth, 3, transitional growth; 4 and 5, early-stationary phase. Total radioactive counts have been reported previously (13), and time points 1, 2, and 5 were treated as having equal (background) counts.
|
E in B. subtilis (26), and it has been shown that in B. subtilis, expression of the acetyl coenzyme A (acetyl-CoA) carboxylase subunits accD and accA is positively regulated by Spo0A (20). Thus, these transcriptional observations may be attributed to the documented accelerated sporulation of 824(pMSPOA) (13). The transcriptional data of Fig. 4A also suggest that these genes may be organized into two operons (one ending at CAC3577 and the other starting at CAC3576) or that there is a secondary regulation of mRNA abundance for some genes within an operon.
![]() View larger version (84K): [in a new window] |
FIG. 3. Expression profiles for all differentially expressed (95% confidence interval) genes; the ratio is 824(pMSPOA)/824(pIMP1). Plots indicate the time course of the average normalized expression ratios within the cluster.
|
![]() View larger version (48K): [in a new window] |
FIG. 4. (A) Structure and comparison of expression profile in a predicted fatty acid synthesis operon (CAC3568 to CAC3579). The gene identified as CAC3580 (annotated as related to 2-nitropropane dioxygenase) was omitted from analysis for missing too many data points. The fatty acid synthesis pathway is based on the work of de Mendoza et al. (8) and the Kyoto Encyclopedia of Genes and Genomes (KEGG; www.genome.ad.jp/kegg/kegg2.html). P, predicted promoter; ACP, acyl carrier protein; R-ScoA, straight- or branched-chain acyl-CoA primer. (B) Structure and comparison of primary metabolism expression profiles in C. acetobutylicum spo0A strains as generated by microarrays. The last column indicates whether differential expression was achieved in each experiment to a 95% confidence interval on at least one array for at least two time points; an asterisk before the slash indicates differential expression in the wild-type (WT)/SKO1 experiment, and an asterisk after the slash indicates differential expression in the 824(pMSPOA)/824(pIMP1) experiment. The colorimetric ratio is on the same scale as in Fig. 3, and the yellow triangle represents the direction of increasing culture time.
|
We also examined if lower expression of glycolytic-pathway genes in 824(pMSPOA) resulted in a corresponding decrease in metabolic flux through glycolysis. Metabolic flux analysis was performed using fermentation data (generated from the same experiment), and the ratio of fluxes [824(pMSPOA)/824(pIMP1)] was calculated (data not shown). At times corresponding to stationary-phase growth, rGLY1, the flux corresponding to the glucose
pyruvate pathway, and rGLY2, the pyruvate
acetyl-CoA pathway flux, were lower in 824(pMSPOA) than in 824(pIMP1). At microarray time point F (t = 22 h) rGLY1 was 38% lower in 824(pMSPOA), and at time point G (t = 25 h), rGLY2 was 42% lower in 824(pMSPOA). Therefore, the downregulation of several glycolytic-pathway genes in 824(pMSPOA) corresponds with lower metabolic fluxes.
A number of heat shock-related genes also had a pattern of higher expression at early times in 824(pMSPOA); these genes included groES (CAC2704), htpG (hsp90 [CAC3315]), and hsp18 (CAC3714) and were included in cluster P3. Also included in this cluster were three genes that are apparently in the cysteine biosynthesis pathway: an ABC sulfate transporter (CAC0107; a cysP homolog), cysD (CAC0109), and cysK (CAC2235). Cluster P0 included a DnaJ-family molecular chaperone (CAC0648), the solventogenic gene adc (CAP0165), a second fabG transcript (CAC3462), an alcohol dehydrogenase gene (CAC1331), and the sigma factor sigG (CAC1696). Cluster P1 included two endoglucanases (CAC0214 and CAC2556).
824(pMSPOA) exhibits increased tolerance and prolonged and enhanced solvent production under butanol stress.
Static flask cultures (n = 3) of 824(pMSPOA) and 824(pIMP1) were grown to exponential phase (A600,
0.6 to 0.7) and then split into four separate flasks, three of which were challenged with 0.2, 0.4, or 0.6% (vol/vol) butanol. The challenged cultures were not typically inhibited within the first 5 h postchallenge (Fig. 5). Butanol-challenged 824(pIMP1) cultures showed 60 to 80% of the unchallenged culture density in stationary phase (t = 25 to 70 h postchallenge). In contrast to the results of bioreactor fermentation experiments (Fig. 1), 824(pMSPOA) grew to significantly lower cell densities than 824(pIMP1) in static flask cultures (n = 3) (Fig. 5). In contrast again to bioreactor experiments (Fig. 1), 824(pMSPOA) produced significantly larger amounts of solvents (despite the lower cell densities) than 824(pIMP1) in static flask cultures (Fig. 6). We have no explanation for these differences, which are presumably due to pH control-associated effects. 824(pMSPOA) cultures challenged with 0.2 and 0.6% butanol typically showed densities 60 to 85% of the unchallenged culture density. In three independent butanol challenge experiments with strain 824(pMSPOA), there was significantly more variability in the cellular response to various butanol levels (0.2, 0.4, and 0.6%) than with other strains [the wild type, 824(pIMP1), and the two strains in reference 31] (data not shown). In some experiments the challenged culture grew better than the unchallenged culture. We believe that this is due to the highly heterogeneous population of the faster-differentiating cells and the very viscous and clumping nature of the 824(pMSPOA) cells (see also reference 13). This apparently results in somewhat different cell populations in the smaller flasks when the larger mother flask is split into four flasks for the butanol challenge experiments (see Materials and Methods). Nevertheless, in all experiments, 824(pMSPOA) cultures showed significantly higher tolerance to butanol challenges than plasmid control cultures or the wild type (data not shown). Unchallenged and challenged 824(pIMP1) cultures stopped utilizing glucose 30 h postchallenge, while 824(pMSPOA) continued to consume glucose over 100 h after the butanol challenge (Fig. 5). 824(pIMP1) cultures challenged with 0.6% butanol could utilize glucose at only 13% of the rate of the unchallenged culture at 9 h postchallenge, and with a 0.2% butanol challenge the glucose utilization rate was 27% that of the unchallenged culture at 18 h postchallenge. 824(pMSPOA) cultures were less inhibited at both the 0.2 and 0.6% challenge levels, with glucose utilization rates at 45 to 75% that of the unchallenged culture up to 70 h postchallenge.
![]() View larger version (48K): [in a new window] |
FIG. 5. Growth and glucose utilization and inhibition profiles for 824(pMSPOA) (solid symbols) and 824(pIMP1) (open symbols) cultures challenged with butanol. Circles, unchallenged cultures; squares and diamonds, cultures challenged with 0.2 and 0.6% butanol, respectively. Arrows on the x axis indicate time of butanol challenge.
|
![]() View larger version (37K): [in a new window] |
FIG. 6. Product profiles for butanol-challenged cultures of 824(pMSPOA) (solid symbols) and 824(pIMP1) (open symbols). Circles, unchallenged cultures; squares and diamonds, cultures challenged with 0.2 and 0.6% butanol, respectively. Arrows on the x axis indicate time of butanol challenge.
|
Transcriptional analysis of butanol-stressed 824(pMSPOA) cultures. RNA samples of the challenged (0.6% [vol/vol] butanol) and unchallenged cultures were taken at 0.17, 0.5, 1, 2, 6, 12, and 24 h after butanol challenge. Two microarray slides (with dye swap) were hybridized for each time point with the exception of the 2- and 6-h samples from 824(pIMP1), where three hybridizations were performed. Data were normalized, and differentially expressed genes were identified to a 95% confidence interval (36). For further analysis, 168 genes that had differential expression at two or more time points in the two experiments were identified (Fig. 7A); these genes were organized using average linkage hierarchical clustering (10).
![]() View larger version (81K): [in a new window] |
FIG. 7. (A) Expression ratio profiles for all 168 differentially expressed genes in butanol-challenged 824(pMSPOA) and 824(pIMP1) cultures. (B) Detail from panel A (location marked by orange vertical bar) identifying genes that show significant upregulation in 824(pMSPOA) (designations printed in red). Samples were taken at 0.17, 0.5, 1, 2, 6, 12, and 24 h post-butanol stress. Yellow triangle represents the direction of increasing time following butanol stress.
|
Additional genes upregulated at 0.17 and/or 2 h (Fig. 7B) included pykA (CAC2660, encoding pyruvate carboxylase), an ABC transporter (CAC3101), a thioredoxin reductase gene (CAC3082), a pyruvate-formate lyase gene (CAC0980), pta (CAC1742, encoding phosphate acetyltransferase), and the gyrase holoenzyme subunit genes gyrB (CAC0006) and gyrA (CAC0007). Several fatty acid synthesis genes were upregulated at 2 h, including fabF and fabG [CAC3573, encoding 3-oxoacyl-(acyl carrier protein) synthase I, and CAC3574, encoding 3-ketoacyl-acyl carrier protein reductase] and accB (CAC3572, encoding the biotin carboxyl carrier protein of acetyl-CoA carboxylase). Genes related to nucleotide and amino acid synthesis in this upregulated cluster include hisC (an amino acid metabolism gene), serA (CAC0015, encoding D-3-phosphoglycerate dehydrogenase, involved in amino acid metabolism), radC (CAC1241, a DNA repair gene), uracil-DNA glycosylase (CAC2997, a DNA repair gene), purQ/purL (CAC1655, a purine metabolism gene), sigA (CAC1300, encoding primary sigma factor
A), spo0J (CAC3729, a DNA packing and segregation gene), and truB (CAC1805, encoding pseudouridine synthase, involved in pyrimidine metabolism). We conclude that early after butanol challenge, strain 824(pMSPOA) showed higher expression of several genes related to butyryl-CoA and butyrate synthesis and assimilation; DNA synthesis, repair, and topology; and fatty acid synthesis. Some of these genes may be related to the enhanced butanol tolerance of this strain.
Comparison of spo0A- and groESL-induced butanol tolerance. We also wanted to identify which genes are part of the general response to butanol stress and which genes might be uniquely associated with the stress response of strain 824(pMSPOA). To do so, we compared our results to data from butanol challenge experiments involving the heat shock protein-overexpressing strain 824(pGROE1) (which overexpresses groESL with a thiolase promoter [31, 32]) at two different butanol challenge levels (0.25 and 0.75% [vol/vol]) and its plasmid control 824(pSOS95del) (0.75% butanol challenge only). Strain 824(pGROE1) also exhibits increased butanol tolerance (31).
We identified 160 genes that were differentially expressed in at least three of the experiments. These were grouped by average linkage hierarchical clustering (Fig. 8A). There is clearly a region of common upregulation in all five experiments (Fig. 8B). Several stress response genes are in this group, including clpC (CAC3189), dnaJK (CAC1283 and CAC1282), groESL (CAC2704 and CAC2703), hsp18 (CAC3714), and hrcA (CAC1280, encoding a regulator of class I heat shock protein transcription). The solvent formation genes aad, ctfA, and ctfB had higher expression at nearly every time point across all five experiments. This is consistent with the observations of Tomas et al. (31) and further suggests the classification of these genes as part of the global stress response. Other strongly upregulated genes include citBC (CAC0971, encoding aconitase A, and CAC0972, encoding isocitrate dehydrogenase; earlier annotations of the C. acetobutylicum genome list these genes as leuCD and leuB) and a pSOL1 megaplasmid membrane protein (CAP0102, homologous to B. subtilis ydgH). While C. acetobutylicum does not appear to utilize the full Krebs cycle in its metabolism, citBC have been found to be necessary for sporulation in B. subtilis. Deletion of citC in B. subtilis blocks sporulation in stage I (16), and Escherichia coli citC mutants have been found to be thermosensitive (6). B. subtilis citB mutants do not sporulate efficiently and are typically blocked in stage 0; the mutation is thought to interfere with the spo0A phosphorelay (7).
![]() ![]() ![]() ![]() View larger version (350K): [in a new window] |
FIG. 8. (A) Comparison of differential expression of genes due to butanol challenge in five experiments. Colored vertical bars indicate regions of interest shown in detail in panels with corresponding bars. (B through D) Detailed views of genes that were universally upregulated (B) (red), upregulated at 0.17 h postchallenge in 824(pMSPOA) but not in 824(pGROE1) (C) (blue), and upregulated at 0.17 and/or 2 h in 824(pMSPOA) but not at all in 824(pGROE1) (D) (green). 824(pMSPOA) and 824(pIMP1) samples were taken at 0.17, 0.5, 1, 2, 6, 12, and 24 h after butanol stress. 824(pGROE1) and 824(pSOS95del) samples were taken at 0.25, 1, 3, 6, 12, and 24 h after butanol stress. Yellow triangle represents the direction of increasing time following butanol stress.
|
|
|
|---|
Due to the highly complex nature of transcriptional regulation, it is expected that either many of the genes differentially expressed due to spo0A overexpression (Fig. 3) would be directly regulated by Spo0A binding or their differential expression would be the product of a secondary effect due to morphological changes or sigma factor binding (11, 20). Previous predictions (30) of 0A binding sites in intergenic regions were compared with the differentially expressed genes diagramed in Fig. 3 in order to search for enrichments within clusters in the number of genes predicted to be Spo0A regulated. By using the 0A box sequence 5'-TGNCGAA-3' and previously described predicted transcriptional units (30), 51% of all genes were predicted to be regulated by Spo0A; the frequencies of predicted Spo0A regulation in the individual clusters of Fig. 3 ranged from 36 to 56%. Because the 0A-binding motif is only 7 bases long, the high number of false positives masks the actual number of genes regulated by Spo0A.
The arrays used in these experiments are targeted microarrays that contain approximately one-fourth of the C. acetobutylicum genome. Genes in this generation of arrays include, among others, nearly all pSOL1 open reading frames, 123 DNA replication and repair genes (90% of the total of such genes as identified by the genome annotation), 97 cell division- and sporulation-related genes (92%), 85 carbohydrate and primary metabolism genes (31%), 67 energy production genes (52%), 63 outer membrane and cell envelope genes (36%), 48 lipid metabolism genes (80%), 42 motility and chemotaxis genes (39%), and all previously identified stress response genes.
The sporulation sigma factor sigF (CAC2306) exhibited nominally higher expression at time points A to F in 824(pMSPOA) [maximum, 23% higher expression than in 824(pIMP1) (data not shown)], but none of these time points were statistically significant. At time point G, the gene was differentially downregulated [higher expression in 824(pIMP1)]. This suggests an accelerated time course similar to the expression of aad/ctfAB (Fig. 2). In B. subtilis, the spoIIA operon (spoIIA-spoIIAB-sigF) is transcribed with the
H sigma factor, a process positively regulated by Spo0A through the presence of four 0A boxes (3). The current prediction (30) of sigma factor and 0A box binding sites of the spoIIA operon in C. acetobutylicum indicates potential
E,
F,
G, and
H binding sites and no 0A boxes. These results suggest that transcriptional regulation of sigF in C. acetobutylicum may not be Spo0A dependent.
In B. subtilis the transitional-phase gene regulator abrB is known to be transcribed with the sigma factor
A; transcription is negatively regulated by Spo0A by binding to one of two 0A boxes (27). The C. acetobutylicum genome (21) contains three genes (CAC0310, CAC1941, and CAC3647) annotated as being abrB-like, and two (CAC1941 and CAC3647) of these three genes were included on the microarrays used in this study. Neither gene was differentially expressed in 824(pMSPOA). Due to the high degree of homology and short lengths (240 and 246 bp) of these genes, the microarray probes for the two abrB genes contained a region with 127 out of 161 bp matching; there was also a comparable amount of homology between these probes and an equally short CAC0310 (246 bp). With this degree of similarity, there was probably a high degree of nonspecific binding. The predicted sigma factor binding sites (30) are different for the genes. CAC0310 has predicted
A,
E,
F,
G, and
H binding sites and a predicted 0A box; CAC3647 has predicted
A and
E binding sites and a predicted 0A box. CAC1941 is predicted to lie on a transcript with an adjacent gene (CAC1940) (the B. subtilis abrB gene occurs as a singleton) with a
G binding site and no 0A box. With the information currently available, it is not clear what the roles of the three abrB genes might be in C. acetobutylicum.
824(pMSPOA) exhibits enhanced butanol tolerance, as demonstrated by prolonged glucose utilization (Fig. 5) and extended solvent production (Fig. 6). It is generally accepted that the main mechanism of solvent toxicity in microorganisms is the chaotropic effect of solvents on the cell membrane, which may inhibit active transport, membrane bound ATPases, and glucose uptake (4), partially or completely abolish the membrane
pH (4, 15, 29) and 
(29), and lower the intracellular pH (4, 15, 29) and ATP concentration (4). It has been shown that clostridial cells increase the saturated-chain content (1, 18, 34) and decrease the fluidity (2) of their lipid membranes in order to adapt to butanol stress. The membrane-related hypothesis of solvent tolerance would suggest that de novo fatty acid synthesis would be required for resistance to solvent stress. Yet the butanol stress experiments indicate that fatty acid metabolism generally decreases upon butanol challenge except in 824(pMSPOA) at 2 h after butanol stress (Fig. 8D).
The observed patterns of differential gene expression in 824(pMSPOA) are likely due to the accelerated differentiation and sporulation of 824(pMSPOA). De novo production of fatty acids has been shown to be required for sporulation in B. subtilis (26), and the accelerated differentiated state in 824(pMSPOA) results in earlier expression of fatty acid synthesis genes than in the plasmid control (Fig. 3). It would be expected that the more robust structures associated with spore development would generally be more resistant to solvents, and this general concept may explain the improved butanol tolerance of 824(pMSPOA). ftsX is found in most bacterial genomes (including B. subtilis and E. coli) and plays a role in cell division and septum formation. Transcription of sigE in B. subtilis (spoIIG) begins near the onset of sporulation, and its product,
E, regulates the transcription of genes involved in spore morphogenesis, genes involved in metabolism within the mother cell, and genes encoding several transcriptional regulators (14). Transcription of such genes apparently leads to enhanced butanol tolerance. Higher expression at early time points of genes related to energy production through butyrate and acetate formation and of genes related to amino acid synthesis, DNA synthesis, and repair may be responsible for, or related to, the improved butanol tolerance of strain 824(pMSPOA).
We thank Christopher Tomas, Seshu Tummala, Stefan Junne, and Carlos Parades for technical assistance and helpful discussions, acknowledge use of the Keck Biophysics Facility and the Center for Genetic Medicine facilities at Northwestern University, and thank Abbott Laboratories for donation of clarithromycin.
|
|
|---|
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»