ASAP: assemble species by automatic partitioning

Examination of the generic concept and species boundaries of the genus Erioscyphella (Lachnaceae, Helotiales, Ascomycota) with the proposal of new species and new combinations based on the Japanese materials

The genus Erioscyphella Kirschst., which was morphologically confused with Lachnum, was herein examined. Based on molecular phylogenetic analyses using a combined dataset of ITS, LSU, mtSSU, and RPB2 and morphological examinations, Erioscyphella was distinguished from Lachnum and redefined by longer ascospores and the presence of apical amorphous materials and/or resinous materials equipped on hairs. Species boundaries recognized by morphology/ecology and phylogenetic analyses were cross-checked using species delimitation analyses based on DNA barcode sequences downloaded from UNITE, resulting in that species' taxonomic problems being uncovered. Six new species (E.boninensis, E.insulae, E.otanii, E.papillaris, E.paralushanensis, and E.sasibrevispora) and two new combinations (E.hainanensis and E.sinensis) were proposed.

Re-evaluating the case for poecilogony in the gastropod Planaxis sulcatus (Cerithioidea, Planaxidae)

Background

Planaxis sulcatus has been touted as a textbook example of poecilogony, with members of this wide-ranging Indo-Pacific marine gastropod said to produce free-swimming veligers as well as brooded juveniles. A recent paper by Wiggering et al. (BMC Evol Biol 20:76, 2020) assessed a mitochondrial gene phylogeny based on partial COI and 16S rRNA sequences for 31 individuals supplemented by observations from the brood pouch of 64 mostly unsequenced individuals. ABGD and bGYMC supported three reciprocally monophyletic clades, with two distributed in the Indo-Pacific, and one restricted to the northern Indian Ocean and Red Sea. Given an apparent lack of correlation between clade membership and morphological differentiation or mode of development, the reported 3.08% maximum K2P model-corrected genetic divergence in COI among all specimens was concluded to represent population structuring. Hence, the hypothesis that phylogenetic structure is evidence of cryptic species was rejected and P. sulcatus was concluded to represent a case of geographic poecilogony.

Results

Our goal was to reassess the case for poecilogony in Planaxis sulcatus with a larger molecular dataset and expanded geographic coverage. We sequenced an additional 55 individuals and included published and unpublished sequence data from other sources, including from Wiggering et al. Our dataset comprised 108 individuals (88 COI, 81 16S rRNA) and included nine countries unrepresented in the previous study. The expanded molecular dataset yielded a maximum K2P model-corrected genetic divergence among all sequenced specimens of 12.09%. The value of 3.08% erroneously reported by Wiggering et al. is the prior maximal distance value that yields a single-species partition in ABGD, and not the maximum K2P intraspecific divergence that can be calculated for the dataset. The bGMYC analysis recognized between two and six subdivisions, while the best-scoring ASAP partitions recognized two, four, or five subdivisions, not all of which were robustly supported in Bayesian and maximum likelihood phylogenetic analyses of the concatenated and single gene datasets. These hypotheses yielded maximum intra-clade genetic distances in COI of 2.56–6.19%, which are more consistent with hypothesized species-level thresholds for marine caenogastropods.

Conclusions

Based on our analyses of a more comprehensive dataset, we conclude that the evidence marshalled by Wiggering et al. in support of Planaxis sulcatus comprising a single widespread, highly variable species with geographic poecilogony is unconvincing and requires further investigation in an integrative taxonomic framework.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-022-01961-7.

The Macrobiotus ariekammensis species complex provides evidence for parallel evolution of claw elongation in macrobiotid tardigrades

The recent integrative revision of the family Macrobiotidae demonstrated monophyly of the genus Macrobiotus and its complex, mosaic morphological evolution. Here, we analyse three Macrobiotus populations that exhibit extraordinary claw morphology characterized by elongated primary branches. Two of these populations, from the Arctic, were initially classified as Macrobiotus ariekammensis, but detailed integrative analyses resulted in splitting them into two subspecies: Macrobiotus ariekammensis ariekammensis and Macrobiotus ariekammensis groenlandicus subsp. nov.. The third population was Macrobiotus kirghizicus from Kyrgyzstan. Given the unusual phenotype of the above-mentioned taxa, we tested whether they constitute a distinct lineage in the family Macrobiotidae and could be delineated as a new genus. Although the phylogenetic investigation showed that the three taxa form a monophyletic group, the clade is nested in the genus Macrobiotus. Therefore, despite their morphological distinctiveness, a new genus cannot be established and we group these taxa in the Macrobiotus ariekammensis species complex instead. The complex includes the three above-mentioned taxa and Macrobiotus ramoli, which is included based on morphological characters. Moreover, our results provide evidence for rapid parallel evolution of long claws in macrobiotid tardigrades inhabiting cold and icy environments. Finally, we discuss the validity of the recent suppression of the genus Xerobiotus, which gathers macrobiotids with reduced claws.

A new species of Hoplias malabaricus species complex (Characiformes: Erythrinidae) from the Crepori River, Amazon basin, Brazil

A new species belonging to the Hoplias malabaricus complex from the Amazon basin, Brazil, is described. The new species is characterized by 15-16 predorsal scales, 37-39 lateral-line scales, 5 scales from dorsal fin to lateral line, 38-39 vertebrae, iii-iv, 7-8 anal-fin rays, ii-iv, 12-15 caudal-fin rays, last vertical series of scales on the base of caudal-fin rays forming a straight line, 6-7 dark bands in anal fin and no distinctive dark bands or blotches on flanks. The new species is also distinguished from other congeners of the H. malabaricus species-group by means of landmark-based morphometrics and DNA Barcoding (Cytochrome c Oxidase I gene). An identification key to species of the H. malabaricus species complex is provided.

New Sets of Primers for DNA Identification of Non-Indigenous Fish Species in the Volga-Kama Basin (European Russia)

