How to fail at species delimitation

Machine learning approaches delimit cryptic taxa in a previously intractable species complex

Cryptic species are not diagnosable via morphological criteria, but can be detected through analysis of DNA sequences. A number of methods have been developed for identifying species based on genetic data; however, these methods are prone to over-splitting taxa with extreme population structure, such as dispersal-limited organisms. Machine learning methodologies have the potential to overcome this challenge. Here, we apply such approaches, using a large dataset generated through hybrid target enrichment of ultraconserved elements (UCEs). Our study taxon is the Aoraki denticulata species complex, a lineage of extremely low-dispersal arachnids endemic to the South Island of Aotearoa New Zealand. This group of mite harvesters has been the subject of previous species delimitation studies using smaller datasets generated through Sanger sequencing and analytical approaches that rely on multispecies coalescent models and barcoding gap discovery. Those analyses yielded a number of putative cryptic species that seems unrealistic and extreme, based on what we know about species' geographic ranges and genetic diversity in non-cryptic mite harvesters. We find that machine learning approaches, on the other hand, identify cryptic species with geographic ranges that are similar to those seen in other morphologically diagnosable mite harvesters in Aotearoa New Zealand's South Island. We performed both unsupervised and supervised machine learning analyses, the latter with training data drawn either from animals broadly (vagile and non-vagile) or from a custom training dataset from dispersal-limited harvesters. We conclude that applying machine learning approaches to the analysis of UCE-derived genetic data is an effective method for delimiting species in complexes of low-vagility cryptic species, and that the incorporation of training data from biologically relevant analogues can be critically informative.

Integrative species delimitation and five new species of lynx spiders (Araneae, Oxyopidae) in Taiwan

An accurate assessment of species diversity is a cornerstone of biology and conservation. The lynx spiders (Araneae: Oxyopidae) represent one of the most diverse and widespread cursorial spider groups, however their species richness in Asia is highly underestimated. In this study, we revised species diversity with extensive taxon sampling in Taiwan and explored species boundaries based on morphological traits and genetic data using a two-step approach of molecular species delimitation. Firstly, we employed a single COI dataset and applied two genetic distance-based methods: ABGD and ASAP, and two topology-based methods: GMYC and bPTP. Secondly, we further analyzed the lineages that were not consistently delimited, and incorporated H3 to the dataset for a coalescent-based analysis using BPP. A total of eight morphological species were recognized, including five new species, Hamataliwa cordivulva sp. nov., Hamat. leporauris sp. nov., Tapponia auriola sp. nov., T. parva sp. nov. and T. rarobulbus sp. nov., and three newly recorded species, Hamadruas hieroglyphica (Thorell, 1887), Hamat. foveata Tang & Li, 2012 and Peucetia latikae Tikader, 1970. All eight morphological species exhibited reciprocally monophyletic lineages. The results of molecular-based delimitation analyses suggested a variety of species hypotheses that did not fully correspond to the eight morphological species. We found that Hamat. cordivulva sp. nov. and Hamat. foveata showed shallow genetic differentiation in the COI, but they were unequivocally distinguishable according to their genitalia. In contrast, T. parva sp. nov. represented a deep divergent lineage, while differences of genitalia were not detected. This study highlights the need to comprehensively employ multiple evidence and methods to delineate species boundaries and the values of diagnostic morphological characters for taxonomic studies in lynx spiders.

Morphological Species Delimitation in The Western Pond Turtle (Actinemys): Can Machine Learning Methods Aid in Cryptic Species Identification?

As the discovery of cryptic species has increased in frequency, there has been an interest in whether geometric morphometric data can detect fine-scale patterns of variation that can be used to morphologically diagnose such species. We used a combination of geometric morphometric data and an ensemble of five supervised machine learning methods (MLMs) to investigate whether plastron shape can differentiate two putative cryptic turtle species, Actinemys marmorata and Actinemys pallida. Actinemys has been the focus of considerable research due to its biogeographic distribution and conservation status. Despite this work, reliable morphological diagnoses for its two species are still lacking. We validated our approach on two datasets, one consisting of eight morphologically disparate emydid species, the other consisting of two subspecies of Trachemys (T. scripta scripta, T. scripta elegans). The validation tests returned near-perfect classification rates, demonstrating that plastron shape is an effective means for distinguishing taxonomic groups of emydids via MLMs. In contrast, the same methods did not return high classification rates for a set of alternative phylogeographic and morphological binning schemes in Actinemys. All classification hypotheses performed poorly relative to the validation datasets and no single hypothesis was unequivocally supported for Actinemys. Two hypotheses had machine learning performance that was marginally better than our remaining hypotheses. In both cases, those hypotheses favored a two-species split between A. marmorata and A. pallida specimens, lending tentative morphological support to the hypothesis of two Actinemys species. However, the machine learning results also underscore that Actinemys as a whole has lower levels of plastral variation than other turtles within Emydidae, but the reason for this morphological conservatism is unclear.

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.

Metazoa-level USCOs as markers in species delimitation and classification

Metazoa-level universal single-copy orthologs (mzl-USCOs) are universally applicable markers for DNA taxonomy in animals that can replace or supplement single-gene barcodes. Previously, mzl-USCOs from target enrichment data were shown to reliably distinguish species. Here, we tested whether USCOs are an evenly distributed, representative sample of a given metazoan genome and therefore able to cope with past hybridization events and incomplete lineage sorting. This is relevant for coalescent-based species delimitation approaches, which critically depend on the assumption that the investigated loci do not exhibit autocorrelation due to physical linkage. Based on 239 chromosome-level assembled genomes, we confirmed that mzl-USCOs are genetically unlinked for practical purposes and a representative sample of a genome in terms of reciprocal distances between USCOs on a chromosome and of distribution across chromosomes. We tested the suitability of mzl-USCOs extracted from genomes for species delimitation and phylogeny in four case studies: Anopheles mosquitos, Drosophila fruit flies, Heliconius butterflies and Darwin's finches. In almost all instances, USCOs allowed delineating species and yielded phylogenies that corresponded to those generated from whole genome data. Our phylogenetic analyses demonstrate that USCOs may complement single-gene DNA barcodes and provide more accurate taxonomic inferences. Combining USCOs from sources that used different versions of ortholog reference libraries to infer marker orthology may be challenging and, at times, impact taxonomic conclusions. However, we expect this problem to become less severe as the rapidly growing number of reference genomes provides a better representation of the number and diversity of organismal lineages.

The genomic diversity of the Eliurus genus in northern Madagascar with a putative new species

Madagascar exhibits extraordinarily high level of species richness and endemism, while being severely threatened by habitat loss and fragmentation (HL&F). In front of these threats to biodiversity, conservation effort can be directed, for instance, in the documentation of species that are still unknown to science, or in investigating how species respond to HL&F. The tufted-tail rats genus (Eliurus spp.) is the most speciose genus of endemic rodents in Madagascar, with 13 described species, which occupy two major habitat types: dry or humid forests. The large species diversity and association to specific habitat types make the Eliurus genus a suitable model for investigating species adaptation to new environments, as well as response to HL&F (dry vs humid). In the present study, we investigated Eliurus spp. genomic diversity across northern Madagascar, a region covered by both dry and humid fragmented forests. From the mitochondrial DNA (mtDNA) and nuclear genomic (RAD-seq) data of 124 Eliurus individuals sampled in poorly studied forests of northern Madagascar, we identified an undescribed Eliurus taxon (Eliurus sp. nova). We tested the hypothesis of a new Eliurus species using several approaches: i) DNA barcoding; ii) phylogenetic inferences; iii) species delimitation tests based on the Multi-Species Coalescent (MSC) model, iv) genealogical divergence index (gdi); v) an ad-hoc test of isolation-by-distance within versus between sister-taxa, vi) comparisons of %GC content patterns and vii) morphological analyses. All analyses support the recognition of the undescribed lineage as a putative distinct species. In addition, we show that Eliurus myoxinus, a species known from the dry forests of western Madagascar, is, surprisingly, found mostly in humid forests in northern Madagascar. In conclusion, we discuss the implications of such findings in the context of Eliurus species evolution and diversification, and use the distribution of northern Eliurus species as a proxy for reconstructing past changes in forest cover and vegetation type in northern Madagascar.

Delimiting the rare, endangered and actively speciating

