An expanding universe of small proteins

A catalog of small proteins from the global microbiome

Small open reading frames (smORFs) shorter than 100 codons are widespread and perform essential roles in microorganisms, where they encode proteins active in several cell functions, including signal pathways, stress response, and antibacterial activities. However, the ecology, distribution and role of small proteins in the global microbiome remain unknown. Here, we construct a global microbial smORFs catalog (GMSC) derived from 63,410 publicly available metagenomes across 75 distinct habitats and 87,920 high-quality isolate genomes. GMSC contains 965 million non-redundant smORFs with comprehensive annotations. We find that archaea harbor more smORFs proportionally than bacteria. We moreover provide a tool called GMSC-mapper to identify and annotate small proteins from microbial (meta)genomes. Overall, this publicly-available resource demonstrates the immense and underexplored diversity of small proteins.

Type I toxin-antitoxin systems in bacteria: from regulation to biological functions

Toxin-antitoxin systems are ubiquitous in the prokaryotic world and widely distributed among chromosomes and mobile genetic elements. Several different toxin-antitoxin system types exist, but what they all have in common is that toxin activity is prevented by the cognate antitoxin. In type I toxin-antitoxin systems, toxin production is controlled by an RNA antitoxin and by structural features inherent to the toxin messenger RNA. Most type I toxins are small membrane proteins that display a variety of cellular effects. While originally discovered as modules that stabilize plasmids, chromosomal type I toxin-antitoxin systems may also stabilize prophages, or serve important functions upon certain stress conditions and contribute to population-wide survival strategies. Here, we will describe the intricate RNA-based regulation of type I toxin-antitoxin systems and discuss their potential biological functions.

Small proteins in bacteria - Big challenges in prediction and identification

Proteins with up to 100 amino acids have been largely overlooked due to the challenges associated with predicting and identifying them using traditional methods. Recent advances in bioinformatics and machine learning, DNA sequencing, RNA and Ribo-seq technologies, and mass spectrometry (MS) have greatly facilitated the detection and characterisation of these elusive proteins in recent years. This has revealed their crucial role in various cellular processes including regulation, signalling and transport, as toxins and as folding helpers for protein complexes. Consequently, the systematic identification and characterisation of these proteins in bacteria have emerged as a prominent field of interest within the microbial research community. This review provides an overview of different strategies for predicting and identifying these proteins on a large scale, leveraging the power of these advanced technologies. Furthermore, the review offers insights into the future developments that may be expected in this field.

The sbiTRS Operon Contributes to Stenobactin-Mediated Iron Utilization in Stenotrophomonas maltophilia

Iron is an essential micronutrient for various bacterial cellular processes. Fur is a global transcriptional regulator participating in iron homeostasis. Stenotrophomonas maltophilia is a ubiquitous environmental bacterium that has emerged as an opportunistic pathogen. To elucidate the novel regulatory mechanism behind iron homeostasis in S. maltophilia, wild-type KJ and KJΔFur, a fur mutant, were subjected to transcriptome assay. A five-gene cluster, sbiBA-sbiTRS, was significantly upregulated in KJΔFur. SbiAB is an ATP type efflux pump, SbiT is an inner membrane protein, and SbiSR is a two-component regulatory system (TCS). The sbiTRS operon organization was verified by reverse transcription-PCR (RT-PCR). Localization prediction and bacterial two-hybrid studies revealed that SbiT resided in the inner membrane and had an intramembrane interaction with SbiS. In iron-replete conditions, SbiT interacted with SbiS and maintained SbiSR TCS in a resting state. In response to iron depletion stress, SbiT no longer interacted with SbiS, leading to SbiSR TCS activation. The iron source utilization assay demonstrated the contribution of SbiSR TCS to stenobactin-mediated ferric iron utilization but notto the utilization of hemin and ferric citrate. Furthermore, SmeDEF and SbiAB pumps, known stenobactin secretion outlets, were members of the SbiSR regulon. Collectively, in an iron-depleted condition, SbiSR activation is regulated by Fur at the transcriptional level and by SbiT at the posttranslational level. Activated SbiSR contributes to stenobactin-mediated ferric iron utilization by upregulating the smeDEF and sbiAB operons. SbiSR is the first TCS found to be involved in iron homeostasis in S. maltophilia. IMPORTANCE Therapeutic options for Stenotrophomonas maltophilia infections are limited because S. maltophilia is intrinsically resistant to several antibiotics. Iron is an essential element for viability, but iron overload is a lethal threat to bacteria. Therefore, disruption of iron homeostasis can be an alternative strategy to cope with S. maltophilia infection. The intricate regulatory networks involved in iron hemostasis have been reported in various pathogens; however, little is known about S. maltophilia. Herein, a novel sbiTRS operon, a member of Fur regulon, was characterized. SbiT, an inner membrane protein, negatively modulated the SbiSR two-component regulatory system by intramembrane protein-protein interaction with SbiS. In response to iron-depleted stress, SbiSR was activated via the regulation of Fur and SbiT. Activated SbiSR upregulated smeDEF and sbiAB, which contributed to stenobactin-mediated ferric iron utilization. A novel fur-sbiT-sbiSR-smeDEF/sbiAB regulatory circuit in S. maltophilia was revealed.

Mapping functional regions of essential bacterial proteins with dominant-negative protein fragments

Massively parallel measurements of dominant-negative inhibition by protein fragments have been used to map protein interaction sites and discover peptide inhibitors. However, the underlying principles governing fragment-based inhibition have thus far remained unclear. Here, we adapted a high-throughput inhibitory fragment assay for use in Escherichia coli, applying it to a set of 10 essential proteins. This approach yielded single amino acid resolution maps of inhibitory activity, with peaks localized to functionally important interaction sites, including oligomerization interfaces and folding contacts. Leveraging these data, we performed a systematic analysis to uncover principles of fragment-based inhibition. We determined a robust negative correlation between susceptibility to inhibition and cellular protein concentration, demonstrating that inhibitory fragments likely act primarily by titrating native protein interactions. We also characterized a series of trade-offs related to fragment length, showing that shorter peptides allow higher-resolution mapping but suffer from lower inhibitory activity. We employed an unsupervised statistical analysis to show that the inhibitory activities of protein fragments are largely driven not by generic properties such as charge, hydrophobicity, and secondary structure, but by the more specific characteristics of their bespoke macromolecular interactions. Overall, this work demonstrates fundamental characteristics of inhibitory protein fragment function and provides a foundation for understanding and controlling protein interactions in vivo.

Exploring the Peptide Potential of Genomes

Recent studies attribute a central role to the noncoding genome in the emergence of novel genes. The widespread transcription of noncoding regions and the pervasive translation of the resulting RNAs offer to the organisms a vast reservoir of novel peptides. Although the majority of these peptides are anticipated as deleterious or neutral, and thereby expected to be degraded right away or short-lived in evolutionary history, some of them can confer an advantage to the organism. The latter can be further subjected to natural selection and be established as novel genes. In any case, characterizing the structural properties of these pervasively translated peptides is crucial to understand (1) their impact on the cell and (2) how some of these peptides, derived from presumed noncoding regions, can give rise to structured and functional de novo proteins. Therefore, we present a protocol that aims to explore the potential of a genome to produce novel peptides. It consists in annotating all the open reading frames (ORFs) of a genome (i.e., coding and noncoding ones) and characterizing the fold potential and other structural properties of their corresponding potential peptides. Here, we apply our protocol to a small genome and show how to apply it to very large genomes. Finally, we present a case study which aims to probe the fold potential of a set of 721 translated ORFs in mouse lncRNAs, identified with ribosome profiling experiments. Interestingly, we show that the distribution of their fold potential is different from that of the nontranslated lncRNAs and more generally from the other noncoding ORFs of the mouse.

CdrS Is a Global Transcriptional Regulator Influencing Cell Division in Haloferax volcanii

