A comparison of key aspects of gene regulation in Streptomyces coelicolor and Escherichia coli using nucleotide-resolution transcription maps produced in parallel by global and differential RNA sequencing

Precise tuning of bacterial translation initiation by non-equilibrium 5'-UTR unfolding observed in single mRNAs

Noncoding, structured 5′-untranslated regions (5′-UTRs) of bacterial messenger RNAs (mRNAs) can control translation efficiency by forming structures that either recruit or repel the ribosome. Here we exploit a 5′-UTR embedded preQ₁-sensing, pseudoknotted translational riboswitch to probe how binding of a small ligand controls recruitment of the bacterial ribosome to the partially overlapping Shine-Dalgarno (SD) sequence. Combining single-molecule fluorescence microscopy with mutational analyses, we find that the stability of 30S ribosomal subunit binding is inversely correlated with the free energy needed to unfold the 5′-UTR during mRNA accommodation into the mRNA binding cleft. Ligand binding to the riboswitch stabilizes the structure to both antagonize 30S recruitment and accelerate 30S dissociation. Proximity of the 5′-UTR and stability of the SD:anti-SD interaction both play important roles in modulating the initial 30S-mRNA interaction. Finally, depletion of small ribosomal subunit protein S1, known to help resolve structured 5′-UTRs, further increases the energetic penalty for mRNA accommodation. The resulting model of rapid standby site exploration followed by gated non-equilibrium unfolding of the 5′-UTR during accommodation provides a mechanistic understanding of how translation efficiency is governed by riboswitches and other dynamic structure motifs embedded upstream of the translation initiation site of bacterial mRNAs.

Ms1 RNA Interacts With the RNA Polymerase Core in Streptomyces coelicolor and Was Identified in Majority of Actinobacteria Using a Linguistic Gene Synteny Search

Bacteria employ small non-coding RNAs (sRNAs) to regulate gene expression. Ms1 is an sRNA that binds to the RNA polymerase (RNAP) core and affects the intracellular level of this essential enzyme. Ms1 is structurally related to 6S RNA that binds to a different form of RNAP, the holoenzyme bearing the primary sigma factor. 6S RNAs are widespread in the bacterial kingdom except for the industrially and medicinally important Actinobacteria. While Ms1 RNA was identified in Mycobacterium, it is not clear whether Ms1 RNA is present also in other Actinobacteria species. Here, using a computational search based on secondary structure similarities combined with a linguistic gene synteny approach, we identified Ms1 RNA in Streptomyces. In S. coelicolor, Ms1 RNA overlaps with the previously annotated scr3559 sRNA with an unknown function. We experimentally confirmed that Ms1 RNA/scr3559 associates with the RNAP core without the primary sigma factor HrdB in vivo. Subsequently, we applied the computational approach to other Actinobacteria and identified Ms1 RNA candidates in 824 Actinobacteria species, revealing Ms1 RNA as a widespread class of RNAP binding sRNAs, and demonstrating the ability of our multifactorial computational approach to identify weakly conserved sRNAs in evolutionarily distant genomes.

Pervasive translation in Mycobacterium tuberculosis

Most bacterial ORFs are identified by automated prediction algorithms. However, these algorithms often fail to identify ORFs lacking canonical features such as a length of >50 codons or the presence of an upstream Shine-Dalgarno sequence. Here, we use ribosome profiling approaches to identify actively translated ORFs in Mycobacterium tuberculosis. Most of the ORFs we identify have not been previously described, indicating that the M. tuberculosis transcriptome is pervasively translated. The newly described ORFs are predominantly short, with many encoding proteins of ≤50 amino acids. Codon usage of the newly discovered ORFs suggests that most have not been subject to purifying selection, and hence are unlikely to contribute to cell fitness. Nevertheless, we identify 90 new ORFs (median length of 52 codons) that bear the hallmarks of purifying selection. Thus, our data suggest that pervasive translation of short ORFs in Mycobacterium tuberculosis serves as a rich source for the evolution of new functional proteins.

The context of the ribosome binding site in mRNAs defines specificity of action of kasugamycin, an inhibitor of translation initiation

Several antibiotics targeting the large ribosomal subunit interfere with translation in a context-specific manner, preventing ribosomes from polymerizing specific amino acid sequences. Here, we reveal kasugamycin as a small ribosomal subunit-targeting antibiotic whose action depends on the sequence context of the untranslated messenger RNA (mRNA) segments. We show that kasugamycin-induced ribosomal arrest at the start codons of the genes and the resulting inhibition of gene expression depend on the nature of the mRNA nucleotide immediately preceding the start codon and on the proximity of the stop codon of the upstream cistron. Our findings underlie the importance of mRNA context for the action of protein synthesis inhibitors and might help to guide the development of better antibiotics.

