Common themes in the genome strategies of pathogens

Systematic comparison of genome information processing and boundary recognition tools used for genomic island detection

Genomic islands are fragments of foreign DNA that are found in bacterial and archaeal genomes, and are typically associated with symbiosis or pathogenesis. While numerous genomic island detection methods have been proposed, there has been limited evaluation of the efficiency of the genome information processing and boundary recognition tools. In this study, we conducted a review of the statistical methods involved in genomic signatures, host signature extraction, informative signature selection, divergence measures, and boundary detection steps in genomic island prediction. We compared the performances of these methods on simulated experiments using alien fragments obtained from both artificial and real genomes. Our results indicate that among the nine genomic signatures evaluated, genomic signature frequency and full probability performed the best. However, their performance declined when normalized to their expectations and variances, such as Z-score and composition vector. Based on our experiments of the E. coli genome, we found that the confidence intervals of the window variances achieved the best performance in the signature extraction of the host, with the best confidence interval being 1.5-2 times the standard error. Ordered kurtosis was most effective in selecting informative signatures from a single genome, without requiring prior knowledge from other datasets. Among the three divergence measures evaluated, the two-sample t-test was the most successful, and a non-overlapping window with a small eye window (size 2) was best suited for identifying compositionally distinct regions. Finally, the maximum of the Markovian Jensen-Shannon divergence score, in terms of GC-content bias, was found to make boundary detection faster while maintaining a similar error rate.

The pangenome of the wheat pathogen Pyrenophora tritici-repentis reveals novel transposons associated with necrotrophic effectors ToxA and ToxB

BACKGROUND

In fungal plant pathogens, genome rearrangements followed by selection pressure for adaptive traits have facilitated the co-evolutionary arms race between hosts and their pathogens. Pyrenophora tritici-repentis (Ptr) has emerged recently as a foliar pathogen of wheat worldwide and its populations consist of isolates that vary in their ability to produce combinations of different necrotrophic effectors. These effectors play vital roles in disease development. Here, we sequenced the genomes of a global collection (40 isolates) of Ptr to gain insights into its gene content and genome rearrangements.

RESULTS

A comparative genome analysis revealed an open pangenome, with an abundance of accessory genes (~ 57%) reflecting Ptr's adaptability. A clear distinction between pathogenic and non-pathogenic genomes was observed in size, gene content, and phylogenetic relatedness. Chromosomal rearrangements and structural organization, specifically around effector coding genes, were detailed using long-read assemblies (PacBio RS II) generated in this work in addition to previously assembled genomes. We also discovered the involvement of large mobile elements associated with Ptr's effectors: ToxA, the gene encoding for the necrosis effector, was found as a single copy within a 143-kb 'Starship' transposon (dubbed 'Horizon') with a clearly defined target site and target site duplications. 'Horizon' was located on different chromosomes in different isolates, indicating mobility, and the previously described ToxhAT transposon (responsible for horizontal transfer of ToxA) was nested within this newly identified Starship. Additionally, ToxB, the gene encoding the chlorosis effector, was clustered as three copies on a 294-kb element, which is likely a different putative 'Starship' (dubbed 'Icarus') in a ToxB-producing isolate. ToxB and its putative transposon were missing from the ToxB non-coding reference isolate, but the homolog toxb and 'Icarus' were both present in a different non-coding isolate. This suggests that ToxB may have been mobile at some point during the evolution of the Ptr genome which is contradictory to the current assumption of ToxB vertical inheritance. Finally, the genome architecture of Ptr was defined as 'one-compartment' based on calculated gene distances and evolutionary rates.

CONCLUSIONS

These findings together reflect on the highly plastic nature of the Ptr genome which has likely helped to drive its worldwide adaptation and has illuminated the involvement of giant transposons in facilitating the evolution of virulence in Ptr.

The Genomics and Cell Biology of Host-Beneficial Intracellular Infections

Microbes gain access to eukaryotic cells as food for bacteria-grazing protists, for host protection by microbe-killing immune cells, or for microbial benefit when pathogens enter host cells to replicate. But microbes can also gain access to a host cell and become an important-often required-beneficial partner. The oldest beneficial microbial infections are the ancient eukaryotic organelles now called the mitochondrion and plastid. But numerous other host-beneficial intracellular infections occur throughout eukaryotes. Here I review the genomics and cell biology of these interactions with a focus on intracellular bacteria. The genomes of host-beneficial intracellular bacteria have features that span a previously unfilled gap between pathogens and organelles. Host cell adaptations to allow the intracellular persistence of beneficial bacteria are found along with evidence for the microbial manipulation of host cells, but the cellular mechanisms of beneficial bacterial infections are not well understood. Expected final online publication date for the Annual Review of Cell and Developmental Biology, Volume 37 is October 2021. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Genomic Island Prediction via Chi-Square Test and Random Forest Algorithm

Genomic islands are related to microbial adaptation and carry different genomic characteristics from the host. Therefore, many methods have been proposed to detect genomic islands from the rest of the genome by evaluating its sequence composition. Many sequence features have been proposed, but many of them have not been applied to the identification of genomic islands. In this paper, we present a scheme to predict genomic islands using the chi-square test and random forest algorithm. We extract seven kinds of sequence features and select the important features with the chi-square test. All the selected features are then input into the random forest to predict the genome islands. Three experiments and comparison show that the proposed method achieves the best performance. This understanding can be useful to design more powerful method for the genomic island prediction.

The key to egress? Babesia bovis perforin-like protein 1 (PLP1) with hemolytic capacity is required for blood stage replication and is involved in the exit of the parasite from the host cell

Bovine babesiosis is a tick-borne disease caused by apicomplexan parasites of the Babesia genus that represents a major constraint to livestock production worldwide. Currently available vaccines are based on live parasites which have archetypal limitations. Our goal is to identify candidate antigens so that new and effective vaccines against Babesia may be developed. The perforin-like protein (PLP) family has been identified as a key player in cell traversal and egress in related apicomplexans and it was also identified in Babesia, but its function in this parasite remains unknown. The aim of this work was to define the PLP family in Babesia and functionally characterize PLP1, a representative member of the family in Babesia bovis. Bioinformatic analyses demonstrate a variable number of plp genes (four to eight) in the genomes of six different Babesia spp. and conservation of the family members at the secondary and tertiary structure levels. We demonstrate here that Babesia PLPs contain the critical domains present in other apicomplexan PLPs to display the lytic capacity. We then focused on the functional characterization of PLP1 of B. bovis, both in vitro and in vivo. PLP1 is expressed and exposed to the host immune system during infection and has high hemolytic capacity under a wide range of conditions in vitro. A B. bovis plp1 knockout line displayed a decreased growth rate in vitro compared with the wild type strain and a peculiar phenotype consisting of multiple parasites within a single red blood cell, although at low frequency. This phenotype suggests that the lack of PLP1 has a negative impact on the mechanism of egression of the parasite and, therefore, on its capacity to proliferate. It is possible that PLP1 is associated with other proteins in the processes of invasion and egress, which were found to have redundant mechanisms in related apicomplexans. Future work will be focused on unravelling the network of proteins involved in these essential parasite functions.

2SigFinder: the combined use of small-scale and large-scale statistical testing for genomic island detection from a single genome

BACKGROUND

Genomic islands are associated with microbial adaptations, carrying genomic signatures different from the host. Some methods perform an overall test to identify genomic islands based on their local features. However, regions of different scales will display different genomic features.

RESULTS

