Genomic signatures of relaxed disruptive selection associated with speciation reversal in whitefish

Background: Speciation reversal: the erosion of species differentiation via an increase in introgressive hybridization due to the weakening of previously divergent selection regimes, is thought to be an important, yet poorly understood, driver of biodiversity loss. Our study system, the Alpine whitefish (Coregonus spp.) species complex is a classic example of a recent postglacial adaptive radiation: forming an array of endemic lake flocks, with the independent origination of similar ecotypes among flocks. However, many of the lakes of the Alpine radiation have been seriously impacted by anthropogenic nutrient enrichment, resulting in a collapse in neutral genetic and phenotypic differentiation within the most polluted lakes. Here we investigate the effects of eutrophication on the selective forces that have shaped this radiation, using population genomics. We studied eight sympatric species assemblages belonging to five independent parallel adaptive radiations, and one species pair in secondary contact. We used AFLP markers, and applied FST outlier (BAYESCAN, DFDIST) and logistic regression analyses (MATSAM), to identify candidate regions for disruptive selection in the genome and their associations with adaptive traits within each lake flock. The number of outlier and adaptive trait associated loci identified per lake were then regressed against two variables (historical phosphorus concentration and contemporary oxygen concentration) representing the strength of eutrophication. Results: Whilst we identify disruptive selection candidate regions in all lake flocks, we find similar trends, across analysis methods, towards fewer disruptive selection candidate regions and fewer adaptive trait/candidate loci associations in the more polluted lakes. Conclusions: Weakened disruptive selection and a concomitant breakdown in reproductive isolating mechanisms in more polluted lakes has lead to increased gene flow between coexisting Alpine whitefish species. We hypothesize that the resulting higher rates of interspecific recombination reduce either the number or extent of genomic islands of divergence surrounding loci evolving under disruptive natural selection. This produces the negative trend seen in the number of selection candidate loci recovered during genome scans of whitefish species flocks, with increasing levels of anthropogenic eutrophication: as the likelihood decreases that AFLP restriction sites will fall within regions of heightened genomic divergence and therefore be classified as FST outlier loci. This study explores for the first time the potential effects of human-mediated relaxation of disruptive selection on heterogeneous genomic divergence between coexisting species.

Widespread phenotypic and genetic divergence along altitudinal gradients in animals

Altitudinal gradients offer valuable study systems to investigate how adap- tive genetic diversity is distributed within and between natural populations and which factors promote or prevent adaptive differentiation. The environ- mental clines along altitudinal gradients tend to be steep relative to the dispersal distance of many organisms, providing an opportunity to study the joint effects of divergent natural selection and gene flow. Temperature is one variable showing consistent altitudinal changes, and altitudinal gradi- ents can therefore provide spatial surrogates for some of the changes antici- pated under climate change. Here, we investigate the extent and patterns of adaptive divergence in animal populations along altitudinal gradients by sur- veying the literature for (i) studies on phenotypic variation assessed under common garden or reciprocal transplant designs and (ii) studies looking for signatures of divergent selection at the molecular level. Phenotypic data show that significant between-population differences are common and taxo- nomically widespread, involving traits such as mass, wing size, tolerance to thermal extremes and melanization. Several lines of evidence suggest that some of the observed differences are adaptively relevant, but rigorous tests of local adaptation or the link between specific phenotypes and fitness are sorely lacking. Evidence for a role of altitudinal adaptation also exists for a number of candidate genes, most prominently haemoglobin, and for anony- mous molecular markers. Novel genomic approaches may provide valuable tools for studying adaptive diversity, also in species that are not amenable to experimentation.

Repeated and predictable patterns of ecotypic differentiation during a biological invasion: lake-stream divergence in parapatric Swiss stickleback

