Evolutionary genomics of Salmonella: gene acquisitions revealed by microarray analysis

Comparative genome analysis of Lactobacillus reuteri and Lactobacillus fermentum reveal a genomic island for reuterin and cobalamin production

Lactobacillus reuteri is a heterofermentative lactic acid bacterium that naturally inhabits the gut of humans and other animals. The probiotic effects of L. reuteri have been proposed to be largely associated with the production of the broad-spectrum antimicrobial compound reuterin during anaerobic metabolism of glycerol. We determined the complete genome sequences of the reuterin-producing L. reuteri JCM 1112(T) and its closely related species Lactobacillus fermentum IFO 3956. Both are in the same phylogenetic group within the genus Lactobacillus. Comparative genome analysis revealed that L. reuteri JCM 1112(T) has a unique cluster of 58 genes for the biosynthesis of reuterin and cobalamin (vitamin B(12)). The 58-gene cluster has a lower GC content and is apparently inserted into the conserved region, suggesting that the cluster represents a genomic island acquired from an anomalous source. Two-dimensional nuclear magnetic resonance (2D-NMR) with (13)C(3)-glycerol demonstrated that L. reuteri JCM 1112(T) could convert glycerol to reuterin in vivo, substantiating the potential of L. reuteri JCM 1112(T) to produce reuterin in the intestine. Given that glycerol is shown to be naturally present in feces, the acquired ability to produce reuterin and cobalamin is an adaptive evolutionary response that likely contributes to the probiotic properties of L. reuteri.

Do Salmonella carry spare tyres?

Salmonellae are enterobacteria that have the unique ability to change their flagellar composition by switching expression among two loci that encode the major flagellin protein. This property is not available to all Salmonella, but is species, subspecies and serotype specific. Curiously, the subsequent loss of the second locus in some lineages of Salmonella has apparently been tolerated and, indeed, has led to considerable success for some lineages. We discuss here an evolutionary model for maintenance of this unique function and the possible evolutionary advantages of loss or preservation of this mechanism. We hypothesize that the second flagellin locus is a genetic 'spare tyre' used in particular environmental circumstances.

RosE represses Std fimbrial expression in Salmonella enterica serotype Typhimurium

The Salmonella enterica serotype Typhimurium (S. typhimurium) genome contains a large repertoire of putative fimbrial operons that remain poorly characterized because they are not expressed in vitro. In this study, insertions that induced expression of the putative stdABCD fimbrial operon were identified from a random bank of transposon mutants by screening with immuno-magnetic particles for ligand expression (SIMPLE). Transposon insertions upstream of csgC and lrhA or within dam, setB and STM4463 (renamed rosE) resulted in expression of StdA and its assembly into fimbrial filaments on the cell surface. RosE is a novel negative regulator of Std fimbrial expression as indicated by its repression of a std::lacZ reporter construct and by binding of the purified protein to a DNA region upstream of the stdA start codon. Expression of Std fimbriae in the rosE mutant resulted in increased attachment of S. typhimurium to human colonic epithelial cell lines (T-84 and CaCo-2). A rosE mutant exhibited a reduced ability to compete with virulent S. typhimurium for colonization of murine organs, while no defect was observed when both competing strains carried a stdAB deletion. These data suggest that a tight control of Std fimbrial expression mediated by RosE is required during host pathogen interaction.

Human Salmonella clinical isolates distinct from those of animal origin

The global trend toward intensive livestock production has led to significant public health risks and industry-associated losses due to an increased incidence of disease and contamination of livestock-derived food products. A potential factor contributing to these health concerns is the prospect that selective pressure within a particular host may give rise to bacterial strain variants that exhibit enhanced fitness in the present host relative to that in the parental host from which the strain was derived. Here, we assessed 184 Salmonella enterica human and animal clinical isolates for their virulence capacities in mice and for the presence of the Salmonella virulence plasmid encoding the SpvB actin cytotoxin required for systemic survival and Pef fimbriae, implicated in adherence to the murine intestinal epithelium. All (21 of 21) serovar Typhimurium clinical isolates derived from animals were virulent in mice, whereas many (16 of 41) serovar Typhimurium isolates derived from human salmonellosis patients lacked this capacity. Additionally, many (10 of 29) serovar Typhimurium isolates derived from gastroenteritis patients did not possess the Salmonella virulence plasmid, in contrast to all animal and human bacteremia isolates tested. Lastly, among serovar Typhimurium isolates that harbored the Salmonella virulence plasmid, 6 of 31 derived from human salmonellosis patients were avirulent in mice, which is in contrast to the virulent phenotype exhibited by all the animal isolates examined. These studies suggest that Salmonella isolates derived from human salmonellosis patients are distinct from those of animal origin. The characterization of these bacterial strain variants may provide insight into their relative pathogenicities as well as into the development of treatment and prophylactic strategies for salmonellosis.

A process for analysis of microarray comparative genomics hybridisation studies for bacterial genomes

Background

Microarray based comparative genomic hybridisation (CGH) experiments have been used to study numerous biological problems including understanding genome plasticity in pathogenic bacteria. Typically such experiments produce large data sets that are difficult for biologists to handle. Although there are some programmes available for interpretation of bacterial transcriptomics data and CGH microarray data for looking at genetic stability in oncogenes, there are none specifically to understand the mosaic nature of bacterial genomes. Consequently a bottle neck still persists in accurate processing and mathematical analysis of these data. To address this shortfall we have produced a simple and robust CGH microarray data analysis process that may be automated in the future to understand bacterial genomic diversity.

Results

The process involves five steps: cleaning, normalisation, estimating gene presence and absence or divergence, validation, and analysis of data from test against three reference strains simultaneously. Each stage of the process is described and we have compared a number of methods available for characterising bacterial genomic diversity, for calculating the cut-off between gene presence and absence or divergence, and shown that a simple dynamic approach using a kernel density estimator performed better than both established, as well as a more sophisticated mixture modelling technique. We have also shown that current methods commonly used for CGH microarray analysis in tumour and cancer cell lines are not appropriate for analysing our data.

Conclusion

After carrying out the analysis and validation for three sequenced Escherichia coli strains, CGH microarray data from 19 E. coli O157 pathogenic test strains were used to demonstrate the benefits of applying this simple and robust process to CGH microarray studies using bacterial genomes.

Salmonellae interplay with host cells

Salmonellae are important causes of enteric diseases in all vertebrates. Characterization of the molecular mechanisms that underpin the interactions of salmonellae with their animal hosts has advanced greatly over the past decade, mainly through the study of Salmonella enterica serovar Typhimurium in tissue culture and animal models of infection. Knowledge of these bacterial processes and host responses has painted a dynamic and complex picture of the interaction between salmonellae and animal cells. This Review focuses on the molecular mechanisms of these host-pathogen interactions, in terms of their context, significance and future perspectives.

Coordinate regulation of Salmonella pathogenicity island 1 (SPI1) and SPI4 in Salmonella enterica serovar Typhimurium

Salmonella enterica serovar Typhimurium harbors five pathogenicity islands (SPI) required for infection in vertebrate hosts. Although the role of SPI1 in promoting epithelial invasion and proinflammatory cell death has been amply documented, SPI4 has only more recently been implicated in Salmonella virulence. SPI4 is a 24-kb pathogenicity island containing six open reading frames, siiA to siiF. Secretion of the 595-kDa SiiE protein requires a type I secretory system encoded by siiC, siiD, and siiF. An operon polarity suppressor (ops) sequence within the 5' untranslated region upstream of siiA is required for optimal SPI4 expression and predicted to bind the antiterminator RfaH. SiiE concentrations are decreased in a SPI1 mutant strain, suggesting that SPI1 and SPI4 may have common regulatory inputs. SPI1 gene expression is positively regulated by the transcriptional activators HilA, HilC, and HilD, encoded within SPI1, and negatively regulated by the regulators HilE and PhoP. Here, we show that mutations in hilA, hilC, or hilD similarly reduce expression of siiE, and mutations in hilE or phoP enhance siiE expression. Individual overexpression of HilA, HilC, or HilD in the absence of SPI1 cannot activate siiE expression, suggesting that these transcriptional regulators act in concert or in combination with additional SPI1-encoded regulatory loci to activate SPI4. HilA is no longer required for siiE expression in an hns mutant strain, suggesting that HilA promotes SPI4 expression by antagonizing the global transcriptional silencer H-NS. Coordinate regulation suggests that SPI1 and SPI4 play complementary roles in the interaction of S. enterica serovar Typhimurium with the host intestinal mucosa.

