Experimental demonstration of operon formation catalyzed by insertion sequence.

Operons are a hallmark of the genomic and regulatory architecture of prokaryotes. However, the mechanism by which two genes placed far apart gradually come close and form operons remains to be elucidated. Here, we propose a new model of the origin of operons: Mobile genetic elements called insertion sequences can facilitate the formation of operons by consecutive insertion-deletion-excision reactions. This mechanism barely leaves traces of insertion sequences and thus difficult to detect in nature. In this study, as a proof-of-concept, we reproducibly demonstrated operon formation in the laboratory. The insertion sequence IS3 and the insertion sequence excision enhancer are genes found in a broad range of bacterial species. We introduced these genes into insertion sequence-less Escherichia coli and found that, supporting our hypothesis, the activity of the two genes altered the expression of genes surrounding IS3, closed a 2.7 kb gap between a pair of genes, and formed new operons. This study shows how insertion sequences can facilitate the rapid formation of operons through locally increasing the structural mutation rates and highlights how coevolution with mobile elements may shape the organization of prokaryotic genomes and gene regulation.

Enhancement of K-strategy evolution in histidine utilization using a container with compartments.

Evolutionary strategies in growth improvement can be classified into r- or K-strategies. The former strategy corresponds to an evolutionary increase in growth rate, whereas the latter corresponds to an increase in the maximum amount of organisms or carrying capacity. What determines the strategies to be adopted during evolution? Spatial structures that compartmentalize the population into small patches are key to inducing the K-strategy. Interestingly, previous evolution experiments using Escherichia coli in a glucose-limited batch culture showed that carrying capacity could improve evolutionally even in the absence of spatial structures. However, it is unclear if the lack of spatial structures can direct evolution toward high carrying capacity for utilization of other resources. To address this question, we established a simplified evolution experiment using histidine-requiring E. coli grown under histidine limitation in a container with compartments. We confirmed the importance of spatial structures in K-strategy evolution in histidine utilization. Whole genome sequencing of the K-adapted strains showed functional variety of the mutated genes during the fitness-increasing period. These results validate the importance of spatial structures and imply that restriction of K-strategy evolution on a sort of nutrients is attributable to a paucity of appropriate selection rather than a paucity of causal mutation.

Molecular Clock of Neutral Mutations in a Fitness-Increasing Evolutionary Process.

The molecular clock of neutral mutations, which represents linear mutation fixation over generations, is theoretically explained by genetic drift in fitness-steady evolution or hitchhiking in adaptive evolution. The present study is the first experimental demonstration for the molecular clock of neutral mutations in a fitness-increasing evolutionary process. The dynamics of genome mutation fixation in the thermal adaptive evolution of Escherichia coli were evaluated in a prolonged evolution experiment in duplicated lineages. The cells from the continuously fitness-increasing evolutionary process were subjected to genome sequencing and analyzed at both the population and single-colony levels. Although the dynamics of genome mutation fixation were complicated by the combination of the stochastic appearance of adaptive mutations and clonal interference, the mutation fixation in the population was simply linear over generations. Each genome in the population accumulated 1.6 synonymous and 3.1 non-synonymous neutral mutations, on average, by the spontaneous mutation accumulation rate, while only a single genome in the population occasionally acquired an adaptive mutation. The neutral mutations that preexisted on the single genome hitchhiked on the domination of the adaptive mutation. The successive fixation processes of the 128 mutations demonstrated that hitchhiking and not genetic drift were responsible for the coincidence of the spontaneous mutation accumulation rate in the genome with the fixation rate of neutral mutations in the population. The molecular clock of neutral mutations to the fitness-increasing evolution suggests that the numerous neutral mutations observed in molecular phylogenetic trees may not always have been fixed in fitness-steady evolution but in adaptive evolution.

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.

Nutrient-dependent growth defects and mutability of mutators in Escherichia coli.

