Differential RNA-seq: the approach behind and the biological insight gained

Whole-genome Identification of Transcriptional Start Sites by Differential RNA-seq in Bacteria

Gene transcription in bacteria often starts some nucleotides upstream of the start codon. Identifying the specific Transcriptional Start Site (TSS) is essential for genetic manipulation, as in many cases upstream of the start codon there are sequence elements that are involved in gene expression regulation. Taken into account the classical gene structure, we are able to identify two kinds of transcriptional start site: primary and secondary. A primary transcriptional start site is located some nucleotides upstream of the translational start site, while a secondary transcriptional start site is located within the gene encoding sequence. Here, we present a step by step protocol for genome-wide transcriptional start sites determination by differential RNA-sequencing (dRNA-seq) using the enteric pathogen Shigella flexneri serotype 5a strain M90T as model. However, this method can be employed in any other bacterial species of choice. In the first steps, total RNA is purified from bacterial cultures using the hot phenol method. Ribosomal RNA (rRNA) is specifically depleted via hybridization probes using a commercial kit. A 5'-monophosphate-dependent exonuclease (TEX)-treated RNA library enriched in primary transcripts is then prepared for comparison with a library that has not undergone TEX-treatment, followed by ligation of an RNA linker adaptor of known sequence allowing the determination of TSS with single nucleotide precision. Finally, the RNA is processed for Illumina sequencing library preparation and sequenced as purchased service. TSS are identified by in-house bioinformatic analysis. Our protocol is cost-effective as it minimizes the use of commercial kits and employs freely available software.

The Transcriptomic Landscape of Cupriavidus metallidurans CH34 Acutely Exposed to Copper

Bacteria are increasingly used for biotechnological applications such as bioremediation, biorecovery, bioproduction, and biosensing. The development of strains suited for such applications requires a thorough understanding of their behavior, with a key role for their transcriptomic landscape. We present a thorough analysis of the transcriptome of Cupriavidus metallidurans CH34 cells acutely exposed to copper by tagRNA-sequencing. C. metallidurans CH34 is a model organism for metal resistance, and its potential as a biosensor and candidate for metal bioremediation has been demonstrated in multiple studies. Several metabolic pathways were impacted by Cu exposure, and a broad spectrum of metal resistance mechanisms, not limited to copper-specific clusters, was overexpressed. In addition, several gene clusters involved in the oxidative stress response and the cysteine-sulfur metabolism were induced. In total, 7500 transcription start sites (TSSs) were annotated and classified with respect to their location relative to coding sequences (CDSs). Predicted TSSs were used to re-annotate 182 CDSs. The TSSs of 2422 CDSs were detected, and consensus promotor logos were derived. Interestingly, many leaderless messenger RNAs (mRNAs) were found. In addition, many mRNAs were transcribed from multiple alternative TSSs. We observed pervasive intragenic TSSs both in sense and antisense to CDSs. Antisense transcripts were enriched near the 5' end of mRNAs, indicating a functional role in post-transcriptional regulation. In total, 578 TSSs were detected in intergenic regions, of which 35 were identified as putative small regulatory RNAs. Finally, we provide a detailed analysis of the main copper resistance clusters in CH34, which include many intragenic and antisense transcripts. These results clearly highlight the ubiquity of noncoding transcripts in the CH34 transcriptome, many of which are putatively involved in the regulation of metal resistance.

Genome-Wide Transcription Start Site Mapping and Promoter Assignments to a Sigma Factor in the Human Enteropathogen Clostridioides difficile

The emerging human enteropathogen Clostridioides difficile is the main cause of diarrhea associated with antibiotherapy. Regulatory pathways underlying the adaptive responses remain understudied and the global view of C. difficile promoter structure is still missing. In the genome of C. difficile 630, 22 genes encoding sigma factors are present suggesting a complex pattern of transcription in this bacterium. We present here the first transcriptional map of the C. difficile genome resulting from the identification of transcriptional start sites (TSS), promoter motifs and operon structures. By 5′-end RNA-seq approach, we mapped more than 1000 TSS upstream of genes. In addition to these primary TSS, this analysis revealed complex structure of transcriptional units such as alternative and internal promoters, potential RNA processing events and 5′ untranslated regions. By following an in silico iterative strategy that used as an input previously published consensus sequences and transcriptomic analysis, we identified candidate promoters upstream of most of protein-coding and non-coding RNAs genes. This strategy also led to refine consensus sequences of promoters recognized by major sigma factors of C. difficile. Detailed analysis focuses on the transcription in the pathogenicity locus and regulatory genes, as well as regulons of transition phase and sporulation sigma factors as important components of C. difficile regulatory network governing toxin gene expression and spore formation. Among the still uncharacterized regulons of the major sigma factors of C. difficile, we defined the SigL regulon by combining transcriptome and in silico analyses. We showed that the SigL regulon is largely involved in amino-acid degradation, a metabolism crucial for C. difficile gut colonization. Finally, we combined our TSS mapping, in silico identification of promoters and RNA-seq data to improve gene annotation and to suggest operon organization in C. difficile. These data will considerably improve our knowledge of global regulatory circuits controlling gene expression in C. difficile and will serve as a useful rich resource for scientific community both for the detailed analysis of specific genes and systems biology studies.

A high-resolution transcriptome map identifies small RNA regulation of metabolism in the gut microbe Bacteroides thetaiotaomicron

Bacteria of the genus Bacteroides are common members of the human intestinal microbiota and important degraders of polysaccharides in the gut. Among them, the species Bacteroides thetaiotaomicron has emerged as the model organism for functional microbiota research. Here, we use differential RNA sequencing (dRNA-seq) to generate a single-nucleotide resolution transcriptome map of B. thetaiotaomicron grown under defined laboratory conditions. An online browser, called ‘Theta-Base’ (www.helmholtz-hiri.de/en/datasets/bacteroides), is launched to interrogate the obtained gene expression data and annotations of ~4500 transcription start sites, untranslated regions, operon structures, and 269 noncoding RNA elements. Among the latter is GibS, a conserved, 145 nt-long small RNA that is highly expressed in the presence of N-acetyl-D-glucosamine as sole carbon source. We use computational predictions and experimental data to determine the secondary structure of GibS and identify its target genes. Our results indicate that sensing of N-acetyl-D-glucosamine induces GibS expression, which in turn modifies the transcript levels of metabolic enzymes.

Bacteroides thetaiotaomicron is a human gut microbe and an emergent model organism. Here, Ryan et al. generate single-nucleotide resolution RNA-seq data for this bacterium and map transcription start sites and noncoding RNAs, one of which modulates expression of metabolic enzymes.

Elucidation of the Biosynthetic Pathway of Vitamin B Groups and Potential Secondary Metabolite Gene Clusters Via Genome Analysis of a Marine Bacterium Pseudoruegeria sp. M32A2M

The symbiotic nature of the relationship between algae and marine bacteria is well-studied among the complex microbial interactions. The mutual profit between algae and bacteria occurs via nutrient and vitamin exchange. It is necessary to analyze the genome sequence of a bacterium to predict its symbiotic relationships. In this study, the genome of a marine bacterium, Pseudoruegeria sp. M32A2M, isolated from the south-eastern isles (GeoJe-Do) of South Korea, was sequenced and analyzed. A draft genome (91 scaffolds) of 5.5 Mb with a DNA G+C content of 62.4% was obtained. In total, 5,101 features were identified from gene annotation, and 4,927 genes were assigned to functional proteins. We also identified transcription core proteins, RNA polymerase subunits, and sigma factors. In addition, full flagella-related gene clusters involving the flagellar body, motor, regulator, and other accessory compartments were detected even though the genus Pseudoruegeria is known to comprise non-motile bacteria. Examination of annotated KEGG pathways revealed that Pseudoruegeria sp. M32A2M has the metabolic pathways for all seven vitamin Bs, including thiamin (vitamin B1), biotin (vitamin B7), and cobalamin (vitamin B12), which are necessary for symbiosis with vitamin B auxotroph algae. We also identified gene clusters for seven secondary metabolites including ectoine, homoserine lactone, beta-lactone, terpene, lasso peptide, bacteriocin, and nonribosomal proteins.

