Sexual isolation promotes divergence between parapatric lake and stream stickleback

Speciation can be initiated by adaptive divergence between populations in ecologically different habitats, but how sexually based reproductive barriers contribute to this process is less well understood. We here test for sexual isolation between ecotypes of threespine stickleback fish residing in adjacent lake and stream habitats in the Lake Constance basin, Central Europe. Mating trials exposing females to pairings of territorial lake and stream males in outdoor mesocosms allowing for natural reproductive behaviour reveal that mating occurs preferentially between partners of the same ecotype. Compared to random mating, this sexual barrier reduces gene flow between the ecotypes by some 36%. This relatively modest strength of sexual isolation is surprising because comparing the males between the two ecotypes shows striking differentiation in traits generally considered relevant to reproductive behaviour (body size, breeding coloration, nest size). Analysing size differences among the individuals in the mating trials further indicates that assortative mating is not related to ecotype differences in body size. Overall, we demonstrate that sexually based reproductive isolation promotes divergence in lake–stream stickleback along with other known reproductive barriers, but we also caution against inferring strong sexual isolation from the observation of strong population divergence in sexually relevant traits.     

Drivers of population genetic differentiation in the wild: isolation by dispersal limitation,isolation by adaptation and isolation by colonization

Empirical population genetic studies have been dominated by a neutralist view, according to which gene flow and drift are the main forces driving population genetic structure in nature. The neutralist view in essence describes a process of isolation by dispersal limitation (IBDL) that generally leads to a pattern of isolation by distance (IBD). Recently, however, conceptual frameworks have been put forward that view local genetic adaptation as an important driver of population genetic structure. Isolation by adaptation (IBA) and monopolization (M) posit that gene flow among natural populations is reduced as a consequence of local genetic adaptation. IBA stresses that effective gene flow is reduced among habitats that show dissimilar ecological characteristics, leading to a pattern of isolation by environment. In monopolization, local genetic adaptation of initial colonizing genotypes results in a reduction in gene flow that fosters the persistence of founder effects. Here, we relate these different processes driving landscape genetic structure to patterns of IBD and isolation by environment (IBE). We propose a method to detect whether IBDL, IBA and M shape genetic differentiation in natural landscapes by studying patterns of variation at neutral and non‐neutral markers as well as at ecologically relevant traits. Finally, we reinterpret a representative number of studies from the recent literature by associating patterns to processes and identify patterns associated with local genetic adaptation to be as common as IBDL in structuring regional genetic variation of populations in the wild. Our results point to the importance of quantifying environmental gradients and incorporating ecology in the analysis of population genetics.     

Isolation by environment

The interactions between organisms and their environments can shape distributions of spatial genetic variation, resulting in patterns of isolation by environment (IBE) in which genetic and environmental distances are positively correlated, independent of geographic distance. IBE represents one of the most important patterns that results from the ways in which landscape heterogeneity influences gene flow and population connectivity, but it has only recently been examined in studies of ecological and landscape genetics. Nevertheless, the study of IBE presents valuable opportunities to investigate how spatial heterogeneity in ecological processes, agents of selection and environmental variables contributes to genetic divergence in nature. New and increasingly sophisticated studies of IBE in natural systems are poised to make significant contributions to our understanding of the role of ecology in genetic divergence and of modes of differentiation both within and between species. Here, we describe the underlying ecological processes that can generate patterns of IBE, examine its implications for a wide variety of disciplines and outline several areas of future research that can answer pressing questions about the ecological basis of genetic diversity.     

The geography of divergence with gene flow facilitates multitrait adaptation and the evolution of pollinator isolation in Mimulus aurantiacus

Ecological adaptation is the driving force during divergence with gene flow and generates reproductive isolation early in speciation. Although gene flow opposes divergence, local adaptation can be facilitated by factors that prevent the breakup of favorable allelic combinations. We investigated how selection, genetic architecture, and geography have contributed to the maintenance of floral trait divergence and pollinator isolation between parapatric ecotypes of Mimulus aurantiacus. Combining greenhouse, field, and genomic studies, we show that sharp clines in floral traits are maintained by spatially varying selection. Although adaptation breaks down where the ecotypes co‐occur, leading to the formation of a hybrid zone, the largely non‐overlapping distributions of the ecotypes shield them from immigrant genes, facilitating divergence across most of the range. In contrast to the sharp genetic discontinuities observed across most hybrid zones, we observed a gradual cline in genome‐wide divergence and a pattern of isolation by distance across the landscape. Thus, contrary to a long period of allopatry followed by recent re‐contact, our data suggest that floral trait divergence in M. aurantiacus may have evolved with locally restricted, but ongoing gene flow. Therefore, our study reveals how the geographic distribution of an organism can contribute to the evolution of premating isolation in the early stages of divergence with gene flow.     

