Application of DNA microarrays to study the evolutionary genomics of Yersinia pestis and Yersinia pseudotuberculosis

A frameshift in Yersinia pestis rcsD alters canonical Rcs signalling to preserve flea-mammal plague transmission cycles

Multiple genetic changes in the enteric pathogen Yersinia pseudotuberculosis have driven the emergence of Yesinia pestis, the arthropod-borne, etiological agent of plague. These include developing the capacity for biofilm-dependent blockage of the flea foregut to enable transmission by flea bite. Previously, we showed that pseudogenization of rcsA, encoding a component of the Rcs signalling pathway, is an important evolutionary step facilitating Y. pestis flea-borne transmission. Additionally, rcsD, another important gene in the Rcs system, harbours a frameshift mutation. Here, we demonstrated that this rcsD mutation resulted in production of a small protein composing the C-terminal RcsD histidine-phosphotransferase domain (designated RcsD-Hpt) and full-length RcsD. Genetic analysis revealed that the rcsD frameshift mutation followed the emergence of rcsA pseudogenization. It further altered the canonical Rcs phosphorylation signal cascade, fine-tuning biofilm production to be conducive with retention of the pgm locus in modern lineages of Y. pestis. Taken together, our findings suggest that a frameshift mutation in rcsD is an important evolutionary step that fine-tuned biofilm production to ensure perpetuation of flea-mammal plague transmission cycles.

Common architecture of Tc toxins from human and insect pathogenic bacteria

Tc toxins use a syringe-like mechanism to penetrate the membrane and translocate toxic enzymes into the host cytosol. They are composed of three components: TcA, TcB, and TcC. Low-resolution structures of TcAs from different bacteria suggest a considerable difference in their architecture and possibly in their mechanism of action. Here, we present high-resolution structures of five TcAs from insect and human pathogens, which show a similar overall composition and domain organization. Essential structural features, including a trefoil protein knot, are present in all TcAs, suggesting a common mechanism of action. All TcAs form functional pores and can be combined with TcB-TcC subunits from other species to form active chimeric holotoxins. We identified a conserved ionic pair that stabilizes the shell, likely operating as a strong latch that only springs open after destabilization of other regions. Our results provide new insights into the architecture and mechanism of the Tc toxin family.

Pan-genomic analyses identify key Helicobacter pylori pathogenic loci modified by carcinogenic host microenvironments

OBJECTIVE

Helicobacter pylori is the strongest risk factor for gastric cancer; however, the majority of infected individuals do not develop disease. Pathological outcomes are mediated by complex interactions among bacterial, host and environmental constituents, and two dietary factors linked with gastric cancer risk are iron deficiency and high salt. We hypothesised that prolonged adaptation of H. pylori to in vivo carcinogenic microenvironments results in genetic modification important for disease.

DESIGN

Whole genome sequencing of genetically related H. pylori strains that differ in virulence and targeted H. pylori sequencing following prolonged exposure of bacteria to in vitro carcinogenic conditions were performed.

RESULTS

A total of 180 unique single nucleotide polymorphisms (SNPs) were identified among the collective genomes when compared with a reference H. pylori genome. Importantly, common SNPs were identified in isolates harvested from iron-depleted and high salt carcinogenic microenvironments, including an SNP within fur (FurR88H). To investigate the direct role of low iron and/or high salt, H. pylori was continuously cultured in vitro under low iron or high salt conditions to assess fur genetic variation. Exposure to low iron or high salt selected for the FurR88H variant after only 5 days. To extend these results, fur was sequenced in 339 clinical H. pylori strains. Among the isolates examined, 17% (40/232) of strains isolated from patients with premalignant lesions harboured the FurR88H variant, compared with only 6% (6/107) of strains from patients with non-atrophic gastritis alone (p=0.0034).

CONCLUSION

These results indicate that specific genetic variation arises within H. pylori strains during in vivo adaptation to conditions conducive for gastric carcinogenesis.

Secretion Systems and Secreted Proteins in Gram-Negative Entomopathogenic Bacteria: Their Roles in Insect Virulence and Beyond

Many Gram-negative bacteria have evolved insect pathogenic lifestyles. In all cases, the ability to cause disease in insects involves specific bacterial proteins exported either to the surface, the extracellular environment, or the cytoplasm of the host cell. They also have several distinct mechanisms for secreting such proteins. In this review, we summarize the major protein secretion systems and discuss examples of secreted proteins that contribute to the virulence of a variety of Gram-negative entomopathogenic bacteria, including Photorhabdus, Xenorhabdus, Serratia, Yersinia, and Pseudomonas species. We also briefly summarize two classes of exported protein complexes, the PVC-like elements, and the Tc toxin complexes that were first described in entomopathogenic bacteria.

A review of methods for subtyping Yersinia pestis: From phenotypes to whole genome sequencing

Numerous subtyping methods have been applied to Yersinia pestis with varying success. Here, we review the various subtyping methods that have been applied to Y. pestis and their capacity for answering questions regarding the population genetics, phylogeography, and molecular epidemiology of this important human pathogen. Methods are evaluated in terms of expense, difficulty, transferability among laboratories, discriminatory power, usefulness for different study questions, and current applicability in light of the advent of whole genome sequencing.

A Yersinia pestis tat mutant is attenuated in bubonic and small-aerosol pneumonic challenge models of infection but not as attenuated by intranasal challenge

