Selecting among three basic fitness landscape models: Additive, multiplicative and stickbreaking.

Fitness landscapes map genotypes to organismal fitness. Their topographies depend on how mutational effects interact - epistasis - andare important for understanding evolutionary processes such as speciation, the rate of adaptation, the advantage of recombination, and the predictability versus stochasticity of evolution. The growing amount of data has made it possible to better test landscape models empirically. We argue that this endeavor will benefit from the development and use of meaningful basic models against which to compare more complex models. Here we develop statistical and computational methods for fitting fitness data from mutation combinatorial networks to three simple models: additive, multiplicative and stickbreaking. We employ a Bayesian framework for doing model selection. Using simulations, we demonstrate that our methods work and we explore their statistical performance: bias, error, and the power to discriminate among models. We then illustrate our approach and its flexibility by analyzing several previously published datasets. An R-package that implements our methods is available in the CRAN repository under the name Stickbreaker.

Positive selection and multiple losses of the LINE-1-derived L1TD1 gene in mammals suggest a dual role in genome defense and pluripotency.

Mammalian genomes comprise many active and fossilized retroelements. The obligate requirement for retroelement integration affords host genomes an opportunity to 'domesticate' retroelement genes for their own purpose, leading to important innovations in genome defense and placentation. While many such exaptations involve retroviruses, the L1TD1 gene is the only known domesticated gene whose protein-coding sequence is almost entirely derived from a LINE-1 (L1) retroelement. Human L1TD1 has been shown to play an important role in pluripotency maintenance. To investigate how this role was acquired, we traced the origin and evolution of L1TD1. We find that L1TD1 originated in the common ancestor of eutherian mammals, but was lost or pseudogenized multiple times during mammalian evolution. We also find that L1TD1 has evolved under positive selection during primate and mouse evolution, and that one prosimian L1TD1 has 'replenished' itself with a more recent L1 ORF1 from the prosimian genome. These data suggest that L1TD1 has been recurrently selected for functional novelty, perhaps for a role in genome defense. L1TD1 loss is associated with L1 extinction in several megabat lineages, but not in sigmodontine rodents. We hypothesize that L1TD1 could have originally evolved for genome defense against L1 elements. Later, L1TD1 may have become incorporated into pluripotency maintenance in some lineages. Our study highlights the role of retroelement gene domestication in fundamental aspects of mammalian biology, and that such domesticated genes can adopt different functions in different lineages.

Integration of molecular cytogenetics, dated molecular phylogeny, and model-based predictions to understand the extreme chromosome reorganization in the Neotropical genus Tonatia (Chiroptera: Phyllostomidae).

BACKGROUND: Defining factors that contributed to the fixation of a high number of underdominant chromosomal rearrangements is a complex task because not only molecular mechanisms must be considered, but also the uniqueness of natural history attributes of each taxon. Ideally, detailed investigation of the chromosome architecture of an organism and related groups, placed within a phylogenetic context, is required. We used multiple approaches to investigate the dynamics of chromosomal evolution in lineages of bats with considerable karyotypic variation, focusing on the different facets contributing to fixation of the exceptional chromosomal changes in Tonatia saurophila. Integration of empirical data with proposed models of chromosome evolution was performed to understand the probable conditions for Tonatia's karyotypic evolution. RESULTS: The trajectory of reorganization of chromosome blocks since the common ancestor of Glossophaginae and Phyllostominae subfamilies suggests that multiple tandem fusions, as well as disruption and fusions of conserved phyllostomid chromosomes were major drivers of karyotypic reshuffling in Tonatia. Considerable variation in the rates of chromosomal evolution between phyllostomid lineages was observed. Thirty-nine unique fusions and fission events reached fixation in Tonatia over a short period of time, followed by ~12 million years of chromosomal stasis. Physical mapping of repetitive DNA revealed an unusual accumulation of LINE-1 sequences on centromeric regions, probably associated with the chromosomal dynamics of this genus. CONCLUSIONS: Multiple rearrangements have reached fixation in a wave-like fashion in phyllostomid bats. Different biological features of Tonatia support distinct models of rearrangement fixation, and it is unlikely that the fixations were a result of solely stochastic processes in small ancient populations. Increased recombination rates were probably facilitated by expansion of repetitive DNA, reinforced by aspects of taxon reproduction and ecology.

