Species delimitation, discovery and conservation in a tiger beetle species complex despite discordant genetic data.

In an age of species declines, delineating and discovering biodiversity is critical for both taxonomic accuracy and conservation. In recent years, there has been a movement away from using exclusively morphological characters to delineate and describe taxa and an increase in the use of molecular markers to describe diversity or through integrative taxonomy, which employs traditional morphological characters, as well as genetic or other data. Tiger beetles are charismatic, of conservation concern, and much work has been done on the morphological delineation of species and subspecies, but few of these taxa have been tested with genetic analyses. In this study, we tested morphologically based taxonomic hypotheses of polymorphic tiger beetles in the Eunota circumpicta (LaFerté-Sénectère, 1841) species complex using multilocus genomic and mtDNA analyses. We find multiple cryptic species within the previous taxonomic concept of Eunota circumpicta, some of which were historically recognized as subspecies. We found that the mtDNA and genomic datasets did not identify the same taxonomic units and that the mtDNA was most at odds with all other genetic and morphological patterns. Overall, we describe new cryptic diversity, which raises important conservation concerns, and provide a working example for testing species and subspecies validity despite discordant data.

A test of trade-offs in dispersal and reproduction within and between a sister species pair of specialist insect herbivores.

Understanding the drivers of trade-offs among traits is vital for comprehending the evolution and maintenance of trait variation. Theoretical frameworks propose that evolutionary mechanisms governing trade-offs frequently exhibit a scale-dependent nature. However, empirical tests of whether trade-offs exhibited across various biological scales (i.e. individuals, populations, species, genera, etc.) remains scarce. In this study, we explore trade-off between dispersal and reproductive effort among sympatric sister species of wasps in the genus Belonocnema (Hymenoptera: Cynipini: Cynipidae) that form galls on live oaks: B. fossoria, which specializes on Quercus geminata, and B. treatae, which specializes on Q. virginiana. Specifically, our results suggest that B. fossoria has evolved reduced flight capability and smaller wings, but a larger abdomen and greater total reproductive effort than B. treatae, which has larger wings and is a stronger flier, but has a smaller abdomen and reduced total reproductive effort. These traits and the relationships among them remain unchanged when B. fossoria and B. treatae are transplanted and reared onto the alternative host plant, suggesting that trait divergence is genetically based as opposed to being a plastic response to the different rearing environments. However, when looking within species, we found no evidence of intraspecific trade-offs between wing length and reproductive traits within either B. fossoria or B. treatae. Overall, our results indicate that observed trade-offs in life history traits between the two gall former species are likely a result of independent adaptations in response to different environments as opposed to the amplified expression of within species intrinsic tradeoffs.

The role of divergent host use and geography in the evolution of habitat isolation and sexual isolation among sister species of Belonocnema gall wasps.

Ecology and geography can play important roles in the evolution of reproductive isolation across the speciation continuum, but few studies address both at the later stages of speciation. This notable gap in knowledge arises from the fact that traditional ecological speciation studies have predominantly focused on the role of ecology in initiating the speciation process, while many studies exploring the effect of geography (e.g., reinforcement) concentrate on species pairs that lack divergent ecological characteristics. We simultaneously examine the strength of habitat isolation and sexual isolation among three closely related species of Belonocnema gall-forming wasps on two species of live oaks, Quercus virginiana and Q. geminata, that experience divergent selection from their host plants and variable rates of migration due to their geographic context. We find that the strength of both habitat isolation and sexual isolation is lowest among allopatric species pairs with the same host plant association, followed by allopatric species with different host plant associations, and highest between sympatric species with different host-plant associations. This pattern suggests that divergent selection due to different host use interacts with geography in the evolution of habitat isolation and sexual isolation during the later stages of speciation of Belonocnema wasps.

A genomic test of subspecies in the Eunota togata species group (Coleoptera: Cicindelidae): Morphology masks evolutionary relationships and taxonomy.