Adequate species’ identification is critical for the detection and monitoring of biological invasions. In this study, we proposed and assessed the efficiency of newly created primer sets for the genetic identification of non-indigenous species (NIS) of fishes in the Volga basin based on: (a) a “long” fragment of cytochrome c oxidase subunit one of the mitochondrial gene (COI) (0.7 kb), used in “classical” DNA barcoding; (b) a short 3’-fragment (0.3 kb) of COI, suitable for use in high-throughput sequencing systems (i.e., for dietary analysis); (c) fragment of 16S mitochondrial rRNA, including those designed to fill the library of reference sequences for work on the metabarcoding of communities and eDNA studies; (d) a fragment of 18S nuclear rRNA, including two hypervariable regions V1-V2, valuable for animal phylogeny. All four sets of primers demonstrated a high amplification efficiency and high specificity for freshwater fish. Also, we proposed the protocols for the cost-effective isolation of total DNA and purification of the PCR product without the use of commercial kits. We propose an algorithm to carry out extremely cheap studies on the assessment of biological diversity without expensive equipment. We also present original data on the genetic polymorphism of all mass NIS fish species in the Volga-Kama region. The high efficiency of DNA identification based on our primers is shown relative to the traditional monitoring of biological invasions.

Oriental freshwater mussels arose in East Gondwana and arrived to Asia on the Indian Plate and Burma Terrane

Freshwater mussels cannot spread through oceanic barriers and represent a suitable model to test the continental drift patterns. Here, we reconstruct the diversification of Oriental freshwater mussels (Unionidae) and revise their taxonomy. We show that the Indian Subcontinent harbors a rather taxonomically poor fauna, containing 25 freshwater mussel species from one subfamily (Parreysiinae). This subfamily most likely originated in East Gondwana in the Jurassic and its representatives arrived to Asia on two Gondwanan fragments (Indian Plate and Burma Terrane). We propose that the Burma Terrane was connected with the Indian Plate through the Greater India up to the terminal Cretaceous. Later on, during the entire Paleogene epoch, these blocks have served as isolated evolutionary hotspots for freshwater mussels. The Burma Terrane collided with mainland Asia in the Late Eocene, leading to the origin of the Mekong's Indochinellini radiation. Our findings indicate that the Burma Terrane had played a major role as a Gondwanan "biotic ferry" alongside with the Indian Plate.

Using Species Groups to Approach the Large and Taxonomically Unresolved Freshwater Fish Family Nemacheilidae (Teleostei: Cypriniformes)

Large animal families with unresolved taxonomy are notoriously difficult to handle with respect to their biodiversity, systematics, and evolutionary history. We approach a large and taxonomically unresolved family of freshwater fishes (Nemacheilidae, >600 species) by proposing, on the basis of morphologic data, a species group within the family and study its phylogeny with conclusions regarding its diversity, taxonomy, and biogeographic history. Phylogenetic analyses of two mitochondrial and three nuclear genes of 139 specimens, representing about 46 species (17 candidate species from the proposed species-group, plus 29 comparative species), revealed that the proposed species group does not form a distinct monophyletic lineage, but that the candidate and comparative species mixed in three different lineages. However, the results revealed more than 20% of undescribed species within the ingroup and showed that species do not cluster according to the presently recognised genera. At least one of the genetic clades shows signs of an eastward range expansion during the second half of Miocene from north India via Myanmar into Laos, western China, and western Thailand. We conclude that the approach of picking monophyletic lineages to study biodiversity, systematics, and evolutionary history helps to open the door to large animal families.

Disentangling the Taxonomic Status of Caprella penantis sensu stricto (Amphipoda: Caprellidae) Using an Integrative Approach

Despite its importance in intertidal and shallow-water marine ecosystems, Caprella penantis continues to be one of the most taxonomically challenging amphipods in the world. A recent molecular study focusing on C. penantis sensu stricto pointed out the existence of three highly divergent lineages, indicating the possible existence of a process of ongoing speciation and, thus, casting doubt on the taxonomic status of this species. In the present study, we used an integrative approach to continue to shed light on the taxonomy and distribution of this caprellid. To this end, we combined morphological and genetic data (COI and 18S) and included, for the first time, populations from its type locality. Our analyses provide strong evidence of the existence of potentially three distinct species, genetically and geographically restricted, within C. penantis sensu stricto, with the distribution of the true C. penantis sensu stricto restricted to the UK (type locality), the northern coast of the Iberian Peninsula, and the Azores. Results show the co-occurrence of two of these species in a locality of northern Portugal and indicate the existence of distinct evolutionary and diversification patterns along the eastern Atlantic region. Overall, our study highlights the use of an integrative approach to properly assess species boundaries and unravel hidden biodiversity in amphipods.

Molecular genetic diversity of seaweeds morphologically related to Ulva rigida at three sites along the French Atlantic coast

Foliose species of the genus Ulva are notoriously difficult to identify due to their variable morphological characteristics and high phenotypic plasticity. We reassessed the taxonomic status of several distromatic foliose Ulva spp., morphologically related to Ulva rigida, using DNA barcoding with the chloroplastic tufA and rbcL (for a subset of taxa) genes for 339 selected attached Ulva specimens collected from three intertidal rocky sites. Two of the collection sites were in Brittany and one site was in Vendée, along the Atlantic coast of France. Molecular analyses included several museum specimens and the holotype of Ulva armoricana Dion, Reviers & Coat. We identified five different tufA haplotypes using a combination of phylogenetic analysis, with the support of several recently sequenced holotypes and lectotypes, and a species delimitation method based on hierarchical clustering. Four haplotypes were supported by validly named species: Ulva australis Areschoug, Ulva fenestrata Postels & Ruprecht, Ulva lacinulata (Kützing) Wittrock and U. rigida C. Agardh. The later was additionally investigated using rbcL. The fifth haplotype represented exact sequence matches to an unnamed species from European Atlantic coasts. Our results support: (1) the synonymy of both U. rigida sensu Bliding non C. Agardh and U. armoricana with U. lacinulata. This finding is based on current genetic analysis of tufA from the U. armoricana holotype and recent molecular characterization of the lectotype of U. laetevirens, which is synonymous to U. australis, (2) the presence of U. australis as a misidentified introduced species in Brittany, and (3) the presence of U. fenestrata and U. rigida in southern Brittany. The taxonomic history of each species is discussed, highlighting issues within distromatic foliose taxa of the genus Ulva and the need to genetically characterize all its available type specimens.