Species delimitation is a contentious topic. The genomics revolution initially brought hope that identifying and classifying species would be easier through better methods and more data, but genomics has also brought complexity and controversy to delimitation. One solution can be to collect a larger sample of individuals at a finer geographic scale. But what if taxa are rare and collecting more samples is difficult or detrimental to the organisms at hand? In this issue of Molecular Ecology Resources, Opatova et al. (2023) tackle the ambiguity of species delimitation in rare and endangered trapdoor spiders (genus Cyclocosmia). The authors propose a framework for delimiting species when samples are hard to come by, such as in these rare and cryptic spiders. The authors combine extensive genomic sampling with statistical approaches that consider both the genetic distinctiveness of each population of spiders and how much gene flow occurs between these populations. Their proposed taxonomy balances two opposing signals, structure and gene flow, to count eight lineages of Cyclocosmia, and to point the way for future taxonomic studies of the rare or difficult to obtain.

Unveiling taxonomic diversity in the deep-sea fish genus Notacanthus (Notacanthiformes: Notacanthidae) with description of Notacanthus arrontei n. sp

Notacanthid fishes constitute a common part of benthopelagic deep-sea fish communities on seamounts and continental slopes around the world. However, their highly conserved morphology and the usual lack of information on deep-water organisms make it difficult to appropriately address their biodiversity. A multidisciplinary approach combining morphological data with a DNA-based species delimitation analyses was used to explore the taxonomy of Notacanthus species. For this purpose, morphological and molecular data were obtained from 43 individuals, and the resulting information was combined with the available data. The results showed the occurrence of Notacanthus arrontei n. sp. from the Iberian Peninsula and highlighted several taxonomic conundrums regarding the Notacanthus genus. For instance, no significant differences were found between Notacanthus indicus and the recently described Notacanthus laccadiviensis, questioning its taxonomic status. Similarly, the result of the species delimitation molecular analysis coincided with previous DNA barcoding studies supporting the snubnosed spiny eel Notacanthus chemnitzii as a species complex that requires further research. Moreover, two unidentified records from the Indian Ocean were confirmed to belong to an unknown species pending formal description, and barcoding data show for the first time the occurrence of the shortfin spiny eel Notacanthus bonaparte in the Australia-New Zealand area. This research confirms the existence of important gaps in the knowledge of notacanthid fishes and represents a step forward toward a better understanding of their biological diversity.

Species delimitation of tea plants (Camellia sect. Thea) based on super-barcodes

BACKGROUND

The era of high throughput sequencing offers new paths to identifying species boundaries that are complementary to traditional morphology-based delimitations. De novo species delimitation using traditional or DNA super-barcodes serve as efficient approaches to recognizing putative species (molecular operational taxonomic units, MOTUs). Tea plants (Camellia sect. Thea) form a group of morphologically similar species with significant economic value, providing the raw material for tea, which is the most popular nonalcoholic caffeine-containing beverage in the world. Taxonomic challenges have arisen from vague species boundaries in this group.

RESULTS

Based on the most comprehensive sampling of C. sect. Thea by far (165 individuals of 39 morphospecies), we applied three de novo species delimitation methods (ASAP, PTP, and mPTP) using plastome data to provide an independent evaluation of morphology-based species boundaries in tea plants. Comparing MOTU partitions with morphospecies, we particularly tested the congruence of MOTUs resulting from different methods. We recognized 28 consensus MOTUs within C. sect. Thea, while tentatively suggesting that 11 morphospecies be discarded. Ten of the 28 consensus MOTUs were uncovered as morphospecies complexes in need of further study integrating other evidence. Our results also showed a strong imbalance among the analyzed MOTUs in terms of the number of molecular diagnostic characters.

CONCLUSION

This study serves as a solid step forward for recognizing the underlying species boundaries of tea plants, providing a needed evidence-based framework for the utilization and conservation of this economically important plant group.

Molecular and morphometric study of Brazilian populations of Psychodopygus davisi

In this study, we analysed the molecular and morphometric differences of several populations of the putative sand fly vector Psychodopygus davisi (Root, 1934) (Diptera, Psychodidae, Phlebotominae) in Brazil. We amplified the 658 base pair fragments of the DNA barcoding region-cytochrome c oxidase subunit 1 (COI) gene-for 57 specimens of P. davisi and three specimens of Psychodopygus claustrei (Abonnenc, Léger & Fauran, 1979). We merged our data with public sequences of the same species available from GenBank. Then, the combined dataset-87 sequences and 20 localities-was analysed using population structure analysis and different species delimitation approaches. Geometric morphometry of wings was performed for 155 specimens of P. davisi populations from the North, Midwest and Southeast Brazilian regions, analysing the differences in centroid sizes and canonical variates. Molecular analysis indicated high intraspecific genetic distance values for P. davisi (maximum p distance = 5.52%). All algorithms identified P. davisi and P. claustrei as distinct molecular taxonomic units, despite the low interspecific distance (p distance to the nearest neighbour = 4.79%). P. davisi sequences were split into four genetic clusters by population structure analysis and at least five genetic lineages using intermediate scenarios of the species delimitation algorithms. The species validation analysis of BPP strongly supported the five-species model in our dataset. We found high genetic diversity in this taxon, which is in agreement with its wide geographic distribution in Brazil. Furthermore, the wing analysis showed that specimens from the Southeast Region of Brazil are different from those in the North and the Midwest. The evolutionary patterns of P. davisi populations in Brazil suggest the presence of candidate species, which need to be validated in future studies using a more comprehensive approach with both genomic data and morphological characters.

A new species of Raorchestes (Anura, Rhacophoridae) from Yunnan Province, China

A new bush frog species is described from Yunnan, China, based on phylogenetic analyses, species delimitation analyses, and morphological comparisons. Raorchesteshekouensissp. nov. is distinguished from all other congeners by a combination of 11 morphological characters. The new species brings the current number of Raorchestes species in China to ten, nine of which are distributed in Yunnan. Molecular analyses supported an unnamed lineage previously recorded as "Raorchestesgryllus" in northern Vietnam. Further studies including additional samples are necessary to clarify the species diversity and boundaries of Raorchestes in China and Indochina.

Species identification through deep learning and geometrical morphology in oaks (Quercus spp.): Pros and cons

Plant phenotypic characteristics, especially leaf morphology of leaves, are an important indicator for species identification. However, leaf shape can be extraordinarily complex in some species, such as oaks. The great variation in leaf morphology and difficulty of species identification in oaks have attracted the attention of scientists since Charles Darwin. Recent advances in discrimination technology have provided opportunities to understand leaf morphology variation in oaks. Here, we aimed to compare the accuracy and efficiency of species identification in two closely related deciduous oaks by geometric morphometric method (GMM) and deep learning using preliminary identification of simple sequence repeats (nSSRs) as a prior. A total of 538 Asian deciduous oak trees, 16 Q. aliena and 23 Q. dentata populations, were firstly assigned by nSSRs Bayesian clustering analysis to one of the two species or admixture and this grouping served as a priori identification of these trees. Then we analyzed the shapes of 2328 leaves from the 538 trees in terms of 13 characters (landmarks) by GMM. Finally, we trained and classified 2221 leaf-scanned images with Xception architecture using deep learning. The two species can be identified by GMM and deep learning using genetic analysis as a priori. Deep learning is the most cost-efficient method in terms of time-consuming, while GMM can confirm the admixture individuals' leaf shape. These various methods provide high classification accuracy, highlight the application in plant classification research, and are ready to be applied to other morphology analysis.

New sand fly (Diptera, Psychodidae) records and COI DNA barcodes in the state of Maranhão, Eastern Amazon, Brazil

