Novel trophic niches drive variable progress towards ecological speciation within an adaptive radiation of pupfishes

Adaptive radiation is recognized by a rapid burst of phenotypic, ecological and species diversification. However, it is unknown whether different species within an adaptive radiation evolve reproductive isolation at different rates. We compared patterns of genetic differentiation between nascent species within an adaptive radiation of Cyprinodon pupfishes using genotyping by sequencing. Similar to classic adaptive radiations, this clade exhibits rapid morphological diversification rates and two species are novel trophic specialists, a scale‐eater and hard‐shelled prey specialist (durophage), yet the radiation is <10 000 years old. Both specialists and an abundant generalist species all coexist in the benthic zone of lakes on San Salvador Island, Bahamas. Based on 13 912 single‐nucleotide polymorphisms (SNPs), we found consistent differences in genetic differentiation between each specialist species and the generalist across seven lakes. The scale‐eater showed the greatest genetic differentiation and clustered by species across lakes, whereas durophage populations often clustered with sympatric generalist populations, consistent with parallel speciation across lakes. However, we found strong evidence of admixture between durophage populations in different lakes, supporting a single origin of this species and genome‐wide introgression with sympatric generalist populations. We conclude that the scale‐eater is further along the speciation‐with‐gene‐flow continuum than the durophage and suggest that different adaptive landscapes underlying these two niche environments drive variable progress towards speciation within the same habitat. Our previous measurements of fitness surfaces in these lakes support this conclusion: the scale‐eating fitness peak may be more distant than the durophage peak on the complex adaptive landscape driving adaptive radiation.     

The genetic architecture of novel trophic specialists: larger effect sizes are associated with exceptional oral jaw diversification in a pupfish adaptive radiation

The genetic architecture of adaptation is fundamental to understanding the mechanisms and constraints governing diversification. However, most case studies focus on loss of complex traits or parallel speciation in similar environments. It is still unclear how the genetic architecture of these local adaptive processes compares to the architecture of evolutionary transitions contributing to morphological and ecological novelty. Here, we identify quantitative trait loci (QTL) between two trophic specialists in an excellent case study for examining the origins of ecological novelty: a sympatric radiation of pupfishes endemic to San Salvador Island, Bahamas, containing a large‐jawed scale‐eater and a short‐jawed molluscivore with a skeletal nasal protrusion. These specialized niches and trophic traits are unique among over 2000 related species. Measurements of the fitness landscape on San Salvador demonstrate multiple fitness peaks and a larger fitness valley isolating the scale‐eater from the putative ancestral intermediate phenotype of the generalist, suggesting that more large‐effect QTL should contribute to its unique phenotype. We evaluated this prediction using an F2 intercross between these specialists. We present the first linkage map for pupfishes and detect significant QTL for sex and eight skeletal traits. Large‐effect QTL contributed more to enlarged scale‐eater jaws than the molluscivore nasal protrusion, consistent with predictions from the adaptive landscape. The microevolutionary genetic architecture of large‐effect QTL for oral jaws parallels the exceptional diversification rates of oral jaws within the San Salvador radiation observed over macroevolutionary timescales and may have facilitated exceptional trophic novelty in this system.     

Context dependence in complex adaptive landscapes: frequency and trait‐dependent selection surfaces within an adaptive radiation of Caribbean pupfishes

The adaptive landscape provides the foundational bridge between micro‐ and macroevolution. One well‐known caveat to this perspective is that fitness surfaces depend on ecological context, including competitor frequency, traits measured, and resource abundance. However, this view is based largely on intraspecific studies. It is still unknown how context‐dependence affects the larger features of peaks and valleys on the landscape which ultimately drive speciation and adaptive radiation. Here, I explore this question using one of the most complex fitness landscapes measured in the wild in a sympatric pupfish radiation endemic to San Salvador Island, Bahamas by tracking survival and growth of laboratory‐reared F2 hybrids. I present new analyses of the effects of competitor frequency, dietary isotopes, and trait subsets on this fitness landscape. Contrary to expectations, decreasing competitor frequency increased survival only among very common phenotypes, whereas less common phenotypes rarely survived despite few competitors, suggesting that performance, not competitor frequency, shapes large‐scale features of the fitness landscape. Dietary isotopes were weakly correlated with phenotype and growth, but did not explain additional survival variation. Nonlinear fitness surfaces varied substantially among trait subsets, revealing one‐, two‐, and three‐peak landscapes, demonstrating the complexity of selection in the wild, even among similar functional traits.     

