Bayesian species delimitation combining multiple genes and traits in a unified framework

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.

Species delimitation 4.0: integrative taxonomy meets artificial intelligence

Although species are central units for biological research, recent findings in genomics are raising awareness that what we call species can be ill-founded entities due to solely morphology-based, regional species descriptions. This particularly applies to groups characterized by intricate evolutionary processes such as hybridization, polyploidy, or asexuality. Here, challenges of current integrative taxonomy (genetics/genomics + morphology + ecology, etc.) become apparent: different favored species concepts, lack of universal characters/markers, missing appropriate analytical tools for intricate evolutionary processes, and highly subjective ranking and fusion of datasets. Now, integrative taxonomy combined with artificial intelligence under a unified species concept can enable automated feature learning and data integration, and thus reduce subjectivity in species delimitation. This approach will likely accelerate revising and unraveling eukaryotic biodiversity.

Species delimitation in Amana (Liliaceae): transcriptomes battle with evolutionary complexity

Species delimitation has long been a subject of controversy, and there are many alternative concepts and approaches used to define species in plants. The genus Amana (Liliaceae), known as "East Asian tulips" has a number of cryptic species and a huge genome size (1C = 21.48-57.35 pg). It also is intriguing how such a spring ephemeral genus thrives in subtropical areas. However, phylogenetic relationships and species delimitation within Amana are challenging. Here we included all species and 84 populations of Amana, which are collected throughout its distribution range. A variety of methods were used to clarify its species relationships based on a combination of morphological, ecological, genetic, evolutionary and phylogenetic species concepts. This evidence supports the recognition of at least 12 species in Amana. Moreover, we explored the complex evolutionary history within the genus and detected several historical hybridization and introgression events based on phylogenetic trees (transcriptomic and plastid), phylonetworks, admixture and ABBA-BABA analyses. Morphological traits have undergone parallel evolution in the genus. This spring ephemeral genus might have originated from a temperate region, yet finally thrives in subtropical areas, and three hypotheses about its adaptive evolution are proposed for future testing. In addition, we propose a new species, Amana polymorpha, from eastern Zhejiang Province, China. This research also demonstrates that molecular evidence at the genome level (such as transcriptomes) has greatly improved the accuracy and reasonability of species delimitation and taxon classification.

Tracing evolutionary trajectories in the presence of gene flow in South American temperate lizards (Squamata: Liolaemus kingii group)

Evolutionary processes behind lineage divergence often involve multidimensional differentiation. However, in the context of recent divergences, the signals exhibited by each dimension may not converge. In such scenarios, incomplete lineage sorting, gene flow, and scarce phenotypic differentiation are pervasive. Here, we integrated genomic (RAD loci of 90 individuals), phenotypic (linear and geometric traits of 823 and 411 individuals, respectively), spatial, and climatic data to reconstruct the evolutionary history of a speciation continuum of liolaemid lizards (Liolaemus kingii group). Specifically, we (a) inferred the population structure of the group and contrasted it with the phenotypic variability; (b) assessed the role of postdivergence gene flow in shaping phylogeographic and phenotypic patterns; and (c) explored ecogeographic drivers of diversification across time and space. We inferred eight genomic clusters exhibiting leaky genetic borders coincident with geographic transitions. We also found evidence of postdivergence gene flow resulting in transgressive phenotypic evolution in one species. Predicted ancestral niches unveiled suitable areas in southern and eastern Patagonia during glacial and interglacial periods. Our study underscores integrating different data and model-based approaches to determine the underlying causes of diversification, a challenge faced in the study of recently diverged groups. We also highlight Liolaemus as a model system for phylogeographic and broader evolutionary studies.

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.

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.

Discerning structure versus speciation in phylogeographic analysis of Seepage Salamanders (Desmognathus aeneus) using demography, environment, geography, and phenotype

Numerous mechanisms can drive speciation, including isolation by adaptation, distance, and environment. These forces can promote genetic and phenotypic differentiation of local populations, the formation of phylogeographic lineages, and ultimately, completed speciation. However, conceptually similar mechanisms may also result in stabilizing rather than diversifying selection, leading to lineage integration and the long-term persistence of population structure within genetically cohesive species. Processes that drive the formation and maintenance of geographic genetic diversity while facilitating high rates of migration and limiting phenotypic differentiation may thereby result in population genetic structure that is not accompanied by reproductive isolation. We suggest that this framework can be applied more broadly to address the classic dilemma of "structure" versus "species" when evaluating phylogeographic diversity, unifying population genetics, species delimitation, and the underlying study of speciation. We demonstrate one such instance in the Seepage Salamander (Desmognathus aeneus) from the southeastern United States. Recent studies estimated up to 6.3% mitochondrial divergence and four phylogenomic lineages with broad admixture across geographic hybrid zones, which could potentially represent distinct species supported by our species-delimitation analyses. However, while limited dispersal promotes substantial isolation by distance, microhabitat specificity appears to yield stabilizing selection on a single, uniform, ecologically mediated phenotype. As a result, climatic cycles promote recurrent contact between lineages and repeated instances of high migration through time. Subsequent hybridization is apparently not counteracted by adaptive differentiation limiting introgression, leaving a single unified species with deeply divergent phylogeographic lineages that nonetheless do not appear to represent incipient species.

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.

Complex Patterns of Diversification in the Gray Zone of Speciation: Model-Based Approaches Applied to Patagonian Liolaemid Lizards (Squamata: Liolaemus kingii clade)

In this study we detangled the evolutionary history of the Patagonian lizard clade Liolaemus kingii, coupling dense geographic sampling and novel computational analytical approaches. We analyzed nuclear and mitochondrial data (restriction site-associated DNA sequencing and cytochrome b) to hypothesize and evaluate species limits, phylogenetic relationships, and demographic histories. We complemented these analyses with posterior predictive simulations to assess the fit of the genomic data to the multispecies coalescent model. We also employed a novel approach to time-calibrate a phylogenetic network. Our results show several instances of mito-nuclear discordance and consistent support for a reticulated history, supporting the view that the complex evolutionary history of the kingii clade is characterized by extensive gene flow and rapid diversification events. We discuss our findings in the contexts of the "gray zone" of speciation, phylogeographic patterns in the Patagonian region, and taxonomic outcomes. [Model adequacy; multispecies coalescent; multispecies network coalescent; phylogenomics; species delimitation.].

Integrative species delimitation in the common ophiuroid Ophiothrix angulata (Echinodermata: Ophiuroidea): insights from COI, ITS2, arm coloration, and geometric morphometrics

Ophiothrix angulata (Say, 1825) is one of the most common and well-known ophiuroids in the Western Atlantic, with a wide geographic and bathymetric range. The taxonomy of this species has been controversial for a century because of its high morphological variability. Here we integrate information from DNA sequence data, color patterns, and geometric morphometrics to assess species delimitation and geographic differentiation in O. angulata. We found three deeply divergent mtDNA-COI clades (K2P 17.0-27.9%). ITS2 nuclear gene and geometric morphometrics of dorsal and ventral arm plates differentiate one of these lineages, as do integrative species delineation analyses, making this a confirmed candidate species.

