Variation in mitochondrial DNA in a cline of allele frequencies and shell phenotype in the dog-whelk Nucella lapillus (L.)

Genetic constitution in the intertidal gastropod Nucella lapillus (L.) influences shell shape, growth rate and physiology. Clinal variation in these traits along a 5 km stretch of coastline in south Devon can be related to environmental variation in temperature and desiccation stress. We have examined mtDNA variation along this shore to investigate whether the cline represents primary or secondary contact. Two distinct mtDNA haplotypes were found which exhibit coincident step clines with karyotypic, allozymic and phenotypic variation and covary with the environmental pressures of temperature and desiccation. These results are interpreted in the context of the wider scale distribution of genetic and phenotypic variation in N. lapillus. It is suggested that the shore studied may represent one of a number of regions of secondary contact within a mosaic hybrid zone in N. lapillus , where coadapted phenotypic variation correlates with habitat and the position of the clines represents an environmental transition.     

Patterns of genetic,phenotypic, and acoustic variation across a chiffchaff (Phylloscopus collybita abietinus/tristis) hybrid zone

Characterizing patterns of evolution of genetic and phenotypic divergence between incipient species is essential to understand how evolution of reproductive isolation proceeds. Hybrid zones are excellent for studying such processes, as they provide opportunities to assess trait variation in individuals with mixed genetic background and to quantify gene flow across different genomic regions. Here, we combine plumage, song, mtDNA and whole‐genome sequence data and analyze variation across a sympatric zone between the European and the Siberian chiffchaff (Phylloscopus collybita abietinus/tristis) to study how gene exchange between the lineages affects trait variation. Our results show that chiffchaff within the sympatric region show more extensive trait variation than allopatric birds, with a large proportion of individuals exhibiting intermediate phenotypic characters. The genomic differentiation between the subspecies is lower in sympatry than in allopatry and sympatric birds have a mix of genetic ancestry indicating extensive ongoing and past gene flow. Patterns of phenotypic and genetic variation also vary between regions within the hybrid zone, potentially reflecting differences in population densities, age of secondary contact, or differences in mate recognition or mate preference. The genomic data support the presence of two distinct genetic clades corresponding to allopatric abietinus and tristis and that genetic admixture is the force underlying trait variation in the sympatric region—the previously described subspecies (“fulvescens”) from the region is therefore not likely a distinct taxon. In addition, we conclude that subspecies identification based on appearance is uncertain as an individual with an apparently distinct phenotype can have a considerable proportion of the genome composed of mixed alleles, or even a major part of the genome introgressed from the other subspecies. Our results provide insights into the dynamics of admixture across subspecies boundaries and have implications for understanding speciation processes and for the identification of specific chiffchaff individuals based on phenotypic characters.     

Geographic structure of genetic and phenotypic variation in the hybrid zone between Quercus affinis and Q. laurina in Mexico

Analyzing the structure of hybrid zones is important for inferring their origin, dynamics and evolutionary significance. We examined the geographic structure of phenotypic and genetic variation in the contact zone between two Mexican red oaks, Quercus affinis and Q. laurina. A total of 105 individuals from seven populations were sampled along a 600‐km latitudinal gradient representing the distribution area of the two species and their contact zone. Individuals were genotyped for nine nuclear and four chloroplast DNA microsatellite loci (ncSSR and cpSSR, respectively), and characterized for several leaf and acorn traits. The cpSSR data revealed extensive haplotype sharing among populations of the two species, while a Bayesian assignment analysis based on ncSSRs identified two main genetic groups, each corresponding to one of the species, and two populations in the contact zone showing evidence of admixture. The proportion of genetic ancestry in the populations was strongly associated with latitude and showed a pattern of variation with the shape of a narrow sigmoidal cline. The variation in three of the seven phenotypic traits was partially congruent with molecular variation, while the other traits did not conform to a geographic cline but instead were correlated with environmental variables. In conclusion, the hybrid zone between the two oak species has some of the characteristics of a tension zone, but heterogeneous variation across traits suggests differential introgression and the action of extrinsic selection.     

Post‐glacial colonization routes coincide with a life‐history breakpoint along a latitudinal gradient

