Microbial Primer: The bacterial flagellum - how bacteria swim

Bacteria swim using membrane-spanning, electrochemical gradient-powered motors that rotate semi-rigid helical filaments. This primer provides a brief overview of the basic synthesis, structure and operation of these nanomachines. Details and variations on the basic system can be found in suggested further reading.

Euo is a developmental regulator that represses late genes and activates midcycle genes in Chlamydia trachomatis

IMPORTANCE

In this study, we developed a correlative approach that combined DNA immunoprecipitation-seq and RNA-seq analyses to define the regulon of the Chlamydia trachomatis transcription factor Euo. We confirmed the proposed role of Euo as a transcriptional repressor of late chlamydial genes but also showed that Euo activates transcription of a subset of midcycle genes and autoregulates its own expression via negative feedback. This study validates and expands the role of Euo as an important developmental regulator in C. trachomatis. In addition, this genome-wide correlative approach can be applied to study transcription factors in other pathogenic bacteria.

The ChvG-ChvI Regulatory Network: A Conserved Global Regulatory Circuit Among the Alphaproteobacteria with Pervasive Impacts on Host Interactions and Diverse Cellular Processes

The ChvG-ChvI two-component system is conserved among multiple Alphaproteobacteria. ChvG is a canonical two-component system sensor kinase with a single large periplasmic loop. Active ChvG directs phosphotransfer to its cognate response regulator ChvI, which controls transcription of target genes. In many alphaproteobacteria, ChvG is regulated by a third component, a periplasmic protein called ExoR, that maintains ChvG in an inactive state through direct interaction. Acidic pH stimulates proteolysis of ExoR, unfettering ChvG-ChvI to control its regulatory targets. Activated ChvI among different alphaproteobacteria controls a broad range of cellular processes, including symbiosis and virulence, exopolysaccharide production, biofilm formation, motility, type VI secretion, cellular metabolism, envelope composition, and growth. Low pH is a virulence signal in Agrobacterium tumefaciens, but in other systems, conditions that cause envelope stress may also generally activate ChvG-ChvI. There is mounting evidence that these regulators influence diverse aspects of bacterial physiology, including but not limited to host interactions.

In-depth mining of single-cell transcriptome reveals the key immune-regulated loops in age-related macular degeneration

Introduction

Age-related macular degeneration (AMD), an ever-increasing ocular disease, has become one of the leading causes of irreversible blindness. Recent advances in single-cell genomics are improving our understanding of the molecular mechanisms of AMD. However, the pathophysiology of this multifactorial disease is complicated and still an ongoing challenge. To better understand disease pathogenesis and identify effective targets, we conducted an in-depth analysis of the single-cell transcriptome of AMD.

Methods

The cell expression specificity of the gene (CESG) was selected as an index to identify the novel cell markers. A computational framework was designed to explore the cell-specific TF regulatory loops, containing the interaction of gene pattern signatures, transcription factors regulons, and differentially expressed genes.

Results

Three potential novel cell markers were DNASE1L3 for endothelial cells, ABCB5 for melanocytes, and SLC39A12 for RPE cells detected. We observed a notable change in the cell abundance and crosstalk of fibroblasts cells, melanocytes, schwann cells, and T/NK cells between AMD and controls, representing a complex cellular ecosystem in disease status. Finally, we identified six cell type related and three disease-associated ternary loops and elaborated on the robust association between key immune-pathway and AMD.

Discussion

In conclusion, this study facilitates the optimization of screening for AMD-related receptor ligand pathways and proposes to further improve the interpretability of disease associations from single-cell data. It illuminated that immune-related regulation paths could be used as potential diagnostic markers for AMD, and in the future, also as therapeutic targets, providing insights into AMD diagnosis and potential interventions.

Deciphering the global roles of Cold shock proteins in Listeria monocytogenes nutrient metabolism and stress tolerance

Listeria monocytogenes (Lm) accounts for serious public health and food safety problems owing to its stress resilience and pathogenicity. Based on their regulatory involvement in global gene expression events, cold-shock domain family proteins (Csps) are crucial in expression of various stress fitness and virulence phenotypes in bacteria. Lm possesses three Csps (CspA, CspB, and CspD) whose regulatory roles in the context of the genetic diversity of this bacterium are not yet fully understood. We examined the impacts of Csps deficiency on Lm nutrient metabolism and stress tolerance using a set of csp deletion mutants generated in different genetic backgrounds. Phenotype microarrays (PM) analysis showed that the absence of Csps in ∆cspABD reduces carbon (C-) source utilization capacity and increases Lm sensitivity to osmotic, pH, various chemical, and antimicrobial stress conditions. Single and double csp deletion mutants in different Lm genetic backgrounds were used to further dissect the roles of individual Csps in these phenotypes. Selected PM-based observations were further corroborated through targeted phenotypic assays, confirming that Csps are crucial in Lm for optimal utilization of various C-sources including rhamnose and glucose as well as tolerance against NaCl, β-phenyethylamine (PEA), and food relevant detergent stress conditions. Strain and genetic lineage background-based differences, division of labour, epistasis, and functional redundancies among the Csps were uncovered with respect to their roles in various processes including C-source utilization, cold, and PEA stress resistance. Finally, targeted transcriptome analysis was performed, revealing the activation of csp gene expression under defined stress conditions and the impact of Csps on expression regulation of selected rhamnose utilization genes. Overall, our study shows that Csps play important roles in nutrient utilization and stress responses in Lm strains, contributing to traits that are central to the public health and food safety impacts of this pathogen.

An ErbB Lineage Co-Regulon Harbors Potentially Co-Druggable Targets for Multimodal Precision Therapy in Head and Neck Squamous Cell Carcinoma

The ErbB lineage of oncogenic receptor tyrosine kinases is frequently overexpressed in head and neck squamous cell carcinomas. A common co-regulon triggered by the ErbB proteins; involving shared signaling circuitries; may harbor co-druggable targets or response biomarkers for potential future multimodal precision therapy in ErbB-driven head and neck squamous cell carcinoma. We here present a cohort-based; genome-wide analysis of 488 head and neck squamous cell carcinomas curated as part of The Cancer Genome Atlas Project to characterize genes that are significantly positively co-regulated with the four ErbB proteins and those that are shared among all ErbBs denoting a common ErbB co-regulon. Significant positive gene correlations involved hundreds of genes that were co-expressed with the four ErbB family members (q < 0.05). A common; overlapping co-regulon consisted of a core set of 268 genes that were uniformly co-regulated with all four ErbB genes and highly enriched for functions in chromatin organization and histone modifications. This high-priority set of genes contained ten putative antineoplastic drug-gene interactions. The nature and directionality of these ten drug-gene associations was an inhibiting interaction for seven (PIK3CB; PIK3C2B; HDAC4; FRK; PRKCE; EPHA4; and DYRK1A) of them in which the drug decreases the biological activity or expression of the gene target. For three (CHD4; ARID1A; and PBRM1) of the associations; the directionality of the interaction was such that the gene predicted sensitivit y to the drug suggesting utility as potential response biomarkers. Drug-gene interactions that predicted the gene product to be reduced by the drug included a variety of potential targeted molecular agent classes. This unbiased genome-wide analysis identified a target-rich environment for multimodal therapeutic approaches in tumors that are putatively ErbB-driven. The results of this study require preclinical validation before ultimately devising lines of combinatorial treatment strategies for ErbB-dependent head and neck squamous cell carcinomas that incorporate these findings.

Human Milk Oligosaccharide Utilization in Intestinal Bifidobacteria Is Governed by Global Transcriptional Regulator NagR

Bifidobacterium longum subsp. infantis is a prevalent beneficial bacterium that colonizes the human neonatal gut and is uniquely adapted to efficiently use human milk oligosaccharides (HMOs) as a carbon and energy source. Multiple studies have focused on characterizing the elements of HMO utilization machinery in B. longum subsp. infantis; however, the regulatory mechanisms governing the expression of these catabolic pathways remain poorly understood. A bioinformatic regulon reconstruction approach used in this study implicated NagR, a transcription factor from the ROK family, as a negative global regulator of gene clusters encoding lacto-N-biose/galacto-N-biose (LNB/GNB), lacto-N-tetraose (LNT), and lacto-N-neotetraose (LNnT) utilization pathways in B. longum subsp. infantis. This conjecture was corroborated by transcriptome profiling upon nagR genetic inactivation and experimental assessment of binding of recombinant NagR to predicted DNA operators. The latter approach also implicated N-acetylglucosamine (GlcNAc), a universal intermediate of LNT and LNnT catabolism, and its phosphorylated derivatives as plausible NagR transcriptional effectors. Reconstruction of NagR regulons in various Bifidobacterium lineages revealed multiple potential regulon expansion events, suggesting evolution from a local regulator of GlcNAc catabolism in ancestral bifidobacteria to a global regulator controlling the utilization of mixtures of GlcNAc-containing host glycans in B. longum subsp. infantis and Bifidobacterium bifidum. IMPORTANCE The predominance of bifidobacteria in the gut of breastfed infants is attributed to the ability of these bacteria to metabolize human milk oligosaccharides (HMOs). Thus, individual HMOs such as lacto-N-tetraose (LNT) and lacto-N-neotetraose (LNnT) are considered promising prebiotics that would stimulate the growth of bifidobacteria and confer multiple health benefits to preterm and malnourished children suffering from impaired (stunted) gut microbiota development. However, the rational selection of HMO-based prebiotics is hampered by the incomplete knowledge of regulatory mechanisms governing HMO utilization in target bifidobacteria. This study describes NagR-mediated transcriptional regulation of LNT and LNnT utilization in Bifidobacterium longum subsp. infantis. The elucidated regulatory network appears optimally adapted to simultaneous utilization of multiple HMOs, providing a rationale to add HMO mixtures (rather than individual components) to infant formulas. The study also provides insights into the evolutionary trajectories of complex regulatory networks controlling carbohydrate metabolism in bifidobacteria.

