This item is only available as the following downloads:
Moreira et al. BMC Genomics 2010, 11 :238 http://www.biomedcentra l.com/1471-2 164/11/238 Open Access RESEARCH ARTICLE 2010 Moreira et al; licensee BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative C ommons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestri cted use, distributi on, and reproductio n in any medium, provided the original work is properly cited.Research articleNovel insights into the genomic basis of citrus canker based on the genome sequences of two strains of Xanthomonas fuscans subsp. aurantifoliiLeandroMMoreira1,2, NalvoFAlmeidaJr13, NehaPotnis3, LucianoADigiampietri12,14, SaidSAdi13, JulioCBortolossi5, AnaCda Silva9, AlineMda Silva2, FabrcioEde Moraes5, JulioCde Oliveira5,6, RobsonFde Souza2, AgdaPFacincani5, AndrLFerraz5, MariaIFerro5, LuizRFurlan7, DanieleFGimenez5, JeffreyBJones3, ElliotWKitajima10, MarceloLLaia5,8, RuiPLeiteJr17, MiltonYNishiyama2, JulioRodrigues Neto11, LetciaANociti5, DavidJNorman18, EricHOstroski14, HaroldoAPereiraJr5, BrianJStaskawicz19, RenataITezza5, JesusAFerro5, BorisAVinatzer4 and JooCSetubal*15,16AbstractBackground: Citrus canker is a disease that ha s severe economic impact on the ci trus industry worldwide. There are three types of canker, called A, B, and C. The three types have different phenotypes and affect different citrus species. The causative agent for type A is Xanthomonas citri subsp. citri whose genome sequence wa s made available in 2002. Xanthomonas fuscans subsp. aurantifolii strain B causes canker B and Xanthomonas fuscans subsp. aurantifolii strain C causes canker C. Results: We have sequenced the genomes of strains B and C to draft status. We have compared their genomic content to X. citri subsp. citri and to other Xanthomonas genomes, with special emphasis on ty pe III secreted effector repertoires. In addition to pthA already known to be present in all three citrus canker strains, two additional effector genes, xopE3 and xopAI are also present in all three strains and are both located on the same putative genomic island. These two effector genes, along with one other effector-like gene in the same region, are thus good candidates for being pathogenicity factors on citrus. Numerous gene content diff erences also exist between the three cankers strains, which can be correlated with their different virulence and host range. Particular attention was placed on the analysis of genes involved in biofilm formation and quorum sensing, ty pe IV secretion, flagellum synthesis and motility, lipopolysacharide synthesis, and on the gene xacPNP which codes for a natriuretic protein. Conclusion: We have uncovered numerous commonalities and diff erences in gene content between the genomes of the pathogenic agents causing citrus canker A, B, and C and other Xanthomonas genomes. Molecular genetics can now be employed to determine the role of these genes in plant-microbe interactions. Th e gained knowledge will be instrumental for improving citrus canker control.BackgroundCitrus canker is a disease with worldwide distribution that has severe economic impact on the citrus industry [1,2]. Disease symptoms consist of water soaked lesions that develop into blisters, then pustules, and, finally, cankers. In severe cases, citrus canker can lead to defoliation and premature fruit drop . Eradication of infected plants is the method of choice to control the disease where it is not yet endemic. When the disease is endemic, control is attempted by planting disease-free trees, limiting the spread between orchards, and using preventive copper sprays [4-6]. However, none of these measures controls citrus canker efficiently. There are three types of citrus canker described in the literature: types A, B and C. Type A originated in Asia, probably in Southern China, Indonesia or India, and it is Correspondence: firstname.lastname@example.org Virginia Bioinformatics Institute, Virginia Poly technic Institute and State University, Blacksburg, VA, USAFull list of author information is available at the end of the article
Moreira et al. BMC Genomics 2010, 11 :238 http://www.biomedcentra l.com/1471-2 164/11/238 Page 2 of 25the type that is most widespread and causes the greatest economic damage [4,7]. The other two types have only been found in South America. Type B (or false canker) was originally identified in Argentina in 1923. This type is present only in Argentina, Paraguay, and Uruguay , whereas type C is limited to the state of So Paulo, Brazil . The causal agent of canker A is Xanthomonas citri subsp. citri (which we abbreviate as XAC for reasons explained below). XAC causes disease on many citrus species, with C. paradisi (grapefruit) and C. aurantifolia (Mexican lime) being most susceptible in the field and C. reticulata (mandarin/tangerine) and C. sinensis (sweet orange) being relatively tolerant [10,11]. Importantly, no citrus species is resistant to XAC after artificial inoculation, suggesting that there is no true genetic resistance against XAC and that field tolerance is mainly due to variation in growth habit . The genome of XAC strain 306 from Brazil was completely sequenced in 2002  and compared to the genomes of Xanthomonas species that are pathogenic in other plants [13,14]. This comparative genomics approach has greatly accelerated the study of the molecular basis of pathogenicity and virulence of XAC. XAC has a hrp/hrc cluster coding for a type III secretion system (T3SS) that is used by the pathogen to inject virulence proteins, called effectors, into host cells. While several genes coding for putative effectors have been identified in the XAC genome, the single most important effector is PthA [3,15]. Even in the absence of the pathogen, PthA induces canker-like symptoms when transiently expressed in plants . Deletion of pthA abolishes the ability of XAC to cause cankers . Intriguingly, PthA induces cankers in all citrus species while it triggers plant immunity in other plant species, thus being the prime determinant of XAC specificity toward citrus [16,18,19]. Two variant forms of canker A have been described. The first was found in Southeast Asia in 1998 infecting C. aurantifolia The pathogen was classified as XAC variant A* . The second variant was isolated in 2003 in Southern Florida in C. aurantifolia and C. macrophyla (alemow), and was named Xanthomonas citri variant AW. AW strains have been shown to be a sub-group within A* . These strains are primarily pathogenic on C. aurantifolia and do not cause disease on C. paradisi even after artificial inoculation [19,20]. A T3SS effector, called AvrGf1, was found to contribute to the exclusion of C. paradisi from the A* host range . A recent study  suggests that A* strains (including AW) have a wider genetic diversity than the strains that cause A-type canker. Canker B is mostly restricted to C. limon (lemon), but has also been found in C. sinensis and in C. paradisi . Its causal agent has been described as X. fuscans subsp. aurantifolii type B (which we abbreviate as XauB for reasons explained below). Even though symptoms are similar to canker A, they take longer to appear, perhaps reflecting the slower growth of XauB in culture when compared to XAC. Canker C has the same symptoms as type A, but, similarly to XAC A* and AW, it is restricted to C. aurantifolia and does not occur in C. paradisi . The causal agent has been described as X. fuscans subsp. aurantifolii type C (which we abbreviate as XauC for reasons explained below). Fig. 1A summarizes observed pheno-Figure 1 Phenotype and Genotype (A) Disease symptoms caused by XAC, XauB, and XauC on leaves of three different citrus species. The lesions caused by XAC and XauB are similar but differ in size. XauC causes a hypersensitive response in C. limonia and C. sinensis with C. aurantifolia as a true host. The pictures were taken 21 days after inoculation. (B) Venn diagram showing the numbers of protei n-coding orthologous genes shared among the three strains.
Moreira et al. BMC Genomics 2010, 11 :238 http://www.biomedcentra l.com/1471-2 164/11/238 Page 3 of 25types in three different citrus species. Recently, what appears to be a new type of X. fuscans subsp. aurantifolii infecting swingle citrumelo ( C. paradisi Macf. Poncirus trifoliata L. Raf.) has been reported in Brazil . XAC, XauB and XauC have been compared phenotypically and analyzed phylogenetically. All three strains present polar flagella with perceptible motility when cultured in semi-solid media . They all grow in the presence of lactose, manitol, and celobiose. However, only XAC is able to grow in the presence of maltose and aspartic acid, and it is also capable of pectate and gel hydrolysis . XauB and XauC have little or no affinity for polyclonal antisera prepared against XAC, and XAC is susceptible to bacteriophages CP1 and CP2 while XauB and XauC are not . It is notable that XauB has fastidious growth in culture media where both XAC and XauC grow well, for example in Agar nutrient and tryptophan-sucrose-agar media. All three grow well in media rich in glutamic acid . Multilocus sequence typing and other molecular analyses [28-30] have shown that XauB and XauC are more closely related to each other than to XAC. Under the rationale that the availability of the genome sequences and annotations of the causative agents of the B and C canker types can substantially improve our understanding of the genomic basis of the disease, we have sequenced the genomes of XauB and XauC to draft status. We have compared them with the genomes of XAC and other xanthomonads. Identified commonalities among the three canker genomes represent candidate genes that may help explain the differences between citrus canker and diseases caused by other xanthomonads. We have also identified numerous gene differences between the three citrus canker genomes. Some of these genes were previously shown to contribute to the virulence of XAC  and are thus primary candidates for explaining the higher virulence of XAC compared to XauB and XauC as well as the host range differences that exist between the three canker types.A note on species abbreviationsThe organisms studied in this work do not have names that are universally accepted. At the time when its genome was sequenced, the accepted name for XAC was Xanthomonas axonopodis subsp. citri In the meantime, Xanthomonas citri subsp. citri has been validly published as a name of this organism [32,33]. However, because the locus tag prefix of XAC genes is 'XAC' we opted for this abbreviation to avoid confusion along the text. Similarly, the causative agents for the B and C cankers are known by different names, but we use XauB and XauC as acronyms so that they are in agreement with their respective locus tag prefixes. For the B species we use the name X. fuscans subsp. aurantifolii type B and for the C species we use the name X. fuscans subsp. aurantifolii type C . We refer to the three organisms collectively as the citrus canker strains (abbreviated by CC strains/genomes ).Results and DiscussionTable 1 shows the genome features of the three canker strains. Even though we do not have complete genome sequences for XauB and XauC, alignments of their scaffolds to the XAC genome suggest that all three chromosomes are highly syntenic (Additional file 1: Fig. S1). Based on these alignments and on total length of contigs we estimate that we have obtained 94% of the XauB genome and 96% of the XauC genome. We estimate that the vast majority of remaining gaps in the two incomplete genomes is under 2 kbp. Fig. 1B shows that XAC shares 74% of its protein-coding genes with the other two strains. The fractions for XauB and XauC are 87% and 84%, respectively. The number of XAC-specific genes is much larger than the anal ogous number for XauB and XauC. Although this difference could be attributed to the incompleteness of the Xau genomes, we have hybridization results (see below) that suggest that many of theseTable 1: General features of XAC, XauB, and XauC genomes.FeatureXACXauB*XauC* genome size (bp) 5,274,1744,877,8085,012,633 # contigs 3239351 %GC 64.764.964.8 protein-coding genes total 4,4273,8043,921 with functional assignment 2,7792,6942,728 conserved hypothetical 1,3869931,009 hypothetical 262117184 RNAs rRNA operons 222 tRNAs 545151 For each strain the numbers are totals for the chromosome plus plasmids. The asterisk (*) denote s that the numbers for XauB and XauC come from their respective incompletely sequenced genomes.
Moreira et al. BMC Genomics 2010, 11 :238 http://www.biomedcentra l.com/1471-2 164/11/238 Page 4 of 25XAC-specific genes are indeed absent from the other two genomes under study. XAC strain 306 has two plasmids, pXAC33 (34 kbp) and pXAC64 (65 kbp). Based on similarities between contig sequences and XAC plasmid sequences it appears that both XauB and XauC have plasmids. At least 46% of plasmid pXAC33 sequence is found in XauB contigs (46% of plasmid pXAC33 sequence is found in XauC contigs); at least 61% of plasmid pXAC64 sequence is found in XauB contigs (55% for XauC). We do not have enough sequence data to ascertain the exact number of plasmids in each Xau genome. The phylogenetic relationship of the three strains with respect to each other and to fully sequenced members of the Xanthomonas and Xylella genera based on their shared protein-coding genes is shown in Fig. 2. The tree shows that the three organisms under study form a welldefined group within the family Xanthomonadaceae This tree is in agreement with trees of the Xanthomonas genus based on multilocus sequence analysis [34,35]. At http://bioinfo.fac om.ufms.br/xanthomonas we provide an interactive tool that allows user-defined gene content comparisons among all sequenced Xanthomonas genomes.Genes shared by all three strains but not present in other Xanthomonadaceae speciesXAC, XauB, and XauC have 65 families of orthologous genes specific to them when compared to all other fully sequenced members of the genera Xanthomonas and Xylella (Additional file 2: Table S2). Among these 65 families we have identified 11 syntenic blocks. Not surprisingly, almost half of these genes code for hypothetical proteins of unknown function. Of the genes with a predicted function the genes encoding the predicted effectors XopE3 (XAC3224) and XopAI (XAC3230) (discussed below) are especially noteworthy. The large number of genes coding for various kinds of transporters are also worth mentioning: four TonB dependent receptors including one with homology to the Escherichia coli receptor FepA, which is involved in transport of siderophores across the bacterial membrane  and seven ABC transporters, which might be used either for translocation of substrates from the citrus apoplast into the bacterial cell to provide nutrients for the pathogen, or, alternatively, for secretion of toxins (either bacterial toxins or expulsion of citrus metabolites toxic to the CC strains). An additional transporter specific to the three CC genomes (XAC3198) is an alkanesulfonate transporter substrate-binding subunit, which is reported to enable E. coli to use sulfonates other than taurine . Besides transporters, several genes encode metabolic enzymes (an amidase, an urea amidolyase, a peptidase, and a nitrilotriacetate monooxygenase). This conspicuous presence of transporters and metabolic enzymes suggests that the CC strains might have adapted to specific metabolites present in the citrus apoplast. However, this will need to be confirmed experimentally by comparing growth of wild type strains and strains mutated in CCspecific genes in apoplastic fluid of citrus species and of other plant species. One other interesting gene present in all CC genomes is a gene co ding for a methyl parathion hydrolase (XAC0726), predicted to degrade the insecticide methyl parathion [38,39]. Orthologs of this gene and other genes that degrade organophosphates are common in soil bacteria and in the soil-borne pathogen Ralstonia solanacearum but have not yet been found in any other foliar plant pathogen. Some of these syntenic CC-specific regions are anomalous in terms of nucleotide composition as determined by the program AlienHunter  and may thus have been acquired by horizontal gene transfer.The three CC genomes have important differences in regard to their repertoires of type III secreted effectorsThe hrp/hrc genes encoding the T3SS are basically the same and found in the same order in all three CC genomes. However, there are notable differences in the three putative T3SS-secreted effector repertoires. A list of twenty-seven T3SS effector genes predicted in the genomes of the CC strains is shown in Table 2. Effectors are important determinants of virulence and host range in many plant pathogenic bacteria, in particular in Xanthomonas sp. and Pseudomonas syringae . Comparison of effector repertoires between the three CC genomes and all other Xanthomonas genomes can thus give us important clues. The effector genes avrBs2, xopL, xopQ and xopX are present in all three CC genomes, in all sequenced genomes of other Xanthomonas species, and in all X. citri and most Xanthomonas strains that were surveyed by PCR and hybridization for these genes by Hajri et al. . These effectors thus belong to the Xanthomonas core set of effectors possibly important for pathogenicity on all plants. The putative effector genes xopK, xopR and xopZ also belong to this group since they can be found in all sequenced Xanthomonas genomes. However, no data exist for these effectors in regard to other Xanthomonas strains . The effector genes xopI, xopV, xopAD and xopAK are present in all three CC genomes and in several, but not all, sequenced Xanthomonas genomes. These effectors, therefore, might contribute to disease in some plant species while they might trigger immunity in others. As already mentioned, PthA is well known to be an important X. citri effector that plays an essential role in citrus canker, while limiting the host range of CC strains to citrus because it triggers immunity in all other tested plant species (see references above). The pthA gene is a
Moreira et al. BMC Genomics 2010, 11 :238 http://www.biomedcentra l.com/1471-2 164/11/238 Page 5 of 25 Figure 2 Maximum likelihood tree of fully sequenced Xanthomonas and Xylella genomes The bootstrap support is 100% for all branches (100 bootstrap runs). Bar, number of amino acid substitutions per site.
Moreira et al. BMC Genomics 2010, 11 :238 http://www.biomedcentr al.com/1471-2164/11/238 Page 6 of 25Table 2: Putative effectors found in the XAC, XauB, and XauC genome sequences.Effector familyXACXauBXauCPfam: functional/structural domain Effectors Candidate effectors common to XAC, XauB, and XAUC AvrBs2 XAC0076XAUB_16770XAUC_23650G lycerophosphoryl diester phosphodiesterase AvrBs2 from X. campestris pv. vesicatoria  AvrBs3 XACa0022 (pthA1) XACa0039 (pthA2) XACb0015 (pthA3) XACb0065 (pthA4) XAUB_40130 XAUB_28490 XAUC_22430 XAUC_24060 XAUC_09900 XAUC_43080 Transcriptional activator, nuclear localization PthA  XopE1 (avrXacE1, hopX, avrPphE) XAC0286XAUB_37010XAUC_37580Putative transglutaminaseAvrXacE1, XopE1 from X. campestris pv. vesicatoria  XopE3 (avrXacE2, hopX, avrPphE) XAC3224XAUB_14680XAUC_00040Putative transglutaminaseAvrXacE2  XopF2 XAC2785 XAUB_07540 XAUC_21000 XopF2  XopI XAC0754XAUB_39080XAUC_07100F-box protein X. campestris pv. vesicatoria  XopK XAC3085XAUB_34090XAUC_12520Identified in Xoo by cya assay  XopL XAC3090XAUB_34130XAUC_02900/12488 LRR protein XopQ (hopQ1) XAC4333XAUB_10220XAUC_14670Inosine uridine nucleoside Nribohydrolase  XopR XAC0277XAUB_36920XAUC_37490Identified in Xoo by cya assay  XopV XAC0601XAUB_23140XAUC_21260Identified in Xoo by cya assay  XopX (HolPsyAE) XAC0543XAUB_14760XAUC_20690 XopZ (HopAS, AWR) XAC2009XAUB_11532/13710 XAUC_25915 XopAD (skwp, RSc3401) XAC4213XAUB_02510XAUC_34870SKWP repeat proteinSkwp from Ralstonia 
Moreira et al. BMC Genomics 2010, 11 :238 http://www.biomedcentr al.com/1471-2164/11/238 Page 7 of 25XopAI (HopO1 (HopPtoO, HopPtoS), HopAI1 (HolPtoAI)) XAC3230XAUB_26830XAUC_23780ADP-ribosyltransferase XopAK (HopAK1 (HopPtoK, HolPtoAB)C terminal domain) XAC3666XAUB_02580XAUC_32490Not confirmed to be effector in Xanthomonas ; homolog of effector in Pseudomonas HrpW (PopW) XAC2922XAUB_19460 (associated with hrp cluster) XAUC_20020 (associated with hrp cluster) Pectate lyase, may not be T3SE Candidate effectors present in XAC and XauB BUT ABSENT in XauC XopE2 (avrXacE3, avrXccE1) XACb0011XAUB_31660-Putative transglutaminaseXopE2 found in another C strain  XopN (hopAU1) XAC2786XAUB_07520-ARM/HEAT repeat XopP XAC1208XAUB_06720- XopAE (HpaF/G/PopC) XAC0393XAUB_19500-LRR proteinXcv8510  Candidate effectors present in XauB and XauC BUT ABSENT from XAC XopB (hopD1, avrPphD1) -XAUB_09070/14842 XAUC_00260 XopE4 (HopX) -XAUB_23330XAUC_31730New class introduced XopJ (AvrXccB) -XAUB_20830XAUC_08850C55-family cysteine protease or Ser/Thr acetyltransferase  XopAF (avrXv3, HopAF1 (HopPtoJ)) -XAUB_02310XAUC_00300 XopAG (AvrGf1, HopG1 (HopPtoG). HolPtoW) -XAUB_03570 XAUC_04910AvrGf1  Candidate effectors present only in XauC XopF1 (Hpa4) --XAUC_20060  Some of the effectors appear to be pseudogenes (indicated by a ). In column 'effector family' we provide the standard effector name along with the names of related proteins belonging to the same family; in column 'effectors' we provide the reference where the effector indicated was characterized.Table 2: Putative effectors found in the XAC, XauB, and XauC genome sequences. (Continued)
Moreira et al. BMC Genomics 2010, 11 :238 http://www.biomedcentra l.com/1471-2 164/11/238 Page 8 of 25member of the avrBS3 family of effector genes, members of which are present in most Xanthomonas genomes and in some R. solanacearum genomes . However, only PthA is known to induce citrus canker. Besides pthA (XACb0065), three paralogs of pthA are also present in the XAC genome (XACa0022, XACa0039, and XACb0015). All four copies are found on plasmids. The three paralogs do not seem to play an important role in citrus canker . We found two pthA homologs in the XauB genome (XAUB_40130 and XAUB_28490) and two in the XauC genome (XAUC_22430 and XAUC_24060/ XAUC_09900 [the latter is a single gene with halves in different contigs]). Not all of these genes have been completely assembled due to the repetitive regions found in avrBS3 family members. However, El Yacoubi et al.  previously assembled a pthA homolog ( pthB [GenBank: 2657482 ]) from the pXcB plasmid [GenBank: NC_005240 ] of a XauB strain with the same repeat copy number (i.e. 17.5) as pthA and Al-Saadi et al.  sequenced and assembled another homolog ( pthC [GenBank: EF473088 ]) from a XauC strain. These genes functionally complemented a pthA deletion in XAC without affecting host range . The XAUC_22430 gene has 99% nucleotide identity to pthC and thus probably corresponds to pthC and would be the functional pthA homolog of XauC. We do not have enough data to confidently report on the repeat copy number of the other three Xau pth homologs, but a phylogenetic analysis (see below) suggests that XAUB_28490 is the functional pthA homolog of XauB.To get insight into the evolution of the pthA/avrBS3 family, a phylogenetic analysis of all available avrBS3 family members in the genus Xanthomonas was performed (Fig. 3). The tree shows that the four XAC pth genes group together, while the XauB and XauC genes form two families, with one member from each strain participating in each family. Moreover, the XauB and XauC pth genes group separately from the XAC genes, with good bootstrap support, which is a surprising result when compared with the species tree (Fig. 2). The genes that flank the pth copies in XauB do not have XAC genes as their best BLAST  hits. In particular, the two genes upstream of XAUB_40130, XAUB_40110 and XAUB_40120, do not match any Xanthomonas genes; instead, their best BLAST hits (e-value 10-102 and 10-99, respectively) are gene sequences from Burkholderia pseudomallei NCTC 13177. (In XauC the pth regions are too fragmented to allow genomic context determination.) We derive the following conclusions from these results: An ancestor of the current XauB and XauC strains already had two different copies of the pth genes; hence the existence of the two noted families. The fact that in XAC the four copies are nearly identical, the fact that the Xau and XAC pth genes group in a way distinctively different from the species grouping, and the fact that the genomic context in which the XAC pth genes are found is different from that of XauB leads us to believe that XAC acquired its pth genes by a different route and that their duplication is more recent when compared to the XauB and XauC pth duplications. The ability to cause citrus canker may thus have independently evolved by XAC and by XauB and XauC, and may represent an example of convergent evolution. Al-Saadi et al. , based on a neighbor-joining phylogeny of pth genes that included the pthB and pthC sequences, have come to the same conclusion.Effectors XopAI and XopE3 may play a role in citrus cankerA comparison of effectors present in all three CC strains with those present in fully sequenced Xanthomonas species, and data from the study by Hajri et al. , suggest that two additional putative effectors may play a special role in citrus canker. These are XopAI and XopE3. Both are present in all three CC genomes. The putative effector xopAI is not found in any other sequenced Xanthomonas species and it was not included in the Hajri et al.  analysis. We do have evidence that it is present in Xanthomonas vesicatoria str. 1111 (Potnis et al., unpublished). Interestingly, the C-terminal region of XopAI has similarity to predicted ADP-ribosyl transferase domains of the effector HopO1-1 of Pseudomonas syringae and of hypothetical proteins in Acidovorax citrulli Ralstonia solanacearum and other bacteria. The Nterminus has high similarity to the N-terminus of the effector XopE2 of X. campestris pv. vesicatoria 85-10 as well the N-termini of a number of other Xanthomonas and Pseudomonas syringae effectors (more on the N-terminal region of xopAI below). XopE3 belongs to the HopX/AvrPphE family of effectors. Effectors belonging to this family have been found in diverse phytopathogenic bacteria including Ralstonia, Pseudomonas, Acidovorax and Xanthomonas suggesting their conserved role in virulence on a wide range of hosts. Sequences from this family have similarity to the transglutaminase superfamily of enzymes, which are responsible for modification of host proteins . The HopX/ AvrPphE effector from Pseudomonas syringae has been shown to be involved in host protein proteolysis, thereby suppressing host defenses [46,47]. In xanthomonads, multiple effectors belonging to this group have been found, such as xopE1, xopE2, xopE3, xopE4 XopE1 and xopE2 have been found in most of the xanthomonads. XopE3 effector gene homologs have been found by PCR and dot-blot hybridization methods in some Xanthomonas axonopodis strains belonging to the alfalfae anacardii glycines phaseoli malvacearum fuscans mangiferae, indicae and citrumelo pathovars . However, sequences of xopE3 from these strains could not be compared against homologs from CC strains since sequence
Moreira et al. BMC Genomics 2010, 11 :238 http://www.biomedcentra l.com/1471-2 164/11/238 Page 9 of 25 Figure 3 Maximum likelihood tree of Xanthomonas pth genes The sequences used to build the tree had their tandem repeat portions masked before alignment. The root position was obtained from a preliminary tree that included as an outgroup a Ralstonia solanacearum pth gene [GenBank: CAD15517.1 ]. Full information about the pth gene sequences used is given in additional file 6 (Table S6). The numbers on branches represent bootstrap support (from 100 bootstrap runs). The bootstrap values obtained fo r the clade indicated by the asterisk are given in a dditional file 7 (Fig. S7). Bar, number of amino acid substitutions per site.
Moreira et al. BMC Genomics 2010, 11 :238 http://www.biomedcentra l.com/1471-2 164/11/238 Page 10 of 25data from the X. axonopodis strains mentioned are not currently available. Phylogenetic analysis of hopX orthologs shows that the xopE3 effector genes found in the CC strains group together with hopX1 effector genes from pseudomonads (data not shown). Although all hopX orthologs show conservation of the catalytic triad (Cys, His, Asp residues) as well as the conserved domain "GRGN" N-terminal to the triad, the region C-terminal to the triad shows high degree of variability. This variable region has been hypothesized to be responsible for targeting different host proteins . In fact, while some AvrPphE (hopX) alleles from P. syringae pv. phaseolicola strains trigger gene for gene disease resistance in some bean cultivars, other alleles were shown to be virulent on these same cultivars. Amino acid differences in the C-terminal region of AvrPphE were identified between alleles . Similarly, comparing XopE3 homologs from different strains at the amino acid level and their corresponding reactions on different hosts might give clues regarding the variable C-terminal domains of XopE3 family members and might determine whether this variability is responsible for targeting different proteins in different host species. Both xopE3 and xopAI belong to an interesting XAC chromosomal region of approximately 15 kbp in size (Fig. 4) that has been hypothesized to be a genomic island . An alignment of the XAC chromosome sequence with the chromosome sequences of X. campestris pv. campestris str. ATCC33913 and X. oryzae pv. oryzae str. PX099A strongly suggests that this region is an insertion (data not shown). The presence of three transposase genes and two phage-related genes in the region provides additional evidence for this hypothesis. The central part of this region (7 kbp) duplicates a region found in XAC plasmid pXAC64 (Fig. 4), suggesting a chromosome-plasmid DNA exchange. In the plasmid we find the effector gene xopE2 (XACb0011), which as described above shares its N-terminal region with xopAI (XAC3230) (Fig. 4). Transposons and phage elements in this region might thus have been responsible for a shuffling process, described as terminal reassortment , resulting in the novel effector gene xopAI Although we can characterize this region completely only in XAC, XauB and XauC contigs contain the most important elements of this region (Fig. 4). Next to xopE3 (XAC3224) we find gene XAC3225, whose product is annotated as tranglycosylase mltB This gene has strong similarity (e-value 10-133, 100% coverage) to hopAJ1 from P. syringae pv. tomato strain DC3000, where it is annotated as a T3SS helper protein. Although the hopAJ1 gene is not itself a T3SS substrate, it contributes to effector translocation . A mutant with a deletion of XAC3225 has reduced ability to cause canker (mutant phenotypes include a reduction in water soaking, hyperplasia, and necrosis compared to wildFigure 4 Region containing citrus-can ker specific effector genes xopE3 and xopAI The XAC region depicted is hypothesized to be a genomic island . The central parts of this region (gray areas) are CCspecific, being present in both Xa uB and XauC as well as in X AC plasmid pXAC64. Gene mltB is not an effector but plays a role in Type III secretion.
Moreira et al. BMC Genomics 2010, 11 :238 http://www.biomedcentra l.com/1471-2 164/11/238 Page 11 of 25type) . We thus conclude that the effector and effector-related genes in this region probably play an important role in citrus canker.Additional differences in effector repertoires among CC genomesIn addition to the pth differences noted above, other effectors that distinguish the XAC genome from the two Xau genomes are XopB, XopE4, XopJ (AvrXccB), XopAF (avrXv3) and XopAG which are all present in both XauFigure 5 Genomotyping and similarity an alysis between the genomes of XAC and XauC ba sed on DNA hybridization and matching of contig sequences (A) The hybridization results showed that 2,486 CDSs (out of a tota l 2,760) gave a hybridization sign al greater than the estimated background noise. Of the 2,486 CDSs, 2,341 (94. 2%) seem to be present in XauC (ratio Cy3/Cy5 between 0 and 1.5). The remaining 145 CDSs (5.8%) seem specific to XAC. (B) Of these 145, most (101 = 70%) belong to regi ons previously described as putative genomic islands in XAC (black ovals: ; gr een ovals: ; yellow ovals: ; the oval numbering corresponds to the original publication numberin g). Pink ovals mark other nonsimilar regions between the genomes of XAC and XauC. (C) The two bottommost graphs (G+C and codon bias) sh ow variation of these tw o metrics along the XAC genome, thus presenting evidence for the pu tative genomic islands denote d by ovals. The numbered triangles correspond to the XA C-specific regions based on DNA hybridization. Blue horizontal bars simply denote regions in these graphs that correspond to ovals not assoc iated with XAC-specific regions. The Heat Map shows gene grou ps that yielded differential hybridization signals, the va st majority of which corre spond to regions marked by ovals. The zoomed-in regions in the upper diagram show genes with differentia l hybridization signals bu t that are in regions shared by the two genomes. The XAC and XauC genes in these regions have less simila rity (from 42 to 68% identity) between them than other shared genes.
Moreira et al. BMC Genomics 2010, 11 :238 http://www.biomedcentra l.com/1471-2 164/11/238 Page 12 of 25genomes but absent from XAC strain 306. (AvrXccB homologs were found in two XAC strains by Hajri et al. .) The absence of these effectors from XAC strain 306 raises the possibility that these effectors might be responsible for limiting the host range of both B and C strains. Interestingly, XauB and XauC strains both contain xopAG an effector gene belonging to the same effector family as avrGf1 from X. citri Aw, which has been shown to be responsible for triggering a hypersensitive defense response in C. paradisi (grapefruit) . The xopAG gene from the B and C genomes shows 44% identity to avrGf1 at the amino acid level. The XauB and and XauC genes are almost identical to each other, with one important difference: in XauB xopAG is interrupted by a transposon. Therefore, the incompatibility between XauC and grapefruit and the ability of XauB to cause disease in grapefruit could be explained by this single gene difference. The xopE4 DNA sequence is identical in the two Xau genomes and has similarity to avrXacE3 but only with 31% identity at the amino acid level; this is why we named this gene xopE4 instead of xopE2 Unlike other XopE family members, XopE4 does not have a predicted myristoylation site, suggesting that it may not be targeted to the cell membrane as the other XopE family members. Presence of an additional effector gene, the avirulence gene avrXccA2 has been shown in some X. aurantifolii B (CFBP3528, CFBP3530) and X. aurantifolii C (CFBP2866) strains by hybridization and PCR analysis . However, this avirulence gene was not found in the two sequenced Xau genomes. A homolog of the effector xopF1 (XAUC_20070) was found only in the XauC strain. It is located in a 5-kbp region that lies between the T3SS genes hrpW (XAUC_20020) and hpa3 (XAUC_20080). The same two genes are adjacent in XauB. Two transposases are present in this region, and the sequence of xopF1 has a frameshift, suggesting that this gene is likely the result of a recent insertion and is not active. There are four effector genes present in the XAC and XauB genomes that have not been found in the genome of XauC: xopE2, xopN, xopP and xopAE These effectors could explain the wider host range of XAC and XauB compared to XauC, assuming a virulence activity of these effectors on citrus species. XopN has been shown to interact with the plant protein TARK1 and to interfere with immunity triggered by pathogen-associated molecular patterns (PAMP-triggered immunity) . Further experiments are required to determine the possible role of XopN in extending host range to lemon, grapefruit and sweet orange. Another effector that could have a similar role is XopAE (a hpaF/PopC homolog) [52,53]. The harpin-like protein HrpW with a pectate lyase domain is present in all CC strains. In the sequenced XAC genome, it is not associated with the T3SS gene cluster, whereas in the genomes of XauB and XauC it is. The role of harpin-like proteins like HrpW as virulence factors or T3SS accessory proteins has not yet been determined in the Xanthomonas genus. Experiments will need to be performed to confirm translocation of the above putative effectors and their putative function as virulence or avirulence genes.XAC-specific genes and genomic regions with respect to XauB and XauCWe compared the XAC genome to the XauB and XauC genomes both computationally and by doing experimental whole genome hybridizations. Results are summarized in Fig. 5. The results of the two methods were consistent. In the following sections we focus on regions and genes specific to XAC with respect to XauB and XauC according to these results. We have identified 25 groups of at least four consecutive genes that we term XAC-specific regions (XACSR) [Fig. 5, additional file 3 (Table S3), and additional file 4 (Fig. S4)]. Nearly all regions contain or are flanked by transposition elements or phage-related genes, suggesting that they could be the result of lateral transfer.Several genomic differences are related to biofilm formation and quorum sensingXanthan gum is an exopolysaccharide that plays an important role in biofilm formation and hence in virulence of pathogens of the Xanthomonadaceae family . Moreover, the synthesis of xanthan gum is regulated by variation in sugar concentration in the culture medium and by the activation of regulatory rpf genes [57,58]. These genes are also responsible for the synthesis of diffusible signal factors, fundamental molecules for quorum sensing processes [59,60]. Both Xau genomes contain an identical xanthan gum operon (XauC: XAUC_26940-27060; XauB: XAUB_007400-007410 and XAUB_10560-10450). The Xau genomes contain gene rpfH (XAUB_10500 and XAUC_27010), but this gene is not found in XAC. Gene rpfI present in the xanthan gum operon of X. campestris pv. campestris strain ATCC33913, is absent from all three CC genomes. XAC and the two Xau genomes contain the xanthomonadin biosynthesis regulon as well as sugar metabolism genes. XauB however presents some important variations that may explain its in planta and in culture fastidious phenotype, when compared to both XAC and XauC. Differences were found in the phosphotransferase system (PTS-Fru), which specializes in internalization of fructose, and in the rpfN gene, a sugar porin, which is regulated by rpf genes, which also regulate xanthan gum synthesis . The importance of the PTS-Fru system and of the sugar porin encoded by the rpfN gene for growth and pathogenicity of certain bacteria has been reported in the literature: PTS-Fru mutants of Spiroplasma citri caus-
Moreira et al. BMC Genomics 2010, 11 :238 http://www.biomedcentra l.com/1471-2 164/11/238 Page 13 of 25ative agent of citrus stubborn disease, are reduced in virulence [61-64], and rpfN mutants of Xanthomonas campestris pv. campestris show an increase in the level of polygalacturonate lyase , an important pathogenicity factor in bacterial plant pathogens. This result shows that the absence of the sugar porin could cause lack of carbohydrate uptake, therefore inducing the synthesis of cell wall degrading enzymes, in order to increase sugar supply [58,66]. In XAC the PTS system, encoded by genes fruBKA and rpfN is organized in one single genomic region (XAC2501-2504) (Fig. 6). XauC presents the same organization (XAUC_01750-01780). In XauB fruBKA corresponds to XAUB_05120, XAUB_05110 and XAUB_05100 respectively, but rpfN was not found. In addition, the fruA gene sequence contains a frameshift, indicating that it may have become a pseudogene. XauBFigure 6 Models for sugar acquisition, quorum sensing modulation, and biofilm formation and the XAC-specific genes and regions related to these processes Region XACSR17 (lower right corner) contains four genes that when mutated (separately) cause decrease in biofilm activity and virulence, as indicated in the leaf photos. Th e photo labels are mutant strain identifiers.
Moreira et al. BMC Genomics 2010, 11 :238 http://www.biomedcentra l.com/1471-2 164/11/238 Page 14 of 25needs a culture medium with glutamic acid , possibly because it is used by the bacterium as an alternative carbon source. Xylella fastidiosa 9a5c, which also lacks the PTS system and the rpfN gene, is also fastidious . Experiments were carried out to compare cellular growth and xanthan gum production in all three organisms under study. Results (Table 3) confirm the fastidiousness of XauB, whereas XAC and XauC have similar cellular mass values. We also determined that XAC produces more than twice as much gum as XauC and almost three times as much as XauB. We believe this variation is due to the XauB deficiencies in the rpfN and PTS-Fru genes discussed above. On the other hand, this result also suggests that the lack of the rpfH and rpfI genes in XAC, caused by the presence of transposition elements , does not influence the gum production capability of XAC. The fastidiousness of XauB would thus be related to the low gum production and the consequent decrease in biofilm formation together with the fact that it needs an additional substrate source (glutamate) to support its growth in culture. In addition to the differences noted above, five additional XAC-specific regions (XACSR7, XACSR9, XACSR10, XACSR14 and XACSR17) may be related to its greater biofilm-formation capability when compared to XauB and XauC. Several of the genes in these regions facilitate adhesion in a process mediated by hemagglutinin . XACSR7 contains two hemin storage system genes, hmsF and hmsH and hemagglutinin coding genes (XAC1811-1816). Genes involved with acquisition and storage of hemin groups ( hmsRFH ) and type I secretion system genes ( fhaC ), and their secreted hemagglutinin ( fhaB ), are found in tandem and flanked by a tRNAR in the XAC genome (Fig. 6). In Yersinia pestis the hms genes are present in a cluster (the pgm cluster) related to temperature-dependent storage of hemin as well as expression of a number of other physiological characteristics . Mutations in these genes cause drastic decrease in Yersinia growth, preventing it from colonizing its point of entrance in infected flies (bucal orifice) [71,72]. These genes also play a role in exopolysacharide synthesis, and reduction in biolfilm formation has also been observed in these mutants in Yersinia . Among all sequenced xanthomonads only XAC and X. oryzae have these genes. In both cases they are similar to (35 to 53% identity at the amino acid level) and syntenic with their homologs in Yersinia, E. coli K-12 and Erwinia carotovora (data not shown). Recent work  describing mutations in the genes that code for hemagglutinin in Xylella fastidiosa strain Temecula (Pierce's disease) has shown that biofilm composition and virulence (adhesion and colonization) were affected in the mutants. This is consistent with results in XAC . This is evidence that the apparent absence of these genes in XauB and XauC might have the same effect (Fig. 6). XACSR9 contains 19 genes. One of them (XAC1918) is a hemolysin-related gene. In enterobacteria hemolysins are an important virulence factor that are associated with proteins related to biofilm formation [75,76]. The protein encoded by this particular gene interacts with VirD4, a Type IV secretion system component , which in turn may play a role in biofilm formation and cell aggregation, as observed for E. coli . XACSR10 contains several noteworthy genes. Gene XAC2151 codes for the YapH protein. Its homolog (XOO2380, 84% identity at the amino acid level) in X. oryzae pv. oryzae KACC10331 when mutated drastically reduced the pathogen adhesion to plant tissue, thus decreasing its virulence ; in addition a homolog of this gene in X. fuscans subsp. fuscans CFBP4834-R (the causative agent of bacterial blight of bean, Phaseolus vulgaris ), was required for adhesion to seed, leaves, and abiotic surfaces . XAC2197 and XAC2198 code for hemolysin-type calcium binding proteins, whereas XAC2201 and XAC2202 code for hemolysin secretion protein D (HlyD) and hemolysin secretion protein B (HlyB), respectively. The latter four genes do not have Xanthomonas matches in the sequence databases; they are similar instead to protein sequences from Acidovorax Xylella and Pseudomonas species. The best hits are from Acidovorax avenae subsp. avenae ATCC19860, which is also a plant pathogen. XACSR14 contains genes related to the type IV pilusdependent system (Fig. 6). This system takes part in several processes, including adhesion, motility, microcolony formation, and protease secretion . In Xylella fastidiosa functional studies of these genes have shown that they are crucial for the host colonization process [82-86]. XACSR17 is also related to biofilm formation. Laia et al.  observed decrease in biofilm activity and virulence in four mutants with changes in this region (XAC324514G01/14G12, XAC3263-10G07/10G09, XAC328510F02 and XAC3294-17B04) (Fig. 6). The only mutated gene with functional assignment is rhsD (XAC3245). This gene has been described as coding for a membrane protein related to adhesion . In Xanthomonas campestrisTable 3: Xanthan gum production.StrainPellet weight mg/ml (sd) Gum production mg/ml (sd) XAC2,57 ( 0,25)3,93 ( 0,83) XauB1,70 ( 0,20)1,62 ( 0,50) XauC2,35 ( 0,92)1,82 ( 0,33) sd Standard deviation
Moreira et al. BMC Genomics 2010, 11 :238 http://www.biomedcentra l.com/1471-2 164/11/238 Page 15 of 25pv. campestris a RhsD protein has been found in the outer membrane vesicle associated with other virulence-associated proteins, such as HrpA/F/X/B4, HrcU, AvrBs1 and AvrBs2 . XAC contains a paralog of this gene (XAC2529), but it was not found in the Xau genomes either. Gene wapA (XAC1305), not part of any XAC-specific region, and which is an adhesion facilitator by way of hemagglutinins , was not found in either of the Xau genomes.XauB contains T4SS gene clusters similar to those found in Ralstonia solanacearum and in Agrobacterium tumefaciensIn addition to the type III secretion system, the type IV secretion system (T4SS) also plays a role in pathogenicity. For example, this system has been shown to contribute to full virulence in the phytopathogen X. campestris pv. campestris strain 8004 . Genome clusters containing T4SS genes are found in several bacterial species but with marked differences in terms of gene presence/absence and organization, which relate to the system's function in the organism where it is present [89,90]. All three genomes under study contain T4SS genes, arrayed in clusters. In order to better understand these gene clusters we have classified the T4SS genes found in the three CC genomes as well as T4SS gene clusters from other organisms into four groups, using as criteria presence/absence of genes and synteny (Fig. 7). We have found important differences between XAC and XauB using this classification; data for XauC was too fragmented to allow a similar general observation. As reported by da Silva et al.  XAC has two T4SS clusters, one in the chromosome and the other in a plasmid (pXAC64). In XauB we found three clusters. Only one of them is similar to a XAC cluster (the pXAC64 cluster), and was therefore placed in group II, along with the T4SS cluster found on plasmid pXcB (already mentioned above in the context of the pth gene discussion). The other two XauB clusters were placed in groups I and III, respectively, while the XAC chromosome cluster was placed in group IV (Fig. 7A). The XauB cluster placed in group III is similar to a cluster found in Ralstonia solanacearum , both in terms of organization as well as individual gene sequence simi-Figure 7 Schematic representation of Type IV Secretion System Gene s in diverse bacteria, grouped by their genomic architecture The colors used for gene arrows correspond to the colors used for their respective protein products in the secretion system representa tion on the right. At: Agrobacterium tumefaciens C58; Lp: Legionella pneumophila ; Xfa: Xylella fastidiosa 9a5c; Rs: Ralstonia solanacearum ; Xcc : Xanthomonas campestris pv. campestris ATCC33913; Xoo: Xanthomonas oryzae pv. oryzae References not cited in the te xt are: Goodner et al. (2001) ; Segal et al. (1999) ; Cazalet et al. (2004) ; Bolland et al. (1990)  ; Marques et al. (2001)  ; Lee et al. (2005) .
Moreira et al. BMC Genomics 2010, 11 :238 http://www.biomedcentra l.com/1471-2 164/11/238 Page 16 of 25larity. This organization is quite different from those found in other bacterial species. This XauB cluster is found in a region containing 45 genes, all of which are XauB-specific when compared to other Xanthomonas genomes. Moreover this cluster is flanked by insertion elements. This evidence suggests that this cluster was likely acquired by lateral transfer. The third XauB cluster was placed in group I, which contains T4SS clusters similar to those found in Agrobacterium tumefaciens C58 plasmid Ti , Xylella fastidiosa 9a5c plasmid pXF51 , and rhizosphere plasmids pIPO2T and pSB102 [94,95]. This XauB cluster is most similar to those found in the rhizosphere plasmids, in particular to pIPO2T. The XauC draft genome sequence did allow the identification of one T4SS gene cluster, and it belongs to this group (Fig. 7A). The XAC chromosomal cluster belongs to group IV (Fig. 7A). It is flanked by regions XACSR12 and XACSR13 (Fig. 8A). In addition, some of the genes in this T4SS cluster had uncertain hybridizations (zoomed-in portion of Fig. 5). This adds evidence that this cluster, as a specific unit, is indeed absent from both Xau genomes. XACSR12 contains genes that code for transporters, alpha-glucosidases, hypothetical genes and transposases, as well as a virB6 -related gene. XACSR13 contains only hypothetical genes and transposases. One of the members of this T4SS cluster is virD4 (XAC2623). Alegria et al.  have identified 12 XAC proteins that interact with VirD4. Six of these (XAC0096, 3266, 0151, 4264, 2609, 1918) are apparently absent from the Xau genomes, and five belong to XAC-specific regions (Fig. 8B). These genes might thus contribute to the increased virulence of XAC when compared to XauB and XauC.Flagellum and motilityThe XAC genome contains all genes known to be needed for flagellum synthesis [12-14]. Three of these clustersFigure 8 The chromosomal T4SS gene cluster in XAC (A) This T4SS cluster is located between regions XACSR12 and XACSR13. Genes shared with XauB and XauC are shown with a black asterisk. Two of the genes are shared only by XAC and XauB (red asterisk). (B) Interactions among T4SS proteins, based on data presented by Alegria et al. . Proteins specific to XAC are represented in yellow, and proteins shared by XAC, XauB and XauC are represented in oran ge. Six of the proteins colored in yellow are in XACspecific regions. Genes XAC0095 and XAC0096, although n ot part of XAC-specific regions, seem to play an especially important role in pa thogenicity. Under this model th e protein coded for by XAC0095 in teracts with HrpG, a protein that participates in the T3SS apparatus , and that causes phenotypic changes when its gene is mutated . The pro duct of gene XAC0096, which is next to XAC0095, interacts wirth VirD4, a component of the T4SS .
Moreira et al. BMC Genomics 2010, 11 :238 http://www.biomedcentra l.com/1471-2 164/11/238 Page 17 of 25(which we call F1, F2, and F3) are close to the terminus of replication and are spread ou t over a 126 kbp region (Fig. 9A). A fourth cluster contains just two genes, coding for flagellar motor proteins A and B ( motA and motB ). All four clusters are found in both XauB and XauC with one exception: in XauB we did not find cluster F2 (XAC1930XAC1955, 31 kbp), which contains 26 genes, coding for flagellum proteins that interact with bacterial membranes and coding for flagellar motor components (Fig. 9CD). On the other hand we have observed (microscopically) that XauB does have a flagellar structure similar to those of XAC and XauC, and motility tests have shown that XauB is able to move (Fig. 9B). Although our genome sequence data for XauB is incomplete, we believe it unlikely that a 31 kbp long fragment would be entirely missed; moreover, the hybridization experiment did notFigure 9 Genes related to flagellum synthesis and regulation (A) Synteny and hybridization results between XAC e XauB/XauC. In the upper left diagram, F1, F2 and F3 denote gene regions related to flagellar functions, whereas ge ne groups XACSR9 and Y appear to contain unrelated genes. Dashed boxes denote regions apparently absent from the XauB and XauC genomes. Horizontal green bars represent hybridization sig nals for probes marked with Cy3 and which are XAC-specific; yellow horizontal ba rs represent regions with positi ve hybridization results. The z oomed-in region XACSR9 shows that it contains four transposases (in blue), twel ve hypothetical genes (in gray), and one gene with assigned func tion (in yellow). That gene (XAC1927) codes for a Fe-S-oxidoreductase that contains a RA DICAL SAM domain. XAC mutants in this particular region (14H02 and 25D11 ) presented decreased virulence phenotypes. (B) Microscopy and motility test of each organism. Th ese results validate flagellum presence and functionality in two different culture media (DYGS and TSA). (C) Flaggellum gene organization in gene clusters F1, F2 and F3, as given in panel A. The color coding is the same for panel D. (D) Representation of the flagellum proteins using the same color coding as in panel C.
Moreira et al. BMC Genomics 2010, 11 :238 http://www.biomedcentra l.com/1471-2 164/11/238 Page 18 of 25indicate presence of this fragment. These results suggest that the absence of genes in cluster F2 is compensated in XauB by other as yet undetermined genes. We did not find in XauB nor in XauC genes that lie between clusters F1 and F2; this region is XACSR9. In this region the gene XAC1927 has been shown to be important for citrus canker since mutations in it significantly decreased virulence  (Fig. 9A).LPS and O-antigen genesTwo of the XAC-specific regions contain genes related to lipopolysaccharides (LPS) that may play a role in the differing phenotypes between XAC and the two Xau strains. Region XACSR18 (XAC3596-3599) (which contains gene rfbC a truncated O-antigen biosynthesis protein) is immediately downstream of a CC-specific region (XAC3588-3595) (Fig. 10). An alignment with X. campestris pv. campestris ATCC33913 shows that these two regions together occupy roughly the same genomic locus as a 25-gene region related to LPS synthesis in X. campes-Figure 10 LPS and O-antigen synthesis (A) Model representing LPS and O-antigen synthesis in Xanthomonas Each of the numbers in this panel correspond to the same numbers in panel B. The model is based on the one described by Vorholter et al.  for X. campestris Genes for products 12, 14, 17 and 18 are absent fr om XAC, XauB and XauC. (B) Genomic context of the LPS synthesis-related genes in XAC and X. campestris The flanking regions (in pink) are conserved. The central region is not conserved between X. campestris on the one hand and XAC, XauB and XauC on the other. Moreover there is a group of genes that is shared and syntenic between the CC-species (in orange) and another (in green) that i s specific to XAC (XACSR18). (C) A XAC mutant for gene wzt showed a different virulence pattern . The wzt homolog in XauB and XauC does not have the C-terminal portion when compared to the XAC wzt gene. (D) Region XACSR1 contains several genes related to LPS biosynthesis. Numbers above gene arrows correspond to the annotations immediately below in this panel.
Moreira et al. BMC Genomics 2010, 11 :238 http://www.biomedcentra l.com/1471-2 164/11/238 Page 19 of 25tris pv. campestris  (Fig. 10). Immediately downstream of XACSR18 we find genes wzt (XAC3600) and wzm (XAC3601), also described as part of the LPS synthesis cluster in X. campestris . Laia et al.  obtained significant phenotypic differences in XAC (less necrosis, more water soaking and more hyperplasia, as compared to the wild type) by mutating gene XAC3600 (which codes for an ABC transporter ATP-binding protein). The wzt gene is present in both Xau genomes (XAUB_16610 and XAUC_09380), although both genes have lost the C-terminus when compared to the XAC wzt gene (Fig. 10). The wzm gene is also present in both Xau genomes (XAUB_ 16600 and XAUC_09370). XACSR1 (XAC0037-0063) contains several genes related to LPS synthesis, including two copies of asnB (XAC0051 and XAC0059), which codes for an asparagine synthase. A homolog of this gene in Pseudomonas aeruginosa has been implicated in O-antigen biosynthesis .The gene xacPNPGottig et al.  have identified a plant natriuretic peptide-like protein in XAC (xacPNP), encoded by gene XAC2654. They have shown that a XAC2654-deletion mutant resulted in more necrotic tissues and earlier bacterial cell death than in the wild type. XAC2654 lies between regions XACSR13 and XACSR14, and is flanked on both sides by phage-related genes. We have experimentally verified that neither XauB nor XauC appears to contain a homolog of XAC2654 (Additional file 5: Fig. S5). Therefore, xacPNP might be another gene contributing to the higher virulence displayed by XAC as compared to XauB and XauC.ConclusionsCitrus canker continues to be an economically important disease. The publication of the XAC strain 306 genome in 2002 opened up new avenues of research, and several important insights into the genetics of canker have been obtained since then, most of them cited here. Yet, understanding the genomic basis of a bacterial plant disease is a complex undertaking. By obtaining the genome sequences of two additional ci trus canker strains we have uncovered several new clues towards a thorough understanding of this disease. We have approached the problem from two basic perspectives. The first was to determine commonalities among the three CC strains that were not found in other Xanthomonas genomes. Such traits are excellent candidates for the general genomic basis of canker and/or adaptation to citrus hosts. The second was to carefully compare the three CC genomes to one another, with special attention to genes that XAC has that the others (apparently) do not, as well as genes present in the Xau genomes but absent in XAC. Because of the draft nature of the Xau genomes here presented, all results concerning gene absence in their sequences are tentative. However, our hybridization platform provided additional evidence for the specificity of XAC genes and regions with respect to the other two strains. Our most important findings are related to presence/ absence of effector genes. In addition to the already known pthA gene, the genes xopE3 and xopAI deserve special attention in future studies. Moreover, we have identified several genes (such as xacPNP ) that differentiate XAC from XauB and/or XauC. These genes or their homologs in other bacterial plant pathogens have demonstrated roles in virulence and/or host specificity. Hypotheses on their role in citrus canker and in host range differences between CC strains can now be tested experimentally. We anticipate that knowledge in regard to CC-specific effectors and other CC-specific genes will be used in the future to engineer citrus species with durable resistance to citrus canker, thus reducing the economic impact of this disease on the citrus industry worldwide. Such knowledge will also be crucial for dealing with new canker variants that may emerge in the field, exemplified by the recent detection of what appears to be a new variant of Xanthomonas fuscans subsp. aurantifolii in swingle citrumelo .MethodsBacterial strains and DNA sequencingThe Xanthomonas fuscans subsp. aurantifolii type B genome sequenced was strain 11122 (B-69), isolated from a Citrus limon tree in Argentina. The Xanthomonas fuscans subsp. aurantifolii type C genome sequenced was strain 10535 (IBSBF338), isolated from a Mexican lime tree in So Paulo state in Brazil. We sequenced the genomes using the Sanger technique as described previously , with sequencers ABI 3700 and ABI 3100. For XauB we generated both shotgun and cosmid libraries; for XauC only shotgun libraries were created. For XauB we obtained 114,874 reads; for XauC we obtained 114,805 reads. We estimate this provided about 15 average coverage for each genome.Microscopy and motility testsInitial attempts to visualize the flagellum of the two Xau strains with a scanning electron microscope by the normal routine (adhesion of the cultured bacteria on a cover slip with the help of the cationic compound poly-L-lysine, fixation, dehydration and critical point drying) failed because apparently attachment of the flagellum is very weak and it tends to fall off easily. To circumvent the problem, a diluted suspension of the bacterial culture was directly transferred onto a cover slip and excess of liquid eliminated. A moist chamber was made in a Petri dish
Moreira et al. BMC Genomics 2010, 11 :238 http://www.biomedcentra l.com/1471-2 164/11/238 Page 20 of 25where the cover slip was inserted together with a small plastic vial containing about 1 ml of 2% aqueous osmium tetroxide. The Petri dish was sealed and wrapped with aluminium foil and left overnight. The next morning the cover slip was removed, air dried and sputter-coated with gold, mounted on the stub and examined in a LEO 435 VP scanning electron microscope. Results are shown in Fig. 9B. For motility tests TSA and DYGS media were used . The Agar concentration of both media was changed (0.7%) so that the media were semisolid. For better growth visualization a phenol dye of 1% was added. Strains were placed in plates with solid media (TSA and DYGS) where isolated colonies were grown. The colonies were then placed in semisolid media. After 96 hours of incubation at 28C, bacterial growth was observed in the test tubes, and results are given in the table in Fig. 9B.Gum production determinationStrains were maintained both in autoclaved tap water at room temperature and at -80C in NA medium (3 g/l meat extract and 5 g/l peptone) containing 25% glycerol. The three strains were picked from -80C stock, streaked on solid TSA medium (10 g/l tryptone, 10 g/l sucrose, 1 g/l sodium glutamate, and 15 g/l agar) and grown overnight at 28C. One single isolated colony from each strain was streaked again on solid TSA medium and grown overnight at 28C. For XAC, a single colony was inoculated into 20 ml of liquid TSA medium (10 g/l tryptone, 10 g/l sucrose, 1 g/l sodium glutamate) in a 125 ml Erlenmeyer flask and incubated at 28C in a rotary shaker at 180 rpm for 17 hours (1.100 OD at 600 nm). A 125 ml Erlenmeyer flask containing 50 ml of liquid TSA medium was inoculated with 1 ml of XAC culture and incubated at 28C in a rotary shaker at 180 rpm for 5 hours (0.300 OD at 600 nm). This XAC culture was used as inoculum for xanthan gum production. For XauB and XauC, an inoculating loop was used to inoculate the bacteria from the solid TSA medium plates into 20 ml of liquid TSA medium in a 50 ml Falcon tube followed by incubation at 28C in a rotary shaker at 180 rpm for 24.5 hours (0.260 and 0.550 OD at 600 nm for XauB and XauC, respectively). One ml of XauB and XauC culture was inoculated into separate 125 ml Erlenmeyer flasks containing 50 ml of liquid TSA medium and incubated at 28C in a rotary shaker at 180 rpm for 15 hours (0.300 OD at 600 nm). The XauB and XauC cultures were used as inoculum for xanthan gum production. For gum production, 2.5 ml of each bacterial strain in liquid TSA medium (0.300 OD at 600 nm) was inoculated in three (triplicate) 250 ml Erlenmeyer flasks containing 100 ml of media for xanthan gum production (25 g/l glucose, 3 g/l yeast extract, 2 g/l K2HPO4, 0.1 g/l MgSO4.7H2O, pH 7.0 with 4 M HCl ) and incubated at 28C in a rotary shaker at 178 rpm for 96 h. The cells from the 96 h culture were centrifuged at 9.666 g for 40 min. The bacterial pellets were stored at -20C and the supernatants were transferred to 500 ml beakers. The gum was recovered from the supernatants by alcohol precipitation. Four grams of KCl were added to each beaker followed by agitation at room temperature for 15 min. Two volumes of cold isopropyl alcohol were added and the gum from each beaker was removed to pre-weighted plastic discs. After 16 h at 37C the discs were weighted again and the gum amount was calculated. The bacterial pellet from each culture was also obtained. For this, each pellet was transferred to a pre-weighted beaker and weighted again after 14 h at 70C.DNA MicroarrayFrom the shotgun libraries made for the sequencing of XAC strain 306 in 2001, 2,653 clones were selected for the design of a glass slide hybridization array. Inserts from selected clones were amplified by PCR with the M13-R and M13-F universal oligonucleotides. Products were purified and placed in duplicate slides, resulting in 6,144 probes, with 768 positive controls and 624 negative controls. Experimental validation was done by hybridization of total XAC DNA probes differentially stained, with 86% of all products with a hybridization signal greater than the average value for the noise signal plus two standard deviations. The final XACarray contains targets for the identification of 2,760 putative coding sequences (61% of all annotated protein-coding genes). The XACarray was used to compare XAC and XauC, and XAC and XauB, with XAC itself as a positive control. The results were similar (128 CDSs considered XAC-specific had similar ratios for both the XauB and XauC experiments, and 101 of these are in XACSRs); we report here detailed results only for the first experiment. Details of the array construction are described elsewhere (Moreira LM, Laia ML, de Souza RF, Zaini PA, da Silva ACR, Ferro JA, da Silva AM: Development and validation of a Xanthomonas citri subsp. citri DNA microarray platform ( XACarray ) generated from the shotgun libraries previously used in the sequencing of this bacterial genome, submitted).XacPNP verificationEvidence for the absence of genes coding for XacPNP homologs in XauB and XauC was obtained by PCR using gene-specific primers (Forward: GGACCAACAACGAATATC; Reverse: ATGGGAATAGTCATGAAAC). XAC was used as positive control.Assembly and genome annotationBase calling, genome assembly and visualization were done with the phred-phrap-consed package [99-101]. Contigs were trimmed to remove low phred quality regions at both ends. For both XauB and XauC all contigs
Moreira et al. BMC Genomics 2010, 11 :238 http://www.biomedcentra l.com/1471-2 164/11/238 Page 21 of 25larger than 1 kbp have average phred quality greater than 20, and nearly all (99% for XauB and 96% for XauC) these contigs have average phred quality greater than or equal to 40 (i.e. accuracy equal to or better than 1 error in 10,000 bp). To improve assembly many transposase sequences were masked, and reads were re-assembled to obtain the final result. We estimated the fraction obtained of the complete genome of XauB by dividing the total length of contigs by the total genome size of XAC, and by adding 1% because of the removed transposases. We did the same with XauC. To estimate the average gap size we divided the number of contigs by the estimated size of un-sequenced genome. Using paired-end reads scaffolds for the XauB and XauC chromosomes were obtained. The overall correctness of each scaffold was validated using a sliding-window GC-skew computation. Each scaffold was aligned against the XAC chromosome using MUMmer . Results are presented in additional file 1 (Fig. S1). The XauB and XauC genomes were automatically annotated with the Genome Reverse Compiler , with a few manual refinements. Ortholog groups were built using OrthoMCL .XAC-specific regionsFor determination of XAC-specific regions we relied on published data about pu tative genomic islands [14,31,105] and on AlienHun ter results . Islands whose genes were not found in the XauB and XauC genomes by BLAST  analysis and with differential hybridization signals in the XACarray were considered XAC-specific regions.Phylogeny reconstruction Species treeWe used a supermatrix approach as in previous work . Protein sequences of eleven Xanthomonas genomes (ingroup) and four Xylella genomes (outgroup) were clustered in 6,375 families using OrthoMCL . We then selected families with one and only one representative from each of the ingroup genomes and at least one outgroup protein, resulting in 1,666 families. Their sequences were aligned using MUSCLE  and the resulting alignments were concatenated. Non-informative columns were removed using Gblocks , resulting in 596,246 positions. RAxML  with the PROTGAMMAWAGF model was used to build the final tree.Pth treeThe same methodology as above was used, with the following differences. Representative pth nucleotide sequences from Xanthomonas species were retrieved from GenBank, and added to the set of XauB and XauC pth nucleotide sequences. A pth gene from Ralstonia solanacearum [GenBank: CAD1557.1 ] was used as an outgroup in a preliminary round of tree construction to ascertain root position. The list of gene sequences used to build the tree is given in additional file 6 (Table S6). The Tandem Repeat Finder program  with parameters 2,7,7,80,10,50,500,1 was used to mask the internal repeats. The masked regions were removed and the resulting sequences were aligned with MUSCLE. A few manual adjustments to the multiple alignment were made before running Gblocks, which yielded 1,679 positions. RAxML with the GTRGAMMA model was used to build the final tree. The bootstrap values obtained for the clade indicated by an asterisk in Fig. 3 are given in additional file 7 (Fig. S7).Effector analysisThe candidate T3SS effectors in the XauB and XauC genomes were identified using tBLASTn  analysis and Pfam domain  searches. For tBLASTn analysis, all known plant and animal pathogen effectors were used as query with an e-value threshold 0.00001. Pfam domains were searched for possible domains found in known effectors in the predicted set of ORFs of draft genome sequences. Candidate effectors were classified according to the nomenclature and classification scheme for effectors in xanthomonads recently described by White et al. .Database submissionThe draft genome sequences of XauB and XauC are available at GenBank under accession numbers ACPX00000000 and ACPY00000000, respectively.Additional materialAbbreviations CC: citrus canker; CDS: coding sequence; ORF: Open Reading Frame; PCR: polymerase chain reaction; T3SS: Type III se cretion system; T4SS: Type IV secretion system; XAC: Xanthomonas citri subsp. citri ; XACSR: XAC-specific region; Xau: Xanthomonas fuscans subsp. aurantifolii ; XauB: Xanthomonas fuscans subsp. aurantifolii strain B; XauC: Xanthomonas fuscans subsp. aurantifolii strain C. Additional file 1 Figure S1: whole chromosome alignments of XauB and XauC scaffolds against XAC. Additional file 2 Table S2: genes shared by XAC, XauB, and XauC but that are not found in other fully sequenced Xanthomonas and Xylella genomes. Additional file 3 Table S3: List of XAC-specific regions. Additional file 4 Figure S4: gene diagrams showing the contents of all 25 XAC-specific regions. Additional file 5 Figure S5: PCR results for gene xacPNP Additional file 6 Table S6: names, organisms, locus tags, and accession numbers for the nucleotide sequences used to build the pth gene phylogeny (Fig. 3 ). Additional file 7 Figure S7: topology of clade in Fig. 3 with bootstrap values.
Moreira et al. BMC Genomics 2010, 11 :238 http://www.biomedcentra l.com/1471-2 164/11/238 Page 22 of 25Authors' contributions JAF and ACRS conceived the project and oversaw genomic sequencing. RPL provided the strains. JCS oversaw all bi oinformatics work, including assembly, annotation and other specialized anal yses. MYN, EHO, and LAD helped with genome assembly. LAD and EHO did comp utational analyses. RIDT maintained the computational infrastructure. LMM conceived the XACarray. LMM and ALJF carried out hybridization experiments. HAP and LASN did the motility tests. EWK did microscopy analysis. DFRJ did the xanthan gum production experiment. JAF oversaw final PCR experiments. APF ran PCR analyses. NFA helped with automated genome annotation and manual refinement, created ortholog families, and generated the GenBank files. NFA, SSA and JCS did the phylogenetic analyses. APF, AMS, HAP, FEM, JCB, JCFO, LMM, LASN, LRF, MLL, MITF, RFS, and RIDT helped with data analyses. DJN and BJS helped create an effector database. NP did the effector analysis. NP, JBJ, BAV, LMM, and JCS interpreted effector analysis. LMM created all figures except Figs. 2 and 3. MLL, LRF, and JRN provided preliminary sections of the manuscript. LMM, BAV, and JCS wrote the final manuscript. All authors approved the final manuscript. Acknowledgements This work was supported by Fundao de Amparo Pesquisa do Estado de So Paulo (FAPESP), by Conselho Nacional de Desenvolvimento Cientfico e Tecnolgico (CNPq), by Coordenao para Ap erfeioamento de Pessoal de Ensino Superior (CAPES), and by FUNDECITRUS, all of them Brazilian funding agencies. We thank the following students for their he lp: Aloisio J. Almeida Jr., Maria L. R. Carrera, Gustavo G. L. Costa, Daniel M. de Faria, Daniel P. Gasparotto, Rafael A. Homem, Vanessa C. Morgan, Julia M. Perdigueiro, Marcelo C. Perez, and Paula F. S. Tabatini. Author Details1Departamento de Cincias Biolgicas, Instituto de Cincias Exatas e Biolgicas, Campus Morro do Cruzeiro, Universidade Federal de Ouro Preto, Ouro Preto, MG, Brazil, 2Departamento de Bioqumic a, Instituto de Qumica, Universidade de So Paulo, So Paulo, SP, Brazil, 3Department of Plant Pathology, University of Florida, Gainesville, FL, USA, 4Department of Plant Pathology, Physiology and Weed Scienc es, Virginia Polytechnic Institute and State University, Blacksburg, VA, USA, 5Departamento de Tecnologia, Faculdade de Cincias Agrrias e Veterinrias de Jaboticabal, UNESP Univ. Estadual Paulista, Jaboticabal, SP, Brazil, 6Departamento de Cincias Biolgicas, Campus de Diadema, Universidade Federal de So Paulo, So Paulo, SP, Brazil, 7Departamento de Melhoramento e Nutr io Animal, Faculdade de Medicina Veterinria e Zootecnia de Botucatu, UNESP Univ. Estadual Paulista, SP, Brazil, 8Departamento de Engenharia Florestal, Centro de Cincias Agroveterinrias, Universidade do Estado de Sant a Catarina, Lages, SC, Brazil, 9Allelyx Applied Genomics, Campinas, SP, Brazil, 10Ncleo de apoio pesquisa em microscopia eletrnica aplicada pesquisa agropecu ria, Escola Superior de Agricultura Luiz de Queiroz, Universidade de So Paulo, Piracicaba, SP, Brazil, 11Laboratrio de Bacteriologia Vegetal, Instituto Biolgico Campinas, Campinas, SP, Brazil, 12Escola de Artes, Cincias, e Humanida des, Universidade de So Paulo, So Paulo, SP, Brazil, 13Faculdade de Computao, Universidade Federal do Mato Grosso do Sul, Campo Grande, MS, Brazil, 14Laboratrio de Bioinformtica, Instituto de Computao, Universidade Estadual de Campinas, Campinas, SP, Brazil, 15Virginia Bioinformatics Institute, Virg inia Polytechnic Institute and State University, Blacksburg, VA, USA, 16Department of Computer Science, Virginia Polytechnic Institute and State University, Blacksburg, VA, USA, 17Instituto Agronmico do Paran, Londrina, PR, Brazil, 18Institute of Food and Agricultural Sciences, Mid-Florida Research & Edu cation Center, University of Florida, Gainesville, FL, USA and 19Department of Plant & Microbial Biology, University of California, Berkeley, Berkeley, CA, USA References 1.Timmer LW, Garnsey SM, Graham JH, (eds): Compendium of Citrus Diseases 2nd edition. St. Paul, MN: American Phytopathological Society; 2000. 2.Gottwald TR, Graham JH, Schubert TS: Citrus canker: The pathogen and its impact Plant Health Progress 2002 [http:// plantmanagementnetwork.org/pub/php/review/citruscanker ]. 3.Brunings AM, Gabriel DW: Xanthomonas citri : breaking the surface Molecular Plant Pathology 2003, 4(3): 141-157. 4.Schubert TS, Miller JW: Bacterial citrus canker In Plant Pathology Circular Volume 377 Florida Department of Agri culture and Consumer Services DoPI, Gainesville, FL; 1996:110-111. 5.Koizumi M: Citrus canker: the world situation University of Florida: Lake Alfred; 1985. 6.Graham JH, Gottwald TR, Cubero J, Achor DS: Xanthomonas axonopodis pv. citri : factors affecting successful eradication of citrus canker Molecular Plant Pathology 2004, 5(1): 1-15. 7.Bitancourt AA: O Cancro Ctrico O Biolgico 1957, 23: 110-111. 8.Civerolo EL: Bacterial canker disease of citrus Journal of the Rio Grande Valley Horticultural Society 1984, 37: 127-145. 9.Malavolta Jnior VA, Yamashiro T, Nogueira EMC, Feichtenberger E: Distribuio do tipo C de Xanthomonas campestris pv. citri no Estado de So Paulo Summa Phytopathologica 1984, 10: 11. 10.Egel DS, Graham JH, Stall RE: Genomic Relatedness of Xanthomonas campestris Strains Causing Diseases of Citrus Applied and Environmental Microbiology 1991, 57(9): 2724-2730. 11.Gabriel D, Hunter G, Kingsley M, Miller J, Lazo G: Clonal Population Structure of Xanthomonas campestris and Genetic Diversity Among Citrus Canker Strains Molecular Plant-Microbe Interactions 1988, 1(2): 59-65. 12.da Silva ACR, Ferro JA, Reinach FC, Farah CS, Furlan LR, Quaggio RB, Monteiro-Vitorello CB, Van Sluys MA, Almeida NF, Alves LMC, do Amaral AM, Bertolini MC, Camargo LE, Cama rotte G, Cannavan F, Cardozo J, Chambergo F, Ciapina LP, Cicarelli RM, Coutinho LL, Cursino-Santos JR, ElDorry H, Faria JB, Ferreira AJ, Ferreira RC, Ferro MI, Formighieri EF, Franco MC, Greggio CC, Gruber A, et al.: Comparison of the genomes of two Xanthomonas pathogens with differing host specificities Nature 2002, 417(6887): 459-463. 13.Moreira LM, de Souza RE, Almeida NF, Setubal JC, Oliveira JCF, Furlan LR, Ferro JA, da Silva ACR: Comparative genomics analyses of citrusassociated bacteria Annual Review of Phytopathology 2004, 42: 163-184. 14.Moreira LM, De Souza RF, Digiampietri LA, Da Silva ACR, Setubal JC: Comparative analyses of Xanthomonas and Xylella complete genomes Omics 2005, 9(1): 43-76. 15.Swarup S, Feyter RD, Brlansky RH, Gabriel DW: A pathogenicity locus from Xanthomonas citri enables strains from several pathovars of X. campestris to elicit cankerlike lesions on citrus Phytopathology 1991, 81: 802-809. 16.Duan YP, Castaneda AL, Zhao G, Erdos G, Gabriel DW: Expression of a single, host-specific, bacterial pathog enicity gene in plant cells elicits division, enlargement and cell death Molecular Plant-Microbe Interactions 1999, 12: 556-560. 17.Al-Saadi A, Reddy JD, Duan YP, Brunings AM, Yuan QP, Gabriel DW: All five host-range variants of Xanthomonas citri carry one pthA homolog with 17.5 repeats that determines pathogenicity on citrus, but none determine host-range variation Molecular Plant-Mi crobe Interactions 2007, 20(8): 934-943. 18.Swarup S, Yang Y, Kingsley MT, Gabriel DW: A Xanthomonas citri pathogenicity gene, pthA pleiotropically encodes gratuitous avirulence on nonhosts Mol Plant Microbe Interact 1992, 5(3): 204-213. 19.Verniere C, Hartung JS, Pruvost OP, Civerolo EL, Alvarez AM, Maestri P, Luisetti J: Characterization of phenotypically distinct strains of Xanthomonas axonopodis pv. citri from Southwest Asia European Journal of Plant Pathology 1998, 104(5): 477-487. 20.Sun XA, Stall RE, Jones JB, Cubero J, Gottwald TR, Graham JH, Dixon WN, Schubert TS, Chaloux PH, Stromberg VK, Lacy GH, Sutton BD: Detection and characterization of a new strain of citrus canker bacteria from key Mexican lime and Alemow in South Florida Plant Disease 2004, 88(11): 1179-1188. 21.Rybak M, Minsavage GV, Stall RE, Jones JB: Identification of Xanthomonas citri ssp citri host specificity genes in a heterologous expression host Molecular Plant Pathology 2009, 10(2): 249-262. 22.Ngoc LBT, Verniere C, Jouen E, Ah-You N, Lefeuvre P, Chiroleu F, Gagnevin L, Pruvost O: Amplified fragment length po lymorphism and multilocus sequence analysis-based genotypi c relatedness among pathogenic variants of Xanthomonas citri pv. citri and Xanthomonas campestris pv. bilvae International Journal of Systematic and Evolutionary Microbiology 2009 in press. 23.Namekata T: Estudos comparativos entre Xanthomonas citri (Hasse) Dow., agente causal do cancro ctrico e Xanthomonas citri (Hasse) Dow., Received: 16 December 2009 Accepted: 13 April 2010 Published: 13 April 2010This article is available from: http: //www.biomedcentral.com/1471-2164/11/238 2010 Moreira et al; licensee BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons. org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. BMC Genomics 2010, 11 :238
Moreira et al. BMC Genomics 2010, 11 :238 http://www.biomedcentra l.com/1471-2 164/11/238 Page 23 of 25n.f.sp. aurantifolia, agente causal da cancrose do limoeiro Galego Piracicaba: University of So Paulo; 1971. 24.Jaciani FJ, Destefano SA, Rodrig ues Neto J, Belasque Jr J: Detection of a New Bacterium Related to Xanthomonas fuscans subsp. aurantifolii Infecting Swingle Ci trumelo in Brazil Plant Disease 2009, 93(10): 1074. 25.Humphries J: Bacteriology London John Murray; 1974. 26.Schaad NW, Jones JB, Lacy GH: Gram-negative bacteria: Xanthomonas In Laboratory Guide for Identification of Plant Pathogenic Bacteria 3rd edition. Edited by: Schaad NW, Jone s JB, Chun W. St. Paul: APS Press; 2001:175-200. 27.Rodrigues Neto J, Mala volta VA Jr, Victor O: Meio simples para isolamento e cultivo de Xanthomonas campestris pv. citri Tipo B [abstract] Summa Phytopathologica 1986, 12(1-2): 16. 28.Cubero J, Graham JH: Genetic relationship among worldwide strains of Xanthomonas causing canker in citrus species and design of new primers for their iden tification by PCR Appl Environ Microbiol 2002, 68(3): 1257-1264. 29.Leite RP Jr, Egel DS, Stall RE: Genetic analysis of hrp -related DNA sequences of Xanthomonas campestris strains causing diseases of citrus Appl Environ Microbiol 1994, 60(4): 1078-1086. 30.Vauterin L, Yang P, Hoste B, Vancanne yt M, Civerolo EL, Swings J, Kersters K: Differentiation of Xanthomonas campestris pv. citri Strains by Sodium Dodecyl Sulfate-Polyacrylamide Gel-Ele ctrophoresis of Proteins, FattyAcid Analysis, and DNA-DNA Hybridization International Journal of Systematic Bacteriology 1991, 41(4): 535-542. 31.Laia ML, Moreira LM, Dezajacomo J, Brigati JB, Ferreira CB, Ferro MI, Silva AC, Ferro JA, Oliveira JC: New genes of Xanthomonas citri subsp. citri involved in pathogenesis and adaptation revealed by a transposonbased mutant library BMC Microbiol 2009, 9: 12. 32.Schaad NW, Postnikova E, Lacy GH, Sechler A, Agarkova I, Stromberg PE, Stromberg VK, Vidaver AK: Reclassification of Xanthomonas campestris pv. citri (ex Hasse 1915) Dye 1978 forms A, B/C/D, and E as X. smithii subsp. citri (ex Hasse) sp. nov. nom. rev. comb. nov., X. fuscans subsp. aurantifolii (ex Gabriel 1989) sp. nov. nom. rev. comb. nov., and X. alfalfae subsp. citrumelo (ex Riker and Jones) Gabriel et al., 1989 sp. nov. nom. rev. comb. nov.; X. campestris pv malvacearum (ex smith 1901) Dye 1978 as X. smithii subsp. smithii nov. comb. nov. nom. nov.; X. campestris pv. alfalfae (ex Riker and Jones, 1935) dye 1978 as X. alfalfae subsp. alfalfae (ex Riker et al., 1935) sp. nov. nom. rev.; and "var. fuscans" of X. campestris pv. phaseoli (ex Smith, 1987) Dye 1978 as X. fuscans subsp. fuscans sp. nov Syst Appl Microbiol 2005, 28(6): 494-518. 33.Schaad NW, Postnikova E, Lacy G, Sechler A, Agarkova I, Stromberg PE, Stromberg VK, Vidaver AK: Emended classification of xanthomonad pathogens on citrus Syst Appl Microbiol 2006, 29(8): 690-695. 34.Young JM, Park DC, Shearman HM, Fargier E: A multilocus sequence analysis of the genus Xanthomonas Syst Appl Microbiol 2008, 31(5): 366-377. 35.Almeida NF, Yan S, Cai R, Clarke CR, Morris CE, Schaad NW, Lacy GH, Jones JB, Castillo JA, Bull CT, Leman S, Guttman DS, Setubal JC, Vinatzer BA: PAMDB, A Multilocus Sequence Typing & Analysis Database and Website for Plant-Associated and Plant-Pathogenic Microorganisms Phytopathology 2010, 100(3): 208-215. 36.Usher KC, Ozkan E, Gardner KH, Deisenhofer J: The plug domain of FepA, a TonB-dependent transport protein from Escherichia coli binds its siderophore in the absence of the transmembrane barrel domain. Proc Natl Acad Sci USA 2001, 98(19): 10676-10681. 37.Uria-Nickelsen MR, Lea dbetter ER, Godchaux W: Comparative aspects of utilization of sulfonate and other sulfur sources by Escherichia coli K12 Arch Microbiol 1994, 161(5): 434-438. 38.Kim KD, Ahn JH, Kim T, Park SC, Seong CN, Song HG, Ka JO: Genetic and phenotypic diversity of fenitrothion-degrading bacteria isolated from soils J Microbiol Biotechnol 2009, 19(2): 113-120. 39.Yang J, Yang C, Jiang H, Qiao C: Overexpression of methyl parathion hydrolase and its application in de toxification of organophosphates Biodegradation 2008, 19(6): 831-839. 40.Vernikos GS, Parkhill J: Interpolated variable order motifs for identification of horizontally acquired DNA: revisiting the Salmonella pathogenicity islands Bioinformatics 2006, 22(18): 2196-2203. 41.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. 42.Hajri A, Brin C, Hunault G, Lardeux F, Lemaire C, Manceau C, Boureau T, Poussier S: A "repertoire for repertoire" hypothesis: repertoires of type three effectors are candidate determinants of host specificity in Xanthomonas PLoS One 2009, 4(8): e6632. 43.Heuer H, Yin YN, Xue QY, Smalla K, Guo JH: Repeat domain diversity of avrBs3-like genes in Ralstonia solanacearum strains and association with host preferences in the field Appl Environ Microbiol 2007, 73(13): 4379-4384. 44.El Yacoubi B, Brunings AM, Yuan Q, Shankar S, Gabriel DW: In planta horizontal transfer of a major pathogenicity effector gene Applied and Environmental Microbiology 2007, 73(5): 1612-1621. 45.Altschul SF, Madden TL, Schaffer AA, Zh ang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs Nucleic Acids Research 1997, 25(17): 3389-3402. 46.Nimchuk ZL, Fisher EJ, Desveaux D, Chang JH, Dangl JL: The HopX (AvrPphE) family of Pseudomonas syringae type III effectors require a catalytic triad and a novel N-terminal domain for function Mol Plant Microbe Interact 2007, 20(4): 346-357. 47.Schechter LM, Roberts KA, Jamir Y, Alfano JR, Collmer A: Pseudomonas syringae type III secretion system targeting signals and novel effectors studied with a Cya translocation reporter J Bacteriol 2004, 186(2): 543-555. 48.Stevens C, Bennett MA, Athanassopoulo s E, Tsiamis G, Ta ylor JD, Mansfield JW: Sequence variations in alleles of the avirulence gene avrPphE. R2 from Pseudomonas syringae pv. phaseolicola lead to loss of recognition of the AvrPphE protein within bean cells and a gain in cultivar-specific virulence Mol Microbiol 1998, 29(1): 165-177. 49.Stavrinides J, Ma W, Guttman DS: Terminal reassortment drives the quantum evolution of type III effectors in bacterial pathogens PLoS Pathog 2006, 2(10): e104. 50.Oh HS, Kvitko BH, Morello JE, Collmer A: Pseudomonas syringae lytic transglycosylases coregulated with the type III secretion system contribute to the translocation of effector proteins into plant cells J Bacteriol 2007, 189(22): 8277-8289. 51.Kim JG, Li XY, Roden JA, Taylor KW, Aakre CD, Su B, Lalonde S, Kirik A, Chen YH, Baranage G, McLane H, Martin GB, Mudgett MB: Xanthomonas T3S Effector XopN Suppresses PAMP-Triggered Immunity and Interacts with a Tomato Atypical Receptor-Like Kinase and TFT1 Plant Cell 2009, 21(4): 1305-1323. 52.Noel L, Thieme F, Nennstiel D, Bonas U: Two novel type III-secreted proteins of Xanthomonas campestris pv. vesicatoria are encoded within the hrp pathogenicity island J Bacteriol 2002, 184(5): 1340-1348. 53.Sugio A, Yang B, White FF: Characterization of the hrpF pathogenicity peninsula of Xanthomonas oryzae pv. oryzae Mol Plant Microbe Interact 2005, 18(6): 546-554. 54.Liakopoulou-Kyriakides M, Psomas SK, Kyriakidis DA: Xanthan gum production by Xanthomonas campestris w.t. fermentation from chestnut extract Applied Biochemistry and Biotechnology 1999, 82(3): 175-183. 55.Katzen F, Ferreiro DU, Oddo CG, Ielm ini MV, Becker A, Puhler A, Ielpi L: Xanthomonas campestris pv. campestris gum mutants: Effects on xanthan biosynthesis and plant virulence Journal of Bacteriology 1998, 180(7): 1607-1617. 56.Boher B, Nicole M, Potin M, Geiger JP: Extracellular polysaccharides from Xanthomonas axonopodis pv. manihotis interact with cassava cell walls during pathogenesis Molecular Plant-Microbe Interactions 1997, 10(7): 803-811. 57.Vojnov AA, Slater H, Daniels MJ, Dow JM: Expression of the gum operon directing xanthan biosynthesis in Xanthomonas campestris and its regulation in planta Molecular Plant-Microbe Interactions 2001, 14(6): 768-774. 58.Tang JL, Liu YN, Barber CE, Dow JM, Wootton JC, Daniels MJ: Genetic and Molecular Analysis of a Cluster of Rpf Genes Involved in Positive Regulation of Synthesis of Extracellular Enzymes and Polysaccharide in Xanthomonas campestris pathovar campestris Molecular & General Genetics 1991, 226(3): 409-417. 59.Newman KL, Almeida RP, Purcell AH, Lindow SE: Cell-cell signaling controls Xylella fastidiosa interactions with both insects and plants Proc Natl Acad Sci USA 2004, 101(6): 1737-1742. 60.He YW, Zhang LH: Quorum sensing and virulence regulation in Xanthomonas campestris FEMS Microbiol Rev 2008, 32(5): 842-857.
Moreira et al. BMC Genomics 2010, 11 :238 http://www.biomedcentra l.com/1471-2 164/11/238 Page 24 of 2561.Andre A, Maucourt M, Moing A, Rolin D, Renaudin J: Sugar import and phytopathogenicity of Spiroplasma citri : Glucose and fructose play distinct roles Molecular Plant-Microbe Interactions 2005, 18(1): 33-42. 62.Gaurivaud P, Danet JL, Laigret F, Garnier M, Bove JM: Fructose utilization and phytopathogenicity of Spiroplasma citri Molecular Plant-Microbe Interactions 2000, 13(10): 1145-1155. 63.Gaurivaud P, Laigret F, Garnier M, Bove JM: Fructose utilization and pathogenicity of Spiroplasma citri : characterization of the fructose operon Gene 2000, 252(1-2): 61-69. 64.Gaurivaud P, Laigret F, Verdin E, Garnier M, Bove JM: Fructose operon mutants of Spiroplasma citri Microbiology-Sgm 2000, 146: 2229-2236. 65.Dow JM, Feng JX, Barber CE, Tang JL, Daniels MJ: Novel genes involved in the regulation of pathogenicity factor production within the rpf gene cluster of Xanthomonas campestris Microbiology-Uk 2000, 146: 885-891. 66.Tang JL, Gough CL, Daniels MJ: Cloning of Genes Involved in Negative Regulation of Production of Extracellular Enzymes and Polysaccharide of Xanthomonas campestris pathovar campestris Molecular & General Genetics 1990, 222(1): 157-160. 67.Lemos EGD, Alves LMC, Campanharo JC: Genomics-based design of defined growth media for the plant pathogen Xylella fastidiosa Fems Microbiology Letters 2003, 219(1): 39-45. 68.da Silva FR, Vettore AL, Kemper EL, Leite A, Arruda P: Fastidian gum: the Xylella fastidiosa exopolysaccharide possibly involved in bacterial pathogenicity Fems Microbiology Letters 2001, 203(2): 165-171. 69.Gottig N, Garavaglia BS, Garofa lo CG, Orellano EG, Ottado J: A Filamentous Hemagglutinin-Like Protein of Xanthomonas axonopodis pv. citri the Phytopathogen Responsible for Citrus Canker, Is Involved in Bacterial Virulence PLoS ONE 2009, 4(2): e4358. 70.Pendrak ML, Perry RD: Characterization of a hemin-storage locus of Yersinia pestis Biol Met 1991, 4(1): 41-47. 71.Hinnebusch BJ, Perry RD, Schwan TG: Role of the Yersinia pestis hemin storage (hms) locus in the transmission of plague by fleas Science 1996, 273(5273): 367-370. 72.Lillard JW Jr, Bearden SW Fetherston JD, Perry RD: The haemin storage (Hms+) phenotype of Yersinia pestis is not essential for the pathogenesis of bubo nic plague in mammals Microbiology 1999, 145(1): 197-209. 73.Jarrett CO, Deak E, Isherwood KE Oyston PC, Fischer ER, Whitney AR, Kobayashi SD, DeLeo FR, Hinnebusch BJ: Transmission of Yersinia pestis from an infectious biofilm in the flea vector Journal of Infectious Diseases 2004, 190(4): 783-792. 74.Guilhabert MR, Kirkpatrick BC: Identification of Xylella fastidiosa antivirulence genes: hemagglutinin adhesins contribute a biofilm maturation to X. fastidiosa and colonization and attenuate virulence Mol Plant Microbe Interact 2005, 18(8): 856-868. 75.Soto SM, Smithson A, Martinez JA, Horcajada JP, Mensa J, Vila J: Biofilm formation in uropathogenic Escherichia coli strains: Relationship with prostatitis, urovirulence factors and antimicrobial resistance Journal of Urology 2007, 177(1): 365-368. 76.Arthur M, Johnson CE, Rubin RH, Arbeit RD, Campanelli C, Kim C, Steinbach S, Agarwal M, Wilkinson R, Goldstein R: Molecular Epidemiology of Adhesin and Hemolysin Virulence Factors among Uropathogenic Escherichia coli Infection and Immunity 1989, 57(2): 303-313. 77.Alegria MC, Souza DP, Andrade MO, Do cena C, Khater L, Ramos CH, da Silva AC, Farah CS: Identification of new protein-protein interactions involving the products of the chromosomeand plasmid-encoded type IV secretion loci of the phytopathogen Xanthomonas axonopodis pv. citri J Bacteriol 2005, 187(7): 2315-2325. 78.Barrios AF, Zuo R, Ren D, Wood TK: Hha, YbaJ, and OmpA regulate Escherichia coli K12 biofilm formation and conjugation plasmids abolish motility Biotechnol Bioeng 2006, 93(1): 188-200. 79.Das A, Rangaraj N, Sonti RV: Multiple adhesin-like functions of Xanthomonas oryzae pv. oryzae are involved in promoting leaf attachment, entry, and virulence on rice Mol Plant Microbe Interact 2009, 22(1): 73-85. 80.Darsonval A, Darrasse A, Durand K, Bureau C, Cesbron S, Jacques MA: Adhesion and fitness in the bean phyllosphere and transmission to seed of Xanthomonas fuscans subsp. fuscans Mol Plant Microbe Interact 2009, 22(6): 747-757. 81.Craig L, Li J: Type IV pili: paradoxes in form and function Curr Opin Struct Biol 2008, 18(2): 267-277. 82.Meng Y, Li Y, Galvani CD, Hao G, Turner JN, Burr TJ, Hoch HC: Upstream migration of Xylella fastidiosa via pilus-driven twitching motility J Bacteriol 2005, 187(16): 5560-5567. 83.De La Fuente L, Montanes E, Meng Y, Li Y, Burr TJ, Hoch HC, Wu M: Assessing adhesion forces of type I and type IV pili of Xylella fastidiosa bacteria by use of a microfluidic flow chamber Appl Environ Microbiol 2007, 73(8): 2690-2696. 84.Li Y, Hao G, Galvani CD, Meng Y, De La Fuente L, Hoch HC, Burr TJ: Type I and type IV pili of Xylella fastidiosa affect twitching motility, biofilm formation and cell-cell aggregation Microbiology 2007, 153(Pt 3): 719-726. 85.De La Fuente L, Burr TJ, Hoch HC: Mutations in type I and type IV pilus biosynthetic genes affect twitching motility rates in Xylella fastidiosa J Bacteriol 2007, 189(20): 7507-7510. 86.De La Fuente L, Burr TJ, Hoch HC: Autoaggregation of Xylella fastidiosa cells is influenced by type I and type IV pili Appl Environ Microbiol 2008, 74(17): 5579-5582. 87.Sidhu VK, Vorholter FJ, Niehaus K, Watt SA: Analysis of outer membrane vesicle associated proteins isolated from the plant pathogenic bacterium Xanthomonas campestris pv. campestris BMC Microbiol 2008, 8: 87. 88.Qian W, Jia Y, Ren SX, He YQ, Feng JX, Lu LF, Sun Q, Ying G, Tang DJ, Tang H, Wu W, Hao P, Wang L, Jiang BL, Zeng S, Gu WY, Lu G, Rong L, Tian Y, Yao Z, Fu G, Chen B, Fang R, Qiang B, Chen Z, Zhao GP, Tang JL, He C: Comparative and functional genomic analyses of the pathogenicity of phytopathogen Xanthomonas campestris pv. campestris Genome Res 2005, 15(6): 757-767. 89.Cascales E, Christie PJ: The versatile bacterial type IV secretion systems Nature Reviews Microbiology 2003, 1(2): 137-149. 90.Backert S, Meyer TF: Type IV secretion systems and their effectors in bacterial pathogenesis Current Opinion in Microbiology 2006, 9(2): 207-217. 91.Salanoubat M, Genin S, Artiguenav e F, Gouzy J, Mangenot S, Arlat M, Billault A, Brottier P, Camus JC, Cattolico L, Ch andler M, Choisne N, Claudel-Renard C, Cunnac S, Demange N, Gaspin C, Lavi e M, Moisan A, Robert C, Saurin W, Schiex T, Siguier P, Thbault P, Whalen M, Wincker P, Levy M, Weissenbach J, Boucher CA: Genome sequence of the plant pathogen Ralstonia solanacearum Nature 2002, 415(6871): 497-502. 92.Wood DW, Setubal JC, Kaul R, Monks DE, Kitajima JP, Okura VK, Zhou Y, Chen L, Wood GE, Almeida NF Jr, Woo L, Chen Y, Paulsen IT, Eisen JA, Karp PD, Bovee D Sr, Chapman P, Clendenning J, Deatherage G, Gillet W, Grant C, Kutyavin T, Levy R, Li MJ, McClella nd E, Palmieri A, Raymond C, Rouse G, Saenphimmachak C, Wu Z, et al. : The genome of the natural genetic engineer Agrobacterium tumefaciens C58 Science 2001, 294(5550): 2317-2323. 93.Simpson AJG, Reinach FC, Arruda P, Abreu FA, Acencio M, Alvarenga R, Alves LMC, Araya JE, Baia GS, Baptista CS, Barros MH, Bonaccorsi ED, Bordin S, Bov JM, Briones MR, Bueno MR, Camargo AA, Camargo LE, Carraro DM, Carrer H, Colauto NB, Colombo C, Costa FF, Costa MC, Costa-Neto CM, Coutinho LL, Cristofani M, Dias-N eto E, Docena C, El-Dorry H, et al. : The genome sequence of the plant pathogen Xylella fastidiosa Nature 2000, 406(6792): 151-157. 94.Schneiker S, Keller M, Droge M, Lanka E, Puhler A, Selbitschka W: The genetic organization and evolution of the broad host range mercury resistance plasmid pSB102 isolated from a microbial population residing in the rhizosphere of alfalfa Nucleic Acids Research 2001, 29(24): 5169-5181. 95.Tauch A, Schneiker S, Selbitschka W, Puhler A, van Overbeek LS, Smalla K, Thomas CM, Bailey MJ, Forney LJ, Weightman Ceglowski P, Pembroke T, Tietze E, Schrder G, Lanka E, van Elsas JD: The complete nucleotide sequence and environmental distribution of the cryptic, conjugative, broad-host-range plasmid pIPO2 isolated from bacteria of the wheat rhizosphere Microbiology 2002, 148(Pt 6): 1637-1653. 96.Vorholter FJ, Niehaus K, Puhler A: Lipopolysaccharide biosynthesis in Xanthomonas campestris pv. campestris : a cluster of 15 genes is involved in the biosynthesis of the LPS O-antigen and the LPS core Mol Genet Genomics 2001, 266(1): 79-95. 97.Rocchetta HL, Burrows LL, Lam JS: Genetics of O-antigen biosynthesis in Pseudomonas aeruginosa Microbiol Mol Biol Rev 1999, 63(3): 523-553. 98.Gottig N, Garavaglia BS, Daurelio LD, Valentine A, Gehring C, Orellano EG, Ottado J: Xanthomonas axonopodis pv. citri uses a plant natriuretic
Moreira et al. BMC Genomics 2010, 11 :238 http://www.biomedcentra l.com/1471-2 164/11/238 Page 25 of 25peptide-like protein to modify host homeostasis Proc Natl Acad Sci USA 2008, 105(47): 18631-18636. 99.Ewing B, Green P: Base-calling of automated sequencer traces using phred. II. Error probabilities Genome Res 1998, 8(3): 186-194. 100.Ewing B, Hillier L, Wendl MC, Green P: Base-calling of automated sequencer traces using phre d. I. Accuracy assessment Genome Res 1998, 8(3): 175-185. 101.Gordon D, Abajian C, Green P: Consed: a graphical tool for sequence finishing Genome Res 1998, 8(3): 195-202. 102.Kurtz S, Phillippy A, Delcher AL Smoot M, Shumway M, Antonescu C, Salzberg SL: Versatile and open software for comparing large genomes Genome Biol 2004, 5(2): R12. 103.Warren AS, Setubal JC: The Genome Reverse Compiler: an explorative annotation tool BMC Bioinformatics 2009, 10: 35. 104.Li L, Stoeckert CJ, Roos DS Jr: OrthoMCL: identification of ortholog groups for eukaryotic genomes Genome Res 2003, 13(9): 2178-2189. 105.Lima WC, Van Sluys MA, Menck CF: Non-gamma-proteobacteria gene islands contribute to the Xanthomonas genome Omics-a Journal of Integrative Biology 2005, 9(2): 160-172. 106.Setubal JC, dos Santos P, Goldman BS Ertesvag H, Espin G, Rubio LM, Valla S, Almeida NF, Balasubramanian D, Crom es L, Curatti L, Du Z, Godsy E, Goodner B, Hellner-Burris K, Hernande z JA, Houmiel K, Imperial J, Kennedy C, Larson TJ, Latreille P, Ligon LS, Lu J, Maerk M, Miller NM, Norton S, O'Carroll IP, Paulsen I, Raulfs EC, Roemer R, et al. : Genome sequence of Azotobacter vinelandii an obligate aerobe specialized to support diverse anaerobic metabolic processes J Bacteriol 2009, 191(14): 4534-4545. 107.Edgar RC: MUSCLE: multiple sequence alignment with high accuracy and high throughput Nucleic Acids Res 2004, 32(5): 1792-1797. 108.Castresana J: Selection of conserved blocks from multiple alignments for their use in phylogenetic analysis Mol Biol Evol 2000, 17(4): 540-552. 109.Stamatakis A: RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models Bioinformatics 2006, 22(21): 2688-2690. 110.Benson G: Tandem repeats finder: a program to analyze DNA sequences Nucleic Acids Res 1999, 27(2): 573-580. 111.Finn RD, Tate J, Mistry J, Coggill PC Sammut SJ, Hotz HR, Ceric G, Forslund K, Eddy SR, Sonnhammer EL, Bateman A: The Pfam protein families database Nucleic Acids Research 2008:D281-288. 112.White FF, Potnis N, Jones JB, Koebnik R: The type III effectors of Xanthomonas Molecular Plant Pathology 2009, 10(6): 749-766. 113.Goodner B, Hinkle G, Gattung S, Mill er N, Blanchard M, Qurollo B, Goldman BS, Cao Y, Askenazi M, Halling C, Mullin L, Houmiel K, Gordon J, Vaudin M, Iartchouk O, Epp A, Liu F, Wollam C, Allinger M, Doughty D, Scott C, Lappas C, Markelz B, Flanagan C, Crow ell C, Gurson J, Lomo C, Sear C, Strub G, Cielo C, et al. : Genome sequence of the plant pathogen and biotechnology agent Agrobacterium tumefaciens C58 Science 2001, 294(5550): 2323-2328. 114.Segal G, Russo JJ, Shuman HA: Relationships between a new type IV secretion system and the icm/dot virulence system of Legionella pneumophila Mol Microbiol 1999, 34(4): 799-809. 115.Cazalet C, Rusniok C, Bruggemann H, Zidane N, Magnier A, Ma L, Tichit M, Jarraud S, Bouchier C, Vandenesch F, Kunst F, Etienne J, Glaser P, Buchrieser C: Evidence in the Legionella pneumophila genome for exploitation of host cell functions and high genome plasticity Nat Genet 2004, 36(11): 1165-1173. 116.Bolland S, Llosa M, Avila P, de la Cruz F: General organization of the conjugal transfer genes of the IncW plasmid R388 and interactions between R388 and IncN and IncP plasmids J Bacteriol 1990, 172(10): 5795-5802. 117.Marques MV, da Silva AM, Gomes SL: Genetic organization of plasmid pXF51 from the plant pathogen Xylella fastidiosa Plasmid 2001, 45(3): 184-199. 118.Lee BM, Park YJ, Park DS, Kang HW, Kim JG, Song ES, Park IC, Yoon UH, Hahn JH, Koo BS, Lee GB, Kim H, Park HS, Yoon KO, Kim JH, Jung CH, Koh NH, Seo JS, Go SJ: The genome sequence of Xanthomonas oryzae pathovar oryzae KACC10331, the bacterial blight pathogen of rice Nucleic Acids Res 2005, 33(2): 577-586. 119.Kearney B, Staskawicz BJ: Widespread Distribution and Fitness Contribution of Xanthomonas campestris Avirulence Gene Avrbs2 Nature 1990, 346(6282): 385-386. 120.Thieme F, Szczesny R, Urban A, Kirchner O, Hause G, Bonas U: New type III effectors from Xanthomonas campestris pv. vesicatoria trigger plant reactions dependent on a conserved N-myristoylation motif Molecular Plant-Microbe Interactions 2007, 20(10): 1250-1261. 121.Dunger G, Pereda R, Farah C, Orellano E, Jorgelina O: Protein-protein interactions identified for effector proteins of the phytopathogen Xanthomonas axonopodis pv. citri [abstract] Proceedings of the V Congreso Argentino de Microbiologia General 2008 [http:// www.conicet.gov.ar/scp/ vista_resumen.php?produccion=637809&id=97&keywords=carri ]. 122.Gurlebeck D, Thieme F, Bonas U: Type III effector proteins from the plant pathogen Xanthomonas and their role in the interaction with the host plant Journal of Plant Physiology 2006, 163(3): 233-255. 123.Furutani A, Takaoka M, Sanada H, Noguchi Y, Oku T, Tsuno K, Ochiai H, Tsuge S: Identification of Novel Type III Secretion Effectors in Xanthomonas oryzae pv. oryzae Molecular Plant-Microbe Interactions 2009, 22(1): 96-106. 124.Roden JA, Belt B, Ross JB, Tach ibana T, Vargas J, Mudgett MB: A genetic screen to isolate type III effectors translocated into pepper cells during Xanthomonas infection Proceedings of the National Academy of Sciences of the United States of America 2004, 101(47): 16624-16629. 125.Metz M, Dahlbeck D, Morales CQ, Al Sady B, Clark ET, Staskawicz BJ: The conserved Xanthomonas campestris pv. vesicatoria effector protein XopX is a virulence factor and suppresses host defense in Nicotiana benthamiana Plant Journal 2005, 41(6): 801-814. 126.Guidot A, Prior P, Schoenfeld J, Carrere S, Genin S, Boucher C: Genomic structure and phylogeny of the plant pathogen Ralstonia solanacearum inferred from gene distribution analysis Journal of Bacteriology 2007, 189(2): 377-387. 127.Park DS, Hyun JW, Park YJ, Kim JS, Kang HW, Hahn JH, Go SJ: Sensitive and specific detection of Xanthomonas axonopodis pv. citri by PCR using pathovar specific primers based on hrpW gene sequences Microbiological Research 2006, 161(2): 145-149. 128.Nol L, Thieme F, Nennstiel D, Bonas U: cDNA-AFLP analysis unravels a genome-wide hrpG -regulon in the plant pathogen Xanthomonas campestris pv. vesicatoria Molecular Microbiology 2001, 41(6): 1271-1281. 129.Xu RQ, Li XZ, Wei HY, Jiang B, Li K, He YQ, Feng JX, Tang JL: Regulation of eight avr genes by hrpG and hrpX in Xanthomonas campestris pv. campestris and their role in pathogenicity Progress in Natural Science 2006, 16(12): 1288-1294. 130.Astua-Monge G, Minsav age GV, Stall RE, Davis MJ, Bonas U, Jones JB: Resistance of Tomato and Pepper to T3 Strains of Xanthomonas campestris pv. vesicatoria Is Specified by a Plant-Inducible Avirulence Gene Molecular Plant-Microbe Interactions 2000, 13(9): 911-921.doi: 10.1186/1471-2164-11-238 Cite this article as: Moreira et al. Novel insights into the genomic basis of citrus canker based on the genome se quences of two strains of Xanthomonas fuscans subsp. aurantifolii BMC Genomics 2010, 11 :238
Xanthomonas citri subsp. citri Xanthomonas fuscans subsp. aurantifolii C Xanthomonas fuscans subsp. aurantifolii B Xanthomonas citri subsp. citri Xanthomonas fuscans subsp. aurantifolii B Xanthomonas fuscans subsp. aurantifolii C
Table S2 Genes shared by XAC, XauB, and XauC but that are not found in other fully sequenced Xanthomonas and Xylella genomes Families were obtained with OrthoMCL (Li, Stoeckert, and Roos, Genome Research 2003, 13:2178 2189) The asterisk indicates th at a nearby gene when mutated caused a phenotype with reduced virulence (Laia et al., BMC Microbiology 2009, 9(12)). Syntenic blocks are indicated by shaded rows Each such block is in a different color. Anomalous regions were determined with AlienHunter results (Vernikos and Parkhill, Bioinformatics 2006 22:2196 2203) on the XAC chromosome. Family XAC Xau B Locus_tag Xau C Locus_tag Function assignment (XAC annotation) Anomalous region ID Strain Start pos. End pos. 4541 XAC0085 103084 104085 XAU B_16700 XAUC_23720 Conserved hypothetical protein YES 4540 XAC0086 104350 104865 XAUB_16690 XAUC_23730 Conserved hypothetical protein YES 4538 XAC0134 157496 158680 XAUB_16030 XAUC_12950 Beta lactamase NO 4537 XAC0283 336818 337519 XAUB_36980 XAUC_37 550 Hydrolase NO 4536 XAC0291 344924 347917 XAUB_37060 XAUC_37630 Oar protein NO 4535 XAC0527 617360 617854 XAUB_33570 XAUC_15720 Conserved hypothetical protein NO 4534 XAC0599 702983 703486 XAUB_23120 XAUC_21280 Conserved hypothetical protein NO 45 32 XAC0726 865373 866395 XAUB_38810 XAUC_30830 Methyl parathion hydrolase NO 4531 XAC0727 866458 867252 XAUB_38820 XAUC_30840 Transcriptional regulator NO 4530 XAC0845 1004256 1005383 XAUB_21650 XAUC_17170 Transcriptional regulator NO 4529 XAC0846 1 005433 1006620 XAUB_21640 XAUC_17180 FMNH2 dependent monooxygenase NO 4528 XAC0847 1006745 1007503 XAUB_21630 XAUC_17190 ABC transporter ATP binding protein NO 4527 XAC0848 1007500 1008351 XAUB_21620 XAUC_17200 ABC transporter permease NO 4526 XAC0849 1008360 1009364 XAUB_21610 XAUC_17210 Sulfonate binding protein NO 4525 XAC0850 1009376 1010539 XAUB_21600 XAUC_17220 Alkanesulfonate monooxygenase NO 4524 XAC0851 1010750 1011313 XAUB_21590 XAUC_17230 NADH dependent FMN reductase NO 4523 XAC0852 1 011721 1013925 XAUB_21580 XAUC_09010 TonB dependent receptor NO 4522 XAC0854 1014835 1016139 XAUB_03080 XAUC_09000 Conserved hypothetical protein NO 4521 XAC0855 1016136 1017596 XAUB_03070 XAUC_08990 Monooxygenase NO 4520 XAC0856 1017611 1019227 XAUB _03060 XAUC_08980 ABC transporter binding protein NO 4519 XAC0858 1020198 1021082 XAUB_18530 XAUC_08960 ABC transporter permease NO 4518 XAC0860 1021286 1022728 XAUB_18540 XAUC_08950 ABC transporter ATP binding protein NO 4516 XAC1388 1595393 1598989 XAUB_20390 XAUC_10530 Conserved hypothetical protein NO
4515 XAC1389 1599002 1599895 XAUB_20400 XAUC_10520 ABC transporter ATP binding protein NO 4512 XAC1613 1859901 1860470 XAUB_24450 XAUC_12090 Conserved hypothetical protein NO 4511 XAC1614 18606 74 1861363 XAUB_24460 XAUC_19550 Conserved hypothetical protein NO 4509 XAC2119 2476483 2477103 XAUB_14670 XAUC_07620 Conserved hypothetical protein NO 4500 XAC2673 3133802 3134287 XAUB_08600 XAUC_12210 Conserved hypothetical protein NO 4499 XAC2860 3355194 3357104 XAUB_06330 XAUC_19060 Conserved hypothetical protein NO 4498 XAC2862 3358551 3359192 XAUB_03690 XAUC_43420 Conserved hypothetical protein NO 4497 XAC3050 3567839 3570760 XAUB_33790 XAUC_34090 TonB dependent receptor YES 4428 XAC3077 3607286 3610177 XAUB_34020 XAUC_12600 TonB dependent receptor NO 4496 XAC3079 3612182 3613549 XAUB_34030 XAUC_12580 Cationic amino acid transporter NO 4495 XAC3080 3613528 3614463 XAUB_34040 XAUC_12570 Ribokinase NO 4494 XAC3081 3614460 3615884 XAUB _34050 XAUC_12560 6 phospho beta glucosidase NO 4493 XAC3088 3625907 3626212 XAUB_34110 XAUC_12500 Conserved hypothetical protein NO 4492 XAC3131 3682724 3683038 XAUB_29300 XAUC_14480 Conserved hypothetical protein NO 4491 XAC3183 3746178 3747353 XAU B_16920 XAUC_01580 Conserved hypothetical protein NO 4490 XAC3197 3762283 3763077 XAUB_10950 XAUC_29340 ABC transporter permease NO 4489 XAC3198 3763196 3764221 XAUB_10940 XAUC_29350 Alkanesulfonate transporter substrate binding NO 4488 XAC3199 3764 320 3765360 XAUB_10930 XAUC_29360 oxidoreductase NO 4487 XAC3200 3765377 3766705 XAUB_10920 XAUC_29370 Nitrilotriacetate monooxygenase component A NO 4486 XAC3224 3797831 3798901 XAUB_14680 XAUC_00040 Avirulence protein (type III SS XopE3 ) NO 4485 XA C3230 3806072 3806962 XAUB_26830 XAUC_23780 Avirulence protein (type III SS XopAI) YES 4484 XAC3233 3808805 3809920 XAUB_03570 XAUC_00450 Transposase YES 4483 XAC3234 3809875 3810390 XAUB_20110 XAUC_00440 Conserved hypothetical protein YES 4482 XAC32 70 3851918 3852886 XAUB_29760 XAUC_00380 Conserved hypothetical protein YES 4477 XAC3337 3933209 3933577 XAUB_15050 XAUC_30350 Conserved hypothetical protein NO 4475 XAC3588 4254271 4255149 XAUB_05370 XAUC_12710 Integral membrane protein BORDER 4426 XAC3589 4255252 4256640 XAUB_05360 XAUC_12700 Conserved hypothetical protein BORDER 4025 XAC3590 4256640 4257941 XAUB_05350 XAUC_12680 Oxidoreductase BORDER 4474 XAC3591 4257944 4258672 XAUB_05340 XAUC_12670 Short chain dehydrogenase BORDER 4473 XAC3 593 4261220 4262164 XAUB_05330 XAUC_12660 NAD dependent epimer /dehydrat /dehydrogen BORDER 4472 XAC3594 4262187 4263473 XAUB_02700 XAUC_12650 Conserved hypothetical protein YES 4471 XAC3595 4263470 4263868 XAUB_02690 XAUC_12640 Conserved hypothetic al protein YES 4470 XAC3611 4283663 4285507 XAUB_16490 XAUC_04090 Peptidase NO 4469 XAC3619 4291954 4292916 XAUB_11810 XAUC_06250 Conserved hypothetical protein NO 4468 XAC3620 4292976 4295219 XAUB_11820 XAUC_06260 Outer membrane receptor FepA NO 44 67 XAC3767 4443316 4443744 XAUB_02890 XAUC_42890 Conserved hypothetical protein NO 4466 XAC3874 4557246 4557638 XAUB_25780 XAUC_27320 Conserved hypothetical protein NO
4465 XAC3968 4665433 4665855 XAUB_28410 XAUC_07270 Conserved hypothetical protein N O 4463 XAC4104 4806727 4807071 XAUB_28690 XAUC_17760 Conserved hypothetical protein NO 4458 XAC4253 5020301 5020702 XAUB_19700 XAUC_34520 Conserved hypothetical protein NO 4457 XAC4326 5111104 5114751 XAUB_03730 XAUC_14760 Urea amidolyas NO 4456 XAC 4327 5114835 5116634 XAUB_10310 XAUC_14750 Amidase NO 4452 XACa0041 99828500 99922440 Partition protein A NO
1 Table S 3. List of XAC specific regions XACSR Region Number; Start LT: locus tag of first protein coding gene in the region; End LT: locus tag of last protein coding gene in the region; Start Pos.: left border of region (nucleotide position); End Pos.: r ight border of region (nucleotide position); NG total number of protein coding genes in region; NG S number of XAC specific genes ( percent of NG ); N G M number o f mobile elements ; NG H number of genes without functional assignment; TGR number of gene s represented in plate (percent of NG ); TGRR>2 number of genes in region with hybridization ratio greater than 2 (percent of TGR)
2 XAC SR Genetic composition of the Regions tRNA Transposase Integrase XACarray Start LT End LT Start Pos. End Pos. N G NG S (%) N G M NG H 3 TGR (%) TGRR>2 (%) 1 0037 0063 46368 77539 27 23 (85) 0 12 ------No No No No 19 (70) 6 (32) 2 0146 0152 176973 184553 7 6 ( 85 ) 3 4 ------No Yes No No 5 (71) 3 (60) 3 1051 1072 1209505 1228677 22 20 (91) 10 11 PRH KL Yes Yes No No 15 (68) 7 (47) 4 1101 1107 1253185 1259423 7 5 (71) 4 2 ---S Yes No No Yes 7 (100) 1 (14) 5 1491 1509 1723743 1744015 20 9 (45) 6 13 ------No No No Yes 12 (60) 4 (30) 6 1657 1664 1908975 1919214 8 6 (75) 4 4 S --No No No Yes 6 (75) 3 (50) 7 1810 1819 2083535 2113580 10 8 (80) 0 2 R ---No No No No 6 (60) 0 (00 8 1866 1872 2164332 2172365 7 5 (71) 2 4 ------No Yes No No 5 (71) 1 (20) 9 1911 1929 2237340 2255111 19 15 (79) 4 5 ------Yes Yes No No 14 (74) 4 (29) 10 2174 2286 2547786 2675160 113 96 (84) 17 58 ------Yes No No Yes 70 (61) 35 (50) 11 2418 2445 2818642 2847010 28 16 (57) 12 13 A ---Yes No No No 19 (68) 10 (53) 12 2601 2609 3068262 3075701 11 9 (69) 2 6 ---V Yes No No No 9 (69) 8 (89) 13 2632 2637 3096409 3103662 6 5 (83) 2 4 ------No No Yes Yes 4 (67) 0 (0) 14 2660 2673 3119404 3134287 14 7 (50) 4 3 ------Yes No No No 10 (71) 2 (20) 15 2901 2904 3406101 3411288 4 4 (100) 1 3 ------No No No Yes 4 ( 100) 0 (0) 16 3014 3025 3524613 3532867 12 12 (100) 0 8 ------No No No No 9 (75) 2 (22) 17 324 5 3298 382 3370 3877600 54 3 6 (67) 7 29 ---G No Yes No Yes 4 3 (81) 9 (21) 18 359 6 3599 4263865 4270247 4 3 (75) 0 1 ------No No No No 3 ( 75 ) 2 ( 66 ) 19 3770 3785 4447486 4454908 16 14 (87) 1 12 ---M No No No Yes 5 (63) 1 (20) 20 3836 3843 4512469 4522281 8 8 (100) 0 8 ------No No No No 5 (31) 1 (20) 21 3932 3955 4624144 4651896 24 16 (66) 5 14 ------No No Yes No 12 (50) 4 (33) 2 2 4061 4064 4756620 4761806 4 4 (100) 0 3 ------No No No No 4 (100) 2 (50) 23 4135 4139 4858635 4865544 5 5 (100) 2 3 ------No Yes No No 2 (50) 0 (0) 24 4261 4266 5030510 5035413 6 6 (100) 0 6 ------No No No No 3 (50) 2 (66) 25 4314 43 38 5099881 5134021 25 19 (76) 4 12 ------No No No No 14 (56) 1 (7)
1 XAC 0037 0063 (27,23,0,12) # # # # PILI 2 XAC 0146 0152 (7,6,3,4) # XAC 1051 1072 PRH KL (22,20,10,11) 3 # # 4 XAC S 1101 1107 (7,5,4,2) # # (20,9,6,13) # 5 XAC 1491 1509 # # # # # # # # # # 6 XAC S 1657 1664 (8,6,4,4) # # XAC X X 1810 1819 R (10,8,0,2) 7 # # T1SS 8 XAC 1861 1883 RPF CLUSTER (7,5,2,4) # # XAC 1911 1929 X (19,15,5,5) 9 # # # # 10 XAC X 2174 2196 XAC XAC XAC X 2197 2213 2214 2249 2250 2274 XAC 2275 2286 (113,96,17,58) # # # # # # # # # # # # # # # # # # T1SS 11 XAC A X 2418 24 45 (28,16,12,13) # # # # # # # # # # # # M obile G enetic E lements (MGE) (transposases/phages) Gene with assigned function H ypotetical or C onserved H ypothetical (HCH) genes Adaptation Pathogenicity Virulence # X Frameshift or Point mutation Low similarity in XauB and XauC Genomes Actual XAC specific region with flanking regions of interest LEGEND (X,Y,Z,W) X total of genes by region Y unique genes by region Z MGE by region W HCH by region
12 XAC X X 2600 2611 V (11,9,2,6) T4SS # # # 14 XAC X 2660 267 3 PILI (14,7,4,3) # # # # # # # 15 XAC 2898 2904 (4,4,1,3) RM I XAC X 2628 2638 2659 (6,5,2,4) 13 N # (28,0,17,10*) 16 XAC 3014 3025 (12,12,0,8) 17 XAC X 3245 3272 (54,35,7,29) XX G VIRULENCE REGULATOR 3273 3298 XAC # # # # # # # # # # # # # # # # # 21 XAC X X A 3930 3960 # # ## # ## # (24,16,5,14) 18 XAC 3596 3599 (4,4,0,3) 19 XAC X X X M 3770 3785 (16,14,1,12) # # 20 XAC 3836 3843 (8,8,0,8) 22 XAC 4061 4064 (4,4,0,3) 25 XAC 4314 4338 (25,19,4,12) RM I # # # X X # # # 23 XAC 4135 4139 (5,5,2,3) 24 XAC 4261 4266 (6,6,0,6) M obile G enetic E lements (MGE) (transposases/phages) Gene with assigned function H ypotetical or C onserved H ypothetical (HCH) genes Adaptation Pathogenicity Virulence # X Frameshift or Point mutation Low similarity in XauB and XauC Genomes Actual XAC specific region with flanking regions of interest LEGEND (X,Y,Z,W) X total of genes by region Y unique genes by region Z MGE by region W HCH by region rhsD
MW Xac XauB XauC Xac XauB XauC Xac XauB XauC MW 100 ng 50 ng 25 ng
Additional file 6 Table S 6 : gene sequences (nucleotides) used to build the pth gene phylogeny. organism Gene NCBI accession number NCBI GeneID label in tree XAC pthA1_XACa0022 1158467 XACa0022 XAC pthA2_XACa0039 1158484 XACa0039 XAC pthA3_XACb00 15 1158502 XACb0015 XAC pthA4_XACb0065 1158552 XACb0065 XauB XAUB_40130 XAUB_40130 XauB XAUB_28490 XAUB_28490 XauC XAUC_22430 XAUC_22430 XauC XAUC_24060 + XAUC_09900 XAUC_24060_09900 X. campestris pthN AF016221 malvacearum pthN X. campes tris pv. malvacearum avrB6 L06634 malvacearum avrB6 X. citri pv. citri strain X0053 PthAW EF473085 citri Aw_X0053 X. gardneri avrHah1 EF436255 gardneri X. oryzae pv. oryzae KACC10331 pthA1 NC_006834.1 3262051 pthA1 Xoo X. oryzae pv. oryzae KACC1 0331 pthA4 NC_006834.1 3265873 pthA4 Xoo X. oryzae pv. oryzae KACC10331 pthA2 NC_006834.1 3264948 pthA2 _Xoo X. oryzae pv. oryzae KACC10331 avrB6 NC_006834.1 3263649 avrB6 _Xoo X. oryzae pv. oryzae KACC10331 pthA3 NC_006834.1 3263519 pthA3 _Xoo X. cit ri apl3 AB021365 apl3 citri X. campestris pv. armoraciae Hax2 AY993937 Hax2_ armoraciae X. citri pv. citri Xc270 pthA* EF473086 citri A* Xc270 X. axonopodis manihotis pBsF2 pthB AF012325 manihotis pBsF2
X. smithii subsp. citri pthA KC21 AB2063 88 pthA KC21 _smithii Ralstonia solanacearum AVRBS3 LIKE PROTEIN CAD15517.1 17428832
malvacearum avrB6 citri A* Xc270 apl3 citri XACb0065 XACa0022 Pth KC21 smithii XACb0015 Citri Aw X0053 XACa0039 Topology and boostrap values only; branch lengths do not reflect phylogenetic distance. 68 56 77 65 72 36 13 16