Search

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. Finally, taking an integrated view, we show that across all phenotypes combined, segregating promoter variants are far more phenotypically similar than would be expected given their genetic divergence. This is the consequence of both stabilizing and directional selection acting on individual phenotypes to minimize differences among segregating variants. 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.

Escherichia coli is commonly considered a host-associated bacterium. However, there is evidence that some strains occupy environmental (non-host-associated) niches. Here, we report the complete genomes of 47 Escherichia coli environmental isolates. These will be useful for understanding the dynamics of plasmids, phages, and other repetitive genetic elements.

Population genetic structure in the marine environment can be influenced by life‐history traits such as developmental mode (biphasic, with distinct adult and larval morphology, and direct development, in which larvae resemble adults) or habitat specificity, as well as geography and selection. Developmental mode is thought to significantly influence dispersal, with direct developers expected to have much lower dispersal potential. However, this prediction can be complicated by the presence of geophysical barriers to dispersal. In this study, we use a panel of 8,020 SNPs to investigate population structure and biogeography over multiple spatial scales for a direct‐developing species, the New Zealand endemic marine isopod Isocladus armatus. Because our sampling range is intersected by two well‐known biogeographic barriers (the East Cape and the Cook Strait), our study provides an opportunity to understand how such barriers influence dispersal in direct developers. On a small spatial scale (20 km), gene flow between locations is extremely high, suggestive of an island model of migration. However, over larger spatial scales (600 km), populations exhibit a clear pattern of isolation‐by‐distance. Our results indicate that I. armatus exhibits significant migration across the hypothesized barriers and suggest that large‐scale ocean currents associated with these locations do not present a barrier to dispersal. Interestingly, we find evidence of a north‐south population genetic break occurring between Māhia and Wellington. While no known geophysical barrier is apparent in this area, it coincides with the location of a proposed border between bioregions. Analysis of loci under selection revealed that both isolation‐by‐distance and adaption may be contributing to the degree of population structure we have observed here. We conclude that developmental life history largely predicts dispersal in the intertidal isopod I. armatus. However, localized biogeographic processes can disrupt this expectation, and this may explain the potential meta‐population detected in the Auckland region., We show, for a New Zealand isopod, a strong relationship between geographic and genetic relationships. Isocladus armatus exhibits both isolation‐by‐distance and geographic concordance with a PCA. Unusually, we find evidence of high and ongoing gene flow across a well studied biogeographic break, which in other species hinders connectivity.

Beatrix, Carthage, Daegal, Dulcie, Fancypants, Fenn, Inca, Naira, and Robyn are newly isolated bacteriophages capable of infecting Mycolicibacterium smegmatis mc2 155. We discovered, sequenced, and annotated these New Zealand bacteriophages., Beatrix, Carthage, Daegal, Dulcie, Fancypants, Fenn, Inca, Naira, and Robyn are newly isolated bacteriophages capable of infecting Mycolicibacterium smegmatis mc2 155. We discovered, sequenced, and annotated these New Zealand bacteriophages. These phages illustrate that New Zealand harbors a selection of the highly diverse and distributed mycobacteriophage clusters found globally.

Seven mycobacteriophages from distinct geographical locations were isolated, using Mycobacterium smegmatis mc 2 155 as the host, and then purified and sequenced. All of the genomes are related to cluster A mycobacteriophages, BobSwaget and Lokk in subcluster A2; Fred313, KADY, Stagni, and StepMih in subcluster A3; and MyraDee in subcluster A18, the first phage to be assigned to that subcluster.

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.

