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

Integrative taxonomy, phylogenetics and historical biogeography of subgenus Aeschyntelus Stål, 1872 (Hemiptera: Heteroptera: Rhopalidae)

The subgenus Aeschyntelus includes six species that show variations in body color and shape, thus making it difficult to identify them based on morphological identification alone. To date, no genetic study has evaluated species within this genus. Herein, we collected 171 individuals from 90 localities of Rhopalus and employed an integrative taxonomic approach that incorporated morphological data, mitochondrial genomic data (COI, whole mitochondrial data) and nuclear genomic data (18S + 28S rRNAs, nuclear genome-wide SNPs) to delineate species boundaries. Our analyses confirmed the status of nine described species of Rhopalus and proposed the recognition of one new species known as Rhopalus qinlinganus sp. nov., which is classified within the subgenus Aeschyntelus. Discrepancies arising from nuclear and mitochondrial data suggest the presence of mito-nuclear discordance. Specifically, mitochondrial data indicated admixture within Clade A, comprising R. kerzhneri and R. latus, whereas genome-wide SNPs unambiguously identified two separate species, aligning with morphological classification. Conversely, mitochondrial data clearly distinguished Clade B- consisting of R. sapporensis into two lineages, whereas genome-wide SNPs unequivocally identified a single species. Our study also provides insights into the evolutionary history of Aeschyntelus, thus indicating that it likely originated in East Asia during the middle Miocene. The development of Aeschyntelus biodiversity in the southwestern mountains of China occurred via an uplift-driven diversification process. Our findings highlight the necessity of integrating both morphological and multiple molecular datasets for precise species identification, particularly when delineating closely related species. Additionally, it reveals the important role of mountain orogenesis on speciation within the southwestern mountains of China.

Phylogeography and paleoclimatic range dynamics explain variable outcomes to contact across a species' range

Replicability of divergence after contact is a poorly characterized process, particularly in the contexts of phylogeography and postglacial range dynamics within species. Using contact zones located at the leading-, mid- and rear-edges of a species' range, we examined variation in outcomes to contact between divergent lineages of Campanula americana. We investigated whether contact zones vary in quantity and directionality of gene flow, how phylogeographic structure differs between contact zones, and how historic range dynamics may affect outcomes to contact. We found that all contact zones formed at similar times via primary contact yet detected significant admixture in only the rear-edge (RE) contact zone. In the northern leading-edge contact zone and the mid-range Virginia contact zone, gene flow was minimal and asymmetric. In the southern RE contact zone, gene flow was strong and symmetric. Asymmetric admixture in the leading-edge and Virginia contact zones matches the directionality of a known cosmopolitan cytonuclear incompatibility between lineages of C. americana. Our results emphasize the dependence of speciation processes on phylogeographic structure, evolutionary history and range dynamics.

"Isolation by Gentes with Asymmetric Migration" shapes the genetic structure of the common cuckoo in China

Amid coevolutionary arms races between brood parasitic birds and their diverse host species, the formation of host-specific races, or gentes, has drawn significant research focus. Nevertheless, numerous questions about gentes evolutionary patterns persist. Here, we investigated the potential for gentes evolution across multiple common cuckoo (Cuculus canorus) populations parasitizing diverse host species in China. Using maternal (mitochondrial and W-linked DNA) and biparental (autosomal and Z-linked DNA) markers, we found consistent clustering of cuckoo gentes (rather than geographical populations) into distinct clades in matrilineal gene trees, indicating robust differentiation. In contrast, biparental markers indicated intermixing of all gentes, suggesting asymmetric gene flow regardless of geography. Unlike the mitonuclear discordance commonly resulting from incomplete lineage sorting, adaptive introgression, or demographic disparities, the observed pattern in brood parasitic cuckoos might reflect biased host preferences between sexes. We hereby present the "Isolation by Gentes with Asymmetric Migration" model. According to this model, the maternal line differentiation of the common cuckoo in China is potentially driven by host preferences in females, whereas males maintained the integrity of the cuckoo species through random mating. To achieve this, cuckoo males could perform flexible migration among gentes or engage in early copulation with females before reaching the breeding sites, allowing female cuckoos to store sperm from various gentes. Future studies collecting additional samples from diverse cuckoo gentes with overlapping distribution and investigating the migratory and copulation patterns of each sex would enhance our understanding of sex-biased differentiation among cuckoo populations in China.

Rapid vertebrate speciation via isolation, bottlenecks, and drift

Speciation is often driven by selective processes like those associated with viability, mate choice, or local adaptation, and "speciation genes" have been identified in many eukaryotic lineages. In contrast, neutral processes are rarely considered as the primary drivers of speciation, especially over short evolutionary timeframes. Here, we describe a rapid vertebrate speciation event driven primarily by genetic drift. The White Sands pupfish (Cyprinodon tularosa) is endemic to New Mexico's Tularosa Basin where the species is currently managed as two Evolutionarily significant units (ESUs) and is of international conservation concern (Endangered). Whole-genome resequencing data from each ESU showed remarkably high and uniform levels of differentiation across the entire genome (global FST ≈ 0.40). Despite inhabiting ecologically dissimilar springs and streams, our whole-genome analysis revealed no discrete islands of divergence indicative of strong selection, even when we focused on an array of candidate genes. Demographic modeling of the joint allele frequency spectrum indicates the two ESUs split only ~4 to 5 kya and that both ESUs have undergone major bottlenecks within the last 2.5 millennia. Our results indicate the genome-wide disparities between the two ESUs are not driven by divergent selection but by neutral drift due to small population sizes, geographic isolation, and repeated bottlenecks. While rapid speciation is often driven by natural or sexual selection, here we show that isolation and drift have led to speciation within a few thousand generations. We discuss these evolutionary insights in light of the conservation management challenges they pose.

Phylogeographic Substructuring in the Southernmost Refugium of the European Common Frog Rana temporaria

