Assessment of evolution of pandemic Vibrio parahaemolyticus by multilocus sequence typing

The genetic relatedness of 81 isolates of Vibrio parahaemolyticus was assessed by multilocus sequence typing. The strain with serotype O3:K6 emerged as a pandemic pathogen in 1996, with subsequent expansion to include strains having serotypes O1:KUT, O4:K68, and O1:K25. Sequence data from gyrB, recA, dnaE, and gnd revealed that 16 distinct serogroups isolated prior to the pandemic were highly variable and only isolates of serogroup O3:K6 shared two alleles with the pandemic strains. The pandemic strains regardless of serotype were clonal, with 51 of 54 isolates having the identical allelic profile (AP). Serotype alone did not adequately define a pandemic strain: among O1:KUT strains tested, seven strains with the identical pandemic AP carried previously described pandemic markers, while five nonpandemic strains had five distinct APs. Our sequence data provide strong molecular support for the clonal origin of pandemic V. parahaemolyticus O3:K6 and suggest that strains within such a clonal group may acquire previously identified serotypes.     

Global dissemination of Vibrio parahaemolyticus serotype O3:K6 and its serovariants

Vibrio parahaemolyticus is recognized as a cause of food-borne gastroenteritis, particularly in the Far East, where raw seafood consumption is high. An unusual increase in admissions of V. parahaemolyticus cases was observed at the Infectious Diseases Hospital in Calcutta, a city in the northeastern part of India, beginning February 1996. Analysis of the strains revealed that a unique serotype, O3:K6, not previously isolated during the surveillance in Calcutta accounted for 50 to 80% of the infections in the following months. After this report, O3:K6 isolates identical to those isolated in Calcutta were reported from food-borne outbreaks and from sporadic cases in Bangladesh, Chile, France, Japan, Korea, Laos, Mozambique, Peru, Russia, Spain, Taiwan, Thailand, and the United States. Other serotypes, such as O4:K68, O1:K25, and O1:KUT (untypeable), that had molecular characteristics identical to that of the O3:K6 serotype were subsequently documented. These serotypes appeared to have diverged from the O3:K6 serotype by alteration of the O:K antigens and were defined as "serovariants" of the O3:K6 isolate. O3:K6 and its serovariants have now spread into Asia, America, Africa, and Europe. This review traces the genesis, virulence features, molecular characteristics, serotype variants, environmental occurrence, and global spread of this unique clone of V. parahaemolyticus.     

Emergence of a unique O3:K6 clone of Vibrio parahaemolyticus in Calcutta, India, and isolation of strains from the same clonal group from Southeast Asian travelers arriving in Japan.

Active surveillance of Vibrio parahaemolyticus infection among hospitalized patients in Calcutta, India, was initiated in January 1994. The incidence of cases of V. parahaemolyticus infection suddenly increased in February 1996 and has remained high since then. One hundred thirty-four strains of V. parahaemolyticus isolated from January 1994 to August 1996 were examined for serovar, the presence of the thermostable direct hemolysin gene (tdh) and tdh-related hemolysin genes (trh1 and trh2), production of urease, and antibiogram. Strains of the O3:K6 serovar appeared for the first time in February 1996. The O3:K6 serovar strains accounted for 50 to 80% of the strains isolated during the high-incidence period (February to August 1996). All of the serovar O3:K6 strains carried the tdh gene but not the trh genes and did not produce urease. All of the isolates except two were sensitive to all of the antibiotics tested. These and the results of analysis by an arbitrarily primed PCR method indicated that the O3:K6 serovar strains belong to a unique clone. When the O3:K6 serovar strains, isolated from travelers arriving in Japan from Southeast Asian countries, were compared by the arbitrarily primed PCR method, the strains isolated between 1982 and 1993 were distinct from Calcutta O3:K6 while the strains isolated in 1995 and 1996 were indistinguishable from the Calcutta O3:K6 strains. The results suggest that this unique O3:K6 clone may have become prevalent not only in Calcutta but also in Southeast Asian countries very recently. Not only the O3:K6 strains but also the non-O3:K6, tdh-bearing strains isolated in 1996 produced thermostable direct hemolysin at high levels, and thus the level of hemolysin produced does not appear to have influenced the high incidence of serovar O3:K6 strains.     

