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.

Simulated climate warming causes asymmetric responses in insect life-history timing potentially disrupting a classic ecological speciation system.

Climate change may alter phenology within populations with cascading consequences for community interactions and on-going evolutionary processes. Here, we measured the response to climate warming in two sympatric, recently diverged (~170 years) populations of Rhagoletis pomonella flies specialized on different host fruits (hawthorn and apple) and their parasitoid wasp communities. We tested whether warmer temperatures affect dormancy regulation and its consequences for synchrony across trophic levels and temporal isolation between divergent populations. Under warmer temperatures, both fly populations developed earlier. However, warming significantly increased the proportion of maladaptive pre-winter development in apple, but not hawthorn, flies. Parasitoid phenology was less affected, potentially generating ecological asynchrony. Observed shifts in fly phenology under warming may decrease temporal isolation, potentially limiting on-going divergence. Our findings of complex sensitivity of life-history timing to changing temperatures predict that coming decades may see multifaceted ecological and evolutionary changes in temporal specialist communities.

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.

Genome-wide variation and transcriptional changes in diverse developmental processes underlie the rapid evolution of seasonal adaptation.

Many organisms enter a dormant state in their life cycle to deal with predictable changes in environments over the course of a year. The timing of dormancy is therefore a key seasonal adaptation, and it evolves rapidly with changing environments. We tested the hypothesis that differences in the timing of seasonal activity are driven by differences in the rate of development during diapause in Rhagoletis pomonella, a fly specialized to feed on fruits of seasonally limited host plants. Transcriptomes from the central nervous system across a time series during diapause show consistent and progressive changes in transcripts participating in diverse developmental processes, despite a lack of gross morphological change. Moreover, population genomic analyses suggested that many genes of small effect enriched in developmental functional categories underlie variation in dormancy timing and overlap with gene sets associated with development rate in Drosophila melanogaster Our transcriptional data also suggested that a recent evolutionary shift from a seasonally late to a seasonally early host plant drove more rapid development during diapause in the early fly population. Moreover, genetic variants that diverged during the evolutionary shift were also enriched in putative cis regulatory regions of genes differentially expressed during diapause development. Overall, our data suggest polygenic variation in the rate of developmental progression during diapause contributes to the evolution of seasonality in R. pomonella We further discuss patterns that suggest hourglass-like developmental divergence early and late in diapause development and an important role for hub genes in the evolution of transcriptional divergence.

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.

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.

Divergent diapause life history timing drives both allochronic speciation and reticulate hybridization in an adaptive radiation of Rhagoletis flies.

Divergent adaptation to new ecological opportunities can be an important factor initiating speciation. However, as niches are filled during adaptive radiations, trait divergence driving reproductive isolation between sister taxa may also result in trait convergence with more distantly related taxa, increasing the potential for reticulated gene flow across the radiation. Here, we demonstrate such a scenario in a recent adaptive radiation of Rhagoletis fruit flies, specialized on different host plants. Throughout this radiation, shifts to novel hosts are associated with changes in diapause life history timing, which act as "magic traits" generating allochronic reproductive isolation and facilitating speciation-with-gene-flow. Evidence from laboratory rearing experiments measuring adult emergence timing and genome-wide DNA-sequencing surveys supported allochronic speciation between summer-fruiting Vaccinium spp.-infesting Rhagoletis mendax and its hypothesized and undescribed sister taxon infesting autumn-fruiting sparkleberries. The sparkleberry fly and R. mendax were shown to be genetically discrete sister taxa, exhibiting no detectable gene flow and allochronically isolated by a 2-month average difference in emergence time corresponding to host availability. At sympatric sites across the southern USA, the later fruiting phenology of sparkleberries overlaps with that of flowering dogwood, the host of another more distantly related and undescribed Rhagoletis taxon. Laboratory emergence data confirmed broadly overlapping life history timing and genomic evidence supported on-going gene flow between sparkleberry and flowering dogwood flies. Thus, divergent phenological adaptation can drive the initiation of reproductive isolation, while also enhancing genetic exchange across broader adaptive radiations, potentially serving as a source of novel genotypic variation and accentuating further diversification.

Genomically correlated trait combinations and antagonistic selection contributing to counterintuitive genetic patterns of adaptive diapause divergence in Rhagoletis flies.

