Divergence hitchhiking and the spread of genomic isolation during ecological speciation-with-gene-flow

In allopatric populations, geographical separation simultaneously isolates the entire genome, allowing genetic divergence to accumulate virtually anywhere in the genome. In sympatric populations, however, the strong divergent selection required to overcome migration produces a genetic mosaic of divergent and non-divergent genomic regions. In some recent genome scans, each divergent genomic region has been interpreted as an independent incidence of migration/selection balance, such that the reduction of gene exchange is restricted to a few kilobases around each divergently selected gene. I propose an alternative mechanism, 'divergence hitchhiking' (DH), in which divergent selection can reduce gene exchange for several megabases around a gene under strong divergent selection. Not all genes/markers within a DH region are divergently selected, yet the entire region is protected to some degree from gene exchange, permitting genetic divergence from mechanisms other than divergent selection to accumulate secondarily. After contrasting DH and multilocus migration/selection balance (MM/SB), I outline a model in which genomic isolation at a given genomic location is jointly determined by DH and genome-wide effects of the progressive reduction in realized migration, then illustrate DH using data from several pairs of incipient species in the wild.     

The genomic consequences of adaptive divergence and reproductive isolation between species of manakins

The processes of adaptation and speciation are expected to shape genomic variation within and between diverging species. Here we analyze genomic heterogeneity of genetic differentiation and introgression in a hybrid zone between two bird species (Manacus candei and M. vitellinus) using 59 100 SNPs, a whole genome assembly, and Bayesian models. Measures of genetic differentiation () and introgression (genomic cline center [α] and rate [β]) were highly heterogeneous among loci. We identified thousands of loci with elevated parameter estimates, some of which are likely to be associated with variation in fitness in Manacus populations. To analyze the genomic organization of differentiation and introgression, we mapped SNPs onto a draft assembly of the M. vitellinus genome. Estimates of , α, and β were autocorrelated at very short physical distances (< 100 bp), but much less so beyond this. In addition, average statistical associations (linkage disequilibrium) between SNPs were generally low and were not higher in admixed populations than in populations of the parental species. Although they did not occur with a constant probability across the genome, loci with elevated , α, and β were not strongly co-localized in the genome. Contrary to verbal models that predict clustering of loci involved in adaptation and isolation in discrete genomic regions, these results are consistent with the hypothesis that genetic regions involved in adaptive divergence and reproductive isolation are scattered throughout the genome. We also found that many loci were characterized by both exceptional genetic differentiation and introgression, consistent with the hypothesis that loci involved in isolation are also often characterized by a history of divergent selection. However, the concordance between isolation and differentiation was only partial, indicating a complex architecture and history of loci involved in isolation.     

Molecular genetic and quantitative trait divergence associated with recent homoploid hybrid speciation: a study of Senecio squalidus (Asteraceae)

Hybridization is increasingly seen as a trigger for rapid evolution and speciation. To quantify and qualify divergence associated with recent homoploid hybrid speciation, we compared quantitative trait (QT) and molecular genetic variation between the homoploid hybrid species Senecio squalidus and its parental species, S. aethnensis and S. chrysanthemifolius, and also their naturally occurring Sicilian hybrids. S. squalidus originated and became invasive in the United Kingdom following the introduction of hybrid plants from Mount Etna, Sicily, about 300 years ago. We recorded considerable molecular genetic differentiation between S. squalidus and its parents and their Sicilian hybrids in terms of both reduced genetic diversity and altered allele frequencies, potentially due to the genetic bottleneck associated with introduction to the United Kingdom. S. squalidus is also distinct from its parents and Sicilian hybrids for QTs, but less so than for molecular genetic markers. We suggest that this is due to resilience of polygenic QTs to changes in allele frequency or lack of selection for hybrid niche divergence in geographic isolation. While S. squalidus is intermediate or parental-like for most QTs, some trangressively distinct traits were observed, which might indicate emerging local adaptation in its invasive range. This study emphasizes the important contribution of founder events and geographic isolation to successful homoploid hybrid speciation.     

Growth form evolution and hybridization in Senecio (Asteraceae) from the high equatorial Andes