Classification of bladder cancer cell lines according to regulon activity

Bladder cancer is one of the most frequent cancers and causes more than 150.000 deaths each year. During the last decade, several studies provided important aspects about genomic characterization, consensus subgroup definition, and transcriptional regulation of bladder cancer. Still, much more research needs to be done to characterize molecular signatures of this cancer in depth. At this point, the use of bladder cancer cell lines is quite useful for the identification and test of new signatures. In this study, we classified the bladder cancer cell lines according to the activities of regulons implicated in the regulation of primary bladder tumors. Our regulon gene expression-based classification revealed three groups, neuronal-basal (NB), luminal-papillary (LP), and basal-squamous (BS). These regulon gene expression-based classifications showed a quite good concordance with the consensus subgroups assigned by the primary bladder cancer classifier. Importantly, we identified FGFR1 regulon to be involved in the characterization of the NB group, where neuroendocrine signature genes were significantly upregulated, and further β-catenin was shown to have significantly higher nuclear localization. LP groups were mainly driven by the regulons ERBB2, FOXA1, GATA3, and PPARG, and they showed upregulation of the genes involved in epithelial differentiation and urogenital development, while the activity of EGFR, FOXM1, STAT3, and HIF1A was implicated for the regulation of BS group. Collectively, our results and classifications may serve as an important guide for the selection and use of bladder cancer cell lines for experimental strategies, which aim to manipulate regulons critical for bladder cancer development.

Overlapping SigH and SigE sigma factor regulons in Corynebacterium glutamicum

The sigma H (σ^Η) and sigma E (σ^E) subunits of Corynebacterium glutamicum RNA polymerase belong to Group 4 of sigma factors, also called extracytoplasmic function (ECF) sigma factors. Genes of the C. glutamicum σ^Η regulon that are involved in heat and oxidative stress response have already been defined, whereas the genes of the σ^E regulon, which is involved in cell surface stress response, have not been explored until now. Using the C. glutamicum RES167 strain and its derivative C. glutamicum ΔcseE with a deletion in the anti-σ^Ε gene, differential gene expression was analyzed by RNA sequencing. We found 296 upregulated and 398 downregulated genes in C. glutamicum ΔcseE compared to C. glutamicum RES167. To confirm the functional link between σ^Ε and the corresponding promoters, we tested selected promoters using the in vivo two-plasmid system with gfpuv as a reporter gene and by in vitro transcription. Analyses with RNAP+σ^Η and RNAP+σ^Ε, which were previously shown to recognize similar promoters, proved that the σ^Η and σ^E regulons significantly overlap. The σ^E-controlled genes were found to be involved for example in protein quality control (dnaK, dnaJ2, clpB, and clpC), the regulation of Clp proteases (clgR), and membrane integrity maintenance. The single-promoter analyses with σ^Η and σ^Ε revealed that there are two groups of promoters: those which are exclusively σ^Η-specific, and the other group of promoters, which are σ^Η/σ^E-dependent. No exclusively σ^E-dependent promoter was detected. We defined the consensus sequences of exclusively σ^Η-regulated promotors to be -35 GGAAt and - 10 GTT and σ^Η/σ^E-regulated promoters to be -35 GGAAC and - 10 cGTT. Fifteen genes were found to belong to the σ^Η/σ^Ε regulon. Homology modeling showed that there is a specific interaction between Met170 in σ^Η and the nucleotides -31 and - 30 within the non-coding strand (AT or CT) of the σ^Η-dependent promoters. In σ^E, Arg185 was found to interact with the nucleotides GA at the same positions in the σ^E-dependent promoters.

Extensive Translational Regulation through the Proliferative Transition of Trypanosoma cruzi Revealed by Multi-Omics

Trypanosoma cruzi is the etiological agent for Chagas disease, a neglected parasitic disease in Latin America. Gene transcription control governs the eukaryotic cell replication but is absent in trypanosomatids; thus, it must be replaced by posttranscriptional regulatory events. We investigated the entrance into the T. cruzi replicative cycle using ribosome profiling and proteomics on G₁/S epimastigote cultures synchronized with hydroxyurea. We identified 1,784 translationally regulated genes (change > 2, false-discovery rate [FDR] < 0.05) and 653 differentially expressed proteins (change > 1.5, FDR < 0.05), respectively. A major translational remodeling accompanied by an extensive proteome change is found, while the transcriptome remains largely unperturbed at the replicative entrance of the cell cycle. The differentially expressed genes comprise specific cell cycle processes, confirming previous findings while revealing candidate cell cycle regulators that undergo previously unnoticed translational regulation. Clusters of genes showing a coordinated regulation at translation and protein abundance share related biological functions such as cytoskeleton organization and mitochondrial metabolism; thus, they may represent posttranscriptional regulons. The translatome and proteome of the coregulated clusters change in both coupled and uncoupled directions, suggesting that complex cross talk between the two processes is required to achieve adequate protein levels of different regulons. This is the first simultaneous assessment of the transcriptome, translatome, and proteome of trypanosomatids, which represent a paradigm for the absence of transcriptional control. The findings suggest that gene expression chronology along the T. cruzi cell cycle is controlled mainly by translatome and proteome changes coordinated using different mechanisms for specific gene groups. IMPORTANCE Trypanosoma cruzi is an ancient eukaryotic unicellular parasite causing Chagas disease, a potentially life-threatening illness that affects 6 to 7 million people, mostly in Latin America. The antiparasitic treatments for the disease have incomplete efficacy and adverse reactions; thus, improved drugs are needed. We study the mechanisms governing the replication of the parasite, aiming to find differences with the human host, valuable for the development of parasite-specific antiproliferative drugs. Transcriptional regulation is essential for replication in most eukaryotes, but in trypanosomatids, it must be replaced by subsequent gene regulation steps since they lack transcription initiation control. We identified the genome-wide remodeling of mRNA translation and protein abundance during the entrance to the replicative phase of the cell cycle. We found that translation is strongly regulated, causing variation in protein levels of specific cell cycle processes, representing the first simultaneous study of the translatome and proteome in trypanosomatids.

Genome-Wide Transcription Factor DNA Binding Sites and Gene Regulatory Networks in Clostridium thermocellum

Clostridium thermocellum is a thermophilic bacterium recognized for its natural ability to effectively deconstruct cellulosic biomass. While there is a large body of studies on the genetic engineering of this bacterium and its physiology to-date, there is limited knowledge in the transcriptional regulation in this organism and thermophilic bacteria in general. The study herein is the first report of a large-scale application of DNA-affinity purification sequencing (DAP-seq) to transcription factors (TFs) from a bacterium. We applied DAP-seq to > 90 TFs in C. thermocellum and detected genome-wide binding sites for 11 of them. We then compiled and aligned DNA binding sequences from these TFs to deduce the primary DNA-binding sequence motifs for each TF. These binding motifs are further validated with electrophoretic mobility shift assay (EMSA) and are used to identify individual TFs’ regulatory targets in C. thermocellum. Our results led to the discovery of novel, uncharacterized TFs as well as homologues of previously studied TFs including RexA-, LexA-, and LacI-type TFs. We then used these data to reconstruct gene regulatory networks for the 11 TFs individually, which resulted in a global network encompassing the TFs with some interconnections. As gene regulation governs and constrains how bacteria behave, our findings shed light on the roles of TFs delineated by their regulons, and potentially provides a means to enable rational, advanced genetic engineering of C. thermocellum and other organisms alike toward a desired phenotype.

Transcriptional Regulation of Plant Biomass Degradation and Carbohydrate Utilization Genes in the Extreme Thermophile Caldicellulosiruptor bescii

Extremely thermophilic bacteria from the genus Caldicellulosiruptor can degrade polysaccharide components of plant cell walls and subsequently utilize the constituting mono- and oligosaccharides. Through metabolic engineering, ethanol and other industrially important end products can be produced. Previous experimental studies identified a variety of carbohydrate-active enzymes in model species Caldicellulosiruptor saccharolyticus and Caldicellulosiruptor bescii, while prior transcriptomic experiments identified their putative carbohydrate uptake transporters. We investigated the mechanisms of transcriptional regulation of carbohydrate utilization genes using a comparative genomics approach applied to 14 Caldicellulosiruptor species. The reconstruction of carbohydrate utilization regulatory network includes the predicted binding sites for 34 mostly local regulators and point to the regulatory mechanisms controlling expression of genes involved in degradation of plant biomass. The Rex and CggR regulons control the central glycolytic and primary redox reactions. The identified transcription factor binding sites and regulons were validated with transcriptomic and transcription start site experimental data for C. bescii grown on cellulose, cellobiose, glucose, xylan, and xylose. The XylR and XynR regulons control xylan-induced transcriptional response of genes involved in degradation of xylan and xylose utilization. The reconstructed regulons informed the carbohydrate utilization reconstruction analysis and improved functional annotations of 51 transporters and 11 catabolic enzymes. Using gene deletion, we confirmed that the shared ATPase component MsmK is essential for growth on oligo- and polysaccharides but not for the utilization of monosaccharides. By elucidating the carbohydrate utilization framework in C. bescii, strategies for metabolic engineering can be pursued to optimize yields of bio-based fuels and chemicals from lignocellulose. IMPORTANCE To develop functional metabolic engineering platforms for nonmodel microorganisms, a comprehensive understanding of the physiological and metabolic characteristics is critical. Caldicellulosiruptor bescii and other species in this genus have untapped potential for conversion of unpretreated plant biomass into industrial fuels and chemicals. The highly interactive and complex machinery used by C. bescii to acquire and process complex carbohydrates contained in lignocellulose was elucidated here to complement related efforts to develop a metabolic engineering platform with this bacterium. Guided by the findings here, a clearer picture of how C. bescii natively drives carbohydrate utilization is provided and strategies to engineer this bacterium for optimal conversion of lignocellulose to commercial products emerge.