So-called mutators emerge when mismatch repair and proofreading mechanisms are defective. Mutators not only accelerate the accumulation of mutations that are beneficial for adaptation but also cause a large number of deleterious mutations that are disadvantageous for cell growth. However, such growth defects may be compensated by nutrient availability. How the growth burden is associated with high mutability in relation to nutritional variation is an intriguing question. To address this question, we constructed a variety of Escherichia coli mutator strains through combinatorial deletions of mismatch repair and proofreading genes and quantitatively evaluated their growth and mutation rates under different nutritional conditions. Growth defects caused by high mutation rates were commonly observed in all mutators, and these defects were alleviated by nutrient supplementation in most mutators. In addition, the mutation rates of the mutators fluctuated greatly in response to nutritional conditions, in contrast to the nearly constant mutation rate of the wild-type strain under varying nutritional conditions. The results showed conditional growth defects and nutrition-sensitive mutability as general features of mutators. This study indicates the importance of modulating mutability in response to changing nutrient conditions to minimize the risk of extinction due to genetic load.

Genomic confirmation of nutrient-dependent mutability of mutators in Escherichia coli.

Mutators with increased mutation rates are prevalent in various environments and have important roles in accelerating adaptive evolution. Previous studies on mutator strains of microorganisms have shown that some mutators have constant mutation rates, whereas others exhibit switchable mutation rates depending on nutritional conditions. This suggests that the contributions of mutators on evolution vary with fluctuating nutritional conditions. However, such conditional mutability has been unclear at the genomic level. In addition, it is still unknown why mutation rates change with nutritional condition. Here, we used two mutator strains of Escherichia coli to explore the nutrient dependence of mutation rates at the genomic level. These strains were transferred repeatedly under different nutritional conditions for hundreds of generations to accumulate mutations. Whole-genome sequencing of the offspring showed that the nutrient dependence of the mutation rates was pervasive at the genomic scale. Neutrality in the mutation accumulation processes and constancy in the mutational bias suggested that nutrient dependence was not derived from conditional selective purges or from shifts of mutational bias. Some mutators could simply switch their mutation rates for both transitions and transversions in response to nutritional shifts.

Directed evolution of cell size in Escherichia coli.

BACKGROUND: In bacteria, cell size affects chromosome replication, the assembly of division machinery, cell wall synthesis, membrane synthesis and ultimately growth rate. In addition, cell size can also be a target for Darwinian evolution for protection from predators. This strong coupling of cell size and growth, however, could lead to the introduction of growth defects after size evolution. An important question remains: can bacterial cell size change and/or evolve without imposing a growth burden? RESULTS: The directed evolution of particular cell sizes, without a growth burden, was tested with a laboratory Escherichia coli strain. Cells of defined size ranges were collected by a cell sorter and were subsequently cultured. This selection-propagation cycle was repeated, and significant changes in cell size were detected within 400 generations. In addition, the width of the size distribution was altered. The changes in cell size were unaccompanied by a growth burden. Whole genome sequencing revealed that only a few mutations in genes related to membrane synthesis conferred the size evolution. CONCLUSIONS: In conclusion, bacterial cell size could evolve, through a few mutations, without growth reduction. The size evolution without growth reduction suggests a rapid evolutionary change to diverse cell sizes in bacterial survival strategies.

A reduced genome decreases the host carrying capacity for foreign DNA.

BACKGROUND: Host-plasmid interactions have been discussed largely in terms of the influences of plasmids, whereas the contributions of variations in host genomes to host interactions with foreign DNA remain unclear. A strain with a so-called "clean genome" (i.e., MDS42) of reduced genome size has recently been generated from the wild-type strain MG1655, a commonly used host strain. A quantitative evaluation of the influence of plasmid burdens in these two Escherichia coli strains can not only provide an understanding of how a reduced genome responds to foreign DNA but also offer insights into the proper application of these strains. RESULTS: The decreases in growth caused by the cost of carrying foreign DNA were similar for the wild-type and clean-genome strains. A negative correlation between the growth rate and the total amount of exogenous DNA was observed in both strains, but a better theoretical fit with a higher statistical significance was found for the strain with the clean genome. Compared to the wild-type strain, the clean-genome strain exhibited a reduced carrying capacity for exogenous DNA, which was largely attributed to its ability to restrict the replication of foreign DNA. A tendency to allocate energy and resources toward gene expression, but not DNA replication, was observed in the strain with the clean genome. CONCLUSIONS: The possession of a clean genome constrained the plasmid copy number to a wild-type-equivalent load. The results indicate that the wild-type strain possesses a greater tolerance for foreign DNA, as in endosymbiosis, and that the use of strains with clean genomes will be favorable in the applications that require precise control and theoretical prediction.

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.

