Defining reproductively isolated units in a cryptic and syntopic species complex using mitochondrial and nuclear markers: the brooding brittle star, Amphipholis squamata (Ophiuroidea)

Genetic variation in the brooding brittle-star: a global hybrid polyploid complex?

The widespread and abundant brooding brittle-star (Amphipholis squamata) is a simultaneous hermaphrodite with a complex mitochondrial phylogeography of multiple divergent overlapping mtDNA lineages, high levels of inbreeding or clonality and unusual sperm morphology. We use exon-capture and transcriptome data to show that the nuclear genome comprises multiple (greater than 3) divergent (π > 6%) expressed components occurring across samples characterized by highly divergent (greater than 20%) mitochondrial lineages, and encompassing several other genera, including diploid dioecious species. We report a massive sperm genome size in A. squamata, an order of magnitude larger than that present in other brittle-stars, and consistent with our SNP-based measure of greatly elevated ploidy. Similarity of these genetic signatures to well-known animal systems suggests that A. squamata (and related taxa) is a hybrid polyploid asexual complex of variable subgenome origins, ploidy and reproductive mode. We discuss enigmatic aspects of A. squamata biology in this light. This putative allopolyploid complex would be the first to be reported from the phylum Echinodermata.

Pseudocryptic diversity and species boundaries in the sea cucumber Stichopus cf. horrens (Echinodermata: Stichopodidae) revealed by mitochondrial and microsatellite markers

Morphologically cryptic and pseudo-cryptic species pose a challenge to taxonomic identification and assessments of species diversity and distributions. Such is the case for the sea cucumber Stichopus horrens, commonly confused with Stichopus monotuberculatus. Here, we used mitochondrial cytochrome oxidase subunit I (COI) and microsatellite markers to examine genetic diversity in Stichopus cf. horrens throughout the Philippine archipelago, to aid species identification and clarify species boundaries. Phylogenetic analysis reveals two recently diverged COI lineages (Clade A and Clade B; c. 1.35-2.54 Mya) corresponding to sequence records for specimens identified as S. monotuberculatus and S. horrens, respectively. Microsatellite markers reveal two significantly differentiated genotype clusters broadly concordant with COI lineages (Cluster 1, Cluster 2). A small proportion of individuals were identified as later-generation hybrids indicating limited contemporary gene flow between genotype clusters, thus confirming species boundaries. Morphological differences in papillae distribution and form are observed for the two species, however tack-like spicules from the dorsal papillae are not a reliable diagnostic character. An additional putative cryptic species was detected within Clade B-Cluster 2 specimens warranting further examination. We propose that these lineages revealed by COI and genotype data be referred to as Stichopus cf. horrens species complex.

Exceptional population genomic homogeneity in the black brittle star Ophiocomina nigra (Ophiuroidea, Echinodermata) along the Atlantic-Mediterranean coast

The Atlantic-Mediterranean marine transition is characterised by strong oceanographic barriers and steep environmental gradients that generally result in connectivity breaks between populations from both basins and may lead to local adaptation. Here, we performed a population genomic study of the black brittle star, Ophiocomina nigra, covering most of its distribution range along the Atlantic-Mediterranean region. Interestingly, O. nigra is extremely variable in its coloration, with individuals ranging from black to yellow-orange, and different colour morphs inhabiting different depths and habitats. In this work, we used a fragment of the mitochondrial COI gene and 2,374 genome-wide ddRADseq-derived SNPs to explore: (a) whether the different colour morphs of O. nigra represent different evolutionary units; (b) the disruptive effects of major oceanographic fronts on its population structure; and (c) genomic signals of local adaptation to divergent environments. Our results revealed exceptional population homogeneity, barely affected by oceanographic fronts, with no signals of local adaptation nor genetic differentiation between colour morphs. This remarkable panmixia likely results from a long pelagic larval duration, a large effective population size and recent demographic expansions. Our study unveils an extraordinary phenotypic plasticity in O. nigra, opening further research questions on the ecological and molecular mechanisms underpinning coloration in Ophiuroidea.

Integrative species delimitation in the common ophiuroid Ophiothrix angulata (Echinodermata: Ophiuroidea): insights from COI, ITS2, arm coloration, and geometric morphometrics

Ophiothrix angulata (Say, 1825) is one of the most common and well-known ophiuroids in the Western Atlantic, with a wide geographic and bathymetric range. The taxonomy of this species has been controversial for a century because of its high morphological variability. Here we integrate information from DNA sequence data, color patterns, and geometric morphometrics to assess species delimitation and geographic differentiation in O. angulata. We found three deeply divergent mtDNA-COI clades (K2P 17.0-27.9%). ITS2 nuclear gene and geometric morphometrics of dorsal and ventral arm plates differentiate one of these lineages, as do integrative species delineation analyses, making this a confirmed candidate species.

