The repertoire of DNA-binding transcriptional regulators in Escherichia coli K-12

The primary transcriptome of Neisseria meningitidis and its interaction with the RNA chaperone Hfq

Neisseria meningitidis is a human commensal that can also cause life-threatening meningitis and septicemia. Despite growing evidence for RNA-based regulation in meningococci, their transcriptome structure and output of regulatory small RNAs (sRNAs) are incompletely understood. Using dRNA-seq, we have mapped at single-nucleotide resolution the primary transcriptome of N. meningitidis strain 8013. Annotation of 1625 transcriptional start sites defines transcription units for most protein-coding genes but also reveals a paucity of classical σ70-type promoters, suggesting the existence of activators that compensate for the lack of −35 consensus sequences in N. meningitidis. The transcriptome maps also reveal 65 candidate sRNAs, a third of which were validated by northern blot analysis. Immunoprecipitation with the RNA chaperone Hfq drafts an unexpectedly large post-transcriptional regulatory network in this organism, comprising 23 sRNAs and hundreds of potential mRNA targets. Based on this data, using a newly developed gfp reporter system we validate an Hfq-dependent mRNA repression of the putative colonization factor PrpB by the two trans-acting sRNAs RcoF1/2. Our genome-wide RNA compendium will allow for a better understanding of meningococcal transcriptome organization and riboregulation with implications for colonization of the human nasopharynx.

The LacI-Family Transcription Factor, RbsR, Is a Pleiotropic Regulator of Motility, Virulence, Siderophore and Antibiotic Production, Gas Vesicle Morphogenesis and Flotation in Serratia

Gas vesicles (GVs) are proteinaceous, gas-filled organelles used by some bacteria to enable upward movement into favorable air/liquid interfaces in aquatic environments. Serratia sp. ATCC39006 (S39006) was the first enterobacterium discovered to produce GVs naturally. The regulation of GV assembly in this host is complex and part of a wider regulatory network affecting various phenotypes, including antibiotic biosynthesis. To identify new regulators of GVs, a comprehensive mutant library containing 71,000 insertion mutants was generated by random transposon mutagenesis and 311 putative GV-defective mutants identified. Three of these mutants were found to have a transposon inserted in a LacI family transcription regulator gene (rbsR) of the putative ribose operon. Each of these rbsR mutants was GV-defective; no GVs were visible by phase contrast microscopy (PCM) or transmission electron microscopy (TEM). GV deficiency was caused by the reduction of gvpA1 and gvrA transcription (the first genes of the two contiguous operons in the GV gene locus). Our results also showed that a mutation in rbsR was highly pleiotropic; the production of two secondary metabolites (carbapenem and prodigiosin antibiotics) was abolished. Interestingly, the intrinsic resistance to the carbapenem antibiotic was not affected by the rbsR mutation. In addition, the production of a siderophore, cellulase and plant virulence was reduced in the mutant, whereas it exhibited increased swimming and swarming motility. The RbsR protein was predicted to bind to regions upstream of at least 18 genes in S39006 including rbsD (the first gene of the ribose operon) and gvrA. Electrophoretic mobility shift assays (EMSA) confirmed that RbsR bound to DNA sequences upstream of rbsD, but not gvrA. The results of this study indicate that RbsR is a global regulator that affects the modulation of GV biogenesis, but also with complex pleiotropic physiological impacts in S39006.

The whole set of the constitutive promoters recognized by four minor sigma subunits of Escherichia coli RNA polymerase

The promoter selectivity of Escherichia coli RNA polymerase (RNAP) is determined by the sigma subunit. The model prokaryote Escherichia coli K-12 contains seven species of the sigma subunit, each recognizing a specific set of promoters. For identification of the "constitutive promoters" that are recognized by each RNAP holoenzyme alone in the absence of other supporting factors, we have performed the genomic SELEX screening in vitro for their binding sites along the E. coli K-12 W3110 genome using each of the reconstituted RNAP holoenzymes and a collection of genome DNA segments of E. coli K-12. The whole set of constitutive promoters for each RNAP holoenzyme was then estimated based on the location of RNAP-binding sites. The first successful screening of the constitutive promoters was achieved for RpoD (σ70), the principal sigma for transcription of growth-related genes. As an extension, we performed in this study the screening of constitutive promoters for four minor sigma subunits, stationary-phase specific RpoS (σ38), heat-shock specific RpoH (σ32), flagellar-chemotaxis specific RpoF (σ28) and extra-cytoplasmic stress-response RpoE (σ24). The total number of constitutive promoters were: 129~179 for RpoS; 101~142 for RpoH; 34~41 for RpoF; and 77~106 for RpoE. The list of constitutive promoters were compared with that of known promoters identified in vivo under various conditions and using varieties of E. coli strains, altogether allowing the estimation of "inducible promoters" in the presence of additional supporting factors.

Inactivation of Transcriptional Regulators during Within-Household Evolution of Escherichia coli

We analyzed the within-household evolution of two household-associated Escherichia coli strains from pandemic clonal group ST131-H30, using isolates recovered from five individuals within two families, each of which had a distinct strain. Family 1's strain was represented by a urine isolate from the index patient (older sister) with recurrent cystitis and a blood isolate from her younger sister with fatal urosepsis. Family 2's strain was represented by a urine isolate from the index patient (father) with pyelonephritis and renal abscesses, blood and kidney drainage isolates from the daughter with emphysematous pyelonephritis, and urine and fecal isolates from the mother with cystitis. Collectively, the several variants of each family's strain had accumulated a total of 8 (family 1) and 39 (family 2) point mutations; no two isolates were identical. Of the 47 total mutations, 36 resulted in amino acid changes or truncation of coded proteins. Fourteen such mutations (39%) targeted genes encoding transcriptional regulators, and 9 (25%) involved DNA-binding transcription factors (TFs), which significantly exceeded the relative contribution of TF genes to the isolates' genomes (∼6%). At least one-half of the transcriptional regulator mutations were inactivating, based on phenotypic and/or transcriptional analysis. In particular, inactivating mutations in the global regulator LrhA (repressor of type 1 fimbriae and flagella) occurred in the blood isolates from both households and increased the virulence of E. coli strains in a murine sepsis model. The results indicate that E. coli undergoes adaptive evolution between and/or within hosts, generating subpopulations with distinctive phenotypes and virulence potential.IMPORTANCE The clonal evolution of bacterial strains associated with interhost transmission is poorly understood. We characterized the genome sequences of clonal descendants of two Escherichia coli strains, recovered at different time points from multiple individuals within two households who had different types of urinary tract infection. We found evidence that the E. coli strains underwent extensive mutational diversification between and within these individuals, driven disproportionately by inactivation of transcriptional regulators. In urosepsis isolates, the mutations observed in the global regulator LrhA increased bacterial virulence in a murine sepsis model. Our findings help in understanding the adaptive dynamics and strategies of E. coli during short-term natural evolution.

Defined chromosome structure in the genome-reduced bacterium Mycoplasma pneumoniae

