Use of representational difference analysis to identify genomic differences between pathogenic strains of Vibrio cholerae

High quality reference genomes for toxigenic and non-toxigenic Vibrio cholerae serogroup O139

Toxigenic Vibrio cholerae of the O139 serogroup have been responsible for several large cholera epidemics in South Asia, and continue to be of clinical and historical significance today. This serogroup was initially feared to represent a new, emerging V. cholerae clone that would lead to an eighth cholera pandemic. However, these concerns were ultimately unfounded. The majority of clinically relevant V. cholerae O139 isolates are closely related to serogroup O1, biotype El Tor V. cholerae, and comprise a single sublineage of the seventh pandemic El Tor lineage. Although related, these V. cholerae serogroups differ in several fundamental ways, in terms of their O-antigen, capsulation phenotype, and the genomic islands found on their chromosomes. Here, we present four complete, high-quality genomes for V. cholerae O139, obtained using long-read sequencing. Three of these sequences are from toxigenic V. cholerae, and one is from a bacterium which, although classified serologically as V. cholerae O139, lacks the CTXφ bacteriophage and the ability to produce cholera toxin. We highlight fundamental genomic differences between these isolates, the V. cholerae O1 reference strain N16961, and the prototypical O139 strain MO10. These sequences are an important resource for the scientific community, and will improve greatly our ability to perform genomic analyses of non-O1 V. cholerae in the future. These genomes also offer new insights into the biology of a V. cholerae serogroup that, from a genomic perspective, is poorly understood.

An improved suppression subtractive hybridization technique to develop species-specific repetitive sequences from Erianthus arundinaceus (Saccharum complex)

BACKGROUND

Sugarcane has recently attracted increased attention for its potential as a source of bioethanol and methane. However, a narrow genetic base has limited germplasm enhancement of sugarcane. Erianthus arundinaceus is an important wild genetic resource that has many excellent traits for improving cultivated sugarcane via wide hybridization. Species-specific repetitive sequences are useful for identifying genome components and investigating chromosome inheritance in noblization between sugarcane and E. arundinaceus. Here, suppression subtractive hybridization (SSH) targeting E. arundinaceus-specific repetitive sequences was performed. The five critical components of the SSH reaction system, including enzyme digestion of genomic DNA (gDNA), adapters, digested gDNA concentrations, primer concentrations, and LA Taq polymerase concentrations, were improved using a stepwise optimization method to establish a SSH system suitable for obtaining E. arundinaceus-specific gDNA fragments.

RESULTS

Specificity of up to 85.42% was confirmed for the SSH method as measured by reverse dot blot (RDB) of an E. arundinaceus subtractive library. Furthermore, various repetitive sequences were obtained from the E. arundinaceus subtractive library via fluorescence in situ hybridization (FISH), including subtelomeric and centromeric regions. EaCEN2-166F/R and EaSUB1-127F/R primers were then designed as species-specific markers to accurately validate E. arundinaceus authenticity.

CONCLUSIONS

This is the first report that E. arundinaceus-specific repetitive sequences were obtained via an improved SSH method. These results suggested that this novel SSH system could facilitate screening of species-specific repetitive sequences for species identification and provide a basis for development of similar applications for other plant species.

Identification and distribution of putative virulence genes in clinical strains of Yersinia enterocolitica biovar 1A by suppression subtractive hybridization

AIMS

To detect putative virulence genes in clinical strains of Yersinia enterocolitica biovar 1A by suppression subtractive hybridization between two closely related strains of clinical and nonclinical origin having the same serotype (O:6,30-6,31).

METHODS AND RESULTS

Suppression Subtractive Hybridization (SSH) was used to identify genomic differences between clinical (serotype O:6,30-6,31, from diarrhoeic human stools) and nonclinical (serotype O:6,30-6,31, from wastewater) strains of Y. enterocolitica biovar 1A. Following genomic subtraction and DNA sequencing, nine DNA sequences that were present only in clinical biovar 1A strains were identified. The sequences identified using SSH showed similarity to conserved hypothetical proteins, proteins related to iron acquisition and haemin storage, type 1 secretion proteins, flagellar hook proteins, exported protein and ABC transport system. All these sequences showed high similarity with Y. enterocolitica 8081 (biovar 1B). The distribution of these genes was further analysed using PCR in 26 clinical strains of Y. enterocolitica biovar 1A. The results revealed that the distribution of these genes was not uniform.

