HaplowebMaker and CoMa: Two web tools to delimit species using haplowebs and conspecificity matrices

Hapsolutely: a user-friendly tool integrating haplotype phasing, network construction, and haploweb calculation

Motivation

Haplotype networks are a routine approach to visualize relationships among alleles. Such visual analysis of single-locus data is still of importance, especially in species diagnosis and delimitation, where a limited amount of sequence data usually are available and sufficient, along with other datasets in the framework of integrative taxonomy. In diploid organisms, this often requires separating (phasing) sequences with heterozygotic positions, and typically separate programs are required for phasing, reformatting of input files, and haplotype network construction. We therefore developed Hapsolutely, a user-friendly program with an ergonomic graphical user interface that integrates haplotype phasing from single-locus sequences with five approaches for network/genealogy reconstruction.

Results

Among the novel options implemented, Hapsolutely integrates phasing and graphical reconstruction steps of haplotype networks, supports input of species partition data in the common SPART and SPART-XML formats, and calculates and visualizes haplowebs and fields for recombination, thus allowing graphical comparison of allele distribution and allele sharing among subsets for the purpose of species delimitation. The new tool has been specifically developed with a focus on the workflow in alpha-taxonomy, where exploring fields for recombination across alternative species partitions may help species delimitation.

Availability and implementation

Hapsolutely is written in Python, and integrates code from Phase, SeqPHASE, and PopART in C++ and Haxe. Compiled stand-alone executables for MS Windows and Mac OS along with a detailed manual can be downloaded from https://www.itaxotools.org; the source code is openly available on GitHub (https://github.com/iTaxoTools/Hapsolutely).

Speciation on the Roof of the World: Parallel Fast Evolution of Cryptic Mole Vole Species in the Pamir-Alay-Tien Shan Region

Speciation is not always accompanied by morphological changes; numerous cryptic closely related species were revealed using genetic methods. In natural populations of Ellobius tancrei (2n = 54-30) and E. alaicus (2n = 52-48) of the Pamir-Alay and Tien Shan, the chromosomal variability due to Robertsonian translocations has been revealed. Here, by comprehensive genetic analysis (karyological analyses as well as sequencing of mitochondrial genes, cytb and COI, and nuclear genes, XIST and IRBP) of E. alaicus and E. tancrei samples from the Inner Tien Shan, the Alay Valley, and the Pamir-Alay, we demonstrated fast and independent diversification of these species. We described an incompletely consistent polymorphism of the mitochondrial and nuclear markers, which arose presumably because of habitat fragmentation in the highlands, rapid karyotype changes, and hybridization of different intraspecific varieties and species. The most intriguing results are a low level of genetic distances calculated from mitochondrial and nuclear genes between some phylogenetic lines of E. tancrei and E. alaicus, as well significant species-specific chromosome variability in both species. The chromosomal rearrangements are what most clearly define species specificity and provide further diversification. The "mosaicism" and inconsistency in polymorphism patterns are evidence of rapid speciation in these mammals.

Surrounded by Kindred: Spermophilus major Hybridization with Other Spermophilus Species in Space and Time

Among the numerous described cases of hybridization in mammals, the most intriguing are (a) cases of introgressive hybridization deeply affecting the evolutionary history of species, and (b) models involving not a pair of species but a multi-species complex. Therefore, the hybridization history of the russet ground squirrel Spermophilus major, whose range has repeatedly changed due to climatic fluctuations and now borders the ranges of four related species, is of great interest. The main aims of this study were to determine the direction and intensity of gene introgression, the spatial depth of the infiltration of extraneous genes into the S. major range, and to refine the hypothesis of the hybridogenic replacement of mitochondrial genomes in the studied group. Using phylogenetic analysis of the variability of mitochondrial (CR, cytb) and nuclear (SmcY, BGN, PRKCI, c-myc, i6p53) markers, we determined the contribution of neighboring species to the S. major genome. We showed that 36% of S. major individuals had extraneous alleles. All peripheral species that were in contact with S. major contributed towards its genetic variability. We also proposed a hypothesis for the sequence and localization of serial hybridization events. Our assessment of the S. major genome implications of introgression highlights the importance of implementing conservation measures to protect this species.

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

Background

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

Methods

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

Results

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

Geographic Mosaic of Extensive Genetic Variations in Subterranean Mole Voles Ellobius alaicus as a Consequence of Habitat Fragmentation and Hybridization

