Research Article |
Corresponding author: Donald J. Colgan ( don.colgan@austmus.gov.au ) Academic editor: Nathalie Yonow
© 2021 Donald J. Colgan, Shane T. Ahyong, Karine Mardon, Ian M. Brereton.
This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Citation:
Colgan DJ, Ahyong ST, Mardon K, Brereton IM (2021) Rare specimen identification in an un-integrated taxonomy: implications of DNA sequences from a Taiwanese Philine (Mollusca, Philinidae). ZooKeys 1060: 93-110. https://doi.org/10.3897/zookeys.1060.28809
|
Many species of the gastropod genus Philine have been named from northeastern Asia but scanty descriptions based predominantly on shells make it difficult to determine which are valid. This, plus the sporadic anatomical and genetic information available for many of these species has led to what may be described as an un-integrated taxonomy. In this situation, it is generally preferable to postpone dissection of rare and unusual specimens until relevant diagnostic characters can be established in broader studies. Micro-CT scanning and DNA sequencing were used to examine such a specimen collected recently from deep waters off northeastern Taiwan. Micro-CT examination of the morphology of the internal shell and gizzard plates suggested that, among named species, the sequenced specimen is most similar to P. otukai. It cannot, however, be definitively referred to P. otukai as that species lacks adequate anatomical description or known DNA sequences. Phylogenetic analyses of newly collected DNA sequences show the specimen to be most closely related to, but distinct from the northern Atlantic Ocean and Mediterranean species, Philine quadripartita. The sequences also confirm genetically that five or more species of Philine occur in northeast Asia, including at least three subject to considerable taxonomic uncertainty.
Gizzard plates, 16S ribosomal RNA, micro-CT scanning, Scaphopoda
Accurate delimitation of species of the highly speciose cephalaspidean genus Philine Ascanius, 1772 is hampered by the brief original descriptions of many of its named taxa. Many such descriptions were made using only the shell, which has limited diagnostic value, and often these have not been subsequently supplemented by studies of internal anatomy, including morphology of the diagnostically important gizzard plates, or DNA sequences.
The difficulty of delimitation is especially marked in northeastern Asia, as exemplified by the situation in Taiwan from which three named Philine species (P. argentata Gould, 1859; P. vitrea Gould, 1859, and P. otukai Habe, 1946) are currently recognised (
Philine species from northeastern Asia with medium- or large-sized shells. Localities of identified type specimens are followed by (T). Information categories are abbreviated as: D, DNA sequence; G, gizzard plate morphology, O, other anatomical.
Species | Locality | Information | Status |
---|---|---|---|
P. acutangula Adams, 1862 | Gulf of Lian-Tung, Hulu-Shan Bay, China [Gulf of Liadong, Hulushan Bay] | No figure or type description | Species inquirendum |
P. argentata Gould, 1859 | Hokkaido Bay, Japan (T) | G, O: |
Synonym of P. orientalis ( |
D: |
|||
P. coreanica Adams, 1855 | Coreen archipelago on mud flats [Korea] |
|
Species inquirendum: “P. aperta” fide |
P. crenata Adams, 1862 | Tsu-Shima, Korea Strait, 46 m | Species inquirendum: “P. aperta” fide |
|
P. habei Valdés, 2008 | Fiji, 17°05'S, 178°55'W, 654–656 m (T) | D: |
Accepted |
G, O: |
|||
P. japonica Lischke, 1872 | Jedo (Tokyo Bay), Japan (T) | D: |
P. orientalis synonym ( |
O: Larvae: |
|||
P. kinglipini Tchang, 1934 | Bays of Tsangkou and Nukukou | G, O: |
Accepted. No type designated. |
P. kurodai Habe, 1946 | Wakayama Prefecture, Japan (T) | Accepted by |
|
P. miyadii Habe, 1946 | Synonym of P. scalpta fide |
||
P. orientalis Adams, 1854 | “Eastern seas” (T) | G, O: |
Accepted |
D: |
|||
P. otukai Habe, 1946 | Penghu Islands, Taiwan (T) | Accepted | |
P. paucipapillata Price, Gosliner & Valdés, 2011 | Kampote and Prek Romeas, Cambodia (T) | G, O: |
|
D: |
|||
P. scalpta A. Adams, 1862 | Tsu-Shima and Korea Strait 50–80 m | D: |
Synonym of P. vitrea ( |
P. striatella Tapparone-Canefri,1874 | Yokohama, Japan (T) | O (radula): |
Synonym of P. japonica fide |
P. vitrea Gould, 1859 | Hong Kong (T) | Accepted |
We used non-destructive morphological analysis of the collected specimen and DNA sequencing of a small external tissue sample to make initial comparisons with named species. For morphological analysis, we used micro-CT scanning (
We determined sequences of 16S ribosomal RNA (16S rRNA), histone H3 (H3), and the D1 expansion region of 28S ribosomal RNA (D1 28S rRNA). We then performed phylogenetic analyses of the sequences, including sequences from previous DNA studies of the Philinidae (
AMS C.559479, off Dasi, Yilan County, north-east Taiwan, 24°54'N, 122°E, 5 xi 2015, coll. S. Ahyong, commercial bottom trawl, 300–400 m depth, fixed and stored in 80% ethanol. The specimen was collected together with other deep-water species including the crustaceans, Bathynomus doederleini Ortmann, 1894, Homola orientalis Henderson, 1888, and Metanephrops thomsoni (Bate, 1888).
Digital callipers were used for external measurements, made after preservation. Measurements of the relative sizes of internal structures, such as gizzard plates, were determined from images using the beta 4.0.3 version of Scion Image (Scion Corporation, Frederick, MD, USA). Sizes of internal structures were estimated by multiplying the ratio of the relative lengths of the relevant structure and the longest shell axis on an image by the length of this axis determined by callipers on the actual specimen.
Micro-CT scanning was performed directly on the specimen after removal from ethanol storage with a Siemens Inveon micro-CT scanner operated at 80 KV energy, 250 µA intensity with 540 projections per 360° and 2200 ms exposure time. The sample was scanned at a nominal isotropic resolution of 27.8 µm. The data were reconstructed using a Feldkamp conebeam back-projection algorithm provided by an Inveon Acquisition Workstation from Siemens.
DNA was extracted from a small amount of mantle tissue using the Bioline Isolate II Kit following the manufacturer’s instructions. The final centrifugation was performed with 100 µL of the kit’s elution buffer. PCR amplifications followed the methods of
16S rRNA: 16Sar CGCCTGTTTATCAAAAACAT (
28S rRNA: 28S D1F ACCCSCTGAAYTTAAGCAT (
H3: H3NF ATGGCTCGTACCAAGCAGAC (
H3NR ATRTCCTTGGGCATGATTGTTAC (
PCR products were checked by electrophoresis on 2% agarose gels run with TAE buffer, including Gel Red (Biotium Inc., Fremont, CA) and visualised with UV-fluorescence. Single-banded products of the appropriate size were purified by Exo-SAP-ITTM (Thermo-Fisher Scientific) and sequenced commercially at Macrogen (Seoul, Korea) with the primers used in the original amplification.
GenBank Accession Numbers for new sequences: MH340050 (D1 28S rRNA), MH340051 (16S rRNA), MH340052 (histone H3).
Phylogenetic analyses were conducted using all Philine sequences for each gene available on 20 September 2020. Sequences from Philinopsis Pease, 1860 were used to root analyses for the 16S rRNA dataset. The sets of sequences for each gene were aligned with ClustalW (
Maximum likelihood analyses (ML) analyses were conducted at the CIPRES Science Gateway (https://www.phylo.org/portal2/home.action,
Animal, external: (Fig.
Shell: (Fig.
Gizzard plates: (Figs
The prey of Philine are predominantly bivalves (
The 16S rRNA alignment comprised 444 positions of which 38 were variable but not parsimony-informative and 169 which were parsimony-informative. The optimal tree found in the ML analysis of the alignment had a ln likelihood of –3619.22 and there were 650 bootstrap replicates conducted under the MRE criterion. After Gblocks filtering, there were 370 positions in the alignment, of which 30 were variable but not parsimony-informative and 133 were parsimony-informative. ML analysis of these data had a ln likelihood of –2831.37 and there were 800 bootstrap replicates conducted under the MRE criterion. The optimal tree found in the ML analysis of the reduced 16S rRNA dataset (Fig.
Phylogenetic relationships of C.559479 based on maximum likelihood analysis of the reduced 16S rRNA dataset. More distant outgroups have been removed and the topology rooted on Philine scabra + P. indistincta. Numbers near nodes refer to bootstrap percentages above 70%. The scale bar indicates 0.05 changes per site. Sequences are identified by accession number and species name or informal designation recorded in GenBank except those labelled P. quadripartita for which the species names have been changed for reasons outlined in the text. Accessions with an sp. number designation followed by a space and “TO” with a one or two digit designation refer to sequences from
In all analyses of the 16S rRNA data, MN326885 collected from P. argentata by
The 28S rRNA D1 alignment (Suppl. material
The histone H3 alignment (Suppl. material
Morphologically, C.559479 has more apparent affinity with P. otukai Habe, 1946, than other Taiwanese species. If the taxa were conspecific, this would represent a very large increase in the depth range of P. otukai which is reported to extend only to 100 metres (
Specimen C.559479 is readily distinguished from P. argentata, which has two slit-like recesses on the outer surfaces of its gizzard plates (
Specimen C.559479 differs from P. vitrea in the small size of the latter (shell dimensions of the type specimen 10 × 8 × 3 mm: Gould 1859). It also differs in the apparently much lower projection of the upper end of the outer lip of P. vitrea, allowing for the damaged condition of this shell section in the specimen studied here.
Sufficient information is available to distinguish C.559479 morphologically from only a few other regional species of Philine. These include P. scalpta in which the gizzard plates are quasi-trapezoidal with elongated pores (
Philine paucipapillata Price, Gosliner & Valdés, 2011 was described as a new species from the South China Sea. Specimen C.559479 is readily distinguished morphologically from this species in coloration and size (P. paucipapillata is uniformly white and ~ 4–5 cm in length). The pores in the gizzard plates of P. paucipapillata are minute (
Philine kurodai Habe, 1946, described from Japan, has been reported from Taiwan although the Taiwanese records are now regarded as being based on specimens of P. otukai (
Three species from northeastern Asia, P. acutangula Adams, 1862, P. crenata Adams, 1862, and P. coreanica Adams, 1855, have received little recent attention, particularly after
Analyses including the DNA sequences of C.559479 collected here highlight the taxonomic complexity of Philine in northeastern Asia. The available 16S rRNA data from the region represent at least five species. One of these species is represented by the C.559479 sequence from Taiwan. The other species are P. paucipapillata, P. argentata sensu
The three sequences of 16S rRNA that were most similar to C.559479 are recorded in GenBank as belonging to Philine aperta. However, the true identity of the specimens from which these sequences were derived appears to be P. quadripartita, a species from the northern Atlantic Ocean and Mediterranean Sea. The sequences form a distinct lineage separate from the two South African specimens from near the type locality of true P. aperta, and were determined from specimens taken within the range of P. quadripartita: AY345016 derives from a specimen from south-eastern Spain (
The “P. orientalis” 16S rRNA sequence accession JQ691684) from Japan (
The difficulty of identifying the unusual deep water Philine specimen studied here emphasises that destructive sampling should be minimised where the taxonomy is un-integrated. We tentatively suggest from examination of the external morphology and micro-CT scanning that it may have affinities with the shallow water P. otukai but any certainty is precluded by the lack of a description of that species, especially one detailing the range of variation it might encompass in characters such as the size of the gizzard plate pores. There is no diagnosis of P. otukai that would guide specimen dissection for relevant anatomical characters. Considerable further research efforts are required to provide a framework in which such characters could be sought. These include DNA studies of a wider range of regional species of Philine, particularly sequences definitely from P. orientalis, P. kurodai, P. otukai, and P. vitrea. Detailed morphological investigations are needed for P. japonica, P. otukai, P. kurodai, and P. vitrea. The status of species such as P. coreanica may remain unknown but it may possibly be resolved if shell fragments from type material ever become available for macromolecular analysis, especially proteomics.
We thank the Australian Museum for funding and Cameron Slatyer for logistical support. We are grateful to Dr W. Rudman, Dr E. Chaban, Dr T. Gosliner, and anonymous reviewers for comments on earlier versions of the manuscript and to the editor, Dr N. Yonow, for very helpful advice. STA gratefully acknowledges support and hospitality of Prof. T.-Y. Chan in Taiwan in 2015. The authors acknowledge the facilities and the scientific and technical assistance of the National Imaging Facility at the Centre for Advanced Imaging, University of Queensland. We thank Professor K.K. Chan for providing the copy of a reference.
Reconstructed micro-CT scan of C.559479
Data type: Mpg file
Explanation note: Reconstructed micro-CT scan of C.559479.
Figure S1. Image from the reconstructed micro-CT scan of C.559479
Data type: Pdf file
Explanation note: Image from the reconstructed micro-CT scan of C.559479. Dorsal view of the shell of C.559479, with boxed area illustrating irregularly cancellate sculpture.
Figure S2. Phylogenetic relationships of C.559479 based on maximum likelihood analysis of the D1 28S rRNA dataset
Data type: Pdf file
Explanation note: Phylogenetic relationships of C.559479 based on maximum likelihood analysis of the D1 28S rRNA dataset. Numbers near nodes refer to bootstrap percentages above 70%. Scale bar indicates 0.01 changes per site. Sequences are identified by accession number and species name or informal designation recorded in GenBank. Accessions with an sp. number designation followed by a space and “TO” with a one or two digit designation refer to sequences from
Figure S3. Phylogenetic relationships of C.559479 based on Maximum Likelihood analysis of the histone H3 dataset.
Data type: Pdf file
Explanation note: Phylogenetic relationships of C.559479 based on Maximum Likelihood analysis of the histone H3 dataset. Numbers near nodes refer to bootstrap percentages above 70%. The scale bar indicates 0.05 changes per site. Sequences are identified by accession number and species name or informal designation recorded in GenBank. Accessions with a sp. number designation followed by a space and “TO” with a one or two digit designation refer to sequences from