CONCLUSIONS

Genes related to iron acquisition and storage, and flagellar proteins might be responsible for virulence of some of the clinical strains of Y. enterocolitica biovar 1A.

SIGNIFICANCE AND IMPACT OF THE STUDY

Genes identified in this study might be useful in understanding the pathogenic potential of clinical strains of Y. enterocolitica biovar 1A.

Application of representational difference analysis to identify genomic differences between Bradyrhizobium elkanii and B. Japonicum species

Bradyrhizobium elkanii is successfully used in the formulation of commercial inoculants and, together with B. japonicum, it fully supplies the plant nitrogen demands. Despite the similarity between B. japonicum and B. elkanii species, several works demonstrated genetic and physiological differences between them. In this work Representational Difference Analysis (RDA) was used for genomic comparison between B. elkanii SEMIA 587, a crop inoculant strain, and B. japonicum USDA 110, a reference strain. Two hundred sequences were obtained. From these, 46 sequences belonged exclusively to the genome of B. elkanii strain, and 154 showed similarity to sequences from B. japonicum genome. From the 46 sequences with no similarity to sequences from B. japonicum, 39 showed no similarity to sequences in public databases and seven showed similarity to sequences of genes coding for known proteins. These seven sequences were divided in three groups: similar to sequences from other Bradyrhizobium strains, similar to sequences from other nitrogen-fixing bacteria, and similar to sequences from non nitrogen-fixing bacteria. These new sequences could be used as DNA markers in order to investigate the rates of genetic material gain and loss in natural Bradyrhizobium strains.

Screening of strain-specific Actinobacillus pleuropneumoniae genes using a combination method

We describe a three-step method designed to identify distinct antigen-coding genes between two related bacterial genomes by: (a) constructing a subtractive library using Representational Difference Analysis (RDA), (b) characterization of gene expression in vitro using a ribosome display system combined with antibody screening and (c) gene recovery and confirmation using RT-PCR and reverse Southern hybridization, respectively. To test the efficacy of this strategy we screened the antigen-coding gene profile of Actinobacillus pleuropneumoniae (APP) strains CCVC259 and CCVC263 that do not elicit cross-protective immunity. This strategy identified six different DNA fragments from CCVC259 and 10 different DNA fragments from CCVC263. Of six sequences identified from CCVC259, 2 were not significantly similar, two were 74% and 87% homologous to the sequences encoding for the Ralstonia eutropha H16 conserved membrane protein and transcriptional regulator respectively, and two were >96% homologous to the Pseudomonas alcaligenes putative transposase subunit genes IS1474 and IS1475. Among ten unique DNA fragments identified from strain CCVC263, eight were homologous to DNA fragments encoding the TBP 1 precursor, ATP-dependent helicase HepA, glycosylase, methyltransferase and GTPase in the APP L20 genome and two genes identified had no significant similarity. Our findings indicated that the three-step method could be utilized to identify unique antigen-coding genes and may be a powerful and efficient technique for serotype-specific identification of pathogens and polyvalent vaccine design.

Microbial ecology in the age of genomics and metagenomics: concepts, tools, and recent advances

Microbial ecology examines the diversity and activity of micro-organisms in Earth's biosphere. In the last 20 years, the application of genomics tools have revolutionized microbial ecological studies and drastically expanded our view on the previously underappreciated microbial world. This review first introduces the basic concepts in microbial ecology and the main genomics methods that have been used to examine natural microbial populations and communities. In the ensuing three specific sections, the applications of the genomics in microbial ecological research are highlighted. The first describes the widespread application of multilocus sequence typing and representational difference analysis in studying genetic variation within microbial species. Such investigations have identified that migration, horizontal gene transfer and recombination are common in natural microbial populations and that microbial strains can be highly variable in genome size and gene content. The second section highlights and summarizes the use of four specific genomics methods (phylogenetic analysis of ribosomal RNA, DNA-DNA re-association kinetics, metagenomics, and micro-arrays) in analysing the diversity and potential activity of microbial populations and communities from a variety of terrestrial and aquatic environments. Such analyses have identified many unexpected phylogenetic lineages in viruses, bacteria, archaea, and microbial eukaryotes. Functional analyses of environmental DNA also revealed highly prevalent, but previously unknown, metabolic processes in natural microbial communities. In the third section, the ecological implications of sequenced microbial genomes are briefly discussed. Comparative analyses of prokaryotic genomic sequences suggest the importance of ecology in determining microbial genome size and gene content. The significant variability in genome size and gene content among strains and species of prokaryotes indicate the highly fluid nature of prokaryotic genomes, a result consistent with those from multilocus sequence typing and representational difference analyses. The integration of various levels of ecological analyses coupled to the application and further development of high throughput technologies are accelerating the pace of discovery in microbial ecology.