We proposed here a novel method "2SigFinder ", the first combined use of small-scale and large-scale statistical testing for genomic island detection. The proposed method was tested by genomic island boundary detection and identification of genomic islands or functional features of real biological data. We also compared the proposed method with the comparative genomics and composition-based approaches. The results indicate that the proposed 2SigFinder is more efficient in identifying genomic islands.

CONCLUSIONS

From real biological data, 2SigFinder identified genomic islands from a single genome and reported robust results across different experiments, without annotated information of genomes or prior knowledge from other datasets. 2SigHunter identified 25 Pathogenicity, 1 tRNA, 2 Virulence and 2 Repeats from 27 Pathogenicity, 1 tRNA, 2 Virulence and 2 Repeats, and detected 101 Phage and 28 HEG out of 130 Phage and 36 HEGs in S. enterica Typhi CT18, which shows that it is more efficient in detecting functional features associated with GIs.

Formylated N-terminal methionine is absent from the Mycoplasma hyopneumoniae proteome: Implications for translation initiation

N-terminal methionine excision (NME) is a proteolytic pathway that cleaves the N-termini of proteins, a process that influences where proteins localise in the cell and their turnover rates. In bacteria, protein biosynthesis is initiated by formylated methionine start tRNA (fMet-tRNA_f^Met). The formyl group is attached by formyltransferase (FMT) and is subsequently removed by peptide deformylase (PDF) in most but not all proteins. Methionine aminopeptidase then cleaves deformylated methionine to complete the process. Components of NME, particularly PDF, are promising therapeutic targets for bacterial pathogens. In Mycoplasma hyopneumoniae, a genome-reduced, major respiratory pathogen of swine, pdf and fmt are absent from its genome. Our bioinformatic analysis uncovered additional enzymes involved in formylated N-terminal methionine (fnMet) processing missing in fourteen mycoplasma species, including M. hyopneumoniae but not in Mycoplasma pneumoniae, a major respiratory pathogen of humans. Consistent with our bioinformatic studies, an analysis of in-house tryptic peptide libraries confirmed the absence of fnMet in M. hyopneumoniae proteins but, as expected fnMet peptides were detected in the proteome of M. pneumoniae. Additionally, computational molecular modelling of M. hyopneumoniae translation initiation factors reveal structural and sequence differences in areas known to interact with fMet-tRNA_f^Met. Our data suggests that some mycoplasmas have evolved a translation process that does not require fnMet.

MTGIpick allows robust identification of genomic islands from a single genome

Genomic islands (GIs) that are associated with microbial adaptations and carry sequence patterns different from that of the host are sporadically distributed among closely related species. This bias can dominate the signal of interest in GI detection. However, variations still exist among the segments of the host, although no uniform standard exists regarding the best methods of discriminating GIs from the rest of the genome in terms of compositional bias. In the present work, we proposed a robust software, MTGIpick, which used regions with pattern bias showing multiscale difference levels to identify GIs from the host. MTGIpick can identify GIs from a single genome without annotated information of genomes or prior knowledge from other data sets. When real biological data were used, MTGIpick demonstrated better performance than existing methods, as well as revealed potential GIs with accurate sizes missed by existing methods because of a uniform standard. Software and supplementary are freely available at http://bioinfo.zstu.edu.cn/MTGI or https://github.com/bioinfo0706/MTGIpick.

An efficient genomic signature ranking method for genomic island prediction from a single genome

Genomic islands that are associated with microbial adaptations and carry genomic signatures different from that of the host, and thus many methods have been proposed to select the informative genomic signatures from a range of organisms and discriminate genomic islands from the rest of the genome in terms of these signature biases. However, they are of limited use when closely related genomes are unavailable. In the present work, we proposed a kurtosis-based ranking method to select the informative genomic signatures from a single genome. In simulations with alien fragments from artificial and real genomes, the proposed kurtosis-based ranking method efficiently selected the informative genomic signatures from a single genome, without annotated information of genomes or prior knowledge from other datasets. This understanding can be useful to design more powerful method for genomic island detection.

Lineage-specific plasmid acquisition and the evolution of specialized pathogens in Bacillus thuringiensis and the Bacillus cereus group

Bacterial plasmids can vary from small selfish genetic elements to large autonomous replicons that constitute a significant proportion of total cellular DNA. By conferring novel function to the cell, plasmids may facilitate evolution but their mobility may be opposed by co-evolutionary relationships with chromosomes or encouraged via the infectious sharing of genes encoding public goods. Here, we explore these hypotheses through large-scale examination of the association between plasmids and chromosomal DNA in the phenotypically diverse Bacillus cereus group. This complex group is rich in plasmids, many of which encode essential virulence factors (Cry toxins) that are known public goods. We characterized population genomic structure, gene content and plasmid distribution to investigate the role of mobile elements in diversification. We analysed coding sequence within the core and accessory genome of 190 B. cereus group isolates, including 23 novel sequences and genes from 410 reference plasmid genomes. While cry genes were widely distributed, those with invertebrate toxicity were predominantly associated with one sequence cluster (clade 2) and phenotypically defined Bacillus thuringiensis. Cry toxin plasmids in clade 2 showed evidence of recent horizontal transfer and variable gene content, a pattern of plasmid segregation consistent with transfer during infectious cooperation. Nevertheless, comparison between clades suggests that co-evolutionary interactions may drive association between plasmids and chromosomes and limit wider transfer of key virulence traits. Proliferation of successful plasmid and chromosome combinations is a feature of specialized pathogens with characteristic niches (Bacillus anthracis, B. thuringiensis) and has occurred multiple times in the B. cereus group.

Comparative Genomics of Pathogenic and Nonpathogenic Beetle-Vectored Fungi in the Genus Geosmithia

Geosmithia morbida is an emerging fungal pathogen which serves as a model for examining the evolutionary processes behind pathogenicity because it is one of two known pathogens within a genus of mostly saprophytic, beetle-associated, fungi. This pathogen causes thousand cankers disease in black walnut trees and is vectored into the host via the walnut twig beetle. Geosmithia morbida was first detected in western United States and currently threatens the timber industry concentrated in eastern United States. We sequenced the genomes of G. morbida in a previous study and two nonpathogenic Geosmithia species in this work and compared these species to other fungal pathogens and nonpathogens to identify genes under positive selection in G. morbida that may be associated with pathogenicity. Geosmithia morbida possesses one of the smallest genomes among the fungal species observed in this study, and one of the smallest fungal pathogen genomes to date. The enzymatic profile in this pathogen is very similar to its nonpathogenic relatives. Our findings indicate that genome reduction or retention of a smaller genome may be an important adaptative force during the evolution of a specialized lifestyle in fungal species that occupy a specificniche, such as beetle vectored tree pathogens. We also present potential genes under selection in G. morbida that could be important for adaptation to a pathogenic lifestyle.

Lack of the PGA exopolysaccharide in Salmonella as an adaptive trait for survival in the host

