Bacterial transcriptomics: what is beyond the RNA horiz-ome?

Exploring the transcription start sites and other genomic features facilitates the accurate identification and annotation of small RNAs across multiple stress conditions in Mycobacterium tuberculosis

Mycobacterium tuberculosis (MTB) is a pathogen that is known for its ability to persist in harsh environments and cause chronic infections. Understanding the regulatory networks of MTB is crucial for developing effective treatments. Small regulatory RNAs (sRNAs) play important roles in gene expression regulation in all kingdoms of life, and their classification based solely on genomic location can be imprecise due to the computational-based prediction of protein-coding genes in bacteria, which often neglects segments of mRNA such as 5'UTRs, 3'UTRs, and intercistronic regions of operons. To address this issue, our study simultaneously discovered genomic features such as TSSs, UTRs, and operons together with sRNAs in the M. tuberculosis H37Rv strain (ATCC 27294) across multiple stress conditions. Our analysis identified 1,376 sRNA candidates and 8,173 TSSs in MTB, providing valuable insights into its complex regulatory landscape. TSS mapping enabled us to classify these sRNAs into more specific categories, including promoter-associated sRNAs, 5'UTR-derived sRNAs, 3'UTR-derived sRNAs, true intergenic sRNAs, and antisense sRNAs. Three of these sRNA candidates were experimentally validated using 3'-RACE-PCR: predictedRNA_0240, predictedRNA_0325, and predictedRNA_0578. Future characterization and validation are necessary to fully elucidate the functions and roles of these sRNAs in MTB. Our study is the first to simultaneously unravel TSSs and sRNAs in MTB and demonstrate that the identification of other genomic features, such as TSSs, UTRs, and operons, allows for more accurate and specific classification of sRNAs.

Decoding the molecular concerto: Toxicotranscriptomic evaluation of microplastic and nanoplastic impacts on aquatic organisms

The pervasive and steadily increasing presence of microplastics/nanoplastics (MPs/NPs) in aquatic environments has raised significant concerns regarding their potential adverse effects on aquatic organisms and their integration into trophic dynamics. This emerging issue has garnered the attention of (eco)toxicologists, promoting the utilization of toxicotranscriptomics to unravel the responses of aquatic organisms not only to MPs/NPs but also to a wide spectrum of environmental pollutants. This review aims to systematically explore the broad repertoire of predicted molecular responses by aquatic organisms, providing valuable intuitions into complex interactions between plastic pollutants and aquatic biota. By synthesizing the latest literature, present analysis sheds light on transcriptomic signatures like gene expression, interconnected pathways and overall molecular mechanisms influenced by various plasticizers. Harmful effects of these contaminants on key genes/protein transcripts associated with crucial pathways lead to abnormal immune response, metabolic response, neural response, apoptosis and DNA damage, growth, development, reproductive abnormalities, detoxification, and oxidative stress in aquatic organisms. However, unique challenge lies in enhancing the fingerprint of MPs/NPs, presenting complicated enigma that requires decoding their specific impact at molecular levels. The exploration endeavors, not only to consolidate existing knowledge, but also to identify critical gaps in understanding, push forward the frontiers of knowledge about transcriptomic signatures of plastic contaminants. Moreover, this appraisal emphasizes the imperative to monitor and mitigate the contamination of commercially important aquatic species by MPs/NPs, highlighting the pivotal role that regulatory frameworks must play in protecting all aquatic ecosystems. This commitment aligns with the broader goal of ensuring the sustainability of aquatic resources and the resilience of ecosystems facing the growing threat of plastic pollutants.

Heterogeneous lineage-specific arginine deiminase expression within dental microbiome species

Arginine catabolism by the bacterial arginine deiminase system (ADS) has anticariogenic properties through the production of ammonia, which modulates the pH of the oral environment. Given the potential protective capacity of the ADS pathway, the exploitation of ADS-competent oral microbes through pre- or probiotic applications is a promising therapeutic target to prevent tooth decay. To date, most investigations of the ADS in the oral cavity and its relation to caries have focused on indirect measures of activity or on specific bacterial groups, yet the pervasiveness and rate of expression of the ADS operon in diverse mixed microbial communities in oral health and disease remain an open question. Here, we use a multivariate approach, combining ultra-deep metatranscriptomic sequencing with paired metataxonomic and in vitro citrulline quantification to characterize the microbial community and ADS operon expression in healthy and late-stage cavitated teeth. While ADS activity is higher in healthy teeth, we identify multiple bacterial lineages with upregulated ADS activity on cavitated teeth that are distinct from those found on healthy teeth using both reference-based mapping and de novo assembly methods. Our dual metataxonomic and metatranscriptomic approach demonstrates the importance of species abundance for gene expression data interpretation and that patterns of differential expression can be skewed by low-abundance groups. Finally, we identify several potential candidate probiotic bacterial lineages within species that may be useful therapeutic targets for the prevention of tooth decay and propose that the development of a strain-specific, mixed-microbial probiotic may be a beneficial approach given the heterogeneity of taxa identified here across health groups.

IMPORTANCE

Tooth decay is the most common preventable chronic disease, affecting more than two billion people globally. The development of caries on teeth is primarily a consequence of acid production by cariogenic bacteria that inhabit the plaque microbiome. Other bacterial strains in the oral cavity may suppress or prevent tooth decay by producing ammonia as a byproduct of the arginine deiminase metabolic pathway, increasing the pH of the plaque biofilm. While the benefits of arginine metabolism on oral health have been extensively documented in specific bacterial groups, the prevalence and consistency of arginine deiminase system (ADS) activity among oral bacteria in a community context remain an open question. In the current study, we use a multi-omics approach to document the pervasiveness of the expression of the ADS operon in both health and disease to better understand the conditions in which ADS activity may prevent tooth decay.

Machine learning identification of Pseudomonas aeruginosa strains from colony image data

When grown on agar surfaces, microbes can produce distinct multicellular spatial structures called colonies, which contain characteristic sizes, shapes, edges, textures, and degrees of opacity and color. For over one hundred years, researchers have used these morphology cues to classify bacteria and guide more targeted treatment of pathogens. Advances in genome sequencing technology have revolutionized our ability to classify bacterial isolates and while genomic methods are in the ascendancy, morphological characterization of bacterial species has made a resurgence due to increased computing capacities and widespread application of machine learning tools. In this paper, we revisit the topic of colony morphotype on the within-species scale and apply concepts from image processing, computer vision, and deep learning to a dataset of 69 environmental and clinical Pseudomonas aeruginosa strains. We find that colony morphology and complexity under common laboratory conditions is a robust, repeatable phenotype on the level of individual strains, and therefore forms a potential basis for strain classification. We then use a deep convolutional neural network approach with a combination of data augmentation and transfer learning to overcome the typical data starvation problem in biological applications of deep learning. Using a train/validation/test split, our results achieve an average validation accuracy of 92.9% and an average test accuracy of 90.7% for the classification of individual strains. These results indicate that bacterial strains have characteristic visual 'fingerprints' that can serve as the basis of classification on a sub-species level. Our work illustrates the potential of image-based classification of bacterial pathogens and highlights the potential to use similar approaches to predict medically relevant strain characteristics like antibiotic resistance and virulence from colony data.

Antisense transcription and its roles in adaption to environmental stress in E. coli