Picks in the Fabric of a Polyploidy Complex: Integrative Species Delimitation in the Tetraploid Leucanthemum Mill. (Compositae, Anthemideae) Representatives

Based on the results of a preceding species-delimitation analysis for the diploid representatives of the genus Leucanthemum (Compositae, Anthemideae), the present study aims at the elaboration of a specific and subspecific taxonomic treatment of the tetraploid members of the genus. Following an integrative taxonomic approach, species-level decisions on eight predefined morphotaxon hypotheses were based on genetic/genealogical, morphological, ecological, and geographical differentiation patterns. ddRADseq fingerprinting and SNP-based clustering revealed genetic integrity for six of the eight morphotaxa, with no clear differentiation patterns observed between the widespread L. ircutianum subsp. ircutianum and the N Spanish (Cordillera Cantábrica) L. cantabricum and the S French L. delarbrei subsp. delabrei (northern Massif Central) and L. meridionale (western Massif Central). The inclusion of differentiation patterns in morphological (leaf dissection and shape), ecological (climatological and edaphic niches), and geographical respects (pair-wise tests of sympatry vs. allopatry) together with the application of a procedural protocol for species-rank decisions (the 'Wettstein tesseract') led to the proposal of an acknowledgement of the eight predefined morphotaxon hypotheses as six species (two of them with two subspecies). Nomenclatural consequences following from these results are drawn and lead to the following new combinations: Leucanthemum delarbrei subsp. meridionale (Legrand) Oberpr., T.Ott & Vogt, comb. nov. and Leucanthemum ruscinonense (Jeanb. & Timb.-Lagr.) Oberpr., T.Ott & Vogt, comb. et stat. nov.

Delimiting species in the taxonomically challenging orchid section Pseudophrys: Bayesian analyses of genetic and phenotypic data

Accurate species delimitation is critical for biodiversity conservation. Integrative taxonomy has been advocated for a long time, yet tools allowing true integration of genetic and phenotypic data have been developed quite recently and applied to few models, especially in plants. In this study, we investigated species boundaries within a group of twelve Pseudophrys taxa from France by analyzing genetic, morphometric and chemical (i.e., floral scents) data in a Bayesian framework using the program integrated Bayesian Phylogenetics and Phylogeography (iBPP). We found that these twelve taxa were merged into four species when only genetic data were used, while most formally described species were recognized as such when only phenotypic (either morphometric or chemical) data were used. The result of the iBPP analysis performed on both genetic and phenotypic data supports the proposal to merge Ophrys bilunulata and O. marmorata on the one hand, and O. funerea and O. zonata on the other hand. Our results show that phenotypic data are particularly informative in the section Pseudophrys and that their integration in a model-based method significantly improves the accuracy of species delimitation. We are convinced that the integrative taxonomic approach proposed in this study holds great promise to conduct taxonomic revisions in other orchid groups.

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.

The Warps and Wefts of a Polyploidy Complex: Integrative Species Delimitation of the Diploid Leucanthemum (Compositae, Anthemideae) Representatives

Species delimitation—owing to the paramount role of the species rank in evolutionary, ecological, and nature conservation studies—is an essential contribution of taxonomy to biodiversity research. In an ‘integrative taxonomy’ approach to species delimitation on the diploid level, we searched for evolutionary significant units (the warps and wefts) that gave rise to the polyploid complex of European ox-eye daisies (Leucanthemum; Compositae-Anthemideae). Species discovery and validation methods based on genetic, ecological, geographical, and morphometric datasets were applied to test the currently accepted diploid morpho-species, i.e., morphologically delimited species, in Leucanthemum. Novel approaches were taken in the analyses of RADseq data (consensus clustering), morphometrics of reconstructed leaf silhouettes from digitized herbarium specimens, and quantification of species-distribution overlaps. We show that 17 of the 20 Leucanthemum morpho-species are supported by genetic evidence. The taxonomic rank of the remaining three morpho-species was resolved by combining genealogic, ecologic, geographic, and morphologic data in the framework of von Wettstein’s morpho-geographical species concept. We herewith provide a methodological pipeline for the species delimitation in an ‘integrative taxonomy’ fashion using sources of evidence from genealogical, morphological, ecological, and geographical data in the philosophy of De Queiroz’s “Unified Species Concept”.

Coalescent-based species delimitation meets deep learning: Insights from a highly fragmented cactus system

Delimiting species boundaries is a major goal in evolutionary biology. An increasing volume of literature has focused on the challenges of investigating cryptic diversity within complex evolutionary scenarios of speciation, including gene flow and demographic fluctuations. New methods based on model selection, such as approximate Bayesian computation, approximate likelihoods, and machine learning are promising tools arising in this field. Here, we introduce a framework for species delimitation using the multispecies coalescent model coupled with a deep learning algorithm based on convolutional neural networks (CNNs). We compared this strategy with a similar ABC approach. We applied both methods to test species boundary hypotheses based on current and previous taxonomic delimitations as well as genetic data (sequences from 41 loci) in Pilosocereus aurisetus, a cactus species complex with a sky-island distribution and taxonomic uncertainty. To validate our method, we also applied the same strategy on data from widely accepted species from the genus Drosophila. The results show that our CNN approach has a high capacity to distinguish among the simulated species delimitation scenarios, with higher accuracy than ABC. For the cactus data set, a splitter hypothesis without gene flow showed the highest probability in both CNN and ABC approaches, a result agreeing with previous taxonomic classifications and in line with the sky-island distribution and low dispersal of P. aurisetus. Our results highlight the cryptic diversity within the P. aurisetus complex and show that CNNs are a promising approach for distinguishing complex evolutionary histories, even outperforming the accuracy of other model-based approaches such as ABC.

A paradigm for the recognition of cryptic trematode species in tropical Indo-west Pacific fishes: the problematic genus Preptetos (Trematoda: Lepocreadiidae)

Molecular data have transformed approaches to trematode taxonomy by providing objective evidence for the delineation of species. However, although the data are objective, the interpretation of these data regarding species boundaries is subjective, especially when different markers conflict. Conserved markers can lead to an underestimation of richness and those used for finer species delineation have the capacity to inflate species recognition, perhaps unrealistically. Here we examine molecular and morphological evidence for species recognition in an especially confusing system, the lepocreadiid genus Preptetos Pritchard, 1960 in acanthuriform fishes of the tropical Indo-west Pacific. We consider species boundaries within this genus based on combined data (ITS2 and 28S rDNA; cox1 mtDNA and morphometrics) for substantial new collections. Delineation of species using only morphological data suggest fewer species than analysis of the sequence data; the latter suggests the presence of potential cryptic species and analysis of different markers suggests the presence of differing numbers of species. We conclude that an integrative interpretation creates the most satisfying taxonomic hypothesis. In the light of the new data, we have chosen and propose a model of trematode species recognition that demands reciprocal monophyly in the most discriminating available molecular marker plus distinction in morphology or host distribution. By invoking these criteria, we distinguish eight species in our new tropical Indo-west Pacific collections. Six of these are new (Preptetos allocaballeroi n. sp., Preptetos paracaballeroi n. sp., Preptetos pearsoni n. sp., Preptetos prudhoei n. sp., Preptetos quandamooka n. sp. and Preptetos zebravaranus n. sp.) and we continue to recognise Preptetos cannoni Barker, Bray & Cribb, 1993 and Preptetos laguncula Bray and Cribb, 1996. Notably; two of the new species, P. allocaballeroi n. sp. and P. paracaballeroi n. sp., are morphologically cryptic relative to each other. Our criteria lead us to recognise, as species, populations with unvarying morphology and similar host relationships but which may have a complex population structure over their range. In our view, this paradigm has the capacity to render tractable the interpretation of the species status of the huge trematode fauna of the tropical Indo-west Pacific.

