Delimiting species of marine gastropods (Turridae, Conoidea) using RAD sequencing in an integrative taxonomy framework

Coupling DNA barcodes and exon-capture to resolve the phylogeny of Turridae (Gastropoda, Conoidea)

Taxon sampling in most phylogenomic studies is often based on known taxa and/or morphospecies, thus ignoring undescribed diversity and/or cryptic lineages. The family Turridae is a group of venomous snails within the hyperdiverse superfamily Conoidea that includes many undescribed and cryptic species. Therefore 'traditional' taxon sampling could constitute a strong risk of undersampling or oversampling Turridae lineages. To minimize potential biases, we establish a robust sampling strategy, from species delimitation to phylogenomics. More than 3,000 cox-1 "barcode" sequences were used to propose 201 primary species hypotheses, nearly half of them corresponding to species potentially new to science, including several cryptic species. A 110-taxa exon-capture tree, including species representatives of the diversity uncovered with the cox-1 dataset, was build using up to 4,178 loci. Our results show the polyphyly of the genus Gemmula, that is split into up to 10 separate lineages, of which half would not have been detected if the sampling strategy was based only on described species. Our results strongly suggest that the use of blind, exploratory and intensive barcode sampling is necessary to avoid sampling biases in phylogenomic studies.

Species delimitation using genomic data to resolve taxonomic uncertainties in a speciation continuum of pelagic seabirds

Speciation is a continuous and complex process shaped by the interaction of numerous evolutionary forces. Despite the continuous nature of the speciation process, the implementation of conservation policies relies on the delimitation of species and evolutionary significant units (ESUs). Puffinus shearwaters are globally distributed and threatened pelagic seabirds. Due to remarkable morphological status the group has been under intense taxonomic debate for the past three decades. Here, we use double digest Restriction-Site Associated DNA sequencing (ddRAD-Seq) to genotype species and subspecies of North Atlantic and Mediterranean Puffinus shearwaters across their entire geographical range. We assess the phylogenetic relationships and population structure among and within the group, evaluate species boundaries, and characterise the genomic landscape of divergence. We find that current taxonomies are not supported by genomic data and propose a more accurate taxonomy by integrating genomic information with other sources of evidence. Our results show that several taxon pairs are at different stages of a speciation continuum. Our study emphasises the potential of genomic data to resolve taxonomic uncertainties, which can help to focus management actions on relevant taxa, even if they do not necessarily coincide with the taxonomic rank of species.

Genomics and ecological modelling clarify species integrity in a confusing group of butterflies

Recent advances in both genomics and ecological modelling present new, multidisciplinary opportunities for resolving species boundaries and understanding the mechanisms that maintain their integrity in regions of contact. Here, we use a combination of high-throughput DNA sequencing and ecological niche modelling to resolve species boundaries and niche divergence within the Speyeria atlantis-hesperis (Lepidoptera: Nymphalidae) complex, a confusing group of North American butterflies. This complex is notorious for its muddled species delimitations, morphological ambiguity, and extensive mito-nuclear discordance. Our admixture and multispecies coalescent-based analyses of single nucleotide polymorphisms identified substantial divergences between S. atlantis and S. hesperis in areas of contact, as well as between distinct northern and southern lineages within S. hesperis. Our results also provide evidence of past introgression relating to another species, S. zerene, which previous work has shown to be more distantly related to the S. atlantis-hesperis complex. We then used ecological models to predict habitat suitability for each of the three recovered genomic lineages in the S. atlantis-hesperis complex and assess their pairwise niche divergence. These analyses resolved that these three lineages are significantly diverged in their respective niches and are not separated by discontinuities in suitable habitat that might present barriers to gene flow. We therefore infer that ecologically-mediated selection resulting in disparate habitat associations is a principal mechanism reinforcing their genomic integrity. Overall, our results unambiguously support significant evolutionary and ecological divergence between the northern and southern lineages of S. hesperis, sufficient to recognize the southern evolutionary lineage as a distinct species, called S. nausicaa based on taxonomic priority.

Evaluating Species Delimitation Methods in Radiations: The Land Snail Albinaria cretensis Complex on Crete

