Algorithmic improvements to species delimitation and phylogeny estimation under the multispecies coalescent

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.

Ranked Subtree Prune and Regraft

Phylogenetic trees are a mathematical formalisation of evolutionary histories between organisms, species, genes, cancer cells, etc. For many applications, e.g. when analysing virus transmission trees or cancer evolution, (phylogenetic) time trees are of interest, where branch lengths represent times. Computational methods for reconstructing time trees from (typically molecular) sequence data, for example Bayesian phylogenetic inference using Markov Chain Monte Carlo (MCMC) methods, rely on algorithms that sample the treespace. They employ tree rearrangement operations such as [Formula: see text] (Subtree Prune and Regraft) and [Formula: see text] (Nearest Neighbour Interchange) or, in the case of time tree inference, versions of these that take times of internal nodes into account. While the classic [Formula: see text] tree rearrangement is well-studied, its variants for time trees are less understood, limiting comparative analysis for time tree methods. In this paper we consider a modification of the classical [Formula: see text] rearrangement on the space of ranked phylogenetic trees, which are trees equipped with a ranking of all internal nodes. This modification results in two novel treespaces, which we propose to study. We begin this study by discussing algorithmic properties of these treespaces, focusing on those relating to the complexity of computing distances under the ranked [Formula: see text] operations as well as similarities and differences to known tree rearrangement based treespaces. Surprisingly, we show the counterintuitive result that adding leaves to trees can actually decrease their ranked [Formula: see text] distance, which may have an impact on the results of time tree sampling algorithms given uncertain "rogue taxa".

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.

Phylogeny of Alisma (Alismataceae) revisited: implications for polyploid evolution and species delimitation

Alisma L. is a genus of aquatic and wetland plants belonging to family Alismataceae. At present, it is thought to contain ten species. Variation in ploidy level is known in the genus, with diploids, tetraploids and hexaploids recorded. Previous molecular phylogenetic studies of Alisma have generated a robust backbone that reveals important aspects of the evolutionary history of this cosmopolitan genus, yet questions remain unresolved about the formation of the polyploid taxa and the taxonomy of one particularly challenging, widely distributed species complex. Here we directly sequenced, or cloned and sequenced, nuclear DNA (nrITS and phyA) and chloroplast DNA (matK, ndhF, psbA-trnH and rbcL) of multiple samples of six putative species and two varieties, and conducted molecular phylogenetic analyses. Alisma canaliculatum and its two varieties known in East Asia and A. rariflorum endemic to Japan possess closely related but heterogeneous genomes, strongly indicating that the two species were generated from two diploid progenitors, and are possibly siblings of one another. This evolutionary event may have occurred in Japan. Alisma canaliculatum var. canaliculatum is segregated into two types, each of which are geographically slightly differentiated in Japan. We reconstructed a single phylogeny based on the multi-locus data using Homologizer and then applied species delimitation analysis (STACEY). This allowed us to discern A. orientale as apparently endemic to the Southeast Asian Massif and distinct from the widespread A. plantago-aquatica. The former species was most likely formed through parapatric speciation at the southern edge of the distribution of the latter.

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.

Molecular Evidence Supports Five Lineages within Chiropotes (Pitheciidae, Platyrrhini)

Pitheciines have unique dental specializations among New World monkeys that allow them to feed on fruits with hard pericarps, thus playing a major role as seed predators. The three extant pitheciine genera, Pithecia, Cacajao and Chiropotes, are all endemic to the Amazon region. Because of the uncertainties about interspecific relationships, we reviewed the systematics and taxonomy of the genus Chiropotes. The phylogenetic analyses were performed based on Maximum Likelihood and Bayesian Inference, while species delimitation analyses were carried out using multispecies coalescent methods. In addition, we estimated genetic distances, divergence time and the probable ancestral distribution of this genus. Our results support five species of Chiropotes that emerged during the Plio-Pleistocene. Biogeographic estimates suggest that the ancestor of the current Chiropotes species occupied the endemism areas from Rondônia and Tapajós. Later, subsequent radiation and founder effects associated with the formation of the Amazonian basins probably determined the speciation events within Chiropotes.

Species delimitation, biogeography, and natural history of dwarf funnel web spiders (Mygalomorphae, Hexurellidae, Hexurella) from the United States / Mexico borderlands

The rarely encountered spider genus Hexurella Gertsch & Platnick, 1979 includes some of the smallest mygalomorph spiders in the world, with four poorly known taxa from central and southeastern montane Arizona, southern California, and northern Baja California Norte. At time of description the genus was known from fewer than 20 individuals, with sparse natural history information suggesting a vagrant, web-building, litter-dwelling natural history. Here the first published taxonomic and natural history information for this taxon is provided in more than 50 years, working from extensive new geographic sampling, consideration of male and female morphology, and sequence capture-based nuclear phylogenomics and mitogenomics. Several new species are easily diagnosed based on distinctive male morphologies, while a complex of populations from central and northern Arizona required an integrative combination of genomic algorithmic species delimitation analyses and morphological study. Four new species are described, including H.ephedrasp. nov., H.uwiiltilsp. nov., H.xericasp. nov., and H.zassp. nov. Females of H.encina Gertsch & Platnick, 1979 are also described for the first time. It is predicted that additional new species will ultimately be found in the mountains of central and northwestern Arizona, northern mainland Mexico, and the Mojave Desert of California.

Stop black and white thinking: Russula subgenus Compactae (Russulaceae, Russulales) in Europe revised

Russula subgenus Compactae is a group of ectomycorrhizal basidiomycetes, usually with large pileate fruitbodies. European members of the group are characterised by the absence of bright colours on the surfaces of their pilei, the context turning grey to black after cutting, the abundance of short lamellulae in the hymenophore, and spores with an inamyloid suprahilar spot and with low reticulate ornamentation. Our multi-locus phylogenetic study confirmed that this morphological delimitation corresponds to a well-supported clade. Within this clade, 16 species are recognised in Europe, of which five belong to the R. albonigra lineage and were described in a previous study, while eleven are fully described in this study. The application of the names R. acrifolia, R. adusta, R. anthracina, R. atramentosa, R. densissima, R. nigricans and R. roseonigra is based on the position of sequences retrieved from types or authentic material. Based on type sequences, R. fuliginosa is synonymised with R. anthracina and two varieties of R. anthracina are considered synonyms of R. atramentosa. The application of the name R. densifolia is based on a morphological match with the traditional species interpretation and the neotype specimen. Three species are described as new, R. marxmuelleriana sp. nov., R. picrophylla sp. nov. and R. thuringiaca sp. nov. This study recognises three major lineages and two species with isolated positions within the European Compactae and a morphological barcode was assigned to the species using an analysis of 23 selected characters. A search of publicly available sequences from the UNITE database revealed that the majority of species are host tree generalists and widely distributed in temperate and Mediterranean areas of Europe. Russula adusta is the only species so far proven to form ectomycorrhiza exclusively with conifers. Citation: De Lange R, Kleine J, Hampe F, et al. 2023. Stop black and white thinking: Russula subgenus Compactae (Russulaceae, Russulales) in Europe revised. Persoonia 51: 152-193. doi: 10.3767/persoonia.2023.51.04.

Diversification of the African legless skinks in the subfamily Acontinae (Family Scincidae)

