Comparative genomic analysis of Vibrio cholerae: genes that correlate with cholera endemic and pandemic disease

How big is the iceberg of which organellar genes in nuclear genomes are but the tip?

As more and more complete bacterial and archaeal genome sequences become available, the role of lateral gene transfer (LGT) in shaping them becomes more and more clear. Over the long term, it may be the dominant force, affecting most genes in most prokaryotes. We review the history of LGT, suggesting reasons why its prevalence and impact were so long dismissed. We discuss various methods purporting to measure the extent of LGT, and evidence for and against the notion that there is a core of never-exchanged genes shared by all genomes, from which we can deduce the "true" organismal tree. We also consider evidence for, and implications of, LGT between prokaryotes and phagocytic eukaryotes.

Bacterial infections: small intestine and colon

Bacterial infections of the small intestine and colon represent a major health problem for developing and developed nations. Recent technological progress has helped research groups to obtain important information on bacterial structure, identify evolutionary relationship between bacterial species, and learn details of the mechanisms involved in the interplay between host and microbes that culminate in disease expression. It is hoped that accumulated knowledge from in vitro experiments and animal models will translate into clinical benefit by means of developing new therapeutic strategies and effective vaccines.

Role of Vibrio cholerae O139 surface polysaccharides in intestinal colonization

Since the first occurrence of O139 Vibrio cholerae as a cause of cholera epidemics, this serogroup has been investigated intensively, and it has been found that its pathogenicity is comparable to that of O1 El Tor strains. O139 isolates express a thin capsule, composed of a polymer of repeating units structurally identical to the lipopolysaccharide (LPS) O side chain. In this study, we investigated the role of LPS O side chain and capsular polysaccharide (CPS) in intestinal colonization by with genetically engineered mutants. We constructed CPS-negative, CPS/LPS O side chain-negative, and CPS-positive/LPS O side chain-negative mutants. Furthermore, we constructed two mutants with defects in LPS core oligosaccharide (OS) assembly. Loss of LPS O side chain or CPS resulted in a approximately 30-fold reduction in colonization of the infant mouse small intestine, indicating that the presence of both LPS O side chain and CPS is important during the colonization process. The strain lacking both CPS and LPS O side chain and a CPS-positive, LPS O side chain-negative core OS mutant were both essentially unable to colonize. To characterize the role of surface polysaccharides in survival in the host intestine, resistance to several antimicrobial substances was investigated in vitro. These investigations revealed that the presence of CPS protects the cell against attack of the complement system and that an intact core OS is necessary for survival in the presence of bile.

Characterization of a novel Vibrio pathogenicity island (VPI-2) encoding neuraminidase (nanH) among toxigenic Vibrio cholerae isolates

Acquisition of virulence genes encoded on mobile genetic elements has played an important role in the emergence of pathogenic isolates of Vibrio cholerae, the causative agent of the diarrhoeal disease cholera. The genes encoding cholera toxin (ctxAB), the main cause of profuse secretory diarrhoea in cholera, are encoded on a filamentous bacteriophage CTXphi. The toxin coregulated pilus (TCP), an essential intestinal colonization factor, was originally designated as part of a pathogenicity island named the Vibrio pathogenicity island (VPI), but this island has more recently been proposed to be the genome of a filamentous phage, VPIphi. In this study, it is shown that nanH, which encodes neuraminidase, maps within a novel pathogenicity island designated VPI-2. The 57.3 kb VPI-2 has all of the characteristic features of a pathogenicity island, including the presence of a bacteriophage-like integrase (int), insertion in a tRNA gene (serine) and the presence of direct repeats at the chromosomal integration sites. Additionally, the G+C content of VPI-2 (42 mol%) is considerably lower than that of the entire genome (47 mol%). VPI-2 encodes several gene clusters, such as a restriction modification system (hsdR and hsdM) and genes required for the utilization of amino sugars (nan-nag region) as well as neuraminidase. To determine the distribution of VPI-2 among V. cholerae, 78 natural isolates were examined using PCR and Southern hybridization analysis for the presence of this region. All toxigenic V. cholerae O1 serogroup isolates examined contained VPI-2, whereas non-toxigenic isolates lacked the island. Of 14 V. cholerae O139 serogroup isolates examined, only one strain, MO2, contained the entire 57.3 kb island, whereas 13 O139 isolates contained only a 20.0 kb region with most of the 5' region of VPI-2 which included nanH deleted in these strains.

