The impact of geographic isolation and host shifts on population divergence of the rare cicada Subpsaltria yangi

The contributions of divergent selection and spatial isolation to population divergence are among the main focuses of evolutionary biology. Here we employed integrated methods to explore genomic divergence, demographic history and calling-song differentiation in the cicada Subpsaltria yangi, and compared the genotype and calling-song phenotype of different populations occurring in distinct habitats. Our results indicate that this species comprises four main lineages with unique sets of haplotypes and calling-song structure, which are distinctly associated with geographic isolation and habitats. The populations occurring on the Loess Plateau underwent substantial expansion at ∼0.130-0.115 Ma during the Last Interglacial. Geographic distance and host shift between pairs of populations predict genomic divergence, with geographic distance and acoustical signal together explaining > 60% of the divergence among populations. Differences in calling songs could reflect adaptation of populations to novel environments with different host plants, habitats and predators, which may have resulted from neutral divergence at the molecular level followed by natural selection. Geomorphic barriers and climate oscillations associated with Pleistocene glaciation may have been primary factors in shaping the population genetic structure of this species. Ultimately this may couple with a host shift in leading toward allopatric speciation in S. yangi, i.e., isolation by distance. Our findings improve understanding of divergence in allopatry of herbivorous insects, and may inform future studies on the molecular mechanisms underlying the association between genetic/phenotypic changes and adaptation of insects to novel niches and host plants.

The role of sexual isolation during rapid ecological divergence: Evidence for a new dimension of isolation in Rhagoletis pomonella

The pace of divergence and likelihood of speciation often depends on how and when different types of reproductive barriers evolve. Questions remain about how reproductive isolation evolves after initial divergence. We tested for the presence of sexual isolation (reduced mating between populations due to divergent mating preferences and traits) in Rhagoletis pomonella flies, a model system for incipient ecological speciation. We measured the strength of sexual isolation between two very recently diverged (~170 generations) sympatric populations, adapted to different host fruits (hawthorn and apple). We found that flies from both populations were more likely to mate within than between populations. Thus, sexual isolation may play an important role in reducing gene flow allowed by early-acting ecological barriers. We also tested how warmer temperatures predicted under climate change could alter sexual isolation and found that sexual isolation was markedly asymmetric under warmer temperatures - apple males and hawthorn females mated randomly while apple females and hawthorn males mated more within populations than between. Our findings provide a window into the early speciation process and the role of sexual isolation after initial ecological divergence, in addition to examining how environmental conditions could shape the likelihood of further divergence.

Cold tolerance and diapause within and across trophic levels: Endoparasitic wasps and their fly host have similar phenotypes

Low temperatures associated with winter can limit the survival of organisms, especially ectotherms whose body temperature is similar to their environment. However, there is a gap in understanding how overwintering may vary among groups of species that interact closely, such as multiple parasitoid species that attack the same host insect. Here, we investigate cold tolerance and diapause phenotypes in three endoparasitoid wasps of the apple maggot fly Rhagoletis pomonella (Diptera: Tephritidae): Utetes canaliculatus, Diachasma alloeum, and Diachasmimorpha mellea (Hymenoptera: Braconidae). Using a combination of respirometry and eclosion tracking, we found that all three wasp species exhibited the same three diapause duration phenotypes as the fly host. Weak (short duration) diapause was rare, with <5 % of all three wasp species prematurely terminating diapause at 21 °C. Most D.mellea (93 %) entered a more intense (longer duration) diapause that did not terminate within 100 d at this warm temperature. The majority of U.canaliculatus (92 %) and D. alloeum (72 %) averted diapause (non-diapause) at 21 °C. There was limited interspecific variation in acute cold tolerance among the three wasp species: wasps and flies had similarly high survival (>87 %) following exposure to extreme low temperatures (-20 °C) as long as their body fluids did not freeze. The three wasp species also displayed little interspecific variation in survival following prolonged exposure to mild chilling of 8 or more weeks at 4 °C. Our study thus documents a remarkable conservation of cold tolerance and diapause phenotypes within and across trophic levels.

Variation in the genomic basis of parallel phenotypic and ecological divergence in benthic and pelagic morphs of Icelandic Arctic charr (Salvelinus alpinus)

Sympatric adaptive phenotypic divergence should be underlain by genomic differentiation between sub-populations. When divergence drives similar patterns of phenotypic and ecological variation within species we expect evolution to draw on common allelic variation. We investigated divergence histories and genomic signatures of adaptive divergence between benthic and pelagic morphs of Icelandic Arctic charr. Divergence histories for each of four populations were reconstructed using coalescent modelling and 14,187 single nucleotide polymorphisms. Sympatric divergence with continuous gene flow was supported in two populations while allopatric divergence with secondary contact was supported in one population; we could not differentiate between demographic models in the fourth population. We detected parallel patterns of phenotypic divergence along benthic-pelagic evolutionary trajectories among populations. Patterns of genomic differentiation between benthic and pelagic morphs were characterized by outlier loci in many narrow peaks of differentiation throughout the genome, which may reflect the eroding effects of gene flow on nearby neutral loci. We then used genome-wide association analyses to relate both phenotypic (body shape and size) and ecological (carbon and nitrogen stable isotopes) variation to patterns of genomic differentiation. Many peaks of genomic differentiation were associated with phenotypic and ecological variation in the three highly divergent populations, suggesting a genomic basis for adaptive divergence. We detected little evidence for a parallel genomic basis of differentiation as most regions and outlier loci were not shared among populations. Our results show that adaptive divergence can have varied genomic consequences in populations with relatively recent common origins, similar divergence histories, and parallel phenotypic divergence.

Pool-GWAS on reproductive dormancy in Drosophila simulans suggests a polygenic architecture

The genetic basis of adaptation to different environments has been of long-standing interest to evolutionary biologists. Dormancy is a well-studied adaptation to facilitate overwintering. In Drosophila melanogaster, a moderate number of genes with large effects have been described, which suggests a simple genetic basis of dormancy. On the other hand, genome-wide scans for dormancy suggest a polygenic architecture in insects. In D. melanogaster, the analysis of the genetic architecture of dormancy is complicated by the presence of cosmopolitan inversions. Here, we performed a genome-wide scan to characterize the genetic basis of this ecologically extremely important trait in the sibling species of D. melanogaster, D. simulans that lacks cosmopolitan inversions. We performed Pool-GWAS in a South African D. simulans population for dormancy incidence at 2 temperature regimes (10 and 12°C, LD 10:14). We identified several genes with SNPs that showed a significant association with dormancy (P-value < 1e-13), but the overall modest response suggests that dormancy is a polygenic trait with many loci of small effect. Our results shed light on controversies on reproductive dormancy in Drosophila and have important implications for the characterization of the genetic basis of this trait.

The Build-Up of Population Genetic Divergence along the Speciation Continuum during a Recent Adaptive Radiation of Rhagoletis Flies

New species form through the evolution of genetic barriers to gene flow between previously interbreeding populations. The understanding of how speciation proceeds is hampered by our inability to follow cases of incipient speciation through time. Comparative approaches examining different diverging taxa may offer limited inferences, unless they fulfill criteria that make the comparisons relevant. Here, we test for those criteria in a recent adaptive radiation of the Rhagoletis pomonella species group (RPSG) hypothesized to have diverged in sympatry via adaptation to different host fruits. We use a large-scale population genetic survey of 1568 flies across 33 populations to: (1) detect on-going hybridization, (2) determine whether the RPSG is derived from the same proximate ancestor, and (3) examine patterns of clustering and differentiation among sympatric populations. We find that divergence of each in-group RPSG taxon is occurring under current gene flow, that the derived members are nested within the large pool of genetic variation present in hawthorn-infesting populations of R. pomonella, and that sympatric population pairs differ markedly in their degree of genotypic clustering and differentiation across loci. We conclude that the RPSG provides a particularly robust opportunity to make direct comparisons to test hypotheses about how ecological speciation proceeds despite on-going gene flow.

Genomic Traces of the Fruit Fly Anastrepha obliqua Associated with Its Polyphagous Nature

Anastrepha obliqua (Macquart) (Diptera: Tephritidae) is an important pest in the neotropical region. It is considered a polyphagous insect, meaning it infests plants of different taxonomic families and readily colonizes new host plants. The change to new hosts can lead to diversification and the formation of host races. Previous studies investigating the effect of host plants on population structure and selection in Anastrepha obliqua have focused on the use of data from the mitochondrial DNA sequence and microsatellite markers of nuclear DNA, and there are no analyses at the genomic level. To better understand this issue, we used a pooled restriction site-associated DNA sequencing (pooled RAD-seq) approach to assess genomic differentiation and population structure across sympatric populations of Anastrepha obliqua that infest three host plants-Spondias purpurea (red mombin), Mangifera indica (mango) of the family Anacardiaceae and Averrhoa carambola (carambola) of the family Oxalidaceae-in sympatric populations of the species Anastrepha obliqua of Inter-Andean Valley of the Cauca River in southwestern Colombia. Our results show genomic differentiation of populations from carambola compared to mango and red mombin populations, but the genetic structure was mainly established by geography rather than by the host plant. On the other hand, we identified 54 SNPs in 23 sequences significantly associated with the use of the host plant. Of these 23 sequences, we identified 17 candidate genes and nine protein families, of which four protein families are involved in the nutrition of these flies. Future studies should investigate the adaptive processes undergone by phytophagous insects in the Neotropics, using fruit flies as a model and state-of-the-art molecular tools.

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.

Contrast in Post-Chill Eclosion Time Strategies Between Two Specialist Braconid Wasps (Hymenoptera: Braconidae) Attacking Rhagoletis Flies (Diptera: Tephritidae) in Western North America

Parasitoids comprise a speciose insect group, displaying a wide array of life history strategies. In the Pacific Northwest of the United States, the tephritid fruit flies Rhagoletis tabellaria (Fitch) and Rhagoletis indifferens Curran infest red osier dogwood, Cornus sericea L. (Cornaceae), and bitter cherry, Prunus emarginata (Douglas ex Hooker) Eaton (Rosaceae), respectively. The flies are parasitized by different braconid wasps at different life stages; Utetes tabellariae (Fischer) oviposits into R. tabellaria eggs, whereas Diachasma muliebre (Muesebeck) oviposits into R. indifferens larvae feeding in cherries. Because Rhagoletis only have one major generation a year and the wasps attack temporally distinct fly life stages, we predicted that eclosion times of U. tabellariae should more closely follow that of its host than the larval-attacking D. muliebre. As predicted, U. tabellariae eclosed on average 6.0-12.5 d later than R. tabellaria, whereas D. muliebre eclosed on average 32.1 d after R. indifferens. Unexpectedly, however, longer chill duration differentially affected the systems; longer overwinters minimally influenced eclosion times of R. tabellaria and U. tabellariae but caused earlier eclosion of both R. indifferens and D. muliebre. Results imply that in temperate regions, diapause timing in braconid wasps evolves in response to both host life stage attacked and fly eclosion characteristics, possibly reflecting differential effects of winter on host plant fruiting phenology. Differences in phenological sensitivity of the lower host plant trophic level to variation in environmental conditions may have cascading effects, sequentially and differentially affecting eclosion times in higher frugivore (fly) and parasitoid (wasp) trophic levels.

Reinventing the wheel? Reassessing the roles of gene flow, sorting and convergence in repeated evolution

Biologists have long been intrigued by apparently predictable and repetitive evolutionary trajectories inferred across a variety of lineages and systems. In recent years, high-throughput sequencing analyses have started to transform our understanding of such repetitive shifts. While researchers have traditionally categorized such shifts as either "convergent" or "parallel," based on relatedness of the lineages involved, emerging genomic insights provide an opportunity to better describe the actual evolutionary mechanisms at play. A synthesis of recent genomic analyses confirms that convergence is the predominant driver of repetitive evolution among species, whereas repeated sorting of standing variation is the major driver of repeated shifts within species. However, emerging data reveal numerous notable exceptions to these expectations, with recent examples of de novo mutations underpinning convergent shifts among even very closely related lineages, while repetitive sorting processes have occurred among even deeply divergent taxa, sometimes via introgression. A number of very recent analyses have found evidence for both processes occurring on different scales within taxa. We suggest that the relative importance of convergent versus sorting processes depends on the interplay between gene flow among populations, and phylogenetic relatedness of the lineages involved.

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.

A Linkage-Based Genome Assembly for the Mosquito Aedes albopictus and Identification of Chromosomal Regions Affecting Diapause

Simple Summary

A genome sequence can provide the basis for understanding a wide range of biological properties of an organism. For vector and pest species, this knowledge can be used to develop novel control strategies based on genome modifications that disrupt traits related to ecological adaptation and disease transmission. Here, we use a genetic mapping experiment to produce an improved version of a previous genome sequence for the invasive vector mosquito, Aedes albopictus. We then use the improved genome sequence to identify regions of the genome that contain candidate genes that may affect the ability of this mosquito to undergo overwintering dormancy, a crucial ecological adaptation.

Abstract

The Asian tiger mosquito, Aedes albopictus, is an invasive vector mosquito of substantial public health concern. The large genome size (~1.19–1.28 Gb by cytofluorometric estimates), comprised of ~68% repetitive DNA sequences, has made it difficult to produce a high-quality genome assembly for this species. We constructed a high-density linkage map for Ae. albopictus based on 111,328 informative SNPs obtained by RNAseq. We then performed a linkage-map anchored reassembly of AalbF2, the genome assembly produced by Palatini et al. (2020). Our reassembled genome sequence, AalbF3, represents several improvements relative to AalbF2. First, the size of the AalbF3 assembly is 1.45 Gb, almost half the size of AalbF2. Furthermore, relative to AalbF2, AalbF3 contains a higher proportion of complete and single-copy BUSCO genes (84.3%) and a higher proportion of aligned RNAseq reads that map concordantly to a single location of the genome (46%). We demonstrate the utility of AalbF3 by using it as a reference for a bulk-segregant-based comparative genomics analysis that identifies chromosomal regions with clusters of candidate SNPs putatively associated with photoperiodic diapause, a crucial ecological adaptation underpinning the rapid range expansion and climatic adaptation of A. albopictus.

Costs of averting or prematurely terminating diapause associated with slow decline of metabolic rates at low temperature

Diapause, a form of insect dormancy, generally facilitates overwintering by increasing cold tolerance and decreasing energy drain at high temperatures via metabolic rate suppression. Averting or terminating diapause prior to winter is generally assumed to be a lethal phenotype. However, low temperature acclimation can also increase cold tolerance and decrease metabolic rates. Here, we tested the hypothesis that non- and post-diapause individuals in a cold-induced quiescence can achieve a diapause-like phenotype, compensating for the potential costs of averting diapause. We tested this in the apple maggot fly Rhagoletis pomonella, which typically overwinters in the soil as a diapause pupa, but can avert diapause (non-diapause) or terminate diapause early ('weak diapause') when reared at warm temperatures. Metabolic rates were initially higher in non- and post-diapause than diapause pupae at high (25 °C) and low (4 °C) temperatures, but quiescent non- and post-diapause pupae achieved diapause-like metabolic rates slowly over time when incubated at 4 °C for several weeks. We found that diapause and quiescent pupae were freeze-avoidant and had similar tolerance of extreme low temperatures (cooling to c. -18 °C) following 8 weeks acclimation at 4 °C. Despite high tolerance of subzero temperatures, quiescent pupae did not survive well when chilled for prolonged periods (8 weeks or more) at 4 °C. We conclude that cold acclimation can only partially compensate for costs associated with aversion or premature termination of diapause, and that energy drain at low (not just high) temperatures likely contributes to chilling mortality in quiescent insects.

Timing strains of the marine insect Clunio marinus diverged and persist with gene flow

Genetic divergence of populations in the presence of gene flow is a central theme in speciation research. Theory predicts that divergence can happen with full range overlap - in sympatry - driven by ecological factors, but there are few empirical examples of how ecologically divergent selection can overcome gene flow and lead to reproductive isolation. In the marine midge Clunio marinus (Diptera: Chironomidae) reproduction is ecologically restricted to the time of the lowest tides, which is ensured through accurate control of development and adult emergence by circalunar and circadian clocks. As tidal regimes differ along the coastline, locally adapted timing strains of C. marinus are found in different sites across Europe. At the same time, ecologically suitable low tides occur at both full and new moon and twice a day, providing C. marinus with four nonoverlapping temporal niches at every geographic location. Along the coast of Brittany, which is characterized by a steep gradient in timing of the tides, we found an unusually large number of differentially adapted timing strains, and the first known instances of sympatric C. marinus strains occupying divergent temporal niches. Analysis of mitochondrial genotypes suggests that these timing strains originated from a single recent colonization event. Nuclear genotypes show strong gene flow, sympatric timing strains being the least differentiated. Even when sympatric strains exist in nonoverlapping temporal niches, timing adaptations do not result in genome-wide genetic divergence, suggesting timing adaptations are maintained by permanent ecological selection. This constitutes a model case for incipient ecological divergence with gene flow.

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.

Divergence, gene flow, and speciation in eight lineages of trans-Beringian birds

Determining how genetic diversity is structured between populations that span the divergence continuum from populations to biological species is key to understanding the generation and maintenance of biodiversity. We investigated genetic divergence and gene flow in eight lineages of birds with a trans-Beringian distribution, where Asian and North American populations have likely been split and reunited through multiple Pleistocene glacial cycles. Our study transects the speciation process, including eight pairwise comparisons in three orders (ducks, shorebirds and passerines) at population, subspecies and species levels. Using ultraconserved elements (UCEs), we found that these lineages represent conditions from slightly differentiated populations to full biological species. Although allopatric speciation is considered the predominant mode of divergence in birds, all of our best divergence models included gene flow, supporting speciation with gene flow as the predominant mode in Beringia. In our eight lineages, three were best described by a split-migration model (divergence with gene flow), three best fit a secondary contact scenario (isolation followed by gene flow), and two showed support for both models. The lineages were not evenly distributed across a divergence space defined by gene flow (M) and differentiation (F_ST ), instead forming two discontinuous groups: one with relatively shallow divergence, no fixed single nucleotide polymorphisms (SNPs), and high rates of gene flow between populations; and the second with relatively deeply divergent lineages, multiple fixed SNPs, and low gene flow. Our results highlight the important role that gene flow plays in avian divergence in Beringia.

Can the genomics of ecological speciation be predicted across the divergence continuum from host races to species? A case study in Rhagoletis

Studies assessing the predictability of evolution typically focus on short-term adaptation within populations or the repeatability of change among lineages. A missing consideration in speciation research is to determine whether natural selection predictably transforms standing genetic variation within populations into differences between species. Here, we test whether and how host-related selection on diapause timing associates with genome-wide differentiation during ecological speciation by comparing ancestral hawthorn and newly formed apple-infesting host races of Rhagoletis pomonella to their sibling species Rhagoletis mendax that attacks blueberries. The associations of 57 857 single nucleotide polymorphisms in a diapause genome-wide-association study (GWAS) on the hawthorn race strongly predicted the direction and magnitude of genomic divergence among the three fly populations at a field site in Fennville, MI, USA. The apple race and R. mendax show parallel changes in the frequencies of putative inversions on three chromosomes associated with the earlier fruiting times of apples and blueberries compared to hawthorns. A diapause GWAS on R. mendax revealed compensatory changes throughout the genome accounting for the earlier eclosion of blueberry, but not apple flies. Thus, a degree of predictability, although not complete, exists in the genomics of diapause across the ecological speciation continuum in Rhagoletis. The generality of this result is placed in the context of other similar systems. This article is part of the theme issue 'Towards the completion of speciation: the evolution of reproductive isolation beyond the first barriers'.

Modular chromosome rearrangements reveal parallel and nonparallel adaptation in a marine fish