Cladogenic diversification is often explained by referring to climatic oscillations and geomorphic shifts that cause allopatric speciation. In this regard, southern Africa retains a high level of landscape heterogeneity in vegetation, geology, and rainfall patterns. The legless skink subfamily Acontinae occurs broadly across the southern African subcontinent and therefore provides an ideal model group for investigating biogeographic patterns associated with the region. A robust phylogenetic study of the Acontinae with comprehensive coverage and adequate sampling of each taxon has been lacking up until now, resulting in unresolved questions regarding the subfamily's biogeography and evolution. In this study, we used multi-locus genetic markers (three mitochondrial and two nuclear) with comprehensive taxon coverage (all currently recognized Acontinae species) and adequate sampling (multiple specimens for most taxa) of each taxon to infer a phylogeny for the subfamily. The phylogeny retrieved four well-supported clades in Acontias and supported the monophyly of Typhlosaurus. Following the General Lineage Concept (GLC), many long-standing phylogenetic enigmas within Acontias occidentalis and the A. kgalagadi, A. lineatus and A. meleagris species complexes, and within Typhlosaurus were resolved. Our species delimitation analyses suggest the existence of hidden taxa in the A. occidentalis, A. cregoi and A. meleagris species groups, but also suggest that some currently recognized species in the A. lineatus and A. meleagris species groups, and within Typhlosaurus, should be synonymised. We also possibly encountered "ghost introgression" in A. occidentalis. Our inferred species tree revealed a signal of gene flow, which implies possible cross-over in some groups. Fossil evidence calibration dating results showed that the divergence between Typhlosaurus and Acontias was likely influenced by cooling and increasing aridity along the southwest coast in the mid-Oligocene caused by the opening of the Drake Passage. Further cladogenesis observed in Typhlosaurus and Acontias was likely influenced by Miocene cooling, expansion of open habitat, uplifting of the eastern Great Escarpment (GE), and variation in rainfall patterns, together with the effect of the warm Agulhas Current since the early Miocene, the development of the cold Benguela Current since the late Miocene, and their co-effects. The biogeographic pattern of the Acontinae bears close resemblance to that of other herpetofauna (e.g., rain frogs and African vipers) in southern Africa.

Geometric Morphometric Versus Genomic Patterns in a Large Polyploid Plant Species Complex

Plant species complexes represent a particularly interesting example of taxonomically complex groups (TCGs), linking hybridization, apomixis, and polyploidy with complex morphological patterns. In such TCGs, mosaic-like character combinations and conflicts of morphological data with molecular phylogenies present a major problem for species classification. Here, we used the large polyploid apomictic European Ranunculus auricomus complex to study relationships among five diploid sexual progenitor species and 75 polyploid apomictic derivate taxa, based on geometric morphometrics using 11,690 landmarked objects (basal and stem leaves, receptacles), genomic data (97,312 RAD-Seq loci, 48 phased target enrichment genes, 71 plastid regions) from 220 populations. We showed that (1) observed genomic clusters correspond to morphological groupings based on basal leaves and concatenated traits, and morphological groups were best resolved with RAD-Seq data; (2) described apomictic taxa usually overlap within trait morphospace except for those taxa at the space edges; (3) apomictic phenotypes are highly influenced by parental subgenome composition and to a lesser extent by climatic factors; and (4) allopolyploid apomictic taxa, compared to their sexual progenitor, resemble a mosaic of ecological and morphological intermediate to transgressive biotypes. The joint evaluation of phylogenomic, phenotypic, reproductive, and ecological data supports a revision of purely descriptive, subjective traditional morphological classifications.

Evaluating species boundaries using coalescent delimitation in pine-killing Monochamus (Coleoptera: Cerambycidae) sawyer beetles

Plant-feeding beetle species are diverse and often individually highly variable. Accurate classifications can be difficult to establish yet are essential for study of evolutionary patterns and processes. Molecular data are key to further characterizing morphologically difficult groups and defining genus and species boundaries. Monochamus Dejean species are ecologically and economically significant, and in coniferous forests they vector the nematode that causes Pine Wilt Disease. This study uses nuclear and mitochondrial genes to test the monophyly and relationships of Monochamus and applies coalescent methods to further delimit the conifer-feeding species. Monochamus has also included approximately 120 Old World species associated with diverse angiosperm tree species. We sample from these additional morphologically diverse species to determine their placement in the Lamiini. Through supermatrix and coalescent methods, the higher-level relationships of Monochamus show that conifer-feeders are a monophyletic group that includes the type species and has split into Nearctic and Palearctic clades. Molecular dating indicates a single dispersal of conifer-feeders to North America over the second Bering Land Bridge circa 5.3 Ma. All other Monochamus sampled fall in different parts of the Lamiini tree. Small-bodied angiosperm-feeding Monochamus group with the monotypic genus Microgoes Casey. The African Monochamus subgenera sampled are distantly related to the conifer-feeding clade. The multispecies coalescent delimitation methods BPP and STACEY delimit 17 conifer-feeding Monochamus species for a total of 18 species, and supports the retention of all current species. An interrogation with nuclear gene allele phasing reveals that unphased data can be unreliable for accurate delimitations and divergence times. The delimited species are discussed with integrative evidence, highlighting real-world challenges in recognizing the completion of speciation trajectories.

Phylogenomics of trans-Andean tetras of the genus Hyphessobrycon Durbin 1908 (Stethaprioninae: Characidae) and colonization patterns of Middle America

Hyphessobrycon is one of the most species rich and widely distributed genera in the family Characidae, with more than 160 species ranging from Veracruz, Mexico to Mar Chiquita Lagoon in Buenos Aires, Argentina. The majority of Hyphessobrycon diversity shows a cis-Andean distribution; only nine species are trans-Andean including H. compressus (Meek 1908). It is well established that Hyphessobrycon is not monophyletic but it has been suggested that natural groups can be identified within the larger Hyphessobrycon species group. In this study, we tested the monophyly of trans-Andean species of Hyphessobrycon and investigated the placement of H. compressus. We inferred the first phylogenomic hypothesis of trans-Andean Hyphessobrycon that includes nearly complete taxonomic sampling (eight of nine valid species) using ultraconserved elements (UCEs). We analyzed 75% (1682 UCEs), 90% (1258 UCEs), and 95% (838 UCEs) complete data matrices, and inferred phylogenomic hypotheses under concatenation and coalescent approaches. In all cases, we recovered the monophyly of trans-Andean Hyphessobrycon inclusive of H. compressus, strong support for three species groups, and evidence of cryptic diversity within the widespread H. compressus and H. condotensis. We used our phylogenomic hypothesis to investigate the biogeographic history of Hyphessobrycon in Middle America. Our ancestral range estimation analysis suggests a single event of cis- to trans-Andean colonization followed by stepwise colonization from the Pacific slope of northwestern South America (Chocó block) to northern Middle America (Maya block). Our work supports the recognition of the trans-Andean species as Hyphessobrycon sensu stricto and provides an evolutionary template to examine morphological characters that will allow us to better understand the diversity of Hyphessobrycon in Middle America.

Phylogenetic relationships of sleeper gobies (Eleotridae: Gobiiformes: Gobioidei), with comments on the position of the miniature genus Microphilypnus

Microphilypnus and Leptophilypnion are miniaturized genera within the family Eleotridae. The evolutionary relationships among these taxa are still poorly understood, and molecular analyses are restricted to mitochondrial genes, which have not been conclusive. We compiled both mitochondrial and nuclear genes to study the phylogenetic position of Microphilypnus and the evolutionary history and relationships of eleotrids. We propose that Microphilypnus and Leptophilypnus (a non-miniature genus) are not sister groups as suggested by previous studies, but rather separate lineages that arose in the early Eocene, with Leptophilypnus recovered as a sister group to the other analyzed eleotrids. In fact, Microphilypnus is currently associated with the Neotropical clade Guavina/Dormitator/Gobiomorus. We also identified a well-supported clade that indicated Gobiomorus and Hemieleotris as paraphyletic groups, besides a close relationship among Calumia godeffroyi, Bunaka gyrinoides, Eleotris and Erotelis species. This is the first comprehensive report about the evolutionary relationships in members of the family Eleotridae, including multiloci and multispecies approaches. Therefore, we provided new insights about the phylogenetic position of some taxa absent in previous studies, such as the miniature genus Microphilypnus and a recently described species of Eleotris from South America.