It has been reported that a highly varying proportion (1% ∼ 93%) of genes in various prokaryotes have antisense RNA (asRNA) transcription. However, the extent of the pervasiveness of asRNA transcription in the well-studied E. coli K12 strain has thus far been an issue of debate. Furthermore, very little is known about the expression patterns and functions of asRNAs under various conditions. To fill these gaps, we determined the transcriptomes and proteomes of E. coli K12 at multiple time points in five culture conditions using strand-specific RNA-seq, differential RNA-seq, and quantitative mass spectrometry methods. To reduce artifacts of possible transcriptional noise, we identified asRNA using stringent criteria with biological replicate verification and transcription start sites (TSSs) information included. We identified a total of 660 asRNAs, which were generally short and largely condition-dependently transcribed. We found that the proportions of the genes which had asRNA transcription highly depended on the culture conditions and time points. We classified the transcriptional activities of the genes in six transcriptional modes according to their relative levels of asRNA to mRNA. Many genes changed their transcriptional modes at different time points of the culture conditions, and such transitions can be described in a well-defined manner. Intriguingly, the protein levels and mRNA levels of genes in the sense-only/sense-dominant mode were moderately correlated, but the same was not true for genes in the balanced/antisense-dominant mode, in which asRNAs were at a comparable or higher level to mRNAs. These observations were further validated by western blot on candidate genes, where an increase in asRNA transcription diminished gene expression in one case and enhanced it in another. These results suggest that asRNAs may directly or indirectly regulate translation by forming duplexes with cognate mRNAs. Thus, asRNAs may play an important role in the bacterium's responses to environmental changes during growth and adaption to different environments.

IMPORTANCE

The cis -antisense RNA (asRNA) is a type of understudied RNA molecules in prokaryotes, which is believed to be important in regulating gene expression. Our current understanding of asRNA is constrained by inconsistent reports about its identification and properties. These discrepancies are partially caused by a lack of sufficient samples, biological replicates, and culture conditions. This study aimed to overcome these disadvantages and identified 660 putative asRNAs using integrated information from strand-specific RNA-seq, differential RNA-seq, and mass spectrometry methods. In addition, we explored the relative expression between asRNAs and sense RNAs and investigated asRNA regulated transcriptional activity changes over different culture conditions and time points. Our work strongly suggests that asRNAs may play a crucial role in bacterium's responses to environmental changes during growth and adaption to different environments.

Development of ONT-cappable-seq to unravel the transcriptional landscape of Pseudomonas phages

RNA sequencing has become the method of choice to study the transcriptional landscape of phage-infected bacteria. However, short-read RNA sequencing approaches generally fail to capture the primary 5' and 3' boundaries of transcripts, confounding the discovery of key transcription initiation and termination events as well as operon architectures. Yet, the elucidation of these elements is crucial for the understanding of the strategy of transcription regulation during the infection process, which is currently lacking beyond a handful of model phages. We developed ONT-cappable-seq, a specialized long-read RNA sequencing technique that allows end-to-end sequencing of primary prokaryotic transcripts using the Nanopore sequencing platform. We applied ONT-cappable-seq to study transcription of Pseudomonas aeruginosa phage LUZ7, obtaining a comprehensive genome-wide map of viral transcription start sites, terminators, and complex operon structures that fine-regulate gene expression. Our work provides new insights in the RNA biology of a non-model phage, unveiling distinct promoter architectures, putative small non-coding viral RNAs, and the prominent regulatory role of terminators during infection. The robust workflow presented here offers a framework to obtain a global, yet fine-grained view of phage transcription and paves the way for standardized, in-depth transcription studies for microbial viruses or bacteria in general.

Visualization of mRNA Expression in Pseudomonas aeruginosa Aggregates Reveals Spatial Patterns of Fermentative and Denitrifying Metabolism

Gaining insight into the behavior of bacteria at the single-cell level is important given that heterogeneous microenvironments strongly influence microbial physiology. The hybridization chain reaction (HCR) is a technique that provides in situ molecular signal amplification, enabling simultaneous mapping of multiple target RNAs at small spatial scales. To refine this method for biofilm applications, we designed and validated new probes to visualize the expression of key catabolic genes in Pseudomonas aeruginosa aggregates. In addition to using existing probes for the dissimilatory nitrate reductase (narG), we developed probes for a terminal oxidase (ccoN1), nitrite reductase (nirS), nitrous oxide reductase (nosZ), and acetate kinase (ackA). These probes can be used to determine gene expression levels across heterogeneous populations such as biofilms. Using these probes, we quantified gene expression across oxygen gradients in aggregate populations grown using the agar block biofilm assay (ABBA). We observed distinct patterns of catabolic gene expression, with upregulation occurring in particular ABBA regions both within individual aggregates and over the aggregate population. Aerobic respiration (ccoN1) showed peak expression under oxic conditions, whereas fermentation (ackA) showed peak expression in the anoxic cores of high metabolic activity aggregates near the air-agar interface. Denitrification genes narG, nirS, and nosZ showed peak expression in hypoxic and anoxic regions, although nirS expression remained at peak levels deeper into anoxic environments than other denitrification genes. These results reveal that the microenvironment correlates with catabolic gene expression in aggregates, and they demonstrate the utility of HCR in unveiling cellular activities at the microscale level in heterogeneous populations. IMPORTANCE To understand bacteria in diverse contexts, we must understand the variations in behaviors and metabolisms they express spatiotemporally. Populations of bacteria are known to be heterogeneous, but the ways this variation manifests can be challenging to characterize due to technical limitations. By focusing on energy conservation, we demonstrate that HCR v3.0 can visualize nuances in gene expression, allowing us to understand how metabolism in Pseudomonas aeruginosa biofilms responds to microenvironmental variation at high spatial resolution. We validated probes for four catabolic genes, including a constitutively expressed oxidase, acetate kinase, nitrite reductase, and nitrous oxide reductase. We showed that the genes for different modes of metabolism are expressed in overlapping but distinct subpopulations according to oxygen concentrations in a predictable fashion. The spatial transcriptomic technique described here has the potential to be used to map microbial activities across diverse environments.

Overlapping genes in natural and engineered genomes

Modern genome-scale methods that identify new genes, such as proteogenomics and ribosome profiling, have revealed, to the surprise of many, that overlap in genes, open reading frames and even coding sequences is widespread and functionally integrated into prokaryotic, eukaryotic and viral genomes. In parallel, the constraints that overlapping regions place on genome sequences and their evolution can be harnessed in bioengineering to build more robust synthetic strains and constructs. With a focus on overlapping protein-coding and RNA-coding genes, this Review examines their discovery, topology and biogenesis in the context of their genome biology. We highlight exciting new uses for sequence overlap to control translation, compress synthetic genetic constructs, and protect against mutation.

The "guiding" principles of noncoding RNA function

The human genome is pervasively transcribed and yet only a small fraction of these RNAs (less than 2%) are known to code for proteins. The vast majority of the RNAs are classified as noncoding RNAs (ncRNAs) and are further subgrouped as small (shorter than 200 bases) and long noncoding RNAs. The ncRNAs have been identified in all three domains of life and regulate diverse cellular processes through transcriptional and posttranscriptional gene regulation. Most of these RNAs work in conjunction with proteins forming a wide array of base pairing interactions. The determinants of these base pairing interactions are now becoming more evident and show striking similarities among the diverse group of ncRNAs. Here we present a mechanistic overview of pairing between RNA-RNA or RNA-DNA that dictates the function of ncRNAs; we provide examples to illustrate that ncRNAs work through shared evolutionary mechanisms that encompasses a guide-target interaction, involving not only classical Watson-Crick but also noncanonical Wobble and Hoogsteen base pairing. We also highlight the similarities in target selection, proofreading, and the ruler mechanism of ncRNA-protein complexes that confers target specificity and target site selection. This article is categorized under: Regulatory RNAs/RNAi/Riboswitches > Regulatory RNAs RNA-Based Catalysis > RNA-Mediated Cleavage RNA Evolution and Genomics > RNA and Ribonucleoprotein Evolution.

Global Analysis of the Zinc Homeostasis Network in Pseudomonas aeruginosa and Its Gene Expression Dynamics