Patterns and Predictors of Small Mammal Phylogenetic and Functional Diversity in Contrasting Elevational Gradients in Kenya

Mountains of the Afrotropics are global biodiversity hotspots and centers of speciation and endemism; however, very few studies have focused on the phylogenetic and functional dimensions of Afromontane small mammals. We investigated the patterns and mechanisms of small mammal phylogenetic and functional diversity and assembly along elevational gradients in Mount Kenya, the second highest mountain in Africa, and a contrasting low mountain range, Chyulu Hills. We sampled 24 200-m interval transects in both sites; 18 in Mt. Kenya (9 each in the windward side, Chogoria, and the leeward side, Sirimon) and 6 in Chyulu. We extracted the mitochondrial Cytochrome b gene to reconstruct a time-calibrated species tree for estimating phylogenetic diversity indices [phylogenetic richness (PD), mean nearest taxon distance (PDMNTD), and nearest taxon index (PDNTI)]. A functional trait data set was compiled from the field-recorded measurements and published data sets for estimating functional diversity indices [functional richness (FD), mean nearest taxon distance (FDMNTD), and nearest taxon index (FDNTI)]. Several environmental variables representing water-energy availability, primary habitat productivity, and topographic heterogeneity were used to estimate the predictive power of abiotic conditions on diversity variances using generalized linear and generalized additive regression models. The PD and FD peaked around mid-elevations in Mt. Kenya, unimodally increased or decreased in Chogoria and Sirimon, and monotonically increased in Chyulu. The divergence and community structure indices—PDMNTD, FDMNTD, and PDNTI and FDNTI—were relatively weakly associated with elevation. Overall, the tendency of assemblages to be phylogenetically and functionally closely related than expected by chance decreased with elevation in Mt. Kenya but increased in Chyulu. Across the indices, the annual precipitation and topographic ruggedness were the strongest predictors in Mt. Kenya, evapotranspiration and temperature seasonality were the strongest predictors in Chyulu, while temperature seasonality and terrain ruggedness overlapped as the strongest predictors in Chogoria and Sirimon in addition to annual precipitation in the latter and normalized difference vegetation index in the former. The observed contrasting trends in diversity distribution and the strongest predictors between elevational gradients are integral to the sustainable management of the high faunal biodiversity in tropical Afromontane ecosystems.

Molecular Identification and Appraisal of the Genetic Variation of Taenia saginata in Central Regions of Vietnam

Taenia saginata is a globally distributed tapeworm responsible for human taeniasis due to the ingestion of raw or undercooked beef. T. saginata is present in several Asian countries, including China, Thailand, Lao PDR, Cambodia, and Vietnam, but little is known about its genetic variation. Studying the tapeworm's phylogeographic patterns is crucial to better understanding their association with the geographic distribution of taeniasis/cysticercosis in human populations. In the present study, 38 specimens of this putative species were collected in central regions of Vietnam and analysed using the mitochondrial gene Cytochrome c Oxidase subunit I (COI) as a molecular marker to assess the correct species identification and investigate the level of genetic variation at different geographic scales. Phylogenetic and phylogeographic analyses were carried out on a dataset that included COI sequences from Vietnamese specimens and from all conspecifics available in GenBank to date. The results showed that the collected Vietnamese specimens belonged to the species T. saginata. In Southeast Asia, signs of a possible founder effect were discovered, with the most common haplotypes frequent and present in many countries, except Lao PDR, which shares its most common haplotype only with individuals from Thailand. Remarkably, a unique taxonomic entity was found worldwide, even though the available COI sequences of T. saginata belonging to non-Asiatic countries are, at present, limited. Therefore, future studies including more COI sequences from a higher number of countries and the use of a combined molecular approach with multiple genetic markers would be useful to provide deeper insight into the global genetic variation of this species.

Multilocus marker-based delimitation of Salicornia persica and its population discrimination assisted by supervised machine learning approach

The Salicornia L. has been considered one of the most taxonomically challenging genera due to high morphological plasticity, intergradation between related species, and lack of diagnostic features in preserved herbarium specimens. In the United Arab Emirates (UAE), only one species of this genus, Salicornia europaea, has been reported, though investigating its identity at the molecular level has not yet been undertaken. Moreover, based on growth form and morphology variation between the Ras-Al-Khaimah (RAK) population and the Umm-Al-Quwain (UAQ) population, we suspect the presence of different species or morphotypes. The present study aimed to initially perform species identification using multilocus DNA barcode markers from chloroplast DNA (cpDNA) and nuclear ribosomal DNA (nrDNA), followed by the genetic divergence between two populations (RAK and UAQ) belonging to two different coastal localities in the UAE. The analysis resulted in high-quality multilocus barcode sequences subjected to species discrimination through the unsupervised OTU picking and supervised learning methods. The ETS sequence data from our study sites had high identity with the previously reported sequences of Salicornia persica using NCBI blast and was further confirmed using OTU picking methods viz., TaxonDNAs Species identifier and Assemble Species by Automatic Partitioning (ASAP). Moreover, matK sequence data showed a non-monophyletic relationship, and significant discrimination between the two populations through alignment-based unsupervised OTU picking, alignment-free Co-Phylog, and alignment & alignment-free supervised learning approaches. Other markers viz., rbcL, trnH-psbA, ITS2, and ETS could not distinguish the two populations individually, though their combination with matK (cpDNA & cpDNA+nrDNA) showed enough population discrimination. However, the ITS2+ETS (nrDNA) exhibited much higher genetic divergence, further splitting both the populations into four haplotypes. Based on the observed morphology, genetic divergence, and the number of haplotypes predicted using the matK marker, it can be suggested that two distinct populations (RAK and UAQ) do exist. Further extensive morpho-taxonomic studies are required to determine the inter-population variability of Salicornia in the UAE. Altogether, our results suggest that S. persica is the species that grow in the present study area in UAE, and do not support previous treatments as S. europaea.

Parasitism, seasonality, and diversity of trombiculid mites (Trombidiformes: Parasitengona, Trombiculidae) infesting bats (Chiroptera) in Poland