Kasugamycin (KSG) is an aminoglycoside antibiotic widely used in agriculture and exhibits considerable medical potential. Previous studies suggested that KSG interferes with translation by blocking binding of canonical messenger RNA (mRNA) and initiator transfer tRNA (tRNA) to the small ribosomal subunit, thereby preventing initiation of protein synthesis. Here, by using genome-wide approaches, we show that KSG can interfere with translation even after the formation of the 70S initiation complex on mRNA, as the extent of KSG-mediated translation inhibition correlates with increased occupancy of start codons by 70S ribosomes. Even at saturating concentrations, KSG does not completely abolish translation, allowing for continuing expression of some Escherichia coli proteins. Differential action of KSG significantly depends on the nature of the mRNA residue immediately preceding the start codon, with guanine in this position being the most conducive to inhibition by the drug. In addition, the activity of KSG is attenuated by translational coupling as genes whose start codons overlap with the coding regions or the stop codons of the upstream cistrons tend to be less susceptible to drug-mediated inhibition. Altogether, our findings reveal KSG as an example of a small ribosomal subunit-targeting antibiotic with a well-pronounced context specificity of action.

Genome-scale analysis of genetic regulatory elements in Streptomyces avermitilis MA-4680 using transcript boundary information

Background

The gram-positive bacterium, Streptomyces avermitilis, holds industrial importance as the producer of avermectin, a widely used anthelmintic agent, and a heterologous expression host of secondary metabolite-biosynthetic gene clusters. Despite its industrial importance, S. avermitilis’ genome organization and regulation of gene expression remain poorly understood. In this study, four different types of Next-Generation Sequencing techniques, including dRNA-Seq, Term-Seq, RNA-Seq and ribosome profiling, were applied to S. avermitilis to determine transcription units of S. avermitilis at a genome-wide level and elucidate regulatory elements for transcriptional and translational control of individual transcription units.

Result

By applying dRNA-Seq and Term-Seq to S. avermitilis MA-4680, a total of 2361 transcription start sites and 2017 transcript 3′-end positions were identified, respectively, leading to determination of 1601 transcription units encoded in S. avermitilis’ genome. Cataloguing the transcription units and integrated analysis of multiple high-throughput data types revealed the presence of diverse regulatory elements for gene expression, such as promoters, 5′-UTRs, terminators, 3′-UTRs and riboswitches. The conserved promoter motifs were identified from 2361 transcription start sites as 5′-TANNNT and 5′-BTGACN for the − 10 and − 35 elements, respectively. The − 35 element and spacer lengths between − 10 and − 35 elements were critical for transcriptional regulation of functionally distinct genes, suggesting the involvement of unique sigma factors. In addition, regulatory sequences recognized by antibiotic regulatory proteins were identified from the transcription start site information. Analysis of the 3′-end of RNA transcript revealed that stem structure formation is a major determinant for transcription termination of most transcription units.

Conclusions

The transcription unit architecture elucidated from the transcripts’ boundary information provides insights for unique genetic regulatory mechanisms of S. avermitilis. Our findings will elevate S. avermitilis’ potential as a production host for a diverse set of secondary metabolites.

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.

Global Chromosome Topology and the Two-Component Systems in Concerted Manner Regulate Transcription in Streptomyces

Bacterial gene expression is controlled at multiple levels, with chromosome supercoiling being one of the most global regulators. Global DNA supercoiling is maintained by the orchestrated action of topoisomerases. In Streptomyces, mycelial soil bacteria with a complex life cycle, topoisomerase I depletion led to elevated chromosome supercoiling, changed expression of a significant fraction of genes, delayed growth, and blocked sporulation. To identify supercoiling-induced sporulation regulators, we searched for Streptomyces coelicolor transposon mutants that were able to restore sporulation despite high chromosome supercoiling. We established that transposon insertion in genes encoding a novel two-component system named SatKR reversed the sporulation blockage resulting from topoisomerase I depletion. Transposition in satKR abolished the transcriptional induction of the genes within the so-called supercoiling-hypersensitive cluster (SHC). Moreover, we found that activated SatR also induced the same set of SHC genes under normal supercoiling conditions. We determined that the expression of genes in this region impacted S. coelicolor growth and sporulation. Interestingly, among the associated products is another two-component system (SitKR), indicating the potential for cascading regulatory effects driven by the SatKR and SitKR two-component systems. Thus, we demonstrated the concerted activity of chromosome supercoiling and a hierarchical two-component signaling system that impacts gene activity governing Streptomyces growth and sporulation.

IMPORTANCEStreptomyces microbes, soil bacteria with complex life cycle, are the producers of a broad range of biologically active compounds (e.g., antibiotics). Streptomyces bacteria respond to various environmental signals using a complex transcriptional regulation mechanism. Understanding regulation of their gene expression is crucial for Streptomyces application as industrial organisms. Here, on the basis of the results of extensive transcriptomics analyses, we describe the concerted gene regulation by global DNA supercoiling and novel two-component system. Our data indicate that regulated genes encode growth and sporulation regulators. Thus, we demonstrate that Streptomyces bacteria link the global regulatory strategies to adjust life cycle to unfavorable conditions.

