Ribosomes bind leaderless mRNA in Escherichia coli through recognition of their 5'-terminal AUG

Unraveling the plasticity of translation initiation in prokaryotes: Beyond the invariant Shine-Dalgarno sequence

Translation initiation in prokaryotes is mainly defined, although not exclusively, by the interaction between the anti-Shine-Dalgarno sequence (antiSD), located at the 3'-terminus of the 16S ribosomal RNA, and a complementary sequence, the ribosome binding site, or Shine-Dalgarno (SD), located upstream of the start codon in prokaryotic mRNAs. The antiSD has a conserved 5'-CCUCC-3' core, but inter-species variations have been found regarding the participation of flanking bases in binding. These variations have been described for certain bacteria and, to a lesser extent, for some archaea. To further analyze these variations, we conducted binding-energy prediction analyses on over 6,400 genomic sequences from both domains. We identified 15 groups of antiSD variants that could be associated with the organisms' phylogenetic origin. Additionally, our findings revealed that certain organisms exhibit variations in the core itself. Importantly, an unaltered core is not necessarily required for the interaction between the 3'-terminus of the rRNA and the region preceding the AUG of the mRNA. In our study, we classified organisms into four distinct categories: i) those possessing a conserved core and demonstrating binding; ii) those with a conserved core but lacking evidence of binding; iii) those exhibiting binding in the absence of a conserved core; and iv) those lacking both a conserved core and evidence of binding. Our results demonstrate the flexibility of organisms in evolving different sequences involved in translation initiation beyond the traditional Shine-Dalgarno sequence. These findings are discussed in terms of the evolution of translation initiation in prokaryotic organisms.

Factors impacting the regulation of nos gene expression in Staphylococcus aureus

Staphylococcus aureus nitric oxide synthase (saNOS) contributes to oxidative stress resistance, antibiotic tolerance, virulence, and modulation of aerobic and nitrate-based cellular respiration. Despite its involvement in these essential processes, the genetic regulation of nos expression has not been well characterized. 5' rapid amplification of cDNA ends on nos RNA isolated from S. aureus UAMS-1 (USA200 strain) and AH1263 (USA300 strain) revealed that the nos transcriptional start site mapped to an adenine nucleotide in the predicted Shine-Dalgarno site located 11 bp upstream of the nos ATG start codon, suggesting that the nos transcript may have a leaderless organization or may be subject to processing. The SrrAB two-component system (TCS) was previously identified as a positive regulator of nos RNA levels, and experiments using a β-galactosidase reporter plasmid confirmed that SrrAB is a positive regulator of nos promoter activity. In addition, the quorum-sensing system Agr was identified as a negative regulator of low-oxygen nos expression in UAMS-1, with activity epistatic to SrrAB. Involvement of Agr was strain dependent, as nos expression remained unchanged in an AH1263 agr mutant, which has higher Agr activity compared to UAMS-1. Furthermore, nos promoter activity and RNA levels were significantly stronger in AH1263 relative to UAMS-1 during late-exponential low-oxygen growth, when nos expression is maximal. Global regulators Rex and MgrA were also implicated as negative regulators of low-oxygen nos promoter activity in UAMS-1. Collectively, these results provide new insight into factors that control nos expression.IMPORTANCEBacterial nitric oxide synthase (bNOS) has recently emerged in several species as a key player in resistance to stresses commonly encountered during infection. Although Staphylococcus aureus (sa)NOS has been suggested to be a promising drug target in S. aureus, an obstacle to this in practice is the existence of mammalian NOS, whose oxygenase domain is like bacterial NOS. Increased understanding of the nos regulatory network in S. aureus could allow targeting of saNOS through its regulators, bypassing the issue of also inhibiting mammalian NOS. Furthermore, the observed strain-dependent differences in S. aureus nos regulation presented in this study reinforce the importance of studying bacterial NOS regulation and function at both the strain and species levels.

Genome-wide transcription start site mapping in the facultative intracellular pathogen Brucella melitensis by Capping-seq

Brucella melitensis (B. melitensis) is an important facultative intracellular bacterium that causes global zoonotic diseases. Continuous intracellular survival and replication are the main obstruction responsible for the accessibility of prevention and treatment of brucellosis. Bacteria respond to complex environment by regulating gene expression. Many regulatory factors function at loci where RNA polymerase initiates messenger RNA synthesis. However, limited gene annotation is a current obstacle for the research on expression regulation in bacteria. To improve annotation and explore potential functional sites, we proposed a novel genome-wide method called Capping-seq for transcription start site (TSS) mapping in B. melitensis. This technique combines capture of capped primary transcripts with Single Molecule Real-Time (SMRT) sequencing technology. We identified 2,369 TSSs at single nucleotide resolution by Capping-seq. TSSs analysis of Brucella transcripts showed a preference of purine on the TSS positions. Our results revealed that -35 and -10 elements of promoter contained consensus sequences of TTGNNN and TATNNN, respectively. The 5' ends analysis showed that 57% genes are associated with more than one TSS and 47% genes contain long leader regions, suggested potential complex regulation at the 5' ends of genes in B. melitensis. Moreover, we identified 52 leaderless genes that are mainly involved in the metabolic processes. Overall, Capping-seq technology provides a unique solution for TSS determination in prokaryotes. Our findings develop a systematic insight into the primary transcriptome characterization of B. melitensis. This study represents a critical basis for investigating gene regulation and pathogenesis of Brucella.

Bistable Expression of a Toxin-Antitoxin System Located in a Cryptic Prophage of Escherichia coli O157:H7

Type II toxin-antitoxin (TA) systems are classically composed of two genes that encode a toxic protein and a cognate antitoxin protein. Both genes are organized in an operon whose expression is autoregulated at the level of transcription by the antitoxin-toxin complex, which binds operator DNA through the antitoxin’s DNA-binding domain. Here, we investigated the transcriptional regulation of a particular TA system located in the immunity region of a cryptic lambdoid prophage in the Escherichia coli O157:H7 EDL933 strain. This noncanonical paaA2-parE2 TA operon contains a third gene, paaR2, that encodes a transcriptional regulator that was previously shown to control expression of the TA. We provide direct evidence that the PaaR2 is a transcriptional regulator which shares functional similarities to the lambda CI repressor. Expression of the paaA2-parE2 TA operon is regulated by two other transcriptional regulators, YdaS and YdaT, encoded within the same region. We argue that YdaS and YdaT are analogous to lambda Cro and CII and that they do not constitute a TA system, as previously debated. We show that PaaR2 primarily represses the expression of YdaS and YdaT, which in turn controls the expression of paaR2-paaA2-parE2 operon. Overall, our results show that the paaA2-parE2 TA is embedded in an intricate lambdoid prophage-like regulation network. Using single-cell analysis, we observed that the entire locus exhibits bistability, which generates diversity of expression in the population. Moreover, we confirmed that paaA2-parE2 is addictive and propose that it could limit genomic rearrangements within the immunity region of the CP-933P cryptic prophage.

A combination of mRNA features influence the efficiency of leaderless mRNA translation initiation

Bacterial translation is thought to initiate by base pairing of the 16S rRNA and the Shine-Dalgarno sequence in the mRNA's 5' untranslated region (UTR). However, transcriptomics has revealed that leaderless mRNAs, which completely lack any 5' UTR, are broadly distributed across bacteria and can initiate translation in the absence of the Shine-Dalgarno sequence. To investigate the mechanism of leaderless mRNA translation initiation, synthetic in vivo translation reporters were designed that systematically tested the effects of start codon accessibility, leader length, and start codon identity on leaderless mRNA translation initiation. Using these data, a simple computational model was built based on the combinatorial relationship of these mRNA features that can accurately classify leaderless mRNAs and predict the translation initiation efficiency of leaderless mRNAs. Thus, start codon accessibility, leader length, and start codon identity combine to define leaderless mRNA translation initiation in bacteria.

