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

Scientific history, biogeography, and biological traits predict presence of cryptic or overlooked species

Genetic data show that many nominal species are composed of more than one biological species, and thus contain cryptic species in the broad sense (including overlooked species). When ignored, cryptic species generate confusion which, beyond biodiversity or vulnerability underestimation, blurs our understanding of ecological and evolutionary processes and may impact the soundness of decisions in conservation or medicine. However, very few hypotheses have been tested about factors that predispose a taxon to contain cryptic or overlooked species. To fill this gap, we surveyed the literature on free-living marine metazoans and built two data sets, one of 187,603 nominal species and another of 83 classes or phyla, to test several hypotheses, correcting for sequence data availability, taxon size and phylogenetic relatedness. We found a strong effect of scientific history: the probability of a taxon containing cryptic species was highest for the earliest described species and varied among time periods potentially consistently with an influence of prevailing scientific theories. The probability of cryptic species being present was also increased for species with large distribution ranges. They were more frequent in the north polar and south polar zones, contradicting previous predictions of more cryptic species in the tropics, and supporting the hypothesis that many cryptic species diverged recently. The number of cryptic species varied among classes, with an excess in hydrozoans and polychaetes, and a deficit in actinopterygians, for example, but precise class ranking was relatively sensitive to the statistical model used. For all models, biological traits, rather than phylum, appeared responsible for the variation among classes: there were fewer cryptic species than expected in classes with hard skeletons (perhaps because they provide good characters for taxonomy) and image-forming vision (in which selection against heterospecific mating may enhance morphological divergence), and more in classes with internal fertilisation. We estimate that among marine free-living metazoans, several thousand additional cryptic species complexes could be identified as more sequence data become available. The factors identified as important for marine animal cryptic species are likely important for other biomes and taxa and should aid many areas in biology that rely on accurate species identification.

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.

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.

Dwarf brooder versus giant broadcaster: combining genetic and reproductive data to unravel cryptic diversity in an Antarctic brittle star

Poecilogony, or multiple developmental modes in a single species, is exceedingly rare. Several species described as poecilogenous were later demonstrated to be multiple (cryptic) species with a different developmental mode. The Southern Ocean is known to harbor a high proportion of brooders (Thorson's Rule) but with an increasing number of counter examples over recent years. Here we evaluated poecilogony vs. crypticism in the brittle star Astrotoma agassizii across the Southern Ocean. This species was initially described from South America as a brooder before some pelagic stages were identified in Antarctica. Reproductive and mitochondrial data were combined to unravel geographic and genetic variation of developmental modes. Our results indicate that A. agassizii is composed of seven well-supported and deeply divergent clades (I: Antarctica and South Georgia; II: South Georgia and Sub-Antarctic locations including Kerguelen, Patagonian shelf, and New Zealand; III-VI-VII: Patagonian shelf, IV-V: South Georgia). Two of these clades demonstrated strong size dimorphism when in sympatry and can be linked to differing developmental modes (Clade V: dwarf brooder vs. Clade I: giant broadcaster). Based on their restricted geographic distributions and on previous studies, it is likely that Clades III-VI-VII are brooders. Clade II is composed of different morphological species, A. agassizii and A. drachi, the latter originally used as the outgroup. By integrating morphology, reproductive, and molecular data we conclude that the variation identified in A. agassizii is best described as crypticism rather than poecilogony.

Abbreviated Development of the Brooding Brittle Star Ophioplocus esmarki

The bilaterally symmetrical, feeding larval stage is an ancestral condition in echinoderms. However, many echinoderms have evolved abbreviated development and form a pentamerous juvenile without a feeding larva. Abbreviated development with a non-feeding vitellaria larva is found in five families of brittle stars, but very little is known about this type of development. In this study, the external anatomy, ciliary bands, neurons, and muscles were examined in the development of the brooded vitellaria larva of Ophioplocus esmarki. The external morphology throughout development shows typical vitellaria features, including morphogenetic movements to set up the vitellaria body plan, an anterior preoral lobe, a posterior lobe, transverse ciliary bands, and development of juvenile structures on the mid-ventral side. An early population of neurons forms at the base of the preoral lobe at the pre-vitellaria stage after the initial formation of the coelomic cavities. These early neurons may be homologous to the apical neurons that develop in echinoderms with feeding larval forms. Neurons form close to the ciliary bands, but the vitellaria larva lacks the tracts of neurons associated with the ciliary bands found in echinoderms with feeding larvae. Additional neurons form in association with the axial complex and persist into the juvenile stage. Juvenile nerves and muscles form with pentamerous symmetry in the late vitellaria stage in a manner similar to their development within the late ophiopluteus larva. Even though O. esmarki is a brooding brittle star, its developmental sequence retains the general vitellaria shape and structure; however, the vitellaria larvae are unable to swim in the water column.

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.

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.

Influence of the larval phase on connectivity: strong differences in the genetic structure of brooders and broadcasters in the Ophioderma longicauda species complex

Closely related species with divergent life history traits are excellent models to infer the role of such traits in genetic diversity and connectivity. Ophioderma longicauda is a brittle star species complex composed of different genetic clusters, including brooders and broadcasters. These species diverged very recently and some of them are sympatric and ecologically syntopic, making them particularly suitable to study the consequences of their trait differences. At the scale of the geographic distribution of the broadcasters (Mediterranean Sea and northeastern Atlantic), we sequenced the mitochondrial marker COI and genotyped an intron (i51) for 788 individuals. In addition, we sequenced 10 nuclear loci newly developed from transcriptome sequences, for six sympatric populations of brooders and broadcasters from Greece. At the large scale, we found a high genetic structure within the brooders (COI: 0.07 < F(ST) < 0.65) and no polymorphism at the nuclear locus i51. In contrast, the broadcasters displayed lower genetic structure (0 < F(ST) < 0.14) and were polymorphic at locus i51. At the regional scale, the multilocus analysis confirmed the contrasting genetic structure between species, with no structure in the broadcasters (global F(ST) < 0.001) and strong structure in the brooders (global F(ST) = 0.49), and revealed a higher genetic diversity in broadcasters. Our study showed that the lecithotrophic larval stage allows on average a 50-fold increase in migration rates, a 280-fold increase in effective size and a threefold to fourfold increase in genetic diversity. Our work, investigating complementary genetic markers on sympatric and syntopic taxa, highlights the strong impact of the larval phase on connectivity and genetic diversity.