Genetic and epigenetic variations associated with adaptation to heterogeneous habitat conditions in a deciduous shrub

Environmentally induced phenotypic plasticity is thought to play an important role in the adaption of plant populations to heterogeneous habitat conditions, and yet the importance of epigenetic variation as a mechanism of adaptive plasticity in natural plant populations still merits further research. In this study, we investigated populations of Vitex negundo var. heterophylla (Chinese chastetree) from adjacent habitat types at seven sampling sites. Using several functional traits, we detected a significant differentiation between habitat types. With amplified fragment length polymorphisms (AFLP) and methylation‐sensitive AFLP (MSAP), we found relatively high levels of genetic and epigenetic diversity but very low genetic and epigenetic differences between habitats within sites. Bayesian clustering showed a remarkable habitat‐related differentiation and more genetic loci associated with the habitat type than epigenetic, suggesting that the adaptation to the habitat is genetically based. However, we did not find any significant correlation between genetic or epigenetic variation and habitat using simple and partial Mantel tests. Moreover, we found no correlation between genetic and ecologically relevant phenotypic variation and a significant correlation between epigenetic and phenotypic variation. Although we did not find any direct relationship between epigenetic variation and habitat environment, our findings suggest that epigenetic variation may complement genetic variation as a source of functional phenotypic diversity associated with adaptation to the heterogeneous habitat in natural plant populations.     

Genetic and morphological divergence between Littorina fabalis ecotypes in Northern Europe

Low dispersal marine intertidal species facing strong divergent selective pressures associated with steep environmental gradients have a great potential to inform us about local adaptation and reproductive isolation. Among these, gastropods of the genus Littorina offer a unique system to study parallel phenotypic divergence resulting from adaptation to different habitats related with wave exposure. In this study, we focused on two Littorina fabalis ecotypes from Northern European shores and compared patterns of habitat‐related phenotypic and genetic divergence across three different geographic levels (local, regional and global). Geometric morphometric analyses revealed that individuals from habitats moderately exposed to waves usually present a larger shell size with a wider aperture than those from sheltered habitats. The phenotypic clustering of L. fabalis by habitat across most locations (mainly in terms of shell size) support an important role of ecology in morphological divergence. A genome scan based on amplified fragment length polymorphisms (AFLPs) revealed a heterogeneous pattern of differentiation across the genome between populations from the two different habitats, suggesting ecotype divergence in the presence of gene flow. The contrasting patterns of genetic structure between nonoutlier and outlier loci, and the decreased sharing of outlier loci with geographic distance among locations are compatible with parallel evolution of phenotypic divergence, with an important contribution of gene flow and/or ancestral variation. In the future, model‐based inference studies based on sequence data across the entire genome will help unravelling these evolutionary hypotheses, improving our knowledge about adaptation and its influence on diversification within the marine realm.     

Genomic differentiation and patterns of gene flow between two long‐tailed tit species (Aegithalos)

Patterns of heterogeneous genomic differentiation have been well documented between closely related species, with some highly differentiated genomic regions (“genomic differentiation islands”) spread throughout the genome. Differential levels of gene flow are proposed to account for this pattern, as genomic differentiation islands are suggested to be resistant to gene flow. Recent studies have also suggested that genomic differentiation islands could be explained by linked selection acting on genomic regions with low recombination rates. Here, we investigate genomic differentiation and gene‐flow patterns for autosomes using RAD‐seq data between two closely related species of long‐tailed tits (Aegithalos bonvaloti and A. fuliginosus) in both allopatric and contact zone populations. The results confirm recent or ongoing gene flow between these two species. However, there is little evidence that the genomic regions that were found to be highly differentiated between the contact zone populations are resistant to gene flow, suggesting that differential levels of gene flow is not the cause of the heterogeneous genomic differentiation. Linked selection may be the cause of genomic differentiation islands between the allopatric populations with no or very limited gene flow, but this could not account for the heterogeneous genomic differentiation between the contact zone populations, which show evidence of recent or ongoing gene flow.     

