Variations in the non-coding transcriptome as a driver of inter-strain divergence and physiological adaptation in bacteria

Adaptation of the Mycobacterium tuberculosis transcriptome to biofilm growth

Mycobacterium tuberculosis (M. tb), the causative agent of tuberculosis (TB), is a leading global cause of death from infectious disease. Biofilms are increasingly recognized as a relevant growth form during M. tb infection and may impede treatment by enabling bacterial drug and immune tolerance. M. tb has a complicated regulatory network that has been well-characterized for many relevant disease states, including dormancy and hypoxia. However, despite its importance, our knowledge of the genes and pathways involved in biofilm formation is limited. Here we characterize the biofilm transcriptomes of fully virulent clinical isolates and find that the regulatory systems underlying biofilm growth vary widely between strains and are also distinct from regulatory programs associated with other environmental cues. We used experimental evolution to investigate changes to the transcriptome during adaptation to biofilm growth and found that the application of a uniform selection pressure resulted in loss of strain-to-strain variation in gene expression, resulting in a more uniform biofilm transcriptome. The adaptive trajectories of transcriptomes were shaped by the genetic background of the M. tb population leading to convergence on a sub-lineage specific transcriptome. We identified widespread upregulation of non-coding RNA (ncRNA) as a common feature of the biofilm transcriptome and hypothesize that ncRNA function in genome-wide modulation of gene expression, thereby facilitating rapid regulatory responses to new environments. These results reveal a new facet of the M. tb regulatory system and provide valuable insight into how M. tb adapts to new environments.

Protein NirP1 regulates nitrite reductase and nitrite excretion in cyanobacteria

When the supply of inorganic carbon is limiting, photosynthetic cyanobacteria excrete nitrite, a toxic intermediate in the ammonia assimilation pathway from nitrate. It has been hypothesized that the excreted nitrite represents excess nitrogen that cannot be further assimilated due to the missing carbon, but the underlying molecular mechanisms are unclear. Here, we identified a protein that interacts with nitrite reductase, regulates nitrogen metabolism and promotes nitrite excretion. The protein, which we named NirP1, is encoded by an unannotated gene that is upregulated under low carbon conditions and controlled by transcription factor NtcA, a central regulator of nitrogen homeostasis. Ectopic overexpression of nirP1 in Synechocystis sp. PCC 6803 resulted in a chlorotic phenotype, delayed growth, severe changes in amino acid pools, and nitrite excretion. Coimmunoprecipitation experiments indicated that NirP1 interacts with nitrite reductase, a central enzyme in the assimilation of ammonia from nitrate/nitrite. Our results reveal that NirP1 is widely conserved in cyanobacteria and plays a crucial role in the coordination of C/N primary metabolism by targeting nitrite reductase.

Adaptation of the Mycobacterium tuberculosis transcriptome to biofilm growth

Mycobacterium tuberculosis ( M. tb ), the causative agent of tuberculosis (TB), is a leading global cause of death from infectious disease. Biofilms are increasingly recognized as a relevant growth form during M. tb infection and may impede treatment by enabling bacterial drug and immune tolerance. M. tb has a complicated regulatory network that has been well-characterized for many relevant disease states, including dormancy and hypoxia. However, despite its importance, our knowledge of the genes and pathways involved in biofilm formation is limited. Here we characterize the biofilm transcriptomes of fully virulent clinical isolates and find that the regulatory systems underlying biofilm growth vary widely between strains and are also distinct from regulatory programs associated with other environmental cues. We used experimental evolution to investigate changes to the transcriptome during adaptation to biofilm growth and found that the application of a uniform selection pressure resulted in loss of strain-to-strain variation in gene expression, resulting in a more uniform biofilm transcriptome. The adaptive trajectories of transcriptomes were shaped by the genetic background of the M. tb population leading to convergence on a sub-lineage specific transcriptome. We identified widespread upregulation of non-coding RNA (ncRNA) as a common feature of the biofilm transcriptome and hypothesize that ncRNA function in genome-wide modulation of gene expression, thereby facilitating rapid regulatory responses to new environments. These results reveal a new facet of the M. tb regulatory system and provide valuable insight into how M. tb adapts to new environments.

Importance

