TETRODOTOXIN RESISTANCE IN GARTER SNAKES: AN EVOLUTIONARY RESPONSE OF PREDATORS TO DANGEROUS PREY

The use of the “arms race” analogy as a conceptualization of evolutionary predator-prey interactions has been criticized because of the lack of evidence that predators can and do adapt to increased antipredator ability of prey. We present evidence that the garter snake Thamnophis sirtalis has evolved resistance to tetrodotoxin (TTX) in response to the toxicity of the newt Taricha granulosa on which the snake feeds. A bioassay (locomotor performance before and after injection of TTX) was used to obtain repeated measures of resistance for individual snakes. We studied interpopulation and interspecific variation by comparing resistance in Thamnophis sirtalis from populations occurring sympatrically and allopatrically with Taricha granulosa, and in Thamnophis ordinoides (which does not feed on the newt) occurring sympatrically with Taricha granulosa. We also examined intrapopulation variation in TTX resistance using snakes from a population known to feed on Taricha granulosa. Resistance differed significantly among individuals and litters; repeatability and heritability estimates of the assay were significantly different from zero, demonstrating the potential for response to selection. The population of Thamnophis sirtalis that occurs with Taricha granulosa exhibited levels of resistance much greater than either of the other groups. These results suggest that the predator-prey arms race analogy may be applicable to this system.     

Patterns of genetic differentiation in Thamnophis and Taricha from the Pacific Northwest

Aim The co‐evolutionary interaction between the common garter snake, Thamnophis sirtalis, and the rough‐skinned newt, Taricha granulosa, takes place throughout much of the Pacific Northwest (North America). The biogeography of the Pacific Northwest has been heavily influenced by the last Pleistocene glaciation, which reached a maximum as late as 14,000 yr bp . We researched: (1) what type of population structure is present for garter snakes and newts, (2) whether the population structure of these species is consistent with a Pleistocene glaciation hypothesis, and (3) how population structure and migration possibly affect co‐evolution between these species. Location The Pacific Northwest of North America, specifically northern California, Oregon and Washington in the USA. Methods We sampled approximately 20 populations for each species from three different transects. Using microsatellite markers and tissue samples from both species, we quantified the population structure for both species. Individual‐based assignment tests were used to estimate contemporary migration rates. Results Both Th. sirtalis and Ta. granulosa exhibited little genetic differentiation among our study sites, even among those separated by large distances. Significant population structure was detected on multiple geographic scales. Differences in population structure were observed among transects and between garter snake and newt transects. Contemporary migration rate estimates indicate high levels of genetic exchange between populations. Main conclusions Prior to this study, little was known about the fine‐scale population structure of either species in this region. Patterns of population structure for garter snakes and newts reflect a shared biogeographical history affected by the Pleistocene glaciation in the Pacific Northwest. Both species apparently migrate frequently between populations, thus potentially retarding the process of adaptive co‐evolution. We find that populations from a northern coastal transect (Washington) are most likely to be locally adapted.     

EVOLUTIONARY RESPONSE OF PREDATORS TO DANGEROUS PREY: PREADAPTATION AND THE EVOLUTION OF TETRODOTOXIN RESISTANCE IN GARTER SNAKES

Coevolutionary interactions typically involve only a few specialized taxa. The factors that cause some taxa and not others to respond evolutionarily to selection by another species are poorly understood. Preadaptation may render some species predisposed for evolutionary response to new pressures, whereas a lack of genetic variation may limit the evolutionary potential of other taxa. We evaluate these factors in the predator-prey interaction between toxic newts (Taricha granulosa) and their resistant garter snake predators (Thamnophis sirtalis). Using a bioassay of resistance to tetrodotoxin (TTX), the primary toxin in the prey, we examined phenotypic evolution in the genus Thamnophis. Reconstruction of ancestral character states suggests that the entire genus Thamnophis, and possibly natricine snakes in general, has slightly elevated TTX resistance compared to other lineages of snakes. While this suggests that T. sirtalis is indeed predisposed to evolving TTX resistance, it also indicates that the potential exists in sympatric congeners not expressing elevated levels of TTX resistance. We also detected significant family level variation for TTX resistance in a species of Thamnophis that does not exhibit elaborated levels of the trait. This finding suggests that evolutionary response in other taxa is not limited by genetic variability. In this predator-prey system, species and population differences in resistance appear to be largely determined by variation in the selective environment rather than preadaptation or constraint.     

