Comparative phylogenomics of pathogenic bacteria by microarray analysis

Microsatellites in Pursuit of Microbial Genome Evolution

Microsatellites or short sequence repeats are widespread genetic markers which are hypermutable 1-6 bp long short nucleotide motifs. Significantly, their applications in genetics are extensive due to their ceaseless mutational degree, widespread length variations and hypermutability skills. These features make them useful in determining the driving forces of evolution by using powerful molecular techniques. Consequently, revealing important questions, for example, what is the significance of these abundant sequences in DNA, what are their roles in genomic evolution? The answers of these important questions are hidden in the ways these short motifs contributed in altering the microbial genomes since the origin of life. Even though their size ranges from 1 -to- 6 bases, these repeats are becoming one of the most popular genetic probes in determining their associations and phylogenetic relationships in closely related genomes. Currently, they have been widely used in molecular genetics, biotechnology and evolutionary biology. However, due to limited knowledge; there is a significant gap in research and lack of information concerning hypermutational mechanisms. These mechanisms play a key role in microsatellite loci point mutations and phase variations. This review will extend the understandings of impacts and contributions of microsatellite in genomic evolution and their universal applications in microbiology.

Role of environmental survival in transmission of Campylobacter jejuni

Campylobacter species are the most common cause of bacterial gastroenteritis, with C. jejuni responsible for the majority of these cases. Although it is clear that livestock, and particularly poultry, are the most common source, it is likely that the natural environment (soil and water) plays a key role in transmission, either directly to humans or indirectly via farm animals. It has been shown using multilocus sequence typing that some clonal complexes (such as ST-45) are more frequently isolated from environmental sources such as water, suggesting that strains vary in their ability to survive in the environment. Although C. jejuni are fastidious microaerophiles generally unable to grow in atmospheric levels of oxygen, C. jejuni can adapt to survival in the environment, exhibiting aerotolerance and starvation survival. Biofilm formation, the viable but nonculturable state, and interactions with other microorganisms can all contribute to survival outside the host. By exploiting high-throughput technologies such as genome sequencing and RNA Seq, we are well placed to decipher the mechanisms underlying the variations in survival between strains in environments such as soil and water and to better understand the role of environmental persistence in the transmission of C. jejuni directly or indirectly to humans.

A beginner's guide to phylogenetics

Metagenomics and the development of high throughput next generation sequencing capabilities have forced significant development in the field of phylogenetics: the study of the evolutionary relatedness of the planet's inhabitants. Herein, I review the major tree-building strategies, challenges and opportunities which exist in this rapidly expanding field of evolutionary biology.

Molecular diagnostics: the changing culture of medical microbiology

Diagnostic molecular biology is arguably the fastest growing area in current laboratory-based medicine. Growth of the so called 'omics' technologies has, over the last decade, led to a gradual migration away from the 'one test, one pathogen' paradigm, toward multiplex approaches to infectious disease diagnosis, which have led to significant improvements in clinical diagnostics and ultimately improved patient care.

Closely related Campylobacter jejuni strains from different sources reveal a generalist rather than a specialist lifestyle

Background

Campylobacter jejuni and Campylobacter coli are human intestinal pathogens of global importance. Zoonotic transmission from livestock animals or animal-derived food is the likely cause for most of these infections. However, little is known about their general and host-specific mechanisms of colonization, or virulence and pathogenicity factors. In certain hosts, Campylobacter species colonize persistently and do not cause disease, while they cause acute intestinal disease in humans.

Results

Here, we investigate putative host-specificity using phenotypic characterization and genome-wide analysis of genetically closely related C. jejuni strains from different sources. A collection of 473 fresh Campylobacter isolates from Germany was assembled between 2006 and 2010 and characterized using MLST. A subset of closely related C. jejuni strains of the highly prevalent sequence type ST-21 was selected from different hosts and isolation sources. PCR typing of strain-variable genes provided evidence that some genes differed between these strains. Furthermore, phenotypic variation of these strains was tested using the following criteria: metabolic variation, protein expression patterns, and eukaryotic cell interaction. The results demonstrated remarkable phenotypic diversity within the ST-21 group, which however did not correlate with isolation source. Whole genome sequencing was performed for five ST-21 strains from chicken, human, bovine, and food sources, in order to gain insight into ST-21 genome diversity. The comparisons showed extensive genomic diversity, primarily due to recombination and gain of phage-related genes. By contrast, no genomic features associated with isolation source or host were identified.

Conclusions

