Comparative genomic and transcriptome analyses of pathotypes of Xanthomonas citri subsp. citri provide insights into mechanisms of bacterial virulence and host range

Multi locus sequencing typing analysis

Multi locus sequence typing (MLST) based phylogenetic analysis was performed for Xcaw12879 and other Xanthomonas spp. using nine housekeeping genes (uvrD, secA, carA, recA, groEL, dnaK, atpD, gyrB, and infB) that are highly conserved in bacteria. The nine protein sequences were aligned and concatenated and then used to construct a maximum-likelihood phylogenetic tree (Figure 1). The results showed that Xcaw12879 is closely related to XccA306. Interestingly, these two citrus canker pathogens form a clade with X. citri pv. mangiferaeindicae LMG 941 and X. axonopodis pv. punicae LMG 859, which cause bacterial black spot in mango and bacterial leaf blight in pomegranates respectively. Both strains were isolated from India [13, 14], a putative origin of XccA. Hence, it is possible that these pathogens have evolved from the same ancestor and evolved to adapt to different hosts. The XccA306 and Xcaw12879 strains share close relationship with the other two citrus canker causing bacteria XauB and XauC (Figure 1). The close relationship between XccA, Xcaw, XauB and XauC agrees with the genome-based phylogeny of the genus Xanthomonas[15].

Chromosome organization and genome plasticity

Whole-genome alignment of Xcaw12879 to closely related XccA306 using MAUVE in progressive mode revealed numerous inversions and translocations (Figure 2). Most of the separated blocks in the alignment are associated with integrases and/or IS elements on at least one of their borders. The IS elements have been known to aid horizontal gene transfer and other genome rearrangements [16].

Xcaw12879 genome contains two plasmids pXcaw19 and pXcaw58 that are significantly different from the plasmids found in XccA306. Plasmid pXcaw19 sequence has no similarity with the plasmids of XccA306, whereas pXcaw58 is only about 35% similar to pXAC64. Plasmid pXcaw58 contains the pthAw2 gene, a homolog of pthA4, which is capable of conferring the ability to cause canker-like symptoms [17]. However, the plasmid pXcaw58 does not contain the Vir like type IV secretion system genes found on pXAC64. The type IV secretion system has been shown to contribute to virulence in X. campestris pv. campestris strain 8004 [18] and absence of these genes from the plasmid could affect virulence of Xcaw12879 strain.

Three clustered regularly interspaced short palindromic repeats (CRISPRs) with short (21–47 bp) direct repeats interspaced with unrelated similarly sized non-repetitive sequences (spacers) are found in Xcaw12879 genome (Additional file 1). The CRISPR1 and CRISPR2 repeats are also present in XccA306. CRISPR2 and CRISPR3 from Xcaw12879 are identical except for a G at the beginning of CRISPR2, indicating that it might be a recent duplication. CRISPR is a bacterial immunity system that helps exclude foreign genetic elements. However the variability in Xcaw12879 and XccA306 suggests that the strains might have had dissimilar exposure to foreign genetic material as suggested in X. oryzae[19].

The TBLASTN analysis of all the proteins from Xcaw12879 and XccA306 revealed various gene clusters specific to each strain. Of the 4,760 proteins from Xcaw12879 and 4,603 (176 not annotated previously [8]) proteins from XccA306, 4,428 proteins are found to be orthologous using the cut-off e-value ≤ 10^-10 and alignments >60% sequence identity, >60% query gene length. Among the 4,428 common proteins, 4,252 were annotated in XccA306 [8] whereas 176 are not annotated. Xcaw12879 has 332 proteins that are either non-orthologous to proteins from XccA306 or unique, whereas XccA306 has 175 such proteins.

The hrp and hrc genes encoding the type 3 secretion system (T3SS) in Xcaw12879 are homologous to the hrp and hrc genes found in XccA306. All the genes are found in similar order with the exception in gene annotation between hrpF and hpaB. The genome of XccA306 contains the annotated gene XAC0395 between the two, which is a hypothetical protein. The annotation in Xcaw12879 in the same region is on the opposite strand and contains hpaI (XCAW_00803) and xopF1 (XCAW_00804/XCAW_00805) which may be nonfunctional due to a frameshift. The nucleotide sequences in both strains are the same and the differences in annotation were confirmed by BLAST similarity of the annotated genes in Xcaw12879 to other xanthomonads.

