Salicaceae detoxification abilities in Florida tiger swallowtail butterflies (Papilio glaucus maynardi Gauthier): Novel ability or Pleistocene holdover?

Abstract Florida populations of the eastern tiger swallowtail butterfly, Papilio glaucus L., have unique morphological features and ecological adaptations that have contributed to their subspecies status (P. g. maynardi Gauthier). We describe geographically unique abilities for detoxification of Carolina willow, Salix caroliniana Michx. (Salicaceae), for several Florida populations of P. g. maynardi. Of all the approximately 570 worldwide species of the Papilionidae, such Salicaceae detoxification abilities exist only in the allopatric North American western and northernmost species (P. rutulus Lucas, P. eurymedon Lucas and P. canadensis Rothschild & Jordan). Females of P. glaucus collected from populations in southeastern USA were examined for oviposition preference in 5‐choice assays, and displayed a low preference for Salicaceae (<5%), but larvae from Florida populations exhibited a high survival (>60%) on these plants. Detoxification abilities have previously shown to be autosomally inherited, and can be transferred via natural or hand‐paired interspecific hybrid introgression. However, these Florida populations are at least 700–1 500 km from the nearest hybrids or the hybrid species, P. appalachiensis Pavulaan & Wright, which possess these detoxification abilities. In any case, the Z ( = X)‐linked oviposition preferences for Salicaceae are lacking in these Florida populations, illustrating genetic independence of oviposition preference determination and larval survival/performance abilities. The origins of detoxification abilities are unlikely to be due to recent climate‐driven introgression, and may represent ancestral trait carry‐overs from interglacial refugium populations of the Pleistocene epoch.     

Assessing ecological and physiological costs of melanism in North American Papilio glaucus females: two decades of dark morph frequency declines

Polymorphisms for melanic form of insects may provide various selective advantages. However, melanic alleles may have significant/subtle pleiotrophic “costs.” Several potential pleiotrophic effects of the W (=Y)‐linked melanism gene in Papilio glaucus L. (Lepidoptera) showed no costs for melanic versus yellow in adult size, oviposition preferences, fecundity, egg viability, larval survival/growth rates, cold stress tolerance, or postdiapause emergence times. Sexual selection (males choosing yellow rather than mimetic dark females) had been suggested to provide a balanced polymorphism in P. glaucus, but spermatophore counts in wild females and direct field tethering studies of size‐matched pairs of virgin females (dark and yellow), show that male preferences are random or frequency‐dependent from Florida to Michigan, providing no yellow counter‐advantages. Recent frequency declines of dark (melanic/mimetic) females in P. glaucus populations are shown in several major populations from Florida (27.3°N latitude) to Ohio (38.5° N). Summer temperatures have increased significantly at all these locations during this time (1999–2018), but whether dark morphs may be more vulnerable (in any stage) to such climate warming remains to be determined. Additional potential reasons for the frequency declines in mimetic females are discussed: (i) genetic introgression of Z‐linked melanism suppressor genes from P. canadensis (R & J) and the hybrid species, P. appalachiensis (Pavulaan & Wright), (ii) differential developmental incompatibilities, or Haldane effects, known to occur in hybrids, (iii) selection against intermediately melanic (“dusty”) females (with the W‐linked melanic gene, b+) which higher temperatures can cause.     

POPULATION STRUCTURE AND HOST USE IN HYBRIDIZING SUBSPECIES OF PAPILIO GLAUCUS (LEPIDOPTERA: PAPILIONIDAE)

Two parapatric subspecies of the eastern tiger swallowtail butterfly, Papilio glaucus glaucus and P. g. canadensis, differ greatly in larval survival and growth on host plants in the Magnoliaceae, Salicaceae, and Betulaceae. In the first part of this study, butterflies were collected from 17 sites along a transect from Georgia to Quebec and used for allozyme electrophoresis. Two X- (or Z-)linked enzyme loci (Ldh and Pgd) showed allele frequency differences that delineated a hybrid zone between the subspecies in northern Pennsylvania and south-central New York. No significant linkage disequilibrium could be detected among allozyme loci within the hybrid zone samples. Genetic differentiation at Ldh and Pgd coincides with subspecies differences in diapause control and female mimicry, which are also sex-linked. Larval offspring of butterflies from 13 of the sites were then tested in the laboratory for survival and growth on Liriodendron tulipifera (Magnoliaceae), Populus deltoïdes (Salicaceae) and Betula lutea (Betulaceae). Steep clines in survival rates matched the position of the hybrid zone. Hybrid zone larvae showed intermediate survival rates and significant heterogeneity among families on all three plants, indicating presence of substantial genetic variation. The results suggest that differential host use by P. g. glaucus and P. g. canadensis is maintained primarily by independent clines in selection intensity for ability to use allopatric sets of host plants, coupled with restricted gene flow across the hybrid zone.     