Evolution of the chaperone/usher assembly pathway: fimbrial classification goes Greek

Many Proteobacteria use the chaperone/usher pathway to assemble proteinaceous filaments on the bacterial surface. These filaments can curl into fimbrial or nonfimbrial surface structures (e.g., a capsule or spore coat). This article reviews the phylogeny of operons belonging to the chaperone/usher assembly class to explore the utility of establishing a scheme for subdividing them into clades of phylogenetically related gene clusters. Based on usher amino acid sequence comparisons, our analysis shows that the chaperone/usher assembly class is subdivided into six major phylogenetic clades, which we have termed alpha-, beta-, gamma-, kappa-, pi-, and sigma-fimbriae. Members of each clade share related operon structures and encode fimbrial subunits with similar protein domains. The proposed classification system offers a simple and convenient method for assigning newly discovered chaperone/usher systems to one of the six major phylogenetic groups.

The complete genome sequence and comparative genome analysis of the high pathogenicity Yersinia enterocolitica strain 8081

The human enteropathogen, Yersinia enterocolitica, is a significant link in the range of Yersinia pathologies extending from mild gastroenteritis to bubonic plague. Comparison at the genomic level is a key step in our understanding of the genetic basis for this pathogenicity spectrum. Here we report the genome of Y. enterocolitica strain 8081 (serotype 0:8; biotype 1B) and extensive microarray data relating to the genetic diversity of the Y. enterocolitica species. Our analysis reveals that the genome of Y. enterocolitica strain 8081 is a patchwork of horizontally acquired genetic loci, including a plasticity zone of 199 kb containing an extraordinarily high density of virulence genes. Microarray analysis has provided insights into species-specific Y. enterocolitica gene functions and the intraspecies differences between the high, low, and nonpathogenic Y. enterocolitica biotypes. Through comparative genome sequence analysis we provide new information on the evolution of the Yersinia. We identify numerous loci that represent ancestral clusters of genes potentially important in enteric survival and pathogenesis, which have been lost or are in the process of being lost, in the other sequenced Yersinia lineages. Our analysis also highlights large metabolic operons in Y. enterocolitica that are absent in the related enteropathogen, Yersinia pseudotuberculosis, indicating major differences in niche and nutrients used within the mammalian gut. These include clusters directing, the production of hydrogenases, tetrathionate respiration, cobalamin synthesis, and propanediol utilisation. Along with ancestral gene clusters, the genome of Y. enterocolitica has revealed species-specific and enteropathogen-specific loci. This has provided important insights into the pathology of this bacterium and, more broadly, into the evolution of the genus. Moreover, wider investigations looking at the patterns of gene loss and gain in the Yersinia have highlighted common themes in the genome evolution of other human enteropathogens.

The goal of this study was to catalogue all the genes encoded within the Y. enterocolitica genome to help us better understand how this bacterium and related bacteria cause different diseases. There are currently genome sequences (complete gene catalogues) available for two other members of this bacterial lineage, which cause dramatically different diseases: Y. pseudotuberculosis, like Y. enterocolitica, is a gut pathogen (enteropathogen) causing gastroenteritis in humans and animals. Yersinia pestis mostly resides within blood (circulating or in fleas following blood meals) and lymph tissue. It causes bubonic plague in humans and animals, and is historically known as “The Black Death.” A three-way comparison of these genomes revealed a patchwork of genes we have defined as being species- or disease-specific and genes that are common to all three Yersinia species. This has provided us with important information on shared gene functions that define the two enteropathogenic yersinias and those that differentiate them. This will help us to connect what we know about the Y. enterocolitica lifestyle within the gut to the disease it causes and its genetic makeup. We have also provided further evidence of gene-loss by Y. pestis as it has evolved from Y. pseudotuberculosis into a more acute systemic pathogen. Similar patterns of gene loss are seen in other important pathogens such as Salmonella enterica serovar Typhi.

Species of environmental mycobacteria differ in their abilities to grow in human, mouse, and carp macrophages and with regard to the presence of mycobacterial virulence genes, as observed by DNA microarray hybridization

There are many species of environmental mycobacteria (EM) that infect animals that are important to the economy and research and that also have zoonotic potential. The genomes of very few of these bacterial species have been sequenced, and little is known about the molecular mechanisms by which most of these opportunistic pathogens cause disease. In this study, 18 isolates of EM isolated from fish and humans (including strains of Mycobacterium avium, Mycobacterium peregrinum, Mycobacterium chelonae, and Mycobacterium salmoniphilum) were examined for their abilities to grow in macrophage lines from humans, mice, and carp. Genomic DNA from 14 of these isolates was then hybridized against DNA from an M. avium reference strain, with a custom microarray containing virulence genes of mycobacteria and a selection of representative genes from metabolic pathways. The strains of EM had different abilities to grow within the three types of cell lines, which grouped largely according to the host from which they were isolated. Genes identified as being putatively absent in some of the strains included those with response regulatory functions, cell wall compositions, and fatty acid metabolisms as well as a recently identified pathogenicity island important to macrophage uptake. Further understanding of the role these genes play in host specificity and pathogenicity will be important to gain insight into the zoonotic potential of certain EM as well as their mechanisms of virulence.

Variation in gene content among geographically diverse Sulfolobus isolates

The ability of competitive (i.e., comparative) genomic hybridization (CGH) to assess similarity across entire microbial genomes suggests that it should reveal diversification within and between natural populations of free-living prokaryotes. We used CGH to measure relatedness of genomes drawn from Sulfolobus populations that had been shown in a previous study to be diversified along geographical lines. Eight isolates representing a wide range of spatial separation were compared with respect to gene-specific tags based on a closely related reference strain (Sulfolobus solfataricus P2). For the purpose of assessing genetic divergence, 232 loci identified as polymorphic were assigned one of two alleles based on the corresponding fluorescence intensities from the arrays. Clustering of these binary genotypes was stable with respect to changes in the threshold and similarity criteria, and most of the groupings were consistent with an isolation-by-distance model of diversification. These results indicate that increasing spatial separation of geothermal sites correlates not only with minor sequence polymorphisms in conserved genes of Sulfolobus (demonstrated in the previous study), but also with the regions of difference (RDs) that occur between genomes of conspecifics. In view of the abundance of RDs in prokaryotic genomes and the relevance that some RDs may have for ecological adaptation, the results further suggest that CGH on microarrays may have advantages for investigating patterns of diversification in other free-living archaea and bacteria.

Cloning and transfer of the Salmonella pathogenicity island 2 type III secretion system for studies of a range of gram-negative genera

The engineering of bacterial strains with specific phenotypes frequently requires the use of blocks or "cassettes" of genes that act together to perform a desired function. The potential benefits of utilizing type III secretion systems in this regard are becoming increasingly realized since these systems can be used to direct interactions with host cells for beneficial purposes such as vaccine development, anticancer therapies, and targeted protein delivery. However, convenient methods to clone and transfer type III secretion systems for studies of a range of different types of bacteria are lacking. In addition to functional applications, such methods would also reveal important information about the evolution of a given type III secretion system, such as its ability to be expressed and functional outside of the strain of origin. We describe here the cloning of the Salmonella enterica serovar Typhimurium pathogenicity island 2 (SPI-2) type III secretion system onto a vector that can be easily transferred to a range of gram-negative bacterial genera. We found that expression of the cloned SPI-2 system in different Gammaproteobacteria and Alphaproteobacteria (as monitored by SseB protein levels) is dependent on the bacterial strain and growth medium. We also demonstrate that the cloned system is functional for secretion, can direct interactions with macrophages, and can be used as a novel tool to analyze the predicted interaction of SseB with host cells. This work provides a foundation for future applications where the cloned SPI-2 region (or other cloned type III systems) can provide a desired function to an engineered gram-negative strain.