DNA-binding proteins are central regulators of chromosome organization; however, in genome-reduced bacteria their diversity is largely diminished. Whether the chromosomes of such bacteria adopt defined three-dimensional structures remains unexplored. Here we combine Hi-C and super-resolution microscopy to determine the structure of the Mycoplasma pneumoniae chromosome at a 10 kb resolution. We find a defined structure, with a global symmetry between two arms that connect opposite poles, one bearing the chromosomal Ori and the other the midpoint. Analysis of local structures at a 3 kb resolution indicates that the chromosome is organized into domains ranging from 15 to 33 kb. We provide evidence that genes within the same domain tend to be co-regulated, suggesting that chromosome organization influences transcriptional regulation, and that supercoiling regulates local organization. This study extends the current understanding of bacterial genome organization and demonstrates that a defined chromosomal structure is a universal feature of living systems.

Tailor-made transcriptional biosensors for optimizing microbial cell factories

Monitoring cellular behavior and eventually properly adapting cellular processes is key to handle the enormous complexity of today's metabolic engineering questions. Hence, transcriptional biosensors bear the potential to augment and accelerate current metabolic engineering strategies, catalyzing vital advances in industrial biotechnology. The development of such transcriptional biosensors typically starts with exploring nature's richness. Hence, in a first part, the transcriptional biosensor architecture and the various modi operandi are briefly discussed, as well as experimental and computational methods and relevant ontologies to search for natural transcription factors and their corresponding binding sites. In the second part of this review, various engineering approaches are reviewed to tune the main characteristics of these (natural) transcriptional biosensors, i.e., the response curve and ligand specificity, in view of specific industrial biotechnology applications, which is illustrated using success stories of transcriptional biosensor engineering.

Functional architecture and global properties of the Corynebacterium glutamicum regulatory network: Novel insights from a dataset with a high genomic coverage

Corynebacterium glutamicum is a Gram-positive, anaerobic, rod-shaped soil bacterium able to grow on a diversity of carbon sources like sugars and organic acids. It is a biotechnological relevant organism because of its highly efficient ability to biosynthesize amino acids, such as l-glutamic acid and l-lysine. Here, we reconstructed the most complete C. glutamicum regulatory network to date and comprehensively analyzed its global organizational properties, systems-level features and functional architecture. Our analyses show the tremendous power of Abasy Atlas to study the functional organization of regulatory networks. We created two models of the C. glutamicum regulatory network: all-evidences (containing both weak and strong supported interactions, genomic coverage=73%) and strongly-supported (only accounting for strongly supported evidences, genomic coverage=71%). Using state-of-the-art methodologies, we prove that power-law behaviors truly govern the connectivity and clustering coefficient distributions. We found a non-previously reported circuit motif that we named complex feed-forward motif. We highlighted the importance of feedback loops for the functional architecture, beyond whether they are statistically over-represented or not in the network. We show that the previously reported top-down approach is inadequate to infer the hierarchy governing a regulatory network because feedback bridges different hierarchical layers, and the top-down approach disregards the presence of intermodular genes shaping the integration layer. Our findings all together further support a diamond-shaped, three-layered hierarchy exhibiting some feedback between processing and coordination layers, which is shaped by four classes of systems-level elements: global regulators, locally autonomous modules, basal machinery and intermodular genes.

Molecular and structural considerations of TF-DNA binding for the generation of biologically meaningful and accurate phylogenetic footprinting analysis: the LysR-type transcriptional regulator family as a study model

BACKGROUND

The goal of most programs developed to find transcription factor binding sites (TFBSs) is the identification of discrete sequence motifs that are significantly over-represented in a given set of sequences where a transcription factor (TF) is expected to bind. These programs assume that the nucleotide conservation of a specific motif is indicative of a selective pressure required for the recognition of a TF for its corresponding TFBS. Despite their extensive use, the accuracies reached with these programs remain low. In many cases, true TFBSs are excluded from the identification process, especially when they correspond to low-affinity but important binding sites of regulatory systems.

RESULTS

We developed a computational protocol based on molecular and structural criteria to perform biologically meaningful and accurate phylogenetic footprinting analyses. Our protocol considers fundamental aspects of the TF-DNA binding process, such as: i) the active homodimeric conformations of TFs that impose symmetric structures on the TFBSs, ii) the cooperative binding of TFs, iii) the effects of the presence or absence of co-inducers, iv) the proximity between two TFBSs or one TFBS and a promoter that leads to very long spurious motifs, v) the presence of AT-rich sequences not recognized by the TF but that are required for DNA flexibility, and vi) the dynamic order in which the different binding events take place to determine a regulatory response (i.e., activation or repression). In our protocol, the abovementioned criteria were used to analyze a profile of consensus motifs generated from canonical Phylogenetic Footprinting Analyses using a set of analysis windows of incremental sizes. To evaluate the performance of our protocol, we analyzed six members of the LysR-type TF family in Gammaproteobacteria.

CONCLUSIONS

The identification of TFBSs based exclusively on the significance of the over-representation of motifs in a set of sequences might lead to inaccurate results. The consideration of different molecular and structural properties of the regulatory systems benefits the identification of TFBSs and enables the development of elaborate, biologically meaningful and precise regulatory models that offer a more integrated view of the dynamics of the regulatory process of transcription.

Carriage of Extended-Spectrum Beta-Lactamase-Plasmids Does Not Reduce Fitness but Enhances Virulence in Some Strains of Pandemic E. coli Lineages

Pathogenic ESBL-producing E. coli lineages occur frequently worldwide, not only in a human health context but in animals and the environment, also in settings with low antimicrobial pressures. This study investigated the fitness costs of ESBL-plasmids and their influence on chromosomally encoded features associated with virulence, such as those involved in the planktonic and sessile behaviors of ST131 and ST648 E. coli. ESBL-plasmid-carrying wild-type E. coli strains, their corresponding ESBL-plasmid-“cured” variants (PCV), and complementary ESBL-carrying transformants were comparatively analyzed using growth curves, Omnilog® phenotype microarray (PM) assays, macrocolony and biofilm formation, swimming motility, and RNA sequence analysis. Growth curves and PM results pointed toward similar growth and metabolic behaviors among the strains. Phenotypic differences in some strains were detected, including enhanced curli fimbriae and/or cellulose production as well as a reduced swimming capacity of some ESBL-carrying strains, as compared to their respective PCVs. RNA sequencing mostly confirmed the phenotypic results, suggesting that the chromosomally encoded csgD pathway is a key factor involved. These results contradict the hypothesis that ESBL-plasmid-carriage leads to a fitness loss in ESBL-carrying strains. Instead, the results indicate an influence of some ESBL-plasmids on chromosomally encoded features associated with virulence in some E. coli strains. In conclusion, apart from antibiotic resistance selective advantages, ESBL-plasmid-carriage may also lead to enhanced virulence or adaption to specific habitats in some strains of pandemic ESBL-producing E. coli lineages.

Transcription profile of Escherichia coli: genomic SELEX search for regulatory targets of transcription factors

