Intraspecific Variation in Mitogenomes of Five Crassostrea Species Provides Insight into Oyster Diversification and Speciation

Morphological and Molecular Analysis Identified a Subspecies of Crassostrea ariakensis (Fujita, 1913) along the Coast of Asia

Crassostrea ariakensis (Fujita, 1913) is one of the most important economic and ecological oysters that is naturally distributed along the coast of Asia, separated by the Yangtze River estuary. They are usually compared as different populations, while there is no consensus on whether C. ariakensis in northern and southern areas should be considered as two species or subspecies. Here, we analyzed morphological characteristics, COI, 16s rRNA, mitogenome sequences, and species delimitation analysis (ASAP and PTP) to resolve the intraspecific taxonomic status of the C. ariakensis. Phylogenetic and ASAP analysis highlight that C. ariakensis was divided into N-type and S-type. PTP was unable to differentiate between the two types of C. ariakensis. The divergence time of N-type and S-type C. ariakinsis is estimated to be 1.6 Mya, using the relaxed uncorrelated lognormal clock method. Additionally, significant morphological differences exist between the two groups in terms of the adductor muscle scar color. Despite these differences, the COI (0.6%) and 16S rRNA (0.6%) genetic distance differences between N-type and S-type C. ariakensis has not yet reached the interspecific level. These results suggest that N-type and S-type C. ariakensis should be treated as different subspecies and renamed as C. ariakensis ariakensis subsp. nov and C. ariakensis meridioyangtzensis subsp. nov.

Effects of salinity on pre- and post-fertilization developmental events in the clam Anomalocardia flexuosa (Linnaeus, 1767)

The knowledge about the effect of salinity on the physiological mechanism of bivalve reproduction is fundamental to improve production strategies in hatcheries. The present work evaluated the influence of different salinity concentrations (15, 20, 25, 30, 35 and 40 g⋅L^-1) on pre- and post-fertilization development processes in the clam, Anomalocardia flexuosa, oocytes obtained by stripping. Salinity directly interfered with the germinal vesicle breakdown (GVBD) rate and in the cellular stability of unfertilized oocytes. Salinity concentrations between 30 and 35 g⋅L^-1 provided better percentages of stable GVBD within 120 min, and incubation of oocytes in the salinity range of 30-35 g⋅L^-1 for a time interval of 80-120 min provided > 80% GVBD. In the post-fertilization analysis, salinity affected the rate of the extrusion of the first and second polar bodies (PB1 and PB2). The release of 50% of the PBs was faster at a salinity of 35 g⋅L^-1, with an estimated time of 10 min for PB1 and 30 min for PB2. Thus, chromosome manipulation methodologies aiming triploids should be applied at 35 g⋅L^-1 salinity, with application of post-fertilization shock before 10 min for PB1 retention or before 30 min for PB2 retention.

Microsatellites as Molecular Markers with Applications in Exploitation and Conservation of Aquatic Animal Populations

A large number of species and taxa has been studied for genetic polymorphism. Microsatellites have been known as hypervariable neutral molecular markers with the highest resolution power in comparison with any other markers. However, the discovery of a new type of molecular marker—single nucleotide polymorphism (SNP) has put the existing applications of microsatellites to the test. To ensure good resolution power in studies of populations and individuals, a number of microsatellite loci from 14 to 20 was often used, which corresponds to about 200 independent alleles. Recently, these numbers have tended to be increased by the application of genomic sequencing of expressed sequence tags (ESTs), and the choice of the most informative loci for genotyping depends on the aims of research. Examples of successful applications of microsatellite molecular markers in aquaculture, fisheries, and conservation genetics in comparison with SNPs have been summarized in this review. Microsatellites can be considered superior markers in such topics as kinship and parentage analysis in cultured and natural populations, the assessment of gynogenesis, androgenesis and ploidization. Microsatellites can be coupled with SNPs for mapping QTL. Microsatellites will continue to be used in research on genetic diversity in cultured stocks, and also in natural populations as an economically advantageous genotyping technique.

A molecular phylogeny of the Persian Gulf and the Gulf of Oman oyster species