The sand fly fauna and the usefulness of the DNA barcoding fragment of the cytochrome c oxidase subunit I (COI) gene were accessed in a forest fragment in the municipality of Governador Newton Bello, state of Maranhão, Brazil. We performed entomological collections in three independent campaigns in May and October 2021, and January 2023. Sand flies were morphologically-identified and then DNA barcoded. Sequences were deposited and analyzed in the BOLD System Database, and various species delimitation algorithms, to assess whether DNA sequences merge into taxonomic units in accordance with nominal species. In total, 1,524 sand flies were collected, comprising 32 nominal species. Nyssomyia antunesi was the most abundant species (31.5 %), followed by Psychodopygus davisi (27 %). We reported for the first time in the state of Maranhão, the presence of Lutzomyia evangelistai, Lutzomyia sherlocki, Pressatia equatorialis, and Psathyromyia barrettoi. We amplified and analyzed 67 COI barcodes of 23 species, which were merged with conspecific sequences extracted from GenBank. The maximum intraspecific p distances ranged from 0.0 % to 14.74 %, while the distances to the nearest neighbor varied from 1.67 % to 13.64 %. The phylogenetic gene tree and species delimitation tools clustered sequences into well-supported clades/clusters for each nominal species, except for Pressatia choti/Pr. equatorialis, which have the lowest interspecific genetic distance (1.67 %). We sequenced for the first time COI barcodes of Brumptomyia brumpti, Evandromyia monstruosa, Micropygomyia rorotaensis, Micropygomyia pilosa, Pintomyia christenseni, Pintomyia pacae, Pr. equatorialis, Pa. barrettoi, and Psathyromyia hermanlenti, which will be useful for further molecular identification and classification proposals of Neotropical species. This study updated the current list of the sand fly fauna for the state of Maranhão to 97, and demonstrated that COI barcodes are useful for specific identification.

Phylogenomics and species delimitation of an abundant and little-studied Amazonian forest spiny rat

Species delimitation studies based on integrating different datasets such as genomic, morphometric, and cytogenetics data are rare in studies focused on Neotropical rodents. As a consequence, the evolutionary history of most of these genera remains poorly understood. Proechimys is a highly diverse and widely distributed genus of Neotropical spiny rats with unique traits like multiple sympatry, micro-habitat segregation, and fuzzy species limits. Here, we applied RAD-Seq to infer the phylogenetic relationships, estimate the species boundaries, and estimate the divergence times for Proechimys, one of the most common and least studied small mammals in the Amazon. We tested whether inferred lineages in the phylogenetic trees could be considered distinct species based on the genomic dataset and morphometric data. Analyses revealed the genus is not monophyletic, with Proechimys hoplomyoides sister to a group of Hoplomys gymnurus + all other Proechimys species, contesting the generic status of Hoplomys. There are five main clades in Proechimys stricto sensu (excluding H. gymnurus and P. hoplomyoides). Species delimitation analyses supported 25 species within the genus Proechimys. The five main clades in Proechimys stricto sensu also showed similar ages for their origins, and two rapid diversification events were identified in the Early Pliocene and in the Early Pleistocene. Most cases of sympatry in Proechimys occur among species from the different main clades, and although Proechimys is an inhabitant of the Amazon, three species occupied the Cerrado biome during the Pleistocene. We could associate available nominal taxon, cytogenetics information, and DNA sequences in Genbank to most of the 25 species we hypothesized from our delimitation analyses. Based on our analyses, we estimate that eight forms represent putative new species that need a taxonomic revision.

Integrative species delimitation uncovers hidden diversity within the Pithecopus hypochondrialis species complex (Hylidae, Phyllomedusinae) and its phylogeography reveals Plio-Pleistocene connectivity among Neotropical savannas

Despite their limited vagility and pronounced habitat heterogeneity in the tropics, many anuran species have unexpectedly extensive geographic ranges. One prominent example of this phenomenon is Pithecopus hypochondrialis, which is found in the Cerrado, Guianan savanna, and Llanos domains, as well as isolated tracts of savanna and open habitat within the Amazon Forest. The present study employs an integrative species delimitation approach to test the hypothesis that P. hypochondrialis is in fact a species complex. We also reconstruct the relationships among the lineages delimited here and other Pithecopus species. In this study, we employ Ecological Niche Modelling (ENM) and spatiotemporal phylogeographic reconstruction approaches to evaluate a multitude of scenarios of connectivity across the Neotropical savannas. We identified three divergent lineages, two of which have been described previously. The lineages were allocated to a lowland Pithecopus clade, although the relationships among these lineages are weakly supported. Both the ENM and the phylogeographic reconstruction highlight the occurrence of periods of connectivity among the Neotropical savannas over the course of the Pliocene and Pleistocene epochs. These processes extended from eastern Amazonia to the northern coast of Brazil. The findings of the present study highlight the presence of hidden diversity within P. hypochondrialis, and reinforce the need for a comprehensive taxonomic review. These findings also indicate intricate and highly dynamic patterns of connectivity across the Neotropical savannas that date back to the Pliocene.

DNA barcoding of the genus Gampsocleis (Orthoptera, Tettigoniidae) from China

DNA barcoding is a useful addition to the traditional morphology-based taxonomy. A ca. 650 bp fragment of the 5' end of mitochondrial cytochrome c oxidase subunit I (hereafter COI-5P) DNA barcoding was sued as a practical tool for Gampsocleis species identification. DNA barcodes from 889 specimens belonging to 8 putative Gampsocleis species was analyzed, including 687 newly generated DNA barcodes. These barcode sequences were clustered/grouped into Operational Taxonomic Units (OTUs) using the criteria of five algorithms, namely Barcode Index Number (BIN) System, Assemble Species by Automatic Partitioning (ASAP), a Java program uses an explicit, determinate algorithm to define Molecular Operational Taxonomic Unit (jMOTU), Generalized Mixed Yule Coalescent (GMYC), and Bayesian implementation of the Poisson Tree Processes model (bPTP). The Taxon ID Tree grouped sequences of morphospecies and almost all MOTUs in distinct nonoverlapping clusters. Both long- and short-winged Gampsocleis species are reciprocally monophyletic in the Taxon ID Tree. In BOLD, 889 barcode sequences are assigned to 17 BINs. The algorithms ASAP, jMOTU, bPTP and GMYC clustered the barcode sequences into 6, 13, 10, and 23 MOTUs, respectively. BIN, ASAP, and bPTP algorithm placed three long-winged species, G. sedakovii, G. sinensis and G. ussuriensis within the same MOTU. All species delimitation algorithms split two short-winged species,G. fletcheri and G. gratiosa into at least two MOTUs each, except for ASAP algorithm. More detailed molecular and morphological integrative studies are required to clarify the status of these MOTUs in the future.

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.

Phylogenomics of the narrowly endemic Eurycheilichthys (Siluriformes: Loricariidae): Sympatric species with non-sister relationships suggest mainly allopatric speciation

Armored catfishes of the genus Eurycheilichthys are endemic to Southern Brazil and Misiones (Argentina) comprising nine species of small size, with a high degree of sympatry and species diversity distributed in two river basins. Here we use new genome-wide data to infer a species phylogeny and test species boundaries for this poorly known group. We estimate 1) the phylogenetic relationships of the species of Eurycheilichthys based on 29,350 loci in 65 individuals of nine species plus outgroups, and 2) the population structure and differentiation based on 43,712 loci and 62 individuals to estimate how geography may have acted on speciation and formation of the sympatric species groups. Analyses support the monophyly of the genus and suggest two species-inclusive clades (East and West) with high support and very recently diverged species. Western clade contains E. limulus (from upper Jacuí River basin) that is sister to Western species of the Taquari-Antas basin plus E. paucidens. The Eastern clade contains E. pantherinus (from Uruguay River basin) sister to the Eastern species of the Taquari-Antas basin E. coryphaenus, plus the central-distributed species E. planus and E. vacariensis, and the more widely-distributed species E. luisae. Eurycheilichthys luisae is not monophyletic and may contain one or more cryptic species or hybrid individuals. A stronger diversity on structure of lineages on the Taquari-Antas, when compared to upper Uruguay and Jacuí River basins, and the fact that most of the sympatrically distributed taxa have non-sister relationships suggest a scenario of mainly allopatric speciation and may indicate a more dynamic landscape with headwater capture events among these tributaries.

Patterns of speciation in a parapatric pair of Saturnia moths as revealed by target capture

