Divergent host-plant use promotes reproductive isolation among cynipid gall wasp populations

Ecological speciation occurs when reproductive isolation evolves as a consequence of divergent natural selection among environments. A direct prediction of this process is that ecologically divergent pairs of populations will exhibit greater reproductive isolation than ecologically similar pairs of populations. By comparing allopatric populations of the cynipid gall wasp Belonocnema treatae infesting Quercus virginiana and Quercus geminata, we tested the role that divergent host use plays in generating ecological divergence and sexual isolation. We found differences in body size and gall structure associated with divergent host use, but no difference in neutral genetic divergence between populations on the same or different host plant. We observed significant assortative mating between populations from alternative host plants but not between allopatric populations on the same host plant. Thus, we provide evidence that divergent host use promotes speciation among gall wasp populations.     

Partitioning of herbivore hosts across time and food plants promotes diversification in the Megastigmus dorsalis oak gall parasitoid complex

Communities of insect herbivores and their natural enemies are rich and ecologically crucial components of terrestrial biodiversity. Understanding the processes that promote their origin and maintenance is thus of considerable interest. One major proposed mechanism is ecological speciation through host‐associated differentiation (HAD), the divergence of a polyphagous species first into ecological host races and eventually into more specialized daughter species. The rich chalcid parasitoid communities attacking cynipid oak gall wasp hosts are structured by multiple host traits, including food plant taxon, host gall phenology, and gall structure. Here, we ask whether the same traits structure genetic diversity within supposedly generalist parasitoid morphospecies. We use mitochondrial DNA sequences and microsatellite genotypes to quantify HAD for Megastigmus (Bootanomyia) dorsalis, a complex of two apparently generalist cryptic parasitoid species attacking oak galls. Ancient Balkan refugial populations showed phenological separation between the cryptic species, one primarily attacking spring galls, and the other mainly attacking autumn galls. The spring species also contained host races specializing on galls developing on different host‐plant lineages (sections Cerris vs. Quercus) within the oak genus Quercus. These results indicate more significant host‐associated structuring within oak gall parasitoid communities than previously thought and support ecological theory predicting the evolution of specialist lineages within generalist parasitoids. In contrast, UK populations of the autumn cryptic species associated with both native and recently invading oak gall wasps showed no evidence of population differentiation, implying rapid recruitment of native parasitoid populations onto invading hosts, and hence potential for natural biological control. This is of significance given recent rapid range expansion of the economically damaging chestnut gall wasp, Dryocosmus kuriphilus, in Europe.     

Evolution of intrinsic reproductive isolation among four North American populations of Rhagoletis pomonella (Diptera: Tephritidae)

Across its range in North America, four geographically separated, ecologically and genetically diverged populations of hawthorn (Crataegus)‐infesting Rhagoletis pomonella (Diptera: Tephritidae) flies inhabit the Eje Volcánico Trans Mexicano (EVTM), the Sierra Madre Oriental (SMO), the Chiapas Highlands (CHIS) and the USA. Here, we tested whether these four populations are reproductively isolated by any intrinsic, nonhost‐related, pre‐ or postmating barriers to gene flow. Crossing experiments suggested that a low level of host‐independent prezygotic isolation may exist between hawthorn flies from EVTM and the three other populations, but only with respect to a slight reduction in copulation duration in EVTM matings. Some evidence for postmating isolation was found, again primarily involving EVTM crossed to SMO, CHIS and US flies. Certain crosses produced no (SMO male × EVTM female) or few (EVTM male × CHIS female; CHIS male × SMO female) F1 hybrid offspring. F2 crosses were generally fertile, except for US male × CHIS female matings. Inherent reproductive isolation therefore appears to be quantitative rather than absolute between populations, as the possibility for gene flow exists through at least some combinations of mating among EVTM, SMO, CHIS and US flies. Our results are consistent with a recently advanced hypothesis that episodic introgression from Mexico into the USA has played a role in providing genetic variation, facilitating sympatric host race formation and the adaptive radiation of the R. pomonella sibling species' complex in the USA. © 2010 The Linnean Society of London, Biological Journal of the Linnean Society, 2010, 100 , 213–223.     

Ecological adaptation and reproductive isolation in sympatry: genetic and phenotypic evidence for native host races of Rhagoletis pomonella

