The evolutionary origin of Xanthomonadales genomes and the nature of the horizontal gene transfer process

Antibiotic Resistance in Plant Pathogenic Bacteria: Recent Data and Environmental Impact of Unchecked Use and the Potential of Biocontrol Agents as an Eco-Friendly Alternative

The economic impact of phytopathogenic bacteria on agriculture is staggering, costing billions of US dollars globally. Pseudomonas syringae is the top most phytopathogenic bacteria, having more than 60 pathovars, which cause bacteria speck in tomatoes, halo blight in beans, and so on. Although antibiotics or a combination of antibiotics are used to manage infectious diseases in plants, they are employed far less in agriculture compared to human and animal populations. Moreover, the majority of antibiotics used in plants are immediately washed away, leading to environmental damage to ecosystems and food chains. Due to the serious risk of antibiotic resistance (AR) and the potential for environmental contamination with antibiotic residues and resistance genes, the use of unchecked antibiotics against phytopathogenic bacteria is not advisable. Despite the significant concern regarding AR in the world today, there are inadequate and outdated data on the AR of phytopathogenic bacteria. This review presents recent AR data on plant pathogenic bacteria (PPB), along with their environmental impact. In light of these findings, we suggest the use of biocontrol agents as a sustainable, eco-friendly, and effective alternative to controlling phytopathogenic bacteria.

Comparative genomics of the black rot pathogen Xanthomonas campestris pv. campestris and non-pathogenic co-inhabitant Xanthomonas melonis from Trinidad reveal unique pathogenicity determinants and secretion system profiles

Black-rot disease caused by the phytopathogen Xanthomonas campestris pv. campestris (Xcc) continues to have considerable impacts on the productivity of cruciferous crops in Trinidad and Tobago and the wider Caribbean region. While the widespread occurrence of resistance of Xcc against bactericidal agrochemicals can contribute to the high disease burdens, the role of virulence and pathogenicity features of local strains on disease prevalence and severity has not been investigated yet. In the present study, a comparative genomic analysis was performed on 6 pathogenic Xcc and 4 co-isolated non-pathogenic Xanthomonas melonis (Xmel) strains from diseased crucifer plants grown in fields with heavy chemical use in Trinidad. Native isolates were grouped into two known and four newly assigned ribosomal sequence types (rST). Mobile genetic elements were identified which belonged to the IS3, IS5 family, Tn3 transposon, resolvases, and tra T4SS gene clusters. Additionally, exogenous plasmid derived sequences with origins from other bacterial species were characterised. Although several instances of genomic rearrangements were observed, native Xcc and Xmel isolates shared a significant level of structural homology with reference genomes, Xcc ATCC 33913 and Xmel CFBP4644, respectively. Complete T1SS hlyDB, T2SS, T4SS vir and T5SS xadA, yapH and estA gene clusters were identified in both species. Only Xmel strains contained a complete T6SS but no T3SS. Both species contained a complex repertoire of extracellular cell wall degrading enzymes. Native Xcc strains contained 37 T3SS and effector genes but a variable and unique profile of 8 avr, 4 xop and 1 hpa genes. Interestingly, Xmel strains contained several T3SS effectors with low similarity to references including avrXccA1 (~89%), hrpG (~73%), hrpX (~90%) and xopAZ (~87%). Furthermore, only Xmel genomes contained a CRISPR-Cas I-F array, but no lipopolysaccharide wxc gene cluster. Xmel strains were confirmed to be non-pathogenic by pathogenicity assays. The results of this study will be useful to guide future research into virulence mechanisms, agrochemical resistance, pathogenomics and the potential role of the co-isolated non-pathogenic Xanthomonas strains on Xcc infections.

Phylogenetic Analyses of Xanthomonads Causing Bacterial Leaf Spot of Tomato and Pepper: Xanthomonas euvesicatoria Revealed Homologous Populations Despite Distant Geographical Distribution

Bacterial leaf spot of tomato and pepper (BLS), an economically important bacterial disease caused by four species of Xanthomonas (X. euvesicatoria (Xe), X. vesicatoria (Xv), X. gardneri (Xg), and X. perforans (Xp)), is a global problem and can cause over 50% crop loss under unfavorable conditions. Among the four species, Xe and Xv are prevalent worldwide. Characterization of the pathogens is crucial for disease management and regulatory purposes. In this study, we performed a multilocus sequence analysis (MLSA) with six genes (hrcN, dnaA gyrB, gapA, pdg, and hmbs) on BLS strains. Other Xanthomonas species were included to determine phylogenetic relationships within and among the tested strains. Four BLS species comprising 76 strains from different serological groups and diverse geographical locations were resolved into three major clades. BLS xanthomonads formed distinct clusters in the phylogenetic analyses. Three other xanthomonads, including X. albilineans, X. sacchari, and X. translucens pv. undolusa revealed less than 85%, 88%, and 89% average nucleotide identity (ANI), respectively, with the other species of Xanthomonas included in this study. Both antibody and MLSA data showed that Xv was clearly separated from Xe and that the latter strains were remarkably clonal, even though they originated from distant geographical locations. The Xe strains formed two separate phylogenetic groups; Xe group A1 consisted only of tomato strains, whereas Xe group A2 included strains from pepper and tomato. In contrast, the Xv group showed greater heterogeneity. Some Xv strains from South America were closely related to strains from California, while others grouped closer to a strain from Indiana and more distantly to a strain from Hawaii. Using this information molecular tests can now be devised to track distribution of clonal populations that may be introduced into new geographic areas through seeds and other infected plant materials.