Understanding mechanisms of resistance and tolerance in Mycobacterium tuberculosis ( M. tb ) can help us develop new treatments that capitalize on M. tb 's vulnerabilities. Here we used transcriptomics to study both the regulation of biofilm formation in clinical isolates as well as how those regulatory systems adapt to new environments. We find that closely related clinical populations have diverse strategies for growth under biofilm conditions, and that genetic background plays a large role in determining the trajectory of evolution. These results have implications for future treatment strategies that may be informed by our knowledge of the evolutionary constraints on strain(s) from an individual infection. This work provides new information about the mechanisms of biofilm formation in M. tb and outlines a framework for population level approaches for studying bacterial adaptation.

Regulation of RNase E during the UV stress response in the cyanobacterium Synechocystis sp. PCC 6803

Endoribonucleases govern the maturation and degradation of RNA and are indispensable in the posttranscriptional regulation of gene expression. A key endoribonuclease in Gram-negative bacteria is RNase E. To ensure an appropriate supply of RNase E, some bacteria, such as Escherichia coli, feedback-regulate RNase E expression via the rne 5'-untranslated region (5' UTR) in cis. However, the mechanisms involved in the control of RNase E in other bacteria largely remain unknown. Cyanobacteria rely on solar light as an energy source for photosynthesis, despite the inherent ultraviolet (UV) irradiation. In this study, we first investigated globally the changes in gene expression in the cyanobacterium Synechocystis sp. PCC 6803 after a brief exposure to UV. Among the 407 responding genes 2 h after UV exposure was a prominent upregulation of rne mRNA level. Moreover, the enzymatic activity of RNase E rapidly increased as well, although the protein stability decreased. This unique response was underpinned by the increased accumulation of full-length rne mRNA caused by the stabilization of its 5' UTR and suppression of premature transcriptional termination, but not by an increased transcription rate. Mapping of RNA 3' ends and in vitro cleavage assays revealed that RNase E cleaves within a stretch of six consecutive uridine residues within the rne 5' UTR, indicating autoregulation. These observations suggest that RNase E in cyanobacteria contributes to reshaping the transcriptome during the UV stress response and that its required activity level is secured at the RNA level despite the enhanced turnover of the protein.

Genome-wide identification and characterization of Fur-binding sites in the cyanobacteria Synechocystis sp. PCC 6803 and PCC 6714

The Ferric uptake regulator (Fur) is crucial to both pathogenic and non-pathogenic bacteria for the maintenance of iron homeostasis as well as the defence against reactive oxygen species. Based on datasets from the genome-wide mapping of transcriptional start sites and transcriptome data, we identified a high confidence regulon controlled by Fur for the model cyanobacterium Synechocystis sp. PCC 6803 and its close relative, strain 6714, based on the conserved strong iron starvation response and Fur-binding site occurrence. This regulon comprises 33 protein-coding genes and the sRNA IsaR1 that are under the control of 16 or 14 individual promoters in strains 6803 and 6714, respectively. The associated gene functions are mostly restricted to transporters and enzymes involved in the uptake and storage of iron ions, with few exceptions or unknown functional relevance. Within the isiABC operon, we identified a previously neglected gene encoding a small cysteine-rich protein, which we suggest calling, IsiE. The regulation of iron uptake, storage, and utilization ultimately results from the interplay between the Fur regulon, several other transcription factors, the FtsH3 protease, and the sRNA IsaR1.

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.

Construction of prokaryotic strand-specific primary-transcripts saturated RNASeq library by controlled heat magnesium-dependent mRNA degradation

The main limiting factors for RNA-Seq analysis are quality and quantity of the isolated mRNA. In prokaryotes, the proportion of messenger RNA to total RNA is rather low. Therefore, the main strategy of library preparation for sequencing is mRNA enrichment. Ribosomal and transfer RNAs, both monophosphorylated at the 5'-ends, are the major fractions of total RNA, while the bulk of primary transcripts is triphosphorylated at the 5'-teminus. Due to its low molecular weight, transfer RNA could be easily removed by a quick precipitation in LiCl solution. Ribosomal RNA may be degraded enzymatically by 5'-end terminal exonuclease XRN-1. These steps allow enriching samples in mRNA during the first stages of RNA-Seq library preparation. The desired level of fragmentation of enriched mRNA necessary for the 2^nd generation sequencing can be controlled by the duration of incubation at elevated temperatures in the presence of Mg²⁺-ions. Here, we describe a simple protocol for construction of the primary prokaryotic mRNA-saturated library without long depletion procedures.

Regulatory RNA at the crossroads of carbon and nitrogen metabolism in photosynthetic cyanobacteria