Transcriptional regulators that integrate cellular and environmental signals to control cell division are well known in bacteria and eukaryotes, but their existence is poorly understood in archaea. We identified a conserved gene (cdrS) that encodes a small protein and is highly transcribed in the model archaeon Haloferax volcanii. The cdrS gene could not be deleted, but CRISPR interference (CRISPRi)-mediated repression of the cdrS gene caused slow growth and cell division defects and changed the expression of multiple genes and their products associated with cell division, protein degradation, and metabolism. Consistent with this complex regulatory network, overexpression of cdrS inhibited cell division, whereas overexpression of the operon encoding both CdrS and a tubulin-like cell division protein (FtsZ2) stimulated division. Chromatin immunoprecipitation-DNA sequencing (ChIP-Seq) identified 18 DNA-binding sites of the CdrS protein, including one upstream of the promoter for a cell division gene, ftsZ1, and another upstream of the essential gene dacZ, encoding diadenylate cyclase involved in c-di-AMP signaling, which is implicated in the regulation of cell division. These findings suggest that CdrS is a transcription factor that plays a central role in a regulatory network coordinating metabolism and cell division. IMPORTANCE Cell division is a central mechanism of life and is essential for growth and development. Members of the Bacteria and Eukarya have different mechanisms for cell division, which have been studied in detail. In contrast, cell division in members of the Archaea is still understudied, and its regulation is poorly understood. Interestingly, different cell division machineries appear in members of the Archaea, with the Euryarchaeota using a cell division apparatus based on the tubulin-like cytoskeletal protein FtsZ, as in bacteria. Here, we identify the small protein CdrS as essential for survival and a central regulator of cell division in the euryarchaeon Haloferax volcanii. CdrS also appears to coordinate other cellular pathways, including synthesis of signaling molecules and protein degradation. Our results show that CdrS plays a sophisticated role in cell division, including regulation of numerous associated genes. These findings are expected to initiate investigations into conditional regulation of division in archaea.

Bottom-up and top-down proteomic approaches for the identification, characterization, and quantification of the low molecular weight proteome with focus on short open reading frame-encoded peptides

The recent discovery of alternative open reading frames creates a need for suitable analytical approaches to verify their translation and to characterize the corresponding gene products at the molecular level. As the analysis of small proteins within a background proteome by means of classical bottom-up proteomics is challenging, method development for the analysis of small open reading frame encoded peptides (SEPs) have become a focal point for research. Here, we highlight bottom-up and top-down proteomics approaches established for the analysis of SEPs in both pro- and eukaryotes. Major steps of analysis, including sample preparation and (small) proteome isolation, separation and mass spectrometry, data interpretation and quality control, quantification, the analysis of post-translational modifications, and exploration of functional aspects of the SEPs by means of proteomics technologies are described. These methods do not exclusively cover the analytics of SEPs but simultaneously include the low molecular weight proteome, and moreover, can also be used for the proteome-wide analysis of proteolytic processing events.

Discovery of a small protein factor involved in the coordinated degradation of phycobilisomes in cyanobacteria

During genome analysis, genes encoding small proteins are frequently neglected. Accordingly, small proteins have remained underinvestigated in all domains of life. Based on a previous systematic search for such genes, we present the functional analysis of the 66 amino acids protein NblD in a photosynthetic cyanobacterium. We show that NblD plays a crucial role during the coordinated dismantling of phycobilisome light-harvesting complexes. This disassembly is triggered when the cells become starved for nitrogen, a condition that frequently occurs in nature. Similar to NblA that tags phycobiliproteins for proteolysis, NblD binds to phycocyanin polypeptides but has a different function. The results show that, even in a well-investigated process, crucial new players can be discovered if small proteins are taken into consideration.

Phycobilisomes are the major pigment–protein antenna complexes that perform photosynthetic light harvesting in cyanobacteria, rhodophyte, and glaucophyte algae. Up to 50% of the cellular nitrogen can be stored in their giant structures. Accordingly, upon nitrogen depletion, phycobilisomes are rapidly degraded following an intricate genetic program. Here, we describe the role of NblD, a cysteine-rich, small protein in this process in cyanobacteria. Deletion of the nblD gene in the cyanobacterium Synechocystis sp. PCC 6803 prevented the degradation of phycobilisomes, leading to a nonbleaching (nbl) phenotype, which could be complemented by a plasmid-localized gene copy. Competitive growth experiments between the ΔnblD and the wild-type strain provided direct evidence for the physiological importance of NblD under nitrogen-limited conditions. Ectopic expression of NblD under nitrogen-replete conditions showed no effect, in contrast to the unrelated proteolysis adaptors NblA1 and NblA2, which can trigger phycobilisome degradation. Transcriptome analysis indicated increased nblA1/2 transcript levels in the ΔnblD strain during nitrogen starvation, implying that NblD does not act as a transcriptional (co)regulator. However, immunoprecipitation and far-western experiments identified the chromophorylated (holo form) of the phycocyanin β-subunit (CpcB) as its target, while apo-CpcB was not bound. The addition of recombinant NblD to isolated phycobilisomes caused a reduction in phycocyanin absorbance and a broadening and shifting of the peak to lower wavelengths, indicating the occurrence of structural changes. These data demonstrate that NblD plays a crucial role in the coordinated dismantling of phycobilisomes and add it as a factor to the genetically programmed response to nitrogen starvation.

smORFer: a modular algorithm to detect small ORFs in prokaryotes

Emerging evidence places small proteins (≤50 amino acids) more centrally in physiological processes. Yet, their functional identification and the systematic genome annotation of their cognate small open-reading frames (smORFs) remains challenging both experimentally and computationally. Ribosome profiling or Ribo-Seq (that is a deep sequencing of ribosome-protected fragments) enables detecting of actively translated open-reading frames (ORFs) and empirical annotation of coding sequences (CDSs) using the in-register translation pattern that is characteristic for genuinely translating ribosomes. Multiple identifiers of ORFs that use the 3-nt periodicity in Ribo-Seq data sets have been successful in eukaryotic smORF annotation. They have difficulties evaluating prokaryotic genomes due to the unique architecture (e.g. polycistronic messages, overlapping ORFs, leaderless translation, non-canonical initiation etc.). Here, we present a new algorithm, smORFer, which performs with high accuracy in prokaryotic organisms in detecting putative smORFs. The unique feature of smORFer is that it uses an integrated approach and considers structural features of the genetic sequence along with in-frame translation and uses Fourier transform to convert these parameters into a measurable score to faithfully select smORFs. The algorithm is executed in a modular way, and dependent on the data available for a particular organism, different modules can be selected for smORF search.

The small DUF1127 protein CcaF1 from Rhodobacter sphaeroides is an RNA-binding protein involved in sRNA maturation and RNA turnover

Many different protein domains are conserved among numerous species, but their function remains obscure. Proteins with DUF1127 domains number >17 000 in current databases, but a biological function has not yet been assigned to any of them. They are mostly found in alpha- and gammaproteobacteria, some of them plant and animal pathogens, symbionts or species used in industrial applications. Bioinformatic analyses revealed similarity of the DUF1127 domain of bacterial proteins to the RNA binding domain of eukaryotic Smaug proteins that are involved in RNA turnover and have a role in development from Drosophila to mammals. This study demonstrates that the 71 amino acid DUF1127 protein CcaF1 from the alphaproteobacterium Rhodobacter sphaeroides participates in maturation of the CcsR sRNAs that are processed from the 3' UTR of the ccaF mRNA and have a role in the oxidative stress defense. CcaF1 binds to many cellular RNAs of different type, several mRNAs with a function in cysteine / methionine / sulfur metabolism. It affects the stability of the CcsR RNAs and other non-coding RNAs and mRNAs. Thus, the widely distributed DUF1127 domain can mediate RNA-binding, affect stability of its binding partners and consequently modulate the bacterial transcriptome, thereby influencing different physiological processes.

The largely unexplored biology of small proteins in pro- and eukaryotes