Bacteria-Killing Type IV Secretion Systems

Bacteria have been constantly competing for nutrients and space for billions of years. During this time, they have evolved many different molecular mechanisms by which to secrete proteinaceous effectors in order to manipulate and often kill rival bacterial and eukaryotic cells. These processes often employ large multimeric transmembrane nanomachines that have been classified as types I-IX secretion systems. One of the most evolutionarily versatile are the Type IV secretion systems (T4SSs), which have been shown to be able to secrete macromolecules directly into both eukaryotic and prokaryotic cells. Until recently, examples of T4SS-mediated macromolecule transfer from one bacterium to another was restricted to protein-DNA complexes during bacterial conjugation. This view changed when it was shown by our group that many Xanthomonas species carry a T4SS that is specialized to transfer toxic bacterial effectors into rival bacterial cells, resulting in cell death. This review will focus on this special subtype of T4SS by describing its distinguishing features, similar systems in other proteobacterial genomes, and the nature of the effectors secreted by these systems and their cognate inhibitors.

Cyanobacterial phylogenetic analysis based on phylogenomics approaches render evolutionary diversification and adaptation: an overview of representative orders

Phylogenetic studies based on a definite set of marker genes usually reconstruct evolutionary relationships among the prokaryotic species. Based on specific target sequences, such studies represent variations and allow identification of similarities or dissimilarities in organisms. With the advent of completely sequenced genomes and accumulation of information on whole prokaryotic genomes, phylogenetic reconstructions should be considered more reliable if they are ideally based on entire genomes to resolve phylogenetic interest. We applied phylogenomics approaches taking into account completely sequenced cyanobacterial genomes to reconstruct underlying species that represented major taxonomic classes and belonged to distinctly different habitats (freshwater, marine, soils, and rocks). We did not rely on describing phylogeny of all representative class of cyanobacterial species on the basis of only ribosomal gene, 16S rDNA gene. In contrast, we analyzed combined molecular marker and phylogenomics approaches (genome alignment, gene content and gene order, composition vector and protein domain content) for accurately inferring phylogenetic relationship of species. We have shown that this approach reflects the impact of evolution on the organisms and considers connects with the ecological adaptation in cyanobacteria in different habitats. Analysis revealed that the members from marine habitat occupy different profile than those from freshwater. Impact of GC content and genomic repetitiveness over the diversification of cyanobacterial species and their possible role in adaptation was also reflected. Members occupying similar habitats cover more evolutionary distance together and also evolve various strategies for adaptation and survival either through genomic repetitiveness or preferences for genes of particular functions or modified GC content. Genomes undergo different changes for their adaptation in diverse habitats.

A highly specific tool for identification of Xanthomonas vasicola pv. musacearum based on five Xvm-specific coding sequences

Xanthomonas vasicola pv. musacearum (Xvm) is a bacterial pathogen responsible for the economically important Xanthomonas wilt disease on banana and enset crops in Sub-Saharan Africa. Given that the symptoms are similar to those of other diseases, molecular diagnosis is essential to unambiguously identify this pathogen and distinguish it from closely related strains not pathogenic on these hosts. Currently, Xvm identification is based on polymerase chain reaction (PCR) with GspDm primers, targeting the gene encoding general secretory protein D. Experimental results and examination of genomic sequences revealed poor specificity of the GspDm PCR. Here, we present and validate five new Xvm-specific primers amplifying only Xvm strains.

Tremblaya phenacola PPER: an evolutionary beta-gammaproteobacterium collage

Many insects rely on bacterial endosymbionts to obtain nutrients that are scarce in their highly specialized diets. The most surprising example corresponds to the endosymbiotic system found in mealybugs from subfamily Pseudococcinae in which two bacteria, the betaproteobacterium 'Candidatus Tremblaya princeps' and a gammaproteobacterium, maintain a nested endosymbiotic consortium. In the sister subfamily Phenacoccinae, however, a single beta-endosymbiont, 'Candidatus Tremblaya phenacola', has been described. In a previous study, we detected a trpB gene of gammaproteobacterial origin in 'Ca. Tremblaya phenacola' from two Phenacoccus species, apparently indicating an unusual case of horizontal gene transfer (HGT) in a bacterial endosymbiont. What we found by sequencing the genome of 'Ca. Tremblaya phenacola' PPER, single endosymbiont of Phenacoccus peruvianus, goes beyond a HGT phenomenon. It rather represents a genome fusion between a beta and a gammaproteobacterium, followed by massive rearrangements and loss of redundant genes, leading to an unprecedented evolutionary collage. Mediated by the presence of several repeated sequences, there are many possible genome arrangements, and different subgenomic sequences might coexist within the same population.

