This Article Abstract Full Text (PDF) Supplemental Research Data Alert me when this article is cited Alert me if a correction is posted Citation Map Services Email this article to a friend Similar articles in this journal Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via HighWire Citing Articles via Google Scholar Google Scholar Articles by Nunes, L. R. Articles by Costa de Oliveira, R. Search for Related Content PubMed PubMed Citation Articles by Nunes, L. R. Articles by Costa de Oliveira, R. Social Bookmarking                   What's this?

LETTER

Microarray Analyses of Xylella fastidiosa Provide Evidence of Coordinated Transcription Control of Laterally Transferred Elements

¹Núcleo Integrado de Biotecnologia, Universidade de Mogi das Cruzes, Mogi das Cruzes, SP 08780-911, Brazil;² Centro de Biologia Molecular e Engenharia Genética, Universidade de Campinas, Caixa Postal 6010, Campinas, SP, Brazil;³ Centro APTA de Citros Sylvio Moreira, Instituto Agronômico de Campinas, Caixa Postal 04, Cordeirópolis, SP 13490, 970, Brazil; ⁴Empresa Brasileira de Pesquisa Agropecuária, EMBRAPA; ⁵Unidade de Biotecnologia, Universidade de Ribeirão Preto, Ribeirão Preto, SP 14096-380, Brazil

	ABSTRACT

Top ABSTRACT RESULTS DISCUSSION METHODS WEB SITE REFERENCES REFERENCES

 
Genetically distinct strains of the plant bacterium Xylella fastidiosa (Xf) are responsible for a variety of plant diseases, accounting for severe economic damage throughout the world. Using as a reference the genome of Xf 9a5c strain, associated with citrus variegated chlorosis (CVC), we developed a microarray-based comparison involving 12 Xf isolates, providing a thorough assessment of the variation in genomic composition across the group. Our results demonstrate that Xf displays one of the largest flexible gene pools characterized to date, with several horizontally acquired elements, such as prophages, plasmids, and genomic islands (GIs), which contribute up to 18% of the final genome. Transcriptome analysis of bacteria grown under different conditions shows that most of these elements are transcriptionally active, and their expression can be influenced in a coordinated manner by environmental stimuli. Finally, evaluation of the genetic composition of these laterally transferred elements identified differences that may help to explain the adaptability of Xf strains to infect such a wide range of plant species.

[Supplemental material is available online at www.genome.org.]