Integrative functional genomics decodes herpes simplex virus 1

The predicted 80 open reading frames (ORFs) of herpes simplex virus 1 (HSV-1) have been intensively studied for decades. Here, we unravel the complete viral transcriptome and translatome during lytic infection with base-pair resolution by computational integration of multi-omics data. We identify a total of 201 transcripts and 284 ORFs including all known and 46 novel large ORFs. This includes a so far unknown ORF in the locus deleted in the FDA-approved oncolytic virus Imlygic. Multiple transcript isoforms expressed from individual gene loci explain translation of the vast majority of ORFs as well as N-terminal extensions (NTEs) and truncations. We show that NTEs with non-canonical start codons govern the subcellular protein localization and packaging of key viral regulators and structural proteins. We extend the current nomenclature to include all viral gene products and provide a genome browser that visualizes all the obtained data from whole genome to single-nucleotide resolution.

Dynamic insights on transcription initiation and RNA processing during bacterial adaptation

Transcription initiation and RNA processing govern gene expression and enable bacterial adaptation by reshaping the RNA landscape. The aim of this study was to simultaneously observe these two fundamental processes in a transcriptome responding to an environmental signal. A controlled σ^E system in E. coli was coupled to our previously described tagRNA-seq method to yield process kinetics information. Changes in transcription initiation frequencies (TIF) and RNA processing frequencies (PF) were followed using 5' RNA tags. Changes in TIF showed a binary increased/decreased pattern that alternated between transcriptionally activated and repressed promoters, providing the bacterial population with transcriptional oscillation. PF variation fell into three categories of cleavage activity: (i) constant and independent of RNA levels, (ii) increased once RNA has accumulated, and (iii) positively correlated to changes in TIF. This work provides a comprehensive and dynamic view of major events leading to transcriptomic reshaping during bacterial adaptation. It unveils an interplay between transcription initiation and the activity of specific RNA cleavage sites. This study utilized a well-known genetic system to analyze fundamental processes and can serve as a blueprint for comprehensive studies that exploit the RNA metabolism to decipher and understand bacterial gene expression control.

Riboregulation in Nitrogen-Fixing Endosymbiotic Bacteria

Small non-coding RNAs (sRNAs) are ubiquitous components of bacterial adaptive regulatory networks underlying stress responses and chronic intracellular infection of eukaryotic hosts. Thus, sRNA-mediated regulation of gene expression is expected to play a major role in the establishment of mutualistic root nodule endosymbiosis between nitrogen-fixing rhizobia and legume plants. However, knowledge about this level of genetic regulation in this group of plant-interacting bacteria is still rather scarce. Here, we review insights into the rhizobial non-coding transcriptome and sRNA-mediated post-transcriptional regulation of symbiotic relevant traits such as nutrient uptake, cell cycle, quorum sensing, or nodule development. We provide details about the transcriptional control and protein-assisted activity mechanisms of the functionally characterized sRNAs involved in these processes. Finally, we discuss the forthcoming research on riboregulation in legume symbionts.

Prevalence of small base-pairing RNAs derived from diverse genomic loci

Small RNAs (sRNAs) that act by base-pairing have been shown to play important roles in fine-tuning the levels and translation of their target transcripts across a variety of model and pathogenic organisms. Work from many different groups in a wide range of bacterial species has provided evidence for the importance and complexity of sRNA regulatory networks, which allow bacteria to quickly respond to changes in their environment. However, despite the expansive literature, much remains to be learned about all aspects of sRNA-mediated regulation, particularly in bacteria beyond the well-characterized Escherichia coli and Salmonella enterica species. Here we discuss what is known, and what remains to be learned, about the identification of regulatory base-pairing RNAs produced from diverse genomic loci including how their expression is regulated. This article is part of a Special Issue entitled: RNA and gene control in bacteria edited by Dr. M. Guillier and F. Repoila.

Navigation through the twists and turns of RNA sequencing technologies: Application to bacterial regulatory RNAs

Discovered in the 1980s, small regulatory RNAs (sRNAs) are now considered key actors in virtually all aspects of bacterial physiology and virulence. Together with transcriptional and translational regulatory proteins, they integrate and often are hubs of complex regulatory networks, responsible for bacterial response/adaptation to various perceived stimuli. The recent development of powerful RNA sequencing technologies has facilitated the identification and characterization of sRNAs (length, structure and expression conditions) and their RNA targets in several bacteria. Nevertheless, it could be very difficult for non-experts to understand the advantages and drawbacks related to each offered option and, consequently, to make an informed choice. Therefore, the main goal of this review is to provide a guide to navigate through the twists and turns of high-throughput RNA sequencing technologies, with a specific focus on those applied to the study of sRNAs. This article is part of a Special Issue entitled: RNA and gene control in bacteria edited by Dr. M. Guillier and F. Repoila.

De novo transcriptome analysis of halotolerant bacterium Staphylococcus sp. strain P-TSB-70 isolated from East coast of India: In search of salt stress tolerant genes

In the present study, we identified salt stress tolerant genes from the marine bacterium Staphylococcus sp. strain P-TSB-70 through transcriptome sequencing. In favour of whole-genome transcriptome profiling of Staphylococcus sp. strain P-TSB-70 (GenBank Accn. No. KP117091) which tolerated upto 20% NaCl stress, the strain was cultured in the laboratory condition with 20% NaCl stress. Transcriptome analyses were performed by SOLiD4.0 sequencing technology from which 10280 and 9612 transcripts for control and treated, respectively, were obtained. The coverage per base (CPB) statistics were analyzed for both the samples. Gene ontology (GO) analysis has been categorized at varied graph levels based on three primary ontology studies viz. cellular components, biological processes, and molecular functions. The KEGG analysis of the assembled transcripts using KAAS showed presumed components of metabolic pathways which perhaps implicated in diverse metabolic pathways responsible for salt tolerance viz. glycolysis/gluconeogenesis, oxidative phosphorylation, glutathione metabolism, etc. further involving in salt tolerance. Overall, 90 salt stress tolerant genes were identified as of 186 salt-related transcripts. Several genes have been found executing normally in the TCA cycle pathway, integral membrane proteins, generation of the osmoprotectants, enzymatic pathway associated with salt tolerance. Recognized genes fit diverse groups of salt stress genes viz. abc transporter, betaine, sodium antiporter, sodium symporter, trehalose, ectoine, and choline, that belong to different families of genes involved in the pathway of salt stress. The control sample of the bacterium showed elevated high proportion of transcript contigs (29%) while upto 20% salt stress treated sample of the bacterium showed a higher percentage of transcript contigs (31.28%). A total of 1,288 and 1,133 transcript contigs were measured entirely as novel transcript contigs in both control and treated samples, respectively. The structure and function of 10 significant salt stress tolerant genes of Staphylococcus sp. have been analyzed in this study. The information acquired in the present study possibly used to recognize and clone the salt stress tolerant genes and support in developing the salt stress-tolerant plant varieties to expand the agricultural productivity in the saline system.

Approaches to study CRISPR RNA biogenesis and the key players involved