The genome information and phenotypic data obtained in vitro and in a chicken infection model provided little evidence of fixed adaptation to a specific host. Instead, the dominant C. jejuni ST-21 appeared to be characterized by phenotypic flexibility and high genetic microdiversity, revealing properties of a generalist. High genetic flexibility might allow generalist variants of C. jejuni to reversibly express diverse fitness factors in changing environments.

Molecular basis of virulence in Staphylococcus aureus mastitis

BACKGROUND

S. aureus is one of the main pathogens involved in ruminant mastitis worldwide. The severity of staphylococcal infection is highly variable, ranging from subclinical to gangrenous mastitis. This work represents an in-depth characterization of S. aureus mastitis isolates to identify bacterial factors involved in severity of mastitis infection.

METHODOLOGY/PRINCIPAL FINDINGS

We employed genomic, transcriptomic and proteomic approaches to comprehensively compare two clonally related S. aureus strains that reproducibly induce severe (strain O11) and milder (strain O46) mastitis in ewes. Variation in the content of mobile genetic elements, iron acquisition and metabolism, transcriptional regulation and exoprotein production was observed. In particular, O11 produced relatively high levels of exoproteins, including toxins and proteases known to be important in virulence. A characteristic we observed in other S. aureus strains isolated from clinical mastitis cases.

CONCLUSIONS/SIGNIFICANCE

Our data are consistent with a dose-dependant role of some staphylococcal factors in the hypervirulence of strains isolated from severe mastitis. Mobile genetic elements, transcriptional regulators, exoproteins and iron acquisition pathways constitute good targets for further research to define the underlying mechanisms of mastitis severity.

Array CGH phylogeny: how accurate are comparative genomic hybridization-based trees?

Background

Array-based Comparative Genomic Hybridization (CGH) data have been used to infer phylogenetic relationships. However, the reliability of array CGH analysis to determine evolutionary relationships has not been well established. In most CGH work, all species and strains are compared to a single reference species, whose genome was used to design the array. In the accompanying work, we critically evaluated CGH-based phylogeny using simulated competitive hybridization data. This work showed that a limited number of conditions, principally the tree topology and placement of the reference taxon in the tree, had a strong effect on the ability to recover the correct tree topology. Here, we add to our simulation study by testing the use of CGH as a phylogenetic tool with experimental CGH data from competitive hybridizations between N. crassa and other Neurospora species. In the discussion, we add to our empirical study of Neurospora by reanalyzing of data from a previous CGH phylogenetic analysis of the yeast sensu stricto complex.

Results

Array ratio data for Neurospora and related species were normalized with loess, robust spline, and linear ratio based methods, and then used to construct Neighbor-Joining and parsimony trees. These trees were compared to published phylogenetic analyses for Neurospora based on multilocus sequence analysis (MLSA). For the Neurospora dataset, the best combination of methods resulted in recovery of the MLSA tree topology less than half the time. Our reanalysis of a yeast dataset found that trees identical to established phylogeny were recovered only by pruning taxa - including the reference taxon - from the analysis.

Conclusion

Our results indicate that CGH data can be problematic for phylogenetic analysis. Success fluctuates based on the methods utilized to construct the tree and the taxa included. Selective pruning of the taxa improves the results - an impractical approach for normal phylogenetic analysis. From the more successful methods we make suggestions on the normalization and post-normalization methods that work best in estimating genetic distance between taxa.

Array Comparative Genomic Hybridizations: assessing the ability to recapture evolutionary relationships using an in silico approach

Background

Comparative Genomic Hybridization (CGH) with DNA microarrays has many biological applications including surveys of copy number changes in tumorogenesis, species detection and identification, and functional genomics studies among related organisms. Array CGH has also been used to infer phylogenetic relatedness among species or strains. Although the use of the entire genome can be seen as a considerable advantage for use in phylogenetic analysis, few such studies have questioned the reliability of array CGH to correctly determine evolutionary relationships. A potential flaw in this application lies in the fact that all comparisons are made to a single reference species. This situation differs from traditional DNA sequence, distance-based phylogenetic analyses where all possible pairwise comparisons are made for the isolates in question. By simulating array data based on the Neurospora crassa genome, we address this potential flaw and other questions regarding array CGH phylogeny.

Results

Our simulation data indicates that having a single reference can, in some cases, be a serious limitation when using this technique. Additionally, the tree building process with a single reference is sensitive to many factors including tree topology, choice of tree reconstruction method, and the distance metric used.