Plotting for change: an analytical framework to aid decisions on which lineages are candidate species in phylogenomic species discovery

A recent study argued that coalescent-based models of species delimitation mostly delineate population structure, not species, and called for the validation of candidate species using biological information additional to the genetic information, such as phenotypic or ecological data. Here, we introduce a framework to interrogate genomic datasets and coalescent-based species trees for the presence of candidate species in situations where additional biological data are unavailable, unobtainable or uninformative. For de novo genomic studies of species boundaries, we propose six steps: (1) visualize genetic affinities among individuals to identify both discrete and admixed genetic groups from first principles and to hold aside individuals involved in contemporary admixture for independent consideration; (2) apply phylogenetic techniques to identify lineages; (3) assess diagnosability of those lineages as potential candidate species; (4) interpret the diagnosable lineages in a geographical context (sympatry, parapatry, allopatry); (5) assess significance of difference or trends in the context of sampling intensity; and (6) adopt a holistic approach to available evidence to inform decisions on species status in the difficult cases of allopatry. We adopt this approach to distinguish candidate species from within-species lineages for a widespread species complex of Australian freshwater fishes (Retropinna spp.). Our framework addresses two cornerstone issues in systematics that are often not discussed explicitly in genomic species discovery: diagnosability and how to determine it, and what criteria should be used to decide whether diagnosable lineages are conspecific or represent different species.

A new species of South America marsh rat (Holochilus, Cricetidae) from northeastern Brazil

Marsh rats of the genus Holochilus are broadly distributed and inhabit several distinct environments throughout South America. As an enigmatic group with a somewhat uncertain taxonomy, the composition and names of taxa have shifted throughout history, including the boundaries of Holochilus sciureus, a species formerly distributed in the lowlands of the Guianas, Peru, Bolivia, and northern, central, and northeastern of Brazil. Based on a combination of morphological and morphometric traits, and genomic sequences, we redefined species boundaries and split this wide concept of H. sciureus into three separated species, including a redefined H. sciureus, the newly erected from synonymy H. nanus, and an unnamed taxon. We describe this unnamed species, and provide emended diagnoses for the two redefined species. The newly named taxon inhabits the northeastern part of Brazil and differs from the other congeners by a unique combination of phenotypic and genomic characters.

Type specimen sequencing, multilocus analyses, and species delimitation methods recognize the cosmopolitan Corallina berteroi and establish the northern Japanese C. yendoi sp. nov. (Corallinaceae, Rhodophyta)

A partial rbcL sequence of the lectotype specimen of Corallina berteroi shows that it is the earliest available name for C. ferreyrae. Multilocus species delimitation analyses (ABGD, SPN, GMYC, bPTP, and BPP) using independent or concatenated COI, psbA, and rbcL sequences recognized one, two, or three species in this complex, but only with weak support for each species hypothesis. Conservatively, we recognize a single worldwide species in this complex of what appears to be multiple, evolving populations. Included in this species, besides C. ferreyrae, are C. caespitosa, the morphologically distinct C. melobesioides, and, based on a partial rbcL sequence of the holotype specimen, C. pinnatifolia. Corallina berteroi, not C. officinalis, is the cosmopolitan temperate species found thus far in the NE Atlantic, Mediterranean Sea, warm temperate NW Atlantic and NE Pacific, cold temperate SW Atlantic (Falkland Islands), cold and warm temperate SE Pacific, NW Pacific and southern Australia. Also proposed is C. yendoi sp. nov. from Hokkaido, Japan, which was recognized as distinct by 10 of the 13 species discrimination analyses, including the multilocus BPP.

Identification of Species by Combining Molecular and Morphological Data Using Convolutional Neural Networks

Integrative taxonomy is central to modern taxonomy and systematic biology, including behaviour, niche preference, distribution, morphological analysis and DNA barcoding. However, decades of use demonstrate that these methods can face challenges when used in isolation, for instance, potential misidentifications due to phenotypic plasticity for morphological methods, and incorrect identifications because of introgression, incomplete lineage sorting and horizontal gene transfer for DNA barcoding. Although researchers have advocated the use of integrative taxonomy, few detailed algorithms have been proposed. Here, we develop a convolutional neural network method (morphology-molecule network (MMNet)) that integrates morphological and molecular data for species identification. The newly proposed method (MMNet) worked better than four currently-available alternative methods when tested with 10 independent datasets representing varying genetic diversity from different taxa. High accuracies were achieved for all groups, including beetles (98.1% of 123 species), butterflies (98.8% of 24 species), fishes (96.3% of 214 species) and moths (96.4% of 150 total species). Further, MMNet demonstrated a high degree of accuracy (>98%) in four datasets including closely related species from the same genus. The average accuracy of two modest sub-genomic (single nucleotide polymorphism) datasets, comprising eight putative subspecies respectively, is 90%. Additional tests show that the success rate of species identification under this method most strongly depends on the amount of training data, and is robust to sequence length and image size. Analyses on the contribution of different data types (image versus gene) indicate that both morphological and genetic data are important to the model, and that genetic data contribute slightly more. The approaches developed here serve as a foundation for the future integration of multi-modal information for integrative taxonomy, such as image, audio, video, 3D scanning and biosensor data, to characterize organisms more comprehensively as a basis for improved investigation, monitoring and conservation of biodiversity.

Phylogenomic Assessment of Biodiversity Using a Reference-Based Taxonomy: An Example With Horned Lizards (Phrynosoma)

Phylogenomic investigations of biodiversity facilitate the detection of fine-scale population genetic structure and the demographic histories of species and populations. However, determining whether or not the genetic divergence measured among populations reflects species-level differentiation remains a central challenge in species delimitation. One potential solution is to compare genetic divergence between putative new species with other closely related species, sometimes referred to as a reference-based taxonomy. To be described as a new species, a population should be at least as divergent as other species. Here, we develop a reference-based taxonomy for Horned Lizards (Phrynosoma; 17 species) using phylogenomic data (ddRADseq data) to provide a framework for delimiting species in the Greater Short-horned Lizard species complex (P. hernandesi). Previous species delimitation studies of this species complex have produced conflicting results, with morphological data suggesting that P. hernandesi consists of five species, whereas mitochondrial DNA support anywhere from 1 to 10 + species. To help address this conflict, we first estimated a time-calibrated species tree for P. hernandesi and close relatives using SNP data. These results support the paraphyly of P. hernandesi; we recommend the recognition of two species to promote a taxonomy that is consistent with species monophyly. There is strong evidence for three populations within P. hernandesi, and demographic modeling and admixture analyses suggest that these populations are not reproductively isolated, which is consistent with previous morphological analyses that suggest hybridization could be common. Finally, we characterize the population-species boundary by quantifying levels of genetic divergence for all 18 Phrynosoma species. Genetic divergence measures for western and southern populations of P. hernandesi failed to exceed those of other Phrynosoma species, but the relatively small population size estimated for the northern population causes it to appear as a relatively divergent species. These comparisons underscore the difficulties associated with putting a reference-based approach to species delimitation into practice. Nevertheless, the reference-based approach offers a promising framework for the consistent assessment of biodiversity within clades of organisms with similar life histories and ecological traits.

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

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