Ecological speciation with gene flow may be an important mode of diversification for phytophagous insects. The recent shift of Rhagoletis pomonella from its native host downy hawthorn (Crataegus mollis) to introduced apple (Malus domestica) in the northeastern United States is a classic example of sympatric host race formation. Here, we test whether R. pomonella has similarly formed host races on four native Crataegus species in the southern United States: western mayhaw (C. opaca), blueberry hawthorn (C. brachyacantha), southern red hawthorn (C. mollis var. texana) and green hawthorn (C. viridis). These four southern hosts differ from each other in their fruiting phenology and in the volatile compounds emitted from the surface of their fruits. These two traits form the basis of ecological reproductive isolation between downy hawthorn and apple flies in the north. We report evidence from microsatellite population surveys and eclosion studies supporting the existence of genetically differentiated and partially reproductively isolated host races of southern hawthorn flies. The results provide an example of host shifting and ecological divergence involving native plants and imply that speciation with gene flow may be commonly initiated in Rhagoletis when ecological opportunity presents itself.     

Testing for host‐associated differentiation in two egg parasitoids of a forest herbivore

Several studies underline the importance of ecological barriers and differential selection in driving sympatric speciation. Host‐associated differentiation (HAD) has been proposed as one of the mechanisms leading to sympatric speciation. However, it is still unclear how common HAD is or which are the factors that could promote it. In particular, not much is known about HAD in predators and parasitoids of herbivorous insects. One of the characteristics postulated to pre‐dispose insects to HAD is parthenogenesis as it may favour adaptive responses to particular environments, amplifying selected gene complexes. In this study, we used amplified fragment length polymorphism (AFLP) markers to determine whether HAD is present in two parthenogenetic egg parasitoids attacking the same herbivore species – the pine processionary moth, Thaumetopoea pityocampa (Denis & Schiffermüller) (Lepidoptera: Notodontidae) – on two host Pinus species. A total of 100 loci for 59 individuals sampled in four populations of Baryscapus servadeii (Domenichini) (Hymenoptera: Eulophidae), a specialist parasitoid, and 106 loci for 117 individuals sampled in six populations of Ooencyrtus pityocampae Mercet (Hymenoptera: Encyrtidae), a generalist parasitoid, were analysed. Levels of genetic differentiation were also assessed with an outlier analysis, checking for alleles associated to host plants. No evidence of HAD was detected in any of the two parasitoid species. We hypothesize that both the lack of strict parthenogenetic reproduction and the ectophagous nature of the insect host could explain the absence of HAD. The genetic variation observed in the generalist parasitoid responded to a pattern of local adaptation, whereas no relationship with either host or geography was found in the specialist parasitoid.     

Seedling response to water stress in valley oak (Quercus lobata) is shaped by different gene networks across populations

Drought is a major stress for plants, creating a strong selection pressure for traits that enable plant growth and survival in dry environments. Many drought responses are conserved species‐wide responses, while others vary among populations distributed across heterogeneous environments. We tested how six populations of the widely distributed California valley oak (Quercus lobata) sampled from contrasting climates would differ in their response to soil drying relative to well‐watered controls in a common environment by measuring ecophysiological traits in 93 individuals and gene expression (RNA‐seq) in 42 individuals. Populations did not differ in their adjustment of turgor loss point during soil drying, suggesting a generalized species‐wide response. Differential expression analysis identified 689 genes with a common response to treatment across populations and 470 genes with population‐specific responses. Weighted gene co‐expression network analysis (WGCNA) identified groups of genes with similar expression patterns that may be regulated together (gene modules). Several gene modules responded differently to water stress among populations, suggesting regional differences in gene network regulation. Populations from sites with a high mean annual temperature responded to the imposed water stress with significantly greater changes in gene module expression, indicating that these populations may be locally adapted to respond to drought. We propose that this variation among valley oak populations provides a mechanism for differential tolerance to the increasingly frequent and severe droughts in California.     

The accumulation of reproductive isolation in early stages of divergence supports a role for sexual selection