Competitive hybridization on spotted microarrays as a tool to conduct comparative genomic analyses of Xylella fastidiosa strains

Xylella fastidiosa strains are responsible for several plant diseases and since such isolates display a broad host range and complex biological behavior, genomic comparisons employing microarray hybridizations may provide an effective method to compare them. Thus, we performed a thorough validation of this type of approach using two recently sequenced strains of this phytopathogen. By matching microarray hybridization results to direct sequence comparisons, we were able to establish precise cutoff ratios for common and exclusive sequences, allowing the identification of exclusive genes involved in important biological traits. This validation will enable the use of microarray-based comparisons across a wide variety of microorganisms

Genomic profiles of clinical and environmental isolates of Vibrio cholerae O1 in cholera-endemic areas of Bangladesh

Diversity, relatedness, and ecological interactions of toxigenic Vibrio cholerae O1 populations in two distinctive habitats, the human intestine and the aquatic environment, were analyzed. Twenty environmental isolates and 42 clinical isolates were selected for study by matching serotype, geographic location of isolation in Bangladesh, and season of isolation. Genetic profiling was done by enterobacterial repetitive intergenic consensus sequence-PCR, optimized for profiling by using the fully sequenced V. cholerae El Tor N16961 genome. Five significant clonal clusters of haplotypes were found from 57 electrophoretic types. Isolates from different areas or habitats intermingled in two of the five significant clusters. Frequencies of haplotypes differed significantly only between the environmental populations (exact test; P < 0.05). Analysis of molecular variance yielded a population genetic structure reflecting the differentiating effects of geographic area, habitat, and sampling time. Although a parameter confounding the latter differences explained 9% of the total molecular variance in the entire population (P < 0.01), the net effect of habitat and time could not be separated because of the small number of environmental isolates included in the study. Five subpopulations from a single area were determined, and from these we were able to estimate a relative differentiating effect of habitat, which was small compared with the effect of temporal change. In conclusion, the resulting population structure supports the hypothesis that spatial and temporal fluctuations in the composition of toxigenic V. cholerae populations in the aquatic environment can cause shifts in the dynamics of the disease.

New variants of Vibrio cholerae O1 biotype El Tor with attributes of the classical biotype from hospitalized patients with acute diarrhea in Bangladesh

The sixth pandemic of cholera and, presumably, the earlier pandemics were caused by the classical biotype of Vibrio cholerae O1, which was progressively replaced by the El Tor biotype representing the seventh cholera pandemic. Although the classical biotype of V. cholerae O1 is extinct, even in southern Bangladesh, the last of the niches where this biotype prevailed, we have identified new varieties of V. cholerae O1, of the El Tor biotype with attributes of the classical biotype, from hospitalized patients with acute diarrhea in Bangladesh. Twenty-four strains of V. cholerae O1 isolated between 1991 and 1994 from hospitalized patients with acute diarrhea in Matlab, a rural area of Bangladesh, were examined for the phenotypic and genotypic traits that distinguish the two biotypes of V. cholerae O1. Standard reference strains of V. cholerae O1 belonging to the classical and El Tor biotypes were used as controls in all of the tests. The phenotypic traits commonly used to distinguish between the El Tor and classical biotypes, including polymyxin B sensitivity, chicken cell agglutination, type of tcpA and rstR genes, and restriction patterns of conserved rRNA genes (ribotypes), differentiated the 24 strains of toxigenic V. cholerae O1 into three types designated the Matlab types. Although all of the strains belonged to ribotypes that have been previously found among El Tor vibrios, type I strains had more traits of the classical biotype while type II and III strains appeared to be more like the El Tor biotype but had some classical biotype properties. These results suggest that, although the classical and El Tor biotypes have different lineages, there are possible naturally occurring genetic hybrids between the classical and El Tor biotypes that can cause cholera and thus provide new insight into the epidemiology of cholera in Bangladesh. Furthermore, the existence of such novel strains may have implications for the development of a cholera vaccine.

