- Research article
- Open Access
Identification of novel growth phase- and media-dependent small non-coding RNAs in Streptococcus pyogenes M49 using intergenic tiling arrays
© Patenge et al.; licensee BioMed Central Ltd. 2012
- Received: 4 July 2012
- Accepted: 10 October 2012
- Published: 13 October 2012
Small non-coding RNAs (sRNAs) have attracted attention as a new class of gene regulators in both eukaryotes and bacteria. Genome-wide screening methods have been successfully applied in Gram-negative bacteria to identify sRNA regulators. Many sRNAs are well characterized, including their target mRNAs and mode of action. In comparison, little is known about sRNAs in Gram-positive pathogens. In this study, we identified novel sRNAs in the exclusively human pathogen Streptococcus pyogenes M49 (Group A Streptococcus, GAS M49), employing a whole genome intergenic tiling array approach. GAS is an important pathogen that causes diseases ranging from mild superficial infections of the skin and mucous membranes of the naso-pharynx, to severe toxic and invasive diseases.
We identified 55 putative sRNAs in GAS M49 that were expressed during growth. Of these, 42 were novel. Some of the newly-identified sRNAs belonged to one of the common non-coding RNA families described in the Rfam database. Comparison of the results of our screen with the outcome of two recently published bioinformatics tools showed a low level of overlap between putative sRNA genes. Previously, 40 potential sRNAs have been reported to be expressed in a GAS M1T1 serotype, as detected by a whole genome intergenic tiling array approach. Our screen detected 12 putative sRNA genes that were expressed in both strains. Twenty sRNA candidates appeared to be regulated in a medium-dependent fashion, while eight sRNA genes were regulated throughout growth in chemically defined medium. Expression of candidate genes was verified by reverse transcriptase-qPCR. For a subset of sRNAs, the transcriptional start was determined by 5′ rapid amplification of cDNA ends-PCR (RACE-PCR) analysis.
In accord with the results of previous studies, we found little overlap between different screening methods, which underlines the fact that a comprehensive analysis of sRNAs expressed by a given organism requires the complementary use of different methods and the investigation of several environmental conditions. Despite a high conservation of sRNA genes within streptococci, the expression of sRNAs appears to be strain specific.
- Streptococcus pyogenes
- Small noncoding RNAs
- Transcriptional regulation
In recent years, the role of small non-coding RNAs (sRNAs) in regulation of bacterial gene expression has become more evident; however, the large number of sRNAs identified in different bacterial species was unexpected[1–3]. Even though sRNAs were conventionally regarded as inhibitory antisense regulators, a significant number of sRNAs that activate bacterial gene expression have been characterized. Furthermore, regulatory mechanisms include both the stabilization and destabilization of target transcripts. Bacterial sRNAs influence the expression of genes involved in processes as diverse as stress response, sugar metabolism, and surface composition[6–10]. With sRNAs representing a whole new level of post-transcriptional regulation, it is no surprise that these molecules play an important role in the tightly controlled expression of virulence factors in many pathogens[11, 12].
We were interested in the regulatory sRNAome of Streptococcus pyogenes (group A streptococci, GAS), a common, exclusively human pathogen that causes a variety of diseases. GAS is responsible for mild superficial infections of the skin (impetigo contagiosa) and mucosal membranes (pharyngitis and tonsillitis). Additionally, there is a high global burden of severe GAS diseases such as post-streptococcal sequelae, and severe systemic (streptococcal toxic-shock-like syndrome) or invasive infections (necrotizing fasciitis), leading to over 500,000 deaths per year[13, 14]. The controlled expression of virulence factors plays a role in GAS infection, persistence in the host, and development of invasive diseases, which makes the investigation of virulence factors and their regulation a research priority. GAS expresses a large number of virulence factor genes coding for a variety of proteins, including surface components, lytic enzymes, proteinases, cytotoxins, superantigens, and immunoprotective proteins, that are controlled at least partially by the 30 stand-alone transcriptional regulators and 13 two-component systems identified to date in GAS. Virulence factor expression in GAS is highly responsive to environmental conditions and greatly depends on the growth phase.
Little is known, however, about the importance of sRNAs for virulence-related gene regulation in GAS. An overview of small RNAs in streptococci is nicely presented by Le Rhun and Charpentier. Several individual sRNAs have been identified in GAS[17, 18], with many more predicted by bioinformatic screens[12, 19]. Previous analysis of sRNA expression in a GAS M1T1 serotype using an intergenic tiling array approach identified 40 potential sRNAs, with a very low predicted overlap with candidate genes. The authors concluded that sRNA expression in GAS is serotype-dependent. The current work focused on sRNA expression in the skin isolate GAS M49. An intergenic tiling array identified 42 novel and 13 known sRNAs. Data from this experiment were compared to the results of the former GAS M1T1 study, and to predictions of two recently published bioinformatics tools[19, 21]. Additionally, we tested the regulation of sRNA expression in correlation to growth media and growth phase. We found very little overlap between the different screening methods, which underlines the importance of using several complementary methods, as well as several environmental conditions, to attain a comprehensive analysis of bacterial sRNAomes.
Identification of sRNAs in the intergenic regions of GAS M49
Custom intergenic tiling arrays representing the genome of S. pyogenes NZ131 (NCBI accession number: NC_011375) were designed to detect the expression of potential sRNAs. A total of 17,823 50-mer probes with an overlap of 15 bp were synthesized to cover the intergenic regions, with 9,082 probes representing the positive strand and 8,741 probes covering the negative strand. Additionally, 174 probes were designed as control probes covering tRNA genes or genes coding for known sRNAs. For example, we probed for fas X, SR914400, and SR1754950, all of which were detected in our experiments. Genedata Selector software was used to integrate genomes, tiling array probe sequences, sRNA predictions, and experimental data. Expression data were analysed in their genomic- and sequence-based contexts (Genedata AG).
Total RNA for the tiling array experiments was isolated from GAS M49 grown in chemically defined medium (CDM). Samples were taken from four biological replicates during both mid-log growth phase (OD600= 0.4–0.6) and stationary growth phase (OD600= 1.2). A signal intensity of >300 was set as a threshold. A positive signal required that a minimum of one probe specific for one strand showed an intensity above the threshold in at least three replicates. Intergenic regions featuring high intensities on both strands were manually removed following the analysis.
S. pyogenes M49 sRNAs, genomic location, Rfam predictions, and conservation across streptococcal genomes
Upstream adjacent gene
Downstream adjacent gene
Rfam and MOSES prediction
In a previous bioinformatic approach (MOSES), 20 probable candidates were predicted. In our array analysis, expression of five of the predicted sRNAs was confirmed under the conditions studied (Table1). Furthermore, we detected 12 streptococcal RNAs that were previously identified by Perez et al. (Table1). The phylogenetic conservation of the putative sRNAs was tested by BLAST analysis, and the taxonomic classification was presented following the nomenclature of Facklam. The tiling array technique with overlapping probes did not allow us to detect an accurate start site for the respective sRNA genes. Thus, we described the nucleotides represented by the active probes as the preliminary start and end (Table1).
Comparison of the sRNA expression data using different sRNA screens
We also compared our tiling array data with a previously published sRNA microarray study in GAS M1T1. The previous study identified 40 putative sRNAs. Even though the sequences of the candidates were conserved across streptococcal genomes, serotype-specific variation of sRNA transcript abundance was observed in northern blot experiments. Screening of GAS M1T1 was conducted with cells grown in complex medium, whereas the expression experiments in this study were performed with GAS M49 grown in CDM. Consequently, we found only a limited overlap between the two microarray screens (Figure3B). Twelve sRNA candidates were detected in both strains.
Analysis of common motifs in the GAS M49 sRNA population
Regulated expression of sRNA genes in GAS M49
Regulation of S. pyogenes M49 sRNA gene expression dependent on culture medium and growth phase
THY vs CDM exponential growth phase
THY vs CDM stationary growth phase
BHI vs CDM exponential growth phase
BHI vs CDM stationary growth phase
Stationary vs exponential in CDM
Validation of sRNA expression by qRT-PCR and northern blot analyses
Expression of sRNA candidates by GAS M49 was tested by gene-specific reverse transcription followed by real-time PCR analysis. Experimental expression validation was performed for the RNAs FASX and sRNASpy490822 (CRISPR1), for the sRNA scan7, which was predicted using sRNAscanner, and for three more candidates identified by tiling arrays in this study (sRNASpy490380c, sRNASpy490957c, and sRNASpy491311c). The expression of the candidates in GAS M49 was verified. Moreover, we confirmed the orientation of the sRNA genes by employing single gene-specific primers for the reverse transcription reaction. Three reactions were performed in parallel: one including the forward primer, one including the reverse primer, and one without any primers. Signals were only detected in samples containing the primer complementary to the coding strand of the respective gene (data not shown). We compared sRNA expression of GAS M49 cultured in CDM medium and THY broth throughout growth (Figure5). 5S RNA was used as an internal control for normalization, and was expressed in comparable amounts under all conditions tested in this experiment. Expression of fasX was equivalent during growth in CDM and THY (Figure5A). FASX was down-regulated in the stationary phase, an observation that confirmed the array data and previously published results from northern blot analyses. We did not detect strong regulation of CRISPR1, sRNASpy490957c, or sRNASpy491311c during growth in CDM (Figure5B). In contrast, scan7 was highly up-regulated in the stationary phase (Figure5B). The expression of sRNASpy490380c was much higher in THY compared to CDM (almost 100-fold, data not shown). During growth in CDM, no changes in the low level expression of sRNASpy490380c were observed (data not shown).
Bacterial gene regulation by sRNAs has gained a lot of attention in recent years, because it plays an important role in many cellular processes, including response to environmental changes, growth, and pathogenesis. There is an intriguingly large diversity of regulatory mechanisms, including cis- and trans-acting sRNAs, untranslated regions, and riboswitches. Some sRNA molecules act as repressors of translation and destabilize mRNA transcripts, but others act by activating and stabilizing target mRNAs[30–32]. One of the best characterized sRNAs in GAS is FASX, which is involved in virulence-related gene regulation[17, 33]. Knock-out mutants of fasX show a reduced expression of secreted virulence factors such as streptokinase and streptolysin S. The mechanism for streptokinase gene (ska) expression control is the stabilization of the ska transcript. Lack of FASX-ska-mRNA-interaction in the fasX deletion mutant decreased transcript levels, and consequently decreased streptokinase protein abundance.
A second example of a regulatory RNA in GAS is the untranslated mRNA of the streptococcal pleiotropic effect locus (pel), which contains sagA, the structural gene for streptolysin S. This region was described as a positive regulator of important streptococcal virulence factors, including M-protein, Sic, and SpeB. Strain specificity of PEL function is indicated by the fact that emm transcription was not affected in a sagA-deficient mutant with a M6 background. Similar results have been obtained in GAS M1 and M18 Tn916 sagA mutant strains. Additionally, pel deletion mutant analysis of four M1T1 GAS isolates did not identify any regulatory function for the pel sRNA in this serotype.
Another, more recently described untranslated RNA with influence on streptococcal virulence is the 4.5S RNA, a component of the bacterial signal recognition particle (SRP). While the 4.5S RNA gene is not essential, mutation impairs bacterial growth, lowers virulence factor secretion, and reduces virulence in a mouse infection model.
Recently, several whole genome sRNA screens in Gram-positive bacteria, employing either tiling array or next generation sequencing approaches, revealed an unexpected number of potential sRNAs in several pathogenic species[38–42]. In this context, it is likely that GAS expresses more sRNAs responsible for virulence gene expression control. One whole-genome intergenic tiling array screen of GAS M1T1 identified approximately 40 sRNAs that were expressed during the exponential growth phase in cells cultivated in THY complex medium. The GAS M49 sRNAome in the present study was determined using cells grown in CDM. From 55 putative sRNAs in GAS M49, only 12 were detected previously in the GAS M1T1 screen (Figure3B). This result is in accord with the concept that sRNA expression is serotype-dependent and regulated by environmental stimuli. Consequently, we detected media- and growth-phase-dependent sRNA gene regulation in the tiling array expression analysis, or by qRT-PCR of selected candidate genes. It would be interesting to monitor sRNA gene expression regulation under infection-relevant conditions.
Clustered, regularly interspaced, short palindromic repeat (CRISPR) loci represent an adaptive RNA-based immune system that protects bacteria and archaea from horizontal transfer of phage and plasmid DNA. Among the putative sRNA genes detected in GAS M49 by the tiling array approach, two sequences were categorised by the Rfam prediction program as CRISPR-related RNAs (Table1). sRNASpy490822 and sRNASpy490827 are encoded by the system II (Nmeni/CASS4 subtype) CRISPR/Cas locus, which was characterized recently by differential RNA sequencing in GAS SF370 (M1 serotype). Our data suggest that this locus is also active in GAS M49. Expression of sRNASpy490822 was confirmed by RT-PCR on the opposite strand of the CRISPR-associated genes under all conditions tested in this study. This transcript corresponds to the trans-activating CRISPR RNA (tracrRNA), which is responsible for the maturation of CRISPR RNA in concert with RNase III and the CRISPR-associated Csn1 protein. A third CRISPR-related RNA detected in our expression screen, sRNASpy491206c, is encoded in the system I-C (Dvulg/CASS subtype) CRISPR/Cas locus, which is also conserved in streptococcal genomes. In contrast to our array data, this locus appeared to be silent in GAS SF370, where no expression was detected in the differential RNA sequencing approach. Even though the CRISPR loci are conserved throughout GAS genomes, the activity of different CRISPR subtypes appears to be serotype-specific.
In the early years of sRNA research, many bioinformatic prediction tools were developed. One of the most prominent programs was the SIPHT tool, which has been used for many bacterial species[45, 46]. However, comparison of the prediction results with the actual in vivo expression of sRNAs often revealed very little overlap between the different screening methods[20, 41, 47]. The reasons for this discrepancy may be the limitations of the prediction programs as well as the fact that not all sRNAs are expressed under all conditions. The development of sRNA prediction software with improved properties is on-going. We compared our tiling array data with the prediction results of two recently published bioinformatics tools, sRNAScanner and MOSES. As depicted in Figure3A, the overlap between the tiling array expression data and the sRNA predictions was low. From the 20 most probable candidates of the MOSES analysis, 25% were expressed in GAS M49, whereas 8% of the predicted sRNAScanner predictions were found in the array analysis. Even the overlap between the two bioinformatics data sets was low. The only sRNA that was detected in all three screens was the previously characterized sRNA FASX. These results strongly suggest that a comprehensive analysis of bacterial genomes requires the combination of mathematical predictions with the collection of expression data. In the long term, testing of different conditions, especially mimicking in vivo situations by employing infection models, might lead to an increased overlap of expression detection and bioinformatics analyses.
We present here the identification of 55 putative sRNAs in GAS M49 by an intergenic tiling array approach. The candidate sRNA genes were expressed during growth in CDM. Forty-two of the RNAs were novel, whereas 13 RNAs have been described previously. The sequences of most of the candidates were conserved over streptococcal genomes. However, comparison of our GAS M49 sRNA expression data to another array analysis of a GAS M1 strain, and to two in silico screening methods, revealed little overlap between the different approaches. Thus, the investigation of several conditions and the combination of screening tools will be necessary to gain a comprehensive understanding of the abundance of sRNAs in GAS. The identification of novel differentially expressed sRNA genes will enhance our understanding of virulence related gene regulation in GAS. To account for specific expression patterns of putative sRNAs, infection relevant conditions combined with next generation RNA sequencing should be employed to investigate sRNA dependent regulatory networks in GAS.
Bacterial strains and culture conditions
GAS serotype M49 strain 591, a clinical isolate from a skin infection, was kindly provided by R. Lütticken (Aachen, Germany). The GAS strain was cultured in chemically defined medium (CDM), Todd-Hewitt broth (Invitrogen) supplemented with 0.5% yeast extract (THY; Invitrogen), or Brain-Heart-Infusion medium (BHI; Oxoid), as indicated, at 37°C under a 5% CO2/20% O2 atmosphere.
Total bacterial RNA from cultures grown to exponential and stationary phase of growth was isolated using the FastRNAProBlue Kit from MP Biomedicals according to manufacturer’s instructions. The purified total RNA was digested with DNaseI (Ambion) to remove remaining traces of chromosomal DNA. The RNA preparation was treated with 10 U of DNase1 for 30 min at 37°C. The enzyme was subsequently heat inactivated at 72°C for 5 min.
Enrichment of small RNAs
Five micrograms of total RNA were fractionated using the Ambion FlashPAGE Fractionator, Ambion FlashPAGE Precast Gels, and the Ambion FlashPAGE Buffer Kit, following manufacturer’s instructions. To collect the fraction of RNA molecules <200 nucleotides in length, the protocol was modified by increasing the running time from 12 min to up to 45 min at 75 V.
The small RNA fraction was ethanol-precipitated overnight at −20°C. The RNA was pelleted by centrifugation, dissolved in nuclease-free water, and labelled using the Ambion mirVana miRNA Labeling Kit following the manufacturer’s instructions. In brief, this kit involves two main steps; the 3′ amine-modified tailing reaction, and labelling with NHS-esters. Poly(A) Polymerase and a mixture of unmodified and amine-modified nucleotides were used to add a 20–50 nucleotide tail to the 3′ end of each RNA molecule in the sample. The amine-modified RNA molecules were purified and coupled to amine-reactive labelled biotin moieties as NHS-esters.
Design and synthesis of microfluidic microarrays
We used a microfluidic biochip (Febit Biomed) consisting of eight independent reaction chambers, the arrays, enclosed in a cartridge for fully automated processing. Each array contains 15,625 features which are synthesized in situ inside the microchannels using the Geniom One technology (febit biomed). The 50mer probes were designed as a whole genome tiling array, covering the intergenic regions of the S. pyogenes NZ131 genome (NCBI accession number: NC_011375). The forward and reverse oriented probes were synthesized in separate arrays. Thus, two arrays per sample were used.
Microarray hybridization and detection
All hybridization and detection steps were carried out using a Geniom RT Analyzer (febit biomed). Hybridizations were performed overnight (16 hours) at 42°C. Subsequently, biotin was detected with streptavidin-phycoerythrin (SAPE). A signal amplification step was added using biotinylated anti-streptavidin antibodies (Vector Laboratories) and a second incubation with SAPE (Invitrogen). Signal detection using the appropriate filter set (Cy3) of the Geniom device employed the auto-exposure function of the Geniom software. The data discussed in this publication have been deposited in the NCBI Gene Expression Omnibus and are accessible through GEO Series accession number GSE31228 (http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE31228).
Microarray data analysis
Raw intensities were analysed and extracted using Geniom Wizard software (febit biomed) as a tab delimited text file. The data were then converted into a matrix, with rows corresponding to the features and columns corresponding to the different samples. Data analysis was performed using GeneSpring GX (version 11) software (Agilent Technologies). The array background was calculated as the median signal intensity of all “blank-control” features on the array. Data were background corrected and then normalized using quantile normalization. Following normalization, a quality control step was performed that removed all data sets with a correlation coefficient less than 0.9 compared to the corresponding biological replicates. Of the original four biological replicate data sets representing cells grown in CDM, at least three were included in the analysis. The remaining probes of the biological replicates required intensity values greater than 300 on all three arrays. Regions that showed signals on probes of both strands were manually removed following the primary analysis. The statistical significance of the determined signals was tested by unpaired student’s t-test with a false discovery rate of 5%. Resulting data were combined with gene information from the flanking coding regions. Terminators and promoters were predicted by TransTermHP (http://transterm.cbcb.umd.edu/tt/Streptococcus_pyogenes_NZ131.tt) and BDGP Neural Network Promoter Prediction, and BProm (http://www.SoftBerry.com), respectively. To investigate sRNA gene regulation, two biological replicates of growth experiments conducted in THY or BHI were included. Following data normalization, three-fold signal intensity differences between various conditions were determined using the GeneSpring GX (version 11) software (Agilent Technologies). A motif search was conducted using MEME Suite, followed by motif analyses using TOMTOM, (http://meme.sdsc.edu/meme/intro.html).
The transcriptional start sites of sRNA candidates were determined using 5′ RACE (Invitrogen) following the manufacturer’s instructions. Briefly, first strand cDNA was synthesized using gene-specific primers ( Additional file2). The original mRNA was enzymatically removed and the 3′ end of the cDNA was tailed with dCTP by terminal deoxynucleotidyl transferase (TdT). PCR amplification was performed with nested, sequence-specific primers and an anchor primer provided by the 5′ RACE system. Primers specific for the sRNA genes tested here are listed in Additional file2. Following amplification, PCR products were cloned into a TOPO-TA vector (Invitrogen) and sequenced (GATC Biotech AG).
Quantitative reverse transcription PCR
Acidic phenol-extracted, DNaseI-treated total RNA was reverse transcribed to generate cDNA using the First-Strand cDNA Synthesis Kit from Invitrogen following the protocol provided by the manufacturer. For gene-specific reverse transcription (RT), three reactions were performed: two strand-specific reactions with either one forward or one reverse primer, and one control reaction without any primer. Primers were designed based on the full genome sequence of S. pyogenes M49 strain NZ131 (NCBI accession number: NC011375) and are listed in Additional file3. Three independent RT experiments were performed and all subsequent PCR reactions were performed in triplicate. Primer efficiency was tested on genomic GAS M49 DNA prior to use in RT reactions. All cDNA products were amplified by PCR with two primers specific for the respective candidate sequences. Real time PCR amplification was performed with SYBR Green (Fermentas) using an ABI PRISM 7000 Sequence Detection system (Applied Biosystems). The level of 5S RNA gene transcription was used for normalization. Relative gene expression was determined by the ΔΔCT method.
Northern blot analyses
Total RNA was isolated during exponential (OD600 = 0.4), transitional (OD600 = 0.8), and stationary (OD600 =1.2) growth phases. RNA samples (10 μg per growth phase) were loaded onto an 8% TBE-Urea polyacrylamide gel and separated by electrophoresis. Size standards (Ultra Low Range Ladder, Fermentas) were loaded on the same gel. RNA was electroblotted onto positively charged nylon membranes (Ambion), UV cross-linked, and probed overnight with a probe complementary to a candidate sRNA. Probes were generated by PCR with the same primers as used for the PCR reaction in qRT-PCR experiments (listed in Additional file3). Probes were labelled with biotin prior to hybridization (Brightstar psoralen-biotin labeling kit, Ambion). A BrightStar BioDetect Kit (Ambion) was used for detection, and autoradiography films were exposed to the luminescent blots.
The work of NP and BK was supported by a BMBF grant in the framework of the ERA-Net PathoGenoMics 2 program (FKZ 0315437B). The work of AB and TH2 was supported by a BMBF grant ERA-NET Pathogenomics Network to the sncRNAomics project 62080061 to TH2.
- Brantl S: Bacterial chromosome-encoded small regulatory RNAs. Future Microbiol. 2009, 4: 85-103. 10.2217/17460922.214.171.124.View ArticlePubMedGoogle Scholar
- Narberhaus F, Vogel J: Regulatory RNAs in prokaryotes: here, there and everywhere. Mol Microbiol. 2009, 74: 261-269. 10.1111/j.1365-2958.2009.06869.x.View ArticlePubMedGoogle Scholar
- Waters LS, Storz G: Regulatory RNAs in bacteria. Cell. 2009, 136: 615-628. 10.1016/j.cell.2009.01.043.PubMed CentralView ArticlePubMedGoogle Scholar
- Frohlich KS, Vogel J: Activation of gene expression by small RNA. Curr Opin Microbiol. 2009, 12: 674-682. 10.1016/j.mib.2009.09.009.View ArticlePubMedGoogle Scholar
- Podkaminski D, Vogel J: Small RNAs promote mRNA stability to activate the synthesis of virulence factors. Mol Microbiol. 2010, 78: 1327-1331. 10.1111/j.1365-2958.2010.07428.x.View ArticlePubMedGoogle Scholar
- Gorke B, Vogel J: Noncoding RNA control of the making and breaking of sugars. Genes Dev. 2008, 22: 2914-2925. 10.1101/gad.1717808.View ArticlePubMedGoogle Scholar
- Heeb S, Valverde C, Gigot-Bonnefoy C, Haas D: Role of the stress sigma factor RpoS in GacA/RsmA-controlled secondary metabolism and resistance to oxidative stress in Pseudomonas fluorescens CHA0. FEMS Microbiol Lett. 2005, 243: 251-258. 10.1016/j.femsle.2004.12.008.View ArticlePubMedGoogle Scholar
- Heidrich N, Chinali A, Gerth U, Brantl S: The small untranslated RNA SR1 from the Bacillus subtilis genome is involved in the regulation of arginine catabolism. Mol Microbiol. 2006, 62: 520-536. 10.1111/j.1365-2958.2006.05384.x.View ArticlePubMedGoogle Scholar
- Gottesman S, McCullen CA, Guillier M, Vanderpool CK, Majdalani N, Benhammou J, Thompson KM, FitzGerald PC, Sowa NA, FitzGerald DJ: Small RNA regulators and the bacterial response to stress. Cold Spring Harb Symp Quant Biol. 2006, 71: 1-11. 10.1101/sqb.2006.71.016.PubMed CentralView ArticlePubMedGoogle Scholar
- Vanderpool CK, Gottesman S: Noncoding RNAs at the membrane. Nat Struct Mol Biol. 2005, 12: 285-286. 10.1038/nsmb0405-285.View ArticlePubMedGoogle Scholar
- Papenfort K, Vogel J: Regulatory RNA in bacterial pathogens. Cell Host Microbe. 2010, 8: 116-127. 10.1016/j.chom.2010.06.008.View ArticlePubMedGoogle Scholar
- Livny J, Brencic A, Lory S, Waldor MK: Identification of 17 Pseudomonas aeruginosa sRNAs and prediction of sRNA-encoding genes in 10 diverse pathogens using the bioinformatic tool sRNAPredict2. Nucleic Acids Res. 2006, 34: 3484-3493. 10.1093/nar/gkl453.PubMed CentralView ArticlePubMedGoogle Scholar
- Carapetis JR, Steer AC, Mulholland EK, Weber M: The global burden of group A streptococcal diseases. Lancet Infect Dis. 2005, 5: 685-694. 10.1016/S1473-3099(05)70267-X.View ArticlePubMedGoogle Scholar
- Cole JN, Barnett TC, Nizet V, Walker MJ: Molecular insight into invasive group A streptococcal disease. Nat Rev Microbiol. 2011, 9: 724-736. 10.1038/nrmicro2648.View ArticlePubMedGoogle Scholar
- Fiedler T, Sugareva V, Patenge N, Kreikemeyer B: Insights into Streptococcus pyogenes pathogenesis from transcriptome studies. Future Microbiol. 2010, 5: 1675-1694. 10.2217/fmb.10.128.View ArticlePubMedGoogle Scholar
- Le RA, Charpentier E: Small RNAs in streptococci. RNA Biol. 2012, 9: 414-426. 10.4161/rna.20104.View ArticleGoogle Scholar
- Kreikemeyer B, Boyle MD, Buttaro BA, Heinemann M, Podbielski A: Group A streptococcal growth phase-associated virulence factor regulation by a novel operon (Fas) with homologies to two-component-type regulators requires a small RNA molecule. Mol Microbiol. 2001, 39: 392-406. 10.1046/j.1365-2958.2001.02226.x.View ArticlePubMedGoogle Scholar
- Roberts SA, Scott JR: RivR and the small RNA RivX: the missing links between the CovR regulatory cascade and the Mga regulon. Mol Microbiol. 2007, 66: 1506-1522.PubMedGoogle Scholar
- Raasch P, Schmitz U, Patenge N, Vera J, Kreikemeyer B, Wolkenhauer O: Non-coding RNA detection methods combined to improve usability, reproducibility and precision. BMC Bioinformatics. 2010, 11: 491-10.1186/1471-2105-11-491.PubMed CentralView ArticlePubMedGoogle Scholar
- Perez N, Trevino J, Liu Z, Ho SC, Babitzke P, Sumby P: A genome-wide analysis of small regulatory RNAs in the human pathogen group A Streptococcus. PLoS One. 2009, 4: e7668-10.1371/journal.pone.0007668.PubMed CentralView ArticlePubMedGoogle Scholar
- Sridhar J, Sambaturu N, Sabarinathan R, Ou HY, Deng Z, Sekar K, Rafi ZA, Rajakumar K: sRNAscanner: a computational tool for intergenic small RNA detection in bacterial genomes. PLoS One. 2010, 5: e11970-10.1371/journal.pone.0011970.PubMed CentralView ArticlePubMedGoogle Scholar
- Gardner PP, Daub J, Tate JG, Nawrocki EP, Kolbe DL, Lindgreen S, Wilkinson AC, Finn RD, Griffiths-Jones S, Eddy SR, et al, et al: Rfam: updates to the RNA families database. Nucleic Acids Res. 2009, 37: D136-D140. 10.1093/nar/gkn766.PubMed CentralView ArticlePubMedGoogle Scholar
- Tillier ER, Collins RA: The contributions of replication orientation, gene direction, and signal sequences to base-composition asymmetries in bacterial genomes. J Mol Evol. 2000, 50: 249-257.PubMedGoogle Scholar
- Rocha EP: The replication-related organization of bacterial genomes. Microbiology. 2004, 150: 1609-1627. 10.1099/mic.0.26974-0.View ArticlePubMedGoogle Scholar
- Rutherford K, Parkhill J, Crook J, Horsnell T, Rice P, Rajandream MA, Barrell B: Artemis: sequence visualization and annotation. Bioinformatics. 2000, 16: 944-945. 10.1093/bioinformatics/16.10.944.View ArticlePubMedGoogle Scholar
- Facklam R: What happened to the streptococci: overview of taxonomic and nomenclature changes. Clin Microbiol Rev. 2002, 15: 613-630. 10.1128/CMR.15.4.613-630.2002.PubMed CentralView ArticlePubMedGoogle Scholar
- Bailey TL, Boden M, Buske FA, Frith M, Grant CE, Clementi L, Ren J, Li WW, Noble WS: MEME SUITE: tools for motif discovery and searching. Nucleic Acids Res. 2009, 37: W202-W208. 10.1093/nar/gkp335.PubMed CentralView ArticlePubMedGoogle Scholar
- Gupta S, Stamatoyannopoulos JA, Bailey TL, Noble WS: Quantifying similarity between motifs. Genome Biol. 2007, 8: R24-10.1186/gb-2007-8-2-r24.PubMed CentralView ArticlePubMedGoogle Scholar
- Deltcheva E, Chylinski K, Sharma CM, Gonzales K, Chao Y, Pirzada ZA, Eckert MR, Vogel J, Charpentier E: CRISPR RNA maturation by trans-encoded small RNA and host factor RNase III. Nature. 2011, 471: 602-607. 10.1038/nature09886.PubMed CentralView ArticlePubMedGoogle Scholar
- Storz G, Vogel J, Wassarman KM: Regulation by small RNAs in bacteria: expanding frontiers. Mol Cell. 2011, 43: 880-891. 10.1016/j.molcel.2011.08.022.PubMed CentralView ArticlePubMedGoogle Scholar
- Thomason MK, Storz G: Bacterial antisense RNAs: how many are there, and what are they doing?. Annu Rev Genet. 2010, 44: 167-188. 10.1146/annurev-genet-102209-163523.PubMed CentralView ArticlePubMedGoogle Scholar
- Gottesman S, Storz G: Bacterial small RNA regulators: versatile roles and rapidly evolving variations. Cold Spring Harb Perspect Biol. 2011, 3: 10.1101/cshperspect.a003798.Google Scholar
- Ramirez-Pena E, Trevino J, Liu Z, Perez N, Sumby P: The group A Streptococcus small regulatory RNA FasX enhances streptokinase activity by increasing the stability of the ska mRNA transcript. Mol Microbiol. 2010, 78: 1332-1347. 10.1111/j.1365-2958.2010.07427.x.PubMed CentralView ArticlePubMedGoogle Scholar
- Mangold M, Siller M, Roppenser B, Vlaminckx BJ, Penfound TA, Klein R, Novak R, Novick RP, Charpentier E: Synthesis of group A streptococcal virulence factors is controlled by a regulatory RNA molecule. Mol Microbiol. 2004, 53: 1515-1527. 10.1111/j.1365-2958.2004.04222.x.View ArticlePubMedGoogle Scholar
- Biswas I, Germon P, McDade K, Scott JR: Generation and surface localization of intact M protein in Streptococcus pyogenes are dependent on sagA. Infect Immun. 2001, 69: 7029-7038. 10.1128/IAI.69.11.7029-7038.2001.PubMed CentralView ArticlePubMedGoogle Scholar
- Betschel SD, Borgia SM, Barg NL, Low DE, De Azavedo JC: Reduced virulence of group A streptococcal Tn916 mutants that do not produce streptolysin S. Infect Immun. 1998, 66: 1671-1679.PubMed CentralPubMedGoogle Scholar
- Trevino J, Perez N, Sumby P: The 4.5S RNA component of the signal recognition particle is required for group A Streptococcus virulence. Microbiology. 2010, 156: 1342-1350. 10.1099/mic.0.036558-0.PubMed CentralView ArticlePubMedGoogle Scholar
- Beaume M, Hernandez D, Docquier M, Delucinge-Vivier C, Descombes P, Francois P: Orientation and expression of methicillin-resistant Staphylococcus aureus small RNAs by direct multiplexed measurements using the nCounter of NanoString technology. J Microbiol Methods. 2011, 84: 327-334. 10.1016/j.mimet.2010.12.025.View ArticlePubMedGoogle Scholar
- Chen Y, Indurthi DC, Jones SW, Papoutsakis ET: Small RNAs in the genus Clostridium. MBio. 2011, 2: e00340-10.PubMed CentralPubMedGoogle Scholar
- Kumar R, Shah P, Swiatlo E, Burgess SC, Lawrence ML, Nanduri B: Identification of novel non-coding small RNAs from Streptococcus pneumoniae TIGR4 using high-resolution genome tiling arrays. BMC Genomics. 2010, 11: 350-10.1186/1471-2164-11-350.PubMed CentralView ArticlePubMedGoogle Scholar
- Mraheil MA, Billion A, Mohamed W, Mukherjee K, Kuenne C, Pischimarov J, Krawitz C, Retey J, Hartsch T, Chakraborty T, et al, et al: The intracellular sRNA transcriptome of Listeria monocytogenes during growth in macrophages. Nucleic Acids Res. 2011, 39: 4235-4248. 10.1093/nar/gkr033.PubMed CentralView ArticlePubMedGoogle Scholar
- Tsui HC, Mukherjee D, Ray VA, Sham LT, Feig AL, Winkler ME: Identification and characterization of noncoding small RNAs in Streptococcus pneumoniae serotype 2 strain D39. J Bacteriol. 2010, 192: 264-279. 10.1128/JB.01204-09.PubMed CentralView ArticlePubMedGoogle Scholar
- Marraffini LA, Sontheimer EJ: Self versus non-self discrimination during CRISPR RNA-directed immunity. Nature. 2010, 463: 568-571. 10.1038/nature08703.PubMed CentralView ArticlePubMedGoogle Scholar
- Haft DH, Selengut J, Mongodin EF, Nelson KE: A guild of 45 CRISPR-associated (Cas) protein families and multiple CRISPR/Cas subtypes exist in prokaryotic genomes. PLoS Comput Biol. 2005, 1: e60-10.1371/journal.pcbi.0010060.PubMed CentralView ArticlePubMedGoogle Scholar
- Livny J, Fogel MA, Davis BM, Waldor MK: sRNAPredict: an integrative computational approach to identify sRNAs in bacterial genomes. Nucleic Acids Res. 2005, 33: 4096-4105. 10.1093/nar/gki715.PubMed CentralView ArticlePubMedGoogle Scholar
- Livny J, Waldor MK: Identification of small RNAs in diverse bacterial species. Curr Opin Microbiol. 2007, 10: 96-101. 10.1016/j.mib.2007.03.005.View ArticlePubMedGoogle Scholar
- Arnvig KB, Young DB: Identification of small RNAs in Mycobacterium tuberculosis. Mol Microbiol. 2009, 73: 397-408. 10.1111/j.1365-2958.2009.06777.x.PubMed CentralView ArticlePubMedGoogle Scholar
- van de Rijn I, Kessler RE: Growth characteristics of group A streptococci in a new chemically defined medium. Infect Immun. 1980, 27: 444-448.PubMed CentralPubMedGoogle Scholar
- Baum M, Bielau S, Rittner N, Schmid K, Eggelbusch K, Dahms M, Schlauersbach A, Tahedl H, Beier M, Guimil R, et al, et al: Validation of a novel, fully integrated and flexible microarray benchtop facility for gene expression profiling. Nucleic Acids Res. 2003, 31: e151-10.1093/nar/gng151.PubMed CentralView ArticlePubMedGoogle Scholar
- Edgar R, Domrachev M, Lash AE: Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002, 30: 207-210. 10.1093/nar/30.1.207.PubMed CentralView ArticlePubMedGoogle Scholar
- Bolstad BM, Irizarry RA, Astrand M, Speed TP: A comparison of normalization methods for high density oligonucleotide array data based on variance and bias. Bioinformatics. 2003, 19: 185-193. 10.1093/bioinformatics/19.2.185.View ArticlePubMedGoogle Scholar
- Kingsford CL, Ayanbule K, Salzberg SL: Rapid, accurate, computational discovery of Rho-independent transcription terminators illuminates their relationship to DNA uptake. Genome Biol. 2007, 8: R22-10.1186/gb-2007-8-2-r22.PubMed CentralView ArticlePubMedGoogle Scholar
- Reese MG: Application of a time-delay neural network to promoter annotation in the Drosophila melanogaster genome. Comput Chem. 2001, 26: 51-56. 10.1016/S0097-8485(01)00099-7.View ArticlePubMedGoogle Scholar
- Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001, 25: 402-408. 10.1006/meth.2001.1262.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.