Taxonomy of Aspergillus series Versicolores: species reduction and lessons learned about intraspecific variability

Aspergillus series Versicolores members occur in a wide range of environments and substrates such as indoor environments, food, clinical materials, soil, caves, marine or hypersaline ecosystems. The taxonomy of the series has undergone numerous re-arrangements including a drastic reduction in the number of species and subsequent recovery to 17 species in the last decade. The identification to species level is however problematic or impossible in some isolates even using DNA sequencing or MALDI-TOF mass spectrometry indicating a problem in the definition of species boundaries. To revise the species limits, we assembled a large dataset of 518 strains. From these, a total of 213 strains were selected for the final analysis according to their calmodulin (CaM) genotype, substrate and geography. This set was used for phylogenetic analysis based on five loci (benA, CaM, RPB2, Mcm7, Tsr1). Apart from the classical phylogenetic methods, we used multispecies coalescence (MSC) model-based methods, including one multilocus method (STACEY) and five single-locus methods (GMYC, bGMYC, PTP, bPTP, ABGD). Almost all species delimitation methods suggested a broad species concept with only four species consistently supported. We also demonstrated that the currently applied concept of species is not sustainable as there are incongruences between single-gene phylogenies resulting in different species identifications when using different gene regions. Morphological and physiological data showed overall lack of good, taxonomically informative characters, which could be used for identification of such a large number of existing species. The characters expressed either low variability across species or significant intraspecific variability exceeding interspecific variability. Based on the above-mentioned results, we reduce series Versicolores to four species, namely A. versicolor, A. creber, A. sydowii and A. subversicolor, and the remaining species are synonymized with either A. versicolor or A. creber. The revised descriptions of the four accepted species are provided. They can all be identified by any of the five genes used in this study. Despite the large reduction in species number, identification based on phenotypic characters remains challenging, because the variation in phenotypic characters is high and overlapping among species, especially between A. versicolor and A. creber. Similar to the 17 narrowly defined species, the four broadly defined species do not have a specific ecology and are distributed worldwide. We expect that the application of comparable methodology with extensive sampling could lead to a similar reduction in the number of cryptic species in other extensively studied Aspergillus species complexes and other fungal genera. Citation: Sklenář F, Glässnerová K, Jurjević Ž, Houbraken J, Samson RA, Visagie CM, Yilmaz N, Gené J, Cano J, Chen AJ, Nováková A, Yaguchi T, Kolařík M, Hubka V (2022). Taxonomy of Aspergillus series Versicolores: species reduction and lessons learned about intraspecific variability. Studies in Mycology 102 : 53-93. doi: 10.3114/sim.2022.102.02.

A monograph of Aspergillus section Candidi

Aspergillus section Candidi encompasses white- or yellow-sporulating species mostly isolated from indoor and cave environments, food, feed, clinical material, soil and dung. Their identification is non-trivial due to largely uniform morphology. This study aims to re-evaluate the species boundaries in the section Candidi and present an overview of all existing species along with information on their ecology. For the analyses, we assembled a set of 113 strains with diverse origin. For the molecular analyses, we used DNA sequences of three house-keeping genes (benA, CaM and RPB2) and employed species delimitation methods based on a multispecies coalescent model. Classical phylogenetic methods and genealogical concordance phylogenetic species recognition (GCPSR) approaches were used for comparison. Phenotypic studies involved comparisons of macromorphology on four cultivation media, seven micromorphological characters and growth at temperatures ranging from 10 to 45 °C. Based on the integrative approach comprising four criteria (phylogenetic and phenotypic), all currently accepted species gained support, while two new species are proposed (A. magnus and A. tenebricus). In addition, we proposed the new name A. neotritici to replace an invalidly described A. tritici. The revised section Candidi now encompasses nine species, some of which manifest a high level of intraspecific genetic and/or phenotypic variability (e.g., A. subalbidus and A. campestris) while others are more uniform (e.g., A. candidus or A. pragensis). The growth rates on different media and at different temperatures, colony colours, production of soluble pigments, stipe dimensions and vesicle diameters contributed the most to the phenotypic species differentiation. Taxonomic novelties: New species: Aspergillus magnus Glässnerová & Hubka; Aspergillus neotritici Glässnerová & Hubka; Aspergillus tenebricus Houbraken, Glässnerová & Hubka. Citation: Glässnerová K, Sklenář F, Jurjević Ž, Houbraken J, Yaguchi T, Visagie CM, Gené J, Siqueira JPZ, Kubátová A, Kolařík M, Hubka V (2022). A monograph of Aspergillus section Candidi. Studies in Mycology 102: 1-51. doi: 10.3114/sim.2022.102.01.

Reducing the number of accepted species in Aspergillus series Nigri

The Aspergillus series Nigri contains biotechnologically and medically important species. They can produce hazardous mycotoxins, which is relevant due to the frequent occurrence of these species on foodstuffs and in the indoor environment. The taxonomy of the series has undergone numerous rearrangements, and currently, there are 14 species accepted in the series, most of which are considered cryptic. Species-level identifications are, however, problematic or impossible for many isolates even when using DNA sequencing or MALDI-TOF mass spectrometry, indicating a possible problem in the definition of species limits or the presence of undescribed species diversity. To re-examine the species boundaries, we collected DNA sequences from three phylogenetic markers (benA, CaM and RPB2) for 276 strains from series Nigri and generated 18 new whole-genome sequences. With the three-gene dataset, we employed phylogenetic methods based on the multispecies coalescence model, including four single-locus methods (GMYC, bGMYC, PTP and bPTP) and one multilocus method (STACEY). From a total of 15 methods and their various settings, 11 supported the recognition of only three species corresponding to the three main phylogenetic lineages: A. niger, A. tubingensis and A. brasiliensis. Similarly, recognition of these three species was supported by the GCPSR approach (Genealogical Concordance Phylogenetic Species Recognition) and analysis in DELINEATE software. We also showed that the phylogeny based on benA, CaM and RPB2 is suboptimal and displays significant differences from a phylogeny constructed using 5 752 single-copy orthologous proteins; therefore, the results of the delimitation methods may be subject to a higher than usual level of uncertainty. To overcome this, we randomly selected 200 genes from these genomes and performed ten independent STACEY analyses, each with 20 genes. All analyses supported the recognition of only one species in the A. niger and A. brasiliensis lineages, while one to four species were inconsistently delimited in the A. tubingensis lineage. After considering all of these results and their practical implications, we propose that the revised series Nigri includes six species: A. brasiliensis, A. eucalypticola, A. luchuensis (syn. A. piperis), A. niger (syn. A. vinaceus and A. welwitschiae), A. tubingensis (syn. A. chiangmaiensis, A. costaricensis, A. neoniger and A. pseudopiperis) and A. vadensis. We also showed that the intraspecific genetic variability in the redefined A. niger and A. tubingensis does not deviate from that commonly found in other aspergilli. We supplemented the study with a list of accepted species, synonyms and unresolved names, some of which may threaten the stability of the current taxonomy. Citation: Bian C, Kusuya Y, Sklenář F, D'hooge E, Yaguchi T, Ban S, Visagie CM, Houbraken J, Takahashi H, Hubka V (2022). Reducing the number of accepted species in Aspergillus series Nigri. Studies in Mycology 102: 95-132. doi: 10.3114/sim.2022.102.03.

Genomic comparisons confirm Giardia duodenalis sub-assemblage AII as a unique species

