Leaderless genes in bacteria: clue to the evolution of translation initiation mechanisms in prokaryotes

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.

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

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.

A Novel Quality Measure and Correction Procedure for the Annotation of Microbial Translation Initiation Sites

The identification of translation initiation sites (TISs) constitutes an important aspect of sequence-based genome analysis. An erroneous TIS annotation can impair the identification of regulatory elements and N-terminal signal peptides, and also may flaw the determination of descent, for any particular gene. We have formulated a reference-free method to score the TIS annotation quality. The method is based on a comparison of the observed and expected distribution of all TISs in a particular genome given prior gene-calling. We have assessed the TIS annotations for all available NCBI RefSeq microbial genomes and found that approximately 87% is of appropriate quality, whereas 13% needs substantial improvement. We have analyzed a number of factors that could affect TIS annotation quality such as GC-content, taxonomy, the fraction of genes with a Shine-Dalgarno sequence and the year of publication. The analysis showed that only the first factor has a clear effect. We have then formulated a straightforward Principle Component Analysis-based TIS identification strategy to self-organize and score potential TISs. The strategy is independent of reference data and a priori calculations. A representative set of 277 genomes was subjected to the analysis and we found a clear increase in TIS annotation quality for the genomes with a low quality score. The PCA-based annotation was also compared with annotation with the current tool of reference, Prodigal. The comparison for the model genome of Escherichia coli K12 showed that both methods supplement each other and that prediction agreement can be used as an indicator of a correct TIS annotation. Importantly, the data suggest that the addition of a PCA-based strategy to a Prodigal prediction can be used to ‘flag’ TIS annotations for re-evaluation and in addition can be used to evaluate a given annotation in case a Prodigal annotation is lacking.

Increasing Avermectin Production in Streptomyces avermitilis by Manipulating the Expression of a Novel TetR-Family Regulator and Its Target Gene Product

Avermectins produced by Streptomyces avermitilis are commercially important anthelmintic agents. The detailed regulatory mechanisms of avermectin biosynthesis remain unclear. Here, we identified SAV3619, a TetR-family transcriptional regulator designated AveT, to be an activator for both avermectin production and morphological differentiation in S. avermitilis. AveT was shown to indirectly stimulate avermectin production by affecting transcription of the cluster-situated activator gene aveR. AveT directly repressed transcription of its own gene (aveT), adjacent gene pepD2 (sav_3620), sav_7490 (designated aveM), and sav_7491 by binding to an 18-bp perfect palindromic sequence (CGAAACGKTKYCGTTTCG, where K is T or G and Y is T or C and where the underlining indicates inverted repeats) within their promoter regions. aveM (which encodes a putative transmembrane efflux protein belonging to the major facilitator superfamily [MFS]), the important target gene of AveT, had a striking negative effect on avermectin production and morphological differentiation. Overexpression of aveT and deletion of aveM in wild-type and industrial strains of S. avermitilis led to clear increases in the levels of avermectin production. In vitro gel-shift assays suggested that C-5-O-B1, the late pathway precursor of avermectin B1, acts as an AveT ligand. Taken together, our findings indicate positive-feedback regulation of aveT expression and avermectin production by a late pathway intermediate and provide the basis for an efficient strategy to increase avermectin production in S. avermitilis by manipulation of AveT and its target gene product, AveM.

Transcriptome analysis of thermophilic methylotrophic Bacillus methanolicus MGA3 using RNA-sequencing provides detailed insights into its previously uncharted transcriptional landscape

Background

Bacillus methanolicus MGA3 is a thermophilic, facultative ribulose monophosphate (RuMP) cycle methylotroph. Together with its ability to produce high yields of amino acids, the relevance of this microorganism as a promising candidate for biotechnological applications is evident. The B. methanolicus MGA3 genome consists of a 3,337,035 nucleotides (nt) circular chromosome, the 19,174 nt plasmid pBM19 and the 68,999 nt plasmid pBM69. 3,218 protein-coding regions were annotated on the chromosome, 22 on pBM19 and 82 on pBM69. In the present study, the RNA-seq approach was used to comprehensively investigate the transcriptome of B. methanolicus MGA3 in order to improve the genome annotation, identify novel transcripts, analyze conserved sequence motifs involved in gene expression and reveal operon structures. For this aim, two different cDNA library preparation methods were applied: one which allows characterization of the whole transcriptome and another which includes enrichment of primary transcript 5′-ends.