Models of speciation by sexual selection propose that male–female coevolution leads to the rapid evolution of behavioural reproductive isolation. Here, we compare the strength of behavioural isolation to ecological isolation, gametic incompatibility and hybrid inviability in a group of dichromatic stream fishes. In addition, we examine whether any of these individual barriers, or a combined measure of total isolation, is predicted by body shape differences, male colour differences, environmental differences or genetic distance. Behavioural isolation reaches the highest values of any barrier and is significantly greater than ecological isolation. No individual reproductive barrier is associated with any of the predictor variables. However, marginally significant relationships between male colour and body shape differences with ecological and behavioural isolation are discussed. Differences in male colour and body shape predict total reproductive isolation between species; hierarchical partitioning of these two variables' effects suggests a stronger role for male colour differences. Together, these results suggest an important role for divergent sexual selection in darter speciation but raise new questions about the mechanisms of sexual selection at play and the role of male nuptial ornaments.     

Genetic structure among greater white‐fronted goose populations of the Pacific Flyway

An understanding of the genetic structure of populations in the wild is essential for long‐term conservation and stewardship in the face of environmental change. Knowledge of the present‐day distribution of genetic lineages (phylogeography) of a species is especially important for organisms that are exploited or utilize habitats that may be jeopardized by human intervention, including climate change. Here, we describe mitochondrial (mtDNA) and nuclear genetic (microsatellite) diversity among three populations of a migratory bird, the greater white‐fronted goose (Anser albifrons), which breeds discontinuously in western and southwestern Alaska and winters in the Pacific Flyway of North America. Significant genetic structure was evident at both marker types. All three populations were differentiated for mtDNA, whereas microsatellite analysis only differentiated geese from the Cook Inlet Basin. In sexual reproducing species, nonrandom mate selection, when occurring in concert with fine‐scale resource partitioning, can lead to phenotypic and genetic divergence as we observed in our study. If mate selection does not occur at the time of reproduction, which is not uncommon in long‐lived organisms, then mechanisms influencing the true availability of potential mates may be obscured, and the degree of genetic and phenotypic diversity may appear incongruous with presumed patterns of gene flow. Previous investigations revealed population‐specific behavioral, temporal, and spatial mechanisms that likely influence the amount of gene flow measured among greater white‐fronted goose populations. The degree of observed genetic structuring aligns well with our current understanding of population differences pertaining to seasonal movements, social structure, pairing behavior, and resource partitioning.     

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

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

Contrasting population genetic structure among freshwater‐resident and anadromous lampreys: the role of demographic history,differential dispersal and anthropogenic barriers to movement

The tendency of many species to abandon migration remains a poorly understood aspect of evolutionary biology that may play an important role in promoting species radiation by both allopatric and sympatric mechanisms. Anadromy inherently offers an opportunity for the colonization of freshwater environments, and the shift from an anadromous to a wholly freshwater life history has occurred in many families of fishes. Freshwater‐resident forms have arisen repeatedly among lampreys (within the Petromyzontidae and Mordaciidae), and there has been much debate as to whether anadromous lampreys, and their derived freshwater‐resident analogues, constitute distinct species or are divergent ecotypes of polymorphic species. Samples of 543 European river lamprey Lampetra fluviatilis (mostly from anadromous populations) and freshwater European brook lamprey Lampetra planeri from across 18 sites, primarily in the British Isles, were investigated for 13 polymorphic microsatellite DNA loci, and 108 samples from six of these sites were sequenced for 829 bp of mitochondrial DNA (mtDNA). We found contrasting patterns of population structure for mtDNA and microsatellite DNA markers, such that low diversity and little structure were seen for all populations for mtDNA (consistent with a recent founder expansion event), while fine‐scale structuring was evident for nuclear markers. Strong differentiation for microsatellite DNA loci was seen among freshwater‐resident L. planeri populations and between L. fluviatilis and L. planeri in most cases, but little structure was evident among anadromous L. fluviatilis populations. We conclude that postglacial colonization founded multiple freshwater‐resident populations with strong habitat fidelity and limited dispersal tendencies that became highly differentiated, a pattern that was likely intensified by anthropogenic barriers.     

The genomic bases of morphological divergence and reproductive isolation driven by ecological speciation in Senecio (Asteraceae)