Mate Choice by Both Sexes Maintains Reproductive Isolation in a Species Flock of Pupfish (Cyprinodon spp) in the Bahamas

Female and male mate choices can reinforce reproductive isolation after sympatric speciation. Using a binary choice design, we examine the importance of visual cues in female mate choice in all three sympatric species of pupfish on San Salvador Island. We also examine the importance of olfactory cues in female choice of the hard‐shelled invertebrate specialist (Cyprinodon brontotheroides). We examine male mate choice in two of the three species, the scale eater (C. desquamator) and the detritivore (C. variegatus). Females of all three species use visual cues and prefer conspecific males. C. brontotheroides females do not use olfactory cues to discriminate between conspecific and heterospecific males. Males of C. desquamator and C. variegatus also preferentially court conspecific females. Thus, mutual mate choice, where both females and males exhibit mate choice, acts as a strong behavioral pre‐mating isolation mechanism in these sympatrically speciated pupfish.     

Male competition fitness landscapes predict both forward and reverse speciation

Speciation is facilitated when selection generates a rugged fitness landscape such that populations occupy different peaks separated by valleys. Competition for food resources is a strong ecological force that can generate such divergent selection. However, it is unclear whether intrasexual competition over resources that provide mating opportunities can generate rugged fitness landscapes that foster speciation. Here we use highly variable male F2 hybrids of benthic and limnetic threespine sticklebacks, Gasterosteus aculeatus Linnaeus, 1758, to quantify the male competition fitness landscape. We find that disruptive sexual selection generates two fitness peaks corresponding closely to the male phenotypes of the two parental species, favouring divergence. Most surprisingly, an additional region of high fitness favours novel hybrid phenotypes that correspond to those observed in a recent case of reverse speciation after anthropogenic disturbance. Our results reveal that sexual selection through male competition plays an integral role in both forward and reverse speciation.     

The Urban Heat Island and its spatial scale dependent impact on survival and development in butterflies of different thermal sensitivity

Climate alteration is one of the most cited ecological consequences of urbanization. However, the magnitude of this impact is likely to vary with spatial scale. We investigated how this alteration affects the biological fitness of insects, which are especially sensitive to ambient conditions and well‐suited organisms to study urbanization‐related changes in phenotypic traits. We monitored temperature and relative air humidity in wooded sites characterized by different levels of urbanization in the surroundings. Using a split‐brood design experiment, we investigated the effect of urbanization at the local (i.e., 200 × 200 m) and landscape (i.e., 3 × 3 km) scale on two key traits of biological fitness in two closely related butterfly species that differ in thermal sensitivity. In line with the Urban Heat Island concept, urbanization led to a 1°C increase in daytime temperature and an 8% decrease in daytime relative humidity at the local scale. The thermophilous species Lasiommata megera responded at the local scale: larval survival increased twofold in urban compared to rural sites. Urbanized sites tended to produce bigger adults, although this was the case for males only. In the woodland species Pararge aegeria, which has recently expanded its ecological niche, we did not observe such a response, neither at the local, nor at the landscape scale. These results demonstrate interspecific differences in urbanization‐related phenotypic plasticity and larval survival. We discuss larval pre‐adaptations in species of different ecological profiles to urban conditions. Our results also highlight the significance of considering fine‐grained spatial scales in urban ecology.     

路径依赖下的物种形成机制

生物科学几乎所有研究都需要物种概念作为基础, 生物多样性研究亦需要可操作的物种概念, 但现有物种概念存在不同程度的人为因素或难操作性, 对物种划分造成不利影响。本文引入“进化路径”这一概念, 说明适合度景观时刻变化着, 物种在每个进化时间点上依据瞬时适合度选择下一时刻的进化状态, 且总是沿着动态适合度景观中适合度增加的方向进化。基于演化博弈的方法, 以随机过程为例模拟物种的进化过程。进而提出路径依赖下的物种形成机制, 并在此基础上给出可操作的物种定义, 即: 针对基因、性状、生态过程等任一状态下两个群体内个体的多个变量做统计分析, 若群体之间同时在两个或多个维度状态下呈现出的不连续性d大于群体内变量呈现出的差异性σ_k, 则拥有相应变量的个体属于不同物种。     

Urbanization drives community shifts towards thermophilic and dispersive species at local and landscape scales

