Identifying spatially concordant evolutionary significant units across multiple species through DNA barcodes: Application to the conservation genetics of the freshwater fishes of Java and Bali

Phylogeography of Microphis retzii (Bleeker, 1856) and Microphis brachyurus (Bleeker, 1854) in the Pacific

Syngnathidae are a charismatic family of teleost fishes, represented by seahorses, seadragons, and pipefishes. Syngnathidae are mainly composed of marine species, but about 30 species of pipefishes inhabit freshwater insular environments of the Indo-Pacific realm. Recent research has shown that some freshwater pipefish species are amphidromous and exhibit high intraspecific divergences across their distribution range, like Microphis brachyurus (Bleeker, 1854) distributed from Sri Lanka to French Polynesia and Microphis retzii (Bleeker, 1856) distributed from Taiwan to Indonesia. In this study, we used the mitochondrial cytochrome oxidase I (COI) partial gene of 91 specimens of M. brachyurus and 30 specimens of M. retzii from localities representative of their respective distribution area to improve knowledge on the population structure of these two widespread species. Genetic species delimitation and phylogeographic analyses were combined to explore spatial patterns of genetic diversity across the distribution ranges of the two species. We have highlighted deep genetic structuring within the two species and relate these results to various biotic and abiotic factors. For M. brachyurus, the population in Polynesia is distinct from those in the West Pacific, suggesting its distinctiveness and recognition as an evolutionary significant unit (ESU). For M. retzii, three lineages are delimited in its range distribution, suggesting the existence of two distinct species in Southeast Asia (Bali/Java/Lombok and China/Taiwan). Pipefish species are particularly vulnerable to anthropogenic pressures (inherent to Syngnathidae and insular environments). The present results, revising species delimitation and geographic distribution, will help implement effective conservation and management measures.

Contrasting genetic and morphological differentiation among geographical lineages of a stenotopic miniature rasborine, Boraras maculatus, in Peninsular Malaysia

The variability in the stenotopic miniature rasborine Boraras maculatus (Cypriniformes: Danionidae: Rasborinae) across acidic-water habitats of Peninsular Malaysia (PM) was investigated using two molecular markers (the mitochondrial cytochrome c oxidase subunit I [COI] gene and the nuclear rhodopsin gene), as well as morphological evidence. Molecular phylogenetic analyses revealed differentiation among populations of B. maculatus in PM with the distinction of four allopatric lineages. Each of them was recognized as a putative species by automatic species delimitation methods. These lineages diverged from each other between 7.4 and 1.9 million years ago. A principal component analysis (PCA) was conducted to examine the multivariate variation in 11 morphometric measurements among three of these lineages. PCA results showed a significant overlap in morphological characteristics among these lineages. Additionally, a photograph-based machine learning approach failed to fully differentiate these lineages, suggesting limited morphological differentiation. B. maculatus represents a case of morphological stasis in a stenotopic miniature species. Strong habitat preference, coupled with long-term habitat fragmentation, may explain why each lineage of B. maculatus has a restricted distribution and did not disperse to other regions within and outside of PM, despite ample possibilities when the Sunda shelf was emerged and drained by large paleodrainages for most of the past 7 million years. The conservation status of B. maculatus and its peat swamp habitats are discussed, and it is concluded that peat swamps comprise several evolutionary units. Each of these units is considered a conservation unit and deserves appropriate protection.

Cryptic diversity impacts model selection and macroevolutionary inferences in diversification analyses

Species persist in landscapes through ecological dynamics but proliferate at wider spatial scales through evolutionary mechanisms. Disentangling the contribution of each dynamic is challenging, but the increasing use of dated molecular phylogenies opened new perspectives. First, the increasing use of DNA sequences in biodiversity inventory shed light on a substantial amount of cryptic diversity in species-rich ecosystems. Second, explicit diversification models accounting for various eco-evolutionary models are now available. Integrating both advances, we explored diversification trajectories among 10 lineages of freshwater fishes in Sundaland, for which time-calibrated and taxonomically rich phylogenies are available. By fitting diversification models to dated phylogenies and incorporating DNA-based species delimitation methods, the impact of cryptic diversity on diversification model selection and related inferences is explored. Eight clades display constant speciation rate model as the most likely if cryptic diversity is accounted, but nine display a signature of diversification slowdowns when cryptic diversity is ignored. Cryptic diversification occurs during the last 5 Myr for most groups, and palaeoecological models received little support. Most cryptic lineages display restricted range distribution, supporting geographical isolation across homogeneous landscapes as the main driver of diversification. These patterns question the persistence of cryptic diversity and its role during species proliferation.

Specific and Intraspecific Diversity of Symphypleona and Neelipleona (Hexapoda: Collembola) in Southern High Appalachia (USA)

