Research article | Open | Published:
Tigecycline challenge triggers sRNA production in Salmonella enterica serovar Typhimurium
BMC Microbiologyvolume 12, Article number: 195 (2012)
Bacteria employ complex transcriptional networks involving multiple genes in response to stress, which is not limited to gene and protein networks but now includes small RNAs (sRNAs). These regulatory RNA molecules are increasingly shown to be able to initiate regulatory cascades and modulate the expression of multiple genes that are involved in or required for survival under environmental challenge. Despite mounting evidence for the importance of sRNAs in stress response, their role upon antibiotic exposure remains unknown. In this study, we sought to determine firstly, whether differential expression of sRNAs occurs upon antibiotic exposure and secondly, whether these sRNAs could be attributed to microbial tolerance to antibiotics.
A small scale sRNA cloning strategy of Salmonella enterica serovar Typhimurium SL1344 challenged with half the minimal inhibitory concentration of tigecycline identified four sRNAs (sYJ5, sYJ20, sYJ75 and sYJ118) which were reproducibly upregulated in the presence of either tigecycline or tetracycline. The coding sequences of the four sRNAs were found to be conserved across a number of species. Genome analysis found that sYJ5 and sYJ118 mapped between the 16S and 23S rRNA encoding genes. sYJ20 (also known as SroA) is encoded upstream of the tbpAyabKyabJ operon and is classed as a riboswitch, whilst its role in antibiotic stress-response appears independent of its riboswitch function. sYJ75 is encoded between genes that are involved in enterobactin transport and metabolism. Additionally we find that the genetic deletion of sYJ20 rendered a reduced viability phenotype in the presence of tigecycline, which was recovered when complemented. The upregulation of some of these sRNAs were also observed when S. Typhimurium was challenged by ampicillin (sYJ5, 75 and 118); or when Klebsiella pneumoniae was challenged by tigecycline (sYJ20 and 118).
Small RNAs are overexpressed as a result of antibiotic exposure in S. Typhimurium where the same molecules are upregulated in a related species or after exposure to different antibiotics. sYJ20, a riboswitch, appears to possess a trans-regulatory sRNA role in antibiotic tolerance. These findings imply that the sRNA mediated response is a component of the bacterial response to antibiotic challenge.
Multiple studies demonstrate that non coding RNAs (or small RNAs (sRNAs)) possess regulatory roles in the bacterial stress response[1–4]. Bacterial sRNA regulators typically range from 50 – 250 nts and are often transcribed from intergenic regions (IGRs), although open reading frames may also encode sRNAs. Most sRNAs act as regulators at the post-transcriptional level by base-pairing with target mRNAs; these sRNA-mRNA binding regions are often short and imperfect and may require an additional RNA chaperone, which in most cases is the Hfq protein[6, 7]. This imperfect binding allows each sRNA molecule to control multiple targets, where either the translation of the target mRNA is upregulated, or more commonly inhibited.
Many sRNA regulators are upregulated when bacteria sense environmental stress: these include oxidative stress, low pH environment, nutrient deprivation and glucose-phosphate stress. Despite overwhelming evidence that sRNAs play a role when bacteria experience physiological stress, no systematic study has been undertaken to ascertain the impact or levels of sRNA production in bacteria when antibiotics are present.
Naturally susceptible pathogens can develop drug resistance when treated with antibiotics. Genetically acquired antibiotic resistance in pathogenic bacteria, via spontaneous / random mutations and horizontal gene transfer, is a significant issue in the treatment of infectious diseases. Intrinsic regulatory networks such as those mediated by the transcriptional regulators MarA, SoxS and RamA are also implicated in the development of antibiotic resistance particularly since these systems control the influx / efflux of antibiotics. Thus far studies that have focused on the intrinsic antibiotic resistome are limited to gene and protein networks mediated by these gene operons or other transcription factors[11–13]. Hence the role of the newly uncovered class of regulatory molecules such as sRNAs in controlling or contributing to the antimicrobial resistance phenotype is largely unknown. Some evidence for the role of sRNAs in mediating antimicrobial resistance already exists: for example, the expression of bacterial outer membrane proteins, OmpF and OmpC, involved in antibiotic import, is controlled by the sRNAs MicF and MicC respectively[14–16]. Additionally, the overexpression of another sRNA (DsrA) was recently found to induce multidrug resistance in Escherichia coli via the MdtEF efflux pump. Nevertheless, whether the functional role of MicF, MicC and DsrA is indeed part of the bacteria’s intrinsic stress response to antibiotic challenge remains unknown.
Tigecycline is a member of the glycylcycline group of antibiotics, and was registered in the EU in April 2006. This bacteriostatic antibiotic acts as a protein synthesis inhibitor by binding to the 30S ribosomal subunit. Tigecycline is active against a broad range of bacteria, with only few naturally resistant exceptions, namely, Proteus spp., Morganella morganii, Providencia spp., and Pseudomonas aeruginosa. Specifically, tigecycline is effective against multidrug resistant bacteria such as Staphylococcus aureus (MRSA), vancomycin-resistant Enterococcus (VRE), extended-spectrum beta-lactamase (ESBL)-expressing Enterobacteriaceae, and carbapenem-resistant strains[20–22]. Reports of resistance to tigecycline have been rare in naturally susceptible pathogens, however in resistant variants efflux pump overexpression has contributed to tigecycline resistance[23–28].
Salmonella, a member of Enterobacteriaceae, encodes both the ramA transcriptional factor and the acrAB efflux pump, which when overexpressed confers tigecycline resistance. Additionally, Salmonella represents a model bacterium for sRNA mining and genome manipulation, making it an ideal system for our study, but more importantly represents a paradigm for other members of Enterobacteriaceae. Hence in this study we used a cloning strategy to determine the sRNA population after tigecycline exposure in Salmonella enterica serovar Typhimurium, and also whether the absence of these sRNAs would render the cells less adaptable to tigecycline challenge.
cDNA library construction and analysis
A cDNA library was constructed from the cells that were challenged by half the minimal inhibitory concentration (MIC) of tigecycline (0.125 μg/ml) at OD600 = 0.6. Approximately ~6000 clones were obtained; from these 200 random candidates were sequenced and analysed. The nature of the cDNA library construction procedure (see Materials and Methods) allowed us to obtain the sequences in an orientation specific manner. The cDNA sequences were mapped to the S. Typhimurium SL1344 genome (FQ312003) using BLAST (http://blast.ncbi.nlm.nih.gov/Blast.cgi). Of the mapped sequences, 31% encoded tRNAs; 6% and 9% matched to rRNAs and protein coding sequences, respectively; 4% partially overlapped with open reading frames (ORFs), and 50% aligned to IGRs. Of all the IGR readings, 90% were located between the 16S and 23S rRNA encoding genes (Figure1).
Analyses of the cDNA sequences encoding partial ORFs indicated which genes were expressed in the presence of tigecycline. As stated above, 9% of the sequences identified matched to rRNAs, in addition to a further sequence which was found to overlap the 30S ribosomal protein and another mapped to elongation factor tu. This is perhaps not surprising, given that the specific target for tigecycline is the ribosome. On the other hand, sequences overlapping known stress-response genes were also captured in the cDNA library, e.g. dinF and a gene encoding a putative outer membrane protein (SL1344_1151). The dinF gene is a member of the SOS response family and encodes an efflux pump which belongs to the multidrug and toxic compound extrusion (MATE) family, and SL1344_1151, encoding a putative outer membrane protein homologous to ycfR in E. coli, which influences biofilm formation through stress response and surface hydrophobicity. The expression of these genes supports our hypothesis that challenge at half the MIC of tigecycline triggers a stress response. Of note, the cDNA library also contained sequences of different lengths that mapped to open reading frames, which we postulate to be a result of mRNA degradation, rather than a representation of bona fide sRNA regulators. Meanwhile, 4% of all sequences that partially overlap ORFs, all do so at the 5’ end of the ORFs. This suggests that these sequences might be 5’ untranslated regions, or encode riboswitches and/or control the expression of the downstream genes.
Northern blot verification
Northern blot analysis was performed on RNA extracted from SL1344 that were either unchallenged or challenged with half the MIC of tigecycline. Since most sRNAs are produced from IGRs, only sequences from these regions (100 out of 200 in total) were selected for further validation by northern blot analysis. As 90% of the IGR sequences are located between 16S and 23S rRNA coding sequences, most of which are identical, there were 20 unique IGR sequences (including those located between 16S and 23S rRNA) that were assayed, of which four (encoding sYJ5, sYJ20, sYJ75 and sYJ118) were found to consistently show elevated expression with tigecycline challenge (Figure2A). The remaining sRNA candidates were either not detectable by northern blots, or did not show differential levels of transcription. Correspondingly all further analyses focused on these four sRNAs. The relative fold increase in sRNA expression was determined by northern blots in challenged versus unchallenged cells. Upon tigecycline exposure, the expression levels of sYJ5, sYJ20, sYJ75 and sYJ118 (performed in triplicate using densitometric analyses), relative to unchallenged cells, were increased to 8, 2, 2 and 8 fold, respectively (Figure2A and B). We also tested the level of the four sRNAs in cells challenged with half the MIC of tetracycline (1 μg/ml). As expected, all of the four sRNAs were also found to be upregulated compared to the control sample (Figure3A). This is possibly due to the fact that tigecycline and tetracycline are related compounds, and they may as well trigger stress response pathways that share a common set of regulatory molecules. Of note and as shown in Figure4A, the level of 5S RNA was not affected by the presence of half the MIC of tigecycline or tetracycline (5Stigecycline: 5Scontrol = 0.88, 5Stetracycline : 5Scontrol = 1.15, average of 4 different experiments).
All four sRNA sequences were searched against S. Typhimurium SL1344 using NCBI BLAST. The sYJ5 encoding sequence is located between the 16S (SL1344_rRNA0001) and 23S rRNA (SL1344_rRNA0002) coding loci on the sense strand (Figure2C (i)). BLAST analysis uncovered two additional identical copies in the genome sequence of SL1344 (one between SL1344_rRNA0014 and SL1344_rRNA0015, the other SL1344_rRNA0017 and SL1344_rRNA0018).
Similar to sYJ5, sYJ118 is also encoded from the IGR between the 16S and 23S rRNA coding sequences, but from a different genetic locus (SL1344_rRNA0009 – SL1344_rRNA0010, Figure2C (iv)). The sequence encoding sYJ118 has an identical copy (SL1344_rRNA0011 – SL1344_rRNA0012) and additionally five other paralogs with 93% - 99% identity on the SL1344 chromosome.
The encoding sequence of sYJ75 is flanked by entC downstream (encoding isochorismate synthase), and fepB upstream (encoding the iron-enterobactin transporter periplasmic binding protein) (Figure2C (iii)). It also has a paralog that shares 90% identity, starting at position 1515629 on the S. Typhimurium SL1344 genome and located between pntB (encoding pyridine nucleotide transhydrogenase β subunit) and an un-annotated gene (encoding a putative membrane protein).
sYJ20 was previously identified by Vogel et al. in E. coli as SroA, encoded by a sequence downstream of yabN (encoding SgrR, a transcriptional regulator in E. coli) and upstream of tbpA (encoding the thiamine-binding periplasmic protein, homologous to thiB in E. coli) (Figures2C (ii) and5A).
sYJ5, sYJ20 (SroA) and sYJ118 are all highly conserved within the different members of Enterobacteriaceae, although the coding sequences of sYJ5, sYJ20 and sYJ118 are also found in other families of bacteria (such as sYJ5 and sYJ118 in Prevotella ruminicola, sYJ20 in Marinobacter aquaeolei VT8), in plants (such as sYJ20 and sYJ118 in Zea mays cultivar line T63) and in animals (sYJ118 in Gryllus bimaculatus). In contrast, sYJ75 is only found in Salmonella, Enterobacter, Photorhabdus and Citrobacter.
sYJ20 (SroA), sYJ5, sYJ75 and sYJ118 in other species and relevance to other drug classes
We proceeded to determine whether the increased expression of these sRNAs would be Salmonella specific or drug-class specific. Hence, we assessed the levels of our sRNA candidates (sYJ5, sYJ20 and sYJ118) in other members of Enterobacteriaceae (Klebsiella pneumoniae and Escherichia coli) when challenged with sub-inhibitory levels of tigecycline (sYJ75 was not included since it is not encoded in the tested species). Additionally, in order to determine whether these sRNAs are upregulated solely as a result of tigecycline challenge or are generally upregulated as a result of sub-inhibitory antibiotic challenge, S. Typhimurium SL1344 was challenged with either half the MIC of ampicillin (1 μg/ml) or ciprofloxacin (0.0156 μg/ml). As shown in Figure3B, none of the four tested sRNAs were upregulated in response to ciprofloxacin exposure, whilst three (sYJ5, sYJ75 and sYJ118) were found to be upregulated in the presence of ampicillin. Interestingly, E. coli did not upregulate the expression of the three candidate sRNAs (sYJ5, sYJ20 and sYJ118) in response to challenge at half the MIC of tigecycline. However, sYJ118 exhibited an elevated level of expression in K. pneumoniae in the presence of tigecycline (Figure3B). Of note, although the sYJ20 (SroA) coding sequence is present in K. pneumoniae, two transcripts were detected after hybridisation. However it was the larger RNA species that appeared upregulated in RNA derived from Klebsiella cells challenged with half the MIC of tigecycline. Hence we surmise that this larger RNA transcript, consistent with the larger intergenic region in K. pneumoniae, is where the sYJ20 homolog coding sequence is located. From these results we show that the upregulation of sRNAs identified in this study are neither species nor drug specific in the presence of unrelated classes of antibiotics.
5’ Rapid Amplifed cDNA Ends (5’ RACE) of sYJ20 (SroA)
To determine the transcriptional start site (TSS) of sYJ20 (shared with the one of tbpA), we performed 5’ RACE analysis. As shown in Figure5, the 5’ RACE result reveals that the TSS of sYJ20 and tbpA lies 129 bases upstream of the start codon of tbpA, consistent with previous findings.
Quantitative real time PCR (qPCR)
sYJ20 (SroA): the upregulation of sYJ20 in S. Typhimurium challenged by half the MIC of tigecycline or tetracycline was quantified with qPCR. As shown in Figure6, compared to the control, cells challenged by tigecycline or tetracycline produced ~3 fold more sYJ20. Interestingly, the transcription level of the downstream gene, tbpA, was hardly affected by the presence of the antibiotics. This suggests that sYJ20, but not the tbpA gene product, is upregulated as a result of tigecycline or tetracycline challenge.
dinF (encoding an efflux pump) and ycfR (encoding a putative outer membrane protein): as mentioned previously, the RNA transcripts of these two stress responsive genes were picked up in the sRNA cloning and is suggestive that half the MIC of tigecycline does induce a stress response in S. Typhimurium. In order to confirm this, we performed a qPCR on S. Typhimurium challenged by half the MIC of tigecycline or tetracycline, and compared the transcriptional levels of dinF and ycfR to the control. As shown in Figure6, the transcriptional level of dinF increased to 7.0 and 2.8 fold when the cells were challenged by half the MIC of tigecycline and tetracycline, respectively; the level of ycfR increased to 390 and 210 fold when the cells were challenged by half the MIC of tigecycline and tetracycline, respectively.
Survival rate assays
Survival rate assays were performed to investigate whether the deletion of sYJ20 (SroA) would highlight any phenotypic deficiencies when challenged with tigecycline. Our initial tests showed that the MICs of the mutant (YJ104) and the wild type strains (SL1344) were identical to tigecycline (MIC: 0.25 μg/ml in RDM). We then performed growth curves in RDM, where both SL1344 and YJ104 exhibited similar growth rates, as determined by OD600 readings, even in the presence of tigecycline (data not shown).
To determine whether sYJ20 confers an advantage to bacterial survival in the presence of tigecycline challenge, the survival frequencies were determined for the wild type SL1344 and YJ104 in the presence of 1 ×, 2 ×, 4 × and 8 × MIC of tigecycline. Both SL1344 and YJ104 failed to form any colonies on 2 ×, 4 × and 8 × MIC plates after overnight incubation at 37°C. The survival rates for SL1344 and YJ104 at 1 × the MIC were ~2.1 × 10-7 and 1.1 × 10-7 respectively (Figure7). Despite this modest decrease, statistical analysis on four biological replicate experiments supports that the reduced survival rate observed in YJ104 is indeed significant (P < 0.05). The survival rate was restored upon complementation where YJ107 (YJ104/pACYC177·sYJ20) yielded a survival frequency close but higher than SL1344 (2.1 × 10-7, Figure7), and as expected the plasmid control YJ110 (YJ104/pACYC177) had a similar survival rate to YJ104 (1.0 × 10-7, Figure7). This reduction in the survival rate of YJ110 compared to the one of YJ107 was also found to be statistically significant (P < 0.05). Overall, it suggests that the absence of sYJ20 could confer a subtle but reduced survival rate in the presence of tigecycline.
Small RNAs are regulatory molecules that enhance a bacterium’s adaptability in a constantly changing environment[1–4]. As regulatory molecules, sRNAs have several advantages over their protein counterparts. Firstly, sRNAs consist of a short nucleotide sequence which does not require translation into a peptide sequence. This ensures that the response from sRNA mediated regulators would be much more rapid than protein mediated factors. Accordingly, modelling studies suggest that due to the rapid kinetics associated with sRNA production, the downstream regulon response is correspondingly prompt when compared to protein based factors, a valuable trait in constantly evolving environments. Moreover, base pairing flexibility presumably allows rapid evolution of sRNAs. Finally, sRNA-mRNA interaction generally lacks specificity and often imperfect binding occurs ensuring that more than one target mRNA is affected, thereby expanding the repertoire of the sRNA regulators. At antibiotic concentrations at or higher than the MIC, cells are likely to halt cellular replication and / or stop growing, or result in the accumulation of advantageous genomic mutations, which may permanently alter the transcriptional profiles of bacteria. Hence we surmised that the sRNAs upregulated in the cells under these conditions may not be a direct result of antibiotic stress response but possibly due to genetic mutations or global perturbations. Therefore, a cDNA library was constructed from the cells that were challenged by half the MIC of tigecycline at mid-log phase.
In support of our hypothesis, our screen identified genes involved in the stress response when the bacterial cells were challenged with half the MIC of tigecycline. These include a SOS response gene, dinF, encoding a MATE family efflux pump, and a gene homologous to ycfR in E. coli, encoding a putative outer membrane protein. QPCR confirms the upregulation of the two genes when S. Typhimurium is challenged with half the MIC of tigecycline or tetracycline (Figure6). Our finding of four sRNAs (sYJ20 (SroA), sYJ5, sYJ75 and sYJ118) that are upregulated in the presence of tigecycline or tetracycline provides the first direct evidence that sRNAs are differentially expressed upon antibiotic exposure. It is known that tetracycline triggers mRNA accumulation in bacteria. However, this is unlikely to be the case as increased transcription was not noted for e.g. tbpA (open reading frame lying downstream of sYJ20, Figure6), and the gene encoding the 5S RNA (Figure4A).
Two of the four sRNAs (sYJ5 and sYJ75) we describe in this study are novel. Additionally, our work shows that these four sRNAs are not species specific as both sYJ20 and sYJ118 are upregulated in K. pneumoniae when challenged with half the MIC of tigecycline, or drug specific as sYJ5, sYJ75 and sYJ118 are upregulated as a result of ampicillin challenge (Figure3B). Both sYJ118, previously identified as StyR-44 in Salmonella, and sYJ5, a novel sRNA discovered in this study, are located between 16S and 23S rRNA coding sequences (Figure2C).
Both tigecycline and tetracycline target the 30S ribosomal subunit in bacterial cells. This might trigger over-production of the 16S-23S rRNA molecules, which also includes sYJ5 and sYJ118. This may raise the possibility that sYJ5 and sYJ118 are “by-products” rather than bona fide sRNA regulators. However, in support of sYJ5 and sYJ118 being classed as sRNAs, not all 16S-23S rRNA intergenic regions identified in our screen were upregulated in the presence of tigecycline when assessed by northern blots (data not shown). Furthermore, only sYJ118, not sYJ5, was upregulated in K. pneumoniae when challenged with tigecycline (Figure3B). Taken together these data lead to the observation that these different inter 16S-23S rRNA regions (including the regions encoding either sYJ5 or sYJ118) may have alternative functions independent of rRNA processing, which could be regulatory sRNA.
In this work we have used the 5S RNA as a loading control for northern blot assays. Given that it is a ribosomal RNA we wondered whether the 5S RNA levels would be affected by either tigecycline or tetracycline exposure. As shown in Figure4A, the 5S RNA expression levels were unaltered when the cells were challenged with half the MIC of tigecycline or tetracycline, and therefore it is a suitable loading control for the northern blot assays.
The four sRNAs (sYJ5, sYJ20, sYJ75 and sYJ118) that were upregulated as a response to tigecycline challenge in S. Typhimurium were also upregulated in tetracycline challenged cells (Figures2A and3A). This is not surprising since both tigecycline and tetracycline target the 30S ribosomal subunit. It is possible that the similar mechanisms of action of tetracycline and tigecycline trigger comparable stress-responsive pathways, which possibly include sYJ5, sYJ20, sYJ75 and sYJ118.
sYJ75 has not been previously described and thus is also a novel sRNA discovered in this study. Its conservation among several species and its upregulation in S. Typhimurium upon challenge with tigecycline and tetracycline, (Figures2A,3A) suggest that sYJ75, combined with its conservation across different species, may represent a common denominator in the response to tigecycline / tetracycline exposure. Interestingly, none of the four sRNAs were found upregulated when S. Typhimurium was exposed to ciprofloxacin, or when E. coli was challenged with tigecycline (Figure3B).
When challenged with tigecycline, both S. Typhimurium and K. pneumoniae upregulated two sRNAs, namely sYJ20 and sYJ118 (Figure3B). Despite encoding these sequences, no upregulation was noted in E. coli cells exposed to tigecycline compared to the unexposed controls (Figure3B). This suggests two possibilities: the first, where the tigecycline stress response involving sRNAs in E. coli is different from that in K. pneumoniae and S. Typhimurium, and the second, where the sRNAs (sYJ20 and sYJ118) may be linked to regulatory networks contributing to tigecycline resistance, i.e. RamA, only found in S. Typhimurium and K.pneumoniae but not in E. coli[40, 41]. However TargetRNA predictions for sYJ20 for cognate mRNA binding partners, using default parameters, yields four mRNA sequences (Table1). Of note, pspB and pspA which are involved in stress-response and the virulence attributes of several bacterial species are potential targets of sYJ20. sYJ20-mediated control of the psp operon may explain the reduced fitness of the sroA (sYJ20) deleted Salmonella strain in a mouse infection model.
A recent study undertaken to map sRNA profiles in SL1344 using massive parallel sequencing technology identified 140 sRNAs. Notably, sYJ5 and sYJ75 were not identified in this large scale study which suggests that firstly, these sRNAs are produced as a result of conditional exposure e.g. tigecycline and secondly that our small scale screen is able to uncover novel sRNAs. The encoding sequences of three sRNAs (sYJ5, sYJ75 and sYJ118) identified in this screen have more than one paralog within S. Typhimurium’s genome, making it difficult to pinpoint their exact roles in the bacterial response against antibiotic challenge through genetic analysis. Due to this reason, only sYJ20 and its associated phenotype were investigated further.
sYJ20, also known as SroA, is encoded immediately upstream of the tbpAyabKyabJ operon (homologous to thiBPQ in E. coli) and contains a THI-box sequence required as a riboswitch for the modulation of the tbpAyabKyabJ operon (Figure5). The deletion of the chromosomal sequence of sYJ20 would have very likely removed the TSS of the downstream gene tbpA (Figure5). However, tbpA transcript levels remained unaltered upon tigecycline / tetracycline exposure (Figure6). Therefore the polar effect of the sYJ20 deletion is considered to be minimal.
When survival rate assays were performed a subtle but reproducible deficiency (P < 0.05) as reflected by a reduction in the viability in the ΔsYJ20 strain (YJ104) compared to the wild type strain (SL1344) (Figure7) was observed. This deficiency was alleviated when a plasmid encoding allele of sYJ20 was transformed in YJ104 (i.e. YJ107), where the vector only control (i.e. YJ110) did not (Figure7). This subtle change of phenotype is not entirely surprising, as it has been observed that sRNA deletions usually have little, if any, effect. In fact, sYJ20, or SroA, has been linked to other phenotypes such as reduced fitness by a ΔsroA S. Typhimurium strain (sroA encodes sYJ20) during competitive infection with the wild type strain in mice. However it is not evident from the work whether the reduction in competitiveness of the ΔsroA S. Typhimurium strain is due to altered tbpA expression.
Previous work suggests that sYJ20 (SroA) may function as a riboswitch for the tbpAyabKyabJ (thiBPQ) operon in E. coli and that this regulatory role does not require Hfq. In our studies, we can show that the wild type strain S. Typhimurium (SL1344) produces sYJ20 (transcript size around 100 nts) in the presence of sub-inhibitory concentration of ciprofloxacin (0.0078 μg/ml) whilst the Δhfq strain produced less (Figure4B). This suggests that sYJ20, apart from its putative riboswitch role, can act as a trans-regulatory sRNA, as Hfq is typically required for functionality and stability by trans-encoded sRNAs. This is further supported by the two facts that A) the mild defect due to the chromosomal deletion of sYJ20 in SL1344 can be complemented by the plasmid-coding allele (YJ107), which cannot be attributed to its role as riboswitch, since the RNA transcripts of sYJ20 (on plasmid) and tbpA (on chromosome) are on separate strands, and B) sYJ20 was upregulated in S. Typhimurium challenged with half the MIC of tigecycline or tetracycline, where the transcriptional level of tbpA remained the same (Figure6). The transcript size of sYJ20, as detected by northern blot analysis, is approximately 100 nts which is consistent with the size reported in E. coli (93 nts). As has been suggested previously, it is possible that sYJ20 is generated by transcription attenuation of tbpAyabKyabJ; and the released short sYJ20 (around 100 nts) functions as a sRNA by regulating alternative targets in trans in the cell.
Our work shows that sRNAs upregulated in response to tigecycline exposure can also be produced in a non drug or species specific manner. The deletion of the sRNA, sYJ20 (SroA) confers a subtle survival disadvantage in the presence of tigecycline, possibly due to its role as a trans-regulatory sRNA after tigecycline exposure. Our results although preliminary, suggest that sRNA levels can be altered upon antibiotic exposure and presumably provide an initial survival advantage under antibiotic challenge. However, ongoing analyses are required to dissect the regulatory impact(s) of sRNA upregulation and its contribution to antibiotic resistance in bacteria.
Bacteria were cultured in Rich Defined Medium (RDM: 1 × M9 salts, 0.4% glucose, 1 × Essential Amino Acids (Gibco), 1 × Nonessential Amino Acids (Sigma-Aldrich, UK), 2 mM MgCl2, 0.1 mM CaCl2) unless otherwise stated. Typically, a strain was grown on a Luria-Bertani (LB) plate from frozen stock prior to experimental manipulations. A 1 in 100 dilution of fresh overnight culture was made in RDM and incubated in a 37°C shaker until OD600 reached 0.6, at which point half the MIC of the selected antibiotic (For SL1344: tigecycline (MIC = 0.25 μg/ml), tetracycline (MIC = 2 μg/ml), ciprofloxacin (MIC = 0.0312 μg/ml), or ampicillin (MIC = 2 μg/ml), for K. pneumoniae: tigecycline (MIC = 0.25 μg/ml), for E. coli: tigecycline (MIC = 0.0625 μg/ml), for JVS-0255: ciprofloxacin (MIC = 0.0156 μg/ml)) was added to the medium. The same volume of sterile water was added to another sample as a control. All strains used in this study are shown in Table2.
Minimum inhibitory concentration (MIC) determination
The MICs of all relevant strains in RDM to tigecycline, (gift from Wyeth Pharmaceuticals, US), tetracycline (Sigma-Aldrich, UK), ciprofloxacin and ampicillin (Sigma-Aldrich, UK) were determined and interpreted according to the BSAC protocols. In order to check whether concentrations at half the MIC would induce stress response rather than kill the cells in liquid medium, half of the MIC of the antibiotic was added to liquid culture at OD600 = 0.6 (sterilised water was added to the control). After growth for an hour or overnight, an aliquot of the culture was taken and spread on plates, to determine colony forming unit per ml (cfu/ml). Additionally growth curves were also generated based on the OD600 readings. The stress response was confirmed by comparison of the antibiotic challenged cells to the control on both the growth curves and the cfu/ml.
Cells were grown to OD600 = 0.6 prior to the addition of the antibiotic. After 1 hour of exposure, cells were harvested by centrifugation. The cell pellet was then resuspended in TRIzol reagent (Invitrogen) and the total RNA was extracted according to Santhakumar et al.. The resulting pellet was washed and resuspended in an appropriate amount of DEPC (Sigma, UK) treated water.
cDNA library construction
The cDNA library was constructed (according to the manufacturer’s instruction) using the Exact START Small RNA Cloning kit from Epicentre (Cambio, UK). Briefly, total RNA was digested with DNase I to remove any contaminating DNA, and small RNAs were enriched with Epicentre enrichment solution by precipitating RNA molecules longer than 200 nts. The enriched RNAs were treated with phosphatase (Cambio, UK) to convert 5’ triphosphate group of RNA molecules to 5’ monophosphate, and a poly-A tail was added to the 3’ end (according to the manufacturer’s instruction). The 5’ end of RNA was ligated with Acceptor Oligo that carries a NotI restriction site. Reverse transcription was performed to yield first cDNA strand, using a primer with poly-T at its 3’ end to cover the poly-A tail of RNA samples, and an AscI restriction site. After RNase digestion, the sample was subject to a PCR with Small RNA PCR Primer 1 and 2. The product was digested by NotI and AscI (New England Biolabs) and was subsequently cloned into the cloning-ready pCDC1-K vector (Cambio, UK). Since the 5’ ligation adaptor differs from the 3’ ligation adaptor, the cloning of these putative small RNA molecules is directional. All oligonucleotides used in this study are listed in Table3.
Ten micrograms of total RNA was separated on 8% polyacrylamide gel containing 8 M urea, and electro-transferred onto Hybond-N nylon membrane (GE Healthcare). The membrane was hybridised and washed according to Vogel et al., and exposed to a phosphor-imager (Fuji). Relative levels of increase in expression were determined by Multi Gauge 2.2 (Fujifilm). The bands were first normalised to the 5S RNA levels prior to calculating the fold increase of challenged versus unchallenged cells. The oligonucleotide probes used in the northern blot experiments are listed in Table3, and were end-labelled with γ32P-ATP using T4-polynucleotide kinase and purified prior to blot hybridisation.
Chromosomal sYJ20 (SroA) inactivation
The chromosomal inactivation of sYJ20 (SroA) was performed according to the manipulation strategy outlined by Datsenko and Wanner. Briefly, primers (sYJ20_Cm_F and sYJ20_Cm_R, sequences listed in Table3) with ~40 bases with 5’ end homology to the flanking regions of the sYJ20 coding sequence were used to amplify the cat locus on pKD3 by PCR. The PCR product was transformed into S. Typhimurium SL1344 carrying the plasmid pKD46. The transformed cells were selected on LB plates supplemented with chloramphenicol. Colonies were picked after an overnight incubation and the replacement of the chromosomal sYJ20 coding sequence with the cat cassette was verified by PCR and sequencing.
Quantitative Real Time PCR (qPCR)
All the primers for qPCR were tested for amplification efficiencies prior to use. cDNA was made with SuperScript® VILOTM cDNA Synthesis Kit (Invitrogen), which was then subject to qPCR with Platinum® SYBR® Green qPCR SuperMix-UDG (Invitrogen). The qPCR was performed using the Mx3005P qPCR system (Agilent/Strategene). Analyses of the QPCR data were undertaken using the MxPro algorithms (Agilent, UK) where the normalisation of the amplification data was to the 5S RNA levels.
The sequence spanning 40 bases upstream and 6 bases downstream up to the sYJ20 sRNA encoding sequence was amplified with primers sYJ20-HF and sYJ20-BR and cloned into pACYC177. The recombinant plasmid carrying the sYJ20 encoding sequence was verified by sequencing before transformation into YJ104 (SL1344 ΔsYJ20) to yield YJ107. Empty pACYC177 was also transformed into YJ104 to yield YJ110, used as a negative control. The levels of sY20 expression were confirmed by northern blots.
In order to determine the TSS of sYJ20 and tbpA, we employed the 5’ RACE System for Rapid Amplification of cDNA Ends (version 2.0, Invitrogen). Briefly, the first strand cDNA was produced using SuperScriptTM II Reverse Transcriptase (Invitrogen) with the GSP1 primer specifically matching to the tbpA RNA transcript. Following purification with the S.N.A.P column (Invitrogen), the 5’ end of the first strand cDNA was tailed with multiple C (cytidines) with dCTP and TdT. A PCR was performed with the Abridged Anchor Primer (Invitrogen) that targets the dC-tailed 5’ cDNA end, and the GSP2 primer attaching to the RNA transcript upstream of the GSP1 matching region. A nested PCR was also performed to increase the specificity with the nested GSP3 primer and the AUAP primer (Invitrogen). The PCR product was ligated onto the pGEM-T EASY vector, and was sequenced with the T7 Forward primer or the SP6 Reverse primer.
Survival rate assay
To assess the fitness of strains challenged with tigecycline, a survival rate assay of the wild type (SL1344), the ΔsYJ20 mutant (YJ104), the plasmid complemented strain (YJ107), and the vector only control (YJ110) was performed. One hundred microlitres of cells from fresh overnight RDM cultures were spread evenly on RDM plates supplemented with tigecycline at the MIC, 2 × MIC, 4 × MIC or 8 × MIC. The same batch of cells was also spread on RDM plates with no antibiotics to establish the baseline levels.
Altuvia S, Weinstein-Fischer D, Zhang A, Postow L, Storz G: A small, stable RNA induced by oxidative stress: role as a pleiotropic regulator and antimutator. Cell. 1997, 90 (1): 43-53. 10.1016/S0092-8674(00)80312-8.
Jin Y, Watt RM, Danchin A, Huang JD: Small noncoding RNA GcvB is a novel regulator of acid resistance in Escherichia coli. BMC Genomics. 2009, 10: 165-10.1186/1471-2164-10-165.
Morita T, Aiba H: Small RNAs making a small protein. Proc Natl Acad Sci U S A. 2007, 104 (51): 20149-20150. 10.1073/pnas.0710634105.
Wassarman KM, Storz G: 6S RNA regulates E. coli RNA polymerase activity. Cell. 2000, 101 (6): 613-623. 10.1016/S0092-8674(00)80873-9.
Vogel J, Bartels V, Tang TH, Churakov G, Slagter-Jager JG, Huttenhofer A, Wagner EG: RNomics in Escherichia coli detects new sRNA species and indicates parallel transcriptional output in bacteria. Nucleic Acids Res. 2003, 31 (22): 6435-6443. 10.1093/nar/gkg867.
Brennan RG, Link TM: Hfq structure, function and ligand binding. Curr Opin Microbiol. 2007, 10 (2): 125-133. 10.1016/j.mib.2007.03.015.
Sittka A, Lucchini S, Papenfort K, Sharma CM, Rolle K, Binnewies TT, Hinton JC, Vogel J: Deep sequencing analysis of small noncoding RNA and mRNA targets of the global post-transcriptional regulator, Hfq. PLoS Genet. 2008, 4 (8): e1000163-10.1371/journal.pgen.1000163.
Gottesman S: Micros for microbes: non-coding regulatory RNAs in bacteria. Trends Genet. 2005, 21 (7): 399-404. 10.1016/j.tig.2005.05.008.
Thi TD, Lopez E, Rodriguez-Rojas A, Rodriguez-Beltran J, Couce A, Guelfo JR, Castaneda-Garcia A, Blazquez J: Effect of recA inactivation on mutagenesis of Escherichia coli exposed to sublethal concentrations of antimicrobials. J Antimicrob Chemother. 2011, 66 (3): 531-538. 10.1093/jac/dkq496.
Wilke MH: Multiresistant bacteria and current therapy - the economical side of the story. Eur J Med Res. 2010, 15 (12): 571-576. 10.1186/2047-783X-15-12-571.
O'Regan E, Quinn T, Pages JM, McCusker M, Piddock L, Fanning S: Multiple regulatory pathways associated with high-level ciprofloxacin and multidrug resistance in Salmonella enterica serovar enteritidis: involvement of RamA and other global regulators. Antimicrob Agents Chemother. 2009, 53 (3): 1080-1087. 10.1128/AAC.01005-08.
Bush K: Alarming β-lactamase-mediated resistance in multidrug-resistant Enterobacteriaceae. Curr Opin Microbiol. 2010, 13 (5): 558-564. 10.1016/j.mib.2010.09.006.
Falagas ME, Rafailidis PI, Matthaiou DK: Resistance to polymyxins: Mechanisms, frequency and treatment options. Drug Resist Updat. 2010, 13 (4–5): 132-138.
Vogel J, Papenfort K: Small non-coding RNAs and the bacterial outer membrane. Curr Opin Microbiol. 2006, 9 (6): 605-611. 10.1016/j.mib.2006.10.006.
Delcour AH: Outer membrane permeability and antibiotic resistance. Biochim Biophys Acta. 2009, 1794 (5): 808-816. 10.1016/j.bbapap.2008.11.005.
Delihas N, Forst S: MicF: an antisense RNA gene involved in response of Escherichia coli to global stress factors. J Mol Biol. 2001, 313 (1): 1-12. 10.1006/jmbi.2001.5029.
Nishino K, Yamasaki S, Hayashi-Nishino M, Yamaguchi A: Effect of overexpression of small non-coding DsrA RNA on multidrug efflux in Escherichia coli. J Antimicrob Chemother. 2010, 66 (2): 291-296.
Hope R, Mushtaq S, James D, Pllana T, Warner M, Livermore DM: Tigecycline activity: low resistance rates but problematic disc breakpoints revealed by a multicentre sentinel survey in the UK. J Antimicrob Chemother. 2010, 65 (12): 2602-2609. 10.1093/jac/dkq370.
Doan TL, Fung HB, Mehta D, Riska PF: Tigecycline: a glycylcycline antimicrobial agent. Clin Ther. 2006, 28 (8): 1079-1106. 10.1016/j.clinthera.2006.08.011.
Kelesidis T, Karageorgopoulos DE, Kelesidis I, Falagas ME: Tigecycline for the treatment of multidrug-resistant Enterobacteriaceae: a systematic review of the evidence from microbiological and clinical studies. J Antimicrob Chemother. 2008, 62 (5): 895-904. 10.1093/jac/dkn311.
Peterson LR: A review of tigecycline–the first glycylcycline. Int J Antimicrob Agents. 2008, 32 (Suppl 4): S215-222.
Stein GE, Craig WA: Tigecycline: a critical analysis. Clin Infect Dis. 2006, 43 (4): 518-524. 10.1086/505494.
Visalli MA, Murphy E, Projan SJ, Bradford PA: AcrAB multidrug efflux pump is associated with reduced levels of susceptibility to tigecycline (GAR-936) in Proteus mirabilis. Antimicrob Agents Chemother. 2003, 47 (2): 665-669. 10.1128/AAC.47.2.665-669.2003.
Ruzin A, Keeney D, Bradford PA: AcrAB efflux pump plays a role in decreased susceptibility to tigecycline in Morganella morganii. Antimicrob Agents Chemother. 2005, 49 (2): 791-793. 10.1128/AAC.49.2.791-793.2005.
Ruzin A, Visalli MA, Keeney D, Bradford PA: Influence of transcriptional activator RamA on expression of multidrug efflux pump AcrAB and tigecycline susceptibility in Klebsiella pneumoniae. Antimicrob Agents Chemother. 2005, 49 (3): 1017-1022. 10.1128/AAC.49.3.1017-1022.2005.
Keeney D, Ruzin A, Bradford PA: RamA, a transcriptional regulator, and AcrAB, an RND-type efflux pump, are associated with decreased susceptibility to tigecycline in Enterobacter cloacae. Microbial drug resistance (Larchmont, NY. 2007, 13 (1): 1-6. 10.1089/mdr.2006.9990.
Keeney D, Ruzin A, McAleese F, Murphy E, Bradford PA: MarA-mediated overexpression of the AcrAB efflux pump results in decreased susceptibility to tigecycline in Escherichia coli. J Antimicrob Chemother. 2008, 61 (1): 46-53.
Hentschke M, Christner M, Sobottka I, Aepfelbacher M, Rohde H: Combined ramR mutation and presence of a Tn1721-associated tet(A) variant in a clinical isolate of Salmonella enterica serovar Hadar resistant to tigecycline. Antimicrob Agents Chemother. 2010, 54 (3): 1319-1322. 10.1128/AAC.00993-09.
Horiyama T, Nikaido E, Yamaguchi A, Nishino K: Roles of Salmonella multidrug efflux pumps in tigecycline resistance. J Antimicrob Chemother. 2010, 66 (1): 105-110.
Vogel J: A rough guide to the non-coding RNA world of Salmonella. Mol Microbiol. 2009, 71 (1): 1-11. 10.1111/j.1365-2958.2008.06505.x.
Brown DG, Swanson JK, Allen C: Two host-induced Ralstonia solanacearum genes, acrA and dinF, encode multidrug efflux pumps and contribute to bacterial wilt virulence. Appl Environ Microbiol. 2007, 73 (9): 2777-2786. 10.1128/AEM.00984-06.
Zhang XS, Garcia-Contreras R, Wood TK: YcfR (BhsA) influences Escherichia coli biofilm formation through stress response and surface hydrophobicity. J Bacteriol. 2007, 189 (8): 3051-3062. 10.1128/JB.01832-06.
Vanderpool CK, Gottesman S: The novel transcription factor SgrR coordinates the response to glucose-phosphate stress. J Bacteriol. 2007, 189 (6): 2238-2248. 10.1128/JB.01689-06.
Kroger C, Dillon SC, Cameron AD, Papenfort K, Sivasankaran SK, Hokamp K, Chao Y, Sittka A, Hebrard M, Handler K: The transcriptional landscape and small RNAs of Salmonella enterica serovar Typhimurium. Proc Natl Acad Sci U S A. 2012, 109 ((20): E1277-E1286.
Waters LS, Storz G: Regulatory RNAs in bacteria. Cell. 2009, 136 (4): 615-628. 10.1016/j.cell.2009.01.043.
Biedenbach DJ, Rhomberg PR, Mendes RE, Jones RN: Spectrum of activity, mutation rates, synergistic interactions, and the effects of pH and serum proteins for fusidic acid (CEM-102). Diagn Microbiol Infect Dis. 2010, 66 (3): 301-307. 10.1016/j.diagmicrobio.2009.10.014.
Kehrenberg C, Cloeckaert A, Klein G, Schwarz S: Decreased fluoroquinolone susceptibility in mutants of Salmonella serovars other than Typhimurium: detection of novel mutations involved in modulated expression of ramA and soxS. J Antimicrob Chemother. 2009, 64 (6): 1175-1180. 10.1093/jac/dkp347.
Chamot D, Owttrim GW: Regulation of cold shock-induced RNA helicase gene expression in the Cyanobacterium anabaena sp. strain PCC 7120. J Bacteriol. 2000, 182 (5): 1251-1256. 10.1128/JB.182.5.1251-1256.2000.
Chinni SV, Raabe CA, Zakaria R, Randau G, Hoe CH, Zemann A, Brosius J, Tang TH, Rozhdestvensky TS: Experimental identification and characterization of 97 novel npcRNA candidates in Salmonella enterica serovar Typhi. Nucleic Acids Res. 2010, 38 (17): 5893-5908. 10.1093/nar/gkq281.
Rosenblum R, Khan E, Gonzalez G, Hasan R, Schneiders T: Genetic regulation of the ramA locus and its expression in clinical isolates of Klebsiella pneumoniae. Int J Antimicrob Agents. 2011, 38 (1): 39-45. 10.1016/j.ijantimicag.2011.02.012.
Horiyama T, Nikaido E, Yamaguchi A, Nishino K: Roles of Salmonella multidrug efflux pumps in tigecycline resistance. J Antimicrob Chemother. 2011, 66 (1): 105-110. 10.1093/jac/dkq421.
Tjaden B: TargetRNA: a tool for predicting targets of small RNA action in bacteria. Nucleic Acids Res. 2008, 36 (Web Server issue): W109-113.
Darwin AJ: The phage-shock-protein response. Mol Microbiol. 2005, 57 (3): 621-628. 10.1111/j.1365-2958.2005.04694.x.
Santiviago CA, Reynolds MM, Porwollik S, Choi SH, Long F, Andrews-Polymenis HL, McClelland M: Analysis of pools of targeted Salmonella deletion mutants identifies novel genes affecting fitness during competitive infection in mice. PLoS Pathog. 2009, 5 (7): e1000477-10.1371/journal.ppat.1000477.
Papenfort K, Pfeiffer V, Lucchini S, Sonawane A, Hinton JC, Vogel J: Systematic deletion of Salmonella small RNA genes identifies CyaR, a conserved CRP-dependent riboregulator of OmpX synthesis. Mol Microbiol. 2008, 68 (4): 890-906. 10.1111/j.1365-2958.2008.06189.x.
Ontiveros-Palacios N, Smith AM, Grundy FJ, Soberon M, Henkin TM, Miranda-Rios J: Molecular basis of gene regulation by the THI-box riboswitch. Mol Microbiol. 2008, 67 (4): 793-803.
Valentin-Hansen P, Eriksen M, Udesen C: The bacterial Sm-like protein Hfq: a key player in RNA transactions. Mol Microbiol. 2004, 51 (6): 1525-1533. 10.1111/j.1365-2958.2003.03935.x.
Hoiseth SK, Stocker BA: Aromatic-dependent Salmonella typhimurium are non-virulent and effective as live vaccines. Nature. 1981, 291 (5812): 238-239. 10.1038/291238a0.
Blattner FR, Plunkett G, Bloch CA, Perna NT, Burland V, Riley M, Collado-Vides J, Glasner JD, Rode CK, Mayhew GF: The complete genome sequence of Escherichia coli K-12. Science New York, NY. 1997, 277 (5331): 1453-1462. 10.1126/science.277.5331.1453.
George AM, Hall RM, Stokes HW: Multidrug resistance in Klebsiella pneumoniae: a novel gene, ramA, confers a multidrug resistance phenotype in Escherichia coli. Microbiology (Reading, England). 1995, 141 (Pt 8): 1909-1920.
Andrews JM: BSAC standardized disc susceptibility testing method (version 8). J Antimicrob Chemother. 2009, 64 (3): 454-489. 10.1093/jac/dkp244.
Santhakumar D, Forster T, Laqtom NN, Fragkoudis R, Dickinson P, Abreu-Goodger C, Manakov SA, Choudhury NR, Griffiths SJ, Vermeulen A: Combined agonist–antagonist genome-wide functional screening identifies broadly active antiviral microRNAs. Proc Natl Acad Sci U S A. 2010, 107 (31): 13830-13835. 10.1073/pnas.1008861107.
Viegas SC, Pfeiffer V, Sittka A, Silva IJ, Vogel J, Arraiano CM: Characterization of the role of ribonucleases in Salmonella small RNA decay. Nucleic Acids Res. 2007, 35 (22): 7651-7664. 10.1093/nar/gkm916.
Vogel J, Wagner EG, Gerhart H: Approaches to identify novel non-messenger RNAs in bacteria and to investigate their biological functions: RNA mining. Handbook of RNA biochemistry. Edited by: Hartmann RK. 2005, Weinheim: Wiley-VCH-Verl, 595-613.
Datsenko KA, Wanner BL: One-step inactivation of chromosomal genes in Escherichia coli K-12 using PCR products. Proc Natl Acad Sci U S A. 2000, 97 (12): 6640-6645. 10.1073/pnas.120163297.
We thank Drs. P. Zucchi and H. Nicoloff for critical comments on the manuscript. Salary (Jing Yu and Thamarai Schneiders) and consumable support for this work were provided by the Department for Employment and Learning (Northern Ireland) through its “Strengthening the all-island Research Base” initiative.
TS and JY designed all the experiments. JY carried out the experiments. TS and JY wrote the manuscript. Both authors read and approved the final manuscript.