The focus of this study has been to understand the evolutionary relationships and taxonomy of a widely distributed parapatric species pair of wild silk moths in Europe: Saturnia pavonia and Saturnia pavoniella (Lepidoptera: Saturniidae). To address species delimitation in these parapatric taxa, target enrichment and mtDNA sequencing was employed alongside phylogenetic, admixture, introgression, and species delimitation analyses. The dataset included individuals from both species close to and farther away from the contact zone as well as two hybrids generated in the lab. Nuclear markers strongly supported both S. pavonia and S. pavoniella as two distinct species, with hybrids forming a sister group to S. pavoniella. However, the Maximum Likelihood (ML) tree generated from mtDNA sequencing data presented a different picture, showing both taxa to be phylogenetically intermixed. This inconsistency is likely attributable to mitonuclear discordance, which can arise from biological factors (e.g., introgressive hybridization and/or incomplete lineage sorting). Our analyses indicate that past introgressions have taken place, but that there is no evidence to suggest an ongoing admixture between the two species, demonstrating that the taxa have reached full postzygotic reproductive isolation and hence represent two distinct biological species. Finally, we discuss our results from an evolutionary point of view taking into consideration the past climatic oscillations that have likely shaped the present dynamics between the two species. Overall, our study demonstrates the effectiveness of the target enrichment approach in resolving shallow phylogenetic relationships under complex evolutionary circumstances and that this approach is useful in establishing robust and well-informed taxonomic delimitations involving parapatric taxa.

Molecular delimitation methods validate morphologically similar species of red snappers (Perciformes: Lutjanidae)

In this study, we tested the taxonomic validation of red snappers species (Southern red snapper Lutjanus purpureus; Silk snapper L. vivanus; Blackfin snapper L. buccanella; and Pacific red snapper L. peru) based on comparative analysis, using four methods for species delimitation. These methods were based on either genetic similarity or phylogenetic trees inferred from two mitochondrial (Cytochrome b and D-loop) and two nuclear (Myostatin and S7 introns) markers. On one hand, the genetic results corroborated the presence of four red snapper species, confirming their taxonomic validation despite their remarkable morphological similarity. On the other hand, few incongruencies in the species delimitation methods were observed according to the phylogenetic reconstruction method (maximum likelihood or Bayesian inference) when using. Based on the phylogenetic results, L. buccanella should represent a more ancient lineage in relation to the clade that encompasses L. purpureus, L. peru and L. vivanus. The single-locus phylogenetic analysis based on Cytb recovered each the red snapper species as a well-supported clade. Overall, this study provided a DNA-based validation of the traditional morphological taxonomy of red snappers.

Species Delimitation in a Polyploid Group of Iberian Jasione (Campanulaceae) Unveils Coherence between Cryptic Speciation and Biogeographical Regionalization

Groups with morphological stasis are an interesting framework to address putative cryptic species that may be hidden behind traditional taxonomic treatments, particularly when distribution ranges suggest disjunct and environmentally heterogeneous biogeographic patterns. New hypotheses of delimitation of evolutionary independent units can lead to the identification of different biogeographic processes, laying the foundation to investigate their historical and ecological significance. Jasione is a plant genus with a distribution centered in the Mediterranean basin, characterized by significant morphological stasis. Within the western Mediterranean J. gr. crispa species complex, J. sessiliflora s.l. and allied taxa form a distinct group, occupying environmentally diverse regions. At least two ploidy levels, diploid and tetraploid, are known to occur in the group. The internal variability is assessed with phylogenetic tools, viz. GMYC and ASAP, for species delimitation. The results are compared with other lines of evidence, including morphology and cytology. The fitting of distribution patterns of the inferred entities to chorological subprovinces is also used as a biogeographical and environmental framework to test the species hypothesis. Despite the scarcity of diagnostic morphological characters in the group, phylogenetic delimitation supports the description of at least one cryptic species, a narrow endemic in the NE Iberian Peninsula. Moreover, the results support the segregation of a thermophilic group of populations in eastern Iberia from J. sessiliflora. Ploidy variation from a wide geographical survey supports the systematic rearrangement suggested by species delimitation. Taxonomic reorganization in J. sessiliflora s.l. would allow ecological interpretations of distribution patterns in great accordance with biogeographical regionalization at the subprovince level, supporting geobotanical boundaries as a framework to interpret species ecological coherence of cryptic lineages. These results suggest that species differentiation, together with geographic isolation and polyploidization, is associated with adaptation to different environments, shifting from more to less thermophilic conditions. Thus, the recognition of concealed evolutionary entities is essential to correctly interpret biogeographical patterns in regions with a complex geologic and evolutionary history, such as the Mediterranean basin, and biogeographical units emerge as biologically sound frameworks to test the species hypothesis.

Plastome phylogenomics and morphological traits analyses provide new insights into the phylogenetic position, species delimitation and speciation of Triplostegia (Caprifoliaceae)

BACKGROUND

The genus Triplostegia contains two recognized species, T. glandulifera and T. grandiflora, but its phylogenetic position and species delimitation remain controversial. In this study, we assembled plastid genomes and nuclear ribosomal DNA (nrDNA) cistrons sampled from 22 wild Triplostegia individuals, each from a separate population, and examined these with 11 recently published Triplostegia plastomes. Morphological traits were measured from herbarium specimens and wild material, and ecological niche models were constructed.

RESULTS

Triplostegia is a monophyletic genus within the subfamily Dipsacoideae comprising three monophyletic species, T. glandulifera, T. grandiflora, and an unrecognized species Triplostegia sp. A, which occupies much higher altitude than the other two. The new species had previously been misidentified as T. glandulifera, but differs in taproot, leaf, and other characters. Triplotegia is an old genus, with stem age 39.96 Ma, and within it T. glandulifera diverged 7.94 Ma. Triplostegia grandiflora and sp. A diverged 1.05 Ma, perhaps in response to Quaternary climate fluctuations. Niche overlap between Triplostegia species was positively correlated with their phylogenetic relatedness.

CONCLUSIONS

Our results provide new insights into the species delimitation of Triplostegia, and indicate that a taxonomic revision of Triplostegia is needed. We also identified that either rpoB-trnC or ycf1 could serve as a DNA barcode for Triplostegia.

Unsupervised machine learning for species delimitation, integrative taxonomy, and biodiversity conservation

Integrative taxonomy, combining data from multiple axes of biologically relevant variation, is a major goal of systematics. Ideally, such taxonomies will derive from similarly integrative species-delimitation analyses. Yet, most current methods rely solely or primarily on molecular data, with other layers often incorporated only in a post hoc qualitative or comparative manner. A major limitation is the difficulty of devising quantitative parametric models linking different datasets in a unified ecological and evolutionary framework. Machine Learning (ML) methods offer flexibility in this arena by easily learning high-dimensional associations between observations (e.g., individual specimens) across a wide array of input features (e.g., genetics, geography, environment, and phenotype) to delimit statistically meaningful clusters. Here, I implement an unsupervised method using Self-Organizing (or "Kohonen") Maps (SOMs) for such purposes. Recent extensions called "SuperSOMs" can integrate multiple layers, each of which exerts independent influence on a two-dimensional output grid via empirically estimated weights. The grid cells are then delimited into K distinct units that can be interpreted as species or other entities. I show empirical examples in salamanders (Desmognathus) and snakes (Storeria) with layers representing alleles, space, climate, and traits. Simulations reveal that the SuperSOM approach can detect K = 1, tends not to over-split, reflects contributions from all layers, and limits large layers (e.g., genetic matrices) from overwhelming other datasets, desirable properties addressing major concerns from previous studies. Finally, I suggest that these and similar methods could integrate conservation-relevant layers such as population trends and human encroachment to delimit management units from an explicitly quantitative framework grounded in the ecology and evolution of species limits and boundaries.

Trichophoromyia auraensis: evidence for cryptic species and first record in the state of Maranhão, Brazil

Trichophoromyia auraensis (Mangabeira, 1942) (Diptera, Psychodidae, Phlebotominae) has a wide geographic distribution in the western region of the Amazon biome, where it is a putative Leishmania vector. Here, we reported for the first time a population of this species in the Brazilian state of Maranhão, in the eastern Amazon, from which we DNA-barcoded and compared with previously processed specimens from Acre State, in the western Amazon. For this, we analyzed the DNA barcoding fragment (658 bp) of the mitochondrial cytochrome c oxidase subunit I (COI) gene and the nuclear internal transcribed spacer 2 (ITS2) of Trichophoromyia species using phylogenetic gene trees, and species delimitation algorithms. The analyses of COI barcodes showed high values of genetic distance (mean K2P = 5.17) and well-supported clades/MOTUs for the eastern and western populations of T. auraensis, which may indicate a possible complex of cryptic species. The western population of this taxon merged with the close-related sand fly Trichophoromyia velezbernali Posada-López, Galvis and Galati, 2018 from Colombia, which may be associated with the recent speciation history and introgression between these populations. These evidences should be evaluated with a more comprehensive sampling in terms of analyzed populations and molecular markers.