Many bacteria build biofilm matrices using a conserved exopolysaccharide named PGA or PNAG (poly-β-1,6-N-acetyl-D-glucosamine). Interestingly, while E. coli and other members of the family Enterobacteriaceae encode the pgaABCD operon responsible for PGA synthesis, Salmonella lacks it. The evolutionary force driving this difference remains to be determined. Here, we report that Salmonella lost the pgaABCD operon after the divergence of Salmonella and Citrobacter clades, and previous to the diversification of the currently sequenced Salmonella strains. Reconstitution of the PGA machinery endows Salmonella with the capacity to produce PGA in a cyclic dimeric GMP (c-di-GMP) dependent manner. Outside the host, the PGA polysaccharide does not seem to provide any significant benefit to Salmonella: resistance against chlorine treatment, ultraviolet light irradiation, heavy metal stress and phage infection remained the same as in a strain producing cellulose, the main biofilm exopolysaccharide naturally produced by Salmonella. In contrast, PGA production proved to be deleterious to Salmonella survival inside the host, since it increased susceptibility to bile salts and oxidative stress, and hindered the capacity of S. Enteritidis to survive inside macrophages and to colonize extraintestinal organs, including the gallbladder. Altogether, our observations indicate that PGA is an antivirulence factor whose loss may have been a necessary event during Salmonella speciation to permit survival inside the host.

During bacterial evolution, specific traits that optimize the organism’s fitness are selected. The production of exopolysaccharides is widespread among bacteria in which they play a protective shielding role as main constituents of biofilms. In contrast to closely related siblings, Salmonella has lost the capacity to produce the exopolysaccharide PGA. Our study reveals that Salmonella lost pga genes, and that the driving force for such a loss may have been the detrimental impact that PGA has during Salmonella invasion of internal organs where it augments the susceptibility to bile salts and oxygen radicals, reducing bacterial survival inside macrophages and rendering Salmonella avirulent. These results suggest that gene-loss has played an important role during Salmonella evolution.

Genome Report: Identification and Validation of Antigenic Proteins from Pajaroellobacter abortibovis Using De Novo Genome Sequence Assembly and Reverse Vaccinology

Epizootic bovine abortion (EBA), or “foothill abortion,” is the leading cause of beef cattle abortion in California and has also been reported in Nevada and Oregon. In the 1970s, the soft-shelled tick Ornithodoros coriaceus, or “pajaroello tick,” was confirmed as the disease-transmitting vector. In 2005, a novel Deltaproteobacterium was discovered as the etiologic agent of EBA (aoEBA), recently named Pajaroellobacter abortibovis. This organism cannot be grown in culture using traditional microbiological techniques; it can only be grown in experimentally-infected severe combined immunodeficient (SCID) mice. The objectives of this study were to perform a de novo genome assembly for P. abortibovis and identify and validate potential antigenic proteins as candidates for future recombinant vaccine development. DNA and RNA were extracted from spleen tissue collected from experimentally-infected SCID mice following exposure to P. abortibovis. This combination of mouse and bacterial DNA was sequenced and aligned to the mouse genome. Mouse sequences were subtracted from the sequence pool and the remaining sequences were de novo assembled at 50x coverage into a 1.82 Mbp complete closed circular Deltaproteobacterial genome containing 2250 putative protein-coding sequences. Phylogenetic analysis of P. abortibovis predicts that this bacterium is most closely related to the organisms of the order Myxococcales, referred to as Myxobacteria. In silico prediction of vaccine candidates was performed using a reverse vaccinology approach resulting in the identification and ranking of the top 10 candidate proteins that are likely to be antigenic. Immunologic testing of these candidate proteins confirmed antigenicity of seven of the nine expressed protein candidates using serum from P. abortibovis immunized mice.

A method for identification of highly conserved elements and evolutionary analysis of superphylum Alveolata

Background

Perfectly or highly conserved DNA elements were found in vertebrates, invertebrates, and plants by various methods. However, little is known about such elements in protists. The evolutionary distance between apicomplexans can be very high, in particular, due to the positive selection pressure on them. This complicates the identification of highly conserved elements in alveolates, which is overcome by the proposed algorithm.

Results

A novel algorithm is developed to identify highly conserved DNA elements. It is based on the identification of dense subgraphs in a specially built multipartite graph (whose parts correspond to genomes). Specifically, the algorithm does not rely on genome alignments, nor pre-identified perfectly conserved elements; instead, it performs a fast search for pairs of words (in different genomes) of maximum length with the difference below the specified edit distance. Such pair defines an edge whose weight equals the maximum (or total) length of words assigned to its ends. The graph composed of these edges is then compacted by merging some of its edges and vertices. The dense subgraphs are identified by a cellular automaton-like algorithm; each subgraph defines a cluster composed of similar inextensible words from different genomes. Almost all clusters are considered as predicted highly conserved elements. The algorithm is applied to the nuclear genomes of the superphylum Alveolata, and the corresponding phylogenetic tree is built and discussed.

Conclusion

We proposed an algorithm for the identification of highly conserved elements. The multitude of identified elements was used to infer the phylogeny of Alveolata.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-016-1257-5) contains supplementary material, which is available to authorized users.

Ammonia-Oligotrophic and Diazotrophic Heavy Metal-Resistant Serratia liquefaciens Strains from Pioneer Plants and Mine Tailings

Mine tailings are man-made environments characterized by low levels of organic carbon and assimilable nitrogen, as well as moderate concentrations of heavy metals. For the introduction of nitrogen into these environments, a key role is played by ammonia-oligotrophic/diazotrophic heavy metal-resistant guilds. In mine tailings from Zacatecas, Mexico, Serratia liquefaciens was the dominant heterotrophic culturable species isolated in N-free media from bulk mine tailings as well as the rhizosphere, roots, and aerial parts of pioneer plants. S. liquefaciens strains proved to be a meta-population with high intraspecific genetic diversity and a potential to respond to these extreme conditions. The phenotypic and genotypic features of these strains reveal the potential adaptation of S. liquefaciens to oligotrophic and nitrogen-limited mine tailings with high concentrations of heavy metals. These features include ammonia-oligotrophic growth, nitrogen fixation, siderophore and indoleacetic acid production, phosphate solubilization, biofilm formation, moderate tolerance to heavy metals under conditions of diverse nitrogen availability, and the presence of zntA, amtB, and nifH genes. The acetylene reduction assay suggests low nitrogen-fixing activity. The nifH gene was harbored in a plasmid of ∼60 kb and probably was acquired by a horizontal gene transfer event from Klebsiella variicola.

The Phylogeny of Bacillus cereus sensu lato

The three main species of the Bacillus cereus sensu lato , B. cereus , B. thuringiensis , and B. anthracis , were recognized and established by the early 1900s because they each exhibited distinct phenotypic traits. B. thuringiensis isolates and their parasporal crystal proteins have long been established as a natural pesticide and insect pathogen. B. anthracis , the etiological agent for anthrax, was used by Robert Koch in the 19th century as a model to develop the germ theory of disease, and B. cereus , a common soil organism, is also an occasional opportunistic pathogen of humans. In addition to these three historical species designations, are three less-recognized and -understood species: B. mycoides , B. weihenstephanensis , and B. pseudomycoides . All of these “species” combined comprise the Bacillus cereus sensu lato group. Despite these apparently clear phenotypic definitions, early molecular approaches to separate the first three by various DNA hybridization and 16S/23S ribosomal sequence analyses led to some “confusion” because there were limited differences to differentiate between these species. These and other results have led to frequent suggestions that a taxonomic change was warranted to reclassify this group to a single species. But the pathogenic properties of B. anthracis and the biopesticide applications of B. thuringiensis appear to “have outweighed pure taxonomic considerations” and the separate species categories are still being maintained. B. cereus sensu lato represents a classic example of a now common bacterial species taxonomic quandary.

Early emergence of Yersinia pestis as a severe respiratory pathogen

