Timeframe of speciation inferred from secondary contact zones in the European tree frog radiation (Hyla arborea group)

Multiple contact zones and karyotypic evolution in a neotropical frog species complex

Previous studies of DNA sequence and karyotypic data have revealed high genetic diversity in the Physalaemus cuvieri - Physalaemus ephippifer species complex-a group of small leptodactylid frogs in South America. To date, seven major genetic lineages have been recognized in this group, with species delimitation tests supporting four to seven of them as valid species. Among these, only P. ephippifer shows heteromorphic sex chromosomes, but the implications of cytogenetic divergence for the evolution of this group are unknown. We analyzed karyotypic, mitochondrial DNA, and 3RAD genomic data to characterize a putative contact zone between P. ephippifer and P. cuvieri Lineage 1, finding evidence for admixture and karyotypic evolution. We also describe preliminary evidence for admixture between two other members of this species complex-Lineage 1 and Lineage 3 of P. cuvieri. Our study sheds new light on evolutionary relationships in the P. cuvieri - P. ephippifer species complex, suggesting an important role of karyotypic divergence in its evolutionary history and underscoring the importance of hybridization as a mechanism of sex chromosome evolution in amphibians.

Subspecies at crossroads: the evolutionary significance of genomic and phenotypic variation in a wide-ranging Australian lizard (Ctenotus pantherinus)

Many subspecies were described to capture phenotypic variation in wide-ranging taxa, with some later being found to correspond to divergent genetic lineages. We investigate whether currently recognized subspecies correspond to distinctive and coherent evolutionary lineages in the widespread Australian lizard Ctenotus pantherinus based on morphological, mitochondrial and genome-wide nuclear variation. We find weak and inconsistent correspondence between morphological patterns and the presumed subspecies ranges, with character polymorphism within regions and broad morphological overlap across regions. Phylogenetic analyses suggest paraphyly of populations assignable to each subspecies, mitonuclear discordance and little congruence between subspecies ranges and the distribution of inferred clades. Genotypic clustering supports admixture across regions. These results undermine the presumed phenotypic and genotypic coherence and distinctiveness of C. pantherinus subspecies. Based on our findings, we comment on the operational and conceptual shortcomings of morphologically defined subspecies and discuss practical challenges in applying the general notion of subspecies as incompletely separated population lineages. We conclude by highlighting a historical asymmetry that has implications for ecology, evolution and conservation: subspecies proposed in the past are difficult to falsify even in the face of new data that challenge their coherence and distinctiveness, whereas modern researchers appear hesitant to propose new subspecies.

SPART: A versatile and standardized data exchange format for species partition information

A wide range of data types can be used to delimit species and various computer-based tools dedicated to this task are now available. Although these formalized approaches have significantly contributed to increase the objectivity of species delimitation (SD) under different assumptions, they are not routinely used by alpha-taxonomists. One obvious shortcoming is the lack of interoperability among the various independently developed SD programs. Given the frequent incongruences between species partitions inferred by different SD approaches, researchers applying these methods often seek to compare these alternative species partitions to evaluate the robustness of the species boundaries. This procedure is excessively time consuming at present, and the lack of a standard format for species partitions is a major obstacle. Here, we propose a standardized format, SPART, to enable compatibility between different SD tools exporting or importing partitions. This format reports the partitions and describes, for each of them, the assignment of individuals to the "inferred species". The syntax also allows support values to be optionally reported, as well as original trees and the full command lines used in the respective SD analyses. Two variants of this format are proposed, overall using the same terminology but presenting the data either optimized for human readability (matricial SPART) or in a format in which each partition forms a separate block (SPART.XML). ABGD, DELINEATE, GMYC, PTP and TR2 have already been adapted to output SPART files and a new version of LIMES has been developed to import, export, merge and split them.

Using replicate hybrid zones to understand the genomic basis of adaptive divergence