Genetic diversity of a short-ranged endemic terrestrial snail

The factors that influence population structure and connectivity are unknown for most terrestrial invertebrates but are of particular interest both for understanding the impacts of disturbance and for determining accurate levels of biodiversity and local endemism. The main objective of this study was to determine the historical patterns of genetic differentiation and contemporary gene flow in the terrestrial snail, Austrochloritis kosciuszkoensis (Shea & O. L. Griffiths, 2010). Snails were collected in the Mt Buffalo and Alpine National Parks in Victoria, in a bid to understand how populations of this species are connected both within continuous habitat and between adjacent, yet separate environments. Utilising both mitochondrial DNA (mtDNA) and single nucleotide polymorphism (SNP) data, the degree of population structure was determined within and between sites. Very high levels of genetic divergence were found between the Mt Buffalo and Alpine snails, with no evidence for genetic exchange detected between the two regions, indicating speciation has possibly occurred between the two regions. Our analyses of the combined mtDNA and nDNA (generated from SNPs) data have revealed patterns of genetic diversity that are consistent with a history of long-term isolation and limited connectivity. This history may be related to past cycles of changes to the climate over hundreds of thousands of years, which have, in part, caused the fragmentation of Australian forests. Within both regions, extremely limited gene flow between separate populations suggests that these land snails have very limited dispersal capabilities across existing landscape barriers, especially at Mt Buffalo: here, populations only 5 km apart from each other are genetically differentiated. The distinct genetic divergences and clearly reduced dispersal ability detected in this data explain the likely existence of at least two previously unnamed cryptic Austrochloritis species within a 30-50 km radius, and highlight the need for more concentrated efforts to understand population structure and gene flow in terrestrial invertebrates.

Integrative approaches to a revision of the liverwort in genus Aneura (Aneuraceae, Marchantiophyta) from Thailand

Background

The genus Aneura Dumort. is a simple thalloid liverwort with cosmopolitan distributions. Species circumscription is problematic in this genus due to a limited number of morphological traits. Two species are currently reported from Thailand, including A. maxima and A. pinguis. At the global scale, A. pinguis is considered a cryptic species, as the species contains several distinct genetic groups without clear morphological differentiation. At the same time, the identity of A. maxima remains unclear. In this work, we examined the level of diversity of Aneura species found in Thailand using both morphological and molecular data.

Methods

We measured the morphological traits and generated the molecular data (four markers: trnL-F, trnH-psbA, rbcL, and ITS2) from the Thai specimens. The concatenated dataset was then used to reconstruct phylogeny. Species delimitation with GMYC, bPTP, ASAP, and ABGD methods was performed to estimate the number of putative species within the genus.

Results

The samples of A. pinguis formed several clades, while A. maxima sequences from Poland were grouped in their clade and nested within another A. pinguis clade. We could not recover a sample of A. maxima from Thailand, even from the reported locality. Two putative species were detected among Thai Aneura samples. However, no morphological trait could distinguish the specimens from the two observed genetic groups.

Discussion

The previously observed paraphyletic nature of A. pinguis globally was also found among Thai samples, including several putative species. However, we could not confirm the identity of A. maxima from Thai specimens. The previous report could result from misidentification and problematic species circumscription within Aneura. The results highlighted the need to include multiple lines of evidence for the future taxonomic investigation of the group.

Genetic diversity analysis in the Brazilian Amazon reveals a new evolutionary lineage and new karyotype for the genus Mesomys (Rodentia, Echimyidae, Eumysopinae)

