Re-annotation of 12,495 prokaryotic 16S rRNA 3' ends and analysis of Shine-Dalgarno and anti-Shine-Dalgarno 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.

Contrasting Genetic Diversity of Listeria Pathogenicity Islands 3 and 4 Harbored by Nonpathogenic Listeria spp

Listeria monocytogenes causes the severe foodborne disease listeriosis. Several clonal groups of L. monocytogenes possess the pathogenicity islands Listeria pathogenicity island 3 (LIPI-3) and LIPI-4. Here, we investigated the prevalence and genetic diversity of LIPI-3 and LIPI-4 among 63 strains of seven nonpathogenic Listeria spp. from the natural environment, i.e., wildlife (black bears [Ursus americanus]) and surface water. Analysis of the whole-genome sequence data suggested that both islands were horizontally acquired but differed considerably in their incidence and genetic diversity. LIPI-3 was identified among half of the L. innocua strains in the same genomic location as in L. monocytogenes (guaA hot spot) in a truncated form, with only three strains harboring full-length LIPI-3, and a highly divergent partial LIPI-3 was observed in three Listeria seeligeri strains, outside the guaA hot spot. Premature stop codons (PMSCs) and frameshifts were frequently noted in the LIPI-3 gene encoding listeriolysin S. On the other hand, full-length LIPI-4 without any PMSCs was found in all Listeria innocua strains, in the same genomic location as L. monocytogenes and with ~85% similarity to the L. monocytogenes counterpart. Our study provides intriguing examples of genetic changes that pathogenicity islands may undergo in nonpathogenic bacterial species, potentially in response to environmental pressures that promote either maintenance or degeneration of the islands. Investigations of the roles that LIPI-3 and LIPI-4 play in nonpathogenic Listeria spp. are warranted to further understand the differential evolution of genetic elements in pathogenic versus nonpathogenic hosts of the same genus. IMPORTANCE Listeria monocytogenes is a serious foodborne pathogen that can harbor the pathogenicity islands Listeria pathogenicity island 3 (LIPI-3) and LIPI-4. Intriguingly, these have also been reported in nonpathogenic L. innocua from food and farm environments, though limited information is available for strains from the natural environment. Here, we analyzed whole-genome sequence data of nonpathogenic Listeria spp. from wildlife and surface water to further elucidate the genetic diversity and evolution of LIPI-3 and LIPI-4 in Listeria. While the full-length islands were found only in L. innocua, LIPI-3 was uncommon and exhibited frequent truncation and genetic diversification, while LIPI-4 was remarkable in being ubiquitous, albeit diversified from L. monocytogenes. These contrasting features demonstrate that pathogenicity islands in nonpathogenic hosts can evolve along different trajectories, leading to either degeneration or maintenance, and highlight the need to examine their physiological roles in nonpathogenic hosts.

Three Phages One Host: Isolation and Characterization of Pantoea agglomerans Phages from a Grasshopper Specimen

The bacterial genus Pantoea comprises species found in a variety of different environmental sources. Pantoea spp. are often recovered from plant material and are capable of both benefitting the plants and acting like phytopathogens. Some species of Pantoea (including P. agglomerans) are considered opportunistic human pathogens capable of causing various infections in immunocompromised subjects. In this study, a strain of P. agglomerans (identified by 16S rRNA gene sequencing) was isolated from a dead specimen of an unidentified Latvian grasshopper species. The retrieved strain of P. agglomerans was then used as a host for the potential retrieval of phages from the same source material. After rounds of plaque purification and propagation, three high-titer lysates corresponding to putatively distinct phages were acquired. Transmission electron microscopy revealed that one of the phages was a myophage with an unusual morphology, while the two others were typical podophages. Whole-genome sequencing (WGS) was performed for each of these isolated phages. Genome de novo assembly and subsequent functional annotation confirmed that three different strictly lytic phages were isolated. Elaborate genomic characterization of the acquired phages was performed to elucidate their place within the so-far-uncovered phage diversity.

Morganella Phage Mecenats66 Utilizes an Evolutionarily Distinct Subtype of Headful Genome Packaging with a Preferred Packaging Initiation Site

Both recognized species from the genus Morganella (M. morganii and M. psychrotolerans) are Gram-negative facultative anaerobic rod-shaped bacteria that have been documented as sometimes being implicated in human disease. Complete genomes of seven Morganella-infecting phages are publicly available today. Here, we report on the genomic characterization of an insect associated Morganella sp. phage, which we named Mecenats66, isolated from dead worker honeybees. Phage Mecenats66 was propagated, purified, and subjected to whole-genome sequencing with subsequent complete genome annotation. After the genome de novo assembly, it was noted that Mecenats66 might employ a headful packaging with a preferred packaging initiation site, although its terminase amino acid sequence did not fall within any of the currently recognized headful packaging subtype employing phage (that had their packaging strategy experimentally verified) with clusters on a terminase sequence phylogenetic tree. The in silico predicted packaging strategy was verified experimentally, validating the packaging initiation site and suggesting that Mecenats66 represents an evolutionarily distinct headful genome packaging with a preferred packaging initiation site strategy subtype. These findings can possibly be attributed to several of the phages already found within the public biological sequence repositories and could aid newly isolated phage packaging strategy predictions in the future.

