Multiple pathways of selected gene amplification during adaptive mutation

Hotspot gene conversion between FKS1 and FKS2 in echinocandin resistant Candida glabrata serial isolates

Candida glabrata (Nakaseomyces glabratus) is the most common cause of drug-resistant candidemia and is associated with a high mortality rate. Only a few mechanisms of drug resistance are known in C. glabrata, predominantly involving recurrent single nucleotide polymorphisms. The importance of structural variation in acquired drug resistance is not understood. We performed comparative phenotypic and genomic analyses of six serial bloodstream isolates of C. glabrata and identified novel mutations associated with resistance to echinocandins. Critically, we identified a novel gene conversion event between the hotspot 2 regions of FKS1 and FKS2 that was associated with increased resistance to micafungin. We further analyzed 621 publicly available C. glabrata genomes and found three additional examples of structural variation involving FKS1 and FKS2. Ultimately, drug resistance in C. glabrata involves structural variants that are missed with current diagnostic methods and need to be considered when designing and implementing more effective antifungal management strategies.

A direct experimental test of Ohno's hypothesis

Gene duplication drives evolution by providing raw material for proteins with novel functions. An influential hypothesis by Ohno (1970) posits that gene duplication helps genes tolerate new mutations and thus facilitates the evolution of new phenotypes. Competing hypotheses argue that deleterious mutations will usually inactivate gene duplicates too rapidly for Ohno's hypothesis to work. We experimentally tested Ohno's hypothesis by evolving one or exactly two copies of a gene encoding a fluorescent protein in Escherichia coli through several rounds of mutation and selection. We analyzed the genotypic and phenotypic evolutionary dynamics of the evolving populations through high-throughput DNA sequencing, biochemical assays, and engineering of selected variants. In support of Ohno's hypothesis, populations carrying two gene copies displayed higher mutational robustness than those carrying a single gene copy. Consequently, the double-copy populations experienced relaxed purifying selection, evolved higher phenotypic and genetic diversity, carried more mutations and accumulated combinations of key beneficial mutations earlier. However, their phenotypic evolution was not accelerated, possibly because one gene copy rapidly became inactivated by deleterious mutations. Our work provides an experimental platform to test models of evolution by gene duplication, and it supports alternatives to Ohno's hypothesis that point to the importance of gene dosage.

Proteostasis modulates gene dosage evolution in antibiotic-resistant bacteria

Evolution of gene expression frequently drives antibiotic resistance in bacteria. We had previously (Patel and Matange, eLife, 2021) shown that, in Escherichia coli, mutations at the mgrB locus were beneficial under trimethoprim exposure and led to overexpression of dihydrofolate reductase (DHFR), encoded by the folA gene. Here, we show that DHFR levels are further enhanced by spontaneous duplication of a genomic segment encompassing folA and spanning hundreds of kilobases. This duplication was rare in wild-type E. coli. However, its frequency was elevated in a lon-knockout strain, altering the mutational landscape early during trimethoprim adaptation. We then exploit this system to investigate the relationship between trimethoprim pressure and folA copy number. During long-term evolution, folA duplications were frequently reversed. Reversal was slower under antibiotic pressure, first requiring the acquisition of point mutations in DHFR or its promoter. Unexpectedly, despite resistance-conferring point mutations, some populations under high trimethoprim pressure maintained folA duplication to compensate for low abundance DHFR mutants. We find that evolution of gene dosage depends on expression demand, which is generated by antibiotic and exacerbated by proteolysis of drug-resistant mutants of DHFR. We propose a novel role for proteostasis as a determinant of copy number evolution in antibiotic-resistant bacteria.

Inter-Relationship Between a Transcriptional Regulator of Flagella Genes cj0440c and Thiamine Metabolic Pathway in Campylobacter jejuni

Campylobacter jejuni is a major cause of gastroenteritis in humans. It has been reported that the pathogenesis of C. jejuni is closely related to the formation, adhesion, and invasion of flagella toxin in host epithelial cells. A putative transcriptional regulator, known as cj0440c, is thought to be involved in the regulation of flagellar synthesis. However, confirmation of this hypothesis requires deep insight into the regulation mechanism of cj0440c and its possible relationship with different antibiotics. Therefore, the study explained here was designed to determine the relationship and function (phenotypically and genotypically) of cj0440c in the flagellar synthesis of C. jejuni NCTC11168. The study determined the mode of expression of cj0440c and flagella-related genes under exposure to various drugs. To verify the involvement of cj0440c protein in the metabolic pathway of thiamine, an enzymatic hydrolysis experiment was performed and analyzed through the application of mass spectrometry. The overexpression vector of C. jejuni NCTC11168 was also constructed to find out whether or not target genes were regulated by cj0440c. The findings of the study showed that cj0440c and other flagella-related genes were expressed differentially under the influence of various antibiotics including erythromycin, tylosin, azithromycin, gentamicin, etimicin, enrofloxacin, gatifloxacin, tetracycline, and tigecycline. The analysis showed that the cj0440c protein did not catalyze the degradation of thiamine. In conclusion, the study aids in the understanding of the inter-relationship between the regulatory mechanism of flagella genes and the thiamine metabolic pathway.

Amplicon remodeling and genomic mutations drive population dynamics after segmental amplification

New enzymes often evolve by duplication and divergence of genes encoding enzymes with promiscuous activities that have become important in the face of environmental opportunities or challenges. Amplifications that increase the copy number of the gene under selection commonly amplify many surrounding genes. Extra copies of these coamplified genes must be removed, either during or after evolution of a new enzyme. Here we report that amplicon remodeling can begin even before mutations occur in the gene under selection. Amplicon remodeling and mutations elsewhere in the genome that indirectly increase fitness result in complex population dynamics, leading to emergence of clones that have improved fitness by different mechanisms. In this work, one of the two most successful clones had undergone two episodes of amplicon remodeling, leaving only four coamplified genes surrounding the gene under selection. Amplicon remodeling in the other clone resulted in removal of 111 genes from the genome, an acceptable solution under these selection conditions, but one that would certainly impair fitness under other environmental conditions.

Evolution of new enzymes by gene duplication and divergence

Thousands of new metabolic and regulatory enzymes have evolved by gene duplication and divergence since the dawn of life. New enzyme activities often originate from promiscuous secondary activities that have become important for fitness due to a change in the environment or a mutation. Mutations that make a promiscuous activity physiologically relevant can occur in the gene encoding the promiscuous enzyme itself, but can also occur elsewhere, resulting in increased expression of the enzyme or decreased competition between the native and novel substrates for the active site. If a newly useful activity is inefficient, gene duplication/amplification will set the stage for divergence of a new enzyme. Even a few mutations can increase the efficiency of a new activity by orders of magnitude. As efficiency increases, amplified gene arrays will shrink to provide two alleles, one encoding the original enzyme and one encoding the new enzyme. Ultimately, genomic rearrangements eliminate co-amplified genes and move newly evolved paralogs to a distant region of the genome.

Selective Inbreeding: Genetic Crosses Drive Apparent Adaptive Mutation in the Cairns-Foster System of Escherichia coli

The Escherichia coli system of Cairns and Foster employs a lac frameshift mutation that reverts rarely (10-9/cell/division) during unrestricted growth. However, when 108 cells are plated on lactose medium, the nongrowing lawn produces ∼50 Lac+ revertant colonies that accumulate linearly with time over 5 days. Revertants carry very few associated mutations. This behavior has been attributed to an evolved mechanism ("adaptive mutation" or "stress-induced mutagenesis") that responds to starvation by preferentially creating mutations that improve growth. We describe an alternative model, "selective inbreeding," in which natural selection acts during intercellular transfer of the plasmid that carries the mutant lac allele and the dinB gene for an error-prone polymerase. Revertant genome sequences show that the plasmid is more intensely mutagenized than the chromosome. Revertants vary widely in their number of plasmid and chromosomal mutations. Plasmid mutations are distributed evenly, but chromosomal mutations are focused near the replication origin. Rare, heavily mutagenized, revertants have acquired a plasmid tra mutation that eliminates conjugation ability. These findings support the new model, in which revertants are initiated by rare pre-existing cells (105) with many copies of the F'lac plasmid. These cells divide under selection, producing daughters that mate. Recombination between donor and recipient plasmids initiates rolling-circle plasmid over-replication, causing a mutagenic elevation of DinB level. A lac+ reversion event starts chromosome replication and mutagenesis by accumulated DinB. After reversion, plasmid transfer moves the revertant lac+ allele into an unmutagenized cell, and away from associated mutations. Thus, natural selection explains why mutagenesis appears stress-induced and directed.