Billions of basepairs of recently expanded, repetitive sequences are eliminated from the somatic genome during copepod development.

BACKGROUND: Chromatin diminution is the programmed deletion of DNA from presomatic cell or nuclear lineages during development, producing single organisms that contain two different nuclear genomes. Phylogenetically diverse taxa undergo chromatin diminution--some ciliates, nematodes, copepods, and vertebrates. In cyclopoid copepods, chromatin diminution occurs in taxa with massively expanded germline genomes; depending on species, germline genome sizes range from 15 - 75 Gb, 12-74 Gb of which are lost from pre-somatic cell lineages at germline--soma differentiation. This is more than an order of magnitude more sequence than is lost from other taxa. To date, the sequences excised from copepods have not been analyzed using large-scale genomic datasets, and the processes underlying germline genomic gigantism in this clade, as well as the functional significance of chromatin diminution, have remained unknown. RESULTS: Here, we used high-throughput genomic sequencing and qPCR to characterize the germline and somatic genomes of Mesocyclops edax, a freshwater cyclopoid copepod with a germline genome of ~15 Gb and a somatic genome of ~3 Gb. We show that most of the excised DNA consists of repetitive sequences that are either 1) verifiable transposable elements (TEs), or 2) non-simple repeats of likely TE origin. Repeat elements in both genomes are skewed towards younger (i.e. less divergent) elements. Excised DNA is a non-random sample of the germline repeat element landscape; younger elements, and high frequency DNA transposons and LINEs, are disproportionately eliminated from the somatic genome. CONCLUSIONS: Our results suggest that germline genome expansion in M. edax reflects explosive repeat element proliferation, and that billions of base pairs of such repeats are deleted from the somatic genome every generation. Thus, we hypothesize that chromatin diminution is a mechanism that controls repeat element load, and that this load can evolve to be divergent between tissue types within single organisms.

Epistasis between beneficial mutations and the phenotype-to-fitness Map for a ssDNA virus.

Epistatic interactions between genes and individual mutations are major determinants of the evolutionary properties of genetic systems and have therefore been well documented, but few quantitative data exist on epistatic interactions between beneficial mutations, presumably because such mutations are so much rarer than deleterious ones. We explored epistasis for beneficial mutations by constructing genotypes with pairs of mutations that had been previously identified as beneficial to the ssDNA bacteriophage ID11 and by measuring the effects of these mutations alone and in combination. We constructed 18 of the 36 possible double mutants for the nine available beneficial mutations. We found that epistatic interactions between beneficial mutations were all antagonistic-the effects of the double mutations were less than the sums of the effects of their component single mutations. We found a number of cases of decompensatory interactions, an extreme form of antagonistic epistasis in which the second mutation is actually deleterious in the presence of the first. In the vast majority of cases, recombination uniting two beneficial mutations into the same genome would not be favored by selection, as the recombinant could not outcompete its constituent single mutations. In an attempt to understand these results, we developed a simple model in which the phenotypic effects of mutations are completely additive and epistatic interactions arise as a result of the form of the phenotype-to-fitness mapping. We found that a model with an intermediate phenotypic optimum and additive phenotypic effects provided a good explanation for our data and the observed patterns of epistatic interactions.

An empirical test of the mutational landscape model of adaptation using a single-stranded DNA virus.

The primary impediment to formulating a general theory for adaptive evolution has been the unknown distribution of fitness effects for new beneficial mutations. By applying extreme value theory, Gillespie circumvented this issue in his mutational landscape model for the adaptation of DNA sequences, and Orr recently extended Gillespie's model, generating testable predictions regarding the course of adaptive evolution. Here we provide the first empirical examination of this model, using a single-stranded DNA bacteriophage related to phiX174, and find that our data are consistent with Orr's predictions, provided that the model is adjusted to incorporate mutation bias. Orr's work suggests that there may be generalities in adaptive molecular evolution that transcend the biological details of a system, but we show that for the model to be useful as a predictive or inferential tool, some adjustments for the biology of the system will be necessary.