Et tu, Brute? Not Even Intracellular Mutualistic Symbionts Escape Horizontal Gene Transfer

Many insect species maintain mutualistic relationships with endosymbiotic bacteria. In contrast to their free-living relatives, horizontal gene transfer (HGT) has traditionally been considered rare in long-term endosymbionts. Nevertheless, meta-omics exploration of certain symbiotic models has unveiled an increasing number of bacteria-bacteria and bacteria-host genetic transfers. The abundance and function of transferred loci suggest that HGT might play a major role in the evolution of the corresponding consortia, enhancing their adaptive value or buffering detrimental effects derived from the reductive evolution of endosymbionts’ genomes. Here, we comprehensively review the HGT cases recorded to date in insect-bacteria mutualistic consortia, and discuss their impact on the evolutionary success of these associations.

Culturable aerobic and facultative bacteria from the gut of the polyphagic dung beetle Thorectes lusitanicus

Unlike other dung beetles, the Iberian geotrupid, Thorectes lusitanicus, exhibits polyphagous behavior; for example, it is able to eat acorns, fungi, fruits, and carrion in addition to the dung of different mammals. This adaptation to digest a wider diet has physiological and developmental advantages and requires key changes in the composition and diversity of the beetle's gut microbiota. In this study, we isolated aerobic, facultative anaerobic, and aerotolerant microbiota amenable to grow in culture from the gut contents of T. lusitanicus and resolved isolate identity to the species level by sequencing 16S rRNA gene fragments. Using BLAST similarity searches and maximum likelihood phylogenetic analyses, we were able to reveal that the analyzed fraction (culturable, aerobic, facultative anaerobic, and aerotolerant) of beetle gut microbiota is dominated by the phyla Proteobacteria, Firmicutes, and Actinobacteria. Among Proteobacteria, members of the order Enterobacteriales (Gammaproteobacteria) were the most abundant. The main functions associated with the bacteria found in the gut of T. lusitanicus would likely include nitrogen fixation, denitrification, detoxification, and diverse defensive roles against pathogens.

Examining marginal sequence similarities between bacterial type III secretion system components and Trypanosoma cruzi surface proteins: horizontal gene transfer or convergent evolution?

The cell invasion mechanism of Trypanosoma cruzi has similarities with some intracellular bacterial taxa especially regarding calcium mobilization. This mechanism is not observed in other trypanosomatids, suggesting that the molecules involved in this type of cell invasion were a product of (1) acquisition by horizontal gene transfer (HGT); (2) secondary loss in the other trypanosomatid lineages of the mechanism inherited since the bifurcation Bacteria-Neomura (1.9 billion to 900 million years ago); or (3) de novo evolution from non-homologous proteins via convergent evolution. Similar to T. cruzi, several bacterial genera require increased host cell cytosolic calcium for intracellular invasion. Among intracellular bacteria, the mechanism of host cell invasion of genus Salmonella is the most similar to T. cruzi. The invasion of Salmonella occurs by contact with the host's cell surface and is mediated by the type III secretion system (T3SS) that promotes the contact-dependent translocation of effector proteins directly into host's cell cytoplasm. Here we provide evidence of distant sequence similarities and structurally conserved domains between T. cruzi and Salmonella spp T3SS proteins. Exhaustive database searches were directed to a wide range of intracellular bacteria and trypanosomatids, exploring sequence patterns for comparison of structural similarities and Bayesian phylogenies. Based on our data we hypothesize that T. cruzi acquired genes for calcium mobilization mediated invasion by ancient HGT from ancestral Salmonella lineages.

Evolutionary history of the plant pathogenic bacterium Xanthomonas axonopodis

Deciphering mechanisms shaping bacterial diversity should help to build tools to predict the emergence of infectious diseases. Xanthomonads are plant pathogenic bacteria found worldwide. Xanthomonas axonopodis is a genetically heterogeneous species clustering, into six groups, strains that are collectively pathogenic on a large number of plants. However, each strain displays a narrow host range. We address the question of the nature of the evolutionary processes--geographical and ecological speciation--that shaped this diversity. We assembled a large collection of X. axonopodis strains that were isolated over a long period, over continents, and from various hosts. Based on the sequence analysis of seven housekeeping genes, we found that recombination occurred as frequently as point mutation in the evolutionary history of X. axonopodis. However, the impact of recombination was about three times greater than the impact of mutation on the diversity observed in the whole dataset. We then reconstructed the clonal genealogy of the strains using coalescent and genealogy approaches and we studied the diversification of the pathogen using a model of divergence with migration. The suggested scenario involves a first step of generalist diversification that spanned over the last 25,000 years. A second step of ecology-driven specialization occurred during the past two centuries. Eventually, secondary contacts between host-specialized strains probably occurred as a result of agricultural development and intensification, allowing genetic exchanges of virulence-associated genes. These transfers may have favored the emergence of novel pathotypes. Finally, we argue that the largest ecological entity within X. axonopodis is the pathovar.

