Process-based species delimitation leads to identification of more biologically relevant species*

Most approaches to species delimitation to date have considered divergence-only models. Although these models are appropriate for allopatric speciation, their failure to incorporate many of the population-level processes that drive speciation, such as gene flow (e.g., in sympatric speciation), places an unnecessary limit on our collective understanding of the processes that produce biodiversity. To consider these processes while inferring species boundaries, we introduce the R-package delimitR and apply it to identify species boundaries in the reticulate taildropper slug (Prophysaon andersoni). Results suggest that secondary contact is an important mechanism driving speciation in this system. By considering process, we both avoid erroneous inferences that can be made when population-level processes such as secondary contact drive speciation but only divergence is considered, and gain insight into the process of speciation in terrestrial slugs. Further, we apply delimitR to three published empirical datasets and find results corroborating previous findings. Finally, we evaluate the performance of delimitR using simulation studies, and find that error rates are near zero when comparing models that include lineage divergence and gene flow for three populations with a modest number of Single Nucleotide Polymorphisms (SNPs; 1500) and moderate divergence times (<100,000 generations). When we apply delimitR to a complex model set (i.e., including divergence, gene flow, and population size changes), error rates are moderate (∼0.15; 10,000 SNPs), and, when present, misclassifications occur among highly similar models.     

Assessing models of speciation under different biogeographic scenarios; an empirical study using multi‐locus and RNA‐seq analyses

Evolutionary biology often seeks to decipher the drivers of speciation, and much debate persists over the relative importance of isolation and gene flow in the formation of new species. Genetic studies of closely related species can assess if gene flow was present during speciation, because signatures of past introgression often persist in the genome. We test hypotheses on which mechanisms of speciation drove diversity among three distinct lineages of desert tortoise in the genus Gopherus. These lineages offer a powerful system to study speciation, because different biogeographic patterns (physical vs. ecological segregation) are observed at opposing ends of their distributions. We use 82 samples collected from 38 sites, representing the entire species' distribution and generate sequence data for mtDNA and four nuclear loci. A multilocus phylogenetic analysis in *BEAST estimates the species tree. RNA‐seq data yield 20,126 synonymous variants from 7665 contigs from two individuals of each of the three lineages. Analyses of these data using the demographic inference package ?a?i serve to test the null hypothesis of no gene flow during divergence. The best‐fit demographic model for the three taxa is concordant with the *BEAST species tree, and the ?a?i analysis does not indicate gene flow among any of the three lineages during their divergence. These analyses suggest that divergence among the lineages occurred in the absence of gene flow and in this scenario the genetic signature of ecological isolation (parapatric model) cannot be differentiated from geographic isolation (allopatric model).     

Tracking the progression of speciation: variable patterns of introgression across the genome provide insights on the species delimitation between progenitor–derivative spruces (Picea mariana × P. rubens)

The genic species concept implies that while most of the genome can be exchanged somewhat freely between species through introgression, some genomic regions remain impermeable to interspecific gene flow. Hence, interspecific differences can be maintained despite ongoing gene exchange within contact zones. This study assessed the heterogeneous patterns of introgression at gene loci across the hybrid zone of an incipient progenitor–derivative species pair, Picea mariana (black spruce) and Picea rubens (red spruce). The spruce taxa likely diverged in geographic isolation during the Pleistocene and came into secondary contact during late Holocene. A total of 300 SNPs distributed across the 12 linkage groups (LG) of black spruce were genotyped for 385 individual trees from 33 populations distributed across the allopatric zone of each species and within the zone of sympatry. An integrative framework combining three population genomic approaches was used to scan the genomes, revealing heterogeneous patterns of introgression. A total of 23 SNPs scattered over 10 LG were considered impermeable to introgression and putatively under diverging selection. These loci revealed the existence of impermeable genomic regions forming the species boundary and are thus indicative of ongoing speciation between these two genetic lineages. Another 238 SNPs reflected selectively neutral diffusion across the porous species barrier. Finally, 39 highly permeable SNPs suggested ancestral polymorphism along with balancing selection. The heterogeneous patterns of introgression across the genome indicated that the speciation process between black spruce and red spruce is young and incomplete, albeit some interspecific differences are maintained, allowing ongoing species divergence even in sympatry. The approach developed in this study can be used to track the progression of ongoing speciation processes.     