Genetic architecture of a feeding adaptation: garter snake (Thamnophis) resistance to tetrodotoxin bearing prey

Detailing the genetic basis of adaptive variation in natural populations is a first step towards understanding the process of adaptive evolution, yet few ecologically relevant traits have been characterized at the genetic level in wild populations. Traits that mediate coevolutionary interactions between species are ideal for studying adaptation because of the intensity of selection and the well-characterized ecological context. We have previously described the ecological context, evolutionary history and partial genetic basis of tetrodotoxin (TTX) resistance in garter snakes (Thamnophis). Derived mutations in a voltage-gated sodium channel gene (Na_v1.4) in three garter snake species are associated with resistance to TTX, the lethal neurotoxin found in their newt prey (Taricha). Here we evaluate the contribution of Na_v1.4 alleles to TTX resistance in two of those species from central coastal California. We measured the phenotypes (TTX resistance) and genotypes (Na_v1.4 and microsatellites) in a local sample of Thamnophis atratus and Thamnophis sirtalis. Allelic variation in Na_v1.4 explains 23 per cent of the variation in TTX resistance in T. atratus while variation in a haphazard sample of the genome (neutral microsatellite markers) shows no association with the phenotype. Similarly, allelic variation in Na_v1.4 correlates almost perfectly with TTX resistance in T. sirtalis, but neutral variation does not. These strong correlations suggest that Na_v1.4 is a major effect locus. The simple genetic architecture of TTX resistance in garter snakes may significantly impact the dynamics of phenotypic coevolution. Fixation of a few alleles of major effect in some garter snake populations may have led to the evolution of extreme phenotypes and an ‘escape’ from the arms race with newts.     

Links between prey assemblages and poison frog toxins: A landscape ecology approach to assess how biotic interactions affect species phenotypes

Ecological studies of species pairs showed that biotic interactions promote phenotypic change and eco‐evolutionary feedbacks. However, it is unclear how phenotypes respond to synergistic interactions with multiple taxa. We investigate whether interactions with multiple prey species explain spatially structured variation in the skin toxins of the neotropical poison frog Oophaga pumilio. Specifically, we assess how dissimilarity (i.e., beta diversity) of alkaloid‐bearing arthropod prey assemblages (68 ant species) and evolutionary divergence between frog populations (from a neutral genetic marker) contribute to frog poison dissimilarity (toxin profiles composed of 230 different lipophilic alkaloids sampled from 934 frogs at 46 sites). We find that models that incorporate spatial turnover in the composition of ant assemblages explain part of the frog alkaloid variation, and we infer unique alkaloid combinations across the range of O. pumilio. Moreover, we find that alkaloid variation increases weakly with the evolutionary divergence between frog populations. Our results pose two hypotheses: First, the distribution of only a few prey species may explain most of the geographic variation in poison frog alkaloids; second, different codistributed prey species may be redundant alkaloid sources. The analytical framework proposed here can be extended to other multitrophic systems, coevolutionary mosaics, microbial assemblages, and ecosystem services.     

Geographic variation in advertisement calls of a Microhylid frog – testing the role of drift and ecology

Acoustic signals for mating are important traits that could drive population differentiation and speciation. Ecology may play a role in acoustic divergence through direct selection (e.g., local adaptation to abiotic environment), constraint of correlated traits (e.g., acoustic traits linked to another trait under selection), and/or interspecific competition (e.g., character displacement). However, genetic drift alone can also drive acoustic divergence. It is not always easy to differentiate the role of ecology versus drift in acoustic divergence. In this study, we tested the role of ecology and drift in shaping geographic variation in the advertisement calls of Microhyla fissipes. We examined three predictions based on ecological processes: (1) the correlation between temperature and call properties across M. fissipes populations; (2) the correlation between call properties and body size across M. fissipes populations; and (3) reproductive character displacement (RCD) in call properties between M. fissipes populations that are sympatric with and allopatric to a congener M. heymonsi. To test genetic drift, we examined correlations among call divergence, geographic distance, and genetic distance across M. fissipes populations. We recorded the advertisement calls from 11 populations of M. fissipes in Taiwan, five of which are sympatrically distributed with M. heymonsi. We found geographic variation in both temporal and spectral properties of the advertisement calls of M. fissipes. However, the call properties were not correlated with local temperature or the callers' body size. Furthermore, we did not detect RCD. By contrast, call divergence, geographic distance, and genetic distance between M. fissipes populations were all positively correlated. The comparisons between phenotypic Q_st (P_st) and F_st values did not show significant differences, suggesting a role of drift. We concluded that genetic drift, rather than ecological processes, is the more likely driver for the geographic variation in the advertisement calls of M. fissipes.     