Rana temporaria is one of the most widespread Palearctic brown frogs. We aimed to clarify distribution pattern of two main genetic clades in the understudied Balkan peninsula by using 16SrRNA and MT-CYTB sequences, already widely applied in analyses of populations from other parts of Europe, while focusing on the broad area along the Morava river (central Balkans) as a known gap in the species distribution. Additionally, we were interested in revealing the extent of haplotype diversity within the main genetic clades in the Balkans, particularly around the supposed suture zone. The results revealed a suture zone between the Western and Eastern Clades in the central part of the Balkan Peninsula. This indicated the existence of a historical barrier between the Balkan Mountain Belt and geographically close mountains surrounding the Vlasina Plateau (Rhodope/Serbian-Macedonian Massif). The overall observed haplotype diversity in populations of R. temporaria from the Balkan Peninsula seems high. Harboring both main genetic clades of R. temporaria qualifies the Balkan Peninsula as another important center of species' genetic diversity, as well as rich in unique haplotypes. This points out the necessity of applying conservation measures focused on the common European frog populations and habitats in this part of the species' distribution area.

Spatial genetic structure in European marbled newts revealed with target enrichment by sequence capture

European marbled newts come in two species that have abutting ranges. The northern species, Triturus marmoratus, is found in France and the northern part of the Iberian Peninsula, whereas the southern species, T. pygmaeus, is found in the southwestern corner of the Iberian Peninsula. We study the intraspecific genetic differentiation of the group because morphological data show geographical variation and because the Iberian Peninsula is a recognized center of speciation and intraspecific genetic diversity for all kinds of organisms, amphibians included. We use target enrichment by sequence capture to generate c. 7 k nuclear DNA markers. We observe limited genetic exchange between the species, which confirms their distinctiveness. Both species show substantial genetic structuring that is only in part mirrored by morphological variation. Genetically differentiated groups are found in the south (T. marmoratus) and west (T. pygmaeus) of the species ranges. Our observations highlight the position of the Iberian Peninsula as a hotspot for genetic differentiation.

Urea and glucose modulation during freezing exposure in three temperate frogs reveals specific targets in relation to climate

Amphibian diversity is most prominent in the warm and humid tropical and subtropical regions across the globe. Nonetheless, amphibians also inhabit high-altitude tropical mountains and regions at medium and high latitudes, exposing them to subzero temperatures and requiring behavioural or physiological adaptations to endure freezing events. While freeze tolerance has been predominantly reported in high-latitude zones where species endure prolonged freezing (several weeks or months), less is known about mid-latitudes amphibians exposed to occasional subzero temperatures. In this study, we employed a controlled ecological protocol, subjecting three frog species from the Iberian Peninsula (Rana parvipalmata, Epidalea calamita, and Pelobates cultripes) to a 2-h exposure to temperatures of -2 °C to investigate the accumulation of urea and glucose as physiological mechanisms associated with survival at freezing temperatures. Our results revealed a moderate response in the production of cryoprotectant metabolites under experimental freezing conditions, particularly urea, with notable findings in R. parvipalmata and E. calamita and no response in P. cultripes. However, no significant alterations in glucose concentrations were observed in any of the studied frog species. This relatively weak freezing tolerance response differs from the strong response exhibited by amphibians inhabiting high latitudes and enduring prolonged freezing conditions, suggesting potential reliance on behavioural adaptations to cope with occasional freezing episodes.

Phylogeography, hybrid zones and contemporary species boundaries in the south-eastern Australian smooth frogs (Anura: Myobatrachidae: Geocrinia)

Paleo-climatic fluctuations have driven episodic changes in species distributions, providing opportunities for populations to diverge in isolation and hybridise following secondary contact. Studies of phylogeographic diversity and patterns of gene flow across hybrid zones can provide insight into contemporary species boundaries and help to inform taxonomic and conservation inferences. Here we explore geographic diversity within the acoustically divergent yet morphologically conserved south-eastern Australian smooth frog complex and assess gene flow across a narrow hybrid zone using mitochondrial nucleotide sequences and nuclear genome-wide single nucleotide polymorphisms. Our analyses reveal the presence of an evolutionarily distinct taxon restricted to the Otway Plains and Ranges, Victoria, which forms a narrow (9-30 km wide), spatiotemporally stable (>50 years) hybrid zone with Geocrinia laevis, which we describe herein as a new species.

Population genetic structure and hybrid zone analyses for species delimitation in the Japanese toad (Bufo japonicus)

Hybridization following secondary contact may produce different outcomes depending on the extent to which genetic diversity and reproductive barriers have accumulated during isolation. The Japanese toad, Bufo japonicus, is distributed on the main islands of Japan. In the present study, we applied multiplexed inter-simple sequence repeat genotyping by sequencing to achieve the fine-scale resolution of the genetic cluster in B. j. japonicus and B. j. formosus. We also elucidated hybridization patterns and gene flow degrees across contact zones between the clusters identified. Using SNP data, we found four genetic clusters in B. j. japonicus and B. j. formosus and three contact zones of the cluster pairs among these four clusters. The two oldest diverged lineages, B. j. japonicus and B. j. formosus, formed a narrow contact zone consistent with species distinctiveness. Therefore, we recommend that these two subspecies be elevated to the species level. In contrast, the less diverged pairs of two clusters in B. j. japonicus and B. j. formosus, respectively, admixed over a hundred kilometers, suggesting that they have not yet developed strong reproductive isolation and need to be treated as conspecifics. These results will contribute to resolving taxonomic confusion in Japanese toads.

Mutations, sex, and genetic diversity: New arguments for partially directed evolution

A review of the theories of the existence of sexes, genetic diversity, and the distribution of mutations among organisms shows that all these concepts are not a product of random evolution and cannot be explained within the framework of Darwinism. Most mutations are the result of the genome acting on itself. This is an organized process that is implemented very differently in different species, in different places in the genome. Because of the fact that it is not random, this process must be directed and regulated, albeit with complex and not fully understood laws. This means that an additional reason must be included in order to model such mutations during evolution. The assumption of directionality must not only be explicitly included in evolutionary theory but must also occupy a central place in it. In this study an updated model of partially directed evolution is constructed, which is capable of qualitatively explaining the indicated features of evolution. Experiments are described that can confirm or disprove the proposed model.

