Considerable genetic diversity and structure despite narrow endemism and limited ecological specialization in the Hayden's ringlet, Coenonympha haydenii

Understanding the processes that underlie the development of population genetic structure is central to the study of evolution. Patterns of genetic structure, in turn, can reveal signatures of isolation by distance (IBD), barriers to gene flow, or even the genesis of speciation. However, it is unclear how severe range restriction might impact the processes that dominate the development of genetic structure. In narrow endemic species, is population structure likely to be adaptive in nature, or rather the result of genetic drift? In this study, we investigated patterns of genetic diversity and structure in the narrow endemic Hayden's ringlet butterfly. Specifically, we asked to what degree genetic structure in the Hayden's ringlet can be explained by IBD, isolation by resistance (IBR) (in the form of geographic or ecological barriers to migration between populations), and isolation by environment (in the form of differences in host plant availability and preference). We employed a genotyping-by-sequencing (GBS) approach coupled with host preference assays, Bayesian modelling, and population genomic analyses to answer these questions. Our results suggest that despite their restricted range, levels of genetic diversity in the Hayden's ringlet are comparable to those seen in more widespread butterfly species. Hayden's ringlets showed a strong preference for feeding on grasses relative to sedges, but neither larval preference nor potential host availability at sampling sites correlated with genetic structure. We conclude that geography, in the form of IBR and simple IBD, was the major driver of contemporary patterns of differentiation in this narrow endemic species.

Genetic variation within a stick insect species associated with community-level traits

Phenotypic variation within species can affect the ecological dynamics of populations and communities. Characterizing the genetic variation underlying such effects can help parse the roles of genetic evolution and plasticity in 'eco-evolutionary dynamics' and inform how genetic variation may shape patterns of evolution. Here we employ genome-wide association (GWA) methods in Timema cristinae stick insects and their co-occurring arthropod communities to identify genetic variation associated with community-level traits. Previous studies have shown that maladaptation (i.e., imperfect crypsis) of T. cristinae can reduce the abundance and species richness of other arthropods due to an increase in bird predation. Whether genetic variation that is independent from crypsis has similar effects is unknown and was tested here using genome-wide genotyping-by-sequencing data of stick insects, arthropod community information, and GWA mapping with Bayesian sparse linear mixed models. We find associations between genetic variation in stick insects and arthropod community traits. However, these associations disappeared when host-plant traits are accounted for. We thus use path analysis to disentangle interrelationships among stick-insect genetic variation, host-plant traits and community traits. This revealed that host-plant size has large effects on arthropod communities, while genetic variation in stick insects has a smaller, but still significant effect. Our findings demonstrate that: (1) genetic variation in a species can be associated with community-level traits, but that (2) interrelationships among multiple factors may need to be analyzed to disentangle whether such associations represent causal relationships. This work helps to build a framework for genomic studies of eco-evolutionary dynamics.

Divergent dynamics of sexual and habitat isolation at the transition between stick insect populations and species

Speciation is often viewed as a continuum along which populations diverge until they become reproductively-isolated species. However, such divergence may be heterogeneous, proceeding in fits and bursts, rather than being uniform and gradual. We show in Timema stick insects that one component of reproductive isolation evolves non-uniformly across this continuum, whereas another does not. Specifically, we use thousands of host-preference and mating trials to study habitat and sexual isolation among 42 pairs of taxa spanning a range of genomic differentiation and divergence time. We find that habitat isolation is uncoupled from genomic differentiation within species, but accumulates linearly with it between species. In contrast, sexual isolation accumulates linearly across the speciation continuum, and thus exhibits similar dynamics to morphological traits not implicated in reproductive isolation. The results show different evolutionary dynamics for different components of reproductive isolation and highlight a special relevance for species status in the process of speciation.

Thirty-six years of butterfly monitoring, snow cover, and plant productivity reveal negative impacts of warmer winters and increased productivity on montane species

Climate change is contributing to declines of insects through rising temperatures, altered precipitation patterns, and an increasing frequency of extreme events. The impacts of both gradual and sudden shifts in weather patterns are realized directly on insect physiology and indirectly through impacts on other trophic levels. Here, we investigated direct effects of seasonal weather on butterfly occurrences and indirect effects mediated by plant productivity using a temporally intensive butterfly monitoring dataset, in combination with high-resolution climate data and a remotely sensed indicator of plant primary productivity. Specifically, we used Bayesian hierarchical path analysis to quantify relationships between weather and weather-driven plant productivity on the occurrence of 94 butterfly species from three localities distributed across an elevational gradient. We found that snow pack exerted a strong direct positive effect on butterfly occurrence and that low snow pack was the primary driver of reductions during drought. Additionally, we found that plant primary productivity had a consistently negative effect on butterfly occurrence. These results highlight mechanisms of weather-driven declines in insect populations and the nuances of climate change effects involving snow melt, which have implications for ecological theories linking topographic complexity to ecological resilience in montane systems.

DNA methylation differences between stick insect ecotypes

Epigenetic mechanisms, such as DNA methylation, can influence gene regulation and affect phenotypic variation, raising the possibility that they contribute to ecological adaptation. Beginning to address this issue requires high-resolution sequencing studies of natural populations to pinpoint epigenetic regions of potential ecological and evolutionary significance. However, such studies are still relatively uncommon, especially in insects, and are mainly restricted to a few model organisms. Here, we characterize patterns of DNA methylation for natural populations of Timema cristinae adapted to two host plant species (i.e. ecotypes). By integrating results from sequencing of whole transcriptomes, genomes and methylomes, we investigate whether environmental, host and genetic differences of these stick insects are associated with methylation levels of cytosine nucleotides in the CpG context. We report an overall genome-wide methylation level for T. cristinae of ~14%, with methylation being enriched in gene bodies and impoverished in repetitive elements. Genome-wide DNA methylation variation was strongly positively correlated with genetic distance (relatedness), but also exhibited significant host-plant effects. Using methylome-environment association analysis, we pinpointed specific genomic regions that are differentially methylated between ecotypes, with these regions being enriched for genes with functions in membrane processes. The observed association between methylation variation and genetic relatedness, and with the ecologically important variable of host plant, suggests a potential role for epigenetic modification in T. cristinae adaptation. To substantiate such adaptive significance, future studies could test whether methylation can be transmitted across generations and the extent to which it responds to experimental manipulation in field and laboratory studies.

Quantitative Analyses of Coupling in Hybrid Zones

