Shifting the paradigm in Evolve and Resequence studies: From analysis of single nucleotide polymorphisms to selected haplotype blocks

For almost a decade the combination of whole genome sequencing with experimental evolution (Evolve and Resequence, E&R; Turner, Stewart, Fields, Rice, & Tarone, 2011) has been used to study adaptation in outcrossing organisms. However, complications caused by inversions and hitchhiking variants have prevented this powerful approach from living up to its potential. In this issue of Molecular Ecology, Michalak, Kang, Schou, Garner, and Loeschke (2018), provide an important step ahead by using a population of Drosophila melanogaster devoid of segregating inversions to identify the genetic basis of resistance to five environmental stressors. They further address the challenge of hitchhiking variants by reconstructing selected haplotype blocks. While it is apparent that the haplotype block reconstruction needs further refinements, their work underpins the potential of E&R studies in Drosophila to address fundamental questions in evolutionary biology.     

Haplotype‐based genotyping‐by‐sequencing in oat genome research

In a de novo genotyping‐by‐sequencing (GBS) analysis of short, 64‐base tag‐level haplotypes in 4657 accessions of cultivated oat, we discovered 164741 tag‐level (TL) genetic variants containing 241224 SNPs. From this, the marker density of an oat consensus map was increased by the addition of more than 70000 loci. The mapped TL genotypes of a 635‐line diversity panel were used to infer chromosome‐level (CL) haplotype maps. These maps revealed differences in the number and size of haplotype blocks, as well as differences in haplotype diversity between chromosomes and subsets of the diversity panel. We then explored potential benefits of SNP vs. TL vs. CL GBS variants for mapping, high‐resolution genome analysis and genomic selection in oats. A combined genome‐wide association study (GWAS) of heading date from multiple locations using both TL haplotypes and individual SNP markers identified 184 significant associations. A comparative GWAS using TL haplotypes, CL haplotype blocks and their combinations demonstrated the superiority of using TL haplotype markers. Using a principal component‐based genome‐wide scan, genomic regions containing signatures of selection were identified. These regions may contain genes that are responsible for the local adaptation of oats to Northern American conditions. Genomic selection for heading date using TL haplotypes or SNP markers gave comparable and promising prediction accuracies of up to r = 0.74. Genomic selection carried out in an independent calibration and test population for heading date gave promising prediction accuracies that ranged between r = 0.42 and 0.67. In conclusion, TL haplotype GBS‐derived markers facilitate genome analysis and genomic selection in oat.     

What paths do advantageous alleles take during short‐term evolutionary change?

Combining experimental evolution with whole‐genome resequencing is a promising new strategy for investigating the dynamics of evolutionary change. Published studies that have resequenced laboratory‐selected populations of sexual organisms have typically focused on populations sampled at the end of an evolution experiment. These studies have attempted to associate particular alleles with phenotypic change and attempted to distinguish between different theoretical models of adaptation. However, neither the population used to initiate the experiment nor multiple time points sampled during the evolutionary trajectory are generally available for examination. In this issue of Molecular Ecology, Orozco‐terWengel et al. (2012) take a significant step forward by estimating genome‐wide allele frequencies at the start, 15 generations into and at the end of a 37‐generation Drosophila experimental evolution study. The authors identify regions of the genome that have responded to laboratory selection and describe the temporal dynamics of allele frequency change. They identify two common trajectories for putatively adaptive alleles: alleles either gradually increase in frequency throughout the entire 37 generations or alleles plateau at a new frequency by generation 15. The identification of complex trajectories of alleles under selection contributes to a growing body of literature suggesting that simple models of adaptation, whereby beneficial alleles arise and increase in frequency unimpeded until they become fixed, may not adequately describe short‐term response to selection.     

Artificial selection on larval growth curves in Tribolium: correlated responses and constraints

Body size is often constrained from evolving. Although artificial selection on body size in insects frequently results in a sizable response, these responses usually bear fitness costs. Further, these experiments tend to select only on size at one landmark age, rather than selecting for patterns of growth over the whole larval life stage. To address whether constraints may be caused by larval growth patterns rather than final size, we implemented a function‐valued (FV) trait method of selection, in which entire larval growth curves from Tribolium were artificially selected. The selection gradient function used was previously predicted to give the maximal response and was implemented using a novel selection index in the FV framework. Results indicated a significant response after one generation of selection, but no response in subsequent generations. Correlated responses included increased mortality, increased critical weight, and decreased development time (DT). The lack of response in size and development time after the first generation was likely caused by increased mortality suffered in selected lines; we demonstrated that the selection criterion caused both increased body size and increased mortality. We conclude that artificial selection on continuous traits using FV methods is very efficient and that the constraint of body size evolution is likely caused by a suite of trade‐offs with other traits.     