Most of the world's biodiversity is described primarily or exclusively using morphological traits that may not always reflect the true evolutionary units. Accurate taxonomy is critical for conservation efforts and re-evaluation of traditional taxonomy may often be warranted since species and subspecies are frequently the focus of conservation and faunistic studies. Here, we test comprehensive taxonomic hypotheses of morphologically defined subspecies in the tiger beetle, Eunota togata (LaFerté-Sénectère, 1841). The four recognized subspecies were delineated based mainly on the dorsal coloration and extent of white markings termed maculations. We combine inferences from mtDNA genealogies and genome-wide multilocus data to elucidate the evolutionary relationships within the group and assess the taxonomic implications. Three of the four subspecific taxa delineated by morphology were not supported by the genomic or mtDNA data. In fact, the species-level diversity in this group was underestimated, as E. togata was found to represent three well-supported distinct species in all genetic analyses. Emerging from these analyses, we also document an intriguing example of convergent evolution in lighter colored E. togata adapting to similar white saline backgrounds. Our collective work underscores the importance of using molecular methods to reevaluate morphological based taxonomy for species and subspecies delimitation and conservation.

Spatial sorting promotes rapid (mal)adaptation in the red-shouldered soapberry bug after hurricane-driven local extinctions.

Predicting future evolutionary change is a critical challenge in the Anthropocene as geographic range shifts and local extinction emerge as hallmarks of planetary change. Hence, spatial sorting-a driver of rapid evolution in which dispersal-associated traits accumulate along expanding range edges and within recolonized habitats-might be of growing importance in ecology and conservation. We report on the results of a natural experiment that monitored recolonization of host plants by the seed-feeding, red-shouldered soapberry bug, Jadera haematoloma, after local extinctions from catastrophic flooding in an extreme hurricane. We tested the contribution of spatial sorting to generate rapid and persistent evolution in dispersal traits, as well as in feeding traits unrelated to dispersal. Long-winged dispersal forms accumulated in recolonized habitats and due to genetic correlation, mouthparts also became longer and this shift persisted across generations. Those longer mouthparts were probably adaptive on one host plant species but maladaptive on two others based on matching the optimum depth of seeds within their host fruits. Moreover, spatial sorting eroded recently evolved adaptive divergence in mouthpart length among all host-associated biotypes, an outcome pointing to profound practical consequences of the extreme weather event for local adaptation, population resilience and evolutionary futures.

Suspended Materials Affect Particle Size Distribution and Removal of Environmental DNA in Flowing Waters.

Environmental DNA (eDNA) in aquatic systems is a complex mixture that includes dissolved DNA, intracellular DNA, and particle-adsorbed DNA. Information about the various components of eDNA and their relative proportions could be used to discern target organism abundance and location. However, a limited knowledge of eDNA adsorption dynamics and interactions with other materials hinders these applications. To address this gap, we used recirculating stream mesocosms to investigate the impact of suspended materials (fine particulate organic matter, plankton, clay, and titanium dioxide) on the eDNA concentration and particle size distribution (PSD) from two fish species in flowing water. Our findings revealed that eDNA rapidly adsorbs to other materials in the water column, affecting its concentration and PSD. Nonetheless, only particulate organic matter affected eDNA removal rate after 30 h. Moreover, we observed that the removal of larger eDNA components (≥10 µm) was more strongly influenced by physical processes, whereas the removal of smaller eDNA components was driven by biological degradation. This disparity in removal mechanisms between larger and smaller eDNA components could explain changes in eDNA composition over time and space, which have implications for modeling the spatial distribution and abundance of target species and optimizing eDNA detection in high turbidity systems.

Environmental DNA (eDNA) removal rates in streams differ by particle size under varying substrate and light conditions.

The use of environmental DNA (eDNA) as a sampling tool offers insights into the detection of invasive and/or rare aquatic species and enables biodiversity assessment without traditional sampling approaches, which are often labor-intensive. However, our understanding of the environmental factors that impact eDNA removal (i.e., how rapidly eDNA is removed from the water column by the combination of decay and physical removal) in flowing waters is limited. This limitation constrains predictions about the location and density of target organisms after positive detection. To address this question, we spiked Common Carp (Cyprinus carpio) eDNA into recirculating mesocosms (n = 24) under varying light (shaded versus open) and benthic substrate conditions (no substrate, bare substrate, and biofilm-colonized substrate). We then collected water samples from each mesocosm at four time points (40 min, 6 h, 18 h, and 48 h), and sequentially filtered the samples through 10, 1.0, and 0.2 µm filters to quantify removal rates for different eDNA particle sizes under varying light and substrate conditions. Combining all size classes, total eDNA removal rates were higher for mesocosms with biofilm-colonized substrate compared to those with no substrate or bare (i.e., no biofilm) substrate, which is consistent with previous findings linking biofilm colonization with increased eDNA removal and degradation. Additionally, when biofilm was present, light availability increased eDNA removal; eDNA levels fell below detection after 6-18 h for open mesocosms versus 18-48 h for shaded mesocosms. Among size classes, larger particles (>10 µm) were removed faster than small particles (1.0-0.2 µm). These results suggest that changes in the distribution of eDNA size classes over time (e.g., with downstream transport) and with differing environmental conditions could be used to predict the location of target organisms in flowing waters, which will advance the use of eDNA as a tool for species monitoring and management.