Characterization of pathogenic Vibrio parahaemolyticus isolates from clinical sources in Spain and comparison with Asian and North American pandemic isolates

In spite of the potential risk involved with contamination of seafood with Vibrio parahaemolyticus, there is a lack of information on the occurrence of pathogenic V. parahaemolyticus in Europe. This organism was isolated in 1999 from a large outbreak (64 cases admitted to a single hospital) associated with raw oyster consumption in Galicia, Spain, one of the most important regions in shellfish production worldwide. Two V. parahaemolyticus isolates from the 1999 Galicia outbreak, three additional clinical isolates obtained in the same period from hospitals in Spain, two reference strains from clinical sources, and five Spanish environmental isolates were examined. Seventeen isolates belonging to the pandemic clone isolated in Asia and North America were included in the study for comparison. All isolates were characterized by serotyping, PCR for virulence-related genes, pulsed-field gel electrophoresis (PFGE), and plasmid analysis. Four of the five clinical isolates from hospitals in Spain belonged to serotype O4:K11; the remaining isolate was O4:K untypeable (KUT). All five isolates were positive for V. parahaemolyticus toxR and tlh (species-specific genes) and tdh and negative for trh and group-specific PCR (a PCR method for detection of the pandemic clone). PFGE analysis with NotI and SfiI discriminated the European isolates in two closely related PFGE types included in a homogeneous cluster, clearly differentiated from the Asian and North American isolates. The five environmental isolates belonged to serotypes O2:K28, O2:KUT, O3:K53, O4:KUT, and O8:K22 and were negative for all virulence genes. The five isolates were discriminated into five different PFGE types unrelated to any other isolate included in the study. While the virulence characteristics (tdh positive, trh negative) of the Spanish clinical isolates matched those of the O3:K6 clone from Asia and North America, they were clearly excluded from this clone by group-specific PCR, PFGE, and serotyping. The results of this study suggest that a unique and specific clone could be related to the V. parahaemolyticus infections in Europe.     

Clonal Diversity among Recently Emerged Strains of Vibrio parahaemolyticus O3:K6 Associated With Pandemic Spread

The genomes of the O3:K6 strains of Vibrio parahaemolyticus which abruptly emerged in Calcutta, India, in February 1996 and which demonstrated an unusual potential to spread and an enhanced propensity to cause infections were examined by different molecular techniques to determine clonality. No restriction fragment length polymorphism (RFLP) in the gene encoding the thermostable direct hemolysin was observed among the O3:K6 isolates of V. parahaemolyticus. Clonal diversity among the O3:K6 strains became evident by examining the RFLPs of the rrn operons and by the use of pulsed-field gel electrophoresis. Five ribotypes were distinguished among the O3:K6 strains examined, with ribotype R4 constituting the major type. Strains of O3:K6 isolated between June and August 1996 showed different pulsotypes compared to the pulsotypes of strains isolated before and after this period, indicating genetic reassortment among these strains, but those isolated between August 1996 and March 1998 showed identical or nearly similar pulsotypes. It is clear that there is a certain degree of genomic reassortment among the O3:K6 clones but that these strains are predominantly one clone.     

Variability of properties of Vibrio parahaemolyticus strains isolated from individual patients

Infections by strains belonging to the O3:K6 pandemic clone of Vibrio parahaemolyticus are prevalent in southern Thailand, and serovariants of these strains have also been detected. V. parahaemolyticus strains lacking important virulence genes (tdh and trh) were isolated from 6.5 to 10.9% of clinical specimens during the period from 2000 to 2003. In order to understand whether changes to the characteristics of V. parahaemolyticus occur during infection, 10 isolates collected from each of 63 patients who presented with diarrhea at the Hat Yai hospital from 2003 to 2004 were examined for the presence of the tdh and trh genes, the O:K serotype, and genetic markers for the pandemic clone. A total of 42 patients (66.7%) yielded identical isolates (homogeneous populations), and 21 of the patients (33.3%) yielded isolates that differed in at least one character from the other isolates (heterogeneous populations). The DNA fingerprints (examined by arbitrarily primed PCR and pulsed-field gel electrophoresis) of some, but not all, of the heterogeneous populations from single patients were indistinguishable. The results indicated that some patients were infected with a unique strain and that in vivo changes (tdh deletion or serotype conversion) might have occurred in certain individuals. It is therefore important to bear in mind that epidemiological studies based on the analysis of a single colony from a single patient might lead to misleading conclusions. Finally, the present study did not rule out the possibility that isolates lacking tdh and trh have unknown virulence mechanisms other than the tdh and trh genes.     

