Riboswitch Discovery by Combining RNA-Seq and Genome-Wide Identification of Transcriptional Start Sites

Extensive Reannotation of the Genome of the Model Streptomycete Streptomyces lividans TK24 Based on Transcriptome and Proteome Information

Streptomyces lividans TK24 is a relevant Gram-positive soil inhabiting bacterium and one of the model organisms of the genus Streptomyces. It is known for its potential to produce secondary metabolites, antibiotics, and other industrially relevant products. S. lividans TK24 is the plasmid-free derivative of S. lividans 66 and a close genetic relative of the strain Streptomyces coelicolor A3(2). In this study, we used transcriptome and proteome data to improve the annotation of the S. lividans TK24 genome. The RNA-seq data of primary 5'-ends of transcripts were used to determine transcription start sites (TSS) in the genome. We identified 5,424 TSS, of which 4,664 were assigned to annotated CDS and ncRNAs, 687 to antisense transcripts distributed between 606 CDS and their UTRs, 67 to tRNAs, and 108 to novel transcripts and CDS. Using the TSS data, the promoter regions and their motifs were analyzed in detail, revealing a conserved -10 (TAnnnT) and a weakly conserved -35 region (nTGACn). The analysis of the 5' untranslated region (UTRs) of S. lividans TK24 revealed 17% leaderless transcripts. Several cis-regulatory elements, like riboswitches or attenuator structures could be detected in the 5'-UTRs. The S. lividans TK24 transcriptome contains at least 929 operons. The genome harbors 27 secondary metabolite gene clusters of which 26 could be shown to be transcribed under at least one of the applied conditions. Comparison of the reannotated genome with that of the strain Streptomyces coelicolor A3(2) revealed a high degree of similarity. This study presents an extensive reannotation of the S. lividans TK24 genome based on transcriptome and proteome analyses. The analysis of TSS data revealed insights into the promoter structure, 5'-UTRs, cis-regulatory elements, attenuator structures and novel transcripts, like small RNAs. Finally, the repertoire of secondary metabolite gene clusters was examined. These data provide a basis for future studies regarding gene characterization, transcriptional regulatory networks, and usage as a secondary metabolite producing strain.

Identification of RNAs bound by Hfq reveals widespread RNA partners and a sporulation regulator in the human pathogen Clostridioides difficile

Noncoding RNAs (ncRNA) have emerged as important components of regulatory networks governing bacterial physiology and virulence. Previous deep-sequencing analysis identified a large diversity of ncRNAs in the human enteropathogen Clostridioides (Clostridium) difficile. Some of them are trans-encoded RNAs that could require the RNA chaperone protein Hfq for their action. Recent analysis suggested a pleiotropic role of Hfq in C. difficile with the most pronounced effect on sporulation, a key process during the infectious cycle of this pathogen. However, a global view of RNAs interacting with C. difficile Hfq is missing. In the present study, we performed RNA immunoprecipitation high-throughput sequencing (RIP-Seq) to identify Hfq-associated RNAs in C. difficile. Our work revealed a large set of Hfq-interacting mRNAs and ncRNAs, including mRNA leaders and coding regions, known and potential new ncRNAs. In addition to trans-encoded RNAs, new categories of Hfq ligands were found including cis-antisense RNAs, riboswitches and CRISPR RNAs. ncRNA-mRNA and ncRNA-ncRNA pairings were postulated through computational predictions. Investigation of one of the Hfq-associated ncRNAs, RCd1, suggests that this RNA contributes to the control of late stages of sporulation in C. difficile. Altogether, these data provide essential molecular basis for further studies of post-transcriptional regulatory network in this enteropathogen.

Conserved and specific features of Streptococcus pyogenes and Streptococcus agalactiae transcriptional landscapes

BACKGROUND

The human pathogen Streptococcus pyogenes, or group A Streptococcus, is responsible for mild infections to life-threatening diseases. To facilitate the characterization of regulatory networks involved in the adaptation of this pathogen to its different environments and their evolution, we have determined the primary transcriptome of a serotype M1 S. pyogenes strain at single-nucleotide resolution and compared it with that of Streptococcus agalactiae, also from the pyogenic group of streptococci.

RESULTS