Cadmium‐induced and trans‐generational changes in the cultivable and total seed endophytic community of Arabidopsis thaliana

Trans‐generational adaptation is important to respond rapidly to environmental challenges and increase overall plant fitness. Besides well‐known mechanisms such as epigenetic modifications, vertically transmitted endophytic bacteria might contribute to this process. The cultivable and total endophytic communities of several generations of Arabidopsis thaliana seeds harvested from plants exposed to cadmium (Cd) or not exposed were investigated. The diversity and richness of the seed endophytic community decreased with an increasing number of generations. Aeromicrobium and Pseudonocardia were identified as indicator species in seeds from Cd‐exposed plants, while Rhizobium was abundantly present in both seed types. Remarkably, Rhizobium was the only genus that was consistently detected in seeds of all generations, which suggests that the phenotypic characteristics were more important as selection criteria for which bacteria are transferred to the next plant generation than the actual genera. Production of IAA was an important trait for endophytes from both seed types, while ACC deaminase activity and Cd tolerance were mainly associated with seed endophytes from Cd‐exposed plants. Understanding how different factors influence the seed endophytic community can help us to improve seed quality and plant growth through different biotechnological applications.     

Mitochondrial DNA and Y‐chromosome diversity in East Adriatic sheep

Variation in mitochondrial DNA (mtDNA) and Y‐chromosome haplotypes was analysed in nine domestic sheep breeds (159 rams) and 21 mouflon ( Ovis musimon) sampled in the East Adriatic. Mitochondrial DNA analyses revealed a high frequency of type B haplotypes, predominantly in European breeds, and a very low frequency of type A haplotypes, which are more frequent in some Asian breeds. Mitochondrial haplotype Hmt‐3 was the most frequent (26.4%), and 37.1%, 20.8% and 7.6% of rams had haplotypes one, two and three mutations remote from Hmt‐3 respectively. In contrast, Y‐chromosome analyses revealed extraordinary paternal allelic richness: HY‐6, 89.3%; HY‐8, 5.0%; HY‐18, 3.1%; HY‐7, 1.3%; and HY‐5, 1.3%. In fact, the number of haplotypes observed is comparable to the number found in Turkish breeds and greater than the number found in European breeds so far. Haplotype HY‐18 (A‐oY1/135‐SRYM18), identified here for the first time, provides a link between the haplotype HY‐12 (A‐oY1/139‐SRYM18) found in a few rams in Turkey and haplotype HY‐9 (A‐oY1/131‐SRYM18) found in one ram in Ethiopia. All mouflons had type B mtDNA haplotypes, including the private haplotype (Hmt‐55), and all were paternally monomorphic for haplotype HY‐6. Our data support a quite homogeneous maternal origin of East Adriatic sheep, which is a characteristic of European breeds. At the same time, the high number of haplotypes found was surprising and intriguing, and it begs for further analysis. Simultaneous analysis of mtDNA and Y‐chromosome information allowed us to detect a large discrepancy between maternal and paternal lineages in some populations. This is most likely the result of breeder efforts to ‘upgrade’ local populations using rams with different paternal origins.     

Death feigning as an adaptive anti‐predator behaviour: Further evidence for its evolution from artificial selection and natural populations

Death feigning is considered to be an adaptive antipredator behaviour. Previous studies on Tribolium castaneum have shown that prey which death feign have a fitness advantage over those that do not when using a jumping spider as the predator. Whether these effects are repeatable across species or whether they can be seen in nature is, however, unknown. Therefore, the present study involved two experiments: (a) divergent artificial selection for the duration of death feigning using a related species T. freemani as prey and a predatory bug as predator, demonstrating that previous results are repeatable across both prey and predator species, and (b) comparison of the death‐feigning duration of T. castaneum populations collected from field sites with and without predatory bugs. In the first experiment, T. freemani adults from established selection regimes with longer durations of death feigning had higher survival rates and longer latency to being preyed on when they were placed with predatory bugs than the adults from regimes selected for shorter durations of death feigning. As a result, the adaptive significance of death‐feigning behaviour was demonstrated in another prey–predator system. In the second experiment, wild T. castaneum beetles from populations with predators feigned death longer than wild beetles from predator‐free populations. Combining the results from these two experiments with those from previous studies provided strong evidence that predators drive the evolution of longer death feigning.     

