Conservation of transcription start sites within genes across a bacterial genus

Role of a single MCP in evolutionary adaptation of Shewanella putrefaciens for swimming in planktonic and structured environments

Bacteria can adapt to their environments by changing phenotypic traits by mutations. However, improving one trait often results in the deterioration of another one, a trade-off that limits the degree of adaptation. The gammaproteobacterium Shewanella putrefaciens CN-32 has an elaborate motility machinery comprising two distinct flagellar systems and an extensive chemotaxis array with 36 methyl-accepting chemotaxis sensor proteins (MCPs). In this study, we performed experimental selection on S. putrefaciens for increased spreading through a porous environment. We readily obtained a mutant that showed a pronounced increase in covered distance. This phenotype was almost completely caused by a deletion of 24 bp from the chromosome, which leads to a moderately enhanced production of a single MCP. Accordingly, chemotaxis assays under free-swimming conditions and cell tracking in soft agar showed that the mutation improved navigation through nutritional gradients. In contrast, further increased levels of the MCP negatively affected spreading. The study demonstrates how moderate differences in the abundance of a single MCP can lead to an efficient upgrade of chemotaxis in specific environments at a low expense of cellular resources.IMPORTANCEExperimental evolution experiments have been used to determine the trade-offs occurring in specific environments. Several studies that have used the spreading behavior of bacteria in structured environments identified regulatory mutants that increase the swimming speed of the cells. While this results in a higher chemotaxis drift, the growth fitness decreases as the higher swimming speed requires substantial cellular resources. Here we show that rapid chemotaxis adaptation can also be achieved by modifying the chemotaxis signal input at a low metabolic cost for the cell.

Non-canonical transcriptional start sites in E. coli O157:H7 EDL933 are regulated and appear in surprisingly high numbers

Analysis of genome wide transcription start sites (TSSs) revealed an unexpected complexity since not only canonical TSS of annotated genes are recognized by RNA polymerase. Non-canonical TSS were detected antisense to, or within, annotated genes as well new intergenic (orphan) TSS, not associated with known genes. Previously, it was hypothesized that many such signals represent noise or pervasive transcription, not associated with a biological function. Here, a modified Cappable-seq protocol allows determining the primary transcriptome of the enterohemorrhagic E. coli O157:H7 EDL933 (EHEC). We used four different growth media, both in exponential and stationary growth phase, replicated each thrice. This yielded 19,975 EHEC canonical and non-canonical TSS, which reproducibly occurring in three biological replicates. This questions the hypothesis of experimental noise or pervasive transcription. Accordingly, conserved promoter motifs were found upstream indicating proper TSSs. More than 50% of 5,567 canonical and between 32% and 47% of 10,355 non-canonical TSS were differentially expressed in different media and growth phases, providing evidence for a potential biological function also of non-canonical TSS. Thus, reproducible and environmentally regulated expression suggests that a substantial number of the non-canonical TSSs may be of unknown function rather than being the result of noise or pervasive transcription.

Transcriptome architecture of the three main lineages of agrobacteria

Agrobacteria are a diverse, polyphyletic group of prokaryotes with multipartite genomes capable of transferring DNA into the genomes of host plants, making them an essential tool in plant biotechnology. Despite their utility in plant transformation, genome-wide transcriptional regulation is not well understood across the three main lineages of agrobacteria. Transcription start sites (TSSs) are a necessary component of gene expression and regulation. In this study, we used differential RNA-seq and a TSS identification algorithm optimized on manually annotated TSS, then validated with existing TSS to identify thousands of TSS with nucleotide resolution for representatives of each lineage. We extend upon the 356 TSSs previously reported in Agrobacterium fabrum C58 by identifying 1,916 TSSs. In addition, we completed genomes and phenotyping of Rhizobium rhizogenes C16/80 and Allorhizobium vitis T60/94, identifying 2,650 and 2,432 TSSs, respectively. Parameter optimization was crucial for an accurate, high-resolution view of genome and transcriptional dynamics, highlighting the importance of algorithm optimization in genome-wide TSS identification and genomics at large. The optimized algorithm reduced the number of TSSs identified internal and antisense to the coding sequence on average by 90.5% and 91.9%, respectively. Comparison of TSS conservation between orthologs of the three lineages revealed differences in cell cycle regulation of ctrA as well as divergence of transcriptional regulation of chemotaxis-related genes when grown in conditions that simulate the plant environment. These results provide a framework to elucidate the mechanistic basis and evolution of pathology across the three main lineages of agrobacteria. IMPORTANCE Transcription start sites (TSSs) are fundamental for understanding gene expression and regulation. Agrobacteria, a group of prokaryotes with the ability to transfer DNA into the genomes of host plants, are widely used in plant biotechnology. However, the genome-wide transcriptional regulation of agrobacteria is not well understood, especially in less-studied lineages. Differential RNA-seq and an optimized algorithm enabled identification of thousands of TSSs with nucleotide resolution for representatives of each lineage. The results of this study provide a framework for elucidating the mechanistic basis and evolution of pathology across the three main lineages of agrobacteria. The optimized algorithm also highlights the importance of parameter optimization in genome-wide TSS identification and genomics at large.

Electrogenetic control of gene expression in Shewanella oneidensis MR-1 using Arc-dependent transcriptional promoters