The study aims to ascertain the diversity of trombiculid species associated with Chiroptera in Poland, and for the first time in the case of research on Central European Trombiculidae, we use both DNA and morphology in an integrative taxonomic approach to determine species identities of trombiculids. The research was carried out from 2015 to 2019. In total, 2725 larvae were collected from 300 specimens of bats belonging to 11 species. Deutonymphs were obtained through laboratory rearing of larvae; few larvae and deutonymphs were collected also from bats' daily roosts. The presence of trombiculid larvae on hosts was observed between July and April of the following year, with the highest numbers recorded in autumn, during bat swarming. Male bats were infested more often than females (16.4 vs. 6.6%). The highest infestation rate was recorded for Barbastella barbastellus, Myotis nattereri and Plecotus auritus, and the highest prevalence of chiggers (> 30%) for Myotis bechsteinii and P. auritus. The larvae found on bats occupied the areas with free access to the host's skin: auricles, tragus, and snout. Morphological identification of specimens to the species level was hindered by the mosaic distribution of diagnostic traits. Morphological analyses indicated the presence of Leptotrombidium russicum and Leptotrombidium spp. in the examined material, whereas molecular analyses additionally suggested three other potential species assigned to the same genus based on the assessed scope of intrageneric variation (ASAP method). We argue that the identification of the parasitic larvae (chiggers) using morphological characters does not address the question of actual species boundaries, which, in turn, affects the inferences about host specificity and host range.

Integrative taxonomy reveals a new species of Cyphocharax (Characiformes: Curimatidae) from the Upper Paraíba do Sul River basin, Brazil

Abstract A new species of Cyphocharax is described from the Upper Paraíba do Sul River basin, São Paulo, Brazil based on integrated morphological and molecular delimitation criteria. It is morphologically distinguished from its congeners by the presence of a round, dark blotch at the midlength of the caudal peduncle not extending to the proximal portions of the median caudal-fin rays, 19–20 circumpeduncular scales, 34–41 perforated lateral-line scales, 6–7 longitudinal scale rows above and below the lateral line, greatest body depth corresponding to 34.7–39.9% of standard length (SL), and the caudal peduncle depth corresponding to 13.3–15.2% of SL. The lowest genetic distances between the new species and other congeners are: 2.5% from C. gilbert, followed by 3.0% from C. santacatarinae, and 3.2% from C. aff. gilbert. All species delimitation criteria employed herein corroborated the recognition of the new species. In addition, comments on its conservation status are provided.

Variation in the Myotis keaysi complex (Chiroptera, Vespertilionidae), with description of a new species from Ecuador

The genus Myotis comprises a diverse group of vespertilionid bats with worldwide distribution. Twenty-eight Neotropical species are currently recognized. In this study, we evaluate molecular and morphological variation in the M. keaysi complex, a high elevation clade of Neotropical myotine bats characterized by complex taxonomy and high morphological variation. A phylogeny inferred with cytochrome-b sequences recovered two clades composed of samples traditionally assigned to M. keaysi, with 9% of genetic divergence between them. These clades were also suggested as putative distinct species by molecular species delimitation methods. Qualitative and quantitative morphological analyses indicated a phenotypic discontinuity between specimens from central Andes (including the holotype of M. keaysi) and western lowlands of Ecuador, showing strong congruence between molecular and morphological approaches. We therefore describe a new species for the Tumbes-Choco-Magdalena region, documenting their external and cranial diagnostic characters by comparing them with other Neotropical species. In addition, we provide an emended diagnosis for our new concept of M. keaysi.

Complete Genome Sequence, Molecular Characterization and Phylogenetic Relationships of a Novel Tern Atadenovirus

Discovery and study of viruses carried by migratory birds are tasks of high importance due to the host's ability to spread infectious diseases over significant distances. With this paper, we present and characterize the first complete genome sequence of atadenovirus from a tern bird (common tern, Sterna hirundo) preliminarily named tern atadenovirus 1 (TeAdV-1). TeAdV-1 genome is a linear double-stranded DNA molecule, 31,334 base pairs which contain 30 methionine-initiated open reading frames with gene structure typical for Atadenovirus genus, and the shortest known inverted terminal repeats (ITRs) within the Atadenovirus genus consisted of 25 bases. The nucleotide composition of the genome is characterized by a low G + C content (33.86%), which is the most AT-rich genome of known avian adenoviruses within Atadenovirus genus. The nucleotide sequence of the TeAdV-1 genome shows high divergence compared to known representatives of the Atadenovirus genus with the highest similarity to the duck atadenovirus 1 (53.7%). Phylogenetic analysis of the protein sequences of core genes confirms the taxonomic affiliation of the new representative to the genus Atadenovirus with the degree of divergence from the known representatives exceeding the interspecies distance within the genus. Thereby we proposed a novel TeAdV-1 to be considered as a separate species.

Phylogenetic analyses of a new freshwater amphipod reveal polyphyly within the Holarctic family Crangonyctidae, with revision of the genus Synurella

A new genus and species of crangonyctid amphipod, Sicifera cahawba gen. & sp. nov., is described from Dallas County, AL, USA, based on both morphological and molecular comparison with similar crangonyctids. These data, with the application of four species delimitation models, identify the taxon as distinct when compared with related species. Nearctic members of the crangonyctid genus Synurella form a separate, well-supported monophyletic lineage when compared with Palaearctic members, which differ considerably in both molecular and morphological markers. Nearctic members, with the exception of the enigmatic Synurella (Eosynurella) johanseni, are placed in the newly erected Sicifera. The separation of these two genera implies that Palaearctic and Nearctic crangonyctid lineages might not be as closely related as once thought, and their evolutionary and biogeographical history requires further review. In addition, a key to Nearctic members of the genera Eosynurella/Sicifera is presented to aid in future identification.

Molecular Analysis Reveals a High Diversity of Anopheline Mosquitoes in Yanomami Lands and the Pantanal Region of Brazil

