MultiFun, a multifunctional classification scheme for Escherichia coli K-12 gene products

Understanding the Genome-Wide Transcription Response To Various cAMP Levels in Bacteria Using Phenomenological Models

Attempts to understand gene regulation by global transcription factors have largely been limited to expression studies under binary conditions of presence and absence of the transcription factor. Studies addressing genome-wide transcriptional responses to changing transcription factor concentration at high resolution are lacking. Here, we create a data set containing the entire Escherichia coli transcriptome in Luria-Bertani (LB) broth as it responds to 10 different cAMP concentrations spanning the biological range. We use the Hill's model to accurately summarize individual gene responses into three intuitively understandable parameters, Emax, n, and k, reflecting the sensitivity, nonlinearity, and midpoint of the dynamic range. Our data show that most cAMP-regulated genes have an n of >2, with their k values centered around the wild-type concentration of cAMP. Additionally, cAMP receptor protein (CRP) affinity to a promoter is correlated with Emax but not k, hinting that a high-affinity CRP promoter need not ensure transcriptional activation at lower cAMP concentrations and instead affects the magnitude of the response. Finally, genes belonging to different functional classes are tuned to have different k, n, and Emax values. We demonstrate that phenomenological models are a better alternative for studying gene expression trends than classical clustering methods, with the phenomenological constants providing greater insights into how genes are tuned in a regulatory network. IMPORTANCE Different genes may follow different trends in response to various transcription factor concentrations. In this study, we ask two questions: (i) what are the trends that different genes follow in response to changing transcription factor concentrations and (ii) what methods can be used to extract information from the gene trends so obtained. We demonstrate a method to analyze transcription factor concentration-dependent genome-wide expression data using phenomenological models. Conventional clustering methods and principal-component analysis (PCA) can be used to summarize trends in data but have limited interpretability. The use of phenomenological models greatly enhances the interpretability and thus utility of conventional clustering. Transformation of dose-response data into phenomenological constants opens up avenues to ask and answer many different kinds of question. We show that the phenomenological constants obtained from the model fits can be used to generate insights about network topology and allows integration of other experimental data such as chromatin immunoprecipitation sequencing (ChIP-seq) to understand the system in greater detail.

Gene regulation in Escherichia coli is commonly selected for both high plasticity and low noise

Bacteria often respond to dynamically changing environments by regulating gene expression. Despite this regulation being critically important for growth and survival, little is known about how selection shapes gene regulation in natural populations. To better understand the role natural selection plays in shaping bacterial gene regulation, here we compare differences in the regulatory behaviour of naturally segregating promoter variants from Escherichia coli (which have been subject to natural selection) to randomly mutated promoter variants (which have never been exposed to natural selection). We quantify gene expression phenotypes (expression level, plasticity and noise) for hundreds of promoter variants across multiple environments and show that segregating promoter variants are enriched for mutations with minimal effects on expression level. In many promoters, we infer that there is strong selection to maintain high levels of plasticity, and direct selection to decrease or increase cell-to-cell variability in expression. Taken together, these results expand our knowledge of how gene regulation is affected by natural selection and highlight the power of comparing naturally segregating polymorphisms to de novo random mutations to quantify the action of selection.

Missing Links Between Gene Function and Physiology in Genomics

Knowledge of biological organisms at the molecular level that has been gathered is now organized into databases, often within ontological frameworks. To enable computational comparisons of annotations across different genomes and organisms, controlled vocabularies have been essential, as is the case in the functional annotation classifications used for bacteria, such as MultiFun and the more widely used Gene Ontology. The function of individual gene products as well as the processes in which collections of them participate constitute a wealth of classes that describe the biological role of gene products in a large number of organisms in the three kingdoms of life. In this contribution, we highlight from a qualitative perspective some limitations of these frameworks and discuss challenges that need to be addressed to bridge the gap between annotation as currently captured by ontologies and databases and our understanding of the basic principles in the organization and functioning of organisms; we illustrate these challenges with some examples in bacteria. We hope that raising awareness of these issues will encourage users of Gene Ontology and similar ontologies to be careful about data interpretation and lead to improved data representation.

Construction of a complete set of Neisseria meningitidis mutants and its use for the phenotypic profiling of this human pathogen

The bacterium Neisseria meningitidis causes life-threatening meningitis and sepsis. Here, we construct a complete collection of defined mutants in protein-coding genes of this organism, identifying all genes that are essential under laboratory conditions. The collection, named NeMeSys 2.0, consists of individual mutants in 1584 non-essential genes. We identify 391 essential genes, which are associated with basic functions such as expression and preservation of genome information, cell membrane structure and function, and metabolism. We use this collection to shed light on the functions of diverse genes, including a gene encoding a member of a previously unrecognised class of histidinol-phosphatases; a set of 20 genes required for type IV pili function; and several conditionally essential genes encoding antitoxins and/or immunity proteins. We expect that NeMeSys 2.0 will facilitate the phenotypic profiling of a major human bacterial pathogen.

A novel knock out strategy to enhance recombinant protein expression in Escherichia coli

Background

The expression of recombinant proteins triggers a stress response which downregulates key metabolic pathway genes leading to a decline in cellular health and feedback inhibition of both growth and protein expression. Instead of individually upregulating these downregulated genes or improving transcription rates by better vector design, an innovative strategy would be to block this stress response thereby ensuring a sustained level of protein expression.

Results

We postulated that the genes which are commonly up-regulated post induction may play the role of signalling messengers in mounting the cellular stress response. We identified those genes which have no known downstream regulatees and created knock outs which were then tested for GFP expression. Many of these knock outs showed significantly higher expression levels which was also sustained for longer periods. The highest product yield (Y_p/x) was observed in a BW25113ΔcysJ knock out (Y_p/x 0.57) and BW25113ΔelaA (Y_p/x 0.49), whereas the Y_p/x of the control W3110 strain was 0.08 and BW25113 was 0.16. Double knock out combinations were then created from the ten best performing single knock outs leading to a further enhancement in expression levels. Out of 45 double knock outs created, BW25113ΔelaAΔyhbC (Y_p/x 0.7) and BW25113ΔcysJΔyhbC (Y_p/x 0.64) showed the highest increase in product yield compared to the single gene mutant strains. We confirmed the improved performance of these knock outs by testing and obtaining higher levels of recombinant asparaginase expression, a system better suited for analysing sustained expression since it gets exported to the extracellular medium.

Conclusion

Creating key knock outs to block the CSR and enhance expression is a radically different strategy that can be synergistically combined with traditional methods of improving protein yields thus helping in the design of superior host platforms for protein expression.

MicroScope: an integrated platform for the annotation and exploration of microbial gene functions through genomic, pangenomic and metabolic comparative analysis