Genotypic analyses of Vibrio parahaemolyticus and development of a pandemic group-specific multiplex PCR assay

A total of 54 Vibrio parahaemolyticus strains including pandemic O3:K6 strains and newly emerged O4:K68, O1:K25, O1:K26, and O1:K untypeable strains (collectively referred to as the "pandemic group") were examined for their pulsed-field gel electrophoresis (PFGE) and arbitrarily primed PCR (AP-PCR) profiles and for the presence or absence of genetic marker DNA sequences, toxRS/new or orf8, that had been reported elsewhere to be specific for the pandemic group. Both PFGE and AP-PCR analyses indicated that all strains of the pandemic group formed a distinct genotypic cluster, suggesting that they originated from the same clone. In addition to the pandemic group, four O3:K6 strains that did not possess the thermostable direct hemolysin (tdh) gene also belonged to this cluster and possessed the toxRS/new sequence. However, three O3:K6 strains that clearly belonged to the pandemic group by PFGE and AP-PCR did not possess the orf8 sequence. The evidence suggests that neither the toxRS/new nor the orf8 sequence is a reliable gene marker for definite identification of the pandemic group. We therefore developed a novel multiplex PCR assay specific for the pandemic group. The assay successfully distinguished pandemic group strains from other V. parahaemolyticus strains by yielding two distinct PCR products for tdh (263 bp) and the toxRS/new sequence (651 bp).     

Isolation of mucoid Vibrio parahaemolyticus strains

Mucoid strains of Vibrio parahaemolyticus were isolated from the stools of two asymptomatic carriers and a patient with gastroenteritis. The strains demonstrated biochemical reactions and antibiotic susceptibility typical of nonmucoid strains of V. parahaemolyticus isolated locally. The slime substance was typed by coagglutination and was antigenically similar to the capsular antigen of the same strain. Three different serotypes (O10:K24, O5:K17, and O5:K15) were involved.     

RS1 element of Vibrio cholerae can propagate horizontally as a filamentous phage exploiting the morphogenesis genes of CTXphi.

In toxigenic Vibrio cholerae, cholera toxin is encoded by the CTX prophage, which consists of a core region carrying ctxAB genes and genes required for CTXPhi morphogenesis, and an RS2 region encoding regulation, replication, and integration functions. Integrated CTXPhi is often flanked by another genetic element known as RS1 which carries all open reading frames (ORFs) found in RS2 and an additional ORF designated rstC. We identified a single-stranded circularized form of the RS1 element, in addition to the CTXPhi genome, in nucleic acids extracted from phage preparations of 32 out of 83 (38.5%) RS1-positive toxigenic V. cholerae strains analyzed. Subsequently, the corresponding double-stranded replicative form (RF) of the RS1 element was isolated from a representative strain and marked with a kanamycin resistance (Km(r)) marker in an intergenic site to construct pRS1-Km. Restriction and PCR analysis of pRS1-Km and sequencing of a 300-bp region confirmed that this RF DNA was the excised RS1 element which formed a novel junction between ig1 and rstC. Introduction of pRS1-Km into a V. cholerae O1 classical biotype strain, O395, led to the production of extracellular Km(r) transducing particles, which carried a single-stranded form of pRS1-Km, thus resembling the genome of a filamentous phage (RS1-KmPhi). Analysis of V. cholerae strains for susceptibility to RS1-KmPhi showed that classical biotype strains were more susceptible to the phage compared to El Tor and O139 strains. Nontoxigenic (CTX(-)) O1 and O139 strains which carried genes encoding the CTXPhi receptor toxin-coregulated pilus (TCP) were also more susceptible (>1,000-fold) to the phage compared to toxigenic El Tor or O139 strains. Like CTXPhi, the RS1Phi genome also integrated into the host chromosomes by using the attRS sequence. However, only transductants of RS1-KmPhi which also harbored the CTXPhi genome produced a detectable level of extracellular RS1-KmPhi. This suggested that the core genes of CTXPhi are also required for the morphogenesis of RS1Phi. The results of this study showed for the first time that RS1 element, which encodes a site-specific recombination system in V. cholerae, can propagate horizontally as a filamentous phage, exploiting the morphogenesis genes of CTXPhi.     

