Genetic dissection of adaptive form and function in rapidly speciating cichlid fishes

Genes of major phenotypic effects and strong genetic correlations can facilitate adaptation, direct selective responses, and potentially lead to phenotypic convergence. However, the preponderance of this type of genetic architecture in repeatedly evolved adaptations remains unknown. Using hybrids between Haplochromis chilotes (thick‐lipped) and Pundamilia nyererei (thin‐lipped) we investigated the genetics underlying hypertrophied lips and elongated heads, traits that evolved repeatedly in cichlids. At least 25 loci of small‐to‐moderate and mainly additive effects were detected. Phenotypic variation in lip and head morphology was largely independent. Although several QTL overlapped for lip and head morphology traits, they were often of opposite effects. The distribution of effect signs suggests strong selection on lips. The fitness implications of several detected loci were demonstrated using a laboratory assay testing for the association between genotype and variation in foraging performance. The persistence of low fitness alleles in head morphology appears to be maintained through antagonistic pleiotropy/close linkage with positive‐effect lip morphology alleles. Rather than being based on few major loci with strong positive genetic correlations, our results indicate that the evolution of the Lake Victoria thick‐lipped ecomorph is the result of selection on numerous loci distributed throughout the genome.     

Reverse evolution of fitness in Drosophila melanogaster

Abstract The evolution of fitness is central to evolutionary theory, yet few experimental systems allow us to track its evolution in genetically and environmentally relevant contexts. Reverse evolution experiments allow the study of the evolutionary return to ancestral phenotypic states, including fitness. This in turn permits well‐defined tests for the dependence of adaptation on evolutionary history and environmental conditions. In the experiments described here, 20 populations of heterogeneous evolutionary histories were returned to their common ancestral environment for 50 generations, and were then compared with both their immediate differentiated ancestors and populations which had remained in the ancestral environment. One measure of fitness returned to ancestral levels to a greater extent than other characters did. The phenotypic effects of reverse evolution were also contingent on previous selective history. Moreover, convergence to the ancestral state was highly sensitive to environmental conditions. The phenotypic plasticity of fecundity, a character directly selected for, evolved during the experimental time frame. Reverse evolution appears to force multiple, diverged populations to converge on a common fitness state through different life‐history and genetic changes.     

The genetic basis of parallel and divergent phenotypic responses in evolving populations of Escherichia coli

Pleiotropy plays a central role in theories of adaptation, but little is known about the distribution of pleiotropic effects associated with different adaptive mutations. Previously, we described the phenotypic effects of a collection of independently arising beneficial mutations in Escherichia coli. We quantified their fitness effects in the glucose environment in which they evolved and their pleiotropic effects in five novel resource environments. Here we use a candidate gene approach to associate the phenotypic effects of the mutations with the underlying genetic changes. Among our collection of 27 adaptive mutants, we identified a total of 21 mutations (18 of which were unique) encompassing five different loci or gene regions. There was limited resolution to distinguish among loci based on their fitness effects in the glucose environment, demonstrating widespread parallelism in the direct response to selection. However, substantial heterogeneity in mutant effects was revealed when we examined their pleiotropic effects on fitness in the five novel environments. Substitutions in the same locus clustered together phenotypically, indicating concordance between molecular and phenotypic measures of divergence.     

THE EVOLUTION OF PHENOTYPIC CORRELATIONS AND “DEVELOPMENTAL MEMORY”