Although adaptive divergence along environmental gradients has repeatedly been demonstrated, the role of post‐glacial colonization routes in determining phenotypic variation along gradients has received little attention. Here, we used a hierarchical Q_ST–F_ST approach to separate the roles of adaptive and neutral processes in shaping phenotypic variation in moor frog (Rana arvalis) larval life histories along a 1,700 km latitudinal gradient across northern Europe. This species has colonized Scandinavia via two routes with a contact zone in northern Sweden. By using neutral SNP and common garden phenotypic data from 13 populations at two temperatures, we showed that most of the variation along the gradient occurred between the two colonizing lineages. We found little phenotypic divergence within the lineages; however, all phenotypic traits were strongly diverged between the southern and northern colonization routes, with higher growth and development rates and larger body size in the north. The Q_ST estimates between the colonization routes were four times higher than F_ST, indicating a prominent role for natural selection. Q_ST within the colonization routes did not generally differ from F_ST, but we found temperature‐dependent adaptive divergence close to the contact zone. These results indicate that lineage‐specific variation can account for much of the adaptive divergence along a latitudinal gradient.     

Adaptation in a variable environment: Phenotypic plasticity and bet‐hedging during egg diapause and hatching in an annual killifish

Two ways in which organisms adapt to variable environments are phenotypic plasticity and bet‐hedging. Theory suggests that bet‐hedging is expected to evolve in unpredictable environments for which reliable cues indicative of future conditions (or season length) are lacking. Alternatively, if reliable cues exist indicating future conditions, organisms will be under selection to produce the most appropriate phenotype —that is, adaptive phenotypic plasticity. Here, we experimentally test which of these modes of adaptation are at play in killifish that have evolved an annual life cycle. These fish persist in ephemeral pools that completely dry each season through the production of eggs that can remain in developmental arrest, or diapause, buried in the soil, until the following rainy season. Consistent with diversified bet‐hedging (a risk spreading strategy), we demonstrate that the eggs of the annual killifish Nothobranchius furzeri exhibit variation at multiple levels—whether or not different stages of diapause are entered, for how long diapause is entered, and the timing of hatching—and this variation persists after controlling for both genetic and environmental sources of variation. However, we show that phenotypic plasticity is also present in that the proportion of eggs that enter diapause is influenced by environmental factors (temperature and light level) that vary seasonally. In nature there is typically a large parameter zone where environmental cues are somewhat correlated with seasonality, but not perfectly so, such that it may be advantageous to have a combination of both bet‐hedging and plasticity.     

Phenotypic effects of the Y chromosome are variable and structured in hybrids among house mouse recombinant lines

Hybrid zones between divergent populations sieve genomes into blocks that introgress across the zone, and blocks that do not, depending on selection between interacting genes. Consistent with Haldane's rule, the Y chromosome has been considered counterselected and hence not to introgress across the European house mouse hybrid zone. However, recent studies detected massive invasion of M. m. musculus Y chromosomes into M. m. domesticus territory. To understand mechanisms facilitating Y spread, we created 31 recombinant lines from eight wild‐derived strains representing four localities within the two mouse subspecies. These lines were reciprocally crossed and resulting F1 hybrid males scored for five phenotypic traits associated with male fitness. Molecular analyses of 51 Y‐linked SNPs attributed ~50% of genetic variation to differences between the subspecies and 8% to differentiation within both taxa. A striking proportion, 21% (frequencies of sperm head abnormalities) and 42% (frequencies of sperm tail dissociations), of phenotypic variation was explained by geographic Y chromosome variants. Our crossing design allowed this explanatory power to be examined across a hierarchical scale from subspecific to local intrastrain effects. We found that divergence and variation were expressed diversely in different phenotypic traits and varied across the whole hierarchical scale. This finding adds another dimension of complexity to studies of Y introgression not only across the house mouse hybrid zone but potentially also in other contact zones.     

Genetic and phenotypic variation across a hybrid zone between ecologically divergent tree squirrels (Tamiasciurus)