Bayesian estimation of genomic clines

We developed a Bayesian genomic cline model to study the genetic architecture of adaptive divergence and reproductive isolation between hybridizing lineages. This model quantifies locus‐specific patterns of introgression with two cline parameters that describe the probability of locus‐specific ancestry as a function of genome‐wide admixture. ‘Outlier’ loci with extreme patterns of introgression relative to most of the genome can be identified. These loci are potentially associated with adaptive divergence or reproductive isolation. We simulated genetic data for admixed populations that included neutral introgression, as well as loci that were subject to directional, epistatic or underdominant selection, and analysed these data using the Bayesian genomic cline model. Under many demographic conditions, underdominance or directional selection had detectable and predictable effects on cline parameters, and ‘outlier’ loci were greatly enriched for genetic regions affected by selection. We also analysed previously published genetic data from two transects through a hybrid zone between Mus domesticus and M. musculus. We found considerable variation in rates of introgression across the genome and particularly low rates of introgression for two X‐linked markers. There were similarities and differences in patterns of introgression between the two transects, which likely reflects a combination of stochastic variability because of genetic drift and geographic variation in the genetic architecture of reproductive isolation. By providing a robust framework to quantify and compare patterns of introgression among genetic regions and populations, the Bayesian genomic cline model will advance our understanding of the genetics of reproductive isolation and the speciation process.     

Hybridization leads to host-use divergence in a polyphagous butterfly sibling species pair

Climate warming has lead to increased genetic introgression across a narrow hybrid zone separating the eastern and Canadian tiger swallowtails (Papilio glaucus and Papilio canadensis). This situation has led to the formation of an allochronically separated hybrid population with a delayed emerging phenotype or “late flight”. Here, we assess how the recombination of the parental genomes that lead to this phenotype may have facilitated another major ecological shift, host-use divergence. We first contrast the ovipositional profiles of the late flight population to that of the parental species P. glaucus and P. canadensis. Subsequently we contrast the larval survival and growth of the late flight, a P. canadensis and a P. glaucus population, and a population from the northern edge of the hybrid zone on five hosts. Our results indicate that the ovipositional preference of this hybrid swarm is identical to that of the introgressing parental species, P. glaucus. Due to the absence of the preferred hosts of P. glaucus (Liriodendron tulipifera L. and Ptelea trifoliata L.) where the late flight occurs, this ovipositional pattern implies a functional specialization onto a secondary host of both parental species, Fraxinus americana L. In contrast, the larval host-use abilities represent a mixture of P. glaucus and P. canadensis, indicating divergence in larval host-use abilities has not taken place. However, high genetic variability (genetic coefficient of variation) is present for growth on F. americana in the late flight hybrid swarm and tradeoffs for larval performance on the preferred hosts of the parental species are evident; indicating a strong potential for future specialization in larval host-use abilities. This current scenario represents an instance where a shift in a major ecological trait, host use, is likely occurring as a byproduct of a shift in an unrelated trait (delayed emergence) leading to partial reproductive isolation.     

GENETICS OF MIMICRY IN THE TIGER SWALLOWTAIL BUTTERFLIES,PAPILIO GLAUCUS AND P. CANADENSIS (LEPIDOPTERA: PAPILIONIDAE)

