FREQUENCY-DEPENDENT SEED PRODUCTION AND HYBRIDIZATION RATES: IMPLICATIONS FOR GENE FLOW BETWEEN LOCALLY ADAPTED PLANT POPULATIONS

Gene flow between genetically distinct plant populations can have significant evolutionary consequences. It can increase genetic diversity, create novel gene combinations, and transfer adaptations from one population to another. This study addresses the roles of frequency-dependent selection and mating system in gene exchange between two subspecies of Gilia capitata (Polemoniaceae). Long-distance migrants are likely to be rare in new habitats, and the importance of immigrant frequency to fitness, gene exchange, and ultimately introgression, has not been explored. To test for the importance of frequency in migration, a field experiment was conducted in which artificial populations (arrays) composed of different mixtures of the two subspecies were placed in the home habitats of both. Female function (seed production) and a portion of male function (hybridization rate) were compared for the two subspecies to assess the potential for gene exchange and introgression between them. Individual fitness (through both hybridization and seed production) for the inland subspecies varied with its frequency as an immigrant at the coastal site. Rare immigrants produced fewer seeds and fathered fewer hybrid offspring. In contrast, both forms of reproductive function were frequency independent for the coastal subspecies when it was an immigrant at the inland site. Seed production was high and insensitive to frequency, and immigrants from the coast never successfully fertilized the inland subspecies' seeds. To control for the effects of frequency-dependent pollinator behavior in the field, hand crosses were performed in the greenhouse using a range of pollen mixtures. The greenhouse experiment demonstrated that cross-fertilization is possible in only one direction, that cross-pollination in the other direction is only partially successful, and that pollen from the coastal subspecies has a strong negative effect on seed production by the inland subspecies. Experimental pollen supplementation in the field verified both the unilateral incompatibility and the negative effect of coastal pollen on inland plant seed production observed in the greenhouse. Contrasts between field array and greenhouse results suggest that pollinator behavior and other ecological factors act to exaggerate reproductive barriers between the two subspecies. In this system, immigrant frequency interacts with reproductive biology and pollinator ecology to enhance gene flow between the populations in one direction, while restricting gene establishment and introgression in the other direction.     

LOCAL ADAPTATION,PHENOTYPIC DIFFERENTIATION,AND HYBRID FITNESS IN DIVERGED NATURAL POPULATIONS OF ARABIDOPSIS LYRATA

Selection for local adaptation results in genetic differentiation in ecologically important traits. In a perennial, outcrossing model plant Arabidopsis lyrata, several differentiated phenotypic traits contribute to local adaptation, as demonstrated by fitness advantage of the local population at each site in reciprocal transplant experiments. Here we compared fitness components, hierarchical total fitness and differentiation in putatively ecologically important traits of plants from two diverged parental populations from different continents in the native climate conditions of the populations in Norway and in North Carolina (NC, U.S.A.). Survival and number of fruits per inflorescence indicated local advantage at both sites and aster life‐history models provided additional evidence for local adaptation also at the level of hierarchical total fitness. Populations were also differentiated in flowering start date and floral display. We also included reciprocal experimental F₁ and F₂ hybrids to examine the genetic basis of adaptation. Surprisingly, the F₂ hybrids showed heterosis at the study site in Norway, likely because of a combination of beneficial dominance effects from different traits. At the NC site, hybrid fitness was mostly intermediate relative to the parental populations. Local cytoplasmic origin was associated with higher fitness, indicating that cytoplasmic genomes also may contribute to the evolution of local adaptation.     

Ecological genetics of Bromus tectorum

Summary By incorporating demographic analyses of fitness components (e.g., survival and reproduction) within a reciprocal sowing design, we tested for 3 consecutive years whether local adaptation has occurred in the alien grass Bromus tectorum (cheatgrass) within 7 habitats along an environmental gradient from arid steppe to subalpine forest in the Intermontain Region of western North America. Patterns of emergence and survival were strongly influenced by the local environment. In terms of survival, expression of significant local adaptation in Tsuga heterophylla habitat varied among years. In contrast, relative differences in flowering time among seed sources were stable across sites and years. Populations from the arid steppe were the earliest to flower; flowering was latest in populations from the mesic Tsuga heterophylla habitat. In terms of net reproductive rate, evidence for local adaptation in B. tectorum was obtained in populations from habitats representing environmental extremes: an arid, saline site dominated by the shrub Sarcobatus vermiculatus and clearings within the cool, mesic Tsuga heterophylla forest habitat. Unlike the plants introduced from other sites, members of the resident population at the Sarcobatus site flowered and produced seeds before soil water became limiting. In contrast, net reproductive rates in other habitats were sometimes the lowest for populations in their home site. This lack of an advantage for local populations within more environmentally moderate sites suggests that limited dispersal may restrict the rate at which superior genotypes are introduced into a particular site.     