Results

Analysis of the primary transcriptome data enabled the detection of 2,167 putative transcription start sites (TSSs) which were categorized into 1,642 TSSs located in the upstream region (5′-UTR) of known protein-coding genes and 525 TSSs of novel antisense, intragenic, or intergenic transcripts. Firstly, 14 wrongly annotated translation start sites (TLSs) were corrected based on primary transcriptome data. Further investigation of the identified 5′-UTRs resulted in the detailed characterization of their length distribution and the detection of 75 hitherto unknown cis-regulatory RNA elements. Moreover, the exact TSSs positions were utilized to define conserved sequence motifs for translation start sites, ribosome binding sites and promoters in B. methanolicus MGA3. Based on the whole transcriptome data set, novel transcripts, operon structures and mRNA abundances were determined. The analysis of the operon structures revealed that almost half of the genes are transcribed monocistronically (940), whereas 1,164 genes are organized in 381 operons. Several of the genes related to methylotrophy had highly abundant transcripts.

Conclusion

The extensive insights into the transcriptional landscape of B. methanolicus MGA3, gained in this study, represent a valuable foundation for further comparative quantitative transcriptome analyses and possibly also for the development of molecular biology tools which at present are very limited for this organism.

Electronic supplementary material

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

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.

Complex and extensive post-transcriptional regulation revealed by integrative proteomic and transcriptomic analysis of metabolite stress response in Clostridium acetobutylicum

BACKGROUND

Clostridium acetobutylicum is a model organism for both clostridial biology and solvent production. The organism is exposed to its own toxic metabolites butyrate and butanol, which trigger an adaptive stress response. Integrative analysis of proteomic and RNAseq data may provide novel insights into post-transcriptional regulation.

RESULTS

The identified iTRAQ-based quantitative stress proteome is made up of 616 proteins with a 15 % genome coverage. The differentially expressed proteome correlated poorly with the corresponding differential RNAseq transcriptome. Up to 31 % of the differentially expressed proteins under stress displayed patterns opposite to those of the transcriptome, thus suggesting significant post-transcriptional regulation. The differential proteome of the translation machinery suggests that cells employ a different subset of ribosomal proteins under stress. Several highly upregulated proteins but with low mRNA levels possessed mRNAs with long 5'UTRs and strong RBS scores, thus supporting the argument that regulatory elements on the long 5'UTRs control their translation. For example, the oxidative stress response rubrerythrin was upregulated only at the protein level up to 40-fold without significant mRNA changes. We also identified many leaderless transcripts, several displaying different transcriptional start sites, thus suggesting mRNA-trimming mechanisms under stress. Downregulation of Rho and partner proteins pointed to changes in transcriptional elongation and termination under stress.

CONCLUSIONS

The integrative proteomic-transcriptomic analysis demonstrated complex expression patterns of a large fraction of the proteome. Such patterns could not have been detected with one or the other omic analyses. Our analysis proposes the involvement of specific molecular mechanisms of post-transcriptional regulation to explain the observed complex stress response.

The coding and noncoding architecture of the Caulobacter crescentus genome

Caulobacter crescentus undergoes an asymmetric cell division controlled by a genetic circuit that cycles in space and time. We provide a universal strategy for defining the coding potential of bacterial genomes by applying ribosome profiling, RNA-seq, global 5′-RACE, and liquid chromatography coupled with tandem mass spectrometry (LC-MS) data to the 4-megabase C. crescentus genome. We mapped transcript units at single base-pair resolution using RNA-seq together with global 5′-RACE. Additionally, using ribosome profiling and LC-MS, we mapped translation start sites and coding regions with near complete coverage. We found most start codons lacked corresponding Shine-Dalgarno sites although ribosomes were observed to pause at internal Shine-Dalgarno sites within the coding DNA sequence (CDS). These data suggest a more prevalent use of the Shine-Dalgarno sequence for ribosome pausing rather than translation initiation in C. crescentus. Overall 19% of the transcribed and translated genomic elements were newly identified or significantly improved by this approach, providing a valuable genomic resource to elucidate the complete C. crescentus genetic circuitry that controls asymmetric cell division.