Pandemic spread of an O3:K6 clone of Vibrio parahaemolyticus and emergence of related strains evidenced by arbitrarily primed PCR and toxRS sequence analyses

Vibrio parahaemolyticus O3:K6 strains responsible for the increase in the number of cases of diarrhea in Calcutta, India, beginning in February 1996 and those isolated from Southeast Asian travelers beginning in 1995 were shown to belong to a unique clone characterized by possession of the tdh gene but not the trh gene and by unique arbitrarily primed PCR (AP-PCR) profiles (J. Okuda, M. Ishibashi, E. Hayakawa, T. Nishino, Y. Takeda, A. K. Mukhopadhyay, S. Garg, S. K. Bhattacharya, G. B. Nair, and M. Nishibuchi, J. Clin. Microbiol. 35:3150-3155, 1997). Evidence supporting a hypothesis that this clone emerged only recently and is spreading to many countries was obtained in this study. Of 227 strains isolated in a hospital in Bangladesh between 1977 and 1998, only 22 strains isolated between 1996 and 1998 belonged to the new O3:K6 clone (defined by the serovar, the tdh and trh typing, and AP-PCR profiles). The O3:K6 strains isolated from clinical sources in Taiwan, Laos, Japan, Thailand, Korea, and the United States between 1997 and 1998 were also shown to belong to the new O3:K6 clone. The clonality of the new O3:K6 strains was also confirmed by analysis of the toxRS sequence, which has been shown to be useful for phylogenetic analysis of the members of the genus Vibrio. The toxRS sequences of the representative strains of the new O3:K6 clone differed from those of the O3:K6 strains isolated before 1995 at least at 7 base positions within a 1,346-bp region. A new PCR method targeted to 2 of the base positions unique to the new O3:K6 clone was developed. This PCR method could clearly differentiate all 172 strains belonging to the new O3:K6 clone from other O3:K6 strains isolated earlier. One hundred sixty-six strains belonging to 28 serovars other than O3:K6 were also examined by the new PCR method. The tdh-positive and trh-lacking strains that belonged to the O4:K68 and O1:K untypeable serovars and were isolated in three countries and from international travelers beginning in 1997 gave positive results. The AP-PCR profiles of these strains were nearly identical to those of the new O3:K6 clone, and their toxRS sequences were 100% identical to that of the new O3:K6 clone. The results suggest that these strains may have diverged from the new O3:K6 clone by alteration of the O:K antigens. In conclusion, this study presents strong evidence for the first pandemicity in the history of V. parahaemolyticus and reports a novel toxRS-targeted PCR method that will be useful in epidemiological investigation of the cases associated with the current pandemic spread.     

The complete genomic sequence analysis of genotype 4 human astrovirus HASTVgz01 strain in Guangzhou