Development of two new sets of PCR primers for eDNA metabarcoding of brittle stars (Echinodermata, Ophiuroidea)

Brittle stars (class Ophiuroidea) are marine invertebrates comprising approximately 2,100 extant species, and are considered to constitute the most diverse taxon of the phylum Echinodermata. As a non-invasive method for monitoring biodiversity, we developed two new sets of PCR primers for metabarcoding environmental DNA (eDNA) from brittle stars. The new primer sets were designed to amplify 2 short regions of the mitochondrial 16S rRNA gene, comprising a conserved region (111–115 bp, 112 bp on average; named “16SOph1”) and a hyper-variable region (180–195 bp, 185 bp on average; named “16SOph2”) displaying interspecific variation. The performance of the primers was tested using eDNA obtained from two sources: a) rearing water of an 2.5 or 170 L aquarium tanks containing 15 brittle star species and b) from natural seawater collected around Misaki, the Pacific coast of central Japan, at depths ranging from shallow (2 m) to deep (> 200 m) sea. To build a reference library, we obtained 16S rRNA sequences of brittle star specimens collected from around Misaki and from similar depths in Japan, and sequences registered in International Nucleotide Sequence Database Collaboration. As a result of comparison of the obtained eDNA sequences with the reference library 37 (including cryptic species) and 26 brittle star species were detected with certain identities by 16SOph1 and 16SOph2 analyses, respectively. In shallow water, the number of species and reads other than the brittle stars detected with 16SOph1 was less than 10% of the total number. On the other hand, the number of brittle star species and reads detected with 16SOph2 was less than half of the total number, and the number of detected non-brittle star metazoan species ranged from 20 to 46 species across 6 to 8 phyla (only the reads at the “Tank” were less than 0.001%). The number of non-brittle star species and reads at 80 m was less than 10% with both of the primer sets. These findings suggest that 16SOph1 is specific to the brittle star and 16SOph2 is suitable for a variety of marine metazoans. It appears, however, that further optimization of primer sequences would still be necessary to avoid possible PCR dropouts from eDNA extracts. Moreover, a detailed elucidation of the brittle star fauna in the examined area, and the accurate identification of brittle star species in the current DNA databank is required.

Mitochondrial phylogeny of the brittle star genus Ophioderma

We reconstructed the mitochondrial phylogeny of the species of the brittle star genus Ophioderma, using sequences of the Cytochrome Oxidase I gene (COI) to address four questions: (i) Are the species of Ophioderma described on morphological evidence reflected in mitochondrial genealogy? (ii) Which species separated from which? (iii) When did speciation events occur? (iv) What is the rate of COI evolution in ophiuroids? We found that most of the 22 described species we sampled coincide with monophyletic clusters of COI sequences, but there are exceptions, most notably in the eastern Pacific, in which three undescribed species were indicated. The COI phylogeny lacks resolution in the deeper nodes, but it does show that there are four species pairs, the members of which are found on either side of the central American Isthmus. Two pairs with a genetic distance of ~ 4% between Atlantic and Pacific members were probably split during the final stages of Isthmus completion roughly 3 million years ago. The rate of divergence provided by these pairs allowed the calibration of a relaxed molecular clock. Estimated dates of divergence indicate that the lineages leading to extant species coalesce at times much older than congeneric species in other classes of echinoderms, suggesting that low extinction rates may be one of the reasons that ophiuroids are species-rich. The mean rate of COI substitution in Ophioderma is three times slower than that of echinoids. Conclusions of previous mitochondrial DNA studies of ophiuroids that relied on echinoid calibrations to determine divergence times need to be revised.

Temporal and spatial variation in population structure among brooding sea stars in the genus Leptasterias

Temporal genetic studies of low-dispersing organisms are rare. Marine invertebrates lacking a planktonic larval stage are expected to have lower dispersal, low gene flow, and a higher potential for local adaptation than organisms with planktonic dispersal. Leptasterias is a genus of brooding sea stars containing several cryptic species complexes. Population genetic methods were used to resolve patterns of fine-scale population structure in central California Leptasterias species using three loci from nuclear and mitochondrial genomes. Historic samples (collected between 1897 and 1998) were compared to contemporary samples (collected between 2008 and 2014) to delineate changes in species distributions in space and time. Phylogenetic analysis of contemporary samples confirmed the presence of a bay-localized clade and revealed the presence of an additional bay-localized and previously undescribed clade of Leptasterias. Analysis of contemporary and historic samples indicates two clades are experiencing a constriction in their southern range limit and suggests a decrease in clade-specific abundance at sites at which they were once prevalent. Historic sampling revealed a dramatically different distribution of diversity along the California coastline compared to contemporary sampling and illustrates the importance of temporal genetic sampling in phylogeographic studies. These samples were collected prior to significant impacts of Sea Star Wasting Disease (SSWD) and represent an in-depth analysis of genetic structure over 117 years prior to the SSWD-associated mass die-off of Leptasterias.

