Global mapping of transcription start sites and promoter motifs in the symbiotic α-proteobacterium Sinorhizobium meliloti 1021

Novel Genes and Regulators That Influence Production of Cell Surface Exopolysaccharides in Sinorhizobium meliloti

Sinorhizobium meliloti is a soil-dwelling alphaproteobacterium that engages in a nitrogen-fixing root nodule symbiosis with leguminous plants. Cell surface polysaccharides are important both for adapting to stresses in the soil and for the development of an effective symbiotic interaction. Among the polysaccharides characterized to date, the acidic exopolysaccharides I (EPS-I; succinoglycan) and II (EPS-II; galactoglucan) are particularly important for protection from abiotic stresses, biofilm formation, root colonization, and infection of plant roots. Previous genetic screens discovered mutants with impaired EPS production, allowing the delineation of EPS biosynthetic pathways. Here we report on a genetic screen to isolate mutants with mucoid colonial morphologies that suggest EPS overproduction. Screening with Tn5-110, which allows the recovery of both null and upregulation mutants, yielded 47 mucoid mutants, most of which overproduce EPS-I; among the 30 unique genes and intergenic regions identified, 14 have not been associated with EPS production previously. We identified a new protein-coding gene, emmD, which may be involved in the regulation of EPS-I production as part of the EmmABC three-component regulatory circuit. We also identified a mutant defective in EPS-I production, motility, and symbiosis, where Tn5-110 was not responsible for the mutant phenotypes; these phenotypes result from a missense mutation in rpoA corresponding to the domain of the RNA polymerase alpha subunit known to interact with transcription regulators.IMPORTANCE The alphaproteobacterium Sinorhizobium meliloti converts dinitrogen to ammonium while inhabiting specialized plant organs termed root nodules. The transformation of S. meliloti from a free-living soil bacterium to a nitrogen-fixing plant symbiont is a complex developmental process requiring close interaction between the two partners. As the interface between the bacterium and its environment, the S. meliloti cell surface plays a critical role in adaptation to varied soil environments and in interaction with plant hosts. We isolated and characterized S. meliloti mutants with increased production of exopolysaccharides, key cell surface components. Our diverse set of mutants suggests roles for exopolysaccharide production in growth, metabolism, cell division, envelope homeostasis, biofilm formation, stress response, motility, and symbiosis.

RNAseq analysis of α-proteobacterium Gluconobacter oxydans 621H

Background

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

Results

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

Conclusions

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

Electronic supplementary material

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

Sinorhizobium meliloti YbeY is an endoribonuclease with unprecedented catalytic features, acting as silencing enzyme in riboregulation

Structural and biochemical features suggest that the almost ubiquitous bacterial YbeY protein may serve catalytic and/or Hfq-like protective functions central to small RNA (sRNA)-mediated regulation and RNA metabolism. We have biochemically and genetically characterized the YbeY ortholog of the legume symbiont Sinorhizobium meliloti (SmYbeY). Co-immunoprecipitation (CoIP) with a FLAG-tagged SmYbeY yielded a poor enrichment in RNA species, compared to Hfq CoIP-RNA uncovered previously by a similar experimental setup. Purified SmYbeY behaved as a monomer that indistinctly cleaved single- and double-stranded RNA substrates, a unique ability among bacterial endoribonucleases. SmYbeY-mediated catalysis was supported by the divalent metal ions Mg²⁺, Mn²⁺ and Ca²⁺, which influenced in a different manner cleavage efficiency and reactivity patterns, with Ca²⁺ specifically blocking activity on double-stranded and some structured RNA molecules. SmYbeY loss-of-function compromised expression of core energy and RNA metabolism genes, whilst promoting accumulation of motility, late symbiotic and transport mRNAs. Some of the latter transcripts are known Hfq-binding sRNA targets and might be SmYbeY substrates. Genetic reporter and in vitro assays confirmed that SmYbeY is required for sRNA-mediated down-regulation of the amino acid ABC transporter prbA mRNA. We have thus discovered a bacterial endoribonuclease with unprecedented catalytic features, acting also as gene silencing enzyme.

Genome-Wide Analysis of the Transcription Start Sites and Promoter Motifs of Phytoplasmas

Phytoplasmas are obligate intracellular parasitic bacteria that infect both plants and insects. We previously identified the sigma factor RpoD-dependent consensus promoter sequence of phytoplasma. However, the genome-wide landscape of RNA transcripts, including non-coding RNAs (ncRNAs) and RpoD-independent promoter elements, was still unknown. In this study, we performed an improved RNA sequencing analysis for genome-wide identification of the transcription start sites (TSSs) and the consensus promoter sequences. We constructed cDNA libraries using a random adenine/thymine hexamer primer, in addition to a conventional random hexamer primer, for efficient sequencing of 5'-termini of AT-rich phytoplasma RNAs. We identified 231 TSSs, which were classified into four categories: mRNA TSSs, internal sense TSSs, antisense TSSs (asTSSs), and orphan TSSs (oTSSs). The presence of asTSSs and oTSSs indicated the genome-wide transcription of ncRNAs, which might act as regulatory ncRNAs in phytoplasmas. This is the first description of genome-wide phytoplasma ncRNAs. Using a de novo motif discovery program, we identified two consensus motif sequences located upstream of the TSSs. While one was almost identical to the RpoD-dependent consensus promoter sequence, the other was an unidentified novel motif, which might be recognized by another transcription initiation factor. These findings are valuable for understanding the regulatory mechanism of phytoplasma gene expression.

RNA silencing in plant symbiotic bacteria: Insights from a protein-centric view

Extensive work in model enterobacteria has evidenced that the RNA chaperone Hfq and several endoribonucleases, such as RNase E or RNase III, serve pivotal roles in small RNA-mediated post-transcriptional silencing of gene expression. Characterization of these protein hubs commonly provide global functional and mechanistic insights into complex sRNA regulatory networks. The legume endosymbiont Sinorhizobium meliloti is a non-classical model bacterium with a very complex lifestyle in which riboregulation is expected to play important adaptive functions. Here, we discuss current knowledge about RNA silencing in S. meliloti from the perspective of the activity of Hfq and a recently discovered endoribonuclease (YbeY) exhibiting unprecedented catalytic versatility for the cleavage of single- and double-stranded RNA molecules.

Transcriptome Analysis of Polyhydroxybutyrate Cycle Mutants Reveals Discrete Loci Connecting Nitrogen Utilization and Carbon Storage in Sinorhizobium meliloti

The ability of bacteria to store carbon and energy as intracellular polymers uncouples cell growth and replication from nutrient uptake and provides flexibility in the use of resources as they are available to the cell. The impact of carbon storage on cellular metabolism would be reflected in global transcription patterns. By investigating the transcriptomic effects of genetically disrupting genes involved in the PHB carbon storage cycle, we revealed a relationship between intracellular carbon storage and nitrogen metabolism. This work demonstrates the utility of combining transcriptome sequencing with metabolic pathway mutations for identifying underlying gene regulatory mechanisms.

Polyhydroxybutyrate (PHB) and glycogen polymers are produced by bacteria as carbon storage compounds under unbalanced growth conditions. To gain insights into the transcriptional mechanisms controlling carbon storage in Sinorhizobium meliloti, we investigated the global transcriptomic response to the genetic disruption of key genes in PHB synthesis and degradation and in glycogen synthesis. Under both nitrogen-limited and balanced growth conditions, transcriptomic analysis was performed with genetic mutants deficient in PHB synthesis (phbA, phbB, phbAB, and phbC), PHB degradation (bdhA, phaZ, and acsA2), and glycogen synthesis (glgA1). Three distinct genomic regions of the pSymA megaplasmid exhibited altered expression in the wild type and the PHB cycle mutants that was not seen in the glycogen synthesis mutant. An Fnr family transcriptional motif was identified in the upstream regions of a cluster of genes showing similar transcriptional patterns across the mutants. This motif was found at the highest density in the genomic regions with the strongest transcriptional effect, and the presence of this motif upstream of genes in these regions was significantly correlated with decreased transcript abundance. Analysis of the genes in the pSymA regions revealed that they contain a genomic overrepresentation of Fnr family transcription factor-encoding genes. We hypothesize that these loci, containing mostly nitrogen utilization, denitrification, and nitrogen fixation genes, are regulated in response to the intracellular carbon/nitrogen balance. These results indicate a transcriptional regulatory association between intracellular carbon levels (mediated through the functionality of the PHB cycle) and the expression of nitrogen metabolism genes.