Development introduces structured correlations among traits that may constrain or bias the distribution of phenotypes produced. Moreover, when suitable heritable variation exists, natural selection may alter such constraints and correlations, affecting the phenotypic variation available to subsequent selection. However, exactly how the distribution of phenotypes produced by complex developmental systems can be shaped by past selective environments is poorly understood. Here we investigate the evolution of a network of recurrent nonlinear ontogenetic interactions, such as a gene regulation network, in various selective scenarios. We find that evolved networks of this type can exhibit several phenomena that are familiar in cognitive learning systems. These include formation of a distributed associative memory that can “store” and “recall” multiple phenotypes that have been selected in the past, recreate complete adult phenotypic patterns accurately from partial or corrupted embryonic phenotypes, and “generalize” (by exploiting evolved developmental modules) to produce new combinations of phenotypic features. We show that these surprising behaviors follow from an equivalence between the action of natural selection on phenotypic correlations and associative learning, well‐understood in the context of neural networks. This helps to explain how development facilitates the evolution of high‐fitness phenotypes and how this ability changes over evolutionary time.     

Fitness recovery and compensatory evolution in natural mutant lines of C. elegans

Deleterious mutation accumulation plays a central role in evolutionary genetics, conservation biology, human health, and evolutionary medicine (e.g., methods of viral attenuation for live vaccines). It is therefore important to understand whether and how quickly populations with accumulated deleterious mutational loads can recover fitness through adaptive evolution. We used laboratory experimental evolution with four long-term mutation-accumulation (MA) lines of Caenorhabditis elegans nematodes to study the dynamics of such fitness evolution. We previously showed that when homozygous mutant populations are evolved in large population sizes, they can rapidly achieve wild-type fitness through the accumulation of new beneficial or compensatory epistatic mutations. Here, we expand this approach to demonstrate that when replicate lineages are initiated from the same mutant genotype, phenotypic evolution is only sometimes repeatable. MA genotypes that recovered ancestral fitness in the previous experiment did not always do so here. Further, the pattern of adaptive evolution in independently evolved replicates was contingent upon the MA genotype and varied among fitness-related traits. Our findings suggest that new beneficial mutations can drive rapid fitness evolution, but that the adaptive process is rendered somewhat unpredictable by its susceptibility to chance events and sensitivity to the evolutionary history of the starting population.     

CONVERGENT AND CORRELATED EVOLUTION OF MAJOR LIFE‐HISTORY TRAITS IN THE ANGIOSPERM GENUS LEUCADENDRON (PROTEACEAE)

Natural selection is expected to cause convergence of life histories among taxa as well as correlated evolution of different life‐history traits. Here, we quantify the extent of convergence of five key life‐history traits (adult fire survival, seed storage, degree of sexual dimorphism, pollination mode, and seed‐dispersal mode) and test hypotheses about their correlated evolution in the genus Leucadendron (Proteaceae) from the fire‐prone South African fynbos. We reconstructed a new molecular phylogeny of this highly diverse genus that involves more taxa and molecular markers than previously. This reconstruction identifies new clades that were not detected by previous molecular study and morphological classifications. Using this new phylogeny and robust methods that account for phylogenetic uncertainty, we show that the five life‐history traits studied were labile during the evolutionary history of the genus. This diversity allowed us to tackle major questions about the correlated evolution of life‐history strategies. We found that species with longer seed‐dispersal distances tended to evolve lower pollen‐dispersal distance, that insect‐pollinated species evolved decreased sexual dimorphism, and that species with a persistent soil seed‐bank evolved toward reduced fire‐survival ability of adults.     

GENE FLOW AND SELECTION ON PHENOTYPIC PLASTICITY IN AN ISLAND SYSTEM OF RANA TEMPORARIA

Gene flow is often considered to be one of the main factors that constrains local adaptation in a heterogeneous environment. However, gene flow may also lead to the evolution of phenotypic plasticity. We investigated the effect of gene flow on local adaptation and phenotypic plasticity in development time in island populations of the common frog Rana temporaria which breed in pools that differ in drying regimes. This was done by investigating associations between traits (measured in a common garden experiment) and selective factors (pool drying regimes and gene flow from other populations inhabiting different environments) by regression analyses and by comparing pairwise F_ST values (obtained from microsatellite analyses) with pairwise Q_ST values. We found that the degree of phenotypic plasticity was positively correlated with gene flow from other populations inhabiting different environments (among‐island environmental heterogeneity), as well as with local environmental heterogeneity within each population. Furthermore, local adaptation, manifested in the correlation between development time and the degree of pool drying on the islands, appears to have been caused by divergent selection pressures. The local adaptation in development time and phenotypic plasticity is quite remarkable, because the populations are young (less than 300 generations) and substantial gene flow is present among islands.     

