PARALLEL EVOLUTION OF LAKE-STREAM PAIRS OF THREESPINE STICKLEBACKS (GASTEROSTEUS) INFERRED FROM MITOCHONDRIAL DNA VARIATION

Adding the third dimension to studies of parallel evolution of morphology and function: An exploration based on parapatric lake-stream stickleback

Recent methodological advances have led to a rapid expansion of evolutionary studies employing three-dimensional landmark-based geometric morphometrics (GM). GM methods generally enable researchers to capture and compare complex shape phenotypes, and to quantify their relationship to environmental gradients. However, some recent studies have shown that the common, inexpensive, and relatively rapid two-dimensional GM methods can distort important information and produce misleading results because they cannot capture variation in the depth (Z) dimension. We use micro-CT scanned threespine stickleback (Gasterosteus aculeatus Linnaeus, 1758) from six parapatric lake-stream populations on Vancouver Island, British Columbia, to test whether the loss of the depth dimension in 2D GM studies results in misleading interpretations of parallel evolution. Using joint locations described with 2D or 3D landmarks, we compare results from separate 2D and 3D shape spaces, from a combined 2D-3D shape space, and from estimates of biomechanical function. We show that, although shape is distorted enough in 2D projections to strongly influence the interpretation of morphological parallelism, estimates of biomechanical function are relatively robust to the loss of the Z dimension.

In the rivers: Multiple adaptive radiations of cyprinid fishes (Labeobarbus) in Ethiopian Highlands

Multiple repeated patterns of adaptive radiation were revealed in cyprinid fish inhabiting the compact geographic region of the Ethiopian Highlands. We found four independently evolved radiations in the evolutionary hexaploid (2n = 150) Labeobarbus lineage based on matrilineal relationships of >800 individuals. Each radiation displayed similar patterns of mouth phenotype diversification, and included ecomorphs/species of the generalized, lipped, scraping (one or two), and large-mouthed (one to three) types. All radiations were detected in geographically isolated rivers, and originated from different ancestral populations. This is the first documented case in which numerous parallel radiations of fishes occurred-via different ways-in a riverine environment. Some radiations are very recent and monophyletic, while others are older and include ecomorphs that originated in separate mini flocks and later combined into one. The diversification bursts among Ethiopian Labeobarbus were detected in the mid-upper reaches of rivers (1050-1550 m above sea level), which likely offer ecological opportunities that include diverse habitats yet poor fauna (i.e. lower competition and relaxed selection). This promising example of parallel evolution of adaptive radiation warrants further investigation.

Repeated Selection of Alternatively Adapted Haplotypes Creates Sweeping Genomic Remodeling in Stickleback

Heterogeneous genetic divergence can accumulate across the genome when populations adapt to different habitats while still exchanging alleles. How long does diversification take and how much of the genome is affected? When divergence occurs in parallel from standing genetic variation, how often are the same haplotypes involved? We explore these questions using restriction site-associated DNA sequencing genotyping data and show that broad-scale genomic repatterning, fueled by copious standing variation, can emerge in just dozens of generations in replicate natural populations of threespine stickleback fish (Gasterosteus aculeatus). After the catastrophic 1964 Alaskan earthquake, marine stickleback colonized newly created ponds on seismically uplifted islands. We find that freshwater fish in these young ponds differ from their marine ancestors across the same genomic segments previously shown to have diverged in much older lake populations. Outside of these core divergent regions the genome shows no population structure across the ocean-freshwater divide, consistent with strong local selection acting in alternative environments on stickleback populations still connected by significant gene flow. Reinforcing this inference, a majority of divergent haplotypes that are at high frequency in ponds are detectable in the sea, even across great geographic distances. Building upon previous population genomics work in this model species, our data suggest that a long history of divergent selection and gene flow among stickleback populations in oceanic and freshwater habitats has maintained polymorphisms of alternatively adapted DNA sequences that facilitate parallel evolution.

Parallel evolution of site-specific changes in divergent caribou lineages

The parallel evolution of phenotypes or traits within or between species provides important insight into the basic mechanisms of evolution. Genetic and genomic advances have allowed investigations into the genetic underpinnings of parallel evolution and the independent evolution of similar traits in sympatric species. Parallel evolution may best be exemplified among species where multiple genetic lineages, descended from a common ancestor, colonized analogous environmental niches, and converged on a genotypic or phenotypic trait. Modern North American caribou (Rangifer tarandus) originated from three ancestral sources separated during the Last Glacial Maximum (LGM): the Beringian-Eurasian lineage (BEL), the North American lineage (NAL), and the High Arctic lineage (HAL). Historical introgression between the NAL and the BEL has been found throughout Ontario and eastern Manitoba. In this study, we first characterized the functional differentiation in the cytochrome-b (cytB) gene by identifying nonsynonymous changes. Second, the caribou lineages were used as a direct means to assess site-specific parallel changes among lineages. There was greater functional diversity within the NAL despite the BEL having greater neutral diversity. The patterns of amino acid substitutions occurring within different lineages supported the parallel evolution of cytB amino acid substitutions suggesting different selective pressures among lineages. This study highlights the independent evolution of identical amino acid substitutions within a wide-ranging mammal species that have diversified from different ancestral haplogroups and where ecological niches can invoke parallel evolution.