Adaptation to novel environments can result in unanticipated genomic responses to selection. Here, we illustrate how multifarious, correlational selection helps explain a counterintuitive pattern of genetic divergence between the recently derived apple- and ancestral hawthorn-infesting host races of Rhagoletis pomonella (Diptera: Tephritidae). The apple host race terminates diapause and emerges as adults earlier in the season than the hawthorn host race, to coincide with the earlier fruiting phenology of their apple hosts. However, alleles at many loci associated with later emergence paradoxically occur at higher frequencies in sympatric populations of the apple compared to the hawthorn race. We present genomic evidence that historical selection over geographically varying environmental gradients across North America generated genetic correlations between two life history traits, diapause intensity and diapause termination, in the hawthorn host race. Moreover, the loci associated with these life history traits are concentrated in genomic regions in high linkage disequilibrium (LD). These genetic correlations are antagonistic to contemporary selection on local apple host race populations that favours increased initial diapause depth and earlier, not later, diapause termination. Thus, the paradox of apple flies appears due, in part, to pleiotropy or linkage of alleles associated with later adult emergence and increased initial diapause intensity, the latter trait strongly selected for by the earlier phenology of apples. Our results demonstrate how understanding of multivariate trait combinations and the correlative nature of selective forces acting on them can improve predictions concerning adaptive evolution and help explain seemingly counterintuitive patterns of genetic diversity in nature.

A reversal in sensory processing accompanies ongoing ecological divergence and speciation in Rhagoletis pomonella.

Changes in behaviour often drive rapid adaptive evolution and speciation. However, the mechanistic basis for behavioural shifts is largely unknown. The tephritid fruit fly Rhagoletis pomonella is an example of ecological specialization and speciation in action via a recent host plant shift from hawthorn to apple. These flies primarily use specific odours to locate fruit, and because they mate only on or near host fruit, changes in odour preference for apples versus hawthorns translate directly to prezygotic reproductive isolation, initiating speciation. Using a variety of techniques, we found a reversal between apple and hawthorn flies in the sensory processing of key odours associated with host fruit preference at the first olfactory synapse, linking changes in the antennal lobe of the brain with ongoing ecological divergence. Indeed, changes to specific neural pathways of any sensory modality may be a broad mechanism for changes in animal behaviour, catalysing the genesis of new biodiversity.

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.

Comparative genome sequencing reveals insights into the dynamics of Wolbachia in native and invasive cherry fruit flies.

Wolbachia is a maternally inherited obligate endosymbiont that can induce a wide spectrum of effects in its host, ranging from mutualism to reproductive parasitism. At the genomic level, recombination within and between strains, transposable elements, and horizontal transfer of strains between host species make Wolbachia an evolutionarily dynamic bacterial system. The invasive cherry fruit fly Rhagoletis cingulata arrived in Europe from North America ~40 years ago, where it now co-occurs with the native cherry pest R. cerasi. This shared distribution has been proposed to have led to the horizontal transfer of different Wolbachia strains between the two species. To better understand transmission dynamics, we performed a comparative genome study of the strain wCin2 in its native United States and invasive European populations of R. cingulata with wCer2 in European R. cerasi. Previous multilocus sequence genotyping (MLST) of six genes implied that the source of wCer2 in R. cerasi was wCin2 from R. cingulata. However, we report genomic evidence discounting the recent horizontal transfer hypothesis for the origin of wCer2. Despite near identical sequences for the MLST markers, substantial sequence differences for other loci were found between wCer2 and wCin2, as well as structural rearrangements, and differences in prophage, repetitive element, gene content, and cytoplasmic incompatibility inducing genes. Our study highlights the need for whole-genome sequencing rather than relying on MLST markers for resolving Wolbachia strains and assessing their evolutionary dynamics.

A rapidly evolved shift in life-history timing during ecological speciation is driven by the transition between developmental phases.

For insect species in temperate environments, seasonal timing is often governed by the regulation of diapause, a complex developmental programme that allows insects to weather unfavourable conditions and synchronize their life cycles with available resources. Diapause development consists of a series of distinct phases including initiation, maintenance, termination and post-diapause development. The evolution of insect seasonal timing depends in part on how these phases of diapause development and post-diapause development interact to affect variation in phenology. Here, we dissect the physiological basis of a recently evolved phenological shift in Rhagoletis pomonella (Diptera: Tephritidae), a model system for ecological divergence. A recently derived population of R. pomonella shifted from specializing on native hawthorn fruit to earlier fruiting introduced apples, resulting in a 3-4 week shift in adult emergence timing. We tracked metabolic rates of individual flies across post-winter development to test which phases of development may act either independently or in combination to contribute to this recently evolved divergence in timing. Apple and hawthorn flies differed in a number of facets of their post-winter developmental trajectories. However, divergent adaptation in adult emergence phenology in these flies was due almost entirely to the end of the pupal diapause maintenance phase, with post-diapause development having a very small effect. The relatively simple underpinnings of variation in adult emergence phenology suggest that further adaptation to seasonal change in these flies for this trait might be largely due to the timing of diapause termination unhindered by strong covariance among different components of post-diapause development.