Collembola, commonly known as springtails, are important detritivores, abundant in leaf litter and soil globally. Springtails are wingless hexapods with many North American species having wide distributions ranging from as far as Alaska to Mexico. Here, we analyze the occurrence and intraspecific diversity of springtails with a globular body shape (Symphypleona and Neelipleona), in southern high Appalachia, a significant biodiversity hotspot. The peaks of high Appalachia represent ‘sky islands’ due to their physical isolation, and they host numerous endemic species in other taxa. We surveyed globular Collembola through COI metabarcoding, assessing geographic and genetic diversity across localities and species. Intraspecific diversity in globular Collembola was extremely high, suggesting that considerable cryptic speciation has occurred. While we were able to associate morphospecies with described species in most of the major families in the region (Dicyrtomidae, Katiannidae, Sminthuridae, and Sminthurididae), other families (Neelidae, and Arrhopalitidae) are in more pressing need of taxonomic revision before species identities can be confirmed. Due to poor representation in databases, and high intraspecific variability, no identifications were accomplished through comparison with available DNA barcodes.

Pleistocene allopatric differentiation followed by recent range expansion explains the distribution and molecular diversity of two congeneric crustacean species in the Palaearctic

Pleistocene glaciations had a tremendous impact on the biota across the Palaearctic, resulting in strong phylogeographic signals of range contraction and rapid postglacial recolonization of the deglaciated areas. Here, we explore the diversity patterns and history of two sibling species of passively dispersing taxa typical of temporary ponds, fairy shrimps (Anostraca). We combine mitochondrial (COI) and nuclear (ITS2 and 18S) markers to conduct a range-wide phylogeographic study including 56 populations of Branchinecta ferox and Branchinecta orientalis in the Palaearctic. Specifically, we investigate whether their largely overlapping ranges in Europe resulted from allopatric differentiation in separate glacial refugia followed by a secondary contact and reconstruct their postglacial recolonization from the inhabited refugia. Our results suggest the existence of distinct refugia for the two species, with genetic divergence among intraspecific lineages consistent with late Pleistocene glacial cycles. While B. ferox lineages originated from Mediterranean refugia, the origin of B. orientalis lineages was possibly located on the Pannonian Plain. We showed that most dispersal events predominantly happened within 100 km, coupled with several recent long-distance events (> 1000 km). Hence the regional habitat density of suitable habitats in Central Europe is possibly a key to the co-existence of the two species. Overall, our study illustrates how isolation in combination with stochastic effects linked to glacial periods are important drivers of the allopatric differentiation of Palaearctic taxa.

To lump, to split or to maintain? Molecular taxonomy of the spotted barb Barbodes binotatus (Cyprinidae) and closely related species in Peninsular Malaysia

The taxonomic status of the Southeast Asian spotted barb, Barbodes binotatus (Teleostei: Cyprinidae), has puzzled researchers because of large but inconsistent geographic variation of its body melanin marking pattern. In this study, the authors appraise the differentiation of B. binotatus and two closely related species, Barbodes rhombeus and saddle barb, Barbodes banksi, in Peninsular Malaysia using mitochondrial and nuclear markers. The results of this study reveal that the Peninsular Malaysia populations of each of the three species form largely reciprocal monophyletic lineages that differ from each other by a minimum of 2.3% p-genetic distance using COI gene. Nonetheless, specimens of B. binotatus in Peninsular Malaysia are only distantly related to specimens of B. binotatus in Java (type locality). The monophyly of B. banksi is not refuted although specimens of Peninsular Malaysia are genetically distinct from those of Sarawak (type locality). The authors discuss alternative hypotheses whether each of these three valid species is a single species or each of the main five genetic lineages revealed in this study represents a distinct species. Preliminary investigations reveal a mito-nuclear discordance at one locality in Peninsular Malaysia where B. binotatus and B. banksi co-occur. Further studies should inform on the extent of reproductive porousness between these two lineages and others.

Aquatic Organisms Research with DNA Barcodes

Since their inception, DNA barcodes have become a powerful tool for understanding the biodiversity and biology of aquatic species, with multiple applications in diverse fields such as food security, fisheries, environmental DNA, conservation, and exotic species detection. Nevertheless, most aquatic ecosystems, from marine to freshwater, are understudied, with many species disappearing due to environmental stress, mostly caused by human activities. Here we highlight the progress that has been made in studying aquatic organisms with DNA barcodes, and encourage its further development in assisting sustainable use of aquatic resources and conservation.

Exploring the vertebrate fauna of the Bird's Head Peninsula (Indonesia, West Papua) through DNA barcodes

