Comparative phylogeography and asymmetric hybridization between cryptic bat species

Cryptic speciation and hybridization are two key processes that affect the origin and maintenance of biodiversity and our ability to understand and estimate it. To determine how these two processes interact, we studied allopatric and sympatric colonies of two cryptic bat species (Eptesicus serotinus and Eptesicus isabellinus) with parapatric distribution in the Iberian Peninsula. These species are the main reservoir for the most commonly rabies virus found in bats in Europe: the European bat Lyssavirus type 1 (EBLV‐1). We used mtDNA and nuclear microsatellite markers to confirm the taxonomic status of both species and to show a more pronounced and geographically based genetic structure in E. isabellinus than in its sibling E. serotinus. Using approximate Bayesian computation (ABC), we inferred rapid range expansion in both species after the Last Glacial Maximum until reaching their present distributions. ABC analysis also supported interspecific differences in genetic diversity and structure, pointing to an earlier expansion of E. isabellinus northward. We found no evidence of mitochondrial introgression between species, but nuclear markers identified a male‐mediated ongoing asymmetric hybridization from E. isabellinus to E. serotinus (28% hybrids in E. serotinus and 5% in E. isabellinus) in the contact zone. Although none of the bats studied tested positive for Lyssavirus RNA, the asymmetric hybridization supports the potential for the recently suggested interspecific transmission of EBLV‐1 from E. isabellinus into E. serotinus.     

Late Pleistocene climate change promoted divergence between Picea asperata and P. crassifolia on the Qinghai–Tibet Plateau through recent bottlenecks

Divergence during the early stage of speciation can be driven by a population bottleneck via reduced gene flow and enhanced lineage sorting. In this study, we aimed to examine whether such bottlenecks occurred during the initial speciation of two closely related spruce species Picea asperata and P. crassifolia occurring on the Qinghai–Tibet Plateau (QTP). We analyzed sequences of three chloroplast, two mitochondrial DNA fragments and a further 13 nuclear loci from 216 individuals of the two species. Both species showed a low level of genetic diversity in contrast to other congeners occurring in the QTP and adjacent regions. The estimated population sizes of P. asperata and P. crassifolia are less than the ancestral population size before splitting. These results together with multiple statistical tests (Tajima's D, Fu and Li's D* and F*) suggest that these two species underwent recent bottlenecks. Based on approximate Bayesian computation (ABC), we also determined that the period of the population shrinkage was consistent with the interspecific divergence during the late Pleistocene. The reduced population sizes and the divergent selection may together have triggered the initial divergence under high gene flow between these two species. Our results therefore highlight the importance of climatic oscillations during the late Pleistocene in promoting speciation through changing demographic sizes of the ancestral species on the QTP and in adjacent regions.     

Contrasted levels of genetic diversity in a benthic Mediterranean octocoral: Consequences of different demographic histories?

Understanding the factors explaining the observed patterns of genetic diversity is an important question in evolutionary biology. We provide the first data on the genetic structure of a Mediterranean octocoral, the yellow gorgonian Eunicella cavolini, along with insights into the demographic history of this species. We sampled populations in four areas of the Mediterranean Sea: continental France, Algeria, Turkey, and the Balearic and Corsica islands. Along French coasts, three sites were sampled at two depths (20 and 40 m). We demonstrated a high genetic structure in this species (overall F_ST = 0.13), and most pairwise differentiation tests were significant. We did not detect any difference between depths at the same site. Clustering analyses revealed four differentiated groups corresponding to the main geographical areas. The levels of allelic richness and heterozygosity were significantly different between regions, with highest diversity in Algeria and lowest levels in Turkey. The highest levels of private allelic richness were observed in Algeria followed by Turkey. Such contrasted patterns of genetic diversity were not observed in other Mediterranean octocorals and could be the result of different evolutionary histories. We also provide new empirical evidence of contrasting results between tests and model‐based studies of demographic history. Our results have important consequences for the management of this species.     

Hybridization and introgression across different ploidy levels in the Neotropical orchids Epidendrum fulgens and E. puniceoluteum (Orchidaceae)