Ecological speciation, driven by adaptation to contrasting environments, provides an attractive opportunity to study the formation of distinct species, and the role of selection and genomic divergence in this process. Here, we focus on a particularly clear‐cut case of ecological speciation to reveal the genomic bases of reproductive isolation and morphological differences between closely related Senecio species, whose recent divergence within the last ~200 000 years was likely driven by the uplift of Mt. Etna (Sicily). These species form a hybrid zone, yet remain morphologically and ecologically distinct, despite active gene exchange. Here, we report a high‐density genetic map of the Senecio genome and map hybrid breakdown to one large and several small quantitative trait loci (QTL). Loci under diversifying selection cluster in three 5 cM regions which are characterized by a significant increase in relative (F_ST), but not absolute (d_XY), interspecific differentiation. They also correspond to some of the regions of greatest marker density, possibly corresponding to ‘cold‐spots’ of recombination, such as centromeres or chromosomal inversions. Morphological QTL for leaf and floral traits overlap these clusters. We also detected three genomic regions with significant transmission ratio distortion (TRD), possibly indicating accumulation of intrinsic genetic incompatibilities between these recently diverged species. One of the TRD regions overlapped with a cluster of high species differentiation, and another overlaps the large QTL for hybrid breakdown, indicating that divergence of these species may have occurred due to a complex interplay of ecological divergence and accumulation of intrinsic genetic incompatibilities.     

Using genomic data to revisit an early example of reproductive character displacement in Haitian Anolis lizards

The pattern of reproductive character displacement (RCD)—in which traits associated with reproductive isolation are more different where two species occur together than where they occur in isolation—is frequently attributed to reinforcement, a process during which natural selection acting against maladaptive mating events leads to enhanced prezygotic isolation between species or incipient species. One of the first studies of RCD to include molecular genetic data was described 40 years ago in a complex of Haitian trunk anole lizards using a small number of allozyme loci. In this example, Anolis caudalis appears to experience divergence in the color and pattern of an extensible throat fan, or dewlap, in areas of contact with closely related species at the northern and southern limits of its range. However, this case study has been largely overlooked for decades; meanwhile, explanations for geographic variation in dewlap color and pattern have focused primarily on adaptation to local signalling environments. We reinvestigate this example using amplified fragment length polymorphism (AFLP) genome scans, mtDNA sequence data, information on dewlap phenotypes and GIS data on environmental variation to test the hypothesis of RCD generated by reinforcement in Haitian trunk anoles. Together, our phenotypic and genetic results are consistent with RCD at the southern and northern limits of the range of A. caudalis. We evaluate the evidence for reinforcement as the explanation for RCD in Haitian trunk anoles, consider alternative explanations and provide suggestions for future work on the relationship between dewlap variation and speciation in Haitian trunk anoles.     

Dietary specialization in mutualistic acacia‐ants affects relative abundance but not identity of host‐associated bacteria

Acacia‐ant mutualists in the genus Pseudomyrmex nest obligately in acacia plants and, as we show through stable isotope analysis, feed at a remarkably low trophic level. Insects with diets such as these sometimes depend on bacterial symbionts for nutritional enrichment. We, therefore, examine the bacterial communities associated with acacia‐ants in order to determine whether they host bacterial partners likely to contribute to their nutrition. Despite large differences in trophic position, acacia‐ants and related species with generalized diets do not host distinct bacterial taxa. However, we find that a small number of previously undescribed bacterial taxa do differ in relative abundance between acacia‐ants and generalists, including several Acetobacteraceae and Nocardiaceae lineages related to common insect associates. Comparisons with an herbivorous generalist, a parasite that feeds on acacias and a mutualistic species with a generalized diet show that trophic level is likely responsible for these small differences in bacterial community structure. While we did not experimentally test for a nutritional benefit to hosts of these bacterial lineages, metagenomic analysis reveals a Bartonella relative with an intact nitrogen‐recycling pathway widespread across Pseudomyrmex mutualists and generalists. This taxon may be contributing to nitrogen enrichment of its ant hosts through urease activity and, concordant with an obligately host‐associated lifestyle, appears to be experiencing genomewide relaxed selection. The lack of distinctiveness in bacterial communities across trophic level in this group of ants shows a remarkable ability to adjust to varied diets, possibly with assistance from these diverse ant‐specific bacterial lineages.     

Microgeographic morphological variation across larval wood frog populations associated with environment despite gene flow