Delimiting species in radiations is notoriously difficult because of the small differences between the incipient species, the star-like tree with short branches between species, incomplete lineage sorting, and the possibility of introgression between several of the incipient species. Next generation sequencing data may help to overcome some of these problems. We evaluated methods for species delimitation based on genome-wide markers in a land snail radiation on Crete. Species delimitation in the Albinaria cretensis group was based exclusively on shell characters until now and resulted in classifications distinguishing 3-9 species. We generated sequences of 4270 loci for 140 specimens of the Albinaria cretensis group from 48 populations by double-digest restriction site-associated DNA sequencing. We evaluated three methods for species discovery. The multispecies coalescent approach implemented in the program Bayesian Phylogenetics and Phylogeography resulted in a drastic overestimating of the number of species, whereas Gaussian clustering resulted in an overlumping. Primary species hypotheses based on the maximum percentage of the genome of the individuals derived from ancestral populations as estimated with the program ADMIXTURE moderately overestimated the number of species, but this was the only approach that provided information about gene flow between groups. Two of the methods for species validation that we applied, BFD* and delimitR, resulted in an acceptance of almost all primary species hypotheses, even such based on arbitrary subdivisions of hypotheses based on ADMIXTURE. In contrast, secondary species hypotheses, resulting from an evaluation of primary species hypotheses based on ADMIXTURE with isolation by distance tests, approached the morphological classification, but also uncovered two cryptic species and indicated that some of the previously delimited units should be combined. Thus, we recommend this combination of approaches that provided more detailed insights in the distinctness of barriers between the taxa of a species complex and the spatial distribution of admixture between them than the other methods. The recognition and delimitation of undersampled species remained a major challenge.

Phylogenomics and species delimitation for effective conservation of manta and devil rays

Practical biodiversity conservation relies on delineation of biologically meaningful units. Manta and devil rays (Mobulidae) are threatened worldwide, yet morphological similarities and a succession of recent taxonomic changes impede the development of an effective conservation strategy. Here, we generate genome-wide single nucleotide polymorphism (SNP) data from a geographically and taxonomically representative set of manta and devil ray samples to reconstruct phylogenetic relationships and evaluate species boundaries under the general lineage concept. We show that nominal species units supported by alternative data sources constitute independently evolving lineages, and find robust evidence for a putative new species of manta ray in the Gulf of Mexico. Additionally, we uncover substantial incomplete lineage sorting indicating that rapid speciation together with standing variation in ancestral populations has driven phylogenetic uncertainty within Mobulidae. Finally, we detect cryptic diversity in geographically distinct populations, demonstrating that management below the species level may be warranted in certain species. Overall, our study provides a framework for molecular genetic species delimitation that is relevant to wide-ranging taxa of conservation concern, and highlights the potential for genomic data to support effective management, conservation and law enforcement strategies.

Lineage Identification Affects Estimates of Evolutionary Mode in Marine Snails

In order to study evolutionary pattern and process, we need to be able to accurately identify species and the evolutionary lineages from which they are derived. Determining the concordance between genetic and morphological variation of living populations, and then directly comparing extant and fossil morphological data, provides a robust approach for improving our identification of lineages through time. We investigate genetic and shell morphological variation in extant species of Penion marine snails from New Zealand, and extend this analysis into deep time using fossils. We find that genetic and morphological variation identify similar patterns and support most currently recognized extant species. However, some taxonomic over-splitting is detected due to shell size being a poor trait for species delimitation, and we identify incorrect assignment of some fossil specimens. We infer that a single evolutionary lineage (Penion sulcatus) has existed for 22 myr, with most aspects of shell shape and shell size evolving under a random walk. However, by removing samples previously classified as the extinct species P. marwicki, we instead detect morphological stasis for one axis of shell shape variation. This result demonstrates how lineage identification can change our perception of evolutionary pattern and process. [Genotyping by sequencing; geometric morphometrics; morphological evolution; Neogastropoda; phenotype; speciation; stasis.].

Where the snails have no name: a molecular phylogeny of Raphitomidae (Neogastropoda: Conoidea) uncovers vast unexplored diversity in the deep seas of temperate southern and eastern Australia

Although raphitomid snails are a dominant component of gastropod communities in deep seas worldwide, their systematics is still largely tentative. We assembled the most complete sampling of Raphitomidae from south-eastern Australia to date. Based on morphological and molecular data from this material, we produced a robust phylogenetic framework and used it to delimit genera. For the focus area, our results show a large proportion of undescribed species- and genus-level taxa, 11 of which are formally described herein. We demonstrate that the examination of purely morphological characters rarely suffices for the purpose of accurate genus delimitation. As a result, some traditionally highly diverse raphitomid genera (such as Gymnobela) turn out to be artificial assemblages of several unrelated, mostly undescribed, genus-level lineages. Our data suggest that comparable configurations of shell and radular features, observed at the genus level, commonly do not reflect true phylogenetic relationships. However, our results are inconclusive as to whether homoplasy or conservatism are the drivers of this phenomenon. Accommodating for the inevitable sampling biases, south-eastern Australia appears as a possible hotspot for both raphitomid diversity and endemism, when compared with adjacent areas.

Diversity of the Pterasteridae (Asteroidea) in the Southern Ocean: a molecular and morphological approach