Contrasting ecosystem CO2 fluxes of inland and coastal wetlands: a meta‐analysis of eddy covariance data

Wetlands play an important role in regulating the atmospheric carbon dioxide (CO₂) concentrations and thus affecting the climate. However, there is still lack of quantitative evaluation of such a role across different wetland types, especially at the global scale. Here, we conducted a meta‐analysis to compare ecosystem CO₂ fluxes among various types of wetlands using a global database compiled from the literature. This database consists of 143 site‐years of eddy covariance data from 22 inland wetland and 21 coastal wetland sites across the globe. Coastal wetlands had higher annual gross primary productivity (GPP), ecosystem respiration (R_e), and net ecosystem productivity (NEP) than inland wetlands. On a per unit area basis, coastal wetlands provided large CO₂ sinks, while inland wetlands provided small CO₂ sinks or were nearly CO₂ neutral. The annual CO₂ sink strength was 93.15 and 208.37 g C m^?2 for inland and coastal wetlands, respectively. Annual CO₂ fluxes were mainly regulated by mean annual temperature (MAT) and mean annual precipitation (MAP). For coastal and inland wetlands combined, MAT and MAP explained 71%, 54%, and 57% of the variations in GPP, R_e, and NEP, respectively. The CO₂ fluxes of wetlands were also related to leaf area index (LAI). The CO₂ fluxes also varied with water table depth (WTD), although the effects of WTD were not statistically significant. NEP was jointly determined by GPP and R_e for both inland and coastal wetlands. However, the NEP/R_e and NEP/GPP ratios exhibited little variability for inland wetlands and decreased for coastal wetlands with increasing latitude. The contrasting of CO₂ fluxes between inland and coastal wetlands globally can improve our understanding of the roles of wetlands in the global C cycle. Our results also have implications for informing wetland management and climate change policymaking, for example, the efforts being made by international organizations and enterprises to restore coastal wetlands for enhancing blue carbon sinks.     

Local Adaptation and the Effects of Isolation and Population Size – the Semelparous Perennial <Emphasis Type="Italic">Carlina vulgaris</Emphasis> as a Study Case

We made a reciprocal transplantation experiment with Carlina vulgaris, among 12 semi-natural grassland sites situated in southwest Sweden. The fate of seeds and seedlings were followed during 2 years. Local adaptation was investigated both by native superiority over non-natives, and by comparing the observed performance of a population to the fitted value of a reduced statistical model that showed the populations’ performance at all sites and the performance of all other populations at its home site. The latter method indicates presence of local adaptation even when natives are inferior to introduced populations as long as the negative difference in fitness between the populations is smaller at the native population’s home site. The strength of local adaptation was measured as the ratio of the observed to the expected performance in reduced statistical models and regressed on the degree of isolation and population size. We found no evidence of local adaptation in terms of native superiority compared to non-natives, but with the relative method we found one of six fitness components, juvenile survival, to be 6% higher for natives at their home sites compared to how they performed at other sites and how others performed at their site. Further, our results indicate that both isolation and population size have a positive effect on the process of local adaptation. Co-ordinating editor: J. Tuomi     

EFFECTS OF HYBRIDIZATION AND INBREEDING ON FITNESS IN PHLOX

The effects of interpopulation hybridization, and self-fertilization and immigration on fitness in Phlox drummondii were analyzed in 5 natural field sites. Germination rates, survivorship to flowering, fecundity, and net reproductive rates (R_o) were determined for planted populations of natives, hybrids, aliens, and the products of one and two generations of self-fertilization. At all sites, seed germination was 36% for natives, 35% for hybrids and 28% for aliens. Survivorship for natives was 50%, compared to 49% for hybrids and 41% for aliens. The mean fecundity was 34, 35, and 39 seeds per plant for natives, hybrids, and aliens, respectively. The R_o of natives averaged 5.2 vs. 6.4 for hybrids and 4.2 aliens. At all the sites, germination averaged 25% for self-1 plants and 19% for self-2 plants compared to 24% for open-pollinated controls. Survivorship progressed from 50% in the controls to 46% in self-1 and 44% in self-2 plants. Seed production per plant averaged 35% in the control vs. 32% in the self-1 and 25% in the self-2 plants. The mean R_o of the control was 3.8 vs. 3.3 in the self-1 and 2.14 in the self-2 plants. Our results demonstrate that the genetic variable may have a substantial effect on plant fitness in the field.     