Mistranslation Reduces Mutation Load in Evolving Proteins through Negative Epistasis with DNA Mutations

Translational errors during protein synthesis cause phenotypic mutations that are several orders of magnitude more frequent than DNA mutations. Such phenotypic mutations may affect adaptive evolution through their interactions with DNA mutations. To study how mistranslation may affect the adaptive evolution of evolving proteins, we evolved populations of green fluorescent protein (GFP) in either high-mistranslation or low-mistranslation Escherichia coli hosts. In both hosts, we first evolved GFP under purifying selection for the ancestral phenotype green fluorescence, and then under directional selection toward the new phenotype yellow fluorescence. High-mistranslation populations evolved modestly higher yellow fluorescence during each generation of evolution than low-mistranslation populations. We demonstrate by high-throughput sequencing that elevated mistranslation reduced the accumulation of deleterious DNA mutations under both purifying and directional selection. It did so by amplifying the fitness effects of deleterious DNA mutations through negative epistasis with phenotypic mutations. In contrast, mistranslation did not affect the incidence of beneficial mutations. Our findings show that phenotypic mutations interact epistatically with DNA mutations. By reducing a population’s mutation load, mistranslation can affect an important determinant of evolvability.     

Parasite‐mediated selection in a natural metapopulation of Daphnia magna

Parasite‐mediated selection varying across time and space in metapopulations is expected to result in host local adaptation and the maintenance of genetic diversity in disease‐related traits. However, nonadaptive processes like migration and extinction‐(re)colonization dynamics might interfere with adaptive evolution. Understanding how adaptive and nonadaptive processes interact to shape genetic variability in life‐history and disease‐related traits can provide important insights into their evolution in subdivided populations. Here we investigate signatures of spatially fluctuating, parasite‐mediated selection in a natural metapopulation of Daphnia magna. Host genotypes from infected and uninfected populations were genotyped at microsatellite markers, and phenotyped for life‐history and disease traits in common garden experiments. Combining phenotypic and genotypic data a Q_ST–F_ST‐like analysis was conducted to test for signatures of parasite mediated selection. We observed high variation within and among populations for phenotypic traits, but neither an indication of host local adaptation nor a cost of resistance. Infected populations have a higher gene diversity (Hs) than uninfected populations and Hs is strongly positively correlated with fitness. These results suggest a strong parasite effect on reducing population level inbreeding. We discuss how stochastic processes related to frequent extinction‐(re)colonization dynamics as well as host and parasite migration impede the evolution of resistance in the infected populations. We suggest that the genetic and phenotypic patterns of variation are a product of dynamic changes in the host gene pool caused by the interaction of colonization bottlenecks, inbreeding, immigration, hybrid vigor, rare host genotype advantage and parasitism. Our study highlights the effect of the parasite in ameliorating the negative fitness consequences caused by the high drift load in this metapopulation.     

JOINT ALLELIC EFFECTS ON FITNESS AND METRIC TRAITS

Theoretical explanations of empirically observed standing genetic variation, mutation, and selection suggest that many alleles must jointly affect fitness and metric traits. However, there are few direct demonstrations of the nature and extent of these pleiotropic associations. We implemented a mutation accumulation (MA) divergence experimental design in Drosophila serrata to segregate genetic variants for fitness and metric traits. By exploiting naturally occurring MA line extinctions as a measure of line‐level total fitness, manipulating sexual selection, and measuring productivity we were able to demonstrate genetic covariance between fitness and standard metric traits, wing size, and shape. Larger size was associated with lower total fitness and male sexual fitness, but higher productivity. Multivariate wing shape traits, capturing major axes of wing shape variation among MA lines, evolved only in the absence of sexual selection, and to the greatest extent in lines that went extinct, indicating that mutations contributing wing shape variation also typically had deleterious effects on both total fitness and male sexual fitness. This pleiotropic covariance of metric traits with fitness will drive their evolution, and generate the appearance of selection on the metric traits even in the absence of a direct contribution to fitness.     