The large abundance of small open reading frames (smORFs) in prokaryotic and eukaryotic genomes and the plethora of smORF-encoded small proteins became only apparent with the constant advancements in bioinformatic, genomic, proteomic, and biochemical tools. Small proteins are typically defined as proteins of < 50 amino acids in prokaryotes and of less than 100 amino acids in eukaryotes, and their importance for cell physiology and cellular adaptation is only beginning to emerge. In contrast to antimicrobial peptides, which are secreted by prokaryotic and eukaryotic cells for combatting pathogens and competitors, small proteins act within the producing cell mainly by stabilizing protein assemblies and by modifying the activity of larger proteins. Production of small proteins is frequently linked to stress conditions or environmental changes, and therefore, cells seem to use small proteins as intracellular modifiers for adjusting cell metabolism to different intra- and extracellular cues. However, the size of small proteins imposes a major challenge for the cellular machinery required for protein folding and intracellular trafficking and recent data indicate that small proteins can engage distinct trafficking pathways. In the current review, we describe the diversity of small proteins in prokaryotes and eukaryotes, highlight distinct and common features, and illustrate how they are handled by the protein trafficking machineries in prokaryotic and eukaryotic cells. Finally, we also discuss future topics of research on this fascinating but largely unexplored group of proteins.

Phycobilisome breakdown effector NblD is required to maintain the cellular amino acid composition during nitrogen starvation

Small proteins are critically involved in the acclimation response of photosynthetic cyanobacteria to nitrogen starvation. NblD is the 66-amino-acid effector of nitrogen-limitation-induced phycobilisome breakdown, which is believed to replenish the cellular amino acid pools. To address the physiological functions of NblD, the concentrations of amino acids, intermediates of the arginine catabolism pathway and several organic acids were measured during the response to nitrogen starvation in the cyanobacterium Synechocystis sp. PCC 6803 wild type and in an nblD deletion strain. A characteristic signature of metabolite pool composition was identified, which shows that NblD-mediated phycobilisome degradation is required to maintain the cellular amino acid and organic acid pools during nitrogen starvation. Specific deviations from the wild type suggest wider-reaching effects that also affect such processes as redox homeostasis via glutathione and tetrapyrrole biosynthesis, both of which are linked to the strongly decreased glutamate pool, and transcriptional reprogramming via an enhanced concentration of 2-oxoglutarate, the metabolite co-regulator of the NtcA transcription factor. The essential role played by NblD in metabolic homeostasis is consistent with the widespread occurrence of NblD throughout the cyanobacterial radiation and the previously observed strong positive selection for the nblD gene under fluctuating nitrogen supply. Importance Cyanobacteria play important roles in the global carbon and nitrogen cycles. In their natural environment, these organisms are exposed to fluctuating nutrient conditions. Nitrogen starvation induces a coordinated nitrogen-saving program that includes the breakdown of nitrogen-rich photosynthetic pigments, particularly phycobiliproteins. The small protein NblD was recently identified as an effector of phycobilisome breakdown in cyanobacteria. In this study, we demonstrate that the NblD-mediated degradation of phycobiliproteins is needed to sustain cellular pools of soluble amino acids and other crucial metabolites. The essential role played by NblD in metabolic homeostasis explains why genes encoding this small protein are conserved in almost all members of cyanobacterial radiation.

Peptides encoded in the Streptococcus mutans RcrRPQ operon are essential for thermotolerance

The MarR-like transcriptional regulator and two ABC transporters encoded by the rcrRPQ operon in the dental caries pathogen Streptococcus mutans have important regulatory roles related to oxidative stress tolerance, genetic competence and (p)ppGpp metabolism. A unique feature of the rcrRPQ operon, when compared to other bacteria, is the presence of two peptides, designated Pep1 and Pep2, encoded in alternative reading frames at the 3' end of rcrQ. Here, we show that the rcrRPQ operon, including Pep1 and 2, is essential for S. mutans to survive and maintain viability at elevated temperatures. No major changes in the levels of the heat shock proteins DnaK or GroEL that could account for the thermosensitivity of rcrRPQ mutants were observed. By introducing a single amino acid substitution into the comX gene that deletes an internally encoded peptide, XrpA, we found that XrpA is a contributing factor to the thermosensitive phenotype of a ΔrcrR strain. Overexpression of XrpA on a plasmid also caused a significant growth defect at 42 °C. Interestingly, loss of the gene for the RelA/SpoT homologue (RSH) enzyme, relA, restored growth of the ΔrcrR strain at 42 °C. During heat stress and when a stringent response was induced, levels of (p)ppGpp were elevated in the ΔrcrR strain. Deletion of relA in the ΔrcrR strain lowered the basal levels of (p)ppGpp to those observed in wild-type S. mutans. Thus, (p)ppGpp pools are dysregulated in ΔrcrR, which likely leads to aberrant control of transcriptional/translational processes and the thermosensitive phenotype. In summary, the genes and peptides encoded in the rcrRPQ operon are critical for thermotolerance, and in some strains these phenotypes are related to altered (p)ppGpp metabolism and increased production of the XrpA peptide.

The Small Toxic Salmonella Protein TimP Targets the Cytoplasmic Membrane and Is Repressed by the Small RNA TimR

Next-generation sequencing (NGS) has enabled the revelation of a vast number of genomes from organisms spanning all domains of life. To reduce complexity when new genome sequences are annotated, open reading frames (ORFs) shorter than 50 codons in length are generally omitted. However, it has recently become evident that this procedure sorts away ORFs encoding small proteins of high biological significance. For instance, tailored small protein identification approaches have shown that bacteria encode numerous small proteins with important physiological functions. As the number of predicted small ORFs increase, it becomes important to characterize the corresponding proteins. In this study, we discovered a conserved but previously overlooked small enterobacterial protein. We show that this protein, which we dubbed TimP, is a potent toxin that inhibits bacterial growth by targeting the cell membrane. Toxicity is relieved by a small regulatory RNA, which binds the toxin mRNA to inhibit toxin synthesis.

Small proteins are gaining increased attention due to their important functions in major biological processes throughout the domains of life. However, their small size and low sequence conservation make them difficult to identify. It is therefore not surprising that enterobacterial ryfA has escaped identification as a small protein coding gene for nearly 2 decades. Since its identification in 2001, ryfA has been thought to encode a noncoding RNA and has been implicated in biofilm formation in Escherichia coli and pathogenesis in Shigella dysenteriae. Although a recent ribosome profiling study suggested ryfA to be translated, the corresponding protein product was not detected. In this study, we provide evidence that ryfA encodes a small toxic inner membrane protein, TimP, overexpression of which causes cytoplasmic membrane leakage. TimP carries an N-terminal signal sequence, indicating that its membrane localization is Sec-dependent. Expression of TimP is repressed by the small RNA (sRNA) TimR, which base pairs with the timP mRNA to inhibit its translation. In contrast to overexpression, endogenous expression of TimP upon timR deletion permits cell growth, possibly indicating a toxicity-independent function in the bacterial membrane.

Posttranslational insertion of small membrane proteins by the bacterial signal recognition particle

Small membrane proteins represent a largely unexplored yet abundant class of proteins in pro- and eukaryotes. They essentially consist of a single transmembrane domain and are associated with stress response mechanisms in bacteria. How these proteins are inserted into the bacterial membrane is unknown. Our study revealed that in Escherichia coli, the 27-amino-acid-long model protein YohP is recognized by the signal recognition particle (SRP), as indicated by in vivo and in vitro site-directed cross-linking. Cross-links to SRP were also observed for a second small membrane protein, the 33-amino-acid-long YkgR. However, in contrast to the canonical cotranslational recognition by SRP, SRP was found to bind to YohP posttranslationally. In vitro protein transport assays in the presence of a SecY inhibitor and proteoliposome studies demonstrated that SRP and its receptor FtsY are essential for the posttranslational membrane insertion of YohP by either the SecYEG translocon or by the YidC insertase. Furthermore, our data showed that the yohP mRNA localized preferentially and translation-independently to the bacterial membrane in vivo. In summary, our data revealed that YohP engages an unique SRP-dependent posttranslational insertion pathway that is likely preceded by an mRNA targeting step. This further highlights the enormous plasticity of bacterial protein transport machineries.