Genomic phylogeography illuminates deep cyto-nuclear discordances in midwife toads (Alytes)

The advent of genomic methods allows us to revisit the evolutionary history of organismal groups for which robust phylogenies are still lacking, particularly in species complexes that frequently hybridize. In this study, we conduct RAD-sequencing (RAD-seq) analyses of midwife toads (genus Alytes), an iconic group of western Mediterranean amphibians famous for their parental care behavior, but equally infamous for the difficulties to reconstruct their evolutionary history. Through admixture and phylogenetic analyses of thousands of loci, we provide the most comprehensive phylogeographic framework for the A. obstetricans complex to date, as well as the first fully resolved phylogeny for the entire genus. As part of this effort, we carefully explore the influence of different sampling schemes and data filtering thresholds on tree reconstruction, showing that several, slightly different, yet robust topologies may be retrieved with small datasets obtained by stringent SNP calling parameters, especially when admixed individuals are included. In contrast, analyses of incomplete but larger datasets converged on the same phylogeny, irrespective of the reconstruction method used or the proportion of missing data. The Alytes tree features three Miocene-diverged clades corresponding to the proposed subgenera Ammoryctis (A. cisternasii), Baleaphryne (A. maurus, A. dickhilleni and A. muletensis), and Alytes (A. obstetricans complex). The latter consists of six evolutionary lineages, grouped into three clades of Pliocene origin, and currently delimited as two species: (1) A. almogavarii almogavarii and A. a. inigoi; (2) A. obstetricans obstetricans and A. o. pertinax; (3) A. o. boscai and an undescribed taxon (A. o. cf. boscai). These results contradict the mitochondrial tree, due to past mitochondrial captures in A. a. almogavarii (central Pyrenees) and A. o. boscai (central Iberia) by A. obstetricans ancestors during the Pleistocene. Patterns of admixture between subspecies appear far more extensive than previously assumed from microsatellites, causing nomenclatural uncertainties, and even underlying the reticulate evolution of one taxon (A. o. pertinax). All Ammoryctis and Baleaphryne species form shallow clades, so their taxonomy should remain stable. Amid the prevalence of cyto-nuclear discordance among terrestrial vertebrates and the usual lack of resolution of conventional nuclear markers, our study advocates for phylogeography based on next-generation sequencing, but also encourages properly exploring parameter space and sampling schemes when building and analyzing genomic datasets.

Environmental and topographic drivers of amphibian phylogenetic diversity and endemism in the Iberian Peninsula

Understanding the ecological and evolutionary processes driving biodiversity patterns and allowing their persistence is of utmost importance. Many hypotheses have been proposed to explain spatial diversity patterns, including water-energy availability, habitat heterogeneity, and historical climatic refugia. The main goal of this study is to identify if general spatial drivers of species diversity patterns of phylogenetic diversity (PD) and phylogenetic endemism (PE) at the global scale are also predictive of PD and PE at regional scales, using Iberian amphibians as a case study. Our main hypothesis assumes that topography along with contemporary and historical climate are drivers of phylogenetic diversity and endemism, but that the strength of these predictors may be weaker at the regional scale than it tends to be at the global scale. We mapped spatial patterns of Iberian amphibians' phylogenetic diversity and endemism, using previously published phylogenetic and distribution data. Furthermore, we compiled spatial data on topographic and climatic variables related to the water-energy availability, topography, and historical climatic instability hypotheses. To test our hypotheses, we used Spatial Autoregressive Models and selected the best model to explain diversity patterns based on Akaike Information Criterion. Our results show that, out of the variables tested in our study, water-energy availability and historical climate instability are the most important drivers of amphibian diversity in Iberia. However, as predicted, the strength of these predictors in our case study is weaker than it tends to be at global scales. Thus, additional drivers should also be investigated and we suggest caution when interpreting these predictors as surrogates for different components of diversity.

Testing assertions of widespread introgressive hybridization in a clade of neotropical toads with low mate selectivity (Rhinella granulosa species group)

Discordance between different genomic regions, often identified through multilocus sequencing of selected markers, presents particular difficulties in identifying historical processes which drive species diversity and boundaries. Mechanisms causing discordance, such as incomplete lineage sorting or introgression due to interspecific hybridization, are better identified based on population-level genomic datasets. In the toads of the Rhinella granulosa species group, patterns of mito-nuclear discordance and potential hybridization have been reported by several studies. However, these patterns were proposed based on few loci, such that alternative mechanisms behind gene-tree heterogeneity cannot be ruled out. Using genome-wide ddRADseq loci from a subset of species within this clade, we found only partial concordance between currently recognized species-level taxon boundaries and patterns of genetic structure. While most taxa within the R. granulosa group correspond to clades, genetic clustering analyses sometimes grouped distinct taxonomic units into a single cluster. Moreover, levels of admixture between inferred clusters were limited and restricted to a single taxon pair which is best explained by incomplete lineage sorting as opposed to introgressive hybridization, according to D-statistics results. These findings contradict previous assertions of widespread cryptic diversity and gene flow within the R. granulosa clade. Lastly, our analyses suggest that diversification events within the Rhinella granulosa group mostly dated back to the early Pliocene, being generally younger than species divergences in other closely related clades that present high levels of cross-species gene flow. This finding uniquely contradicts common assertions that this young clade of toads exhibits interspecific hybridization.

The nuanced nature of mesic refugia in arid landscapes: a tale of two peas

BACKGROUND AND AIMS