Combining hybrid zone analysis with genomic data is a promising approach to understanding the genomic basis of adaptive divergence. It allows for the identification of genomic regions underlying barriers to gene flow. It also provides insights into spatial patterns of allele frequency change, informing about the interplay between environmental factors, dispersal and selection. However, when only a single hybrid zone is analysed, it is difficult to separate patterns generated by selection from those resulting from chance. Therefore, it is beneficial to look for repeatable patterns across replicate hybrid zones in the same system. We applied this approach to the marine snail Littorina saxatilis, which contains two ecotypes, adapted to wave-exposed rocks vs. high-predation boulder fields. The existence of numerous hybrid zones between ecotypes offered the opportunity to test for the repeatability of genomic architectures and spatial patterns of divergence. We sampled and phenotyped snails from seven replicate hybrid zones on the Swedish west coast and genotyped them for thousands of single nucleotide polymorphisms. Shell shape and size showed parallel clines across all zones. Many genomic regions showing steep clines and/or high differentiation were shared among hybrid zones, consistent with a common evolutionary history and extensive gene flow between zones, and supporting the importance of these regions for divergence. In particular, we found that several large putative inversions contribute to divergence in all locations. Additionally, we found evidence for consistent displacement of clines from the boulder-rock transition. Our results demonstrate patterns of spatial variation that would not be accessible without continuous spatial sampling, a large genomic data set and replicate hybrid zones.

A reassessment of the diversity of green toads (Bufotes) in the circum-Aegean region

The dynamic biogeography of glacial refugia may cause complex patterns of genetic admixture between parapatric taxa, which in turn can mislead their systematics, diversity, and distributions. We investigated this issue for green toads (Bufotes) inhabiting the circum-Aegean region, a biodiversity hotspot of the Eastern Mediterranean. A previous phylogeographic study based on mitochondrial and microsatellite loci identified the hybrid zone between the European (viridis) and Anatolian (sitibundus) lineages of B. viridis all over the Balkan Peninsula, but subsequent range-wide genomic analyses (>1000 SNPs) located this transition in Turkey, a thousand kilometers eastwards. In order to clarify the diversity and taxonomy of the circum-Aegean populations, we reconciled these conflicting findings by integrating previous data with pure sitibundus individuals. Our results confirmed that the viridis/sitibundus hybrid zone extends in Western Anatolia, but that southeastern European populations feature cytonuclear discordances and a high and structured microsatellite diversity. This remarkable situation may stem from a massive geographic displacement of the hybrid zone during the last glacial fluctuations, an underappreciated yet seemingly common feature among the herpetofauna of the region. Our study thus contributes to the rising view that mitochondrial DNA can be a poor predictor of current phylogeographic structure, hence the need for genomic data, especially for narrowly distributed taxa. Finally, the analyses unambiguously support the distinction of a micro-endemic clade of green toads unique to some Cyclades islands, for which we provide a formal taxonomic description.

Dispersal Predicts Hybrid Zone Widths across Animal Diversity: Implications for Species Borders under Incomplete Reproductive Isolation

Hybrid zones occur as range boundaries for many animal taxa. One model for how hybrid zones form and stabilize is the tension zone model, a version of which predicts that hybrid zone widths are determined by a balance between random dispersal into hybrid zones and selection against hybrids. Here, we examine whether random dispersal and proxies for selection against hybrids (genetic distances between hybridizing pairs) can explain variation in hybrid zone widths across 131 hybridizing pairs of animals. We show that these factors alone can explain ∼40% of the variation in zone width among animal hybrid zones, with dispersal explaining far more of the variation than genetic distances. Patterns within clades were idiosyncratic. Genetic distances predicted hybrid zone widths particularly well for reptiles, while this relationship was opposite tension zone predictions in birds. Last, the data suggest that dispersal and molecular divergence set lower bounds on hybrid zone widths in animals, indicating that there are geographic restrictions on hybrid zone formation. Overall, our analyses reinforce the fundamental importance of dispersal in hybrid zone formation and more generally in the ecology of range boundaries.

The effect of phylogeographic history on species boundaries: a comparative framework in Hyla tree frogs