Changes in growth forms frequently accompany plant adaptive radiations, including páramo–a high‐elevation treeless habitat type of the northern Andes. We tested whether diverse group of Senecio inhabiting montane forests and páramo represented such growth form changes. We also investigated the role of Andean geography and environment in structuring genetic variation of this group. We sampled 108 populations and 28 species of Senecio (focusing on species from former genera Lasiocephalus and Culcitium) and analyzed their genetic relationships and patterns of intraspecific variation using DNA fingerprinting (AFLPs) and nuclear DNA sequences (ITS). We partitioned genetic variation into environmental and geographical components. ITS‐based phylogeny supported monophyly of a Lasiocephalus‐Culcitium clade. A grade of herbaceous alpine Senecio species subtended the Lasiocephalus‐Culcitium clade suggesting a change from the herbaceous to the woody growth form. Both ITS sequences and the AFLPs separated a group composed of the majority of páramo subshrubs from other group(s) comprising both forest and páramo species of various growth forms. These morphologically variable group(s) further split into clades encompassing both the páramo subshrubs and forest lianas, indicating independent switches among the growth forms and habitats. The finest AFLP genetic structure corresponded to morphologically delimited species except in two independent cases in which patterns of genetic variation instead reflected geography. Several morphologically variable species were genetically admixed, which suggests possible hybrid origins. Latitude and longitude accounted for 5%–8% of genetic variation in each of three AFLP groups, while the proportion of variation attributed to environment varied between 8% and 31% among them. A change from the herbaceous to the woody growth form is suggested for species of high‐elevation Andean Senecio. Independent switches between habitats and growth forms likely occurred within the group. Hybridization likely played an important role in species diversification.     

Integrating Bayesian genomic cline analyses and association mapping of morphological and ecological traits to dissect reproductive isolation and introgression in a Louisiana Iris hybrid zone

Hybrid zones provide unique opportunities to examine reproductive isolation and introgression in nature. We utilized 45,384 single nucleotide polymorphism (SNP) loci to perform association mapping of 14 floral, vegetative and ecological traits that differ between Iris hexagona and Iris fulva, and to investigate, using a Bayesian genomic cline (BGC) framework, patterns of genomic introgression in a large and phenotypically diverse hybrid zone in southern Louisiana. Many loci of small effect size were consistently found to be associated with phenotypic variation across all traits, and several individual loci were revealed to influence phenotypic variation across multiple traits. Patterns of genomic introgression were quite heterogeneous throughout the Louisiana Iris genome, with I. hexagona alleles tending to be favoured over those of I. fulva. Loci that were found to have exceptional patterns of introgression were also found to be significantly associated with phenotypic variation in a small number of morphological traits. However, this was the exception rather than the rule, as most loci that were associated with morphological trait variation were not significantly associated with excess ancestry. These findings provide insights into the complexity of the genomic architecture of phenotypic differences and are a first step towards identifying loci that are associated with both trait variation and reproductive isolation in nature.     

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.     

Change in sexual signalling traits outruns morphological divergence across an ecological gradient in the post‐glacial radiation of the songbird genus Junco

The relative roles of natural and sexual selection in promoting evolutionary lineage divergence remains controversial and difficult to assess in natural systems. Local adaptation through natural selection is known to play a central role in promoting evolutionary divergence, yet secondary sexual traits can vary widely among species in recent radiations, suggesting that sexual selection may also be important in the early stages of speciation. Here, we compare rates of divergence in ecologically relevant traits (morphology) and sexually selected signalling traits (coloration) relative to neutral structure in genome‐wide molecular markers and examine patterns of variation in sexual dichromatism to explore the roles of natural and sexual selection in the diversification of the songbird genus Junco (Aves: Passerellidae). Juncos include divergent lineages in Central America and several dark‐eyed junco (J. hyemalis) lineages that diversified recently as the group recolonized North America following the last glacial maximum (ca. 18,000 years ago). We found an accelerated rate of divergence in sexually selected characters relative to ecologically relevant traits. Moreover, sexual dichromatism measurements suggested a positive relationship between the degree of colour divergence and the strength of sexual selection when controlling for neutral genetic distance. We also found a positive correlation between dichromatism and latitude, which coincides with the geographic axis of decreasing lineage age in juncos but also with a steep ecological gradient. Finally, we found significant associations between genome‐wide variants linked to functional genes and proxies of both sexual and natural selection. These results suggest that the joint effects of sexual and ecological selection have played a prominent role in the junco radiation.     

Major structural genomic alterations can be associated with hybrid speciation in Aegilops markgrafii (Triticeae)

