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

Tracing the serendipitous genesis of radiation resistance

Free-living organisms frequently encounter unfavorable abiotic environmental factors. Those who adapt and cope with sudden changes in the external environment survive. Desiccation is one of the most common and frequently encountered stresses in nature. On the contrary, ionizing radiations are limited to high local concentrations of naturally occurring radioactive materials and related anthropogenic activities. Yet, resistance to high doses of ionizing radiation is evident across the tree of life. The evolution of desiccation resistance has been linked to the evolution of ionizing radiation resistance, although, evidence to support the idea that the evolution of desiccation tolerance is a necessary precursor to ionizing radiation resistance is lacking. Moreover, the presence of radioresistance in hyperthermophiles suggests multiple paths lead to radiation resistance. In this minireview, we focus on the molecular aspects of damage dynamics and damage response pathways comprising protective and restorative functions with a definitive survival advantage, to explore the serendipitous genesis of ionizing radiation resistance.

DNA repair and oxidative stress defense systems in radiation-resistant Deinococcus murrayi

Deinococcus murrayi is a bacterium isolated from hot springs in Portugal, and named after Dr. Robert G.E. Murray in recognition of his research on the genus Deinococcus. Like other Deinococcus species, D. murrayi is extremely resistant to ionizing radiation. Repair of massive DNA damage and limitation of oxidative protein damage are two important factors contributing to the robustness of Deinococcus bacteria. Here, we identify, among others, the DNA repair and oxidative stress defense proteins in D. murrayi, and highlight special features of D. murrayi. For DNA repair, D. murrayi does not contain a standalone uracil-DNA glycosylase (Ung), but it encodes a protein in which Ung is fused to a DNA photolyase domain (PhrB). UvrB and UvrD contain large insertions corresponding to inteins. One of its endonuclease III enzymes lacks a [4Fe-4S] cluster. Deinococcus murrayi possesses a homolog of the error-prone DNA polymerase IV. Concerning oxidative stress defense, D. murrayi encodes a manganese catalase in addition to a heme catalase. Its organic hydroperoxide resistance protein Ohr is atypical because the redox active cysteines are present in a CXXC motif. These and other characteristics of D. murrayi show further diversity among Deinococcus bacteria with respect to resistance-associated mechanisms.

Coping with stress: How bacteria fine-tune protein synthesis and protein transport

Maintaining a functional proteome under different environmental conditions is challenging for every organism, in particular for unicellular organisms, such as bacteria. In order to cope with changing environments and stress conditions, bacteria depend on strictly coordinated proteostasis networks that control protein production, folding, trafficking, and degradation. Regulation of ribosome biogenesis and protein synthesis are cornerstones of this cellular adaptation in all domains of life, which is rationalized by the high energy demand of both processes and the increased resistance of translationally silent cells against internal or external poisons. Reduced protein synthesis ultimately also reduces the substrate load for protein transport systems, which are required for maintaining the periplasmic, inner, and outer membrane subproteomes. Consequences of impaired protein transport have been analyzed in several studies and generally induce a multifaceted response that includes the upregulation of chaperones and proteases and the simultaneous downregulation of protein synthesis. In contrast, generally less is known on how bacteria adjust the protein targeting and transport machineries to reduced protein synthesis, e.g., when cells encounter stress conditions or face nutrient deprivation. In the current review, which is mainly focused on studies using Escherichia coli as a model organism, we summarize basic concepts on how ribosome biogenesis and activity are regulated under stress conditions. In addition, we highlight some recent developments on how stress conditions directly impair protein targeting to the bacterial membrane. Finally, we describe mechanisms that allow bacteria to maintain the transport of stress-responsive proteins under conditions when the canonical protein targeting pathways are impaired.

Mycobacterium smegmatis: The Vanguard of Mycobacterial Research

The genus Mycobacterium contains several slow-growing human pathogens, including Mycobacterium tuberculosis, Mycobacterium leprae, and Mycobacterium avium. Mycobacterium smegmatis is a nonpathogenic and fast growing species within this genus. In 1990, a mutant of M. smegmatis, designated mc²155, that could be transformed with episomal plasmids was isolated, elevating M. smegmatis to model status as the ideal surrogate for mycobacterial research. Classical bacterial models, such as Escherichia coli, were inadequate for mycobacteria research because they have low genetic conservation, different physiology, and lack the novel envelope structure that distinguishes the Mycobacterium genus. By contrast, M. smegmatis encodes thousands of conserved mycobacterial gene orthologs and has the same cell architecture and physiology. Dissection and characterization of conserved genes, structures, and processes in genetically tractable M. smegmatis mc²155 have since provided previously unattainable insights on these same features in its slow-growing relatives. Notably, tuberculosis (TB) drugs, including the first-line drugs isoniazid and ethambutol, are active against M. smegmatis, but not against E. coli, allowing the identification of their physiological targets. Furthermore, Bedaquiline, the first new TB drug in 40 years, was discovered through an M. smegmatis screen. M. smegmatis has become a model bacterium, not only for M. tuberculosis, but for all other Mycobacterium species and related genera. With a repertoire of bioinformatic and physical resources, including the recently established Mycobacterial Systems Resource, M. smegmatis will continue to accelerate mycobacterial research and advance the field of microbiology.