Because it is indicative of reproductive isolation, the amount of genetic introgression across secondary contact zones is increasingly considered in species delimitation. However, patterns of admixture at range margins can be skewed by the regional dynamics of hybrid zones. In this context, we posit an important role for phylogeographic history: hybrid zones located within glacial refugia (putatively formed during the Late-Pleistocene) should be better defined than those located in post-glacial or introduced ranges (putatively formed during the Holocene and the Anthropocene). We test this hypothesis in a speciation continuum of tree frogs from the Western Palearctic (Hyla), featuring ten identified contacts between species spanning Plio-Pleistocene to Miocene divergences. We review the rich phylogeographic literature of this group and examine the overlooked transition between H. arborea and H. molleri in Western France using a multilocus dataset. Our comparative analysis supports a trend that contacts zones resulting from post-glacial expansions and human translocations feature more extensive introgression than those established within refugial areas. Integrating the biogeographic history of incipient species, i.e. their age since first contact together with their genetic divergence, thus appears timely to draw sound evolutionary and taxonomic inferences from patterns of introgression across hybrid zones.

Are glacial refugia hotspots of speciation and cytonuclear discordances? Answers from the genomic phylogeography of Spanish common frogs

Subdivided Pleistocene glacial refugia, best known as "refugia within refugia", provided opportunities for diverging populations to evolve into incipient species and/or to hybridize and merge following range shifts tracking the climatic fluctuations, potentially promoting extensive cytonuclear discordances and "ghost" mtDNA lineages. Here, we tested which of these opposing evolutionary outcomes prevails in northern Iberian areas hosting multiple historical refugia of common frogs (Rana cf. temporaria), based on a genomic phylogeography approach (mtDNA barcoding and RAD-sequencing). We found evidence for both incipient speciation events and massive cytonuclear discordances. On the one hand, populations from northwestern Spain (Galicia and Asturias, assigned to the regional endemic R. parvipalmata), are deeply-diverged at mitochondrial and nuclear genomes (~4 My of independent evolution), and barely admix with northeastern populations (assigned to R. temporaria sensu stricto) across a narrow hybrid zone (~25 km) located in the Cantabrian Mountains, suggesting that they represent distinct species. On the other hand, the most divergent mtDNA clade, widespread in Cantabria and the Basque country, shares its nuclear genome with other R. temporaria s. s. lineages. Patterns of population expansions and isolation-by-distance among these populations are consistent with past mitochondrial capture and/or drift in generating and maintaining this ghost mitochondrial lineage. This remarkable case study emphasizes the complex evolutionary history that shaped the present genetic diversity of refugial populations, and stresses the need to revisit their phylogeography by genomic approaches, in order to make informed taxonomic inferences.

Discordant patterns of introgression across a narrow hybrid zone between two cryptic lineages of an Iberian endemic newt

The study of natural hybrid zones can illuminate aspects of lineage divergence and speciation in morphologically cryptic taxa. We studied a hybrid zone between two highly divergent but morphologically similar lineages (south-western and south-eastern) of the Iberian endemic Bosca's newt (Lissotriton boscai) in SW Iberia with a multilocus dataset (microsatellites, nuclear and mitochondrial genes). STRUCTURE and NEWHYBRIDS analyses retrieved few admixed individuals, which classified as backcrosses involving parental individuals of the south-western lineage. Our results show asymmetric introgression of mtDNA beyond the contact from this lineage into the south-eastern lineage. Analysis of nongeographic introgression patterns revealed asymmetries in the direction of introgression, but except for mtDNA, we did not find evidence for nonconcordant introgression patterns across nuclear loci. Analysis of a 150-km transect across the hybrid zone showed broadly coincident cline widths (ca. 3.2-27.9 km), and concordant cline centres across all markers, except for mtDNA that is displaced ca. 60 km northward. Results from ecological niche modelling show that the hybrid zone is in a climatically homogenous area with suitable habitat for the species, suggesting that contact between the two lineages is unlikely to occur further south as their distributions are currently separated by an extensive area of unfavourable habitat. Taken together, our findings suggest the genetic structure of this hybrid zone results from the interplay of historical (biogeographic) and population-level processes. The narrowness and coincidence of genetic clines can be explained by weak selection against hybrids and reflect a degree of reproductive isolation that is consistent with cryptic speciation.