In hybrid zones, whether barrier loci experience selection mostly independently or as a unit depends on the ratio of selection to recombination as captured by the coupling coefficient. Theory predicts a sharper transition between an uncoupled and coupled system when more loci affect hybrid fitness. However, the extent of coupling in hybrid zones has rarely been quantified. Here, we use simulations to characterize the relationship between the coupling coefficient and variance in clines across genetic loci. We then reanalyze 25 hybrid zone data sets and find that cline variances and estimated coupling coefficients form a smooth continuum from high variance and weak coupling to low variance and strong coupling. Our results are consistent with low rates of hybridization and a strong genome-wide barrier to gene flow when the coupling coefficient is much greater than 1, but also suggest that this boundary might be approached gradually and at a near constant rate over time.

A stabilizing eco-evolutionary feedback loop in the wild

There is increasing evidence that evolutionary and ecological processes can operate on the same timescale1,2 (i.e., contemporary time). As such, evolution can be sufficiently rapid to affect ecological processes such as predation or competition. Thus, evolution can influence population, community, and ecosystem-level dynamics. Indeed, studies have now shown that evolutionary dynamics can alter community structure3,4,5,6 and ecosystem function.7,8,9,10 In turn, shifts in ecological dynamics driven by evolution might feed back to affect the evolutionary trajectory of individual species.11 This feedback loop, where evolutionary and ecological changes reciprocally affect one another, is a central tenet of eco-evolutionary dynamics.1,12 However, most work on such dynamics in natural populations has focused on one-way causal associations between ecology and evolution.13 Hence, direct empirical evidence for eco-evolutionary feedback is rare and limited to laboratory or mesocosm experiments.13,14,15,16 Here, we show in the wild that eco-evolutionary dynamics in a plant-feeding arthropod community involve a negative feedback loop. Specifically, adaptation in cryptic coloration in a stick-insect species mediates bird predation, with local maladaptation increasing predation. In turn, the abundance of arthropods is reduced by predation. Here, we experimentally manipulate arthropod abundance to show that these changes at the community level feed back to affect the stick-insect evolution. Specifically, low-arthropod abundance increases the strength of selection on crypsis, increasing local adaptation of stick insects in a negative feedback loop. Our results suggest that eco-evolutionary feedbacks are able to stabilize complex systems by preventing consistent directional change and therefore increasing resilience.

Complex evolutionary processes maintain an ancient chromosomal inversion

Genome re-arrangements such as chromosomal inversions are often involved in adaptation. As such, they experience natural selection, which can erode genetic variation. Thus, whether and how inversions can remain polymorphic for extended periods of time remains debated. Here we combine genomics, experiments, and evolutionary modeling to elucidate the processes maintaining an inversion polymorphism associated with the use of a challenging host plant (Redwood trees) in Timema stick insects. We show that the inversion is maintained by a combination of processes, finding roles for life-history trade-offs, heterozygote advantage, local adaptation to different hosts, and gene flow. We use models to show how such multi-layered regimes of balancing selection and gene flow provide resilience to help buffer populations against the loss of genetic variation, maintaining the potential for future evolution. We further show that the inversion polymorphism has persisted for millions of years and is not a result of recent introgression. We thus find that rather than being a nuisance, the complex interplay of evolutionary processes provides a mechanism for the long-term maintenance of genetic variation.

Early-warning signals of impending speciation

Species formation is a central topic in biology and a large body of theoretical work has explored the conditions under which speciation occurs, including whether speciation dynamics are gradual or abrupt. In some cases of abrupt speciation, differentiation slowly builds up until it reaches a threshold, at which point linkage disequilibrium (LD) and divergent selection enter a positive feedback loop that triggers accelerated change. Notably, such abrupt transitions powered by a positive feedback have also been observed in a range of other systems. Efforts to anticipate abrupt transitions have led to the development of 'early warning signals' (EWS), i.e. specific statistical patterns preceding abrupt transitions. Examples of EWS are rising autocorrelation and variance in time-series data, due to the reduction of the ability of the system to recover from disturbances. Here, we investigate whether speciation dynamics in theoretical models also exhibit EWS. Using a model of genetic divergence between two populations, we search for EWS before gradual and abrupt speciation events. We do so using six different metrics of differentiation: the effective migration rate, the number of selected loci, the mean fitness of our studied population, LD, FST and Dabs, a metric analogous to DXY . We find evidence for EWS, but with a heterogeneity in their strength among differentiation metrics. We specifically identify FST and the effective migration rate as the most reliable EWS of upcoming abrupt speciation events. Our results provide initial insights into potential EWS of impending speciation and contribute to efforts to generalize the mechanisms underlying EWS.

Disruptive selection and the evolution of discrete color morphs in Timema stick insects

A major unresolved issue in biology is why phenotypic and genetic variation is sometimes continuous, yet other times packaged into discrete units of diversity, such as morphs, ecotypes, and species. In theory, ecological discontinuities can impose strong disruptive selection that promotes the evolution of discrete forms, but direct tests of this hypothesis are lacking. Here, we show that Timema stick insects exhibit genetically determined color morphs that range from weakly to strongly discontinuous. Color data from nature and a manipulative field experiment demonstrate that greater morph differentiation is associated with shifts from host plants exhibiting more continuous color variation to those exhibiting greater coloration distance between green leaves and brown stems, the latter of which generates strong disruptive selection. Our results show how ecological factors can promote discrete variation, and we further present results on how this can have variable effects on the genetic differentiation that promotes speciation.

Population genomic evidence of selection on structural variants in a natural hybrid zone

Structural variants (SVs) can promote speciation by directly causing reproductive isolation or by suppressing recombination across large genomic regions. Whereas examples of each mechanism have been documented, systematic tests of the role of SVs in speciation are lacking. Here, we take advantage of long-read (Oxford nanopore) whole-genome sequencing and a hybrid zone between two Lycaeides butterfly taxa (L. melissa and Jackson Hole Lycaeides) to comprehensively evaluate genome-wide patterns of introgression for SVs and relate these patterns to hypotheses about speciation. We found >100,000 SVs segregating within or between the two hybridizing species. SVs and SNPs exhibited similar levels of genetic differentiation between species, with the exception of inversions, which were more differentiated. We detected credible variation in patterns of introgression among SV loci in the hybrid zone, with 562 of 1419 ancestry-informative SVs exhibiting genomic clines that deviated from null expectations based on genome-average ancestry. Overall, hybrids exhibited a directional shift towards Jackson Hole Lycaeides ancestry at SV loci, consistent with the hypothesis that these loci experienced more selection on average than SNP loci. Surprisingly, we found that deletions, rather than inversions, showed the highest skew towards excess ancestry from Jackson Hole Lycaeides. Excess Jackson Hole Lycaeides ancestry in hybrids was also especially pronounced for Z-linked SVs and inversions containing many genes. In conclusion, our results show that SVs are ubiquitous and suggest that SVs in general, but especially deletions, might disproportionately affect hybrid fitness and thus contribute to reproductive isolation.