Caulobacter crescentus is a model system for studying asymmetric cell division, a fundamental process that, through differential gene expression in the two daughter cells, enables the generation of cells with different fates. To explore how the genome directs and maintains asymmetry upon cell division, we performed a coordinated analysis of multiple genomic and proteomic datasets to identify the RNA and protein coding features in the C. crescentus genome. Our integrated analysis identifies many new genetic regulatory elements, adding significant regulatory complexity to the C. crescentus genome. Surprisingly, 75.4% of protein coding genes lack a canonical translation initiation sequence motif (the Shine-Dalgarno site) which hybridizes to the 3′ end of the ribosomal RNA allowing translation initiation. We find Shine-Dalgarno sites primarily inside of genes where they cause translating ribosomes to pause, possibly allowing nascent proteins to correctly fold. With our detailed map of genomic transcription and translation elements, a systems view of the genetic network that controls asymmetric cell division is within reach.

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.

The adnAB locus, encoding a putative helicase-nuclease activity, is essential in Streptomyces

Homologous recombination is a crucial mechanism that repairs a wide range of DNA lesions, including the most deleterious ones, double-strand breaks (DSBs). This multistep process is initiated by the resection of the broken DNA ends by a multisubunit helicase-nuclease complex exemplified by Escherichia coli RecBCD, Bacillus subtilis AddAB, and newly discovered Mycobacterium tuberculosis AdnAB. Here we show that in Streptomyces, neither recBCD nor addAB homologues could be detected. The only putative helicase-nuclease-encoding genes identified were homologous to M. tuberculosis adnAB genes. These genes are conserved as a single copy in all sequenced genomes of Streptomyces. The disruption of adnAB in Streptomyces ambofaciens and Streptomyces coelicolor could not be achieved unless an ectopic copy was provided, indicating that adnAB is essential for growth. Both adnA and adnB genes were shown to be inducible in response to DNA damage (mitomycin C) and to be independently transcribed. Introduction of S. ambofaciens adnAB genes in an E. coli recB mutant restored viability and resistance to UV light, suggesting that Streptomyces AdnAB could be a functional homologue of RecBCD and be involved in DNA damage resistance.

Predicting the minimal translation apparatus: lessons from the reductive evolution of mollicutes

Mollicutes is a class of parasitic bacteria that have evolved from a common Firmicutes ancestor mostly by massive genome reduction. With genomes under 1 Mbp in size, most Mollicutes species retain the capacity to replicate and grow autonomously. The major goal of this work was to identify the minimal set of proteins that can sustain ribosome biogenesis and translation of the genetic code in these bacteria. Using the experimentally validated genes from the model bacteria Escherichia coli and Bacillus subtilis as input, genes encoding proteins of the core translation machinery were predicted in 39 distinct Mollicutes species, 33 of which are culturable. The set of 260 input genes encodes proteins involved in ribosome biogenesis, tRNA maturation and aminoacylation, as well as proteins cofactors required for mRNA translation and RNA decay. A core set of 104 of these proteins is found in all species analyzed. Genes encoding proteins involved in post-translational modifications of ribosomal proteins and translation cofactors, post-transcriptional modifications of t+rRNA, in ribosome assembly and RNA degradation are the most frequently lost. As expected, genes coding for aminoacyl-tRNA synthetases, ribosomal proteins and initiation, elongation and termination factors are the most persistent (i.e. conserved in a majority of genomes). Enzymes introducing nucleotides modifications in the anticodon loop of tRNA, in helix 44 of 16S rRNA and in helices 69 and 80 of 23S rRNA, all essential for decoding and facilitating peptidyl transfer, are maintained in all species. Reconstruction of genome evolution in Mollicutes revealed that, beside many gene losses, occasional gains by horizontal gene transfer also occurred. This analysis not only showed that slightly different solutions for preserving a functional, albeit minimal, protein synthetizing machinery have emerged in these successive rounds of reductive evolution but also has broad implications in guiding the reconstruction of a minimal cell by synthetic biology approaches.

Mapping of genotype-phenotype diversity among clinical isolates of mycobacterium tuberculosis by sequence-based transcriptional profiling