The tiger swallowtail butterfly, Papilio glaucus, exhibits a female-limited polymorphism for Batesian mimicry; the Canadian tiger swallowtail, Papilio canadensis, lacks the mimetic (dark) form entirely. The species hybridize to a limited extent where their ranges overlap. Field collections and censuses indicate that mimetic females occur throughout the range of P. glaucus but at lowest frequencies in populations at the latitudinal edges of its geographic range such as the southernmost part of Florida and along the entire northern edge of its distribution from Massachusetts to Minnesota. Frequencies of mimetic females have remained relatively stable over time. Inheritance of the mimetic form is controlled primarily by two interacting sex-linked loci. The typical matrilineal pattern of inheritance in P. glaucus can be explained by polymorphism at a Y-linked locus, b. Analysis of P. glaucus × P. canadensis crosses has also revealed an X-linked locus, s, which controls the expression of the mimetic phenotype. The P. canadensis allele, s^can, suppresses the mimetic phenotype in hybrid and backcross females. Results from more than 12 yr of rearing tiger swallowtails, including interspecies hybrids, indicate that the absence of mimetic P. canadensis females is due to both a high frequency of the “suppressing” allele s^can and low frequency of the black-pigment-determining b + allele. The frequency of s^can (or other suppressing alleles of s) in P. glaucus populations outside the hybrid zone is low. Some males heterozygous at the s locus and some suppressed mimetic females occur within the hybrid zone. A simple genetic model predicts the frequency of daughters that differ in phenotype from their mothers.     

Steep clines within a highly permeable genome across a hybrid zone between two subspecies of the European rabbit

Maintenance of genetic distinction in the face of gene flow is an important aspect of the speciation process. Here, we provide a detailed spatial and genetic characterization of a hybrid zone between two subspecies of the European rabbit. We examined patterns of allele frequency change for 22 markers located on the autosomes, X‐chromosome, Y‐chromosome and mtDNA in 1078 individuals sampled across the hybrid zone. While some loci revealed extremely wide clines (w ≥ 300 km) relative to an estimated dispersal of 1.95–4.22 km/generation, others showed abrupt transitions (w ≈ 10 km), indicating localized genomic regions of strong selection against introgression. The subset of loci showing steep clines had largely coincident centers and stepped changes in allele frequency that did not co‐localize with any physical barrier or ecotone, suggesting that the rabbit hybrid zone is a tension zone. The steepest clines were for X‐ and Y‐chromosome markers. Our results are consistent with previous inference based on DNA sequence variation of individuals sampled in allopatry in suggesting that a large proportion of each genome has escaped the overall barrier to gene flow in the middle of the hybrid zone. These results imply an old history of hybridization and high effective gene flow and anticipate that isolation factors should often localize to small genomic regions.     

The role of hybridization in the origin and spread of asexuality in Daphnia

The molecular mechanisms leading to asexuality remain little understood despite their substantial bearing on why sexual reproduction is dominant in nature. Here, we examine the role of hybridization in the origin and spread of obligate asexuality in Daphnia pulex, arguably the best‐documented case of contagious asexuality. Obligately parthenogenetic (OP) clones of D. pulex have traditionally been separated into ‘hybrid’ (Ldh SF) and ‘nonhybrid’ (Ldh SS) forms because the lactase dehydrogenase (Ldh) locus distinguishes the cyclically parthenogenetic (CP) lake dwelling Daphnia pulicaria (Ldh FF) from its ephemeral pond dwelling sister species D. pulex (Ldh SS). The results of our population genetic analyses based on microsatellite loci suggest that both Ldh SS and SF OP individuals can originate from the crossing of CP female F₁ (D. pulex × D. pulicaria) and backcross with males from OP lineages carrying genes that suppress meiosis specifically in female offspring. In previous studies, a suite of diagnostic markers was found to be associated with OP in Ldh SS D. pulex lineages. Our association mapping supports a similar genetic mechanism for the spread of obligate parthenogenesis in Ldh SF OP individuals. Interestingly, our study shows that CP D. pulicaria carry many of the diagnostic microsatellite alleles associated with obligate parthenogenesis. We argue that the assemblage of mutations that suppress meiosis and underlie obligate parthenogenesis in D. pulex originated due to a unique historical hybridization and introgression event between D. pulex and D. pulicaria.     

The role of latitudinal,genetic and temperature variation in the induction of diapause of Papilio glaucus (Lepidoptera: Papilionidae)