Controlling Recombination to Evolve Bacteriophages.

Recombination among different phages sometimes facilitates their ability to grow on new hosts. Protocols to direct the evolution of phage host range, as might be used in the application of phage therapy, would then benefit from including steps to enable recombination. Applying mathematical and computational models, in addition to experiments using phages T3 and T7, we consider ways that a protocol may influence recombination levels. We first address coinfection, which is the first step to enabling recombination. The multiplicity of infection (MOI, the ratio of phage to cell concentration) is insufficient for predicting (co)infection levels. The force of infection (the rate at which cells are infected) is also critical but is more challenging to measure. Using both a high force of infection and high MOI (>1) for the different phages ensures high levels of coinfection. We also apply a four-genetic-locus model to study protocol effects on recombinant levels. Recombinants accumulate over multiple generations of phage growth, less so if one phage outgrows the other. Supplementing the phage pool with the low-fitness phage recovers some of this 'lost' recombination. Overall, fine tuning of phage recombination rates will not be practical with wild phages, but qualitative enhancement can be attained with some basic procedures.

Genetic and functional basis of the reduction effect in bacteriophage ΦX174.

The ΦX174 reduction effect describes a plasmid-based inhibitory phenomenon that mimics the superinfection inhibition found in wild phage populations. In this effect, when a portion of the ΦX174 genome - the 3' end of the pilot protein gene (H), the 5' end of the replication gene (A), and the H-A intergenic region - is present on a plasmid in the host cell, almost complete protection from phage infection occurs. Here we demonstrate that only the phage pilot protein H portion of the plasmid is sufficient for the observed inhibition, that protein synthesis is necessary for inhibition to occur, that inserting the entire H gene in the plasmid may also impart a blocking effect, and that partial to complete recovery from this inhibition is possible with minimal viral evolution.

Retrofitting the genome: L1 extinction follows endogenous retroviral expansion in a group of muroid rodents.

Long interspersed nuclear element 1 (LINE-1; L1) retrotransposons are the most common retroelements in mammalian genomes. Unlike individual families of endogenous retroviruses (ERVs), they have remained active throughout the mammalian radiation and are responsible for most of the retroelement movement and much genome rearrangement within mammals. They can be viewed as occupying a substantial niche within mammalian genomes. Our previous demonstration that L1s and B1 short interspersed nuclear elements (SINEs) are inactive in a group of South American rodents led us to ask if other elements have amplified to fill the empty niche. We identified a novel and highly active family of ERVs (mysTR). To determine whether loss of L1 activity was correlated with expansion of mysTR, we examined mysTR activity in four South American rodent species that have lost L1 and B1 activity and four sister species with active L1s. The copy number of recent mysTR insertions was extremely high, with an average of 4,200 copies per genome. High copy numbers exist in both L1-active and L1-extinct species, so the mysTR expansion appears to have preceded the loss of both SINE and L1 activity rather than to have filled an empty niche created by their loss. It may be coincidental that two unusual genomic events--loss of L1 activity and massive expansion of an ERV family--occur in the same group of mammals. Alternatively, it is possible that this large ERV expansion set the stage for L1 extinction.

Modeling the Directed Evolution of Broad Host Range Phages.