Design of 5′-UTR to Enhance Keratinase Activity in Bacillus subtilis

Keratinase is an important industrial enzyme, but its application performance is limited by its low activity. A rational design of 5′-UTRs that increases translation efficiency is an important approach to enhance protein expression. Herein, we optimized the 5′-UTR of the recombinant keratinase KerZ1 expression element to enhance its secretory activity in Bacillus subtilis WB600 through Spacer design, RBS screening, and sequence simplification. First, the A/U content in Spacer was increased by the site-directed saturation mutation of G/C bases, and the activity of keratinase secreted by mutant strain B. subtilis WB600-SP was 7.94 times higher than that of KerZ1. Subsequently, the keratinase activity secreted by the mutant strain B. subtilis WB600-SP-R was further increased to 13.45 times that of KerZ1 based on the prediction of RBS translation efficiency and the multi-site saturation mutation screening. Finally, the keratinase activity secreted by the mutant strain B. subtilis WB600-SP-R-D reached 204.44 KU mL−1 by reducing the length of the 5′ end of the 5′-UTR, which was 19.70 times that of KerZ1. In a 5 L fermenter, the keratinase activity secreted by B. subtilis WB600-SP-R-D after 25 h fermentation was 797.05 KU mL−1, which indicated its high production intensity. Overall, the strategy of this study and the obtained keratinase mutants will provide a good reference for the expression regulation of keratinase and other industrial enzymes.

Features of smaller ribosomes in candidate phyla radiation (CPR) bacteria revealed with a molecular evolutionary analysis

The candidate phyla radiation (CPR) is a large bacterial group consisting mainly of uncultured lineages. They have small cells and small genomes, and they often lack ribosomal proteins uL1, bL9, and/or uL30, which are basically ubiquitous in non-CPR bacteria. Here, we comprehensively analyzed the genomic information on CPR bacteria and identified their unique properties. The distribution of protein lengths in CPR bacteria peaks at around 100-150 amino acids, whereas the position of the peak varies in the range of 100-300 amino acids in free-living non-CPR bacteria, and at around 100-200 amino acids in most symbiotic non-CPR bacteria. These results show that the proteins of CPR bacteria are smaller, on average, than those of free-living non-CPR bacteria, like those of symbiotic non-CPR bacteria. We found that ribosomal proteins bL28, uL29, bL32, and bL33 have been lost in CPR bacteria in a taxonomic lineage-specific manner. Moreover, the sequences of approximately half of all ribosomal proteins of CPR differ, in part, from those of non-CPR bacteria, with missing regions or specifically added regions. We also found that several regions in the 16S, 23S, and 5S rRNAs of CPR bacteria are lacking, which presumably caused the total predicted lengths of the three rRNAs of CPR bacteria to be smaller than those of non-CPR bacteria. The regions missing in the CPR ribosomal proteins and rRNAs are located near the surface of the ribosome, and some are close to one another. These observations suggest that ribosomes are smaller in CPR bacteria than those in free-living non-CPR bacteria, with simplified surface structures.

RiboGrove: a database of full-length prokaryotic 16S rRNA genes derived from completely assembled genomes

16S rRNA gene is frequently used for the identification of prokaryotic organisms and for phylogeny inference. Several specialized public databases exist that contain complete and partial sequences of 16S rRNA genes. In this paper, we present RiboGrove: the first publicly available database that comprises only full-length sequences of 16S rRNA genes originating from completely assembled prokaryotic genomes deposited in RefSeq. Despite being strongly biased towards frequently sequenced genomes, RiboGrove is a useful complement to existing 16S rRNA resources and allows for analyses that would not be possible using amplicon-derived gene sequences. For instance, the absence of partial gene sequences in RiboGrove allowed us to make a summary of prokaryotic organisms, which lack core anti-Shine-Dalgarno sequence in their 16S rRNA genes. In this study, we describe the collected sequence data and present the results of exploratory data analysis of 16S rRNA gene sequences.

Genome Characterization of Nocturne116, Novel Lactococcus lactis-Infecting Phage Isolated from Moth