Habitat quality as a predictor of spatial variation in blue tit reproductive performance: a multi-plot analysis in a heterogeneous landscape

Vertebrate studies have rarely investigated the influence of spatial variation in habitat richness on both short-term (breeding) and long-term (offspring recruitment) reproductive performance using simultaneously multi-patch, multi-habitat type and multi-year approaches at landscape level. Here we present results of such an approach using the influence of two oak tree (Quercus ilex, Q. humilis) species on reproductive performance in Corsican blue tits (Parus caeruleus ogliastrae) as a model system. We found that blue tits breeding in rich broad-leaved deciduous patches consistently laid eggs earlier in the season, and produced larger clutches and more fledglings of higher quality, than those breeding in poor evergreen patches. Also, parents, especially males, were in better physical condition in the broad-leaved deciduous than in the evergreen patches. Surprisingly, estimates of long-term effects of reproduction, such as recruitment rates of locally born offspring, did not differ between the two habitat types. Our results suggest that short-term breeding performance and phenotypic quality of both chicks and parents do not necessarily provide reliable information about contributions to following generations at a scale larger than that of the local study plot. Differences in reproductive performance between the two oak habitat types could not be attributed to density-dependent effects, differences in levels of nest predation, or differences in age structure of the birds. We suggest that habitats that are optimal for breeding are not necessarily optimal for survival after the breeding season.     

Personality‐related differences in response to habitat in Mediterranean blue tits

Habitat‐specific selection pressures have been widely recognized, but whether selection favours different personality types in different habitats has rarely been evaluated. This study aimed to test whether personality‐related differences in annual reproductive success differed between two populations of blue tits (Cyanistes caeruleus) occupying different Mediterranean habitats (oakwood and pinewood). We measured exploration and parental provisioning behaviours and used a path analysis to ask how the interplay between these two behavioural traits affected reproductive success in each of these two habitats. We found that blue tits breeding in the pinewood were slow‐exploring compared to blue tits breeding in the oakwood, suggesting the occurrence of personality‐related differences in settlement, or behavioural plasticity in response to habitat. Exploration behaviour and feeding rates were positively associated, suggesting that they affect each other or that there is an environmental factor affecting both traits simultaneously. Finally, fast explorers were favoured in the pinewood, while there was no selection acting on exploration behaviour in the oak habitat. These findings emphasize the importance of integrating habitat selection, plasticity and personality in the study of behavioural evolution.     

Evolution along the Great Rift Valley: phenotypic and genetic differentiation of East African white‐eyes (Aves,Zosteropidae)

The moist and cool cloud forests of East Africa represent a network of isolated habitats that are separated by dry and warm lowland savannah, offering an opportunity to investigate how strikingly different selective regimes affect species diversification. Here, we used the passerine genus Zosterops (white‐eyes) from this region as our model system. Species of the genus occur in contrasting distribution settings, with geographical mountain isolation driving diversification, and savannah interconnectivity preventing differentiation. We analyze (1) patterns of phenotypic and genetic differentiation in high‐ and lowland species (different distribution settings), (2) investigate the potential effects of natural selection and temporal and spatial isolation (evolutionary drivers), and (3) critically review the taxonomy of this species complex. We found strong phenotypic and genetic differentiation among and within the three focal species, both in the highland species complex and in the lowland taxa. Altitude was a stronger predictor of phenotypic patterns than the current taxonomic classification. We found longitudinal and latitudinal phenotypic gradients for all three species. Furthermore, wing length and body weight were significantly correlated with altitude and habitat type in the highland species Z. poliogaster. Genetic and phenotypic divergence showed contrasting inter‐ and intraspecific structures. We suggest that the evolution of phenotypic characters is mainly driven by natural selection due to differences in the two macro‐habitats, cloud forest and savannah. In contrast, patterns of neutral genetic variation appear to be rather driven by geographical isolation of the respective mountain massifs. Populations of the Z. poliogaster complex, as well as Z. senegalensis and Z. abyssinicus, are not monophyletic based on microsatellite data and have higher levels of intraspecific differentiation compared to the currently accepted species.     