Inference of transcriptome signatures of Escherichia coli in long-term stationary phase.

"Non-growing" is a dominant life form of microorganisms in nature, where available nutrients and resources are limited. In laboratory culture systems, Escherichia coli can survive for years under starvation, denoted as long-term stationary phase, where a small fraction of cells manages to survive by recycling resources released from nonviable cells. Although the physiology by which viable cells in long-term stationary phase adapt to prolonged starvation is of great interest, their genome-wide response has not been fully understood. In this study, we analyzed transcriptional profiles of cells exposed to the supernatant of 30-day long-term stationary phase culture and found that their transcriptome profiles displayed several similar responses to those of cells in the 16-h short-term stationary phase. Nevertheless, our results revealed that cells in long-term stationary phase supernatant exhibit higher expressions of stress-response genes such as phage shock proteins (psp), and lower expressions of growth-related genes such as ribosomal proteins than those in the short-term stationary phase. We confirmed that the mutant lacking the psp operon showed lower survival and growth rate in the long-term stationary phase culture. This study identified transcriptional responses for stress-resistant physiology in the long-term stationary phase environment.

Synthetic symbiosis between a cyanobacterium and a ciliate toward novel chloroplast-like endosymbiosis.

Chloroplasts are thought to have co-evolved through endosymbiosis, after a cyanobacterial-like prokaryote was engulfed by a eukaryotic cell; however, it is impossible to observe the process toward chloroplasts. In this study, we constructed an experimental symbiosis model to observe the initial stage in the process from independent organisms to a chloroplast-like organelle. Our system of synthetic symbiosis is capable of long-term coculture of two model organisms: a cyanobacterium (Synechocystis sp. PCC6803) as a symbiont and a ciliate (Tetrahymena thermophila) as a host with endocytic ability. The experimental system was clearly defined, because we used a synthetic medium and the cultures were shaken to avoid spatial complexity. We determined the experimental conditions for sustainable coculture, by analyzing population dynamics using a mathematical model. We experimentally demonstrated that the coculture was sustainable for at least 100 generations, through serial transfers. Moreover, we found that cells isolated after the serial transfer improved the probability of coexistence of both species without extinction in re-coculture. The constructed system will be useful for understanding the initial stage of primary endosymbiosis from cyanobacteria to chloroplasts, i.e., the origin of algae and plants.

Adaptation by stochastic switching of a monostable genetic circuit in Escherichia coli.

Stochastic switching is considered as a cost-saving strategy for adaptation to environmental challenges. We show here that stochastic switching of a monostable circuit can mediate the adaptation of the engineered OSU12-hisC Escherichia coli strain to histidine starvation. In this strain, the hisC gene was deleted from the His operon and placed under the control of a monostable foreign promoter. In response to histidine depletion, the OSU12-hisC population shifted to a higher HisC expression level, which is beneficial under starving conditions but is not favoured by the monostable circuit. The population shift was accompanied by growth recovery and was reversible upon histidine addition. A weak directionality in stochastic switching of hisC was observed in growing microcolonies under histidine-free conditions. Directionality and fate decision were in part dependent on the initial cellular status. Finally, microarray analysis indicated that OSU12-hisC reorganized its transcriptome to reach the appropriate physiological state upon starvation. These findings suggest that bacteria do not necessarily need to evolve signalling mechanisms to control gene expression appropriately, even for essential genes.

Purifying selection enduringly acts on the sequence evolution of highly expressed proteins in Escherichia coli.