Zinc is one of the most important trace elements for life and its deficiency, like its excess, can be fatal. In the bacterial opportunistic pathogen Pseudomonas aeruginosa, Zn homeostasis is not only required for survival, but also for virulence and antibiotic resistance. Thus, the bacterium possesses multiple Zn import/export/storage systems. In this work, we determine the expression dynamics of the entire P. aeruginosa Zn homeostasis network at both transcript and protein levels. Precisely, we followed the switch from a Zn-deficient environment, mimicking the initial immune strategy to counteract bacterial infections, to a Zn-rich environment, representing the phagocyte metal boost used to eliminate an engulfed pathogen. Thanks to the use of the NanoString technology, we timed the global silencing of Zn import systems and the orchestrated induction of Zn export systems. We show that the induction of Zn export systems is hierarchically organized as a function of their impact on Zn homeostasis. Moreover, we identify PA2807 as a novel Zn resistance component in P. aeruginosa and highlight new regulatory links among Zn-homeostasis systems. Altogether, this work unveils a sophisticated and adaptive homeostasis network, which complexity is key in determining a pathogen spread in the environment and during host-colonization.

Heterotrophic bacterial diazotrophs are more abundant than their cyanobacterial counterparts in metagenomes covering most of the sunlit ocean

Biological nitrogen fixation contributes significantly to marine primary productivity. The current view depicts few cyanobacterial diazotrophs as the main marine nitrogen fixers. Here, we used 891 Tara Oceans metagenomes derived from surface waters of five oceans and two seas to generate a manually curated genomic database corresponding to free-living, filamentous, colony-forming, particle-attached, and symbiotic bacterial and archaeal populations. The database provides the genomic content of eight cyanobacterial diazotrophs including a newly discovered population related to known heterocystous symbionts of diatoms, as well as 40 heterotrophic bacterial diazotrophs that considerably expand the known diversity of abundant marine nitrogen fixers. These 48 populations encapsulate 92% of metagenomic signal for known nifH genes in the sunlit ocean, suggesting that the genomic characterization of the most abundant marine diazotrophs may be nearing completion. Newly identified heterotrophic bacterial diazotrophs are widespread, express their nifH genes in situ, and also occur in large planktonic size fractions where they might form aggregates that provide the low-oxygen microenvironments required for nitrogen fixation. Critically, we found heterotrophic bacterial diazotrophs to be more abundant than cyanobacterial diazotrophs in most metagenomes from the open oceans and seas, emphasizing the importance of a wide range of heterotrophic populations in the marine nitrogen balance.

Global distribution patterns of marine nitrogen-fixers by imaging and molecular methods

Nitrogen fixation has a critical role in marine primary production, yet our understanding of marine nitrogen-fixers (diazotrophs) is hindered by limited observations. Here, we report a quantitative image analysis pipeline combined with mapping of molecular markers for mining >2,000,000 images and >1300 metagenomes from surface, deep chlorophyll maximum and mesopelagic seawater samples across 6 size fractions (<0.2–2000 μm). We use this approach to characterise the diversity, abundance, biovolume and distribution of symbiotic, colony-forming and particle-associated diazotrophs at a global scale. We show that imaging and PCR-free molecular data are congruent. Sequence reads indicate diazotrophs are detected from the ultrasmall bacterioplankton (<0.2 μm) to mesoplankton (180–2000 μm) communities, while images predict numerous symbiotic and colony-forming diazotrophs (>20 µm). Using imaging and molecular data, we estimate that polyploidy can substantially affect gene abundances of symbiotic versus colony-forming diazotrophs. Our results support the canonical view that larger diazotrophs (>10 μm) dominate the tropical belts, while unicellular cyanobacterial and non-cyanobacterial diazotrophs are globally distributed in surface and mesopelagic layers. We describe co-occurring diazotrophic lineages of different lifestyles and identify high-density regions of diazotrophs in the global ocean. Overall, we provide an update of marine diazotroph biogeographical diversity and present a new bioimaging-bioinformatic workflow.

Nitrogen fixation by diazotrophs is critical for marine primary production. Using Tara Oceans datasets, this study combines a quantitative image analysis pipeline with metagenomic mining to provide an improved global overview of diazotroph abundance, diversity and distribution.

Elucidating the Regulatory Elements for Transcription Termination and Posttranscriptional Processing in the Streptomyces clavuligerus Genome

Identification of transcriptional regulatory elements in the GC-rich Streptomyces genome is essential for the production of novel biochemicals from secondary metabolite biosynthetic gene clusters (smBGCs). Despite many efforts to understand the regulation of transcription initiation in smBGCs, information on the regulation of transcription termination and posttranscriptional processing remains scarce. In this study, we identified the transcriptional regulatory elements in β-lactam antibiotic-producing Streptomyces clavuligerus ATCC 27064 by determining a total of 1,427 transcript 3'-end positions (TEPs) using the term-seq method. Termination of transcription was governed by three classes of TEPs, of which each displayed unique sequence features. The data integration with transcription start sites and transcriptome data generated 1,648 transcription units (TUs) and 610 transcription unit clusters (TUCs). TU architecture showed that the transcript abundance in TU isoforms of a TUC was potentially affected by the sequence context of their TEPs, suggesting that the regulatory elements of TEPs could control the transcription level in additional layers. We also identified TU features of a xenobiotic response element (XRE) family regulator and DUF397 domain-containing protein, particularly showing the abundance of bidirectional TEPs. Finally, we found that 189 noncoding TUs contained potential cis- and trans-regulatory elements that played a major role in regulating the 5' and 3' UTR. These findings highlight the role of transcriptional regulatory elements in transcription termination and posttranscriptional processing in Streptomyces sp.IMPORTANCE Streptomyces sp. is a great source of bioactive secondary metabolites, including antibiotics, antifungal agents, antiparasitic agents, immunosuppressant compounds, and other drugs. Secondary metabolites are synthesized via multistep conversions of the precursor molecules from primary metabolism, governed by multicomplex enzymes from secondary metabolite biosynthetic gene clusters. As their production is closely related with the growth phase and dynamic cellular status in response to various intra- and extracellular signals, complex regulatory systems tightly control the gene expressions related to secondary metabolism. In this study, we determined genome-wide transcript 3'-end positions and transcription units in the β-lactam antibiotic producer Streptomyces clavuligerus ATCC 27064 to elucidate the transcriptional regulatory elements in transcription termination and posttranscriptional processing by integration of multiomics data. These unique features, such as transcript 3'-end sequence, potential riboregulators, and potential 3'-untranslated region (UTR) cis-regulatory elements, can be potentially used to design engineering tools that can regulate the transcript abundance of genes for enhancing secondary metabolite production.

Protein dosage of the lldPRD operon is correlated with RNase E-dependent mRNA processing

The ability of Escherichia coli to grow on L-lactate as a sole carbon source depends on the expression of the lldPRD operon. A striking feature of this operon is that the transcriptional regulator (LldR) encoding gene is located between the permease (LldP) and the dehydrogenase (LldD) encoding genes. In this study we report that dosage of the LldP, LldR, and LldD proteins is not modulated on the transcriptional level. Instead, modulation of protein dosage is primarily correlated with RNase E-dependent mRNA processing events that take place within the lldR mRNA, leading to the immediate inactivation of lldR, to differential segmental stabilities of the resulting cleavage products, and to differences in the translation efficiencies of the three cistrons. A model for the processing events controlling the molar quantities of the proteins in the lldPRD operon is presented and discussed.ImportanceAdjustment of gene expression is critical for proper cell function. For the case of polycistronic transcripts, posttranscriptional regulatory mechanisms can be used to fine-tune the expression of individual cistrons. Here, we elucidate how protein dosage of the Escherichia coli lldPRD operon, which presents the paradox of having the gene encoding a regulator protein located between genes that code for a permease and an enzyme, is regulated. Our results demonstrate that the key event in this regulatory mechanism involves the RNase E-dependent cleavage of the primary lldPRD transcript at internal site(s) located within the lldR cistron, resulting in a drastic decrease of intact lldR mRNA, to differential segmental stabilities of the resulting cleavage products, and to differences in the translation efficiencies of the three cistrons.