Autorepressor PsrA is required for optimal Legionella pneumophila growth in Acanthamoeba castellanii protozoa

Legionella pneumophila possesses a unique intracellular lifecycle featuring distinct morphological stages that include replicative forms and transmissive cyst forms. Expression of genes associated with virulence traits and cyst morphogenesis is concomitant, and governed by a complex stringent response based-regulatory network and the stationary phase sigma factor RpoS. In Pseudomonas spp., rpoS expression is controlled by the autorepressor PsrA, and orthologs of PsrA and RpoS are required for cyst formation in Azotobacter. Here we report that the L. pneumophila psrA ortholog, expressed as a leaderless monocistronic transcript, is also an autorepressor, but is not a regulator of rpoS expression. Further, the binding site sequence recognized by L. pneumophila PsrA is different from that of Pseudomonas PsrA, suggesting a repertoire of target genes unique to L. pneumophila. While PsrA was dispensable for growth in human U937-derived macrophages, lack of PsrA affected bacterial intracellular growth in Acanthamoeba castellanii protozoa, but also increased the quantity of poly-3-hydroxybutyrate (PHB) inclusions in matured transmissive cysts. Interestingly, overexpression of PsrA increased the size and bacterial load of the replicative vacuole in both host cell types. Taken together, we report that PsrA is a host-specific requirement for optimal temporal progression of L. pneumophila intracellular lifecycle in A. castellanii.

The diversity of Shine-Dalgarno sequences sheds light on the evolution of translation initiation

Shine-Dalgarno (SD) sequences, the core element of prokaryotic ribosome-binding sites, facilitate mRNA translation by base-pair interaction with the anti-SD (aSD) sequence of 16S rRNA. In contrast to this paradigm, an inspection of thousands of prokaryotic species unravels tremendous SD sequence diversity both within and between genomes, whereas aSD sequences remain largely static. The pattern has led many to suggest unidentified mechanisms for translation initiation. Here we review known translation-initiation pathways in prokaryotes. Moreover, we seek to understand the cause and consequence of SD diversity through surveying recent advances in biochemistry, genomics, and high-throughput genetics. These findings collectively show: (1) SD:aSD base pairing is beneficial but nonessential to translation initiation. (2) The 5' untranslated region of mRNA evolves dynamically and correlates with organismal phylogeny and ecological niches. (3) Ribosomes have evolved distinct usage of translation-initiation pathways in different species. We propose a model portraying the SD diversity shaped by optimization of gene expression, adaptation to environments and growth demands, and the species-specific prerequisite of ribosomes to initiate translation. The model highlights the coevolution of ribosomes and mRNA features, leading to functional customization of the translation apparatus in each organism.

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

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

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

Leaderless mRNAs in the Spotlight: Ancient but Not Outdated!

Previously, leaderless mRNAs (lmRNAs) were perceived to make up only a minor fraction of the transcriptome in bacteria. However, advancements in RNA sequencing technology are uncovering vast numbers of lmRNAs, particularly in archaea, Actinobacteria , and extremophiles and thus underline their significance in cellular physiology and regulation. Due to the absence of conventional ribosome binding signals, lmRNA translation initiation is distinct from canonical mRNAs and can therefore be differentially regulated. The ribosome’s inherent ability to bind a 5′-terminal AUG can stabilize and protect the lmRNA from degradation or allow ribosomal loading for downstream initiation events. As a result, lmRNAs remain translationally competent during a variety of physiological conditions, allowing them to contribute to multiple regulatory mechanisms. Furthermore, the abundance of lmRNAs can increase during adverse conditions through the upregulation of lmRNA transcription from alternative promoters or by the generation of lmRNAs from canonical mRNAs cleaved by an endonucleolytic toxin. In these ways, lmRNA translation can continue during stress and contribute to regulation, illustrating their importance in the cell. Due to their presence in all domains of life and their ability to be translated by heterologous hosts, lmRNAs appear further to represent ancestral transcripts that might allow us to study the evolution of the ribosome and the translational process.

RNAseq analysis of α-proteobacterium Gluconobacter oxydans 621H

Background

The acetic acid bacterium Gluconobacter oxydans 621H is characterized by its exceptional ability to incompletely oxidize a great variety of carbohydrates in the periplasm. The metabolism of this α-proteobacterium has been characterized to some extent, yet little is known about its transcriptomes and related data. In this study, we applied two different RNAseq approaches. Primary transcriptomes enriched for 5′-ends of transcripts were sequenced to detect transcription start sites, which allow subsequent analysis of promoter motifs, ribosome binding sites, and 5´-UTRs. Whole transcriptomes were sequenced to identify expressed genes and operon structures.

Results

Sequencing of primary transcriptomes of G. oxydans revealed 2449 TSSs, which were classified according to their genomic context followed by identification of promoter and ribosome binding site motifs, analysis of 5´-UTRs including validation of predicted cis-regulatory elements and correction of start codons. 1144 (41%) of all genes were found to be expressed monocistronically, whereas 1634 genes were organized in 571 operons. Together, TSSs and whole transcriptome data were also used to identify novel intergenic (18), intragenic (328), and antisense transcripts (313).

Conclusions

This study provides deep insights into the transcriptional landscapes of G. oxydans. The comprehensive transcriptome data, which we made publicly available, facilitate further analysis of promoters and other regulatory elements. This will support future approaches for rational strain development and targeted gene expression in G. oxydans. The corrections of start codons further improve the high quality genome reference and support future proteome analysis.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4415-x) contains supplementary material, which is available to authorized users.

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.

Genome-wide primary transcriptome analysis of H2-producing archaeon Thermococcus onnurineus NA1

In spite of their pivotal roles in transcriptional and post-transcriptional processes, the regulatory elements of archaeal genomes are not yet fully understood. Here, we determine the primary transcriptome of the H2-producing archaeon Thermococcus onnurineus NA1. We identified 1,082 purine-rich transcription initiation sites along with well-conserved TATA box, A-rich B recognition element (BRE), and promoter proximal element (PPE) motif in promoter regions, a high pyrimidine nucleotide content (T/C) at the -1 position, and Shine-Dalgarno (SD) motifs (GGDGRD) in 5' untranslated regions (5' UTRs). Along with differential transcript levels, 117 leaderless genes and 86 non-coding RNAs (ncRNAs) were identified, representing diverse cellular functions and potential regulatory functions under the different growth conditions. Interestingly, we observed low GC content in ncRNAs for RNA-based regulation via unstructured forms or interaction with other cellular components. Further comparative analysis of T. onnurineus upstream regulatory sequences with those of closely related archaeal genomes demonstrated that transcription of orthologous genes are initiated by highly conserved promoter sequences, however their upstream sequences for transcriptional and translational regulation are largely diverse. These results provide the genetic information of T. onnurineus for its future application in metabolic engineering.

Transcription Profiling of Bacillus subtilis Cells Infected with AR9, a Giant Phage Encoding Two Multisubunit RNA Polymerases

Bacteriophage AR9 is a recently sequenced jumbo phage that encodes two multisubunit RNA polymerases. Here we investigated the AR9 transcription strategy and the effect of AR9 infection on the transcription of its host, Bacillus subtilis. Analysis of whole-genome transcription revealed early, late, and continuously expressed AR9 genes. Alignment of sequences upstream of the 5′ ends of AR9 transcripts revealed consensus sequences that define early and late phage promoters. Continuously expressed AR9 genes have both early and late promoters in front of them. Early AR9 transcription is independent of protein synthesis and must be determined by virion RNA polymerase injected together with viral DNA. During infection, the overall amount of host mRNAs is significantly decreased. Analysis of relative amounts of host transcripts revealed notable differences in the levels of some mRNAs. The physiological significance of up- or downregulation of host genes for AR9 phage infection remains to be established. AR9 infection is significantly affected by rifampin, an inhibitor of host RNA polymerase transcription. The effect is likely caused by the antibiotic-induced killing of host cells, while phage genome transcription is solely performed by viral RNA polymerases.