Genome sequencing has identified an extensive repertoire of single nucleotide polymorphisms among clinical isolates of Mycobacterium tuberculosis, but the extent to which these differences influence phenotypic properties of the bacteria remains to be elucidated. To determine whether these polymorphisms give rise to phenotypic diversity, we have integrated genome data sets with RNA sequencing to assess their impact on the comparative transcriptome profiles of strains belonging to M. tuberculosis Lineages 1 and 2. We observed clear correlations between genotype and transcriptional phenotype. These arose by three mechanisms. First, lineage-specific changes in amino acid sequence of transcriptional regulators were associated with alterations in their ability to control gene expression. Second, changes in nucleotide sequence were associated with alteration of promoter activity and generation of novel transcriptional start sites in intergenic regions and within coding sequences. We show that in some cases this mechanism is expected to generate functionally active truncated proteins involved in innate immune recognition. Finally, genes showing lineage-specific patterns of differential expression not linked directly to primary mutations were characterized by a striking overrepresentation of toxin–antitoxin pairs. Taken together, these findings advance our understanding of mycobacterial evolution, contribute to a systems level understanding of this important human pathogen, and more broadly demonstrate the application of state-of-the-art techniques to provide novel insight into mechanisms by which intergenic and silent mutations contribute to diversity.

Comparative transcriptomics of two environmentally relevant cyanobacteria reveals unexpected transcriptome diversity

Prochlorococcus is a genus of abundant and ecologically important marine cyanobacteria. Here, we present a comprehensive comparison of the structure and composition of the transcriptomes of two Prochlorococcus strains, which, despite their similarities, have adapted their gene pool to specific environmental constraints. We present genome-wide maps of transcriptional start sites (TSS) for both organisms, which are representatives of the two most diverse clades within the two major ecotypes adapted to high- and low-light conditions, respectively. Our data suggest antisense transcription for three-quarters of all genes, which is substantially more than that observed in other bacteria. We discovered hundreds of TSS within genes, most notably within 16 of the 29 prochlorosin genes, in strain MIT9313. A direct comparison revealed very little conservation in the location of TSS and the nature of non-coding transcripts between both strains. We detected extremely short 5' untranslated regions with a median length of only 27 and 29 nt for MED4 and MIT9313, respectively, and for 8% of all protein-coding genes the median distance to the start codon is only 10 nt or even shorter. These findings and the absence of an obvious Shine-Dalgarno motif suggest that leaderless translation and ribosomal protein S1-dependent translation constitute alternative mechanisms for translation initiation in Prochlorococcus. We conclude that genome-wide antisense transcription is a major component of the transcriptional output from these relatively small genomes and that a hitherto unrecognized high degree of complexity and variability of gene expression exists in their transcriptional architecture.

Haloferax volcanii, a prokaryotic species that does not use the Shine Dalgarno mechanism for translation initiation at 5'-UTRs

It was long assumed that translation initiation in prokaryotes generally occurs via the so-called Shine Dalgarno (SD) mechanism. Recently, it became clear that translation initiation in prokaryotes is more heterogeneous. In the haloarchaeon Haloferax volcanii, the majority of transcripts is leaderless and most transcripts with a 5′-UTR lack a SD motif. Nevertheless, a bioinformatic analysis predicted that 20–30% of all genes are preceded by a SD motif in haloarchaea. To analyze the importance of the SD mechanism for translation initiation in haloarchaea experimentally the monocistronic sod gene was chosen, which contains a 5′-UTR with an extensive SD motif of seven nucleotides and a length of 19 nt, the average length of 5′UTRs in this organism. A translational fusion of part of the sod gene with the dhfr reporter gene was constructed. A mutant series was generated that matched the SD motif from zero to eight positions, respectively. Surprisingly, there was no correlation between the base pairing ability between transcripts and 16S rRNA and translational efficiency in vivo under several different growth conditions. Furthermore, complete replacement of the SD motif by three unrelated sequences did not reduce translational efficiency. The results indicate that H. volcanii does not make use of the SD mechanism for translation initiation in 5′-UTRs. A genome analysis revealed that while the number of SD motifs in 5′-UTRs is rare, their fraction within open reading frames is high. Possible biological functions for intragenic SD motifs are discussed, including re-initiation of translation at distal genes in operons.