Genetic diversity among major endemic strains of Leptospira interrogans in China

BACKGROUND

Leptospirosis is a world-widely distributed zoonosis. Humans become infected via exposure to pathogenic Leptospira spp. from contaminated water or soil. The availability of genomic sequences of Leptospira interrogans serovar Lai and serovar Copenhageni opened up opportunities to identify genetic diversity among different pathogenic strains of L. interrogans representing various kinds of serotypes (serogroups and serovars).

RESULTS

Comparative genomic hybridization (CGH) analysis was used to compare the gene content of L. interrogans serovar Lai strain Lai with that of other 10 L. interrogans strains prevailed in China and one identified from Brazil using a microarray spotted with 3,528 protein coding sequences (CDSs) of strain Lai. The cutoff ratio of sample/reference (S/R) hybridization for detecting the absence of genes from one tested strain was set by comparing the ratio of S/R hybridization and the in silico sequence similarities of strain Lai and serovar Copenhageni strain Fiocruz L1-130. Among the 11 strains tested, 275 CDSs were found absent from at least one strain. The common backbone of the L. interrogans genome was estimated to contain about 2,917 CDSs. The genes encoding fundamental cellular functions such as translation, energy production and conversion were conserved. While strain-specific genes include those that encode proteins related to either cell surface structures or carbohydrate transport and metabolism. We also found two genomic islands (GIs) in strain Lai containing genes divergently absent in other strains. Because genes encoding proteins with potential pathogenic functions are located within GIs, these elements might contribute to the variations in disease manifestation. Differences in genes involved in O-antigen biosynthesis were also identified for strains belonging to different serogroups, which offers an opportunity for future development of genomic typing tools for serological classification.

CONCLUSION

CGH analyses for pathogenic leptospiral strains prevailed in China against the L. interrogans serovar Lai strain Lai CDS-spotted microarrays revealed 2,917 common backbone CDSs and strain specific genes encoding proteins mainly related to cell surface structures and carbohydrated transport/metabolism. Of the 275 CDSs considered absent from at least one of the L. interrogans strains tested, most of them were clustered in the rfb gene cluster and two putative genomic islands (GI A and B) in strain Lai. The strain-specific genes detected via this work will provide a knowledge base for further investigating the pathogenesis of L interrogans and/or for the development of effective vaccines and/or diagnostic tools.

Comparative genomics of Streptomyces avermitilis, Streptomyces cattleya, Streptomyces maritimus and Kitasatospora aureofaciens using a Streptomyces coelicolor microarray system

DNA/DNA microarray hybridization was used to compare the genome content of Streptomyces avermitilis, Streptomyces cattleya, Streptomyces maritimus and Kitasatospora aureofaciens with that of Streptomyces coelicolor A3(2). The array data showed an about 93% agreement with the genome sequence data available for S. avermitilis and also showed a number of trends in the genome structure for Streptomyces and closely related Kitasatospora. A core central region was well conserved, which might be predicted from previous research and this was linked to a low degree of gene conservation in the terminal regions of the linear chromosome across all four species. Between these regions there are two areas of intermediate gene conservation by microarray analysis where gene synteny is still detectable in S. avermitilis. Nonetheless, a range of conserved genes could be identified within the terminal regions. Variation in the genes involved in differentiation, transcription, DNA replication, etc. provides interesting insights into which genes in these categories are generally conserved and which are not. The results also provide target priorities for possible gene knockouts in a group of bacteria with a very large numbers of genes with unknown functions compared to most bacterial species.

A Hidden Markov model web application for analysing bacterial genomotyping DNA microarray experiments

Whole genome DNA microarray genomotyping experiments compare the gene content of different species or strains of bacteria. A statistical approach to analysing the results of these experiments was developed, based on a Hidden Markov model (HMM), which takes adjacency of genes along the genome into account when calling genes present or absent. The model was implemented in the statistical language R and applied to three datasets. The method is numerically stable with good convergence properties. Error rates are reduced compared with approaches that ignore spatial information. Moreover, the HMM circumvents a problem encountered in a conventional analysis: determining the cut-off value to use to classify a gene as absent. An Apache Struts web interface for the R script was created for the benefit of users unfamiliar with R. The application may be found at http://hmmgd.cryst.bbk.ac.uk/hmmgd. The source code illustrating how to run R scripts from an Apache Struts-based web application is available from the corresponding author on request. The application is also available for local installation if required.

Insertion sequence-driven diversification creates a globally dispersed emerging multiresistant subspecies of E. faecium

Enterococcus faecium, an ubiquous colonizer of humans and animals, has evolved in the last 15 years from an avirulent commensal to the third most frequently isolated nosocomial pathogen among intensive care unit patients in the United States. E. faecium combines multidrug resistance with the potential of horizontal resistance gene transfer to even more pathogenic bacteria. Little is known about the evolution and virulence of E. faecium, and genomic studies are hampered by the absence of a completely annotated genome sequence. To further unravel its evolution, we used a mixed whole-genome microarray and hybridized 97 E. faecium isolates from different backgrounds (hospital outbreaks (n = 18), documented infections (n = 34) and asymptomatic carriage of hospitalized patients (n = 15), and healthy persons (n = 15) and animals (n = 21)). Supported by Bayesian posterior probabilities (PP = 1.0), a specific clade containing all outbreak-associated strains and 63% of clinical isolates was identified. Sequencing of 146 of 437 clade-specific inserts revealed mobile elements (n = 74), including insertion sequence (IS) elements (n = 42), phage genes (n = 6) and plasmid sequences (n = 26), hypothetical (n = 58) and membrane proteins (n = 10), and antibiotic resistance (n = 9) and regulatory genes (n = 11), mainly located on two contigs of the unfinished E. faecium DO genome. Split decomposition analysis, varying guanine cytosine content, and aberrant codon adaptation indices all supported acquisition of these genes through horizontal gene transfer with IS16 as the predicted most prominent insert (98% sensitive, 100% specific). These findings suggest that acquisition of IS elements has facilitated niche adaptation of a distinct E. faecium subpopulation by increasing its genome plasticity. Increased genome plasticity was supported by higher diversity indices (ratio of average genetic similarities of pulsed-field gel electrophoresis and multi locus sequence typing) for clade-specific isolates. Interestingly, the previously described multi locus sequence typing-based clonal complex 17 largely overlapped with this clade. The present data imply that the global emergence of E. faecium, as observed since 1990, represents the evolution of a subspecies with a presumably better adaptation than other E. faecium isolates to the constraints of a hospital environment.

Adaptation to the host environment: regulation of the SPI1 type III secretion system in Salmonella enterica serovar Typhimurium

Salmonella enterica invades the intestinal epithelium of the host using a type III secretion system encoded on Salmonella pathogenicity island 1 (SPI1). The bacteria integrate environmental signals from a variety of global regulatory systems to precisely induce transcription of SPI1. The regulatory circuit converges on expression of HilA, which directly regulates transcription of the SPI1 apparatus genes. Transcription of hilA is controlled by a complex feed-forward loop. Regulatory signals feed into the system through post-transcriptional and post-translational control of HilD, which in turn activates HilC and RtsA. These three regulators act in concert to control hilA transcription. The system acts as a switch, ensuring that SPI1 is fully on at the appropriate time.

Validation of mixed-genome microarrays as a method for genetic discrimination