The existence of different disease symptoms observed in a wide range of plant hosts (Pooler and Hartung 1995; Chen et al. 2000), possibly associated with genetically distinct Xf strains, has led researchers to hypothesize that total genome comparisons among such strains may uncover important information regarding genes involved in the interaction with specific hosts. These results could provide the basis for further epidemiological studies and more efficient control of these phytopathogens (Hendson et al. 2001). Thus, the genomes of two other strains isolated from oleander and almond trees were partially sequenced and annotated (http://www.jgi.doe.gov/jgi_microbial/html), and a fourth strain, isolated from grapevines and responsible for PD in California, has also been sequenced to completion (http://www.lbi.ic.unicamp.br/world/xf-grape). Total genome comparisons have been performed among these strains, indicating the presence of specific gene sets for each isolate (Bhattacharyya et al. 2002a,b; Costa de Oliveira et al. 2002). However, no attempt has been made to gather information about specific gene expression modulation under controlled experimental conditions. This problem becomes more relevant since the annotated genomes point to the existence of a high percentage of open reading frames (ORFs) whose functions are still unknown (over 50% of the genome from strain 9a5c), whereas genes assigned with putative functions may be associated with several, nonexclusive, possible pathogenicity mechanisms.

In this report, a spotted microarray carrying approximately 2200 ORFs from the genome of Xf 9a5c strain (Simpson et al. 2000) was hybridized with DNA from 11 Xf strains, obtained from different plant hosts and geographic origins. The results revealed the existence of a highly conserved core gene pool, containing all genes involved in the major biochemical pathways and cellular functions, which are likely to determine the adaptation of Xf to life in the plant xylem. In contrast, we found evidence that Xf displays a large flexible gene pool, characterized by the presence of several clusters containing foreign genes. Some of these clusters fulfilled the designation of genomic islands (GI), whereas others were considered putative GIs (Hacker and Kaper 2000; Hacker and Carniel 2001). The major GIs were consistently found in the citrus and coffee strains, but not in those obtained from other hosts. Moreover, transcriptome analysis revealed coordinated transcription of most genes within GIs and prophages, indicating that gene expression in these horizontally transferred elements can be controlled by growth conditions and, perhaps, other environmental stimuli. The number of genes contained in GIs, prophages, and plasmids indicates that Xf carries one of the largest flexible gene pools described to date, and dissection of its composition helps to explain the efficient adaptability of the group to infect such a wide range of plant species.

	RESULTS

Top ABSTRACT RESULTS DISCUSSION METHODS WEB SITE REFERENCES REFERENCES

 
Genomic Comparison of Xf Strains
Evidence for the existence of laterally transferred genes in Xf was first reported by Simpson et al. (2000), who mapped the presence of two ORFs, originally found in Dichelobacter nodosus and Ralstonia eutropha, within cryptic prophages found in the genome of Xf strain 9a5c. Recent genome comparisons involving Xf strains obtained from citrus, almond, and oleander plants broadened this scenario, suggesting that a series of conjugation-related ORFs found in Xf strains could have been acquired from soil-inhabiting bacteria (Bhattacharyya et al. 2002a). Moreover, a large (67-kb) element present in the genome of Xf strain 9a5c was identified as a new prophage, probably derived from the Siphophage group of double-stranded bacteriophages (Bhattacharyya et al. 2002b).

However, these comparisons employed only in silico approaches and therefore, relied solely on information derived from completely or partially sequenced genomes, limiting the scope of the observations to a small number of strains. In the present work, a microarray hybridization approach allowed comparison of the genome from strain 9a5c with a significantly larger set of Xf strains, obtained from different hosts: two strains from citrus, plum, and grapevine and one strain from almond, coffee, elm, mulberry, and ragweed (Table 1). The results obtained allowed detection of missing and duplicated genes across the 12 Xf strains, generating the genomic profiles shown in Figure 1B (all genetic differences observed in the analyzed strains can be found in the supplemental material). Although the precise genomic composition seems to be unique for each strain, one can clearly visualize an overall pattern that differentiates the citrus/coffee group (group 1) from the isolates obtained from plum, ragweed, grape, mulberry, elm, and almond (group 2).

View this table: [in this window] [in a new window]   Table 1. Xylella fastidiosaStrains Analyzed

View larger version (51K): [in this window] [in a new window]   Figure 1. A comparative and functional analysis of X. fastidiosa genomes. (A) The relative chromosomal position of genes that showed transcription upregulation in PW (red bars) and XDM (green bars) media. Most PW-overexpressed genes are located in clusters that coincide with prophages or GI regions. (B) The genomic profiles among 12 Xf isolates obtained from different hosts and disparate geographical origins. DNA from these isolates was labeled with Cy5-dCTP and submitted to competitive microarray hybridization with Cy3-labeled DNA from Xf strain 9a5c. ORFs were categorized as common to the two strains, exclusive to Xf strain 9a5c, or duplicated in the test strain, depending on the specific Cy5/Cy3 ratio. The figure shows the presence/absence of 2200 ORFs from Xf strain 9a5c across the genome of all tested strains. A red bar indicates the presence of a gene, a black bar, absence of the gene, and duplicated sequences are represented by a green bar. (C,D) The CG content and codon bias variations in strain 9a5c, respectively. The figure also indicates the relative positions of five genomic islands (GI1–GI5), six putative GIs (pGI1–pGI6), four integrated prophages (XfP1–XfP4), one plasmid integration region (pI1), and pXf51 plasmid, identified in the genome of Xf strain 9a5c.

View this table: [in this window] [in a new window]   Table 2. General Characteristics of Putative Horizontally Transferred Elements Found in the Genome of Xf Strain 9a5c

A second large GI (named GI₂) spans almost 67 kb, from ORF XF1718 to XF1792. It also displays typical GI characteristics, and ORF1718 is very similar to an integrase associated with a self-transmissible 105-kb element that carries the clcRABDE genes, which encode catechol-degrading enzymes in Pseudomonas putida (Ravatn et al. 1998). In Xf, this element is integrated at the 3' end of tRNA-G and carries a series of 12 ORFs encoding different ketoreductases/dehydrogenases and eight transcription regulators. A few other ORFs have assigned functions associated with outer membrane structure (including one fimbrillin gene), DNA replication, phage-related metabolism, and others. However, the majority of ORFs found in this element also encode hypothetical and conserved hypothetical proteins (42 ORFs). As mentioned above, an independent study by Bhattacharyya et al. (2002b) identified this element as another cryptic prophage. However, we see no strong evidence to support this hypothesis, because two independent genome annotations were not capable of identifying a complete set of genes encoding prohead components, phage regulatory factors, and a full DNA replication machinery, as expected from a typical integrated prophage (Simpson et al. 2000; Battacharyya et al. 2002a,b). We shall, therefore, continue to refer to this element as a GI.

Three smaller GIs were identified, spanning from ORFs XF0631 to XF0641 (GI₃), XF1859 to XF1885 (GI₄), and XF2108 to XF2132 (GI₅). Due to their reduced size (9, 15, and 16 kb, respectively), they should be referred to as genomic islets (Hacker and Carniel 2001). GI₃ is the only GI identified that does not display significant variation in codon bias and GC content, but other than that, GI₃ has all expected properties of a typical GI: It seems to be integrated at the 3' end of tRNA-S and carries a phage-related integrase (ORF XF0631), along with nine hypothetical proteins and one DNA modifying enzyme. It is conserved across group 1 and in a few group 2 isolates (elm, plum I, and plum II), but missing in the remaining strains. The small variations in codon bias and GC content and the presence of GI₃ in genomes other than citrus and coffee suggest that this islet might be more ancient than the others. GI₄ is inserted next to tRNA-K and carries only 27 ORFs, five with a known function and all encoding phage-related proteins, probably derived from bacteriophage Pf3. This islet is likely missing in all group 2 isolates, with the exception of the mulberry strain, where most of its structure seems to be preserved. Similarly, this GI is conserved in the group 1 strains, although some degree of genomic instability could be found in the coffee strain and both citrus isolates. Surprisingly, no known integrase was found within GI₄, and it is the only identified GI that does not display transcription upregulation when cells are grown in Periwinkle Wilt (PW) medium (see Methods). GI₅, on the other hand, bears a phage-related integrase (ORF XF2132) and is inserted next to tRNA-C. It also carries a great number of ORFs with unknown function, and another copy of the D. nodosus vapE gene is present in its structure. This element displays a partial inverted duplication of the vapE-containing region of GI₁ and seems to be conserved in the coffee strain, whereas it is duplicated in the two citrus isolates and partly missing in all group 2 representatives.

A series of other regions with altered codon bias and GC content were encountered throughout the genome of strain 9a5c, but because most of them did not seem to be associated with areas of genomic instability in the tested isolates, we decided to name them pGIs (for putative GIs). The general characteristics of these pGIs are listed in Table 2. They all resemble, to some extent, typical GIs, but the most striking observation is that most of them also carry PW-upregulated genes, either along their entire structure or associated with subsets of genes (see below).

Finally, a large 52-kb element located between ORFs XF2025 and XF2087 seems to be the result of plasmid-mediated gene integration. Although the majority of genes from this region are missing in the group 2 isolates, no alterations in GC content or codon bias were observed. Moreover, this element does not carry any integrase-like ORF and is not inserted in the vicinity of tRNA-like genes. A series of ORFs encoding bacterial conjugation factors (TraCDEF and TrbBCDEFGHIJLN) are present in this element. These tra/trb genes are typically found in widespread conjugationally transmitted plasmids, originally isolated from E. coli and other enterobacteria (Motallebi-Veshareh et al. 1992). This element also bears two plasmid stabilization proteins, ParD/E, reinforcing its possible plasmid-mediated origin. This element was named PI₁ (for Plasmid Integration) and also carries a particularly noteworthy multidrug efflux system (ORFs XF2083 to XF2085), along with three transcription regulators (ORFs XF2038, XF2062, and XF2086) and a series of 34 hypothetical/conserved hypothetical protein-coding ORFs.

Transcriptome Analysis of Xf Strain 9a5c Grown in Different Media
Lateral gene transfer is currently viewed as an important mediator of bacterial evolution. However, acquisition of foreign genes present in GIs does not necessarily lead to the development of new phenotypes, because these newly acquired elements will only be expressed in the host cell if recognized by the transcriptional machinery and their expression control is efficiently coordinated with the rest of the genome (Hacker and Carniel 2001). Thus, the presence of so many horizontally transferred elements in Xf strain 9a5c raised questions about the functional role(s) played by their resident genes, because there is no evidence for their effective transcription in vivo. Moreover, because most GIs are lacking in some of the strains under analysis, it is possible that these elements are nonessential for bacterial survival.

Thus, to verify whether Xf can efficiently regulate the expression of genes present in any of the transposable elements mentioned above, microarray hybridizations were carried out with RNA samples obtained from bacteria grown in two very distinct media: Periwinkle Wilt (PW; Davis et al. 1981) and Xylella Defined Medium (XDM; Lemos and Alves 2001). PW contains amino acids, BSA, and peptone as sources of carbon, nitrogen, and other nutrients, whereas XDM contains only glycerol and glutamic acid. Due to the marked differences in the composition of these two media, they are able to induce variations in bacterial physiology and growth kinetics. Indeed, when cells were harvested for RNA extraction after several days of incubation, Xf grown in PW showed a twofold higher cellular density compared to cells grown in XDM, even though both cultures had reached early to mid-exponential growth phase at this point, as verified by OD₆₀₀ measurements (data not shown).

As expected, the two cultures revealed striking differences in their respective transcriptomes (see supplemental materials for a complete list of modulated genes). Surprisingly, only 55 genes were found to be overexpressed in XDM, and their transcription modulation seems to be much smaller compared to what is observed in PW-overexpressed genes (Fig. 2). Among these, 17 ORFs encode hypothetical and conserved hypothetical proteins. ORFs encoding stress-related factors, such as GroEL, GroES (ORFs XF0615 and XF0616), and SodA (ORF XF1921) were also found under these conditions, as were several ribosomal proteins. These findings are only partly consistent with the observations made by Paustian et al. (2002), who used microarray hybridizations to evaluate the transcriptional response of Pasteurella multocida to nutrient limitation, comparing bacteria grown in rich and minimal media. The Paustian studies showed that one-third of the bacterial genome appeared to have its transcription modulated over the course of the experiment. A greater number of genes, including those involved in energy metabolism, transport, protein synthesis, and binding were overexpressed in rich medium, whereas growth in minimal medium led to upregulation of a smaller number of genes, encoding enzymes associated with amino acid biosynthesis, membrane transport systems, and heat shock proteins. Some of the differences observed between the Paustian study and the findings reported herein might be due to the composition of XDM, which, although limited in nutrient availability compared to PW, cannot be considered minimal medium, which is yet to be developed for X. fastidiosa. Furthermore, it is important to note that most genes involved in major biosynthetic pathways were not modulated in any of the growth conditions. Although a few such genes were found to be upregulated in PW (see below), their expression could also be detected after hybridization with RNA obtained from cells grown in XDM, suggesting that X. fastidiosa might have most of these ORFs under the control of poorly induced and/or constitutive promoters, which may represent an adaptation to life in a naturally nutrient-poor environment such as the plant xylem.

View larger version (54K): [in this window] [in a new window]   Figure 2. Results from competitive hybridizations between RNA samples obtained from cells grown in XDM medium (labeled with Cy3) and PW medium (labeled with Cy5). Each spot on the array is represented by a dot in the scatter plot shown above. The graph displays the specific, background-subtracted, average intensities for each gene on the array. Gray lines represent fold change values comparing the Cy3 and Cy5 intensities. Lines below the diagonal represent XDM/PW values of two-, three-, five- and 10-fold, respectively. Lines above the diagonal represent XDM/PW ratios of 1/2, 1/3, 1/5, and 1/10.

Surprisingly however, chromosomal mapping of PW-overexpressed genes showed that most of them were present within prophages, GIs, and pGIs detected in Xf strain 9a5c (Fig. 1A). Moreover, regulation of gene expression within these elements seems to occur in a coordinated manner, because large sets of contiguous genes were simultaneously upregulated. For example, in the large, 67-kb GI₂, upregulation of nearly all genes could be observed. This coordinated control of transcription was also observed in the other GIs (with the exception of GI₄, which bears no integrase gene), as well as in prophages XfP1, XfP2, and XfP3 (Fig. 1A). Although there is no evidence of phage particle release by Xf cells, the coordinated expression pattern observed across XfP1, XfP2, and XfP3 indicates that these elements might have conserved a relatively independent mechanism of transcription control. It was surprising, however, to see that the same phenomenon was associated with genes present in a genomic island, because there is no precedent for the existence of a similar control mechanism in GIs.

	DISCUSSION

Top ABSTRACT RESULTS DISCUSSION METHODS WEB SITE REFERENCES REFERENCES

 
The recent advances in the field of comparative genomics resulted in the development of a new paradigm to explain bacterial evolution. It is now clear that prokaryotic genomes are composed of a core gene pool and a flexible gene pool (Hacker and Kaper 2000; Hacker and Carniel 2001). The majority of genes in the core gene pool encode proteins that play roles in basic cellular functions and exhibit a rather homogeneous GC content and codon usage preference. Elements from the flexible gene pool, on the other hand, often encode additional functions that are not essential for bacterial growth, but provide advantages under particular conditions, such as changes in the environment and adaptation to a new host. These genes are normally associated with phages, plasmids, integrons, and other transposable elements. Their contribution to the overall gene content has been estimated to vary among bacteria (Hacker and Carniel 2001). Some microorganisms, such as E. coli and Synechocystis sp., have larger flexible gene pools, corresponding to 18% and 16% of the total genome, respectively (Lawrence and Ochman 1998; Ochman et al. 2000). This percentage is around 10% in Methanobacterium thermoautrophicum and Mycoplasma pneumoniae, 5% in Helicobacter pylori and Deinococus radiodurans, and virtually 0% in Rickettsia prowasekii and Mycobacterium genitalium (Ochman et al. 2000). In Xf 9a5c, the number of genes carried by plasmid pXf51, prophages, pI₁, and GIs add up to 510, which corresponds to 18% of the total number of ORFs annotated by Simpson et al. (2000). Thus, the flexible gene pool in this bacterium seems to be one of the largest described to date. Moreover, different Xf strains display variable genome sizes and carry strain-specific gene sets, which mainly include ORFs encoding conjugal transfer factors, restriction-modification systems, methyltransferases, hypothetical/conserved hypothetical proteins, and other genes typically found in horizontally transferred elements (Bhattacharyya et al. 2002a,b). Thus, strain-specific flexible gene pools are likely to be found within the Xf group, a theory that is reinforced by the presence of different plasmids and a putative plasmid-derived Pathogenicity Island in some North American isolates (Pooler et al. 1997; Bhattacharyya et al. 2002a,b).

Several examples of GIs have been observed in different microorganisms associated with symbiotic and saprophytic behaviors, as well as other aspects of environmental adaptation and pathogenicity (Sullivan and Ronson 1998; Hacker and Kaper 2000; Hacker and Carniel 2001). Efforts have been made to understand the control of gene expression in GIs, which are believed to be coordinately regulated with the core genetic pool, allowing fully functional integration of GIs and other components of the flexible gene pool to the global regulatory network of the cell. For instance, the global regulator PhoP-PhoQ of Salmonella typhimurium is located in the ancestral chromosome (Deiwick et al. 1999), but can positively regulate expression of the mgtC genes located in a GI. Other examples are the ToxT and YbtA transcriptional activators of V. cholerae and Yersinia sp., respectively, which are located in GIs, but regulate the expression of genes located outside these regions (Carniel 2002; Sarkar et al. 2002). However, these analyses always addressed the study of gene expression in PAIs and other GIs through the individual evaluation of a few pathogenicity genes. Our array-based global analysis of gene expression in Xf demonstrated a situation where the entire length of several elements of the flexible gene pool have been subjected to upregulation. This phenomenon seems consistent with relatively recent acquisition of such horizontally transferred elements by Xf, with insufficient evolutionary time to guarantee proper gene expression coordination with the core genetic pool.

Interestingly, this upregulation does not seem to be related to differences in cell density or growth kinetics, because hybridizations comparing RNA samples obtained from cells harvested at early and late exponential growth phases did not show distinct expression profiles (data not shown). Thus, it is possible that these laterally acquired genes may be responding collectively to environmental stimuli, perhaps due to differences in media composition. This raises the possibility that GIs might not only be capable of transferring large gene sets among bacteria, but may also control their transcriptional activation in a rather independent manner. If that is the case, the ability of GIs to mediate the acquisition of new characteristics might have a much deeper impact on bacterial evolution than previously imagined, and a better assessment of the mechanisms that control gene expression in transposable elements might be of extreme importance to understand the processes that mediate pathogenicity and environmental adaptation of microorganisms.

In the case of Xf, these genes are likely to hold the key to explain the broad host range of this phytopathogen, because the overall genomic diversity observed among the strains is mainly due to elements from the flexible gene pool. As a result, most variations involve ORFs associated with phage metabolism, plasmid maintenance, bacterial conjugation, and enzymes associated with environmental response and adaptation, such as transcriptional regulators, drug-resistance enzymes, and virulence factors. Cell wall- and outer membrane-related structural components also accounted for significant differences among the tested strains. Particularly noteworthy is the observation of differences among fimbrillin genes, whose products may establish direct contacts with host cells and other bacteria (Park et al. 2001). One fimbrillin (ORF XF0539) was not found in the strains obtained from elm and plum, whereas the almond, grape, mulberry, and ragweed isolates lacked four of the fimbrillin genes from strain 9a5c (ORFs XF0487, XF0538, XF0539, and XF1791). These factors are adhesins of Type IV fimbriae found in numerous bacterial species that infect humans, animals, and plants and are important players in the process of host colonization and disease development (Kennan et al. 2001; van Doorn 2001). Another enzyme involved in modification of outer membrane carbohydrates (OafA–ORF XF0778) is missing from all North American strains except grape-FL. Virulence factors, such as bacteriocin, VapD, VapE and the virulence regulator XrvA (ORFs XF2007, XFa0052, XF0506, and XF0749, respectively) also displayed similar profiles, as did ORF XF1744, a hypothetical oxidoreductase, possibly associated with invasive capability, which is also missing in all group 2 strains (Gaillard et al. 1991). Although such differences among pathogenicity factors might help to explain the disparate pathogenies caused by Xf strains, recent evidence points to host-specific factors as key mediators of disease manifestation by these bacteria (Li et al. 2002).

The subdivision of Xf strains into at least two distinct groups is in accord with previous observations (Chen et al. 2000; Hendson et al. 2001; Mehta and Rosato 2001), but even inside each group, a considerable degree of variation can be observed. The citrus group bears duplicated regions not observed in the coffee strain. The latter isolate does not carry an 32-kb-long region of pI₁—possibly excised by homologous recombination between two copies of the traC gene (ORFs XF2025 and XF2061). Within group 2, the almond/grape strains display a similar profile, whereas elm and plum II isolates are nearly identical to each other, contrasting with plum I. Independent studies found strain plum I to have a closer phylogenetic relationship with the citrus-coffee strains than with others (da Costa et al. 2000; Mehta and Rosato 2001). This finding is in accord with geographic distribution, as plum I was originally obtained from a diseased tree in the southern region of Brazil, being the first Xf isolate obtained in that country (French and Kitajima 1978). Interestingly, strain plum I was the only group 2 isolate to bear nearly all ORFs from pXf51, the large 51-kb plasmid found in CVC strain 9a5c. However, considering the overall genomic profile, it is likely that strain plum I derives directly from the North American isolates from group 2. It is tempting to speculate that this isolate might have been brought to Brazil as a graft or seedling contaminant and acquired pXF51 sequences through lateral gene transfer after coming into contact with other South American Xf strains, raising interesting possibilities concerning genetic interactions among Xf isolates in the wild and genomic plasticity within the group. In this sense, the observation that GI₂ seems to be present only in group 1 suggests that it may have been acquired after geographical separation between North and South American strains. The remaining horizontally transferred elements, however, appear to display a more complex distribution pattern, hindering more conclusive evaluation of their evolutionary history. For instance, the mulberry strain seems to carry most sequences from GI₄, which is lacking in other group 2 representatives. On the other hand, this very same strain does not carry sequences from GI₃, which is present in elm, plum I, and plum II. Several sequences from prophage XfP1 are lacking in elm and plum II strains, but present in other group 2 isolates. These observations reinforce the idea that the flexible gene pool in Xf might display a complex, strain-specific evolutionary pattern.

Nonetheless, despite all of the differences categorized and discussed above, it must be noted that the overall similarity among Xf genomes was very high, with more than 90% ORF conservation between strain 9a5c and the group 2 representatives. The conservation of ORFs increased to 97% when strain 9a5c was compared to the coffee strain and to 99% when the comparison involved the other citrus isolates. Essentially, all genes involved in the major biosynthetic pathways and cellular functions are present in all analyzed Xf strains, suggesting that adaptation of the entire group to life in the plant xylem is likely to depend on genes from the highly conserved core gene pool, responsible for the major physiological characteristics already discussed by Simpson et al. (2000), which involve a mechanism for energy generation based on the utilization of carbohydrates, including enzymes for cellulose degradation, but with no alternative energy sources based on catabolism of fatty acids or amino acids. Fully equipped biochemical pathways guarantee biosynthesis of all amino acids, purines, pyrimidines, and other nutritional requirements, which are also met through the presence of an elaborated panel of membrane transporters for the efficient uptake of different ions and organic compounds. On the other hand, up to 18% of the genomic composition of strain 9a5c may be associated with elements from the flexible gene pool, which have been shown to carry important adaptation and putative pathogenicity factors. Moreover, we provide evidence that gene expression within these elements can be efficiently modulated by environmental stimuli, and that Xf strains are likely to exchange genetic material in the wild, contributing to the considerable genomic differences observed among strains from groups 1 and 2. Thus, an efficient mechanism for exchanging mobile elements in a large flexible gene pool may be of capital importance in mediating adaptation of Xf to new hosts.

	METHODS

Top ABSTRACT RESULTS DISCUSSION METHODS WEB SITE REFERENCES REFERENCES

 
Origin of the Isolates and Growth Conditions
Xf strains from different hosts and geographic origins (Table 1) were grown in either Periwinkle Wilt (PW) medium (Davis et al. 1981), which contained 4.0 g phytone peptone, 1.0 g tryptcase peptone, 1.2 g K₂HPO_4, 1.0 g KH₂PO_4, 0.4 g MgSO₄, 0.1% hemin chloride, and 10 mL phenol red per L, or in Xylella Defined Medium (XDM), which contained 0.23 g K₂HPO_4, 0.10 g MgSO₄.7H₂O, 1.1 g glutamic acid, 4 mL glycerol, and 0.5 g ferric pyrophosphate per L (Lemos and Alves 2001). The bacteria were initially grown on solid medium and then transferred consecutively to increasing volumes of liquid medium (5, 20, 50, and 100 mL). In each step, the cells were incubated for 4 d in a rotatory shaker at 28°C. Because the bacteria grew adhering to the flask walls, the cellular mass was disrupted by pipetting and vortexing before proceeding to the next transfer. After the last step, the cellular mass was transferred to 200 mL of PW or XDM liquid medium and incubated for another 4 d. Bacterial growth in both PW and XDM media was evaluated through measurement of optical density (OD₆₀₀). Aliquots (1.0 mL) were sampled at 24-h intervals over a 10-d period. The samples were centrifuged, and the resulting pellet was disrupted in 1.0 mL of saline (0.9% NaCl) and OD₆₀₀ immediately read in a spectrophotometer. Three independent samples were analyzed for each time point.

Microarray Fabrication
Specific primer pairs for amplification of 2838 ORFs found in the genome of Xf strain 9a5c were designed employing the software Primer3 (kindly provided by Dr. S. Rozen at http://www.genome.wi.mit.edu/genome_software/other/primer3.html). Approximately 2500 primer pairs, representing 2205 different ORFs from the main chromosomal DNA and from the pXf51 plasmid, yielded specific PCR amplification products. Most of the amplified fragments ranged in size between 350 and 1000 base pairs. The PCR products obtained were purified with the QiaQuick 96 PCR purification system (QIAGEN) and spotted onto CMT-GAPS silane-coated slides (Corning), using an Affymetrix 417 Arrayer, according to the manufacturer's instructions. The printed microarrays were used for competitive hybridization experiments using DNA samples obtained from Xf strains. This approach allowed a global comparison of approximately 2200 ORFs from Xf strain 9a5c with the genome of the other 11 Xf isolates. Although the strain 9a5c genome is composed of more than 2800 ORFs, many of them were, actually, paralogous duplications scattered throughout the chromosome and/or small ORFs (<250 bp) encoding hypothetical proteins. Many of these sequences could not be efficiently and/or specifically amplified, due to their reduced size and chromosomal location (closely flanked by larger ORFs). We estimate, however, that these analyses cover more than 90% of the most relevant ORFs present in the genome from Xf strain 9a5c. The microarrays were also used to monitor the global transcription profile of Xf strain 9a5c under two different growth conditions.

Nucleic Acid Labeling and Hybridization Conditions
DNA samples were extracted according to the protocol described by Dungan (1989). Labeling reactions and purification were performed as described by Behr et al. (1999). Total RNA was isolated from Xf cells using the RNeasy RNA extraction kit (QIAGEN) and labeled by reverse transcription. Briefly, 30 µg of total bacterial RNA was mixed with 16 µg of random hexamers (Invitrogen) in a 30 µL final volume. Annealing was accomplished by incubation for 2 min at 75°C, 2 min at 55°C, 2 min at 45°C, 2 min at 37°C, and 2 min at 22°C, followed by addition of 6 µL SuperScript II reaction buffer (Invitrogen), 3 µL of 0.1 M dithiothreitol, 1 µL dNTP mix (10 mM dATP, 10 mM dGTP, 10 mM dTTP, 5 mM dCTP), 2 µL of 1.0mM Cy3- or Cy5-labeled dCTP (Amersham Biosciences), and 2 µL of SuperScript II reverse transcriptase (200 U/µL; Invitrogen) into the reaction. The cDNA synthesis was carried out at 42°C for 2 h. After labeling, the RNA was hydrolyzed with 0.1 N NaOH treatment and neutralized by the addition of 0.1 N HCl. The labeled cDNA was diluted to 500 µL with TE and concentrated using a Microcon-30 (Amicon) to 10 µL.

Arrays were hybridized overnight (42°C) in a GeneTac Hybridization Station (Genomic Solutions), in 6 x SSC, 5 x Denhardt's solution, 0.25 mg/mL sheared salmon sperm DNA, 0.5% SDS, and 2 µg of each labeled DNA or cDNA sample. After hybridization, slides were washed twice (42°C) in 0.5 x SSC, 0.01% SDS, followed by two washes in 0.06 x SSC, 0.01% SDS and two final washes in 0.06 x SSC. All washing steps consisted of 1 min of flow, followed by 5 min of incubation. Slides were then dried and submitted to fluorescence detection.

Image Acquisition and Analysis
Hybridized arrays were scanned in an Affymetrix 418 Array Scanner, and images were analyzed with Affymetrix Jaguar v 2.0. Quality control of the hybridized spots is automatically performed by the software, based on spot morphology and local signal-to-background ratio, using the Easy Threshold and Variable Circle Size Algorithm. In all experiments, reliable hybridization signals could be obtained for more than 90% of the arrayed probes (see Costa de Oliveira et al. 2002). For DNA–DNA hybridizations, normalization between the intensities in the two channels was achieved with the Jaguar Control Spots option, using a list of 30 control ORFs that had 100% sequence identity in the genomes of strains 9a5c and Grape CA. For each pair of strains, hybridizations were performed in triplicate. These data were consolidated into a GATC database with Affymetrix MicroDB v 2.0, and the averages from all six readings were submitted to scatter plot visualization with Affymetrix Data Mining Tool v 2.0. Statistical validation of fold change variations was performed with the aid of the Significance Analysis of Microarrays (SAM) method proposed by Tusher et al. (2001), using the software developed by B. Narasimham and available at http://www-stat.stanford.edu/tibs/SAM/index.html. Spots that showed a Reference/Test ratio <1:2 were considered to be present in greater copy number in the test over the reference strain, as proposed by Smoot et al. (2002), whereas spots that showed an average Reference/Test ratio >5:1 were considered to be missing in the test strain. The application of these criteria in a direct sequence comparison between strains 9a5c and Grape CA, which have been completely sequenced, provided an estimated error rate below 0.3% (Costa de Oliveira et al. 2002).

For cDNA–cDNA hybridizations, three independent experiments were performed with different RNA preparations. Hybridizations were performed in triplicate, generating nine intensity readings for each gene. Normalization was achieved through the sum of all intensities in each channel, using the Jaguar All Spots option. Normalized data were consolidated into a GATC database with Affymetrix MicroDB v 2.0, and genes that displayed statistically significant variations in gene expression were identified with the aid of the SAM method. A fold change cutoff of at least 2.0 was used in these analyses.

Genome comparisons and transcription profiles were viewed using the programs Cluster and TreeView, available at http://www.microarrays.org/software (Eisen et al. 1998). For the visualization of comparative profiles from the 12 analyzed strains, we applied the method proposed by Smoot et al. (2002), where ORFs shared by the reference and each test strain was labeled 1, whereas ORFs exclusive to the reference strain, or present in greater copy number in the test over the reference strain were labeled 0 and –1, respectively.

	WEB SITE REFERENCES

Top ABSTRACT RESULTS DISCUSSION METHODS WEB SITE REFERENCES REFERENCES

 
http://www.jgi.doe.gov/jgi_microbial/html; Accession to the partial genome sequences from Xf strains isolated from Almond and Oleander trees.

http://www.lbi.ic.unicamp.br/world/xf-grape; Accession to the complete genome sequence from Xf strains isolated from grapevines.

http://www-genome.wi.mit.edu/genome_software/other/primer3.html; Accession to software PRIMER3.

http://www-stat.stanford.edu/tibs/SAM/index.html; Accession to software SAM.

http://www.microarrays.org/software; Accession to softwares CLUSTER and TREEVIEW.

	Acknowledgements

 
This work was supported by Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP), under the X. fastidiosa Functional Genome Program.

The publication costs of this article were defrayed in part by payment of page charges. This article must therefore be hereby marked "advertisement" in accordance with 18 USC section 1734 solely to indicate this fact.

	Footnotes

 
⁶ Corresponding authors.

E-MAIL Lnunes{at}umc.br; FAX 55-11-4798-7104.

E-MAIL reginaco{at}umc.br; FAX 55-11-4798-7106.

Article and publication are at http://www.genome.org/cgi/doi/10.1101/gr.930803.

	REFERENCES

Top ABSTRACT RESULTS DISCUSSION METHODS WEB SITE REFERENCES REFERENCES