Inferring the evolutionary history of divergence despite gene flow in a lizard species,Scincella lateralis (Scincidae), composed of cryptic lineages

Although recent radiations are fruitful for studying the process of speciation, they are difficult to characterize and require the use of multiple loci and analytical methods that account for processes such as gene flow and genetic drift. Using multilocus sequence data, we combine hierarchical cluster analysis, coalescent species tree inference, and isolation‐with‐migration analysis to investigate evolutionary relationships among cryptic lineages of North American ground skinks. We also estimate the extent that gene flow has accompanied or followed diversification, and also attempt to account for and minimize the influence of gene flow when reconstructing relationships. The data best support seven largely parapatric populations that are broadly concordant with mitochondrial (mt)DNA phylogeography throughout most of the species range, although they fail to fully represent extensive mtDNA divergence along the Gulf Coast. Relationships within and among three broad geographical groups are well supported, despite evidence of gene flow among them. Rejection of an allopatric divergence model partially depends on the inclusion of samples from near parapatric boundaries in the analyses, suggesting that allopatric divergence followed by recent migration may best explain migration rate estimates. Accounting for geographical variation in patterns of gene flow can improve estimates of migration–divergence parameters and minimize the influence of contemporary gene flow on phylogenetic inference. © 2012 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, ?? , ??–??.     

Phylogeographic structure,demographic history and morph composition in a colour polymorphic lizard

In polymorphic species, population divergence in morph composition and frequency has the potential to promote speciation. We assessed the relationship between geographic variation in male throat colour polymorphism and phylogeographic structure in the tawny dragon lizard, Ctenophorus decresii. We identified four genetically distinct lineages, corresponding to two polymorphic lineages in the Northern Flinders Ranges and Southern Flinders Ranges/Olary Ranges regions respectively, and a monomorphic lineage in the Mt Lofty Ranges/Kangaroo Island region. The degree of divergence between these three lineages was consistent with isolation to multiple refugia during Pleistocene glacial cycles, whereas a fourth, deeply divergent (at the interspecific level) and monomorphic lineage was restricted to western New South Wales. The same four morphs occurred in both polymorphic lineages, although populations exhibited considerable variation in the frequency of morphs. By contrast, male throat coloration in the monomorphic lineages differed from each other and from the polymorphic lineages. Our results suggest that colour polymorphism has evolved once in the C. decresii species complex, with subsequent loss of polymorphism in the Mt Lofty Ranges/Kangaroo Island lineage. However, an equally parsimonious scenario, that polymorphism arose independently twice within C. decresii, could not be ruled out. We also detected evidence of a narrow contact zone with limited genotypic admixture between the polymorphic Olary Ranges and monomorphic Mt Lofty Ranges regions, yet no individuals of intermediate colour phenotype. Such genetic divergence and evidence for barriers to gene flow between lineages suggest incipient speciation between populations that differ in morph composition.     

Consequences of divergence and introgression for speciation in Andean cloud forest birds

Divergence with gene flow is well documented and reveals the influence of ecological adaptation on speciation. Yet, it remains intuitive that gene exchange inhibits speciation in many scenarios, particularly among ecologically similar populations. The influence of gene flow on the divergence of populations facing similar selection pressures has received less empirical attention than scenarios where differentiation is coupled with local environmental adaptation. I used a paired study design to test the influence of genomic divergence and introgression on plumage differentiation between ecologically similar allopatric replacements of Andean cloud forest birds. Through analyses of short‐read genome‐wide sequences from over 160 individuals in 16 codistributed lineages, I found that plumage divergence is associated with deep genetic divergence, implicating a prominent role of geographic isolation in speciation. By contrast, lineages that lack plumage divergence across the same geographic barrier are more recently isolated or exhibit a signature of secondary genetic introgression, indicating a negative relationship between gene flow and divergence in phenotypic traits important to speciation. My results suggest that the evolutionary outcomes of cycles of isolation and divergence in this important theatre of biotic diversification are sensitive to time spent in the absence of gene flow.     