Small membrane proteins represent a largely unexplored yet abundant class of proteins, but how they are inserted into the bacterial membrane is unknown. This study identifies a novel posttranslational protein transport pathway that relies on the signal recognition particle and the SecYEG translocon/YidC insertase.

Heterogeneous resistance to colistin in Enterobacter cloacae complex due to a new small transmembrane protein

BACKGROUND

Enterobacter strains can display heterogeneous resistance (heteroresistance) to colistin but the mechanisms remain largely unknown. We investigated potential mechanisms of colistin heteroresistance in an Enterobacter clinical strain, WCHECl-1060, and found a new mechanism.

METHODS

Strain WCHECl-1060 was subjected to WGS to identify known colistin resistance mechanisms. Tn5 insertional mutagenesis, gene knockout and complementation and shotgun cloning were employed to investigate unknown colistin heteroresistance mechanisms. RNA sequencing was performed to link the newly identified mechanism with known ones.

RESULTS

We showed that the phoP gene [encoding part of the PhoP-PhoQ two-component system (TCS)], the dedA(Ecl) gene (encoding an inner membrane protein of the DedA family) and the tolC gene (encoding part of the AcrAB-TolC efflux pump) are required for colistin heteroresistance. We identified a new gene, ecr, encoding a 72 amino acid transmembrane protein, which was able to mediate colistin heteroresistance. We then performed RNA sequencing and transcriptome analysis and found that in the presence of ecr the expression of phoP and the arnBCADTEF operon, which synthesizes and transfers l-Ara4N to lipid A, was increased significantly.

CONCLUSIONS

The small protein encoded by ecr represents a new colistin heteroresistance mechanism and is likely to mediate colistin heteroresistance via the PhoP-PhoQ TCS to act on the arnBCADTEF operon.

SR7 - a dual-function antisense RNA from Bacillus subtilis

Here, we describe SR7, a dual-function antisense RNA encoded on the Bacillus subtilis chromosome. This RNA was earlier described as SigB-dependent regulatory RNA S1136 and reported to reduce the amount of the small ribosomal subunit under ethanol stress. We found that the 5ʹ portion of SR7 encodes a small protein composed of 39 amino acids which we designated SR7P. It is translated from a 185 nt SigB-dependent mRNA under five different stress conditions and a longer SigB-independent RNA constitutively. About three-fold higher amounts of SR7P were detected in B. subtilis cells exposed to salt, ethanol, acid or heat stress. Co-elution experiments with SR7P_C-FLAG and Far-Western blotting demonstrated that SR7P interacts with the glycolytic enzyme enolase. Enolase is a scaffolding component of the B. subtilis degradosome where it interacts with RNase Y and phosphofructokinase PfkA. We found that SR7P increases the amount of RNase Y bound to enolase without affecting PfkA. RNA does not bridge the SR7P-enolase-RNase Y interaction. In vitro-degradation assays with the known RNase Y substrates yitJ and rpsO mRNA revealed enhanced enzymatic activity of enolase-bound RNase Y in the presence of SR7P. Northern blots showed a major effect of enolase and a minor effect of SR7P on the half-life of rpsO mRNA indicating a fine-tuning role of SR7P in RNA degradation.

Polycysteine-encoding leaderless short ORFs function as cysteine-responsive attenuators of operonic gene expression in mycobacteria

Genome-wide transcriptomic analyses have revealed abundant expressed short open reading frames (ORFs) in bacteria. Whether these short ORFs, or the small proteins they encode, are functional remains an open question. One quarter of mycobacterial mRNAs are leaderless, beginning with a 5'-AUG or GUG initiation codon. Leaderless mRNAs often encode unannotated short ORFs as the first gene of a polycistronic transcript. Here, we show that polycysteine-encoding leaderless short ORFs function as cysteine-responsive attenuators of operonic gene expression. Detailed mutational analysis shows that one polycysteine short ORF controls expression of the downstream genes. Our data indicate that ribosomes stalled in the polycysteine tract block mRNA structures that otherwise sequester the ribosome-binding site of the 3'gene. We assessed endogenous proteomic responses to cysteine limitation in Mycobacterium smegmatis using mass spectrometry. Six cysteine metabolic loci having unannotated polycysteine-encoding leaderless short ORF architectures responded to cysteine limitation, revealing widespread cysteine-responsive attenuation in mycobacteria. Individual leaderless short ORFs confer independent operon-level control, while their shared dependence on cysteine ensures a collective response mediated by ribosome pausing. We propose the term ribulon to classify ribosome-directed regulons. Regulon-level coordination by ribosomes on sensory short ORFs illustrates one utility of the many unannotated short ORFs expressed in bacterial genomes.

AcnSP - A Novel Small Protein Regulator of Aconitase Activity in the Cyanobacterium Synechocystis sp. PCC 6803

Synechocystis sp. PCC 6803 is a widely used model cyanobacterium whose genome has been well annotated. However, several additional small protein coding sequences (sORFs) have been recently identified, which might play important roles, for example in the regulation of cellular metabolism. Here, we analyzed the function of a sORF encoding a 44 amino acid peptide showing high similarity to the N-terminal part of aconitase (AcnB). The expression of the gene, which probably originated from a partial gene duplication of chromosomal acnB into the plasmid pSYSA, was verified and it was designated as acnSP. The protein-coding part of acnSP was inactivated by interposon mutagenesis. The obtained mutant displayed slower growth under photoautotrophic conditions with light exceeding 100 μmol photons m^–2 s^–1 and showed significant changes in the metabolome compared to wild type, including alterations in many metabolites associated to the tricarboxylic acid (TCA) cycle. To analyze a possible direct impact of AcnSP on aconitase, the recombinant Synechocystis enzyme was generated and biochemically characterized. Biochemical analysis revealed that addition of equimolar amounts of AcnSP resulted in an improved substrate affinity (lower K_m) and lowered V_max of aconitase. These results imply that AcnSP can regulate aconitase activity, thereby impacting the carbon flow into the oxidative branch of the cyanobacterial TCA cycle, which is mainly responsible for the synthesis of carbon skeletons needed for ammonia assimilation.

OCCAM: prediction of small ORFs in bacterial genomes by means of a target-decoy database approach and machine learning techniques

Small open reading frames (ORFs) have been systematically disregarded by automatic genome annotation. The difficulty in finding patterns in tiny sequences is the main reason that makes small ORFs to be overlooked by computational procedures. However, advances in experimental methods show that small proteins can play vital roles in cellular activities. Hence, it is urgent to make progress in the development of computational approaches to speed up the identification of potential small ORFs. In this work, our focus is on bacterial genomes. We improve a previous approach to identify small ORFs in bacteria. Our method uses machine learning techniques and decoy subject sequences to filter out spurious ORF alignments. We show that an advanced multivariate analysis can be more effective in terms of sensitivity than applying the simplistic and widely used e-value cutoff. This is particularly important in the case of small ORFs for which alignments present higher e-values than usual. Experiments with control datasets show that the machine learning algorithms used in our method to curate significant alignments can achieve average sensitivity and specificity of 97.06% and 99.61%, respectively. Therefore, an important step is provided here toward the construction of more accurate computational tools for the identification of small ORFs in bacteria.

Genome-wide discovery of structured noncoding RNAs in bacteria

Background

Structured noncoding RNAs (ncRNAs) play essential roles in many biological processes such as gene regulation, signaling, RNA processing, and protein synthesis. Among the most common groups of ncRNAs in bacteria are riboswitches. These cis-regulatory, metabolite-binding RNAs are present in many species where they regulate various metabolic and signaling pathways. Collectively, there are likely to be hundreds of novel riboswitch classes that remain hidden in the bacterial genomes that have already been sequenced, and potentially thousands of classes distributed among various other species in the biosphere. The vast majority of these undiscovered classes are proposed to be exceedingly rare, and so current bioinformatics search techniques are reaching their limits for differentiating between true riboswitch candidates and false positives.