Bacterial proteins destined for the Tat pathway are folded before crossing the inner membrane and are typically identified by an N-terminal signal peptide containing a twin arginine motif. Translocation by the Tat pathway is dependent on the products of genes which encode proteins possessing the binding site of the signal peptide and mediating the actual translocation event. In the fully virulent CO92 strain of Yersinia pestis, the tatA gene was deleted. The mutant was assayed for loss of virulence through various in vitro and in vivo assays. Deletion of the tatA gene resulted in several consequences for the mutant as compared to wild-type. Cell morphology of the mutant bacteria was altered and demonstrated a more elongated form. In addition, while cultures of the mutant strain were able to produce a biofilm, we observed a loss of adhesion of the mutant biofilm structure compared to the biofilm produced by the wild-type strain. Immuno-electron microscopy revealed a partial disruption of the F1 antigen on the surface of the mutant. The virulence of the ΔtatA mutant was assessed in various murine models of plague. The mutant was severely attenuated in the bubonic model with full virulence restored by complementation with the native gene. After small-particle aerosol challenge in a pneumonic model of infection, the mutant was also shown to be attenuated. In contrast, when mice were challenged intranasally with the mutant, very little difference in the LD₅₀ was observed between wild-type and mutant strains. However, an increased time-to-death and delay in bacterial dissemination was observed in mice infected with the ΔtatA mutant as compared to the parent strain. Collectively, these findings demonstrate an essential role for the Tat pathway in the virulence of Y. pestis in bubonic and small-aerosol pneumonic infection but less important role for intranasal challenge.

A mutant with aberrant extracellular LcrV-YscF interactions fails to form pores and translocate Yop effector proteins but retains the ability to trigger Yop secretion in response to host cell contact

The plasmid-encoded type three secretion system (TTSS) of Yersinia spp. is responsible for the delivery of effector proteins into cells of the innate immune system, where these effectors disrupt the target cells' activity. Successful translocation of effectors into mammalian cells requires Yersinia to both insert a translocon into the host cell membrane and sense contact with host cells. To probe the events necessary for translocation, we investigated protein-protein interactions among TTSS components of the needle-translocon complex using a chemical cross-linking-based approach. We detected extracellular protein complexes containing YscF, LcrV, and YopD that were dependent upon needle formation. The formation of these complexes was evaluated in a secretion-competent but translocation-defective mutant, the YscFD28AD46A strain (expressing YscF with the mutations D28A and D46A). We found that one of the YscF and most of the LcrV and YopD cross-linked complexes were nearly absent in this mutant. Furthermore, the YscFD28AD46A strain did not support YopB insertion into mammalian membranes, supporting the idea that the LcrV tip complex is required for YopB insertion and translocon formation. However, the YscFD28AD46A strain did secrete Yops in the presence of host cells, indicating that a translocation-competent tip complex is not required to sense contact with host cells to trigger Yop secretion. In conclusion, in the absence of cross-linkable LcrV-YscF interactions, translocon insertion is abolished, but Yersinia still retains the ability to sense cell contact.

Yersinia pestis versus Yersinia pseudotuberculosis: effects on host macrophages

Yersinia pestis, the causative agent of plague, is proved to be a recently emerged clone from Y. pseudotuberculosis. However, the diseases they cause and their patterns of transmission are very different. People always focus on the genetic changes between Y. pestis and Y. pseudotuberculosis to reveal their pathogenic differences, and little is known about host defence differences to these two Yersinia. In this study, the effects of Y. pestis and Y. pseudotuberculosis on macrophages were analysed. Cell apoptosis showed significant difference after the macrophages were infected by these two strains, and caspase-3 activity also demonstrated a similar tendency. Further, macrophage function activities were evaluated. We found during the early infection of Y. pestis, several basic functions of macrophages, including phagocytosis, secretion of cytokine tumour necrosis factor-α and nitric oxide, macrophage polarity and antigen presenting, were significantly interrupted. In comparison, Y. pseudotuberculosis infection showed lower inhibition on macrophages. Especially, Y. pestis infection might cause macrophage to polarize to M2 macrophages in the early phase, compared with Y. pseudotuberculosis infection, which was different from the common acute infection. These results clearly indicated even in the early stage of infection, different host macrophage defence patterns could help us to understand the obvious virulence differences between Y. pestis and Y. pseudotuberculosis.

Transgenic approaches to western corn rootworm control

The western corn rootworm, Diabrotica virgifera virgifera LeConte (Coleoptera: Chrysomelidae) is a significant corn pest throughout the United States corn belt. Rootworm larvae feed on corn roots causing yield losses and control expenditures that are estimated to exceed US$1 billion annually. Traditional management practices to control rootworms such as chemical insecticides or crop rotation have suffered reduced effectiveness due to the development of physiological and behavioral resistance. Transgenic maize expressing insecticidal proteins are very successful in protecting against rootworm damage and preserving corn yield potential. However, the high rate of grower adoption and early reliance on hybrids expressing a single mode of action and low-dose traits threatens the durability of commercialized transgenic rootworm technology for rootworm control. A summary of current transgenic approaches for rootworm control and the corresponding insect resistance management practices is included. An overview of potential new modes of action based on insecticidal proteins, and especially RNAi targeting mRNA coding for essential insect proteins is provided.

A multi-omic systems approach to elucidating Yersinia virulence mechanisms

