Post-transcriptional regulons coordinate the expression of groups of genes in eukaryotic cells, yet relatively few have been characterized. Parasitic trypanosomatids are particularly good models for studies on such mechanisms because they exhibit almost exclusive polycistronic, and unregulated, transcription. Here, we identify the Trypanosoma brucei ZC3H39/40 RNA-binding proteins as regulators of the respiratome; the mitochondrial electron transport chain (complexes I-IV) and the FoF1-ATP synthase (complex V). A high-throughput RNAi screen initially implicated both ZC3H proteins in variant surface glycoprotein (VSG) gene silencing. This link was confirmed and both proteins were shown to form a cytoplasmic ZC3H39/40 complex. Transcriptome and mRNA-interactome analyses indicated that the impact on VSG silencing was indirect, while the ZC3H39/40 complex specifically bound and stabilized transcripts encoding respiratome-complexes. Quantitative proteomic analyses revealed specific positive control of >20 components from complexes I, II and V. Our findings establish a link between the mitochondrial respiratome and VSG gene silencing in bloodstream form T. brucei. They also reveal a major respiratome regulon controlled by the conserved trypanosomatid ZC3H39/40 RNA-binding proteins.

Fine-tuning of inflammatory responses by microRNAs (miRNAs) is complex, as they can both enhance and repress expression of pro-inflammatory mediators. In this study, we investigate inflammatory responses following global miRNA depletion, to better define the overall contribution of miRNAs to inflammation. We demonstrate that miRNAs positively regulate Toll-like receptor signaling using inducible Dicer1 deletion and global miRNA depletion. We establish an important contribution of miR-19b in this effect, which potentiates nuclear factor-κB (NF-κB) activity in human and mouse cells. Positive regulation of NF-κB signaling by miR-19b involves the coordinated suppression of a regulon of negative regulators of NF-κB signaling (including A20/Tnfaip3, Rnf11, Fbxl11/Kdm2a and Zbtb16). Transfection of miR-19b mimics exacerbated the inflammatory activation of rheumatoid arthritis primary fibroblast-like synoviocytes, demonstrating its physiological importance in the pathology of this disease. This study constitutes, to our knowledge, the first description of a miR-19 regulon that controls NF-κB signaling, and suggests that targeting this miRNA and linked family members could regulate the activity of NF-κB signaling in inflammation.