Experimental evolution of competing bean beetle species reveals long‐term reversals of short‐term evolution,but no consistent character displacement

Interspecific competition for shared resources should select for evolutionary divergence in resource use between competing species, termed character displacement. Many purported examples of character displacement exist, but few completely rule out alternative explanations. We reared genetically diverse populations of two species of bean beetles, Callosobruchus maculatus and Callosobruchus chinensis, in allopatry and sympatry on a mixture of adzuki beans and lentils, and assayed oviposition preference and other phenotypic traits after four, eight, and twelve generations of (co)evolution. C. maculatus specializes on adzuki beans; the generalist C. chinensis uses both beans. C. chinensis growing in allopatry emerged equally from both bean species. In sympatry, the two species competing strongly and coexisted via strong realized resource partitioning, with C. chinensis emerging almost exclusively from lentils and C. maculatus emerging almost exclusively from adzuki beans. However, oviposition preferences, larval survival traits, and larval development rates in both beetle species did not vary consistently between allopatric versus sympatric treatments. Rather, traits evolved in treatment‐independent fashion, with several traits exhibiting reversals in their evolutionary trajectories. For example, C. chinensis initially evolved a slower egg‐to‐adult development rate on adzuki beans in both allopatry and sympatry, then subsequently evolved back toward the faster ancestral development rate. Lack of character displacement is consistent with a previous similar experiment in bean beetles and may reflect lack of evolutionary trade‐offs in resource use. However, evolutionary reversals were unexpected and remain unexplained. Together with other empirical and theoretical work, our results illustrate the stringency of the conditions for character displacement.     

Loss of adaptation following reversion suggests trade‐offs in host use by a seed beetle

Experimental evolution has provided little support for the hypothesis that the narrow diets of herbivorous insects reflect trade‐offs in performance across hosts; selection lines can sometimes adapt to an inferior novel host without a decline in performance on the ancestral host. An alternative approach for detecting trade‐offs would be to measure adaptation decay after selection is relaxed, that is, when populations newly adapted to a novel host are reverted to the ancestral one. Lines of the seed beetle Callosobruchus maculatus rapidly adapted to a poor host (lentil); survival in lentil seeds increased from 2% to > 90% in < 30 generations. After the lines had reached a plateau with respect to survival in lentil, sublines were reverted to the ancestral host, mung bean. Twelve generations of reversion had little effect on performance in lentil, but after 25–35 generations, the reverted lines exhibited lower survival, slower development and smaller size. The most divergent pair of lines was then assayed on both lentil and mung bean. Performance on lentil was again much poorer in the reverted line than in the nonreverted one, but the lines performed equally well on mung bean. Moreover, the performance of the nonreverted line on mung bean remained comparable to that of the original mung‐bean population. Our results thus present a paradox: loss of adaptation to lentil following reversion implies a trade‐off, but the continued strong performance of lentil‐adapted lines on mung bean does not. Genomic comparisons of the reverted, nonreverted and ancestral lines may resolve this paradox and determine the importance of selection vs. drift in causing a loss of adaptation following reversion.     

Modelling as a tool for analysing the temperature‐dependent future of the Colorado potato beetle in Europe

A warmer climate may increase the risk of attacks by insect pests on agricultural crops, and questions on how to adapt management practice have created a need for impact models. Phenological models driven by climate data can be used for assessing the potential distribution and voltinism of different insect species, but the quality of the simulations is influenced by a range of uncertainties. In this study, we model the temperature‐dependent activity and development of the Colorado potato beetle, and analyse the influence of uncertainty associated with parameterization of temperature and day length response. We found that the developmental threshold has a major impact on the simulated number of generations per year. Little is known about local adaptations and individual variations, but the use of an upper and a lower developmental threshold gave an indication on the potential variation. The day length conditions triggering diapause are known only for a few populations. We used gridded observed temperature data to estimate local adaptations, hypothesizing that cold autumns can leave a footprint in the population genetics by low survival of individuals not reaching the adult stage before winter. Our study indicated that the potential selection pressure caused by climate conditions varies between European regions. Provided that there is enough genetic variation, a local adaption at the northern distribution limit would reduce the number of unsuccessful initiations and thereby increase the potential for spreading to areas currently not infested. The simulations of the impact model were highly sensitive to biases in climate model data, i.e. systematic deviations in comparison with observed weather, highlightening the need of improved performance of regional climate models. Even a moderate temperature increase could change the voltinism of Leptinotarsa decemlineata in Europe, but knowledge on agricultural practice and strategies for countermeasures is needed to evaluate changes in risk of attacks.     