Yersinia pestis causes the fatal respiratory disease pneumonic plague. Y. pestis recently evolved from the gastrointestinal pathogen Y. pseudotuberculosis; however, it is not known at what point Y. pestis gained the ability to induce a fulminant pneumonia. Here we show that the acquisition of a single gene encoding the protease Pla was sufficient for the most ancestral, deeply rooted strains of Y. pestis to cause pneumonic plague, indicating that Y. pestis was primed to infect the lungs at a very early stage in its evolution. As Y. pestis further evolved, modern strains acquired a single amino-acid modification within Pla that optimizes protease activity. While this modification is unnecessary to cause pneumonic plague, the substitution is instead needed to efficiently induce the invasive infection associated with bubonic plague. These findings indicate that Y. pestis was capable of causing pneumonic plague before it evolved to optimally cause invasive infections in mammals.

Distinct mutations led to inactivation of type 1 fimbriae expression in Shigella spp

Shigella spp. are responsible for bacillary dysentery in humans. The acquisition or the modification of the virulence plasmid encoding factors promoting entry of bacteria into and dissemination within epithelial cells was a critical step in the evolution of these bacteria from their Escherichia coli ancestor(s). Incorporation of genomic islands (GI) and gene inactivation also shaped interactions between these pathogens and their human host. Sequence analysis of the GI inserted next to the leuX tRNA gene in S. boydii, S. dysenteriae, S. flexneri, S. sonnei and enteroinvasive E. coli (EIEC) suggests that this region initially carried the fec, yjhATS and fim gene clusters. The fim cluster encoding type I fimbriae is systematically inactivated in both reference strains and clinical isolates and distinct mutations are responsible for this inactivation in at least three phylogenetic groups. To investigate consequences of the presence of fimbriae on the outcome of the interaction of Shigella with host cells, we used a S. flexneri strain harboring a plasmid encoding the E. coli fim operon. Production of fimbriae by this recombinant strain increased the ability of bacteria to adhere to and enter into epithelial cells and had no effect on their ability to disseminate from cell to cell. The observations that production of type I fimbriae increases invasion of epithelial cells and that independent mutations abolish fimbriae production in Shigella suggest that these mutations correspond to pathoadaptive events.

Prevalence of virulence determinants and antimicrobial resistance among commensal Escherichia coli derived from dairy and beef cattle

Cattle is a reservoir of potentially pathogenic E. coli, bacteria that can represent a significant threat to public health, hence it is crucial to monitor the prevalence of the genetic determinants of virulence and antimicrobial resistance among the E. coli population. The aim of this study was the analysis of the phylogenetic structure, distribution of virulence factors (VFs) and prevalence of antimicrobial resistance among E. coli isolated from two groups of healthy cattle: 50 cows housed in the conventional barn (147 isolates) and 42 cows living on the ecological pasture (118 isolates). The phylogenetic analysis, identification of VFs and antimicrobial resistance genes were based on either multiplex or simplex PCR. The antimicrobial susceptibilities of E. coli were examined using the broth microdilution method. Two statistical approaches were used to analyse the results obtained for two groups of cattle. The relations between the dependent (VFs profiles, antibiotics) and the independent variables were described using the two models. The mixed logit model was used to characterise the prevalence of the analysed factors in the sets of isolates. The univariate logistic regression model was used to characterise the prevalence of these factors in particular animals. Given each model, the odds ratio (OR) and the 95% confidence interval for the population were estimated. The phylogroup B1 was predominant among isolates from beef cattle, while the phylogroups A, B1 and D occurred with equal frequency among isolates from dairy cattle. The frequency of VFs-positive isolates was significantly higher among isolates from beef cattle. E. coli from dairy cattle revealed significantly higher resistance to antibiotics. Some of the tested resistance genes were present among isolates from dairy cattle. Our study showed that the habitat and diet may affect the genetic diversity of commensal E. coli in the cattle. The results suggest that the ecological pasture habitat is related to the increased spreading rate of the VFs, while the barn habitat is characterised by the higher levels of antimicrobial resistance among E. coli.

Turnover and accumulation of genetic diversity across large time-scale cycles of isolation and connection of populations

Major climatic and geological events but also population history (secondary contacts) have generated cycles of population isolation and connection of long and short periods. Recent empirical and theoretical studies suggest that fast evolutionary processes might be triggered by such events, as commonly illustrated in ecology by the adaptive radiation of cichlid fishes (isolation and reconnection of lakes and watersheds) and in epidemiology by the fast adaptation of the influenza virus (isolation and reconnection in hosts). We test whether cyclic population isolation and connection provide the raw material (standing genetic variation) for species evolution and diversification. Our analytical results demonstrate that population isolation and connection can provide, to populations, a high excess of genetic diversity compared with what is expected at equilibrium. This excess is either cyclic (high allele turnover) or cumulates with time depending on the duration of the isolation and the connection periods and the mutation rate. We show that diversification rates of animal clades are associated with specific periods of climatic cycles in the Quaternary. We finally discuss the importance of our results for macroevolutionary patterns and for the inference of population history from genomic data.

Investigating the path of plastid genome degradation in an early-transitional clade of heterotrophic orchids, and implications for heterotrophic angiosperms

Parasitic organisms exemplify morphological and genomic reduction. Some heterotrophic, parasitic plants harbor drastically reduced and degraded plastid genomes resulting from relaxed selective pressure on photosynthetic function. However, few studies have addressed the initial stages of plastome degradation in groups containing both photosynthetic and nonphotosynthetic species. Corallorhiza is a genus of leafless, heterotrophic orchids that contains both green, photosynthetic species and nongreen, putatively nonphotosynthetic species, and represents an ideal system in which to assess the beginning of the transition to a "minimal plastome." Complete plastomes were generated for nine taxa of Corallorhiza using Illumina paired-end sequencing of genomic DNA to assess the degree of degradation among taxa, and for comparison with a general model of degradation among angiosperms. Quantification of total chlorophyll suggests that nongreen Corallorhiza still produce chlorophyll, but at 10-fold lower concentrations than green congeners. Complete plastomes and partial nuclear rDNA cistrons yielded a fully resolved tree for Corallorhiza, with at least two independent losses of photosynthesis, evidenced by gene deletions and pseudogenes in Co. striata and nongreen Co. maculata. All Corallorhiza show some evidence of degradation in genes of the NAD(P)H dehydrogenase complex. Among genes with open reading frames, photosynthesis-related genes displayed evidence of neutral evolution in nongreen Corallorhiza, whereas genes of the ATP synthase complex displayed some evidence of positive selection in these same groups, though for reasons unknown. Corallorhiza spans the early stages of a general model of plastome degradation and has added critical insight for understanding the process of plastome evolution in heterotrophic angiosperms.

Association of clustered regularly interspaced short palindromic repeat (CRISPR) elements with specific serotypes and virulence potential of shiga toxin-producing Escherichia coli

