Two virulent sRNAs identified by genomic sequencing target the type III secretion system in rice bacterial blight pathogen
BMC Plant Biology volume 18, Article number: 237 (2018)
Small non-coding RNA (sRNA) short sequences regulate various biological processes in all organisms, including bacteria that are animal or plant pathogens. Virulent or pathogenicity-associated sRNAs have been increasingly elucidated in animal pathogens but little is known about similar category of sRNAs in plant-pathogenic bacteria. This is particularly true regarding rice bacterial blight pathogen Xanthomonas oryzae pathovar oryzae (Xoo) as studies on the virulent role of Xoo sRNAs is very limited at present.
The number and genomic distribution of sRNAs in Xoo were determined by bioinformatics analysis based on high throughput sequencing (sRNA-Seq) of the bacterial cultures from virulence-inducing and standard growth media, respectively. A total of 601 sRNAs were identified in the Xoo genome and ten virulent sRNA candidates were screened out based on significant differences of their expression levels between the culture conditions. In addition, trans3287 and trans3288 were also selected as candidates due to high expression levels in both media. The differential expression of 12 sRNAs evidenced by the sRNA-Seq data was confirmed by a convincing quantitative method. Based on genetic analysis of Xoo ΔsRNA mutants generated by deletion of the 12 single sRNAs, trans217 and trans3287 were characterized as virulent sRNAs. They are essential not only for the formation of bacterial blight in a susceptible rice variety Nipponbare but also for the induction of hypersensitive response (HR) in nonhost plant tobacco. Xoo Δtrans217 and Δtrans3287 mutants fail to induce bacterial blight in Nipponbare and also fail to induce the HR in tobacco, whereas, genetic complementation restores both mutants to the wild type in the virulent performance and HR induction. Similar effects of gene knockout and complementation were found in the expression of hrpG and hrpX genes, which encode regulatory proteins of the type III secretion system. Consistently, secretion of a type III effector, PthXo1, is blocked in Δtrans217 or Δtrans3287 bacterial cultures but retrieved by genetic complementation to both mutants.
The genetic analysis characterizes trans217 and trans3287 as pathogenicity-associated sRNAs essential for the bacterial virulence on the susceptible rice variety and for the HR elicitation in the nonhost plant. The molecular evidence suggests that both virulent sRNAs regulate the bacterial virulence by targeting the type III secretion system.
Small non-coding RNAs (sRNAs) are characteristic of unique nucleotide (nt) sequences comprising 50–500 nt under most circumstances and constitute a universal group of post-transcriptional regulators for gene expression [1,2,3]. They function by base pairing with target mRNAs via limited and extended complementarity, which are used to distinguish two categories, namely cis-sRNAs and trans-sRNAs, respectively . In addition, many sRNAs require the RNA binding protein Hfq to maintain stability and co-regulate target transcripts [4, 5]. Because of the primary function in post-transcription of target genes, sRNAs are able to regulate numerous bioprocesses in all living organisms [2, 3, 6]. In particular, animal and plant pathogens including bacteria employ certain sRNAs to regulate pathogenicity or virulence on their hosts [3, 7,8,9].
Animal-pathogenic bacteria in the Escherichia and Salmonella genera possess approximately a half of bacterial sRNAs as identified to date, but only a few have been associated with bacterial pathogenicity and/or virulence [10,11,12]. E. coli RyhB is a member of small antisense regulatory RNA family, with requirement for RNA chaperone Hfq to maintain stability and perform function [13,14,15]. RyhB works with the ferric uptake regulator to regulate cellular iron homeostasis, the clbA gene transcription, and the colibactin production . These responses affect pathogenicity not only in E. coli , but also in Shigella flexneri , S. dysenteriae , and Listeria monocytogenes . One more example is from the animal pathogen S. typhimurium, which has at least 280 sRNAs, as identified by deep sequencing of Hfq-bound transcripts . The sRNA InvS controls the bacterial invasion by coordinating the production of PrgH, a type III secretion apparatus protein, and FimZ, a negative regulator of invasion-related gene expression . In Staphylococcus aureus, moreover, RNAIII is a pathogenicity-associated sRNA reported as the first case [22, 25] and a secreted effector protein with multiple functions in the bacterial virulence [22,23,24]. While RNAIII directly targets hla, spa, and rot genes involved in quorum sensing [22, 24], its noncoding parts act as antisense RNAs to regulate translation and stability of related transcription factors, major virulence effectors, and cell wall metabolism enzymes [24, 25]. RNAIII also governs the expression of SA1000 mRNA, which hypothetically encodes a fibrinogen binding protein indispensable for bacterial adhesion to host surfaces and invasion of host tissues . These demonstrations suggest that bacterial RNAs use distinct mechanisms while cooperating with different functional partners to regulate pathogenicity and virulence.
In contrast to extensive studies on sRNAs in animal-pathogenic bacteria, little is known about plant pathogens regarding functions of sRNAs in relevance to virulence or pathogenicity . Limited information comes from major group of plant-pathogenic bacteria, especially Xanthomonas species and pathovar (pv.), with some sRNAs already identified in X. campestris pv. vesicatoria (Xcv), X. campestris pv. campestris (Xcc), X. oryzae pv. oryzae (Xoo) [26,27,28,29,30,31]. sX12 and sX13 from Xcv are the earliest and best studied sRNAs in the bacterial genus. Both sRNAs contribute to the bacterial virulence on host plant pepper [28, 32]. Deletion of sX13 not only eliminates Xcv virulence but also impinges the type III secretion system (T3SS) , which is ubiquitously present in Gram-negative bacteria and functions to deliver effectors en route to the host cytoplasm . Knockout of sX13 results in decreased expression of genes encoding structural protein HrcJ, regulators HrpG and HrpX, and translocator HrpF of the T3SS . In essence, HrpF serves as a type III translocator indispensable for effector translocation from bacterial cells into the cytosol of plant cells, where effectors execute their pathological functions . In agreement with the effect of sX13 on T3SS, sRNA-Xcc1 transcription is subject to positive control of HrpG and HrpX . Therefore, structural, regulatory, and secretory components of the T3SS are potential targets of sRNAs in Xcc.
It is unclear whether Xoo also deploys the functional mechanism of sRNAs targeting the T3SS due to paucity in related studies. Xoo is an important plant pathogenic bacterium, causing bacterial blight of rice, which is a major disease in China southern rice cultivation area and also in the southeast of Asia. The disease occurs mostly in rice leaves at the later stage of rice growth and is a severe threat to the high and stable yield of rice. Like in all plant pathogenic Gram-negative bacteria, Xoo uses the T3SS to secrete decades of effector proteins and translocate them into the cytosol of rice cells, where effectors execute their pathological roles [33, 36]. As a well-studied example, the Xoo strain PXO99A confers virulence to the susceptible rice variety Nipponbare through the function of PthXo1, a transcription activator-like (TAL) protein, in manipulating rice transcriptome . Evidently, an unappreciated functional relationship exists between certain sRNAs and effector secretion by the T3SS. However, this notion is not examined until now due to limited knowledge on sRNAs of plant pathogenic bacteria. To date, only eight sRNAs have been identified and shown to associate with various processes characterized by proteomics profiling analysis , but their pathological functions remain unclear. In conjecture, the genetic repertoire of Xoo is not likely to possess only eight sRNAs since hundreds of sRNAs exist in animal-pathogenic bacteria . Clearly, it is a great challenge to determine the number and biological functions of sRNAs in Xoo.
This study is devised to determine the repertoire and virulence-associated members of Xoo sRNAs by the genomic sRNA-Seq technique performed with Xoo bacterial cultures under standard growth conditions in contrast to the medium that induces virulence factor production. High-throughput sequencing and bioinformatics analysis identify a total of 601 sRNAs in the bacterial genome and 12 out of them alter expression levels with virulence induction in contrast to the standard growth condition. Genetic analysis of mutagenesis and complementation suggests that two sRNAs, trans217 and trans3287, are essential for the bacterial virulence. We present evidence that trans217 and trans3287 act on the T3SS and effector secretion, providing the first case study of Xoo sRNAs targeting the effector secretion mechanism.
Aggregate information on the genomic profile of sRNAs
The genomic cDNA libraries of Xoo used in sRNA-Seq were prepared from cultures of the international standard strain PXO99A grown on the standard growth medium polypeptone sucrose agar (PSA) or X. oryzae growth medium recipe 2 (XOM2), which induces gene expression related to the bacterial virulence . The sRNA cDNA libraries were normalized to the size ranging from 50 to 500 nt and analyzed by a high-throughput sequencing technique to yield sRNA reads . The number of sRNA reads from XOM2 was found to be less than that of PSA, suggesting that virulence induction impairs the expression of potential sRNAs in the bacterial genome. The number of sRNAs was determined as 601 in total, including 337 cis-sRNAs and 264 trans-sRNAs (Fig. 1; Table 1). Based on the different mechanisms by which cis- or trans-sRNAs function on targets , this sRNA classification information was used subsequently in the study to elucidate pathological functions of any of those Xoo sRNAs.
Identification of virulence-associated sRNA candidates
In the genomic profiling of Xoo sRNAs from PSA and XOM2 cultures, expression levels of 10 sRNAs were found to alter between both cultures and quantified as XOM2-to-PSA fold changes in base-2 logarithms (Fig. 2; Additional file 1: Table S1). These sRNAs were subjected to deletion mutagenesis and virulence test. Meanwhile, trans3287 and trans3288 with high counts under both culture conditions (Additional file 1: Table S1) were also selected in the future studies to test whether or not they associate with virulence. The profiling values of XOM2-to-PSA fold changes (Fig. 2) were utilized to estimate degrees by which the tested sRNAs possibly associate with the bacterial virulence, which is induced by the XOM2 medium . According to this criterion, trans198, trans217, and trans202, showing greater than 3-fold XOM2-to-PSA increases, were most positively relevant to virulence. Relatively less relevance to virulence was found in trans238, trans1513, and trans3288, which had XOM2-to-PSA fold changes of smaller than 2. The intermediate sRNAs were trans2192 and trans191, with fold increase values being 2.69 and 2.52, respectively. By contrast, four sRNAs were negatively associated with the bacterial virulence. In particular, the expression level of trans3747 was highly reduced by XOM2, showing the lowest XOM2-to-PSA fold decrease of 2 (Fig. 2).
These results from the sRNA-Seq profiling were confirmed by RT-qPCR, namely real-time quantitative reverse transcriptase polymerase chain reaction (RT-qPCR). The relative level of sRNA gene expression was quantified as the ratio of sRNA to 16S rRNA, a gene used as a reference owing to its characteristics of constitutive expression . As shown in Fig. 3, RT-qPCR data indicated that relative amounts of trans191, trans198, trans202, trans217, and trans238 transcripts from XOM2 cultures were increased by dozen to thousand times compared to those from PSA cultures. Several fold increases or decreases were detected for the other 7 sRNAs (Fig. 3), in agreement with Log2 XOM2-to-PSA fold changes found in the genomic profile (Fig. 2).
Both the sRNA-Seq and RT-qPCR data suggests that trans191, trans198, trans202, and trans217, and trans238 are highly responsive to the culture condition that induces gene expression related to bacterial virulence. These sRNAs are regarded as candidates of virulence-associated sRNAs.
Trans217 and trans3287 determine the bacterial virulence on rice
By NCBI blasting against adjacent sequences of published Xoo genes, length of the 12 sRNAs identified in this study varies from 291 to 487 nt and locate between or on known gene sequences at the bacterial chromosome (Fig. 4). This information was used to construct the bacterial sRNA-knockout mutants by dual-exchange of homologous arms.
For use in analyzing whether any of the 12 virulent sRNA candidates is certainly involved in the bacterial virulence on the susceptible variety Nipponbare, every sRNA sequence was deleted from the PXO99A genome by dual-exchange of homologous arms. Every sRNA deletion mutant was compared with the wild-type (WT) PXO99A strain in virulence assessments by inoculation experiments. Great variations were found between the different sRNAs with respect to the effects of gene deletion on the bacterial virulence (Additional file 2: Figure S1). Virulence was shown as induction of the bacterial blight symptom (Fig. 5a), the disease severities (Fig. 5b), and Xoo populations propagated in tissues of inoculated leaves before the symptom development (Fig. 5c). Based on the genetic analysis of 12 sRNAs, only trans217 and trans3287 had relevance to the bacterial virulence, as evidenced by the virulence level reduction by gene deletion, shown as Δtrans217 and Δtrans3287. All the other 10 sRNAs did not show evident relevance to the bacterial virulence tested in comparative experiments.
The virulence compromise was attributed to the gene knockout (Additional file 2: Figure S1), which was confirmed by PCR analysis of the bacterial genomic DNA. Consistently, transcripts of both trans217 and trans3287 were not detected in RT-qPCR analysis of RNAs from cultures either from PSA or from XOM2 in corresponding mutant (Fig. 6a, b). Both PCR and RT-qPCR analyses also conform gene deletion of virulence-unrelated sRNAs, such as trans191 (Fig. 6c). Moreover, trans217 or trans3287 gene deletion did not affect bacterial multiplication on the growth medium (Fig. 7). These data offer additional evidence that regulation of the bacterial virulence is a predominant function of trans217 and trans3287 during infection of the susceptible rice variety.
This notion was corroborated by genetic complementation experiments. The original sequences of trans217 and trans3287 cloned from the WT PXO99A genome was returned back into the genome of Δtrans217 and Δtrans3287, respectively. Resulting complement strains CΔtrans217 and CΔtrans3287 were restored to the WT not only in gene expression extents (Fig. 6), but also in virulence levels (Fig. 5a–c). In the parallel experiments, Δtrans191 was also complemented and tested similarly as for trans217 and trans3287 as trans191 is representative of virulence-unrelated Xoo sRNAs, with expression highly induced by XOM2 (Figs. 2 and 3) and without evident role in the bacterial virulence (Fig. 5a–c). The performances of trans191 in these experiments were constant, the virulence level of CΔtrans191 bacteria was similar to PXO99A. This was evidenced by equivalent severities of bacterial blight symptom (Fig. 5a), as well as equivalent quantities of blight lesion length (Fig. 5b) and bacterial population (Fig. 5c) in Nipponbare leaves irrespectively of inoculation with the different bacterial strains. These analyses strongly suggest that Xoo sRNAs showing high expression levels with the virulence induction are not necessarily related to the bacterial virulence but both trans217 and trans3287 definitely play a predominant role in the bacterial virulence.
Trans217 and trans3287 function for the HR induction in tobacco
In agreement with the original definition about the hypersensitive response and pathogenicity (hrp) gene cluster present in plant-pathogenic Gram-negative bacteria , the presence of an hrp cluster enables the WT PXO99A strain to induce the HR in tobacco N. benthamiana leaves . However, not all hrp genes contribute to the HR induction . Therefore, we investigated whether the HR induction involves trans217 and trans3287 by using the conventional method of infiltrating a bacterial suspension into the apoplastic space of tobacco leaves. We found that deleting trans217 or trans3287 nullified the bacterial capability of inducing the HR, whereas, genetic complementation of Δtrans217 and Δtrans3287 returned the HR induction (Fig. 8). Once again, the Δtrans191 or CΔtrans191 strain performed like the WT strain, eliciting a full development of the HR in tobacco leaves (Fig. 8). This result rules out the functional association of trans191 with the HR induction and suggests that both trans217 and trans3287 function in the way of certain hrp components to induce the HR.
Trans217 and trans3287 are needed for hrpG and hrpX expression
Once the hrpG and hrpX genes are expressed, their encoding products, regulatory proteins HrpG and HrpX, are essential for the subsequent expression of other hrp genes that encode structural and functional components of the T3SS machinery [32, 44]. To infer the functional relationship between the virulent sRNAs and hrp genes, we carried out RT-qPCR analysis to compare expression levels of hrpG and hrpX genes in PSA and XOM2 cultures of Xoo strains with and without trans191, trans217, and trans3287, respectively. Based on the ratios of sRNA to 16S rRNA transcript quantities, both hrpG and hrpX were expressed to equivalent degrees in the WT, Δtrans191, and CΔtran191 strains irrespectively of culture in PSA or XOM2 (Fig. 9a). However, significant decreases of hrpG and hrpX expression levels were caused by Δtrans217 or Δtrans3287 compared to the WT or complementary strains (Fig. 9b, c). Evidently, trans217 and trans3287 take part in hrpG and hrpX expression but trans191 does not.
Trans217 and trans3287 participate in PthXo1 secretion
Since HrpG and HrpX regulate the expression of other hrp genes that encode T3SS accessory and secreted proteins , trans217- and trans3287-dependent expression of hrpG and hrpX gene expression may have subsequent effect on T3SS effector secretion. To test the hypothesis, we analyzed the effects of both sRNAs on the PthXo1 secretion, which was monitored by using calmodulin-dependent adenylate cyclase (Cya), a eukaryotic cytoplasmic import marker . As shown in Fig. 10, PthXo1-Cya was highly secreted out of cells in the Xoo WT strain, whereas, the fusion protein secretion was totally cancelled in Δtrans217 and Δtrans3287 mutants, which were created by deletion of trans217 and trans3287 accordingly from the bacterial genome. By contrast, when the Δtrans217 and Δtrans3287 mutants were complemented, PthXo1 secretion out of bacterial cells of complementation CΔtrans217 and CΔtrans3287 strains was brought back to the WT level (Fig. 10). Clearly, trans217 and trans3287 are critical to PthXo1 secretion, providing the first case that connects bacterial sRNAs with effector secretion.
By the genomic sRNA-Seq profiling and experimental genetic analysis, one of the main findings in this study is the importance of trans217 and trans3287 to the virulence of Xoo strain PXO99A on the susceptible rice variety Nipponbare. We show that both sRNAs function as key components of the bacterial virulence. The nullification and restoration of the bacterial virulence correspondingly by deleting trans217 or trans3287 and by the genetic complementation suggest that trans217 and trans3287 are virulent sRNAs present in the Xoo genome. This information adds useful insight into the current understanding of Xoo sRNAs, which are paid attention only in recent years . To date, merely a few sRNAs have been characterized in the Xoo genome while none of them was definitely related to the bacterial pathogenicity or virulence .
To infer the molecular mechanism that underpins the pathological function of trans217 and trans3287, we have correlated the virulent role of both sRNAs with their positive effects on the expression of hrpG and hrpX genes. Both genes reduce expression quantities in the bacterial Δtrans217 and Δtrans3287 mutants but are highly expressed in the genetically complemented bacteria as in the WT strain. This result sheds a light into exploring molecular targets of trans217 and trans3287 while we do not have evidence at present to demonstrate how the sRNAs affects hrpG or hrpX expression. There is no homology present in trans217 and trans3287 sequences based on Blast Search. Analysis with TargetRNA2 (http://cs.wellesley.edu/~btjaden/TargetRNA2/) suggests that hrpG and hrpX are not targets of trans217 and trans3287. Instead, several predicted targets of trans217 and trans3287 are related to the T3SS apparatus formation or biochemical reactions. Thus, the effects of both sRNAs on hrpG and hrpX are indirect and involved mechanisms remain to be studied.
In the protein form, HrpG and HrpX are essential regulators of T3SS and execute the regulatory function at the transcription level [48, 49]. HrpG regulates the expression of hrpX gene while subsequently produced HrpX protein acts in turn to regulate the expression of other hrp genes, which encode T3SS accessory proteins like Hpa1 [50, 51] and effector proteins like PthXo1 secreted through the T3SS pathway . It was recently found that the Xcv sRNAs sX13 took part in the expression of hrpX in correlation with the virulent function of sX13 . Thus, targeting T3SS regulators is likely to be a common mechanism conserved in the virulence function of certain sRNAs from plant-pathogenic bacteria of the Xanthomonas genus. It is apparent that trans217 and trans3287 differ from sX13 in terms of direct molecular targets. Irrespective of the direct gene targets, however, the effect of sRNAs on T3SS components establishes a mechanistic connection between the sRNAs and bacterial virulence performance though particular effectors .
In agreement with the assumed functional connection, we demonstrate that trans217 and trans3287 are required for PthXo1 secretion from the bacterial cells outward to the culture environment. This finding indicates that targeting effector secretion may represent one of the molecular mechanisms by which bacterial sRNAs perform their virulent functions. An alternative mechanism has been characterized to be the suppression of immune responses in plants . Because PthXo1 determines the virulence of Xoo strain PXO99A on the susceptible rice variety Nipponbare , the apparent role of sRNAs in the control of PthXo1 secretion is obviously indirect but crucial for the effector to be translocated into the plant cell [33, 37]. Translocation from bacterial cells into the cytosol of plant cells is a key step for effectors to fulfil their pathological roles [33, 53, 54]. Effector translocation indispensably needs direct mediation by type III translocators, such as Hpa1 . Like HrpG and HrpX, type III translocators also belong to T3SS components but are closer than both Hrp proteins in the spatial and temporal patterns of functions from secretion to translocation of effectors . Therefore, sRNAs regulate bacterial virulence by complicated functional networks, which remain to be studied in the future. In the future, it is especially necessary to elucidate whether trans217 and trans3287 function during the Xoo-rice interaction process or during rice infection by the pathogen.
Based on the bacterial performance in the absence and presence of trans217 and trans3287, trans217 and trans3287 serve as pathogenicity-associated sRNAs essential for the bacterial virulence on the susceptible rice variety and for the HR elicitation in the nonhost plant. Gene expression and protein secretion data offer the molecular evidence suggesting that both virulent sRNAs regulate the bacterial virulence by targeting the type III secretion system.
Plant growth and bacteria cultures
Seeds of rice Oryza sativa L. japonica variety Nipponbare were sown in pots filled with a mixture of sand and potting soil (1:1 v/v) and seedlings were grown in a plant growth chamber under 25–26 °C, 85% humidity, and a 14-h light cycle at 250 μE/m2/s. Seeds of tobacco Nicotiana benthamiana were sown in the soil mixture as for rice and plants were grown in a different plant growth chamber with controlled environment as for rice growth. The WT Xoo strain PXO99A and different sRNA-related strains (Additional file 3: Table S2) are maintained in this lab. Bacteria were cultured in liquid nutrient broth (NB) medium in a 28 °C shaking incubator for preparation of inoculum while bacteria isolated from plant leaves were incubated in nutrient broth agar (NB) medium. Alternatively, bacteria were cultured separately on PSA and XOM2 media  for use in sRNA sequencing.
RNA isolation and quality optimization
Nursery bacterial PSA or XOM2 culture suspension with a desirable density (OD600 ≈ 0.5) was supplied to fresh medium and incubated in a 28 °C shaker for 6 h. Bacterial cells were collected by centrifuge and employed to isolate total RNA using the trizol reagent kit (Invitrogen) as per the manufacturer’s protocol. In order to eliminate possible variations between bacterial cultures, 2 biological repetitions were performed. In each repetition, RNA was isolated from a mixture of five PSA or XOM2 cultures and purity was assessed using the Nanodrop Spectrophotometer ND-1000 (Nanodrop Technologies). Each RNA sample had an A260:A280 ratio greater than 1.8 and A260:A230 ratio higher than 2.0. RNA integrity was evaluated using the Agilent 2200 TapeStation (Agilent) and each sample was quarantined to have the RIN (RNA Integrity Number) greater than 7.0. Then, RNAs were ligated with 3’-RNA adapter, followed by 3′-adapter ligation. Subsequently, the adapter-ligated RNAs were subjected to RT-PCR and amplified with a low-cycle program (10 cycles). To obtain 50–500 nt RNA, the PCR products were fractioned by PAGE according to instructions of NEBNext® Multiplex Small RNA Library Prep Set for Illumina® (Illumina). Finally, the purified library products were evaluated using the Agilent 2200 TapeStation and diluted to 10 pM for cluster generation in situ on the HiSeq2500 single-end flow cell, followed by sequencing (1 × 50 bp) with the HiSeq 2500 facility.
Xoo sRNA sequencing and sRNA identification
The PSA or XOM2-derived cDNA libraries were amplified by PCR and products were sequenced with the IlluminaHiSeq™2500 sequencer in a commercial company. Image data output from the sequencing device were converted into raw reads and stored in the FASTQ format. Clean reads were obtained after removing reads that contained adaptor sequences, reads containing more than 10% of unknown bases, and reads in which more than half of the quality values of the bases were less than 5. Sequences of the clear reads were mapped to the public PXO99A genome (NC_010717.2) and the overlap reads were selected to splice sRNA fragments if the coverage depth at each base was higher than 50X. This also yield the information on sRNA sequence, length and orientation. These sRNA fragments were named as sRNA1, sRNA2, and sequential codes for more RNAs . All sRNAs from the different samples (PSA and XOM2 cultures) were collectively analyzed to identify candidates of trans-sRNAs and cis-sRNAs according to their origin from the PXO99A genome. A sRNA shorter than 500 nt and transcribed from the opposite strand of an mRNA or another sRNA was identified as cis-sRNAs while a sRNA shorter than 500 nt and transcribed from the intergenic region was identified as a trans-sRNA [39, 55].
Verification of sRNAs differentially expressed in PSA and XOM2 cultures
Expression levels of sRNAs from both cultures quantified as average transcript levels as previously described . Log 2 fold change ratios of XOM2 to PSA transcript quantities were calculated by the edgeR method  and assessed for significant differences at P < 0.01 when the ratio was greater than 1. This type of sRNAs were analyzed by RT-qPCR to confirm expression levels in the Seq data, using specific primers (Additional file 4: Table S3). In every RT-qPCR protocol, the 25 μl reaction mixture was composed of 1 μl first-strand cDNA diluted 1:10, 2.5 μM primer, and 1 × SYBR Premix Ex Taq (TaKaRa). All reactions were performed in triplicate with null-template controls in which cDNA was absent. The constitutively expressed 16S rRNA gene was used as a reference, and relative expression level of an sRNA was quantified as the ratio of sRNA to 16S rRNA transcript amounts.
The sRNA-directed mutation and complementation
To generate sRNA deletion PXO99A mutants, flanking sequences of up- and downstream of sRNA were amplified by PCR from PXO99A using primers (Additional file 4: Table S3). PCR products were digested with respective restriction enzyme and ligated into the suicide vector pK18mabSacB (Additional file 3: Table S2). The recombinant vectors were transferred into PXO99A competent cells by electroporation. Transformed bacterial cells were incubated in sugar-absent NA plates containing 100 μg/mL kanamycin under 28 °C for 3 days and then single clones were verified by PCR. The suspension of positive clones was smeared on plates of NA medium containing 10% sucrose, incubated for 3 days, and subjected to PCR analysis to identify sRNA deletion mutants. Complemented strains were created by sRNA PCR fragment ligating into the expressive vector pHM1 via appropriate restriction enzyme (Additional file 4: Table S3). Every recombinant vector was transferred into the competent cells of the corresponding sRNA deletion mutant. Transformants were identified similarly as for mutants.
In plant inoculation assay, the leaves of two-month old rice were clipped using sterile scissors which were dipped in the bacterial cultures (OD600 ≈ 0.5) at a distance of about 2 cm from the leaves edge. Each strain was inoculated in five plants and three leaves of each plant were tested. Lesion length was scored 2 weeks after inoculation and the average was calculated. To determine bacterial population in leaf tissues, leaves were segmented, sterilized with 75% (v/v) ethanol, and homogenized in sterile water; bacteria were recovered from the resulting homogenates by culturing on NA medium. Bacterial number in leaves was gives as colony formation unit (cfu). At ultimate, virulence levels of different Xoo strains were evaluated by bacterial blight severity (lesion length) and bacterial populations in leaf tissues.
Bacterial growth assessment
Overnight Xoo strains PSA cultures (OD600 ≈ 1.0, 10 uL) were corresponding pipetted to a new bottle of 20 mL medium, cultured in a shaking incubator under 28 °C. Bacterial populations were measured every 2 h during 20–34 h.
The HR induction
The activated bacterial cultures were pipetted to NB medium until the optical density (OD600 ≈ 0.5). Bacterial suspensions of the different strains were separately infiltrated by using needleless syringes into intercellular spaces of at the 5th and 6th tobacco leaves. Water was applied similarly in control. Two days later, infiltrated leaves were photographed.
PthXo1 secretion assay
The recombinant PXO99A sRNA pthXo1-cya strains were grown under 28 °C in liquid NB with 50 μg/mL spectinomycin to logarithmic phase. Bacterial cells were harvested by centrifugation. The precipitated bacterial pellet was washed twice with sterile water and resuspended in 100 mL XOM2 liquid media  with 50 μg/mL spectinomycin to OD600 = 0.6 and cultured in a 28 °C shaker at 220 rpm for 16 h. XOM2 cultures were then separated into cell pellet and supernatant fractions by centrifugation. The proteins in the two fractions were extracted by sonication and by precipitation with 12.5% trichloroacetic acid, respectively . Proteins were separated by 6% SDS-PAGE and transferred to Immobilon-P membranes (Millipore) for immunoblotting analyses using a Cya antibody (Santa Cruz) or β-lactamase antibody (Abcam). The ampicillin resistance protein β-lactamase protein is encoded by the pHM1-cya vector, remains cell-bound unless non-specific cell leakage occurred, and was used as a control for nonspecific cell lysis. Protein blots were incubated with the specific antibody and hybridized to horseradish peroxidase-conjugated goat antimouse immunoglobulin G from the BeyoECL Plus kit (Beyotime).
All experiments were conducted at least three times with similar results. Quantitative data were analyzed by using the commercial IBM SPSS19.0 software package . Homogeneity-of-variance in data was determined by using the Levene test, and the formal distribution pattern of the data was confirmed by using the Kolmogorov-Smirnov test and P-P plots. Data were subjected to ANOVA, along with Fisher’s least significant difference test and Tukey-Kramer’s test, respectively. Significance was tested for differences in pair or multiple comparisons of different bacterial strains.
Reverse transcriptase-polymerase chain reaction
Quantitative real-time RT-PCR
- Xoo :
Xanthomonas oryzae pv. oryzae
Gottesman S, Storz G. Bacterial small RNA regulators: versatile roles and rapidly evolving variations. Cold Spring Harb Perspect Biol. 2011. https://doi.org/10.1101/cshperspect.a003798.
Storz G, Vogel J, Wassarman KM. Regulation by small RNAs in bacteria: expanding frontiers. Mol Cell. 2011;43:880–91.
Papenfort K, Vogel J. Regulatory RNA in bacterial pathogens. Cell Host Microbe. 2010;8:116–27.
Gottesman S. The small RNA regulators of Escherichia coli: roles and mechanisms. Annu Rev Microbiol. 2004;58:303–28.
Christiansen JK, Nielsen JS, Ebersbach T, Valentin-Hansen P, Søgaard-Andersen L, Kallipolitis BH. Identification of small Hfq-binding RNAs in Listeria monocytogenes. RNA. 2006;12:1383–96.
Bejerano-Sagie M, Xavier KB. The role of small RNAs in quorum sensing. Curr Opin Microbiol. 2007;10:189–98.
Toledo-Arana A, Dussurget O, Nikitas G, Sesto N, Guet-Revillet H, Balestrino D, Loh E, Gripenland J, Tiensuu T, Vaitkevicius K, Barthelemy M, Vergassola M, Nahori MA, Soubigou G, Régnault B, Coppée JY, Lecuit M, Johansson J, Cossart P. The Listeria transcriptional landscape from saprophytism to virulence. Nature. 2009;459:950–6.
Chabelskaya S, Gaillot O, Felden B. A Staphylococcus aureus small RNA is required for bacterial virulence and regulates the expression of an immune-evasion molecule. PLoS Pathog. 2010. https://doi.org/10.1371/journal.ppat.1000927.
Cui Y, Chatterjee A, Yang H, Chatterjee AK. Regulatory network controlling extracellular proteins in Erwinia carotovora subsp. carotovora: FlhDC, the master regulator of flagellar genes, activates rsmB regulatory RNA production by affecting gacA and hexA (lrhA) expression. J Bacteriol. 2008;190:4610–23.
Livny J, Waldor MK. Identification of small RNAs in diverse bacterial species. Curr Opin Microbiol. 2007;10:96–101.
Altuvia S. Identification of bacterial small non-coding RNAs: experimental approaches. Curr Opin Microbiol. 2007;10:257–61.
Gottesman S. Micros for microbes: non-coding regulatory RNAs in bacteria. Trends Genet. 2005;21:399–404.
Massé E, Escorcia FE, Gottesman S. Coupled degradation of a small regulatory RNA and its mRNA targets in Escherichia coli. Genes Dev. 2014;17:2374–83.
Zhang A, Wassarman KM, Rosenow C, Tjaden BC, Storz G, Gottesman S. Global analysis of small RNA and mRNA targets of Hfq. Mol Microbiol. 2003;50:1111–24.
Tronnet S, Garcie C, Brachmann AO, Piel J, Oswald E, Martin P. High iron supply inhibits the synthesis of the genotoxin colibactin by pathogenic Escherichia coli through a non-canonical Fur/RyhB-mediated pathway. Pathog Dis. 2017. https://doi.org/10.1093/femspd/ftx066.
Porcheron G, Habib R, Houle S, Caza M, Lépine F, Daigle F, Massé E. The small RNA RyhB contributes to siderophore production and virulence of uropathogenic Escherichia coli. Infect Immun. 2014;82:5056–68.
Oglesby AG, Murphy ER, Iyer VR, Payne SM. Fur regulates acid resistance in Shigella flexneri via RyhB and ydeP. Mol Microbiol. 2005;58:1354–67.
Murphy ER, Payne SM. RyhB, an iron-responsive small RNA molecule regulates Shigella dysenteriae virulence. Infect Immun. 2007;75:3470–7.
Sesto N, Koutero M, Cossart P. Bacterial and cellular RNAs at work during Listeria infection. Future Microbiol. 2014;9:1025–37.
Chao Y, Papenfort K, Reinhardt R, Sharma CM, Vogel J. An atlas of Hfq-bound transcripts reveals 3’ UTRs as a genomic reservoir of regulatory small RNAs. EMBO J. 2012;31:4005–19.
Wang L, Cai X, Wu S, Bomjan R, Nakayasu ES, Händler K, Hinton JCD, Zhou DG. InvS coordinates expression of PrgH and FimZ, and is required for invasion of epithelial Cells by Salmonella enterica serovar Typhimurium. J Bacteriol. 2017. https://doi.org/10.1128/JB.00824-16.
Chevalier C, Boisset S, Romilly C, Masquida B, Fechter P, Geissmann T, Vandenesch F, Romby P. Staphylococcus aureus RNAIII Binds to two distant regions of coa mRNA to arrest translation and promote mRNA degradation. PLoS Pathog. 2010. https://doi.org/10.1371/journal.ppat.1000809.
Liang X, Yu C, Sun J, Liu H, Landwehr C, Holmes D, Ji YD. Inactivation of a two-component signal transduction system, SaeRS, eliminates adherence and attenuates virulence of Staphylococcus aureus. Infect Immun. 2006;74:4655–65.
Bronesky D, Wu Z, Marzi S, Walter P, Geissmann T, Moreau K, Vandenesch F, Caldelari I, Romby P. Staphylococcus aureus RNAIII and its regulon link quorum sensing, stress responses, metabolic adaptation, and regulation of virulence gene expression. Annu Rev Microbiol. 2016;70:299–316.
Gupta RK, Luong TT, Lee CY. RNAIII of the Staphylococcus aureus agr system activates global regulator MgrA by stabilizing mRNA. Proc Natl Acad Sci U S A. 2015;112:14036–41.
Bonas U, Schulte R, Fenselau S, Minsavage GV, Staskawicz BJ, Stall RE. Isolation of a gene cluster from Xanthomonas campestris pv. vesicatoria that determines pathogenicity and the hypersensitive response on pepper and tomato. Mol Plant-Microbe Interact. 1991;4:81–8.
Findeiss S, Schmidtke C, Stadler PF, Bonas U. A novel family of plasmid-transferred anti-sense ncRNAs. RNA Biol. 2010;7:120–4.
Jiang R, Tang D, Chen X, He Y, Feng J, Jiang B, Lu G, Lin M, Tang J. Identification of four novel small non-coding RNAs from Xanthomonas campestris pathovar campestris. BMC Genomics. 2010. https://doi.org/10.1186/1471-2164-11-316.
Liang H, Zhao Y, Zhang J, Wang X, Fang R, Jia Y. Identification and functional characterization of small non-coding RNAs in Xanthomonas oryzae pathovar oryzae. BMC Genomics. 2011. https://doi.org/10.1186/1471-2164-12-87.
Schmidtke C, Findeiss S, Sharma CM, Kuhfuss J, Hoffmann S, Vogel J, Stadler PF, Bonas U. Genome-wide transcriptome analysis of the plant pathogen Xanthomonas identifies sRNAs with putative virulence functions. Nucleic Acids Res. 2012;40:2020–31.
An S, Febrer M, McCarthy Y, Tang D, Clissold L, Kaithakottil G, Swarbreck D, Tang J, Rogers J, Dow JM, Ryan RP. High-resolution transcriptional analysis of the regulatory influence of cell-to-cell signalling reveals novel genes that contribute to Xanthomonas phytopathogenesis. Mol Microbiol. 2013;88:1058–69.
Schmidtke C, Abendroth U, Brock J, Serrania J, Becker A, Bonas U. Small RNA sX13: a multifaceted regulator of virulence in the plant pathogen Xanthomonas. PLoS Pathog. 2013. https://doi.org/10.1371/journal.ppat.1003626.
Ji H, Dong H. Key steps in type III secretion system (T3SS) towards translocon assembly with potential sensorat plant plasma membrane. Mol Plant Pathol. 2015;16:762–73.
Sugio A, Yang B, White FF. Characterization of the hrpF pathogenicity peninsula of Xanthomonas oryzae pv. oryzae. Mol Plant-Microbe Interact. 2005;18:546–54.
Chen X, Tang D, Jiang R, He Y, Jiang B, Lu G, Tang J. sRNA-Xcc1, an integron-encoded transposon- and plasmid-transferred trans-acting sRNA, is under the positive control of the key virulence regulators HrpG and HrpX of Xanthomonas campestris pathovar campestris. RNA Biol. 2011;8:947–53.
Thomas NA, Brett Finlay B. Establishing order for type III secretion substrates-a hierarchical process. Trends Microbiol. 2003;11:398–403.
Yang B, Zhu W, Johnson LB, White FF. The virulence factor AvrXa7 of Xanthomonas oryzae pv. oryzae is a type III secretion pathway-dependent nuclear-localized double-stranded DNA-binding protein. Proc Natl Acad Sci U S A. 2000;97:9807–12.
Tsuge S, Furutani A, Fukunaka R, Oku T, Tsuno K, Ochiai H, Inoue Y, Kaku H, Kubo Y. Expression of Xanthomonas oryzae pv. oryzae hrp genes in XOM2, a novel synthetic medium. J Gen Plant Pathol. 2002;68:363–71.
Gómez-Lozano M, Marvig RL, Molin S, Long KS. Genome-wide identification of novel small RNAs in Pseudomonas aeruginosa. Environ Microbiol. 2012;14:2006–16.
Vogel J, Luisi BF. Hfq and its constellation of RNA. Nat Rev Microbiol. 2011;9:578–89.
Rashid MM, Ikawa Y, Tsuge S. GamR, the LysR-type galactose metabolism regulator, regulates hrp gene expression via transcriptional activation of two key hrp regulators, HrpG and HrpX, in Xanthomonas oryzae pv. oryzae. Appl Environ Microbiol. 2016;82:3947–58.
Bauer DW, Wei ZM, Beer SV, Collmer A. Erwinia chrysanthemi harpinEch: an elicitor of the hypersensitive response that contributes to soft-rot pathogenesis. Mol Plant-Microbe Interact. 1995;8:484–91.
He SY, Nomuraa K, Whittam TS. Type III protein secretion mechanism in mammalian and plant pathogens. Biochim Biophys Acta. 2004;1694:181–206.
Wengelnik K, Bonas U. HrpXv, an AraC-type regulator, activates expression of five of the six loci in the hrp cluster of Xanthomonas campestris pv. vesicatoria. J Bacteriol. 1996;178:3462–9.
Wengelnik K, Rossier O, Bonas U. Mutations in the regulatory gene hrpG of Xanthomonas campestris pv. vesicatoria result in constitutive expression of all hrp. J Bacteriol. 1999;181:6828–31.
Kvitko BH, Ramos AR, Morello JE, Oh HS, Collmer A. Identification of harpins in Pseudomonas syringae pv. tomato DC3000, which are functionally similar to HrpK1 in promoting translocation of type III secretion system effectors. J Bacteriol. 2007;189:8059–72.
Zhang H, Tao Z, Hong H, Chen Z, Wu C, Li X, Xiao J, Wang S. Transposon-derived small RNA is responsible for modified function of WRKY45 locus. Nat Plants. 2016. https://doi.org/10.1038/nplants.2016.16.
Wengelnik K, Van den Ackerveken G, Bonas U. HrpG, a key hrp regulatory protein of Xanthomonas campestris pv. vesicatoria is homologous to two component response regulators. Mol Plant-Microbe Interact. 1996;9:704–12.
Guo Y, Figueiredo F, Jones J, Wang N. HrpG and HrpX play global roles in coordinating different virulence traits of Xanthomonas axonopodis pv. citri. Mol Plant-Microbe Interact. 2010;24:649–61.
Liu Y, Long J, Shen D, Song C. Xanthomonas oryzae pv. oryzae requires H-NS-family protein XrvC to regulate virulence during rice infection. FEMS Microbiol Lett. 2016. https://doi.org/10.1093/femsle/fnw067.
Ikawa Y, Tsuge S. The quantitative regulation of the hrp regulator HrpX is involved in sugar-source-dependent hrp gene expression in Xanthomonas oryzae pv. oryzae. FEMS Microbiol Lett. 2016. http://doi.org/10.1093/femsle/fnw071.
Okan NA, Bliska JB, Karzai AW. A role for the SmpB-SsrA system in Yersinia pseudotuberculosis pathogenesis. PLoS Pathog. 2006;2:50–62.
Fu ZQ, Guo M, Alfano JR. Pseudomonas syringae HrpJ is a type III secreted protein that is required for plant pathogenesis, injection of effectors, and secretion of the HrpZ1 Harpin. J Bacteriol. 2006;188:6060–9.
Wei ZM, Laby RJ, Zumoff CH, Bauer DW, He SY, Collmer A, Beer SV. Harpin, elicitor of the hypersensitive response produced by the plant pathogen Erwinia amylovora. Science. 2004;257:85–8.
Bronsard J, Pascreau G, Sassi M, Mauro T, Augagneur Y, Felden B. sRNA and cis-antisense sRNA identification in Staphylococcus aureus highlights an unusualsRNA gene cluster with one encoding a secreted peptide. Sci Rep. 2017. https://doi.org/10.1038/s41598-017-04786-3.
Park W, Rengaraj D, Kil DY, Kim H, Lee HK, Song KD. RNA-Seq analysis of the kidneys of broiler chickens fed diets containing different concentrations of calcium. Sci Rep. 2017. https://doi.org/10.1038/s41598-017-11379-7.
Li YR, Che YZ, Zou HS, Cui YP, Guo W, Zou LF, Biddle EM, Yang CH, Chen GY. Hpa2 required by HrpF to translocate Xanthomonas oryzae transcriptional activator-like effectors into rice for pathogenicity. Appl Environ Microbiol. 2011;77:3809–18.
Shi LW. SPSS19.0 statistical analysis from accidence to conversance. Beijing: Tsinghua University Press; 2012. p. 109–43.
We thank our former colleague Dr. Dan Shen for initiating this study.
This study was supported by Natural Science Foundation of China (grant number 31772247), the national key research and development plan (grant number 2017YFD0200901) and the Jiangsu Provincial Scientific Plan Program (grant number BK20150668).
Availability of data and materials
All data supporting the findings in this manuscript is included here or in the Supporting Information. The raw reads files have been deposited in the National Center for Biotechnology Information’s Sequence Read Archive (SRA), and are accessible via the SRA accession SRR7867914 (https://www.ncbi.nlm.nih.gov/sra/?term=SRR7867914).
Ethics approval and consent to participate
Not applicable. Both of the plant materials and bacterial strains were obtained from our lab and the materials were cultivated by the authors. Our study did not involve endangered or protected species. Therefore, there is no need to get a permit for collection of plant material from natural habitats.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Information on 12 sRNA selected from the Xanthomonas oryzae pv. oryzae sRNA-Seq profiling. (DOCX 14 kb)
Figure S1. Bacterial blight severities of the susceptible rice variety Nipponbare inoculated with sRNA-related strains. (DOCX 80 kb)
Table S2. Strains and plasmids used and created in this study. (DOCX 15 kb)
Table S3. Information on genes tested and primers used in this study. (DOCX 17 kb)
About this article
Cite this article
Hu, Y., Zhang, L., Wang, X. et al. Two virulent sRNAs identified by genomic sequencing target the type III secretion system in rice bacterial blight pathogen. BMC Plant Biol 18, 237 (2018). https://doi.org/10.1186/s12870-018-1470-7