Giardia duodenalis is a parasitic flagellated protozoan which infects a wide range of mammalian hosts, including humans, and is subdivided into at least eight genetic assemblages commonly thought to represent cryptic species. Molecular studies have shown that G. duodenalis assemblage A, which parasitizes humans and animals, contains several phylogenetically distinct groupings known as sub-assemblages. Molecular studies employing poor phylogenetic-resolution markers routinely recover these sub-assemblages, implying that they represent evolutionarily distinct clades and possibly cryptic species, a hypothesis which is supported by epidemiologic trends. Here, we further tested this hypothesis by using available data from 41 whole genomes to characterize sub-assemblages and coalescent techniques for statistical estimation of species boundaries coupled to functional gene content analysis, thereby assessing the stability and distinctiveness of clades. Our analysis revealed two new sub-assemblage clades as well as novel signatures of gene content geared toward differential host adaptation and population structuring via vertical inheritance rather than recombination or panmixia. We formally propose sub-assemblage AII as a new species, Giardia hominis, while preserving the name Giardia duodenalis for sub-assemblage AI. Additionally, our bioinformatic methods broadly address the challenges of identifying cryptic microbial species to advance our understanding of emerging disease epidemiology, which should be broadly applicable to other lower eukaryotic taxa of interest. Giardia hominis n. sp. Zoobank LSID: urn:lsid: zoobank.org:pub:4298F3E1-E3EF-4977-B9DD-5CC59378C80E.

Integrative species delimitation reveals fine-scale allopatric speciation in a good-flying insect: a case study on

Alpha taxonomy is fundamental for many biological fields. Delineation of the species boundary, however, can be challenging in a species complex, where different species share a similar morphology and diagnostic characters may not be available. In this context, integrative approaches that incorporate molecular and morphological data sets, and account for speciation history can be helpful to alpha taxonomy. Different approaches to species delimitation based on different assumptions are complementary and by integrating the results from multiple approaches we can generate a more reliable and objective taxonomic decision. In this study, we applied three molecular approaches to species delimitation and inferred the demographic history based on an isolation with migration model to test a morphologically based taxonomic hypothesis for the Cylindera pseudocylindriformis complex. We discuss the association between genetic divergence and microhabitat specialisation, and further corroborate that C. subtilis sp. nov. is a valid new species by integrating the results from model-based species delimitation and the genealogical divergence index. We argue that genetic endemism can occur at a small geographic scale, even in a winged insect like tiger beetles. Our results also indicated that there may still be undocumented species diversity of Taiwanese Cylindera remaining to be discovered. ZooBank LSID: urn:lsid:zoobank.org:pub:9DEC1432-365C-4872-8D06-73B95F30624F

A first complete phylogenomic hypothesis for diploid blueberries (Vaccinium section Cyanococcus)

PREMISE

The true blueberries (Vaccinium sect. Cyanococcus; Ericaceae), endemic to North America, have been intensively studied for over a century. However, with species estimates ranging from nine to 24 and much confusion regarding species boundaries, this ecologically and economically valuable group remains inadequately understood at a basic evolutionary and taxonomic level. As a first step toward understanding the evolutionary history and taxonomy of this species complex, we present the first phylogenomic hypothesis of the known diploid blueberries.

METHODS

We used flow cytometry to verify the ploidy of putative diploid taxa and a target-enrichment approach to obtain a genomic data set for phylogenetic analyses.

RESULTS

Despite evidence of gene flow, we found that a primary phylogenetic signal is present. Monophyly for all morphospecies was recovered, with two notable exceptions: one sample of V. boreale was consistently nested in the V. myrtilloides clade and V. caesariense was nested in the V. fuscatum clade. One diploid taxon, Vaccinium pallidum, is implicated as having a homoploid hybrid origin.

CONCLUSIONS

This foundational study represents the first attempt to elucidate evolutionary relationships of the true blueberries of North America with a phylogenomic approach and sets the stage for multiple avenues of future study such as a taxonomic revision of the group, the verification of a homoploid hybrid taxon, and the study of polyploid lineages within the context of a diploid phylogeny.

Untying Gordian knots: unraveling reticulate polyploid plant evolution by genomic data using the large Ranunculus auricomus species complex

Speciation via hybridization and polyploidization is a major evolutionary force in plant evolution but is still poorly understood for neopolyploid groups. Challenges are attributed to high heterozygosity, low genetic divergence, and missing information on progenitors, ploidy, and reproduction. We study the large Eurasian Ranunculus auricomus species complex and use a comprehensive workflow integrating reduced-representation sequencing (RRS) genomic data to unravel reticulate evolution, genome diversity and composition of polyploids. We rely on 97 312 restriction site-associated DNA sequencing (RAD-Seq) loci, 576 targeted nuclear genes (48 phased), and 71 plastid regions derived from 78 polyploid apomictic taxa and four diploid and one tetraploid putative sexual progenitor species. We applied (phylo)genomic structure, network, and single nucleotide polymorphism (SNP)-origin analyses. Results consistently showed only 3-5 supported and geographically structured polyploid genetic groups, each containing extant sexual and one unknown progenitor species. Combined analyses demonstrated predominantly allopolyploid origins, each involving 2-3 different diploid sexual progenitor species. Young allotetraploids were characterized by subgenome dominance and nonhybrid SNPs, suggesting substantial post-origin but little lineage-specific evolution. The biodiversity of neopolyploid complexes can result from multiple hybrid origins involving different progenitors and substantial post-origin evolution (e.g. homoeologous exchanges, hybrid segregation, gene flow). Reduced-representation sequencing genomic data including multi-approach information is efficient to delimit shallow reticulate relationships.

Consolidation of Chloridium: new classification into eight sections with 37 species and reinstatement of the genera Gongromeriza and Psilobotrys