A combination of mRNA features influence the efficiency of leaderless mRNA translation initiation

Bacterial translation is thought to initiate by base pairing of the 16S rRNA and the Shine-Dalgarno sequence in the mRNA's 5' untranslated region (UTR). However, transcriptomics has revealed that leaderless mRNAs, which completely lack any 5' UTR, are broadly distributed across bacteria and can initiate translation in the absence of the Shine-Dalgarno sequence. To investigate the mechanism of leaderless mRNA translation initiation, synthetic in vivo translation reporters were designed that systematically tested the effects of start codon accessibility, leader length, and start codon identity on leaderless mRNA translation initiation. Using these data, a simple computational model was built based on the combinatorial relationship of these mRNA features that can accurately classify leaderless mRNAs and predict the translation initiation efficiency of leaderless mRNAs. Thus, start codon accessibility, leader length, and start codon identity combine to define leaderless mRNA translation initiation in bacteria.

Protein synthesis in Mycobacterium tuberculosis as a potential target for therapeutic interventions

Mycobacterium tuberculosis (Mtb) causes one of humankind's deadliest diseases, tuberculosis. Mtb protein synthesis machinery possesses several unique species-specific features, including its ribosome that carries two mycobacterial specific ribosomal proteins, bL37 and bS22, and ribosomal RNA segments. Since the protein synthesis is a vital cellular process that occurs on the ribosome, a detailed knowledge of the structure and function of mycobacterial ribosomes is essential to understand the cell's proteome by translation regulation. Like in many bacterial species such as Bacillus subtilis and Streptomyces coelicolor, two distinct populations of ribosomes have been identified in Mtb. Under low-zinc conditions, Mtb ribosomal proteins S14, S18, L28, and L33 are replaced with their non-zinc binding paralogues. Depending upon the nature of physiological stress, species-specific modulation of translation by stress factors and toxins that interact with the ribosome have been reported. In addition, about one-fourth of messenger RNAs in mycobacteria have been reported to be leaderless, i.e., without 5' UTR regions. However, the mechanism by which they are recruited to the Mtb ribosome is not understood. In this review, we highlight the mycobacteria-specific features of the translation apparatus and propose exploiting these features to improve the efficacy and specificity of existing antibiotics used to treat tuberculosis.

Combining transposon mutagenesis and reporter genes to identify novel regulators of the topA promoter in Streptomyces

BACKGROUND

Identifying the regulatory factors that control transcriptional activity is a major challenge of gene expression studies. Here, we describe the application of a novel approach for in vivo identification of regulatory proteins that may directly or indirectly control the transcription of a promoter of interest in Streptomyces.

RESULTS

A method based on the combination of Tn5 minitransposon-driven random mutagenesis and lux reporter genes was applied for the first time for the Streptomyces genus. As a proof of concept, we studied the topA supercoiling-sensitive promoter, whose activity is dependent on unknown regulatory factors. We found that the sco4804 gene product positively influences topA transcription in S. coelicolor, demonstrating SCO4804 as a novel player in the control of chromosome topology in these bacteria.

CONCLUSIONS

Our approach allows the identification of novel Streptomyces regulators that may be critical for the regulation of gene expression in these antibiotic-producing bacteria.

Systems and synthetic biology to elucidate secondary metabolite biosynthetic gene clusters encoded in Streptomyces genomes

Covering: 2010 to 2020 Over the last few decades, Streptomyces have been extensively investigated for their ability to produce diverse bioactive secondary metabolites. Recent advances in Streptomyces research have been largely supported by improvements in high-throughput technology 'omics'. From genomics, numerous secondary metabolite biosynthetic gene clusters were predicted, increasing their genomic potential for novel bioactive compound discovery. Additional omics, including transcriptomics, translatomics, interactomics, proteomics and metabolomics, have been applied to obtain a system-level understanding spanning entire bioprocesses of Streptomyces, revealing highly interconnected and multi-layered regulatory networks for secondary metabolism. The comprehensive understanding derived from this systematic information accelerates the rational engineering of Streptomyces to enhance secondary metabolite production, integrated with the exploitation of the highly efficient 'Design-Build-Test-Learn' cycle in synthetic biology. In this review, we describe the current status of omics applications in Streptomyces research to better understand the organism and exploit its genetic potential for higher production of valuable secondary metabolites and novel secondary metabolite discovery.

Genome-Scale Transcription-Translation Mapping Reveals Features of Zymomonas mobilis Transcription Units and Promoters

Efforts to rationally engineer synthetic pathways in Zymomonas mobilis are impeded by a lack of knowledge and tools for predictable and quantitative programming of gene regulation at the transcriptional, posttranscriptional, and posttranslational levels. With the detailed functional characterization of the Z. mobilis genome presented in this work, we provide crucial knowledge for the development of synthetic genetic parts tailored to Z. mobilis. This information is vital as researchers continue to develop Z. mobilis for synthetic biology applications. Our methods and statistical analyses also provide ways to rapidly advance the understanding of poorly characterized bacteria via empirical data that enable the experimental validation of sequence-based prediction for genome characterization and annotation.

