Long-term climatic stability drives accumulation and maintenance of divergent freshwater fish lineages in a temperate biodiversity hotspot

Anthropogenic climate change is forecast to drive regional climate disruption and instability across the globe. These impacts are likely to be exacerbated within biodiversity hotspots, both due to the greater potential for species loss but also to the possibility that endemic lineages might not have experienced significant climatic variation in the past, limiting their evolutionary potential to respond to rapid climate change. We assessed the role of climatic stability on the accumulation and persistence of lineages in an obligate freshwater fish group endemic to the southwest Western Australia (SWWA) biodiversity hotspot. Using 19,426 genomic (ddRAD-seq) markers and species distribution modelling, we explored the phylogeographic history of western (Nannoperca vittata) and little (Nannoperca pygmaea) pygmy perches, assessing population divergence and phylogenetic relationships, delimiting species and estimating changes in species distributions from the Pliocene to 2100. We identified two deep phylogroups comprising three divergent clusters, which showed no historical connectivity since the Pliocene. We conservatively suggest these represent three isolated species with additional intraspecific structure within one widespread species. All lineages showed long-term patterns of isolation and persistence owing to climatic stability but with significant range contractions likely under future climate change. Our results highlighted the role of climatic stability in allowing the persistence of isolated lineages in the SWWA. This biodiversity hotspot is under compounding threat from ongoing climate change and habitat modification, which may further threaten previously undetected cryptic diversity across the region.

Congruent patterns of cryptic cladogenesis revealed using RADseq and Sanger sequencing in a velvet worm species complex (Onychophora: Peripatopsidae: Peripatopsis sedgwicki)

In the present study, first generation DNA sequencing (mitochondrial cytochrome c oxidase subunit one, COI) and reduced-representative genomic RADseq data were used to understand the patterns and processes of diversification of the velvet worm, Peripatopsis sedgwicki species complex across its distribution range in South Africa. For the RADseq data, three datasets (two primary and one supplementary) were generated corresponding to 1,259-11,468 SNPs, in order to assess the diversity and phylogeography of the species complex. Tree topologies for the two primary datasets were inferred using maximum likelihood and Bayesian inferences methods. Phylogenetic analyses using the COI datasets retrieved four distinct, well-supported clades within the species complex. Five species delimitation methods applied to the COI data (ASAP, bPTP, bGMYC, STACEY and iBPP) all showed support for the distinction of the Fort Fordyce Nature Reserve specimens. In the main P. sedgwicki species complex, the species delimitation methods revealed a variable number of operational taxonomic units and overestimated the number of putative taxa. Divergence time estimates coupled with the geographic exclusivity of species and phylogeographic results suggest recent cladogenesis during the Plio/Pleistocene. The RADseq data were subjected to a principal components analysis and a discriminant analysis of principal components, under a maximum-likelihood framework. The latter results corroborate the four main clades observed using the COI data, however, applying additional filtering revealed additional diversity. The high overall congruence observed between the RADseq data and COI data suggest that first generation sequence data remain a cheap and effective method for evolutionary studies, although RADseq does provide a far greater resolution of contemporary temporo-spatial patterns.

Establishing species boundaries in Bornean geckos

Species delimitation using mitochondrial DNA (mtDNA) remains an important and accessible approach for discovering and delimiting species. However, delimiting species with a single locus (e.g. DNA barcoding) is biased towards overestimating species diversity. The highly diverse gecko genus Cyrtodactylus is one such group where delimitation using mtDNA remains the paradigm. In this study, we use genomic data to test putative species boundaries established using mtDNA within three recognized species of Cyrtodactylus on the island of Borneo. We predict that multi-locus genomic data will estimate fewer species than mtDNA, which could have important ramifications for the species diversity within the genus. We aim to (i) investigate the correspondence between species delimitations using mtDNA and genomic data, (ii) infer species trees for each target species, and (iii) quantify gene flow and identify migration patterns to assess population connectivity. We find that species diversity is overestimated and that species boundaries differ between mtDNA and nuclear data. This underscores the value of using genomic data to reassess mtDNA-based species delimitations for taxa lacking clear species boundaries. We expect the number of recognized species within Cyrtodactylus to continue increasing, but, when possible, genomic data should be included to inform more accurate species boundaries.

Endless forms most frustrating: disentangling species boundaries in the Ramalina decipiens group (Lecanoromycetes, Ascomycota), with the description of six new species and a key to the group

Oceanic islands have been recognized as natural laboratories in which to study a great variety of evolutionary processes. One such process is evolutionary radiations, the diversification of a single ancestor into a number of species that inhabit different environments and differ in the traits that allow them to exploit those environments. The factors that drive evolutionary radiations have been studied for decades in charismatic organisms such as birds or lizards, but are lacking in lichen-forming fungi, despite recent reports of some lineages showing diversification patterns congruent with radiation. Here we propose the Ramalina decipiens group as a model system in which to carry out such studies. This group is currently thought to be comprised of five saxicolous species, all of them endemic to the Macaronesian region (the Azores, Madeira, Selvagens, Canary and Cape Verde islands). Three species are single-island endemics (a rare geographic distribution pattern in lichens), whereas two are widespread and show extreme morphological variation. The latter are suspected to harbor unrecognized species-level lineages. In order to use the Ramalina decipiens group as a model system it is necessary to resolve the group's phylogeny and to clarify its species boundaries. In this study we attempt to do so following an integrative taxonomy approach. We constructed a phylogenetic tree based on six molecular markers, four of which are newly developed and generated competing species hypotheses based on molecular (species discovery strategies based on both single locus and multilocus datasets) and phenotypic data (unsupervised clustering algorithms based on morphology, secondary chemistry and geographic origin). We found that taxonomic diversity in the Ramalina decipiens group has been highly underestimated in previous studies. In consequence, we describe six new species, most of them single-island endemics and provide a key to the group. Phylogenetic relationships among species have been reconstructed with almost full support which, coupled with the endemic character of the group, makes it an excellent system for the study of island radiations in lichen-forming fungi. Citation: Blázquez M, Pérez-Vargas I, Garrido-Benavent I, et al. 2024. Endless forms most frustrating: disentangling species boundaries in the Ramalina decipiens group (Lecanoromycetes, Ascomycota), with the description of six new species and a key to the group. Persoonia 52: 44-93. https://doi.org/10.3767/persoonia.2024.52.03 .

DNA barcoding reveals cryptic species in the sea slater Ligia italica (Crustacea, Isopoda) from Tunisia