Morphological, molecular and chromosomal studies in the genera Lonchothrix and Mesomys have contributed to a better understanding of taxonomic design, phylogenetic relationships and karyotypic patterns. Recent molecular investigations have shown a yet undescribed diversity, suggesting that these taxa are even more diverse than previously assumed. Furthermore, some authors have questioned the limits of geographic distribution in the Amazon region for the species M. hispidus and M. stimulax. In this sense, the current study sought to understand the karyotypic evolution and geographic limits of the genus Mesomys, based on classical (G- and C-banding) and molecular cytogenetic analysis (FISH using rDNA 18S and telomeric probes) and through the sequencing of mitochondrial genes Cytochrome b (Cytb) and Cytochrome Oxidase-Subunit I (CO using phylogeny, species delimitation and time of divergence, from samples of different locations in the Brazilian Amazon. The species M. stimulax and Mesomys sp. presented 2n = 60/FN = 110, while M. hispidus presented 2n = 60/FN = 112, hitherto unpublished. Molecular dating showed that Mesomys diversification occurred during the Plio-Pleistocene period, with M. occultus diverging at around 5.1 Ma, followed by Mesomys sp. (4.1 Ma) and, more recently, the separation between M. hispidus and M. stimulax (3.5 Ma). The ABGD and ASAP species delimiters support the formation of 7 and 8 potential species of the genus Mesomys, respectively. Furthermore, in both analyzes Mesomys sp. was recovered as a valid species. Our multidisciplinary approach involving karyotypic, molecular and biogeographic analysis is the first performed in Mesomys, with the description of a new karyotype for M. hispidus, a new independent lineage for the genus and new distribution data for M. hispidus and M. stimulax.

How we study cryptic species and their biological implications: A case study from marine shelled gastropods

Methodological and biological considerations are intertwined when studying cryptic species. A potentially large component of modern biodiversity, the frequency of cryptic species among taxonomic groups is not well documented. The term "cryptic species" is imprecisely used in scientific literature, causing ambiguity when interpreting their evolutionary and ecological significance. This study reviews how cryptic species have been defined, discussing implications for taxonomy and biology, and explores these implications with a case study based on recently published literature on extant shelled marine gastropods. Reviewed gastropods were recorded by species. Records of cryptic gastropods were presented by authors with variable levels of confidence but were difficult to disentangle from inherent biases in the study effort. These complexities notwithstanding, most gastropod species discussed were not cryptic. To the degree that this review's sample represents extinct taxa, the results suggest that a high proportion of shelled marine gastropod species are identifiable for study in the fossil record. Much additional work is needed to provide a more adequate understanding of the relative frequency of cryptic species in shelled marine gastropods, which should start with more explicit definitions and targeted case studies.

Penicillium rhizophilum, a novel species in the section Exilicaulis isolated from the rhizosphere of sugarcane in Southwest Iran

During a survey of species diversity of Penicillium and Talaromyces in sugarcane (Saccharum officinarum) rhizosphere in the Khuzestan province of Iran [1], 195 strains were examined, from which 187 belonged to Penicillium (11 species) and eight to Talaromyces (one species). In the present study, three strains of Penicillium belonging to section Exilicaulis series Restricta, identified as P. restrictum by Ansari et al. [1], were subjected to a phylogenetic study. The multilocus phylogeny of partial β-tubulin, calmodulin and RNA polymerase II second largest subunit genes enabled the recognition of one new phylogenetic species that is here formally described as Penicillium rhizophilum sp. nov. This species is phylogenetically distinct in series Restricta, but it does not show significant morphological differences from other species previously classified in the series. Therefore, we here placed bias on the phylogenetic species concept. The holotype of Penicillium rhizophilum sp. nov. is IRAN 18169F and the ex-type culture is LA30T (=IRAN 4042CT=CBS 149737T).

Losing the forest for the tree? On the wisdom of subpopulation management

Animal habitats are changing around the world in many ways, presenting challenges to the survival of species. Zoo animal populations are also challenged by small population sizes and limited genetic diversity. Some ex situ populations are managed as subpopulations based on presumed subspecies or geographic locality and related concerns over genetic purity or taxonomic integrity. However, these decisions can accelerate the loss of genetic diversity and increase the likelihood of population extinction. Here I challenge the wisdom of subpopulation management, pointing out significant concerns in the literature with delineation of species, subspecies, and evolutionarily significant units. I also review literature demonstrating the value of gene flow for preserving adaptive potential, the often-misunderstood role of hybridization in evolution, and the likely overstated concerns about outbreeding depression, and preservation of local adaptations. I argue that the most effective way to manage animal populations for the long term be they in human care, in the wild, or if a captive population is being managed for reintroduction, is to manage for maximum genetic diversity rather than managing subpopulations focusing on taxonomic integrity, genetic purity, or geographic locale because selection in the future, rather than the past, will determine what genotypes and phenotypes are the most fit. Several case studies are presented to challenge the wisdom of subpopulation management and stimulate thinking about the preservation of genomes rather than species, subspecies, or lineages because those units evolved in habitats that are likely very different from those habitats today and in the future.

Abundant, diverse, unknown: Extreme species richness and turnover despite drastic undersampling in two closely placed tropical Malaise traps

Arthropods account for a large proportion of animal biomass and diversity in terrestrial systems, making them crucial organisms in our environments. However, still too little is known about the highly abundant and megadiverse groups that often make up the bulk of collected samples, especially in the tropics. With molecular identification techniques ever more evolving, analysis of arthropod communities has accelerated. In our study, which was conducted within the Global Malaise trap Program (GMP) framework, we operated two closely placed Malaise traps in Padang, Sumatra, for three months. We analyzed the samples by DNA barcoding and sequenced a total of more than 70,000 insect specimens. For sequence clustering, we applied three different delimitation techniques, namely RESL, ASAP, and SpeciesIdentifier, which gave similar results. Despite our (very) limited sampling in time and space, our efforts recovered more than 10,000 BINs, of which the majority are associated with "dark taxa". Further analysis indicates a drastic undersampling of both sampling sites, meaning that the true arthropod diversity at our sampling sites is even higher. Regardless of the close proximity of both Malaise traps (< 360 m), we discovered significantly distinct communities.

Consilience Across Multiple, Independent Genomic Data Sets Reveals Species in a Complex with Limited Phenotypic Variation

Species delimitation in the genomic era has focused predominantly on the application of multiple analytical methodologies to a single massive parallel sequencing (MPS) data set, rather than leveraging the unique but complementary insights provided by different classes of MPS data. In this study, we demonstrate how the use of two independent MPS data sets, a sequence capture data set and a single-nucleotide polymorphism (SNP) data set generated via genotyping-by-sequencing, enables the resolution of species in three complexes belonging to the grass genus Ehrharta, whose strong population structure and subtle morphological variation limit the effectiveness of traditional species delimitation approaches. Sequence capture data are used to construct a comprehensive phylogenetic tree of Ehrharta and to resolve population relationships within the focal clades, while SNP data are used to detect patterns of gene pool sharing across populations, using a novel approach that visualizes multiple values of K. Given that the two genomic data sets are independent, the strong congruence in the clusters they resolve provides powerful ratification of species boundaries in all three complexes studied. Our approach is also able to resolve a number of single-population species and a probable hybrid species, both of which would be difficult to detect and characterize using a single MPS data set. Overall, the data reveal the existence of 11 and five species in the E. setacea and E. rehmannii complexes, with the E. ramosa complex requiring further sampling before species limits are finalized. Despite phenotypic differentiation being generally subtle, true crypsis is limited to just a few species pairs and triplets. We conclude that, in the absence of strong morphological differentiation, the use of multiple, independent genomic data sets is necessary in order to provide the cross-data set corroboration that is foundational to an integrative taxonomic approach. [Species delimitation; genotyping-by-sequencing; population structure; integrative taxonomy; cryptic species; Ehrharta (Poaceae).].

Species delimitation and integrative taxonomy of the Reithrodontomys mexicanus (Rodentia: Cricetidae) cryptic complex

Species boundaries are difficult to establish in groups with very similar morphology. As an alternative, it has been suggested to integrate multiple sources of data to clarify taxonomic problems in taxa where cryptic speciation processes have been reported. This is the case of the harvest mouse Reithrodontomys mexicanus, which has a problematic taxonomy history as it is considered a complex species. Here, we evaluate the cryptic diversity of R. mexicanus using an integrative taxonomy approach in order to detect candidate lineages at the species level. The molecular analysis used one mitochondrial (cytb) and two nuclear (Fgb-I7 and IRBP) genes. Species hypotheses were suggested based on three molecular delimitation methods (mPTP, bGMYC, and STACEY) and cytb genetic distance values. Skull and environmental space differences between the delimited species were also tested to complement the discrimination of candidate species. Based on the consensus across the delimitation methods and genetic distance values, four species were proposed, which were mostly supported by morphometric and ecological data: R. mexicanus clade I, R. mexicanus clade IIA, R. mexicanus clade IIIA, and R. mexicanus clade IIIB. In addition, the evolutionary relationships between the species that comprise the R. mexicanus group were discussed from a phylogenetic approach. Our findings present important taxonomic implications for Reithrodontomys, as the number of known species for this genus increases. Furthermore, we highlight the importance of the use of multiple sources of data in systematic studies to establish robust delimitations between species considered taxonomically complex.

A Custom Regional DNA Barcode Reference Library for Lichen-Forming Fungi of the Intermountain West, USA, Increases Successful Specimen Identification

DNA barcoding approaches provide powerful tools for characterizing fungal diversity. However, DNA barcoding is limited by poor representation of species-level diversity in fungal sequence databases. Can the development of custom, regionally focused DNA reference libraries improve species-level identification rates for lichen-forming fungi? To explore this question, we created a regional ITS database for lichen-forming fungi (LFF) in the Intermountain West of the United States. The custom database comprised over 4800 sequences and represented over 600 formally described and provisional species. Lichen communities were sampled at 11 sites throughout the Intermountain West, and LFF diversity was characterized using high-throughput ITS2 amplicon sequencing. We compared the species-level identification success rates from our bulk community samples using our regional ITS database and the widely used UNITE database. The custom regional database resulted in significantly higher species-level assignments (72.3%) of candidate species than the UNITE database (28.3-34.2%). Within each site, identification of candidate species ranged from 72.3-82.1% using the custom database; and 31.5-55.4% using the UNITE database. These results highlight that developing regional databases may accelerate a wide range of LFF research by improving our ability to characterize species-level diversity using DNA barcoding.

Mitochondrial introgression and mito-nuclear discordance obscured the closely related species boundaries in Cletus Stål from China (Heteroptera: Coreidae)

Accurate taxonomy and delimitation are of great importance for pest control strategies and management programs. Here, we focus on Cletus (Insecta: Hemiptera: Coreidae), which includes many crop pests. The species boundaries still conflict and only cytochrome c oxidase subunit I (COI) barcoding has been previously used for molecular studies. We generated new mitochondrial genome and nuclear genome-wide SNPs to explore the species boundaries of 46 Cletus samples from China using multiple species delimitation approaches. All results recovered a monophyly with high support, except for two closely related species in clade I - C. punctiger and C. graminis. Mitochondrial data demonstrated admixture in clade I, while genome-wide SNPs unambiguously identified two separate species, which were confirmed by morphological classification. Inconsistent nuclear and mitochondrial data indicated mito-nuclear discordance. Mitochondrial introgression is the most likely explanation, and more extensive sampling and more comprehensive data are needed to ascertain a pattern. Accurate species delimitation will shed light on species status; thus, an accurate taxonomy is of particular concern, as there is a pressing need to implement precise control of agricultural pests and to perform further research on diversification.

Delimitation and species discovery in the Profundulidae fish family: Using genetic, environmental and morphologic data to address taxonomic uncertainty

The family Profundulidae includes some one of the most enigmatic freshwater fishes of Mesoamerica: despite many attempts, a robust phylogenetic framework to delimit species is lacking, mainly due to limited morphological variation within the group. The accumulation of molecular data of profundulid fishes has led to advances in the description of new taxa, but relatively less progress has been made estimating evolutionary and phylogenetic relationships for this fish family. Here, we adopt an integrative taxonomy approach including the use of nuclear and mitochondrial DNA sequences, morphometric and ecological data, to test species boundaries in profundulid fishes in the westernmost area of their known distribution range in the states of Guerrero and Oaxaca, Mexico. Using a combination of methods for species discovery and validation based on Bayesian gene tree topologies, our analyses support the delimitation of 15 valid species of profundulid fishes - a combination of previously described species validated by this study, the synonymy of unsupported taxa, and the description of two new species. Using species delimitation methods, examination of phenotypic variation, and ecological niche characterization, we also identify five potentially new lineages which require further evidence to be erected as new species. We demonstrate that the use of an integrative taxonomy approach provides a robust methodology to delimit species in a taxonomically complex group like Profundulidae. Accurate taxonomic and ecological information is crucial for the conservation of these microendemic fishes, as several species are endangered.

Improving the COI DNA barcoding library for Neotropical phlebotomine sand flies (Diptera: Psychodidae)

Sand fly species are traditionally identified using morphological traits, though this method is hampered by the presence of cryptic species. DNA barcoding is a widely used tool in the case of insects of medical importance, where it is necessary to know quickly which species are present in a transmission area. Here, we assess the usefulness of mitochondrial cytochrome c oxidase subunit I (COI) DNA barcoding as a practical tool for species identification, correct assignment of isomorphic females, and to evaluate the detection of cryptic diversity that occurs in the same species. A fragment of the COI gene was used to generate 156 new barcode sequences for sand flies from different countries of the Neotropical region, mainly Colombia, which had been identified morphologically as 43 species. The sequencing of the COI gene allowed the detection of cryptic diversity within species and correctly associated isomorphic females with males identified by morphology. The maximum intraspecific genetic distances ranged from 0 to 8.32% and 0 to 8.92% using uncorrected p distances and the Kimura 2-parameter (K2P) model, respectively. The minimum interspecific distance (nearest neighbor) for each species ranged from 1.5 to 14.14% and 1.51 to 15.7% using p and K2P distances, respectively. Three species had more than 3% maximum intraspecific distance: Psychodopygus panamensis, Micropygomyia cayennensis cayennensis, and Pintomyia evansi. They also were split into at least two molecular operational taxonomic units (MOTUs) each, using different species delimitation algorithms. Regarding interspecific genetic distances, the species of the genera Nyssomyia and Trichophoromyia generated values lower than 3% (except Nyssomyia ylephiletor and Ny. trapidoi). However, the maximum intraspecific distances did not exceed these values, indicating the presence of a barcode gap despite their proximity. Also, nine sand fly species were DNA barcoded for the first time: Evandromyia georgii, Lutzomyia sherlocki, Ny. ylephiletor, Ny. yuilli pajoti, Psathyromyia punctigeniculata, Sciopemyia preclara, Trichopygomyia triramula, Trichophoromyia howardi, and Th. velezbernali. The COI DNA barcode analysis enabled the correct delimitation of several Neotropical sand fly species from South and Central America and raised questions about the presence of cryptic species for some taxa, which should be further assessed.

Integrative species delimitation helps to find the hidden diversity of the leaf-litter frog Ischnocnema manezinho (Garcia, 1996) (Anura, Brachycephalidae), endemic to the southern Atlantic Forest

Background

The delimitation of cryptic species is a challenge for biodiversity conservation. Anurans show high cryptic diversity levels, and molecular species delimitation methods could help identify putative new species. Additionally, species delimitation approaches can provide important results for cryptic species conservation, with integrative methods adding robustness to results. Ischnocnema manezinho was described from Santa Catarina Island (SCI), southern Brazil. More recently, some inventories indicated continental populations supposedly similar in morphology to it. If these records are confirmed as I. manezinho, it would likely change its endangered status on National Red List, removing the species from conservation agendas. We investigated the threatened frog Ischnocnema manezinho, to evaluate if the continental populations belong to this species or if they form an undescribed species complex.

Methods

We used coalescent, distance, and allele-sharing-based species delimitation methods and integrative analyses of morphometric and bioacoustics traits to test evolutionary independence between I. manezinho from SCI, Arvoredo Island, and continental populations.

Results

Ischnocnema manezinho is restricted to Santa Catarina Island, while the five remaining lineages should be further investigated through a taxonomic review. Our results point to a small geographic range of Ischnocnema manezinho. Additionally, the species occurs in isolated fragments of forest in SCI surrounded by expanding urban areas, confirming its status as Endangered. Thus, the protection and monitoring of I. manezinho and the taxonomic description of the continental and Arvoredo Island candidate species should be priorities.

Biogeography and evolution of social parasitism in Australian Myrmecia bulldog ants revealed by phylogenomics

Studying the historical biogeography and life history transitions from eusocial colony life to social parasitism contributes to our understanding of the evolutionary mechanisms generating biodiversity in eusocial insects. The ants in the genus Myrmecia are a well-suited system for testing evolutionary hypotheses about how their species diversity was assembled through time because the genus is endemic to Australia with the single exception of the species M. apicalis inhabiting the Pacific Island of New Caledonia, and because at least one social parasite species exists in the genus. However, the evolutionary mechanisms underlying the disjunct biogeographic distribution of M. apicalis and the life history transition(s) to social parasitism remain unexplored. To study the biogeographic origin of the isolated, oceanic species M. apicalis and to reveal the origin and evolution of social parasitism in the genus, we reconstructed a comprehensive phylogeny of the ant subfamily Myrmeciinae. We utilized Ultra Conserved Elements (UCEs) as molecular markers to generate a comprehensive molecular genetic dataset consisting of 2,287 loci per taxon on average for 66 out of the 93 known Myrmecia species as well as for the sister lineage Nothomyrmecia macrops and selected outgroups. Our time-calibrated phylogeny inferred that: (i) stem Myrmeciinae originated during the Paleocene ∼58 Ma ago; (ii) the current disjunct biogeographic distribution of M. apicalis was driven by long-distance dispersal from Australia to New Caledonia during the Miocene ∼14 Ma ago; (iii) the single social parasite species, M. inquilina, evolved directly from one of the two known host species, M. nigriceps, in sympatry via the intraspecific route of social parasite evolution; and (iv) 5 of the 9 previously established taxonomic species groups are non-monophyletic. We suggest minor changes to reconcile the molecular phylogenetic results with the taxonomic classification. Our study enhances our understanding of the evolution and biogeography of Australian bulldog ants, contributes to our knowledge about the evolution of social parasitism in ants, and provides a solid phylogenetic foundation for future inquiries into the biology, taxonomy, and classification of Myrmeciinae.

Betula mcallisteri sp. nov. (sect. Acuminatae, Betulaceae), a new diploid species overlooked in the wild and in cultivation, and its relation to the widespread B. luminifera

Taxa are traditionally identified using morphological proxies for groups of evolutionarily isolated populations. These proxies are common characters deemed by taxonomists as significant. However, there is no general rule on which character or sets of characters are appropriate to circumscribe taxa, leading to discussions and uncertainty. Birch species are notoriously hard to identify due to strong morphological variability and factors such as hybridization and the existence of several ploidy levels. Here, we present evidence for an evolutionarily isolated line of birches from China that are not distinguishable by traditionally assumed taxon recognition proxies, such as fruit or leaf characters. We have discovered that some wild material in China and some cultivated in the Royal Botanic Gardens Edinburgh, formerly recognized as Betula luminifera, differ from other individuals by having a peeling bark and a lack of cambial fragrance. We use restriction site-associated DNA sequencing and flow cytometry to study the evolutionary status of the unidentified Betula samples to assess the extent of hybridization between the unidentified Betula samples and typical B. luminifera in natural populations. Molecular analyses show the unidentified Betula samples as a distinct lineage and reveal very little genetic admixture between the unidentified samples and B. luminifera. This may also be facilitated by the finding that B. luminifera is tetraploid, while the unidentified samples turned out to be diploid. We therefore conclude that the samples represent a yet unrecognized species, which is here described as Betula mcallisteri.

Improving species delimitation for effective conservation: a case study in the endemic maple-leaf oak (Quercus acerifolia)

Species delimitation is challenging in lineages that exhibit both high plasticity and introgression. This challenge can be compounded by collection biases, which may downweight specimens morphologically intermediate between traditional species. Additionally, mismatch between named species and observable phenotypes can compromise species conservation. We studied the species boundaries of Quercus acerifolia, a tree endemic to Arkansas, U.S. We performed morphometric analyses of leaves and acorns from 527 field and 138 herbarium samples of Q. acerifolia and its close relatives, Q. shumardii and Q. rubra. We employed two novel approaches: sampling ex situ collections to detect phenotypic plasticity caused by environmental variation and comparing random field samples with historical herbarium samples to identify collection biases that might undermine species delimitation. To provide genetic evidence, we also performed molecular analyses on genome-wide SNPs. Quercus acerifolia shows distinctive morphological, ecological, and genomic characteristics, rejecting the hypothesis that Q. acerifolia is a phenotypic variant of Q. shumardii. We found mismatches between traditional taxonomy and phenotypic clusters. We detected underrepresentation of morphological intermediates in herbarium collections, which may bias species discovery and recognition. Rare species conservation requires considering and addressing taxonomic problems related to phenotypic plasticity, mismatch between taxonomy and morphological clusters, and collection biases.

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.

The phylogenetic structure and coalescent species delimitation of an endemic trapdoor spider genus, Stasimopus (Araneae, Mygalomorphae, Stasimopidae) in the Karoo region of South Africa

The Karoo region of South Africa is a unique and sensitive ecosystem which is facing pressure for development due to economic incentives such as mining, farming and shale gas exploration. The species diversity of many taxa in the area is largely unknown. A phylogenetic analysis of the cork-lid trapdoor spider genus, Stasimopus (Stasimopidae) was undertaken in order to gain insight into the relationships between the species that may be present in the area. The species within Stasimopus are challenging to identify and define using traditional morphological methods due to a high degree of morphological conservatism within the genus. For this reason, multiple coalescent based species delimitation methods were used to attempt to determine the species present for Stasimopus in the region which was tested against the morphological identifications and genetic clades (based on CO1, 16S and EF-1ɣ). We tested single-locus methods Automatic Barcode Gap Discovery (ABGD), Bayesian implementation of Poisson Tree Processes (bPTP) and General Mixed Yule- Coalescent (GMYC), as well as multi-locus Brownie. The phylogenetic analysis of Stasimopus in the Karoo showed that there is a high degree of genetic diversity within the genus. The species delimitation results proved unfruitful for the genus, as they appear to delimit population structure rather than species for most methods. Alternative methods should be investigated to aid in the identification of the species in order truly understand the species diversity of the genus.

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.

A systematic revision of the Shovel-nosed Salamander (Plethodontidae: Desmognathus marmoratus), with re-description of the related D. aureatus and D. intermedius

Shovel-nosed Salamanders, Desmognathus marmoratus (Moore, 1899), were long thought to represent a single species from the southern Appalachian Mountains of the eastern United States, ranging from northeastern Georgia to extreme southwestern Virginia. These populations have a highly derived ecomorphology, being fully aquatic with a specialized flattened and elongated phenotype adapted to rocky riffle zones in fast-flowing, high-gradient mountain streams. Because of this, they were originally described in a separate genus, Leurognathus Moore, 1899. Four additional species or subspecies were described from 1928-1956 based on regional geographic variation in phenotype before being synonymized with L. marmoratus in 1962, which was reassigned to Desmognathus in 1996. Molecular analyses subsequently revealed four distinct candidate lineages in two distantly related clades, which were recently re-delimited into three species. These are D. aureatus (Martof, 1956) from northeastern Georgia, D. intermedius (Pope, 1928) from western North Carolina, and D. marmoratus from northwestern North Carolina. We provide a systematic revision of these taxa, which do not represent a natural group but instead exhibit convergent phenotypes across multiple species, potentially driven by ancient episodes of adaptive introgression between ancestral lineages. Our recent fieldwork revealed an astonishingly disjunct and morphologically distinct population of D. marmoratus in the New River Gorge of West Virginia, which were previously confused with D. kanawha Pyron and Beamer, 2022. This locality is ~120 airline km away from the nearest populations of D. marmoratus in Virginia. No Shovel-nosed Salamanders have ever been found in the New River drainage during our extensive previous explorations or credibly reported in museum specimens or the literature. Additional cryptic populations of these taxa may remain.

Single-specimen systematics resolves the phylogeny and diversity conundrum of enigmatic crustacean y-larvae

Resolving the evolutionary history of organisms is a major goal in biology. Yet for some taxa the diversity, phylogeny, and even adult stages remain unknown. The enigmatic crustacean "y-larvae" (Facetotecta) is one particularly striking example. Here we use extensive video-imaging and single-specimen molecular sequencing of >200 y-larval specimens to comprehensively explore for the first time their evolutionary history and diversity. This integrative approach revealed five major clades of Facetotecta, four of which encompass a considerable larval diversity. Whereas morphological analyses recognized 35 y-naupliar "morphospecies", molecular species delimitation analyses suggested the existence of between 88 and 127 species. The phenotypic and genetic diversity between the morphospecies suggests that a more elaborate classification than the current one-genus approach is needed. Morphology and molecular data were highly congruent at shallower phylogenetic levels, but no morphological synapomorphies could be unambiguously identified for major clades, which mostly comprise both planktotrophic and lecithotrophic y-nauplii. We argue that lecithotrophy arose several times independently whereas planktotrophic y-nauplii, which are structurally more similar across clades, most likely display the ancestral feeding mode of Facetotecta. We document a remarkably complex and highly diverse phylogenetic backbone for a taxon of marine crustaceans, the full life cycle of which remains a mystery.

High Andean Steppes of Southern Chile Contain Little-Explored Peltigera Lichen Symbionts

Peltigera lichens can colonize extreme habitats, such as high-elevation ecosystems, but their biodiversity is still largely unknown in these environments, especially in the southern hemi- sphere. We assessed the genetic diversity of mycobionts and cyanobionts of 60 Peltigera lichens collected in three high Andean steppes of southern Chile using LSU, β-tubulin, COR3 and ITS loci for mycobionts, and SSU and rbcLX loci for cyanobionts. We obtained 240 sequences for the different mycobiont markers and 118 for the cyanobiont markers, including the first report of β-tubulin sequences of P. patagonica through modifying a previously designed primer. Phylogenetic analyses, ITS scrutiny and variability of haplotypes were used to compare the sequences with those previously reported. We found seven mycobiont species and eleven cyanobiont haplotypes, including considerable novel symbionts. This was reflected by ~30% of mycobionts and ~20% of cyanobionts haplotypes that yielded less than 99% BLASTn sequence identity, 15 new sequences of the ITS1-HR, and a putative new Peltigera species associated with 3 Nostoc haplotypes not previously reported. Our results suggest that high Andean steppe ecosystems are habitats of unknown or little-explored lichen species and thus valuable environments to enhance our understanding of global Peltigera biodiversity.

Gene flow assessment helps to distinguish strong genomic structure from speciation in an Iberian ant-eating spider

Although genomic data is boosting our understanding of evolution, we still lack a solid framework to perform reliable genome-based species delineation. This problem is especially critical in the case of phylogeographically structured organisms, with allopatric populations showing similar divergence patterns as species. Here, we assess the species limits and phylogeography of Zodarion alacre, an ant-eating spider widely distributed across the Iberian Peninsula. We first performed species delimitation based on genome-wide data and then validated these results using additional evidence. A commonly employed species delimitation strategy detected four distinct lineages with almost no admixture, which present allopatric distributions. These lineages showed ecological differentiation but no clear morphological differentiation, and evidence of introgression in a mitochondrial barcode. Phylogenomic networks found evidence of substantial gene flow between lineages. Finally, phylogeographic methods highlighted remarkable isolation by distance and detected evidence of range expansion from south-central Portugal to central-north Spain. We conclude that despite their deep genomic differentiation, the lineages of Z. alacre do not show evidence of complete speciation. Our results likely shed light on why Zodarion is among the most diversified spider genera despite its limited distribution and support the use of gene flow evidence to inform species boundaries.

To explore strange new worlds - The diversification in Tremella caloplacae was linked to the adaptive radiation of the Teloschistaceae

Lichenicolous fungi are a heterogeneous group of organisms that grow exclusively on lichens, forming obligate associations with them. It has often been assumed that cospeciation has occurred between lichens and lichenicolous fungi, but this has been seldom analysed from a macroevolutionary perspective. Many lichenicolous species are rare or are rarely observed, which results in frequent and large gaps in the knowledge of the diversity of many groups. This, in turn, hampers evolutionary studies that necessarily are based on a reasonable knowledge of this diversity. Tremella caloplacae is a heterobasidiomycete growing on various hosts from the lichen-forming family Teloschistaceae, and evidence suggests that it may represent a species complex. We combine an exhaustive sampling with molecular and ecological data to study species delimitation, cophylogenetic events and temporal concordance of this association. Tremella caloplacae is here shown to include at least six distinct host-specific lineages (=putative species). Host switch is the dominant and most plausible event influencing diversification and explaining the coupled evolutionary history in this system, although cospeciation cannot be discarded. Speciation in T. caloplacae would therefore have occurred coinciding with the rapid diversification - by an adaptive radiation starting in the late Cretaceous - of their hosts. New species in T. caloplacae would have developed as a result of specialization on diversifying lichen hosts that suddenly offered abundant new ecological niches to explore or adapt to.