Heritable gene expression differences between lake and stream stickleback include both parallel and antiparallel components

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.

Divergent Macroparasite Infections in Parapatric Swiss Lake-Stream Pairs of Threespine Stickleback (Gasterosteus aculeatus)

Spatial heterogeneity in diversity and intensity of parasitism is a typical feature of most host-parasite interactions, but understanding of the evolutionary implications of such variation is limited. One possible outcome of infection heterogeneities is parasite-mediated divergent selection between host populations, ecotypes or species which may facilitate the process of ecological speciation. However, very few studies have described infections in population-pairs along the speciation continuum from low to moderate or high degree of genetic differentiation that would address the possibility of parasite-mediated divergent selection in the early stages of the speciation process. Here we provide an example of divergent parasitism in freshwater fish ecotypes by examining macroparasite infections in threespine stickleback (Gasterosteus aculeatus) of four Swiss lake systems each harbouring parapatric lake-stream ecotype pairs. We demonstrate significant differences in infections within and between the pairs that are driven particularly by the parasite taxa transmitted to fish from benthic invertebrates. The magnitude of the differences tended to correlate positively with the extent of neutral genetic differentiation between the parapatric lake and stream populations of stickleback, whereas no such correlation was found among allopatric populations from similar or contrasting habitats. This suggests that genetic differentiation is unrelated to the magnitude of parasite infection contrasts when gene flow is constrained by geographical barriers while in the absence of physical barriers, genetic differentiation and the magnitude of differences in infections tend to be positively correlated.

Range dynamics, rather than convergent selection, explain the mosaic distribution of red-winged blackbird phenotypes

Geographic distributions of genetic and phenotypic characters can illuminate historical evolutionary processes. In particular, mosaic distributions of phenotypically similar populations can arise from parallel evolution or from irregular patterns of dispersal and colonization by divergent forms. Two phenotypically divergent forms of the red-winged blackbird (Agelaius phoeniceus) show a mosaic phenotypic distribution, with a "bicolored" form occurring disjunctly in California and Mexico. We analyzed the relationships among these bicolored populations and neighboring typical populations, using ∼600 bp of mitochondrial DNA sequence data and 10 nuclear short tandem repeat loci. We find that bicolored populations, although separated by ∼3000 km, are genetically more similar to one other than they are to typical populations separated by ∼400 km. We also find evidence of ongoing gene flow among populations, including some evidence of asymmetric gene flow. We conclude that the current distribution of bicolored forms represents incomplete speciation, where recent asymmetric hybridization with typical A. phoeniceus is dividing the range of a formerly widespread bicolored form. This hypothesis predicts that bicolored forms may suffer extinction by hybridization. Future work will use fine-scaled geographical sampling and nuclear sequence data to test for hybrid origins of currently typical populations and to more precisely quantify the directionality of gene flow.

Repeated lake-stream divergence in stickleback life history within a Central European lake basin

Life history divergence between populations inhabiting ecologically distinct habitats might be a potent source of reproductive isolation, but has received little attention in the context of speciation. We here test for life history divergence between threespine stickleback inhabiting Lake Constance (Central Europe) and multiple tributary streams. Otolith analysis shows that lake fish generally reproduce at two years of age, while their conspecifics in all streams have shifted to a primarily annual life cycle. This divergence is paralleled by a striking and consistent reduction in body size and fecundity in stream fish relative to lake fish. Stomach content analysis suggests that life history divergence might reflect a genetic or plastic response to pelagic versus benthic foraging modes in the lake and the streams. Microsatellite and mitochondrial markers further reveal that life history shifts in the different streams have occurred independently following the colonization by Lake Constance stickleback, and indicate the presence of strong barriers to gene flow across at least some of the lake-stream habitat transitions. Given that body size is known to strongly influence stickleback mating behavior, these barriers might well be related to life history divergence.

Population genetic dynamics of three-spined sticklebacks (Gasterosteus aculeatus) in anthropogenic altered habitats