The taxonomy of oysters along the northern coasts of the Persian Gulf and the Gulf of Oman is not well recognized. We present a phylogenetic analysis of oyster species in these regions. We combined morphological and molecular techniques to obtain the identity of oysters to the lowest taxonomic levels. Analysis of partial nucleotide sequences from mitochondrial cytochrome c oxidase subunit I (COI) was used for the phylogenetic evaluation. Based on our findings, Iranian samples nested within the genus Saccostrea and belonged to Saccostrea mordax and Saccostrea palmula clades. The shell morphology of the studied samples was variable, as in other rock oyster species. The examination of morphological features was in line with the molecular outcomes, but despite some similarities, Iranian S. palmula had well-developed and elongated chomata. The results also showed that S. mordax and S. palmula possessed significant relative abundance as dominant oysters in the Persian Gulf and the Gulf of Oman, respectively. Phylogenetic analysis revealed that Iranian samples of S. palmula formed a separate subclade from the Gulf of California and Panama samples, with large genetic distances (6–7%). Iranian specimens differed morphologically and genetically, suggesting that they could be a new species, although more research is needed.

Cryptic Marine Diversity in the Northern Arabian Gulf: An Integrative Approach Uncovers a New Species of Oyster (Bivalvia: Ostreidae), Ostrea oleomargarita

Animal biodiversity is greatly underestimated in nontemperate marine regions, especially for intertidal benthic organisms such as oysters. Recent surveys in the northern Arabian Gulf suggest the presence of numerous unidentified species, some of which form shallow reef ecosystems while others are cryptic and found under rocks. In this study, we focused on small oysters from Kuwait, which show typical characteristics in common with the genus Ostrea except for the presence of lophine chomata that would link it to the genera Lopha, Dendostrea, and Alectryonella. Phylogenetic analyses based on mitochondrial and nuclear markers unambiguously placed the Kuwait oyster within the Ostreinae as a sister to the Japanese species Ostrea futamiensis. The hypothesis that the Kuwait oyster represents a new species was assessed with phylogenetic and species delimitation methods combined with a morphological assessment. Results corroborated the Kuwait oysters as a new species herein described as Ostrea oleomargarita Oliver, Salvi, and Al-Kandari, sp. nov. The phylogeny of the Ostreinae shows extensive disagreement between morphology-based genera and phylogenetic clades. The genus Ostrea is polyphyletic, and the form and distribution of taxonomic characters such as chomata are not as definitive as suggested in previous studies. This study, along with other recent investigations, confirmed the Arabian Gulf as a key region for discovering marine animal diversity and suggested a possible biogeographic divide between the Eastern and Western Indo-Pacific. A pattern that has been documented in a growing number of taxa and that warrants further research attention.

Population Genomics, Transcriptional Response to Heat Shock, and Gut Microbiota of the Hong Kong Oyster Magallana hongkongensis

The Hong Kong oyster Magallana hongkongensis, previously known as Crassostrea hongkongensis, is a true oyster species native to the estuarine-coast of the Pearl River Delta in southern China. The species—with scientific, ecological, cultural, and nutritional importance—has been farmed for hundreds of years. However, there is only limited information on its genetics, stress adaptation mechanisms, and gut microbiota, restricting the sustainable production and use of oyster resources. Here, we present population structure analysis on M. hongkongensis oysters collected from Deep Bay and Lantau Island in Hong Kong, as well as transcriptome analysis on heat shock responses and the gut microbiota profile of M. hongkongensis oysters collected from Deep Bay. Single nucleotide polymorphisms (SNPs), including those on the homeobox genes and heat shock protein genes, were revealed by the whole genome resequencing. Transcriptomes of oysters incubated at 25 °C and 32 °C for 24 h were sequenced which revealed the heat-induced regulation of heat shock protein pathway genes. Furthermore, the gut microbe community was detected by 16S rRNA sequencing which identified Cyanobacteria, Proteobacteria and Spirochaetes as the most abundant phyla. This study reveals the molecular basis for the adaptation of the oyster M. hongkongensis to environmental conditions.

The complete mitochondrial genome of Crassostrea hongkongensis from East China Sea indicates species' range may extend northward

BACKGROUND

Crassostrea hongkongensis is an important mariculture shellfish with a relatively narrow distribution range. Recently, larger wild oysters were identified as C. hongkongensis from Sanmen bay in East China Sea. No natural distribution had been reported for this species here, and its origin remains unknown.

METHODS AND RESULTS