Quantifying the constraining influence of gene flow on adaptive divergence in the lake-stream threespine stickleback system

The constraining effect of gene flow on adaptive divergence is often inferred but rarely quantified. We illustrate ways of doing so using stream populations of threespine stickleback (Gasterosteus aculeatus) that experience different levels of gene flow from a parapatric lake population. In the Misty Lake watershed (British Columbia, Canada), the inlet stream population is morphologically divergent from the lake population, and presumably experiences little gene flow from the lake. The outlet stream population, however, shows an intermediate phenotype and may experience more gene flow from the lake. We first used microsatellite data to demonstrate that gene flow from the lake is low into the inlet but high into the outlet, and that gene flow from the lake remains relatively constant with distance along the outlet. We next combined gene flow data with morphological and habitat data to quantify the effect of gene flow on morphological divergence. In one approach, we assumed that inlet stickleback manifest well-adapted phenotypic trait values not constrained by gene flow. We then calculated the deviation between the observed and expected phenotypes for a given habitat in the outlet. In a second approach, we parameterized a quantitative genetic model of adaptive divergence. Both approaches suggest a large impact of gene flow, constraining adaptation by 80-86% in the outlet (i.e., only 14-20% of the expected morphological divergence in the absence of gene flow was observed). Such approaches may be useful in other taxa to estimate how important gene flow is in constraining adaptive divergence in nature.     

Quantifying the roles of ecology and geography in spatial genetic divergence

Investigating the properties of ecological landscapes that influence gene flow among populations can provide key insights into the earliest stages of biological divergence. Both ecological and geographical factors can reduce gene flow, which can lead to population divergence, but we know little of the relative strengths of these phenomena in nature. Here, we use a novel application of structural equation modelling to quantify the contributions of ecological and geographical isolation to spatial genetic divergence in 17 species of Anolis lizards. Our comparative analysis shows that although both processes contributed significantly, geographical isolation explained substantially more genetic divergence than ecological isolation (36.3 vs. 17.9% of variance respectively), suggesting that despite the proposed ubiquity of ecological divergence, non‐ecological factors play the dominant role in the evolution of spatial genetic divergence.     

Rapid sexual and genomic isolation in sympatric Drosophila without reproductive character displacement

The rapid evolution of sexual isolation in sympatry has long been associated with reinforcement (i.e., selection to avoid maladaptive hybridization). However, there are many species pairs in sympatry that have evolved rapid sexual isolation without known costs to hybridization. A major unresolved question is what evolutionary processes are involved in driving rapid speciation in such cases. Here, we focus on one such system; the Drosophila athabasca species complex, which is composed of three partially sympatric and interfertile semispecies: WN, EA, and EB. To study speciation in this species complex, we assayed sexual and genomic isolation within and between these semispecies in both sympatric and allopatric populations. First, we found no evidence of reproductive character displacement (RCD) in sympatric zones compared to distant allopatry. Instead, semispecies were virtually completely sexually isolated from each other across their entire ranges. Moreover, using spatial approaches and coalescent demographic simulations, we detected either zero or only weak heterospecific gene flow in sympatry. In contrast, within each semispecies we found only random mating and little population genetic structure, except between highly geographically distant populations. Finally, we determined that speciation in this system is at least an order of magnitude older than previously assumed, with WN diverging first, around 200K years ago, and EA and EB diverging 100K years ago. In total, these results suggest that these semispecies should be given full species status and we adopt new nomenclature: WN—D. athabasca, EA—D. mahican, and EB—D. lenape. While the lack of RCD in sympatry and interfertility do not support reinforcement, we discuss what additional evidence is needed to further decipher the mechanisms that caused rapid speciation in this species complex.     

Rapid ecological speciation in three‐spined stickleback Gasterosteus aculeatus from Middleton Island,Alaska: the roles of selection and geographic isolation

