Microbial Biotechnology
Electronic Journal of Biotechnology ISSN: 0717-3458 Electron. J. Biotechnol.
© 2012 by Pontificia Universidad Católica de Valparaíso -- Chile Vol. 15 No. 3, Issue of May 15, 2012
DOI: 10.2225/vol15-issue3-fulltext-4 Received November 28, 2011 / Accepted March 28, 2012
RESEARCH ARTICLE

Xylella fastidiosa: An in vivo system to study possible survival strategies within citrus xylem vessels based on global gene expression analysis

María Teresa Federici*1,2 · Jackson A. Marcondes2 · Simone Cristina Picchi2,6 · Eduardo S. Stuchi3 · André Luiz Fadel4 · Marcelo L. Laia5 · Manoel Victor F. Lemos7 · Eliana Gertrudes Macedo Lemos2

1 Instituto Nacional de Investigación Agropecuaria (INIA), Estación Experimental Wilson Ferreira Aldunate, Unidad de Biotecnología, Rincón del Colorado, Uruguay
2 Universidade Estadual Paulista “Júlio de Mesquita Filho”, Faculdade de Ciências Agrárias e Veterinárias, Departamento de Tecnologia, Laboratório de Bioquímica de Microrganismos e Plantas, Jaboticabal, Brazil
3 Embrapa Mandioca e Fruticultura, Núcleo Tecnológico de Recursos Genéticos e Desenvolvimento de Variedades - Nugene, Cruz das Almas, BA, Brasil and “Estação Experimental de Citricultura”, Bebedouro, São Paulo, Brazil
4 Universidade de São Paulo, Escola Superior de Agricultura "Luiz de Queiroz", Piracicaba, São Paulo, Brazil
5 Universidade Federal dos Vales do Jequitinhonha e Mucuri, Faculdade de Ciências Agrárias, Departamento de Engenharia Florestal, Diamantina, MG, Brasil
6 Centro de Citricultura Sylvio Moreira - APTA, Cordeirópolis, SP, Brasil
7 Universidade Estadual Paulista “Júlio de Mesquita Filho”, Faculdade de Ciências Agrárias e Veterinárias. Departamento de Biologia Aplicada à Agropecuária, Jaboticabal, Brazil

*Corresponding author: mfederici@lb.inia.org.uy

Financial support: This research was supported by the Ministry of Culture of Uruguay through the Program for Technical Development (PDT), the Agency CAPES (“Coordinadoria de Aperfeiçõamento de Pessoal de Nível Superior”) and the National Institute of Agronomic Research of Uruguay (INIA).

Keywords: citrus, microarray, Xylella fastidiosa

Abstract    

Xylella fastidiosa inhabits the plant xylem, a nutrient-poor environment, so that mechanisms to sense and respond to adverse environmental conditions are extremely important for bacterial survival in the plant host. Although the complete genome sequences of different Xylella strains have been determined, little is known about stress responses and gene regulation in these organisms. In this work, a DNA microarray was constructed containing 2,600 ORFs identified in the genome sequencing project of Xylella fastidiosa 9a5c strain, and used to check global gene expression differences in the bacteria when it is infecting a symptomatic and a tolerant citrus tree. Different patterns of expression were found in each variety, suggesting that bacteria are responding differentially according to each plant xylem environment. The global gene expression profile was determined and several genes related to bacterial survival in stressed conditions were found to be differentially expressed between varieties, suggesting the involvement of different strategies for adaptation to the environment. The expression pattern of some genes related to the heat shock response, toxin and detoxification processes, adaptation to atypical conditions, repair systems as well as some regulatory genes are discussed in this paper. DNA microarray proved to be a powerful technique for global transcriptome analyses. This is one of the first studies of Xylella fastidiosa gene expression in vivo which helped to increase insight into stress responses and possible bacterial survival mechanisms in the nutrient-poor environment of xylem vessels.

Introduction

Xylella fastidiosa is a gram-negative, xylem-limited bacterium, phylogenetically most closely related to Xanthomonas spp., and is transmitted from plant to plant by several xylem-feeding insect vectors (Hopkins, 1989). Strains of X. fastidiosa cause diseases in many economically important plants, including grapevines, citrus, peach, plum, oleander, elm, sycamore, oak, maple, and coffee (Purcel, 1997; De Lima et al. 1998). The major symptoms of most X. fastidiosa diseases are associated with water stress, due to reduced xylem flow, which is thought to result from occlusion of the xylem vessels by bacterial aggregates that likely contain extracellular polysaccharides (EPS) (Da Silva et al. 2001), gums, and tyloses (Hopkins, 1989). Diseases caused by Xylella fastidiosa have attained great importance worldwide as the pathogen and its insect vectors have been disseminated.

The genome sequence determination of the pathogenic strain 9a5c, associated with citrus variegated chlorosis (CVC) was reported by Simpson et al. (2000), revealing many genes associated with pathogenicity and virulence, involving toxins, antibiotics and ion sequestration systems, as well as bacterium-bacterium and bacterium-host interactions. Also, at least 83 genes were bacteriophage-derived and included virulence-associated genes from other bacteria of the genome, providing direct evidence of phage-mediated horizontal transfer (Simpson et al. 2000).

X. fastidiosa is a member of the gamma subdivision of Proteobacteria which inhabits the plant xylem, a nutrient-poor environment. Mechanisms to sense and respond to adverse environmental conditions are certainly extremely important for bacterial survival in the plant host (Koide et al. 2006). Survival can be regarded in a broader context as the fitness of an organism to quickly adapt to particular environmental changes, and to ensure adaptation, bacteria rely on invasive DNA (transposons, integrons, plasmids and prophages) brought in by horizontal gene transfer mechanisms to generate the genetic diversity needed. Such elements can carry genes that code for a specific ecological adaptation, for new features associated with pathogenicity or bacterium-host interaction (Sullivan and Ronson, 1998; Wong and Golding, 2003; Monteiro-Vittorello et al. 2005). Also, a wide range of signal molecules, environmental stresses, and other factors trigger changes in gene expression in human, animal and plant pathogens, resulting in new phenotypes that increase the fitness of individuals (Camilli and Bassler, 2006; Shild et al. 2007). Phenotypic changes are especially important for pathogens with complex life histories, which must survive in multiple environments under different habitat-dependent stresses. Therefore, many pathogens have phenotypic plasticity that permits the exploration of multiple environments, which is likely regulated by signals or cues provided by their habitat (Killiny and Almeida, 2009). Therefore, an in vivo transcriptomic analysis could provide many clues to the understanding of adaptation and survival of X. fastidiosa.

In this work we characterized the global gene expression profile of Xylella fastidiosa in vivo, without culturing bacteria, in order to achieve a deeper understanding of the stress responses and adaptability of this pathogen. We focused on the transcriptional activity of some genes related to biochemical processes involved in heat shock and antioxidant responses, toxins and detoxification processes as well as other possible mechanisms involved in bacterial survival. For this purpose, we used microarray analysis to identify global changes in gene expression patterns using a DNA microarray previously standardized by the Laboratory of Biochemistry of Microorganisms and Plants, by Travensolo et al. (2008) and Travensolo et al. (2009), including all genes identified in the X. fastidiosa genome.

We analyzed the differential expression of genes when the bacteria infected either a susceptible or a tolerant sweet orange plant. The susceptible variety used in this study, which rapidly developed severe symptoms, was Pera (Citrus sinensis (L.) Osbeck) and the tolerant was the cultivar Navelina ISA 315 (commonly known as “Bahianinha”). This cultivar belongs to the navel group of sweet oranges. This cultivar was originated in Italy in 1976 as the result of a clone that was recovered in vitro by immature ovule culture. It was introduced in Brazil and established in the field by the year 2000 to initiate CVC resistance studies. In preliminary studies, Navelina ISA 315 was found to be an asymptomatic host for Xylella fastidiosa (Stuchi et al. 2007), since other sweet oranges are all susceptible to CVC. The use of resistant or tolerant varieties is one of the main options for managing the disease. Fadel et al. (2010) confirmed that Navelina ISA 315 is a tolerant cultivar using visual observation of symptoms in greenhouse and field, as well as quantification by RT qPCR analysis.

