Comparative genomic analysis of translation initiation mechanisms for genes lacking the Shine-Dalgarno sequence in prokaryotes

Exploring the roles of ribosomal peptides in prokaryote-phage interactions through deep learning-enabled metagenome mining

BACKGROUND

Microbial secondary metabolites play a crucial role in the intricate interactions within the natural environment. Among these metabolites, ribosomally synthesized and post-translationally modified peptides (RiPPs) are becoming a promising source of therapeutic agents due to their structural diversity and functional versatility. However, their biosynthetic capacity and ecological functions remain largely underexplored.

RESULTS

Here, we aim to explore the biosynthetic profile of RiPPs and their potential roles in the interactions between microbes and viruses in the ocean, which encompasses a vast diversity of unique biomes that are rich in interactions and remains chemically underexplored. We first developed TrRiPP to identify RiPPs from ocean metagenomes, a deep learning method that detects RiPP precursors in a hallmark gene-independent manner to overcome the limitations of classic methods in processing highly fragmented metagenomic data. Applying this method to metagenomes from the global ocean microbiome, we uncover a diverse array of previously uncharacterized putative RiPP families with great novelty and diversity. Through correlation analysis based on metatranscriptomic data, we observed a high prevalence of antiphage defense-related and phage-related protein families that were co-expressed with RiPP families. Based on this putative association between RiPPs and phage infection, we constructed an Ocean Virus Database (OVD) and established a RiPP-involving host-phage interaction network through host prediction and co-expression analysis, revealing complex connectivities linking RiPP-encoding prokaryotes, RiPP families, viral protein families, and phages. These findings highlight the potential of RiPP families involved in prokaryote-phage interactions and coevolution, providing insights into their ecological functions in the ocean microbiome.

CONCLUSIONS

This study provides a systematic investigation of the biosynthetic potential of RiPPs from the ocean microbiome at a global scale, shedding light on the essential insights into the ecological functions of RiPPs in prokaryote-phage interactions through the integration of deep learning approaches, metatranscriptomic data, and host-phage connectivity. This study serves as a valuable example of exploring the ecological functions of bacterial secondary metabolites, particularly their associations with unexplored microbial interactions. Video Abstract.

Analysis of programmed frameshifting during translation of prfB in Flavobacterium johnsoniae

Ribosomes of Bacteroidia fail to recognize Shine-Dalgarno (SD) sequences due to sequestration of the 3' tail of the 16S rRNA on the 30S platform. Yet in these organisms, the prfB gene typically contains the programmed +1 frameshift site with its characteristic SD sequence. Here, we investigate prfB autoregulation in Flavobacterium johnsoniae, a member of the Bacteroidia. We find that the efficiency of prfB frameshifting in F. johnsoniae is low (∼7%) relative to that in Escherichia coli (∼50%). Mutation or truncation of bS21 in F. johnsoniae increases frameshifting substantially, suggesting that anti-SD (ASD) sequestration is responsible for the reduced efficiency. The frameshift site of certain Flavobacteriales, such as Winogradskyella psychrotolerans, has no SD. In F. johnsoniae, this W. psychrotolerans sequence supports frameshifting as well as the native sequence, and mutation of bS21 causes no enhancement. These data suggest that prfB frameshifting normally occurs without SD-ASD pairing, at least under optimal laboratory growth conditions. Chromosomal mutations that remove the frameshift or ablate the SD confer subtle growth defects in the presence of paraquat or streptomycin, respectively, indicating that both the autoregulatory mechanism and the SD element contribute to F. johnsoniae cell fitness. Analysis of prfB frameshift sites across 2686 representative bacteria shows loss of the SD sequence in many clades, with no obvious relationship to genome-wide SD usage. These data reveal unexpected variation in the mechanism of frameshifting and identify another group of organisms, the Verrucomicrobiales, that globally lack SD sequences.

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.

The distinct translational landscapes of gram-negative Salmonella and gram-positive Listeria

Translational control in pathogenic bacteria is fundamental to gene expression and affects virulence and other infection phenotypes. We used an enhanced ribosome profiling protocol coupled with parallel transcriptomics to capture accurately the global translatome of two evolutionarily distant pathogenic bacteria-the Gram-negative bacterium Salmonella and the Gram-positive bacterium Listeria. We find that the two bacteria use different mechanisms to translationally regulate protein synthesis. In Salmonella, in addition to the expected correlation between translational efficiency and cis-regulatory features such as Shine-Dalgarno (SD) strength and RNA secondary structure around the initiation codon, our data reveal an effect of the 2nd and 3rd codons, where the presence of tandem lysine codons (AAA-AAA) enhances translation in both Salmonella and E. coli. Strikingly, none of these features are seen in efficiently translated Listeria transcripts. Instead, approximately 20% of efficiently translated Listeria genes exhibit 70 S footprints seven nt upstream of the authentic start codon, suggesting that these genes may be subject to a novel translational initiation mechanism. Our results show that SD strength is not a direct hallmark of translational efficiency in all bacteria. Instead, Listeria has evolved additional mechanisms to control gene expression level that are distinct from those utilised by Salmonella and E. coli.

Chloroplast gene expression: Recent advances and perspectives

