Redefining fundamental concepts of transcription initiation in bacteria.

Despite enormous progress in understanding the fundamentals of bacterial gene regulation, our knowledge remains limited when compared with the number of bacterial genomes and regulatory systems to be discovered. Derived from a small number of initial studies, classic definitions for concepts of gene regulation have evolved as the number of characterized promoters has increased. Together with discoveries made using new technologies, this knowledge has led to revised generalizations and principles. In this Expert Recommendation, we suggest precise, updated definitions that support a logical, consistent conceptual framework of bacterial gene regulation, focusing on transcription initiation. The resulting concepts can be formalized by ontologies for computational modelling, laying the foundation for improved bioinformatics tools, knowledge-based resources and scientific communication. Thus, this work will help researchers construct better predictive models, with different formalisms, that will be useful in engineering, synthetic biology, microbiology and genetics.

Structural Insights into the Substrate Range of a Bacterial Monoamine Oxidase.

Monoamine oxidases (MAOs) play a key role in the breakdown of primary and secondary amines. In eukaryotic organisms, these enzymes are vital to the regulation of monoamine neurotransmitters and the degradation of dietary monoamines. MAOs have also been identified in prokaryotic species, although their role in these organisms is not well understood. Here, we report the biophysical and structural properties of a promiscuous, bacterial MAO from Corynebacterium ammoniagenes (caMAO). caMAO catalyzes the oxidation of a number of monoamine substrates including dopamine and norepinephrine, as well as exhibiting some activity with polyamine substrates such as cadaverine. The X-ray crystal structures of Michaelis complexes with seven substrates show that conserved hydrophobic interactions and hydrogen-bonding pattern (for polar substrates) allow the broad specificity range. The structure of caMAO identifies an unusual cysteine (Cys424) residue in the so-called "aromatic cage", which flanks the flavin isoalloxazine ring in the active site. Site-directed mutagenesis, steady-state kinetics in air-saturated buffer, and UV-vis spectroscopy revealed that Cys424 plays a role in the pH dependence and modulation of electrostatics within the caMAO active site. Notably, bioinformatic analysis shows a propensity for variation at this site within the "aromatic cage" of the flavin amine oxidase (FAO) superfamily. Structural analysis also identified the conservation of a secondary substrate inhibition site, present in a homologous member of the superfamily. Finally, genome neighborhood diagram analysis of caMAO in the context of the FAO superfamily allows us to propose potential roles for these bacterial MAOs in monoamine and polyamine degradation and catabolic pathways related to scavenging of nitrogen.

Transcription-Factor-Induced Aggregation of Biomimetic Oligonucleotide-b-Protein Micelles.

Polymeric micelles and especially those based on natural diblocks are of particular interest due to their advantageous properties in terms of molecular recognition, biocompatibility, and biodegradability. We herein report a facile and straightforward synthesis of thermoresponsive elastin-like polypeptide (ELP) and oligonucleotide (ON) diblock bioconjugates, ON-b-ELP, through copper-catalyzed azide-alkyne cycloaddition. The resulting thermosensitive diblock copolymer self-assembles above its critical micelle temperature (CMT ∼30 °C) to form colloidally stable micelles of ∼50 nm diameter. The ON-b-ELP micelles hybridize with an ON complementary strand and maintain their size and stability. Next, we describe the capacity of these micelles to bind proteins, creating more complex structures using the classic biotin-streptavidin pairing and the specific recognition between a transcription factor protein and the ON strand. In both instances, the micelles are intact, form larger structures, and retain their sensitivity to temperature.

Epigenetic competition reveals density-dependent regulation and target site plasticity of phosphorothioate epigenetics in bacteria.