A hybrid zone along an environmental gradient should contain a clinal pattern of genetic and phenotypic variation. This occurs because divergent selection in the two parental habitats is typically strong enough to overcome the homogenizing effects of gene flow across the environmental transition. We studied hybridization between two parapatric tree squirrels (Tamiasciurus spp.) across a forest gradient over which the two species vary in coloration, cranial morphology and body size. We sampled 397 individuals at 29 locations across a 600‐km transect to seek genetic evidence for hybridization; upon confirming hybridization, we examined levels of genetic admixture in relation to maintenance of phenotypic divergence despite potentially homogenizing gene flow. Applying population assignment analyses to microsatellite data, we found that Tamiasciurus douglasii and T. hudsonicus form two distinct genetic clusters but also hybridize, mostly within transitional forest habitat. Overall, based on this nuclear analysis, 48% of the specimens were characterized as T. douglasii, 9% as hybrids and 43% as T. hudsonicus. Hybrids appeared to be reproductively viable, as evidenced by the presence of later‐generation hybrid genotypes. Observed clines in ecologically important phenotypic traits—fur coloration and cranial morphology—were sharper than the cline of putatively neutral mtDNA, which suggests that divergent selection may maintain phenotypic distinctiveness. The relatively recent divergence of these two species (probably late Pleistocene), apparent lack of prezygotic isolating mechanisms and geographic coincidence of cline centres for both genetic and phenotypic variation suggest that environmental factors play a large role in maintaining the distinctiveness of these two species across the hybrid zone.     

The effect of phenotypic variation on metapopulation persistence

Demographic stochasticity (due to the probabilistic nature of the birth–death process) and demographic heterogeneity (between-individual differences in demographic parameters) have long been seen as factors affecting extinction risk. While demographic stochasticity can be independent of underlying species traits, demographic heterogeneity may strongly depend on phenotypic variation. However, how phenotypic variation can affect extinction risk is largely unknown. Here, I develop a stochastic metapopulation model that takes into account the effects of demographic stochasticity and phenotypic variation in the traits controlling colonization rates to assess what the effect of phenotypic variation may be on the persistence of the metapopulation. Although phenotypic variation can lead to a decrease in metapopulation persistence under some conditions, it also may lead to an increase in persistence whenever phenotypic mismatch—or the distance between the optimal trait value and the population mean—is large. This mismatch can in turn arise from a variety of ecological and evolutionary reasons, including weak selection or a recent history of invasion. Last, the effect of phenotypic variation has a deterministic component on colonization rates, and a stochastic component on persistence through colonization rates, but both are important to understand the overall effect. These results have important implications for the conservation of threatened species and management practices that may historically have overlooked phenotypic variation as unimportant noise around mean values of interest.     

Fine-scale phenotypic change across a species transition zone in the genus neotoma: disentangling independent evolution from phylogenetic history

Zones of secondary contact between closely related species provide a rare opportunity to examine evidence of evolutionary processes that reinforce species boundaries and/or promote diversification. Here, we report on genetic and morphological variation in two sister species of woodrats, Neotoma fuscipes and N. macrotis, across a 30-km transition zone in the Sierra Nevada of California. We assessed whether these lineages readily hybridize, and whether their morphology suggests ecological interactions favoring phenotypic diversification. We combined measurements of body size and 11 craniodental traits from nine populations with genetic data to examine patterns of variation within and between species. We used phylogenetic autocorrelation methods to estimate the degree to which phenotypic variation in our dataset arose from independent evolution within populations versus phylogenetic history. Although no current sympatry or hybridization was evident, craniodental morphology diverged in both lineages near their distributional limits, whereas body size converged. The shift in craniodental morphology arose independently within populations whereas body size retained a strong phylogenetic signal, yet both patterns are consistent with expectations of phenotypic change based on different models of resource competition. Our findings demonstrate the importance of examining a suite of morphological traits across contact zones to provide a more complete picture of potential ecological interactions: competition may drive both diversification and convergence in different phenotypic traits.     

The genetic basis of color‐related local adaptation in a ring‐like colonization around the Mediterranean