The hypothesis of gene flow between species with large differences in chromosome numbers has rarely been tested in the wild, mainly because species of different ploidy are commonly assumed to be reproductively isolated from each other because of instantaneous and strong postzygotic barriers. In this study, a broad‐scale survey of molecular variation was carried out between two orchid species with different ploidy levels: Epidendrum fulgens (2n = 2x = 24 chromosomes) and Epidendrum puniceoluteum (2n = 4x = 52 chromosomes). To test the strength of their reproductive barriers, we investigated the distribution of genetic variation in sympatric and allopatric populations of these two species and conducted crossing experiments. Nuclear and plastid microsatellite loci were used to genotype 463 individuals from eight populations across the geographical range of both species along the Brazilian coastal plain. All six sympatric populations analysed presented hybrid zones, indicating that hybridization between E. fulgens and E. puniceoluteum is a common phenomenon. Bayesian assignment analysis detected the presence of F₁ and F₂ individuals and also signs of introgression, demonstrating a high potential for interspecific gene flow. Introgression occurs preferentially from E. fulgens to E. puniceoluteum. Pure parental individuals of both species display strong genotype–habitat associations, indicating that environment‐dependent selection could be acting in all hybrid zones. This study suggests that hybridization and introgression are evolutionary processes playing a role in the diversification of Epidendrum and indicates the importance of investigations of hybrid zones in understanding reproductive barriers and speciation processes in Neotropical orchid species.     

Evaluation of the utility of COI and ITS markers as tools for population genetic studies of temperate gorgonians

The aim of this study was to analyse the polymorphism of two molecular markers, mitochondrial COI and nuclear ITS 2, to determine the genetic structure of Mediterranean gorgonian populations. Four gorgonians, Eunicella cavolinii, Eunicella singularis, Paramuricea clavata and Corallium rubrum were analysed for COI, while only E. cavolinii was analysed for ITS. These species have been recently affected by mass mortality events probably linked to climate change. In the context of these strong disturbances, reliable data on the genetic structure of affected populations could be critical in order to propose management and conservation strategies for these species. A hierarchical sampling plan encompassing different spatial scales (from tenths to hundreds of kilometres) was designed to test these markers. Both markers failed to make evident any genetic structure of these populations. The evolutionary features and utility of these markers are critically discussed and future research directions are proposed.     

Evidence for hybridization and introgression within a species-rich oak (<Emphasis Type="Italic">Quercus</Emphasis> spp.) community

Background  

Analysis of interspecific gene flow is crucial for the understanding of speciation processes and maintenance of species integrity. Oaks (genus Quercus, Fagaceae) are among the model species for the study of hybridization. Natural co-occurrence of four closely related oak species is a very rare case in the temperate forests of Europe. We used both morphological characters and genetic markers to characterize hybridization in a natural community situated in west-central Romania and which consists of Quercus robur, Q. petraea, Q. pubescen s, and Q. frainetto, respectively.     

Postglacial climate changes and rise of three ecotypes of harbour porpoises,Phocoena phocoena,in western Palearctic waters

Despite no obvious barriers to gene flow in the marine realm, environmental variation and ecological specializations can lead to genetic differentiation in highly mobile predators. Here, we investigated the genetic structure of the harbour porpoise over the entire species distribution range in western Palearctic waters. Combined analyses of 10 microsatellite loci and a 5085 base‐pair portion of the mitochondrial genome revealed the existence of three ecotypes, equally divergent at the mitochondrial genome, distributed in the Black Sea (BS), the European continental shelf waters, and a previously overlooked ecotype in the upwelling zones of Iberia and Mauritania. Historical demographic inferences using approximate Bayesian computation (ABC) suggest that these ecotypes diverged during the last glacial maximum (c. 23–19 kilo‐years ago, kyrbp ). ABC supports the hypothesis that the BS and upwelling ecotypes share a more recent common ancestor (c. 14 kyrbp ) than either does with the European continental shelf ecotype (c. 28 kyrbp ), suggesting they probably descended from the extinct populations that once inhabited the Mediterranean during the glacial and post‐glacial period. We showed that the two Atlantic ecotypes established a narrow admixture zone in the Bay of Biscay during the last millennium, with highly asymmetric gene flow. This study highlights the impacts that climate change may have on the distribution and speciation process in pelagic predators and shows that allopatric divergence can occur in these highly mobile species and be a source of genetic diversity.     

Oxygen consumption in Mediterranean octocorals under different temperatures