Environment and shipping drive environmental DNA beta-diversity among commercial ports.

The spread of nonindigenous species by shipping is a large and growing global problem that harms coastal ecosystems and economies and may blur coastal biogeographical patterns. This study coupled eukaryotic environmental DNA (eDNA) metabarcoding with dissimilarity regression to test the hypothesis that ship-borne species spread homogenizes port communities. We first collected and metabarcoded water samples from ports in Europe, Asia, Australia and the Americas. We then calculated community dissimilarities between port pairs and tested for effects of environmental dissimilarity, biogeographical region and four alternative measures of ship-borne species transport risk. We predicted that higher shipping between ports would decrease community dissimilarity, that the effect of shipping would be small compared to that of environment dissimilarity and shared biogeography, and that more complex shipping risk metrics (which account for ballast water and stepping-stone spread) would perform better. Consistent with our hypotheses, community dissimilarities increased significantly with environmental dissimilarity and, to a lesser extent, decreased with ship-borne species transport risks, particularly if the ports had similar environments and stepping-stone risks were considered. Unexpectedly, we found no clear effect of shared biogeography, and that risk metrics incorporating estimates of ballast discharge did not offer more explanatory power than simpler traffic-based risks. Overall, we found that shipping homogenizes eukaryotic communities between ports in predictable ways, which could inform improvements in invasive species policy and management. We demonstrated the usefulness of eDNA metabarcoding and dissimilarity regression for disentangling the drivers of large-scale biodiversity patterns. We conclude by outlining logistical considerations and recommendations for future studies using this approach.

Particle size influences decay rates of environmental DNA in aquatic systems.

Environmental DNA (eDNA) analysis is a powerful tool for remote detection of target organisms. However, obtaining quantitative and longitudinal information from eDNA data is challenging, requiring a deep understanding of eDNA ecology. Notably, if the various size components of eDNA decay at different rates, and we can separate them within a sample, their changing proportions could be used to obtain longitudinal dynamics information on targets. To test this possibility, we conducted an aquatic mesocosm experiment in which we separated fish-derived eDNA components using sequential filtration to evaluate the decay rate and changing proportion of various eDNA particle sizes over time. We then fit four alternative mathematical decay models to the data, building towards a predictive framework to interpret eDNA data from various particle sizes. We found that medium-sized particles (1-10 µm) decayed more slowly than other size classes (i.e., <1 and > 10 µm), and thus made up an increasing proportion of eDNA particles over time. We also observed distinct eDNA particle size distribution (PSD) between our Common carp and Rainbow trout samples, suggesting that target-specific assays are required to determine starting eDNA PSDs. Additionally, we found evidence that different sizes of eDNA particles do not decay independently, with particle size conversion replenishing smaller particles over time. Nonetheless, a parsimonious mathematical model where particle sizes decay independently best explained the data. Given these results, we suggest a framework to discern target distance and abundance with eDNA data by applying sequential filtration, which theoretically has both metabarcoding and single-target applications.

Speciation in Nearctic oak gall wasps is frequently correlated with changes in host plant, host organ, or both.