The underlying mechanisms that lead to dramatic differences between closely related pathogens are not always readily apparent. For example, the genomes of Yersinia pestis (YP) the causative agent of plague with a high mortality rate and Yersinia pseudotuberculosis (YPT) an enteric pathogen with a modest mortality rate are highly similar with some species specific differences; however the molecular causes of their distinct clinical outcomes remain poorly understood. In this study, a temporal multi-omic analysis of YP and YPT at physiologically relevant temperatures was performed to gain insights into how an acute and highly lethal bacterial pathogen, YP, differs from its less virulent progenitor, YPT. This analysis revealed higher gene and protein expression levels of conserved major virulence factors in YP relative to YPT, including the Yop virulon and the pH6 antigen. This suggests that adaptation in the regulatory architecture, in addition to the presence of unique genetic material, may contribute to the increased pathogenecity of YP relative to YPT. Additionally, global transcriptome and proteome responses of YP and YPT revealed conserved post-transcriptional control of metabolism and the translational machinery including the modulation of glutamate levels in Yersiniae. Finally, the omics data was coupled with a computational network analysis, allowing an efficient prediction of novel Yersinia virulence factors based on gene and protein expression patterns.

Evolution and virulence contributions of the autotransporter proteins YapJ and YapK of Yersinia pestis CO92 and their homologs in Y. pseudotuberculosis IP32953

Yersinia pestis, the causative agent of plague, evolved from the gastrointestinal pathogen Yersinia pseudotuberculosis. Both species have numerous type Va autotransporters, most of which appear to be highly conserved. In Y. pestis CO92, the autotransporter genes yapK and yapJ share a high level of sequence identity. By comparing yapK and yapJ to three homologous genes in Y. pseudotuberculosis IP32953 (YPTB0365, YPTB3285, and YPTB3286), we show that yapK is conserved in Y. pseudotuberculosis, while yapJ is unique to Y. pestis. All of these autotransporters exhibit >96% identity in the C terminus of the protein and identities ranging from 58 to 72% in their N termini. By extending this analysis to include homologous sequences from numerous Y. pestis and Y. pseudotuberculosis strains, we determined that these autotransporters cluster into a YapK (YPTB3285) class and a YapJ (YPTB3286) class. The YPTB3286-like gene of most Y. pestis strains appears to be inactivated, perhaps in favor of maintaining yapJ. Since autotransporters are important for virulence in many bacterial pathogens, including Y. pestis, any change in autotransporter content should be considered for its impact on virulence. Using established mouse models of Y. pestis infection, we demonstrated that despite the high level of sequence identity, yapK is distinct from yapJ in its contribution to disseminated Y. pestis infection. In addition, a mutant lacking both of these genes exhibits an additive attenuation, suggesting nonredundant roles for yapJ and yapK in systemic Y. pestis infection. However, the deletion of the homologous genes in Y. pseudotuberculosis does not seem to impact the virulence of this organism in orogastric or systemic infection models.

Evolutionary insights into scleractinian corals using comparative genomic hybridizations

Background

Coral reefs belong to the most ecologically and economically important ecosystems on our planet. Yet, they are under steady decline worldwide due to rising sea surface temperatures, disease, and pollution. Understanding the molecular impact of these stressors on different coral species is imperative in order to predict how coral populations will respond to this continued disturbance. The use of molecular tools such as microarrays has provided deep insight into the molecular stress response of corals. Here, we have performed comparative genomic hybridizations (CGH) with different coral species to an Acropora palmata microarray platform containing 13,546 cDNA clones in order to identify potentially rapidly evolving genes and to determine the suitability of existing microarray platforms for use in gene expression studies (via heterologous hybridization).

Results

Our results showed that the current microarray platform for A. palmata is able to provide biological relevant information for a wide variety of coral species covering both the complex clade as well the robust clade. Analysis of the fraction of highly diverged genes showed a significantly higher amount of genes without annotation corroborating previous findings that point towards a higher rate of divergence for taxonomically restricted genes. Among the genes with annotation, we found many mitochondrial genes to be highly diverged in M. faveolata when compared to A. palmata, while the majority of nuclear encoded genes maintained an average divergence rate.

Conclusions

The use of present microarray platforms for transcriptional analyses in different coral species will greatly enhance the understanding of the molecular basis of stress and health and highlight evolutionary differences between scleractinian coral species. On a genomic basis, we show that cDNA arrays can be used to identify patterns of divergence. Mitochondrion-encoded genes seem to have diverged faster than nuclear encoded genes in robust corals. Accordingly, this needs to be taken into account when using mitochondrial markers for scleractinian phylogenies.

Cluster analysis of host cytokine responses to biodefense pathogens in a whole blood ex vivo exposure model (WEEM)

Background

Rapid detection and therapeutic intervention for infectious and emerging diseases is a major scientific goal in biodefense and public health. Toward this end, cytokine profiles in human blood were investigated using a human whole blood ex vivo exposure model, called WEEM.

Results

Samples of whole blood from healthy volunteers were incubated with seven pathogens including Yersinia pseudotuberculosis, Yersinia enterocolitica, Bacillus anthracis, and multiple strains of Yersinia pestis, and multiplexed protein expression profiling was conducted on supernatants of these cultures with an antibody array to detect 30 cytokines simultaneously. Levels of 8 cytokines, IL-1α, IL-1β, IL-6, IL-8, IL-10, IP-10, MCP-1 and TNFα, were significantly up-regulated in plasma after bacterial exposures of 4 hours. Statistical clustering was applied to group the pathogens based on the host response protein expression profiles. The nearest phylogenetic neighbors clustered more closely than the more distant pathogens, and all seven pathogens were clearly differentiated from the unexposed control. In addition, the Y. pestis and Yersinia near neighbors were differentiated from the B. anthracis strains.

Conclusions

Cluster analysis, based on host response cytokine profiles, indicates that distinct patterns of immunomodulatory proteins are induced by the different pathogen exposures and these patterns may enable further development into biomarkers for diagnosing pathogen exposure.

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.

Expression during host infection and localization of Yersinia pestis autotransporter proteins

