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 of 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 disappear 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.

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.

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), that is, 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, 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.

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.

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.

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.

Additive genetic effects in interacting species jointly determine the outcome of caterpillar herbivory.

Plant-insect interactions are common and important in basic and applied biology. Trait and genetic variation can affect the outcome and evolution of these interactions, but the relative contributions of plant and insect genetic variation and how these interact remain unclear and are rarely subject to assessment in the same experimental context. Here, we address this knowledge gap using a recent host-range expansion onto alfalfa by the Melissa blue butterfly. Common garden rearing experiments and genomic data show that caterpillar performance depends on plant and insect genetic variation, with insect genetics contributing to performance earlier in development and plant genetics later. Our models of performance based on caterpillar genetics retained predictive power when applied to a second common garden. Much of the plant genetic effect could be explained by heritable variation in plant phytochemicals, especially saponins, peptides, and phosphatidyl cholines, providing a possible mechanistic understanding of variation in the species interaction. We find evidence of polygenic, mostly additive effects within and between species, with consistent effects of plant genotype on growth and development across multiple butterfly species. Our results inform theories of plant-insect coevolution and the evolution of diet breadth in herbivorous insects and other host-specific parasites.

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.