Parallel evolution of similar traits by independent populations in similar environments is considered strong evidence for adaptation by natural selection. Often, however, replicate populations in similar environments do not all evolve in the same way, thus deviating from any single, predominant outcome of evolution. This variation might arise from non-adaptive, population-specific effects of genetic drift, gene flow or limited genetic variation. Alternatively, these deviations from parallel evolution might also reflect predictable adaptation to cryptic environmental heterogeneity within discrete habitat categories. Here, we show that deviations from parallel evolution are the consequence of environmental variation within habitats combined with variation in gene flow. Threespine stickleback (Gasterosteus aculeatus) in adjoining lake and stream habitats (a lake-stream 'pair') diverge phenotypically, yet the direction and magnitude of this divergence is not always fully parallel among 16 replicate pairs. We found that the multivariate direction of lake-stream morphological divergence was less parallel between pairs whose environmental differences were less parallel. Thus, environmental heterogeneity among lake-stream pairs contributes to deviations from parallel evolution. Additionally, likely genomic targets of selection were more parallel between environmentally more similar pairs. In contrast, variation in the magnitude of lake-stream divergence (independent of direction) was better explained by differences in lake-stream gene flow; pairs with greater lake-stream gene flow were less morphologically diverged. Thus, both adaptive and non-adaptive processes work concurrently to generate a continuum of parallel evolution across lake-stream stickleback population pairs.

This study asks whether the spatial scale of sampling alters structural properties of food webs and whether any differences are attributable to changes in species richness and connectance with scale. Understanding how different aspects of sampling effort affect ecological network structure is important for both fundamental ecological knowledge and the application of network analysis in conservation and management. Using a highly resolved food web for the marine intertidal ecosystem of the Sanak Archipelago in the Eastern Aleutian Islands, Alaska, we assess how commonly studied properties of network structure differ for 281 versions of the food web sampled at five levels of spatial scale representing six orders of magnitude in area spread across the archipelago. Species (S) and link (L) richness both increased by approximately one order of magnitude across the five spatial scales. Links per species (L/S) more than doubled, while connectance (C) decreased by approximately two-thirds. Fourteen commonly studied properties of network structure varied systematically with spatial scale of sampling, some increasing and others decreasing. While ecological network properties varied systematically with sampling extent, analyses using the niche model and a power-law scaling relationship indicate that for many properties, this apparent sensitivity is attributable to the increasing S and decreasing C of webs with increasing spatial scale. As long as effects of S and C are accounted for, areal sampling bias does not have a special impact on our understanding of many aspects of network structure. However, attention does need be paid to some properties such as the fraction of species in loops, which increases more than expected with greater spatial scales of sampling.

The evolutionary consequences of temporal variation in selection remain hotly debated. We explored these consequences by studying threespine stickleback in a set of bar-built estuaries along the central California coast. In most years, heavy rains induce water flow strong enough to break through isolating sand bars, connecting streams to the ocean. New sand bars typically re-form within a few weeks or months, thereby re-isolating populations within the estuaries. These breaching events cause severe and often extremely rapid changes in abiotic and biotic conditions, including shifts in predator abundance. We investigated whether this strong temporal environmental variation can maintain within-population variation while eroding adaptive divergence among populations that would be caused by spatial variation in selection. We used neutral genetic markers to explore population structure and then analysed how stickleback armor traits, the associated genes Eda and Pitx1 and elemental composition (%P) varies within and among populations. Despite strong gene flow, we detected evidence for divergence in stickleback defensive traits and Eda genotypes associated with predation regime. However, this among-population variation was lower than that observed among other stickleback populations exposed to divergent predator regimes. In addition, within-population variation was very high as compared to populations from environmentally stable locations. Elemental composition was strongly associated with armor traits, Eda genotype and the presence of predators, thus suggesting that spatiotemporal variation in armor traits generates corresponding variation in elemental phenotypes. We conclude that gene flow, and especially temporal environmental variation, can maintain high levels of within-population variation while reducing, but not eliminating, among-population variation driven by spatial environmental variation.