Phylogenomics and molecular signatures for species from the plant pathogen-containing order xanthomonadales

The species from the order Xanthomonadales, which harbors many important plant pathogens and some human pathogens, are currently distinguished primarily on the basis of their branching in the 16S rRNA tree. No molecular or biochemical characteristic is known that is specific for these bacteria. Phylogenetic and comparative analyses were conducted on 26 sequenced Xanthomonadales genomes to delineate their branching order and to identify molecular signatures consisting of conserved signature indels (CSIs) in protein sequences that are specific for these bacteria. In a phylogenetic tree based upon sequences for 28 proteins, Xanthomonadales species formed a strongly supported clade with Rhodanobacter sp. 2APBS1 as its deepest branch. Comparative analyses of protein sequences have identified 13 CSIs in widely distributed proteins such as GlnRS, TypA, MscL, LysRS, LipA, Tgt, LpxA, TolQ, ParE, PolA and TyrB that are unique to all species/strains from this order, but not found in any other bacteria. Fifteen additional CSIs in proteins (viz. CoxD, DnaE, PolA, SucA, AsnB, RecA, PyrG, LigA, MutS and TrmD) are uniquely shared by different Xanthomonadales except Rhodanobacter and in a few cases by Pseudoxanthomonas species, providing further support for the deep branching of these two genera. Five other CSIs are commonly shared by Xanthomonadales and 1–3 species from the orders Chromatiales, Methylococcales and Cardiobacteriales suggesting that these deep branching orders of Gammaproteobacteria might be specifically related. Lastly, 7 CSIs in ValRS, CarB, PyrE, GlyS, RnhB, MinD and X001065 are commonly shared by Xanthomonadales and a limited number of Beta- or Gamma-proteobacteria. Our analysis indicates that these CSIs have likely originated independently and they are not due to lateral gene transfers. The Xanthomonadales-specific CSIs reported here provide novel molecular markers for the identification of these important plant and human pathogens and also as potential targets for development of drugs/agents that specifically target these bacteria.

Multilocus sequence analysis and type III effector repertoire mining provide new insights into the evolutionary history and virulence of Xanthomonas oryzae

Multilocus sequence analysis (MLSA) and type III effector (T3E) repertoire mining were performed to gain new insights into the genetic relatedness of Xanthomonas oryzae pv. oryzae (Xoo) and Xanthomonas oryzae pv. oryzicola (Xoc), two major bacterial pathogens of rice. Based on a collection of 45 African and Asian strains, we first sequenced and analysed three housekeeping genes by MLSA, Bayesian clustering and a median-joining network approach. Second, we investigated the distribution of 32 T3E genes, which are known to be major virulence factors of plant pathogenic bacteria, in all selected strains, by polymerase chain reaction and dot-blot hybridization methods. The diversity observed within housekeeping genes, as well as within T3E repertoires, clearly showed that both pathogens belong to closely related, but distinct, phylogenetic groups. Interestingly, these evolutionary groups are differentiated according to the geographical origin of the strains, suggesting that populations of Xoo and Xoc might be endemic in Africa and Asia, and thus have evolved separately. We further revealed that T3E gene repertoires of both pathogens comprise core and variable gene suites that probably have distinct roles in pathogenicity and different evolutionary histories. In this study, we carried out a functional analysis of xopO, a differential T3E gene between Xoo and Xoc, to determine the involvement of this gene in tissue specificity. Altogether, our data contribute to a better understanding of the evolutionary history of Xoo and Xoc in Africa and Asia, and provide clues for functional studies aiming to understand the virulence, host and tissue specificity of both rice pathogens.

ClaMS: A Classifier for Metagenomic Sequences

ClaMS - "Classifier for Metagenomic Sequences" - is a Java application for binning assembled contigs in metagenomes using user-specified training sets and initial parameters. Since ClaMS trains on sequence composition-based genomic signatures, it is much faster than binning tools that rely on alignments to homologs; ClaMS can bin ~20,000 sequences in 3 minutes on a laptop with a 2.4 GH× Intel Core 2 Duo processor and 2 GB RAM. ClaMS is meant to be a desktop application for biologists and can be run on any machine under any Operating System on which the Java Runtime Environment can be installed.

Pathogenomics of Xanthomonas: understanding bacterium-plant interactions

Xanthomonas is a large genus of Gram-negative bacteria that cause disease in hundreds of plant hosts, including many economically important crops. Pathogenic species and pathovars within species show a high degree of host plant specificity and many exhibit tissue specificity, invading either the vascular system or the mesophyll tissue of the host. In this Review, we discuss the insights that functional and comparative genomic studies are providing into the adaptation of this group of bacteria to exploit the extraordinary diversity of plant hosts and different host tissues.

Evolutionary placement of Xanthomonadales based on conserved protein signature sequences