Yersinia pestis CO92 has 12 open reading frames encoding putative conventional autotransporters (yaps), nine of which appear to produce functional proteins. Here, we demonstrate the ability of the Yap proteins to localize to the cell surface of both Escherichia coli and Yersinia pestis and show that a subset of these proteins undergoes processing by bacterial surface omptins to be released into the supernatant. Numerous autotransporters have been implicated in pathogenesis, suggesting a role for the Yaps as virulence factors in Y. pestis. Using the C57BL/6 mouse models of bubonic and pneumonic plague, we determined that all of these genes are transcribed in the lymph nodes during bubonic infection and in the lungs during pneumonic infection, suggesting a role for the Yaps during mammalian infection. In vitro transcription studies did not identify a particular environmental stimulus responsible for transcriptional induction. The primary sequences of the Yaps reveal little similarity to any characterized autotransporters; however, two of the genes are present in operons, suggesting that the proteins encoded in these operons may function together. Further work aims to elucidate the specific functions of the Yaps and clarify the contributions of these proteins to Y. pestis pathogenesis.

Identification and characterisation of a novel adhesin Ifp in Yersinia pseudotuberculosis

Background

In order to identify new virulence determinants in Y. pseudotuberculosis a comparison between its genome and that of Yersinia pestis was undertaken. This reveals dozens of pseudogenes in Y. pestis, which are still putatively functional in Y. pseudotuberculosis and may be important in the enteric lifestyle. One such gene, YPTB1572 in the Y. pseudotuberculosis IP32953 genome sequence, encodes a protein with similarity to invasin, a classic adhesion/invasion protein, and to intimin, the attaching and effacing protein from enteropathogenic (EPEC) and enterohaemorraghic (EHEC) Escherichia coli.

Results

We termed YPTB1572 Ifp (Intimin family protein) and show that it is able to bind directly to human HEp-2 epithelial cells. Cysteine and tryptophan residues in the C-terminal region of intimin that are essential for function in EPEC and EHEC are conserved in Ifp. Protein binding occurred at distinct foci on the HEp-2 cell surface and can be disrupted by mutation of a single cysteine residue at the C-terminus of the protein. Temporal expression analysis using lux reporter constructs revealed that ifp is expressed at late log phase at 37°C in contrast to invasin, suggesting that Ifp is a late stage adhesin. An ifp defined mutant showed a reduction in adhesion to HEp-2 cells and was attenuated in the Galleria mellonella infection model.

Conclusion

A new Y. pseudotuberculosis adhesin has been identified and characterised. This Ifp is a new member in the family of invasin/intimin outer membrane adhesins.

The diversity of insect-bacteria interactions and its applications for disease control

Prokaryotic microorganisms are widespread in all environments on Earth, establishing diverse interactions with many eukaryotic taxa, including insects. These associations may be symbiotic, pathogenic and vectoring. Independently of the type of interaction, each association starts with the adhesion of the microorganism to the host, entry and "invasion" of the host, then progresses to establishment and dissemination within the host, by avoiding host immune responses, and concludes with transmission back to the environment or to a new host. Advances in genomics and genetics have allowed the dissection of these processes and provided important information on the elements driving the shaping of the members of each association. Furthermore, many mechanisms involved in the establishment of the associations have been scrutinised, along with the development of new methods for the management of insect populations.

Targeting of the actin cytoskeleton by insecticidal toxins from Photorhabdus luminescens

Photorhabdus luminescens produces several types of protein toxins, which are essential for participation in a trilateral symbiosis with nematodes and insects. The nematodes, carrying the bacteria, invade insect larvae and release the bacteria, which kill the insects with their toxins. Recently, the molecular mechanisms of the toxin complexes PTC3 and PTC5 have been elucidated. The biologically active components of the toxin complexes are ADP-ribosyltransferases, which modify actin and Rho GTPases, respectively. The actions of the toxins are described and compared with other bacterial protein toxins acting on the cytoskeleton.

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.

Clostridium difficile using DNA microarrays

Clostridium difficile is a pathogen on the move, as evidenced by the rapid transcontinental spread of the so-called hypervirulent 027 strains, followed by the emergence of further PCR ribotypes such as 017, 078 and 106. This provides a rare opportunity to study the evolution of virulence in action. However, to fully exploit this opportunity, robust phylogenetic methods on a diverse set of characterised strains are required to provide a reference evolutionary framework to study C. difficile epidemiology, ecology and virulence. Traditional phylogenetic classification of bacteria to study evolutionary relatedness is based on the characterisation of a limited number of genes, rRNA or signature sequences. However, due to the acquisition of DNA through lateral gene transfer, the differences between closely related bacterial strains can be vast. By contrast, whole genome sequencing comparisons allow all genes to be compared. Nevertheless, whole-scale genome sequencing remains an expensive endeavour and such comparisons are limited to only a handful of strains. DNA microarrays represent an alternative technology for whole genome comparisons enabling a "birds eye view" of all the genes absent or present in a given genome as compared to the reference genome on the microarray.

Accurate, rapid and high-throughput detection of strain-specific polymorphisms in Bacillus anthracis and Yersinia pestis by next-generation sequencing

Background

In the event of biocrimes or infectious disease outbreaks, high-resolution genetic characterization for identifying the agent and attributing it to a specific source can be crucial for an effective response. Until recently, in-depth genetic characterization required expensive and time-consuming Sanger sequencing of a few strains, followed by genotyping of a small number of marker loci in a panel of isolates at or by gel-based approaches such as pulsed field gel electrophoresis, which by necessity ignores most of the genome. Next-generation, massively parallel sequencing (MPS) technology (specifically the Applied Biosystems sequencing by oligonucleotide ligation and detection (SOLiD™) system) is a powerful investigative tool for rapid, cost-effective and parallel microbial whole-genome characterization.