Hybridization monitor: A method for identifying differences between complex genomes

We have developed a method to identify and amplify differential fragments between two complex genomes. This technique, named hybridization-monitored genome differential analysis (HMDA), incorporates a monitor system into a PCR-based solid subtraction hybridization that tracks the entire hybridization process. This is achieved by monitoring the subtraction progress using PCR analysis of the conserved sequence of 18S rDNA in the tester sample after each round of subtraction. Homologous fragments can then be eliminated when bound to the driver DNA immobilized on a solid membrane. The hybridization continues until the conserved DNA sequence of 18S rDNA can no longer be detected, and most of the unbound DNA fragments left in the liquid were mainly the tester-specific fragments, thus greatly decreasing the complexity of DNA template of PCR amplification, increasing the amplification efficiency of differences accordingly, and ensuring high positive efficiency and coverage across the tester genome. We have applied the technique in a comparison between the genomes of Saccharomyces cerevisiae and Schizosaccharomyces pombe, which are two completely sequenced organisms. Results indicated that 95% of the subtracted clones have been confirmed to be different to the driver analyzed using the BLASTN homology alignment. With this technique, 240-fold enrichment of differences is obtained, and the coverage of the difference is up to 79%. These results indicate that HMDA can efficiently identify sequences that differ between two complex genomes.

Novel DNA sequences from natural strains of the nitrogen-fixing symbiotic bacterium Sinorhizobium meliloti

Variation in genome size and content is common among bacterial strains. Identifying these naturally occurring differences can accelerate our understanding of bacterial attributes, such as ecological specialization and genome evolution. In this study, we used representational difference analysis to identify potentially novel sequences not present in the sequenced laboratory strain Rm1021 of the nitrogen-fixing bacterium Sinorhizobium meliloti. Using strain Rm1021 as the driver and the type strain of S. meliloti ATCC 9930, which has a genome size approximately 370 kilobases bigger than that of strain Rm1021, as the tester, we identified several groups of sequences in the ATCC 9930 genome not present in strain Rm1021. Among the 85 novel DNA fragments examined, 55 showed no obvious homologs anywhere in the public databases. Of the remaining 30 sequences, 24 contained homologs to the Rm1021 genome as well as unique segments not found in Rm1021, 3 contained sequences homologous to those published for another S. meliloti strain but absent in Rm1021, 2 contained sequences homologous to other symbiotic nitrogen-fixing bacteria (Rhizobium etli and Bradyrhizobium japonicum), and 1 contained a sequence homologous to a gene in a non-nitrogen-fixing species, Pseudomonas sp. NK87. Using PCR, we assayed the distribution of 12 of the above 85 novel sequences in a collection of 59 natural S. meliloti strains. The distribution varied widely among the 12 novel DNA fragments, from 1.7% to 72.9%. No apparent correlation was found between the distribution of these novel DNA sequences and their genotypes obtained using multilocus enzyme electrophoresis. Our results suggest potentially high rates of gene gain and loss in S. meliloti genomes.

Identification of Vibrio parahaemolyticus pandemic group-specific DNA sequence by genomic subtraction

A genomic subtraction between a pandemic Vibrio parahaemolyticus and a nonpandemic strain that seemed to be clonally related was performed. A subtractive DNA fragment was identified to be a part of a 16-kbp insertion sequence which was present in almost all pandemic strains but not in nonpandemic strains tested.

Microplate subtractive hybridization to enrich for bacteroidales genetic markers for fecal source identification