Xanthomonadales comprises one of the largest phytopathogenic bacterial groups, and is currently classified within the gamma-proteobacteria. However, the phylogenetic placement of this group is not clearly resolved, and the results of different studies contradict one another. In this work, the evolutionary position of Xanthomonadales was determined by analyzing the presence of shared insertions and deletions (INDELs) in highly conserved proteins. Several distinctive insertions found in most of the members of the gamma-proteobacteria are absent in Xanthomonadales and groups such as Legionelalles, Chromatiales, Methylococcales, Thiotrichales and Cardiobacteriales. These INDELs were most likely introduced after the branching of Xanthomonadales from most of the gamma-proteobacteria and provide evidence for the phylogenetic placement of the early gamma-proteobacteria. Moreover, other proteins contain insertions exclusive to the Xanthomonadales order, confirming that this is a monophyletic group and provide important specific genetic markers. Thus, the data presented clearly support the Xanthomonadales group as an independent subdivision, and constitute one of the deepest branching lineage within the gamma-proteobacteria clade.

A "repertoire for repertoire" hypothesis: repertoires of type three effectors are candidate determinants of host specificity in Xanthomonas

Background

The genetic basis of host specificity for animal and plant pathogenic bacteria remains poorly understood. For plant pathogenic bacteria, host range is restricted to one or a few host plant species reflecting a tight adaptation to specific hosts.

Methodology/Principal Findings

Two hypotheses can be formulated to explain host specificity: either it can be explained by the phylogenetic position of the strains, or by the association of virulence genes enabling a pathological convergence of phylogenically distant strains. In this latter hypothesis, host specificity would result from the interaction between repertoires of bacterial virulence genes and repertoires of genes involved in host defences. To challenge these two hypotheses, we selected 132 Xanthomonas axonopodis strains representative of 18 different pathovars which display different host range. First, the phylogenetic position of each strain was determined by sequencing the housekeeping gene rpoD. This study showed that many pathovars of Xanthomonas axonopodis are polyphyletic. Second, we investigated the distribution of 35 type III effector genes (T3Es) in these strains by both PCR and hybridization methods. Indeed, for pathogenic bacteria T3Es were shown to trigger and to subvert host defences. Our study revealed that T3E repertoires comprise core and variable gene suites that likely have distinct roles in pathogenicity and different evolutionary histories. Our results showed a correspondence between composition of T3E repertoires and pathovars of Xanthomonas axonopodis. For polyphyletic pathovars, this suggests that T3E genes might explain a pathological convergence of phylogenetically distant strains. We also identified several DNA rearrangements within T3E genes, some of which correlate with host specificity of strains.

Conclusions/Significance

These data provide insight into the potential role played by T3E genes for pathogenic bacteria and support a “repertoire for repertoire” hypothesis that may explain host specificity. Our work provides resources for functional and evolutionary studies aiming at understanding host specificity of pathogenic bacteria, functional redundancy between T3Es and the driving forces shaping T3E repertoires.

Phylogenetic evidence for extensive horizontal gene transfer of type III secretion system genes among enterobacterial plant pathogens

This study uses sequences from four genes, which are involved in the formation of the type III secretion apparatus, to determine the role of horizontal gene transfer in the evolution of virulence genes for the enterobacterial plant pathogens. Sequences of Erwinia, Brenneria, Pectobacterium, Dickeya and Pantoea were compared (a) with one another, (b) with sequences of enterobacterial animal pathogens, and (c) with sequences of plant pathogenic gamma and beta proteobacteria, to evaluate probable paths of lateral exchange leading to the current distribution of virulence determinants among these micro-organisms. Phylogenies were reconstructed based on hrcC, hrcR, hrcJ and hrcV gene sequences using parsimony and maximum-likelihood algorithms. Virulence gene phylogenies were also compared with several housekeeping gene loci in order to evaluate patterns of lateral versus vertical acquisition. The resulting phylogenies suggest that multiple horizontal gene transfer events have occurred both within and among the enterobacterial plant pathogens and plant pathogenic gamma and beta proteobacteria. hrcJ sequences are the most similar, exhibiting anywhere from 2 to 50 % variation at the nucleotide level, with the highest degree of variation present between plant and animal pathogen sequences. hrcV sequences are conserved among plant and animal pathogens at the N terminus. The C-terminal domain is conserved only among the enterobacterial plant pathogens, as are the hrcC and hrcR sequences. Additionally, hrcJ and hrcV sequence phylogenies suggest that at least some type III secretion system virulence genes from enterobacterial plant pathogens are related more closely to those of the genus Pseudomonas, a conclusion neither supported nor refuted by hrcC or hrcR.

Dealing with incongruence in phylogenomic analyses

