EXPERIMENTAL EVIDENCE FOR THE EVOLUTION OF INDIRECT GENETIC EFFECTS: CHANGES IN THE INTERACTION EFFECT COEFFICIENT,PSI (Ψ), DUE TO SEXUAL SELECTION

Indirect genetics effects (IGEs)—when the genotype of one individual affects the phenotypic expression of a trait in another—may alter evolutionary trajectories beyond that predicted by standard quantitative genetic theory as a consequence of genotypic evolution of the social environment. For IGEs to occur, the trait of interest must respond to one or more indicator traits in interacting conspecifics. In quantitative genetic models of IGEs, these responses (reaction norms) are termed interaction effect coefficients and are represented by the parameter psi (Ψ). The extent to which Ψ exhibits genetic variation within a population, and may therefore itself evolve, is unknown. Using an experimental evolution approach, we provide evidence for a genetic basis to the phenotypic response caused by IGEs on sexual display traits in Drosophila serrata. We show that evolution of the response is affected by sexual but not natural selection when flies adapt to a novel environment. Our results indicate a further mechanism by which IGEs can alter evolutionary trajectories—the evolution of interaction effects themselves.     

Rapid genome‐wide evolution in Brassica rapa populations following drought revealed by sequencing of ancestral and descendant gene pools