Chloridium is a little-studied group of soil- and wood-inhabiting dematiaceous hyphomycetes that share a rare mode of phialidic conidiogenesis on multiple loci. The genus has historically been divided into three morphological sections, i.e. Chloridium, Gongromeriza, and Psilobotrys. Sexual morphs have been placed in the widely perceived genus Chaetosphaeria, but unlike their asexual counterparts, they show little or no morphological variation. Recent molecular studies have expanded the generic concept to include species defined by a new set of morphological characters, such as the collar-like hyphae, setae, discrete phialides, and penicillately branched conidiophores. The study is based on the consilience of molecular species delimitation methods, phylogenetic analyses, ancestral state reconstruction, morphological hypotheses, and global biogeographic analyses. The multilocus phylogeny demonstrated that the classic concept of Chloridium is polyphyletic, and the original sections are not congeneric. Therefore, we abolish the existing classification and propose to restore the generic status of Gongromeriza and Psilobotrys. We present a new generic concept and define Chloridium as a monophyletic, polythetic genus comprising 37 species distributed in eight sections. In addition, of the taxa earlier referred to Gongromeriza, two have been redisposed to the new genus Gongromerizella. Analysis of published metabarcoding data showed that Chloridium is a common soil fungus representing a significant (0.3 %) proportion of sequence reads in environmental samples deposited in the GlobalFungi database. The analysis also showed that they are typically associated with forest habitats, and their distribution is strongly influenced by climate, which is confirmed by our data on their ability to grow at different temperatures. We demonstrated that Chloridium forms species-specific ranges of distribution, which is rarely documented for microscopic soil fungi. Our study shows the feasibility of using the GlobalFungi database to study the biogeography and ecology of fungi. Taxonomic novelties: New genus: Gongromerizella Réblová; New sections: Chloridium section Cryptogonytrichum Réblová, Hern.-Restr., M. Kolařík & F. Sklenar, Chloridium section Gonytrichopsis Réblová, Hern.-Restr., M. Kolařík & F. Sklenar, Chloridium section Metachloridium Réblová, Hern.-Restr., M. Kolařík & F. Sklenar, Chloridium section Volubilia Réblová, Hern.-Restr., M. Kolařík & F. Sklenar; New species: Chloridium bellum Réblová & Hern.-Restr., Chloridium biforme Réblová & Hern.-Restr., Chloridium detriticola Réblová & Hern.-Restr., Chloridium gamsii Réblová & Hern.-Restr., Chloridium guttiferum Réblová & Hern.-Restr., Chloridium moratum Réblová & Hern.-Restr., Chloridium peruense Réblová & Hern.-Restr., Chloridium novae-zelandiae Réblová & Hern.-Restr., Chloridium elongatum Réblová & Hern.-Restr., Chloridium volubile Réblová & Hern.-Restr.; New varieties: Chloridium bellum var. luteum Réblová & Hern.-Restr., Chloridium detriticola var. effusum Réblová & Hern.-Restr., Chloridium chloridioides var. convolutum Réblová & Hern.-Restr.; New combinations: Chloridium section Gonytrichum (Nees & T. Nees) Réblová, Hern.-Restr., M. Kolařík & F. Sklenar, Chloridium section Mesobotrys (Sacc.) Réblová, Hern.-Restr., M. Kolařík & F. Sklenar, Chloridium section Pseudophialocephala (M.S. Calabon et al.) Réblová, Hern.-Restr., M. Kolařík & F. Sklenar, Chloridium simile (W. Gams & Hol.-Jech.) Réblová & Hern.-Restr., Chloridium chloridioides (W. Gams & Hol.-Jech.) Réblová & Hern.-Restr., Chloridium subglobosum (W. Gams & Hol.-Jech.) Réblová & Hern.-Restr., Chloridium fuscum (Corda) Réblová & Hern.-Restr., Chloridium ypsilosporum (Hol.-Jech.) Réblová & Hern.-Restr., Chloridium costaricense (G. Weber et al.) Réblová & Hern.-Restr., Chloridium cuneatum (N.G. Liu et al.) Réblová & Hern.-Restr., Fusichloridium cylindrosporum (W. Gams & Hol.-Jech.) Réblová, Gongromeriza myriocarpa (Fr.) Réblová, Gongromeriza pygmaea (P. Karst.) Réblová, Gongromerizella lignicola (F. Mangenot) Réblová, Gongromerizella pachytrachela (W. Gams & Hol.-Jech) Réblová, Gongromerizella pini (Crous & Akulov) Réblová; New name: Chloridium pellucidum Réblová & Hern.-Restr.; Epitypifications (basionyms): Chaetopsis fusca Corda, Gonytrichum caesium var. subglobosum W. Gams & Hol.-Jech.; Lectotypification (basionym): Gonytrichum caesium Nees & T. Nees. Citation: Réblová M, Hernández-Restrepo M, Sklenář F, Nekvindová J, Réblová K, Kolařík M (2022). Consolidation of Chloridium: new classification into eight sections with 37 species and reinstatement of the genera Gongromeriza and Psilobotrys. Studies in Mycology 103: 87-212. doi: 10.3114/sim.2022.103.04.

The evolutionary history of the relict scorpion family Iuridae of the eastern Mediterranean

Iuridae is a family of scorpions that exhibits a highly complex biogeographic and taxonomic history. Iuridae taxa are mainly found in Turkey and Greece, whereas a single species is found in northern Iraq. Several taxonomic revisions have been conducted on this family that initially comprised two genera. The latest taxonomic review, based on morphological and anatomical features, raised the number of Iuridae genera to four, and the number of species to 14. Sequence data from three molecular markers (COX1, 16S rDNA, ITS1) originating from numerous Iuridae taxa were analyzed within a phylogenetic framework. Divergence time-estimate analyses, species delimitation approaches and estimation of ancestral areas were implemented in order to: 1) reconstruct the phylogenetic relationships of the Iuridae taxa, 2) evaluate the morphological classifications, and 3) obtain insights into the biogeographic history of the family in the East Mediterranean. The multi-locus phylogeny clearly confirms an ancient division into two clades, Calchinae and Iurinae. Ancient patterns of isolation and dispersal are revealed. Both subfamilies are largely confined to the Anatolian peninsula and its few coastal islands; only the most derived genus Iurus has dispersed westward to Crete and Peloponnese. Based on our findings, three new genera of Iurinae (Metaiurus, Anatoliurus, and Letoiurus) are established. The genus Neocalchas emerges as one of the most ancient scorpion clades, with divergence time about 27 mya.

A tale of two bellies: systematics of the oval frogs (Anura: Microhylidae: Elachistocleis Parker, 1927)

Oval frogs (Elachistocleis) have a broad geographic distribution covering nearly all of South America and parts of Central America. They also have a large inter- and intraspecific variation of the few morphological characters commonly used as diagnostic traits among species of the genus. Based on molecular data, we provide the most complete phylogeny of Elachistocleis to date, and explore its genetic diversity using distance-based and tree-based methods for putative species delimitation. Our results show that at least two of the most relevant traditional characters used in the taxonomy of this group (belly pattern and dorsal median white line) carry less phylogenetic information than previously thought. Based on our results, we propose some synonymizations and some candidate new species. This study is a first major step in disentangling the current systematics of Elachistocleis. Yet, a comprehensive review of morphological data is needed before any new species descriptions can be properly made.

Assessing diversity of King Crab Lithodes spp. in the south‐eastern pacific using phylogeny and molecular species delimitation methods

The purpose of this study was to test the hypothesis that the genetic diversity of commercially significant species of King Crabs (Lithodes spp.) along the south‐eastern Pacific (SEP) comprises different independent evolutionary units (IEUs) with spatially isolated distribution. Nine localities from inner and open waters along the SEP Chilean coast (39°S‐55°S) were sampled. We analyzed sequences from 173 individuals for the mitochondrial gene Cytochrome oxidase I (COX‐I), 151 individuals for the Internal Transcribed Spacer 1 (ITS) and 135 for the structural ribosomal RNA (28S). Genetic delimitation was performed through three analytical methods: ABGD, GMYC, and its Bayesian implementation, bGMYC. Bayesian phylogenetic analyses and haplotype networks were also performed. Divergence time between clades was assessed for the COX‐I marker and estimated from known evolutionary rates for this marker in other crustacean species and fossil calibration from other Anomuran species. Delimitation analyses, phylogenetic analyses, and mitochondrial haplotype networks suggested the presence of two deeply divergent mitochondrial lineages of Lithodes in the SEP, referred to as Clade1 and Clade 2. Nuclear markers showed low phylogenetic resolution and therefore were unsuitable for molecular species delimitation. Divergence time analysis of the mitochondrial lineages suggests a separation between Clades of approximately 2.3 Mya. The divergence time obtained suggested that Pliocene glaciations and deglaciations cycles could be involved in hybridization events between Lithodes IEUs at southern tip of South American coasts. The different frequencies of Lithodes haplotypes in inner and open water environments along SEP coasts could be explained by events such as the last glacial maximum or by differences in the adaptation of each clade to different environments. These findings support the necessity of evaluating the taxonomic status of Lithodes individuals found along SEP coasts under an integrative taxonomy approach or through markers with other evolution rates than those already used.

Bayesian phylogenetic inference using relaxed-clocks and the multispecies coalescent

The multispecies coalescent (MSC) model accommodates both species divergences and within-species coalescent and provides a natural framework for phylogenetic analysis of genomic data when the gene trees vary across the genome. The MSC model implemented in the program bpp assumes a molecular clock and the Jukes–Cantor model, and is suitable for analyzing genomic data from closely related species. Here we extend our implementation to more general substitution models and relaxed clocks to allow the rate to vary among species. The MSC-with-relaxed-clock model allows the estimation of species divergence times and ancestral population sizes using genomic sequences sampled from contemporary species when the strict clock assumption is violated, and provides a simulation framework for evaluating species tree estimation methods. We conducted simulations and analyzed two real datasets to evaluate the utility of the new models. We confirm that the clock-JC model is adequate for inference of shallow trees with closely related species, but it is important to account for clock violation for distant species. Our simulation suggests that there is valuable phylogenetic information in the gene-tree branch lengths even if the molecular clock assumption is seriously violated, and the relaxed-clock models implemented in bpp are able to extract such information. Our Markov chain Monte Carlo algorithms suffer from mixing problems when used for species tree estimation under the relaxed clock and we discuss possible improvements. We conclude that the new models are currently most effective for estimating population parameters such as species divergence times when the species tree is fixed.