The relative contribution of drift and selection to phenotypic divergence: A test case using the horseshoe bats Rhinolophus simulator and Rhinolophus swinnyi

Natural selection and drift can act on populations individually, simultaneously or in tandem and our understanding of phenotypic divergence depends on our ability to recognize the contribution of each. According to the quantitative theory of evolution, if an organism has diversified through neutral evolutionary processes (mutation and drift), variation of phenotypic characteristics between different geographic localities (B) should be directly proportional to the variation within localities (W), that is, B ∝ W. Significant deviations from this null model imply that non‐neutral forces such as natural selection are acting on a phenotype. We investigated the relative contributions of drift and selection to intraspecific diversity using southern African horseshoe bats as a test case. We characterized phenotypic diversity across the distributional range of Rhinolophus simulator (n = 101) and Rhinolophus swinnyi (n = 125) using several traits associated with flight and echolocation. Our results suggest that geographic variation in both species was predominantly caused by disruptive natural selection (B was not directly proportional to W). Evidence for correlated selection (co‐selection) among traits further confirmed that our results were not compatible with drift. Selection rather than drift is likely the predominant evolutionary process shaping intraspecific variation in traits that strongly impact fitness.     

Genomic variation across two barn swallow hybrid zones reveals traits associated with divergence in sympatry and allopatry

Hybrid zones are geographic regions where isolating barriers between divergent populations are challenged by admixture. Identifying factors that facilitate or inhibit hybridization in sympatry can illuminate the processes that maintain those reproductive barriers. We analysed patterns of hybridization and phenotypic variation across two newly discovered hybrid zones between three subspecies of barn swallow (Hirundo rustica). These subspecies differ in ventral coloration and wing length, traits that are targets of sexual and natural selection, respectively, and are associated with genome‐wide differentiation in allopatry. We tested the hypothesis that the degree of divergence in these traits is associated with the extent of hybridization in secondary contact. We applied measures of population structure based on >23,000 SNPs to confirm that named subspecies correspond to distinct genomic clusters, and assessed coincidence between geographic clines for ancestry and phenotype. Although gene flow was ongoing across both hybrid zones and pairwise F_ST between subspecies was extremely low, we found striking differences in the extent of hybridization. In the more phenotypically differentiated subspecies pair, clines for ancestry, wing length and ventral coloration were steep and coincident, suggestive of strong isolation and, potentially, selection associated with phenotype. In the less phenotypically differentiated pair, gene flow and phenotypic variation occurred over a wide geographic span, indicative of weaker isolation. Traits associated with genome‐wide differentiation in allopatry may thus also contribute to isolation in sympatry. We discuss potentially important additional roles for evolutionary history and ecology in shaping variation in the extent hybridization between closely related pairs of subspecies.     

Microgeographic morphological variation across larval wood frog populations associated with environment despite gene flow

Gene flow has historically been thought to constrain local adaptation; yet, recent research suggests that populations can diverge despite exchanging genes. Here I use a common garden experiment to assess the combined effects of gene flow and natural selection on morphological variation of 16 wood frog (Rana sylvatica) populations, a species known to experience divergent selection pressures in open‐ and closed‐canopy ponds across relatively small geographic scales. Wood frog tadpoles from different ponds showed significant morphological variation associated with canopy type with a trade‐off between tail length and body depth consistent with previous research. In contrast, neutral genetic differentiation of nine microsatellite loci as measured by Jost's D was not associated with canopy type, indicating no pattern of isolation by environment. Genetic structure analyses indicated some substructure across the 16 ponds (K = 4); however, three out of four assigned clusters included both open‐ and closed‐canopy ponds. Together, these results suggest that morphological divergence among these wood frog populations is occurring despite gene flow and that selection within these environments is strong. Furthermore, morphological variation among ponds differed across two sampling periods during larval development, demonstrating the importance of evaluating phenotypic divergence over multiple time periods and at a time relevant to the processes being studied.     

The evolutionary response of predators to dangerous prey: hotspots and coldspots in the geographic mosaic of coevolution between garter snakes and newts