Results

To demonstrate the utility of MPS for whole-genome typing of monomorphic pathogens, four Bacillus anthracis and four Yersinia pestis strains were sequenced in parallel. Reads were aligned to complete reference genomes, and genomic variations were identified. Resequencing of the B. anthracis Ames ancestor strain detected no false-positive single-nucleotide polymorphisms (SNPs), and mapping of reads to the Sterne strain correctly identified 98% of the 133 SNPs that are not clustered or associated with repeats. Three geographically distinct B. anthracis strains from the A branch lineage were found to have between 352 and 471 SNPs each, relative to the Ames genome, and one strain harbored a genomic amplification. Sequencing of four Y. pestis strains from the Orientalis lineage identified between 20 and 54 SNPs per strain relative to the CO92 genome, with the single Bolivian isolate having approximately twice as many SNPs as the three more closely related North American strains. Coverage plotting also revealed a common deletion in two strains and an amplification in the Bolivian strain that appear to be due to insertion element-mediated recombination events. Most private SNPs (that is, a, variant found in only one strain in this set) selected for validation by Sanger sequencing were confirmed, although rare false-positive SNPs were associated with variable nucleotide tandem repeats.

Conclusions

The high-throughput, multiplexing capability, and accuracy of this system make it suitable for rapid whole-genome typing of microbial pathogens during a forensic or epidemiological investigation. By interrogating nearly every base of the genome, rare polymorphisms can be reliably discovered, thus facilitating high-resolution strain tracking and strengthening forensic attribution.

Insights into the infective properties of YpfΦ, the Yersinia pestis filamentous phage

YpfΦ is a filamentous phage that infected Yersinia pestis, the plague bacillus, during its emergence. Using an experimental transduction approach, we show here that this phage has the capacity to infect with variable efficiencies, all three pathogenic Yersinia species as well as Escherichia coli. Like other Inovirus phages, its genetic organization comprises three functional modules necessary for the production of infectious virions. Upon infection, YpfΦ integrates into the chromosomal dif site, but extrachromosomal forms are also frequently observed. Several pieces of evidence suggest that the absence of chromosomal YpfΦ in natural non-Orientalis Y. pestis isolates results from a higher chromosomal excision rate rather than from a defective integration machinery. A resident YpfΦ confers some protection against a superinfection. In contrast to other filamentous phages, the incoming YpfΦ genome inserts itself between two copies of the resident prophage. This analysis thus unravels infective properties specific to YpfΦ.

Delineation and analysis of chromosomal regions specifying Yersinia pestis

Yersinia pestis, the causative agent of plague, has recently diverged from the less virulent enteropathogen Yersinia pseudotuberculosis. Its emergence has been characterized by massive genetic loss and inactivation and limited gene acquisition. The acquired genes include two plasmids, a filamentous phage, and a few chromosomal loci. The aim of this study was to characterize the chromosomal regions acquired by Y. pestis. Following in silico comparative analysis and PCR screening of 98 strains of Y. pseudotuberculosis and Y. pestis, we found that eight chromosomal loci (six regions [R1pe to R6pe] and two coding sequences [CDS1pe and CDS2pe]) specified Y. pestis. Signatures of integration by site specific or homologous recombination were identified for most of them. These acquisitions and the loss of ancestral DNA sequences were concentrated in a chromosomal region opposite to the origin of replication. The specific regions were acquired very early during Y. pestis evolution and were retained during its microevolution, suggesting that they might bring some selective advantages. Only one region (R3pe), predicted to carry a lambdoid prophage, is most likely no longer functional because of mutations. With the exception of R1pe and R2pe, which have the potential to encode a restriction/modification and a sugar transport system, respectively, no functions could be predicted for the other Y. pestis-specific loci. To determine the role of the eight chromosomal loci in the physiology and pathogenicity of the plague bacillus, each of them was individually deleted from the bacterial chromosome. None of the deletants exhibited defects during growth in vitro. Using the Xenopsylla cheopis flea model, all deletants retained the capacity to produce a stable and persistent infection and to block fleas. Similarly, none of the deletants caused any acute flea toxicity. In the mouse model of infection, all deletants were fully virulent upon subcutaneous or aerosol infections. Therefore, our results suggest that acquisition of new chromosomal materials has not been of major importance in the dramatic change of life cycle that has accompanied the emergence of Y. pestis.

The sudden dominance of blaCTX-M harbouring plasmids in Shigella spp. Circulating in Southern Vietnam

BACKGROUND

Plasmid mediated antimicrobial resistance in the Enterobacteriaceae is a global problem. The rise of CTX-M class extended spectrum beta lactamases (ESBLs) has been well documented in industrialized countries. Vietnam is representative of a typical transitional middle income country where the spectrum of infectious diseases combined with the spread of drug resistance is shifting and bringing new healthcare challenges.

METHODOLOGY

We collected hospital admission data from the pediatric population attending the hospital for tropical diseases in Ho Chi Minh City with Shigella infections. Organisms were cultured from all enrolled patients and subjected to antimicrobial susceptibility testing. Those that were ESBL positive were subjected to further investigation. These investigations included PCR amplification for common ESBL genes, plasmid investigation, conjugation, microarray hybridization and DNA sequencing of a bla(CTX-M) encoding plasmid.

PRINCIPAL FINDINGS

We show that two different bla(CTX-M) genes are circulating in this bacterial population in this location. Sequence of one of the ESBL plasmids shows that rather than the gene being integrated into a preexisting MDR plasmid, the bla(CTX-M) gene is located on relatively simple conjugative plasmid. The sequenced plasmid (pEG356) carried the bla(CTX-M-24) gene on an ISEcp1 element and demonstrated considerable sequence homology with other IncFI plasmids.