During evolutionary history many grasses from the tribe Triticeae have undergone interspecific hybridization, resulting in allopolyploidy; whereas homoploid hybrid speciation was found only in rye. Homoeologous chromosomes within the Triticeae preserved cross‐species macrocolinearity, except for a few species with rearranged genomes. Aegilops markgrafii, a diploid wild relative of wheat (2n = 2x = 14), has a highly asymmetrical karyotype that is indicative of chromosome rearrangements. Molecular cytogenetics and next‐generation sequencing were used to explore the genome organization. Fluorescence in situ hybridization with a set of wheat cDNAs allowed the macrostructure and cross‐genome homoeology of the Ae. markgrafii chromosomes to be established. Two chromosomes maintained colinearity, whereas the remaining were highly rearranged as a result of inversions and inter‐ and intrachromosomal translocations. We used sets of barley and wheat orthologous gene sequences to compare discrete parts of the Ae. markgrafii genome involved in the rearrangements. Analysis of sequence identity profiles and phylogenic relationships grouped chromosome blocks into two distinct clusters. Chromosome painting revealed the distribution of transposable elements and differentiated chromosome blocks into two groups consistent with the sequence analyses. These data suggest that introgressive hybridization accompanied by gross chromosome rearrangements might have had an impact on karyotype evolution and homoploid speciation in Ae. markgrafii.     

Insight into the roles of selection in speciation from genomic patterns of divergence and introgression in secondary contact in venomous rattlesnakes

Investigating secondary contact of historically isolated lineages can provide insight into how selection and drift influence genomic divergence and admixture. Here, we studied the genomic landscape of divergence and introgression following secondary contact between lineages of the Western Diamondback Rattlesnake (Crotalus atrox) to determine whether genomic regions under selection in allopatry also contribute to reproductive isolation during introgression. We used thousands of nuclear loci to study genomic differentiation between two lineages that have experienced recent secondary contact following isolation, and incorporated sampling from a zone of secondary contact to identify loci that are resistant to gene flow in hybrids. Comparisons of patterns of divergence and introgression revealed a positive relationship between allelic differentiation and resistance to introgression across the genome, and greater‐than‐expected overlap between genes linked to lineage‐specific divergence and loci that resist introgression. Genes linked to putatively selected markers were related to prominent aspects of rattlesnake biology that differ between populations of Western Diamondback rattlesnakes (i.e., venom and reproductive phenotypes). We also found evidence for selection against introgression of genes that may contribute to cytonuclear incompatibility, consistent with previously observed biased patterns of nuclear and mitochondrial alleles suggestive of partial reproductive isolation due to cytonuclear incompatibilities. Our results provide a genome‐scale perspective on the relationships between divergence and introgression in secondary contact that is relevant for understanding the roles of selection in maintaining partial isolation of lineages, causing admixing lineages to not completely homogenize.     

Multiple forms of selection shape reproductive isolation in a primate hybrid zone

Speciation occurs when populations diverge and become reproductively isolated from each other. Natural selection is commonly accepted to play a large role in this process, and it has been widely assumed that reproductive isolation often results as a by‐product of divergence driven by adaptation in allopatry. When such populations come into secondary contact, reinforcement can act to strengthen reproductive isolation, but the frequency and importance of this process are still unknown. Here, we explored genomic signatures of selection in allopatry and sympatry for loci associated with reproductive isolation using a natural primate hybrid zone. By analysing reduced‐representation sequencing data, we quantified admixture and population structure across a howler monkey hybrid zone and examined the relationship between locus‐specific differentiation and introgression. We detected extensive admixture that was mostly limited to the narrow contact zone. Loci with reduced introgression into the heterospecific genomic background (the pattern expected for loci associated with reproductive isolation due to selection against hybrids) were significantly more differentiated between allopatric parental populations than loci with neutral and increased introgression, supporting the hypothesis that reproductive isolation is a by‐product of divergence in allopatry. Further, loci with reduced introgression showed greater differentiation in sympatry than in allopatry, suggesting a role for reinforcement. Thus, our results reflect multiple forms of selection that have shaped reproductive isolation in this system. We conclude that reproductive isolation may have initially been driven by divergence in allopatry, but later reinforced by divergent selection in sympatry.     

Life history affects the evolution of reproductive isolation among species of Coreopsis (Asteraceae)