Spatiotemporal and ontogenetic variation, microbial selection, and predicted Bd-inhibitory function in the skin-associated microbiome of a Rocky Mountain amphibian

Host-associated microbiomes play important roles in host health and pathogen defense. In amphibians, the skin-associated microbiota can contribute to innate immunity with potential implications for disease management. Few studies have examined season-long temporal variation in the amphibian skin-associated microbiome, and the interactions between bacteria and fungi on amphibian skin remain poorly understood. We characterize season-long temporal variation in the skin-associated microbiome of the western tiger salamander (Ambystoma mavortium) for both bacteria and fungi between sites and across salamander life stages. Two hundred seven skin-associated microbiome samples were collected from salamanders at two Rocky Mountain lakes throughout the summer and fall of 2018, and 127 additional microbiome samples were collected from lake water and lake substrate. We used 16S rRNA and ITS amplicon sequencing with Bayesian Dirichlet-multinomial regression to estimate the relative abundances of bacterial and fungal taxa, test for differential abundance, examine microbial selection, and derive alpha diversity. We predicted the ability of bacterial communities to inhibit the amphibian chytrid fungus Batrachochytrium dendrobatidis (Bd), a cutaneous fungal pathogen, using stochastic character mapping and a database of Bd-inhibitory bacterial isolates. For both bacteria and fungi, we observed variation in community composition through time, between sites, and with salamander age and life stage. We further found that temporal trends in community composition were specific to each combination of salamander age, life stage, and lake. We found salamander skin to be selective for microbes, with many taxa disproportionately represented relative to the environment. Salamander skin appeared to select for predicted Bd-inhibitory bacteria, and we found a negative relationship between the relative abundances of predicted Bd-inhibitory bacteria and Bd. We hope these findings will assist in the conservation of amphibian species threatened by chytridiomycosis and other emerging diseases.

Laplace's demon in biology: Models of evolutionary prediction

Our ability to predict natural phenomena can be limited by incomplete information. This issue is exemplified by "Laplace's demon," an imaginary creature proposed in the 18th century, who knew everything about everything, and thus could predict the full nature of the universe forward or backward in time. Quantum mechanics, among other things, has cast doubt on the possibility of Laplace's demon in the full sense, but the idea still serves as a useful metaphor for thinking about the extent to which prediction is limited by incomplete information on deterministic processes versus random factors. Here, we use simple analytical models and computer simulations to illustrate how data limits can be captured in a Bayesian framework, and how they influence our ability to predict evolution. We show how uncertainty in measurements of natural selection, or low predictability of external environmental factors affecting selection, can greatly reduce predictive power, often swamping the influence of intrinsic randomness caused by genetic drift. Thus, more accurate knowledge concerning the causes and action of natural selection is key to improving prediction. Fortunately, our analyses and simulations show quantitatively that reasonable improvements in data quantity and quality can meaningfully increase predictability.

Climatic similarity and genomic background shape the extent of parallel adaptation in Timema stick insects

Evolution can repeat itself, resulting in parallel adaptations in independent lineages occupying similar environments. Moreover, parallel evolution sometimes, but not always, uses the same genes. Two main hypotheses have been put forth to explain the probability and extent of parallel evolution. First, parallel evolution is more likely when shared ecologies result in similar patterns of natural selection in different taxa. Second, parallelism is more likely when genomes are similar because of shared standing variation and similar mutational effects in closely related genomes. Here we combine ecological, genomic, experimental and phenotypic data with Bayesian modelling and randomization tests to quantify the degree of parallelism and its relationship with ecology and genetics. Our results show that the extent to which genomic regions associated with climate are parallel among species of Timema stick insects is shaped collectively by shared ecology and genomic background. Specifically, the extent of genomic parallelism decays with divergence in climatic conditions (that is, habitat or ecological similarity) and genomic similarity. Moreover, we find that climate-associated loci are likely subject to selection in a field experiment, overlap with genetic regions associated with cuticular hydrocarbon traits and are not strongly shaped by introgression between species. Our findings shed light on when evolution is most expected to repeat itself.

Evidence that genomic incompatibilities and other multilocus processes impact hybrid fitness in a rattlesnake hybrid zone

Hybrid zones provide valuable opportunities to understand the genomic mechanisms that promote speciation by providing insight into factors involved in intermediate stages of speciation. Here, we investigate introgression in a hybrid zone between two rattlesnake species (Crotalus viridis and Crotalus oreganus concolor) that have undergone historical allopatric divergence and recent range expansion and secondary contact. We use Bayesian genomic cline models to characterize genomic patterns of introgression between these lineages and identify loci potentially subject to selection in hybrids. We find evidence for a large number of genomic regions with biased ancestry that deviate from the genomic background in hybrids (i.e., excess ancestry loci), which tend to be associated with genomic regions with higher recombination rates. We also identify suites of excess ancestry loci that show highly correlated allele frequencies (including conspecific and heterospecific combinations) across physically unlinked genomic regions in hybrids. Our findings provide evidence for multiple multilocus evolutionary processes impacting hybrid fitness in this system.

A continental-scale survey of Wolbachia infections in blue butterflies reveals evidence of interspecific transfer and invasion dynamics

Infections by maternally inherited bacterial endosymbionts, especially Wolbachia, are common in insects and other invertebrates but infection dynamics across species ranges are largely under studied. Specifically, we lack a broad understanding of the origin of Wolbachia infections in novel hosts, and the historical and geographical dynamics of infections that are critical for identifying the factors governing their spread. We used Genotype-by-Sequencing data from previous population genomics studies for range-wide surveys of Wolbachia presence and genetic diversity in North American butterflies of the genus Lycaeides. As few as one sequence read identified by assembly to a Wolbachia reference genome provided high accuracy in detecting infections in host butterflies as determined by confirmatory PCR tests, and maximum accuracy was achieved with a threshold of only 5 sequence reads per host individual. Using this threshold, we detected Wolbachia in all but 2 of the 107 sampling localities spanning the continent, with infection frequencies within populations ranging from 0% to 100% of individuals, but with most localities having high infection frequencies (mean = 91% infection rate). Three major lineages of Wolbachia were identified as separate strains that appear to represent 3 separate invasions of Lycaeides butterflies by Wolbachia. Overall, we found extensive evidence for acquisition of Wolbachia through interspecific transfer between host lineages. Strain wLycC was confined to a single butterfly taxon, hybrid lineages derived from it, and closely adjacent populations in other taxa. While the other 2 strains were detected throughout the rest of the continent, strain wLycB almost always co-occurred with wLycA. Our demographic modeling suggests wLycB is a recent invasion. Within strain wLycA, the 2 most frequent haplotypes are confined almost exclusively to separate butterfly taxa with haplotype A1 observed largely in Lycaeides melissa and haplotype A2 observed most often in Lycaeides idas localities, consistent with either cladogenic mode of infection acquisition from a common ancestor or by hybridization and accompanying mutation. More than 1 major Wolbachia strain was observed in 15 localities. These results demonstrate the utility of using resequencing data from hosts to quantify Wolbachia genetic variation and infection frequency and provide evidence of multiple colonizations of novel hosts through hybridization between butterfly lineages and complex dynamics between Wolbachia strains.