Morphological analysis of three‐spined stickleback Gasterosteus aculeatus collected in Middleton Island, Alaska, was conducted in order to study how gene flow and selection interact during divergence. Middleton Island was uplifted by 3·4 m during the Great Alaska Earthquake in 1964; this event formed a series of new freshwater sites, triggering rapid evolution, and probably rapid speciation, in G. aculeatus populations that colonized them. The level of hybridization between the anadromous and the resident freshwater populations is reflected by the level of morphological variance of the resident freshwater G. aculeatus. Therefore, geographic isolation of the sites from the sea (approximating gene flow) and ionic concentration of the water (reflecting selection pressures) were correlated with morphological variance of the resident freshwater populations. Geographic isolation was negatively correlated with morphological variance in a majority of the analysed traits. Both selection and gene flow surrogates were found to be important influences on variance in morphology, though selection had a larger effect, especially on armour traits. It was concluded that gene flow appeared to constrain ecological speciation, but even in the presence of gene flow the strong selection in the freshwater environment was apparently leading to rapid divergence.     

Rapid genetic divergence in postglacial populations of threespine stickleback (Gasterosteus aculeatus): the role of habitat type, drainage and geographical proximity

The goal of this study was to evaluate the role of common ancestry, and of geographical or reproductive isolation, in genetic divergence in populations of threespine sticklebacks (Gasterosteus aculeatus). Using seven DNA microsatellite loci we compared the effects of habitat type, drainage system and geographical proximity on genetic distance among 16 populations situated in an area in Schleswig-Holstein (Germany) that became deglaciated approximately 12 000 years ago. Stickleback population structure correlated only weakly with drainage system, whereas the primary divergence was among habitat types. Phylogenetic analysis revealed that lake (n = 7) and river (n = 5) populations formed two distinct clades (Cavalli-Sforza's and Edwards' chord distance, 82-100% bootstrap support) at approximately equal genetic distances to a third clade, comprising putative estuarine (n = 4) ancestors. Allele frequencies in lake and river populations represented different subsets of the genetically more diverse estuarine populations. In nested amovas approximately twice the genetic variance was distributed among lake vs. river vs. estuarine populations as compared with the combined effects of drainage system and geographical distance. Limited gene flow between habitat types must have been established after postglacial colonization, suggesting ecological hybrid inferiority or behavioural mating barriers between ecotypes. Within estuarine and lake populations, population differentiation followed an isolation-by-distance model. Given the high observed heterozygosities within the 16 study populations (HO = 0.65-0.87), the mean divergence between lake and river population pairs (FST = 0.18 +/- 0.007) would be reached after 300-6000 generations in a stepwise mutation model, depending on the size of N(e). This demonstrates both the utility of hypervariable microsatellites for detecting recent population divergences and the danger of operating at temporal or spatial scales which are beyond their resolution.     

Ecological divergence and habitat isolation between two migratory forms of Japanese threespine stickleback (Gasterosteus aculeatus)

When two closely related species migrate to divergent spawning sites, divergent use of spawning habitats can directly reduce heterospecific mating. Furthermore, adaptations to divergent spawning habitats can promote speciation as a by‐product of ecological divergence. Here, we investigated habitat isolation and ecological divergence between two anadromous forms of threespine stickleback (Gasterosteus aculeatus), the Japan Sea and Pacific Ocean forms. In several coastal regions of eastern Hokkaido, Japan, these forms migrate to the same watershed to spawn. Our field surveys in a single watershed revealed that segregation of distinct spawning sites between the two forms was maintained within the watershed across multiple years. These spawning sites diverged in salinity and predator composition. Morphological and physiological divergence between the forms also occurs in the direction predicted by ecological differences between the spawning sites. Our data indicate that migration into divergent spawning habitats can be an important mechanism contributing to speciation and phenotypic divergence in anadromous fishes.     

The influence of landscape configuration and environment on population genetic structure in a sedentary passerine: insights from loci located in different genomic regions

The study of the factors structuring genetic variation can help to infer the neutral and adaptive processes shaping the demographic and evolutionary trajectories of natural populations. Here, we analyse the role of isolation by distance (IBD), isolation by resistance (IBR, defined by landscape composition) and isolation by environment (IBE, estimated as habitat and elevation dissimilarity) in structuring genetic variation in 25 blue tit (Cyanistes caeruleus) populations. We typed 1385 individuals at 26 microsatellite loci classified into two groups by considering whether they are located into genomic regions that are actively (TL; 12 loci) or not (NTL; 14 loci) transcribed to RNA. Population genetic differentiation was mostly detected using the panel of NTL. Landscape genetic analyses showed a pattern of IBD for all loci and the panel of NTL, but genetic differentiation estimated at TL was only explained by IBR models considering high resistance for natural vegetation and low resistance for agricultural lands. Finally, the absence for IBE suggests a lack of divergent selection pressures associated with differences in habitat and elevation. Overall, our study shows that markers located in different genomic regions can yield contrasting inferences on landscape‐level patterns of realized gene flow in natural populations.     