Barcoding studies have provided significant insights into phylogenetic relationships among species belonging to the genus Ligia (Crustacea, Isopoda). Herein the diversity of the Italian sea slater Ligia italica from Tunisia is studied for the first time. Samples were collected from 18 localities in Tunisia, and the analysis included previously published sequences from Italy and Greece available in GenBank. Bayesian and Maximum Likelihood phylogenetic analyses were carried out using a fragment of the mitochondrial COI gene. Putative cryptic species were explored using the 'barcode gap' approach in the software ASAP. A genetic landscape shape analysis was carried out using the program Alleles in Space. The analyses revealed highly divergent and well-supported clades of L. italica dispersed across Tunisia (Clades A1 and A2), Greece (Clade B) and Italy (Clades C1 and C2). High genetic dissimilarity among clades suggested that L. italica constitute a cryptic species complex. Divergence among different L. italica lineages (Clades A, B and C) occurred around 7-4.5 Ma. The detected high genetic distances among clades did not result from atypical mitochondrial DNAs or intracellular infection by Wolbachia bacteria. The complex history of the Mediterranean Sea appears to have played a significant role in shaping the phylogeographic pattern of Ligia italica. Additional morphological and molecular studies are needed to confirm the existence of cryptic species in Ligia italica in Mediterranean.

Pleistocene glaciation advances the cryptic speciation of Stellera chamaejasme L. in a major biodiversity hotspot

The mountains of Southwest China comprise a significant large mountain range and biodiversity hotspot imperiled by global climate change. The high species diversity in this mountain system has long been attributed to a complex set of factors, and recent large-scale macroevolutionary investigations have placed a broad timeline on plant diversification that stretches from 10 million years ago (Mya) to the present. Despite our increasing understanding of the temporal mode of speciation, finer-scale population-level investigations are lacking to better refine these temporal trends and illuminate the abiotic and biotic influences of cryptic speciation. This is largely due to the dearth of organismal sampling among closely related species and populations, spanning the incredible size and topological heterogeneity of this region. Our study dives into these evolutionary dynamics of speciation using genomic and eco-morphological data of Stellera chamaejasme L. We identified four previously unrecognized cryptic species having indistinct morphological traits and large metapopulation of evolving lineages, suggesting a more recent diversification (~2.67-0.90 Mya), largely influenced by Pleistocene glaciation and biotic factors. These factors likely influenced allopatric speciation and advocated cyclical warming-cooling episodes along elevational gradients during the Pleistocene. The study refines the evolutionary timeline to be much younger than previously implicated and raises the concern that projected future warming may influence the alpine species diversity, necessitating increased conservation efforts.

DNA metabarcoding reveals seasonal changes in diet composition across four arthropod-eating lizard species (Phrynosomatidae: Sceloporus)

Diet composition and its ecological drivers are rarely investigated in coexisting closely related species. We used a molecular approach to characterize the seasonal variation in diet composition in four spiny lizard species inhabiting a mountainous ecosystem. DNA metabarcoding revealed that the lizards Sceloporus aeneus, S. bicanthalis, S. grammicus, and S. spinosus mostly consumed arthropods of the orders Hemiptera, Araneae, Hymenoptera, and Coleoptera. The terrestrial lizards S. aeneus and S. bicanthalis mostly predated ants and spiders, whereas the arboreal-saxicolous S. grammicus and saxicolous S. spinosus largely consumed grasshoppers and leafhoppers. The taxonomic and phylogenetic diversity of the prey was higher during the dry season than the rainy season, likely because reduced prey availability in the dry season forced lizards to diversify their diets to meet their nutritional demands. Dietary and phylogenetic composition varied seasonally depending on the species, but only dietary composition varied with altitude. Seasonal dietary turnover was greater in S. spinosus than in S. bicanthalis, suggesting site-specific seasonal variability in prey availability; no other differences among species were observed. S. bicanthalis, which lives at the highest altitude in our study site, displayed interseasonal variation in diet breadth. Dietary differences were correlated with the species' feeding strategies and elevational distribution, which likely contributed to the coexistence of these lizard species in the studied geographic area and beyond.

Species delimitation with limited sampling: An example from rare trapdoor spider genus Cyclocosmia (Mygalomorphae, Halonoproctidae)

The outcome of species delimitation depends on many factors, including conceptual framework, study design, data availability, methodology employed and subjective decision making. Obtaining sufficient taxon sampling in endangered or rare taxa might be difficult, particularly when non-lethal tissue collection cannot be utilized. The need to avoid overexploitation of the natural populations may thus limit methodological framework available for downstream data analyses and bias the results. We test species boundaries in rare North American trapdoor spider genus Cyclocosmia Ausserer (1871) inhabiting the Southern Coastal Plain biodiversity hotspot with the use of genomic data and two multispecies coalescent model methods. We evaluate the performance of each methodology within a limited sampling framework. To mitigate the risk of species over splitting, common in taxa with highly structured populations, we subsequently implement a species validation step via genealogical diversification index (gdi), which accounts for both genetic isolation and gene flow. We delimited eight geographically restricted lineages within sampled North American Cyclocosmia, suggesting that major river drainages in the region are likely barriers to dispersal. Our results suggest that utilizing BPP in the species discovery step might be a good option for datasets comprising hundreds of loci, but fewer individuals, which may be a common scenario for rare taxa. However, we also show that such results should be validated via gdi, in order to avoid over splitting.

Species Diagnosis and DNA Taxonomy

The use of DNA has helped to improve and speed up species identification and delimitation. However, it also provides new challenges to taxonomists. Incongruence of outcome from various markers and delimitation methods, bias from sampling and skewed species distribution, implemented models, and the choice of methods/priors may mislead results and also may, in conclusion, increase elements of subjectivity in species taxonomy. The lack of direct diagnostic outcome from most contemporary molecular delimitation approaches and the need for a reference to existing and best sampled trait reference systems reveal the need for refining the criteria of species diagnosis and diagnosability in the current framework of nomenclature codes and good practices to avoid nomenclatorial instability, parallel taxonomies, and consequently more and new taxonomic impediment.

Unraveling phylogenetic relationships and species boundaries in the arid adapted Gerbillus rodents (Muridae: Gerbillinae) by RAD-seq data

Gerbillus is one of the most speciose genera among rodents, with ca. 51 recognized species. Previous attempts to reconstruct the evolutionary history of Gerbillus mainly relied on the mitochondrial cyt-b marker as a source of phylogenetic information. In this study, we utilize RAD-seq genomic data from 37 specimens representing 11 species to reconstruct the phylogenetic tree for Gerbillus, applying concatenation and coalescence methods. We identified four highly supported clades corresponding to the traditionally recognized subgenera: Dipodillus, Gerbillus, Hendecapleura and Monodia. Only two uncertain branches were detected in the resulting trees, with one leading to diversification of the main lineages in the genus, recognized by quartet sampling analysis as uncertain due to possible introgression. We also examined species boundaries for four pairs of sister taxa, including potentially new species from Morocco, using SNAPP. The results strongly supported a speciation model in which all taxa are treated as separate species. The dating analyses confirmed the Plio-Pleistocene diversification of the genus, with the uncertain branch coinciding with the beginning of aridification of the Sahara at the the Plio-Pleistocene boundary. This study aligns well with the earlier analyses based on the cyt-b marker, reaffirming its suitability as an adequate marker for estimating genetic diversity in Gerbillus.