IMPORTANCE The ability of bacteria to store carbon and energy as intracellular polymers uncouples cell growth and replication from nutrient uptake and provides flexibility in the use of resources as they are available to the cell. The impact of carbon storage on cellular metabolism would be reflected in global transcription patterns. By investigating the transcriptomic effects of genetically disrupting genes involved in the PHB carbon storage cycle, we revealed a relationship between intracellular carbon storage and nitrogen metabolism. This work demonstrates the utility of combining transcriptome sequencing with metabolic pathway mutations for identifying underlying gene regulatory mechanisms.

Author Video: An author video summary of this article is available.

A Key Regulator of the Glycolytic and Gluconeogenic Central Metabolic Pathways in Sinorhizobium meliloti

The order Rhizobiales contains numerous agriculturally, biotechnologically, and medically important bacteria, including the rhizobia, and the genera Agrobacterium, Brucella, and Methylobacterium, among others. These organisms tend to be metabolically versatile, but there has been relatively little investigation into the regulation of their central carbon metabolic pathways. Here, RNA-sequencing and promoter fusion data are presented to show that the PckR protein is a key regulator of central carbon metabolism in Sinorhizobium meliloti; during growth with gluconeogenic substrates, PckR represses expression of the complete Entner-Doudoroff glycolytic pathway and induces expression of the pckA and fbaB gluconeogenic genes. Electrophoretic mobility shift assays indicate that PckR binds an imperfect palindromic sequence that overlaps the promoter or transcriptional start site in the negatively regulated promoters, or is present in tandem upstream the promoter motifs in the positively regulated promoters. Genetic and in vitro electrophoretic mobility shift assay experiments suggest that elevated concentrations of a PckR effector ligand results in the dissociation of PckR from its target binding site, and evidence is presented that suggests phosphoenolpyruvate may function as the effector. Characterization of missense pckR alleles identified three conserved residues important for increasing the affinity of PckR for its cognate effector molecule. Bioinformatics analyses illustrates that PckR is limited to a narrow phylogenetic range consisting of the Rhizobiaceae, Phyllobacteriaceae, Brucellaceae, and Bartonellaceae families. These data provide novel insights into the regulation of the core carbon metabolic pathways of this pertinent group of α-proteobacteria.

The Irr and RirA Proteins Participate in a Complex Regulatory Circuit and Act in Concert To Modulate Bacterioferritin Expression in Ensifer meliloti 1021

In this work we found that the bfr gene of the rhizobial species Ensifer meliloti, encoding a bacterioferritin iron storage protein, is involved in iron homeostasis and the oxidative stress response. This gene is located downstream of and overlapping the smc03787 open reading frame (ORF). No well-predicted RirA or Irr boxes were found in the region immediately upstream of the bfr gene although two presumptive RirA boxes and one presumptive Irr box were present in the putative promoter of smc03787 We demonstrate that bfr gene expression is enhanced under iron-sufficient conditions and that Irr and RirA modulate this expression. The pattern of bfr gene expression as well as the response to Irr and RirA is inversely correlated to that of smc03787 Moreover, our results suggest that the small RNA SmelC759 participates in RirA- and Irr-mediated regulation of bfr expression and that additional unknown factors are involved in iron-dependent regulation.IMPORTANCEE. meliloti belongs to the Alphaproteobacteria, a group of bacteria that includes several species able to associate with eukaryotic hosts, from mammals to plants, in a symbiotic or pathogenic manner. Regulation of iron homeostasis in this group of bacteria differs from that found in the well-studied Gammaproteobacteria In this work we analyzed the effect of rirA and irr mutations on bfr gene expression. We demonstrate the effect of an irr mutation on iron homeostasis in this bacterial genus. Moreover, results obtained indicate a complex regulatory circuit where multiple regulators, including RirA, Irr, the small RNA SmelC759, and still unknown factors, act in concert to balance bfr gene expression.

The NtrY/NtrX System of Sinorhizobium meliloti GR4 Regulates Motility, EPS I Production, and Nitrogen Metabolism but Is Dispensable for Symbiotic Nitrogen Fixation

Sinorhizobium meliloti can translocate over surfaces. However, little is known about the regulatory mechanisms that control this trait and its relevance for establishing symbiosis with alfalfa plants. To gain insights into this field, we isolated Tn5 mutants of S. meliloti GR4 with impaired surface motility. In mutant strain GRS577, the transposon interrupted the ntrY gene encoding the sensor kinase of the NtrY/NtrX two-component regulatory system. GRS577 is impaired in flagella synthesis and overproduces succinoglycan, which is responsible for increased biofilm formation. The mutant also shows altered cell morphology and higher susceptibility to salt stress. GRS577 induces nitrogen-fixing nodules in alfalfa but exhibits decreased competitive nodulation. Complementation experiments indicate that both ntrY and ntrX account for all the phenotypes displayed by the ntrY::Tn5 mutant. Ectopic overexpression of VisNR, the motility master regulator, was sufficient to rescue motility and competitive nodulation of the transposant. A transcriptome profiling of GRS577 confirmed differential expression of exo and flagellar genes, and led to the demonstration that NtrY/NtrX allows for optimal expression of denitrification and nifA genes under microoxic conditions in response to nitrogen compounds. This study extends our knowledge of the complex role played by NtrY/NtrX in S. meliloti.

A conserved α-proteobacterial small RNA contributes to osmoadaptation and symbiotic efficiency of rhizobia on legume roots

Small non-coding RNAs (sRNAs) are expected to have pivotal roles in the adaptive responses underlying symbiosis of nitrogen-fixing rhizobia with legumes. Here, we provide primary insights into the function and activity mechanism of the Sinorhizobium meliloti trans-sRNA NfeR1 (Nodule Formation Efficiency RNA). Northern blot probing and transcription tracking with fluorescent promoter-reporter fusions unveiled high nfeR1 expression in response to salt stress and throughout the symbiotic interaction. The strength and differential regulation of nfeR1 transcription are conferred by a motif, which is conserved in nfeR1 promoter regions in α-proteobacteria. NfeR1 loss-of-function compromised osmoadaptation of free-living bacteria, whilst causing misregulation of salt-responsive genes related to stress adaptation, osmolytes catabolism and membrane trafficking. Nodulation tests revealed that lack of NfeR1 affected competitiveness, infectivity, nodule development and symbiotic efficiency of S. meliloti on alfalfa roots. Comparative computer predictions and a genetic reporter assay evidenced a redundant role of three identical NfeR1 unpaired anti Shine-Dalgarno motifs for targeting and downregulation of translation of multiple mRNAs from transporter genes. Our data provide genetic evidence of the hyperosmotic conditions of the endosymbiotic compartments. NfeR1-mediated gene regulation in response to this cue could contribute to coordinate nutrient uptake with the metabolic reprogramming concomitant to symbiotic transitions.

Deep sequencing approaches for the analysis of prokaryotic transcriptional boundaries and dynamics

The identification of the protein-coding regions of a genome is straightforward due to the universality of start and stop codons. However, the boundaries of the transcribed regions, conditional operon structures, non-coding RNAs and the dynamics of transcription, such as pausing of elongation, are non-trivial to identify, even in the comparatively simple genomes of prokaryotes. Traditional methods for the study of these areas, such as tiling arrays, are noisy, labour-intensive and lack the resolution required for densely-packed bacterial genomes. Recently, deep sequencing has become increasingly popular for the study of the transcriptome due to its lower costs, higher accuracy and single nucleotide resolution. These methods have revolutionised our understanding of prokaryotic transcriptional dynamics. Here, we review the deep sequencing and data analysis techniques that are available for the study of transcription in prokaryotes, and discuss the bioinformatic considerations of these analyses.

Heterologous Complementation Reveals a Specialized Activity for BacA in the Medicago-Sinorhizobium meliloti Symbiosis

The bacterium Sinorhizobium meliloti Rm2011 forms N2-fixing root nodules on alfalfa and other leguminous plants. The pSymB chromid contains a 110-kb region (the ETR region) showing high synteny to a chromosomally located region in Sinorhizobium fredii NGR234 and related rhizobia. We recently introduced the ETR region from S. fredii NGR234 into the S. meliloti chromosome. Here, we report that, unexpectedly, the S. fredii NGR234 ETR region did not complement deletion of the S. meliloti ETR region in symbiosis with Medicago sativa. This phenotype was due to the bacA gene of NGR234 not being functionally interchangeable with the S. meliloti bacA gene during M. sativa symbiosis. Further analysis revealed that, whereas bacA genes from S. fredii or Rhizobium leguminosarum bv. viciae 3841 failed to complement the Fix- phenotype of a S. meliloti bacA mutant with M. sativa, they allowed for further developmental progression prior to a loss of viability. In contrast, with Melilotus alba, bacA from S. fredii and R. leguminosarum supported N2 fixation by a S. meliloti bacA mutant. Additionally, the S. meliloti bacA gene can support N2 fixation of a R. leguminosarum bacA mutant during symbiosis with Pisum sativum. A phylogeny of BacA proteins illustrated that S. meliloti BacA has rapidly diverged from most rhizobia and has converged toward the sequence of pathogenic genera Brucella and Escherichia. These data suggest that the S. meliloti BacA has evolved toward a specific interaction with Medicago and highlights the limitations of using a single model system for the study of complex biological topics.