Is local adaptation in Mimulus guttatus caused by trade‐offs at individual loci?

Local adaptation is considered to be the result of fitness trade‐offs for particular phenotypes across different habitats. However, it is unclear whether such phenotypic trade‐offs exist at the level of individual genetic loci. Local adaptation could arise from trade‐offs of alternative alleles at individual loci or by complementary sets of loci with different fitness effects of alleles in one habitat but selective neutrality in the alternative habitat. To evaluate the genome‐wide basis of local adaptation, we performed a field‐based quantitative trait locus (QTL) mapping experiment on recombinant inbred lines (RILs) created from coastal perennial and inland annual races of the yellow monkeyflower (Mimulus guttatus) grown reciprocally in native parental habitats. Overall, we detected 19 QTLs affecting one or more of 16 traits measured in two environments, most of small effect. We identified 15 additional QTL effects at two previously identified candidate QTLs [DIV ERGENCE (DIV)]. Significant QTL by environment interactions were detected at the DIV loci, which was largely attributable to genotypic differences at a single field site. We found no detectable evidence for trade‐offs for any one component of fitness, although DIV2 showed a trade‐off involving different fitness traits between sites, suggesting that local adaptation is largely controlled by non‐overlapping loci. This is surprising for an outcrosser, implying that reduced gene flow prevents the evolution of individuals adapted to multiple environments. We also determined that native genotypes were not uniformly adaptive, possibly reflecting fixed mutational load in one of the populations.     

Complex trait divergence contributes to environmental niche differentiation in ecological speciation of Boechera stricta

Ecological factors may contribute to reproductive isolation if differential local adaptation causes immigrant or hybrid fitness reduction. Because local adaptation results from the interaction between natural selection and adaptive traits, it is crucial to investigate both to understand ecological speciation. Previously, we used niche modelling to identify local water availability as an environmental correlate of incipient ecological speciation between two subspecies in Boechera stricta, a close relative of Arabidopsis. Here, we performed several large‐scale glasshouse experiments to investigate the divergence of various physiological, phenological and morphological traits. Although we found no significant difference in physiological traits, the Western subspecies has significantly faster growth rate, larger leaf area, less succulent leaves, delayed reproductive time and longer flowering duration. These trait differences are concordant with previous results that habitats of the Western genotypes have more consistent water availability, while Eastern genotypes inhabit locations with more ephemeral water supplies. In addition, by comparing univariate and multivariate divergence of complex traits (Q_ST) to the genomewide distribution of SNP F_ST, we conclude that the aspects of phenology and morphology (but not physiology) are under divergent selection. In addition, we also identified several highly diverged traits without obvious water‐related functions.     

GENETIC DIFFERENTIATION AND SELECTION AGAINST MIGRANTS IN EVOLUTIONARILY REPLICATED EXTREME ENVIRONMENTS

We investigated mechanisms of reproductive isolation in livebearing fishes (genus Poecilia) inhabiting sulfidic and nonsulfidic habitats in three replicate river drainages. Although sulfide spring fish convergently evolved divergent phenotypes, it was unclear if mechanisms of reproductive isolation also evolved convergently. Using microsatellites, we found strongly reduced gene flow between adjacent populations from different habitat types, suggesting that local adaptation to sulfidic habitats repeatedly caused the emergence of reproductive isolation. Reciprocal translocation experiments indicate strong selection against immigrants into sulfidic waters, but also variation among drainages in the strength of selection against immigrants into nonsulfidic waters. Mate choice experiments revealed the evolution of assortative mating preferences in females from nonsulfidic but not from sulfidic habitats. The inferred strength of sexual selection against immigrants (RI_s) was negatively correlated with the strength of natural selection (RI_m), a pattern that could be attributed to reinforcement, whereby natural selection strengthens behavioral isolation due to reduced hybrid fitness. Overall, reproductive isolation and genetic differentiation appear to be replicated and direct consequences of local adaptation to sulfide spring environments, but the relative contributions of different mechanisms of reproductive isolation vary across these evolutionarily independent replicates, highlighting both convergent and nonconvergent evolutionary trajectories of populations in each drainage.     

Local adaptation and ecological differentiation under selection,migration, and drift in Arabidopsis lyrata*

