Research Article
Print
Research Article
A new species of the genus Yoldiella (Bivalvia, Protobranchia, Yoldiidae) from Haima Cold Seep, South China Sea, China
expand article infoQi Gao§, Yan Tang§, Junlong Zhang§|
‡ Qingdao Agricultural University, Qingdao, China
§ Institute of Oceanology, Chinese Academy of Sciences, Qingdao, China
| Marine Biological Museum, Chinese Academy of Sciences, Qingdao, China
¶ University of Chinese Academy of Sciences, Beijing, China
Open Access

Abstract

In present study, a previously unidentified but frequently encountered species of deep-sea protobranch, Yoldiella haimaensis sp. nov., is described new to science from the Haima Cold Seep on the northwestern slope of the South China Sea. A morphological analysis confirmed that this species belongs to a previously undescribed species of the genus Yoldiella A.E. Verrill & K.J. Bush, 1897. It differs morphologically from other known species within the genus in its shell shape, degree of inflation, beaks, and number of hinge teeth. Furthermore, we sequenced three gene segments of Y. haimaensis sp. nov., comprising a nuclear ribosomal gene (18S rRNA), a nuclear protein-coding gene (histone H3), and a mitochondrial gene (cytochrome c oxidase subunit I, COI). Our phylogenetic analysis performed on the superfamily Nuculanoidea and family Yoldiidae indicates that the genus Yoldiella is non-monophyletic, and the widely recognized families within the superfamily Nuculanoidea are also not monophyletic. Our results provide molecular insights into the Protobranchia and highlight the necessity for further samples and data to revise the classification of families and genera within the superfamily using an integrative approach that combines morphological analysis and molecular data.

Key words

Anatomy, chemosynthetic ecosystems, deep-sea, molecular analysis, morphology

Introduction

Protobranchia, with a significant evolutionary history dating back to the Cambrian, represents an ancestral and basal group of Bivalvia. The protobranchs are primarily found in the subtidal zone, especially in the deep sea, and are generally deposit feeders that bury themselves in the soft sediment (Allen 1978). They have limited presence in the intertidal zone. So, it is difficult to collect specimens of this group (Xu 1999). The highly conserved and distinctive morphology and anatomy, including gill structure, hinge conformation, shell microstructure, as well as the pericalymma larval development, small in size and lifestyle of the group make it a difficult but fascinating taxon of Bivalvia (Zardus 2002). The simplicity in the form of protobranch bivalves veils the complexity of their phylogeny. The monophyly of Protobranchia has been discussed intensively and become a subject of controversy due to the extensive use of molecular methods (Smith et al. 2011; Sharma et al. 2012; Bieler et al. 2014; González et al. 2015; Combosch et al. 2017; Lemer et al. 2019). Phylogenetic analyses using four nuclear genes (Sharma et al. 2012) and an exemplar-based approach combining Sanger-based sequences and an extensive morphological data matrix (Bieler et al. 2014), as well as the phylogenomic analysis using genomes and transcriptomes (González et al. 2015), recovered the monophyly of Protobranchia. However, the subsequent analysis by Combosch et al. (2017), which utilized five genes and included more taxa, supported the polyphyly of Protobranchia. The latest research has divided Protobranchia into three orders and five superfamilies: Nuculanoidea and Sareptoidea in Nuculanida, Solemyoidea and Manzanelloidea in Solemyida, and Nuculoidea in the order Nuculida (Sharma et al. 2013; Sato et al. 2020). The reconstruction of the phylogeny has indicated that eight families (i.e. Nuculanidae, Bathyspinulidae, Malletiidae, Neilonellidae, Phaseolidae, Siliculidae, Tindariidae, and Yoldiidae) within the superfamily Nuculanoidea are all non-monophyletic (Sato et al. 2020). Shell microstructure of protobranchs play a crucial role in their classification at the superfamily level. Moreover, subtle differences in shell microstructure can aid in distinguishing similar species (Sato et al. 2020). Extensive research has revealed a multitude of unknown or cryptic protobranch species awaiting description and classification under the integrative taxonomy framework (e.g. Neulinger et al. 2006; Zardus et al. 2006). Additional samples and data are required to revise the families and genera within the superfamily Nuculanoidea through a combination of morphological diagnosis and molecular analysis.

Cold seeps are natural phenomena widely distributed across the globe. On the seafloor of these areas, the hydrocarbon-rich fluids and gases leak from cracks and enter the water column through sediment, forming a distinctive habitat (Dong et al. 2021). In March 2015, a newly active cold seep was discovered using the Haima remotely operated vehicle (ROV) on the northwestern slope of the South China Sea (SCS) (Zhao et al. 2020). Dong et al. (2021) documented 34 epibenthic macrofauna species collected from Haima Cold Seep with 24 species being identified. Yao et al. (2022) identified 12 macrobenthic species from five phyla and 12 families in the Haima Cold Seep, including two species first found in this location. Seven phyla, 14 classes, and 65 species were identified by He et al. (2023) in this cold seep. To date, more than 80 species of macrobenthic organisms have been collected from the Haima cold seeps (Wang et al. 2022).

This study presents the description of a new Yoldiella species, Y. haimaensis sp. nov., from the Haima Cold Seep in the SCS. This species had been identified as Malletia sp. or Yoldiella sp. (Dong et al. 2021; Ke et al. 2022; He et al. 2023). Additionally, we provide the sequences of three gene segments of the new species, including a nuclear ribosomal gene (18S rRNA), a nuclear protein-coding gene (histone H3), and a mitochondrial gene (cytochrome c oxidase subunit 1, COI). A phylogenetic analysis was conducted on the superfamily Nuculanoidea and family Yoldiidae, providing molecular data for the study of Protobranchia and enhancing understanding of macrobenthos at the Haima Cold Seep.