Multiple genes clustered into ten groups were identified in Xcaw12879 but not in XccA306 (Table 2). Many genes of these clusters present in Xcaw12879 but not in XccA306 have homologs in other Xanthomonas species. All these regions contain transposase, integrase or phage related genes indicating possible acquisition by horizontal gene transfer. An interesting difference noted in the above-mentioned regions is in cluster 5, which encodes for lipopolysaccharide (LPS) biosynthetic pathway. Interestingly, the LPS cluster in Xcaw12879 contains regions orthologous to both XccA306 and X. oryzae pv. oryzicola BLS256 as shown in Figure 3, which indicates that there has been horizontal gene transfer. Cluster 4 is almost 100 kb long and parts of cluster 4 are syntenic with regions from X. campestris pv. campestris 8004, a black spot pathogen of cabbage (Table 2). A MUMmer comparison between cluster 4 and X. campestris pv. campestris 8004 shows high synteny (Additional file 2). Three transcriptional regulators (XCAW_01037, XCAW_01129, XCAW_01131) and one two-component system (TCS) sensor kinase (XCAW_01148) and its response regulator (XCAW_01150) are present in Xcaw12879, but are absent in XccA306.

Cluster number	Locus tag	Homologs in other genomes	Function
1	XCAW_01029 to XCAW_01069		hypothetical proteins, RhsA family protein, transcriptional regulator, integrase, adenine specific DNA methylase, type III restriction enzyme: res subunit, ATP dependent exoDNAse, thermonulease
2	XCAW_01117 to XCAW_01151	Some present in X. campestris pv. campestris ATCC 33913	transcriptional regulator, phage-related tail proteins, TCS response sensor and regulator, chitinase, Zn peptidase, transcriptional repressor, protein -glutamate methylesterase
3	XCAW_01571 to XCAW_01582	Some present in X. campestris pv. campestris str. 8004	phage related proteins, hypothetical protein
4	XCAW_01620 to XCAW_01726	Homologous to Acidovorax sp. JS42 and X. campestris pv. campestris str. 8004	transposases, hypothetical proteins, type II restriction enzyme: methylase subunit, phage related regulatory proteins, chromosome partitioning related protein, soluble lytic murein transglycosylase, VirB6 protein
5	XCAW_04292 to XCAW_04303	Homologous to X. oryzae pv. oryzicola BLS256	lipopolysaccharide biosynthesis genes
6	XCAW_04482 to XCAW_04496		transposases, hypothetical proteins, transcriptional repressor, polymerase V subunit
7	XCAW_04519 to XCAW_04545	Homologous to X. oryzae pv. oryzae PXO99A	phage related proteins, transcriptional regulator, transposases, hypothetical proteins
8	XCAW_a00001 to XCAW_a00017		Plasmid partition protein, conjugal transfer protein, hypothetical proteins
9	XCAW_b00002 to XACW_b00018		Transposases, plasmid stability proteins, avirulence protein, hypothetical proteins
10	XCAW_b00048 to XACW_b00056		Plasmid mobilization proteins, transposases, hypothetical proteins

PthA and homologs