Fail-safe genetic codes designed to intrinsically contain engineered organisms

One challenge in engineering organisms is taking responsibility for their behavior over many generations. Spontaneous mutations arising before or during use can impact heterologous genetic functions, disrupt system integration, or change organism phenotype. Here, we propose restructuring the genetic code itself such that point mutations in protein-coding sequences are selected against. Synthetic genetic systems so-encoded should fail more safely in response to most spontaneous mutations. We designed fail-safe codes and simulated their expected effects on the evolution of so-encoded proteins. We predict fail-safe codes supporting expression of 20 or 15 amino acids could slow protein evolution to ∼30% or 0% the rate of standard-encoded proteins, respectively. We also designed quadruplet-codon codes that should ensure all single point mutations in protein-coding sequences are selected against while maintaining expression of 20 or more amino acids. We demonstrate experimentally that a reduced set of 21 tRNAs is capable of expressing a protein encoded by only 20 sense codons, whereas a standard 64-codon encoding is not expressed. Our work suggests that biological systems using rationally depleted but otherwise natural translation systems should evolve more slowly and that such hypoevolvable organisms may be less likely to invade new niches or outcompete native populations.

Mutations that improve efficiency of a weak-link enzyme are rare compared to adaptive mutations elsewhere in the genome

New enzymes often evolve by gene amplification and divergence. Previous experimental studies have followed the evolutionary trajectory of an amplified gene, but have not considered mutations elsewhere in the genome when fitness is limited by an evolving gene. We have evolved a strain of Escherichia coli in which a secondary promiscuous activity has been recruited to serve an essential function. The gene encoding the ‘weak-link’ enzyme amplified in all eight populations, but mutations improving the newly needed activity occurred in only one. Most adaptive mutations occurred elsewhere in the genome. Some mutations increase expression of the enzyme upstream of the weak-link enzyme, pushing material through the dysfunctional metabolic pathway. Others enhance production of a co-substrate for a downstream enzyme, thereby pulling material through the pathway. Most of these latter mutations are detrimental in wild-type E. coli, and thus would require reversion or compensation once a sufficient new activity has evolved.

Experimental Evolution of Innovation and Novelty

How does novelty, a new, genetically based function, evolve? A compelling answer has been elusive because there are few model systems where both the genetic mechanisms generating novel functions and the ecological conditions that govern their origin and spread can be studied in detail. This review article considers what we have learned about the evolution of novelty from microbial selection experiments. This work reveals that the genetic routes to novelty can be more highly variable than standard models have led us to believe and underscores the importance of considering both genetics and ecology in this process.

Selection-Enhanced Mutagenesis of lac Genes Is Due to Their Coamplification with dinB Encoding an Error-Prone DNA Polymerase

To test whether growth limitation induces mutations, Cairns and Foster constructed an Escherichia coli strain whose mutant lac allele provides 1-2% of normal ability to use lactose. This strain cannot grow on lactose, but produces ∼50 Lac⁺ revertant colonies per 10⁸ plated cells over 5 days. About 80% of revertants carry a stable lac⁺ mutation made by the error-prone DinB polymerase, which may be induced during growth limitation; 10% of Lac⁺ revertants are stable but form without DinB; and the remaining 10% grow by amplifying their mutant lac allele and are unstably Lac⁺ Induced DinB mutagenesis has been explained in two ways: (1) upregulation of dinB expression in nongrowing cells ("stress-induced mutagenesis") or (2) selected local overreplication of the lac and dinB⁺ genes on lactose medium (selected amplification) in cells that are not dividing. Transcription of dinB is necessary but not sufficient for mutagenesis. Evidence is presented that DinB enhances reversion only when encoded somewhere on the F'lac plasmid that carries the mutant lac gene. A new model will propose that rare preexisting cells (1 in a 1000) have ∼10 copies of the F'lac plasmid, providing them with enough energy to divide, mate, and overreplicate their F'lac plasmid under selective conditions. In these clones, repeated replication of F'lac in nondividing cells directs opportunities for lac reversion and increases the copy number of the dinB⁺ gene. Amplification of dinB⁺ increases the error rate of replication and increases the number of lac⁺ revertants. Thus, reversion is enhanced in nondividing cells not by stress-induced mutagenesis, but by selected coamplification of the dinB and lac genes, both of which happen to lie on the F'lac plasmid.

The role of recombination in evolutionary adaptation of Escherichia coli to a novel nutrient

The benefits and detriments of recombination for adaptive evolution have been studied both theoretically and experimentally, with conflicting predictions and observations. Most pertinent experiments examine recombination's effects in an unchanging environment and do not study its genomewide effects. Here, we evolved six replicate populations of either highly recombining R⁺ or lowly recombining R^- E. coli strains in a changing environment, by introducing the novel nutrients L-arabinose or indole into the environment. The experiment's ancestral strains are not viable on these nutrients, but 130 generations of adaptive evolution were sufficient to render them viable. Recombination conferred a more pronounced advantage to populations adapting to indole. To study the genomic changes associated with this advantage, we sequenced the genomes of 384 clones isolated from selected replicates at the end of the experiment. These genomes harbour complex changes that range from point mutations to large-scale DNA amplifications. Among several candidate adaptive mutations, those in the tryptophanase regulator tnaC stand out, because the tna operon in which it resides has a known role in indole metabolism. One of the highly recombining populations also shows a significant excess of large-scale segmental DNA amplifications that include the tna operon. This lineage also shows a unique and potentially adaptive combination of point mutations and DNA amplifications that may have originated independently from one another, to be joined later by recombination. Our data illustrate that the advantages of recombination for adaptive evolution strongly depend on the environment and that they can be associated with complex genomic changes.

Pathoadaptive Mutations in Uropathogenic Escherichia coli

Uropathogenic Escherichia coli (UPEC) are opportunistic human pathogens that primarily circulate as part of commensal intestinal microbiota. Though they have the ability to survive and proliferate in various urinary tract compartments, the urinary tract is a transient, occasional habitat for UPEC. Because of this, most of the UPEC traits have originally evolved to serve in intestinal colonization and transmission. Some of these bacterial traits serve as virulence factors - they are critical to or assist in survival of UPEC as pathogens, and the structure and/or function may be specialized for the infection. Other traits could serve as anti-virulence factors - they represent liability in the urinary tract and are under selection to be lost or inactivated during the infection. Inactivation, variation, or other changes of the bacterial genes that increase the pathogen's fitness during the infection are called pathoadaptive mutations. This chapter describes examples of pathoadaptive mutations in UPEC and provides rationale for their further in-depth study.

Cj0440c Affects Flagella Formation and In Vivo Colonization of Erythromycin-Susceptible and -Resistant Campylobacter jejuni