How the balance between selection, migration, and drift influences the evolution of local adaptation has been under intense theoretical scrutiny. Yet, empirical studies that relate estimates of local adaptation to quantification of gene flow and effective population sizes have been rare. Here, we conducted a reciprocal transplant trial, a common garden trial, and a whole‐genome‐based demography analysis to examine these effects among Arabidopsis lyrata populations from two altitudinal gradients in Norway. Demography simulations indicated that populations within the two gradients are connected by gene flow (0.1 < 4N_em < 11) and have small effective population sizes (N_e < 6000), suggesting that both migration and drift can counteract local selection. However, the three‐year field experiments showed evidence of local adaptation at the level of hierarchical multiyear fitness, attesting to the strength of differential selection. In the lowland habitat, local superiority was associated with greater fecundity, while viability accounted for fitness differences in the alpine habitat. We also demonstrate that flowering time differentiation has contributed to adaptive divergence between these locally adapted populations. Our results show that despite the estimated potential of gene flow and drift to hinder differentiation, selection among these A. lyrata populations has resulted in local adaptation.     

Phytophagous insect oviposition shifts in response to probability of flower abortion owing to the presence of basal fruits

Phytophagous insects use a wide range of indicators or associated cues to avoid laying eggs in sites where offspring survival is low. For insects that lay eggs in flowers, these unsuitable sites may be created by the host plant's resource allocation to flowers. In the sequentially flowering host plant, Yucca glauca, late‐opening distal flowers are more likely to be aborted in the presence of already‐initiated basal fruits because they are strong resource sinks. If flowers are aborted, all eggs of the phytophagous insect, Tegeticula yuccasella, within the flower die. We used the phytophagous insect T. yuccasella that lays eggs in and pollinates host plant Y. glauca flowers to test the hypothesis that phytophagous insect females are less likely to invest eggs in host plant flowers if basal fruits are present because they are more likely to be aborted. We also investigated potential predictors of arrival of T. yuccasella at inflorescences at the onset of flowering. These factors may influence a phytophagous insect's decisions to select oviposition sites. We carried out a behavioral experiment using wild‐caught T. yuccasella females on manipulated inflorescences with distal flowers with basal fruits and without fruits. As potential predictors of T. yuccasella arriving at inflorescences, we used floral display size and day of onset of flowering. In support of our hypothesis, our experimental results showed that T. yuccasella was significantly less likely to oviposit in distal flowers on inflorescences with basal fruits. We also found that T. yuccasella arrival was higher at inflorescences with larger floral display size and earlier in the flowering season. These findings uncover a novel indicator of unsuitable oviposition sites—the presence of basal fruits, that phytophagous insects use to make oviposition decisions. Further, our study contributes to the growing body of evidence that shows that females prefer sites that increase the probability of survival of their offspring.     

Flowering time QTL in natural populations of Arabidopsis thaliana and implications for their adaptive value

The genetic basis of phenotypic traits is of great interest to evolutionary biologists, but their contribution to adaptation in nature is often unknown. To determine the genetic architecture of flowering time in ecologically relevant conditions, we used a recombinant inbred line population created from two locally adapted populations of Arabidopsis thaliana from Sweden and Italy. Using these RILs, we identified flowering time QTL in growth chambers that mimicked the natural temperature and photoperiod variation across the growing season in each native environment. We also compared the genomic locations of flowering time QTL to those of fitness (total fruit number) QTL from a previous three‐year field study. Ten total flowering time QTL were found, and in all cases, the Italy genotype caused early flowering regardless of the conditions. Two QTL were consistent across chamber environments, and these had the largest effects on flowering time. Five of the fitness QTL colocalized with flowering time QTL found in the Italy conditions, and in each case, the local genotype was favoured. In contrast, just two flowering time QTL found in the Sweden conditions colocalized with fitness QTL and in only one case was the local genotype favoured. This implies that flowering time may be more important for adaptation in Italy than Sweden. Two candidate genes (FLC and VIN3) underlying the major flowering time QTL found in the current study are implicated in local adaptation.     

Effects of Competition and Life History Stage on the Expression of Local Adaptation in Two Native Bunchgrasses