The relative importance of ecological selection and geographical isolation in promoting and constraining genetic and phenotypic differentiation among populations is not always obvious. Interacting with divergent selection, restricted opportunity for gene flow may in some cases be as much a cause as a consequence of adaptation, with the latter being a hallmark of ecologi- cal speciation. Ecological speciation is well studied in parts of the native range of the three-spined stickleback. Here, we study this process in a recently invaded part of its range. Switzerland was colonized within the past 140 years from at least three different colonization events involving differ- ent stickleback lineages. They now occupy diverse habitats, ranging from small streams to the pelagic zone of large lakes. We use replicated systems of parapatric lake and stream populations, some of which trace their origins to different invasive lineages, to ask (i) whether phenotypic divergence occurred among populations inhabiting distinct habitats, (ii) whether trajec- tories of phenotypic divergence follow predictable parallel patterns and (iii) whether gene flow constrains divergent adaptation or vice versa. We find consistent phenotypic divergence between populations occupying distinct habitats. This involves parallel evolution in several traits with known eco- logical relevance in independent evolutionary lineages. Adaptive divergence supersedes homogenizing gene flow even at a small spatial scale. We find evidence that adaptive phenotypic divergence places constraints on gene flow over and above that imposed by geographical distance, signalling the early onset of ecological speciation.

A test of genetic association among male nuptial coloration, female mating preference, and male aggression bias within a polymorphic population of cichlid fish

Both inter- and intrasexual selection have been implicated in the origin and maintenance of species-rich taxa with diverse sexual traits. Simultaneous disruptive selection by female mate choice and male-male competition can, in theory, lead to speciation without geographical isolation if both act on the same male trait. Female mate choice can generate discontinuities in gene flow, while male-male competition can generate negative frequency-dependent selection stabilizing the male trait polymorphism. Speciation may be facilitated when mating preference and/or aggression bias are physically linked to the trait they operate on. We tested for genetic associations among female mating preference, male aggression bias and male coloration in the Lake Victoria cichlid Pundamilia. We crossed females from a phenotypically variable population with males from both extreme ends of the phenotype distribution in the same population (blue or red). Male offspring of a red sire were significantly redder than males of a blue sire, indicating that intra-population variation in male coloration is heritable. We tested mating preferences of female offspring and aggression biases of male offspring using binary choice tests. There was no evidence for associations at the family level between female mating preferences and coloration of sires, but dam identity had a significant effect on female mate preference. Sons of the red sire directed significantly more aggression to red than blue males, whereas sons of the blue sire did not show any bias. There was a positive correlation among individuals between male aggression bias and body coloration, possibly due to pleiotropy or physical linkage, which could facilitate the maintenance of color polymorphism.

Hidden Biodiversity in an Ecologically Important Freshwater Amphipod: Differences in Genetic Structure between Two Cryptic Species

Cryptic species, i.e. species that are morphologically hard to distinguish, have been detected repeatedly in various taxa and ecosystems. In order to evaluate the importance of this finding, we have to know in how far cryptic species differ in various aspects of their biology. The amphipod Gammarus fossarum is a key invertebrate in freshwater streams and contains several cryptic species. We examined the population genetic structure, genetic diversity and demographic history of two of them (type A and type B) using microsatellite markers and asked whether they show significant differences. We present results of population genetic analyses based on a total of 37 populations from the headwaters of two major European drainages, Rhine and Rhone. We found that, in both species, genetic diversity was geographically structured among and within drainages. For type A in the Rhine and type B in the Rhone, we detected significant patterns of isolation by distance. The increase of genetic differentiation with geographical distance, however, was much higher in type A than in type B. This result indicates substantial interspecific differences in population history and/or the extent of current gene flow between populations. In the Rhine, type B does not show evidence of isolation by distance, and population differentiation is relatively low across hundreds of kilometres. The majority of these populations also show signatures of recent bottlenecks. These patterns are consistent with a recent expansion of type B into the Rhine drainage. In summary, our results suggest considerable and previously unrecognized interspecific differences in the genetic structure of these cryptic keystone species.

Phylogeography of silver pheasant (Lophura nycthemera L.) across China: aggregate effects of refugia, introgression, and riverine barriers