Desert‐like badlands and surrounding (semi‐)dry grasslands of Central Germany promote small‐scale phenotypic and genetic differentiation in Thymus praecox

Environmental heterogeneity among sites can generate phenotypic and genetic variation facilitating differentiation and microevolution of plant populations. Badlands are desert‐like, predominantly vegetation‐poor habitats often embedded in (semi‐)dry grasslands. The desert‐like conditions of badlands demand extreme adaptation of plants, that is, phenotypic modifications in short‐term and/or natural adaptation in long‐term. However, detailed knowledge is missing about both plant phenotypic and genetic differentiation in this unique and widely occurring habitat type. The present study focused on the largest known badlands systems in Central Europe located in the “Drei Gleichen” region, a designated nature conservation area in Central Germany. Locations were suitable for this study in terms of having co‐occurring badlands and (semi‐)dry grassland habitats (sites) occupied by the pioneer plant Thymus praecox. Here, we studied the environmental preferences, morphological and functional trait variation, and genetic variation using microsatellite markers of T. praecox. Results revealed significant, mainly site‐dependent environmental, phenotypic, and genetic differentiation. In general, individuals in badlands are shorter in height and have lower patch sizes (length × width), relative growth rates, and smaller stomata. The PCA additionally unveiled slightly increased leaf robustness, trichome density, decreased stomatal conductance, fewer females, and earlier phenology in badlands. We interpret differentiation patterns as adaptive responses to light, temperature, drought, and nutrient stress conditions supported by reviewed literature. Genetic differentiation was strongest between local badlands and grassland sites, and clearly weaker among locations and between sites (in total) as indicated by G_ST, AMOVA, PCoA, and population structure. Our study supports the importance of small‐scale microhabitat conditions as a driver of microevolutionary processes, and the population's need for sufficient phenotypic variation and genetic resources to deal with environmental changes. We demonstrated that badlands are an appropriate model system for testing plant response to extreme habitats and that more research is needed on these fascinating landscapes.     

Parallel pattern of differentiation at a genomic island shared between clinal and mosaic hybrid zones in a complex of cryptic seahorse lineages

Diverging semi‐isolated lineages either meet in narrow clinal hybrid zones, or have a mosaic distribution associated with environmental variation. Intrinsic reproductive isolation is often emphasized in the former and local adaptation in the latter, although both reduce gene flow between groups. Rarely are these two patterns of spatial distribution reported in the same study system. Here, we report that the long‐snouted seahorse Hippocampus guttulatus is subdivided into discrete panmictic entities by both types of hybrid zones. Along the European Atlantic coasts, a northern and a southern lineage meet in the southwest of France where they coexist in sympatry—i.e., in the same geographical zone—with little hybridization. In the Mediterranean Sea, two lineages have a mosaic distribution, associated with lagoon‐like and marine habitats. A fifth lineage was identified in the Black Sea. Genetic homogeneity over large spatial scales contrasts with isolation maintained in sympatry or close parapatry at a fine scale. A high variation in locus‐specific introgression rates provides additional evidence that partial reproductive isolation must be maintaining the divergence. We find that fixed differences between lagoon and marine populations in the Mediterranean Sea belong to the most differentiated SNPs between the two Atlantic lineages, against the genome‐wide pattern of structure that mostly follow geography. These parallel outlier SNPs cluster on a single chromosome‐wide island of differentiation. Since Atlantic lineages do not map to lagoon‐sea habitat variation, genetic parallelism at the genomic island suggests a shared genetic barrier contributes to reproductive isolation in contrasting contexts–i.e., spatial versus ecological. We discuss how a genomic hotspot of parallel differentiation could have evolved and become associated both with space and with a patchy environment in a single study system.     

Parallel genetic divergence among coastal–marine ecotype pairs of European anchovy explained by differential introgression after secondary contact

