First molecular phylogeny of the circumtropical bivalve family Pinnidae (Mollusca, Bivalvia): evidence for high levels of cryptic species diversity

Reconstructing the Evolutionary History of Pinna nobilis: New Genetic Signals from the Past of a Species on the Brink of Extinction

Pinna nobilis, commonly known as the noble pen shell, is a marine bivalve endemic to the Mediterranean Sea. Unfortunately, due to a multifactorial disease that began affecting its populations in 2016, the species is currently facing the threat of extinction. To gain insights into the evolutionary history of P. nobilis before the mass mortality event (MME), and to obtain a comprehensive understanding of how evolutionary processes led to the adaptation of the species into the Mediterranean Sea, phylogenetic and phylogeographic analyses were carried out. The dataset analysed includes 469 sequences of COI gene fragment both from GenBank and the present study (100). The analysis performed evidenced that P. nobilis diverged about 2.5 mya, after the entrance of its ancestor into the Mediterranean Sea following the Zanclean flood (5.33 mya). Moreover, our results suggest that the starting point of colonisation was the central part of the western Mediterranean basin, with the eastern basin being populated subsequently. From a conservational viewpoint, these results provide important hints for present and future restocking plans, helping to reconstruct the pre-existing genetic variability in sites where the species became extinct.

Replacing mechanical protection with colorful faces-twice: parallel evolution of the non-operculate marine worm-snail genera Thylacodes (Guettard, 1770) and Cayo n. gen. (Gastropoda: Vermetidae)

Vermetid worm-snails are sessile and irregularly coiled marine mollusks common in warmer nearshore and coral reef environments that are subject to high predation pressures by fish. Often cryptic, some have evolved sturdy shells or long columellar muscles allowing quick withdrawal into better protected parts of the shell tube, and most have variously developed opercula that protect and seal the shell aperture trapdoor-like. Members of Thylacodes (previously: Serpulorbis) lack such opercular protection. Its species often show polychromatic head-foot coloration, and some have aposematic coloration likely directed at fish predators. A new polychromatic species, Thylacodes bermudensis n. sp., is described from Bermuda and compared morphologically and by DNA barcode markers to the likewise polychromatic western Atlantic species T. decussatus (Gmelin, 1791). Operculum loss, previously assumed to be an autapomorphy of Thylacodes, is shown to have occurred convergently in a second clade of the family, for which a new genus Cayo n. gen. and four new western Atlantic species are introduced: C. margarita n. sp. (type species; with type locality in the Florida Keys), C. galbinus n. sp., C. refulgens n. sp., and C. brunneimaculatus n. sp. (the last three with type locality in the Belizean reef) (all new taxa authored by Bieler, Collins, Golding & Rawlings). Cayo n. gen. differs from Thylacodes in morphology (e.g., a protoconch that is wider than tall), behavior (including deep shell entrenchment into the substratum), reproductive biology (fewer egg capsules and eggs per female; an obliquely attached egg capsule stalk), and in some species, a luminous, "neon-like", head-foot coloration. Comparative investigation of the eusperm and parasperm ultrastructure also revealed differences, with a laterally flattened eusperm acrosome observed in two species of Cayo n. gen. and a spiral keel on the eusperm nucleus in one, the latter feature currently unique within the family. A molecular phylogenetic analysis based on mitochondrial and nuclear rRNA gene sequences (12SrRNA, trnV, 16SrRNA, 28SrRNA) strongly supports the independent evolution of the two non-operculate lineages of vermetids. Thylacodes forms a sister grouping to a clade comprising Petaloconchus, Eualetes, and Cupolaconcha, whereas Cayo n. gen is strongly allied with the small-operculate species Vermetus triquetrus and V. bieleri. COI barcode markers provide support for the species-level status of the new taxa. Aspects of predator avoidance/deterrence are discussed for these non-operculate vermetids, which appear to involve warning coloration, aggressive behavior when approached by fish, and deployment of mucous feeding nets that have been shown, for one vermetid in a prior study, to contain bioactive metabolites avoided by fish. As such, non-operculate vermetids show characteristics similar to nudibranch slugs for which the evolution of warning coloration and chemical defenses has been explored previously.