Shedding Light on the Ghost Proteome

Conventionally, eukaryotic mRNAs were thought to be monocistronic, leading to the translation of a single protein. However, large-scale proteomics has led to the identification of proteins translated from alternative open reading frames (AltORFs) in mRNAs. AltORFs are found in addition to predicted reference ORFs and noncoding RNA. Alternative proteins are not represented in the conventional protein databases, and this 'Ghost proteome' was not considered until recently. Some of these proteins are functional, and there is growing evidence that they are involved in central functions in physiological and physiopathological contexts. Here, we review how this Ghost proteome fills the gap in our understanding of signaling pathways, establishes new markers of pathologies, and highlights therapeutic targets.

Correlated chromosomal periodicities according to the growth rate and gene expression

Linking genetic information to population fitness is crucial to understanding living organisms. Despite the abundant knowledge of the genetic contribution to growth, the overall patterns/features connecting genes, their expression, and growth remain unclear. To reveal the quantitative and direct connections, systematic growth assays of single-gene knockout Escherichia coli strains under both rich and poor nutritional conditions were performed; subsequently, the resultant growth rates were associated with the original expression levels of the knockout genes in the parental genome. Comparative analysis of growth and the transcriptome identified not only the nutritionally differentiated fitness cost genes but also a significant correlation between the growth rates of the single-gene knockout strains and the original expression levels of these knockout genes in the parental strain, regardless of the nutritional variation. In addition, the coordinated chromosomal periodicities of the wild-type transcriptome and the growth rates of the strains lacking the corresponding genes were observed. The common six-period periodicity was somehow attributed to the essential genes, although the underlying mechanism remains to be addressed. The correlated chromosomal periodicities associated with the gene expression-growth dataset were highly valuable for bacterial growth prediction and discovering the working principles governing minimal genetic information.

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.

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.

Trends in Symbiont-Induced Host Cellular Differentiation

Bacteria participate in a wide diversity of symbiotic associations with eukaryotic hosts that require precise interactions for bacterial recognition and persistence. Most commonly, host-associated bacteria interfere with host gene expression to modulate the immune response to the infection. However, many of these bacteria also interfere with host cellular differentiation pathways to create a hospitable niche, resulting in the formation of novel cell types, tissues, and organs. In both of these situations, bacterial symbionts must interact with eukaryotic regulatory pathways. Here, we detail what is known about how bacterial symbionts, from pathogens to mutualists, control host cellular differentiation across the central dogma, from epigenetic chromatin modifications, to transcription and mRNA processing, to translation and protein modifications. We identify four main trends from this survey. First, mechanisms for controlling host gene expression appear to evolve from symbionts co-opting cross-talk between host signaling pathways. Second, symbiont regulatory capacity is constrained by the processes that drive reductive genome evolution in host-associated bacteria. Third, the regulatory mechanisms symbionts exhibit correlate with the cost/benefit nature of the association. And, fourth, symbiont mechanisms for interacting with host genetic regulatory elements are not bound by native bacterial capabilities. Using this knowledge, we explore how the ubiquitous intracellular Wolbachia symbiont of arthropods and nematodes may modulate host cellular differentiation to manipulate host reproduction. Our survey of the literature on how infection alters gene expression in Wolbachia and its hosts revealed that, despite their intermediate-sized genomes, different strains appear capable of a wide diversity of regulatory manipulations. Given this and Wolbachia's diversity of phenotypes and eukaryotic-like proteins, we expect that many symbiont-induced host differentiation mechanisms will be discovered in this system.

Large-scale and significant expression from pseudogenes in Sodalis glossinidius - a facultative bacterial endosymbiont

The majority of bacterial genomes have high coding efficiencies, but there are some genomes of intracellular bacteria that have low gene density. The genome of the endosymbiont Sodalis glossinidius contains almost 50 % pseudogenes containing mutations that putatively silence them at the genomic level. We have applied multiple 'omic' strategies, combining Illumina and Pacific Biosciences Single-Molecule Real-Time DNA sequencing and annotation, stranded RNA sequencing and proteome analysis to better understand the transcriptional and translational landscape of Sodalis pseudogenes, and potential mechanisms for their control. Between 53 and 74 % of the Sodalis transcriptome remains active in cell-free culture. The mean sense transcription from coding domain sequences (CDSs) is four times greater than that from pseudogenes. Comparative genomic analysis of six Illumina-sequenced Sodalis isolates from different host Glossina species shows pseudogenes make up ~40 % of the 2729 genes in the core genome, suggesting that they are stable and/or that Sodalis is a recent introduction across the genus Glossina as a facultative symbiont. These data shed further light on the importance of transcriptional and translational control in deciphering host-microbe interactions. The combination of genomics, transcriptomics and proteomics gives a multidimensional perspective for studying prokaryotic genomes with a view to elucidating evolutionary adaptation to novel environmental niches.

Full-length RNA profiling reveals pervasive bidirectional transcription terminators in bacteria

The ability to determine full-length nucleotide composition of individual RNA molecules is essential for understanding the architecture and function of a transcriptome. However, experimental approaches capable of capturing the sequences of both 5' and 3' termini of the same transcript remain scarce. In the present study, simultaneous 5' and 3' end sequencing (SEnd-seq)-a high-throughput and unbiased method that simultaneously maps transcription start and termination sites with single-nucleotide resolution-is presented. Using this method, a comprehensive view of the Escherichia coli transcriptome was obtained, which displays an unexpected level of complexity. SEnd-seq notably expands the catalogue of transcription start sites and termination sites, defines unique transcription units and detects prevalent antisense RNA. Strikingly, the results of the present study unveil widespread overlapping bidirectional terminators located between opposing gene pairs. Furthermore, it has been shown that convergent transcription is a major contributor to highly efficient bidirectional termination both in vitro and in vivo. This finding highlights an underappreciated role of RNA polymerase conflicts in shaping transcript boundaries and suggests an evolutionary strategy for modulating transcriptional output by arranging gene orientation.

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.

A crash course in sequencing for a microbiologist

For the last 40 years, "Sanger sequencing" allowed to unveil crucial secrets of life. However, this method of sequencing has been time-consuming, laborious and remains expensive even today. Human Genome Project was a huge impulse to improve sequencing technologies, and unprecedented financial and human effort prompted the development of cheaper high-throughput technologies and strategies called next-generation sequencing (NGS) or whole genome sequencing (WGS). This review will discuss applications of high-throughput methods to study bacteria in a much broader context than simply their genomes. The major goal of next-generation sequencing for a microbiologist is not really resolving another circular genomic sequence. NGS started its infancy from basic structural and functional genomics, to mature into the molecular taxonomy, phylogenetic and advanced comparative genomics. Today, the use of NGS expended capabilities of diagnostic microbiology and epidemiology. The use of RNA sequencing techniques allows studying in detail the complex regulatory processes in the bacterial cells. Finally, NGS is a key technique to study the organization of the bacterial life-from complex communities to single cells. The major challenge in understanding genomic and transcriptomic data lies today in combining it with other sources of global data such as proteome and metabolome, which hopefully will lead to the reconstruction of regulatory networks within bacterial cells that allow communicating with the environment (signalome and interactome) and virtual cell reconstruction.

The Toxin-Antidote Model of Cytoplasmic Incompatibility: Genetics and Evolutionary Implications