Clustered regularly interspaced short palindromic repeats (CRISPR) and CRISPR-associated (Cas) proteins provide an inheritable and adaptive immune system against phages and foreign genetic elements in many bacteria and archaea. The three stages of CRISPR-Cas immunity comprise adaptation, CRISPR RNA (crRNA) biogenesis and interference. The maturation of the pre-crRNA into mature crRNAs, short guide RNAs that target invading nucleic acids, is crucial for the functionality of CRISPR-Cas defense systems. Mature crRNAs assemble with Cas proteins into the ribonucleoprotein (RNP) effector complex and guide the Cas nucleases to the cognate foreign DNA or RNA target. Experimental approaches to characterize these crRNAs, the specific steps toward their maturation and the involved factors, include RNA-seq analyses, enzyme assays, methods such as cryo-electron microscopy, the crystallization of proteins, or UV-induced protein-RNA crosslinking coupled to mass spectrometry analysis. Complex and multiple interactions exist between CRISPR-cas-encoded specific riboendonucleases such as Cas6, Cas5d and Csf5, endonucleases with dual functions in maturation and interference such as the enzymes of the Cas12 and Cas13 families, and nucleases belonging to the cell's degradosome such as RNase E, PNPase and RNase J, both in the maturation as well as in interference. The results of these studies have yielded a picture of unprecedented diversity of sequences, enzymes and biochemical mechanisms.

Methylation guide RNA evolution in archaea: structure, function and genomic organization of 110 C/D box sRNA families across six Pyrobaculum species

Archaeal homologs of eukaryotic C/D box small nucleolar RNAs (C/D box sRNAs) guide precise 2′-O-methyl modification of ribosomal and transfer RNAs. Although C/D box sRNA genes constitute one of the largest RNA gene families in archaeal thermophiles, most genomes have incomplete sRNA gene annotation because reliable, fully automated detection methods are not available. We expanded and curated a comprehensive gene set across six species of the crenarchaeal genus Pyrobaculum, particularly rich in C/D box sRNA genes. Using high-throughput small RNA sequencing, specialized computational searches and comparative genomics, we analyzed 526 Pyrobaculum C/D box sRNAs, organizing them into 110 families based on synteny and conservation of guide sequences which determine methylation targets. We examined gene duplications and rearrangements, including one family that has expanded in a pattern similar to retrotransposed repetitive elements in eukaryotes. New training data and inclusion of kink-turn secondary structural features enabled creation of an improved search model. Our analyses provide the most comprehensive, dynamic view of C/D box sRNA evolutionary history within a genus, in terms of modification function, feature plasticity, and gene mobility.

The protein regulator ArgR and the sRNA derived from the 3'-UTR region of its gene, ArgX, both regulate the arginine deiminase pathway in Lactococcus lactis

Small regulatory RNAs (sRNAs) and their enormous potential and versatility have provided us with an astounding insight in the complexity of bacterial transcriptomes. sRNAs have been shown to be involved in a variety of cellular processes that range from stress to general metabolism. Here we report that the gene encoding the transcriptional regulator ArgR is immediately followed by the gene of the small regulatory RNA ArgX. The latter is transcribed from its own promoter. The production of ArgX is induced by increasing arginine concentrations and repressed by CcpA. Previously, ArgR was shown to act as a transcriptional repressor of the catabolic arginine deiminase pathway (arc operon) by binding in the promoter region of arcA. Here we demonstrate that ArgX downregulates arc mRNA levels. Furthermore, ArgX putatively blocks the translation of one of the genes in the operon, arcC1, a process that would redirect an intermediate in arginine degradation, carbamoyl phosphate, towards pyrimidine synthesis. Our findings exemplify, for the first time, the combinatorial power of a transcription factor and a small regulatory RNA derived from the 3’-UTR region. The regulators ArgR and ArgX share a common target, but act on transcription and on RNA level, respectively.

Expanding the Vocabulary of Peptide Signals in Streptococcus mutans

Streptococci, including the dental pathogen Streptococcus mutans, undergo cell-to-cell signaling that is mediated by small peptides to control critical physiological functions such as adaptation to the environment, control of subpopulation behaviors and regulation of virulence factors. One such model pathway is the regulation of genetic competence, controlled by the ComRS signaling system and the peptide XIP. However, recent research in the characterization of this pathway has uncovered novel operons and peptides that are intertwined into its regulation. These discoveries, such as cell lysis playing a critical role in XIP release and importance of bacterial self-sensing during the signaling process, have caused us to reevaluate previous paradigms and shift our views on the true purpose of these signaling systems. The finding of new peptides such as the ComRS inhibitor XrpA and the peptides of the RcrRPQ operon also suggests there may be more peptides hidden in the genomes of streptococci that could play critical roles in the physiology of these organisms. In this review, we summarize the recent findings in S. mutans regarding the integration of other circuits into the ComRS signaling pathway, the true mode of XIP export, and how the RcrRPQ operon controls competence activation. We also look at how new technologies can be used to re-annotate the genome to find new open reading frames that encode peptide signals. Together, this summary of research will allow us to reconsider how we perceive these systems to behave and lead us to expand our vocabulary of peptide signals within the genus Streptococcus.

Characterization of the transcriptome of Haloferax volcanii, grown under four different conditions, with mixed RNA-Seq

Haloferax volcanii is a well-established model species for haloarchaea. Small scale RNomics and bioinformatics predictions were used to identify small non-coding RNAs (sRNAs), and deletion mutants revealed that sRNAs have important regulatory functions. A recent dRNA-Seq study was used to characterize the primary transcriptome. Unexpectedly, it was revealed that, under optimal conditions, H. volcanii contains more non-coding sRNAs than protein-encoding mRNAs. However, the dRNA-Seq approach did not contain any length information. Therefore, a mixed RNA-Seq approach was used to determine transcript length and to identify additional transcripts, which are not present under optimal conditions. In total, 50 million paired end reads of 150 nt length were obtained. 1861 protein-coding RNAs (cdRNAs) were detected, which encoded 3092 proteins. This nearly doubled the coverage of cdRNAs, compared to the previous dRNA-Seq study. About 2/3 of the cdRNAs were monocistronic, and 1/3 covered more than one gene. In addition, 1635 non-coding sRNAs were identified. The highest fraction of non-coding RNAs were cis antisense RNAs (asRNAs). Analysis of the length distribution revealed that sRNAs have a median length of about 150 nt. Based on the RNA-Seq and dRNA-Seq results, genes were chosen to exemplify characteristics of the H. volcanii transcriptome by Northern blot analyses, e.g. 1) the transcript patterns of gene clusters can be straightforward, but also very complex, 2) many transcripts differ in expression level under the four analyzed conditions, 3) some genes are transcribed into RNA isoforms of different length, which can be differentially regulated, 4) transcripts with very long 5'-UTRs and with very long 3'-UTRs exist, and 5) about 30% of all cdRNAs have overlapping 3'-ends, which indicates, together with the asRNAs, that H. volcanii makes ample use of sense-antisense interactions. Taken together, this RNA-Seq study, together with a previous dRNA-Seq study, enabled an unprecedented view on the H. volcanii transcriptome.

The Primary Antisense Transcriptome of Halobacterium salinarum NRC-1