The evolutionary speed of a protein sequence is constrained by its expression level, with highly expressed proteins evolving relatively slowly. This negative correlation between expression levels and evolutionary rates (known as the E-R anticorrelation) has already been widely observed in past macroevolution between species from bacteria to animals. However, it remains unclear whether this seemingly general law also governs recent evolution, including past and de novo, within a species. However, the advent of genomic sequencing and high-throughput phenotyping, particularly for bacteria, has revealed fundamental gaps between the 2 evolutionary processes and has provided empirical data opposing the possible underlying mechanisms which are widely believed. These conflicts raise questions about the generalization of the E-R anticorrelation and the relevance of plausible mechanisms. To explore the ubiquitous impact of expression levels on molecular evolution and test the relevance of the possible underlying mechanisms, we analyzed the genome sequences of 99 strains of Escherichia coli for evolution within species in nature. We also analyzed genomic mutations accumulated under laboratory conditions as a model of de novo evolution within species. Here, we show that E-R anticorrelation is significant in both past and de novo evolution within species in E. coli. Our data also confirmed ongoing purifying selection on highly expressed genes. Ongoing selection included codon-level purifying selection, supporting the relevance of the underlying mechanisms. However, the impact of codon-level purifying selection on the constraints in evolution within species might be smaller than previously expected from evolution between species.

Bacterial cells carrying synthetic dual-function operon survived starvation.

A synthetic dual-function operon with a bistable structure was designed and successfully integrated into the bacterial genome. Bistability was generated by the mutual inhibitory structure comprised of the promoters P(tet) and P(lac) and the repressors LacI and TetR. Dual function essential for cell growth was introduced by replacing the genes (i.e., hisC and leuB) encoding proteins involved in the biosynthesis of histidine and leucine from their native chromosomal locations to the synthetic operon. Both colony formation and population dynamics of the cells carrying this operon showed that the cells survived starvation and the newly formed population transited between the two stable states, representing the induced hisC and leuB levels, in accordance with the nutritional status. The results strongly suggested that the synthetic design of proto-operons sensitive to external perturbations is practical and functional in native cells.

Construction of Escherichia coli gene expression level perturbation collection.

We generated 61 strains of Escherichia coli in which the expression level of a specific single gene can be changed continuously over a physiologically significant range. In each strain, one auxotrophic gene was deleted from its original position and reinserted at a specific position on the chromosome under the control of the tetA promoter. Therefore, the level of expression of the target gene can be controlled easily by altering the concentrations of inducers, e.g., anhydrotetracycline and doxycycline, in the medium. Protein and mRNA levels and changes in proliferation rate were examined in some of the strains in our collection to determine the ability to control the level of target gene expression over a physiologically significant range. These strains will be useful for extracting omics data sets and for the construction of genome-scale mathematical models, because causality between perturbations in gene expression level and their consequences can be clearly determined.

Noisy cell growth rate leads to fluctuating protein concentration in bacteria.

The present study discusses a prime cause of fluctuating protein concentrations, which play a significant role in generating phenotypic diversity in bacteria. A genetic circuit integrated in a bacterial genome was used to evaluate the cell-to-cell variation in protein concentration. A simple dynamic model, comprising terms for synthesis and dilution, was used to elucidate the contributions of distinct noises to the fluctuation in cell protein concentration. Experimental and theoretical results demonstrated that noise in the rate of increase in cell volume (cell growth rate) serves as a source of extrinsic noise that accounts for dozens of percent of the total noise, whereas intrinsic noise in protein synthesis makes only a moderate contribution to the fluctuation in protein concentration. This suggests that such external noise in the cell growth rate has a global effect on cellular components, resulting in a large fluctuation in protein concentration in bacterial cells.

Development of an Automated UV Irradiation Device for Microbial Cell Culture.

Ultraviolet (UV) mutagenesis is a widely used technique to increase bacterial mutation rates in laboratory experiments. UV mutagenesis requires fine regulation of UV dose, because the number of dead cells increases exponentially as the dose increases. Ignoring this hazard can cause extinction of UV-exposed populations. Therefore, an automated system that cooperatively conducts both growth measurement and UV irradiation is needed for efficient UV mutagenesis experiments. To address this task, we constructed an automated UV irradiation device for microbial cell culture. This device can measure cell density and irradiate the bacterial cells with UV light automatically according to the state of cell growth. We demonstrated that this growth feedback control avoided extinction and enabled accumulation of mutations in bacterial genomes at a rapid rate for a long period. Whole-genome sequencing revealed the high accumulation rate, neutrality, and spectrum of UV-induced mutations. These characteristics were all consistent with those obtained by manual UV irradiation. These results indicate that our automated device is useful in accelerating mutation accumulation over a long duration.