Alterations in rRNA-mRNA interaction during plastid evolution

Translation initiation depends on the recognition of mRNA by a ribosome. For this to occur, prokaryotes primarily use the Shine-Dalgarno (SD) interaction, where the 3'-tail of small subunit rRNA (core motif: 3'CCUCC) forms base pairs with a complementary signal sequence in the 5'-untranslated region of mRNA. Here, we examined what happened to SD interactions during the evolution of a cyanobacterial endosymbiont into modern plastids (including chloroplasts). Our analysis of available complete plastid genome sequences revealed that the majority of plastids retained SD interactions but with varying levels of usage. Parallel losses of SD interactions took place in plastids of Chlorophyta, Euglenophyta, and Chromerida/Apicomplexa lineages, presumably related to their extensive reductive evolution. Interestingly, we discovered that the classical SD interaction (3'CCUCC/5'GGAGG [rRNA/mRNA]) was replaced by an altered SD interaction (3'CCCU/5'GGGA or 3'CUUCC/5'GAAGG) through coordinated changes in the sequences of the core rRNA motif and its paired mRNA signal. These changes in plastids of Chlorophyta and Euglenophyta proceeded through intermediate stages that allowed both the classical and altered SD interactions. This coevolution between the rRNA motif and the mRNA signal demonstrates unexpected plasticity in the translation initiation machinery.

Insertion sequence distribution bias in Archaea

Insertion sequences (IS) are common transposable elements in Archaea. Intergenic IS elements are usually less harmful than intragenic ISs, simply because they are less likely to disrupt host gene function. However, because regulatory sequences are intergenic and upstream of genes, we hypothesized that not all intergenic regions are selectively equivalent for IS insertion. We tested this hypothesis by analyzing the distributions of intergenic IS elements within 155 fully sequenced archaeal genomes. Of the 22 genomes with enough IS elements for statistical analysis, five have significantly fewer ISs between divergently oriented neighboring genes than expected by chance, and seven have significantly more ISs between convergently oriented genes. Furthermore, of the 85 genomes with at least one expected IS within each of the three possible neighboring gene orientations (i.e., divergent, convergent, and tandem), 73 genomes have fewer ISs between divergently oriented genes than expected, and 60 have more ISs between convergently oriented genes than expected (both values deviate significantly from binomial probabilities of random distribution). We suspect that these non-random IS distributions are molded by natural selection resulting from differential disruption of neighboring gene regulation, and that this selective pressure has affected transposable element distributions in prokaryotes for billions of years.

A sporulation-specific, sigF-dependent protein, SspA, affects septum positioning in Streptomyces coelicolor

The RNA polymerase sigma factor SigF controls late development during sporulation in the filamentous bacterium Streptomyces coelicolor. The only known SigF-dependent gene identified so far, SCO5321, is found in the biosynthetic cluster encoding spore pigment synthesis. Here we identify the first direct target for SigF, the gene sspA, encoding a sporulation-specific protein. Bioinformatic analysis suggests that SspA is a secreted lipoprotein with two PepSY signature domains. The sspA deletion mutant exhibits irregular sporulation septation and altered spore shape, suggesting that SspA plays a role in septum formation and spore maturation. The fluorescent translational fusion protein SspA–mCherry localized first to septum sites, then subsequently around the surface of the spores. Both SspA protein and sspA transcription are absent from the sigF null mutant. Moreover, in vitro transcription assay confirmed that RNA polymerase holoenzyme containing SigF is sufficient for initiation of transcription from a single sspA promoter. In addition, in vivo and in vitro experiments showed that sspA is a direct target of BldD, which functions to repress sporulation genes, including whiG, ftsZ and ssgB, during vegetative growth, co-ordinating their expression during sporulation septation.

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 comprehensive analysis of the importance of translation initiation factors for Haloferax volcanii applying deletion and conditional depletion mutants

