Discordance between morphological and molecular species boundaries among Caribbean species of the reef sponge Callyspongia

Sponges are among the most species‐rich and ecologically important taxa on coral reefs, yet documenting their diversity is difficult due to the simplicity and plasticity of their morphological characters. Genetic attempts to identify species are hampered by the slow rate of mitochondrial sequence evolution characteristic of sponges and some other basal metazoans. Here we determine species boundaries of the Caribbean coral reef sponge genus Callyspongia using a multilocus, model‐based approach. Based on sequence data from one mitochondrial (COI), one ribosomal (28S), and two single‐copy nuclear protein‐coding genes, we found evolutionarily distinct lineages were not concordant with current species designations in Callyspongia. While C. fallax, C. tenerrima, and C. plicifera were reciprocally monophyletic, four taxa with different morphologies (C. armigera, C. longissima, C. eschrichtii, and C. vaginalis) formed a monophyletic group and genetic distances among these taxa overlapped distances within them. A model‐based method of species delimitation supported collapsing these four into a single evolutionary lineage. Variation in spicule size among these four taxa was partitioned geographically, not by current species designations, indicating that in Callyspongia, these key taxonomic characters are poor indicators of genetic differentiation. Taken together, our results suggest a complex relationship between morphology and species boundaries in sponges.     

Integrative species delimitation reveals cryptic diversity in the southern Appalachian Antrodiaetus unicolor (Araneae: Antrodiaetidae) species complex

Although species delimitation can be highly contentious, the development of reliable methods to accurately ascertain species boundaries is an imperative step in cataloguing and describing Earth's quickly disappearing biodiversity. Spider species delimitation remains largely based on morphological characters; however, many mygalomorph spider populations are morphologically indistinguishable from each other yet have considerable molecular divergence. The focus of our study, the Antrodiaetus unicolor species complex containing two sympatric species, exhibits this pattern of relative morphological stasis with considerable genetic divergence across its distribution. A past study using two molecular markers, COI and 28S, revealed that A. unicolor is paraphyletic with respect to A. microunicolor. To better investigate species boundaries in the complex, we implement the cohesion species concept and use multiple lines of evidence for testing genetic exchangeability and ecological interchangeability. Our integrative approach includes extensively sampling homologous loci across the genome using a RADseq approach (3RAD), assessing population structure across their geographic range using multiple genetic clustering analyses that include structure , principal components analysis and a recently developed unsupervised machine learning approach (Variational Autoencoder). We evaluate ecological similarity by using large‐scale ecological data for niche‐based distribution modelling. Based on our analyses, we conclude that this complex has at least one additional species as well as confirm species delimitations based on previous less comprehensive approaches. Our study demonstrates the efficacy of genomic‐scale data for recognizing cryptic species, suggesting that species delimitation with one data type, whether one mitochondrial gene or morphology, may underestimate true species diversity in morphologically homogenous taxa with low vagility.     

A genomic evaluation of taxonomic trends through time in coast horned lizards (genus Phrynosoma)

Determining the boundaries between species and deciding when to describe new species are challenging practices that are particularly difficult in groups with high levels of geographic variation. The coast horned lizards (Phrynosoma blainvillii, Phrynosoma cerroense and P. coronatum) have an extensive geographic distribution spanning many distinctive ecological regions ranging from northern California to the Cape Region of Baja California, Mexico, and populations differ substantially with respect to external morphology across much of this range. The number of taxa recognized in the group has been reevaluated by herpetologists over 20 times during the last 180 years, and typically without the aid of explicit species delimitation methods, resulting in a turbulent taxonomy containing anywhere from one to seven taxa. In this study, we evaluate taxonomic trends through time by ranking 15 of these species delimitation models (SDMs) using coalescent analyses of nuclear loci and SNPs in a Bayesian model comparison framework. Species delimitation models containing more species were generally favoured by Bayesian model selection; however, several three‐species models outperformed some four‐ and five‐species SDMs, and the top‐ranked model, which contained five species, outperformed all SDMs containing six species. Model performance peaked in the 1950s based on marginal likelihoods estimated from nuclear loci and SNPs. Not surprisingly, SDMs based on genetic data outperformed morphological taxonomies when using genetic data alone to evaluate models. The de novo estimation of population structure favours a three‐population model that matches the currently recognized integrative taxonomy containing three species. We discuss why Bayesian model selection might favour models containing more species, and why recognizing more than three species might be warranted.     