Density-Dependent Recycling Promotes the Long-Term Survival of Bacterial Populations during Periods of Starvation.

The amount of natural resources in the Earth's environment is in flux, which can trigger catastrophic collapses of ecosystems. How populations survive under nutrient-poor conditions is a central question in ecology. Curiously, some bacteria persist for a long time in nutrient-poor environments. Although this survival may be accomplished through cell death and the recycling of dead cells, the importance of these processes and the mechanisms underlying the survival of the populations have not been quantitated. Here, we use microbial laboratory experiments and mathematical models to demonstrate that death and recycling are essential activities for the maintenance of cell survival. We also show that the behavior of the survivors is governed by population density feedback, wherein growth is limited not only by the available resources but also by the population density. The numerical simulations suggest that population density-dependent recycling could be an advantageous behavior under starvation conditions. IMPORTANCE: How organisms survive after exhaustion of resources is a central question in ecology. Starving Escherichia coli constitute a model system to understand survival mechanisms during long-term starvation. Although death and the recycling of dead cells might play a key role in the maintenance of long-term survival, their mechanisms and importance have not been quantitated. Here, we verified the significance of social recycling of dead cells for long-term survival. We also show that the survivors restrained their recycling and did not use all available nutrients released from dead cells, which may be advantageous under starvation conditions. These results indicate that not only the utilization of dead cells but also restrained recycling coordinate the effective utilization of limited resources for long-term survival under starvation.

Mutation accumulation under UV radiation in Escherichia coli.

Mutations are induced by not only intrinsic factors such as inherent molecular errors but also by extrinsic mutagenic factors such as UV radiation. Therefore, identifying the mutational properties for both factors is necessary to achieve a comprehensive understanding of evolutionary processes both in nature and in artificial situations. Although there have been extensive studies on intrinsic factors, the mutational profiles of extrinsic factors are poorly understood on a genomic scale. Here, we explored the mutation profiles of UV radiation, a ubiquitous mutagen, in Escherichia coli on the genomic scale. We performed an evolution experiment under periodic UV radiation for 28 days. The accumulation speed of the mutations was found to increase so that it exceeded that of a typical mutator strain with deficient mismatch repair processes. The huge contribution of the extrinsic factors to all mutations consequently increased the risk of the destruction of inherent error correction systems. The spectrum of the UV-induced mutations was broader than that of the spontaneous mutations in the mutator. The broad spectrum and high upper limit of the frequency of occurrence suggested ubiquitous roles for UV radiation in accelerating the evolutionary process.

Stochasticity in gene expression in a cell-sized compartment.

The gene expression in a clonal cell population fluctuates significantly, and its relevance to various cellular functions is under intensive debate. A fundamental question is whether the fluctuation is a consequence of the complexity and redundancy in living cells or an inevitable attribute of the minute microreactor nature of cells. To answer this question, we constructed an artificial cell, which consists of only necessary components for the gene expression (in vitro transcription and translation system) and its boundary as a microreactor (cell-sized lipid vesicle), and investigated the gene expression noise. The variation in the expression of two fluorescent proteins was decomposed into the components that were correlated and uncorrelated between the two proteins using a method similar to the one used by Elowitz and co-workers to analyze the expression noise in E. coli. The observed fluctuation was compared with a theoretical model that expresses the amplitude of noise as a function of the average number of intermediate molecules and products. With the assumption that the transcripts are partly active, the theoretical model was able to well describe the noise in the artificial system. Furthermore, the same measurement for E. coli cells harboring an identical plasmid revealed that the E. coli exhibited a similar level of expression noise. Our results demonstrated that the level of fluctuation found in bacterial cells is mostly an intrinsic property that arises even in a primitive form of the cell.