Antisense RNAs (asRNAs) are present in diverse organisms and play important roles in gene regulation. In this work, we mapped the primary antisense transcriptome in the halophilic archaeon Halobacterium salinarum NRC-1. By reanalyzing publicly available data, we mapped antisense transcription start sites (aTSSs) and inferred the probable 3′ ends of these transcripts. We analyzed the resulting asRNAs according to the size, location, function of genes on the opposite strand, expression levels and conservation. We show that at least 21% of the genes contain asRNAs in H. salinarum. Most of these asRNAs are expressed at low levels. They are located antisense to genes related to distinctive characteristics of H. salinarum, such as bacteriorhodopsin, gas vesicles, transposases and other important biological processes such as translation. We provide evidence to support asRNAs in type II toxin–antitoxin systems in archaea. We also analyzed public Ribosome profiling (Ribo-seq) data and found that ~10% of the asRNAs are ribosome-associated non-coding RNAs (rancRNAs), with asRNAs from transposases overrepresented. Using a comparative transcriptomics approach, we found that ~19% of the asRNAs annotated in H. salinarum belong to genes with an ortholog in Haloferax volcanii, in which an aTSS could be identified with positional equivalence. This shows that most asRNAs are not conserved between these halophilic archaea.

De novo whole transcriptome profiling of Edwardsiella tarda isolated from infected fish (Labeo catla)

Edwardsiella tarda belongs to the genera of Gram negative bacterium mainly associated with edwardsiellosis, the most commonly found infectious fish disease throughout the globe. E. tarda is also a widespread pathogen which cause infections such as cellulitis or gas gangrene and generalized infections in humans. To control the escalating infection of E. trada on various species, it is essential to decoded the mysterious mechanism behind the bacterial infection at transcript level. In this present study, we carry out a de novo E. tarda Whole transcriptome sequencing, isolated from infected fish intestine using SOLiD sequencing platform. RNA-Seq data analysis was performed using various bioinformatics pipelines. Protein-protein interaction study for pathway enrichment and gene ontology study were executed for further investigation. Assembly statistics for E. tarda dataset showed that the number of transcript contigs was 9657 out of which 6749 were GO annotated whereas 1528 were not assigned any GO terms. GO analysis showed that the expressed genes were enhanced with molecular function, cellular component and biological process. A KEGG enrichment study showed that pathway's that are directly linked with immune diseases like Rheumatoid arthritis (0.2%), Tuberculosis (0.3%) Endocytosis (0.6%) was considerably enriched. Protein-protein interaction study showed that most of the expressed proteins were involved in metabolic pathways, flagellar assembly, Propanoate metabolism, Microbial metabolism in diverse environments, Butanoate metabolism and Carbon. The present study provides novel E. tarda transcriptome sequence data, allowing us to identify biologically significant genes and their functional relationship with fish diseases, and will be useful in recognize the reliable therapeutic targets in near feature.

An integrative overview of genomic, transcriptomic and proteomic analyses in organohalide respiration research

Organohalide respiration (OHR) is a crucial process in the global halogen cycle and of interest for bioremediation. However, investigations on OHR are hampered by the restricted genetic accessibility and the poor growth yields of many organohalide-respiring bacteria (OHRB). Therefore, genomics, transcriptomics and proteomics are often used to investigate OHRB. In general, these gene expression studies are more useful when the data of the different 'omics' approaches are integrated and compared among a wide range of cultivation conditions and ideally involve several closely related OHRB. Despite the availability of a couple of proteomic and transcriptomic datasets dealing with OHRB, such approaches are currently not covered in reviews. Therefore, we here present an integrative and comparative overview of omics studies performed with the OHRB Sulfurospirillum multivorans, Dehalococcoides mccartyi, Desulfitobacterium spp. and Dehalobacter restrictus. Genes, transcripts, proteins and the regulatory and biochemical processes involved in OHR are discussed, and a comprehensive view on the unusual metabolism of D. mccartyi, which is one of the few bacteria possibly using a quinone-independent respiratory chain, is provided. Several 'omics'-derived theories on OHRB, e.g. the organohalide-respiratory chain, hydrogen metabolism, corrinoid biosynthesis or one-carbon metabolism are critically discussed on the basis of this integrative approach.

Analysis of post-transcriptional RNA metabolism in prokaryotes

Post-transcriptional RNA metabolic pathways play important roles in permitting prokaryotes to operate under a variety of environmental conditions. Although significant progress has been made during the last decade in deciphering RNA processing pathways in a number of bacteria, a complete understanding of post-transcriptional RNA metabolism in any single microorganism is far from reality. Here we describe multiple experimental approaches that can be used to study mRNA stability, tRNA and rRNA processing, sRNA metabolism, and polyadenylation in prokaryotes. The methods described here can be readily utilized in both Gram-negative and Gram-positive bacteria with simple modifications.

Genome-scale analysis of Acetobacterium bakii reveals the cold adaptation of psychrotolerant acetogens by post-transcriptional regulation

Acetogens synthesize acetyl-CoA via CO2 or CO fixation, producing organic compounds. Despite their ecological and industrial importance, their transcriptional and post-transcriptional regulation has not been systematically studied. With completion of the genome sequence of Acetobacterium bakii (4.28-Mb), we measured changes in the transcriptome of this psychrotolerant acetogen in response to temperature variations under autotrophic and heterotrophic growth conditions. Unexpectedly, acetogenesis genes were highly up-regulated at low temperatures under heterotrophic, as well as autotrophic, growth conditions. To mechanistically understand the transcriptional regulation of acetogenesis genes via changes in RNA secondary structures of 5'-untranslated regions (5'-UTR), the primary transcriptome was experimentally determined, and 1379 transcription start sites (TSS) and 1100 5'-UTR were found. Interestingly, acetogenesis genes contained longer 5'-UTR with lower RNA-folding free energy than other genes, revealing that the 5'-UTRs control the RNA abundance of the acetogenesis genes under low temperature conditions. Our findings suggest that post-transcriptional regulation via RNA conformational changes of 5'-UTRs is necessary for cold-adaptive acetogenesis.

Regulatory non-coding sRNAs in bacterial metabolic pathway engineering

Non-coding RNAs (ncRNAs) are versatile and powerful controllers of gene expression that have been increasingly linked to cellular metabolism and phenotype. In bacteria, identified and characterized ncRNAs range from trans-acting, multi-target small non-coding RNAs to dynamic, cis-encoded regulatory untranslated regions and riboswitches. These native regulators have inspired the design and construction of many synthetic RNA devices. In this work, we review the design, characterization, and impact of ncRNAs in engineering both native and exogenous metabolic pathways in bacteria. We also consider the opportunities afforded by recent high-throughput approaches for characterizing sRNA regulators and their corresponding networks to showcase their potential applications and impact in engineering bacterial metabolism.

Capping-RACE: a simple, accurate, and sensitive 5' RACE method for use in prokaryotes

Rapid amplification of cDNA ends (RACE) is a prevalent technique used to obtain the 5' ends of transcripts. Several different 5' RACE methods have been developed, and one particularly simple and efficient approach called CapFinder relies on the 5' cap-dependent template-switching that occurs in eukaryotes. However, most prokaryotic transcripts lack a 5' cap structure. Here, we report a procedure to capture primary transcripts based on capping the 5' triphosphorylated RNA in prokaryotes. Primary transcripts were first treated with vaccinia capping enzyme to add a 5' cap structure. First-strand cDNA was then synthesized using Moloney murine leukaemia virus reverse transcriptase. Finally, a template-switching oligonucleotide with a tail containing three ribonucleic acid guanines was hybridized to the cDNA 3' poly(C) and further used as template for reverse transcriptase. It is oligonucleotide sequence independent and is more sensitive compared to RLM-RACE. This approach specifically identified the transcription start sites of ompA, sodB and shiA in Escherichia coli and of ompA, rne and rppH in Brucella melitensis. Furthermore, we also successfully identified the transcription start sites of small noncoding genes ryhB and micC in E. coli and bsnc135 and bsnc149 in B. melitensis. Our findings suggest that Capping-RACE is a simple, accurate, and sensitive 5' RACE method for use in prokaryotes.