Phylogenomic species delimitation dramatically reduces species diversity in an Antarctic adaptive radiation

Application of genetic data to species delimitation often builds confidence in delimitations previously hypothesized using morphological, ecological, and geographic data and frequently yields recognition of previously-undescribed cryptic diversity. However, a recent critique of genomic data-based species delimitation approaches is that they have the potential to conflate population structure with species diversity, resulting in taxonomic oversplitting. The need for an integrative approach to species delimitation, in which molecular, morphological, ecological, and geographic lines of evidence are evaluated together, is becoming increasingly apparent. Here, we integrate phylogenetic, population genetic, and coalescent analyses of genome-wide sequence data with investigation of variation in multiple morphological traits to delimit species within the Antarctic barbeled plunderfishes (Artedidraconidae: Pogonophryne). Pogonophryne currently comprises 29 valid species, most of which are distinguished solely by variation in ornamentation of the mental barbel that projects from the lower jaw, a structure previously shown to vary widely within a single species. However, our genomic and phenotypic analyses result in a dramatic reduction in the number of distinct species recognized within the clade, providing evidence to support the recognition of no more than six species. We propose to synonymize 24 of the currently recognized species with five species of Pogonophryne. We find genomic and phenotypic evidence for a new species of Pogonophryne from specimens collected in the Ross Sea. Our findings represent a rare example in which application of molecular data provides evidence of taxonomic oversplitting on the basis of morphology, clearly demonstrating the utility of an integrative species delimitation framework.

A taxonomist's nightmare - Cryptic diversity in Caribbean intertidal arthropods (Arachnida, Acari, Oribatida)

There has been a long controversy about what defines a species and how to delimitate them which resulted in the existence of more than two dozen different species concepts. Recent research on so-called "cryptic species" heated up this debate as some scientists argue that these cryptic species are only a result of incompatible species concepts. While this may be true, we should keep in mind that all concepts are nothing more than human constructs and that the phenomenon of high phenotypic similarity despite reproductive isolation is real. To investigate and understand this phenomenon it is important to classify and name cryptic species as it allows to communicate them with other fields of science that use Linnaean binomials. To provide a common framework for the description of cryptic species, we propose a possible protocol of how to formally name and describe these taxa in practice. The most important point of this protocol is to explain which species concept was used to delimitate the cryptic taxon. As a model, we present the case of the allegedly widespread Caribbean intertidal mite Thalassozetes barbara, which in fact consists of seven phenotypically very similar but genetically distinct species. All species are island or short-range endemics with poor dispersal abilities that have evolved in geographic isolation. Stabilizing selection caused by the extreme conditions of the intertidal environment is suggested to be responsible for the morphological stasis of this cryptic species complex.

Phylogeography, historical demography and systematics of the world's smallest pythons (Pythonidae, Antaresia)

Advances from empirical studies in phylogeography, systematics and species delimitation highlight the importance of integrative approaches for quantifying taxonomic diversity. Genomic data have greatly improved our ability to discern both systematic diversity and evolutionary history. Here we combine analyses of mitochondrial DNA sequences, thousands of genome-wide SNPs and linear and geometric morphometrics on Antaresia, a clade of four currently recognised dwarf pythons from Australia and New Guinea (Antaresia childreni, A. stimsoni, A. maculosa and A. perthensis). Our integrative analyses of phylogenetics, population structure, species delimitation, historical demography and morphometrics revealed that the true evolutionary diversity is not well reflected in the current appraisal of the diversity of the group. We find that Antaresia childreni and A. stimsoni comprise a widespread network of populations connected by gene flow and without evidence of species-level divergence among them. However, A. maculosa shows considerable genetic structuring which leads us to recognise two subspecies in northeastern Australia and a new species in Torres Strait and New Guinea. These two contrasting cases of over and under estimation of diversity, respectively, illustrate the power of thorough integrative approaches into understanding evolution of biodiversity. Furthermore, our analyses of historical demographic patterns highlight the importance of the Kimberley, Pilbara and Cape York as origins of biodiversity in Australia.

The number of species of degus (genus Octodon) is currently underestimated: An appraisal of species limits and their phylogenetic relationships (Rodentia: Hystricomorpha: Octodontidae)

As currently understood, the genus Octodon contains five species degus, lunatus, bridgesii, pacificus, and ricardojeda. Previous phylogenetic studies suggest that genus specific diversity is underestimated. In order to evaluate the taxonomic diversity of Octodon, we implemented unilocus (cytochrome-b) and multilocus (cytochrome-b + 4 nuclear genes) species delimitation methods. Octodon degus was recovered as a sister of the other species of the genus. The unilocus bGMYC and mPTP methods, based on cytochrome-b sequences, delimits 11 and 7 candidate species respectively, and both methods fail to recognize O. pacificus from O. ricardojeda. Results of the multilocus analysis (BPP) vary as a function of the dataset used. When the five genes are used 11 species are delimited, while eight species are delimited when only the nuclear genes are used. Octodon bridgesii is shown as comprising at least two species (one on the Pacific coast and the typical form found on the Andean slopes), while O. ricardojeda may comprise two species (one on the Chilean side of the Andes and the other in Argentina). Likewise, both multilocus matrices recover O. pacificus as a distinct species. This shows that species diversity of Octodon is underestimated. Remarkably, many of the delimited species based on genetic data are morphologically differentiated in cranio-dental characteristics. However, a pair of species has not achieved morphological differentiation, being cryptic species. Finally, the incongruence between mitochondrial and nuclear phylogenies suggests that processes such as incomplete lineage sorting and/or introgression have been present during the radiation of the genus.

Investigating Morphological Complexes Using Informational Dissonance and Bayes Factors: A Case Study in Corbiculate Bees

It is widely recognized that different regions of a genome often have different evolutionary histories and that ignoring this variation when estimating phylogenies can be misleading. However, the extent to which this is also true for morphological data is still largely unknown. Discordance among morphological traits might plausibly arise due to either variable convergent selection pressures or else phenomena such as hemiplasy. Here, we investigate patterns of discordance among 282 morphological characters, which we scored for 50 bee species particularly targeting corbiculate bees, a group that includes the well-known eusocial honeybees and bumblebees. As a starting point for selecting the most meaningful partitions in the data, we grouped characters as morphological modules, highly integrated trait complexes that as a result of developmental constraints or coordinated selection we expect to share an evolutionary history and trajectory. In order to assess conflict and coherence across and within these morphological modules, we used recently developed approaches for computing Bayesian phylogenetic information allied with model comparisons using Bayes factors. We found that despite considerable conflict among morphological complexes, accounting for among-character and among-partition rate variation with individual gamma distributions, rate multipliers, and linked branch lengths can lead to coherent phylogenetic inference using morphological data. We suggest that evaluating information content and dissonance among partitions is a useful step in estimating phylogenies from morphological data, just as it is with molecular data. Furthermore, we argue that adopting emerging approaches for investigating dissonance in genomic datasets may provide new insights into the integration and evolution of anatomical complexes. [Apidae; entropy; morphological modules; phenotypic integration; phylogenetic information.].

