Department of Environmental Toxicology, University of California, Santa Cruz, 269 Jack Baskin Engineering Bldg., Santa Cruz, California 95064,1 Howard Hughes Medical Institute and Department of Molecular Biology and Microbiology, Tufts University School of Medicine, 136 Harrison Avenue, Boston, Massachusetts 021112
Received 10 November 2005/ Accepted 15 February 2006
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
Recent studies have indicated the importance in these processes of an intracellular signaling molecule termed 3',5'-cyclic diguanylic acid (c-di-GMP), which acts as a second messenger for integrating environmental signals that affect cellular physiology (11, 25). These studies have shown the cellular levels of c-di-GMP are controlled by diguanylate cyclases (DGCs) and phosphodiesterases (PDEAs), which are proteins that contain conserved domains termed GGDEF and EAL, respectively (8, 36, 44-48).
Proteins containing DGC and/or PDEA domains are widespread among bacteria (16, 17). Genome analysis of V. cholerae revealed that it has 31 genes that encode DGC proteins, 12 that encode PDEA proteins, and 10 that encode proteins with both the DGC and PDEA domains on the same polypeptide (17). The prevalence of genes encoding DGC- and PDEA-containing proteins in V. cholerae and the evidence that c-di-GMP is an essential modulator of biofilm formation in many microorganisms prompted us to examine the c-di-GMP regulon and associated phenotypic changes in V. cholerae (3, 11, 25, 26, 39, 48).
To analyze the consequences of elevated c-di-GMP concentration in the cytoplasm of V. cholerae, we modulated the c-di-GMP level by overexpressing proteins with DGC or PDEA domains from an arabinose-inducible promoter. We then determined the whole-genome expression profiles of these cells. The study revealed that 4.5% of V. cholerae genes are differentially expressed in response to an increased level of c-di-GMP. The main responses were an increase in the expression of vps genes that are required for Vibrio polysaccharide (VPS) biosynthesis, eps genes involved in the extracellular protein secretion system (EPS), and msh genes required for the mannose-sensitive hemagglutinin (MSHA) type IV pilus biogenesis as well as a decrease in the expression of fla genes required for flagellum biogenesis.
To further understand the response of pathogenic V. cholerae to c-di-GMP, we performed transcriptome analysis in two V. cholerae O1 biotypes: El Tor and classical. The classification of V. cholerae strains is based on the lipopolysaccharide O antigen, which identifies over 200 serogroups in the species. Among these serogroups, only O1 and O139 cause epidemic cholera, as they harbor genes encoding cholera toxin and the colonization factor toxin-coregulated pilus in their genomes. The O1 serogroup can be further divided into the classical and El Tor biotypes on the basis of biochemical properties and phage susceptibility (14).
Although the responses to changes in c-di-GMP level were similar between the classical and El Tor biotypes with respect to the expression of the vps, eps, msh, and fla genes, there were significant differences in the expression patterns of other genes of unknown function. These results demonstrate that V. cholerae responds to changes in the level of c-di-GMP by modulating production of its cell surface properties, which in turn leads to altered motility and biofilm formation. The functions of many potentially important c-di-GMP-regulated genes have yet to be identified.
| MATERIALS AND METHODS |
|---|
|
|
|---|
|
pir strain was used to maintain the pGP704sac28 plasmid carrying the deletion construct pFY-331. Biparental mating was carried out between recipient wild-type V. cholerae and the E. coli S17-1
pir strain harboring the deletion construct. Ampicillin- and rifampin-resistant V. cholerae strains resulting from a single homologous recombination were identified, grown without ampicillin, and then subjected to sucrose selection. Strains with properties of a double recombination event (ampicillin-sensitive and sucrose-resistant) were identified and further analyzed by PCR to confirm that the VCA0956 gene was deleted. In-frame deletion mutants of mshA and the vpsI operon (vpsI) were generated by following the same method described above, using p28-mshA and pAJH5 plasmids (Table 1), respectively.
GFP tagging of V. cholerae strains.
The insertion of the green fluorescence protein (GFP) gene into the V. cholerae genome was performed by triparental mating among V. cholerae, E. coli S17-1
pir carrying the helper plasmid pUX-BF13, and E. coli S17-1
pir carrying pMCM11, which encodes a Tn7 cassette with the GFP gene. Exconjugants were selected on thiosulfate-citrate-bile salts-sucrose agar plates containing gentamicin (50 µg/ml). The mini Tn7 cassette inserts at a specific site in the V. cholerae genome, between the genes encoding VC0487 and VC0488 (nucleotides 519257 to 519261), causing a 5-bp duplication. The confirmation of GFP insertion was performed by colony PCR using primers GFP_1 (5'-CGATGTTGACCAGCCTCGTA-3') and GFP_2 (5'-CGGATCTCGACACAAGCGTA-3'). The growth rate of each GFP-tagged strain was compared to the corresponding parental background to ensure that the GFP insertion did not interfere with growth. Arabinose-inducible vectors were then introduced to GFP-tagged V. cholerae strains via electroporation.
RNA isolation. V. cholerae strains were grown on LB agar plates containing chloramphenicol (5 µg/ml) at 30°C overnight. For the VCA0956 overexpression studies, colonies from semiconfluent plates of each strain were resuspended in 2 ml MOPS medium. The optical density of each culture at 600 nm (OD600) was determined, and cultures were then inoculated to a starting OD600 of 0.05 in MOPS medium. Cultures were grown with aeration at 30°C, and c-di-GMP production was induced by the addition of arabinose to a final concentration of 0.2% when the OD600 of the cultures reached 0.30 to 0.35. Samples were harvested, and RNA was isolated as described previously (55).
Gene expression profiling. Whole-genome expression analysis was performed by using a common reference RNA that was isolated from uninduced V. cholerae AC-1724 cells grown to an OD600 of 0.30 to 0.35. RNA samples from test and reference samples were used in a reverse transcription reaction as follows: 3 µg of RNA and 5 µg of random hexamers were hydrolyzed at 80°C for 8 min and then chilled on ice for 5 min. The reaction mixture, containing first-strand buffer, 0.01 M dithiothreitol, aminoallyl (aa) dNTP labeling mix with aa-dUTP in a 3 (dTTP)-to-2 (aa-dUTP) ratio, and 400 U SuperScript III (Invitrogen), was added to cooled samples and incubated at 42°C for 3 to 4 h. In order to hydrolyze the RNA in the cDNA-RNA mixture, samples were incubated at 65°C for 10 min in the presence of 100 mM NaOH and 10 mM EDTA. After hydrolysis, reactions were neutralized by the addition of 1 M HEPES [pH 7.0], at a final concentration of 500 mM HEPES. Reactions were cleaned up with the QIAquick PCR purification kit (QIAGEN) by using phosphate wash and phosphate elution buffers and dried. The samples were indirectly labeled by covalent coupling of the N-hydroxysuccinimide esters of the cyanine fluorophores to the aminoallyl-labeled cDNAs as follows. Cleaned and dried samples were resuspended in freshly prepared 50 mM Na-bicarbonate buffer, pH 9.0. One vial of Cy3/Cy5 fluor (Amersham CyDye postlabeling reactive dye pack) was mixed with each sample and incubated for 1 h in the dark at room temperature. Unincorporated dye molecules were removed by using the QIAquick PCR purification kit (QIAGEN). Dye incorporation and cDNA concentrations were determined by measuring optical density at 260, 550, and 650 nm. The amount of nucleotides in each sample was determined by using following formula: amount of nucleotides (pmol) = (OD260 · volume [µl] · 37 ng/µl · 1,000 pg/ng)/324.5 pg/pmol. The amount of incorporated dye molecules was determined by using following formula: pmol Cy3 = (OD550 · volume [µl])/0.15; pmol Cy5 = (OD650 · volume [µl])/0.25. The nucleotides-to-dye ratio was then calculated by dividing cDNA picomoles by Cy dye picomoles. Samples with >100 pmol dye incorporation and <10 nucleotides per dye ratio were used in the hybridizations (this method is modified from the Institute for Genomic Research's microbial RNA aminoallyl labeling for microarrays protocol). For hybridizations, dye-coupled test samples were mixed with the corresponding reference sample and dried. Amine-silane slides containing 70-mer oligonucleotides representing most of the open reading frames present in the V. cholerae genome were UV cross-linked at 250 mJoules and prehybridized in a solution containing 5x SSC (1x SSC is 0.15 M NaCl plus 0.015 M sodium citrate), 0.1% sodium dodecyl sulfate, and 1% bovine serum albumin for at least 45 min at 42°C. The hybridization mix was prepared by resuspending dried samples in a solution containing 15 µg salmon sperm DNA, 15 µg tRNA, 3x SSC, and 0.1% sodium dodecyl sulfate. Concentrated hybridization mix was applied to the slides and incubated at 65°C for 12 to 20 h. After hybridization, slides were washed, dried, and scanned by using an Axon scanner to determine the fluorescence in each open reading frame-specific spot. The raw data was obtained by using the software package GenePix 4.1 (Axon). Normalized signal ratios were obtained with LOWESS print-tip normalization by using Bioconductor packages (http://www.bioconductor.org) (18) in an R environment. Differentially regulated genes were determined (with two biological and two technical replicates for each data point) with the significance analysis of microarrays (SAM) package (50) using twofold differences in gene expression and 1% false discovery rate (FDR) as a cutoff value.
Biofilm assays. The biofilm-forming capacities of the V. cholerae strains were determined by crystal violet staining assay after 2, 4, 6, and 12 h of incubating cultures in polyvinyl chloride microtiter plates under static conditions in LB and MOPS medium at 30°C (55, 57). For flow cell experiments, biofilms were grown at room temperature in flow chambers (individual channel dimensions of 1 by 4 by 40 mm) supplied with half-strength LB (0.5% peptone, 0.25% yeast extract, and 1% NaCl) containing chloramphenicol (3 µg/ml) at a flow rate of 13.5 ml h1. Arabinose was added at 0.2% final concentration to induce expression from the ParaBAD promoter. The flow cell system was assembled and prepared as described previously (22). V. cholerae cultures used to inoculate flow channels were prepared by inoculating 6 to 10 single colonies from a plate into flasks containing LB medium and growing them with aeration at 30°C for 20 h. Cultures were diluted to an OD600 of 0.1 in LB and used for inoculation as described previously (55). Image acquisition was done with a Zeiss Axiovert 200 M laser scanning microscope. IMARIS personal software and Adobe Photoshop were used for image processing.
Motility assays. LB and MOPS soft agar plates (0.3% agar) were used to determine the motility of bacterial strains. The diameter of the zone of motility was measured after 20 and 48 h of incubation at 30°C for LB- and MOPS-grown cells, respectively.
Two-dimensional TLC. The amount of c-di-GMP was quantified by using two-dimensional thin-layer chromatography (TLC) as previously described (48).
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
|
2.0-fold change in gene expression and an FDR of
0.01 as criteria. The analysis identified 79 differentially expressed genes (63 genes were induced and 13 were repressed in the E-DGC strain relative to the control pBAD33 strain) after 15 min of arabinose induction. After 30 min, we identified 173 differentially regulated genes (151 induced and 22 repressed). The abundance of the VCA0956 transcript, whose expression was modulated by arabinose to increase cellular c-di-GMP level, increased 59- and 38-fold after 15 and 30 min of induction, respectively. This observation, along with that of the increased cellular c-di-GMP concentration, validates our experimental design. Based on annotated functions provided by the V. cholerae genome sequencing project, the genes regulated by c-di-GMP are predicted to participate in various cellular functions (Fig. 2). Below we discuss a selected subset of the differentially regulated genes and their significance. The complete list of differentially expressed genes is provided as supplemental material.
|
|
Taken together, our results show a positive correlation between the intracellular concentration of c-di-GMP and VPS production in V. cholerae. This finding is congruent with our previous observation that the V. cholerae classical biotype vieA deletion mutant, which accumulates more c-di-GMP than does the wild-type strain, expresses vpsR, vpsA, vpsL, and VC0928 at higher levels (48). In contrast, the same vieA deletion mutation in the V. cholerae El Tor biotype causes neither an increase in c-di-GMP concentration nor an altered vps gene expression (data not shown).
Expression of other factors required for biofilm formation is induced by c-di-GMP. The V. cholerae EPS system (also known as the type II secretion system), is responsible for the secretion of proteins through the outer membrane, including cholera toxin, hemagglutinin/protease, chitinase, neuraminidase, and lipase. Whether eps genes are coordinately regulated with those encoding proteins that are secreted by the EPS system has yet to be determined. Coordinate expression of the vps and eps genes was first observed in the rugose variant of the V. cholerae El Tor biotype, which is capable of producing high levels of VPS (55). In our study, the expression of 3 eps genes located in the 12-gene operon harboring genes epsC-epsN was induced in response to increased c-di-GMP (Fig. 3A) (43). Two products of the EPS pathway, EpsD and EpsE, are predicted to be important for the secretion of either VPS itself or a protein(s) that is involved in the transport or assembly of VPS (1). Our findings suggest that the expression of the eps and vps genes that are required for VPS production and extracellular protein secretion processes may be transcriptionally coordinated by the signaling molecule c-di-GMP.
Expression of the genes mshF, mshG, mshI, mshJ, and mshN, which encode proteins required for the production and function of the MSHA type IV pilus, increased by 1.4- to 1.8-fold in response to elevated c-di-GMP concentrations. MSHA facilitates the attachment of V. cholerae to both biotic and abiotic surfaces, and the biofilm development dynamics of an mshA mutant are markedly different from those of the wild type (7, 33, 51, 52). Since the expression of msh genes increases in response to an increased cellular concentration of c-di-GMP, we show for the first time a correlation between vps and msh gene expression.
An increase in the cytoplasmic c-di-GMP level alters biofilm development dynamics. MSHA and VPS are required for the development of mature biofilms by V. cholerae (27, 51, 52, 54, 57, 59). To determine the physiological consequences of increased vps and msh gene transcription, we compared the biofilm-forming capacities of the E-DGC and E-pBAD33 strains. Initially, we quantified the ability of cells grown both in LB and MOPS at 30°C under static conditions to form biofilms on the surface of microtiter plates. The results presented in Fig. 3B (in LB) and C (in MOPS) reveal that there is a significant increase in the capacity of E-DGC to form biofilms upon induction with arabinose. In LB-grown cultures, at 4 h of biofilm formation, the biofilm-forming capacity of induced E-DGC increased fourfold compared to that of uninduced culture. Similarly, at 6 h and 12 h into biofilm development, E-DGC biofilms were increased by 3.7- and 35-fold, respectively. In MOPS-grown cultures, at 2, 4, 6, and 12 h of biofilm formation, the biofilm-forming capacity of induced E-DGC increased 11.5-, 40-, 43-, and 7.6-fold, respectively.
Transcriptome analysis revealed that the expression of vps genes increased in E-DGC within 15 min of induction with arabinose. However, an increase in the cells' capacity to form biofilms was not observed as rapidly. Biofilms are multicellular structures requiring the synthesis and export of large polymers. The time required for both the production of the exopolysaccharide and synthesis/assembly of the transport system may be one of the reasons for the delay in biofilm formation. Similarly, the differences in physiology and basal transcriptome states of free-living and attached bacteria could result in different phenotypic consequences.
We also attempted to modulate the c-di-GMP level in the strain V. cholerae O1 El Tor by overexpressing the PDEA domain of the VieA protein, which has c-di-GMP phosphodiesterase function, under the control of an arabinose-inducible promoter (47, 48). This strain, referred to as E-PDEA, showed decreased biofilm formation when grown in LB at 2 and 4 h after induction with arabinose (Fig. 3B). However, no significant change in biofilm formation was observed at later time points (6 and 12 h). When grown in MOPS, E-PDEA showed a decrease in biofilm formation after 12 h of induction with arabinose (Fig. 3C). The modest effect on biofilm formation upon induction of VieA expression is consistent with the fact that c-di-GMP levels are already very low (undetectable by two-dimensional TLC [data not shown]) in cells under the conditions used in this experiment. We did not observe a significant increase in biofilm formation upon arabinose induction in the control strain, E-pBAD33 (Fig. 3B and C).
Although these experiments showed that increasing the intracellular c-di-GMP concentration in V. cholerae O1 El Tor alters its biofilm-forming capacity, they did not provide any information about the architecture of the biofilms formed. To investigate the influence of c-di-GMP on biofilm morphology, we first generated a V. cholerae O1 El Tor strain harboring the gene encoding GFP on the large chromosome (E-GFP). The plasmids pBAD33 and pBAD33::VCA0956 were then introduced into E-GFP to obtain strains E-GFP-pBAD33 and E-GFP-DGC, respectively. Each strain formed biofilms in a flow cell system, and images were acquired by confocal scanning laser microscopy over the course of biofilm development. COMSTAT software was used to quantify biofilm features, including average thickness, maximum thickness, and total biomass (22).
A remarkable change in the dynamics of biofilm formation was observed in E-GFP-DGC upon arabinose induction. During the initial stages of biofilm development, there was no significant difference between induced and uninduced E-GFP-DGC cultures (data not shown). However, at 7 and 11 h postinduction with arabinose, E-GFP-DGC formed larger microcolonies than uninduced cultures, presumably because it accumulates a higher level of c-di-GMP. At 25 h, E-GFP-DGC biofilms exhibited full-surface coverage and formed well-developed biofilms, reaching average and maximum thicknesses of 38 µm and 49 µm, respectively. In contrast, the biofilms of uninduced cells were unable to fully cover the substratum and had average and maximum thicknesses of 32 µm and 43 µm, respectively. After 45 h of biofilm development, biofilm detachment characteristics were also markedly different between induced and uninduced cells. Whereas biofilms of the induced cells did not detach, those of uninduced cells showed a marked detachment from the substratum (Fig. 4). In addition, total biomass significantly increased at the 25- and 45-h time points when cells were grown in the presence of 0.2% arabinose. The E-GFP-pBAD33 did not exhibit significant changes in biofilm formation capacity following induction with arabinose (data not shown). These results indicate that an increase in the intracellular c-di-GMP level leads to the development of more compact biofilms, alters biofilm development kinetics, and prevents biofilm detachment.
|
mshA and S
vpsI mutant strains and examined the biofilm development dynamics in a flow cell system. The results shown in Fig. 5 revealed that the surface colonization capacity of the S
mshA mutant and the mature biofilm-forming capacity of the S
vpsI mutant were drastically reduced relative to those of the wild type, and this phenotype did not change in the presence of arabinose. This observation suggests that both MSHA and VPS are required for c-di-GMP-induced biofilms, and their requirement for biofilm formation cannot be compensated for by solely increasing c-di-GMP levels in the cell.
|
|
In addition to genes that encode components of the flagellum, certain genes predicted to encode proteins involved in chemotaxis were differentially regulated in response to increased c-di-GMP levels. The V. cholerae genome contains 43 genes encoding methyl-accepting chemotaxis proteins (MCP). The amount of messages corresponding to two genes that encode MCPs, VCA1034 and VC1405, increased by 2.0- and 1.7-fold, respectively. Conversely, the expression of VC1898, VCA0773, VC2439, VC1298, VCA0176, and VC2161 decreased by 1.4- to 2.9-fold. Thus, it is possible that the decreased motility behavior discussed above could be due to both altered chemotaxis and flagellar motility. We believe that c-di-GMP plays an important role in regulating flagellar motility rather than chemotaxis because the induction of DGC expression caused a reduction in the percentage of motile cells from
80% to
10% (data not shown).
The processes discussed thus far, VPS production and MSHA pili and flagellar motility, are critical to biofilm development (51, 52, 57). Biofilm formation by aquatic microorganisms is well documented (10). It has been proposed that the attachment of bacteria to a surface and subsequent biofilm growth exemplify a survival strategy. Although the steps involved in biofilm formation and some of the necessary structural and regulatory genes involved have been identified, relatively little is known about the signals and genetic networks that control gene expression during biofilm development. Our study highlights the importance of the c-di-GMP second messenger in linking the processes required for V. cholerae biofilm formation, although how the c-di-GMP level changes during biofilm development remains to be determined. Furthermore, in this study, we determined that the expression of VC1888, a protein recently shown to be involved in biofilm development (34), increases in response to elevated c-di-GMP levels. Thus, it is likely that we have identified additional genes that are critical for biofilm development in V. cholerae.
As the overexpression of VCA0956 altered the dynamics of biofilm development, we tested whether a mutation in the gene that encodes VCA0956 (termed cdgF for cyclic-diguanylate protein F) would cause a decrease in the biofilm-forming capacity of the El Tor biotype. The cdgF null mutant did not exhibit any phenotypic differences relative to the wild-type strain with respect to colony morphology, motility, and biofilm formation (data not shown). As the El Tor strain utilized in this study does not have a level of c-di-GMP detectable by two-dimensional TLC (B. Lim, S. Beyhan, J. Meir, and F. H. Yildiz, unpublished data), further decreasing the already low c-di-GMP levels by deleting VCA0956 may not contribute to any phenotypic changes under the conditions tested. Alternatively, since there are 40 other genes in V. cholerae that encode proteins with predicted DGC activity, it is possible that these proteins can compensate for the loss of VCA0956.
The increased concentration of c-di-GMP alters the expression of purine biosynthesis genes. The VCA0956 protein is a DGC that condenses two molecules of GTP to produce c-di-GMP. Hence, we surmised that the overexpression of the DGC-domain of VCA0956 could lead to changes in purine nucleotide biosynthetic pathways. Indeed, this was the case. Following VCA0956 DGC induction in the E-DGC strain expression of VC1190 (purC), which encodes a putative phosphoribosylaminoimidazole-succinocarboxamide (SAICAR) synthetase, decreased by 1.7-fold. SAICAR synthetase catalyzes the ATP-dependent conversion of 5'-phosphoribosyl-5-aminoimidazole-4-carboxylic acid and aspartic acid to SAICAR. In contrast, the expression of VC0767 (guaB), which encodes IMP dehydrogenase, required for the synthesis of GMP from IMP, increased in response to a higher concentration of c-di-GMP. Similarly, the expression of VC1992 (purU), which encodes 10-formyl-THF hydrolase, an enzyme that provides the formate for the synthesis of the purine ring used by 5'-phosphoribosylglycinamide transformylase, also increased.
Purine nucleotides can be synthesized either by de novo or salvage pathways. Synthesis of 5'-phosphoribosyl-1-pyrophosphate (PRPP), which is a precursor of purine and pyrimidine biosynthesis, is catalyzed by PRPP synthase using ATP and ribose-5'-phosphate as substrates. Synthesis of the nicotinamide nucleotides histidine and tryptophan requires 10 to 15% of PRPP for cells growing in minimal medium (58). In our transcriptome analysis, seven of the histidine biosynthesis genes (which are located in an eight-gene operon [VC1132 to VC1139]), were expressed 1.8-fold higher, on average, in response to increased c-di-GMP levels. Four of the genes involved in thiamine biosynthesis (thiC, thiE, thiF, and thiS), which uses an intermediate of the de novo purine biosynthesis pathway (5'-aminoimidazole ribonucleotide) as a reactant, were also expressed at higher levels. These observations may suggest that fluxes in the levels of precursors of the metabolic pathways using either GTP or its intermediates may alter the expression of the genes indicated above.
Increased concentrations of c-di-GMP alter expression of a large number of genes encoding hypothetical proteins. The most striking finding was that 50.9% (88 out of 173) of the genes that were differentially regulated in response to increased c-di-GMP concentration are predicted to encode hypothetical or conserved hypothetical proteins. These genes are organized into 74 different transcriptional units, and 68% of them are located on the small chromosome, while the remaining 32% are on the large chromosome. The identification of a large set of hypothetical proteins regulated by c-di-GMP underlines how little we know of c-di-GMP-regulated physiology and behavior. It is of particular interest that one of these genes, VC0928, which encodes a hypothetical protein whose expression was increased by 4.5-fold after 15 min of VCA0956 DGC induction, is required for the development and stability of biofilm structures in V. cholerae (J. C. N. Fong and F. H. Yildiz, unpublished data).
Analysis of early and late responses to an increased c-di-GMP level. To better understand the immediate and late responses of V. cholerae to c-di-GMP, we compared the transcriptome profiles acquired 15 and 30 min after VCA0956 DGC induction. Differences in transcription profiles between these two time points were determined by using scatter plot analysis. In the graph shown in Fig. 7, the x and y axes denote the changes (n-fold) in gene expression after 30 min and 15 min of VCA0956 DGC induction, respectively, whereby each spot represents a gene that was differentially regulated by twofold in at least one experiment. Genes in E-DGC that were differentially regulated at only 30 min, at only 15 min, and at both time points are indicated. The graph allowed us to divide genes according to their regulation state into four quadrants (I through IV) as follows: quadrant I represents genes whose expression increased at both 15 and 30 min; II represents genes expressed at higher levels at the 15-min time point but lower levels at the 30-min time point; III represents genes expressed at lower levels at both time points; and IV represents genes expressed at higher levels at the 30-min time point but lower levels at the 15-min time point. The genes that were differentially regulated at both time points are found in only quadrants I and III, indicating that their expression is not temporally regulated in response to increased c-di-GMP concentration. Quadrant I contains mostly vps genes (40%) and genes encoding hypothetical proteins (40%), while quadrant III contains the flgC gene and those encoding additional hypothetical proteins. Genes differentially expressed at both time points comprised 16.4% of all the differentially expressed genes, whereas 20.4% and 63.2% were differentially expressed at only 15 or 30 min, respectively. Altogether the results suggest that there are initial and secondary responses to an increase in cellular concentrations of c-di-GMP. The detailed molecular mechanism of the timing of these responses is being investigated.
|
Transcriptome changes due to increased c-di-GMP levels were determined by using the growth conditions described above. SAM analysis revealed genes exhibiting changes in expression of at least twofold and FDRs less than 0.01. Using these criteria, 139 genes were up-regulated and 43 were down-regulated in C-DGC relative to C-pBAD33 after 30 min of induction with arabinose. The complete list of differentially expressed genes is provided as supplemental material. Below we discuss selected cellular processes that are altered by increased c-di-GMP concentration in the classical biotype.
One of the main cellular processes that was differentially regulated in the classical biotype in response to c-di-GMP is VPS production. Twenty-two of the 24 genes found in the "vps region" were up-regulated by 1.7- to 26-fold. The transcript abundance of vpsR and vpsT, the genes encoding two positive regulators of other vps genes, was increased by 7.7- and 9.1-fold, respectively. The magnitude of the increase in vps gene expression, including the positive regulators, vpsT and vpsR, in response to increased c-di-GMP concentration was greater in C-DGC than E-DGC (Fig. 3A). This observation suggests that, in the classical biotype, there is either a higher level of expression of a positive regulator of vpsT and/or vpsR or a lower level of expression of a negative regulator of these genes. We favor the latter explanation, since the classical strain harbors a frame shift mutation in the hapR gene, which encodes a quorum-sensing transcriptional regulator that negatively controls vps gene expression. However, it is also possible that yet-to-be-identified components of the network that regulates vps gene expression contribute to the observed difference.
In C-DGC, as in E-DGC, the expression of genes encoding MSHA (mshD, mshE, mshI, and mshJ) and EPS (epsC, epsJ, and epsN) increased on average by 1.5-fold and 1.4-fold, respectively, further evidence of a correlation between vps, eps, and msh gene expression and increased c-di-GMP concentration (Fig. 3A). Similarly, the expression of a set of flagellar biosynthesis genes (flgF, flgG, flhF, fliH, fliE, and motY) decreased two- to threefold in C-DGC following induction with arabinose. When expression of all the genes involved in chemotaxis and motility was analyzed, we observed a 1.7-fold average decrease in the expression of most of the flagellar biosynthesis genes (flaA, flaB, flaC, flaG, flrB, flrC, flgB, flgC, flgD, flgI, flhA, fliA, fliF, fliG, and fliN) (Fig. 6D). Some of the chemotaxis genes (cheA-2, cheY-3, cheW-2, and cheD) and four MCPs (VC1289, VC2161, VCA0658, and VCA0773) also exhibited decreased gene expression (1.5- to 1.7-fold). In contrast, the expression of two chemotaxis genes (cheV-3 and a putative cheW) and five MCPs (VC0098, VC1859, VCA0864, VCA0988, and VCA1034) increased by 1.4- to 4-fold.
The expression of eight genes encoding DGC-family proteins (VC0703, VC1029, VC1185, VC1216, VC1934, VCA0074, VCA0697, and VCA0785) increased in C-DGC following induction, indicating that the c-di-GMP signal can lead to increased transcription of genes capable of further modulating the c-di-GMP concentration. Among these DGC-family proteins, VC0703, VC1934, and VCA0785 also contain a PDEA domain. In E-DGC, the expression of two genes with both DGC and PDEA domains (VC1934 and VCA0785) increased upon DCG induction. The role of each of the differentially regulated genes encoding DGC- and/or PDEA-family proteins in modulating the intracellular c-di-GMP concentration needs to be determined in order to understand differences in the c-di-GMP response in the two O1 biotypes. At present, only mbaA, which encodes VC0703, has been characterized in detail and found to be a negative regulator of biofilm formation in V. cholerae (4).
The expression of some of the pathogenesis genes (acfB, tagE-1, tcpE, tcpJ, tcpT, and VC1621) located in the Vibrio pathogenicity island increased by 1.6- to 2.7-fold in the classical biotype in response to an increase in the c-di-GMP level. However, the expression of Vibrio pathogenicity island genes and ctx genes was markedly reduced in the classical biotype vieA deletion mutant (data not shown). This mutant produces a high level of c-di-GMP relative to the wild type when grown under a condition that induces virulence gene expression in the classical biotype, in M9 minimal medium supplemented with asparagine, arginine, glutamate, and serine (32). In light of these apparently contradictory findings, the exact mechanism by which c-di-GMP controls the expression of virulence genes in V. cholerae remains to be determined. It should be noted, however, that the observed transcriptional differences between the classical and El Tor biotypes could not be attributed to differences in the nucleotide sequence of VCA0956 between the two strains.
The ability of an organism to modulate the transcriptome in response to extracellular and intracellular signals will greatly influence its overall fitness. In this study, we determined genome-wide transcriptional responses and resulting phenotypic responses of the V. cholerae O1 El Tor biotype to changes in the concentration of the intracellular second messenger c-di-GMP. The results indicated that significant changes to the transcriptome occur in response to an increased c-di-GMP level, which in turn alter cell physiology and lead to alterations in biofilm development and motility behavior. Similar phenotypical changes were observed when the gene adrA (which regulates c-di-GMP concentration) was overexpressed in Salmonella enterica serovar Typhimurium (45), suggesting that different microorganisms respond to cellular increases in c-di-GMP in similar manners. Recent studies in P. aeruginosa indicate that an increased cellular c-di-GMP level can cause large changes in gene expression patterns (23), further suggesting that, in addition to being an allosteric activator (40, 41), c-di-GMP performs a general role as a second messenger (11, 25). The present study also revealed that the classical and El Tor biotypes of V. cholerae exhibit markedly different transcriptional profiles in response to increased c-di-GMP concentration. This finding may be important in light of the reported differences in the regulation of virulence factor production between the two biotypes.
| ACKNOWLEDGMENTS |
|---|
We thank members of the Yildiz laboratory for their suggestions, Lindsay Odell and Vanessa L. Soliven for construction of mshA and vpsI null mutants, and M. Miller and G. Schoolnik for providing plasmids pGP704-sac28, pMCM11, and pUX-BF13. We also thank Lily Shiue for her help in printing FY-Vc array slides.
| FOOTNOTES |
|---|
Supplemental material for this article may be found at http://jb.asm.org/. ![]()
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles: