Host plant and latitude-related diapause variation in Rhagoletis pomonella: a test for multifaceted life history adaptation on different stages of diapause development

Variation in Thermal Sensitivity of Diapause Development among Individuals and over Time Predicts Life History Timing in a Univoltine Insect

AbstractPhysiological time is important for understanding the development and seasonal timing of ectothermic animals but has largely been applied to developmental processes that occur during spring and summer, such as morphogenesis. There is a substantial knowledge gap in the relationship between temperature and development during winter, a season that is increasingly impacted by climate change. Most temperate insects overwinter in diapause, a developmental process with little obvious morphological change. We used principles from the physiological time literature to measure and model the thermal sensitivity of diapause development rate in the apple maggot fly Rhagoletis pomonella, a univoltine fly whose diapause duration varies substantially within and among populations. We show that diapause duration can be predicted by modeling a relationship between temperature and development rate that is shifted toward lower temperatures compared with typical models of morphogenic, nondiapause development. However, incorporating interindividual variation and ontogenetic variation in the temperature-to-development rate relationship was critical for accurately predicting fly emergence, as diapause development proceeded more quickly at high temperatures later in diapause. We conclude that the conceptual framework may be flexibly applied to other insects and discuss possible mechanisms of diapause timers and implications for phenology with warming winters.

Crowd-sourced observations of a polyphagous moth reveal evidence of allochronic speciation varying along a latitudinal gradient

Allochronic speciation, where reproductive isolation between populations of a species is facilitated by a difference in reproductive timing, depends on abiotic factors such as seasonality and biotic factors such as diapause intensity. These factors are strongly influenced by latitudinal trends in climate, so we hypothesized that there is a relationship between latitude and divergence among populations separated by life history timing. Hyphantria cunea (the fall webworm), a lepidopteran defoliator with red and black colour morphs, is hypothesized to be experiencing an incipient allochronic speciation. However, given their broad geographic range, the strength of allochronic speciation may vary across latitude. We annotated >11,000 crowd-sourced observations of fall webworm to model geographic distribution, phenology, and differences in colour phenotype between morphs across North America. We found that red and black morph life history timing differs across North America, and the phenology of morphs diverges more in warmer climates at lower latitudes. We also found some evidence that the colour phenotype of morphs also diverges at lower latitudes, suggesting reduced gene flow between colour morphs. Our results demonstrate that seasonality in lower latitudes may increase the strength of allochronic speciation in insects, and that the strength of sympatric speciation can vary along a latitudinal gradient. This has implications for our understanding of broad-scale speciation events and trends in global biodiversity.

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.

Rapid brain development and reduced neuromodulator titres correlate with host shifts in Rhagoletis pomonella

Host shifts are considered a key generator of insect biodiversity. For insects, adaptation to new host plants often requires changes in larval/pupal development and adult behavioural preference toward new hosts. Neurochemicals play key roles in both development and behaviour and therefore provide a potential source for such synchronization. Here, we correlated life-history timing, brain development and corresponding levels of 14 neurochemicals in Rhagoletis pomonella (Diptera: Tephritidae), a species undergoing ecological speciation through an ongoing host shift from hawthorn to apple fruit. These races exhibit differences in pupal diapause timing as well as adult behavioural preference with respect to their hosts. This difference in behavioural preference is coupled with differences in neurophysiological response to host volatiles. We found that apple race pupae exhibited adult brain morphogenesis three weeks faster after an identical simulated winter than the hawthorn race, which correlated with significantly lower titres of several neurochemicals. In some cases, particularly biogenic amines, differences in titres were reflected in the mature adult stage, when host preference is exhibited. In summary, life-history timing, neurochemical titre and brain development can be coupled in this speciating system, providing new hypotheses for the origins of new species through host shifts.

Genotoxicity Set Up in Artemia franciscana Nauplii and Adults Exposed to Phenanthrene, Naphthalene, Fluoranthene, and Benzo(k)fluoranthene

Polycyclic aromatic hydrocarbons (PAHs) consist of a group of over 100 different organic compounds mainly generated by anthropogenic activities. Because of their low water solubility, they tend to be accumulated in sediment, where their degradation rate is very low. Few studies have been carried out so far to investigate the effects of PAHs on Artemia franciscana. Artemia is easy to manage at laboratory scale, but it is not a sensitive biological model considering the traditional endpoints (i.e., mortality). In addition to evaluating the lethality on nauplii and adults of A. franciscana after 24 and 48 h, we focused on the genotoxicity to investigate the potential effects of phenanthrene (PHE), naphthalene (NAP), fluoranthene (FLT), and benzo(k)fluoranthene (BkF). Results showed that FLT was the most toxic both for nauplii and adults after 48 h of exposure. Real-time qPCR showed that all toxicants, including BkF, which had no negative effects on the survival of the crustacean, were able to switch the gene expression of all nine genes. This work has important ecological implications, especially on contaminated sediment assessment considering that PAHs represent the most abundant organic group of compounds in marine environment, opening new perspectives in understanding the molecular pathways activated by crustaceans.