Genomic evidence sheds new light on phylogeny of Rhabdophis nuchalis (sensu lato) complex (Serpentes: Natricidae)

Stable taxonomy and robust phylogeny are essential for the evolution and conservation of organisms. The Rhabdophis nuchalis (sensu lato) complex presently contains three species (R. nuchalis, R. chiwen, R. pentasupralabialis). Although several studies have explored the diversity and phylogeography of this group, certain issues related to systematics and taxonomy remain unresolved. Here, based on genome-wide data, including single nucleotide polymorphisms (SNPs) generated from ddRAD-seq and mitochondrial DNA (mtDNA), we re-evaluated the phylogenetic relationships and cryptic diversity of this species group. Our results are generally consistent with previous studies but provide some new insights. Phylogenetic relationship reconstruction based on SNPs and mtDNA revealed that three species in the R. nuchalis (sensu lato) complex did not form a monophyly but each species is well supported as monophyletic lineage in SNP-based analyses. Population structure analyses showed genetic admixture between several species pairs. Additionally, the population in eastern Yunnan, China, was identified as a potential cryptic species and thus described as a new species based on morphological data. From our results and previous studies, we redefined the distribution boundary for each species in the R. nuchalis (sensu lato) species complex.

Whole-genome sequencing reveals adaptations of hairy-footed jerboas (Dipus, Dipodidae) to diverse desert environments

BACKGROUND

Environmental conditions vary among deserts across the world, spanning from hyper-arid to high-elevation deserts. However, prior genomic studies on desert adaptation have focused on desert and non-desert comparisons overlooking the complexity of conditions within deserts. Focusing on the adaptation mechanisms to diverse desert environments will advance our understanding of how species adapt to extreme desert environments. The hairy-footed jerboas are well adapted to diverse desert environments, inhabiting high-altitude arid regions, hyper-arid deserts, and semi-deserts, but the genetic basis of their adaptation to different deserts remains unknown.

RESULTS

Here, we sequenced the whole genome of 83 hairy-footed jerboas from distinct desert zones in China to assess how they responded under contrasting conditions. Population genomics analyses reveal the existence of three species in hairy-footed jerboas distributed in China: Dipus deasyi, Dipus sagitta, and Dipus sowerbyi. Analyses of selection between high-altitude desert (elevation ≥ 3000m) and low-altitude desert (< 500m) populations identified two strongly selected genes, ATR and HIF1AN, associated with intense UV radiation and hypoxia in high-altitude environments. A number of candidate genes involved in energy and water homeostasis were detected in the comparative genomic analyses of hyper-arid desert (average annual precipitation < 70mm) and arid desert (< 200mm) populations versus semi-desert (> 360mm) populations. Hyper-arid desert animals also exhibited stronger adaptive selection in energy homeostasis, suggesting water and resource scarcity may be the main drivers of desert adaptation in hairy-footed jerboas.

CONCLUSIONS

Our study challenges the view of deserts as homogeneous environments and shows that distinct genomic adaptations can be found among desert animals depending on their habitats.

Disparate genetic divergence patterns in three corals across a pan-Pacific environmental gradient highlight species-specific adaptation

Tropical coral reefs are among the most affected ecosystems by climate change and face increasing loss in the coming decades. Effective conservation strategies that maximize ecosystem resilience must be informed by the accurate characterization of extant genetic diversity and population structure together with an understanding of the adaptive potential of keystone species. Here we analyzed samples from the Tara Pacific Expedition (2016-2018) that completed an 18,000 km longitudinal transect of the Pacific Ocean sampling three widespread corals-Pocillopora meandrina, Porites lobata, and Millepora cf. platyphylla-across 33 sites from 11 islands. Using deep metagenomic sequencing of 269 colonies in conjunction with morphological analyses and climate variability data, we can show that despite a targeted sampling the transect encompasses multiple cryptic species. These species exhibit disparate biogeographic patterns and, most importantly, distinct evolutionary patterns in identical environmental regimes. Our findings demonstrate on a basin scale that evolutionary trajectories are species-specific and can only in part be predicted from the environment. This highlights that conservation strategies must integrate multi-species investigations to discern the distinct genomic footprints shaped by selection as well as the genetic potential for adaptive change.

Integrative taxonomy in a rapid speciation group associated with mating system transition: A case study in the Primula cicutariifolia complex

Accurate species delimitation is the key to biodiversity conservation and is fundamental to most branches of biology. However, species delimitation remains challenging in those evolutionary radiations associated with mating system transition from outcrossing to self-fertilization, which have frequently occurred in angiosperms and are usually accompanied by rapid speciation. Here, using the Primula cicutariifolia complex as a case, we integrated molecular, morphological and reproductive isolation evidence to test and verify whether its outcrossing (distylous) and selfing (homostylous) populations have developed into independent evolutionary lineages. Phylogenetic trees based on whole plastomes and SNPs of the nuclear genome both indicated that the distylous and homostylous populations grouped into two different clades. Multispecies coalescent, gene flow and genetic structure analyses all supported such two clades as two different genetic entities. In morphology, as expected changes in selfing syndrome, homostylous populations have significantly fewer umbel layers and smaller flower and leaf sizes compared to distylous populations, and the variation range of some floral traits, such as corolla diameter and umbel layers, show obvious discontinuity. Furthermore, hand-pollinated hybridization between the two clades produced almost no seeds, indicating that well post-pollination reproductive isolation has been established between them. Therefore, the distylous and homostylous populations in this studied complex are two independent evolutionary lineages, and thus these distylous populations should be treated as a distinct species, here named Primula qiandaoensis W. Zhang & J.W. Shao sp. nov.. Our empirical study of the P. cicutariifolia complex highlights the importance of applying multiple lines of evidence, in particular genomic data, to delimit species in pervasive evolutionary plant radiations associated with mating system transition.

Ocean crossers: a tale of disjunctions and speciation in the dwarf-fruticose Lichina (lichenized Ascomycota)