Quantitatively defining species boundaries with more efficiency and more biological realism

We introduce a widely applicable species delimitation method based on the multispecies coalescent model that is more efficient and more biologically realistic than existing methods. We extend a threshold-based method to allow the ancestral speciation rate to vary through time as a smooth piecewise function. Furthermore, we introduce the cutting-edge proposal kernels of StarBeast3 to this model, thus enabling rapid species delimitation on large molecular datasets and allowing the use of relaxed molecular clock models. We validate these methods with genomic sequence data and SNP data, and show they are more efficient than existing methods at achieving parameter convergence during Bayesian MCMC. Lastly, we apply these methods to two datasets (Hemidactylus and Galagidae) and find inconsistencies with the published literature. Our methods are powerful for rapid quantitative testing of species boundaries in large multilocus datasets and are implemented as an open source BEAST 2 package called SPEEDEMON.

Introducing SPEEDEMON, a package for BEAST 2 that better defines species boundaries based on molecular data demonstrated on gecko and loris datasets.

Phylogeny of Telmatobius marmoratus complex (Anura, Telmatobiidae) reveals high cryptic diversity in the Andean Altiplano

Telmatobius is the most diverse group of anurans in the Andean Altiplano (highlands) Morphologically, these amphibians have a generally conserved morphology but in turn present large intraspecific variation, which has led to a complex taxonomy and systematics. T. marmoratus has the widest distribution of the genus and forms a complex composed of at least two Telmatobius species. Partial systematic studies based on molecular evidence reveal the existence of three lineages with a complex spatial distribution. However, these studies did not include the entire distribution of T. marmoratus. Our study aims to reassess the current systematic scenario including the complete distribution of the complex. For this, we used a multilocus approach based on mitochondrial (16S, Cytb) and nuclear (RAG1-1, BFIB) DNA sequences to build a phylogenetic hypothesis based on Bayesian inference, maximum likelihood and maximum parsimony. Subsequently, we performed single-locus (ABGD and PTP) and multilocus (STACEY) species delimitation analyses to verify the diversity of nominal species within the complex. The analyses suggest seven non-sibling lineages and 6-10 candidate species within the marmoratus complex. Only one of the two lineages restricted to the central northern plateau correspond to T. marmoratus sensu stricto. South-central marbled water frogs belong to completely new lineages closer to T. gigas and T. culeus, evidencing the polyphyletic condition of the marmoratus complex. The findings of several sympatric lineages in some localities reveal a complex history of ancient water connections in south-central Altiplano.

Phylogenetic assessment of the Metamasius hemipterus species complex (Coleoptera, Curculionidae, Dryophthorinae)

Metamasius is a large genus of dryophthorine weevils, with nearly 85 species. Among the economically important pests in the genus, M. hemipterus is currently separated in three subspecies, based largely on color patterns of the elytra, pronotum, and sternum. The tenuous limits of M. hemipterus subspecies were created over fifty years ago and never tested under a phylogenetic framework. Here, for the first time, we address the M. hemipterus species boundaries applying a molecular approach. We constructed a reduced genome representation of a few species using restriction site-associated DNA sequencing (RADseq). Phylogenetic analysis using either a complete supermatrix or only SNPs revealed a clear separation of Metamasius species. We suggest that M. h. carbonarius syn. nov. and M. h. sericeus be treated as the same species, M. sericeus (Oliver) stat. n., and elevate M. h. hemipterus as a separate species M. hemipterus (Linnaeus). We updated Vaurie's identification key to reflect the new species status. This systematic reassessment reflects a more natural classification for these remarkable and economically significant weevils.

Convergent Morphological Evolution in Silene Sect. Italicae (Caryophyllaceae) in the Mediterranean Basin

Recent divergence can obscure species boundaries among closely related taxa. Silene section Italicae (Caryophyllaceae) has been taxonomically controversial, with about 30 species described. We investigate species delimitation within this section using 500 specimens sequenced for one nuclear and two plastid markers. Despite the use of a small number of genes, the large number of sequenced samples allowed confident delimitation of 50% of the species. The delimitation of other species (e.g., Silene nemoralis, S. nodulosa and S. andryalifolia) was more challenging. We confirmed that seven of the ten chasmophyte species in the section are not related to each other but are, instead, genetically closer to geographically nearby species belonging to Italicae yet growing in open habitats. Adaptation to chasmophytic habitats therefore appears to have occurred independently, as a result of convergent evolution within the group. Species from the Western Mediterranean Basin showed more conflicting species boundaries than species from the Eastern Mediterranean Basin, where there are fewer but better-delimited species. Significant positive correlations were found between an estimation of the effective population size of the taxa and their extent of occurrence (EOO) or area of occupancy (AOO), and negative but non-significant correlations between the former and the posterior probability (PP) of the corresponding clades. These correlations might suggest a lower impact of incomplete lineage sorting in species with low effective population sizes and small distributional ranges compared with that in species inhabiting large areas. Finally, we confirmed that S. italica and S. nemoralis are distinct species, that S. nemoralis might furthermore include two different species and that S. velutina from Corsica and S. hicesiae from the Lipari Islands are sister species.

DNA barcoding, multilocus phylogeny, and morphometry reveal phenotypic plasticity in the Chinese freshwater mussel Lamprotula caveata (Bivalvia: Unionidae)

Accurate species identification is crucial for developing conservation strategies for freshwater mussels, one of the most imperiled faunas in the world. Traditionally, mussel species description primarily relied on conchological characters. However, shell morphology has great variability, which leads to the complexity of species delimitation. As endemic species to China, Lamprotula caveata was originally described by Heude (1877). Lamprotula quadrangulosus and Lamprotula contritus were considered for synonymization of L. caveata based on shell variants in the early 20th century, which has been long debated due to lack of rigorous molecular analysis. Moreover, great morphological variation caused doubt whether there are cryptic species. In this study, we used a combined phylogenetic and morphometric approach to verify the validity of the synonymization of L. caveata. The results of molecular species delimitation showed that two molecular operational taxonomic units (MOTUs) were identified in Lamprotula spp., including the L. leaii lineage and the complex lineage (L. quadrangulosa, L. cornuumlunae, L. contritus, and L. caveata). Phylogenetic analyses revealed that L. cornuumlunae formed a basal monophyletic clade, whose divergence time was relatively recent (4.26 Ma [95% HPD = 1.91–7.22 Ma]), and L. contritus, L. caveata, and L. quadrangulosa formed a large polytomy group with very shallow branches. In the previous study, we have demonstrated the validity of L. cornuumlunae. The molecular evidences supported that the complex (L. quadrangulosa + L. contritus + L. caveata) was a valid species; L. quadrangulosa and L. contritus were synonyms of L. caveata. In addition, three morphospecies (L. quadrangulosa, L. contritus, and L. caveata) were aggregated without clear differentiation based on shell morphometric analysis. We confirmed multiple phenotypes in L. caveata for species identification and presumed that the phenotypic plasticity was a response to specific habitats. This study clarified the diversity and phylogeny of the Lamprotula group, which is a crucial step for developing new conservation and management strategies for this imperiled group.

Multilocus Data Analysis Reveal the Diversity of Cryptic Species in the Tillandsia ionantha (Bromeliaceae: Tillansiodeae) Complex