Electrochemically active bacteria (EAB) are capable of electrically interacting with electrodes, enabling their application in bioelectrochemical systems (BESs). As the performance of BES is related to the metabolic activities of EAB, the development of methods to control their metabolic activities is important to facilitate BES applications. A recent study found that the EAB Shewanella oneidensis MR-1 uses the Arc system to regulate the expression of catabolic genes in response to electrode potentials, suggesting that a methodology for electrical control of gene expression in EAB, referred to as electrogenetics, can be developed by using electrode potential-responsive, Arc-dependent transcriptional promoters. Here, we explored Arc-dependent promoters in the genomes of S. oneidensis MR-1 and Escherichia coli to identify electrode potential-responsive promoters that are differentially activated in MR-1 cells exposed to high- and low-potential electrodes. LacZ reporter assays using electrode-associated cells of MR-1 derivatives revealed that the activities of promoters located upstream of the E. coli feo gene (P_feo) and the MR-1 nqrA2 (SO_0902) gene (P_nqr2) were significantly increased when S. oneidensis cells were exposed to electrodes poised at +0.7 V and -0.4 V (versus the standard hydrogen electrode), respectively. Additionally, we developed a microscopic system for in situ monitoring of promoter activity in electrode-associated cells and found that P_nqr2 activity was persistently induced in MR-1 cells associated with an electrode poised at -0.4 V. Our results indicate that these electrode potential-responsive promoters enable efficient regulation of gene expression in EAB, providing a molecular basis for the development of electrogenetics.

Sense and antisense RNA products of the uxuR gene can affect motility and chemotaxis acting independent of the UxuR protein

Small non-coding and antisense RNAs are widespread in all kingdoms of life, however, the diversity of their functions in bacteria is largely unknown. Here, we study RNAs synthesised from divergent promoters located in the 3'-end of the uxuR gene, encoding transcription factor regulating hexuronate metabolism in Escherichia coli. These overlapping promoters were predicted in silico with rather high scores, effectively bound RNA polymerase in vitro and in vivo and were capable of initiating transcription in sense and antisense directions. The genome-wide correlation between in silico promoter scores and RNA polymerase binding in vitro and in vivo was higher for promoters located on the antisense strands of the genes, however, sense promoters within the uxuR gene were more active. Both regulatory RNAs synthesised from the divergent promoters inhibited expression of genes associated with the E. coli motility and chemotaxis independent of a carbon source on which bacteria had been grown. Direct effects of these RNAs were confirmed for the fliA gene encoding σ²⁸ subunit of RNA polymerase. In addition to intracellular sRNAs, promoters located within the uxuR gene could initiate synthesis of transcripts found in the fraction of RNAs secreted in the extracellular medium. Their profile was also carbon-independent suggesting that intragenic uxuR transcripts have a specific regulatory role not directly related to the function of the protein in which gene they are encoded.

FlrA-independent production of flagellar proteins is required for proper flagellation in Shewanella putrefaciens

Flagella are multiprotein complexes whose assembly and positioning require complex spatiotemporal control. Flagellar assembly is thought to be controlled by several transcriptional tiers, which are mediated through various master regulators. Here, we revisited the regulation of flagellar genes in polarly flagellated gammaproteobacteria by the regulators FlrA, RpoN (σ⁵⁴ ) and FliA (σ²⁸ ) in Shewanella putrefaciens CN-32 at the transcript and protein level. We found that a number of regulatory and structural proteins were present in the absence of the main regulators, suggesting that initiation of flagella assembly and motor activation relies on the abundance control of only a few structural key components that are required for the formation of the MS- and C-ring and the flagellar type III secretion system. We identified FlrA-independent promoters driving expression of the regulators of flagellar number and positioning, FlhF and FlhG. Reduction of the gene expression levels from these promoters resulted in the emergence of hyperflagellation. This finding indicates that basal expression is required to adjust the flagellar counter in Shewanella. This is adding a deeper layer to the regulation of flagellar synthesis and assembly.

Precise Regulation of Differential Transcriptions of Various Catabolic Genes by OdcR via a Single Nucleotide Mutation in the Promoter Ensures the Safety of Metabolic Flux

Synergistic regulation of the expression of various genes in a catabolic pathway is crucial for the degradation, survival, and adaptation of microorganisms in polluted environments. However, how a single regulator accurately regulates and controls differential transcriptions of various catabolic genes to ensure metabolic safety remains largely unknown. Here, a LysR-type transcriptional regulator (LTTR), OdcR, encoded by the regulator gene odcR, was confirmed to be essential for 3,5-dibromo-4-hydroxybenozate (DBHB) catabolism and simultaneously activated the transcriptions of a gene with unknown function, orf419, and three genes, odcA, odcB, and odcC, involved in the DBHB catabolism in Pigmentiphaga sp. strain H8. OdcB further metabolized the highly toxic intermediate 2,6-dibromohydroquinone, which was produced from DBHB by OdcA. The upregulated transcriptional level of odcB was 7- to 9-fold higher than that of orf419, odcA, or odcC in response to DBHB. Through an electrophoretic mobility shift assay and DNase I footprinting assay, DBHB was found to be the effector and essential for OdcR binding to all four promoters of orf419, odcA, odcB, and odcC. A single nucleotide mutation in the regulatory binding site (RBS) of the promoter of odcB (TAT-N₁₁-ATG), compared to those of odcA/orf419 (CAT-N₁₁-ATG) and odcC (CAT-N₁₁-ATT), was identified and shown to enable the significantly higher transcription of odcB. The precise regulation of these genes by OdcR via a single nucleotide mutation in the promoter avoided the accumulation of 2,6-dibromohydroquinone, ensuring the metabolic safety of DBHB. IMPORTANCE Prokaryotes use various mechanisms, including improvement of the activity of detoxification enzymes, to cope with toxic intermediates produced during catabolism. However, studies on how bacteria accurately regulate differential transcriptions of various catabolic genes via a single regulator to ensure metabolic safety are scarce. This study revealed a LysR-type transcriptional activator, OdcR, which strongly activated odcB transcription for the detoxification of the toxic intermediate 2,6-dibromohydroquinone and slightly activated the transcriptions of other genes (orf419, odcA, and odcC) for 3,5-dibromo-4-hydroxybenozate (DBHB) catabolism in Pigmentiphaga sp. strain H8. Interestingly, the differential transcription/expression of the four genes, which ensured the metabolic safety of DBHB in cells, was determined by a single nucleotide mutation in the regulatory binding sites of the four promoters. This study describes a new and ingenious regulatory mode of ensuring metabolic safety in bacteria, expanding our understanding of synergistic transcriptional regulation in prokaryotes.