Chloroplasts evolved from an ancient cyanobacterial endosymbiont more than 1.5 billion years ago. During subsequent coevolution with the nuclear genome, the chloroplast genome has remained independent, albeit strongly reduced, with its own transcriptional machinery and distinct features, such as chloroplast-specific innovations in gene expression and complicated post-transcriptional processing. Light activates the expression of chloroplast genes via mechanisms that optimize photosynthesis, minimize photodamage, and prioritize energy investments. Over the past few years, studies have moved from describing phases of chloroplast gene expression to exploring the underlying mechanisms. In this review, we focus on recent advances and emerging principles that govern chloroplast gene expression in land plants. We discuss engineering of pentatricopeptide repeat proteins and its biotechnological effects on chloroplast RNA research; new techniques for characterizing the molecular mechanisms of chloroplast gene expression; and important aspects of chloroplast gene expression for improving crop yield and stress tolerance. We also discuss biological and mechanistic questions that remain to be answered in the future.

Identification and validation of ferroptosis-related genes in patients infected with dengue virus: implication in the pathogenesis of DENV

Ferroptosis, an iron-dependent form of regulated cell death, has been associated with many virus infections. However, the role of ferroptosis in dengue virus (DENV) infection remains to be clarified. In our study, a dengue fever microarray dataset (GSE51808) of whole blood samples was downloaded from the Gene Expression Omnibus (GEO), and a list of ferroptosis related genes (FRGs) was extracted from the FerrDb. We identified 37 ferroptosis-related differentially expressed genes (FR-DEGs) in DENV-infected patient blood samples compared to healthy individuals. Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses as well as protein-protein interaction (PPI) network of FR-DEGs revealed that these 37 FR-DEGs were mainly related to the C-type lectin receptor and p53 signaling pathway. Nine out of the 37 FR-DEGs (HSPA5, CAV1, HRAS, PTGS2, JUN, IL6, ATF3, XBP1, and CDKN2A) were hub genes, of which 5 were validated by qRT-PCR in DENV-infected HepG2 cells. Finally, using miRNA-mRNA regulatory network, we identified has-miR-124-3p and has-miR-16-5p as the most critical miRNAs in regulating the expression of these hub genes. In conclusion, our findings demonstrated that 5 FR-DEGs, JUN, IL6, ATF3, XBP1, and CDKN2A, and two miRNAs, has-miR-124-3p and has-miR-16-5p may implicate an essential role of ferroptosis in DENV infection, and further studies are warranted to explore the underlying mechanisms.

Ultradeep characterisation of translational sequence determinants refutes rare-codon hypothesis and unveils quadruplet base pairing of initiator tRNA and transcript

Translation is a key determinant of gene expression and an important biotechnological engineering target. In bacteria, 5'-untranslated region (5'-UTR) and coding sequence (CDS) are well-known mRNA parts controlling translation and thus cellular protein levels. However, the complex interaction of 5'-UTR and CDS has so far only been studied for few sequences leading to non-generalisable and partly contradictory conclusions. Herein, we systematically assess the dynamic translation from over 1.2 million 5'-UTR-CDS pairs in Escherichia coli to investigate their collective effect using a new method for ultradeep sequence-function mapping. This allows us to disentangle and precisely quantify effects of various sequence determinants of translation. We find that 5'-UTR and CDS individually account for 53% and 20% of variance in translation, respectively, and show conclusively that, contrary to a common hypothesis, tRNA abundance does not explain expression changes between CDSs with different synonymous codons. Moreover, the obtained large-scale data provide clear experimental evidence for a base-pairing interaction between initiator tRNA and mRNA beyond the anticodon-codon interaction, an effect that is often masked for individual sequences and therefore inaccessible to low-throughput approaches. Our study highlights the indispensability of ultradeep sequence-function mapping to accurately determine the contribution of parts and phenomena involved in gene regulation.

Ribosomes lacking bS21 gain function to regulate protein synthesis in Flavobacterium johnsoniae

Ribosomes of Bacteroidia (formerly Bacteroidetes) fail to recognize Shine-Dalgarno (SD) sequences even though they harbor the anti-SD (ASD) of 16S rRNA. Inhibition of SD-ASD pairing is due to sequestration of the 3' tail of 16S rRNA in a pocket formed by bS21, bS18, and bS6 on the 30S platform. Interestingly, in many Flavobacteriales, the gene encoding bS21, rpsU, contains an extended SD sequence. In this work, we present genetic and biochemical evidence that bS21 synthesis in Flavobacterium johnsoniae is autoregulated via a subpopulation of ribosomes that specifically lack bS21. Mutation or depletion of bS21 in the cell increases translation of reporters with strong SD sequences, such as rpsU'-gfp, but has no effect on other reporters. Purified ribosomes lacking bS21 (or its C-terminal region) exhibit higher rates of initiation on rpsU mRNA and lower rates of initiation on other (SD-less) mRNAs than control ribosomes. The mechanism of autoregulation depends on extensive pairing between mRNA and 16S rRNA, and exceptionally strong SD sequences, with predicted pairing free energies of < -13 kcal/mol, are characteristic of rpsU across the Bacteroidota. This work uncovers a clear example of specialized ribosomes in bacteria.

Initiation at AUGUG and GUGUG sequences can lead to translation of overlapping reading frames in E. coli