Cyanobacteria are photosynthetic bacteria that populate widely different habitats. Accordingly, cyanobacteria exhibit a wide spectrum of lifestyles, physiologies, and morphologies and possess genome sizes and gene numbers which may vary by up to a factor of ten within the phylum. Consequently, large differences exist between individual species in the size and complexity of their regulatory networks. Several non-coding RNAs have been identified that play crucial roles in the acclimation responses of cyanobacteria to changes in the environment. Some of these regulatory RNAs are conserved throughout the cyanobacterial phylum, while others exist only in a few taxa. Here we give an overview on characterized regulatory RNAs in cyanobacteria, with a focus on regulators of photosynthesis, carbon and nitrogen metabolism. However, chances are high that these regulators represent just the tip of the iceberg.

The Primary Antisense Transcriptome of Halobacterium salinarum NRC-1

Antisense RNAs (asRNAs) are present in diverse organisms and play important roles in gene regulation. In this work, we mapped the primary antisense transcriptome in the halophilic archaeon Halobacterium salinarum NRC-1. By reanalyzing publicly available data, we mapped antisense transcription start sites (aTSSs) and inferred the probable 3′ ends of these transcripts. We analyzed the resulting asRNAs according to the size, location, function of genes on the opposite strand, expression levels and conservation. We show that at least 21% of the genes contain asRNAs in H. salinarum. Most of these asRNAs are expressed at low levels. They are located antisense to genes related to distinctive characteristics of H. salinarum, such as bacteriorhodopsin, gas vesicles, transposases and other important biological processes such as translation. We provide evidence to support asRNAs in type II toxin–antitoxin systems in archaea. We also analyzed public Ribosome profiling (Ribo-seq) data and found that ~10% of the asRNAs are ribosome-associated non-coding RNAs (rancRNAs), with asRNAs from transposases overrepresented. Using a comparative transcriptomics approach, we found that ~19% of the asRNAs annotated in H. salinarum belong to genes with an ortholog in Haloferax volcanii, in which an aTSS could be identified with positional equivalence. This shows that most asRNAs are not conserved between these halophilic archaea.

Small RNAs Involved in Regulation of Nitrogen Metabolism

Global (metabolic) regulatory networks allow microorganisms to survive periods of nitrogen starvation or general nutrient stress. Uptake and utilization of various nitrogen sources are thus commonly tightly regulated in Prokarya (Bacteria and Archaea) in response to available nitrogen sources. Those well-studied regulations occur mainly at the transcriptional and posttranslational level. Surprisingly, and in contrast to their involvement in most other stress responses, small RNAs (sRNAs) involved in the response to environmental nitrogen fluctuations are only rarely reported. In addition to sRNAs indirectly affecting nitrogen metabolism, only recently it was demonstrated that three sRNAs were directly involved in regulation of nitrogen metabolism in response to changes in available nitrogen sources. All three trans-acting sRNAs are under direct transcriptional control of global nitrogen regulators and affect expression of components of nitrogen metabolism (glutamine synthetase, nitrogenase, and PII-like proteins) by either masking the ribosome binding site and thus inhibiting translation initiation or stabilizing the respective target mRNAs. Most likely, there are many more sRNAs and other types of noncoding RNAs, e.g., riboswitches, involved in the regulation of nitrogen metabolism in Prokarya that remain to be uncovered. The present review summarizes the current knowledge on sRNAs involved in nitrogen metabolism and their biological functions and targets.

The Distinctive Regulation of Cyanobacterial Glutamine Synthetase

Glutamine synthetase (GS) features prominently in bacterial nitrogen assimilation as it catalyzes the entry of bioavailable nitrogen in form of ammonium into cellular metabolism. The classic example, the comprehensively characterized GS of enterobacteria, is subject to exquisite regulation at multiple levels, among them gene expression regulation to control GS abundance, as well as feedback inhibition and covalent modifications to control enzyme activity. Intriguingly, the GS of the ecologically important clade of cyanobacteria features fundamentally different regulatory systems to those of most prokaryotes. These include the interaction with small proteins, the so-called inactivating factors (IFs) that inhibit GS linearly with their abundance. In addition to this protein interaction-based regulation of GS activity, cyanobacteria use alternative elements to control the synthesis of GS and IFs at the transcriptional level. Moreover, cyanobacteria evolved unique RNA-based regulatory mechanisms such as glutamine riboswitches to tightly tune IF abundance. In this review, we aim to outline the current knowledge on the distinctive features of the cyanobacterial GS encompassing the overall control of its activity, sensing the nitrogen status, transcriptional and post-transcriptional regulation, as well as strain-specific differences.

