Transcriptomic dissection of the rice – Burkholderia glumae interaction
© Magbanua et al.; licensee BioMed Central Ltd. 2014
Received: 5 May 2014
Accepted: 26 August 2014
Published: 3 September 2014
Bacterial panicle blight caused by the bacterium Burkholderia glumae is an emerging disease of rice in the United States. Not much is known about this disease, the disease cycle or any source of disease resistance. To understand the interaction between rice and Burkholderia glumae, we used transcriptomics via next-generation sequencing (RNA-Seq) and bioinformatics to identify differentially expressed transcripts between resistant and susceptible interactions and formulate a model for rice resistance to the disease.
Using inoculated young seedlings as sample tissues, we identified unique transcripts involved with resistance to bacterial panicle blight, including a PIF-like ORF1 and verified differential expression of some selected genes using qRT-PCR. These transcripts, which include resistance genes of the NBS-LRR type, kinases, transcription factors, transporters and expressed proteins with functions that are not known, have not been reported in other pathosystems including rice blast or bacterial blight. Further, functional annotation analysis reveals enrichment of defense response and programmed cell death (biological processes); ATP and protein binding (molecular functions); and mitochondrion-related (cell component) transcripts in the resistant interaction.
Taken together, we formulated a model for rice resistance to bacterial panicle blight that involves an activation of previously unknown resistance genes and their activation partners upon challenge with B. glumae. Other interesting findings are that 1) though these resistance transcripts were up-regulated upon inoculation in the resistant interaction, some of them were already expressed in the water-inoculated control from the resistant genotype, but not in the water- and bacterium-inoculated samples from the susceptible genotype; 2) rice may have co-opted an ORF that was previously a part of a transposable element to aid in the resistance mechanism; and 3) resistance may have existed immediately prior to rice domestication.
KeywordsBacterial panicle blight Disease resistance Next-generation sequencing
Bacterial panicle blight (BPB) is an emerging disease on rice in the United States [1–4]. BPB was first reported in Japan as the cause of grain rotting and seedling blight on rice . The disease can also cause sheath rot and panicle blight . Diseased panicles are characterized by having florets with a darker base and a reddish-brown margin, and frequently upright due to poor filling. Epidemics of the disease occurred during 1995 and 2000 with yield losses in some fields estimated up to 40% [3, 6]; however, the disease can be worse in years with hot summers . With global temperature on the rise, it is expected that the disease will pose a threat to rice production worldwide . Currently, rice that is highly resistant to the disease is not available commercially and some rice varieties, such as Clearfield 161 showed promising, but not complete resistance to the disease . Major practices to control the disease are use of pathogen free seeds and application of antibiotics .
The bacterium Burkholderia glumae is the major causal agent of BPB of rice . Pathogenicity and virulence factors of B. glumae include the phytotoxin toxoflavin , lipase activity  and the motility driven by flagella . Type III effectors that are delivered into host cells by a type III secretion system play critical roles in bacterial pathogenicity and these effectors challenge host resistance mechanisms in many bacterial pathosystems [13, 14]. Thus, mutation in any of the genes encoding a type III secretion system could significantly attenuate the virulence of B. glumae . On the host side, extensive research has been performed on plant resistance against rice blast (caused by the fungus Magnaporthe grisea) and bacterial blight (caused by the bacterium Xanthomonas oryzae pv. oryzae) [16–35]. Although weakly virulent strains of X. oryzae pv. oryzae used to be found in 1989  bacterial blight has not been identified in the U.S. since then . Because it is caused by a different pathogen, we believe BPB needs to be investigated especially the resistance mechanism against it.
The goal of this research was to understand the rice – B. glumae pathosystem using transcriptomics via next-generation sequencing and bioinformatics. Here, we present data that the mechanism for resistance against BPB is independent of resistance mechanisms against other known rice diseases and other novel findings related to this pathosystem.
Responses of moderately resistant and susceptible rice varieties to inoculation with B. glumae
Illumina libraries and preliminary analysis of Illumina sequence reads
Number of transcripts expressed at each sample point
Differentially expressed transcripts upon inoculation of rice by B. glumae
Number of differentially expressed transcripts between conditions tested
We checked for the expression of known disease resistance genes against blast and bacterial blight in rice such as Pi-2, Pi-36, Pi5-1, Pi9, Pib, Pi-d2, Pikm1-TS, Pikm2-TS, Pikp-2, Pit, Pita, Piz-t, Xa21[26, 27], Xa26[28, 29], Xa1, Xa27, xa5; and xa13[33–35]. Of these, one putative xa1 transcript (LOC_Os02g16330) was detected in the control S and in the inoculated R conditions, while another (LOC_Os02g16260) was shown in the control R. In these three cases, the levels were very low, implying background levels. The rest of the other genes were not expressed in the inoculated R, which suggests that the resistance genes against blast and bacterial blight are not involved in BPB. Further, this suggests that the resistance pathways in rice against B. glumae are different from blast and bacterial blight. Of note, these known genes were not expressed in the control R, control or inoculated S either.
Gene Ontology annotation of differentially expressed transcripts
Validation of selected differentially expressed transcripts
Rice is one of the most studied plants for several reasons, one of which is its economic significance. Major pathogens to rice include Magnaporthe grisea, the causative agent of rice blast and Xanthomonas oryzae pv. oryzae, the causative agent of bacterial blight. Because these diseases inflict rice, extensive research has been done to understand host-pathogen interactions and develop disease management strategies. As a result, most rice resistance genes have been cloned against these diseases (Results Section). In this study, we embarked to understand BPB, a rice disease caused by another pathogen, B. glumae. To obtain a broad representation of the interaction between rice and B. glumae, we used RNA-Seq to identify transcripts that were differentially expressed between a resistant and a susceptible interactions. With young stems from seedlings as test tissues, we observed symptom development after 24, 48, 72 and 96 hours of inoculation. Stems over panicles were chosen for better distinction of symptoms between a susceptible and a resistant interaction over a short period of time. Reports have shown that symptoms on the panicles usually appear after two weeks of inoculation [3, 4]. At this time period, gene expression may not be indicative of what happens early in the interaction, but earlier time points may not show significant difference in the symptoms in the panicles. In the stem tissues, browning and lesion formation progressed at a remarkable rate in the susceptible genotype, especially when compared to the control and the resistant interaction. However, we found that after 48 hours was the earliest time point that showed the most difference in the responses between the two genotypes. This was our basis for choosing 48 hours post inoculation for the RNA-Seq experiment. We opted for replicates over deeper sequence coverage to provide statistical measures on comparative analysis between any two sample points as we were ultimately interested in possible sources of resistance genes from rice. Haas et al. expressed the same sentiment when they said that for some systems, the tradeoff for having replicates rather than sequence depth will provide better biological insight and statistical confidence. We had three biological replications and each replicate was loaded into two lanes. We pooled the reads from two lanes to constitute the read counts for a specific replicate (multiplexed with other samples). The application of the appropriate bioinformatic tools was equally necessary. In this work, we used Bowtie and rpkmforgenes.py script using default parameters to pre-analyze our data, and DESeq to calculate for the differentially expressed transcripts. High-throughput data analysis requires accurate prediction of variability within the dynamic range of values and a suitable error model and DESeq attempts to achieve them by using the negative binomial distribution with the variance and mean linked by local regression. DESeq was preferred over other programs because it provides better statistics for high-throughput data with few numbers of replicates such as RNA-Seq and it addresses the issue of data normalization in a more robust way compared to other available systems . As an initial filter after alignment into the rice genome, we eliminated reads that mapped to more than one locus, rRNAs and other repetitive sequences such as transposable elements. DESeq calculates padj (p-values adjusted for false discovery rate) to correct for multiple testing. We applied an FDR cutoff of 5% and used fold change values of 2 or greater to select for differentially expressed transcripts in the host (rice). These stringent conditions nevertheless generated sufficient transcripts for analysis.
The rice genes that were differentially expressed between the resistant and susceptible interaction were disease resistant-types or related, different enzymes, transcription factors, expressed and hypothetical proteins as well as proteins of unknown function. An enumeration of the transcript-types did not show significant distinction between the two types of interactions. However, a GO analysis, which classifies genes into biological process, molecular function and cellular component, demonstrated a clear distinction between them. In the R interaction, we saw an up-regulation of transcripts that are enriched for defense response, programmed cell death and generalized cell death transcripts under biological process; ATP, nucleoside, nucleotide and protein binding transcripts under molecular function; and mitochondrion supporting transcripts for cell component. By contrast, the S interaction displayed an up-regulation of transcripts enriched in lipid metabolic process, defense response and programmed cell death under biological processes; signal transducer, receptor activities and carbohydrate binding under molecular functions; and plasma membrane and membrane parts under cell component. This side by side comparison of ontologies presented that although disease resistance transcripts and most likely proteins were also up-regulated in the S interaction, other constituents appeared to play important roles in the resistance mechanism. For example, transcripts supporting molecular function and cell component were different between the two interaction types. Of note, a QTL for bacterial seedling rot, another rice disease caused by B. glumae is qRBS1, was mapped in chromosome 10 . This locus was not differentially expressed in any of the comparisons made, suggesting different pathways for resistance in these two rice diseases even if they were caused by the same pathogen. When we looked closer into specific disease-related transcripts, none of the previously cloned rice genes were differentially expressed, suggesting the known resistance genes were not involved in this interaction. Most were not even expressed at the tested conditions. The differentially expressed transcripts may represent genes that are unique to the rice - B. glumae interaction, indicating that resistance of rice to BPB may be conferred by a different set of genes and their roles in the interaction need to be further investigated. Of characterized resistance-related transcripts, those of the NBS-LRR-type and in some cases, sub-families such as NB-ARC and RPM1, were found to be both up-regulated and down-regulated in the R vs S comparison (down-regulation means up-regulation in the S vs R comparison), though actual transcripts associated with each group were not shared. We selected a few of these transcripts to verify their expression using qRT-PCR. All six disease-related transcripts (four RPM1 and two NBS-LRR types) were co-expressed in both the control and inoculated R conditions but not expressed at all in the control or inoculated S. A PIF-like ORF1 that is mapped in chromosome 8 follows the same trend. These results suggest that the R genotype maybe keeping a constitutive level of resistance arsenal to help it combat future B. glumae attacks.
It has been documented that the nucleotide-binding site or NBS (also NB-ARC) is a conserved domain for ATP binding and hydrolysis and sequences at the amino terminus are required for protein-protein interaction [48, 49]. The leucine rich repeats (LRRs) vary in number and the amino terminal domain seems to regulate activation while the carboxy terminal domain appears to function in recognition . It implies then that where NBS-LRR resistance genes are involved, so does specificity in the interaction. More so, it also suggests the involvement of an effector protein, which initiates the cascade of events that will lead to resistance. RPM1, a type of NBS-LRR resistance gene was originally cloned from Arabidopsis in response to the bacterium Pseudomonas syringae. Prior studies had shown that the NBS-LRR gene family is constitutively negatively regulated [63–66] and gets activated in the presence of pathogens through a mechanism that is not clearly understood. However, their activation needs to be precise (in space and time) for resistance to ensue. In addition, earlier studies have demonstrated indirectly that the NBS motif binds to ATP or GTP for activation [48, 49, 64]. If this gene family functions in the same manner in this pathosystem, then processes necessary for their activation should be up-regulated or activated as well. Our GO annotation results suggest an enrichment of ATP binding activities under molecular function, supporting the premise that an activation of NBS type motifs occurs in the R genotype during B. glumae challenge. No evidence for NBS activation was shown in the inoculated S genotype, despite the up-regulation of this type motif.
The rice genome has been annotated with more than 500 NBS-LRR-type genes although more than a hundred were predicted to be pseudogenes [46, 47]. Available literature shows that they cluster where mapped  and high sequence diversity exists in both the NBS and the LRR domains . It has been proposed that this gene family arose by several independent events of gene duplications all throughout rice evolution . Prior research also demonstrated that diversifying selection has shaped the evolution of the family, giving rise to the diversity that has been observed among its members . Hence, it is conceivable that of those members that are functional, the mechanisms that they provide may not necessarily be similar. Our results showed a clustering of up-regulated transcripts in chromosomes 11 and 8 in the R vs S comparison. Although previous work showed a bias clustering of disease-related genes in chromosomes 11 and 12 , our results suggested that resistance against BPB was not a direct result of clustering alone as none of the resistance genes in chromosome 12 were differentially expressed between the conditions tested. Furthermore, there are other loci where resistance genes are clustered, though to a lesser degree . Our results suggest that rice may have utilized the clusters of resistance genes together with another factor/s to devise a resistance mechanism against BPB. These factors may include NBS-LRR activation partners like ATP binding. Another interesting finding was the co-expression in the water- and pathogen-inoculated states, and probably, constitutive expression of a PIF-like ORF1 transcript in the R and almost none in the S genotype. Because the filtering method for the reads that we used involved elimination of those mapped to more than one locus and known repeats like rRNA and transposable elements, we can state that the reads only mapped to the PIF-like ORF1 transcript in the genotypes that we studied. When we further investigated it from the Rice Genome Annotation Project , it showed only one match (located in chromosome 8), suggesting that only one copy exists per haploid genome. P instability factor or PIF[52–56], a family of Class 2 transposable element is widely distributed in plants and other metazoans [55, 57, 58]. Jiao and Deng  performed a genome-wide survey of transcriptional activity of transposable element-related genes in 15 developmental stages and stress conditions in rice and found no expression of PIF-like transcripts in their test plants, suggesting that the PIF-like ORF1 is not expressed in all rice genotypes. PIF has two open reading frames, ORF1 and ORF2, of which ORF2 is most likely the transposase TPase . The function of ORF1 is still unknown, but its predicted protein sequence has significant homology to the Myb/SANT domain. The Myb domain is involved in DNA binding  while the SANT domain, although shares a strong homology with Myb sequences, is involved in protein-protein interactions . When we searched for homologs of the transposase in our transcriptomes, we found out that the reads mapped to several loci and were eliminated from the analysis. Whether the homologs were truly repetitive or this was an artifact of high-throughput sequencing analysis remains to be explored and is beyond the scope of this study. Based on the GO result that the R genotype did not show significant (p ≤ 0.05) enrichment for signal transduction, it appears that the PIF-like ORF1 may have been recruited to behave as a transcriptional regulator through DNA binding and not as a participant in signal transduction processes in this pathosystem. The transcripts that were constitutively expressed in the R genotype as quantified using qRT-PCR were all mapped in chromosome 11, suggesting that the PIF-like ORF1 (chromosome 8) may be acting in trans on the genes that it regulates.
All things considered, we propose a resistance mechanism in rice against BPB that existed early in rice domestication and that is not shared with other diseases including rice blast and bacterial blight. This was supported by the recent occurrence of this disease  and several observations that we noted in this work that are linked to resistance. We propose that shortly before it is domesticated, encounters between rice and B. glumae are limited. The genome of rice along with the prevailing environment at that time may have supported resistance. Specifically, the cluster of resistance genes that include NBS-LRR and related types in chromosomes 11 and 8, the up-regulation of the PIF-like ORF1 and the enrichment for ATP binding all contribute to this resistance. Because they are available, rice may have co-opted them as resistance contributors against BPB. The involvement of NBS-LRR-type transcripts and activation partner ATP binding suggests that the resistance mechanism consists of an effector molecule, probably from the pathogen, that is recognized by the host. The effector activates a cascade of events that will eventually lead to resistance in the host. It is possible that the PIF-like ORF1 may have been recruited to participate in the activation of the NBS-LRR genes. However, changes in global weather patterns, specifically gradual warming, favored the breaking of the resistance originally held by wild rice species. This is not outrageous as an increase in new or previously insignificant plant diseases caused by pathogens that grow optimally at higher temperatures has been observed with the increase in global temperatures [1, 71]. B. glumae is one of these pathogens.
Alternatively, because it may have been a part of a DNA transposon, the PIF-like ORF1 may be performing a more active role in the resistance pathway. This remains to be tested but is not a part of this study. Furthermore, we do not exclude that other processes may be occurring in parallel. The list of differentially expressed transcripts includes proteins of unknown functions and other disease related proteins. Their roles in the resistance pathway need to be uncovered in order to paint a complete picture of the resistance mechanism.
The main objective of this study was to understand the interaction between rice and B. glumae, the bacterium that causes bacterial panicle blight, using transcriptomics via next-generation sequencing technology and bioinformatics. This is a timely study of BPB as it is an emerging disease in the rice growing regions of the United States. With our strategy, we were able to provide a model for the resistance mechanism and present hypotheses to be tested. Of note is the hypothesis that resistance existed just prior to rice domestication suggesting that sources of resistance may be found in wild species. A good example is Xa21, a resistance gene that confers resistance to rice bacterial blight (caused by the bacterium X. oryzae pv. oryzae) was originally identified in the wild rice relative Oryza longistaminata and introgressed into Oryza sativa. We have generated candidate loci that may play major roles in conferring resistance against BPB. Along with phenotypic studies on the response of wild species against BPB, these loci may be used as molecular markers as well as foundation to build evolutionary history studies of this disease. However, in order to complete the model, functional assays need to be performed on these candidate loci/genes and similar studies with time series and on other developmental stages of rice especially the panicle stage need to be carried out.
A diagram of the Materials and Methods is given in Additional file 4.
Generation of samples for RNA-Seq and qRT-PCR
B. glumae strain 189gr-4, requested from Dr. Jong Hyun Ham, Louisiana State University, was used to inoculate rice seedlings. Originally isolated from Texas, it is a wild type strain with confirmed pathogenicity on rice seedlings and panicles . To prepare the inoculum, the bacterium was grown on NBY plates , incubated at 30°C for 24 h, harvested with a sterile cotton swab and suspended in a vial containing sterile distilled water to OD420 of 0.3 (108 CFU/mL) .
Clearfield rice varieties CL 151 and CL 161 were used for this study . Five to six seedlings were planted in individual pots grown in a greenhouse at Mississippi State University. Each seedling was injected with 0.5 mL of inoculum into the main stem and greenhouse conditions were maintained at 35-40°C during the day with >75% relative humidity. Water-inoculated plants were used as negative control. Three biological replications were completed for each sample point. Symptom development was observed 24, 48, 72 and 96 hours after inoculation. Inoculated parts were cut 1.0 cm above and below the inoculation point and frozen immediately in liquid nitrogen. Total RNA was extracted from tissues using a modified hot borate method  or Trizol reagent (Life Technologies, Carlsbad CA).
Sequencing using Illumina Hi-Seq
Total RNA was treated with DNAse I (Promega Madison, WI) and subsequently cleaned using the RNeasy Mini Kit (Qiagen, Valencia, CA) as described by the manufacturer. mRNA were selected using primers designed for polyA tails. The selected mRNA was subjected to chemical fragmentation that resulted in sizes ranging from 100–300 bp. cDNA construction and attachment of unique barcodes to each sample ensued (Illumina Inc, San Diego CA). The sequencing libraries with covalently attached unique barcodes were constructed from the cDNA using Illumina chemistry. After RT-PCR, a portion of each library (35ul) was loaded onto a Pippin Size Selection System gel cassette (3%) that was programmed to select the final sequencing library (~148 bp). Sample was eluted in 40ul and concentrated to 10ul using AmpureXP beads (Agencourt). Bioanalyzer, Qubit fluorometer and qPCR were used to quantify them. The sequencing pools were composed of twelve individually bar coded samples. Equimolar amounts (10 μl of 10 nM solution) from each library sample were pooled to constitute the samples that were loaded into a lane in a flowcell. Concentrations measured by using qPCR, which is more sensitive, were used as basis in pooling the samples. Each biological replicate of the samples was loaded into two different Hi-Seq lanes, together with eleven other biological replicates. Sequences of 1 × 50 bp single reads were generated. The raw reads that passed quality control were submitted to NCBI (SRA SRP015433 BioProject: PRJNA174463).
High-throughput sequence analysis
After the reads were grouped into their specific sample points, adaptors and indices were trimmed. Sequence reads from rice were aligned using Bowtie  against the most recent version of the Rice Genome available from the Rice Genome Annotation Project  using default parameters. The alignment quality threshold is not more than two mismatches within the first 28 bases of the read and the sum of the Phred quality at all mismatches may not exceed 70. Read count transcript expression values of every annotated gene region were calculated using the rpkmforgenes.py script . This created a list of gene IDs, loci where they were mapped and their corresponding read counts for each replicate of each sample point. Using gene IDs as reference, preliminary analysis was performed to eliminate those with zero read count values and those mapped to more than one locus, transposable elements and rRNA. This initial result was used as basis for all succeeding analyses. Pairwise analyses between genotypes within an inoculation state and between inoculation states within a genotype were performed to obtain differential gene expression using DESeq . DESeq utilized read counts from the rpkmforgenes.py output and normalized them initially for the analysis. Consequently, this step created read count values that were not whole numbers. For any two conditions being compared, DESeq accepted the read counts from all replicates from these conditions and processed the read counts to derive the fold change values between the two conditions. The differential expression analysis was performed on gene IDs that mapped to one and only one locus. To determine differential expression, significance was set to a padj value of at most 0.05 and a fold change value of 2 was used as cut-off score. Gene Ontology (GO) annotations for the differentially expressed transcripts were identified by searching the primary GO databases [75–77] using the associated GO retrieving tools [78, 79]. Furthermore, GO enrichment for the differentially expressed transcripts was analyzed using agriGO enrichment tools . Default parameters were used in all cases.
Validation of differential gene expression using qRT-PCR
Primers used in this study
Primer sequence (5′ → 3′)
18S rRNA-1 (forward)
18S rRNA-2 (reverse)
This research did not involve any human or animal subjects, materials or data and therefore did not require any ethics oversight or approval in these respects.
Bacterial panicle blight
Resistant genotype or CL 161
Susceptible genotype or CL 151
RNA integrity number
Flse discovery rate
The authors thank Mary Rougeau and Ryan Chen for their technical assistance. This research was funded by the MAFES - SRI research grant to ZVM (USDA-NIFA, MIS- 241030) and SL (USDA-NIFA, MIS-401170) and to DGP from the US Department of Agriculture – NIFA (Award Number 2010-34609-21206 10080747) and the US Department of Agriculture – ARS (Award Number 58-6402-1-645 12030374).
- Ham JH, Melanson RA, Rush MC: Burkholderia glumae: next major pathogen of rice?. Mol Plant Pathol. 2011, 12: 329-339. 10.1111/j.1364-3703.2010.00676.x.PubMedView ArticleGoogle Scholar
- Nandakumar R, Bollich P, Groth D, Rush MC: Confirmation of the partial resistance of Jupiter rice to bacterial panicle blight caused by Burkholderia glumae through reduced disease and yield loss in inoculated field tests. Phytopathology. 2007, 97: S82-S83.Google Scholar
- Nandakumar R, Shahjahan AKM, Yuan XL, Dickstein ER, Groth DE, Clark CA, Cartwright RD, Rush MC: Burkholderia glumae and B. gladioli cause bacterial panicle blight in rice in the southern United States. Plant Dis. 2009, 93: 896-905. 10.1094/PDIS-93-9-0896.View ArticleGoogle Scholar
- Sayler RJ, Cartwright RD, Yang YN: Genetic characterization and real-time PCR detection of Burkholderia glumae, a newly emerging bacterial pathogen of rice in the United States. Plant Dis. 2006, 90: 603-610. 10.1094/PD-90-0603.View ArticleGoogle Scholar
- Goto K, Ohata K: New bacterial disease of rice (brown stripe and grain rot). Annual Review of Phytopathological Society of Japan. 1956, 21: 46-47.Google Scholar
- Shahjahan AKM, Rush MC, Groth D, Clark C: Panicle blight. Rice Journal. 2000, 15: 26-29.Google Scholar
- Wamishe Y: Bacterial Panicle Blight of Rice. [https://www.agriskmanagementforum.org/content/bacterial-panicle-blight-rice]
- Lu S, Allen T: Causal Agents of Bacterial Panicle Blight of Rice and Evaluation of Disease Resistance of Rice Cultivars in Mississippi. Proceedings of Thirty-Four Rice Technical Working Group Meeting: 27February-1 March 2012. 2012, Louisianan State University Agricultural Center and Louisianan Agricultural Experiment Station: Hot Springs, AR. by Salassi ME, Wilson Jr CE, Walker TW, 7-Google Scholar
- Stockwell VO, Duffy B: Use of antibiotics in plant agriculture. Rev Sci Tech. 2012, 31: 199-210.PubMedGoogle Scholar
- Kim J, Kim J, Kang Y, Jang J, Jog G, Lim J, Kim S, Suga H, Nagamatsu T, Hwang I: Quorum sensing and the LysR-type transcriptional activator ToxR regulate toxoflavin biosynthesis and transport in Burkholderia glumae. Mol Microbiol. 2004, 54 (4): 921-934. 10.1111/j.1365-2958.2004.04338.x.PubMedView ArticleGoogle Scholar
- Devescovi G, Bigirimana J, Degrassi G, Cabrio L, LiPuma JJ, Kim J, Hwang I, Venturi V: Involvement of a quorum-sensing-regulated lipase secreted by a clinical isolate of Burkholderia glumae in severe disease symptoms in rice. Appl Environ Microbiol. 2007, 73 (15): 4950-4958. 10.1128/AEM.00105-07.PubMed CentralPubMedView ArticleGoogle Scholar
- Kim J, Kang Y, Choi O, Jeong Y, Jeong JE, Lim JY, Kim M, Moon JS, Suga H, Hwang I: Regulation of polar flagellum genes is mediated by quorum sensing and FlhDC in Burkholderia glumae. Mol Microbiol. 2007, 64 (1): 165-179. 10.1111/j.1365-2958.2007.05646.x.PubMedView ArticleGoogle Scholar
- Buttner D, Bonas U: Getting across - bacterial type III effector proteins on their way to the plant cell. EMBO J. 2002, 21: 5313-5322. 10.1093/emboj/cdf536.PubMed CentralPubMedView ArticleGoogle Scholar
- Cornelis GR: Type III secretion: a bacterial device for close combat with cells of their eukaryotic host. Philos Trans R Soc Lond B. 2000, 355: 681-693. 10.1098/rstb.2000.0608.View ArticleGoogle Scholar
- Kang Y, Kim J, Kim S, Kim H, Lim JY, Kim M, Kwak J, Moon JS, Hwang I: Proteomic analysis of the proteins regulated by HrpB from the plant pathogenic bacterium Burkholderia glumae. Proteomics. 2008, 8: 106-121. 10.1002/pmic.200700244.PubMedView ArticleGoogle Scholar
- Zhou B, Qu S, Liu G, Dolan M, Sakai H, Lu G, Bellizzi M, Wang G: The eight amino-acid differences within three leucine-rich repeats between Pi2 and Piz-t resistance proteins determine the resistance specificity to magnaporthe grisea. Mol Plant Microbe Interact. 2006, 19 (11): 1216-1228. 10.1094/MPMI-19-1216.PubMedView ArticleGoogle Scholar
- Liu X, Lin F, Wang L, Pan Q: The in silico map-based cloning of Pi36, a rice coiled-coil nucleotide-binding site leucine-rich repeat gene that confers race-specific resistance to the blast fungus. Genetics. 2007, 176 (4): 2541-2549. 10.1534/genetics.107.075465.PubMed CentralPubMedView ArticleGoogle Scholar
- Lee S, Song M, Seo Y, Kim H, Ko S, Cao P, Suh J, Yi G, Roh J, Lee S, An G, Hahn T, Wang G, Ronald P, Jeon J: Rice Pi5-mediated resistance to Magnaporthe oryzae requires the presence of two coiled-coil-nucleotide-binding-leucine-rich repeat genes. Genetics. 2009, 181 (4): 1627-1638. 10.1534/genetics.108.099226.PubMed CentralPubMedView ArticleGoogle Scholar
- Qu S, Liu G, Zhou B, Bellizzi M, Zeng L, Dai L, Han B, Wang G-L: The broad-spectrum blast resistance gene Pi9 encodes a nucleotide-binding site–leucine-rich repeat protein and is a member of a multigene family in rice. Genetics. 2006, 172 (3): 1901-1914.PubMed CentralPubMedView ArticleGoogle Scholar
- Wang Z, Yano M, Yamanouchi U, Iwamoto M, Monna L, Hayasaka H, Katayose Y, Sasaki T: The Pib gene for rice blast resistance belongs to the nucleotide binding and leucine-rich repeat class of plant disease resistance genes. Plant J. 1999, 19 (1): 55-64. 10.1046/j.1365-313X.1999.00498.x.PubMedView ArticleGoogle Scholar
- Chen X, Shang J, Chen D, Lei C, Zou Y, Zhai W, Liu G, Xu J, Ling Z, Cao G, Ma B, Wang Y, Zhao X, Li S, Zhu L: A B-lectin receptor kinase gene conferring rice blast resistance. Plant J. 2006, 46 (5): 794-804. 10.1111/j.1365-313X.2006.02739.x.PubMedView ArticleGoogle Scholar
- Ashikawa I, Hayashi N, Yamane H, Kanamori H, Wu J, Matsumoto T, Ono K, Yano M: Two adjacent nucleotide-binding site-leucine-rich repeat class genes are required to confer Pikm-specific rice blast resistance. Genetics. 2008, 180 (4): 2267-2276. 10.1534/genetics.108.095034.PubMed CentralPubMedView ArticleGoogle Scholar
- Yuan B, Zhai C, Wang W, Zeng X, Xu X, Hu H, Lin F, Wang L, Pan Q: The Pik-p resistance to Magnaporthe oryzae in rice is mediated by a pair of closely linked CC-NBS-LRR genes. Theor Appl Genet. 2011, 122 (5): 1017-1028. 10.1007/s00122-010-1506-3.PubMedView ArticleGoogle Scholar
- Hayashi K, Yoshida H: Refunctionalization of the ancient rice blast disease resistance gene Pit by the recruitment of a retrotransposon as a promoter. Plant J. 2009, 57 (3): 413-425. 10.1111/j.1365-313X.2008.03694.x.PubMedView ArticleGoogle Scholar
- Bryan G, Wu K-S, Farrall L, Jia Y, Hershey H, McAdams S, Faulk K, Donaldson G, Tarchini R, Valent B: A single amino acid difference distinguishes resistant and susceptible alleles of the rice blast resistance gene Pi-ta. Plant Cell. 2000, 12 (11): 2033-2045. 10.1105/tpc.12.11.2033.PubMed CentralPubMedView ArticleGoogle Scholar
- Song W, Wang G, Chen L, Kim H, Pi L, Holsten T, Gardner J, Wang B, Zhai W, Zhu L: A receptor kinase-like protein encoded by the rice disease resistance gene, Xa21. Science. 1995, 270: 1804-1806. 10.1126/science.270.5243.1804.PubMedView ArticleGoogle Scholar
- da Silva FG, Shen Y, Dardick C, Burdman S, Yadav RC, de Leon AL, Ronald PC: Bacterial genes involved in type I secretion and sulfation are required to elicit the rice Xa21-mediated innate immune response. Mol Plant Microbe Interact. 2004, 17: 593-601. 10.1094/MPMI.2004.17.6.593.PubMedView ArticleGoogle Scholar
- Chun H, Choi O, Goo E, Kim N, Kim H, Kang Y, Kim J, Moon JS, Hwang I: The quorum sensing-dependent gene katG of Burkholderia glumae is important for protection from visible light. J Bacteriol. 2009, 191 (13): 4152-4157. 10.1128/JB.00227-09.PubMed CentralPubMedView ArticleGoogle Scholar
- Sun X, Cao Y, Yang Z, Xu C, Li X, Wang S, Zhang Q: Xa26, a gene conferring resistance to Xanthomonas oryzae pv. oryzae in rice, encodes an LRR receptor kinase-like protein. Plant J. 2004, 37: 517-527. 10.1046/j.1365-313X.2003.01976.x.PubMedView ArticleGoogle Scholar
- Yoshimura S, Yamanouchi U, Katayose Y, Toki S, Wang Z, Kono I, Yano M, Iwata N, Sasaki T: Expression of Xa1, a bacterial blight-resistance gene in rice, is induced by bacterial inoculation. Proc Natl Acad Sci. 1998, 95: 1663-1668. 10.1073/pnas.95.4.1663.PubMed CentralPubMedView ArticleGoogle Scholar
- Gu K, Yang B, Tian D, Wu L, Wang D, Sreekala C, Yang F, Chu Z, Wang GL, White FF: R gene expression induced by a type-III effector triggers disease resistance in rice. Nature. 2005, 435: 1122-1125. 10.1038/nature03630.PubMedView ArticleGoogle Scholar
- Iyer AS, McCouch SR: The rice bacterial blight resistance gene xa5 encodes a novel form of disease resistance. Mol Plant Microbe Interact. 2004, 17: 1348-1354. 10.1094/MPMI.2004.17.12.1348.PubMedView ArticleGoogle Scholar
- Jiang GH, Xia ZH, Zhou YL, Wan J, Li DY, Chen RS, Zhai WX, Zhu LH: Testifying the rice bacterial blight resistance gene xa5 by genetic complementation and further analyzing xa5 (Xa5) in comparison with its homolog TFIIAgamma1. Mol Genet Genomics. 2006, 275: 354-366. 10.1007/s00438-005-0091-7.PubMedView ArticleGoogle Scholar
- Yang B, Sugio A, White FF: Os8N3 is a host disease-susceptibility gene for bacterial blight of rice. Proc Natl Acad Sci U S A. 2006, 103: 10503-10508. 10.1073/pnas.0604088103.PubMed CentralPubMedView ArticleGoogle Scholar
- Chu Z, Yuan M, Yao J, Ge X, Yuan B, Xu C, Li X, Fu B, Li Z, Bennetzen J, Zhang Q, Wang S: Promoter mutations of an essential gene for pollen development result in disease resistance in rice. Genes Dev. 2006, 20 (10): 1250-1255. 10.1101/gad.1416306.PubMed CentralPubMedView ArticleGoogle Scholar
- Jones RK, Barnes LW, Gonzales CF, Leach JE, Alvaerez AM, Benedict AA: Identification of Low-Virulence Strains of Xanthomonas campestries pv. oryzae from Rice in the United States. Phytopathology. 1989, 79 (9): 984-990. 10.1094/Phyto-79-984.View ArticleGoogle Scholar
- Corral R, Leach JE, Verdier V, Vera Cruz C: Recovery Plan for Xanthomonas Oryzae Causing Bacterial Blight and Bacterial Leaf Streak of Rice. 2013, http://www.lopdf.net/preview/02L5ox4SQ6qqX7O3y7eK1VAH_vqDEX7H3sZekKw_8Ew,/Recovery-Plan-for-USDA.html?query=Toward-the-Identification-of-Features-of-Effective,Google Scholar
- Groth D, Hollier C: Rice Diseases of Louisiana. 2009, [https://www.lsuagcenter.com/NR/rdonlyres/B80795F5-1ECB-4497-9079-D78A61287FEF/54719/pub3084RiceDiseasesLOWRES.pdf]Google Scholar
- Langmead B, Trapnel C, Pop M, Salzberg SL: Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009, 10 (3): R25-10.1186/gb-2009-10-3-r25.PubMed CentralPubMedView ArticleGoogle Scholar
- Ramsköld D, Wang ET, Burge CB, Sandberg R: An abundance of ubiquitously expressed genes revealed by tissue transcriptome sequence data. PLoS Comput Biol. 2009, 5 (12): e1000598-10.1371/journal.pcbi.1000598.PubMed CentralPubMedView ArticleGoogle Scholar
- Anders S, Huber W: Differential expression analysis for sequence count data. Genome Biol. 2010, 11: R106-10.1186/gb-2010-11-10-r106.PubMed CentralPubMedView ArticleGoogle Scholar
- Mizobuchi R, Sato H, Fukuoka S, Tsushima S, Imbe T, Yano M: Identification of qRBS1, a QTL involved in resistance to bacterial seedling rot in rice. Theor Appl Genet. 2013, 126: 2417-2425. 10.1007/s00122-013-2145-2.PubMed CentralPubMedView ArticleGoogle Scholar
- Dai L-Y, Liu X-L, Xiao Y-H, Wang G-L: Recent advances in cloning and characterization of disease resistance genes in rice. J Integr Plant Biol. 2007, 49 (1): 112-119. 10.1111/j.1744-7909.2006.00413.x.View ArticleGoogle Scholar
- Jones JDG: Plant disease resistance genes: structure, function and evolution. Curr Opin Biotechnol. 1996, 7: 155-160. 10.1016/S0958-1669(96)80006-1.View ArticleGoogle Scholar
- Meyers BC, Kaushik S, Nandety RS: Evolving disease resistance genes. Curr Opin Plant Biol. 2005, 8: 129-134. 10.1016/j.pbi.2005.01.002.PubMedView ArticleGoogle Scholar
- Monosi B, Wisser RJ, Pennill L, Hulbert SH: Full-genome analysis of resistance gene homologues in rice. Theor Appl Genet. 2004, 109: 1434-1447. 10.1007/s00122-004-1758-x.PubMedView ArticleGoogle Scholar
- Zhou T, Wang Y, Chen J-Q, Araki H, Jing Z, Jiang K, Shen J, Tian D: Genome-wide identification of NBS genes in japonica rice reveals significant expansion of divergent non-TIR NBS-LRR genes. Mol Gen Genomics. 2004, 271: 402-415. 10.1007/s00438-004-0990-z.View ArticleGoogle Scholar
- Belkhadir Y, Subramaniam R, Dangl JL: Plant disease resistance protein signaling: NBS–LRR proteins and their partners. Curr Opin Plant Biol. 2004, 7: 391-399. 10.1016/j.pbi.2004.05.009.PubMedView ArticleGoogle Scholar
- van Ooijen G, Mayr G, Kasiem MMA, Albrecht M, Cornelissen BJC, Takken FLW: Structure–function analysis of the NB-ARC domain of plant disease resistance proteins. J Exp Bot. 2008, 59 (6): 1383-1397. 10.1093/jxb/ern045.PubMedView ArticleGoogle Scholar
- Boyes DC, Nam J, Dangl JL: The Arabidopsis thaliana RPM1 disease resistance gene product is a peripheral plasma membrane protein that is degraded coincident with the hypersensitive response. Proc Natl Acad Sci U S A. 1998, 95: 15849-15854. 10.1073/pnas.95.26.15849.PubMed CentralPubMedView ArticleGoogle Scholar
- Kawasaki T, Nam J, Boyes DC, Holt BF, Hubert DA, Wiig A, Dangl JL: A duplicated pair of Arabidopsis RING-finger E3 ligases contribute to the RPM1- and RPS2-mediated hypersensitive response. Plant J. 2005, 44: 258-270. 10.1111/j.1365-313X.2005.02525.x.PubMedView ArticleGoogle Scholar
- Alexandrov NN, Brover VV, Freidin S, Troukhan ME, Tatarinova TV, Zhang H, Swaller TJ, Lu YP, Bouck J, Flavell RB, Feldmann KA: PIF-like orf1 [Zea mays]. GenBank. 2008, http://www.ncbi.nlm.nih.gov/protein/ACG44320.1,Google Scholar
- Alexandrov NN, Vyacheslav V, Brover VV, Freidin S, Troukhan ME, Tatarinova TV, Zhang H, Swaller TJ, Lu Y-P, Bouck J, Flavell RB, Feldmann KA: Insights into corn genes derived from large-scale cDNA sequencing. Plant Mol Biol. 2009, 69: 179-194. 10.1007/s11103-008-9415-4.PubMed CentralPubMedView ArticleGoogle Scholar
- Zhang X, Jiang N, Feschotte C, Wessler SR: PIF- and pong-like transposable elements: distribution, evolution and relationship with tourist-like miniature inverted-repeat transposable elements. Genetics. 2004, 166 (2): 971-986. 10.1534/genetics.166.2.971.PubMed CentralPubMedView ArticleGoogle Scholar
- Grzebelus D, Lasota S, Gambin T, Kucherov G, Gambin A: Diversity and structure of PIF/Harbinger-like elements in the genome of Medicago truncatula. BMC Genomics. 2007, 8: 409-10.1186/1471-2164-8-409.PubMed CentralPubMedView ArticleGoogle Scholar
- Jiang N, Feschotte C, Zhang X, Wessler SR: Using rice to understand the origin and amplification of miniature inverted repeat transposable elements (MITEs). Curr Opin Plant Biol. 2004, 7: 115-119. 10.1016/j.pbi.2004.01.004.PubMedView ArticleGoogle Scholar
- Casola C, Lawing AM, Betrán E, Feschotte C: PIF-like transposons are common in drosophila and have been repeatedly domesticated to generate New host genes. Mol Biol Evol. 2007, 24 (8): 1872-1888. 10.1093/molbev/msm116.PubMedView ArticleGoogle Scholar
- Jiao Y, Deng XW: A genome-wide transcriptional activity survey of rice transposable element-related genes. Genome Biol. 2007, 8: R28-10.1186/gb-2007-8-2-r28.PubMed CentralPubMedView ArticleGoogle Scholar
- International Rice Genome Sequencing Project: The map-based sequence of the rice genome. Nature. 2005, 436: 793-800. 10.1038/nature03895.View ArticleGoogle Scholar
- Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, Harris MA, Hill DP, Issel-Tarver L, Kasarskis A, Lewis S, Matese JC, Richardson JE, Ringwald M, Rubin GM, Sherlock G: Gene ontology: tool for the unification of biology. The Gene Ontology Consortium Nat Genet. 2000, 25 (1): 25-29.Google Scholar
- Haas BJ, Chin M, Nusbaum C, Birren BW, Livny J: How deep is deep enough for RNA-Seq profiling of bacterial transcriptomes?. BMC Genomics. 2012, 13: 734-10.1186/1471-2164-13-734.PubMed CentralPubMedView ArticleGoogle Scholar
- Dillies M-A, Rau A, Aubert J, Hennequet-Antier C, Jeanmougin M, Servant N, Keime C, Marot G, Castel D, Estelle J, Guernec G, Jagla B, Jouneau L, Laloe D, Le Gall C, Schaeffer B, Le Crom S, Guedj M, Jaffrezic F: A comprehensive evaluation of normalization methods for Illumina high-throughput RNA sequencing data analysis. Brief Bioinform. 2012, 14: 1-13.Google Scholar
- Bendahmane A, Farnham G, Moffett P, Baulcombe DC: Constitutive gain-of-function mutants in a nucleotide binding site–leucine rich repeat protein encoded at the Rx locus of potato. Plant J. 2002, 32: 195-204. 10.1046/j.1365-313X.2002.01413.x.PubMedView ArticleGoogle Scholar
- Shirano Y, Kachroo P, Shah J, Klessig DF: A gain-of-function mutation in an Arabidopsis toll Interleukin1 receptor–nucleotide binding site–leucine-rich repeat type R gene triggers defense responses and results in enhanced disease resistance. Plant Cell. 2002, 14: 3149-3162. 10.1105/tpc.005348.PubMed CentralPubMedView ArticleGoogle Scholar
- Hwang CF, Williamson VM: Leucine-rich repeat-mediated intramolecular interactions in nematode recognition and cell death signaling by the tomato resistance protein Mi. Plant J. 2003, 34: 585-593. 10.1046/j.1365-313X.2003.01749.x.PubMedView ArticleGoogle Scholar
- Zhang Y, Goritschnig S, Dong X, Li X: A gain-of-function mutation in a plant disease resistance gene leads to constitutive activation of downstream signal transduction pathways in suppressor of npr1-1, constitutive 1. Plant Cell. 2003, 11: 2636-2646.View ArticleGoogle Scholar
- The Rice Chromosomes 11 and 12 Sequencing Consortia: The sequence of rice chromosomes 11 and 12, rich in disease resistance genes and recent gene duplications. BMC Biol. 2005, 3: 20-PubMed CentralView ArticleGoogle Scholar
- Ouyang S, Zhu W, Hamilton J, Lin H, Campbell M, Childs K, Thibaud-Nissen F, Malek RL, Lee Y, Zheng L, Orvis J, Haas B, Wortman J, Buell CR: The TIGR rice genome annotation resource: improvements and new features. Nucleic Acids Res. 2007, 35: D883-D887. 10.1093/nar/gkl976.PubMed CentralPubMedView ArticleGoogle Scholar
- Klempnauer KH, Sippel AE: The highly conserved amino-terminal region of the protein encoded by the v-myb oncogene functions as a DNA-binding domain. EMBO J. 1987, 6 (9): 2719-2725.PubMed CentralPubMedGoogle Scholar
- Grüne T, Brzeski J, Eberharter A, Clapier CR, Corona DF, Becker PB, Müller CW: Crystal structure and functional analysis of a nucleosome recognition module of the remodeling factor ISWI. Mol Cell. 2003, 12 (2): 449-460. 10.1016/S1097-2765(03)00273-9.PubMedView ArticleGoogle Scholar
- Schaad NW: Emerging Plant Pathogenic Bacteria and Global Warming. Pseudomonas Syringae Pathovars and Related Pathogens—Identification, Epidemiology and Genomics. Edited by: Fatmi M, Collmer A, Iacobellis NS, Mansfield JW, Murillo J, Schaad NW, Ullrich M. 2008, New York, NY: Springer, 369-379.View ArticleGoogle Scholar
- Vidaver AK: Synthetic and complex media for the rapid detection of fluorescence of phytopathogenic pseudomonads: Effect of the carbon source. Appl Microbiol. 1967, 15: 1523-1524.PubMed CentralPubMedGoogle Scholar
- Buehring NW: CLEARFIELD® Rice. Mississippi State University Extension Service. 2008Google Scholar
- Wan CY, Wilkins TA: A modified hot borate method significantly enhances the yield of high-quality RNA from cotton (Gossypium hirsutum L.). Anal Biochem. 1994, 223 (1): 7-12. 10.1006/abio.1994.1538.PubMedView ArticleGoogle Scholar
- Barrell D, Dimmer E, Huntley RP, Binns D, O’Donovan C, Apweiler R: The GOA database in 2009–an integrated Gene Ontology Annotation resource. Nucleic Acids Res. 2009, 37 (Database issue): D396-D403.PubMed CentralPubMedView ArticleGoogle Scholar
- McCarthy FM, Gresham CR, Buza TJ, Chouvarine P, Pillai LR, Kumar R, Ozkan S, Wang H, Manda P, Arick T, Bridges SM, Burgess SC: AgBase: supporting functional modeling in agricultural organisms. Nucleic Acids Res. 2011, 39 (Database issue): D497-D506.PubMed CentralPubMedView ArticleGoogle Scholar
- Carbon S, Ireland A, Mungall CJ, Shu S, Marshall B, Lewis S: AmiGO: online access to ontology and annotation data. Bioinformatics. 2009, 25 (2): 288-289. 10.1093/bioinformatics/btn615.PubMed CentralPubMedView ArticleGoogle Scholar
- Binns D, Dimmer E, Huntley R, Barrell D, O’Donovan C, Apweiler R: QuickGO: a web-based tool for Gene Ontology searching. Bioinformatics. 2009, 25 (22): 3045-3046. 10.1093/bioinformatics/btp536.PubMed CentralPubMedView ArticleGoogle Scholar
- Du Z, Zhou X, Ling Y, Zhang Z, Su Z: agriGO: a GO analysis toolkit for the agricultural community. Nucleic Acids Res. 2010, 38 (Web Server issue): W64-W70.PubMed CentralPubMedView ArticleGoogle Scholar
- Pfaffl MW: A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acid Res. 2001, 29: e45-10.1093/nar/29.9.e45.PubMed CentralPubMedView ArticleGoogle Scholar
- Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2-ΔΔCT method. Methods Mol Biol. 2001, 25: 402-408.Google 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.