Background: The host ranges of individual phages tend to be narrow, yet many applications of phages would benefit from expanded host ranges. Empirical methods have been developed to direct the evolution of phages to attack new strains, but the methods have not been evaluated or compared for their consequences. In particular, how do different methods favor generalist (broad host range) phages over specialist phages? All methods involve exposing phages to two or more novel bacterial strains, but the methods differ in the order in which those hosts are presented through time: Parallel presentation, Sequential presentation, and Mixed presentation. Methods: We use a combination of simple analytical methods and numerical analyses to study the effect of these different protocols on the selection of generalist versus specialist phages. Results: The three presentation protocols have profoundly different consequences for the evolution of generalist versus specialist phages. Sequential presentation favors generalists almost to the exclusion of specialists, whereas Parallel presentation does the least so. However, other protocol attributes (the nature of dilution between transfers of phages to new cultures) also have effects on selection and phage maintenance. It is also noted that protocols can be designed to enhance recombination to augment evolution and to reduce stochastic loss of newly arisen mutants.

Ecological Approach to Understanding Superinfection Inhibition in Bacteriophage.

In microbial communities, viruses compete with each other for host cells to infect. As a consequence of competition for hosts, viruses evolve inhibitory mechanisms to suppress their competitors. One such mechanism is superinfection exclusion, in which a preexisting viral infection prevents a secondary infection. The bacteriophage ΦX174 exhibits a potential superinfection inhibition mechanism (in which secondary infections are either blocked or resisted) known as the reduction effect. In this auto-inhibitory phenomenon, a plasmid containing a fragment of the ΦX174 genome confers resistance to infection among cells that were once permissive to ΦX174. Taking advantage of this plasmid system, we examine the inhibitory properties of the ΦX174 reduction effect on a range of wild ΦX174-like phages. We then assess how closely the reduction effect in the plasmid system mimics natural superinfection inhibition by carrying out phage-phage competitions in continuous culture, and we evaluate whether the overall competitive advantage can be predicted by phage fitness or by a combination of fitness and reduction effect inhibition. Our results show that viral fitness often correctly predicts the winner. However, a phage's reduction sequence also provides an advantage to the phage in some cases, modulating phage-phage competition and allowing for persistence where competitive exclusion was expected. These findings provide strong evidence for more complex dynamics than were previously thought, in which the reduction effect may inhibit fast-growing viruses, thereby helping to facilitate coexistence.

ΦX174 Attenuation by Whole-Genome Codon Deoptimization.

Natural selection acting on synonymous mutations in protein-coding genes influences genome composition and evolution. In viruses, introducing synonymous mutations in genes encoding structural proteins can drastically reduce viral growth, providing a means to generate potent, live-attenuated vaccine candidates. However, an improved understanding of what compositional features are under selection and how combinations of synonymous mutations affect viral growth is needed to predictably attenuate viruses and make them resistant to reversion. We systematically recoded all nonoverlapping genes of the bacteriophage ΦX174 with codons rarely used in its Escherichia coli host. The fitness of recombinant viruses decreases as additional deoptimizing mutations are made to the genome, although not always linearly, and not consistently across genes. Combining deoptimizing mutations may reduce viral fitness more or less than expected from the effect size of the constituent mutations and we point out difficulties in untangling correlated compositional features. We test our model by optimizing the same genes and find that the relationship between codon usage and fitness does not hold for optimization, suggesting that wild-type ΦX174 is at a fitness optimum. This work highlights the need to better understand how selection acts on patterns of synonymous codon usage across the genome and provides a convenient system to investigate the genetic determinants of virulence.

Genic incompatibilities in two hybrid bacteriophages.

Horizontal gene transfer and recombination play a major role in microbial evolution and have been detected in diverse groups, including many of medical relevance such as HIV and dengue virus. In the absence of mechanistic barriers, the evolutionary success of a particular recombination event is determined by whether the recombinant genotype suffers a fitness cost through the disruption of favorable epistatic interactions within the genome, and if so, the extent to which this fitness cost might be mitigated by subsequent compensatory evolution. To investigate the importance of epistatic interactions between genes and the evolutionary viability of a homologous recombination event between diverged ancestral genotypes, we constructed two recombinant microvirid bacteriophages by exchanging their alleles of the gene encoding the coat protein. The coding sequences for this gene differ by approximately 8% at the amino acid level and were interchanged between two ancestral phages related to varphiX174 and well adapted to their culture conditions. Because the recombinant phages showed drastically reduced fitnesses, we further explored their evolutionary viability by subjecting replicate lines of each of them to selection. We found that all four lineages achieved fitnesses commensurate with ancestral fitnesses in as few as 60 generations, and on average, the first substitution accounted for more than half of the total fitness recovery. Fitness recovery required three to five substitutions in each lineage, and overall eight of the nine essential phage genes were involved, suggesting extensive epistatic interactions throughout the genome. Interestingly, the proteins with the most extensive and apparent physical interactions with the exchanged protein in the viral capsid did not appear to have much of a role in fitness recovery. This result appears to be a consequence of the conservation of the amino acid residues involved in the interactions. It suggests that strong epistatic interactions are less important than weaker, transient ones in producing genic incompatibilities because they preclude variability in the interacting regions of the proteins.