Ecosystem resilience to climate anomalies is related to the physiological plasticity of organisms. To characterize the physiological response of some common Mediterranean gorgonians to fluctuations in temperature, four species (Paramuricea clavata, Eunicella singularis, Eunicella cavolinii and Corallium rubrum) were maintained in aquaria, in which the temperature was increased every ten days with increments of 2-3 °C, starting at 14 °C, ending at 25 °C. Oxygen consumption, number of open/closed polyps and percentage of necrotic tissue were monitored. All species showed similar activity patterns with increasing temperature. P. clavata and E. singularis showed the highest respiration rate at 18 °C, E. cavolinii and C. rubrum at 20 °C. Above these temperatures, both oxygen consumption and polyp reactivity decreased in all species. The present data confirm a reduction of the metabolic activity in Mediterranean gorgonians during periods of high temperature. At temperatures above 18 °C, the percentage of open polyps (considered as a parameter to evaluate polyps reactivity) decreased, thus mirroring the trend of oxygen consumption. The average values of Q₁₀ indicated that gorgonians have a definite temperature limit over which the metabolism (oxygen consumption) stop to follow the temperature increase. After three days at 25 °C, metabolic activity in E. cavolinii, C. rubrum and P. clavata further decreased and the first signs of necrosis were observed. At this temperature, activity remained unchanged in E. singularis. This species seems to more resistant to thermal stress. The symbiotic zooxanthellae present in this species are likely to provide an alternative source of energy when polyps reduce their feeding activity.     

Spatial and host related genomic variation in partially sympatric cactophagous moth species

Surveys of patterns of genetic variation in natural sympatric and allopatric populations of recently diverged species are necessary to understand the processes driving intra- and interspecific diversification. The South American moths Cactoblastis cactorum, Cactoblastis doddi and Cactoblastis bucyrus are specialized in the use of cacti as host plants. These species have partially different geographic ranges and differ in patterns of host plant use. However, there are areas that overlap, particularly, in northwestern Argentina, where they are sympatric. Using a combination of genome-wide SNPs and mitochondrial data we assessed intra and interspecific genetic variation and investigated the relative roles of geography and host plants on genetic divergence. We also searched for genetic footprints of hybridization between species. We identified three well delimited species and detected signs of hybridization in the area of sympatry. Our results supported a hypothetical scenario of allopatric speciation in the generalist C. cactorum and genetic interchange during secondary geographic contact with the pair of specialists C. bucyrus and C. doddi that probably speciated sympatrically. In both cases, adaptation to new host plants probably played an important role in speciation. The results also suggested the interplay of geography and host plant use as drivers of divergence and limiting gene flow at intra and interspecific levels.     

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.     

Combining hydrodynamic modelling with genetics: can passive larval drift shape the genetic structure of Baltic Mytilus populations?

While secondary contact between Mytilus edulis and Mytilus trossulus in North America results in mosaic hybrid zone formation, both species form a hybrid swarm in the Baltic. Despite pervasive gene flow, Baltic Mytilus species maintain substantial genetic and phenotypic differentiation. Exploring mechanisms underlying the contrasting genetic composition in Baltic Mytilus species will allow insights into processes such as speciation or adaptation to extremely low salinity. Previous studies in the Baltic indicated that only weak interspecific reproductive barriers exist and discussed the putative role of adaptation to environmental conditions. Using a combination of hydrodynamic modelling and multilocus genotyping, we investigate how oceanographic conditions influence passive larval dispersal and hybrid swarm formation in the Baltic. By combining our analyses with previous knowledge, we show a genetic transition of Baltic Mytilus species along longitude 12°‐13°E, that is a virtual line between Malmö (Sweden) and Stralsund (Germany). Although larval transport only occurs over short distances (10–30 km), limited larval dispersal could not explain the position of this genetic transition zone. Instead, the genetic transition zone is located at the area of maximum salinity change (15–10 psu). Thus, we argue that selection results in weak reproductive barriers and local adaptation. This scenario could maintain genetic and phenotypic differences between Baltic Mytilus species despite pervasive introgressive hybridization.     

Differences in gene flow in a twofold secondary contact zone of pond turtles in southern Italy (Testudines: Emydidae: Emys orbicularis galloitalica,E. o. hellenica,E. trinacris)