Regulation of Leaderless mRNA Translation in Bacteria

In bacteria, the translation of genetic information can begin through at least three different mechanisms: canonical or Shine-Dalgarno-led initiation, readthrough or 70S scanning initiation, or leaderless initiation. Here, we discuss the main features and regulation of the last, which is characterized mainly by the ability of 70S ribosomal particles to bind to AUG located at or near the 5′ end of mRNAs to initiate translation. These leaderless mRNAs (lmRNAs) are rare in enterobacteria, such as Escherichia coli, but are common in other bacteria, such as Mycobacterium tuberculosis and Deinococcus deserti, where they may represent more than 20% and even up to 60% of the genes. Given that lmRNAs are devoid of a 5′ untranslated region and the Shine-Dalgarno sequence located within it, the mechanism of translation regulation must depend on molecular strategies that are different from what has been observed in the Shine-Dalgarno-led translation. Diverse regulatory mechanisms have been proposed, including the processing of ribosomal RNA and changes in the abundance of translation factors, but all of them produce global changes in the initiation of lmRNA translation. Thus, further research will be required to understand how the initiation of the translation of particular lmRNA genes is regulated.

Thiol Reductases in Deinococcus Bacteria and Roles in Stress Tolerance

Deinococcus species possess remarkable tolerance to extreme environmental conditions that generate oxidative damage to macromolecules. Among enzymes fulfilling key functions in metabolism regulation and stress responses, thiol reductases (TRs) harbour catalytic cysteines modulating the redox status of Cys and Met in partner proteins. We present here a detailed description of Deinococcus TRs regarding gene occurrence, sequence features, and physiological functions that remain poorly characterised in this genus. Two NADPH-dependent thiol-based systems are present in Deinococcus. One involves thioredoxins, disulfide reductases providing electrons to protein partners involved notably in peroxide scavenging or in preserving protein redox status. The other is based on bacillithiol, a low-molecular-weight redox molecule, and bacilliredoxin, which together protect Cys residues against overoxidation. Deinococcus species possess various types of thiol peroxidases whose electron supply depends either on NADPH via thioredoxins or on NADH via lipoylated proteins. Recent data gained on deletion mutants confirmed the importance of TRs in Deinococcus tolerance to oxidative treatments, but additional investigations are needed to delineate the redox network in which they operate, and their precise physiological roles. The large palette of Deinococcus TR representatives very likely constitutes an asset for the maintenance of redox homeostasis in harsh stress conditions.

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.

A small RNA regulates pprM, a modulator of pleiotropic proteins promoting DNA repair, in Deinococcus radiodurans under ionizing radiation

Networks of transcriptional and post-transcriptional regulators are critical for bacterial survival and adaptation to environmental stressors. While transcriptional regulators provide rapid activation and/or repression of a wide-network of genes, post-transcriptional regulators, such as small RNAs (sRNAs), are also important to fine-tune gene expression. However, the mechanisms of sRNAs remain poorly understood, especially in less-studied bacteria. Deinococcus radiodurans is a gram-positive bacterium resistant to extreme levels of ionizing radiation (IR). Although multiple unique regulatory systems (e.g., the Radiation and Desiccation Response (RDR)) have been identified in this organism, the role of post-transcriptional regulators has not been characterized within the IR response. In this study, we have characterized an sRNA, PprS (formerly Dsr2), as a post-transcriptional coordinator of IR recovery in D. radiodurans. PprS showed differential expression specifically under IR and knockdown of PprS resulted in reduced survival and growth under IR, suggesting its importance in regulating post-radiation recovery. We determined a number of potential RNA targets involved in several pathways including translation and DNA repair. Specifically, we confirmed that PprS binds within the coding region to stabilize the pprM (DR_0907) transcript, a RDR modulator. Overall, these results are the first to present an additional layer of sRNA-based control in DNA repair pathways associated with bacterial radioresistance.

Characterization of the DdrD protein from the extremely radioresistant bacterium Deinococcus radiodurans

Here, we report the in vitro and in vivo characterization of the DdrD protein from the extraordinary stress-resistant bacterium, D. radiodurans. DdrD is one of the most highly induced proteins following cellular irradiation or desiccation. We confirm that DdrD belongs to the Radiation Desiccation Response (RDR) regulon protein family whose expression is regulated by the IrrE/DdrO proteins after DNA damage. We show that DdrD is a DNA binding protein that binds to single-stranded DNA In vitro, but not to duplex DNA unless it has a 5' single-stranded extension. In vivo, we observed no significant effect of the absence of DdrD on the survival of D. radiodurans cells after exposure to γ-rays or UV irradiation in different genetic contexts. However, genome reassembly is affected in a ∆ddrD mutant when cells recover from irradiation in the absence of nutrients. Thus, DdrD likely contributes to genome reconstitution after irradiation, but only under starvation conditions. Lastly, we show that the absence of the DdrD protein partially restores the frequency of plasmid transformation of a ∆ddrB mutant, suggesting that DdrD could also be involved in biological processes other than the response to DNA damage.

Coexistence of SOS-Dependent and SOS-Independent Regulation of DNA Repair Genes in Radiation-Resistant Deinococcus Bacteria