Phages regulate the timing of the expression of their own genes to coordinate processes in the infected cell and maximize the release of viral progeny. Phages also alter the levels of host transcripts. Here we present the results of a temporal analysis of the host and viral transcriptomes of Bacillus subtilis infected with a giant phage, AR9. We identify viral promoters recognized by two virus-encoded RNA polymerases that are a unique feature of the phiKZ-related group of phages to which AR9 belongs. Our results set the stage for future analyses of highly unusual RNA polymerases encoded by AR9 and other phiKZ-related phages.

Effects of Kasugamycin on the Translatome of Escherichia coli

It is long known that Kasugamycin inhibits translation of canonical transcripts containing a 5’-UTR with a Shine Dalgarno (SD) motif, but not that of leaderless transcripts. To gain a global overview of the influence of Kasugamycin on translation efficiencies, the changes of the translatome of Escherichia coli induced by a 10 minutes Kasugamycin treatment were quantified. The effect of Kasugamycin differed widely, 102 transcripts were at least twofold more sensitive to Kasugamycin than average, and 137 transcripts were at least twofold more resistant, and there was a more than 100-fold difference between the most resistant and the most sensitive transcript. The 5’-ends of 19 transcripts were determined from treated and untreated cultures, but Kasugamycin resistance did neither correlate with the presence or absence of a SD motif, nor with differences in 5’-UTR lengths or GC content. RNA Structure Logos were generated for the 102 Kasugamycin-sensitive and for the 137 resistant transcripts. For both groups a short Shine Dalgarno (SD) motif was retrieved, but no specific motifs associated with resistance or sensitivity could be found. Notably, this was also true for the region -3 to -1 upstream of the start codon and the presence of an extended SD motif, which had been proposed to result in Kasugamycin resistance. Comparison of the translatome results with the database RegulonDB showed that the transcript with the highest resistance was leaderless, but no further leaderless transcripts were among the resistant transcripts. Unexpectedly, it was found that translational coupling might be a novel feature that is associated with Kasugamycin resistance. Taken together, Kasugamycin has a profound effect on translational efficiencies of E. coli transcripts, but the mechanism of action is different than previously described.

Negative Regulation of Ectoine Uptake and Catabolism in Sinorhizobium meliloti: Characterization of the EhuR Gene

Ectoine has osmoprotective effects on Sinorhizobium meliloti that differ from its effects in other bacteria. Ectoine does not accumulate in S. meliloti cells; instead, it is degraded. The products of the ehuABCD-eutABCDE operon were previously discovered to be responsible for the uptake and catabolism of ectoine in S. meliloti However, the mechanism by which ectoine is involved in the regulation of the ehuABCD-eutABCDE operon remains unclear. The ehuR gene, which is upstream of and oriented in the same direction as the ehuABCD-eutABCDE operon, encodes a member of the MocR/GntR family of transcriptional regulators. Quantitative reverse transcription-PCR and promoter-lacZ reporter fusion experiments revealed that EhuR represses transcription of the ehuABCD-eutABCDE operon, but this repression is inhibited in the presence of ectoine. Electrophoretic mobility shift assays and DNase I footprinting assays revealed that EhuR bound specifically to the DNA regions overlapping the -35 region of the ehuA promoter and the +1 region of the ehuR promoter. Surface plasmon resonance assays further demonstrated direct interactions between EhuR and the two promoters, although EhuR was found to have higher affinity for the ehuA promoter than for the ehuR promoter. In vitro, DNA binding by EhuR could be directly inhibited by a degradation product of ectoine. Our work demonstrates that EhuR is an important negative transcriptional regulator involved in the regulation of ectoine uptake and catabolism and is likely regulated by one or more end products of ectoine catabolism.

IMPORTANCE

Sinorhizobium meliloti is an important soil bacterium that displays symbiotic interactions with legume hosts. Ectoine serves as a key osmoprotectant for S. meliloti However, ectoine does not accumulate in the cells; rather, it is degraded. In this study, we characterized the transcriptional regulation of the operon responsible for ectoine uptake and catabolism in S. meliloti We identified and characterized the transcription repressor EhuR, which is the first MocR/GntR family member found to be involved in the regulation of compatible solute uptake and catabolism. More importantly, we demonstrated for the first time that an ectoine catabolic end product could modulate EhuR DNA-binding activity. Therefore, this work provides new insights into the unique mechanism of ectoine-induced osmoprotection in S. meliloti.

5'-Terminal AUGs in Escherichia coli mRNAs with Shine-Dalgarno Sequences: Identification and Analysis of Their Roles in Non-Canonical Translation Initiation

Analysis of the Escherichia coli transcriptome identified a unique subset of messenger RNAs (mRNAs) that contain a conventional untranslated leader and Shine-Dalgarno (SD) sequence upstream of the gene’s start codon while also containing an AUG triplet at the mRNA’s 5’- terminus (5’-uAUG). Fusion of the coding sequence specified by the 5’-terminal putative AUG start codon to a lacZ reporter gene, as well as primer extension inhibition assays, reveal that the majority of the 5’-terminal upstream open reading frames (5’-uORFs) tested support some level of lacZ translation, indicating that these mRNAs can function both as leaderless and canonical SD-leadered mRNAs. Although some of the uORFs were expressed at low levels, others were expressed at levels close to that of the respective downstream genes and as high as the naturally leaderless cI mRNA of bacteriophage λ. These 5’-terminal uORFs potentially encode peptides of varying lengths, but their functions, if any, are unknown. In an effort to determine whether expression from the 5’-terminal uORFs impact expression of the immediately downstream cistron, we examined expression from the downstream coding sequence after mutations were introduced that inhibit efficient 5’-uORF translation. These mutations were found to affect expression from the downstream cistrons to varying degrees, suggesting that some 5’-uORFs may play roles in downstream regulation. Since the 5’-uAUGs found on these conventionally leadered mRNAs can function to bind ribosomes and initiate translation, this indicates that canonical mRNAs containing 5’-uAUGs should be examined for their potential to function also as leaderless mRNAs.

Alternative Translation Initiation of a Haloarchaeal Serine Protease Transcript Containing Two In-Frame Start Codons

Recent studies have shown that haloarchaea employ leaderless and Shine-Dalgarno (SD)-less mechanisms for translation initiation of leaderless transcripts with a 5' untranslated region (5' UTR) of <10 nucleotides (nt) and leadered transcripts with a 5' UTR of ≥10 nt, respectively. However, whether the two mechanisms can operate on the same naturally occurring haloarchaeal transcript carrying multiple potential start codons is unknown. In this study, the transcript of the sptA gene (encoding an extracellular serine protease of Natrinema sp. strain J7-2) was experimentally determined and found to contain two potential in-frame AUG codons (AUG(1) and AUG(2)) located 5 and 29 nt, respectively, downstream of the transcription start site. Mutational analysis revealed that both AUGs can function as the translation start codon for production of active SptA, although AUG(1) is more efficient than AUG(2) for translation initiation. Insertion of a stable stem-loop structure between the two AUGs completely abolished initiation at AUG(1) but did not affect initiation at AUG(2), indicating that AUG(2)-initiated translation does not involve ribosome scanning from the 5' end of the transcript. Furthermore, the efficiency of AUG(2)-initiated translation was not influenced by an upstream SD-like sequence. In addition, both AUG(1) and AUG(2) contribute to transcript stability, probably by recruiting ribosomes to protect the transcript against degradation. These data suggest that depending on which of two in-frame start codons is used, the sptA transcript can act as either a leaderless or a leadered transcript for SptA production in haloarchaea.