Integrating genomic and phenotypic data to evaluate alternative phylogenetic and species delimitation hypotheses in a recent evolutionary radiation of grasshoppers

Although resolving phylogenetic relationships and establishing species limits are primary goals of systematics, these tasks remain challenging at both conceptual and analytical levels. Here, we integrated genomic and phenotypic data and employed a comprehensive suite of coalescent‐based analyses to develop and evaluate competing phylogenetic and species delimitation hypotheses in a recent evolutionary radiation of grasshoppers (Chorthippus binotatus group) composed of two species and eight putative subspecies. To resolve the evolutionary relationships within this complex, we first evaluated alternative phylogenetic hypotheses arising from multiple schemes of genomic data processing and contrasted genetic‐based inferences with different sources of phenotypic information. Second, we examined the importance of number of loci, demographic priors, number and kind of phenotypic characters and sex‐based trait variation for developing alternative species delimitation hypotheses. The best‐supported topology was largely compatible with phenotypic data and showed the presence of two clades corresponding to the nominative species groups, one including three well‐resolved lineages and the other comprising a four‐lineage polytomy and a well‐differentiated sister taxon. Integrative species delimitation analyses indicated that the number of employed loci had little impact on the obtained inferences but revealed the higher power provided by an increasing number of phenotypic characters and the usefulness of assessing their phylogenetic information content and differences between sexes in among‐taxa trait variation. Overall, our study highlights the importance of integrating multiple sources of information to test competing phylogenetic hypotheses and elucidate the evolutionary history of species complexes representing early stages of divergence where conflicting inferences are more prone to appear.     

Sky island diversification meets the multispecies coalescent – divergence in the spruce‐fir moss spider (Microhexura montivaga,Araneae, Mygalomorphae) on the highest peaks of southern Appalachia

Microhexura montivaga is a miniature tarantula‐like spider endemic to the highest peaks of the southern Appalachian mountains and is known only from six allopatric, highly disjunct montane populations. Because of severe declines in spruce‐fir forest in the late 20th century, M. montivaga was formally listed as a US federally endangered species in 1995. Using DNA sequence data from one mitochondrial and seven nuclear genes, patterns of multigenic genetic divergence were assessed for six montane populations. Independent mitochondrial and nuclear discovery analyses reveal obvious genetic fragmentation both within and among montane populations, with five to seven primary genetic lineages recovered. Multispecies coalescent validation analyses [guide tree and unguided Bayesian Phylogenetics and Phylogeography (BPP), Bayes factor delimitation (BFD)] using nuclear‐only data congruently recover six or seven distinct lineages; BFD analyses using combined nuclear plus mitochondrial data favour seven or eight lineages. In stark contrast to this clear genetic fragmentation, a survey of secondary sexual features for available males indicates morphological conservatism across montane populations. While it is certainly possible that morphologically cryptic speciation has occurred in this taxon, this system may alternatively represent a case where extreme population genetic structuring (but not speciation) leads to an oversplitting of lineage diversity by multispecies coalescent methods. Our results have clear conservation implications for this federally endangered taxon and illustrate a methodological issue expected to become more common as genomic‐scale data sets are gathered for taxa found in naturally fragmented habitats.     

A new integrative framework for large‐scale assessments of biodiversity and community dynamics,using littoral gastropods and crabs of British Columbia,Canada