Zymomonas mobilis is an ethanologenic alphaproteobacterium with promise for the industrial conversion of renewable plant biomass into fuels and chemical bioproducts. Limited functional annotation of the Z. mobilis genome is a current barrier to both fundamental studies of Z. mobilis and its development as a synthetic biology chassis. To gain insight, we collected sample-matched multiomics data, including RNA sequencing (RNA-seq), transcription start site (TSS) sequencing (TSS-seq), termination sequencing (term-seq), ribosome profiling, and label-free shotgun proteomic mass spectrometry, across different growth conditions and used these data to improve annotation and assign functional sites in the Z. mobilis genome. Proteomics and ribosome profiling informed revisions of protein-coding genes, which included 44 start codon changes and 42 added proteins. We developed statistical methods for annotating transcript 5′ and 3′ ends, enabling the identification of 3,940 TSSs and their corresponding promoters and 2,091 transcription termination sites, which were distinguished from RNA processing sites by the lack of an adjacent RNA 5′ end. Our results revealed that Z. mobilis σ^A −35 and −10 promoter elements closely resemble canonical Escherichia coli −35 and −10 elements, with one notable exception: the Z. mobilis −10 element lacks the highly conserved −7 thymine observed in E. coli and other previously characterized σ^A promoters. The σ^A promoters of another alphaproteobacterium, Caulobacter crescentus, similarly lack the conservation of −7 thymine in their −10 elements. Our results anchor the development of Z. mobilis as a platform for synthetic biology and establish strategies for empirical genome annotation that can complement purely computational methods.

IMPORTANCE Efforts to rationally engineer synthetic pathways in Zymomonas mobilis are impeded by a lack of knowledge and tools for predictable and quantitative programming of gene regulation at the transcriptional, posttranscriptional, and posttranslational levels. With the detailed functional characterization of the Z. mobilis genome presented in this work, we provide crucial knowledge for the development of synthetic genetic parts tailored to Z. mobilis. This information is vital as researchers continue to develop Z. mobilis for synthetic biology applications. Our methods and statistical analyses also provide ways to rapidly advance the understanding of poorly characterized bacteria via empirical data that enable the experimental validation of sequence-based prediction for genome characterization and annotation.

Primary transcriptome and translatome analysis determines transcriptional and translational regulatory elements encoded in the Streptomyces clavuligerus genome

Determining transcriptional and translational regulatory elements in GC-rich Streptomyces genomes is essential to elucidating the complex regulatory networks that govern secondary metabolite biosynthetic gene cluster (BGC) expression. However, information about such regulatory elements has been limited for Streptomyces genomes. To address this limitation, a high-quality genome sequence of β-lactam antibiotic-producing Streptomyces clavuligerus ATCC 27 064 is completed, which contains 7163 newly annotated genes. This provides a fundamental reference genome sequence to integrate multiple genome-scale data types, including dRNA-Seq, RNA-Seq and ribosome profiling. Data integration results in the precise determination of 2659 transcription start sites which reveal transcriptional and translational regulatory elements, including −10 and −35 promoter components specific to sigma (σ) factors, and 5′-untranslated region as a determinant for translation efficiency regulation. Particularly, sequence analysis of a wide diversity of the −35 components enables us to predict potential σ-factor regulons, along with various spacer lengths between the −10 and −35 elements. At last, the primary transcriptome landscape of the β-lactam biosynthetic pathway is analyzed, suggesting temporal changes in metabolism for the synthesis of secondary metabolites driven by transcriptional regulation. This comprehensive genetic information provides a versatile genetic resource for rational engineering of secondary metabolite BGCs in Streptomyces.

Multi-level regulation of coelimycin synthesis in Streptomyces coelicolor A3(2)

Despite being a yellow pigment visible to the human eye, coelimycin (CPK) remained to be an undiscovered secondary metabolite for over 50 years of Streptomyces research. Although the function of this polyketide is still unclear, we now know that its "cryptic" nature is attributed to a very complex and precise mechanism of cpk gene cluster regulation in the model actinomycete S. coelicolor A3(2). It responds to the stringent culture density and timing of the transition phase by the quorum-sensing butanolide system and to the specific nutrient availability/uptake signals mediated by the global (pleiotropic) regulators; many of which are two-component signal transduction systems. The final effectors of this regulation cascade are predicted to be two cluster-situated Streptomyces antibiotic regulatory proteins (SARPs) putatively activating the expression of type I polyketide synthase (PKS I) genes. After its synthesis, unstable, colorless antibiotic coelimycin A reacts with specific compounds in the medium losing its antibacterial properties and giving rise to yellow coelimycins P1 and P2. Here we review the current knowledge on coelimycin synthesis regulation in Streptomyces coelicolor A3(2). We focus on the regulatory feedback loop which interconnects the butanolide system with other cpk cluster-situated regulators. We also present the effects exerted on cpk genes expression by the global, pleiotropic regulators, and the regulatory connections between cpk and other biosynthetic gene clusters.