Bacterial genomes are transcribed by DNA-dependent RNA polymerase (RNAP), which achieves gene selectivity through interaction with sigma factors that recognize promoters, and transcription factors (TFs) that control the activity and specificity of RNAP holoenzyme. To understand the molecular mechanisms of transcriptional regulation, the identification of regulatory targets is needed for all these factors. We then performed genomic SELEX screenings of targets under the control of each sigma factor and each TF. Here we describe the assembly of 156 SELEX patterns of a total of 116 TFs performed in the presence and absence of effector ligands. The results reveal several novel concepts: (i) each TF regulates more targets than hitherto recognized; (ii) each promoter is regulated by more TFs than hitherto recognized; and (iii) the binding sites of some TFs are located within operons and even inside open reading frames. The binding sites of a set of global regulators, including cAMP receptor protein, LeuO and Lrp, overlap with those of the silencer H-NS, suggesting that certain global regulators play an anti-silencing role. To facilitate sharing of these accumulated SELEX datasets with the research community, we compiled a database, ‘Transcription Profile of Escherichia coli’ (www.shigen.nig.ac.jp/ecoli/tec/).

Transcription Factors Exhibit Differential Conservation in Bacteria with Reduced Genomes

The description of transcriptional regulatory networks has been pivotal in the understanding of operating principles under which organisms respond and adapt to varying conditions. While the study of the topology and dynamics of these networks has been the subject of considerable work, the investigation of the evolution of their topology, as a result of the adaptation of organisms to different environmental conditions, has received little attention. In this work, we study the evolution of transcriptional regulatory networks in bacteria from a genome reduction perspective, which manifests itself as the loss of genes at different degrees. We used the transcriptional regulatory network of Escherichia coli as a reference to compare 113 smaller, phylogenetically-related γ-proteobacteria, including 19 genomes of symbionts. We found that the type of regulatory action exerted by transcription factors, as genomes get progressively smaller, correlates well with their degree of conservation, with dual regulators being more conserved than repressors and activators in conditions of extreme reduction. In addition, we found that the preponderant conservation of dual regulators might be due to their role as both global regulators and nucleoid-associated proteins. We summarize our results in a conceptual model of how each TF type is gradually lost as genomes become smaller and give a rationale for the order in which this phenomenon occurs.

Genetic Regulation of Yersinia pestis

Y. pestis exhibits dramatically different traits of pathogenicity and transmission, albeit their close genetic relationship with its ancestor-Y. pseudotuberculosis, a self-limiting gastroenteric pathogen. Y. pestis is evolved into a deadly pathogen and transmitted to mammals and/or human beings by infected flea biting or directly contacting with the infected animals. Various kinds of environmental changes are implicated into its complex life cycle and pathogenesis. Dynamic regulation of gene expression is critical for environmental adaptation or survival, primarily reflected by genetic regulation mediated by transcriptional factors and small regulatory RNAs at the transcriptional and posttranscriptional level, respectively. The effects of genetic regulation have been shown to profoundly influence Y. pestis physiology and pathogenesis such as stress resistance, biofilm formation, intracellular survival, and replication. In this chapter, we mainly summarize the progresses on popular methods of genetic regulation and on regulatory patterns and consequences of many key transcriptional and posttranscriptional regulators, with a particular emphasis on how genetic regulation influences the biofilm and virulence of Y. pestis.

Collation and analyses of DNA-binding protein domain families from sequence and structural databanks

DNA-protein interactions govern several high fidelity cellular processes like DNA-replication, transcription, DNA repair, etc. Proteins that have the ability to recognise and bind DNA sequences can be classified either according to their DNA-binding motif or based on the sequence of the target nucleotides. We have collated the DNA-binding families by integrating information from both protein sequence family and structural databases. This resulted in a dataset of 1057 DNA-binding protein domain families. Their family properties (the number of members, percent identity distribution and length of members) and domain architectures were examined. Further, sequence domain families were mapped to structures in the protein databank (PDB) and the protein domain structure classification database (SCOP). The DNA-binding families, with no structural information, were clustered together into potential superfamilies based on sequence associations. On the basis of functions attributed to DNA-binding protein folds, we observe that a majority of the DNA-binding proteins follow divergent evolution. This study can serve as a basis for annotation and distribution of DNA-binding proteins in genome(s) of interest. The entire collated set of DNA-binding protein domains is available for download as Hidden Markov Models.

Novel gene sets improve set-level classification of prokaryotic gene expression data

Background

Set-level classification of gene expression data has received significant attention recently. In this setting, high-dimensional vectors of features corresponding to genes are converted into lower-dimensional vectors of features corresponding to biologically interpretable gene sets. The dimensionality reduction brings the promise of a decreased risk of overfitting, potentially resulting in improved accuracy of the learned classifiers. However, recent empirical research has not confirmed this expectation. Here we hypothesize that the reported unfavorable classification results in the set-level framework were due to the adoption of unsuitable gene sets defined typically on the basis of the Gene ontology and the KEGG database of metabolic networks. We explore an alternative approach to defining gene sets, based on regulatory interactions, which we expect to collect genes with more correlated expression. We hypothesize that such more correlated gene sets will enable to learn more accurate classifiers.

Methods

We define two families of gene sets using information on regulatory interactions, and evaluate them on phenotype-classification tasks using public prokaryotic gene expression data sets. From each of the two gene-set families, we first select the best-performing subtype. The two selected subtypes are then evaluated on independent (testing) data sets against state-of-the-art gene sets and against the conventional gene-level approach.

Results

The novel gene sets are indeed more correlated than the conventional ones, and lead to significantly more accurate classifiers. The novel gene sets are indeed more correlated than the conventional ones, and lead to significantly more accurate classifiers.

Conclusion

Novel gene sets defined on the basis of regulatory interactions improve set-level classification of gene expression data. The experimental scripts and other material needed to reproduce the experiments are available at http://ida.felk.cvut.cz/novelgenesets.tar.gz.

Electronic supplementary material

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

GntR Family of Bacterial Transcription Factors and Their DNA Binding Motifs: Structure, Positioning and Co-Evolution

The GntR family of transcription factors (TFs) is a large group of proteins present in diverse bacteria and regulating various biological processes. Here we use the comparative genomics approach to reconstruct regulons and identify binding motifs of regulators from three subfamilies of the GntR family, FadR, HutC, and YtrA. Using these data, we attempt to predict DNA-protein contacts by analyzing correlations between binding motifs in DNA and amino acid sequences of TFs. We identify pairs of positions with high correlation between amino acids and nucleotides for FadR, HutC, and YtrA subfamilies and show that the most predicted DNA-protein interactions are quite similar in all subfamilies and conform well to the experimentally identified contacts formed by FadR from E. coli and AraR from B. subtilis. The most frequent predicted contacts in the analyzed subfamilies are Arg-G, Asn-A, Asp-C. We also analyze the divergon structure and preferred site positions relative to regulated genes in the FadR and HutC subfamilies. A single site in a divergon usually regulates both operons and is approximately in the middle of the intergenic area. Double sites are either involved in the co-operative regulation of both operons and then are in the center of the intergenic area, or each site in the pair independently regulates its own operon and tends to be near it. We also identify additional candidate TF-binding boxes near palindromic binding sites of TFs from the FadR, HutC, and YtrA subfamilies, which may play role in the binding of additional TF-subunits.

The functional landscape bound to the transcription factors of Escherichia coli K-12