Viral adaptation to an antiviral protein enhances the fitness level to above that of the uninhibited wild type.

Viruses often evolve resistance to antiviral agents. While resistant strains are able to replicate in the presence of the agent, they generally exhibit lower fitness than the wild-type strain in the absence of the inhibitor. In some cases, resistant strains become dependent on the antiviral agent. However, the agent rarely, if ever, elevates dependent strain fitness above the uninhibited wild-type level. This would require an adaptive mechanism to convert the antiviral agent into a beneficial growth factor. Using an inhibitory scaffolding protein that specifically blocks phiX174 capsid assembly, we demonstrate that such mechanisms are possible. To obtain the quintuple-mutant resistant strain, the wild-type virus was propagated for approximately 150 viral life cycles in the presence of increasing concentrations of the inhibitory protein. The expression of the inhibitory protein elevated the strain's fitness significantly above the uninhibited wild-type level. Thus, selecting for resistance coselected for dependency, which was characterized and found to operate on the level of capsid nucleation. To the best of our knowledge, this is the first report of a virus evolving a mechanism to productively utilize an antiviral agent to stimulate its fitness above the uninhibited wild-type level. The results of this study may be predictive of the types of resistant phenotypes that could be selected by antiviral agents that specifically target capsid assembly.

Saturation Mutagenesis Genome Engineering of Infective ΦX174 Bacteriophage via Unamplified Oligo Pools and Golden Gate Assembly.

Here we present a novel protocol for the construction of saturation single-site-and massive multisite-mutant libraries of a bacteriophage. We segmented the ΦX174 genome into 14 nontoxic and nonreplicative fragments compatible with Golden Gate assembly. We next used nicking mutagenesis with oligonucleotides prepared from unamplified oligo pools with individual segments as templates to prepare near-comprehensive single-site mutagenesis libraries of genes encoding the F capsid protein (421 amino acids scanned) and G spike protein (172 amino acids scanned). Libraries possessed greater than 99% of all 11 860 programmed mutations. Golden Gate cloning was then used to assemble the complete ΦX174 mutant genome and generate libraries of infective viruses. This protocol will enable reverse genetics experiments for studying viral evolution and, with some modifications, can be applied for engineering therapeutically relevant bacteriophages with larger genomes.

The genetics of adaptation for eight microvirid bacteriophages.

Theories of adaptive molecular evolution have recently experienced significant expansion, and their predictions and assumptions have begun to be subjected to rigorous empirical testing. However, these theories focus largely on predicting the first event in adaptive evolution, the fixation of a single beneficial mutation. To address long-term adaptation it is necessary to include new assumptions, but empirical data are needed for guidance. To empirically characterize the general properties of adaptive walks, eight recently isolated relatives of the single-stranded DNA (ssDNA) bacteriophage phiX174 (family Microviridae) were adapted to identical selective conditions. Three of the eight genotypes were adapted in replicate, for a total of 11 adaptive walks. We measured fitness improvement and identified the genetic changes underlying the observed adaptation. Nearly all phages were evolvable; nine of the 11 lineages showed a significant increase in fitness. However, fitness plateaued quickly, and adaptation was achieved through only three substitutions on average. Parallel evolution was rampant, both across replicates of the same genotype as well as across different genotypes, yet adaptation of replicates never proceeded through the exact same set of mutations. Despite this, final fitnesses did not vary significantly among replicates. Final fitnesses did vary significantly across genotypes but not across phylogenetic groupings of genotypes. A positive correlation was found between the number of substitutions in an adaptive walk and the magnitude of fitness improvement, but no correlation was found between starting and ending fitness. These results provide an empirical framework for future adaptation theory.