Gene flow has historically been thought to constrain local adaptation; yet, recent research suggests that populations can diverge despite exchanging genes. Here I use a common garden experiment to assess the combined effects of gene flow and natural selection on morphological variation of 16 wood frog (Rana sylvatica) populations, a species known to experience divergent selection pressures in open‐ and closed‐canopy ponds across relatively small geographic scales. Wood frog tadpoles from different ponds showed significant morphological variation associated with canopy type with a trade‐off between tail length and body depth consistent with previous research. In contrast, neutral genetic differentiation of nine microsatellite loci as measured by Jost's D was not associated with canopy type, indicating no pattern of isolation by environment. Genetic structure analyses indicated some substructure across the 16 ponds (K = 4); however, three out of four assigned clusters included both open‐ and closed‐canopy ponds. Together, these results suggest that morphological divergence among these wood frog populations is occurring despite gene flow and that selection within these environments is strong. Furthermore, morphological variation among ponds differed across two sampling periods during larval development, demonstrating the importance of evaluating phenotypic divergence over multiple time periods and at a time relevant to the processes being studied.     

Effects of quantitative and qualitative differences in volatiles from host‐ and non‐host‐infested maize on the attraction of the larval parasitoid Cotesia kariyai

Cotesia kariyai Watanabe (Hymenoptera: Braconidae) is a specialist larval parasitoid of Mythimna separata Walker (Lepidoptera: Noctuidae). Cotesia kariyai wasps use herbivore‐induced plant volatiles (HIPVs) to locate hosts. However, complex natural habitats are full of volatiles released by both herbivorous host‐ and non‐host‐infested plants at various levels of intensity. Therefore, the presence of non‐hosts may affect parasitoid decisions while foraging. Here, the host‐finding efficiency of naive C. kariyai from HIPVs influenced by host‐ and non‐host‐infested maize [Zea mays L. (Poaceae)] plants was investigated with a four‐arm olfactometer. Ostrinia furnacalis Guenée (Lepidoptera: Crambidae) was selected as a non‐host species. One unit (1 U) of host‐ or non‐host‐infested plant was prepared by infesting a potted plant with five host or seven non‐host larvae. In two‐choice bioassays, host‐infested plants fed upon by different numbers of larvae, and various units of host‐ and non‐host‐infested plants (infestation units; 1 U, 2 U, and 3 U) were arranged to examine the effects of differences in volatile quantity and quality on the olfactory responses of C. kariyai with the assumption that volatile quantity and quality changes with differences in numbers of insects and plants. Cotesia kariyai was found to perceive quantitative differences in volatiles from host‐infested plants, preferring larger quantities of volatiles from larger numbers of larvae or plants. Also, the parasitoids discriminated between healthy plants, host‐infested plants, and non‐host‐infested plants by recognising volatiles released from those plants. Cotesia kariyai showed a reduced preference for host‐induced volatiles, when larger numbers of non‐host‐infested plants were present. Therefore, quantitative and qualitative differences in volatiles from host‐ and non‐host‐infested plants appear to affect the decision of C. kariyai during host‐habitat searching in multiple tritrophic systems.     

Cascading reproductive isolation: Plant phenology drives temporal isolation among populations of a host‐specific herbivore

All organisms exist within a complex network of interacting species, thus evolutionary change may have reciprocal effects on multiple taxa. Here, we demonstrate “cascading reproductive isolation,” whereby ecological differences that reduce gene flow between populations at one trophic level affect reproductive isolation (RI) among interacting species at the next trophic level. Using a combination of field, laboratory and common‐garden studies and long‐term herbaria records, we estimate and evaluate the relative contribution of temporal RI to overall prezygotic RI between populations of Belonocnema treatae, a specialist gall‐forming wasp adapted to sister species of live oak (Quercus virginiana and Q. geminata). We link strong temporal RI between host‐associated insect populations to differences between host plant budbreak phenology. Budbreak initiates flowering and the production of new leaves, which are an ephemeral resource critical to insect reproduction. As flowering time is implicated in RI between plant species, budbreak acts as a “multitrophic multi‐effect trait,” whereby differences in budbreak phenology contribute to RI in plants and insects. These sister oak species share a diverse community of host‐specific gall‐formers and insect natural enemies similarly dependent on ephemeral plant tissues. Thus, our results set the stage for testing for parallelism in a role of plant phenology in driving temporal cascading RI across multiple species and trophic levels.     