Understanding how genetic diversity is distributed and maintained within species is a central tenet of evolutionary and conservation biology, yet is understudied in arid regions of the globe. In temperate, glaciated environments, high genetic diversity in plant species is frequently found in refugial areas, which are often associated with southern non-glaciated landscapes. In arid, unglaciated environments, landscape features providing mesic conditions are likely to be refugia, although our understanding needs more refinement in these biomes. We test whether refugia and nuclear diversity hotspots occur in high-elevation, topographically complex areas for co-distributed shrubs (Petalostylis labicheoides and Indigofera monophylla; Fabaceae) in the ancient, arid Pilbara bioregion of north-western Australia.

METHODS

We conducted extensive sampling of the Pilbara (>1400 individuals from 62 widespread populations) to detect patterns in nuclear diversity and structure based on 13-16 microsatellite loci. Evidence of historical refugia was investigated based on patterns of diversity in three non-coding chloroplast (cp) sequence regions for approx. 240 individuals per species. Haplotype relationships were defined with median-joining networks and maximum likelihood phylogenetic trees.

KEY RESULTS

We found cpDNA evidence for a high-elevation refugium in P. labicheoides but not for I. monophylla that instead exhibited extraordinary haplotype diversity and evidence for persistence across a widespread area. Nuclear diversity hotspots occurred in, but were not exclusive to, high-elevation locations and extended to adjacent, low-elevation riparian areas in both species.

CONCLUSIONS

Phylogeographic refugia in arid environments may occur in high-elevation areas for some species but not all, and may be influenced by species-specific traits: a mesic montane refugium in P. labicheoides could be related to its preference for growth in water-gaining areas, while a lack of such evidence in I. monophylla could be related to maintenance of cpDNA diversity in a large soil seed bank and dynamic evolutionary history. Mesic environments created by the intersection of topographically complex landscapes with riparian zones can be contemporary reservoirs of genetic diversity in arid landscapes.

New insights on patterns of genetic admixture and phylogeographic history in Iberian high mountain populations of midwife toads

Multiple Quaternary glacial refugia in the Iberian Peninsula, commonly known as "refugia within refugia", allowed diverging populations to come into contact and admix, potentially boosting substantial mito-nuclear discordances. In this study, we employ a comprehensive set of mitochondrial and nuclear markers to shed light onto the drivers of geographical differentiation in Iberian high mountain populations of the midwife toads Alytes obstetricans and A. almogavarii from the Pyrenees, Picos de Europa and Guadarrama Mountains. In the three analysed mountain regions, we detected evidence of extensive mito-nuclear discordances and/or admixture between taxa. Clustering analyses identified three major divergent lineages in the Pyrenees (corresponding to the eastern, central and central-western Pyrenees), which possibly recurrently expanded and admixed during the succession of glacial-interglacial periods that characterised the Late Pleistocene, and that currently follow a ring-shaped diversification pattern. On the other hand, populations from the Picos de Europa mountains (NW Iberian Peninsula) showed a mitochondrial affinity to central-western Pyrenean populations and a nuclear affinity to populations from the central Iberian Peninsula, suggesting a likely admixed origin for Picos de Europa populations. Finally, populations from the Guadarrama Mountain Range (central Iberian Peninsula) were depleted of genetic diversity, possibly as a consequence of a recent epidemic of chytridiomycosis. This work highlights the complex evolutionary history that shaped the current genetic composition of high mountain populations, and underscores the importance of using a multilocus approach to better infer the dynamics of population divergence.

Phenology and plasticity can prevent adaptive clines in thermal tolerance across temperate mountains: The importance of the elevation-time axis

Critical thermal limits (CT_max and CT_min) decrease with elevation, with greater change in CT_min, and the risk to suffer heat and cold stress increasing at the gradient ends. A central prediction is that populations will adapt to the prevailing climatic conditions. Yet, reliable support for such expectation is scant because of the complexity of integrating phenotypic, molecular divergence and organism exposure. We examined intraspecific variation of CT_max and CT_min, neutral variation for 11 microsatellite loci, and micro- and macro-temperatures in larvae from 11 populations of the Galician common frog (Rana parvipalmata) across an elevational gradient, to assess (1) the existence of local adaptation through a P_ST-F_ST comparison, (2) the acclimation scope in both thermal limits, and (3) the vulnerability to suffer acute heat and cold thermal stress, measured at both macro- and microclimatic scales. Our study revealed significant microgeographic variation in CT_max and CT_min, and unexpected elevation gradients in pond temperatures. However, variation in CT_max and CT_min could not be attributed to selection because critical thermal limits were not correlated to elevation or temperatures. Differences in breeding phenology among populations resulted in exposure to higher and more variable temperatures at mid and high elevations. Accordingly, mid- and high-elevation populations had higher CT_max and CT_min plasticities than lowland populations, but not more extreme CT_max and CT_min. Thus, our results support the prediction that plasticity and phenological shifts may hinder local adaptation, promoting thermal niche conservatism. This may simply be a consequence of a coupled variation of reproductive timing with elevation (the "elevation-time axis" for temperature variation). Mid and high mountain populations of R. parvipalmata are more vulnerable to heat and cool impacts than lowland populations during the aquatic phase. All of this contradicts some of the existing predictions on adaptive thermal clines and vulnerability to climate change in elevational gradients.

Population genomics reveal deep divergence and strong geographical structure in gentians in the Hengduan Mountains

Understanding the evolutionary and ecological processes driving population differentiation and speciation can provide critical insights into the formation of biodiversity. Here, we examine the link between population genetic processes and biogeographic history underlying the generation of diversity in the Hengduan Mountains (HM), a region harboring a rich and dynamic flora. We used restriction site-associated DNA sequencing to generate 1,907 single-nucleotide polymorphisms (SNPs) and four-kb of plastid sequence in species of the Gentiana hexaphylla complex (Gentianaceae). We performed genetic clustering with spatial and non-spatial models, phylogenetic reconstructions, and ancestral range estimation, with the aim of addressing the processes influencing diversification of G. hexaphylla in the HM. We find the G. hexaphylla complex is characterized by geographic genetic structure with clusters corresponding to the South, North and the central HM. Phylogenetic reconstruction and pairwise F _ST analyses showed deep differentiation between Southern and Northern populations in the HM. The population in Mount Taibai exhibited the highest genetic similarity to the North HM. Ancestral range estimation indicated that the G. hexaphylla complex originated in the central HM and then diverged in the Pliocene and the Early Pleistocene, before dispersing widely, resulting in the current distinct lineages. Overall, we found deep genomic differentiation in the G. hexaphylla complex corresponds to geographic barriers to dispersal in the HM and highlights a critical role of the uplift of the Daxue Mountains and subsequent climatic fluctuations underlying diversification. The colonization of G. hexaphylla in the Mount Taibai region suggests directional dispersal between the alpine flora of the Qinling Mountains and the HM.