Hidden species diversity and mito-nuclear discordance within the Mediterranean cone snail, Lautoconus ventricosus

The Mediterranean cone snail, Lautoconus ventricosus, is currently considered a single species inhabiting the whole Mediterranean basin and the adjacent Atlantic coasts. Yet, no population genetic study has assessed its taxonomic status. Here, we collected 245 individuals from 75 localities throughout the Mediterranean Sea and used cox1 barcodes, complete mitochondrial genomes, and genome skims to test whether L. ventricosus represents a complex of cryptic species. The maximum likelihood phylogeny based on complete mitochondrial genomes recovered six main clades (hereby named blue, brown, green, orange, red, and violet) with sufficient sequence divergence to be considered putative species. On the other hand, phylogenomic analyses based on 437 nuclear genes only recovered four out of the six clades: blue and orange clades were thoroughly mixed and the brown one was not recovered. This mito-nuclear discordance revealed instances of incomplete lineage sorting and introgression, and may have caused important differences in the dating of main cladogenetic events. Species delimitation tests proposed the existence of at least three species: green, violet, and red+blue+orange (i.e., cyan). Green plus cyan (with sympatric distributions) and violet, had West and East Mediterranean distributions, respectively, mostly separated by the Siculo-Tunisian biogeographical barrier. Morphometric analyses of the shell using species hypotheses as factor and shell length as covariate showed that the discrimination power of the studied parameters was only 70.2%, reinforcing the cryptic nature of the uncovered species, and the importance of integrative taxonomic approaches considering morphology, ecology, biogeography, and mitochondrial and nuclear population genetic variation.

Blue mussels of the Mytilus edulis species complex from South America: The application of species delimitation models to DNA sequence variation

Smooth-shelled blue mussels, Mytilus spp., have a worldwide antitropical distribution and are ecologically and economically important. Mussels of the Mytilus edulis species complex have been the focus of numerous taxonomic and biogeographical studies, in particular in the Northern hemisphere, but the taxonomic classification of mussels from South America remains unclear. The present study analysed 348 mussels from 20 sites in Argentina, Chile, Uruguay and the Falkland Islands on the Atlantic and Pacific coasts of South America. We sequenced two mitochondrial locus, Cytochrome c Oxidase subunit I (625 bp) and 16S rDNA (443 bp), and one nuclear gene, ribosomal 18S rDNA (1770 bp). Mitochondrial and nuclear loci were analysed separately and in combination using maximum likelihood and Bayesian inference methods to identify the combination of the most informative dataset and model. Species delimitation using five different models (GMYC single, bGMYC, PTP, bPTP and BPP) revealed that the Mytilus edulis complex in South America is represented by three species: native M. chilensis, M. edulis, and introduced Northern Hemisphere M. galloprovincialis. However, all models failed to delimit the putative species Mytilus platensis. In contrast, however, broad spatial scale genetic structure in South America using Geneland software to analyse COI sequence variation revealed a group of native mussels (putatively M. platensis) in central Argentina and the Falkland Islands. We discuss the scope of species delimitation methods and the use of nuclear and mitochondrial genetic data to the recognition of species within the Mytilus edulis complex at regional and global scales.

Phylogeography, colouration, and cryptic speciation across the Indo-Pacific in the sea urchin genus Echinothrix