In industrialized and/or agriculturally used landscapes, inhabiting species are exposed to a variety of anthropogenic changes in their environments. Genetic diversity may be reduced if populations encounter founder events, bottlenecks, or isolation. Conversely, genetic diversity may increase if populations adapt to changes in selective regimes in newly created habitats. With the present study, genetic variability of 918 sticklebacks from 43 samplings (21.3 ± 3.8 per sample) at 36 locations from cultivated landscapes in Northwest Germany was analyzed at nine neutral microsatellite loci. To test if differentiation is influenced by habitat alterations, sticklebacks were collected from ancient running waters and adjacent artificial stagnant waters, from brooks with salt water inflow of anthropogenic and natural origin and adjacent freshwater sites. Overall population structure was dominated by isolation by distance (IBD), which was significant across all populations, and analysis of molecular variance (AMOVA) revealed that 10.6% of the variation was explained by river catchment area. Populations in anthropogenic modified habitats deviated from the general IBD structure and in the AMOVA, grouping by habitat type running/stagnant water explained 4.9% of variation and 1.4% of the variation was explained by salt-/freshwater habitat. Sticklebacks in salt-polluted water systems seem to exhibit elevated migratory activity between fresh- and saltwater habitats, reducing IBD. In other situations, populations showed distinct signs of genetic isolation, which in some locations was attributed to mechanical migration barriers, but in others to potential anthropogenic induced bottleneck or founder effects. The present study shows that anthropogenic habitat alterations may have diverse effects on the population genetic structure of inhabiting species. Depending on the type of habitat change, increased genetic differentiation, diversification, or isolation are possible consequences.

Independent axes of genetic variation and parallel evolutionary divergence of opercle bone shape in threespine stickleback

Evolution of similar phenotypes in independent populations is often taken as evidence of adaptation to the same fitness optimum. However, the genetic architecture of traits might cause evolution to proceed more often toward particular phenotypes, and less often toward others, independently of the adaptive value of the traits. Freshwater populations of Alaskan threespine stickleback have repeatedly evolved the same distinctive opercle shape after divergence from an oceanic ancestor. Here we demonstrate that this pattern of parallel evolution is widespread, distinguishing oceanic and freshwater populations across the Pacific Coast of North America and Iceland. We test whether this parallel evolution reflects genetic bias by estimating the additive genetic variance-covariance matrix (G) of opercle shape in an Alaskan oceanic (putative ancestral) population. We find significant additive genetic variance for opercle shape and that G has the potential to be biasing, because of the existence of regions of phenotypic space with low additive genetic variation. However, evolution did not occur along major eigenvectors of G, rather it occurred repeatedly in the same directions of high evolvability. We conclude that the parallel opercle evolution is most likely due to selection during adaptation to freshwater habitats, rather than due to biasing effects of opercle genetic architecture.

Parallel and nonparallel aspects of ecological, phenotypic, and genetic divergence across replicate population pairs of lake and stream stickleback

Parallel (or convergent) evolution provides strong evidence for a deterministic role of natural selection: similar phenotypes evolve when independent populations colonize similar environments. In reality, however, independent populations in similar environments always show some differences: some nonparallel evolution is present. It is therefore important to explicitly quantify the parallel and nonparallel aspects of trait variation, and to investigate the ecological and genetic explanations for each. We performed such an analysis for threespine stickleback (Gasterosteus aculeatus) populations inhabiting lake and stream habitats in six independent watersheds. Morphological traits differed in the degree to which lake-stream divergence was parallel across watersheds. Some aspects of this variation were correlated with ecological variables related to diet, presumably reflecting the strength and specifics of divergent selection. Furthermore, a genetic scan revealed some markers that diverged between lakes and streams in many of the watersheds and some that diverged in only a few watersheds. Moreover, some of the lake-stream divergence in genetic markers was associated within some of the lake-stream divergence in morphological traits. Our results suggest that parallel evolution, and deviations from it, are primarily the result of natural selection, which corresponds in only some respects to the dichotomous habitat classifications frequently used in such studies.

Population genomics of parallel adaptation in threespine stickleback using sequenced RAD tags

Next-generation sequencing technology provides novel opportunities for gathering genome-scale sequence data in natural populations, laying the empirical foundation for the evolving field of population genomics. Here we conducted a genome scan of nucleotide diversity and differentiation in natural populations of threespine stickleback (Gasterosteus aculeatus). We used Illumina-sequenced RAD tags to identify and type over 45,000 single nucleotide polymorphisms (SNPs) in each of 100 individuals from two oceanic and three freshwater populations. Overall estimates of genetic diversity and differentiation among populations confirm the biogeographic hypothesis that large panmictic oceanic populations have repeatedly given rise to phenotypically divergent freshwater populations. Genomic regions exhibiting signatures of both balancing and divergent selection were remarkably consistent across multiple, independently derived populations, indicating that replicate parallel phenotypic evolution in stickleback may be occurring through extensive, parallel genetic evolution at a genome-wide scale. Some of these genomic regions co-localize with previously identified QTL for stickleback phenotypic variation identified using laboratory mapping crosses. In addition, we have identified several novel regions showing parallel differentiation across independent populations. Annotation of these regions revealed numerous genes that are candidates for stickleback phenotypic evolution and will form the basis of future genetic analyses in this and other organisms. This study represents the first high-density SNP-based genome scan of genetic diversity and differentiation for populations of threespine stickleback in the wild. These data illustrate the complementary nature of laboratory crosses and population genomic scans by confirming the adaptive significance of previously identified genomic regions, elucidating the particular evolutionary and demographic history of such regions in natural populations, and identifying new genomic regions and candidate genes of evolutionary significance.