Regulation of Polyhydroxybutyrate Accumulation in Sinorhizobium meliloti by the Trans-Encoded Small RNA MmgR

Riboregulation has a major role in the fine-tuning of multiple bacterial processes. Among the RNA players, trans-encoded untranslated small RNAs (sRNAs) regulate complex metabolic networks by tuning expression from multiple target genes in response to numerous signals. In Sinorhizobium meliloti, over 400 sRNAs are expressed under different stimuli. The sRNA MmgR (standing for Makes more granules Regulator) has been of particular interest to us since its sequence and structure are highly conserved among the alphaproteobacteria and its expression is regulated by the amount and quality of the bacterium's available nitrogen source. In this work, we explored the biological role of MmgR in S. meliloti 2011 by characterizing the effect of a deletion of the internal conserved core of mmgR (mmgR^Δ33-51). This mutation resulted in larger amounts of polyhydroxybutyrate (PHB) distributed into more intracellular granules than are found in the wild-type strain. This phenotype was expressed upon cessation of balanced growth owing to nitrogen depletion in the presence of surplus carbon (i.e., at a carbon/nitrogen molar ratio greater than 10). The normal PHB accumulation was complemented with a wild-type mmgR copy but not with unrelated sRNA genes. Furthermore, the expression of mmgR limited PHB accumulation in the wild type, regardless of the magnitude of the C surplus. Quantitative proteomic profiling and quantitative reverse transcription-PCR (qRT-PCR) revealed that the absence of MmgR results in a posttranscriptional overexpression of both PHB phasin proteins (PhaP1 and PhaP2). Together, our results indicate that the widely conserved alphaproteobacterial MmgR sRNA fine-tunes the regulation of PHB storage in S. melilotiIMPORTANCE High-throughput RNA sequencing has recently uncovered an overwhelming number of trans-encoded small RNAs (sRNAs) in diverse prokaryotes. In the nitrogen-fixing alphaproteobacterial symbiont of alfalfa root nodules Sinorhizobium meliloti, only four out of hundreds of identified sRNA genes have been functionally characterized. Thus, uncovering the biological role of sRNAs currently represents a major issue and one that is particularly challenging because of the usually subtle quantitative regulation contributed by most characterized sRNAs. Here, we have characterized the function of the broadly conserved alphaproteobacterial sRNA gene mmgR in S. meliloti Our results strongly suggest that mmgR encodes a negative regulator of the accumulation of polyhydroxybutyrate, the major carbon and reducing power storage polymer in S. meliloti cells growing under conditions of C/N overbalance.

Caulobacter crescentus CdnL is a non-essential RNA polymerase-binding protein whose depletion impairs normal growth and rRNA transcription

CdnL is an essential RNA polymerase (RNAP)-binding activator of rRNA transcription in mycobacteria and myxobacteria but reportedly not in Bacillus. Whether its function and mode of action are conserved in other bacteria thus remains unclear. Because virtually all alphaproteobacteria have a CdnL homolog and none of these have been characterized, we studied the homolog (CdnL_Cc) of the model alphaproteobacterium Caulobacter crescentus. We show that CdnL_Cc is not essential for viability but that its absence or depletion causes slow growth and cell filamentation. CdnL_Cc is degraded in vivo in a manner dependent on its C-terminus, yet excess CdnL_Cc resulting from its stabilization did not adversely affect growth. We find that CdnL_Cc interacts with itself and with the RNAP β subunit, and localizes to at least one rRNA promoter in vivo, whose activity diminishes upon depletion of CdnL_Cc. Interestingly, cells expressing CdnL_Cc mutants unable to interact with the RNAP were cold-sensitive, suggesting that CdnL_Cc interaction with RNAP is especially required at lower than standard growth temperatures in C. crescentus. Our study indicates that despite limited sequence similarities and regulatory differences compared to its myco/myxobacterial homologs, CdnL_Cc may share similar biological functions, since it affects rRNA synthesis, probably by stabilizing open promoter-RNAP complexes.

Massive parallel insertion site sequencing of an arrayed Sinorhizobium meliloti signature-tagged mini-Tn 5 transposon mutant library

Transposon mutagenesis in conjunction with identification of genomic transposon insertion sites is a powerful tool for gene function studies. We have implemented a protocol for parallel determination of transposon insertion sites by Illumina sequencing involving a hierarchical barcoding method that allowed for tracking back insertion sites to individual clones of an arrayed signature-tagged transposon mutant library. This protocol was applied to further characterize a signature-tagged mini-Tn 5 mutant library comprising about 12,000 mutants of the symbiotic nitrogen-fixing alphaproteobacterium Sinorhizobium meliloti (Pobigaylo et al., 2006; Appl. Environ. Microbiol. 72, 4329-4337). Previously, insertion sites have been determined for 5000 mutants of this library. Combining an adapter-free, inverse PCR method for sequencing library preparation with next generation sequencing, we identified 4473 novel insertion sites, increasing the total number of transposon mutants with known insertion site to 9562. The number of protein-coding genes that were hit at least once by a transposon increased by 1231 to a total number of 3673 disrupted genes, which represents 59% of the predicted protein-coding genes in S. meliloti.

If it transcribes, we can sequence it: mining the complexities of host-pathogen-environment interactions using RNA-seq

Host-pathogen interactions are exceedingly complex because they involve multiple host tissues, often occur in the context of normal microflora, and can span diverse microenvironments. Although decades of gene expression studies have provided detailed insights into infection processes, technical challenges have restricted experiments to single pathogenic species or host tissues. RNA-sequencing (RNA-seq) has revolutionized the study of gene expression because in addition to quantifying transcriptional output, it allows detection and characterization of all transcripts in a genome. Here, we review how refined approaches to RNA-seq are used to map the transcriptional networks that control host-pathogen interactions. These enhanced techniques include dRNA-seq and term-seq for the fine-scale mapping of transcriptional start and termination sites, and dual RNA-seq for simultaneous sequencing of host and bacterial pathogen transcriptomes. Dual RNA-seq experiments are currently limited to in vitro infection systems that do not fully reflect the complexities of the in vivo environment, thus a challenge is to develop in vivo model systems and experimental approaches that address the biological heterogeneity of host environments, followed by the integration of RNA-seq with other genome-scale datasets to identify the transcriptional networks that mediate host-pathogen interactions.

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

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

Global repositioning of transcription start sites in a plant-fermenting bacterium

Bacteria respond to their environment by regulating mRNA synthesis, often by altering the genomic sites at which RNA polymerase initiates transcription. Here, we investigate genome-wide changes in transcription start site (TSS) usage by Clostridium phytofermentans, a model bacterium for fermentation of lignocellulosic biomass. We quantify expression of nearly 10,000 TSS at single base resolution by Capp-Switch sequencing, which combines capture of synthetically capped 5' mRNA fragments with template-switching reverse transcription. We find the locations and expression levels of TSS for hundreds of genes change during metabolism of different plant substrates. We show that TSS reveals riboswitches, non-coding RNA and novel transcription units. We identify sequence motifs associated with carbon source-specific TSS and use them for regulon discovery, implicating a LacI/GalR protein in control of pectin metabolism. We discuss how the high resolution and specificity of Capp-Switch enables study of condition-specific changes in transcription initiation in bacteria.

Negative Regulation of Ectoine Uptake and Catabolism in Sinorhizobium meliloti: Characterization of the EhuR Gene

Ectoine has osmoprotective effects on Sinorhizobium meliloti that differ from its effects in other bacteria. Ectoine does not accumulate in S. meliloti cells; instead, it is degraded. The products of the ehuABCD-eutABCDE operon were previously discovered to be responsible for the uptake and catabolism of ectoine in S. meliloti However, the mechanism by which ectoine is involved in the regulation of the ehuABCD-eutABCDE operon remains unclear. The ehuR gene, which is upstream of and oriented in the same direction as the ehuABCD-eutABCDE operon, encodes a member of the MocR/GntR family of transcriptional regulators. Quantitative reverse transcription-PCR and promoter-lacZ reporter fusion experiments revealed that EhuR represses transcription of the ehuABCD-eutABCDE operon, but this repression is inhibited in the presence of ectoine. Electrophoretic mobility shift assays and DNase I footprinting assays revealed that EhuR bound specifically to the DNA regions overlapping the -35 region of the ehuA promoter and the +1 region of the ehuR promoter. Surface plasmon resonance assays further demonstrated direct interactions between EhuR and the two promoters, although EhuR was found to have higher affinity for the ehuA promoter than for the ehuR promoter. In vitro, DNA binding by EhuR could be directly inhibited by a degradation product of ectoine. Our work demonstrates that EhuR is an important negative transcriptional regulator involved in the regulation of ectoine uptake and catabolism and is likely regulated by one or more end products of ectoine catabolism.