Phylogeography of a cryptic speciation continuum in Eurasian spadefoot toads (Pelobates)

Cryptic phylogeographic diversifications provide unique models to examine the role of phylogenetic divergence on the evolution of reproductive isolation, without extrinsic factors such as ecological and behavioural differentiation. Yet, to date very few comparative studies have been attempted within such radiations. Here, we characterize a new speciation continuum in a group of widespread Eurasian amphibians, the Pelobates spadefoot toads, by conducting multilocus (restriction site associated DNA sequencing and mitochondrial DNA) phylogenetic, phylogeographic and hybrid zone analyses. Within the P. syriacus complex, we discovered species-level cryptic divergences (>5 million years ago [My]) between populations distributed in the Near-East (hereafter P. syriacus sensu stricto [s.s.]) and southeastern Europe (hereafter P. balcanicus), each featuring deep intraspecific lineages. Altogether, we could scale hybridizability to divergence time along six different stages, spanning from sympatry without gene flow (P. fuscus and P. balcanicus, >10 My), parapatry with highly restricted hybridization (P. balcanicus and P. syriacus s.s., >5 My), narrow hybrid zones (~15 km) consistent with partial reproductive isolation (P. fuscus and P. vespertinus, ~3 My), to extensive admixture between Pleistocene and refugial lineages (≤2 My). This full spectrum empirically supports a gradual build up of reproductive barriers through time, reversible up until a threshold that we estimate at ~3 My. Hence, cryptic phylogeographic lineages may fade away or become reproductively isolated species simply depending on the time they persist in allopatry, and without definite ecomorphological divergence.

A tree of tree frogs around the Black Sea

Speciation, the process by which one species evolves into two or more, is a major focus of ongoing debate, particularly regarding the geographic context in which it occurs. Geographic models of speciation tend to fall into discrete categories, typically referred to as allopatric, parapatric and sympatric speciation, according to whether two groups evolve reproductive isolation while geographically isolated, differentiated but connected by gene flow, or completely co-occurring. Yet molecular studies indicate that full development of reproductive isolation can take very long compared with the timescale at which climatic oscillations occur, such that the geographic context of differentiating forms might change often during the long process to full species. Studies of genetic relationships across the ranges of organisms with low-dispersal distances have the potential to reveal these complex histories. In a particularly elegant example in this issue, Dufresnes et al. () use genetic variation and ecological niche modelling to show that a ring of populations of the eastern tree frog (Hyla orientalis) surrounding the Black Sea had a complex history of geographic differentiation. Alternating phases of geographic fragmentation and phases of gene flow between neighbouring populations have produced a pattern of gradual genetic change connecting the western, southern and eastern sides of the ring, with the northwestern and northeastern forms being most differentiated. In the north, a population in Crimea appears to have been produced through mixture of the two extreme forms. The overall genetic relationships are reminiscent of those found in ring species, which have been used as prime demonstrations of the process of speciation. The difference, however, is that the terminal forms appear to have mixed rather than be reproductively isolated, although more research is needed to infer whether there might be some reproductive isolation on the northern side of the ring.

A Maladaptive Combination of Traits Contributes to the Maintenance of a Drosophila Hybrid Zone

Drosophila teissieri and D. yakuba diverged approximately 3 mya and are thought to share a large, ancestral, African range [1-3]. These species now co-occur in parts of continental Africa and in west Africa on the island of Bioko [1, 4]. While D. yakuba is a human commensal, D. teissieri seems to be associated with Parinari fruits, restricting its range to forests [4-6]. Genome data indicate introgression, despite no evidence of contemporary hybridization. Here we report the discovery of D. yakuba-D. teissieri hybrids at the interface of secondary forests and disturbed, open habitats on Bioko. We demonstrate that hybrids are the F₁ progeny of D. yakuba females and D. teissieri males. At high temperatures like those found on Bioko, D. teissieri females are generally less receptive to mating, and in combination with temperature effects on egg lay and egg-to-adult viability, this decreases the potential for gene flow between female D. teissieri and male D. yakuba relative to the reciprocal cross. Field and laboratory experiments demonstrate that F₁ hybrids have a maladaptive combination of D. yakuba behavior and D. teissieri physiology, generating additional barriers to gene flow. Nevertheless, analysis of introgressed and non-introgressed regions of the genome indicate that, while rare, gene flow is relatively recent. Our observations identify precise intrinsic and extrinsic factors that, along with hybrid male sterility, limit gene flow and maintain these species. These data contribute to a growing body of literature that suggests the Gulf of Guinea may be a hotspot for hybridization.