There is increasing evidence that evolution can occur rapidly in response to selection. Recent advances in sequencing suggest the possibility of documenting genetic changes as they occur in populations, thus uncovering the genetic basis of evolution, particularly if samples are available from both before and after selection. Here, we had a unique opportunity to directly assess genetic changes in natural populations following an evolutionary response to a fluctuation in climate. We analysed genome‐wide differences between ancestors and descendants of natural populations of Brassica rapa plants from two locations that rapidly evolved changes in multiple phenotypic traits, including flowering time, following a multiyear late‐season drought in California. These ancestor‐descendant comparisons revealed evolutionary shifts in allele frequencies in many genes. Some genes showing evolutionary shifts have functions related to drought stress and flowering time, consistent with an adaptive response to selection. Loci differentiated between ancestors and descendants (F_ST outliers) were generally different from those showing signatures of selection based on site frequency spectrum analysis (Tajima's D), indicating that the loci that evolved in response to the recent drought and those under historical selection were generally distinct. Very few genes showed similar evolutionary responses between two geographically distinct populations, suggesting independent genetic trajectories of evolution yielding parallel phenotypic changes. The results show that selection can result in rapid genome‐wide evolutionary shifts in allele frequencies in natural populations, and highlight the usefulness of combining resurrection experiments in natural populations with genomics for studying the genetic basis of adaptive evolution.     

Adaptive differentiation in floral traits in the presence of high gene flow in scarlet gilia (Ipomopsis aggregata)

Plant–pollinator interactions are thought to be major drivers of floral trait diversity. However, the relative importance of divergent pollinator‐mediated selection vs. neutral processes in floral character evolution has rarely been explored. We tested for adaptive floral trait evolution by comparing differentiation at neutral genetic loci to differentiation at quantitative floral traits in a putative Ipomopsis aggregata hybrid zone. Typical I. aggregata subsp. candida displays slender white tubular flowers that are typical of flowers pollinated by hawkmoths, and subsp. collina displays robust red tubular flowers typical of flowers pollinated by hummingbirds; yet, hybrid flower morphs are abundant across the East Slope of the Colorado Rockies. We estimated genetic differentiation (F_ST) for nuclear and chloroplast microsatellite loci and used a half‐sib design to calculate quantitative trait divergence (Q_ST) from collection sites across the morphological hybrid zone. We found little evidence for population structure and estimated mean F_ST to be 0.032. Q_ST values for several floral traits including corolla tube length and width, colour, and nectar volume were large and significantly greater than mean F_ST. We performed multivariate comparisons of neutral loci to genetic correlations within and between populations and found a strong signal for divergent selection, suggesting that specific combinations of floral display and reward traits may be the targets of selection. Our results show little support for historical subspecies categories, yet floral traits are more diverged than expected due to drift alone. Non‐neutral divergence for multivariate quantitative traits suggests that selection by pollinators is maintaining a correlation between display and reward traits.     

Condition dependence and the maintenance of genetic variance in a sexually dimorphic black scavenger fly

The maintenance of genetic variation in traits under strong sexual selection is a longstanding problem in evolutionary biology. The genic capture model proposes that this problem can be explained by the evolution of condition dependence in exaggerated male traits. We tested the predictions that condition dependence should be more pronounced in male sexual traits and that genetic variance in expression of these traits should increase under stress as among‐genotype variation in overall condition is exposed. Genetic variance in female and nonsexual traits should, by contrast, be similar across environments as a result of stabilizing selection on trait expression. The relationship between the degree of sexual dimorphism, condition dependence and additive genetic variance (V_a) was assessed for two morphological traits (body size and relative fore femur width) affecting male mating success in the black scavenger fly Sepsis punctum (Diptera: Sepsidae) and for development time (a nonsexual trait often correlated with body size). We compared trait expression between the sexes for two cross‐continental populations that differ in degree of sexual dimorphism (Ottawa and Zurich). Condition dependence was indeed most pronounced in males of the strongly dimorphic Zurich population (males larger), and V_a was similar for males and females unless the trait was strongly sex specific and condition dependent. Contrary to prediction, however, V_a primarily increased under food limitation in both sexes, and genetic variance in fore femur width was low to nil, perhaps depleted by putatively strong sexual selection. Solely for body size of Zurich males, V_a increased more in males than females at limited food, in accordance with the predictions of the genic capture model. Overall therefore, quantitative genetic evidence in support of the model was inconsistent and weak at best.     

Change in sexual signalling traits outruns morphological divergence across an ecological gradient in the post‐glacial radiation of the songbird genus Junco

The relative roles of natural and sexual selection in promoting evolutionary lineage divergence remains controversial and difficult to assess in natural systems. Local adaptation through natural selection is known to play a central role in promoting evolutionary divergence, yet secondary sexual traits can vary widely among species in recent radiations, suggesting that sexual selection may also be important in the early stages of speciation. Here, we compare rates of divergence in ecologically relevant traits (morphology) and sexually selected signalling traits (coloration) relative to neutral structure in genome‐wide molecular markers and examine patterns of variation in sexual dichromatism to explore the roles of natural and sexual selection in the diversification of the songbird genus Junco (Aves: Passerellidae). Juncos include divergent lineages in Central America and several dark‐eyed junco (J. hyemalis) lineages that diversified recently as the group recolonized North America following the last glacial maximum (ca. 18,000 years ago). We found an accelerated rate of divergence in sexually selected characters relative to ecologically relevant traits. Moreover, sexual dichromatism measurements suggested a positive relationship between the degree of colour divergence and the strength of sexual selection when controlling for neutral genetic distance. We also found a positive correlation between dichromatism and latitude, which coincides with the geographic axis of decreasing lineage age in juncos but also with a steep ecological gradient. Finally, we found significant associations between genome‐wide variants linked to functional genes and proxies of both sexual and natural selection. These results suggest that the joint effects of sexual and ecological selection have played a prominent role in the junco radiation.     

Polymorphisms in a desaturase 2 ortholog associate with cuticular hydrocarbon and male mating success variation in a natural population of Drosophila serrata

Elucidating the nature of genetic variation underlying both sexually selected traits and the fitness components of sexual selection is essential to understanding the broader consequences of sexual selection as an evolutionary process. To date, there have been relatively few attempts to connect the genetic variance in sexually selected traits with segregating DNA sequence polymorphisms. We set out to address this in a well‐characterized sexual selection system – the cuticular hydrocarbons (CHCs) of Drosophila serrata – using an indirect association study design that allowed simultaneous estimation of the genetic variance in CHCs, sexual fitness and single nucleotide polymorphism (SNP) effects in an outbred population. We cloned and sequenced an ortholog of the D. melanogaster desaturase 2 gene, previously shown to affect CHC biosynthesis in D. melanogaster, and associated 36 SNPs with minor allele frequencies > 0.02 with variance in CHCs and sexual fitness. Three SNPs had significant multivariate associations with CHC phenotype (q‐value < 0.05). At these loci, minor alleles had multivariate effects on CHCs that were weakly associated with the multivariate direction of sexual selection operating on these traits. Two of these SNPs had pleiotropic associations with male mating success, suggesting these variants may underlie responses to sexual selection due to this locus. There were 15 significant male mating success associations (q‐value < 0.1), and interestingly, we detected a nonrandom pattern in the relationship between allele frequency and direction of effect on male mating success. The minor‐frequency allele usually reduced male mating success, suggesting a positive association between male mating success and total fitness at this locus.     

Selection on life‐history traits and genetic population divergence in rotifers

A combination of founder effects and local adaptation – the Monopolization hypothesis – has been proposed to reconcile the strong population differentiation of zooplankton dwelling in ponds and lakes and their high dispersal abilities. The role genetic drift plays in genetic differentiation of zooplankton is well documented, but the impact of natural selection has received less attention. Here, we compare differentiation in neutral genetic markers (F_ST) and in quantitative traits (Q_ST) in six natural populations of the rotifer Brachionus plicatilis to assess the importance of natural selection in explaining genetic differentiation of life‐history traits. Five life‐history traits were measured in four temperature × salinity combinations in common‐garden experiments. Population differentiation for neutral genetic markers – 11 microsatellite loci – was very high (F_ST = 0.482). Differentiation in life‐history traits was higher in traits related to sexual reproduction than in those related to asexual reproduction. Q_ST values for diapausing egg production (a trait related to sexual reproduction) were higher than their corresponding F_ST in some pairs of populations. Our results indicate the importance of divergent natural selection in these populations and suggest local adaptation to the unpredictability of B. plicatilis habitats.     

Evolutionary changes in plant tolerance against herbivory through a resurrection experiment

Both theoretical and empirical works have highlighted the difference in the evolutionary implications of host resistance and tolerance against their enemies. However, it has been difficult to show evolutionary changes in host defences in natural populations; thus, evaluating theoretical predictions of simultaneous evolution of defences remains a challenge. We studied the evolutionary changes in traits related to resistance and tolerance against herbivory in a natural plant population using seeds from two collections made in a period of 20 years. In a common garden experiment, we compared defensive traits of ancestral (1987) and descendant (2007) subpopulations of the annual plant Datura stramonium that shows genetic variation for tolerance and to which the specialist herbivore Lema daturaphila is locally adapted. We also examined the effects of different plant genotypes on the herbivore for testing the plant genetic variation in resistance. Based on the response to the contemporary herbivore populations, results revealed a nonsignificant response in plant resistance traits (herbivore consumption, foliar trichomes and tropane alkaloids), but a significant one in tolerance. The survival of herbivores in laboratory experiments depended on the plant genotype, which suggests genetic variation in plant resistance. Although we cannot identify the selective agent for the change nor exclude genetic drift, the results are consistent with the expectation that when resistance fails to control herbivory, tolerance should play a more important role in the evolution of the interaction.     

An empirical test for a zone of canalization in thermal reaction norms

Theoretical models on the evolution of phenotypic plasticity predict a zone of canalization where reaction norms cross, and genetic variation is minimized in the environment a population most frequently encounter. Empirical tests of this prediction are largely missing, in particular for life‐history traits. We addressed this prediction by quantifying thermal reaction norms of three life‐history traits (somatic growth rate, age and size at maturation) of a Norwegian population of Daphnia magna and testing for the occurrence of an intermediate temperature (T_m) at which genetic variance in the traits is minimized. Size at maturation changed relatively little with temperature compared to the other traits, and there was no genetic variance in the shape of the reaction norm. Consequently, age at maturation and somatic growth rate were strongly negatively correlated. Both traits showed a strong genotype–environment interaction, and the estimated T_m was 14 °C for both age at maturation and growth rate. This value of T_m corresponds well with mean summer temperatures experienced by the population and suggests that the population has evolved under stabilizing selection in temperatures that fluctuate around this mean temperature. These results suggest local adaptation to temperature in the studied population and allow predicting evolutionary trajectories of thermal reaction norms under changing thermal regimes.     

Parasite‐mediated selection in a natural metapopulation of Daphnia magna

Parasite‐mediated selection varying across time and space in metapopulations is expected to result in host local adaptation and the maintenance of genetic diversity in disease‐related traits. However, nonadaptive processes like migration and extinction‐(re)colonization dynamics might interfere with adaptive evolution. Understanding how adaptive and nonadaptive processes interact to shape genetic variability in life‐history and disease‐related traits can provide important insights into their evolution in subdivided populations. Here we investigate signatures of spatially fluctuating, parasite‐mediated selection in a natural metapopulation of Daphnia magna. Host genotypes from infected and uninfected populations were genotyped at microsatellite markers, and phenotyped for life‐history and disease traits in common garden experiments. Combining phenotypic and genotypic data a Q_ST–F_ST‐like analysis was conducted to test for signatures of parasite mediated selection. We observed high variation within and among populations for phenotypic traits, but neither an indication of host local adaptation nor a cost of resistance. Infected populations have a higher gene diversity (Hs) than uninfected populations and Hs is strongly positively correlated with fitness. These results suggest a strong parasite effect on reducing population level inbreeding. We discuss how stochastic processes related to frequent extinction‐(re)colonization dynamics as well as host and parasite migration impede the evolution of resistance in the infected populations. We suggest that the genetic and phenotypic patterns of variation are a product of dynamic changes in the host gene pool caused by the interaction of colonization bottlenecks, inbreeding, immigration, hybrid vigor, rare host genotype advantage and parasitism. Our study highlights the effect of the parasite in ameliorating the negative fitness consequences caused by the high drift load in this metapopulation.     

The evolutionary consequences of habitat fragmentation: Body morphology and coloration differentiation among brook trout populations of varying size

A reduction in population size due to habitat fragmentation can alter the relative roles of different evolutionary mechanisms in phenotypic trait differentiation. While deterministic (selection) and stochastic (genetic drift) mechanisms are expected to affect trait evolution, genetic drift may be more important than selection in small populations. We examined relationships between mature adult traits and ecological (abiotic and biotic) variables among 14 populations of brook trout. These naturally fragmented populations have shared ancestry but currently exhibit considerable variability in habitat characteristics and population size (49 < N_c < 10,032; 3 < N_b < 567). Body size, shape, and coloration differed among populations, with a tendency for more variation among small populations in both trait means and CV when compared to large populations. Phenotypic differences were more frequently and directly linked to habitat variation or operational sex ratio than to population size, suggesting that selection may overcome genetic drift at small population size. Phenotype–environment associations were also stronger in females than males, suggesting that natural selection due to abiotic conditions may act more strongly on females than males. Our results suggest that natural and sexual‐selective pressures on phenotypic traits change during the process of habitat fragmentation, and that these changes are largely contingent upon existing habitat conditions within isolated fragments. Our study provides an improved understanding of the ecological and evolutionary consequences of habitat fragmentation and lends insight into the ability of some small populations to respond to selection and environmental change.     

In situ marker‐based assessment of leaf trait evolutionary potential in a marginal European beech population

Evolutionary processes are expected to be crucial for the adaptation of natural populations to environmental changes. In particular, the capacity of rear edge populations to evolve in response to the species limiting conditions remains a major issue that requires to address their evolutionary potential. In situ quantitative genetic studies based on molecular markers offer the possibility to estimate evolutionary potentials manipulating neither the environment nor the individuals on which phenotypes are measured. The goal of this study was to estimate heritability and genetic correlations of a suite of leaf functional traits involved in climate adaptation for a natural population of the tree Fagus sylvatica, growing at the rear edge of the species range. Using two marker‐based quantitative genetics approaches, we obtained consistent and significant estimates of heritability for leaf phenological (phenology of leaf flush), morphological (mass, area, ratio mass/area) and physiological (δ¹³C, nitrogen content) traits. Moreover, we found only one significant positive genetic correlation between leaf area and leaf mass, which likely reflected mechanical constraints. We conclude first that the studied population has considerable genetic diversity for important ecophysiological traits regarding drought adaptation and, second, that genetic correlations are not likely to impose strong genetic constraints to future population evolution. Our results bring important insights into the question of the capacity of rear edge populations to evolve.     

A multivariate test of evolutionary constraints for thermal tolerance in Drosophila melanogaster

Exposure to extreme temperatures is increasingly likely to impose strong selection on many organisms in their natural environments. The ability of organisms to adapt to such selective pressures will be determined by patterns of genetic variation and covariation. Despite increasing interest in thermal adaptation, few studies have examined the extent to which the genetic covariance between traits might constrain thermal responses. Furthermore, it remains unknown whether sex‐specific genetic architectures will constrain responses to climatic selection. We used a paternal half‐sibling breeding design to examine whether sex‐specific genetic architectures and genetic covariances between traits might constrain evolutionary responses to warming climates in a population of Drosophila melanogaster. Our results suggest that the sexes share a common genetic underpinning for heat tolerance as indicated by a strong positive inter‐sexual genetic correlation. Further, we found no evidence in either of the sexes that genetic trade‐offs between heat tolerance and fitness will constrain responses to thermal selection. Our results suggest that neither trade‐offs, nor sex‐specific genetics, will significantly constrain an evolutionary response to climatic warming, at least in this population of D. melanogaster.     

The role of immigration and local adaptation on fine‐scale genotypic and phenotypic population divergence in a less mobile passerine

Dispersal and local patterns of adaptation play a major role on the ecological and evolutionary trajectory of natural populations. In this study, we employ a combination of genetic (25 microsatellite markers) and field‐based information (seven study years) to analyse the impact of immigration and local patterns of adaptation in two nearby (< 7 km) blue tit (Cyanistes caeruleus) populations. We used genetic assignment analyses to identify immigrant individuals and found that dispersal rate is female‐biased (72%). Data on lifetime reproductive success indicated that immigrant females produced fewer local recruits than their philopatric counterparts whereas immigrant males recruited more offspring than those that remained in their natal location. In spite of the considerably higher immigration rates of females, our results indicate that, in absolute terms, their demographic and genetic impact in the receiving populations is lower than that in immigrant males. Immigrants often brought novel alleles into the studied populations and a high proportion of them were transmitted to their recruits, indicating that the genetic impact of immigrants is not ephemeral. Although only a few kilometres apart, the two study populations were genetically differentiated and showed strong divergence in different phenotypic and life‐history traits. An almost absent inter‐population dispersal, together with the fact that both populations receive immigrants from different source populations, is probably the main cause of the observed pattern of genetic differentiation. However, phenotypic differentiation (P_ST) for all the studied traits greatly exceeded neutral genetic differentiation (F_ST), indicating that divergent natural selection is the prevailing factor determining the evolutionary trajectory of these populations. Our study highlights the importance of integrating individual‐ and population‐based approaches to obtain a comprehensive view about the role of dispersal and natural selection on structuring the genotypic and phenotypic characteristics of natural populations.     

Genome scans for the contemporary response to selection in quantitative traits

Genome scans have been an important approach for discovering historical signatures of selection in both model and nonmodel species. An exciting new experimental design for genome scans is to measure the change in allele frequency before and after contemporary selection within a generation, from a single population. The most widely‐used methods, however, have two major limitations: they are based on testing one locus at a time, and they only have power to uncover loci that have evolved under relatively strong selection. On the other hand, complex quantitative traits are common in nature and are caused by several loci of small effect. Selection on a quantitative trait at the phenotypic level is predicted to be accompanied by subtle allele frequency changes in many loci that covary (a polygenic soft sweep), rather than a large, single‐effect allele (a selective sweep). In this issue of Molecular Ecology, Bourret et al. (2014) measure the contemporary response to natural selection across the genome in multiple cohorts of Atlantic salmon during their first year at sea. They introduce a multilocus framework based on groups of markers that covary in their genotypic distribution. While the traditional, single‐locus approach did not find evidence for repeated patterns of selection, the multivariate approach found that a group of covarying SNPs was selected for in different cohorts at one site. Their multilocus framework has potential to be a more fruitful approach for uncovering the genomic basis of adaptation in quantitative traits, although caution should be applied as the framework has yet to be validated with simulated data.     

Evolutionary stasis despite selection on a heritable trait in an invasive zooplankton

Invasive species are one of the greatest threats to ecosystems, and there is evidence that evolution plays an important role in the success or failure of invasions. Yet, few studies have measured natural selection and evolutionary responses to selection in invasive species, particularly invasive animals. We quantified the strength of natural selection on the defensive morphology (distal spine) of an invasive zooplankton, Bythotrephes longimanus, in Lake Michigan across multiple months during three growing seasons. We used multiple lines of evidence, including historic and contemporary wild‐captured individuals and palaeoecology of retrieved spines, to assess phenotypic change in distal spine length since invasion. We found evidence of temporally variable selection, with selection for decreased distal spine length early in the growing season and selection for increased distal spine length later in the season. This trend in natural selection is consistent with seasonal changes in the relative strength of non‐gape‐limited and gape‐limited fish predation. Yet, despite net selection for increased distal spine length and a known genetic basis for distal spine length, we observed little evidence of an evolutionary response to selection. Multiple factors likely limit an evolutionary response to selection, including genetic correlations, trade‐offs between components of fitness, and phenotypic plasticity.     

Sexual selection and population divergence II. Divergence in different sexual traits and signal modalities in field crickets (Teleogryllus oceanicus)

Sexual selection can target many different types of traits. However, the relative influence of different sexually selected traits during evolutionary divergence is poorly understood. We used the field cricket Teleogryllus oceanicus to quantify and compare how five traits from each of three sexual signal modalities and components diverge among allopatric populations: male advertisement song, cuticular hydrocarbon (CHC) profiles and forewing morphology. Population divergence was unexpectedly consistent: we estimated the among‐population (genetic) variance‐covariance matrix, D , for all 15 traits, and D_max explained nearly two‐thirds of its variation. CHC and wing traits were most tightly integrated, whereas song varied more independently. We modeled the dependence of among‐population trait divergence on genetic distance estimated from neutral markers to test for signatures of selection versus neutral divergence. For all three sexual trait types, phenotypic variation among populations was largely explained by a neutral model of divergence. Our findings illustrate how phenotypic integration across different types of sexual traits might impose constraints on the evolution of mating isolation and divergence via sexual selection.     

An evolutionary limit to male mating success

The well-known phenotypic diversity of male sexual displays, and the high levels of genetic variation reported for individual display traits have generated the expectation that male display traits, and consequently male mating success, are highly evolvable. It has not been shown however that selection for male mating success, exerted by female preferences in an unmanipulated population, results in evolutionary change. Here, we tested the expectation that male mating success is highly evolvable in Drosophila bunnanda using an experimental evolution approach. Female D. bunnanda exhibit a strong, consistent preference for a specific combination of male cuticular hydrocarbons (CHCs). We used female preference to select for male mating success by propagating replicate populations from either attractive or unattractive males over 10 generations. Neither the combination of CHCs under sexual selection (the sexual signal) nor male mating success itself evolved. The lack of a response to selection was consistent with previous quantitative genetic experiments in D. bunnanda that demonstrated the virtual absence of genetic variance in the combination of CHCs under sexual selection. Persistent directional selection, such as applied by female mate choice, may erode genetic variance, resulting in multitrait evolutionary limits.     

Population genomic,olfactory, dietary,and gut microbiota analyses demonstrate the unique evolutionary trajectory of feral pigs

Domestication is an intriguing evolutionary process. Many domestic populations are subjected to strong human-mediated selection, and when some individuals return to the wild, they are again subjected to selective forces associated with new environments. Generally, these feral populations evolve into something different from their wild predecessors and their members typically possess a combination of both wild and human selected traits. Feralisation can manifest in different forms on a spectrum from a wild to a domestic phenotype. This depends on how the rewilded domesticated populations can readapt to natural environments based on how much potential and flexibility the ancestral genome retains after its domestication signature. Whether feralisation leads to the evolution of new traits that do not exist in the wild or to convergence with wild forms, however, remains unclear. To address this question, we performed population genomic, olfactory, dietary, and gut microbiota analyses on different populations of Sus scrofa (wild boar, hybrid, feral and several domestic pig breeds). Porcine single nucleotide polymorphisms (SNPs) analysis shows that the feral population represents a cluster distinctly separate from all others. Its members display signatures of past artificial selection, as demonstrated by values of F_ST in specific regions of the genome and bottleneck signature, such as the number and length of runs of homozygosity. Generalised F_ST values, reacquired olfactory abilities, diet, and gut microbiota variation show current responses to natural selection. Our results suggest that feral pigs are an independent evolutionary unit which can persist so long as levels of human intervention remain unchanged.     

Multiple aspects of the selfing syndrome of the morning glory Ipomoea lacunosa evolved in response to selection: A Qst‐Fst comparison

The frequent transition from outcrossing to selfing in flowering plants is often accompanied by changes in multiple aspects of floral morphology, termed the “selfing syndrome.” While the repeated evolution of these changes suggests a role for natural selection, genetic drift may also be responsible. To determine whether selection or drift shaped different aspects of the pollination syndrome and mating system in the highly selfing morning glory Ipomoea lacunosa, we performed multivariate and univariate Qst‐Fst comparisons using a wide sample of populations of I. lacunosa and its mixed‐mating sister species Ipomoea cordatotriloba. The two species differ in early growth, floral display, inflorescence traits, corolla size, nectar, and pollen number. Our analyses support a role for natural selection driving trait divergence, specifically in corolla size and nectar traits, but not in early growth, display size, inflorescence length, or pollen traits. We also find evidence of selection for reduced herkogamy in I. lacunosa, consistent with selection driving both the transition in mating system and the correlated floral changes. Our research demonstrates that while some aspects of the selfing syndrome evolved in response to selection, others likely evolved due to drift or correlated selection, and the balance between these forces may vary across selfing species.