An integrative approach is crucial in discrimination of species, especially for taxa that are difficult to identify based on morphological characters. In this study, we combine genetics and morphology to assess the diversity of Pterasteridae, a sea star family diversified in deep-sea and polar environments. Because of their derived anatomy and the frequent loss of characters during preservation, Pterasteridae are a suitable case for an integrative study. The molecular identification of 191 specimens (mostly from the Southern Ocean) suggests 26–33 species in three genera (Diplopteraster, Hymenaster and Pteraster), which match the morphological identification in 54–62% of cases. The mismatches are either different molecular units that are morphologically indistinguishable (e.g. Pteraster stellifer units 2 and 4) or, conversely, nominal species that are genetically identical (e.g. Hymenaster coccinatus/densus/praecoquis). Several species are shared between the Northern and Southern Hemispheres (e.g. Pteraster jordani/affinis). In conclusion, the taxonomic status of some groups is confirmed, but for others we find the need to re-evaluate the taxonomy at both genus and species levels. This work significantly increases the DNA barcode library of the Southern Ocean species and merges taxonomic information into an identification key that could become a baseline for future studies (pterasteridae-so.identificationkey.org).

Just the once will not hurt: DNA suggests species lumping over two oceans in deep-sea snails (Cryptogemma)

The practice of species delimitation using molecular data commonly leads to the revealing of species complexes and an increase in the number of delimited species. In a few instances, however, DNA-based taxonomy has led to lumping together of previously described species. Here, we delimit species in the genus Cryptogemma (Gastropoda: Conoidea: Turridae), a group of deep-sea snails with a wide geographical distribution, primarily by using the mitochondrial COI gene. Three approaches of species delimitation (ABGD, mPTP and GMYC) were applied to define species partitions. All approaches resulted in eight species. According to previous taxonomic studies and shell morphology, 23 available names potentially apply to the eight Cryptogemma species that were recognized herein. Shell morphometrics, radular characters and geographical and bathymetric distributions were used to link type specimens to these delimited species. In all, 23 of these available names are here attributed to seven species, resulting in 16 synonymizations, and one species is described as new: Cryptogemma powelli sp. nov. We discuss the possible reasons underlying the apparent overdescription of species within Cryptogemma, which is shown here to constitute a rare case of DNA-based species lumping in the hyper-diversified superfamily Conoidea.

Venom Diversity and Evolution in the Most Divergent Cone Snail Genus Profundiconus

Profundiconus is the most divergent cone snail genus and its unique phylogenetic position, sister to the rest of the family Conidae, makes it a key taxon for examining venom evolution and diversity. Venom gland and foot transcriptomes of Profundiconus cf. vaubani and Profundiconus neocaledonicus were de novo assembled, annotated, and analyzed for differential expression. One hundred and thirty-seven venom components were identified from P. cf. vaubani and 82 from P. neocaledonicus, with only four shared by both species. The majority of the transcript diversity was composed of putative peptides, including conotoxins, profunditoxins, turripeptides, insulin, and prohormone-4. However, there were also a significant percentage of other putative venom components such as chymotrypsin and L-rhamnose-binding lectin. The large majority of conotoxins appeared to be from new gene superfamilies, three of which are highly different from previously reported venom peptide toxins. Their low conotoxin diversity and the type of insulin found suggested that these species, for which no ecological information are available, have a worm or molluscan diet associated with a narrow dietary breadth. Our results indicate that Profundiconus venom is highly distinct from that of other cone snails, and therefore important for examining venom evolution in the Conidae family.

Genomic and Mitochondrial Data Identify Different Species Boundaries in Aposematically Polymorphic Eniclases Net-Winged Beetles (Coleoptera: Lycidae)

Species delineation is essential for any evolutionary and biodiversity research, and recent advances in genomic sequencing have made it possible to robustly define species boundaries and detect hidden diversity. Here, we studied 14 species of aposematically colored New Guinean Eniclases (Coleoptera: Lycidae) whose conventional morphology- and single-locus mtDNA-based taxonomy has been contentious. We analyzed mitochondrial and restriction site associated DNA fragments to obtain a phylogenetic hypothesis and compared relationships recovered by the RAD analysis with species limits based on other information. The results show the presence of cryptic diversity and common mitonuclear discordance when over 30% of individuals were incorrectly assigned to species if only mitogenomic markers were considered. Nuclear data falsified the species rank of one species and identified one earlier unrecognized lineage deserving species rank. Further, our analyses demonstrate a highly variable phenotypic differentiation, with several pairs of cryptic species standing in contrast with genetically close but phenotypically highly divergent lineages. We show that morphological and mitogenomic analyses produce reliable information for taxonomy in most cases. Nevertheless, the species boundaries among closely related species should be based on all lines of evidence, including nuclear markers.