Climate-driven genetic divergence of limpets with different life histories across a southeast African marine biogeographic disjunction: different processes, same outcome

Genetic divergence among populations of marine broadcast spawners in the absence of past geological barriers presents an intriguing challenge to understanding speciation in the sea. To determine how differences in life history affect genetic divergence and demographic histories across incomplete dispersal barriers, we conducted a comparative phylogeographic study of three intertidal limpets (Siphonaria spp.) represented on either side of a biogeographic disjunction separating tropical and subtropical marine provinces in southeastern Africa. Using a combination of mitochondrial and nuclear sequence data, we identified two distinct evolutionary lineages each in both Siphonaria concinna (a planktonic disperser) and S. nigerrima (a direct developer), and panmixia in a second planktonic disperser, S. capensis. Although phylogeographic breaks were present in two species, how these became established differed depending on their life histories. In the direct developer, lack of gene flow following divergence, and demographic expansion from a small initial size in the species' subtropical population, point to a single colonisation event. In contrast, the evolutionary lineages of the planktonic disperser split into two genetic lineages with much larger initial population sizes and southward gene flow continued at least periodically, indicating that divergence in this species may have been driven by a combination of reduced larval dispersal and divergent selection. These findings help explain why the presence or absence of phylogeographic breaks often appears to be independent of species' dispersal potential.     

Increased differentiation and reduced gene flow in sex chromosomes relative to autosomes between lineages of the brown creeper Certhia americana

The properties of sex chromosomes, including patterns of inheritance, reduced levels of recombination, and hemizygosity in one of the sexes may result in the faster fixation of new mutations via drift and natural selection. Due to these patterns and processes, the two rules of speciation to describe the genetics of postzygotic isolation, Haldane's rule and the large‐X effect, both explicitly include quicker evolution on sex chromosomes relative to autosomes. Because sex‐linked mutations may be the first to become fixed in the speciation process, and appear to be due to stronger genetic drift (in birds), we may identify pronounced genetic differentiation in sex chromosomes in taxa experiencing recent speciation and diverging mainly via genetic drift. Here, we use nine sex‐linked and 21 autosomal genetic markers to investigate differential divergence and introgression between marker types in Certhia americana. We identified increased levels of genetic differentiation and reduced levels of gene flow on sex chromosomes relative to autosomes. This pattern is similar to those observed in other recently‐divergent avian species, providing another case study of the earlier role of sex chromosomes in divergence, relative to autosomes. Additionally, we identify three markers that may be under selection between Certhia americana lineages.     

Plio-Pleistocene climatic oscilations, Holarctic biogeography and speciation in an avian subfamily

Aim In this paper, I discuss the temporal and spatial aspects of historical biogeography and speciation in a widely distributed Holarctic subfamily of birds (Tetraoninae). Location Northern Holarctic. Results Using dated fossils, I calibrated the molecular clock for the mitochondrial control region at 7.23 ± 1.58% nucleotide divergence (maximum likelihood corrected) per million years. The data suggest that grouse (Tetraoninae) originated in the Middle Pliocene, 6.3 Ma. Grouse apparently originated in the northern part of western Nearctic, and Palearctic was colonized independently three times, first by the ancestor of all grouse in the Middle Pliocene, then by the ancestor of forest (Falcipennis, Tetrao and Lyrurus) and prairie (Centrocercus, Dendragapus and Tympanuchus) grouse in the Late Pliocene, and finally by the ancestral Lagopus in the Early Pleistocene. Only once Nearctic was colonized from Palearctic by a common ancestor of forest grouse. Sympatry and range symmetry were positively correlated with molecular divergence. These correlations suggest that peripatric isolation was the predominant mode of speciation throughout grouse history. Main conclusions Speciation events in grouse were driven by climatic oscillations of the Pliocene and Pleistocene. Isolation of small peripheral populations from widely distributed ancestors was the dominant mode of speciation in grouse. Isolations during interglacials both across Beringia, and in southern mountain areas when boreal habitats were restricted to high elevations, suggest an important role for vicariance in grouse speciation.     