Regulation of antibiotic production in Actinobacteria: new perspectives from the post-genomic era

Covering: 2000 to 2018 The antimicrobial activity of many of their natural products has brought prominence to the Streptomycetaceae, a family of Gram-positive bacteria that inhabit both soil and aquatic sediments. In the natural environment, antimicrobial compounds are likely to limit the growth of competitors, thereby offering a selective advantage to the producer, in particular when nutrients become limited and the developmental programme leading to spores commences. The study of the control of this secondary metabolism continues to offer insights into its integration with a complex lifecycle that takes multiple cues from the environment and primary metabolism. Such information can then be harnessed to devise laboratory screening conditions to discover compounds with new or improved clinical value. Here we provide an update of the review we published in NPR in 2011. Besides providing the essential background, we focus on recent developments in our understanding of the underlying regulatory networks, ecological triggers of natural product biosynthesis, contributions from comparative genomics and approaches to awaken the biosynthesis of otherwise silent or cryptic natural products. In addition, we highlight recent discoveries on the control of antibiotic production in other Actinobacteria, which have gained considerable attention since the start of the genomics revolution. New technologies that have the potential to produce a step change in our understanding of the regulation of secondary metabolism are also described.

Transposable element insertions shape gene regulation and melanin production in a fungal pathogen of wheat

BACKGROUND

Fungal plant pathogens pose major threats to crop yield and sustainable food production if they are highly adapted to their host and the local environment. Variation in gene expression contributes to phenotypic diversity within fungal species and affects adaptation. However, very few cases of adaptive regulatory changes have been reported in fungi and the underlying mechanisms remain largely unexplored. Fungal pathogen genomes are highly plastic and harbor numerous insertions of transposable elements, which can potentially contribute to gene expression regulation. In this work, we elucidated how transposable elements contribute to variation in melanin accumulation, a quantitative trait in fungi that affects survival under stressful conditions.

RESULTS

We demonstrated that differential transcriptional regulation of the gene encoding the transcription factor Zmr1, which controls expression of the genes in the melanin biosynthetic gene cluster, is responsible for variation in melanin accumulation in the fungal plant pathogen Zymoseptoria tritici. We show that differences in melanin levels between two strains of Z. tritici are due to two levels of transcriptional regulation: (1) variation in the promoter sequence of Zmr1 and (2) an insertion of transposable elements upstream of the Zmr1 promoter. Remarkably, independent insertions of transposable elements upstream of Zmr1 occurred in 9% of Z. tritici strains from around the world and negatively regulated Zmr1 expression, contributing to variation in melanin accumulation.

CONCLUSIONS

Our studies identified two levels of transcriptional control that regulate the synthesis of melanin. We propose that these regulatory mechanisms evolved to balance the fitness costs associated with melanin production against its positive contribution to survival in stressful environments.

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.

Sporulation-specific cell division defects in ylmE mutants of Streptomyces coelicolor are rescued by additional deletion of ylmD

Cell division during the reproductive phase of the Streptomyces life-cycle requires tight coordination between synchronous formation of multiple septa and DNA segregation. One remarkable difference with most other bacterial systems is that cell division in Streptomyces is positively controlled by the recruitment of FtsZ by SsgB. Here we show that deletion of ylmD (SCO2081) or ylmE (SCO2080), which lie in operon with ftsZ in the dcw cluster of actinomycetes, has major consequences for sporulation-specific cell division in Streptomyces coelicolor. Electron and fluorescence microscopy demonstrated that ylmE mutants have a highly aberrant phenotype with defective septum synthesis, and produce very few spores with low viability and high heat sensitivity. FtsZ-ring formation was also highly disturbed in ylmE mutants. Deletion of ylmD had a far less severe effect on sporulation. Interestingly, the additional deletion of ylmD restored sporulation to the ylmE null mutant. YlmD and YlmE are not part of the divisome, but instead localize diffusely in aerial hyphae, with differential intensity throughout the sporogenic part of the hyphae. Taken together, our work reveals a function for YlmD and YlmE in the control of sporulation-specific cell division in S. coelicolor, whereby the presence of YlmD alone results in major developmental defects.

Non-canonical Escherichia coli transcripts lacking a Shine-Dalgarno motif have very different translational efficiencies and do not form a coherent group