Sensitivities to Chill Durations and No-Chill Temperatures Regulating Eclosion Responses Differ Between Rhagoletis zephyria (Diptera: Tephritidae) and its Braconid Parasitoids (Hymenoptera: Braconidae)

Seasonal temperatures select for eclosion timing of temperate insects and their parasitoids. In western North America, the fruit fly Rhagoletis zephyria Snow (Diptera: Tephritidae) is parasitized by the hymenopterous wasps Utetes lectoides (Gahan), an egg parasite, and Opius downesi Gahan, a larval parasite (both Braconidae). Eclosion of wasps should be timed with the presence of susceptible fly stages, but reports indicate U. lectoides ecloses in the absence of flies under no-chill conditions. Based on this, we tested the hypotheses that chill durations and no-chill temperatures both differentially regulate eclosion times of R. zephyria and its parasitic wasps. When fly puparia were chilled at ~3°C for 130-180 d, U. lectoides and O. downesi always eclosed on average later than flies. However, after 180-d chill, flies eclosed on average earlier than after 130- and 150-d chill, whereas eclosion times of U. lectoides and O. downesi were less or not affected by chill duration. When fly puparia were exposed to 20-22°C (no chill), U. lectoides eclosed before flies, with 88.9% of U. lectoides versus only 0.61% of flies eclosing. Taken together, findings show that eclosion times of flies are more sensitive to changes in chill duration than those of wasps. Flies are less sensitive than wasps to no-chill in that most flies do not respond by eclosing after no-chill while most wasps do. Our results suggest that shorter winters and longer summers due to climate change could cause mismatches in eclosion times of flies and wasps, with potentially significant evolutionary consequences.

Host plant stimuli effects on survival and duration of dormancy of Rhagoletis zoqui, Rhagoletis completa (Diptera: Tephritidae) and associated parasitoids

Tephritid fruit flies in the genus Rhagoletis bridge between predictable periods of fruit availability by becoming dormant. To cope with acyclic unpredictable events (e.g., frost, mast seeding, etc), a proportion of the population can undergo prolonged dormancy. In the case of walnut infesting Rhagoletis, host plant-derived cues such as juglone soil concentration vary seasonally in predictable patterns. Here, we examined the effects of host plant parts and derived compounds on emergence rates and dormancy duration of Rhagoletis completa (Cresson), Rhagoletis zoqui (Bush) (Diptera: Tephritidae), and associated parasitoids. Pupae of both species were exposed to walnut leaves, fruit, or fruit and leaves and compared to a control. In a second experiment, R. zoqui were exposed to 10 mg l-1 of juglone applied to pupation medium during four consecutive 4-week time periods under variable combinations of temperature and frequency of exposure. Overall, the presence of fruit resulted in greater overwintering survival of R. completa but had no effect on the duration of dormancy of either fly species. Application of juglone over two consecutive periods produced greater mortality of R. zoqui than the control. Three parasitoid species emerged from R. completa and one from R. zoqui. Duration of dormancy for parasitoids was longer than that of fly hosts. Regardless of treatment, 13.3-18.4% of R. completa pupae and 1.3-2.8% R. zoqui engaged in prolonged (>year) dormancy. Our results indicate that host plant derived cues have little or no effect on survival and duration of dormancy of walnut infesting Rhagoletis, and at the tested concentration juglone is toxic to R. zoqui pupae. Testing the effect of juglone at lower concentrations is necessary to rule out its role as an environmental cue for regulation of dormancy. So far, host plant fruiting phenology appears to play a greater role than host plant derived cues in selecting for fly eclosion times.

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.

Diapause differentially modulates the transcriptomes of fat body and flight muscle in the Colorado potato beetle

Many temperate insects, such as the Colorado potato beetle, enter diapause in winter, during which they arrest their development, suppress their metabolic rate and have high stress tolerance. Diapause phenotypes can be transcriptionally regulated, however many studies to date report only whole animal gene expression rather than tissue-specific processes during diapause. We used RNA-seq to measure gene expression in fat body and flight muscle of diapausing and non-diapausing beetles. We used differential expression and GO enrichment analyses to evaluate longstanding hypotheses about the mechanisms that drive arrested development, changes in energy metabolism, and increased stress tolerance during diapause. We found evidence of G2/M cell cycle arrest, juvenile hormone catabolism, increased antioxidant metabolism, epigenetic modification, transposable element regulation, and cytoskeletal remodeling in both the fat body and flight muscle of diapausing beetles. Beetles differentially modulated the fat body and flight muscle transcriptomes during diapause with fat body playing a larger role in the hypoxia response and immunity, whereas flight muscle had higher abundance of transcripts related to the chaperone response and proteostasis. Our transcriptome provides evidence for distinct roles and responses of fat body and flight muscle during diapause in the Colorado potato beetle, and we provide testable hypotheses for biological processes that appear to drive diapause phenotypes in insects.

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.

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.

Does Drosophila sechellia escape parasitoid attack by feeding on a toxic resource?

Host shifts can drastically change the selective pressures that animals experience from their environment. Drosophila sechellia is a species restricted to the Seychelles islands, where it specializes on the fruit Morinda citrifolia (noni). This fruit is known to be toxic to closely related Drosophila species, including D. melanogaster and D. simulans, releasing D. sechellia from interspecific competition when breeding on this substrate. Previously, we showed that larvae of D. sechellia are unable to mount an effective immunological response against wasp attack, while larvae of closely-related species can defend themselves from parasitoid attack by melanotic encapsulation. We hypothesized that this inability constitutes a trait loss due to a reduced risk of parasitoid attack in noni. Here we present a lab experiment and field survey aimed to test the hypothesis that specialization on noni has released D. sechellia from the antagonistic interaction with its larval parasitoids. Our results from the lab experiment suggest that noni may be harmful to parasitoid wasps. Our results from the field survey indicate that D. sechellia was found in ripe noni, whereas another Drosophila species, D. malerkotliana, was present in unripe and overripe stages. Parasitic wasps of the species Leptopilina boulardi emerged from overripe noni, where D. malerkotliana was the most abundant host, but not from ripe noni. These results indicate that the specialization of D. sechellia on noni has indeed drastically altered its ecological interactions, leading to a relaxation in the selection pressure to maintain parasitoid resistance.

Effects of photoperiod and relative humidity on diapause termination and post-winter development of Rhagoletis cerasi pupae

The European cherry fruit fly, Rhagoletis cerasi (Diptera: Tephritidae), is a univoltine species that undergoes obligatory summer-winter diapause at pupal stage in the soil (2-5 cm) beneath host trees. To study the effects of photoperiod and relative humidity on diapause termination and post-winter developmental duration of R. cerasi, pupae collected from Dossenheim (Germany) were exposed to different photoperiod or relative humidity regimes during a chilling period ranging from 2 to 8.5 months. Specifically, pupae were exposed to four photoperiod regimes: (a) light conditions (24L:00D), (b) dark conditions (00L:24D), (c) short photoperiod (08L:16D) and (d) long photoperiod (16L:08D), as well as to three relative humidity regimes: (a) low (40% RH), (b) medium (60% RH) and (c) high (70-80% RH). Data revealed that relative humidity is not a significant predictor of diapause termination, but it affects the post-winter developmental period. Higher relative humidity promotes post-winter pupae development. On the other hand, photoperiod significantly affected both diapause termination and post-winter development of R. cerasi pupae. Light conditions (24L:00D) accelerate adult emergence, particularly for females. Regardless of the photoperiod (24L:00D, 00L:24D, 08L:16D), rates of adult emergence were high (>75%) for chilling intervals longer than 6.5 months. Nonetheless, exposure to a long day photoperiod (16L:08D), during chilling, dramatically reduced the proportion of adult emergence following 6 months exposure to chilling. Our findings broaden the understanding of factors regulating diapause responses in European cherry fruit fly, local adaptation and synchronization of adult emergence with the ripening period of major hosts.

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.

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

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.

Effects of Pre-Diapause Temperature and Body Weight on the Diapause Intensity of the Overwintering Generation of Bactrocera minax (Diptera: Tephritidae)

The Chinese citrus fruit fly, Bactrocera minax (Enderlein), is an economically important pest of citrus. The fly has an obligatory pupal diapause in soil from November to March. However, techniques for predicting or determining the emergence of the adult have, thus far, not been well documented. In this study, we investigated the effects of different pre-diapause temperatures (8, 12, 16, and 20°C) and pupal body weight (five groups according to pupal weight: G-58, 55.0-61.0 mg; G-68, 65-71 mg; G-78, 75-81 mg; G-88, 85-91 mg; G-95, 92-98 mg) on pupal period (the indicator of diapause intensity). The pupal period of B. minax larvae pupated at 8°C was 193.41 d, which was significantly shorter than that of larvae incubated at higher temperatures, suggesting that there was a lower diapause intensity for larvae pupated at lower pre-diapause temperatures. There were also significant differences in the pupal periods at different pupal body weights. The pupal period of G-58 was significantly shorter than that of the heavier groups (G-88 and G-95), and the pupal period increased with increasing pupal body weight in the five groups. Moreover, the pupal period of B. minax significantly and positively correlated to pupal body weight. These findings demonstrate that the pre-diapause temperature and pupal body weight are suitable indicators for predicting the pupal period of overwintering individuals, and the results of this study will contribute to the development of new and effective strategies for predicting the occurrence and population dynamics of B. minax adult.

Strong divergent selection at multiple loci in two closely related species of ragworts adapted to high and low elevations on Mount Etna

Recently diverged species present particularly informative systems for studying speciation and maintenance of genetic divergence in the face of gene flow. We investigated speciation in two closely related Senecio species, S. aethnensis and S. chrysanthemifolius, which grow at high and low elevations, respectively, on Mount Etna, Sicily and form a hybrid zone at intermediate elevations. We used a newly generated genome-wide single nucleotide polymorphism (SNP) dataset from 192 individuals collected over 18 localities along an elevational gradient to reconstruct the likely history of speciation, identify highly differentiated SNPs, and estimate the strength of divergent selection. We found that speciation in this system involved heterogeneous and bidirectional gene flow along the genome, and species experienced marked population size changes in the past. Furthermore, we identified highly-differentiated SNPs between the species, some of which are located in genes potentially involved in ecological differences between species (such as photosynthesis and UV response). We analysed the shape of these SNPs' allele frequency clines along the elevational gradient. These clines show significantly variable coincidence and concordance, indicative of the presence of multifarious selective forces. Selection against hybrids is estimated to be very strong (0.16-0.78) and one of the highest reported in literature. The combination of strong cumulative selection across the genome and previously identified intrinsic incompatibilities probably work together to maintain the genetic and phenotypic differentiation between these species - pointing to the importance of considering both intrinsic and extrinsic factors when studying divergence and speciation.

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.

Changes in Energy Metabolism Trigger Pupal Diapause Transition of Bactrocera minax After 20-Hydroxyecdysone Application

Correct timing of diapause entry and exit is critical for a species' survival. While many aspects of insect diapause are well-studied, the mechanisms underlying diapause termination remain largely unknown. The Chinese citrus fly, Bactrocera minax, is a univoltine insect with an obligatory pupal diapause. The application of 20-hydroxyecdysone (20E) is known to terminate diapause in B. minax, and we used this approach, along with isobaric tags for relative and absolute quantitation technology, to determine the proteins associated with diapause termination in this fly. Among 2,258 identified proteins, 1,169 proteins significantly differed at 1, 2, and 5 days post-injection of 20E, compared with the solvent-injected control group. Functional annotation revealed that the majority of differentially expressed proteins were enriched in the core energy metabolism of amino acids, proteins, lipids, and carbohydrates as well as in signal transduction pathways including PPAR signaling, Calcium signaling, Glucagon signaling, VEGF signaling, Ras signaling, cGMP-PKG signaling, and cAMP signaling. A combined transcriptomic and proteomic analysis suggested the involvement of energy metabolism in the response of diapause transition. RNA interference experiments disclosed that a 20E injection triggers diapause termination probably through non-genomic actions, rather than nuclear receptor mediated genomic actions. Our results provide extensive proteomic resources for insect diapause transition and offer a potential for pest control by incapacitating the regulation of diapause termination either by breaking diapause prematurely or by delaying diapause termination to render diapausing individuals at a high risk of mortality.

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.

Differential Response to Photoperiod by Diverging Rhagoletis pomonella (Diptera: Tephritidae) Populations Exploiting Host Plants with Different Fruiting Phenology

Photoperiod has been found to influence the proportion of non-dormant individuals and the duration of dormancy among North American populations of Rhagoletis pomonella (Walsh) (Diptera: Tephritidae). In particular, long days combined with high temperatures can produce a 100% non-dormant generation. There are several genetically distinct populations of R. pomonella in the highlands of Mexico, two of which occur at the same latitude, at different elevations, that exploit hawthorn plants with different fruiting phenology (early (summer-fall) and late (winter)). Flies exploiting such hosts might use day length in different ways as a cue to match adult emergence with fruit availability. Here, we compared responses of pupae from two Mexican populations of R. pomonella to a 12/12 L/D photoperiod, a long-day photoperiod 17/7 L/D, and continuous illumination. Experiments were performed under warm conditions (27°C). Day length had no effect on the proportion of adults emerged from the early-fruiting hawthorn population, while pupae extracted from late-fruiting hawthorns and exposed to long days emerged as adults in lower proportions and engaged in prolonged dormancy in greater proportions than those exposed to a short day. Photoperiod had no effect on the proportion from both Mexican populations foregoing dormancy. Dormancy duration was affected by photoperiod and was longer than previous reports for North American populations. Host plant fruiting phenology may be the driver of these differences. These findings, coupled with previous reports of genetic and biological differences, suggest that the Eje Volcanico Trans-Mexicano R. pomonella population may actually be a distinct species.

First Insights into the Intrapuparial Development of Bactrocera dorsalis (Hendel): Application in Predicting Emergence Time for Tephritid Fly Control

Intrapuparial development is a special pattern of metamorphosis in cyclorrhaphous flies, in which the pupa forms in an opaque, barrel-like puparium. This has been well studied in forensic insects for age estimations. In this study, the intrapuparial development of a quarantine agricultural pest, Bactrocera dorsalis (Hendel), was studied under a constant temperature of 27 ± 1 °C and 70 ± 5% relative humidity. Results showed that intrapuparial development could be divided into five stages: Larval-pupal apolysis, cryptocephalic pupa, phanerocephalic pupa, pharate adult, and emergent adult. It lays a morphology-based foundation for molecular mechanism studies and enhances the understanding of the physiological basis for changes in intrapuparial development. More importantly, the chronology of intrapuparial development can be used to predict the emergence time of tephritid flies, indicating when to spray insecticides to control these phytophagous agricultural pests. This may be an effective approach to reduce the use of insecticides and slow down the evolution of insecticidal resistance.

The Effect of Winter Length on Duration of Dormancy and Survival of Rhagoletis completa (Diptera: Tephritidae) and Associated Parasitoids From Northeastern Mexico

The walnut husk fly Rhagoletis completa (Cresson), native to the Midwestern United States and Mexico, is invasive in California and Europe. It is one of the most important pests of walnuts in areas gathering 30% of the world production. Knowledge of life-history regulation is important for the design of management strategies. Research on dormancy has been performed on invasive populations, and not on populations at the southern extreme of its native range. Here, we examined the effect of winter length on fly and parasitoid emergence, survival, and duration of dormancy. Percent emergence was higher for chill periods at 5°C ranging from 8 to 20 wk. No or insufficient chill resulted in low emergence and a significant proportion of individuals in prolonged dormancy (>1 yr). Duration of dormancy was longer for pupae at constant temperatures and a 4-wk chill period than longer winter durations. Dormancy was longer for Mexican than that reported for U.S. populations, suggesting the existence of a latitudinal cline where populations at southern latitudes have evolved slower metabolic rates. Three parasitoid species were found associated with R. completa (Aganaspis alujai (Wharton and Ovruski) (Hymenoptera: Figitidae), Diachasmimorpha juglandis Muesebeck, and Diachasmimorpha mellea Gahan) (Hymenoptera: Braconidae). Results suggest that rearing of R. completa is possible by subjecting pupae to chill periods between 8 and 20 wk. Overwintering mortality of flies and A. alujai could be further reduced above 5°C. Our findings can contribute for the accurate development of predictive models on invasion potential, development, fly and parasitoid rearing, and biological control.

Mate choice and gene expression signatures associated with nutritional adaptation in the medfly (Ceratitis capitata)

Evolutionary responses to nutrition are key to understanding host shifts and the resulting potential for reproductive isolation. Experimental evolution has previously been used to describe the responses of the medfly (Ceratitis capitata) to larval diets with different nutritional properties. Within 30 generations this led to divergence in larval development time, egg to adult survival and adaptation in adult body size. Here we used mRNA-seq to identify differences in gene expression patterns in these same populations, using males from the 60^th generation of nutritional selection. We validated differential expression by using qRT-PCR and found that genes linked to metabolism, oxidative phosphorylation and proteolysis were significantly over-represented among the differentially expressed genes. The results provide the first genome-wide survey of the putative mechanisms underpinning evolved responses to nutritional adaptation. In addition, we tested the hypothesis that nutritional adaptation can alter mating patterns. We found evidence for assortative mating by diet at generation 60, but not 90. Hence, the pattern was variable across generations and there was no evidence overall for any isolating mating divergence between the lines. Overall, the results provide insight into the mechanisms underpinning dietary adaptation and extend our knowledge of which traits represent core responses to nutritional selection.

Estimating the impact of divergent mating phenology between residents and migrants on the potential for gene flow

Gene flow between populations can allow the spread of beneficial alleles and genetic diversity between populations, with importance to conservation, invasion biology, and agriculture. Levels of gene flow between populations vary not only with distance, but also with divergence in reproductive phenology. Since phenology is often locally adapted, arriving migrants may be reproductively out of synch with residents, which can depress realized gene flow. In flowering plants, the potential impact of phenological divergence on hybridization between populations can be predicted from overlap in flowering schedules-the daily count of flowers capable of pollen exchange-between a resident and migrant population. The accuracy of this prospective hybridization estimate, based on parental phenotypes, rests upon the assumptions of unbiased pollen transfer between resident and migrant active flowers. We tested the impact of phenological divergence on resident-migrant mating frequencies in experiments that mimicked a single large gene flow event. We first prospectively estimated mating frequencies two lines of Brassica rapaselected or early and late flowering. We then estimated realized mating frequencies retrospectively through progeny testing. The two estimates strongly agreed in a greenhouse experiment, where procedures ensured saturating, unbiased pollination. Under natural pollination in the field, the rate of resident-migrant mating, was lower than estimated by phenological divergence alone, although prospective and retrospective estimates were correlated. In both experiments, differences between residents and migrants in flowering schedule shape led to asymmetric hybridization. Results suggest that a prospective estimate of hybridization based on mating schedules can be a useful, although imperfect, tool for evaluating potential gene flow. They also illustrate the impact of mating phenology on the magnitude and symmetry of reproductive isolation.

Sources of Variation in the Adult Flight of Walnut Husk Fly (Diptera: Tephritidae): A Phenology Model for California Walnut Orchards

A phenology model of the walnut husk fly, Rhagoletis completa Cresson, was developed to more accurately predict the timing of the flight period and optimize management decisions. A data set of 153 orchard years in which adults were trapped throughout the season was used for the development and validation of this model. Data from California Irrigation Management Information System (CIMIS) weather stations were used to match orchard-year datasets with historical climatic data on degree-day (DD) accumulation, winter chill, and winter rainfall. A cumulative Weibull distribution was used to model the relationship between cumulative trap catch and DD accumulation for R. completa in California. The model was used to predict thermal requirements for the start (5% cumulative trap catch) and mid-point (50% cumulative trap catch) of the flight period, which were 1,670 and 2,179 DDs, respectively. The prediction for 50% cumulative trap catch of R. completa in California was much higher than the thermal requirement estimated in Oregon previously (1,751 DDs). Linear mixed effects models were used to evaluate other environmental and orchard-specific factors which could explain the large variation between predicted and observed thermal requirements for both the start and mid-point of the flight period. Latitude, walnut cultivar leaf-out time, orchard age and year, as a continuous variable, all contributed significantly to explain deviations from the predictions of the DD model for individual orchard years. Such factors can be used both to adjust predicted thermal requirements for these two specific and informative stages of the flight period, and to provide a basis for ecological and evolutionary hypotheses.

Urban heat islands advance the timing of reproduction in a social insect

For many species, the timing of life cycle events is advancing under contemporary global climate change. However, much less is known regarding phenological shifts as a result of other sources of anthropogenic change, such as urban warming. In both cases, progress has been hampered by a focus on phenological traits such as the timing of emergence, rather than the phenology of more directly related fitness traits such as the timing of reproduction. Here we explore how urban heat island effects shape the timing of reproduction in an acorn-dwelling ant. We used a common garden experiment with acorn ants collected from three cities in the eastern United States along a latitudinal gradient and reared long-term in the laboratory under five temperature treatments. This allowed us to quantify the effects of temperature on reproductive phenology across three scales-a biogeographic temperature cline, three urban vs. rural temperature comparisons, and five laboratory rearing temperatures. At our northernmost and southernmost cities (spanning 6° of latitude), we found both urbanization and warmer laboratory rearing temperature significantly advanced reproductive phenology; ants from the lowest latitude city also had earlier reproductive phenology compared with the higher latitude cities. In the field, the differences in urban versus rural acorn ant reproductive phenology translate to approximately one month earlier reproduction in the urban populations. For insects with synchronous mating events, such as ants, shifts in the already short window of time to reproduce could limit mating across environments, potentially leading to reproductive isolation between urban and rural populations.

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.

Divergence in larval diapause induction between the rice and water-oat populations of the striped stem borer, Chilo suppressalis (Walker) (Lepidoptera: Crambidae)

Differences in diapause traits can result in the seasonal reproductive isolation of host plant-associated insect populations and thereby facilitate the population divergence. The striped stem borer, Chilo suppressalis, has two host plant-associated populations: rice population and water-oat population. Several studies have found evidence that seasonal reproductive isolation between these populations is at least partially due to interpopulation differences in diapause. However, there still lack unambiguous evidence comparing characteristics of diapause induction for both populations. We compared the photoperiodic response and the age of peak photoperiod sensitivity of these populations and used RNA-Seq to compare the molecular response of diapause induction between populations. The photoperiodic response of the two populations differed at 25 °C; the critical night length of larvae from the rice population was 11 h and 20 min, whereas no obvious critical night length was in those from the water-oat population. In rice population, larvae were most sensitive to photoperiod at 9-12 days of age, whereas in water-oat population, larvae were the most sensitive to photoperiod at 9-10 days of age. The RNA-Seq results indicated that there were several differences in the molecular response of diapause induction and small overlap in differentially expressed genes (DEGs) between populations. Furthermore, GO analysis indicated that both rice and water-oat population's DEGs were significantly enriched in heme and iron binding. Besides, water-oat population's DEGs were significantly enriched in metabolizing nutrients but rice population's DEGs do not. Thus, our results described differences in diapause induction between rice and water-oat populations of C. suppressalis which could affect the timing of diapause and thereby contribute to the seasonal reproductive isolation of these host plant-associated populations. In conclusion, this work suggests that difference in diapause induction could promote the population divergence in insects associated with different host plants.

The effect of winter length on duration of dormancy and survival of specialized herbivorous Rhagoletis fruit flies from high elevation environments with acyclic climatic variability