Genome-Wide TSS Distribution in Three Related Clostridia with Normalized Capp-Switch Sequencing

Transcription initiation is a tightly regulated process that is crucial for many aspects of prokaryotic physiology. High-throughput transcription start site (TSS) mapping can shed light on global and local regulation of transcription initiation, which in turn may help us understand and predict microbial behavior. In this study, we used Capp-Switch sequencing to determine the TSS positions in the genomes of three model solventogenic clostridia: Clostridium acetobutylicum ATCC 824, C. beijerinckii DSM 6423, and C. beijerinckii NCIMB 8052. We first refined the approach by implementing a normalization pipeline accounting for gene expression, yielding a total of 12,114 mapped TSSs across the species. We further compared the distributions of these sites in the three strains. Results indicated similar distribution patterns at the genome scale, but also some sharp differences, such as for the butyryl-CoA synthesis operon, particularly when comparing C. acetobutylicum to the C. beijerinckii strains. Lastly, we found that promoter structure is generally poorly conserved between C. acetobutylicum and C. beijerinckii. A few conserved promoters across species are discussed, showing interesting examples of how TSS determination and comparison can improve our understanding of gene expression regulation at the transcript level. IMPORTANCE Solventogenic clostridia have been employed in industry for more than a century, initially being used in the acetone-butanol-ethanol (ABE) fermentation process for acetone and butanol production. Interest in these bacteria has recently increased in the context of green chemistry and sustainable development. However, our current understanding of their genomes and physiology limits their optimal use as industrial solvent production platforms. The gene regulatory mechanisms of solventogenesis are still only partly understood, impeding efforts to increase rates and yields. Genome-wide mapping of transcription start sites (TSSs) for three model solventogenic Clostridium strains is an important step toward understanding mechanisms of gene regulation in these industrially important bacteria.

Promotech: a general tool for bacterial promoter recognition

Promoters are genomic regions where the transcription machinery binds to initiate the transcription of specific genes. Computational tools for identifying bacterial promoters have been around for decades. However, most of these tools were designed to recognize promoters in one or few bacterial species. Here, we present Promotech, a machine-learning-based method for promoter recognition in a wide range of bacterial species. We compare Promotech's performance with the performance of five other promoter prediction methods. Promotech outperforms these other programs in terms of area under the precision-recall curve (AUPRC) or precision at the same level of recall. Promotech is available at https://github.com/BioinformaticsLabAtMUN/PromoTech .

Design and in silico analysis of a whole-cell biosensor able to kill methicillin-resistant Staphylococcus aureus

The rise of methicillin-resistant Staphylococcus aureus (MRSA) infections has gained concern throughout the world over the past decades. Alternative therapeutic agents to antibiotics are rapidly growing to impede the proliferation of MRSA-caused infections. Lately, synthetic biology techniques have developed whole-cell biosensors by designing gene circuitry capable of sensing quorum-sensing (QS) molecules of pathogens and triggering expression of an antimicrobial moiety that kills MRSA and therefore prevents its further proliferation. Here, an E. coli was engineered in silico to act as a whole-cell biosensor that senses QS molecules from MRSA and triggers the expression of a bacteriocin that kills MRSA. To achieve this functionality, biosensor and bacteriocin modules were constructed and assembled into a vector. Both modules were codon-optimized to increase the yield production of the recombinant proteins. We then demonstrate in silico that the construction of a dual biosensor-killer plasmid, which holds two genetical modules known as biosensor and bacteriocin modules, enables the recombinant host to sense QS molecules from MRSA. Our designed whole-cell biosensor demonstrates in silico its ability to produce and secrete the bacteriocin as a function of the external concentration of autoinducer peptide from MRSA. These in silico results unravel the possibility of designing antimicrobial smarter therapeutics against resistant pathogens.

Redefining fundamental concepts of transcription initiation in bacteria

Despite enormous progress in understanding the fundamentals of bacterial gene regulation, our knowledge remains limited when compared with the number of bacterial genomes and regulatory systems to be discovered. Derived from a small number of initial studies, classic definitions for concepts of gene regulation have evolved as the number of characterized promoters has increased. Together with discoveries made using new technologies, this knowledge has led to revised generalizations and principles. In this Expert Recommendation, we suggest precise, updated definitions that support a logical, consistent conceptual framework of bacterial gene regulation, focusing on transcription initiation. The resulting concepts can be formalized by ontologies for computational modelling, laying the foundation for improved bioinformatics tools, knowledge-based resources and scientific communication. Thus, this work will help researchers construct better predictive models, with different formalisms, that will be useful in engineering, synthetic biology, microbiology and genetics.

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.

The Signal and the Noise: Characteristics of Antisense RNA in Complex Microbial Communities

This study systematically evaluated the global patterns of microbial antisense expression across various environments and provides a bird’s-eye view of general patterns observed across data sets, which can provide guidelines in our understanding of antisense expression as well as interpretation of metatranscriptomic data in general. This analysis highlights that in some environments, antisense expression from microbial communities can dominate over regular gene expression. We explored some potential drivers of antisense transcription, but more importantly, this study serves as a starting point, highlighting topics for future research and providing guidelines to include antisense expression in generic bioinformatic pipelines for metatranscriptomic data.