IMPORTANCE

Sinorhizobium meliloti is an important soil bacterium that displays symbiotic interactions with legume hosts. Ectoine serves as a key osmoprotectant for S. meliloti However, ectoine does not accumulate in the cells; rather, it is degraded. In this study, we characterized the transcriptional regulation of the operon responsible for ectoine uptake and catabolism in S. meliloti We identified and characterized the transcription repressor EhuR, which is the first MocR/GntR family member found to be involved in the regulation of compatible solute uptake and catabolism. More importantly, we demonstrated for the first time that an ectoine catabolic end product could modulate EhuR DNA-binding activity. Therefore, this work provides new insights into the unique mechanism of ectoine-induced osmoprotection in S. meliloti.

Cyclic mononucleotide- and Clr-dependent gene regulation in Sinorhizobium meliloti

To identify physiological processes affected by cAMP in the plant-symbiotic nitrogen-fixing α-proteobacterium Sinorhizobium meliloti Rm2011, cAMP levels were artificially increased by overexpression of its cognate adenylate/guanylate cyclase gene cyaJ. This resulted in high accumulation of cAMP in the culture supernatant, decreased swimming motility and increased production of succinoglycan, an exopolysaccharide involved in host invasion. Weaker, similar phenotypic changes were induced by overexpression of cyaB and cyaG1. Effects on swimming motility and succinoglycan production were partially dependent on clr encoding a cyclic AMP receptor-like protein. Transcriptome profiling of an cyaJ-overexpressing strain identified 72 upregulated and 82 downregulated genes. A considerable number of upregulated genes are related to polysaccharide biosynthesis and osmotic stress response. These included succinoglycan biosynthesis genes, genes of the putative polysaccharide synthesis nodP2-exoF3 cluster and feuN, the first gene of the operon encoding the FeuNPQ regulatory system. Downregulated genes were mostly related to respiration, central metabolism and swimming motility. Promoter-probe studies in the presence of externally added cAMP revealed 18 novel Clr-cAMP-regulated genes. Moreover, the addition of cGMP into the growth medium also promoted clr-dependent gene regulation. In vitro binding of Clr-cAMP and Clr-cGMP to the promoter regions of SMc02178, SMb20906,SMc04190, SMc00925, SMc01136 and cyaF2 required the DNA motif (A/C/T)GT(T/C)(T/C/A)C (N4) G(G/A)(T/A)ACA. Furthermore, SMc02178, SMb20906,SMc04190and SMc00653 promoters were activated by Clr-cAMP/cGMP in Escherichia coli as heterologous host. These findings suggest direct activation of these 7 genes by Clr-cAMP/cGMP.

A Sinorhizobium meliloti RpoH-Regulated Gene Is Involved in Iron-Sulfur Protein Metabolism and Effective Plant Symbiosis under Intrinsic Iron Limitation

In Sinorhizobium meliloti, RpoH-type sigma factors have a global impact on gene expression during heat shock and play an essential role in symbiosis with leguminous plants. Using mutational analysis of a set of genes showing highly RpoH-dependent expression during heat shock, we identified a gene indispensable for effective symbiosis. This gene, designated sufT, was located downstream of the sufBCDS homologs that specify the iron-sulfur (Fe/S) cluster assembly pathway. The identified transcription start site was preceded by an RpoH-dependent promoter consensus sequence. SufT was related to a conserved protein family of unknown molecular function, of which some members are involved in Fe/S cluster metabolism in diverse organisms. A sufT mutation decreased bacterial growth in both rich and minimal media, tolerance to stresses such as iron starvation, and activities of some Fe/S cluster-dependent enzymes. These results support the involvement of SufT in SUF (sulfur mobilization) system-mediated Fe/S protein metabolism. Furthermore, we isolated spontaneous pseudorevertants of the sufT mutant with partially recovered growth; each of them had a mutation in rirA This gene encodes a global iron regulator whose loss increases the intracellular iron content. Deletion of rirA in the original sufT mutant improved growth and restored Fe/S enzyme activities and effective symbiosis. These results suggest that enhanced iron availability compensates for the lack of SufT in the maintenance of Fe/S proteins.

IMPORTANCE

Although RpoH-type sigma factors of the RNA polymerase are present in diverse proteobacteria, their role as global regulators of protein homeostasis has been studied mainly in the enteric gammaproteobacterium Escherichia coli In the soil alphaproteobacterium Sinorhizobium meliloti, the rpoH mutations have a strong impact on symbiosis with leguminous plants. We found that sufT is a unique member of the S. meliloti RpoH regulon; sufT contributes to Fe/S protein metabolism and effective symbiosis under intrinsic iron limitation exerted by RirA, a global iron regulator. Our study provides insights into the RpoH regulon function in diverse proteobacteria adapted to particular ecological niches and into the mechanism of conserved Fe/S protein biogenesis.

The early events underlying genome evolution in a localized Sinorhizobium meliloti population

BACKGROUND

Population genetic analyses based on genome-wide sequencing data have been carried out for Sinorhizobium medicae and S. meliloti, two closely related bacterial species forming nitrogen-fixing symbioses with plants of the genus Medicago. However, genome coverage was low or the isolates had a broad geographic distribution, making it difficult to interpret the estimated diversity and to unravel the early events underlying population genetic variations and ecological differentiation.

RESULTS

Here, to gain insight into the early genome level variation and diversification within S. meliloti populations, we first used Illumina paired-end reads technology to sequence a new clone of S. meliloti strain GR4, a highly competitive strain for alfalfa nodulation. The Illumina data and the GR4 genome sequence previously obtained with 454 technology were used to generate a high-quality reference genome sequence. We then used Illumina technology to sequence the genomes of 13 S. meliloti isolates representative of the genomic variation within the GR4-type population, obtained from a single field site with a high degree of coverage. The genome sequences obtained were analyzed to determine nucleotide diversity, divergence times, polymorphism and genomic variation. Similar low levels of nucleotide diversity were observed for the chromosome, pSymB and pSymA replicons. The isolates displayed other types of variation, such as indels, recombination events, genomic island excision and the transposition of mobile elements.

CONCLUSIONS

Our results suggest that the GR4-type population has experienced a process of demographic expansion and behaves as a stable genotypic cluster of genome-wide similarity, with most of the genome following a clonal pattern of evolution. Although some of genetic variation detected within the GR4-type population is probably due to genetic drift, others might be important in diversification and environmental adaptation.

LDSS-P: an advanced algorithm to extract functional short motifs associated with coordinated gene expression

Identifying functional elements in promoter sequences is a major goal in computational and experimental genome biology. Here, we describe an algorithm, Local Distribution of Short Sequences for Prokaryotes (LDSS-P), to identify conserved short motifs located at specific positions in the promoters of co-expressed prokaryotic genes. As a test case, we applied this algorithm to a symbiotic nitrogen-fixing bacterium, Sinorhizobium meliloti. The LDSS-P profiles that overlap with the 5′ section of the extracytoplasmic function RNA polymerase sigma factor RpoE2 consensus sequences displayed a sharp peak between -34 and -32 from TSS positions. The corresponding genes overlap significantly with RpoE2 targets identified from previous experiments. We further identified several groups of genes that are co-regulated with characterized marker genes. Our data indicate that in S. meliloti, and possibly in other Rhizobiaceae species, the master cell cycle regulator CtrA may recognize an expanded motif (AACCAT), which is positionally shifted from the previously reported CtrA consensus sequence in Caulobacter crescentus. Bacterial one-hybrid experiments showed that base substitution in the expanded motif either increase or decrease the binding by CtrA. These results show the effectiveness of LDSS-P as a method to delineate functional promoter elements.

Evolution of a Sigma Factor: An All-In-One of Gene Duplication, Horizontal Gene Transfer, Purifying Selection, and Promoter Differentiation