Improving our understanding of species responses to environmental changes is an important contribution ecologists can make to facilitate effective management decisions. Novel synthetic approaches to assessing biodiversity and ecosystem integrity are needed, ideally including all species living in a community and the dynamics defining their ecological relationships. Here, we present and apply an integrative approach that links high‐throughput, multicharacter taxonomy with community ecology. The overall purpose is to enable the coupling of biodiversity assessments with investigations into the nature of ecological interactions in a community‐level data set. We collected 1195 gastropods and crabs in British Columbia. First, the General mixed Yule‐coalescent (GMYC) and the Poisson Tree Processes (PTP) methods for proposing primary species‐hypotheses based on cox1 sequences were evaluated against an integrative taxonomic framework. We then used data on the geographic distribution of delineated species to test species co‐occurrence patterns for nonrandomness using community‐wide and pairwise approaches. Results showed that PTP generally outperformed GMYC and thus constitutes a more effective option for producing species‐hypotheses in community‐level data sets. Nonrandom species co‐occurrence patterns indicative of ecological relationships or habitat preferences were observed for grazer gastropods, whereas assemblages of carnivorous gastropods and crabs appeared influenced by random processes. Species‐pair associations were consistent with current ecological knowledge, thus suggesting that applying community assembly within a large taxonomical framework constitutes a valuable tool for assessing ecological interactions. Combining phylogenetic, morphological and co‐occurrence data enabled an integrated view of communities, providing both a conceptual and pragmatic framework for biodiversity assessments and investigations into community dynamics.     

Delimiting species in the genus Otospermophilus (Rodentia: Sciuridae), using genetics,ecology, and morphology

We apply an integrative taxonomy approach to delimit species of ground squirrels in the genus Otospermophilus because the diverse evolutionary histories of organisms shape the existence of taxonomic characters. Previous studies of mitochondrial DNA from this group recovered three divergent lineages within Otospermophilus beecheyi separated into northern, central, and southern geographical populations, with Otospermophilus atricapillus nested within the southern lineage of O. beecheyi. To further evaluate species boundaries within this complex, we collected additional genetic data (one mitochondrial locus, 11 microsatellite markers, and 11 nuclear loci), environmental data (eight bioclimatic variables), and morphological data (23 skull measurements). We used the maximum number of possible taxa (O. atricapillus, Northern O. beecheyi, Central O. beecheyi, and Southern O. beecheyi) as our operational taxonomic units (OTUs) and examined patterns of divergence between these OTUs. Phenotypic measures (both environmental and morphological) showed little differentiation among OTUs. By contrast, all genetic datasets supported the evolutionary independence of Northern O. beecheyi, although they were less consistent in their support for other OTUs as distinct species. Based on these data, we support the conclusions from a previous study that synonymized O. atricapillus with O. beecheyi, and we elevate the northern lineage of O. beecheyi to a separate species. © 2014 The Linnean Society of London, Biological Journal of the Linnean Society, 2014, 113 , 1136–1151.     

Multivariate species boundaries and conservation of harlequin poison frogs

In this study, we present an iterative method for delimiting species under the general lineage concept (GLC) based on the multivariate clustering of morphological, ecological and genetic data. Our rationale is that distinct multivariate groups correspond to evolutionarily independent metapopulation lineages because they reflect the common signal of different secondary defining properties (environmental and genetic distinctiveness, phenotypic diagnosability, etc.) that imply the existence of barriers preventing or limiting gene exchange. We applied this method to study a group of endangered poison frogs, the Oophaga histrionica complex. In our study case, we used next‐generation targeted amplicon sequencing to obtain a robust genetic data set that we combined with patterns of morphological and ecological features. Our analyses revealed the existence of at least five different species in the histrionica complex (three, new to science), some of them, occurring in small isolated populations outside any protected areas. The lineage delimitation proposed here has important conservation implications as it revealed that some of the Oophaga species should be considered among the most vulnerable of the Neotropical frogs. More broadly, our study exemplifies how multiple‐amplicon and multivariate statistical techniques can be integrated to successfully identify species and their boundaries.     

Phylogeny,phylogeography and geographical variation in the Crocidura monax (Soricidae) species complex from the montane islands of Tanzania,with descriptions of three new species

We assess morphological and multilocus genetic variation among 11 isolated montane populations of white‐toothed shrews from Tanzania that have been referred to either Crocidura monax Thomas or C. montis Thomas. The montane sites we sampled represent ‘sky‐islands’ from two geologically distinct archipelagos (Northern Highlands and the Eastern Arc Mountains) and are a significant component of the Eastern Afromontane Biodiversity Hotspot. We used multivariate analyses of morphometric traits and phylogenetic and species‐delimitation analyses of multilocus DNA sequence data to assess species‐level diversity. Our species delimitation analyses included a novel, pairwise validation approach that avoids potential biases associated with specifying a guide tree. These analyses reveal several distinct lineages, which we treat as six allopatric species. Each species is restricted to one, two or four mountains. We use available names to recognize C. monax, C. tansaniana Hutterer and C. usambarae Dippenaar, while naming and describing three new species. Our results demonstrate the effectiveness of combining morphological and genetic data to uncover and describe hidden diversity in a cryptic mammalian system. © 2015 The Linnean Society of London     