High-throughput sequencing has allowed unprecedented insight into the composition and function of complex microbial communities. With metatranscriptomics, it is possible to interrogate the transcriptomes of multiple organisms simultaneously to get an overview of the gene expression of the entire community. Studies have successfully used metatranscriptomics to identify and describe relationships between gene expression levels and community characteristics. However, metatranscriptomic data sets contain a rich suite of additional information that is just beginning to be explored. Here, we focus on antisense expression in metatranscriptomics, discuss the different computational strategies for handling it, and highlight the strengths but also potentially detrimental effects on downstream analysis and interpretation. We also analyzed the antisense transcriptomes of multiple genomes and metagenome-assembled genomes (MAGs) from five different data sets and found high variability in the levels of antisense transcription for individual species, which were consistent across samples. Importantly, we challenged the conceptual framework that antisense transcription is primarily the product of transcriptional noise and found mixed support, suggesting that the total observed antisense RNA in complex communities arises from the combined effect of unknown biological and technical factors. Antisense transcription can be highly informative, including technical details about data quality and novel insight into the biology of complex microbial communities.

IMPORTANCE This study systematically evaluated the global patterns of microbial antisense expression across various environments and provides a bird’s-eye view of general patterns observed across data sets, which can provide guidelines in our understanding of antisense expression as well as interpretation of metatranscriptomic data in general. This analysis highlights that in some environments, antisense expression from microbial communities can dominate over regular gene expression. We explored some potential drivers of antisense transcription, but more importantly, this study serves as a starting point, highlighting topics for future research and providing guidelines to include antisense expression in generic bioinformatic pipelines for metatranscriptomic data.

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.

Defining the Transcriptional and Post-transcriptional Landscapes of Mycobacterium smegmatis in Aerobic Growth and Hypoxia

The ability of Mycobacterium tuberculosis to infect, proliferate, and survive during long periods in the human lungs largely depends on the rigorous control of gene expression. Transcriptome-wide analyses are key to understanding gene regulation on a global scale. Here, we combine 5′-end-directed libraries with RNAseq expression libraries to gain insight into the transcriptome organization and post-transcriptional mRNA cleavage landscape in mycobacteria during log phase growth and under hypoxia, a physiologically relevant stress condition. Using the model organism Mycobacterium smegmatis, we identified 6,090 transcription start sites (TSSs) with high confidence during log phase growth, of which 67% were categorized as primary TSSs for annotated genes, and the remaining were classified as internal, antisense, or orphan, according to their genomic context. Interestingly, over 25% of the RNA transcripts lack a leader sequence, and of the coding sequences that do have leaders, 53% lack a strong consensus Shine-Dalgarno site. This indicates that like M. tuberculosis, M. smegmatis can initiate translation through multiple mechanisms. Our approach also allowed us to identify over 3,000 RNA cleavage sites, which occur at a novel sequence motif. To our knowledge, this represents the first report of a transcriptome-wide RNA cleavage site map in mycobacteria. The cleavage sites show a positional bias toward mRNA regulatory regions, highlighting the importance of post-transcriptional regulation in gene expression. We show that in low oxygen, a condition associated with the host environment during infection, mycobacteria change their transcriptomic profiles and endonucleolytic RNA cleavage is markedly reduced, suggesting a mechanistic explanation for previous reports of increased mRNA half-lives in response to stress. In addition, a number of TSSs were triggered in hypoxia, 56 of which contain the binding motif for the sigma factor SigF in their promoter regions. This suggests that SigF makes direct contributions to transcriptomic remodeling in hypoxia-challenged mycobacteria. Taken together, our data provide a foundation for further study of both transcriptional and posttranscriptional regulation in mycobacteria.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Transcriptional noise and exaptation as sources for bacterial sRNAs

Understanding how new genes originate and integrate into cellular networks is key to understanding evolution. Bacteria present unique opportunities for both the natural history and experimental study of gene origins, due to their large effective population sizes, rapid generation times, and ease of genetic manipulation. Bacterial small non-coding RNAs (sRNAs), in particular, many of which operate through a simple antisense regulatory logic, may serve as tractable models for exploring processes of gene origin and adaptation. Understanding how and on what timescales these regulatory molecules arise has important implications for understanding the evolution of bacterial regulatory networks, in particular, for the design of comparative studies of sRNA function. Here, we introduce relevant concepts from evolutionary biology and review recent work that has begun to shed light on the timescales and processes through which non-functional transcriptional noise is co-opted to provide regulatory functions. We explore possible scenarios for sRNA origin, focusing on the co-option, or exaptation, of existing genomic structures which may provide protected spaces for sRNA evolution.

Roles of d-Lactate Dehydrogenases in the Anaerobic Growth of Shewanella oneidensis MR-1 on Sugars

Shewanella oneidensis MR-1 is a facultative anaerobe that respires using a variety of electron acceptors. Although this organism is incapable of fermentative growth in the absence of electron acceptors, its genome encodes LdhA (a putative fermentative NADH-dependent d-lactate dehydrogenase [d-LDH]) and Dld (a respiratory quinone-dependent d-LDH). However, the physiological roles of LdhA in MR-1 are unclear. Here, we examined the activity, transcriptional regulation, and traits of deletion mutants to gain insight into the roles of LdhA in the anaerobic growth of MR-1. Analyses of d-LDH activity in MR-1 and the ldhA deletion mutant confirmed that LdhA functions as an NADH-dependent d-LDH that catalyzes the reduction of pyruvate to d-lactate. In vivo and in vitro assays revealed that ldhA expression was positively regulated by the cyclic-AMP receptor protein, a global transcription factor that regulates anaerobic respiratory pathways in MR-1, suggesting that LdhA functions in coordination with anaerobic respiration. Notably, we found that a deletion mutant of all four NADH dehydrogenases (NDHs) in MR-1 (ΔNDH mutant) retained the ability to grow on N-acetylglucosamine under fumarate-respiring conditions, while an additional deletion of ldhA or dld deprived the ΔNDH mutant of this growth ability. These results indicate that LdhA-Dld serves as a bypass of NDH in electron transfer from NADH to quinones. Our findings suggest that the LdhA-Dld system manages intracellular redox balance by utilizing d-lactate as a temporal electron sink under electron acceptor-limited conditions.IMPORTANCE NADH-dependent LDHs are conserved among diverse organisms and contribute to NAD⁺ regeneration in lactic acid fermentation. However, this type of LDH is also present in nonfermentative bacteria, including members of the genus Shewanella, while their physiological roles in these bacteria remain unknown. Here, we show that LdhA (an NADH-dependent d-LDH) works in concert with Dld (a quinone-dependent d-LDH) to transfer electrons from NADH to quinones during sugar catabolism in S. oneidensis MR-1. Our results indicate that d-lactate acts as an intracellular electron mediator to transfer electrons from NADH to membrane quinones. In addition, d-lactate serves as a temporal electron sink when respiratory electron acceptors are not available. Our study suggests novel physiological roles for d-LDHs in providing nonfermentative bacteria with catabolic flexibility under electron acceptor-limited conditions.

