Speciation in Western Scrub-Jays, Haldane’s rule, and genetic clines in secondary contact
© Gowen et al.; licensee BioMed Central Ltd. 2014
Received: 11 March 2014
Accepted: 4 June 2014
Published: 17 June 2014
Haldane’s Rule, the tendency for the heterogametic sex to show reduced fertility in hybrid crosses, can obscure the signal of gene flow in mtDNA between species where females are heterogametic. Therefore, it is important when studying speciation and species limits in female-heterogametic species like birds to assess the signature of gene flow in the nuclear genome as well. We studied introgression of microsatellites and mtDNA across a secondary contact zone between coastal and interior lineages of Western Scrub-Jays (Aphelocoma californica) to test for a signature of Haldane’s Rule: a narrower cline of introgression in mtDNA compared to nuclear markers.
Our initial phylogeographic analysis revealed that there is only one major area of contact between coastal and interior lineages and identified five genetic clusters with strong spatial structuring: Pacific Slope, Interior US, Edwards Plateau (Texas), Northern Mexico, and Southern Mexico. Consistent with predictions from Haldane’s Rule, mtDNA showed a narrower cline than nuclear markers across a transect through the hybrid zone. This result is not being driven by female-biased dispersal because neutral diffusion analysis, which included estimates of sex-specific dispersal rates, also showed less diffusion of mtDNA. Lineage-specific plumage traits were associated with nuclear genetic profiles for individuals in the hybrid zone, indicating that these differences are under genetic control.
This study adds to a growing list of studies that support predictions of Haldane’s Rule using cline analysis of multiple loci of differing inheritance modes, although alternate hypotheses like selection on different mtDNA types cannot be ruled out. That Haldane’s Rule appears to be operating in this system suggests a measure of reproductive isolation between the Pacific Slope and interior lineages. Based on a variety of evidence from the phenotype, ecology, and genetics, we recommend elevating three lineages to species level: A. californica (Pacific Slope); A. woodhouseii (Interior US plus Edwards Plateau plus Northern Mexico); A. sumichrasti (Southern Mexico). The distinctive Edwards Plateau population in Texas, which was monophyletic in mtDNA except for one individual, should be studied in greater detail given habitat threat.
KeywordsBirds Speciation Gene flow Postzygotic reproductive isolation Phylogeography
The study of speciation has been influenced during the last decade by increasing awareness that different parts of the genome tell different evolutionary stories [1–3]. Previously, speciation studies relied heavily on mitochondrial DNA (mtDNA), given the numerous benefits offered by these data for inferring evolutionary history at diverse timescales (e.g., high substitution rates) . Yet, along with these benefits come drawbacks , among them that mtDNA is a single recombinational unit that can present a biased snapshot of evolutionary history when mtDNA data are not balanced by independent evidence from nuclear markers.
One case in which mtDNA data can be particularly misleading occurs during gene flow between populations in secondary contact. Haldane was the first to note that when two differentiated lineages come back into reproductive contact, the heterogametic sex (i.e., the one with two types of sex chromosomes) often shows reduced fertility . The mechanistic reasons behind Haldane’s Rule, as it came to be known, are under intense study [7–10]. While we may not understand why Haldane’s Rule occurs, the consequences for female-heterogametic species like birds and butterflies are clear: if female hybrids show reduced fertility, then mtDNA, which is inherited as a single recombinational unit from the female parent, is less likely to be passed on to the next generation. As a result, mtDNA might show little or no evidence of hybridization, despite the fact that hybridization occurs and paternal, nuclear alleles move between lineages.
The potential mismatch between gene flow estimated from mtDNA and gene flow estimated using nuclear markers is not just theoretical; it has been demonstrated empirically by studies of female-heterogametic organisms [11–18]. This phenomenon is of interest to biologists, among other reasons, because of the continuing influence of the Biological Species Concept (BSC) and its emphasis on cessation of gene flow as a necessary criterion for species delimitation . While the BSC may not be applied as commonly in practice as is claimed , it is nevertheless evident that many attempts to split species taxonomically fail due to the existence of gene flow and the perceived importance of the BSC. As empirical evidence continues to mount demonstrating mito-nuclear discordance [21, 22], the use of nuclear markers has become essential to species delimitation and the study of speciation.
These two lineages (‘coastal’ and ‘interior’) are thought to have diverged from each other anywhere from 1 to 4 million years ago (Ma) based on divergence times calibrated with a fossil and estimated molecular substitution rates . Coastal individuals are easily distinguished from those in the interior by their brighter blue plumage, bolder white supercilium, and more pronounced blue collar. Coastal individuals also have stout, hooked bills, which are better for husking acorns than the smaller, more pointed bills of interior individuals, which are better suited for extracting pine seeds from cones [27, 28].
The hybrid zone was originally described by Pitelka  who traversed the zone and analyzed numerous museum specimens. He found intermingling of coastal and interior phenotypes in western Nevada, as well as individuals that appeared to be F1 hybrids. The hybrid zone has never been studied in detail using genetic techniques, although previous work documented a relatively rapid turnover of mtDNA types and mismatches of taxonomy and genotype of birds from this zone [23, 24]. Nevertheless, a prior proposal to elevate the coastal and interior lineages to species status  failed because of doubts concerning the extent of gene flow in the contact zone.
The purpose of this study was to test for a prediction of Haldane’s Rule that there will be a mismatch between mtDNA and nuclear gene flow across the contact zone. Specifically, we predict that if Haldane’s Rule is operating in this system, we will see a narrower (i.e., steeper) cline in the transition of mtDNA haplotypes compared to nuclear markers. Before addressing this hypothesis, however, we needed to determine the exact nature of the contact zone in western Nevada and whether there were any undescribed contact zones nearby that might confound our results. We therefore first assessed broad-scale patterns of divergence and gene flow among Western Scrub-Jays across their entire geographic range, using deep sampling of populations and a suite of nuclear markers. Previous molecular studies sampled only one or a few individuals of different populations, focusing largely on allozymes  or mtDNA . Our study design includes nearly 700 individuals sampled from throughout the range of the species, thus providing a comprehensive portrait of genetic divergence and allowing us to locate potentially undescribed contact zones. Finally, we correlated proportion of hybrid ancestry with phenotypic traits for a subset of individuals in the contact zone, as a way to control for the influence of environment and therefore test for a genetic basis of phenotypic differences between the coastal and interior lineages. If phenotypic traits are genetically based, we expect to see a correlation between neutral genetic assignment to coastal or interior lineages and the phenotypic differences known to characterize those lineages.
Geographic structure of coastal and interior mtDNA types
To assess gene flow throughout the range of Western Scrub-Jays, we conducted simple mtDNA typing of all 689 individuals using restriction digest, which showed that most populations are composed entirely of coastal or interior mtDNA types, with very few populations showing mixing of mtDNA types. The four populations showing mixing were (1) the previously described contact zone in western Nevada, which featured several populations having mixed types 100 km east and southeast of Reno and Carson City (sites inside the inset square in Figure 1); (2) a single interior type (among n = 16 coastal types) found at the base of the western Sierra 100 km east of Stockton, CA (site 6); (3) an even mix of coastal and interior types in the eastern Sierra Nevada to the west of the Owens Valley near Independence, CA (site 18); and (4) one remarkable interior type (among n = 16 coastal types) in the Santa Ana Mountains in Orange County, Los Angeles (site 23).
Broad-scale geographic structure of microsatellite variation
To complement the mtDNA typing results with information from nuclear markers, we genotyped all 689 individuals at 14 microsatellite loci. No pairs of loci were physically linked based on linkage tests. Of the 14 loci, one locus (MJG8) showed significant deviations from HWE across many populations following Bonferroni correction. We removed this locus from the data set. Four other loci showed deviations from HWE across some populations. We ran downstream analyses with and without these loci to investigate whether they were driving any of the patterns of genetic structure we observed. Results among analyses including and excluding loci were qualitatively similar, so we retained the larger group of 13 loci. No loci showed consistent signs of other microsatellite artifacts like large-allele drop-out or null alleles.
Phylogeny from mtDNA
Geographic cline analysis
Correlation between phenotype and genotype
Regression ( R 2 ) values between hybrid index and phenotypic traits within and outside the hybrid zone
Haldane’s Rule and the mismatch in cline width between mtDNA and microsatellites
Our results support a key prediction from Haldane’s Rule that mtDNA will show a steeper cline across a hybrid zone than nuclear markers. Although Haldane’s Rule may explain the clinal pattern we observed, this pattern could also result from divergent selection on mtDNA types  or sex-biased dispersal, wherein females (and the mitochondrial genomes they pass on) tend to disperse less than males. Although we have not looked at selection, what little is known about natal dispersal in Western Scrub-Jays suggests that, to the contrary, females disperse about three times farther than males on average . The results from our time-since-contact analyses, which incorporated estimates of dispersal distances, also support predictions from Haldane’s Rule, with less diffusion of mtDNA compared to nuclear markers.
Our results showing steeper clines in mtDNA compared to nuclear markers add to a gathering corpus of research showing discordant rates of DNA transmission across hybrid zones among different marker types [11–18, 22], a pattern that has often been attributed to Haldane’s Rule . Whatever the root cause, these results caution that conclusions about individual-level gene flow made on the basis of mtDNA alone may be underestimated. In our study, gene flow is generally restricted in both mtDNA and nuclear markers, due to the low dispersal distances of Western Scrub-Jays. In other species, such as long-distance dispersers with no sex-biased dispersal, mismatches could be much more pronounced.
Recent work on the genetic basis for Haldane’s Rule suggests it results from incompatibilities among genes that are exposed to selection on hybrid female Z chromosomes, causing sterility , whereas only those incompatibilities that are dominant are visible to selection in hybrid ZZ males. If Haldane’s Rule is indeed the cause of the clinal mismatch we observed, then our results suggest that the coastal and interior lineages of Western Scrub-Jay have begun the process of reproductive isolation. Ultimately, Haldane’s Rule as the cause of the mismatch in marker introgression rates could be determined through field studies or captive mating experiments that document the reproductive viability of female offspring resulting from hybrid crosses .
Clarifying the hybrid zone between coastal and interior Western Scrub-Jays
Our results also help clarify the nature of the hybrid zone between coastal and interior lineages of Western Scrub-Jays. The fact that our hybrid transect populations could be fit to a cline suggests that Pitelka  was correct that the major axis of gene flow is from the northwest to southeast, as coastal populations continue around the north side of the Sierra Nevada and meet with interior populations in western Nevada. However, the coastal influence of the Woodford site (45 in Figure 2) warrants further study. In addition to gene flow from the north, some gene flow may link the southern end of the Sierra Nevada to the Woodford site, with coastal individuals inhabiting canyons of the southeastern Sierra that were not sampled for this study.
Our results also support Pitelka’s  conclusions that the hybrid zone is narrow and the result of secondary contact. Even the oldest time estimates from our time-since-contact analysis across the hybrid zone (~437,000 years) are younger than divergence estimates of the original split between coastal and interior lineages [1 to 4 Ma; 24], supporting the hypothesis that the hybrid zone results from secondary contact and not in situ divergence. In terms of broader biogeographic evidence, several other species have secondary contact zones where coastal and interior lineages meet in this general region of the Great Basin [34, 35].
Our results also confirm that the contact zone in western Nevada is the only place where coastal and interior lineages come into contact to any great degree. Outside of western Nevada, two sites (6 and 23 in Figure 1) appear to contain examples of rare, long-distance dispersers from pure interior populations because the nuclear DNA profile of the two individuals matches their mtDNA type (the interior mtDNA individual at site 6 was 97% interior in nuclear DNA, and the interior mtDNA individual at site 23 was 95% interior in nuclear DNA). In contrast, at a third site in the eastern Sierra Nevada (18 in Figure 1), the three individuals with interior mtDNA had strongly coastal nuclear DNA profiles (91%, 95%, and 89% coastal), suggesting that they were advanced backcrosses from older hybridization events (i.e., nuclear DNA profiles of F1s would be closer to 50/50). Thus, although there appears to be one major contact zone, the exceptions show that the reality of dispersal and gene flow among jays in and around the Sierras – and across the Mojave Desert farther south – is more complex. Fine-scale sampling, combined with habitat-based methods of identifying dispersal corridors , could provide a more complete picture of gene flow across and around the Sierra Nevada for this and other species.
Previous studies documented phenotypic differences between coastal and interior Western Scrub-Jays [25, 28], but ours is the first to provide evidence for a genetic basis for some traits. There was a correlation between nuclear genetic variation and phenotypic traits, such as bill shape and plumage traits, when considering all individuals, even those from pure populations (Table 1). However, including individuals from different geographic locations potentially confounds environmental and genetic sources of phenotypic variation. Still, many of these relationships remained, for the plumage traits at least, within the contact zone where there are hybrid individuals and all variety of backcrosses. The correlation between proportion of interior vs. coastal genetic ancestry and phenotypic traits for individuals that were fledged in the hybrid zone and therefore experienced the same environment since fledging provides evidence for a genetic basis of plumage traits, at the very least. Why we did not observe similar relationships for bill traits is not clear, but could be due to low statistical power, stabilizing selection, or greater plasticity in bill traits. Further studies with more individuals and better genomic coverage will help discriminate between these alternate hypotheses, as will studies that attempt to identify the genetic loci or environmental pressures responsible for the phenotypic differences.
Genetic differentiation across North America’s scrub woodland
Previous genetic studies of species inhabiting North America’s pine-oak and scrub woodland have uncovered unexpectedly deep phylogenetic structure within species [37, 38]. Results of the current study, which are the most comprehensive in terms of geographic, individual, and genetic sampling of Western Scrub-Jays, indicate at least five genetic clusters with strong spatial structuring (colored dots in Figure 1). The Pacific Coast and Southern Mexico groups have been recovered previously in studies of phenotype, allozymes, and mtDNA [23–26]. The Edwards Plateau group and the split between Interior US and Northern Mexico individuals are novel, although the inferred population clusters align with current subspecies designations .
This study is the first to show clear evidence of nuclear and mtDNA differentiation of Western Scrub-Jays on the Edwards Plateau in Texas (currently a separate subspecies, A. c. texana). All individuals in the Edwards Plateau population except one possess a unique mtDNA haplotype differing from other interior US haplotypes by two substitutions. Nuclear and mtDNA differentiation on the Edwards Plateau is interesting given the history of isolation of this region and its unique flora and fauna. The oak savanna and Ashe juniper woodland of the Edwards Plateau are separated from other scrub woodland to the west by an extension of the Chihuahuan Desert in the Pecos Valley. Fossil and pollen data suggest, however, that pygmy woodland was continuous across this region, and connected the Edwards Plateau to western forests throughout the last glacial maximum until about 11,000 years ago . At least one other distinct bird species with a similar geographic range occurs on the Edwards Plateau, the Golden-cheeked Warbler (Setophaga chrysoparia), and other organisms in this region have distinct populations or subspecies [40, 41]. Our results lend additional support to the distinctiveness of the biota in this region, where native habitat is under persistent conservation threat from urbanization.
The nuclear genetic break between the Interior US and Northern Mexico populations, occurring roughly where the Sierra Madre Occidental and Oriental turn into a network of sky islands near the US border, could be explained either by a break in suitable habitat or a sampling artifact. Addressing the latter point, observations posted on eBird (ebird.org) confirm observations of Western Scrub-Jays in southeastern Arizona and northern Sonora (in the vicinity of the genetic break), but these records could represent wandering individuals. Likewise few observations exist in the vicinity of the genetic break in southwestern Texas, and virtually no observations exist from northern Coahuila, where congeneric Mexican Jays (Aphelocoma wollweberi) occur at lower elevations, seemingly occupying the Western Scrub-Jay niche [25, 42–44]. Thus, while available evidence suggests that the genetic break between Interior US and Mexico populations relates to a real distributional break near the US border, the exact causes of this break remain unclear. Future niche modeling could be used to determine if there are cryptic breaks in suitable habitat.
Lack of genetic differentiation along the Pacific Coast
In contrast to the highly structured genetic portrait of the interior US and Mexico, Western Scrub-Jays along the Pacific Slope show a remarkable lack of genetic structure despite 2,500 km of nearly unbroken sampling from northern Oregon to the southern tip of Baja California. Lack of genetic structure over large distances suggests two possible scenarios: (1) high levels of dispersal and gene flow link coastal populations; or (2) recently expanded coastal populations have had insufficient time to accumulate genetic differences. Available evidence suggests that Western Scrub-Jays generally have low dispersal , notwithstanding a few rare, long-distance dispersal events discussed above. Furthermore, the original split between coastal and interior lineages occurred between 1–4 million years ago , providing ample time for genetic structure to emerge in coastal populations, absent other demographic events. Thus, the hypothesis that seems most likely is a recent expansion following a population bottleneck. However, Pleistocene niche models of Western Scrub-Jay niches suggest ample habitat along the Pacific Slope during the last glacial maximum . Further study is needed to determine the cause of the apparent genetic homogeneity of Pacific Slope populations.
We recommend splitting the current concept of the Western Scrub-Jay into three species, as has already been suggested previously . The Pacific Slope lineage (A. californica, including subspecies californica, oocleptica, caurina, obscura, hypoleuca, superciliosa, immanis, and cactophila) is phenotypically distinct , possesses a bill morphology adapted to local resources [27, 28], is monophyletic in mtDNA, and is sister to another recognized species, the Island Scrub-Jay [23, 24]. This study shows that the Pacific Slope lineage is also well-differentiated in nuclear markers and has only a narrow zone of contact with interior populations. The genetic clines across the hybrid zone are generally steep, and the mismatch between mtDNA and nuclear clines suggests a measure of reproductive isolation. These results would appear to meet criteria for species distinction under even the most stringent species concept.
Southern Mexico populations (A. sumichrasti, including subspecies sumichrasti and remota) also meet criteria for species recognition. They are phenotypically distinct , with behavioral evidence that some southern populations are cooperative breeders , whereas other Western Scrub-Jays breed in territorial pairs. This study demonstrates that these southern populations are reciprocally monophyletic to other interior populations in mtDNA and well-differentiated in nuclear markers.
Finally, we recommend that the remaining interior populations (A. woodhouseii including subspecies texana, woodhouseii, nevadae, grisea, and cyanotis) be elevated to species level, given their diagnostic phenotype , adaptive bill morphology , high support for monophyly in mtDNA in previous work [23, 24] and this study (0.94 PP), and differentiation in nuclear markers (this study). Within this group, we do not recommend elevating the genetically distinctive Edwards Plateau population to species level at this time because our study sampled only one population. More information is needed on potential gene flow with other interior US populations, as well as a modern multivariate analysis of their phenotypic differentiation and studies of their ecological and behavioral divergence. Likewise, we do not currently recommend splitting Interior US from Northern Mexico populations, unless, at minimum, further sampling near the US-Mexico confirms that the break seen in microsatellite variation is not the result of a sampling artifact.
We report a mismatch in cline widths between mtDNA and nuclear markers across a hybrid zone in Western Scrub-Jays. This result is consistent with predictions based on Haldane’s Rule, where there is less mtDNA introgression between lineages in female-heterogametic species. This study adds to a growing list of studies that support predictions of Haldane’s Rule using cline analysis of multiple loci of differing inheritance modes. Structure results of microsatellite variation confirm that there is only one major area of contact between coastal and interior lineages. Furthermore, Structure identified at least five clusters of individuals with strong spatial structuring: Pacific Slope, Interior US, Edwards Plateau, Northern Mexico, and Southern Mexico. Based on a variety of evidence from the phenotype, ecology, and genetics, we recommend elevating three lineages to species level: A. californica (Pacific Slope); A. woodhouseii (Interior US plus Edwards Plateau plus Northern Mexico); A. sumichrasti (Southern Mexico). The distinctive Edwards Plateau population in Texas, which was also monophyletic in mtDNA except for one individual, should be studied in greater detail given conservation concerns.
Specimens and DNA extraction
Locality information for Western Scrub-Jays sampled for this study
Vida, Oregon, USA
Alturas, California, USA
Lassen Co., California, USA
Douglas City, California, USA
Beckwourth, California, USA
Poole Station Rd, Calaveras Co., California, USA
Martinez, Contra Costa Co., California, USA
El Sobrante, Contra Costa Co., California, USA
Alhambra Valley Rd., Contra Costa Co., California, USA
Berkeley, California, USA
Carmel Valley Rd., Alameda Co., California, USA
Moraga 1, Contra Costa Co., California, USA
Moraga 2, Contra Costa Co., California, USA
Moraga 3, Contra Costa Co., California, USA
Orinda, Contra Costa Co., California, USA
Danville, Contra Costa, Co., California, USA
San Ramon, Contra Costa Co., California, USA
Dublin, Alameda Co., California, USA
Purisima Creek, San Mateo Co., California, USA
North Fork, Madera Co., California, USA
Independence, Inyo Co., California, USA
Bradley, Monterey Co., California, USA
Bodfish, Kern Co., California, USA
Big Bear City, San Bernadino Co., California, USA
Pioneertown, San Bernadino Co., California, USA
Santa Ana Mountains, Orange Co., California, USA
Love Valley, San Diego Co., California, USA
Pinyon Flat, Riverside Co., California, USA
Vallecito Mountains, San Bernadino Co., California, USA
Mountain Springs Pass, Imperial Co., California, USA
Jacumba, San Diego Co., California, USA
La Rumorosa, Baja California, Mexico
Potrero, San Diego Co., California, USA
La Rosa de Castilla, Baja California, Mexico
San Lucas, Baja California Sur, Mexico
Bahia Magdalena, Baja California Sur, Mexico
La Burrera, Baja California Sur, Mexico
Virginia Mountains, Storey Co., Nevada, USA
superciliosa x nevadae
Lousetown, Storey Co., Nevada, USA
superciliosa x nevadae
Brunswick Cyn, Pine Nut Mountains, Douglas Co., Nevada, USA
superciliosa x nevadae
Sunrise Pass Rd. 1, Pine Nut Mountains, Douglas Co., Nevada, USA
superciliosa x nevadae
Sunrise Pass Rd. 2, Pine Nut Mountains, Douglas Co., Nevada, USA
superciliosa x nevadae
Lebo Spring, Pine Nut Mountains, Douglas Co., Nevada, USA
superciliosa x nevadae
Sunrise Pass Rd. 3, Pine Nut Mountains, Douglas Co., Nevada, USA
superciliosa x nevadae
Hot Springs Mountain, Douglas Co., Nevada, USA
superciliosa x nevadae
Pine Nut Creek, Douglas Co., Nevada, USA
superciliosa x nevadae
Gardnerville, Douglas Co., Nevada, USA
superciliosa x nevadae
Woodford, Alpine Co., California, USA
superciliosa x nevadae
Nye Canyon, Pine Grove Hills, Lyon Co., Nevada, USA
Pine Grove Hills, Lyon Co., Nevada, USA
Lucky Boy Pass Rd, Wassuk Range, Mineral Co., Nevada, USA
North Canyon, Wassuk Range, Mineral Co., Nevada, USA
Excelsior Mountains, Mineral Co., Nevada, USA
White Mountains, Inyo Co., California, USA
Mt Charleston, Clark Co., Nevada, USA
Toiyabe Mountains, Lander Co., Nevada, USA
Baker, Millard Co., Utah, USA
Hualapi Mountains, Mohave Co., Arizona, USA
Drake, Yavapai Co., Arizona, USA
Stansbury Mountains, Tooele Co., Utah, USA
Uinta Mountains, Duschesne Co., Utah, USA
Duschesne Co., Utah, USA
Bonanza, Uitah Co., Utah, USA
Whitewater, Mesa Co., Colorado, USA
Monument, Teller Co., Colorado, USA
Gardner, Huerfano Co., Colorado, USA
Manzano, Valencia Co., New Mexico, USA
Fort Davis, Jeff Davis Co., Texas, USA
Carta Valley, Edwards Co., Texas, USA
Villa Ocampo, Durango, Mexico
Sombrerete, Zacatecas, Mexico
Rancho Santa Rita, Jalisco, Mexico
El Diamante Pass, Coahuila, Mexico
Bledos, San Luis Potosí, Mexico
Taxco, Guerrero, Mexico
Xocomanantlán, Guerrero, Mexico
San Lorenzo de Abarrados, Oaxaca, Mexico
Assigning coastal vs. interior mtDNA haplotypes
Phylogenies based on mtDNA have been characterized in previous studies of Western Scrub-Jays using a limited numbers of individuals [23, 24]. Our goal was to conduct deep mtDNA sampling across the range of the species and assign coastal and interior types to determine the extent of mixing within populations, and to determine whether rare long-distance dispersal events have occurred.
Instead of direct sequencing, we screened individuals for coastal versus interior mtDNA type with a rapid, diagnostic assay. We used previously sequenced individuals  to identify a rare-cutting restriction enzyme that differed between coastal individuals, which had the cut site, and interior individuals, which did not. Then, for new samples of unknown haplotype, we amplified the cyt b gene using previously published primers and PCR conditions . Once completed, we added 0.5 μL of BSR-DI (2000 U/mL, New England Biolabs, Ipswitch, MA) to each reaction and incubated for one hour at 65°C, followed by visualization on an agarose gel. We scored individuals having two bands as coastal and individuals having one band as interior, employing both positive and negative controls. We ran samples where mtDNA type conflicted with geographic (sampling) location a second time to confirm the result.
Nuclear microsatellite loci
We used microsatellites to gain insight into differentiation and gene flow of the nuclear genome and to score individuals as potential hybrids (given that mtDNA is haploid and provides no information as to the hybrid ancestry of individuals). Previously characterized microsatellite loci have been isolated from the Florida Scrub-Jay . From these loci, we chose a panel of 14 loci that amplified well and were variable among Western Scrub-Jays. To reduce processing time and cost while maintaining data quality, we pooled up to three loci and amplified them together (multiplexing) with a single dye , after first confirming that allele size classes for pooled loci did not overlap.
After amplification, we determined fragment size using an ABI PRISM 3730 capillary sequencer and analyzed the resulting data using Geneious v. 6.0.4 (Biomatters). We called alleles objectively by creating bins based on known microsatellite repeat motifs , and assessed artifacts after genotyping (e.g., stuttering, null alleles, big-allele drop-out) using Microchecker . To ensure that microsatellite loci were evolving under neutral processes, we assessed each locus for deviation from Hardy-Weinberg equilibrium (HWE) using GenePop on the Web (http://genepop.curtin.edu.au). Strong population structure can lead to false positives for deviations from HWE when populations are pooled (Wahlund Effect), so we assessed HWE in each population at each locus individually, applying Bonferroni correction to control for false positives resulting from the effect of conducting multiple simultaneous tests. We also performed a test of linkage disequilibrium to confirm independent sorting of loci.
Analysis of microsatellite variation for genetic clustering across the geographic range
We assessed genetic structure of microsatellite data using Structure v. 2.3.4 , which uses Bayesian analysis to infer the number of populations (K) from a group of individuals based only on their genetic variation, without prior information on where the individuals originated. We did not use the method of Evanno et al.  to detect the ‘true K’ because this method is known to underestimate genetic structure in all but cases of very strong genetic differentiation . Initial tests of the Evanno et al.  method suggested K = 2 as the best model for our data; however, analysis at higher K revealed strong structure that was highly justified based on biology and geography. Thus, because one of our goals was discovery of the basic units of spatially structured genetic diversity, we instead used Structure to analyze successively smaller clusters of individuals at K = 2 until Structure could no longer find any clustering within groups, as demonstrated previously . For example, when two clear geographic population clusters with strong biological justification emerged in the complete data set, we ran each of these clusters in its own analysis until K = 2 revealed no remaining geographical clusters. After the first run, we excluded hybrid populations (identified as populations with both coastal and interior mtDNA types) from the analysis. We ran Structure with an admixture model, correlated allele frequencies, and a burn-in period set to 100,000 generations, followed by 500,000 generations, which was sufficient for each run to reach stationarity.
Analysis of microsatellite variation across the hybrid zone
To assess population assignment and potential hybrid ancestry of individuals across the hybrid zone, we delimited 13 populations along a transect running northwest to southeast from the Sierra Nevada in northeastern California into west-central Nevada (Figure 2). This line conforms to prior descriptions of the contact zone [25, 55, 56] and our own field observations. To obtain estimates of hybrid ancestry (i.e., Q-scores) within the hybrid zone, individuals from these 13 populations were analyzed in Structure separately from all other individuals in the study using K = 2; their resulting Q-scores were recorded for later use in cline analysis. The Q-score values provide an overall estimate of nuclear variation that is expected to suffer less from random processes affecting individual loci. Moreover, cline analysis requires binary data or frequency values varying from 0 to 1. Due to large number of alleles at a given locus, individual microsatellite loci do not usually conform to this standard.
Testing for phylogenetic structure of microsatellite groups using mtDNA
Previous studies that generated mtDNA phylogenies of Western Scrub-Jays did not sample deeply within populations [23, 24]. Thus, to determine if nuclear DNA groups discovered in this study were discernible using mtDNA, we sequenced multiple individuals from all major groups suggested by our microsatellite results for the mtDNA cyt b gene, using previous protocols . We combined these new data with existing data for Aphelocoma and close outgroups (Calocitta formosa, Gymnorhinus cyanocephalus, and Cyanocitta stelleri) from GenBank, and we generated a phylogeny from these sequences using BEAST v. 1.7.5 . We ran the analysis using the best-fit substitution model for the cyt b gene determined through Akaike Information Criterion (AIC), which was HKY + I + G with six gamma categories and base frequencies estimated as determined by MrModeltest, v2 . We used a strict clock with a random starting tree. We ran the chain to 1.0 × 106 iterations, sampling every 1,000. The analysis converged quickly as determined by observing the plot of likelihood scores and ESS scores > 200 with Tracer. We discarded the first 1,000 of 10,000 trees as burn-in and produced a consensus tree from the remaining posterior trees using Geneious . The maximum clade credibility produced with TreeAnnotator was almost identical, with key nodes mentioned in the Results having the same PP values.
Fitting geographic clines to genetic data across the hybrid zone
We used HZAR v2.5 , a statistical package implemented in R, to determine the width of geographic clines estimated from mtDNA and nuclear markers, and the extent of mismatch, if any. We measured the distance between points from the coastal northwest end of the hybrid zone (site 2 in Figure 1), through the contact zone, continuing to pure interior populations to the southeast (site 52 in Figure 1). These were then compressed to a single line using HZAR. This transect likely captures the major axis of gene flow across the hybrid zone , with minor caveats covered in the Discussion. We estimated clines for mtDNA haplotype frequency and average Q-score per site as determined in our Structure analyses. We used these cline analyses to estimate changes of the molecular characters in local mean frequencies. We modeled the cline shape using three equations [61, 62] describing a sigmoid shape at the center of the transition with two exponential decay curves on either side of the transition. We estimated several parameters, including width (w), center (c), delta (d, distance between the center and the tail), and tau (t, slope of the tail). We also incorporated the possibility that Pmax and Pmin (the “top” and “bottom” of the cline) were either fixed or free to vary. We fit three sets of five cline models using the Metropolis-Hasting algorithm in R. Model set 1 has no scaling (Pmin = 0, Pmax =1), model set 2 has fixed scaling (Pmin = observed minimum, Pmax = observed maximum), and model set 3 allows Pmin and Pmax to vary. Within each model set, scaling and tails are fixed or free to vary. We compared these models to a null model of no clinal transition.
For the two genetic datasets, we set up a covariance matrix by running each model for 1.0 × 106 generations. We ran three independent chains for 9.0 × 106 generations and assessed for convergence and stability by visualizing the MCMC traces. We performed model selection using corrected AIC (AICc). We estimated two log-likelihood confidence intervals around each parameter estimate. If the parameter estimate from one dataset did not overlap with the confidence interval of the parameter estimate from the other dataset, we considered them to be significantly different. We also used estimates of cline width to determine the time since contact under a neutral diffusion model solving for t in the equation w = 2.51σ√t, where t is the time since contact and σ is the root mean square dispersal distance from estimates in . Because females and males have different dispersal distances, we used the female dispersal distance (1.5 – 5 km) to estimate time since contact using mtDNA, and we used the average dispersal estimate of males and females (0.5 – 3.0) to estimate time since contact for the nuclear dataset.
Correlating phenotype with hybrid ancestry
We collected phenotypic data on a subset of 124 coastal and interior Western Scrub-Jay specimens within and outside the hybrid zone. Unfortunately we could not perform a cline analysis along the transect because we did not have sufficient sample sizes from pure populations at both ends of the transect. However, to determine the extent to which phenotype is correlated with nuclear DNA variation, we conducted regression analysis of phenotypic traits (plumage and morphology) with Q-scores at several spatial scales. The set of 124 specimens includes the full ranges of coastal and interior phenotypes, many hybrids, and pure individuals of each lineage with no evidence for hybrid genetic ancestry (Q-score at or close to 1 or 0). Because environment can influence phenotype when considering allopatric lineages, a more robust test of the genotype-phenotype link occurs within the hybrid zone where individuals share the same environment and there is more continuous variation in Q-scores. Here, we looked specifically at the 68 specimens that are part of the hybrid transect, as well as a more restrictive sample of 50 individuals within the core hybrid zone defined as sites 35–45 in Figure 1.
To ensure consistent measurements, one author (TCW) measured wing, tail, tarsus, bill length, bill depth, and bill width to the nearest 0.1 mm with digital calipers on all 124 birds after first verifying high repeatability scores for all traits. On a subset of 66 of these birds, another author (JEM) assessed three qualitative plumage traits known to vary between coastal and interior lineages: the amount of blue collar (reduced in interior birds); size of the eyestripe (reduced in interior birds); and color of the undertail coverts (blue-tinged in interior birds as opposed to white in coastal birds). We determined variation in plumage traits by lining up all 66 specimens at the same time and arranging them in order of trait appearance (e.g., very white undertail coverts to very blue undertail coverts). We then divided the arranged specimens into six even categories and assigned a value of 1 to 6. We repeated these steps for each trait. We analyzed morphological data with principal components analysis on the correlation matrix in Stata 10 to identify axes defining the most variation in the data. We assessed relationships between Q-scores and univariate and multivariate phenotypic traits with linear regression at each of the three spatial extents described above (full range, hybrid transect, core hybrid zone) using Stata.
Availability of supporting data
New DNA sequences have been deposited in Genbank under accession numbers KJ835799-KJ835861. Museum catalog numbers, localities, raw microsatellite data, phenotypic data, and Q-scores for each individual are available through Dryad http://dx.doi.org/10.5061/dryad.57f48.
FCG conducted part of this research for her Master’s thesis at the Moore Laboratory of Zoology at Occidental College; JMM is the Collections Manager at the Moore Laboratory of Zoology; CC is Staff Curator of Birds at the Museum of Vertebrate Zoology at Berkeley; ATP is Distinguished Professor in the Ecology and Evolutionary Biology Department and Curator of Ornithology at the Biodiversity Institute at Kansas University; BCF is an Assistant Research Scientist at UCLA; TCW was an undergraduate at Louisiana State University while conducting this research; JEM is Curator of Birds and Mammals and Director of the Moore Laboratory of Zoology and an Assistant Professor in the Biology Department at Occidental College.
We thank Hemani Wijesuriya and Uma Dandekar at the UCLA Sequencing Core. Nathan Jackson helped set up the microsatellite protocol at LSU. Liz Derryberry helped with the sequencer at LSU. Many Occidental College undergraduates helped collect microsatellite data: Dan Crowley, Elise Hanson, Andrew Hurley, Duncan Marks, Olivia Reed, Sean Sazaki, Leela Thapa, Rebecca Tribelhorn, and Christina Turner. Dave Willard at the Field Museum of Natural History issued a massive tissue loan and, along with Steve Cardiff at LSU, assisted with study skin loans. Adolfo Navarro at Universidad Nacional Autónoma de México facilitated field collecting. Whitney Tsai assisted with the lab work and field work. We thank John Klicka for assistance with obtaining specimens in Nevada. We thank two reviewers and Frank Gill for comments on the manuscript. The illustration of the coastal Scrub-Jay in Figure 1 was modified from a photo by Dawn Beattie provided under Creative Commons license (https://creativecommons.org/licenses/by/2.0/legalcode). The illustration of the interior Scrub-Jay was modified from a photo ©www.briansmallphoto.com and has been used here with permission. This project was partially funded by an NSF EAGER grant DEB-1244739 to JEM, CC, and BCF and the Borestone Fund.
- Edwards SV: Is a new and general theory of molecular systematics emerging?. Evolution. 2009, 63: 1-19. 10.1111/j.1558-5646.2008.00549.x.PubMedView Article
- Maddison WP: Gene trees in species trees. Syst Biol. 1997, 46: 523-536. 10.1093/sysbio/46.3.523.View Article
- Wu C-I: The genic view of the process of speciation. J Evol Biol. 2001, 14: 851-865. 10.1046/j.1420-9101.2001.00335.x.View Article
- Zink RM, Barrowclough GF: Mitochondrial DNA under siege in avian phylogeography. Mol Ecol. 2008, 17: 2107-2121. 10.1111/j.1365-294X.2008.03737.x.PubMedView Article
- Ballard JWO, Whitlock MC: The incomplete natural history of mitochondria. Mol Ecol. 2004, 13: 729-744. 10.1046/j.1365-294X.2003.02063.x.PubMedView Article
- Haldane JBS: Sex ratio and unisexual sterility in hybrid animals. J Genet. 1922, 12: 101-109. 10.1007/BF02983075.View Article
- Schilthuizen M, Giesbers M, Beukeboom L: Haldane’s rule in the 21st century. Heredity. 2011, 107: 95-102. 10.1038/hdy.2010.170.PubMedPubMed CentralView Article
- Wu C-I, Johnson NA, Palopoli MF: Haldane’s rule and its legacy: why are there so many sterile males?. Trends Ecol Evol. 1996, 11: 281-284. 10.1016/0169-5347(96)10033-1.PubMedView Article
- Orr HA: Haldane’s rule. Annu Rev Ecol Syst. 1997, 28: 195-218. 10.1146/annurev.ecolsys.28.1.195.View Article
- Turelli M, Orr HA: The dominance theory of Haldane’s rule. Genetics. 1995, 140: 389-402.PubMedPubMed Central
- Carling MD, Brumfield RT: Haldane’s rule in an avian system: using cline theory and divergence population genetics to test for differential introgression of mitochondrial, autosomal, and sex-linked loci across the Passerina bunting hybrid zone. Evolution. 2008, 62: 2600-2615. 10.1111/j.1558-5646.2008.00477.x.PubMedView Article
- Naisbit RE, Jiggins CD, Linares M, Salazar C, Mallet J: Hybrid sterility, Haldane’s rule and speciation in Heliconius cydno and H melpomene. Genetics. 2002, 161: 1517-1526.PubMedPubMed Central
- Tegelström H, Gelter HP: Haldane’s rule and sex biased gene flow between two hybridizing flycatcher species (Ficedula albicollis and F. hypoleuca, Aves: Muscicapidae). Evolution. 1990, 44: 2012-2021. 10.2307/2409611.View Article
- Irwin DE, Brelsford A, Toews DPL, MacDonald C, Phinney M: Extensive hybridization in a contact zone between MacGillivray’s warblers Oporornis tolmiei and mourning warblers O. philadelphia detected using molecular and morphological analyses. J Avian Biol. 2009, 40: 539-552. 10.1111/j.1600-048X.2009.04687.x.View Article
- Carling MD, Lovette IJ, Brumfield RT: Historical divergence and gene flow: coalescent analyses of mitochondrial, autosomal and sex-linked loci in Passerina buntings. Evolution. 2010, 64: 1762-1772. 10.1111/j.1558-5646.2010.00954.x.PubMedView Article
- Kingston SE, Jernigan RW, Fagan WF, Braun D, Braun MJ: Genomic variation in cline shape across a hybrid zone. Ecol Evol. 2012, 2: 2737-2748. 10.1002/ece3.375.PubMedPubMed CentralView Article
- Sætre GP, Borge T, Lindell J, Moum T, Primmer CR, Sheldon BC, Haavie J, Johnsen A, Ellegren H: Speciation, introgressive hybridization and nonlinear rate of molecular evolution in flycatchers. Mol Ecol. 2001, 10: 737-749.PubMedView Article
- Sætre G-P, Borge T, Lindroos K, Haavie J, Sheldon BC, Primmer C, Syvänen A-C: Sex chromosome evolution and speciation in Ficedula flycatchers. Proc R Soc Lond Ser B Biol Sci. 2003, 270: 53-59. 10.1098/rspb.2002.2204.View Article
- Mayr E: Systematics and the Origin of Species: From the Viewpoint of a Zoologist. 1942, Cambridge, Massachusetts: Harvard University Press
- Sangster G: The application of species criteria in avian taxonomy and its implications for the debate over species concepts. Biol Rev. 2014, 89: 199-214.PubMedView Article
- Toews DPL, Brelsford A: The biogeography of mitochondrial and nuclear discordance in animals. Mol Ecol. 2012, 21: 3907-3930. 10.1111/j.1365-294X.2012.05664.x.PubMedView Article
- Rheindt FE, Edwards SV: Genetic introgression: an integral but neglected component of speciation in birds. Auk. 2011, 128: 620-632. 10.1525/auk.2011.128.4.620.View Article
- Delaney KS, Zafar S, Wayne RK: Genetic divergence and differentiation within the Western Scrub-Jay (Aphelocoma californica). Auk. 2008, 125: 839-849. 10.1525/auk.2008.07088.View Article
- McCormack JE, Heled J, Delaney KS, Peterson AT, Knowles LL: Calibrating divergence times on species trees versus gene trees: implications for speciation history of Aphelocoma jays. Evolution. 2011, 65: 184-202. 10.1111/j.1558-5646.2010.01097.x.PubMedView Article
- Pitelka FA: Speciation and Ecologic Distribution in American Jays of the Genus Aphelocoma. 1951, Berkeley and Los Angeles, California: University of California Press
- Peterson AT: Phylogeny and rates of molecular evolution in the Aphelocoma jays (Corvidae). Auk. 1992, 109: 133-147. 10.2307/4088274.View Article
- Bardwell E, Benkman CW, Gould WR: Adaptive geographic variation in Western Scrub-Jays. Ecology. 2001, 82: 2617-2627. 10.1890/0012-9658(2001)082[2617:AGVIWS]2.0.CO;2.View Article
- Peterson AT: Adaptive geographical variation in bill shape of scrub jays (Aphelocoma coerulescens). Am Nat. 1993, 142: 508-527. 10.1086/285552.View Article
- Chesser RT, Banks RC, Barker FK, Cicero C, Dunn JL, Kratter AW, Lovette IJ, Rasmussen PC, Remsen JV, Rising JD: Fifty-first supplement to the American Ornithologists’ Union check-list of North American birds. Auk. 2010, 127: 726-744. 10.1525/auk.2010.127.3.726.View Article
- Gowen F, Maley J, Cicero C, Peterson A, Faircloth B, Warr T, McCormack J: Data from: Speciation in Western Scrub-Jays, Haldane’s rule, and genetic clines in secondary contact. Dryad Digital Repository. 2014, http://dx.doi.org/10.5061/dryad.57f48,
- Curry RL, Peterson AT, Langen TA: Western Scrub-Jay (Aphelocoma californica). The Birds of North America. Edited by: Poole A. Ithaca, NY: Cornell Lab of Ornithology
- Yuri T, Jernigan RW, Brumfield RT, Bhagabati NK, Braun MJ: The effect of marker choice on estimated levels of introgression across an avian (Pipridae: Manacus) hybrid zone. Mol Ecol. 2009, 18: 4888-4903. 10.1111/j.1365-294X.2009.04381.x.PubMedView Article
- Sætre G-P, Saether SA: Ecology and genetics of speciation in Ficedula flycatchers. Mol Ecol. 2010, 19: 1091-1106. 10.1111/j.1365-294X.2010.04568.x.PubMedView Article
- Spellman GM, Riddle B, Klicka J: Phylogeography of the Mountain Chickadee (Poecile gambeli): diversification, introgression, and expansion in response to Quaternary climate change. Mol Ecol. 2007, 16: 1055-1068. 10.1111/j.1365-294X.2007.03199.x.PubMedView Article
- Cicero C: Sibling species of titmice in the Parus inornatus complex (Aves: Paridae). Univ Calif Press Publ Zoology. 1996, 128:
- Zellmer AJ, Knowles LL: Disentangling the effects of historic vs. contemporary landscape structure on population genetic divergence. Mol Ecol. 2009, 18: 3593-3602. 10.1111/j.1365-294X.2009.04305.x.PubMedView Article
- Klicka J, Spellman GM, Winker K, Chua V, Smith BT: A phylogeographic and population genetic analysis of a widespread, sedentary North American bird: the Hairy Woodpecker (Picoides villosus). Auk. 2011, 128: 346-362. 10.1525/auk.2011.10264.View Article
- Spellman GM, Klicka J: Phylogeography of the White-breasted Nuthatch (Sitta carolinensis): diversification in North American pine and oak woodlands. Mol Ecol. 2007, 16: 1729-1740. 10.1111/j.1365-294X.2007.03237.x.PubMedView Article
- Betancourt JL, Van Devender TR, Martin PPS: Packrat Middens: The Last 40,000 Years of Biotic Change. 1990, Tucson, Arizona: University of Arizona Press
- Lucas LK, Gompert Z, Ott JR, Nice CC: Geographic and genetic isolation in spring-associated Eurycea salamanders endemic to the Edwards Plateau region of Texas. Conserv Genet. 2009, 10: 1309-1319. 10.1007/s10592-008-9710-2.View Article
- Bowles DE, Arsuffi TL: Karst aquatic ecosystems of the Edwards Plateau region of central Texas, USA: a consideration of their importance, threats to their existence, and efforts for their conservation. Aquat Conserv Mar Freshwat Ecosyst. 1993, 3: 317-329. 10.1002/aqc.3270030406.View Article
- McCormack JE, Smith TB: Niche expansion leads to small-scale adaptive divergence along an elevation gradient in a medium-sized passerine bird. Proc R Soc B. 2008, 275: 2155-2164. 10.1098/rspb.2008.0470.PubMedPubMed CentralView Article
- Peterson AT: New distributional information on the Aphelocoma jays. Bull Br Ornithol Club. 1991, 111: 28-33.
- Miller AH: The avifauna of the Sierra del Carmen of Coahuila, Mexico. Condor. 1955, 57: 154-178. 10.2307/1364864.View Article
- Peterson AT, Martinez-Meyer E, Gonzalez-Salazar C: Reconstructing the Pleistocene geography of the Aphelocoma jays (Corvidae). Divers Distrib. 2004, 10: 237-246. 10.1111/j.1366-9516.2004.00097.x.View Article
- Navarro-Sigüenza AG, Peterson AT: An alternative species taxonomy of the birds of Mexico. Biota Neotrop. 2004, 4: 1-32.View Article
- Burt DB, Peterson AT: Biology of cooperative-breeding Scrub Jays (Aphelocoma coerulescens) of Oaxaca, Mexico. Auk. 1993, 110: 207-214.
- Stenzler LM, Fitzpatrick JW: Isolation of microsatellite loci in the Florida scrub-jay Aphelocoma coerulescens. Mol Ecol Notes. 2002, 2: 547-550. 10.1046/j.1471-8286.2002.00312.x.View Article
- McCormack JE, Peterson AT, Bonaccorso E, Smith TB: Speciation in the highlands of Mexico: genetic and phenotypic divergence in the Mexican jay (Aphelocoma ultramarina). Mol Ecol. 2008, 17: 2505-2521. 10.1111/j.1365-294X.2008.03776.x.PubMedView Article
- van Oosterhout C, Hutchinson WF, Wills DPM, Shipley P: Micro-checker: software for identifying and correcting genotyping errors in microsatellite data. Mol Ecol Notes. 2004, 4: 535-538. 10.1111/j.1471-8286.2004.00684.x.View Article
- Pritchard JK, Stephens M, Donnelly P: Inference of population structure using multilocus genotype data. Genetics. 2000, 155: 945-959.PubMedPubMed Central
- Evanno G, Regnaut S, Goudet J: Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol Ecol. 2005, 14: 2611-2620. 10.1111/j.1365-294X.2005.02553.x.PubMedView Article
- Waples RS, Gaggiotti O: Invited review: what is a population? An empirical evaluation of some genetic methods for identifying the number of gene pools and their degree of connectivity. Mol Ecol. 2006, 15: 1419-1439. 10.1111/j.1365-294X.2006.02890.x.PubMedView Article
- Brown DM, Brenneman RA, Koepfli K-P, Pollinger JP, Milá B, Georgiadis NJ, Louis EE, Grether GF, Jacobs DK, Wayne RK: Extensive population genetic structure in the giraffe. BMC Biol. 2007, 5: 57-10.1186/1741-7007-5-57.PubMedPubMed CentralView Article
- Peterson AT: Evolutionary Relationships of the Aphelocoma Jays. 1990, University of Chicago, Committee on Evolutionary Biology
- Johnson NK: Patterns of avian geography and speciation in the Intermountain Region. Great Basin Nat Mem. 1978, 2: 137-159.
- Drummond AJ, Suchard MA, Xie D, Rambaut A: Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol Biol Evol. 2012, 29: 1969-1973. 10.1093/molbev/mss075.PubMedPubMed CentralView Article
- Nylander J: MrModeltest v2. Evolutionary Biology Centre, Uppsala University, Program distributed by the author. 2004
- Drummond AJ, Ashton B, Buxton S, Cheung M, Cooper A, Duran C, Field M, Heled J, Kearse M, Markowitz S: Geneious v5. 4. 2010, Available from http://www.geneious.com/
- Derryberry EP, Derryberry GE, Maley JM, Brumfield RT: HZAR: hybrid zone analysis using an R software package. Mol Ecol Resour. 2014, 14: 652-66361. 10.1111/1755-0998.12209.PubMedView Article
- Szymura JM, Barton NH: Genetic analysis of a hybrid zone between the fire-bellied toads, Bombina bombina and B. variegate, near Cracow in southern Poland. Evolution. 1986, 40: 1141-1159. 10.2307/2408943.View Article
- Szymura JM, Barton NH: The genetic structure of the hybrid zone between the fire-bellied toads Bombina bombina and B. variegata: comparisons between transects and between loci. Evolution. 1991, 45: 237-261. 10.2307/2409660.View Article
- Barton NH, Gale KS: Genetic Analysis of Hybrid Zones. Hybrid Zones and the Evolutionary Process. Edited by: Harrison RG. 1993, New York, NY: Oxford University Press
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.