Concerns about the use of genetically appropriate material in restoration often focus on questions of local adaptation. Many reciprocal transplant studies have demonstrated local adaptation in native plant species, but very few have examined how interspecific competition affects the expression of adaptive variation. Our study examined regional scales of adaptation between foothill and coastal populations of two California native bunchgrasses (Elymus glaucus and Nassella pulchra). By combining competitive manipulations with reciprocal transplants, we examined the importance of the vegetation at a site as a selective factor in the process of local adaptation. By monitoring survival and reproduction of reciprocally transplanted populations over the course of 3 years, we also studied the effect of life history stage on the expression of local adaptation. For most of the fitness components we measured, local adaptation was detected and interspecific competition consistently amplified its expression. Expression of local adaptation was especially apparent in the more inbreeding species E. glaucus and suggests that with weaker gene flow, selection may be more effective in creating ecotypes within this species. Local adaptation was detected at all life history stages but was most strongly expressed in traits associated with adult reproduction and the viability of seeds produced by the transplants. Taken together, our results indicate that the importance of local adaptation will become more apparent in the later stages of a restoration project as the plants at a site begin to reproduce and as they experience greater interspecific competition from the maturing vegetation at the site.     

Divergent selection on flowering time contributes to local adaptation in Mimulus guttatus populations

The timing of when to initiate reproduction is an important transition in any organism's life cycle. There is much variation in flowering time among populations, but we do not know to what degree this variation contributes to local adaptation. Here we use a reciprocal transplant experiment to examine the presence of divergent natural selection for flowering time and local adaptation between two distinct populations of Mimulus guttatus. We plant both parents and hybrids (to tease apart differences in suites of associated parental traits) between these two populations into each of the two native environments and measure floral, vegetative, life-history, and fitness characters to assess which traits are under selection at each site. Analysis of fitness components indicates that each of these plant populations is locally adapted. We obtain striking evidence for divergent natural selection on date of first flower production at these two sites. Early flowering is favored at the montane site, which is inhabited by annual plants and characterized by dry soils in midsummer, whereas intermediate (though later) flowering dates are selectively favored at the temperate coastal site, which is inhabited by perennial plants and is almost continually moist. Divergent selection on flowering time contributes to local adaptation between these two populations of M. guttatus, suggesting that genetic differentiation in the timing of reproduction may also serve as a partial reproductive isolating barrier to gene flow among populations.     

Patterns of Differentiation in Wiregrass (Aristida beyrichiana): Implications for Restoration Efforts

Aristida beyrichiana (wiregrass) is increasingly being planted in restoration projects across the southeastern coastal plain, with little focus on genetic differences among populations across the region. Local and regional population differentiation for establishment and growth traits were examined in common garden and reciprocal transplant experiments. Seeds from up to 20 plants from each of seven populations were collected in northern and central Florida sites that encompassed gradients of soils, hydrology, and temperature. Reciprocal seed transplants using three of the common garden populations were conducted in two consecutive years. In the common garden, significant population differences were seen in seed weight, seedling emergence and survival, tiller height, number of tillers, the relationship between tiller number and tiller height, and flowering. Variation among maternal families was seen in tiller number and in the relationship between tiller number and tiller height. The reciprocal transplant study did not detect either local adaptation to sites of origin or consistent superiority of one source population or planting site in seedling establishment. These results suggest that the probability of seedling establishment is primarily dependent on environmental conditions rather than genetic differences. Genetic variation for traits related to fitness (e.g., tiller number) may be retained within populations because phenotypically plastic growth responses of seedlings to environmental variation buffer genetic variation against the action of selection. But despite the lack of evidence for genetic influences on initial establishment in wiregrass, our common garden study suggests genetic differences among populations. This result, when combined with previous results indicating local adaptation in later life stages of wiregrass, suggests that restoration efforts involving this species should use local seed sources from sites with similar soil and hydrological conditions.     

Genomic adaptation of flowering‐time genes during the expansion of rice cultivation area

The diversification of flowering time in response to natural environments is critical for the spread of crops to diverse geographic regions. In contrast with recent advances in understanding the molecular basis of photoperiodic flowering in rice (Oryza sativa), little is known about how flowering‐time diversification is structured within rice subspecies. By analyzing genome sequencing data and a set of 429 chromosome segment substitution lines (CSSLs) originating from 10 diverse rice accessions with wide distributions, we revealed diverse effects of allelic variations for common flowering‐time quantitative trait loci in the recipient's background. Although functional variations associated with a few loci corresponded to standing variations among subspecies, the identified functional nucleotide polymorphisms occurred recently after rice subgroup differentiation, indicating that the functional diversity of flowering‐time gene sequences was not particularly associated with phylogenetic relationship between rice subspecies. Intensive analysis of the Hd1 genomic region identified the signature of an early introgression of the Hd1 with key mutation(s) in aus and temperate japonica accessions. Our data suggested that, after such key introgressions, new mutations were selected and accelerated the flowering‐time diversity within subspecies during the expansion of rice cultivation area. This finding may imply that new genome‐wide changes for flowering‐time adaptation are one of the critical determinants for establishing genomic architecture of local rice subgroups. In‐depth analyses of various rice genomes coupling with the genetically confirmed phenotypic changes in a large set of CSSLs enabled us to demonstrate how rice genome dynamics has coordinated with the adaptation of cultivated rice during the expansion of cultivation area.     