The role of Pleistocene glacial cycles in forming the contemporary genetic structure of organisms has been well studied in China with a particular focus on the Tibetan Plateau. However, China has a complex topography and diversity of local climates, and how glacial cycles may have shaped the subtropical and tropical biota of the region remains mostly unaddressed. To investigate the factors that affected the phylogeography and population history of a widely distributed and nondeciduous forest species, we analysed morphological characters, mitochondrial DNA sequences and nuclear microsatellite loci in the Silver Pheasant (Lophura nycthemera). In a pattern generally consistent with phenotypic clusters, but not nominal subspecies, deeply divergent mitochondrial lineages restricted to different geographic regions were detected. Coalescent simulations indicated that the time of main divergence events corresponded to major glacial periods in the Pleistocene and gene flow was only partially lowered by drainage barriers between some populations. Intraspecific cytonuclear discordance was revealed in mitochondrial lineages from Hainan Island and the Sichuan Basin with evidence of nuclear gene flow from neighbouring populations into the latter. Unexpectedly, hybridization was revealed in Yingjiang between the Silver Pheasant and Kalij Pheasant (Lophura leucomelanos) with wide genetic introgression at both the mtDNA and nuclear levels. Our results highlight a novel phylogeographic pattern in a subtropical area generated from the combined effects of climate oscillation, partial drainage barriers and interspecific hybridization. Cytonuclear discordance combined with morphological differentiation implies that complex historical factors shaped the divergence process in this biodiversity hot spot area of southern China.

Sex-specific clines support incipient speciation in a common European mammal

Hybrid zones provide excellent opportunities to study processes and mechanisms underlying reproductive isolation and speciation. Here we investigated sex-specific clines of molecular markers in hybrid zones of morphologically cryptic yet genetically highly-diverged evolutionary lineages of the European common vole (Microtus arvalis). We analyzed the position and width of four secondary contact zones along three independent transects in the region of the Alps using maternally (mitochondrial DNA) and paternally (Y-chromosome) inherited genetic markers. Given male-biased dispersal in the common vole, a selectively neutral secondary contact would show broader paternal marker clines than maternal ones. In a selective case, for example, involving a form of Haldane’s rule, Y-chromosomal clines would not be expected to be broader than maternal markers because they are transmitted by the heterogametic sex and thus gene flow would be restricted. Consistent with the selective case, paternal clines were significantly narrower or at most equal in width to maternal clines in all contact zones. In addition, analyses using maximum likelihood cline-fitting detected a shift of paternal relative to maternal clines in three of four contact zones. These patterns suggest that processes at the contact zones in the common vole are not selectively neutral, and that partial reproductive isolation is already established between these evolutionary lineages. We conclude that hybrid zone movement, sexual selection and/or genetic incompatibilities are likely associated with an unusual unidirectional manifestation of Haldane’s rule in this common European mammal.

Female preference for male color is necessary and sufficient for assortative mating in 2 cichlid sister species

A critical step for speciation in the face of gene flow is the origination of reproductive isolation. The evolution of assortative mating greatly facilitates this process. Assortative mating can be mediated by one or multiple cues across an array of sensory modalities. We here explore possible cues that may underlie female mate choice in a sympatric species pair of cichlid fish from Lake Victoria, Pundamilia pundamilia and Pundamilia nyererei. Previous studies identified species-specific female preferences for male coloration, but effects of other cues could not be ruled out. Therefore, we assessed female choice in a series of experiments in which we manipulated visual (color) and chemical cues. We show that the visibility of differences in nuptial hue (i.e., either blue or red) between males of the 2 species is necessary and sufficient for assortative mating by female mate choice. Such assortment mediated by a single cue may evolve relatively quickly, but could make reproductive isolation vulnerable to environmental changes. These findings confirm the important role of female mate choice for male nuptial hue in promoting the explosive speciation of African haplochromine cichlids.

Recent speciation between sympatric Tanganyikan cichlid colour morphs