Phylogeography of Kisaura Ross (Trichoptera: Philopotamidae) of the Japanese Archipelago and the character displacement evolution observed in a secondary contact area between genetically differentiated intra-specific lineages

In this study we present evidence of an ongoing speciation event in the Japanese Archipelago. The Kisaura nozakii species complex (Philopotamidae) is a group of a small trichopteran insects, and it has been reported that this species has large differences between populations in the copulatory organs (genitalia) of males. Therefore, we conducted molecular phylogenetic analyses based on the mtDNA COI region (658 bp) and genome-wide nDNA single nucleotide polymorphisms (SNPs) (16 254 loci) using the K. nozakii species complex and a sister-species, K. borealis. We also conducted principal component analyses (PCA) based on the forewing length and seven morphological characteristics of the male genitalia. It was revealed that the K. nozakii species complex is composed of multiple allopatrically differentiated genetic lineages. However, two genetic lineages were distributed sympatrically/parapatrically in the western area of the Japanese Archipelago, and the occurrence of ‘character displacement’ in the male genitalia was detected in this area. We consider that the evolutionary event in this species complex has high potential to become a model case for elucidating the processes of speciation and species diversification.

Diversity of Land Snail Tribe Helicini (Gastropoda: Stylommatophora: Helicidae): Where Do We Stand after 20 Years of Sequencing Mitochondrial Markers?

Sequences of mitochondrial genes revolutionized the understanding of animal diversity and continue to be an important tool in biodiversity research. In the tribe Helicini, a prominent group of the western Palaearctic land snail fauna, mitochondrial data accumulating since the 2000s helped to newly delimit genera, inform species-level taxonomy and reconstruct past range dynamics. We combined the published data with own unpublished sequences and provide a detailed overview of what they revealed about the diversity of the group. The delimitation of Helix is revised by placing Helix godetiana back in the genus and new synonymies are suggested within the genera Codringtonia and Helix. The spatial distribution of intraspecific mitochondrial lineages of several species is shown for the first time. Comparisons between species reveal considerable variation in distribution patterns of intraspecific lineages, from broad postglacial distributions to regions with a fine-scale pattern of allopatric lineage replacement. To provide a baseline for further research and information for anyone re-using the data, we thoroughly discuss the gaps in the current dataset, focusing on both taxonomic and geographic coverage. Thanks to the wealth of data already amassed and the relative ease with which they can be obtained, mitochondrial sequences remain an important source of information on intraspecific diversity over large areas and taxa.

Diversity, distribution and molecular species delimitation in frogs and toads from the Eastern Palaearctic

Biodiversity analyses can greatly benefit from coherent species delimitation schemes and up-to-date distribution data. In this article, we have made the daring attempt to delimit and map described and undescribed lineages of anuran amphibians in the Eastern Palaearctic (EP) region in its broad sense. Through a literature review, we have evaluated the species status considering reproductive isolation and genetic divergence, combined with an extensive occurrence dataset (nearly 85k localities). Altogether 274 native species from 46 genera and ten families were retrieved, plus eight additional species introduced from other realms. Independent hotspots of species richness were concentrated in southern Tibet (Medog County), the circum-Sichuan Basin region, Taiwan, the Korean Peninsula and the main Japanese islands. Phylogeographic breaks responsible for recent in situ speciation events were shared around the Sichuan Mountains, across Honshu and between the Ryukyu Island groups, but not across shallow water bodies like the Yellow Sea and the Taiwan Strait. Anuran compositions suggested to restrict the zoogeographical limits of the EP to East Asia. In a rapidly evolving field, our study provides a checkpoint to appreciate patterns of species diversity in the EP under a single, spatially explicit, species delimitation framework that integrates phylogeographic data in taxonomic research.

Next-generation phylogeography of the banded newts (Ommatotriton): A phylogenetic hypothesis for three ancient species with geographically restricted interspecific gene flow and deep intraspecific genetic structure

Technological developments now make it possible to employ many markers for many individuals in a phylogeographic setting, even for taxa with large and complex genomes such as salamanders. The banded newt (genus Ommatotriton) from the Near East has been proposed to contain three species (O. nesterovi, O. ophryticus and O. vittatus) with unclear phylogenetic relationships, apparently limited interspecific gene flow and deep intraspecific geographic mtDNA structure. We use parallel tagged amplicon sequencing to obtain 177 nuclear DNA markers for 35 banded newts sampled throughout the range. We determine population structure (with Bayesian clustering and principal component analysis), interspecific gene flow (by determining the distribution of species-diagnostic alleles) and phylogenetic relationships (by maximum likelihood inference of concatenated sequence data and based on a summary-coalescent approach). We confirm that the three proposed species are genetically distinct. A sister relationship between O. nesterovi and O. ophryticus is suggested. We find evidence for introgression between O. nesterovi and O. ophryticus, but this is geographically limited. Intraspecific structuring is extensive, with the only recognized banded newt subspecies, O. vittatus cilicensis, representing the most distinct lineage below the species level. While mtDNA mostly mirrors the pattern observed in nuclear DNA, all banded newt species show mito-nuclear discordance as well.

Complex and divergent histories gave rise to genome-wide divergence patterns amongst European whitefish (Coregonus lavaretus)