Large-scale genome sequencing and the increasingly massive use of high-throughput approaches produce a vast amount of new information that completely transforms our understanding of thousands of microbial species. However, despite the development of powerful bioinformatics approaches, full interpretation of the content of these genomes remains a difficult task. Launched in 2005, the MicroScope platform (https://www.genoscope.cns.fr/agc/microscope) has been under continuous development and provides analysis for prokaryotic genome projects together with metabolic network reconstruction and post-genomic experiments allowing users to improve the understanding of gene functions. Here we present new improvements of the MicroScope user interface for genome selection, navigation and expert gene annotation. Automatic functional annotation procedures of the platform have also been updated and we added several new tools for the functional annotation of genes and genomic regions. We finally focus on new tools and pipeline developed to perform comparative analyses on hundreds of genomes based on pangenome graphs. To date, MicroScope contains data for >11 800 microbial genomes, part of which are manually curated and maintained by microbiologists (>4500 personal accounts in September 2019). The platform enables collaborative work in a rich comparative genomic context and improves community-based curation efforts.

MicroScope-an integrated resource for community expertise of gene functions and comparative analysis of microbial genomic and metabolic data

The overwhelming list of new bacterial genomes becoming available on a daily basis makes accurate genome annotation an essential step that ultimately determines the relevance of thousands of genomes stored in public databanks. The MicroScope platform (http://www.genoscope.cns.fr/agc/microscope) is an integrative resource that supports systematic and efficient revision of microbial genome annotation, data management and comparative analysis. Starting from the results of our syntactic, functional and relational annotation pipelines, MicroScope provides an integrated environment for the expert annotation and comparative analysis of prokaryotic genomes. It combines tools and graphical interfaces to analyze genomes and to perform the manual curation of gene function in a comparative genomics and metabolic context. In this article, we describe the free-of-charge MicroScope services for the annotation and analysis of microbial (meta)genomes, transcriptomic and re-sequencing data. Then, the functionalities of the platform are presented in a way providing practical guidance and help to the nonspecialists in bioinformatics. Newly integrated analysis tools (i.e. prediction of virulence and resistance genes in bacterial genomes) and original method recently developed (the pan-genome graph representation) are also described. Integrated environments such as MicroScope clearly contribute, through the user community, to help maintaining accurate resources.

Forecasting cell fate during antibiotic exposure using stochastic gene expression

Antibiotic killing does not occur at a single, precise time for all cells within a population. Variability in time to death can be caused by stochastic expression of genes, resulting in differences in endogenous stress-resistance levels between individual cells in a population. Here we investigate whether single-cell differences in gene expression prior to antibiotic exposure are related to cell survival times after antibiotic exposure for a range of genes of diverse function. We quantified the time to death of single cells under antibiotic exposure in combination with expression of reporters. For some reporters, including genes involved in stress response and cellular processes like metabolism, the time to cell death had a strong relationship with the initial expression level of the genes. Our results highlight the single-cell level non-uniformity of antibiotic killing and also provide examples of key genes where cell-to-cell variation in expression is strongly linked to extended durations of antibiotic survival.

Exploring the transcriptome of luxI- and ΔainS mutants and the impact of N-3-oxo-hexanoyl-L- and N-3-hydroxy-decanoyl-L-homoserine lactones on biofilm formation in Aliivibrio salmonicida

Background

Bacterial communication through quorum sensing (QS) systems has been reported to be important in coordinating several traits such as biofilm formation. In Aliivibrio salmonicida two QS systems the LuxI/R and AinS/R, have been shown to be responsible for the production of eight acyl-homoserine lactones (AHLs) in a cell density dependent manner. We have previously demonstrated that inactivation of LitR, the master regulator of the QS system resulted in biofilm formation, similar to the biofilm formed by the AHL deficient mutant ΔainSluxI⁻. In this study, we aimed to investigate the global gene expression patterns of luxI and ainS autoinducer synthases mutants using transcriptomic profiling. In addition, we examined the influence of the different AHLs on biofilm formation.

Results

The transcriptome profiling of ΔainS and luxI⁻ mutants allowed us to identify genes and gene clusters regulated by QS in A. salmonicida. Relative to the wild type, the ΔainS and luxI⁻ mutants revealed 29 and 500 differentially expressed genes (DEGs), respectively. The functional analysis demonstrated that the most pronounced DEGs were involved in bacterial motility and chemotaxis, exopolysaccharide production, and surface structures related to adhesion. Inactivation of luxI, but not ainS genes resulted in wrinkled colony morphology. While inactivation of both genes (ΔainSluxI⁻) resulted in strains able to form wrinkled colonies and mushroom structured biofilm. Moreover, when the ΔainSluxI⁻ mutant was supplemented with N-3-oxo-hexanoyl-L-homoserine lactone (3OC6-HSL) or N-3-hydroxy-decanoyl-L-homoserine lactone (3OHC10-HSL), the biofilm did not develop. We also show that LuxI is needed for motility and for repression of EPS production, where repression of EPS is likely operated through the RpoQ-sigma factor.

Conclusion

These findings imply that the LuxI and AinS autoinducer synthases play a critical role in the regulation of biofilm formation, EPS production, and motility.

Tracing the phylogenetic history of the Crl regulon through the Bacteria and Archaea genomes

Background

Crl, identified for curli production, is a small transcription factor that stimulates the association of the σ^S factor (RpoS) with the RNA polymerase core through direct and specific interactions, increasing the transcription rate of genes during the transition from exponential to stationary phase at low temperatures, using indole as an effector molecule. The lack of a comprehensive collection of information on the Crl regulon makes it difficult to identify a dominant function of Crl and to generate any hypotheses concerning its taxonomical distribution in archaeal and bacterial organisms.

Results

In this work, based on a systematic literature review, we identified the first comprehensive dataset of 86 genes under the control of Crl in the bacterium Escherichia coli K-12; those genes correspond to 40% of the σ^S regulon in this bacterium. Based on an analysis of orthologs in 18 archaeal and 69 bacterial taxonomical divisions and using E. coli K-12 as a framework, we suggest three main events that resulted in this regulon’s actual form: (i) in a first step, rpoS, a gene widely distributed in bacteria and archaea cellular domains, was recruited to regulate genes involved in ancient metabolic processes, such as those associated with glycolysis and the tricarboxylic acid cycle; (ii) in a second step, the regulon recruited those genes involved in metabolic processes, which are mainly taxonomically constrained to Proteobacteria, with some secondary losses, such as those genes involved in responses to stress or starvation and cell adhesion, among others; and (iii) in a posterior step, Crl might have been recruited in Enterobacteriaceae; because its taxonomical pattern constrained to this bacterial order, however further analysis are necessary.

Conclusions

Therefore, we suggest that the regulon Crl is highly flexible for phenotypic adaptation, probably as consequence of the diverse growth environments associated with all organisms in which members of this regulatory network are present.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5619-z) contains supplementary material, which is available to authorized users.

Stochastic modelling reveals mechanisms of metabolic heterogeneity

Phenotypic variation is a hallmark of cellular physiology. Metabolic heterogeneity, in particular, underpins single-cell phenomena such as microbial drug tolerance and growth variability. Much research has focussed on transcriptomic and proteomic heterogeneity, yet it remains unclear if such variation permeates to the metabolic state of a cell. Here we propose a stochastic model to show that complex forms of metabolic heterogeneity emerge from fluctuations in enzyme expression and catalysis. The analysis predicts clonal populations to split into two or more metabolically distinct subpopulations. We reveal mechanisms not seen in deterministic models, in which enzymes with unimodal expression distributions lead to metabolites with a bimodal or multimodal distribution across the population. Based on published data, the results suggest that metabolite heterogeneity may be more pervasive than previously thought. Our work casts light on links between gene expression and metabolism, and provides a theory to probe the sources of metabolite heterogeneity.

Differential expression profiling of ΔlitR and ΔrpoQ mutants reveals insight into QS regulation of motility, adhesion and biofilm formation in Aliivibrio salmonicida

BACKGROUND

The coordination of group behaviors in bacteria is achieved by a cell-cell signaling process called quorum sensing (QS). QS is an intercellular communication system, which synchronously controls expression of a vast range of genes in response to changes in cell density and is mediated by autoinducers that act as extracellular signals. Aliivibrio salmonicida, the causative agent of cold-water vibrosis in marine aquacultures, uses QS to regulate several activities such as motility, biofilm formation, adhesion and rugose colony morphology. However, little is known about either genes or detailed mechanisms involved in the regulation of these phenotypes.

RESULTS

Differential expression profiling allowed us to define the genes involved in controlling phenotypes related to QS in A. salmonicida LFI1238. RNA sequencing data revealed that the number of expressed genes in A. salmonicida, ΔlitR and ΔrpoQ mutants were significantly altered due to changes in cell density. These included genes that were distributed among the 21 functional groups, mainly presented in cell envelope, cell processes, extrachromosomal/foreign DNA and transport-binding proteins functional groups. The comparative transcriptome of A. salmonicida wild-type at high cell density relative to low cell density revealed 1013 genes to be either up- or downregulated. Thirty-six downregulated genes were gene clusters encoding biosynthesis of the flagellar and chemotaxis genes. Additionally we identified significant expression for genes involved in acyl homoserine lactone (AHL) synthesis, adhesion and early colonization. The transcriptome profile of ΔrpoQ compared to the wild-type revealed 384 differensially expressed genes (DEGs) that allowed us to assign genes involved in regulating motility, adhesion and colony rugosity. Indicating the importance of RpoQ in controlling several QS related activities. Furthermore, the comparison of the transcriptome profiles of ΔlitR and ΔrpoQ mutants, exposed numerous overlapping DEGs that were essential for motility, exopolysaccharide production via syp operon and genes associated with tad operon.

CONCLUSION

Our findings indicate previously unexplained functional roles for LitR and RpoQ in regulation of different phenotypes related to QS. Our transcriptome data provide a better understanding of the regulation cascade of motility, wrinkling colony morphology and biofilm formation and will offer a major source for further research and analysis on this important field.

The EcoCyc Database

EcoCyc is a bioinformatics database available at EcoCyc.org that describes the genome and the biochemical machinery of Escherichia coli K-12 MG1655. The long-term goal of the project is to describe the complete molecular catalog of the E. coli cell, as well as the functions of each of its molecular parts, to facilitate a system-level understanding of E. coli. EcoCyc is an electronic reference source for E. coli biologists and for biologists who work with related microorganisms. The database includes information pages on each E. coli gene product, metabolite, reaction, operon, and metabolic pathway. The database also includes information on E. coli gene essentiality and on nutrient conditions that do or do not support the growth of E. coli. The website and downloadable software contain tools for analysis of high-throughput data sets. In addition, a steady-state metabolic flux model is generated from each new version of EcoCyc and can be executed via EcoCyc.org. The model can predict metabolic flux rates, nutrient uptake rates, and growth rates for different gene knockouts and nutrient conditions. This review outlines the data content of EcoCyc and of the procedures by which this content is generated.

Mining Novel Constitutive Promoter Elements in Soil Metagenomic Libraries in Escherichia coli

Although functional metagenomics has been widely employed for the discovery of genes relevant to biotechnology and biomedicine, its potential for assessing the diversity of transcriptional regulatory elements of microbial communities has remained poorly explored. Here, we experimentally mined novel constitutive promoter sequences in metagenomic libraries by combining a bi-directional reporter vector, high-throughput fluorescence assays and predictive computational methods. Through the expression profiling of fluorescent clones from two independent soil sample libraries, we have analyzed the regulatory dynamics of 260 clones with candidate promoters as a set of active metagenomic promoters in the host Escherichia coli. Through an in-depth analysis of selected clones, we were able to further explore the architecture of metagenomic fragments and to report the presence of multiple promoters per fragment with a dominant promoter driving the expression profile. These approaches resulted in the identification of 33 novel active promoters from metagenomic DNA originated from very diverse phylogenetic groups. The in silico and in vivo analysis of these individual promoters allowed the generation of a constitutive promoter consensus for exogenous sequences recognizable by E. coli in metagenomic studies. The results presented here demonstrates the potential of functional metagenomics for exploring environmental bacterial communities as a source of novel regulatory genetic parts to expand the toolbox for microbial engineering.

Construction of a fur null mutant and RNA-sequencing provide deeper global understanding of the Aliivibrio salmonicida Fur regulon

BACKGROUND

The ferric uptake regulator (Fur) is a transcription factor and the main regulator of iron acquisition in prokaryotes. When bound to ferric iron, Fur recognizes its DNA binding site and generally executes its function by repressing transcription of its target genes. Due to its importance in virulence, the Fur regulon is well studied for several model bacteria. In our previous work, we used computational predictions and microarray to gain insights into Fur-regulation in Aliivibrio salmonicida, and have identified a number of genes and operons that appear to be under direct control of Fur. To provide a more accurate and deeper global understanding of the biological role of Fur we have now generated an A. salmonicida fur knock-out strain and used RNA-sequencing to compare gene expression between the wild-type and fur null mutant strains.

RESULTS

An A. salmonicida fur null mutant strain was constructed. Biological assays demonstrate that deletion of fur results in loss of fitness, with reduced growth rates, and reduced abilities to withstand low-iron conditions, and oxidative stress. When comparing expression levels in the wild-type and the fur null mutant we retrieved 296 differentially expressed genes distributed among 18 of 21 functional classes of genes. A gene cluster encoding biosynthesis of the siderophore bisucaberin represented the highest up-regulated genes in the fur null mutant. Other highly up-regulated genes all encode proteins important for iron acquisition. Potential targets for the RyhB sRNA was predicted from the list of down-regulated genes, and significant complementarities were found between RyhB and mRNAs of the fur, sodB, cysN and VSAL_I0422 genes. Other sRNAs with potential functions in iron homeostasis were identified.

CONCLUSION

The present work provides by far the most comprehensive and deepest understanding of the Fur regulon in A. salmonicida to date. Our data also contribute to a better understanding of how Fur plays a key role in iron homeostasis in bacteria in general, and help to show how Fur orchestrates iron uptake when iron levels are extremely low.

Genomewide Mutational Diversity in Escherichia coli Population Evolving in Prolonged Stationary Phase

Prolonged stationary phase in bacteria, contrary to its name, is highly dynamic, with extreme nutrient limitation as a predominant stress. Stationary-phase cultures adapt by rapidly selecting a mutation(s) that confers a growth advantage in stationary phase (GASP). The phenotypic diversity of starving E. coli populations has been studied in detail; however, only a few mutations that accumulate in prolonged stationary phase have been described. This study documented the spectrum of mutations appearing in Escherichia coli during 28 days of prolonged starvation. The genetic diversity of the population increases over time in stationary phase to an extent that cannot be explained by random, neutral drift. This suggests that prolonged stationary phase offers a great model system to study adaptive evolution by natural selection.

Prolonged stationary phase is an approximation of natural environments presenting a range of stresses. Survival in prolonged stationary phase requires alternative metabolic pathways for survival. This study describes the repertoire of mutations accumulating in starving Escherichia coli populations in lysogeny broth. A wide range of mutations accumulates over the course of 1 month in stationary phase. Single nucleotide polymorphisms (SNPs) constitute 64% of all mutations. A majority of these mutations are nonsynonymous and are located at conserved loci. There is an increase in genetic diversity in the evolving populations over time. Computer simulations of evolution in stationary phase suggest that the maximum frequency of mutations observed in our experimental populations cannot be explained by neutral drift. Moreover, there is frequent genetic parallelism across populations, suggesting that these mutations are under positive selection. Finally, functional analysis of mutations suggests that regulatory mutations are frequent targets of selection.

IMPORTANCE Prolonged stationary phase in bacteria, contrary to its name, is highly dynamic, with extreme nutrient limitation as a predominant stress. Stationary-phase cultures adapt by rapidly selecting a mutation(s) that confers a growth advantage in stationary phase (GASP). The phenotypic diversity of starving E. coli populations has been studied in detail; however, only a few mutations that accumulate in prolonged stationary phase have been described. This study documented the spectrum of mutations appearing in Escherichia coli during 28 days of prolonged starvation. The genetic diversity of the population increases over time in stationary phase to an extent that cannot be explained by random, neutral drift. This suggests that prolonged stationary phase offers a great model system to study adaptive evolution by natural selection.

A genome-wide screen of bacterial mutants that enhance dauer formation in C. elegans

Molecular pathways involved in dauer formation, an alternate larval stage that allows Caenorhabditis elegans to survive adverse environmental conditions during development, also modulate longevity and metabolism. The decision to proceed with reproductive development or undergo diapause depends on food abundance, population density, and temperature. In recent years, the chemical identities of pheromone signals that modulate dauer entry have been characterized. However, signals derived from bacteria, the major source of nutrients for C. elegans, remain poorly characterized. To systematically identify bacterial components that influence dauer formation and aging in C. elegans, we utilized the individual gene deletion mutants in E. coli (K12). We identified 56 diverse E. coli deletion mutants that enhance dauer formation in an insulin-like receptor mutant (daf-2) background. We describe the mechanism of action of a bacterial mutant cyaA, that is defective in the production of cyclic AMP, which extends lifespan and enhances dauer formation through the modulation of TGF-β (daf-7) signaling in C. elegans. Our results demonstrate the importance of bacterial components in influencing developmental decisions and lifespan in C. elegans. Furthermore, we demonstrate that C. elegans is a useful model to study bacterial-host interactions.

Screen for genes involved in radiation survival of Escherichia coli and construction of a reference database

A set of 3907 single-gene knockout (Keio collection) strains of Escherichia coli K-12 was examined for strains with increased susceptibility to killing by X- or UV-radiation. After screening with a high-throughput resazurin-based assay and determining radiation survival with triplicate clonogenic assays, we identified 76 strains (and associated deleted genes) showing statistically-significant increased radiation sensitivity compared to a control strain. To determine gene novelty, we constructed a reference database comprised of genes found in nine similar studies including ours. This database contains 455 genes comprised of 103 common genes (found 2-7 times), and 352 uncommon genes (found once). Our 76 genes includes 43 common genes and 33 uncommon (potentially novel) genes, i.e., appY, atoS, betB, bglJ, clpP, cpxA, cysB, cysE, ddlA, dgkA, dppF, dusB, elfG, eutK, fadD, glnA, groL, guaB, intF, prpR, queA, rplY, seqA, sufC,yadG, yagJ, yahD, yahO, ybaK, ybfA, yfaL, yhjV, and yiaL. Of our 33 uncommon gene mutants, 4 (12%) were sensitive only to UV-radiation, 10 (30%) only to X-radiation, and 19 (58%) to both radiations. Our uncommon mutants vs. our common mutants showed more radiation specificity, i.e., 12% vs. 9% (sensitive only to UV-); 30% vs. 16% (X-) and 58% vs. 74% (both radiations). Considering just our radiation-sensitive mutants, the median UV-radiation survival (75Jm^-2) for 23 uncommon mutants was 6.84E-3 compared to 1.85E-3 for 36 common mutants (P=0.025). Similarly, the average X-radiation survival for 29 uncommon mutants was 1.08E-2, compared to 6.19E-3 for 39 common mutants (P=0.010). Comparing gene functions using MultiFun terms, uncommon genes tended to show less involvement in DNA repair-relevant categories (information transfer and cell processes), but greater involvement in seven other categories. Our analysis of 455 genes suggests cell survival and DNA repair processes are more complex than previously understood, and may be compromised by deficiencies in other processes.

Combining Shigella Tn-seq data with gold-standard E. coli gene deletion data suggests rare transitions between essential and non-essential gene functionality

BACKGROUND

Gene essentiality - whether or not a gene is necessary for cell growth - is a fundamental component of gene function. It is not well established how quickly gene essentiality can change, as few studies have compared empirical measures of essentiality between closely related organisms.

RESULTS

Here we present the results of a Tn-seq experiment designed to detect essential protein coding genes in the bacterial pathogen Shigella flexneri 2a 2457T on a genome-wide scale. Superficial analysis of this data suggested that 481 protein-coding genes in this Shigella strain are critical for robust cellular growth on rich media. Comparison of this set of genes with a gold-standard data set of essential genes in the closely related Escherichia coli K12 BW25113 revealed that an excessive number of genes appeared essential in Shigella but non-essential in E. coli. Importantly, and in converse to this comparison, we found no genes that were essential in E. coli and non-essential in Shigella, implying that many genes were artefactually inferred as essential in Shigella. Controlling for such artefacts resulted in a much smaller set of discrepant genes. Among these, we identified three sets of functionally related genes, two of which have previously been implicated as critical for Shigella growth, but which are dispensable for E. coli growth.

CONCLUSIONS

The data presented here highlight the small number of protein coding genes for which we have strong evidence that their essentiality status differs between the closely related bacterial taxa E. coli and Shigella. A set of genes involved in acetate utilization provides a canonical example. These results leave open the possibility of developing strain-specific antibiotic treatments targeting such differentially essential genes, but suggest that such opportunities may be rare in closely related bacteria.

Shifts in Buchnera aphidicola density in soybean aphids (Aphis glycines) feeding on virus-infected soybean

Vertically transmitted bacterial symbionts are common in arthropods. Aphids undergo an obligate symbiosis with Buchnera aphidicola, which provides essential amino acids to its host and contributes directly to nymph growth and reproduction. We previously found that newly adult Aphis glycines feeding on soybean infected with the beetle-transmitted Bean pod mottle virus (BPMV) had significantly reduced fecundity. We hypothesized that the reduced fecundity was attributable to detrimental impacts of the virus on the aphid microbiome, namely Buchnera. To test this, mRNA sequencing and quantitative real-time PCR were used to assay Buchnera transcript abundance and titre in A. glycines feeding on Soybean mosaic virus-infected, BPMV-infected, and healthy soybean for up to 14 days. Our results indicated that Buchnera density was lower and ultimately suppressed in aphids feeding on virus-infected soybean. While the decreased Buchnera titre may be associated with reduced aphid fecundity, additional mechanisms are probably involved. The present report begins to describe how interactions among insects, plants, and plant pathogens influence endosymbiont population dynamics.

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

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

Time-series analysis of the transcriptome and proteome of Escherichia coli upon glucose repression

Time-series transcript- and protein-profiles were measured upon initiation of carbon catabolite repression in Escherichia coli, in order to investigate the extent of post-transcriptional control in this prototypical response. A glucose-limited chemostat culture was used as the CCR-free reference condition. Stopping the pump and simultaneously adding a pulse of glucose, that saturated the cells for at least 1h, was used to initiate the glucose response. Samples were collected and subjected to quantitative time-series analysis of both the transcriptome (using microarray analysis) and the proteome (through a combination of 15N-metabolic labeling and mass spectrometry). Changes in the transcriptome and corresponding proteome were analyzed using statistical procedures designed specifically for time-series data. By comparison of the two sets of data, a total of 96 genes were identified that are post-transcriptionally regulated. This gene list provides candidates for future in-depth investigation of the molecular mechanisms involved in post-transcriptional regulation during carbon catabolite repression in E. coli, like the involvement of small RNAs.

Global coordination in adaptation to gene rewiring

Gene rewiring is a common evolutionary phenomenon in nature that may lead to extinction for living organisms. Recent studies on synthetic biology demonstrate that cells can survive genetic rewiring. This survival (adaptation) is often linked to the stochastic expression of rewired genes with random transcriptional changes. However, the probability of adaptation and the underlying common principles are not clear. We performed a systematic survey of an assortment of gene-rewired Escherichia coli strains to address these questions. Three different cell fates, designated good survivors, poor survivors and failures, were observed when the strains starved. Large fluctuations in the expression of the rewired gene were commonly observed with increasing cell size, but these changes were insufficient for adaptation. Cooperative reorganizations in the corresponding operon and genome-wide gene expression largely contributed to the final success. Transcriptome reorganizations that generally showed high-dimensional dynamic changes were restricted within a one-dimensional trajectory for adaptation to gene rewiring, indicating a general path directed toward cellular plasticity for a successful cell fate. This finding of global coordination supports a mechanism of stochastic adaptation and provides novel insights into the design and application of complex genetic or metabolic networks.

Proteome-wide subcellular topologies of E. coli polypeptides database (STEPdb)

Cell compartmentalization serves both the isolation and the specialization of cell functions. After synthesis in the cytoplasm, over a third of all proteins are targeted to other subcellular compartments. Knowing how proteins are distributed within the cell and how they interact is a prerequisite for understanding it as a whole. Surface and secreted proteins are important pathogenicity determinants. Here we present the STEP database (STEPdb) that contains a comprehensive characterization of subcellular localization and topology of the complete proteome of Escherichia coli. Two widely used E. coli proteomes (K-12 and BL21) are presented organized into thirteen subcellular classes. STEPdb exploits the wealth of genetic, proteomic, biochemical, and functional information on protein localization, secretion, and targeting in E. coli, one of the best understood model organisms. Subcellular annotations were derived from a combination of bioinformatics prediction, proteomic, biochemical, functional, topological data and extensive literature re-examination that were refined through manual curation. Strong experimental support for the location of 1553 out of 4303 proteins was based on 426 articles and some experimental indications for another 526. Annotations were provided for another 320 proteins based on firm bioinformatic predictions. STEPdb is the first database that contains an extensive set of peripheral IM proteins (PIM proteins) and includes their graphical visualization into complexes, cellular functions, and interactions. It also summarizes all currently known protein export machineries of E. coli K-12 and pairs them, where available, with the secretory proteins that use them. It catalogs the Sec- and TAT-utilizing secretomes and summarizes their topological features such as signal peptides and transmembrane regions, transmembrane topologies and orientations. It also catalogs physicochemical and structural features that influence topology such as abundance, solubility, disorder, heat resistance, and structural domain families. Finally, STEPdb incorporates prediction tools for topology (TMHMM, SignalP, and Phobius) and disorder (IUPred) and implements the BLAST2STEP that performs protein homology searches against the STEPdb.

The EcoCyc Database

EcoCyc is a bioinformatics database available at EcoCyc.org that describes the genome and the biochemical machinery of Escherichia coli K-12 MG1655. The long-term goal of the project is to describe the complete molecular catalog of the E. coli cell, as well as the functions of each of its molecular parts, to facilitate a system-level understanding of E. coli. EcoCyc is an electronic reference source for E. coli biologists and for biologists who work with related microorganisms. The database includes information pages on each E. coli gene, metabolite, reaction, operon, and metabolic pathway. The database also includes information on E. coli gene essentiality and on nutrient conditions that do or do not support the growth of E. coli. The website and downloadable software contain tools for analysis of high-throughput data sets. In addition, a steady-state metabolic flux model is generated from each new version of EcoCyc. The model can predict metabolic flux rates, nutrient uptake rates, and growth rates for different gene knockouts and nutrient conditions. This review provides a detailed description of the data content of EcoCyc and of the procedures by which this content is generated.

Data visualization in environmental proteomics

From raw data to gene expression profiles, from single cultures to complex microbial communities, environmental proteomics works with data of different complexity levels that need to be interpreted in detail or in its entirety. Although data visualization is closely connected with data analysis approaches, this work will solely focus on data visualization. Complementing traditional tools such as bar charts or line graphs, scientists and visualization professionals have been provided sophisticated visualization tools. Many rules and concerns regarding the display of single but also complex data will be reviewed and discussed. Visual approaches such as microcharts, heat maps, stream graphs, and tree maps will be brought to the reader's attention and demonstrated by utilizing real data sets.

Annotation enrichment analysis: an alternative method for evaluating the functional properties of gene sets

Gene annotation databases (compendiums maintained by the scientific community that describe the biological functions performed by individual genes) are commonly used to evaluate the functional properties of experimentally derived gene sets. Overlap statistics, such as Fishers Exact test (FET), are often employed to assess these associations, but don't account for non-uniformity in the number of genes annotated to individual functions or the number of functions associated with individual genes. We find FET is strongly biased toward over-estimating overlap significance if a gene set has an unusually high number of annotations. To correct for these biases, we develop Annotation Enrichment Analysis (AEA), which properly accounts for the non-uniformity of annotations. We show that AEA is able to identify biologically meaningful functional enrichments that are obscured by numerous false-positive enrichment scores in FET, and we therefore suggest it be used to more accurately assess the biological properties of gene sets.

The predictability of molecular evolution during functional innovation

Determining the molecular changes that give rise to functional innovations is a major unresolved problem in biology. The paucity of examples has served as a significant hindrance in furthering our understanding of this process. Here we used experimental evolution with the bacterium Escherichia coli to quantify the molecular changes underlying functional innovation in 68 independent instances ranging over 22 different metabolic functions. Using whole-genome sequencing, we show that the relative contribution of regulatory and structural mutations depends on the cellular context of the metabolic function. In addition, we find that regulatory mutations affect genes that act in pathways relevant to the novel function, whereas structural mutations affect genes that act in unrelated pathways. Finally, we use population genetic modeling to show that the relative contributions of regulatory and structural mutations during functional innovation may be affected by population size. These results provide a predictive framework for the molecular basis of evolutionary innovation, which is essential for anticipating future evolutionary trajectories in the face of rapid environmental change.

Growth rate-coordinated transcriptome reorganization in bacteria

Background

Cell growth rate reflects an organism’s physiological state and largely relies on the ability of gene expression to respond to the environment. The relationship between cellular growth rate and gene expression remains unknown.

Results

Growth rate-coordinated changes in gene expression were discovered by analyzing exponentially growing Escherichia coli cells cultured under multiple defined environments, in which osmotic pressure, temperature and starvation status were varied. Gene expression analyses showed that all 3,740 genes in the genome could be simply divided into three clusters (C1, C2 and C3), which were accompanied by a generic trend in the growth rate that was coordinated with transcriptional changes. The direction of transcriptional change in C1 indicated environmental specificity, whereas those in C2 and C3 were correlated negatively and positively with growth rates, respectively. The three clusters exhibited differentiated gene functions and gene regulation task division.

Conclusions

We identified three gene clusters, exhibiting differential gene functions and distinct directions in their correlations with growth rates. Reverses in the direction of the growth rate correlated transcriptional changes and the distinguished duties of the three clusters indicated how transcriptome homeostasis is maintained to balance the total expression cost for sustaining life in new habitats.

Genome sequence of Xanthomonas fuscans subsp. fuscans strain 4834-R reveals that flagellar motility is not a general feature of xanthomonads

BACKGROUND

Xanthomonads are plant-associated bacteria responsible for diseases on economically important crops. Xanthomonas fuscans subsp. fuscans (Xff) is one of the causal agents of common bacterial blight of bean. In this study, the complete genome sequence of strain Xff 4834-R was determined and compared to other Xanthomonas genome sequences.

RESULTS

Comparative genomics analyses revealed core characteristics shared between Xff 4834-R and other xanthomonads including chemotaxis elements, two-component systems, TonB-dependent transporters, secretion systems (from T1SS to T6SS) and multiple effectors. For instance a repertoire of 29 Type 3 Effectors (T3Es) with two Transcription Activator-Like Effectors was predicted. Mobile elements were associated with major modifications in the genome structure and gene content in comparison to other Xanthomonas genomes. Notably, a deletion of 33 kbp affects flagellum biosynthesis in Xff 4834-R. The presence of a complete flagellar cluster was assessed in a collection of more than 300 strains representing different species and pathovars of Xanthomonas. Five percent of the tested strains presented a deletion in the flagellar cluster and were non-motile. Moreover, half of the Xff strains isolated from the same epidemic than 4834-R was non-motile and this ratio was conserved in the strains colonizing the next bean seed generations.

CONCLUSIONS

This work describes the first genome of a Xanthomonas strain pathogenic on bean and reports the existence of non-motile xanthomonads belonging to different species and pathovars. Isolation of such Xff variants from a natural epidemic may suggest that flagellar motility is not a key function for in planta fitness.

Transcription factors in Escherichia coli prefer the holo conformation

The transcriptional regulatory network of Escherichia coli K-12 is among the best studied gene networks of any living cell. Transcription factors bind to DNA either with their effector bound (holo conformation), or as a free protein (apo conformation) regulating transcription initiation. By using RegulonDB, the functional conformations (holo or apo) of transcription factors, and their mode of regulation (activator, repressor, or dual) were exhaustively analyzed. We report a striking discovery in the architecture of the regulatory network, finding a strong under-representation of the apo conformation (without allosteric metabolite) of transcription factors when binding to their DNA sites to activate transcription. This observation is supported at the level of individual regulatory interactions on promoters, even if we exclude the promoters regulated by global transcription factors, where three-quarters of the known promoters are regulated by a transcription factor in holo conformation. This genome-scale analysis enables us to ask what are the implications of these observations for the physiology and for our understanding of the ecology of E. coli. We discuss these ideas within the framework of the demand theory of gene regulation.

Proteomic analysis of an unculturable bacterial endosymbiont (Blochmannia) reveals high abundance of chaperonins and biosynthetic enzymes

Many insect groups have coevolved with bacterial endosymbionts that live within specialized host cells. As a salient example, ants in the tribe Camponotini rely on Blochmannia, an intracellular bacterial mutualist that synthesizes amino acids and recycles nitrogen for the host. We performed a shotgun, label-free, LC/MS/MS quantitative proteomic analysis to investigate the proteome of Blochmannia associated with Camponotus chromaiodes. We identified more than 330 Blochmannia proteins, or 54% coverage of the predicted proteome, as well as 244 Camponotus proteins. Using the average intensity of the top 3 "best flier" peptides along with spiking of a surrogate standard at a known concentration, we estimated the concentration (fmol/μg) of those proteins with confident identification. The estimated dynamic range of Blochmannia protein abundance spanned 3 orders of magnitude and covered diverse functional categories, with particularly high representation of metabolism, information transfer, and chaperones. GroEL, the most abundant protein, totaled 6% of Blochmannia protein abundance. Biosynthesis of essential amino acids, fatty acids, and nucleotides, and sulfate assimilation had disproportionately high coverage in the proteome, further supporting a nutritional role of the symbiosis. This first quantitative proteomic analysis of an ant endosymbiont illustrates a promising approach to study the functional basis of intimate symbioses.

Patterns of evolutionary conservation of essential genes correlate with their compensability

Essential genes code for fundamental cellular functions required for the viability of an organism. For this reason, essential genes are often highly conserved across organisms. However, this is not always the case: orthologues of genes that are essential in one organism are sometimes not essential in other organisms or are absent from their genomes. This suggests that, in the course of evolution, essential genes can be rendered nonessential. How can a gene become non-essential? Here we used genetic manipulation to deplete the products of 26 different essential genes in Escherichia coli. This depletion results in a lethal phenotype, which could often be rescued by the overexpression of a non-homologous, non-essential gene, most likely through replacement of the essential function. We also show that, in a smaller number of cases, the essential genes can be fully deleted from the genome, suggesting that complete functional replacement is possible. Finally, we show that essential genes whose function can be replaced in the laboratory are more likely to be non-essential or not present in other taxa. These results are consistent with the notion that patterns of evolutionary conservation of essential genes are influenced by their compensability-that is, by how easily they can be functionally replaced, for example through increased expression of other genes.

Membrane transport metabolons

In this review evidence from a wide variety of biological systems is presented for the genetic, functional, and likely physical association of membrane transporters and the enzymes that metabolize the transported substrates. This evidence supports the hypothesis that the dynamic association of transporters and enzymes creates functional membrane transport metabolons that channel substrates typically obtained from the extracellular compartment directly into their cellular metabolism. The immediate modification of substrates on the inner surface of the membrane prevents back-flux through facilitated transporters, increasing the efficiency of transport. In some cases products of the enzymes are themselves substrates for the transporters that efflux the products in an exchange or antiport mechanism. Regulation of the binding of enzymes to transporters and their mutual activities may play a role in modulating flux through transporters and entry of substrates into metabolic pathways. Examples showing the physical association of transporters and enzymes are provided, but available structural data is sparse. Genetic and functional linkages between membrane transporters and enzymes were revealed by an analysis of Escherichia coli operons encoding polycistronic mRNAs and provide a list of predicted interactions ripe for further structural studies. This article supports the view that membrane transport metabolons are important throughout Nature in organisms ranging from bacteria to humans.

A genome-wide analysis of promoter-mediated phenotypic noise in Escherichia coli

Gene expression is subject to random perturbations that lead to fluctuations in the rate of protein production. As a consequence, for any given protein, genetically identical organisms living in a constant environment will contain different amounts of that particular protein, resulting in different phenotypes. This phenomenon is known as “phenotypic noise.” In bacterial systems, previous studies have shown that, for specific genes, both transcriptional and translational processes affect phenotypic noise. Here, we focus on how the promoter regions of genes affect noise and ask whether levels of promoter-mediated noise are correlated with genes' functional attributes, using data for over 60% of all promoters in Escherichia coli. We find that essential genes and genes with a high degree of evolutionary conservation have promoters that confer low levels of noise. We also find that the level of noise cannot be attributed to the evolutionary time that different genes have spent in the genome of E. coli. In contrast to previous results in eukaryotes, we find no association between promoter-mediated noise and gene expression plasticity. These results are consistent with the hypothesis that, in bacteria, natural selection can act to reduce gene expression noise and that some of this noise is controlled through the sequence of the promoter region alone.

Many biological processes in a cell involve small numbers of molecules and therefore fluctuate over time. As a consequence, genetically identical cells that live in the same environment differ from each other in many phenotypic traits, including the expression level of different genes. Our aim was to identify types of genes with particularly low or high levels of variation (“noise”) and to understand molecular and evolutionary factors that determine noise level. Working with the bacterium E. coli, we analyzed the expression—at the single cell level—of more than 1,500 different genes. We found particularly low levels of noise in genes that E. coli needs to live and genes that this bacterium shares with many related taxa. This suggests that cellular functions that are particularly important for this organism evolved towards low levels of variation. In contrast to previous results with yeast, we find that genes that change their expression levels in response to environmental signals do not have high levels of noise. This suggests that there may be fundamental differences in how noise is controlled in bacteria and eukaryotes.

Stringent response of Escherichia coli: revisiting the bibliome using literature mining

Background

Understanding the mechanisms responsible for cellular responses depends on the systematic collection and analysis of information on the main biological concepts involved. Indeed, the identification of biologically relevant concepts in free text, namely genes, tRNAs, mRNAs, gene products and small molecules, is crucial to capture the structure and functioning of different responses.

Results

In this work, we review literature reports on the study of the stringent response in Escherichia coli. Rather than undertaking the development of a highly specialised literature mining approach, we investigate the suitability of concept recognition and statistical analysis of concept occurrence as means to highlight the concepts that are most likely to be biologically engaged during this response. The co-occurrence analysis of core concepts in this stringent response, i.e. the (p)ppGpp nucleotides with gene products was also inspected and suggest that besides the enzymes RelA and SpoT that control the basal levels of (p)ppGpp nucleotides, many other proteins have a key role in this response. Functional enrichment analysis revealed that basic cellular processes such as metabolism, transcriptional and translational regulation are central, but other stress-associated responses might be elicited during the stringent response. In addition, the identification of less annotated concepts revealed that some (p)ppGpp-induced functional activities are still overlooked in most reviews.

Conclusions

In this paper we applied a literature mining approach that offers a more comprehensive analysis of the stringent response in E. coli. The compilation of relevant biological entities to this stress response and the assessment of their functional roles provided a more systematic understanding of this cellular response. Overlooked regulatory entities, such as transcriptional regulators, were found to play a role in this stress response. Moreover, the involvement of other stress-associated concepts demonstrates the complexity of this cellular response.

Reduced selective constraint in endosymbionts: elevation in radical amino acid replacements occurs genome-wide

As predicted by the nearly neutral model of evolution, numerous studies have shown that reduced N(e) accelerates the accumulation of slightly deleterious changes under genetic drift. While such studies have mostly focused on eukaryotes, bacteria also offer excellent models to explore the effects of N(e). Most notably, the genomes of host-dependent bacteria with small N(e) show signatures of genetic drift, including elevated K(a)/K(s). Here, I explore the utility of an alternative measure of selective constraint: the per-site rate of radical and conservative amino acid substitutions (D(r)/D(c)). I test the hypothesis that purifying selection against radical amino acid changes is less effective in two insect endosymbiont groups (Blochmannia of ants and Buchnera of aphids), compared to related gamma-Proteobacteria. Genome comparisons demonstrate a significant elevation in D(r)/D(c) in endosymbionts that affects the majority (66-79%) of shared orthologs examined. The elevation of D(r)/D(c) in endosymbionts affects all functional categories examined. Simulations indicate that D(r)/D(c) estimates are sensitive to codon frequencies and mutational parameters; however, estimation biases occur in the opposite direction as the patterns observed in genome comparisons, thereby making the inference of elevated D(r)/D(c) more conservative. Increased D(r)/D(c) and other signatures of genome degradation in endosymbionts are consistent with strong effects of genetic drift in their small populations, as well as linkage to selected sites in these asexual bacteria. While relaxed selection against radical substitutions may contribute, genome-wide processes such as genetic drift and linkage best explain the pervasive elevation in D(r)/D(c) across diverse functional categories that include basic cellular processes. Although the current study focuses on a few bacterial lineages, it suggests D(r)/D(c) is a useful gauge of selective constraint and may provide a valuable alternative to K(a)/K(s) when high sequence divergences preclude estimates of K(s). Broader application of D(r)/D(c) will benefit from approaches less prone to estimation biases.

The association of DNA damage response and nucleotide level modulation with the antibacterial mechanism of the anti-folate drug trimethoprim

BACKGROUND

Trimethoprim is a widely prescribed antibiotic for a variety of bacterial infections. It belongs to a class of anti-metabolites - antifolates - which includes drugs used against malarial parasites and in cancer therapy. However, spread of bacterial resistance to the drug has severely hampered its clinical use and has necessitated further investigations into its mechanism of action and treatment regimen. Trimethoprim selectively starves bacterial cells for tetrahydrofolate, a vital cofactor necessary for the synthesis of several metabolites. The outcome (bacteriostatic or bactericidal) of such starvation, however, depends on the availability of folate-dependent metabolites in the growth medium. To characterize this dependency, we investigated in detail the regulatory and structural components of Escherichia coli cellular response to trimethoprim in controlled growth and supplementation conditions.

RESULTS

We surveyed transcriptional responses to trimethoprim treatment during bacteriostatic and bactericidal conditions and analyzed associated gene sets/pathways. Concurrent starvation of all folate dependent metabolites caused growth arrest, and this was accompanied by induction of general stress and stringent responses. Three gene sets were significantly associated with the bactericidal effect of TMP in different media including LB: genes of the SOS regulon, genes of the pyrimidine nucleotide biosynthetic pathway and members of the multiple antibiotic resistance (mar) regulon controlled by the MarR repressor. However, the SOS response was identified as the only universal transcriptional signature associated with the loss of viability by direct thymine starvation or by folate stress. We also used genome-wide gene knock-out screen to uncover means of sensitization of bacteria to the drug. We observed that among a number of candidate genes and pathways, the effect of knock-outs in the deoxyribose nucleotide salvage pathway, encoded by the deoCABD operon and under the control of the DeoR repressor, was most informative.

CONCLUSION

Transcriptional induction of DNA damage response is an essential feature of the bactericidal effect of trimethoprim. Either the observation of the transcriptional response or DNA damage itself, or both, is made possible by thymine starvation when other folate-dependent metabolites are not limited. The effect of DNA damage by the drug takes place prior to its bactericidal effect, at the beginning of the lag stage of the treatment. Mutations in the deoxyribose nucleotide salvage pathway can affect duration of the lag as well as the rate of killing. This information can be used to postulate certain mechanistic differences between direct thymine starvation in thymidylate synthase deficient mutants and thymine starvation by anti-folate inhibitors.

SubtiWiki--a comprehensive community resource for the model organism Bacillus subtilis

In the post-genomic era, most components of a cell are known and they can be quantified by large-scale functional genomics approaches. However, genome annotation is the bottleneck that hampers our understanding of living cells and organisms. Up-to-date functional annotation is of special importance for model organisms that provide a frame of reference for studies with other relevant organisms. We have generated a Wiki-type database for the Gram-positive model bacterium Bacillus subtilis, SubtiWiki (http://subtiwiki.uni-goettingen.de/). This Wiki is centered around the individual genes and gene products of B. subtilis and provides information on each aspect of gene function and expression as well as protein activity and its control. SubtiWiki is accompanied by two companion databases SubtiPathways and SubtInteract that provide graphical representations of B. subtilis metabolism and its regulation and of protein-protein interactions, respectively. The diagrams of both databases are easily navigatable using the popular Google maps API, and they are extensively linked with the SubtiWiki gene pages. Moreover, each gene/gene product was assigned to one or more functional categories and transcription factor regulons. Pages for the specific categories and regulons provide a rapid overview of functionally related genes/proteins. Today, SubtiWiki can be regarded as one of the most complete inventories of knowledge on a living organism in one single resource.

Transcriptional profiling of the iron starvation response in Bordetella pertussis provides new insights into siderophore utilization and virulence gene expression

Serological studies of patients with pertussis and the identification of antigenic Bordetella pertussis proteins support the hypothesis that B. pertussis perceives an iron starvation cue and expresses multiple iron source utilization systems in its natural human host environment. Furthermore, previous studies using a murine respiratory tract infection model showed that several of these B. pertussis iron systems are required for colonization and persistence and are differentially expressed over the course of infection. The present study examined genome-wide changes in B. pertussis gene transcript abundance in response to iron starvation in vitro. In addition to known iron source utilization genes, we identified a previously uncharacterized iron-repressed cytoplasmic membrane transporter system, fbpABC, that is required for the utilization of multiple structurally distinct siderophores including alcaligin, enterobactin, ferrichrome, and desferrioxamine B. Expression of type III secretion system genes was also found to be upregulated during iron starvation in both B. pertussis strain Tohama I and Bordetella bronchiseptica strain RB50. In a survey of type III secretion system protein production by an assortment of B. pertussis laboratory-adapted and low-passage clinical isolate strains, iron limitation increased the production and secretion of the type III secretion system-specific translocation apparatus tip protein Bsp22 in all Bvg-proficient strains. These results indicate that iron starvation in the infected host is an important environmental cue influencing not only Bordetella iron transport gene expression but also the expression of other important virulence-associated genes.

Large-scale label-free quantitative proteomics of the pea aphid-Buchnera symbiosis

Many insects are nutritionally dependent on symbiotic microorganisms that have tiny genomes and are housed in specialized host cells called bacteriocytes. The obligate symbiosis between the pea aphid Acyrthosiphon pisum and the γ-proteobacterium Buchnera aphidicola (only 584 predicted proteins) is particularly amenable for molecular analysis because the genomes of both partners have been sequenced. To better define the symbiotic relationship between this aphid and Buchnera, we used large-scale, high accuracy tandem mass spectrometry (nanoLC-LTQ-Orbtrap) to identify aphid and Buchnera proteins in the whole aphid body, purified bacteriocytes, isolated Buchnera cells and the residual bacteriocyte fraction. More than 1900 aphid and 400 Buchnera proteins were identified. All enzymes in amino acid metabolism annotated in the Buchnera genome were detected, reflecting the high (68%) coverage of the proteome and supporting the core function of Buchnera in the aphid symbiosis. Transporters mediating the transport of predicted metabolites were present in the bacteriocyte. Label-free spectral counting combined with hierarchical clustering, allowed to define the quantitative distribution of a subset of these proteins across both symbiotic partners, yielding no evidence for the selective transfer of protein among the partners in either direction. This is the first quantitative proteome analysis of bacteriocyte symbiosis, providing a wealth of information about molecular function of both the host cell and bacterial symbiont.

Sequence information on simple sequence repeats and single nucleotide polymorphisms through transcriptome analysis of mungbean

Mungbean (Vigna radiata (L.) Wilczek) is a unique species in its ability to fix atmospheric nitrogen, with early maturity, and relatively good drought resistance. We used 454 sequencing technology for transcriptome sequencing. A total of 150 159 and 142 993 reads produced 5 254 and 6 374 large contigs (≥ 500 bp) with an average length of 833 and 853 for Sunhwa and Jangan, respectively. Functional annotation to known sequences yielded 41.34% and 41.74% unigenes for Jangan and Sunhwa. A higher number of simple sequence repeat (SSR) motifs was identified in Jangan (1 630) compared with that of Sunhwa (1 334). A similar SSR distribution pattern was observed in both varieties. A total of 8 249 single nucleotide polymorphisms (SNPs) and indels with 2 098 high-confidence candidates were identified in the two mungbean varieties. The average distance between individual SNPs was approximately 860 bp. Our report demonstrates the utility of transcriptomic data for implementing a functional annotation and development of genetic markers. We also provide large resource sequence data for mungbean improvement programs.

Array-based synthetic genetic screens to map bacterial pathways and functional networks in Escherichia coli

Cellular processes are carried out through a series of molecular interactions. Various experimental approaches can be used to investigate these functional relationships on a large-scale. Recently, the power of investigating biological systems from the perspective of genetic (gene-gene or epistatic) interactions has been evidenced by the ability to elucidate novel functional relationships. Examples of functionally related genes include genes that buffer each other's function or impinge on the same biological process. Genetic interactions have traditionally been investigated in bacteria by combining pairs of mutations (e.g., gene deletions) and assessing deviation of the phenotype of each double mutant from an expected neutral (or no interaction) phenotype. Fitness is a particularly convenient phenotype to measure: when the double mutant grows faster or slower than expected, the two mutated genes are said to show alleviating or aggravating interactions, respectively. The most commonly used neutral model assumes that the fitness of the double mutant is equal to the product of individual single mutant fitness. A striking genetic interaction is exemplified by the loss of two nonessential genes that buffer each other in performing an essential biological function: deleting only one of these genes produces no detectable fitness defect; however, loss of both genes simultaneously results in systems failure, leading to synthetic sickness or lethality. Systematic large-scale genetic interaction screens have been used to generate functional maps for model eukaryotic organisms, such as yeast, to describe the functional organization of gene products into pathways and protein complexes within a cell. They also reveal the modular arrangement and cross talk of pathways and complexes within broader functional neighborhoods (Dixon et al., Annu Rev Genet 43:601-625, 2009). Here, we present a high-throughput quantitative Escherichia coli Synthetic Genetic Array (eSGA) screening procedure, which we developed to systematically infer genetic interactions by scoring growth defects among large numbers of double mutants in a classic Gram-negative bacterium. The eSGA method exploits the rapid colony growth, ease of genetic manipulation, and natural efficient genetic exchange via conjugation of laboratory E. coli strains. Replica pinning is used to grow and mate arrayed sets of single gene mutant strains and to select double mutants en masse. Strain fitness, which is used as the eSGA readout, is quantified by the digital imaging of the plates and subsequent measuring and comparing single and double mutant colony sizes. While eSGA can be used to screen select mutants to probe the functions of individual genes, using eSGA more broadly to collect genetic interaction data for many combinations of genes can help reconstruct a functional interaction network to reveal novel links and components of biological pathways as well as unexpected connections between pathways. A variety of bacterial systems can be investigated, wherein the genes impinge on a essential biological process (e.g., cell wall assembly, ribosome biogenesis, chromosome replication) that are of interest from the perspective of drug development (Babu et al., Mol Biosyst 12:1439-1455, 2009). We also show how genetic interactions generated by high-throughput eSGA screens can be validated by manual small-scale genetic crosses and by genetic complementation and gene rescue experiments.

Mobile genetic element proliferation and gene inactivation impact over the genome structure and metabolic capabilities of Sodalis glossinidius, the secondary endosymbiont of tsetse flies

BACKGROUND

Genome reduction is a common evolutionary process in symbiotic and pathogenic bacteria. This process has been extensively characterized in bacterial endosymbionts of insects, where primary mutualistic bacteria represent the most extreme cases of genome reduction consequence of a massive process of gene inactivation and loss during their evolution from free-living ancestors. Sodalis glossinidius, the secondary endosymbiont of tsetse flies, contains one of the few complete genomes of bacteria at the very beginning of the symbiotic association, allowing to evaluate the relative impact of mobile genetic element proliferation and gene inactivation over the structure and functional capabilities of this bacterial endosymbiont during the transition to a host dependent lifestyle.

RESULTS

A detailed characterization of mobile genetic elements and pseudogenes reveals a massive presence of different types of prophage elements together with five different families of IS elements that have proliferated across the genome of Sodalis glossinidius at different levels. In addition, a detailed survey of intergenic regions allowed the characterization of 1501 pseudogenes, a much higher number than the 972 pseudogenes described in the original annotation. Pseudogene structure reveals a minor impact of mobile genetic element proliferation in the process of gene inactivation, with most of pseudogenes originated by multiple frameshift mutations and premature stop codons. The comparison of metabolic profiles of Sodalis glossinidius and tsetse fly primary endosymbiont Wiglesworthia glossinidia based on their whole gene and pseudogene repertoires revealed a novel case of pathway inactivation, the arginine biosynthesis, in Sodalis glossinidius together with a possible case of metabolic complementation with Wigglesworthia glossinidia for thiamine biosynthesis.

CONCLUSIONS

The complete re-analysis of the genome sequence of Sodalis glossinidius reveals novel insights in the evolutionary transition from a free-living ancestor to a host-dependent lifestyle, with a massive proliferation of mobile genetic elements mainly of phage origin although with minor impact in the process of gene inactivation that is taking place in this bacterial genome. The metabolic analysis of the whole endosymbiotic consortia of tsetse flies have revealed a possible phenomenon of metabolic complementation between primary and secondary endosymbionts that can contribute to explain the co-existence of both bacterial endosymbionts in the context of the tsetse host.