All the citrus canker causing variants (XccA, Xcaw, XccA*, XauB, and XauC) contain PthA or its functional homologs. Thus, PthA or its functional homologs is likely the pathogenicity determinant of citrus canker pathogen as suggested in a previous study [17] that linked the strains of Xanthomonas with different host range together. Al-Saadi et al. [17] have shown that all the variants carry one pthA homolog with 17.5 repeats which determines pathogenicity on citrus and triggers immunity in various other plant species [42]. The avrBs3/pthA family of effectors includes various pth genes but only PthA [42] is known to induce canker. The functional homolog of this gene in XccA strain 306 is pthA4, which also has three other paralogs on its two plasmids (Table 1). We found two homologs pthAw1 and pthAw2 in Xcaw12879 genome, both located on plasmid pXcaw58. The pthAw2 gene is 99% identical to pthA4 from XccA and also to pthAw sequenced from another Wellington strain 0053 that is able to complement a knockout mutant of pthA in XccA strain 3213 [17], indicating that PthAw2 is the functional homolog of pthA in Xcaw. PthAw2 has the same repeat number (17.5) as the functional homologs PthA4, PthB and PthC from the three respective citrus canker causing strains XccA, XauB, and XauC [10]. The other homolog PthAw1 in Xcaw has 18.5 tandem repeats, which is different from PthA homologs found in XccA that have either 15.5 or 16.5 tandem repeats. The AvrBs3/PthA family of effectors are known as transcription activator-like (TAL) effectors since they reprogram host cells by specifically binding to the promoters of plant genes recognized by the central domain of tandem repeats [43]. Comparing the DNA binding TAL effector codes for PthA from XccA as predicted by Boch et al. [44] to PthAw indicate that the codes for PthA4 and PthAw2 are quite divergent (Additional file 3). Al-Saadi et al. [17] predicted that the well conserved sequence of the 17^th repeat in functional PthA might be important for pathogenicity on citrus, and this sequence is present in PthAw2. The rest of PthAw2 sequence however potentially encodes a DNA binding code that is only about 79% similar to the one encoded by PthA4 of XccA306 (Additional file 3). This may result in recognition of different target genes in host plant or differences in strength of induction of plant genes and thus affect virulence of Xcaw and XccA.

Pathogenicity and growth assays

All three host limited strains of Xanthomonas affecting citrus, Xcaw, XauB and XauC had avrGf1 and xopAF genes in their genomes. The gene avrGf1 has been previously studied in Xcaw and is known to be responsible for HR in grapefruit [7]. However its effect on other varieties of citrus such as sweet orange is unknown. Also, since xopAF is the other putative effector gene, its effect on host limitation was further characterized by pathogenicity and growth assays of XcawΔxopAF and XcawΔavrGf1 ΔxopAF.

Pathogenicity assays indicated that Xcaw12879 did not elicit any reaction on Valencia or Hamlin oranges at our test conditions while wild type XccA306 caused necrotic raised lesions, typical of citrus canker on the leaves at a high bacterial inoculation concentration of 10⁸ cfu/ml (Figure 4). Xcaw12879 showed a HR on grapefruit leaves that was abolished by deleting avrGf1 gene (XcawΔavrGf1), however the growth of the mutant was visibly reduced compared to XccA306 strain. XcawΔavrGf1 did not show any symptoms or reaction on either Valencia or Hamlin (Figure 4).

To check whether mutation of xopAF affects Xcaw12879 growth in planta, the wild-type strains of XccA306 and Xcaw12879, XcawΔxopAF, XcawΔxopAF-53:xopAF (complementary strain), XcawΔavrGf1 and XcawΔxopAF ΔavrGf1 mutant strains were inoculated into grapefruit, Mexican lime and Valencia leaves. As shown in Figure 5A, the population of Xcaw12879 was much lower compared to XccA306 in grapefruit. This population of XcawΔavrGf1 was increased compared to the wild type Xcaw12879 and XcawΔavrGf1 caused symptoms on grapefruit. However, the populations of XcawΔxopAF and XcawΔxopAF ΔavrGf1 mutants were one order magnitude lower than that of Xcaw12879 and XcawΔavrGf1 respectively, indicating that mutation of xopAF gene decreased the growth of Xcaw12879 in planta. A similar trend was observed in Mexican lime where the populations of xopAF single and xopAF, avrGf1 double mutants were lower compared to Xcaw12879 and XcawΔavrGf1 respectively (Figure 5B). The growth of XcawΔxopAF in grapefruit and Mexican lime was restored to similar levels as Xcaw12879 by the complementation (Figure 5). No significant changes were observed in Valencia leaves as neither Xcaw12879 nor any of its mutants grew well in the sweet orange variety as compared to XccA306 (Figure 5C).

Transcriptome analyses of Xcaw12879 and XccA306 under nutrient rich (NB) and hrp gene expression inducing (XVM2) conditions