By using a combination of differential RNA-sequencing and oriented RNA-sequencing we have identified 892 transcription start sites (TSS) and 885 promoters in the S. pyogenes M1 strain S119. 8.6% of S. pyogenes mRNAs were leaderless, among which 81% were also classified as leaderless in S. agalactiae. 26% of S. pyogenes transcript 5' untranslated regions (UTRs) were longer than 60 nt. Conservation of long 5' UTRs with S. agalactiae allowed us to predict new potential regulatory sequences. In addition, based on the mapping of 643 transcript ends in the S. pyogenes strain S119, we constructed an operon map of 401 monocistrons and 349 operons covering 81.5% of the genome. One hundred fifty-six operons and 254 monocistrons retained the same organization, despite multiple genomic reorganizations between S. pyogenes and S. agalactiae. Genomic reorganization was found to more often go along with variable promoter sequences and 5' UTR lengths. Finally, we identified 117 putative regulatory RNAs, among which nine were regulated in response to magnesium concentration.

CONCLUSIONS

Our data provide insights into transcriptome evolution in pyogenic streptococci and will facilitate the analysis of genetic polymorphisms identified by comparative genomics in S. pyogenes.

Synthetic Biology of Small RNAs and Riboswitches

In bacteria and archaea, small RNAs (sRNAs) regulate complex networks through antisense interactions with target mRNAs in trans, and riboswitches regulate gene expression in cis based on the ability to bind small-molecule ligands. Although our understanding and characterization of these two important regulatory RNA classes is far from complete, these RNA-based mechanisms have proven useful for a wide variety of synthetic biology applications. Besides classic and contemporary applications in the realm of metabolic engineering and orthogonal gene control, this review also covers newer applications of regulatory RNAs as biosensors, logic gates, and tools to determine RNA-RNA interactions. A separate section focuses on critical insights gained and challenges posed by fundamental studies of sRNAs and riboswitches that should aid future development of synthetic regulatory RNAs.

Iron Regulation in Clostridioides difficile

The response to iron limitation of several bacteria is regulated by the ferric uptake regulator (Fur). The Fur-regulated transcriptional, translational and metabolic networks of the Gram-positive, pathogen Clostridioides difficile were investigated by a combined RNA sequencing, proteomic, metabolomic and electron microscopy approach. At high iron conditions (15 μM) the C. difficile fur mutant displayed a growth deficiency compared to wild type C. difficile cells. Several iron and siderophore transporter genes were induced by Fur during low iron (0.2 μM) conditions. The major adaptation to low iron conditions was observed for the central energy metabolism. Most ferredoxin-dependent amino acid fermentations were significantly down regulated (had, etf, acd, grd, trx, bdc, hbd). The substrates of these pathways phenylalanine, leucine, glycine and some intermediates (phenylpyruvate, 2-oxo-isocaproate, 3-hydroxy-butyryl-CoA, crotonyl-CoA) accumulated, while end products like isocaproate and butyrate were found reduced. Flavodoxin (fldX) formation and riboflavin biosynthesis (rib) were enhanced, most likely to replace the missing ferredoxins. Proline reductase (prd), the corresponding ion pumping RNF complex (rnf) and the reaction product 5-aminovalerate were significantly enhanced. An ATP forming ATPase (atpCDGAHFEB) of the F₀F₁-type was induced while the formation of a ATP-consuming, proton-pumping V-type ATPase (atpDBAFCEKI) was decreased. The [Fe-S] enzyme-dependent pyruvate formate lyase (pfl), formate dehydrogenase (fdh) and hydrogenase (hyd) branch of glucose utilization and glycogen biosynthesis (glg) were significantly reduced, leading to an accumulation of glucose and pyruvate. The formation of [Fe-S] enzyme carbon monoxide dehydrogenase (coo) was inhibited. The fur mutant showed an increased sensitivity to vancomycin and polymyxin B. An intensive remodeling of the cell wall was observed, Polyamine biosynthesis (spe) was induced leading to an accumulation of spermine, spermidine, and putrescine. The fur mutant lost most of its flagella and motility. Finally, the CRISPR/Cas and a prophage encoding operon were downregulated. Fur binding sites were found upstream of around 20 of the regulated genes. Overall, adaptation to low iron conditions in C. difficile focused on an increase of iron import, a significant replacement of iron requiring metabolic pathways and the restructuring of the cell surface for protection during the complex adaptation phase and was only partly directly regulated by Fur.

Tracking gene expression and oxidative damage of O2-stressed Clostridioides difficile by a multi-omics approach