Phosphorothioate (PT) DNA modifications-in which a nonbonding phosphate oxygen is replaced with sulfur-represent a widespread, horizontally transferred epigenetic system in prokaryotes and have a highly unusual property of occupying only a small fraction of available consensus sequences in a genome. Using Salmonella enterica as a model, we asked a question of fundamental importance: How do the PT-modifying DndA-E proteins select their GPSAAC/GPSTTC targets? Here, we applied innovative analytical, sequencing, and computational tools to discover a novel behavior for DNA-binding proteins: The Dnd proteins are "parked" at the G6mATC Dam methyltransferase consensus sequence instead of the expected GAAC/GTTC motif, with removal of the 6mA permitting extensive PT modification of GATC sites. This shift in modification sites further revealed a surprising constancy in the density of PT modifications across the genome. Computational analysis showed that GAAC, GTTC, and GATC share common features of DNA shape, which suggests that PT epigenetics are regulated in a density-dependent manner partly by DNA shape-driven target selection in the genome.

Development and Clinical Validation of Iso-IMRS: A Novel Diagnostic Assay for P. falciparum Malaria.

In many countries targeting malaria elimination, persistent malaria infections can have parasite loads significantly below the lower limit of detection (LLOD) of standard diagnostic techniques, making them difficult to identify and treat. The most sensitive diagnostic methods involve amplification and detection of Plasmodium DNA by polymerase chain reaction (PCR), which requires expensive thermal cycling equipment and is difficult to deploy in resource-limited settings. Isothermal DNA amplification assays have been developed, but they require complex primer design, resulting in high nonspecific amplification, and show a decrease in sensitivity than PCR methods. Here, we have used a computational approach to design a novel isothermal amplification assay with a simple primer design to amplify P. falciparum DNA with analytical sensitivity comparable to PCR. We have identified short DNA sequences repeated throughout the parasite genome to be used as primers for DNA amplification and demonstrated that these primers can be used, without modification, to isothermally amplify P. falciparum parasite DNA via strand displacement amplification. Our novel assay shows a LLOD of ∼1 parasite/µL within a 30 min amplification time. The assay was demonstrated with clinical samples using patient blood and saliva. We further characterized the assay using direct amplicon next-generation sequencing and modified the assay to work with a visual readout. The technique developed here achieves similar analytical sensitivity to current gold standard PCR assays requiring a fraction of time and resources for PCR. This highly sensitive isothermal assay can be more easily adapted to field settings, making it a potentially useful tool for malaria elimination.

RegulonDB v 10.5: tackling challenges to unify classic and high throughput knowledge of gene regulation in E. coli K-12.

RegulonDB, first published 20 years ago, is a comprehensive electronic resource about regulation of transcription initiation of Escherichia coli K-12 with decades of knowledge from classic molecular biology experiments, and recently also from high-throughput genomic methodologies. We curated the literature to keep RegulonDB up to date, and initiated curation of ChIP and gSELEX experiments. We estimate that current knowledge describes between 10% and 30% of the expected total number of transcription factor- gene regulatory interactions in E. coli. RegulonDB provides datasets for interactions for which there is no evidence that they affect expression, as well as expression datasets. We developed a proof of concept pipeline to merge binding and expression evidence to identify regulatory interactions. These datasets can be visualized in the RegulonDB JBrowse. We developed the Microbial Conditions Ontology with a controlled vocabulary for the minimal properties to reproduce an experiment, which contributes to integrate data from high throughput and classic literature. At a higher level of integration, we report Genetic Sensory-Response Units for 200 transcription factors, including their regulation at the metabolic level, and include summaries for 70 of them. Finally, we summarize our research with Natural language processing strategies to enhance our biocuration work.

A Förster Resonance Energy Transfer-Based Ratiometric Sensor with the Allosteric Transcription Factor TetR.

A recent description of an antibody-free assay is significantly extended for small molecule analytes using allosteric transcription factors (aTFs) and Förster resonance energy transfer (FRET). The FRET signal indicates the differential binding of an aTF-DNA pair with a dose-dependent response to its effector molecule, i.e., the analyte. The new sensors described here, based on the well-characterized aTF TetR, demonstrate several new features of the assay approach: 1) the generalizability of the sensors to additional aTF-DNA-analyte systems, 2) sensitivity modulation through the choice of donor fluorophore (quantum dots or fluorescent proteins, FPs), and 3) sensor tuning using aTF variants with differing aTF-DNA binding affinities. While all of these modular sensors self-assemble, the design reported here based on a recombinant aTF-FP chimera with commercially available dye-labeled DNA uses readily accessible sensor components to facilitate easy adoption of the sensing approach by the broader community.