Osmoregulation of the transcriptome of the hypothalamic supraoptic nucleus: A resource for the community

The hypothalamic supraoptic nucleus (SON) is a core osmoregulatory control centre that deciphers information about the metabolic state of the organism and orchestrates appropriate homeostatic (endocrine) and allostatic (behavioural) responses. We have used RNA sequencing to describe the polyadenylated transcriptome of the SON of the male Wistar Han rat. These data have been mined to generate comprehensive catalogues of functional classes of genes (enzymes, transcription factors, endogenous peptides, G protein coupled receptors, transporters, catalytic receptors, channels and other pharmacological targets) expressed in this nucleus in the euhydrated state, and that together form the basal substrate for its physiological interactions. We have gone on to show that fluid deprivation for 3 days (dehydration) results in changes in the expression levels of 2247 RNA transcripts, which have similarly been functionally catalogued, and further mined to describe enriched gene categories and putative regulatory networks (Regulons) that may have physiological importance in SON function related plasticity. We hope that the revelation of these genes, pathways and networks, most of which have no characterised roles in the SON, will encourage the neuroendocrine community to pursue new investigations into the new 'known-unknowns' reported in the present study.

The AraC-Type Transcriptional Regulator GliR (PA3027) Activates Genes of Glycerolipid Metabolism in Pseudomonas aeruginosa

Pseudomonas aeruginosa encodes a large set of transcriptional regulators (TRs) that modulate and manage cellular metabolism to survive in variable environmental conditions including that of the human body. The AraC family regulators are an abundant group of TRs in bacteria, mostly acting as gene expression activators, controlling diverse cellular functions (e.g., carbon metabolism, stress response, and virulence). The PA3027 protein from P. aeruginosa has been classified in silico as a putative AraC-type TR. Transcriptional profiling of P. aeruginosa PAO1161 overexpressing PA3027 revealed a spectacular increase in the mRNA levels of PA3026-PA3024 (divergent to PA3027), PA3464, and PA3342 genes encoding proteins potentially involved in glycerolipid metabolism. Concomitantly, chromatin immunoprecipitation-sequencing (ChIP-seq) analysis revealed that at least 22 regions are bound by PA3027 in the PAO1161 genome. These encompass promoter regions of PA3026, PA3464, and PA3342, showing the major increase in expression in response to PA3027 excess. In Vitro DNA binding assay confirmed interactions of PA3027 with these regions. Furthermore, promoter-reporter assays in a heterologous host showed the PA3027-dependent activation of the promoter of the PA3026-PA3024 operon. Two motifs representing the preferred binding sites for PA3027, one localized upstream and one overlapping with the -35 promoter sequence, were identified in PA3026p and our data indicate that both motifs are required for full activation of this promoter by PA3027. Overall, the presented data show that PA3027 acts as a transcriptional regulator in P. aeruginosa, activating genes likely engaged in glycerolipid metabolism. The GliR name, from a glycerolipid metabolism regulator, is proposed for PA3027 of P. aeruginosa.

Gene Regulation Network Analysis on Human Prostate Orthografts Highlights a Potential Role for the JMJD6 Regulon in Clinical Prostate Cancer

Simple Summary

Prostate cancer is a very common malignancy worldwide. Treatment resistant prostate cancer poses a big challenge to clinicians and is the second most common cause of premature death in men with cancer. Gene expression analysis has been performed on clinical tumours but to date none of the gene expression-based biomarkers for prostate cancer have been successfully integrated to into clinical practice to improve patient management and treatment choice. We applied a novel laboratory prostate cancer model to mimic clinical hormone responsive and resistant prostate cancer and tested whether a network of genes similarly regulated by transcription factors (gene products that control the expression of target genes) are associated with patient outcome. We identified regulons (networks of genes similarly regulated) from our preclinical prostate cancer models and further evaluated the top ranked JMJD6 gene related regulated network in three independent clinical patient cohorts.

Abstract

Background: Prostate cancer (PCa) is the second most common tumour diagnosed in men. Tumoral heterogeneity in PCa creates a significant challenge to develop robust prognostic markers and novel targets for therapy. An analysis of gene regulatory networks (GRNs) in PCa may provide insight into progressive PCa. Herein, we exploited a graph-based enrichment score to integrate data from GRNs identified in preclinical prostate orthografts and differentially expressed genes in clinical resected PCa. We identified active regulons (transcriptional regulators and their targeted genes) associated with PCa recurrence following radical prostatectomy. Methods: The expression of known transcription factors and co-factors was analysed in a panel of prostate orthografts (n = 18). We searched for genes (as part of individual GRNs) predicted to be regulated by the highest number of transcriptional factors. Using differentially expressed gene analysis (on a per sample basis) coupled with gene graph enrichment analysis, we identified candidate genes and associated GRNs in PCa within the UTA cohort, with the most enriched regulon being JMJD6, which was further validated in two additional cohorts, namely EMC and ICGC cohorts. Cox regression analysis was performed to evaluate the association of the JMJD6 regulon activity with disease-free survival time in the three clinical cohorts as well as compared to three published prognostic gene signatures (TMCC11, BROMO-10 and HYPOXIA-28). Results: 1308 regulons were correlated to transcriptomic data from the three clinical prostatectomy cohorts. The JMJD6 regulon was identified as the top enriched regulon in the UTA cohort and again validated in the EMC cohort as the top-ranking regulon. In both UTA and EMC cohorts, the JMJD6 regulon was significantly associated with cancer recurrence. Active JMJD6 regulon also correlated with disease recurrence in the ICGC cohort. Furthermore, Kaplan–Meier analysis confirmed shorter time to recurrence in patients with active JMJD6 regulon for all three clinical cohorts (UTA, EMC and ICGC), which was not the case for three published prognostic gene signatures (TMCC11, BROMO-10 and HYPOXIA-28). In multivariate analysis, the JMJD6 regulon status significantly predicted disease recurrence in the UTA and EMC, but not ICGC datasets, while none of the three published signatures significantly prognosticate for cancer recurrence. Conclusions: We have characterised gene regulatory networks from preclinical prostate orthografts and applied transcriptomic data from three clinical cohorts to evaluate the prognostic potential of the JMJD6 regulon.

So alike yet so different. Differential expression of the long non-coding RNAs NORAD and HCG11 in breast cancer subtypes

Breast cancer (BC) is a heterogeneous disease, and it is the leading cause of death among women. NORAD and HCG11 are highly similar lncRNAs that present binding sites for PUMILIO proteins. PUMILIO acts on hundreds of mRNA targets, contributing to the modulation of gene expression. We analyzed the expression levels of NORAD and HCG11 in the BC subtypes luminal A (LA) and basal-like (BL), and the regulatory networks associated with these lncRNAs. In the analysis of TCGA cohort (n=329) and Brazilian BC samples (n=44), NORAD was up-regulated in LA while HCG11 was up-regulated in BL subtype. An increased expression of NORAD is associated with reduced disease-free survival in basal-like patients (p = 0.002), which suggests that its prognostic value could be different in specific subtypes. The biological pathways observed for the HCG11 network are linked to the epithelial-to-mesenchymal transition; while NORAD associated pathways appear to be related to luminal epithelial cell transformation. NORAD and HCG11 regulons respectively present 36% and 21.5% of PUMILIO targets, which suggests that these lncRNAs act as a decoy for PUMILIO. These lncRNAs seem to work as players in the differentiation process that drives breast cells to acquire distinct phenotypes related to a specific BC subtype.

Genome-Wide Mapping Reveals Complex Regulatory Activities of BfmR in Pseudomonas aeruginosa

BfmR is a response regulator that modulates diverse pathogenic phenotypes and induces an acute-to-chronic virulence switch in Pseudomonas aeruginosa, an important human pathogen causing serious nosocomial infections. However, the mechanisms of action of BfmR remain largely unknown. Here, using chromatin immunoprecipitation followed by high-throughput sequencing (ChIP-seq), we showed that 174 chromosomal regions of P. aeruginosa MPAO1 genome were highly enriched by coimmunoprecipitation with a C-terminal Flag-tagged BfmR. Integration of these data with global transcriptome analyses revealed that 172 genes in 106 predicted transcription units are potential targets for BfmR. We determined that BfmR binds to and modulates the promoter activity of genes encoding transcriptional regulators CzcR, ExsA, and PhoB. Intriguingly, BfmR bound to the promoters of a number of genes belong to either CzcR or PhoB regulon, or both, indicating that CzcRS and PhoBR two-component systems (TCSs) deeply feed into the BfmR-mediated regulatory network. In addition, we demonstrated that phoB is required for BfmR to promote the biofilm formation by P. aeruginosa. These results delineate the direct BfmR regulon and exemplify the complexity of BfmR-mediated regulation of cellular functions in P. aeruginosa.

Regulatory Effects of CsrA in Vibrio cholerae

Vibrio cholerae, a Gram-negative bacterium, is a natural inhabitant of the aqueous environment. However, once ingested, this bacterium can colonize the human host and cause the disease cholera.