Bacteriophages are the most numerous entities in the biosphere. Despite this numerical dominance, the genetic structure of bacteriophage populations is poorly understood. Here, we present a biogeography study involving 25 previously undescribed bacteriophages from the Cystoviridae clade, a group characterized by a dsRNA genome divided into three segments. Previous laboratory manipulation has shown that, when multiple Cystoviruses infect a single host cell, they undergo ( i ) rare intrasegment recombination events and ( ii ) frequent genetic reassortment between segments. Analyzing linkage disequilibrium (LD) within segments, we find no significant evidence of intrasegment recombination in wild populations, consistent with ( i ). An extensive analysis of LD between segments supports frequent reassortment, on a time scale similar to the genomic mutation rate. The absence of LD within this group of phages is consistent with expectations for a completely sexual population, despite the fact that some segments have >50% nucleotide divergence at 4-fold degenerate sites. This extraordinary rate of genetic exchange between highly unrelated individuals is unprecedented in any taxa. We discuss our results in light of the biological species concept applied to viruses.

Although it is often tacitly assumed that gene regulatory interactions are finely tuned, how accurate gene regulation could evolve from a state without regulation is unclear. Moreover, gene expression noise would seem to impede the evolution of accurate gene regulation, and previous investigations have provided circumstantial evidence that natural selection has acted to lower noise levels. By evolving synthetic Escherichia coli promoters de novo, we here show that, contrary to expectations, promoters exhibit low noise by default. Instead, selection must have acted to increase the noise levels of highly regulated E. coli promoters. We present a general theory of the interplay between gene expression noise and gene regulation that explains these observations. The theory shows that propagation of expression noise from regulators to their targets is not an unwanted side-effect of regulation, but rather acts as a rudimentary form of regulation that facilitates the evolution of more accurate regulation. DOI: http://dx.doi.org/10.7554/eLife.05856.001, eLife digest Genes are stretches of DNA that contain the instructions needed to make proteins and other molecules. By changing how much protein is produced from each gene (i.e., its expression), many organisms—including humans—can produce a wide variety of cell types with very different behaviors. Similarly, single-celled organisms, such as bacteria, can adapt to survive and grow in different environments by changing gene expression levels. It is thus thought that gene expression must be precisely controlled. However, the molecular processes involved in gene expression are subject to random fluctuations, and so gene expression is inherently ‘noisy’. This means that even groups of identical cells in identical environments will show variation in their gene expression patterns. Furthermore, different genes show different levels of noise. The DNA sequence of a part of each gene, called the promoter, has a big effect on these noise levels. Consequently, gene expression noise is a genetically encoded trait, and can therefore be shaped by natural selection. But it remains largely unclear how natural selection has affected gene expression noise. Now, Wolf et al. have carefully measured the gene expression noise of hundreds of synthetic promoters that were evolved in the laboratory from random DNA sequences, and a similar number of natural promoters in a bacterium called E. coli. These experiments revealed that, contrary to expectation, most lab-evolved promoters had low levels of noise. On the other hand, many natural promoters had high levels of noise. Wolf et al. also found that noisy promoters tend to be highly regulated by transcription factors: the proteins that control gene expression by binding to promoter regions. Together, these results imply that unregulated promoters start by having low noise as a default state. Selection pressures must then have caused some E. coli promoters to become regulated by transcription factors and raise their noise levels. But, what might these selection pressures have been? Many genes need to be expressed at different levels in different conditions, and it is generally accepted that regulation by transcription factors evolves to ‘satisfy’ these requirements. However, transcription factors are themselves noisy, and this noise necessarily propagates to their target genes. Wolf et al. have now developed a general theory showing that this noise-propagation can often benefit an organism. This explains why natural selection can favor an increase in noise levels for regulated genes. Importantly, by showing that the main role of a transcription factor can be to increase the noise of its targets, it suddenly becomes very easy to see how new gene regulatory interactions can evolve from scratch. The next steps in understanding of how gene expression noise evolves will involve manipulating the expression noise of a gene, and measuring how selection acts on such changes. DOI: http://dx.doi.org/10.7554/eLife.05856.002