To determine the differential gene expression amongst the strains of X. citri subsp. citri, we grew Xcaw12879 and XccA306 under nutrient rich condition in Nutrient Broth (NB) and in XVM2 [12]. Three biological replicates of the strains were used for RNA-Seq. Over 45 million reads were obtained on average for each sample. After trimming and mapping, approximately 96% of the reads were mapped to the genomes (data not shown) indicating that RNA-Seq provides high quality reads suitable for Xanthomonas transcriptomics. Of all the reads, over 6.5 to 14 million reads could be mapped from each sample to mRNA specifically (Additional file 4). This gave an enrichment of mRNA from 11.3% up to 28.5% for each sample. It has been suggested that 5–10 million non-rRNA fragments enable profiling of the vast majority of transcriptional activity in diverse species including E. coli grown under diverse culture conditions [45]. It was also found that when RNA-Seq data from biological replicates is available, differential expression of numerous genes can be detected with high statistical significance even when the number of fragments per sample is reduced to 2–3 million [45]. Thus our RNA-Seq data is likely sufficient for the transcriptome analysis of XccA306 and Xcaw12879.

To quantify the expression of each gene, the reads aligned to each gene were pooled and normalized for gene size by calculating the Reads Per Kb per Million reads (RPKM) values. The values for each gene from all the replicates were further quantile normalized to test them statistically. The resulting values were log₂ transformed and t-test was performed on these expression values to compare differential gene expression (DGE) between XccA306 and Xcaw12879 under the same growth conditions or between the same strains in NB or XVM2 growth conditions. High correlation was observed between differential expression values of biological replicates (Additional file 5), signifying that the method was reproducible. Principal component analysis indicates that the biological replicates of XccA formed a separate cluster from Xcaw in both growth conditions (Additional file 6).

qRT-PCR was used to validate the RNA-Seq data. Eight genes were chosen (Additional file 7) that were differentially expressed in Xcaw as compared to XccA under both NB and XVM2 growth conditions to compare data obtained from the two methods. The resulting transcriptional ratio from qRT-PCR analysis was log₂ transformed to compare with the DGE values obtained by RNA-Seq (Additional file 8). Although the scale of fold changes between the two techniques is different, high correlation coefficient of 0.87 verifies that the general trend of gene expression is consistent for both data sets.

We studied the expression profile of Xcc strains in XVM2 as compared to NB. At the cut-off of │fold change│ = 3, FDR < 0.05, 292 genes showed differential expression (173 up-regulated and 119 down-regulated in XVM2 compared to NB) in XccA (Additional file 9) and 281 genes (129 up-regulated and 152 down-regulated in XVM2 compared to NB) for Xcaw (Additional file 10). The entire T3SS cluster consisting of twenty-five genes except one gene (XAC0395) was up-regulated in XVM2 for both XccA and Xcaw strains (Additional files 9 and 10). Among all the effectors, sixteen were induced for XccA whereas nineteen effectors were overexpressed for Xcaw in XVM2 compared to in NB. As identified in this study, the effector genes avrBs2, xopA, xopE1, xopE3, xopI, xopX, xopZ1, xopAD, xopAP, xopAQ, hpaA, xopN and xopP were up-regulated in XVM2 in both strains, while pthA1, pthA2, avrXacE3 and xopK were induced only in XccA and xopL, xopR, xopAI, xopAK, xopAF and xopAG only in Xcaw strain.

The 11-gene xps cluster encodes for type 2 secretion system (T2SS) in Xanthomonas secreting various enzymes including pectate lyase, cellulase, and xylanase. The xps genes were down-regulated in XVM2 as compared to in NB for Xcaw, with xpsE being the most significantly down-regulated. For XccA, the xps genes were not down-regulated. Besides the T2SS genes, at least 22 genes encoding T2SS substrates in XccA were overexpressed in XVM2 as compared to only 12 in Xcaw. To the contrary 11 genes for Xcaw and 8 for XccA were down-regulated in XVM2 compared to in NB (Additional files 9 and 10).

Our analysis showed that all the flagella biosynthesis genes encoded by flg and fli, motility by mot and chemotaxis by mcp, che and tsr were repressed in XVM2 for XccA and Xcaw except cheY (XAC3284 in XccA and XCAW_03412 in Xcaw) and tar (XCAW_03417, XCAW_04009 and XCAW_02497). The genes encoding LPS were down-regulated in both Xcaw and XccA, whereas the xanthan gum (EPS) genes were overexpressed in both except gumP in XccA. A few genes encoding outer membrane proteins, which are involved in adhesion, including ompW, blc and hms were up-regulated in XVM2 as compared to in NB for both strains while xadA and yapH were induced in XccA but down-regulated in Xcaw. The Type IV pili genes encoded by pil and fim genes except pilB and filamentous haemagglutinin related genes (fhaB, XAC1816) were down-regulated in both Xcaw and XccA (Additional files 9 and 10).