Different expression patterns in symptomatic and asymptomatic plants would be expected, due to differences in plant defenses, xylem conditions or endophytic microorganism competition. This information can provide insights into bacterial survival mechanisms within susceptible and tolerant plants, as well as into plant-pathogen interactions.

Materials and Methods

Microarray construction

For genomic DNA extraction, X. fastidiosa 9a5c was cultivated in Petri dishes containing BCYE medium (Wells et al. 1981) at 28ºC for six days. Genomic DNA was extracted according to Ausubel et al. (1987) with a modification that includes RNAse treatment (200 µg/ml) for 1 hr 30 min at 37ºC. Specific primer pairs were designed using PRIMER3 software and used to amplify the 2,600 ORFs of the complete genome of X. fastidiosa. A melting temperature ranging from 48 to 57ºC was determined for all primers. PCR reactions were carried out in 1X PCR buffer (50 mM KCl, 20 mM Tris-HCl, pH 8.4), 2 mM MgCl2, 10 mM dNTP, 2 U Taq DNA polymerase, 5 pmols of each primer, 60 ng of genomic DNA and pure sterile water to a total volume of 10 µl. The following reaction conditions were used: 94ºC/2 min; 35 cycles of 94ºC/1 min, 58ºC/1 min, 72ºC/1 min 30 sec; followed by a final extension of 72ºC/5 min. All products were analyzed and quantified by electrophoresis in 1.5% agarose gel in 1X TBE buffer with 0.5 µg/ml ethidium bromide. Amplified products were suspended in 50% v/v DMSO at a final concentration of 10 to 30 ng/µl and arranged in duplicate at a distance of 25 µm on glass slides treated with aminosilane (Corning). Printing of microarrays was done by a robot model GMS 417 Arrayer (Affymetrix Inc.). After printing, the DNAs were re-hydrated (42ºC for 1 sec), dried (70ºC for 1 min) and fixed by UV cross-linking (130 x 10 µJ/cm2). The slides were kept at 70ºC for 2 hrs and then stored under vacuum at room temperature. Genetically distant negative controls were included in this array, consisting of human genes (pHUM1 and pHUM7) and plant genes (707050B11-Rubisco).

Experimental design and controls

From each host, three individual plants were collected. RNA was extracted, and fluorescent labeled cDNA was prepared and used to hybridize the slides. The slides contain PCR products from the 2,600 Xylella fastidiosa ORFs spotted in duplicate as well as the genetically distant controls. Xylella fastidiosa within asymptomatic cultivar Navelina ISA 315 was considered in this analysis as the experimental condition (labeled with Cy5 fluorophore) whereas susceptible variety Pera was considered the control condition (labeled with Cy3 fluorophore). Three control hybridizations of healthy samples of the cultivar Navelina ISA 315 and healthy Pera variety were performed in order to detect non-specific hybridization with the plant genome. Both healthy varieties produced by micrografting were collected within the greenhouse which was protected from any contamination by insects.

Plant material used for total RNA extraction

Plants used in this study were collected from the Citrus Experimental Station of Bebedouro (SP-Brazil), including symptomatic leaves of the susceptible variety Pera and asymptomatic leaves from the cultivar Navelina ISA 315. Samples from three different trees of each variety were collected in order to minimize variability due to individual plants. Navelina ISA 315 plants had been grafted over symptomatic Pera to ensure a high pressure of Xylella fastidosa inoculum. After 10 months, bacterial concentration was evaluated by observation of symptoms in the field based on a diagrammatic scale developed by Amorim et al. (1993), and by real time PCR (RT qPCR) in two periods, January 2009 and August 2009. These plants were found to be asymptomatic in the field, however, RT qPCR revealed high concentrations of bacteria at the first evaluation, in January 2009. Bacterial concentration in these plants did not decrease between evaluations as revealed by RT qPCR, nor were symptoms observed in the field in August 2009 (Fadel et al. 2010). RT qPCR analyses were performed based on the standard curve for Xylella fastidiosa developed by Oliveira et al. (2002). Samples of symptomatic Pera were collected from an orchard highly affected with CVC disease, located at about 100 meters from the Navelina ISA 315 trees, and from the greenhouse in the case of the healthy plant controls.

Isolation of total RNA

As Xylella fastidiosa is a xylem limited pathogen, the central veins (midribs) were cut from about 33 to 35 leaves of each plant and were macerated with liquid nitrogen. An RNA extraction methodology was used which involved a monophasic solution of phenol and guanidine isothiocyanate (TRizol®, Invitrogen). RNA integrity was verified in a 1.2% (w/v) agarose gel in 1X buffer (20 mM MOPS, 5 mM sodium acetate; 1 mM EDTA), DEPC-treated H2O and 6.7% w/v formaldehyde. The RNA molecular marker used was “RNA ladder” (Life Technologies) with six fragments visualized from 0.24 to 9.49 kb. The same procedure was performed for samples from both symptomatic and asymptomatic plants, as well as healthy control plants.

Synthesis of fluorescent cDNA from total RNA

Synthesis of fluorescent cDNA was carried out using 100 µg of total RNA, 2 µg of random primers (GE Healthcare, Piscataway, NJ, USA), and 4 U of RNAsin (Invitrogen). This mixture was incubated at 70ºC for 5 min and cooled to 4ºC for 5 min to allow the primers and mRNA template to anneal. Then, a solution containing 6 µl of ImProm 5X reaction buffer; 3.6 µl of MgCl2 25 mM, 3 µl of dNTPs 5 mM A/C/G, 2 mM T; 2 µl of ImPromII Reverse transcriptase (Promega), and 2 µl of dUTP-Cy3 or dUTP-Cy5 was added to the RNA/primer mixture. For the synthesis of fluorescent cDNA, the mixture was incubated at 25ºC for 5 min and 40ºC for 3 hrs in a programmable thermocycler PC-100 Programmable Thermal Controller (MJ Research Inc). Three different reverse transcription reactions were carried out for the three biological samples, which were subsequently hybridized on different slides (triplicates) to minimize variation in gene expression not related to infection. After cDNA synthesis, a step of RNA degradation was included by adding 2.5 µl of EDTA 0.5M (pH 8.0) and 5.0 µl of NaOH 1M and incubating at 37ºC for 4 min. All the reagents and primers were removed by using Microcon YM-3 columns (Ambion). The labeled cDNA was stored at 4ºC for subsequent hybridization.

Hybridization and washing

Arrays were hybridized and washed in a Gene-Tac Hybridization Station (Genomic Solutions, Ann Arbor, MI, USA). Initially, the arrays were denatured at 65ºC for 5 min. A solution containing 8 µl of blocking reagent (GE Healthcare, Piscataway, NJ, USA RPN 3601), 1 µl of SSC 2x, 5,5 µl of 2% w/v SDS, 10 pmoles of each cDNA labeled with the fluorescent dyes Cy3 or Cy5, in a final volume of 11 µl, was heated at 95ºC for 2 min. This solution was injected into the hybridization chamber to cover the arrays and the hybridization was performed for 16 hrs at 42ºC. Afterwards, the slides were washed at 25ºC with the following solutions: 0.5% (w/v) SDS in SSC of different concentrations: 2x SSC (high stringency), 0.5x SSC (medium stringency) and 0.05x SSC (low stringency). All the washing steps consisted of 1 cycle of solution flow (1s) and incubation (2s). The slides were dried for 15 min and subjected to fluorescence detection.

Image acquisition and data analyses