Dormancy can be defined as a state of suppressed development allowing insects to cope with adverse conditions and plant phenology. Among specialized herbivorous insects exploiting seasonal resources, diapause frequently evolves as a strategy to adjust to predictable plant seasonal cycles. To cope with acyclic and unpredictable climatic events, it has been found for some insects that a proportion of the population undergoes prolonged dormancy. We compared the response of three species in the Rhagoletis cingulata species group exploiting plants differing in fruiting phenology from environments varying in frequency and timing of acyclic climatic catastrophic events (frost during flowering and fruit set) and varying also in the time of the onset of the rainy season. Small proportions (10 months), and large proportions of pupae died without emerging as adults. The number of days elapsed from the end of artificial winter and adult eclosion was longer for R. cingulata exploiting late fruiting Prunus serotina in Northeastern Mexico than for flies recovered from earlier fruiting plants in the central Altiplano. Rhagoletis turpiniae and northeastern R. cingulata pupae suffered high proportions of parasitism. Large proportions of R. cingulata from central Mexico engaging in prolonged dormancy may be explained by the fact that flowering and fruit set for its host, P. serotina var capuli, driven by the timing of maximum precipitation, matches a period of highest probability of frost often resulting in large areas with fruitless trees at unpredictable time intervals. As a consequence of differences in host plant fruiting phenology, central and northeastern Mexican R. cingulata were found to be allochronically isolated. Prolonged dormancy may have resulted in escape from parasitism.

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.

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.

Phenological mismatch and the effectiveness of assisted gene flow

The persistence of narrowly adapted species under climate change will depend on their ability to migrate apace with their historical climatic envelope or to adapt in place to maintain fitness. This second path to persistence can only occur if there is sufficient genetic variance for response to new selection regimes. Inadequate levels of genetic variation can be remedied through assisted gene flow (AGF), that is the intentional introduction of individuals genetically adapted to localities with historic climates similar to the current or future climate experienced by the resident population. However, the timing of reproduction is frequently adapted to local conditions. Phenological mismatch between residents and migrants can reduce resident × migrant mating frequencies, slowing the introgression of migrant alleles into the resident genetic background and impeding evolutionary rescue efforts. Focusing on plants, we devised a method to estimate the frequency of resident × migrant matings based on flowering schedules and applied it in an experiment that mimicked the first generation of an AGF program with Chamaecrista fasciculata, a prairie annual, under current and expected future temperature regimes. Phenological mismatch reduced the potential for resident × migrant matings by 40-90%, regardless of thermal treatment. The most successful migrant sires were the most resident like in their flowering time, further biasing the genetic admixture between resident and migrant populations. Other loci contributing to local adaptation-heat-tolerance genes, for instance-may be in linkage disequilibrium with phenology when residents and migrants are combined into a single mating pool. Thus, introgression of potentially adaptive migrant alleles into the resident genetic background is slowed when selection acts against migrant phenology. Successful AGF programs may require sustained high immigration rates or preliminary breeding programs when phenologically matched migrant source populations are unavailable.

Identification of Host Fruit Volatiles from Snowberry (Symphoricarpos albus), Attractive to Rhagoletis zephyria Flies from the Western United States

A mixture of behaviorally active volatiles was identified from the fruit of snowberry, Symphoricarpos albus laevigatus, for Rhagoletis zephyria flies reared from snowberry fruit. A nine-component blend containing 3-methylbutan-1-ol (3%), dimethyl trisulfide (1%), 1-octen-3-ol (40%), myrcene (8%), nonanal (9%), linalool (13%), (3E)-4,8-dimethyl-1,3,7-nonatriene (DMNT, 6%), decanal (15%), and β-caryophyllene (5%) was identified that gave consistent electroantennogram activity and was behaviorally active in flight tunnel tests. In other flight tunnel assays, snowberry flies from two sites in Washington state, USA, displayed significantly greater levels of upwind oriented flight to sources with the snowberry volatile blend compared with previously identified volatile blends from domestic apple (Malus domestica) and downy hawthorn (Crataegus mollis) fruit from the eastern USA, and domestic apple, black hawthorn (C. douglasii) and ornamental hawthorn (C. monogyna) from Washington state. Selected subtraction assays showed that whereas removal of DMNT or 1-octen-3-ol significantly reduced the level of upwind flight, removal of myrcene and β-caryophyllene, or dimethyl trisulfide alone did not significantly affect the proportion of upwind flights. Our findings add to previous studies showing that populations of Rhagoletis flies infesting different host fruit are attracted to unique mixtures of volatile compounds specific to their respective host plants. Taken together, the results support the hypothesis that differences among flies in their behavioral responses to host fruit odors represent key adaptations involved in sympatric host plant shifts, contributing to host specific mating and generating prezygotic reproductive isolation among members of the R. pomonella sibling species complex.