Integrative taxonomy of the primitively segmented spider genus Ganthela (Araneae: Mesothelae: Liphistiidae): DNA barcoding gap agrees with morphology

Species delimitation is difficult for taxa in which the morphological characters are poorly known because of the rarity of adult morphs or sexes, and in cryptic species. In primitively segmented spiders, family Liphistiidae, males are often unknown, and female genital morphology – usually species‐specific in spiders – exhibits considerable intraspecific variation. Here, we report on an integrative taxonomic study of the liphistiid genus Ganthela Xu & Kuntner, 2015, endemic to south‐east China, where males are only available for two of the seven morphological species (two known and five undescribed). We obtained DNA barcodes (cytochrome c oxidase subunit I gene, COI) for 51 newly collected specimens of six morphological species and analysed them using five species‐delimitation methods: DNA barcoding gap, species delimitation plugin [P ID(Liberal)], automatic barcode gap discovery (ABGD), generalized mixed Yule‐coalescent model (GMYC), and statistical parsimony (SP). Whereas the first three agreed with the morphology, GMYC and SP indicate several additional species. We used the consensus results to delimit and diagnose six Ganthela species, which in addition to the type species Ganthela yundingensis Xu, 2015, completes the revision of the genus. Although multi‐locus phylogenetic approaches may be needed for complex taxonomic delimitations, our results indicate that even single‐locus analyses based on the COI barcodes, if integrated with morphological and geographical data, may provide sufficiently reliable species delimitation. © 2015 The Linnean Society of London     

Assessing the diversity and distribution of Cephalothrix species (Nemertea: Palaeonemertea) in European waters by comparing different species delimitation methods

Soft‐bodied marine taxa, like ribbon worms (Nemertea), often lack clear diagnostic morphological characters impeding traditional species delimitation. Therefore, recent studies concentrated on molecular genetic methods to solve taxonomic issues. Different delimitation methods were employed to explore species boundaries and the presence of cryptic species. However, the performance of the different delimitation methods needs to be tested. A particularly promising nemertean genus in this regard is the palaeonemertean genus Cephalothrix that is commonly found in European waters. In order to gain information on the number and distribution of European cephalotrichids and to test different tree‐based and non‐tree‐based delimitation methods, we analyzed a dataset comprising the barcoding region of the mitochondrial cytochrome c oxidase subunit I (COI) of 215 European Cephalothrix specimens, of which 78 were collected for this study. Our results show the presence of 12–13 European lineages of which several can be assigned to known European species. Analyzing a second dataset comprising 74 additional sequences from the Pacific and the Atlantic Oceans helped identify some of the unassigned European specimens. One resulting clade seems to represent a non‐native introduced Cephalothrix species, while another has never been recorded from Europe before. In our analysis, especially the tree‐based methods and the phylogenetic analysis proved to be a useful tool when delimiting species. It remains unclear whether the different identified clades result from cryptic speciation or from a high genetic variability of the COI gene.     

Delimiting the rare,endangered and actively speciating

Species delimitation is a contentious topic. The genomics revolution initially brought hope that identifying and classifying species would be easier through better methods and more data, but genomics has also brought complexity and controversy to delimitation. One solution can be to collect a larger sample of individuals at a finer geographic scale. But what if taxa are rare and collecting more samples is difficult or detrimental to the organisms at hand? In this issue of Molecular Ecology Resources, Opatova et al. (2023) tackle the ambiguity of species delimitation in rare and endangered trapdoor spiders (genus Cyclocosmia). The authors propose a framework for delimiting species when samples are hard to come by, such as in these rare and cryptic spiders. The authors combine extensive genomic sampling with statistical approaches that consider both the genetic distinctiveness of each population of spiders and how much gene flow occurs between these populations. Their proposed taxonomy balances two opposing signals, structure and gene flow, to count eight lineages of Cyclocosmia, and to point the way for future taxonomic studies of the rare or difficult to obtain.     