Pleistocene glaciations dramatically affected species distribution in regions that were impacted by ice cover and subsequent postglacial range expansion impacted contemporary biodiversity in complex ways. The European whitefish, Coregonus lavaretus, is a widely distributed salmonid fish species on mainland Europe, but in Britain it has only seven native populations, all of which are found on the western extremes of the island. The origins and colonization routes of the species into Britain are unknown but likely contributed to contemporary genetic patterns and regional uniqueness. Here, we used up to 25,751 genome‐wide polymorphic loci to reconstruct the history and to discern the demographic and evolutionary forces underpinning divergence between British populations. Overall, we found lower genetic diversity in Scottish populations but high differentiation (F_ST = 0.433–0.712) from the English/Welsh and other European populations. Differentiation was elevated genome‐wide rather than in particular genomic regions. Demographic modelling supported a postglacial colonization into western Scotland from northern refugia and a separate colonization route for the English/Welsh populations from southern refugia, with these two groups having been separated for more than ca. 50 Ky. We found cyto‐nuclear discordance at a European scale, with the Scottish populations clustering closely with Baltic population in the mtDNA analysis but not in the nuclear data, and with the Norwegian and Alpine populations displaying the same mtDNA haplotype but being distantly related in the nuclear tree. These findings suggest that neutral processes, primarily drift and regionally distinct pre‐glacial evolutionary histories, are important drivers of genomic divergence in British populations of European whitefish. This sheds new light on the establishment of the native British freshwater fauna after the last glacial maximum.

Phylogeographic analysis of shrubby beardtongues reveals range expansions during the Last Glacial Maximum and implicates the Klamath Mountains as a hotspot for hybridization

Quaternary glacial cycles often altered species' geographic distributions, which in turn altered the geographic structure of species' genetic diversity. In many cases, glacial expansion forced species in temperate climates to contract their ranges and reside in small pockets of suitable habitat (refugia), where they were likely to interact closely with other species, setting the stage for potential gene exchange. These introgression events, in turn, would have degraded species boundaries, making the inference of phylogenetic relationships challenging. Using high-throughput sequence data, we employed a combination of species distribution models and hybridization tests to assess the effect of glaciation on the geographic distributions, phylogenetic relationships, and patterns of gene flow of five species of Penstemon subgenus Dasanthera, long-lived shrubby angiosperms distributed throughout the Pacific Northwest of North America. Surprisingly, we found that rather than reducing their ranges to small refugia, most Penstemon subgenus Dasanthera species experienced increased suitable habitat during the Last Glacial Maximum relative to the present day. We also found substantial evidence for gene exchange between species, with the bulk of introgression events occurring in or near the Klamath Mountains of southwestern Oregon and northwestern California. Subsequently, our phylogenetic inference reveals blurred taxonomic boundaries in the Klamath Mountains, where introgression is most prevalent. Our results question the classical paradigm of temperate species' responses to glaciation and highlight the importance of contextualizing phylogenetic inference with species' histories of introgression.

Landscape resistance constrains hybridization across contact zones in a reproductively and morphologically polymorphic salamander

Explicitly accounting for phenotypic differentiation together with environmental heterogeneity is crucial to understand the evolutionary dynamics in hybrid zones. Species showing intra-specific variation in phenotypic traits that meet across environmentally heterogeneous regions constitute excellent natural settings to study the role of phenotypic differentiation and environmental factors in shaping the spatial extent and patterns of admixture in hybrid zones. We studied three environmentally distinct contact zones where morphologically and reproductively divergent subspecies of Salamandra salamandra co-occur: the pueriparous S. s. bernardezi that is mostly parapatric to its three larviparous subspecies neighbours. We used a landscape genetics framework to: (i) characterise the spatial location and extent of each contact zone; (ii) assess patterns of introgression and hybridization between subspecies pairs; and (iii) examine the role of environmental heterogeneity in the evolutionary dynamics of hybrid zones. We found high levels of introgression between parity modes, and between distinct phenotypes, thus demonstrating the evolution to pueriparity alone or morphological differentiation do not lead to reproductive isolation between these highly divergent S. salamandra morphotypes. However, we detected substantial variation in patterns of hybridization across contact zones, being lower in the contact zone located on a topographically complex area. We highlight the importance of accounting for spatial environmental heterogeneity when studying evolutionary dynamics of hybrid zones.

Stability in the Balkans: phylogeography of the endemic Greek stream frog, Rana graeca

We still have little knowledge concerning the phylogeography of amphibians and reptiles from the Balkan Peninsula compared with the other two Mediterranean peninsulas. This raises concerns for endemic taxa from these peninsulas, because it might interfere with further conservation efforts. Here we focus on the endemic Greek stream frog (Rana graeca) and reconstruct its biogeography and evolutionary history. Using four genetic markers (Cytb, 16S, COI and BDNF) in > 350 sequences covering the whole distribution range, we conducted phylogenetic, demographic and ecological niche analyses, which revealed the phylogeography of this species. Surprisingly, this examination of R. graeca reveals a very shallow level of intraspecific genetic variability through the Balkans, with two main, statistically supported lineages having a partly sympatric distribution. The most variable marker was Cytb, which showed 19 haplotypes in 123 analysed sequences in the whole species distribution area. Here presented genetic data, together with the environmental niche projection and demographic analyses suggest that R. graeca was probably affected only marginally by climatic oscillations, with the Hellenides as the most suitable area for the occurrence of the species in different geological periods. This is consistent with the observed genetic diversity, which is mostly related to these mountains. Although the species shows a certain level of phenotypic variability and ecological preferences, this might be related to species plasticity affected by the micro-climatic conditions in small areas, which merits further research. Comparing phylogeography of other amphibian and reptile species in the Balkans, we showed that the observed pattern represents a new view on the phylogeography of the Balkan herpetofauna.

Climatic Refugia and Geographical Isolation Contribute to the Speciation and Genetic Divergence in Himalayan-Hengduan Tree Peonies (Paeonia delavayi and Paeonia ludlowii)

