- Research article
- Open Access
Dynamics in the resistant and susceptible peanut (Arachis hypogaea L.) root transcriptome on infection with the Ralstonia solanacearum
© Chen et al.; licensee BioMed Central Ltd. 2014
- Received: 1 May 2014
- Accepted: 5 November 2014
- Published: 7 December 2014
Bacterial wilt caused by Ralstonia solanacearum is a serious soil-borne disease of peanut (Arachis hypogaea L). The molecular basis of peanut response to R. solanacearum remains unknown. To understand the resistance mechanism behind peanut resistance to R. solanacearum, we used RNA-Seq to perform global transcriptome profiling on the roots of peanut resistant (R) and susceptible (S) genotypes under R. solanacearum infection.
A total of 4.95 x 108 raw sequence reads were generated and subsequently assembled into 271, 790 unigenes with an average length of 890 bp and a N50 of 1, 665 bp. 179, 641 unigenes could be annotated by public protein databases. The pairwise transcriptome comparsions of time course (6, 12, 24, 48 and 72 h post inoculation) were conducted 1) between inoculated and control samples of each genotype, 2) between inoculated samples of R and S genotypes. The linear dynamics of transcriptome profile was observed between adjacent samples for each genotype, two genotypes shared similar transcriptome pattern at early time points with most significant up regulation at 12 hour, and samples from R genotype at 24 h and S genotype at 48 h showed similar transcriptome pattern, significant differences of transcriptional profile were observed in pairwise comparisons between R and S genotypes. KEGG analysis showed that the primary metabolisms were inhibited in both genotypes and stronger inhibition in R genotype post inoculation. The defense related genes (R gene, LRR-RLK, cell wall genes, etc.) generally showed a genotype-specific down regulation and different expression between both genotypes.
This transcriptome profiling provided the largest data set that explores the dynamic in crosstalk between peanut and R. solanacearum. The results suggested that the down-regulation of primary metabolism is contributed to the resistance difference between R and S genotypes. The genotype-specific expression pattern of defense related DEGs also contributed to the resistance difference between R and S genotype. This study will strongly contribute to better understand the molecular interaction between plant and R. solanacearum.
- Arachis hypogaea L
- Ralstonia solanacearum
Peanut (Arachis hypogaea L.) is an important crop for oil and protein production in the tropical and subtropical regions of the world, Asia is the top region of peanut planting with the growing area of 11 million ha (http://faostat.fao.org/site/339/default.aspx). The bacterial wilt caused by Ralstonia solanacearum is a destructive soil borne peanut disease in Asia. It could cause peanut production reduction at least of 10% even a mass extinction. R. solanacearum has a wide host range expanding over more than 200 plant species . The process of R. solanacearum infecting plant had been well characterized in model crop, briefly, the R. solanacearum penetrated into cortical tissue of host roots, colonized and exploded in numbers, caused a sudden deadly wilt of plant [2–4]. The resistant breeding is the most ideal strategy for controlling bacterial wilt with great benefit of economy and environmental protection . However, even in the resistance varieties, R. solanacearum can multiply in a high level of number and caused the symptoms of stunted growth, weak wilting and reduced resistance to other pathogens, finally resulted in a potential crop failure . Up to now, the molecular basis of peanut resistance to R. solanacearum is poorly understood.
Understanding the complexity of disease resistance will contribute to the development of peanut resistance to bacterial wilt. In the past few decades, the molecular cross-talk between plants and pathogens had been characterized, the intrinsic mechanism of plant resistance to pathogens had been well documented [6–8]. During plant-pathogen interactions, plant evolved a two-tiered innate immunity system to defend against pathogens attack. The host cell surface localized pattern-recognition receptors (PRRs) recognizes pathogen-associated molecular patterns (PAMPs), then activated the PAMP-triggered immunity (PTI) followed by pathogen effector-triggered immunity (ETI). In PTI and ETI, a set of defense response on transcriptome level were activated, and resulted in the arrest of pathogen clone [9–12]. However, the molecular reaction between plant and R. solanacearum have received far less attentions.
Efforts had been made in discovering the molecular mechanisms underlying interactions between Arabidopsis–, tomato–, potato– R. solanacearum, several resistance related genes and enzymes had been well characterized [13–18]. The significant changes on level of transcriptome and proteome were also observed in interaction between plant and R. solanacearum [19–24]. Especially, the mechanism of silicon in priming tomato resistance to R. solanacearum has been systematically studied [18, 21, 24], roles of cell wall proteins in tomato defend against R. solanacearum were well discussed [25–30]. These have provided preliminary understanding of molecular mechanism of plant response to R. solanacearum. Up to now, the resistant mechanism of plant to R. solanacearum is obscure, the literature and molecular resources available for plant resistance to R. solanacearum remain to be enriched. About 68, 094 ESTs differentially expressed in plant after R. solanacearum challenge were identified in previous study (http://www.ncbi.nlm.nih.gov/nucest/?term=ralstonia). The present informations are too poor to clearly shed light on the mechanism behind plant resistance to R. solanacearum. The genetic control factors that determined the consequence between plant and R. solanacearum are not yet fully identified, and changes in the global transcriptome of plant resistance to R. solanacearum is yet to be explored. It is still a great challenge to isolate genes by map-based cloning for the huge genome size of peanut. The traditional methods of sequencing cDNA clones resulted in the loss of rare transcript with inefficient cost, low throughput and lack of quantitation of the transcripts. Recent years, RNA sequencing technologies were developed, it is a more comprehensive and efficient way to carry study of transcriptome level on detecting the expression pattern, explore new exons and novel genes [31, 32]. Especially, the application of this technology is not limited to the prior knowledge of genomic sequence, it had been successfully applied in peanut transcriptome study on development and response to stress [33–36].
In the present study, we invested globally and compared the transcriptome profile in the roots of peanut resistant (R) and susceptible (S) genotypes under R. solanacearum infection. The dynamic differences of transcriptome profiles in peanut roots under R. solanacearum infection were investigated. The specific transcripts related to peanut response to R. solanacearum were identified. The possible roles of differentially expressed unigenes (DEGs) were discussed and the resistant mechanism of peanut to R. solanacearum was also deduced. In addition, a better understanding of peanut resistance to R. solanacearum could be a reference for exploring the resistance to bacterial pathogen in other crop plants. This study also provided a significant transcriptome resource in systemic plant-R. solanacearum interactions.
Observation of the bacterial number in peanut roots post-inoculation
Transcriptome sequencing and de novoassembly
The sumuary of transcriptome
Number of reads
High quality length(bp)
Number of unigenes
To characterize the unigenes, they were searched against Nr protein database using Blastx. For all 271, 790 unigenes, 92, 147 (33.91%) unigenes had no match with Nr record, 179, 643 (66.09%) unigenes matched protein accessions in Nr database, only 100, 336 of 271, 790 unigenes (36.92%) were assigned valuable Nr annotations (Additional file 2: Table S1, Additional file 3: File S1). For sequence similarity analyses by Blast against Nt database, 109, 845 (40.42%) unigenes had matched significant similar sequences (Additional file 2: Table S1, Additional file 4: File S2). The 45.46% sequences have strong homology at the E-value of ≤1.0e-50 of top hit in the Nt database. (Additional file 4: File S2). To get a more overall view of present sequences, they were searched against special database of Arachis species, Glycine max and Arabidopsis with a cutoff of E-value of 1.0E-10. For the increasing of peanut transcripts from Sanger and RNA-sequencing, a present largest database for Arachis species database (peanut DB) were built (http://bioinfolab.muohio.edu/txid3818v1/) . Up to 109, 968 (47.03%) of 271, 790 unigenes matched 27, 694 contigs of Arachis species. 47, 201 unigenes matched to 21, 576 sequence from Glycine max and 107, 399 unigenes matched to 8, 299 genes from Arabidopsis, respectively (Additional file 5: File S3). It also observed that 1, 936 unigenes matched 1, 063 genes of R. solanacearum genome (Additional file 6: File S4). This observation proved a reliable representation of the comprehensive transcriptome in which even R. solanacearum transcripts were detected.
To further characterized sequence annotation, all unigenes were searched against COG database to analyze phylogenetically widespread domain families. From present unigenes set, 92, 326 of 271, 790 unigenes have significant homologies in COG database (Additional file 2: Table S1, Additional file 7: File S5). A total of 59, 472 unigenes were assigned to the 25 functions categories (Additional file 8: Figure S2). Among this, the top six categories included “Replication, recombination and repair” (13, 226, 14.36%), “translation, ribosomal structure and biogenesis” (13, 194, 14.29%), “Transcription” (12, 944, 14.01%), “posttranslational modification, protein turnover, chaperones” (10, 646, 11.53%) and “Signal transduction mechanisms” (10, 105, 11.96%). 1, 596 (1.72%) of unigenes were grouped into “defense mechanism”.
Following the Nr annotations, all unigenes were mapped into the records of GO database. 111, 663 of them retrieved 333, 350 annotations from 52 sub-categories of three GO categories. “Biological Process” took up the majority of GO annotations (93, 914, 28.17%), followed by “Molecular Function” (200, 315, 60.09%) and Cellular Component (39, 183, 11.75%) (Additional file 2: Table S1, Additional file 9: Figure S3, Additional file 10: File S6). Due to the limited descriptions, a small proportion of unigenes with GO distributions was determined.
Cluster analysis of transcriptome
Identification of DEGs
The Co-expression pattern of DEGs
To identify genes shared same expression pattern, the co-expression patterns of DEGs were further investigated in RD and SD data sets (Additional file 15: Table S4) and D data sets (Additional file 16: Table S5). It uncovered various and genotype-specific expression patterns of DEGs, and the down-regulation and less share of co-expressed DEGs were dominated in both genotypes. The most popular expression patterns shared by both genotypes were down-regulated meanwhile at 6, 12 and 24 h, there were 1, 979 and 1, 911 DEGs in R and S genotype respectively, 852 DEGs were shared by both genotypes.
Experimental verification of DEGs
To verify the reliability of sequencing-based approach in identifying R. solanacearum -responsive genes, we monitored the expression pattern of twelve candidate DEGs at the five time points post inoculation for both genotypes using qRT-PCR. These candidate DEGs included genes that were proved to be related to defend response in other plant species, such as R genes and genes involving in ethylene signal transduction pathways. Their expression showed basically linear correlation to the results of the RNA- sequencing (Additional file 17: Table S6, Additional file 18: Figure S4).
The annotation of DEGs
Similar results were observed in metabolism of energy and lipid, DEGs with down-regulation were more than those with up-regulation at 6, 12 and 24 h in S genotype and all time points in R genotype (Additional file 20: Figure S5 and Additional file 21: Figure S6). For D series data sets, the more DEGs with down-regulation involved these two pathways than those with up-regulation at 12, 24 and 48 h. The main related pathway was “Oxidative phosphorylation” and “Fatty acid metabolism”, the expressions of DEGs these two pathways were decreased hundreds of times than those in S genotype (Additional file 19: File S9).
Interestingly, it observed that the numbers of the up-regulated DEGs were close to or more than those of down-regulated DEGs at 12 and 48 h in S genotype. However, the DEGs with up regulation involved in enrichment analysis were significantly less than those with down regulation, the number declined disproportionately. For example of the DEGs involved in metabolism of carbohydrate, energy and lipid at 12 h, the percentages of up-regulation DEGs were almost half or less than those of down-regulated DEGs, they were 1.90% (193/10, 159) and 3.91% (388/9, 916) for carbohydrate metabolism, 1.75% (178/10, 159) and 3.0% (300/9, 916), 1.22% (124/10, 159) and 2.68% (266/9, 916). Further analysis showed that, compared with down-regulated DEGs, low percentage of KEGG annotations occured on up-regulated DEGs in R genotype at 12 h and in S genotype at 12 and 48 h.
The DEGs with always up-regulation in D data set
Zinc finger, BED-type predicted
D-binding pseudobarrel domain
Domain of unknown function
Domain of unknown function
Proteise inhibitor I25,
Myb transcription factor
The DEGs with always down-regulation in D data set
Ubiquitin thiolesterase 1-like
FAR1 D binding domain
Sugar nucleotide epimerase
Domain of unknown function
Disease resistance protein
AAA + ATPase
Myb transcription factor
Ribosomal protein L7
11-S seed storage protein
The expression patterns of defense-related genes
The expression patterns of genes taking key roles in plant resistance to pathogens were investigated, including those of the NBS-LRR (R) genes, the genes encoding leucine-rich repeat receptor-like protein (LRR-RLK), MAP kinase and WRKY factors, the DEGs related to ADP ribosylation, the genes involved in synthesis of phytoalexins including lectin, terpenoid and chalcone, and the genes involved in ethylene, salicylic acid and cytokinkins signaling pathways, and the factors involved in defense pathways such as reactive oxygen species (ROS), cell wall proteins and pathogenesis-related proteins (PR proteins).
Total 168 NBS-LRR type R genes were identified (Additional file 22: File S10). They were mainly down-regulated continuously at 12, 24 and 48 in R genotype and intermittently down-regulated at 12 and 48 h in S genotype (Figure 8C, D), and showed slightly up- or down-regulation between R and S genotype at single time points. Compared to S genotype, only one (Ahy163031) was up-regulated at five time points in R genotype, the DEGs encoding Glycine max R gene RPM1 (Ahy090925) and one R gene (Ahy077818) were significantly up-regulated (above 8 fold change) at 12 h in R genotype. Interestingly, total 35 DEGs encoding Medicago truncatula disease resistance protein RPS1, twenty-eight of them were shared by both genotypes.
Total 30 DEGs encoding WRKY factor (IPR003657) were identified (Additional file 22: File S10). Two DEGs, Ahy114509 and Ahy109642 were down-regulated at all five time points in R genotype and without expression difference in S genotype, 11 DEGs were meanwhile up-regulated at 12 h in both genotypes including DEGs encoding WRKY3, WRKY4, WRKY33, WRKY40, WRKY42, WRKY65 and WRKY70 of Glycine max. Four DEGs were significantly down-regulated, WRKY50 (Ahy108250) at 12 h and WRKY20 (Ahy066190) at 48 h in R genotype, WRKY17 (Ahy012174) at 12 h and WRKY48 (Ahy099173) at 48 h in S genotype. Compared with the expression pattern of R genes, the expression of WRKY showed an opposite pattern with up-regulation at 12 h when the R genes were down-regulated (Figure 8E,F).
Total 22 DEGs encoding MAP kinase (IPR003527) were identified, most of them were down-regulated at 12 and 24 h in R genotype and 6, 12, 24 and 48 h in S genotype (Additional file 22: File S10). DEGs encoding MEK2 were down-regulated in both genotype, DEGs encoding MEK4 was up-regulated at 6 h in R genotype. Similar expression pattern were observed in calcium-dependent protein kinases (IPR020636), they were mainly down-regulated in both genotypes (Additional file 22: File S10).
Five DEGs encoding ADP-ribosylation factor (IPR005502) were identified with down regulation in R or S genotype (Additional file 22: File S10). And one (Ahy245540) was up-regulated post inoculation and with no detection in control samples in both genotypes. Compared with S genotype, the expression of Ahy245540 was down-regulated at 48 h and then significantly up-regulated at 72 h in R genotype.
For the DEGs involving in the biosynthesis of phytoalexins, most of them were down-regulated at both genotypes with different expression patterns at later stages. For DEGs encoding domain of lectin (235 DEGs), compared with control, the down-regulation is mainstream in both genotypes (Additional file 23: File S11). Interestingly, all eleven DEGs with ricin B lectin domain were down-regulated in R genotype meanwhile five of them were significantly up-regulated in S genotype.
Twenty-one DEGs encoding chalcone synthase were significant down-regulated at 6 h in R genotype and at 12 h in S genotype (Additional file 23: File S11). For DEGs encoding terpenoid synthase, significant down-regulated occured 12 h in both genotypes (Additional file 23: File S11).
In ROS, the DEGs encoding peroxidase (POD) were down regulated with distinct expression pattern in both genotypes and showed significant up-regulation in S genotype when compared to R genotype (Additional file 23: File S11). The expression pattern of cell wall protein also showed differential regulation (Additional file 22: File S11). DEGs encoding callose synthase (CalS) were down regulated in both genotypes and with higher expression in S genotype when compared to R genotype. On the contrary, DEGs encoding pectinesterase (PE) and Polygalacturonase (PG) were down-regulated in S genotype, one PG (Ahy107854) were significant down-regulated at 48 h in S genotype when compared to R genotype, both are key enzymes in cell wall degradation (Additional file 23: File S11). Interestingly, in the COG annotation, the DEGs involved in cell wall biosynthesis generally showed significant up regulation in S than in R genotype (Additional file 24: Figure S7). DEGs encoding PR protein showed different regulation. Compared to R genotype, DEGs encoding PR-2 (β-1,3-glucanase, Ahy269290) were no detected in all samples (except SC1) in R genotypes with a significant up-regulation at 72 h in S genotype, PR-1 (Ahy193681) and PR-5 (Thaumatin-like protein, Ahy119583) were significantly down-regulated at 48 and 72 h respectively in S genotype, one PR-10 (Ahy120907) was up-regulated at 6 and 48 h respectively in S genotype, one PR-10 (Ahy249270) was significantly down regulated in both genotypes with no difference between them.
Expression patterns of the genes in plant hormone signaling pathway
Total 67 ethylene responsive factors (IPR001471) were identified and showed distinct expression pattern (Additional file 25: File S12). DEGs encoding ERF4 and ERF6 were up-regulated in both genotypes. Three DEGs encoding salicylic acid carboxyl methyltransferase (SAMT), all were down-regulated at R and two were down-regulated in S genotype (Additional file 25: File S12).
The distinct expression patterns of DEGs involving in plant hormone occured on DEGs involved in biosynthesis of zeatin. Compared with control, most DEGs encoding zeatin showed up-regulation but no expression changes at 6 h in both genotypes (Additional file 25: File S12). There were 10 DEGs involving in zeatin biosynthesis in both genotypes and another four were detected only in S genotype, and only one was down-regulated in both genotypes and two were down-regulated in S genotype.
Assembly and annotation of Transcriptome
RNA-seq technologies are powerful for de novo transcripts assembly, identification of novel genes and analysis of transcriptional profiles [38, 39]. In present study, RNA-seq enabled the generation of huge amounts of sequence reads, for characterization of the transcriptional changes of peanut roots under R. solanacearum challenges. More than 90 Gbp were generated. The parameters of present transcriptome were far more than those of previous reports in the number of total bp and reads, value of N50, average length of unigenes and coverage of genome sequence. The characters of other transcriptome researches were shown here: covering 211 Mbp of both genomes and 21, 714 contigs with average length of 494 bp for Arachis wild species , 27.2 million paired-end reads and 59, 077 unigenes with average length of 619 bp and N50 value of 823 , 114.87 Mbp of genome sequence, 74, 974 contigs with average length of 859 bp and N50 value of 974 , and the coverage of genome sequence were also more than the data from published database which represented 34.41 Mb (http://www.ncbi.nlm.nih.gov/UniGene/UGOrg.cgi?TAXID=3818)) and 33.97 Mb (http://bioinfolab.muohio.edu/txid3818v1/) genome sequences, respectively. It is well known that due to the lack of available annotations for the limited basic research in peanut, the huge number of unigenes (more than 33.91% of total unigenes) had no any matches in Nr database, it added the difficulty of evaluation and understanding of peanut transcriptome.
Similar transcriptome profiles between R and S genotypes
We performed a time course transcriptome analysis to investigate the roots global response between peanut R and S genotypes during R. solanacearum colonization. The previous reported that the plant resistance to R. solanacearum were mainly focused on stem (18,24,29) as well as our studies (5), however, the inhibition of root development was also observed in apparently healthy peanut (5), our primary study had revealed a significant difference on root bacterial population at 96 h post inoculation between R and S genotype (Figure 1). So study on the roots response to R. solanacearum could contribute to comprehensively understand the mechanism of plant resistance to R. solanacearum. The cluster analyses suggested a high synchronization in root development of both genotypes (RC1 Vs SC1, RC2 Vs SC2) and early response to R. solanacearum infection (R6 Vs S6, R12 Vs S12). Unexpectedly, the cluster analysis failed to build similarity between the transcriptome profiles of RC1 Vs RC2 and SC1 Vs SC2, showed a dramatic transcriptome change even in the developmental stage of root during only in 72 hours. Interestingly, the similarity patterns of transcriptomes were shared by R24 and S48, which indicated an earlier and more rapid molecular response to R. solanacearum infection in R genotype. Although the great difference of an order of magnitude of the bacterial population was observed at 96 hour post inoculation, however, the different transcriptomes were observed as early as 24 h post inoculation, it could be deduced that the resistance to R. solanacearum was determined earlier in R genotype before the multiplication of bacteria was inhibited in roots.
The primary metabolisms were inhibited after inoculation in both genotypes
For RD, SD and D series datasets, the enrichment of DEGs acted in primary metabolism was observed by KEGG analyses. The collective behavior of DEGs indicated that the primary metabolisms were inhibited at early stage in S and whole stage in R genotype, and more evident inhibition of primary metabolism occured on R genotype. Up to now, the roles of primary metabolism in plant defense are not clear and less characterized [40, 41]. Although a number of studies showed that the genes acted in primary metabolism were up-regulated and identified as positive roles in plant resistance response . Interestingly, in the cases of wheat response to powdery mildew and Arabidopsis thaliana response to Pseudomonas syringae pv. Tomato, the genes acting in primary metabolism (glycolysis, the Krebs cycle and the pentose phosphate pathway) were down-regulated, the switchs from primary metabolism to defense metabolism were observed [42, 43]. The change of primary metabolism in host plant caused by pathogen had been documented [40, 44, 45]. It is a consensus that the host plant defense response to pathogen and inhibition growth of pathogen is cost-intensive, the nutrition and energy were transported from primary metabolism to the defense interaction . However, it also showed a very complicated regulation of genes involving in the primary metabolism for their up- or down-regulation even in the same pathways [44, 45]. Interestingly, in present study, the expression pattern of DEGs encoding sugar transporter (IPR005828 ) were opposite to those were previous reported. The sugar transporter were induced to supply nutrient for the pathogen colonization in the cases of Arabidopsis to bacterium Pseudomonas syringae pv. tomato DC3000 and powdery mildew fungus G. cichoracearum, rice to Xanthomonas oryzae pv. Oryzae [46, 47]. It seemed that the peanut reduced the nutrition supply to R. solanacearum by down-regulation of sugar transport. In the case of peanut bacterial wilt, even in the incompatible interactions between peanut R genotype and R. solanacearum, the concentration of R. solanacearum is up to the 108cfu/ml even in the R genotype during a long period post inoculation. It could be assumed that the more nutrition and energy was switched to the defense pathways by more intensive and longer-lasting down-regulation of primary metabolism in R genotype, and the down-regulation of primary metabolism is contributed to the resistance difference between R and S genotypes.
The PTI and ETI were partly suppressed in both genotypes
The molecular event occured in ETI and PTI had been documented and the key components of PTI and ETI were characterized well [48–50]. The sequential molecular events caused by pathogens were perception of PAMP by PRR, and then activation of the MAPK cascade followed the activation of WRKY-type transcription and R gene. And the leucine-rich repeat receptor kinases (LRR-RK) were identified as one of PRR. In Arabidopsis, the LRR-RLK, named FLS2, can bind bacterial flagellin peptides and contribute to resistance . FLS2 is a component of the preexisting recognition system, its expression was not affected by flagellin. However, in present study, the expressions of DEGs encoding LRR-RLK including FLS2-like kinase were mainly down-regulated, and one LRR-RLK (Ahy086504) was not detected in all samples from S genotype, the role of Ahy086504 deserves to be further studied. After perception of PAMP, the MAPK cascade was activated and MKK4 acted positively and MKK2 acted negatively in PTI in Arabidopsis [52, 53]. The expression pattern indicated that a similar regulation of MAPK component occured on the interaction of peanut to R. solanacearum like that in Arabidopsis to bacterial pathogen. The MAPK cascade then activated WRKY-type transcription factors which act in a complex defense response network as both positive and negative regulators in plant immunity . The WRKY27 took negative role in Arabidopsis response to R. solanacearum, on the contrary, the CaWRKY27 promoted tobacco resistance to R. solanacearum [14, 55]. The CaWRKY40 took a positive role in pepper resistance to R. solanacearum . The WRKY70 were positively regulated in Brachypodium distachyon response to pathogens . In this study, the identified DEGS encoding WRKY showed distinct expression pattern in both genotypes, WRKY33, -40, -70 were up-regulated and WRKY17, -20, -48 and -50 were down-regulated, it could be deduced that the WRKYs acted positively and negatively in peanut resistance to R. solanacearum, the exact roles of WRKY still need further study. Based on the expression pattern analysis of components in PTI, with consideration of the increased bacterial populations in both genotypes (Figure 1), it seemed that the components of PTI were activated, however, the R. solanacearum partly inhibited PTI in both genotypes.
In the followed ETI, the two types of NBS-LRR genes, CC-NBS-LRR and TIR-NBS-LRR were identified and with general down-regulation post-inoculation, and three R genes showed higher expression in R than in S genotype, it argued that the NBS-LRR type of R genes were responsible for peanut resistance to R. solanacearum. The R gene was the central component in ETI by detecting and binding to bacterial effectors and triggered the subsequent defense response. More interestingly, one CC-NBS-LRR type of R gene, the RPS1 of Medicago truncatula was represented by 35 unigenes with down-regulation in both genotypes. The soybean RPS1 (CC-NBS-LRR) was identified as the resistance gene to Phytophthora sojae which cause root and stem disease in soybean . It raised the suspicion that the RPS1 mediated partially the resistance of peanut to R. solanacearum. The DEGs encoding Glycine max R gene RPm1 (Ahy090925) were also identified with down-regulation in both genotypes and higher regulation in R genotype. The Arabidopsis RPm1 is a resistance gene responding to Pseudomonas syringae pv. tomato DC3000, it associates physically with bacterial FLS2 which destabilizes host ADP ribosylation factor AtMIN7 which is required for both PTI and ETI in Arabidopsis . The expression pattern of Ahy245540 encoding ADP ribosylation factor indicated that the degradation of ADP-ribosylation was more effectively blocked in R genotype, Ahy090925 and Ahy245540 maybe act like RPM1 and AtMIN7 respectively in Arabidopsis response to Pseudomonas syringae pv. tomato DC3000, thus the PTI in R genotype was rescued by ETI.
In PTI and ETI, the salicylic acid (SA) and ET signaling pathway involves in the activation of defense-related genes that contribute to resistance to pathogen [50, 60]. ERFs involved in SA and ET signaling pathways and take positive or negative roles in plant immunity response [60, 61]. In present study, many types of ERF were identified, ERF1, 4 and 6 were up-regulated and ERF3 was down-regulated in both genotypes, this reproduced the results of our previous report, and confirmed the biosynthesis of ET may contribute to peanut resistance to R. solanacearum . The plant cytokinins involved in the resistance to pathogen and promoted the accumulation of salicylic acid in plant defense response . In Arabidopsis, zeatin involved in the defense response of R gene RPM1 and caused the accumulation of salicylic acid [62, 63]. Recent research showed that SA positively regulated Arabidopsis reistance to R. solanacearum . Given the up-regulation of DEGs involved in zeatin biosynthesis and down-regulation of SAMT, which balances the level of salicylic acid by methylating SA , it appeared the promotion of zeatin lead to the accumulation of salicylic acid, both cytokinins and SA synergistically acted in peanut resistance to R. solanacearum in both genotypes, but did not contribute to the resistance difference between R and S genotype.
As the result of R gene-mediated resistance, the ROS and cell wall proteins as well as PR proteins contributed to a final immune response to pathogen [50, 60]. In present study, POD were up-regulated in S than in R genotype, which involved in production of ROS and reinforcement of the cell wall, increased tomato resistance to R. solanacearum [18, 24, 29]. Interestingly, compared to R genotype, DEGs involved in cell wall biosynthesis, such as CalS, POD and DEGs with COG’s annotation of cell wall biosynthesis were generally up-regulated in S genotype, meanwhile those involved in cell wall degradation (PE and PG) were down-regulated in S genotype. It could be understand that the DEGs involving in cell walls were inhibited in both genotypes under a high inoculation pressures, however, a more intense defense response in cell wall biosynthesis occurred in S genotype. The induction of PR protein indicates activation of systematic acquired resistance to pathogen, PR protein including PR1, PR2, PR5 and PR10 were identified as positive roles in tomato resistance to R. solanacearum (21,24,29). In present study, PR-2 were detected only in inoculated samples of S genotype and RC1, PR-5 were down regulated in S genotype meanwhile one PR-10 showed up regulation in S genotype, their expression were different from those in tomato, indicated a potential different role of them in peanut resistance to R. solanacearum.
In present study, we have investigated the transcriptome in roots of peanut R and S genotypes response to R. solanacearum infection at five time points. The dynamics profiles of transcriptome during the infection were roughly characterized. Both genotypes showed a general similarity on transcriptome level at early time points, meanwhile resistance-related DEGs showed distinct genotype-specific expression pattern in response to R. solanacearum challenge. The enrichment analyses for the DEGs showed genotype-specific patterns of transcriptome remodelling under R. solanacearum challenge. The present work provided the largest unigene dataset for cultivated peanut response to R. solanacearum infection, greatly expanded the vision of the genetic basis and improved our understanding on the molecular mechanisms underlying peanut resistance to R. solanacearum.
Plant inoculation and tissue harvest
Two peanut genotypes, R genotype J04 and S genotype J62, were provided by Oil Crop Research Institute, Chinese Academy of Agricultural Sciences. Two genotypes greatly differ in resistance to Ralstonia solanacearum and are close in genetic background . The seeds were sterilized and germinated on wet double filters paper in a culture dish and transplanted to a porcelain pot full of water when the roots were about 2 cm long, under 16 h photoperiod at 26°C with the relative humidity of 65-80%. The virulent R. solanacearum strain (Race 1, biovar 3) was isolated from soil of disease nurseries of Hongan County. Bacterial cells were cultured in CPG medium (pH 7.0) containing 0.1% casamino acids, 1% peptone and 0.5% glucose, at 200 rpm and 28°C in a shaker for 48 hours.
When seedlings were at 3-leaf stage, the root tips were cut to produce wound and then cultivated in a cultivation cabinet. After 24 hours when the root tips were cut, the roots were inoculated with R. solanacearum strain race 1 suspension (107cfu/ml) in a culture dishes for 30 min, meanwhile the mock inoculation seedlings were used as control. The roots of seven individual seedlings were sampled at 0, 6, 12, 24, 48 and 72 hours post inoculation. The samples were frozen in liquid nitrogen immediately and stored at -80°C. About 15 individuals of each line were kept for more than 10 days to confirm successful inoculation (Additional file 26: Figure S8).
Detection of population dynamics of R. solanacearumin roots
The roots of infected samples at each time point were vigorously grinded in a mortar for 1 g of roots per 1 ml distilled water, and the serial dilutions of extracted suspensions were made with dilution multiple from 103 to 108, and 0.1 ml aliquots were spread on the surface of a CPG medium then cultured in a CPG medium. The medium is made of 10 g bacto peptone (Difco), 5 g Glucose, 1 g casamino acid (Difco), 15 g bacto agar (Difco), and 1 L distilled water. This medium was supplemented with 1% polymyxin B sulphate (Sigma), 1% crystal violet, 1% bacitracin (Sigma), and 1% cycloheximide (Sigma). After incubating plates at 28°C for 48 hours, colonies of R. solanacearum were counted and cfu were calculated per gram of roots. Three replicates were prepared for each sample.
Construction of the transcriptome
For samples of two genotypes, RNAs were isolated from pooled roots using Trizol reagent (Invitrogen, US), then analysed with Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA). The RNA samples were numbered with the sampling time point of 6, 12, 24, 48 and 72 h. For inoculated samples, they were denoted as S6, S12, S24, S48 and S72 for S genotype; and R6, R12, R24, R48 and R72 for R genotype. For control samples, the equal amount of RNA from time points 6, 12 and 24 h were pooled and denoted as RC1 and SC1 for R and S genotype, the pools of RNAs from 48 and 72 h were denoted as RC2 and SC2 for R and for S genotype. Rna-sequencing was carried by Beijing Genomics Institute for Illumina sequencing on a HiSeq2000 system. The raw sequence from fourteen libraries were assembled into comprehensive unigenes using Trinity and TGICL [27, 29].
Transcriptome functional annotation
The assignment of sequence orientations and functional annotations of the all-unigenes were determined by BLASTx against the followed database, the NCBI nonredundant (Nr) protein database, the Swiss-Prot protein database, the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway database and the Cluster of Orthologous Groups (COG) database with an E-value cut-off of 1.0e–5. The all-unigenes were assigned GO annotations using Blast2GO (http://www.blast2go.rog/). In addition, unigenes were aligned with the NCBI nucleotide (Nt) databases using BLASTn with an E value of 1.0e–5.
First, All reads of each libraries were respectively mapped onto the unigenes using the default parameters by SOAP . Second, the uniquely mapped reads were extracted for quantifing the abundance. Third, Unigene expression was normalized using the value of RPKM (reads per kilobase per million reads). Multiple pairwise comparisons were carried between the data sets of different samples (Additional file 12: Table S2). The absolute value of log2(Ratio) ≥ 1 (under the criterion of P ≤ 0.01 and FDR ≤0.001) were used as threshold to assess the significance of gene expression difference.
The reverse transcripts were performed using an Invitrogen SuperScript Reagent Kit. The primer designed using the Oligo6 software. For RT-PCR, the SYBR® Premix Ex Taq™ (TAKARA) was used on a Bio-Rad CFX96 real-time PCR detection system (Bio-Rad, Hercules, CA). Gene expression was analyzed for samples at 6, 12, 24, 48 and 72 hours post-inoculation of R and S genotypes. All reactions for each gene were performed in triplicate. The relative expression level of each gene among samples was calculated using the 2-△△Ct method with normalization to the internal reference actin gene. The parameters of thermal cycle were 95°C for 30 s, followed by 40 cycles of 95°C for 10 s, 50–56°C for 25 s at a volume of 20 μl.
This research was financed by funds from the National Key Basic Research Special Foundation of China “973” Project (2011CB109304), the National Natural Science Foundation of China (31000724), the National Natural Science Foundation of China (31271764), the Plant Germplasm Resources Sharing Platform [NB2011-2130135-29(2)] and Agricultural Well Breeds Project of Shandong Province.
- Hayward AC: Diagnosis, Distribution and Status of Groundnut Bacterial Wilt. Proceedings of an ACIAR/ICRISAT Collaborative Research planning meeting: 18–19 March 1990; Genting Highlands, Malaysia. Edited by: Milddleton and Hayward, ACIAR Proc No 31, 58. 1990Google Scholar
- Van Overbeek LS, Cassidy M, Kozdroj J, Trevors JT, Van Elsas JD: A polyphasic approach for studying the interaction between Ralstonia solanacearum and potential control agents in the tomato phytosphere. J Microbiol Meth. 2002, 48: 69-86. 10.1016/S0167-7012(01)00347-5.View ArticleGoogle Scholar
- Götz M, Gomes NCM, Dratwinski A, Costa R, Berg G, Peixoto R, Mendonça-Hagler L, Smalla K: Survival of gfp-tagged antagonistic bacteria in the rhizosphere of tomato plants and their effects on the indigenous bacterial community. Fems Microbiol Ecol. 2006, 56: 207-218. 10.1111/j.1574-6941.2006.00093.x.PubMedView ArticleGoogle Scholar
- Fujie M, Takamoto H, Kawasaki T, Fujiwara A, Yamada T: Monitoring growth and movement of Ralstonia solanacearum cells harboring plasmid pRSS12 derived from bacteriophage ϕRSS1. J Biosc Bioeng. 2010, 109: 153-158. 10.1016/j.jbiosc.2009.07.012.View ArticleGoogle Scholar
- Liao B, Shan Z, Lei Y, Tan Y, Li D, Duan N: Reaction to latent infection by Ralstonia solanacearum in groundnut. Chin J Oil Crop Sci. 1998, 20: 61-65.Google Scholar
- Dangl JL, Jones JDG: Plant pathogens and integrated defence responses to infection. Nature. 2001, 411: 826-833. 10.1038/35081161.PubMedView ArticleGoogle Scholar
- McDowell JM, Simon SA: Molecular diversity at the plant–pathogen interface. Dev Comp Immunol. 2008, 32: 736-744. 10.1016/j.dci.2007.11.005.PubMedView ArticleGoogle Scholar
- Zhang J, Zhou J: Plant immunity triggered by microbial molecular signatures. Mol Plant. 2010, 3: 783-793. 10.1093/mp/ssq035.PubMedView ArticleGoogle Scholar
- Dong X: SA, JA, ethylene, and disease resistance in plants. Curr Opin Plant Biol. 1998, 1: 316-323. 10.1016/1369-5266(88)80053-0.PubMedView ArticleGoogle Scholar
- Wise RP, Moscou MJ, Bogdanove AJ, Whitham SA: Transcript profiling in host–pathogen interactions. Annu Rev Phytopathol. 2007, 45: 329-369. 10.1146/annurev.phyto.45.011107.143944.PubMedView ArticleGoogle Scholar
- Bonas U, Lahaye T: Plant disease resistance triggered by pathogen-derived molecules: refined models of specific recognition. Curr Opin Microbiol. 2002, 5: 44-50. 10.1016/S1369-5274(02)00284-9.PubMedView ArticleGoogle Scholar
- Göhre V, Robatzek S: Breaking the barriers: microbial effector molecules subvert plant immunity. Annu Rev Phytopathol. 2008, 46: 189-215. 10.1146/annurev.phyto.46.120407.110050.PubMedView ArticleGoogle Scholar
- Deslandes L, Olivier J, Theulières F, Hirsch J, Feng D, Bittner-Eddy P, Beynon J, Marco Y: Resistance to Ralstonia solanacearum in Arabidopsis thaliana is confered by the recessive RRS1-R gene, a member of a novel family of resistance genes. Proc Nalt Acad Sci. 2002, 99: 2404-2409. 10.1073/pnas.032485099.View ArticleGoogle Scholar
- Mukhtar MS, Deslandes L, Auriac MC, Marco Y, Somssich IE: The Arabidopsis transcription factor WRKY27 influences wilt disease symptom development caused by Ralstonia solanacearu. Plant J. 2008, 56: 935-947. 10.1111/j.1365-313X.2008.03651.x.PubMedView ArticleGoogle Scholar
- Gao G, Jin L, Xie K, Qu D: The potato StLTPa7 gene displays a complex Ca2+-associated patternof expression during the early stage of potato–Ralstonia solanacearum interaction. Mol Plant Pathol. 2009, 10: 15-27. 10.1111/j.1364-3703.2008.00508.x.PubMedView ArticleGoogle Scholar
- Zhang H, Zhang D, Chen J, Yang Y, Huang Z, Huang D, Wang X, Huang R: Tomato stress-responsive factor TSRF1 interacts with ethylene responsive element GCC box and regulates pathogen resistance to Ralstonia solanacearum. Plant Mol Biol. 2004, 55: 825-834. 10.1007/s11103-005-2140-3.PubMedView ArticleGoogle Scholar
- Gao G, Jin L, Xie K, Qu D: Differential space-time expression of StLTPb1 gene between resistant and susceptible potato genotypes in response to Ralstonia solanacearum. Agr Sci China. 2008, 7: 713-725. 10.1016/S1671-2927(08)60106-6.View ArticleGoogle Scholar
- Kurabachew H, Wydra K: Induction of systemic resistance and defense-related enzymes after elicitation of resistance by rhizobacteria and silicon application against Ralstonia solanacearum in tomato (Solanum lycopersicum). Crop Prot. 2014, 57: 1-7.View ArticleGoogle Scholar
- Hu J, Barlet X, Deslandes L, Hirsch J, Feng D, Somssich I, Marco Y: Transcriptional responses of Arabidopsis thaliana during wilt disease caused by the soil-borne phytopathogenic bacterium, Ralstonia solanacearum. PLoS One. 2008, 3: e2589-10.1371/journal.pone.0002589.PubMed CentralPubMedView ArticleGoogle Scholar
- Li G, Jin L, Wang X, Xie K, Yang Y, Van Der Vossen EAG, Huanga S, Qu D: Gene transcription analysis during interaction between potato and Ralstonia solanacearum. Russ J Plant Physl. 2010, 57: 685-695. 10.1134/S1021443710050122.View ArticleGoogle Scholar
- Ghareeb H, Bozsó Z, Ott PG, Repenning C, Stahl F, Wydra K: Transcriptome of silicon-induced resistance against Ralstonia solanacearum in the silicon non accumulator tomato implicates priming effect. Physiol Mol Plant Pathol. 2011, 75: 83-89. 10.1016/j.pmpp.2010.11.004.View ArticleGoogle Scholar
- Huang J, Yan L, Lei Y, Jiang H, Ren X, Liao B: Expressed sequence tags in cultivated peanut (Arachis hypogaea): discovery of genes in seed development and response to Ralstonia solanacearum challenge. J Plant Res. 2012, 125: 755-769. 10.1007/s10265-012-0491-9.PubMedView ArticleGoogle Scholar
- Kurabachew H, Stahl F, Wydra K: Global gene expression of rhizobacteria-silicon mediated induced systemic resistance in tomato (Solanum lycopersicum) against Ralstonia solanacearum. Physiol Mol Plant P. 2013, 84: 44-52.View ArticleGoogle Scholar
- Kiirika L, Stahl F, Wydra K: Phenotypic and molecular characterization of resistance induction by single and combined application of chitosan and silicon in tomato against Ralstonia solanacearum. Physiol Mol Plant P. 2013, 81: 1-12.View ArticleGoogle Scholar
- Wydra K, Beri H: Structural changes of homogalacturonan, rhamnogalacturonan I and arabinogalactan protein in xylem cell walls of tomato genotypes in reaction to Ralstonia solanacearum. Physiol Mol Plant P. 2006, 68: 41-50. 10.1016/j.pmpp.2006.06.001.View ArticleGoogle Scholar
- Diogo R, Wydra K: Silicon-induced basal resistance in tomato against Ralstonia solanacearum is related to modification of pectic cell wall polysaccharide structure. Physiol Mol Plant P. 2007, 70: 120-129. 10.1016/j.pmpp.2007.07.008.View ArticleGoogle Scholar
- Wydra K, Beri H: Immunohistochemical changes in methyl-ester distribution of homogalacturonan and side chain composition of rhamnogalacturonan I as possible components of basal resistance in tomato inoculated with Ralstonia solanacearum. Physiol Mol Plant P. 2007, 70: 13-24. 10.1016/j.pmpp.2007.05.006.View ArticleGoogle Scholar
- Dahal D, Heintz D, Van Dorsselaer A, Braun HP, Wydra K: Pathogenesis and stress related, as well as metabolic proteins are regulated in tomato stems infected with Ralstonia solanacearum. Plant Physiol Bioch. 2009, 47: 838-846. 10.1016/j.plaphy.2009.05.001.View ArticleGoogle Scholar
- Dahal D, Pich A, Braun HP, Wydra K: Analysis of cell wall proteins regulated in stems of susceptible and resistant tomato genotypes after inoculation with Ralstonia solanacearum: a proteomics approach. Plant Mol Biol. 2010, 73: 643-658. 10.1007/s11103-010-9646-z.PubMed CentralPubMedView ArticleGoogle Scholar
- Schacht T, Unger C, Pich A, Wydra K: Endo- and exopolygalacturonases of Ralstonia solanacearum are inhibited by polygalacturonase-inhibiting protein (PGIP) activity in tomato stem extracts. Plant Physiol Bioch. 2011, 49: 377-387. 10.1016/j.plaphy.2011.02.001.View ArticleGoogle Scholar
- Li C, Shao J, Wang Y, Li W, Guo D, Yan B, Xia Y, Peng M: Analysis of banana transcriptome and global gene expression profiles in banana roots in response to infection by race 1 and tropical race 4 of Fusarium oxysporum f. sp. cubense. BMC Genomics. 2013, 14: 851-10.1186/1471-2164-14-851.PubMed CentralPubMedView ArticleGoogle Scholar
- Varshney RK, Nayak SN, May GD, Jackson SA: Next generation sequencing technologies and their implications for crop genetics and breeding. Trends Biotechnol. 2009, 27: 522-530. 10.1016/j.tibtech.2009.05.006.PubMedView ArticleGoogle Scholar
- Pandey MK, Monyo E, Ozias-Akins P, Liang X, Guimarães P, Nigam SN, Upadhyaya HD, Janila P, Zhang X, Guo B: Advances in Arachis genomics for peanut improvement. Biotechnol Adv. 2012, 30: 639-651. 10.1016/j.biotechadv.2011.11.001.PubMedView ArticleGoogle Scholar
- Zhang J, Liang S, Duan J, Wang J, Chen S, Cheng Z, Zhang Q, Liang X, Li Y: De novo assembly and characterisation of the transcriptome during seed development, and generation of genic-SSR markers in Peanut (Arachis hypogaea L.). BMC Genomics. 2012, 13: 90-10.1186/1471-2164-13-90.PubMed CentralPubMedView ArticleGoogle Scholar
- Guimarães PM, Brasileiro ACM, Morgante CV, Martins AC, Pappas G, Silva OB, Togawa R, Leal-Bertioli SCM, Araujo ACG, Moretzsohn MC: Global transcriptome analysis of two wild relatives of peanut under drought and fungi infection. BMC Genomics. 2012, 13: 387-10.1186/1471-2164-13-387.PubMed CentralPubMedView ArticleGoogle Scholar
- Chen X, Zhu W, Azam S, Li H, Zhu F, Li H, Hong Y, Liu H, Zhang E, Wu H: Deep sequencing analysis of the transcriptomes of peanut aerial and subterranean young pods identifies candidate genes related to early embryo abortion. Plant Biotechnol J. 2013, 11: 115-127. 10.1111/pbi.12018.PubMedView ArticleGoogle Scholar
- Duan X, Schmidt E, Li P, Lenox D, Liu L, Shu C, Zhang J, Liang C: PeanutDB: an integrated bioinformatics web portal for Arachis hypogaea transcriptomics. BMC Plant Biol. 2012, 12: 94-10.1186/1471-2229-12-94.PubMed CentralPubMedView ArticleGoogle Scholar
- Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, Adiconis X, Fan L, Raychowdhury R, Zeng Q: Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011, 29: 644-652. 10.1038/nbt.1883.PubMed CentralPubMedView ArticleGoogle Scholar
- Zhao Q, Wang Y, Kong Y, Luo D, Xuan L, Hao P: Optimizing de novo transcriptome assembly from short-read RNA-Seq data: a comparative study. BMC Bioinformatics. 2011, 12 (Suppl 14): S2-10.1186/1471-2105-12-S14-S2.View ArticleGoogle Scholar
- Berger S, Sinha AK, Roitsch T: Plant physiology meets phytopathology: plant primary metabolism and plant-pathogen interactions. J Exp Bot. 2007, 58: 4019-4026. 10.1093/jxb/erm298.PubMedView ArticleGoogle Scholar
- Bolton MD, Thomma BPHJ: The complexity of nitrogen metabolism and nitrogen-regulated gene expression in plant pathogenic fungi. Physiol Mol Plant Pathol. 2008, 72: 104-110. 10.1016/j.pmpp.2008.07.001.View ArticleGoogle Scholar
- Scheideler M, Schlaich NL, Fellenberg K, Beissbarth T, Hauser NC, Vingron M, Slusarenko AJ, Hoheisel JD: Monitoring the switch from housekeeping to pathogen defense metabolism in Arabidopsis thaliana using cDNA arrays. J Biol Chem. 2002, 277: 10555-10561. 10.1074/jbc.M104863200.PubMedView ArticleGoogle Scholar
- Xin M, Wang X, Peng H, Yao Y, Xie C, Han Y, Ni Z, Sun Q: Transcriptome comparison of susceptible and resistant wheat in response to powdery mildew infection gluconeogenesis process. Genomics Proteomics Bioinformatics. 2012, 10: 94-106. 10.1016/j.gpb.2012.05.002.PubMedView ArticleGoogle Scholar
- Chou HM, Bundock N, Rolfe SA, Scholes JD: Infection of Arabidopsis thaliana leaves with Albugo candida causes a reprogramming of host metabolism. Mol Plant Pathol. 2000, 1: 99-113. 10.1046/j.1364-3703.2000.00013.x.PubMedView ArticleGoogle Scholar
- Charte J, Scho¨n H, Weis E: Photosynthesis and carbohydrate metabolism in tobacco leaves during an incompatible interaction with Phytophthora nicotianae. Plant Cell Environ. 2005, 28: 1421-1435. 10.1111/j.1365-3040.2005.01380.x.View ArticleGoogle Scholar
- Chen LQ, Hou BH, Lalonde S, Takanaga H, Hartung ML, Qu XQ, Guo WJ, Kim JG, Underwood W, Chaudhuri B, Chermak D, Antony G, White FF, Somerville SC, Mudgett MB, Frommer WB: Sugar transporters for intercellular exchange and nutrition of pathogens. Nature. 2010, 468: 527-532. 10.1038/nature09606.PubMed CentralPubMedView ArticleGoogle Scholar
- Chen LQ, Qu XQ, Hou BH, Sosso D, Osorio S, Fernie AR, Frommer WB: Sucrose efflux mediated by SWEET proteins as a key step for phloem transport. Science. 2012, 335: 207-211. 10.1126/science.1213351.PubMedView ArticleGoogle Scholar
- Jones JDG, Dangl JL: The plant immune system. Nature. 2006, 444: 323-329. 10.1038/nature05286.PubMedView ArticleGoogle Scholar
- Schwessinger B, Zipfel C: News from the frontline: recent insights into PAMP-triggered immunity in plants. Curr Opin Plant Biol. 2008, 11: 389-395. 10.1016/j.pbi.2008.06.001.PubMedView ArticleGoogle Scholar
- Zhang Y, Lubberstedt T, Xu M: The genetic and molecular basis of plant resistance to pathogens. J Genet Genomics. 2013, 40: 23-25. 10.1016/j.jgg.2012.11.003.PubMedView ArticleGoogle Scholar
- Gomez-Gomez L, Boller T: FLS2: an LRR receptor-like kinase involved in the perception of the bacterial elicitor flagellin in Arabidopsis. Mol Cell. 2000, 5: 1003-1011. 10.1016/S1097-2765(00)80265-8.PubMedView ArticleGoogle Scholar
- Asai T, Tena G, Plotnikova J, Willmann MR, Chiu WL, Gomez-Gomez L, Boller T, Ausubel FM, Sheen J: MAP kinase signalling cascade in Arabidopsis innate immunity. Nature. 2002, 415: 977-983. 10.1038/415977a.PubMedView ArticleGoogle Scholar
- Nicaise V, Roux M, Zipfel C: Recent advances in PAMP-triggered immunity against bacteria pattern recognition receptors watch over and raise the alarm. Plant Physiol. 2009, 150: 1638-1647. 10.1104/pp.109.139709.PubMed CentralPubMedView ArticleGoogle Scholar
- Pandey SP, Somssich IE: The role of WRKY transcription factors in plant immunity. Plant Physiol. 2009, 150: 1648-1655. 10.1104/pp.109.138990.PubMed CentralPubMedView ArticleGoogle Scholar
- Dang F, Wang Y, She J, Lei Y, Liu Z, Eulgem T, Lai Y, Lin J, Yu L: Overexpression of CaWRKY27, a subgroup IIe WRKY transcription factor of Capsicum annuum, positively regulates tobacco resistance to Ralstonia solanacearum infection. Physiol Plantarum. 2014, 150: 397-411. 10.1111/ppl.12093.View ArticleGoogle Scholar
- Dang F, Wang Y, Lu Y, Eulgem T, Lai Y, Liu Z, Wang X, Qiu A, Zhang T, Lin J: CaWRKY40, a WRKY protein of pepper, plays an important role in the regulation of tolerance to heat stress and resistance to Ralstonia solanacearum infection. Plant Cell Environ. 2013, 36: 757-774. 10.1111/pce.12011.PubMedView ArticleGoogle Scholar
- Wen F, Zhu H, Li P, Jiang M, Mao W, Ong C, Chu Z: Genome-wide evolutionary characterization and expression analyses of WRKY family genes in Brachypodium distachyon. DNA Res. in pressGoogle Scholar
- Bhattacharyya MK, Narayanan NN, Gao H, Santra DK, Salimath SS, Kasuga T, Liu Y, Espinosa B, Ellison L, Marek L: Identification of a large cluster of coiled coil-nucleotide binding site–leucine rich repeat-type genes from the Rps1 region containing Phytophthora resistance genes in soybean. Theor Appl Genet. 2005, 111: 75-86. 10.1007/s00122-005-1993-9.PubMedView ArticleGoogle Scholar
- Nomura K, Mecey C, Lee YN, Imboden LA, Chang J, He S: Effector-triggered immunity blocks pathogen degradation of an immunity-associated vesicle traffic regulator in Arabidopsis. Proc Natl Acad Sci. 2011, 108: 10774-10779. 10.1073/pnas.1103338108.PubMed CentralPubMedView ArticleGoogle Scholar
- Robert-Seilaniantz A, Grant M, Jones JD: Hormone crosstalk in plant disease and defense: more than just jasmonate-salicylate antagonism. Annu Rev Phytopathol. 2011, 49: 317-343. 10.1146/annurev-phyto-073009-114447.PubMedView ArticleGoogle Scholar
- Berrocal-Lobo M, Molina A, Solano R: Constitutive expression of Ethylene-Response-Factor1 in Arabidopsis confers resistance to several necrotrophic fungi. Plant J. 2002, 29: 23-32. 10.1046/j.1365-313x.2002.01191.x.PubMedView ArticleGoogle Scholar
- Choi J, Huh SU, Kojima M, Sakakibara H, Paek KH, Hwang I: The cytokinin-activated transcription factor ARR2: promotes plant immunity via TGA3/NPR1-dependent salicylic acid signaling in Arabidopsis. Dev Cell. 2010, 19: 284-295. 10.1016/j.devcel.2010.07.011.PubMedView ArticleGoogle Scholar
- Igari K, Endo S, Hibara K, Aida M, Sakakibara H, Kawasaki T, Tasaka M: Constitutive activation of a CC-NB-LRR protein alters morphogenesis through the cytokinin pathway in Arabidopsis. Plant J. 2008, 55: 14-27. 10.1111/j.1365-313X.2008.03466.x.PubMedView ArticleGoogle Scholar
- Denance N, Ranocha P, Oria N, Barlet X, Rivière MP, Yadeta KA, Hoffmann L, Perreau F, Clément G, Maia-Grondard A: Arabidopsis wat1 (walls are thin1)-mediated resistance to the bacterial vascular pathogen, Ralstonia solanacearum, is accompanied by cross-regulation of salicylic acid and tryptophan metabolism. Plant J. 2013, 73: 225-239. 10.1111/tpj.12027.PubMedView ArticleGoogle Scholar
- Ross JR, Nam KH, D'Auria JC, Pichersky E: S-adenosyl-L-methionine: salicylic acid carboxyl methyltransferase, an enzyme involved in floral scent production and plant defense, represents a new class of plant methyltransferases. Arch Biochem Biophys. 1999, 367: 9-16. 10.1006/abbi.1999.1255.PubMedView ArticleGoogle Scholar
- Jiang H, Chen B, Ren X, Liao B, Lei Y, Fu T, Ma C, Mace E, Crouch JH: Identification of SSR markers linked to bacterial wilt resistance of peanut with RILs. Chin J Oil Crop Sci. 2007, 29: 26-30.Google Scholar
- Li R, Li Y, Kristiansen K, Wang J: SOAP: short oligonucleotide alignment program. Bioinformatics. 2008, 24: 713-714. 10.1093/bioinformatics/btn025.PubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.