SODA: Multi-locus species delimitation using quartet frequencies

MOTIVATION

Species delimitation, the process of deciding how to group a set of organisms into units called species, is one of the most challenging problems in evolutionary computational biology. While many methods exist for species delimitation, most based on the coalescent theory, few are scalable to very large datasets, and methods that scale tend to be not accurate. Species delimitation is closely related to species tree inference from discordant gene trees, a problem that has enjoyed rapid advances in recent years.

RESULTS

In this paper, we build on the accuracy and scalability of recent quartet-based methods for species tree estimation and propose a new method called SODA for species delimitation. SODA relies heavily on a recently developed method for testing zero branch length in species trees. In extensive simulations, we show that SODA can easily scale to very large datasets while maintaining high accuracy.

AVAILABILITY

The code and data presented here are available on https://github.com/maryamrabiee/SODA.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Minimalist revision and description of 403 new species in 11 subfamilies of Costa Rican braconid parasitoid wasps, including host records for 219 species

Three new genera are described: Michener (Proteropinae), Bioalfa (Rogadinae), and Hermosomastax (Rogadinae). Keys are given for the New World genera of the following braconid subfamilies: Agathidinae, Braconinae, Cheloninae, Homolobinae, Hormiinae, Ichneutinae, Macrocentrinae, Orgilinae, Proteropinae, Rhysipolinae, and Rogadinae. In these subfamilies 416 species are described or redescribed. Most of the species have been reared and all but 13 are new to science. A consensus sequence of the COI barcodes possessed by each species is employed to diagnose the species, and this approach is justified in the introduction. Most descriptions consist of a lateral or dorsal image of the holotype, a diagnostic COI consensus barcode, the Barcode Index Number (BIN) code with a link to the Barcode of Life Database (BOLD), and the holotype specimen information required by the International Code of Zoological Nomenclature. The following species are treated and those lacking authorship are newly described here with authorship attributable to Sharkey except for the new species of Macrocentrinae which are by Sharkey & van Achterberg: AGATHIDINAE: Aerophiluspaulmarshi, Mesocoelusdavidsmithi, Neothlipsisbobkulai, Plesiocoelusvanachterbergi, Pneumagathiserythrogastra (Cameron, 1905), Therophilusbobwhartoni, T.donaldquickei, T.gracewoodae, T.maetoi, T.montywoodi, T.penteadodiasae, Zacremnopsbrianbrowni, Z.coatlicue Sharkey, 1990, Zacremnopscressoni (Cameron, 1887), Z.ekchuah Sharkey, 1990, Z.josefernandezi, Zelomorphasarahmeierottoae. BRACONINAE: Braconalejandromarini, B.alejandromasisi, B.alexamasisae, B.andresmarini, B.andrewwalshi, B.anniapicadoae, B.anniemoriceae, B.barryhammeli, B.bernardoespinozai, B.carlossanabriai, B.chanchini, B.christophervallei, B.erasmocoronadoi, B.eugeniephillipsae, B.federicomatarritai, B.frankjoycei, B.gerardovegai, B.germanvegai, B.isidrochaconi, B.jimlewisi, B.josejaramilloi, B.juanjoseoviedoi, B.juliodiazi, B.luzmariaromeroae, B.manuelzumbadoi, B.marialuisariasae, B.mariamartachavarriae, B.mariorivasi, B.melissaespinozae, B.nelsonzamorai, B.nicklaphami, B.ninamasisae, B.oliverwalshi, B.paulamarinae, B.rafamoralesi, B.robertofernandezi, B.rogerblancoi, B.ronaldzunigai, B.sigifredomarini, B.tihisiaboshartae, B.wilberthbrizuelai, Digonogastramontylloydi, D.montywoodi, D.motohasegawai, D.natwheelwrighti, D.nickgrishini. CHELONINAE: Adeliusadrianguadamuzi, A.gauldi Shimbori & Shaw, 2019, A.janzeni Shimbori & Shaw, 2019, Ascogastergloriasihezarae, A.grettelvegae, A.guillermopereirai, A.gustavoecheverrii, A.katyvandusenae, A.luisdiegogomezi, Chelonusalejandrozaldivari, C.gustavogutierrezi, C.gustavoinduni, C.harryramirezi, C.hartmanguidoi, C.hazelcambroneroae, C.iangauldi, C.isidrochaconi, C.janecheverriae, C.jeffmilleri, C.jennyphillipsae, C.jeremydewaardi, C.jessiehillae, C.jesusugaldei, C.jimlewisi, C.jimmilleri, C.jimwhitfieldi, C.johanvalerioi, C.johnburnsi, C.johnnoyesi, C.jorgebaltodanoi, C.jorgehernandezi, C.josealfredohernandezi, C.josefernandeztrianai, C.josehernandezcortesi, C.josemanuelperezi, C.josephinerodriguezae, C.juanmatai, C.junkoshimurae, C.kateperezae, C.luciariosae, C.luzmariaromeroae, C.manuelpereirai, C.manuelzumbadoi, C.marianopereirai, C.maribellealvarezae, C.markmetzi, C.markshawi, C.martajimenezae, C.mayrabonillae, C.meganmiltonae, C.melaniamunozae, C.michaelstroudi, C.michellevanderbankae, C.mingfangi, C.minorcarmonai, C.monikaspringerae, C.moniquegilbertae, C.motohasegawai, C.nataliaivanovae, C.nelsonzamorai, C.normwoodleyi, C.osvaldoespinozai, C.pamelacastilloae, C.paulgoldsteini, C.paulhansoni, C.paulheberti, C.petronariosae, C.ramyamanjunathae, C.randallgarciai, C.rebeccakittelae, C.robertoespinozai, C.robertofernandezi, C.rocioecheverriae, C.rodrigogamezi, C.ronaldzunigai, C.rosibelelizondoae, C.rostermoragai, C.ruthfrancoae, C.scottmilleri, C.scottshawi, C.sergioriosi, C.sigifredomarini, C.stevearonsoni, C.stevestroudi, C.sujeevanratnasinghami, C.sureshnaiki, C.torbjornekremi, C.yeimycedenoae, Leptodrepanaalexisae, L.erasmocoronadoi, L.felipechavarriai, L.freddyquesadai, L.gilbertfuentesi, L.manuelriosi, Phanerotomaalmasolisae, P.alvaroherrerai, P.anacordobae, P.anamariamongeae, P.andydeansi, P.angelagonzalezae, P.angelsolisi, P.barryhammeli, P.bernardoespinozai, P.calixtomoragai, P.carolinacanoae, P.christerhanssoni, P.christhompsoni, P.davesmithi, P.davidduthiei, P.dirksteinkei, P.donquickei, P.duniagarciae, P.duvalierbricenoi, P.eddysanchezi, P.eldarayae, P.eliethcantillanoae, P.jenopappi, Pseudophanerotomaalanflemingi, Ps.albanjimenezi, Ps.alejandromarini, Ps.alexsmithi, Ps.allisonbrownae, Ps.bobrobbinsi. HOMOLOBINAE: Exasticolusjennyphillipsae, E.randallgarciai, E.robertofernandezi, E.sigifredomarini, E.tomlewinsoni. HORMIINAE: Hormiusanamariamongeae, H.angelsolisi, H.anniapicadoae, H.arthurchapmani, H.barryhammeli, H.carmenretanae, H.carloswalkeri, H.cesarsuarezi, H.danbrooksi, H.eddysanchezi, H.erikframstadi, H.georgedavisi, H.grettelvegae, H.gustavoinduni, H.hartmanguidoi, H.hectoraritai, H.hesiquiobenitezi, H.irenecanasae, H.isidrochaconi, H.jaygallegosi, H.jimbeachi, H.jimlewisi, H.joelcracrafti, H.johanvalerioi, H.johnburleyi, H.joncoddingtoni, H.jorgecarvajali, H.juanmatai, H.manuelzumbadoi, H.mercedesfosterae, H.modonnellyae, H.nelsonzamorai, H.pamelacastilloae, H.raycypessi, H.ritacolwellae, H.robcolwelli, H.rogerblancosegurai, H.ronaldzunigai, H.russchapmani, H.virginiaferrisae, H.warrenbrighami, H.willsflowersi. ICHNEUTINAE: Oligoneuruskriskrishtalkai, O.jorgejimenezi, Paroligoneuruselainehoaglandae, P.julianhumphriesi, P.mikeiviei. MACROCENTRINAE: Austrozelejorgecampabadali, A.jorgesoberoni, Dolichozelegravitarsis (Muesebeck, 1938), D.josefernandeztrianai, D.josephinerodriguezae, Hymenochaoniakalevikulli, H.kateperezae, H.katherinebaillieae, H.katherineellisonae, H.katyvandusenae, H.kazumifukunagae, H.keithlangdoni, H.keithwillmotti, H.kenjinishidai, H.kimberleysheldonae, H.krisnorvigae, H.lilianamadrigalae, H.lizlangleyae, Macrocentrusfredsingeri, M.geoffbarnardi, M.gregburtoni, M.gretchendailyae, M.grettelvegae, M.gustavogutierrezi, M.hannahjamesae, M.harisridhari, M.hillaryrosnerae, M.hiroshikidonoi, M.iangauldi, M.jennyphillipsae, M.jesseausubeli, M.jessemaysharkae, M.jimwhitfieldi, M.johnbrowni, M.johnburnsi, M.jonathanfranzeni, M.jonathanrosenbergi, M.jorgebaltodanoi, M.lucianocapelli. ORGILINAE: Orgilusamyrossmanae, O.carrolyoonae, O.christhompsoni, O.christinemcmahonae, O.dianalipscombae, O.ebbenielsoni, O.elizabethpennisiae, O.evertlindquisti, O.genestoermeri, O.jamesriegeri, O.jeanmillerae, O.jeffmilleri, O.jerrypowelli, O.jimtiedjei, O.johnlundbergi, O.johnpipolyi, O.jorgellorentei, O.larryspearsi, O.marlinricei, O.mellissaespinozae, O.mikesmithi, O.normplatnicki, O.peterrauchi, O.richardprimacki, O.sandraberriosae, O.sarahmirandae, O.scottmilleri, O.scottmorii, Stantoniabillalleni, S.brookejarvisae, S.donwilsoni, S.erikabjorstromae, S.garywolfi, S.henrikekmani, S.luismirandai, S.miriamzunzae, S.quentinwheeleri, S.robinkazmierae, S.ruthtifferae. PROTEROPINAE: Hebichneutestricolor Sharkey & Wharton, 1994, Proteropsiangauldi, P.vickifunkae, Michenercharlesi. RHYSIPOLINAE: Pseudorhysipolisluisfonsecai, P. mailyngonzalezaeRhysipolisjulioquirosi. ROGADINAE: Aleiodesadrianaradulovae, A.adrianforsythi, A.agnespeelleae, A.alaneaglei, A.alanflemingi, A.alanhalevii, A.alejandromasisi, A.alessandracallejae, A.alexsmithi, A.alfonsopescadori, A.alisundermieri, A.almasolisae, A.alvarougaldei, A.alvaroumanai, A.angelsolisi, A.annhowdenae, A.bobandersoni, A.carolinagodoyae, A.charlieobrieni, A.davefurthi, A.donwhiteheadi, A.doylemckeyi, A.frankhovorei, A.henryhowdeni, A.inga Shimbori & Shaw, 2020, A.johnchemsaki, A.johnkingsolveri, A.gonodontovorus Shimbori & Shaw, 2020, A.manuelzumbadoi, A.mayrabonillae, A.michelledsouzae, A.mikeiviei, A.normwoodleyi, A.pammitchellae, A.pauljohnsoni, A.rosewarnerae, A.steveashei, A.terryerwini, A.willsflowersi, Bioalfapedroleoni, B.alvarougaldei, B.rodrigogamezi, Choreborogasandydeansi, C.eladiocastroi, C.felipechavarriai, C.frankjoycei, Clinocentrusandywarreni, Cl.angelsolisi, Cystomastaxalexhausmanni, Cy.angelagonzalezae, Cy.ayaigarashiae, Hermosomastaxclavifemorus Quicke sp. nov., Heterogamusdonstonei, Pseudoyeliconesbernsweeneyi, Stiropiusbencrairi, S.berndkerni, S.edgargutierrezi, S.edwilsoni, S.ehakernae, Triraphisbillfreelandi, T.billmclarneyi, T.billripplei, T.bobandersoni, T.bobrobbinsi, T.bradzlotnicki, T.brianbrowni, T.brianlaueri, T.briannestjacquesae, T.camilocamargoi, T.carlosherrerai, T.carolinepalmerae, T.charlesmorrisi, T.chigiybinellae, T.christerhanssoni, T.christhompsoni, T.conniebarlowae, T.craigsimonsi, T.defectus Valerio, 2015, T.danielhubi, T.davidduthiei, T.davidwahli, T.federicomatarritai, T.ferrisjabri, T.mariobozai, T.martindohrni, T.matssegnestami, T.mehrdadhajibabaei, T.ollieflinti, T.tildalauerae, Yeliconesdirksteinkei, Y.markmetzi, Y.monserrathvargasae, Y.tricolor Quicke, 1996. Y.woldai Quicke, 1996.