To analyze the genomic molecular structure and genotype of human astrovirus isolated from infant in Guangzhou of China, the primers were designed based on the genomic sequence of astrovirus from the GenBank and the target sequence were amplified by RT-PCR. Then the PCR-products were cloned to T vector and sequenced. The genomic nucleotide sequences were analyzed by the programs CLUSTAL W and DNASTAR. It was found that the full genomic length of HASTVgz01 strain was 6721 bp and the ORFs were 6558 bp. The 5' and 3' UTR were 82 and 81 nucleotides. The genome included 3 open reading frames (ORFs): ORF1a, ORF1b and ORF2. The 5'-terminal ORF1a started at nucleotide 83 and extended to nucleotide 2845. ORFlb (nt 2785 to nt 4332) overlaped ORFla by 61 nucleotides. The 3'-terminal ORF2 began at nucleotide 4325 and terminated at nucleotide 6640. ORF2 had 2316 nucleotides. Compared with other astrovirus sequences in GenBank, the homology of the amino acid sequence of ORF2 of HASTVgz01 strain with that of serotype 4 was 93% . Homology with other serotypes ranged from 61% to 70% . The complete nucleotide sequence of astrovirus HASTVgz01 strain isolated from Guangzhou in China was 6721 bp in length, GenBank accession NO. DQ344027. Comparing the ORF2 of astrovirus HASTVgz01 with the known sequences of types 1-8 the highest homology was serotype 4 (93%). Comparative sequence analysis of the HASTVgz01 ORF2 with the reported human astrovirus sequences revealed that the isolated astrovirus belongs to genotype (serotype) 4.     

Virulence gene- and pandemic group-specific marker profiling of clinical Vibrio parahaemolyticus isolates

Vibrio parahaemolyticus is a halophilic bacterium capable of causing food- and waterborne gastroenteritis, wound infections, and septicemia in humans. The organism has recently received increasing attention, as the emergence of a new clone, V. parahaemolyticus O3:K6, has resulted in the first documented pandemic spread of V. parahaemolyticus. We used microarray analyses to explore the presence of known virulence factors and genetic markers thought to be specific for V. parahaemolyticus O3:K6 and its clonal derivatives. Analyses of 48 human clinical isolates collected between 1997 and 2005 revealed that the V. parahaemolyticus chromosome 2 type III secretion system is not specifically associated with pandemic strains and can be found in tdh-negative (i.e., Kanagawa-negative) clinical isolates. These results highlight the genetic dynamism of V. parahaemolyticus and aid in refining the genetic definition of the pandemic group members.     

Complete genomic sequence of the lytic bacteriophage DT1 of Streptococcus thermophilus

Streptococcus thermophilus lytic bacteriophage DT1, isolated from a mozzarella whey, was characterized at the microbiological and molecular levels. Phage DT1 had an isometric head of 60 nm and a noncontractile tail of 260 x 8 nm, two major structural proteins of 26 and 32 kDa, and a linear double-stranded DNA genome with cohesive ends at its extremities. The host range of phage DT1 was limited to 5 of the 21 S. thermophilus strains tested. Using S. thermophilus SMQ-301 as a host, phage DT1 had a burst size of 276 +/- 36 and a latent period of 25 min. The genome of phage DT1 contained 34,820 bp with a GC content of 39.1%. Forty-six open reading frames (ORFs) of more than 40 codons were found and putative functions were assigned to 20 ORFs, mostly in the late region of phage DT1. Comparative genomic analysis of DT1 with the completely sequenced S. thermophilus temperate phage O1205 revealed two large homologous regions interspersed by two heterologous segments. The homologous regions consisted of the early replication genes, the late morphogenesis genes, and the lysis cassette. The divergent segments contained the DNA packaging machinery, the major structural proteins, and remnants of a lysogeny module.     

The complete DNA sequence and analysis of the virulence plasmid and of five additional plasmids carried by Shiga toxin-producing Escherichia coli O26:H11 strain H30

Shiga toxin-producing Escherichia coli (STEC) strains belonging to serogroup O26 have been associated with sporadic cases and outbreaks of hemorrhagic colitis and hemolytic uremic syndrome. In addition to chromosomal virulence genes, STEC strains usually harbor a large plasmid that carries genes associated with pathogenicity. The complete nucleotide sequence and genetic organization of 6 plasmids carried by STEC O26:H11 strain H30 were determined. The large virulence plasmid (pO26-Vir) was approximately 168 kb in size and contained 196 open reading frames (ORFs). pO26-Vir possesses a mosaic structure and shows similarity to the virulence plasmids in locus of enterocyte effacement (LEE)-negative STEC O113:H21 EH41 (pO113), in E. coli clinical strain C1096 (pSERB1), and in E. coli O157:H7 RIMD 0509952 (pO157). Plasmid pO26-Vir shares several highly conserved regions with pO157 and carries important virulence genes, including toxB, katP, espP, and the hly gene cluster. In addition, pO26-Vir possesses genes encoding for type IV pili (pilL-V). The second largest plasmid, pO26-L (73 kb) contains 101 ORFs. pO26-L carries the tetracycline resistance gene and has regions that show similarity to the E. coli conjugative resistance plasmid NR1. The third largest plasmid, pO26-S4 (5.8 kb), is homologous to the ColE2 colicinogenic plasmid that encodes for colicin E2. The remaining 3 plasmids, pO26-S1 (1.5 kb), pO26-S2 (3.1 kb), and pO26-S3 (4.2 kb), carry very little genetic information except for putative proteins involved in plasmid replication and DNA maintenance. The data presented underscore the diversity among the STEC virulence plasmids and provide insights into the evolution of these plasmids in STEC strains that cause serious human illness.     