Widespread Antisense Transcription in Prokaryotes

Although bacterial genomes are usually densely protein-coding, genome-wide mapping approaches of transcriptional start sites revealed that a significant fraction of the identified promoters drive the transcription of noncoding RNAs. These can be trans -acting RNAs, mainly originating from intergenic regions and, in many studied examples, possessing regulatory functions. However, a significant fraction of these noncoding RNAs consist of natural antisense transcripts (asRNAs), which overlap other transcriptional units. Naturally occurring asRNAs were first observed to play a role in bacterial plasmid replication and in bacteriophage λ more than 30 years ago. Today’s view is that asRNAs abound in all three domains of life. There are several examples of asRNAs in bacteria with clearly defined functions. Nevertheless, many asRNAs appear to result from pervasive initiation of transcription, and some data point toward global functions of such widespread transcriptional activity, explaining why the search for a specific regulatory role is sometimes futile. In this review, we give an overview about the occurrence of antisense transcription in bacteria, highlight particular examples of functionally characterized asRNAs, and discuss recent evidence pointing at global relevance in RNA processing and transcription-coupled DNA repair.

Ethylene causes transcriptomic changes in Synechocystis during phototaxis

Ethylene is well known as a plant hormone, but its role in bacteria is poorly studied. We recently showed that Synechocystis sp. Strain PCC 6803 has a functional receptor for ethylene, ethylene response 1 (Etr1), that is involved in various processes such as phototaxis in response to directional light and biofilm formation. Here, we use RNA sequencing to examine the changes in gene transcripts caused by ethylene under phototaxis conditions. Over 500 gene transcripts across many functional categories, of approximately 3700 protein-encoding genes, were altered by application of ethylene. In general, ethylene caused both up- and downregulation of genes within a functional category. However, the transcript levels of amino acid metabolism genes were mainly upregulated and cell envelope genes were mostly downregulated by ethylene. The changes in cell envelope genes correlate with our prior observation that ethylene affects cell surface properties to alter cell motility. Ethylene caused a twofold or more change in 62 transcripts with the largest category of upregulated genes annotated as transporters and the largest category of downregulated genes annotated as glycosyltransferases which sometimes are involved in changing the composition of sugars on the cell surface. Consistent with changes in cell envelope, glycosyltransferase, and transporter gene transcripts, application of ethylene altered the levels of specific sugar moieties on the surface of cells. Light signaling from Etr1 involves two proteins (Slr1213 and Slr1214) and a small, noncoding RNA, carbon stress-induced RNA1 (csiR1). Application of ethylene caused a rapid, but transient, decrease in the transcript levels of etr1, slr1213, and slr1214 and a rapid and prolonged decrease in csiR1 transcript. Deletion of Slr1214 caused a large increase in csiR1 transcript levels and ethylene lowered csiR1 transcript. These data combined with prior reports indicate that ethylene functions as a signal to affect a variety of processes altering the physiology of Synechocystis cells.

Benefit from decline: the primary transcriptome of Alteromonas macleodii str. Te101 during Trichodesmium demise

Interactions between co-existing microorganisms deeply affect the physiology of the involved organisms and, ultimately, the function of the ecosystem as a whole. Copiotrophic Alteromonas are marine gammaproteobacteria that thrive during the late stages of phytoplankton blooms in the marine environment and in laboratory co-cultures with cyanobacteria such as Trichodesmium. The response of this heterotroph to the sometimes rapid and transient changes in nutrient supply when the phototroph crashes is not well understood. Here, we isolated and sequenced the strain Alteromonas macleodii str. Te101 from a laboratory culture of Trichodesmium erythraeum IMS101, yielding a chromosome of 4.63 Mb and a single plasmid of 237 kb. Increasing salinities to ≥43 ppt inhibited the growth of Trichodesmium but stimulated growth of the associated Alteromonas. We characterized the transcriptomic responses of both microorganisms and identified the complement of active transcriptional start sites in Alteromonas at single-nucleotide resolution. In replicate cultures, a similar set of genes became activated in Alteromonas when growth rates of Trichodesmium declined and mortality was high. The parallel activation of fliA, rpoS and of flagellar assembly and growth-related genes indicated that Alteromonas might have increased cell motility, growth, and multiple biosynthetic activities. Genes with the highest expression in the data set were three small RNAs (Aln1a-c) that were identified as analogs of the small RNAs CsrB-C in E. coli or RsmX-Z in pathogenic bacteria. Together with the carbon storage protein A (CsrA) homolog Te101_05290, these RNAs likely control the expression of numerous genes in responding to changes in the environment.