Sigma factors are an essential part of bacterial gene regulation and have been extensively studied as far as their molecular mechanisms and protein structure are concerned. However, their molecular evolution, especially for the alternative sigma factors, is poorly understood. Here, we analyze the evolutionary forces that have shaped the rpoH sigma factors within the alphaproteobacteria. We found that an ancient duplication gave rise to two major groups of rpoH sigma factors and that after this event horizontal gene transfer (HGT) occurred in rpoH₁ group. We also noted that purifying selection has differentially affected distinct parts of the gene; singularly, the gene segment that encodes the region 4.2, which interacts with the −35 motif of the RpoH-dependent genes, has been under relaxed purifying selection. Furthermore, these two major groups are clearly differentiated from one another regarding their promoter selectivity, as rpoH₁ is under the transcriptional control of σ⁷⁰ and σ³², whereas rpoH₂ is under the transcriptional control of σ²⁴. Our results suggest a scenario in which HGT, gene loss, variable purifying selection and clear promoter specialization occurred after the ancestral duplication event. More generally, our study offers insights into the molecular evolution of alternative sigma factors and highlights the importance of analyzing not only the coding regions but also the promoter regions.

Genome-wide transcription start site mapping of Bradyrhizobium japonicum grown free-living or in symbiosis - a rich resource to identify new transcripts, proteins and to study gene regulation

BACKGROUND

Differential RNA-sequencing (dRNA-seq) is indispensable for determination of primary transcriptomes. However, using dRNA-seq data to map transcriptional start sites (TSSs) and promoters genome-wide is a bioinformatics challenge. We performed dRNA-seq of Bradyrhizobium japonicum USDA 110, the nitrogen-fixing symbiont of soybean, and developed algorithms to map TSSs and promoters.

RESULTS

A specialized machine learning procedure for TSS recognition allowed us to map 15,923 TSSs: 14,360 in free-living bacteria, 4329 in symbiosis with soybean and 2766 in both conditions. Further, we provide proteomic evidence for 4090 proteins, among them 107 proteins corresponding to new genes and 178 proteins with N-termini different from the existing annotation (72 and 109 of them with TSS support, respectively). Guided by proteomics evidence, previously identified TSSs and TSSs experimentally validated here, we assign a score threshold to flag 14 % of the mapped TSSs as a class of lower confidence. However, this class of lower confidence contains valid TSSs of low-abundant transcripts. Moreover, we developed a de novo algorithm to identify promoter motifs upstream of mapped TSSs, which is publicly available, and found motifs mainly used in symbiosis (similar to RpoN-dependent promoters) or under both conditions (similar to RpoD-dependent promoters). Mapped TSSs and putative promoters, proteomic evidence and updated gene annotation were combined into an annotation file.

CONCLUSIONS

The genome-wide TSS and promoter maps along with the extended genome annotation of B. japonicum represent a valuable resource for future systems biology studies and for detailed analyses of individual non-coding transcripts and ORFs. Our data will also provide new insights into bacterial gene regulation during the agriculturally important symbiosis between rhizobia and legumes.

Contributions of Sinorhizobium meliloti Transcriptional Regulator DksA to Bacterial Growth and Efficient Symbiosis with Medicago sativa

The stringent response, mediated by the (p)ppGpp synthetase RelA and the RNA polymerase-binding protein DksA, is triggered by limiting nutrient conditions. For some bacteria, it is involved in regulation of virulence. We investigated the role of two DksA-like proteins from the Gram-negative nitrogen-fixing symbiont Sinorhizobium meliloti in free-living culture and in interaction with its host plant Medicago sativa The two paralogs, encoded by the genes SMc00469 and SMc00049, differ in the constitution of two major domains required for function in canonical DksA: the DXXDXA motif at the tip of a coiled-coil domain and a zinc finger domain. Using mutant analyses of single, double, and triple deletions for SMc00469(designated dksA),SMc00049, and relA, we found that the ΔdksA mutant but not the ΔSMc00049 mutant showed impaired growth on minimal medium, reduced nodulation on the host plant, and lower nitrogen fixation activity in early nodules, while its nod gene expression was normal. The ΔrelA mutant showed severe pleiotropic phenotypes under all conditions tested. Only S. meliloti dksA complemented the metabolic defects of an Escherichia coli dksA mutant. Modifications of the DXXDXA motif in SMc00049 failed to establish DksA function. Our results imply a role for transcriptional regulator DksA in the S. meliloti-M. sativa symbiosis.

IMPORTANCE

The stringent response is a bacterial transcription regulation process triggered upon nutritional stress.Sinorhizobium meliloti, a soil bacterium establishing agriculturally important root nodule symbioses with legume plants, undergoes constant molecular adjustment during host interaction. Analyzing the components of the stringent response in this alphaproteobacterium helps understand molecular control regarding the development of plant interaction. Using mutant analyses, we describe how the lack of DksA influences symbiosis with Medicago sativa and show that a second paralogous S. meliloti protein cannot substitute for this missing function. This work contributes to the field by showing the similarities and differences of S. meliloti DksA-like proteins to orthologs from other species, adding information to the diversity of the stringent response regulatory system.

Expression of the Sinorhizobium meliloti small RNA gene mmgR is controlled by the nitrogen source

Small non-coding regulatory RNAs (sRNAs) are key players in post-transcriptional regulation of gene expression. Hundreds of sRNAs have been identified in Sinorhizobium meliloti, but their biological function remains unknown for most of them. In this study, we characterized the expression pattern of the gene encoding the 77-nt sRNA MmgR in S. meliloti strain 2011. A chromosomal transcriptional reporter fusion (PmmgR-gfp) showed that the mmgR promoter is active along different stages of the interaction with alfalfa roots. In pure cultures, PmmgR-gfp activity paralleled the sRNA abundance indicating that mmgR expression is primarily controlled at the level of transcriptional initiation. PmmgR-gfp activity was higher during growth in rhizobial defined medium (RDM) than in TY medium. Furthermore, PmmgR-gfp was induced at 60 min after shifting growing cells from TY to RDM medium, i.e. shorter than the cell doubling time. In defined RDM medium containing NO3 (-), both PmmgR-gfp and MmgR level were repressed by the addition of tryptone or single amino acids, suggesting that mmgR expression depends on the cellular nitrogen (N) status. In silico analysis failed to detect conserved motifs upstream the promoter RNA polymerase binding site, but revealed a strongly conserved motif centered at -28 that may be linked to the observed regulatory pattern by the N source.

Thermodynamic matchers for the construction of the cuckoo RNA family

RNA family models describe classes of functionally related, non-coding RNAs based on sequence and structure conservation. The most important method for modeling RNA families is the use of covariance models, which are stochastic models that serve in the discovery of yet unknown, homologous RNAs. However, the performance of covariance models in finding remote homologs is poor for RNA families with high sequence conservation, while for families with high structure but low sequence conservation, these models are difficult to built in the first place. A complementary approach to RNA family modeling involves the use of thermodynamic matchers. Thermodynamic matchers are RNA folding programs, based on the established thermodynamic model, but tailored to a specific structural motif. As thermodynamic matchers focus on structure and folding energy, they unfold their potential in discovering homologs, when high structure conservation is paired with low sequence conservation. In contrast to covariance models, construction of thermodynamic matchers does not require an input alignment, but requires human design decisions and experimentation, and hence, model construction is more laborious. Here we report a case study on an RNA family that was constructed by means of thermodynamic matchers. It starts from a set of known but structurally different members of the same RNA family. The consensus secondary structure of this family consists of 2 to 4 adjacent hairpins. Each hairpin loop carries the same motif, CCUCCUCCC, while the stems show high variability in their nucleotide content. The present study describes (1) a novel approach for the integration of the structurally varying family into a single RNA family model by means of the thermodynamic matcher methodology, and (2) provides the results of homology searches that were conducted with this model in a wide spectrum of bacterial species.

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

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

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

The complexity of bacterial transcriptomes

For eukaryotes there seems to be no doubt that differences on the trancriptomic level substantially contribute to the process of species diversification, whereas for bacteria this is thought to be less important. Recent years saw a significant increase in full transcriptome studies for bacteria, which provided deep insight into the architecture of bacterial transcriptomes. Most notably, it became evident that, in contrast to previous scientific consensus, bacterial transcriptomes are quite complex. There exist a large number of cis-antisense RNAs, non-coding RNAs, overlapping transcripts and RNA elements that regulate transcription, such as riboswitches. Furthermore, processing and degradation of RNA has gained interest, because it has a significant impact on the composition of the transcriptome. In this review, we summarize recent findings and put them into a broader context with respect to the complexity of bacterial transcriptomes and its putative biological meanings.

Evolution of Intra-specific Regulatory Networks in a Multipartite Bacterial Genome