Phylogeography of Aegean green toads (Bufo viridis subgroup): continental hybrid swarm vs. insular diversification with discovery of a new island endemic

Background

Debated aspects in speciation research concern the amount of gene flow between incipient species under secondary contact and the modes by which post-zygotic isolation accumulates. Secondary contact zones of allopatric lineages, involving varying levels of divergence, provide natural settings for comparative studies, for which the Aegean (Eastern Mediterranean) geography offers unique scenarios. In Palearctic green toads (Bufo viridis subgroup or Bufotes), Plio-Pleistocene (~ 2.6 Mya) diverged species show a sharp transition without contemporary gene flow, while younger lineages, diverged in the Lower-Pleistocene (~ 1.9 Mya), admix over tens of kilometers. Here, we conducted a fine-scale multilocus phylogeographic analysis of continental and insular green toads from the Aegean, where a third pair of taxa, involving Mid-Pleistocene diverged (~ 1.5 Mya) mitochondrial lineages, earlier tentatively named viridis and variabilis, (co-)occurs.

Results

We discovered a new lineage, endemic to Naxos (Central Cyclades), while coastal islands and Crete feature weak genetic differentiation from the continent. In continental Greece, both lineages, viridis and variabilis, form a hybrid swarm, involving massive mitochondrial and nuclear admixture over hundreds of kilometers, without obvious selection against hybrids.

Conclusions

The genetic signatures of insular Aegean toads appear governed by bathymetry and Quaternary sea level changes, resulting in long-term isolation (Central Cyclades: Naxos) and recent land-bridges (coastal islands). Conversely, Crete has been isolated since the end of the Messinian salinity crisis (5.3 My) and Cretan populations thus likely result from human-mediated colonization, at least since Antiquity, from Peloponnese and Anatolia. Comparisons of green toad hybrid zones support the idea that post-zygotic hybrid incompatibilities accumulate gradually over the genome. In this radiation, only one million years of divergence separate a scenario of complete reproductive isolation, from a secondary contact resulting in near panmixia.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1179-0) contains supplementary material, which is available to authorized users.

The unusual case of the widely distributed fiddler crab Minuca rapax (Smith, 1870) from the western Atlantic: an exemplary polytypic species

A classic dilemma in taxonomy is distinguishing intraspecific from interspecific variation. In order to better comprehend the process of divergence and speciation, we examine morphological, genetic, developmental and behavioural variation among related fiddler crab populations from eastern North America, the Caribbean and South America. We chose geographically remote populations that appear related to Minuca rapax (Smith, 1870). First, using females from across the range of the species, we use geometric morphometric techniques to identify regional differences in carapace shape. Second, in the northern portion of the range, the Caribbean into the Gulf of Mexico, we report variation in the relationship between corporal size and cheliped length in males. Third, we examine the major components of the courtship waves produced by males from several locations in the western Gulf of Mexico. Fourth, we compare the structure of the gastric mill between different populations in the Gulf of Mexico, the Caribbean and the Atlantic Ocean. And, fifth, we use mitochondrial 16S rDNA and cytochrome oxidase subunit I as genetic markers to define the phylogeographic relationship among specimens from more than 20 populations. From these studies, we find discrete, distinct populations across the original range of the species. In particular, populations in the northern Gulf of Mexico appear to represent a lineage that has resulted from limited gene flow and sustained selection pressures. On the basis of the observed degree of divergence, it is apparent that some separated populations in M. rapax should be recognised as evolutionary significant units. The geographic range of these populations is consistent with the historical range for Minuca virens (Salmon & Atsaides, 1968), a putative species that otherwise cannot be consistently distinguished from M. rapax based on discrete external morphological characters. This study provides evidence for M. virens as an emergent but possibly not completely isolated subclade of the M. rapax species complex.