The primary transcriptome of the fast-growing cyanobacterium Synechococcus elongatus UTEX 2973

BACKGROUND

Cyanobacteria have shown promising potential for the production of various biofuels and chemical feedstocks. Synechococcus elongatus UTEX 2973 is a fast-growing strain with pronounced tolerance to high temperatures and illumination. Hence, this strain appears to be ideal for the development of photosynthetic biotechnology. However, molecular insights on how this strain can rapidly accumulate biomass and carbohydrates under high-light and high-temperature conditions are lacking.

RESULTS

Differential RNA-Sequencing (dRNA-Seq) enabled the genome-wide identification of 4808 transcription start sites (TSSs) in S. elongatus UTEX 2973 using a background reduction algorithm. High light promoted the transcription of genes associated with central metabolic pathways, whereas the highly induced small RNA (sRNA) PsrR1 likely contributed to the repression of phycobilisome genes and the accelerated glycogen accumulation rates measured under this condition. Darkness caused transcriptome remodeling with a decline in the expression of genes for carbon fixation and other major metabolic pathways and an increase in the expression of genes for glycogen catabolism and Calvin cycle inhibitor CP12. Two of the identified TSSs drive the transcription of highly abundant sRNAs in darkness. One of them is widely conserved throughout the cyanobacterial phylum. Its gene is fused to a protein-coding gene in some species, illustrating the evolutionary origin of sRNAs from an mRNA 3'-end.

CONCLUSIONS

Our comprehensive set of genome-wide mapped TSSs, sRNAs and promoter activities will be valuable for projects requiring precise information about the control of transcription aimed at metabolic engineering and the elucidation of stress acclimation mechanisms in this promising strain.

Termination factor Rho: From the control of pervasive transcription to cell fate determination in Bacillus subtilis

In eukaryotes, RNA species originating from pervasive transcription are regulators of various cellular processes, from the expression of individual genes to the control of cellular development and oncogenesis. In prokaryotes, the function of pervasive transcription and its output on cell physiology is still unknown. Most bacteria possess termination factor Rho, which represses pervasive, mostly antisense, transcription. Here, we investigate the biological significance of Rho-controlled transcription in the Gram-positive model bacterium Bacillus subtilis. Rho inactivation strongly affected gene expression in B. subtilis, as assessed by transcriptome and proteome analysis of a rho-null mutant during exponential growth in rich medium. Subsequent physiological analyses demonstrated that a considerable part of Rho-controlled transcription is connected to balanced regulation of three mutually exclusive differentiation programs: cell motility, biofilm formation, and sporulation. In the absence of Rho, several up-regulated sense and antisense transcripts affect key structural and regulatory elements of these differentiation programs, thereby suppressing motility and biofilm formation and stimulating sporulation. We dissected how Rho is involved in the activity of the cell fate decision-making network, centered on the master regulator Spo0A. We also revealed a novel regulatory mechanism of Spo0A activation through Rho-dependent intragenic transcription termination of the protein kinase kinB gene. Altogether, our findings indicate that distinct Rho-controlled transcripts are functional and constitute a previously unknown built-in module for the control of cell differentiation in B. subtilis. In a broader context, our results highlight the recruitment of the termination factor Rho, for which the conserved biological role is probably to repress pervasive transcription, in highly integrated, bacterium-specific, regulatory networks.

Global Transcriptional Responses to Osmotic, Oxidative, and Imipenem Stress Conditions in Pseudomonas putida