Deinococcus bacteria are extremely resistant to radiation and able to repair a shattered genome in an essentially error-free manner after exposure to high doses of radiation or prolonged desiccation. An efficient, SOS-independent response mechanism to induce various DNA repair genes such as recA is essential for radiation resistance. This pathway, called radiation/desiccation response, is controlled by metallopeptidase IrrE and repressor DdrO that are highly conserved in Deinococcus. Among various Deinococcus species, Deinococcus radiodurans has been studied most extensively. Its genome encodes classical DNA repair proteins for error-free repair but no error-prone translesion DNA polymerases, which may suggest that absence of mutagenic lesion bypass is crucial for error-free repair of massive DNA damage. However, many other radiation-resistant Deinococcus species do possess translesion polymerases, and radiation-induced mutagenesis has been demonstrated. At least dozens of Deinococcus species contain a mutagenesis cassette, and some even two cassettes, encoding error-prone translesion polymerase DnaE2 and two other proteins, ImuY and ImuB-C, that are probable accessory factors required for DnaE2 activity. Expression of this mutagenesis cassette is under control of the SOS regulators RecA and LexA. In this paper, we review both the RecA/LexA-controlled mutagenesis and the IrrE/DdrO-controlled radiation/desiccation response in Deinococcus.

Redox signaling through zinc activates the radiation response in Deinococcus bacteria

Deinococcus bacteria are extremely resistant to radiation and other DNA damage- and oxidative stress-generating conditions. An efficient SOS-independent response mechanism inducing expression of several DNA repair genes is essential for this resistance, and is controlled by metalloprotease IrrE that cleaves and inactivates transcriptional repressor DdrO. Here, we identify the molecular signaling mechanism that triggers DdrO cleavage. We show that reactive oxygen species (ROS) stimulate the zinc-dependent metalloprotease activity of IrrE in Deinococcus. Sudden exposure of Deinococcus to zinc excess also rapidly induces DdrO cleavage, but is not accompanied by ROS production and DNA damage. Further, oxidative treatment leads to an increase of intracellular free zinc, indicating that IrrE activity is very likely stimulated directly by elevated levels of available zinc ions. We conclude that radiation and oxidative stress induce changes in redox homeostasis that result in IrrE activation by zinc in Deinococcus. We propose that a part of the zinc pool coordinated with cysteine thiolates is released due to their oxidation. Predicted regulation systems involving IrrE- and DdrO-like proteins are present in many bacteria, including pathogens, suggesting that such a redox signaling pathway including zinc as a second messenger is widespread and participates in various stress responses.

Lost and Found: Re-searching and Re-scoring Proteomics Data Aids Genome Annotation and Improves Proteome Coverage

Delineation of open reading frames (ORFs) causes persistent inconsistencies in prokaryote genome annotation. We demonstrate that by advanced (re)analysis of omics data, a higher proteome coverage and sensitive detection of unannotated ORFs can be achieved, which can be exploited for conditional bacterial genome (re)annotation, which is especially relevant in view of annotating the wealth of sequenced prokaryotic genomes obtained in recent years.

Prokaryotic genome annotation is heavily dependent on automated gene annotation pipelines that are prone to propagate errors and underestimate genome complexity. We describe an optimized proteogenomic workflow that uses ribosome profiling (ribo-seq) and proteomic data for Salmonella enterica serovar Typhimurium to identify unannotated proteins or alternative protein forms. This data analysis encompasses the searching of cofragmenting peptides and postprocessing with extended peptide-to-spectrum quality features, including comparison to predicted fragment ion intensities. When this strategy is applied, an enhanced proteome depth is achieved, as well as greater confidence for unannotated peptide hits. We demonstrate the general applicability of our pipeline by reanalyzing public Deinococcus radiodurans data sets. Taken together, our results show that systematic reanalysis using available prokaryotic (proteome) data sets holds great promise to assist in experimentally based genome annotation.

IMPORTANCE Delineation of open reading frames (ORFs) causes persistent inconsistencies in prokaryote genome annotation. We demonstrate that by advanced (re)analysis of omics data, a higher proteome coverage and sensitive detection of unannotated ORFs can be achieved, which can be exploited for conditional bacterial genome (re)annotation, which is especially relevant in view of annotating the wealth of sequenced prokaryotic genomes obtained in recent years.

System-level understanding of gene expression and regulation for engineering secondary metabolite production in Streptomyces

The gram-positive bacterium, Streptomyces, is noticed for its ability to produce a wide array of pharmaceutically active compounds through secondary metabolism. To discover novel bioactive secondary metabolites and increase the production, Streptomyces species have been extensively studied for the past decades. Among the cellular components, RNA molecules play important roles as the messengers for gene expression and diverse regulations taking place at the RNA level. Thus, the analysis of RNA-level regulation is critical to understanding the regulation of Streptomyces’ metabolism and secondary metabolite production. A dramatic advance in Streptomyces research was made recently, by exploiting high-throughput technology to systematically understand RNA levels. In this review, we describe the current status of the system-wide investigation of Streptomyces in terms of RNA, toward expansion of its genetic potential for secondary metabolite synthesis.

Recent Advances in Archaeal Translation Initiation