Lichens thrive in rocky coastal areas in temperate and cold regions of both hemispheres. Species of the genus Lichina, which form characteristic black fruiting thalli associated with cyanobacteria, often create distinguishable bands in the intertidal and supralittoral zones. The present study uses a comprehensive specimen dataset and four gene loci to (1) delineate and discuss species boundaries in this genus, (2) assess evolutionary relationships among species, and (3) infer the most likely causes of their current geographic distribution in the Northern and Southern hemispheres. A dated phylogeny describes the time frame in which extant disjunctions of species and populations were established. The results showed that the genus is integrated by four species, with Lichina pygmaea, L. confinis and the newly described L. canariensis from rocky seashores in the Canary Islands, occurring in the Northern Hemisphere, whereas L. intermedia is restricted to the Southern Hemisphere. Lichina intermedia hosted a much higher intraspecific genetic diversity than the other species, with subclades interpreted as species-level lineages by the different species delimitation approaches. However, a conservative taxonomic approach was adopted. This species showed a striking disjunct distribution between Australasia and southern South America. The timing for the observed interspecific and intraspecific divergences and population disjunctions postdated continental plate movements, suggesting that long-distance dispersal across body waters in the two hemispheres played a major role in shaping the current species distributions. Such ocean crossings were, as in L. canariensis, followed by speciation. New substitution rates for the nrITS of the genus Lichina were inferred using a tree spanning the major Ascomycota lineages calibrated using fossils. In conclusion, this work lays the foundation for a better understanding of the evolution through time and space of maritime lichens.

From genomics to integrative species delimitation? The case study of the Indo-Pacific Pocillopora corals

With the advent of genomics, sequencing thousands of loci from hundreds of individuals now appears feasible at reasonable costs, allowing complex phylogenies to be resolved. This is particularly relevant for cnidarians, for which insufficient data is available due to the small number of currently available markers and obscures species boundaries. Difficulties in inferring gene trees and morphological incongruences further blur the study and conservation of these organisms. Yet, can genomics alone be used to delimit species? Here, focusing on the coral genus Pocillopora, whose colonies play key roles in Indo-Pacific reef ecosystems but have challenged taxonomists for decades, we explored and discussed the usefulness of multiple criteria (genetics, morphology, biogeography and symbiosis ecology) to delimit species of this genus. Phylogenetic inferences, clustering approaches and species delimitation methods based on genome-wide single-nucleotide polymorphisms (SNP) were first used to resolve Pocillopora phylogeny and propose genomic species hypotheses from 356 colonies sampled across the Indo-Pacific (western Indian Ocean, tropical southwestern Pacific and south-east Polynesia). These species hypotheses were then compared to other lines of evidence based on genetic, morphology, biogeography and symbiont associations. Out of 21 species hypotheses delimited by genomics, 13 were strongly supported by all approaches, while six could represent either undescribed species or nominal species that have been synonymised incorrectly. Altogether, our results support (1) the obsolescence of macromorphology (i.e., overall colony and branches shape) but the relevance of micromorphology (i.e., corallite structures) to refine Pocillopora species boundaries, (2) the relevance of the mtORF (coupled with other markers in some cases) as a diagnostic marker of most species, (3) the requirement of molecular identification when species identity of colonies is absolutely necessary to interpret results, as morphology can blur species identification in the field, and (4) the need for a taxonomic revision of the genus Pocillopora. These results give new insights into the usefulness of multiple criteria for resolving Pocillopora, and more widely, scleractinian species boundaries, and will ultimately contribute to the taxonomic revision of this genus and the conservation of its species.

A New Assessment of Robust Capuchin Monkey (Sapajus) Evolutionary History Using Genome-Wide SNP Marker Data and a Bayesian Approach to Species Delimitation

Robust capuchin monkeys, Sapajus genus, are among the most phenotypically diverse and widespread groups of primates in South America, with one of the most confusing and often shifting taxonomies. We used a ddRADseq approach to generate genome-wide SNP markers for 171 individuals from all putative extant species of Sapajus to access their evolutionary history. Using maximum likelihood, multispecies coalescent phylogenetic inference, and a Bayes Factor method to test for alternative hypotheses of species delimitation, we inferred the phylogenetic history of the Sapajus radiation, evaluating the number of discrete species supported. Our results support the recognition of three species from the Atlantic Forest south of the São Francisco River, with these species being the first splits in the robust capuchin radiation. Our results were congruent in recovering the Pantanal and Amazonian Sapajus as structured into three monophyletic clades, though new morphological assessments are necessary, as the Amazonian clades do not agree with previous morphology-based taxonomic distributions. Phylogenetic reconstructions for Sapajus occurring in the Cerrado, Caatinga, and northeastern Atlantic Forest were less congruent with morphology-based phylogenetic reconstructions, as the bearded capuchin was recovered as a paraphyletic clade, with samples from the Caatinga biome being either a monophyletic clade or nested with the blond capuchin monkey.

Genome-wide data reveal cryptic diversity and hybridization in a group of tree ferns

Discovery of cryptic diversity is essential to understanding both the process of speciation and the conservation of species. Determining species boundaries in fern lineages represents a major challenge due to lack of morphologically diagnostic characters and frequent hybridization. Genomic data has substantially enhanced our understanding of the speciation process, increased the resolution of species delimitation studies, and led to the discovery of cryptic diversity. Here, we employed restriction-site-associated DNA sequencing (RAD-seq) and integrated phylogenomic and population genomic analyses to investigate phylogenetic relationships and evolutionary history of 16 tree ferns with marginate scales (Cyatheaceae) from China and Vietnam. We conducted multiple species delimitation analyses using the multispecies coalescent (MSC) model and novel approaches based on genealogical divergence index (gdi) and isolation by distance (IBD). In addition, we inferred species trees using concatenation and several coalescent-based methods, and assessed hybridization patterns and rate of gene flow across the phylogeny. We obtained highly supported and generally congruent phylogenies inferred from concatenated and summary-coalescent methods, and the monophyly of all currently recognized species were strongly supported. Our results revealed substantial evidence of cryptic diversity in three widely distributed Gymnosphaera species, each of which was composite of two highly structure lineages that may correspond to cryptic species. We found that hybridization was fairly common between not only closely related species, but also distantly related species. Collectively, it appears that scaly tree ferns may contain cryptic diversity and hybridization has played an important role throughout the evolutionary history of this group.

Taxonomic review of Saguinus mystax (Spix, 1823) (Primates, Callitrichidae), and description of a new species

Although the Amazon has the greatest diversity of primates, there are still taxonomic uncertainties for many taxa, such as the species of the Saguinus mystax group. The most geographically broadly distributed and phenotypically diverse species in this group is S. mystax, and its phenotypic diversity has been recognized as three subspecies-S. mystax mystax, S. mystax pileatus and S. mystax pluto-with non-overlapping geographic distributions. In this sense, we carried out an extensive field survey in their distribution areas and used a framework of taxonomic hypothesis testing of genomic data combined with an integrative taxonomic decision-making framework to carry out a taxonomic revision of S. mystax. Our tests supported the existence of three lineages/species. The first species corresponds to Saguinus mystax mystax from the left bank of the Juruá River, which was raised to the species level, and we also discovered and described animals from the Juruá-Tefé interfluve previously attributed to S. mystax mystax as a new species. The subspecies S. m. pileatus and S. m. pluto are recognized as a single species, under a new nomenclatural combination. However, given their phenotypic distinction and allopatric distribution, they potentially are a manifestation of an early stage of speciation, and therefore we maintain their subspecific designations.