Comparative genomic hybridizations have been used to examine genetic relationships among bacteria. The microarrays used in these experiments may have open reading frames from one or more reference strains (whole-genome microarrays), or they may be composed of random DNA fragments from a large number of strains (mixed-genome microarrays [MGMs]). In this work both experimental and virtual arrays are analyzed to assess the validity of genetic inferences from these experiments with a focus on MGMs. Empirical data are analyzed from an Enterococcus MGM, while a virtual MGM is constructed in silico using sequenced genomes (Streptococcus). On average, a small MGM is capable of correctly deriving phylogenetic relationships between seven species of Enterococcus with accuracies of 100% (n=100 probes) and 95% (n=46 probes); more probes are required for intraspecific differentiation. Compared to multilocus sequence methods and whole-genome microarrays, MGMs provide additional discrimination between closely related strains and offer the possibility of identifying unique strain or lineage markers. Representational bias can have mixed effects. Microarrays composed of probes from a single genome can be used to derive phylogenetic relationships, although branch length can be exaggerated for the reference strain. We describe a case where disproportional representation of different strains used to construct an MGM can result in inaccurate phylogenetic inferences, and we illustrate an algorithm that is capable of correcting this type of bias. The bias correction algorithm automatically provides bootstrap confidence values and can provide multiple bias-corrected trees with high confidence values.

Determination of the gene content of Salmonella genomes by microarray analysis

Microarray technology provides a convenient and relatively inexpensive way of investigating the genetic content of bacterial genomes by comparative genomic hybridization. In this method, genomic DNA of an unknown bacterial strain of interest and that of a closely related sequenced isolate are hybridized to the same array. Hybridization signals are subsequently translated into gene absence and presence predictions for the experimental strain. Our nonredundant microarray of PCR products representing almost all genes from a number of the sequenced Salmonella enterica serovars (including Typhimurium, Typhi, Paratyphi A, and Enteritidis) allows accurate predictions of gene presence and absence in hundreds of Salmonella isolates on whole genome scale, for a fraction of the cost of complete genome sequencing, or resequencing using tiled oligo-arrays.

Studying bacterial genome dynamics using microarray-based comparative genomic hybridization

Genome sequencing has revealed the remarkable amount of genetic diversity that can be encountered in bacterial genomes. In particular, the comparison of genome sequences from closely related strains has uncovered significant differences in gene content, hinting at the dynamic nature of bacterial genomes. The study of these genome dynamics is crucial to leveraging genomic information because the genome sequence of a single bacterial strain may not accurately represent the genome of the species. The dynamic nature of bacterial genome content has required us to apply the concepts of comparative genomics (CG) at the species level. Although direct genome sequence comparisons are an ideal method of performing CG, one current constraint is the limited availability of multiple genome sequences from a given bacterial species. DNA microarray-based comparative genomic hybridization (MCGH), which can be used to determine the presence or absence of thousands of genes in a single hybridization experiment, provides a powerful alternative for determining genome content and has been successfully used to investigate the genome dynamics of a wide number of bacterial species. Although MCGH-based studies have already provided a new vista on bacterial genome diversity, original methods for MCGH have been limited by the absence of novel gene sequences included in the microarray. New applications of the MCGH platform not only promise to accelerate the pace of novel gene discovery but will also help provide an integrated microarray-based approach to the study of bacterial CG.

Identification of putative ancestors of the multidrug-resistant Salmonella enterica serovar typhimurium DT104 clone harboring the Salmonella genomic island 1

The origin of multidrug-resistant Salmonella enterica serovar typhimurium (S. typhimurium) harboring the Salmonella Genomic Island 1 (SGI1), which was detected for the first time in the mid-1980s is unknown. In this study, we performed microarray genomotyping of four multidrug-resistant SGI1 positive strains and found that unlike the S. typhimurium LT2 strain, the multidrug-resistant strains lacked genes STM0517-0529 allowing the utilization of allantoin as a sole nitrogen source. We extended this observation by PCR screening of additional 120 S. typhimurium field strains and found that this locus was absent in all SGI1 positive and also in 24% of SGI1 negative strains, which were proposed to be the original recipients of SGI1. To prove this hypothesis, we compared the STM0517-0529 negative strains (with or without the SGI1) by PFGE and PCR prophage typing and found that 8 out of 11 of the SGI1 negative strains and 17 out of 22 SGI1 positive strains were of identical PFGE pattern and PCR prophage pattern, while this specific pattern was never observed among STM0517-0529 positive strains. We therefore propose that a lineage of the S. typhimurium DT104 sensitive strain first lost the ability to metabolize allantoin and then acquired SGI1.

Flagellin gene sequence evolution in Salmonella

Salmonella exhibits 70 serologically distinct flagellins, used internationally to diagnose and track infections. The terminal sequences of flagellin protein subunits are conserved in a range of bacteria and are here used as evolutionary markers to reveal how new serotypes arise. Terminal sequences of flagellins that exhibit factors g or m (G-group) were distinct from other Salmonella antigens (Non-G-group) and cluster more closely with Escherichia coli. It is postulated that G-group flagellins were inherited from a common ancestor of E. coli and Salmonella and that these antigens were among the original set in Salmonella. Sequence differences at the 5' termini may prevent recombination between co-infecting strains. Evidence of increased variation of flagellin in rare biphasic G-group serotypes suggests that the presence of a second flagellin locus allows mutation of the G-group flagellin. FljB probably arose from a single duplication of a Non-G gene, since which synonymous mutations resulted in the fljB-specific sequence at the 5' termini.

Altered levels of Salmonella DNA adenine methylase are associated with defects in gene expression, motility, flagellar synthesis, and bile resistance in the pathogenic strain 14028 but not in the laboratory strain LT2

Comparative genomic analysis has revealed limited strain diversity between Salmonella pathogenic and nonpathogenic isolates. Thus, some of the relative virulence and host-immune response disparities may be credited to differential gene regulation rather than gross differences in genomic content. Here we show that altered levels of Salmonella DNA adenine methylase (Dam) resulted in acute defects in virulence-associated gene expression, motility, flagellin synthesis, and bile resistance in the Salmonella pathogenic strain 14028 but not in avirulent laboratory strain LT2. The defects in motility exhibited by 14028 in response to altered Dam levels was not dependent on the presence of the regulatory protein, RpoS. The transitioning between flagellar types (phase variation) was also differentially regulated in 14028 versus LT2 in response to dam levels, resulting in distinct differences in flagellin expression states. These data suggest that differential gene regulation may contribute to the relative virulence disparities observed between Salmonella serovars that are closely related at the DNA level.

SseL is a salmonella-specific translocated effector integrated into the SsrB-controlled salmonella pathogenicity island 2 type III secretion system

Bacterial pathogens use horizontal gene transfer to acquire virulence factors that influence host colonization, alter virulence traits, and ultimately shape the outcome of disease following infection. One hallmark of the host-pathogen interaction is the prokaryotic type III secretion system that translocates virulence factors into host cells during infection. Salmonella enterica possesses two type III secretion systems that are utilized during host colonization and intracellular replication. Salmonella pathogenicity island 2 (SPI2) is a genomic island containing approximately 30 contiguous genes required to assemble a functional secretion system including the two-component regulatory system called SsrA-SsrB that positively regulates transcription of the secretion apparatus. We used transcriptional profiling with DNA microarrays to search for genes that coregulate with the SPI2 type III secretion machinery in an SsrB-dependent manner. Here we report the identification of a Salmonella-specific translocated effector called SseL that is required for full virulence during murine typhoid-like disease. Analysis of infected macrophages using fluorescence-activated cell sorting revealed that sseL is induced inside cells and requires SsrB for expression. SseL is retained predominantly in the cytoplasm of infected cells following translocation by the type III system encoded in SPI2. Animal infection experiments with sseL mutant bacteria indicate that integration of SseL into the SsrB response regulatory system contributes to systemic virulence of this pathogen.

Similarity and differences in the Lactobacillus acidophilus group identified by polyphasic analysis and comparative genomics