During initiation, the ribosome is tasked to efficiently recognize open reading frames (ORFs) for accurate and fast translation of mRNAs. A critical step is start codon recognition, which is modulated by initiation factors, mRNA structure, a Shine Dalgarno (SD) sequence and the start codon itself. Within the Escherichia coli genome, we identified more than 50 annotated initiation sites harboring AUGUG or GUGUG sequence motifs that provide two canonical start codons, AUG and GUG, in immediate proximity. As these sites may challenge start codon recognition, we studied if and how the ribosome is accurately guided to the designated ORF, with a special focus on the SD sequence as well as adenine at the fourth coding sequence position (A4). By in vitro and in vivo experiments, we characterized key requirements for unambiguous start codon recognition, but also discovered initiation sites that lead to the translation of both overlapping reading frames. Our findings corroborate the existence of an ambiguous translation initiation mechanism, implicating a multitude of so far unrecognized ORFs and translation products in bacteria.

INRI-seq enables global cell-free analysis of translation initiation and off-target effects of antisense inhibitors

Ribosome profiling (Ribo-seq) is a powerful method for the transcriptome-wide assessment of protein synthesis rates and the study of translational control mechanisms. Yet, Ribo-seq also has limitations. These include difficulties with the analysis of translation-modulating molecules such as antibiotics, which are often toxic or challenging to deliver into living cells. Here, we have developed in vitro Ribo-seq (INRI-seq), a cell-free method to analyze the translational landscape of a fully customizable synthetic transcriptome. Using Escherichia coli as an example, we show how INRI-seq can be used to analyze the translation initiation sites of a transcriptome of interest. We also study the global impact of direct translation inhibition by antisense peptide nucleic acid (PNA) to analyze PNA off-target effects. Overall, INRI-seq presents a scalable, sensitive method to study translation initiation in a transcriptome-wide manner without the potentially confounding effects of extracting ribosomes from living cells.

Airway microecology in rifampicin-resistant and rifampicin-sensitive pulmonary tuberculosis patients

BACKGROUND

Pulmonary tuberculosis is a chronic infectious disease of the respiratory system. It is still one of the leading causes of death from a single infectious disease, but it has been stuck in the study of a single pathogen. Recent studies have shown that many diseases are associated with disruption of the native microbiota. In this study we investigated the occurrence of tuberculosis and the correlation between drug resistance and respiratory flora. High-throughput 16 S rRNA gene sequencing was used to characterize the respiratory microbiota composition of 30 tuberculosis (TB) affected patients and compared with 30 healthy (H) controls. According to their Gene Xpert results, 30 pulmonary tuberculosis patients were divided into 12 persons in the drug-sensitive group (DS0) and 18 persons in the drug-resistant group (DR0). The microbial flora of the two were compared with the H group.

RESULTS

The data generated by sequencing showed that Firmicutes, Proteus, Bacteroides, Actinomyces and Fusobacterium were the five main bacterial phyla detected, and they constituted more than 96% of the microbial community. The relative abundances of Fusobacterium, Haemophilus, Porphyromonas, Neisseria, TM7, Spirochetes, SR1, and Tenericutes in the TB group was lower than that of the H group, and Granulicatella was higher than the H group. The PcoA diagrams of the two groups had obvious clustering differences. The Alpha diversity of the TB group was lower than that of the H group, and the Beta diversity was higher than that of the H group (P < 0.05). The relative abundance of Streptococcus in the DS0 group was significantly higher than that in the DR0 group (P < 0.05).

CONCLUSION

Pulmonary tuberculosis can cause disorders of the respiratory tract microbial flora, in which the relative abundance of Streptococcus was significantly different between rifampicin-sensitive and rifampicin-resistant patients.

Mechanisms governing codon usage bias and the implications for protein expression in the chloroplast of Chlamydomonas reinhardtii

Chloroplasts possess a considerably reduced genome that is decoded via an almost minimal set of tRNAs. These features make an excellent platform for gaining insights into fundamental mechanisms that govern protein expression. Here, we present a comprehensive and revised perspective of the mechanisms that drive codon selection in the chloroplast of Chlamydomonas reinhardtii and the functional consequences for protein expression. In order to extract this information, we applied several codon usage descriptors to genes with different expression levels. We show that highly expressed genes strongly favor translationally optimal codons, while genes with lower functional importance are rather affected by directional mutational bias. We demonstrate that codon optimality can be deduced from codon-anticodon pairing affinity and, for a small number of amino acids (leucine, arginine, serine, and isoleucine), tRNA concentrations. Finally, we review, analyze, and expand on the impact of codon usage on protein yield, secondary structures of mRNA, translation initiation and termination, and amino acid composition of proteins, as well as cotranslational protein folding. The comprehensive analysis of codon choice provides crucial insights into heterologous gene expression in the chloroplast of C. reinhardtii, which may also be applicable to other chloroplast-containing organisms and bacteria.

Role of aIF5B in archaeal translation initiation