CsrA is a posttranscriptional global regulator in Vibrio cholerae. Although CsrA is critical for V. cholerae survival within the mammalian host, the regulatory targets of CsrA remain mostly unknown. To identify pathways controlled by CsrA, RNA-seq transcriptome analysis was carried out by comparing the wild type and the csrA mutant grown to early exponential, mid-exponential, and stationary phases of growth. This enabled us to identify the global effects of CsrA-mediated regulation throughout the V. cholerae growth cycle. We found that CsrA regulates 22% of the V. cholerae transcriptome, with significant regulation within the gene ontology (GO) processes that involve amino acid transport and metabolism, central carbon metabolism, lipid metabolism, iron uptake, and flagellum-dependent motility. Through CsrA-RNA coimmunoprecipitation experiments, we found that CsrA binds to multiple mRNAs that encode regulatory proteins. These include transcripts encoding the major sigma factors RpoS and RpoE, which may explain how CsrA regulation affects such a large proportion of the V. cholerae transcriptome. Other direct targets include flrC, encoding a central regulator in flagellar gene expression, and aphA, encoding the virulence gene transcription factor AphA. We found that CsrA binds to the aphA mRNA both in vivo and in vitro, and CsrA significantly increases AphA protein synthesis. The increase in AphA was due to increased translation, not transcription, in the presence of CsrA, consistent with CsrA binding to the aphA transcript and enhancing its translation. CsrA is required for the virulence of V. cholerae and this study illustrates the central role of CsrA in virulence gene regulation.

monocytogenes PrfA regulon identifies six novel putative PrfA-regulated genes

The transcriptional activator Positive Regulatory Factor A (PrfA) regulates expression of genes essential for virulence in Listeria monocytogenes. To define the PrfA regulon, the 10403S wildtype (WT) strain, a constitutively active prfA* mutant, and an isogenic ∆prfA mutant were grown under PrfA-inducing conditions in a medium containing glucose-1-phosphate and pre-treated with 0.2% activated charcoal. RNA-seq-generated transcript levels were compared as follows: (i) prfA* and WT; (ii) WT and ∆prfA and (iii) prfA* and ∆prfA. Significantly higher transcript levels in the induced WT or constitutively active PrfA* were identified for 18 genes and 2 ncRNAs in at least one of the three comparisons. These genes included: (i) 10/12 of the genes previously identified as directly PrfA-regulated; (ii) 2 genes previously identified as PrfA-regulated, albeit likely indirectly; and (iii) 6 genes newly identified as PrfA-regulated, including one (LMRG_0 2046) with a σA-dependent promoter and PrfA box located within an upstream open reading frame. LMRG_0 2046, which encodes a putative cyanate permease, is reported to be downregulated by a σB-dependent anti-sense RNA. This newly identified overlap between the σB and PrfA regulons highlights the complexity of regulatory networks important for fine-tuning bacterial gene expression in response to the rapidly changing environmental conditions associated with infection.

Identification of Qk as a Glial Precursor Cell Marker that Governs the Fate Specification of Neural Stem Cells to a Glial Cell Lineage

During brain development, neural stem cells (NSCs) initially produce neurons and change their fate to generate glias. While the regulation of neurogenesis is well characterized, specific markers for glial precursor cells (GPCs) and the master regulators for gliogenesis remain unidentified. Accumulating evidence suggests that RNA-binding proteins (RBPs) have significant roles in neuronal development and function, as they comprehensively regulate the expression of target genes in a cell-type-specific manner. We systematically investigated the expression profiles of 1,436 murine RBPs in the developing mouse brain and identified quaking (Qk) as a marker of the putative GPC population. Functional analysis of the NSC-specific Qk-null mutant mouse revealed the key role of Qk in astrocyte and oligodendrocyte generation and differentiation from NSCs. Mechanistically, Qk upregulates gliogenic genes via quaking response elements in their 3′ untranslated regions. These results provide crucial directions for identifying GPCs and deciphering the regulatory mechanisms of gliogenesis from NSCs.

•

Differential expression analysis identified Qk as a glial precursor cell marker

•

Loss of Qk ablated both astrocyte and OL production from neural stem cells

•

Qk^−/− NSCs failed to become glia and aberrantly expressed neural genes

•

Qk comprehensively upregulates essential genes for gliogenesis as regulons via QREs

Redefining the Clostridioides difficile σB Regulon: σB Activates Genes Involved in Detoxifying Radicals That Can Result from the Exposure to Antimicrobials and Hydrogen Peroxide

In many Gram-positive bacteria, the general stress response is regulated at the transcriptional level by the alternative sigma factor sigma B (σ^B). In C. difficile, σ^B has been implicated in protection against stressors such as reactive oxygen species (ROS) and antimicrobial compounds. Here, we used an anti-σ^B antibody to demonstrate time-limited overproduction of σ^B in C. difficile despite its toxicity at higher cellular concentrations. This toxicity eventually led to the loss of the plasmid used for anhydrotetracycline-induced σ^B gene expression. Inducible σ^B overproduction uncouples σ^B expression from its native regulatory network and allows for the refinement of the previously proposed σ^B regulon. At least 32% of the regulon was found to consist of genes involved in the response to reactive radicals. Direct gene activation by C. difficile σ^B was demonstrated through in vitro runoff transcription of specific target genes (cd0350, cd3614, cd3605, and cd2963). Finally, we demonstrated that different antimicrobials and hydrogen peroxide induce these genes in a manner dependent on this sigma factor, using a plate-based luciferase reporter assay. Together, our work suggests that lethal exposure to antimicrobials may result in the formation of toxic radicals that lead to σ^B-dependent gene activation.IMPORTANCE Sigma B is the alternative sigma factor governing stress response in many Gram-positive bacteria. In C. difficile, a sigB mutant shows pleiotropic transcriptional effects. Here, we determine genes that are likely direct targets of σ^B by evaluating the transcriptional effects of σ^B overproduction, provide biochemical evidence of direct transcriptional activation by σ^B, and show that σ^B-dependent genes can be activated by antimicrobials. Together, our data suggest that σ^B is a key player in dealing with toxic radicals.

Sigma 54-Regulated Transcription Is Associated with Membrane Reorganization and Type III Secretion Effectors during Conversion to Infectious Forms of Chlamydia trachomatis

The factors that control the growth and infectious processes for Chlamydia are still poorly understood. This study used recently developed genetic tools to determine the regulon for one of the key transcription factors encoded by Chlamydia, sigma 54. Surrogate and computational analyses provide additional support for the hypothesis that sigma 54 plays a key role in controlling the expression of many components critical to converting and enabling the infectious capability of Chlamydia. These components include those that remodel the membrane for the extracellular environment and incorporation of an arsenal of type III secretion effectors in preparation for infecting new cells.

Chlamydia bacteria are obligate intracellular organisms with a phylum-defining biphasic developmental cycle that is intrinsically linked to its ability to cause disease. The progression of the chlamydial developmental cycle is regulated by the temporal expression of genes predominantly controlled by RNA polymerase sigma (σ) factors. Sigma 54 (σ⁵⁴) is one of three sigma factors encoded by Chlamydia for which the role and regulon are unknown. CtcC is part of a two-component signal transduction system that is requisite for σ⁵⁴ transcriptional activation. CtcC activation of σ⁵⁴ requires phosphorylation, which relieves inhibition by the CtcC regulatory domain and enables ATP hydrolysis by the ATPase domain. Prior studies with CtcC homologs in other organisms have shown that expression of the ATPase domain alone can activate σ⁵⁴ transcription. Biochemical analysis of CtcC ATPase domain supported the idea of ATP hydrolysis occurring in the absence of the regulatory domain, as well as the presence of an active-site residue essential for ATPase activity (E242). Using recently developed genetic approaches in Chlamydia to induce expression of the CtcC ATPase domain, a transcriptional profile was determined that is expected to reflect the σ⁵⁴ regulon. Computational evaluation revealed that the majority of the differentially expressed genes were preceded by highly conserved σ⁵⁴ promoter elements. Reporter gene analyses using these putative σ⁵⁴ promoters reinforced the accuracy of the model of the proposed regulon. Investigation of the gene products included in this regulon supports the idea that σ⁵⁴ controls expression of genes that are critical for conversion of Chlamydia from replicative reticulate bodies into infectious elementary bodies.

Beyond Trees: Regulons and Regulatory Motif Characterization

Trees and their seeds regulate their germination, growth, and reproduction in response to environmental stimuli. These stimuli, through signal transduction, trigger transcription factors that alter the expression of various genes leading to the unfolding of the genetic program. A regulon is conceptually defined as a set of target genes regulated by a transcription factor by physically binding to regulatory motifs to accomplish a specific biological function, such as the CO-FT regulon for flowering timing and fall growth cessation in trees. Only with a clear characterization of regulatory motifs, can candidate target genes be experimentally validated, but motif characterization represents the weakest feature of regulon research, especially in tree genetics. I review here relevant experimental and bioinformatics approaches in characterizing transcription factors and their binding sites, outline problems in tree regulon research, and demonstrate how transcription factor databases can be effectively used to aid the characterization of tree regulons.

Systematic Reconstruction of the Complete Two-Component Sensorial Network in Staphylococcus aureus

Bacteria are able to sense environmental conditions and respond accordingly. Their sensorial system relies on pairs of sensory and regulatory proteins, known as two-component systems (TCSs). The majority of bacteria contain dozens of TCSs, each of them responsible for sensing and responding to a different range of signals. Traditionally, the function of each TCS has been determined by analyzing the changes in gene expression caused by the absence of individual TCSs. Here, we used a bacterial strain deprived of the complete TC sensorial system to introduce, one by one, the active form of every TCS. This gain-of-function strategy allowed us to identify the changes in gene expression conferred by each TCS without interference of other members of the family.