Studies of microbial evolutionary dynamics are being transformed by the availability of affordable high-throughput sequencing technologies, which allow whole-genome sequencing of hundreds of related taxa in a single study. Reconstructing a phylogenetic tree of these taxa is generally a crucial step in any evolutionary analysis. Instead of constructing genome assemblies for all taxa, annotating these assemblies, and aligning orthologous genes, many recent studies 1) directly map raw sequencing reads to a single reference sequence, 2) extract single nucleotide polymorphisms (SNPs), and 3) infer the phylogenetic tree using maximum likelihood methods from the aligned SNP positions. However, here we show that, when using such methods to reconstruct phylogenies from sets of simulated sequences, both the exclusion of nonpolymorphic positions and the alignment to a single reference genome, introduce systematic biases and errors in phylogeny reconstruction. To address these problems, we developed a new method that combines alignments from mappings to multiple reference sequences and show that this successfully removes biases from the reconstructed phylogenies. We implemented this method as a web server named REALPHY (Reference sequence Alignment-based Phylogeny builder), which fully automates phylogenetic reconstruction from raw sequencing reads., Molecular Biology and Evolution, 31 (5), ISSN:0737-4038, ISSN:1537-1719

X-linked lymphoproliferative disease (XLP) is a primary immunodeficiency characterized by extreme susceptibility to Epstein-Barr virus. The XLP disease gene product SH2D1A (SAP) interacts via its SH2 domain with a motif (TIYXXV) present in the cytoplasmic tail of the cell-surface receptors CD150/SLAM, CD84, CD229/Ly-9, and CD244/2B4. Characteristically, the SH2D1A three-pronged interaction with Tyr(281) of CD150 can occur in absence of phosphorylation. Here we analyze the effect of SH2D1A protein missense mutations identified in 10 XLP families. Two sets of mutants were found: (i) mutants with a marked decreased protein half-life (e.g. Y7C, S28R, Q99P, P101L, V102G, and X129R) and (ii) mutants with structural changes that differently affect the interaction with the four receptors. In the second group, mutations that disrupt the interaction between the SH2D1A hydrophobic cleft and Val +3 of its binding motif (e.g. T68I) and mutations that interfere with the SH2D1A phosphotyrosine-binding pocket (e.g. C42W) abrogated SH2D1A binding to all four receptors. Surprisingly, a mutation in SH2D1A able to interfere with Thr -2 of the CD150 binding motif (mutant T53I) severely impaired non-phosphotyrosine interactions while preserving unaffected the binding of SH2D1A to phosphorylated CD150. Mutant T53I, however, did not bind to CD229 and CD224, suggesting that SH2D1A controls several critical signaling pathways in T and natural killer cells. Because no correlation is present between identified types of mutations and XLP patient clinical presentation, additional unidentified genetic or environmental factors must play a strong role in XLP disease manifestations.

Figure S6 is incorrectly a duplicate of Figure S7. The correct Figure S6 can be viewed here: Click here for additional data file.(244K, pdf) [^]

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., PLoS Genetics, 8 (1), ISSN:1553-7390, ISSN:1553-7404

The genetic structure of natural bacteriophage populations is poorly understood. Recent metagenomic studies suggest that phage biogeography is characterized by frequent migration. Using virus samples mostly isolated in Southern California, we recently showed that very little population structure exists in segmented RNA phage of the Cystoviridae family due to frequent segment reassortment (sexual genetic mixis) between unrelated virus individuals. Here we use a larger genetic dataset to examine the structure of Cystoviridae phage isolated from three geographic locations in Southern New England. We document extensive natural variation in the physical sizes of RNA genome segments for these viruses. In addition, consistent with earlier findings, our phylogenetic analyses and calculations of linkage disequilibrium (LD) show no evidence of within-segment recombination in wild populations. However, in contrast to the prior study, our analysis finds that reassortment of segments between individual phage plays a lesser role among cystoviruses sampled in New England, suggesting that the evolutionary importance of genetic mixis in Cystoviridae phage may vary according to geography. We discuss possible explanations for these conflicting results across the studies, such as differing local ecology and its impact on phage growth, and geographic differences in selection against hybrid phage genotypes.