SIGNIFICANCE

The rapid dissemination, spread of antimicrobial resistance and changing population of Shigella spp. concurrent with economic growth are pertinent to many other countries undergoing similar development. Third generation cephalosporins are commonly used empiric antibiotics in Ho Chi Minh City. We recommend that these agents should not be considered for therapy of dysentery in this setting.

A critical assessment of cross-species detection of gene duplicates using comparative genomic hybridization

Background

Comparison of genomic DNA among closely related strains or species is a powerful approach for identifying variation in evolutionary processes. One potent source of genomic variation is gene duplication, which is prevalent among individuals and species. Array comparative genomic hybridization (aCGH) has been successfully utilized to detect this variation among lineages. Here, beyond the demonstration that gene duplicates among species can be quantified with aCGH, we consider the effect of sequence divergence on the ability to detect gene duplicates.

Results

Using the X chromosome genomic content difference between male D. melanogaster and female D. yakuba and D. simulans, we describe a decrease in the ability to accurately measure genomic content (copy number) for orthologs that are only 90% identical. We demonstrate that genome characteristics (e.g. chromatin environment and non-orthologous sequence similarity) can also affect the ability to accurately measure genomic content. We describe a normalization strategy and statistical criteria to be used for the identification of gene duplicates among any species group for which an array platform is available from a closely related species.

Conclusions

Array CGH can be used to effectively identify gene duplication and genome content; however, certain biases are present due to sequence divergence and other genome characteristics resulting from the divergence between lineages. Highly conserved gene duplicates will be more readily recovered by aCGH. Duplicates that have been retained for a selective advantage due to directional selection acting on many loci in one or both gene copies are likely to be under-represented. The results of this study should inform the interpretation of both previously published and future work that employs this powerful technique.

Using comparative genomic hybridization to survey genomic sequence divergence across species: a proof-of-concept from Drosophila

Background

Genome-wide analysis of sequence divergence among species offers profound insights into the evolutionary processes that shape lineages. When full-genome sequencing is not feasible for a broad comparative study, we propose the use of array-based comparative genomic hybridization (aCGH) in order to identify orthologous genes with high sequence divergence. Here we discuss experimental design, statistical power, success rate, sources of variation and potential confounding factors. We used a spotted PCR product microarray platform from Drosophila melanogaster to assess sequence divergence on a gene-by-gene basis in three fully sequenced heterologous species (D. sechellia, D. simulans, and D. yakuba). Because complete genome assemblies are available for these species this study presents a powerful test for the use of aCGH as a tool to measure sequence divergence.

Results

We found a consistent and linear relationship between hybridization ratio and sequence divergence of the sample to the platform species. At higher levels of sequence divergence (< 92% sequence identity to D. melanogaster) ~84% of features had significantly less hybridization to the array in the heterologous species than the platform species, and thus could be identified as "diverged". At lower levels of divergence (≥ 97% identity), only 13% of genes were identified as diverged. While ~40% of the variation in hybridization ratio can be accounted for by variation in sequence identity of the heterologous sample relative to D. melanogaster, other individual characteristics of the DNA sequences, such as GC content, also contribute to variation in hybridization ratio, as does technical variation.

Conclusions

Here we demonstrate that aCGH can accurately be used as a proxy to estimate genome-wide divergence, thus providing an efficient way to evaluate how evolutionary processes and genomic architecture can shape species diversity in non-model systems. Given the increased number of species for which microarray platforms are available, comparative studies can be conducted for many interesting lineages in order to identify highly diverged genes that may be the target of natural selection.

Insights in Bartonella host specificity

The genus Bartonella comprises a unique group of emerging gram-negative, intracellular bacteria that can cause a long-lasting intraerythrocytic bacteremia in their reservoir hosts. In recent years, the widespread occurrence and diversity of these bacteria has been increasingly recognized. This has resulted in a dramatic expansion of the genus Bartonella to 24 currently described species or subspecies, among which at least half have been associated with human disease. Bartonella infections have been observed in virtually all species examined, extending from humans to carnivores, ungulates, rodents, lagomorphs, insectivores, and bats. Adaptation by Bartonellae to such a diverse range of mammals has resulted in host specificity, and all validated Bartonella species described to date are capable of parasitizing only a limited number of animal species. In this review, the possible mechanisms explaining the specificity of each Bartonella species for its reservoir host are discussed.

Genotyping and phylogenetic analysis of Yersinia pestis by MLVA: insights into the worldwide expansion of Central Asia plague foci

BACKGROUND

The species Yersinia pestis is commonly divided into three classical biovars, Antiqua, Medievalis, and Orientalis, belonging to subspecies pestis pathogenic for human and the (atypical) non-human pathogenic biovar Microtus (alias Pestoides) including several non-pestis subspecies. Recent progress in molecular typing methods enables large-scale investigations in the population structure of this species. It is now possible to test hypotheses about its evolution which were proposed decades ago. For instance the three classical biovars of different geographical distributions were suggested to originate from Central Asia. Most investigations so far have focused on the typical pestis subspecies representatives found outside of China, whereas the understanding of the emergence of this human pathogen requires the investigation of strains belonging to subspecies pestis from China and to the Microtus biovar.

METHODOLOGY/PRINCIPAL FINDINGS

Multi-locus VNTR analysis (MLVA) with 25 loci was performed on a collection of Y. pestis isolates originating from the majority of the known foci worldwide and including typical rhamnose-negative subspecies pestis as well as rhamnose-positive subspecies pestis and biovar Microtus. More than 500 isolates from China, the Former Soviet Union (FSU), Mongolia and a number of other foci around the world were characterized and resolved into 350 different genotypes. The data revealed very close relationships existing between some isolates from widely separated foci as well as very high diversity which can conversely be observed between nearby foci.