In bacteria, adaptation to changes in the environment is mainly controlled through two-component signal transduction systems (TCSs). Most bacteria contain dozens of TCSs, each of them responsible for sensing a different range of signals and controlling the expression of a repertoire of target genes (regulon). Over the years, identification of the regulon controlled by each individual TCS in different bacteria has been a recurrent question. However, limitations associated with the classical approaches used have left our knowledge far from complete. In this report, using a pioneering approach in which a strain devoid of the complete nonessential TCS network was systematically complemented with the constitutively active form of each response regulator, we have reconstituted the regulon of each TCS of S. aureus in the absence of interference between members of the family. Transcriptome sequencing (RNA-Seq) and proteomics allowed us to determine the size, complexity, and insulation of each regulon and to identify the genes regulated exclusively by one or many TCSs. This gain-of-function strategy provides the first description of the complete TCS regulon in a living cell, which we expect will be useful to understand the pathobiology of this important pathogen.

IMPORTANCE Bacteria are able to sense environmental conditions and respond accordingly. Their sensorial system relies on pairs of sensory and regulatory proteins, known as two-component systems (TCSs). The majority of bacteria contain dozens of TCSs, each of them responsible for sensing and responding to a different range of signals. Traditionally, the function of each TCS has been determined by analyzing the changes in gene expression caused by the absence of individual TCSs. Here, we used a bacterial strain deprived of the complete TC sensorial system to introduce, one by one, the active form of every TCS. This gain-of-function strategy allowed us to identify the changes in gene expression conferred by each TCS without interference of other members of the family.

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

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

paratuberculosis (MAP)

Intracellular iron concentration is tightly regulated to maintain cell viability. Iron plays important roles in electron transport, nucleic acid synthesis, and oxidative stress. A Mycobacterium avium subsp. paratuberculosis (MAP)-specific genomic island carries a putative metal transport operon that includes MAP3773c, which encodes a Fur-like protein. Although well characterized as a global regulator of iron homeostasis in multiple bacteria, the function of Fur (ferric uptake regulator) in MAP is unknown as this organism also carries IdeR (iron dependent regulator), a native iron regulatory protein specific to mycobacteria. Computational analysis using PRODORIC identified 23 different pathways involved in respiration, metabolism, and virulence that were likely regulated by MAP3773c. Thus, chromatin immunoprecipitation followed by high-throughput sequencing (ChIP-seq) was performed to confirm the putative regulon of MAP3773c (Fur-like protein) in MAP. ChIP-Seq revealed enriched binding to 58 regions by Fur under iron-replete and -deplete conditions, located mostly within open reading frames (ORFs). Three ChIP peaks were identified in genes that are directly related to iron regulation: MAP3638c (hemophore-like protein), MAP3736c (Fur box), and MAP3776c (ABC transporter). Fur box consensus sequence was identified, and binding specificity and dependence on Mn2+ availability was confirmed by a chemiluminescent electrophoresis mobility shift assay (EMSA). The results confirmed that MAP3773c is a Fur ortholog that recognizes a 19 bp DNA sequence motif (Fur box) and it is involved in metal homeostasis. This work provides a regulatory network of MAP Fur binding sites during iron-replete and -deplete conditions, highlighting unique properties of Fur regulon in MAP.

Analysis of HrpG regulons and HrpG-interacting proteins by ChIP-seq and affinity proteomics in Xanthomonas campestris

Gamma-proteobacteria Xanthomonas spp. cause at least 350 different plant diseases among important agricultural crops, which result in serious yield losses. Xanthomonas spp. rely mainly on the type III secretion system (T3SS) to infect their hosts and induce a hypersensitive response in nonhosts. HrpG, the master regulator of the T3SS, plays the dominant role in bacterial virulence. In this study, we used chromatin immunoprecipitation followed by sequencing (ChIP-seq) and tandem affinity purification (TAP) to systematically characterize the HrpG regulon and HrpG interacting proteins in vivo. We obtained 186 candidate HrpG downstream genes from the ChIP-seq analysis, which represented the genomic-wide regulon spectrum. A consensus HrpG-binding motif was obtained and three T3SS genes, hpa2, hrcU, and hrpE, were confirmed to be directly transcriptionally activated by HrpG in the inducing medium. A total of 273 putative HrpG interacting proteins were identified from the TAP data and the DNA-binding histone-like HU protein of Xanthomonas campestris pv. campestris (HU_xcc ) was proved to be involved in bacterial virulence by increasing the complexity and intelligence of the bacterial signalling pathways in the T3SS.

Hierarchical and Dynamic Regulation of Defense-Responsive Specialized Metabolism by WRKY and MYB Transcription Factors

The plant kingdom produces hundreds of thousands of specialized bioactive metabolites, some with pharmaceutical and biotechnological importance. Their biosynthesis and function have been studied for decades, but comparatively less is known about how transcription factors with overlapping functions and contrasting regulatory activities coordinately control the dynamics and output of plant specialized metabolism. Here, we performed temporal studies on pathogen-infected intact host plants with perturbed transcription factors. We identified WRKY33 as the condition-dependent master regulator and MYB51 as the dual functional regulator in a hierarchical gene network likely responsible for the gene expression dynamics and metabolic fluxes in the camalexin and 4-hydroxy-indole-3-carbonylnitrile (4OH-ICN) pathways. This network may have also facilitated the regulatory capture of the newly evolved 4OH-ICN pathway in Arabidopsis thaliana by the more-conserved transcription factor MYB51. It has long been held that the plasticity of plant specialized metabolism and the canalization of development should be differently regulated; our findings imply a common hierarchical regulatory architecture orchestrated by transcription factors for specialized metabolism and development, making it an attractive target for metabolic engineering.

The Staphylococcus aureus ArlRS two-component system regulates virulence factor expression through MgrA

The Gram-positive bacterium, Staphylococcus aureus, is a versatile pathogen that can sense and adapt to a wide variety of environments within the human host, in part through its 16 two-component regulatory systems. The ArlRS two-component system has been shown to affect many cellular processes in S. aureus, including autolysis, biofilm formation, capsule synthesis and virulence. Yet the molecular details of this regulation remained largely unknown. We used RNA sequencing to identify the ArlRS regulon, and found 70% overlap with that of the global regulator MgrA. These genes included cell wall-anchored adhesins (ebh, sdrD), polysaccharide and capsule synthesis genes, cell wall remodeling genes (lytN, ddh), the urease operon, genes involved in metal transport (feoA, mntH, sirA), anaerobic metabolism genes (adhE, pflA, nrdDG) and a large number of virulence factors (lukSF, lukAB, nuc, gehB, norB, chs, scn and esxA). We show that ArlR directly activates expression of mgrA and identify a probable ArlR-binding site (TTTTCTCAT-N4 -TTTTAATAA). A highly similar sequence is also found in the spx P2 promoter, which was recently shown to be regulated by ArlRS. We also demonstrate that ArlS has kinase activity toward ArlR in vitro, although it has slower kinetics than other similar histidine kinases.

PAX8 regulon in human ovarian cancer links lineage dependency with epigenetic vulnerability to HDAC inhibitors

PAX8 is a prototype lineage-survival oncogene in epithelial ovarian cancer. However, neither its underlying pro-tumorigenic mechanisms nor potential therapeutic implications have been adequately elucidated. Here, we identified an ovarian lineage-specific PAX8 regulon using modified cancer outlier profile analysis, in which PAX8-FGF18 axis was responsible for promoting cell migration in an autocrine fashion. An image-based drug screen pinpointed that PAX8 expression was potently inhibited by small-molecules against histone deacetylases (HDACs). Mechanistically, HDAC blockade altered histone H3K27 acetylation occupancies and perturbed the super-enhancer topology associated with PAX8 gene locus, resulting in epigenetic downregulation of PAX8 transcripts and related targets. HDAC antagonists efficaciously suppressed ovarian tumor growth and spreading as single agents, and exerted synergistic effects in combination with standard chemotherapy. These findings provide mechanistic and therapeutic insights for PAX8-addicted ovarian cancer. More generally, our analytic and experimental approach represents an expandible paradigm for identifying and targeting lineage-survival oncogenes in diverse human malignancies.

Combined Transcriptome and Proteome Analysis of RpoS Regulon Reveals Its Role in Spoilage Potential of Pseudomonas fluorescens

Microbial contamination is considered the main cause of food spoilage. Pseudomonas fluorescens is a typical spoilage bacterium contributing to a large extent to the spoilage process of proteinaceous foods. RpoS is known as an alternative sigma factor controlling stress resistance and virulence in many pathogens. Our previous work revealed that RpoS contributes to the spoilage activities of P. fluorescens by regulating resistance to different stress conditions, extracellular acylated homoserine lactone (AHL) levels, extracellular protease and total volatile basic nitrogen (TVB-N) production. However, RpoS-dependent genes in P. fluorescens remained undefined. RNA-seq transcriptomics analysis combined with quantitative proteomics analysis based on multiplexed isobaric tandem mass tag (TMT) labeling was performed in the P. fluorescens wild-type strain UK4 and its derivative carrying an rpoS mutation. A total of 375 differentially expressed coding sequences (DECs) and 212 differentially expressed proteins (DEPs) were identified. The DECs were further verified by qRT-PCR. The combined transcriptome and proteome analyses revealed the involvement of this regulator in several cellular processes, mainly including polysaccharide metabolism, intracellular secretion, extracellular structures, cell wall biogenesis, stress responses, and amino acid and biogenic amine metabolism, which may contribute to the biofilm formation, stress resistance, and spoilage activities of P. fluorescens. Moreover, we indeed observed that RpoS contributed to the production of the macrocolony biofilm's matrix. Our results provide insights into the regulatory network of RpoS and expand the knowledge about the role of RpoS in the functioning of P. fluorescens in food spoilage.