A set of lactobacilli were investigated by polyphasic analysis. Multilocus sequence analysis, DNA typing, microarray analysis, and in silico whole-genome alignments provided a remarkably consistent pattern of similarity within the Lactobacillus acidophilus complex. On microarray analysis, 17 and 5% of the genes from Lactobacillus johnsonii strain NCC533 represented variable and strain-specific genes, respectively, when tested against four independent isolates of L. johnsonii. When projected on the NCC533 genome map, about 10 large clusters of variable genes were identified, and they were enriched around the terminus of replication. A quarter of the variable genes and two-thirds of the strain-specific genes were associated with mobile DNA. Signatures for horizontal gene transfer and modular evolution were found in prophages and in DNA from the exopolysaccharide biosynthesis cluster. On microarray hybridizations, Lactobacillus gasseri strains showed a shift to significantly lower fluorescence intensities than the L. johnsonii test strains, and only genes encoding very conserved cellular functions from L. acidophilus hybridized to the L. johnsonii array. In-silico comparative genomics showed extensive protein sequence similarity and genome synteny of L. johnsonii with L. gasseri, L. acidophilus, and Lactobacillus delbrueckii; moderate synteny with Lactobacillus casei; and scattered X-type sharing of protein sequence identity with the other sequenced lactobacilli. The observation of a stepwise decrease in similarity between the members of the L. acidophilus group suggests a strong element of vertical evolution in a natural phylogenetic group. Modern whole-genome-based techniques are thus a useful adjunct to the clarification of taxonomical relationships in problematic bacterial groups.

Comparison of Salmonella enterica serovar Typhimurium LT2 and non-LT2 salmonella genomic sequences, and genotyping of salmonellae by using PCR

Genes of Salmonella enterica serovar Typhimurium LT2 expected to be specifically present in Salmonella were selected using the Basic Local Alignment Search Tool (BLAST) program. The 152 selected genes were compared with 11 genomic sequences of Salmonella serovars, including Salmonella enterica subsp. I and IIIb and Salmonella bongori (V), and were clustered into 17 groups by their comparison patterns. A total of 38 primer pairs were constructed to represent each of the 17 groups, and PCR was performed with various Salmonella subspecies including Salmonella enterica subsp. I, II, IIIa, IIIb, IV, VI, and V to evaluate a comprehensive DNA-based scheme for identification of Salmonella subspecies and the major disease-causing Salmonella serovars. Analysis of PCR results showed that Salmonella enterica subsp. I was critically divided from other subspecies, and Salmonella strains belonging to S. enterica subsp. I were clustered based on their serovars. In addition, genotypic relationships within S. enterica subsp. I by PCR results were investigated. Also, Salmonella signature genes, Salmonella enterica serovar Typhimurium signature genes, and Salmonella enterica subsp. I signature genes were demonstrated based on their PCR results. The described PCR method suggests a rapid and convenient method for identification of Salmonella serovars that can be used by nonspecialized laboratories. Genome sequence comparison can be a useful tool in epidemiologic and taxonomic studies of Salmonella.

Identification of specific gene sequences conserved in contemporary epidemic strains of Salmonella enterica

Genetic elements specific to recent and contemporary epidemic strains of Salmonella enterica were identified using comparative genomic analysis. Two epidemic multidrug-resistant (MDR) strains, MDR Salmonella enterica serovar Typhimurium definitive phage type 104 (DT104) and cephalosporin-resistant MDR Salmonella enterica serovar Newport, and an epidemic pansusceptible strain, Salmonella serovar Typhimurium DT160, were subjected to Salmonella gene microarray and suppression subtractive hybridization analyses. Their genome contents were compared with those of coexisting sporadic strains matched by serotype, geographic and temporal distribution, and host species origin. These paired comparisons revealed that epidemic strains of S. enterica had specific genes and gene regions that were shared by isolates of the same subtype. Most of these gene sequences are related to mobile genetic elements, including phages, plasmids, and plasmid-like and transposable elements, and some genes may encode proteins conferring growth or survival advantages. The emergence of epidemic MDR strains may therefore be associated with the presence of fitness-associated genetic factors in addition to their antimicrobial resistance genes.

Multiplex PCR-based method for identification of common clinical serotypes of Salmonella enterica subsp. enterica

A multiplex PCR method has been developed to differentiate between the most common clinical serotypes of Salmonella enterica subsp. enterica encountered in Washington State and the United States in general. Six genetic loci from S. enterica serovar Typhimurium and four from S. enterica serovar Typhi were used to create an assay consisting of two five-plex PCRs. The assays gave reproducible results with 30 different serotypes that represent the most common clinical isolates of S. enterica subsp. enterica. Of these, 22 serotypes gave unique amplification patterns compared with each other and the other 8 serotypes were grouped into four pairs. These were further resolved by two additional PCRs. We compared the data from PCR serotyping with conventional serotyping and found that PCR serotyping was nearly as discriminatory as conventional serotyping was. The results from a blind test screening 111 clinical isolates revealed that 97% were correctly identified using the multiplex PCR assay. The assay can be easily performed on multiple samples with final results in less than 5 h and, in conjunction with pulsed-field gel electrophoresis, forms a very robust test method for the molecular subtyping of Salmonella enterica subsp. enterica.

Comparative genetics of the rdar morphotype in Salmonella

The Salmonella rdar morphotype is a distinct, rough and dry colony morphology formed by the extracellular interaction of thin aggregative fimbriae (Tafi or curli), cellulose, and other polysaccharides. Cells in rdar colonies are more resistant to desiccation and exogenous stresses, which is hypothesized to aid in the passage of pathogenic Salmonella spp. between hosts. Here we analyzed the genetic and phenotypic conservation of the rdar morphotype throughout the entire Salmonella genus. The rdar morphotype was conserved in 90% of 80 isolates representing all 7 Salmonella groups; however, the frequency was only 31% in a reference set of 16 strains (Salmonella reference collection C [SARC]). Comparative gene expression analysis was used to separate cis- and trans-acting effects on promoter activity for the 16 SARC strains, focusing on the 780-bp intergenic region containing divergent promoters for the master regulator of the rdar morphotype (agfD) and the Tafi structural genes (agfB). Surprisingly, promoter functionality was conserved in most isolates, and loss of the phenotype was due primarily to defects in trans-acting regulatory factors. We hypothesize that trans differences have been caused by domestication, whereas cis differences, detected for Salmonella enterica subsp. arizonae isolates, may reflect an evolutionary change in lifestyle. Our results demonstrate that the rdar morphotype is conserved throughout the salmonellae, but they also emphasize that regulation is an important source of variability among isolates.

Molecular characterisation of Salmonella strains by an oligonucleotide multiprobe microarray

A DNA microarray has been developed for the simultaneous characterisation and typing of Salmonella enterica subsp. enterica isolates. One-hundred and nine 35-40 mer oligonucleotides probes detect flagellar and somatic antigen encoding genes (serogroup or serotype specific), important virulence genes located within or outside the pathogenicity islands, phage-associated genes and antibiotic resistance determinants. The probes were printed on glass slides and whole genomic Cy5-labelled Salmonella DNA was hybridised to the substrate. A set of 19 different Salmonella strains and one Escherichia coli strain has been selected as positive and negative controls for each probe. The validity of the results is confirmed by gene-specific PCRs or phenotypic methods (serotyping, MIC determination for various antimicrobial agents). Of 2071 data points generated, an agreement of 97.4% has been obtained between microarray and PCR/phenotypic results. Twenty-six data points (1.3%) were classified as uncertain and, similarly, 1.3% showed a discordant result. The microarray described here is a new tool to study the epidemiology of Salmonella strains on the genotypic level and might become a powerful method in risk assessment studies.

The use of comparative genomic hybridization to characterize genome dynamics and diversity among the serotypes of Shigella

Background