In eukaryotes and in archaea late steps of translation initiation involve the two initiation factors e/aIF5B and e/aIF1A. In eukaryotes, the role of eIF5B in ribosomal subunit joining is established and structural data showing eIF5B bound to the full ribosome were obtained. To achieve its function, eIF5B collaborates with eIF1A. However, structural data illustrating how these two factors interact on the small ribosomal subunit have long been awaited. The role of the archaeal counterparts, aIF5B and aIF1A, remains to be extensively addressed. Here, we study the late steps of Pyrococcus abyssi translation initiation. Using in vitro reconstituted initiation complexes and light scattering, we show that aIF5B bound to GTP accelerates subunit joining without the need for GTP hydrolysis. We report the crystallographic structures of aIF5B bound to GDP and GTP and analyze domain movements associated to these two nucleotide states. Finally, we present the cryo-EM structure of an initiation complex containing 30S bound to mRNA, Met-tRNAiMet, aIF5B and aIF1A at 2.7 Å resolution. Structural data shows how archaeal 5B and 1A factors cooperate to induce a conformation of the initiator tRNA favorable to subunit joining. Archaeal and eukaryotic features of late steps of translation initiation are discussed.

Identification and Cross-Characterisation of Artificial Promoters and 5' Untranslated Regions in Vibrio natriegens

Vibrio natriegens has recently gained attention as a novel fast-growing bacterium in synthetic biology applications. Currently, a limited set of genetic elements optimised for Escherichia coli are used in V. natriegens due to the lack of DNA parts characterised in this novel host. In this study, we report the identification and cross-characterisation of artificial promoters and 5' untranslated regions (artificial regulatory sequence, ARES) that lead to production of fluorescent proteins with a wide-range of expression levels. We identify and cross-characterise 52 constructs in V. natriegens and E. coli. Furthermore, we report the DNA sequence and motif analysis of the ARESs using various algorithms. With this study, we expand the pool of characterised genetic DNA parts that can be used for different biotechnological applications using V. natriegens as a host microorganism.

Comparative Analysis of anti-Shine- Dalgarno Function in Flavobacterium johnsoniae and Escherichia coli

The anti-Shine-Dalgarno (ASD) sequence of 16S rRNA is highly conserved across Bacteria, and yet usage of Shine-Dalgarno (SD) sequences in mRNA varies dramatically, depending on the lineage. Here, we compared the effects of ASD mutagenesis in Escherichia coli, a Gammaproteobacteria which commonly employs SD sequences, and Flavobacterium johnsoniae, a Bacteroidia which rarely does. In E. coli, 30S subunits carrying any single substitution at positions 1,535–1,539 confer dominant negative phenotypes, whereas subunits with mutations at positions 1,540–1,542 are sufficient to support cell growth. These data suggest that CCUCC (1,535–1,539) represents the functional core of the element in E. coli. In F. johnsoniae, deletion of three ribosomal RNA (rrn) operons slowed growth substantially, a phenotype largely rescued by a plasmid-borne copy of the rrn operon. Using this complementation system, we found that subunits with single mutations at positions 1,535–1,537 are as active as control subunits, in sharp contrast to the E. coli results. Moreover, subunits with quadruple substitution or complete replacement of the ASD retain substantial, albeit reduced, activity. Sedimentation analysis revealed that these mutant subunits are overrepresented in the subunit fractions and underrepresented in polysome fractions, suggesting some defect in 30S biogenesis and/or translation initiation. Nonetheless, our collective data indicate that the ASD plays a much smaller role in F. johnsoniae than in E. coli, consistent with SD usage in the two organisms.

The Role of the Universally Conserved ATPase YchF/Ola1 in Translation Regulation during Cellular Stress

The ability to respond to metabolic or environmental changes is an essential feature in all cells and involves both transcriptional and translational regulators that adjust the metabolic activity to fluctuating conditions. While transcriptional regulation has been studied in detail, the important role of the ribosome as an additional player in regulating gene expression is only beginning to emerge. Ribosome-interacting proteins are central to this translational regulation and include universally conserved ribosome interacting proteins, such as the ATPase YchF (Ola1 in eukaryotes). In both eukaryotes and bacteria, the cellular concentrations of YchF/Ola1 determine the ability to cope with different stress conditions and are linked to several pathologies in humans. The available data indicate that YchF/Ola1 regulates the stress response via controlling non-canonical translation initiation and via protein degradation. Although the molecular mechanisms appear to be different between bacteria and eukaryotes, increased non-canonical translation initiation is a common consequence of YchF/Ola1 regulated translational control in E. coli and H. sapiens. In this review, we summarize recent insights into the role of the universally conserved ATPase YchF/Ola1 in adapting translation to unfavourable conditions.

StartLink and StartLink+: Prediction of Gene Starts in Prokaryotic Genomes

State-of-the-art algorithms of ab initio gene prediction for prokaryotic genomes were shown to be sufficiently accurate. A pair of algorithms would agree on predictions of gene 3'ends. Nonetheless, predictions of gene starts would not match for 15-25% of genes in a genome. This discrepancy is a serious issue that is difficult to be resolved due to the absence of sufficiently large sets of genes with experimentally verified starts. We have introduced StartLink that infers gene starts from conservation patterns revealed by multiple alignments of homologous nucleotide sequences. We also have introduced StartLink+ combining both ab initio and alignment-based methods. The ability of StartLink to predict the start of a given gene is restricted by the availability of homologs in a database. We observed that StartLink made predictions for 85% of genes per genome on average. The StartLink+ accuracy was shown to be 98-99% on the sets of genes with experimentally verified starts. In comparison with database annotations, we observed that the annotated gene starts deviated from the StartLink+ predictions for ∼5% of genes in AT-rich genomes and for 10-15% of genes in GC-rich genomes on average. The use of StartLink+ has a potential to significantly improve gene start annotation in genomic databases.