Contrasting Phylogeographic Patterns of Mitochondrial and Genome-Wide Variation in the Groundwater Amphipod Crangonyx islandicus That Survived the Ice Age in Iceland

The analysis of phylogeographic patterns has often been based on mitochondrial DNA variation, but recent analyses dealing with nuclear DNA have in some instances revealed mito-nuclear discordances and complex evolutionary histories. These enigmatic scenarios, which may involve stochastic lineage sorting, ancestral hybridization, past dispersal and secondary contacts, are increasingly scrutinized with a new generation of genomic tools such as RADseq, which also poses additional analytical challenges. Here, we revisited the previously inconclusive phylogeographic history, showing the mito-nuclear discordance of an endemic groundwater amphipod from Iceland, Crangonyx islandicus, which is the only metazoan known to have survived the Pleistocene beneath the glaciers. Previous studies based on three DNA markers documented a mitochondrial scenario with the main divergence occurring between populations in northern Iceland and an ITS scenario with the main divergence between the south and north. We used double digest restriction-site-associated DNA sequencing (ddRADseq) to clarify this mito-nuclear discordance by applying several statistical methods while estimating the sensitivity to different analytical approaches (data-type, differentiation indices and base call uncertainty). A majority of nuclear markers and methods support the ITS divergence. Nevertheless, a more complex scenario emerges, possibly involving introgression led by male-biased dispersal among northern locations or mitochondrial capture, which may have been further strengthened by natural selection.

Cryptic species diversity in a widespread neotropical tree genus: The case of Cedrela odorata

PREMISE

Reconciling the use of taxonomy to partition morphological variation and describe genetic divergence within and among closely related species is a persistent challenge in phylogenetics. We reconstructed phylogenetic relationships among Cedrela odorata (Meliaceae) and five closely allied species to test the genetic basis for the current model of species delimitation in this economically valuable and threatened genus.

METHODS

We prepared a nuclear species tree with the program SNPhylo and 16,000 single-nucleotide polymorphisms from 168 Cedrela specimens. Based on clades present and ancestral patterns ADMIXTURE, we designed nine species delimitation models and compared each model to current taxonomy with Bayes factor delimitation. Timing of major lineage divergences was estimated with the program SNAPP.

RESULTS

The resulting analysis revealed that modern C. odorata evolved from two genetically distinct ancestral sources. All species delimitation models tested better fit the data than the model representing current taxonomic delimitation. Models with the greatest marginal likelihoods separated Mesoamerican C. odorata and South American C. odorata into two species and lumped C. angustifolia and C. montana as a single species. We estimated that Cedrela diversified in South America within the last 19 million years following one or more dispersal events from Mesoamerican lineages.

CONCLUSIONS

Our analyses show that the present taxonomic understanding within the genus obscures divergent lineages in C. odorata due in part to morphological differentiation and taxonomic distinctions that are not predictably associated with genetic divergence. A more accurate application of taxonomy to C. odorata and related species may aid in its conservation, management, and restoration efforts.

Lineage and role in integrative taxonomy of a heterotrophic orchid complex

Lineage-based species definitions applying coalescent approaches to species delimitation have become increasingly popular. Yet, the application of these methods and the recognition of lineage-only definitions have recently been questioned. Species delimitation criteria that explicitly consider both lineages and evidence for ecological role shifts provide an opportunity to incorporate ecologically meaningful data from multiple sources in studies of species boundaries. Here, such criteria were applied to a problematic group of mycoheterotrophic orchids, the Corallorhiza striata complex, analysing genomic, morphological, phenological, reproductive-mode, niche, and fungal host data. A recently developed method for generating genomic polymorphism data-ISSRseq-demonstrates evidence for four distinct lineages, including a previously unidentified lineage in the Coast Ranges and Cascades of California and Oregon, USA. There is divergence in morphology, phenology, reproductive mode, and fungal associates among the four lineages. Integrative analyses, conducted in population assignment and redundancy analysis frameworks, provide evidence of distinct genomic lineages and a similar pattern of divergence in the extended data, albeit with weaker signal. However, none of the extended data sets fully satisfy the condition of a significant role shift, which requires evidence of fixed differences. The four lineages identified in the current study are recognized at the level of variety, short of comprising different species. This study represents the most comprehensive application of lineage + role to date and illustrates the advantages of such an approach.

Evolutionary conservation genomics reveals recent speciation and local adaptation in threatened takins

Incorrect species delimitation will lead to inappropriate conservation decisions, especially for threatened species. The takin (Budorcas taxicolor) is a large artiodactyl endemic to the Himalayan–Hengduan–Qinling Mountains and is well known for its threatened status and peculiar appearance. However, the speciation, intraspecies taxonomy, evolutionary history, and adaptive evolution of this species still remain unclear, which greatly hampers its scientific conservation. Here, we de novo assembled a high-quality chromosome-level genome of takin and resequenced the genomes of 75 wild takins. Phylogenomics revealed that takin was positioned at the root of Caprinae. Population genomics based on the autosome, X chromosome, and Y chromosome SNPs and mitochondrial genomes consistently revealed the existence of two phylogenetic species and recent speciation in takins: the Himalayan takin (B. taxicolor) and the Chinese takin (B. tibetana), with the support of morphological evidence. Two genetically divergent subspecies were identified in both takin species, rejecting three previously proposed taxonomical viewpoints. Furthermore, their distribution boundaries were determined, suggesting that large rivers play important roles in shaping the genetic partition. Compared with the other subspecies, the Qinling subspecies presented the lowest genomic diversity, higher linkage disequilibrium, inbreeding, and genetic load, thus is in urgent need of genetic management and protection. Moreover, coat color gene (PMEL) variation may be responsible for the adaptive coat color difference between the two species following Gloger’s rule. Our findings provide novel insights into the recent speciation, local adaptation, scientific conservation of takins, and biogeography of the Himalaya–Hengduan biodiversity hotspot.

Lack of Statistical Rigor in DNA Barcoding Likely Invalidates the Presence of a True Species' Barcode Gap

