Ecological genomics in Xanthomonas: the nature of genetic adaptation with homologous recombination and host shifts
- Chao-Li Huang†1,
- Pei-Hua Pu†1,
- Hao-Jen Huang†1,
- Huang-Mo Sung†1,
- Hung-Jiun Liaw†1,
- Yi-Min Chen†2,
- Chien-Ming Chen3,
- Ming-Ban Huang1,
- Naoki Osada4,
- Takashi Gojobori1, 4, 5,
- Tun-Wen Pai3,
- Yu-Tin Chen1Email author,
- Chi-Chuan Hwang6Email author and
- Tzen-Yuh Chiang1Email author
© Huang et al.; licensee BioMed Central. 2015
Received: 3 July 2014
Accepted: 20 February 2015
Published: 15 March 2015
Comparative genomics provides insights into the diversification of bacterial species. Bacterial speciation usually takes place with lasting homologous recombination, which not only acts as a cohering force between diverging lineages but brings advantageous alleles favored by natural selection, and results in ecologically distinct species, e.g., frequent host shift in Xanthomonas pathogenic to various plants.
Using whole-genome sequences, we examined the genetic divergence in Xanthomonas campestris that infected Brassicaceae, and X. citri, pathogenic to a wider host range. Genetic differentiation between two incipient races of X. citri pv. mangiferaeindicae was attributable to a DNA fragment introduced by phages. In contrast to most portions of the genome that had nearly equivalent levels of genetic divergence between subspecies as a result of the accumulation of point mutations, 10% of the core genome involving with homologous recombination contributed to the diversification in Xanthomonas, as revealed by the correlation between homologous recombination and genomic divergence. Interestingly, 179 genes were under positive selection; 98 (54.7%) of these genes were involved in homologous recombination, indicating that foreign genetic fragments may have caused the adaptive diversification, especially in lineages with nutritional transitions. Homologous recombination may have provided genetic materials for the natural selection, and host shifts likely triggered ecological adaptation in Xanthomonas. To a certain extent, we observed positive selection nevertheless contributed to ecological divergence beyond host shifting.
Altogether, mediated with lasting gene flow, species formation in Xanthomonas was likely governed by natural selection that played a key role in helping the deviating populations to explore novel niches (hosts) or respond to environmental cues, subsequently triggering species diversification.
KeywordsXanthomonas Adaptive diversification Host shift Parapatric speciation Homologous recombination Comparative genomics
Diversification of prokaryotes attracts much attention from ecologists because it plays a critical role in ecosystem equilibrium and dynamics. Speciation of bacteria is distinct from eukaryotes, especially given the predominantly asexual reproduction  and the rare occurrence of geographic barriers. Bacterial speciation is often triggered by adaptive divergence [2,3], while homologous recombination, which leads to gene flow cohering diverging populations, simultaneously occurs, a model approximating parapatric or sympatric speciation [4,5]. Accordingly, diversifying selection plays a key role in differentiating sister species. On the other hand, favorable genes could be brought by homologous recombination enable the recipient to explore a new niche [6,7], an example of adaptive processes in bacterial diversification under gene flow [8-11].
Many of the bacteria in the genus Xanthomonas, of the γ-subdivision of Proteobacteria, cause plant diseases, e.g., bacterial spots and blights in leaves and fruits . These plant pathogens often display a high degree of host specificity , e.g., X. citri pv. citri exclusively infecting citrus, with various genetic mechanisms associated with the host specificity . Host shifts occurred in pv. mangiferaeindicae that causes bacterial black spot in mango, and in pv. vesicatoria that attacks pepper and tomato, displaying a wide phylogenetic range of hosts in X. citri [15,16]. In contrast, X. campestris pathovars predominantly infect crucifers . It has been known that favorable traits enable deviating populations to explore novel niches in an ecosystem . During habitat shifts such as host transformation, evolutionary footprints of adaptation were often reserved in the genomes . Lu et al. found that six pathogenicity-related gene clusters were associated with the genomic divergences in Xanthomonas . In this study, via a comparative genomic analysis, we comprehensively investigated genomic divergence and adaptation in Xanthomonas, and the contributions of host shifts and homologous recombination to the adaptive diversification between species.
To date, whole genomes of 8 Xanthomonas taxa, including X. campestris pv. campestris (3 strains) [21-23], X. campestris pv. raphani , X. citri pv. vesicatoria , X. citri pv. citrumelo , X. citri pv. citri, X. albilineans , X. oryzae pv. oryzae (3 strains) [27-29], and X. oryzicola , have been sequenced with conventional shotgun sequencing. The size of their genomes varies from 4.9 to 5.4 Mb, with a high GC content of 63.3–69.7% in the chromosome . Phylogenomic analysis revealed a close affinity among X. citri pv. mangiferaeindicae, X. citri pv. citri, and X. citri pv. citrumelo, three taxa that are closely related to X. citri. pv. vesicatoria [25,31]. Recently, a 5.1-Mb genome of X. citri pv. mangiferaeindicae strain LMG 941 from India with 195 contigs was obtained by pyrosequencing . In order to comprehend the genetic divergence among recently diverging strains, we sequenced the genome of BCRC 13182, a local mangiferaeindicae strain from Taiwan. As the two strains LMG 941 and BCRC 13182 of X. citri pv. mangiferaeindicae diverged only recently, how they became genetically differentiated is quite fascinating. In this study, for exploring the nature of adaptive diversification in Xanthomonas, comparative genomic analyses were conducted on two nonsister species complexes of Xanthomonas, i.e., X. citri and X. campestris. The patterns of genomic divergence, homologous recombination, and genes under positive selection were examined to elucidate the ecological interactions among Xanthomonas taxa.
General features of the Xanthomonas citri pv. mangiferaeindicae genome
Genome features of X. citri pv. mangiferaeindicae strains BCRC 13182 (XCM-B) and LMG 941 (XCM-L)
General feature of the genomes
N50 size (bp)
Average length (bp)
Genome size (bp)
G + C content (%)
With assigned to COGs
Shared protein-coding genes
Identical amino acid sequences
Sequences with indels
Phylogenetic relationships and homologous recombination among Xanthomonas species
Recombination can cohere the taxa within a species, but dissociates sister taxa when it occurs between long-split species. Here, we combined the alignment-based programs geneconv and PhiPack to detect genes affected by recombination between or within species. In total, 283 core genes (9.9% of the genome) were detected with genetic recombination, which included 83 genes with cross-overs occurring between X. citri and X. campestris, 172 genes only detected within species, while 28 have likely derived from other distant species (Additional file 3: Table S3). Single gene trees were generated with the neighbor-joining method. For 2786 of the 2851 trees (97.7%), the topologies agreed with the division between citri and campestris, suggesting a long split between the 2 species.
Genome-wide variations and recombination events
Genes under positive selection
Distributions of the genes with recombination and the genes under positive selection in Xanthomonas
Fisher’s exact test
P < 0.0001
P < 0.0001
P < 0.0001
Distributions of positively selected genes and tree topologies in Xanthomonas
Separation of species groups
Fisher’s exact test
P = 1
P = 0.3503
P = 1
Ecological effects of positively selected genes
Positively selected genes
Related to host shift
Related to homologous recombination
Putative ecological effects
RpfB, RpfC (XCM0813, XCM0815)
X. citri pv. citri
Altering pathogenicity to favor exploration of different host
Increasing copper tolerance in escaping from copper-containing bactericides
X. citri pv. vesicatoria
Reducing AMF colonization of host plant and increasing infectivity
I1 branch (X. citri pv. mangiferaeindicae)
Increasing iron uptake with host-shifting from citrus to mango
I2, I3 branches
Controlling invasive routes and host specificity of X. citri
Helping pili assembly and associating with virulence
Breaking down xylan and causing virulence
Speciation processes in Xanthomonas
Bacterial diversity, which often results from species diversification through ecological interactions, greatly influences the ecosystem health. Xanthomonas species, which can cause serious diseases and economic loss of crops, are an excellent group for examining speciation models. The genomes of 8 taxa that belong to X. citri and X. campestris have been completely sequenced, providing a large data set as a sound foundation for investigating the extent of genetic divergence between species. Their host specificity makes the group even more useful as a system for examining adaptation over the host shifts because the host range of the group as a whole is quite broad. In this study, two strains of X. citri pv. mangiferaeindicae from Taiwan (XCM-B) and India (XCM-L) and two strains of X. campestris (XCCA and XCC8) displayed slight genetic differentiation, providing a window for investigating incipient genetic divergence. As expected, given a shorter time for isolation, strains had lower levels of genetic divergence than that of species/pathovars. It is noticeable that prophage-introduced DNA fragments were detected in X. citri pv. mangiferaeindicae and pv. citri , an event that exclusively occurred in X. citri. Several salient features were observed in these inserted DNA fragments. First, long DNA fragments of 28,787 bp were lost in genomes of the two pathovars (Additional file 5: Figure S2). Second, a high proportion of genes (18 of 71 between strains of mangiferaeindicae, and 33 of 71 between pathovars) became eliminated along with the divergence time. Third, higher nucleotide substitution rates were detected in these inserted DNAs than in the host bacterial genome. These facts indicate that the foreign genes introduced by the phage tended to lose functions and were likely to be removed from the genome eventually (e.g., gene XCM0436), although gene residues may still remain .
On the other hand, homologous recombination had started to contribute noticeably to the global diverging patterns among Xanthomonas strains at the incipient stage, as shown by the positive and significant correlation with genomic divergence (Figure 3). Homologous recombination is well documented in prokaryotes with foreign genes introgressed into the recipient species . Recombination often occurs between closely related bacteria, resulting in convergence of two bacterial populations; recombination with divergent bacteria brings variations into the recipient population, facilitating the differentiation among incipient species [42,43].
Furthermore, 1 and 6 genes under positive selection were detected in the lineages leading to XCM-L and XCM-B, respectively. Also, of the 7 genes involved in diversifying selection of pv. mangiferaeindicae, 3 genes were involved in homologous recombination. Similarly, 6 of 9 genes under positive selection were detected with recombination for XCCA and XCC8. At early stages of speciation or diversification, the interplay between natural selection and homologous recombination clearly played a key role in differentiating incipient species/races.
Homologous recombination in Xanthomonas
In this study, we found frequent genetic recombination in Xanthomonas, with about 10% of the core genome mediated with homologous recombination, whereas only 2% of the core genes displayed phylogenies deviating from the species tree. This gap to the expectations simply came from intraspecific crossovers, which constitute the most recombination events (about 88%) with short DNA fragments. Apparently, those genes displaying trees inconsistent with the species tree represented footprints of gene flow between species. This result suggests that the divergence of Xanthomonas tends to follow parapatric speciation, in which gene flow between species can be long lasting and nonstop . This is a pattern well documented in the Archaea and Bacteria [1,7]. It is intriguing that most of the foreign genes were favored by natural selection (Table 2), especially those associated with plant infection (Figure 4B), reflecting a fact that homologous recombination may create advantageous new alleles for a so-called “hopeful monster”  in bacteria. Altogether, the topology of the putative species tree in Xanthomonas reflects a scenario of deep divergence between species, mediated by recurrent gene flow at the same time .
In this study, 98 out of 283 genes with recombination were under positive selection (Table 2). Via homologous recombination, genes were able to be interchanged between strains, and advantageous alleles that helped bacterial colonization would be favored by natural selection [3,45]. Of these genes, it is noticeable that RpfB (XCM0815) and RpfC (XCM0813) of the Rpf gene cluster, which regulates the virulence of Xanthomonas [20,46,47], were positively selected in X. citri pv. citri lineage (Table 4). Functionally, RpfB encodes a fatty acyl-CoA ligase, which catalyzes the synthesis of an important signal molecule that regulates the expression of virulence genes [48,49]. RpfB signaling is perceived by the two component system of RpfC and RpfG , and the mutations of RpfC were found to reduce the pathogenicity of X. citri pv. citri . Accordingly, existence of favorable RpfB alleles suggests that homologous recombination may have helped the adaptation of Xanthomonas. Besides, DsbD (XCM1788) has been reported to be associated with copper tolerance . Copper compounds are frequently used as bactericides in controlling the leaf infection of Xanthomonas species, including X. citri [52,53]. A previous experiment showed that a DsbD knockout mutant was highly sensitive to environmental copper . An effective DsbD allele was capable of reducing the impacts of the copper-containing bactericides, thus enabling the bacterium to escape from the agricultural control.
Host shifts likely triggering genetic adaptation
Like the findings of previous studies, compared with the vast portions of the genome that are usually shaped by purifying selection, adaptive genes only constitute a relatively smaller proportion of the chromosome instead . Furthermore, in this study, genes under positive selection (only 6.3% of the core genome) were unequally distributed along lineages, with significantly more genes located at lineages leading to nodes coupled with host shifts (Figure 4A). For example, at the lineages from X. citri pv. citri leading to X. citri pv. mangiferaeindicae and X. citri pv. vescatoria with hosts shifts among citrus, mango and pepper, respectively, 18 and 25 positively selected genes were detected, revealing faster accumulation of positively selected genes than other lineages without host shifts. However, the lineages of X. campestris pv. raphani and pv. campestris, both infecting the Brassicaceae without host shifting, did not show facilitated positive selection. This contrasting pattern suggested that favorable genes may have helped the bacteria in exploring new niches. The tight association between positively selected genes and host shifts further suggests that Xanthomonas likely followed ecological speciation, which describes species arising from ecological diversification .
We found that some positively selected genes were associated with the host shifting events. It is notable that the PotG gene (XCM0993), that encodes ATPase involved in putrescine transport, was shaped by positive selection in the X. citri pv. vesicatoria lineage. Putrescine is a ubiquitous polyamine that enhances the interaction between plants and arbuscular mycorrhizal fungi (AMF) , which in turn have been shown to activate plant defense against infectious pathogens such as Xanthomonas . Alternatively, the removal of exogenous putrescine from the soil may impede host colonization by AMF; thus, an effective PotG gene that decreases AMF colonization of host plant will be favored by natural selection, subsequently helping X. citri pv. vesicatoria to explore a new niche (host). In addition, BtuB gene (XCM3148) was positively selected in X. citri pv. mangiferaeindicae lineage. BtuB is responsible for the iron uptake, while iron dependent superoxide dismutases are vital in inhibiting the reactive oxygen species (ROS) responses in host cells, thus increasing the infection rates . Interestingly, previous studies revealed that the iron contents in mango leaves might be hundreds times lower than that in citrus leaves (0.27 vs. 41.7 mg kg-1 dry weight) [58,59]. The sharp difference suggested that an effective iron transporter was needed when the host shifted from citrus to mango.
In addition to the genes related to the host shifting, an Hrp regulon gene XCM1612 that encodes components for the type III secretion system (T3SS)  was detected as loci under positive selection in the X. citri. The pathogenicity of Xanthomonas mainly depends on the T3SS, which is highly conserved among plant and animal pathogenic bacteria [61,62]. Curiously, the amino acid sequences of the major subunit of Hrp pili are hypervariable in different subspecies of bacterial pathogens [63,64]. The rapidly evolving Hrp pili provide evidence for strong positive selection in Xanthomonas spp. and diversifying selection in Pseudomonas syringae [65,66]. The acquisition of the Hrp gene cluster has been found to be associated with the adaptation and plant pathogenicity in Xanthomonas . Thus, we hypothesized that the positive selection on HrpF gene may be responsible for the invasion routes and the host specificity of Xanthomonas. Furthermore, two genes involved in type IV secretion system were detected with positive selections exclusively in X. campestris, seemingly agreeing that type IV secretion system may not be involved in the infection of X. citri on citrus . Moreover, EcpD (XCM0369) that encodes an adhesion protein to help pili assembly was positively selected in XCM-L. EcpD has been shown to facilitate the polymerization of EcpA to form pilus in Escherichia coli and involve in host cell recognition or biofilm formation . Xanthomonas species only possessed EcpD while lacks other members of Ecp operon. We therefore hypothesized that EcpD might participate the assembly of the other pili in association with the virulence in Xanthomonas. On the other hand, a xylanase gene XynA (XCM1323) was positively selected at the branch leading to X. campestris pv. campestris (Figure 4B, Table 4 ). Xylan is a major component in the cell walls of land plants and exists in all plant tissues . Previous studies showed xylanases are responsible for the virulence of X. citri pv. vesicatoria and X. oryzae pv. oryzae [70,71]. Two gene clusters xcs and xps of the Type II secretion system (T2SS) have been shown to control the secretion of xylanases in XCV and are associated with the virulence . Nevertheless, no member of these two gene clusters was detected under positive selection in this study (Additional file 3: Table S3).
In this study, we sequenced a genome of X. citri pv. mangiferaeindicae and conducted genomic analyses of 9 taxa of X. citri and campestris. Between the 2 strains of XCM, only the prophage region displayed sharp differentiation, while gradually losing the constituting genes. In addition, we found homologous recombination frequently occurring in the Xanthomonas genomes, which likely represented footprints of gene flow between species, thus most likely suggesting parapatric speciation. It is noticeable that facilitated accumulations of positively selected genes occurred along the lineages with host shifts. Interestingly, most of the favored genes were acquired from homologous recombination. Taken together, the genes with recombination enabled Xanthomonas species to explore novel niches and respond to environmental stresses, subsequently resulting in adaptive diversification in this pathogenic genus.
Whole sequencing, assembly, and annotation of Xanthomonas citri pv. mangiferaeindicae
Bacterial strain BCRC 13182 (NCHUPP Xma1) from Taiwan of X. citri pv. mangiferaeindicae was purchased from the Food Industry Research and Development Institute and sequenced in this study. The whole genome was sequenced using a high-throughput sequencing technique with Illumina GA IIx and HiSeq sequencers . Paired-end sequence reads of 60-bp were obtained, with an average distance of 150 bp between pair reads. Mate-paired sequence reads of 100 bp were obtained with an average distance of 3000 bp between pair reads. All raw sequences were stored in NCBI SRA database (SRP049288). A total of 110,671,234 and 15,543,685 high-quality sequences of 60 and 100 bp, respectively, were assembled into contigs using the de novo assembler ABySS version 1.2.7 . After contigs shorter than 200 bp were removed, 221 contigs were generated. Scaffolding was performed with these contigs using SSPACE v2.0 , and 43 scaffolds were built based on the distance information between paired reads.
The draft genome was analyzed using an integrated annotation pipeline glued by the perl programming language. Glimmer version 3.02  was used to predict the protein-coding regions, which were annotated with BLAST  against the nonredundant protein database (http://ncbi.nlm.nih.gov) and the Clusters of Orthologous Groups (COG) database  with a cutoff of E-value < 1 × 10-5. All protein-coding regions were manually curated with the EMBOSS analysis package according to the BLAST results . We used tRNAscan-SE and RNAmmer version 1.2 to predict the prokaryotic transfer RNA genes and the ribosomal RNA genes [78,79].
Comparative genomics analysis and orthologous protein identification
For comparative analyses, we downloaded 8 bacteria genomes from the NCBI GenBank (http://www.ncbi.nlm.nih.gov/genbank/) in addition to the newly assembled genomes of X. citri pv. mangiferaeindicae BCRC 13182, including X. citri pv. citri 306, X. citri pv. citrumelo F1, X. citri pv. vesicatoria 85–10, X. campestris pv. campestris ATCC 33913, X. campestris pv. campestris 8004, X. campestris pv. campestris B100, X. campestris pv. raphani 756C, and a draft genome of X. citri pv. mangiferaeindicae LMG 941 (Additional file 1: Table S1). Orthologous genes in these genomes were identified using bidirectional best hits (BBH) in a BLAST search  based on the criteria set at E-value < 1 × 10-5, identity > 60%, and a threshold of 70% of orthologous alignment length. The core genome with 2,851 orthologous genes and common proteins shared by Xanthomonas species were recognized. In total, 43 scaffolds of the X. citri pv. mangiferaeindicae strain BCRC 13182 were aligned against the strain LMG 941, which has 195 contigs. We used BLAST for pairwise nucleotide alignment of the two draft genomes.
Phylogeny of Xanthomonas
For aligning nucleotide sequences of protein-encoding genes without interrupting codons, the protein sequences were first aligned by ClustalW2 , and the corresponding nucleotide sequences were aligned accordingly using PAL2NAL . A maximum-likelihood (ML) phylogeny based on the concatenated 2851 genes was generated with the program RAxML . We used the GTRGAMMA substitution model suggested by jModelTest 2.1.4 [84,85]. Confidence in the internal nodes of the ML tree was tested with 1000 rapid bootstrap replicates. For each gene phylogeny, we used ClustalW2 to generate a neighbor-joining tree with default settings.
Detection of genetic recombination
Putative genetic recombination events occurring in Xanthomonas genomes were examined by using the software geneconv  and PhiPack . First, the alignments of the 2851 orthologous genes were tested with geneconv, and a recombination event can be recognized when a Bonferroni corrected KA p-value is less than 0.05. Notably, the geneconv program is able to distinguish a recombination occurring within or between species groups. In addition, to reduce the probability of detecting false positive when merely using a single method, we used the PhiPack program, which performs three different methods (Neighbor similarity score, MaxChi, and Phi) for identifying genetic recombination, to confirm the results of geneconv. Using the PhiPack program, recombination was identified for p-values lower than 0.05. Taken together, a gene with recombination was recognized when recombination was detected with geneconv and at least two additional methods implemented in PhiPack.
Numbers of substitutions per synonymous site (K s) and per nonsynonymous site (K a) were calculated using codeml program implemented in the PAML package with the codon-based model . All core orthologous genes were mapped onto the XCC genome to perform window-sliding analyses. The genome-wide comparisons of K s values and the number of genes with recombination were made with a window size of 50 genes and an overlap of 10 genes. Spearman’s rho coefficient was used to examine the correlations between Ks values and the number of the genes with recombination. In addition, the significant of difference of Ks distributions between genes with and without recombination were evaluated with Mann–Whitney U test. Both tests were performed by using SPSS Statistics 17.0.
Detection of genes under positive selection
For detecting genes under positive selection, we applied the codeml program from the PAML package to all core genes. The multiple alignments were fit to the F3X4 model of codon frequencies, and the branch-site alternative model (model = 2, NSsites = 2) was adopted . Three independent simulations were performed with initial omega (K a/K s) values of 0.5, 1, and 1.5, respectively, while the null model was fixed with the omega value at 1. Likelihood ratio tests were used to assess the significance between the test and null models, and the P values were adjusted with a false discovery rate (FDR) of 0.01 for multiple testing corrections. We performed codeml analyses on each branch of individual gene trees. For above three tests, branches detected with significance were deemed to be affected by positive selection. Assuming the number of positively selected genes was proportional to the branch length, we used a chi-square test to examine the homogeneity in the distribution of positively selected genes between branches with and without host shifts. Fisher’s exact tests were also used to assess the differences of the number of positively selected genes between outlier and background genes, as well as between the genes with and without homologous recombination.
Availability of supporting data
Alignments and phylogenetic trees of the core orthologous genes were deposited in TreeBASE ( http://treebase.org/treebase-web/home.html ).
- Cadillo-Quiroz H, Didelot X, Held NL, Herrera A, Darling A, Reno ML, et al. Patterns of gene flow define species of thermophilic Archaea. PLoS Biol. 2012;10:e1001265.View ArticlePubMed CentralPubMedGoogle Scholar
- Kopac S, Wang Z, Wiedenbeck J, Sherry J, Wu M, Cohan FM. Genomic heterogeneity and ecological speciation within one subspecies of Bacillus subtilis. Appl Environ Microbiol. 2014;80:4842–53.View ArticlePubMed CentralPubMedGoogle Scholar
- Vos M. A species concept for bacteria based on adaptive divergence. Trends Microbiol. 2011;19:1–7.View ArticlePubMedGoogle Scholar
- Oakley BB, Carbonero F, van der Gast CJ, Hawkins RJ, Purdy KJ. Evolutionary divergence and biogeography of sympatric niche-differentiated bacterial populations. ISME J. 2010;4:488–97.View ArticlePubMedGoogle Scholar
- Sikorski J, Nevo E. Adaptation and incipient sympatric speciation of Bacillus simplex under microclimatic contrast at “Evolution Canyons” I and II, Israel. Proc Natl Acad Sci U S A. 2005;102:15924–9.View ArticlePubMed CentralPubMedGoogle Scholar
- Cohan FM, Koeppel AF. The origins of ecological diversity in prokaryotes. Curr Biol. 2008;18:R1024–34.View ArticlePubMedGoogle Scholar
- Vos M, Didelot X. A comparison of homologous recombination rates in bacteria and archaea. ISME J. 2008;3:199–208.View ArticlePubMedGoogle Scholar
- Luo C, Walk ST, Gordon DM, Feldgarden M, Tiedje JM, Konstantinidis KT. Genome sequencing of environmental Escherichia coli expands understanding of the ecology and speciation of the model bacterial species. Proc Natl Acad Sci U S A. 2011;108:7200–5.View ArticlePubMed CentralPubMedGoogle Scholar
- Cohan FM. Genetic exchange and evolutionary divergence in prokaryotes. Trends Ecol Evol. 1994;9:175–80.View ArticlePubMedGoogle Scholar
- Nogueira T, Rankin DJ, Touchon M, Taddei F, Brown SP, Rocha EP. Horizontal gene transfer of the secretome drives the evolution of bacterial cooperation and virulence. Curr Biol. 2009;19:1683–91.View ArticlePubMed CentralPubMedGoogle Scholar
- Wiedenbeck J, Cohan FM. Origins of bacterial diversity through horizontal genetic transfer and adaptation to new ecological niches. FEMS Microbiol Rev. 2011;35:957–76.View ArticlePubMedGoogle Scholar
- Vauterin L, Hoste B, Kersters K, Swings J. Reclassification of xanthomonas. Int J Syst Bacteriol. 1995;45:472–89.View ArticleGoogle Scholar
- Dye D, Bradbury J, Goto M, Hayward A, Lelliott R, Schroth MN. International standards for naming pathovars of phytopathogenic bacteria and a list of pathovar names and pathotype strains. Rev Plant Pathol. 1980;59:153–68.Google Scholar
- Siciliano F, Torres P, Sendín L, Bermejo C, Filippone P, Vellice G, et al. Analysis of the molecular basis of Xanthomonas axonopodis pv. citri pathogenesis in Citrus limon. Electron J Biotechnol. 2006;9:200–4.View ArticleGoogle Scholar
- Gagnevin L, Pruvost O. Epidemiology and control of mango bacterial black spot. Plant Dis. 2001;85:928–35.View ArticleGoogle Scholar
- Mhedbi-Hajri N, Hajri A, Boureau T, Darrasse A, Durand K, Brin C, et al. Evolutionary history of the plant pathogenic bacterium Xanthomonas axonopodis. PLoS One. 2013;8:e58474.View ArticlePubMed CentralPubMedGoogle Scholar
- Bogdanove AJ, Koebnik R, Lu H, Furutani A, Angiuoli SV, Patil PB, et al. Two new complete genome sequences offer insight into host and tissue specificity of plant pathogenic Xanthomonas spp. J Bacteriol. 2011;193:5450–64.View ArticlePubMed CentralPubMedGoogle Scholar
- Huang C-L, Ho C-W, Chiang Y-C, Shigemoto Y, Hsu T-W, Hwang C-C, et al. Adaptive divergence with gene flow in incipient speciation of Miscanthus floridulus/sinensis complex (Poaceae). Plant J. 2014;80:834–47.View ArticlePubMedGoogle Scholar
- Schmidt H, Greshake B, Feldmeyer B, Hankeln T, Pfenninger M. Genomic basis of ecological niche divergence among cryptic sister species of non-biting midges. BMC Genomics. 2013;14:384.View ArticlePubMed CentralPubMedGoogle Scholar
- Lu H, Patil P, Van Sluys M-A, White FF, Ryan RP, Dow JM, et al. Acquisition and evolution of plant pathogenesis–associated gene clusters and candidate determinants of tissue-specificity in Xanthomonas. PLoS One. 2008;3:e3828.View ArticlePubMed CentralPubMedGoogle Scholar
- da Silva AR, Ferro JA, Reinach F, Farah C, Furlan L, Quaggio R, et al. Comparison of the genomes of two Xanthomonas pathogens with differing host specificities. Nature. 2002;417:459–63.View ArticlePubMedGoogle Scholar
- Qian W, Jia Y, Ren S-X, He Y-Q, Feng J-X, Lu L-F, et al. Comparative and functional genomic analyses of the pathogenicity of phytopathogen Xanthomonas campestris pv. campestris. Genome Res. 2005;15:757–67.View ArticlePubMed CentralPubMedGoogle Scholar
- Vorhölter F-J, Schneiker S, Goesmann A, Krause L, Bekel T, Kaiser O, et al. The genome of Xanthomonas campestris pv. campestris B100 and its use for the reconstruction of metabolic pathways involved in xanthan biosynthesis. J Biotechnol. 2008;134:33–45.View ArticlePubMedGoogle Scholar
- Thieme F, Koebnik R, Bekel T, Berger C, Boch J, Büttner D, et al. Insights into genome plasticity and pathogenicity of the plant pathogenic bacterium Xanthomonas campestris pv. vesicatoria revealed by the complete genome sequence. J Bacteriol. 2005;187:7254–66.View ArticlePubMed CentralPubMedGoogle Scholar
- Jalan N, Aritua V, Kumar D, Yu F, Jones JB, Graham JH, et al. Comparative genomic analysis of xanthomonas axonopodis pv. citrumelo f1, which causes citrus bacterial spot disease, and related strains provides insights into virulence and host specificity. J Bacteriol. 2011;193:6342–57.View ArticlePubMed CentralPubMedGoogle Scholar
- Pieretti I, Royer M, Barbe V, Carrere S, Koebnik R, Cociancich S, et al. The complete genome sequence of Xanthomonas albilineans provides new insights into the reductive genome evolution of the xylem-limited Xanthomonadaceae. BMC Genomics. 2009;10:616.View ArticlePubMed CentralPubMedGoogle Scholar
- Ochiai H, Inoue Y, Takeya M, Sasaki A, Kaku H. Genome sequence of Xanthomonas oryzae pv. oryzae suggests contribution of large numbers of effector genes and insertion sequences to its race diversity. JARQ-Jpn Agr Res Q. 2005;39:275.View ArticleGoogle Scholar
- Lee B-M, Park Y-J, Park D-S, Kang H-W, Kim J-G, Song E-S, et al. The genome sequence of Xanthomonas oryzae pathovar oryzae KACC10331, the bacterial blight pathogen of rice. Nucleic Acids Res. 2005;33:577–86.View ArticlePubMed CentralPubMedGoogle Scholar
- Salzberg SL, Sommer DD, Schatz MC, Phillippy AM, Rabinowicz PD, Tsuge S, et al. Genome sequence and rapid evolution of the rice pathogen Xanthomonas oryzae pv. oryzae PXO99A. BMC Genomics. 2008;9:204.View ArticlePubMed CentralPubMedGoogle Scholar
- Fleischmann RD, Adams MD, White O, Clayton RA, Kirkness EF, Kerlavage AR, et al. Whole-genome random sequencing and assembly of Haemophilus influenzae Rd. Science. 1995;269:496–512.View ArticlePubMedGoogle Scholar
- Midha S, Ranjan M, Sharma V, Pinnaka AK, Patil PB. Genome sequence of Xanthomonas citri pv. mangiferaeindicae strain LMG 941. J Bacteriol. 2012;194:3031.View ArticlePubMed CentralPubMedGoogle Scholar
- Ah-You N, Gagnevin L, Chiroleu F, Jouen E, Neto JR, Pruvost O. Pathological variations within Xanthomonas campestris pv. mangiferaeindicae support its separation into three distinct pathovars that can be distinguished by amplified fragment length polymorphism. Phytopathology. 2007;97:1568–77.View ArticlePubMedGoogle Scholar
- Du P, Yang Y, Wang H, Liu D, Gao G, Chen C. A large scale comparative genomic analysis reveals insertion sites for newly acquired genomic islands in bacterial genomes. BMC Microbiol. 2011;11:135.View ArticlePubMed CentralPubMedGoogle Scholar
- Caza M, Kronstad JW. Shared and distinct mechanisms of iron acquisition by bacterial and fungal pathogens of humans. Front Cell Infect Microbiol. 2013;3:80.View ArticlePubMed CentralPubMedGoogle Scholar
- Gotoh Y, Eguchi Y, Watanabe T, Okamoto S, Doi A, Utsumi R. Two-component signal transduction as potential drug targets in pathogenic bacteria. Curr Opin Microbiol. 2010;13:232–9.View ArticlePubMedGoogle Scholar
- Alfano JR, Collmer A. Type III secretion system effector proteins: double agents in bacterial disease and plant defense. Annu Rev Phytopathol. 2004;42:385–414.View ArticlePubMedGoogle Scholar
- Christie PJ, Vogel JP. Bacterial type IV secretion: conjugation systems adapted to deliver effector molecules to host cells. Trends Microbiol. 2000;8:354–60.View ArticlePubMedGoogle Scholar
- Brown JS, Holden DW. Iron acquisition by gram-positive bacterial pathogens. Microb Infect. 2002;4:1149–56.View ArticleGoogle Scholar
- Marraffini LA, Sontheimer EJ. CRISPR interference limits horizontal gene transfer in staphylococci by targeting DNA. Science. 2008;322:1843–5.View ArticlePubMed CentralPubMedGoogle Scholar
- Chiang T-Y, Chiang YC, Pan C-H, Wang W-K, Chen C-N, Hsu T-W, et al. Comparative genomics of horizontal transfer between chloroplast and nuclear genomes in rice and Arabidopsis. In: Plant evolutionary genetics and biology of weeds. Nantou, Taiwan: Endemic Species Research Institute; 2004. p. 1–10.Google Scholar
- Mallet J. Hybridization, ecological races and the nature of species: empirical evidence for the ease of speciation. Phil Trans R Soc B. 2008;363:2971–86.View ArticlePubMed CentralPubMedGoogle Scholar
- Didelot X, Maiden MCJ. Impact of recombination on bacterial evolution. Trends Microbiol. 2010;18:315–22.View ArticlePubMed CentralPubMedGoogle Scholar
- Guttman DS, Dykhuizen DE. Clonal divergence in Escherichia coli as a result of recombination, not mutation. Science. 1994;266:1380–3.View ArticlePubMedGoogle Scholar
- Goldschmidt R. The material basis of evolution. New Haven, CT: Yale University Press; 1940.Google Scholar
- Koeppel A, Perry EB, Sikorski J, Krizanc D, Warner A, Ward DM, et al. Identifying the fundamental units of bacterial diversity: a paradigm shift to incorporate ecology into bacterial systematics. Proc Natl Acad Sci U S A. 2008;105:2504–9.View ArticlePubMed CentralPubMedGoogle Scholar
- Roohani R, Najafabadi MS, Alavi SM, Farrokhi N, Shams-bakhsh M. Transcript analysis of some pathogenicity-related elements in an Iranian A* isolate of Xanthomonas citri subsp. citri. J Crop Prot. 2012;1:337–47.Google Scholar
- Büttner D, Bonas U. Regulation and secretion of Xanthomonas virulence factors. FEMS Microbiol Rev. 2010;34:107–33.View ArticlePubMedGoogle Scholar
- Ryan RP, Fouhy Y, Lucey JF, Jiang B-L, He Y-Q, Feng J-X, et al. Cyclic di-GMP signalling in the virulence and environmental adaptation of Xanthomonas campestris. Mol Microbiol. 2007;63:429–42.View ArticlePubMedGoogle Scholar
- Slater H, Alvarez-Morales A, Barber CE, Daniels MJ, Dow JM. A two-component system involving an HD-GYP domain protein links cell–cell signalling to pathogenicity gene expression in Xanthomonas campestris. Mol Microbiol. 2000;38:986–1003.View ArticlePubMedGoogle Scholar
- Ryan RP, Fouhy Y, Lucey JF, Crossman LC, Spiro S, He YW, et al. Cell-cell signaling in Xanthomonas campestris involves an HD-GYP domain protein that functions in cyclic di-GMP turnover. Proc Natl Acad Sci U S A. 2006;103:6712–7.View ArticlePubMed CentralPubMedGoogle Scholar
- Gupta SD, Wu HC, Rick PD. A Salmonella typhimurium genetic locus which confers copper tolerance on copper-sensitive mutants of Escherichia coli. J Bacteriol. 1997;179:4977–84.PubMed CentralPubMedGoogle Scholar
- Palmieri ACB, Amaral AM, Homem RA, Machado MA. Differential expression of pathogenicity- and virulence-related genes of Xanthomonas axonopodis pv. citri under copper stress. Genet Mol Biol. 2010;33:348–53.View ArticlePubMed CentralPubMedGoogle Scholar
- Ritchie D. Copper-containing fungicides/bactericides and their use in management of bacterial spot on peaches. Southeast Reg Newsl. 2004;4:1.Google Scholar
- Katzen F, Beckwith J. Role and location of the unusual redox-active cysteines in the hydrophobic domain of the transmembrane electron transporter DsbD. Proc Natl Acad Sci U S A. 2003;100:10471–6.View ArticlePubMed CentralPubMedGoogle Scholar
- Sella G, Petrov DA, Przeworski M, Andolfatto P. Pervasive natural selection in the Drosophila genome? PLoS Genet. 2009;5:e1000495.View ArticlePubMed CentralPubMedGoogle Scholar
- Wu Q-S, Zou Y-N, He X-H. Exogenous putrescine, not spermine or spermidine, enhances root mycorrhizal development and plant growth of trifoliate orange (Poncirus trifoliata) seedlings. Int J Agr Biol. 2010;12:576–80.Google Scholar
- Erb M, Lenk C, Degenhardt J, Turlings TC. The underestimated role of roots in defense against leaf attackers. Trends Plant Sci. 2009;14:653–9.View ArticlePubMedGoogle Scholar
- Pestana M, de Varennes A, Goss M, Abadía J, Faria E. Floral analysis as a tool to diagnose iron chlorosis in orange trees. Plant Soil. 2004;259:287–95.View ArticleGoogle Scholar
- Hossain T, Alam Z, Absar N. Changes in different nutrients and enzyme contents in mango leaves infected with Colletotrichum glaeosorioides. Ind Phytopayhol. 1999;52:75–6.Google Scholar
- Cornelis GR. The type III secretion injectisome. Nat Rev Microbiol. 2006;4:811–25.View ArticlePubMedGoogle Scholar
- Cornelis GR, Van Gijsegem F. Assembly and function of type III secretory systems. Annu Rev Microbiol. 2000;54:735–74.View ArticlePubMedGoogle Scholar
- Tampakaki A, Fadouloglou V, Gazi A, Panopoulos N, Kokkinidis M. Conserved features of type III secretion. Cell Microbiol. 2004;6:805–16.View ArticlePubMedGoogle Scholar
- Lee YH, Kolade OO, Nomura K, Arvidson DN, He SY. Use of dominant-negative HrpA mutants to dissect Hrp pilus assembly and type III secretion in Pseudomonas syringae pv. tomato. J Biol Chem. 2005;280:21409–17.View ArticlePubMedGoogle Scholar
- Weber E, Koebnik R. Domain structure of HrpE, the Hrp pilus subunit of Xanthomonas campestris pv. vesicatoria. J Bacteriol. 2005;187:6175–86.View ArticlePubMed CentralPubMedGoogle Scholar
- Weber E, Koebnik R. Positive selection of the Hrp pilin HrpE of the plant pathogen Xanthomonas. J Bacteriol. 2006;188:1405–10.View ArticlePubMed CentralPubMedGoogle Scholar
- Guttman DS, Gropp SJ, Morgan RL, Wang PW. Diversifying selection drives the evolution of the type III secretion system pilus of Pseudomonas syringae. Mol Biol Evol. 2006;23:2342–54.View ArticlePubMedGoogle Scholar
- Jacob TR, de Laia ML, Moreira LM, Goncalves JF, de Souza Carvalho FM, Ferro MIT, et al. Type IV secretion system is not involved in infection process in citrus. Int J Microbiol. 2014;2014:9.View ArticleGoogle Scholar
- Garnett JA, Martínez-Santos VI, Saldaña Z, Pape T, Hawthorne W, Chan J, et al. Structural insights into the biogenesis and biofilm formation by the Escherichia coli common pilus. Proc Natl Acad Sci U S A. 2012;109:3950–5.View ArticlePubMed CentralPubMedGoogle Scholar
- Wong KKY, Tan LUL, Saddler JN. Multiplicity of beta-1,4-xylanase in microorganisms: functions and applications. Microbiol Rev. 1988;52:305–17.PubMed CentralPubMedGoogle Scholar
- Szczesny R, Jordan M, Schramm C, Schulz S, Cogez V, Bonas U, et al. Functional characterization of the Xcs and Xps type II secretion systems from the plant pathogenic bacterium Xanthomonas campestris pv vesicatoria. New Phytol. 2010;187:983–1002.View ArticlePubMedGoogle Scholar
- Rajeshwari R, Jha G, Sonti RV. Role of an in planta-expressed xylanase of Xanthomonas oryzae pv. oryzae in promoting virulence on rice. Mol Plant-Microbe Interact. 2005;18:830–7.View ArticlePubMedGoogle Scholar
- Bentley DR. Whole-genome re-sequencing. Curr Opin Genet Dev. 2006;16:545–52.View ArticlePubMedGoogle Scholar
- Simpson JT, Wong K, Jackman SD, Schein JE, Jones SJM, Birol İ. ABySS: a parallel assembler for short read sequence data. Genome Res. 2009;19:1117–23.View ArticlePubMed CentralPubMedGoogle Scholar
- Boetzer M, Henkel CV, Jansen HJ, Butler D, Pirovano W. Scaffolding pre-assembled contigs using SSPACE. Bioinforma. 2011;27:578–9.View ArticleGoogle Scholar
- Delcher AL, Harmon D, Kasif S, White O, Salzberg SL. Improved microbial gene identification with GLIMMER. Nucleic Acids Res. 1999;27:4636–41.View ArticlePubMed CentralPubMedGoogle Scholar
- Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215:403–10.View ArticlePubMedGoogle Scholar
- Tatusov RL, Fedorova ND, Jackson JD, Jacobs AR, Kiryutin B, Koonin EV, et al. The COG database: an updated version includes eukaryotes. BMC Bioinforma. 2003;4:41.View ArticleGoogle Scholar
- Lowe TM, Eddy SR. tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 1997;25:0955–64.View ArticleGoogle Scholar
- Lagesen K, Hallin P, Rødland EA, Stærfeldt HH, Rognes T, Ussery DW. RNAmmer: consistent and rapid annotation of ribosomal RNA genes. Nucleic Acids Res. 2007;35:3100–8.View ArticlePubMed CentralPubMedGoogle Scholar
- Fuchsman CA, Rocap G. Whole-genome reciprocal BLAST analysis reveals that planctomycetes do not share an unusually large number of genes with Eukarya and Archaea. Appl Environ Microbiol. 2006;72:6841–4.View ArticlePubMed CentralPubMedGoogle Scholar
- Larkin M, Blackshields G, Brown N, Chenna R, McGettigan P, McWilliam H, et al. Clustal W and Clustal X version 2.0. Bioinformatics. 2007;23:2947–8.View ArticlePubMedGoogle Scholar
- Suyama M, Torrents D, Bork P. PAL2NAL: robust conversion of protein sequence alignments into the corresponding codon alignments. Nucleic Acids Res. 2006;34:W609–12.View ArticlePubMed CentralPubMedGoogle Scholar
- Stamatakis A. RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinforma. 2006;22:2688–90.View ArticleGoogle Scholar
- Darriba D, Taboada GL, Doallo R, Posada D. jModelTest 2: more models, new heuristics and parallel computing. Nat Methods. 2012;9:772.View ArticlePubMedGoogle Scholar
- Guindon S, Gascuel O. A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2003;52:696–704.View ArticlePubMedGoogle Scholar
- Padidam M, Sawyer S, Fauquet CM. Possible emergence of new geminiviruses by frequent recombination. Virol. 1999;265:218–25.View ArticleGoogle Scholar
- Bruen TC, Philippe H, Bryant D. A simple and robust statistical test for detecting the presence of recombination. Genet. 2006;172:2665–81.View ArticleGoogle Scholar
- Yang Z. PAML 4: phylogenetic analysis by maximum likelihood. Mol Biol Evol. 2007;24:1586–91.View ArticlePubMedGoogle Scholar
- Yang Z, Wong WS, Nielsen R. Bayes empirical Bayes inference of amino acid sites under positive selection. Mol Biol Evol. 2005;22:1107–18.View ArticlePubMedGoogle Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.