Primer on medical genomics. Part III: Microarray experiments and data analysis

Genomics has been defined as the comprehensive study of whole sets of genes, gene products, and their interactions as opposed to the study of single genes or proteins. Microarray technology is one of many novel tools that are allowing global and high-throughput analysis of genes and gene products. In addition to an introduction on underlying principles, the current review focuses on the use of both complementary DNA and oligodeoxynucleotide microarrays in gene expression analysis. Genome-wide experiments generate a massive amount of data points that require systematic methods of analysis to extract biologically useful information. Accordingly, the current educational communication discusses different methods of data analysis, including supervised and unsupervised clustering algorithms. Illustrative clinical examples show clinical applications, including (1) identification of candidate genes or pathological pathways (ie, elucidation of pathogenesis); (2) identification of "new" molecular classes of diseases that may be relevant in disease reclassification, prognostication, and treatment selection (ie, class discovery); and (3) use of expression profiles of known disease classes to predict diagnosis and classification of unknown samples (ie, class prediction). The current review should serve as an introduction to the subject for clinician investigators, physicians and medical scientists in training, practicing clinicians, and other students of medicine.

Genome sequence of a serotype M3 strain of group A Streptococcus: phage-encoded toxins, the high-virulence phenotype, and clone emergence

Genome sequences are available for many bacterial strains, but there has been little progress in using these data to understand the molecular basis of pathogen emergence and differences in strain virulence. Serotype M3 strains of group A Streptococcus (GAS) are a common cause of severe invasive infections with unusually high rates of morbidity and mortality. To gain insight into the molecular basis of this high-virulence phenotype, we sequenced the genome of strain MGAS315, an organism isolated from a patient with streptococcal toxic shock syndrome. The genome is composed of 1,900,521 bp, and it shares approximately 1.7 Mb of related genetic material with genomes of serotype M1 and M18 strains. Phage-like elements account for the great majority of variation in gene content relative to the sequenced M1 and M18 strains. Recombination produces chimeric phages and strains with previously uncharacterized arrays of virulence factor genes. Strain MGAS315 has phage genes that encode proteins likely to contribute to pathogenesis, such as streptococcal pyrogenic exotoxin A (SpeA) and SpeK, streptococcal superantigen (SSA), and a previously uncharacterized phospholipase A(2) (designated Sla). Infected humans had anti-SpeK, -SSA, and -Sla antibodies, indicating that these GAS proteins are made in vivo. SpeK and SSA were pyrogenic and toxic for rabbits. Serotype M3 strains with the phage-encoded speK and sla genes increased dramatically in frequency late in the 20th century, commensurate with the rise in invasive disease caused by M3 organisms. Taken together, the results show that phage-mediated recombination has played a critical role in the emergence of a new, unusually virulent clone of serotype M3 GAS.

Evolutionary genomics of Salmonella: gene acquisitions revealed by microarray analysis

The presence of homologues of Salmonella enterica sv. Typhimurium LT2 genes was assessed in 22 other Salmonella including members of all seven subspecies and Salmonella bongori. Genomes were hybridized to a microarray of over 97% of the 4,596 annotated ORFs in the LT2 genome. A phylogenetic tree based on homologue content, relative to LT2, was largely concordant with previous studies using sequence information from several loci. Based on the topology of this tree, homologues of genes in LT2 acquired by various clades were predicted including 513 homologues acquired by the ancestor of all Salmonella, 111 acquired by S. enterica, 105 by diphasic Salmonella, and 216 by subspecies 1, most of which are of unknown function. Because this subspecies is responsible for almost all Salmonella infections of mammals and birds, these genes will be of particular interest for further mechanistic studies. Overall, a high level of gene gain, loss, or rapid divergence was predicted along all lineages. For example, at least 425 close homologues of LT2 genes may have been laterally transferred into Salmonella and then between Salmonella lineages.