Genomics and the origin of species.

Speciation is a fundamental evolutionary process, the knowledge of which is crucial for understanding the origins of biodiversity. Genomic approaches are an increasingly important aspect of this research field. We review current understanding of genome-wide effects of accumulating reproductive isolation and of genomic properties that influence the process of speciation. Building on this work, we identify emergent trends and gaps in our understanding, propose new approaches to more fully integrate genomics into speciation research, translate speciation theory into hypotheses that are testable using genomic tools and provide an integrative definition of the field of speciation genomics.

Comparing Adaptive Radiations Across Space, Time, and Taxa.

Adaptive radiation plays a fundamental role in our understanding of the evolutionary process. However, the concept has provoked strong and differing opinions concerning its definition and nature among researchers studying a wide diversity of systems. Here, we take a broad view of what constitutes an adaptive radiation, and seek to find commonalities among disparate examples, ranging from plants to invertebrate and vertebrate animals, and remote islands to lakes and continents, to better understand processes shared across adaptive radiations. We surveyed many groups to evaluate factors considered important in a large variety of species radiations. In each of these studies, ecological opportunity of some form is identified as a prerequisite for adaptive radiation. However, evolvability, which can be enhanced by hybridization between distantly related species, may play a role in seeding entire radiations. Within radiations, the processes that lead to speciation depend largely on (1) whether the primary drivers of ecological shifts are (a) external to the membership of the radiation itself (mostly divergent or disruptive ecological selection) or (b) due to competition within the radiation membership (interactions among members) subsequent to reproductive isolation in similar environments, and (2) the extent and timing of admixture. These differences translate into different patterns of species accumulation and subsequent patterns of diversity across an adaptive radiation. Adaptive radiations occur in an extraordinary diversity of different ways, and continue to provide rich data for a better understanding of the diversification of life.

Host plant-related genomic differentiation in the European cherry fruit fly, Rhagoletis cerasi.

Elucidating the mechanisms and conditions facilitating the formation of biodiversity are central topics in evolutionary biology. A growing number of studies imply that divergent ecological selection may often play a critical role in speciation by counteracting the homogenising effects of gene flow. Several examples involve phytophagous insects, where divergent selection pressures associated with host plant shifts may generate reproductive isolation, promoting speciation. Here, we use ddRADseq to assess the population structure and to test for host-related genomic differentiation in the European cherry fruit fly, Rhagoletis cerasi (L., 1758) (Diptera: Tephritidae). This tephritid is distributed throughout Europe and western Asia, and has adapted to two different genera of host plants, Prunus spp. (cherries) and Lonicera spp. (honeysuckle). Our data imply that geographic distance and geomorphic barriers serve as the primary factors shaping genetic population structure across the species range. Locally, however, flies genetically cluster according to host plant, with consistent allele frequency differences displayed by a subset of loci between Prunus and Lonicera flies across four sites surveyed in Germany and Norway. These 17 loci display significantly higher FST values between host plants than others. They also showed high levels of linkage disequilibrium within and between Prunus and Lonicera flies, supporting host-related selection and reduced gene flow. Our findings support the existence of sympatric host races in R. cerasi embedded within broader patterns of geographic variation in the fly, similar to the related apple maggot, Rhagoletis pomonella, in North America.

Identification of a New Blend of Host Fruit Volatiles from Red Downy Hawthorn, Crataegus mollis, Attractive to Rhagoletis pomonella Flies from the Northeastern United States.

A new blend of volatiles was identified for the fruit of downy red hawthorn, Crataegus mollis, that is attractive to Rhagoletis pomonella flies infesting this host in the northeastern USA. The new blend was as attractive as the previously identified mixture but is more complex in the number of odorants (six in the old versus ten in the new) and differs significantly in the ratio of three volatiles, 3-methylbutan-1-ol, butyl hexanoate, and dihydro-ß-ionone, that are common to both blends and exerted agonist or antagonist effects on behavior in a flight tunnel assay. However, behavioral results with the old and new northern hawthorn blends, as well as modified blends with substituted ratios of 3-methylbutan-1-ol, butyl hexanoate, dihydro-ß-ionone, indicated that the 'agonist' or 'antagonist' effects of these volatiles depended on the ratio, or balance of compounds within the blend. In addition, the new blend contains a number of esters identified from the headspace of domesticated apple, Malus domestica, that are attractive to apple-origin R. pomonella, and present in the five other blends from southern hawthorns, including the southern C. mollis var. texana blend, but are not part of the previously identified blend from northern C. mollis fruit. This finding supports the hypothesis that in addition to providing specificity to the odor blends of the northern and southern hawthorn populations, the presence of the significant amounts of ester compounds in the new northern hawthorn blend might have provided a source of standing variation that could help explain the shift in host preference by C. mollis-infesting flies to introduced apple in the mid-1800's.