Genetically identical populations of unicellular organisms often show marked variation in some phenotypic traits. To investigate the molecular causes and possible biological functions of this phenotypic noise, it would be useful to have a method to identify genes whose expression varies stochastically on a certain time scale. Here, we developed such a method and used it for identifying genes with high levels of phenotypic noise in Salmonella enterica ssp. I serovar Typhimurium (S. Typhimurium). We created a genomic plasmid library fused to a green fluorescent protein (GFP) reporter and subjected replicate populations harboring this library to fluctuating selection for GFP expression using fluorescent-activated cell sorting (FACS). After seven rounds of fluctuating selection, the populations were strongly enriched for promoters that showed a high amount of noise in gene expression. Our results indicate that the activity of some promoters of S. Typhimurium varies on such a short time scale that these promoters can absorb rapid fluctuations in the direction of selection, as imposed during our experiment. The genomic fragments that conferred the highest levels of phenotypic variation were promoters controlling the synthesis of flagella, which are associated with virulence and host–pathogen interactions. This confirms earlier reports that phenotypic noise may play a role in pathogenesis and indicates that these promoters have among the highest levels of noise in the S. Typhimurium genome. This approach can be applied to many other bacterial and eukaryotic systems as a simple method for identifying genes with noisy expression., PLoS Genetics, 4, ISSN:1553-7390, ISSN:1553-7404

The most consistent result in more than two decades of experimental evolution is that the fitness of populations adapting to a constant environment does not increase indefinitely, but reaches a plateau. Using experimental evolution with bacteriophage, we show here that the converse is also true. In populations small enough such that drift overwhelms selection and causes fitness to decrease, fitness declines down to a plateau. We demonstrate theoretically that both of these phenomena must be due either to changes in the ratio of beneficial to deleterious mutations, the size of mutational effects, or both. We use mutation accumulation experiments and molecular data from experimental evolution to show that the most significant change in mutational effects is a drastic increase in the rate of beneficial mutation as fitness decreases. In contrast, the size of mutational effects changes little even as organismal fitness changes over several orders of magnitude. These findings have significant implications for the dynamics of adaptation., PLoS Biology, 5 (4), ISSN:1544-9173, ISSN:1545-7885

Epistasis and the evolution of recombination are closely intertwined: epistasis generates linkage disequilibria (i.e. statistical associations between alleles), whereas recombination breaks them up. The mutational deterministic hypothesis (MDH) states that high recombination rates are maintained because the breaking up of linkage disequilibria generated by negative epistasis enables more efficient purging of deleterious mutations. However, recent theoretical and experimental work challenges the MDH. Experimental evidence suggests that negative epistasis, required by the MDH, is relatively uncommon. On the theoretical side, population genetic models suggest that, compared with the combined effects of drift and selection, epistasis generates a negligible amount of linkage disequilibria. Here, we assess these criticisms and discuss to what extent they invalidate the MDH as an explanation for the evolution of recombination.

Background. Orthologous genes are frequently presumed to perform similar functions. However, outside of model organisms, this is rarely tested. One means of inferring changes in function is if there are changes in the level of gene conservation and selective constraint. Here we compare levels of gene conservation across three bacterial groups to test for changes in gene functionality. Findings. The level of gene conservation for different orthologous genes is highly correlated across clades, even for highly divergent groups of bacteria. These correlations do not arise from broad differences in gene functionality (e.g. informational genes vs. metabolic genes), but instead seem to result from very specific differences in gene function. Furthermore, these functional differences appear to be maintained over very long periods of time. Conclusion. These results suggest that even over broad time scales, most bacterial genes are under a nearly constant level of purifying selection, and that bacterial evolution is thus dominated by selective and functional stasis., BMC Research Notes, 2 (1), ISSN:1756-0500