The "geographic mosaic" approach to understanding coevolution is predicated on the existence of variable selection across the landscape of an interaction between species. A range of ecological factors, from differences in resource availability to differences in community composition, can generate such a mosaic of selection among populations, and thereby differences in the strength of coevolution. The result is a mixture of hotspots, where reciprocal selection is strong, and coldspots, where reciprocal selection is weak or absent, throughout the ranges of species. Population subdivision further provides the opportunity for nonadaptive forces, including gene flow, drift, and metapopulation dynamics, to influence the coevolutionary interaction between species. Some predicted results of this geographic mosaic of coevolution include maladapted or mismatched phenotypes, maintenance of high levels of polymorphism, and prevention of stable equilibrium trait combinations. To evaluate the potential for the geographic mosaic to influence predator-prey coevolution, we investigated the geographic pattern of genetically determined TTX resistance in the garter snake Thamnophis sirtalis over much of the range of its ecological interaction with toxic newts of genus Taricha. We assayed TTX resistance in over 2900 garter snakes representing 333 families from 40 populations throughout western North America. Our results provide dramatic evidence that geographic structure is an important component in coevolutionary interactions between predators and prey. Resistance levels vary substantially (over three orders of magnitude) among populations and over short distances. The spatial array of variation is consistent with two areas of intense evolutionary response by predators ("hotspots") surrounded by clines of decreasing resistance. Some general predictions of the geographic mosaic process are supported, including clinal variation in phenotypes, polymorphism in some populations, and divergent outcomes of the interaction between predator and prey. Conversely, our data provide little support for one of the major predictions, mismatched values of interacting traits. Two lines of evidence suggest selection is paramount in determining population variation in resistance. First, phylogenetic information indicates that two hotspots of TTX resistance have evolved independently. Second, in the one region that TTX levels in prey have been quantified, resistance and toxicity levels match almost perfectly over a wide phenotypic and geographic range. However, these results do not preclude the role the nonadaptive forces in generating the overall geographic mosaic of TTX resistance. Much work remains to fill in the geographic pattern of variation among prey populations and, just as importantly, to explore the variation in the ecology of the interaction that occurs within populations.     

High degree of genetic differentiation in marine three‐spined sticklebacks (Gasterosteus aculeatus)

Populations of widespread marine organisms are typically characterized by a low degree of genetic differentiation in neutral genetic markers, but much less is known about differentiation in genes whose functional roles are associated with specific selection regimes. To uncover possible adaptive population divergence and heterogeneous genomic differentiation in marine three‐spined sticklebacks (Gasterosteus aculeatus), we used a candidate gene‐based genome‐scan approach to analyse variability in 138 microsatellite loci located within/close to (<6 kb) functionally important genes in samples collected from ten geographic locations. The degree of genetic differentiation in markers classified as neutral or under balancing selection—as determined with several outlier detection methods—was low (F_ST = 0.033 or 0.011, respectively), whereas average F_ST for directionally selected markers was significantly higher (F_ST = 0.097). Clustering analyses provided support for genomic and geographic heterogeneity in selection: six genetic clusters were identified based on allele frequency differences in the directionally selected loci, whereas four were identified with the neutral loci. Allelic variation in several loci exhibited significant associations with environmental variables, supporting the conjecture that temperature and salinity, but not optic conditions, are important drivers of adaptive divergence among populations. In general, these results suggest that in spite of the high degree of physical connectivity and gene flow as inferred from neutral marker genes, marine stickleback populations are strongly genetically structured in loci associated with functionally relevant genes.     

Experimental evolution of competing bean beetle species reveals long‐term reversals of short‐term evolution,but no consistent character displacement

Interspecific competition for shared resources should select for evolutionary divergence in resource use between competing species, termed character displacement. Many purported examples of character displacement exist, but few completely rule out alternative explanations. We reared genetically diverse populations of two species of bean beetles, Callosobruchus maculatus and Callosobruchus chinensis, in allopatry and sympatry on a mixture of adzuki beans and lentils, and assayed oviposition preference and other phenotypic traits after four, eight, and twelve generations of (co)evolution. C. maculatus specializes on adzuki beans; the generalist C. chinensis uses both beans. C. chinensis growing in allopatry emerged equally from both bean species. In sympatry, the two species competing strongly and coexisted via strong realized resource partitioning, with C. chinensis emerging almost exclusively from lentils and C. maculatus emerging almost exclusively from adzuki beans. However, oviposition preferences, larval survival traits, and larval development rates in both beetle species did not vary consistently between allopatric versus sympatric treatments. Rather, traits evolved in treatment‐independent fashion, with several traits exhibiting reversals in their evolutionary trajectories. For example, C. chinensis initially evolved a slower egg‐to‐adult development rate on adzuki beans in both allopatry and sympatry, then subsequently evolved back toward the faster ancestral development rate. Lack of character displacement is consistent with a previous similar experiment in bean beetles and may reflect lack of evolutionary trade‐offs in resource use. However, evolutionary reversals were unexpected and remain unexplained. Together with other empirical and theoretical work, our results illustrate the stringency of the conditions for character displacement.     