Fine‐grained adaptive divergence in an amphibian: genetic basis of phenotypic divergence and the role of nonrandom gene flow in restricting effective migration among wetlands

Adaptive ecological differentiation among sympatric populations is promoted by environmental heterogeneity, strong local selection and restricted gene flow. High gene flow, on the other hand, is expected to homogenize genetic variation among populations and therefore prevent local adaptation. Understanding how local adaptation can persist at the spatial scale at which gene flow occurs has remained an elusive goal, especially for wild vertebrate populations. Here, we explore the roles of natural selection and nonrandom gene flow (isolation by breeding time and habitat choice) in restricting effective migration among local populations and promoting generalized genetic barriers to neutral gene flow. We examined these processes in a network of 17 breeding ponds of the moor frog Rana arvalis, by combining environmental field data, a common garden experiment and data on variation in neutral microsatellite loci and in a thyroid hormone receptor (TRβ) gene putatively under selection. We illustrate the connection between genotype, phenotype and habitat variation and demonstrate that the strong differences in larval life history traits observed in the common garden experiment can result from adaptation to local pond characteristics. Remarkably, we found that haplotype variation in the TRβ gene contributes to variation in larval development time and growth rate, indicating that polymorphism in the TRβ gene is linked with the phenotypic variation among the environments. Genetic distance in neutral markers was correlated with differences in breeding time and environmental differences among the ponds, but not with geographical distance. These results demonstrate that while our study area did not exceed the scale of gene flow, ecological barriers constrained gene flow among contrasting habitats. Our results highlight the roles of strong selection and nonrandom gene flow created by phenological variation and, possibly, habitat preferences, which together maintain genetic and phenotypic divergence at a fine‐grained spatial scale.     

Distinct haplotype structure at the innate immune receptor Toll‐like receptor 2 across bank vole populations and lineages in Europe

Parasite‐mediated selection may contribute to the maintenance of genetic variation at host immune genes over long time scales. To date, the best evidence for the long‐term maintenance of immunogenetic variation in natural populations comes from studies on the major histocompatibility complex (MHC) genes, whereas evidence for such processes from other immune genes remains scarce. In the present study, we show that, despite pronounced population differentiation and the occurrence of numerous private alleles within populations, the innate immune gene Toll‐like receptor 2 (TLR2) displays a distinct haplotype structure in 21 bank vole (Myodes glareolus) populations across Europe. Haplotypes from all populations grouped in four clearly differentiated clusters, with the three main clusters co‐occurring in at least three previously described mitochondrial lineages. This pattern indicates that the distinct TLR2 haplotype structure may precede the split of the mitochondrial lineages 0.19–0.56 Mya and suggests that haplotype clusters at this innate immune receptor are maintained over prolonged time in wild bank vole populations.     

Long‐term evolution of the natural isolate of Escherichia coli 536 in the mouse gut colonized after maternal transmission reveals convergence in the constitutive expression of the lactose operon

In vitro experimental evolution has taught us many lessons on the molecular bases of adaptation. To move towards more natural settings, evolution in the mice gut has been successfully performed. Yet, these experiments suffered from the use of laboratory strains as well as the use of axenic or streptomycin‐treated mice to maintain the inoculated strains. To circumvent these limitations, we conducted a one‐year experimental evolution in vivo using a natural isolate of E. coli, strain 536, in conditions mimicking as much as possible natural environment with mother‐to‐offspring microbiota transmission. Mice were then distributed in 24 independent cages and separated into two different diets: a regular one (chow diet, CD) and high‐fat and high‐sugar one (Western Diet, WD). Genome sequences revealed an early and rapid selection during the breastfeeding period that selected the constitutive expression of the well‐characterized lactose operon. E. coli was lost significantly more in CD than WD; however, we could not detect any genomic signature of selection, nor any diet specificities during the later part of the experiments. The apparently neutral evolution presumably due to low population size maintained nevertheless at high frequency the early selected mutations affecting lactose regulation. The rapid loss of lactose operon regulation challenges the idea that plastic gene expression is both optimal and stable in the wild.     