Translation initiation in 50-70 % of transcripts in Escherichia coli requires base pairing between the Shine-Dalgarno (SD) motif in the mRNA and the anti-SD motif at the 3' end of the 16S rRNA. However, 30-50 % of E. coli transcripts are non-canonical and are not preceded by an SD motif. The 5' ends of 44 E. coli transcripts were determined, all of which contained a 5'-UTR (no leaderless transcripts), but only a minority contained an SD motif. The 5'-UTR lengths were compared with those listed in RegulonDB and reported in previous publications, and the identities and differences were obtained in all possible combinations. We aimed to quantify the translational efficiencies of non-canonical 5'-UTRs using GusA reporter gene assays and Northern blot analyses. Ten non-canonical 5'-UTRs and two control 5'-UTRs with an SD motif were cloned upstream of the gusA gene. The translational efficiencies were quantified under five different conditions (different growth rates via two different temperatures and two different carbon sources, and heat shock). The translational efficiencies of the non-canonical 5'-UTRs varied widely, from 5 to 384 % of the positive control. In addition, the non-canonical transcripts did not exhibit a common regulatory pattern with changing environmental parameters. No correlation could be observed between the translational efficiencies of the non-canonical 5'-UTRs and their lengths, sequences, GC content, or predicted secondary structures. The introduction of an SD motif enhanced the translational efficiency of a poorly translated non-canonical transcript, while the efficiency of a well-translated non-canonical transcript remained unchanged. Taken together, the mechanisms of translation initiation at non-canonical transcripts in E. coli still need to be elucidated.

Genome-wide mRNA processing in methanogenic archaea reveals post-transcriptional regulation of ribosomal protein synthesis

Unlike stable RNAs that require processing for maturation, prokaryotic cellular mRNAs generally follow an 'all-or-none' pattern. Herein, we used a 5΄ monophosphate transcript sequencing (5΄P-seq) that specifically captured the 5΄-end of processed transcripts and mapped the genome-wide RNA processing sites (PSSs) in a methanogenic archaeon. Following statistical analysis and stringent filtration, we identified 1429 PSSs, among which 23.5% and 5.4% were located in 5΄ untranslated region (uPSS) and intergenic region (iPSS), respectively. A predominant uridine downstream PSSs served as a processing signature. Remarkably, 5΄P-seq detected overrepresented uPSS and iPSS in the polycistronic operons encoding ribosomal proteins, and the majority upstream and proximal ribosome binding sites, suggesting a regulatory role of processing on translation initiation. The processed transcripts showed increased stability and translation efficiency. Particularly, processing within the tricistronic transcript of rplA-rplJ-rplL enhanced the translation of rplL, which can provide a driving force for the 1:4 stoichiometry of L10 to L12 in the ribosome. Growth-associated mRNA processing intensities were also correlated with the cellular ribosomal protein levels, thereby suggesting that mRNA processing is involved in tuning growth-dependent ribosome synthesis. In conclusion, our findings suggest that mRNA processing-mediated post-transcriptional regulation is a potential mechanism of ribosomal protein synthesis and stoichiometry.

Global Mapping of Small RNA-Target Interactions in Bacteria

Small RNAs (sRNAs) associated with the RNA chaperon protein Hfq are key posttranscriptional regulators of gene expression in bacteria. Deciphering the sRNA-target interactome is an essential step toward understanding the roles of sRNAs in the cellular networks. We developed a broadly applicable methodology termed RIL-seq (RNA interaction by ligation and sequencing), which integrates experimental and computational tools for in vivo transcriptome-wide identification of interactions involving Hfq-associated sRNAs. By applying this methodology to Escherichia coli we discovered an extensive network of interactions involving RNA pairs showing sequence complementarity. We expand the ensemble of targets for known sRNAs, uncover additional Hfq-bound sRNAs encoded in various genomic regions along with their trans encoded targets, and provide insights into binding and possible cycling of RNAs on Hfq. Comparison of the sRNA interactome under various conditions has revealed changes in the sRNA repertoire as well as substantial re-wiring of the network between conditions.

•

A widely applicable method for in vivo global mapping of small RNA interactome

•

Substantial re-wiring of the network upon changes in cellular conditions

•

Regulatory circuits involving two regulators derived from the same transcript

•

sRNAs acting in trans are encoded within almost every possible genomic element

dRNA-seq transcriptional profiling of the FK506 biosynthetic gene cluster in Streptomyces tsukubaensis NRRL18488 and general analysis of the transcriptome

FK506 (tacrolimus) is a valuable immunosuppressant produced by several Streptomyces strains. In the genome of the wild type producer Streptomyces tsukubaensis NRRL18488, FK506 biosynthesis is encoded by a gene cluster that spans 83.5 (kb). A whole transcriptome differential shotgun sequencing (dRNA-seq) of S. tsukubaensis was performed to analyze transcription at 2 different time points; before and during active FK506 production. In total, 8,914 transcription start sites were identified in either condition, which enabled precise determination of the 5'-UTR length of the corresponding transcripts as well as the identification of 2 consensus sequence motifs in the promoter regions. The transcription start sites of all gene operons within the FK506 cluster were identified, including 3 examples of leaderless RNA transcripts. These data provide detailed insight into the transcription of the FK506 biosynthetic gene cluster to support future regulatory studies, genetic manipulation, and industrial production.

Landscape of RNA polyadenylation in E. coli