Genomic insights into tuberculosis.

Prevalent since pre-history, human tuberculosis - caused by the pathogen Mycobacterium tuberculosis - remains a major source of death worldwide. Moreover, increasing drug resistance poses the threat of disease resurgence. However, the expanding application of genomic techniques is providing new avenues for combating this old foe. Whole-genome sequencing, comparative genomics and systems biology are generating new insights into the origins and ongoing evolution of M. tuberculosis, as well as the molecular basis for its pathogenicity. These have important implications for our perspective of the disease, development of new drugs and vaccines, and treatment of patients using existing therapeutics.

Transcription Factor Based Small-Molecule Sensing with a Rapid Cell Phone Enabled Fluorescent Bead Assay.

Recently, allosteric transcription factors (TFs) were identified as a novel class of biorecognition elements for in vitro sensing, whereby an indicator of the differential binding affinity between a TF and its cognate DNA exhibits dose-dependent responsivity to an analyte. Described is a modular bead-based biosensor design that can be applied to such TF-DNA-analyte systems. DNA-functionalized beads enable efficient mixing and spatial separation, while TF-labeled semiconductor quantum dots serve as bright fluorescent indicators of the TF-DNA bound (on bead) and unbound states. The prototype sensor for derivatives of the antibiotic tetracycline exhibits nanomolar sensitivity with visual detection of bead fluorescence. Facile changes to the sensor enable sensor response tuning without necessitating changes to the biomolecular affinities. Assay components self-assemble, and readout by eye or digital camera is possible within 5 minutes of analyte addition, making sensor use facile, rapid, and instrument-free.

The Mycobacterium tuberculosis regulatory network and hypoxia.

We have taken the first steps towards a complete reconstruction of the Mycobacterium tuberculosis regulatory network based on ChIP-Seq and combined this reconstruction with system-wide profiling of messenger RNAs, proteins, metabolites and lipids during hypoxia and re-aeration. Adaptations to hypoxia are thought to have a prominent role in M. tuberculosis pathogenesis. Using ChIP-Seq combined with expression data from the induction of the same factors, we have reconstructed a draft regulatory network based on 50 transcription factors. This network model revealed a direct interconnection between the hypoxic response, lipid catabolism, lipid anabolism and the production of cell wall lipids. As a validation of this model, in response to oxygen availability we observe substantial alterations in lipid content and changes in gene expression and metabolites in corresponding metabolic pathways. The regulatory network reveals transcription factors underlying these changes, allows us to computationally predict expression changes, and indicates that Rv0081 is a regulatory hub.

A unified resource for transcriptional regulation in Escherichia coli K-12 incorporating high-throughput-generated binding data into RegulonDB version 10.0.

BACKGROUND: Our understanding of the regulation of gene expression has benefited from the availability of high-throughput technologies that interrogate the whole genome for the binding of specific transcription factors and gene expression profiles. In the case of widely used model organisms, such as Escherichia coli K-12, the new knowledge gained from these approaches needs to be integrated with the legacy of accumulated knowledge from genetic and molecular biology experiments conducted in the pre-genomic era in order to attain the deepest level of understanding possible based on the available data. RESULTS: In this paper, we describe an expansion of RegulonDB, the database containing the rich legacy of decades of classic molecular biology experiments supporting what we know about gene regulation and operon organization in E. coli K-12, to include the genome-wide dataset collections from 32 ChIP and 19 gSELEX publications, in addition to around 60 genome-wide expression profiles relevant to the functional significance of these datasets and used in their curation. Three essential features for the integration of this information coming from different methodological approaches are: first, a controlled vocabulary within an ontology for precisely defining growth conditions; second, the criteria to separate elements with enough evidence to consider them involved in gene regulation from isolated transcription factor binding sites without such support; and third, an expanded computational model supporting this knowledge. Altogether, this constitutes the basis for adequately gathering and enabling the comparisons and integration needed to manage and access such wealth of knowledge. CONCLUSIONS: This version 10.0 of RegulonDB is a first step toward what should become the unifying access point for current and future knowledge on gene regulation in E. coli K-12. Furthermore, this model platform and associated methodologies and criteria can be emulated for gathering knowledge on other microbial organisms.