Identifying the species of the subfamily Anophelinae that are Plasmodium vectors is important to vector and malaria control. Despite the increase in cases, vector mosquitoes remain poorly known in Brazilian indigenous communities. This study explores Anophelinae mosquito diversity in the following areas: (1) a Yanomami reserve in the northwestern Amazon Brazil biome and (2) the Pantanal biome in southwestern Brazil. This is carried out by analyzing cytochrome c oxidase (COI) gene data using Refined Single Linkage (RESL), Assemble Species by Automatic Partitioning (ASAP), and tree-based multi-rate Poisson tree processes (mPTP) as species delimitation approaches. A total of 216 specimens collected from the Yanomami and Pantanal regions were sequenced and combined with 547 reference sequences for species delimitation analyses. The mPTP analysis for all sequences resulted in the delimitation of 45 species groups, while the ASAP analysis provided the partition of 48 groups. RESL analysis resulted in 63 operational taxonomic units (OTUs). This study expands our scant knowledge of anopheline species in the Yanomami and Pantanal regions. At least 18 species of Anophelinae mosquitoes were found in these study areas. Additional studies are now required to determine the species that transmit Plasmodium spp. in these regions.

DNA barcodes evidence the contact zone of eastern and western caddisfly lineages in the Western Carpathians

The region of the Western Carpathians is, among other aspects, very important for survival and diversity of European freshwater fauna due to the presence of a large number of (sub)mountain springs and streams. However, these ecologically and faunistically diversified habitats are still understudied in the context of genetic diversity and population structure of their inhabitants. This study focuses on genetic diversity and distribution patterns of the caddisfly Rhyacophila tristis, common and widespread representative of mountain freshwater fauna. Analysis of the COI mitochondrial marker revealed presence of the western and eastern lineages, with samples from both lineages being grouped in BOLD (Barcode of Life Data System) into separate BINs (Barcode Index Numbers). Our data indicates that eastern lineage (BIN_E) is more closely related to the Balkan populations than to co-occurring western lineage (BIN_W), and that the contact zone of the lineages passes through the W Carpathians. The study revealed phylogeographic and demographic differences between lineages, supporting hypothesis of their evolutionary independence and specific ecological preferences. The obtained genetic data of the R. tristis population from W Carpathians improved our knowledge about population genetics of this aquatic species and can contribute to understanding the state and evolution of biodiversity of freshwater ecosystems in Europe.

Distribution of Stenasellidae in Africa and description of a new species of Metastenasellus from Cameroonian groundwaters

During recent investigations of the groundwater fauna of Cameroon, specimens of a new species of the stygobitic genus Metastenasellus, M. boutini sp. nov. were collected in wells of the city of Douala. The new species can be easily distinguished from the other species of the genus by its relatively large size (up to 11 mm), pleonite 1 and 2 half the length of pereonite 7, the shape of pleopod 2 in males (presence of an external lobe on the protopodite, distal part of the spermatic duct slightly protruding out of the second article, lack of a distal seta on the exopodite), and uropod half the length of the pleotelson. Ecological data and a key to Metastenasellus species are provided. We also performed an exhaustive analysis of the literature on Stenasellidae in Africa to study the geographical distribution of the family in this continent and discuss some hypotheses about the origin of African species.

G. C. Genetic identification of marine eels (Anguilliformes: Congroidei) through DNA barcoding from Kasimedu fishing harbour

Along with the mysteries of their body's shape like snakes, marine eels have fascinated biologists for centuries. Information on the molecular taxonomy of marine eels is scarce from the Southeast Indian region and hence, the present study aimed to barcode marine eels collected from Kasimedu fishing harbor, Chennai, Tamil Nadu. A total of 44 specimens were collected and DNA barcoding was done with a COI marker. The evolutionary history was inferred using the BA method. We observed 17 species, 10 genera, 4 families from the suborder Congroidei of which the genus Ariosoma and Conger were found to be predominant. The species of the family Muraenesocidae and Congridae are highly variable. The average Kimura two-parameter (K2P) distances within species, genera, and families were 3.08%, 6.80%, 13.80%, respectively. Maximum genetic distance (0.307) was observed between the species Muraenesox cinereus and Ariosoma sp.1. BA tree topology revealed distinct clusters in concurrence with the taxonomic status of the species. A deeper split was observed in Uroconger lepturus. We sequenced for the first-time barcode of Sauromuraenesox vorax and a new species Ophichthus chennaiensis is the gap-filling in identifying this taxon in the Indian context. We found a correct match between morphological and genetic identification of the species analyzed, depending on the cluster analysis performed (BINs and ASAP). This demonstrates that the COI gene sequence is suitable for phylogenetic analysis and species identification.

Species boundaries in the messy middle-A genome-scale validation of species delimitation in a recently diverged lineage of coastal fog desert lichen fungi

Species delimitation among closely related species is challenging because traditional phenotype-based approaches, for example, using morphology, ecological, or chemical characteristics, may not coincide with natural groupings. With the advent of high-throughput sequencing, it has become increasingly cost-effective to acquire genome-scale data which can resolve previously ambiguous species boundaries. As the availability of genome-scale data has increased, numerous species delimitation analyses, such as BPP and SNAPP+Bayes factor delimitation (BFD*), have been developed to delimit species boundaries. However, even empirical molecular species delimitation approaches can be biased by confounding evolutionary factors, for example, hybridization/introgression and incomplete lineage sorting, and computational limitations. Here, we investigate species boundaries and the potential for micro-endemism in a lineage of lichen-forming fungi, Niebla Rundel & Bowler, in the family Ramalinaceae by analyzing single-locus and genome-scale data consisting of (a) single-locus species delimitation analysis using ASAP, (b) maximum likelihood-based phylogenetic tree inference, (c) genome-scale species delimitation models, e.g., BPP and SNAPP+BFD, and (d) species validation using the genealogical divergence index (gdi). We specifically use these methods to cross-validate results between genome-scale and single-locus datasets, differently sampled subsets of genomic data and to control for population-level genetic divergence. Our species delimitation models tend to support more speciose groupings that were inconsistent with traditional taxonomy, supporting a hypothesis of micro-endemism, which may include morphologically cryptic species. However, the models did not converge on robust, consistent species delimitations. While the results of our analysis are somewhat ambiguous in terms of species boundaries, they provide a valuable perspective on how to use these empirical species delimitation methods in a nonmodel system. This study thus highlights the challenges inherent in delimiting species, particularly in groups such as Niebla, with complex, relatively recent phylogeographic histories.