Bacteria cope with and adapt to stress by modulating gene expression in response to specific environmental cues. In this study, the transcriptional response of Pseudomonas putida KT2440 to osmotic, oxidative, and imipenem stress conditions at two time points was investigated via identification of differentially expressed mRNAs and small RNAs (sRNAs). A total of 440 sRNA transcripts were detected, of which 10% correspond to previously annotated sRNAs, 40% to novel intergenic transcripts, and 50% to novel transcripts antisense to annotated genes. Each stress elicits a unique response as far as the extent and dynamics of the transcriptional changes. Nearly 200 protein-encoding genes exhibited significant changes in all stress types, implicating their participation in a general stress response. Almost half of the sRNA transcripts were differentially expressed under at least one condition, suggesting possible functional roles in the cellular response to stress conditions. The data show a larger fraction of differentially expressed sRNAs than of mRNAs with >5-fold expression changes. The work provides detailed insights into the mechanisms through which P. putida responds to different stress conditions and increases understanding of bacterial adaptation in natural and industrial settings.IMPORTANCE This study maps the complete transcriptional response of P. putida KT2440 to osmotic, oxidative, and imipenem stress conditions at short and long exposure times. Over 400 sRNA transcripts, consisting of both intergenic and antisense transcripts, were detected, increasing the number of identified sRNA transcripts in the strain by a factor of 10. Unique responses to each type of stress are documented, including both the extent and dynamics of the gene expression changes. The work adds rich detail to previous knowledge of stress response mechanisms due to the depth of the RNA sequencing data. Almost half of the sRNAs exhibit significant expression changes under at least one condition, suggesting their involvement in adaptation to stress conditions and identifying interesting candidates for further functional characterization.

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.

mdRNA-Seq analysis of marine microbial communities from the northern Red Sea

Metatranscriptomic differential RNA-Seq (mdRNA-Seq) identifies the suite of active transcriptional start sites at single-nucleotide resolution through enrichment of primary transcript 5′ ends. Here we analyzed the microbial community at 45 m depth at Station A in the northern Gulf of Aqaba, Red Sea, during 500 m deep mixing in February 2012 using mdRNA-Seq and a parallel classical RNA-Seq approach. We identified promoters active in situ for five different pico-planktonic genera (the SAR11 clade of Alphaproteobacteria, Synechococcus of Cyanobacteria, Euryarchaeota, Thaumarchaeota, and Micromonas as an example for picoeukaryotic algae), showing the applicability of this approach to highly diverse microbial communities. 16S rDNA quantification revealed that 24% of the analyzed community were group II marine Euryarchaeota in which we identified a highly abundant non-coding RNA, Tan1, and detected very high expression of genes encoding intrinsically disordered proteins, as well as enzymes for the synthesis of specific B vitamins, extracellular peptidases, carbohydrate-active enzymes, and transport systems. These results highlight previously unknown functions of Euryarchaeota with community-wide relevance. The complementation of metatranscriptomic studies with mdRNA-Seq provides substantial additional information regarding transcriptional start sites, promoter activities, and the identification of non-coding RNAs.

Natural antisense RNAs as mRNA regulatory elements in bacteria: a review on function and applications

Naturally occurring antisense RNAs are small, diffusible, untranslated transcripts that pair to target RNAs at specific regions of complementarity to control their biological function by regulating gene expression at the post-transcriptional level. This review focuses on known cases of antisense RNA control in prokaryotes and provides an overview of some natural RNA-based mechanisms that bacteria use to modulate gene expression, such as mRNA sensors, riboswitches and antisense RNAs. We also highlight recent advances in RNA-based technology. The review shows that studies on both natural and synthetic systems are reciprocally beneficial.

Type II Toxin-Antitoxin Systems in the Unicellular Cyanobacterium Synechocystis sp. PCC 6803

Bacterial toxin–antitoxin (TA) systems are genetic elements, which are encoded by plasmid as well as chromosomal loci. They mediate plasmid and genomic island maintenance through post-segregational killing mechanisms but may also have milder effects, acting as mobile stress response systems that help certain cells of a population in persisting adverse growth conditions. Very few cyanobacterial TA system have been characterized thus far. In this work, we focus on the cyanobacterium Synechocystis 6803, a widely used model organism. We expand the number of putative Type II TA systems from 36 to 69 plus seven stand-alone components. Forty-seven TA pairs are located on the chromosome and 22 are plasmid-located. Different types of toxins are associated with various antitoxins in a mix and match principle. According to protein domains and experimental data, 81% of all toxins in Synechocystis 6803 likely exhibit RNase activity, suggesting extensive potential for toxicity-related RNA degradation and toxin-mediated transcriptome remodeling. Of particular interest is the Ssr8013–Slr8014 system encoded on plasmid pSYSG, which is part of a larger defense island or the pSYSX system Slr6056–Slr6057, which is linked to a bacterial ubiquitin-like system. Consequently, Synechocystis 6803 is one of the most prolific sources of new information about these genetic elements.

A Feed-Forward Loop Consisting of the Response Regulator RpaB and the Small RNA PsrR1 Controls Light Acclimation of Photosystem I Gene Expression in the Cyanobacterium Synechocystis sp. PCC 6803