Quantifying the frequency of shifts to new host plants within diverse clades of specialist herbivorous insects is critically important to understand whether and how host shifts contribute to the origin of species. Oak gall wasps (Hymenoptera: Cynipidae: Cynipini) comprise a tribe of ∼1000 species of phytophagous insects that induce gall formation on various organs of trees in the family Fagacae-primarily the oaks (genus Quercus; ∼435 sp.). The association of oak gall wasps with oaks is ancient (∼50 my), and most oak species are galled by one or more gall wasp species. Despite the diversity of both gall wasp species and their plant associations, previous phylogenetic work has not identified the strong signal of host plant shifting among oak gall wasps that has been found in other phytophagous insect systems. However, most emphasis has been on the Western Palearctic and not the Nearctic where both oaks and oak gall wasps are considerably more species rich. We collected 86 species of Nearctic oak gall wasps from most of the major clades of Nearctic oaks and sequenced >1000 Ultraconserved Elements (UCEs) and flanking sequences to infer wasp phylogenies. We assessed the relationships of Nearctic gall wasps to one another and, by leveraging previously published UCE data, to the Palearctic fauna. We then used phylogenies to infer historical patterns of shifts among host tree species and tree organs. Our results indicate that oak gall wasps have moved between the Palearctic and Nearctic at least four times, that some Palearctic wasp clades have their proximate origin in the Nearctic, and that gall wasps have shifted within and between oak tree sections, subsections, and organs considerably more often than previous data have suggested. Given that host shifts have been demonstrated to drive reproductive isolation between host-associated populations in other phytophagous insects, our analyses of Nearctic gall wasps suggest that host shifts are key drivers of speciation in this clade, especially in hotspots of oak diversity. Although formal assessment of this hypothesis requires further study, two putatively oligophagous gall wasp species in our dataset show signals of host-associated genetic differentiation unconfounded by geographic distance, suggestive of barriers to gene flow associated with the use of alternative host plants.

Delimiting the cryptic diversity and host preferences of Sycophila parasitoid wasps associated with oak galls using phylogenomic data.

Cryptic species diversity is a major challenge regarding the species-rich community of parasitoids attacking oak gall wasps due to a high degree of sexual dimorphism, morphological plasticity, small size and poorly known biology. As such, we know very little about the number of species present, nor the evolutionary forces responsible for generating this diversity. One hypothesis is that trait diversity in the gall wasps, including the morphology of the galls they induce, has evolved in response to selection imposed by the parasitoid community, with reciprocal selection driving diversification of the parasitoids. Using a rare, continental-scale data set of Sycophila parasitoid wasps reared from 44 species of cynipid galls from 18 species of oak across the USA, we combined mitochondrial DNA barcodes, ultraconserved elements (UCEs), morphological and natural history data to delimit putative species. Using these results, we generate the first large-scale assessment of ecological specialization and host association in this species-rich group, with implications for evolutionary ecology and biocontrol. We find most Sycophila target specific subsets of available cynipid host galls with similar morphologies, and generally attack larger galls. Our results suggest that parasitoid wasps such as Sycophila have adaptations allowing them to exploit particular host trait combinations, while hosts with contrasting traits are resistant to attack. These findings support the tritrophic niche concept for the structuring of plant-herbivore-parasitoid communities.

Novel Field-Based Protein Detection Method Using Light Transmission Spectroscopy and Antibody Functionalized Gold Nanoparticles.

Protein detection is a universal tool critical to many applications in medicine, agriculture, and biotechnology. We developed a novel protein detection method combining light transmission spectroscopy and particle-size analysis of gold nanospheres monovalently functionalized with polyclonal antibodies and applied it to an emerging challenge for such technologies─the monitoring of environmental proteins (eProteins) present in natural aquatic systems. These are an underreported source of pollution and include the pseudopersistent Cry toxins that enter aquatic ecosystems from surrounding genetically engineered crops. The assay is capable of detecting proteins in complex matrices, such as water samples collected in the field, making it a competitive assay for eProtein detection. It is sensitive, reaching 1.25 ng mL-1, and we demonstrate its application to the detection of Cry1Ab from subsurface tile-drain and streamwater samples from agricultural waterways. The assay can also be quickly adapted for other protein detection applications in the future.

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.

Ultraconserved elements-based systematics reveals evolutionary patterns of host-plant family shifts and phytophagy within the predominantly parasitoid braconid wasp subfamily Doryctinae.