The sea urchins Echinothrix calamaris and Echinothrix diadema have sympatric distributions throughout the Indo-Pacific. Diverse colour variation is reported in both species. To reconstruct the phylogeny of the genus and assess gene flow across the Indo-Pacific we sequenced mitochondrial 16S rDNA, ATPase-6, and ATPase-8, and nuclear 28S rDNA and the Calpain-7 intron. Our analyses revealed that E. diadema formed a single trans-Indo-Pacific clade, but E. calamaris contained three discrete clades. One clade was endemic to the Red Sea and the Gulf of Oman. A second clade occurred from Malaysia in the West to Moorea in the East. A third clade of E. calamaris was distributed across the entire Indo-Pacific biogeographic region. A fossil calibrated phylogeny revealed that the ancestor of E. diadema diverged from the ancestor of E. calamaris ~ 16.8 million years ago (Ma), and that the ancestor of the trans-Indo-Pacific clade and Red Sea and Gulf of Oman clade split from the western and central Pacific clade ~ 9.8 Ma. Time since divergence and genetic distances suggested species level differentiation among clades of E. calamaris. Colour variation was extensive in E. calamaris, but not clade or locality specific. There was little colour polymorphism in E. diadema.

Natural hybridization between pen shell species: Pinna rudis and the critically endangered Pinna nobilis may explain parasite resistance in P. nobilis

Recently, Pinna nobilis pen shells population in Mediterranean Sea has plummeted due to a Mass Mortality Event caused by an haplosporidian parasite. In consequence, this bivalve species has been included in the IUCN Red List as "Critically Endangered". In the current scenario, several works are in progress to protect P. nobilis from extinction, being identification of hybrids (P. nobilis x P. rudis) among survivors extremely important for the conservation of the species.Morphological characteristics and molecular analyses were used to identify putative hybrids. A total of 10 individuals of each species (P. nobilis and P. rudis) and 3 doubtful individuals were considered in this study. The putative hybrids showed shell morphology and mantle coloration intermingled exhibiting both P. nobilis and P. rudis traits. Moreover, the analyses of 1150 bp of the 28S gene showed 9 diagnostic sites between P. rudis and P. nobilis, whereas hybrids showed both parental diagnostic alleles at the diagnostic loci. Regarding the multilocus genotypes from the 8 microsatellite markers, the segregation of two Pinna species was clearly detected on the PCoA plot and the 3 hybrids showed intermediate positions.This is the first study evidencing the existence of hybrids P. nobilis x P. rudis, providing molecular methodology for a proper identification of new hybrids. Further studies testing systematically all parasite-resisting isolated P. nobilis should be undertaken to determine if the resistance is resulting from introgression of P. rudis into P. nobilis genome and identifying aspects related to resistance.

Survival and behavior patterns associated with hypoxia at different life stages of the pen shell Atrina cf. japonica

Pen shell (Atrina cf. japonica) resources have been devastated in Ariake Bay, Japan, and to facilitate the recovery of this species, there is an urgent need to fully understand the factors contributing to its high levels of mortality. Pen shells living in natural waters grow through successive life stages, and environmental factors may affect these bivalves differently at different stages. Accordingly, to elucidate the causes of mortality in natural waters, it is necessary to gain an understanding of the quantitative effects of environmental factors on bivalves at each life stage. In this study, we sought to determine the differential effects of hypoxic conditions on 1-year-old (average shell length: 93.9 mm) and 2-year-old (146.5 mm) pen shells bred under artificial conditions. We exposed shells of each age group to six different dissolved oxygen (DO) concentrations for 96 h and monitored their behavior and survival rate. Based on the survival status, we estimated the lethal DO concentrations that induced 50%, 95%, and 5% mortality (LC₅₀, LC₉₅, and LC₅, respectively) at each age. We found that for 1-year-old shells, the LC₅₀ values at 48, 72, and 96 h were 0.51, 0.74, and 0.84 mg/L, respectively, whereas the corresponding values for 2-year-old shells were 0.74, 1.27, and 1.80 mg/L. Furthermore, we found that for 1- and 2-year-old shells, the estimated ranges from LC₉₅ to LC₅ at 48, 72, and 96 h were 0.39-0.68, 0.62-0.88, and 0.64-1.12 mg/L, and 0.31-1.75, 0.77-2.09, and 1.29-2.53 mg/L, respectively. Under low DO concentrations (0.47 to 1.93 mg/L and 0.49 to 3.30 mg/L for 1- and 2-year-old shells, respectively), we observed pen shells with more than half of their shell length protruding above the substrate. In addition to age and body size, the 1- and 2-year-old pen shells used in the present study also differed with respect to reproductive status, with 7.6% of 1-year-old and 96.7% of 2-year-old shells considered to be fully ripe. Collectively, our observations indicate that 2-year-old pen shells are less tolerant to hypoxic conditions than are 1-year-old pen shells, and we suspect that the differences in hypoxic tolerance could be attributable to differences in the physiological status of the pen shells during gonadal development. We believe the findings of this study will make an important contribution to enhancing our understanding of the effects of hypoxia on the viability of A. cf. japonica in natural waters.