Clostridioides difficile is the major pathogen causing diarrhea following antibiotic treatment. It is considered to be a strictly anaerobic bacterium, however, previous studies have shown a certain and strain-dependent oxygen tolerance. In this study, the model strain C. difficile 630Δerm was shifted to micro-aerobiosis and was found to stay growing to the same extent as anaerobically growing cells with only few changes in the metabolite pattern. However, an extensive change in gene expression was determined by RNA-Seq. The most striking adaptation strategies involve a change in the reductive fermentation pathways of the amino acids proline, glycine and leucine. But also a far-reaching restructuring in the carbohydrate metabolism was detected with changes in the phosphotransferase system (PTS) facilitated uptake of sugars and a repression of enzymes of glycolysis and butyrate fermentation. Furthermore, a temporary induction in the synthesis of cofactor riboflavin was detected possibly due to an increased demand for flavin mononucleotid (FMN) and flavin adenine dinucleotide (FAD) in redox reactions. However, biosynthesis of the cofactors thiamin pyrophosphate and cobalamin were repressed deducing oxidation-prone enzymes and intermediates in these pathways. Micro-aerobically shocked cells were characterized by an increased demand for cysteine and a thiol redox proteomics approach revealed a dramatic increase in the oxidative state of cysteine in more than 800 peptides after 15 min of micro-aerobic shock. This provides not only a catalogue of oxidation-prone cysteine residues in the C. difficile proteome but also puts the amino acid cysteine into a key position in the oxidative stress response. Our study suggests that tolerance of C. difficile towards O₂ is based on a complex and far-reaching adjustment of global gene expression which leads to only a slight change in phenotype.

Identification of c-di-GMP-Responsive Riboswitches

Bis-(3'-5')-cyclic dimeric guanosine monophosphate (c-di-GMP) is an important signaling molecule for community behavior control, cell morphogenesis, and virulence in bacteria. In addition to protein effectors, this second messenger binds RNA molecules that act as riboswitches to control target gene expression. In this chapter, we describe a method for experimental validation of the functionality of c-di-GMP-responsive riboswitches and the analysis of c-di-GMP control of target gene expression by qRT-PCR and Northern blot. This procedure can be used for the studies of in silico-predicted riboswitch candidates, as well as a targeted experimental approach for exploring the data from next-generation sequencing. The examples on the analysis of type I and type II c-di-GMP-responsive riboswitches in Clostridium difficile are provided to illustrate the application of the method.

Transcriptome Landscape of Mycobacterium smegmatis

The non-pathogenic bacterium Mycobacterium smegmatis mc²155 has been widely used as a model organism in mycobacterial research, yet a detailed study about its transcription landscape remains to be established. Here we report the transcriptome, expression profiles and transcriptional structures through growth-phase-dependent RNA sequencing (RNA-seq) as well as other related experiments. We found: (1) 2,139 transcriptional start sites (TSSs) in the genome-wide scale, of which eight samples were randomly selected and further verified by 5′-RACE; (2) 2,233 independent monocistronic or polycistronic mRNAs in the transcriptome within the operon/sub-operon structures which are classified into five groups; (3) 47.50% (1016/2139) genes were transcribed into leaderless mRNAs, with the TSSs of 41.3% (883/2139) mRNAs overlapping with the first base of the annotated start codon. Initial amino acids of MSMEG_4921 and MSMEG_6422 proteins were identified by Edman degradation, indicating the presence of distinctive widespread leaderless features in M. smegmatis mc²155. (4) 150 genes with potentially wrong structural annotation, of which 124 proposed genes have been corrected; (5) eight highly active promoters, with their activities further determined by β-galactosidase assays. These data integrated the transcriptional landscape to genome information of model organism mc²155 and lay a solid foundation for further works in Mycobacterium.

Genome improvement of the acarbose producer Actinoplanes sp. SE50/110 and annotation refinement based on RNA-seq analysis

Actinoplanes sp. SE50/110 is the natural producer of acarbose, which is used in the treatment of diabetes mellitus type II. However, until now the transcriptional organization and regulation of the acarbose biosynthesis are only understood rudimentarily. The genome sequence of Actinoplanes sp. SE50/110 was known before, but was resequenced in this study to remove assembly artifacts and incorrect base callings. The annotation of the genome was refined in a multi-step approach, including modern bioinformatic pipelines, transcriptome and proteome data. A whole transcriptome RNA-seq library as well as an RNA-seq library enriched for primary 5'-ends were used for the detection of transcription start sites, to correct tRNA predictions, to identify novel transcripts like small RNAs and to improve the annotation through the correction of falsely annotated translation start sites. The transcriptome data sets were also applied to identify 31 cis-regulatory RNA structures, such as riboswitches or RNA thermometers as well as three leaderless transcribed short peptides found in putative attenuators upstream of genes for amino acid biosynthesis. The transcriptional organization of the acarbose biosynthetic gene cluster was elucidated in detail and fourteen novel biosynthetic gene clusters were suggested. The accurate genome sequence and precise annotation of the Actinoplanes sp. SE50/110 genome will be the foundation for future genetic engineering and systems biology studies.

