Development and validation of a Xanthomonas axonopodis pv. citri DNA microarray platform (XACarray) generated from the shotgun libraries previously used in the sequencing of this bacterial genome
© Moreira et al; licensee BioMed Central Ltd. 2010
Received: 3 February 2010
Accepted: 27 May 2010
Published: 27 May 2010
From shotgun libraries used for the genomic sequencing of the phytopathogenic bacterium Xanthomonas axonopodis pv. citri (XAC), clones that were representative of the largest possible number of coding sequences (CDSs) were selected to create a DNA microarray platform on glass slides (XACarray). The creation of the XACarray allowed for the establishment of a tool that is capable of providing data for the analysis of global genome expression in this organism.
The inserts from the selected clones were amplified by PCR with the universal oligonucleotide primers M13R and M13F. The obtained products were purified and fixed in duplicate on glass slides specific for use in DNA microarrays. The number of spots on the microarray totaled 6,144 and included 768 positive controls and 624 negative controls per slide. Validation of the platform was performed through hybridization of total DNA probes from XAC labeled with different fluorophores, Cy3 and Cy5. In this validation assay, 86% of all PCR products fixed on the glass slides were confirmed to present a hybridization signal greater than twice the standard deviation of the deviation of the global median signal-to-noise ration.
Our validation of the XACArray platform using DNA-DNA hybridization revealed that it can be used to evaluate the expression of 2,365 individual CDSs from all major functional categories, which corresponds to 52.7% of the annotated CDSs of the XAC genome. As a proof of concept, we used this platform in a previously work to verify the absence of genomic regions that could not be detected by sequencing in related strains of Xanthomonas.
Citrus canker, or cancrosis, is a disease that affects most species of the Citrus genus . Its symptoms can be observed on leaves, fruits or branches and are characterized by small, pointed and spongy pustules that are surrounded by a chlorotic halo on both sides of the leaf. This halo tends to spread quickly through the tissue, increasing its diameter and becoming brownish or darker in color with a rough and salient appearance . The causal agent of cancrosis is the bacterium Xanthomonas axonopodis pv. citri (XAC) , more recently renamed as Xanthomonas citri subsp. citri [4, 5] whose genome was completely sequenced in 2002  and compared with others organisms [6–8]. Data from this sequencing project revealed that the XAC genome contains a circular chromosome of 5.2 Mbp and two plasmids (33 and 64 Kbp), containing a total of 4,489 coding regions .
To perform the complete sequencing of this phytopathogen, shotgun libraries containing 46,462 clones were created and distributed over approximately 500 96-well plates, representing 98% of the genome, with an estimated coverage of seven times its size . The remaining 2% of the genome was sequenced using a cosmid library.
Upon completion of the sequencing project, the clones were stored with the expectation that they would be useful for future functional genomics studies. We here explore one of the possible applications of this library, namely the use of these clones to build a DNA microarray platform for use in studies of gene expression and genotyping . Building on the data provided by the XAC genome project, this tool will enable the exploration of the physiological and biochemical machinery of this organism in distinct environmental situations, both in vitro and in vivo, thus providing clues to the working of Xanthomonas metabolism and the mechanisms of infection and disease.
Clone selection from shotgun libraries
Analysis of clone selection grouped according to the overlap models.
Type of clones
Total of selected clones
Average size of overlap Insert/CDS
Average size of non-overlap Insert/CDS
Amplification, purification and identification of PCR products
The selected clone inserts were amplified through PCR after DNA extraction by the boiling preparation technique . The amplification reactions were performed with an initial denaturation at 95°C for 5 minutes, followed by 40 denaturing cycles at 95°C for 45 seconds, annealing at 50°C for 30 seconds and extension at 72°C for 1 minute in a final volume of 50 μL of amplification solution containing 1.0 μL boiling prep product, 1 U Taq DNA Polymerase™ (Invitrogen), 40 mM dNTPs (Invitrogen) and 10 pmols of each universal oligonucleotide M13 (forward 5'-GTTTTCCCAGTCACGAC-3' and reverse 5'-CAGGAAACAGCTATGAC-3'). In an attempt to increase the amount of the generated product, 1.0 μL of the amplification reaction was added to new amplification solution with the same volume and concentration of reagents and submitted to the amplification protocol described above.
The amplification products were purified with Multiscreen Millipore™ (cat. # MAFB NOB 50) filters plates and diluted in 50 μL of 10 mM Tris pH 8.0 directly into 96-well plates. The purification products were then manually combined in 8 384-well plates; one of these plates contained the hybridization controls (ScoreCard ).
To ensure that each amplified PCR product was in fact the clone selected for each CDS, an end of the insert in each of the clones was sequenced using the universal oligonucleotide primer M13R and the BigDye terminator v3™ kit (Applied Biosystems). These sequences were then compared to the database of sequencing reads and the XAC genome using the BLASTn algorithm .
Fixation of PCR products onto glass slides
Equal volumes of DMSO were added to the purified PCR products in 384-well plates so that the final concentration was approximately 200 to 400 fmol/μL of DNA . Afterwards, DNA was deposited onto Type 7 mirrored glass slides (Amersham Biosciences/GE Healthcare) using a robot (Generation III Microarrays Spotter™ - Amersham Biosciences/GE Healthcare). This robot allowed for the deposition of up to 4,608 DNA samples that were organized into 12 subsets (subarrays) of 384 spots (12 × 32) (Figure 2A). The set of 3,072 probes, which corresponded to the products on the 384-well plates and included positive and negative controls, was placed in duplicate in the two longitudinal halves of the slide. The duplicates were called set A and set B, respectively, and constituted a technical replica on each slide (Figure 2A). After DNA deposition, the slides were submitted to 50 mJ UV light, which promoted the fixation of the DNA onto the support. The slides were stored in desiccators (relative moisture ~5%), ready to be used for any experimental purpose.
XAC cultivation and DNA extraction
The XAC 306 isolate was extracted from Citrus limonia L. Osbeck (rangpur lime)-infected leaves and stored in 30% glycerol at room temperature until use. The strain was sowed to reactivate its growth in NA solid medium (3 g/L of meat extract, 5/L g of peptone, 5 g/L of sodium chloride, 15 g/L of agar) and was incubated at 28°C for 48 hours. Samples of the reactivated strain were minced in 12 mL of NA liquid medium in 25-mL Falcon tubes and maintained under gentle agitation at 28°C for 48 hours. Appropriated antibiotic was used when necessary. After this period, the samples were centrifuged for 30 minutes at 3,000 g at room temperature. The supernatant was discarded, and 12.5 mL of wash buffer (10 mM Tris-HCl, pH 8.8, 3.0 mM KCl, 1.25 mM NaCl) was added to the bacterial pellet. The samples were shaken vigorously using a vortex and then centrifuged for 15 minutes at 3,000 g at room temperature. The bacterial pellet was resuspended again, now using 10 mL of D buffer (25 mM sodium citrate, pH 7.0, 5.0 g/L sarkosyl, 4 M guanidine isothiocyanate), and incubated in a 65°C water bath for 1 hour to promote cell lysis. After the lysis step, 5.0 mL of P buffer (667 mM Tris-HCl, 833 mM NaCl, 83 mM EDTA, pH 8.0) was added. The solution was shaken by inversion and centrifuged for 30 minutes at 3,000 g at room temperature. Aliquots (1.0 mL) of the supernatant were transferred to tubes with a 2.0 mL capacity, 1.0 mL of isopropanol was added and the samples were kept at -20°C. The samples were then homogenized by inversion and centrifuged for 45 minutes at 3,000 g at 4°C. The DNA pellet was resuspended, washed twice with 70% ethanol to promote precipitation and centrifuged for 15 minutes at 3,000 g. Finally, the DNA pellet was resuspended in 200 μL TE buffer (10 mM Tris-HCl, pH 8.0, 0.1 mM EDTA) containing 1 μg/mL of RNaseA and incubated for 1 hour at 37°C. The quality and concentration of the extracted DNA was evaluated on a 0.8% agarose gel stained with ethidium bromide ; the profile of the band obtained for the DNA samples was compared to a molecular weight standard with known sizes and concentrations (Figure 2A).
Preparation of labeled DNA and microarray hybridization
Double-stranded DNA labeled with fluorophore was produced by Cy3-dCTP or Cy5-dCTP™ (Amersham Biosciences/GE Healthcare) incorporation during polymerization from 2 μg XAC total DNA (chromosomal + plasmid) that had been fragmented through the shearing technique (~14 times through a syringe with 8 cm and 18 gauge needles) . The reaction was performed using 1 μL of the DNA polymerase Klenow Fragment™ (Gibco) at a high concentration (40 U/μL) and 500 μg of random nonamer primers . Afterwards, the labeled DNA was purified with Multiscreen™ filter plates (Millipore), and the total amount of incorporated fluorophore was indirectly quantified through absorbance detection at 550 nm for Cy3 and 650 nm for Cy5. To ensure accuracy in the analysis, equivalent quantities of fluorescent DNA were used in the microarray hybridizations, which were performed for 16 hours at 42°C in hybridization buffer (Amersham Biosciences/GE Healthcare) with 50% formamide. According to protocols described previously by Koide and coauthors , after washing, the slides were dried by applying a jet of compressed nitrogen gas. Microarray data was aquired using a Generation III Scanner™ (Amersham Biosciences/GE Healthcare) in order to obtain images of each fluorescent channel, Cy3 and Cy5, for each of the microarray probes (Figure 2A).
Detection, quantification and normalization of fluorescence signals
After the microarrays were scanned, the crude signal intensity of the images was determined using Array Vision 6.0™ software (Image Research/Molecular Dynamics). For each spot on the microarray, which represents a specific probe, the foreground and background intensity were measured.
Data preparation included subtraction of the local noise for each spot, and calculation of average fluorescece density and artifacts removal (MTM Dens). Array Vision software allowed the exclusion of bad pixels (pixels within a section that showed signal intensity either above or under four median absolute deviations (MADs) of the signal intensity of all of the pixels in the section).
Construction of XAC DNA microarrays (XACarray)
The selection of clones from the XAC shotgun libraries allowed for the identification of 4,421 clones, which represented 3,084 individual CDSs of the XAC genome (Figure 1A). The remaining 1,405 CDSs, which would complete the genome, were not selected because they represented very small CDSs, smaller than 200 bp for example (Figure 1B), or because they were not physically found in the shotgun library .
Analysis of the CDSs deposited on the XACarray based on the hybridization signal during the experimental validation of the platform.
Type of clones
Physical analysis (Fixed PCR products)
After the in silico selection, the clones were rearranged in 47, 96 well plates using a Q-Bot robot™ (Genetix, UK), and the inserts were amplified by PCR (see methods). Approximately 5% of the clones were shown to be not viable, and 10% presented an insert with a size different from that which was predicted (data not shown), suggesting the existence of probable errors in clone labeling, a common occurrence in clone library construction, even in commercial libraries . Moreover, another 15% presented negative results with respect to their amplification. For these reasons, it was necessary to re-sequence one end of the insert of each of the previously selected clones in order to confirm the identity of all of the clones. The sequences were compared to the XAC genome, and the identity of the clones was corrected when necessary. Unfortunately, at the end of all the steps described above, approximately 30% of the initially selected clones (1,326) could not be used. Thus, at the end of these steps, expected identity of 2,653 clones were confirmed and this set was used. They were rearranged again, now in 28 96-well plates. Although the number of clones was severely reduced, the microarray contained a set of clones with a reliable identity. However, after the PCR was performed to increase the amount of the probe to be immobilized on the slides, the number of clones was reduced to 2,639, which represented 58.8% of the XAC genome (Table 2). To this set of 2,639 PCR products, 121 products derived from amplifications of genomic DNA sequences using specific oligonucleotides that were previously acquired for other experimental proposals were added; these products represented the sequences of genes of particular interest, including the genes vir, hrp, rpf and gum. Thus, the 2,760 PCR products were purified (probes), rearranged in 8 384-well plates, and used in the construction of microarrays as described above.
It is necessary to highlight that the molecular mass pattern obtained from the pUC19 plasmid not only facilitated the analysis of the amplified product, but also demonstrated that this approach was correct. The banding pattern quantification of some randomly selected PCR products corresponded to the visual quantifications that were performed using only the agarose gels (data not shown).
For an initial validation of the quality of the XACarray, hybridization with XAC genomic DNA labeled with distinct fluorophores was performed. Figure 2A shows the fluxogram of the XACarray protocol and the hybridization images obtained using this probe. Note that the spots on the slide are practically all yellow in color, which visually represents the hybridization of probes labeled with Cy3 and Cy5 in similar quantities. This result can be observed more rigorously from the monochromatic images obtained from laser scanning at 550 and 650 nm, respectively.
The analysis of these images revealed a noticeable signal above the background noise for 88.5% of the CDSs represented in the XACarray, a percentage that is equivalent to 2,365 individual CDSs of the XAC genome (52.7%). Table 2 summarizes the physical features of the 2,670 probes deposited on the XACarray. The probes are represented in Figure 1C in relation to the type of overlap of the insert and the average size of these sequences. The classification of the CDSs regarding their annotation characteristics is shown in Table 2. The CDS are classified into two major groups: CDSs previously annotated with a determined putative function, which represent 63.5% of the XACarray; and CDSs that compose the hypothetical group, in this case exclusive hypothetical and conserved hypothetical in other organisms, which represent the other 36.5% of the composition of the array.
In this study, we developed a new and low-cost method to generate a DNA microarray. Previously, Astua-Monge et al.  used a DNA macroarray, implemented on nylon membranes and containing probes for 279 XAC genes associated with pathogenicity and virulence, to investigate genes differentially expressed in a synthetic medium (XVM2) that simulates in planta conditions. In contrast, our plataform is the first DNA microarray platform for XAC that representes more than 50% of all XAC genome CDSs.
The platform was constructed from shotgun libraries, previously generated for the genomic sequencing project. Although the clones of type II are ideal, because they are restricted to include just an internal fragment of the target gene, we have made an effort to select clones that would not overlap with other genes. In our array, there are at least 1,290 clones (~30% of the genome, not counting the ones for which specific primers were designed) that contain genome fragments that do not overlap any gene besides the target gene they were selected to represent. Also, many genes that have fragments in the same clone will correspond to operons and could, in a first crude analysis, be regarded as a single expression unit, with additional experiments clarifying which gene(s) are indeed relevant. Therefore, our platform can be used to generate lists of putative differentially expressed genes in XAC under different physiological conditions.
On the other hand, the presence of overlapping probes makes it necessary to subsequently validate individual candidate genes using methods with higher specificity and sensitivity, like qRT-PCR. Moreover, under most conditions, some more expensive methods, like designed-probe chips and RNA-Seq will always provide better specificity/sensitivity rates and better reproducibility and such methods should be preferred whenever researchers have access to them. Functional results using this platform were described in a recent work that use this tool to validated putative insertion/deletion regions between different and incomplete genomes sequences of other Xanthomonas species .
We acknowledge the Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP) for granting PhD scholarships to LMM, RFS, PZ and MLM and for funding the project (CAGE).
- Gottwald TR, Graham JH, Schubert TS: Citrus Canker: The Pathogen and Its Impact. APS Net: Peer-Reviewed by Plant Health Progress. 2002Google Scholar
- Brunings AM, Gabriel DW: Xanthomonas citri: breaking the surface. Molecular Plant Pathology. 2003, 4 (3): 141-157. 10.1046/j.1364-3703.2003.00163.x.PubMedView ArticleGoogle Scholar
- Bitancourt AA: O Cancro Cítrico. Biológico. 1957, 23: 110-111.Google Scholar
- 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. 10.1016/j.syapm.2005.03.017.PubMedView ArticleGoogle Scholar
- 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. 10.1016/j.syapm.2006.08.001.PubMedView ArticleGoogle Scholar
- da Silva AC, Ferro JA, Reinach FC, Farah CS, Furlan LR, Quaggio RB, Monteiro-Vitorello CB, Van Sluys MA, Almeida NF, Alves LM: Comparison of the genomes of two Xanthomonas pathogens with differing host specificities. Nature. 2002, 417 (6887): 459-463. 10.1038/417459a.PubMedView ArticleGoogle Scholar
- Moreira LM, de Souza RE, Almeida NF, Setubal JC, Oliveira JCF, Furlan LR, Ferro JA, da Silva ACR: Comparative genomics analyses of citrus-associated bacteria. Annual Review of Phytopathology. 2004, 42: 163-184. 10.1146/annurev.phyto.42.040803.140310.PubMedView ArticleGoogle Scholar
- 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.9. 10.1089/omi.2005.9.43.PubMedView ArticleGoogle Scholar
- Hayward RE, Derisi JL, Alfadhli S, Kaslow DC, Brown PO, Rathod PK: Shotgun DNA microarrays and stage-specific gene expression in Plasmodium falciparum malaria. Mol Microbiol. 2000, 35 (1): 6-14. 10.1046/j.1365-2958.2000.01730.x.PubMedView ArticleGoogle Scholar
- Fleischmann RD, Adams MD, White O, Clayton RA, Kirkness EF, Kerlavage AR, Bult CJ, Tomb JF, Dougherty BA, Merrick JM: Whole-genome random sequencing and as sembly of Haem ophilus influenzae Rd. Science. 1995, 269 (5223): 496-512. 10.1126/science.7542800.PubMedView ArticleGoogle Scholar
- Holmes DS, Quigley M: A rapid boiling method for the preparation of bacterial plasm ids. Anal Biochem. 1981, 114 (1): 193-197. 10.1016/0003-2697(81)90473-5.PubMedView ArticleGoogle Scholar
- Raffelsberger W, Dembele D, Neubauer MG, Gottardis MM, Gronemeyer H: Quality indicators increase the reliability of microarray data. Genomics. 2002, 80 (4): 385-394. 10.1006/geno.2002.6848.PubMedView ArticleGoogle Scholar
- Ausubel FM, Brent R, Kingston RE, Moore DD, Seidman JG, Smith JA, Struhl K: Current protocols in Molecular Biology. 1990, New York:: Greene Publishing Associates and Wiley-Interscience, 3Google Scholar
- Sambrook J, Fritsch EF, Maniatis T: Molecular Cloning - A Laboratory Manual. 1989, New York: Cold Spring Harbor Laboratory Press, 2Google Scholar
- Altschul S, Madden T, Schaffer A, Zhang J, Zhang Z, Miller W, Lipman D: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Research. 1997, 25 (17): 3389-3402. 10.1093/nar/25.17.3389.PubMed CentralPubMedView ArticleGoogle Scholar
- Diehl F, Grahlmann S, Beier M, Hoheisel JD: Manufacturing DNA microarrays of high spot homogeneity and reduced background signal. Nucleic Acids Res. 2001, 29 (7): E38-10.1093/nar/29.7.e38.PubMed CentralPubMedView ArticleGoogle Scholar
- Schriefer LA, Gebauer BK, Qui LQ, Waterston RH, Wilson RK: Low pressure DNA shearing: a method for random DNA sequence analysis. Nucleic Acids Res. 1990, 18 (24): 7455-7456. 10.1093/nar/18.24.7455.PubMed CentralPubMedView ArticleGoogle Scholar
- Koide T, Zaini PA, Moreira LM, Vencio RZ, Matsukuma AY, Durham AM, Teixeira DC, El-Dorry H, Monteiro PB, da Silva AC: DNA microarray-based genome comparison of a pathogenic and a nonpathogenic strain of Xylella fastidiosa delineates genes important for bacterial virulence. J Bacteriol. 2004, 186 (16): 5442-5449. 10.1128/JB.186.16.5442-5449.2004.PubMed CentralPubMedView ArticleGoogle Scholar
- Halgren RG, Fielden MR, Fong CJ, Zacharewski TR: Assessment of clone identity and sequence fidelity for 1189 IMAGE cDNA clones. Nucleic Acids Res. 2001, 29 (2): 582-588. 10.1093/nar/29.2.582.PubMed CentralPubMedView ArticleGoogle Scholar
- Astua-Monge G, Freitas-Astua J, Bacocina G, Roncoletta J, Carvalho SA, Machado MA: Expression Profiling of Virulence and Pathogenicity Genes of Xanthomonas axonopodis pv. citri. Journal of Bacteriology. 2005, 187 (3): 1201-1205. 10.1128/JB.187.3.1201-1205.2005.PubMed CentralPubMedView ArticleGoogle Scholar
- Moreira LM, Almeida NF, Potnis N, Digiampietri LA, Adi SS, Bortolossi LC, da Silva AC, da Silva AM, de Moraes FE, Oliveira JC, de Souza RF, Fancincani AP: Novel insights into the genomic basis of citrus canker based on the genome sequences of two strains of Xanthomonas fuscans subsp. aurantifolii. BMC Genomics. 2010, 11 (1): 238-10.1186/1471-2164-11-238.PubMed CentralPubMedView ArticleGoogle Scholar
This article is published under license to 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.