A key adaptation in insects for dealing with variable environmental conditions is the ability to diapause. The tiger swallowtail butterflies, Papilio glaucus and P. canadensis are ideal species to explore the genetic causes and population genetic consequences of diapause because divergence in this trait is believed to be a salient factor in maintaining a hybrid zone between these species. Yet little is known about the factors that influence diapause induction in this system. Here we explored how spatial (latitudinal), environmental (temperature) and genetic (hybridization) factors affect diapause induction in this system. Specifically, a series of growth chamber experiments using wild caught individuals from across the eastern United States were performed to: (1) evaluate how critical photoperiod varies with latitude, (2) isolate the stage in which induction occurs, (3) test whether changes in temperature affected rates of diapause induction, and (4) explore how the incidence of diapause is affected in hybrid offspring. We find that induction occurs in the larval stage, is not sensitive to a relatively broad range of temperatures, appears to have a complex genetic basis (i.e., is not simply a dominant trait following a Mendelian inheritance pattern) and that the critical photoperiod increases by 0.4 h with each increasing degree in latitude. This work deepens our understanding of how spatial, environmental and genetic variation influences a key seasonal adaptation (diapause induction) in a well‐developed ecological model system and will make possible future studies that explore how climatic variation affects the population dynamics and genetics of this system.     

Differential effects of short term winter thermal stress on diapausing tiger swallowtail butterflies (Papilio spp.)

Abstract It is generally thought that insects inhabiting lower latitudes are more severely impacted by changes in their thermal environment than are high latitude species. This is attributed to the wider range of temperatures to which high‐latitude species are exposed. By contrast, low‐latitude species have typically evolved in more thermally stable environments with a narrower range of temperature variation. However, deviation from this pattern can occur and here we report that under variable winter conditions a higher latitude species may be more sensitive to thermal variation than its lower latitude sister species. Using split broods, we examined the survival and adult emergence success of diapausing pupae of Papilio canadensis and P. glaucus, as well as a unique, recombinant hybrid population (“late‐flight”) to short periods of mid‐winter cold and heat stress. Our results indicate that the higher latitude, univoltine populations (P. canadensis and late‐flights) exhibit lower pupal survival than the lower latitude, facultative diapauser (P. glaucus) for all mid‐winter thermal stress treatments, both high and low. Size differences alone do not appear to account for the observed differences in survival or metabolic costs in these three phenotypes, as late‐flight individuals are similar in size to P. glaucus. We attribute the observed differences in survival and weight loss to potential metabolic differences and variation in the intensity of diapause, in addition to divergent adaptation to winter precipitation levels (e.g. snow cover) and the influences this may have on microhabitat temperature moderation.     

Spatiotemporally consistent genomic signatures of reproductive isolation in a moving hybrid zone

Studies of hybrid zone dynamics often investigate a single sampling period and draw conclusions from that temporal snapshot. Stochasticity can, however, result in loci with spurious outlier patterns, which is exacerbated by limited temporal or geographic sampling. Comparing admixed populations from different geographic regions is one way to detect repeatedly divergent genomic regions potentially involved in reproductive isolation. Temporal comparisons also allow us to control partially for the role of stochasticity, but the power of temporal sampling has not yet been adequately explored. In North America, black‐capped (Poecile atricapillus) and Carolina (P. carolinensis) chickadees hybridize in a contact zone extending from New Jersey to Kansas. The hybrid zone is likely maintained by strong intrinsic selection against hybrids, and it is moving north. We used a reduced representation genomic approach and temporally spaced sampling—two samples of ～80 individuals separated by a decade—to determine the pattern and consistency of selection and genomic introgression in the chickadee hybrid zone. We report consistently low introgression for highly divergent loci between P. atricapillus and P. carolinensis in this moving hybrid zone. This is strong evidence that these loci may be linked to genomic regions involved in reproductive isolation between chickadees.     

Reinforcement selection acting on the European house mouse hybrid zone