The ability to identify sources of fecal pollution plays a key role in the analysis of human health risk and the implementation of water resource management strategies. One approach to this problem involves the identification of bacterial lineages or gene sequences that are found exclusively in a particular host species or group. We used subtractive hybridization to enrich for target host-specific fecal Bacteroidales rRNA gene fragments that were different from those of very closely related reference (subtracter) host sources. Target host rRNA gene fragments were hybridized to subtracter rRNA gene fragments immobilized in a microplate well, and target sequences that did not hybridize were cloned and sequenced for PCR primer design. The use of microplates for DNA immobilization resulted in a one-step subtractive hybridization in which the products could be directly amplified with PCR. The new host-specific primers designed from subtracted target fragments differentiated among very closely related Bacteroidales rRNA gene sequences and distinguished between similar fecal sources, such as elk and cow or human and domestic pet (dog).

Subtractive hybridization reveals a high genetic diversity in the fish pathogen Photobacterium damselae subsp. piscicida: evidence of a SXT-like element

Photobacterium damselae subsp. piscicida is the causative agent of fish pasteurellosis, a severe disease affecting cultured marine fish worldwide. In this study, suppression subtractive hybridization was used to identify DNA fragments present in the virulent strain PC554.2, but absent in the avirulent strain EPOY 8803-II. Twenty-one genomic regions of this type (that included twenty-six distinct putative ORFs) were analysed by DNA sequencing. Twenty ORFs encoded proteins with homology to proteins in other bacteria, including four homologues involved in siderophore biosynthesis, and four homologues related to mobile elements; three of these were putative transposases and one was a putative conjugative transposon related to the Vibrio cholerae SXT element. This sequence was shown to be integrated into a prfC gene homologue. Six ORFs showed no significant homology to known bacterial proteins. Among the 21 DNA fragments specific to strain PC554.2, 5 DNA fragments (representing 7 ORFs) were also absent in the avirulent strain ATCC 29690. The analysis of these differential regions, as well as the screening of their presence in a collection of strains, demonstrated the high genetic heterogeneity of this pathogen.

A glimpse into the expanded genome content of Vibrio cholerae through identification of genes present in environmental strains

Vibrio cholerae has multiple survival strategies which are reflected both in its broad distribution in many aquatic environments and its high genotypic diversity. To obtain additional information regarding the content of the V. cholerae genome, suppression subtractive hybridization (SSH) was used to prepare libraries of DNA sequences from two southern California coastal isolates which are divergent or absent in the clinical strain V. cholerae O1 El Tor N16961. More than 1,400 subtracted clones were sequenced. This revealed the presence of novel sequences encoding functions related to cell surface structures, transport, metabolism, signal transduction, luminescence, mobile elements, stress resistance, and virulence. Flanking sequence information was determined for loci of interest, and the distribution of these sequences was assessed for a collection of V. cholerae strains obtained from southern California and Mexican environments. This led to the surprising observation that sequences related to the toxin genes toxA, cnf1, and exoY are widespread and more common in these strains than those of the cholera toxin genes which are a hallmark of the pandemic strains of V. cholerae. Gene transfer among these strains could be facilitated by a 4.9-kbp plasmid discovered in one isolate, which possesses similarity to plasmids from other environmental vibrios. By investigating some of the nucleotide sequence basis for V. cholerae genotypic diversity, DNA fragments have been uncovered which could promote survival in coastal environments. Furthermore, a set of genes has been described which could be involved in as yet undiscovered interactions between V. cholerae and eukaryotic organisms.

Representational difference analysis: critical appraisal and method development for the identification of unique DNA sequences from prokaryotes

Representational difference analysis (RDA) has great potential for preferential amplification of unique but uncharacterised DNA sequences present in one source such as a whole genome, but absent from a related genome or other complex population of sequences. While a few examples of its successful exploitation have been published, the method has not been well dissected and robust, detailed published protocols are lacking. Here we examine the method in detail, suggest improvements and provide a protocol that has yielded key unique sequences from a pathogenic bacterial genome.

Novel virulence-associated type II secretion system unique to high-pathogenicity Yersinia enterocolitica