Prevalence of the Pandemic Genotype of Vibrio parahaemolyticus in Dhaka, Bangladesh, and Significance of Its Distribution across Different Serotypes

Sixty-six strains of Vibrio parahaemolyticus belonging to 14 serotypes were isolated from hospitalized patients in Dhaka, Bangladesh, from January 1998 to December 2000. Among these, 48 strains belonging to four serotypes had the pandemic genotype and possessed the tdh gene. A marker (open reading frame ORF8) for a filamentous phage previously thought to correspond to the pandemic genotype was found to have a poor correlation with the pandemic genotype.     

Genetic characterization of DNA region containing the trh and ure genes of Vibrio parahaemolyticus

We have demonstrated that possession of the gene for thermostable direct hemolysin-related hemolysin (trh) coincides with the presence of the urease gene among clinical Vibrio parahaemolyticus strains and that the location of the two genes are in close proximity on the chromosome. Here, we cloned and sequenced the 15,754-bp DNA region containing the trh gene and the gene cluster for urease production from the chromosome of clinical V. parahaemolyticus (TH3996). We found 16 open reading frames (ORFs) and a lower G+C content (41%) compared with the total genome of this bacterium (46 to 47%). The ure cluster consisted of eight genes, namely, ureDABCEFG and ureR. ureR was located 5.2 kb upstream of the other seven genes in the opposite direction. The genetic organization and sequences of the ure genes resembled those found in Proteus mirabilis. Between ureR and the other ure genes, there were five ORFs, which are homologous with the nickel transport operon (nik) of Escherichia coli. We disrupted each of the ureR, ureC, and nikD genes in TH3996 by homologous recombination and analyzed the phenotype of the mutants. In the presence of urea these mutant strains had dramatically less urease activity than the strain they were derived from. Disruption of ureR, nikD, or ureC, however, had no effect on TRH production. The DNA region containing the trh, nik, and ure genes was found in only trh-positive strains and not in Kanagawa phenomenon-positive and environmental V. parahaemolyticus strains. At the end of the region, an insertion sequence-like element existed. These results suggest that the DNA region was introduced into V. parahaemolyticus in the past through a mechanism mediated by insertion sequences. This is the first reported case that the genes for an ATP-binding cassette-type nickel transport system, which may play a role in nickel transport through bacterial cytoplasmic membrane, are located adjacent to the ure cluster on the genome of an organism.     

Genomic analysis of bacteriophage PBECO4 infecting Escherichia coli O157:H7

Escherichia coli O157:H7 is a human pathogen. We isolated a novel bacteriophage infecting this bacterium from a sewage water treatment facility. Phage PBECO4 belongs to the family Myoviridae, having an isometric head and a contractile tail. It has a linear double-stranded DNA genome of 348,113 base pairs in length with a GC content of 34.09 %. Whole-genome sequencing revealed that PBECO4 is distantly related to enterobacteria phage vB_KleM_RaK2, with 10 % similarity, and Cronobacter phage vB_CsaM_GAP32 with 6 % similarity. Five hundred fifty-one putative open reading frames (ORFs) and six tRNA genes were found. Eight ORFs are related to genes encoding structural proteins, nine to DNA packaging, two to DNA lysis activity, and 42 to replication and regulation. Four hundred ninety ORFs have not been functionally annotated.