Pervasive,yet idiosyncratic,epistatic pleiotropy during adaptation in a behaviourally complex microbe

Understanding how multiple mutations interact to jointly impact multiple ecologically important traits is critical for creating a robust picture of organismal fitness and the process of adaptation. However, this is complicated by both environmental heterogeneity and the complexity of genotype‐to‐phenotype relationships generated by pleiotropy and epistasis. Moreover, little is known about how pleiotropic and epistatic relationships themselves change over evolutionary time. The soil bacterium Myxococcus xanthus employs several distinct social traits across a range of environments. Here, we use an experimental lineage of M. xanthus that evolved a novel form of social motility to address how interactions between epistasis and pleiotropy evolve. Specifically, we test how mutations accumulated during selection on soft agar pleiotropically affect several other social traits (hard agar motility, predation and spore production). Relationships between changes in swarming rate in the selective environment and the four other traits varied greatly over time in both direction and magnitude, both across timescales of the entire evolutionary lineage and individual evolutionary time steps. We also tested how a previously defined epistatic interaction is pleiotropically expressed across these traits. We found that phenotypic effects of this epistatic interaction were highly correlated between soft and hard agar motility, but were uncorrelated between soft agar motility and predation, and inversely correlated between soft agar motility and spore production. Our results show that ‘epistatic pleiotropy’ varied greatly in magnitude, and often even in sign, across traits and over time, highlighting the necessity of simultaneously considering the interacting complexities of pleiotropy and epistasis when studying the process of adaptation.     

Genotype‐by‐environment interactions due to antibiotic resistance and adaptation in Escherichia coli

Mutations that are beneficial in one environment can have different fitness effects in other environments. In the context of antibiotic resistance, the resulting genotype‐by‐environment interactions potentially make selection on resistance unpredictable in heterogeneous environments. Furthermore, resistant bacteria frequently fix additional mutations during evolution in the absence of antibiotics. How do these two types of mutations interact to determine the bacterial phenotype across different environments? To address this, I used Escherichia coli as a model system, measuring the effects of nine different rifampicin resistance mutations on bacterial growth in 31 antibiotic‐free environments. I did this both before and after approximately 200 generations of experimental evolution in antibiotic‐free conditions (LB medium), and did the same for the antibiotic‐sensitive wild type after adaptation to the same environment. The following results were observed: (i) bacteria with and without costly resistance mutations adapted to experimental conditions and reached similar levels of competitive fitness; (ii) rifampicin resistance mutations and adaptation to LB both indirectly altered growth in other environments; and (iii) resistant‐evolved genotypes were more phenotypically different from the ancestor and from each other than resistant‐nonevolved and sensitive‐evolved genotypes. This suggests genotype‐by‐environment interactions generated by antibiotic resistance mutations, observed previously in short‐term experiments, are more pronounced after adaptation to other types of environmental variation, making it difficult to predict long‐term selection on resistance mutations from fitness effects in a single environment.     

Apparent directional selection by biased pleiotropic mutation