A comprehensive phylogeny of dwarf geckos of the genus Lygodactylus, with insights into their systematics and morphological variation

The 71 currently known species of dwarf geckos of the genus Lygodactylus are a clade of biogeographic interest due to their occurrence in continental Africa, Madagascar, and South America. Furthermore, because many species are morphologically cryptic, our knowledge of species-level diversity within this genus is incomplete, as indicated by numerous unnamed genetic lineages revealed in previous molecular studies. Here we provide an extensive multigene phylogeny covering 56 of the named Lygodactylus species, four named subspecies, and 34 candidate species of which 19 are newly identified in this study. Phylogenetic analyses, based on ∼10.1 kbp concatenated sequences of eight nuclear-encoded and five mitochondrial gene fragments, confirm the monophyly of 14 Lygodactylus species groups, arranged in four major clades. We recover two clades splitting from basal nodes, one comprising exclusively Malagasy species groups, and the other containing three clades. In the latter, there is a clade with only Madagascar species, which is followed by a clade containing three African and one South American species groups, and its sister clade containing six African and two Malagasy species groups. Relationships among species groups within these latter clades remain weakly supported. We reconstruct a Lygodactylus timetree based on a novel fossil-dated phylotranscriptomic tree of squamates, in which we included data from two newly sequenced Lygodactylus transcriptomes. We estimate the crown diversification of Lygodactylus started at 46 mya, and the dispersal of Lygodactylus among the main landmasses in the Oligocene and Miocene, 35-22 mya, but emphasize the wide confidence intervals of these estimates. The phylogeny suggests an initial out-of-Madagascar dispersal as most parsimonious, but accounting for poorly resolved nodes, an out-of-Africa scenario may only require one extra dispersal step. More accurate inferences into the biogeographic history of these geckos will likely require broader sampling of related genera and phylogenomic approaches to provide better topological support. A survey of morphological characters revealed that most of the major clades and species groups within Lygodactylus cannot be unambiguously characterized by external morphology alone, neither by unique character states nor by a diagnostic combination of character states. Thus, any future taxonomic work will likely benefit from integrative, phylogenomic approaches.

The importance of being integrative: a remarkable case of synonymy in the genus Viridiscus (Heterotardigrada: Echiniscidae)

There are two predominant sources of taxonomically useful morphological variability in the diverse tardigrade family Echiniscidae: the internal structure and surface sculpture of the cuticular plates covering the dorsum (sculpturing) and the arrangement and morphology of the trunk appendages (chaetotaxy). However, since the appendages often exhibit intraspecific variation (they can be reduced or can develop asymmetrically), sculpturing has been considered more stable at the species level and descriptions of new echiniscid species based solely on morphology are still being published. Here, we present a case study in which a detailed analysis of the morphology and multiple genetic markers of several species of the genus Viridiscus shows that cuticular sculpture may also exhibit considerable intraspecific variation and lead to false taxonomic conclusions. In a population collected from the eastern Nearctic, in the type locality of the recently described species V. miraviridis, individuals with transitional morphotypes between those reported for V. viridissimus and V. miraviridis were found. Importantly, all morphotypes within the viridissimus-miraviridis spectrum were grouped in a single monospecific clade according to rapidly evolving markers (ITS-1, ITS-2 and COI). Given the morphological and genetic evidence, we establish V. miraviridis as a junior synonym of V. viridissimus. This study explicitly demonstrates that a lack of DNA data associated with morphological descriptions of new taxa jeopardizes the efforts to unclutter tardigrade systematics. Additionally, V. perviridis and V. viridissimus are reported from Lâm Đồng Province in southern Vietnam, which considerably broadens their known geographic ranges.

Elevational Diversity Patterns of Green Lacewings (Neuroptera: Chrysopidae) Uncovered With DNA Barcoding in a Biodiversity Hotspot of Southwest China

Elevational diversity patterns can reflect the responses of biodiversity to climate change spatially. We investigate the species diversity patterns of green lacewings (an important predatory group of insects) along the gradient of elevation from the Shaluli Mountains (Mts. Shaluli), which belong to the Hengduan Mountains in southwestern China, one of the important hotspots of global biodiversity. We combined multiple approaches, including Automatic Barcode Gap Discovery (ABGD), Assemble Species by Automatic Partitioning analysis (ASAP), General Mixed Yule Coalescent (GMYC), Poisson tree processes (bPTP), multi-rate Poisson tree processes (mPTP), to delimit the green lacewings species based on the standard barcoding region of cytochrome c oxidase subunit I (COI). The α-diversity and β-diversity patterns of green lacewings from the Mts. Shaluli along the gradient of elevation were analyzed, with further exploration on how the temperature effect elevational-diversity pattern on broad-scale (county scale) elevational gradients. The DNA barcoding reference library consisted of 40 green lacewing species from the Mts. Shaluli. The α-diversity of green lacewings decreased with the increasing elevation. The temperature was found to have a significant effect on the abundance and Shannon-Wiener diversity index but not on the species richness. Nestedness replaced turnover as the main component of Sørensen’s dissimilarity with the increasing elevation, and greater nestedness occurred at low temperature areas. The combination of a reliable DNA barcoding database could improve the accuracy and efficiency to investigate the species diversity patterns of green lacewings. Temperature, resource, and resultant interspecific competitions may have important roles in explaining the species diversity patterns of green lacewings from the Mts. Shaluli. Priority of conservation should be given to the species at low elevation, middle elevation, and relatively high temperature regions under the background of global climate warming.

First report of mitochondrial COI in foraminifera and implications for DNA barcoding

Foraminifera are a species-rich phylum of rhizarian protists that are highly abundant in many marine environments and play a major role in global carbon cycling. Species recognition in Foraminifera is mainly based on morphological characters and nuclear 18S ribosomal RNA barcoding. The 18S rRNA contains variable sequence regions that allow for the identification of most foraminiferal species. Still, some species show limited variability, while others contain high levels of intragenomic polymorphisms, thereby complicating species identification. The use of additional, easily obtainable molecular markers other than 18S rRNA will enable more detailed investigation of evolutionary history, population genetics and speciation in Foraminifera. Here we present the first mitochondrial cytochrome c oxidase subunit 1 (COI) gene sequences ("barcodes") of Foraminifera. We applied shotgun sequencing to single foraminiferal specimens, assembled COI, and developed primers that allow amplification of COI in a wide range of foraminiferal species. We obtained COI sequences of 49 specimens from 17 species from the orders Rotaliida and Miliolida. Phylogenetic analysis showed that the COI tree is largely congruent with previously published 18S rRNA phylogenies. Furthermore, species delimitation with ASAP and ABGD algorithms showed that foraminiferal species can be identified based on COI barcodes.