DNA barcoding has been largely successful in satisfactorily exposing levels of standing genetic diversity for a wide range of taxonomic groups through the employment of only one or a few universal gene markers. However, sufficient coverage of geographically-broad intra-specific haplotype variation within genomic databases like the Barcode of Life Data Systems (BOLD) and GenBank remains relatively sparse. As reference sequence libraries continue to grow exponentially in size, there is now the need to identify novel ways of meaningfully analyzing vast amounts of available DNA barcode data. This is an important issue to address promptly for the routine tasks of specimen identification and species discovery, which have seen broad adoption in areas as diverse as regulatory forensics and resource conservation. Here, it is demonstrated that the interpretation of DNA barcoding data is lacking in statistical rigor. To highlight this, focus is set specifically on one key concept that has become a household name in the field: the DNA barcode gap. Arguments outlined herein specifically center on DNA barcoding in animal taxa and stem from three angles: (1) the improper allocation of specimen sampling effort necessary to capture adequate levels of within-species genetic variation, (2) failing to properly visualize intra-specific and interspecific genetic distances, and (3) the inconsistent, inappropriate use, or absence of statistical inferential procedures in DNA barcoding gap analyses. Furthermore, simple statistical solutions are outlined which can greatly propel the use of DNA barcoding as a tool to irrefutably match unknowns to knowns on the basis of the barcoding gap with a high degree of confidence. Proposed methods examined herein are illustrated through application to DNA barcode sequence data from Canadian Pacific fish species as a case study.

Genome and cuticular hydrocarbon-based species delimitation shed light on potential drivers of speciation in a Neotropical ant species complex

Geographic separation that leads to the evolution of reproductive isolation between populations generally is considered the most common form of speciation. However, speciation may also occur in the absence of geographic barriers due to phenotypic and genotypic factors such as chemical cue divergence, mating signal divergence, and mitonuclear conflict. Here, we performed an integrative study based on two genome-wide techniques (3RAD and ultraconserved elements) coupled with cuticular hydrocarbon (CHC) and mitochondrial (mt) DNA sequence data, to assess the species limits within the Ectatomma ruidum species complex, a widespread and conspicuous group of Neotropical ants for which heteroplasmy (i.e., presence of multiple mtDNA variants in an individual) has been recently discovered in some populations from southeast Mexico. Our analyses indicate the existence of at least five distinct species in this complex: two widely distributed across the Neotropics, and three that are restricted to southeast Mexico and that apparently have high levels of heteroplasmy. We found that species boundaries in the complex did not coincide with geographic barriers. We therefore consider possible roles of alternative drivers that may have promoted the observed patterns of speciation, including mitonuclear incompatibility, CHC differentiation, and colony structure. Our study highlights the importance of simultaneously assessing different sources of evidence to disentangle the species limits of taxa with complicated evolutionary histories.

Genomics and ecological modelling clarify species integrity in a confusing group of butterflies

Recent advances in both genomics and ecological modelling present new, multidisciplinary opportunities for resolving species boundaries and understanding the mechanisms that maintain their integrity in regions of contact. Here, we use a combination of high-throughput DNA sequencing and ecological niche modelling to resolve species boundaries and niche divergence within the Speyeria atlantis-hesperis (Lepidoptera: Nymphalidae) complex, a confusing group of North American butterflies. This complex is notorious for its muddled species delimitations, morphological ambiguity, and extensive mito-nuclear discordance. Our admixture and multispecies coalescent-based analyses of single nucleotide polymorphisms identified substantial divergences between S. atlantis and S. hesperis in areas of contact, as well as between distinct northern and southern lineages within S. hesperis. Our results also provide evidence of past introgression relating to another species, S. zerene, which previous work has shown to be more distantly related to the S. atlantis-hesperis complex. We then used ecological models to predict habitat suitability for each of the three recovered genomic lineages in the S. atlantis-hesperis complex and assess their pairwise niche divergence. These analyses resolved that these three lineages are significantly diverged in their respective niches and are not separated by discontinuities in suitable habitat that might present barriers to gene flow. We therefore infer that ecologically-mediated selection resulting in disparate habitat associations is a principal mechanism reinforcing their genomic integrity. Overall, our results unambiguously support significant evolutionary and ecological divergence between the northern and southern lineages of S. hesperis, sufficient to recognize the southern evolutionary lineage as a distinct species, called S. nausicaa based on taxonomic priority.

Using natural history to guide supervised machine learning for cryptic species delimitation with genetic data

The diversity of biological and ecological characteristics of organisms, and the underlying genetic patterns and processes of speciation, makes the development of universally applicable genetic species delimitation methods challenging. Many approaches, like those incorporating the multispecies coalescent, sometimes delimit populations and overestimate species numbers. This issue is exacerbated in taxa with inherently high population structure due to low dispersal ability, and in cryptic species resulting from nonecological speciation. These taxa present a conundrum when delimiting species: analyses rely heavily, if not entirely, on genetic data which over split species, while other lines of evidence lump. We showcase this conundrum in the harvester Theromaster brunneus, a low dispersal taxon with a wide geographic distribution and high potential for cryptic species. Integrating morphology, mitochondrial, and sub-genomic (double-digest RADSeq and ultraconserved elements) data, we find high discordance across analyses and data types in the number of inferred species, with further evidence that multispecies coalescent approaches over split. We demonstrate the power of a supervised machine learning approach in effectively delimiting cryptic species by creating a "custom" training data set derived from a well-studied lineage with similar biological characteristics as Theromaster. This novel approach uses known taxa with particular biological characteristics to inform unknown taxa with similar characteristics, using modern computational tools ideally suited for species delimitation. The approach also considers the natural history of organisms to make more biologically informed species delimitation decisions, and in principle is broadly applicable for taxa across the tree of life.

From Gondwana to the Yellow Sea, evolutionary diversifications of true toads Bufo sp. in the Eastern Palearctic and a revisit of species boundaries for Asian lineages

Taxa with vast distribution ranges often display unresolved phylogeographic structures and unclear taxonomic boundaries resulting in hidden diversity. This hypothesis-driven study reveals the evolutionary history of Bufonidae, covering the phylogeographic patterns found in Holarctic bufonids from the West Gondwana to the phylogenetic taxonomy of Asiatic true toads in the Eastern Palearctic. We used an integrative approach relying on fossilized birth-death calibrations, population dynamics, gene-flow, species distribution, and species delimitation modeling to resolve the biogeography of the clade and highlight cryptic lineages. We verified the near-simultaneous Miocene radiations within Western and Eastern Palearctic Bufo, c. 14.49-10.00 Mya, temporally matching with the maximum dust outflows in Central Asian deserts. Contrary to earlier studies, we demonstrated that the combined impacts of long dispersal and ice-age refugia equally contributed to the current genetic structure of Bufo in East Asia. Our findings reveal a climate-driven adaptation in septentrional Eastern Asian Bufo, explaining its range shifts toward northern latitudes. We resolve species boundaries within the Eastern Palearctic Bufo, and redefine the taxonomic and conservation units of the northeastern species: B. sachalinensis and its subspecies.

Variation in size and shape sexual dimorphism in the Sceloporus scalaris species group (Squamata: Phrynosomatidae) from the Transvolcanic Belt of Mexico