Ecophenotypic differentiation among replicate ecotype pairs within a species complex is often attributed to independent outcomes of parallel divergence driven by adaptation to similar environmental contrasts. However, the extent to which parallel phenotypic and genetic divergence patterns have emerged independently is increasingly questioned by population genomic studies. Here, we document the extent of genetic differentiation within and among two geographic replicates of the coastal and marine ecotypes of the European anchovy (Engraulis encrasicolus) gathered from Atlantic and Mediterranean locations. Using a genome‐wide data set of RAD‐derived SNPs, we show that habitat type (marine vs. coastal) is the most important component of genetic differentiation among populations of anchovy. By analysing the joint allele frequency spectrum of each coastal–marine ecotype pair, we show that genomic divergence patterns between ecotypes can be explained by a postglacial secondary contact following a long period of allopatric isolation (c. 300 kyrs). We found strong support for a model including heterogeneous migration among loci, suggesting that secondary gene flow has eroded past differentiation at different rates across the genome. Markers experiencing reduced introgression exhibited strongly correlated differentiation levels among Atlantic and Mediterranean regions. These results support that partial reproductive isolation and parallel genetic differentiation among replicate pairs of anchovy ecotypes are largely due to a common divergence history prior to secondary contact. They moreover provide comprehensive insights into the origin of a surprisingly strong fine‐scale genetic structuring in a high gene flow marine fish, which should improve stock management and conservation actions.     

Recent local adaptation of sockeye salmon to glacial spawning habitats

Salmonids spawn in highly diverse habitats, exhibit strong genetic population structuring, and can quickly colonize newly created habitats with few founders. Spawning traits often differ among populations, but it is largely unknown if these differences are adaptive or due to genetic drift. To test if sockeye salmon (Oncorhynchus nerka) populations are adapted to glacial, beach, and tributary spawning habitats, we examined variation in heritable phenotypic traits associated with spawning in 13 populations of wild sockeye salmon in Lake Clark, Alaska. These populations were commonly founded between 100 and 400 hundred sockeye salmon generations ago and exhibit low genetic divergence at 11 microsatellite loci (F _ST < 0.024) that is uncorrelated with spawning habitat type. We found that mean P _ST (phenotypic divergence among populations) exceeded neutral F _ST for most phenotypic traits measured, indicating that phenotypic differences among populations could not be explained by genetic drift alone. Phenotypic divergence among populations was associated with spawning habitat differences, but not with neutral genetic divergence. For example, female body color was lighter and egg color was darker in glacial than non-glacial habitats. This may be due to reduced sexual selection for red spawning color in glacial habitats and an apparent trade-off in carotenoid allocation to body and egg color in females. Phenotypic plasticity is an unlikely source of phenotypic differences because Lake Clark sockeye salmon spend nearly all their lives in a common environment. Our data suggest that Lake Clark sockeye salmon populations are adapted to spawning in glacial, beach and tributary habitats and provide the first evidence of a glacial spawning ecotype in salmonids. Glacial spawning habitats are often young (i.e., <200 years old) and ephemeral. Thus, local adaptation of sockeye salmon to glacial habitats appears to have occurred recently.     

Contemporary ecotypic divergence during a recent range expansion was facilitated by adaptive introgression

Although rapid phenotypic evolution during range expansion associated with colonization of contrasting habitats has been documented in several taxa, the evolutionary mechanisms that underlie such phenotypic divergence have less often been investigated. A strong candidate for rapid ecotype formation within an invaded range is the three‐spine stickleback in the Lake Geneva region of central Europe. Since its introduction only about 140 years ago, it has undergone a significant expansion of its range and its niche, now forming phenotypically differentiated parapatric ecotypes that occupy either the pelagic zone of the large lake or small inlet streams, respectively. By comparing museum collections from different times with contemporary population samples, we here reconstruct the evolution of parapatric phenotypic divergence through time. Using genetic data from modern samples, we infer the underlying invasion history. We find that parapatric habitat‐dependent phenotypic divergence between the lake and stream was already present in the first half of the twentieth century, but the magnitude of differentiation increased through time, particularly in antipredator defence traits. This suggests that divergent selection between the habitats occurred and was stable through much of the time since colonization. Recently, increased phenotypic differentiation in antipredator defence traits likely results from habitat‐dependent selection on alleles that arrived through introgression from a distantly related lineage from outside the Lake Geneva region. This illustrates how hybridization can quickly promote phenotypic divergence in a system where adaptation from standing genetic variation was constrained.     