Conclusions

Without prior knowledge of the topology and placement of the reference taxon in the topology, the outcome is likely to be wrong and the error undetected. Given these limitations, using CGH to reveal phylogeny based on sequence divergence does not offer a robust alternative to traditional phylogenetic analysis.

A novel genotyping scheme for Vibrio parahaemolyticus with combined use of large variably-presented gene clusters (LVPCs) and variable-number tandem repeats (VNTRs)

A total of 18 variably-presented gene clusters (LVPCs) and nine previously characterized variable-number tandem repeats (VNTRs), and all known virulence markers were screened for their frequency and/or copy number in 251 global strains of Vibrio parahaemolyticus using PCR and gel or capillary electrophoresis. A two-step genotyping approach combining the use of LVPCs and VNTRs was established accordingly. The frequency profiles of LVPCs and virulence markers were primarily used to group the strains into six distinct complexes with different potential pathogenicity natures. The strains from each of these complexes were further analyzed with VNTRs to give a much more detailed discrimination of the strains. A genetic fingerprint-like database of a large collection of strains established with this two-stage approach would be very useful for identification, genotyping, origin tracing, and risk estimation of V. parahaemolyticus.

Genomic variations define divergence of water/wildlife-associated Campylobacter jejuni niche specialists from common clonal complexes

Although the major food-borne pathogen Campylobacter jejuni has been isolated from diverse animal, human and environmental sources, our knowledge of genomic diversity in C. jejuni is based exclusively on human or human food-chain-associated isolates. Studies employing multilocus sequence typing have indicated that some clonal complexes are more commonly associated with particular sources. Using comparative genomic hybridization on a collection of 80 isolates representing diverse sources and clonal complexes, we identified a separate clade comprising a group of water/wildlife isolates of C. jejuni with multilocus sequence types uncharacteristic of human food-chain-associated isolates. By genome sequencing one representative of this diverse group (C. jejuni 1336), and a representative of the bank-vole niche specialist ST-3704 (C. jejuni 414), we identified deletions of genomic regions normally carried by human food-chain-associated C. jejuni. Several of the deleted regions included genes implicated in chicken colonization or in virulence. Novel genomic insertions contributing to the accessory genomes of strains 1336 and 414 were identified. Comparative analysis using PCR assays indicated that novel regions were common but not ubiquitous among the water/wildlife group of isolates, indicating further genomic diversity among this group, whereas all ST-3704 isolates carried the same novel accessory regions. While strain 1336 was able to colonize chicks, strain 414 was not, suggesting that regions specifically absent from the genome of strain 414 may play an important role in this common route of Campylobacter infection of humans. We suggest that the genomic divergence observed constitutes evidence of adaptation leading to niche specialization.

High-throughput microarray technology in diagnostics of enterobacteria based on genome-wide probe selection and regression analysis

Background

The Enterobacteriaceae comprise a large number of clinically relevant species with several individual subspecies. Overlapping virulence-associated gene pools and the high overall genome plasticity often interferes with correct enterobacterial strain typing and risk assessment. Array technology offers a fast, reproducible and standardisable means for bacterial typing and thus provides many advantages for bacterial diagnostics, risk assessment and surveillance. The development of highly discriminative broad-range microbial diagnostic microarrays remains a challenge, because of marked genome plasticity of many bacterial pathogens.

Results

We developed a DNA microarray for strain typing and detection of major antimicrobial resistance genes of clinically relevant enterobacteria. For this purpose, we applied a global genome-wide probe selection strategy on 32 available complete enterobacterial genomes combined with a regression model for pathogen classification. The discriminative power of the probe set was further tested in silico on 15 additional complete enterobacterial genome sequences. DNA microarrays based on the selected probes were used to type 92 clinical enterobacterial isolates. Phenotypic tests confirmed the array-based typing results and corroborate that the selected probes allowed correct typing and prediction of major antibiotic resistances of clinically relevant Enterobacteriaceae, including the subspecies level, e.g. the reliable distinction of different E. coli pathotypes.

Conclusions

Our results demonstrate that the global probe selection approach based on longest common factor statistics as well as the design of a DNA microarray with a restricted set of discriminative probes enables robust discrimination of different enterobacterial variants and represents a proof of concept that can be adopted for diagnostics of a wide range of microbial pathogens. Our approach circumvents misclassifications arising from the application of virulence markers, which are highly affected by horizontal gene transfer. Moreover, a broad range of pathogens have been covered by an efficient probe set size enabling the design of high-throughput diagnostics.