Space, time, and host evolution facilitate coexistence of competing bacteriophages: theory and experiment.

We present a joint experimental/theoretical investigation into the roles of spatial structure and time in the competition between two pathogens for a single host. We suggest a natural mechanism by which competing pathogens can coexist when host evolution and competitive dynamics occur on similar timescales. Our experimental system consisted of a single bacterial host species and two competing bacteriophage strains grown on agar plates, with a serial transfer of samples of the bacteriophage population to fresh host populations after each incubation cycle. The experiments included two incubation times and two transfer protocols that either maintained or disrupted the spatial structure of the viruses at each transfer. The same bacteriophage acted as the dominant competitor under both transfer protocols. A striking difference between the treatments is that the weak competitor was able to persist in the long-incubation experiments but not in the short-incubation experiments. Mathematical and experimental evidence suggest that coexistence is due to the appearance of resistant mutant host cells that provide a transient "spatiotemporal refuge" for the weaker competitor. Our mathematical model is individual based, captures the stochastic spatial dynamics down to the level of individual cells, and helps to explain the differences in behavior under the various experimental conditions.

Loss of LINE-1 activity in the megabats.

LINE-1 (L1) retrotransposons are the most abundant type of mammalian retroelement. They have profound effects on genome plasticity and have been proposed to fulfill essential host functions, yet it remains unclear where they lie on the spectrum from parasitism to mutualism. Their ubiquity makes it difficult to determine the extent of their effects on genome evolution and gene expression because of the relative dearth of animal models lacking L1 activity. We have isolated L1 sequences from 11 megabat species by a method that enriches for recently inserted L1s and have done a bioinformatic examination of L1 sequences from a 12th species whose genome was recently shotgun sequenced. An L1 extinction event appears to have occurred at least 24 million years ago (MYA) in an ancestor of the megabats. The ancestor was unusual in having maintained two highly divergent long-term L1 lineages with different levels of activity, which appear, on an evolutionary scale, to have simultaneously lost that activity. These megabat species can serve as new animal models to ask what effect loss of L1 activity has on mammalian genome evolution and gene expression.

Tracing the history of LINE and SINE extinction in sigmodontine rodents.

BACKGROUND: L1 retrotransposons have co-evolved with their mammalian hosts for the entire history of mammals and currently compose ~ 20% of a mammalian genome. B1 retrotransposons are dependent on L1 for retrotransposition and span the evolutionary history of rodents since their radiation. L1s were found to have lost their activity in a group of South American rodents, the Sigmodontinae, and B1 inactivation preceded the extinction of L1 in the same group. Consequently, a basal group of sigmodontines have active L1s but inactive B1s and a derived clade have both inactive L1s and B1s. It has been suggested that B1s became extinct during a long period of L1 quiescence and that L1s subsequently reemerged in the basal group. RESULTS: Here we investigate the evolutionary histories of L1 and B1 in the sigmodontine rodents and show that L1 activity continued until after the L1-extinct clade and the basal group diverged. After the split, L1 had a small burst of activity in the former group, followed by extinction. In the basal group, activity was initially low but was followed by a dramatic increase in L1 activity. We found the last wave of B1 retrotransposition was large and probably preceded the split between the two rodent clades. CONCLUSIONS: Given that L1s had been steadily retrotransposing during the time corresponding to B1 extinction and that the burst of B1 activity preceding B1 extinction was large, we conclude that B1 extinction was not a result of L1 quiescence. Rather, the burst of B1 activity may have contributed to L1 extinction both by competition with L1 and by putting strong selective pressure on the host to control retrotransposition.