The increasing conversion of agricultural and natural areas to human‐dominated urban landscapes is predicted to lead to a major decline in biodiversity worldwide. Two conditions that typically differ between urban environments and the surrounding landscape are increased temperature, and high patch isolation and habitat turnover rates. However, the extent and spatial scale at which these altered conditions shape biotic communities through selection and/or filtering on species traits are currently poorly understood. We sampled carabid beetles at 81 sites in Belgium using a hierarchically nested sampling design wherein three local‐scale (200 × 200 m) urbanization levels were repeatedly sampled across three landscape‐scale (3 × 3 km) urbanization levels. First, we showed that communities sampled in the most urbanized locations and landscapes displayed a distinct species composition at both local and landscape scale. Second, we related community means of species‐specific thermal preferences and dispersal capacity (based on European distribution and wing morphology, respectively) to the urbanization gradients. We showed that urban communities consisted on average of species with a preference for higher temperatures and with better dispersal capacities compared to rural communities. These shifts were caused by an increased number of species tolerating higher temperatures, a decreased richness of species with low thermal preference, and an almost complete depletion of species with very low‐dispersal capacity in the most urbanized localities. Effects of urbanization were most clearly detected at the local scale, although more subtle effects could also be found at the scale of entire landscapes. Our results demonstrate that urbanization may fundamentally and consistently alter species composition by exerting a strong filtering effect on species dispersal characteristics and favouring replacement by warm‐dwelling species.     

The rank ordering of genotypic fitness values predicts genetic constraint on natural selection on landscapes lacking sign epistasis

Sewall Wright's genotypic fitness landscape makes explicit one mechanism by which epistasis for fitness can constrain evolution by natural selection. Wright distinguished between landscapes possessing multiple fitness peaks and those with only a single peak and emphasized that the former class imposes substantially greater constraint on natural selection. Here I present novel formalism that more finely partitions the universe of possible fitness landscapes on the basis of the rank ordering of their genotypic fitness values. In this report I focus on fitness landscapes lacking sign epistasis (i.e., landscapes that lack mutations the sign of whose fitness effect varies epistatically), which constitute a subset of Wright's single peaked landscapes. More than one fitness rank ordering lacking sign epistasis exists for L > 2 (where L is the number of interacting loci), and I find that a highly statistically significant effect exists between landscape membership in fitness rank-ordering partition and two different proxies for genetic constraint, even within this subset of landscapes. This statistical association is robust to population size, permitting general inferences about some of the characteristics of fitness rank orderings responsible for genetic constraint on natural selection.     

Multilevel landscape utilization of the Siberian flying squirrel: Scale effects on species habitat use

Animals use and select habitat at multiple hierarchical levels and at different spatial scales within each level. Still, there is little knowledge on the scale effects at different spatial levels of species occupancy patterns. The objective of this study was to examine nonlinear effects and optimal‐scale landscape characteristics that affect occupancy of the Siberian flying squirrel, Pteromys volans, in South‐ and Mid‐Finland. We used presence–absence data (n = 10,032 plots of 9 ha) and novel approach to separate the effects on site‐, landscape‐, and regional‐level occupancy patterns. Our main results were: landscape variables predicted the placement of population patches at least twice as well as they predicted the occupancy of particular sites; the clear optimal value of preferred habitat cover for species landscape‐level abundance is a surprisingly low value (10% within a 4 km buffer); landscape metrics exert different effects on species occupancy and abundance in high versus low population density regions of our study area. We conclude that knowledge of regional variation in landscape utilization will be essential for successful conservation of the species. The results also support the view that large‐scale landscape variables have high predictive power in explaining species abundance. Our study demonstrates the complex response of species occurrence at different levels of population configuration on landscape structure. The study also highlights the need for data in large spatial scale to increase the precision of biodiversity mapping and prediction of future trends.     

Fine‐scale geographic patterns of gene flow and reproductive character displacement in Drosophila subquinaria and Drosophila recens