Rapid postglacial diversification and long-term stasis within the songbird genus Junco: phylogeographic and phylogenomic evidence

Natural systems composed of closely related taxa that vary in the degree of phenotypic divergence and geographic isolation provide an opportunity to investigate the rate of phenotypic diversification and the relative roles of selection and drift in driving lineage formation. The genus Junco (Aves: Emberizidae) of North America includes parapatric northern forms that are markedly divergent in plumage pattern and colour, in contrast to geographically isolated southern populations in remote areas that show moderate phenotypic divergence. Here, we quantify patterns of phenotypic divergence in morphology and plumage colour and use mitochondrial DNA genes, a nuclear intron, and genomewide SNPs to reconstruct the demographic and evolutionary history of the genus to infer relative rates of evolutionary divergence among lineages. We found that geographically isolated populations have evolved independently for hundreds of thousands of years despite little differentiation in phenotype, in sharp contrast to phenotypically diverse northern forms, which have diversified within the last few thousand years as a result of the rapid postglacial recolonization of North America. SNP data resolved young northern lineages into reciprocally monophyletic lineages, indicating low rates of gene flow even among closely related parapatric forms, and suggesting a role for strong genetic drift or multifarious selection acting on multiple loci in driving lineage divergence. Juncos represent a compelling example of speciation in action, where the combined effects of historical and selective factors have produced one of the fastest cases of speciation known in vertebrates.

The effect of winter length on survival and duration of dormancy of four sympatric species of Rhagoletis exploiting plants with different fruiting phenology

Dormancy has been thoroughly studied for several species of economic importance in the genus Rhagoletis in temperate areas of North America and Europe. Much less is known on life history regulation for species inhabiting high-elevation areas in the subtropics at the southern extreme of their geographical range. Host plant phenology has been found to play a key role in generating allochronic isolation among sibling species and host races of Rhagoletis in the course of sympatric speciation, and has important implications for pest management. We compare the effect of winter length on survival to adult eclosion and dormancy duration among four species of Rhagoletis (three of them sympatric) exploiting hosts with different fruiting phenology in subtropical isolated highlands. Survival and duration of dormancy was found to be different among the four species. At 24°C, a very small proportion (<1%) of R. pomonella, R. turpiniae and R. zoqui completed development without becoming dormant, while in the case of R. solanophaga the majority of the population emerged after development within 40 days of pupation. Also, a large proportion of braconid parasitoids infesting Rhagoletis eggs and larvae emerged as adults without becoming dormant. Greatest survival after artificial winter was obtained for R. pomonella (50-60%) and R. zoqui (30%) after only four weeks at 5°C (a third of the time reported for studies on northern R. pomonella), while R. turpiniae, under identical environmental conditions experienced low adult emergence, and highest survival (11%) was recorded for flies exposed to 5°C during 10 and 12 weeks. For R. pomonella, there was a strong positive relationship between winter length and time to post-winter adult eclosion that was not observed for R. zoqui. In sum, for R. pomonella, mild winters in highland subtropical areas appear to select for flies better able to withstand longer periods of warm temperature before winter than flies exploiting late fruiting hosts and inhabiting northern latitudes. In the case of R. turpiniae and R. zoqui environmental cues such as fluctuations in humidity and/or different temperature thresholds (5°C) may play a more important role than winter length in life history regulation. Continuous host availability for R. solanophaga appears to have selected for non-diapausing flies. From an applied perspective our results are useful for handling flies in the laboratory to conduct research and suggest that non-diapausing strains of flies and parasitoids may be selected for SIT and innundative biological control programs.

Population structure of a vector-borne plant parasite

Parasites are among the most diverse groups of life on Earth, yet complex natural histories often preclude studies of their speciation processes. The biology of parasitic plants facilitates in situ collection of data on both genetic structure and the mechanisms responsible for that structure. Here, we studied the role of mating, dispersal and establishment in host race formation of a parasitic plant. We investigated the population genetics of a vector-borne desert mistletoe (Phoradendron californicum) across two legume host tree species (Senegalia greggii and Prosopis velutina) in the Sonoran desert using microsatellites. Consistent with host race formation, we found strong host-associated genetic structure in sympatry, little genetic variation due to geographic site and weak isolation by distance. We hypothesize that genetic differentiation results from differences in the timing of mistletoe flowering by host species, as we found initial flowering date of individual mistletoes correlated with genetic ancestry. Hybrids with intermediate ancestry were detected genetically. Individuals likely resulting from recent, successful establishment events following dispersal between the host species were detected at frequencies similar to hybrids between host races. Therefore, barriers to gene flow between the host races may have been stronger at mating than at dispersal. We also found higher inbreeding and within-host individual relatedness values for mistletoes on the more rare and isolated host species (S. greggii). Our study spanned spatial scales to address how interactions with both vectors and hosts influence parasitic plant structure with implications for parasite virulence evolution and speciation.