Uncovering the genetic basis of phenotypic variation and the population history under which it established is key to understand the trajectories along which local adaptation evolves. Here, we investigated the genetic basis and evolutionary history of a clinal plumage color polymorphism in European barn owls (Tyto alba). Our results suggest that barn owls colonized the Western Palearctic in a ring‐like manner around the Mediterranean and meet in secondary contact in Greece. Rufous coloration appears to be linked to a recently evolved nonsynonymous‐derived variant of the melanocortin 1 receptor (MC1R) gene, which according to quantitative genetic analyses evolved under local adaptation during or following the colonization of Central Europe. Admixture patterns and linkage disequilibrium between the neutral genetic background and color found exclusively within the secondary contact zone suggest limited introgression at secondary contact. These results from a system reminiscent of ring species provide a striking example of how local adaptation can evolve from derived genetic variation.     

Morphological and genetic divergence of intralacustrine stickleback morphs in Iceland: a case for selective differentiation?

The evolutionary processes involved in population divergence and local adaptation are poorly understood. Theory predicts that divergence of adjacent populations is possible but depends on several factors including gene flow, divergent selection, population size and the number of genes involved in divergence and their distribution on the genome. We analyse variation in neutral markers, markers linked to putative quantitative trait loci and morphological traits in a recent (<10000 years) zone of primary divergence between stickleback morphs in Lake Thingvallavatn, Iceland. Environmental factors, especially predation, are clearly implicated in reducing gene flow between morphs. There is continuous morphological and genetic variation between habitats with a zone centre similar to secondary contact zones. Individual microsatellite loci are implicated as being linked to adaptive variation by direct tests as well as by differences in cline shape. Patterns of linkage disequilibria indicate that the morphs have diverged at several loci. This divergence shows parallels and differences with the well-studied limnetic-benthic stickleback morphs, both in phenotypic divergence and at the genomic level.     

Parallel plumage colour evolution and introgressive hybridization in wheatears

Genetic and phenotypic mosaics, in which various phenotypes and different genomic regions show discordant patterns of species or population divergence, offer unique opportunities to study the role of ancestral and introgressed genetic variation in phenotypic evolution. Here, we investigated the evolution of discordant phenotypic and genetic divergence in a monophyletic clade of four songbird taxa—pied wheatear (O. pleschanka), Cyprus wheatear (Oenanthe cypriaca), and western and eastern subspecies of black‐eared wheatear (O. h. hispanica and O. h. melanoleuca). Phenotypically, black back and neck sides distinguish pied and Cyprus wheatears from the white‐backed/necked black‐eared wheatears. Meanwhile, mitochondrial variation only distinguishes western black‐eared wheatear. In the absence of nuclear genetic data, and given frequent hybridization among eastern black‐eared and pied wheatear, it remains unclear whether introgression is responsible for discordance between mitochondrial divergence patterns and phenotypic similarities, or whether plumage coloration evolved in parallel. Multispecies coalescent analyses of about 20,000 SNPs obtained from RAD data mapped to a draft genome assembly resolve the species tree, provide evidence for the parallel evolution of colour phenotypes and establish western and eastern black‐eared wheatears as independent taxa that should be recognized as full species. The presence of the entire admixture spectrum in the Iranian hybrid zone and the detection of footprints of introgression from pied into eastern black‐eared wheatear beyond the hybrid zone despite strong geographic structure of ancestry proportions furthermore suggest a potential role for introgression in parallel plumage colour evolution. Our results support the importance of standing heterospecific and/or ancestral variation in phenotypic evolution.     

Phenotypic and genetic diversity in aposematic Malagasy poison frogs (genus Mantella)

Intraspecific color variation has long fascinated evolutionary biologists. In species with bright warning coloration, phenotypic diversity is particularly compelling because many factors, including natural and sexual selection, contribute to intraspecific variation. To better understand the causes of dramatic phenotypic variation in Malagasy poison frogs, we quantified genetic structure and color and pattern variation across three closely related species, Mantella aurantiaca, Mantella crocea, and Mantella milotympanum. Although our restriction site‐associated DNA (RAD) sequencing approach identified clear genetic clusters, they do not align with current species designations, which has important conservation implications for these imperiled frogs. Moreover, our results suggest that levels of intraspecific color variation within this group have been overestimated, while species diversity has been underestimated. Within major genetic clusters, we observed distinct patterns of variation including: populations that are phenotypically similar yet genetically distinct, populations where phenotypic and genetic breaks coincide, and populations that are genetically similar but have high levels of within‐population phenotypic variation. We also detected admixture between two of the major genetic clusters. Our study suggests that several mechanisms—including hybridization, selection, and drift—are contributing to phenotypic diversity. Ultimately, our work underscores the need for a reevaluation of how polymorphic and polytypic populations and species are classified, especially in aposematic organisms.     