Mito‐nuclear discord in six congeneric lineages of Holarctic ducks (genus Anas)

Many species have Holarctic distributions that extend across Europe, Asia and North America. Most genetics research on these species has examined only mitochondrial (mt) DNA, which has revealed wide variance in divergence between Old World (OW) and New World (NW) populations, ranging from shallow, unstructured genealogies to deeply divergent lineages. In this study, we sequenced 20 nuclear introns to test for concordant patterns of OW–NW differentiation between mtDNA and nuclear (nu) DNA for six lineages of Holarctic ducks (genus Anas). Genetic differentiation for both marker types varied widely among these lineages (idiosyncratic population histories), but mtDNA and nuDNA divergence within lineages was not significantly correlated. Moreover, compared with the association between mtDNA and nuDNA divergence observed among different species, OW–NW nuDNA differentiation was generally lower than mtDNA divergence, at least for lineages with deeply divergent mtDNA. Furthermore, coalescent estimates indicated significantly higher rates of gene flow for nuDNA than mtDNA for four of the six lineages. Thus, Holarctic ducks show prominent mito‐nuclear discord between OW and NW populations, and we reject differences in sorting rates as the sole cause of the within‐species discord. Male‐mediated intercontinental gene flow is likely a leading contributor to this discord, although selection could also cause increased mtDNA divergence relative to weak nuDNA differentiation. The population genetics of these ducks contribute to growing evidence that mtDNA can be an unreliable indicator of stage of speciation and that more holistic approaches are needed for species delimitation.     

Population genomics of local adaptation versus speciation in coral reef fishes (Hypoplectrus spp,Serranidae)

Are the population genomic patterns underlying local adaptation and the early stages of speciation similar? Addressing this question requires a system in which (i) local adaptation and the early stages of speciation can be clearly identified and distinguished, (ii) the amount of genetic divergence driven by the two processes is similar, and (iii) comparisons can be repeated both taxonomically (for local adaptation) and geographically (for speciation). Here, we report just such a situation in the hamlets (Hypoplectrus spp), brightly colored reef fishes from the wider Caribbean. Close to 100,000 SNPs genotyped in 126 individuals from three sympatric species sampled in three repeated populations provide genome‐wide levels of divergence that are comparable among allopatric populations (F_st estimate = 0.0042) and sympatric species (F_st estimate = 0.0038). Population genetic, clustering, and phylogenetic analyses reveal very similar patterns for local adaptation and speciation, with a large fraction of the genome undifferentiated (F_st estimate ≈ 0), a very small proportion of F_st outlier loci (0.05–0.07%), and remarkably few repeated outliers (1–3). Nevertheless, different loci appear to be involved in the two processes in Hypoplectrus, with only 7% of the most differentiated SNPs and outliers shared between populations and species comparisons. In particular, a tropomyosin (Tpm4) and a previously identified hox (HoxCa) locus emerge as candidate loci (repeated outliers) for local adaptation and speciation, respectively. We conclude that marine populations may be locally adapted notwithstanding shallow levels of genetic divergence, and that from a population genomic perspective, this process does not appear to differ fundamentally from the early stages of speciation.     

Genetic variability,local selection and demographic history: genomic evidence of evolving towards allopatric speciation in Asian seabass