Importance of the 5' regulatory region to bacterial synthetic biology applications

The field of synthetic biology is evolving at a fast pace. It is advancing beyond single-gene alterations in single hosts to the logical design of complex circuits and the development of integrated synthetic genomes. Recent breakthroughs in deep learning, which is increasingly used in de novo assembly of DNA components with predictable effects, are also aiding the discipline. Despite advances in computing, the field is still reliant on the availability of pre-characterized DNA parts, whether natural or synthetic, to regulate gene expression in bacteria and make valuable compounds. In this review, we discuss the different bacterial synthetic biology methodologies employed in the creation of 5' regulatory regions - promoters, untranslated regions and 5'-end of coding sequences. We summarize methodologies and discuss their significance for each of the functional DNA components, and highlight the key advances made in bacterial engineering by concentrating on their flaws and strengths. We end the review by outlining the issues that the discipline may face in the near future.

Translation of a Leaderless Reporter Is Robust During Exponential Growth and Well Sustained During Stress Conditions in Mycobacterium tuberculosis

Mycobacterium tuberculosis expresses a large number of leaderless mRNA transcripts; these lack the 5' leader region, which usually contains the Shine-Dalgarno sequence required for translation initiation in bacteria. In M. tuberculosis, transcripts encoding proteins like toxin-antitoxin systems are predominantly leaderless and the overall ratio of leaderless to Shine-Dalgarno transcripts significantly increases during growth arrest, suggesting that leaderless translation might be important during persistence in the host. However, whether these two types of transcripts are translated with differing efficiencies during optimal growth conditions and during stress conditions that induce growth arrest, is unclear. Here, we have used the desA1 (Rv0824c) and desA2 (Rv1094) gene pair as representative for Shine-Dalgarno and leaderless transcripts in M. tuberculosis respectively; and used them to construct bioluminescent reporter strains. We detect robust leaderless translation during exponential in vitro growth, and we show that leaderless translation is more stable than Shine-Dalgarno translation during adaptation to stress conditions. These changes are independent from transcription, as transcription levels did not significantly change following quantitative real-time PCR analysis. Upon entrance into nutrient starvation and after nitric oxide exposure, leaderless translation is significantly less affected by the stress than Shine-Dalgarno translation. Similarly, during the early stages of infection of macrophages, the levels of leaderless translation are transiently more stable than those of Shine-Dalgarno translation. These results suggest that leaderless translation may offer an advantage in the physiology of M. tuberculosis. Identification of the molecular mechanisms underlying this translational regulation may provide insights into persistent infection.

Light-Dependent Translation Change of Arabidopsis psbA Correlates with RNA Structure Alterations at the Translation Initiation Region

mRNA secondary structure influences translation. Proteins that modulate the mRNA secondary structure around the translation initiation region may regulate translation in plastids. To test this hypothesis, we exposed Arabidopsis thaliana to high light, which induces translation of psbA mRNA encoding the D1 subunit of photosystem II. We assayed translation by ribosome profiling and applied two complementary methods to analyze in vivo RNA secondary structure: DMS-MaPseq and SHAPE-seq. We detected increased accessibility of the translation initiation region of psbA after high light treatment, likely contributing to the observed increase in translation by facilitating translation initiation. Furthermore, we identified the footprint of a putative regulatory protein in the 5' UTR of psbA at a position where occlusion of the nucleotide sequence would cause the structure of the translation initiation region to open up, thereby facilitating ribosome access. Moreover, we show that other plastid genes with weak Shine-Dalgarno sequences (SD) are likely to exhibit psbA-like regulation, while those with strong SDs do not. This supports the idea that changes in mRNA secondary structure might represent a general mechanism for translational regulation of psbA and other plastid genes.

PPD: A Manually Curated Database for Experimentally Verified Prokaryotic Promoters

As a key region, promoter plays a key role in transcription regulation. A eukaryotic promoter database called EPD has been constructed to store eukaryotic POL II promoters. Although there are some promoter databases for specific prokaryotic species or specific promoter type, such as RegulonDB for Escherichia coli K-12, DBTBS for Bacillus subtilis and Pro54DB for sigma 54 promoter, because of the diversity of prokaryotes and the development of sequencing technology, huge amounts of prokaryotic promoters are scattered in numerous published articles, which is inconvenient for researchers to explore the process of gene regulation in prokaryotes. In this study, we constructed a Prokaryotic Promoter Database (PPD), which records the experimentally validated promoters in prokaryotes, from published articles. Up to now, PPD has stored 129,148 promoters across 63 prokaryotic species manually extracted from published papers. We provided a friendly interface for users to browse, search, blast, visualize, submit and download data. The PPD will provide relatively comprehensive resources of prokaryotic promoter for the study of prokaryotic gene transcription. The PPD is freely available and easy accessed at http://lin-group.cn/database/ppd/.

Structural basis of sequestration of the anti-Shine-Dalgarno sequence in the Bacteroidetes ribosome