Evolutionary history of a dispersal‐associated locus across sympatric and allopatric divergent populations of a wing‐polymorphic beetle across Atlantic Europe

Studying the evolutionary history of trait divergence, in particular those related to dispersal capacity, is of major interest for the process of local adaptation and metapopulation dynamics. Here, we reconstruct the evolution of different alleles at the nuclear‐encoded mitochondrial NADP⁺‐dependent isocitrate dehydrogenase (mtIdh) locus of the ground beetle Pogonus chalceus that are differentially and repeatedly selected in short‐ and long‐winged populations in response to different hydrological regimes at both allopatric and sympatric scales along the Atlantic European coasts. We sequenced 2788 bp of the mtIdh locus spanning a ~7‐kb genome region and compared its variation with that of two supposedly neutral genes. mtIdh sequences show (i) monophyletic clustering of the short‐winged associated mtIDH‐DE haplotypes within the long‐winged associated mtIDH‐AB haplotypes, (ii) a more than tenfold lower haplotype diversity associated with the mtIDH‐DE alleles compared to the mtIDH‐AB alleles and (iii) a high number of fixed nucleotide differences between both mtIDH haplotype clusters. Coalescent simulations suggest that this observed sequence variation in the mtIdh locus is most consistent with a singular origin in a partially isolated subpopulation, followed by a relatively recent spread of the mtIDH‐DE allele in short‐winged populations along the Atlantic coast. These results demonstrate that even traits associated with decreased dispersal capacity can rapidly spread and that reuse of adaptive alleles plays an important role in the adaptive potential within this sympatric mosaic of P. chalceus populations.     

Genetic changes in a novel breeding population of Brassica napus synthesized from hundreds of crosses between B. rapa and B. carinata

Introgression of genomic variation between and within related crop species is a significant evolutionary approach for population differentiation, genome reorganization and trait improvement. Using the Illumina Infinium Brassica 60K SNP array, we investigated genomic changes in a panel of advanced generation new‐type Brassica napus breeding lines developed from hundreds of interspecific crosses between 122 Brassica rapa and 74 Brassica carinata accessions, and compared them with representative accessions of their three parental species. The new‐type B. napus population presented rich genetic diversity and abundant novel genomic alterations, consisting of introgressions from B. rapa and B. carinata, novel allelic combinations, reconstructed linkage disequilibrium patterns and haplotype blocks, and frequent deletions and duplications (nonrandomly distributed), particularly in the C subgenome. After a much shorter, but very intensive, selection history compared to traditional B. napus, a total of 15 genomic regions with strong selective sweeps and 112 genomic regions with putative signals of selective sweeps were identified. Some of these regions were associated with important agronomic traits that were selected for during the breeding process, while others were potentially associated with restoration of genome stability and fertility after interspecific hybridization. Our results demonstrate how a novel method for population‐based crop genetic improvement can lead to rapid adaptation, restoration of genome stability and positive responses to artificial selection.     

Rapid genome‐wide evolution in Brassica rapa populations following drought revealed by sequencing of ancestral and descendant gene pools

There is increasing evidence that evolution can occur rapidly in response to selection. Recent advances in sequencing suggest the possibility of documenting genetic changes as they occur in populations, thus uncovering the genetic basis of evolution, particularly if samples are available from both before and after selection. Here, we had a unique opportunity to directly assess genetic changes in natural populations following an evolutionary response to a fluctuation in climate. We analysed genome‐wide differences between ancestors and descendants of natural populations of Brassica rapa plants from two locations that rapidly evolved changes in multiple phenotypic traits, including flowering time, following a multiyear late‐season drought in California. These ancestor‐descendant comparisons revealed evolutionary shifts in allele frequencies in many genes. Some genes showing evolutionary shifts have functions related to drought stress and flowering time, consistent with an adaptive response to selection. Loci differentiated between ancestors and descendants (F_ST outliers) were generally different from those showing signatures of selection based on site frequency spectrum analysis (Tajima's D), indicating that the loci that evolved in response to the recent drought and those under historical selection were generally distinct. Very few genes showed similar evolutionary responses between two geographically distinct populations, suggesting independent genetic trajectories of evolution yielding parallel phenotypic changes. The results show that selection can result in rapid genome‐wide evolutionary shifts in allele frequencies in natural populations, and highlight the usefulness of combining resurrection experiments in natural populations with genomics for studying the genetic basis of adaptive evolution.