IMPORTANCE

In eukaryotes and bacteria, alternative translation start sites contribute to proteome complexity and can be used as a functional mechanism to increase translation efficiency. However, little is known about alternative translation initiation in archaea. Our results demonstrate that leaderless and SD-less mechanisms can be used for translation initiation of the sptA transcript from two in-frame start codons, raising the possibility that in haloarchaea, alternative translation initiation on one transcript functions to increase translation efficiency and/or contribute to proteome complexity.

Genome wide transcription start sites analysis of Xanthomonas campestris pv. campestris B100 with insights into the gum gene cluster directing the biosynthesis of the exopolysaccharide xanthan

Xanthomonas campestris pv. campestris (Xcc) is the major producer of the exopolysaccharide xanthan, the commercially most important natural polysaccharide of microbial origin. The current work provides deeper insights into the yet uncharacterized transcriptomic features of the xanthan producing strain Xcc-B100. Towards this goal, RNA sequencing of a library based on the selective enrichment of the 5' ends of native transcripts was performed. This approach resulted in the genome wide identification of 3067 transcription start sites (TSSs) that were further classified based on their genomic positions. Among them, 1545 mapped upstream of an actively transcribed CDS and 1363 were classified as novel TSSs representing antisense, internal, and TSSs belonging to previously unidentified genomic features. Analyzing the transcriptional strength of primary and antisense TSSs revealed that in some instances antisense transcription seemed to be initiated at a higher level than its sense counterpart. Mapping the exact positions of TSSs aided in the identification of promoter consensus motifs, ribosomal binding sites, and enhanced the genome annotation of 159 in silico predicted translational start (TLS) sites. The global view on length distribution of the 5' untranslated regions (5'-UTRs) deduced from the data pointed to the occurrence of leaderless transcripts and transcripts with unusually long 5'-UTRs, in addition to identifying seven putative riboswitch elements for Xcc-B100. Concerning the biosynthesis of xanthan, we focused on the transcriptional organization of the gum gene cluster. Under the conditions tested, we present evidence for a complex transcription pattern of the gum genes with multiple TSSs and an obvious considerable role of antisense transcription. The gene gumB, encoding an outer membrane xanthan exporter, is presented here as an example for genes that possessed a strong antisense TSS.

Transcriptome landscape of Lactococcus lactis reveals many novel RNAs including a small regulatory RNA involved in carbon uptake and metabolism

RNA sequencing has revolutionized genome-wide transcriptome analyses, and the identification of non-coding regulatory RNAs in bacteria has thus increased concurrently. Here we reveal the transcriptome map of the lactic acid bacterial paradigm Lactococcus lactis MG1363 by employing differential RNA sequencing (dRNA-seq) and a combination of manual and automated transcriptome mining. This resulted in a high-resolution genome annotation of L. lactis and the identification of 60 cis-encoded antisense RNAs (asRNAs), 186 trans-encoded putative regulatory RNAs (sRNAs) and 134 novel small ORFs. Based on the putative targets of asRNAs, a novel classification is proposed. Several transcription factor DNA binding motifs were identified in the promoter sequences of (a)sRNAs, providing insight in the interplay between lactococcal regulatory RNAs and transcription factors. The presence and lengths of 14 putative sRNAs were experimentally confirmed by differential Northern hybridization, including the abundant RNA 6S that is differentially expressed depending on the available carbon source. For another sRNA, LLMGnc_147, functional analysis revealed that it is involved in carbon uptake and metabolism. L. lactis contains 13% leaderless mRNAs (lmRNAs) that, from an analysis of overrepresentation in GO classes, seem predominantly involved in nucleotide metabolism and DNA/RNA binding. Moreover, an A-rich sequence motif immediately following the start codon was uncovered, which could provide novel insight in the translation of lmRNAs. Altogether, this first experimental genome-wide assessment of the transcriptome landscape of L. lactis and subsequent sRNA studies provide an extensive basis for the investigation of regulatory RNAs in L. lactis and related lactococcal species.

Leaderless Transcripts and Small Proteins Are Common Features of the Mycobacterial Translational Landscape

RNA-seq technologies have provided significant insight into the transcription networks of mycobacteria. However, such studies provide no definitive information on the translational landscape. Here, we use a combination of high-throughput transcriptome and proteome-profiling approaches to more rigorously understand protein expression in two mycobacterial species. RNA-seq and ribosome profiling in Mycobacterium smegmatis, and transcription start site (TSS) mapping and N-terminal peptide mass spectrometry in Mycobacterium tuberculosis, provide complementary, empirical datasets to examine the congruence of transcription and translation in the Mycobacterium genus. We find that nearly one-quarter of mycobacterial transcripts are leaderless, lacking a 5’ untranslated region (UTR) and Shine-Dalgarno ribosome-binding site. Our data indicate that leaderless translation is a major feature of mycobacterial genomes and is comparably robust to leadered initiation. Using translational reporters to systematically probe the cis-sequence requirements of leaderless translation initiation in mycobacteria, we find that an ATG or GTG at the mRNA 5’ end is both necessary and sufficient. This criterion, together with our ribosome occupancy data, suggests that mycobacteria encode hundreds of small, unannotated proteins at the 5’ ends of transcripts. The conservation of small proteins in both mycobacterial species tested suggests that some play important roles in mycobacterial physiology. Our translational-reporter system further indicates that mycobacterial leadered translation initiation requires a Shine Dalgarno site in the 5’ UTR and that ATG, GTG, TTG, and ATT codons can robustly initiate translation. Our combined approaches provide the first comprehensive view of mycobacterial gene structures and their non-canonical mechanisms of protein expression.

The current paradigm for bacterial translation is based on an mRNA that includes an untranslated leader sequence containing the ribosome-binding site upstream of the initiation codon. We applied genome-scale approaches to map the protein-coding regions in the genomes of Mycobacterium smegmatis and Mycobacterium tuberculosis. We found that nearly one-quarter of mycobacterial transcripts are leaderless in mycobacterial species, thus indicating that ribosomes must recognize these mRNAs by a novel mechanism and suggesting that there are alternative modes of bacterial translation beyond the Escherichia coli paradigm. Our translational profiling showed that many mycobacterial proteins are mis-annotated, and also found many new genes encoding small proteins that had been previously overlooked, which are likely to play novel roles in diverse cellular processes. We also developed a new reporter system that provides mechanistic insights into translation initiation through deep sequencing. Our data show that leaderless translation is a robust process that is conserved in mycobacteria, that leaderless translation only requires that the mRNA begin with a start codon, and predict that mycobacteria encode hundreds of small proteins. This work will help us understand gene structure, genome organization and protein expression in bacteria, and how the translational machinery differs in different organisms.

The abundant and essential HU proteins in Deinococcus deserti and Deinococcus radiodurans are translated from leaderless mRNA

HU proteins have an important architectural role in nucleoid organization in bacteria. Compared with HU of many bacteria, HU proteins from Deinococcus species possess an N-terminal lysine-rich extension similar to the eukaryotic histone H1 C-terminal domain involved in DNA compaction. The single HU gene in Deinococcus radiodurans, encoding DrHU, is required for nucleoid compaction and cell viability. Deinococcus deserti contains three expressed HU genes, encoding DdHU1, DdHU2 and DdHU3. Here, we show that either DdHU1 or DdHU2 is essential in D. deserti. DdHU1 and DdHU2, but not DdHU3, can substitute for DrHU in D. radiodurans, indicating that DdHU3 may have a non-essential function different from DdHU1, DdHU2 and DrHU. Interestingly, the highly abundant DrHU and DdHU1 proteins, and also the less expressed DdHU2, are translated in Deinococcus from leaderless mRNAs, which lack a 5'-untranslated region and, hence, the Shine-Dalgarno sequence. Unexpectedly, cloning the DrHU or DdHU1 gene under control of a strong promoter in an expression plasmid, which results in leadered transcripts, strongly reduced the DrHU and DdHU1 protein level in D. radiodurans compared with that obtained from the natural leaderless gene. We also show that the start codon position for DrHU and DdHU1 should be reannotated, resulting in proteins that are 15 and 4 aa residues shorter than initially reported. The expression level and start codon correction were crucial for functional characterization of HU in Deinococcus.

Distribution and diversity of ribosome binding sites in prokaryotic genomes

BACKGROUND

Prokaryotic translation initiation involves the proper docking, anchoring, and accommodation of mRNA to the 30S ribosomal subunit. Three initiation factors (IF1, IF2, and IF3) and some ribosomal proteins mediate the assembly and activation of the translation initiation complex. Although the interaction between Shine-Dalgarno (SD) sequence and its complementary sequence in the 16S rRNA is important in initiation, some genes lacking an SD ribosome binding site (RBS) are still well expressed. The objective of this study is to examine the pattern of distribution and diversity of RBS in fully sequenced bacterial genomes. The following three hypotheses were tested: SD motifs are prevalent in bacterial genomes; all previously identified SD motifs are uniformly distributed across prokaryotes; and genes with specific cluster of orthologous gene (COG) functions differ in their use of SD motifs.

RESULTS

Data for 2,458 bacterial genomes, previously generated by Prodigal (PROkaryotic DYnamic programming Gene-finding ALgorithm) and currently available at the National Center for Biotechnology Information (NCBI), were analyzed. Of the total genes examined, ~77.0% use an SD RBS, while ~23.0% have no RBS. Majority of the genes with the most common SD motifs are distributed in a manner that is representative of their abundance for each COG functional category, while motifs 13 (5'-GGA-3'/5'-GAG-3'/5'-AGG-3') and 27 (5'-AGGAGG-3') appear to be predominantly used by genes for information storage and processing, and translation and ribosome biogenesis, respectively.

CONCLUSION

These findings suggest that an SD sequence is not obligatory for translation initiation; instead, other signals, such as the RBS spacer, may have an overarching influence on translation of mRNAs. Subsequent analyses of the 5' secondary structure of these mRNAs may provide further insight into the translation initiation mechanism.

Initiation of mRNA translation in bacteria: structural and dynamic aspects

Initiation of mRNA translation is a major checkpoint for regulating level and fidelity of protein synthesis. Being rate limiting in protein synthesis, translation initiation also represents the target of many post-transcriptional mechanisms regulating gene expression. The process begins with the formation of an unstable 30S pre-initiation complex (30S pre-IC) containing initiation factors (IFs) IF1, IF2 and IF3, the translation initiation region of an mRNA and initiator fMet-tRNA whose codon and anticodon pair in the P-site following a first-order rearrangement of the 30S pre-IC produces a locked 30S initiation complex (30SIC); this is docked by the 50S subunit to form a 70S complex that, following several conformational changes, positional readjustments of its ligands and ejection of the IFs, becomes a 70S initiation complex productive in initiation dipeptide formation. The first EF-G-dependent translocation marks the beginning of the elongation phase of translation. Here, we review structural, mechanistic and dynamical aspects of this process.

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

Background

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

Results

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

Conclusions

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

Electronic supplementary material

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

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

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

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

RNA sequencing and proteogenomics reveal the importance of leaderless mRNAs in the radiation-tolerant bacterium Deinococcus deserti

Deinococcus deserti is a desiccation- and radiation-tolerant desert bacterium. Differential RNA sequencing (RNA-seq) was performed to explore the specificities of its transcriptome. Strikingly, for 1,174 (60%) mRNAs, the transcription start site was found exactly at (916 cases, 47%) or very close to the translation initiation codon AUG or GUG. Such proportion of leaderless mRNAs, which may resemble ancestral mRNAs, is unprecedented for a bacterial species. Proteomics showed that leaderless mRNAs are efficiently translated in D. deserti. Interestingly, we also found 173 additional transcripts with a 5′-AUG or 5′-GUG that would make them competent for ribosome binding and translation into novel small polypeptides. Fourteen of these are predicted to be leader peptides involved in transcription attenuation. Another 30 correlated with new gene predictions and/or showed conservation with annotated and nonannotated genes in other Deinococcus species, and five of these novel polypeptides were indeed detected by mass spectrometry. The data also allowed reannotation of the start codon position of 257 genes, including several DNA repair genes. Moreover, several novel highly radiation-induced genes were found, and their potential roles are discussed. On the basis of our RNA-seq and proteogenomics data, we propose that translation of many of the novel leaderless transcripts, which may have resulted from single-nucleotide changes and maintained by selective pressure, provides a new explanation for the generation of a cellular pool of small peptides important for protection of proteins against oxidation and thus for radiation/desiccation tolerance and adaptation to harsh environmental conditions.

Transcriptome sequencing revealed the transcriptional organization at ribosome-mediated attenuation sites in Corynebacterium glutamicum and identified a novel attenuator involved in aromatic amino acid biosynthesis

The Gram-positive bacterium Corynebacterium glutamicum belongs to the order Corynebacteriales and is used as a producer of amino acids at industrial scales. Due to its economic importance, gene expression and particularly the regulation of amino acid biosynthesis has been investigated extensively. Applying the high-resolution technique of transcriptome sequencing (RNA-seq), recently a vast amount of data has been generated that was used to comprehensively analyze the C. glutamicum transcriptome. By analyzing RNA-seq data from a small RNA cDNA library of C. glutamicum, short transcripts in the known transcriptional attenuators sites of the trp operon, the ilvBNC operon and the leuA gene were verified. Furthermore, whole transcriptome RNA-seq data were used to elucidate the transcriptional organization of these three amino acid biosynthesis operons. In addition, we discovered and analyzed the novel attenuator aroR, located upstream of the aroF gene (cg1129). The DAHP synthase encoded by aroF catalyzes the first step in aromatic amino acid synthesis. The AroR leader peptide contains the amino acid sequence motif F-Y-F, indicating a regulatory effect by phenylalanine and tyrosine. Analysis by real-time RT-PCR suggests that the attenuator regulates the transcription of aroF in dependence of the cellular amount of tRNA loaded with phenylalanine when comparing a phenylalanine-auxotrophic C. glutamicum mutant fed with limiting and excess amounts of a phenylalanine-containing dipeptide. Additionally, the very interesting finding was made that all analyzed attenuators are leaderless transcripts.

Mitochondrial protein synthesis: figuring the fundamentals, complexities and complications, of mammalian mitochondrial translation

Mitochondrial protein synthesis is essential for all mammals, being responsible for providing key components of the oxidative phosphorylation complexes. Although only thirteen different polypeptides are made, the molecular details of this deceptively simple process remain incomplete. Central to this process is a non-canonical ribosome, the mitoribosome, which has evolved to address its unique mandate. In this review, we integrate the current understanding of the molecular aspects of mitochondrial translation with recent advances in structural biology. We identify numerous key questions that we will need to answer if we are to increase our knowledge of the molecular mechanisms underlying mitochondrial protein synthesis.

Comprehensive analysis of the Corynebacterium glutamicum transcriptome using an improved RNAseq technique

Background

The use of RNAseq to resolve the transcriptional organization of an organism was established in recent years and also showed the complexity and dynamics of bacterial transcriptomes. The aim of this study was to comprehensively investigate the transcriptome of the industrially relevant amino acid producer and model organism Corynebacterium glutamicum by RNAseq in order to improve its genome annotation and to describe important features for transcription and translation.

Results

RNAseq data sets were obtained by two methods, one that focuses on 5′-ends of primary transcripts and another that provides the overall transcriptome with an improved resolution of 3′-ends of transcripts. Subsequent data analysis led to the identification of more than 2,000 transcription start sites (TSSs), the definition of 5′-UTRs (untranslated regions) for annotated protein-coding genes, operon structures and many novel transcripts located between or in antisense orientation to protein-coding regions. Interestingly, a high number of mRNAs (33%) is transcribed as leaderless transcripts. From the data, consensus promoter and ribosome binding site (RBS) motifs were identified and it was shown that the majority of genes in C. glutamicum are transcribed monocistronically, but operons containing up to 16 genes are also present.

Conclusions

The comprehensive transcriptome map of C. glutamicum established in this study represents a major step forward towards a complete definition of genetic elements (e.g. promoter regions, gene starts and stops, 5′-UTRs, RBSs, transcript starts and ends) and provides the ideal basis for further analyses on transcriptional regulatory networks in this organism. The methods developed are easily applicable for other bacteria and have the potential to be used also for quantification of transcriptomes, replacing microarrays in the near future.

Genome-wide mapping of transcriptional start sites defines an extensive leaderless transcriptome in Mycobacterium tuberculosis

Deciphering physiological changes that mediate transition of Mycobacterium tuberculosis between replicating and nonreplicating states is essential to understanding how the pathogen can persist in an individual host for decades. We have combined RNA sequencing (RNA-seq) of 5′ triphosphate-enriched libraries with regular RNA-seq to characterize the architecture and expression of M. tuberculosis promoters. We identified over 4,000 transcriptional start sites (TSSs). Strikingly, for 26% of the genes with a primary TSS, the site of transcriptional initiation overlapped with the annotated start codon, generating leaderless transcripts lacking a 5′ UTR and, hence, the Shine-Dalgarno sequence commonly used to initiate ribosomal engagement in eubacteria. Genes encoding proteins with active growth functions were markedly depleted from the leaderless transcriptome, and there was a significant increase in the overall representation of leaderless mRNAs in a starvation model of growth arrest. The high percentage of leaderless genes may have particular importance in the physiology of nonreplicating M. tuberculosis.

•

A resource for the identification of in vitro active promoters in M. tuberculosis

•

A quarter of all genes in M. tuberculosis are expressed as leaderless mRNAs

•

Leaderless mRNAs are differentially associated with toxin-antitoxin modules

•

Abundance of leaderless mRNAs increases during starvation-induced growth arrest

A combination of improved differential and global RNA-seq reveals pervasive transcription initiation and events in all stages of the life-cycle of functional RNAs in Propionibacterium acnes, a major contributor to wide-spread human disease

BACKGROUND

Sequencing of the genome of Propionibacterium acnes produced a catalogue of genes many of which enable this organism to colonise skin and survive exposure to the elements. Despite this platform, there was little understanding of the gene regulation that gives rise to an organism that has a major impact on human health and wellbeing and causes infections beyond the skin. To address this situation, we have undertaken a genome-wide study of gene regulation using a combination of improved differential and global RNA-sequencing and an analytical approach that takes into account the inherent noise within the data.

RESULTS

We have produced nucleotide-resolution transcriptome maps that identify and differentiate sites of transcription initiation from sites of stable RNA processing and mRNA cleavage. Moreover, analysis of these maps provides strong evidence for 'pervasive' transcription and shows that contrary to initial indications it is not biased towards the production of antisense RNAs. In addition, the maps reveal an extensive array of riboswitches, leaderless mRNAs and small non-protein-coding RNAs alongside vegetative promoters and post-transcriptional events, which includes unusual tRNA processing. The identification of such features will inform models of complex gene regulation, as illustrated here for ribonucleotide reductases and a potential quorum-sensing, two-component system.

CONCLUSIONS

The approach described here, which is transferable to any bacterial species, has produced a step increase in whole-cell knowledge of gene regulation in P. acnes. Continued expansion of our maps to include transcription associated with different growth conditions and genetic backgrounds will provide a new platform from which to computationally model the gene expression that determines the physiology of P. acnes and its role in human disease.

Parallel evolution of genome structure and transcriptional landscape in the Epsilonproteobacteria

Background

Gene reshuffling, point mutations and horizontal gene transfer contribute to bacterial genome variation, but require the genome to rewire its transcriptional circuitry to ensure that inserted, mutated or reshuffled genes are transcribed at appropriate levels. The genomes of Epsilonproteobacteria display very low synteny, due to high levels of reshuffling and reorganisation of gene order, but still share a significant number of gene orthologs allowing comparison. Here we present the primary transcriptome of the pathogenic Epsilonproteobacterium Campylobacter jejuni, and have used this for comparative and predictive transcriptomics in the Epsilonproteobacteria.

Results

Differential RNA-sequencing using 454 sequencing technology was used to determine the primary transcriptome of C. jejuni NCTC 11168, which consists of 992 transcription start sites (TSS), which included 29 putative non-coding and stable RNAs, 266 intragenic (internal) TSS, and 206 antisense TSS. Several previously unknown features were identified in the C. jejuni transcriptional landscape, like leaderless mRNAs and potential leader peptides upstream of amino acid biosynthesis genes. A cross-species comparison of the primary transcriptomes of C. jejuni and the related Epsilonproteobacterium Helicobacter pylori highlighted a lack of conservation of operon organisation, position of intragenic and antisense promoters or leaderless mRNAs. Predictive comparisons using 40 other Epsilonproteobacterial genomes suggests that this lack of conservation of transcriptional features is common to all Epsilonproteobacterial genomes, and is associated with the absence of genome synteny in this subdivision of the Proteobacteria.

Conclusions

Both the genomes and transcriptomes of Epsilonproteobacteria are highly variable, both at the genome level by combining and division of multicistronic operons, but also on the gene level by generation or deletion of promoter sequences and 5′ untranslated regions. Regulatory features may have evolved after these species split from a common ancestor, with transcriptome rewiring compensating for changes introduced by genomic reshuffling and horizontal gene transfer.

Gene expression control by selective RNA processing and stabilization in bacteria

RNA maturation is a key event regulating genes at post-transcriptional level. In bacteria, it is employed to adjust the amounts of proteins and functional RNAs, often in response to environmental constraints. During the process of RNA maturation, enzymes and factors that would otherwise promote RNA degradation convert a labile RNA into a stable and biologically functional molecule.

Global mapping of transcription start sites and promoter motifs in the symbiotic α-proteobacterium Sinorhizobium meliloti 1021

BACKGROUND

Sinorhizobium meliloti is a soil-dwelling α-proteobacterium that possesses a large, tripartite genome and engages in a nitrogen fixing symbiosis with its plant hosts. Although much is known about this important model organism, global characterization of genetic regulatory circuits has been hampered by a lack of information about transcription and promoters.

RESULTS

Using an RNAseq approach and RNA populations representing 16 different growth and stress conditions, we comprehensively mapped S. meliloti transcription start sites (TSS). Our work identified 17,001 TSS that we grouped into six categories based on the genomic context of their transcripts: mRNA (4,430 TSS assigned to 2,657 protein-coding genes), leaderless mRNAs (171), putative mRNAs (425), internal sense transcripts (7,650), antisense RNA (3,720), and trans-encoded sRNAs (605). We used this TSS information to identify transcription factor binding sites and putative promoter sequences recognized by seven of the 15 known S. meliloti σ factors σ70, σ54, σH1, σH2, σE1, σE2, and σE9). Altogether, we predicted 2,770 new promoter sequences, including 1,302 located upstream of protein coding genes and 722 located upstream of antisense RNA or trans-encoded sRNA genes. To validate promoter predictions for targets of the general stress response σ factor, RpoE2 (σE2), we identified rpoE2-dependent genes using microarrays and confirmed TSS for a subset of these by 5' RACE mapping.