Regulation of PI-2b Pilus Expression in Hypervirulent Streptococcus agalactiae ST-17 BM110

The widely spread Streptococcus agalactiae (also known as Group B Streptococcus, GBS) “hypervirulent” ST17 clone is strongly associated with neonatal meningitis. The PI-2b locus is mainly found in ST17 strains but is also present in a few non ST17 human isolates such as the ST-7 prototype strain A909. Here, we analysed the expression of the PI-2b pilus in the ST17 strain BM110 as compared to the non ST17 A909. Comparative genome analyses revealed the presence of a 43-base pair (bp) hairpin-like structure in the upstream region of PI-2b operon in all 26 ST17 genomes, which was absent in the 8 non-ST17 strains carrying the PI-2b locus. Deletion of this 43-bp sequence in strain BM110 resulted in a 3- to 5-fold increased transcription of PI-2b. Characterization of PI-2b promoter region in A909 and BM110 strains was carried out by RNAseq, primer extension, qRT-PCR and transcriptional fusions with gfp as reporter gene. Our results indicate the presence of a single promoter (Ppi2b) with a transcriptional start site (TSS) mapped 37 bases upstream of the start codon of the first PI-2b gene. The large operon of 16 genes located upstream of PI-2b codes for the group B carbohydrate (also known as antigen B), a major constituent of the bacterial cell wall. We showed that the hairpin sequence located between antigen B and PI-2b operons is a transcriptional terminator. In A909, increased expression of PI-2b probably results from read-through transcription from antigen B operon. In addition, we showed that an extended 5’ promoter region is required for maximal transcription of gfp as a reporter gene in S. agalactiae from Ppi2b promoter. Gene reporter assays performed in Lactococcus lactis strain NZ9000, a related non-pathogenic Gram-positive species, revealed that GBS-specific regulatory factors are required to drive PI-2b transcription. PI-2b expression is up-regulated in the BM110ΔcovR mutant as compared to the parental BM110 strain, but this effect is probably indirect. Collectively, our results indicate that PI-2b expression is regulated in GBS ST17 strains, which may confer a selective advantage in the human host either by reducing host immune responses and/or increasing their dissemination potential.

Opportunities in the design and application of RNA for gene expression control

The past decade of synthetic biology research has witnessed numerous advances in the development of tools and frameworks for the design and characterization of biological systems. Researchers have focused on the use of RNA for gene expression control due to its versatility in sensing molecular ligands and the relative ease by which RNA can be modeled and designed compared to proteins. We review the recent progress in the field with respect to RNA-based genetic devices that are controlled through small molecule and protein interactions. We discuss new approaches for generating and characterizing these devices and their underlying components. We also highlight immediate challenges, future directions and recent applications of synthetic RNA devices in engineered biological systems.

Function of the CRISPR-Cas System of the Human Pathogen Clostridium difficile

Clostridium difficile is the cause of most frequently occurring nosocomial diarrhea worldwide. As an enteropathogen, C. difficile must be exposed to multiple exogenous genetic elements in bacteriophage-rich gut communities. CRISPR (clustered regularly interspaced short palindromic repeats)-Cas (CRISPR-associated) systems allow bacteria to adapt to foreign genetic invaders. Our recent data revealed active expression and processing of CRISPR RNAs from multiple type I-B CRISPR arrays in C. difficile reference strain 630. Here, we demonstrate active expression of CRISPR arrays in strain R20291, an epidemic C. difficile strain. Through genome sequencing and host range analysis of several new C. difficile phages and plasmid conjugation experiments, we provide evidence of defensive function of the CRISPR-Cas system in both C. difficile strains. We further demonstrate that C. difficile Cas proteins are capable of interference in a heterologous host, Escherichia coli. These data set the stage for mechanistic and physiological analyses of CRISPR-Cas-mediated interactions of important global human pathogen with its genetic parasites.