Genomic studies have indicated that certain bacterial lineages such as the Bacteroidetes lack Shine-Dalgarno (SD) sequences, and yet with few exceptions ribosomes of these organisms carry the canonical anti-SD (ASD) sequence. Here, we show that ribosomes purified from Flavobacterium johnsoniae, a representative of the Bacteroidetes, fail to recognize the SD sequence of mRNA in vitro. A cryo-electron microscopy structure of the complete 70S ribosome from F. johnsoniae at 2.8 Å resolution reveals that the ASD is sequestered by ribosomal proteins bS21, bS18 and bS6, explaining the basis of ASD inhibition. The structure also uncovers a novel ribosomal protein—bL38. Remarkably, in F. johnsoniae and many other Flavobacteriia, the gene encoding bS21 contains a strong SD, unlike virtually all other genes. A subset of Flavobacteriia have an alternative ASD, and in these organisms the fully complementary sequence lies upstream of the bS21 gene, indicative of natural covariation. In other Bacteroidetes classes, strong SDs are frequently found upstream of the genes for bS21 and/or bS18. We propose that these SDs are used as regulatory elements, enabling bS21 and bS18 to translationally control their own production.

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.

Integrative characterization of the near-minimal bacterium Mesoplasma florum

The near-minimal bacterium Mesoplasma florum is an interesting model for synthetic genomics and systems biology due to its small genome (~ 800 kb), fast growth rate, and lack of pathogenic potential. However, fundamental aspects of its biology remain largely unexplored. Here, we report a broad yet remarkably detailed characterization of M. florum by combining a wide variety of experimental approaches. We investigated several physical and physiological parameters of this bacterium, including cell size, growth kinetics, and biomass composition of the cell. We also performed the first genome-wide analysis of its transcriptome and proteome, notably revealing a conserved promoter motif, the organization of transcription units, and the transcription and protein expression levels of all protein-coding sequences. We converted gene transcription and expression levels into absolute molecular abundances using biomass quantification results, generating an unprecedented view of the M. florum cellular composition and functions. These characterization efforts provide a strong experimental foundation for the development of a genome-scale model for M. florum and will guide future genome engineering endeavors in this simple organism.

Translational regulation of environmental adaptation in bacteria

Bacteria must rapidly respond to both intracellular and environmental changes to survive. One critical mechanism to rapidly detect and adapt to changes in environmental conditions is control of gene expression at the level of protein synthesis. At each of the three major steps of translation-initiation, elongation, and termination-cells use stimuli to tune translation rate and cellular protein concentrations. For example, changes in nutrient concentrations in the cell can lead to translational responses involving mechanisms such as dynamic folding of riboswitches during translation initiation or the synthesis of alarmones, which drastically alter cell physiology. Moreover, the cell can fine-tune the levels of specific protein products using programmed ribosome pausing or inducing frameshifting. Recent studies have improved understanding and revealed greater complexity regarding long-standing paradigms describing key regulatory steps of translation such as start-site selection and the coupling of transcription and translation. In this review, we describe how bacteria regulate their gene expression at the three translational steps and discuss how translation is used to detect and respond to changes in the cellular environment. Finally, we appraise the costs and benefits of regulation at the translational level in bacteria.

The Effect of Translation Promoting Site (TPS) on Protein Expression in E. coli Cells

The study of translation initiation in prokaryotes assumes that there should be a mechanism different from the canonical model, which postulates the formation of the pre-initiation complex through the interaction of the Shine-Dalgarno sequence (SD) at the 5'-end of mRNA and the anti-Shine-Dalgarno site at the 3'-end of 16S rRNA. In this paper we've studied the effect of TPS (Translation-initiation Promoting Site) on β-glucuronidase expression in E. coli cells at different cultivation temperatures. The examined leader sequences were cloned into the pET23c plasmid upstream the β-glucuronidase gene; protein expression was performed in E. coli BL21 (DE3) cells. β-glucuronidase activity was measured in bacterial cell extracts via paranitrophenyl b-D-glucuronide assay. The quantity of expressed protein was measured by Western blotting with following densitometry. It was shown that TPS increases the level of protein expression at stressful conditions (10 °C and 44 °C) 5-8 times compared to control. The combination of TPS and SD sites in the 5'-leader sequence of the mRNA created an enhancer that increased the expression level 2-3.6 times compared to a single SD-sequence. Based on the obtained data and the computer modeling of interaction between 16S rRNA and TPS, we proposed an alternative variation of prokaryotic translation initiation.

Unique Shine-Dalgarno Sequences in Cyanobacteria and Chloroplasts Reveal Evolutionary Differences in Their Translation Initiation