Restricted mobility, sociality, and high inbreeding—characteristic for subterranean mammals—lead to rapid changes in their genome structure. Up to now, the Alay mole vole Ellobius alaicus was a data-deficient species; its spatial and phylogenetic relationships with a sibling species, E. tancrei, were not clarified. We carried out a genetic analysis including differential G-banding of chromosomes and mitochondrial (cytb) and nuclear gene (XIST and IRBP) sequencing. The phylogenetic reconstruction based on cytb represented the expected phylogenetic relationships of two species. Using the XIST, we revealed two new lineages among E. alaicus from the Alay Valley (Southern Kyrgyzstan). Analysis of IRBP demonstrated presence of the specific genotype in most of E. alaicus specimens, but also revealed the haplotype, typical for E. tancrei, in some Alay mole voles. The results may be explained as persistence of ancestral gene polymorphism in E. alaicus or limited interspecific hybridization with E. tancrei. Several chromosomal forms were revealed in E. alaicus in the Alay Valley. We propose that ‘mosaic’ genetic polymorphism might appear in E. alaicus due to fragmentation of their habitats in highlands of the Alay Valley, Tien Shan, and Pamir-Alay as well as due to hybridization with E. tancrei or persistence of ancestral polymorphisms.

Morphological and molecular diversity patterns of the genus Tropodiaptomus Kiefer, 1932 (Copepoda, Calanoida, Diaptomidae) in Thailand

Tropodiaptomus is one of the most specious genera in the family Diaptomidae, but it is often rare in terms of distribution and abundance. Moreover, Tropodiaptomus species show a noteworthy variability in some of the morphological characters considered of prime importance in diaptomid taxonomy, and the presence of cryptic or pseudocryptic species is likely. Thus, through a geographically-wide sampling in Thailand, we aimed to investigate the local diversity of the genus and to compare the morphological and molecular diversity pattern based on mitochondrial and nuclear genes sequences. DNA taxonomy was also implemented in order to check whether the Tropodiaptomus lineages were independent species according to the “evolutionary genetic species concept”. Six Tropodiaptomus morphospecies were found, three of which are putative species new to Science pending a formal description. The finding of such a high incidence of undescribed species stresses the existence of a significant “Linnean shortfall” affecting Thai diaptomids. The molecular results showed that most of the studied species could be identified consistently with their morphology-based taxonomy. However, Tropodiaptomus vicinus and T. cf. lanaonus showed a high level of genetic diversity, suggesting that traditional morphological techniques might be inadequate for correctly assessing their taxonomical status.

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).

Rediscovery and redescription of Niphargus enslini Karaman, 1932 (Amphipoda, Niphargidae) in southern Germany

Niphargus enslini Karaman, 1932 was collected only once in 1905 from the Falkensteiner Höhle (Baden-Württemberg, Germany). Two years after its description, the species was synonymized with Niphargus virei and not studied any more. During recent surveys on German niphargids, further samples collected in this cave did not yield N. enslini specimens, but this species was collected in the Blätterteighöhle and in the Schwarzer Brunnen, two caves located in Baden-Württemberg and intercepting the same karstic aquifer feeding Falkensteiner Höhle. In an integrative taxonomic approach, we carefully studied the morphology of the newly collected specimens and sequenced two molecular markers (fragments of the cytochrome c oxidase subunit I (COI) and of the nuclear 28S rRNA gene) to test for possible conspecificity of N. enslini with N. virei. Morphological analysis confirmed that N. enslini is distinct from the N. virei species complex. We provide a redescription of newly collected material, together with new drawings of a more than 100 years old topotypic female. We briefly discuss the putative origin of N. enslini and the age of its split from the N. virei species complex.

Solving the Coral Species Delimitation Conundrum

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

Evolutionary history of Sundaland shrews (Eulipotyphla: Soricidae: Crocidura) with a focus on Borneo

The hyperdiverse shrew genus Crocidura is one of few small mammal genera distributed across Sundaland and all of its boundaries. This represents a rare opportunity to study the geological history of this region through the evolutionary history of these shrews. We generate a phylogeny of all recognized species of Sundaland Crocidura and show that most speciation events took place during the Pleistocene, prior to the inundation of the Sunda Shelf around 400 000 years ago. We find east–west differentiation within two separate lineages on Borneo, and that the current taxonomy of its two endemic species does not reflect evolutionary history, but ecophenotypic variation of plastic traits related to elevation. Sulawesi shrews are monophyletic, with a single notable exception: the black-footed shrew (C. nigripes). We show that the black-footed shrew diverged from its relatives on Borneo recently, suggesting a human-assisted breach of Wallace’s line. Overall, the number of Crocidura species, especially on Borneo, probably remains an underestimate.

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.

A Revised Phylogeny of the Mentha spicata Clade Reveals Cryptic Species