Is there more than one way to skin a newt? Convergent toxin resistance in snakes is not due to a common genetic mechanism

Convergent evolution of tetrodotoxin (TTX) resistance, at both the phenotypic and genetic levels, characterizes coevolutionary arms races between amphibians and their snake predators around the world, and reveals remarkable predictability in the process of adaptation. Here we examine the repeatability of the evolution of TTX resistance in an undescribed predator–prey relationship between TTX-bearing Eastern Newts (Notophthalmus viridescens) and Eastern Hog-nosed Snakes (Heterodon platirhinos). We found that that local newts contain levels of TTX dangerous enough to dissuade most predators, and that Eastern Hog-nosed Snakes within newt range are highly resistant to TTX. In fact, these populations of Eastern Hog-nosed Snakes are so resistant to TTX that the potential for current reciprocal selection might be limited. Unlike all other cases of TTX resistance in vertebrates, H. platirhinos lacks the adaptive amino acid substitutions in the skeletal muscle sodium channel that reduce TTX binding, suggesting that physiological resistance in Eastern Hog-nosed Snakes is conferred by an alternate genetic mechanism. Thus, phenotypic convergence in this case is not due to parallel molecular evolution, indicating that there may be more than one way for this adaptation to arise, even among closely related species.     

Genome‐wide differentiation in closely related populations: the roles of selection and geographic isolation

Population divergence in geographic isolation is due to a combination of factors. Natural and sexual selection may be important in shaping patterns of population differentiation, a pattern referred to as ‘isolation by adaptation’ (IBA). IBA can be complementary to the well‐known pattern of ‘isolation by distance’ (IBD), in which the divergence of closely related populations (via any evolutionary process) is associated with geographic isolation. The barn swallow Hirundo rustica complex comprises six closely related subspecies, where divergent sexual selection is associated with phenotypic differentiation among allopatric populations. To investigate the relative contributions of selection and geographic distance to genome‐wide differentiation, we compared genotypic and phenotypic variation from 350 barn swallows sampled across eight populations (28 pairwise comparisons) from four different subspecies. We report a draft whole‐genome sequence for H. rustica, to which we aligned a set of 9493 single nucleotide polymorphisms (SNPs). Using statistical approaches to control for spatial autocorrelation of phenotypic variables and geographic distance, we find that divergence in traits related to migratory behaviour and sexual signalling, as well as geographic distance, together explain over 70% of genome‐wide divergence among populations. Controlling for IBD, we find 42% of genomewide divergence is attributable to IBA through pairwise differences in traits related to migratory behaviour and sexual signalling alone. By (i) combining these results with prior studies of how selection shapes morphological differentiation and (ii) accounting for spatial autocorrelation, we infer that morphological adaptation plays a large role in shaping population‐level differentiation in this group of closely related populations.     

Müllerian mimicry of a quantitative trait despite contrasting levels of genomic divergence and selection

Hybrid zones, where distinct populations meet and interbreed, give insight into how differences between populations are maintained despite gene flow. Studying clines in genetic loci and adaptive traits across hybrid zones is a powerful method for understanding how selection drives differentiation within a single species, but can also be used to compare parallel divergence in different species responding to a common selective pressure. Here, we study parallel divergence of wing colouration in the butterflies Heliconius erato and H. melpomene, which are distantly related Müllerian mimics which show parallel geographic variation in both discrete variation in pigmentation, and quantitative variation in structural colour. Using geographic cline analysis, we show that clines in these traits are positioned in roughly the same geographic region for both species, which is consistent with direct selection for mimicry. However, the width of the clines varies markedly between species. This difference is explained in part by variation in the strength of selection acting on colour traits within each species, but may also be influenced by differences in the dispersal rate and total strength of selection against hybrids between the species. Genotyping‐by‐sequencing also revealed weaker population structure in H. melpomene, suggesting the hybrid zones may have evolved differently in each species, which may also contribute to the patterns of phenotypic divergence in this system. Overall, we conclude that multiple factors are needed to explain patterns of clinal variation within and between these species, although mimicry has probably played a central role.     