Polyadenylation is thought to be involved in the degradation and quality control of bacterial RNAs but relatively few examples have been investigated. We used a combination of 5΄-tagRACE and RNA-seq to analyze the total RNA content from a wild-type strain and from a poly(A)polymerase deleted mutant. A total of 178 transcripts were either up- or down-regulated in the mutant when compared to the wild-type strain. Poly(A)polymerase up-regulates the expression of all genes related to the FliA regulon and several previously unknown transcripts, including numerous transporters. Notable down-regulation of genes in the expression of antigen 43 and components of the type 1 fimbriae was detected. The major consequence of the absence of poly(A)polymerase was the accumulation of numerous sRNAs, antisense transcripts, REP sequences and RNA fragments resulting from the processing of entire transcripts. A new algorithm to analyze the position and composition of post-transcriptional modifications based on the sequence of unencoded 3΄-ends, was developed to identify polyadenylated molecules. Overall our results shed new light on the broad spectrum of action of polyadenylation on gene expression and demonstrate the importance of poly(A) dependent degradation to remove structured RNA fragments.

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.

Toxins MazF and MqsR cleave Escherichia coli rRNA precursors at multiple sites

The endoribonuclease toxins of the E. coli toxin-antitoxin systems arrest bacterial growth and protein synthesis by targeting cellular mRNAs. As an exception, E. coli MazF was reported to cleave also 16S rRNA at a single site and separate an anti-Shine-Dalgarno sequence-containing RNA fragment from the ribosome. We noticed extensive rRNA fragmentation in response to induction of the toxins MazF and MqsR, which suggested that these toxins can cleave rRNA at multiple sites. We adapted differential RNA-sequencing to map the toxin-cleaved 5'- and 3'-ends. Our results show that the MazF and MqsR cleavage sites are located within structured rRNA regions and, therefore, are not accessible in assembled ribosomes. Most of the rRNA fragments are located in the aberrant ribosomal subunits that accumulate in response to toxin induction and contain unprocessed rRNA precursors. We did not detect MazF- or MqsR-cleaved rRNA in stationary phase bacteria and in assembled ribosomes. Thus, we conclude that MazF and MqsR cleave rRNA precursors before the ribosomes are assembled and potentially facilitate the decay of surplus rRNA transcripts during stress.

The Coordinated Positive Regulation of Topoisomerase Genes Maintains Topological Homeostasis in Streptomyces coelicolor

Maintaining an optimal level of chromosomal supercoiling is critical for the progression of DNA replication and transcription. Moreover, changes in global supercoiling affect the expression of a large number of genes and play a fundamental role in adapting to stress. Topoisomerase I (TopA) and gyrase are key players in the regulation of bacterial chromosomal topology through their respective abilities to relax and compact DNA. Soil bacteria such as Streptomyces species, which grow as branched, multigenomic hyphae, are subject to environmental stresses that are associated with changes in chromosomal topology. The topological fluctuations modulate the transcriptional activity of a large number of genes and in Streptomyces are related to the production of antibiotics. To better understand the regulation of topological homeostasis in Streptomyces coelicolor, we investigated the interplay between the activities of the topoisomerase-encoding genes topA and gyrBA. We show that the expression of both genes is supercoiling sensitive. Remarkably, increased chromosomal supercoiling induces the topA promoter but only slightly influences gyrBA transcription, while DNA relaxation affects the topA promoter only marginally but strongly activates the gyrBA operon. Moreover, we showed that exposure to elevated temperatures induces rapid relaxation, which results in changes in the levels of both topoisomerases. We therefore propose a unique mechanism of S. coelicolor chromosomal topology maintenance based on the supercoiling-dependent stimulation, rather than repression, of the transcription of both topoisomerase genes. These findings provide important insight into the maintenance of topological homeostasis in an industrially important antibiotic producer.

IMPORTANCE We describe the unique regulation of genes encoding two topoisomerases, topoisomerase I (TopA) and gyrase, in a model Streptomyces species. Our studies demonstrate the coordination of topoisomerase gene regulation, which is crucial for maintenance of topological homeostasis. Streptomyces species are producers of a plethora of biologically active secondary metabolites, including antibiotics, antitumor agents, and immunosuppressants. The significant regulatory factor controlling the secondary metabolism is the global chromosomal topology. Thus, the investigation of chromosomal topology homeostasis in Streptomyces strains is crucial for their use in industrial applications as producers of secondary metabolites.

Cosmid based mutagenesis causes genetic instability in Streptomyces coelicolor, as shown by targeting of the lipoprotein signal peptidase gene

Bacterial lipoproteins are extracellular proteins tethered to cell membranes by covalently attached lipids. Deleting the lipoprotein signal peptidase (lsp) gene in Streptomyces coelicolor results in growth and developmental defects that cannot be restored by reintroducing lsp. This led us to hypothesise that lsp is essential and that the lsp mutant we isolated previously had acquired compensatory secondary mutations. Here we report resequencing of the genomes of wild-type M145 and the cis-complemented ∆lsp mutant (BJT1004) to map and identify these secondary mutations but we show that they do not increase the efficiency of disrupting lsp and are not lsp suppressors. We provide evidence that they are induced by introducing the cosmid St4A10∆lsp, as part of ReDirect PCR mutagenesis protocol, which transiently duplicates a number of important cell division genes. Disruption of lsp using a suicide vector (which does not result in gene duplication) still results in growth and developmental delays and we conclude that loss of Lsp function results in developmental defects due to the loss of all lipoproteins from the cell membrane. Significantly, our results also indicate the use of cosmid libraries for the genetic manipulation of bacteria can lead to phenotypes not necessarily linked to the gene(s) of interest.