Pleiotropic effects of deleterious mutations are considered to be among the factors responsible for genetic constraints on evolution by long-term directional selection acting on a quantitative trait. If pleiotropic phenotypic effects are biased in a particular direction, mutations generate apparent directional selection, which refers to the covariance between fitness and the trait owing to a linear association between the number of mutations possessed by individuals and the genotypic values of the trait. The present analysis has shown how the equilibrium mean value of the trait is determined by a balance between directional selection and biased pleiotropic mutations. Assuming that genes act additively both on the trait and on fitness, the total variance-standardized directional selection gradient was decomposed into apparent and true components. Experimental data on mutation bias from the bristle traits of Drosophila and life history traits of Daphnia suggest that apparent selection explains a small but significant fraction of directional selection pressure that is observed in nature; the data suggest that changes induced in a trait by biased pleiotropic mutation (i.e., by apparent directional selection) are easily compensated for by (true) directional selection.     

GENETIC DRIFT IN ANTAGONISTIC GENES LEADS TO DIVERGENCE IN SEX‐SPECIFIC FITNESS BETWEEN EXPERIMENTAL POPULATIONS OF DROSOPHILA MELANOGASTER

Males and females differ in their reproductive roles and as a consequence are often under diverging selection pressures on shared phenotypic traits. Theory predicts that divergent selection can favor the invasion of sexually antagonistic alleles, which increase the fitness of one sex at the detriment of the other. Sexual antagonism can be subsequently resolved through the evolution of sex‐specific gene expression, allowing the sexes to diverge phenotypically. Although sexual dimorphism is very common, recent evidence also shows that antagonistic genetic variation continues to segregate in populations of many organisms. Here we present empirical data on the interaction between sexual antagonism and genetic drift in populations that have independently evolved under standardized conditions. We demonstrate that small experimental populations of Drosophila melanogaster have diverged in male and female fitness, with some populations showing high male, but low female fitness while other populations show the reverse pattern. The between‐population patterns are consistent with the differentiation in reproductive fitness being driven by genetic drift in sexually antagonistic alleles. We discuss the implications of our results with respect to the maintenance of antagonistic variation in subdivided populations and consider the wider implications of drift in fitness‐related genes.     

Phenotypic integration between claw and toepad traits promotes microhabitat specialization in the Anolis adaptive radiation

The performance of an organism in its environment frequently depends more on its composite phenotype than on individual phenotypic traits. Thus, understanding environmental adaptation requires investigating patterns of covariation across functionally related traits. The replicated adaptive radiations of Greater Antillean Anolis lizards are characterized by ecological and morphological convergence, thus, providing an opportunity to examine the role of multiple phenotypes in microhabitat adaptation. Here, we examine integrated claw and toepad morphological evolution in relation to habitat partitioning across the adaptive radiations of Greater Antillean anoles. Based on analysis of 428 specimens from 57 species, we found that different aspects of claw morphology were associated with different perch dimensions, with claw height positively associated with perch diameter and claw curvature positively associated with perch height. Patterns of integration also varied across claw and toepad traits, likely driven by correlative selection for performance on smoother and rougher substrates. Finally, rates of evolution differed between claw and toepad traits, with claw length evolving faster than all other traits despite having no predicted functional importance. Our results highlight the multivariate nature of phenotypic adaptation and suggest that phenotypic integration across Greater Antillean anoles is driven by fine‐scale correlative selection based on structural habitat specialization.     

Evolutionary genomics of host adaptation in vesicular stomatitis virus

Populations experiencing similar selection pressures can sometimes diverge in the genetic architectures underlying evolved complex traits. We used RNA virus populations of large size and high mutation rate to study the impact of historical environment on genome evolution, thus increasing our ability to detect repeatable patterns in the evolution of genetic architecture. Experimental vesicular stomatitis virus populations were evolved on HeLa cells, on MDCK cells, or on alternating hosts. Turner and Elena (2000. Cost of host radiation in an RNA virus. Genetics. 156:1465-1470.) previously showed that virus populations evolved in single-host environments achieved high fitness on their selected hosts but failed to increase in fitness relative to their ancestor on the unselected host and that alternating-host-evolved populations had high fitness on both hosts. Here we determined the complete consensus sequence for each evolved population after 95 generations to gauge whether the parallel phenotypic changes were associated with parallel genomic changes. We also analyzed the patterns of allele substitutions to discern whether differences in fitness across hosts arose through true pleiotropy or the presence of not only a mutation that is beneficial in both hosts but also 1 or more mutations at other loci that are costly in the unselected environment (mutation accumulation [MA]). We found that ecological history may influence to what extent pleiotropy and MA contribute to fitness asymmetries across environments. We discuss the degree to which current genetic architecture is expected to constrain future evolution of complex traits, such as host use by RNA viruses.     