Incongruence between gene trees is the main challenge faced by phylogeneticists in the genomic era. Incongruence can occur for artefactual reasons, when we fail to recover the correct gene trees, or for biological reasons, when true gene trees are actually distinct from each other, and from the species tree. Horizontal gene transfers (HGTs) between genomes are an important process of bacterial evolution resulting in a substantial amount of phylogenetic conflicts between gene trees. We argue that the (bacterial) species tree is still a meaningful scientific concept even in the case of HGTs, and that reconstructing it is still a valid goal. We tentatively assess the amount of phylogenetic incongruence caused by HGTs in bacteria by comparing bacterial datasets to a metazoan dataset in which transfers are presumably very scarce or absent.We review existing phylogenomic methods and their ability to return to the user, both the vertical (speciation/extinction history) and horizontal (gene transfers) phylogenetic signals.

Horizontal gene transfer between Ralstonia solanacearum strains detected by comparative genomic hybridization on microarrays

The plant pathogenic Betaproteobacterium Ralstonia solanacearum is a complex species in that most of the strains share the common characteristic of being naturally transformable. In this study, we used a new approach based on comparative genomic hybridization (CGH) on microarrays to investigate the extent of horizontal gene transfers (HGTs) between different strains of R. solanacearum. Recipient strains from phylotypes I, II and III were naturally transformed in vitro by genomic DNA from the GMI1000 reference strain (phylotype I) and the resulting DNAs were hybridized on a microarray representative of the 5120 predicted genes from the GMI1000 strain. In addition to transfer of the antibiotic resistance marker, in 8 of the 16 tested transformants, CGH on microarrays detected other transferred GMI1000 genes and revealed their number, category, function and localization along the genome. We showed that DNA blocks up to 30 kb and 33 genes could be integrated during a single event. Most of these blocks flanked the marker gene DNA but, interestingly, multiple DNA acquisitions along the genome also occurred in a single recombinant clone in one transformation experiment. The results were confirmed by PCR amplification, cloning and sequencing and Southern blot hybridization. This represents the first comprehensive identification of gene acquisitions and losses along the genome of the recipient bacterial strain during natural transformation experiments. In future studies, this strategy should help to answer many questions related to HGT mechanisms.

Acquisition and evolution of plant pathogenesis-associated gene clusters and candidate determinants of tissue-specificity in xanthomonas

Background

Xanthomonas is a large genus of plant-associated and plant-pathogenic bacteria. Collectively, members cause diseases on over 392 plant species. Individually, they exhibit marked host- and tissue-specificity. The determinants of this specificity are unknown.

Methodology/Principal Findings

To assess potential contributions to host- and tissue-specificity, pathogenesis-associated gene clusters were compared across genomes of eight Xanthomonas strains representing vascular or non-vascular pathogens of rice, brassicas, pepper and tomato, and citrus. The gum cluster for extracellular polysaccharide is conserved except for gumN and sequences downstream. The xcs and xps clusters for type II secretion are conserved, except in the rice pathogens, in which xcs is missing. In the otherwise conserved hrp cluster, sequences flanking the core genes for type III secretion vary with respect to insertion sequence element and putative effector gene content. Variation at the rpf (regulation of pathogenicity factors) cluster is more pronounced, though genes with established functional relevance are conserved. A cluster for synthesis of lipopolysaccharide varies highly, suggesting multiple horizontal gene transfers and reassortments, but this variation does not correlate with host- or tissue-specificity. Phylogenetic trees based on amino acid alignments of gum, xps, xcs, hrp, and rpf cluster products generally reflect strain phylogeny. However, amino acid residues at four positions correlate with tissue specificity, revealing hpaA and xpsD as candidate determinants. Examination of genome sequences of xanthomonads Xylella fastidiosa and Stenotrophomonas maltophilia revealed that the hrp, gum, and xcs clusters are recent acquisitions in the Xanthomonas lineage.

Conclusions/Significance

Our results provide insight into the ancestral Xanthomonas genome and indicate that differentiation with respect to host- and tissue-specificity involved not major modifications or wholesale exchange of clusters, but subtle changes in a small number of genes or in non-coding sequences, and/or differences outside the clusters, potentially among regulatory targets or secretory substrates.

The global landscape of sequence diversity

BACKGROUND

Systematic comparisons between genomic sequence datasets have revealed a wide spectrum of sequence specificity from sequences that are highly conserved to those that are specific to individual species. Due to the limited number of fully sequenced eukaryotic genomes, analyses of this spectrum have largely focused on prokaryotes. Combining existing genomic datasets with the partial genomes of 193 eukaryotes derived from collections of expressed sequence tags, we performed a quantitative analysis of the sequence specificity spectrum to provide a global view of the origins and extent of sequence diversity across the three domains of life.

RESULTS

Comparisons with prokaryotic datasets reveal a greater genetic diversity within eukaryotes that may be related to differences in modes of genetic inheritance. Mapping this diversity within a phylogenetic framework revealed that the majority of sequences are either highly conserved or specific to the species or taxon from which they derive. Between these two extremes, several evolutionary landmarks consisting of large numbers of sequences conserved within specific taxonomic groups were identified. For example, 8% of sequences derived from metazoan species are specific and conserved within the metazoan lineage. Many of these sequences likely mediate metazoan specific functions, such as cell-cell communication and differentiation.

CONCLUSION