Translation is an important step in gene expression. The initiation of translation is phylogenetically diverse, since currently five different initiation mechanisms are known. For bacteria the three initiation factors IF1 – IF3 are described in contrast to archaea and eukaryotes, which contain a considerably higher number of initiation factor genes. As eukaryotes and archaea use a non-overlapping set of initiation mechanisms, orthologous proteins of both domains do not necessarily fulfill the same function. The genome of Haloferax volcanii contains 14 annotated genes that encode (subunits of) initiation factors. To gain a comprehensive overview of the importance of these genes, it was attempted to construct single gene deletion mutants of all genes. In 9 cases single deletion mutants were successfully constructed, showing that the respective genes are not essential. In contrast, the genes encoding initiation factors aIF1, aIF2γ, aIF5A, aIF5B, and aIF6 were found to be essential. Factors aIF1A and aIF2β are encoded by two orthologous genes in H. volcanii. Attempts to generate double mutants failed in both cases, indicating that also these factors are essential. A translatome analysis of one of the single aIF2β deletion mutants revealed that the translational efficiency of the second ortholog was enhanced tenfold and thus the two proteins can replace one another. The phenotypes of the single deletion mutants also revealed that the two aIF1As and aIF2βs have redundant but not identical functions. Remarkably, the gene encoding aIF2α, a subunit of aIF2 involved in initiator tRNA binding, could be deleted. However, the mutant had a severe growth defect under all tested conditions. Conditional depletion mutants were generated for the five essential genes. The phenotypes of deletion mutants and conditional depletion mutants were compared to that of the wild-type under various conditions, and growth characteristics are discussed.

RNA-Seq of Bacillus licheniformis: active regulatory RNA features expressed within a productive fermentation

BACKGROUND

The production of enzymes by an industrial strain requires a complex adaption of the bacterial metabolism to the conditions within the fermenter. Regulatory events within the process result in a dynamic change of the transcriptional activity of the genome. This complex network of genes is orchestrated by proteins as well as regulatory RNA elements. Here we present an RNA-Seq based study considering selected phases of an industry-oriented fermentation of Bacillus licheniformis.

RESULTS

A detailed analysis of 20 strand-specific RNA-Seq datasets revealed a multitude of transcriptionally active genomic regions. 3314 RNA features encoded by such active loci have been identified and sorted into ten functional classes. The identified sequences include the expected RNA features like housekeeping sRNAs, metabolic riboswitches and RNA switches well known from studies on Bacillus subtilis as well as a multitude of completely new candidates for regulatory RNAs. An unexpectedly high number of 855 RNA features are encoded antisense to annotated protein and RNA genes, in addition to 461 independently transcribed small RNAs. These antisense transcripts contain molecules with a remarkable size range variation from 38 to 6348 base pairs in length. The genome of the type strain B. licheniformis DSM13 was completely reannotated using data obtained from RNA-Seq analyses and from public databases.

CONCLUSION

The hereby generated data-sets represent a solid amount of knowledge on the dynamic transcriptional activities during the investigated fermentation stages. The identified regulatory elements enable research on the understanding and the optimization of crucial metabolic activities during a productive fermentation of Bacillus licheniformis strains.

RNA-Seq analysis of Mycobacterium avium non-coding transcriptome

Deep sequencing was implemented to study the transcriptional landscape of Mycobacterium avium. High-resolution transcriptome analysis identified the transcription start points for 652 genes. One third of these genes represented leaderless transcripts, whereas the rest of the transcripts had 5' UTRs with the mean length of 83 nt. In addition, the 5' UTRs of 6 genes contained SAM-IV and Ykok types of riboswitches. 87 antisense RNAs and 10 intergenic small RNAs were mapped. 6 intergenic small RNAs, including 4.5S RNA and rnpB, were transcribed at extremely high levels. Although several intergenic sRNAs are conserved in M. avium and M. tuberculosis, both of these species have unique intergenic sRNAs. Moreover, we demonstrated that even conserved small RNAs are regulated differently in these species. Different sets of intergenic sRNAs may underlie differences in physiology between conditionally pathogenic M. avium and highly specialized pathogen M. tuberculosis.

Active and accurate trans-translation requires distinct determinants in the C-terminal tail of SmpB protein and the mRNA-like domain of transfer messenger RNA (tmRNA)