Integrating genetic and environmental forces that shape the evolution of geographic variation in a marine snail

Temporal and spatial patterns of phenotypic variation have traditionally been thought to reflect genetic differentiation produced by natural selection. Recently, however, there has been growing interest in how natural selection may shape the genetics of phenotypic plasticity to produce patterns of geographic variation and phenotypic evolution. Because the covariance between genetic and environmental influences can modulate the expression of phenotypic variation, a complete understanding of geographic variation requires determining whether these influences covary in the same (cogradient variation) or in opposing (countergradient variation) directions. We focus on marine snails from rocky intertidal shores as an ideal system to explore how genetic and plastic influences contribute to geographic and historical patterns of phenotypic variation. Phenotypic plasticity in response to predator cues, wave action, and water temperature appear to exert a strong influence on small and large-scale morphological variation in marine snails. In particular, plasticity in snail shell thickness: (i) may contribute to phenotypic evolution, (ii) appears to have evolved across small and large spatial scales, and (iii) may be driven by life history trade-offs tied to architectural constraints imposed by the shell. The plasticity exhibited by these snails represents an important adaptive strategy to the pronounced heterogeneity of the intertidal zone and undoubtedly has played a key role in their evolution.     

Interspecific territoriality in gibbons (Hylobates lar and H. pileatus) and its effects on the dynamics of interspecies contact zones

We investigated the ecology and interspecific interactions of the two gibbon species (Hylobates lar and H. pileatus) that overlap in distribution within a narrow zone of contact in the headwaters of the Takhong River at Khao Yai National Park in central Thailand. The zone is about 10-km wide, with phenotypic hybrids comprising 6.5% of the adult population. We compared species with respect to diet, territory size, intra- and interspecific group encounters, and territory quality. The two gibbon species exploited the same types of resources within their territories despite variation in the relative abundance of food-plant species between territories. The gibbons were interspecifically territorial, and males of both species displayed aggressive behaviors at both intraspecific and interspecific territorial boundaries. There was no marked difference in the amount of overlap between territories of conspecific and heterospecific pairs of groups. Although the habitat was not homogeneous, territory quality did not vary significantly between species. The species have not diverged in habitat preference or in diet. Neither species dominated in interspecific encounters, and both were reproducing well in the contact zone. We analyzed the potential consequences of several types of interspecific interactions on individual dispersal options and on the structure of the contact zone. Interference competition through interspecific territoriality affects the dispersal of individuals into the range of the other species. In general, territorial competition coupled with limited hybridization leads to predictions of a narrow contact zone or parapatry between species; thus, behavioral and ecological interactions between species need to be considered as potential factors in explaining range borders of primate species.     

Selection on outlier loci and their association with adaptive phenotypes in Littorina saxatilis contact zones

A fundamental issue in speciation research is to evaluate phenotypic variation and the genomics driving the evolution of reproductive isolation between sister taxa. Above all, hybrid zones are excellent study systems for researchers to examine the association of genetic differentiation, phenotypic variation and the strength of selection. We investigated two contact zones in the marine gastropod Littorina saxatilis and utilized landmark‐based geometric morphometric analysis together with amplified fragment length polymorphism (AFLP) markers to assess phenotypic and genomic divergence between ecotypes under divergent selection. From genetic markers, we calculated the cline width, linkage disequilibrium and the average effective selection on a locus. Additionally, we conducted an association analysis linking the outlier loci and phenotypic variation between ecotypes and show that a proportion of outlier loci are associated with key adaptive phenotypic traits.     

Sexual signal evolution and patterns of assortative mating across an intraspecific contact zone