Wolbachia bacteria inhabit the cells of about half of all arthropod species, an unparalleled success stemming in large part from selfish invasive strategies. Cytoplasmic incompatibility (CI), whereby the symbiont makes itself essential to embryo viability, is the most common of these and constitutes a promising weapon against vector-borne diseases. After decades of theoretical and experimental struggle, major recent advances have been made toward a molecular understanding of this phenomenon. As pieces of the puzzle come together, from yeast and Drosophila fly transgenesis to CI diversity patterns in natural mosquito populations, it becomes clearer than ever that the CI induction and rescue stem from a toxin-antidote (TA) system. Further, the tight association of the CI genes with prophages provides clues to the possible evolutionary origin of this phenomenon and the levels of selection at play.

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.

Emergence of Bias During the Synthesis and Amplification of cDNA for scRNA-seq

The advent of single-cell omics technology has promoted our understanding of the genomic, epigenomic, and transcriptomic heterogeneity in individual cells. Compared to traditional sequencing studies using bulk cells, single-cell transcriptome technology is naturally more dynamic for in depth analysis of genomic variation resulting from cell division and is useful in unraveling the regulatory mechanisms of gene networks in many diseases. However, there are still some limitations of current single-cell RNA sequencing (scRNA-seq) protocols. Biases that arise during the RNA reverse transcription and cDNA pre-amplification steps are the most common problems and play pivotal roles in limiting the quantitative accuracy of scRNA-seq. In this review, we will describe how these biases emerge and impact scRNA-seq protocols. Moreover, we will introduce several current and convenient modified scRNA-seq methods that allow for bias to be decreased and estimated.

Extensive reshaping of bacterial operons by programmed mRNA decay

Bacterial operons synchronize the expression of multiple genes by placing them under the control of a shared promoter. It was previously shown that polycistronic transcripts can undergo differential RNA decay, leaving some genes within the polycistron more stable than others, but the extent of regulation by differential mRNA decay or its evolutionary conservation remains unknown. Here, we find that a substantial fraction of E. coli genes display non-uniform mRNA stoichiometries despite being coded from the same operon. We further show that these altered operon stoichiometries are shaped post-transcriptionally by differential mRNA decay, which is regulated by RNA structures that protect specific regions in the transcript from degradation. These protective RNA structures are generally coded within the protein-coding regions of the regulated genes and are frequently evolutionarily conserved. Furthermore, we provide evidence that differences in ribosome densities across polycistronic transcript segments, together with the conserved structural RNA elements, play a major role in the differential decay process. Our results highlight a major role for differential mRNA decay in shaping bacterial transcriptomes.

Bacteria utilize operonic transcription to synchronize the expression of multiple consecutive genes. However, this strategy lacks the ability to fine-tune the expression of specific operon members, which is often biologically important. In this report, we integrate multiple transcriptome-wide RNA-sequencing methods to show that bacteria commonly employ differential mRNA decay rates for genes residing within the same operon, generating differential transcript abundances for equally transcribed operon members, at steady state. By comparing the transcriptomes of different bacteria, we show that differential decay not only regulates the expression levels of hundreds of genes but also often evolutionarily conserved, providing support for its biological importance. By mapping the RNA termini positions at steady-state, we show that stabilized operon segments are protected from different RNases through a combination of protective RNA structures, which surprisingly, are often encoded within protein-coding regions and are evolutionarily conserved. In addition, we provide evidence that differential ribosome densities over the regulated operons guide the initial events in the differential decay mechanism. Our results highlight differential mRNA decay as a major shaping force of bacterial transcriptomes and gene regulatory programs.

Genome-wide mRNA processing in methanogenic archaea reveals post-transcriptional regulation of ribosomal protein synthesis

Unlike stable RNAs that require processing for maturation, prokaryotic cellular mRNAs generally follow an 'all-or-none' pattern. Herein, we used a 5΄ monophosphate transcript sequencing (5΄P-seq) that specifically captured the 5΄-end of processed transcripts and mapped the genome-wide RNA processing sites (PSSs) in a methanogenic archaeon. Following statistical analysis and stringent filtration, we identified 1429 PSSs, among which 23.5% and 5.4% were located in 5΄ untranslated region (uPSS) and intergenic region (iPSS), respectively. A predominant uridine downstream PSSs served as a processing signature. Remarkably, 5΄P-seq detected overrepresented uPSS and iPSS in the polycistronic operons encoding ribosomal proteins, and the majority upstream and proximal ribosome binding sites, suggesting a regulatory role of processing on translation initiation. The processed transcripts showed increased stability and translation efficiency. Particularly, processing within the tricistronic transcript of rplA-rplJ-rplL enhanced the translation of rplL, which can provide a driving force for the 1:4 stoichiometry of L10 to L12 in the ribosome. Growth-associated mRNA processing intensities were also correlated with the cellular ribosomal protein levels, thereby suggesting that mRNA processing is involved in tuning growth-dependent ribosome synthesis. In conclusion, our findings suggest that mRNA processing-mediated post-transcriptional regulation is a potential mechanism of ribosomal protein synthesis and stoichiometry.

A comprehensive benchmark of RNA-RNA interaction prediction tools for all domains of life

Motivation

The aim of this study is to assess the performance of RNA-RNA interaction prediction tools for all domains of life.

Results

Minimum free energy (MFE) and alignment methods constitute most of the current RNA interaction prediction algorithms. The MFE tools that include accessibility (i.e. RNAup, IntaRNA and RNAplex) to the final predicted binding energy have better true positive rates (TPRs) with a high positive predictive values (PPVs) in all datasets than other methods. They can also differentiate almost half of the native interactions from background. The algorithms that include effects of internal binding energies to their model and alignment methods seem to have high TPR but relatively low associated PPV compared to accessibility based methods.

Availability and Implementation

We shared our wrapper scripts and datasets at Github (github.com/UCanCompBio/RNA_Interactions_Benchmark). All parameters are documented for personal use.

Contact

sinan.umu@pg.canterbury.ac.nz.

Supplementary information

Supplementary data are available at Bioinformatics online.

Malonate degradation in Acinetobacter baylyi ADP1: operon organization and regulation by MdcR

Transcriptional regulators in the LysR or GntR families are typically encoded in the genomic neighbourhood of bacterial genes for malonate degradation. While these arrangements have been evaluated using bioinformatics methods, experimental studies demonstrating co-transcription of predicted operons were lacking. Here, transcriptional regulation was characterized for a cluster of mdc genes that enable a soil bacterium, Acinetobacter baylyi ADP1, to use malonate as a carbon source. Despite previous assumptions that the mdc-gene set forms one operon, our studies revealed distinct promoters in two different regions of a nine-gene cluster. Furthermore, a single promoter is insufficient to account for transcription of mdcR, a regulatory gene that is convergent to other mdc genes. MdcR, a LysR-type transcriptional regulator, was shown to bind specifically to a site where it can activate mdc-gene transcription. Although mdcR deletion prevented growth on malonate, a 1 nt substitution in the promoter of mdcA enabled MdcR-independent growth on this carbon source. Regulation was characterized by methods including transcriptional fusions, quantitative reverse transcription PCR, reverse transcription PCR, 5'-rapid amplification of cDNA ends and gel shift assays. Moreover, a new technique was developed for transcriptional characterization of low-copy mRNA by increasing the DNA copy number of specific chromosomal regions. MdcR was shown to respond to malonate, in the absence of its catabolism. These studies contribute to ongoing characterization of the structure and function of a set of 44 LysR-type transcriptional regulators in A. baylyi ADP1.

Deep sequencing approaches for the analysis of prokaryotic transcriptional boundaries and dynamics

The identification of the protein-coding regions of a genome is straightforward due to the universality of start and stop codons. However, the boundaries of the transcribed regions, conditional operon structures, non-coding RNAs and the dynamics of transcription, such as pausing of elongation, are non-trivial to identify, even in the comparatively simple genomes of prokaryotes. Traditional methods for the study of these areas, such as tiling arrays, are noisy, labour-intensive and lack the resolution required for densely-packed bacterial genomes. Recently, deep sequencing has become increasingly popular for the study of the transcriptome due to its lower costs, higher accuracy and single nucleotide resolution. These methods have revolutionised our understanding of prokaryotic transcriptional dynamics. Here, we review the deep sequencing and data analysis techniques that are available for the study of transcription in prokaryotes, and discuss the bioinformatic considerations of these analyses.