Campylobacter jejuni is one of the most common foodborne pathogen worldwide. A putative transcriptional regulator, Cj0440c, was up-regulated in the erythromycin-resistant C. jejuni, however, the precise role of Cj0440c is yet to be determined. The aim of this study was to determine the biological functions of Cj0440c. The Cj0440c isogenic mutants were constructed from erythromycin-susceptible C. jejuni NCTC 11168 (S) and -resistant C. jejuni 68-ER (R), designating as SM and RM, respectively. The isogenic Cj0440c mutants (SM and RM) and parental strains (S and R) were subjected to microarray and qRT-PCR analysis to examine the transcriptional profile changes contributed by Cj0440c. The antimicrobial susceptibility, flagellar morphology, in vitro growth and in vivo colonization in chickens were carried out to analyze the biological function of Cj0440c. The results showed that 17 genes were down-regulated in SM compared to S, while 9 genes were down-regulated in RM compared to R. The genes with transcriptional change were mainly involved in flagella biosynthesis and assembly. Using transmission electron microscopy, we found that the filaments were impaired in SM and lost in RM. The chicken colonization experiments showed that Cj0440c mutants (SM and RM) had reduced colonization ability in chickens when compared with corresponding parental strains (S and R). In conclusion, Cj0440c regulates flagella biosynthesis and assembly, and consequently affect the in vivo colonization of erythromycin-susceptible and -resistant C. jejuni.

Competing for Iron: Duplication and Amplification of the isd Locus in Staphylococcus lugdunensis HKU09-01 Provides a Competitive Advantage to Overcome Nutritional Limitation

Staphylococcus lugdunensis is a coagulase negative bacterial pathogen that is particularly associated with severe cases of infectious endocarditis. Unique amongst the coagulase-negative staphylococci, S. lugdunensis harbors an iron regulated surface determinant locus (isd). This locus facilitates the acquisition of heme as a source of nutrient iron during infection and allows iron limitation caused by “nutritional immunity” to be overcome. The isd locus is duplicated in S. lugdunensis HKU09-01 and we show here that the duplication is intrinsically unstable and undergoes accordion-like amplification and segregation leading to extensive isd copy number variation. Amplification of the locus increased the level of expression of Isd proteins and improved binding of hemoglobin to the cell surface of S. lugdunensis. Furthermore, Isd overexpression provided an advantage when strains were competing for a limited amount of hemoglobin as the sole source of iron. Gene duplications and amplifications (GDA) are events of fundamental importance for bacterial evolution and are frequently associated with antibiotic resistance in many species. As such, GDAs are regarded as evolutionary adaptions to novel selective pressures in hostile environments pointing towards a special importance of isd for S. lugdunensis. For the first time we show an example of a GDA that involves a virulence factor of a Gram-positive pathogen and link the GDA directly to a competitive advantage when the bacteria were struggling with selective pressures mimicking “nutritional immunity”.

Sometimes changing environmental conditions force bacteria to boost protein expression above the level that can be achieved by transcriptional or translational control. Gene duplication and amplification (GDA) represents a simple and effective means to augment protein expression and is frequently associated with the development of resistance against antibacterial agents. As such GDAs can be seen as recent evolutionary adaptions towards novel selective pressures within the environment. We describe herein the GDA of the isd locus in a clinical isolate of Staphylococcus lugdunensis a bacterial species frequently associated with infectious endocarditis. This amplification made individual lineages significantly more successful when competing for the trace element iron supplied in the form of hemoglobin (hb). Human body fluids are actively depleted of iron to prevent bacterial proliferation and hb represents the most abundant iron source during infection. This deepens our understanding about adaption processes of pathogens and demonstrates how selective pressures drive the evolution of pathogens to become more successful. This study suggests that investigating GDAs in clinical isolates can help to identify chromosomal regions that are of special importance for bacterial fitness under changing environmental conditions. The gene products can therefore represent attractive targets for novel antibiotics supporting the natural immune defenses by targeting the vulnerable spots of the bacterial anti-immunity strategies.

Host Mitochondrial Association Evolved in the Human Parasite Toxoplasma gondii via Neofunctionalization of a Gene Duplicate

In Toxoplasma gondii, an intracellular parasite of humans and other animals, host mitochondrial association (HMA) is driven by a gene family that encodes multiple mitochondrial association factor 1 (MAF1) proteins. However, the importance of MAF1 gene duplication in the evolution of HMA is not understood, nor is the impact of HMA on parasite biology. Here we used within- and between-species comparative analysis to determine that the MAF1 locus is duplicated in T. gondii and its nearest extant relative Hammondia hammondi, but not another close relative, Neospora caninum. Using cross-species complementation, we determined that the MAF1 locus harbors multiple distinct paralogs that differ in their ability to mediate HMA, and that only T. gondii and H. hammondi harbor HMA⁺ paralogs. Additionally, we found that exogenous expression of an HMA⁺ paralog in T. gondii strains that do not normally exhibit HMA provides a competitive advantage over their wild-type counterparts during a mouse infection. These data indicate that HMA likely evolved by neofunctionalization of a duplicate MAF1 copy in the common ancestor of T. gondii and H. hammondi, and that the neofunctionalized gene duplicate is selectively advantageous.

Reinterpreting Long-Term Evolution Experiments: Is Delayed Adaptation an Example of Historical Contingency or a Consequence of Intermittent Selection?