Results

Herein, we exploit a computational search pipeline that can efficiently identify intergenic regions most likely to encode structured ncRNAs. Application of this method to five bacterial genomes yielded nearly 70 novel genetic elements including 30 novel candidate ncRNA motifs. Among the riboswitch candidates identified is an RNA motif involved in the regulation of thiamin biosynthesis.

Conclusions

Analysis of other genomes will undoubtedly lead to the discovery of many additional novel structured ncRNAs, and provide insight into the range of riboswitches and other kinds of ncRNAs remaining to be discovered in bacteria and archaea.

Electronic supplementary material

The online version of this article (10.1186/s12866-019-1433-7) contains supplementary material, which is available to authorized users.

Mining for missed sORF-encoded peptides

INTRODUCTION

Small open reading frames (sORFs) with potential protein-coding capacity have been disclosed in various transcripts, including long noncoding RNAs (LncRNAs), mRNAs (5'-upstream, coding domain, and 3'-downstream), circular RNAs, pri-miRNAs, and ribosomal RNAs (rRNAs). Recent characterization of several sORF-encoded peptides (SEPs or micropeptides) revealed their important roles in many fundamental biological processes in a broad range of species from yeast to human. The success in the mining of micropeptides attributes to the advanced bioinformatics and high-throughput sequencing techniques. Areas covered: sORFs and SEPs were overlooked for their tiny size and the difficulty of identification by bioinformatics analyses. With more and more sORFs and SEPs have been identified, this field has attracted more attention. This review covers recent advances in the strategies for the detection and identification of sORFs and SEPs. Expert commentary: The advantages and drawbacks of the strategies for detection and identification of sORFs and SEPs are discussed, as well as the techniques that are used to decipher the roles of micropeptides in organisms are described.

Regulatory Tools for Controlling Gene Expression in Cyanobacteria

Cyanobacteria are attractive hosts for converting carbon dioxide and sunlight into desirable chemical products. To engineer these organisms and manipulate their metabolic pathways, the biotechnology community has developed genetic tools to control gene expression. Many native cyanobacterial promoters and related sequence elements have been used to regulate genes of interest, and heterologous tools that use non-native small molecules to induce gene expression have been demonstrated. Overall, IPTG-based induction systems seem to be leaky and initially demonstrate small dynamic ranges in cyanobacteria. Consequently, a variety of other induction systems have been optimized to enable tighter control of gene expression. Tools require significant optimization because they function quite differently in cyanobacteria when compared to analogous use in model heterotrophs. We hypothesize that these differences are due to fundamental differences in physiology between organisms. This review is not intended to summarize all known products made in cyanobacteria nor the performance (titer, rate, yield) of individual strains, but instead will focus on the genetic tools and the inherent aspects of cellular physiology that influence gene expression in cyanobacteria.

A Hydrophobic Small Protein, BpOF4_01690, Is Critical for Alkaliphily of Alkaliphilic Bacillus pseudofirmus OF4

A monocistronic small protein, BpOF4_01690, was annotated in alkaliphilic Bacillus pseudofirmus OF4. It comprises 59 amino acids and is hydrophobic. Importantly, homologs of this protein were identified only in alkaliphiles. In this study, a mutant with a BpOF4_01690 gene deletion (designated Δ01690) exhibited weaker growth than that of the wild type in both malate-based defined and glucose-based defined media under low-sodium conditions at pH 10.5. Additionally, the enzymatic activity of the respiratory chain of Δ01690 was much lower than that of the wild type. These phenotypes were similar to those of a ctaD deletion mutant and an atpB-F deletion mutant. Therefore, we hypothesize that BpOF4_01690 plays a critical role in oxidative phosphorylation under highly alkaline conditions.

Application of targeted mass spectrometry in bottom-up proteomics for systems biology research

The enormous diversity of proteoforms produces tremendous complexity within cellular proteomes, facilitates intricate networks of molecular interactions, and constitutes a formidable analytical challenge for biomedical researchers. Currently, quantitative whole-proteome profiling often relies on non-targeted liquid chromatography-mass spectrometry (LC-MS), which samples proteoforms broadly, but can suffer from lower accuracy, sensitivity, and reproducibility compared with targeted LC-MS. Recent advances in bottom-up proteomics using targeted LC-MS have enabled previously unachievable identification and quantification of target proteins and posttranslational modifications within complex samples. Consequently, targeted LC-MS is rapidly advancing biomedical research, especially systems biology research in diverse areas that include proteogenomics, interactomics, kinomics, and biological pathway modeling. With the recent development of targeted LC-MS assays for nearly the entire human proteome, targeted LC-MS is positioned to enable quantitative proteomic profiling of unprecedented quality and accessibility to support fundamental and clinical research. Here we review recent applications of bottom-up proteomics using targeted LC-MS for systems biology research. SIGNIFICANCE: Advances in targeted proteomics are rapidly advancing systems biology research. Recent applications include systems-level investigations focused on posttranslational modifications (such as phosphoproteomics), protein conformation, protein-protein interaction, kinomics, proteogenomics, and metabolic and signaling pathways. Notably, absolute quantification of metabolic and signaling pathway proteins has enabled accurate pathway modeling and engineering. Integration of targeted proteomics with other technologies, such as RNA-seq, has facilitated diverse research such as the identification of hundreds of "missing" human proteins (genes and transcripts that appear to encode proteins but direct experimental evidence was lacking).

Small bacterial and phagic proteins: an updated view on a rapidly moving field

Small proteins, that is, polypeptides of 50 amino acids (aa) or less, are increasingly recognized as important regulators in bacteria. Secreted or not, their small size make them versatile proteins, involved in a wide range of processes. They may allow bacteria to sense and to respond to stresses, to send signals and communicate, and to modulate infections. Bacteriophages also produce small proteins to influence lysogeny/lysis decisions. In this review, we update the present view on small proteins functions, and discuss their possible applications.

Functional analysis of a biosynthetic cluster essential for production of 4-formylaminooxyvinylglycine, a germination-arrest factor from Pseudomonas fluorescens WH6

Rhizosphere-associated Pseudomonas fluorescens WH6 produces the germination-arrest factor 4-formylaminooxyvinylglycine (FVG). FVG has previously been shown to both arrest the germination of weedy grasses and inhibit the growth of the bacterial plant pathogen Erwinia amylovora. Very little is known about the mechanism by which FVG is produced. Although a previous study identified a region of the genome that may be involved in FVG biosynthesis, it has not yet been determined which genes within that region are sufficient and necessary for FVG production. In the current study, we explored the role of each of the putative genes encoded in that region by constructing deletion mutations. Mutant strains were assayed for their ability to produce FVG with a combination of biological assays and TLC analyses. This work defined the core FVG biosynthetic gene cluster and revealed several interesting characteristics of FVG production. We determined that FVG biosynthesis requires two small ORFs of less than 150 nucleotides and that multiple transporters have overlapping but distinct functionality. In addition, two genes in the centre of the biosynthetic gene cluster are not required for FVG production, suggesting that additional products may be produced from the cluster. Transcriptional analysis indicated that at least three active promoters play a role in the expression of genes within this cluster. The results of this study enrich our knowledge regarding the diversity of mechanisms by which bacteria produce non-proteinogenic amino acids like vinylglycines.

The Conserved Spore Coat Protein SpoVM Is Largely Dispensable in Clostridium difficile Spore Formation

The spore-forming obligate anaerobe Clostridium difficile is the leading cause of antibiotic-associated diarrheal disease in the United States. When C. difficile spores are ingested by susceptible individuals, they germinate within the gut and transform into vegetative, toxin-secreting cells. During infection, C. difficile must also induce spore formation to survive exit from the host. Since spore formation is essential for transmission, understanding the basic mechanisms underlying sporulation in C. difficile could inform the development of therapeutic strategies targeting spores. In this study, we determine the requirement of the C. difficile homolog of SpoVM, a protein that is essential for spore formation in Bacillus subtilis due to its regulation of coat and cortex formation. We observed that SpoVM plays a minor role in C. difficile spore formation, in contrast with B. subtilis, indicating that this protein would not be a good target for inhibiting spore formation.