Translation initiation (TI) allows accurate selection of the initiation codon on a messenger RNA (mRNA) and defines the reading frame. In all domains of life, translation initiation generally occurs within a macromolecular complex made up of the small ribosomal subunit, the mRNA, a specialized methionylated initiator tRNA, and translation initiation factors (IFs). Once the start codon is selected at the P site of the ribosome and the large subunit is associated, the IFs are released and a ribosome competent for elongation is formed. However, even if the general principles are the same in the three domains of life, the molecular mechanisms are different in bacteria, eukaryotes, and archaea and may also vary depending on the mRNA. Because TI mechanisms have evolved lately, their studies bring important information about the evolutionary relationships between extant organisms. In this context, recent structural data on ribosomal complexes and genome-wide studies are particularly valuable. This review focuses on archaeal translation initiation highlighting its relationships with either the eukaryotic or the bacterial world. Eukaryotic features of the archaeal small ribosomal subunit are presented. Ribosome evolution and TI mechanisms diversity in archaeal branches are discussed. Next, the use of leaderless mRNAs and that of leadered mRNAs having Shine-Dalgarno sequences is analyzed. Finally, the current knowledge on TI mechanisms of SD-leadered and leaderless mRNAs is detailed.

Crystal structure of the transcriptional repressor DdrO: insight into the metalloprotease/repressor-controlled radiation response in Deinococcus

Exposure to harmful conditions such as radiation and desiccation induce oxidative stress and DNA damage. In radiation-resistant Deinococcus bacteria, the radiation/desiccation response is controlled by two proteins: the XRE family transcriptional repressor DdrO and the COG2856 metalloprotease IrrE. The latter cleaves and inactivates DdrO. Here, we report the biochemical characterization and crystal structure of DdrO, which is the first structure of a XRE protein targeted by a COG2856 protein. DdrO is composed of two domains that fold independently and are separated by a flexible linker. The N-terminal domain corresponds to the DNA-binding domain. The C-terminal domain, containing three alpha helices arranged in a novel fold, is required for DdrO dimerization. Cleavage by IrrE occurs in the loop between the last two helices of DdrO and abolishes dimerization and DNA binding. The cleavage site is hidden in the DdrO dimer structure, indicating that IrrE cleaves DdrO monomers or that the interaction with IrrE induces a structural change rendering accessible the cleavage site. Predicted COG2856/XRE regulatory protein pairs are found in many bacteria, and available data suggest two different molecular mechanisms for stress-induced gene expression: COG2856 protein-mediated cleavage or inhibition of oligomerization without cleavage of the XRE repressor.

Full-Length Transcript-Based Proteogenomics of Rice Improves Its Genome and Proteome Annotation

Rice (Oryza sativa) molecular breeding has gained considerable attention in recent years, but inaccurate genome annotation hampers its progress and functional studies of the rice genome. In this study, we applied single-molecule long-read RNA sequencing (lrRNA_seq)-based proteogenomics to reveal the complexity of the rice transcriptome and its coding abilities. Surprisingly, approximately 60% of loci identified by lrRNA_seq are associated with natural antisense transcripts (NATs). The high-density genomic arrangement of NAT genes suggests their potential roles in the multifaceted control of gene expression. In addition, a large number of fusion and intergenic transcripts have been observed. Furthermore, 906,456 transcript isoforms were identified, and 72.9% of the genes can generate splicing isoforms. A total of 706,075 posttranscriptional events were subsequently categorized into 10 subtypes, demonstrating the interdependence of posttranscriptional mechanisms that contribute to transcriptome diversity. Parallel short-read RNA sequencing indicated that lrRNA_seq has a superior capacity for the identification of longer transcripts. In addition, over 190,000 unique peptides belonging to 9,706 proteoforms/protein groups were identified, expanding the diversity of the rice proteome. Our findings indicate that the genome organization, transcriptome diversity, and coding potential of the rice transcriptome are far more complex than previously anticipated.

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

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

Conservation and diversity of radiation and oxidative stress resistance mechanisms in Deinococcus species

Deinococcus bacteria are famous for their extreme resistance to ionising radiation and other DNA damage- and oxidative stress-generating agents. More than a hundred genes have been reported to contribute to resistance to radiation, desiccation and/or oxidative stress in Deinococcus radiodurans. These encode proteins involved in DNA repair, oxidative stress defence, regulation and proteins of yet unknown function or with an extracytoplasmic location. Here, we analysed the conservation of radiation resistance-associated proteins in other radiation-resistant Deinococcus species. Strikingly, homologues of dozens of these proteins are absent in one or more Deinococcus species. For example, only a few Deinococcus-specific proteins and radiation resistance-associated regulatory proteins are present in each Deinococcus, notably the metallopeptidase/repressor pair IrrE/DdrO that controls the radiation/desiccation response regulon. Inversely, some Deinococcus species possess proteins that D. radiodurans lacks, including DNA repair proteins consisting of novel domain combinations, translesion polymerases, additional metalloregulators, redox-sensitive regulator SoxR and manganese-containing catalase. Moreover, the comparisons improved the characterisation of several proteins regarding important conserved residues, cellular location and possible protein–protein interactions. This comprehensive analysis indicates not only conservation but also large diversity in the molecular mechanisms involved in radiation resistance even within the Deinococcus genus.