The genus Mentha is taxonomically and phylogenetically challenging due to complex genomes, polyploidization and an extensive historical nomenclature, potentially hiding cryptic taxa. A straightforward interpretation of phylogenetic relationships within the section Mentha is further hindered by dominant but outdated concepts on historically identified hybrid taxa. Mentha spicata is traditionally considered to be of hybrid origin, but the evidence for this is weak. Here, we aim to understand the phylogenetic relationships within the section Mentha using large sample sizes and to revisit the hybrid status and identity of M. spicata. We show that two of three traditional species in the subsection Spicatae are polyphyletic, as is the subsection as a whole, while the real number of cryptic species was underestimated. Compared to previous studies we present a fundamentally different phylogeny, with a basal split between M. spicata s.s. and M. longifolia s.s. Cluster analyses of morphological and genotypic data demonstrate that there is a dissociation between morphologically and genotypically defined groups of samples. We did not find any evidence that M. spicata is of hybrid origin, and we conclude its taxonomic status should be revised. The combination of genetic and phenotypic information is essential when evaluating hyperdiverse taxonomic groups.

Five animal phyla in glacier ice reveal unprecedented biodiversity in New Zealand's Southern Alps

Glacier ice is an extreme environment in which most animals cannot survive. Here we report the colonization of high elevation, climate-threatened glaciers along New Zealand's southwestern coast by species of Arthropoda, Nematoda, Platyhelminthes, Rotifera and Tardigrada. Based on DNA barcoding and haplotype-inferred evidence for deep genetic variability, at least 12 undescribed species are reported, some of which have persisted in this niche habitat throughout the Pleistocene. These findings identify not only an atypical biodiversity hotspot but also highlight the adaptive plasticity of microinvertebrate Animalia.

Population genetics of the brooding coral Seriatopora hystrix reveals patterns of strong genetic differentiation in the Western Indian Ocean

Coral reefs provide essential goods and services but are degrading at an alarming rate due to local and global anthropogenic stressors. The main limitation that prevents the implementation of adequate conservation measures is that connectivity and genetic structure of populations are poorly known. Here, the genetic diversity and connectivity of the brooding scleractinian coral Seriatopora hystrix were assessed at two scales by genotyping ten microsatellite markers for 356 individual colonies. S. hystrix showed high differentiation, both at large scale between the Red Sea and the Western Indian Ocean (WIO), and at smaller scale along the coast of East Africa. As such high levels of differentiation might indicate the presence of more than one species, a haploweb analysis was conducted with the nuclear marker ITS2, confirming that the Red Sea populations are genetically distinct from the WIO ones. Based on microsatellite analyses three groups could be distinguished within the WIO: (1) northern Madagascar, (2) south-west Madagascar together with one site in northern Mozambique (Nacala) and (3) all other sites in northern Mozambique, Tanzania and Kenya. These patterns of restricted connectivity could be explained by the short pelagic larval duration of S. hystrix, and/or by oceanographic factors, such as eddies in the Mozambique Channel (causing larval retention in northern Madagascar but facilitating dispersal from northern Mozambique towards south-west Madagascar). This study provides an additional line of evidence supporting the conservation priority status of the Northern Mozambique Channel and should inform coral reef management decisions in the region.

Testing species hypotheses for Fridericia magna, an enchytraeid worm (Annelida: Clitellata) with great mitochondrial variation

BACKGROUND

Deep mitochondrial divergences were observed in Scandinavian populations of the terrestrial to semi-aquatic annelid Fridericia magna (Clitellata: Enchytraeidae). This raised the need for testing whether the taxon is a single species or a complex of cryptic species.

RESULTS

A total of 62 specimens from 38 localities were included in the study, 44 of which were used for species delimitation. First, the 44 specimens were divided into clusters using ABGD (Automatic Barcode Gap Discovery) on two datasets, consisting of sequences of the mitochondrial markers COI and 16S. For each dataset, the worms were divided into six not completely congruent clusters. When they were combined, a maximum of seven clusters, or species hypotheses, were obtained, and the seven clusters were used as input in downstream analyses. We tested these hypotheses by constructing haplowebs for two nuclear markers, H3 and ITS, and in both haplowebs the specimens appeared as a single species. Multi-locus species delimitation analyses performed with the Bayesian BPP program also mainly supported a single species. Furthermore, no apparent morphological differences were found between the clusters. Two of the clusters were partially separated from each other and the other clusters, but not strongly enough to consider them as separate species. All 62 specimens were used to visualise the Scandinavian distribution, of the species, and to compare with published COI data from other Fridericia species.

CONCLUSION

We show that the morphospecies Fridericia magna is a single species, harbouring several distinct mitochondrial clusters. There is partial genetic separation between some of them, which may be interpreted as incipient speciation. The study shows the importance of rigorous species delimitation using several independent markers when deep mitochondrial divergences might give the false impression of cryptic speciation.