Molecular species delimitation refines the taxonomy of native and nonnative physinine snails in North America

Being able to associate an organism with a scientific name is fundamental to our understanding of its conservation status, ecology, and evolutionary history. Gastropods in the subfamily Physinae have been especially troublesome to identify because morphological variation can be unrelated to interspecific differences and there have been widespread introductions of an unknown number of species, which has led to a speculative taxonomy. To resolve uncertainty about species diversity in North America, we targeted an array of single-locus species delimitation methods at publically available specimens and new specimens collected from the Snake River basin, USA to generate species hypotheses, corroborated using nuclear analyses of the newly collected specimens. A total-evidence approach delineated 18 candidate species, revealing cryptic diversity within recognized taxa and a lack of support for other named taxa. Hypotheses regarding certain local endemics were confirmed, as were widespread introductions, including of an undescribed taxon likely belonging to a separate genus in southeastern Idaho for which the closest relatives are in southeast Asia. Overall, single-locus species delimitation was an effective first step toward understanding the diversity and distribution of species in Physinae and to guiding future investigation sampling and analyses of species hypotheses.

Fuzzy species borders of glacial survivalists in the Carpathian biodiversity hotspot revealed using a multimarker approach

The Carpathians are one of the key biodiversity hotspots in Europe. The mountain chain uplifted during Alpine orogenesis and is characterised by a complex geological history. Its current biodiversity was highly influenced by Pleistocene glaciations. The goal of the current study was to examine the phylogenetic and demographic history of Gammarus balcanicus species complex in the Carpathians using multiple markers as well as to delimit, using an integrative approach, and describe new species hidden so far under the name G. balcanicus. Results showed that divergence of the studied lineages reaches back to the Miocene, which supports the hypothesis of their survival in multiple micro refugia. Moreover, the increase of their diversification rate in the Pleistocene suggests that glaciation was the driving force of their speciation. The climatic changes during and after the Pleistocene also played a major role in the demography of the local Carpathian lineages. Comparison of diversity patterns and phylogenetic relationships of both, the mitochondrial and nuclear markers, provide evidence of putative hybridisation and retention of ancient polymorphism (i.e., incomplete lineage sorting). The morphological examination supported the existence of two morphological types; one we describe as a G. stasiuki sp. nov. and another we redescribe as a G. tatrensis (S. Karaman, 1931).

Sea surface temperature, rather than land mass or geographic distance, may drive genetic differentiation in a species complex of highly dispersive seabirds

Seabirds, particularly Procellariiformes, are highly mobile organisms with a great capacity for long dispersal, though simultaneously showing high philopatry, two conflicting life-history traits that may lead to contrasted patterns of genetic population structure. Landmasses were suggested to explain differentiation patterns observed in seabirds, but philopatry, isolation by distance, segregation between breeding and nonbreeding zones, and oceanographic conditions (sea surface temperatures) may also contribute to differentiation patterns. To our knowledge, no study has simultaneously contrasted the multiple factors contributing to the diversification of seabird species, especially in the gray zone of speciation. We conducted a multilocus phylogeographic study on a widespread seabird species complex, the little shearwater complex, showing highly homogeneous morphology, which led to considerable taxonomic debate. We sequenced three mitochondrial and six nuclear markers on all extant populations from the Atlantic (lherminieri) and Indian Oceans (bailloni), that is, five nominal lineages from 13 populations, along with one population from the eastern Pacific Ocean (representing the dichrous lineage). We found sharp differentiation among populations separated by the African continent with both mitochondrial and nuclear markers, while only mitochondrial markers allowed characterizing the five nominal lineages. No differentiation could be detected within these five lineages, questioning the strong level of philopatry showed by these shearwaters. Finally, we propose that Atlantic populations likely originated from the Indian Ocean. Within the Atlantic, a stepping-stone process accounts for the current distribution. Based on our divergence time estimates, we suggest that the observed pattern of differentiation mostly resulted from historical and current variation in sea surface temperatures.

Genetic structure of the protist Physarum albescens (Amoebozoa) revealed by multiple markers and genotyping by sequencing

Myxomycetes are terrestrial protists with many presumably cosmopolitan species dispersing via airborne spores. A truly cosmopolitan species would suffer from outbreeding depression hampering local adaptation, while locally adapted species with limited distribution would be at a higher risk of extinction in changing environments. Here, we investigate intraspecific genetic diversity and phylogeography of Physarum albescens over the entire Northern Hemisphere. We sequenced 324 field collections of fruit bodies for 1-3 genetic markers (SSU, EF1A, COI) and analysed 98 specimens with genotyping by sequencing. The structure of the three-gene phylogeny, SNP-based phylogeny, phylogenetic networks, and the observed recombination pattern of three independently inherited gene markers can be best explained by the presence of at least 18 reproductively isolated groups, which can be seen as cryptic species. In all intensively sampled regions and in many localities, members of several phylogroups coexisted. Some phylogroups were found to be abundant in only one region and completely absent in other well-studied regions, and thus may represent regional endemics. Our results demonstrate that the widely distributed myxomycete species Ph. albescens represents a complex of at least 18 cryptic species, and some of these seem to have a limited geographical distribution. In addition, the presence of groups of presumably clonal specimens suggests that sexual and asexual reproduction coexist in natural populations of myxomycetes.

DNA barcoding for species delimitation of the freshwater leech genus Glossiphonia from the Western Balkan (Hirudinea, Glossiphoniidae)