Identification, Functional Characterization, and Regulon Prediction of the Zinc Uptake Regulator (zur) of Bacillus anthracis - An Insight Into the Zinc Homeostasis of the Pathogen

Zinc has an abounding occurrence in the prokaryotes and plays paramount roles including catalytic, structural, and regulatory. Zinc uptake regulator (Zur), a Fur family transcriptional regulator, is connoted in maintaining zinc homeostasis in the pathogenic bacteria by binding to zinc and regulating the genes involved in zinc uptake and mobilization. Zinc homeostasis has been marginally scrutinized in Bacillus anthracis, the top-rated bio-terror agent, with no decipherment of the role of Zur. Of the three Fur family regulators in B. anthracis, BAS4181 is annotated as a zinc-specific transcriptional regulator. This annotation was further substantiated by our stringent computational and experimental analyses. The residues critical for zinc and DNA binding were delineated by homology modeling and sequence/structure analysis. ba zur existed as a part of a three-gene operon. Purified BaZur prodigiously existed in the dimeric form, indicated by size exclusion chromatography and blue native-polyacrylamide gel electrophoresis (PAGE). Computational and manual strategies were employed to decipher the putative regulon of ba zur, comprising of 11 genes, controlled by six promoters, each harboring at least one Zur box. The DNA binding capability of the purified BaZur to the upstream regions of the ba zur operon, yciC, rpmG, znuA, and genes encoding a GTPase cobalamine synthesis protein and a permease was ascertained by electrophoretic mobility shift assays. The regulon genes, implicated in zinc uptake and mobilization, were mostly negatively regulated by BaZur. The ba zur expression was downregulated upon exposure of cells to an excess of zinc. Conversely, it exhibited a marked upregulation under N, N, N', N'-Tetrakis (2-pyridylmethyl) ethylenediamine (TPEN) mediated zinc-depleted environment, adding credence to its negative autoregulation. Moreover, an increase in the transcript levels of the regulon genes znuA, rpmG, and yciC upon exposure of cells to TPEN connoted their role in combating hypo-zincemic conditions by bringing about zinc uptake and mobilization. Thus, this study functionally characterizes Zur of B. anthracis and elucidates its role in maintaining zinc homeostasis.

The Recruitment Model of Metabolic Evolution: Jasmonate-Responsive Transcription Factors and a Conceptual Model for the Evolution of Metabolic Pathways

Plants produce a vast array of structurally diverse specialized metabolites with various biological activities, including medicinal alkaloids and terpenoids, from relatively simple precursors through a series of enzymatic steps. Massive metabolic flow through these pathways usually depends on the transcriptional coordination of a large set of metabolic, transport, and regulatory genes known as a regulon. The coexpression of genes involved in certain metabolic pathways in a wide range of developmental and environmental contexts has been investigated through transcriptomic analysis, which has been successfully exploited to mine the genes involved in various metabolic processes. Transcription factors are DNA-binding proteins that recognize relatively short sequences known as cis-regulatory elements residing in the promoter regions of target genes. Transcription factors have positive or negative effects on gene transcription mediated by RNA polymerase II. Evolutionarily conserved transcription factors of the APETALA2/ETHYLENE RESPONSE FACTOR (AP2/ERF) and basic helix-loop-helix (bHLH) families have been identified as jasmonate (JA)-responsive transcriptional regulators of unrelated specialized pathways in distinct plant lineages. Here, I review the current knowledge and propose a conceptual model for the evolution of metabolic pathways, termed "recruitment model of metabolic evolution." According to this model, structural genes are repeatedly recruited into regulons under the control of conserved transcription factors through the generation of cognate cis-regulatory elements in the promoters of these genes. This leads to the adjustment of catalytic activities that improve metabolic flow through newly established passages.

RECTA: Regulon Identification Based on Comparative Genomics and Transcriptomics Analysis