Genomewide analysis of genetic divergence is critically important in understanding the genetic processes of allopatric speciation. We sequenced RAD tags of 131 Asian seabass individuals of six populations from South‐East Asia and Australia/Papua New Guinea. Using 32 433 SNPs, we examined the genetic diversity and patterns of population differentiation across all the populations. We found significant evidence of genetic heterogeneity between South‐East Asian and Australian/Papua New Guinean populations. The Australian/Papua New Guinean populations showed a rather lower level of genetic diversity. F_ST and principal components analysis revealed striking divergence between South‐East Asian and Australian/Papua New Guinean populations. Interestingly, no evidence of contemporary gene flow was observed. The demographic history was further tested based on the folded joint site frequency spectrum. The scenario of ancient migration with historical population size changes was suggested to be the best fit model to explain the genetic divergence of Asian seabass between South‐East Asia and Australia/Papua New Guinea. This scenario also revealed that Australian/Papua New Guinean populations were founded by ancestors from South‐East Asia during mid‐Pleistocene and were completely isolated from the ancestral population after the last glacial retreat. We also detected footprints of local selection, which might be related to differential ecological adaptation. The ancient gene flow was examined and deemed likely insufficient to counteract the genetic differentiation caused by genetic drift. The observed genomic pattern of divergence conflicted with the ‘genomic islands’ scenario. Altogether, Asian seabass have likely been evolving towards allopatric speciation since the split from the ancestral population during mid‐Pleistocene.     

Repeated and predictable patterns of ecotypic differentiation during a biological invasion: lake–stream divergence in parapatric Swiss stickleback

The relative importance of ecological selection and geographical isolation in promoting and constraining genetic and phenotypic differentiation among populations is not always obvious. Interacting with divergent selection, restricted opportunity for gene flow may in some cases be as much a cause as a consequence of adaptation, with the latter being a hallmark of ecological speciation. Ecological speciation is well studied in parts of the native range of the three‐spined stickleback. Here, we study this process in a recently invaded part of its range. Switzerland was colonized within the past 140 years from at least three different colonization events involving different stickleback lineages. They now occupy diverse habitats, ranging from small streams to the pelagic zone of large lakes. We use replicated systems of parapatric lake and stream populations, some of which trace their origins to different invasive lineages, to ask (i) whether phenotypic divergence occurred among populations inhabiting distinct habitats, (ii) whether trajectories of phenotypic divergence follow predictable parallel patterns and (iii) whether gene flow constrains divergent adaptation or vice versa. We find consistent phenotypic divergence between populations occupying distinct habitats. This involves parallel evolution in several traits with known ecological relevance in independent evolutionary lineages. Adaptive divergence supersedes homogenizing gene flow even at a small spatial scale. We find evidence that adaptive phenotypic divergence places constraints on gene flow over and above that imposed by geographical distance, signalling the early onset of ecological speciation.     

Mitochondrial DNA data imply a stepping‐stone colonization of Beringia by arctic warbler Phylloscopus borealis

Arctic warbler Phylloscopus borealis is one of several high‐latitude Passerines which are widely distributed across one northern continent but restricted to the Beringian part of the other. Most species with such asymmetric intercontinental ranges are monomorphic across Beringia, suggesting either recent colonization of the second continent or considerable gene flow across the Bering Strait. Arctic warbler is the only migratory species in this group that has three different subspecies in Beringia: Ph. b. borealis (Scandinavia to western Beringia, south to Mongolia), Ph. b. xanthodryas (Japan, Sakhalin, Kamchatka, western Beringia), and Ph. b. kennicotti (Alaska). This polymorphism may indicate that Arctic warbler has a unique and complex phylogeographic history that differs significantly from other species with similar ranges. Our analyses of complete mtDNA ND2 sequences of 88 Arctic warblers collected across the species range showed that the clade comprised of birds breeding on Sakhalin Island and Kamchatka Peninsula diverged from the Palearctic/Beringian clade by 3.8% in ND2 sequence. Beringian birds formed a recently derived clade embedded within the Palearctic clade. Nucleotide diversity declined sharply eastward from Palearctic to western Beringia and then to eastern Beringia. Our data provided no support for currently recognized subspecies. They suggested that the barrier at the western edge of Beringia was crossed by Arctic warbler earlier than the Bering Strait resulting in a stepping‐stone colonization of Beringia by this species. Gene flow appears to be restricted across the western border of Beringia but not the Bering Strait.     

The hidden history of the snowshoe hare,Lepus americanus: extensive mitochondrial DNA introgression inferred from multilocus genetic variation