The short-term, genome-wide effects of indirect selection deserve study: A response to Charlesworth and Jensen (2022)

We recently published a paper quantifying the genome-wide consequences of natural selection, including the effects of indirect selection due to the correlation of genetic regions (neutral or selected) with directly selected regions (Gompert et al., 2022). In their critique of our paper, Charlesworth and Jensen (2022) make two main points: (i) indirect selection is equivalent to hitchhiking and thus well documented (i.e., our results are not novel) and (ii) that we do not demonstrate the source of linkage disequilibrium (LD) between SNPs and the Mel-Stripe locus in the Timema cristinae experiment we analyse. As we discuss in detail below, neither of these are substantial criticisms of our work.

Testing for fitness epistasis in a transplant experiment identifies a candidate adaptive locus in Timema stick insects

Identifying the genetic basis of adaptation is a central goal of evolutionary biology. However, identifying genes and mutations affecting fitness remains challenging because a large number of traits and variants can influence fitness. Selected phenotypes can also be difficult to know a priori, complicating top-down genetic approaches for trait mapping that involve crosses or genome-wide association studies. In such cases, experimental genetic approaches, where one maps fitness directly and attempts to infer the traits involved afterwards, can be valuable. Here, we re-analyse data from a transplant experiment involving Timema stick insects, where five physically clustered single-nucleotide polymorphisms associated with cryptic body coloration were shown to interact to affect survival. Our analysis covers a larger genomic region than past work and revealed a locus previously not identified as associated with survival. This locus resides near a gene, Punch (Pu), involved in pteridine pigments production, implying that it could be associated with an unmeasured coloration trait. However, by combining previous and newly obtained phenotypic data, we show that this trait is not eye or body coloration. We discuss the implications of our results for the discovery of traits, genes and mutations associated with fitness in other systems, as well as for supergene evolution. This article is part of the theme issue 'Genetic basis of adaptation and speciation: from loci to causative mutations'.

Frequency dependence and the predictability of evolution in a changing environment

Frequency‐dependent (FD) selection, whereby fitness and selection depend on the genetic or phenotypic composition of the population, arises in numerous ecological contexts (competition, mate choice, crypsis, mimicry, etc.) and can strongly impact evolutionary dynamics. In particular, negative frequency‐dependent selection (NFDS) is well known for its ability to potentially maintain stable polymorphisms, but it has also been invoked as a source of persistent, predictable frequency fluctuations. However, the conditions under which such fluctuations persist are not entirely clear. In particular, previous work rarely considered that FD is unlikely to be the sole driver of evolutionary dynamics when it occurs, because most environments are not static but instead change dynamically over time. Here, we investigate how FD interacts with a temporally fluctuating environment to shape the dynamics of population genetic change. We show that a simple metric introduced by Lewontin, the slope of frequency change against frequency near equilibrium, works as a key criterion for distinguishing microevolutionary outcomes, even in a changing environment. When this slope D is between 0 and –2 (consistent with the empirical examples we review), substantial fluctuations would not persist on their own in a large population occupying a constant environment, but they can still be maintained indefinitely as quasi‐cycles fueled by environmental noise or genetic drift. However, such moderate NFDS buffers and temporally shifts evolutionary responses to periodic environments (e.g., seasonality). Stronger FD, with slope D < –2, can produce self‐sustained cycles that may overwhelm responses to a changing environment, or even chaos that fundamentally limits predictability. This diversity of expected outcomes, together with the empirical evidence for both FD and environment‐dependent selection, suggests that the interplay of internal dynamics with external forcing should be investigated more systematically to reach a better understanding and prediction of evolution.

Natural selection drives genome-wide evolution via chance genetic associations

Understanding selection's impact on the genome is a major theme in biology. Functionally neutral genetic regions can be affected indirectly by natural selection, via their statistical association with genes under direct selection. The genomic extent of such indirect selection, particularly across loci not physically linked to those under direct selection, remains poorly understood, as does the time scale at which indirect selection occurs. Here, we use field experiments and genomic data in stick insects, deer mice and stickleback fish to show that widespread statistical associations with genes known to affect fitness cause many genetic loci across the genome to be impacted indirectly by selection. This includes regions physically distant from those directly under selection. Then, focusing on the stick insect system, we show that statistical associations between SNPs and other unknown, causal variants result in additional indirect selection in general and specifically within genomic regions of physically linked loci. This widespread indirect selection necessarily makes aspects of evolution more predictable. Thus, natural selection combines with chance genetic associations to affect genome-wide evolution across linked and unlinked loci and even in modest-sized populations. This process has implications for the application of evolutionary principles in basic and applied science.

Biodiversity, resilience and the stability of evolutionary systems

Various macro-evolutionary phenomena, such as long-term stability punctuated by bursts of evolution, are difficult to explain via the micro-evolutionary process of weak selection acting steadily on individual mutations. In contrast, bursts of change are expected if evolutionary systems are complex and balanced, with occasional disruption of balance. Such disruption represents the collapse of resilience, akin to the snapping of an elastic band. It can be driven by external factors, or by self-propagating feedback loops internal to a system. Thus, evolutionary resilience could help explain how evolution generates broader patterns of biodiversity. We outline evidence and tests for this hypothesis, which emphasizes the processes balancing evolution, as urged fifty years ago in ecological genetics and via modern results in a range of systems.

Disparate genetic variants associated with distinct components of cowpea resistance to the seed beetle Callosobruchus maculatus