Characterization of a cAMP responsive transcription factor, Cmr (Rv1675c), in TB complex mycobacteria reveals overlap with the DosR (DevR) dormancy regulon.

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

Mycobacterium tuberculosis Whole Genome Sequences From Southern India Suggest Novel Resistance Mechanisms and the Need for Region-Specific Diagnostics.

BACKGROUND.: India is home to 25% of all tuberculosis cases and the second highest number of multidrug resistant cases worldwide. However, little is known about the genetic diversity and resistance determinants of Indian Mycobacterium tuberculosis, particularly for the primary lineages found in India, lineages 1 and 3. METHODS.: We whole genome sequenced 223 randomly selected M. tuberculosis strains from 196 patients within the Tiruvallur and Madurai districts of Tamil Nadu in Southern India. Using comparative genomics, we examined genetic diversity, transmission patterns, and evolution of resistance. RESULTS.: Genomic analyses revealed (11) prevalence of strains from lineages 1 and 3, (11) recent transmission of strains among patients from the same treatment centers, (11) emergence of drug resistance within patients over time, (11) resistance gained in an order typical of strains from different lineages and geographies, (11) underperformance of known resistance-conferring mutations to explain phenotypic resistance in Indian strains relative to studies focused on other geographies, and (11) the possibility that resistance arose through mutations not previously implicated in resistance, or through infections with multiple strains that confound genotype-based prediction of resistance. CONCLUSIONS.: In addition to substantially expanding the genomic perspectives of lineages 1 and 3, sequencing and analysis of M. tuberculosis whole genomes from Southern India highlight challenges of infection control and rapid diagnosis of resistant tuberculosis using current technologies. Further studies are needed to fully explore the complement of diversity and resistance determinants within endemic M. tuberculosis populations.

Decoding ChIP-seq with a double-binding signal refines binding peaks to single-nucleotides and predicts cooperative interaction.

The comprehension of protein and DNA binding in vivo is essential to understand gene regulation. Chromatin immunoprecipitation followed by sequencing (ChIP-seq) provides a global map of the regulatory binding network. Most ChIP-seq analysis tools focus on identifying binding regions from coverage enrichment. However, less work has been performed to infer the physical and regulatory details inside the enriched regions. This research extends a previous blind-deconvolution approach to develop a post-peak-calling algorithm that improves binding site resolution and predicts cooperative interactions. At the core of our new method is a physically motivated model that characterizes the binding signal as an extreme value distribution. This model suggests a mathematical framework to study physical properties of DNA shearing from the ChIP-seq coverage. The model explains the ChIP-seq coverage with two signals: The first considers DNA fragments with only a single binding event, whereas the second considers fragments with two binding events (a double-binding signal). The model incorporates motif discovery and is able to detect multiple sites in an enriched region with single-nucleotide resolution, high sensitivity, and high specificity. Our method improves peak caller sensitivity, from less than 45% up to 94%, at a false positive rate < 11% for a set of 47 experimentally validated prokaryotic sites. It also improves resolution of highly enriched regions of large-scale eukaryotic data sets. The double-binding signal provides a novel application in ChIP-seq analysis: the identification of cooperative interaction. Predictions of known cooperative binding sites show a 0.85 area under an ROC curve.

Genetic Determinants of Drug Resistance in Mycobacterium tuberculosis and Their Diagnostic Value.