CONCLUSIONS

By identifying TSS and promoters on a global scale, our work provides a firm foundation for the continued study of S. meliloti gene expression with relation to gene organization, σ factors and other transcription factors, and regulatory RNAs.

Rhamnosyltransferase genes migA and wapR are regulated in a differential manner to modulate the quantities of core oligosaccharide glycoforms produced by Pseudomonas aeruginosa

migA and wapR are rhamnosyltransferase genes involved in the biosynthesis of Pseudomonas aeruginosa lipopolysaccharide core oligosaccharide. Here, we show that preferential expression of migA and wapR correlated with the levels of uncapped and O polysaccharide-capped core, respectively. wapR is negatively regulated, while migA is positively regulated by RhlR/RhlI quorum sensing.

Selective detection of peptide-oligonucleotide heteroconjugates utilizing capillary HPLC-ICPMS

A method for the selective detection and quantification of peptide:oligonucleotide heteroconjugates, such as those generated by protein:nucleic acid cross-links, using capillary reversed-phase high performance liquid chromatography (cap-RPHPLC) coupled with inductively coupled plasma mass spectrometry detection (ICPMS) is described. The selective detection of phosphorus as (31)P(+), the only natural isotope, in peptide-oligonucleotide heteroconjugates is enabled by the elemental detection capabilities of the ICPMS. Mobile phase conditions that allow separation of heteroconjugates while maintaining ICPMS compatibility were investigated. We found that trifluoroacetic acid (TFA) mobile phases, used in conventional peptide separations, and hexafluoroisopropanol/triethylamine (HFIP/TEA) mobile phases, used in conventional oligonucleotide separations, both are compatible with ICPMS and enable heteroconjugate separation. The TFA-based separations yielded limits of detection (LOD) of ~40 ppb phosphorus, which is nearly seven times lower than the LOD for HFIP/TEA-based separations. Using the TFA mobile phase, 1-2 pmol of a model heteroconjugate were routinely separated and detected by this optimized capLC-ICPMS method.