Using virtually range‐wide sampling for three pond turtle taxa (Emys orbicularis galloitalica, E. o. hellenica, E. trinacris), we analyse gene flow across their southern Italian contact zone. Based on population genetic analyses of 15 highly polymorphic microsatellite loci and a mitochondrial marker, we show that the general genetic pattern matches well with the current taxon delimitation. Yet, single individuals with conflicting genetic identity suggest translocation of turtles by humans. In addition, we identify in south‐western France and the vicinity of Rome populations being heavily impacted by introduced turtles. Cline analyses reveal that the major genetic break between E. o. galloitalica and E. o. hellenica corresponds well with the currently accepted intergradation zone in southern Italy. However, introgression is largely unidirectional from E. o. galloitalica into E. o. hellenica. In the distribution range of the latter subspecies, genetic footprints of E. o. galloitalica are evident along most of the Italian east coast. Our results corroborate that E. o. galloitalica was introduced long ago in Corsica and Sardinia and naturalized there. Gene flow between E. orbicularis and E. trinacris is negligible, with the Strait of Messina matching well with the narrow cline centre between the two species. This contrasts with other Mediterranean freshwater turtle species with extensive transoceanic gene flow. Compared to the two subspecies of E. orbicularis, the Sicilian E. trinacris shows an unexpectedly strong population structuring, a finding also of some relevance for conservation. The differences between the two taxon pairs E. orbicularis/E. trinacris and E. o. galloitalica/E. o. hellenica support their current taxonomic classification and make them attractive objects for follow‐up studies to elucidate the underlying mechanisms of speciation by comparing their properties.     

Genetic divergence with ongoing gene flow is maintained by the use of different hosts in phytophagous ladybird beetles genus Henosepilachna

Adaptation to different environments can promote population divergence via natural selection even in the presence of gene flow – a phenomenon that typically occurs during ecological speciation. To elucidate how natural selection promotes and maintains population divergence during speciation, we investigated the population genetic structure, degree of gene flow and heterogeneous genomic divergence in three closely related Japanese phytophagous ladybird beetles: Henosepilachna pustulosa, H. niponica and H. yasutomii. These species act as a generalist, a wild thistle (Cirsium spp.) specialist and a blue cohosh (Caulophyllum robustum) specialist, respectively, and their ranges differ accordingly. The two specialist species widely co‐occur but are reproductively isolated solely due to their high specialization to a particular host plant. Genomewide amplified fragment‐length polymorphism (AFLP) markers and mitochondrial cytochrome c oxidase subunit I (COI) gene sequences demonstrated obvious genomewide divergence associated with both geographic distance and ecological divergence. However, a hybridization assessment for both AFLP loci and the mitochondrial sequences revealed a certain degree of unidirectional gene flow between the two sympatric specialist species. Principal coordinates analysis (PCoA) based on all of the variable AFLP loci demonstrated that there are genetic similarities between populations from adjacent localities irrespective of the species (i.e. host range). However, a further comparative genome scan identified a few fractions of loci representing approximately 1% of all loci as different host‐associated outliers. These results suggest that these three species had a complex origin, which could be obscured by current gene flow, and that ecological divergence can be maintained with only a small fraction of the genome is related to different host use even when there is a certain degree of gene flow between sympatric species pairs.     

Phylogenetics and biogeography of the two‐wing flyingfish (Exocoetidae: Exocoetus)

Two‐wing flyingfish (Exocoetus spp.) are widely distributed, epipelagic, mid‐trophic organisms that feed on zooplankton and are preyed upon by numerous predators (e.g., tunas, dolphinfish, tropical seabirds), yet an understanding of their speciation and systematics is lacking. As a model of epipelagic fish speciation and to investigate mechanisms that increase biodiversity, we studied the phylogeny and biogeography of Exocoetus, a highly abundant holoepipelagic fish taxon of the tropical open ocean. Morphological and molecular data were used to evaluate the phylogenetic relationships, species boundaries, and biogeographic patterns of the five putative Exocoetus species. We show that the most widespread species (E. volitans) is sister to all other species, and we find no evidence for cryptic species in this taxon. Sister relationship between E. monocirrhus (Indo‐Pacific) and E. obtusirostris (Atlantic) indicates the Isthmus of Panama and/or Benguela Barrier may have played a role in their divergence via allopatric speciation. The sister species E. peruvianus and E. gibbosus are found in different regions of the Pacific Ocean; however, our molecular results do not show a clear distinction between these species, indicating recent divergence or ongoing gene flow. Overall, our phylogeny reveals that the most spatially restricted species are more recently derived, suggesting that allopatric barriers may drive speciation, but subsequent dispersal and range expansion may affect the distributions of species.     