Motivated by the experimental evidences accumulated in the last ten years and based on information deposited in RegulonDB, literature look up, and sequence analysis, we analyze the repertoire of 304 DNA-binding Transcription factors (TFs) in Escherichia coli K-12. These regulators were grouped in 78 evolutionary families and are regulating almost half of the total genes in this bacterium. In structural terms, 60% of TFs are composed by two-domains, 30% are monodomain, and 10% three- and four-structural domains. As previously noticed, the most abundant DNA-binding domain corresponds to the winged helix-turn-helix, with few alternative DNA-binding structures, resembling the hypothesis of successful protein structures with the emergence of new ones at low scales. In summary, we identified and described the characteristics associated to the DNA-binding TF in E. coli K-12. We also identified twelve functional modules based on a co-regulated gene matrix. Finally, diverse regulons were predicted based on direct associations between the TFs and potential regulated genes. This analysis should increase our knowledge about the gene regulation in the bacterium E. coli K-12, and provide more additional clues for comprehensive modelling of transcriptional regulatory networks in other bacteria.

Regulation of Streptococcus mutans PTS Bio by the transcriptional repressor NigR

Streptococcus mutans is implicated in human dental caries, and the carbohydrate metabolism of this organism plays an important role in the formation of this disease. Carbohydrate transport and metabolism are essential for the survival of S. mutans in the oral cavity. It is known that a unique phosphoenolpyruvate-sugar phosphotransferase system PTS(B) (io) of S. mutans UA159 is expressed in sucrose-grown biofilms (Mol Oral Microbiol 28: 2013; 114). In this study we analyzed the transcriptional regulation of the operon (O(B) (io) ) encoding the PTS(B) (io) and showed that it was repressed by NigR, a LacI-like transcriptional regulator. Using electro-mobility shift assay, we described two operators to which NigR bound with different affinities. We also identified the transcriptional start site and showed that one of the operators overlaps with the promoter and presumably represses initiation of transcription. Mutational analyses revealed the key nucleotides in the operators required for high-affinity binding of NigR. PTS(B) (io) is expressed in S. mutans biofilms so understanding its regulation may provide improved strategies for caries treatment and prevention.

The Subtleties and Contrasts of the LeuO Regulator in Salmonella Typhi: Implications in the Immune Response

Salmonella are facultative intracellular pathogens. Salmonella infection occurs mainly by expression of two Salmonella pathogenicity Islands (SPI-1 and SPI-2). SPI-1 encodes transcriptional factors that participate in the expression of virulence factors encoded in the island. However, there are transcriptional factors encoded outside the island that also participate in the expression of SPI-1-encoded genes. Upon infection, bacteria are capable of avoiding the host immune response with several strategies that involve several virulence factors under the control of transcriptional regulators. Interestingly, LeuO a transcriptional global regulator which is encoded outside of any SPI, is proposed to be part of a complex regulatory network that involves expression of several genes that help bacteria to survive stress conditions and, also, induces the expression of porins that have been shown to be immunogens and can thus be considered as antigenic candidates for acellular vaccines. Hence, the understanding of the LeuO regulon implies a role of bacterial genetic regulation in determining the host immune response.

Prediction of DNA binding motifs from 3D models of transcription factors; identifying TLX3 regulated genes

Proper cell functioning depends on the precise spatio-temporal expression of its genetic material. Gene expression is controlled to a great extent by sequence-specific transcription factors (TFs). Our current knowledge on where and how TFs bind and associate to regulate gene expression is incomplete. A structure-based computational algorithm (TF2DNA) is developed to identify binding specificities of TFs. The method constructs homology models of TFs bound to DNA and assesses the relative binding affinity for all possible DNA sequences using a knowledge-based potential, after optimization in a molecular mechanics force field. TF2DNA predictions were benchmarked against experimentally determined binding motifs. Success rates range from 45% to 81% and primarily depend on the sequence identity of aligned target sequences and template structures, TF2DNA was used to predict 1321 motifs for 1825 putative human TF proteins, facilitating the reconstruction of most of the human gene regulatory network. As an illustration, the predicted DNA binding site for the poorly characterized T-cell leukemia homeobox 3 (TLX3) TF was confirmed with gel shift assay experiments. TLX3 motif searches in human promoter regions identified a group of genes enriched in functions relating to hematopoiesis, tissue morphology, endocrine system and connective tissue development and function.

Expression levels of transcription factors in Escherichia coli: growth phase- and growth condition-dependent variation of 90 regulators from six families

The expression pattern of the genome in Escherichia coli is controlled by regulating the utilization of a limited number of RNA polymerases between a total of 4600 genes on its genome. The distribution pattern of RNA polymerase on the genome changes after two steps of protein–protein interaction with seven sigma subunits and about 300 transcription factors (TFs). Based on a systematic search for the regulation target promoters recognized by each TF, we propose two novel concepts: each TF regulates a number of target promoters; and each promoter is regulated by many TFs. In parallel, attempts have been made to determine the intracellular concentrations of all TFs using two systems: quantitative immunoblot analysis using TF-specific antibodies; and reporter assay of TF promoter activities. The direct measurement of TF protein level has so far been published for a set of 60 regulators with known functions. This study describes the determination of growth phase-dependent expression levels of 90 TFs using the reporter assay system. The translational fusion vector was constructed from the TF promoter sequence including an N-terminal proximal TF segment and the reporter GFP. At the beginning of cell growth, high-level expression was observed only for a small number of TFs. In the exponential phase, approximately 80 % TFs are expressed, but the expressed TF species change upon transfer to the stationary phase. Significant changes in the pattern of TF expression were observed between aerobic and anaerobic conditions. The list of intracellular levels of TFs provides further understanding to the transcription regulation of the E. coli genome under various stressful conditions.

Intracellular concentrations of 65 species of transcription factors with known regulatory functions in Escherichia coli

The expression pattern of the Escherichia coli genome is controlled in part by regulating the utilization of a limited number of RNA polymerases among a total of its approximately 4,600 genes. The distribution pattern of RNA polymerase changes from modulation of two types of protein-protein interactions: the interaction of core RNA polymerase with seven species of the sigma subunit for differential promoter recognition and the interaction of RNA polymerase holoenzyme with about 300 different species of transcription factors (TFs) with regulatory functions. We have been involved in the systematic search for the target promoters recognized by each sigma factor and each TF using the newly developed Genomic SELEX system. In parallel, we developed the promoter-specific (PS)-TF screening system for identification of the whole set of TFs involved in regulation of each promoter. Understanding the regulation of genome transcription also requires knowing the intracellular concentrations of the sigma subunits and TFs under various growth conditions. This report describes the intracellular levels of 65 species of TF with known function in E. coli K-12 W3110 at various phases of cell growth and at various temperatures. The list of intracellular concentrations of the sigma factors and TFs provides a community resource for understanding the transcription regulation of E. coli under various stressful conditions in nature.

The TetR family of regulators

