Natural selection can drive the repeated evolution of reproductive isolation, but the genomic basis of parallel speciation remains poorly understood. We analyzed whole-genome divergence between replicate pairs of stick insect populations that are adapted to different host plants and undergoing parallel speciation. We found thousands of modest-sized genomic regions of accentuated divergence between populations, most of which are unique to individual population pairs. We also detected parallel genomic divergence across population pairs involving an excess of coding genes with specific molecular functions. Regions of parallel genomic divergence in nature exhibited exceptional allele frequency changes between hosts in a field transplant experiment. The results advance understanding of biological diversification by providing convergent observational and experimental evidence for selection's role in driving repeatable genomic divergence.

Predicting evolution remains difficult. We studied the evolution of cryptic body coloration and pattern in a stick insect using 25 years of field data, experiments, and genomics. We found that evolution is more difficult to predict when it involves a balance between multiple selective factors and uncertainty in environmental conditions than when it involves feedback loops that cause consistent back-and-forth fluctuations. Specifically, changes in color-morph frequencies are modestly predictable through time ( = 0.14) and driven by complex selective regimes and yearly fluctuations in climate. In contrast, temporal changes in pattern-morph frequencies are highly predictable due to negative frequency-dependent selection ( = 0.86). For both traits, however, natural selection drives evolution around a dynamic equilibrium, providing some predictability to the process.

potentially more complex, as the homogenising effects of gene flow must be countered [1][2][3] . The 49 genic model of speciation proposes that specific genetic regions subject to strong divergent 50

Understanding natural selection's effect on genetic variation is a major goal in biology, but the genome-scale consequences of contemporary selection are not well known. In a release and recapture field experiment we transplanted stick insects to native and novel host plants and directly measured allele frequency changes within a generation at 186 576 genetic loci. We observed substantial, genome-wide allele frequency changes during the experiment, most of which could be attributed to random mortality (genetic drift). However, we also documented that selection affected multiple genetic loci distributed across the genome, particularly in transplants to the novel host. Host-associated selection affecting the genome acted on both a known colour-pattern trait as well as other (unmeasured) phenotypes. We also found evidence that selection associated with elevation affected genome variation, although our experiment was not designed to test this. Our results illustrate how genomic data can identify previously underappreciated ecological sources and phenotypic targets of selection.

Diverse geographical modes and mechanisms of speciation are known, and individual speciation genes have now been identified. Despite this progress, genome-wide outcomes of different evolutionary processes during speciation are less understood. Here, we integrate ecological and spatial information, mating trials, transplantation data and analysis of 86 130 single nucleotide polymorphisms (SNPs) in eight populations (28 pairwise comparisons) of Timema cristinae stick insects to test the effects of different factors on genomic divergence in a system undergoing ecological speciation. We find patterns consistent with effects of numerous factors, including geographical distance, gene flow, divergence in host plant use and climate, and selection against maladaptive hybridization (i.e. reinforcement). For example, the number of highly differentiated ‘outlier loci’, allele-frequency clines and the overall distribution of genomic differentiation were recognizably affected by these factors. Although host use has strong effects on phenotypic divergence and reproductive isolation, its effects on genomic divergence were subtler and other factors had pronounced effects. The results demonstrate how genomic data can provide new insights into speciation and how genomic divergence can be complex, yet predictable. Future work could adopt experimental, mapping and functional approaches to directly test which genetic regions are affected by selection and determine their physical location in the genome.

The interplay between selection and aspects of the genetic architecture of traits (such as linkage, dominance, and epistasis) can either drive or constrain speciation [1-3]. Despite accumulating evidence that speciation can progress to "intermediate" stages-with populations evolving only partial reproductive isolation-studies describing selective mechanisms that impose constraints on speciation are more rare than those describing drivers. The stick insect Timema cristinae provides an example of a system in which partial reproductive isolation has evolved between populations adapted to different host plant environments, in part due to divergent selection acting on a pattern polymorphism [4, 5]. Here, we demonstrate how selection on a green/melanistic color polymorphism counteracts speciation in this system. Specifically, divergent selection between hosts does not occur on color phenotypes because melanistic T. cristinae are cryptic on the stems of both host species, are resistant to a fungal pathogen, and have a mating advantage. Using genetic crosses and genome-wide association mapping, we quantify the genetic architecture of both the pattern and color polymorphism, illustrating their simple genetic control. We use these empirical results to develop an individual-based model that shows how the melanistic phenotype acts as a "genetic bridge" that increases gene flow between populations living on different hosts. Our results demonstrate how variation in the nature of selection acting on traits, and aspects of trait genetic architecture, can impose constraints on both local adaptation and speciation.

The results demonstrate a pervasive effect of ongoing evolution in a spatial context on population and community dynamics. The eco-evolutionary model makes testable predictions about the influence of the spatial configuration of the patch network on metapopulation size and the spatial scale of adaptation.

Effects of plant traits on herbivore-carnivore interactions are well documented in component communities but are not well understood at the level of large, complex communities. We report on a 2-year field experiment testing mechanisms by which variation in food quality among eight temperate forest tree species alters avian suppression of an assemblage of dietary generalist caterpillars. Plant quality and bird effects varied dramatically among tree species; high-quality plants yielded herbivores of 50% greater mass than those on low-quality plants, and bird effects ranged from near 0% to 97% reductions in caterpillar density. We also find evidence for two mechanisms linking host plant quality to bird effects. If caterpillar density was statistically controlled for, birds had relatively strong effects on the herbivores of low-quality plants, as predicted by the slow-growth/high-mortality hypothesis. At the same time, caterpillar density increased with plant quality, and bird effects were density dependent. Consequently, the net effect of birds was strongest on the herbivores of high-quality plants, a dynamic we call the high-performance/high-mortality hypothesis. Host plant quality thus changes highly generalized herbivore-carnivore interactions by two complementary but opposing mechanisms. These results highlight the interrelatedness of plant-herbivore and herbivore-carnivore interactions and thus the importance of a tritrophic perspective.