Reconstruction of the regulatory network is an important step in understanding how organisms control the expression of gene products and therefore phenotypes. Recent studies have pointed out the importance of regulatory network plasticity in bacterial adaptation and evolution. The evolution of such networks within and outside the species boundary is however still obscure. Sinorhizobium meliloti is an ideal species for such study, having three large replicons, many genomes available and a significant knowledge of its transcription factors (TF). Each replicon has a specific functional and evolutionary mark; which might also emerge from the analysis of their regulatory signatures. Here we have studied the plasticity of the regulatory network within and outside the S. meliloti species, looking for the presence of 41 TFs binding motifs in 51 strains and 5 related rhizobial species. We have detected a preference of several TFs for one of the three replicons, and the function of regulated genes was found to be in accordance with the overall replicon functional signature: house-keeping functions for the chromosome, metabolism for the chromid, symbiosis for the megaplasmid. This therefore suggests a replicon-specific wiring of the regulatory network in the S. meliloti species. At the same time a significant part of the predicted regulatory network is shared between the chromosome and the chromid, thus adding an additional layer by which the chromid integrates itself in the core genome. Furthermore, the regulatory network distance was found to be correlated with both promoter regions and accessory genome evolution inside the species, indicating that both pangenome compartments are involved in the regulatory network evolution. We also observed that genes which are not included in the species regulatory network are more likely to belong to the accessory genome, indicating that regulatory interactions should also be considered to predict gene conservation in bacterial pangenomes.

The influence of transcriptional regulatory networks on the evolution of bacterial pangenomes has not yet been elucidated, even though the role of transcriptional regulation is widely recognized. Using the model symbiont Sinorhizobium meliloti we have predicted the regulatory targets of 41 transcription factors in 51 strains and 5 other rhizobial species, showing a correlation between regulon diversity and pangenome evolution, through upstream sequence diversity and accessory genome composition. We have also shown that genes not wired to the regulatory network are more likely to belong to the accessory genome, thus suggesting that inclusion in the regulatory circuits may be an indicator of gene conservation. We have also highlighted a series of transcription factors that preferentially regulate genes belonging to one of the three replicons of this species, indicating the presence of replicon-specific regulatory modules, with peculiar functional signatures. At the same time the chromid shares a significant part of the regulatory network with the chromosome, indicating an additional way by which this replicon integrates itself in the pangenome.

Transcriptomic Analysis of Sinorhizobium meliloti and Medicago truncatula Symbiosis Using Nitrogen Fixation-Deficient Nodules

The bacterium Sinorhizobium meliloti interacts symbiotically with legume plant hosts such as Medicago truncatula to form nitrogen-fixing root nodules. During symbiosis, plant and bacterial cells differentiate in a coordinated manner, resulting in specialized plant cells that contain nitrogen-fixing bacteroids. Both plant and bacterial genes are required at each developmental stage of symbiosis. We analyzed gene expression in nodules formed by wild-type bacteria on six plant mutants with defects in nitrogen fixation. We observed differential expression of 482 S. meliloti genes with functions in cell envelope homeostasis, cell division, stress response, energy metabolism, and nitrogen fixation. We simultaneously analyzed gene expression in M. truncatula and observed differential regulation of host processes that may trigger bacteroid differentiation and control bacterial infection. Our analyses of developmentally arrested plant mutants indicate that plants use distinct means to control bacterial infection during early and late symbiotic stages.

Genomic context drives transcription of insertion sequences in the bacterial endosymbiont Wolbachia wVulC

Transposable elements (TEs) are DNA pieces that are present in almost all the living world at variable genomic density. Due to their mobility and density, TEs are involved in a large array of genomic modifications. In eukaryotes, TE expression has been studied in detail in several species. In prokaryotes, studies of IS expression are generally linked to particular copies that induce a modification of neighboring gene expression. Here we investigated global patterns of IS transcription in the Alphaproteobacterial endosymbiont Wolbachia wVulC, using both RT-PCR and bioinformatic analyses. We detected several transcriptional promoters in all IS groups. Nevertheless, only one of the potentially functional IS groups possesses a promoter located upstream of the transposase gene, that could lead up to the production of a functional protein. We found that the majority of IS groups are expressed whatever their functional status. RT-PCR analyses indicate that the transcription of two IS groups lacking internal promoters upstream of the transposase start codon may be driven by the genomic environment. We confirmed this observation with the transcription analysis of individual copies of one IS group. These results suggest that the genomic environment is important for IS expression and it could explain, at least partly, copy number variability of the various IS groups present in the wVulC genome and, more generally, in bacterial genomes.

Cell Cycle Control by the Master Regulator CtrA in Sinorhizobium meliloti

In all domains of life, proper regulation of the cell cycle is critical to coordinate genome replication, segregation and cell division. In some groups of bacteria, e.g. Alphaproteobacteria, tight regulation of the cell cycle is also necessary for the morphological and functional differentiation of cells. Sinorhizobium meliloti is an alphaproteobacterium that forms an economically and ecologically important nitrogen-fixing symbiosis with specific legume hosts. During this symbiosis S. meliloti undergoes an elaborate cellular differentiation within host root cells. The differentiation of S. meliloti results in massive amplification of the genome, cell branching and/or elongation, and loss of reproductive capacity. In Caulobacter crescentus, cellular differentiation is tightly linked to the cell cycle via the activity of the master regulator CtrA, and recent research in S. meliloti suggests that CtrA might also be key to cellular differentiation during symbiosis. However, the regulatory circuit driving cell cycle progression in S. meliloti is not well characterized in both the free-living and symbiotic state. Here, we investigated the regulation and function of CtrA in S. meliloti. We demonstrated that depletion of CtrA cause cell elongation, branching and genome amplification, similar to that observed in nitrogen-fixing bacteroids. We also showed that the cell cycle regulated proteolytic degradation of CtrA is essential in S. meliloti, suggesting a possible mechanism of CtrA depletion in differentiated bacteroids. Using a combination of ChIP-Seq and gene expression microarray analysis we found that although S. meliloti CtrA regulates similar processes as C. crescentus CtrA, it does so through different target genes. For example, our data suggest that CtrA does not control the expression of the Fts complex to control the timing of cell division during the cell cycle, but instead it negatively regulates the septum-inhibiting Min system. Our findings provide valuable insight into how highly conserved genetic networks can evolve, possibly to fit the diverse lifestyles of different bacteria.

In order to propagate, all living cells must ensure that their genetic material is faithfully copied and properly partitioned into the daughter cells before division. These coordinated processes of DNA replication and cell division are termed the “cell cycle” and are controlled by a complex network of regulatory proteins in all organisms. In the class Alphaproteobacteria, the regulation of the cell cycle is closely linked to cellular differentiation processes that are vital for survival in the environment. In these bacteria, the cell cycle regulator CtrA is thought to serve as the primary link between the coordination of the cell cycle and cellular differentiation. The alphaproteobacterium, Sinorhizobium meliloti, an important model symbiont of alfalfa plants, undergoes a striking cellular differentiation that is vital to the formation of an efficient symbiosis dedicated to the conversion of atmospheric nitrogen to biologically available organic nitrogen. However, the link between cellular differentiation and cell cycle control in S. meliloti has not been made. In this study, we showed that S. meliloti cells without CtrA are similar to the symbiotic form. By the identification of the genes whose expression is directly and indirectly controlled by CtrA, we found that CtrA regulates vital cell cycle processes, including DNA replication and cell division, but through different genetic pathways than in other alphaproteobacteria. We importantly show that the levels of CtrA protein are governed by an essential cell cycle regulated proteolysis, which may also be an important mode of CtrA down-regulation during symbiosis to drive cellular differentiation.

An empirical strategy to detect bacterial transcript structure from directional RNA-seq transcriptome data

Background

As sequencing costs are being lowered continuously, RNA-seq has gradually been adopted as the first choice for comparative transcriptome studies with bacteria. Unlike microarrays, RNA-seq can directly detect cDNA derived from mRNA transcripts at a single nucleotide resolution. Not only does this allow researchers to determine the absolute expression level of genes, but it also conveys information about transcript structure. Few automatic software tools have yet been established to investigate large-scale RNA-seq data for bacterial transcript structure analysis.

Results

In this study, 54 directional RNA-seq libraries from Salmonella serovar Typhimurium (S. Typhimurium) 14028s were examined for potential relationships between read mapping patterns and transcript structure. We developed an empirical method, combined with statistical tests, to automatically detect key transcript features, including transcriptional start sites (TSSs), transcriptional termination sites (TTSs) and operon organization. Using our method, we obtained 2,764 TSSs and 1,467 TTSs for 1331 and 844 different genes, respectively. Identification of TSSs facilitated further discrimination of 215 putative sigma 38 regulons and 863 potential sigma 70 regulons. Combining the TSSs and TTSs with intergenic distance and co-expression information, we comprehensively annotated the operon organization in S. Typhimurium 14028s.