Himalaya and Hengduan Mountains (HHM) is a biodiversity hotspot, and very rich in endemic species. Previous phylogeographical studies proposed different hypotheses (vicariance and climate-driven speciation) in explaining diversification and the observed pattern of extant biodiversity, but it is likely that taxa are forming in this area in species-specific ways. Here, we reexplored the phylogenetic relationship and tested the corresponding hypotheses within Paeonia subsect. Delavayanae composed of one widespread species (Paeonia delavayi) and the other geographically confined species (Paeonia ludlowii). We gathered genetic variation data at three chloroplast DNA fragments and one nuclear gene from 335 individuals of 34 populations sampled from HHM. We performed a combination of population genetic summary statistics, isolation-with-migration divergence models, isolation by environment, and demographic history analyses. We found evidence for the current taxonomic treatment that P. ludlowii and P. delavayi are two different species with significant genetic differentiation. The significant isolation by environment was revealed within all sampled populations but genetic distances only explained by geographical distances within P. delavayi populations. The results of population divergence models and demographic history analyses indicated a progenitor–derivative relationship and the Late Quaternary divergence without gene flow between them. The coalescence of all sampled cpDNA haplotypes could date to the Late Miocene, and P. delavayi populations probably underwent a severe bottleneck in population size during the last glacial period. Genetic variation in Paeonia subsect. Delavayanae is associated with geographical and environmental distances. These findings point to the importance of geological and climatic changes as causes of the speciation event and lineage diversification within Paeonia subsect. Delavayanae.

Slow and steady wins the race: contrasted phylogeographic signatures in two Alpine amphibians

A deeper phylogeographic structure is expected for slow-dispersing habitat specialists compared to widespread adaptable species, especially in topographically complex regions. We tested this classic assumption by comparing the genomic (RAD-sequencing) phylogeographies of two amphibians inhabiting the Swiss Alps: the mobile, cosmopolitan common frog (Rana temporaria) against the stationary, mountain endemic Alpine salamander (Salamandra atra). Our results ran opposite of predictions: the frog displayed significantly higher genetic divergences and lower within-population variation compared to the salamander. This implies a prominent role for their distinctive glacial histories in shaping intraspecific diversity and structure: diversification and recolonization from several circum-Alpine micro-refugia for the frog versus a single refugium for the salamander, potentially combined with better population connectivity and stability. These striking differences emphasize the great variability of phylogeographic responses to the Quaternary glaciations, hence the complexity to predict general patterns of genetic diversity at the regional scale, and the forces that underlie them.

Phylogenomic inference of species and subspecies diversity in the Palearctic salamander genus Salamandra

The salamander genus Salamandra is widespread across Europe, North Africa, and the Near East and is renowned for its conspicuous and polymorphic colouration and diversity of reproductive modes. The phylogenetic relationships within the genus, and especially in the highly polymorphic species S. salamandra, have been very challenging to elucidate, leaving its real evolutionary history and classification at species and subspecies levels a topic of debate and contention. However, the distribution of diversity and species delimitation within the genus are critically important for identifying evolutionarily significant units for conservation and management, especially in light of threats posed by the pathogenic chytrid fungus Batrachochytrium salamandrivorans that is causing massive declines of S. salamandra populations in central Europe. Here, we conducted a phylogenomic analysis from across the taxonomic and geographic breadth of the genus Salamandra in its entire range. Bayesian, maximum likelihood and network-based phylogenetic analyses of up to 4905 ddRADseq-loci (294,300 nucleotides of sequence) supported the distinctiveness of all currently recognised species (Salamandra algira, S. atra, S. corsica, S. infraimmaculata, S. lanzai, and S. salamandra), and all five species for which we have multiple exemplars were confirmed as monophyletic. Within S. salamandra, two main clades can be distinguished: one clade with the Apenninic subspecies S. s. gigliolii nested within the Iberian S. s. bernardezi/fastuosa; and a second clade comprising all other Iberian, Central and East European subspecies. Our analyses revealed that some of the currently recognized subspecies of S. salamandra are paraphyletic and may require taxonomic revision, with the Central- and Eastern-European subspecies all being poorly differentiated at the analysed genomic markers. Salamandra s. longirostris - sometimes considered a separate species - was nested within S. salamandra, consistent with its subspecies status. The relationships identified within and between Salamandra species provide valuable context for future systematic and biogeographic studies, and help elucidate critical evolutionary units for conservation and taxonomy.

Pleistocene climate fluctuations drove demographic history of African golden wolves (Canis lupaster)

Pleistocene climate change impacted entire ecosystems throughout the world. In the northern hemisphere, the distribution of Arctic species expanded during glacial periods, while more temperate and mesic species contracted into climatic refugia, where isolation drove genetic divergence. Cycles of local cooling and warming in the Sahara region of northern Africa caused repeated contractions and expansions of savannah-like environments which connected mesic species isolated in refugia during interglacial times, possibly driving population expansions and contractions; divergence and geneflow in the associated fauna. Here, we use whole genome sequences of African golden wolves (Canis lupaster), a generalist mesopredator with a wide distribution in northern Africa to estimate their demographic history and past episodes of geneflow. We detect a correlation between divergence times and cycles of increased aridity-associated Pleistocene glacial cycles. A complex demographic history with responses to local climate change in different lineages was found, including a relict lineage north of the High Atlas Mountains of Morocco that has been isolated for more than 18,000 years, possibly a distinct ecotype.

Genetic Divergence Across Glacial Refugia Despite Interglacial Gene Flow in a Crested Newt