When two species are incompletely isolated, strengthening premating isolation barriers in response to the production of low fitness hybrids may complete the speciation process. Here, we use the sister species Drosophila subquinaria and Drosophila recens to study the conditions under which this reinforcement of species boundaries occurs in natural populations. We first extend the region of known sympatry between these species, and then we conduct a fine‐scale geographic survey of mate discrimination coupled with estimates of gene flow within and admixture between species. Within D. subquinaria, reinforcement is extremely effective: we find variation in mate discrimination both against D. recens males and against conspecific allopatric males on the scale of a few kilometres and in the face of gene flow both from conspecific populations and introgression from D. recens. In D. recens, we do not find evidence for increased mate discrimination in sympatry, even where D. recens is rare, consistent with substantial gene flow throughout the species’ range. Finally, we find that introgression between species is asymmetric, with more from D. recens into D. subquinaria than vice versa. Within each species, admixture is highest in the geographic region where it is rare relative to the other species, suggesting that when hybrids are produced they are of low fitness. In sum, reinforcement within D. subquinaria is effective at maintaining species boundaries, but even when reinforcing selection is strong it may not always result in a pattern of strong reproductive character displacement due to variation in the frequency of hybridization and gene flow from neighbouring populations.     

Large moths captures by a pest monitoring system depend on farmland heterogeneity

As intensive farmlands cover increasing areas of the world, associated biotic richness is crucial for the biodiversity of entire regions. Using data on non‐target Macrolepidopteran moths captured by a crop pest monitoring system, we compared local (100 m perimeter)‐ and landscape‐scale (1000 m perimeter) predictors of the numbers of moth individuals and moth species richness. During a single year (2009), eighteen light traps captured 91 726 individuals of 564 moths species. Typically for biotically impoverished habitats, the catches were dominated by a few superabundant species. Even in these impoverished assemblages, numbers of species increased with increasing herb and woody plants diversity (100 m around the traps), crop diversity (1000 m perimeter), landscape composition (1000 m) and configuration (100 and 1000 m). Abundance of the catches increased only with woody plants diversity in 100 m perimeters. In separate analyses of two species‐rich families, the presumably less mobile Geometridae increased with landscape configuration (i.e. density of edges) within 100 m perimeters around the traps, whereas the species richness of more mobile Noctuidae also reflected the landscape composition (i.e. proportional representation of land covers) at 100 m perimeters. Proportional representation of pest species decreased with increasing richness of herb and woody plants. Taken together, farmland heterogeneity increases moth species richness, whereas abundance of the catches mainly depends on local factors in the vicinity of light traps, and the local factors are more important for presumably less mobile Geometridae than for more mobile Noctuidae.     

Speciation and dietary specialization in Drupa,a genus of predatory marine snails (Gastropoda: Muricidae)

Claremont, M., Reid, D.G. & Williams, S.T. (2012) Speciation and dietary specialization in Drupa, a genus of predatory marine snails (Gastropoda: Muricidae). —Zoologica Scripta, 41, 137–149. We test the competing predictions of allopatric speciation and of ecological speciation by dietary specialization in Drupa, an Indo‐Pacific genus of carnivorous marine gastropods in the family Muricidae. We use a well‐resolved molecular phylogeny (reconstructed from one nuclear and two mitochondrial genes) to show the validity of the traditional species D. elegans, D. rubusidaeus, D. clathrata, D. morum and D. speciosa.‘Drupa ricinus’ is shown to consist of three species: D. ricinus s. s., D. albolabris and a new species, possibly endemic to Japan. ‘Purpura’aperta is transferred to Drupa. Despite potential widespread dispersal and a high degree of range overlap among sister species, range sizes between sister species are highly asymmetric, suggesting that speciation has been predominately peripatric. The exception is the sister pair D. ricinus s. s. and D. albolabris, which have symmetric range sizes and are sympatric over broad Indo‐Pacific ranges. Such symmetry and extensive sympatry are contrary to the predictions of the (peripatric) allopatric model of speciation. Nevertheless, contrary to the predictions of an ecological speciation model based upon dietary specialization, broad dietary range appears to be identical between the species. Small differences in microhabitat preferences (or hypothetical dietary specialization at a fine taxonomic scale) may have been significant in the speciation process or, if initial divergence was allopatric, in permitting subsequent sympatry. Broad dietary shifts appear to have accompanied more ancient divergences within the genus Drupa.     

Modeling the two‐locus architecture of divergent pollinator adaptation: how variation in SAD paralogs affects fitness and evolutionary divergence in sexually deceptive orchids