The repeated phenotypic patterns that characterize populations undergoing parallel evolution provide support for a deterministic role of adaptation by natural selection. Determining the level of parallelism also at the genetic level is thus central to our understanding of how natural selection works. Many studies have looked for repeated genomic patterns in natural populations, but work on gene expression is less common. The studies that have examined gene expression have found some support for parallelism, but those studies almost always used samples collected from the wild that potentially confounds the effects of plasticity with heritable differences. Here we use two independent pairs of lake and stream threespine stickleback (Gasterosteus aculeatus) raised in common garden conditions to assess both parallel and antiparallel (that is, similar versus different directions of lake–stream expression divergence in the two watersheds) heritable gene expression differences as measured by total RNA sequencing. We find that more genes than expected by chance show either parallel (22 genes, 0.18% of expressed genes) or antiparallel (24 genes, 0.20% of expressed genes) lake–stream expression differences. These results correspond well with previous genomic studies in stickleback ecotype pairs that found similar levels of parallelism. We suggest that parallelism might be similarly constrained at the genomic and transcriptomic levels.

A population model was constructed to simulate the development of imazamox-resistant jointed goatgrass (AEGCY) in imazamox-resistant (Clearfield™) wheat. The model computed changes in the surface and in the buried AEGCY seed banks for both resistant and susceptible biotypes. Simulations started with an initial density of 1,000 susceptible and zero resistant seeds/m2 in each seed bank. Simulation of continuous, no-till Clearfield wheat production resulted in rapid development of resistant AEGCY without hybridization with wheat and in extremely rapid resistance development with hybridization. In less than 10 yr, the resistant population was growing exponentially in both simulations. Adding a fallow year with tillage into the simulated rotation did not substantially slow down the appearance of resistance but did delay the rate of resistant population increase by several orders of magnitude over 10 yr. Alternating Clearfield and a nonresistant winter wheat in combination with fallowing prevented the establishment of a significant resistant AEGCY population and prevented the susceptible seed population from increasing exponentially. These projections suggest that imazamox-resistant wheat can be a tool for managing AEGCY populations especially if combined with rotations that include fallow and crops other than Clearfield winter wheat.

No abstract

Haminoea japonica is a species of opisthobranch sea slug native to Japan and Korea. Non-native populations have spread unnoticed for decades due to difficulties in the taxonomy of Haminoea species. Haminoea japonica is associated with a schistosome parasite in San Francisco Bay, thus further spread could have consequence to human health and economies. Anecdotal evidence suggests that H. japonica has displaced native species of Haminoea in North America and Europe, becoming locally dominant in estuaries and coastal lagoons. In this paper we study the population genetics of native and non-native populations of H. japonica based on mt-DNA data including newly discovered populations in Italy and France. The conclusions of this study further corroborate a Northeastern Japan origin for the non-native populations and suggest possible independent introductions into North America and Europe. Additionally, the data obtained revealed possible secondary introductions within Japan. Although non-native populations have experienced severe genetic bottlenecks they have colonized different regions with a broad range of water temperatures and other environmental conditions. The environmental tolerance of this species, along with its ability to become dominant in invaded areas and its association with a schistosome parasite, suggest H. japonica could be a dangerous invasive species.

Populations receiving high maladaptive gene flow are expected to experience strong directional selection-because gene flow pulls mean phenotypes away from local fitness peaks. We tested this prediction by means of a large and replicated mark-recapture study of threespine stickleback (Gasterosteus aculeatus) in two stream populations. One of the populations (outlet) experiences high gene flow from the lake population and its morphology is correspondingly poorly adapted. The other population (inlet) experiences very low gene flow from the lake population and its morphology is correspondingly well adapted. Contrary to the above prediction, selection was not stronger in the outlet than in the inlet, a result that forced us to consider potential reasons for why maladaptive gene flow might not increase selection. Of particular interest, we show by means of a simple population genetic model that maladaptive gene flow can-under reasonable conditions-decrease the strength of directional selection. This outcome occurs when immigrants decrease mean fitness in the resident population, which decreases the strength of selection against maladapted phenotypes. We argue that this previously unrecognized effect of gene flow deserves further attention in theoretical and empirical studies.