Leaderless mRNAs in the Spotlight: Ancient but Not Outdated!

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

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.

RNAseq analysis of α-proteobacterium Gluconobacter oxydans 621H

Background

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

Results

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

Conclusions

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

Electronic supplementary material

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

A Genome-Wide Search for Ionizing-Radiation-Responsive Elements in Deinococcus radiodurans Reveals a Regulatory Role for the DNA Gyrase Subunit A Gene's 5' Untranslated Region in the Radiation and Desiccation Response

Tight regulation of gene expression is important for the survival of Deinococcus radiodurans, a model bacterium of extreme stress resistance. Few studies have examined the use of regulatory RNAs as a possible contributing mechanism to ionizing radiation (IR) resistance, despite their proffered efficient and dynamic gene expression regulation under IR stress. This work presents a transcriptome-based approach for the identification of stress-responsive regulatory 5' untranslated region (5'-UTR) elements in D. radiodurans R1 that can be broadly applied to other bacteria. Using this platform and an in vivo fluorescence screen, we uncovered the presence of a radiation-responsive regulatory motif in the 5' UTR of the DNA gyrase subunit A gene. Additional screens under H2O2-induced oxidative stress revealed the specificity of the response of this element to IR stress. Further examination of the sequence revealed a regulatory motif of the radiation and desiccation response (RDR) in the 5' UTR that is necessary for the recovery of D. radiodurans from high doses of IR. Furthermore, we suggest that it is the preservation of predicted RNA structure, in addition to DNA sequence consensus of the motif, that permits this important regulatory ability.IMPORTANCEDeinococcus radiodurans is an extremely stress-resistant bacterium capable of tolerating up to 3,000 times more ionizing radiation than human cells. As an integral part of the stress response mechanism of this organism, we suspect that it maintains stringent control of gene expression. However, understanding of its regulatory pathways remains incomplete to date. Untranslated RNA elements have been demonstrated to play crucial roles in gene regulation throughout bacteria. In this work, we focus on searching for and characterizing responsive RNA elements under radiation stress and propose that multiple levels of gene regulation work simultaneously to enable this organism to efficiently recover from exposure to ionizing radiation. The model we propose serves as a generic template to investigate similar mechanisms of gene regulation under stress that have likely evolved in other bacterial species.

In vivo and in vitro characterization of DdrC, a DNA damage response protein in Deinococcus radiodurans bacterium

The bacterium Deinococcus radiodurans possesses a set of Deinococcus-specific genes highly induced after DNA damage. Among them, ddrC (dr0003) was recently re-annotated, found to be in the inverse orientation and called A2G07_00380. Here, we report the first in vivo and in vitro characterization of the corrected DdrC protein to better understand its function in irradiated cells. In vivo, the ΔddrC null mutant is sensitive to high doses of UV radiation and the ddrC deletion significantly increases UV-sensitivity of ΔuvrA or ΔuvsE mutant strains. We show that the expression of the DdrC protein is induced after γ-irradiation and is under the control of the regulators, DdrO and IrrE. DdrC is rapidly recruited into the nucleoid of the irradiated cells. In vitro, we show that DdrC is able to bind single- and double-stranded DNA with a preference for the single-stranded DNA but without sequence or shape specificity and protects DNA from various nuclease attacks. DdrC also condenses DNA and promotes circularization of linear DNA. Finally, we show that the purified protein exhibits a DNA strand annealing activity. Altogether, our results suggest that DdrC is a new DNA binding protein with pleiotropic activities. It might maintain the damaged DNA fragments end to end, thus limiting their dispersion and extensive degradation after exposure to ionizing radiation. DdrC might also be an accessory protein that participates in a single strand annealing pathway whose importance in DNA repair becomes apparent when DNA is heavily damaged.

Conservation and diversity of the IrrE/DdrO-controlled radiation response in radiation-resistant Deinococcus bacteria

The extreme radiation resistance of Deinococcus bacteria requires the radiation‐stimulated cleavage of protein DdrO by a specific metalloprotease called IrrE. DdrO is the repressor of a predicted radiation/desiccation response (RDR) regulon, composed of radiation‐induced genes having a conserved DNA motif (RDRM) in their promoter regions. Here, we showed that addition of zinc ions to purified apo‐IrrE, and short exposure of Deinococcus cells to zinc ions, resulted in cleavage of DdrO in vitro and in vivo, respectively. Binding of IrrE to RDRM‐containing DNA or interaction of IrrE with DNA‐bound DdrO was not observed. The data are in line with IrrE being a zinc peptidase, and indicate that increased zinc availability, caused by oxidative stress, triggers the in vivo cleavage of DdrO unbound to DNA. Transcriptomics and proteomics of Deinococcus deserti confirmed the IrrE‐dependent regulation of predicted RDR regulon genes and also revealed additional members of this regulon. Comparative analysis showed that the RDR regulon is largely well conserved in Deinococcus species, but also showed diversity in the regulon composition. Notably, several RDR genes with an important role in radiation resistance in Deinococcus radiodurans, for example pprA, are not conserved in some other radiation‐resistant Deinococcus species.