Small Noncoding Regulatory RNAs from Pseudomonas aeruginosa and Burkholderia cepacia Complex

Cystic fibrosis (CF) is the most life-limiting autosomal recessive disorder in Caucasians. CF is characterized by abnormal viscous secretions that impair the function of several tissues, with chronic bacterial airway infections representing the major cause of early decease of these patients. Pseudomonas aeruginosa and bacteria from the Burkholderia cepacia complex (Bcc) are the leading pathogens of CF patients’ airways. A wide array of virulence factors is responsible for the success of infections caused by these bacteria, which have tightly regulated responses to the host environment. Small noncoding RNAs (sRNAs) are major regulatory molecules in these bacteria. Several approaches have been developed to study P. aeruginosa sRNAs, many of which were characterized as being involved in the virulence. On the other hand, the knowledge on Bcc sRNAs remains far behind. The purpose of this review is to update the knowledge on characterized sRNAs involved in P. aeruginosa virulence, as well as to compile data so far achieved on sRNAs from the Bcc and their possible roles on bacteria virulence.

A Small Regulatory RNA Controls Cell Wall Biosynthesis and Antibiotic Resistance

Small regulatory RNAs play an important role in the adaptation to changing conditions. Here, we describe a differentially expressed small regulatory RNA (sRNA) that affects various cellular processes in the plant pathogen Agrobacterium tumefaciens Using a combination of bioinformatic predictions and comparative proteomics, we identified nine targets, most of which are positively regulated by the sRNA. According to these targets, we named the sRNA PmaR for peptidoglycan biosynthesis, motility, and ampicillin resistance regulator. Agrobacterium spp. are long known to be naturally resistant to high ampicillin concentrations, and we can now explain this phenotype by the positive PmaR-mediated regulation of the beta-lactamase gene ampC Structure probing revealed a spoon-like structure of the sRNA, with a single-stranded loop that is engaged in target interaction in vivo and in vitro Several riboregulators have been implicated in antibiotic resistance mechanisms, such as uptake and efflux transporters, but PmaR represents the first example of an sRNA that directly controls the expression of an antibiotic resistance gene.IMPORTANCE The alphaproteobacterium Agrobacterium tumefaciens is able to infect various eudicots causing crown gall tumor formation. Based on its unique ability of interkingdom gene transfer, Agrobacterium serves as a crucial biotechnological tool for genetic manipulation of plant cells. The presence of hundreds of putative sRNAs in this organism suggests a considerable impact of riboregulation on A. tumefaciens physiology. Here, we characterized the biological function of the sRNA PmaR that controls various processes crucial for growth, motility, and virulence. Among the genes directly targeted by PmaR is ampC coding for a beta-lactamase that confers ampicillin resistance, suggesting that the sRNA is crucial for fitness in the competitive microbial composition of the rhizosphere.

Discovering RNA-Based Regulatory Systems for Yersinia Virulence

The genus Yersinia includes three human pathogenic species, Yersinia pestis, the causative agent of the bubonic and pneumonic plague, and enteric pathogens Y. enterocolitica and Y. pseudotuberculosis that cause a number of gut-associated diseases. Over the past years a large repertoire of RNA-based regulatory systems has been discovered in these pathogens using different RNA-seq based approaches. Among them are several conserved or species-specific RNA-binding proteins, regulatory and sensory RNAs as well as various RNA-degrading enzymes. Many of them were shown to control the expression of important virulence-relevant factors and have a very strong impact on Yersinia virulence. The precise targets, the molecular mechanism and their role for Yersinia pathogenicity is only known for a small subset of identified genus- or species-specific RNA-based control elements. However, the ongoing development of new RNA-seq based methods and data analysis methods to investigate the synthesis, composition, translation, decay, and modification of RNAs in the bacterial cell will help us to generate a more comprehensive view of Yersinia RNA biology in the near future.

SuhB, a small non-coding RNA involved in catabolite repression of tetralin degradation genes in Sphingopyxis granuli strain TFA

Global dRNA-seq analysis of transcription start sites combined with in silico annotation using Infernal software revealed the expression of 91 putative non-coding sRNA in Sphingopyxis granuli TFA cells grown on different carbon sources. Excluding housekeeping sRNAs, only one additional sRNA, which belongs to the Rfam SuhB family (RF00519), was detected by Infernal but with an incorrect size according to the experimental results. SuhB is highly conserved across the Sphingopyxis genus. Expression data revealed that SuhB is present in rapidly growing TFA cells. A suhB deletion mutant exhibited de-repression of tetralin degradation (thn) gene expression and higher amounts of their LysR-type activator, ThnR, under conditions of carbon catabolite repression (CCR). Interaction between SuhB and the 5'UTR of thnR mRNA was demonstrated in vitro. Moreover, co-immunoprecipitation experiments, combined with fluorescence measurements of gfp fusions to the 5'UTR of thnR mRNA and the phenotype of an hfq deletion mutant, suggest the involvement of Hfq in this interaction. Taken together, these data support an Hfq-mediated repressive role for SuhB, on ThnR mRNA translation that prevents thn gene induction. SuhB, which is a highly conserved sRNA in the Sphingopyxis genus, is the first identified element directly involved in CCR of thn gene expression in S. granuli strain TFA.

Neisseria gonorrhoeae Exposed to Sublethal Levels of Hydrogen Peroxide Mounts a Complex Transcriptional Response

Neisseria gonorrhoeae mounts a substantial transcriptional program in response to hydrogen peroxide (HP), a prominent reactive oxygen species (ROS) encountered during infection. We tested which strain FA1090 genes show differential transcript abundance in response to sublethal amounts of HP to differentiate HP-responsive signaling from widespread cellular death and dysregulation. RNA sequencing (RNA-Seq) revealed that 150 genes were significantly upregulated and 143 genes downregulated following HP exposure. We annotated HP-responsive operons and all transcriptional start sites (TSSs) and identified which TSSs responded to HP treatment. We compared the HP responses and other previously reported genes and found only partial overlapping of other regulatory networks, indicating that the response to HP involves multiple biological functions. Using a representative subset of responsive genes, we validated the RNA-Seq results and found that the HP transcriptome was similar to that of sublethal organic peroxide. None of the genes in the representative subset, however, responded to sublethal levels of HOCl or O2 -. These results support the idea that N. gonorrhoeae may use variations in HP levels as a signal for different stages of infection. IMPORTANCE The strict human pathogen Neisseria gonorrhoeae is the only causative agent of the sexually transmitted disease gonorrhea. This bacterium encounters hydrogen peroxide produced from host cells during infection, but the organism survives in the presence of this antimicrobial agent. This work shows that the bacterium responds to hydrogen peroxide by regulating the expression of many genes involved in multiple processes.

ANNOgesic: a Swiss army knife for the RNA-seq based annotation of bacterial/archaeal genomes

To understand the gene regulation of an organism of interest, a comprehensive genome annotation is essential. While some features, such as coding sequences, can be computationally predicted with high accuracy based purely on the genomic sequence, others, such as promoter elements or noncoding RNAs, are harder to detect. RNA sequencing (RNA-seq) has proven to be an efficient method to identify these genomic features and to improve genome annotations. However, processing and integrating RNA-seq data in order to generate high-resolution annotations is challenging, time consuming, and requires numerous steps. We have constructed a powerful and modular tool called ANNOgesic that provides the required analyses and simplifies RNA-seq-based bacterial and archaeal genome annotation. It can integrate data from conventional RNA-seq and differential RNA-seq and predicts and annotates numerous features, including small noncoding RNAs, with high precision. The software is available under an open source license (ISCL) at https://pypi.org/project/ANNOgesic/.