The emotion system promotes diversity and evolvability

Studies on the relationship between the optimal phenotype and its environment have had limited focus on genotype-to-phenotype pathways and their evolutionary consequences. Here, we study how multi-layered trait architecture and its associated constraints prescribe diversity. Using an idealized model of the emotion system in fish, we find that trait architecture yields genetic and phenotypic diversity even in absence of frequency-dependent selection or environmental variation. That is, for a given environment, phenotype frequency distributions are predictable while gene pools are not. The conservation of phenotypic traits among these genetically different populations is due to the multi-layered trait architecture, in which one adaptation at a higher architectural level can be achieved by several different adaptations at a lower level. Our results emphasize the role of convergent evolution and the organismal level of selection. While trait architecture makes individuals more constrained than what has been assumed in optimization theory, the resulting populations are genetically more diverse and adaptable. The emotion system in animals may thus have evolved by natural selection because it simultaneously enhances three important functions, the behavioural robustness of individuals, the evolvability of gene pools and the rate of evolutionary innovation at several architectural levels.     

Rapid,broad‐scale gene expression evolution in experimentally harvested fish populations

Gene expression changes potentially play an important role in adaptive evolution under human‐induced selection pressures, but this has been challenging to demonstrate in natural populations. Fishing exhibits strong selection pressure against large body size, thus potentially inducing evolutionary changes in life history and other traits that may be slowly reversible once fishing ceases. However, there is a lack of convincing examples regarding the speed and magnitude of fisheries‐induced evolution, and thus, the relevant underlying molecular‐level effects remain elusive. We use wild‐origin zebrafish (Danio rerio) as a model for harvest‐induced evolution. We experimentally demonstrate broad‐scale gene expression changes induced by just five generations of size‐selective harvesting, and limited genetic convergence following the cessation of harvesting. We also demonstrate significant allele frequency changes in genes that were differentially expressed after five generations of size‐selective harvesting. We further show that nine generations of captive breeding induced substantial gene expression changes in control stocks likely due to inadvertent selection in the captive environment. The large extent and rapid pace of the gene expression changes caused by both harvest‐induced selection and captive breeding emphasizes the need for evolutionary enlightened management towards sustainable fisheries.     

Evolution of Drosophila resistance against different pathogens and infection routes entails no detectable maintenance costs

Pathogens exert a strong selective pressure on hosts, entailing host adaptation to infection. This adaptation often affects negatively other fitness‐related traits. Such trade‐offs may underlie the maintenance of genetic diversity for pathogen resistance. Trade‐offs can be tested with experimental evolution of host populations adapting to parasites, using two approaches: (1) measuring changes in immunocompetence in relaxed‐selection lines and (2) comparing life‐history traits of evolved and control lines in pathogen‐free environments. Here, we used both approaches to examine trade‐offs in Drosophila melanogaster populations evolving for over 30 generations under infection with Drosophila C Virus or the bacterium Pseudomonas entomophila, the latter through different routes. We find that resistance is maintained after up to 30 generations of relaxed selection. Moreover, no differences in several classical life‐history traits between control and evolved populations were found in pathogen‐free environments, even under stresses such as desiccation, nutrient limitation, and high densities. Hence, we did not detect any maintenance costs associated with resistance to pathogens. We hypothesize that extremely high selection pressures commonly used lead to the disproportionate expression of costs relative to their actual occurrence in natural systems. Still, the maintenance of genetic variation for pathogen resistance calls for an explanation.