Species detection and individual assignment in species delimitation: can integrative data increase efficacy?

Statistical species delimitation usually relies on singular data, primarily genetic, for detecting putative species and individual assignment to putative species. Given the variety of speciation mechanisms, singular data may not adequately represent the genetic, morphological and ecological diversity relevant to species delimitation. We describe a methodological framework combining multivariate and clustering techniques that uses genetic, morphological and ecological data to detect and assign individuals to putative species. Our approach recovers a similar number of species recognized using traditional, qualitative taxonomic approaches that are not detected when using purely genetic methods. Furthermore, our approach detects groupings that traditional, qualitative taxonomic approaches do not. This empirical test suggests that our approach to detecting and assigning individuals to putative species could be useful in species delimitation despite varying levels of differentiation across genetic, phenotypic and ecological axes. This work highlights a critical, and often overlooked, aspect of the process of statistical species delimitation—species detection and individual assignment. Irrespective of the species delimitation approach used, all downstream processing relies on how individuals are initially assigned, and the practices and statistical issues surrounding individual assignment warrant careful consideration.     

Detecting cryptic species in sympatry and allopatry: analysis of hidden diversity in Polyommatus (Agrodiaetus) butterflies (Lepidoptera: Lycaenidae)

Modern multilocus molecular techniques are a powerful tool in the detection and analysis of cryptic taxa. However, its shortcoming is that with allopatric populations it reveals phylogenetic lineages, not biological species. The increasing power of coalescent multilocus analysis leads to the situation in which nearly every geographically isolated or semi‐isolated population can be identified as a lineage and therefore raised to species rank. It leads to artificial taxonomic inflation and as a consequence creates an unnecessary burden on the conservation of biodiversity. To solve this problem, we suggest combining modern lineage delimitation techniques with the biological species concept. We discuss several explicit principles on how genetic markers can be used to detect cryptic entities that have properties of biological species (i.e. of actually or potentially reproductively isolated taxa). Using these principles we rearranged the taxonomy of the butterfly species close to Polyommatus (Agrodiaetus) ripartii. The subgenus Agrodiaetus is a model system in evolutionary research, but its taxonomy is poorly elaborated because, as a rule, most of its species are morphologically poorly differentiated. The taxon P. (A.) valiabadi has been supposed to be one of the few exceptions from this rule due to its accurately distinguishable wing pattern. We discovered that in fact traditionally recognized P. valiabadi is a triplet of cryptic species, strongly differentiated by their karyotypes and mitochondrial haplotypes.     

Molecules,ecology, morphology,and songs in concert: how many species is Arremon torquatus (Aves: Emberizidae)?

The acceptance of the generalized or unified concept of species (i.e. that species are segments of population lineages) implies that an important task for systematists is to focus on identifying lineages and on testing hypotheses about the acquisition of properties such as phenotypic diagnosability, reciprocal monophyly, or mechanisms of reproductive isolation. However, delimiting species objectively remains one of the most challenging problems faced by biologists. In the present study, we begin to tackle the thorny issue of species delimitation in a complicated group of Neotropical passerine birds (the Arremon torquatus complex, Emberizidae) in which sets of characters vary substantially across space, but do not obviously vary in a concerted fashion. To earlier discussions of species limits in the group, we add a historical perspective offered by a recent molecular phylogeny, present quantitative analyses of morphological and vocal variation, and incorporate ecological niche models as a new tool that aids species delimitation by highlighting cases of ecological distinctiveness and cases where populations appear to be in independent evolutionary trajectories, despite being connected by environments unlikely to represent barriers to gene flow. We demonstrate that at least one pair of taxa (and likely another) currently treated as conspecific are, in fact, distinct lineages that merit species status under essentially any species criterion. However, other pairwise comparisons are not as straightforward owing to nonconcordant patterns of variation in different traits and to the impossibility of distinguishing which characters are causes and which are consequences of reproductive (and evolutionary) isolation. After considering several alternatives, we propose a provisional classification of the complex recognizing eight tentative species‐level taxa. Although this classification is likely to change as more detailed work is conducted, it provides a better foundation for studying the biology of these birds and helps to better describe their diversity, which is obscured when all taxa are subsumed into a single species name. The present study highlights several outstanding challenges, both practical and conceptual, for future studies. © 2010 The Linnean Society of London, Biological Journal of the Linnean Society, 2010, 99 , 152–176.     