The evolution of reproductive isolation within Coreopsis is investigated by integrating phylogenetic data with estimates of pollen viability of plants from inter- and intraspecific crosses. Three different models that predict F1 fitness are compared. The first uses ITS pairwise distances, the second is based on phylogenetic branch lengths derived from DNA sequences, and the third elaborates on the second model by dividing branch length according to reconstructions of the evolution of life history. This is the first study to use phylogenetic branch-length estimates for predicting levels of reproductive isolation. Estimated branch lengths (model 2) predict hybrid fitness more accurately than simply genetic distance (model 1) but only very slightly. This is probably because the two variables are strongly correlated in Coreopsis. Prediction is substantially improved by allowing evolutionary rates to differ between annual and perennial branches (model 3). A bootstrapping procedure indicates that the life-history effect is statistically significant. The more rapid evolution of reproductive isolation within annual species of Coreopsis may be due to differing mechanisms of reproductive isolation, that is, chromosomal rearrangements rather than genetic incompatibilities.     

Genomic divergence across ecological gradients in the Central African rainforest songbird (Andropadus virens)

The little greenbul, a common rainforest passerine from sub‐Saharan Africa, has been the subject of long‐term evolutionary studies to understand the mechanisms leading to rainforest speciation. Previous research found morphological and behavioural divergence across rainforest–savannah transition zones (ecotones), and a pattern of divergence with gene flow suggesting divergent natural selection has contributed to adaptive divergence and ecotones could be important areas for rainforests speciation. Recent advances in genomics and environmental modelling make it possible to examine patterns of genetic divergence in a more comprehensive fashion. To assess the extent to which natural selection may drive patterns of differentiation, here we investigate patterns of genomic differentiation among populations across environmental gradients and regions. We find compelling evidence that individuals form discrete genetic clusters corresponding to distinctive environmental characteristics and habitat types. Pairwise F_ST between populations in different habitats is significantly higher than within habitats, and this differentiation is greater than what is expected from geographic distance alone. Moreover, we identified 140 SNPs that showed extreme differentiation among populations through a genomewide selection scan. These outliers were significantly enriched in exonic and coding regions, suggesting their functional importance. Environmental association analysis of SNP variation indicates that several environmental variables, including temperature and elevation, play important roles in driving the pattern of genomic diversification. Results lend important new genomic evidence for environmental gradients being important in population differentiation.     

Patterns of reproductive isolation within and between two Lygaeus species characterized by sexual conflicts over mating

Theory suggests that, under some circumstances, sexual conflict over mating can lead to divergent sexually antagonistic coevolution among populations for traits associated with mating, and that this can promote reproductive isolation and hence speciation. However, sexual conflict over mating may also select for traits (e.g. male willingness to mate) that enhance gene flow between populations, limiting population divergence. In the present study, we compare pre‐ and post‐mating isolation within and between two species characterized by male–female conflict over mating rate. We quantify sexual isolation among five populations of the seed bug Lygaeus equestris collected from Italy and Sweden, and two replicates of a population of the sister‐species Lygaeus simulans, also collected from Italy. We find no evidence of reproductive isolation amongst populations of L. equestris, suggesting that sexual conflict over mating has not led to population divergence in relevant mating traits in L. equestris. However, there was strong asymmetric pre‐mating isolation between L. equestris and L. simulans: male L. simulans were able to mate successfully with female L. equestris, whereas male L. equestris were largely unable to mate with female L. simulans. We found little evidence for strong post‐mating isolation between the two species, however, with hybrid F₂ offspring being produced. Our results suggest that sexual conflict over mating has not led to population divergence, and indeed perhaps supports the contrary theoretical prediction that male willingness to mate may retard speciation by promoting gene flow.     

Maternal–foetal genomic conflict and speciation: no evidence for hybrid placental dysplasia in crosses between two house mouse subspecies

Interspecific hybridization between closely related mammalian species, including various species of the genus Mus, is commonly associated with abnormal growth of the placenta and hybrid foetuses, a phenomenon known as hybrid placental dysplasia (HPD). The role of HPD in speciation is anticipated but still poorly understood. Here, we studied placental and foetal growth in F₁ crosses between four inbred mouse strains derived from two house mouse subspecies, Mus musculus musculus and Mus musculus domesticus. These subspecies are in the early stage of speciation and still hybridize in nature. In accordance with the maternal–foetal genomic conflict hypothesis, we found different parental influences on placental and foetal development, with placental weight most affected by the father's body weight and foetal weight by the mother's body weight. After removing the effects of parents’ body weight, we did not find any significant differences in foetal or placental weights between intra‐subspecific and inter‐subspecific F₁ crosses. Nevertheless, we found that the variability in placental weight in inter‐subspecific crosses is linked to the X chromosome, similarly as for HPD in interspecific mouse crosses. Our results suggest that maternal–foetal genomic conflict occurs in the house mouse system, but has not yet diverged sufficiently to cause abnormalities in placental and foetal growth in inter‐subspecific crosses. HPD is thus unlikely to contribute to speciation in the house mouse system. However, we cannot rule out that it might have contributed to other speciation events in the genus Mus, where differences in the levels of polyandry exist between the species.     

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.     