Conclusions

Our results show that directional RNA-seq can be used to detect transcriptional borders at an acceptable resolution of ±10-20 nucleotides. Technical limitations of the RNA-seq procedure may prevent single nucleotide resolution. The automatic transcript border detection methods, statistical models and operon organization pipeline that we have described could be widely applied to RNA-seq studies in other bacteria. Furthermore, the TSSs, TTSs, operons, promoters and unstranslated regions that we have defined for S. Typhimurium 14028s may constitute valuable resources that can be used for comparative analyses with other Salmonella serotypes.

Electronic supplementary material

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

Small RNA Deep-Sequencing Analyses Reveal a New Regulator of Virulence in Agrobacterium fabrum C58

Novel ways of regulating Ti plasmid functions were investigated by studying small RNAs (sRNAs) that are known to act as posttranscriptional regulators in plant pathogenic bacteria. sRNA-seq analyses of Agrobacterium fabrum C58 allowed us to identify 1,108 small transcripts expressed in several growth conditions that could be sRNAs. A quarter of them were confirmed by bioinformatics or by biological experiments. Antisense RNAs represent 24% of the candidates and they are over-represented on the pTi (with 62% of pTi sRNAs), suggesting differences in the regulatory mechanisms between the essential and accessory replicons. Moreover, a large number of these pTi antisense RNAs are transcribed opposite to those genes involved in virulence. Others are 5'- and 3'-untranslated region RNAs and trans-encoded RNAs. We have validated, by rapid amplification of cDNA ends polymerase chain reaction, the transcription of 14 trans-encoded RNAs, among which RNA1111 is expressed from the pTiC58. Its deletion decreased the aggressiveness of A. fabrum C58 on tomatoes, tobaccos, and kalanchoe, suggesting that this sRNA activates virulence. The identification of its putative target mRNAs (6b gene, virC2, virD3, and traA) suggests that this sRNA may coordinate two of the major pTi functions, the infection of plants and its dissemination among bacteria.

A stress-induced small RNA modulates alpha-rhizobial cell cycle progression

Mechanisms adjusting replication initiation and cell cycle progression in response to environmental conditions are crucial for microbial survival. Functional characterization of the trans-encoded small non-coding RNA (trans-sRNA) EcpR1 in the plant-symbiotic alpha-proteobacterium Sinorhizobium meliloti revealed a role of this class of riboregulators in modulation of cell cycle regulation. EcpR1 is broadly conserved in at least five families of the Rhizobiales and is predicted to form a stable structure with two defined stem-loop domains. In S. meliloti, this trans-sRNA is encoded downstream of the divK-pleD operon. ecpR1 belongs to the stringent response regulon, and its expression was induced by various stress factors and in stationary phase. Induced EcpR1 overproduction led to cell elongation and increased DNA content, while deletion of ecpR1 resulted in reduced competitiveness. Computationally predicted EcpR1 targets were enriched with cell cycle-related mRNAs. Post-transcriptional repression of the cell cycle key regulatory genes gcrA and dnaA mediated by mRNA base-pairing with the strongly conserved loop 1 of EcpR1 was experimentally confirmed by two-plasmid differential gene expression assays and compensatory changes in sRNA and mRNA. Evidence is presented for EcpR1 promoting RNase E-dependent degradation of the dnaA mRNA. We propose that EcpR1 contributes to modulation of cell cycle regulation under detrimental conditions.

Microorganisms frequently encounter adverse conditions unfavorable for cell proliferation. They have evolved diverse mechanisms, including transcriptional control and targeted protein degradation, to adjust cell cycle progression in response to environmental cues. Non-coding RNAs are widespread regulators of various cellular processes in all domains of life. In prokaryotes, trans-encoded small non-coding RNAs (trans-sRNAs) contribute to a rapid cellular response to changing environments, but so far have not been directly related to cell cycle regulation. Here, we report the first example of a trans-sRNA (EcpR1) with two experimentally confirmed targets in the core of cell cycle regulation and demonstrate that in the plant-symbiotic alpha-proteobacterium Sinorhizobium meliloti the regulatory mechanism involves base-pairing of this sRNA with the dnaA and gcrA mRNAs. Most trans-sRNAs are restricted to closely related species, but the stress-induced EcpR1 is broadly conserved in the order of Rhizobiales suggesting an evolutionary advantage conferred by ecpR1. It broadens the functional diversity of prokaryotic sRNAs and adds a new regulatory level to the mechanisms that contribute to interlinking stress responses with the cell cycle machinery.

Transcriptomic profiling of Yersinia pseudotuberculosis reveals reprogramming of the Crp regulon by temperature and uncovers Crp as a master regulator of small RNAs

One hallmark of pathogenic yersiniae is their ability to rapidly adjust their life-style and pathogenesis upon host entry. In order to capture the range, magnitude and complexity of the underlying gene control mechanisms we used comparative RNA-seq-based transcriptomic profiling of the enteric pathogen Y. pseudotuberculosis under environmental and infection-relevant conditions. We identified 1151 individual transcription start sites, multiple riboswitch-like RNA elements, and a global set of antisense RNAs and previously unrecognized trans-acting RNAs. Taking advantage of these data, we revealed a temperature-induced and growth phase-dependent reprogramming of a large set of catabolic/energy production genes and uncovered the existence of a thermo-regulated ‘acetate switch’, which appear to prime the bacteria for growth in the digestive tract. To elucidate the regulatory architecture linking nutritional status to virulence we also refined the CRP regulon. We identified a massive remodelling of the CRP-controlled network in response to temperature and discovered CRP as a transcriptional master regulator of numerous conserved and newly identified non-coding RNAs which participate in this process. This finding highlights a novel level of complexity of the regulatory network in which the concerted action of transcriptional regulators and multiple non-coding RNAs under control of CRP adjusts the control of Yersinia fitness and virulence to the requirements of their environmental and virulent life-styles.

Many bacterial pathogens cycle between environmental sources and mammalian hosts. Adaptation to the different natural habitats and host niches is achieved through complex regulatory networks which adjust synthesis of the large repertoire of crucial virulence factors and fitness determinants. To uncover underlying control circuits, we determined the first in-depth single-nucleotide resolution transcriptome of Yersinia. This revealed important novel genetic information, such as global locations of transcriptional start sites, non-coding RNAs, potential riboswitches and provided a set of virulence-relevant expression profiles, which constitute a valuable tool for the research community. The analysis further uncovered a temperature-induced global reprogramming of central metabolic functions, likely to support intestinal colonization of the pathogen. This is accompanied by a major reorganization of the CRP regulon, which involves a multitude of regulatory RNAs. The primary consequence is a fine-tuned, coordinated control of metabolism and virulence through a plethora of environmentally controlled regulatory RNAs allowing rapid adaptation and high flexibility during life-style changes.

The Sinorhizobium meliloti SyrM regulon: effects on global gene expression are mediated by syrA and nodD3

In Sinorhizobium meliloti, three NodD transcriptional regulators activate bacterial nodulation (nod) gene expression. NodD1 and NodD2 require plant compounds to activate nod genes. The NodD3 protein does not require exogenous compounds to activate nod gene expression; instead, another transcriptional regulator, SyrM, activates nodD3 expression. In addition, NodD3 can activate syrM expression. SyrM also activates expression of another gene, syrA, which when overexpressed causes a dramatic increase in exopolysaccharide production. In a previous study, we identified more than 200 genes with altered expression in a strain overexpressing nodD3. In this work, we define the transcriptomes of strains overexpressing syrM or syrA. The syrM, nodD3, and syrA overexpression transcriptomes share similar gene expression changes; analyses imply that nodD3 and syrA are the only targets directly activated by SyrM. We propose that most of the gene expression changes observed when nodD3 is overexpressed are due to NodD3 activation of syrM expression, which in turn stimulates SyrM activation of syrA expression. The subsequent increase in SyrA abundance results in broad changes in gene expression, most likely mediated by the ChvI-ExoS-ExoR regulatory circuit.

IMPORTANCE

Symbioses with bacteria are prevalent across the animal and plant kingdoms. Our system of study, the rhizobium-legume symbiosis (Sinorhizobium meliloti and Medicago spp.), involves specific host-microbe signaling, differentiation in both partners, and metabolic exchange of bacterial fixed nitrogen for host photosynthate. During this complex developmental process, both bacteria and plants undergo profound changes in gene expression. The S. meliloti SyrM-NodD3-SyrA and ChvI-ExoS-ExoR regulatory circuits affect gene expression and are important for optimal symbiosis. In this study, we defined the transcriptomes of S. meliloti overexpressing SyrM or SyrA. In addition to identifying new targets of the SyrM-NodD3-SyrA regulatory circuit, our work further suggests how it is linked to the ChvI-ExoS-ExoR regulatory circuit.

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