Uncovering the shell game with barcodes: diversity of meiofaunal Caecidae snails (Truncatelloidea, Caenogastropoda) from Central America

Caecidae is a species-rich family of microsnails with a worldwide distribution. Typical for many groups of gastropods, caecid taxonomy is largely based on overt shell characters. However, identification of species using shell characteristics is problematic due to their rather uniform, tubular shells, the presence of different growth stages, and a high degree of intraspecific variability. In the present study, a first integrative approach to caecid taxonomy is provided using light-microscopic investigation with microsculptural analyses and multi-marker barcoding, in conjunction with molecular species delineation analyses (ABGD, haplotype networks, GMYC, and bPTP). In total 132 specimens of Caecum and Meioceras collected during several sampling trips to Central America were analyzed and delineated into a minimum of 19 species to discuss putative synonyms, and supplement the original descriptions. Molecular phylogenetic analyses suggest Meiocerasnitidum and M.cubitatum should be reclassified as Caecum, and the genus Meioceras might present a junior synonym of Caecum. Meiofaunal caecids morphologically resembling C.glabrum from the Northeast Atlantic are a complex of cryptic species with independent evolutionary origins, likely associated with multiple habitat shifts to the mesopsammic environment. Caecuminvisibile Egger & Jörger, sp. nov. is formally described based on molecular diagnostic characters. This first integrative approach towards the taxonomy of Caecidae increases the known diversity, reveals the need for a reclassification of the genus Caecum and serves as a starting point for a barcoding library of the family, thereby enabling further reliable identifications of these taxonomically challenging microsnails in future studies.

Development of 26 highly polymorphic microsatellite markers for the highly endangered fan mussel Pinna nobilis and cross-species amplification

The fan mussel, Pinna nobilis is a highly endangered bivalve species endemic to the Mediterranean Sea. During the last few decades, populations have been greatly reduced due to anthropic impacts and they are now under strict protection in most Mediterranean countries. Today, the species is facing a major crisis following the introduction of an haplosporidan parasite which is driving mass mortality in almost all P. nobilis populations throughout the Mediterranean Sea. Gathering additional knowledge regarding dynamics and connectivity patterns of P. nobilis populations is now more than ever critical. Here, we describe the development of 26 highly polymorphic microsatellite markers. Average allelic diversity of 10.9 alleles per locus was reported and heterozygosity ranged from 0.0294 to 0.9737. We tested cross-species amplification in four Pinna species for the new markers together with 10 already published markers, and analysed its success according to the genetic distances among species. Cross-species transferability success ranged from 3 to 38% and had a negative relationship with the genetic distance between the target species and the tested species. The establishment of this new set of high-resolution markers provides a useful tool to understand processes driving gene flow and genetic diversity in P. nobilis populations and the closest congeneric species.

A mycobacterial disease is associated with the silent mass mortality of the pen shell Pinna nobilis along the Tyrrhenian coastline of Italy