Lake Tanganyika, Africa’s oldest lake, harbours an impressive diversity of cichlid fishes. Although diversification in its radiating groups is thought to have been initially rapid, cichlids from Lake Tanganyika show little evidence for ongoing speciation. In contrast, examples of recent divergence among sympatric colour morphs are well known in haplochromine cichlids from Lakes Malawi and Victoria. Here, we report genetic evidence for recent divergence between two sympatric Tanganyikan cichlid colour morphs. These Petrochromis morphs share mitochondrial haplotypes, yet microsatellite loci reveal that their sympatric populations form distinct genetic groups. Nuclear divergence between the two morphs is equivalent to that which arises geographically within one of the morphs over short distances and is substantially smaller than that among other sympatric species in this genus. These patterns suggest that these morphs diverged only recently, yet that barriers to gene flow exist which prevent extensive admixture despite their sympatric distribution. The morphs studied here provide an unusual example of active diversification in Lake Tanganyika’s generally ancient cichlid fauna and enable comparisons of speciation processes between Lake Tanganyika and other African lakes.

Tree growth response along an elevational gradient: climate or genetics?

Environment and genetics combine to influence tree growth and should therefore be jointly considered when evaluating forest responses in a warming climate. Here, we combine dendroclimatology and population genetic approaches with the aim of attributing climatic influences on growth of European larch (Larix decidua) and Norway spruce (Picea abies). Increment cores and genomic DNA samples were collected from populations along a ~900-m elevational transect where the air temperature gradient encompasses a ~4 °C temperature difference. We found that low genetic differentiation among populations indicates gene flow is high, suggesting that migration rate is high enough to counteract the selective pressures of local environmental variation. We observed lower growth rates towards higher elevations and a transition from negative to positive correlations with growing season temperature upward along the elevational transect. With increasing elevation there was also a clear increase in the explained variance of growth due to summer temperatures. Comparisons between climate sensitivity patterns observed along this elevational transect with those from Larix and Picea sites distributed across the Alps reveal good agreement, and suggest that tree-ring width (TRW) variations are more climate-driven than genetics-driven at regional and larger scales. We conclude that elevational transects are an extremely valuable platform for understanding climatic-driven changes over time and can be especially powerful when working within an assessed genetic framework.

Hybridization and adaptive radiation

Whether interspecific hybridization is important as a mechanism that generates biological diversity is a matter of controversy. Whereas some authors focus on the potential of hybridization as a source of genetic variation, functional novelty and new species, others argue against any important role, because reduced fitness would typically render hybrids an evolutionary dead end. By drawing on recent developments in the genetics and ecology of hybridization and on principles of ecological speciation theory, I develop a concept that reconciles these views and adds a new twist to this debate. Because hybridization is common when populations invade new environments and potentially elevates rates of response to selection, it predisposes colonizing populations to rapid adaptive diversification under disruptive or divergent selection. I discuss predictions and suggest tests of this hybrid swarm theory of adaptive radiation and review published molecular phylogenies of adaptive radiations in light of the theory. Some of the confusion about the role of hybridization in evolutionary diversification stems from the contradiction between a perceived necessity for cessation of gene flow to enable adaptive population differentiation on the one hand [1], and the potential of hybridization for generating adaptive variation, functional novelty and new species 2, 3 and 4 on the other. Much progress in the genetics 5, 6, 7, 8 and 9 and ecology of hybridization 9, 10 and 11, and in our understanding of the role of ecology in speciation (see Glossary) 12, 13 and 14 make a re-evaluation timely. Whereas botanists traditionally stressed the diversity-generating potential of hybridization 2, 3 and 14, zoologists traditionally saw it as a process that limits diversification [1] and refer to it mainly in the contexts of hybrid zones (Box 1) and reinforcement of reproductive isolation [15]. Judging by the wide distribution of allopolyploidy among plants, many plant species might be of direct hybrid origin or descended from a hybrid species in the recent past [16]. The ability to reproduce asexually might explain why allopolyploid hybrid species are more common in plants than in animals. Allopolyploidy arises when meiotic mismatch of parental chromosomes or karyotypes causes hybrid sterility. Mitotic error, duplicating the karyotype, can restore an asexually maintained hybrid line to fertility. Although bisexual allopolyploid hybrid species are not uncommon in fish [17] and frogs [18], the difficulty with which allopolyploid animals reproduce, typically requiring gynogenesis[19], makes establishment and survival of allopolyploid animal species difficult.