The following new combinations are proposed: Neothlipsissmithi (Ashmead), new combination for Microdussmithi Ashmead, 1894; Neothlipsispygmaeus (Enderlein), new combination for Microduspygmaeus Enderlein, 1920; Neothlipsisunicinctus (Ashmead), new combination for Microdusunicinctus Ashmead, 1894; Therophilusanomalus (Bortoni and Penteado-Dias) new combination for Plesiocoelusanomalus Bortoni and Penteado-Dias, 2015; Aerophilusareolatus (Bortoni and Penteado-Dias) new combination for Plesiocoelusareolatus Bortoni and Penteado-Dias, 2015; Pneumagathiserythrogastra (Cameron) new combination for Agathiserythrogastra Cameron, 1905. Dolichozelecitreitarsis (Enderlein), new combination for Paniscozelecitreitarsis Enderlein, 1920. Dolichozelefuscivertex (Enderlein) new combination for Paniscozelefuscivertex Enderlein, 1920. Finally, Bassusbrooksi Sharkey, 1998 is synonymized with Agathiserythrogastra Cameron, 1905; Paniscozelegriseipes Enderlein, 1920 is synonymized with Dolichozelekoebelei Viereck, 1911; Paniscozelecarinifrons Enderlein, 1920 is synonymized with Dolichozelefuscivertex (Enderlein, 1920); and Paniscozelenigricauda Enderlein,1920 is synonymized with Dolichozelequaestor (Fabricius, 1804). (originally described as Ophionquaestor Fabricius, 1804).