In order to further understand the molecular mechanisms determining the differences in virulence and host range of Xcaw and XccA, we compared the expression profile of common genes of Xcaw and XccA. When expression of orthologous genes in Xcaw was compared to XccA, 603 genes (426 overexpressed and 177 down-regulated) in NB (Additional file 11) and 450 genes (319 overexpressed and 131 down-regulated) genes in XVM2 (Additional file 12) conditions were significantly differentially regulated at cut-off value of │fold change│ = 3 and FDR < 0.05. On comparing the differentially expressed genes in both conditions, 126 genes were differentially regulated in Xcaw as compared to XccA, irrespective of the growth conditions (Figure 6). Of these 87 were overexpressed in Xcaw and 39 genes were repressed as compared to XccA (Additional file 13). Of the 87 genes overexpressed in Xcaw, 35 were virulence-related genes including hrpX, hrpG, phoP-phoQ regulatory genes, and T2SS substrate genes (XAC2537, XAC2763, XAC2999, XAC4004) (Additional file 13). Of the 39 genes overexpressed in XccA, 21 were virulence-related genes including cellulase genes (XAC0028, XAC0029 and engXCA), reactive oxygen species scavenging enzyme genes, e.g., superoxide dismutase gene sodC2, and genes encoding heat shock protein GrpE and heat stress protein Muc.

Since the gene expression of T2SS substrate genes was different, we compared the protease and pectate lyase activities of XccA306 and Xcaw12879. Xcaw12879 showed higher protease activity than XccA306 (Figure 7A). Xcaw12879 showed lower pectate lyase activity compared to XccA306 (Figure 7B).

Phylogenetic and comparative analysis

The deduced protein sequences of nine housekeeping genes (uvrD, secA, carA, recA, groEL, dnaK, atpD, gyrB and infB) from 13 completely sequenced and 10 draft Xanthomonas spp., and three Xylella fastidiosa strains (out-group species) were used to construct the phylogenetic tree. Amino acid sequences were aligned using ClustalW 2.1 [60]. A phylogenetic tree was constructed from the concatenated sequences using CLC Genomics Workbench v6.0 (CLC Bio, Aarhus, Denmark) by the maximum likelihood method. Comparative analyses of XccA306 and Xcaw12879 was conducted by, a two-way BLAST of the nucleotide sequences to identify unique genes in each strain using the standalone blast + software (ncbi-blast-2.2.4). The genes were considered orthologous if reciprocal TBLASTN hits were found between two genes with e-value less than or equal to 10^-10 and alignments exceeding 60% sequence identity and 60% query gene length. A gene was considered singleton or unique to each strain if it had no hits or with an e-values less than or equal to 10^-5[61, 62]. The CRISPRfinder server [63] was used to identify CRISPRs. Only confirmed structures are reported here. Alignment between whole chromosomes was done using the script Promer from the MUMmer package [64]. Promer does alignments between translated nucleotide sequences.

Preparation of RNA samples for transcriptome analysis

RNA sample preparation and cDNA library generation were performed according to procedures outlined by Filiatrault et al. [65] with some modifications. RNA samples were extracted from XccA306 and Xcaw12879 grown to OD₆₀₀ of 0.4 in XVM2 medium and NB medium at 28°C on shaker at 200 rpm. The starting OD₆₀₀ for each culture was 0.03. Three biological replicates of each strain in each medium were used for RNA extraction. When the OD₅₆₀ reached 0.4 for each condition, RNA was stabilized immediately by mixing 10 ml of the culture with two volumes of RNAprotect bacterial reagent (Qiagen, Valencia, CA). The cells were centrifuged at 5000 × g at 4°C and cell pellets were treated with lysozyme and RNA extractions were performed using RiboPure bacteria kit (Ambion, Austin, TX) per manufacturers” instructions. Genomic DNA was removed by treatment with TURBO DNA-free kit (Ambion, Austin, TX). Total RNA samples were quantified using spectrophotometry (Nanodrop ND-1000, NanoDrop Tech. Inc.). RNA quality was assessed using the Agilent 2100 bioanalyzer (Agilent Technologies, Palo Alto, CA).