The spore-forming bacterial pathogen Clostridium difficile is a leading cause of health care-associated infections in the United States. In order for this obligate anaerobe to transmit infection, it must form metabolically dormant spores prior to exiting the host. A key step during this process is the assembly of a protective, multilayered proteinaceous coat around the spore. Coat assembly depends on coat morphogenetic proteins recruiting distinct subsets of coat proteins to the developing spore. While 10 coat morphogenetic proteins have been identified in Bacillus subtilis, only two of these morphogenetic proteins have homologs in the Clostridia: SpoIVA and SpoVM. C. difficile SpoIVA is critical for proper coat assembly and functional spore formation, but the requirement for SpoVM during this process was unknown. Here, we show that SpoVM is largely dispensable for C. difficile spore formation, in contrast with B. subtilis. Loss of C. difficile SpoVM resulted in modest decreases (~3-fold) in heat- and chloroform-resistant spore formation, while morphological defects such as coat detachment from the forespore and abnormal cortex thickness were observed in ~30% of spoVM mutant cells. Biochemical analyses revealed that C. difficile SpoIVA and SpoVM directly interact, similarly to their B. subtilis counterparts. However, in contrast with B. subtilis, C. difficile SpoVM was not essential for SpoIVA to encase the forespore. Since C. difficile coat morphogenesis requires SpoIVA-interacting protein L (SipL), which is conserved exclusively in the Clostridia, but not the more broadly conserved SpoVM, our results reveal another key difference between C. difficile and B. subtilis spore assembly pathways.

IMPORTANCE The spore-forming obligate anaerobe Clostridium difficile is the leading cause of antibiotic-associated diarrheal disease in the United States. When C. difficile spores are ingested by susceptible individuals, they germinate within the gut and transform into vegetative, toxin-secreting cells. During infection, C. difficile must also induce spore formation to survive exit from the host. Since spore formation is essential for transmission, understanding the basic mechanisms underlying sporulation in C. difficile could inform the development of therapeutic strategies targeting spores. In this study, we determine the requirement of the C. difficile homolog of SpoVM, a protein that is essential for spore formation in Bacillus subtilis due to its regulation of coat and cortex formation. We observed that SpoVM plays a minor role in C. difficile spore formation, in contrast with B. subtilis, indicating that this protein would not be a good target for inhibiting spore formation.

Common and phylogenetically widespread coding for peptides by bacterial small RNAs

Background

While eukaryotic noncoding RNAs have recently received intense scrutiny, it is becoming clear that bacterial transcription is at least as pervasive. Bacterial small RNAs and antisense RNAs (sRNAs) are often assumed to be noncoding, due to their lack of long open reading frames (ORFs). However, there are numerous examples of sRNAs encoding for small proteins, whether or not they also have a regulatory role at the RNA level.

Methods

Here, we apply flexible machine learning techniques based on sequence features and comparative genomics to quantify the prevalence of sRNA ORFs under natural selection to maintain protein-coding function in 14 phylogenetically diverse bacteria. Importantly, we quantify uncertainty in our predictions, and follow up on them using mass spectrometry proteomics and comparison to datasets including ribosome profiling.

Results

A majority of annotated sRNAs have at least one ORF between 10 and 50 amino acids long, and we conservatively predict that 409±191.7 unannotated sRNA ORFs are under selection to maintain coding (mean estimate and 95% confidence interval), an average of 29 per species considered here. This implies that overall at least 10.3±0.5% of sRNAs have a coding ORF, and in some species around 20% do. 165±69 of these novel coding ORFs have some antisense overlap to annotated ORFs. As experimental validation, many of our predictions are translated in published ribosome profiling data and are identified via mass spectrometry shotgun proteomics. B. subtilis sRNAs with coding ORFs are enriched for high expression in biofilms and confluent growth, and S. pneumoniae sRNAs with coding ORFs are involved in virulence. sRNA coding ORFs are enriched for transmembrane domains and many are predicted novel components of type I toxin/antitoxin systems.

Conclusions

We predict over two dozen new protein-coding genes per bacterial species, but crucially also quantified the uncertainty in this estimate. Our predictions for sRNA coding ORFs, along with predicted novel type I toxins and tools for sorting and visualizing genomic context, are freely available in a user-friendly format at http://disco-bac.web.pasteur.fr. We expect these easily-accessible predictions to be a valuable tool for the study not only of bacterial sRNAs and type I toxin-antitoxin systems, but also of bacterial genetics and genomics.

Electronic supplementary material

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

Construction of a fur null mutant and RNA-sequencing provide deeper global understanding of the Aliivibrio salmonicida Fur regulon

BACKGROUND

The ferric uptake regulator (Fur) is a transcription factor and the main regulator of iron acquisition in prokaryotes. When bound to ferric iron, Fur recognizes its DNA binding site and generally executes its function by repressing transcription of its target genes. Due to its importance in virulence, the Fur regulon is well studied for several model bacteria. In our previous work, we used computational predictions and microarray to gain insights into Fur-regulation in Aliivibrio salmonicida, and have identified a number of genes and operons that appear to be under direct control of Fur. To provide a more accurate and deeper global understanding of the biological role of Fur we have now generated an A. salmonicida fur knock-out strain and used RNA-sequencing to compare gene expression between the wild-type and fur null mutant strains.

RESULTS

An A. salmonicida fur null mutant strain was constructed. Biological assays demonstrate that deletion of fur results in loss of fitness, with reduced growth rates, and reduced abilities to withstand low-iron conditions, and oxidative stress. When comparing expression levels in the wild-type and the fur null mutant we retrieved 296 differentially expressed genes distributed among 18 of 21 functional classes of genes. A gene cluster encoding biosynthesis of the siderophore bisucaberin represented the highest up-regulated genes in the fur null mutant. Other highly up-regulated genes all encode proteins important for iron acquisition. Potential targets for the RyhB sRNA was predicted from the list of down-regulated genes, and significant complementarities were found between RyhB and mRNAs of the fur, sodB, cysN and VSAL_I0422 genes. Other sRNAs with potential functions in iron homeostasis were identified.

CONCLUSION

The present work provides by far the most comprehensive and deepest understanding of the Fur regulon in A. salmonicida to date. Our data also contribute to a better understanding of how Fur plays a key role in iron homeostasis in bacteria in general, and help to show how Fur orchestrates iron uptake when iron levels are extremely low.

Pleiotropic Regulation of Virulence Genes in Streptococcus mutans by the Conserved Small Protein SprV

Streptococcus mutans, an oral pathogen associated with dental caries, colonizes tooth surfaces as polymicrobial biofilms known as dental plaque. S. mutans expresses several virulence factors that allow the organism to tolerate environmental fluctuations and compete with other microorganisms. We recently identified a small hypothetical protein (90 amino acids) essential for the normal growth of the bacterium. Inactivation of the gene, SMU.2137, encoding this protein caused a significant growth defect and loss of various virulence-associated functions. An S. mutans strain lacking this gene was more sensitive to acid, temperature, osmotic, oxidative, and DNA damage-inducing stresses. In addition, we observed an altered protein profile and defects in biofilm formation, bacteriocin production, and natural competence development, possibly due to the fitness defect associated with SMU.2137 deletion. Transcriptome sequencing revealed that nearly 20% of the S. mutans genes were differentially expressed upon SMU.2137 deletion, thereby suggesting a pleiotropic effect. Therefore, we have renamed this hitherto uncharacterized gene as sprV (streptococcal pleiotropic regulator of virulence). The transcript levels of several relevant genes in the sprV mutant corroborated the phenotypes observed upon sprV deletion. Owing to its highly conserved nature, inactivation of the sprV ortholog in Streptococcus gordonii also resulted in poor growth and defective UV tolerance and competence development as in the case of S. mutans Our experiments suggest that SprV is functionally distinct from its homologs identified by structure and sequence homology. Nonetheless, our current work is aimed at understanding the importance of SprV in the S. mutans biology.IMPORTANCEStreptococcus mutans employs several virulence factors and stress resistance mechanisms to colonize tooth surfaces and cause dental caries. Bacterial pathogenesis is generally controlled by regulators of fitness that are critical for successful disease establishment. Sometimes these regulators, which are potential targets for antimicrobials, are lost in the genomic context due to the lack of annotated homologs. This work outlines the regulatory impact of a small, highly conserved hypothetical protein, SprV, encoded by S. mutans We show that SprV affects the transcript levels of various virulence factors required for normal growth, biofilm formation, stress tolerance, genetic competence, and bacteriocin production.