Yersinia enterocolitica strains comprise an important group of bacterial enteropathogens that cause a broad range of gastrointestinal syndromes. Three groups are distinguishable within this bacterial species, namely, the nonpathogenic group (biotype 1A strains), the low-pathogenicity, non-mouse-lethal group (biotypes 2 to 5), and the high-pathogenicity, mouse-lethal group (biotype 1B). To date, the presence of the high-pathogenicity island (HPI), a chromosomal locus that encodes the yersiniabactin system (involved in iron uptake), defines essentially the difference between low-pathogenicity and high-pathogenicity Y. enterocolitica strains, with the low-pathogenicity strains lacking the HPI. Using the powerful tool of representational difference analysis between the nonpathogenic 1A strain, NF-O, and its high-pathogenicity 1B counterpart, WA-314, we have identified a novel type II secretion gene cluster (yts1C-S) occurring exclusively in the high-pathogenicity group. The encoded secreton, designated Yts1 (for Yersinia type II secretion 1) was shown to be important for virulence in mice. A close examination of the almost completed genome sequence of another high-pathogenicity representative, Y. enterocolitica 8081, revealed a second putative type II secretion cluster uniformly distributed among all Y. enterocolitica isolates. This putative species-specific cluster (designated yts2) differed significantly from yts1, while resembling more closely the putative type II cluster present on the genome of Y. pestis. The Yts1 secreton thus appears to have been additionally acquired by the high-pathogenicity assemblage for a virulence-associated function.

Peptide ligands in antibacterial drug discovery: use as inhibitors in target validation and target-based screening

There is an urgent need to develop novel classes of antibiotics to counter the inexorable rise of resistant bacterial pathogens. Modern antibacterial drug discovery is focused on the identification and validation of novel protein targets that may have a suitable therapeutic index. In combination with assays for function, the advent of microbial genomics has been invaluable in identifying novel antibacterial drug targets. The major challenge in this field is the implementation of methods that validate protein targets leading to the discovery of new chemical entities. Ligand-directed drug discovery has the distinct advantage of having a concurrent analysis of both the importance of a target in the disease process and its amenability to functional modulation by small molecules. VITA is a process that enables a target-based paradigm by using peptide ligands for direct in vitro and in vivo validation of antibacterial targets and the implementation of high-throughput assays to identify novel inhibitory molecules. This process can establish sufficient levels of confidence indicating that the target is relevant to the disease process and inhibition of the target will lead to effective disease treatment.

Identification of a type III secretion system in uropathogenic Escherichia coli

To determine virulence-related genes in uropathogenic Escherichia coli (UPEC) showing invasiveness to T-24 bladder cancer cells, genomic subtractive hybridization was performed between a highly invasive and a less invasive strain. Forty-nine DNA fragments were isolated from the invasive strain. One of them showed homology with Salmonella invA gene. By chromosomal walking of the strain, a type III secretion system that has been described in E. coli O157:H7 was identified on the genome of the invasive strains. Three strains out of 100 UPEC isolates had a type III secretion system inserted at 64 min of the chromosome, corresponding to E. coli K-12 MG1655. This finding suggested that the type III secretion system could play a part in uropathogenicity of UPEC.

Identification of gene expression differences between Globodera pallida and G.'mexicana' by suppression subtractive hybridization

Summary Globodera pallida and G.'mexicana' are closely related nematode species that can mate and form viable hybrids on tomato but usually develop on different Solanaceous plants. Identification of nematode genes involved in parasitism is important for elucidation of disease resistance mechanisms in plants. In this study, we have used suppression subtractive hybridization (SSH) to investigate differences between the transcriptomes of G. pallida and G.'mexicana' J2s. This provides a basis for further studies characterizing pathogenicity factors in these nematodes. None of the cDNA fragments isolated in the SSH experiments appeared to be completely absent from the other transcriptome. Differences in expression levels of some of the isolated cDNAs between the two species were detected. Sequence analysis revealed that nearly 85% of the cloned sequences are nematode specific and a high proportion were pioneer genes for which no putative homologues were present in the databases. However, homologues of a cellulase and a putative pathogenicity factor previously described from G. rostochiensis were isolated. The putative roles of these sequences in parasitism are discussed.

Identification of virulence genes in a pathogenic strain of Pseudomonas aeruginosa by representational difference analysis