Glossiphoniid leeches are a diverse group and sometimes abundant elements of the aquatic fauna inhabiting various types of freshwater habitats. In this study, we sampled leeches of the genus Glossiphonia from the Western Balkan in order to test the suitability of the mitochondrial cytochrome c oxidase subunit 1 (COI) marker sequence for species delimitation. Morphological analysis revealed the presence of four taxa, G.complanata with two subspecies, G.c.complanata and G.c.maculosa, the latter an endemic of Ohrid Lake, G.nebulosa and endemic G.balcanica. In total, 29 new barcodes of Glossiphonia were sequenced in the course of this study and compared with the available molecular dataset of the latter genus from GenBank/BOLD databases. The applied ASAP distance-based species delimitation method for the analysed dataset revealed an interspecific threshold between 4-8% K2P distance as suitable for species identification purposes of the Western Balkan Glossiphonia species. Our study revealed that morphologically identified taxa as G.nebulosa and G.concolor each consists of more than one clearly different phylogenetic clade. This study contributes to a better knowledge of the taxonomy of glossiphoniid leeches and emphasises future work on the revision of this genus using a standard molecular COI marker in species identification.

The dark side of pseudoscorpion diversity: The German Barcode of Life campaign reveals high levels of undocumented diversity in European false scorpions

DNA barcoding is particularly useful for identification and species delimitation in taxa with conserved morphology. Pseudoscorpions are arachnids with high prevalence of morphological crypsis. Here, we present the first comprehensive DNA barcode library for Central European Pseudoscorpiones, covering 70% of the German pseudoscorpion fauna (35 out of 50 species). For 21 species, we provide the first publicly available COI barcodes, including the rare Anthrenochernes stellae Lohmander, a species protected by the FFH Habitats Directive. The pattern of intraspecific COI variation and interspecific COI variation (i.e., presence of a barcode gap) generally allows application of the DNA barcoding approach, but revision of current taxonomic designations is indicated in several taxa. Sequences of 36 morphospecies were assigned to 74 BINs (barcode index numbers). This unusually high number of intraspecific BINs can be explained by the presence of overlooked cryptic species and by the accelerated substitution rate in the mitochondrial genome of pseudoscorpions, as known from previous studies. Therefore, BINs may not be an appropriate proxy for species numbers in pseudoscorpions, while partitions built with the ASAP algorithm (Assemble Species by Automatic Partitioning) correspond well with putative species. ASAP delineated 51 taxonomic units from our data, an increase of 42% compared with the present taxonomy. The Neobisium carcionoides complex, currently considered a polymorphic species, represents an outstanding example of cryptic diversity: 154 sequences from our dataset were allocated to 23 BINs and 12 ASAP units.

A re-analysis of the data in Sharkey et al.'s (2021) minimalist revision reveals that BINs do not deserve names, but BOLD Systems needs a stronger commitment to open science

Halting biodiversity decline is one of the most critical challenges for humanity, but monitoring biodiversity is hampered by taxonomic impediments. One impediment is the large number of undescribed species (here called "dark taxon impediment") whereas another is caused by the large number of superficial species descriptions, that can only be resolved by consulting type specimens ("superficial description impediment"). Recently, Sharkey et al. (2021) proposed to address the dark taxon impediment for Costa Rican braconid wasps by describing 403 species based on COI barcode clusters ("BINs") computed by BOLD Systems. More than 99% of the BINs (387 of 390) were converted into species by assigning binominal names (e.g. BIN "BOLD:ACM9419" becomes Bracon federicomatarritai) and adding a minimal diagnosis (consisting only of a consensus barcode for most species). We here show that many of Sharkey et al.'s species are unstable when the underlying data are analyzed using different species delimitation algorithms. Add the insufficiently informative diagnoses, and many of these species will become the next "superficial description impediment" for braconid taxonomy because they will have to be tested and redescribed after obtaining sufficient evidence for confidently delimiting species. We furthermore show that Sharkey et al.'s approach of using consensus barcodes as diagnoses is not functional because it cannot be applied consistently. Lastly, we reiterate that COI alone is not suitable for delimiting and describing species, and voice concerns over Sharkey et al.'s uncritical use of BINs because they are calculated by a proprietary algorithm (RESL) that uses a mixture of public and private data. We urge authors, reviewers and editors to maintain high standards in taxonomy by only publishing new species that are rigorously delimited with open-access tools and supported by publicly available evidence.

Alpine lichen diversity in an isolated sky island in the Colorado Plateau, USA-Insight from an integrative biodiversity inventory

Lichens are major components of high altitude/latitude ecosystems. However, accurately characterizing their biodiversity is challenging because these regions and habitats are often underexplored, there are numerous poorly known taxonomic groups, and morphological variation in extreme environments can yield conflicting interpretations. Using an iterative taxonomic approach based on over 800 specimens and incorporating both traditional morphology-based identifications and information from the standard fungal DNA barcoding marker, we compiled a voucher-based inventory of biodiversity of lichen-forming fungi in a geographically limited and vulnerable alpine community in an isolated sky island in the Colorado Plateau, USA-the La Sal Mountains. We used the newly proposed Assemble Species by Automatic Partitioning (ASAP) approach to empirically delimit candidate species-level lineages from family-level multiple sequence alignments. Specimens comprising DNA-based candidate species were evaluated using traditional taxonomically diagnostic phenotypic characters to identify specimens to integrative species hypotheses and link these, where possible, to currently described species. Despite the limited alpine habitat (ca. 3,250 ha), we document the most diverse alpine lichen community known to date from the southern Rocky Mountains, with up to 240 candidate species/species-level lineages of lichen-forming fungi. 139 species were inferred using integrative taxonomy, plus an additional 52 candidate species within 29 different putative species complexes. Over 68% of sequences could not be assigned to species-level rank with statistical confidence, corroborating the limited utility of current sequence repositories for species-level DNA barcoding of lichen-forming fungi. By integrating vouchered specimens, DNA sequence data, and photographic documentation, we provide an important baseline of lichen-forming fungal diversity for the limited alpine habitat in the Colorado Plateau. These data provide an important resource for subsequent research in the ecology and evolution of lichens alpine habitats, including DNA barcodes for most putative species/species-level lineages occurring in the La Sal Mountains, and vouchered collections representing any potentially undescribed species that can be used for future taxonomic studies.

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.