Multiple uprising invasions of Pelophylax water frogs, potentially inducing a new hybridogenetic complex

The genetic era has revolutionized our perception of biological invasions. Yet, it is usually too late to understand their genesis for efficient management. Here, we take the rare opportunity to reconstruct the scenario of an uprising invasion of the famous water frogs (Pelophylax) in southern France, through a fine-scale genetic survey. We identified three different taxa over less than 200 km²: the autochthonous P. perezi, along with the alien P. ridibundus and P. kurtmuelleri, which have suddenly become invasive. As a consequence, the latter hybridizes and may now form a novel hybridogenetic complex with P. perezi, which could actively promote its replacement. This exceptional situation makes a textbook application of genetics to early-detect, monitor and understand the onset of biological invasions before they pose a continental-wide threat. It further emphasizes the alarming rate of amphibian translocations, both at global and local scales, as well as the outstanding invasive potential of Pelophylax aliens.

Integrative inference of population history in the Ibero-Maghrebian endemic Pleurodeles waltl (Salamandridae)

Inference of population histories from the molecular signatures of past demographic processes is challenging, but recent methodological advances in species distribution models and their integration in time-calibrated phylogeographic studies allow detailed reconstruction of complex biogeographic scenarios. We apply an integrative approach to infer the evolutionary history of the Iberian ribbed newt (Pleurodeles waltl), an Ibero-Maghrebian endemic with populations north and south of the Strait of Gibraltar. We analyzed an extensive multilocus dataset (mitochondrial and nuclear DNA sequences and ten polymorphic microsatellite loci) and found a deep east-west phylogeographic break in Iberian populations dating back to the Plio-Pleistocene. This break is inferred to result from vicariance associated with the formation of the Guadalquivir river basin. In contrast with previous studies, North African populations showed exclusive mtDNA haplotypes, and formed a monophyletic clade within the Eastern Iberian lineage in the mtDNA genealogy. On the other hand, microsatellites failed to recover Moroccan populations as a differentiated genetic cluster. This is interpreted to result from post-divergence gene flow based on the results of IMA2 and Migrate analyses. Thus, Moroccan populations would have originated after overseas dispersal from the Iberian Peninsula in the Pleistocene, with subsequent gene flow in more recent times, implying at least two trans-marine dispersal events. We modeled the distribution of the species and of each lineage, and projected these models back in time to infer climatically favourable areas during the mid-Holocene, the last glacial maximum (LGM) and the last interglacial (LIG), to reconstruct more recent population dynamics. We found minor differences in climatic favourability across lineages, suggesting intraspecific niche conservatism. Genetic diversity was significantly correlated with the intersection of environmental favourability in the LIG and LGM, indicating that populations of P. waltl are genetically more diverse in regions that have remained environmentally favourable through the last glacial cycle, particularly southern Iberia and northern Morocco. This study provides novel insights into the relative roles of geology and climate on the biogeography of a biodiversity hotspot.

Limited genomic consequences of hybridization between two African clawed frogs, Xenopus gilli and X. laevis (Anura: Pipidae)

The Cape platanna, Xenopus gilli, an endangered frog, hybridizes with the African clawed frog, X. laevis, in South Africa. Estimates of the extent of gene flow between these species range from pervasive to rare. Efforts have been made in the last 30 years to minimize hybridization between these two species in the west population of X. gilli, but not the east populations. To further explore the impact of hybridization and the efforts to minimize it, we examined molecular variation in one mitochondrial and 13 nuclear genes in genetic samples collected recently (2013) and also over two decades ago (1994). Despite the presence of F₁ hybrids, none of the genomic regions we surveyed had evidence of gene flow between these species, indicating a lack of extensive introgression. Additionally we found no significant effect of sampling time on genetic diversity of populations of each species. Thus, we speculate that F₁ hybrids have low fitness and are not backcrossing with the parental species to an appreciable degree. Within X. gilli, evidence for gene flow was recovered between eastern and western populations, a finding that has implications for conservation management of this species and its threatened habitat.