Since cyanobacteria need to decrease PSI content to avoid absorption of excess light energy, down-regulation of PSI gene expression is one of the key characteristics of the high-light (HL) acclimation response. The transcriptional regulator RpaB and the small RNA PsrR1 (photosynthesis regulatory RNA1) have been suggested to be the two most critical factors for this response in Synechocystis sp. PCC 6803. In this study, we found that the HLR1 DNA-binding motif, the recognition sequence for RpaB, is highly conserved in the core promoter region of the psrR1 gene among cyanobacterial species. Gel mobility shift assay revealed that RpaB binds to the HLR1 sequence of psrR1 in vitro. RNA gel blot analysis together with chromatin affinity purification (ChAP) analysis suggested that PSI genes are activated and the psrR1 gene is repressed by the binding of RpaB under low-light (LL) conditions. A decrease in DNA binding affinity of RpaB occurs within 5 min after the shift from LL to HL conditions, leading to the prompt decrease in PSI promoter activity together with derepression of psrR1 gene expression. Accumulating PsrR1 molecules then prevent translation from pre-existing PSI transcripts. By this dual repression at transcriptional and post-transcriptional levels, rapid and strict down-regulation of PSI expression under HL is secured. Our findings suggest that RpaB and PsrR1 constitute a feed-forward loop for the regulation of PSI gene expression to achieve a rapid acclimation response to the damaging HL conditions.

Advances in genomics, transcriptomics and proteomics of toxin-producing cyanobacteria

A common misconception persists that the genomes of toxic and non-toxic cyanobacterial strains are largely conserved with the exception of the presence or absence of the genes responsible for toxin production. Implementation of -omics era technologies has challenged this paradigm, with comparative analyses providing increased insight into the differences between strains of the same species. The implementation of genomic, transcriptomic and proteomic approaches has revealed distinct profiles between toxin-producing and non-toxic strains. Further, metagenomics and metaproteomics highlight the genomic potential and functional state of toxic bloom events over time. In this review, we highlight how these technologies have shaped our understanding of the complex relationship between these molecules, their producers and the environment at large within which they persist.

Integrated Transcriptomic and Metabolomic Characterization of the Low-Carbon Response Using an ndhR Mutant of Synechocystis sp. PCC 6803

The acquisition and assimilation of inorganic carbon (Ci) represents the largest flux of inorganic matter in photosynthetic organisms; hence, this process is tightly regulated. We examined the Ci-dependent transcriptional and metabolic regulation in wild-type Synechocystis sp. PCC 6803 compared with a mutant defective in the main transcriptional repressor for Ci acquisition genes, the NAD(P)H dehydrogenase transcriptional regulator NdhR. The analysis revealed that many protein-coding transcripts that are normally repressed in the presence of high CO2 (HC) concentrations were strongly expressed in ∆ndhR, whereas other messenger RNAs were strongly down-regulated in mutant cells, suggesting a potential activating role for NdhR. A conserved NdhR-binding motif was identified in the promoters of derepressed genes. Interestingly, the expression of some NdhR-regulated genes remained further inducible under low-CO2 conditions, indicating the involvement of additional NdhR-independent Ci-regulatory mechanisms. Intriguingly, we also observed that the abundance of 52 antisense RNAs and 34 potential noncoding RNAs was affected by Ci supply, although most of these molecules were not regulated through NdhR. Thus, antisense and noncoding RNAs could contribute to NdhR-independent carbon regulation. In contrast to the transcriptome, the metabolome in ∆ndhR cells was similar to that of wild-type cells under HC conditions. This observation and the delayed metabolic responses to the low-CO2 shift in ∆ndhR, specifically the lack of transient increases in the photorespiratory pathway intermediates 2-phosphoglycolate, glycolate, and glycine, suggest that the deregulation of gene expression in the ΔndhR mutant successfully preacclimates cyanobacterial cells to lowered Ci supply under HC conditions.

The sRNA NsiR4 is involved in nitrogen assimilation control in cyanobacteria by targeting glutamine synthetase inactivating factor IF7