Pseudomonas aeruginosa is an opportunistic pathogen that may cause severe infections in humans and other vertebrates. In addition, a human clinical isolate of P. aeruginosa, strain PA14, also causes disease in a variety of nonvertebrate hosts, including plants, Caenorhabditis elegans, and the greater wax moth, Galleria mellonella. This has led to the development of a multihost pathogenesis system in which plants, nematodes, and insects have been used as adjuncts to animal models for the identification of P. aeruginosa virulence factors. Another approach to identifying virulence genes in bacteria is to take advantage of the natural differences in pathogenicity between isolates of the same species and to use a subtractive hybridization technique to recover relevant genomic differences. The sequenced strain of P. aeruginosa, strain PAO1, has substantial differences in virulence from strain PA14 in several of the multihost models of pathogenicity, and we have utilized the technique of representational difference analysis (RDA) to directly identify genomic differences between P. aeruginosa strains PA14 and PAO1. We have found that the pilC, pilA, and uvrD genes in strain PA14 differ substantially from their counterparts in strain PAO1. In addition, we have recovered a gene homologous to the ybtQ gene from Yersinia, which is specifically present in strain PA14 but absent in strain PAO1. Mutation of the ybtQ homolog in P. aeruginosa strain PA14 significantly attenuates the virulence of this strain in both G. mellonella and a burned mouse model of sepsis to levels comparable to those seen with PAO1. This suggests that the increased virulence of P. aeruginosa strain PA14 compared to PAO1 may relate to specific genomic differences identifiable by RDA.

A subtractive hybridisation analysis of genomic differences between the uropathogenic E. coli strain 536 and the E. coli K-12 strain MG1655

Suppression subtractive hybridisation (SSH) was performed to identify genomic differences between the uropathogenic Escherichia coli strain 536 and the non-pathogenic E. coli K-12 strain MG1655. In total, 22 DNA fragments were isolated which were specific for strain 536. Five of these fragments showed homology to known virulence determinants and four fragments matched genes for lipopolysaccharide (LPS) or capsule biosynthesis and a siderophore receptor. Seven fragments did not show any homology to known genes. These fragments may represent parts of putative pathogenicity islands (PAIs). Whereas two fragments were highly specific for uropathogenic E. coli (UPEC), the other fragments could also be detected among the other tested wild-type strains.

Use of representational difference analysis to identify Escherichia coli O157-specific DNA sequences

Whereas several important virulence factors in Escherichia coli O157 have been identified, studies suggest they are not always essential and are probably insufficient to account for the severe clinical manifestation of E. coli O157 infection. Identification of putative virulence determinants is crucial to the understanding of bacterial pathogenesis and genomic comparison analysis may aid the characterisation of unidentified virulence attributes. In this study, representational difference analysis (RDA) was used for genomic comparison of E. coli O157 with the proposed ancestral strain, E. coli O55. Unique E. coli O157 gene sequences were isolated and one, termed RDA-1, taken forward for further analysis. Southern blotting with labelled RDA-1 as a probe showed it to be present in 77% of E. coli O157 isolates and absent in all non-E. coli O157 screened. Sequence flanking RDA-1 was obtained from a genomic clone identified by hybridisation, and contained an open reading frame predicted to encode a novel iron-regulated outer membrane protein.

Representational differential analysis detects amplification of satellite sequences in postweaning multisystemic wasting syndrome of pigs

Representational difference analysis (RDA) was used as a molecular approach to identify unique sequences associated with postweaning multisystemic wasting syndrome (PMWS) in pigs. Three rounds of subtractive hybridization and amplification between driver DNA extracted from normal pigs and tester DNA from PMWS-affected animals were performed. The final product corresponding to sequences associated with PMWS in pigs was analyzed using agarose gel electrophoresis, and 9 fragments were visualized after staining with ethidium bromide. Eight recombinants were successively cloned and sequenced, and the results were then compared with existing databases. Most of the PMWS clones isolated were satellite sequences from pig centrometric regions and 1 was a microsatellite sequence. One clone represented a microsatellite sequence, and 2 clones showed no homology with any gene found in the databases. The sequence comparison data did not reveal any homology with an infectious agent such as a virus or a bacterium. In the present experimental setting, it was concluded that PMWS in pigs triggers molecular changes such as an amplification of genomic regions containing repeated sequences.