Unproductive ribosome stalling in eubacteria is resolved by the actions of SmpB protein and transfer messenger (tm) RNA. We examined the functional significance of conserved regions of SmpB and tmRNA to the trans-translation process. Our investigations reveal that the N-terminal 20 residues of SmpB, which are located near the ribosomal decoding center, are dispensable for all known SmpB activities. In contrast, a set of conserved residues that reside at the junction between the tmRNA-binding core and the C-terminal tail of SmpB play an important role in tmRNA accommodation. Our data suggest that the highly conserved glycine 132 acts as a flexible hinge that enables movement of the C-terminal tail, thus permitting proper positioning and establishment of the tmRNA open reading frame (ORF) as the surrogate template. To gain further insights into the function of the SmpB C-terminal tail, we examined the tagging activity of hybrid variants of tmRNA and the SmpB protein, in which the tmRNA ORF or the SmpB C-terminal tail was substituted with the equivalent but highly divergent sequences from Francisella tularensis. We observed that the hybrid tmRNA was active but resulted in less accurate selection of the resume codon. Cognate hybrid SmpB was necessary to restore activity. Furthermore, accurate tagging was observed when the identity of the resume codon was reverted from GGC to GCA. Taken together, these data suggest that the engagement of the tmRNA ORF and the selection of the correct translation resumption point are distinct activities that are influenced by independent tmRNA and SmpB determinants.

Isolation of fully synthetic promoters for high-level gene expression in Corynebacterium glutamicum

Corynebacterium glutamicum is an important industrial organism that is widely used in the production of amino acids, nucleotides and vitamins. To extend its product spectrum and improve productivity, C. glutamicum needs to undergo further engineering, including the development of applicable promoter system. Here, we isolated new promoters from the fully synthetic promoter library consisting of 70-bp random sequences in C. glutamicum. Using green fluorescent protein (GFP) as a reporter, highly fluorescent cells were screened from the library by fluorescent activated cell sorting (FACS). Twenty potential promoters of various strengths were isolated and characterized through extensive analysis of DNA sequences and mRNA transcripts. Among 20 promoters, 6 promoters which have different strengths were selected and their activities were successfully demonstrated using two model proteins (antibody fragment and endoxylanase). Finally, the strongest promoter (P(H36)) was employed for the secretory production of endoxylanase in fed-batch cultivation, achieving production levels of 746 mg/L in culture supernatant. This is the first report of synthetic promoters constructed in C. glutamicum, and our screening strategy together with the use of synthetic promoters of various strengths will contribute to the future engineering of C. glutamicum.

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.

Histidine biosynthesis, its regulation and biotechnological application in Corynebacterium glutamicum

l-Histidine biosynthesis is an ancient metabolic pathway present in bacteria, archaea, lower eukaryotes, and plants. For decades l-histidine biosynthesis has been studied mainly in Escherichia coli and Salmonella typhimurium, revealing fundamental regulatory processes in bacteria. Furthermore, in the last 15 years this pathway has been also investigated intensively in the industrial amino acid-producing bacterium Corynebacterium glutamicum, revealing similarities to E. coli and S. typhimurium, as well as differences. This review summarizes the current knowledge of l-histidine biosynthesis in C. glutamicum. The genes involved and corresponding enzymes are described, in particular focusing on the imidazoleglycerol-phosphate synthase (HisFH) and the histidinol-phosphate phosphatase (HisN). The transcriptional organization of his genes in C. glutamicum is also reported, including the four histidine operons and their promoters. Knowledge of transcriptional regulation during stringent response and by histidine itself is summarized and a translational regulation mechanism is discussed, as well as clues about a histidine transport system. Finally, we discuss the potential of using this knowledge to create or improve C. glutamicum strains for the industrial l-histidine production.

Gene prediction in metagenomic fragments based on the SVM algorithm

Background

Metagenomic sequencing is becoming a powerful technology for exploring micro-ogranisms from various environments, such as human body, without isolation and cultivation. Accurately identifying genes from metagenomic fragments is one of the most fundamental issues.

Results

In this article, we present a novel gene prediction method named MetaGUN for metagenomic fragments based on a machine learning approach of SVM. It implements in a three-stage strategy to predict genes. Firstly, it classifies input fragments into phylogenetic groups by a k-mer based sequence binning method. Then, protein-coding sequences are identified for each group independently with SVM classifiers that integrate entropy density profiles (EDP) of codon usage, translation initiation site (TIS) scores and open reading frame (ORF) length as input patterns. Finally, the TISs are adjusted by employing a modified version of MetaTISA. To identify protein-coding sequences, MetaGun builds the universal module and the novel module. The former is based on a set of representative species, while the latter is designed to find potential functionary DNA sequences with conserved domains.