Biodiversity knowledge is widely heterogeneous across the Earth's biomes. Some areas, due to their remoteness and difficult access, present large taxonomic knowledge gaps. Mostly located in the tropics, these areas have frequently experienced a fast development of anthropogenic activities during the last decades and are therefore of high conservation concerns. The biodiversity hotspots of Southeast Asia exemplify the stakes faced by tropical countries. While the hotspots of Sundaland (Java, Sumatra, Borneo) and Wallacea (Sulawesi, Moluccas) have long attracted the attention of biologists and conservationists alike, extensive parts of the Sahul area, in particular the island of New Guinea, have been much less explored biologically. Here, we describe the results of a DNA-based inventory of aquatic and terrestrial vertebrate communities, which was the objective of a multidisciplinary expedition to the Bird's Head Peninsula (West Papua, Indonesia) conducted between 17 October and 20 November 2014. This expedition resulted in the assembly of 1005 vertebrate DNA barcodes. Based on the use of multiple species-delimitation methods (GMYC, PTP, RESL, ABGD), 264 molecular operational taxonomic units (MOTUs) were delineated, among which 75 were unidentified and an additional 48 were considered cryptic. This study suggests that the diversity of vertebrates of the Bird's Head is severely underestimated and considerations on the evolutionary origin and taxonomic knowledge of these biotas are discussed.

Conservation genetics of the yellow-bellied toad (Bombina variegata): population structure, genetic diversity and landscape effects in an endangered amphibian

Revealing patterns of genetic diversity and barriers for gene flow are key points for successful conservation in endangered species. Methods based on molecular markers are also often used to delineate conservation units such as evolutionary significant units and management units. Here we combine phylo-geographic analyses (based on mtDNA) with population and landscape genetic analyses (based on microsatellites) for the endangered yellow-bellied toad Bombina variegata over a wide distribution range in Germany. Our analyses show that two genetic clusters are present in the study area, a northern and a southern/central one, but that these clusters are not deeply divergent. The genetic data suggest high fragmentation among toad occurrences and consequently low genetic diversity. Genetic diversity and genetic connectivity showed a negative relationship with road densities and urban areas surrounding toad occurrences, indicating that these landscape features act as barriers to gene flow. To preserve a maximum of genetic diversity, we recommend considering both genetic clusters as management units, and to increase gene flow among toad occurrences with the aim of restoring and protecting functional meta-populations within each of the clusters. Several isolated populations with especially low genetic diversity and signs of inbreeding need particular short-term conservation attention to avoid extinction. We also recommend to allow natural gene flow between both clusters but not to use individuals from one cluster for translocation or reintroduction into the other. Our results underscore the utility of molecular tools for species conservation, highlight outcomes of habitat fragmentation onto the genetic structure of an endangered amphibian and reveal particularly threatened populations in need for urgent conservation efforts.

Impact of Pleistocene Eustatic Fluctuations on Evolutionary Dynamics in Southeast Asian Biodiversity Hotspots

Pleistocene Climatic Fluctuations (PCF) are frequently highlighted as important evolutionary engines that triggered cycles of biome expansion and contraction. While there is ample evidence of the impact of PCF on biodiversity of continental biomes, the consequences in insular systems depend on the geology of the islands and the ecology of the taxa inhabiting them. The idiosyncratic aspects of insular systems are exemplified by the islands of the Sunda Shelf in Southeast Asia (Sundaland), where PCF-induced eustatic fluctuations had complex interactions with the geology of the region, resulting in high species diversity and endemism. Emergent land in Southeast Asia varied drastically with sea level fluctuations during the Pleistocene. Climate-induced fluctuations in sea level caused temporary connections between insular and continental biodiversity hotspots in Southeast Asia. These exposed lands likely had freshwater drainage systems that extended between modern islands: the Paleoriver Hypothesis. Built upon the assumption that aquatic organisms are among the most suitable models to trace ancient river boundaries and fluctuations of landmass coverage, the present study aims to examine the evolutionary consequences of PCF on the dispersal of freshwater biodiversity in Southeast Asia. Time-calibrated phylogenies of DNA-delimited species were inferred for six species-rich freshwater fish genera in Southeast Asia (Clarias, Channa, Glyptothorax, Hemirhamphodon, Dermogenys, Nomorhamphus). The results highlight rampant cryptic diversity and the temporal localization of most speciation events during the Pleistocene, with 88% of speciation events occurring during this period. Diversification analyses indicate that sea level-dependent diversification models poorly account for species proliferation patterns for all clades excepting Channa. Ancestral area estimations point to Borneo as the most likely origin for most lineages, with two waves of dispersal to Sumatra and Java during the last 5 Myr. Speciation events are more frequently associated with boundaries of the paleoriver watersheds, with 60%, than islands boundaries, with 40%. In total, one-third of speciation events are inferred to have occured within paleorivers on a single island, suggesting that habitat heterogeneity and factors other than allopatry between islands substantially affected diversification of Sundaland fishes. Our results suggest that species proliferation in Sundaland is not wholly reliant on Pleistocene sea-level fluctuations isolating populations on different islands.