We assembled the complete 18,617 bp circular mitochondrial genome of C. hongkongensis from Sanmen bay by next generation sequencing. It included 12 protein-coding genes, 23 tRNAs, and two rRNAs. The A/T content of the mitogenome was higher than its G/C content. Similar values and features were previously found for five other specimens of C. hongkongensis, and were comparable to those of other congeneric species. A phylogenetic analysis based on the 12 protein-coding genes and complete mitochondrial sequence indicated that the six specimens of C. hongkongensis formed a monophyletic group and shared a sister group relationship with C. ariakensis, C. nippona, C. sikamea, C. angulata, C. gigas, and C. iredalei, whereas specimens from the Sanmen bay area clustered later with the five other C. hongkongensis individuals, sharing a sub-clade. The newly sequenced mitogenome had more singleton sites than previously published C. hongkongensis mitogenomes.

CONCLUSIONS

Crassostrea hongkongensis may be a native species, and the species' range extends further to the north than previously known. Our data may therefore contribute to a better understanding of the species diversity and conservation of Crassostrea oysters.

Chromosome-level genome and population genomic analysis provide insights into the evolution and environmental adaptation of Jinjiang oyster Crassostrea ariakensis

The Jinjiang oyster Crassostrea ariakensis, naturally distributing in estuarine regions with low salinity, is an important economic and ecological species in China. However, studies on its genomics and population genetics remain lacking. Here, we assembled the chromosome-level genome of a female C. ariakensis and re-sequenced 261 individuals from five locations in China representing three typical habitats. The C. ariakensis genome was 662.9 Mb with contig N50 length of 5.9 Mb using PacBio HiFi-CCS long reads, and 99.83% sequences were anchored onto 10 pseudochromosomes using Hi-C data. A total of 26,354 protein-coding genes were predicted. We identified three significantly expanded gene families which are closely associated with osmotic pressure regulation, including CDO, SLC13 and SDR. Population structure analysis revealed that the C. ariakensis from five locations were clustered into three typical groups (northern, southern and Shanghai) (K = 3) and their phylogenetic relationship was consistently correlated to their geographical distribution. Furtherly, the differentiation between northern and southern groups was clearly demonstrated by estimated population differentiation coefficient (F_ST  = 0.1154), and the PSMC distribution showed the two groups of effective population size separated at 0.1 Ma. Meanwhile gene flow from southern to Shanghai was detected. Selective sweep analysis between northern and southern group detected genes associated with heat response and salinity adaptation. This study could provide valuable genomic resources and information for further research on the molecular evolution, genetic breeding, biological function and evolutionary adaptation of C. ariakensis.

The queen conch mitogenome: intra- and interspecific mitogenomic variability in Strombidae and phylogenetic considerations within the Hypsogastropoda

Aliger gigas is an economically important and vulnerable marine species. We present a new mitogenome of A. gigas from the Mexican Caribbean and use the eight publicly available Strombidae mitogenomes to analyze intra- and interspecific variation. We present the most complete phylogenomic understanding of Hypsogastropoda to date (17 superfamilies, 39 families, 85 genera, 109 species) to revisit the phylogenetic position of the Stromboidea and evaluate divergence times throughout the phylogeny. The A. gigas mitogenome comprises 15,460 bp including 13 PCGs, 22 tRNAs, and two rRNAs. Nucleotide diversity suggested divergence between the Mexican and Colombian lineages of A. gigas. Interspecific divergence showed high differentiation among Strombidae species and demonstrated a close relationship between A. gigas and Strombus pugilis, between Lambis lambis and Harpago chiragra, and among Tridentarius dentatus/Laevistrombus canarium/Ministrombus variabilis. At the intraspecific level, the gene showing the highest differentiation is ATP8 and the lowest is NAD4L, whereas at the interspecific level the NAD genes show the highest variation and the COX genes the lowest. Phylogenomic analyses confirm that Stromboidea belongs in the non-Latrogastropoda clade and includes Xenophoridea. The phylogenomic position of other superfamilies, including those of previously uncertain affiliation, is also discussed. Finally, our data indicated that Stromboidea diverged into two principal clades in the early Cretaceous while Strombidae diversified in the Paleocene, and lineage diversification within A. gigas took place in the Pleistocene.

Cloning and expression of the ChGstα and ChGstκ genes in the gills of Crassostrea hongkongensis under nanoparticulate and ionic Zn stress