MRNA enrichment and library construction

mRNA was enriched from total RNA using MicrobExpress kit (Ambion) to remove the 23S and 16S ribosomal RNAs (rRNAs). Removal of rRNAs was assessed using an Agilent Bioanalyzer. Double stranded cDNA synthesis was performed following the Illumina mRNA Sequencing sample preparation guide (Cat. No. RS-930-1001) in accordance with the manufacturer’s standard protocol. Enriched mRNA was fragmented via incubation for 5 min at 94°C with the Illumina-supplied fragmentation buffer. The first strand of cDNA was synthesized by reverse transcription using random oligo primers. Second-strand synthesis was conducted by incubation with RNAse H and DNA polymerase I. The resulting dsDNA fragments were further end-repaired, and A-nucleotide overhangs were added. After the ligation of Illumina adaptors, the samples were run on a denaturing gel and the band correlating to 200 (±25) base pairs on the denatured DNA ladder was selected. The selected DNA constructs were amplified by PCR using the primers provided in the Illumina library kit. The amplified constructs were purified and the library was validated using Agilent 2100 bioanalyzer.

Illumina sequencing and alignment

Paired-end, 75-cycle sequencing of the libraries was performed using an Illumina GAIIx at Yale Center for Genomic Analysis. The raw sequencing reads were further analyzed using CLC Genomics Workbench v6.0 (CLC Bio, Aarhus, Denmark). The reads were trimmed using the quality score limit of 0.08 and maximum limit of 2 ambiguous nucleotides. The trimmed reads were mapped to the genome and the protein-coding genes of XccA306 (GenBank accession no. NC_003919, NC_003921.3 and NC_003922.1) and Xcaw12879, with the parameters allowing mapping of reads to the genome with up to 2 mismatches. The reads mapped to rRNA and the reads not uniquely mapped were removed from further analysis. The expression levels were evaluated by RPKM method as described by Mortazavi et al. [66].

Differential gene expression analysis

The differential gene expression of the pooled samples from each condition was analyzed using CLC Genomics Workbench v6.0 (CLC Bio, Aarhus, Denmark). RPKM values were normalized using quantile normalization and further log₂ transformed for statistical analysis. Box plots, hierarchical clustering of samples and principal component analysis were done to examine data quality and comparability. A t-test was performed on log₂-transformed data to identify the genes with significant changes in expression between the two growth conditions and between the two strains. The p-values were adjusted for the false discovery rate (FDR) using the Benjamini and Hochberg method [67].

Quantitative reverse transcription - PCR (qRT-PCR)

To verify the RNA-Seq result, qRT-PCR assays were carried out using the same sets of RNA for RNA-Seq analysis. Gene specific primers listed in Additional file 7 were designed to generate sequences of 100–250 bp in length from the XccA306 genome. qRT-PCR was performed for all 3 biological replicates of XccA306 and Xcaw12879 grown in NB and XVM2 on a 7500 fast real-time PCR system (Applied Biosystems) using QuantiTect™ SYBR® Green RT-PCR kit (Qiagen) following the manufacturers’ instructions. 16S rRNA was used as an endogenous control. The fold change of gene expression was calculated by using the formula 2^-ΔΔC_T[68]. The fold change was further log₂ transformed to compare with the RNA-Seq data.

Generation of the xopAF mutant and xopAF, avrGf1 double mutant