RATIONALE: The development of molecular diagnostics that detect both the presence of Mycobacterium tuberculosis in clinical samples and drug resistance-conferring mutations promises to revolutionize patient care and interrupt transmission by ensuring early diagnosis. However, these tools require the identification of genetic determinants of resistance to the full range of antituberculosis drugs. OBJECTIVES: To determine the optimal molecular approach needed, we sought to create a comprehensive catalog of resistance mutations and assess their sensitivity and specificity in diagnosing drug resistance. METHODS: We developed and validated molecular inversion probes for DNA capture and deep sequencing of 28 drug-resistance loci in M. tuberculosis. We used the probes for targeted sequencing of a geographically diverse set of 1,397 clinical M. tuberculosis isolates with known drug resistance phenotypes. We identified a minimal set of mutations to predict resistance to first- and second-line antituberculosis drugs and validated our predictions in an independent dataset. We constructed and piloted a web-based database that provides public access to the sequence data and prediction tool. MEASUREMENTS AND MAIN RESULTS: The predicted resistance to rifampicin and isoniazid exceeded 90% sensitivity and specificity but was lower for other drugs. The number of mutations needed to diagnose resistance is large, and for the 13 drugs studied it was 238 across 18 genetic loci. CONCLUSIONS: These data suggest that a comprehensive M. tuberculosis drug resistance diagnostic will need to allow for a high dimension of mutation detection. They also support the hypothesis that currently unknown genetic determinants, potentially discoverable by whole-genome sequencing, encode resistance to second-line tuberculosis drugs.

Role of intragenic binding of cAMP responsive protein (CRP) in regulation of the succinate dehydrogenase genes Rv0249c-Rv0247c in TB complex mycobacteria.

Bacterial pathogens adapt to changing environments within their hosts, and the signaling molecule adenosine 3', 5'-cyclic monophosphate (cAMP) facilitates this process. In this study, we characterized in vivo DNA binding and gene regulation by the cAMP-responsive protein CRP in M. bovis BCG as a model for tuberculosis (TB)-complex bacteria. Chromatin immunoprecipitation followed by deep-sequencing (ChIP-seq) showed that CRP associates with ∼900 DNA binding regions, most of which occur within genes. The most highly enriched binding region was upstream of a putative copper transporter gene (ctpB), and crp-deleted bacteria showed increased sensitivity to copper toxicity. Detailed mutational analysis of four CRP binding sites upstream of the virulence-associated Rv0249c-Rv0247c succinate dehydrogenase genes demonstrated that CRP directly regulates Rv0249c-Rv0247c expression from two promoters, one of which requires sequences intragenic to Rv0250c for maximum expression. The high percentage of intragenic CRP binding sites and our demonstration that these intragenic DNA sequences significantly contribute to biologically relevant gene expression greatly expand the genome space that must be considered for gene regulatory analyses in mycobacteria. These findings also have practical implications for an important bacterial pathogen in which identification of mutations that affect expression of drug target-related genes is widely used for rapid drug resistance screening.

Analysis of clock-regulated genes in Neurospora reveals widespread posttranscriptional control of metabolic potential.

Neurospora crassa has been for decades a principal model for filamentous fungal genetics and physiology as well as for understanding the mechanism of circadian clocks. Eukaryotic fungal and animal clocks comprise transcription-translation-based feedback loops that control rhythmic transcription of a substantial fraction of these transcriptomes, yielding the changes in protein abundance that mediate circadian regulation of physiology and metabolism: Understanding circadian control of gene expression is key to understanding eukaryotic, including fungal, physiology. Indeed, the isolation of clock-controlled genes (ccgs) was pioneered in Neurospora where circadian output begins with binding of the core circadian transcription factor WCC to a subset of ccg promoters, including those of many transcription factors. High temporal resolution (2-h) sampling over 48 h using RNA sequencing (RNA-Seq) identified circadianly expressed genes in Neurospora, revealing that from ∼10% to as much 40% of the transcriptome can be expressed under circadian control. Functional classifications of these genes revealed strong enrichment in pathways involving metabolism, protein synthesis, and stress responses; in broad terms, daytime metabolic potential favors catabolism, energy production, and precursor assembly, whereas night activities favor biosynthesis of cellular components and growth. Discriminative regular expression motif elicitation (DREME) identified key promoter motifs highly correlated with the temporal regulation of ccgs. Correlations between ccg abundance from RNA-Seq, the degree of ccg-promoter activation as reported by ccg-promoter-luciferase fusions, and binding of WCC as measured by ChIP-Seq, are not strong. Therefore, although circadian activation is critical to ccg rhythmicity, posttranscriptional regulation plays a major role in determining rhythmicity at the mRNA level.