Behavioural isolation may lead to complete speciation when partial postzygotic isolation acts in the presence of divergent‐specific mate‐recognition systems. These conditions exist where Mus musculus musculus and M. m. domesticus come into contact and hybridize. We studied two mate‐recognition signal systems, based on urinary and salivary proteins, across a Central European portion of the mouse hybrid zone. Introgression of the genomic regions responsible for these signals: the major urinary proteins (MUPs) and androgen binding proteins (ABPs), respectively, was compared to introgression at loci assumed to be nearly neutral and those under selection against hybridization. The preference of individuals taken from across the zone regarding these signals was measured in Y mazes, and we develop a model for the analysis of the transition of such traits under reinforcement selection. The strongest assortative preferences were found in males for urine and females for ABP. Clinal analyses confirm nearly neutral introgression of an Abp locus and two loci closely linked to the Abp gene cluster, whereas two markers flanking the Mup gene region reveal unexpected introgression. Geographic change in the preference traits matches our reinforcement selection model significantly better than standard cline models. Our study confirms that behavioural barriers are important components of reproductive isolation between the house mouse subspecies.     

Diapause strategy in the parasitoid Psyllaephagus pistaciae

The processes involved in the induction and termination of diapause in the parasitoid Psyllaephagus pistaciae Ferrière (Hymenoptera: Encyrtidae) were investigated under controlled laboratory conditions. Results indicated that the parasitoid P. pistaciae is able to overwinter successfully as a pupa within the mummified psyllid host Agonoscena pistaciae Burckhardt and Lauterer (Homoptera: Psylloidea), and the parasitoid clearly responded to a short photoperiod by entering diapause. The incidence of diapause increased significantly with decreasing temperature, indicating the existence of an interaction between low temperature and short‐day photoperiod that resulted in a 100% diapause. However, low temperature alone partly stimulated P. pistaciae to enter diapause. The incidence of diapause was likely to be independent from psyllid nymphal instar and the condition of its mother. This investigation showed that the termination of diapause in the overwintering parasitoid pupae and the duration of diapause completion or adult emergence is photoperiod independent, but dependent on temperature. However, chilling treatment was found to be unnecessary for diapause termination in P. pistaciae, although it decreased the duration of diapause maturation. The positive responses of parasitoids to environmental factors appears to act as a basic strategy to induce diapause in the parasitoid P. pistaciae against extreme climatic conditions.     

Genetics of diapause in the comma butterfly Polygonia c‐album

The processes of local adaptation and ecological speciation can be better understood by studying the genetic background of life‐history decisions. The sex chromosomes host genes for many population differences in the Lepidoptera and therefore the inheritance of diapause determination in the butterfly Polygonia c‐album may be hypothesized to be sex‐linked. In the present study, Polygonia c‐album (L.) from Spain and Sweden and hybrid offspring are raised under an LD 17 : 7 h photocycle that induces most pure Swedish individuals to develop into the diapausing dark morph and most pure Spanish individuals into the light and directly‐developing morph. If inheritance of the daylength threshold for diapause is X‐linked, as is known to be the case for host‐plant preferences, females should follow the developmental path of their male parents' populations. However, female hybrids instead have a diapause propensity intermediate to that of their parental stocks and, consequentially, diapause determination is not X‐linked. However, male hybrids eclose as the diapausing morph to a higher extent than females and, moreover, this pattern is more pronounced in the Spanish female × Swedish male cross than in the reciprocal cross. Hence, it is concluded that the genetic determination of the critical daylength for diapause is mainly autosomal but with some influence of sex‐linked genes and/or parental effects, possibly as an effect of the importance of protandry for males. Such sex effects could provide a starting point for the evolution of population differences inherited on the sex chromosomes.     

Physiological adaptation along environmental gradients and replicated hybrid zone structure in swordtails (Teleostei: Xiphophorus)

Local adaptation is often invoked to explain hybrid zone structure, but empirical evidence of this is generally rare. Hybrid zones between two poeciliid fishes, Xiphophorus birchmanni and X. malinche, occur in multiple tributaries with independent replication of upstream‐to‐downstream gradients in morphology and allele frequencies. Ecological niche modelling revealed that temperature is a central predictive factor in the spatial distribution of pure parental species and their hybrids and explains spatial and temporal variation in the frequency of neutral genetic markers in hybrid populations. Among populations of parentals and hybrids, both thermal tolerance and heat‐shock protein expression vary strongly, indicating that spatial and temporal structure is likely driven by adaptation to local thermal environments. Therefore, hybrid zone structure is strongly influenced by interspecific differences in physiological mechanisms for coping with the thermal environment.     