Effects of genetics and light environment on colour expression in threespine sticklebacks

The genetic basis of traits that are under sexual selection and that are involved in recognizing conspecific mates is poorly known, even in systems in which the phenotypic basis of these traits has been well studied. In the present study, we investigate genetic and environmental influences on nuptial colour, which plays important roles in sexual selection and sexual isolation in species pairs of limnetic and benthic threespine sticklebacks ( Gasterosteus aculeatus species complex). Previous work demonstrated that colour differences among species correlate to differences in the ambient light prevalent in their mating habitat. Red fish are found in clear water and black fish in red-shifted habitats. We used a paternal half-sib split-clutch design to investigate the genetic and environmental basis of nuptial colour. We found genetic differences between a red and a black stickleback population in the expression of both red and black nuptial colour. In addition, the light environment influenced colour expression, and genotype by environment interactions were also present. We found evidence for both phenotypic and genetic correlations between our colour traits; some of these correlations are in opposite directions for our red and black populations. These results suggest that both genetic change and phenotypic plasticity underlie the correlation of male colour with light environment.  © 2008 The Linnean Society of London, Biological Journal of the Linnean Society , 2008, 94 , 663–673.     

Historical divergence of mechanical isolation agents in the ground beetle Carabus arrowianus as revealed by phylogeographical analyses

In the carabid genus Carabus subgenus Ohomopterus , diverged body size and genital morphology serve as mechanical reproductive barriers. To elucidate the diverging process of body and genital sizes in Carabus arrowianus , which exhibits marked morphological diversity among geographical populations and may represent an early stage of speciation, we analysed a mitochondrial gene sequence for 1051 individuals from 63 populations and male morphology for 359 individuals from 47 populations. Two discrete morphological groups segregated by geographical barriers were distinguished, one of which possessed smaller bodies and shorter genitalia ( S group) than the other ( L group), which exhibited larger bodies and exaggerated genitalia. Genetic divergence between the two groups was significant but not large. Phylogeographical and population genetic analyses indicated that the L group was derived from the S group, and a coalescent simulation revealed that the two groups diverged during the latest middle Pleistocene (0.13 million years ago), with a much larger effective population size in the L group than the S group. Because the body size divergence could not be explained by adaptation to climatic conditions and genital morphology is considered to be subject to sexual selection, we postulated that a population division and colonization in favourable habitats caused by the Pleistocene climatic and geographical change might facilitate natural and sexual selection for enlarged body and genital sizes in the L group.     

Isolation by distance in the Atlantic cod, Gadus morhua, at large and small geographic scales

Genetic isolation by distance (IBD) has rarely been described in marine species with high potential for dispersal at both the larval and adult life-history stages. Here, we report significant relationships between inferred levels of gene flow and geographic distance in the Atlantic cod, Gadus morhua, at 10 nuclear restriction-fragment-length-polymorphism (RFLP) loci at small regional scales in the western north Atlantic region (< 1,600 km) that mirror those previously detected over its entire geographic range (up to 7,300 km). Highly significant allele frequency differences were observed among eight northwestern Atlantic populations, although the mean FST for all 10 loci was only 0.014. Despite this weak population structuring, the distance separating populations explained between 54% and 62% of the variation in gene flow depending on whether nine or 10 loci were used to estimate Nm. Across the species' entire geographic range, highly significant differences were observed among six regional populations at nine of the 10 loci (mean FST = 0.068) and seven loci exhibited significant negative relationships between gene flow and distance. At this large geographic scale, natural selection acting in the vicinity of one RFLP locus (GM798) had a significant effect on the correlation between gene flow and distance, and eliminating it from the analysis caused the coefficient of determination to increase from 17% to 62%. The role of vicariance was assessed by sequentially removing populations from the analysis and was found to play a minor role in contributing to the relationship between gene flow and distance at either geographic scale. The correlation between gene flow and distance detected in G. morhua at small and large spatial scales suggests that dispersal distances and effective population sizes are much smaller than predicted for the species and that the recent age of populations, rather than extensive gene flow, may be responsible for its weak population structure. Our results suggest that interpreting limited genetic differences among populations as reflecting high levels of ongoing gene flow should be made with caution.