The role of environment and core‐margin effects on range‐wide phenotypic variation in a montane grasshopper

The integration of genetic information with ecological and phenotypic data constitutes an effective approach to gain insight into the mechanisms determining interpopulation variability and the evolutionary processes underlying local adaptation and incipient speciation. Here, we use the Pyrenean Morales grasshopper (Chorthippus saulcyi moralesi) as study system to (i) analyse the relative role of genetic drift and selection in range‐wide patterns of phenotypic differentiation and (ii) identify the potential selective agents (environment, elevation) responsible for variation. We also test the hypothesis that (iii) the development of dispersal‐related traits is associated with different parameters related to population persistence/turnover, including habitat suitability stability over the last 120 000 years, distance to the species distribution core and population genetic variability. Our results indicate that selection shaped phenotypic differentiation across all the studied morphological traits (body size, forewing length and shape). Subsequent analyses revealed that among‐population differentiation in forewing length was significantly explained by a temperature gradient, suggesting an adaptive response to thermoregulation or flight performance under contrasting temperature regimes. We found support for our hypothesis predicting a positive association between the distance to the species distribution core and the development of dispersal‐related morphology, which suggests an increased dispersal capability in populations located at range edges that, in turn, exhibit lower levels of genetic variability. Overall, our results indicate that range‐wide patterns of phenotypic variation are partially explained by adaptation in response to local environmental conditions and differences in habitat persistence between core and peripheral populations.     

Evidence for morphological and adaptive genetic divergence between lake and stream habitats in European minnows (Phoxinus phoxinus, Cyprinidae)

Natural selection drives local adaptation, potentially even at small temporal and spatial scales. As a result, adaptive genetic and phenotypic divergence can occur among populations living in different habitats. We investigated patterns of differentiation between contrasting lake and stream habitats in the cyprinid fish European minnow (Phoxinus phoxinus) at both the morphological and genomic levels using geometric morphometrics and AFLP markers, respectively. We also used a spatial correlative approach to identify AFLP loci associated with environmental variables representing potential selective forces responsible for adaptation to divergent habitats. Our results identified different morphologies between lakes and streams, with lake fish presenting a deeper body and caudal peduncle compared to stream fish. Body shape variation conformed to a priori predictions concerning biomechanics and swimming performance in lakes vs. streams. Moreover, morphological differentiation was found to be associated with several environmental variables, which could impose selection on body and caudal peduncle shape. We found adaptive genetic divergence between these contrasting habitats in the form of 'outlier' loci (2.9%) whose genetic divergence exceeded neutral expectations. We also detected additional loci (6.6%) not associated with habitat type (lake vs. stream), but contributing to genetic divergence between populations. Specific environmental variables related to trophic dynamics, landscape topography and geography were associated with several neutral and outlier loci. These results provide new insights into the morphological divergence and genetic basis of adaptation to differentiated habitats.     

Adaptive inter-population differences in blue tit life-history traits on Corsica

Few studies of natural populations have investigated how phenotypic variation across populations relates to key factors in the environment and landscape structure. In the blue tits of southern France, inter-population differences in reproductive life-history traits (e.g. laying date and clutch size) are small, whatever the timing of maximum caterpillar availability, a key factor for offspring survival in tits. These small differences are attributed to gene flow between local populations occupying different habitat types. In contrast, in blue tits on the island of Corsica, we noted large differences in reproductive life-history traits between two populations, where each population is synchronized with the peak-date of caterpillar abundance. These occur over a short geographical distance (25km). Considering our study within a framework of long-term population studies in tits, our results support the hypothesis that different blue tit populations on Corsica show adaptive differences in life-history traits, and suggest that landscape structure at a small spatial scale can have profound effects on adaptive between-population differentiation in life-history traits that are closely linked with fitness.     