The most common prokaryotic signal transduction mechanisms are the one-component systems in which a single polypeptide contains both a sensory domain and a DNA-binding domain. Among the >20 classes of one-component systems, the TetR family of regulators (TFRs) are widely associated with antibiotic resistance and the regulation of genes encoding small-molecule exporters. However, TFRs play a much broader role, controlling genes involved in metabolism, antibiotic production, quorum sensing, and many other aspects of prokaryotic physiology. There are several well-established model systems for understanding these important proteins, and structural studies have begun to unveil the mechanisms by which they bind DNA and recognize small-molecule ligands. The sequences for more than 200,000 TFRs are available in the public databases, and genomics studies are identifying their target genes. Three-dimensional structures have been solved for close to 200 TFRs. Comparison of these structures reveals a common overall architecture of nine conserved α helices. The most important open question concerning TFR biology is the nature and diversity of their ligands and how these relate to the biochemical processes under their control.

Metabolic flux analysis for Escherichia coli by flux balance analysis

Conventional metabolic flux analysis (MFA) of Escherichia coli wild type and of pathway gene knockout mutants cultivated under anaerobic condition is explained in detail in this chapter. To place the MFA results into the context of the literature, the regulation of central carbon metabolism in terms of catabolite regulation by the phosphotransferase system (PTS) and the response to oxygen limitations via global regulators is reviewed. The effects of gene deletions such as pflA, pta, ppc, pykF, adhE, and ldhA on the metabolic network are presented. Moreover, for the pflA mutant the effects of various carbon sources were quantified. The chapter thereby contributes to the discussion of metabolic network function and the design of microbial cell factories.

Analysis of the transcriptional regulator GlpR, promoter elements, and posttranscriptional processing involved in fructose-induced activation of the phosphoenolpyruvate-dependent sugar phosphotransferase system in Haloferax mediterranei