Independent evolutionary lineages or species that lack phenotypic variation as an operative criterion for their delimitation are known as cryptic species. However, these have been delimited using other data sources and analysis. The aims of this study are: (1) to evaluate the divergence of the populations of the T. ionantha complex; and (2) to delimit the species using multilocus data, phylogenetic analysis and the coalescent model. Phylogenetic analyses, genetic diversity and population structure, and isolation by distance analysis were performed. A multispecies coalescent analysis to delimit the species was conducted. Phylogenetic analysis showed that T. ionantha is polyphyletic composed of eight evolutionary lineages. Haplotype distribution and genetic differentiation analysis detected strong population structure and high values of genetic differentiation among populations. The positive correlation between genetic differences with geographic distance indicate that the populations are evolving under the model of isolation by distance. The coalescent multispecies analysis performed with starBEAST supports the recognition of eight lineages as different species. Only three out of the eight species have morphological characters good enough to recognize them as different species, while five of them are cryptic species. Tillandsia scaposa and T. vanhyningii are corroborated as independent lineages, and T. ionantha var. stricta changed status to the species level.

Reticulate Evolution in the Western Mediterranean Mountain Ranges: The Case of the Leucanthemopsis Polyploid Complex

Polyploidization is one of the most common speciation mechanisms in plants. This is particularly relevant in high mountain environments and/or in areas heavily affected by climatic oscillations. Although the role of polyploidy and the temporal and geographical frameworks of polyploidization have been intensively investigated in the alpine regions of the temperate and arctic biomes, fewer studies are available with a specific focus on the Mediterranean region. Leucanthemopsis (Asteraceae) consists of six to ten species with several infraspecific entities, mainly distributed in the western Mediterranean Basin. It is a polyploid complex including montane, subalpine, and strictly alpine lineages, which are locally distributed in different mountain ranges of Western Europe and North Africa. We used a mixed approach including Sanger sequencing and (Roche-454) high throughput sequencing of amplicons to gather information from single-copy nuclear markers and plastid regions. Nuclear regions were carefully tested for recombinants/PCR artifacts and for paralogy. Coalescent-based methods were used to infer the number of polyploidization events and the age of formation of polyploid lineages, and to reconstruct the reticulate evolution of the genus. Whereas the polyploids within the widespread Leucanthemopsis alpina are autopolyploids, the situation is more complex among the taxa endemic to the western Mediterranean. While the hexaploid, L. longipectinata, confined to the northern Moroccan mountain ranges (north-west Africa), is an autopolyploid, the Iberian polyploids are clearly of allopolyploid origins. At least two different polyploidization events gave rise to L. spathulifolia and to all other tetraploid Iberian taxa, respectively. The formation of the Iberian allopolyploids took place in the early Pleistocene and was probably caused by latitudinal and elevational range shifts that brought into contact previously isolated Leucanthemopsis lineages. Our study thus highlights the importance of the Pleistocene climatic oscillations and connected polyploidization events for the high plant diversity in the Mediterranean Basin.

Unveiling the Mycodrosophila projectans (Diptera, Drosophilidae) species complex: Insights into the evolution of three Neotropical cryptic and syntopic species

The Zygothrica genus group has been shown to be speciose, with a high number of cryptic species. DNA barcoding approaches have been a valuable tool to uncover cryptic diversity in this lineage, as recently suggested for the Neotropical Mycodrosophila projectans complex, which seems to comprise at least three different species. The aim of this study was to confirm the subdivision of the M. projectans complex while shedding some light on the patterns and processes related to its diversification. In this sense, the use of single and multi-locus datasets under phylogenetic, distance, coalescence, and diagnostic nucleotide approaches confirmed the presence of at least three species under the general morphotype previously described as M. projectans. Only a few subtle morphological differences were found for the three species in terms of aedeagus morphology and abdominal color patterns. Ecologically, sympatry and syntopy seem to be recurrent for these three cryptic species, which present widely overlapping niches, implying niche conservatism. This morphological and ecological similarity has persisted though cladogenesis within the complex, which dates back to the Miocene, providing an interesting example of morphological conservation despite ancient divergence. These results, in addition to contrasting patterns of past demographic fluctuations, allowed us to hypothesize patterns of allopatric or parapatric diversification with secondary contact in Southern Brazil. Nevertheless, genetic diversity was generally high within species, suggesting that migration may encompass an adaptive response to the restrictions imposed by the ephemerality of resources.

Phylogeography and evolutionary history of the Panamic Clingfish Gobiesox adustus in the Tropical Eastern Pacific

The Panamic Clingfish Gobiesox adustus is widely distributed in the Tropical Eastern Pacific (TEP), from the central Gulf of California, Mexico to Ecuador, including the oceanic Revillagigedo Archipelago, and Isla del Coco. This cryptobenthic species is restricted to very shallow rocky-reef habitats. Here, we used one mitochondrial and three nuclear DNA markers from 155 individuals collected across the distribution range of the species in order to evaluate if geographically structured populations exist and to elucidate its evolutionary history. Phylogenetic analyses recovered a monophyletic group, with four well-supported, allopatric subgroups. Each subgroup corresponded to one of the following well-known biogeographic regions/provinces: 1) the Revillagigedo Archipelago, 2) the Cortez + Mexican provinces (Mexico), 3) the Panamic province (from El Salvador to Ecuador), and 4) Isla del Coco. A molecular-clock analysis showed a mean date for the divergence between clade I (the Revillagigedos and Cortez + Mexican provinces) and clade II (Panamic province and Isla del Coco) in the Pliocene, at ca. 5.33 Mya. Within clade I, the segregation between the Revillagigedos and Cortez + Mexican province populations was dated at ca. 1.18 Mya, during the Pleistocene. Within clade II, the segregation between samples of Isla del Coco and the Panamic province samples was dated at ca. 0.77 Mya, during the Pleistocene. The species tree, Bayesian species delimitation tests (BPP and STACEY), the Φ_ST, AMOVA, and the substantial genetic distances that exist between those four subgroups, indicate that they are independent evolutionary units. These cladogenetic events seem to be related to habitat discontinuities, and oceanographic and geological processes that produce barriers to gene flow for G. adustus, effects of which are enhanced by the intrinsic ecological characteristics of this species.

Molecular data confirm Triatoma pallidipennis Stål, 1872 (Hemiptera: Reduviidae: Triatominae) as a novel cryptic species complex

Triatoma pallidipennis constitues one of the most important Chagas disease vector in Mexico. Previous studies based on molecular data suggest T. pallidipennis as a complex of cryptic species. For that reason, we analyzed the phylogenetic relationships of T. pallidipennis using DNA sequences from the mitochondrial ND4 gene and the ITS-2 gene. In addition, the divergence times were estimated, and possible new taxa were delimited with three species delimitation methods. Finally, genetic distances and possible connectivity routes based on shared haplotypes were obtained among the T. pallidipennis populations. Five haplogroups (possible cryptic species) were found, based on delimitation methods and genetic distances. Haplogroup divergence began about 3 Ma, in the Pleistocene. Moreover, none of the haplogroups showed potential connectivity routes between them, evidencing lack of gene flow. Our results suggest the existence of a new cryptic species complex within what is currently recognized as a T. pallidipennis.

Molecular phylogeny of the genus Chondrina (Gastropoda, Panpulmonata, Chondrinidae) in the Iberian Peninsula