While looking for novel insect-associated phages, a unique siphophage, Nocturne116, was isolated from a deceased local moth specimen along with its host, which was identified by 16S rRNA gene sequencing as a strain of Lactococcus lactis. Next-generation sequencing and the subsequent genome annotation elaborated on herein revealed that the genome of Nocturne116 is a 25,554 bp long dsDNA molecule with 10 bp long 3' cos overhangs and a GC content of 37.99%, comprising 52 predicted open reading frames. The complete nucleotide sequence of phage Nocturne116 genome is dissimilar to any of the already sequenced phages, save for a distant link with Lactococcus phage Q54. Functions for only 15/52 of Nocturne116 gene products could be reliably predicted using contemporary comparative genomics approaches, while 22 of its gene products do not yet have any homologous entries in the public biological sequence repositories. Despite the public availability of nearly 350 elucidated Lactococcus phage complete genomes as of now, Nocturne116 firmly stands out as a sole representative of novel phage genus.

Simplification of Ribosomes in Bacteria with Tiny Genomes

The ribosome is an essential cellular machine performing protein biosynthesis. Its structure and composition are highly conserved in all species. However, some bacteria have been reported to have an incomplete set of ribosomal proteins. We have analyzed ribosomal protein composition in 214 small bacterial genomes (<1 Mb) and found that although the ribosome composition is fairly stable, some ribosomal proteins may be absent, especially in bacteria with dramatically reduced genomes. The protein composition of the large subunit is less conserved than that of the small subunit. We have identified the set of frequently lost ribosomal proteins and demonstrated that they tend to be positioned on the ribosome surface and have fewer contacts to other ribosome components. Moreover, some proteins are lost in an evolutionary correlated manner. The reduction of ribosomal RNA is also common, with deletions mostly occurring in free loops. Finally, the loss of the anti-Shine-Dalgarno sequence is associated with the loss of a higher number of ribosomal proteins.

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.

Simplification of Ribosomes in Bacteria with Tiny Genomes

The ribosome is an essential cellular machine performing protein biosynthesis. Its structure and composition are highly conserved in all species. However, some bacteria have been reported to have an incomplete set of ribosomal proteins. We have analyzed ribosomal protein composition in 214 small bacterial genomes (<1 Mb) and found that although the ribosome composition is fairly stable, some ribosomal proteins may be absent, especially in bacteria with dramatically reduced genomes. The protein composition of the large subunit is less conserved than that of the small subunit. We have identified the set of frequently lost ribosomal proteins and demonstrated that they tend to be positioned on the ribosome surface and have fewer contacts to other ribosome components. Moreover, some proteins are lost in an evolutionary correlated manner. The reduction of ribosomal RNA is also common, with deletions mostly occurring in free loops. Finally, the loss of the anti-Shine-Dalgarno sequence is associated with the loss of a higher number of ribosomal proteins.

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.

Species-specific genomic sequences for classification of bacteria