The population genomic signature of environmental association and gene flow in an ecologically divergent tree species Metrosideros polymorpha (Myrtaceae)

Genomewide markers enable us to study genetic differentiation within a species and the factors underlying it at a much higher resolution than before, which advances our understanding of adaptation in organisms. We investigated genomic divergence in Metrosideros polymorpha, a woody species that occupies a wide range of ecological habitats across the Hawaiian Islands and shows remarkable phenotypic variation. Using 1659 single nucleotide polymorphism (SNP) markers annotated with the genome assembly, we examined the population genetic structure and demographic history of nine populations across five elevations and two ages of substrates on Mauna Loa, the island of Hawaii. The nine populations were differentiated into two genetic clusters distributed on the lower and higher elevations and were largely admixed on the middle elevation. Demographic modelling revealed that the two genetic clusters have been maintained in the face of gene flow, and the effective population size of the high‐altitude cluster was much smaller. A F_ST‐based outlier search among the 1659 SNPs revealed that 34 SNPs (2.05%) were likely to be under divergent selection and the allele frequencies of 21 of them were associated with environmental changes along elevations, such as temperature and precipitation. This study shows a genomic mosaic of M. polymorpha, in which contrasting divergence patterns were found. While most genomic polymorphisms were shared among populations, a small fraction of the genome was significantly differentiated between populations in diverse environments and could be responsible for the dramatic adaptation to a wide range of environments.     

Population genomics reveals a candidate gene involved in bumble bee pigmentation

Variation in bumble bee color patterns is well‐documented within and between species. Identifying the genetic mechanisms underlying such variation may be useful in revealing evolutionary forces shaping rapid phenotypic diversification. The widespread North American species Bombus bifarius exhibits regional variation in abdominal color forms, ranging from red‐banded to black‐banded phenotypes and including geographically and phenotypically intermediate forms. Identifying genomic regions linked to this variation has been complicated by strong, near species level, genome‐wide differentiation between red‐ and black‐banded forms. Here, we instead focus on the closely related black‐banded and intermediate forms that both belong to the subspecies B. bifarius nearcticus. We analyze an RNA sequencing (RNAseq) data set and identify a cluster of single nucleotide polymorphisms (SNPs) within one gene, Xanthine dehydrogenase/oxidase‐like, that exhibit highly unusual differentiation compared to the rest of the sequenced genome. Homologs of this gene contribute to pigmentation in other insects, and results thus represent a strong candidate for investigating the genetic basis of pigment variation in B. bifarius and other bumble bee mimicry complexes.     

On humans and wildlife in Mediterranean islands