Evaluation of one- and two-color gene expression arrays for microbial comparative genome hybridization analyses in routine applications

DNA microarray technology has already revolutionized basic research in infectious diseases, and whole-genome sequencing efforts have allowed for the fabrication of tailor-made spotted microarrays for an increasing number of bacterial pathogens. However, the application of microarrays in diagnostic microbiology is currently hampered by the high costs associated with microarray experiments and the specialized equipment needed. Here, we show that a thorough bioinformatic postprocessing of the microarray design to reduce the amount of unspecific noise also allows the reliable use of spotted gene expression microarrays for gene content analyses. We further demonstrate that the use of only single-color labeling to halve the costs for dye-labeled nucleotides results in only a moderate decrease in overall specificity and sensitivity. Therefore, gene expression microarrays using only single-color labeling can also reliably be used for gene content analyses, thus reducing the costs for potential routine applications such as genome-based pathogen detection or strain typing.

Genome dynamics of Bartonella grahamii in micro-populations of woodland rodents

Background

Rodents represent a high-risk reservoir for the emergence of new human pathogens. The recent completion of the 2.3 Mb genome of Bartonella grahamii, one of the most prevalent blood-borne bacteria in wild rodents, revealed a higher abundance of genes for host-cell interaction systems than in the genomes of closely related human pathogens. The sequence variability within the global B. grahamii population was recently investigated by multi locus sequence typing, but no study on the variability of putative host-cell interaction systems has been performed.

Results

To study the population dynamics of B. grahamii, we analyzed the genomic diversity on a whole-genome scale of 27 B. grahamii strains isolated from four different species of wild rodents in three geographic locations separated by less than 30 km. Even using highly variable spacer regions, only 3 sequence types were identified. This low sequence diversity contrasted with a high variability in genome content. Microarray comparative genome hybridizations identified genes for outer surface proteins, including a repeated region containing the fha gene for filamentous hemaggluttinin and a plasmid that encodes a type IV secretion system, as the most variable. The estimated generation times in liquid culture medium for a subset of strains ranged from 5 to 22 hours, but did not correlate with sequence type or presence/absence patterns of the fha gene or the plasmid.

Conclusion

Our study has revealed a geographic microstructure of B. grahamii in wild rodents. Despite near-identity in nucleotide sequence, major differences were observed in gene presence/absence patterns that did not segregate with host species. This suggests that genetically similar strains can infect a range of different hosts.

Changes in the genomic content of circulating Bordetella pertussis strains isolated from the Netherlands, Sweden, Japan and Australia: adaptive evolution or drift?

Background

Bordetella pertussis is the causative agent of human whooping cough (pertussis) and is particularly severe in infants. Despite worldwide vaccinations, whooping cough remains a public health problem. A significant increase in the incidence of whooping cough has been observed in many countries since the 1990s. Several reasons for the re-emergence of this highly contagious disease have been suggested. A particularly intriguing possibility is based on evidence indicating that pathogen adaptation may play a role in this process. In an attempt to gain insight into the genomic make-up of B. pertussis over the last 60 years, we used an oligonucleotide DNA microarray to compare the genomic contents of a collection of 171 strains of B. pertussis isolates from different countries.

Results

The CGH microarray analysis estimated the core genome of B. pertussis, to consist of 3,281 CDSs that are conserved among all B. pertussis strains, and represent 84.8% of all CDSs found in the 171 B. pertussis strains. A total of 64 regions of difference consisting of one or more contiguous CDSs were identified among the variable genes. CGH data also revealed that the genome size of B. pertussis strains is decreasing progressively over the past 60 years. Phylogenetic analysis of microarray data generated a minimum spanning tree that depicted the phylogenetic structure of the strains. B. pertussis strains with the same gene content were found in several different countries. However, geographic specificity of the B. pertussis strains was not observed. The gene content was determined to highly correlate with the ptxP-type of the strains.

Conclusions

An overview of genomic contents of a large collection of isolates from different countries allowed us to derive a core genome and a phylogenetic structure of B. pertussis. Our results show that B. pertussis is a dynamic organism that continues to evolve.

Genomic Approaches to the Host Response to Pathogens