Sequential divergence and the multiplicative origin of community diversity.

Phenotypic and genetic variation in one species can influence the composition of interacting organisms within communities and across ecosystems. As a result, the divergence of one species may not be an isolated process, as the origin of one taxon could create new niche opportunities for other species to exploit, leading to the genesis of many new taxa in a process termed "sequential divergence." Here, we test for such a multiplicative effect of sequential divergence in a community of host-specific parasitoid wasps, Diachasma alloeum, Utetes canaliculatus, and Diachasmimorpha mellea (Hymenoptera: Braconidae), that attack Rhagoletis pomonella fruit flies (Diptera: Tephritidae). Flies in the R. pomonella species complex radiated by sympatrically shifting and ecologically adapting to new host plants, the most recent example being the apple-infesting host race of R. pomonella formed via a host plant shift from hawthorn-infesting flies within the last 160 y. Using population genetics, field-based behavioral observations, host fruit odor discrimination assays, and analyses of life history timing, we show that the same host-related ecological selection pressures that differentially adapt and reproductively isolate Rhagoletis to their respective host plants (host-associated differences in the timing of adult eclosion, host fruit odor preference and avoidance behaviors, and mating site fidelity) cascade through the ecosystem and induce host-associated genetic divergence for each of the three members of the parasitoid community. Thus, divergent selection at lower trophic levels can potentially multiplicatively and rapidly amplify biodiversity at higher levels on an ecological time scale, which may sequentially contribute to the rich diversity of life.

Multilocus approaches for the measurement of selection on correlated genetic loci.

The study of ecological speciation is inherently linked to the study of selection. Methods for estimating phenotypic selection within a generation based on associations between trait values and fitness (e.g. survival) of individuals are established. These methods attempt to disentangle selection acting directly on a trait from indirect selection caused by correlations with other traits via multivariate statistical approaches (i.e. inference of selection gradients). The estimation of selection on genotypic or genomic variation could also benefit from disentangling direct and indirect selection on genetic loci. However, achieving this goal is difficult with genomic data because the number of potentially correlated genetic loci (p) is very large relative to the number of individuals sampled (n). In other words, the number of model parameters exceeds the number of observations (p â‰« n). We present simulations examining the utility of whole-genome regression approaches (i.e. Bayesian sparse linear mixed models) for quantifying direct selection in cases where p â‰« n. Such models have been used for genome-wide association mapping and are common in artificial breeding. Our results show they hold promise for studies of natural selection in the wild and thus of ecological speciation. But we also demonstrate important limitations to the approach and discuss study designs required for more robust inferences.

A test of genomic modularity among life-history adaptations promoting speciation with gene flow.

Speciation with gene flow may require adaptive divergence of multiple traits to generate strong ecologically based reproductive isolation. Extensive negative pleiotropy or physical linkage of genes in the wrong phase affecting these diverging traits may therefore hinder speciation, while genetic independence or "modularity" among phenotypic traits may reduce constraints and facilitate divergence. Here, we test whether the genetics underlying two components of diapause life history, initial diapause intensity and diapause termination timing, constrain differentiation between sympatric hawthorn and apple-infesting host races of the fly Rhagoletis pomonella through analysis of 10,256 SNPs measured via genotyping-by-sequencing (GBS). Loci genetically associated with diapause termination timing were mainly observed for SNPs mapping to chromosomes 1-3 in the genome, most notably for SNPs displaying higher levels of linkage disequilibrium (LD), likely due to inversions. In contrast, selection on initial diapause intensity affected loci on all five major chromosomes of the genome, specifically those showing low levels of LD. This lack of overlap in genetically associated loci suggests that the two diapause phenotypes are largely modular. On chromosome 2, however, intermediate level LD loci and a subgroup of high LD loci displayed significant negative relationships between initial diapause intensity and diapause termination time. These gene regions on chromosome 2 therefore affected both traits, while most regions were largely independent. Moreover, loci associated with both measured traits also tended to exhibit highly divergent allele frequencies between the host races. Thus, the presence of nonoverlapping genetic modules likely facilitates simultaneous, adaptive divergence for the measured life-history components.