Speciation across the Tree of Life

Much of what we know about speciation comes from detailed studies of well-known model systems. Although there have been several important syntheses on speciation, few (if any) have explicitly compared speciation among major groups across the Tree of Life. Here, we synthesize and compare what is known about key aspects of speciation across taxa, including bacteria, protists, fungi, plants, and major animal groups. We focus on three main questions. Is allopatric speciation predominant across groups? How common is ecological divergence of sister species (a requirement for ecological speciation), and on what niche axes do species diverge in each group? What are the reproductive isolating barriers in each group? Our review suggests the following patterns. (i) Based on our survey and projected species numbers, the most frequent speciation process across the Tree of Life may be co-speciation between endosymbiotic bacteria and their insect hosts. (ii) Allopatric speciation appears to be present in all major groups, and may be the most common mode in both animals and plants, based on non-overlapping ranges of sister species. (iii) Full sympatry of sister species is also widespread, and may be more common in fungi than allopatry. (iv) Full sympatry of sister species is more common in some marine animals than in terrestrial and freshwater ones. (v) Ecological divergence of sister species is widespread in all groups, including ~70% of surveyed species pairs of plants and insects. (vi) Major axes of ecological divergence involve species interactions (e.g. host-switching) and habitat divergence. (vii) Prezygotic isolation appears to be generally more widespread and important than postzygotic isolation. (viii) Rates of diversification (and presumably speciation) are strikingly different across groups, with the fastest rates in plants, and successively slower rates in animals, fungi, and protists, with the slowest rates in prokaryotes. Overall, our study represents an initial step towards understanding general patterns in speciation across all organisms.

Mitogenome phylogenetics in the genus Palaemon (Crustacea: Decapoda) sheds light on species crypticism in the rockpool shrimp P. elegans

The genus Palaemon comprises worldwide marine and freshwater shrimps and prawns, and some of them are ecologically or commercially important species. Palaemon is not currently a monophyletic group, so phylogenetics and systematics are constantly changing. Species crypticism has been pointed out in several Palaemon species, being the clearest evidence in the European rockpool shrimp P. elegans. Here we sequenced and described seven European Palaemon mitochondrial genomes. The mitochondrial protein-coding genes were used, along with those of three other Palaemon species, to perform mitogenome phylogenetic analyses to clarify the evolutionary relationships within the genus, and particularly to shed light on the cryptic species found within P. elegans. The Messinian Salinity Crisis (5.3-5.9 Ma, late Miocene) was proposed to be the origin of this cryptic species and it was used as aged constraint for calibration analysis. We provide the largest and the first time-calibrated mitogenome phylogeny of the genus Palaemon and mitogenome substitution rate was estimated (1.59% per million years) in Decapoda for the first time. Our results highlighted the need for future systematics changes in Palaemon and crypticism in P. elegans was confirmed. Mitochondrial genome and cox1 (1.41%) substitution rate estimates matched those published elsewhere, arguing that the Messinian Salinity Crisis was a plausible event driving the split between P. elegans and its cryptic species. Molecular dating suggested that Pleistocene glaciations were likely involved in the differentiation between the Atlantic and Mediterranean populations of P. elegans. On the contrary, the divergence between the Atlantic and Mediterranean populations of the common littoral shrimp P. serratus was greater and dated to be much older (4.5-12.3 Ma, Plio-Miocene), so we considered that they could represent two separated species. Therefore, species crypticism in the genus Palaemon seems to be a common phenomenon.

Genetic diversity and population structure of the rockpool shrimp Palaemon elegans based on microsatellites: evidence for a cryptic species and differentiation across the Atlantic-Mediterranean transition

The rockpool shrimp Palaemon elegans is an ecologically important crustacean species within the European coastline fauna. In the present study, genetic diversity and population structure and connectivity were assessed by examining 21 polymorphic microsatellite loci at 13 sampling sites located along the natural distribution range. All localities showed similar levels of genetic variability. Significant deficits of heterozygosity were recorded, most likely due to the presence of null alleles. Genetic structure analyses revealed two clearly genetically distinct groups within P. elegans but without following any geographical or oceanographic basis. Thus, our results provided nuclear evidence for the existence of a Mediterranean cryptic species within P. elegans, highlighting the need to revise its taxonomic status. Regarding P. elegans sensu stricto, population structuring was reported across the Atlantic-Mediterranean transition area, where the Almería-Orán Front restricts the gene flow between the Atlantic and the Mediterranean population. Moreover, while population connectivity was suggested between all Mediterranean localities, some substructure was found within the Atlantic group. Canary Islands exhibited a weak but significant genetic differentiation from all Atlantic mainland localities, consistent with the isolation-by-distance pattern detected throughout the Atlantic population. Overall, all these findings provided new insights into the population biology of P. elegans complex.