Microorganisms require efficient translation to grow and replicate rapidly, and translation is often rate-limited by initiation. A prominent feature that facilitates translation initiation in bacteria is the Shine-Dalgarno (SD) sequence. However, there is much debate over its conservation in Cyanobacteria and in chloroplasts which presumably originated from endosymbiosis of ancient Cyanobacteria. Elucidating the utilization of SD sequences in Cyanobacteria and in chloroplasts is therefore important to understand whether 1) SD role in Cyanobacterial translation has been reduced prior to chloroplast endosymbiosis or 2) translation in Cyanobacteria and in plastid has been subjected to different evolutionary pressures. To test these alternatives, we employed genomic, proteomic, and transcriptomic data to trace differences in SD usage among Synechocystis species, Microcystis aeruginosa, cyanophages, Nicotiana tabacum chloroplast, and Arabidopsis thaliana chloroplast. We corrected their mis-annotated 16S rRNA 3' terminus using an RNA-Seq-based approach to determine their SD/anti-SD locational constraints using an improved measurement DtoStart. We found that cyanophages well-mimic Cyanobacteria in SD usage because both have been under the same selection pressure for SD-mediated initiation. Whereas chloroplasts lost this similarity because the need for SD-facilitated initiation has been reduced in plastids having much reduced genome size and different ribosomal proteins as a result of host-symbiont coevolution. Consequently, SD sequence significantly increases protein expression in Cyanobacteria but not in chloroplasts, and only Cyanobacterial genes compensate for a lack of SD sequence by having weaker secondary structures at the 5' UTR. Our results suggest different evolutionary pressures operate on translation initiation in Cyanobacteria and in chloroplast.

Global analysis of protein synthesis in Flavobacterium johnsoniae reveals the use of Kozak-like sequences in diverse bacteria

In all cells, initiation of translation is tuned by intrinsic features of the mRNA. Here, we analyze translation in Flavobacterium johnsoniae, a representative of the Bacteroidetes. Members of this phylum naturally lack Shine-Dalgarno (SD) sequences in their mRNA, and yet their ribosomes retain the conserved anti-SD sequence. Translation initiation is tuned by mRNA secondary structure and by the identities of several key nucleotides upstream of the start codon. Positive determinants include adenine at position -3, reminiscent of the Kozak sequence of Eukarya. Comparative analysis of Escherichia coli reveals use of the same Kozak-like sequence to enhance initiation, suggesting an ancient and widespread mechanism. Elimination of contacts between A-3 and the conserved β-hairpin of ribosomal protein uS7 fails to diminish the contribution of A-3 to initiation, suggesting an indirect mode of recognition. Also, we find that, in the Bacteroidetes, the trinucleotide AUG is underrepresented in the vicinity of the start codon, which presumably helps compensate for the absence of SD sequences in these organisms.

AssessORF: combining evolutionary conservation and proteomics to assess prokaryotic gene predictions

MOTIVATION

A core task of genomics is to identify the boundaries of protein coding genes, which may cover over 90% of a prokaryote's genome. Several programs are available for gene finding, yet it is currently unclear how well these programs perform and whether any offers superior accuracy. This is in part because there is no universal benchmark for gene finding and, therefore, most developers select their own benchmarking strategy.

RESULTS

Here, we introduce AssessORF, a new approach for benchmarking prokaryotic gene predictions based on evidence from proteomics data and the evolutionary conservation of start and stop codons. We applied AssessORF to compare gene predictions offered by GenBank, GeneMarkS-2, Glimmer and Prodigal on genomes spanning the prokaryotic tree of life. Gene predictions were 88-95% in agreement with the available evidence, with Glimmer performing the worst but no clear winner. All programs were biased towards selecting start codons that were upstream of the actual start. Given these findings, there remains considerable room for improvement, especially in the detection of correct start sites.

AVAILABILITY AND IMPLEMENTATION

AssessORF is available as an R package via the Bioconductor package repository.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Identification of common highly expressed genes of Salmonella Enteritidis by in silico prediction of gene expression and in vitro transcriptomic analysis

Chickens are the reservoir host of Salmonella Enteritidis. Salmonella Enteritidis colonizes the gastro-intestinal tract of chickens and replicates within macrophages without causing clinically discernable illness. Persistence of S. Enteritidis in the hostile environments of intestinal tract and macrophages allows it to disseminate extra-intestinally to liver, spleen, and reproductive tract. Extra-intestinal dissemination into reproductive tract leads to contamination of internal contents of eggs, which is a major risk factor for human infection. Understanding the genes that contribute to S. Enteritidis persistence in the chicken host is central to elucidate the genetic basis of the unique pathobiology of this public health pathogen. The aim of this study was to identify a succinct set of genes associated with infection-relevant in vitro environments to provide a rational foundation for subsequent biologically-relevant research. We used in silico prediction of gene expression and RNA-seq technology to identify a core set of 73 S. Enteritidis genes that are consistently highly expressed in multiple S. Enteritidis strains cultured at avian physiologic temperature under conditions that represent intestinal and intracellular environments. These common highly expressed (CHX) genes encode proteins involved in bacterial metabolism, protein synthesis, cell-envelope biogenesis, stress response, and a few proteins with uncharacterized functions. Further studies are needed to dissect the contribution of these CHX genes to the pathobiology of S. Enteritidis in the avian host. Several of the CHX genes could serve as promising targets for studies towards the development of immunoprophylactic and novel therapeutic strategies to prevent colonization of chickens and their environment with S. Enteritidis.

Control of Translation at the Initiation Phase During Glucose Starvation in Yeast