Hybridization drives the evolutionary trajectory of many species or local populations, and assessing the geographic extent and genetic impact of interspecific gene flow may provide invaluable clues to understand population divergence or the adaptive relevance of admixture. In North America, hares (Lepus spp.) are key species for ecosystem dynamics and their evolutionary history may have been affected by hybridization. Here we reconstructed the speciation history of the three most widespread hares in North America – the snowshoe hare (Lepus americanus), the white‐tailed jackrabbit (L. townsendii) and the black‐tailed jackrabbit (L. californicus) – by analysing sequence variation at eight nuclear markers and one mitochondrial DNA (mtDNA) locus (6240 bp; 94 specimens). A multilocus–multispecies coalescent‐based phylogeny suggests that L. americanus diverged ~2.7 Ma and that L. californicus and L. townsendii split more recently (~1.2 Ma). Within L. americanus, a deep history of cryptic divergence (~2.0 Ma) was inferred, which coincides with major speciation events in other North American species. While the isolation‐with‐migration model suggested that nuclear gene flow was generally rare or absent among species or major genetic groups, coalescent simulations of mtDNA divergence revealed historical mtDNA introgression from L. californicus into the Pacific Northwest populations of L. americanus. This finding marks a history of past reticulation between these species, which may have affected other parts of the genome and influence the adaptive potential of hares during climate change.     

Introgression and selection shaped the evolutionary history of sympatric sister‐species of coral reef fishes (genus: Haemulon)

Closely related marine species with large overlapping ranges provide opportunities to study mechanisms of speciation, particularly when there is evidence of gene flow between such lineages. Here, we focus on a case of hybridization between the sympatric sister‐species Haemulon maculicauda and H. flaviguttatum, using Sanger sequencing of mitochondrial and nuclear loci, as well as 2422 single nucleotide polymorphisms (SNPs) obtained via restriction site‐associated DNA sequencing (RADSeq). Mitochondrial markers revealed a shared haplotype for COI and low divergence for CytB and CR between the sister‐species. On the other hand, complete lineage sorting was observed at the nuclear loci and most of the SNPs. Under neutral expectations, the smaller effective population size of mtDNA should lead to fixation of mutations faster than nDNA. Thus, these results suggest that hybridization in the recent past (0.174–0.263 Ma) led to introgression of the mtDNA, with little effect on the nuclear genome. Analyses of the SNP data revealed 28 loci potentially under divergent selection between the two species. The combination of mtDNA introgression and limited nuclear DNA introgression provides a mechanism for the evolution of independent lineages despite recurrent hybridization events. This study adds to the growing body of research that exemplifies how genetic divergence can be maintained in the presence of gene flow between closely related species.     

Demographic Divergence History of Pied Flycatcher and Collared Flycatcher Inferred from Whole-Genome Re-sequencing Data

Profound knowledge of demographic history is a prerequisite for the understanding and inference of processes involved in the evolution of population differentiation and speciation. Together with new coalescent-based methods, the recent availability of genome-wide data enables investigation of differentiation and divergence processes at unprecedented depth. We combined two powerful approaches, full Approximate Bayesian Computation analysis (ABC) and pairwise sequentially Markovian coalescent modeling (PSMC), to reconstruct the demographic history of the split between two avian speciation model species, the pied flycatcher and collared flycatcher. Using whole-genome re-sequencing data from 20 individuals, we investigated 15 demographic models including different levels and patterns of gene flow, and changes in effective population size over time. ABC provided high support for recent (mode 0.3 my, range <0.7 my) species divergence, declines in effective population size of both species since their initial divergence, and unidirectional recent gene flow from pied flycatcher into collared flycatcher. The estimated divergence time and population size changes, supported by PSMC results, suggest that the ancestral species persisted through one of the glacial periods of middle Pleistocene and then split into two large populations that first increased in size before going through severe bottlenecks and expanding into their current ranges. Secondary contact appears to have been established after the last glacial maximum. The severity of the bottlenecks at the last glacial maximum is indicated by the discrepancy between current effective population sizes (20,000–80,000) and census sizes (5–50 million birds) of the two species. The recent divergence time challenges the supposition that avian speciation is a relatively slow process with extended times for intrinsic postzygotic reproductive barriers to evolve. Our study emphasizes the importance of using genome-wide data to unravel tangled demographic histories. Moreover, it constitutes one of the first examples of the inference of divergence history from genome-wide data in non-model species.     