Differentiation of ncRNAs from small mRNAs in Escherichia coli O157:H7 EDL933 (EHEC) by combined RNAseq and RIBOseq - ryhB encodes the regulatory RNA RyhB and a peptide, RyhP

Background

While NGS allows rapid global detection of transcripts, it remains difficult to distinguish ncRNAs from short mRNAs. To detect potentially translated RNAs, we developed an improved protocol for bacterial ribosomal footprinting (RIBOseq). This allowed distinguishing ncRNA from mRNA in EHEC. A high ratio of ribosomal footprints per transcript (ribosomal coverage value, RCV) is expected to indicate a translated RNA, while a low RCV should point to a non-translated RNA.

Results

Based on their low RCV, 150 novel non-translated EHEC transcripts were identified as putative ncRNAs, representing both antisense and intergenic transcripts, 74 of which had expressed homologs in E. coli MG1655. Bioinformatics analysis predicted statistically significant target regulons for 15 of the intergenic transcripts; experimental analysis revealed 4-fold or higher differential expression of 46 novel ncRNA in different growth media. Out of 329 annotated EHEC ncRNAs, 52 showed an RCV similar to protein-coding genes, of those, 16 had RIBOseq patterns matching annotated genes in other enterobacteriaceae, and 11 seem to possess a Shine-Dalgarno sequence, suggesting that such ncRNAs may encode small proteins instead of being solely non-coding. To support that the RIBOseq signals are reflecting translation, we tested the ribosomal-footprint covered ORF of ryhB and found a phenotype for the encoded peptide in iron-limiting condition.

Conclusion

Determination of the RCV is a useful approach for a rapid first-step differentiation between bacterial ncRNAs and small mRNAs. Further, many known ncRNAs may encode proteins as well.

Electronic supplementary material

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

Small proteins in cyanobacteria provide a paradigm for the functional analysis of the bacterial micro-proteome

Background

Despite their versatile functions in multimeric protein complexes, in the modification of enzymatic activities, intercellular communication or regulatory processes, proteins shorter than 80 amino acids (μ-proteins) are a systematically underestimated class of gene products in bacteria. Photosynthetic cyanobacteria provide a paradigm for small protein functions due to extensive work on the photosynthetic apparatus that led to the functional characterization of 19 small proteins of less than 50 amino acids. In analogy, previously unstudied small ORFs with similar degrees of conservation might encode small proteins of high relevance also in other functional contexts.

Results

Here we used comparative transcriptomic information available for two model cyanobacteria, Synechocystis sp. PCC 6803 and Synechocystis sp. PCC 6714 for the prediction of small ORFs. We found 293 transcriptional units containing candidate small ORFs ≤80 codons in Synechocystis sp. PCC 6803, also including the known mRNAs encoding small proteins of the photosynthetic apparatus. From these transcriptional units, 146 are shared between the two strains, 42 are shared with the higher plant Arabidopsis thaliana and 25 with E. coli. To verify the existence of the respective μ-proteins in vivo, we selected five genes as examples to which a FLAG tag sequence was added and re-introduced them into Synechocystis sp. PCC 6803. These were the previously annotated gene ssr1169, two newly defined genes norf1 and norf4, as well as nsiR6(nitrogen stress-induced RNA 6) and hliR1(high light-inducible RNA 1) , which originally were considered non-coding. Upon activation of expression via the Cu^2+.responsive petE promoter or from the native promoters, all five proteins were detected in Western blot experiments.

Conclusions

The distribution and conservation of these five genes as well as their regulation of expression and the physico-chemical properties of the encoded proteins underline the likely great bandwidth of small protein functions in bacteria and makes them attractive candidates for functional studies.

Mutually exclusive RNA secondary structures regulate translation initiation of DinQ in Escherichia coli

Protein translation can be affected by changes in the secondary structure of mRNA. The dinQ gene in Escherichia coli encodes a primary transcript (+1) that is inert to translation. Ribonucleolytic removal of the 44 first nucleotides converts the +1 transcript into a translationally active form, but the mechanism behind this structural change is unknown. Here we present experimental evidence for a mechanism where alternative RNA secondary structures in the two dinQ mRNA variants affect translation initiation by mediating opening or closing of the ribosome binding sequence. This structural switch is determined by alternative interactions of four sequence elements within the dinQ mRNA and also by the agrB antisense RNA. Additionally, the structural conformation of +1 dinQ suggests a locking mechanism comprised of an RNA stem that both stabilizes and prevents translation initiation from the full-length dinQ transcript. BLAST search and multiple sequence alignments define a new family of dinQ-like genes widespread in Enterobacteriaceae with close RNA sequence similarities in their 5' untranslated regions. Thus, it appears that a whole new family of genes is regulated by the same mechanism of alternative secondary RNA structures.

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

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

Transcriptomic analysis of staphylococcal sRNAs: insights into species-specific adaption and the evolution of pathogenesis

Next-generation sequencing technologies have dramatically increased the rate at which new genomes are sequenced. Accordingly, automated annotation programs have become adept at identifying and annotating protein coding regions, as well as common and conserved RNAs. Additionally, RNAseq techniques have advanced our ability to identify and annotate regulatory RNAs (sRNAs), which remain significantly understudied. Recently, our group catalogued and annotated all previously known and newly identified sRNAs in several Staphylococcus aureus strains. These complete annotation files now serve as tools to compare the sRNA content of S. aureus with other bacterial strains to investigate the conservation of their sRNomes. Accordingly, in this study we performed RNAseq on two staphylococcal species, Staphylococcus epidermidis and Staphylococcus carnosus, identifying 118 and 89 sRNAs in these organisms, respectively. The sRNA contents of all three species were then compared to elucidate their common and species-specific sRNA content, identifying a core set of between 53 and 36 sRNAs encoded in each organism. In addition, we determined that S. aureus has the largest set of unique sRNAs (137) while S. epidermidishas the fewest (25). Finally, we identify a highly conserved sequence and structural motif differentially represented within, yet common to, both S. aureus and S. epidermidis. Collectively, in this study, we uncover the sRNome common to three staphylococcal species, shedding light on sRNAs that are likely to be involved in basic physiological processes common to the genus. More significantly, we have identified species-specific sRNAs that are likely to influence the individual lifestyle and behaviour of these diverse staphylococcal strains.

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.

Towards the complete small RNome of Acinetobacter baumannii

In recent years, the Gram-negative bacterium Acinetobacter baumannii has garnered considerable attention for its unprecedented capacity to rapidly develop resistance to antibacterial therapeutics. This is coupled with the seemingly epidemic emergence of new hyper-virulent strains. Although strain-specific differences for A. baumannii isolates have been well described, these studies have primarily focused on proteinaceous factors. At present, only limited publications have investigated the presence and role of small regulatory RNA (sRNA) transcripts. Herein, we perform such an analysis, describing the RNA-seq-based identification of 78 A. baumannii sRNAs in the AB5075 background. Together with six previously identified elements, we include each of these in a new genome annotation file, which will serve as a tool to investigate regulatory events in this organism. Our work reveals that the sRNAs display high expression, accounting for >50 % of the 20 most strongly expressed genes. Through conservation analysis we identified six classes of similar sRNAs, with one found to be particularly abundant and homologous to regulatory, C4 antisense RNAs found in bacteriophages. These elements appear to be processed from larger transcripts in an analogous manner to the phage C4 molecule and are putatively controlled by two further sRNAs that are strongly antisense to them. Collectively, this study offers a detailed view of the sRNA content of A. baumannii, exposing sequence and structural conservation amongst these elements, and provides novel insight into the potential evolution, and role, of these understudied regulatory molecules.