Internal RNAs overlapping coding sequences can drive the production of alternative proteins in archaea

Prokaryotic genomes show a high level of information compaction often with different molecules transcribed from the same locus. Although antisense RNAs have been relatively well studied, RNAs in the same strand, internal RNAs (intraRNAs), are still poorly understood. The question of how common is the translation of overlapping reading frames remains open. We address this question in the model archaeon Halobacterium salinarum. In the present work we used differential RNA-seq (dRNA-seq) in H. salinarum NRC-1 to locate intraRNA signals in subsets of internal transcription start sites (iTSS) and establish the open reading frames associated to them (intraORFs). Using C-terminally flagged proteins, we experimentally observed isoforms accurately predicted by intraRNA translation for kef1, acs3 and orc4 genes. We also recovered from the literature and mass spectrometry databases several instances of protein isoforms consistent with intraRNA translation such as the gas vesicle protein gene gvpC1. We found evidence for intraRNAs in horizontally transferred genes such as the chaperone dnaK and the aerobic respiration related cydA in both H. salinarum and Escherichia coli. Also, intraRNA translation evidence in H. salinarum, E. coli and yeast of a universal elongation factor (aEF-2, fusA and eEF-2) suggests that this is an ancient phenomenon present in all domains of life.

Widespread Antisense Transcription in Prokaryotes

Although bacterial genomes are usually densely protein-coding, genome-wide mapping approaches of transcriptional start sites revealed that a significant fraction of the identified promoters drive the transcription of noncoding RNAs. These can be trans-acting RNAs, mainly originating from intergenic regions and, in many studied examples, possessing regulatory functions. However, a significant fraction of these noncoding RNAs consist of natural antisense transcripts (asRNAs), which overlap other transcriptional units. Naturally occurring asRNAs were first observed to play a role in bacterial plasmid replication and in bacteriophage λ more than 30 years ago. Today's view is that asRNAs abound in all three domains of life. There are several examples of asRNAs in bacteria with clearly defined functions. Nevertheless, many asRNAs appear to result from pervasive initiation of transcription, and some data point toward global functions of such widespread transcriptional activity, explaining why the search for a specific regulatory role is sometimes futile. In this review, we give an overview about the occurrence of antisense transcription in bacteria, highlight particular examples of functionally characterized asRNAs, and discuss recent evidence pointing at global relevance in RNA processing and transcription-coupled DNA repair.

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.

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 Small RNA RyhB Is a Regulator of Cytochrome Expression in Shewanella oneidensis

Shewanella oneidensis produces an extensive electron transfer network that results in metabolic flexibility. A large number of c-type cytochromes are expressed by S. oneidensis and these function as the fundamental electron transport chain proteins. Although several S. oneidensis cytochromes have been well-characterized, little is known about how their expression is regulated. In this study, we investigate the role of the ferric uptake regulator (Fur) and the sRNA RyhB in regulation. Our results demonstrate that loss of Fur leads to diminished growth and an apparent decrease in heme-containing proteins. Remarkably, deleting the Fur-repressed ryhB gene almost completely reverses these physiological changes, indicating that the phenotypes resulting from loss of Fur are (at least partially) dependent on RyhB. RNA sequencing identified a number of possible RyhB repressed genes. A large fraction of these encode c-type cytochromes, among them two of the most abundant periplasmic cytochromes CctA (also known as STC) and ScyA. We show that RyhB destabilizes the mRNA of four of its target genes, cctA, scyA, omp35, and nrfA and this requires the presence of the RNA chaperone Hfq. Iron limitation decreases the expression of the RyhB target genes cctA and scyA and this regulation relies on the presence of both Fur and RyhB. Overall, this study suggests that controlling cytochrome expression is of importance to maintain iron homeostasis and that sRNAs molecules are important players in the regulation of fundamental processes in S. oneidensis MR-1.

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.

Nitrogen cost minimization is promoted by structural changes in the transcriptome of N-deprived Prochlorococcus cells

Prochlorococcus is a globally abundant marine cyanobacterium with many adaptations that reduce cellular nutrient requirements, facilitating growth in its nutrient-poor environment. One such genomic adaptation is the preferential utilization of amino acids containing fewer N-atoms, which minimizes cellular nitrogen requirements. We predicted that transcriptional regulation might further reduce cellular N budgets during transient N limitation. To explore this, we compared transcription start sites (TSSs) in Prochlorococcus MED4 under N-deprived and N-replete conditions. Of 64 genes with primary and internal TSSs in both conditions, N-deprived cells initiated transcription downstream of primary TSSs more frequently than N-replete cells. Additionally, 117 genes with only an internal TSS demonstrated increased internal transcription under N-deprivation. These shortened transcripts encode predicted proteins with an average of 21% less N content compared to full-length transcripts. We hypothesized that low translation rates, which afford greater control over protein abundances, would be beneficial to relatively slow-growing organisms like Prochlorococcus. Consistent with this idea, we found that Prochlorococcus exhibits greater usage of glycine–glycine motifs, which causes translational pausing, when compared to faster growing microbes. Our findings indicate that structural changes occur within the Prochlorococcus MED4 transcriptome during N-deprivation, potentially altering the size and structure of proteins expressed under nutrient limitation.