Deinococcus deserti is a desiccation- and radiation-tolerant desert bacterium. Differential RNA sequencing (RNA-seq) was performed to explore the specificities of its transcriptome. Strikingly, for 1,174 (60%) mRNAs, the transcription start site was found exactly at (916 cases, 47%) or very close to the translation initiation codon AUG or GUG. Such proportion of leaderless mRNAs, which may resemble ancestral mRNAs, is unprecedented for a bacterial species. Proteomics showed that leaderless mRNAs are efficiently translated in D. deserti. Interestingly, we also found 173 additional transcripts with a 5′-AUG or 5′-GUG that would make them competent for ribosome binding and translation into novel small polypeptides. Fourteen of these are predicted to be leader peptides involved in transcription attenuation. Another 30 correlated with new gene predictions and/or showed conservation with annotated and nonannotated genes in other Deinococcus species, and five of these novel polypeptides were indeed detected by mass spectrometry. The data also allowed reannotation of the start codon position of 257 genes, including several DNA repair genes. Moreover, several novel highly radiation-induced genes were found, and their potential roles are discussed. On the basis of our RNA-seq and proteogenomics data, we propose that translation of many of the novel leaderless transcripts, which may have resulted from single-nucleotide changes and maintained by selective pressure, provides a new explanation for the generation of a cellular pool of small peptides important for protection of proteins against oxidation and thus for radiation/desiccation tolerance and adaptation to harsh environmental conditions.

Genetic redundancy is prevalent within the 6.7 Mb Sinorhizobium meliloti genome

Biological pathways are frequently identified via a genetic loss-of-function approach. While this approach has proven to be powerful, it is imperfect as illustrated by well-studied pathways continuing to have missing steps. One potential limiting factor is the masking of phenotypes through genetic redundancy. The prevalence of genetic redundancy in bacterial species has received little attention, although isolated examples of functionally redundant gene pairs exist. Here, we made use of a strain of Sinorhizobium meliloti whose genome was reduced by 45 % through the complete removal of a megaplasmid and a chromid (3 Mb of the 6.7 Mb genome was removed) to begin quantifying the level of genetic redundancy within a large bacterial genome. A mutagenesis of the strain with the reduced genome identified a set of transposon insertions precluding growth of this strain on minimal medium. Transfer of these mutations to the wild-type background revealed that 10-15 % of these chromosomal mutations were located within duplicated genes, as they did not prevent growth of cells with the full genome. The functionally redundant genes were involved in a variety of metabolic pathways, including central carbon metabolism, transport, and amino acid biosynthesis. These results indicate that genetic redundancy may be prevalent within large bacterial genomes. Failing to account for redundantly encoded functions in loss-of-function studies will impair our understanding of a broad range of biological processes and limit our ability to use synthetic biology in the construction of designer cell factories.

The global regulatory architecture of transcription during the Caulobacter cell cycle

Each Caulobacter cell cycle involves differentiation and an asymmetric cell division driven by a cyclical regulatory circuit comprised of four transcription factors (TFs) and a DNA methyltransferase. Using a modified global 5′ RACE protocol, we globally mapped transcription start sites (TSSs) at base-pair resolution, measured their transcription levels at multiple times in the cell cycle, and identified their transcription factor binding sites. Out of 2726 TSSs, 586 were shown to be cell cycle-regulated and we identified 529 binding sites for the cell cycle master regulators. Twenty-three percent of the cell cycle-regulated promoters were found to be under the combinatorial control of two or more of the global regulators. Previously unknown features of the core cell cycle circuit were identified, including 107 antisense TSSs which exhibit cell cycle-control, and 241 genes with multiple TSSs whose transcription levels often exhibited different cell cycle timing. Cumulatively, this study uncovered novel new layers of transcriptional regulation mediating the bacterial cell cycle.

The generation of diverse cell types occurs through two fundamental processes; asymmetric cell division and cell differentiation. Cells progress through these developmental changes guided by complex and layered genetic programs that lead to differential expression of the genome. To explore how a genetic program directs cell cycle progression, we examined the global activity of promoters at distinct stages of the cell cycle of the bacterium Caulobacter crescentus, that undergoes cellular differentiation and divides asymmetrically at each cell division. We found that approximately 21% of transcription start sites are cell cycle-regulated, driving the transcription of both mRNAs and non-coding and antisense RNAs. In addition, 102 cell cycle-regulated genes are transcribed from multiple promoters, allowing multiple regulatory inputs to control the logic of gene activation. We found combinatorial control by the five master transcription regulators that provide the core regulation for the genetic circuitry controlling the cell cycle. Much of this combinatorial control appears to be directed at refinement of temporal expression of various genes over the cell cycle, and at tighter control of asymmetric gene expression between the swarmer and stalked daughter cells.

Complex and extensive post-transcriptional regulation revealed by integrative proteomic and transcriptomic analysis of metabolite stress response in Clostridium acetobutylicum

BACKGROUND

Clostridium acetobutylicum is a model organism for both clostridial biology and solvent production. The organism is exposed to its own toxic metabolites butyrate and butanol, which trigger an adaptive stress response. Integrative analysis of proteomic and RNAseq data may provide novel insights into post-transcriptional regulation.

RESULTS

The identified iTRAQ-based quantitative stress proteome is made up of 616 proteins with a 15 % genome coverage. The differentially expressed proteome correlated poorly with the corresponding differential RNAseq transcriptome. Up to 31 % of the differentially expressed proteins under stress displayed patterns opposite to those of the transcriptome, thus suggesting significant post-transcriptional regulation. The differential proteome of the translation machinery suggests that cells employ a different subset of ribosomal proteins under stress. Several highly upregulated proteins but with low mRNA levels possessed mRNAs with long 5'UTRs and strong RBS scores, thus supporting the argument that regulatory elements on the long 5'UTRs control their translation. For example, the oxidative stress response rubrerythrin was upregulated only at the protein level up to 40-fold without significant mRNA changes. We also identified many leaderless transcripts, several displaying different transcriptional start sites, thus suggesting mRNA-trimming mechanisms under stress. Downregulation of Rho and partner proteins pointed to changes in transcriptional elongation and termination under stress.

CONCLUSIONS

The integrative proteomic-transcriptomic analysis demonstrated complex expression patterns of a large fraction of the proteome. Such patterns could not have been detected with one or the other omic analyses. Our analysis proposes the involvement of specific molecular mechanisms of post-transcriptional regulation to explain the observed complex stress response.

Analysis of a taurine-dependent promoter in Sinorhizobium meliloti that offers tight modulation of gene expression

Background

Genetic models have been developed in divergent branches of the class Alphaproteobacteria to help answer a wide spectrum of questions regarding bacterial physiology. For example, Sinorhizobium meliloti serves as a useful representative for investigating rhizobia-plant symbiosis and nitrogen fixation, Caulobacter crescentus for studying cell cycle regulation and organelle biogenesis, and Zymomonas mobilis for assessing the potentials of metabolic engineering and biofuel production. A tightly regulated promoter that enables titratable expression of a cloned gene in these different models is highly desirable, as it can facilitate observation of phenotypes that would otherwise be obfuscated by leaky expression.

Results

We compared the functionality of four promoter regions in S. meliloti (P_araA, P_tauA, P_rhaR, and P_melA) by constructing strains carrying fusions to the uidA reporter in their genomes and measuring beta-glucuronidase activities when they were induced by arabinose, taurine, rhamnose, or melibiose. P_tauA was chosen for further study because it, and, to a lesser extent, P_melA, exhibited characteristics suitable for efficient modulation of gene expression. The levels of expression from P_tauA depended on the concentrations of taurine, in both complex and defined media, in S. meliloti as well as C. crescentus and Z. mobilis. Moreover, our analysis indicated that TauR, TauC, and TauY are each necessary for taurine catabolism and substantiated their designated roles as a transcriptional activator, the permease component of an ABC transporter, and a major subunit of the taurine dehydrogenase, respectively. Finally, we demonstrated that P_tauA can be used to deplete essential cellular factors in S. meliloti, such as the PleC histidine kinase and TatB, a component of the twin-arginine transport machinery.

Conclusions

The P_tauA promoter of S. meliloti can control gene expression with a relatively inexpensive and permeable inducer, taurine, in diverse alpha-proteobacteria. Regulated expression of the same gene in different hosts can be achieved by placing both tauR and P_tauA on appropriate vectors, thus facilitating inspection of conservation of gene function across species.

Electronic supplementary material

The online version of this article (doi:10.1186/s12866-014-0295-2) contains supplementary material, which is available to authorized users.