Aim To investigate the effects of human‐induced landscape changes in Mediterranean islands on the ecological and evolutionary responses of bird communities and populations. The combination of mass extinction of large mammals and massive deforestation by humans was hypothesized to produce new selection regimes to which organisms were likely to respond. Habitat selection and niche breadth have been investigated at the scale of species, and phenotypic variation at the scale of local populations. Location The study was carried out along habitat gradients and in habitat mosaics at different spatial scales on the island of Corsica and in areas of similar size and structure in continental France. Methods Two sets of gradients have been used for investigating habitat selection and niche breadth: gradients of altitude, and gradients of vegetation structure. Population studies focused on the blue tit, Cyanistes caeruleus. Large samples of breeding attempts by this species in 10 habitats provided detailed data on phenotypic variation of fitness‐related traits both on Corsica and on the mainland. Results The extent of niche space used by birds differed substantially depending on which habitat gradient was considered. Many species have been found to contract their habitat niche along the elevation gradient on Corsica compared with the mainland, whereas all species in the vegetation gradient broadened their niche on the island. Breeding patterns of the blue tit differed considerably depending on whether they settle in deciduous oaks (Quercus humilis) or in evergreen sclerophyllous oaks (Quercus ilex). Phenotypic variation of breeding traits was much higher on the island, where more populations were correctly timed for the best breeding period than on the mainland, a pattern that is likely to result from lower dispersal of organisms on the island. Main conclusions The differences in observed niche breadth between the two series of habitat gradients is explained both by the species‐specific ecology of the species and the human‐induced environmental history of Corsica. Large‐scale landscape changes provided new opportunities for island colonization by non‐forest species, which are isolated as small, ‘fugitive’ local populations. In both gradients, forest species that are typical components of the Corsican bird fauna definitely expanded their niche and occupied a wider range of habitats on Corsica than on the mainland. At the population scale, landscapes included habitat patches with contrasted selection regimes, which resulted in high phenotypic variation for many fitness‐related traits. Reduced dispersal of birds on the island resulted in a much higher degree of local differentiation on Corsica than on the mainland.     

Nonrandom dispersal drives phenotypic divergence within a bird population

Gene flow through dispersal has traditionally been thought to function as a force opposing evolutionary differentiation. However, directional gene flow may actually reinforce divergence of populations in close proximity. This study documents the phenotypic differentiation over more than two decades in body size (tarsus length) at a very short spatial scale (1.1 km) within a population of pied flycatchers Ficedula hypoleuca inhabiting deciduous and coniferous habitats. Unlike females, males breeding in the deciduous forest were consistently larger than those from the managed coniferous forest. This assortment by size is likely explained by preset habitat preferences leading to dominance of the largest males and exclusion of the smallest ones toward the nonpreferred coniferous forest coupled with directional dispersal. Movements of males between forests were nonrandom with respect to body size and flow rate, which might function to maintain the phenotypic variation in this heritable trait at such a small spatial scale. However, a deeply rooted preference for the deciduous habitat might not be in line with its quality due to the increased levels of breeding density of hole‐nesting competitors therein. These results illustrate how eco‐evolutionary scenarios can develop under directional gene flow over surprisingly small spatial scales. Our findings come on top of recent studies concerning new ways in which dispersal and gene flow can influence microevolution.     

Fine‐grained adaptive divergence in an amphibian: genetic basis of phenotypic divergence and the role of nonrandom gene flow in restricting effective migration among wetlands

Adaptive ecological differentiation among sympatric populations is promoted by environmental heterogeneity, strong local selection and restricted gene flow. High gene flow, on the other hand, is expected to homogenize genetic variation among populations and therefore prevent local adaptation. Understanding how local adaptation can persist at the spatial scale at which gene flow occurs has remained an elusive goal, especially for wild vertebrate populations. Here, we explore the roles of natural selection and nonrandom gene flow (isolation by breeding time and habitat choice) in restricting effective migration among local populations and promoting generalized genetic barriers to neutral gene flow. We examined these processes in a network of 17 breeding ponds of the moor frog Rana arvalis, by combining environmental field data, a common garden experiment and data on variation in neutral microsatellite loci and in a thyroid hormone receptor (TRβ) gene putatively under selection. We illustrate the connection between genotype, phenotype and habitat variation and demonstrate that the strong differences in larval life history traits observed in the common garden experiment can result from adaptation to local pond characteristics. Remarkably, we found that haplotype variation in the TRβ gene contributes to variation in larval development time and growth rate, indicating that polymorphism in the TRβ gene is linked with the phenotypic variation among the environments. Genetic distance in neutral markers was correlated with differences in breeding time and environmental differences among the ponds, but not with geographical distance. These results demonstrate that while our study area did not exceed the scale of gene flow, ecological barriers constrained gene flow among contrasting habitats. Our results highlight the roles of strong selection and nonrandom gene flow created by phenological variation and, possibly, habitat preferences, which together maintain genetic and phenotypic divergence at a fine‐grained spatial scale.