Concordant morphological and molecular clines in a contact zone of the Common and Spined toad (Bufo bufo and B. spinosus) in the northwest of France

BACKGROUND

Hybrid zones are regions where individuals of two species meet and produce hybrid progeny, and are often regarded as natural laboratories to understand the process of species formation. Two microevolutionary processes can take place in hybrid zones, with opposing effects on population differentiation. Hybridization tends to produce genetic homogenization, reducing species differences, whereas the presence of mechanisms of reproductive isolation result in barriers to gene flow, maintaining or increasing differences between taxa.

RESULTS

Here we study a contact zone between two hybridizing toad species, Bufo bufo and B. spinosus, through a combination of molecular (12 polymorphic microsatellites, four nuclear and two mitochondrial SNP markers) and morphological data in a transect in the northwest of France. The results show largely concordant clines across markers, defining a narrow hybrid zone of ca. 30 km wide. Most hybrids in the centre of the contact zone are classified as F₂ or backcrossed individuals, with no individuals assigned to the F₁ hybrid class.

CONCLUSIONS

We discuss the implications of these results for our understanding of the evolutionary history of these species. We anticipate that the toad contact zone here described will become an important asset in the study of hybrid zone dynamics and evolutionary biology because of its easy access and the abundance of the species involved.

Phylogeography reveals an ancient cryptic radiation in East-Asian tree frogs (Hyla japonica group) and complex relationships between continental and island lineages

Background

In contrast to the Western Palearctic and Nearctic biogeographic regions, the phylogeography of Eastern-Palearctic terrestrial vertebrates has received relatively little attention. In East Asia, tectonic events, along with Pleistocene climatic conditions, likely affected species distribution and diversity, especially through their impact on sea levels and the consequent opening and closing of land-bridges between Eurasia and the Japanese Archipelago. To better understand these effects, we sequenced mitochondrial and nuclear markers to determine phylogeographic patterns in East-Asian tree frogs, with a particular focus on the widespread H. japonica.

Results

We document several cryptic lineages within the currently recognized H. japonica populations, including two main clades of Late Miocene divergence (~5 Mya). One occurs on the northeastern Japanese Archipelago (Honshu and Hokkaido) and the Russian Far-East islands (Kunashir and Sakhalin), and the second one inhabits the remaining range, comprising southwestern Japan, the Korean Peninsula, Transiberian China, Russia and Mongolia. Each clade further features strong allopatric Plio-Pleistocene subdivisions (~2-3 Mya), especially among continental and southwestern Japanese tree frog populations. Combined with paleo-climate-based distribution models, the molecular data allowed the identification of Pleistocene glacial refugia and continental routes of postglacial recolonization. Phylogenetic reconstructions further supported genetic homogeneity between the Korean H. suweonensis and Chinese H. immaculata, suggesting the former to be a relic population of the latter that arose when the Yellow Sea formed, at the end of the last glaciation.

Conclusions

Patterns of divergence and diversity were likely triggered by Miocene tectonic activities and Quaternary climatic fluctuations (including glaciations), causing the formation and disappearance of land-bridges between the Japanese islands and the continent. Overall, this resulted in a ring-like diversification of H. japonica around the Sea of Japan. Our findings urge for important taxonomic revisions in East-Asian tree frogs. First, they support the synonymy of H. suweonensis (Kuramoto, 1980) and H. immaculata (Boettger, 1888). Second, the nominal H. japonica (Günther, 1859) represents at least two species: an eastern (new taxon A) on the northern Japanese and Russian Far East islands, and a southwestern species (n. t. B) on southern Japanese islands and possibly also forming continental populations. Third, these continental tree frogs may also represent an additional entity, previously described as H. stepheni Boulenger, 1888 (senior synonym of H. ussuriensis Nikolskii, 1918). A complete revision of this group requires further taxonomic and nomenclatural analyses, especially since it remains unclear to which taxon the species-epitheton japonica corresponds to.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0814-x) contains supplementary material, which is available to authorized users.