Polygenic genome-wide association mapping identified two regions of the cowpea genome associated with different components of resistance to its major post-harvest pest, the seed beetle Callosobruchus maculatus. Cowpea (Vigna unguiculata) is an important grain and fodder crop in arid and semi-arid regions of Africa, Asia, and South America, where the cowpea seed beetle, Callosobruchus maculatus, is a serious post-harvest pest. Development of cultivars resistant to C. maculatus population growth in storage could increase grain yield and quality and reduce reliance on insecticides. Here, we use a MAGIC (multi-parent, advanced-generation intercross) population of cowpea consisting of 305 recombinant inbred lines (RILs) to identify genetic variants associated with resistance to seed beetles. Because inferences regarding the genetic basis of resistance may depend on the source of the pest or the assay protocol, we used two divergent geographic populations of C. maculatus and two complementary assays to measure several aspects of resistance. Using polygenic genome-wide association mapping models, we found that the cowpea RILs harbor substantial additive-genetic variation for most resistance measures. Variation in several components of resistance, including larval development time and survival, was largely explained by one or several linked loci on chromosome 5. A second region on chromosome 8 explained increased seed resistance via the induction of early-exiting larvae. Neither of these regions contained genes previously associated with resistance to insects that infest grain legumes. We found some evidence of gene-gene interactions affecting resistance, but epistasis did not contribute substantially to resistance variation in this mapping population. The combination of mostly high heritabilities and a relatively consistent and simple genetic architecture increases the feasibility of breeding for enhanced resistance to C. maculatus.

Genomic time-series data show that gene flow maintains high genetic diversity despite substantial genetic drift in a butterfly species

Effective population size affects the efficacy of selection, rate of evolution by drift, and neutral diversity levels. When species are subdivided into multiple populations connected by gene flow, evolutionary processes can depend on global or local effective population sizes. Theory predicts that high levels of diversity might be maintained by gene flow, even very low levels of gene flow, consistent with species long-term effective population size, but tests of this idea are mostly lacking. Here, we show that Lycaeides buttery populations maintain low contemporary (variance) effective population sizes (e.g., ~200 individuals) and thus evolve rapidly by genetic drift. In contrast, populations harbored high levels of genetic diversity consistent with an effective population size several orders of magnitude larger. We hypothesized that the differences in the magnitude and variability of contemporary versus long-term effective population sizes were caused by gene flow of sufficient magnitude to maintain diversity but only subtly affect evolution on generational time scales. Consistent with this hypothesis, we detected low but non-trivial gene flow among populations. Furthermore, using short-term population-genomic time-series data, we documented patterns consistent with predictions from this hypothesis, including a weak but detectable excess of evolutionary change in the direction of the mean (migrant gene pool) allele frequencies across populations, and consistency in the direction of allele frequency change over time. The documented decoupling of diversity levels and short-term change by drift in Lycaeides has implications for our understanding of contemporary evolution and the maintenance of genetic variation in the wild.

Inversion breakpoints and the evolution of supergenes

The coexistence of discrete morphs that differ in multiple traits is common within natural populations of many taxa. Such morphs are often associated with chromosomal inversions, presumably because the recombination suppressing effects of inversions help maintain alternate adaptive combinations of alleles across the multiple loci affecting these traits. However, inversions can also harbour selected mutations at their breakpoints, leading to their rise in frequency in addition to (or independent from) their role in recombination suppression. In this review, we first describe the different ways that breakpoints can create mutations. We then critically examine the evidence for the breakpoint-mutation and recombination suppression hypotheses for explaining the existence of discrete morphs associated with chromosomal inversions. We find that the evidence that inversions are favoured due to recombination suppression is often indirect. The evidence that breakpoints harbour mutations that are adaptive is also largely indirect, with the characterization of inversion breakpoints at the sequence level being incomplete in most systems. Direct tests of the role of suppressed recombination and breakpoint mutations in inversion evolution are thus needed. Finally, we emphasize how the two hypotheses of recombination suppression and breakpoint mutation can act in conjunction, with implications for understanding the emergence of supergenes and their evolutionary dynamics. We conclude by discussing how breakpoint characterization could improve our understanding of complex, discrete phenotypic forms in nature.

Standing genetic variation in laboratory populations of insecticide-susceptible Phlebotomus papatasi and Lutzomyia longipalpis (Diptera: Psychodidae: Phlebotominae) for the evolution of resistance

Insecticides can exert strong selection on insect pest species, including those that vector diseases, and have led to rapid evolution of resistance. Despite such rapid evolution, relatively little is known about standing genetic variation for resistance in insecticide-susceptible populations of many species. To help fill this knowledge gap, we generated genotyping-by-sequencing data from insecticide-susceptible Phlebotomus papatasi and Lutzomyia longipalpis sand flies that survived or died from a sub-diagnostic exposure to either permethrin or malathion using a modified version of the Centers for Disease Control and Prevention bottle bioassay. Multi-locus genome-wide association mapping methods were used to quantify standing genetic variation for insecticide resistance in these populations and to identify specific alleles associated with insecticide survival. For each insecticide treatment, we estimated the proportion of the variation in survival explained by the genetic data (i.e., "chip" heritability) and the number and contribution of individual loci with measurable effects. For all treatments, survival to an insecticide exposure was heritable with a polygenic architecture. Both P. papatasi and L. longipalpis had alleles for survival that resided within many genes throughout their genomes. The implications for resistance conferred by many alleles, as well as inferences made about the utility of laboratory insecticide resistance association studies compared to field observations, are discussed.

A population-genomic approach for estimating selection on polygenic traits in heterogeneous environments