MtDNA-based phylogeography has illuminated the impact of the Pleistocene Ice Age on species distribution dynamics and the build-up of genetic divergence. The well-known shortcomings of mtDNA in biogeographical inference can be compensated by integrating multilocus data and species distribution modelling into phylogeography. We re-visit the phylogeography of the Italian crested newt (Triturus carnifex), a species distributed in two of Europe’s main glacial refugia, the Balkan and Italian Peninsulas. While a new 51 nuclear DNA marker dataset supports the existence of three lineages previously suggested by mtDNA (Balkan, northern Italy and southern Italy), the nuclear DNA dataset also provides improved resolution where these lineages have obtained secondary contact. We observe geographically restricted admixture at the contact between the Balkan and northern Italy gene pools and identify a potential mtDNA ghost lineage here. At the contact between the northern and southern Italy gene pools we find admixture over a broader area, as well as asymmetric mtDNA introgression. Our species distribution model is in agreement with a distribution restricted to distinct refugia during Pleistocene glacial cycles and postglacial expansion with secondary contact. Our study supports: (1) the relevance of the north-western Balkan Peninsula as a discrete glacial refugium; (2) the importance of north-eastern Italy and the northern Apennine as suture zones; and (3) the applicability of a refugia-within-refugia scenario within the Italian Peninsula.

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.

Italian odonates in the Pandora's box: A comprehensive DNA barcoding inventory shows taxonomic warnings at the Holarctic scale

The Odonata are considered among the most endangered freshwater faunal taxa. Their DNA-based monitoring relies on validated reference data sets that are often lacking or do not cover important biogeographical centres of diversification. This study presents the results of a DNA barcoding campaign on Odonata, based on the standard 658-bp 5' end region of the mitochondrial COI gene, involving the collection of 812 specimens (409 of which barcoded) from peninsular Italy and its main islands (328 localities), belonging to all the 88 species (31 Zygoptera and 57 Anisoptera) known from the country. Additional BOLD and GenBank data from Holarctic samples expanded the data set to 1,294 DNA barcodes. A multi-approach species delimitation analysis involving two distance (OT and ABGD) and four tree-based (PTP, MPTP, GMYC and bGMYC) methods was used to explore these data. Of the 88 investigated morphospecies, 75 (85%) unequivocally corresponded to distinct molecular operational units, whereas the remaining ones were classified as 'warnings' (i.e. showing a mismatch between morphospecies assignment and DNA-based species delimitation). These results are in contrast with other DNA barcoding studies on Odonata showing up to 95% of identification success. The species causing warnings were grouped into three categories depending on if they showed low, high or mixed genetic divergence patterns. The analysis of haplotype networks revealed unexpected intraspecific complexity at the Italian, Palearctic and Holarctic scale, possibly indicating the occurrence of cryptic species. Overall, this study provides new insights into the taxonomy of odonates and a valuable basis for future DNA and eDNA-based monitoring studies.

The toad's warts: Discordance creates bumpy expectations of mitochondrial-nuclear evolution between species

Discordance between the mitochondrial and nuclear genomes is a prevalent phenomenon in nature, in which the underlying processes responsible are considered to be important in shaping genetic variation in natural populations. Among the evolutionary processes that best explain such genomic mismatches incomplete lineage sorting and introgression are commonly identified, however, many studies are unable to distinguish between these hypotheses, which has become a major challenge in the field. In this issue of Molecular Ecology, Firneno et al. (2020) present an elegant exploration of mitochondrial-nuclear discordance in Mesoamerican toads. Integrating genome-scale and spatial data to test between these hypotheses within an empirical model testing framework, they find strong support that incomplete lineage sorting explains the observed discordance. Their work, along with many previous articles in Molecular Ecology, highlights the commonality of mito-nuclear discordance among species despite the expectations of tightly concerted mitochondrial and nuclear genome evolution. It is increasingly clear that the nuclear genomes of many species are (at least for short periods of evolutionary time) functionally compatible with multiple, divergent mitochondrial haplotypes. As such, we suggest future research not only seeks to understand the processes causing spatial mito-nuclear discordance (e.g. incomplete lineage sorting, introgression), but also explores those that maintain discordance through time and space (e.g. relaxed selection on mito-nuclear interactions, heterozygosity, population demographics). We also discuss the vital role that taxonomy plays in interpreting patterns of mito-nuclear discordance when data-consistent yet differing taxonomies are used, such as treating allopatrically distributed taxa as multiple isolated populations versus multiple micro-endemic species.

Reconstruction of past distribution for the Mongolian toad, Strauchbufo raddei (Anura: Bufonidae) using environmental modeling

The use of ecological models enables determining the current distribution of species, but also their past distribution when matching climatic conditions are available. In specific cases, they can also be used to determine the likelihood of fossils to belong to the same species-under the hypothesis that all individuals of a species have the same ecological requirements. Here, using environmental modeling, we reconstructed the distribution of the Mongolian toad, Strauchbufo raddei, since the Last Glacial Maximum and thus covering the time period between the Late Pleistocene and the Holocene. We found the range of the species to have shifted over time, with the LGM population clustered around the current southern range of the species, before expanding east and north during the Pleistocene, and reaching the current range since the mid-Holocene. Finally, we determined that the ecological conditions during the life-time of the mid-Pleistocene fossils attributed to the species in Europe were too different from the one of the extant species or fossils occurring at the same period in Asia to belong to the same species.

Extensive cytonuclear discordance in a crested newt from the Balkan Peninsula glacial refugium

Integration of multilocus data and species distribution modelling into phylogeography allows mitochondrial DNA (mtDNA)-based scenarios to be fine-tuned. We address the question of whether extensive mtDNA substructuring in the crested newt Triturus macedonicus from the Balkan Peninsula is matched in the nuclear genome. We determine the intraspecific population structure based on 52 nuclear DNA markers and project a species distribution model on climate layers for the Last Glacial Maximum. We show that T. macedonicus accumulated nuclear DNA population structure in an area predicted to have been climatically stable during the Pleistocene, with four nuclear DNA groups in the western part of the species range. The distribution of these nuclear DNA groups shows little agreement with that of mtDNA structuring, which shows three highly distinct species-specific clades and a fourth one introgressed from another crested newt species. This cytonuclear discordance conveys that historical biogeographical scenarios based on mtDNA exclusively should be interpreted with caution. Our findings further highlight the important role the Balkan Peninsula has played in the evolution and preservation of European biodiversity.