Clostridium difficile is the major cause of nosocomial infections associated with antibiotic therapy worldwide. To survive in bacteriophage-rich gut communities, enteropathogens must develop efficient systems for defense against foreign DNA elements. CRISPR-Cas systems have recently taken center stage among various anti-invader bacterial defense systems. We provide experimental evidence for the function of the C. difficile CRISPR system against plasmid DNA and bacteriophages. These data demonstrate the original features of active C. difficile CRISPR system and bring important insights into the interactions of this major enteropathogen with foreign DNA invaders during its infection cycle.

Single nucleotide resolution RNA-seq uncovers new regulatory mechanisms in the opportunistic pathogen Streptococcus agalactiae

Background

Streptococcus agalactiae, or Group B Streptococcus, is a leading cause of neonatal infections and an increasing cause of infections in adults with underlying diseases. In an effort to reconstruct the transcriptional networks involved in S. agalactiae physiology and pathogenesis, we performed an extensive and robust characterization of its transcriptome through a combination of differential RNA-sequencing in eight different growth conditions or genetic backgrounds and strand-specific RNA-sequencing.

Results

Our study identified 1,210 transcription start sites (TSSs) and 655 transcript ends as well as 39 riboswitches and cis-regulatory regions, 39 cis-antisense non-coding RNAs and 47 small RNAs potentially acting in trans. Among these putative regulatory RNAs, ten were differentially expressed in response to an acid stress and two riboswitches sensed directly or indirectly the pH modification. Strikingly, 15% of the TSSs identified were associated with the incorporation of pseudo-templated nucleotides, showing that reiterative transcription is a pervasive process in S. agalactiae. In particular, 40% of the TSSs upstream genes involved in nucleotide metabolism show reiterative transcription potentially regulating gene expression, as exemplified for pyrG and thyA encoding the CTP synthase and the thymidylate synthase respectively.

Conclusions

This comprehensive map of the transcriptome at the single nucleotide resolution led to the discovery of new regulatory mechanisms in S. agalactiae. It also provides the basis for in depth analyses of transcriptional networks in S. agalactiae and of the regulatory role of reiterative transcription following variations of intra-cellular nucleotide pools.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1583-4) contains supplementary material, which is available to authorized users.

Transcriptomic profiling of Yersinia pseudotuberculosis reveals reprogramming of the Crp regulon by temperature and uncovers Crp as a master regulator of small RNAs

One hallmark of pathogenic yersiniae is their ability to rapidly adjust their life-style and pathogenesis upon host entry. In order to capture the range, magnitude and complexity of the underlying gene control mechanisms we used comparative RNA-seq-based transcriptomic profiling of the enteric pathogen Y. pseudotuberculosis under environmental and infection-relevant conditions. We identified 1151 individual transcription start sites, multiple riboswitch-like RNA elements, and a global set of antisense RNAs and previously unrecognized trans-acting RNAs. Taking advantage of these data, we revealed a temperature-induced and growth phase-dependent reprogramming of a large set of catabolic/energy production genes and uncovered the existence of a thermo-regulated ‘acetate switch’, which appear to prime the bacteria for growth in the digestive tract. To elucidate the regulatory architecture linking nutritional status to virulence we also refined the CRP regulon. We identified a massive remodelling of the CRP-controlled network in response to temperature and discovered CRP as a transcriptional master regulator of numerous conserved and newly identified non-coding RNAs which participate in this process. This finding highlights a novel level of complexity of the regulatory network in which the concerted action of transcriptional regulators and multiple non-coding RNAs under control of CRP adjusts the control of Yersinia fitness and virulence to the requirements of their environmental and virulent life-styles.

Many bacterial pathogens cycle between environmental sources and mammalian hosts. Adaptation to the different natural habitats and host niches is achieved through complex regulatory networks which adjust synthesis of the large repertoire of crucial virulence factors and fitness determinants. To uncover underlying control circuits, we determined the first in-depth single-nucleotide resolution transcriptome of Yersinia. This revealed important novel genetic information, such as global locations of transcriptional start sites, non-coding RNAs, potential riboswitches and provided a set of virulence-relevant expression profiles, which constitute a valuable tool for the research community. The analysis further uncovered a temperature-induced global reprogramming of central metabolic functions, likely to support intestinal colonization of the pathogen. This is accompanied by a major reorganization of the CRP regulon, which involves a multitude of regulatory RNAs. The primary consequence is a fine-tuned, coordinated control of metabolism and virulence through a plethora of environmentally controlled regulatory RNAs allowing rapid adaptation and high flexibility during life-style changes.