Nanoparticulate and ionic Zn have potential impacts on the detoxification systems of organisms, and Gst genes play key roles in the detoxification of xenobiotics. In this study, we cloned the ChGstα and ChGstκ genes of C. hongkongensis, and studied their expression in gills under nanoparticulate and ionic Zn stress. The results showed that the coding sequences of the ChGstα and ChGstκ genes were 684 and 675 bp, respectively, and had no signal peptide; ChGstα was cytoplasmic, while ChGstκ was mitochondrial. The two genes were expressed in all 8 tested samples, with the most abundant expression observed in hemocytes for ChGstα and digestive glands for ChGstκ. After ZnCl₂ or ZnoNP challenge, the expression of ChGstα decreased significantly in the ZnCl₂ groups, and its expression was higher in the ZnoNP groups than in the ZnCl₂ groups. The expression of ChGstκ was significantly decreased in the ZnCl₂ and ZnoNP groups, and its expression was higher in the ZnoNP groups than in the ZnCl₂ groups except at 3 h post metal Zn stress, which suggested that ChGstα and ChGstκ were more sensitive to ZnoNP than ZnCl₂.

Molecular Species Delimitation of the Genus Reishia (Mollusca: Gastropoda) along the Coasts of China and Korea

Species of the predatory gastropod genus Reishia Kuroda and Habe, 1971 (Muricidae) inhabit intertidal rocky shores in East Asia. Due to their highly variable external shell morphology, the taxonomy of this genus at species-level is still in need of re-evaluation. Using DNA-based delimitation methods, we aimed to ascertain the number of species of Reishia along the coasts of China and adjacent Asian areas. Also, we looked for diagnostic traits using morphology-based statistical approaches. Our genetic data suggest that the studied individuals comprised two separate species of a Reishia complex in this region, in contrast to the previously proposed four or more taxa. This conclusion is further supported by statistical analyses of shell morphological characteristics. The morphospecies R. bronni (Dunker, 1860), R. jubilaea (Tan and Sigurdsson, 1990), and R. luteostoma (Holten, 1803) were assigned to a single taxon, indicating that they might be synonyms of the same species. The morphospecies R. clavigera (Küster, 1860) singly formed one group, suggesting that it is likely a valid name. The estimated divergence time of the two identified taxa indicates that speciation might have been associated with the sea level and temperature fluctuations during the Plio-Pleistocene period. Our study on Reishia species provides crucial information for further research on the ecology, evolutionary biology, and conservation of this genus.

Population genetics and comparative mitogenomic analyses reveal cryptic diversity of Amphioctopus neglectus (Cephalopoda: Octopodidae)

This study presented 96 cox1 and 76 cox3 genes of Amphioctopus neglectus populations. Three distinct lineages were formed from phylogenetic trees and networks constructed using haplotypes. Mitogenomes of A. neglectus-a and A. neglectus-b as the representatives of two lineages separated from population genetics were sequenced to compare with A. neglectus at the genome-level. Amphioctopus neglectus-a showed significant differences with A. neglectus, mainly reflected in gene length, intergenic regions and the secondary structure of tandem repeat motifs. Notably, two sequence deletions in mitogenomes of the two representative species were detected in different positions of major non-coding regions, which were the most distinct differences with A. neglectus. Pairwise genetic distances and the phylogenetic analysis supported the relationship of (A. neglectus-a + (A. neglectus + A. neglectus-b)). This study suggested that A. neglectus-a should be considered as a potential cryptic species of this complex, while A. neglectus-b needed further verification to be defined.

Chromosome-level analysis of the Crassostrea hongkongensis genome reveals extensive duplication of immune-related genes in bivalves

Crassostrea hongkongensis is a popular and important native oyster species that is cultured mainly along the coast of the South China Sea. However, the absence of a reference genome has restricted genetic studies and the development of molecular breeding schemes for this species. Here, we combined PacBio and 10 × Genomics technologies to create a C. hongkongensis genome assembly, which has a size of 610 Mb, and is close to that estimated by flow cytometry (~650 Mb). Contig and scaffold N50 are 2.57 and 4.99 Mb, respectively, and BUSCO analysis indicates that 95.8% of metazoan conserved genes are completely represented. Using a high-density linkage map of its closest related species, C. gigas, a total of 521 Mb (85.4%) was anchored to 10 haploid chromosomes. Comparative genomic analyses with other molluscs reveal that several immune- or stress response-related genes extensively expanded in bivalves by tandem duplication, including C1q, Toll-like receptors and Hsp70, which are associated with their adaptation to filter-feeding and sessile lifestyles in shallow sea and/or deep-sea ecosystems. Through transcriptome sequencing, potential genes and pathways related to sex determination and gonad development were identified. The genome and transcriptome of C. hongkongensis provide valuable resources for future molecular studies, genetic improvement and genome-assisted breeding of oysters.

