Fortunella margarita Transcriptional Reprogramming Triggered by Xanthomonas citri subsp. citri
© Khalaf et al; licensee BioMed Central Ltd. 2011
Received: 22 June 2011
Accepted: 11 November 2011
Published: 11 November 2011
Citrus canker disease caused by the bacterial pathogen Xanthomonas citri subsp. citri (Xcc) has become endemic in areas where high temperature, rain, humidity, and windy conditions provide a favourable environment for the dissemination of the bacterium. Xcc is pathogenic on many commercial citrus varieties but appears to elicit an incompatible reaction on the citrus relative Fortunella margarita Swing (kumquat), in the form of a very distinct delayed necrotic response. We have developed subtractive libraries enriched in sequences expressed in kumquat leaves during both early and late stages of the disease. The isolated differentially expressed transcripts were subsequently sequenced. Our results demonstrate how the use of microarray expression profiling can help assign roles to previously uncharacterized genes and elucidate plant pathogenesis-response related mechanisms. This can be considered to be a case study in a citrus relative where high throughput technologies were utilized to understand defence mechanisms in Fortunella and citrus at the molecular level.
cDNAs from sequenced kumquat libraries (ESTs) made from subtracted RNA populations, healthy vs. infected, were used to make this microarray. Of 2054 selected genes on a customized array, 317 were differentially expressed (P < 0.05) in Xcc challenged kumquat plants compared to mock-inoculated ones. This study identified components of the incompatible interaction such as reactive oxygen species (ROS) and programmed cell death (PCD). Common defence mechanisms and a number of resistance genes were also identified. In addition, there were a considerable number of differentially regulated genes that had no homologues in the databases. This could be an indication of either a specialized set of genes employed by kumquat in response to canker disease or new defence mechanisms in citrus.
Functional categorization of kumquat Xcc-responsive genes revealed an enhanced defence-related metabolism as well as a number of resistant response-specific genes in the kumquat transcriptome in response to Xcc inoculation. Gene expression profile(s) were analyzed to assemble a comprehensive and inclusive image of the molecular interaction in the kumquat/Xcc system. This was done in order to elucidate molecular mechanisms associated with the development of the hypersensitive response phenotype in kumquat leaves. These data will be used to perform comparisons among citrus species to evaluate means to enhance the host immune responses against bacterial diseases.
Citrus trees are susceptible to a number of diseases with different degrees of economic impact. One of the most severe in terms of economic losses is citrus canker disease (sometimes referred to as Asiatic citrus canker) caused by Xanthomonas citri subsp. citri, (synonym, Xanthomonas axonopodis pv. citri strain A; Xac-A). Xcc is a biotrophic bacterial phytopathogen that belongs to the genus Xanthomonas of the α-subdivision v within Proteobacteria. Susceptibility to citrus canker disease varies among citrus types and relatives, but most of the commercially grown citrus types are susceptible hosts to Xcc . Disease symptoms include canker lesions on the green aerial parts of the plant as well as fruit; infections can result in both foliar and fruit abscission, thereby decreasing the productivity of affected trees. In addition there can be reduced profitability as a result of blemished fruit that can be harvested but not sold in the fresh market.
Plants have evolved multiple defence mechanisms to survive pathogen attacks . The first branch of the indispensable plant innate immunity system is triggered by pathogen-associated molecular patterns (PAMPs) such as the lipopolysaccharides (LPSs), peptidoglycan and bacterial flagellin, as well as the chitin and glucan from fungi. The second branch utilizes the nucleotide-binding site-leucine-rich repeat (NBS-LRR) encoded by R (resistance) genes named the effector-triggered immunity (ETI) . The Xanthomonas spp. phytobacterial pathogens have evolved unique pathogenesis mechanisms to avoid host recognition and suppress host defences [4, 5]. Bacterial effector proteins are delivered via the bacterial type III secretion system (TTSS) into the plant cell to evade recognition by the different plant surveillance systems . These effectors in general contribute to host resistance or susceptibility as well as to modifying host responses. A fundamental element of the ETI in resistant plants is a localized cell collapse or a hypersensitive response (HR) at infection sites in an attempt to restrict the growth of the pathogen [7, 8]. This is a common feature of disease resistant responses in incompatible plant-pathogen, and occasionally some non-host, interactions [9, 10]. Some of the Xanthomonas spp. effector proteins, for instance PthA/AvrBs3, are essential to elicit citrus canker symptoms and if expressed by itself inside host cells, pthA is sufficient to cause symptoms of citrus canker disease [11–15]. In the meantime however, other recent studies show that other types of proteins are injected through the Xcc TTSS and do not necessarily alter the physiological and transcriptional responses to the pathogen in citrus [8, 10, 16, 17].
While certain genes involved in systemic acquired resistance (SAR) have been characterized and used as markers for studying plant defence mechanisms , crosstalk between signals and hormone pathways has also been proposed [19–21]. Consequently, plant resistance is correlated with the activation of a complex network of defence pathways and the response of the host plant to a microbial assault is therefore expected to result in drastic changes in the patterns of gene expression throughout the plant [22, 23].
Kumquats (Fortunella spp.), close relatives to citrus species, are reported to have high levels of field resistance to citrus canker . Previously, we have shown a sharply contrasting phenotype in grapefruit and kumquat when both plants were challenged with a high concentration of Xccr (OD600 nm = 0.3) . Grapefruit (Citrus paradisi Macf. cv. Duncan), considered to be highly susceptible to the bacterium, showed the characteristic sequence of canker lesion development. Initially lesions appeared as water soaking, followed by the development of a raised corky form; each such lesion is a reservoir of new bacterial inoculum. Bacterial exudates were visible between 10 and 21 days post-inoculation. In contrast, PCD was observed in kumquat leaves in the form of a HR 3-5 days after inoculation with the canker-causing bacterium. Only necrotic lesions were observed and the bacterial population over time was shown to have an 'avirulent' incompatible growth pattern where bacterial multiplication ceased upon the development of necrosis [8, 25].
New tools have been developed in recent years through advances in genomics, proteomics, and bioinformatics that have particular utility for examining pathogen: host interaction complexities [22, 26–28]. The purpose of this study was to examine simultaneous changes in expression profiles for genes differentially expressed in the early stages (6-72 hpi) of citrus canker infection in kumquat, particularly those previously implicated in PCD-related responses such as HR.
Results and Discussion
In this study, identification of differentially expressed kumquat genes during its interaction with Xcc was pursued in an attempt to unravel the nature of the resistance mechanism(s) employed by the plant. Previously, kumquat suppression subtractive hybridization (SSH) cDNA libraries were constructed from Xcc-inoculated vs. mock inoculated leaves . Since SSH allows differential amplification of rare target sequences due to the elimination of more abundant house-keeping cDNA transcripts found in common from both samples, the technique has the potential of uncovering pertinent cDNA sequences. Subtraction was done in both directions, forward (inoculated-mock) and reverse (mock-inoculated) and the resulting cDNAs were subsequently sequenced. Preliminary screening macroarrays were used to confirm enrichment of the subtracted libraries with differentially expressed genes (data not shown).
Microarray experimental design
Functional annotation and an overview of global gene expression
The B2GO program  was used to assign GO (Gene Ontology) terms for hits obtained through eBLAST homology searches in NCBI. A general view of the similarity of the query set with the NCBI database, the distribution of the cut off for the e-value as well as the distribution of species with similar sequences are shown in Additional files 1, 2, and 3. The GO annotation score is considered to be more intuitive than regular blast e-values since GO annotation is carried out by applying an annotation rule (AR) on the ontology terms. Additionally, query sequence descriptions are obtained by applying a language processing algorithm that extracts informative names and avoids low-content terms such as "hypothetical protein" or "expressed protein". Using Blast2Go suite default parameters, 1042 probes were provided with GO annotations (Additional file 4). Approximately 25% of the transcripts on the array do not show similarity to proteins present in public databases. Some of these could represent exclusive genes of the citrus or kumquat lineages, but a fraction of these uncharacterized sequences may possibly represent low quality or 3'UTR sequences. Similar percentages of unknown sequences have been reported in other small-scale EST projects [30–32] and therefore this pattern can be considered characteristic of this approach. Since a citrus genome sequence is now available, future studies will have a wealth of citrus genomic sequence information that can be utilized to identify kumquat-specific as well as novel citrus genes involved in diverse defence mechanisms .
Kumquat transcriptional changes in response to Xcc infection
Functional categorization of transcripts underlines key elements in kumquat response to Xcc infection
Functional categorization of cDNAs identified from microarray analysis.
OXIDATIVE BURST/STRESS, APOPTOSIS
*Class III peroxidase
*lon protease homologue
*Ubiquitin-conjugating enzyme ubc7
Catalase (EC 22.214.171.124)CAT-2
Benzoic acid salicylic acid methyltransferase
Chlorophyll ab binding protein
Chloroplast photosystemII 22kda
*Pathogenesis-related protein 1a
Pathogenesis-related protein 4-1
*Class IV chitinase
*Zinc finger protein
Kumquat transcriptional changes in response to Xcc infection
ROS vs ROS scavenging
In order to maintain homeostasis and overcome the damaging effects of ROS (reactive oxygen species), a balance between SODs (superoxide dismutase) and the different H2O2-scavenging enzymes is considered to be critical in determining the levels of O2- and H2O2 in plant cells [41, 42]. Accordingly, there is a constant interplay between the antioxidant state and processes generating ROS. ROS are produced in chloroplasts, peroxisomes, and mitochondria in response to biotic as well as abiotic stresses [43, 44]. Accordingly, the expression of different enzymes that produce ROS were evidently stringently controlled and coordinated during the kumquat/Xcc interaction. For instance, while formate dehydrogenase (FDH; KLLRI2-G05), a mitochondrial NAD dependent enzyme, was 1.5 fold upregulated by 6 hpi, amine oxidase (KSLFI3-G05) that contributes to the synthesis of H2O2 and secondary metabolites was downregulated by 1.5 and 1.6 fold at 6 and 24 hpi respectively in response to Xcc challenge (Table 1). Concurrently, Xcc-inoculated kumquat plants overexpressed genes related to ROS scavenging to restrict damage to the inoculated parts of the plant, in this case the leaves. For instance, CuZnSOD (KLLF13-A03) expression in kumquats was increased ~1.5 fold at 6 hpi and was stabilized at 24 hpi and 72 dpi (Table 1). The same phenomenon was observed previously in tomato infected with Botrytis cinerea, a sign of increased ROS production by the host as part of the defence response to infection . Furthermore, while the expression of some of the genes linked to protease inhibitors and endopeptidase activities such as protease inhibitor homologue (KLLFI2-D02) was suppressed by the bacteria, other serine-type endopeptidase inhibitors such as an ATP-dependent ion protease (KSLFIV1-H05) was >2 fold up-regulated as early as 6 hpi subsequent to Xcc inoculation. In the same context, the redox coupling ascorbate-glutathione cycle, known to be responsible for peroxide detoxification , was repressed by 6 hpi in the kumquat dataset; examples include dehydroascorbate reductase (KSLRI1-F02) and glutathione peroxidase (KLLFII3-G07). Ascorbate and glutathione are non-enzymatic antioxidant molecules that have a role in other cycles, including those that synthesize and in some cases modulate flavonoids, alkaloids, phenolic compounds, α-tocopherol and carotenoids, all of which contribute in scavenging ROS . Dehydroascorbate reductase activity is indispensible when the ascorbate peroxidase (APX) levels are higher than normal under certain conditions to ensure preservation of the reduced form of ascorbate. Both proteins in addition to certain types of trypsin inhibitors might also catalyze a plant response . A similar study to investigate the Xanthomonas-grapefruit compatible interaction might present a platform to compare gene expression profiles of some genes of interest in both plants.
Accumulating evidence indicates that protein ubiquitination and degradation, last steps in protein turnover, are involved in plant defence responses. A number of recent studies have investigated a possible role of U-box E3 ubiquitin ligases in PTI (PAMPS-triggered immunity), ETI (effector-triggered immunity), as well as plant cell death and defence [49, 50]. In the present study, 6 ubiquitination pathway-related genes, for example ubiquitin-conjugating enzyme ucb7 (CSL1-A02), were isolated in the kumquat forward subtracted libraries; more investigation of their expression levels after infection will follow. Other induced genes that are involved in the proteolysis process are present in clusters A, B and C.
Genes involved in photosynthesis
A distinct down-regulation in the expression of ribulose1,5-bisphosphate carboxylase/oxygenase at 6 hpi, followed by an increase in expression that reaches maximum expression at 24 hpi, was observed in the microarray dataset (Rubisco small and large subunits; for example KLLFIII3-G09 and KSLRII2-F01) (Table 1). Rubisco, the most abundant protein in leaves, is the main source of energy production in plant cells. A decrease in photosynthesis was previously shown in Arabidopsis leaves as early as 3 h after challenge with the P. syringae avirulent strain, while after 48 h the rate of photosynthesis was lower with the virulent strain . Most of the photosynthetic machinery in challenged kumquat leaves was repressed at 6 hpi, including chlorophyll A/B binding protein (KLLFIII3-A06 and KLLFIII3-E08). Three photosynthesis-related genes were differentially down-regulated during the first 24 hours. In Pto-mediated resistance, 30 photosynthesis-related genes and 12 genes encoding chloroplast-associated proteins were suppressed . These results show that plants reduce photosynthetic potential to induce HR following pathogen attack. Further, Quirino et al.  suggested that HR and senescence are two programs that involve biochemical similarities as well as an overlap. The research reinforced the idea of a connection between defence response and senescence. Evidently, the down-regulation of genes involved in photosynthesis during the Xcc/kumquat interaction represents a cost for the plant fitness where energy resources were redirected to defence response.
Cell wall remodelling
Xcc inoculation of kumquat was followed by the down-regulation of various genes related to cell wall remodeling and rapid expansion such as endoglucanases. The expression level of a kumquat homologue of this wall loosening protein (KLLFI2-C10) was insignificant. On the other hand, genes related to cell wall reorganization, for example xyloglucan endotransglycosylase/hydrolase (XET,-- an enzyme involved in cell wall elongation and restructuring), were significantly up-regulated by 24 hpi (KSLFIII1-H08). In 'Pera' sweet orange, a major difference in the response to inoculation of the two bacterial strains was that Xcc strongly upregulated several cell wall remodelling enzymes, while Xaa upregulated genes related to endoglucanase inhitors and lignin biosynthesis. A phenomenon that we observed in kumquat plants is the development of a few minute necrotic flecks on the leaves when inoculated with low concentrations of the bacterium (Xcc). Neither leaf abscission nor water soaked lesions were observed on the leaves later under our conditions. It is also worth mentioning that although Cernadas et al. used a relatively high concentration of Xaa (OD600 nm = 0.6, ~double what we used for Xcc with kumquat) only pustles were recorded in sweet orange inoculated with Xaa that were not followed by necrosis . Using light microscopy, we have previously shown mesophyl collapse in kumquat leaves which was followed by leaf abscission 72 hrs post inoculation with Xcc. Alternatively, grapefruit mesophyl cells from inoculated leaves showed enlargement (hypertrophy) and division (hyperplasia) followed by raised circular lesions that became raised and developed into white or yellow spongy pustules. These pustules then darkened and thickened into brown corky canker lesions . Pustule formation and hypertrophy were linked previously to the PthA effector in Nicotiana benthamiana . Alternatively, accumulation of the tomato XTH (xyloglucan endotransglucosylase/hydrolase LeXTH1) protein 6 hours after attachment of the parasite has provided evidence for a role of XTH in defence reactions associated with the incompatible tomato-Cuscuta interaction as was presented in Albert et al. .
Resistance genes and related proteins
qRT-PCR analysis of genes expressed in response to Xcc inoculation a (5 × 108cfu/ml) concentration of the Miami strain X04-59.
Basic leucine zipper transcription factor
Class I chitinase (CHI1)
A disease resistance leucine-rich repeat protein
Receptor-like serine threonine kinase
A putative beta-galactosidase BG1
protein kinase 3
Interestingly, the expression levels of the kumquat PR1 (Pathogenesis -related gene1) gene homologue, normally a marker of salicylic acid-induced systemic acquired resistance (SAR) that is usually up-regulated after pathogen infections, was lower in the infected samples as early as 6 hpi compared to control mock inoculated samples (KLLFII2-A04) (according to the microarray results and non-published data). The region upstream of the PR1 promoter, W-box sequences, was shown previously to act as a negative cis-acting element in the expression of defence related genes . This implies a different basal defence response and SAR regulation mechanism from that of Arabidopsis and other dicotyledonous plants in kumquat after Xcc infection. PR genes that were identified in citrus vary in their responses to different pathogens as shown in this study and others where variation in expression levels of different members of the PR gene family was dependent on the nature of different elicitors [60–62]. The biological activity of a large majority of PR genes in plants during biotic stress is yet to be revealed. More interestingly, according to the microarray results, expression levels of a number of homologues to other defence related genes such as the NDR1 gene (KLLFII2-E03) were also repressed 24 hpi after Xcc infection.
LRR proteins are known to be a part of the early signal transduction cascade involved in the recognition of pathogen Avr products . Sequences for a number of homologues known to be part of different hormonal defence pathways (for instance transcription factors, receptor like and receptor-like kinases) were found to be differentially expressed in kumquat after Xcc inoculation.
Key molecular features of kumquat PCD
KSLFI3-C10 is homologous to hsr203J , a carboxylesterase (CXE) gene implicated previously in the incompatible interactions between tobacco and the bacterial pathogen Ralstonia solanacearum. Its promoter is highly, rapidly, and specifically activated in response to HR inducing bacterial inoculation, does not respond to various stress conditions, and is strongly dependent on hrp (hypersensitive response and pathogenicity) genes of the pathogenic bacterium . It has been proposed that its expression should be a useful marker for programmed cell death occurring in response to diverse pathogens.
KLLRI2-B05 shares homology with DND1 (DEFENCE NO DEATH 1), which encodes a cyclic nucleotide-gated ion channel that allows passage of Ca2+, K+ and other cations. The Arabidopsis thaliana dnd1 mutant failed to produce HR cell death in response to an avirulent pathogen infection .
During programmed cell death or apoptosis cytochrome c is released to the cytoplasm from the intermembrane space of the mitochondrion [66, 67]. Once in the cytoplasm, it activates caspases (cysteine aspartate-specific proteases), killer proteins that dismantle the cell . Two key proteins known to be core components of the apoptic machinery in animals, caspase and Bax-Inhibitor1 gene homologues, were identified in our dataset. According to our data, a homologue that has caspase activity KLLRI2_A12 was slightly up-regulated by 6 hpi in kumquat challenged leaves. Bax is a member of the Bcl2 family that plays a regulatory role preventing apoptosis by inhibiting adapters needed for the activation of caspases . A kumquat homologue of the Bax-Inhibitor1 gene (KLLFIII2-E06) was shown to be slightly up regulated 6 hpi in response to Xcc challenge as was previously shown with Arabidopsis thaliana Bax Inhibitor-1 (AtBI-1), isolated during a differential screen of plants challenged with the phytopathogen Pseudomonas syringae . In the same context, Bax inhibitor has been shown to trigger cytochrome c release from mitochondria both in vitro and in vivo in animals.
Endopeptidase inhibitors are often part of an inducible, jasmonic acid associated defence pathway that accumulates upon wounding, pathogen, or herbivore damage in leaves . The antagonistic interaction between proteases and endopeptidase inhibitors is considered to be a cell death control mechanism . Li et al., 2008 demonstrated that a serine protease (Kunitz trypsin) inhibitor (KTI1) of Arabidopsis is involved in modulating PCD in plant-pathogen interactions . RNAi silencing of the AtKTI1 gene resulted in enhanced lesion development after infiltration of leaf tissue with the PCD-eliciting fungal toxin fumonisin B1 (FB1) or the avirulent bacterial pathogen Pseudomonas syringae pv tomato DC3000 carrying avrB (Pst avrB). Trypsin inhibitor (KSLFIII1-H12 and KLLFIII3-F03) and a miraculin serine type endopeptidase inhibitor (FI2-A05) sequences were found in original early subtraction libraries representing transcripts expressed during early infection (30 min.pi-24 hpi). While KLLFIII3-F03 (trypsin homologue) gene expression was significantly suppressed 6 hpi, the expression of KSLIII1-H12 was slightly suppressed and then 2 fold upregulated 24 hpi. Further analysis should be done to study the difference between the mechanism(s) of action of these two genes.
Suppression of defence responses
A very evident down-regulation of a considerable number of genes was recorded by 6 hpi which may be caused by defence suppression imposed by Xcc effectors (Clusters A and C; Figure 6A and 6B). It has been shown previously that Xcc exploits the Type III secretion system (T3SS) to inject different effector proteins into citrus plants in order to avoid host recognition and subsequently MAMPS/PAMP-triggered immunity. The bacterial effector proteins suppress plant defences including basal defence, gene-for-gene resistance, and nonhost resistance. There was no accumulation of any SAR gene transcripts including PR1, a marker for enhanced defence; in addition some other key elements in the SA defence pathway were suppressed. On the other hand, the S-adenosyl-l-methionine:benzoic acid salicylic acid carboxyl methyltransferase gene (KLLRI2-C03) was at least 1-fold upregulated in kumquat leaves by 6 hpi in response to Xcc inoculation; the gene is known to play a role in plant defence responses . In addition, the SA-binding protein 2 (SABP2 KLLRI2-G01), a lipase protein that belongs to the hydrolase super family, was found to be up-regulated at 6 hpi by at least 2 fold; the gene was previously found to be required for the plant immune response in tobacco .
Realtime Quantitative Polymerase Chain Reaction Validation
Validation of the presented microarray dataset was carried out using TaqMan gene expression assay for a number of homologues on the array. Genes that were implicated in plant defence including basic leucine zipper transcription factor (KLLFII2-C05), a putative chitinase protein (CHI1)(KLLRI2-D05), a putative disease resistance leucine rich protein (CSL1-D05), a receptor like protein kinase (KLLRI2-H10), a beta galactosidase like protein (KLLFII1-B11) and a putative mitogen-activated protein (CSL2-A02) were selected for validation using q RT-PCR. As summarized in Table 2 and Additional file 10, the qRT-PCR data correlated with the microarray results confirming the up-or down-regulation of all analyzed genes although as expected the qRT-PCR was more sensitive.
In this study, a F. margarita custom microarray representing 1024 unigenes was used to study the response to inoculation with X. axonopodis pv. citri. A very distinct though delayed HR was observed in Xcc-inoculated kumquat plants where initially the bacterium grew exponentially, followed by a sudden leaf tissue collapse (necrosis with no canker lesions) 2-5 days after inoculation . A comparable delayed HR was observed in tomato resistance response to race T3 mediated by AvrXv3 effector and RxvT3 R protein . The current kumquat analysis allowed simultaneous investigation of the expression of more than one group of genes known to be linked to more than one biological process and cellular compartment in relation to the HR caused by Xcc infection. A large number of genes were found to be differentially expressed after infection. Most of the genes involved in defence mechanisms in kumquat appear to be associated with the phenomena that precede the HR including oxidative burst, protein degradation, and regulation of photosynthesis as well as the production of ROS that is associated with the oxidative burst. One very distinct observation was that some of the defence genes such as PR1 and NDR1 were down-regulated in kumquats in response to Xcc inoculation as early as 6 hpi, a phenomenon currently under further examination. What clearly appears to be a resistant response and a drastic decrease in the bacterial population, in addition to the activation of genes involved in ROS production as well as and programmed cell death, seems to be a common mechanism that is pursued by more than one citrus bacterial pathogen with no associated-resistance genes yet identified [24, 36]. Future work will compare differences in gene response in both resistant and susceptible citrus types.
Plant material and inoculation with bacteria
Fortunella margarita (Lour.) Swingle (Nagami kumquat) plants were used in all of the experiments described in this study. Plants were approximately 2 years old at the time of the experiment and were maintained in the quarantine greenhouse facility at the Division of Plant Industry, Florida Department of Agriculture (Gainesville, FL, USA) under controlled conditions. Leaves from a set of six kumquat plants were infiltrated with bacterial cultures according to Lund et al. . The bacterial strain used was Xanthomonas citri subsp. citri A; Miami X04-59 (Xcc). The inoculum was adjusted to 5 × 108 cfu/ml. A similar set of plants was mock-inoculated using sterile tap water as controls. Leaves from the two sets of plants were used in subsequent experiments.
The kumquat microarray chip was developed and printed at the University of Florida (Gainesville, FL, USA). The array included ESTs chosen from 4 previously constructed Nagami kumquat forward and reverse leaf subtraction cDNA libraries. The cDNA libraries were constructed using RNA extracted from leaf tissue collected at different intervals post inoculation with Xcc (see below) and pooled into early and late library sets to capture a wide spectrum of differentially expressed transcripts . Random DNA sequencing was performed from 5'and 3' ends of randomly selected clones using universal primers, generating sequence information from 2788 and 1655 clones from the early and late leaf subtraction libraries respectively.
The initial dataset was reduced to a total of 2304 transcripts that were selected according to sequence alignment similarities with proteins in the Genbank database. Sequences were selected based on quality and length. The dataset included 2254 kumquat ESTs comprising 738 contigs and 1516 singletons, in addition to 50 cDNA control elements. Each probe was printed in 3 locations on the array using the Omnigrid Microarrayer (Gene Machines, San Carlos, CA, USA) so that all clones had 3 technical replicates on each slide, generating a total of 6912 spots. Post-printing slide processing was performed as described in Heller et al. (1997)  with some modifications. In brief, a combination of sequential baking and UV crosslinking was implemented where slides were baked for 80 min at 80°C in a drying oven without vacuum. The slides were then washed twice in 0.1% SDS for 5 minutes each to remove any unbound DNA.
To identify genes that are considered to be differentially expressed in kumquat, a time-course experiment was designed utilizing the kumquat/Xcc pathosystem. Six independent Xcc or mock-inoculated kumquat plants were used; each plant was considered an independent biological replicate. Since citrus canker is a non-systemic disease, 6-10 leaves per treated kumquat plant were independently infiltrated using 5 × 108 cfu/ml Xcc. All RNA samples isolated from healthy mock- or Xcc inoculated (infected) leaves were processed independently. It is unlikely that differential gene expression observed was caused by the pressure infiltration inoculation method used, since this factor was normalized by treating the mock inoculated plants in the exact same way as infected plants.
Individual leaves were harvested from the inoculated and mock-inoculated plants at specific time-points post-inoculation (pi) according to designated conditions for each experiment; there were 3 time points for the microarray experiment and 5 for subsequent real-time PCR assays. For the microarray experiment, the three time points (6 hpi, 24 hpi, and 72 hpi) were chosen based upon the internal bacterial populations previously detected at these times following inoculation and the knowledge that kumquat leaves abscised 3-5 days after inoculation. In addition, previous experiments revealed that there were some transcripts differentially expressed as early as 30 min post-inoculation with Xcc. Finally, the RNA yield and the abundance of cellular transcripts decreased as the leaves approached total PCD, as has been shown previously by others . Time points of 0 and 120 hpi were added for the real time PCR assays. The healthy mock-inoculated and Xcc-inoculated leaf samples were immediately frozen in liquid nitrogen. For each respective time point, total RNA was extracted using RNeasy columns (QIAGEN, Valencia, CA, U.S.A.) according to the manufacturer's protocol. RNA purity, concentration and quality were assessed using a spectrophotometer and a BioAnalyzer 2100 (Agilent Technologies, Palo Alto, CA).
Fluorescent probe, hybridization, and scanning
Prior to slide hybridization with probe, slides were pre-hybridized in a solution containing 5× SSC, 0.1% SDS and 1% bovine serum albumin at 42°C for 45 min to eliminate nonspecific binding of the probe to the slide. Slides were washed using MilliQ RNase free water, then isopropanol, once each, and air-dried. Slides were maintained at the hybridization temperature until loaded with probe. cDNA labeling was performed using the Genisphere (Hatfield, PA, USA) 3DNA Array50® Expression Array Detection Kit according to the manufacturer's protocol for total RNA. For each time-point 125 μg of DNA-free total RNA isolated from an independent plant (biological replicate) per slide was reverse transcribed for each of the mock inoculated and the infected leaf samples separately using Ambion reverse transcriptase (Ambion; LaJolla, CA) in the presence of Genisphere dT primers. Two-step hybridization was performed as follows. The first hybridization, carried out at 48°C overnight, contained 10 μl of the concentrated cDNA (heat denatured probe) made using either the Cy5-RT primer capture sequence or the Cy3-RT primer capture sequence, in Genisphere 2× formamide-based hybridization buffer. Three successive post hybridization washes were performed, first in 2X SSC, 0.2% SDS at 55°C for 10 min, then 2× SSC for 10 min and finally in 0.2× SSC for 10 min at room temperature. The second hybridization (for addition of the dye) was performed using 2.5 μl of either Cy3 or Cy5 dendrimer, 2 μl of high-end differential buffer and 58 μl of hybridization buffer. For each time point, three mock-inoculated samples from three individual plants were labeled with Cy-3 and the Xcc infected samples were labeled with Cy-5 and both were hybridized to the same array. Post-hybridization washes were conducted as performed earlier following the primary hybridization, with the addition of 0.1 ml dithiothreitol (DTT) into the first and second wash solutions to reduce oxidation of fluorescent dyes.
cDNA microarray setup and quality control
Control measures such as the detection sensitivity level were determined using internal control probes and non-specific control elements. Human genomic DNA, the green fluorescent protein gene, and the lambda control template DNA fragment were included as negative controls. Additionally, cDNAs previously implicated in pathogen defence such as PR1 and NPR1 from Arabidopsis, NDR1 from citrus, were printed 3 times on the array to test the ability of the microarray method to detect changes in gene expression. These were considered to be specific positive controls. In addition, the microarray ratio for each gene analyzed was normalized against the microarray ratio obtained for 18 S.
Transcriptome data analysis
Agilent's Feature Extraction Software (Agilent Corp., Palo Alto, CA) was used to analyze the microarray data. Data were uploaded into the statistical platform R  for statistical analysis and the Limma package was used for pre-processing. Data were Lowes transformed followed by scaling between arrays . Two array slides were chosen for each time-point experiment based on the consistency of the signal across the replicates. The fold difference in expression was computed as: 2٨average ratio (2 to the power of the average ratio). cDNAs with an average ratio of 1.0 or higher were considered differentially expressed, which represents a 1.5 fold or higher difference in expression. Statistical analysis was performed using two different approaches. Time-dependent gene expression changes were analyzed by the maSigPro methodology . Data were subsequently subjected to ASCA-gene analysis that combines ANOVA and multivariate methods to identify main and secondary patterns of gene expression associated with different experimental factors . Statistical analysis identified a number of selected genes that were further grouped into clusters.
Functional information about the ESTs represented in the array was obtained by Blast2GO analysis using default parameters . Blast2GO uses Blast and an elaborated annotation algorithm to assign Gene Ontology (GO), Enzyme Code and InterPro functional labels to a set of uncharacterized sequences . The functional characterization of these clusters was done by applying the Functional Enrichment (FE) method included in Blast2GO which implements the Gossip algorithm. FE methods assess which functional categories are over-represented within a group of genes in relation to a broader list, in this case the whole kumquat array. Finally, the major induced transcriptional changes considered functional classes as a whole were studied with the PCA-maSigFun method . This method combines Principal Component Analysis and maSigPro to characterize the "expression profiles" associated with cellular functions. Sequence data from this work have been deposited in the NCBI Genbank database libraries (http://www.ncbi.nlm.nih.gov/GenBank/index.html), using the BankIt dbEST database, and accession numbers were obtained. [Genbank: GW687757to GW690680]. (See Additional File 11).
Quantitative Real-time Quantitative PCR
Kumquat leaves were infiltrated with Xcc (5 × 108 cfu per milliliter), then total RNA was isolated from inoculated leaves 0,6,24,72, 120 hpi for both the microarray and the quantitative real-time PCR experiment as previously stated in the 'experimental plan'.
Total RNA was isolated separately for each respective time point using the TRIzol reagent (Invitrogen, Carlsbad,CA, U.S.A.) according to manufacturer's instructions. RNA samples were further purified using the RNeasy Plant Mini Kit (Qiagen, Valencia, CA) including DNaseI TURBO DNase (Ambion, Austin, TX, U.S.A.) treatment according to the manufacturer's instructions. RNA quality and quantity was then assessed using microspectrophotometry (Nanodrop Technologies, Inc., Wilmington, DE). cDNA was synthesized from 1 μg of total RNA using Applied Biosystems (High Capacity cDNA Reverse Transcription Kit PN 4368813, 4374966). A TaqMan gene expression assay was then used to validate the transcript accumulation levels of a specific subset of genes from the kumquat microarrays. Reactions were performed in the ABI Prism7900 HT sequence detector (Applied Biosystems, Foster City, CA, U.S.A.). Primers for qRT-PCR were designed using the Primer Express 3.0 software (Applied Biosystems), and data were normalized using a Taqman ribosomal RNA control, in addition to the kumquat 18S ribosomal gene that served as an internal control where each real-time PCR reaction was done in parallel with the 18S primers. For internal controls, a number of genes, for example actin showed inconsistencies. The 18s surprisingly showed coherency throughout the interaction, and this was noticed during the microarrays and was subsequently confirmed with the Realtime PCR (RT-qPCR) 18S expression curve. qRT-PCR was carried out at 50°C for 2 minutes, 95°C for 10 minutes, followed by 40 cycles at 95°C for 15 seconds and 60°C for 1 minute.
Acknowledgements and Funding
This research was supported in part by the Citrus Research and Education Center (CREC), Dr. Frederick Gmitter, Jr. and by a grant from the Florida Citrus Production Research Advisory Council, Florida Department of Citrus (FDACS Contract Number 67), "Peter McClure, Chairman".
We gratefully acknowledge the generous financial and intellectual support from Florida citrus growers, Mrs. Hariot and Mr.Barnette E. "Barney" Greene Jr. that helped complete this work.
We would also like to thank Dr. Henry Baker, Professor and Interim Chair and Dr. Hassan Badrane (Medicine: Molecular Genetics and Microbiology Dept at UF) for providing help with microarrays equipments. We would also like to thank Dr. Jeffrey Jones and Dr. Vicente Febres for their objective critical review of the manuscript.
- Gottwald TR, Graham JH, Schubert TS: Citrus canker: The pathogen and its impact. Online Plant Health Progress. 2002.
- Głowacki S, Macioszek V, Kononowicz A: R proteins as fundamentals of plant innate immunity. Cellular & Molecular Biology Letters. 2010, 1-24.
- Jones JD, Dangl JL: The plant immune system. Nature. 2006, 444 (7117): 323-329. 10.1038/nature05286.PubMedView Article
- Fujikawa T, Ishihara H, Leach JE, Tuyumu S: Suppression of defense response in plants by avrBs3/pthA gene family of Xanthomonas spp. MolPlant-Microbe Interact. 2006, 19: 342-349. 10.1094/MPMI-19-0342.View Article
- Yang B, Sugio A, White FF: Avoidance of host recognition by alterations in the repetitive and C-terminal regions of AvrXa7, a type III effector of Xanthomonas oryzae pv. oryzae. Mol Plant Microbe Interact. 2005, 18 (2): 142-149. 10.1094/MPMI-18-0142.PubMedView Article
- Alfano JR, Collmer A: The type III (Hrp) secretion pathway of plant pathogenic bacteria: trafficking harpins, Avr proteins, and death. J Bacteriol. 1997, 179 (18): 5655-5662.PubMedPubMed Central
- Sanabria N, Goring D, Nurnberger T, Dubery I: Self/nonself perception and recognition mechanisms in plants: a comparison of self-incompatibility and innate immunity. New Phytol. 2008, 178 (3): 503-514. 10.1111/j.1469-8137.2008.02403.x.PubMedView Article
- Gottig N, Garavaglia BS, Daurelio LD, Valentine A, Gehring C, Orellano EG, Ottado J: Xanthomonas axonopodis pv. citri uses a plant natriuretic peptide-like protein to modify host homeostasis. Proc Natl Acad Sci USA. 2008, 105 (47): 18631-18636. 10.1073/pnas.0810107105.PubMedPubMed CentralView Article
- Bestwick CS, Brown IR, Bennett MH, Mansfield JW: Localization of hydrogen peroxide accumulation during the hypersensitive reaction of lettuce cells to Pseudomonas syringae pv phaseolicola. Plant Cell. 1997, 9 (2): 209-221.PubMedPubMed CentralView Article
- Ali R, Ma W, Lemtiri-Chlieh F, Tsaltas D, Leng Q, von Bodman S, Berkowitz GA: Death don't have no mercy and neither does calcium: Arabidopsis CYCLIC NUCLEOTIDE GATED CHANNEL2 and innate immunity. Plant Cell. 2007, 19 (3): 1081-1095. 10.1105/tpc.106.045096.PubMedPubMed CentralView Article
- Shiotani H, Fujikawa T, Ishihara H, Tsuyumu S, Ozaki K: A pthA homolog from Xanthomonas axonopodis pv. citri responsible for host-specific suppression of virulence. J Bacteriol. 2007, 189 (8): 3271-3279. 10.1128/JB.01790-06.PubMedPubMed CentralView Article
- Brunings AM, Gabriel DW: Xanthomonas citri: breaking the surface. Molecular Plant Pathology. 2003, 4 (3): 141-157. 10.1046/j.1364-3703.2003.00163.x.PubMedView Article
- Romer P, Hahn S, Jordan T, Strauss T, Bonas U, Lahaye T: Plant pathogen recognition mediated by promoter activation of the pepper Bs3 resistance gene. Science. 2007, 318 (5850): 645-648. 10.1126/science.1144958.PubMedView Article
- Yang B, White FF: Diverse members of the AvrBs3/PthA family of type III effectors are major virulence determinants in bacterial blight disease of rice. Mol Plant Microbe Interact. 2004, 17 (11): 1192-1200. 10.1094/MPMI.2004.17.11.1192.PubMedView Article
- Duan YP, A C, G Z, G E, Gabriel DW: Expression of a single, host-specific, bacterial pathogenicity gene in plant cells elicits division, enlargement, and cell death. Molecular Plant-Microbe Interactions. 1999, 12 (6): 556-560. 10.1094/MPMI.19126.96.36.1996.View Article
- Tasic L, Borin PF, Khater LC, Ramos CH: Cloning and characterization of three hypothetical secretion chaperone proteins from Xanthomonas axonopodis pv. citri. Protein Expr Purif. 2007, 53 (2): 363-369. 10.1016/j.pep.2007.01.011.PubMedView Article
- White FF, Potnis N, Jones JB, Koebnik R: The type III effectors of Xanthomonas. Mol Plant Pathol. 2009, 10 (6): 749-766. 10.1111/j.1364-3703.2009.00590.x.PubMedView Article
- Schenk PM, Kazan K, Manners JM, Anderson JP, Simpson RS, Wilson IW, Somerville SC, Maclean DJ: Systemic gene expression in Arabidopsis during an incompatible interaction with Alternaria brassicicola. Plant Physiol. 2003, 132 (2): 999-1010. 10.1104/pp.103.021683.PubMedPubMed CentralView Article
- Spoel SH, Dong X: Making sense of hormone crosstalk during plant immune responses. Cell Host Microbe. 2008, 3 (6): 348-351. 10.1016/j.chom.2008.05.009.PubMedView Article
- Lorenzo O, Piqueras R, Sanchez-Serrano JJ, Solano R: ETHYLENE RESPONSE FACTOR1 integrates signals from ethylene and jasmonate pathways in plant defense. Plant Cell. 2003, 15 (1): 165-178. 10.1105/tpc.007468.PubMedPubMed CentralView Article
- Nemhauser JL, Hong F, Chory J: Different plant hormones regulate similar processes through largely nonoverlapping transcriptional responses. Cell. 2006, 126 (3): 467-475. 10.1016/j.cell.2006.05.050.PubMedView Article
- Glazebrook J: Use of microarray analysis to dissect the plant defense response. Methods Mol Biol. 2007, 354: 121-130.PubMed
- Leon-Reyes A, Spoel SH, De Lange ES, Abe H, Kobayashi M, Tsuda S, Millenaar FF, Welschen RA, Ritsema T, Pieterse CM: Ethylene Modulates the Role of NPR1 in Cross-Talk Between Salicylate and Jasmonate Signaling. Plant Physiol. 2009.
- Khalaf AA, Moore GA, Jones JB, Gmitter FG: New insights into the resistance of Nagami kumquat to canker disease. Physiological and Molecular Plant Pathology. 2007, 71: 240-250. 10.1016/j.pmpp.2008.03.001.View Article
- Alfano JR, Collmer A: Bacterial pathogens in plants: Life up against the wall. Plant Cell. 1996, 8 (10): 1683-1698.PubMedPubMed CentralView Article
- Alba R, Fei Z, Payton P, Liu Y, Moore SL, Debbie P, Cohn J, D'Ascenzo M, Gordon JS, Rose JK, et al: ESTs, cDNA microarrays, and gene expression profiling: tools for dissecting plant physiology and development. Plant J. 2004, 39 (5): 697-714. 10.1111/j.1365-313X.2004.02178.x.PubMedView Article
- Gotz S, Garcia-Gomez JM, Terol J, Williams TD, Nagaraj SH, Nueda MJ, Robles M, Talon M, Dopazo J, Conesa A: High-throughput functional annotation and data mining with the Blast2GO suite. Nucleic Acids Res. 2008, 36 (10): 3420-3435. 10.1093/nar/gkn176.PubMedPubMed CentralView Article
- Talon M, Gmitter FG: Citrus genomics. Int J Plant Genomics. 2008, 2008: 528361.PubMedPubMed CentralView Article
- Conesa A, Gotz S, Garcia-Gomez JM, Terol J, Talon M, Robles M: Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005, 21 (18): 3674-3676. 10.1093/bioinformatics/bti610.PubMedView Article
- Cohn JR, Martin GB: Pseudomonas syringae pv. tomato type III effectors AvrPto and AvrPtoB promote ethylene-dependent cell death in tomato. Plant J. 2005, 44 (1): 139-154. 10.1111/j.1365-313X.2005.02516.x.PubMedView Article
- Desikan R, S AH-M, Hancock JT, Neill SJ: Regulation of the Arabidopsis transcriptome by oxidative stress. Plant Physiol. 2001, 127 (1): 159-172. 10.1104/pp.127.1.159.PubMedPubMed CentralView Article
- Gang DR, Wang J, Dudareva N, Nam KH, Simon JE, Lewinsohn E, Pichersky E: An investigation of the storage and biosynthesis of phenylpropenes in sweet basil. Plant Physiol. 2001, 125 (2): 539-555. 10.1104/pp.125.2.539.PubMedPubMed CentralView Article
- Conesa A, Gotz S: Blast2GO: A comprehensive suite for functional analysis in plant genomics. Int J Plant Genomics. 2008, 2008: 619832.PubMedPubMed CentralView Article
- Nueda MJ, Conesa A, Westerhuis JA, Hoefsloot HC, Smilde AK, Talon M, Ferrer A: Discovering gene expression patterns in time course microarray experiments by ANOVA-SCA. Bioinformatics. 2007, 23 (14): 1792-1800. 10.1093/bioinformatics/btm251.PubMedView Article
- Kay S, Hahn S, Marois E, Hause G, Bonas U: A Bacterial Effector Acts as a Plant Transcription Factor and Induces a Cell Size Regulator. Science. 2007, 318 (5850): 648-651. 10.1126/science.1144956.PubMedView Article
- Cernadas RA, Camillo LR, Benedetti CE: Transcriptional analysis of the sweet orange interaction with the citrus canker pathogens Xanthomonas axonopodis pv. citri and Xanthomonas axonopodis pv. aurantifolii. Mol Plant Pathol. 2008, 9 (5): 609-631. 10.1111/j.1364-3703.2008.00486.x.PubMedView Article
- Albrecht U, Bowman KD: Transcriptional response of susceptible and tolerant citrus to infection with Candidatus Liberibacter asiaticus. Plant Science. 2011.
- Ross J, Li Y, Lim E, Bowles DJ: Higher plant glycosyltransferases. Genome Biol. 2001, 2 (2): Reviews 3004.
- Dixon RA: Natural products and plant disease resistance. Nature. 2001, 411 (6839): 843-847. 10.1038/35081178.PubMedView Article
- Mugford ST, Qi X, Bakht S, Hill L, Wegel E, Hughes RK, Papadopoulou K, Melton R, Philo M, Sainsbury F, et al: A serine carboxypeptidase-like acyltransferase is required for synthesis of antimicrobial compounds and disease resistance in oats. Plant Cell. 2009, 21 (8): 2473-2484. 10.1105/tpc.109.065870.PubMedPubMed CentralView Article
- Nanda AK, Andrio E, Marino D, Pauly N, Dunand C: Reactive Oxygen Species during plant-microorganism early interactions. Journal of Integrative Plant Biology. 2010, 52 (2): 195-204. 10.1111/j.1744-7909.2010.00933.x.PubMedView Article
- Torres MA: ROS in biotic interactions. Physiologia Plantarum. 2010, 138 (4): 414-429. 10.1111/j.1399-3054.2009.01326.x.PubMedView Article
- de Pinto MC, Paradiso A, Leonetti P, De Gara L: Hydrogen peroxide, nitric oxide and cytosolic ascorbate peroxidase at the crossroad between defence and cell death. Plant J. 2006, 48 (5): 784-795. 10.1111/j.1365-313X.2006.02919.x.PubMedView Article
- Vranova E, Inze D, Van Breusegem F: Signal transduction during oxidative stress. J Exp Bot. 2002, 53 (372): 1227-1236. 10.1093/jexbot/53.372.1227.PubMedView Article
- Kuzniak E, Sklodowska M: The effect of Botrytis cinerea infection on the antioxidant profile of mitochondria from tomato leaves. J Exp Bot. 2004, 55 (397): 605-612. 10.1093/jxb/erh076.PubMedView Article
- Halliwell B, Foyer CH: Ascorbic acid, metal ions and the superoxide radical. Biochem J. 1976, 155 (3): 697-700.PubMedPubMed CentralView Article
- Foyer CH, Noctor G: Redox homeostasis and antioxidant signaling: a metabolic interface between stress perception and physiological responses. Plant Cell. 2005, 17 (7): 1866-1875. 10.1105/tpc.105.033589.PubMedPubMed CentralView Article
- Trumper S, Follmann H, Haberlein I: A novel-dehydroascorbate reductase from spinach chloroplasts homologous to plant trypsin inhibitor. FEBS Lett. 1994, 352 (2): 159-162. 10.1016/0014-5793(94)00947-3.PubMedView Article
- Yang C-W, Gonzalez-Lamothe R, Ewan RA, Rowland O, Yoshioka H, Shenton M, Ye H, O'Donnell E, Jones JDG, Sadanandom A: The E3 ubiquitin ligase activity of Arabidopsis PLANT U-BOX17 and Its functional tobacco Homolog ACRE276 are required for cell death and defense. Plant Cell. 2006, 18 (4): 1084-1098. 10.1105/tpc.105.039198.PubMedPubMed CentralView Article
- Trujillo M, Ichimura K, Casais C, Shirasu K: Negative regulation of PAMP-triggered immunity by an E3 ubiquitin ligase triplet in Arabidopsis. Curr Biol. 2008, 18 (18): 1396-1401. 10.1016/j.cub.2008.07.085.PubMedView Article
- Bonfig KB, Schreiber U, Gabler A, Roitsch T, Berger S: Infection with virulent and avirulent P. syringae strains differentially affects photosynthesis and sink metabolism in Arabidopsis leaves. Planta. 2006, 225 (1): 1-12. 10.1007/s00425-006-0303-3.PubMedView Article
- Mysore KS, D'Ascenzo MD, He X, Martin GB: Overexpression of the disease resistance gene Pto in tomato induces gene expression changes similar to immune responses in human and fruitfly. Plant Physiol. 2003, 132 (4): 1901-1912. 10.1104/pp.103.022731.PubMedPubMed CentralView Article
- Quirino BF, Noh YS, Himelblau E, Amasino RM: Molecular aspects of leaf senescence. Trends Plant Sci. 2000, 5 (7): 278-282. 10.1016/S1360-1385(00)01655-1.PubMedView Article
- Albert M, Werner M, Proksch P, Fry SC, Kaldenhoff R: The cell wall-modifying xyloglucan endotransglycosylase/hydrolase LeXTH1 is expressed during the defence reaction of tomato against the plant parasite Cuscuta reflexa. Plant Biol (Stuttg). 2004, 6 (4): 402-407. 10.1055/s-2004-817959.View Article
- Zhu-Salzman K, Salzman RA, Ahn JE, Koiwa H: Transcriptional regulation of sorghum defense determinants against a phloem-feeding aphid. Plant Physiol. 2004, 134 (1): 420-431. 10.1104/pp.103.028324.PubMedPubMed CentralView Article
- Lee S, Kim SY, Chung E, Joung YH, Pai HS, Hur CG, Choi D: EST and microarray analyses of pathogen-responsive genes in hot pepper (Capsicum annuum L.) non-host resistance against soybean pustule pathogen (Xanthomonas axonopodis pv. glycines). Funct Integr Genomics. 2004, 4 (3): 196-205.PubMedView Article
- Brunner F, Stintzi A, Fritig B, Legrand M: Substrate specificities of tobacco chitinases. Plant J. 1998, 14 (2): 225-234. 10.1046/j.1365-313X.1998.00116.x.PubMedView Article
- Kasprzewska A: Plant chitinases--regulation and function. Cell Mol Biol Lett. 2003, 8 (3): 809-824.PubMed
- Lebel E, Heifetz P, Thorne L, Uknes S, Ryals J, Ward E: Functional analysis of regulatory sequences controlling PR-1 gene expression in Arabidopsis. Plant J. 1998, 16 (2): 223-233. 10.1046/j.1365-313x.1998.00288.x.PubMedView Article
- Mitsuhara I, Iwai T, Seo S, Yanagawa Y, Kawahigasi H, Hirose S, Ohkawa Y, Ohashi Y: Characteristic expression of twelve rice PR1 family genes in response to pathogen infection, wounding, and defense-related signal compounds (121/180). Molecular genetics and genomics: MGG. 2008, 279 (4): 415-427. 10.1007/s00438-008-0322-9.PubMedPubMed CentralView Article
- Campos MA, Rosa DD, Teixeira JÉC, Targon MLPN, Souza AA, Paiva LV, Stach-Machado DR, Machado MA: PR gene families of citrus: their organ specific-biotic and abiotic inducible expression profiles based on ESTs approach. Genetics and Molecular Biology. 2007, 30: 917-930. 10.1590/S1415-47572007000500020.View Article
- Li C-W, Su R-C, Cheng C-P, You S-J, Hsieh T-H, Chao T-C, Chan M-T: Tomato RAV Transcription Factor Is a Pivotal Modulator Involved in the AP2/EREBP-Mediated Defense Pathway. Plant Physiology. 2011, 156 (1): 213-227. 10.1104/pp.111.174268.PubMedPubMed CentralView Article
- Dodds PN, Lawrence GJ, Ellis JG: Six amino acid changes confined to the leucine-rich repeat beta-strand/beta-turn motif determine the difference between the P and P2 rust resistance specificities in flax. Plant Cell. 2001, 13 (1): 163-178.PubMedPubMed Central
- Pontier D, Godiard L, Marco Y, Roby D: hsr203J, a tobacco gene whose activation is rapid, highly localized and specific for incompatible plant/pathogen interactions. Plant J. 1994, 5 (4): 507-521. 10.1046/j.1365-313X.1994.5040507.x.PubMedView Article
- Clough SJ, Fengler KA, Yu IC, Lippok B, Smith RK, Bent AF: The Arabidopsis dnd1 "defense, no death" gene encodes a mutated cyclic nucleotide-gated ion channel. Proc Natl Acad Sci USA. 2000, 97 (16): 9323-9328. 10.1073/pnas.150005697.PubMedPubMed CentralView Article
- Green DR, Reed JC: Mitochondria and apoptosis. Science. 1998, 281 (5381): 1309-1312.PubMedView Article
- Logan DC: The dynamic plant chondriome. Seminars in Cell & Developmental Biology. 2010, 21 (6): 550-557. 10.1016/j.semcdb.2009.12.010.View Article
- Martinou J-C, Desagher S, Antonsson B: Cytochrome c release from mitochondria: all or nothing. Nat Cell Biol. 2000, 2 (3): E41-E43. 10.1038/35004069.PubMedView Article
- Sanchez P, de Torres Zabala M, Grant M: AtBI-1, a plant homologue of Bax inhibitor-1, suppresses Bax-induced cell death in yeast and is rapidly upregulated during wounding and pathogen challenge. Plant J. 2000, 21 (4): 393-399. 10.1046/j.1365-313x.2000.00690.x.PubMedView Article
- Hoeberichts FA, Woltering EJ: Multiple mediators of plant programmed cell death: Interplay of conserved cell death mechanisms and plant-specific regulators. BioEssays. 2003, 25 (1): 47-57. 10.1002/bies.10175.PubMedView Article
- Li J, Brader G, Palva ET: Kunitz trypsin inhibitor: an antagonist of cell death triggered by phytopathogens and fumonisin b1 in Arabidopsis. Mol Plant. 2008, 1 (3): 482-495. 10.1093/mp/ssn013.PubMedView Article
- Xu R, Song F, Zheng Z: OsBISAMT1, a gene encoding S-adenosyl-L-methionine: salicylic acid carboxyl methyltransferase, is differentially expressed in rice defense responses. Mol Biol Rep. 2006, 33 (3): 223-231. 10.1007/s11033-005-4823-x.PubMedView Article
- Kumar D, Klessig DF: High-affinity salicylic acid-binding protein 2 is required for plant innate immunity and has salicylic acid-stimulated lipase activity. Proc Natl Acad Sci USA. 2003, 100 (26): 16101-16106. 10.1073/pnas.0307162100.PubMedPubMed CentralView Article
- Gibly A, Bonshtien A, Balaji V, Debbie P, Martin GB, Sessa G: Identification and expression profiling of tomato genes differentially regulated during a resistance response to Xanthomonas campestris pv. vesicatoria. Mol Plant Microbe Interact. 2004, 17 (11): 1212-1222. 10.1094/MPMI.2004.17.11.1212.PubMedView Article
- Lund ST, Stall RE, Klee HJ: Ethylene regulates the susceptible response to pathogen infection in tomato. Plant Cell. 1998, 10 (3): 371-382.PubMedPubMed CentralView Article
- Heller RA, Schena M, Chai A, Shalon D, Bedilion T, Gilmore J, Woolley DE, Davis RW: Discovery and analysis of inflammatory disease-related genes using cDNA microarrays. Proc Natl Acad Sci USA. 1997, 94 (6): 2150-2155. 10.1073/pnas.94.6.2150.PubMedPubMed CentralView Article
- Smyth GK, (ed): Limma: Linear models for microarray data. New York:Springer; 2005.
- Yang YH, Dudoit S, Luu P, Lin DM, Peng V, Ngai J, Speed TP: Normalization for cDNA microarray data: a robust composite method addressing single and multiple slide systematic variation. Nucleic Acids Res. 2002, 30 (4): e15-10.1093/nar/30.4.e15.PubMedPubMed CentralView Article
- Conesa A, Nueda MJ, Ferrer A, Talon M: maSigPro: a method to identify significantly differential expression profiles in time-course microarray experiments. Bioinformatics. 2006, 22 (9): 1096-1102. 10.1093/bioinformatics/btl056.PubMedView Article
- Nueda MJ, Sebastian P, Tarazona S, Garcia-Garcia F, Dopazo J, Ferrer A, Conesa A: Functional assessment of time course microarray data. BMC Bioinformatics. 2009, 10 (Suppl 6): S9-10.1186/1471-2105-10-S6-S9.PubMedPubMed CentralView Article
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.