An integrative approach to address species limits in the southernmost members of the Liolaemus kingii group (Squamata: Liolaemini)

Recent conceptual and methodological advances have enabled an increasing number of studies to address the problem of species delimitation in a comprehensive manner. This is of particular interest in cases of species whose divergence times are recent and/or effective population sizes are large, where the conclusions obtained from a single source of evidence may lead to erroneous estimations of true species numbers or incorrect assignment of individuals to species. Iguanian lizards of the Liolaemus kingii group (13 species) comprise an important component of the endemic fauna of Patagonia. The southernmost species of this group (namely L. baguali, L. escarchadosi, L. sarmientoi, and L. tari) show widely overlapping distributions across southern Patagonia, also, their phylogenetic relationships are ambiguous and species boundaries have not been explicitly tested. Here we use a comprehensive approach to assess species limits through the use of molecular and morphological information (mitochondrial cytb, nuclear sequences collected by ddRADseq, and linear, meristic and landmark-based morphometrics). We found support for the current taxonomy given that the different analyses recognized the nominal species (4 entities), also a candidate species was supported by mitochondrial and morphological data. In addition, we detected signs of admixture between some of the species. Our results indicate that the L. kingii group can serve as a model system in studies of diversification accompanied by hybridization in nature, which in turn might have been promoted by past climatic oscillations and generalist morphologies. We emphasize the importance of using multiple lines of evidence in order to solve evolutionary stories, and minimizing potential erroneous results that may arise when relying on a single source of information.

Influence of niche breadth and position on the historical biogeography of seafaring scincid lizards

Niche breadth and position can influence diversification among closely related species or populations, yet limited empirical data exist concerning the predictability of the outcomes. We explored the effects of these factors on the evolution of the Emoia atrocostata species group, an insular radiation of lizards in the western Pacific Ocean and Indo-Australasia composed of both endemic and widespread species that differ in niche occupancy. We used molecular data and phylogeographical diffusion models to estimate the timing and patterns of range expansion, and ancestral reconstruction methods to infer shifts in ecology. We show evidence of multidirectional spread from a centre of origin in western Micronesia, and that the phyletic diversity of the group is derived from a putative habitat specialist that survives in the littoral zone. This species is composed of paraphyletic lineages that represent stages or possible endpoints in the continuum toward speciation. Several descendant species have transitioned to either strand or interior forest habitat, but only on remote islands with depauperate terrestrial faunas. Our results suggest that the atrocostata group might be in the early phases of a Wilsonian taxon cycle and that the capacity to tolerate salt stress has promoted dispersal and colonization of remote oceanic islands. Divergence itself, however, is largely driven by geographical isolation rather than shifts in ecology.

Repositories for Taxonomic Data: Where We Are and What is Missing

Natural history collections are leading successful large-scale projects of specimen digitization (images, metadata, DNA barcodes), thereby transforming taxonomy into a big data science. Yet, little effort has been directed towards safeguarding and subsequently mobilizing the considerable amount of original data generated during the process of naming 15,000-20,000 species every year. From the perspective of alpha-taxonomists, we provide a review of the properties and diversity of taxonomic data, assess their volume and use, and establish criteria for optimizing data repositories. We surveyed 4113 alpha-taxonomic studies in representative journals for 2002, 2010, and 2018, and found an increasing yet comparatively limited use of molecular data in species diagnosis and description. In 2018, of the 2661 papers published in specialized taxonomic journals, molecular data were widely used in mycology (94%), regularly in vertebrates (53%), but rarely in botany (15%) and entomology (10%). Images play an important role in taxonomic research on all taxa, with photographs used in >80% and drawings in 58% of the surveyed papers. The use of omics (high-throughput) approaches or 3D documentation is still rare. Improved archiving strategies for metabarcoding consensus reads, genome and transcriptome assemblies, and chemical and metabolomic data could help to mobilize the wealth of high-throughput data for alpha-taxonomy. Because long-term-ideally perpetual-data storage is of particular importance for taxonomy, energy footprint reduction via less storage-demanding formats is a priority if their information content suffices for the purpose of taxonomic studies. Whereas taxonomic assignments are quasifacts for most biological disciplines, they remain hypotheses pertaining to evolutionary relatedness of individuals for alpha-taxonomy. For this reason, an improved reuse of taxonomic data, including machine-learning-based species identification and delimitation pipelines, requires a cyberspecimen approach-linking data via unique specimen identifiers, and thereby making them findable, accessible, interoperable, and reusable for taxonomic research. This poses both qualitative challenges to adapt the existing infrastructure of data centers to a specimen-centered concept and quantitative challenges to host and connect an estimated $ \le $2 million images produced per year by alpha-taxonomic studies, plus many millions of images from digitization campaigns. Of the 30,000-40,000 taxonomists globally, many are thought to be nonprofessionals, and capturing the data for online storage and reuse therefore requires low-complexity submission workflows and cost-free repository use. Expert taxonomists are the main stakeholders able to identify and formalize the needs of the discipline; their expertise is needed to implement the envisioned virtual collections of cyberspecimens. [Big data; cyberspecimen; new species; omics; repositories; specimen identifier; taxonomy; taxonomic data.].

Multiple species delimitation approaches applied to the avian lark genus Alaudala

Species delimitation has advanced from a purely phenotypic exercise to a branch of science that integrates multiple sources of data to identify independently evolving lineages that can be treated as species. We here test species limits in the avian Lesser Short-toed Lark Alaudala rufesens-Sand Lark A. raytal complex, which has an intricate taxonomic history, ranging from a single to three recognised species, with different inclusiveness in different treatments. Our integrative taxonomic approach is based on a combination of DNA sequences, plumage, biometrics, songs, song-flights, geographical distributions, habitat, and bioclimatic data, and using various methods including a species delimitation program (STACEY) based on the multispecies coalescent model. We propose that four species should be recognised: Lesser Short-toed Lark A. rufescens (sensu stricto), Heine's Short-toed Lark A. heinei, Asian Short-toed Lark A. cheleensis and Sand Lark A. raytal. There is also some evidence suggesting lineage separation within A. cheleensis and A. raytal, but additional data are required to evaluate this. The species delimitation based on STACEY agrees well with the non-genetic data. Although computer-based species delimitation programs can be useful in identifying independently evolving lineages, we stress that whenever possible, species hypotheses proposed by these programs should be tested by independent, non-genetic data. Our results highlight the difficulty and subjectivity of delimiting lineages and species, especially at early stages in the speciation process.