DNA barcoding reveals cryptic diversity in the underestimated genus Triplophysa (Cypriniformes: Cobitidae, Nemacheilinae) from the northeastern Qinghai-Tibet Plateau

Background

The northeastern part of the Qinghai-Tibet Plateau (QTP) presents a high number of plateau loach species. As one of the three major groups of fishes distributed on the QTP, plateau loach has high ecological value. However, the taxonomy and systematics of these fish are still controversial, and a large number of new species have been reported. The reason for this phenomenon is that the degree of morphological variation is low, the phylogenetic information provided by morphological and anatomical features used for species identification is relatively poor, and many cryptic species are observed. Based on the high-density sampling points from the biodiversity hotspots surveyed, this study aims to evaluate the biodiversity of plateau loach in the northeastern part of the QTP and reveal the hidden diversity by comparing morphological species with molecular operational taxonomic units (MOTUs).

Results

After careful identification and comparison of the morphology and DNA barcoding of 1630 specimens, 22 species were identified, with 20 considered valid local species and two identified as new species that had not been previously described. Based on the combination of morphological and molecular methods, a total of 24 native species were found, two of which were cryptic species: Triplophysa robusta sp1 and Triplophysa minxianensis sp1. Fourteen of the 24 species form clusters of barcodes that allow them to be reliably identified. The remaining cases involved 10 closely related species, including rapidly differentiated species and species that seemed to have experienced incomplete lineage sorting or showed introgressions.

Conclusions

The results highlight the need to combine traditional taxonomies with molecular methods to correctly identify species, especially closely related species, such as the plateau loach. This study provides a basis for protecting the biodiversity of plateau loach.

Niche divergence of evolutionarily significant units with implications for repopulation programs of the world's largest amphibians

The niche divergence and potential climate change-induced loss of evolutionarily significant units (ESUs) of flagship amphibian species in China, the Chinese giant salamander clade, were investigated. We tested niche-related ecological hypotheses and identified suitable habitats that are essential for the conservation of ESUs in response to future climate change according to ecological niche models (ENMs). We predicted the localized habitat loss crisis of ESUs induced by global climate heating using the predicted climate derived from two representative concentration pathway (RCP) scenarios 2.6 and 8.5, respectively. In our study, a niche conservatism pattern was found between the two distinctive northern and southern ESUs with sufficient distributional records, but their niches were not equivalent. Furthermore, there was neither abrupt environmental change in nor remarkable biogeographic barriers between the suitable habitats of the species, as indicated by random linear, blob and ribbon range-breaking tests. Under the low-emission scenario RCP2.6, the northern ESU had a moderate loss of suitable range, while the southern ESU had range expansion in the 2070s. The climatic velocities were low in the ranges of both ESUs. However, under the high-emission scenario RCP8.5, the climatic velocities were found to become larger in the suitable ranges of both ESUs. Moreover, the northern ESU had severe habitat loss, bringing it to the edge of extinction, while the southern ESU also had intensified range loss. Considering this, climatic velocity can be an effective indicator of range loss. We argued conclusively that conservation prioritization of ESUs should effectively take into account the underlying geographic and ecological mechanisms driving the speciation process. The conservation of ESUs should consider the conservation of both evolutionary potential and ecological adaptation capacity of each lineage. The present study provided practical guidelines for repopulation programs for endangered species and the conservation of evolutionary diversity.

Disentangling the taxonomy of the subfamily Rasborinae (Cypriniformes, Danionidae) in Sundaland using DNA barcodes

Sundaland constitutes one of the largest and most threatened biodiversity hotspots; however, our understanding of its biodiversity is afflicted by knowledge gaps in taxonomy and distribution patterns. The subfamily Rasborinae is the most diversified group of freshwater fishes in Sundaland. Uncertainties in their taxonomy and systematics have constrained its use as a model in evolutionary studies. Here, we established a DNA barcode reference library of the Rasborinae in Sundaland to examine species boundaries and range distributions through DNA-based species delimitation methods. A checklist of the Rasborinae of Sundaland was compiled based on online catalogs and used to estimate the taxonomic coverage of the present study. We generated a total of 991 DNA barcodes from 189 sampling sites in Sundaland. Together with 106 previously published sequences, we subsequently assembled a reference library of 1097 sequences that covers 65 taxa, including 61 of the 79 known Rasborinae species of Sundaland. Our library indicates that Rasborinae species are defined by distinct molecular lineages that are captured by species delimitation methods. A large overlap between intraspecific and interspecific genetic distance is observed that can be explained by the large amounts of cryptic diversity as evidenced by the 166 Operational Taxonomic Units detected. Implications for the evolutionary dynamics of species diversification are discussed.