Genetic diversity, phenotypic variation and local adaptation in the alpine landscape: case studies with alpine plant species

Plant survival in alpine landscapes is constantly challenged by the harsh and often unpredictable environmental conditions. Steep environmental gradients and patchy distribution of habitats lead to small size and spatial isolation of populations and restrict gene flow. Agricultural land use has further increased the diversity of habitats below and above the treeline. We studied the consequences of the highly structured alpine landscape for evolutionary processes in four study plants: Epilobium fleischeri, Geum reptans, Campanula thyrsoides and Poa alpina. The main questions were: (1) How is genetic diversity distributed within and among populations and is it affected by altitude, population size or land use? (2) Do reproductive traits such as allocation to sexual or vegetative reproduction vary with altitude or land use? Furthermore, we studied if seed weight increases with altitude. Within-population genetic diversity of the four species was high and mostly not related to altitude and population size. Nevertheless, genetic differentiation among populations was pronounced and strongly increasing with distance. In Poa alpina genetic diversity was affected by land use. Results suggest considerable genetic drift among populations of alpine plants. Reproductive allocation was affected by altitude and land use in Poa alpina and by succession in Geum reptans. Seed weight was usually higher in alpine species than in related lowland species. We conclude that the evolutionary potential to respond to global change is mostly intact in alpine plants, even at high altitude. Phenotypic variability is shaped by adaptive as well as by random evolutionary processes; moreover plastic responses to growth conditions seem to be crucial for survival of plants in the alpine landscape.

Genomics of Divergence along a Continuum of Parapatric Population Differentiation

The patterns of genomic divergence during ecological speciation are shaped by a combination of evolutionary forces. Processes such as genetic drift, local reduction of gene flow around genes causing reproductive isolation, hitchhiking around selected variants, variation in recombination and mutation rates are all factors that can contribute to the heterogeneity of genomic divergence. On the basis of 60 fully sequenced three-spined stickleback genomes, we explore these different mechanisms explaining the heterogeneity of genomic divergence across five parapatric lake and river population pairs varying in their degree of genetic differentiation. We find that divergent regions of the genome are mostly specific for each population pair, while their size and abundance are not correlated with the extent of genome-wide population differentiation. In each pair-wise comparison, an analysis of allele frequency spectra reveals that 25–55% of the divergent regions are consistent with a local restriction of gene flow. Another large proportion of divergent regions (38–75%) appears to be mainly shaped by hitchhiking effects around positively selected variants. We provide empirical evidence that alternative mechanisms determining the evolution of genomic patterns of divergence are not mutually exclusive, but rather act in concert to shape the genome during population differentiation, a first necessary step towards ecological speciation.

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.

Structure and dynamics of hybrid zones at different stages of speciation in the common vole (Microtus arvalis)

The genetic structure and dynamics of hybrid zones provide crucial information for understanding the processes and mechanisms of evolutionary divergence and speciation. In general, higher levels of evolutionary divergence between taxa are more likely to be associated with reproductive isolation and may result in suppressed or strongly restricted hybridization. In this study, we examined two secondary contact zones between three deep evolutionary lineages in the common vole (Microtus arvalis). Differences in divergence times between the lineages can shed light on different stages of reproductive isolation and thus provide information on the ongoing speciation process in M. arvalis. We examined more than 800 individuals for mitochondrial (mtDNA), Y-chromosome and autosomal markers and used assignment and cline analysis methods to characterize the extent and direction of gene flow in the contact zones. Introgression of both autosomal and mtDNA markers in a relatively broad area of admixture indicates selectively neutral hybridization between the least-divergent lineages (Central and Eastern) without evidence for partial reproductive isolation. In contrast, a very narrow area of hybridization, shifts in marker clines and the quasi-absence of Y-chromosome introgression support a moving hybrid zone and unidirectional selection against male hybrids between the lineages with older divergence (Central and Western). Data from a replicate transect further support non-neutral processes in this hybrid zone and also suggest a role for landscape history in the movement and shaping of geneflow profiles.