“Forms” of water mites (Acari: Hydrachnidia): intraspecific variation or valid species?

In many groups of organisms, especially in the older literature, it has been common practice to recognize sympatrically occurring phenotypic variants of a species as “forms”. However, what these forms really represent often remains unclear, especially in poorly studied groups. With new algorithms for DNA‐based species delimitation, the status of forms can be explicitly tested with molecular data. In this study, we test a number of what is now recognized as valid species of water mites (Hydrachnidia), but have in the past been treated as forms sympatrically occurring with their nominate species. We also test a form without prior taxonomical status, using DNA and morphometrics. The barcoding fragment of COI, nuclear 28S and quantitative analyses of morphological data were used to test whether these taxa merit species status, as suggested by several taxonomists. Our results confirm valid species. Genetic distances between the form and nominate species (Piona dispersa and Piona variabilis, COI 11%), as well as likelihood ratio tests under the general mixed‐Yule coalescent model, supported that these are separately evolving lineages as defined by the unified species concept. In addition, they can be diagnosed with morphological characters. The study also reveals that some taxa genetically represent more than one species. We propose that P. dispersa are recognized as valid taxa at the species level. Unionicola minor (which may consist of several species), Piona stjordalensis, P. imminuta s. lat., and P. rotundoides are confirmed as species using this model. The results also imply that future studies of other water mite species complexes are likely to reveal many more genetically and morphologically distinct species.     

Tracking climate change in a dispersal-limited species: reduced spatial and genetic connectivity in a montane salamander

Tropical montane taxa are often locally adapted to very specific climatic conditions, contributing to their lower dispersal potential across complex landscapes. Climate and landscape features in montane regions affect population genetic structure in predictable ways, yet few empirical studies quantify the effects of both factors in shaping genetic structure of montane-adapted taxa. Here, we considered temporal and spatial variability in climate to explain contemporary genetic differentiation between populations of the montane salamander, Pseudoeurycea leprosa. Specifically, we used ecological niche modelling (ENM) and measured spatial connectivity and gene flow (using both mtDNA and microsatellite markers) across extant populations of P. leprosa in the Trans-Mexican Volcanic Belt (TVB). Our results indicate significant spatial and genetic isolation among populations, but we cannot distinguish between isolation by distance over time or current landscape barriers as mechanisms shaping population genetic divergences. Combining ecological niche modelling, spatial connectivity analyses, and historical and contemporary genetic signatures from different classes of genetic markers allows for inference of historical evolutionary processes and predictions of the impacts future climate change will have on the genetic diversity of montane taxa with low dispersal rates. Pseudoeurycea leprosa is one montane species among many endemic to this region and thus is a case study for the continued persistence of spatially and genetically isolated populations in the highly biodiverse TVB of central Mexico.     

Phylogeny and species delimitation of the C‐genome diploid species in Oryza

Abstract The diploid Oryza species with C‐genome type possesses abundant genes useful for rice improvement and provides parental donors of many tetraploid species with the C‐genome (BBCC, CCDD). Despite extensive studies, the phylogenetic relationship among the C‐genome species and the taxonomic status of some taxa remain controversial. In this study, we reconstructed the phylogeny of three diploid species with C‐genome (Oryza officinalis, O. rhizomatis, and O. eichingeri) based on sequences of 68 nuclear single‐copy genes. We obtained a fully resolved phylogenetic tree, clearly indicating the sister relationship of O. officinalis and O. rhizomatis, with O. eichingeri being the more divergent lineage. Incongruent phylogenies of the C‐genome species found in previous studies might result from lineage sorting, introgression/hybridization and limited number of genetic markers used. We further applied a recently developed Bayesian species delimitation method to investigate the species status of the Sri Lankan and African O. eichingeri. Analyses of two datasets (68 genes with a single sample, and 10 genes with multiple samples) support the distinct species status of the Sri Lankan and African O. eichingeri. In addition, we evaluated the impact of the number of sampled individuals and loci on species delimitation. Our simulation suggests that sampling multiple individuals is critically important for species delimitation, particularly for closely related species.