Among all known archaeal strains, the phosphoenolpyruvate-dependent phosphotransferase system (PTS) for fructose utilization is used primarily by haloarchaea, which thrive in hypersaline environments, whereas the molecular details of the regulation of the archaeal PTS under fructose induction remain unclear. In this study, we present a comprehensive examination of the regulatory mechanism of the fructose PTS in the haloarchaeon Haloferax mediterranei. With gene knockout and complementation, microarray analysis, and chromatin immunoprecipitation-quantitative PCR (ChIP-qPCR), we revealed that GlpR is the indispensable activator, which specifically binds to the PTS promoter (PPTS) during fructose induction. Further promoter-scanning mutation indicated that three sites located upstream of the H. mediterranei PPTS, which are conserved in most haloarchaeal PPTSs, are involved in this induction. Interestingly, two PTS transcripts (named T8 and T17) with different lengths of 5' untranslated region (UTR) were observed, and promoter or 5' UTR swap experiments indicated that the shorter 5' UTR was most likely generated from the longer one. Notably, the translation efficiency of the transcript with this shorter 5' UTR was significantly higher and the ratio of T8 (with the shorter 5' UTR) to T17 increased during fructose induction, implying that a posttranscriptional mechanism is also involved in PTS activation. With these insights into the molecular regulation of the haloarchaeal PTS, we have proposed a working model for haloarchaea in response to environmental fructose.

SdrA, a new DeoR family regulator involved in Streptomyces avermitilis morphological development and antibiotic production

The SAV3339 (SdrA) protein of Streptomyces avermitilis, a member of the DeoR family of regulators, was assessed to determine its in vivo function by gene knockdown through the use of cis-encoded noncoding RNA and knockout of the sdrA gene. These analyses revealed that SdrA represents another class of Streptomyces regulator that controls morphological development and antibiotic production.

Quantitative influence of macromolecular crowding on gene regulation kinetics

We introduce macromolecular crowding quantitatively into the model for kinetics of gene regulation in Escherichia coli. We analyse and compute the specific-site searching time for 180 known transcription factors (TFs) regulating 1300 operons. The time is between 160 s (e.g. for SoxS M_w = 12.91 kDa) and 1550 s (e.g. for PepA₆ of M_w = 329.28 kDa). Diffusion coefficients for one-dimensional sliding are between for large proteins up to for small monomers or dimers. Three-dimensional diffusion coefficients in the cytoplasm are 2 orders of magnitude larger than 1D sliding coefficients, nevertheless the sliding enhances the binding rates of TF to specific sites by 1–2 orders of magnitude. The latter effect is due to ubiquitous non-specific binding. We compare the model to experimental data for LacI repressor and find that non-specific binding of the protein to DNA is activation- and not diffusion-limited. We show that the target location rate by LacI repressor is optimized with respect to microscopic rate constant for association to non-specific sites on DNA. We analyse the effect of oligomerization of TFs and DNA looping effects on searching kinetics. We show that optimal searching strategy depends on TF abundance.

How reliable is the linear noise approximation of gene regulatory networks?

Background

The linear noise approximation (LNA) is commonly used to predict how noise is regulated and exploited at the cellular level. These predictions are exact for reaction networks composed exclusively of first order reactions or for networks involving bimolecular reactions and large numbers of molecules. It is however well known that gene regulation involves bimolecular interactions with molecule numbers as small as a single copy of a particular gene. It is therefore questionable how reliable are the LNA predictions for these systems.

Results

We implement in the software package intrinsic Noise Analyzer (iNA), a system size expansion based method which calculates the mean concentrations and the variances of the fluctuations to an order of accuracy higher than the LNA. We then use iNA to explore the parametric dependence of the Fano factors and of the coefficients of variation of the mRNA and protein fluctuations in models of genetic networks involving nonlinear protein degradation, post-transcriptional, post-translational and negative feedback regulation. We find that the LNA can significantly underestimate the amplitude and period of noise-induced oscillations in genetic oscillators. We also identify cases where the LNA predicts that noise levels can be optimized by tuning a bimolecular rate constant whereas our method shows that no such regulation is possible. All our results are confirmed by stochastic simulations.

Conclusion

The software iNA allows the investigation of parameter regimes where the LNA fares well and where it does not. We have shown that the parametric dependence of the coefficients of variation and Fano factors for common gene regulatory networks is better described by including terms of higher order than LNA in the system size expansion. This analysis is considerably faster than stochastic simulations due to the extensive ensemble averaging needed to obtain statistically meaningful results. Hence iNA is well suited for performing computationally efficient and quantitative studies of intrinsic noise in gene regulatory networks.

Current challenges in bacterial transcriptomics

Over the past decade or so, dramatic developments in our ability to experimentally determine the content and function of genomes have taken place. In particular, next-generation sequencing technologies are now inspiring a new understanding of bacterial transcriptomes on a global scale. In bacterial cells, whole-transcriptome studies have not received attention, owing to the general view that bacterial genomes are simple. However, several recent RNA sequencing results are revealing unexpected levels of complexity in bacterial transcriptomes, indicating that the transcribed regions of genomes are much larger and complex than previously anticipated. In particular, these data show a wide array of small RNAs, antisense RNAs, and alternative transcripts. Here, we review how current transcriptomics are now revolutionizing our understanding of the complexity and regulation of bacterial transcriptomes.

Identification of Salmonella enterica serovar Pullorum antigenic determinants expressed in vivo

Salmonella enterica serovar Pullorum affecting poultry causes pullorum disease and results in severe economic loss in the poultry industry. Currently, it remains a major threat in countries with poor poultry surveillance and no efficient control measures. As S. Pullorum could induce strong humoral immune responses, we applied an immunoscreening technique, the in vivo-induced antigen technology (IVIAT), to identify immunogenic bacterial proteins expressed or upregulated during S. Pullorum infection. Convalescent-phase sera from chickens infected with S. Pullorum were pooled, adsorbed against antigens expressed in vitro, and used to screen an S. Pullorum genomic expression library. Forty-five proteins were screened out, and their functions were implicated in molecular biosynthesis and degradation, transport, metabolism, regulation, cell wall synthesis and antibiotic resistance, environmental adaptation, or putative functions. In addition, 11 of these 45 genes were assessed for their differential expression by quantitative real-time reverse transcription-PCR (RT-PCR), revealing that 9 of 11 genes were upregulated to different degrees under in vivo conditions, especially the regulator of virulence determinants, phoQ. Then, four in vivo-induced proteins (ShdA, PhoQ, Cse3, and PbpC) were tested for their immunoreactivity in 28 clinical serum samples from chickens infected with S. Pullorum. The rate of detection of antibodies against ShdA reached 82% and was the highest among these proteins. ShdA is a host colonization factor known to be upregulated in vivo and related to the persistence of S. Typhimurium in the intestine. Furthermore, these antigens identified by IVIAT warrant further evaluation for their contributions to pathogenesis, and more potential roles, such as diagnostic, therapeutic, and preventive uses, need to be developed in future studies.

Characterization of the transcriptional stimulatory properties of the Pseudomonas putida RapA protein

RNA polymerase-associated factors can significantly affect its performance at specific promoters. Here we identified a Pseudomonas putida RNA polymerases-associated protein as a homolog of Escherichia coli RapA. We found that P. putida RapA stimulates the transcription from promoters dependent on a variety of σ-factors (σ(70), σ(S), σ(54), σ(32), σ(E)) in vitro. The level of stimulation varied from 2- to 10-fold, with the maximal effect observed with the σ(E)-dependent PhtrA promoter. Stimulation by RapA was apparent in the multi-round reactions and was modulated by salt concentration in vitro. However, in contrast to findings with E. coli RapA, P. putida RapA-mediated stimulation of transcription was also evident using linear templates. These properties of P. putida RapA were apparent using either E. coli- or P. putida-derived RNA polymerases. Analysis of individual steps of transcription revealed that P. putida RapA enhances the stability of competitor-resistant open-complexes formed by RNA polymerase at promoters. In vivo, P. putida RapA can complement the inhibitory effect of high salt on growth of an E. coli RapA null strain. However, a P. putida RapA null mutant was not sensitive to high salt. The in vivo effects of lack of RapA were only detectable for the σ(E)-PhtrA promoter where the RapA-deficiency resulted in lower activity. The presented characteristics of P. putida RapA indicate that its functions may extend beyond a role in facilitating RNA polymerase recycling to include a role in transcription initiation efficiency.

Genome, functional gene annotation, and nuclear transformation of the heterokont oleaginous alga Nannochloropsis oceanica CCMP1779

Unicellular marine algae have promise for providing sustainable and scalable biofuel feedstocks, although no single species has emerged as a preferred organism. Moreover, adequate molecular and genetic resources prerequisite for the rational engineering of marine algal feedstocks are lacking for most candidate species. Heterokonts of the genus Nannochloropsis naturally have high cellular oil content and are already in use for industrial production of high-value lipid products. First success in applying reverse genetics by targeted gene replacement makes Nannochloropsis oceanica an attractive model to investigate the cell and molecular biology and biochemistry of this fascinating organism group. Here we present the assembly of the 28.7 Mb genome of N. oceanica CCMP1779. RNA sequencing data from nitrogen-replete and nitrogen-depleted growth conditions support a total of 11,973 genes, of which in addition to automatic annotation some were manually inspected to predict the biochemical repertoire for this organism. Among others, more than 100 genes putatively related to lipid metabolism, 114 predicted transcription factors, and 109 transcriptional regulators were annotated. Comparison of the N. oceanica CCMP1779 gene repertoire with the recently published N. gaditana genome identified 2,649 genes likely specific to N. oceanica CCMP1779. Many of these N. oceanica-specific genes have putative orthologs in other species or are supported by transcriptional evidence. However, because similarity-based annotations are limited, functions of most of these species-specific genes remain unknown. Aside from the genome sequence and its analysis, protocols for the transformation of N. oceanica CCMP1779 are provided. The availability of genomic and transcriptomic data for Nannochloropsis oceanica CCMP1779, along with efficient transformation protocols, provides a blueprint for future detailed gene functional analysis and genetic engineering of Nannochloropsis species by a growing academic community focused on this genus.

YdfH identified as a repressor of rspA by the use of reduced genome Escherichia coli MGF-01

We have reported the construction of 1 Mb reduced genome Escherichia coli MGF-01 by a 28-step operation. This time, transcriptome analysis of MGF-01 was performed. Although the transcriptome profiles of the exponential phase in parental strain W3110red were well-conserved in MGF-01, the rspAB operon was highly expressed. A LacZ reporter assay of a series of stepwise deletion strains prepared in the course of MGF-01 construction indicated that rspA was highly expressed after the 5th step. Further analysis indicated that Δ29, one of the deleted regions at the 5th step, relates to an increase in rspA expression, and that transcriptional regulator ydfH, in the Δ29 region, is responsible for the expression of rspA, gel shift assay indicated that YdfH bound directly to the upstream region of rspA. Based on these results, it was concluded that YdfH is a transcriptional repressor of the rspAB operon.

Hierarchical control on polyene macrolide biosynthesis: PimR modulates pimaricin production via the PAS-LuxR transcriptional activator PimM

Control of polyene macrolide production in Streptomyces natalensis is mediated by the transcriptional activator PimR. This regulator combines an N-terminal domain corresponding to the Streptomyces antibiotic regulatory protein (SARP) family of transcriptional activators with a C-terminal half homologous to guanylate cyclases and large ATP-binding regulators of the LuxR family. The PimR SARP domain (PimR(SARP)) was expressed in Escherichia coli as a glutathione S-transferase (GST)-fused protein. Electrophoretic mobility shift assays showed that GST-PimR(SARP) binds a single target, the intergenic region between the regulatory genes pimR and pimMs in the pimaricin cluster. The PimR(SARP)-binding site was investigated by DNaseI protection studies, revealing that it contains three heptameric direct repeats adjusting to the consensus 5'-CGGCAAG-3'. Transcription start points of pimM and pimR promoters were identified by 5'-RACE, revealing that unlike other SARPs, PimR(SARP) does not interact with the -35 region of its target promoter. Quantitative transcriptional analysis of these regulatory genes on mutants on each of them has allowed the identification of the pimM promoter as the transcriptional target for PimR. Furthermore, the constitutive expression of pimM restored pimaricin production in a pimaricin-deficient strain carrying a deletion mutant of pimR. These results reveal that PimR exerts its positive effect on pimaricin production by controlling pimM expression level, a regulator whose gene product activates transcription from eight different promoters of pimaricin structural genes directly.

Mycobacterium tuberculosis RbpA protein is a new type of transcriptional activator that stabilizes the σ A-containing RNA polymerase holoenzyme

RbpA is an RNA polymerase (RNAP)-binding protein whose presence increases the tolerance levels of Mycobacteria to the first-line anti-tuberculosis drug rifampicin by an unknown mechanism. Here, we show that the role of Mycobacterium tuberculosis RbpA in resistance is indirect because it does not affect the sensitivity of RNAP to rifampicin while it stimulates transcription controlled by the housekeeping σ^A-factor. The transcription regulated by the stress-related σ^F was not affected by RbpA. The binding site of RbpA maps to the RNAP β subunit Sandwich-Barrel Hybrid Motif, which has not previously been described as an activator target and does not overlap the rifampicin binding site. Our data suggest that RbpA modifies the structure of the core RNAP, increases its affinity for σ^A and facilitates the assembly of the transcriptionally competent promoter complexes. We propose that RbpA is an essential partner which advantages σ^A competitiveness for core RNAP binding with respect to the alternative σ factors. The RbpA-driven stimulation of the housekeeping gene expression may help Mycobacteria to tolerate high rifampicin levels and to adapt to the stress conditions during infection.

Prokaryotic genome regulation: a revolutionary paradigm

After determination of the whole genome sequence, the research frontier of bacterial molecular genetics has shifted to reveal the genome regulation under stressful conditions in nature. The gene selectivity of RNA polymerase is modulated after interaction with two groups of regulatory proteins, 7 sigma factors and 300 transcription factors. For identification of regulation targets of transcription factors in Escherichia coli, we have developed Genomic SELEX system and subjected to screening the binding sites of these factors on the genome. The number of regulation targets by a single transcription factor was more than those hitherto recognized, ranging up to hundreds of promoters. The number of transcription factors involved in regulation of a single promoter also increased to as many as 30 regulators. The multi-target transcription factors and the multi-factor promoters were assembled into complex networks of transcription regulation. The most complex network was identified in the regulation cascades of transcription of two master regulators for planktonic growth and biofilm formation.

Extracting regulatory networks of Escherichia coli from RegulonDB

RegulonDB contains the largest and currently best-known data set on transcriptional regulation in a single free-living organism, that of Escherichia coli K-12 (Gama-Castro et al. Nucleic Acids Res 36:D120-D124, 2008). This organized knowledge has been the gold standard for the implementation of bioinformatic predictive methods on gene regulation in bacteria (Collado-Vides et al. J Bacteriol 191:23-31, 2009). Given the complexity of different types of interactions, the difficulty of visualizing in a single figure of the whole network, and the different uses of this knowledge, we are making available different views of the genetic network. This chapter describes case studies about how to access these views, via precomputed files, web services and SQL, including sigma-gene relationships corresponding to transcription of alternative RNA polymerase holoenzyme promoters; as well as, transcription factor (TF)-genes, TF-operons, TF-TF, and TF-regulon interactions. 17.

Role of the coiled-coil tip of Escherichia coli DksA in promoter control

Escherichia coli DksA works in conjunction with the small-molecule ppGpp to regulate transcription initiation negatively or positively, depending on the identity of the promoter. DksA is in a class of transcription factors that do not bind directly to DNA such as classical repressors or activators but rather bind in the RNA polymerase (RNAP) secondary channel such as the transcription elongation factors GreA and GreB in E. coli and TFIIS in eukaryotes. We found that substitution for either of two residues in its coiled-coil tip, D74 or A76, eliminates DksA function without affecting its apparent affinity for RNAP. The properties of DksA-Gre factor chimeras indicated that the coiled-coil tip is responsible for the DksA-specific effects on open complex formation. A conservative substitution at position 74, D74E, resulted in a loss of DksA function in both negative and positive control, and an E44D substitution at the analogous position in GreA resulted in a gain of function in both negative and positive control. That a single methylene group has such an extraordinary effect on these transcription factors highlights the critical nature of the identity of coiled-coil tip interactions with RNAP for open complex formation.

Analyzing the metabolic stress response of recombinant Escherichia coli cultures expressing human interferon-beta in high cell density fed batch cultures using time course transcriptomic data

Fed batch cultures expressing recombinant interferon beta under the T7 promoter were run with different exponential feeding rates of a complex substrate and induced at varying cell densities. Post-induction profiles of the specific product formation rates showed a strong dependence on the specific growth rate with the maximum product yield obtained at 0.2 h(-1). A study of the relative transcriptomic profiles as a function of pre-induction μ was therefore done to provide insight into the role of cellular physiology in enhancing recombinant protein expression. Hierarchical clustering analysis of the significantly regulated genes allowed us to identify biologically important groups of genes which fall under specific master regulators. The groups were: rpoH, ArcB, CreB, Lrp, RelA, Fis and Hfq. The response of these regulators, which exert a feedback control on the growth and product formation rates correlated well with the expression levels obtained. Thus at the optimum pre-induction μ, the alternative sigma factors and ribosomal machinery genes did not get depressed till the 6th hour post-induction unlike at other specific growth rates, demonstrating a critical role for the genes in sustaining recombinant protein expression.

Comparative transcriptomic profile analysis of fed-batch cultures expressing different recombinant proteins in Escherichia coli

There is a need to elucidate the product specific features of the metabolic stress response of the host cell to the induction of recombinant protein synthesis. For this, the method of choice is transcriptomic profiling which provides a better insight into the changes taking place in complex global metabolic networks. The transcriptomic profiles of three fed-batch cultures expressing different proteins viz. recombinant human interferon-beta (rhIFN-β), Xylanase and Green Fluorescence Protein (GFP) were compared post induction. We observed a depression in the nutrient uptake and utilization pathways, which was common for all the three expressed proteins. Thus glycerol transporters and genes involved in ATP synthesis as well as aerobic respiration were severely down-regulated. On the other hand the amino acid uptake and biosynthesis genes were significantly repressed only when soluble proteins were expressed under different promoters, but not when the product was expressed as an inclusion body (IB). High level expression under the T7 promoter (rhIFN-β and xylanase) triggered the cellular degradation machinery like the osmoprotectants, proteases and mRNA degradation genes which were highly up-regulated, while this trend was not true with GFP expression under the comparatively weaker ara promoter. The design of a better host platform for recombinant protein production thus needs to take into account the specific nature of the cellular response to protein expression.

New insights on gene regulation in archaea

Archaea represent an important and vast domain of life. This cellular domain includes a large diversity of organisms characterized as prokaryotes with basal transcriptional machinery similar to eukarya. In this work we explore the most recent findings concerning the transcriptional regulatory organization in archaeal genomes since the perspective of the DNA-binding transcription factors (TFs), such as the high proportion of archaeal TFs homologous to bacteria, the apparent deficit of TFs, only comparable to the proportion of TFs in parasites or intracellular pathogenic bacteria, suggesting a deficit in this class of proteins. We discuss an appealing hypothesis to explain the apparent deficit of TFs in archaea, based on their characteristics, such as their small length sizes. The hypothesis suggests that a large fraction of these small-sized TFs could supply the deficit of TFs in archaea, by forming different combinations of monomers similar to that observed in eukaryotic transcriptional machinery, where a wide diversity of protein-protein interactions could act as mediators of regulatory feedback, indicating a chimera of bacterial and eukaryotic TFs' functionality. Finally, we discuss how global experiments can help to understand in a global context the role of TFs in these organisms.

Recognition of prokaryotic promoters based on a novel variable-window Z-curve method

Transcription is the first step in gene expression, and it is the step at which most of the regulation of expression occurs. Although sequenced prokaryotic genomes provide a wealth of information, transcriptional regulatory networks are still poorly understood using the available genomic information, largely because accurate prediction of promoters is difficult. To improve promoter recognition performance, a novel variable-window Z-curve method is developed to extract general features of prokaryotic promoters. The features are used for further classification by the partial least squares technique. To verify the prediction performance, the proposed method is applied to predict promoter fragments of two representative prokaryotic model organisms (Escherichia coli and Bacillus subtilis). Depending on the feature extraction and selection power of the proposed method, the promoter prediction accuracies are improved markedly over most existing approaches: for E. coli, the accuracies are 96.05% (σ⁷⁰ promoters, coding negative samples), 90.44% (σ⁷⁰ promoters, non-coding negative samples), 92.13% (known sigma-factor promoters, coding negative samples), 92.50% (known sigma-factor promoters, non-coding negative samples), respectively; for B. subtilis, the accuracies are 95.83% (known sigma-factor promoters, coding negative samples) and 99.09% (known sigma-factor promoters, non-coding negative samples). Additionally, being a linear technique, the computational simplicity of the proposed method makes it easy to run in a matter of minutes on ordinary personal computers or even laptops. More importantly, there is no need to optimize parameters, so it is very practical for predicting other species promoters without any prior knowledge or prior information of the statistical properties of the samples.

Identification and functional analysis of CT069 as a novel transcriptional regulator in Chlamydia

Only a small number of transcription factors have been predicted in Chlamydia spp., which are obligate intracellular bacteria that include a number of important human pathogens. We used a bioinformatics strategy to identify novel transcriptional regulators from the Chlamydia trachomatis genome by predicting proteins with the general structure and characteristic functional domains of a bacterial transcription factor. With this approach, we identified CT069 as a candidate transcription factor with sequence similarity at its C terminus to Treponema pallidum TroR. Like TroR, the gene for CT069 belongs to an operon that encodes components of a putative ABC transporter for importing divalent metal cations. However, CT069 has been annotated as YtgC because of sequence similarity at its N terminus to TroC, a transmembrane component of this metal ion transporter. Instead, CT069 appears to be a fusion protein composed of YtgC and a TroR ortholog that we have called YtgR. Although it has not been previously reported, a similar YtgC-YtgR fusion protein is predicted to be encoded by other Chlamydia spp. and several other bacteria, including Bacillus subtilis. We show that recombinant YtgR polypeptide bound specifically to an operator sequence upstream of the ytg operon and that binding was enhanced by Zn(2+). We also demonstrate that YtgR repressed transcription from the ytg promoter in a heterologous in vivo reporter assay. These results provide evidence that CT069 is a negative regulator of the ytg operon, which encodes a putative metal ion transporter in C. trachomatis.

Intramolecular signal transmission in a tetrameric repressor of the IclR family

Members of the IclR family control bacterial genes involved in a number of physiological processes. The IclR-family member TtgV crystallizes as a tetramer, with each TtgV monomer consisting of two domains--a DNA binding domain and an effector recognition domain, which are interconnected by an extended α-helix. When bound to DNA, a kink is introduced so that the extended helix is split in two α-helices (helix-4 and -5). Differential scanning calorimetry studies revealed that TtgV unfolds in a single event, suggesting that the two domains unfold cooperatively. When mutations are introduced in helix-5 that disrupt interactions between Arg98 and Glu102, the thermal unfolding of the TtgV domains becomes uncoupled without compromising effector binding. Two of these mutants (TtgVE102R and TtgVE102A) showed impaired release from target DNA, suggesting that these mutations alter signal transmission. By combining various mutants, we found that the mutations in the connecting α-helix exhibited a dominant effect over mutations in DNA binding and effector binding domains. We propose a model in which the loss of cooperativity of unfolding of TtgV reflects perturbed interdomain communication, and that the transition from the continuous to discontinuous helix may mediate interdomain communication necessary for the proper functioning of TtgV.

Deciphering the regulatory codes in bacterial genomes

Interactions between cis-regulatory elements and trans-acting factors are fundamental for cellular functions such as transcription. With the revolution in microarrays and sequencing technologies, genome-wide binding locations of trans-acting factors are being determined in large numbers. The richness of the genome-scale information has revealed that the nature of the bacterial transcriptome and regulome are considerably more complex than previously expected. In addition, the emerging view of the bacterial transcriptome is revising the concept of the operon organization of the genome. This review describes current advances in the genome-scale analysis of the interaction between cis-regulatory elements and trans-acting factors in microorganisms.

Bacterial transcriptomics: what is beyond the RNA horiz-ome?

Over the past 3 years, bacterial transcriptomics has undergone a massive revolution. Increased sequencing capacity and novel tools have made it possible to explore the bacterial transcriptome to an unprecedented depth, which has revealed that the transcriptome is more complex and dynamic than expected. Alternative transcripts within operons challenge the classic operon definition, and many small RNAs involved in the regulation of transcription, translation and pathogenesis have been discovered. Furthermore, mRNAs may localize to specific areas in the cell, and the spatial organization and dynamics of the chromosome have been shown to be important for transcription. Epigenetic modifications of DNA also affect transcription, and RNA processing affects translation. Therefore, transcription in bacteria resembles that in eukaryotes in terms of complexity more closely than was previously thought. Here we will discuss the contribution of 'omics' approaches to these discoveries as well as the possible impact that they are expected to have in the future.

Genome wide identification of Acidithiobacillus ferrooxidans (ATCC 23270) transcription factors and comparative analysis of ArsR and MerR metal regulators

Acidithiobacillus ferrooxidans is a chemolithoautotrophic acidophilic bacterium that obtains its energy from the oxidation of ferrous iron, elemental sulfur, or reduced sulfur minerals. This capability makes it of great industrial importance due to its applications in biomining. During the industrial processes, A. ferrooxidans survives to stressing circumstances in its environment, such as an extremely acidic pH and high concentration of transition metals. In order to gain insight into the organization of A. ferrooxidans regulatory networks and to provide a framework for further studies in bacterial growth under extreme conditions, we applied a genome-wide annotation procedure to identify 87 A. ferrooxidans transcription factors. We classified them into 19 families that were conserved among diverse prokaryotic phyla. Our annotation procedure revealed that A. ferrooxidans genome contains several members of the ArsR and MerR families, which are involved in metal resistance and detoxification. Analysis of their sequences revealed known and potentially new mechanism to coordinate gene-expression in response to metal availability. A. ferrooxidans inhabit some of the most metal-rich environments known, thus transcription factors identified here seem to be good candidates for functional studies in order to determine their physiological roles and to place them into A. ferrooxidans transcriptional regulatory networks.