Validating regulatory predictions from diverse bacteria with mutant fitness data

Although transcriptional regulation is fundamental to understanding bacterial physiology, the targets of most bacterial transcription factors are not known. Comparative genomics has been used to identify likely targets of some of these transcription factors, but these predictions typically lack experimental support. Here, we used mutant fitness data, which measures the importance of each gene for a bacterium’s growth across many conditions, to test regulatory predictions from RegPrecise, a curated collection of comparative genomics predictions. Because characterized transcription factors often have correlated fitness with one of their targets (either positively or negatively), correlated fitness patterns provide support for the comparative genomics predictions. At a false discovery rate of 3%, we identified significant cofitness for at least one target of 158 TFs in 107 ortholog groups and from 24 bacteria. Thus, high-throughput genetics can be used to identify a high-confidence subset of the sequence-based regulatory predictions.

Global repositioning of transcription start sites in a plant-fermenting bacterium

Bacteria respond to their environment by regulating mRNA synthesis, often by altering the genomic sites at which RNA polymerase initiates transcription. Here, we investigate genome-wide changes in transcription start site (TSS) usage by Clostridium phytofermentans, a model bacterium for fermentation of lignocellulosic biomass. We quantify expression of nearly 10,000 TSS at single base resolution by Capp-Switch sequencing, which combines capture of synthetically capped 5' mRNA fragments with template-switching reverse transcription. We find the locations and expression levels of TSS for hundreds of genes change during metabolism of different plant substrates. We show that TSS reveals riboswitches, non-coding RNA and novel transcription units. We identify sequence motifs associated with carbon source-specific TSS and use them for regulon discovery, implicating a LacI/GalR protein in control of pectin metabolism. We discuss how the high resolution and specificity of Capp-Switch enables study of condition-specific changes in transcription initiation in bacteria.

Genome-wide identification of transcriptional start sites in the haloarchaeon Haloferax volcanii based on differential RNA-Seq (dRNA-Seq)

BACKGROUND

Differential RNA-Seq (dRNA-Seq) is a recently developed method of performing primary transcriptome analyses that allows for the genome-wide mapping of transcriptional start sites (TSSs) and the identification of novel transcripts. Although the transcriptomes of diverse bacterial species have been characterized by dRNA-Seq, the transcriptome analysis of archaeal species is still rather limited. Therefore, we used dRNA-Seq to characterize the primary transcriptome of the model archaeon Haloferax volcanii.

RESULTS

Three independent cultures of Hfx. volcanii grown under optimal conditions to the mid-exponential growth phase were used to determine the primary transcriptome and map the 5'-ends of the transcripts. In total, 4749 potential TSSs were detected. A position weight matrix (PWM) was derived for the promoter predictions, and the results showed that 64 % of the TSSs were preceded by stringent or relaxed basal promoters. Of the identified TSSs, 1851 belonged to protein-coding genes. Thus, fewer than half (46 %) of the 4040 protein-coding genes were expressed under optimal growth conditions. Seventy-two percent of all protein-coding transcripts were leaderless, which emphasized that this pathway is the major pathway for translation initiation in haloarchaea. A total of 2898 of the TSSs belonged to potential non-coding RNAs, which accounted for an unexpectedly high fraction (61 %) of all transcripts. Most of the non-coding TSSs had not been previously described (2792) and represented novel sequences (59 % of all TSSs). A large fraction of the potential novel non-coding transcripts were cis-antisense RNAs (1244 aTSSs). A strong negative correlation between the levels of antisense transcripts and cognate sense mRNAs was found, which suggested that the negative regulation of gene expression via antisense RNAs may play an important role in haloarchaea. The other types of novel non-coding transcripts corresponded to internal transcripts overlapping with mRNAs (1153 iTSSs) and intergenic small RNA (sRNA) candidates (395 TSSs).

CONCLUSION

This study provides a comprehensive map of the primary transcriptome of Hfx. volcanii grown under optimal conditions. Fewer than half of all protein-coding genes have been transcribed under these conditions. Unexpectedly, more than half of the detected TSSs belonged to several classes of non-coding RNAs. Thus, RNA-based regulation appears to play a more important role in haloarchaea than previously anticipated.

Comparative transcriptomics across the prokaryotic tree of life

Whole-transcriptome sequencing studies from recent years revealed an unexpected complexity in transcriptomes of bacteria and archaea, including abundant non-coding RNAs, cis-antisense transcription and regulatory untranslated regions (UTRs). Understanding the functional relevance of the plethora of non-coding RNAs in a given organism is challenging, especially since some of these RNAs were attributed to ‘transcriptional noise’. To allow the search for conserved transcriptomic elements we produced comparative transcriptome maps for multiple species across the microbial tree of life. These transcriptome maps are detailed in annotations, comparable by gene families, and BLAST-searchable by user provided sequences. Our transcriptome collection includes 18 model organisms spanning 10 phyla/subphyla of bacteria and archaea that were sequenced using standardized RNA-seq methods. The utility of the comparative approach, as implemented in our web server, is demonstrated by highlighting genes with exceptionally long 5′UTRs across species, which correspond to many known riboswitches and further suggest novel putative regulatory elements. Our study provides a standardized reference transcriptome to major clinically and environmentally important microbial phyla. The viewer is available at http://exploration.weizmann.ac.il/TCOL, setting a framework for comparative studies of the microbial non-coding genome.