Materials and methods

Specimen collection and identification

The specimens were collected on July 4–12, 2022 from the Haima Cold Seep at a water depth of 1390 m in the SCS (16°43'N, 110°28'E) (Fig. 1) using the ROV and TV Grab of Research Vessel Kexue. All specimens were fixed in 100% ethanol and deposited at the Marine Biological Museum (MBM), Chinese Academy of Sciences, Qingdao, China.

Figure 1. 

Map of samples site (triangle) of Yoldiella haimaensis sp. nov.

DNA extraction, amplification, and sequencing

Two specimens with tissues were randomly selected from the specimens collected. DNA was extracted from the muscle tissues using the TIANamp Marine Animals DNA Kit. Three gene fragments were amplified, including a nuclear ribosomal gene (18S rRNA), a nuclear protein-coding gene (histone H3), and a mitochondrial gene (cytochrome c oxidase subunit 1, COI), which were subsequently utilized for phylogenetic analyses. The PCR program for the mitochondrial genes was as follows: initial denaturation for 180 s at 94 °C, followed by 35 cycles of denaturation for 30 s at 94 °C, annealing for 45 s at 46 °C and elongation for 60 s at 72 °C. The final elongation step was conducted for 10 min at 72 °C. The PCR conditions for the 18S rRNA and H3 genes were performed the same as above, with the only difference being in the annealing step. The annealing temperature for H3 was 55 °C, and the annealing step for 18S rRNA was performed in a touch-down manner (Don et al. 1991) with an initial annealing temperature of 70 °C followed by a reduction of 1 °C per cycle until 65 °C. The primer sequences are listed in Suppl. material 1.

The PCR products were separated by electrophoresis using 1.0% agarose gels, purified, and then sent for sequencing to Sangon Biotech Co. Ltd.

Phylogenetic analyses