Genomic architecture and standing variation can play a key role in ecological adaptation and contribute to the predictability of evolution. In Atlantic cod (Gadus morhua), four large chromosomal rearrangements have been associated with ecological gradients and migratory behavior in regional analyses. However, the degree of parallelism, the extent of independent inheritance, and functional distinctiveness of these rearrangements remain poorly understood. Here, we use a 12K single nucleotide polymorphism (SNP) array to demonstrate extensive individual variation in rearrangement genotype within populations across the species range, suggesting that local adaptation to fine-scale ecological variation is enabled by rearrangements with independent inheritance. Our results demonstrate significant association of rearrangements with migration phenotype and environmental gradients across the species range. Individual rearrangements exhibit functional modularity, but also contain loci showing multiple environmental associations. Clustering in genetic distance trees and reduced differentiation within rearrangements across the species range are consistent with shared variation as a source of contemporary adaptive diversity in Atlantic cod. Conversely, we also find that haplotypes in the LG12 and LG1 rearranged region have diverged across the Atlantic, despite consistent environmental associations. Exchange of these structurally variable genomic regions, as well as local selective pressures, has likely facilitated individual diversity within Atlantic cod stocks. Our results highlight the importance of genomic architecture and standing variation in enabling fine-scale adaptation in marine species.

Evolutionary and functional genetics of insect diapause: a call for greater integration

Diapause in response to seasonality is an important model for rapid evolutionary adaptation that is highly genetically variable, and experiences strong natural selection. Forward genetic methods using various genomic and transcriptomic approaches have begun to characterize the genetic architecture and candidate genes underlying diapause evolution. Largely in parallel, reverse genetic studies have identified functional roles for candidate genes that may or may not be genetically variable. We illustrate the disconnect between the evolutionary and physiological literature using a suite of studies of the role of the circadian clock in diapause regulation. These extensive studies in two different disciplines provide excellent opportunities for integration, which should facilitate rapid progress in understanding both the regulation and evolution of diapause.

The genetic architecture underlying diapause termination in a planktonic crustacean

Diapause is a feature of the life cycle of many invertebrates by which unfavourable environmental conditions can be outlived. The seasonal timing of diapause allows organisms to adapt to seasonal changes in habitat suitability and thus is key to their fitness. In the planktonic crustacean Daphnia, various cues can induce the production of diapause stages that are resistant to heat, drought or freezing and contain one to two embryos in developmental arrest. Daphnia is a keystone species of many freshwater ecosystems, where it acts as the main link between phytoplankton and higher trophic levels. The correct seasonal timing of diapause termination is essential to maintain trophic interactions and is achieved via a genetically based interpretation of environmental cues like photoperiod and temperature. Field monitoring and modelling studies raised concerns on whether populations can advance their seasonal release from diapause to advances in spring phenology under global change, or if a failure to adapt will cause trophic mismatches negatively affecting ecosystem functioning. Our capacity to understand and predict the evolution of diapause timing requires information about the genetic architecture underlying this trait. In this study, we identified eight quantitative trait loci (QTLs) and four epistatic interactions that together explained 66.5% of the variation in diapause termination in Daphnia magna using QTL mapping. Our results suggest that the most significant QTL is modulating diapause termination dependent on photoperiod and is involved in three of the four detected epistatic interactions. Candidate genes at this QTL could be identified through the integration with genome data and included the presynaptic active zone protein bruchpilot. Our findings contribute to understanding the genomic control of seasonal diapause timing in an ecological relevant species.

Rapid and repeatable host plant shifts drive reproductive isolation following a recent human-mediated introduction of the apple maggot fly, Rhagoletis pomonella

Ecological speciation via host-shifting is often invoked as a mechanism for insect diversification, but the relative importance of this process is poorly understood. The shift of Rhagoletis pomonella in the 1850s from the native downy hawthorn, Crataegus mollis, to introduced apple, Malus pumila, is a classic example of sympatric host race formation, a hypothesized early stage of ecological speciation. The accidental human-mediated introduction of R. pomonella into the Pacific Northwest (PNW) in the late 1970s allows us to investigate how novel ecological opportunities may trigger divergent adaptation and host race formation on a rapid timescale. Since the introduction, the fly has spread in the PNW, where in addition to apple, it now infests native black hawthorn, Crataegus douglasii, and introduced ornamental hawthorn, Crataegus monogyna. We use this "natural experiment" to test for genetic differentiation among apple, black, and ornamental hawthorn flies co-occurring at three sympatric sites. We report evidence that populations of all three host-associations are genetically differentiated at the local level, indicating that partial reproductive isolation has evolved in this novel habitat. Our results suggest that conditions suitable for initiating host-associated divergence may be common in nature, allowing for the rapid evolution of new host races when ecological opportunity arises.

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 F_ST 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.