CONCLUSIONS/SIGNIFICANCE

The results obtained are in full agreement with the view that the Y. pestis subsp. pestis pathogenic for humans emerged in the Central Asia region between China, Kazakhstan, Russia and Mongolia, only three clones of which spread out of Central Asia. The relationships among the strains in China, Central Asia and the rest of the world based on the MLVA25 assay provide an unprecedented view on the expansion and microevolution of Y. pestis.

GO4genome: a prokaryotic phylogeny based on genome organization

Determining the phylogeny of closely related prokaryotes may fail in an analysis of rRNA or a small set of sequences. Whole-genome phylogeny utilizes the maximally available sample space. For a precise determination of genome similarity, two aspects have to be considered when developing an algorithm of whole-genome phylogeny: (1) gene order conservation is a more precise signal than gene content; and (2) when using sequence similarity, failures in identifying orthologues or the in situ replacement of genes via horizontal gene transfer may give misleading results. GO4genome is a new paradigm, which is based on a detailed analysis of gene function and the location of the respective genes. For characterization of genes, the algorithm uses gene ontology enabling a comparison of function independent of evolutionary relationship. After the identification of locally optimal series of gene functions, their length distribution is utilized to compute a phylogenetic distance. The outcome is a classification of genomes based on metabolic capabilities and their organization. Thus, the impact of effects on genome organization that are not covered by methods of molecular phylogeny can be studied. Genomes of strains belonging to Escherichia coli, Shigella, Streptococcus, Methanosarcina, and Yersinia were analyzed. Differences from the findings of classical methods are discussed.

Genomic comparison of Bradyrhizobium japonicum strains with different symbiotic nitrogen-fixing capabilities and other Bradyrhizobiaceae members

Comparative genomic hybridization (CGH) was performed with nine strains of Bradyrhizobium japonicum (a symbiotic nitrogen-fixing bacterium associated with soybean) and eight other members of the Bradyrhizobiaceae by DNA macroarray of B. japonicum USDA110. CGH clearly discriminated genomic variations in B. japonicum strains, but similar CGH patterns were observed in other members of the Bradyrhizobiaceae. The most variable regions were 14 genomic islands (4-97 kb) and low G+C regions on the USDA110 genome, some of which were missing in several strains of B. japonicum and other members of the Bradyrhizobiaceae. The CGH profiles of B. japonicum were classified into three genome types: 110, 122 and 6. Analysis of DNA sequences around the boundary regions showed that at least seven genomic islands were missing in genome type 122 as compared with type 110. Phylogenetic analysis for internal transcribed sequences revealed that strains belonging to genome types 110 and 122 formed separate clades. Thus genomic islands were horizontally inserted into the ancestor genome of type 110 after divergence of the type 110 and 122 strains. To search for functional relationships of variable genomic islands, we conducted linear models of the correlation between the existence of genomic regions and the parameters associated with symbiotic nitrogen fixation in soybean. Variable genomic regions including genomic islands were associated with the enhancement of symbiotic nitrogen fixation in B. japonicum USDA110.

Evolution, population structure, and phylogeography of genetically monomorphic bacterial pathogens

Genetically monomorphic bacteria contain so little sequence diversity that sequencing a few gene fragments yields little or no information. As a result, our understanding of their evolutionary patterns presents greater technical challenges than exist for genetically diverse microbes. These challenges are now being met by analyses at the genomic level for diverse types of genetic variation, the most promising of which are single nucleotide polymorphisms. Many of the most virulent bacterial pathogens are genetically monomorphic, and understanding their evolutionary and phylogeographic patterns will help our understanding of the effects of infectious disease on human history.

YtxR acts as an overriding transcriptional off switch for the Yersinia enterocolitica Ysc-Yop type 3 secretion system

The Yersinia enterocolitica YtxR protein is a LysR-type transcriptional regulator that induces expression of the ytxAB locus, which encodes a putative ADP-ribosylating toxin. The ytxR and ytxAB genes are not closely linked in the Y. enterocolitica chromosome, and whereas ytxR is present in all sequenced Yersinia spp., the ytxAB locus is not. These observations suggested that there might be other YtxR-regulon members besides ytxAB and prompted us to investigate coregulated genes and gene products by using transcriptional and proteomic approaches. Microarray and reverse transcription-PCR analysis showed that YtxR strongly activates expression of the yts2 locus, which encodes a putative type 2 secretion system, as well as several uncharacterized genes predicted to encode extracytoplasmic proteins. Strikingly, we also discovered that under Ysc-Yop type 3 secretion system-inducing conditions, YtxR prevented the appearance of Yop proteins in the culture supernatant. Microarray and lacZ operon fusion analysis showed that this was due to specific repression of ysc-yop gene expression. YtxR was also able to repress VirF-dependent Phi(yopE-lacZ) and Phi(yopH-lacZ) expression in a strain lacking the virulence plasmid, which suggested a direct repression mechanism. This was supported by DNase I footprinting, which showed that YtxR interacted with the yopE and yopH control regions. Therefore, YtxR is a newly identified regulator of the ysc-yop genes that can act as an overriding off switch for this critical virulence system.

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.

Development and application of the active surveillance of pathogens microarray to monitor bacterial gene flux

Background

Human and animal health is constantly under threat by emerging pathogens that have recently acquired genetic determinants that enhance their survival, transmissibility and virulence. We describe the construction and development of an Active Surveillance of Pathogens (ASP) oligonucleotide microarray, designed to 'actively survey' the genome of a given bacterial pathogen for virulence-associated genes.