Through the use of partial genome datasets, this study provides a unique perspective of sequence conservation across the three domains of life. The provision of taxon restricted sequences should prove valuable for future computational and biochemical analyses aimed at understanding evolutionary and functional relationships.

Modular networks and cumulative impact of lateral transfer in prokaryote genome evolution

Lateral gene transfer is an important mechanism of natural variation among prokaryotes, but the significance of its quantitative contribution to genome evolution is debated. Here, we report networks that capture both vertical and lateral components of evolutionary history among 539,723 genes distributed across 181 sequenced prokaryotic genomes. Partitioning of these networks by an eigenspectrum analysis identifies community structure in prokaryotic gene-sharing networks, the modules of which do not correspond to a strictly hierarchical prokaryotic classification. Our results indicate that, on average, at least 81 +/- 15% of the genes in each genome studied were involved in lateral gene transfer at some point in their history, even though they can be vertically inherited after acquisition, uncovering a substantial cumulative effect of lateral gene transfer on longer evolutionary time scales.

Uncovering rate variation of lateral gene transfer during bacterial genome evolution

BACKGROUND

Large scale genome arrangement, such as whole gene insertion/deletion, plays an important role in bacterial genome evolution. Various methods have been employed to study the dynamic process of gene insertions and deletions, such as parsimony methods and maximum likelihood methods. Previous maximum likelihood studies have assumed that the rate of gene insertions/deletions is constant over different genes. This assumption is unrealistic. For instance, it has been shown that informational genes are less likely to be laterally transferred than non-informational genes. However, how much of the variation in gene transfer rates is due to the difference between informational genes and non-informational genes is unclear. In this study, a Gamma-distribution was incorporated in the likelihood estimation by considering rate variation for gene insertions/deletions between genes. This makes it possible to address whether a difference between informational genes and non-informational genes is the main contributor to rate variation of lateral gene transfers.

RESULTS

The results show that models incorporating rate variation fit the data better than do constant rate models in many phylogenetic groups. Even though informational genes are less likely to be laterally transferred than non-informational genes, the degree of rate variation for insertions/deletions did not change dramatically and remained high even when informational genes were excluded from the study. This suggests that the variation in rate of insertions/deletions is not due mainly to the simple difference between informational genes and non-informational genes. Among genes that are not classified as informational and among the informational genes themselves, there are still large differences in the rates that these genes are inserted and deleted.

CONCLUSION

While the difference in informational gene rates contributes to rate variation, it is only a small fraction of the variation present; instead, a substantial amount of rate variation for insertions/deletions remains among both informational genes and among non-informational genes.

Estimating the extent of horizontal gene transfer in metagenomic sequences

BACKGROUND

Although the extent of horizontal gene transfer (HGT) in complete genomes has been widely studied, its influence in the evolution of natural communities of prokaryotes remains unknown. The availability of metagenomic sequences allows us to address the study of global patterns of prokaryotic evolution in samples from natural communities. However, the methods that have been commonly used for the study of HGT are not suitable for metagenomic samples. Therefore it is important to develop new methods or to adapt existing ones to be used with metagenomic sequences.

RESULTS

We have created two different methods that are suitable for the study of HGT in metagenomic samples. The methods are based on phylogenetic and DNA compositional approaches, and have allowed us to assess the extent of possible HGT events in metagenomes for the first time. The methods are shown to be compatible and quite precise, although they probably underestimate the number of possible events. Our results show that the phylogenetic method detects HGT in between 0.8% and 1.5% of the sequences, while DNA compositional methods identify putative HGT in between 2% and 8% of the sequences. These ranges are very similar to these found in complete genomes by related approaches. Both methods act with a different sensitivity since they probably target HGT events of different ages: the compositional method mostly identifies recent transfers, while the phylogenetic is more suitable for the detections of older events. Nevertheless, the study of the number of HGT events in metagenomic sequences from different communities shows a consistent trend for both methods: the lower amount is found for the sequences of the Sargasso Sea metagenome, while the higher quantity is found in the whale fall metagenome from the bottom of the ocean. The significance of these observations is discussed.

CONCLUSION

The computational approaches that are used to find possible HGT events in complete genomes can be adapted to work with metagenomic samples, where a level of high performance is shown in different metagenomic samples. The percentage of possible HGT events that were observed is close to that found for complete genomes, and different microbiomes show diverse ratios of putative HGT events. This is probably related with both environmental factors and the composition in the species of each particular community.

Assessment of the genetic diversity of Xanthomonas axonopodis pv. phaseoli and Xanthomonas fuscans subsp. fuscans as a basis to identify putative pathogenicity genes and a type III secretion system of the SPI-1 family by multiple suppression subtractive hybridizations

Fluorescent amplified fragment length polymorphism revealed that strains of Xanthomonas axonopodis pv. phaseoli and Xanthomonas fuscans subsp. fuscans are genetically distinct and can be grouped into four genetic lineages. Four suppression subtractive hybridizations were then performed to isolate DNA fragments present in these bean pathogens and absent from closely related xanthomonads. Virulence gene candidates were identified such as homologs of hemagglutinins, TonB-dependent receptors, zinc-dependent metalloproteases, type III effectors, and type IV secretion system components. Unexpectedly, homologs of the type III secretion apparatus components (SPI-1 family), usually reported in animal pathogens and insect symbionts, were also detected.