物种形成过程中的分化基因组岛及其形成机制

理解物种形成机制是生态和进化领域的重要任务。得益于测序技术的快速发展, 越来越多研究发现分化种群(亚种、物种)间的基因组常呈现异质性分化景观, 存在分化基因组岛, 这被认为是基因流存在下的歧化选择引起的, 支持基因流存在下的成种假说。然而, 基因渐渗、祖先多态性的差异分选、连锁选择等其他进化过程也可导致分化基因组岛的形成。现有实证研究在解析分化基因组岛的形成机制时, 往往忽略了上述其他进化过程的作用。为此, 本文在辨析分化基因组岛相关概念的基础上, 总结了利用种群基因组数据鉴定分化基因组岛的方法, 对比了不同进化过程形成分化基因组岛的特征, 指出在区分不同机制时联用基因渐渗程度、绝对分化指数(d_XY)、相对节点深度(RND)、重组率等多个指标的必要性, 归纳了物种形成过程中分化基因组岛形成机制解析的研究思路, 并对未来在生殖隔离机制上的深入探索以及实证研究的整合分析等方面进行了展望。     

A high frequency of allopolyploid speciation in the gymnospermous genus Ephedra and its possible association with some biological and ecological features

The origin and evolution of polyploids have been studied extensively in angiosperms and ferns but very rarely in gymnosperms. With the exception of three species of conifers, all natural polyploid species of gymnosperms belong to Ephedra, in which more than half of the species show polyploid cytotypes. Here, we investigated the origin and evolution of polyploids of Ephedra distributed in the Qinghai–Tibetan Plateau (QTP) and neighbouring areas. Flow cytometry (FCM) was used to measure the ploidy levels of the sampled species that are represented by multiple individuals from different populations, and then, two single‐copy nuclear genes (LFY and DDB2) and two chloroplast DNA fragments were used to unravel the possible origins and maternal donors of the polyploids. The results indicate that the studied polyploid species are allopolyploids, and suggest that allotetraploidy is a dominant mode of speciation in Ephedra. The high percentage of polyploids in the genus could be related to some of its biological attributes such as vegetative propagation, a relatively high rate of unreduced gamete formation, and a small genome size relative to most other gymnosperms. Significant ecological divergences between allotetraploids and their putative progenitors were detected by PCAs and anova and Tukey's tests, with the exception of E. saxatilis. The overlap of geographical distributions and ecological niches of some diploid species could have provided opportunities for interspecific hybridization and allopolyploid speciation.     

The origin of ecological divergence in Helianthus paradoxus (Asteraceae): selection on transgressive characters in a novel hybrid habitat

Diploid hybrid speciation in plants is often accompanied by rapid ecological divergence between incipient neospecies and their parental taxa. One plausible means by which novel adaptation in hybrid lineages may arise is transgressive segregation, that is, the generation of extreme phenotypes that exceed those of the parental lines. Early generation (BC2) hybrids between two wild, annual sunflowers, Helianthus annuus and Helianthus petiolaris, were used to study directional selection on transgressive characters associated with the origin of Helianthus paradoxus, a diploid hybrid species adapted to extremely saline marshes. The BC2 plants descended from a single F1 hybrid backcrossed toward H. petiolaris. The strength of selection on candidate adaptive traits in the interspecific BC2 was measured in natural H. paradoxus salt marsh habitat. Positive directional selection was detected for leaf succulence and Ca uptake, two traits that are known to be important in salt stress response in plants. Strong negative directional selection operated on uptake of Na and correlated elements. A significant decrease in trait correlations over time was observed in the BC2 population for Na and Ca content, suggesting an adaptive role for increased Ca uptake coupled with increased net exclusion of Na from leaves. Patterns of directional selection in BC2 hybrids were concordant with character expression in the natural hybrid species, H. paradoxus, transplanted into the wild. Moreover, the necessary variation for generating the H. paradoxus phenotype existed only in the BC2 population, but not in samples of the two parental species, H. annuus and H. petiolaris. These results are consistent with the hypothesis that transgressive segregation of elemental uptake and leaf succulence contributed to the origin of salt adaptation in the diploid hybrid species H. paradoxus.     

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.