Hybridization and mitochondrial genome introgression between Rana chensinensis and R. kukunoris

Mitochondrial genome (mito‐genome) introgression among metazoans is commonplace, and several biological processes may promote such introgression. We examined two proposed processes for the mito‐genome introgression between Rana chensinensis and R. kukunoris: natural hybridization and sex‐biased dispersal. We sampled 477 individuals from 28 sites in the potential hybrid zone in the western Tsinling Mountains. Mitochondrial gene (cyt‐b) trees were used to examine the introgression events. Microsatellite DNA loci, cyt‐b and morphological data were used to identify hybrids and to examine the extent of natural hybridization. We detected rampant bidirectional introgressions, both ancient and recent, between the two species. Furthermore, we found a wide hybrid zone, and frequent and asymmetric hybridization. The hybrid zone cline analysis revealed a clear mitochondrial–nuclear discordance; while most nuclear markers displayed similar and steep clines, cyt‐b had a displaced cline centre and a more gradual and wider cline. We also detected strong and asymmetric historical maternal gene flow across the hybrid zone. This widespread hybridization and detected low mito‐nuclear conflicts may, at least partially, explain the high frequency of introgression. Lastly, microsatellite data and population genetic methods were used to assess sex‐biased dispersal. A weak pattern of female‐biased dispersal was detected in both species, suggesting it may not play an important role in the observed introgression. Our data are consistent with the hybridization hypothesis, but support for the sex‐biased dispersal hypothesis is weak. We further suggest that selective advantages of the R. kukunoris‐type mito‐genome in thermal adaptation may also contribute to the introgression between the two species.     

Diversification of host use in two polyphagous butterflies: differences in oviposition specificity or host rank hierarchy?

Novel host usage may represent an initial step towards diversification or radiation onto novel hosts within an evolutionary lineage, particularly if a shift in host plant preference ranking takes place. Polyphagous stages of evolutionary lineages may represent transitional states in which novel host associations are more likely to develop, but may be more difficult to detect experimentally. The polyphagous sister species Papilio glaucus L. and Papilio canadensis (Lepidoptera: Papilionidae; these Papilio = Pterourus) are known to exhibit differences in host‐plant use, despite significant overlap in host‐use abilities, providing an opportunity to examine how host shifts in polyphagous species may occur and what the implications for future divergence may be. In particular, we were interested in (i) determining whether differences in oviposition behavior of these species were due to changes in specificity or shifts in host‐plant hierarchy, (ii) whether the varying preference for primary hosts also affected the preference for secondary hosts, and (iii) what the oviposition preferences of a new hybrid swarm population are. We examined more than 40 000 oviposition bouts from more than 400 P. glaucus, P. canadensis, and hybrid females placed in seven‐, three‐, or two‐choice assays. In each of the choice assays, leaves from plants in different plant families of varying suitability for P. glaucus and P. canadensis larvae were used. We found the primary difference between P. glaucus and P. canadensis to be limited to a Z‐linked shift in host rank hierarchy due to an acceptance of Populus tremuloides Michx. (Salicaceae) and reduced specificity for Liriodendron tulipifera L. (Magnoliaceae) in P. canadensis. In addition, we found the absence of the Z‐linked oviposition acceptance of P. tremuloides in a recently formed allochronically separated hybrid swarm population found in P. canadensis territory at the northern border of the P. glaucus and P. canadensis hybrid zone.     

THE PONTIA DAPLIDICE-ED USA HYBRID ZONE IN NORTHWESTERN ITALY