Shiga toxin-producing Escherichia coli (STEC) strains (n = 194) representing 43 serotypes and E. coli K-12 were examined for clustered regularly interspaced short palindromic repeat (CRISPR) arrays to study genetic relatedness among STEC serotypes. A subset of the strains (n = 81) was further analyzed for subtype I-E cas and virulence genes to determine a possible association of CRISPR elements with potential virulence. Four types of CRISPR arrays were identified. CRISPR1 and CRISPR2 were present in all strains tested; 1 strain also had both CRISPR3 and CRISPR4, whereas 193 strains displayed a short, combined array, CRISPR3-4. A total of 3,353 spacers were identified, representing 528 distinct spacers. The average length of a spacer was 32 bp. Approximately one-half of the spacers (54%) were unique and found mostly in strains of less common serotypes. Overall, CRISPR spacer contents correlated well with STEC serotypes, and identical arrays were shared between strains with the same H type (O26:H11, O103:H11, and O111:H11). There was no association identified between the presence of subtype I-E cas and virulence genes, but the total number of spacers had a negative correlation with potential pathogenicity (P < 0.05). Fewer spacers were found in strains that had a greater probability of causing outbreaks and disease than in those with lower virulence potential (P < 0.05). The relationship between the CRISPR-cas system and potential virulence needs to be determined on a broader scale, and the biological link will need to be established.

A unique megaplasmid contributes to stress tolerance and pathogenicity of an emergent Salmonella enterica serovar Infantis strain

Of all known Salmonella enterica serovars, S. Infantis is one of the most commonly isolated and has been recently emerging worldwide. To understand the recent emergence of S. Infantis in Israel, we performed extensive comparative analyses between pre-emergent and the clonal emergent S. Infantis populations. We demonstrate the fixation of adaptive mutations in the DNA gyrase (gyrA) and nitroreductase (nfsA) genes, conferring resistance to quinolones and nitrofurans, respectively, and the carriage of an emergent-specific plasmid, designated pESI. This self-transferred episome is a mosaic megaplasmid (∼280 kb), which increases bacterial tolerance to environmental mercury (mer operon) and oxidative stress, and provides further resistance to tetracycline, sulfamethoxazole and trimethoprim, most likely due to the presence of tetRA, sulI and dfrA genes respectively. Moreover, pESI carries the yersiniabactin siderophore system and two novel chaperone-usher fimbriae. In vitro studies established that pESI conjugation into a plasmidless S. Infantis strain results in superior biofilm formation, adhesion and invasion into avian and mammalian host cells. In vivo mouse infections demonstrated higher pathogenicity and increased intestinal inflammation caused by an S. Infantis strain harboring pESI compared with the plasmidless parental strain. Our results indicate that the presence of pESI that was found only in the emergent population of S. Infantis in Israel contributes significantly to antimicrobials tolerance and pathogenicity of its carrier. It is highly likely that pESI plays a key role in the successful spread of the emergent clone that replaced the local S. Infantis community in the short time of only 2-3 years.

Sequencing of bacterial genomes: principles and insights into pathogenesis and development of antibiotics

The impact of bacterial diseases on public health has become enormous, and is partly due to the increasing trend of antibiotic resistance displayed by bacterial pathogens. Sequencing of bacterial genomes has significantly improved our understanding about the biology of many bacterial pathogens as well as identification of novel antibiotic targets. Since the advent of genome sequencing two decades ago, about 1,800 bacterial genomes have been fully sequenced and these include important aetiological agents such as Streptococcus pneumoniae, Mycobacterium tuberculosis, Escherichia coli O157:H7, Vibrio cholerae, Clostridium difficile and Staphylococcus aureus. Very recently, there has been an explosion of bacterial genome data and is due to the development of next generation sequencing technologies, which are evolving so rapidly. Indeed, the field of microbial genomics is advancing at a very fast rate and it is difficult for researchers to be abreast with the new developments. This highlights the need for regular updates in microbial genomics through comprehensive reviews. This review paper seeks to provide an update on bacterial genome sequencing generally, and to analyze insights gained from sequencing in two areas, including bacterial pathogenesis and the development of antibiotics.

Virulence gene profiles and population genetic analysis for exploration of pathogenic serogroups of Shiga toxin-producing Escherichia coli

Infection with Shiga toxin (Stx)-producing Escherichia coli (STEC) is a serious public health concern, causing severe diarrhea and hemolytic-uremic syndrome. Patient symptoms are varied among STEC strains, possibly implying the presence of markers for STEC virulence other than Stx. To reveal the genotypic traits responsible for STEC virulence, we investigated 282 strains of various serogroups for the presence of 17 major virulence genes, i.e., stx1, stx2a, stx2c, stx2d, stx2e, stx2f, eae, tir, espB, espD, iha, saa, subA, ehxA, espP, katP, and stcE. Next, we examined the prevalence of virulence genes according to the seropathotypes in which serotypes were classified (seropathotypes A through E) based on the reported frequencies in human illness, as well as known associations with outbreaks and with severe disease. Our results demonstrate that the presence of both katP and stcE in STEC, in addition to the genes located in the locus of enterocyte effacement (LEE), including eae, tir, espB, and espD, may indicate the most pathogenic genotype of STEC. A population structure analysis of the profiles of virulence genes statistically supported the pathogenic genotype and, furthermore, revealed that there are serogroups with potentially higher pathogenicity than previously thought. Some strains in serogroups O26, O145, and O165 may have high virulence equivalent to that of serogroup O157. Several serogroups, including O14, O16, O45, O63, O74, 119, O128, and O untypeable, also may be potentially pathogenic, although rarely in humans.

Virulence gene profiling and pathogenicity characterization of non-typhoidal Salmonella accounted for invasive disease in humans

Human infection with non-typhoidal Salmonella serovars (NTS) infrequently causes invasive systemic disease and bacteremia. To understand better the nature of invasive NTS (iNTS), we studied the gene content and the pathogenicity of bacteremic strains from twelve serovars (Typhimurium, Enteritidis, Choleraesuis, Dublin, Virchow, Newport, Bredeney, Heidelberg, Montevideo, Schwarzengrund, 9,12:l,v:- and Hadar). Comparative genomic hybridization using a Salmonella enterica microarray revealed a core of 3233 genes present in all of the iNTS strains, which include the Salmonella pathogenicity islands 1-5, 9, 13, 14; five fimbrial operons (bcf, csg, stb, sth, sti); three colonization factors (misL, bapA, sinH); and the invasion gene, pagN. In the iNTS variable genome, we identified 16 novel genomic islets; various NTS virulence factors; and six typhoid-associated virulence genes (tcfA, cdtB, hlyE, taiA, STY1413, STY1360), displaying a wider distribution among NTS than was previously known. Characterization of the bacteremic strains in C3H/HeN mice showed clear differences in disease manifestation. Previously unreported characterization of serovars Schwarzengrund, 9,12:l,v:-, Bredeney and Virchow in the mouse model showed low ability to elicit systemic disease, but a profound and elongated shedding of serovars Schwarzengrund and 9,12:l,v:- (as well as Enteritidis and Heidelberg) due to chronic infection of the mouse. Phenotypic comparison in macrophages and epithelial cell lines demonstrated a remarkable intra-serovar variation, but also showed that S. Typhimurium bacteremic strains tend to present lower intracellular growth than gastroenteritis isolates. Collectively, our data demonstrated a common core of virulence genes, which might be required for invasive salmonellosis, but also an impressive degree of genetic and phenotypic heterogeneity, highlighting that bacteremia is a complex phenotype, which cannot be attributed merely to an enhanced invasion or intracellular growth of a particular strain.

Identification of accessory genome regions in poultry Clostridium perfringens isolates carrying the netB plasmid