Genomic transitions during host race and species formation

Darwin recognized species as discontinuous, yet considered them to be formed by an incremental process of natural selection. Recent theoretical work on 'genome-wide congealing' is bridging this gap between the gradualism of divergent selection and rapid genome-wide divergence, particularly during ecological speciation-with-gene-flow. Host races and species of phytophagous insects, displaying a spectrum of divergence and gene flow among member taxa, provide model systems for testing predicted non-linear transitions from 'genic' divergence at a few uncoupled loci to 'genomic' divergence with genome-wide coupling of selected loci and strong reproductive isolation. Integrating across natural history, genomics, and evolutionary theory, emerging research suggests a tipping point from 'genic' to 'genomic' divergence between host races and species, during both sympatric speciation and secondary contact.

Standing geographic variation in eclosion time and the genomics of host race formation in Rhagoletis pomonella fruit flies

Taxa harboring high levels of standing variation may be more likely to adapt to rapid environmental shifts and experience ecological speciation. Here, we characterize geographic and host-related differentiation for 10,241 single nucleotide polymorphisms in Rhagoletis pomonella fruit flies to infer whether standing genetic variation in adult eclosion time in the ancestral hawthorn (Crataegus spp.)-infesting host race, as opposed to new mutations, contributed substantially to its recent shift to earlier fruiting apple (Malus domestica). Allele frequency differences associated with early vs. late eclosion time within each host race were significantly related to geographic genetic variation and host race differentiation across four sites, arrayed from north to south along a 430-km transect, where the host races co-occur in sympatry in the Midwest United States. Host fruiting phenology is clinal, with both apple and hawthorn trees fruiting earlier in the North and later in the South. Thus, we expected alleles associated with earlier eclosion to be at higher frequencies in northern populations. This pattern was observed in the hawthorn race across all four populations; however, allele frequency patterns in the apple race were more complex. Despite the generally earlier eclosion timing of apple flies and corresponding apple fruiting phenology, alleles on chromosomes 2 and 3 associated with earlier emergence were paradoxically at lower frequency in the apple than hawthorn host race across all four sympatric sites. However, loci on chromosome 1 did show higher frequencies of early eclosion-associated alleles in the apple than hawthorn host race at the two southern sites, potentially accounting for their earlier eclosion phenotype. Thus, although extensive clinal genetic variation in the ancestral hawthorn race exists and contributed to the host shift to apple, further study is needed to resolve details of how this standing variation was selected to generate earlier eclosing apple fly populations in the North.

Habitat or temporal isolation: Unraveling herbivore-parasitoid speciation patterns using double digest RADseq

Ecological speciation is often observed in phytophagous insects and their parasitoids due to divergent selection caused by host-associated or temporal differences. Most previous studies have utilized limited genetic markers or distantly related species to look for reproductive barriers of speciation. In our study, we focus on closely related species of Lygus bugs and two sister species of Peristenus parasitoid wasps. Using mitochondrial DNA COI and genomewide SNPs generated using ddRADseq, we tested for potential effects of host-associated differentiation (HAD) or temporal isolation in this system. While three species of Lygus are clearly delineated with both COI and SNPs, no evidence of HAD or temporal differentiation was detected. Two Peristenus sister species were supported by both sets of markers and separated temporally, with P. mellipes emerging early in June and attacking the first generation of Lygus, while P. howardi emerging later in August and attacking the second generation of their hosts. This is one of the few studies to examine closely related hosts and parasitoids to examine drivers of diversification. Given the results of this study, the Lygus-Peristenus system demonstrates temporal isolation as a potential barrier to reproductive isolation for parasitoids, which could indicate higher parasitoid diversity in regions of multivoltine hosts. This study also demonstrates that incorporating systematics improves studies of parasitoid speciation, particularly by obtaining accurate host records through rearing, carefully delimiting cryptic species and examining population-level differences with genomic-scale data among closely related taxa.