The need to better understand host–pathogen interactions has risen with the expansion in genomics and related technologies. This chapter focuses on two aspects of the host response to pathogens where major advances are being made using genomic approaches. The availability of complete genomic sequences of an expanding number of pathogens, the human and mouse genome sequences, and the advent of genome-wide genotyping and gene expression profiling has opened up new avenues of investigation in the field. The genotype of the pathogen plays a major role in the response of the host to infection with more virulent pathogenic strains possessing the capability to interfere with the host immune response. In addition, different individuals in a population can have very different responses to a genetically identical pathogen. Part of the differential response is governed by the underlying genetic differences between individuals. The advent of genome-wide genotyping using single nucleotide polymorphisms or microsatellite markers is leading to major advances in molecular epidemiology. The future impact of genomic approaches on the development of diagnostics and therapeutics is discussed for infectious diseases. This includes defining the basis of genetic susceptibility to infection and system-wide molecular response to a pathogen.

Comparative genomic analysis of Acidithiobacillus ferrooxidans strains using the A. ferrooxidans ATCC 23270 whole-genome oligonucleotide microarray

Acidithiobacillus ferrooxidans is an important microorganism used in biomining operations for metal recovery. Whole-genomic diversity analysis based on the oligonucleotide microarray was used to analyze the gene content of 12 strains of A. ferrooxidans purified from various mining areas in China. Among the 3100 open reading frames (ORFs) on the slides, 1235 ORFs were absent in at least 1 strain of bacteria and 1385 ORFs were conserved in all strains. The hybridization results showed that these strains were highly diverse from a genomic perspective. The hybridization results of 4 major functional gene categories, namely electron transport, carbon metabolism, extracellular polysaccharides, and detoxification, were analyzed. Based on the hybridization signals obtained, a phylogenetic tree was built to analyze the evolution of the 12 tested strains, which indicated that the geographic distribution was the main factor influencing the strain diversity of these strains. Based on the hybridization signals of genes associated with bioleaching, another phylogenetic tree showed an evolutionary relationship from which the co-relation between the clustering of specific genes and geochemistry could be observed. The results revealed that the main factor was geochemistry, among which the following 6 factors were the most important: pH, Mg, Cu, S, Fe, and Al.

On the origin of prokaryotic species

The notion that all prokaryotes belong to genomically and phenomically cohesive clusters that we might legitimately call "species" is a contentious one. At issue are (1) whether such clusters actually exist; (2) what species definition might most reliably identify them, if they do; and (3) what species concept -- by which is meant a genetic and ecological theory of speciation -- might best explain species existence and rationalize a species definition, if we could agree on one. We review existing theories and some relevant data. We conclude that microbiologists now understand in some detail the various genetic, population, and ecological processes that effect the evolution of prokaryotes. There will be on occasion circumstances under which these, working together, will form groups of related organisms sufficiently like each other that we might all agree to call them "species," but there is no reason that this must always be so. Thus, there is no principled way in which questions about prokaryotic species, such as how many there are, how large their populations are, or how globally they are distributed, can be answered. These questions can, however, be reformulated so that metagenomic methods and thinking will meaningfully address the biological patterns and processes whose understanding is our ultimate target.

Genomics-based molecular epidemiology of Campylobacter jejuni isolates from feedlot cattle and from people in Alberta, Canada

Feedlot cattle in Alberta, Canada, have been identified as reservoirs for Campylobacter jejuni, an important human pathogen. Oligonucleotide DNA microarrays were used as a platform to compare C. jejuni isolates from feedlot cattle and human clinical cases from Alberta. Comparative genomic hybridization (CGH) analysis was performed on 87 isolates (46 bovine, 41 human) obtained within the same geographical regions and time frame. Thirteen CGH clusters were obtained based on overall comparative genomic profile similarity. Nine CGH clusters contained human and cattle isolates, three contained only human isolates, and one contained only cattle isolates. The study isolates clustered regardless of temporal or geographical frameworks. In addition, array genes (n = 1,399) were investigated on a gene-by-gene basis to see if any were unequally distributed between human and cattle sources or between clusters dominated by either human or cattle isolates ("human enriched" versus "cattle enriched"). Using Fisher's exact test with the Westfall and Young correction for these comparisons, a small number of differentially distributed genes were identified. Our findings suggest that feedlot cattle and human C. jejuni strains are very similar and may be endemic within Alberta. Further, the common distribution of human clinical and bovine C. jejuni isolates within the same genetically based clusters suggests that dynamic and important transmission routes between cattle and human populations may exist.

Applications of microarrays in pathogen detection and biodefence