Divergent selection by pollinators can bring about strong reproductive isolation via changes at few genes of large effect. This has recently been demonstrated in sexually deceptive orchids, where studies (1) quantified the strength of reproductive isolation in the field; (2) identified genes that appear to be causal for reproductive isolation; and (3) demonstrated selection by analysis of natural variation in gene sequence and expression. In a group of closely related Ophrys orchids, specific floral scent components, namely n‐alkenes, are the key floral traits that control specific pollinator attraction by chemical mimicry of insect sex pheromones. The genetic basis of species‐specific differences in alkene production mainly lies in two biosynthetic genes encoding stearoyl–acyl carrier protein desaturases (SAD) that are associated with floral scent variation and reproductive isolation between closely related species, and evolve under pollinator‐mediated selection. However, the implications of this genetic architecture of key floral traits on the evolutionary processes of pollinator adaptation and speciation in this plant group remain unclear. Here, we expand on these recent findings to model scenarios of adaptive evolutionary change at SAD2 and SAD5, their effects on plant fitness (i.e., offspring number), and the dynamics of speciation. Our model suggests that the two‐locus architecture of reproductive isolation allows for rapid sympatric speciation by pollinator shift; however, the likelihood of such pollinator‐mediated speciation is asymmetric between the two orchid species O. sphegodes and O. exaltata due to different fitness effects of their predominant SAD2 and SAD5 alleles. Our study not only provides insight into pollinator adaptation and speciation mechanisms of sexually deceptive orchids but also demonstrates the power of applying a modeling approach to the study of pollinator‐driven ecological speciation.     

Scale dependent drivers of wild bee diversity in tropical heterogeneous agricultural landscapes

Factors associated with agricultural intensification, for example, loss of seminatural vegetation and pesticide use has been shown to adversely affect the bee community. These factors may impact the bee community differently at different landscape scales. The scale dependency is expected to be more pronounced in heterogeneous landscapes. However, the scale‐dependent response of the bee community to drivers of its decline is relatively understudied, especially in the tropics where the agricultural landscape is often heterogeneous. This study looked at effects of agricultural intensification on bee diversity at patch and landscape scales in a tropical agricultural landscape. Wild bees were sampled using 12 permanent pan trap stations. Patch and landscape characteristics were measured within a 100 m (patch scale) and a 500 m (landscape scale) radius of pan trap stations. Information on pesticide input was obtained from farmer surveys. Data on vegetation cover, productivity, and percentage of agricultural and fallow land (FL) were collected using satellite imagery. Intensive areas in a bee‐site network were less specialized in terms of resources to attract rare bee species while the less intensive areas, which supported more rare species, were more vulnerable to disturbance. A combination of patch quality and diversity as well as pesticide use regulates species diversity at the landscape scale (500 m), whereas pesticide quantity drove diversity at the patch scale (100 m). At the landscape scale, specialization of each site in terms of resources for bees increased with increasing patch diversity and FL while at the patch scale specialization declined with increased pesticide use. Bee functional groups responded differentially to landscape characteristics as well as pesticide use. Wood nesting bees were negatively affected by the number of pesticides used but other bee functional groups were not sensitive to pesticides. Synthesis and Applications: Different factors affect wild bee diversity at the scale of landscape and patch in heterogeneous tropical agricultural systems. The differential response of bee functional groups to agricultural intensification underpins the need for guild‐specific management strategies for wild bee conservation. Less intensively farmed areas support more rare species and are vulnerable to disturbance; consequently, these areas should be prioritized for conservation to maintain heterogeneity in the landscape. It is important to conserve and restore seminatural habitats to maintain complexity in the landscapes through participatory processes and to regulate synthetic chemical pesticides in farm operations to conserve the species and functional diversity of wild bees.     

Ecological divergence and speciation between lemur (Eulemur) sister species in Madagascar

Understanding ecological niche evolution over evolutionary timescales is crucial to elucidating the biogeographic history of organisms. Here, we used, for the first time, climate‐based ecological niche models (ENMs) to test hypotheses about ecological divergence and speciation processes between sister species pairs of lemurs (genus Eulemur) in Madagascar. We produced ENMs for eight species, all of which had significant validation support. Among the four sister species pairs, we found nonequivalent niches between sisters, varying degrees of niche overlap in ecological and geographic space, and support for multiple divergence processes. Specifically, three sister‐pair comparisons supported the null model that niches are no more divergent than the available background region. These findings are consistent with an allopatric speciation model, and for two sister pairs (E. collaris–E. cinereiceps and E. rufus–E. rufifrons), a riverine barrier has been previously proposed for driving allopatric speciation. However, for the fourth sister pair E. flavifrons–E. macaco, we found support for significant niche divergence, and consistent with their parapatric distribution on an ecotone and the lack of obvious geographic barriers, these findings most strongly support a parapatric model of speciation. These analyses thus suggest that various speciation processes have led to diversification among closely related Eulemur species.     