Mitogenomic Perspectives on the Adaptation to Extreme Alkaline Environment of Amur ide (Leuciscus waleckii)

Amur ide (Leuciscus waleckii, Family Cyprinidae) is widely distributed in Northeast Asia. L. waleckii usually inhabits freshwater environments but can also survive in the Lake Dali Nur, one of the most extreme aquatic environments on the earth, with an alkalinity up to 50 mmol/L (pH 9.6). To investigate mechanisms of mitogenomic evolution underlying adaptation to extreme environments, we determined 30 complete mitogenomes that included Lake Dali Nur (alkaline environment, AL) population and Amur basin (freshwater environment, FW) population. Through phylogenetic and divergence time analysis, we found that AL and FW populations forming distinct two groups which were consistent with geographic divergence (the formation of Lake Dali Nur). In addition, we found that almost of the windows exhibited higher nucleotide diversity in FW population (avg 0.0046) than AL population (avg 0.0012). This result indicated that severe environment selection had remarkably reduced the genetic diversity of mitogenome in AL population and suggested that severe environment selection had remarkably reduced the genetic diversity of mitogenome in the AL population. Compared with the FW population (ω = 0.064), the AL population (ω = 0.092) had a larger mean ω (dN/dS ratios) value for the 13 concatenated mitochondrial protein-coding genes, indicating that the high alkaline tolerated group had accumulated more nonsynonymous mutations. These nonsynonymous mutations had resulted in slightly beneficial amino acid changes that allowed adaption to the severe conditions. This study provides an additional view to decipher the adaptive mitogenome evolution of L. waleckii of the high alkaline environment.

Thermotolerance Divergence Revealed by the Physiological and Molecular Responses in Two Oyster Subspecies of Crassostrea gigas in China

Investigating the physiological mechanisms of closely related species that exhibit distinct geographic distributions and thermal niches is essential for understanding their thermal tolerance capacities and local adaptations in view of climate warming. The variations in upper thermal limits (LT50) under acute heat shock and cardiac activity, standard metabolic rate (SMR), anaerobic metabolite production and molecular responses (expression of molecular chaperones and glycolysis metabolism genes) under increasing temperatures in two oyster subspecies were studied. The populations of two oyster subspecies, Crassostrea gigas gigas and C. gigas angulata, exhibit different latitudinal distributions along the northern and southern coastlines of China, respectively, which experience different environmental conditions. The LT50 was significantly higher, by ∼1°C, in the southern than in the northern oysters. In both subspecies, temperature increases had powerful effects on heart rate, SMR and gene expression. The southern oysters had the highest Arrhenius breakpoint temperatures for heart rate (31.4 ± 0.17°C) and SMR (33.09°C), whereas the heart rate (28.86 ± 0.3°C) and SMR (29.22°C) of the northern oysters were lower. The same patterns were observed for the Q 10 coefficients. More thermal sensitivity was observed in the northern oysters than in their southern counterparts, as the heat-shock proteins (HSPs) in the northern oysters were expressed first and had a higher induction at a lower temperature than those of southern oysters. Furthermore, different expression patterns of energetic metabolism genes (HK, PK, and PEPCK) were observed. In the northern oysters, increasing anaerobic glycolysis genes (PEPCK) and end products (succinate) were found at 36-43°C, indicating a transition from aerobic to anaerobic metabolism and a lower aerobic scope compared with the southern oysters. These two subspecies experience different environmental conditions, and their physiological performances suggested species-specific thermal tolerance windows in which the southern oysters, with mild physiological flexibility, had a higher potential capability to withstand heat stress. Overall, our results indicate that comparing and unifying physiological and molecular mechanisms can provide a framework for understanding the likely effects of global warming on marine ectotherms in intertidal regions.

Computational Characterization of the mtORF of Pocilloporid Corals: Insights into Protein Structure and Function in Stylophora Lineages from Contrasting Environments