Modern bacterial classification relies on genomic relatedness. Genetic variation in bacterial populations present a big challenge for taxonomic classification and recently several bacterial species have been reclassified based on the intra-species genome comparison. These were facilitated by next generation sequencing technologies and advances in genome comparison approaches which led to the rearrangement of diverse bacterial species and revolution in the microbial classification system. One of the outcome of these studies is the development of suitable DNA barcodes as reliable and cost-effective method for identifying various bacterial genera. Towards refining this further, we have applied a genome comparison approach in 1104 bacterial genome assemblies (excluding plasmids) to identify unique genomic segments among intra-species genome assemblies. Using extensive bioinformatics analysis, we have identified species-specific genomic regions and designed unique primers for 100 different species (belonging to 62 genera) which includes 62 pathogenic and 13 opportunistic pathogenic bacterial species and built a database (http://slsdb.manipal.edu/bact/). These species-specific genomic regions will have a major impact on in silico and molecular methods aimed at bacterial classification and identification. These may also serve as better DNA barcodes than the markers currently used for delineation of bacteria and may also find application in various translational research programs.

Recommendations for bacterial ribosome profiling experiments based on bioinformatic evaluation of published data

Ribosome profiling (RIBO-Seq) has improved our understanding of bacterial translation, including finding many unannotated genes. However, protocols for RIBO-Seq and corresponding data analysis are not yet standardized. Here, we analyzed 48 RIBO-Seq samples from nine studies of Escherichia coli K12 grown in lysogeny broth medium and particularly focused on the size-selection step. We show that for conventional expression analysis, a size range between 22 and 30 nucleotides is sufficient to obtain protein-coding fragments, which has the advantage of removing many unwanted rRNA and tRNA reads. More specific analyses may require longer reads and a corresponding improvement in rRNA/tRNA depletion. There is no consensus about the appropriate sequencing depth for RIBO-Seq experiments in prokaryotes, and studies vary significantly in total read number. Our analysis suggests that 20 million reads that are not mapping to rRNA/tRNA are required for global detection of translated annotated genes. We also highlight the influence of drug-induced ribosome stalling, which causes bias at translation start sites. The resulting accumulation of reads at the start site may be especially useful for detecting weakly expressed genes. As different methods suit different questions, it may not be possible to produce a "one-size-fits-all" ribosome profiling data set. Therefore, experiments should be carefully designed in light of the scientific questions of interest. We propose some basic characteristics that should be reported with any new RIBO-Seq data sets. Careful attention to the factors discussed should improve prokaryotic gene detection and the comparability of ribosome profiling data sets.

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.

Identification of a circular code periodicity in the bacterial ribosome: origin of codon periodicity in genes?

Three-base periodicity (TBP), where nucleotides and higher order n-tuples are preferentially spaced by 3, 6, 9, etc. bases, is a well-known intrinsic property of protein-coding DNA sequences. However, its origins are still not fully understood. One hypothesis is that the periodicity reflects a primordial coding system that was used before the emergence of the modern standard genetic code (SGC). Recent evidence suggests that the X circular code, a set of 20 trinucleotides allowing the reading frames in genes to be retrieved locally, represents a possible ancestor of the SGC. Motifs from the X circular code have been found in the reading frame of protein-coding regions in extant organisms from bacteria to eukaryotes, in many transfer RNA (tRNA) genes and in important functional regions of the ribosomal RNA (rRNA), notably in the peptidyl transferase centre and the decoding centre. Here, we have used a powerful correlation function to search for periodicity patterns involving the 20 trinucleotides of the X circular code in a large set of bacterial protein-coding genes, as well as in the translation machinery, including rRNA and tRNA sequences. As might be expected, we found a strong circular code periodicity 0 modulo 3 in the protein-coding genes. More surprisingly, we also identified a similar circular code periodicity in a large region of the 16S rRNA. This region includes the 3' major domain corresponding to the primordial proto-ribosome decoding centre and containing numerous sites that interact with the tRNA and messenger RNA (mRNA) during translation. Furthermore, 3D structural analysis shows that the periodicity region surrounds the mRNA channel that lies between the head and the body of the SSU. Our results support the hypothesis that the X circular code may constitute an ancestral translation code involved in reading frame retrieval and maintenance, traces of which persist in modern mRNA, tRNA and rRNA despite their long evolution and adaptation to the SGC.

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.

Dengue and Zika Virus 5' Untranslated Regions Harbor Internal Ribosomal Entry Site Functions

Members of the genus Flavivirus of Flaviviridae are important human pathogens of great concern because they cause serious diseases, sometimes death, in human populations living in tropical, subtropical (dengue virus [DENV], Zika virus [ZIKV], and yellow fever virus), or moderate climates (West Nile virus). Flaviviruses are known to control their translation by a cap-dependent mechanism. We have observed, however, that the uncapped genomes of DENV or ZIKV can initiate infection of mammalian and insect cells. We provide evidence that the short 5′ untranslated region (5′-UTR) of DENV or ZIKV genomes can fulfill the function of an internal ribosomal entry site (IRES). This strategy frees these organisms from the cap-dependent mechanism of gene expression at an as yet unknown stage of proliferation. The data raise new questions about the biology and evolution of flaviviruses, possibly leading to new controls of flavivirus disease.

The Flavivirus genus of the Flaviviridae family encompasses numerous enveloped plus-strand RNA viruses. Dengue virus (DENV), a flavivirus, is the leading cause of serious arthropod-borne disease globally. The genomes of DENV, like the genomes of yellow fever virus (YFV), West Nile fever virus (WNV), or Zika virus (ZIKV), control their translation by a 5′-terminal capping group. Three other genera of Flaviviridae are remarkable because their viruses use internal ribosomal entry sites (IRESs) to control translation, and they are not arthropod transmitted. In 2006, E. Harris’ group published work suggesting that DENV RNA does not stringently need a cap for translation. They proposed that instead DENV translation is controlled by an interplay between 5′ and 3′ termini. Here we present evidence that the DENV or ZIKV 5′ untranslated regions (5′-UTRs) alone have IRES competence. This conclusion is based, first, on the observation that uncapped monocistronic mRNAs 5′ terminated with the DENV or ZIKV 5′-UTRs can efficiently direct translation of a reporter gene in BHK and C6/36 cells and second, that either 5′-UTR placed between two reporter genes can efficiently induce expression of the downstream gene in BHK cells but not in C6/36 cells. These experiments followed observations that uncapped DENV/ZIKV genomic transcripts, 5′ terminated with pppAN… or GpppAN…, can initiate infections of mammalian (BHK) or mosquito (C6/36) cells. IRES competence of the 5′-UTRs of DENV/ZIKV raises many open questions regarding the biology and control, as well as the evolution, of insect-borne flaviviruses.