A genome-wide map of diversity in Plasmodium falciparum.

Genetic variation allows the malaria parasite Plasmodium falciparum to overcome chemotherapeutic agents, vaccines and vector control strategies and remain a leading cause of global morbidity and mortality. Here we describe an initial survey of genetic variation across the P. falciparum genome. We performed extensive sequencing of 16 geographically diverse parasites and identified 46,937 SNPs, demonstrating rich diversity among P. falciparum parasites (pi = 1.16 x 10(-3)) and strong correlation with gene function. We identified multiple regions with signatures of selective sweeps in drug-resistant parasites, including a previously unidentified 160-kb region with extremely low polymorphism in pyrimethamine-resistant parasites. We further characterized 54 worldwide isolates by genotyping SNPs across 20 genomic regions. These data begin to define population structure among African, Asian and American groups and illustrate the degree of linkage disequilibrium, which extends over relatively short distances in African parasites but over longer distances in Asian parasites. We provide an initial map of genetic diversity in P. falciparum and demonstrate its potential utility in identifying genes subject to recent natural selection and in understanding the population genetics of this parasite.

The role of selection in shaping diversity of natural M. tuberculosis populations.

Mycobacterium tuberculosis (M.tb), the cause of tuberculosis (TB), is estimated to infect a new host every second. While analyses of genetic data from natural populations of M.tb have emphasized the role of genetic drift in shaping patterns of diversity, the influence of natural selection on this successful pathogen is less well understood. We investigated the effects of natural selection on patterns of diversity in 63 globally extant genomes of M.tb and related pathogenic mycobacteria. We found evidence of strong purifying selection, with an estimated genome-wide selection coefficient equal to -9.5 × 10(-4) (95% CI -1.1 × 10(-3) to -6.8 × 10(-4)); this is several orders of magnitude higher than recent estimates for eukaryotic and prokaryotic organisms. We also identified different patterns of variation across categories of gene function. Genes involved in transport and metabolism of inorganic ions exhibited very low levels of non-synonymous polymorphism, equivalent to categories under strong purifying selection (essential and translation-associated genes). The highest levels of non-synonymous variation were seen in a group of transporter genes, likely due to either diversifying selection or local selective sweeps. In addition to selection, we identified other important influences on M.tb genetic diversity, such as a 25-fold expansion of global M.tb populations coincident with explosive growth in human populations (estimated timing 1684 C.E., 95% CI 1620-1713 C.E.). These results emphasize the parallel demographic histories of this obligate pathogen and its human host, and suggest that the dominant effect of selection on M.tb is removal of novel variants, with exceptions in an interesting group of genes involved in transportation and defense. We speculate that the hostile environment within a host imposes strict demands on M.tb physiology, and thus a substantial fitness cost for most new mutations. In this respect, obligate bacterial pathogens may differ from other host-associated microbes such as symbionts.

Evolutionary roles of upstream open reading frames in mediating gene regulation in fungi.

Upstream open reading frames (uORFs) are frequently present in the 5'-leader regions of fungal mRNAs. They can affect translation by controlling the ability of ribosomes that scan from the mRNA 5' end to reach the downstream genic reading frame. The translation of uORFs can also affect mRNA stability. For several genes, including Saccharomyces cerevisiae GCN4, S. cerevisiae CPA1, and Neurospora crassa arg-2, regulation by uORFs controls expression in response to specific physiological signals. The roles of many uORFs that are identified by genome-level approaches, as have been initiated for Saccharomyces, Aspergillus, and Cryptococcus species, remain to be determined. Some uORFs may have regulatory roles, while others may exist to insulate the genic reading frame from the negative impacts of upstream translation start sites in the mRNA 5' leader.