The slides were subjected to fluorescence reading in an Axon Gene Pix 4000B Microarray Scanner (Molecular Devices) under different wavelengths: 550 nm (Cy3) and 650 nm (Cy5). The location and identity of each gene on the slide was defined in a text file, created with the aid of the Clone-Tracker2 program (Biodiscovery) by Travensolo et al. (2008) and Travensolo et al. (2009). The signal was quantified through ImaGene software (v.4.1, BioDisco-very), in which two images from the Cy3 and Cy5 fluorescent dyes were overlapped and the spots classified according to morphology and intensity. The computer displays an electronic symbol as a false-color image where a red or green spot corresponds to expression of a gene in sample 1 or 2, respectively, while a yellow-orange spot indicates that the gene was expressed at similar levels in both samples.

Normalization and statistical analysis were carried out using the R Program (www.r-project.org) (R. Development core team, 2009) and the limma package (http// bioconductor.org/packages/2.4/bioc/html/limma.html) (Smyth, 2005). The “Array Weights” software was run to check the quality of the arrays and the background signal was discounted from the signal of each spot using the method normexp with an offset equal to 1 (Ritchie et al. 2006) and normalized by the “loess” method (Yang et al. 2001; Yang et al. 2002; Smyth and Speed, 2003; Ritchie et al. 2007). After normalization, the Students Q-Q test was applied and 400 genes were identified as differentially expressed (Smyth, 2004; Smyth et al. 2005). Genes above a threshold value of LogFC=1 were considered differentially expressed in Navelina, while those below LogFC=-1 were considered differentially expressed in Pera (Figure 1).

Microarray validation

The genes used for microchip validation by RT-qPCR were all related to pathogenicity functions with high and medium values of LogFC (in Pera or Navelina ISA 315).

Fasta sequences were downloaded from the Gene Bank of the Xylella fastidiosa Genome Project (http://aeg.lbi.ic.unicamp.br/xf/) and primer sequences were designed using Primer Express Software for Real Time PCR 3.0, provided by the real time PCR equipment of Applied Biosystems, model 7500. A total of fifty candidate sequences were analyzed for each gene, in order to exclude secondary structures like “hairpins” and possible primer dimer formation. Table 1 presents the sequences of the primers used in the RT-qPCR analysis. All primers were validated for amplification efficiency and specificity following the “MIQE guidelines: minimum information for publication of quantitative real-time PCR experiments” (Bustin et al. 2009).

Reverse transcription (RT) was carried out with 2 µg of X. fastidiosa total RNA from the three biological replicates primed with 1.5 µl. of random hexamers pdN6 (Invitrogen) using the reverse transcription enzyme Improm II (Promega), following kit instructions. Total RNA used in reverse transcription was composed of an equal amount of the RNA from each individual plant used in microarray analysis. Two hundred ng of the resulting cDNA was subjected to quantitative relative expression PCR (RT qPCR) on the Applied Biosystems model 7500 System, using 400 nM (each) of the forward and reverse primers and 6.25 µl of SYBR Green PCR Master Mix (Applied Biosystems) in a final volume of 25 µl. The default thermocycler program was used for all genes and qrPCR assays were performed in triplicate for each primer pair. Threshold values were normalized according to the threshold cycle of ORF XF0305 (NuoA gene) which is expressed at similar levels under all conditions tested according to our microarray data. The change in the expression of each gene was calculated using the 2-∆∆Ct method. The ORF used as the endogenous control, NuoA, codes for the NADH-ubiquinone oxidoreductase (subunit NQ07) which belongs to an operon with many subunits of NADH dehydrogenase I, a multimeric protein. This enzyme participates in aerobic respiration, and would therefore be expressed in a constant and invariable way in bacteria. For this reason it is used as a standard to normalize variations in the gene expression of other genes. It had no significant variation under the experimental conditions of this study in the microarray (LogFC:-0.2)

Results and Discussion

Global expression profile

Among the genes analyzed, approximately 15.4% (400 genes) were detected as differentially expressed within the plants studied. From these, 24 genes (6%) showed higher expression in susceptible Pera, and 151 genes (37.75%) in asymptomatic Navelina. Differentially expressed genes were distributed among Xylella´s functional categories (http://unicamp.lbi.ic.unicamp.br/xf) (Figure 2). According to the results of the microarray analysis, bacteria infecting the susceptible variety Pera expressed a higher number of genes than in the asymptomatic plant, the cultivar Navelina ISA315. We found differentially expressed genes within the eight Xylella fastidiosa gene categories in both hosts. Just one bacterial gene differentially expressed in Pera was included in the ninth category, “ORFs with undefined category”, but, as it is described as an integral membrane protein, we included it in category IV “Cell Structure” for a better graphic visualization. The highest percentage of differentially expressed genes in both varieties was found for hypothetical proteins (47.5%) corresponding to 190 genes (Figure 3). All genes described in “Results” are described in Table 2, including gene ID, name of gene, category, gene product and P number. Table 3 shows the abbreviations corresponding to all gene categories presented in Table 2.

Many genes distributed among the above mentioned categories are involved in bacterial survival, and influence pathogenicity. A complex network of signals between plant, bacteria and the environment, activates or represses different genes. As can be observed in Figure 2 and Figure 3, hypothetical and conserved hypothetical proteins represented 47% of the total of genes found differentially expressed in both varieties, including a percentage of 29.25% in Pera and a percentage of 18.25% in Navelina ISA 315 (corresponding to 117 and 73 genes respectively). It is important to mention that a significant number of hypothetical proteins differentially expressed in Pera presented a very high value of LogFC. This fact suggests a possible role, still unexplored, for many of them in pathogenicity. One of these proteins, ORF XFa0032, still classified as a conserved hypothetical protein at the Xylella fastidiosa Genome Project site (http://www.lbi.ic.unicamp.br/xf/), was identified as a hydoxynitrile lyase related to the process of cyanogenesis by Caruso et al. (2009). Its possible role in Navelina ISA 315 cultivar is further discussed below.

Differentially expressed genes involved in the response to stress conditions

Heat shock response. Three genes coding chaperons were induced in Navelina ISA 315, DNAK, DNA J and HTPG (with LogFC values of 3.20, 1.32 and 1.71 respectively) whereas only one was induced in Pera, the gene HTRA (LogFC -1.72). The complete genome sequences of different Xylella strains have been determined (Simpson et al. 2000; Bhattacharyya et al. 2002; Weiner III et al. 2003), revealing the presence of genes encoding different HSPs as well as the heat shock sigma factor σ32. However, little is known about stress responses and gene regulation in these organisms (Koide et al. 2006) and a more direct assignment of gene function and evaluation of its association with pathogenesis require the analysis of mutants but only a few were obtained for CVC strains by homologous recombination (Da Silva Neto et al. 2002; Da Silva Neto et al. 2006; Gaurivaud et al. 2002). Koide et al. (2006) carried out a global analysis of heat shock response in vitro and determined genes involved in the citrus strain 9a5c, after exposure to 40ºC. In addition to the induction of classical Hsp genes, they observed the up-regulation of several genes related to pathogenesis and adaptation, as well as the repression of fimbria-related genes and genes involved in energy metabolism, revealing a complex network of genes that work together in response to heat stress in Xylella fastidiosa.

DNAK, the Escherichia coli Hsp70 homologue, is an ATP-dependent molecular chaperone that acts in conjunction with the co-chaperones, DNAJ and GrpE, to mediate protein folding and remodeling reactions in the cell. The DNAK and Hsp70 chaperone systems participate in a wide variety of cellular processes in both normal and stressed cells, including nascent protein folding, protein trafficking across intracellular membranes, proteolysis, assembly of multiprotein structures, disassembly of protein aggregates, cell division, DNA replication of several phages and plasmids, and regulation of the heat shock response (Bukau and Horwich, 1998; Bukau et al. 2000; Frydman, 2001; Kim et al. 2002). The Xylella fastidiosa HTPG gene, also diferentially expressed in Navelina ISA 315, has the strongest similarity with Escherichia coli HTPG high temperature protein G, which is a molecular chaperone, has ATPase activity and belongs to the heat shock protein 90 family. The gene HTRA, differentially expressed in Pera, codes for a protein very similar to the E. coli heat shock protein htrA, a periplasmic serine protease, which is required at high temperatures and is involved in degradation of damaged proteins. This protein is essential for bacterial survival at temperatures higher than 42ºC (Xylella fastidiosa Genome Project http://www.lbi.ic.unicamp.br/xf/). It is possible that this activity helps the bacteria to survive at high temperatures within Pera degrading damaged polypeptides.

The fact that more genes coding for chaperones were induced in Navelina ISA 315 suggested that bacteria could be under more stress within asymptomatic plants, maybe because of plant defenses or competition from endophytic microorganisms. Experiments had shown that Xylella fastidiosa growth was stimulated in vitro by Methylobacterium extorquens and inhibited by Curtobacterium flaccumaciens, which demonstrated that real interaction occurs between endophytic microorganisms and the pathogen (Lacava et al. 2004). The DNAK/DNAJ/GrpE system could be more active in the function of re-assembling disaggregated polypeptides in Navelina ISA 315. Defense mechanisms in the Pera variety could be less effective, generating fewer damaged polypeptides, and the DNAK/DNAJ/GrpE system may therefore not be essential for bacterial survival. In contrast, HTRA chaperone is more active in Xylella fastidiosa within Navelina ISA 315. A different heat shock response could be activated in each variety depending on environmental conditions within the plant, specific plant defenses or endophytic competition.

Antioxidant response, toxins and detoxification proteins. During X. fastidiosa infection, elicitors of plant defense responses synthesized by the bacterium and/or released from partial plant cell wall degradation may induce the generation of reactive oxygen species (ROS), such as superoxide and hydrogen peroxide which may inhibit bacterial growth and disease development (Wojtaszec, 1997). Several genes involved in ROS detoxification have been identified in X. fastidiosa including genes that code for catalase, superoxide dismutase, glutathione peroxidase, and genes oxyR and Ohr, suggesting that X. fastidiosa can respond to oxidative stresses (Simpson et al. 2000).

Four genes were differentially expressed in the cultivar Navelina ISA 315, within subcategory “Toxin production and detoxification” of the category “Pathogenicity, virulence and adaptation”: the genes KSGA, cpeB, st3F7.11 and pbp, coding for a dimethyladenosine transferase, catalase/peroxidase, peptide synthase and the protein beta-lactamase, respectively (Log FC values: 1.3; 1.5; 1.5 and 2.2). In Pera, the genes FRPC, AHPC, cutC, czcA and TolC were differentially expressed with LogFC values of -3.4, -2.0 and -1.4, -1.49 and -1.48. The last 3 mentioned genes, cutC, czcA and TolC, involved in copper resistance were also differentially expressed in Pera. Rodrigues et al. (2008) also found a higher expression of these genes in cells forming biofilms, which also suggests a higher resistance to copper related to biofilm formation in the symptomatic variety.

The gene coding for catalase/peroxidase was differentially expressed in Navelina ISA 315 whereas subunit C of alkyl-hydroperoxide reductase was differentially expressed in Pera (-2.0 LogFC). Alkyl hydroperoxide reductase subunit C (AhpC) is the catalytic subunit responsible for alkyl-peroxide metabolism in Xanthomonas. The increase in catalase activity was a compensatory response to lack of AhpC, through a complex regulation response which differs from other bacteria (Mongkolsuk et al. 2000). Higher expression of catalase/peroxidase in Navelina ISA 315 could compensate for lower expression of AhpC, whereas in Pera AhpC is more expressed. Probably, Xylella fastidiosa is reacting in a different way within symptomatic and asymptomatic plants according to the ROS plant response which is reflected by differentially expressed antioxidant proteins.

Among the toxin-encoding genes identified in X. fastidiosa, four are hemolysin-like genes belonging to the RTX family (Simpson et al. 2000). Related to this family, the FRPC gene, which codes for a hemolysin-type calcium binding protein, is induced in Pera with a high LogFC value (-3.47). This gene’s product is similar to other hemolysins of the RTX family, being located in the outer membrane and secreted. It has a GLY-rich region probably involved in binding calcium, which is required for target cell binding or cytolytic activity. This type of toxins forms pores that interact with the plasma membrane of host cells, causing leakage of cellular substances (Lally et al. 1999). Interestingly, one conserved hypothetical protein (XFa0032) of 27.1 kDa recently characterized as a α-hydroxynitrile lyase, was differentially expressed in cultivar Navelina ISA 315 with a high logFC value of 2.4. Alpha-hydroxynitrile lyases (HNLs) are plant defense enzymes that protect plant cells from hervibores or microbial attack, producing HCN through a process called cyanogenesis (Caruso et al. 2009). Xylella fastidiosa could be producing this enzyme in Navelina ISA 315, responding to plant defenses. This plant, or maybe an endophyte, could be secreting more HCN, which is a toxic compound to defend itself from Xylella fastidiosa, or from other microorganisms within the plant. Xylella fastidiosa in turn, codes for more HNLs to detoxify the environment and survive.

Differentially expressed genes involved in survival in atypical conditions and related mechanisms

The genes MDOH (XF1623) and MDOG (XF2682) included in subcategory “Adaptation to atypical conditions” of the category “Pathogenicity, Virulence and Adaptation” are responsible for the synthesis of osmoregulated periplasmic glycan proteins (OPGs). These genes were differentially expressed in cultivar Navelina ISA 315 and Pera with LogFC values of 1.55 and -1.55, respectively. Osmoregulated periplasmic glycan proteins (OPGs) are important membrane components of Gram negative cells which are also involved in survival at low nutrient concentrations as well as inadequate osmotic conditions, especially during the process of infection of eukaryotic hosts. Normally they are associated with phytopathogenic processes (Bohin, 2000). There is probably an opposite regulation at the transcriptional level of those genes within the varieties Pera and Navelina ISA 315, leading to different concentrations of the proteins, which then interfere in host pathogen interactions and determine the degree of pathogenicity and symptom development.

The gene SspB, highly induced in Pera variety (-3.69 LogFC) encodes the stringent starvation protein B. It was recently demonstrated that in Escherichia coli, the SspA gene has an important function in stress response (including acid tolerance) through negative regulation of the global regulator H-NS, which regulates multiple defense systems (Cowtan, 1994; Vassylyev et al. 2002). Williams et al. (2006) reported the starvation-induced expression of SspA and SspB in Escherichia coli and their functions for survival during growth and prolonged starvation. The inactivation of the SspA gene blocked SspB expression. In our microarray experiment, high expression of SspB in Pera could indicate that the bacterium is surviving in the absence of nutrients in the stationary phase and inducing acid tolerance within the xylem vessels of the symptomatic plant.

Many regulatory genes, involved in many metabolic cascades, were differentially expressed in this experiment. The gene rpfC was differentially expressed in Pera variety with a value of -1.59 LogFC. The RpfC gene codes for Rpf C protein, the sensor component of the two-component system RpfF/RpfC, which is part of the Rpf (regulation of pathogenicity factors) cluster. Chatterjee et al. (2008) proposed a model in which the same sensor molecule (rpfC) is required both to develop virulence and insect transmission. The induction of the RpfC gene in the susceptible variety is related to a gene regulation cascade which involves biofilm formation. Perhaps different biofilm maturation or denser composition could occur within the Pera variety, which would be consistent with the severe symptoms observed, and probably with insect transmission. However, this hypothesis should be confirmed with further studies.

As shown in Table 2, other regulatory genes like ALGR, colR, algZ and VapD were differentially expressed in variety Pera, while colR and ColS were induced in Navelina ISA 315. This fact suggests changes in gene expression of sensor and regulatory molecules related to the different conditions within each plant.

Differentially expressed genes involved in iron homeostasis

Among some putative ORFs identified as possibly related to iron transport (Silva-Stenico et al. 2005), only the ORFs XF0599 (YBIL gene) and XF1038 (St3F7.11 gene) were differentially expressed in this experiment in Pera (-1.7 LogFC) and Navelina ISA315 (1.5 LogFC), respectively. The gene YBIL codes for a TonB dependent receptor and St3F7.11 for a peptide synthase enzyme. The TonB gene product is involved in iron-siderophore transport across the membrane.

When pathogenic bacteria are challenged with limiting iron concentrations, they can mantain homeostasis of this metal by releasing iron from intracellular reservoirs and enhancing the expression of the systems involved in iron uptake (Ratledge and Dover, 2000). This system allows bacterial survival in nutrient-poor environments, as iron is an essential cofactor for many proteins involved in electron transfer and redox reactions. There may be differences in iron availability between the varieties we studied which could influence bacterial response related to gene expression. Levels of expression of genes related to iron uptake may be induced depending on limiting iron conditions within the xylem vessels. As the YBIL gene was induced in Pera, probably iron is a more limited factor within this variety than in Navelina ISA 315. That would be consistent with the symptoms observed in Pera. Simpson et al. (2000) proposed that the absorption of iron and probably other metallic ions like manganese causes the depletion of essential micronutrients in the xylem, leading to the typical symptoms of variegated clorosis.

Zaini et al. (2009) reported a constant level of expression of Xylella fastidiosa transcripts in limiting iron conditions using microarray analysis. They observed the modulation of many CDS related to regulatory functions, pathogenicity and cellular structure, toxins and pili/ fimbriae. We found that many genes involved in pathogenicity were also induced in Pera. The complex TonB-ExB-ExbD is also likely to be involved in greater uptake of iron.

Genes related to mobile genetic elements: phages, prophages and plasmids

The genes VapD, differentially expressed in Pera and SCJ2116, in Navelina ISA 315, are related to virulence and detoxification, respectively, and are both included in mobile genetic elements. The presence of virulence genes within prophage sequences is strong evidence of the important role of horizontal transfer of genes via bacteriophages in the definition of bacterial phenotype (Hendrix et al. 1999). Many of the differentially expressed ORFs of this experiment were classified as hypothetical proteins, suggesting a high transcriptional activity in this region, perhaps related to pathogenicity or adaptability, that is still unexplored (Table 2). Figure 2 shows that more genes related to mobile genetic elements are differentially expressed in Pera than in Navelina ISA 315. Nunes et al. (2003) showed that Xylella fastidiosa has the biggest “pool” of laterally transmitted elements characterized until then (like prophages, plasmids and genomic islands) which contributed up to 18% of the total genome. Transcriptome analysis confirmed that these elements are transcriptionally active, possibly responding to environmental signals. The differentially expressed ORFs within these elements, found in this experiment, also confirmed this activity. A deeper study of laterally transferred elements would help to explain the adaptability and capacity of the bacteria to infect such a wide host range.

DNA repair, DNA restriction and modification

More genes related to DNA repair, restriction and modification systems were expressed at higher levels in Pera variety (Table 1 and Table 3), which could suggest a higher response to DNA damage caused by plant defense systems. Xylella fastidiosa could be more actively replicating in the susceptible variety, therefore overcoming the stress generated by plant defense mechanisms within this variety. Different survival strategies involving different processes may be adopted by the bacteria depending on plant environment.

Validation of the microarray

Primers for RT-qPCR were all validated for amplification efficiency and specificity following the “MIQE guidelines: minimum information for publication of quantitative real-time PCR experiments” (Bustin et al. 2009). Efficiencies of about 100% ± 20% were found for all of them. The RT-qPCR analysis of the genes ChpA, ColS, DNAJ, HTPG, chpA, Pil U, Vap D confirmed the results obtained by the microarray analysis in the system studied. Figure 4 shows relative levels of increased gene expression, given as the multiplication factor (number of times) in Navelina ISA 315 compared to Pera, for each gene, with both microarray and RT-qPCR analysis.

The inherent variation of both microarray and RT-qrPCR due to processing of material and normalization methods can influence the correlation between results obtained. It is more common to find in the literature the simple statement that results were validated, often with no, or extremely low, reported correlations (Freeman et al. 1999; Bustin, 2002; Chyaqui et al. 2002; Wurumbach et al. 2003). For that reason, when the same direction of expression in the same variety was found with both methods, this was considered validation.

Healthy Pera and Navelina plants did not reveal any fluorescence due to possible non-specific hybridization with plant genome. Also, one slide without hybridization was scanned in order to detect any non-specific hybridization signal, but similar background fluorescence to the healthy control plants was obtained. These results, together with microarray internal controls and qrtPCR results, confirmed that the hybridization signal from bacteria infecting both Navelina ISA 315 and Pera, is the result of X. fastidiosa genes hybridization in the system studied.

Concluding Remarks

The finding that more bacterial genes were differentially expressed in the susceptible variety was to be expected since the bacteria are probably more actively migrating through vessels, blocking water and nutrient flow and degrading plant tissues while colonizing different parts of the susceptible variety. Even though bacterial concentration was found to be high in cultivar Navelina ISA315, it was probably not blocking the xylem flow in vessels as it was in Pera, since it did not cause any symptoms. Plant defenses or competition from other microorganisms may be more efficient in this variety, thus generating higher stress in the bacteria, suggested by differentially expressed “heat shock” genes, as well as the above-mentioned detoxification systems. Perhaps in Navelina ISA 315, cells are not so embedded in biofilms, and move more freely through vessels, therefore having more plankton-like characteristics and probably being exposed to much more stress without biofilm protection. Also, biofilms could have different compositions (Martins et al. 2010) in each variety.

Different bacterial gene expression patterns probably reflect different bacterial survival strategies within each variety involving differences in biofilm composition, sensing systems and migrating patterns, as well as different responses to stress conditions like detoxification proteins and chaperones. Different levels of gene expression could be the key to understanding the ability of X. fastidiosa to colonize different hosts and survive in such adverse conditions as xylem vessels. Moreover, bacterial non-coding small RNAs (sRNA) have attracted considerable attention due to their ubiquitous nature and roles in controlling numerous cellular processes including survival, adaptation and pathogenesis. However, little is known about small RNA in this bacterium, in spite of the fact that several whole genome sequences of X. fastidiosa have been published. Results from BLAST analysis showed that 34 small RNA genes were shared by all four X. fastidiosa strains. Species-, subspecies- and pathotype-specific small RNAs (size from 40 to 350 pb) were also identified (Chen and Huang, 2011). A clearer idea of regulation processes and survival strategies could be gained by including these promising newly identified small RNA sequences in future transcriptome analyses.

It is important to consider the adaptation potential of X. fastidiosa in CVC control strategies and breeding programs when using the cultivar Navelina ISA 315 because, even without symptoms, it could be a source of infection due to the concentration of bacteria inside the plant.

Spotted microrrays proved to be a powerful technique to analyze changes in gene expression for Xylella fastidiosa infecting citrus plants (in vivo system). The technique could also be used to monitor gene expression changes in time course experiments in natural habitats. Still more research is needed to understand the disease process and elucidate the environmental signals which could coordinate bacterial behaviour in terms of gene expression and regulation of cell density, as well as plant-pathogen interactions. Nevertheless, transcriptomic analysis could give many clues about the different processes involved in pathogenicity, like the heat shock response and some of the survival mechanisms presented in this paper. Many other pathogenicity and regulatory genes were differentially expressed in this experiment within each variety which has not been discussed in this article.

Other techniques and studies could be applied in order to achieve a better and deeper understanding of this system, like gene complementation analyses, assessing the diversity of other microorganisms within the plant, as well as electron microscopy. In any case, transcriptomic analysis gave a general view of different biochemical processes activated or repressed in natural conditions, considering the inherent variables of an in vivo system. On the other hand, it would be interesting to complement this study with transcriptomic analysis of the plant hosts studied here, which could reveal much more information about plant-pathogen interactions.

To our knowledge, this is one of the first gene expression analyses of X. fastidiosa in vivo within the plant, which gave a more accurate image of this pathogen in nature and its ability to live in different environments (hosts), without the usual bias introduced by in vitro culture conditions.

Acknowledgements

In first place, I would like to thank Professors PhD Eliana G. M. Lemos and PhD Jackson Marcondes, my thesis advisor and co-advisor. I am especially grateful to Ing. Agr. MSc Carmen Goñi, Director of the INIA Citrus Program of Uruguay at the time I started my PhD studies and to the current Director Ing. Agr. PhD Fernando Rivas for supporting this research. I would also like to thank the current and former Directors of the Biotechnology Unit of INIA, Ing. Agr. PhD Marco Dalla Rizza and Ing. Agr. PhD Dan Piestun, as well as Ing. Agr. MSc. Elena Perez for supporting my work at INIA during my PhD studies, and for their valuable suggestions about this research.

References

AMORIM, L.; BERGAMIN FILHO, A.; PALAZZO, D.; BASSANEZ, R.B.; GODOY, C.V. and TORRES, G.A.M. (1993). Clorose variegada dos citros: Uma escala diagrammática para avaliação da doença. Fitopatologia Brasileira, vol. 17, no. 1, p. 42-48.

AUSUBEL, F.M.; BRENT, R.; KINGSTON, R.E.; MOORE, D.D.; SEIDMAN, J.G.; SMITH, J.A. and STRUHL, K. (1987). Current Protocols in Molecular Biology. Greene Publishing Associates/Wiley Interscience, New York, 5300 p. ISBN 9780471142720.

BHATTACHARYYA, A.; STILWAGEN, S.; IVANOVA, N.; D’SOUZA, M.; BERNAL, A.; LYKIDIS, A.; KAPATRAL, V.; ANDERSON, I.; LARSEN, N.; LOS, T.; REZNIK, G.; SELKOV JR., E.; WALUNAS, T.L.; FEIL, H.; FEIL, W.S.; PURCELL, A.; LASSEZ, J.-L.; HAWKINS, T.L.; HASELKORN, R.; OVERBEEK, R.; PREDKI, P.F. and KYRPIDES, N.C. (2002). Whole-genome comparative analysis of three phytopathogenic Xylella fastidiosa strains. Proceedings of the National Academy of Sciences of the United States of America, vol. 99, no. 19, p. 12403-12408. [CrossRef]

BOHIN, J.-P. (2000). Osmoregulated periplasmic glucans in proteobacteria. FEMS Microbiology Letters, vol. 186, no. 1, p. 11-19. [CrossRef]

BUSTIN, S.A. (2002). Quantification of mRNA using real-time reverse transcription PCR (RT-PCR): Trends and problems. Journal of Molecular Endocrinology, vol. 29, p. 23-39. [CrossRef]

BUSTIN, S.; BENES, V.; GARSON, J.A.; HELLEMANS, J.; HUGGETT, J.; KUBISTA, M.; MUELLER, R.; NOLAN, T.; PFAFFL, M.W.; SHIPLEY, G.L.; VANDESOMPELE, J. and WITTWER, C.T. (2009). The MIQE guidelines: Minimum information for publication of quantitative real-time PCR experiments. Clinical Chemistry, vol. 55, no. 4, p. 611-622. [CrossRef]

BUKAU, B. and HORWICH, A.L. (1998). The Hsp70 and Hsp60 chaperone machines. Cell, vol. 92, no. 3, p. 351-366. [CrossRef]

BUKAU, B.; DEUERLING, E.; PFUND, C. and CRAIG, E.A. (2000). Getting newly synthesized proteins into shape. Cell, vol. 101, no. 2, p. 119-122. [CrossRef]

CAMILLI, A. and BASSLER, B.L. (2006). Bacterial small-molecule signaling pathways. Science, vol. 311, no. 5764, p. 1133-1116. [CrossRef]

CARUSO, C.S.; DE FATIMA TRAVENSOLO, R.; DE CAMPUS BICUDO, R.; DE MACEDO LEMOS, E.G.; ULIAN DE ARAÚJO, A.P. and CARRILHO, E. (2009). α-Hydroxynitrile lyase protein from Xylella fastidiosa: Cloning, expression, and characterization. Microbial Pathogenesis, vol. 47, no. 3, p. 118-127. [CrossRef]

CHATTERJEE, S.; WISTROM, C. and LINDOW, S.E. (2008). A cell-cell signaling sensor is required for virulence and insect transmission of Xylella fastidiosa. Proceedings of the National Academy of Sciences of the United States of America, vol. 105, no. 7, p. 2670-2675. [CrossRef]

CHEN, J. and HUANG, H. (2011). Searching for small RNA genes in Xylella fastidiosa genomes. American Phytopathological Society Abstracts, vol. 101, no. S34.

CHYAQUI, R.F.; BONNER, R.F.; BEST, C.J.M; GILLESPIE, J.W.; FLAIG, M.J.; HEWITT, S.M.; PHILLIPS, J.L.; KRIZMAN, D.B.; TANGREA, M.A.; AHRAM, M.; LINEHAN, W.M.; KNEZEVIC, V. and EMMERT-BUCK, M.R. (2002). Post-analysis follow-up and validation of microarray experiments. Nature Genetics, vol. 32, no. 4, p. 509-514. [CrossRef]

COWTAN, K. (1994). dm: An automated procedure for phase improvement by density modification. Joint CCP4 and ESF-EACBM Newsletter on Protein Crystallography, vol. 31, p. 34-38.

DA SILVA, F.R.; VETTORE, A.L.; KEMPER, E.L.; LEITE, A. and ARRUDA, P. (2001). Fastidian gum: The Xylella fastidiosa exopolysaccharide possibly involved in bacterial pathogenicity. FEMS Microbiology Letters, vol. 203, no. 2, p. 165-171. [CrossRef]

DA SILVA NETO, J.F.; KOIDE, T.; GOMES, S.L. and MARQUES, M.V. (2002). Site-directed gene disruption in Xylella fastidiosa. FEMS Microbiology Letters, vol. 210, no. 1, p. 105-110. [CrossRef]

DA SILVA NETO, J.F.; KOIDE, T.; GOMES, S.L. and MARQUES, M.V. (2006). The single extracytoplasmic-function sigma factor of Xylella fastidiosa is involved in the heat shock response and presents an unusual regulatory mechanism. Journal of Bacteriology, vol. 189, no. 2, p. 551-560. [CrossRef]

DE LIMA, J.E.O.; MIRANDA, V.S.; HARTUNG, J.S.; BRLANSKY, R.H.; COUTINHO, A.; ROBERTO, S.R. and CARLOS, E.F. (1998). Coffee leaf scorch bacterium: Axenic culture, pathogenicity, and comparison with Xylella fastidiosa of citrus. Plant Disease, vol. 82, no. 1, p. 94-97. [CrossRef]

FADEL, A.L.; COLETTA FILHO, H.D.; RODRIGUES, C.M.; FEDERICI, M.T. and STUCHI, E.S. (2010). Reação da cultivar Navelina ISA 315 (Citrus sinensis L. Osb.), à clorose variegada dos citros em condições de campo. In: Proceedings of the XXI Brazilian Fruit Congress. Frutas: saúde, inovação e responsabilidade: anais. (17th - 22 th October, 2010, Natal, RN, Brazil). Sociedade Brasileira de Fruticultura.

FREEMAN, W.M.; WALKER, S.J. and VRANA, K.E. (1999). Quantitative RT-PCR: Pitfalls and potential. BioTechniques, vol. 26, no. 1, p. 112-125.

FRYDMAN, J. (2001). Folding of newly translated proteins in vivo: The role of molecular chaperones. Annual Review of Biochemistry, vol. 70, p. 603-647. [CrossRef]

GAURIVAUD, P.; SOUZA, L.C.A.; VIRGILIO, A.C.D.; MARIANO, A.G.; PALMA, R.R. and MONTEIRO, P.B. (2002). Gene disruption by homologous recombination in the Xylella fastidiosa citrus variegated chlorosis strain. Applied and Environmental Microbiology, vol. 68, no. 9, p. 4658-4665. [CrossRef]

HENDRIX, R.W.; SMITH, M.C.M.; BURNS, R.N.; FORD, M.E. and HATFULL, G.F. (1999). Evolutionary relationships among diverse bacteriophages and prophages: All the world’s a phage. Proceedings of the National Academy of Sciences of the United States of America, vol. 96, no. 5, p. 2192-2197. [CrossRef]

HOPKINS, D.L. (1989). Xylella fastidiosa: Xylem-limited bacterial pathogen of plants. Annual Review of Phytopathology, vol. 27, p. 271-290. [CrossRef]

KILLINY, N. and ALMEIDA, R.P.P. (2009). Host structural carbohydrate induces vector transmission of a bacterial plant pathogen. Proceedings of the National Academy of Sciences of the United States of America, vol. 106, no. 52, p. 22416-22420. [CrossRef]

KIM, S.-Y.; SHARMA, S.; HOSKINS, J.R. and WICKNER, S. (2002). Interaction of the DnaK and DnaJ chaperone system with a native substrate, P1 RepA. The Journal of Biological Chemistry, vol. 277; no. 47, p. 44778-44783. [CrossRef]

KOIDE, T.; VÊNCIO, R.Z.N. and GOMES, S.L. (2006). Global gene expression analysis of the heat shock response in the phytopathogen Xylella fastidiosa. Journal of Bacteriology, vol. 118, no. 16, p. 5821-5830. [CrossRef]

LACAVA, P.T.; ARAUJO, W.L.; MARCON, J.; MACCHERONI JR, W. and AZEVEDO. J.L. (2004). Interaction between endophytic bacteria from citrus plants and the phytopathogenic bacteria Xylella fastidiosa, causal agent of citrus-variegated chlorosis. Letters in Applied Microbiology, vol. 39, no. 1, p. 55-59. [CrossRef]

LALLY, E.T.; HILL, R.B.; KIEBA, I.R. and KOROSTOFF, J. (1999). The interaction between RTX toxins and target cells. Trends in Microbiology, vol. 7, no. 9, p. 356-361. [CrossRef]

MARTINS, L.S.; PICCHI, S.C. and LEMOS, M.V.F. (2010). Monosaccharides composition of biofilm produced by Xylella fastidiosa wild type and rpf mutants. Research Journal of Agriculture and Biological Sciences, vol. 6, no. 6, p. 1029-1035.

MONGKOLSUK, S.; WHANGSUK, W.; VATTANAVIBOON, P.; LOPRASERT, S. and FUANGTHONG, M. (2000). A Xanthomonas alkyl hydroperoxide reductase subunit C (ahpC) mutant showed an altered peroxide stress response and complex regulation of the compensatory response of peroxide detoxification enzymes. Journal of Bacteriology, vol. 182, no. 23, p. 6845-6849. [CrossRef]

MONTEIRO-VITTORELLO, C.V.; DE OLIVEIRA, M.; ZERILLO, M.M.; VARANI, A.M.; CIVEROLO, E. and VAN SLUYS, M. (2005). Xylella and Xanthomonas mobil'omics. OMICS: A Journal of Integrative Biology, vol. 9, no. 2, p. 146-158. [CrossRef]

NUNES, L.R.; ROSATO, Y.B.; MUTO, N.H.; YANAI, G.M.; DA SILVA, V.S.; LEITE, D.B.; GONÇALVES, E.R.; DE SOUZA, A.A.; COLETTA-FILHO, H.D.; MACHADO, M.A.; LOPES, S.A. and DE OLIVEIRA, R.C. (2003). Microarray analyses of Xylella fastidiosa provide evidence of coordinated transcription control of laterally transferred elements. Genome Research, vol. 13, p. 570-578. [CrossRef]

OLIVEIRA, A.C.; VALLIM, M.A.; SEMIGHINI, C.P.; ARAÚJO, W.L.; GOLDMAN, G.H. and MACHADO, M.A. (2002). Quantification of Xylella fastidiosa from citrus trees by real-time polymerase chain reaction assay. Phytopatology, vol. 92, no. 10, p. 1048-1054. [CrossRef]

PURCEL, A.H. (1997). Xylella fastidiosa, a regional problem or global threat? Journal of Plant Pathology, vol. 79, no. 2, p. 99-105.

RATLEDGE, C. and DOVER, L.G. (2000). Iron metabolism in pathogenic bacteria. Annual Review of Microbiology, vol. 54, p. 881-941. [CrossRef]

R. DEVELOPMENT CORE TEAM (2009). A language and environment for statistical computing. Vienna, Austria. R Foundation for Statistical Computing.

RITCHIE, M.E.; DIYAGAMA, D.; NEILSON, J.; VAN LAAR, R.; DOBROVIC, A.; HOLLOWAY, A. and SMYTH, G.K. (2006). Empirical array quality weights in the analysis of microarray data. BMC Bioinformatics, vol. 7, p. 261. [CrossRef]

RITCHIE, M.E.; SILVER, J.; OSHLACK, A.; HOLMES, M.; DIYAGAMA, D.; HOLLOWAY, A. and SMYTH, G.K. (2007). A comparison of background correction methods for two-colour microarrays. Bioinformatics, vol. 23, no. 20, p. 2700-2707. [CrossRef]

RODRIGUES, C.M.; TAKITA, M.A.; COLETTA-FILHO, H.D.; OLIVATO, J.C.; CASERTA, R.; MACHADO, M.A. and DE SOUZA, A.A. (2008). Copper resistance of biofilm cells of the plant pathogen Xylella fastidiosa. Applied Microbiology and Biotechnology, vol. 77, no. 5, p. 1145-1157. [CrossRef]

SHILD, S.; TAMAYO, R.; NELSON, E.J.; QADRI, F.; CALDERWOOD, S.B. and CAMILLI, A. (2007). Genes induced late in infection increase fitness of Vibrio cholerae after release into the environment. Cell Host & Microbe, vol. 2, no. 4, p. 264-277. [CrossRef]

SILVA-STENICO, M.E.; PACHECO, F.T.H.; RODRIGUES, J.L.M.; CARRILHO, E. and TSAI, S.M. (2005). Growth and siderophore production of Xylella fastidiosa under iron-limited conditions. Microbiological Research, vol. 160, no. 4, p. 429-436. [CrossRef]

SIMPSON, A.J.G.; REINACH, F.C.; ARRUDA, P.; ABREU, F.A.; ACENCIO, M.; ALVARENGA, R.; ALVES, L.M.C.; ARAYA, J.E.; BAIA, G.S.; BAPTISTA, C.S.; BARROS, M.H.; BONACCORSI, E.D.; BORDIN, S.; BOVÉ, J.M.; BRIONES, M.R.S.; BUENO, M.R.P.; CAMARGO, A.A.; CAMARGO, L.E.A.; CARRARO, D.M.; CARRER, H.; COLAUTO, N.B.; COLOMBO, C.; COSTA, F.F.; COSTA, M.C.R.; COSTA-NETO, C.M.; COUTINHO, L.L.; CRISTOFANI, M.; DIAS-NETO, E.; DOCENA, C.; EL-DORRY, H.; FACINCANI, A.P.; FERREIRA, A.J.S.; FERREIRA, V.C.A.; FERRO, J.A.; FRAGA, J.S.; FRANÇA, S.C.; FRANCO, M.C.; FROHME, M.; FURLAN, L.R.; GARNIER, M.; GOLDMAN, G.H.; GOLDMAN, M.H.S.; GOMES, S.L.; GRUBER, A.; HO, P.L.; HOHEISEL, J.D.; JUNQUEIRA, M.L.; KEMPER, E.L.; KITAJIMA, J.P.; KRIEGER, J.E.; KURAMAE, E.E.; LAIGRET, F.; LAMBAIS, M.R.; LEITE, L.C.C.; LEMOS, E.G.M.; LEMOS, M.V.F.; LOPES, S.A.; LOPES, C.R.; MACHADO, J.A.; MACHADO, M.A.; MADEIRA, A.M.B.N.; MADEIRA, H.M.F.; MARINO, C.L.; MARQUES, M.V.; MARTINS, E.A.L.; MARTINS, E.M.F.; MATSUKUMA, A.Y.; MENCK, C.F.M.; MIRACCA, E.C.; MIYAKI, C.Y.; MONTEIRO-VITORELLO, C.B.; MOON, D.H.; NAGAI, D.H.; NASCIMENTO, A.L.T.O.; NETTO, L.E.S.; NHANI, JR, A.; NOBREGA, F.G.; NUNES, L.R.; OLIVEIRA, M.A.; DE OLIVEIRA, M.C.; DE OLIVEIRA, R.C.; PALMIERI, D.A.; PARIS, A.; PEIXOTO, B.R.; PEREIRA, G.A.G.; PEREIRA, JR, H.A.; PESQUERO, J.B.; QUAGGIO, R.B.; ROBERTO, P.G.; RODRIGUES, V.; ROSA, A.J.M.; DE ROSA, JR, V.E.; DE SÁ, R.G.; SANTELLI, R.V.; SAWASAKI, H.E.; DA SILVA, A.C.R.; DA SILVA, A.M.; DA SILVA, F.R.; SILVA, W.A.; DA SILVEIRA, J.F.; SILVESTRI, M.L.Z.; SIQUEIRA, W.J.; DE SOUZA, A.A.; DE SOUZA, A.P.; TERENZI, M.F.; TRUFFI, D.; TSAI, S.M.; TSUHAKO, M.H.; VALLADA, H.; VAN SLUYS, M.A.; VERJOVSKI-ALMEIDA, S.; VETTORE, A.L.; ZAGO, M.A.; ZATZ, M.; MEIDANIS, J. and SETUBAL, J.C. (2000). The genome sequence of the plant pathogen Xylella fastidiosa. Nature, vol. 406, no. 6792, p. 151-157. [CrossRef]

SMYTH, G.K. and SPEED, T.P. (2003). Normalization of cDNA microarray data. Methods, vol. 31, no. 4, p. 265-273. [CrossRef]

SMYTH, G.K. (2004). Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Statistical Applications in Genetics and Molecular Biology, vol. 3, no. 1.

SMYTH, G.K. (2005). Limma: Linear models for microarray data. In: GENTLEMAN, R.; CAREY, V.; HUBER, W.; IRIZARRY, R. and DUDOIT, S. eds. Bioinformatics and Computational Biology Solutions using R and Bioconductor. Springer, New York, p. 397-420.

SMYTH, G.K.; MICHAUD, J. and SCOTT, H.S. (2005). Use of within-array replicate spots for assessing differential expression in microarray experiments. Bioinformatics, vol. 21, no. 9, p. 2067-2075. [CrossRef]

STUCHI, E.S.; SILVA, S.R.; COLLETTA-FILHO, H.D.; FRANCO, D.; CARVALHO, S.A.; SEMPIONATO, O.R.; DONADIO, L.C. and ALVES, K.C.S. (2007). Navelina ISA 315 sweet orange: A citrus variegated chlorosis (CVC) resistant cultivar. In: Conference of the International Organization of Citrus Virologists. (17th, Adana, Turkey). Program & Abstracs, p. 89.

SULLIVAN, J.T. and RONSON, C.W. (1998). Evolution of rhizobia by acquisition of a 500-kb symbiosis island that integrates into a phe-tRNA gene. Proceedings of the National Academy of Sciences of the United States of America, vol. 95, no. 9, p. 5145-5149. [CrossRef]

TRAVENSOLO, R.F.; CARARETO-ALVES, L.M.; COSTA, M.V.C.G.; LOPES, T.J.S.; CARRILHO, E. and LEMOS, E.G.M. (2008). Xylella fastidiosagene expression analysis by DNA microarrays. Genetics and Molecular Biology, vol. 32, no. 2. [CrossRef]

TRAVENSOLO, R.F.; COSTA, M.V.C.G.; CARARETO-ALVES, L.M.; CARRILHO, E. and LEMOS, E.G.M. (2009). Production of DNA microarray and expression analysis of genes from Xylella fastidiosa in different culture media. Brazilian Archives of Biology and Technology, vol. 52, no. 3, p. 555-566. [CrossRef]

VASSYLYEV, D.G.; SEKINE, S.-I.; LAPTENKO, O.; LEE, J.; VASSYLYEVA, M.N.; BORUKHOV, S. and YOKOYAMA, S. (2002). Crystal structure of a bacterial RNA polymerase holoenzyme at 2.6 Å resolution. Nature, vol. 417, no. 6890, p. 712-719. [CrossRef]

WEINER III, J.; ZIMMERMAN, C.-U.; GÖHLMANN, H.W.H. and HERRMANN, R. (2003). Transcription profiles of the bacterium Mycoplasma pneumoniae grown at different temperatures. Nucleic Acids Research, vol. 31, no. 21, p. 6306-6320. [CrossRef]

WELLS, J.M.; RAJU, B.C.; NYLAND, G. and LOWE, S.K. (1981). Medium for isolation and growth of bacteria associated with plum leaf scald and phony peach diseases. Applied and Environmental Microbiology, vol. 42, no. 2, p. 357-363.

WILLIAMS, M.D.; OUYANG, T.X. and FLICKINGER, M.C. (2006). Starvation-induced expression of SspA and SspB: The effects of a null mutation in sspA on Escherichia coli protein synthesis and survival during growth and prolonged starvation. Molecular Microbiology, vol. 11, no. 6, p. 1029-1043. [CrossRef]

WOJTASZEC, P. (1997). Oxidative burst: An early plant response to pathogen infection. Biochemical Journal, vol. 322, p. 681-692.

WONG, K. and GOLDING, G.B. (2003). A phylogenetic analysis of the pSymB replicon from the Sinorhizobium meliloti genome reveals a complex evolutionary history. Canadian Journal of Microbiology, vol. 49, no. 4, p. 269-280. [CrossRef]

WURUMBACH, E.; YUEN, T. and SEALFON, S.C. (2003). Focused microarray analysis. Methods, vol. 31, no. 4, p. 306-316. [CrossRef]

YANG, Y.H.; DUDOIT, S.; LUU, P. and SPEED, T.P. (2001). Normalization for cDNA microarray data. In: BITTNER, M.L.; CHEN, Y.; DORSEL, A.N. and DOUGHERTY, E.R. eds. Microarrays: Optical technologies and informatics. Proceedings of SPIE, vol. 4266, p. 141-152.

YANG, Y.H.; DUDOIT, S.; LUU, P.; LIN, D.M.; PENG, V.; NGAI, J. and SPEED, T.P. (2002). Normalization for cDNA microarray data: A robust composite method addressing single and multiple slide systematic variation. Nucleic Acids Research, vol. 30, no. 4, p. e15. [CrossRef]

ZAINI, P.A.; DE LA FUENTE, L.; HOCH, H.C. and BURR, T.J. (2009). Grapevine xylem sap enhances biofilm development by Xylella fastidiosa. FEMS Microbiology Letters, vol. 295, no. 1, p. 129-134. [CrossRef]

Note: Electronic Journal of Biotechnology is not responsible if on-line references cited on manuscripts are not available any more after the date of publication.