We attempted to identify the factors influencing size and shape dimorphism between sexes, as well as among populations and species in the Sceloporus scalaris group (Sceloporus aeneus, S. scalaris, S. bicanthalis and S. subniger). Our analysis focused on five morphological characteristics: snout–vent length, head length, head width, forearm length and tibia length. The effect of environmental variables (precipitation and temperature) on these variables was also tested. We found differences in morphological traits between sexes, and among populations of the same species. The oviparous species (S. aeneus and S. scalaris) were larger in overall body size than the viviparous species (S. bicanthalis and S. subniger). Differences in overall body size among populations were recorded only in S. aeneus and S. scalaris. Male-biased sexual size dimorphism occurred in oviparous but not viviparous lizards (except for one population of S. bicanthalis). An absence of sexual size dimorphism was also recorded in S. subniger and some populations of the remaining species. Two different shape patterns were found; the first was female-biased with larger relative body length in almost all populations, which could be explained by fecundity, and the second was male-biased with relatively larger head and limbs in a few populations, which may be explained by sexual selection. The patterns of sexual size and shape dimorphism show that environment, rather than phylogeny, may be determining the extent of sexual dimorphism. These types of studies show the importance of an integrated evaluation of interpopulation and interspecies variation to determine the factors that generate sexual dimorphism.

OUP accepted manuscript

As genomic sequence data become increasingly available, inferring the phylogeny of the species as that of concatenated genomic data can be enticing. However, this approach makes for a biased estimator of branch lengths and substitution rates and an inconsistent estimator of tree topology. Bayesian multispecies coalescent (MSC) methods address these issues. This is achieved by constraining a set of gene trees within a species tree and jointly inferring both under a Bayesian framework. However, this approach comes at the cost of increased computational demand. Here, we introduce StarBeast3—a software package for efficient Bayesian inference under the MSC model via Markov chain Monte Carlo. We gain efficiency by introducing cutting-edge proposal kernels and adaptive operators, and StarBeast3 is particularly efficient when a relaxed clock model is applied. Furthermore, gene-tree inference is parallelized, allowing the software to scale with the size of the problem. We validated our software and benchmarked its performance using three real and two synthetic data sets. Our results indicate that StarBeast3 is up to one-and-a-half orders of magnitude faster than StarBeast2, and therefore more than two orders faster than *BEAST, depending on the data set and on the parameter, and can achieve convergence on large data sets with hundreds of genes. StarBeast3 is open-source and is easy to set up with a friendly graphical user interface. [Adaptive; Bayesian inference; BEAST 2; effective population sizes; high performance; multispecies coalescent; parallelization; phylogenetics.]

Molecular species delimitation of shrub frogs of the genus Pseudophilautus (Anura, Rhacophoridae)

Sri Lanka is an amphibian hotspot of global significance. Its anuran fauna is dominated by the shrub frogs of the genus Pseudophilautus. Except for one small clade of four species in Peninsular India, these cool-wet adapted frogs, numbering some 59 extant species, are distributed mainly across the montane and lowland rain forests of the island. With species described primarily by morphological means, the diversification has never yet been subjected to a molecular species delimitation analysis, a procedure now routinely applied in taxonomy. Here we test the species boundaries of Pseudophilautus in the context of the phylogenetic species concept (PSC). We use all the putative species for which credible molecular data are available (nDNA-Rag-1; mt-DNA- 12S rRNA, 16S rRNA) to build a well resolved phylogeny, which is subjected to species delimitation analyses. The ABGD, bPTP, mPTP and bGMYC species delimitation methods applied to the 16S rRNA frog barcoding gene (for all species), 12S rRNA and Rag-1 nDNA grouped P. procax and P. abundus; P. hallidayi and P. fergusonianus; P. reticulatus and P. pappilosus; P. pleurotaenia and P. hoipolloi; P. hoffmani and P. asankai; P. silvaticus and P. limbus; P. dilmah and P. hankeni; P. fulvus and P. silus.. Surprisingly, all analyses recovered 14 unidentified potential new species as well. The geophylogeny affirms a distribution across the island's aseasonal 'wet zone' and its three principal hill ranges, suggestive of allopatric speciation playing a dominant role, especially between mountain masses. Among the species that are merged by the delimitation analyses, a pattern leading towards a model of parapatric speciation emerges-ongoing speciation in the presence of gene flow. This delimitation analysis reinforces the species hypotheses, paving the way to a reasonable understanding of Sri Lankan Pseudophilautus, enabling both deeper analyses and conservation efforts of this remarkable diversification. http://zoobank.org/urn:lsid:zoobank.org:pub:DA869B6B-870A-4ED3-BF5D-5AA3F69DDD27.

Solving the Coral Species Delimitation Conundrum

Distinguishing coral species is not only crucial for physiological, ecological and evolutionary studies, but also to enable effective management of threatened reef ecosystems. However, traditional hypotheses that delineate coral species based on morphological traits from the coral skeleton are frequently at odds with tree-based molecular approaches. Additionally, a dearth of species-level molecular markers has made species delimitation particularly challenging in species-rich coral genera, leading to the widespread assumption that inter-specific hybridization might be responsible for this apparent conundrum. Here, we used three lines of evidence - morphology, breeding trials and molecular approaches - to identify species boundaries in a group of ecologically important tabular Acropora corals. In contrast to previous studies, our morphological analyses yielded groups that were congruent with experimental crosses as well as with coalescent-based and allele sharing-based multilocus approaches to species delimitation. Our results suggest that species of the genus Acropora are reproductively isolated and independently evolving units that can be distinguished morphologically. These findings not only pave the way for a taxonomic revision of coral species, but also outline an approach that can provide a solid basis to address species delimitation and provide conservation support to a wide variety of keystone organisms.

A species of the genus Panophrys (Anura, Megophryidae) from southeastern Guizhou Province, China

Herein, we describe Panophryscongjiangensissp. nov. obtained from the Yueliangshan Nature Reserve, Congjiang County, Guizhou Province, China. Phylogenetic analyses based on the mitochondrial genes 16S rRNA and COI indicated that this new species represented an independent lineage, closely related to P.leishanensis. The uncorrected genetic distances between the new species and its closest congener, P.leishanensis, were 3.0% for 16S rRNA and 8.4% for COI. The new species is distinguished from its congeners by a combination of the following morphological characteristics (1) medium body size (SVL 28.6–33.4 mm in males and 38.4–40.2 mm in females); (2) a small horn-like tubercle at the edge of each upper eyelid; (3) the tympanum distinctly visible (TD/ED ratio 0.47–0.66); (4) vomerine teeth absent; (5) the tongue not notched behind; (6) a narrow and unobvious lateral fringe on toes; (7) relative finger lengths II < I < V < III; (8) rudimentary webs on toes; (9) hindlimbs slender, heels overlapping when thighs are positioned at right angles to the body; (10) two metacarpal tubercles on the palm, with the inner metatarsal tubercle long and oval-shaped; (11) the tibiotarsal articulation reaching the nostril when the leg is adpressed and stretched forward; (12) dorsal skin rough with numerous orange–red granules, ventral surface smooth; (13) a single internal subgular vocal sac present in males; and (14) in breeding males, weak gray-black nuptial pads with black nuptial spines present on the dorsal surface of the bases of the first and second fingers. To date, the new species is only known from the type locality.