Modeling leaderless transcription and atypical genes results in more accurate gene prediction in prokaryotes

In a conventional view of the prokaryotic genome organization, promoters precede operons and ribosome binding sites (RBSs) with Shine-Dalgarno consensus precede genes. However, recent experimental research suggesting a more diverse view motivated us to develop an algorithm with improved gene-finding accuracy. We describe GeneMarkS-2, an ab initio algorithm that uses a model derived by self-training for finding species-specific (native) genes, along with an array of precomputed "heuristic" models designed to identify harder-to-detect genes (likely horizontally transferred). Importantly, we designed GeneMarkS-2 to identify several types of distinct sequence patterns (signals) involved in gene expression control, among them the patterns characteristic for leaderless transcription as well as noncanonical RBS patterns. To assess the accuracy of GeneMarkS-2, we used genes validated by COG (Clusters of Orthologous Groups) annotation, proteomics experiments, and N-terminal protein sequencing. We observed that GeneMarkS-2 performed better on average in all accuracy measures when compared with the current state-of-the-art gene prediction tools. Furthermore, the screening of ∼5000 representative prokaryotic genomes made by GeneMarkS-2 predicted frequent leaderless transcription in both archaea and bacteria. We also observed that the RBS sites in some species with leadered transcription did not necessarily exhibit the Shine-Dalgarno consensus. The modeling of different types of sequence motifs regulating gene expression prompted a division of prokaryotic genomes into five categories with distinct sequence patterns around the gene starts.

Small RNAs Involved in Regulation of Nitrogen Metabolism

Global (metabolic) regulatory networks allow microorganisms to survive periods of nitrogen starvation or general nutrient stress. Uptake and utilization of various nitrogen sources are thus commonly tightly regulated in Prokarya (Bacteria and Archaea) in response to available nitrogen sources. Those well-studied regulations occur mainly at the transcriptional and posttranslational level. Surprisingly, and in contrast to their involvement in most other stress responses, small RNAs (sRNAs) involved in the response to environmental nitrogen fluctuations are only rarely reported. In addition to sRNAs indirectly affecting nitrogen metabolism, only recently it was demonstrated that three sRNAs were directly involved in regulation of nitrogen metabolism in response to changes in available nitrogen sources. All three trans-acting sRNAs are under direct transcriptional control of global nitrogen regulators and affect expression of components of nitrogen metabolism (glutamine synthetase, nitrogenase, and PII-like proteins) by either masking the ribosome binding site and thus inhibiting translation initiation or stabilizing the respective target mRNAs. Most likely, there are many more sRNAs and other types of noncoding RNAs, e.g., riboswitches, involved in the regulation of nitrogen metabolism in Prokarya that remain to be uncovered. The present review summarizes the current knowledge on sRNAs involved in nitrogen metabolism and their biological functions and targets.

Functional Genomics Approaches to Studying Symbioses between Legumes and Nitrogen-Fixing Rhizobia

Biological nitrogen fixation gives legumes a pronounced growth advantage in nitrogen-deprived soils and is of considerable ecological and economic interest. In exchange for reduced atmospheric nitrogen, typically given to the plant in the form of amides or ureides, the legume provides nitrogen-fixing rhizobia with nutrients and highly specialised root structures called nodules. To elucidate the molecular basis underlying physiological adaptations on a genome-wide scale, functional genomics approaches, such as transcriptomics, proteomics, and metabolomics, have been used. This review presents an overview of the different functional genomics approaches that have been performed on rhizobial symbiosis, with a focus on studies investigating the molecular mechanisms used by the bacterial partner to interact with the legume. While rhizobia belonging to the alpha-proteobacterial group (alpha-rhizobia) have been well studied, few studies to date have investigated this process in beta-proteobacteria (beta-rhizobia).

Epitranscriptomics: RNA Modifications in Bacteria and Archaea

The increasingly complex functionality of RNA is contrasted by its simple chemical composition. RNA is generally built from only four different nucleotides (adenine, guanine, cytosine, and uracil). To date, >160 chemical modifications are known to decorate RNA molecules and thereby alter their function or stability. Many RNA modifications are conserved throughout bacteria, archaea, and eukaryotes, while some are unique to each branch of life. Most known modifications occur at internal positions, while there is limited diversity at the termini. The dynamic nature of RNA modifications and newly discovered regulatory functions of some of these RNA modifications gave birth to a new field, now often referred to as "epitranscriptomics." This review highlights the major developments in this field and summarizes detection principles for internal as well as 5'-terminal mRNA modifications in prokaryotes and archaea to investigate their biological significance.

Post-transcriptional control of virulence gene expression in Staphylococcus aureus

Opportunistic pathogens have to be ready to change life-style whenever the occasion arises, and therefore need to keep tight control over the expression of their virulence factors. Doubly so for commensal bacteria, such as Staphylococcus aureus, which should avoid harming their hosts when they are in a state of peaceful co-existence. S. aureus carries very few sigma factors to help define the transcriptional programs, but instead uses a plethora of small RNA molecules and RNA-RNA interactions to regulate gene expression post-transcriptionally. The endoribonucleases RNase III and RNase Y contribute to this regulatory diversity, and provide a link to RNA-decay and intra-cellular spatiotemporal control of expression. In this review we describe some of these post-transcriptional mechanisms as well as some of the novel transcriptomic approaches that have been used to find and to study them.

The evolutionary impact of intragenic FliA promoters in proteobacteria

In Escherichia coli, one sigma factor recognizes the majority of promoters, and six 'alternative' sigma factors recognize specific subsets of promoters. The alternative sigma factor FliA (σ28 ) recognizes promoters upstream of many flagellar genes. We previously showed that most E. coli FliA binding sites are located inside genes. However, it was unclear whether these intragenic binding sites represent active promoters. Here, we construct and assay transcriptional promoter-lacZ fusions for all 52 putative FliA promoters previously identified by ChIP-seq. These experiments, coupled with integrative analysis of published genome-scale transcriptional datasets, strongly suggest that most intragenic FliA binding sites are active promoters that transcribe highly unstable RNAs. Additionally, we show that widespread intragenic FliA-dependent transcription may be a conserved phenomenon, but that specific promoters are not themselves conserved. We conclude that intragenic FliA-dependent promoters and the resulting RNAs are unlikely to have important regulatory functions. Nonetheless, one intragenic FliA promoter is broadly conserved and constrains evolution of the overlapping protein-coding gene. Thus, our data indicate that intragenic regulatory elements can influence bacterial protein evolution and suggest that the impact of intragenic regulatory sequences on genome evolution should be considered more broadly.

The Non-Coding Regulatory RNA Revolution in Archaea

Small non-coding RNAs (sRNAs) are ubiquitously found in the three domains of life playing large-scale roles in gene regulation, transposable element silencing and defense against foreign elements. While a substantial body of experimental work has been done to uncover function of sRNAs in Bacteria and Eukarya, the functional roles of sRNAs in Archaea are still poorly understood. Recently, high throughput studies using RNA-sequencing revealed that sRNAs are broadly expressed in the Archaea, comprising thousands of transcripts within the transcriptome during non-challenged and stressed conditions. Antisense sRNAs, which overlap a portion of a gene on the opposite strand (cis-acting), are the most abundantly expressed non-coding RNAs and they can be classified based on their binding patterns to mRNAs (3' untranslated region (UTR), 5' UTR, CDS-binding). These antisense sRNAs target many genes and pathways, suggesting extensive roles in gene regulation. Intergenic sRNAs are less abundantly expressed and their targets are difficult to find because of a lack of complete overlap between sRNAs and target mRNAs (trans-acting). While many sRNAs have been validated experimentally, a regulatory role has only been reported for very few of them. Further work is needed to elucidate sRNA-RNA binding mechanisms, the molecular determinants of sRNA-mediated regulation, whether protein components are involved and how sRNAs integrate with complex regulatory networks.