The transcriptional landscape and small RNAs of Salmonella enterica serovar Typhimurium

More than 50 y of research have provided great insight into the physiology, metabolism, and molecular biology of Salmonella enterica serovar Typhimurium (S. Typhimurium), but important gaps in our knowledge remain. It is clear that a precise choreography of gene expression is required for Salmonella infection, but basic genetic information such as the global locations of transcription start sites (TSSs) has been lacking. We combined three RNA-sequencing techniques and two sequencing platforms to generate a robust picture of transcription in S. Typhimurium. Differential RNA sequencing identified 1,873 TSSs on the chromosome of S. Typhimurium SL1344 and 13% of these TSSs initiated antisense transcripts. Unique findings include the TSSs of the virulence regulators phoP, slyA, and invF. Chromatin immunoprecipitation revealed that RNA polymerase was bound to 70% of the TSSs, and two-thirds of these TSSs were associated with σ(70) (including phoP, slyA, and invF) from which we identified the -10 and -35 motifs of σ(70)-dependent S. Typhimurium gene promoters. Overall, we corrected the location of important genes and discovered 18 times more promoters than identified previously. S. Typhimurium expresses 140 small regulatory RNAs (sRNAs) at early stationary phase, including 60 newly identified sRNAs. Almost half of the experimentally verified sRNAs were found to be unique to the Salmonella genus, and <20% were found throughout the Enterobacteriaceae. This description of the transcriptional map of SL1344 advances our understanding of S. Typhimurium, arguably the most important bacterial infection model.