Immune system gene regulation perturbation has been found to be a major cause of the development of various types of cancer. Numbers of mechanisms contribute to gene expression regulation, thus, systematically identification of potential regulons of immune-related pathways is critical to cancer immunotherapy. Here, we comprehensively chart the landscape of transcription factors, microRNAs, RNA binding proteins and long noncoding RNAs regulation in 17 immune-related pathways across 33 cancers. The potential immunology regulons are likely to exhibit higher expressions in immune cells, show expression perturbations in cancer, and are significantly correlated with immune cell infiltrations. We also identify a panel of clinically relevant immunology regulons across cancers. Moreover, the regulon atlas of immune-related pathways helps prioritizing cancer-related genes (i.e. ETV7, miR-146a-5p, ZFP36 and HCP5). We further identified two molecular subtypes of glioma (cold and hot tumour phenotypes), which were characterized by differences in immune cell infiltrations, expression of checkpoints, and prognosis. Finally, we developed a user-friendly resource, ImmReg (http://bio-bigdata.hrbmu.edu.cn/ImmReg/), with multiple modules to visualize, browse, and download immunology regulation. Our study provides a comprehensive landscape of immunology regulons, which will shed light on future development of RNA-based cancer immunotherapies.

The plant-specific, B3 domain-containing transcription factor ABSCISIC ACID INSENSITIVE3 (ABI3) is an essential component of the regulatory network controlling the development and maturation of the Arabidopsis thaliana seed. Genome-wide chromatin immunoprecipitation (ChIP-chip), transcriptome analysis, quantitative reverse transcriptase-polymerase chain reaction and a transient promoter activation assay have been combined to identify a set of 98 ABI3 target genes. Most of these presumptive ABI3 targets require the presence of abscisic acid for their activation and are specifically expressed during seed maturation. ABI3 target promoters are enriched for G-box-like and RY-like elements. The general occurrence of these cis motifs in non-ABI3 target promoters suggests the existence of as yet unidentified regulatory signals, some of which may be associated with epigenetic control. Several members of the ABI3 regulon are also regulated by other transcription factors, including the seed-specific, B3 domain-containing FUS3 and LEC2. The data strengthen and extend the notion that ABI3 is essential for the protection of embryonic structures from desiccation and raise pertinent questions regarding the specificity of promoter recognition.

HrdB in streptomycetes is a principal sigma factor whose deletion is lethal. This is also the reason why its regulon has not been investigated so far. To overcome experimental obstacles, for investigating the HrdB regulon, we constructed a strain whose HrdB protein was tagged by an HA epitope. ChIP-seq experiment, done in 3 repeats, identified 2137 protein-coding genes organized in 337 operons, 75 small RNAs, 62 tRNAs, 6 rRNAs and 3 miscellaneous RNAs. Subsequent kinetic modeling of regulation of protein-coding genes with HrdB alone and with a complex of HrdB and a transcriptional cofactor RbpA, using gene expression time series, identified 1694 genes that were under their direct control. When using the HrdB-RbpA complex in the model, an increase of the model fidelity was found for 322 genes. Functional analysis revealed that HrdB controls the majority of gene groups essential for the primary metabolism and the vegetative growth. Particularly, almost all ribosomal protein-coding genes were found in the HrdB regulon. Analysis of promoter binding sites revealed binding motif at the -10 region and suggested the possible role of mono- or di-nucleotides upstream of the -10 element.

The genome-wide mapping of gene-regulatory motifs remains a major goal that will facilitate the modelling of gene-regulatory networks and their evolution. The repressor element 1 is a long, conserved transcription factor-binding site which recruits the transcriptional repressor REST to numerous neuron-specific target genes. REST plays important roles in multiple biological processes and disease states. To map RE1 sites and target genes, we created a position specific scoring matrix representing the RE1 and used it to search the human and mouse genomes. We identified 1301 and 997 RE1s inhuman and mouse genomes, respectively, of which >40% are novel. By employing an ontological analysis we show that REST target genes are significantly enriched in a number of functional classes. Taking the novel REST target gene CACNA1A as an experimental model, we show that it can be regulated by multiple RE1s of different binding affinities, which are only partially conserved between human and mouse. A novel BLAST methodology indicated that many RE1s belong to closely related families. Most of these sequences are associated with transposable elements, leading us to propose that transposon-mediated duplication and insertion of RE1s has led to the acquisition of novel target genes by REST during evolution.

RNA-binding proteins (RBPs) are essential to fine-tune gene expression. RBPs containing the cold-shock domain are RNA chaperones that have been extensively studied. However, the RNA targets and specific functions for many of them remain elusive. Here, combining comparative proteomics and RBP-immunoprecipitation-microarray profiling, we have determined the regulon of the RNA chaperone CspA of Staphylococcus aureus. Functional analysis revealed that proteins involved in carbohydrate and ribonucleotide metabolism, stress response and virulence gene expression were affected by cspA deletion. Stress-associated phenotypes such as increased bacterial aggregation and diminished resistance to oxidative-stress stood out. Integration of the proteome and targetome showed that CspA post-transcriptionally modulates both positively and negatively the expression of its targets, denoting additional functions to the previously proposed translation enhancement. One of these repressed targets was its own mRNA, indicating the presence of a negative post-transcriptional feedback loop. CspA bound the 5'UTR of its own mRNA disrupting a hairpin, which was previously described as an RNase III target. Thus, deletion of the cspA 5'UTR abrogated mRNA processing and auto-regulation. We propose that CspA interacts through a U-rich motif, which is located at the RNase III cleavage site, portraying CspA as a putative RNase III-antagonist.

RcsB is a transcriptional regulator that controls expression of numerous genes in enteric bacteria. RcsB accomplishes this role alone or in combination with auxiliary transcriptional factors independently or dependently of phosphorylation. To understand the mechanisms by which RcsB regulates such large number of genes, we performed structural studies as well as in vitro and in vivo functional studies with different RcsB variants. Our structural data reveal that RcsB binds promoters of target genes such as rprA and flhDC in a dimeric active conformation. In this state, the RcsB homodimer docks the DNA-binding domains into the major groove of the DNA, facilitating an initial weak read-out of the target sequence. Interestingly, comparative structural analyses also show that DNA binding may stabilize an active conformation in unphosphorylated RcsB. Furthermore, RNAseq performed in strains expressing wild-type or several RcsB variants provided new insights into the contribution of phosphorylation to gene regulation and assign a potential role of RcsB in controlling iron metabolism. Finally, we delimited the RcsB box for homodimeric active binding to DNA as the sequence TN(G/A)GAN4TC(T/C)NA. This RcsB box was found in promoter, intergenic and intragenic regions, facilitating both increased or decreased gene transcription.

A comparative genomic approach was used to reconstruct transcriptional regulation of NAD biosynthesis in bacteria containing orthologs of Bacillus subtilis gene yrxA, a previously identified niacin-responsive repressor of NAD de novo synthesis. Members of YrxA family (re-named here NiaR) are broadly conserved in the Bacillus/Clostridium group and in the deeply branching Fusobacteria and Thermotogales lineages. We analyzed upstream regions of genes associated with NAD biosynthesis to identify candidate NiaR-binding DNA motifs and assess the NiaR regulon content in these species. Representatives of the two distinct types of candidate NiaR-binding sites, characteristic of the Firmicutes and Thermotogales, were verified by an electrophoretic mobility shift assay. In addition to transcriptional control of the nadABC genes, the NiaR regulon in some species extends to niacin salvage (the pncAB genes) and includes uncharacterized membrane proteins possibly involved in niacin transport. The involvement in niacin uptake proposed for one of these proteins (re-named NiaP), encoded by the B. subtilis gene yceI, was experimentally verified. In addition to bacteria, members of the NiaP family are conserved in multicellular eukaryotes, including human, pointing to possible NaiP involvement in niacin utilization in these organisms. Overall, the analysis of the NiaR and NrtR regulons (described in the accompanying paper) revealed mechanisms of transcriptional regulation of NAD metabolism in nearly a hundred diverse bacteria.

While global transcription factors (TFs) have been studied extensively in Escherichia coli model strains, conservation and diversity in TF regulation between strains is still unknown. Here we use a combination of ChIP-exo-to define ferric uptake regulator (Fur) binding sites-and differential gene expression-to define the Fur regulon in nine E. coli strains. We then define a pan-regulon consisting of 469 target genes that includes all Fur target genes in all nine strains. The pan-regulon is then divided into the core regulon (target genes found in all the strains, n = 36), the accessory regulon (target found in two to eight strains, n = 158) and the unique regulon (target genes found in one strain, n = 275). Thus, there is a small set of Fur regulated genes common to all nine strains, but a large number of regulatory targets unique to a particular strain. Many of the unique regulatory targets are genes unique to that strain. This first-established pan-regulon reveals a common core of conserved regulatory targets and significant diversity in transcriptional regulation amongst E. coli strains, reflecting diverse niche specification and strain history.

Flexible adaptation to environmental stress is vital for bacteria. An energy-efficient post-transcriptional stress response mechanism in Escherichia coli is governed by the toxin MazF. After stress-induced activation the endoribonuclease MazF processes a distinct subset of transcripts as well as the 16S ribosomal RNA in the context of mature ribosomes. As these 'stress-ribosomes' are specific for the MazF-processed mRNAs, the translational program is changed. To identify this 'MazF-regulon' we employed Poly-seq (polysome fractionation coupled with RNA-seq analysis) and analyzed alterations introduced into the transcriptome and translatome after mazF overexpression. Unexpectedly, our results reveal that the corresponding protein products are involved in all cellular processes and do not particularly contribute to the general stress response. Moreover, our findings suggest that translational reprogramming serves as a fast-track reaction to harsh stress and highlight the so far underestimated significance of selective translation as a global regulatory mechanism in gene expression. Considering the reported implication of toxin-antitoxin (TA) systems in persistence, our results indicate that MazF acts as a prime effector during harsh stress that potentially introduces translational heterogeneity within a bacterial population thereby stimulating persister cell formation.

In the absence of arabinose, the AraR transcription factor represses the expression of genes involved in the utilization of arabinose, xylose and galactose in Bacillus subtilis. AraR exhibits a chimeric organization: the N-terminal DNA-binding region belongs to the GntR family and the C-terminal effector-binding domain is homologous to the GalR/LacI family. Here, the AraR-DNA-binding interactions were characterized in vivo and in vitro. The effect of residue substitutions in the AraR N-terminal domain and of base-pair exchanges into an AraR-DNA-binding operator site were examined by assaying for AraR-mediated regulatory activity in vivo and DNA-binding activity in vitro. The results showed that residues K4, R45 and Q61, located in or near the winged-helix DNA-binding motif, were the most critical amino acids required for AraR function. In addition, the analysis of the various mutations in an AraR palindromic operator sequence indicated that bases G9, A11 and T16 are crucial for AraR binding. Moreover, an AraR mutant M34T was isolated that partially suppressed the effect of mutations in the regulatory cis-elements. Together, these findings extend the knowledge on the nature of AraR nucleoprotein complexes and provide insight into the mechanism that underlies the mode of action of AraR and its orthologues.

The Mycoplasma genitalium MG428 protein shows homology to members of the sigma-70 family of sigma factors. Herein, we found that MG428 activates transcription of recA, ruvA and ruvB as well as several genes with unknown function. Deletion of MG_428 or some of the up-regulated unknown genes led to severe recombination defects. Single cell analyses revealed that activation of the MG428-regulon is a rare event under laboratory growth conditions. A conserved sequence with sigma-70 promoter architecture (TTGTCA-N(18/19)-ATTWAT) was identified in the upstream region of all of the MG428-regulated genes or operons. Primer extension analyses demonstrated that transcription initiates immediately downstream of this sigma70-type promoter in a MG428-dependent manner. Furthermore, mutagenesis of the conserved -10 and -35 elements corroborated the requirement of these regions for promoter function. Therefore, a new mycoplasma promoter directs transcription of a unique recombination regulon. Additionally, MG428 was found to interact with the RNAP core enzyme, reinforcing the predicted role of this protein as an alternative sigma factor. Finally, our results indicate that MG428 contributes to the generation of genetic diversity in this model organism. Since recombination is an important mechanism to generate antigenic variation, MG428 emerges as a novel factor contributing to M. genitalium virulence.

The LysR-family regulator MexT modulates the expression of the MexEF-OprN efflux system in the human pathogen Pseudomonas aeruginosa. Recently, we demonstrated that MexT regulates certain virulence phenotypes, including the type-three secretion system and early attachment independent of its role in regulating MexEF-OprN. In this study, transcriptome profiling was utilized to investigate the global nature of MexT regulation in P. aeruginosa PAO1 and an isogenic mexEF mutant. Twelve genes of unknown function were highly induced by overexpressing MexT independent of MexEF-OprN. A well-conserved DNA motif was identified in the upstream regulatory region of nine of these genes and upstream of mexE. Reporter fusion analysis demonstrated that the expression of the genes was significantly induced by MexT in P. aeruginosa and a heterogenous Escherichia coli strain and that the conserved sequence was required for this induction. The conserved DNA motif was further characterized as the MexT binding site by site-directed mutagenesis and electrophoretic mobility shift assays. Genes containing this conserved regulatory sequence were identified across other Pseudomonas species, and their expression was activated by MexT. Thus, a novel regulon directly modulated by MexT, that includes but is independent of mexEF-oprN, has been identified.

Bacterial flagellar motility and chemotaxis help cells to reach the most favorable environments and to successfully compete with other micro-organisms in response to external stimuli. Escherichia coli is a motile gram-negative bacterium, and the flagellar regulon in E. coli is controlled by a master regulator FlhDC as well as a second regulator, flagellum-specific sigma factor, sigma(F). To define the physiological role of these two regulators, we carried out transcription profiling experiments to identify, on a genome-wide basis, genes under the control of these two regulators. In addition, the synchronized pattern of increasing CRP activity causing increasing FlhDC expression with decreasing carbon source quality, together with the apparent coupling of motility activity with the activation of motility and chemotaxis genes in poor quality carbon sources, highlights the importance of CRP activation in allowing E. coli to devote progressively more of its limited reserves to search out better conditions. In adaptation to a variety of carbon sources, the motile bacteria carry out tactical responses by increasing flagellar operation but restricting costly flagellar synthesis, indicating its capability of strategically using the precious energy in nutrient-poor environments for maximizing survival.

Bacteria core RNA polymerase (RNAP) must associate with a sigma factor to recognize promoter sequences. Promoters recognized by the sigma(54) (or sigma(N)) associated RNA polymerase are unique in having conserved positions around -24 and -12 nucleotides upstream from the transcriptional start site. Using DNA microarrays representing the entire Escherichia coli genome and promoter validation approaches, we identify 40 in vivo targets of sigma(54), the nitrogen assimilation sigma factor, and estimate that there are 70 sigma(54) promoters in total. Immunoprecipitation assays have been performed to further evaluate the efficiency of our approaches. In addition, promoter consensus binding search and primer extension assay helped us to identify a new sigma(54) promoter carried by insB-5 in the upstream of flhDC operon. The involvement of sigma(54) in flagellar biosynthesis in sequenced E. coli strain MG1655 indicates a fluid gene regulation phenomenon carried by some mobile elements in bacteria genome.

The alternative sigma factor RpoN is an essential colonization factor of Vibrio cholerae and controls important cellular functions including motility and type VI secretion (T6SS). The RpoN regulon has yet to be clearly defined in T6SS-active V. cholerae isolates, which use T6SS to target both bacterial competitors and eukaryotic cells. We hypothesize that T6SS-dependent secreted effectors are co-regulated by RpoN. To systemically identify RpoN-controlled genes, we used chromatin immunoprecipitation coupled with sequencing (ChIP-Seq) and transcriptome analysis (RNA-Seq) to determine RpoN-binding sites and RpoN-controlled gene expression. There were 68 RpoN-binding sites and 82 operons positively controlled by RpoN, among which 37 operons had ChIP-identified binding sites. A consensus RpoN-binding motif was identified with a highly conserved thymine (-14) and an AT-rich region in the middle between the hallmark RpoN-recognized motif GG(-24)/GC(-12). There were seven new RpoN-dependent promoters in the flagellar regions. We identified a small RNA, flaX, downstream of the major flagellin gene flaA. Mutation of flaX substantially reduced motility. In contrast to previous results, we report that RpoN positively regulates the expression of hcp operons and vgrG3 that encode T6SS secreted proteins but has no effect on the expression of the main T6SS cluster encoding sheath and other structural components.

In cyanobacteria, nitrogen homeostasis is maintained by an intricate regulatory network around transcription factor NtcA. Although mechanisms controlling NtcA activity appear to be well understood, its regulon remains poorly defined. To determine the NtcA regulon during the early stages of nitrogen starvation for the model cyanobacterium Synechocystis sp. PCC 6803, we performed chromatin immunoprecipitation, followed by sequencing (ChIP-seq), in parallel with transcriptome analysis (RNA-seq). Through combining these methods, we determined 51 genes activated and 28 repressed directly by NtcA. In addition to genes associated with nitrogen and carbon metabolism, a considerable number of genes without current functional annotation were among direct targets providing a rich reservoir for further studies. The NtcA regulon also included eight non-coding RNAs, of which Ncr1071, Syr6 and NsiR7 were experimentally validated, and their putative targets were computationally predicted. Surprisingly, we found substantial NtcA binding associated with delayed expression changes indicating that NtcA can reside in a poised state controlled by other factors. Indeed, a role of PipX as modulating factor in nitrogen regulation was confirmed for selected NtcA-targets. We suggest that the indicated poised state of NtcA enables a more differentiated response to nitrogen limitation and can be advantageous in native habitats of Synechocystis.

Mycobacterium tuberculosis (Mtb) Cmr (Rv1675c) is a CRP/FNR family transcription factor known to be responsive to cAMP levels and during macrophage infections. However, Cmr's DNA binding properties, cellular targets and overall role in tuberculosis (TB) complex bacteria have not been characterized. In this study, we used experimental and computational approaches to characterize Cmr's DNA binding properties and identify a putative regulon. Cmr binds a 16-bp palindromic site that includes four highly conserved nucleotides that are required for DNA binding. A total of 368 binding sites, distributed in clusters among ~200 binding regions throughout the Mycobacterium bovis BCG genome, were identified using ChIP-seq. One of the most enriched Cmr binding sites was located upstream of the cmr promoter, and we demonstrated that expression of cmr is autoregulated. cAMP affected Cmr binding at a subset of DNA loci in vivo and in vitro, including multiple sites adjacent to members of the DosR (DevR) dormancy regulon. Our findings of cooperative binding of Cmr to these DNA regions and the regulation by Cmr of the DosR-regulated virulence gene Rv2623 demonstrate the complexity of Cmr-mediated gene regulation and suggest a role for Cmr in the biology of persistent TB infection.

Pathogenicity of Pseudomonas aeruginosa, a major cause of many acute and chronic human infections, is determined by tightly regulated expression of multiple virulence factors. Quorum sensing (QS) controls expression of many of these pathogenic determinants. Previous microarray studies have shown that the AmpC β-lactamase regulator AmpR, a member of the LysR family of transcription factors, also controls non-β-lactam resistance and multiple virulence mechanisms. Using RNA-Seq and complementary assays, this study further expands the AmpR regulon to include diverse processes such as oxidative stress, heat shock and iron uptake. Importantly, AmpR affects many of these phenotypes, in part, by regulating expression of non-coding RNAs such as rgP32, asRgsA, asPrrF1 and rgRsmZ. AmpR positively regulates expression of the major QS regulators LasR, RhlR and MvfR, and genes of the Pseudomonas quinolone system. Chromatin immunoprecipitation (ChIP)-Seq and ChIP-quantitative real-time polymerase chain reaction studies show that AmpR binds to the ampC promoter both in the absence and presence of β-lactams. In addition, AmpR directly binds the lasR promoter, encoding the QS master regulator. Comparison of the AmpR-binding sequences from the transcriptome and ChIP-Seq analyses identified an AT-rich consensus-binding motif. This study further attests to the role of AmpR in regulating virulence and physiological processes in P. aeruginosa.