Genome-wide transcription start site mapping of Bradyrhizobium japonicum grown free-living or in symbiosis - a rich resource to identify new transcripts, proteins and to study gene regulation

BACKGROUND

Differential RNA-sequencing (dRNA-seq) is indispensable for determination of primary transcriptomes. However, using dRNA-seq data to map transcriptional start sites (TSSs) and promoters genome-wide is a bioinformatics challenge. We performed dRNA-seq of Bradyrhizobium japonicum USDA 110, the nitrogen-fixing symbiont of soybean, and developed algorithms to map TSSs and promoters.

RESULTS

A specialized machine learning procedure for TSS recognition allowed us to map 15,923 TSSs: 14,360 in free-living bacteria, 4329 in symbiosis with soybean and 2766 in both conditions. Further, we provide proteomic evidence for 4090 proteins, among them 107 proteins corresponding to new genes and 178 proteins with N-termini different from the existing annotation (72 and 109 of them with TSS support, respectively). Guided by proteomics evidence, previously identified TSSs and TSSs experimentally validated here, we assign a score threshold to flag 14 % of the mapped TSSs as a class of lower confidence. However, this class of lower confidence contains valid TSSs of low-abundant transcripts. Moreover, we developed a de novo algorithm to identify promoter motifs upstream of mapped TSSs, which is publicly available, and found motifs mainly used in symbiosis (similar to RpoN-dependent promoters) or under both conditions (similar to RpoD-dependent promoters). Mapped TSSs and putative promoters, proteomic evidence and updated gene annotation were combined into an annotation file.

CONCLUSIONS

The genome-wide TSS and promoter maps along with the extended genome annotation of B. japonicum represent a valuable resource for future systems biology studies and for detailed analyses of individual non-coding transcripts and ORFs. Our data will also provide new insights into bacterial gene regulation during the agriculturally important symbiosis between rhizobia and legumes.

Staphylococcus aureus Transcriptome Architecture: From Laboratory to Infection-Mimicking Conditions

Staphylococcus aureus is a major pathogen that colonizes about 20% of the human population. Intriguingly, this Gram-positive bacterium can survive and thrive under a wide range of different conditions, both inside and outside the human body. Here, we investigated the transcriptional adaptation of S. aureus HG001, a derivative of strain NCTC 8325, across experimental conditions ranging from optimal growth in vitro to intracellular growth in host cells. These data establish an extensive repertoire of transcription units and non-coding RNAs, a classification of 1412 promoters according to their dependence on the RNA polymerase sigma factors SigA or SigB, and allow identification of new potential targets for several known transcription factors. In particular, this study revealed a relatively low abundance of antisense RNAs in S. aureus, where they overlap only 6% of the coding genes, and only 19 antisense RNAs not co-transcribed with other genes were found. Promoter analysis and comparison with Bacillus subtilis links the small number of antisense RNAs to a less profound impact of alternative sigma factors in S. aureus. Furthermore, we revealed that Rho-dependent transcription termination suppresses pervasive antisense transcription, presumably originating from abundant spurious transcription initiation in this A+T-rich genome, which would otherwise affect expression of the overlapped genes. In summary, our study provides genome-wide information on transcriptional regulation and non-coding RNAs in S. aureus as well as new insights into the biological function of Rho and the implications of spurious transcription in bacteria.

The major human pathogen Staphylococcus aureus can survive under a wide range of conditions, both inside and outside the human body. The goal of this study was to determine how S. aureus adapts to such different conditions and, additionally, we wanted to identify general factors governing the staphylococcal transcriptome architecture. Therefore, we performed a precise analysis of all RNA transcripts of S. aureus across experimental conditions ranging from in vitro growth in different media to internalization by eukaryotic host cells. We systematically mapped all transcription units, annotated non-coding RNAs, and assigned promoters controlled by particular RNA polymerase sigma factors and transcription factors. By a comparison with data available for the related Gram-positive bacterium Bacillus subtilis, we made key observations concerning the abundance and origin of antisense RNAs. Intriguingly, these findings support the view that many antisense RNAs in a bacterium like B. subtilis could be byproducts of spurious promoter recognition by condition-specific alternative sigma factors. We also report that the transcription termination factor Rho prevents widespread antisense transcription, presumably caused by pervasive transcription initiation in the A+T-rich genome of S. aureus. Altogether our study presents new perspectives on the biological significance of antisense and pervasive transcription in bacteria.

Sense overlapping transcripts in IS1341-type transposase genes are functional non-coding RNAs in archaea

The existence of sense overlapping transcripts that share regulatory and coding information in the same genomic sequence shows an additional level of prokaryotic gene expression complexity. Here we report the discovery of ncRNAs associated with IS1341-type transposase (tnpB) genes, at the 3'-end of such elements, with examples in archaea and bacteria. Focusing on the model haloarchaeon Halobacterium salinarum NRC-1, we show the existence of sense overlapping transcripts (sotRNAs) for all its IS1341-type transposases. Publicly available transcriptome compendium show condition-dependent differential regulation between sotRNAs and their cognate genes. These sotRNAs allowed us to find a UUCA tetraloop motif that is present in other archaea (ncRNA family HgcC) and in a H. salinarum intergenic ncRNA derived from a palindrome associated transposable elements (PATE). Overexpression of one sotRNA and the PATE-derived RNA harboring the tetraloop motif improved H. salinarum growth, indicating that these ncRNAs are functional.

The complexity of bacterial transcriptomes