A 5'-terminal phosphate is required for stable ternary complex formation and translation of leaderless mRNA in Escherichia coli

The bacteriophage λ's cI mRNA was utilized to examine the importance of the 5'-terminal phosphate on expression of leadered and leaderless mRNA in Escherichia coli. A hammerhead ribozyme was used to produce leadered and leaderless mRNAs, in vivo and in vitro, that contain a 5'-hydroxyl. Although these mRNAs may not occur naturally in the bacterial cell, they allow for the study of the importance of the 5'-phosphorylation state in ribosome binding and translation of leadered and leaderless mRNAs. Analyses with mRNAs containing either a 5'-phosphate or a 5'-hydroxyl indicate that leaderless cI mRNA requires a 5'-phosphate for stable ribosome binding in vitro as well as expression in vivo. Ribosome-binding assays show that 30S subunits and 70S ribosomes do not bind as strongly to 5'-hydroxyl as they do to 5'-phosphate containing leaderless mRNA and the tRNA-dependent ternary complex is less stable. Additionally, filter-binding assays revealed that the 70S ternary complex formed with a leaderless mRNA containing a 5'-hydroxyl has a dissociation rate (k(off)) that is 4.5-fold higher compared with the complex formed with a 5'-phosphate leaderless mRNA. Fusion to a lacZ reporter gene revealed that leaderless cI mRNA expression with a 5'-hydroxyl was >100-fold lower than the equivalent mRNA with a 5'-phosphate. These data indicate that a 5'-phosphate is an important feature of leaderless mRNA for stable ribosome binding and expression.

The architecture and ppGpp-dependent expression of the primary transcriptome of Salmonella Typhimurium during invasion gene expression

Background

Invasion of intestinal epithelial cells by Salmonella enterica serovar Typhimurium (S. Typhimurium) requires expression of the extracellular virulence gene expression programme (ST^EX), activation of which is dependent on the signalling molecule guanosine tetraphosphate (ppGpp). Recently, next-generation transcriptomics (RNA-seq) has revealed the unexpected complexity of bacterial transcriptomes and in this report we use differential RNA sequencing (dRNA-seq) to define the high-resolution transcriptomic architecture of wild-type S. Typhimurium and a ppGpp null strain under growth conditions which model ST^EX. In doing so we show that ppGpp plays a much wider role in regulating the S. Typhimurium ST^EX primary transcriptome than previously recognised.

Results

Here we report the precise mapping of transcriptional start sites (TSSs) for 78% of the S. Typhimurium open reading frames (ORFs). The TSS mapping enabled a genome-wide promoter analysis resulting in the prediction of 169 alternative sigma factor binding sites, and the prediction of the structure of 625 operons. We also report the discovery of 55 new candidate small RNAs (sRNAs) and 302 candidate antisense RNAs (asRNAs). We discovered 32 ppGpp-dependent alternative TSSs and determined the extent and level of ppGpp-dependent coding and non-coding transcription. We found that 34% and 20% of coding and non-coding RNA transcription respectively was ppGpp-dependent under these growth conditions, adding a further dimension to the role of this remarkable small regulatory molecule in enabling rapid adaptation to the infective environment.

Conclusions

The transcriptional architecture of S. Typhimurium and finer definition of the key role ppGpp plays in regulating Salmonella coding and non-coding transcription should promote the understanding of gene regulation in this important food borne pathogen and act as a resource for future research.

Translation Initiation

Selection of correct start codons on messenger RNAs is a key step required for faithful translation of the genetic message. Such a selection occurs in a complex process, during which a translation-competent ribosome assembles, eventually having in its P site a specialized methionyl-tRNAMet base-paired with the start codon on the mRNA. This chapter summarizes recent advances describing at the molecular level the successive steps involved in the process. Special emphasis is put on the roles of the three initiation factors and of the initiator tRNA, which are crucial for the efficiency and the specificity of the process. In particular, structural analyses concerning complexes containing ribosomal subunits, as well as detailed kinetic studies, have shed new light on the sequence of events leading to faithful initiation of protein synthesis in Bacteria.