Population genomics of parallel adaptation in threespine stickleback using sequenced RAD tags

Next-generation sequencing technology provides novel opportunities for gathering genome-scale sequence data in natural populations, laying the empirical foundation for the evolving field of population genomics. Here we conducted a genome scan of nucleotide diversity and differentiation in natural populations of threespine stickleback (Gasterosteus aculeatus). We used Illumina-sequenced RAD tags to identify and type over 45,000 single nucleotide polymorphisms (SNPs) in each of 100 individuals from two oceanic and three freshwater populations. Overall estimates of genetic diversity and differentiation among populations confirm the biogeographic hypothesis that large panmictic oceanic populations have repeatedly given rise to phenotypically divergent freshwater populations. Genomic regions exhibiting signatures of both balancing and divergent selection were remarkably consistent across multiple, independently derived populations, indicating that replicate parallel phenotypic evolution in stickleback may be occurring through extensive, parallel genetic evolution at a genome-wide scale. Some of these genomic regions co-localize with previously identified QTL for stickleback phenotypic variation identified using laboratory mapping crosses. In addition, we have identified several novel regions showing parallel differentiation across independent populations. Annotation of these regions revealed numerous genes that are candidates for stickleback phenotypic evolution and will form the basis of future genetic analyses in this and other organisms. This study represents the first high-density SNP-based genome scan of genetic diversity and differentiation for populations of threespine stickleback in the wild. These data illustrate the complementary nature of laboratory crosses and population genomic scans by confirming the adaptive significance of previously identified genomic regions, elucidating the particular evolutionary and demographic history of such regions in natural populations, and identifying new genomic regions and candidate genes of evolutionary significance.

Can gene flow have negative demographic consequences? Mixed evidence from stream threespine stickleback

Dispersal and gene flow can have both positive and negative effects on population size, but little empirical support from nature exists for the negative effects. We test for such effects in a stream population of threespine stickleback (Gasterosteus aculeatus L.) that is subject to high gene flow from a lake and is thus maladapted to stream conditions. In this system, maladaptation increases with distance along the stream, and this increase is associated with decreasing population densities until stickleback are no longer present (2.5 km from the lake). We conducted field experiments to inform whether this association might reflect a negative role for gene flow in constraining population size and therefore causing a local range limit. We specifically tested predictions deriving from theory: peripheral populations should show partial local adaptation, be under strong selection and not simply be maintained by dispersal. First, a transplant experiment suggested a weak home-site advantage in the peripheral population. Second, a mark-recapture study showed directional selection for a stream-adapted phenotype in 1 of 2 years. Third, another mark-recapture experiment showed that dispersal is limited to the point that positive demographic effects of dispersal are probably minimal. We conclude that, although gene flow does constrain morphological maladaptation in the outlet stream population, the evidence for its contribution to population size and range limits is mixed. We discuss the implications of our work for the study of factors influencing the evolution of species' ranges.

Historical contingency and ecological determinism interact to prime speciation in sticklebacks, Gasterosteus

Historical contingency and determinism are often cast as opposing paradigms under which evolutionary diversification operates. It may be, however, that both factors act together to promote evolutionary divergence, although there are few examples of such interaction in nature. We tested phylogenetic predictions of an explicit historical model of divergence (double invasions of freshwater by marine ancestors) in sympatric species of three-spined sticklebacks (Gasterosteus aculeatus) where determinism has been implicated as an important factor driving evolutionary novelty. Microsatellite DNA variation at six loci revealed relatively low genetic variation in freshwater populations, supporting the hypothesis that they were derived by colonization of freshwater by more diverse marine ancestors. Phylogenetic and genetic distance analyses suggested that pairs of sympatric species have evolved multiple times, further implicating determinism as a factor in speciation. Our data also supported predictions based on the hypothesis that the evolution of sympatric species was contingent upon 'double invasions' of postglacial lakes by ancestral marine sticklebacks. Sympatric sticklebacks, therefore, provide an example of adaptive radiation by determinism contingent upon historical conditions promoting unique ecological interactions, and illustrate how contingency and determinism may interact to generate geographical variation in species diversity