Microsatellite markers and cytoplasmic sequences reveal contrasting pattern of spatial genetic structure in the red algae species complex Mazzaella laminarioides

Mazzaella laminarioides is a common haploid-diploid red alga that forms dense beds. This alga has a wide distributional range, covering 3,500 km of the Chilean coast, but is restricted to high rocky intertidal zones. Recently, the existence of three highly divergent genetic lineages was demonstrated for this taxon, and two cytoplasmic markers were used to determine that these lineages are distributed in strict parapatry. Here, using 454 next-generation sequencing, we developed polymorphic microsatellite loci that cross amplify in all three cytoplasmic lineages. Six sites (i.e., two sites within each lineage) were analyzed using nine microsatellite loci. Our work shows that, although substantial cytoplasmic differentiation occurs within M. laminarioides, the microsatellite loci did not retrieve three nuclear genetic clusters as expected. Indeed, while the northernmost and southernmost cytoplasmic lineages form two strongly divergent nuclear groups characterized by diagnostic alleles, the third cytoplasmic lineage did not form a third nuclear independent group. It is possible that inter-lineage gene exchange has occurred, particularly at sites along the contact zone between the different cytoplasmic lineages. This nuclear-cytoplasmic incongruence in M. laminarioides could be explained by incomplete lineage sorting of the nuclear genes or asymmetric introgressive hybridization between the lineages. Finally, highly significant heterozygote deficiencies (suggesting occurrence of intergametophytic selfing) were observed in the three small northernmost sites while the large southernmost sites generally approached panmixia.

Evolutionary melting pots: a biodiversity hotspot shaped by ring diversifications around the Black Sea in the Eastern tree frog (Hyla orientalis)

Hotspots of intraspecific genetic diversity, which are of primary importance for the conservation of species, have been associated with glacial refugia, that is areas where species survived the Quaternary climatic oscillations. However, the proximate mechanisms generating these hotspots remain an open issue. Hotspots may reflect the long-term persistence of large refugial populations; alternatively, they may result from allopatric differentiation between small and isolated populations, that later admixed. Here, we test these two scenarios in a widely distributed species of tree frog, Hyla orientalis, which inhabits Asia Minor and southeastern Europe. We apply a fine-scale phylogeographic survey, combining fast-evolving mitochondrial and nuclear markers, with a dense sampling throughout the range, as well as ecological niche modelling, to understand what shaped the genetic variation of this species. We documented an important diversity centre around the Black Sea, composed of multiple allopatric and/or parapatric diversifications, likely driven by a combination of Pleistocene climatic fluctuations and complex regional topography. Remarkably, this diversification forms a ring around the Black Sea, from the Caucasus through Anatolia and eastern Europe, with terminal forms coming into contact and partially admixing in Crimea. Our results support the view that glacial refugia generate rather than host genetic diversity and can also function as evolutionary melting pots of biodiversity. Moreover, we report a new case of ring diversification, triggered by a large, yet cohesive dispersal barrier, a very rare situation in nature. Finally, we emphasize the Black Sea region as an important centre of intraspecific diversity in the Palearctic with implications for conservation.

Empirical evidence for large X-effects in animals with undifferentiated sex chromosomes

Reproductive isolation is crucial for the process of speciation to progress. Sex chromosomes have been assigned a key role in driving reproductive isolation but empirical evidence from natural population processes has been restricted to organisms with degenerated sex chromosomes such as mammals and birds. Here we report restricted introgression at sex-linked compared to autosomal markers in a hybrid zone between two incipient species of European tree frog, Hyla arborea and H. orientalis, whose homologous X and Y sex chromosomes are undifferentiated. This large X-effect cannot result from the dominance or faster-X aspects of Haldane's rule, which are specific to degenerated sex chromosomes, but rather supports a role for faster-heterogametic-sex or faster-male evolutionary processes. Our data suggest a prominent contribution of undifferentiated sex chromosomes to speciation.