Necrotic enteritis (NE) is an economically important disease of poultry caused by certain Clostridium perfringens type A strains. NE pathogenesis involves the NetB toxin, which is encoded on a large conjugative plasmid within a 42-kb pathogenicity locus. Recent multilocus sequence type (MLST) studies have identified two predominant NE-associated clonal groups, suggesting that host genes are also involved in NE pathogenesis. We used microarray comparative genomic hybridization (CGH) to assess the gene content of 54 poultry isolates from birds that were healthy or that suffered from NE. A total of 400 genes were variably present among the poultry isolates and nine nonpoultry strains, many of which had putative functions related to nutrient uptake and metabolism and cell wall and capsule biosynthesis. The variable genes were organized into 142 genomic regions, 49 of which contained genes significantly associated with netB-positive isolates. These regions included three previously identified NE-associated loci as well as several apparent fitness-related loci, such as a carbohydrate ABC transporter, a ferric-iron siderophore uptake system, and an adhesion locus. Additional loci were related to plasmid maintenance. Cluster analysis of the CGH data grouped all of the netB-positive poultry isolates into two major groups, separated according to two prevalent clonal groups based on MLST analysis. This study identifies chromosomal loci associated with netB-positive poultry strains, suggesting that the chromosomal background can confer a selective advantage to NE-causing strains, possibly through mechanisms involving iron acquisition, carbohydrate metabolism, and plasmid maintenance.

An evolutionary perspective on the kinome of malaria parasites

Malaria parasites belong to an ancient lineage that diverged very early from the main branch of eukaryotes. The approximately 90-member plasmodial kinome includes a majority of eukaryotic protein kinases that clearly cluster within the AGC, CMGC, TKL, CaMK and CK1 groups found in yeast, plants and mammals, testifying to the ancient ancestry of these families. However, several hundred millions years of independent evolution, and the specific pressures brought about by first a photosynthetic and then a parasitic lifestyle, led to the emergence of unique features in the plasmodial kinome. These include taxon-restricted kinase families, and unique peculiarities of individual enzymes even when they have homologues in other eukaryotes. Here, we merge essential aspects of all three malaria-related communications that were presented at the Evolution of Protein Phosphorylation meeting, and propose an integrated discussion of the specific features of the parasite's kinome and phosphoproteome.

Evolution of variation in presence and absence of genes in bacterial pathways

Background

Bacterial genomes exhibit a remarkable degree of variation in the presence and absence of genes, which probably extends to the level of individual pathways. This variation may be a consequence of the significant evolutionary role played by horizontal gene transfer, but might also be explained by the loss of genes through mutation. A challenge is to understand why there would be variation in gene presence within pathways if they confer a benefit only when complete.

Results

Here, we develop a mathematical model to study how variation in pathway content is produced by horizontal transfer, gene loss and partial exposure of a population to a novel environment.

Conclusions

We discuss the possibility that variation in gene presence acts as cryptic genetic variation on which selection acts when the appropriate environment occurs. We find that a high level of variation in gene presence can be readily explained by decay of the pathway through mutation when there is no longer exposure to the selective environment, or when selection becomes too weak to maintain the genes. In the context of pathway variation the role of horizontal gene transfer is probably the initial introduction of a complete novel pathway rather than in building up the variation in a genome without the pathway.

Polyphyly of non-bioluminescent Vibrio fischeri sharing a lux-locus deletion

This study reports the first description and molecular characterization of naturally occurring, non-bioluminescent strains of Vibrio fischeri. These 'dark' V. fischeri strains remained non-bioluminescent even after treatment with both autoinducer and aldehyde, substrate additions that typically maximize light production in dim strains of luminous bacteria. Surprisingly, the entire lux locus (eight genes) was absent in over 97% of these dark V. fischeri strains. Although these strains were all collected from a Massachusetts (USA) estuary in 2007, phylogenetic reconstructions allowed us to reject the hypothesis that these newly described non-bioluminescent strains exhibit monophyly within the V. fischeri clade. These dark strains exhibited a competitive disadvantage against native bioluminescent strains when colonizing the light organ of the model V. fischeri host, the Hawaiian bobtail squid Euprymna scolopes. Significantly, we believe that the data collected in this study may suggest the first observation of a functional, parallel locus-deletion event among independent lineages of a non-pathogenic bacterial species.

Structural and evolutionary divergence of eukaryotic protein kinases in Apicomplexa

Background

The Apicomplexa constitute an evolutionarily divergent phylum of protozoan pathogens responsible for widespread parasitic diseases such as malaria and toxoplasmosis. Many cellular functions in these medically important organisms are controlled by protein kinases, which have emerged as promising drug targets for parasitic diseases. However, an incomplete understanding of how apicomplexan kinases structurally and mechanistically differ from their host counterparts has hindered drug development efforts to target parasite kinases.

Results

We used the wealth of sequence data recently made available for 15 apicomplexan species to identify the kinome of each species and quantify the evolutionary constraints imposed on each family of apicomplexan kinases. Our analysis revealed lineage-specific adaptations in selected families, namely cyclin-dependent kinase (CDK), calcium-dependent protein kinase (CDPK) and CLK/LAMMER, which have been identified as important in the pathogenesis of these organisms. Bayesian analysis of selective constraints imposed on these families identified the sequence and structural features that most distinguish apicomplexan protein kinases from their homologs in model organisms and other eukaryotes. In particular, in a subfamily of CDKs orthologous to Plasmodium falciparum crk-5, the activation loop contains a novel PTxC motif which is absent from all CDKs outside Apicomplexa. Our analysis also suggests a convergent mode of regulation in a subset of apicomplexan CDPKs and mammalian MAPKs involving a commonly conserved arginine in the αC helix. In all recognized apicomplexan CLKs, we find a set of co-conserved residues involved in substrate recognition and docking that are distinct from metazoan CLKs.

Conclusions

We pinpoint key conserved residues that can be predicted to mediate functional differences from eukaryotic homologs in three identified kinase families. We discuss the structural, functional and evolutionary implications of these lineage-specific variations and propose specific hypotheses for experimental investigation. The apicomplexan-specific kinase features reported in this study can be used in the design of selective kinase inhibitors.

The influence of rickettsiologists on post-modern microbiology

Many of the definitions in microbiology are currently false. We have reviewed the great denominations of microbiology and attempted to free microorganisms from the theories of the twentieth century. The presence of compartmentation and a nucleoid in Planctomycetes clearly calls into question the accuracy of the definitions of eukaryotes and prokaryotes. Archaea are viewed as prokaryotes resembling bacteria. However, the name archaea, suggesting an archaic origin of lifestyle, is inconsistent with the lifestyle of this family. Viruses are defined as small, filterable infectious agents, but giant viruses challenge the size criteria used for the definition of a virus. Pathogenicity does not require the acquisition of virulence factors (except for toxins), and in many cases, gene loss is significantly inked to the emergence of virulence. Species classification based on 16S rRNA is useless for taxonomic purposes of human pathogens, as a 2% divergence would classify all Rickettsiae within the same species and would not identify bacteria specialized for mammal infection. The use of metagenomics helps us to understand evolution and physiology by elucidating the structure, function, and interactions of the major microbial communities, but it neglects the minority populations. Finally, Darwin’s descent with modification theory, as represented by the tree of life, no longer matches our current genomic knowledge because genomics has revealed the occurrence of de novo-created genes and the mosaic structure of genomes, the Rhizome of life is therefore more appropriate.

Targeted enrichment of ancient pathogens yielding the pPCP1 plasmid of Yersinia pestis from victims of the Black Death