In situ real-time evaluation of radiation-responsive promoters in the extremely radioresistant microbe Deinococcus radiodurans

A third generation promoter probe shuttle vector pKG was constructed, using the green fluorescent protein as a reporter, for in situ evaluation of Deinococcal promoter activity in Escherichia coli or Deinococcus radiodurans. The construct yielded zero background fluorescence in both the organisms, in the absence of promoter sequences. Fifteen Deinococcal promoters, either harbouring Radiation and Desiccation Response Motif (RDRM) or not, were cloned in vector pKG. Only the RDRM-promoter constructs displayed (i) gamma radiation inducible GFP expression in D. radiodurans, following gamma irradiation, (ii) DdrO-mediated repression of GFP expression in heterologous E. coli, or (iii) abolition in GFP induction following gamma irradiation, in pprI mutant of D. radiodurans. Utility of pKG vector for real-time in situ assessment of Deinococcal promoter function was, thus, successfully demonstrated.

Genome-Wide Transcriptional Start Site Mapping and sRNA Identification in the Pathogen Leptospira interrogans

Leptospira are emerging zoonotic pathogens transmitted from animals to humans typically through contaminated environmental sources of water and soil. Regulatory pathways of pathogenic Leptospira spp. underlying the adaptive response to different hosts and environmental conditions remains elusive. In this study, we provide the first global Transcriptional Start Site (TSS) map of a Leptospira species. RNA was obtained from the pathogen Leptospira interrogans grown at 30°C (optimal in vitro temperature) and 37°C (host temperature) and selectively enriched for 5′ ends of native transcripts. A total of 2865 and 2866 primary TSS (pTSS) were predicted in the genome of L. interrogans at 30 and 37°C, respectively. The majority of the pTSSs were located between 0 and 10 nucleotides from the translational start site, suggesting that leaderless transcripts are a common feature of the leptospiral translational landscape. Comparative differential RNA-sequencing (dRNA-seq) analysis revealed conservation of most pTSS at 30 and 37°C. Promoter prediction algorithms allow the identification of the binding sites of the alternative sigma factor sigma 54. However, other motifs were not identified indicating that Leptospira consensus promoter sequences are inherently different from the Escherichia coli model. RNA sequencing also identified 277 and 226 putative small regulatory RNAs (sRNAs) at 30 and 37°C, respectively, including eight validated sRNAs by Northern blots. These results provide the first global view of TSS and the repertoire of sRNAs in L. interrogans. These data will establish a foundation for future experimental work on gene regulation under various environmental conditions including those in the host.

Effects of Kasugamycin on the Translatome of Escherichia coli

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

Bacterial host and reporter gene optimization for genetically encoded whole cell biosensors

Whole-cell biosensors based on reporter genes allow detection of toxic metals in water with high selectivity and sensitivity under laboratory conditions; nevertheless, their transfer to a commercial inline water analyzer requires specific adaptation and optimization to field conditions as well as economical considerations. We focused here on both the influence of the bacterial host and the choice of the reporter gene by following the responses of global toxicity biosensors based on constitutive bacterial promoters as well as arsenite biosensors based on the arsenite-inducible P_ars promoter. We observed important variations of the bioluminescence emission levels in five different Escherichia coli strains harboring two different lux-based biosensors, suggesting that the best host strain has to be empirically selected for each new biosensor under construction. We also investigated the bioluminescence reporter gene system transferred into Deinococcus deserti, an environmental, desiccation- and radiation-tolerant bacterium that would reduce the manufacturing costs of bacterial biosensors for commercial water analyzers and open the field of biodetection in radioactive environments. We thus successfully obtained a cell survival biosensor and a metal biosensor able to detect a concentration as low as 100 nM of arsenite in D. deserti. We demonstrated that the arsenite biosensor resisted desiccation and remained functional after 7 days stored in air-dried D. deserti cells. We also report here the use of a new near-infrared (NIR) fluorescent reporter candidate, a bacteriophytochrome from the magnetotactic bacterium Magnetospirillum magneticum AMB-1, which showed a NIR fluorescent signal that remained optimal despite increasing sample turbidity, while in similar conditions, a drastic loss of the lux-based biosensors signal was observed.

Using Galaxy for Proteomics

The area of informatics for mass spectrometry (MS)-based proteomics data has steadily grown over the last two decades. Numerous, effective software programs now exist for various aspects of proteomic informatics. However, many researchers still have difficulties in using these software. These difficulties arise from problems with running and integrating disparate software programs, scalability issues when dealing with large data volumes, and lack of ability to share and reproduce workflows comprised of different software. The Galaxy framework for bioinformatics provides an attractive option for solving many of these current issues in proteomic informatics. Originally developed as a workbench to enable genomic data analysis, numerous researchers are now turning to Galaxy to implement software for MS-based proteomics applications. Here, we provide an introduction to Galaxy and its features, and describe how software tools are deployed, published and shared via the scalable framework. We also describe some of the existing tools in Galaxy for basic MS-based proteomics data analysis and informatics. Finally, we describe how proteomics tools in Galaxy can be combined with other existing tools for genomic and transcriptomic data analysis to enable powerful multi-omic data analysis applications.