Phytophagy has promoted species diversification in many insect groups, including Hymenoptera, one of the most diverse animal orders on Earth. In the predominantly parasitoid family Braconidae, an association with insect-induced, plant galls in angiosperms have been reported in three subfamilies, but in particular in the Doryctinae, where it has been recorded to occur in species of ten genera. Allorhogas Gahan is the most species-rich of these genera, with its species having different phytophagous strategies. Here we conducted a comprehensive phylogenomic study for the doryctine gall-associated genera, with an emphasis on Allorhogas, using ultraconserved elements (UCEs). Based on this estimate of phylogeny we: (1) evaluated their taxonomic composition, (2) estimated the timing of origin of the gall-associated clade and divergence of its main subclades, and (3) performed ancestral state reconstruction analyses for life history traits related to their host-plant association. Our phylogenetic hypothesis confirmed Allorhogas as polyphyletic, with most of its members being nested in a main clade composed of various subclades, each comprising species with a particular host-plant family and herbivorous feeding habit. The origin of gall-association was estimated to have occurred during the late Oligocene to early Miocene, with a subsequent diversification of subclades during the middle to late Miocene and Pliocene. Overlap in divergence timing appears to occur between some taxa and their host-associated plant lineages. Evolution of the feeding strategies in the group shows "inquilinism-feeding" as the likely ancestral state, with gall-formation in different plant organs and seed predation having independently evolved on multiple occasions.

Tissue-specific gene expression shows a cynipid wasp repurposes oak host gene networks to create a complex and novel parasite-specific organ.

Every organism on Earth depends on interactions with other organisms to survive. In each of these interactions, an organism must utilize the limited toolbox of genes and proteins it possesses to successfully manipulate or cooperate with another species, but it can also co-opt the genome machinery of its partner to expand its available tools. Insect-induced plant galls are an extreme example of this, wherein an insect hijacks the plant's genome to direct the initiation and development of galls consisting of plant tissue. However, previous transcriptomic studies have not evaluated individual tissues within a gall to determine the full extent to which a galling insect manipulates its host plant. Here we demonstrate that the cynipid wasp Dryocosmus quercuspalustris creates a complex parasite-specific organ from red oak tissue via massive changes in host gene expression. Our results show that the gall wasp is not merely modifying oak leaf tissue but creating extensive changes in gene expression between galled and ungalled tissue (differential expression in 28% of genes) and distinct gall tissue types (20% of genes). The outer gall tissue shows increases in various plant defence systems, which is consistent with its predicted functional role of protecting the wasp larva. The inner larval capsule shows suppression of large parts of the plant innate immune system and evidence for the wasp utilizing the plant's RNA interference mechanisms, which may be a potential mechanism for the wasp's control on gall growth.

The Arthropod Associates of 155 North American Cynipid Oak Galls.

The identities of most arthropod associates of cynipid-induced oak galls in the western Palearctic are generally known. However, a comprehensive accounting of associates has been performed for only a small number of the galls induced by the estimated 700 species of cynipid gall wasps in the Nearctic. This gap in knowledge stymies many potential studies of diversity, coevolution, and community ecology, for which oak gall systems are otherwise ideal models. We report rearing records of insects and other arthropods from more than 527,306 individual galls representing 201 different oak gall types collected from 32 oak tree species in North America. Of the 201 gall types collected, 155 produced one or more arthropods. A total of 151,075 arthropods were found in association with these 155 gall types, and of these 61,044 (40.4%) were gall wasps while 90,031 (59.6%) were other arthropods. We identified all arthropods to superfamily, family, or, where possible, to genus. We provide raw numbers and summaries of collections, alongside notes on natural history, ecology, and previously published associations for each taxon. For eight common gall-associated genera (Synergus, Ceroptres, Euceroptres, Ormyrus, Torymus, Eurytoma, Sycophila, and Euderus), we also connect rearing records to gall wasp phylogeny, geography, and ecology -including host tree and gall location (host organ), and their co-occurrence with other insect genera. Though the diversity of gall wasps and the large size of these communities is such that many Nearctic oak gall-associated insects still remain undescribed, this large collection and identification effort should facilitate the testing of new and varied ecological and evolutionary hypotheses in Nearctic oak galls.

Context-Dependent Reproductive Isolation: Host Plant Variability Drives Fitness of Hybrid Herbivores.

AbstractThe role of divergent selection between alternative environments in promoting reproductive isolation (RI) between lineages is well recognized. However, most studies view each divergent environment as homogenous, thereby overlooking the potential role within-environment variation plays in RI between differentiating lineages. Here, we test the importance of microenvironmental variation in RI by using individual trees of two host plants, each harboring locally adapted populations of the cynipid wasp Belonocnema treatae. We compared the fitness surrogate (survival) of offspring from hybrid crosses with resident crosses across individual trees on each of two primary host plants, Quercus virginiana and Q. geminata. We found evidence of weak hybrid inviability between host-associated lineages of B. treatae despite strong genomic differentiation. However, averaging across environments masked great variation in hybrid fitness on individual trees, where hybrids performed worse than, equal to, or better than residents. Thus, considering the environmental context of hybridization is critical to improving the predictability of divergence under variable selection.