Marine Bacterial and Archaeal Ion-Pumping Rhodopsins: Genetic Diversity, Physiology, and Ecology

The recognition of a new family of rhodopsins in marine planktonic bacteria, proton-pumping proteorhodopsin, expanded the known phylogenetic range, environmental distribution, and sequence diversity of retinylidene photoproteins. At the time of this discovery, microbial ion-pumping rhodopsins were known solely in haloarchaea inhabiting extreme hypersaline environments. Shortly thereafter, proteorhodopsins and other light-activated energy-generating rhodopsins were recognized to be widespread among marine bacteria. The ubiquity of marine rhodopsin photosystems now challenges prior understanding of the nature and contributions of "heterotrophic" bacteria to biogeochemical carbon cycling and energy fluxes. Subsequent investigations have focused on the biophysics and biochemistry of these novel microbial rhodopsins, their distribution across the tree of life, evolutionary trajectories, and functional expression in nature. Later discoveries included the identification of proteorhodopsin genes in all three domains of life, the spectral tuning of rhodopsin variants to wavelengths prevailing in the sea, variable light-activated ion-pumping specificities among bacterial rhodopsin variants, and the widespread lateral gene transfer of biosynthetic genes for bacterial rhodopsins and their associated photopigments. Heterologous expression experiments with marine rhodopsin genes (and associated retinal chromophore genes) provided early evidence that light energy harvested by rhodopsins could be harnessed to provide biochemical energy. Importantly, some studies with native marine bacteria show that rhodopsin-containing bacteria use light to enhance growth or promote survival during starvation. We infer from the distribution of rhodopsin genes in diverse genomic contexts that different marine bacteria probably use rhodopsins to support light-dependent fitness strategies somewhere between these two extremes.

Insights into the Mechanisms of Basal Coordination of Transcription Using a Genome-Reduced Bacterium

Coordination of transcription in bacteria occurs at supra-operonic scales, but the extent, specificity, and mechanisms of such regulation are poorly understood. Here, we tackle this problem by profiling the transcriptome of the model organism Mycoplasma pneumoniae across 115 growth conditions. We identify three qualitatively different levels of co-expression corresponding to distinct relative orientations and intergenic properties of adjacent genes. We reveal that the degree of co-expression between co-directional adjacent operons, and more generally between genes, is tightly related to their capacity to be transcribed en bloc into the same mRNA. We further show that this genome-wide pervasive transcription of adjacent genes and operons is specifically repressed by DNA regions preferentially bound by RNA polymerases, by intrinsic terminators, and by large intergenic distances. Taken together, our findings suggest that the basal coordination of transcription is mediated by the physical entities and mechanical properties of the transcription process itself, and that operon-like behaviors may strongly vary from condition to condition.

Term-seq reveals abundant ribo-regulation of antibiotics resistance in bacteria

Riboswitches and attenuators are cis-regulatory RNA elements, most of which control bacterial gene expression via metabolite-mediated, premature transcription termination. We developed an unbiased experimental approach for genome-wide discovery of such ribo-regulators in bacteria. We also devised an experimental platform that quantitatively measures the in vivo activity of all such regulators in parallel and enables rapid screening for ribo-regulators that respond to metabolites of choice. Using this approach, we detected numerous antibiotic-responsive ribo-regulators that control antibiotic resistance genes in pathogens and in the human microbiome. Studying one such regulator in Listeria monocytogenes revealed an attenuation mechanism mediated by antibiotic-stalled ribosomes. Our results expose broad roles for conditional termination in regulating antibiotic resistance and provide a tool for discovering riboswitches and attenuators that respond to previously unknown ligands.

Transcriptome Sequencing Reveals Wide Expression Reprogramming of Basal and Unknown Genes in Leptospira biflexa Biofilms

In this work, we describe the first transcriptome based on RNA-seq technology focused on studying transcriptional changes associated with biofilm growth in a member of the genus Leptospira. As many pathogenic species of this genus can survive inside the host but also persist in environmental water, mostly forming biofilms, identifying the molecular basis of this capacity can impact the understanding of how leptospires are able to fulfill a complete life cycle that alternates between adaptation to the host and adaptation to hostile external environmental conditions. We identified several genes and regulatory networks that can be the kickoff for deepening understanding of the molecular mechanisms involving bacterial persistence via biofilm formation; understanding this is important for the future development of tools for controlling leptospirosis.

The genus Leptospira is composed of pathogenic and saprophytic spirochetes. Pathogenic Leptospira is the etiological agent of leptospirosis, a globally spread neglected disease. A key ecological feature of some pathogenic species is their ability to survive both within and outside the host. For most leptospires, the ability to persist outside the host is associated with biofilm formation, a most important bacterial strategy to face and overcome hostile environmental conditions. The architecture and biochemistry of leptospiral biofilms are rather well understood; however, the genetic program underpinning biofilm formation remains mostly unknown. In this work, we used the saprophyte Leptospira biflexa as a model organism to assess over- and underrepresented transcripts during the biofilm state, using transcriptome sequencing (RNA-seq) technology. Our results showed that some basal biological processes like DNA replication and cell division are downregulated in the mature biofilm. Additionally, we identified significant expression reprogramming for genes involved in motility, sugar/lipid metabolism, and iron scavenging, as well as for outer membrane-encoding genes. A careful manual annotation process allowed us to assign molecular functions to many previously uncharacterized genes that are probably involved in biofilm metabolism. We also provided evidence for the presence of small regulatory RNAs in this species. Finally, coexpression networks were reconstructed to pinpoint functionally related gene clusters that may explain how biofilm maintenance is regulated. Beyond elucidating some genetic aspects of biofilm formation, this work reveals a number of pathways whose functional dissection may impact our understanding of leptospiral biology, in particular how these organisms adapt to environmental changes.

IMPORTANCE In this work, we describe the first transcriptome based on RNA-seq technology focused on studying transcriptional changes associated with biofilm growth in a member of the genus Leptospira. As many pathogenic species of this genus can survive inside the host but also persist in environmental water, mostly forming biofilms, identifying the molecular basis of this capacity can impact the understanding of how leptospires are able to fulfill a complete life cycle that alternates between adaptation to the host and adaptation to hostile external environmental conditions. We identified several genes and regulatory networks that can be the kickoff for deepening understanding of the molecular mechanisms involving bacterial persistence via biofilm formation; understanding this is important for the future development of tools for controlling leptospirosis.

Genome-wide transcriptional responses of Alteromonas naphthalenivorans SN2 to contaminated seawater and marine tidal flat sediment

A genome-wide transcriptional analysis of Alteromonas naphthalenivorans SN2 was performed to investigate its ecophysiological behavior in contaminated tidal flats and seawater. The experimental design mimicked these habitats that either added naphthalene or pyruvate; tidal flat-naphthalene (TF-N), tidal flat-pyruvate (TF-P), seawater-naphthalene (SW-N), and seawater-pyruvate (SW-P). The transcriptional profiles clustered by habitat (TF-N/TF-P and SW-N/SW-P), rather than carbon source, suggesting that the former may exert a greater influence on genome-wide expression in strain SN2 than the latter. Metabolic mapping of cDNA reads from strain SN2 based on KEGG pathway showed that metabolic and regulatory genes associated with energy metabolism, translation, and cell motility were highly expressed in all four test conditions, probably highlighting the copiotrophic properties of strain SN2 as an opportunistic marine r-strategist. Differential gene expression analysis revealed that strain SN2 displayed specific cellular responses to environmental variables (tidal flat, seawater, naphthalene, and pyruvate) and exhibited certain ecological fitness traits -- its notable PAH degradation capability in seasonally cold tidal flat might be reflected in elevated expression of stress response and chaperone proteins, while fast growth in nitrogen-deficient and aerobic seawater probably correlated with high expression of glutamine synthetase, enzymes utilizing nitrite/nitrate, and those involved in the removal of reactive oxygen species.