Compelling evidence indicates that Shigella species, the etiologic agents of bacillary dysentery, as well as enteroinvasive Escherichia coli, are derived from multiple origins of Escherichia coli and form a single pathovar. To further understand the genome diversity and virulence evolution of Shigella, comparative genomic hybridization microarray analysis was employed to compare the gene content of E. coli K-12 with those of 43 Shigella strains from all lineages.

Results

For the 43 strains subjected to CGH microarray analyses, the common backbone of the Shigella genome was estimated to contain more than 1,900 open reading frames (ORFs), with a mean number of 726 undetectable ORFs. The mosaic distribution of absent regions indicated that insertions and/or deletions have led to the highly diversified genomes of pathogenic strains.

Conclusion

These results support the hypothesis that by gain and loss of functions, Shigella species became successful human pathogens through convergent evolution from diverse genomic backgrounds. Moreover, we also found many specific differences between different lineages, providing a window into understanding bacterial speciation and taxonomic relationships.

Comparison of different approaches for comparative genetic analysis using microarray hybridization

A robust analysis of comparative genomic microarray data is critical for meaningful genomic comparison studies. In this paper, we compare our method (implemented in a new software tool, GENCOM, freely available at http://www.ifr.ac.uk/safety/gencom ) with three commonly used analysis methods: GACK (freely available at http://falkow.stanford.edu ), an empirical cut-off value of twofold difference between the fluorescence intensities after LOWESS normalization or after AVERAGE normalization in which the fluorescence intensity is divided by the average fluorescence intensity of the entire data set. Each method was tested using data sets from real experiments with prior knowledge of conserved and divergent genes. GENCOM and GACK were superior when a high proportion of genes were divergent. GENCOM was the most suitable method for the data set in which the relationship between the fluorescence intensities was not linear. GENCOM has proved robust in an analysis of all the data sets tested.

DNA microarray technology: a new tool for the epidemiological typing of bacterial pathogens?

Genomic hybridization on whole genome arrays detects the presence or absence of similar DNA regions in sufficiently related microorganisms, allowing genome-wide comparison of their genetic contents. A whole genome array is based on a sequenced bacterial isolate, and is a collection of DNA probes fixed on a solid support. In a single hybridization experiment, the absence/presence status of all genes of the sequenced microbe in the queried isolate can be examined. The objective of this minireview is to summarize the past usage of DNA microarray technology for microbial strain characterizations, and to estimate its future utilization in epidemiological studies and molecular typing of bacterial pathogens. The studies reviewed here confirm the usefulness of microarray technology for the detection of genetic polymorphisms. However, the construction or purchase of DNA microarrays and the performance of strain to strain hybridization experiments are still prohibitively expensive for routine application. Future use of arrays in epidemiology is likely to depend on the development of more cost-effective protocols, more robust and simplified formats, and the adequate evaluation of their performance (efficacy) and convenience (efficiency) compared with other genotyping methods. It seems more likely that a more focused assay, concentrating on genomic regions of variability previously detected by genome-wide microarrays, will find broad application in routine bacterial epidemiology.

Application of DNA microarray technology for detection, identification, and characterization of food-borne pathogens

DNA microarrays represent the latest advance in molecular technology. In combination with bioinformatics, they provide unparalleled opportunities for simultaneous detection of thousands of genes or target DNA sequences and offer tremendous potential for studying food-borne microorganisms. This review provides an up-to-date look at the application of DNA microarray technology to detect food-borne pathogenic bacteria, viruses, and parasites. In addition, it covers the advantages of using microarray technology to further characterize microorganisms by providing information for specific identification of isolates, to understand the pathogenesis based on the presence of virulence genes, and to indicate how new pathogenic strains evolved epidemiologically and phylogenetically.

Genetic diversity of the Q fever agent, Coxiella burnetii, assessed by microarray-based whole-genome comparisons

Coxiella burnetii, a gram-negative obligate intracellular bacterium, causes human Q fever and is considered a potential agent of bioterrorism. Distinct genomic groups of C. burnetii are revealed by restriction fragment-length polymorphisms (RFLP). Here we comprehensively define the genetic diversity of C. burnetii by hybridizing the genomes of 20 RFLP-grouped and four ungrouped isolates from disparate sources to a high-density custom Affymetrix GeneChip containing all open reading frames (ORFs) of the Nine Mile phase I (NMI) reference isolate. We confirmed the relatedness of RFLP-grouped isolates and showed that two ungrouped isolates represent distinct genomic groups. Isolates contained up to 20 genomic polymorphisms consisting of 1 to 18 ORFs each. These were mostly complete ORF deletions, although partial deletions, point mutations, and insertions were also identified. A total of 139 chromosomal and plasmid ORFs were polymorphic among all C. burnetii isolates, representing ca. 7% of the NMI coding capacity. Approximately 67% of all deleted ORFs were hypothetical, while 9% were annotated in NMI as nonfunctional (e.g., frameshifted). The remaining deleted ORFs were associated with diverse cellular functions. The only deletions associated with isogenic NMI variants of attenuated virulence were previously described large deletions containing genes involved in lipopolysaccharide (LPS) biosynthesis, suggesting that these polymorphisms alone are responsible for the lower virulence of these variants. Interestingly, a variant of the Australia QD isolate producing truncated LPS had no detectable deletions, indicating LPS truncation can occur via small genetic changes. Our results provide new insight into the genetic diversity and virulence potential of Coxiella species.

Screening for rapidly evolving genes in the ectomycorrhizal fungus Paxillus involutus using cDNA microarrays

We have examined the variations in gene content and sequence divergence that could be associated with symbiotic adaptations in the ectomycorrhizal fungus Paxillus involutus and the closely related species Paxillus filamentosus. Strains with various abilities to form mycorrhizae were analysed by comparative genomic hybridizations using a cDNA microarray containing 1076 putative unique genes of P. involutus. To screen for genes diverging at an enhanced and presumably non-neutral rate, we implemented a simple rate test using information from both the variations in hybridizations signal and data on sequence divergence of the arrayed genes relative to the genome of Coprinus cinereus. C. cinereus is a free-living saprophyte and is the closest evolutionary relative to P. involutus that has been fully sequenced. Approximately 17% of the genes investigated were detected as rapidly diverging within Paxillus. Furthermore, 6% of the genes varied in copy numbers between the analysed strains. Genome rearrangements associated with this variation including duplications and deletions may also play a role in adaptive evolution. The cohort of divergent and duplicated genes showed an over-representation of either orphans, genes whose products are located at membranes, or genes encoding for components of stress/defence reactions. Some of the identified genomic changes may be associated with the variation in host specificity of ectomycorrhizal fungi. The proposed procedure could be generally applicable to screen for rapidly evolving genes in closely related strains or species where at least one has been sequenced or characterized by expressed sequence tag analysis.

Comparative phylogenomics of pathogenic bacteria by microarray analysis

DNA microarrays represent a powerful technology that enables whole-scale comparison of bacterial genomes. This, coupled with new methods to model DNA microarray data, is facilitating the development of robust comparative phylogenomics analyses. Such studies have dramatically increased our ability to differentiate between bacteria, highlighting previously undetected genetic differences and population structures and providing new insight into virulence and evolution of bacterial pathogens. Recent results from such studies have generated insights into the evolution of bacterial pathogens, the levels of diversity and plasticity in the genome of a species, as well as the differences in virulence amongst pathogenic bacteria.

Identification of core and variable components of the Salmonella enterica subspecies I genome by microarray

We have performed microarray hybridization studies on 40 clinical isolates from 12 common serovars within Salmonella enterica subspecies I to identify the conserved chromosomal gene pool. We were able to separate the core invariant portion of the genome by a novel mathematical approach using a decision tree based on genes ranked by increasing variance. All genes within the core component were confirmed using available sequence and microarray information for S. enterica subspecies I strains. The majority of genes within the core component had conserved homologues in Escherichia coli K-12 strain MG1655. However, many genes present in the conserved set which were absent or highly divergent in K-12 had close homologues in pathogenic bacteria such as Shigella flexneri and Pseudomonas aeruginosa. Genes within previously established virulence determinants such as SPI1 to SPI5 were conserved. In addition several genes within SPI6, all of SPI9, and three fimbrial operons (fim, bcf, and stb) were conserved within all S. enterica strains included in this study. Although many phage and insertion sequence elements were missing from the core component, approximately half the pseudogenes present in S. enterica serovar Typhi were conserved. Furthermore, approximately half the genes conserved in the core set encoded hypothetical proteins. Separation of the core and variant gene sets within S.enterica subspecies I has offered fundamental biological insight into the genetic basis of phenotypic similarity and diversity across S. enterica subspecies I and shown how the core genome of these pathogens differs from the closely related E. coli K-12 laboratory strain.

A new experimental approach for studying bacterial genomic island evolution identifies island genes with bacterial host-specific expression patterns

Background

Genomic islands are regions of bacterial genomes that have been acquired by horizontal transfer and often contain blocks of genes that function together for specific processes. Recently, it has become clear that the impact of genomic islands on the evolution of different bacterial species is significant and represents a major force in establishing bacterial genomic variation. However, the study of genomic island evolution has been mostly performed at the sequence level using computer software or hybridization analysis to compare different bacterial genomic sequences. We describe here a novel experimental approach to study the evolution of species-specific bacterial genomic islands that identifies island genes that have evolved in such a way that they are differentially-expressed depending on the bacterial host background into which they are transferred.

Results

We demonstrate this approach by using a "test" genomic island that we have cloned from the Salmonella typhimurium genome (island 4305) and transferred to a range of Gram negative bacterial hosts of differing evolutionary relationships to S. typhimurium. Systematic analysis of the expression of the island genes in the different hosts compared to proper controls allowed identification of genes with genera-specific expression patterns. The data from the analysis can be arranged in a matrix to give an expression "array" of the island genes in the different bacterial backgrounds. A conserved 19-bp DNA site was found upstream of at least two of the differentially-expressed island genes. To our knowledge, this is the first systematic analysis of horizontally-transferred genomic island gene expression in a broad range of Gram negative hosts. We also present evidence in this study that the IS200 element found in island 4305 in S. typhimurium strain LT2 was inserted after the island had already been acquired by the S. typhimurium lineage and that this element is likely not involved in the integration or excision of island 4305.

Conclusion

The "clone-and-transfer" approach of evolutionary study identifies genes whose expression patterns indicate the existence of genera-specific regulatory mechanisms that influence the expression of horizontally-transferred DNA sections. The results provide key information that can be used to facilitate the identification of these regulatory mechanisms.

Differences in gene content between Salmonella enterica serovar Enteritidis isolates and comparison to closely related serovars Gallinarum and Dublin

Salmonella enterica serovar Enteritidis is often transmitted into the human food supply through eggs of hens that appear healthy. This pathogen became far more prevalent in poultry following eradication of the fowl pathogen S. enterica serovar Gallinarum in the mid-20th century. To investigate whether changes in serovar Enteritidis gene content contributed to this increased prevalence, and to evaluate genetic heterogeneity within the serovar, comparative genomic hybridization was performed on eight 60-year-old and nineteen 10- to 20-year-old serovar Enteritidis strains from various hosts, using a Salmonella-specific microarray. Overall, almost all the serovar Enteritidis genomes were very similar to each other. Excluding two rare strains classified as serovar Enteritidis in the Salmonella reference collection B, only eleven regions of the serovar Enteritidis phage type 4 (PT4) chromosome (sequenced at the Sanger Center) were absent or divergent in any of the other serovar Enteritidis strains tested. The more recent isolates did not have consistent differences from 60-year-old field isolates, suggesting that no large genomic additions on a whole-gene scale were needed for serovar Enteritidis to become more prevalent in domestic fowl. Cross-hybridization of phage genes on the array with related genes in the examined genomes grouped the serovar Enteritidis isolates into two major lineages. Microarray comparisons of the sequenced serovar Enteritidis PT4 to isolates of the closely related serovars Dublin and Gallinarum (biovars Gallinarum and Pullorum) revealed several genomic areas that distinguished them from serovar Enteritidis and from each other. These differences in gene content could be useful in DNA-based typing and in understanding the different phenotypes of these related serovars.

Characterization of the genome composition of Bartonella koehlerae by microarray comparative genomic hybridization profiling

Bartonella henselae is present in a wide range of wild and domestic feline hosts and causes cat-scratch disease and bacillary angiomatosis in humans. We have estimated here the gene content of Bartonella koehlerae, a novel species isolated from cats that was recently identified as an agent of human endocarditis. The investigation was accomplished by comparative genomic hybridization (CGH) to a microarray constructed from the sequenced 1.93-Mb genome of B. henselae. Control hybridizations of labeled DNA from the human pathogen Bartonella quintana with a reduced genome of 1.58 Mb were performed to evaluate the accuracy of the array for genes with known levels of sequence divergence. Genome size estimates of B. koehlerae by pulsed-field gel electrophoresis matched that calculated by the CGH, indicating a genome of 1.7 to 1.8 Mb with few unique genes. As in B. quintana, sequences in the prophage and the genomic islands were reported absent in B. koehlerae. In addition, sequence variability was recorded in the chromosome II-like region, where B. koehlerae showed an intermediate retention pattern of both coding and noncoding sequences. Although most of the genes missing in B. koehlerae are also absent from B. quintana, its phylogenetic placement near B. henselae suggests independent deletion events, indicating that host specificity is not solely attributed to genes in the genomic islands. Rather, the results underscore the instability of the genomic islands even within bacterial populations adapted to the same host-vector system, as in the case of B. henselae and B. koehlerae.

Design, validation, and application of a seven-strain Staphylococcus aureus PCR product microarray for comparative genomics

Bacterial comparative genomics has been revolutionized by microarrays, but the power of any microarray is dependent on the number and diversity of gene reporters it contains. Staphylococcus aureus is an important human pathogen causing a wide range of invasive and toxin-mediated diseases, and more than 20% of the genome of any isolate consists of variable genes. Seven whole-genome sequences of S. aureus are available, and we exploited this rare opportunity to design, build, and validate a comprehensive, nonredundant PCR product microarray carrying reporters that represent every predicted open reading frame (3,623 probes). Such a comprehensive microarray necessitated a novel design strategy. Validation with the seven sequenced strains showed correct identification of 93.9% of genes present or absent/divergent but was dependent on the method of analysis chosen. Microarray data were highly reproducible, reducing the need for many replicate slides. Interpretation of microarray data was enhanced by focusing on the major areas of variation--the presence or absence of mobile genetic elements (MGEs). We compiled "composite genomes" of every individual MGE and visualized their distribution. This allowed the sensitive discrimination of related isolates, including the first clear description of how isolates of the same clone of epidemic methicillin-resistant S. aureus differ substantially in their carriage of MGEs. These MGEs carry virulence and resistance genes, suggesting differences in pathogenic potential. The novel methods of design and interpretation of data generated from this microarray will enable further studies of S. aureus evolution, epidemiology, and pathogenesis.

Detection of Vi-negative Salmonella enterica serovar typhi in the peripheral blood of patients with typhoid fever in the Faisalabad region of Pakistan

The synthesis and transportation proteins of the Vi capsular polysaccharide of Salmonella enterica serovar Typhi (serovar Typhi) are encoded by the viaB operon, which resides on a 134-kb pathogenicity island known as SPI-7. In recent years, Vi-negative strains of serovar Typhi have been reported in regions where typhoid fever is endemic. However, because Vi negativity can arise during in vitro passage, the clinical significance of Vi-negative serovar Typhi is not clear. To investigate the loss of Vi expression at the genetic level, 60 stored strains of serovar Typhi from the Faisalabad region of Pakistan were analyzed by PCR for the presence of SPI-7 and two genes essential for Vi production: tviA and tviB. Nine of the sixty strains analyzed (15%) tested negative for both tviA and tviB; only two of these strains lacked SPI-7. In order to investigate whether this phenomenon occurred in vivo, blood samples from patients with the clinical symptoms of typhoid fever were also investigated. Of 48 blood samples tested, 42 tested positive by fliC PCR for serovar Typhi; 4 of these were negative for tviA and tviB. Three of these samples tested positive for SPI-7. These results demonstrate that viaB-negative, SPI-7-positive serovar Typhi is naturally occurring and can be detected by PCR in the peripheral blood of typhoid patients in this region. The method described here can be used to monitor the incidence of Vi-negative serovar Typhi in regions where the Vi vaccine is used.

Pilot study to evaluate microarray hybridization as a tool for Salmonella enterica serovar Typhimurium strain differentiation

In developed countries, Salmonella enterica subspecies 1 serovars Enteritidis and Typhimurium range among the most common causes of bacterial food-borne infections. The surveillance and typing of epidemic Salmonella strains are important tools in epidemiology. Usually, Salmonella enterica subspecies 1 serovars are differentiated by serotyping for diagnostic purposes. Further differentiation is done by phage typing as well as molecular typing techniques. Here we have designed and evaluated a prototype DNA microarray as a tool for serovar Typhimurium strain differentiation. It harbors 83 serovar Typhimurium probes obtained by differential subtractive hybridization and from the public database. The microarray yielded reproducible hybridization patterns in repeated hybridizations with chromosomal DNA of the same strain and could differentiate five serovar Typhimurium reference strains (DT204, DT104, DT208, DT36, and LT2). Furthermore, the microarray identified two distinct groups among 13 serovar Typhimurium DT104 strains. This correlated with observations from pulsed-field gel electrophoresis analysis. Twenty-three further serovar Typhimurium strains were analyzed to explore future directions for optimization of the simple 83-probe DNA microarray. The data presented here demonstrate that DNA microarrays harboring small numbers of selected probes are promising tools for serovar Typhimurium strain typing.

Cell size and nucleoid organization of engineered Escherichia coli cells with a reduced genome

The minimization of a genome is necessary to identify experimentally the minimal gene set that contains only those genes that are essential and sufficient to sustain a functioning cell. Recent developments in genetic techniques have made it possible to generate bacteria with a markedly reduced genome. We developed a simple system for formation of markerless chromosomal deletions, and constructed and characterized a series of large-scale chromosomal deletion mutants of Escherichia coli that lack between 2.4 and 29.7% of the parental chromosome. Combining deletion mutations changes cell length and width, and the mutant cells with larger deletions were even longer and wider than the parental cells. The nucleoid organization of the mutants is also changed: the nucleoids occur as multiple small nucleoids and are localized peripherally near the envelope. Inhibition of translation causes them to condense into one or two packed nucleoids, suggesting that the coupling of transcription and translation of membrane proteins peripherally localizes chromosomes. Because these phenotypes are similar to those of spherical cells, those may be a consequence of the morphological change. Based on the nucleoid localization observed with these mutants, we discuss the cellular nucleoid dynamics.

Correlation of phenotype with the genotype of egg-contaminating Salmonella enterica serovar Enteritidis

The genotype of Salmonella enterica serovar Enteritidis was correlated with the phenotype using DNA-DNA microarray hybridization, ribotyping, and Phenotype MicroArray analysis to compare three strains that differed in colony morphology and phage type. No DNA hybridization differences were found between two phage type 13A (PT13A) strains that varied in biofilm formation; however, the ribotype patterns were different. Both PT13A strains had DNA sequences similar to that of bacteriophage Fels2, whereas the PT4 genome to which they were compared, as well as a PT4 field isolate, had a DNA sequence with some similarity to the bacteriophage ST64b sequence. Phenotype MicroArray analysis indicated that the two PT13A strains and the PT4 field isolate had similar respiratory activity profiles at 37 degrees C. However, the wild-type S. enterica serovar Enteritidis PT13A strain grew significantly better in 20% more of the 1,920 conditions tested when it was assayed at 25 degrees C than the biofilm-forming PT13A strain grew. Statistical analysis of the respiratory activity suggested that S. enterica serovar Enteritidis PT4 had a temperature-influenced dimorphic metabolism which at 25 degrees C somewhat resembled the profile of the biofilm-forming PT13A strain and that at 37 degrees C the metabolism was nearly identical to that of the wild-type PT13A strain. Although it is possible that lysogenic bacteriophage alter the balance of phage types on a farm either by lytic competition or by altering the metabolic processes of the host cell in subtle ways, the different physiologies of the S. enterica serovar Enteritidis strains correlated most closely with minor, rather than major, genomic changes. These results strongly suggest that the pandemic of egg-associated human salmonellosis that came into prominence in the 1980s is primarily an example of bacterial adaptive radiation that affects the safety of the food supply.

Generation of subspecies level-specific microbial diagnostic microarrays using genes amplified from subtractive suppression hybridization as microarray probes

The generation of microarray probes with specificity below the species level is an ongoing challenge, not least because the high-throughput detection of microorganisms would be an efficient means of identifying environmentally relevant microbes. Here, we describe how suppression subtractive hybridization (SSH) can be applied to the production of microarray probes that are useful for microbial differentiation at the subspecies level. SSH was used to initially isolate unique genomic sequences of nine Salmonella strains, and these were validated in quadruplicate by microarray analysis. The results obtained indicate that a large group of genes subtracted by SSH could serve together, as one probe, for detecting a microbial subspecies. Similarly, the whole microbial genome (not subjected to SSH) can be used as a species-specific probe. The detailed methods described herein could be used and adapted for the estimation of any cultivable bacteria from different environments.

Chips and SNPs, bugs and thugs: a molecular sleuthing perspective

Recent events both here and abroad have focused attention on the need for ensuring a safe and secure food supply. Although much has been written about the potential of particular select agents in bioterrorism, we must consider seriously the more mundane pathogens, especially those that have been implicated previously in foodborne outbreaks of human disease, as possible agents of bioterrorism. Given their evolutionary history, the enteric pathogens are more diverse than agents such as Bacillus anthracis, Francisella tularensis, or Yersinia pestis. This greater diversity, however, is a double-edged sword; although diversity affords the opportunity for unequivocal identification of an organism without the need for whole-genome sequencing, the same diversity can confound definitive forensic identification if boundaries are not well defined. Here, we discuss molecular approaches used for the identification of Salmonella enterica, Escherichia coli, and Shigella spp. and viral pathogens and discuss the utility of these approaches to the field of microbial molecular forensics.

Genetic and Molecular Analysis of GogB, a Phage-encoded Type III-secreted Substrate in Salmonella enterica Serovar Typhimurium with Autonomous Expression from its Associated Phage

Salmonella enterica serovar Typhimurium is lysogenized by several temperate bacteriophages that encode lysogenic conversion genes, which can act as virulence factors during infection and contribute to the genetic diversity and pathogenic potential of the lysogen. We have investigated the temperate bacteriophage called Gifsy-1 in S.enterica serovar Typhimurium and show here that the product of the gogB gene encoded within this phage shares similarity with proteins from other Gram-negative pathogens. The amino-terminal portion of GogB shares similarity with leucine-rich repeat-containing virulence-associated proteins from other Gram-negative pathogens, whereas the carboxyl-terminal portion of GogB shares similarity with uncharacterized proteins in other pathogens. We show that GogB is secreted by both type III secretion systems encoded in Salmonella Pathogenicity Island-1 (SPI-1) and SPI-2 but translocation into host cells is a SPI-2-mediated process. Once translocated, GogB localizes to the cytoplasm of infected host cells. The genetic regulation of gogB in Salmonella is influenced by the transcriptional activator, SsrB, under SPI-2-inducing conditions, but the modular nature of the gogB gene allows for autonomous expression and type III secretion following horizontal gene transfer into a heterologous pathogen. These data define the first autonomously expressed lysogenic conversion gene within Gifsy-1 that acts as a modular and promiscuous type III-secreted substrate of the infection process.