Results

The microarray consists of 4958 reporters from 151 bacterial species and include genes for the identification of individual bacterial species as well as mobile genetic elements (transposons, plasmid and phage), virulence genes and antibiotic resistance genes. The ASP microarray was validated with nineteen bacterial pathogens species, including Francisella tularensis, Clostridium difficile, Staphylococcus aureus, Enterococcus faecium and Stenotrophomonas maltophilia. The ASP microarray identified these bacteria, and provided information on potential antibiotic resistance (eg sufamethoxazole resistance and sulfonamide resistance) and virulence determinants including genes likely to be acquired by horizontal gene transfer (e.g. an alpha-haemolysin).

Conclusion

The ASP microarray has potential in the clinic as a diagnostic tool, as a research tool for both known and emerging pathogens, and as an early warning system for pathogenic bacteria that have been recently modified either naturally or deliberately.

Characterization of chromosomal regions conserved in Yersinia pseudotuberculosis and lost by Yersinia pestis

The transformation of the enteropathogenic bacterium Yersinia pseudotuberculosis into the plague bacillus, Yersinia pestis, has been accompanied by extensive genetic loss. This study focused on chromosomal regions conserved in Y. pseudotuberculosis and lost during its transformation into Y. pestis. An extensive PCR screening of 78 strains of the two species identified five regions (R1 to R5) and four open reading frames (ORFs; orf1 to orf4) that were conserved in Y. pseudotuberculosis and absent from Y. pestis. Their conservation in Y. pseudotuberculosis suggests a positive selective pressure and a role during the life cycle of this species. Attempts to delete two ORFs (orf3 and orf4) from the chromosome of strain IP32953 were unsuccessful, indicating that they are essential for its viability. The seven remaining loci were individually deleted from the IP32953 chromosome, and the ability of each mutant to grow in vitro and to kill mice upon intragastric infection was evaluated. Four loci (orf1, R2, R4, and R5) were not required for optimal growth or virulence of Y. pseudotuberculosis. In contrast, orf2, encoding a putative pseudouridylate synthase involved in RNA stability, was necessary for the optimal growth of IP32953 at 37 degrees C in a chemically defined medium (M63S). Deletion of R1, a region predicted to encode the methionine salvage pathway, altered the mutant pathogenicity, suggesting that the availability of free methionine is severely restricted in vivo. R3, a region composed mostly of genes of unknown functions, was necessary for both optimal growth of Y. pseudotuberculosis at 37 degrees C in M63S and for virulence. Therefore, despite their loss in Y. pestis, five of the nine Y. pseudotuberculosis-specific chromosomal loci studied play a role in the survival, growth, or virulence of this species.

Different region analysis for genotyping Yersinia pestis isolates from China

Background

DFR (different region) analysis has been developed for typing Yesinia pestis in our previous study, and in this study, we extended this method by using 23 DFRs to investigate 909 Chinese Y. pestis strains for validating DFR-based genotyping method and better understanding adaptive microevolution of Y. pestis.

Methodology/Principal Findings

On the basis of PCR and Bionumerics data analysis, 909 Y. pestis strains were genotyped into 32 genomovars according to their DFR profiles. New terms, Major genomovar and Minor genomovar, were coined for illustrating evolutionary relationship between Y. pestis strains from different plague foci and different hosts. In silico DFR profiling of the completed or draft genomes shed lights on the evolutionary scenario of Y. pestis from Y. pseudotuberculosis. Notably, several sequenced Y. pestis strains share the same DFR profiles with Chinese strains, providing data for revealing the global plague foci expansion.

Conclusions/significance

Distribution of Y. pestis genomovars is plague focus-specific. Microevolution of biovar Orientalis was deduced according to DFR profiles. DFR analysis turns to be an efficient and inexpensive method to portrait the genome plasticity of Y. pestis based on horizontal gene transfer (HGT). DFR analysis can also be used as a tool in comparative and evolutionary genomic research for other bacteria with similar genome plasticity.

Characterisation of Yersinia pestis isolates from natural foci of plague in the Republic of Georgia, and their relationship to Y. pestis isolates from other countries

Forty Yersinia pestis isolates from endemic foci of plague in the Republic of Georgia, and six Y. pestis isolates from neighbouring former Soviet Union countries, were analysed for their biochemical and phenotypic properties, and their genetic relatedness was compared with Y. pestis strains KIM and CO92 by pulsed-field gel electrophoresis (PFGE). In addition, 11 Y. pestis isolates from the USA, together with published nucleotide sequences from Y. pestis strains KIM, CO92 and 91001, were compared with the 46 isolates in the present collection using multilocus sequence typing (MLST), based on sequence data for the 16S rRNA, hsp60, glnA, gyrB, recA, manB, thrA and tmk loci. Four virulence gene loci (caf1, lcrV, psaA and pla) were also sequenced and analysed. Two sequence types (ST1 and ST2), which differed by a single nucleotide, were identified by MLST. With the exception of a single isolate (771G), all of the Georgian Y. pestis isolates belonged to ST2. PFGE also grouped the Georgian Y. pestis isolates separately from the non-Georgian isolates. Overall, PFGE discriminated the Y. pestis isolates more effectively than MLST. The sequences of three of the four virulence genes (lcrV, psaA and pla) were identical in all Georgian and non-Georgian isolates, but the caf1 locus was represented by two allele types, with caf1 NT1 being associated with the non-Georgian isolates and caf1 NT2 being associated with the Georgian isolates. These results suggest that Georgian Y. pestis isolates are of clonal origin.