Genome‐wide single‐nucleotide polymorphism data reveal cryptic species within cryptic freshwater snail species—The case of the Ancylus fluviatilis species complex

DNA barcoding utilizes short standardized DNA sequences to identify species and is increasingly used in biodiversity assessments. The technique has unveiled an unforeseeably high number of morphologically cryptic species. However, if speciation has occurred relatively recently and rapidly, the use of single gene markers, and especially the exclusive use of mitochondrial markers, will presumably fail in delimitating species. Therefore, the true number of biological species might be even higher. One mechanism that can result in rapid speciation is hybridization of different species in combination with polyploidization, that is, allopolyploid speciation. In this study, we analyzed the population genetic structure of the polyploid freshwater snail Ancylus fluviatilis, for which allopolyploidization was postulated as a speciation mechanism. DNA barcoding has already revealed four cryptic species within A. fluviatilis (i.e., A. fluviatilis s. str., Ancylus sp. A–C), but early allozyme data even hint at the presence of additional cryptic lineages in Central Europe. We combined COI sequencing with high‐resolution genome‐wide SNP data (ddRAD data) to analyze the genetic structure of A. fluviatilis populations in a Central German low mountain range (Sauerland). The ddRAD data results indicate the presence of three cryptic species within A. fluviatilis s. str. occurring in sympatry and even syntopy, whereas mitochondrial sequence data only support the existence of one species, with shared haplotypes between species. Our study hence points to the limitations of DNA barcoding when dealing with organismal groups where speciation is assumed to have occurred rapidly, for example, through the process of allopolyploidization. We therefore emphasize that single marker DNA barcoding can underestimate the true species diversity and argue in strong favor of using genome‐wide data for species delimitation in such groups.     

Rapid sexual and genomic isolation in sympatric Drosophila without reproductive character displacement

The rapid evolution of sexual isolation in sympatry has long been associated with reinforcement (i.e., selection to avoid maladaptive hybridization). However, there are many species pairs in sympatry that have evolved rapid sexual isolation without known costs to hybridization. A major unresolved question is what evolutionary processes are involved in driving rapid speciation in such cases. Here, we focus on one such system; the Drosophila athabasca species complex, which is composed of three partially sympatric and interfertile semispecies: WN, EA, and EB. To study speciation in this species complex, we assayed sexual and genomic isolation within and between these semispecies in both sympatric and allopatric populations. First, we found no evidence of reproductive character displacement (RCD) in sympatric zones compared to distant allopatry. Instead, semispecies were virtually completely sexually isolated from each other across their entire ranges. Moreover, using spatial approaches and coalescent demographic simulations, we detected either zero or only weak heterospecific gene flow in sympatry. In contrast, within each semispecies we found only random mating and little population genetic structure, except between highly geographically distant populations. Finally, we determined that speciation in this system is at least an order of magnitude older than previously assumed, with WN diverging first, around 200K years ago, and EA and EB diverging 100K years ago. In total, these results suggest that these semispecies should be given full species status and we adopt new nomenclature: WN—D. athabasca, EA—D. mahican, and EB—D. lenape. While the lack of RCD in sympatry and interfertility do not support reinforcement, we discuss what additional evidence is needed to further decipher the mechanisms that caused rapid speciation in this species complex.     

Maintenance of species boundaries in a Neotropical radiation of Begonia

A major goal of evolutionary biology is to determine the mechanisms generating biodiversity. In Begonia, one of the largest plant genera (1900+ species), it has been postulated that the high number of endemic species is a by‐product of low gene flow among populations, which predisposes the group to speciation. However, this model of divergence requires that reproductive barriers accumulate rapidly among diverging species that overlap in their geographic ranges, otherwise speciation will be opposed by homogenizing gene flow in zones of secondary contact. Here, we test the outcomes of secondary contact in Begonia by genotyping multiple sympatric sites with 12 nuclear and seven plastid loci. We show that three sites of secondary contact between B. heracleifolia and B. nelumbiifolia are highly structured, mostly containing parental genotypes, with few F1 hybrids. A sympatric site between B. heracleifolia and B. sericoneura contains a higher proportion of F1s, but little evidence of introgression. The lack of later‐generation hybrids contrasts with that documented in many other plant taxa, where introgression is extensive. Our results, in conjunction with previous genetic work, show that Begonia demonstrate properties making them exceptionally prone to speciation, at multiple stages along the divergence continuum. Not only are populations weakly connected by gene flow, promoting allopatric speciation, but species often show strong reproductive barriers in secondary contact. Whether similar mechanisms contribute to diversification in other large genera remains to be tested.     