Disease is an increasing threat for marine bivalves worldwide. Recently, a mass mortality event (MME) impacting the bivalve Pinna nobilis was detected across a wide geographical area of the Spanish Mediterranean Sea and linked to a haplosporidian parasite. In 2017–2018, mass mortality events affecting the pen shell Pinna nobilis were recorded in two different regions of Italy, Campania and Sicily, in the Tyrrhenian Sea (Mediterranean Sea). Histopathological and molecular examinations of specimens showed the presence of Haplosporidium sp. in only one specimen in one area. Conversely, in all of the surveyed moribund animals, strong inflammatory lesions at the level of connective tissue surrounding the digestive system and gonads and linked to the presence of intracellular Zhiel-Neelsen-positive bacteria were observed. Molecular analysis of all of the diseased specimens (13) confirmed the presence of a Mycobacterium. Blast analysis of the sequences from all of the areas revealed that they were grouped together with the human mycobacterium M. sherrisii close to the group including M. shigaense, M. lentiflavum and M. simiae. Based on pathological and molecular findings, it is proposed that a mycobacterial disease is associated with the mortality episodes of Pinna nobilis, indicating that, at this time, Haplosporidium sp. is not responsible for these events in Campanian and Sicilian waters.

Morphological and textural evolution of the prismatic ultrastructure in mollusc shells: A comparative study of Pinnidae species

Molluscan shells, exhibiting a variety of complex three-dimensional architectures, are an exemplar model system to study biogenic mineral formation by living organisms. Recent studies have demonstrated that the deposition process of some shell ultrastructures can be described using classical analytical models borrowed from materials physics, which were developed to predict the structural evolution of man-made and geological polycrystalline composite assemblies. In the current study, we use this newly developed capacity to quantitatively describe the morphogenesis of the prismatic ultrastructure in three shells from the bivalve family Pinnidae towards establishing a correlation between structure, texture, growth kinetics, topology and phylogeny of the species. Using data collected by electron microscopy, synchrotron-based microtomography, electron backscatter diffraction analysis (EBSD) and X-ray diffraction we demonstrate that the prismatic ultrastructures in Pinnidae are formed following either ideal or triple-junction-controlled kinetics, which are shown to be closely linked to the morphological and topological characteristics, as well as crystallographic texture of these biocomposites. The experimental and analytical framework presented in this comparative study can serve as an additional tool for classifying molluscan shell ultrastructures on the levels of structural and textural morphogenesis. STATEMENT OF SIGNIFICANCE: The ability to quantitatively describe the structural evolution of the prismatic architecture in mollusc shells is used for the first time to derive and compare between analytical parameters that define the growth kinetics and morphological and topological evolution during the growth of three shells from the family Pinnidae from two different genera. Furthermore, these parameters are linked to the evolution of crystallographic texture in the studied architectures. The developed experimental and analytical framework not only enables us to quantitatively describe species-specific growth mechanisms but also suggests a direct correlation between the evolution of morphology and texture.

Speciation among sympatric lineages in the genus Palythoa (Cnidaria: Anthozoa: Zoantharia) revealed by morphological comparison, phylogenetic analyses and investigation of spawning period

Zoantharians are sessile marine invertebrates and colonial organisms possessing sexual and asexual reproductive ability. The zooxanthellate zoantharian genus Palythoa is widely distributed in coral reef ecosystems. In the Ryukyu Archipelago, Japan, sympatric Palythoa tuberculosa and P. mutuki are the dominant species of this genus in the intertidal zone. Previous phylogenetic analyses have shown that these two species are closely related, and additionally revealed a putative sympatric hybrid species (designated as Palythoa sp. yoron). In this study, we attempted to delineate Palythoa species boundaries and to clarify the relationships among these three groups plus another additional putative sympatric species (P. aff. mutuki) by multiple independent criteria. The morphology of these four lineages was clearly different; for example the number of tentacles was significantly different for each species group in all pairwise comparisons. From observations of gonadal development conducted in 2010 and 2011, P. sp. yoron and P. aff. mutuki appear to be reproductively isolated from P. tuberculosa. In the phylogenetic tree resulting from maximum likelihood analyses of the ITS-rDNA sequence alignment, P. tuberculosa and P. sp. yoron formed a very well supported monophyletic clade (NJ = 100%, ML = 95%, Bayes = 0.99). This study demonstrates that despite clear morphological and/or reproductive differences, P. tuberculosa and P. sp. yoron are phylogenetically entangled and closely related to each other, as are P. mutuki and P. aff. mutuki. Additionally, no single molecular marker was able to divide these four lineages into monophyletic clades by themselves, and a marker that has enough resolution to solve this molecular phylogenetic species complex is required. In summary, the morphological and reproductive results suggest these lineages are four separate species, and that incomplete genetic lineage sorting may prevent the accurate phylogenetic detection of distinct species with the DNA markers utilized in this study, demonstrating the value of morphological and reproductive data when examining closely related lineages.