To construct the xopAF deletion mutant, the 1096-bp fragment containing entire xopAF gene was amplified using genomic DNA of Xcaw12879 as template and primers xopAFF1 and xopAFR. This resulted in F1, containing a Bam HI restriction site within the xopAF gene. A 422 bp fragment containing 337 bp of xopAF gene and its downstream region was amplified further from F1 using primers xopAFF2-Bam HI and xopAFR (Additional file 7), resulting in F2. Both F1 and F2 were digested with Bam HI and fragments F3 (414 bp) and F4 (500 bp) were gel purified. The fragments were ligated and cloned into pGEM-T easy vector, resulting in the construct named pGEM-ΔxopAF that was confirmed by PCR and sequencing. From pGEM-ΔxopAF, an Apa I-Pst I fragment containing xopAF gene with 192 bp internal deletion was transferred into Apa I-Pst I digested suicide vector pNTPS138, resulting in pNTPS-ΔxopAF. The construct pNTPS-ΔxopAF was transformed into E. coli DH5αλPIR. The construct was purified from E. coli and subsequently transferred into Xcaw12879 and Xcaw12879ΔavrGf1 generated in a previous study [7] by electroporation. Transformants were selected on NA medium supplemented with 50 μg/μl kanamycin. Positive colonies were replicated on both NA plates supplemented with 5% (w/v) sucrose and kanamycin, and only NA and kanamycin. The sucrose sensitive colonies were selected from NA plus kanamycin plate and grown in NB medium overnight at 28°C. The culture was then dilution-plated on NA containing 5% sucrose to select for resolution of the construct by a second cross-over event. The resulting deletion mutant of xopAF and double mutant of xopAF and avrGf1 was confirmed by PCR (data not shown). The complete xopAF and avrGf1 genes were complemented in the single and double mutants using pUFR053 and pUFR034 respectively. The resulting complement strains were Xcaw12879ΔxopAF-53:xopAF, Xcaw12879ΔavrGf1-34:avrGf1 and Xcaw12879ΔavrGf ΔxopAF-34:avrGf1-53:xopAF were used in this study.

Pathogenicity assay

Pathogenicity assays were conducted in a quarantine greenhouse facility at Citrus Research and Education Center, Lake Alfred, FL. XccA306, Xcaw12879, and XcawΔavrGf1 strains were grown with shaking overnight at 28°C in NB, centrifuged down and suspended in sterile tap water and the concentrations were adjusted to 10⁸ cfu/ml. The bacterial solutions were infiltrated into fully expanded, immature leaves of Duncan grapefruit, Valencia sweet orange and Hamlin sweet orange, with needleless syringes [54]. The test was repeated three times with similar results. Disease symptoms were photographed 10 days post inoculation.

Growth assay in planta

XccA306, Xcaw12879, XcawΔxopAF, XcawΔxopAF-53:xopAF, XcawΔavrGf1 and XcawΔxopAF ΔavrGf1 strains were grown with shaking overnight at 28°C in NB, centrifuged down and suspended in sterile tap water and the concentrations we re adjusted to 10⁶ cfu/ml. The bacterial solutions were infiltrated into fully expanded, immature leaves of Duncan grapefruit, Mexican lime and Valencia sweet orange with needleless syringes [54]. To evaluate the growth of various Xcc strains and mutants in these plants 2 inoculated leaves were collected from each plant at 0, 2, 4, 7, 10, 14 and 21 days. 1 cm² leaf disks from inoculated leaves were cut with a cork borer and then ground in 1 ml sterile water. These were serially diluted and plated on NA plates. The bacterial colonies were counted after 3-day incubation at 28°C. The test was repeated three times independently.

Pectate lyase and proteinase assay

XccA306 and Xcaw12879 were grown on nutrient agar at 28°C, then suspended in sterile deionized water to the O.D. of 0.3 at 560 nm. Hildebrand’s medium A, B and C were used to test for pectolytic activity [69]. In short the medium contained bromothymol blue dye, calcium chloride, 2% sodium polypectate and 0.4% agar. The pH was adjusted to 4.5, 7.0 and 8.5 for the medium A, B and C. One μl of the cultures were inoculated onto the plates and incubated at 28°C for 6 days before confirming pitting due to pectate lyase production. 10% skim milk agar was used to test the bacterial protease activity. The cultures were grown and suspended in sterile water as explained above. One μl of the cultures were inoculated onto the skim-milk plates and cultured at 28°C for 6 days to observe protease activity.

Availability of supporting data

The genome sequences of Xanthomonas citri subsp. citri strain A^w12879 are available at GenBank under the accession numbers CP003778, CP003779 and CP003780. The RNA-Seq data from this study are available in the NCBI’s Gene Expression Omnibus database under the accession number GSE41519.