Reannotation of Genomes by Means of Proteomics Data

Omics approaches have become popular in biology as powerful discovery tools, and currently gain in interest for diagnostic applications. Establishing the accurate genome sequence of any organism is easy, but the outcome of its annotation by means of automatic pipelines remains imprecise. Some protein-encoding genes may be missed as soon as they are specific and poorly conserved in a given taxon, while important to explain the specific traits of the organism. Translational starts are also poorly predicted in a relatively important number of cases, thus impacting the protein sequence database used in proteomics, comparative genomics, and systems biology. The use of high-throughput proteomics data to improve genome annotation is an attractive option to obtain a more comprehensive molecular picture of a given organism. Here, protocols for reannotating prokaryote genomes are described based on shotgun proteomics and derivatization of protein N-termini with a positively charged reagent coupled to high-resolution tandem mass spectrometry.

Proteogenomics for the Comprehensive Analysis of Human Cellular and Serum Antibody Repertoires

The vast repertoire of immunoglobulins produced by the immune system is a consequence of the huge amount of antigens to which we are exposed every day. The diversity of these immunoglobulins is due to different mechanisms (including VDJ recombination, somatic hypermutation, and antigen selection). Understanding how the immune system is capable of generating this diversity and which are the molecular bases of the composition of immunoglobulins are key challenges in the immunological field. During the last decades, several techniques have emerged as promising strategies to achieve these goals, but it is their combination which appears to be the fruitful solution for increasing the knowledge about human cellular and serum antibody repertoires.In this chapter, we address the diverse strategies focused on the analysis of immunoglobulin repertoires as well as the characterization of the genomic and peptide sequences. Moreover, the advantages of combining various -omics approaches are discussed through review different published studies, showing the benefits in clinical areas.

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

BACKGROUND

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

RESULTS

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

CONCLUSION

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

A novel enrichment strategy reveals unprecedented number of novel transcription start sites at single base resolution in a model prokaryote and the gut microbiome

Background

The initiating nucleotide found at the 5’ end of primary transcripts has a distinctive triphosphorylated end that distinguishes these transcripts from all other RNA species. Recognizing this distinction is key to deconvoluting the primary transcriptome from the plethora of processed transcripts that confound analysis of the transcriptome. The currently available methods do not use targeted enrichment for the 5′end of primary transcripts, but rather attempt to deplete non-targeted RNA.

Results

We developed a method, Cappable-seq, for directly enriching for the 5' end of primary transcripts and enabling determination of transcription start sites at single base resolution. This is achieved by enzymatically modifying the 5′ triphosphorylated end of RNA with a selectable tag. We first applied Cappable-seq to E. coli, achieving up to 50 fold enrichment of primary transcripts and identifying an unprecedented 16539 transcription start sites (TSS) genome-wide at single base resolution. We also applied Cappable-seq to a mouse cecum sample and identified TSS in a microbiome.

Conclusions

Cappable-seq allows for the first time the capture of the 5′ end of primary transcripts. This enables a unique robust TSS determination in bacteria and microbiomes.  In addition to and beyond TSS determination, Cappable-seq depletes ribosomal RNA and reduces the complexity of the transcriptome to a single quantifiable tag per transcript enabling digital profiling of gene expression in any microbiome.

Electronic supplementary material

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

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

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

Proteogenomics of rare taxonomic phyla: A prospective treasure trove of protein coding genes

Sustainable innovations in sequencing technologies have resulted in a torrent of microbial genome sequencing projects. However, the prokaryotic genomes sequenced so far are unequally distributed along their phylogenetic tree; few phyla contain the majority, the rest only a few representatives. Accurate genome annotation lags far behind genome sequencing. While automated computational prediction, aided by comparative genomics, remains a popular choice for genome annotation, substantial fraction of these annotations are erroneous. Proteogenomics utilizes protein level experimental observations to annotate protein coding genes on a genome wide scale. Benefits of proteogenomics include discovery and correction of gene annotations regardless of their phylogenetic conservation. This not only allows detection of common, conserved proteins but also the discovery of protein products of rare genes that may be horizontally transferred or taxonomy specific. Chances of encountering such genes are more in rare phyla that comprise a small number of complete genome sequences. We collated all bacterial and archaeal proteogenomic studies carried out to date and reviewed them in the context of genome sequencing projects. Here, we present a comprehensive list of microbial proteogenomic studies, their taxonomic distribution, and also urge for targeted proteogenomics of underexplored taxa to build an extensive reference of protein coding genes.

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

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

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

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

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

Differential RNA-seq (dRNA-seq) for annotation of transcriptional start sites and small RNAs in Helicobacter pylori

The global mapping of transcription boundaries is a key step in the elucidation of the full complement of transcriptional features of an organism. It facilitates the annotation of operons and untranslated regions as well as novel transcripts, including cis- and trans-encoded small RNAs (sRNAs). So called RNA sequencing (RNA-seq) based on deep sequencing of cDNAs has greatly facilitated transcript mapping with single nucleotide resolution. However, conventional RNA-seq approaches typically cannot distinguish between primary and processed transcripts. Here we describe the recently developed differential RNA-seq (dRNA-seq) approach, which facilitates the annotation of transcriptional start sites (TSS) based on deep sequencing of two differentially treated cDNA library pairs, with one library being enriched for primary transcripts. Using the human pathogen Helicobacter pylori as a model organism, we describe the application of dRNA-seq together with an automated TSS annotation approach for generation of a genome-wide TSS map in bacteria. Besides a description of transcriptome and regulatory features that can be identified by this approach, we discuss the impact of different library preparation protocols and sequencing platforms as well as manual and automated TSS annotation. Moreover, we have set up an easily accessible online browser for visualization of the H. pylori transcriptome data from this and our previous H. pylori dRNA-seq study.

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