More than a decade ago, a new mitochondrial Open Reading Frame (mtORF) was discovered in corals of the family Pocilloporidae and has been used since then as an effective barcode for these corals. Recently, mtORF sequencing revealed the existence of two differentiated Stylophora lineages occurring in sympatry along the environmental gradient of the Red Sea (18.5°C to 33.9°C). In the endemic Red Sea lineage RS_LinB, the mtORF and the heat shock protein gene hsp70 uncovered similar phylogeographic patterns strongly correlated with environmental variations. This suggests that the mtORF too might be involved in thermal adaptation. Here, we used computational analyses to explore the features and putative function of this mtORF. In particular, we tested the likelihood that this gene encodes a functional protein and whether it may play a role in adaptation. Analyses of full mitogenomes showed that the mtORF originated in the common ancestor of Madracis and other pocilloporids, and that it encodes a transmembrane protein differing in length and domain architecture among genera. Homology-based annotation and the relative conservation of metal-binding sites revealed traces of an ancient hydrolase catalytic activity. Furthermore, signals of pervasive purifying selection, lack of stop codons in 1830 sequences analyzed, and a codon-usage bias similar to that of other mitochondrial genes indicate that the protein is functional, i.e., not a pseudogene. Other features, such as intrinsically disordered regions, tandem repeats, and signals of positive selection particularly in StylophoraRS_LinB populations, are consistent with a role of the mtORF in adaptive responses to environmental changes.

Phylogeography of bivalve Meretrix petechialis in the Northwestern Pacific indicated by mitochondrial and nuclear DNA data

The marine clam Meretrix petechialis is an important economic shellfish species in Northwestern Pacific, but little is known about its phylogeographical pattern. Here, we analyzed 311 samples from 22 locations along the northwestern Pacific using combined profiling of one mitochondrial gene (the first subunit of cytochrome coxidase, COI) and one nuclear DNA marker (the internal transcribed spacer region 1, ITS-1) to investigate contemporary genetic structure and reconstruct phylogenetic history of this species. The results revealed that two distinct phylogeographic lineages dominated marginal seas-the East China Sea (ECS) and the South China Sea (SCS) respectively. The estimation of divergence time between two lineages was 2.1-3.8 Ma, corresponding to a period of the early Pleistocene to late Pliocene. The vicariance of the two lineages was connected to the historical isolation of marginal seas and sea surface temperature (SST) gradient, pointing that SST might play an important role in maintaining phylogeographical patterns of M. petachialis. Significant overlaps between two lineages were observed in 23° to 29° N, located at the adjacent area of the ECS and SCS, which might be promoted by the connectivity of China Coast Current. However, the influence of ocean currents on mixings between two lineages was limited. In comparison, significant relationships were found between genetic distances and geographic distances if the North and South populations were analyzed separately, result of which might be due to some small reciprocal, rotating flows along coastal areas and special geographical conditions.

DNA barcoding reveals that the common cupped oyster in Taiwan is the Portuguese oyster Crassostrea angulata (Ostreoida; Ostreidae), not C. gigas

The Pacific cupped oyster, Crassostrea gigas, is one of the major aquacultural shellfish species that has been introduced to Europe and America from its native source in the West Pacific. In Taiwan, the cultivated cupped oysters along the west coast have been identified as C. gigas for over centuries; however, several molecular phylogenetic studies have cast doubt upon the existence of this species in Taiwan and adjacent waters. Indeed, our analyses of mitochondrial cytochrome oxidase I (COI) sequences from 313 Crassostrea collected from 12 locations along Taiwanese and southern Chinese coastlines confirm that all samples were the Portuguese oyster, C. angulata, rather than C. gigas. Multiple lines of evidence, including haplotypic and nucleotide diversity of the COI gene, demographic history, and population genetics, suggest that Taiwanese C. angulata is unique, probably experienced a sudden population expansion after the Last Glacial Maxima around 20,000 years ago, and has a significantly limited genetic connectivity across the Taiwan Strait. Our study applies an extended sampling and DNA barcoding to confirm the absence of C. gigas in natural and cultivated populations in Taiwan and southern China, where we only found C. angulata. We highlight the importance of conserving the gene pool of the C. angulata population in Taiwan, particularly considering the current threats by large-scale environmental disturbances such as marine pollution, habitat destruction, and climate change.