Conclusions

Comparisons on artificial shotgun fragments with multiple current metagenomic gene finders show that MetaGUN predicts better results on both 3' and 5' ends of genes with fragments of various lengths. Especially, it makes the most reliable predictions among these methods. As an application, MetaGUN was used to predict genes for two samples of human gut microbiome. It identifies thousands of additional genes with significant evidences. Further analysis indicates that MetaGUN tends to predict more potential novel genes than other current metagenomic gene finders.

Regulation of pluripotency and self- renewal of ESCs through epigenetic-threshold modulation and mRNA pruning

Embryonic stem cell (ESC) pluripotency requires bivalent epigenetic modifications of key developmental genes regulated by various transcription factors and chromatin-modifying enzymes. How these factors coordinate with one another to maintain the bivalent chromatin state so that ESCs can undergo rapid self-renewal while retaining pluripotency is poorly understood. We report that Utf1, a target of Oct4 and Sox2, is a bivalent chromatin component that buffers poised states of bivalent genes. By limiting PRC2 loading and histone 3 lysine-27 trimethylation, Utf1 sets proper activation thresholds for bivalent genes. It also promotes nuclear tagging of messenger RNAs (mRNAs) transcribed from insufficiently silenced bivalent genes for cytoplasmic degradation through mRNA decapping. These opposing functions of Utf1 promote coordinated differentiation. The mRNA degradation function also ensures rapid cell proliferation by blocking the Myc-Arf feedback control. Thus, Utf1 couples the core pluripotency factors with Myc and PRC2 networks to promote the pluripotency and proliferation of ESCs.

Kinetic control of translation initiation in bacteria

Translation initiation is a crucial step of protein synthesis which largely defines how the composition of the cellular transcriptome is converted to the proteome and controls the response and adaptation to environmental stimuli. The efficiency of translation of individual mRNAs, and hence the basal shape of the proteome, is defined by the structures of the mRNA translation initiation regions. Initiation efficiency can be regulated by small molecules, proteins, or antisense RNAs, underscoring its importance in translational control. Although initiation has been studied in bacteria for decades, many aspects remain poorly understood. Recent evidence has suggested an unexpected diversity of pathways by which mRNAs can be recruited to the bacterial ribosome, the importance of structural dynamics of initiation intermediates, and the complexity of checkpoints for mRNA selection. In this review, we discuss how the ribosome shapes the landscape of translation initiation by non-linear kinetic processing of the transcriptome information. We summarize the major pathways by which mRNAs enter the ribosome depending on the structure of their 5' untranslated regions, the assembly and the structure of initiation intermediates, the individual and synergistic roles of initiation factors, and the mechanisms of mRNA and initiator tRNA selection.

The Pseudomonas aeruginosa transcriptome in planktonic cultures and static biofilms using RNA sequencing

In this study, we evaluated how gene expression differs in mature Pseudomonas aeruginosa biofilms as opposed to planktonic cells by the use of RNA sequencing technology that gives rise to both quantitative and qualitative information on the transcriptome. Although a large proportion of genes were consistently regulated in both the stationary phase and biofilm cultures as opposed to the late exponential growth phase cultures, the global biofilm gene expression pattern was clearly distinct indicating that biofilms are not just surface attached cells in stationary phase. A large amount of the genes found to be biofilm specific were involved in adaptation to microaerophilic growth conditions, repression of type three secretion and production of extracellular matrix components. Additionally, we found many small RNAs to be differentially regulated most of them similarly in stationary phase cultures and biofilms. A qualitative analysis of the RNA-seq data revealed more than 3000 putative transcriptional start sites (TSS). By the use of rapid amplification of cDNA ends (5′-RACE) we confirmed the presence of three different TSS associated with the pqsABCDE operon, two in the promoter of pqsA and one upstream of the second gene, pqsB. Taken together, this study reports the first transcriptome study on P. aeruginosa that employs RNA sequencing technology and provides insights into the quantitative and qualitative transcriptome including the expression of small RNAs in P. aeruginosa biofilms.