Evolution of interspecies unilateral incompatibility in the relatives of Arabidopsis thaliana

The evolutionary concurrence of intraspecies self‐incompatibility (SI) and explosive angiosperm radiation in the Cretaceous have led to the hypothesis that SI was one of the predominant drivers of rapid speciation in angiosperms. Interspecies unilateral incompatibility (UI) usually occurs when pollen from a self‐compatible (SC) species is rejected by the pistils of a SI species, while the reciprocal pollination is compatible (UC). Although this SI × SC type UI is most prevalent and viewed as a prezygotic isolation barrier to promote incipient speciation of angiosperms, comparative evidence to support such a role is lacking. We show that SI × SI type UI in SI species pairs is also common in the well‐characterized accessions representing the four major lineages of the Arabidopsis genus and is developmentally regulated. This allowed us to reveal a strong correlation between UI strength and species divergence in these representative accessions. In addition, analyses of a SC accession and the pseudo‐self‐compatible (psc) spontaneous mutant of Arabidopsis lyrata indicate that UI shares, at least, common pollen rejection pathway with SI. Furthermore, genetic and genomic analyses of SI × SI type UI in A. lyrata × A. arenosa species pair showed that two major‐effect quantitative trait loci are the stigma and pollen‐side determinant of UI, respectively, which could be involved in heterospecies pollen discrimination. By revealing a close link between UI and SI pathway, particularly between UI and species divergence in these representative accessions, our findings establish a connection between SI and speciation. Thus, the pre‐existence of SI system would have facilitated the evolution of UI and accordingly promote speciation.     

Differential responses of cryptic bat species to the urban landscape

Urbanization is a key global driver in the modification of land use and has been linked to population declines even in widespread and relatively common species. Cities comprise a complex assortment of habitat types yet we know relatively little about the effects of their composition and spatial configuration on species distribution. Although many bat species exploit human resources, the majority of species are negatively impacted by urbanization. Here, we use data from the National Bat Monitoring Programme, a long‐running citizen science scheme, to assess how two cryptic European bat species respond to the urban landscape. A total of 124 × 1 km² sites throughout Britain were surveyed. The landscape surrounding each site was mapped and classified into discrete biotope types (e.g., woodland). Generalized linear models were used to assess differences in the response to the urban environment between the two species, and which landscape factors were associated with the distributions of P. pipistrellus and P. pygmaeus. The relative prevalence of P. pygmaeus compared to P. pipistrellus was greater in urban landscapes with a higher density of rivers and lakes, whereas P. pipistrellus was frequently detected in landscapes comprising a high proportion of green space (e.g., parklands). Although P. pipistrellus is thought to be well adapted to the urban landscape, we found a strong negative response to urbanization at a relatively local scale (1 km), whilst P. pygmaeus was detected more regularly in wooded urban landscapes containing freshwater. These results show differential habitat use at a landscape scale of two morphologically similar species, indicating that cryptic species may respond differently to anthropogenic disturbance. Even species considered relatively common and well adapted to the urban landscape may respond negatively to the built environment highlighting the future challenges involved in maintaining biodiversity within an increasingly urbanized world.     

Spatial and ecological population genetic structures within two island‐endemic Aeonium species of different niche width

The Crassulacean genus Aeonium is a well‐known example for plant species radiation on oceanic archipelagos. However, while allopatric speciation among islands is documented for this genus, the role of intra‐island speciation due to population divergence by topographical isolation or ecological heterogeneity has not yet been addressed. The aim of this study was to investigate intraspecific genetic structures and to identify spatial and ecological drivers of genetic population differentiation on the island scale. We analyzed inter simple sequence repeat variation within two island‐endemic Aeonium species of La Palma: one widespread generalist that covers a large variety of different habitat types (Ae. davidbramwellii) and one narrow ecological specialist (Ae. nobile), in order to assess evolutionary potentials on this island. Gene pool differentiation and genetic diversity patterns were associated with major landscape structures in both species, with phylogeographic implications. However, overall levels of genetic differentiation were low. For the generalist species, outlier loci detection and loci–environment correlation approaches indicated moderate signatures of divergent selection pressures linked to temperature and precipitation variables, while the specialist species missed such patterns. Our data point to incipient differentiation among populations, emphasizing that ecological heterogeneity and topographical structuring within the small scales of an island can foster evolutionary processes. Very likely, such processes have contributed to the radiation of Aeonium on the Canary Islands. There is also support for different evolutionary mechanisms between generalist and specialist species.