The microarray is a platform with wide-ranging potential in biodefence. Owing to the high level of throughput attainable through miniaturization, microarrays have accelerated the ability to respond in an epidemic or crisis. Extending beyond diagnostics, recent studies have applied microarrays as a research tool towards understanding the etiology and pathogenicity of dangerous pathogens, as well as in vaccine development. The original emphasis was on DNA microarrays, but the range now includes protein, antibody and carbohydrate microarrays, and research groups have exploited this diversity to further extend microarray applications in the area of biodefence. Here, we discuss the impact and contributions of the growing range of microarrays and emphasize the concepts that might shape the future of biodefence research.

Characterization of probiotic Escherichia coli isolates with a novel pan-genome microarray

A high-density microarray has been designed that covers the genomes of 24 Escherichia coli and 8 Shigella strains. As a proof-of-principle the genomes of four probiotic E. coli strains were analyzed and their phylogenetic relationship to other E.coli strains investigated.

Background

Microarrays have recently emerged as a novel procedure to evaluate the genetic content of bacterial species. So far, microarrays have mostly covered single or few strains from the same species. However, with cheaper high-throughput sequencing techniques emerging, multiple strains of the same species are rapidly becoming available, allowing for the definition and characterization of a whole species as a population of genomes - the 'pan-genome'.

Results

Using 32 Escherichia coli and Shigella genome sequences we estimate the pan- and core genome of the species. We designed a high-density microarray in order to provide a tool for characterization of the E. coli pan-genome. Technical performance of this pan-genome microarray based on control strain samples (E. coli K-12 and O157:H7) demonstrated a high sensitivity and relatively low false positive rate. A single-channel analysis approach is robust while allowing the possibility for deriving presence/absence predictions for any gene included on our pan-genome microarray. Moreover, the array was highly sufficient to investigate the gene content of non-pathogenic isolates, despite the strong bias towards pathogenic E. coli strains that have been sequenced so far.

Conclusion

This high-density microarray provides an excellent tool for characterizing the genetic makeup of unknown E. coli strains and can also deliver insights into phylogenetic relationships. Its design poses a considerably larger challenge and involves different considerations than the design of single strain microarrays. Here, lessons learned and future directions will be discussed in order to optimize design of microarrays targeting entire pan-genomes.

Comparative genomic characterization of citrus-associated Xylella fastidiosa strains

Background

The xylem-inhabiting bacterium Xylella fastidiosa (Xf) is the causal agent of Pierce's disease (PD) in vineyards and citrus variegated chlorosis (CVC) in orange trees. Both of these economically-devastating diseases are caused by distinct strains of this complex group of microorganisms, which has motivated researchers to conduct extensive genomic sequencing projects with Xf strains. This sequence information, along with other molecular tools, have been used to estimate the evolutionary history of the group and provide clues to understand the capacity of Xf to infect different hosts, causing a variety of symptoms. Nonetheless, although significant amounts of information have been generated from Xf strains, a large proportion of these efforts has concentrated on the study of North American strains, limiting our understanding about the genomic composition of South American strains – which is particularly important for CVC-associated strains.

Results

This paper describes the first genome-wide comparison among South American Xf strains, involving 6 distinct citrus-associated bacteria. Comparative analyses performed through a microarray-based approach allowed identification and characterization of large mobile genetic elements that seem to be exclusive to South American strains. Moreover, a large-scale sequencing effort, based on Suppressive Subtraction Hybridization (SSH), identified 290 new ORFs, distributed in 135 Groups of Orthologous Elements, throughout the genomes of these bacteria.

Conclusion

Results from microarray-based comparisons provide further evidence concerning activity of horizontally transferred elements, reinforcing their importance as major mediators in the evolution of Xf. Moreover, the microarray-based genomic profiles showed similarity between Xf strains 9a5c and Fb7, which is unexpected, given the geographical and chronological differences associated with the isolation of these microorganisms. The newly identified ORFs, obtained by SSH, represent an approximately 10% increase in our current knowledge of the South American Xf gene pool and include new putative virulence factors, as well as novel potential markers for strain identification. Surprisingly, this list of novel elements include sequences previously believed to be unique to North American strains, pointing to the necessity of revising the list of specific markers that may be used for identification of distinct Xf strains.

Bacterial pathogenomics

Genomes from all of the crucial bacterial pathogens of humans, plants and animals have now been sequenced, as have genomes from many of the important commensal, symbiotic and environmental microorganisms. Analysis of these sequences has revealed the forces that shape pathogen evolution and has brought to light unexpected aspects of pathogen biology. The finding that horizontal gene transfer and genome decay have key roles in the evolution of bacterial pathogens was particularly surprising. It has also become evident that even the definitions for 'pathogen' and 'virulence factor' need to be re-evaluated.