Drift versus selection as drivers of phenotypic divergence at small spatial scales: The case of Belgjarskógur threespine stickleback

Divergence in phenotypic traits is facilitated by a combination of natural selection, phenotypic plasticity, gene flow, and genetic drift, whereby the role of drift is expected to be particularly important in small and isolated populations. Separating the components of phenotypic divergence is notoriously difficult, particularly for multivariate phenotypes. Here, we assessed phenotypic divergence of threespine stickleback (Gasterosteus aculeatus) across 19 semi‐interconnected ponds within a small geographic region (~7.5 km²) using comparisons of multivariate phenotypic divergence (PST), neutral genetic (FST), and environmental (EST) variation. We found phenotypic divergence across the ponds in a suite of functionally relevant phenotypic traits, including feeding, defense, and swimming traits, and body shape (geometric morphometric). Comparisons of PSTs with FSTs suggest that phenotypic divergence is predominantly driven by neutral processes or stabilizing selection, whereas phenotypic divergence in defensive traits is in accordance with divergent selection. Comparisons of population pairwise PSTs with ESTs suggest that phenotypic divergence in swimming traits is correlated with prey availability, whereas there were no clear associations between phenotypic divergence and environmental difference in the other phenotypic groups. Overall, our results suggest that phenotypic divergence of these small populations at small geographic scales is largely driven by neutral processes (gene flow, drift), although environmental determinants (natural selection or phenotypic plasticity) may play a role.     

Deep sympatric mtDNA divergence in the autumnal moth (Epirrita autumnata)

Deep sympatric intraspecific divergence in mtDNA may reflect cryptic species or formerly distinct lineages in the process of remerging. Preliminary results from DNA barcoding of Scandinavian butterflies and moths showed high intraspecific sequence variation in the autumnal moth, Epirrita autumnata. In this study, specimens from different localities in Norway and some samples from Finland and Scotland, with two congeneric species as outgroups, were sequenced with mitochondrial and nuclear markers to resolve the discrepancy found between mtDNA divergence and present species‐level taxonomy. We found five COI sub‐clades within the E. autumnata complex, most of which were sympatric and with little geographic structure. Nuclear markers (ITS2 and Wingless) showed little variation and gave no indications that E. autumnata comprises more than one species. The samples were screened with primers for Wolbachia outer surface gene (wsp) and 12% of the samples tested positive. Two Wolbachia strains were associated with different mtDNA sub‐clades within E. autumnata, which may indicate indirect selection/selective sweeps on haplotypes. Our results demonstrate that deep mtDNA divergences are not synonymous with cryptic speciation and this has important implications for the use of mtDNA in species delimitation, like in DNA barcoding.     

Genetic variation in two populations of the rough‐skinned newt (Taricha granulosa) assessed using novel tetranucleotide microsatellite loci

Seven novel tetranucleotide microsatellite loci were identified from a partial genomic DNA library, enriched for GATA‐motif microsatellites, from the rough‐skinned newt (Taricha granulosa). All loci were polymorphic, and one displayed a high frequency null allele. A related species, T. rivularis, displays strong site fidelity and detectable population structure over small spatial scales, so we assessed genetic variation in two samples of T. granulosa separated by 16 km. Distributions of allele frequencies differ significantly between our two sites, but small values of F_ST and Rho_ST suggest that the populations are linked by a large amount of gene flow.     

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.     

Evidence for ecological speciation via a host shift in the holly leaf miner,Phytomyza glabricola (Diptera: Agromyzidae)

Evolutionary radiations have been well documented in plants and insects, and natural selection may often underly these radiations. If radiations are adaptive, the diversity of species could be due to ecological speciation in these lineages. Agromyzid flies exhibit patterns of repeated host‐associated radiations. We investigated whether host‐associated population divergence and evidence of divergent selection exist in the leaf miner Phytomyza glabricola on its sympatric host plants, the holly species, Ilex coriacea and I. glabra. Using AFLPs and nuclear sequence data, we found substantial genetic divergence between host‐associated populations of these flies throughout their geographic range. Genome scans using the AFLP data identified 13 loci under divergent selection, consistent with processes of ecological speciation. EF‐1α data suggest that I. glabra is the original host of P. glabricola and that I. coriacea is the novel host, but the AFLP data are ambiguous with regard to directionality of the host shift.