The pierid butterflies Pontia daplidice and P. edusa, parapatrically distributed in southern Europe, have very similar morphologies and life histories, but show fixed differences at four allozyme markers. We sampled these allozymes in a 28-population transect north of Genoa in Italy, through the hybrid zone where these taxa meet. We used the numerical techniques developed for hybrid zone analysis to study the patterns of genetic differentiation and their underlying evolutionary causes. The hybrid zone is characterized by a very short and steep central region, flanked by broad tails of introgression extended up to 100 km in either direction. From mean two-locus disequilibium of D = 0.148 (maximum-likelihood two-unit support limits 0.139-0.153), and after accounting for minor differences in the center locations of the single-locus clines, which act to bias the dispersal estimate, we estimated a dispersal rate of σ = 4.4 (3.7-5.5) km/gen^1/2. The effective selection needed to maintain the steep central portion is strong, 0.47 < s? < 0.64, when combined over potential intrinsic (genetic background) and extrinsic (ecological) sources of selection. The clines in allozyme loci showed variation that was significantly different between the most divergent shapes, and the differences are attributable to different degrees of introgression on the edusa side of the zone. The average selection acting on individual allozyme loci was high at s_???e  1.5%, but because of the narrowness of the central region of the cline, we suspect that this estimate is somewhat biased by selection on loci closely linked to the allozyme markers. A common question for taxa that show fixed allozyme differences in parapatry is whether or not they are genetically isolated. A fairly general measure of genetic isolation across hybrid zones is the time, T, that it takes a neutral allele to cross the hybrid zone and recombine into the opposite genetic background, given by T = (β/σ)², where β is the barrier strength of the hybrid zone. Genetic isolation in the Pontia zone is weak, with T  25 generations for most allozyme markers. By this measure, populations of daplidice and edusa on opposite sides of the hybrid zone share more identical-by-descent alleles than do populations of phenotypically pure daplidice in, say, France and Morocco. Accordingly, we think it best for systematists to consider edusa as a well-marked subspecies of P. daplidice.     

A phylogeographic investigation of the hybrid origin of a species of swordtail fish from Mexico

Hybrid speciation may contribute significantly to generating biodiversity, but only a few well‐documented examples for it exist so far that do not involve polyploidization as a mechanism. The swordtail fish, Xiphophorus clemenciae, shows common hallmarks of a hybrid origin and still overlaps in its current geographic distribution with its putative ancestral species (Xiphophorus hellerii and Xiphophorus maculatus). Xiphophorus clemenciae provides an ideal system for investigating the possible continued genetic interactions between a hybrid and its parental species. Here, we use microsatellite and mitochondrial markers to investigate the population structure of these species of swordtails and search for signs of recent hybridization. Individuals were sampled from 21 localities across the known range of X. clemenciae– the Isthmus of Tehuantepec (IT) Mexico, and several environmental parameters that might represent barriers to dispersal were recorded. The hybridization event that gave rise to X. clemenciae appears to be rather ancient, and a single origin is likely. We find negligible evidence for ongoing hybridization and introgression between the putative ancestral species, because they now occupy distinct ecological niches, and a common haplotype is shared by most populations of X. clemenciae. The population structure within these species shows an isolation‐by‐distance (IBD) pattern and genetic differentiation between most populations is significant and high. We infer that tectonic evolution in the Isthmus has greatly restricted gene flow between the southern and central IT populations of X. clemenciae and X. helleriii and provide preliminary information to aid in conservation management of this geographically restricted hybrid species, X. clemenciae.     

Hybridization and introgression in two ecologically dissimilar Fundulus hybrid zones

Hybridization and introgression appear more common in rapidly evolving groups, suggesting an important role in the evolutionary process. Detailed studies of how extrinsic or intrinsic forces regulate hybridization and introgression have the potential for broadening our understanding of mechanisms generating diversity. Species in the Fundulus notatus species complex have broad overlapping ranges and occur in replicated hybrid zones along predictable stream gradients. Typical hybrid zone structure has Fundulus olivaceus in headwaters, F. notatus downstream, and hybrid zones near confluences or abrupt shifts in habitat. Rarely, the typical upstream‐downstream orientation is reversed raising questions as to how hybrid zones are formed and maintained. We used next‐generation sequencing data to study hybridization and introgression in hybrid zones in neighboring drainages that differ in orientation (typical and reversed). We predicted extrinsic forces linked to stream gradients would result in noticeable differences between the two. Contrary to predictions, the data indicate the hybrid zones are remarkably similar. We used individual‐based simulations to explore the potential role of intrinsic and extrinsic forces in generating and maintaining typical and reversed hybrid zones. Simulation results were consistent with reversed hybrid zones being formed from stochastic processes combined with strong intrinsic forces and weak extrinsic forces.