Loss of Conserved Noncoding RNAs in Genomes of Bacterial Endosymbionts

The genomes of intracellular symbiotic or pathogenic bacteria, such as of Buchnera, Mycoplasma, and Rickettsia, are typically smaller compared with their free-living counterparts. Here we showed that noncoding RNA (ncRNA) families, which are conserved in free-living bacteria, frequently could not be detected by computational methods in the small genomes. Statistical tests demonstrated that their absence is not an artifact of low GC content or small deletions in these small genomes, and thus it was indicative of an independent loss of ncRNAs in different endosymbiotic lineages. By analyzing the synteny (conservation of gene order) between the reduced and nonreduced genomes, we revealed instances of protein-coding genes that were preserved in the reduced genomes but lost cis-regulatory elements. We found that the loss of cis-regulatory ncRNA sequences, which regulate the expression of cognate protein-coding genes, is characterized by the reduction of secondary structure formation propensity, GC content, and length of the corresponding genomic regions.

Rescuing discarded spectra: Full comprehensive analysis of a minimal proteome

A common problem encountered when performing large-scale MS proteome analysis is the loss of information due to the high percentage of unassigned spectra. To determine the causes behind this loss we have analyzed the proteome of one of the smallest living bacteria that can be grown axenically, Mycoplasma pneumoniae (729 ORFs). The proteome of M. pneumoniae cells, grown in defined media, was analyzed by MS. An initial search with both Mascot and a species-specific NCBInr database with common contaminants (NCBImpn), resulted in around 79% of the acquired spectra not having an assignment. The percentage of non-assigned spectra was reduced to 27% after re-analysis of the data with the PEAKS software, thereby increasing the proteome coverage of M. pneumoniae from the initial 60% to over 76%. Nonetheless, 33,413 spectra with assigned amino acid sequences could not be mapped to any NCBInr database protein sequence. Approximately, 1% of these unassigned peptides corresponded to PTMs and 4% to M. pneumoniae protein variants (deamidation and translation inaccuracies). The most abundant peptide sequence variants (Phe-Tyr and Ala-Ser) could be explained by alterations in the editing capacity of the corresponding tRNA synthases. About another 1% of the peptides not associated to any protein had repetitions of the same aromatic/hydrophobic amino acid at the N-terminus, or had Arg/Lys at the C-terminus. Thus, in a model system, we have maximized the number of assigned spectra to 73% (51,453 out of the 70,040 initial acquired spectra). All MS data have been deposited in the ProteomeXchange with identifier PXD002779 (http://proteomecentral.proteomexchange.org/dataset/PXD002779).

Changes in protein abundance are observed in bacterial isolates from a natural host

Bacterial proteomic studies frequently use strains cultured in synthetic liquid media over many generations. It is uncertain whether bacterial proteins expressed under these conditions will be the same as the repertoire found in natural environments, or when bacteria are infecting a host organism. Thus, genomic and proteomic characterization of bacteria derived from the host environment in comparison to reference strains grown in the lab, should aid understanding of pathogenesis. Isolates of Corynebacterium pseudotuberculosis were obtained from the lymph nodes of three naturally infected sheep and compared to a laboratory reference strain using bottom-up proteomics, after whole genome sequencing of each of the field isolates. These comparisons were performed following growth in liquid media that allowed us to reach the required protein amount for proteomic analysis. Over 1350 proteins were identified in the isolated strains, from which unique proteome features were revealed. Several of the identified proteins demonstrated a significant abundance difference in the field isolates compared to the reference strain even though there were no obvious differences in the DNA sequence of the corresponding gene or in nearby non-coding DNA. Higher abundance in the field isolates was observed for proteins related to hypoxia and nutrient deficiency responses as well as to thiopeptide biosynthesis.

Accelerating Discovery and Functional Analysis of Small RNAs with New Technologies

Over the past decade, bacterial small RNAs (sRNAs) have gone from a biological curiosity to being recognized as a major class of regulatory molecules. High-throughput methods for sampling the transcriptional output of bacterial cells demonstrate that sRNAs are universal features of bacterial transcriptomes, are plentiful, and appear to vary extensively over evolutionary time. With ever more bacteria coming under study, the question becomes how can we accelerate the discovery and functional characterization of sRNAs in diverse organisms. New technologies built on high-throughput sequencing are emerging that can rapidly provide global insight into the numbers and functions of sRNAs in bacteria of interest, providing information that can shape hypotheses and guide research. In this review, we describe recent developments in transcriptomics (RNA-seq) and functional genomics that we expect to help us develop an integrated, systems-level view of sRNA biology in bacteria.

Genome-wide transcription start site profiling in biofilm-grown Burkholderia cenocepacia J2315

BACKGROUND

Burkholderia cenocepacia is a soil-dwelling Gram-negative Betaproteobacterium with an important role as opportunistic pathogen in humans. Infections with B. cenocepacia are very difficult to treat due to their high intrinsic resistance to most antibiotics. Biofilm formation further adds to their antibiotic resistance. B. cenocepacia harbours a large, multi-replicon genome with a high GC-content, the reference genome of strain J2315 includes 7374 annotated genes. This study aims to annotate transcription start sites and identify novel transcripts on a whole genome scale.

METHODS

RNA extracted from B. cenocepacia J2315 biofilms was analysed by differential RNA-sequencing and the resulting dataset compared to data derived from conventional, global RNA-sequencing. Transcription start sites were annotated and further analysed according to their position relative to annotated genes.

RESULTS

Four thousand ten transcription start sites were mapped over the whole B. cenocepacia genome and the primary transcription start site of 2089 genes expressed in B. cenocepacia biofilms were defined. For 64 genes a start codon alternative to the annotated one was proposed. Substantial antisense transcription for 105 genes and two novel protein coding sequences were identified. The distribution of internal transcription start sites can be used to identify genomic islands in B. cenocepacia. A potassium pump strongly induced only under biofilm conditions was found and 15 non-coding small RNAs highly expressed in biofilms were discovered.

CONCLUSIONS

Mapping transcription start sites across the B. cenocepacia genome added relevant information to the J2315 annotation. Genes and novel regulatory RNAs putatively involved in B. cenocepacia biofilm formation were identified. These findings will help in understanding regulation of B. cenocepacia biofilm formation.

Comparative "-omics" in Mycoplasma pneumoniae Clinical Isolates Reveals Key Virulence Factors

The human respiratory tract pathogen M. pneumoniae is one of the best characterized minimal bacterium. Until now, two main groups of clinical isolates of this bacterium have been described (types 1 and 2), differing in the sequence of the P1 adhesin gene. Here, we have sequenced the genomes of 23 clinical isolates of M. pneumoniae. Studying SNPs, non-synonymous mutations, indels and genome rearrangements of these 23 strains and 4 previously sequenced ones, has revealed new subclasses in the two main groups, some of them being associated with the country of isolation. Integrative analysis of in vitro gene essentiality and mutation rates enabled the identification of several putative virulence factors and antigenic proteins; revealing recombination machinery, glycerol metabolism and peroxide production as possible factors in the genetics and physiology of these pathogenic strains. Additionally, the transcriptomes and proteomes of two representative strains, one from each of the two main groups, have been characterized to evaluate the impact of mutations on RNA and proteins levels. This study has revealed that type 2 strains show higher expression levels of CARDS toxin, a protein recently shown to be one of the major factors of inflammation. Thus, we propose that type 2 strains could be more toxigenic than type 1 strains of M. pneumoniae.

Virulence and Fitness Determinants of Uropathogenic Escherichia coli

Urinary tract infection (UTI) caused by uropathogenic Escherichia coli (UPEC) is a major global public health concern. Increasing antibiotic resistance found in clinical UPEC isolates underscores the immediate need for development of novel therapeutics against this pathogen. Better understanding of the fitness and virulence mechanisms that are integral to the pathogenesis of UTI will facilitate identification of novel strategies to prevent and treat infection with UPEC. Working towards that goal, the global UPEC research community has made great strides at unraveling various virulence and fitness genes. Here, we summarize major findings on virulence and fitness determinants that enable UPEC to successfully survive and colonize the urinary tract of mammalian hosts. Major sections of this chapter are devoted to the role of iron acquisition systems, metabolic pathways, fimbriae, flagella, toxins, biofilm formation, capsule, and strain-specific genes in the initiation and progression of UTIs. Transcriptomes of UPEC during experimental UTI in a murine model and naturally occurring UTI in women are compared to elucidate virulence mechanisms specifically involved in human UTI. Capitalizing on the advances in molecular pathogenesis research by translating these findings will help develop better clinical strategies for prevention and management of UTIs.

A Post-Genomic View of the Ecophysiology, Catabolism and Biotechnological Relevance of Sulphate-Reducing Prokaryotes

Dissimilatory sulphate reduction is the unifying and defining trait of sulphate-reducing prokaryotes (SRP). In their predominant habitats, sulphate-rich marine sediments, SRP have long been recognized to be major players in the carbon and sulphur cycles. Other, more recently appreciated, ecophysiological roles include activity in the deep biosphere, symbiotic relations, syntrophic associations, human microbiome/health and long-distance electron transfer. SRP include a high diversity of organisms, with large nutritional versatility and broad metabolic capacities, including anaerobic degradation of aromatic compounds and hydrocarbons. Elucidation of novel catabolic capacities as well as progress in the understanding of metabolic and regulatory networks, energy metabolism, evolutionary processes and adaptation to changing environmental conditions has greatly benefited from genomics, functional OMICS approaches and advances in genetic accessibility and biochemical studies. Important biotechnological roles of SRP range from (i) wastewater and off gas treatment, (ii) bioremediation of metals and hydrocarbons and (iii) bioelectrochemistry, to undesired impacts such as (iv) souring in oil reservoirs and other environments, and (v) corrosion of iron and concrete. Here we review recent advances in our understanding of SRPs focusing mainly on works published after 2000. The wealth of publications in this period, covering many diverse areas, is a testimony to the large environmental, biogeochemical and technological relevance of these organisms and how much the field has progressed in these years, although many important questions and applications remain to be explored.

Experimental approaches to identify small RNAs and their diverse roles in bacteria--what we have learnt in one decade of MicA research

Nowadays the identification of small RNAs (sRNAs) and characterization of their role within regulatory networks takes a prominent place in deciphering complex bacterial phenotypes. Compared to the study of other components of bacterial cells, this is a relatively new but fast-growing research field. Although reports on new sRNAs appear regularly, some sRNAs are already subject of research for a longer time. One of such sRNAs is MicA, a sRNA best described for its role in outer membrane remodeling, but probably having a much broader function than anticipated. An overview of what we have learnt from MicA led to the conclusion that even for this well-described sRNA, we still do not have the overall picture. More general, the story of MicA might become an experimental lead for unraveling the many sRNAs with unknown functions. In this review, three important topics in the sRNA field are covered, exemplified from the perspective of MicA: (i) identification of new sRNAs, (ii) target identification and unraveling the biological function, (iii) structural analysis. The complex mechanisms of action of MicA deliver some original insights in the sRNA field which includes the existence of dimer formation or simultaneous cis and trans regulation, and might further inspire the understanding of the function of other sRNAs.

An empirical strategy to detect bacterial transcript structure from directional RNA-seq transcriptome data

Background

As sequencing costs are being lowered continuously, RNA-seq has gradually been adopted as the first choice for comparative transcriptome studies with bacteria. Unlike microarrays, RNA-seq can directly detect cDNA derived from mRNA transcripts at a single nucleotide resolution. Not only does this allow researchers to determine the absolute expression level of genes, but it also conveys information about transcript structure. Few automatic software tools have yet been established to investigate large-scale RNA-seq data for bacterial transcript structure analysis.

Results

In this study, 54 directional RNA-seq libraries from Salmonella serovar Typhimurium (S. Typhimurium) 14028s were examined for potential relationships between read mapping patterns and transcript structure. We developed an empirical method, combined with statistical tests, to automatically detect key transcript features, including transcriptional start sites (TSSs), transcriptional termination sites (TTSs) and operon organization. Using our method, we obtained 2,764 TSSs and 1,467 TTSs for 1331 and 844 different genes, respectively. Identification of TSSs facilitated further discrimination of 215 putative sigma 38 regulons and 863 potential sigma 70 regulons. Combining the TSSs and TTSs with intergenic distance and co-expression information, we comprehensively annotated the operon organization in S. Typhimurium 14028s.

Conclusions

Our results show that directional RNA-seq can be used to detect transcriptional borders at an acceptable resolution of ±10-20 nucleotides. Technical limitations of the RNA-seq procedure may prevent single nucleotide resolution. The automatic transcript border detection methods, statistical models and operon organization pipeline that we have described could be widely applied to RNA-seq studies in other bacteria. Furthermore, the TSSs, TTSs, operons, promoters and unstranslated regions that we have defined for S. Typhimurium 14028s may constitute valuable resources that can be used for comparative analyses with other Salmonella serotypes.

Electronic supplementary material

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

Distinguishing between productive and abortive promoters using a random forest classifier in Mycoplasma pneumoniae

Distinguishing between promoter-like sequences in bacteria that belong to true or abortive promoters, or to those that do not initiate transcription at all, is one of the important challenges in transcriptomics. To address this problem, we have studied the genome-reduced bacterium Mycoplasma pneumoniae, for which the RNAs associated with transcriptional start sites have been recently experimentally identified. We determined the contribution to transcription events of different genomic features: the –10, extended –10 and –35 boxes, the UP element, the bases surrounding the –10 box and the nearest-neighbor free energy of the promoter region. Using a random forest classifier and the aforementioned features transformed into scores, we could distinguish between true, abortive promoters and non-promoters with good –10 box sequences. The methods used in this characterization of promoters can be extended to other bacteria and have important applications for promoter design in bacterial genome engineering.

Proteomic profile of dormancy within Staphylococcus epidermidis biofilms using iTRAQ and label-free strategies

Staphylococcus epidermidis is an important nosocomial bacterium among carriers of indwelling medical devices, since it has a strong ability to form biofilms. The presence of dormant bacteria within a biofilm is one of the factors that contribute to biofilm antibiotic tolerance and immune evasion. Here, we provide a detailed characterization of the quantitative proteomic profile of S. epidermidis biofilms with different proportions of dormant bacteria. A total of 427 and 409 proteins were identified by label-free and label-based quantitative methodologies, respectively. From these, 29 proteins were found to be differentially expressed between S. epidermidis biofilms with prevented and induced dormancy. Proteins overexpressed in S. epidermidis with prevented dormancy were associated with ribosome synthesis pathway, which reflects the metabolic state of dormant bacteria. In the opposite, underexpressed proteins were related to catalytic activity and ion binding, with involvement in purine, arginine, and proline metabolism. Additionally, GTPase activity seems to be enhanced in S. epidermidis biofilm with induced dormancy. The role of magnesium in dormancy modulation was further investigated with bioinformatics tool based in predicted interactions. The main molecular function of proteins, which strongly interact with magnesium, was nucleic acid binding. Different proteomic strategies allowed to obtain similar results and evidenced that prevented dormancy led to an expression of a markedly different repertoire of proteins in comparison to the one of dormant biofilms.