Although investigations of medieval plague victims have identified Yersinia pestis as the putative etiologic agent of the pandemic, methodological limitations have prevented large-scale genomic investigations to evaluate changes in the pathogen's virulence over time. We screened over 100 skeletal remains from Black Death victims of the East Smithfield mass burial site (1348-1350, London, England). Recent methods of DNA enrichment coupled with high-throughput DNA sequencing subsequently permitted reconstruction of ten full human mitochondrial genomes (16 kb each) and the full pPCP1 (9.6 kb) virulence-associated plasmid at high coverage. Comparisons of molecular damage profiles between endogenous human and Y. pestis DNA confirmed its authenticity as an ancient pathogen, thus representing the longest contiguous genomic sequence for an ancient pathogen to date. Comparison of our reconstructed plasmid against modern Y. pestis shows identity with several isolates matching the Medievalis biovar; however, our chromosomal sequences indicate the victims were infected with a Y. pestis variant that has not been previously reported. Our data reveal that the Black Death in medieval Europe was caused by a variant of Y. pestis that may no longer exist, and genetic data carried on its pPCP1 plasmid were not responsible for the purported epidemiological differences between ancient and modern forms of Y. pestis infections.

Endodontic microbiology

Root canal therapy has been practiced ever since 1928 and the success rate has tremendously increased over the years owing to various advancements in the field. One main reason is the complete understanding of the microbiology involved in the endodontic pathology. This has helped us to modify the conventional treatment plans and effectively combat the microorganisms. Now, studies are aiming to explore the characteristics of the "most" resistant organism and the methods to eliminate them. This article gives an insight of the microbiology involved in endodontic pathology and discusses its role in our treatment procedure. Information from original reviews listed in PubMed, published from 1995 to 2010, has been mainly included in this review.

The evolution of virulence in non-o157 shiga toxin-producing Escherichia coli

Shiga toxin-producing Escherichia coli (STEC) are zoonotic foodborne and waterborne pathogens that are a serious public health concern because they cause outbreaks and the potentially fatal hemolytic uremic syndrome (HUS). The most common STEC serotype associated with human disease is O157:H7, but there is a growing recognition of over 100 non-O157 serotypes that also may result in human illness. Some of these non-O157 STEC strains cause outbreaks and severe disease such as HUS and hemorrhagic colitis, whereas others are associated with only mild diarrhea or with no human disease at all. The relative scarceness of whole genome sequence data for non-O157 STEC has limited the scientific discovery into the genetic basis of these differences in virulence. Uncovering the scope of genetic diversity and phylogeny of the non-O157 STEC through targeting sequencing of clinically relevant isolates will offer new biological insight to the pathogenic behavior of these emerging pathogens. These approaches would also enable molecular risk assessment strategies to rapidly identify and respond to emerging non-O157 STEC that pose a serious public health risk to humans.

Accurate localization of the integration sites of two genomic islands at single-nucleotide resolution in the genome of Bacillus cereus ATCC 10987

We have identified two genomic islands, that is, BCEGI-1 and BCEGI-2, in the genome of Bacillus cereus ATCC 10987, based on comparative analysis with Bacillus cereus ATCC 14579. Furthermore, by using the cumulative GC profile and performing homology searches between the two genomes, the integration sites of the two genomic islands were determined at single-nucleotide resolution. BCEGI-1 is integrated between 159705 bp and 198000 bp, whereas BCEGI-2 is integrated between the end of ORF BCE4594 and the start of the intergenic sequence immediately following BCE4626, that is, from 4256803 bp to 4285534 bp. BCEGI-1 harbors two bacterial Tn7 transposons, which have two sets of genes encoding TnsA, B, C, and D. It is generally believed that unlike the TnsABC+E pathway, the TnsABC+D pathway would only promote vertical transmission to daughter cells. The evidence presented in this paper, however, suggests a role of the TnsABC+D pathway in the horizontal transfer of some genomic islands.

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.

Evolutionary diversification of an ancient gene family (rhs) through C-terminal displacement

Background

Rhs genes are prominent features of bacterial genomes that have previously been implicated in genomic rearrangements in E. coli. By comparing rhs repertoires across the Enterobacteriaceae, this study provides a robust explanation of rhs diversification and evolution, and a mechanistic model of how rhs diversity is gained and lost.

Results

Rhs genes are ubiquitous and comprise six structurally distinct lineages within the Enterobacteriaceae. There is considerable intergenomic variation in rhs repertoire; for instance, in Salmonella enterica, rhs are restricted to mobile elements, while in Escherichia coli one rhs lineage has diversified through transposition as older lineages have been deleted. Overall, comparative genomics reveals frequent, independent gene gains and losses, as well as occasional lateral gene transfer, in different genera. Furthermore, we demonstrate that Rhs 'core' domains and variable C-termini are evolutionarily decoupled, and propose that rhs diversity is driven by homologous recombination with circular intermediates. Existing C-termini are displaced by laterally acquired alternatives, creating long arrays of dissociated 'tips' that characterize the appearance of rhs loci.

Conclusion

Rhs repertoires are highly dynamic among Enterobacterial genomes, due to repeated gene gains and losses. In contrast, the primary structures of Rhs genes are evolutionarily conserved, indicating that rhs sequence diversity is driven, not by rapid mutation, but by the relatively slow evolution of novel core/tip combinations. Hence, we predict that a large pool of dissociated rhs C-terminal tips exists episomally and these are potentially transmitted across taxonomic boundaries.

Gene transfer and diversification of microbial eukaryotes

The importance of lateral gene transfer in genome evolution of microbial eukaryotes is slowly being appreciated. Acquisitions of genes have led to metabolic adaptation in diverse eukaryotic lineages. In most cases the metabolic genes have originated from prokaryotes, often followed by sequential transfers between eukaryotes. However, the knowledge of gene transfer in eukaryotes is still mainly based on anecdotal evidence. Some of the observed patterns may be biases in experimental approaches and sequence databases rather than evolutionary trends. Rigorous systematic studies of gene acquisitions that allow for the possibility of exchanges of all categories of genes from all sources are needed to get a more objective view of gene transfer in eukaryote evolution. It may be that the role of gene transfer in the diversification process of microbial eukaryotes currently is underestimated.

Bordetella pertussis and vaccination: the persistence of a genetically monomorphic pathogen

Before childhood vaccination was introduced in the 1950s, pertussis or whooping cough was a major cause of infant death worldwide. Widespread vaccination of children was successful in significantly reducing morbidity and mortality. However, despite vaccination, pertussis has persisted and, in the 1990s, resurged in a number of countries with highly vaccinated populations. Indeed, pertussis has become the most prevalent vaccine-preventable disease in developed countries with estimated infection frequencies of 1-6%. Recently vaccinated children are well protected against pertussis disease and its increase is mainly seen in adolescents and adults in which disease symptoms are often mild. The etiologic agent of pertussis, Bordetella pertussis, is extremely monomorphic and its ability to persist in the face of intensive vaccination is intriguing. Numerous studies have shown that B. pertussis populations changed after the introduction of vaccination suggesting adaptation. These adaptations did not involve the acquisition of novel genes but small genetic changes, mainly SNPs, and occurred in successive steps in a period of 40 years. The earliest adaptations resulted in antigenic divergence with vaccine strains. More recently, strains emerged with increased pertussis toxin (Ptx) production. Here I argue that the resurgence of pertussis is the compound effect of pathogen adaptation and waning immunity. I propose that the removal by vaccination of naïve infants as the major source for transmission was the crucial event which has driven the changes in B. pertussis populations. This has selected for strains which are more efficiently transmitted by primed hosts in which immunity has waned. The adaptation of B. pertussis to primed hosts involved delaying an effective immune response by antigenic divergence with vaccine strains and by increasing immune suppression through higher levels of Ptx production. Higher levels of Ptx may also benefit transmission by enhancing clinical symptoms. The study of B. pertussis populations has not only increased our understanding of pathogen evolution, but also suggests way to improve pertussis vaccines, underlining the public health significance of population-based studies of pathogens.