Chondrina Reichenbach, 1828 is a highly diverse genus of terrestrial molluscs currently including 44 species with about 28 subspecific taxa. It is distributed through North Africa, central and southern Europe, from Portugal in the West to the Caucasus and Asia Minor in the East. Approximately 70% of the species are endemic to the Iberian Peninsula constituting its main center of speciation with 34 species. This genus includes many microendemic taxa, some of them not yet described, confined to limestone habitats (being strictly rock-dwelling species). They are distributed on rocky outcrops up to 2000 m.a.s.l. It is a genus of conical-fusiform snails that differ mainly in shell characters and in the number and position of teeth in their aperture. So far, molecular studies on Chondrina have been based exclusively on the mitochondrial Cytochrome Oxidase subunit I region (COI). These studies gave a first view of the phylogeny of the genus but many inner nodes were not statistically supported. The main objective of the study is to obtain a better understanding of the phylogeny and systematics of the genus Chondrina on the Iberian Peninsula, using multilocus molecular analysis. Partial sequences of the COI and 16S rRNA genes, as well as of the nuclear Internal Transcribed Spacer 1 (ITS1-5.8S) and Internal Transcribed Spacer 2 (5.8S-ITS2-28S) were obtained from individuals of all the extant Chondrina species known from the Iberian Peninsula. In addition to this, the newly obtained COI sequences were combined with those previously published in the GenBank. Phylogenetic relationships were inferred using maximum likelihood and Bayesian methods. The reconstructed phylogenies showed high values of support for more recent branches and basal nodes. Moreover, molecular species delimitation allowed to better definethe studied species and check the presence of new taxa.

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

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

An Efficient Coalescent Epoch Model for Bayesian Phylogenetic Inference

We present a two-headed approach called Bayesian Integrated Coalescent Epoch PlotS (BICEPS) for efficient inference of coalescent epoch models. Firstly, we integrate out population size parameters, and secondly, we introduce a set of more powerful Markov chain Monte Carlo (MCMC) proposals for flexing and stretching trees. Even though population sizes are integrated out and not explicitly sampled through MCMC, we are still able to generate samples from the population size posteriors. This allows demographic reconstruction through time and estimating the timing and magnitude of population bottlenecks and full population histories. Altogether, BICEPS can be considered a more muscular version of the popular Bayesian skyline model. We demonstrate its power and correctness by a well-calibrated simulation study. Furthermore, we demonstrate with an application to SARS-CoV-2 genomic data that some analyses that have trouble converging with the traditional Bayesian skyline prior and standard MCMC proposals can do well with the BICEPS approach. BICEPS is available as open-source package for BEAST 2 under GPL license and has a user-friendly graphical user interface.[Bayesian phylogenetics; BEAST 2; BICEPS; coalescent model.].

Sea–land transition drove terrestrial amphipod diversification in East Asia, with a description of a new species

Sea–land transition caused by marine incursion and regression is hypothesized to be a major driving force in terrestrial biodiversity, providing opportunities for marine ancestors to colonize terrestrial habitats and driving vicariant speciation in distinct geographical regions. Here, we test this hypothesis in East Asia using amphipods of the Morinoia japonica complex. We constructed a dataset from 269 individuals covering all known ranges of this species complex. Phylogenetic and biogeographic analyses revealed that sea–land transition during the Miocene drove the coastal ancestor to invade terrestrial habitats in East Asia and subsequently split into eight biogeographic lineages in eastern China, Korea, Japan and some Pacific islands. Stepping-stone dispersal resulted in a relatively wide distribution of M. japonica, and long-term geographical isolation led to the diversification of the M. japonica complex. Species delimitation analysis suggests that this complex contains eight species. We describe the geographical group from eastern China as a new species, Morinoia aosen sp. nov., based on genetic and morphological comparisons with other geographical groups. Type specimens are deposited in the Institute of Zoology, Chinese Academy of Sciences (IZCAS) in Beijing, China.

Molecular systematics of the Reithrodontomys tenuirostris group (Rodentia: Cricetidae) highlighting the Reithrodontomys microdon species complex

The Reithrodontomys tenuirostris species group is considered "the most specialized" within the genus Reithrodontomys from morphological and ecological perspectives. Previous studies based on molecular data recommended changes in the taxonomy of the group. In particular, R. microdon has been the most taxonomically questioned, with the suggestion that it constitutes a complex of cryptic species. We analyzed the phylogenetic relationships of the R. tenuirostris species group using DNA sequences from the mitochondrial Cytochrome b gene and Intron 7 of the nuclear beta fibrinogen gene. In addition, divergence times were estimated, and possible new taxa delimited with three widely used species delimitation methods. Finally, possible connectivity routes based on shared haplotypes were tested among the R. microdon populations. All species were recovered as monophyletic with the exception of R. microdon, whose individuals were grouped into four different haplogroups, one of which included specimens of R. bakeri. Diversification within the R. tenuirostris species group began about 3 Ma, in the Pleistocene. The bGMYC and STACEY delimitation methods were congruent with each other, delimiting at the species-level each haplogroup within R. microdon, while the mPTP suggested a greater number of species. Moreover, none of the haplogroups showed potential connectivity routes between them, evidencing lack of gene flow. Our results suggest the existence of a higher number of species in the R. tenuirostris group, because we show that there are four species within what is currently recognized as R. microdon.

Polyozellus vs. Pseudotomentella: generic delimitation with a multi-gene dataset

Polyozellus and Pseudotomentella are two genera of closely related, ectomycorrhizal fungi in the order Thelephorales; the former stipitate and the latter corticioid. Both are widespread in the Northern Hemisphere and many species from both genera seem to be restricted to old growth forest. This study aimed to: a) identify genetic regions useful in inferring the phylogenetic relationship between Polyozellus and Pseudotomentella, b) infer this relationship with the regions identified and c) make any taxonomic changes warranted by the result. RPB2, mtSSU and nearly full-length portions of nrLSU and nrSSU were found to be comparatively easy to sequence and provide a strong phylogenetic signal. A STACEY species tree of these three regions revealed that Polyozellus makes Pseudotomentella paraphyletic. As a result, nearly all species currently placed in Pseudotomentella were recombined to Polyozellus. Pseudotomentella larsenii was found to be closer to Tomentellopsis than Polyozellus, but its placement needs further study and it was hence not recombined.

Multispecies coalescent and its applications to infer species phylogenies and cross-species gene flow

Multispecies coalescent (MSC) is the extension of the single-population coalescent model to multiple species. It integrates the phylogenetic process of species divergences and the population genetic process of coalescent, and provides a powerful framework for a number of inference problems using genomic sequence data from multiple species, including estimation of species divergence times and population sizes, estimation of species trees accommodating discordant gene trees, inference of cross-species gene flow and species delimitation. In this review, we introduce the major features of the MSC model, discuss full-likelihood and heuristic methods of species tree estimation and summarize recent methodological advances in inference of cross-species gene flow. We discuss the statistical and computational challenges in the field and research directions where breakthroughs may be likely in the next few years.

Dendrodacrys: a new genus for species with branched hyphidia in Dacrymyces s.l., with the description of four new species

A new genus named Dendrodacrys is proposed for a monophyletic group in Dacrymycetaceae, containing species with pulvinate to depressed basidiocarps, distinctly branched hymenial hyphidia, and up to 3-septate mature basidiospores. Four taxa in this group, occurring in Europe, are proposed as new species, viz. De. ciprense, De. concrescens, De. ellipsosporum, and De. oblongisporum, based both on morphological and DNA data (nrDNA, RPB1, RPB2, TEF-1α, 12S). These new species are all described in detail, illustrated, and compared with other published taxa that with which they can be confounded. The new combination De. paraphysatum is proposed after revising the type material of Dacrymyces paraphysatus, but other combinations or potentially new non-European species descriptions are postponed pending further studies of additional specimens.

Citation: Zamora JC, Savchenko A, González-Cruz Á, Prieto-García F, Olariaga I, Ekman S (2022). Dendrodacrys: a new genus for species with branched hyphidia in Dacrymyces s.l., with the description of four new species. Fungal Systematics and Evolution9: 27–42. doi: 10.3114/fuse.2022.09.04

OUP accepted manuscript

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