Comparative genomic hybridizations reveal absence of large Streptomyces coelicolor genomic islands in Streptomyces lividans

Background

The genomes of Streptomyces coelicolor and Streptomyces lividans bear a considerable degree of synteny. While S. coelicolor is the model streptomycete for studying antibiotic synthesis and differentiation, S. lividans is almost exclusively considered as the preferred host, among actinomycetes, for cloning and expression of exogenous DNA. We used whole genome microarrays as a comparative genomics tool for identifying the subtle differences between these two chromosomes.

Results

We identified five large S. coelicolor genomic islands (larger than 25 kb) and 18 smaller islets absent in S. lividans chromosome. Many of these regions show anomalous GC bias and codon usage patterns. Six of them are in close vicinity of tRNA genes while nine are flanked with near perfect repeat sequences indicating that these are probable recent evolutionary acquisitions into S. coelicolor. Embedded within these segments are at least four DNA methylases and two probable methyl-sensing restriction endonucleases. Comparison with S. coelicolor transcriptome and proteome data revealed that some of the missing genes are active during the course of growth and differentiation in S. coelicolor. In particular, a pair of methylmalonyl CoA mutase (mcm) genes involved in polyketide precursor biosynthesis, an acyl-CoA dehydrogenase implicated in timing of actinorhodin synthesis and bldB, a developmentally significant regulator whose mutation causes complete abrogation of antibiotic synthesis belong to this category.

Conclusion

Our findings provide tangible hints for elucidating the genetic basis of important phenotypic differences between these two streptomycetes. Importantly, absence of certain genes in S. lividans identified here could potentially explain the relative ease of DNA transformations and the conditional lack of actinorhodin synthesis in S. lividans.

Use of suppression subtractive hybridisation to extend our knowledge of genome diversity in Campylobacter jejuni

BACKGROUND

Previous studies have sought to identify a link between the distribution of variable genes amongst isolates of Campylobacter jejuni and particular host preferences. The genomic sequence data available currently was obtained using only isolates from human or chicken hosts. In order to identify variable genes present in isolates from alternative host species, five subtractions between C. jejuni isolates from different sources (rabbit, cattle, wild bird) were carried out, designed to assess genomic variability within and between common multilocus sequence type (MLST) clonal complexes (ST-21, ST-42, ST-45 and ST-61).

RESULTS

The vast majority (97%) of the 195 subtracted sequences identified had a best BLASTX match with a Campylobacter protein. However, there was considerable variation within and between the four clonal complexes included in the subtractions. The distributions of eight variable sequences, including four with putative roles in the use of alternative terminal electron acceptors, amongst a panel of C. jejuni isolates representing diverse sources and STs, were determined.

CONCLUSION

There was a clear correlation between clonal complex and the distribution of the metabolic genes. In contrast, there was no evidence to support the hypothesis that the distribution of such genes may be related to host preference. The other variable genes studied were also generally distributed according to MLST type. Thus, we found little evidence for widespread horizontal gene transfer between clonal complexes involving these genes.

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.

Bacillus subtilis genome diversity

Microarray-based comparative genomic hybridization (M-CGH) is a powerful method for rapidly identifying regions of genome diversity among closely related organisms. We used M-CGH to examine the genome diversity of 17 strains belonging to the nonpathogenic species Bacillus subtilis. Our M-CGH results indicate that there is considerable genetic heterogeneity among members of this species; nearly one-third of Bsu168-specific genes exhibited variability, as measured by the microarray hybridization intensities. The variable loci include those encoding proteins involved in antibiotic production, cell wall synthesis, sporulation, and germination. The diversity in these genes may reflect this organism's ability to survive in diverse natural settings.

Differentiation of Escherichia coli pathotypes by oligonucleotide spotted array

To accurately determine the pathotypes of Escherichia coli strains, a comprehensive assessment of each strain that targets multiple genes is required. A new approach to the identification and characterization of E. coli pathotypes was developed by constructing gene-specific probes (70-mers) for not only the virulence genes associated with each E. coli pathotype but also the O157-, CFT073-, and K-12-specific and common genes of each pathotype. Analysis of oligonucleotide probes with reference and clinical isolates of E. coli pathotypes indicated that the array could differentiate the pathotypes on the basis of their virulence and specific gene patterns. Probes targeting common genes of E. coli were present in all the reference and clinical strains. Salmonella enterica subsp. enterica-specific genes and Salmonella core genes were used as negative controls. The entire E. coli pathotype showed reactivity to only 4 of the 81 Salmonella-specific gene probes. Characterization of the genetic and virulence profiles of a single strain by using probes for virulence factors and specific and common genes in the spotted array is an ideal diagnostic tool for determination of E. coli pathotypes and could also have a significant impact on the epidemiological analysis of E. coli infections.