Loci under selection and markers associated with host plant and host-related strains shape the genetic structure of Brazilian populations of Spodoptera frugiperda (Lepidoptera, Noctuidae)

We applied the ddRAD genotyping-by-sequencing technique to investigate the genetic distinctiveness of Brazilian populations of the noctuid moth Spodoptera frugiperda, the fall armyworm (FAW), and the role of host-plant association as a source of genetic diversification. By strain-genotyping all field-collected individuals we found that populations collected from corn were composed primarily of corn-strain individuals, while the population collected from rice was composed almost entirely of rice-strain individuals. Outlier analyses indicated 1,184 loci putatively under selection (ca. 15% of the total) related to 194 different Gene Ontologies (GOs); the most numerous GOs were nucleotide binding, ATP binding, metal-ion binding and nucleic-acid binding. The association analyses indicated 326 loci associated with the host plant, and 216 loci associated with the individual strain, including functions related to Bacillus thuringiensis and insecticide resistance. The genetic-structure analyses indicated a moderate level of differentiation among all populations, and lower genetic structure among populations collected exclusively from corn, which suggests that the population collected from rice has a strong influence on the overall genetic structure. Populations of S. frugiperda are structured partially due to the host plant, and pairs of populations using the same host plant are more genetically similar than pairs using different hosts. Loci putatively under selection are the main factors responsible for the genetic structure of these populations, which indicates that adaptive selection on important traits, including the response to control tactics, is acting in the genetic differentiation of FAW populations in Brazil.

Genomic Differentiation during Speciation-with-Gene-Flow: Comparing Geographic and Host-Related Variation in Divergent Life History Adaptation in Rhagoletis pomonella

A major goal of evolutionary biology is to understand how variation within populations gets partitioned into differences between reproductively isolated species. Here, we examine the degree to which diapause life history timing, a critical adaptation promoting population divergence, explains geographic and host-related genetic variation in ancestral hawthorn and recently derived apple-infesting races of Rhagoletis pomonella. Our strategy involved combining experiments on two different aspects of diapause (initial diapause intensity and adult eclosion time) with a geographic survey of genomic variation across four sites where apple and hawthorn flies co-occur from north to south in the Midwestern USA. The results demonstrated that the majority of the genome showing significant geographic and host-related variation can be accounted for by initial diapause intensity and eclosion time. Local genomic differences between sympatric apple and hawthorn flies were subsumed within broader geographic clines; allele frequency differences within the races across the Midwest were two to three-fold greater than those between the races in sympatry. As a result, sympatric apple and hawthorn populations displayed more limited genomic clustering compared to geographic populations within the races. The findings suggest that with reduced gene flow and increased selection on diapause equivalent to that seen between geographic sites, the host races may be recognized as different genotypic entities in sympatry, and perhaps species, a hypothesis requiring future genomic analysis of related sibling species to R. pomonella to test. Our findings concerning the way selection and geography interplay could be of broad significance for many cases of earlier stages of divergence-with-gene flow, including (1) where only modest increases in geographic isolation and the strength of selection may greatly impact genetic coupling and (2) the dynamics of how spatial and temporal standing variation is extracted by selection to generate differences between new and discrete units of biodiversity.

Sympatric Speciation in the Genomic Era

Sympatric speciation has been of key interest to biologists investigating how natural and sexual selection drive speciation without the confounding variable of geographic isolation. The advent of the genomic era has provided a more nuanced and quantitative understanding of the different and often complex modes of speciation by which sympatric sister taxa arose, and a reassessment of some of the most compelling empirical case studies of sympatric speciation. However, I argue that genomic studies based on contemporary populations may never be able to provide unequivocal evidence of true primary sympatric speciation, and there is a need to incorporate palaeogenomic studies into this field. This inability to robustly distinguish cases of primary and secondary 'divergence with gene flow' may be inconsequential, as both are useful for understanding the role of large effect barrier loci in the progression from localised genic isolation to genome-wide reproductive isolation. I argue that they can be of equivalent interest due to shared underlying mechanisms driving divergence and potentially leaving similar patterns of coalescence.