Use of the Salmonella MgtR peptide as an antagonist of the Mycobacterium MgtC virulence factor

Background: The MgtC virulence factor has been proposed as an attractive target for antivirulence strategies because it is shared by several important bacterial pathogens, including Salmonella enterica and Mycobacterium tuberculosis (Mtb). Aim: A natural antagonistic peptide, MgtR, which interacts with MgtC and modulates its stability, has been identified in Salmonella, and we investigated its efficiency to target MgtC in another pathogen. Materials & methods: We evaluated the interaction between Salmonella MgtR peptide and the Mtb MgtC protein using an in vivo bacterial two-hybrid system and we addressed the effect of exogenously added synthetic MgtR and endogenously expressed peptide. Results: MgtR peptide strongly interacted with Mtb MgtC protein and exogenously added synthetic MgtR peptide-reduced Mtb MgtC level and interfered with the dimerization of Mtb MgtC. Importantly, heterologous expression of MgtR in Mycobacterium bovis BCG resulted in increased phagocytosis and reduced intramacrophage survival. Conclusion: MgtR peptide can target Mtb MgtC protein and reduce mycobacterial macrophage resistance, thus providing a promising new scaffold for the development of antivirulence compounds.

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.

Comparative Life Cycle Transcriptomics Revises Leishmania mexicana Genome Annotation and Links a Chromosome Duplication with Parasitism of Vertebrates

Leishmania spp. are protozoan parasites that have two principal life cycle stages: the motile promastigote forms that live in the alimentary tract of the sandfly and the amastigote forms, which are adapted to survive and replicate in the harsh conditions of the phagolysosome of mammalian macrophages. Here, we used Illumina sequencing of poly-A selected RNA to characterise and compare the transcriptomes of L. mexicana promastigotes, axenic amastigotes and intracellular amastigotes. These data allowed the production of the first transcriptome evidence-based annotation of gene models for this species, including genome-wide mapping of trans-splice sites and poly-A addition sites. The revised genome annotation encompassed 9,169 protein-coding genes including 936 novel genes as well as modifications to previously existing gene models. Comparative analysis of gene expression across promastigote and amastigote forms revealed that 3,832 genes are differentially expressed between promastigotes and intracellular amastigotes. A large proportion of genes that were downregulated during differentiation to amastigotes were associated with the function of the motile flagellum. In contrast, those genes that were upregulated included cell surface proteins, transporters, peptidases and many uncharacterized genes, including 293 of the 936 novel genes. Genome-wide distribution analysis of the differentially expressed genes revealed that the tetraploid chromosome 30 is highly enriched for genes that were upregulated in amastigotes, providing the first evidence of a link between this whole chromosome duplication event and adaptation to the vertebrate host in this group. Peptide evidence for 42 proteins encoded by novel transcripts supports the idea of an as yet uncharacterised set of small proteins in Leishmania spp. with possible implications for host-pathogen interactions.

Complete Genome Sequence of Bifidobacterium longum GT15: Identification and Characterization of Unique and Global Regulatory Genes

In this study, we report the first completely annotated genome sequence of the Russia origin Bifidobacterium longum subsp. longum strain GT15. Comparative genomic analysis of this genome with other available completely annotated genome sequences of B. longum strains isolated from other countries has revealed a high degree of conservation and synteny across the entire genomes. However, it was discovered that the open reading frames to 35 genes were detected only from the B. longum GT15 genome and absent from other genomes B. longum strains (not of Russian origin). These so-called unique genes (UGs) represent a total length of 39,066 bp, with G + C content ranging from 37 to 65 %. Interestingly, certain genes were detected in other B. longum strains of Russian origin. In our analysis, we examined genes for global regulatory systems: proteins of toxin-antitoxin (TA) systems type II, serine/threonine protein kinases (STPKs) of eukaryotic type, and genes of the WhiB-like family proteins. In addition, we have made in silico analysis of all the most significant probiotic genes and considered genes involved in epigenetic regulation and genes responsible for producing various neuromediators. This genome sequence may elucidate the biology of this probiotic strain as a promising candidate for practical (pharmaceutical) applications.

Small membrane proteins - elucidating the function of the needle in the haystack

Membrane proteins are important mediators between the cell and its environment or between different compartments within a cell. However, much less is known about the structure and function of membrane proteins compared to water-soluble proteins. Moreover, until recently a subset of membrane proteins, those shorter than 100 amino acids, have almost completely evaded detection as a result of technical difficulties. These small membrane proteins (SMPs) have been underrepresented in most genomic and proteomic screens of both pro- and eukaryotic cells and, hence, we know much less about their functions in both. Currently, through a combination of bioinformatics, ribosome profiling, and more sensitive proteomics, large numbers of SMPs are being identified and characterized. Herein we describe recent advances in identifying SMPs from genomic and proteomic datasets and describe examples where SMPs have been successfully characterized biochemically. Finally we give an overview of identified functions of SMPs and speculate on the possible roles SMPs play in the cell.

Defining a minimal cell: essentiality of small ORFs and ncRNAs in a genome-reduced bacterium

Identifying all essential genomic components is critical for the assembly of minimal artificial life. In the genome-reduced bacterium Mycoplasma pneumoniae, we found that small ORFs (smORFs; < 100 residues), accounting for 10% of all ORFs, are the most frequently essential genomic components (53%), followed by conventional ORFs (49%). Essentiality of smORFs may be explained by their function as members of protein and/or DNA/RNA complexes. In larger proteins, essentiality applied to individual domains and not entire proteins, a notion we could confirm by expression of truncated domains. The fraction of essential non-coding RNAs (ncRNAs) non-overlapping with essential genes is 5% higher than of non-transcribed regions (0.9%), pointing to the important functions of the former. We found that the minimal essential genome is comprised of 33% (269,410 bp) of the M. pneumoniae genome. Our data highlight an unexpected hidden layer of smORFs with essential functions, as well as non-coding regions, thus changing the focus when aiming to define the minimal essential genome.

Identifying (non-)coding RNAs and small peptides: challenges and opportunities

Over the past decade, high-throughput studies have identified many novel transcripts. While their existence is undisputed, their coding potential and functionality have remained controversial. Recent computational approaches guided by ribosome profiling have indicated that translation is far more pervasive than anticipated and takes place on many transcripts previously assumed to be non-coding. Some of these newly discovered translated transcripts encode short, functional proteins that had been missed in prior screens. Other transcripts are translated, but it might be the process of translation rather than the resulting peptides that serves a function. Here, we review annotation studies in zebrafish to discuss the challenges of placing RNAs onto the continuum that ranges from functional protein-encoding mRNAs to potentially non-functional peptide-producing RNAs to non-coding RNAs. As highlighted by the discovery of the novel signaling peptide Apela/ELABELA/Toddler, accurate annotations can give rise to exciting opportunities to identify the functions of previously uncharacterized transcripts.

Transcriptional analysis of porcine circovirus-like virus P1

Background

Recently identified porcine circovirus-like virus P1 has the smallest DNA viral genome. In this study, we identified the viral genes and their corresponding mRNA transcripts.

Results

The RNAs of P1, synthesized in porcine kidney cells, were examined with northern blotting and PCR analyses.

Eight virus-specific RNAs were detected. Four mRNAs (open reading frames (ORFs) 1, 2, 4, and 5) are encoded by the viral (−) strand and four (ORFs 3, 6, 7, and 8) are encoded by the viral (+) strand. All proteins encoded by the ORFs of the P1 virus are less than 50 amino acids in length, except that encoded by ORF1 (113 amino acids).

Conclusions

We show a very complex viral transcription pattern in P1-infected cells.