Van Hofwegen et al. demonstrated that Escherichia coli rapidly evolves the ability to use citrate when long selective periods are provided (D. J. Van Hofwegen, C. J. Hovde, and S. A. Minnich, J Bacteriol 198:1022-1034, 2016, http://dx.doi.org/10.1128/JB.00831-15). This contrasts with the extreme delay (15 years of daily transfers) seen in the long-term evolution experiments of Lenski and coworkers. Their idea of "historical contingency" may require reinterpretation. Rapid evolution seems to involve selection for duplications of the whole cit locus that are too unstable to contribute when selection is provided in short pulses.

Rapid Evolution of Citrate Utilization by Escherichia coli by Direct Selection Requires citT and dctA

The isolation of aerobic citrate-utilizing Escherichia coli (Cit(+)) in long-term evolution experiments (LTEE) has been termed a rare, innovative, presumptive speciation event. We hypothesized that direct selection would rapidly yield the same class of E. coli Cit(+) mutants and follow the same genetic trajectory: potentiation, actualization, and refinement. This hypothesis was tested with wild-type E. coli strain B and with K-12 and three K-12 derivatives: an E. coli ΔrpoS::kan mutant (impaired for stationary-phase survival), an E. coli ΔcitT::kan mutant (deleted for the anaerobic citrate/succinate antiporter), and an E. coli ΔdctA::kan mutant (deleted for the aerobic succinate transporter). E. coli underwent adaptation to aerobic citrate metabolism that was readily and repeatedly achieved using minimal medium supplemented with citrate (M9C), M9C with 0.005% glycerol, or M9C with 0.0025% glucose. Forty-six independent E. coli Cit(+) mutants were isolated from all E. coli derivatives except the E. coli ΔcitT::kan mutant. Potentiation/actualization mutations occurred within as few as 12 generations, and refinement mutations occurred within 100 generations. Citrate utilization was confirmed using Simmons, Christensen, and LeMaster Richards citrate media and quantified by mass spectrometry. E. coli Cit(+) mutants grew in clumps and in long incompletely divided chains, a phenotype that was reversible in rich media. Genomic DNA sequencing of four E. coli Cit(+) mutants revealed the required sequence of mutational events leading to a refined Cit(+) mutant. These events showed amplified citT and dctA loci followed by DNA rearrangements consistent with promoter capture events for citT. These mutations were equivalent to the amplification and promoter capture CitT-activating mutations identified in the LTEE.IMPORTANCE E. coli cannot use citrate aerobically. Long-term evolution experiments (LTEE) performed by Blount et al. (Z. D. Blount, J. E. Barrick, C. J. Davidson, and R. E. Lenski, Nature 489:513-518, 2012, http://dx.doi.org/10.1038/nature11514 ) found a single aerobic, citrate-utilizing E. coli strain after 33,000 generations (15 years). This was interpreted as a speciation event. Here we show why it probably was not a speciation event. Using similar media, 46 independent citrate-utilizing mutants were isolated in as few as 12 to 100 generations. Genomic DNA sequencing revealed an amplification of the citT and dctA loci and DNA rearrangements to capture a promoter to express CitT, aerobically. These are members of the same class of mutations identified by the LTEE. We conclude that the rarity of the LTEE mutant was an artifact of the experimental conditions and not a unique evolutionary event. No new genetic information (novel gene function) evolved.

The Origin of Mutants under Selection: Interactions of Mutation, Growth, and Selection

The classical experiments of Luria and Delbrück showed convincingly that mutations exist before selection and do not contribute to the creation of mutations when selection is lethal. In contrast, when nonlethal selections are used,measuring mutation rates and separating the effects of mutation and selection are difficult and require methods to fully exclude growth after selection has been applied. Although many claims of stress-induced mutagenesis have been made, it is difficult to exclude the influence of growth under nonlethal selection conditions in accounting for the observed increases in mutant frequency. Instead, for many of the studied experimental systems the increase in mutant frequency can be explainedbetter by the ability of selection to detect small differences in growth rate caused by common small effect mutations. A verycommon mutant class,found in response to many different types of selective regimensin which increased gene dosage can resolve the problem, is gene amplification. In the well-studiedlac system of Cairns and Foster, the apparent increase in Lac+revertants can be explained by high-level amplification of the lac operon and the increased probability for a reversion mutation to occur in any one of the amplified copies. The associated increase in general mutation rate observed in revertant cells in that system is an artifact caused by the coincidental co-amplification of the nearby dinB gene (encoding the error-prone DNA polymerase IV) on the particular plasmid used for these experiments. Apart from the lac system, similar gene amplification processes have been described for adaptation to toxic drugs, growth in host cells, and various nutrient limitations.

Experimental Evolution Identifies Vaccinia Virus Mutations in A24R and A35R That Antagonize the Protein Kinase R Pathway and Accompany Collapse of an Extragenic Gene Amplification

Most new human infectious diseases emerge from cross-species pathogen transmissions; however, it is not clear how viruses adapt to productively infect new hosts. Host restriction factors represent one species-specific barrier that viruses may initially have little ability to inhibit in new hosts. For example, viral antagonists of protein kinase R (PKR) vary in their ability to block PKR-mediated inhibition of viral replication, in part due to PKR allelic variation between species. We previously reported that amplification of a weak PKR antagonist encoded by rhesus cytomegalovirus, rhtrs1, improved replication of a recombinant poxvirus (VVΔEΔK+RhTRS1) in several resistant primate cells. To test whether amplification increases the opportunity for mutations that improve virus replication with only a single copy of rhtrs1 to evolve, we passaged rhtrs1-amplified viruses in semipermissive primate cells. After passage, we isolated two viruses that contained only a single copy of rhtrs1 yet replicated as well as the amplified virus. Surprisingly, rhtrs1 was not mutated in these viruses; instead, we identified mutations in two vaccinia virus (VACV) genes, A24R and A35R, either of which was sufficient to improve VVΔEΔK+RhTRS1 replication. Neither of these genes has previously been implicated in PKR antagonism. Furthermore, the mutation in A24R, but not A35R, increased resistance to the antipoxviral drug isatin-β-thiosemicarbazone, suggesting that these mutations employ different mechanisms to evade PKR. This study supports our hypothesis that gene amplification may provide a "molecular foothold," broadly improving replication to facilitate rapid adaptation, while subsequent mutations maintain this efficient replication in the new host without requiring gene amplification.

IMPORTANCE

Understanding how viruses adapt to a new host may help identify viruses poised to cross species barriers before an outbreak occurs. Amplification of rhtrs1, a weak viral antagonist of the host antiviral protein PKR, enabled a recombinant vaccinia virus to replicate in resistant cells from humans and other primates. After serial passage of rhtrs1-amplified viruses, there arose in two vaccinia virus genes mutations that improved viral replication without requiring rhtrs1 amplification. Neither of these genes has previously been associated with inhibition of the PKR pathway. These data suggest that gene amplification can improve viral replication in a resistant host species and facilitate the emergence of novel adaptations that maintain the foothold needed for continued replication and spread in the new host.

The Origin of Mutants Under Selection: How Natural Selection Mimics Mutagenesis (Adaptive Mutation)

Selection detects mutants but does not cause mutations. Contrary to this dictum, Cairns and Foster plated a leaky lac mutant of Escherichia coli on lactose medium and saw revertant (Lac(+)) colonies accumulate with time above a nongrowing lawn. This result suggested that bacteria might mutagenize their own genome when growth is blocked. However, this conclusion is suspect in the light of recent evidence that revertant colonies are initiated by preexisting cells with multiple copies the conjugative F'lac plasmid, which carries the lac mutation. Some plated cells have multiple copies of the simple F'lac plasmid. This provides sufficient LacZ activity to support plasmid replication but not cell division. In nongrowing cells, repeated plasmid replication increases the likelihood of a reversion event. Reversion to lac(+) triggers exponential cell growth leading to a stable Lac(+) revertant colony. In 10% of these plated cells, the high-copy plasmid includes an internal tandem lac duplication, which provides even more LacZ activity—sufficient to support slow growth and formation of an unstable Lac(+) colony. Cells with multiple copies of the F'lac plasmid have an increased mutation rate, because the plasmid encodes the error-prone (mutagenic) DNA polymerase, DinB. Without DinB, unstable and stable Lac(+) revertant types form in equal numbers and both types arise with no mutagenesis. Amplification and selection are central to behavior of the Cairns-Foster system, whereas mutagenesis is a system-specific side effect or artifact caused by coamplification of dinB with lac. Study of this system has revealed several broadly applicable principles. In all populations, gene duplications are frequent stable genetic polymorphisms, common near-neutral mutant alleles can gain a positive phenotype when amplified under selection, and natural selection can operate without cell division when variability is generated by overreplication of local genome subregions.

Mechanisms of gene duplication and amplification

Changes in gene copy number are among the most frequent mutational events in all genomes and were among the mutations for which a physical basis was first known. Yet mechanisms of gene duplication remain uncertain because formation rates are difficult to measure and mechanisms may vary with position in a genome. Duplications are compared here to deletions, which seem formally similar but can arise at very different rates by distinct mechanisms. Methods of assessing duplication rates and dependencies are described with several proposed formation mechanisms. Emphasis is placed on duplications formed in extensively studied experimental situations. Duplications studied in microbes are compared with those observed in metazoan cells, specifically those in genomes of cancer cells. Duplications, and especially their derived amplifications, are suggested to form by multistep processes often under positive selection for increased copy number.

Experimental design to evaluate directed adaptive mutation in Mammalian cells

BACKGROUND

We describe the experimental design for a methodological approach to determine whether directed adaptive mutation occurs in mammalian cells. Identification of directed adaptive mutation would have profound practical significance for a wide variety of biomedical problems, including disease development and resistance to treatment. In adaptive mutation, the genetic or epigenetic change is not random; instead, the presence and type of selection influences the frequency and character of the mutation event. Adaptive mutation can contribute to the evolution of microbial pathogenesis, cancer, and drug resistance, and may become a focus of novel therapeutic interventions.

OBJECTIVE

Our experimental approach was designed to distinguish between 3 types of mutation: (1) random mutations that are independent of selective pressure, (2) undirected adaptive mutations that arise when selective pressure induces a general increase in the mutation rate, and (3) directed adaptive mutations that arise when selective pressure induces targeted mutations that specifically influence the adaptive response. The purpose of this report is to introduce an experimental design and describe limited pilot experiment data (not to describe a complete set of experiments); hence, it is an early report.

METHODS

An experimental design based on immortalization of mouse embryonic fibroblast cells is presented that links clonal cell growth to reversal of an inactivating polyadenylation site mutation. Thus, cells exhibit growth only in the presence of both the countermutation and an inducing agent (doxycycline). The type and frequency of mutation in the presence or absence of doxycycline will be evaluated. Additional experimental approaches would determine whether the cells exhibit a generalized increase in mutation rate and/or whether the cells show altered expression of error-prone DNA polymerases or of mismatch repair proteins.

RESULTS

We performed the initial stages of characterizing our system and have limited preliminary data from several pilot experiments. Cell growth and DNA sequence data indicate that we have identified a cell clone that exhibits several suitable characteristics, although further study is required to identify a more optimal cell clone.

CONCLUSIONS

The experimental approach is based on a quantum biological model of basis-dependent selection describing a novel mechanism of adaptive mutation. This project is currently inactive due to lack of funding. However, consistent with the objective of early reports, we describe a proposed study that has not produced publishable results, but is worthy of report because of the hypothesis, experimental design, and protocols. We outline the project's rationale and experimental design, with its strengths and weaknesses, to stimulate discussion and analysis, and lay the foundation for future studies in this field.

Large chromosomal rearrangements during a long-term evolution experiment with Escherichia coli

Large-scale rearrangements may be important in evolution because they can alter chromosome organization and gene expression in ways not possible through point mutations. In a long-term evolution experiment, twelve Escherichia coli populations have been propagated in a glucose-limited environment for over 25 years. We used whole-genome mapping (optical mapping) combined with genome sequencing and PCR analysis to identify the large-scale chromosomal rearrangements in clones from each population after 40,000 generations. A total of 110 rearrangement events were detected, including 82 deletions, 19 inversions, and 9 duplications, with lineages having between 5 and 20 events. In three populations, successive rearrangements impacted particular regions. In five populations, rearrangements affected over a third of the chromosome. Most rearrangements involved recombination between insertion sequence (IS) elements, illustrating their importance in mediating genome plasticity. Two lines of evidence suggest that at least some of these rearrangements conferred higher fitness. First, parallel changes were observed across the independent populations, with ~65% of the rearrangements affecting the same loci in at least two populations. For example, the ribose-utilization operon and the manB-cpsG region were deleted in 12 and 10 populations, respectively, suggesting positive selection, and this inference was previously confirmed for the former case. Second, optical maps from clones sampled over time from one population showed that most rearrangements occurred early in the experiment, when fitness was increasing most rapidly. However, some rearrangements likely occur at high frequency and may have simply hitchhiked to fixation. In any case, large-scale rearrangements clearly influenced genomic evolution in these populations.

Bacterial chromosomes are dynamic structures shaped by long histories of evolution. Among genomic changes, large-scale DNA rearrangements can have important effects on the presence, order, and expression of genes. Whole-genome sequencing that relies on short DNA reads cannot identify all large-scale rearrangements. Therefore, deciphering changes in the overall organization of genomes requires alternative methods, such as optical mapping. We analyzed the longest-running microbial evolution experiment (more than 25 years of evolution in the laboratory) by optical mapping, genome sequencing, and PCR analyses. We found multiple large genome rearrangements in all 12 independently evolving populations. In most cases, it is unclear whether these changes were beneficial themselves or, alternatively, hitchhiked to fixation with other beneficial mutations. In any case, many genome rearrangements accumulated over decades of evolution, providing these populations with genetic plasticity reminiscent of that observed in some pathogenic bacteria.

Plasmid copy number underlies adaptive mutability in bacteria

The origin of mutations under selection has been intensively studied using the Cairns-Foster system, in which cells of an Escherichia coli lac mutant are plated on lactose and give rise to 100 Lac+ revertants over several days. These revertants have been attributed variously to stress-induced mutagenesis of nongrowing cells or to selective improvement of preexisting weakly Lac+ cells with no mutagenesis. Most revertant colonies (90%) contain stably Lac+ cells, while others (10%) contain cells with an unstable amplification of the leaky mutant lac allele. Evidence is presented that both stable and unstable Lac+ revertant colonies are initiated by preexisting cells with multiple copies of the F'lac plasmid, which carries the mutant lac allele. The tetracycline analog anhydrotetracycline (AnTc) inhibits growth of cells with multiple copies of the tetA gene. Populations with tetA on their F'lac plasmid include rare cells with an elevated plasmid copy number and multiple copies of both the tetA and lac genes. Pregrowth of such populations with AnTc reduces the number of cells with multiple F'lac copies and consequently the number of Lac+ colonies appearing under selection. Revertant yield is restored rapidly by a few generations of growth without AnTc. We suggest that preexisting cells with multiple F'lac copies divide very little under selection but have enough energy to replicate their F'lac plasmids repeatedly until reversion initiates a stable Lac+ colony. Preexisting cells whose high-copy plasmid includes an internal lac duplication grow under selection and produce an unstable Lac+ colony. In this model, all revertant colonies are initiated by preexisting cells and cannot be stress induced.

Copy-number changes in evolution: rates, fitness effects and adaptive significance

Gene copy-number differences due to gene duplications and deletions are rampant in natural populations and play a crucial role in the evolution of genome complexity. Per-locus analyses of gene duplication rates in the pre-genomic era revealed that gene duplication rates are much higher than the per nucleotide substitution rate. Analyses of gene duplication and deletion rates in mutation accumulation lines of model organisms have revealed that these high rates of copy-number mutations occur at a genome-wide scale. Furthermore, comparisons of the spontaneous duplication and deletion rates to copy-number polymorphism data and bioinformatic-based estimates of duplication rates from sequenced genomes suggest that the vast majority of gene duplications are detrimental and removed by natural selection. The rate at which new gene copies appear in populations greatly influences their evolutionary dynamics and standing gene copy-number variation in populations. The opportunity for mutations that result in the maintenance of duplicate copies, either through neofunctionalization or subfunctionalization, also depends on the equilibrium frequency of additional gene copies in the population, and hence on the spontaneous gene duplication (and loss) rate. The duplication rate may therefore have profound effects on the role of adaptation in the evolution of duplicated genes as well as important consequences for the evolutionary potential of organisms. We further discuss the broad ramifications of this standing gene copy-number variation on fitness and adaptive potential from a population-genetic and genome-wide perspective.

Recombination and annealing pathways compete for substrates in making rrn duplications in Salmonella enterica

Tandem genetic duplications arise frequently between the seven directly repeated 5.5-kb rrn loci that encode ribosomal RNAs in Salmonella enterica. The closest rrn genes, rrnB and rrnE, flank a 40-kb region that includes the purHD operon. Duplications of purHD arise by exchanges between rrn loci and form at a high rate (10(-3)/cell/division) that remains high in strains blocked for early steps in recombination (recA, recB, and/or recF), but drops 30-fold in mutants blocked for later Holliday junction resolution (ruvC recG). The duplication defect of a ruvC recG mutant was fully corrected by an added mutation in any one of the recA, recB, or recF genes. To explain these results, we propose that early recombination defects activate an alternative single-strand annealing pathway for duplication formation. In wild-type cells, rrn duplications form primarily by the action of RecFORA on single-strand gaps. Double-strand breaks cannot initiate rrn duplications because rrn loci lack Chi sites, which are essential for recombination between two separated rrn sequences. A recA or recF mutation allows unrepaired gaps to accumulate such that different rrn loci can provide single-strand rrn sequences that lack the RecA coating that normally inhibits annealing. A recB mutation activates annealing by allowing double-strand ends within rrn to avoid digestion by RecBCD and provide a new source of rrn ends for use in annealing. The equivalent high rates of rrn duplication by recombination and annealing pathways may reflect a limiting economy of gaps and breaks arising in heavily transcribed, palindrome-rich rrn sequences.

A case of adaptation through a mutation in a tandem duplication during experimental evolution in Escherichia coli

Background

DNA duplications constitute important precursors for genome variation. Here we analyzed an unequal duplication harboring a beneficial mutation that may provide alternative evolutionary outcomes.

Results

We characterized this evolutionary event during experimental evolution for only 100 generations of an Escherichia coli strain under glucose limitation within chemostats. By combining Insertion Sequence based Restriction Length Polymorphism experiments, pulsed field gel electrophoresis and two independent genome re-sequencing experiments, we identified an evolved lineage carrying a 180 kb duplication of the 46’ region of the E. coli chromosome. This evolved duplication revealed a heterozygous state, with one copy harboring a 2668 bp deletion that included part of the ogrK gene and both the yegR and yegS genes. By genetically manipulating ancestral and evolved strains, we showed that the single yegS inactivation was sufficient to confer a frequency dependent fitness increase under the chemostat selective conditions in both the ancestor and evolved genetic contexts, implying that the duplication itself was not a direct fitness contributor. Nonetheless, the heterozygous duplicated state was relatively stable in the conditions prevailing during evolution in chemostats, in striking contrast to non selective conditions in which the duplication resolved at high frequency into either its ancestral or deleted copy.

Conclusions

Our results suggest that the duplication state may constitute a second order selection process providing higher evolutionary potential. Moreover, its heterozygous nature may provide differential evolutionary opportunities in alternating environments. Our results also highlighted how careful analyses of whole genome data are needed to identify such complex rearrangements.

Quantum biology at the cellular level--elements of the research program

Quantum biology is emerging as a new field at the intersection between fundamental physics and biology, promising novel insights into the nature and origin of biological order. We discuss several elements of QBCL (quantum biology at cellular level) - a research program designed to extend the reach of quantum concepts to higher than molecular levels of biological organization. We propose a new general way to address the issue of environmentally induced decoherence and macroscopic superpositions in biological systems, emphasizing the 'basis-dependent' nature of these concepts. We introduce the notion of 'formal superposition' and distinguish it from that of Schroedinger's cat (i.e., a superposition of macroscopically distinct states). Whereas the latter notion presents a genuine foundational problem, the former one contradicts neither common sense nor observation, and may be used to describe cellular 'decision-making' and adaptation. We stress that the interpretation of the notion of 'formal superposition' should involve non-classical correlations between molecular events in a cell. Further, we describe how better understanding of the physics of Life can shed new light on the mechanism driving evolutionary adaptation (viz., 'Basis-Dependent Selection', BDS). Experimental tests of BDS and the potential role of synthetic biology in closing the 'evolvability mechanism' loophole are also discussed.

Gross chromosomal rearrangement mediated by DNA replication in stressed cells: evidence from Escherichia coli

Gross chromosomal rearrangements (GCRs), or changes in chromosome structure, play central roles in evolution and are central to cancer formation and progression. GCRs underlie copy number variation (CNV), and therefore genomic disorders that stem from CNV. We study amplification in Escherichia coli as a model system to understand mechanisms and circumstances of GCR formation. Here, we summarize observations that led us to postulate that GCR occurs by a replicative mechanism as part of activated stress responses. We report that we do not find RecA to be downregulated by stress on a population basis and that constitutive expression of RecA does not inhibit amplification, as would be expected if downregulation of RecA made cells permissive for nonhomologous recombination. Strains deleted for the genes for three proteins that inhibit RecA activity, psiB, dinI, and recX, all show unaltered amplification, suggesting that if they do downregulate RecA indirectly, this activity does not promote amplification.

Expression of a cryptic secondary sigma factor gene unveils natural competence for DNA transformation in Staphylococcus aureus

It has long been a question whether Staphylococcus aureus, a major human pathogen, is able to develop natural competence for transformation by DNA. We previously showed that a novel staphylococcal secondary sigma factor, SigH, was a likely key component for competence development, but the corresponding gene appeared to be cryptic as its expression could not be detected during growth under standard laboratory conditions. Here, we have uncovered two distinct mechanisms allowing activation of SigH production in a minor fraction of the bacterial cell population. The first is a chromosomal gene duplication rearrangement occurring spontaneously at a low frequency [≤10⁻⁵], generating expression of a new chimeric sigH gene. The second involves post-transcriptional regulation through an upstream inverted repeat sequence, effectively suppressing expression of the sigH gene. Importantly, we have demonstrated for the first time that S. aureus cells producing active SigH become competent for transformation by plasmid or chromosomal DNA, which requires the expression of SigH-controlled competence genes. Additionally, using DNA from the N315 MRSA strain, we successfully transferred the full length SCCmecII element through natural transformation to a methicillin-sensitive strain, conferring methicillin resistance to the resulting S. aureus transformants. Taken together, we propose a unique model for staphylococcal competence regulation by SigH that could help explain the acquisition of antibiotic resistance genes through horizontal gene transfer in this important pathogen.

Staphylococcus aureus is a major human pathogen responsible for a broad spectrum of infections, emphasized by the emergence of multiple antibiotic-resistant strains with up to 60% of strains worldwide resistant to methicillin (Methicillin Resistant Staphylococcus aureus or MRSA). Indeed, MRSA-related infections are now one of the leading causes of death in the USA, highlighting the growing threat this bacterium poses to human health. Many bacteria have the ability to acquire novel genetic characteristics, including antibiotic resistance, through the uptake of extracellular DNA, a phenomenon known as natural genetic transformation or competence. We have shown that the SigH staphylococcal sigma factor is a likely key component for competence development, but that its gene is not expressed under standard laboratory conditions. Here, we have uncovered two distinct mechanisms allowing activation of SigH production in S. aureus: a chromosomal gene duplication rearrangement and post-transcriptional regulation through an upstream inverted repeat sequence. Importantly, we have demonstrated for the first time that S. aureus cells producing active SigH become competent for natural transformation by plasmid or chromosomal DNA, and we were able to confer methicillin resistance to a methicillin-sensitive strain by transformation with chromosomal DNA. SigH-dependent competence development in S. aureus could help explain the acquisition of antibiotic resistance genes and the rise of the so-called “superbug."

Multiple pathways of duplication formation with and without recombination (RecA) in Salmonella enterica

Duplications are often attributed to "unequal recombination" between separated, directly repeated sequence elements (>100 bp), events that leave a recombinant element at the duplication junction. However, in the bacterial chromosome, duplications form at high rates (10(-3)-10(-5)/cell/division) even without recombination (RecA). Here we describe 1800 spontaneous lac duplications trapped nonselectively on the low-copy F'(128) plasmid, where lac is flanked by direct repeats of the transposable element IS3 (1258 bp) and by numerous quasipalindromic REP elements (30 bp). Duplications form at a high rate (10(-4)/cell/division) that is reduced only about 11-fold in the absence of RecA. With and without RecA, most duplications arise by recombination between IS3 elements (97%). Formation of these duplications is stimulated by IS3 transposase (Tnp) and plasmid transfer functions (TraI). Three duplication pathways are proposed. First, plasmid dimers form at a high rate stimulated by RecA and are then modified by deletions between IS3 elements (resolution) that leave a monomeric plasmid with an IS3-flanked lac duplication. Second, without RecA, duplications occur by single-strand annealing of DNA ends generated in different sister chromosomes after transposase nicks DNA near participating IS3 elements. The absence of RecA may stimulate annealing by allowing chromosome breaks to persist. Third, a minority of lac duplications (3%) have short (0-36 bp) junction sequences (SJ), some of which are located within REP elements. These duplication types form without RecA, Tnp, or Tra by a pathway in which the palindromic junctions of a tandem inversion duplication (TID) may stimulate deletions that leave the final duplication.

Analysis of IS1236-mediated gene amplification events in Acinetobacter baylyi ADP1

Recombination between insertion sequence copies can cause genetic deletion, inversion, or duplication. However, it is difficult to assess the fraction of all genomic rearrangements that involve insertion sequences. In previous gene duplication and amplification studies of Acinetobacter baylyi ADP1, an insertion sequence was evident in approximately 2% of the characterized duplication sites. Gene amplification occurs frequently in all organisms and has a significant impact on evolution, adaptation, drug resistance, cancer, and various disorders. To understand the molecular details of this important process, a previously developed system was used to analyze gene amplification in selected mutants. The current study focused on amplification events in two chromosomal regions that are near one of six copies of the only transposable element in ADP1, IS1236 (an IS3 family member). Twenty-one independent mutants were analyzed, and in contrast to previous studies of a different chromosomal region, IS1236 was involved in 86% of these events. IS1236-mediated amplification could occur through homologous recombination between insertion sequences on both sides of a duplicated region. However, this mechanism presupposes that transposition generates an appropriately positioned additional copy of IS1236. To evaluate this possibility, PCR and Southern hybridization were used to determine the chromosomal configurations of amplification mutants involving IS1236. Surprisingly, the genomic patterns were inconsistent with the hypothesis that intramolecular homologous recombination occurred between insertion sequences following an initial transposition event. These results raise a novel possibility that the gene amplification events near the IS1236 elements arise from illegitimate recombination involving transposase-mediated DNA cleavage.

What limits the efficiency of double-strand break-dependent stress-induced mutation in Escherichia coli?

Stress-induced mutation is a collection of molecular mechanisms in bacterial, yeast and human cells that promote mutagenesis specifically when cells are maladapted to their environment, i.e. when they are stressed. Here, we review one molecular mechanism: double-strand break (DSB)-dependent stress-induced mutagenesis described in starving Escherichia coli. In it, the otherwise high-fidelity process of DSB repair by homologous recombination is switched to an error-prone mode under the control of the RpoS general stress response, which licenses the use of error-prone DNA polymerase, DinB, in DSB repair. This mechanism requires DSB repair proteins, RpoS, the SOS response and DinB. This pathway underlies half of spontaneous chromosomal frameshift and base substitution mutations in starving E. coli [Proc Natl Acad Sci USA 2011;108:13659-13664], yet appeared less efficient in chromosomal than F' plasmid-borne genes. Here, we demonstrate and quantify DSB-dependent stress-induced reversion of a chromosomal lac allele with DSBs supplied by I-SceI double-strand endonuclease. I-SceI-induced reversion of this allele was previously studied in an F'. We compare the efficiencies of mutagenesis in the two locations. When we account for contributions of an F'-borne extra dinB gene, strain background differences, and bypass considerations of rates of spontaneous DNA breakage by providing I-SceI cuts, the chromosome is still ∼100 times less active than F. We suggest that availability of a homologous partner molecule for recombinational break repair may be limiting. That partner could be a duplicated chromosomal segment or sister chromosome.

Ecological and temporal constraints in the evolution of bacterial genomes

Studies on the experimental evolution of microorganisms, on their in vivo evolution (mainly in the case of bacteria producing chronic infections), as well as the availability of multiple full genomic sequences, are placing bacteria in the playground of evolutionary studies. In the present article we review the differential contribution to the evolution of bacterial genomes that processes such as gene modification, gene acquisition and gene loss may have when bacteria colonize different habitats that present characteristic ecological features. In particular, we review how the different processes contribute to evolution in microbial communities, in free-living bacteria or in bacteria living in isolation. In addition, we discuss the temporal constraints in the evolution of bacterial genomes, considering bacterial evolution from the perspective of processes of short-sighted evolution and punctual acquisition of evolutionary novelties followed by long stasis periods.

Global chromosomal structural instability in a subpopulation of starving Escherichia coli cells

Copy-number variations (CNVs) constitute very common differences between individual humans and possibly all genomes and may therefore be important fuel for evolution, yet how they form remains elusive. In starving Escherichia coli, gene amplification is induced by stress, controlled by the general stress response. Amplification has been detected only encompassing genes that confer a growth advantage when amplified. We studied the structure of stress-induced gene amplification in starving cells in the Lac assay in Escherichia coli by array comparative genomic hybridization (aCGH), with polymerase chain reaction (pcr) and DNA sequencing to establish the structures generated. About 10% of 300 amplified isolates carried other chromosomal structural change in addition to amplification. Most of these were inversions and duplications associated with the amplification event. This complexity supports a mechanism similar to that seen in human non-recurrent copy number variants. We interpret these complex events in terms of repeated template switching during DNA replication. Importantly, we found a significant occurrence (6 out of 300) of chromosomal structural changes that were apparently not involved in the amplification event. These secondary changes were absent from 240 samples derived from starved cells not carrying amplification, suggesting that amplification happens in a differentiated subpopulation of stressed cells licensed for global chromosomal structural change and genomic instability. These data imply that chromosomal structural changes occur in bursts or showers of instability that may have the potential to drive rapid evolution.

Much of the difference between individual humans is in the number of copies of genes and lengths of genome. The mechanisms by which copy number variation arises are not well understood. We sought information on copy number change mechanisms by extensive use of array comparative genomic hybridization of whole genomes in bacteria selected for amplification of part of the genome. We report that about 10% of amplified isolates carried other chromosomal structural changes associated with the amplification, a result comparable to that seen in human copy number variants. Importantly, we found a significant occurrence of structural changes that were not involved in the amplification event. These were not seen in a control sample of stressed cells not carrying amplification. This establishes that chromosomal structural change happens in a subpopulation of cells apparently licensed to undergo these changes. Because the changes occur under the stress of starvation and require two of the cells' stress-response systems, we propose that licensing for cell-wide structural change in this subpopulation is a component of response to stress. This idea has implications for the mechanisms of evolution and cancer progression, suggesting that changes occur in a shower of events rather than as isolated random events.

Effect of growth under selection on appearance of chromosomal mutations in Salmonella enterica

Populations adapt physiologically using regulatory mechanisms and genetically by means of mutations that improve growth. During growth under selection, genetic adaptation can be rapid. In several genetic systems, the speed of adaptation has been attributed to cellular mechanisms that increase mutation rates in response to growth limitation. An alternative possibility is that growth limitation serves only as a selective agent but acts on small-effect mutations that are common under all growth conditions. The genetic systems that initially suggested stress-induced mutagenesis have been analyzed without regard for multistep adaptation and some include features that make such analysis difficult. To test the selection-only model, a simpler system is examined, whose behavior was originally attributed to stress-induced mutagenesis (Yang et al. 2001, 2006). A population with a silent chromosomal lac operon gives rise to Lac+ revertant colonies that accumulate over 6 days under selection. Each colony contains a mixture of singly and doubly mutant cells. Evidence is provided that the colonies are initiated by pre-existing single mutants with a weak Lac+ phenotype. Under selection, these cells initiate slow-growing clones, in which a second mutation arises and improves growth of the resulting double mutant. The system shows no evidence of general mutagenesis during selection. Selection alone may explain rapid adaptation in this and other systems that give the appearance of mutagenesis.

Genome-wide detection of chromosomal rearrangements, indels, and mutations in circular chromosomes by short read sequencing

Whole-genome sequencing (WGS) with new short-read sequencing technologies has recently been applied for genome-wide identification of mutations. Genomic rearrangements have, however, often remained undetected by WGS, and additional analyses are required for their detection. Here, we have applied a combination of WGS and genome copy number analysis, for the identification of mutations that suppress the growth deficiency imposed by excessive initiations from the Escherichia coli origin of replication, oriC. The E. coli chromosome, like the majority of bacterial chromosomes, is circular, and DNA replication is initiated by assembling two replication complexes at the origin, oriC. These complexes then replicate the chromosome bidirectionally toward the terminus, ter. In a population of growing cells, this results in a copy number gradient, so that origin-proximal sequences are more frequent than origin-distal sequences. Major rearrangements in the chromosome are, therefore, readily identified by changes in copy number, i.e., certain sequences become over- or under-represented. Of the eight mutations analyzed in detail here, six were found to affect a single gene only, one was a large chromosomal inversion, and one was a large chromosomal duplication. The latter two mutations could not be detected solely by WGS, validating the present approach for identification of genomic rearrangements. We further suggest the use of copy number analysis in combination with WGS for validation of newly assembled bacterial chromosomes.

Horizontal transfer, not duplication, drives the expansion of protein families in prokaryotes

Gene duplication followed by neo- or sub-functionalization deeply impacts the evolution of protein families and is regarded as the main source of adaptive functional novelty in eukaryotes. While there is ample evidence of adaptive gene duplication in prokaryotes, it is not clear whether duplication outweighs the contribution of horizontal gene transfer in the expansion of protein families. We analyzed closely related prokaryote strains or species with small genomes (Helicobacter, Neisseria, Streptococcus, Sulfolobus), average-sized genomes (Bacillus, Enterobacteriaceae), and large genomes (Pseudomonas, Bradyrhizobiaceae) to untangle the effects of duplication and horizontal transfer. After removing the effects of transposable elements and phages, we show that the vast majority of expansions of protein families are due to transfer, even among large genomes. Transferred genes--xenologs--persist longer in prokaryotic lineages possibly due to a higher/longer adaptive role. On the other hand, duplicated genes--paralogs--are expressed more, and, when persistent, they evolve slower. This suggests that gene transfer and gene duplication have very different roles in shaping the evolution of biological systems: transfer allows the acquisition of new functions and duplication leads to higher gene dosage. Accordingly, we show that paralogs share most protein-protein interactions and genetic regulators, whereas xenologs share very few of them. Prokaryotes invented most of life's biochemical diversity. Therefore, the study of the evolution of biology systems should explicitly account for the predominant role of horizontal gene transfer in the diversification of protein families.

Evolvability and Speed of Evolutionary Algorithms in Light of Recent Developments in Biology

Biological and artificial evolutionary systems exhibit varying degrees of evolvability and different rates of evolution. Such quantities can be affected by various factors. Here, we review some evolutionary mechanisms and discuss new developments in biology that can potentially improve evolvability or accelerate evolution in artificial systems. Biological notions are discussed to the degree they correspond to notions in Evolutionary Computation. We hope that the findings put forward here can be used to design computational models of evolution that produce significant gains in evolvability and evolutionary speed.

Metagenomic sequencing of an in vitro-simulated microbial community

Background

Microbial life dominates the earth, but many species are difficult or even impossible to study under laboratory conditions. Sequencing DNA directly from the environment, a technique commonly referred to as metagenomics, is an important tool for cataloging microbial life. This culture-independent approach involves collecting samples that include microbes in them, extracting DNA from the samples, and sequencing the DNA. A sample may contain many different microorganisms, macroorganisms, and even free-floating environmental DNA. A fundamental challenge in metagenomics has been estimating the abundance of organisms in a sample based on the frequency with which the organism's DNA was observed in reads generated via DNA sequencing.

Methodology/Principal Findings

We created mixtures of ten microbial species for which genome sequences are known. Each mixture contained an equal number of cells of each species. We then extracted DNA from the mixtures, sequenced the DNA, and measured the frequency with which genomic regions from each organism was observed in the sequenced DNA. We found that the observed frequency of reads mapping to each organism did not reflect the equal numbers of cells that were known to be included in each mixture. The relative organism abundances varied significantly depending on the DNA extraction and sequencing protocol utilized.

Conclusions/Significance

We describe a new data resource for measuring the accuracy of metagenomic binning methods, created by in vitro-simulation of a metagenomic community. Our in vitro simulation can be used to complement previous in silico benchmark studies. In constructing a synthetic community and sequencing its metagenome, we encountered several sources of observation bias that likely affect most metagenomic experiments to date and present challenges for comparative metagenomic studies. DNA preparation methods have a particularly profound effect in our study, implying that samples prepared with different protocols are not suitable for comparative metagenomics.

The tandem inversion duplication in Salmonella enterica: selection drives unstable precursors to final mutation types

During growth under selection, mutant types appear that are rare in unselected populations. Stress-induced mechanisms may cause these structures or selection may favor a series of standard events that modify common preexisting structures. One such mutation is the short junction (SJ) duplication with long repeats separated by short sequence elements: AB*(CD)*(CD)*E (* = a few bases). Another mutation type, described here, is the tandem inversion duplication (TID), where two copies of a parent sequence flank an inverse-order segment: AB(CD)(E'D'C'B')(CD)E. Both duplication types can amplify by unequal exchanges between direct repeats (CD), and both are rare in unselected cultures but common after prolonged selection for amplification. The observed TID junctions are asymmetric (aTIDs) and may arise from a symmetrical precursor (sTID)-ABCDE(E'D'C'B'A')ABCDE-when sequential deletions remove each palindromic junction. Alternatively, one deletion can remove both sTID junctions to generate an SJ duplication. It is proposed that sTID structures form frequently under all growth conditions, but are usually lost due to their instability and fitness cost. Selection for increased copy number helps retain the sTID and favors deletions that remodel junctions, improve fitness, and allow higher amplification. Growth improves with each step in formation of an SJ or aTID amplification, allowing selection to favor completion of the mutation process.

Alterations in DNA replication and histone levels promote histone gene amplification in Saccharomyces cerevisiae

Gene amplification, a process that increases the copy number of a gene or a genomic region to two or more, is utilized by many organisms in response to environmental stress or decreased levels of a gene product. Our previous studies in Saccharomyces cerevisiae identified the amplification of a histone H2A-H2B gene pair, HTA2-HTB2, in response to the deletion of the other H2A-H2B gene pair, HTA1-HTB1. This amplification arises from a recombination event between two flanking Ty1 elements to form a new, stable circular chromosome and occurs at a frequency higher than has been observed for other Ty1-Ty1 recombination events. To understand the regulation of this amplification event, we screened the S. cerevisiae nonessential deletion set for mutations that alter the amplification frequency. Among the deletions that increase HTA2-HTB2 amplification frequency, we identified those that either decrease DNA replication fork progression (rrm3Delta, dpb3Delta, dpb4Delta, and clb5Delta) or that reduce histone H3-H4 levels (hht2-hhf2Delta). These two classes are related because reduced histone H3-H4 levels increase replication fork pauses, and impaired replication forks cause a reduction in histone levels. Consistent with our mutant screen, we found that the introduction of DNA replication stress by hydroxyurea induces the HTA2-HTB2 amplification event. Taken together, our results suggest that either reduced histone levels or slowed replication forks stimulate the HTA2-HTB2 amplification event, contributing to the restoration of normal chromatin structure.

Duplication frequency in a population of Salmonella enterica rapidly approaches steady state with or without recombination

Tandem duplications are among the most common mutation events. The high loss rate of duplication suggested that the frequency of duplications in a bacterial population (1/1000) might reflect a steady state dictated by relative rates of formation (k(F)) and loss (k(L)). This possibility was tested for three genetic loci. Without homologous recombination (RecA), duplication loss rate dropped essentially to zero, but formation rate decreased only slightly and a steady state was still reached rapidly. Under all conditions, steady state was reached faster than predicted by formation and loss rates alone. A major factor in determining steady state proved to be the fitness cost, which can exceed 40% for some genomic regions. Depending on the region tested, duplications reached 40-98% of the steady-state frequency within 30 generations-approximately the growth required for a single cell to produce a saturated overnight culture or form a large colony on solid medium (10(9) cells). Long-term bacterial populations are stably polymorphic for duplications of every region of their genome. These polymorphisms contribute to rapid genetic adaptation by providing frequent preexisting mutations that are beneficial whenever imposed selection favors increases in some gene activity. While the reported results were obtained with the bacterium Salmonella enterica, the genetic implications seem likely to be of broader biological relevance.

The transcription elongation factor NusA is required for stress-induced mutagenesis in Escherichia coli

Stress-induced mutagenesis describes the accumulation of mutations that occur in nongrowing cells, in contrast to mutagenesis that occurs in actively dividing populations, and has been referred to as stationary-phase or adaptive mutagenesis. The most widely studied system for stress-induced mutagenesis involves monitoring the appearance of Lac(+) revertants of the strain FC40 under starvation conditions in Escherichia coli. The SOS-inducible translesion DNA polymerase DinB plays an important role in this phenomenon. Loss of DinB (DNA pol IV) function results in a severe reduction of Lac(+) revertants. We previously reported that NusA, an essential component of elongating RNA polymerases, interacts with DinB. Here we report our unexpected observation that wild-type NusA function is required for stress-induced mutagenesis. We present evidence that this effect is unlikely to be due to defects in transcription of lac genes but rather is due to an inability to adapt and mutate in response to environmental stress. Furthermore, we extended our analysis to the formation of stress-induced mutants in response to antibiotic treatment, observing the same striking abolition of mutagenesis under entirely different conditions. Our results are the first to implicate NusA as a crucial participant in the phenomenon of stress-induced mutagenesis.