Phylogeny and species delimitation of the genus Longgenacris and Fruhstorferiola viridifemorata species group (Orthoptera: Acrididae: Melanoplinae) based on molecular evidence

Phylogenetic positions of the genus Longgenacris and one of its members, i.e. L. rufiantennus are controversial. The species boundaries within both of L. rufiantennus+Fruhstorferiola tonkinensis and F. viridifemorata species groups are unclear. In this study, we explored the phylogenetic positions of the genus Longgenacris and the species L. rufiantennus and the relationships among F. viridifemorata group based on the 658-base fragment of the mitochondrial gene cytochrome c oxidase subunit I (COI) barcode and the complete sequences of the internal transcribed spacer regions (ITS1 and ITS2) of the nuclear ribosomal DNA. The phylogenies were reconstructed in maximum likelihood framework using IQ-TREE. K2P distances were used to assess the overlap range between intraspecific variation and interspecific divergence. Phylogenetic species concept and NJ tree, K2P distance, the statistical parsimony network as well as the generalized mixed Yule coalescent model (GMYC) were employed to delimitate the species boundaries in L. rufiantennus+F. tonkinensis and F. viridifemorata species groups. The results demonstrated that the genus Longgenacris should be placed in the subfamily Melanoplinae but not Catantopinae, and L. rufiantennus should be a member of the genus Fruhstorferiola but not Longgenacris. Species boundary delimitation confirmed the presence of oversplitting in L. rufiantennus+F. tonkinensis and F. viridifemorata species groups and suggested that each group should be treated as a single species.

Shining a light on species delimitation in the tree genus Engelhardia Leschenault ex Blume (Juglandaceae)

Enhanced efficacy in species delimitation is critically important in biology given the pending biodiversity crisis under global warming and anthropogenic activity. In particular, delineation of traditional classifications in view of the complexity of species requires an integrative approach to effectively define species boundaries, and this is a major focus of systematic biology. Here, we explored species delimitation of Engelhardia in tropical and subtropical Asia. In total, 716 individuals in 71 populations were genotyped using five chloroplast regions, one nuclear DNA region (nrITS), and 11 nuclear simple sequence repeats (nSSR). Phylogenetic trees were constructed and relationships among species were assessed. Molecular analyses were then combined with 14 morphological characteristics of 720 specimens to further explore the species boundaries of Engelhardia. Integrating phylogenetic and morphological clusters provided well-resolved relationships to delineate seven species. The results suggested that: first, that E. fenzelii, E. roxburghiana, E. hainanensis, E. apoensis, and E. serrata are distinct species; second, E. spicata var. spicata, E. spicata var. aceriflora, E. spicata var. colebrookeana, and E. rigida should be combined under E. spicata and treated as a species complex; third, E. serrata var. cambodica should be raised to species level and named E. villosa. We illuminated that bias thresholds determining the cluster number for delimiting species boundaries were substantially reduced when morphological data were incorporated. Our results urge caution when using the concepts of subspecies and varieties in order to prevent confusion, particularly with respect to species delimitation for tropical and subtropical species. In some cases, re-ranking or combining subspecies and/or varieties may enable more accurate species delimitation.

Phylogenetic and spatial distribution of evolutionary diversification, isolation, and threat in turtles and crocodilians (non-avian archosauromorphs)

Background

The origin of turtles and crocodiles and their easily recognized body forms dates to the Triassic and Jurassic. Despite their long-term success, extant species diversity is low, and endangerment is extremely high compared to other terrestrial vertebrate groups, with ~ 65% of ~ 25 crocodilian and ~ 360 turtle species now threatened by exploitation and habitat loss. Here, we combine available molecular and morphological evidence with statistical and machine learning algorithms to present a phylogenetically informed, comprehensive assessment of diversification, threat status, and evolutionary distinctiveness of all extant species.

Results

In contrast to other terrestrial vertebrates and their own diversity in the fossil record, the recent extant lineages of turtles and crocodilians have not experienced any global mass extinctions or lineage-wide shifts in diversification rate or body-size evolution over time. We predict threat statuses for 114 as-yet unassessed or data-deficient species and identify a concentration of threatened turtles and crocodilians in South and Southeast Asia, western Africa, and the eastern Amazon. We find that unlike other terrestrial vertebrate groups, extinction risk increases with evolutionary distinctiveness: a disproportionate amount of phylogenetic diversity is concentrated in evolutionarily isolated, at-risk taxa, particularly those with small geographic ranges. Our findings highlight the important role of geographic determinants of extinction risk, particularly those resulting from anthropogenic habitat-disturbance, which affect species across body sizes and ecologies.

Conclusions

Extant turtles and crocodilians maintain unique, conserved morphologies which make them globally recognizable. Many species are threatened due to exploitation and global change. We use taxonomically complete, dated molecular phylogenies and various approaches to produce a comprehensive assessment of threat status and evolutionary distinctiveness of both groups. Neither group exhibits significant overall shifts in diversification rate or body-size evolution, or any signature of global mass extinctions in recent, extant lineages. However, the most evolutionarily distinct species tend to be the most threatened, and species richness and extinction risk are centered in areas of high anthropogenic disturbance, particularly South and Southeast Asia. Range size is the strongest predictor of threat, and a disproportionate amount of evolutionary diversity is at risk of imminent extinction.

Genomic nucleotide-based distance analysis for delimiting old world monkey derived herpes simplex virus species

BACKGROUND

Herpes simplex viruses form a genus within the alphaherpesvirus subfamily, with three identified viral species isolated from Old World monkeys (OWM); Macacine alphaherpesvirus 1 (McHV-1; herpes B), Cercopithecine alphaherpesvirus 2 (SA8), and Papiine alphaherpesvirus 2 (PaHV-2; herpes papio). Herpes B is endemic to macaques, while PaHV-2 and SA8 appear endemic to baboons. All three viruses are genetically and antigenically similar, with SA8 and PaHV-2 thought to be avirulent in humans, while herpes B is a biosafety level 4 pathogen. Recently, next-generation sequencing (NGS) has resulted in an increased number of published OWM herpes simplex genomes, allowing an encompassing phylogenetic analysis.

RESULTS

In this study, phylogenetic networks, in conjunction with a genome-based genetic distance cutoff method were used to examine 27 OWM monkey herpes simplex isolates. Genome-based genetic distances were calculated, resulting in distances between lion and pig-tailed simplex viruses themselves, and versus herpes B core strains that were higher than those between PaHV-2 and SA8 (approximately 14 and 10% respectively). The species distance cutoff was determined to be 8.94%, with the method recovering separate species status for PaHV-2 and SA8 and showed that lion and pig-tailed simplex viruses (vs core herpes B strains) were well over the distance species cutoff.

CONCLUSIONS

We propose designating lion and pig-tailed simplex viruses as separate, individual viral species, and that this may be the first identification of viral cryptic species.