Bacterial RNA Biology on a Genome Scale

Bacteria are an exceedingly diverse group of organisms whose molecular exploration is experiencing a renaissance. While the classical view of bacterial gene expression was relatively simple, the emerging view is more complex, encompassing extensive post-transcriptional control involving riboswitches, RNA thermometers, and regulatory small RNAs (sRNAs) associated with the RNA-binding proteins CsrA, Hfq, and ProQ, as well as CRISPR/Cas systems that are programmed by RNAs. Moreover, increasing interest in members of the human microbiota and environmental microbial communities has highlighted the importance of understudied bacterial species with largely unknown transcriptome structures and RNA-based control mechanisms. Collectively, this creates a need for global RNA biology approaches that can rapidly and comprehensively analyze the RNA composition of a bacterium of interest. We review such approaches with a focus on RNA-seq as a versatile tool to investigate the different layers of gene expression in which RNA is made, processed, regulated, modified, translated, and turned over.

The Small RNA ncS35 Regulates Growth in Burkholderia cenocepacia J2315

Small RNAs play an important role in the survival of bacteria in diverse environments. We explored the physiological role of ncS35, a small RNA expressed in B. cenocepacia J2315, an opportunistic pathogen in cystic fibrosis patients. In cystic fibrosis patients, infections can lead to “cepacia syndrome,” a rapidly progressing and often fatal pneumonia. Infections with Burkholderia spp. are difficult to threat with antibiotics because of their high intrinsic resistance and ability to form biofilms. We show that ncS35 attenuates the growth and reduces the metabolic rate of B. cenocepacia and influences biofilm structure. This demonstrates that as-yet-uncharacterized small RNAs with regulatory function can influence physiological traits of B. cenocepacia that are relevant for infection.

Burkholderia cenocepacia J2315 is a member of the B. cepacia complex. It has a large genome with three replicons and one plasmid; 7,261 genes code for annotated proteins, while 113 code for functional RNAs. Small regulatory RNAs of B. cenocepacia have not yet been functionally characterized. We investigated a small regulatory RNA, designated ncS35, that was discovered by differential RNA sequencing. Its expression under various conditions was quantified, and a deletion mutant, ΔncS35, was constructed. Compared to planktonic growth in a rich medium, the expression of ncS35 was elevated when B. cenocepacia J2315 was grown in biofilms and in minimal medium. Cells of the deletion mutant showed increased aggregation, higher metabolic activity, a higher growth rate, and an increased susceptibility to tobramycin. A transcriptomic analysis revealed upregulation of the phenylacetic acid and tryptophan degradation pathways in ΔncS35. Computational target prediction indicated that ncS35 likely interacts with the first gene of the tryptophan degradation pathway. Overall, we demonstrated that small RNA ncS35 is a noncoding RNA with an attenuating effect on the metabolic rate and growth. It is possible that slower growth protects B. cenocepacia J2315 against stressors acting on fast-dividing cells and enhances survival under unfavorable conditions.

IMPORTANCE Small RNAs play an important role in the survival of bacteria in diverse environments. We explored the physiological role of ncS35, a small RNA expressed in B. cenocepacia J2315, an opportunistic pathogen in cystic fibrosis patients. In cystic fibrosis patients, infections can lead to “cepacia syndrome,” a rapidly progressing and often fatal pneumonia. Infections with Burkholderia spp. are difficult to threat with antibiotics because of their high intrinsic resistance and ability to form biofilms. We show that ncS35 attenuates the growth and reduces the metabolic rate of B. cenocepacia and influences biofilm structure. This demonstrates that as-yet-uncharacterized small RNAs with regulatory function can influence physiological traits of B. cenocepacia that are relevant for infection.

The 2011 German Enterohemorrhagic Escherichia Coli O104:H4 Outbreak-The Danger Is Still Out There

Enterohemorrhagic Escherichia coli (EHEC) are Shiga toxin (Stx) producing bacteria causing a disease characterized by bloody (or non-bloody) diarrhea, which might progress to hemolytic uremic syndrome (HUS). EHEC O104:H4 caused the largest ever recorded EHEC outbreak in Germany in 2011, which in addition showed the so far highest incidence rate of EHEC-related HUS worldwide. The aggressive outbreak strain carries an unusual combination of virulence traits characteristic to both EHEC-a chromosomally integrated Stx-encoding bacteriophage, and enteroaggregative Escherichia coli-pAA plasmid-encoded aggregative adherence fimbriae mediating its tight adhesion to epithelia cells. There are currently still open questions regarding the 2011 EHEC outbreak, e.g., with respect to the exact molecular mechanisms resulting in the hypervirulence of the strain, the natural reservoir of EHEC O104:H4, and suitable therapeutic strategies. Nevertheless, our knowledge on these issues has substantially expanded since 2011. Here, we present an overview of the epidemiological, clinical, microbiological, and molecular biological data available on the 2011 German EHEC O104:H4 outbreak.

RNA Sequencing and Analysis in Microorganisms for Metabolic Network Reconstruction

There is a strict interplay between metabolic networks and transcriptional regulation in bacteria; indeed, the transcriptome regulation, affecting the expression of large gene sets, can be used to predict the likely "on" or "off" state of metabolic genes as a function of environmental factors. Up to date, many bacterial transcriptomes have been studied by RNAseq, hundreds of experiments have been performed, and Giga bases of sequences have been produced. All this transcriptional information could potentially be integrated into metabolic networks in order to obtain a more comprehensive view of their regulation and to increase their prediction power.To get high-quality transcriptomic data, to be integrated into metabolic networks, it is paramount to clearly know how to produce highly informative RNA sequencing libraries and how to manage RNA sequencing data.In this chapter, we will get across the main steps of an RNAseq experiment: from removal of ribosomal RNAs, to strand-specific library preparation, till data analysis and integration. We will try to share our experience and know-how, to give you a precise protocol to follow, and some useful recommendations or tips and tricks to adopt in order to go straightforward toward a successful RNAseq experiment.

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.

Identification of small RNAs abundant in Burkholderia cenocepacia biofilms reveal putative regulators with a potential role in carbon and iron metabolism

Small RNAs play a regulatory role in many central metabolic processes of bacteria, as well as in developmental processes such as biofilm formation. Small RNAs of Burkholderia cenocepacia, an opportunistic pathogenic beta-proteobacterium, are to date not well characterised. To address that, we performed genome-wide transcriptome structure analysis of biofilm grown B. cenocepacia J2315. 41 unannotated short transcripts were identified in intergenic regions of the B. cenocepacia genome. 15 of these short transcripts, highly abundant in biofilms, widely conserved in Burkholderia sp. and without known function, were selected for in-depth analysis. Expression profiling showed that most of these sRNAs are more abundant in biofilms than in planktonic cultures. Many are also highly abundant in cells grown in minimal media, suggesting they are involved in adaptation to nutrient limitation and growth arrest. Their computationally predicted targets include a high proportion of genes involved in carbon metabolism. Expression and target genes of one sRNA suggest a potential role in regulating iron homoeostasis. The strategy used for this study to detect sRNAs expressed in B. cenocepacia biofilms has successfully identified sRNAs with a regulatory function.