Redefining bacterial populations: a post-genomic reformation

Sexual reproduction and recombination are essential for the survival of most eukaryotic populations. Until recently, the impact of these processes on the structure of bacterial populations has been largely overlooked. The advent of large-scale whole-genome sequencing and the concomitant development of molecular tools, such as microarray technology, facilitate the sensitive detection of recombination events in bacteria. These techniques are revealing that bacterial populations are comprised of isolates that show a surprisingly wide spectrum of genetic diversity at the DNA level. Our new awareness of this genetic diversity is increasing our understanding of population structures and of how these affect host pathogen relationships.

Quorum-sensing regulators control virulence gene expression in Vibrio cholerae

The production of virulence factors including cholera toxin and the toxin-coregulated pilus in the human pathogen Vibrio cholerae is strongly influenced by environmental conditions. The well-characterized ToxR signal transduction cascade is responsible for sensing and integrating the environmental information and controlling the virulence regulon. We show here that, in addition to the known components of the ToxR signaling circuit, quorum-sensing regulators are involved in regulation of V. cholerae virulence. We focused on the regulators LuxO and HapR because homologues of these two proteins control quorum sensing in the closely related luminous marine bacterium Vibrio harveyi. Using an infant mouse model, we found that a luxO mutant is severely defective in colonization of the small intestine. Gene arrays were used to profile transcription in the V. cholerae wild type and the luxO mutant. These studies revealed that the ToxR regulon is repressed in the luxO mutant, and that this effect is mediated by another negative regulator, HapR. We show that LuxO represses hapR expression early in log-phase growth, and constitutive expression of hapR blocks ToxR-regulon expression. Additionally, LuxO and HapR regulate a variety of other cellular processes including motility, protease production, and biofilm formation. Together these data suggest a role for quorum sensing in modulating expression of blocks of virulence genes in a reciprocal fashion in vivo.

Improved analytical methods for microarray-based genome-composition analysis

BACKGROUND

Whereas genome sequencing has given us high-resolution pictures of many different species of bacteria, microarrays provide a means of obtaining information on genome composition for many strains of a given species. Genome-composition analysis using microarrays, or 'genomotyping', can be used to categorize genes into 'present' and 'divergent' categories based on the level of hybridization signal. This typically involves selecting a signal value that is used as a cutoff to discriminate present (high signal) and divergent (low signal) genes. Current methodology uses empirical determination of cutoffs for classification into these categories, but this methodology is subject to several problems that can result in the misclassification of many genes.

RESULTS

We describe a method that depends on the shape of the signal-ratio distribution and does not require empirical determination of a cutoff. Moreover, the cutoff is determined on an array-to-array basis, accounting for variation in strain composition and hybridization quality. The algorithm also provides an estimate of the probability that any given gene is present, which provides a measure of confidence in the categorical assignments.

CONCLUSIONS

Many genes previously classified as present using static methods are in fact divergent on the basis of microarray signal; this is corrected by our algorithm. We have reassigned hundreds of genes from previous genomotyping studies of Helicobacter pylori and Campylobacter jejuni strains, and expect that the algorithm should be widely applicable to genomotyping data.

Upstream – News in Genomics

In recent months a bumper crop of genomes has been completed, including the fission yeast (Schizosaccharomyces pombe) and rice (Oryza sativa). Two large-scale studies ofSaccharomyces cerevisiaeprotein complexes provided a picture of the eukaryotic proteome as a network of complexes. Amongst the other stories of interest was a demonstration that proteomic analysis of blood samples can be used to detect ovarian cancer, perhaps even as early as stage I.