Strong selection can cause rapid evolutionary change, but temporal fluctuations in the form, direction and intensity of selection can limit net evolutionary change over longer time periods. Fluctuating selection could affect molecular diversity levels and the evolution of plasticity and ecological specialization. Nonetheless, this phenomenon remains understudied, in part because of analytical limitations and the general difficulty of detecting selection that does not occur in a consistent manner. Herein, I fill this analytical gap by presenting an approximate Bayesian computation (ABC) method to detect and quantify fluctuating selection on polygenic traits from population genomic time-series data. I propose a model for environment-dependent phenotypic selection. The evolutionary genetic consequences of selection are then modelled based on a genotype-phenotype map. Using simulations, I show that the proposed method generates accurate and precise estimates of selection when the generative model for the data is similar to the model assumed by the method. The performance of the method when applied to an evolve-and-resequence study of host adaptation in the cowpea seed beetle (Callosobruchus maculatus) was more idiosyncratic and depended on specific analytical choices. Despite some limitations, these results suggest the proposed method provides a powerful approach to connect the causes of (variable) selection to traits and genome-wide patterns of evolution. Documentation and open-source computer software (fsabc) implementing this method are available from github (https://github.com/zgompert/fsabc.git).

Model-based genotype and ancestry estimation for potential hybrids with mixed-ploidy

Non-random mating among individuals can lead to spatial clustering of genetically similar individuals and population stratification. This deviation from panmixia is commonly observed in natural populations. Consequently, individuals can have parentage in single populations or involving hybridization between differentiated populations. Accounting for this mixture and structure is important when mapping the genetics of traits and learning about the formative evolutionary processes that shape genetic variation among individuals and populations. Stratified genetic relatedness among individuals is commonly quantified using estimates of ancestry that are derived from a statistical model. Development of these models for polyploid and mixed-ploidy individuals and populations has lagged behind those for diploids. Here, we extend and test a hierarchical Bayesian model, called entropy, which can use low-depth sequence data to estimate genotype and ancestry parameters in autopolyploid and mixed-ploidy individuals (including sex chromosomes and autosomes within individuals). Our analysis of simulated data illustrated the trade-off between sequencing depth and genome coverage and found lower error associated with low-depth sequencing across a larger fraction of the genome than with high-depth sequencing across a smaller fraction of the genome. The model has high accuracy and sensitivity as verified with simulated data and through analysis of admixture among populations of diploid and tetraploid Arabidopsis arenosa.

Measuring the effect of environmental stress on inbreeding depression alone obscures the relative importance of inbreeding-stress interactions on overall fitness in Callosobruchus maculatus

Environmental stress can have a profound effect on inbreeding depression. Quantifying this effect is of particular importance in threatened populations, which are often simultaneously subject to both inbreeding and environmental stress. But while the prevalence of inbreeding-stress interactions is well known, the importance and broader applicability of such interactions in conservation are not clearly understood. We used seed beetles, Callosobruchus maculatus, as a model system to quantify how environmental stressors (here host quality and temperature stress) interact with inbreeding as measured by changes in the magnitude of inbreeding depression, δ, as well as the relative importance of inbreeding-stress interactions to overall fitness. We found that while both environmental stressors caused substantial inbreeding-stress interactions as measured by change in δ, the relative importance of these interactions to overall survival was modest. This suggests that assessing inbreeding-stress interactions within the framework of δ alone may give an inaccurate representation of the relevance of interactions to population persistence. Furthermore, we found that the effect of environmental stress on fitness, but not inbreeding depression, varied strongly among populations. These results suggest that the outcomes of inbreeding-stress interactions are not easily generalized, an important consideration in conservation settings.

Predicting patch occupancy reveals the complexity of host range expansion

Specialized plant-insect interactions are a defining feature of life on earth, yet we are only beginning to understand the factors that set limits on host ranges in herbivorous insects. To better understand the recent adoption of alfalfa as a host plant by the Melissa blue butterfly, we quantified arthropod assemblages and plant metabolites across a wide geographic region while controlling for climate and dispersal inferred from population genomic variation. The presence of the butterfly is successfully predicted by direct and indirect effects of plant traits and interactions with other species. Results are consistent with the predictions of a theoretical model of parasite host range in which specialization is an epiphenomenon of the many barriers to be overcome rather than a consequence of trade-offs in developmental physiology.

Genetic and phenotypic evidence of a contact zone between divergent colour morphs of the iconic red-eyed treefrog

Hybrid zones act as natural laboratories where divergent genomes interact, providing powerful systems for examining the evolutionary processes underlying biological diversity. In this study, we characterized patterns of genomic and phenotypic variation resulting from hybridization between divergent intraspecific lineages of the Neotropical red-eyed treefrog (Agalychnis callidryas). We found genetic evidence of a newly discovered contact zone and phenotypic novelty in leg colour-a trait suspected to play a role in mediating assortative mating in this species. Analysis of hybrid ancestry revealed an abundance of later-generation F_n individuals, suggesting persistence of hybrids in the contact zone. Hybrids are predominantly of southern ancestry but are phenotypically more similar to northern populations. Genome-wide association mapping revealed QTL with measurable effects on leg-colour variation, but further work is required to dissect the architecture of this trait and establish causal links. Further, genomic cline analyses indicated substantial variation in patterns of introgression across the genome. Directional introgression of loci associated with different aspects of leg colour are inherited from each parental lineage, creating a distinct hybrid colour pattern. We show that hybridization can generate new phenotypes, revealing the evolutionary processes that potentially underlie patterns of phenotypic diversity in this iconic polytypic frog. Our study is consistent with a role of hybridization and sexual selection in lineage diversification, evolutionary processes that have been implicated in accelerating divergence in the most phenotypically diverse species.

Colonization of Marginal Host Plants by Seed Beetles (Coleoptera: Chrysomelidae): Effects of Geographic Source and Genetic Admixture

The ability to adapt to a novel host plant may vary among insect populations with different genetic histories, and colonization of a marginal host may be facilitated by genetic admixture of disparate populations. We assembled populations of the seed beetle, Callosobruchus maculatus (F.), from four continents, and compared their ability to infest two hosts, lentil and pea. We also formed two cross-continent hybrids (Africa × N.A. and Africa × S.A.). In pre-selection assays, survival was only ~3% in lentil and ~40% in pea. For three replicate populations per line, colonization success on lentil was measured as cumulative exit holes after 75-175 d. On pea, we estimated the change in larval survival after five generations of selection. Females in all lines laid few eggs on lentil, and survival of F1 larvae was uniformly <5%. Subsequently, however, the lines diverged considerably in population growth. Performance on lentil was highest in the Africa × N.A. hybrid, which produced far more adults (mean > 11,000) than either parental line. At the other extreme, Asian populations on lentil appeared to have gone extinct. The Africa × N.A. line also exhibited the highest survival on pea, and again performed better than either parent line. However, no line displayed a rapid increase in survival on pea, as is sometimes observed on lentil. Our results demonstrate that geographic populations can vary substantially in their responses to the same novel resource. In addition, genetic admixtures (potentially caused by long-distance transport of infested seeds) may facilitate colonization of an initially poor host.

Recent hybrids recapitulate ancient hybrid outcomes

Genomic outcomes of hybridization depend on selection and recombination in hybrids. Whether these processes have similar effects on hybrid genome composition in contemporary hybrid zones versus ancient hybrid lineages is unknown. Here we show that patterns of introgression in a contemporary hybrid zone in Lycaeides butterflies predict patterns of ancestry in geographically adjacent, older hybrid populations. We find a particularly striking lack of ancestry from one of the hybridizing taxa, Lycaeides melissa, on the Z chromosome in both the old and contemporary hybrids. The same pattern of reduced L. melissa ancestry on the Z chromosome is seen in two other ancient hybrid lineages. More generally, we find that patterns of ancestry in old or ancient hybrids are remarkably predictable from contemporary hybrids, which suggests selection and recombination affect hybrid genomes in a similar way across disparate time scales and during distinct stages of speciation and species breakdown.

Combining Experimental Evolution and Genomics to Understand How Seed Beetles Adapt to a Marginal Host Plant

Genes that affect adaptive traits have been identified, but our knowledge of the genetic basis of adaptation in a more general sense (across multiple traits) remains limited. We combined population-genomic analyses of evolve-and-resequence experiments, genome-wide association mapping of performance traits, and analyses of gene expression to fill this knowledge gap and shed light on the genomics of adaptation to a marginal host (lentil) by the seed beetle Callosobruchus maculatus. Using population-genomic approaches, we detected modest parallelism in allele frequency change across replicate lines during adaptation to lentil. Mapping populations derived from each lentil-adapted line revealed a polygenic basis for two host-specific performance traits (weight and development time), which had low to modest heritabilities. We found less evidence of parallelism in genotype-phenotype associations across these lines than in allele frequency changes during the experiments. Differential gene expression caused by differences in recent evolutionary history exceeded that caused by immediate rearing host. Together, the three genomic datasets suggest that genes affecting traits other than weight and development time are likely to be the main causes of parallel evolution and that detoxification genes (especially cytochrome P450s and beta-glucosidase) could be especially important for colonization of lentil by C. maculatus.

Contrasting genomic and phenotypic outcomes of hybridization between pairs of mimetic butterfly taxa across a suture zone

Hybrid zones, whereby divergent lineages come into contact and eventually hybridize, can provide insights on the mechanisms involved in population differentiation and reproductive isolation, and ultimately speciation. Suture zones offer the opportunity to compare these processes across multiple species. In this paper we use reduced-complexity genomic data to compare the genetic and phenotypic structure and hybridization patterns of two mimetic butterfly species, Ithomia salapia and Oleria onega (Nymphalidae: Ithomiini), each consisting of a pair of lineages differentiated for their wing colour pattern and that come into contact in the Andean foothills of Peru. Despite similarities in their life history, we highlight major differences, both at the genomic and phenotypic level, between the two species. These differences include the presence of hybrids, variations in wing phenotype, and genomic patterns of introgression and differentiation. In I. salapia, the two lineages appear to hybridize only rarely, whereas in O. onega the hybrids are not only more common, but also genetically and phenotypically more variable. We also detected loci statistically associated with wing colour pattern variation, but in both species these loci were not over-represented among the candidate barrier loci, suggesting that traits other than wing colour pattern may be important for reproductive isolation. Our results contrast with the genomic patterns observed between hybridizing lineages in the mimetic Heliconius butterflies, and call for a broader investigation into the genomics of speciation in Ithomiini - the largest radiation of mimetic butterflies.

The role of structural genomic variants in population differentiation and ecotype formation in Timema cristinae walking sticks

Theory predicts that structural genomic variants such as inversions can promote adaptive diversification and speciation. Despite increasing empirical evidence that adaptive divergence can be triggered by one or a few large inversions, the degree to which widespread genomic regions under divergent selection are associated with structural variants remains unclear. Here we test for an association between structural variants and genomic regions that underlie parallel host-plant-associated ecotype formation in Timema cristinae stick insects. Using mate-pair resequencing of 20 new whole genomes we find that moderately sized structural variants such as inversions, deletions and duplications are widespread across the genome, being retained as standing variation within and among populations. Using 160 previously published, standard-orientation whole genome sequences we find little to no evidence that the DNA sequences within inversions exhibit accentuated differentiation between ecotypes. In contrast, a formerly described large region of reduced recombination that harbours genes controlling colour-pattern exhibits evidence for accentuated differentiation between ecotypes, which is consistent with differences in the frequency of colour-pattern morphs between host-associated ecotypes. Our results suggest that some types of structural variants (e.g., large inversions) are more likely to underlie adaptive divergence than others, and that structural variants are not required for subtle yet genome-wide genetic differentiation with gene flow.

Components of Cowpea Resistance to the Seed Beetle Callosobruchus maculatus (Coleoptera: Chrysomelidae: Bruchinae)

Cowpea, Vigna unguiculata (L.) Walp., serves as a major source of dietary protein in many tropical and subtropical regions around the world. To identify loci associated with agronomically desirable traits, eight elite cowpea cultivars were systematically inter-crossed for eight generations to yield 305 recombinant inbred lines. Here, we investigated whether these founder parents also possess resistance to the seed beetle Callosobruchus maculatus (F.), a highly destructive post-harvest pest. We estimated larval survival in seeds, egg-to-adult development time, adult mass at emergence, and seed acceptance for oviposition. Survival varied significantly among cowpea cultivars, but the pattern was complicated by an unexpected source of mortality; on three cultivars, mature larvae in a substantial fraction of seeds (20-36%) exited seeds prematurely, and consequently failed to molt into viable adults. Even if such seeds were eliminated from the analysis, survival in the remaining seeds varied from 49 to 92% across the eight parents. Development time and body mass also differed among hosts, with particularly slow larval development on three closely related cultivars. Egg-laying females readily accepted all cultivars except one with a moderately rugose seed coat. Overall, suitability ranks of the eight cultivars depended on beetle trait; a cultivar that received the most eggs (IT82E-18) also conferred low survival. However, one cultivar (IT93K-503-1) was a relatively poor host for all traits. Given the magnitude of variation among parental cultivars, future assays of genotyped recombinant progeny can identify genomic regions and candidate genes associated with resistance to seed beetles.

Extreme heterogeneity of population response to climatic variation and the limits of prediction

Certain general facets of biotic response to climate change, such as shifts in phenology and geographic distribution, are well characterized; however, it is not clear whether the observed similarity of responses across taxa will extend to variation in other population-level processes. We examined population response to climatic variation using long-term incidence data (collected over 42 years) encompassing 149 butterfly species and considerable habitat diversity (10 sites along an elevational gradient from sea level to over 2,700 m in California). Population responses were characterized by extreme heterogeneity that was not attributable to differences in species composition among sites. These results indicate that habitat heterogeneity might be a buffer against climate change and highlight important questions about mechanisms maintaining interpopulation differences in responses to weather. Despite overall heterogeneity of response, population dynamics were accurately predicted by our model for many species at each site. However, the overall correlation between observed and predicted incidence in a cross validation analysis was moderate (Pearson's r = 0.23, SE 0.01), and 97% of observed data fell within the predicted 95% credible intervals. Prediction was most successful for more abundant species as well as for sites with lower annual turnover. Population-level heterogeneity in response to climate variation and the limits of our predictive power highlight the challenges for a future of increasing climatic variability.

Vertical differentiation in tropical forest butterflies: a novel mechanism generating insect diversity?

Many tropical fruit-feeding nymphalid butterflies are associated with either the forest canopy or the understorey; however, the exceptions offer insights into the origins of tropical diversity. As it occurs in both habitats of tropical forests in Ecuador and Peru, Archaeoprepona demophon is one such exception. We compared patterns of occurrence of A. demophon in the canopy and understorey and population genomic variation for evidence of ecological and genetic differentiation between habitats. We found that butterfly occurrences in the canopy were largely uncorrelated with occurrences in the understorey at both localities, indicating independent demographic patterns in the two habitats. We also documented modest, significant genome-level differentiation at both localities. Genetic differentiation between habitat types (separated by approx. 20 m in elevation) was comparable to levels of differentiation between sampling locations (approx. 1500 km). We conclude that canopy and understorey populations of A. demophon represent incipient independent evolutionary units. These findings support the hypothesis that divergence between canopy and understorey-associated populations might be a mechanism generating insect diversity in the tropics.

Long-term experimental hybridisation results in the evolution of a new sex chromosome in swordtail fish

The remarkable diversity of sex determination mechanisms known in fish may be fuelled by exceptionally high rates of sex chromosome turnovers or transitions. However, the evolutionary causes and genomic mechanisms underlying this variation and instability are yet to be understood. Here we report on an over 30-year evolutionary experiment in which we tested the genomic consequences of hybridisation and selection between two Xiphophorus fish species with different sex chromosome systems. We find that introgression and imposing selection for pigmentation phenotypes results in the retention of an unexpectedly large maternally derived genomic region. During the hybridisation process, the sex-determining region of the X chromosome from one parental species was translocated to an autosome in the hybrids leading to the evolution of a new sex chromosome. Our results highlight the complexity of factors contributing to patterns observed in hybrid genomes, and we experimentally demonstrate that hybridisation can catalyze rapid evolution of a new sex chromosome.

Fish have a high diversity of sex-determining systems, but the mechanisms responsible for this are not well understood. Here, Franchini et al. show how hybridization and backcrossing have led to the evolution of a new sex chromosome in swordtail fish during 30 years of experimental evolution.

Variable Responses to Novel Hosts by Populations of the Seed Beetle Callosobruchus maculatus (Coleoptera: Chrysomelidae: Bruchinae)

Cosmopolitan pests can consist of geographic populations that differ in their current host ranges or in their ability to colonize a novel host. We compared the responses of cowpea-adapted, seed-beetle populations (Callosobruchus maculatus [F.] (Coleoptera: Chrysomelidae: Bruchinae)) from Africa, North America, and South America to four novel legumes: chickpea, lentil, mung bean, and pea. We also qualitatively compared these results to those obtained earlier for an Asian population. For each host, we measured larval survival to adult emergence and used both no-choice and choice tests to estimate host acceptance. The pattern of larval survival was similar among populations: high or moderately high survival on cowpea, mung bean, and chickpea, intermediate survival on pea, and very low survival on lentil. One exception was unusually high survival of African larvae on pea, and there was modest variation among populations for survival on lentil. The African population was also an outlier with respect to host acceptance; under no-choice conditions, African females showed a much greater propensity to accept the two least preferred hosts, chickpea and lentil. However, greater acceptance of these hosts by African females was not evident in choice tests. Inferences about population differences in host acceptance can thus strongly depend on experimental protocol. Future selection experiments can be used to determine whether the observed population differences in initial performance will affect the probability of producing self-sustaining populations on a marginal crop host.

Transitions from Single- to Multi-Locus Processes during Speciation with Gene Flow

During speciation-with-gene-flow, a transition from single-locus to multi-locus processes can occur, as strong coupling of multiple loci creates a barrier to gene flow. Testing predictions about such transitions with empirical data requires building upon past theoretical work and the continued development of quantitative approaches. We simulated genomes under several evolutionary scenarios of gene flow and divergent selection, extending previous work with the additions of neutral sites and coupling statistics. We used these simulations to investigate, in a preliminary way, if and how selected and neutral sites differ in the conditions they require for transitions during speciation. For the parameter combinations we explored, as the per-locus strength of selection grew and/or migration decreased, it became easier for selected sites to show divergence-and thus to rise in linkage disequilibrium (LD) with each other as a statistical consequence-farther in advance of the conditions under which neutral sites could diverge. Indeed, even very low rates of effective gene flow were sufficient to prevent differentiation at neutral sites. However, once strong enough, coupling among selected sites eventually reduced gene flow at neutral sites as well. To explore whether similar transitions might be detectable in empirical data, we used published genome resequencing data from three taxa of Heliconius butterflies. We found that fixation index ( F S T ) outliers and allele-frequency outliers exhibited stronger patterns of within-deme LD than the genomic background, as expected. The statistical characteristics of within-deme LD-likely indicative of the strength of coupling of barrier loci-varied between chromosomes and taxonomic comparisons. Qualitatively, the patterns we observed in the empirical data and in our simulations suggest that selection drives rapid genome-wide transitions to multi-locus coupling, illustrating how divergence and gene flow interact along the speciation continuum.

Admixture, evolution, and variation in reproductive isolation in the Boechera puberula clade

BACKGROUND

Hybridization is very common in plants, and the incorporation of new alleles into existing lineages (i.e. admixture) can blur species boundaries. However, admixture also has the potential to increase standing genetic variation. With new sequencing methods, we can now study admixture and reproductive isolation at a much finer scale than in the past. The genus Boechera is an extraordinary example of admixture, with over 400 hybrid derivates of varying ploidy levels. Yet, few studies have assessed admixture in this genus on a genomic scale.

RESULTS

In this study, we used Genotyping-by-Sequencing (GBS) to clarify the evolution of the Boechera puberula clade, whose six members are scattered across the western United States. We further assessed patterns of admixture and reproductive isolation within the group, including two additional species (B. stricta and B. retrofracta) that are widespread across North America. Based on 14,815 common genetic variants, we found evidence for some cases of hybridization. We find evidence of both recent and more ancient admixture, and that levels of admixture vary across species.

CONCLUSIONS

We present evidence for a monophyletic origin of the B. puberula group, and a split of B. puberula into two subspecies. Further, when inferring reproductive isolation on the basis of presence and absence of admixture, we found that the accumulation of reproductive isolation between species does not seem to occur linearly with time since divergence in this system. We discuss our results in the context of sexuality and asexuality in Boechera.