Genetic diversity of Vibrio cholerae in Chesapeake Bay determined by amplified fragment length polymorphism fingerprinting

Vibrio cholerae is indigenous to the aquatic environment, and serotype non-O1 strains are readily isolated from coastal waters. However, in comparison with intensive studies of the O1 group, relatively little effort has been made to analyze the population structure and molecular evolution of non-O1 V. cholerae. In this study, high-resolution genomic DNA fingerprinting, amplified fragment length polymorphism (AFLP), was used to characterize the temporal and spatial genetic diversity of 67 V. cholerae strains isolated from Chesapeake Bay during April through July 1998, at four different sampling sites. Isolation of V. cholerae during the winter months (January through March) was unsuccessful, as observed in earlier studies (J. H. L. Kaper, R. R. Colwell, and S. W. Joseph, Appl. Environ. Microbiol. 37:91-103, 1979). AFLP fingerprints subjected to similarity analysis yielded a grouping of isolates into three large clusters, reflecting time of the year when the strains were isolated. April and May isolates were closely related, while July isolates were genetically diverse and did not cluster with the isolates obtained earlier in the year. The results suggest that the population structure of V. cholerae undergoes a shift in genotype that is linked to changes in environmental conditions. From January to July, the water temperature increased from 3 degrees C to 27.5 degrees C, bacterial direct counts increased nearly an order of magnitude, and the chlorophyll a concentration tripled (or even quadrupled at some sites). No correlation was observed between genetic similarity among isolates and geographical source of isolation, since isolates found at a single sampling site were genetically diverse and genetically identical isolates were found at several of the sampling sites. Thus, V. cholerae populations may be transported by surface currents throughout the entire Bay, or, more likely, similar environmental conditions may be selected for a specific genotype. The dynamic nature of the population structure of this bacterial species in Chesapeake Bay provides new insight into the ecology and molecular evolution of V. cholerae in the natural environment.

Identification of a vibrio cholerae RTX toxin gene cluster that is tightly linked to the cholera toxin prophage

We identify and characterize a gene cluster in El Tor Vibrio cholerae that encodes a cytotoxic activity for HEp-2 cells in vitro. This gene cluster contains four genes and is physically linked to the cholera toxin (CTX) element in the V. cholerae genome. We demonstrate by using insertional mutagenesis that this gene cluster is required for the cytotoxic activity. The toxin, RtxA, resembles members of the RTX (repeats in toxin) toxin family in that it contains a GD-rich repeated motif. Like other RTX toxins, its activity depends on an activator, RtxC, and an associated ABC transporter system, RtxB and RtxD. In V. cholerae strains of the classical biotype, a deletion within the gene cluster removes rtxC and eliminates cytotoxic activity. Other strains, including those of the current cholera pandemic, contain a functional gene cluster and display cytotoxic activity. Thus, the RTX gene cluster in El Tor O1 and O139 strains might have contributed significantly to their emergence. Furthermore, the RTX toxin of V. cholerae may be associated with residual adverse properties displayed by certain live, attenuated cholera vaccines.

PCR-based subtractive hybridization and differences in gene content among strains of Helicobacter pylori

Genes that are characteristic of only certain strains of a bacterial species can be of great biologic interest. Here we describe a PCR-based subtractive hybridization method for efficiently detecting such DNAs and apply it to the gastric pathogen Helicobacter pylori. Eighteen DNAs specific to a monkey-colonizing strain (J166) were obtained by subtractive hybridization against an unrelated strain whose genome has been fully sequenced (26695). Seven J166-specific clones had no DNA sequence match to the 26695 genome, and 11 other clones were mixed, with adjacent patches that did and did not match any sequences in 26695. At the protein level, seven clones had homology to putative DNA restriction-modification enzymes, and two had homology to putative metabolic enzymes. Nine others had no database match with proteins of assigned function. PCR tests of 13 unrelated H. pylori strains by using primers specific for 12 subtracted clones and complementary Southern blot hybridizations indicated that these DNAs are highly polymorphic in the H. pylori population, with each strain yielding a different pattern of gene-specific PCR amplification. The search for polymorphic DNAs, as described here, should help identify previously unknown virulence genes in pathogens and provide new insights into microbial genetic diversity and evolution.