A parapatric propensity for breeding precludes the completion of speciation in common teal (Anas crecca, sensu lato)

Speciation is a process in which genetic drift and selection cause divergence over time. However, there is no rule dictating the time required for speciation, and even low levels of gene flow hinder divergence, so that taxa may be poised at the threshold of speciation for long periods of evolutionary time. We sequenced mitochondrial DNA (mtDNA) and eight nuclear introns (nuDNA) to estimate genomic levels of differentiation and gene flow between the Eurasian common teal (Anas crecca crecca) and the North American green‐winged teal (Anas crecca carolinensis). These ducks come into contact in Beringia (north‐eastern Asia and north‐western North America) and have probably done so, perhaps cyclically, since the Pliocene–Pleistocene transition, ～2.6 Ma, when they apparently began diverging. They have diagnosable differences in male plumage and are 6.9% divergent in the mtDNA control region, with only 1 of 58 crecca and 2 of 86 carolinensis having haplotypes grouping with the other. Two nuclear loci were likewise strongly structured between these teal (Φ_st ≥ 0.35), but six loci were undifferentiated or only weakly structured (Φ_st = 0.0–0.06). Gene flow between crecca and carolinensis was ～1 individual per generation in both directions in mtDNA, but was asymmetrical in nuDNA, with ～1 and ～20 individuals per generation immigrating into crecca and carolinensis, respectively. This study illustrates that species delimitation using a single marker oversimplifies the complexity of the speciation process, and it suggests that even with divergent selection, moderate levels of gene flow may stall the speciation process short of completion.     

Hybridization between two parapatric ranid frog species in the northern Sierra Nevada,California, USA

Contact zones between species provide a unique opportunity to test whether taxa can hybridize or not. Cross‐breeding or hybridization between closely related taxa can promote gene flow (introgression) between species, adaptation, or even speciation. Though hybridization events may be short‐lived and difficult to detect in the field, genetic data can provide information about the level of introgression between closely related taxa. Hybridization can promote introgression between species, which may be an important evolutionary mechanism for either homogenization (reversing initial divergence between species) or reproductive isolation (potentially leading to speciation). Here, we used thousands of genetic markers from nuclear DNA to detect hybridization between two parapatric frog species (Rana boylii and Rana sierrae) in the Sierra Nevada of California. Based on principal components analysis, admixture, and analysis of heterozygosity at species diagnostic SNPs, we detected two F1 hybrid individuals in the Feather River basin, as well as a weak signal of introgression and gene flow between the frog species compared with frog populations from two other adjacent watersheds. This study provides the first documentation of hybridization and introgression between these two species, which are of conservation concern.     

Isolation‐driven divergence: speciation in a widespread North American songbird (Aves: Certhiidae)

Lineage, or true ‘species’, trees may differ from gene trees because of stochastic processes in molecular evolution leading to gene‐tree heterogeneity. Problems with inferring species trees because of excessive incomplete lineage sorting may be exacerbated in lineages with rapid diversification or recent divergences necessitating the use of multiple loci and individuals. Many recent multilocus studies that investigate divergence times identify lineage splitting to be more recent than single‐locus studies, forcing the revision of biogeographic scenarios driving divergence. Here, we use 21 nuclear loci from regional populations to re‐evaluate hypotheses identified in an mtDNA phylogeographic study of the Brown Creeper (Certhia americana), as well as identify processes driving divergence. Nuclear phylogeographic analyses identified hierarchical genetic structure, supporting a basal split at approximately 32°N latitude, splitting northern and southern populations, with mixed patterns of genealogical concordance and discordance between data sets within the major lineages. Coalescent‐based analyses identify isolation, with little to no gene flow, as the primary driver of divergence between lineages. Recent isolation appears to have caused genetic bottlenecks in populations in the Sierra Madre Oriental and coastal mountain ranges of California, which may be targets for conservation concerns.