The biogeography of Dromiciops in southern South America: Middle Miocene transgressions, speciation and associations with Nothofagus

The current distribution of the flora and fauna of southern South America is the result of drastic geological events that occurred during the last 20 million years, including marine transgressions, glaciations and active vulcanism. All these have been associated with fragmentation, isolation and subsequent expansion of the biota, south of 35°S, such as the temperate rainforest. This forest is mostly dominated by Nothofagus trees and is the habitat of the relict marsupial monito del monte, genus Dromiciops, sole survivor of the order Microbiotheria. Preliminary analyses using mtDNA proposed the existence of three main Dromiciops lineages, distributed latitudinally, whose divergence was initially attributed to recent Pleistocene glaciations. Using fossil-calibrated dating on nuclear and mitochondrial genes, here we reevaluate this hypothesis and report an older (Miocene) biogeographic history for the genus. We performed phylogenetic reconstructions using sequences from two mitochondrial DNA and four nuclear DNA genes in 159 specimens from 31 sites across Chile and Argentina. Our phylogenetic analysis resolved three main clades with discrete geographic distributions. The oldest and most differentiated clade corresponds to that of the northern distribution (35.2°S to 39.3°S), which should be considered a distinct species (D. bozinovici, sensu D'Elía et al. 2016). According to our estimations, this species shared a common ancestor with D. gliroides (southern clades) about ~13 million years ago. Divergence time estimates for the southern clades (39.6°S to 42.0°S) ranged from 9.57 to 6.5 Mya. A strong genetic structure was also detected within and between clades. Demographic analyses suggest population size stability for the northern clade (D. bozinovici), and recent demographic expansions for the central and southern clades. All together, our results suggest that the diversification of Dromiciops were initiated by the Middle Miocene transgression (MMT), the massive marine flooding that covered several lowlands of the western face of Los Andes between 37 and 48°S. The MMT resulted from an increase in global sea levels at the Miocene climatic optimum, which shaped the biogeographic origin of several species, including Nothofagus forests, the habitat of Dromiciops.

Introgression and Species Delimitation in the Longear Sunfish Lepomis megalotis (Teleostei: Percomorpha: Centrarchidae)

Introgression and hybridization are major impediments to genomic-based species delimitation because many implementations of the multispecies coalescent framework assume no gene flow among species. The sunfish genus Lepomis, one of the world's most popular groups of freshwater sport fish, has a complicated taxonomic history. The results of ddRAD phylogenomic analyses do not provide support for the current taxonomy that recognizes two species, L. megalotis and L. peltastes, in the L. megalotis complex. Instead, evidence from phylogenomics and phenotype warrants recognizing six relatively ancient evolutionary lineages in the complex. The introgressed and hybridizing populations in the L. megalotis complex are localized and appear to be the result of secondary contact or rare hybridization events between non-sister species. Segregating admixed populations from our multispecies coalescent analyses identifies six species with moderate to high genealogical divergence, whereas including admixed populations drives all but one lineage below the species threshold of genealogical divergence. Segregation of admixed individuals also helps reveal phenotypic distinctiveness among the six species in morphological traits used by ichthyologists to discover and delimit species over the last two centuries. Our protocols allow for the identification and accommodation of hybridization and introgression in species delimitation. Genomic-based species delimitation validated with multiple lines of evidence provides a path towards the discovery of new biodiversity and resolving long-standing taxonomic problems.

Sareomycetes: more diverse than meets the eye

Since its resurrection, the resinicolous discomycete genus Sarea has been accepted as containing two species, one with black apothecia and pycnidia, and one with orange. We investigate this hypothesis using three ribosomal (nuITS, nuLSU, mtSSU) regions from and morphological examination of 70 specimens collected primarily in Europe and North America. The results of our analyses support separation of the traditional Sarea difformis s.lat. and Sarea resinae s.lat. into two distinct genera, Sarea and Zythia. Sarea as circumscribed is shown to conservatively comprise three phylospecies, with one corresponding to Sarea difformis s.str. and two, morphologically indistinguishable, corresponding to the newly combined Sarea coeloplata. Zythia is provisionally maintained as monotypic, containing only a genetically and morphologically variable Z. resinae. The new genus Atrozythia is erected for the new species A. klamathica. Arthrographis lignicola is placed in this genus on molecular grounds, expanding the concept of Sareomycetes by inclusion of a previously unknown type of asexual morph. Dating analyses using additional marker regions indicate the emergence of the Sareomycetes was roughly concurrent with the diversification of the genus Pinus, suggesting that this group of fungi emerged to exploit the newly-available resinous ecological niche supplied by Pinus or another, extinct group of conifers. Our phylogeographic studies also permitted us to study the introductions of these fungi to areas where they are not native, including Antarctica, Cape Verde, and New Zealand and are consistent with historical hypotheses of introduction.

Single-locus species delimitation and ecological niche modelling provide insights into the evolution, historical distribution and taxonomy of the Pacific chorus frogs

The Pacific chorus frogs are a complex of three wide-ranging species (i.e. Hyliola hypochondriaca, Hyliola regilla, Hyliola sierra) whose current taxonomy remains unresolved. We conducted species delimitation analyses of these taxa using fragments of the cytochrome b and 12S–16S mtDNA genes to assess the species diversity. Importantly, we included samples from new locations throughout the range to better understand species distributions and identify potential contact zones among clades. Our analyses revealed three slightly parapatric but distinct species-level clades. Molecular dating revealed that these species began diverging in the Pleistocene c. 1.4 Mya with H. hypochondriaca and H. sierra diverging more recently c. 0.8 Mya. We found that populations from western Montana and Idaho originated recently from populations to the southwest that belong to H. sierra, rather than from H. regilla populations directly to the west. Population sizes of each species expanded c. 130–80 Kya with H. hypochondriaca exhibiting a more pronounced expansion beginning c. 100 Kya than the more gradual expansion of the other two species. The climatic niche models suggest that distributions of the three species were similar during the last interglacial (LIG) as they are today. During the Last Glacial Maximum (LGM), H. hypochondriaca and H. sierra occupied a larger range than they do today whereas H. regilla occupied a smaller refugium, shifted south from the current distribution. This study highlights the continued effectiveness of utilizing single-locus data sets for species delimitation and biogeographic analyses.