On the distribution of the Slow Worm,Anguis fragilis,and the European Glass Lizard,Ophisaurus apodus,in Turkey

Abstract

The two subspecies of Ophisaurus apodus represent different ecotypes: thracicus is a coastal type, which occurs in Turkey along all coasts (but only local along the Black Sea coast), whereas nominate apodus is an inland type which is found in Turkey only in the extreme east. Anguis fragilis is confined to areas in which the monthly mean precipitation is above 20 mm.     

Asymmetric isolating barriers between different microclimatic environments caused by low immigrant survival

Spatially variable selection has the potential to result in local adaptation unless counteracted by gene flow. Therefore, barriers to gene flow will help facilitate divergence between populations that differ in local selection pressures. We performed spatially and temporally replicated reciprocal field transplant experiments between inland and coastal habitats using males of the common blue damselfly (Enallagma cyathigerum) as our study organism. Males from coastal populations had lower local survival rates than resident males at inland sites, whereas we detected no differences between immigrant and resident males at coastal sites, suggesting asymmetric local adaptation in a source–sink system. There were no intrinsic differences in longevity between males from the different environments suggesting that the observed differences in male survival are environment-dependent and probably caused by local adaptation. Furthermore, the coastal environment was found to be warmer and drier than the inland environment, further suggesting local adaptation to microclimatic factors has lead to differential survival of resident and immigrant males. Our results suggest that low survival of immigrant males mediates isolation between closely located populations inhabiting different microclimatic environments.     

Comparing demographic parameters for philopatric and immigrant individuals in a long-lived bird adapted to unstable habitats

Theoretical models about the benefits of philopatry predict that immigrant fitness can be higher, lower or similar to that of philopatrics depending on habitat heterogeneity, dispersal costs, distance between patches or population densities. In this study, we compared transience rates, local survival and recruitment among philopatric and immigrant individuals of Audouin’s gull Larus audouinii, a long-lived bird with high dispersal capacities. Several previous studies have shown that these capacities were probably the result of adaptation to unstable and highly discrete habitats; hence, we tested the hypothesis that fitness components for philopatrics and immigrants were similar. During 1988–2006, ca. 27,800 chicks were marked in 31 colonies in the western Mediterranean metapopulation, and more than 52,000 resightings were made in a single, high-quality colony, to estimate local demographic parameters by capture–recapture analyses. Results suggest that, even though parameters related to site-tenacity (e.g. recapture rates) were higher for philopatrics than for immigrants, survival and recruitment were fundamentally similar. Small differences between philopatrics and immigrants were probably influenced by a highly suitable habitat at the study site, which reduced dispersal costs for immigrants; furthermore, the similarities in most fitness components were also probably the result of a life-history strategy of a species living in unpredictable, unstable habitats with high emigration rates among local populations, and with a relatively low cost of dispersal.     

Potential adaptive divergence between subspecies and populations of snapdragon plants inferred from QST–FST comparisons

Phenotypic divergence among natural populations can be explained by natural selection or by neutral processes such as drift. Many examples in the literature compare putatively neutral (F_ST) and quantitative genetic (Q_ST) differentiation in multiple populations to assess their evolutionary signature and identify candidate traits involved with local adaptation. Investigating these signatures in closely related or recently diversified species has the potential to shed light on the divergence processes acting at the interspecific level. Here, we conducted this comparison in two subspecies of snapdragon plants (eight populations of Antirrhinum majus pseudomajus and five populations of A. m. striatum) in a common garden experiment. We also tested whether altitude was involved with population phenotypic divergence. Our results identified candidate phenological and morphological traits involved with local adaptation. Most of these traits were identified in one subspecies but not the other. Phenotypic divergence increased with altitude for a few biomass‐related traits, but only in A. m. striatum. These traits therefore potentially reflect A. m. striatum adaptation to altitude. Our findings imply that adaptive processes potentially differ at the scale of A. majus subspecies.