Fate of Environmental Proteins (eProteins) from Genetically Engineered Crops in Streams is Controlled by Water pH and Ecosystem Metabolism.

Environmental proteins (eProteins), such as Cry proteins associated with genetically engineered (GE) organisms, are present in ecosystems worldwide, but only rarely reach concentrations with detectable ecosystem-level impacts. Despite their ubiquity, the degradation and fate of Cry and other eProteins are mostly unknown. Here, we report the results of an experiment where we added Cry proteins leached from GE Bt maize to a suite of 19 recirculating experimental streams. We found that Cry exhibited a biphasic degradation with an initial phase of rapid and variable degradation within 1 h, followed by a slow and steady phase of degradation with traces of protein persisting after 48 h. The initial degradation was correlated with heterotrophic respiration and water column dissolved oxygen, confirming a previously documented association with stream metabolism. However, protein degradation persisted even with no biofilm and was faster at a more acidic pH, suggesting that water chemistry is also a critical factor in both degradation and subsequent detection. We suggest that Cry, as well as other eProteins, will have a rapid degradation caused by denaturation of proteins and pH changes, which confirms that the detection of Cry proteins in natural streams must be the result of steady and consistent leaching into the environment.

Asymmetric, but opposing reductions in immigrant viability and fecundity promote reproductive isolation among host-associated populations of an insect herbivore.

Immigrant inviability can contribute to reproductive isolation (RI) during ecological speciation by reducing the survival of immigrants in non-native environments. However, studies that assess the fitness consequence of immigrants moving from native to non-native environments typically fail to explore the potential role of concomitant reductions in immigrant fecundity despite recent evidence suggesting its prominent role during local adaptation. Here, we evaluate the directionality and magnitude of both immigrant viability and fecundity to RI in a host-specific gall-forming wasp, Belonocnema treatae. Using reciprocal transplant experiments replicated across sites, we measure immigrant viability and fecundity by comparing differences in the incidence of gall formation (viability) and predicted the number of eggs per female (fecundity) between residents and immigrants in each of two host-plant environments. Reduced immigrant viability was found in one host environment while reduced immigrant fecundity was found in the other. Such habitat-dependent barriers resulted in asymmetric RI between populations. By surveying recent literature on local adaptation, we find that asymmetry in immigrant viability and fecundity are widespread across disparate taxa, which highlights the need to combine estimates of both common and overlooked barriers in cases of potential bi-directional gene flow to create a more comprehensive view of the evolution of RI.

Evidence for spatial clines and mixed geographic modes of speciation for North American cherry-infesting Rhagoletis (Diptera: Tephritidae) flies.

An important criterion for understanding speciation is the geographic context of population divergence. Three major modes of allopatric, parapatric, and sympatric speciation define the extent of spatial overlap and gene flow between diverging populations. However, mixed modes of speciation are also possible, whereby populations experience periods of allopatry, parapatry, and/or sympatry at different times as they diverge. Here, we report clinal patterns of variation for 21 nuclear-encoded microsatellites and a wing spot phenotype for cherry-infesting Rhagoletis (Diptera: Tephritidae) across North America consistent with these flies having initially diverged in parapatry followed by a period of allopatric differentiation in the early Holocene. However, mitochondrial DNA (mtDNA) displays a different pattern; cherry flies at the ends of the clines in the eastern USA and Pacific Northwest share identical haplotypes, while centrally located populations in the southwestern USA and Mexico possess a different haplotype. We hypothesize that the mitochondrial difference could be due to lineage sorting but more likely reflects a selective sweep of a favorable mtDNA variant or the spread of an endosymbiont. The estimated divergence time for mtDNA suggests possible past allopatry, secondary contact, and subsequent isolation between USA and Mexican fly populations initiated before the Wisconsin glaciation. Thus, the current genetics of cherry flies may involve different mixed modes of divergence occurring in different portions of the fly's range. We discuss the need for additional DNA sequencing and quantification of prezygotic and postzygotic reproductive isolation to verify the multiple mixed-mode hypothesis for cherry flies and draw parallels from other systems to assess the generality that speciation may commonly involve complex biogeographies of varying combinations of allopatric, parapatric, and sympatric divergence.