Background

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

Results

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

Conclusions

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

Electronic supplementary material

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

Expression and mutational analysis of DinB-like protein DR0053 in Deinococcus radiodurans

In order to understand the mechanism governing radiation resistance in Deinococcus radiodurans, current efforts are aimed at identifying potential candidates from a large repertoire of unique Deinococcal genes and protein families. DR0053 belongs to the DinB/YfiT protein family, which is an over-represented protein family in D. radiodurans. We observed that dr0053 transcript levels were highly induced in response to gamma radiation (γ-radiation) and mitomycin C (MMC) exposure depending on PprI, RecA and the DrtR/S two-component signal transduction system. Protein profiles demonstrated that DR0053 is a highly induced protein in cultures exposed to 10 kGy γ-radiation. We were able to determine the transcriptional start site of dr0053, which was induced upon irradiation, and to assign the 133-bp promoter region of dr0053 as essential for radiation responsiveness through primer extension and promoter deletion analyses. A dr0053 mutant strain displayed sensitivity to γ-radiation and MMC exposure, but not hydrogen peroxide, suggesting that DR0053 helps cells recover from DNA damage. Bioinformatic analyses revealed that DR0053 is similar to the Bacillus subtilis protein YjoA, which is a substrate of bacterial protein-tyrosine kinases. Taken together, the DNA damage-inducible (din) gene dr0053 may be regulated at the transcriptional and post-translational levels.

Proteogenomics in microbiology: taking the right turn at the junction of genomics and proteomics

High-accuracy and high-throughput proteomic methods have completely changed the way we can identify and characterize proteins. MS-based proteomics can now provide a unique supplement to genomic data and add a new level of information to the interpretation of genomic sequences. Proteomics-driven genome annotation has become especially relevant in microbiology where genomes are sequenced on a daily basis and limitations of an in silico driven annotation process are well recognized. In this review paper, we outline different strategies on how one can design a proteogenomic experiment, for example on genome-sequenced (synonymous proteogenomics) versus unsequenced organisms (ortho-proteogenomics) or with the aid of other "omic" data such as RNA-seq. We touch upon many challenges that are encountered during a typical proteogenomic study, mostly concerning bioinformatics methods and downstream data analysis, but also related to creation and use of sequence databases. A large list of proteogenomic case studies of different microorganisms is provided to illustrate the mapping of MS/MS-derived peptide spectra to genomic DNA sequences. These investigations have led to accurate determination of translational initiation sites, pointed out eventual read-throughs or programmed frameshifts, detected signal peptide processing or other protein maturation events, removed questionable annotation assignments, and provided evidence for predicted hypothetical proteins.

N-terminomics and proteogenomics, getting off to a good start

Proteogenomics consists of the annotation or reannotation of protein-coding nucleic acid sequences based on the empirical observation of their gene products. While functional annotation of predicted genes is increasingly feasible given the multiplicity of genomes available for many branches of the tree of life, the accurate annotation of the translational start sites is still a point of contention. Extensive coverage of the proteome, including specifically the N-termini, is now possible, thanks to next-generation mass spectrometers able to record data from thousands of proteins at once. Efforts to increase the peptide coverage of protein sequences and to detect low abundance proteins are important to make proteomic and proteogenomic studies more comprehensive. In this review, we present the panoply of N-terminus-oriented strategies that have been developed over the last decade.

Radiation response in Deinococcus deserti: IrrE is a metalloprotease that cleaves repressor protein DdrO

Deinococcus bacteria are famous for their extreme radiation tolerance. The IrrE protein was shown to be essential for radiation tolerance and, in an unelucidated manner, for induction of a number of genes in response to radiation, including recA and other DNA repair genes. Earlier studies indicated that IrrE could be a zinc peptidase, but proteolytic activity was not demonstrated. Here, using several in vivo and in vitro experiments, IrrE from Deinococcus deserti was found to interact with DdrO, a predicted regulator encoded by a radiation-induced gene that is, like irrE, highly conserved in Deinococcus. Moreover, IrrE was found to cleave DdrO in vitro and when the proteins were coexpressed in Escherichia coli. This cleavage was not observed in the presence of metal chelator EDTA or when IrrE contains a mutation in the conserved active-site motif of metallopeptidases. In D. deserti, IrrE-dependent cleavage of DdrO was observed after exposure to radiation. Furthermore, DdrO-dependent repression of the promoter of a radiation-induced gene was shown. These results demonstrate that IrrE is a metalloprotease and we propose that IrrE-mediated cleavage inactivates repressor protein DdrO, leading to transcriptional induction of various genes required for repair and survival after exposure of Deinococcus to radiation.