Glutamine synthetase (GS), a key enzyme in biological nitrogen assimilation, is regulated in multiple ways in response to varying nitrogen sources and levels. Here we show a small regulatory RNA, NsiR4 (nitrogen stress-induced RNA 4), which plays an important role in the regulation of GS in cyanobacteria. NsiR4 expression in the unicellular Synechocystis sp. PCC 6803 and in the filamentous, nitrogen-fixing Anabaena sp. PCC 7120 is stimulated through nitrogen limitation via NtcA, the global transcriptional regulator of genes involved in nitrogen metabolism. NsiR4 is widely conserved throughout the cyanobacterial phylum, suggesting a conserved function. In silico target prediction, transcriptome profiling on pulse overexpression, and site-directed mutagenesis experiments using a heterologous reporter system showed that NsiR4 interacts with the 5'UTR of gifA mRNA, which encodes glutamine synthetase inactivating factor (IF)7. In Synechocystis, we observed an inverse relationship between the levels of NsiR4 and the accumulation of IF7 in vivo. This NsiR4-dependent modulation of gifA (IF7) mRNA accumulation influenced the glutamine pool and thus [Formula: see text] assimilation via GS. As a second target, we identified ssr1528, a hitherto uncharacterized nitrogen-regulated gene. Competition experiments between WT and an ΔnsiR4 KO mutant showed that the lack of NsiR4 led to decreased acclimation capabilities of Synechocystis toward oscillating nitrogen levels. These results suggest a role for NsiR4 in the regulation of nitrogen metabolism in cyanobacteria, especially for the adaptation to rapid changes in available nitrogen sources and concentrations. NsiR4 is, to our knowledge, the first identified bacterial sRNA regulating the primary assimilation of a macronutrient.

Transcriptomic and Proteomic Profiling of Anabaena sp. Strain 90 under Inorganic Phosphorus Stress

Inorganic phosphorus (Pi) is one of the main growth-limiting factors of diazotrophic cyanobacteria. Due to human activity, the availability of Pi has increased in water bodies, resulting in eutrophication and the formation of massive cyanobacterial blooms. In this study, we examined the molecular responses of the cyanobacterium Anabaena sp. strain 90 to phosphorus deprivation, aiming at the identification of candidate genes to monitor the Pi status in cyanobacteria. Furthermore, this study increased the basic understanding of how phosphorus affects diazotrophic and bloom-forming cyanobacteria as a major growth-limiting factor. Based on RNA sequencing data, we identified 246 differentially expressed genes after phosphorus starvation and 823 differentially expressed genes after prolonged Pi limitation, most of them related to central metabolism and cellular growth. The transcripts of the genes related to phosphorus transport and assimilation (pho regulon) were most upregulated during phosphorus depletion. One of the most increased transcripts encodes a giant protein of 1,869 amino acid residues, which contains, among others, a phytase-like domain. Our findings predict its crucial role in phosphorus starvation, but future studies are still needed. Using two-dimensional difference in gel electrophoresis (2D-DIGE) and liquid chromatography-tandem mass spectrometry (LC-MS/MS), we found 43 proteins that were differentially expressed after prolonged phosphorus stress. However, correlation analysis unraveled an association only to some extent between the transcriptomic and proteomic abundances. Based on the present results, we suggest that the method used for monitoring the Pi status in cyanobacterial bloom should contain wider combinations of pho regulon genes (e.g., PstABCS transport systems) in addition to the commonly used alkaline phosphatase gene alone.

Regulatory RNAs in photosynthetic cyanobacteria

Regulatory RNAs play versatile roles in bacteria in the coordination of gene expression during various physiological processes, especially during stress adaptation. Photosynthetic bacteria use sunlight as their major energy source. Therefore, they are particularly vulnerable to the damaging effects of excess light or UV irradiation. In addition, like all bacteria, photosynthetic bacteria must adapt to limiting nutrient concentrations and abiotic and biotic stress factors. Transcriptome analyses have identified hundreds of potential regulatory small RNAs (sRNAs) in model cyanobacteria such as Synechocystis sp. PCC 6803 or Anabaena sp. PCC 7120, and in environmentally relevant genera such as Trichodesmium, Synechococcus and Prochlorococcus. Some sRNAs have been shown to actually contain μORFs and encode short proteins. Examples include the 40-amino-acid product of the sml0013 gene, which encodes the NdhP subunit of the NDH1 complex. In contrast, the functional characterization of the non-coding sRNA PsrR1 revealed that the 131 nt long sRNA controls photosynthetic functions by targeting multiple mRNAs, providing a paradigm for sRNA functions in photosynthetic bacteria. We suggest that actuatons comprise a new class of genetic elements in which an sRNA gene is inserted upstream of a coding region to modify or enable transcription of that region.