Regulons, which serve as co-regulated gene groups contributing to the transcriptional regulation of microbial genomes, have the potential to aid in understanding of underlying regulatory mechanisms. In this study, we designed a novel computational pipeline, regulon identification based on comparative genomics and transcriptomics analysis (RECTA), for regulon prediction related to the gene regulatory network under certain conditions. To demonstrate the effectiveness of this tool, we implemented RECTA on Lactococcus lactis MG1363 data to elucidate acid-response regulons. A total of 51 regulons were identified, 14 of which have computational-verified significance. Among these 14 regulons, five of them were computationally predicted to be connected with acid stress response. Validated by literature, 33 genes in Lactococcus lactis MG1363 were found to have orthologous genes which were associated with six regulons. An acid response related regulatory network was constructed, involving two trans-membrane proteins, eight regulons (llrA, llrC, hllA, ccpA, NHP6A, rcfB, regulons #8 and #39), nine functional modules, and 33 genes with orthologous genes known to be associated with acid stress. The predicted response pathways could serve as promising candidates for better acid tolerance engineering in Lactococcus lactis. Our RECTA pipeline provides an effective way to construct a reliable gene regulatory network through regulon elucidation, and has strong application power and can be effectively applied to other bacterial genomes where the elucidation of the transcriptional regulation network is needed.

Gene Enrichment Analysis Reveals Major Regulators of Mycobacterium tuberculosis Gene Expression in Two Models of Antibiotic Tolerance

The development of antibiotic tolerance is believed to be a major factor in the lengthy duration of current tuberculosis therapies. In the current study, we have modeled antibiotic tolerance in vitro by exposing Mycobacterium tuberculosis to two distinct stress conditions: progressive hypoxia and nutrient starvation [phosphate-buffered saline (PBS)]. We then studied the bacterial transcriptional response using RNA-seq and employed a bioinformatics approach to identify important transcriptional regulators, which was facilitated by a novel Regulon Enrichment Test (RET). A total of 17 transcription factor (TF) regulons were enriched in the hypoxia gene set and 16 regulons were enriched in the nutrient starvation, with 12 regulons enriched in both conditions. Using the same approach to analyze previously published gene expression datasets, we found that three M. tuberculosis regulons (Rv0023, SigH, and Crp) were commonly induced in both stress conditions and were also among the regulons enriched in our data. These regulators are worthy of further study to determine their potential role in the development and maintenance of antibiotic tolerance in M. tuberculosis following stress exposure.

Transcriptome Profiling Reveals the EanI/R Quorum Sensing Regulon in Pantoea Ananatis LMG 2665T

Pantoea ananatis LMG 2665T synthesizes and utilizes acyl homoserine lactones (AHLs) for signalling. The complete set of genes regulated by the EanI/R quorum sensing (QS) system in this strain is still not fully known. In this study, RNA-sequencing (RNA-seq) was used to identify the EanI/R regulon in LMG 2665T. Pairwise comparisons of LMG 2665T in the absence of AHLs (Optical density (OD)600 = 0.2) and in the presence of AHLs (OD600 = 0.5) were performed. Additionally, pairwise comparisons of LMG 2665T and its QS mutant at OD600 = 0.5 were undertaken. In total, 608 genes were differentially expressed between LMG 2665T at OD600 = 0.5 versus the same strain at OD600 = 0.2 and 701 genes were differentially expressed between LMG 2665T versus its QS mutant at OD600 = 0.5. A total of 196 genes were commonly differentially expressed between the two approaches. These constituted approximately 4.5% of the whole transcriptome under the experimental conditions used in this study. The RNA-seq data was validated by reverse transcriptase quantitative polymerase chain reaction (RT-qPCR). Genes found to be regulated by EanI/R QS were those coding for redox sensing, metabolism, flagella formation, flagella dependent motility, cell adhesion, biofilm formation, regulators, transport, chemotaxis, methyl accepting proteins, membrane proteins, cell wall synthesis, stress response and a large number of hypothetical proteins. The results of this study give insight into the genes that are regulated by the EanI/R system in LMG 2665T. Functional characterization of the QS regulated genes in LMG 2665T could assist in the formulation of control strategies for this plant pathogen.

Comprehensive Molecular Characterization of Muscle-Invasive Bladder Cancer

We report a comprehensive analysis of 412 muscle-invasive bladder cancers characterized by multiple TCGA analytical platforms. Fifty-eight genes were significantly mutated, and the overall mutational load was associated with APOBEC-signature mutagenesis. Clustering by mutation signature identified a high-mutation subset with 75% 5-year survival. mRNA expression clustering refined prior clustering analyses and identified a poor-survival "neuronal" subtype in which the majority of tumors lacked small cell or neuroendocrine histology. Clustering by mRNA, long non-coding RNA (lncRNA), and miRNA expression converged to identify subsets with differential epithelial-mesenchymal transition status, carcinoma in situ scores, histologic features, and survival. Our analyses identified 5 expression subtypes that may stratify response to different treatments.

P-Body Purification Reveals the Condensation of Repressed mRNA Regulons

Within cells, soluble RNPs can switch states to coassemble and condense into liquid or solid bodies. Although these phase transitions have been reconstituted in vitro, for endogenous bodies the diversity of the components, the specificity of the interaction networks, and the function of the coassemblies remain to be characterized. Here, by developing a fluorescence-activated particle sorting (FAPS) method to purify cytosolic processing bodies (P-bodies) from human epithelial cells, we identified hundreds of proteins and thousands of mRNAs that structure a dense network of interactions, separating P-body from non-P-body RNPs. mRNAs segregating into P-bodies are translationally repressed, but not decayed, and this repression explains part of the poor genome-wide correlation between RNA and protein abundance. P-bodies condense thousands of mRNAs that strikingly encode regulatory processes. Thus, we uncovered how P-bodies, by condensing and segregating repressed mRNAs, provide a physical substrate for the coordinated regulation of posttranscriptional mRNA regulons.

Regulation of gene expression by translation factor eIF5A: Hypusine-modified eIF5A enhances nonsense-mediated mRNA decay in human cells

Nonsense-mediated mRNA decay (NMD) couples protein synthesis to mRNA turnover. It eliminates defective transcripts and controls the abundance of certain normal mRNAs. Our study establishes a connection between NMD and the translation factor eIF5A (eukaryotic initiation factor 5A) in human cells. eIF5A modulates the synthesis of groups of proteins (the eIF5A regulon), and undergoes a distinctive two-step post-translational modification (hypusination) catalyzed by deoxyhypusine synthase and deoxyhypusine hydroxylase. We show that expression of NMD-susceptible constructs was increased by depletion of the major eIF5A isoform, eIF5A1. NMD was also attenuated when hypusination was inhibited by RNA interference with either of the two eIF5A modifying enzymes, or by treatment with the drugs ciclopirox or deferiprone which inhibit deoxyhypusine hydroxylase. Transcriptome analysis by RNA-Seq identified human genes whose expression is coordinately regulated by eIF5A1, its modifying enzymes, and the pivotal NMD factor, Upf1. Transcripts encoding components of the translation system were highly represented, including some encoding ribosomal proteins controlled by alternative splicing coupled to NMD (AS-NMD). Our findings extend and strengthen the association of eIF5A with NMD, previously inferred in yeast, and show that hypusination is important for this function of human eIF5A. In addition, they advance drug-mediated NMD suppression as a therapeutic opportunity for nonsense-associated diseases. We propose that regulation of mRNA stability contributes to eIF5A's role in selective gene expression.

A model for evolution and regulation of nicotine biosynthesis regulon in tobacco

In tobacco, the defense alkaloid nicotine is produced in roots and accumulates mainly in leaves. Signaling mediated by jasmonates (JAs) induces the formation of nicotine via a series of structural genes that constitute a regulon and are coordinated by JA-responsive transcription factors of the ethylene response factor (ERF) family. Early steps in the pyrrolidine and pyridine biosynthesis pathways likely arose through duplication of the polyamine and nicotinamide adenine dinucleotide (NAD) biosynthetic pathways, respectively, followed by recruitment of duplicated primary metabolic genes into the nicotine biosynthesis regulon. Transcriptional regulation of nicotine biosynthesis by ERF and cooperatively-acting MYC2 transcription factors is implied by the frequency of cognate cis-regulatory elements for these factors in the promoter regions of the downstream structural genes. Indeed, a mutant tobacco with low nicotine content was found to have a large chromosomal deletion in a cluster of closely related ERF genes at the nicotine-controlling NICOTINE2 (NIC2) locus.

Novel DNA Binding and Regulatory Activities for σ54 (RpoN) in Salmonella enterica Serovar Typhimurium 14028s

The variable sigma (σ) subunit of the bacterial RNA polymerase (RNAP) holoenzyme, which is responsible for promoter specificity and open complex formation, plays a strategic role in the response to environmental changes. Salmonella enterica serovar Typhimurium utilizes the housekeeping σ⁷⁰ and five alternative sigma factors, including σ⁵⁴ The σ⁵⁴-RNAP differs from other σ-RNAP holoenzymes in that it forms a stable closed complex with the promoter and requires ATP hydrolysis by an activated cognate bacterial enhancer binding protein (bEBP) to transition to an open complex and initiate transcription. In S. Typhimurium, σ⁵⁴-dependent promoters normally respond to one of 13 different bEBPs, each of which is activated under a specific growth condition. Here, we utilized a constitutively active, promiscuous bEBP to perform a genome-wide identification of σ⁵⁴-RNAP DNA binding sites and the transcriptome of the σ⁵⁴ regulon of S. Typhimurium. The position and context of many of the identified σ⁵⁴ RNAP DNA binding sites suggest regulatory roles for σ⁵⁴-RNAP that connect the σ⁵⁴ regulon to regulons of other σ factors to provide a dynamic response to rapidly changing environmental conditions.IMPORTANCE The alternative sigma factor σ⁵⁴ (RpoN) is required for expression of genes involved in processes with significance in agriculture, bioenergy production, bioremediation, and host-microbe interactions. The characterization of the σ⁵⁴ regulon of the versatile pathogen S. Typhimurium has expanded our understanding of the scope of the σ⁵⁴ regulon and how it links to other σ regulons within the complex regulatory network for gene expression in bacteria.

Regulation of Sporangium Formation by BldD in the Rare Actinomycete Actinoplanes missouriensis

The rare actinomycete Actinoplanes missouriensis forms sporangia, including hundreds of flagellated spores that start swimming as zoospores after their release. Under conditions suitable for vegetative growth, zoospores stop swimming and germinate. A comparative proteome analysis between zoospores and germinating cells identified 15 proteins that were produced in larger amounts in germinating cells. They include an orthologue of BldD (herein named AmBldD [BldD of A. missouriensis]), which is a transcriptional regulator involved in morphological development and secondary metabolism in Streptomyces AmBldD was detected in mycelia during vegetative growth but was barely detected in mycelia during the sporangium-forming phase, in spite of the constant transcription of AmbldD throughout growth. An AmbldD mutant started to form sporangia much earlier than the wild-type strain, and the resulting sporangia were morphologically abnormal. Recombinant AmBldD bound a palindromic sequence, the AmBldD box, located upstream from AmbldD 3',5'-Cyclic di-GMP significantly enhanced the in vitro DNA-binding ability of AmBldD. A chromatin immunoprecipitation-sequencing analysis and an in silico search for AmBldD boxes revealed that AmBldD bound 346 genomic loci that contained the 19-bp inverted repeat 5'-NN(G/A)TNACN(C/G)N(G/C)NGTNA(C/T)NN-3' as the consensus AmBldD-binding sequence. The transcriptional analysis of 27 selected AmBldD target gene candidates indicated that AmBldD should repress 12 of the 27 genes, including bldM, ssgB, whiD, ddbA, and wblA orthologues. These genes are involved in morphological development in Streptomyces coelicolor A3(2). Thus, AmBldD is a global transcriptional regulator that seems to repress the transcription of tens of genes during vegetative growth, some of which are likely to be required for sporangium formation.IMPORTANCE The rare actinomycete Actinoplanes missouriensis undergoes complex morphological differentiation, including sporangium formation. However, almost no molecular biological studies have been conducted on this bacterium. BldD is a key global regulator involved in the morphological development of streptomycetes. BldD orthologues are highly conserved among sporulating actinomycetes, but no BldD orthologues, except one in Saccharopolyspora erythraea, have been studied outside the streptomycetes. Here, it was revealed that the BldD orthologue AmBldD is essential for normal developmental processes in A. missouriensis The AmBldD regulon seems to be different from the BldD regulon in Streptomyces coelicolor A3(2), but they share four genes that are involved in morphological differentiation in S. coelicolor A3(2).

Quorum sensing and RsaM regulons of the rice pathogen Pseudomonas fuscovaginae

Pseudomonas fuscovaginae (Pfv) is an emerging plant pathogen causing sheath brown rot in rice, as well as diseases in other gramineae food crops including maize, sorghum and wheat. Pfv possesses two conserved N-acyl homoserine lactone (AHL) quorum sensing (QS) systems called PfvI/R and PfsI/R, which are repressed by RsaL and RsaM, respectively. The two systems are not hierarchically organized and are involved in plant virulence. In this study the AHL QS PfsI/R, PfvI/R and RsaM regulons were determined by transcriptomic analysis. The PfsI/R system regulates 98 genes, whereas 26 genes are regulated by the PfvI/R AHL QS system; only two genes are regulated by both systems. RsaM, on the other hand, regulates over 400 genes: 206 are negatively regulated and 260 are positively regulated. More than half of the genes controlled by the PfsI/R system and 65 % by the PfvI/R system are also part of the RsaM regulon; this is due to RsaM being involved in the regulation of both systems. It is concluded that the two QS systems regulate a unique set of genes and that RsaM is a global regulator mediating the expression of different genes through the two QS systems as well as genes independently of QS.

Conservation and diversity of the IrrE/DdrO-controlled radiation response in radiation-resistant Deinococcus bacteria

The extreme radiation resistance of Deinococcus bacteria requires the radiation‐stimulated cleavage of protein DdrO by a specific metalloprotease called IrrE. DdrO is the repressor of a predicted radiation/desiccation response (RDR) regulon, composed of radiation‐induced genes having a conserved DNA motif (RDRM) in their promoter regions. Here, we showed that addition of zinc ions to purified apo‐IrrE, and short exposure of Deinococcus cells to zinc ions, resulted in cleavage of DdrO in vitro and in vivo, respectively. Binding of IrrE to RDRM‐containing DNA or interaction of IrrE with DNA‐bound DdrO was not observed. The data are in line with IrrE being a zinc peptidase, and indicate that increased zinc availability, caused by oxidative stress, triggers the in vivo cleavage of DdrO unbound to DNA. Transcriptomics and proteomics of Deinococcus deserti confirmed the IrrE‐dependent regulation of predicted RDR regulon genes and also revealed additional members of this regulon. Comparative analysis showed that the RDR regulon is largely well conserved in Deinococcus species, but also showed diversity in the regulon composition. Notably, several RDR genes with an important role in radiation resistance in Deinococcus radiodurans, for example pprA, are not conserved in some other radiation‐resistant Deinococcus species.

Components of the SNARE-containing regulon are co-regulated in root cells undergoing defense

The term regulon has been coined in the genetic model plant Arabidopsis thaliana, denoting a structural and physiological defense apparatus defined genetically through the identification of the penetration (pen) mutants. The regulon is composed partially by the soluble N-ethylmaleimide-sensitive fusion protein attachment protein receptor (SNARE) syntaxin PEN1. PEN1 has homology to a Saccharomyces cerevisae gene that regulates a Secretion (Sec) protein, Suppressor of Sec 1 (Sso1p). The regulon is also composed of the β-glucosidase (PEN2) and an ATP binding cassette (ABC) transporter (PEN3). While important in inhibiting pathogen infection, limited observations have been made regarding the transcriptional regulation of regulon genes until now. Experiments made using the model agricultural Glycine max (soybean) have identified co-regulated gene expression of regulon components. The results explain the observation of hundreds of genes expressed specifically in the root cells undergoing the natural process of defense. Data regarding additional G. max genes functioning within the context of the regulon are presented here, including Sec 14, Sec 4 and Sec 23. Other examined G. max homologs of membrane fusion genes include an endosomal bromo domain-containing protein1 (Bro1), syntaxin6 (SYP6), SYP131, SYP71, SYP8, Bet1, coatomer epsilon (ϵ-COP), a coatomer zeta (ζ-COP) paralog and an ER to Golgi component (ERGIC) protein. Furthermore, the effectiveness of biochemical pathways that would function within the context of the regulon ave been examined, including xyloglucan xylosyltransferase (XXT), reticuline oxidase (RO) and galactinol synthase (GS). The experiments have unveiled the importance of the regulon during defense in the root and show how the deposition of callose relates to the process.

The late-annotated small ORF LSO1 is a target gene of the iron regulon of Saccharomyces cerevisiae

We have identified a new downstream target gene of the Aft1/2‐regulated iron regulon in budding yeast Saccharomyces cerevisiae, the late‐annotated small open reading frame LSO1. LSO1 transcript is among the most highly induced from a transcriptome analysis of a fet3‐1 mutant grown in the presence of the iron chelator bathophenanthrolinedisulfonic acid. LSO1 has a paralog, LSO2, which is constitutively expressed and not affected by iron availability. In contrast, we find that the LSO1 promoter region contains three consensus binding sites for the Aft1/2 transcription factors and that an LSO1‐lacZ reporter is highly induced under low‐iron conditions in a Aft1‐dependent manner. The expression patterns of the Lso1 and Lso2 proteins mirror those of their mRNAs. Both proteins are localized to the nucleus and cytoplasm, but become more cytoplasmic upon iron deprivation consistent with a role in iron transport. LSO1 and LSO2 appear to play overlapping roles in the cellular response to iron starvation since single lso1 and lso2 mutants are sensitive to iron deprivation and this sensitivity is exacerbated when both genes are deleted.

Identification of sdiA-regulated genes in a mouse commensal strain of Enterobacter cloacae

Many bacteria determine their population density using quorum sensing. The most intensively studied mechanism of quorum sensing utilizes proteins of the LuxI family to synthesize a signaling molecule of the acylhomoserine lactone (AHL) type, and a protein of the LuxR family to bind AHL and regulate transcription. Genes regulated by quorum sensing often encode functions that are most effective when a group of bacteria are working cooperatively (e.g., luminescence, biofilm formation, host interactions). Bacteria in the Escherichia, Salmonella, Klebsiella, and Enterobacter genera do not encode an AHL synthase but they do encode an AHL receptor of the LuxR family, SdiA. Instead of detecting their own AHL synthesis, these organisms use SdiA to detect the AHLs synthesized by other bacterial species. In this study, we used a genetic screen to identify AHL-responsive genes in a commensal Enterobacter cloacae strain that was isolated from a laboratory mouse. The genes include a putative type VI secretion system, copA (a copper transporter), and fepE (extends O-antigen chain length). A new transposon mutagenesis strategy and suicide vectors were used to construct an sdiA mutant of E. cloacae. The AHL-responsiveness of all fusions was entirely sdiA-dependent, although some genes were regulated by sdiA in the absence of AHL.

The mRNAs associated to a zinc finger protein from Trypanosoma cruzi shift during stress conditions

Trypanosome gene expression is regulated almost exclusively at the posttranscriptional level, through mRNA stability, storage and degradation. Here, we characterize the ribonucleoprotein complex (mRNPs) corresponding to the zinc finger protein TcZC3H39 from T. cruzi comparing cells growing in normal conditions and under nutritional stress. The nutritional stress is a key step during T. cruzi differentiation from epimastigote form to human infective metacyclic trypomastigote form. The mechanisms by which the stress, altogether with other stimuli, triggers differentiation is not well understood. This work aims to characterize the TcZC3H39 protein during stress response. Using cells cultured in normal and stress conditions, we observed a dynamic change in TcZC3H39 granule distribution, which appeared broader in stressed epimastigotes. The protein core of the TcZC3H39-mRNP is composed of ribosomes, translation factors and RBPs. The TcZC3H39-mRNP could act sequestering highly expressed mRNAs and their associated ribosomes, potentially slowing translation in stress conditions. A shift were observed in the mRNAs associated with TcZC3H39: the number of targets in unstressed epimastigotes was smaller than that in stressed parasites, with no clear functional clustering in normal conditions. By contrast, in stressed parasites, the targets of TcZC3H39 were mRNAs encoding ribosomal proteins and a remarkable enrichment in mRNAs for the cytochrome c complex (COX), highly expressed mRNAs in the replicative form. This identification of a new component of RNA granules in T. cruzi, the TcZC3H39 protein, provides new insight into the mechanisms involved in parasite stress responses and the regulation of gene expression during T. cruzi differentiation.

Regulation of galactomannan biosynthesis in coffee seeds

The seed of Coffea arabica accumulates large amounts of cell wall storage polysaccharides (CWSPs) of the mannan family in the cell walls of the endosperm. The variability induced by the growing environment and extensive pairwise correlation analysis with stringent significance thresholds was used to investigate transcript-transcript and transcript-metabolite relationships among 26 sugar-related genes, and the amount of CWSPs and seven soluble low molecular weight carbohydrates in the developing coffee endosperm. A dense module of nine quantitatively co-expressed genes was detected at the mid-developmental stage when CWSPs accumulate. This module included the five genes of the core galactomannan synthetic machinery, namely genes coding for the enzymes needed to assemble the mannan backbone (mannan synthase, ManS), and genes that introduce the galactosyl side chains (galactosyltransferase, GMGT), modulate the post-depositional degree of galactose substitution (α-galactosidase), and produce the nucleotide sugar building blocks GDP-mannose and UDP-galactose (mannose-1P guanyltransferase and UDP-glucose 4'-epimerase, respectively). The amount of CWSPs stored in the endosperm at the onset of their accumulation was primarily and quantitatively modulated at the transcriptional level (i.e. positively correlated with the expression level of these key galactomannan biosynthetic genes). This analysis also suggests a role for sorbitol and raffinose family oligosaccharides as transient auxiliary sources of building blocks for galactomannan synthesis. Finally, a microarray-based analysis of the developing seed transcriptome revealed that all genes of the core galactomannan synthesis machinery grouped in a single cluster of 209 co-expressed genes. Analysis of the gene composition of this cluster revealed remarkable functional coherence and identified transcription factors that putatively control galactomannan biosynthesis in coffee.

Comparative genomics and functional analysis of rhamnose catabolic pathways and regulons in bacteria

L-rhamnose (L-Rha) is a deoxy-hexose sugar commonly found in nature. L-Rha catabolic pathways were previously characterized in various bacteria including Escherichia coli. Nevertheless, homology searches failed to recognize all the genes for the complete L-Rha utilization pathways in diverse microbial species involved in biomass decomposition. Moreover, the regulatory mechanisms of L-Rha catabolism have remained unclear in most species. A comparative genomics approach was used to reconstruct the L-Rha catabolic pathways and transcriptional regulons in the phyla Actinobacteria, Bacteroidetes, Chloroflexi, Firmicutes, Proteobacteria, and Thermotogae. The reconstructed pathways include multiple novel enzymes and transporters involved in the utilization of L-Rha and L-Rha-containing polymers. Large-scale regulon inference using bioinformatics revealed remarkable variations in transcriptional regulators for L-Rha utilization genes among bacteria. A novel bifunctional enzyme, L-rhamnulose-phosphate aldolase (RhaE) fused to L-lactaldehyde dehydrogenase (RhaW), which is not homologous to previously characterized L-Rha catabolic enzymes, was identified in diverse bacteria including Chloroflexi, Bacilli, and Alphaproteobacteria. By using in vitro biochemical assays we validated both enzymatic activities of the purified recombinant RhaEW proteins from Chloroflexus aurantiacus and Bacillus subtilis. Another novel enzyme of the L-Rha catabolism, L-lactaldehyde reductase (RhaZ), was identified in Gammaproteobacteria and experimentally validated by in vitro enzymatic assays using the recombinant protein from Salmonella typhimurium. C. aurantiacus induced transcription of the predicted L-Rha utilization genes when L-Rha was present in the growth medium and consumed L-Rha from the medium. This study provided comprehensive insights to L-Rha catabolism and its regulation in diverse Bacteria.