The sequences newly acquired in this study and obtained from NCBI are listed in Suppl. materials 2, 3. All sequences were aligned using MAFFT v. 7 software (Katoh and Standley 2013) based on the amino acid sequences and employing the Auto and G-INS-I algorithms. Subsequently, the aligned sequences in each dataset were manually trimmed to the same length. Ambiguously aligned sites in the ribosomal gene were removed using GBLOCKS v. 0.91b (Castresana 2000) with the least stringent settings. Because of the limited molecular data of Protobranchia, different combined gene datasets were used for phylogenetic analyses of the family Yoldiidae and superfamily Nuculanoidea. The final dataset for Yoldiidae phylogenetic trees comprised 1059 bp of 18S rRNA (88% of 1192 bp before Gblock), and 594 bp of COI. The final dataset for Nuculanoidea phylogenetic trees included 1101 bp of 18S rRNA (81% of 1352 bp before Gblock), 588 bp of COI, and 305 bp of H3. The COI, 18S rRNA, and H3 gene sequences from the same individuals were concatenated using SequenceMatrix software (Vaidya et al. 2011) to form a combined gene dataset. The Maximum likelihood (ML) and Bayesian inference (BI) analysis based on concatenated datasets of COI and 18S rRNA were used for phylogenetic analyses. The ML tree was conducted using IQ-TREE v. 2.2.0-Linux (Nguyen et al. 2015). The most suitable evolution model was found by ModelFinder (Kalyaanamoorthy et al. 2017) and adopted automatically to infer the ML tree. Bootstrap supports (BS) were calculated with 1,000 replicates to assess branch supports. Gene partition models chosen for IQ-TREE were 18S rRNA, K2P+I; COI, TIM+F+I+G4; and combined gene, TIM3+F+I+G4. jModelTest v. 2.1.10 (Darriba et al. 2012) was used to evaluate the best-fitting nucleotide substitution model and derive the optimal model of the Bayesian phylogenetic tree. According to the Akaike information criterion (AIC), the best model for each isolated gene was 18S rRNA, HKY+I; COI, GTR+I+G; and combined gene, HKY+I+G. A Bayesian-inference (BI) analysis was performed using MrBayes v. 3.2.6 (Sharma et al. 2013) and the best model of each dataset. The posterior probability (PP) was estimated using four chains running 1 million generations and sampled every 100 generations. The first 25% of sampled trees were discarded as burn-in. The results of ML and BI trees were visualized and rendered using Figtree v. 1.4.4 (http://tree.bio.ed.ac.uk/software/figtree/).

Species delimitation

A variety of species delimitation methods were employed to determine that the species described here is not conspecific with another, already known species of Yoldiidae. COI data were analyzed using the program Automated Barcode Gap Discovery (ABGD; Puillandre et al. 2012) and Assemble Species by Automatic Partitioning (ASAP; Puillandre et al. 2021), a method to build species partitions from single-locus sequence alignments. Single gene trees were analyzed by applying the Bayesian implementation of the Poisson Tree Processes model (bPTP; Zhang et al. 2013) at the web server of the Heidelberg Institute for Theoretical Studies, Germany (http://species.h-its.org/).

Results

Systematics

Subclass Protobranchia Pelseneer, 1889

Order Nuculanida J.G. Carter, D.C. Campbell & M.R. Campbell, 2000

Superfamily Nuculanoidea H. Adams & A. Adams, 1858 (1854)

Family Yoldiidae Dall, 1908

Yoldiella A.E. Verrill & K.J. Bush, 1897

Type species

Yoldiella lucida (Lovén, 1846) (International Commission on Zoological Nomenclature 1985: Opinion 1306) by original designation (Recent, North Atlantic).

Yoldiella haimaensis sp. nov.

Figs 2, 3, 4

Malletia sp.: Dong et al. 2021: 4, fig. 5b; Ke et al. 2022: 4, fig. 2h.

Yoldiella sp.: He et al. 2023: 6, fig. 2Q.

Type specimens

Holotype : MBM 229041: length 7.5 mm, width 3.2, height 5.1 mm. Paratypes: MBM 229042: length 7.4 mm, width 3.3 mm, height 5.0 mm; MBM 229043: length 6.5 mm, width 2.6 mm, height 4.4 mm; MBM229044: length 6.3 mm, width 2.5 mm, height 4.1 mm; MBM229045: length 7.6 mm, width 3.0 mm, height 4.9 mm.

Type locality

Haima Cold Seep (depth 1390 m, 16°43.00'N, 110°28.00'E), off southern Hainan Island, South China Sea.

Diagnosis

Yoldiella haimaensis sp. nov. differs morphologically from other known species within the genus in shell shape, degree of inflation, beaks, and number of hinge teeth. Diagnostic characteristics: shell small, ovate, inflated medially. Posterior end slight produced. Resilifer triangular, projecting. Beak rather lower than other species, suborthogyrate, and easily worn. Hinge plate narrow; posterior hinge plate smaller than anterior one, with taxodont teeth in two series; 17–19 anterior and 15–16 posterior teeth on hinge plate.

Figure 2. 

Yoldiella haimaensis sp. nov. A–D holotype, MBM229041. Scale bar: 2 mm (all at the same scale).

Description

Shell small, elongate, ovate in outline, moderately inflated, opaque, fragile, 2.2–8.2 mm long, W/L about 0.40; H/L about 0.66, usually subequivalve, inequilateral. Shell surface smooth, with numerous very fine, regular, and nearly isometric growth lines, without radial stria. Periostracum light brown and flaky. Umbo slightly posterior to middle, low, large, obscure, opisthogyrate, and easily worn. Antero-dorsal margin convex; anterior end broadly rounded, merging smoothly to ventral margin. Ventral margin slightly convex, with very shallow sinus at postero-ventral corner. Postero-dorsal margin oblique and then convex, descending to blunt posterior margin. Posterior end slight produced. Escutcheon and lunule obscure. ligament amphidetic, thin, short.

Figure 3. 

Yoldiella haimaensis sp. nov. A–D paratype 1, MBM229042 E–H paratype 2, MBM229043 I paratype 3, MBM229044 J paratype 4 MBM229045. Scale bar: 2 mm (all at the same scale).

Internal surface porcelaneous white. Hinge plate moderately broad, narrow below umbo, moderately long, and rather strong, with two chevron-shaped columns and moderately sized taxodont lateral teeth, about 17–19 anterior teeth, about 15–16 posterior teeth, interrupted by a triangular, projecting resilifer, and not extend beyond the inner limit of adductor muscles. Angle of about 140° between anterior and posterior hinge plates. Posterior hinge plate usually smaller than anterior. Resilium oblique and often obscure in dry preserved specimens. Adductor scar obscure to evident; triangular anterior adductor scar larger than droplet-shaped posterior adductor scar. Pallial sinus obsolete; pallial line usually entire.

Figure 4. 

Anatomy features of Yoldiella haimaensis sp. nov. aa, anterior adductor; c, ctenidium; f, foot; lp, labial palp; m, mantle; pa, posterior adductor; pp, Palp proboscis; dg, digestive gland. Scale bar: 2 mm.

Mantle large, thin, and opaque; anterior adductor crescent-shaped, twice or three times size of posterior. Ctenidium structure simple and lamellar, at posterior side parallel to the postero-dorsal shell margin. Labial palp size moderate, consisting of flat, paired lamellae on each side, with appendages of elongated palp proboscis. Foot muscular and large, with a regular series of nearly rectangular protrusions at margins, partially covered by labial palp. Siphons combined posteriorly.

Etymology

The species epithet “haimaensis” is Latin and means “from Haima”, which refers to the name of the cold seep where the specimens were collected.

Distribution

Currently known only from the Haima Cold Seep on the northwestern slope of the South China Sea.

Remarks

Yoldiella haimaensis sp. nov. differs morphologically from other known species of Yoldiella in its shell shape, degree of inflation, beak characteristics, and number of hinge teeth. Its beaks are lower than those of other species, suborthogyrate, and prone to wearing easily. This new species resembles the type species of Yoldiella, Y. lucida. However, Y. haimaensis sp. nov. differs from Y. lucida in having lower, suborthogyrate beaks that wear easily and a slightly rounded posterior end. Yoldiella haimaensis sp. nov. has more teeth (17–19 anterior teeth, 15–16 posterior teeth) than Y. lucida (8 teeth on each end), with more anterior teeth than posterior teeth. Another species closely resembling Y. haimaensis sp. nov. in outline is Yoldiella sagamiana T. Okutani & K. Fujikura, 2022 from Sagami Bay, but Y. sagamiana has a larger W/L ratio and fewer teeth (15 anterior teeth, 10 posterior teeth) than the new species, and Y. sagamiana also has more pointed beaks and finer commarginal cords and lines. The outline of the new species is similar to Yoldiella biguttata Allen, H. L. Sanders & F. Hannah, 1995 from the Guyana Basin, but Y. biguttata has the more prominent umbo and the anterior and posterior series are either equal or with the anterior series having one additional tooth (5–6 in the largest specimen).

Species delimitation

Our species delimitation using the ABGD, ASAP, and bPTP methods show slightly different results, but they all support that our two samples are the same new species, and distinct from others in the family Yoldiidae (Fig. 5). ABGD and bPTP delineated the data into 19 species, while ASAP delineated the data into 12 species. ASAP grouped Yoldiella nana (M. Sars, 1865), Yoldiella frigida (Torrell, 1859), Yoldiella inconspicua A.E. Verrill & K.J. Bush, 1898, and Yoldiella orcia (Dall, 1916) into one species; ASAP and bPTP grouped Yoldiella philippiana (Nyst, 1845) and Yoldiella propinqua (Leche, 1878), which clustered together, into one species; ABGD and bPTP grouped Yoldia sp. and Yoldia scissurata Dall, 1897 into one species; ABGD grouped Yoldia notabilis Yokoyama, 1922 and Yoldia johanni Dall, 1925 into one species; bPTP delineated two species for the two individuals of Y. notabilis; ASAP grouped Yoldia hyperborea (A. Gould, 1841), Yoldia sp., Y. scissurata, Y. notabilis, and Y. johanni into one species; ABGD delineated two species for the seven individuals of Portlandia arctica (Gray, 1824).

Figure 5. 

Phylogenetic tree inferred by Bayesian inference (BI) and maximum likelihood (ML), and species delimitation based on COI gene. Numbers adjacent to nodes refer to ML bootstrap scores and BI posterior probability (left, and right, respectively). The results of three species delimitation methods are shown on the right of phylogenetic tree (each rectangle represents one species).

Molecular phylogenies

Due to the limitation of molecular data of Protobranchia, only COI and 18S rRNA could meet the requirements of the phylogenetic tree construction of the family Yoldiidae, and their data were different, so these genes were used to construct two separate phylogenetic trees (Figs 5, 6). Sarepta speciosa A. Adams, 1860 from the closely related superfamily Sareptoidea was selected as the outgroup in our phylogenetic analyses. BI and ML analyses yielded similar results. Phylogenetic analyses based on COI and 18S rRNA both showed that Yoldiella is not monophyletic. However, there were some differences. In the COI phylogenetic analysis, Y. nana and Y. frigida clustered into a sister taxon and has a close relationship with Y. inconspicua and Y. orcia, successively, while in 18S rRNA phylogenetic analysis, Y. nana clustered with species of the genus Megayoldia A.E. Verrill & K.J. Bush, 1897, and was distantly related to Y. inconspicua and Y. orcia. In addition, the COI analysis revealed that species of the genus Yoldia Møller, 1842 were nested among species of Yoldiella. Conversely, in the phylogenetic analysis constructed using 18S rRNA, Yoldia clustered as a distinct branch. The phylogenetic status of Y. haimaensis sp. nov. was different owing to the incongruent results of separate molecular analyses. The phylogenetic analysis based on COI showed that the new species was positioned at the base of the clade including species of Yoldia (Fig. 5), while it was at the basal position of the clade containing all Yoldiella, Yoldia, and Megayoldia species in the 18S rRNA phylogenetic tree with low support values (Fig. 6).

Figure 6. 

Phylogenetic tree inferred by Bayesian inference (BI) and maximum likelihood (ML) based on 18S rRNA gene. Numbers adjacent to nodes refer to ML bootstrap scores and BI posterior probability (left, and right, respectively).

Phylogenetic analysis constructed using a combined gene dataset (COI+18S+H3) did not improve resolution as more taxa were added within the superfamily Nuculanoidea (Suppl. materials 4, 5). The BI and ML analyses recovered different topologies, and most nodes at high levels received very low support values. Furthermore, the results indicate that the family Yoldiidae is non-monophyletic, and that none of the widely recognized families (Nuculanidae, Neilonellidae, Malletiidae, Siliculidae, Phaseolidae, Tindariidae, and Bathyspinulidae) form monophyletic groups in this reconstruction.

Discussion

The genus Yoldiella is believed to be ubiquitous in all the world’s oceans, with a particularly high abundance in deep waters (Benaim and Absalão 2011a). Environmental differences are likely to have a significant impact on the weight of their shell and the rate of their growth (Reed et al. 2013). However, their small size and offshore habitat contribute to their rarity in collections. Prior to this study, no species of this genus had been identified in the South China Sea. The discovery of Y. haimaensis sp. nov. bridges this geographical gap. In fact, this new species has been encountered more than once during the investigation and research on cold-seep fauna in the South China Sea, but in these studies, it has been misidentified as Malletia sp. or as Yoldiella not to species (Dong et al. 2021; Ke et al. 2022; He et al. 2023). The detailed description of Y. haimaensis sp. nov. contributes to the understanding of macrobenthos in the Haima Cold Seep.

Yoldiella is a difficult taxon to define because the morphological differences within this genus are mostly subtle and there are many closely related species. This genus was established by Verrill and Bush (1897), and it is usually characterized by ovate or wedge-shaped shells, which always have a slight postero-ventral sinuosity. The internal cartilage is often relatively large and occupies a simple notch. The external ligament is weak, and the pallial sinus is usually indistinct. Since then, its classification and composition have been controversial. Allen and Hannah (1986) and Allen et al. (1995) redefined the genus Yoldiella considering shell shape, hinge morphology, musculature, and the extent and course taken by the hindgut, which have been widely accepted by researchers. The redescription of Yoldiella by Coan and Valentich Scott (2012) limited it to forms having an elongate, amphidetic ligament with an internal section, an obscure or absent escutcheon and lunule, and a small pallial sinus.

It is widely acknowledged that soft-tissue analysis plays a crucial role in contemporary malacology. However, there is a significant lack of detailed morphological descriptions for this taxon. Yonge and Calman (1939) presented a comprehensive summary of the anatomical characteristics of Protobranchia. Purchon (1956) surveyed the structure and function of the stomach throughout the Protobranchia to establish evidence for phylogenetic relationships between the families and orders within the subclass. Xu (1999) highlighted a significant distinction in the anatomical features of Protobranchia compared to other bivalves, emphasizing the muscular nature of their feet for crawling. The ctenidium assumes a lamellar shape primarily for respiratory purposes, while the well-developed labial palp functions as a distinct feeding organ. Okutani and Fujiwara (2005) illustrated the soft part of Yoldiella kaikonis Okutani & Fujiwara, 2005. Reed et al. (2014) described and deliberated on the gonad morphology and oocyte size of four species of Yoldiella. Yonge and Calman (1939) and Allen et al. (1995) provided detailed accounts of the anatomical characteristics of Y. lucida, the type species of Yoldiella, noting that the organs in this species are more compactly arranged due to its relatively shorter shell, which exhibits a greatly abbreviated posterior rostrum. The labial palps are moderately large and extended, with well-developed, long, and muscular palp proboscides. The foot is relatively larger, as are the ctenidia, the filaments of which are broader and deeper. The shell and anatomical features of the new species align with the aforementioned characteristics of Y. lucida.

The genus Yoldiella may encompass species from other genera due to its small size and the potential for confusion with immature specimens. La Perna (2004) observed that, based on the traditional description at the time, Yoldiella was regarded as a provisional “pigeon-hole”, where numerous species were temporarily allocated in a kind of waiting list, rather than a natural group. It was proposed that within this genus, some clusters of morphologically similar species could be recognized, with the morphological differences among these clusters suggesting distinct systematic ranks. Benaim and Absalão (2011a) followed this conjecture by analyzing the morphological characteristics of some Atlantic species using empty shell specimens and proposed three distinct clusters, grouping together Y. nana, Y. inconspicua, and Y. americana Allen, H.L. Sanders & F. Hannah, 1995. However, in our molecular phylogenetic analysis of these three species, they did not cluster in a single branch. The feasibility of this hypothesis requires validation through the combination of anatomy, morphometry, molecular data, and other methods. In addition, hinge plate features are considered a significant diagnostic feature in descriptions of Yoldiella species, especially the width of the posterior hinge plate (Benaim and Absalão 2011b; Benaim et al. 2011). Further meticulous examination of all taxa is essential before achieving confidence in the accuracy of classification.

Based on anatomical and morphological characteristics, Protobranchia has traditionally been regarded as a monophyletic group (Pelseneer 1889; Xu 1999), but the widespread application of molecular methods in phylogenetic analysis has not provided support for the monophyletism of the Protobranchia (Combosch et al. 2017). Recent studies have supported the monophyly of Protobranchia and five superfamilies (Nuculoidea, Manzanelloidea, Solemyoidea, Sareptoidea, and Nuculanoidea) (Smith et al. 2011; Sharma et al. 2012; Bieler et al. 2014; González et al. 2015; Lemer et al. 2019; Sato et al. 2020). Notably, the monophyly of lower taxa (family and below) is still uncertain, particularly within the superfamily Nuculanoidea. In this study, all currently recognized families were found to be non-monophyletic, which is consistent with previous research findings (Sharma et al. 2013; Sato et al. 2020). In addition, the limited support for the majority of nodes within the Nuculanoidea has resulted in largely ambiguous internal relationships. This phenomenon may be the result of a combination of factors, possibly the lack of genetic information, deficient taxon sampling, or both, or even the possibility that this taxon may have undergone rapid radiation (Lemer et al. 2019).

Phylogenetic analysis of the Yoldiidae revealed differences between the evolutionary relationships inferred from the COI and the 18S rRNA genes. Specifically, the phylogenetic tree based on COI gene shows that Yoldia is polyphyletic, in which species were interspersed among other genera rather than forming a distinct clade. In contrast, the phylogenetic tree derived from the 18S rRNA gene shows that Yoldia is a monophyletic group within a single clade. Such incongruities may be attributed to variations in the evolutionary rates of the two genes. The genetic markers employed in this study, such as COI and H3, exhibit rapid evolutionary rates, particularly at the third codon position, which are almost certainly saturated when applied across multiple taxonomic families. In addition, Mallatt et al. (2010) explored the metazoan tree using almost complete rRNA genes (18S and 28S), which suggested that non-monophyly of Mollusca in analyses based on 18S rRNA might be explained by these non-homologous forms, and indeed paraphyly/polyphyly of Mollusca has been found in many studies of metazoan phylogeny based on 18S rRNA. In addition, Mallatt et al. (2010) investigated the metazoan phylogenetic tree by analyzing nearly complete 18S rRNA genes. Their findings did not support the monophyly of Mollusca, indicating that the non-monophyly of Mollusca in 18S rRNA-based analyses could be attributed to non-homologous forms. Indeed, the paraphyly/polyphyly of Mollusca has been observed in numerous studies on metazoan phylogeny utilizing 18S rRNA (e.g. Winnepenninckx et al. 1996; Mallatt et al. 2010, 2012). Based on the findings derived from the analysis of 18S rRNA, the potential scenario cannot be ruled out. Consequently, great care needed to be taken in the application of 18S rRNA for phylogenetics to avoid the incorporation of non-homologous variants of 18S rRNA in the assessments.

There may be more problems with traditional identification that relies solely on anatomy and morphology. Challenges persist due to the difficulty in obtaining protobranch samples, resulting in relatively limited molecular data with insufficient resolution. As shown in our phylogenetic analysis, the bootstrap values are low which may be unreliable. Nonetheless, this work represented a minor step forward and demonstrated the complexity of the taxonomic classification system within this taxon and made it clear that a greater coverage of taxa and more informative genetic markers (not limited to a single gene, but based on high-throughput sequencing technologies for massive orthologous genes from transcriptomes or even genomes) will provide great potential for improving the resolution of phylogenetic classification in this group in the future.

Additional information

Conflict of interest

The authors have declared that no competing interests exist.

Ethical statement

No ethical statement was reported.

Funding

This work was supported by the National Key Research and Development Program of China (2021YFE0193700), Qingdao New Energy Shandong Laboratory Open Project (QNESL OP202306), the Science and Technology Innovation Project of Laoshan Laboratory (LSKJ202203100), the Strategic Priority Research Program of the Chinese Academy of Sciences (XDB42000000), and the Taishan Scholars Program (tsqn202306280).

Author contributions

Junlong Zhang: research conceptualization, funding provision, writing, and editing; Qi Gao: morphological study, molecular experiments, methodology, data analysis and interpretation, original draft writing. Yan Tang: methodology, data analysis and interpretation, writing. All authors contributed critically to the drafts and gave final approval for publication.

Author ORCIDs

Yan Tang https://orcid.org/0000-0003-0564-7279

Junlong Zhang https://orcid.org/0000-0003-3831-4673

Data availability

All of the data that support the findings of this study are available in the main text or Supplementary Information.

References

  • Allen JA (1978) Evolution of the deep sea protobranch bivalves. Philosophical Transactions of the Royal Society of London. Series B, Biological Sciences 284: 387–401. https://doi.org/10.1098/rstb.1978.0076
  • Allen JA, Hannah FJ (1986) A reclassification of the Recent genera of the subclass Protobranchia Mollusca Bivalvia. Journal of Conchology 324: 225–249.
  • Allen JA, Sanders HL, Hannah F (1995) Studies on the deep-sea Protobranchia (Bivalvia); the subfamily Yoldiellinae. Bulletin of the Natural History Museum Zoology Series 61: 11–90.
  • Benaim NP, Absalão RS (2011a) Deep sea Yoldiella (Pelecypoda: Protobranchia: Yoldiidae) from Campos Basin, Rio de Janeiro, Brazil. Journal of the Marine Biological Association of the United Kingdom 91(2): 513–529. https://doi.org/10.1017/S0025315410001372
  • Benaim NP, Absalão RS (2011b) Discriminating among similar deep-sea Yoldiella (Pelecypoda: Protobranchia) species with a morphometric approach. Journal of the Marine Biological Association of the United Kingdom 91(8): 1665–1672. https://doi.org/10.1017/S0025315411000063
  • Benaim NP, Viegas D, Absalão RS (2011) How features of the hinge plate aid in discriminating among three Yoldiella (Pelecypoda, Protobranchia) species from the Campos Basin, Brazil. Zootaxa 2883(1): 39–51. https://doi.org/10.11646/zootaxa.2883.1.3
  • Bieler R, Mikkelsen PM, Collins TM, Glover EA, González VL, Graf DL, Harper EM, Healy J, Kawauchi GY, Sharma PP, Staubach S, Strong EE, Taylor JD, Tëmkin I, Zardus JD, Clark S, Guzmán A, McIntyre E, Sharp P, Giribet G (2014) Investigating the bivalve tree of life—An exemplar-based approach combining molecular and novel morphological characters. Invertebrate Systematics 28(1): 32. https://doi.org/10.1071/IS13010
  • Coan EV, Valentich Scott P (2012) Bivalve Seashells of Tropical West America: Marine Bivalve Mollusks from Baja California to Northern Perú (Part I). Santa Barbara Museum of Natural History, Santa Barbara, 598 pp.
  • Combosch DJ, Collins TM, Glover EA, Graf DL, Harper EM, Healy JM, Kawauchi GY, Lemer S, McIntyre E, Strong EE, Taylor JD, Zardus JD, Mikkelsen PM, Giribet G, Bieler R (2017) A family-level Tree of Life for bivalves based on a Sanger-sequencing approach. Molecular Phylogenetics and Evolution 107: 191–208. https://doi.org/10.1016/j.ympev.2016.11.003
  • Dall WH (1908) Reports on the dredging operations off the west coast of Central America to the Galapagos, to the west coast of Mexico, and in the Gulf of California, in charge of Alexander Agassiz, carried on by the U.S. Fish Commission steamer “Albatross,” during 1891, Lieut.-Commander Z.L. Tanner, U.S.N., commanding. XXXVII. Reports on the scientific results of the expedition to the eastern tropical Pacific, in charge of Alexander Agassiz, by the U.S. Fish Commission steamer “Albatross”, from October, 1904 to March, 1905, Lieut.-Commander L.M. Garrett, U.S.N., commanding. XIV. The Mollusca and Brachiopoda. Bulletin of the Museum of Comparative Zoology 43(6): 205–487, [pls 1–22].
  • Darriba D, Taboada GL, Doallo R, Posada D (2012) jModelTest 2: More models, new heuristics and parallel computing. Nature Methods 9(8): 772–772. https://doi.org/10.1038/nmeth.2109
  • Don RH, Cox PT, Wainwright BJ, Baker K, Mattick JS (1991) ‘Touchdown’ PCR to circumvent spurious priming during gene amplification. Nucleic Acids Research 19(14): 4008–4008. https://doi.org/10.1093/nar/19.14.4008
  • Dong D, Li X, Yang M, Gong L, Li Y, Sui J, Gan Z, Kou Q, Xiao N, Zhang J (2021) Report of epibenthic macrofauna found from Haima cold seeps and adjacent deep-sea habitats, South China Sea. Marine Life Science & Technology 3(1): 1–12. https://doi.org/10.1007/s42995-020-00073-9
  • González VL, Andrade SCS, Bieler R, Collins TM, Dunn CW, Mikkelsen PM, Taylor JD, Giribet G (2015) A phylogenetic backbone for Bivalvia: An RNA-seq approach. Proceedings of the Royal Society B: Biological Sciences 282(1801): 20142332. https://doi.org/10.1098/rspb.2014.2332
  • He X, Xu T, Chen C, Liu X, Li Y-X, Zhong Z, Gu X, Lin Y-T, Lan Y, Yan G, Sun Y, Qiu J-W, Qian P-Y, Sun J (2023) Same (sea) bed different dreams: Biological community structure of the Haima seep reveals distinct biogeographic affinities. The Innovation Geoscience 1(2): 100019. https://doi.org/10.59717/j.xinn-geo.2023.100019
  • International Commission on Zoological Nomenclature (1985) Opinion 1306. Ledella bushae Warén, 1978 is the type species of Ledella Verrill & Bush, 1897 (Mollusca, Bivalvia). Bulletin of Zoological Nomenclature 42: 146–147.
  • Kalyaanamoorthy S, Minh BQ, Wong TKF, Von Haeseler A, Jermiin LS (2017) ModelFinder: Fast model selection for accurate phylogenetic estimates. Nature Methods 14(6): 587–589. https://doi.org/10.1038/nmeth.4285
  • Katoh K, Standley DM (2013) MAFFT Multiple Sequence Alignment Software version 7: Improvements in performance and usability. Molecular Biology and Evolution 30(4): 772–780. https://doi.org/10.1093/molbev/mst010
  • Ke Z, Li R, Chen Y, Chen D, Chen Z, Lian X, Tan Y (2022) A preliminary study of macrofaunal communities and their carbon and nitrogen stable isotopes in the Haima cold seeps, South China Sea. Deep-sea Research. Part I, Oceanographic Research Papers 184: 103774. https://doi.org/10.1016/j.dsr.2022.103774
  • La Perna R (2004) The identity of Yoldia micrometrica Seguenza, 1877 and three new deep-sea protobranchs from the Mediterranean (Bivalvia). Journal of Natural History 38(8): 1045–1057. https://doi.org/10.1080/0022293021000056936
  • Lemer S, Bieler R, Giribet G (2019) Resolving the relationships of clams and cockles: Dense transcriptome sampling drastically improves the bivalve tree of life. Proceedings. Biological Sciences 286(1896): 20182684. https://doi.org/10.1098/rspb.2018.2684
  • Lovén SL (1846) Index Molluscorum litora Scandinaviae occidentalia habitantium. Kongliga Vetenskaps Akademiens Förhandlingar 1846: 134–160, 182–204.
  • Mallatt J, Craig CW, Yoder MJ (2010) Nearly complete rRNA genes assembled from across the metazoan animals: Effects of more taxa, a structure-based alignment, and paired-sites evolutionary models on phylogeny reconstruction. Molecular Phylogenetics and Evolution 55(1): 1–17. https://doi.org/10.1016/j.ympev.2009.09.028
  • Mallatt J, Craig CW, Yoder MJ (2012) Nearly complete rRNA genes from 371 Animalia: Updated structure-based alignment and detailed phylogenetic analysis. Molecular Phylogenetics and Evolution 64(3): 603–617. https://doi.org/10.1016/j.ympev.2012.05.016
  • Neulinger S, Sahling H, Süling J, Imhoff J (2006) Presence of two phylogenetically distinct groups in the deep-sea mussel Acharax (Mollusca: Bivalvia: Solemyidae). Marine Ecology Progress Series 312: 161–168. https://doi.org/10.3354/meps312161
  • Nguyen L-T, Schmidt HA, Von Haeseler A, Minh BQ (2015) IQ-TREE: A fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Molecular Biology and Evolution 32(1): 268–274. https://doi.org/10.1093/molbev/msu300
  • Okutani T, Fujiwara Y (2005) Four protobranch bivalves collected by the ROV Kaiko from hadal depths in the Japan Trench. Venus (Tokyo) 63(3–4): 87–94.
  • Pelseneer P (1889) Sur la classification phylogenetique de pelecypods. Bulletin Scientifique de la France et de la Belgique (série 3) 20 (1–4): 27–52.
  • Reed AJ, Morris JP, Linse K, Thatje S (2013) Plasticity in shell morphology and growth among deep-sea protobranch bivalves of the genus Yoldiella (Yoldiidae) from contrasting Southern Ocean regions. Deep-sea Research. Part I, Oceanographic Research Papers 81: 14–24. https://doi.org/10.1016/j.dsr.2013.07.006
  • Reed AJ, Morris JP, Linse K, Thatje S (2014) Reproductive morphology of the deep-sea protobranch bivalves Yoldiella ecaudata, Yoldiella sabrina, and Yoldiella valettei (Yoldiidae) from the Southern Ocean. Polar Biology 37(10): 1383–1392. https://doi.org/10.1007/s00300-014-1528-4
  • Sato K, Kano Y, Setiamarga DHE, Watanabe HK, Sasaki T (2020) Molecular phylogeny of protobranch bivalves and systematic implications of their shell microstructure. Zoologica Scripta 49(4): 458–472. https://doi.org/10.1111/zsc.12419
  • Sharma PP, González VL, Kawauchi GY, Andrade SCS, Guzmán A, Collins TM, Glover EA, Harper EM, Healy JM, Mikkelsen PM, Taylor JD, Bieler R, Giribet G (2012) Phylogenetic analysis of four nuclear protein-encoding genes largely corroborates the traditional classification of Bivalvia (Mollusca). Molecular Phylogenetics and Evolution 65(1): 64–74. https://doi.org/10.1016/j.ympev.2012.05.025
  • Sharma PP, Zardus JD, Boyle EE, González VL, Jennings RM, McIntyre E, Wheeler WC, Etter RJ, Giribet G (2013) Into the deep: A phylogenetic approach to the bivalve subclass Protobranchia. Molecular Phylogenetics and Evolution 69(1): 188–204. https://doi.org/10.1016/j.ympev.2013.05.018
  • Smith SA, Wilson NG, Goetz FE, Feehery C, Andrade SCS, Rouse GW, Giribet G, Dunn CW (2011) Resolving the evolutionary relationships of molluscs with phylogenomic tools. Nature 480(7377): 364–367. https://doi.org/10.1038/nature10526
  • Vaidya G, Lohman DJ, Meier R (2011) SequenceMatrix: Concatenation software for the fast assembly of multi-gene datasets with character set and codon information. Cladistics 27(2): 171–180. https://doi.org/10.1111/j.1096-0031.2010.00329.x
  • Verrill AE, Bush KJ (1897) Revision of the genera of Ledidae and Nuculidae of the Atlantic Coast of the United States. American Journal of Science 3(13): 51–63. https://doi.org/10.2475/ajs.s4-3.13.51
  • Wang H, Liu H, Wang X, Zhang J, Sirenko BI, Liu C, Dong D, Li X (2022) Stirring the deep, disentangling the complexity: Report on the third species of Thermochiton (Mollusca: Polyplacophora) From Haima Cold Seeps. Frontiers in Marine Science 9: 889022. https://doi.org/10.3389/fmars.2022.889022
  • Xu F (1999) Fauna Sinica. Phylum Mollusca. Class Bivalvia. Subclasses Protobranchia and Anomalodesmata. Science Press, Beijing, 244 pp.
  • Yao G, Zhang H, Xiong P, Jia H, Shi Y, He M (2022) Community characteristics and genetic diversity of macrobenthos in Haima Cold Seep. Frontiers in Marine Science 9: 920327. https://doi.org/10.3389/fmars.2022.920327
  • Yonge M, Calman WT (1939) The protobranchiate mollusca; a functional interpretation of their structure and evolution. Philosophical Transactions of the Royal Society of London. Series B, Biological Sciences 230(566): 79–148. https://doi.org/10.1098/rstb.1939.0005
  • Zardus JD, Etter RJ, Chase MR, Rex MA, Boyle EE (2006) Bathymetric and geographic population structure in the pan‐Atlantic deep‐sea bivalve Deminucula atacellana (Schenck, 1939). Molecular Ecology 15(3): 639–651. https://doi.org/10.1111/j.1365-294X.2005.02832.x
  • Zhao J, Liang Q, Wei J, Tao J, Yang S, Liang J, Lu J, Wang J, Fang Y, Gong Y, He X (2020) Seafloor geology and geochemistry characteristic of methane seepage of the “Haima” cold seep, northwestern slope of the South China Sea. Geochimica 49: 108–118.

Supplementary materials

Supplementary material 1 

Primers used in this study

Qi Gao , Yan Tang , Junlong Zhang

Data type: xlsx

This dataset is made available under the Open Database License (http://opendatacommons.org/licenses/odbl/1.0/). The Open Database License (ODbL) is a license agreement intended to allow users to freely share, modify, and use this Dataset while maintaining this same freedom for others, provided that the original source and author(s) are credited.
Download file (11.89 kb)
Supplementary material 2 

The sequences of COI, 18S and H3 genes used for phylogenetic analyses of the superfamily Nuculanoidea

Qi Gao , Yan Tang , Junlong Zhang

Data type: xlsx

This dataset is made available under the Open Database License (http://opendatacommons.org/licenses/odbl/1.0/). The Open Database License (ODbL) is a license agreement intended to allow users to freely share, modify, and use this Dataset while maintaining this same freedom for others, provided that the original source and author(s) are credited.
Download file (13.91 kb)
Supplementary material 3 

The sequences of COI and 18S genes used for phylogenetic analyses of the family Yoldiidae

Qi Gao , Yan Tang , Junlong Zhang

Data type: xlsx

This dataset is made available under the Open Database License (http://opendatacommons.org/licenses/odbl/1.0/). The Open Database License (ODbL) is a license agreement intended to allow users to freely share, modify, and use this Dataset while maintaining this same freedom for others, provided that the original source and author(s) are credited.
Download file (10.61 kb)
Supplementary material 4 

Phylogenetic tree inferred by Bayesian inference (BI) based on combined gene dataset (COI+18S+H3)

Qi Gao , Yan Tang , Junlong Zhang

Data type: pdf

This dataset is made available under the Open Database License (http://opendatacommons.org/licenses/odbl/1.0/). The Open Database License (ODbL) is a license agreement intended to allow users to freely share, modify, and use this Dataset while maintaining this same freedom for others, provided that the original source and author(s) are credited.
Download file (209.59 kb)
Supplementary material 5 

Phylogenetic tree inferred by Maximum likelihood (ML) based on combined gene dataset (COI+18S +H3)

Qi Gao , Yan Tang , Junlong Zhang

Data type: pdf

This dataset is made available under the Open Database License (http://opendatacommons.org/licenses/odbl/1.0/). The Open Database License (ODbL) is a license agreement intended to allow users to freely share, modify, and use this Dataset while maintaining this same freedom for others, provided that the original source and author(s) are credited.
Download file (206.46 kb)
login to comment