The dynamic transcriptional and translational landscape of the model antibiotic producer Streptomyces coelicolor A3(2)

Individual Streptomyces species have the genetic potential to produce a diverse array of natural products of commercial, medical and veterinary interest. However, these products are often not detectable under laboratory culture conditions. To harness their full biosynthetic potential, it is important to develop a detailed understanding of the regulatory networks that orchestrate their metabolism. Here we integrate nucleotide resolution genome-scale measurements of the transcriptome and translatome of Streptomyces coelicolor, the model antibiotic-producing actinomycete. Our systematic study determines 3,570 transcription start sites and identifies 230 small RNAs and a considerable proportion (∼21%) of leaderless mRNAs; this enables deduction of genome-wide promoter architecture. Ribosome profiling reveals that the translation efficiency of secondary metabolic genes is negatively correlated with transcription and that several key antibiotic regulatory genes are translationally induced at transition growth phase. These findings might facilitate the design of new approaches to antibiotic discovery and development.

Insights into the Stress Response Triggered by Kasugamycin in Escherichia coli

The bacteriostatic aminoglycoside antibiotic kasugamycin inhibits protein synthesis at an initial step without affecting translation elongation. It binds to the mRNA track of the ribosome and prevents formation of the translation initiation complex on canonical mRNAs. In contrast, translation of leaderless mRNAs continues in the presence of the drug in vivo. Previously, we have shown that kasugamycin treatment in E. coli stimulates the formation of protein-depleted ribosomes that are selective for leaderless mRNAs. Here, we provide evidence that prolonged kasugamycin treatment leads to selective synthesis of specific proteins. Our studies indicate that leaderless and short-leadered mRNAs are generated by different molecular mechanisms including alternative transcription and RNA processing. Moreover, we provide evidence for ribosome heterogeneity in response to kasugamycin treatment by alteration of the modification status of the stalk proteins bL7/L12.

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.

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

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.

Transcriptional analysis of the cell division-related ssg genes in Streptomyces coelicolor reveals direct control of ssgR by AtrA

SsgA-like proteins are a family of actinomycete-specific regulatory proteins that control cell division and spore maturation in streptomycetes. SsgA and SsgB together activate sporulation-specific cell division by controlling the localization of FtsZ. Here we report the identification of novel regulators that control the transcription of the ssgA-like genes. Transcriptional regulators controlling ssg gene expression were identified using a DNA-affinity capture assay. Supporting transcriptional and DNA binding studies showed that the ssgA activator gene ssgR is controlled by the TetR-family regulator AtrA, while the γ-butyrolactone-responsive AdpA (SCO2792) and SlbR (SCO0608) and the metabolic regulator Rok7B7 (SCO6008) were identified as candidate regulators for the cell division genes ssgA, ssgB and ssgG. Transcription of the cell division gene ssgB depended on the sporulation genes whiA and whiH, while ssgR, ssgA and ssgD were transcribed independently of the whi genes. Our work sheds new light on the mechanisms by which sporulation-specific cell division is controlled in Streptomyces.

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

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

A MarR Family Transcriptional Regulator, DptR3, Activates Daptomycin Biosynthesis and Morphological Differentiation in Streptomyces roseosporus

Daptomycin produced by Streptomyces roseosporus is an important lipopeptide antibiotic used to treat human infections caused by Gram-positive pathogenic bacteria, including drug-resistant strains. The genetic basis for regulatory mechanisms of daptomycin production is poorly known. Here, we characterized the dptR3 gene, which encodes a MarR family transcriptional regulator located adjacent to the known daptomycin biosynthetic (dpt) genes. Deletion of dptR3 reduced daptomycin production significantly and delayed aerial mycelium formation and sporulation on solid media. Dissection of the mechanism underlying the function of DptR3 in daptomycin production revealed that it stimulates daptomycin production indirectly by altering the transcription of dpt structural genes. DptR3 directly activated the transcription of its own gene, dptR3, but repressed the transcription of the adjacent, divergent gene orf16 (which encodes a putative ABC transporter ATP-binding protein). A 66-nucleotide DptR3-binding site in the intergenic region of dptR3-orf16 was determined by DNase I footprinting, and the palindromic sequence TCATTGTTACCTATGCTCACAATGA (underlining indicates inverted repeats) in the protected region was found to be essential for DptR3 binding. orf16, the major target gene of DptR3, exerted a positive effect on daptomycin biosynthesis. Our findings indicate that DptR3 functions as a global regulator that positively controls daptomycin production and morphological development in S. roseosporus.