Impact of pre‐ and post‐pollination barriers on pollen transfer and reproductive isolation among three sympatric Pedicularis (Orobanchaceae) species

• Pedicularis is the largest genus in the Orobanchaceae (>300) with many species co‐occurring and co‐blooming in subalpine to alpine meadows in the Himalayas. Although it is well known that different Pedicularis species place pollen on different parts of the same bumblebee's body, thus reducing interspecific pollen transfer, it is not known whether post‐pollination components also contribute to reproductive isolation (RI).
• In this study, we quantified the individual strengths and absolute contributions of six pre‐ and post‐pollination components of RI between three sympatric species in two pairs; Pedicularis gruina × Pedicularis tenuisecta (gru × ten) and Pedicularis comptoniifolia × Pedicularis tenuisecta (com × ten).
• All three Pedicularis species shared the same Bombus species. Individual foragers showed a high, but incomplete, floral constancy for each species. Therefore, pre‐pollination barriers were potentially ‘leaky’ as Bombus species showed a low but consistent frequency of interspecific visitation. The RI strength of pre‐pollination was lower in com × ten than in gru × ten. In contrast, post‐pollination barriers completely blocked gene flow between both sets of species pairs. Two post‐pollination recognition sites were identified. Late acting rejection of interspecific pollen tube growth occurred in com♀ × ten♂, while seeds produced in bi‐directional crosses of gru × ten failed to germinate.
• We propose that, although floral isolation based on pollen placement on pollinators in the genus Pedicularis is crucial to avoid interspecific pollen transfer, the importance of this mode of interspecific isolation may be exaggerated. Post‐pollination barriers may play even larger roles for currently established populations of co‐blooming and sympatric species in this huge genus in the Himalayas.

     

Female reproductive morphology of coral‐inhabiting gall crabs (Crustacea: Decapoda: Brachyura: Cryptochiridae)

Gall crabs are obligate associates of stony corals in which they induce skeletal modifications. In some cryptochirid species, females live in open depressions accessible to males; while in others, females are rather isolated in semiclosed galls, which necessitates elaborate sperm storage capabilities by the female. In this study we investigate the female gross morphology and reproductive systems of Fungicola syzygia lodged in semiclosed flattened pits, Opecarcinus cathyae with semiopen pits and Pseudocryptochirus viridis from shallow open depressions using line drawings and histological methods. The general morphology of the cryptochirids' reproductive systems is uniform and conforms to other thoracotreme brachyurans: paired muscular vaginae of the concave pattern lead from the sternal gonopores into paired seminal receptacles where sperm is stored. The seminal receptacle is internally lined by distinct types of epithelia: a cuticle underlined by a columnar epithelium ventrally, a monolayered secretory epithelium dorsally and a multilayered transfer tissue where the oviducts enter the seminal receptacle. In all studied specimens, the seminal receptacle contained free spermatozoa; however, in specimens of Pseudocryptochirus viridis it also contained spermatophores, indicating a recent insemination. In contrast to most other brachyurans ovaries of the investigated cryptochirids extend into the pleon. The specific degree of ovary extension differs between the studied species and is closely related to female body shape.     

Experimental evidence of genome‐wide impact of ecological selection during early stages of speciation‐with‐gene‐flow

Theory predicts that speciation‐with‐gene‐flow is more likely when the consequences of selection for population divergence transitions from mainly direct effects of selection acting on individual genes to a collective property of all selected genes in the genome. Thus, understanding the direct impacts of ecologically based selection, as well as the indirect effects due to correlations among loci, is critical to understanding speciation. Here, we measure the genome‐wide impacts of host‐associated selection between hawthorn and apple host races of Rhagoletis pomonella (Diptera: Tephritidae), a model for contemporary speciation‐with‐gene‐flow. Allele frequency shifts of 32 455 SNPs induced in a selection experiment based on host phenology were genome wide and highly concordant with genetic divergence between co‐occurring apple and hawthorn flies in nature. This striking genome‐wide similarity between experimental and natural populations of R. pomonella underscores the importance of ecological selection at early stages of divergence and calls for further integration of studies of eco‐evolutionary dynamics and genome divergence.