Glucose is one of the most important sources of carbon across all life. Glucose starvation is a key stress relevant to all eukaryotic cells. Glucose starvation responses have important implications in diseases, such as diabetes and cancer. In yeast, glucose starvation causes rapid and dramatic effects on the synthesis of proteins (mRNA translation). Response to glucose deficiency targets the initiation phase of translation by different mechanisms and with diverse dynamics. Concomitantly, translationally repressed mRNAs and components of the protein synthesis machinery may enter a variety of cytoplasmic foci, which also form with variable kinetics and may store or degrade mRNA. Much progress has been made in understanding these processes in the last decade, including with the use of high-throughput/omics methods of RNA and RNA:protein detection. This review dissects the current knowledge of yeast reactions to glucose starvation systematized by the stage of translation initiation, with the focus on rapid responses. We provide parallels to mechanisms found in higher eukaryotes, such as metazoans, for the most critical responses, and point out major remaining gaps in knowledge and possible future directions of research on translational responses to glucose starvation.

RNA-Seq-Based Analysis Reveals Heterogeneity in Mature 16S rRNA 3' Termini and Extended Anti-Shine-Dalgarno Motifs in Bacterial Species

We present an RNA-Seq based approach to map 3′ end sequences of mature 16S rRNA (3′ TAIL) in bacteria with single-base specificity. Our results show that 3′ TAILs are heterogeneous among species; they contain the core CCUCC anti-Shine-Dalgarno motif, but vary in downstream lengths. Importantly, our findings rectify the mis-annotated 16S rRNAs in 11 out of 13 bacterial species studied herein (covering Cyanobacteria, Deinococcus-Thermus, Firmicutes, Proteobacteria, Tenericutes, and Spirochaetes). Furthermore, our results show that species-specific 3′ TAIL boundaries are retained due to their high complementarity with preferred Shine-Dalgarno sequences, suggesting that 3′ TAIL bases downstream of the canonical CCUCC motif play a more important role in translation initiation than previously reported.

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

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

Comparative genetic and genomic analysis of the novel fusellovirus Sulfolobus spindle-shaped virus 10

Viruses that infect thermophilic Archaea are unique in both their structure and genetic makeup. The lemon-shaped fuselloviruses—which infect members of the order Sulfolobales, growing optimally at 80 °C and pH 3—are some of the most ubiquitous and best studied viruses of the thermoacidophilic Archaea. Nonetheless, much remains to be learned about these viruses. In order to investigate fusellovirus evolution, we have isolated and characterized a novel fusellovirus, Sulfolobus spindle-shaped virus 10 (formerly SSV-L1). Comparative genomic analyses highlight significant similarity with both SSV8 and SSV9, as well as conservation of promoter elements within the Fuselloviridae. SSV10 encodes five ORFs with no homology within or outside of the Fuselloviridae, as well as a putatively functional Cas4-like ORF, which may play a role in evading CRISPR-mediated host defenses. Moreover, we demonstrate the ability of SSV10 to withstand mutation in a fashion consistent with mutagenesis in SSV1.

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

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

Chloroplast Translation: Structural and Functional Organization, Operational Control, and Regulation

Chloroplast translation is essential for cellular viability and plant development. Its positioning at the intersection of organellar RNA and protein metabolism makes it a unique point for the regulation of gene expression in response to internal and external cues. Recently obtained high-resolution structures of plastid ribosomes, the development of approaches allowing genome-wide analyses of chloroplast translation (i.e., ribosome profiling), and the discovery of RNA binding proteins involved in the control of translational activity have greatly increased our understanding of the chloroplast translation process and its regulation. In this review, we provide an overview of the current knowledge of the chloroplast translation machinery, its structure, organization, and function. In addition, we summarize the techniques that are currently available to study chloroplast translation and describe how translational activity is controlled and which cis-elements and trans-factors are involved. Finally, we discuss how translational control contributes to the regulation of chloroplast gene expression in response to developmental, environmental, and physiological cues. We also illustrate the commonalities and the differences between the chloroplast and bacterial translation machineries and the mechanisms of protein biosynthesis in these two prokaryotic systems.

Shine-Dalgarno Sequences Play an Essential Role in the Translation of Plastid mRNAs in Tobacco

In prokaryotic systems, the translation initiation of many, though not all, mRNAs depends on interaction between a sequence element upstream of the start codon (the Shine-Dalgarno sequence [SD]) and a complementary sequence in the 3' end of the 16S rRNA (anti-Shine-Dalgarno sequence [aSD]). Although many chloroplast mRNAs harbor putative SDs in their 5' untranslated regions and the aSD displays strong conservation, the functional relevance of SD-aSD interactions in plastid translation is unclear. Here, by generating transplastomic tobacco (Nicotiana tabacum) mutants with point mutations in the aSD coupled with genome-wide analysis of translation by ribosome profiling, we provide a global picture of SD-dependent translation in plastids. We observed a pronounced correlation between weakened predicted SD-aSD interactions and reduced translation efficiency. However, multiple lines of evidence suggest that the strength of the SD-aSD interaction is not the only determinant of the translational output of many plastid mRNAs. Finally, the translation efficiency of mRNAs with strong secondary structures around the start codon is more dependent on the SD-aSD interaction than weakly structured mRNAs. Thus, our data reveal the importance of the aSD in plastid translation initiation, uncover chloroplast genes whose translation is influenced by SD-aSD interactions, and provide insights into determinants of translation efficiency in plastids.