Determination of the Genome and Primary Transcriptome of Syngas Fermenting Eubacterium limosum ATCC 8486

Autotrophic conversion of CO₂ to value-added biochemicals has received considerable attention as a sustainable route to replace fossil fuels. Particularly, anaerobic acetogenic bacteria are naturally capable of reducing CO₂ or CO to various metabolites. To fully utilize their biosynthetic potential, an understanding of acetogenesis-related genes and their regulatory elements is required. Here, we completed the genome sequence of the syngas fermenting Eubacterium limosum ATCC 8486 and determined its transcription start sites (TSS). We constructed a 4.4 Mb long circular genome with a GC content of 47.2% and 4,090 protein encoding genes. To understand the transcriptional and translational regulation, the primary transcriptome was augmented, identifying 1,458 TSSs containing a high pyrimidine (T/C) and purine nucleotide (A/G) content at the −1 and +1 position, respectively, along with 1,253 5′-untranslated regions, and principal promoter elements such as −10 (TATAAT) and −35 (TTGACA), and Shine-Dalgarno motifs (GGAGR). Further analysis revealed 93 non-coding RNAs, including one for potential transcriptional regulation of the hydrogenase complex via interaction with molybdenum or tungsten cofactors, which in turn controls formate dehydrogenase activity of the initial step of Wood-Ljungdahl pathway. Our results provide comprehensive genomic information for strain engineering to enhance the syngas fermenting capacity of acetogenic bacteria.

The primary transcriptome of Neisseria meningitidis and its interaction with the RNA chaperone Hfq

Neisseria meningitidis is a human commensal that can also cause life-threatening meningitis and septicemia. Despite growing evidence for RNA-based regulation in meningococci, their transcriptome structure and output of regulatory small RNAs (sRNAs) are incompletely understood. Using dRNA-seq, we have mapped at single-nucleotide resolution the primary transcriptome of N. meningitidis strain 8013. Annotation of 1625 transcriptional start sites defines transcription units for most protein-coding genes but also reveals a paucity of classical σ70-type promoters, suggesting the existence of activators that compensate for the lack of −35 consensus sequences in N. meningitidis. The transcriptome maps also reveal 65 candidate sRNAs, a third of which were validated by northern blot analysis. Immunoprecipitation with the RNA chaperone Hfq drafts an unexpectedly large post-transcriptional regulatory network in this organism, comprising 23 sRNAs and hundreds of potential mRNA targets. Based on this data, using a newly developed gfp reporter system we validate an Hfq-dependent mRNA repression of the putative colonization factor PrpB by the two trans-acting sRNAs RcoF1/2. Our genome-wide RNA compendium will allow for a better understanding of meningococcal transcriptome organization and riboregulation with implications for colonization of the human nasopharynx.

Systematic analysis of transcription start sites in avian development

Cap Analysis of Gene Expression (CAGE) in combination with single-molecule sequencing technology allows precision mapping of transcription start sites (TSSs) and genome-wide capture of promoter activities in differentiated and steady state cell populations. Much less is known about whether TSS profiling can characterize diverse and non-steady state cell populations, such as the approximately 400 transitory and heterogeneous cell types that arise during ontogeny of vertebrate animals. To gain such insight, we used the chick model and performed CAGE-based TSS analysis on embryonic samples covering the full 3-week developmental period. In total, 31,863 robust TSS peaks (>1 tag per million [TPM]) were mapped to the latest chicken genome assembly, of which 34% to 46% were active in any given developmental stage. ZENBU, a web-based, open-source platform, was used for interactive data exploration. TSSs of genes critical for lineage differentiation could be precisely mapped and their activities tracked throughout development, suggesting that non-steady state and heterogeneous cell populations are amenable to CAGE-based transcriptional analysis. Our study also uncovered a large set of extremely stable housekeeping TSSs and many novel stage-specific ones. We furthermore demonstrated that TSS mapping could expedite motif-based promoter analysis for regulatory modules associated with stage-specific and housekeeping genes. Finally, using Brachyury as an example, we provide evidence that precise TSS mapping in combination with Clustered Regularly Interspaced Short Palindromic Repeat (CRISPR)-on technology enables us, for the first time, to efficiently target endogenous avian genes for transcriptional activation. Taken together, our results represent the first report of genome-wide TSS mapping in birds and the first systematic developmental TSS analysis in any amniote species (birds and mammals). By facilitating promoter-based molecular analysis and genetic manipulation, our work also underscores the value of avian models in unravelling the complex regulatory mechanism of cell lineage specification during amniote development.

Plasticity first: molecular signatures of a complex morphological trait in filamentous cyanobacteria

Background

Filamentous cyanobacteria that differentiate multiple cell types are considered the peak of prokaryotic complexity and their evolution has been studied in the context of multicellularity origins. Species that form true-branching filaments exemplify the most complex cyanobacteria. However, the mechanisms underlying the true-branching morphology remain poorly understood despite of several investigations that focused on the identification of novel genes or pathways. An alternative route for the evolution of novel traits is based on existing phenotypic plasticity. According to that scenario – termed genetic assimilation – the fixation of a novel phenotype precedes the fixation of the genotype.

Results

Here we show that the evolution of transcriptional regulatory elements constitutes a major mechanism for the evolution of new traits. We found that supplementation with sucrose reconstitutes the ancestral branchless phenotype of two true-branching Fischerella species and compared the transcription start sites (TSSs) between the two phenotypic states. Our analysis uncovers several orthologous TSSs whose transcription level is correlated with the true-branching phenotype. These TSSs are found in genes that encode components of the septosome and elongasome (e.g., fraC and mreB).

Conclusions

The concept of genetic assimilation supplies a tenable explanation for the evolution of novel traits but testing its feasibility is hindered by the inability to recreate and study the evolution of present-day traits. We present a novel approach to examine transcription data for the plasticity first route and provide evidence for its occurrence during the evolution of complex colony morphology in true-branching cyanobacteria. Our results reveal a route for evolution of the true-branching phenotype in cyanobacteria via modification of the transcription level of pre-existing genes. Our study supplies evidence for the ‘plasticity-first’ hypothesis and highlights the importance of transcriptional regulation in the evolution of novel traits.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-1053-5) contains supplementary material, which is available to authorized users.

A combinatorial strategy of alternative promoter use during differentiation of a heterocystous cyanobacterium

Heterocystous cyanobacteria such as Nostoc sp. are filamentous photosynthetic organisms that, in response to nitrogen deficiency, undergo a differentiation process transforming certain, semi-regularly spaced cells into heterocysts, devoted to nitrogen fixation. During transition to a nitrogen-fixing regime, growth of most vegetative cells in the filament is temporarily arrested due to nutritional deprivation, but developing heterocysts require intense transcriptional activity. Therefore, the coexistence of arrested vegetative cells and actively developing prospective heterocysts relies on the simultaneous operation of somewhat opposite transcriptional programs. We have identified genes with multiple nitrogen-responsive transcriptional starts appearing in seemingly paradoxical combinations. For instance, sigA, encoding the RNA polymerase housekeeping sigma factor, is transcribed from one major nitrogen stress-repressed promoter and from a second, nitrogen stress-induced promoter. Here, we show that both promoters are expressed with complementary temporal dynamics. Using a gfp reporter we also show that transcription from the inducible promoter takes place exclusively in differentiating heterocysts and is already detected before any morphological or fluorescence signature of differentiation is observed. Tandem promoters with opposite dynamics could operate a compensatory mechanism in which repression of transcription from the major promoter operative in vegetative cells is offset by transcription from a new promoter only in developing heterocyst.