Genetic and oceanographic tools reveal high population connectivity and diversity in the endangered pen shell Pinna nobilis

For marine meta-populations with source-sink dynamics knowledge about genetic connectivity is important to conserve biodiversity and design marine protected areas (MPAs). We evaluate connectivity of a Mediterranean sessile species, Pinna nobilis. To address a large geographical scale, partial sequences of cytochrome oxidase I (COI, 590 bp) were used to evaluate phylogeographical patterns in the Western Mediterranean, and in the whole basin using overlapping sequences from the literature (243 bp). Additionally, we combined (1) larval trajectories based on oceanographic currents and early life-history traits and (2) 10 highly polymorphic microsatellite loci collected in the Western Mediterranean. COI results provided evidence for high diversity and low inter-population differentiation. Microsatellite genotypes showed increasing genetic differentiation with oceanographic transport time (isolation by oceanographic distance (IBD) set by marine currents). Genetic differentiation was detected between Banyuls and Murcia and between Murcia and Mallorca. However, no genetic break was detected between the Balearic populations and the mainland. Migration rates together with numerical Lagrangian simulations showed that (i) the Ebro Delta is a larval source for the Balearic populations (ii) Alicante is a sink population, accumulating allelic diversity from nearby populations. The inferred connectivity can be applied in the development of MPA networks in the Western Mediterranean.

Unifying latitudinal gradients in range size and richness across marine and terrestrial systems

Many marine and terrestrial clades show similar latitudinal gradients in species richness, but opposite gradients in range size-on land, ranges are the smallest in the tropics, whereas in the sea, ranges are the largest in the tropics. Therefore, richness gradients in marine and terrestrial systems do not arise from a shared latitudinal arrangement of species range sizes. Comparing terrestrial birds and marine bivalves, we find that gradients in range size are concordant at the level of genera. Here, both groups show a nested pattern in which narrow-ranging genera are confined to the tropics and broad-ranging genera extend across much of the gradient. We find that (i) genus range size and its variation with latitude is closely associated with per-genus species richness and (ii) broad-ranging genera contain more species both within and outside of the tropics when compared with tropical- or temperate-only genera. Within-genus species diversification thus promotes genus expansion to novel latitudes. Despite underlying differences in the species range-size gradients, species-rich genera are more likely to produce a descendant that extends its range relative to the ancestor's range. These results unify species richness gradients with those of genera, implying that birds and bivalves share similar latitudinal dynamics in net species diversification.

Human Stressors Are Driving Coastal Benthic Long-Lived Sessile Fan Mussel Pinna nobilis Population Structure More than Environmental Stressors

Coastal degradation and habitat disruption are severely compromising sessile marine species. The fan shell Pinna nobilis is an endemic, vulnerable species and the largest bivalve in the Mediterranean basin. In spite of species legal protection, fan shell populations are declining. Models analyzed the contributions of environmental (mean depth, wave height, maximum wave height, period of waves with high energy and mean direction of wave source) versus human-derived stressors (anchoring, protection status, sewage effluents, fishing activity and diving) as explanatory variables depicting Pinna nobilis populations at a mesoscale level. Human stressors were explaining most of the variability in density spatial distribution of fan shell, significantly disturbing benthic communities. Habitat protection affected P. nobilis structure and physical aggression by anchoring reveals a high impact on densities. Environmental variables instead played a secondary role, indicating that global change processes are not so relevant in coastal benthic communities as human-derived impacts.