Genome-wide transcriptome analysis of the plant pathogen Xanthomonas identifies sRNAs with putative virulence functions

The Gram-negative plant-pathogenic bacterium Xanthomonas campestris pv. vesicatoria (Xcv) is an important model to elucidate the mechanisms involved in the interaction with the host. To gain insight into the transcriptome of the Xcv strain 85–10, we took a differential RNA sequencing (dRNA-seq) approach. Using a novel method to automatically generate comprehensive transcription start site (TSS) maps we report 1421 putative TSSs in the Xcv genome. Genes in Xcv exhibit a poorly conserved −10 promoter element and no consensus Shine-Dalgarno sequence. Moreover, 14% of all mRNAs are leaderless and 13% of them have unusually long 5′-UTRs. Northern blot analyses confirmed 16 intergenic small RNAs and seven cis-encoded antisense RNAs in Xcv. Expression of eight intergenic transcripts was controlled by HrpG and HrpX, key regulators of the Xcv type III secretion system. More detailed characterization identified sX12 as a small RNA that controls virulence of Xcv by affecting the interaction of the pathogen and its host plants. The transcriptional landscape of Xcv is unexpectedly complex, featuring abundant antisense transcripts, alternative TSSs and clade-specific small RNAs.

Mycobacterium tuberculosis H37Rv sigC is expressed from two promoters but is not auto-regulatory

The extracytoplasmic function (ECF) sigma factor SigC has been implicated in the pathogenesis of Mycobacterium tuberculosis but control of its expression and activity is poorly understood. No proteins that interact with SigC have been detected leading to the suggestion that this sigma factor may be primarily controlled at the level of transcription. It has been suggested that SigC may be autoregulatory and a role has also been proposed for SigF in the expression of sigC. In this study we identified two promoters that were active under standard growth conditions by a combination of transcript start site mapping and promoter-lacZ fusion assays. The dominant promoter, P1, closely resembled mycobacterial SigA-dependent promoters, and introduction of a single base change at the conserved A of the -10 region eliminated promoter activity. Although the sequence of the other, P2, closely resembled the reported SigC consensus motifs, expression directed by this promoter was unaltered in a ΔsigC mutant strain, or in strains defective in other ECF sigma factors for which some similarity in consensus sequences was apparent. Comparison of the effects of different changes in the -10 region suggested that the P2 promoter was most likely recognised by SigA.

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

Over the past 3 years, bacterial transcriptomics has undergone a massive revolution. Increased sequencing capacity and novel tools have made it possible to explore the bacterial transcriptome to an unprecedented depth, which has revealed that the transcriptome is more complex and dynamic than expected. Alternative transcripts within operons challenge the classic operon definition, and many small RNAs involved in the regulation of transcription, translation and pathogenesis have been discovered. Furthermore, mRNAs may localize to specific areas in the cell, and the spatial organization and dynamics of the chromosome have been shown to be important for transcription. Epigenetic modifications of DNA also affect transcription, and RNA processing affects translation. Therefore, transcription in bacteria resembles that in eukaryotes in terms of complexity more closely than was previously thought. Here we will discuss the contribution of 'omics' approaches to these discoveries as well as the possible impact that they are expected to have in the future.

Identification of prokaryotic small proteins using a comparative genomic approach

MOTIVATION

Accurate prediction of genes encoding small proteins (on the order of 50 amino acids or less) remains an elusive open problem in bioinformatics. Some of the best methods for gene prediction use either sequence composition analysis or sequence similarity to a known protein coding sequence. These methods often fail for small proteins, however, either due to a lack of experimentally verified small protein coding genes or due to the limited statistical significance of statistics on small sequences. Our approach is based upon the hypothesis that true small proteins will be under selective pressure for encoding the particular amino acid sequence, for ease of translation by the ribosome and for structural stability. This stability can be achieved either independently or as part of a larger protein complex. Given this assumption, it follows that small proteins should display conserved local protein structure properties much like larger proteins. Our method incorporates neural-net predictions for three local structure alphabets within a comparative genomic approach using a genomic alignment of 22 closely related bacteria genomes to generate predictions for whether or not a given open reading frame (ORF) encodes for a small protein.

RESULTS

We have applied this method to the complete genome for Escherichia coli strain K12 and looked at how well our method performed on a set of 60 experimentally verified small proteins from this organism. Out of a total of 11 407 possible ORFs, we found that 6 of the top 10 and 27 of the top 100 predictions belonged to the set of 60 experimentally verified small proteins. We found 35 of all the true small proteins within the top 200 predictions. We compared our method to Glimmer, using a default Glimmer protocol and a modified small ORF Glimmer protocol with a lower minimum size cutoff. The default Glimmer protocol identified 16 of the true small proteins (all in the top 200 predictions), but failed to predict on 34 due to size cutoffs. The small ORF Glimmer protocol made predictions for all the experimentally verified small proteins but only contained 9 of the 60 true small proteins within the top 200 predictions.

CONTACT

jsamayoa@jhu.edu

Enhanced expression of recX in Mycobacterium tuberculosis owing to a promoter internal to recA

RecX is a small protein that interacts with, and modulates the activity of, RecA protein. In mycobacteria the recX gene is located immediately downstream of the recA gene, and the coding regions overlap. It has previously been shown that these two genes are co-transcribed in Mycobacterium smegmatis. In this study we examine the expression of recX in Mycobacterium tuberculosis. In addition to being co-transcribed with recA from the DNA-damage inducible recA promoters, we identify a constitutive recX promoter located within the recA coding sequence that is strong enough to make a significant contribution to the expression level of recX in the absence of DNA damage. Intriguingly, this promoter is inactivated in M. smegmatis by a critical base change in the -10 promoter motif, which probably accounts for the lower level of expression of recX relative to recA that we observed in that species. It is possible that this difference in relative expression influences RecA functions including the response to DNA damage of LexA-regulated genes.

Proximity of the start codon to a leaderless mRNA's 5' terminus is a strong positive determinant of ribosome binding and expression in Escherichia coli

An AUG start codon is an important determinant of ribosome binding and expression of leaderless mRNAs in Escherichia coli. Using reporter constructs encoding mRNAs where the AUG start codon is preceded by untranslated leaders of various length and sequence, we find that close proximity of the start codon to the 5' terminus and the leader sequence are strong determinants of both ribosome binding and expression.

Preferential selection of the 5'-terminal start codon on leaderless mRNAs by mammalian mitochondrial ribosomes

Mammalian mitochondrial mRNAs are basically leaderless, having few or no untranslated nucleotides prior to the 5'-start codon. We demonstrate here that mammalian mitochondrial 55 S ribosomes preferentially form initiation complexes at a 5'-terminal AUG codon over an internal AUG. The preferential use of the 5'-start codon is also seen on mitochondrial 28 S small subunits, which suggests that mitochondrial translation initiation on leaderless mRNAs does not require the large ribosomal subunit. The selection of the 5'-AUG is dependent on the presence of fMet-tRNA and is enhanced by the presence of the mitochondrial initiation factor IF2(mt). In prokaryotes, IF3 is believed to antagonize initiation on leaderless mRNAs. However, IF3(mt) stimulates initiation complex formation on leaderless mRNAs when tested with 55 S ribosomes. The addition of even a few nucleotides 5' to the AUG codon significantly reduces the efficiency of initiation, highlighting the importance of the leaderless or nearly leaderless nature of mitochondrial mRNAs. In addition, very few initiation complexes could form on a hybrid mRNA construct consisting of tRNA(Met) attached at the 5'-end of a mitochondrial protein-coding sequence. This observation demonstrates that post-transcriptional processing must occur prior to translation in mammalian mitochondria.