Aeons of distress: an evolutionary perspective on the bacterial SOS response

The SOS response of bacteria is a global regulatory network targeted at addressing DNA damage. Governed by the products of the lexA and recA genes, it co-ordinates a comprehensive response against DNA lesions and its description in Escherichia coli has stood for years as a textbook paradigm of stress-response systems in bacteria. In this paper we review the current state of research on the SOS response outside E. coli. By retracing research on the identification of multiple diverging LexA-binding motifs across the Bacteria Domain, we show how this work has led to the description of a minimum regulon core, but also of a heterogeneous collection of SOS regulatory networks that challenges many tenets of the E. coli model. We also review recent attempts at reconstructing the evolutionary history of the SOS network that have cast new light on the SOS response. Exploiting the newly gained knowledge on LexA-binding motifs and the tight association of LexA with a recently described mutagenesis cassette, these works put forward likely evolutionary scenarios for the SOS response, and we discuss their relevance on the ultimate nature of this stress-response system and the evolutionary pressures driving its evolution.

Population structure and recombination in environmental isolates of Legionella pneumophila

Legionella pneumophila is a water-borne bacteria responsible for most cases of legionellosis, an emerging disease with an increasing incidence in industrialized countries. Although early analysis based on multilocus enzyme electrophoresis (MLEE) described the population structure of this species as clonal, more recent reports have suggested that recombination also contributes to shaping variation across its genome. We report here the results of analysing the nucleotide sequences of 19 loci in 31 environmental samples of L. pneumophila from a small Spanish region (near Alcoi, province of Alicante) where legionellosis has become almost endemic. We analysed the six loci currently incorporated to the sequence-based typing scheme developed by European Working Group for Legionella Infections (EWGLI) for L. pneumophila and 13 intergenic regions, for which we developed primers anchored in flanking, conserved genes. Our results show that recombination among natural isolates of this species is a common phenomenon, as 20 of the 31 isolates contained at least one locus in which recombination was revealed by at least three different methods. The mapping of the recombination events on the maximum likelihood tree of the concatenate sequence of the 19 loci indicated that at least nine independent recombination events might explain the observed distribution of recombinant loci among isolates. In consequence, we have shown that recombination in L. pneumophila is much more frequent than previously considered and that it does not seem to be restricted to already described pathogenicity islands or other genome constituents which provide it with a high plasticity.

From phylogenetics to phylogenomics: the evolutionary relationships of insect endosymbiotic gamma-Proteobacteria as a test case

The increasing availability of complete genome sequences and the development of new, faster methods for phylogenetic reconstruction allow the exploration of the set of evolutionary trees for each gene in the genome of any species. This has led to the development of new phylogenomic methods. Here, we have compared different phylogenetic and phylogenomic methods in the analysis of the monophyletic origin of insect endosymbionts from the gamma-Proteobacteria, a hotly debated issue with several recent, conflicting reports. We have obtained the phylogenetic tree for each of the 579 identified protein-coding genes in the genome of the primary endosymbiont of carpenter ants, Blochmannia floridanus, after determining their presumed orthologs in 20 additional Proteobacteria genomes. A reference phylogeny reflecting the monophyletic origin of insect endosymbionts was further confirmed with different approaches, which led us to consider it as the presumed species tree. Remarkably, only 43 individual genes produced exactly the same topology as this presumed species tree. Most discrepancies between this tree and those obtained from individual genes or by concatenation of different genes were due to the grouping of Xanthomonadales with beta-Proteobacteria and not to uncertainties over the monophyly of insect endosymbionts. As previously noted, operational genes were more prone to reject the presumed species tree than those included in information-processing categories, but caution should be exerted when selecting genes for phylogenetic inference on the basis of their functional category assignment. We have obtained strong evidence in support of the monophyletic origin of gamma-Proteobacteria insect endosymbionts by a combination of phylogenetic and phylogenomic methods. In our analysis, the use of concatenated genes has shown to be a valuable tool for analyzing primary phylogenetic signals coded in the genomes. Nevertheless, other phylogenomic methods such as supertree approaches were useful in revealing alternative phylogenetic signals and should be included in comprehensive phylogenomic studies.

Characterization of ISXax1, a novel insertion sequence restricted to Xanthomonas axonopodis pv. phaseoli (variants fuscans and non-fuscans) and Xanthomonas axonopodis pv. vesicatoria

ISXax1 is a novel insertion sequence belonging to the IS256 and Mutator families. Dot blot, Southern blot, and PCR analyses revealed that ISXax1 is restricted to Xanthomonas axonopodis pv. phaseoli (variants fuscans and non-fuscans) and X. axonopodis pv. vesicatoria strains. Directed AFLP also showed that a high degree of polymorphism is associated with ISXax1 insertion in these strains.