Comparative genomic analyses of the human fungal pathogens Coccidioides and their relatives

While most Ascomycetes tend to associate principally with plants, the dimorphic fungi Coccidioides immitis and Coccidioides posadasii are primary pathogens of immunocompetent mammals, including humans. Infection results from environmental exposure to Coccidiodies, which is believed to grow as a soil saprophyte in arid deserts. To investigate hypotheses about the life history and evolution of Coccidioides, the genomes of several Onygenales, including C. immitis and C. posadasii; a close, nonpathogenic relative, Uncinocarpus reesii; and a more diverged pathogenic fungus, Histoplasma capsulatum, were sequenced and compared with those of 13 more distantly related Ascomycetes. This analysis identified increases and decreases in gene family size associated with a host/substrate shift from plants to animals in the Onygenales. In addition, comparison among Onygenales genomes revealed evolutionary changes in Coccidioides that may underlie its infectious phenotype, the identification of which may facilitate improved treatment and prevention of coccidioidomycosis. Overall, the results suggest that Coccidioides species are not soil saprophytes, but that they have evolved to remain associated with their dead animal hosts in soil, and that Coccidioides metabolism genes, membrane-related proteins, and putatively antigenic compounds have evolved in response to interaction with an animal host.

Genetic identity and differential gene expression between Trichomonas vaginalis and Trichomonas tenax

BACKGROUND

Trichomonas vaginalis is a human urogenital pathogen responsible for trichomonosis, the number-one, non-viral sexually transmitted disease (STD) worldwide, while T. tenax is a commensal of the human oral cavity, found particularly in patients with poor oral hygiene and advanced periodontal disease. The extent of genetic identity between T. vaginalis and its oral commensal counterpart is unknown.

RESULTS

Genes that were differentially expressed in T. vaginalis were identified by screening three independent subtraction cDNA libraries enriched for T. vaginalis genes. The same thirty randomly selected cDNA clones encoding for proteins with specific functions associated with colonization were identified from each of the subtraction cDNA libraries. In addition, a T. vaginalis cDNA expression library was screened with patient sera that was first pre-adsorbed with an extract of T. tenax antigens, and seven specific cDNA clones were identified from this cDNA library. Interestingly, some of the clones identified by the subtraction cDNA screening were also obtained from the cDNA expression library with the pre-adsorbed sera. Moreover and noteworthy, clones identified by both the procedures were found to be up-regulated in expression in T. vaginalis upon contact with vaginal epithelial cells, suggesting a role for these gene products in host colonization. Semi-quantitative RT-PCR analysis of select clones showed that the genes were not unique to T. vaginalis and that these genes were also present in T. tenax, albeit at very low levels of expression.

CONCLUSION

These results suggest that T. vaginalis and T. tenax have remarkable genetic identity and that T. vaginalis has higher levels of gene expression when compared to that of T. tenax. The data may suggest that T. tenax could be a variant of T. vaginalis.

Pathogenic adaptation of intracellular bacteria by rewiring a cis-regulatory input function

The acquisition of DNA by horizontal gene transfer enables bacteria to adapt to previously unexploited ecological niches. Although horizontal gene transfer and mutation of protein-coding sequences are well-recognized forms of pathogen evolution, the evolutionary significance of cis-regulatory mutations in creating phenotypic diversity through altered transcriptional outputs is not known. We show the significance of regulatory mutation for pathogen evolution by mapping and then rewiring a cis-regulatory module controlling a gene required for murine typhoid. Acquisition of a binding site for the Salmonella pathogenicity island-2 regulator, SsrB, enabled the srfN gene, ancestral to the Salmonella genus, to play a role in pathoadaptation of S. typhimurium to a host animal. We identified the evolved cis-regulatory module and quantified the fitness gain that this regulatory output accrues for the bacterium using competitive infections of host animals. Our findings highlight a mechanism of pathogen evolution involving regulatory mutation that is selected because of the fitness advantage the new regulatory output provides the incipient clones.

Horizontal gene transfer in eukaryotic parasites: a case study of Entamoeba histolytica and Trichomonas vaginalis

Over the past few years it has become apparent that horizontal gene transfer (HGT) has played an important role in the evolution of pathogenic prokaryotes. What is less clear is the exact role that HGT has played in shaping the metabolism of eukaryotic organisms. The main problems are the reliable inference of HGT on a genomic scale as well as the functional assignment of genes in these poorly studied organisms. We have screened the completed genomes of the protists Entamoeba histolytica and Trichomonas vaginalis for cases of HGT from prokaryotes. Using a fast primary screen followed by a conservative phylogenetic approach, we found 68 and 153 recent cases of HGT in the respective organisms. The majority of transferred genes that fall into functional categories code for enzymes involved in metabolism. We found a broad range of prokaryotic lineages represented among the donors, but organisms that share similar environmental niches with E. histolytica and T. vaginalis, such as the gut and the vaginal mucosa, dominate.

Comparative genomic profiling of Dutch clinical Bordetella pertussis isolates using DNA microarrays: identification of genes absent from epidemic strains

Background

Whooping cough caused by Bordetella pertussis in humans, is re-emerging in many countries despite vaccination. Several studies have shown that significant shifts have occurred in the B. pertussis population resulting in antigenic divergence between vaccine strains and circulating strains and suggesting pathogen adaptation. In the Netherlands, the resurgence of pertussis is associated with the rise of B. pertussis strains with an altered promoter region for pertussis toxin (ptxP3).

Results

We used Multi-Locus Sequence Typing (MLST), Multiple-Locus Variable Number of Tandem Repeat Analysis (MLVA) and microarray-based comparative genomic hybridization (CGH) to characterize the ptxP3 strains associated with the Dutch epidemic. For CGH analysis, we developed an oligonucleotide (70-mers) microarray consisting of 3,581 oligonucleotides representing 94% of the gene repertoire of the B. pertussis strain Tohama I. Nine different MLST profiles and 38 different MLVA types were found in the period 1993 to 2004. Forty-three Dutch clinical isolates were analyzed with CGH, 98 genes were found to be absent in at least one of the B. pertussis strains tested, these genes were clustered in 8 distinct regions of difference.

Conclusion

The presented MLST, MLVA and CGH-analysis identified distinctive characteristics of ptxP3 B. pertussis strains -the most prominent of which was a genomic deletion removing about 23,000 bp. We propose a model for the emergence of ptxP3 strains.

Completely sequenced genomes of pathogenic bacteria: a review

Six out of ten completely sequenced bacterial genomes are pathogenic or opportunistic bacteria. The genome sequence of at least one strain of all the principal pathogenic bacteria will soon be available. This information should enable us to identify genes that encode virulence factors. As these genes are potential targets for drugs and vaccines, their identification should have considerable repercussions on prevention, diagnosis, and treatment of the main bacterial infectious diseases. Comparison of genome sequences of several strains of the same species should allow identification of the genetic clues responsible for the differing behavior of related bacterial pathogens. This article reviews the genomes from pathogenic bacteria that have been or are currently being sequenced, describes the main tasks to be accomplished after a genome sequence becomes available, and discusses the benefits of having the genome sequence of bacterial pathogens.