The role of hybridization and introgression in maintaining species integrity and cohesion in naturally isolated inselberg bromeliad populations

• Hybridization is a widespread phenomenon present in numerous lineages across the tree of life. Its evolutionary consequences range from effects on the origin and maintenance, to the loss of biodiversity.
• We studied genetic diversity and intra‐ and interspecific gene flow between two sympatric populations of closely‐related species, Pitcairnia flammea and P. corcovadensis (Bromeliaceae), which are adapted to naturally fragmented Neotropical inselbergs, based on nuclear and plastidial DNA.
• Our main results indicate a strong reproductive isolation barrier, although low levels of interspecific gene flow were observed in both sympatric populations. The low rates of intraspecific gene flow observed for both P. corcovadensis and P. flammea populations corroborate the increasing body of evidence that inselberg bromeliad species are maintained as discrete evolutionary units despite the presence of low genetic connectivity. Nuclear patterns of genetic diversity and gene flow revealed that hybridization and introgression might not cause species extinction via genetic assimilation of the rare P. corcovadensis.
• In the face of reduced intraspecific gene exchange, hybridization and introgression may be important aspects of the Pitcairnia diversification process, with a positive evolutionary impact at the bromeliad community level, and thus contribute to increasing and maintaining genetic diversity in local isolated inselberg populations.

     

Deep mitochondrial introgression and hybridization among ecologically divergent vole species

The completion of speciation is typically difficult to ascertain in rapidly diverging taxa but the amount of hybridization and gene flow in sympatry or parapatry contains important information about the level of reproductive isolation achieved. Here, we examined the progress in speciation between the Mediterranean (Microtus duodecimcostatus) and the Lusitanian pine vole (M. lusitanicus), which are part of the most rapid radiation of species known in mammals. These two Iberian pine voles are classified as separate species because of differences in morphology and ecology, but relatively many ambiguous individuals can be found in sympatric conditions. Our phylogenetic analyses of rangewide data from the mitochondrial cytochrome b gene (mtDNA) demonstrated high levels of diversity and a basal separation in two parapatric lineages. However, mtDNA affiliation was at odds with morphological classification or geographical distribution of the taxa. In contrast, statistical analyses of microsatellites (nucDNA) showed two clear genetic clusters in allopatry and sympatry generally matching morphological classification. This cytonuclear discordance over a large geographic area suggests historical introgression of mtDNA from M. duodecimcostatus to M. lusitanicus. There was statistical evidence for at least two recent hybrids in the sympatry zone but gene flow is apparently low given clear‐cut differences in nucDNA. Our results indicate a relatively advanced speciation process in these Iberian pine voles without fully established reproductive isolation. This situation enables use of combined population genomic and experimental approaches for the separation of patterns and mechanisms in the ongoing explosive diversification of these and other Arvicoline rodents in the future.     

AN EVALUATION OF THE HYBRID SPECIATION HYPOTHESIS FOR XIPHOPHORUS CLEMENCIAE BASED ON WHOLE GENOME SEQUENCES

Once thought rare in animal taxa, hybridization has been increasingly recognized as an important and common force in animal evolution. In the past decade, a number of studies have suggested that hybridization has driven speciation in some animal groups. We investigate the signature of hybridization in the genome of a putative hybrid species, Xiphophorus clemenciae, through whole genome sequencing of this species and its hypothesized progenitors. Based on analysis of this data, we find that X. clemenciae is unlikely to have been derived from admixture between its proposed parental species. However, we find significant evidence for recent gene flow between Xiphophorus species. Although we detect genetic exchange in two pairs of species analyzed, the proportion of genomic regions that can be attributed to hybrid origin is small, suggesting that strong behavioral premating isolation prevents frequent hybridization in Xiphophorus. The direction of gene flow between species is potentially consistent with a role for sexual selection in mediating hybridization.