The Echinoderm Fauna of the Azores (NE Atlantic Ocean)

In more than 150 years of research in the waters surrounding the Azores, several publications on the fauna of echinoderms of the archipelago have been produced, in the form of papers, notes, reports, reviews, and monographs. This work attempts to summarize the present knowledge on this marine group in the Azorean exclusive economic zone (i.e., waters within 200 nautical miles of the archipelago's shores). A short review of the history of the species' taxonomy is given, with key references, geographical distribution, ecology, additional notes and, when possible, figures. We herein report 172 species of echinoderms (6 crinoids, 55 ophiuroids, 45 asteroids, 36 holothurians, and 30 echinoids) from the Azores Archipelago, most of them inhabiting deep waters (>200 m). Only 29 shallow-water species were recorded locally (≤50 m depth). In general, the echinoderm species present in the Azores are characterized by a wide geographical distribution in the Atlantic Ocean. Only nine taxa (all deep-water species, >840 m) appear to be restricted to the Azorean waters. Overall, the knowledge of the echinoderm fauna of the Azores is out-dated, with many species last collected in the archipelago over 100 years ago. A recent interest in the Azorean Mid-Atlantic waters has brought oceanographic cruises back to the archipelago, thus providing new opportunities for the renewal of 150 years of echinoderm studies in the area.

Problems and Questions Posed by Cryptic Species. A Framework to Guide Future Studies

Species are the currency of biology and important units of biodiversity, thus errors in species delimitations potentially have important consequences. During the last decades, owing to the use of genetic markers, many nominal species appeared to consist of several reproductively isolated entities called cryptic species (hereafter CS). In this chapter we explain why CS are important for practical reasons related to community and ecosystem monitoring, and for biological knowledge, particularly for understanding ecological and evolutionary processes. To find solutions to practical problems and to correct biological errors, a thorough analysis of the distinct types of CS reported in the literature is necessary and some general rules have to be identified. Here we explain how to identify CS, and we propose a rational and practical classification of CS (and putative CS), based on the crossing of distinct levels of genetic isolation with distinct levels of morphological differentiation. We also explain how to identify likely explanations for a given CS (either inherent to taxonomic processes or related to taxon biology, ecology and geography) and how to build a comprehensive database aimed at answering these practical and theoretical questions. Our pilot review of the literature in marine animals established that half of the reported cases are not CS sensu stricto (i.e. where morphology cannot distinguish the entities) and just need taxonomic revision. It also revealed significant associations between CS features, such as a higher proportion of diagnostic morphological differences in sympatric than in allopatric CS and more frequent ecological differentiation between sympatric than allopatric CS, both observations supporting the competitive exclusion theory, thus suggesting that ignoring CS causes not only species diversity but also functional diversity underestimation.

A multispecies approach reveals hot spots and cold spots of diversity and connectivity in invertebrate species with contrasting dispersal modes

Genetic diversity is crucial for species' maintenance and persistence, yet is often overlooked in conservation studies. Species diversity is more often reported due to practical constraints, but it is unknown if these measures of diversity are correlated. In marine invertebrates, adults are often sessile or sedentary and populations exchange genes via dispersal of gametes and larvae. Species with a larval period are expected to have more connected populations than those without larval dispersal. We assessed the relationship between measures of species and genetic diversity, and between dispersal ability and connectivity. We compiled data on genetic patterns and life history traits in nine species across five phyla. Sampling sites spanned 600 km in the northwest Mediterranean Sea and focused on a 50-km area near Marseilles, France. Comparative population genetic approaches yielded three main results. (i) Species without larvae showed higher levels of genetic structure than species with free-living larvae, but the role of larval type (lecithotrophic or planktotrophic) was negligible. (ii) A narrow area around Marseilles, subject to offshore advection, limited genetic connectivity in most species. (iii) We identified sites with significant positive contributions to overall genetic diversity across all species, corresponding with areas near low human population densities. In contrast, high levels of human activity corresponded with a negative contribution to overall genetic diversity. Genetic diversity within species was positively and significantly linearly related to local species diversity. Our study suggests that local contribution to overall genetic diversity should be taken into account for future conservation strategies.

DNA barcoding of reef brittle stars (Ophiuroidea, Echinodermata) from the southwestern Indian Ocean evolutionary hot spot of biodiversity

In anticipation of the current biodiversity crisis, it has become critical to rapidly and accurately assess biodiversity. DNA barcoding has proved efficient in facilitating the discovery and description of thousands of species and also provides insight into the dynamics of biodiversity. Here, we sequenced a portion of the mitochondrial cytochrome c oxidase subunit I (COI) gene from all morphospecies of reef brittle stars collected during a large‐scale biodiversity survey in the southwestern Indian Ocean (SWIO). Three methods of species delineation (Automatic Barcode Gap Discovery, Generalized Mixed Yule Coalescent model, and Bayesian Poisson Tree Processes) showed concordant results and revealed 51 shallow reef species in the region. Mean intraspecific genetic distances (0.005–0.064) and mean interspecific genetic distances within genera (0.056–0.316) were concordant with previous echinoderm studies. This study revealed that brittle‐star biodiversity is underestimated by 20% within SWIO and by >40% when including specimens from the Pacific Ocean. Results are discussed in terms of endemism, diversification processes, and conservation implications for the Indo‐West Pacific marine biodiversity. We emphasize the need to further our knowledge on biodiversity of invertebrate groups in peripheral areas.

One species for one island? Unexpected diversity and weak connectivity in a widely distributed tropical hydrozoan

Isolation by distance (IBD) is one of the main modes of differentiation in marine species, above all in species presenting low dispersal capacities. This study reports the genetic structuring in the tropical hydrozoan Lytocarpia brevirostris α (sensu Postaire et al, 2016b), a brooding species, from 13 populations in the Western Indian Ocean (WIO) and one from New Caledonia (Tropical Southwestern Pacific). At the local scale, populations rely on asexual propagation at short distance, which was not found at larger scales; identical genotypes were restricted to single populations. After the removal of repeated genotypes, all populations presented significant positive F_IS values (between 0.094*** and 0.335***). Gene flow was extremely low at all spatial scales, between sites within islands (<10 km distance) and among islands (100 to>11 000 km distance), with significant pairwise F_ST values (between 0.012*** and 0.560***). A general pattern of IBD was found at the Indo-Pacific scale, but also within sampled ecoregions of the WIO province. Clustering analyses identified each sampled island as an independent population, whereas analysis of molecular variance indicated that population genetic differentiation was significant at small (within island) and intermediate (among islands within province) spatial scales. The high population differentiation might reflect the life cycle of this brooding hydrozoan, possibly preventing regular dispersal at distances more than a few kilometres and probably leading to high cryptic diversity, each island housing an independent evolutionary lineage.

Positive selection on sperm ion channels in a brooding brittle star: consequence of life-history traits evolution

Closely related species are key models to investigate mechanisms leading to reproductive isolation and early stages of diversification, also at the genomic level. The brittle star cryptic species complex Ophioderma longicauda encompasses the sympatric broadcast-spawning species C3 and the internal brooding species C5. Here, we used de novo transcriptome sequencing and assembly in two closely related species displaying contrasting reproductive modes to compare their genetic diversity and to investigate the role of natural selection in reproductive isolation. We reconstructed 20 146 and 22 123 genes for C3 and C5, respectively, and characterized a set of 12 229 orthologs. Genetic diversity was 1.5-2 times higher in C3 compared to C5, confirming that species with low parental investment display higher levels of genetic diversity. Forty-eight genes were the targets of positive diversifying selection during the evolution of the two species. Notably, two genes (NHE and TetraKCNG) are sperm-specific ion channels involved in sperm motility. Ancestral sequence reconstructions show that natural selection targeted the two genes in the brooding species. This may result from an adaptation to the novel environmental conditions surrounding sperm in the brooding species, either directly affecting sperm or via an increase in male/female conflict. This phenomenon could have promoted prezygotic reproductive isolation between C3 and C5. Finally, the sperm receptors to egg chemoattractants differed between C3 and C5 in the ligand-binding region. We propose that mechanisms of species-specific gamete recognition in brittle stars occur during sperm chemotaxis (sperm attraction towards the eggs), contrary to other marine invertebrates where prezygotic barriers to interspecific hybridization typically occur before sperm-egg fusion.

Contrasted phylogeographic patterns on mitochondrial DNA of shallow and deep brittle stars across the Atlantic-Mediterranean area

Previous studies on Ophiothrix in European waters demonstrated the existence of two distinct species, Ophiothrix fragilis and Ophiothrix sp. II. Using phylogenetic and species delimitation techniques based on two mitochondrial genes (cytochrome c oxidase I and 16S rRNA) we prove the existence of a new congeneric species (Ophiothrix sp. III), occurring in the deep Atlantic coast of the Iberian Peninsula and the Alboran Sea. We compared phylogeographic patterns of these three Ophiothrix species to test whether closely related species are differentially affected by past demographic events and current oceanographic barriers. We used 432 sequences (137 of O. fragilis, 215 of Ophiothrix sp. II, and 80 of Ophiothrix sp. III) of the 16S rRNA from 23 Atlantic-Mediterranean locations for the analyses. We observed different geographic and bathymetric distributions, and contrasted phylogeography among species. Ophiothrix fragilis appeared genetically isolated between the Atlantic and Mediterranean basins, attributed to past vicariance during Pleistocene glaciations and a secondary contact associated to demographic expansion. This contrasts with the panmixia observed in Ophiothrix sp. II across the Atlantic-Mediterranean area. Results were not conclusive for Ophiothrix sp. III due to the lack of a more complete sampling within the Mediterranean Sea.

The effect of copper and temperature on juveniles of the eurybathic brittle star Amphipholis squamata--exploring responses related to motility and the water vascular system

The limited availability of test organisms that represent tropical and deeper water environments is a significant concern when assessing the risk of contaminants in these environments. Amphipholis squamata (Delle Chiaje 1828) is a widely distributed brittle star with many phylogenetic clades reported from different latitudes, and it also occurs from the intertidal zone to a depth of ∼1300 m. In the present study, the effect of copper on four behavioural responses and mortality of A. squamata were quantified at four different temperatures including 25, 20, 15 and 10°C. At 25°C the four behavioural responses and mortality were relatively sensitive to copper, with 96 h EC50 values of 25 (confidence interval 18-44), 24 (7-26), 32 (24-41), 29 (9-41) μg L(-1) for the measured ability to turn from the oral surface up to oral surface down, curling behaviour, tube foot movement, and tube foot retraction respectively. The average 96-h LC50 value for copper at 25°C was 46 μg L(-1). Some endpoints investigated showed significant effects of reduced temperature compared to the optimal temperature. These effects were enhanced with increasing copper concentrations and significant differences in copper toxicity between temperature treatments were most notable when measuring the ability to turn from the oral surface up to oral surface down where the EC50 changed from 25 (18 to 44) to 6 (-18 to 14) μg L(-1) with a reduction of temperature from 25 to 15°C. The results showed that A. squamata is relatively sensitive to copper and that further investigation into the effects of other stressors on these endpoints is warranted.

Genetic data, reproduction season and reproductive strategy data support the existence of biological species in Ophioderma longicauda

Cryptic species are numerous in the marine environment. The brittle star Ophioderma longicauda is composed of six mitochondrial lineages, encompassing brooders, which form a monophyletic group, and broadcasters, from which the brooders are derived. To clarify the species limits within O. longicauda, we compared the reproductive status of the sympatric lineages L1 and L3 (defined after sequencing a portion of the mitochondrial gene COI) during the month of May in Greece. In addition, we genotyped a nuclear marker, intron i51. Each L3 female was brooding, whereas all L1 specimens displayed full gonads, suggesting temporal pre-zygotic isolation between brooders and broadcasters. Statistical differences were found among lineages in morphology and bathymetric distribution. Finally, the intron i51 was polymorphic in L1 (60 individuals), but monomorphic in L3 (109 individuals), confirming the absence of gene flow between brooders and broadcasters. In conclusion, the broadcasting lineage L1 and the brooding lineage L3 are different biological species.

Hidden diversity in sardines: genetic and morphological evidence for cryptic species in the goldstripe sardinella, Sardinella gibbosa (Bleeker, 1849)

Cryptic species continue to be uncovered in many fish taxa, posing challenges for fisheries conservation and management. In Sardinella gibbosa, previous investigations revealed subtle intra-species variations, resulting in numerous synonyms and a controversial taxonomy for this sardine. Here, we tested for cryptic diversity within S. gibbosa using genetic data from two mitochondrial and one nuclear gene regions of 248 individuals of S. gibbosa, collected from eight locations across the Philippine archipelago. Deep genetic divergence and subsequent clustering was consistent across both mitochondrial and nuclear markers. Clade distribution is geographically limited: Clade 1 is widely distributed in the central Philippines, while Clade 2 is limited to the northernmost sampling site. In addition, morphometric analyses revealed a unique head shape that characterized each genetic clade. Hence, both genetic and morphological evidence strongly suggests a hidden diversity within this common and commercially-important sardine.

Cryptic speciation or global spread? The case of a cosmopolitan marine invertebrate with limited dispersal capabilities

The existence of globally-distributed species with low dispersal capabilities is a paradox that has been explained as a result of human-mediated transport and by hidden diversity in the form of unrecognized cryptic species. Both factors are not mutually exclusive, but relatively few studies have demonstrated the presence of both. Here we analyse the genetic patterns of the colonial ascidian Diplosoma listerianum, a species nowadays distributed globally. The study of a fragment of a mitochondrial gene in localities worldwide revealed the existence of multiple cryptic species. In addition, we found a complex geographic structure and multiple clades occurred in sympatry. One of the species showed strong population structure irrespective of geographical distances, which is coherent with stochastic dispersal linked to human transport. The present study shows the complexity of discerning the role of cryptic diversity from human-driven range shifts worldwide, as well as disentangling the effects of natural and artificial dispersal.

Design of phylum-specific hybrid primers for DNA barcoding: addressing the need for efficient COI amplification in the Echinodermata

Recent research has shown the usefulness of the Folmer region of the cytochrome oxidase I (COI) as a genetic barcode to assist in species delimitation of echinoderms. However, amplification of COI is often challenging in echinoderms (low success or pseudogenes). We present a method that allows the design of phylum-specific hybrid primers, and use this to develop COI primers for the Echinodermata. We aligned COI sequences from 310 echinoderm species and designed all possible primers along the consensus sequence with two methods (standard degenerate and hybrid). We found much lower degeneracy for hybrid primers (4-fold degeneracy) than for standard degenerate primers (≥48-fold degeneracy). We then designed the most conserved hybrid primers to amplify a >500-bp region within COI. These primers successfully amplified this gene region in all tested taxa (123 species across all echinoderm classes). Sequencing of 30 species among these confirmed both the quality of the sequences (>500 bp, no pseudogenes) and their utility as a DNA barcode. This method should be useful for developing primers for other mitochondrial genes and other phyla. The method will also be of interest for the development of future projects involving both community-based genetic assessments on macroorganisms and biodiversity assessment of environmental samples using high-throughput sequencing.

Large-scale species delimitation method for hyperdiverse groups

Accelerating the description of biodiversity is a major challenge as extinction rates increase. Integrative taxonomy combining molecular, morphological, ecological and geographical data is seen as the best route to reliably identify species. Classic molluscan taxonomic methodology proposes primary species hypotheses (PSHs) based on shell morphology. However, in hyperdiverse groups, such as the molluscan family Turridae, where most of the species remain unknown and for which homoplasy and plasticity of morphological characters is common, shell-based PSHs can be arduous. A four-pronged approach was employed to generate robust species hypotheses of a 1000 specimen South-West Pacific Turridae data set in which: (i) analysis of COI DNA Barcode gene is coupled with (ii) species delimitation tools GMYC (General Mixed Yule Coalescence Method) and ABGD (Automatic Barcode Gap Discovery) to propose PSHs that are then (iii) visualized using Klee diagrams and (iv) evaluated with additional evidence, such as nuclear gene rRNA 28S, morphological characters, geographical and bathymetrical distribution to determine conclusive secondary species hypotheses (SSHs). The integrative taxonomy approach applied identified 87 Turridae species, more than doubling the amount previously known in the Gemmula genus. In contrast to a predominantly shell-based morphological approach, which over the last 30 years proposed only 13 new species names for the Turridae genus Gemmula, the integrative approach described here identified 27 novel species hypotheses not linked to available species names in the literature. The formalized strategy applied here outlines an effective and reproducible protocol for large-scale species delimitation of hyperdiverse groups.

Genetic population structure in the Antarctic benthos: insights from the widespread amphipod, Orchomenella franklini

Currently there is very limited understanding of genetic population structure in the Antarctic benthos. We conducted one of the first studies of microsatellite variation in an Antarctic benthic invertebrate, using the ubiquitous amphipod Orchomenella franklini (Walker, 1903). Seven microsatellite loci were used to assess genetic structure on three spatial scales: sites (100 s of metres), locations (1-10 kilometres) and regions (1000 s of kilometres) sampled in East Antarctica at Casey and Davis stations. Considerable genetic diversity was revealed, which varied between the two regions and also between polluted and unpolluted sites. Genetic differentiation among all populations was highly significant (F(ST) = 0.086, R(ST) = 0.139, p<0.001) consistent with the brooding mode of development in O. franklini. Hierarchical AMOVA revealed that the majority of the genetic subdivision occurred across the largest geographical scale, with N(e)m≈1 suggesting insufficient gene flow to prevent independent evolution of the two regions, i.e., Casey and Davis are effectively isolated. Isolation by distance was detected at smaller scales and indicates that gene flow in O. franklini occurs primarily through stepping-stone dispersal. Three of the microsatellite loci showed signs of selection, providing evidence that localised adaptation may occur within the Antarctic benthos. These results provide insights into processes of speciation in Antarctic brooders, and will help inform the design of spatial management initiatives recently endorsed for the Antarctic benthos.

Global diversity of brittle stars (Echinodermata: Ophiuroidea)

This review presents a comprehensive overview of the current status regarding the global diversity of the echinoderm class Ophiuroidea, focussing on taxonomy and distribution patterns, with brief introduction to their anatomy, biology, phylogeny, and palaeontological history. A glossary of terms is provided. Species names and taxonomic decisions have been extracted from the literature and compiled in The World Ophiuroidea Database, part of the World Register of Marine Species (WoRMS). Ophiuroidea, with 2064 known species, are the largest class of Echinodermata. A table presents 16 families with numbers of genera and species. The largest are Amphiuridae (467), Ophiuridae (344 species) and Ophiacanthidae (319 species). A biogeographic analysis for all world oceans and all accepted species was performed, based on published distribution records. Approximately similar numbers of species were recorded from the shelf (n = 1313) and bathyal depth strata (1297). The Indo-Pacific region had the highest species richness overall (825 species) and at all depths. Adjacent regions were also relatively species rich, including the North Pacific (398), South Pacific (355) and Indian (316) due to the presence of many Indo-Pacific species that partially extended into these regions. A secondary region of enhanced species richness was found in the West Atlantic (335). Regions of relatively low species richness include the Arctic (73 species), East Atlantic (118), South America (124) and Antarctic (126).

Did vicariance and adaptation drive cryptic speciation and evolution of brooding in Ophioderma longicauda (Echinodermata: Ophiuroidea), a common Atlanto-Mediterranean ophiuroid?

Over the last decade, cryptic speciation has been discovered in an increasing number of taxa. Species complexes are useful models for the understanding of speciation processes. Motivated by the discovery of brooding specimens in the common Atlanto-Mediterranean broadcast spawning brittle star, Ophioderma longicauda, a recent study revealed the occurrence of divergent mitochondrial lineages. We analysed 218 specimens from 23 locations spread over the geographic range of the species with partial Cytochrome c Oxidase subunit I (COI) sequences. A subset of this sample was also surveyed with the internal transcribed spacer of the ribosomal DNA cluster (nuclear ITS-1). Our study revealed six highly divergent mitochondrial lineages, and the ITS-1 data confirmed that they most likely represent a species complex. Geographic ranges, abundances and genetic structures are contrasted among the putative cryptic species. Lineages in which brooding specimens have been found form a monophyletic group and are restricted to the Eastern Mediterranean basin, an oligotrophic zone. A phylogeny-trait association analysis revealed a phylogenetic signal for low 'chlorophyll a' values (our proxy for oligotrophy). An ecological shift related to the hyper oligotrophy of the Eastern Mediterranean region is therefore likely to have played a role in the evolution of brooding. This study revealed that a complex mixture of vicariance, population expansion, adaptive divergence and possibly high local diversification rates resulting from brooding has shaped the evolution of this species complex. The dating analysis showed that these events probably occurred in the Pleistocene epoch.

An efficient method to find potentially universal population genetic markers, applied to metazoans

Background

Despite the impressive growth of sequence databases, the limited availability of nuclear markers that are sufficiently polymorphic for population genetics and phylogeography and applicable across various phyla restricts many potential studies, particularly in non-model organisms. Numerous introns have invariant positions among kingdoms, providing a potential source for such markers. Unfortunately, most of the few known EPIC (Exon Primed Intron Crossing) loci are restricted to vertebrates or belong to multigenic families.

Results

In order to develop markers with broad applicability, we designed a bioinformatic approach aimed at avoiding multigenic families while identifying intron positions conserved across metazoan phyla. We developed a program facilitating the identification of EPIC loci which allowed slight variation in intron position. From the Homolens databases we selected 29 gene families which contained 52 promising introns for which we designed 93 primer pairs. PCR tests were performed on several ascidians, echinoderms, bivalves and cnidarians. On average, 24 different introns per genus were amplified in bilaterians. Remarkably, five of the introns successfully amplified in all of the metazoan genera tested (a dozen genera, including cnidarians). The influence of several factors on amplification success was investigated. Success rate was not related to the phylogenetic relatedness of a taxon to the groups that most influenced primer design, showing that these EPIC markers are extremely conserved in animals.

Conclusions

Our new method now makes it possible to (i) rapidly isolate a set of EPIC markers for any phylum, even outside the animal kingdom, and thus, (ii) compare genetic diversity at potentially homologous polymorphic loci between divergent taxa.