Comparative phylogenomics of Clostridium difficile reveals clade specificity and microevolution of hypervirulent strains

Clostridium difficile is the most frequent cause of nosocomial diarrhea worldwide, and recent reports suggested the emergence of a hypervirulent strain in North America and Europe. In this study, we applied comparative phylogenomics (whole-genome comparisons using DNA microarrays combined with Bayesian phylogenies) to model the phylogeny of C. difficile, including 75 diverse isolates comprising hypervirulent, toxin-variable, and animal strains. The analysis identified four distinct statistically supported clusters comprising a hypervirulent clade, a toxin A(-) B(+) clade, and two clades with human and animal isolates. Genetic differences among clades revealed several genetic islands relating to virulence and niche adaptation, including antibiotic resistance, motility, adhesion, and enteric metabolism. Only 19.7% of genes were shared by all strains, confirming that this enteric species readily undergoes genetic exchange. This study has provided insight into the possible origins of C. difficile and its evolution that may have implications in disease control strategies.

Unique DNA sequences of avian pathogenicEscherichia coliisolates as determined by genomic suppression subtractive hybridization

Genomic suppression subtractive hybridization was carried out between an avian pathogenic Escherichia coli (APEC) and an E. coli isolate from apparently healthy poultry (avian fecal E. coli or AFEC) to understand the genomic differences between APEC and AFEC and to identify potential virulence genes. In an effort to ensure that the genes identified were novel, tester and driver strains were chosen by their similarities in certain known virulence genes. In total, 24 subtractive hybridization fragments specific for APEC were obtained. Eight of the fragments were phage or plasmid related, 13 matched to chromosomal regions and three did not show homology to any known genes. We further determined the association of APEC-specific fragments with virulence by screening a collection of APEC, AFEC and human uropathogenic E. coli (UPEC) using the PCR. Chromosomally located fragments that matched to the ireA iron-responsive element gene, the tia invasion determinant gene, the ibeA invasion gene, genes encoding P pili and a gene encoding aldo/keto reductase were present at significantly higher frequencies in APEC and UPEC than in AFEC, suggesting their contribution to the virulence of extraintestinal pathogenic E. coli. No genes were found to be specific to APEC alone.

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 comparative phylogenomics to study the evolution of Yersinia enterocolitica and to identify genetic differences relating to pathogenicity

Yersinia enterocolitica, an important cause of human gastroenteritis generally caused by the consumption of livestock, has traditionally been categorized into three groups with respect to pathogenicity, i.e., nonpathogenic (biotype 1A), low pathogenicity (biotypes 2 to 5), and highly pathogenic (biotype 1B). However, genetic differences that explain variation in pathogenesis and whether different biotypes are associated with specific nonhuman hosts are largely unknown. In this study, we applied comparative phylogenomics (whole-genome comparisons of microbes with DNA microarrays combined with Bayesian phylogenies) to investigate a diverse collection of 94 strains of Y. enterocolitica consisting of 35 human, 35 pig, 15 sheep, and 9 cattle isolates from nonpathogenic, low-pathogenicity, and highly pathogenic biotypes. Analysis confirmed three distinct statistically supported clusters composed of a nonpathogenic clade, a low-pathogenicity clade, and a highly pathogenic clade. Genetic differences revealed 125 predicted coding sequences (CDSs) present in all highly pathogenic strains but absent from the other clades. These included several previously uncharacterized CDSs that may encode novel virulence determinants including a hemolysin, a metalloprotease, and a type III secretion effector protein. Additionally, 27 CDSs were identified which were present in all 47 low-pathogenicity strains and Y. enterocolitica 8081 but absent from all nonpathogenic 1A isolates. Analysis of the core gene set for Y. enterocolitica revealed that 20.8% of the genes were shared by all of the strains, confirming this species as highly heterogeneous, adding to the case for the existence of three subspecies of Y. enterocolitica. Further analysis revealed that Y. enterocolitica does not cluster according to source (host).