For eukaryotes there seems to be no doubt that differences on the trancriptomic level substantially contribute to the process of species diversification, whereas for bacteria this is thought to be less important. Recent years saw a significant increase in full transcriptome studies for bacteria, which provided deep insight into the architecture of bacterial transcriptomes. Most notably, it became evident that, in contrast to previous scientific consensus, bacterial transcriptomes are quite complex. There exist a large number of cis-antisense RNAs, non-coding RNAs, overlapping transcripts and RNA elements that regulate transcription, such as riboswitches. Furthermore, processing and degradation of RNA has gained interest, because it has a significant impact on the composition of the transcriptome. In this review, we summarize recent findings and put them into a broader context with respect to the complexity of bacterial transcriptomes and its putative biological meanings.

Catabolic and regulatory systems in Shewanella oneidensis MR-1 involved in electricity generation in microbial fuel cells

Shewanella oneidensis MR-1 is a facultative anaerobe that respires using a variety of inorganic and organic compounds. MR-1 is also capable of utilizing extracellular solid materials, including anodes in microbial fuel cells (MFCs), as electron acceptors, thereby enabling electricity generation. As MFCs have the potential to generate electricity from biomass waste and wastewater, MR-1 has been extensively studied to identify the molecular systems that are involved in electricity generation in MFCs. These studies have demonstrated the importance of extracellular electron-transfer (EET) pathways that electrically connect the quinone pool in the cytoplasmic membrane to extracellular electron acceptors. Electricity generation is also dependent on intracellular catabolic pathways that oxidize electron donors, such as lactate, and regulatory systems that control the expression of genes encoding the components of catabolic and electron-transfer pathways. In addition, recent findings suggest that cell-surface polymers, e.g., exopolysaccharides, and secreted chemicals, which function as electron shuttles, are also involved in electricity generation. Despite these advances in our knowledge on the EET processes in MR-1, further efforts are necessary to fully understand the underlying intra- and extracellular molecular systems for electricity generation in MFCs. We suggest that investigating how MR-1 coordinates these systems to efficiently transfer electrons to electrodes and conserve electrochemical energy for cell proliferation is important for establishing the biological basis for MFCs.

Transcriptional mechanisms for differential expression of outer membrane cytochrome genes omcA and mtrC in Shewanella oneidensis MR-1

Background

Shewanella oneidensis MR-1 is capable of reducing extracellular electron acceptors, such as metals and electrodes, through the Mtr respiratory pathway, which consists of the outer membrane cytochromes OmcA and MtrC and associated proteins MtrA and MtrB. These proteins are encoded in the mtr gene cluster (omcA-mtrCAB) in the MR-1 chromosome.

Results

Here, we investigated the transcriptional mechanisms for the mtr genes and demonstrated that omcA and mtrC are transcribed from two upstream promoters, P_omcA and P_mtrC, respectively. In vivo transcription and in vitro electrophoretic mobility shift assays revealed that a cAMP receptor protein (CRP) positively regulates the expression of the mtr genes by binding to the upstream regions of P_omcA and P_mtrC. However, the expression of omcA and mtrC was differentially regulated in response to culture conditions; specifically, the expression from P_mtrC was higher under aerobic conditions than that under anaerobic conditions with fumarate as an electron acceptor, whereas expression from P_omcA exhibited the opposite trend. Deletion of the region upstream of the CRP-binding site of P_omcA resulted in a significant increase in promoter activity under aerobic conditions, demonstrating that the deleted region is involved in the negative regulation of P_omcA.

Conclusions

Taken together, the present results indicate that transcription of the mtr genes is regulated by multiple promoters and regulatory systems, including the CRP/cAMP-dependent regulatory system and yet-unidentified negative regulators.

Electronic supplementary material

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

Quantitative bacterial transcriptomics with RNA-seq

RNA sequencing has emerged as the premier approach to study bacterial transcriptomes. While the earliest published studies analyzed the data qualitatively, the data are readily digitized and lend themselves to quantitative analysis. High-resolution RNA sequence (RNA-seq) data allows transcriptional features (promoters, terminators, operons, among others) to be pinpointed on any bacterial transcriptome. Once the transcriptome is mapped, the activity of transcriptional features can be quantified. Here we highlight how quantitative transcriptome analysis can reveal biological insights and briefly discuss some of the challenges to be faced by the field of bacterial transcriptomics in the near future.

Where to begin? Mapping transcription start sites genome-wide in Escherichia coli

Recent genome-wide studies of bacterial transcription have revealed large numbers of promoters located inside genes. In this issue of the Journal of Bacteriology, Thomason and colleagues (J. Bacteriol. 197:18-28, 2015, doi:10.1128/JB.02096-14) map transcription start sites in Escherichia coli on an unprecedented scale. This work provides important insights into the regulation of transcripts that initiate inside genes and sources of variability between studies aimed at identifying these RNAs.

Global transcriptional start site mapping using differential RNA sequencing reveals novel antisense RNAs in Escherichia coli

While the model organism Escherichia coli has been the subject of intense study for decades, the full complement of its RNAs is only now being examined. Here we describe a survey of the E. coli transcriptome carried out using a differential RNA sequencing (dRNA-seq) approach, which can distinguish between primary and processed transcripts, and an automated prediction algorithm for transcriptional start sites (TSS). With the criterion of expression under at least one of three growth conditions examined, we predicted 14,868 TSS candidates, including 5,574 internal to annotated genes (iTSS) and 5,495 TSS corresponding to potential antisense RNAs (asRNAs). We examined expression of 14 candidate asRNAs by Northern analysis using RNA from wild-type E. coli and from strains defective for RNases III and E, two RNases reported to be involved in asRNA processing. Interestingly, nine asRNAs detected as distinct bands by Northern analysis were differentially affected by the rnc and rne mutations. We also compared our asRNA candidates with previously published asRNA annotations from RNA-seq data and discuss the challenges associated with these cross-comparisons. Our global transcriptional start site map represents a valuable resource for identification of transcription start sites, promoters, and novel transcripts in E. coli and is easily accessible, together with the cDNA coverage plots, in an online genome browser.