Contact zones provide important insights into the evolutionary processes that underlie lineage divergence and speciation. Here, we use a contact zone to ascertain speciation potential in the red-eyed treefrog (Agalychnis callidryas), a brightly coloured and polymorphic frog that exhibits unusually high levels of intraspecific variation. Populations of A. callidryas differ in a number of traits, several of which are known sexual signals that mediate premating reproductive isolation in allopatric populations. Along the Caribbean coast of Costa Rica, a ~100 km contact zone, situated between two phenotypically and genetically divergent parent populations, contains multiple colour pattern phenotypes and late-generation hybrids. This contact zone provides the opportunity to examine processes that are important in the earliest stages of lineage divergence. We performed analyses of colour pattern variation in five contact zone sites and six parental sites and found complex, continuous colour variation along the contact zone. We found discordance between the geographic distribution of colour pattern and previously described genomic population structure. We then used a parental site and contact zone site to measure assortative mating and directional selection from naturally-occurring amplectant mating pairs. We found assortative mating in a parental population, but no assortative mating in the contact zone. Furthermore, we uncovered evidence of directional preference towards the adjacent parental phenotype in the contact zone population, but no directional preference in the parent population. Combined, these data provide insights into potential dynamics at the contact zone borders and indicate that incipient speciation between parent populations will be slowed.     

Where the Lake Meets the Sea: Strong Reproductive Isolation Is Associated with Adaptive Divergence between Lake Resident and Anadromous Three-Spined Sticklebacks

Contact zones between divergent forms of the same species are often characterised by high levels of phenotypic diversity over small geographic distances. What processes are involved in generating such high phenotypic diversity? One possibility is that introgression and recombination between divergent forms in contact zones results in greater phenotypic and genetic polymorphism. Alternatively, strong reproductive isolation between forms may maintain distinct phenotypes, preventing homogenisation by gene flow. Contact zones between divergent freshwater-resident and anadromous stickleback (Gasterosteus aculeatus L.) forms are numerous and common throughout the species distribution, offering an opportunity to examine these contrasting hypotheses in greater detail. This study reports on an interesting new contact zone located in a tidally influenced lake catchment in western Ireland, characterised by high polymorphism for lateral plate phenotypes. Using neutral and QTL-linked microsatellite markers, we tested whether the high diversity observed in this contact zone arose as a result of introgression or reproductive isolation between divergent forms: we found strong support for the latter hypothesis. Three phenotypic and genetic clusters were identified, consistent with two divergent resident forms and a distinct anadromous completely plated population that migrates in and out of the system. Given the strong neutral differentiation detected between all three morphotypes (mean F_ST = 0.12), we hypothesised that divergent selection between forms maintains reproductive isolation. We found a correlation between neutral genetic and adaptive genetic differentiation that support this. While strong associations between QTL linked markers and phenotypes were also observed in this wild population, our results support the suggestion that such associations may be more complex in some Atlantic populations compared to those in the Pacific. These findings provide an important foundation for future work investigating the dynamics of gene flow and adaptive divergence in this newly discovered stickleback contact zone.     

Introgressive hybridization and morphological transgression in the contact zone between two Mediterranean Solea species

Hybrid zones provide natural experiments where new combinations of genotypes and phenotypes are produced. Studying the reshuffling of genotypes and remodeling of phenotypes in these zones is of particular interest to document the building of reproductive isolation and the possible emergence of transgressive phenotypes that can be a source of evolutionary novelties. Here, we specifically investigate the morphological variation patterns associated with introgressive hybridization between two species of sole, Solea senegalensis and Solea aegyptiaca. The relationship between genetic composition at nuclear loci and individual body shape variation was studied in four populations sampled across the hybrid zone located in northern Tunisia. A strong correlation between genetic and phenotypic variation was observed among all individuals but not within populations, including the two most admixed ones. Morphological convergence between parental species was observed close to the contact zone. Nevertheless, the samples taken closest to the hybrid zone also displayed deviant segregation of genotypes and phenotypes, as well as transgressive phenotypes. In these samples, deviant body shape variation could be partly attributed to a reduced condition index, and the distorted genetic composition was most likely due to missing allelic combinations. These results were interpreted as an indication of hybrid breakdown, which likely contributes to postmating reproductive isolation between the two species.