Previous Article | Next Article ![]()
Applied and Environmental Microbiology, August 2007, p. 5183-5189, Vol. 73, No. 16
0099-2240/07/$08.00+0 doi:10.1128/AEM.02776-06
Copyright © 2007, American Society for Microbiology. All Rights Reserved.

Laboratoire de Microbiologie, Département Environnement Microbiologie et Phycotoxines, IFREMER, 29280 Plouzané, France
Received 28 November 2006/ Accepted 30 May 2007
|
|
|---|
|
|
|---|
To determine if cells induce the expression of virulence genes after the recovery of culturability, a relative quantification of expression rates was performed with standardization by using several genes of reference. The literature shows that, in many cases, using only one gene of reference is not sufficient because the expression can vary considerably depending on experimental conditions (4, 34). The limited concentration of iron in surface seawater suggests that iron is one of the factors that affect the growth of bacteria in the oceans (32). To survive under conditions of limited iron, several genes of the iron uptake pathway are involved, including the pvsA gene required for the biosynthesis and transport of the siderophore vibrioferrin in V. parahaemolyticus (31) and the pvuA gene encoding the ferric vibrioferrin receptor (12). The gene fur, encoding a factor of transcription implicated in the siderophore-mediated iron acquisition system, was also chosen for this study (36). These genes were tested as potential control genes under our experimental conditions.
The objective was to determine (i) if V. parahaemolyticus, after the recovery of culturability under environmental conditions, could be a danger to public health by inducing the expression of virulence genes, such as tdh2 encoding TDH and VPA1354 (escU), VPA1342 (spa24), and VPA1346 (vopP) encoding cytosolic, inner membrane, and effector proteins of TTSS2, and (ii) if real-time reverse transcription (RT)-PCR could be used in a field survey to determine the potential virulence of this pathogenic bacterium. A real-time RT-PCR method was developed to quantify the relative expression ratios of the virulence genes and to detect the viable states of V. parahaemolyticus bacteria, culturable or not, by commonly used techniques based on mRNA amplification. The validity of this approach for the detection of viability was tested by calculating the half-lives of several mRNA species at a low temperature (4°C).
|
|
|---|
VBNC cells in artificial seawater (ASW) microcosms at 4°C were obtained by a method described previously (7). To measure the recovery of culturability, 50-ml volumes harvested from ASW microcosms 1 day after the VBNC state was reached were incubated at 20 or 37°C with shaking (70 rpm) in the dark. The CFU counts on HI (2% NaCl) agar plates after incubation at 37°C for 24 h were determined.
Strains of V. parahaemolyticus (tdh negative, tdh negative, tdh+, and/or trh+), non-O1/non-O139 V. cholerae, V. mimicus, V. alginolyticus, V. vulnificus, Aeromonas hydrophila/caviae, A. salmonicida, Pseudomonas aeruginosa, Shewanella putrefaciens, and Escherichia coli O157:H7 were selected to test the specificity of the primers used in this study.
Nucleic acid extractions.
DNA and RNA extractions were performed as previously described (7).
Probes and primers.
PCR and real-time RT-PCR were performed with the primers described in Table 1. In this study, probes and primers specific for virulence and housekeeping genes were designed based on gene sequences from 5 and 15 strains of V. parahaemolyticus, respectively. The design of the primers and probes was performed as previously described (22).
|
View this table: [in a new window] |
TABLE 1. Primers used in PCR and real-time RT-PCR and probes used in slot blot hybridization for V. parahaemolyticus
|
Sequencing.
The PCR product obtained with each set of primers for the DNA from strain Vp5 of V. parahaemolyticus was purified using a column system (QIAquick PCR purification kit) according to the recommendations of the manufacturer (QIAGEN). The nucleotide sequence of the purified PCR product was determined as previously described (9).
Slot blot hybridization.
The slot blot hybridization approach was used to determine the specificity of the PCR amplicons obtained with newly designed primers. Slot blots were prepared using positively charged membrane (Roche Diagnostics) and processed by the method of Sambrook et al. (28). Briefly, the PCR-amplified products were denatured (100°C for 10 min) and transferred onto nylon membrane by using the slot blot apparatus BIO-DOT SF (Bio-Rad, Marnes la Coquette, France). The DNA was cross-linked to the membrane with UV light for 3 min. Oligonucleotide probes specifically designed for this assay (Table 1) were labeled using the Dig oligonucleotide tailing kit (Roche Diagnostics) and hybridized at the annealing temperature in a solution containing 5x SSC (0.75 M NaCl, 0.075 M sodium citrate [pH 7.0]), 1.0% (wt/vol) blocking reagent (Roche Diagnostics), 0.1% (wt/vol) N-lauroylsarcosine (Sigma), and 0.02% (wt/vol) sodium dodecyl sulfate (SDS; Sigma). Hybridized blots were washed twice in 2x SSC-0.1% SDS for 5 min and three times in 0.5x SSC-0.1% SDS for 5 min at the annealing temperature. The membrane was incubated with buffer containing the antidigoxigenin-alkaline phosphatase conjugate for 30 min at room temperature. The immunoblots were washed with buffer (0.1 M maleic acid, 0.15 M NaCl [pH 7.5], 0.3% Tween 20) for 15 min and developed using p-nitroblue tetrazolium chloride and BCIP (5-bromo-4-chloro-3-indolylphosphate) toluidine substrates (Roche Diagnostics) at room temperature.
cDNA synthesis.
RT was performed with the SuperScript III first-strand synthesis system for RT-PCR as described by the manufacturer (Invitrogen, Cergy-Pontoise, France). The first cDNA strand was obtained using random hexamers with 800 ng of DNase I-treated RNA. To determine the potential inhibition during the RT and/or PCR, 1 ng of MS2 RNA (Roche Diagnostics, Meylan, France) was added to each RT tube. The primers Oco-1 and Oco-2 (Table 1) were used in real-time RT-PCR to detect MS2 RNA. Each population of cDNA was also used for the construction of a standard curve specific to the corresponding targeted gene.
Real-time RT-PCR conditions.
Real-time RT-PCR was performed with viable culturable and VBNC cells (between day 0 and day 106 of the study) and VBNC samples subjected to temperature upshifts to 20 and 37°C. Detection by real-time RT-PCR was performed with the brilliant SYBR green quantitative PCR core reagent kit (Stratagene, Amsterdam, The Netherlands). Each PCR was performed with a 25-µl final volume containing 2 µl of cDNA with 1x core PCR buffer, 3.5 mM MgCl2, 200 µM (each) deoxynucleoside triphosphates, 300 nM (each) primers, 3% dimethyl sulfoxide, 0.167x SYBR green I, and 1.25 U of SureStart Taq DNA polymerase. The reactions were run on an Mx3000P quantitative PCR system (Stratagene). The following thermal cycling conditions were used: a denaturation program (95°C for 10 min), an amplification program repeated 40 times (95°C for 15 s and 60°C for 1 min), and a melting-curve program (55 to 95°C with warming of 0.2°C per s). RT and PCR positive controls (RNA and DNA, respectively) and negative controls (deionized water) were included in each run.
Construction of specific standard curves.
Specific PCR products were identified based on predicted amplicon sizes and extracted using the QIAquick gel extraction kit (QIAGEN, Courtaboeuf, France). Extracted PCR products were cloned using a TOPO TA cloning kit (pCR2.1-TOPO vector and One Shot TOP10F' competent E. coli cells; Invitrogen). After the selection of bacterial colonies containing plasmids with the correctly sized inserts, plasmids were extracted and purified using the HiSpeed plasmid midi kit (QIAGEN). The presence of the targeted sequence was verified by PCR with specific primers. In a total volume of 40 µl, 2 µg of plasmid DNA was digested by 4 U of NotI enzyme (Roche Diagnostics) for 3 h at 37°C. For each gene, plasmid DNA dilution curves were generated and used to calculate the efficiency (E) of the real-time PCR [E = 10(–1/slope)].
Relative quantification and statistical analysis.
The relative expression ratio was calculated for each gene of interest by using a mathematical model described by Pfaffl (26):
![]() |
CT(target) and
CT(ref) represent the CT deviation of control (untreated) minus sample (treated) of the target and reference gene transcripts, respectively. The stability of the transcription of the potential reference genes was evaluated using the freely distributed Microsoft Excel application geNorm (34). The genetic expression in culturable V. parahaemolyticus Vp5 cells, those maintained in ASW at 4°C (day 0) and those subjected to temperature upshifts to 20 and 37°C (responsible for the recovery of the culturability of VBNC cells), was quantified. Samples were taken at two different times, 2 and 6 days after the temperature upshifts. The real-time RT-PCR analysis was performed twice in triplicate. The most stable reference genes were identified, and a geometric mean was calculated and used as a normalization factor. The latter allowed the normalization of the relative expression ratio of each gene of interest. The expression ratios of the genes of interest were tested for significance by a permutation test implemented in the relative expression software tool REST which allows the user to determine if observed effects are the result of the treatment or random allocation (27).
Procedures for measurement of RNA decay.
V. parahaemolyticus Vp5 was inoculated as previously described into a final 2-liter volume of ASW in a 5-liter flask and stored at 4°C. Rifampin (Sigma Aldrich Chimie, St. Quentin Fallavier, France) was added to the ASW microcosm at a final concentration of 100 µg/ml to block transcription. After rifampin treatment, 50-ml volumes were harvested from the ASW microcosm at serial time points (0, 1, 2, 3, 4, 5, 6, and 7 days). The RNA extractions were performed as described above, and the decay of RNA was analyzed by real-time RT-PCR. Additionally, CFU were enumerated after the passage of 10-ml volumes through a 0.22-µm-pore-size nitrocellulose filter (Millipore, St. Quentin en Yvelines, France). The filters were placed on HI agar (2% NaCl) in duplicate and incubated at 37°C for 48 h.
Nucleotide sequence accession numbers.
The GenBank accession numbers for the V. parahaemolyticus Vp5 sequences obtained in this study are EF139123 (rpoS), EF139124 (pvsA), EF139126 (fur), EF139125 (pvuA), EF139129 (tdh2), EF139128 (escU), EF139127 (vopP), and EF139122 (spa24).
|
|
|---|
![]() View larger version (27K): [in a new window] |
FIG. 1. Specificity of primers for V. parahaemolyticus used in real-time RT-PCR. PCR products were hybridized with digoxigenin-labeled probes P2-rpoS (A) and P3-tdh2 (B) specific for rpoS and tdh2 gene amplicons, respectively. The specificity was tested with various strains of V. parahaemolyticus, the tdh1+ tdh2+ (positions 1 to 4), tdh2+ trh+ (positions 5 and 6), tdh1 tdh2 trh negative (position 7), and trh+ (positions 8 to 14) strains, and strains of V. cholerae (position 15), V. mimicus (positions 16 and 17), V. alginolyticus (positions 18 and 19), V. vulnificus (positions 20 and 21), A. hydrophila/caviae (position 22), A. salmonicida (position 23), P. aeruginosa (position 24), S. putrefaciens (position 25), and E. coli O157:H7 (position 26). The no-template control (NTC) contained deionized water.
|
|
View this table: [in a new window] |
TABLE 2. Standard curve evaluation of putative control and virulence genes in V. parahaemolyticus
|
|
View this table: [in a new window] |
TABLE 3. geNorm stability measurements and rankings of expression of selected housekeeping genes
|
![]() View larger version (16K): [in a new window] |
FIG. 2. Mean relative expression ratios for five genes in ASW after temperature upshift. The temperature of ASW was increased from 4°C to 20°C (A and C) or 37°C (B and D). The levels of expression after 2 days (D2) and 6 days (D6) of temperature upshift were computed and were normalized according to the geometric mean of results obtained for two control genes, pvsA and pvuA. The data are the normalized means ± standard deviations of the results for two runs, each with three replicate samples. *, significant results (P < 0.05).
|
![]() View larger version (10K): [in a new window] |
FIG. 3. Transcriptional activity in V. parahaemolyticus Vp5 in ASW at 4°C. Quantification of the mRNA copy numbers for the rpoS, fur, pvsA, pvuA, tdh2, escU, vopP, and spa24 genes was performed by real-time RT-PCR. The number of days necessary to decrease the mRNA population by half (T ) is indicated for each gene. Sample graphs depicting the reduction in the mRNA for a housekeeping gene (rpoS; panel A) and a virulence gene (tdh2; panel B) are shown.
|
![]() View larger version (12K): [in a new window] |
FIG. 4. Half-lives of mRNA species of V. parahaemolyticus Vp5 in ASW at 4°C. The table summarizes the half-lives calculated for mRNA for eight genes, i.e., rpoS, fur, pvsA, pvuA, tdh2, escU, vopP, and spa24. Sample graphs representing the quantities of mRNA for a housekeeping gene (rpoS; panel A) and a virulence gene (tdh2; panel B) are shown. , MS2 RNA was used as an internal control for detecting potential inhibitors of quantitative RT-PCR in RNA samples. , mRNA copies were quantified on days 1, 2, 3, 4, 5, 6, and 7 after a rifampin treatment which blocked the transcriptional activity in the clinical strain V. parahaemolyticus Vp5. The data are expressed as the means ± standard deviations (SD) of results from experiments performed with three replicate samples.
|
|
|
|---|
In vitro, we observed no significant up-regulation of the virulence gene tdh2 and of the genes for TTSS2 in the bacteria after the recovery of culturability. These results are in agreement with those of Baffone et al. (1). Indeed, these authors observed that virulence characteristics disappeared after resuscitation in the mouse and were subsequently reactivated by means of two consecutive passages of V. parahaemolyticus in the rat ileal loop model. The pathogenicity of VBNC cells in vivo has been studied often. Indeed, VBNC V. cholerae cells have previously proven to maintain virulence in human volunteers after ingestion (6), VBNC V. vulnificus cells maintain virulence in a murine model after intraperitoneal injections (23), VBNC V. parahaemolyticus is virulent in a murine model after ingestion (1), and VBNC Campylobacter jejuni is virulent against embryonated eggs (5). However, these studies did not investigate if this virulence was the result of the recovery of culturability and the growth of the cells, with pathogenicity depending also on the infective dose, or if it was caused by an up-regulation of the virulence genes. Previous studies have demonstrated that VBNC bacteria maintain the expression of virulence genes, e.g., a hemolysin gene (vvhA) in V. vulnificus (9, 30, 35), cholera toxin genes (ctxAB) in V. cholerae, and a TDH gene (tdh) in V. parahaemolyticus (35). In our study, we confirmed in vitro that the virulence genes were expressed in VBNC bacteria but not induced in these cells after the recovery of culturability. However, short transitory inductions may be possible, as demonstrated with the yopH and yopE genes in Yersinia spp. after temperature upshifts (10, 17, 19, 20). Nevertheless, contact with eukaryotic cells appears to be necessary for a strong induction of the transcriptional level. This point was made previously for P. aeruginosa, in which the following environmental signals induce the TTSS: (i) the in vitro removal of calcium from the medium, (ii) in vivo contact with eukaryotic host cells, and (iii) the presence of serum (11, 33, 37). The transcriptional induction of the TTSS genes appears also to be a mechanism down-regulated by the transcriptional regulator sigma S, encoded by the rpoS gene, in enterohemorrhagic E. coli and P. aeruginosa (15, 29). Our results confirm that when the rpoS gene is highly expressed, no transcription of the TTSS genes is observed. We also observed a repression of rpoS expression before a decrease in the rate of survival of V. parahaemolyticus cells in ASW at 37°C. Thus, the RpoS factor may be involved in the maintenance of viability of the bacteria in an oligotrophic environment, as described previously for E. coli, Salmonella enterica serovar Typhimurium, and V. cholerae (21, 39).
RT-PCR, which targets short-lived mRNA molecules, has become an increasingly used molecular method for assessing the viability of bacteria, culturable or not. We developed a real-time RT-PCR for the detection of pathogenic strains of V. parahaemolyticus in a VBNC state. In this study, mRNA specific to V. parahaemolyticus was detected in ASW at 4°C and was confirmed to be the result of de novo RNA synthesis. Indeed, the amounts of time taken to decrease the quantity of mRNA by half were four to eight times higher than the mRNA half-lives, indicating transcriptional activity in V. parahaemolyticus in ASW at 4°C. In our study, we found that the rpoS mRNA half-life was 3.5 days. These results are different from previous findings, which indicated that the half-life of rpoS mRNA in V. vulnificus is less than 60 min (30). However, a number of experimental and biological variations can explain this discrepancy. To our knowledge, this is the first time that half-lives have been determined using multiple mRNA populations under conditions which induce the VBNC state. Because of the observed differences in mRNA half-lives, caution is needed when measuring the expression of more than one gene. In previous studies, the average half-life of E. coli mRNA was determined to be 1 min, with a range of 40 s to 20 min and with 80% of all mRNA half-lives computed to be between 3 and 8 min (2, 3). The mRNA turnover appears to have a large impact on transcript levels (16). However, the half-lives were not assessed under conditions which produce the VBNC state, i.e., in ASW at 4°C. Thus, the viability and the gene expression appear to be difficult to confirm because mRNA detected in the VBNC samples may result from residual mRNA and not from de novo RNA synthesis. Thus, one should take this into account when reporting the results obtained in previous studies on the expression of housekeeping (7, 13, 38) and pathogenicity (9, 30, 35) genes in VBNC populations. Indeed, the detection of mRNA in an early population of VBNC cells should be associated with a study of the mRNA half-lives under the same conditions to ensure the viability of the bacterial population.
In conclusion, a real-time RT-PCR method using virulence and housekeeping genes of V. parahaemolyticus was developed and used to study the expression of these genes. After the recovery of culturability, the VBNC cells of V. parahaemolyticus Vp5 did not significantly induce the transcription of the main virulence genes in the environment. mRNAs proved to be good viability markers, but their half-lives, a few days at 4°C, must be considered for the development of experiments and for the interpretation of real-time RT-PCR data. These first results suggest that V. parahaemolyticus needs specific conditions to induce virulence factors. In the future, verification at the translational level could be performed and the possible induction of the virulence genes in the presence of either nutrients or host cells at a favorable temperature could be investigated.
We gratefully acknowledge Jonathan Porter (Environment Agency, Starcross, United Kingdom) and Françoise S. Le Guyader (IFREMER, Nantes, France) for critical reading of the manuscript and Jacques Le Pendu (INSERM, Nantes, France) for his advice on scientific aspects of the project.
Published ahead of print on 8 June 2007. ![]()
|
|
|---|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»