Next Article in Journal
Active Coral Restoration to Preserve the Biodiversity of a Highly Impacted Reef in the Maldives
Next Article in Special Issue
Fish of Low Commercial Value in Lakes of Different Trophic Status (Poland)
Previous Article in Journal
Transcriptomic Analysis of Medicago truncatula under Long-Day Conditions
Previous Article in Special Issue
Exploring the Drivers of Spatiotemporal Patterns in Fish Community in a Non-Fed Aquaculture Reservoir
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Diversity
Volume 15
Issue 9
10.3390/d15091021
diversity-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Population Structure and Phylogeography of Marine Gastropods Monodonta labio and M. confusa (Trochidae) along the Northwestern Pacific Coast

by
Yuh-Wen Chiu
1,*,
Hor Bor
2,
Jin-Xian Wu
3,
Bao-Sen Shieh
4 and
Hung-Du Lin
5
1
Department of Biological Resources, National Chiayi University, Chiayi 600, Taiwan
2
Environmental Protection Department, New Taipei City Government, New Taipei 220, Taiwan
3
Guangzhou Key Laboratory of Subtropical Biodiversity and Biomonitoring, School of Life Science, South China Normal University, Guangzhou 510631, China
4
Department of Biomedical Science and Environmental Biology, Kaohsiung Medical University, Kaohsiung 807, Taiwan
5
The Affiliated School of National Tainan First Senior High School, Tainan 701, Taiwan
*
Author to whom correspondence should be addressed.
Diversity 2023, 15(9), 1021; https://doi.org/10.3390/d15091021
Submission received: 11 May 2023 / Revised: 12 September 2023 / Accepted: 14 September 2023 / Published: 17 September 2023
(This article belongs to the Special Issue Evaluating the Influence of Environmental Variables on Fish Ecology and Diversity)
Browse Figures
Review Reports Versions Notes

Abstract

:
The genetic structure and demographic history of marine organisms are influenced by biological and ecological features, oceanic currents, and the paleo-geological effects of sea-level fluctuations. In this study, we used mitochondrial COI + 16S gene analysis to investigate the phylogeographic pattern and demography of Monodonta labio and M. confusa in Taiwan, the Ryukyu Islands, Japan, mainland China, and Korea. Our genetic analysis identified two major lineages that were not evident from the distribution patterns of different populations. The Taiwan Strait, which acted as a land bridge during Pleistocene glaciations, and the lack of strong dispersal barriers to gene flow between ocean basins after glaciations shaped the phylogeographic pattern. The genetic differentiation in the Ryukyu Islands was influenced by the specialist-generalist variation hypothesis and the Kuroshio Current. Bayesian skyline plot analyses suggested that the effective population size of M. labio and M. confusa rapidly increased approximately 0.1 and 0.075 million years ago, respectively. Our approximate Bayesian computation analysis suggested that all M. labio and M. confusa populations experienced a decline in population size following a recent population expansion and constant size, respectively. Our study provides a baseline for future investigations of the biogeographical patterns of marine gastropods in the Northwest Pacific and offers valuable insights for the management, sustainable resource utilization, and conservation of this species.

1. Introduction

Global climate conditions have changed during glacial-to-interglacial periods, leading to the oscillation of the sea level, oceanic current pattern, and life history traits (e.g., dispersal ability and life cycle), which could have an unavoidable effect on the phylogeography, demographic history, and population genetic diversity of coastal marine systems [1,2]. The Northwest Pacific, mainly including the East China Sea (ECS), the South China Sea (SCS), the Yellow Sea, and the Sea of Japan, is an epicenter of marine biodiversity and provides an excellent model system to investigate as a major glacial refuge, driving population differentiation at different spatial and temporal scales in some marine species during glacial periods of the Plio-Pleistocene [3]. During major glaciations, when there is a sea level decline of approximately 120–140 m below the present level, the northwestern (NW) Pacific is thought to form a continuous land mass, effectively closing the sea passage (the Korean Strait and Taiwan Strait). Because of their geographical isolation, the East China Sea (ECS) and the South China Sea (SCS) are partially enclosed, with connections to the Pacific occurring through the Okinawa Trough and the Bashi Strait in a restricted manner. The emergence of the Strait as a barrier is one of the significant factors that have shaped the current phylogeographical patterns and the presence of deeply divergent lineages in marine species, such as marine gastropods (Reticunassa festiva) [4], grey mullet (Mugil cephalus) [5], spotted scat (Scatophagus argus) [6], and cutlassfish (Lepturacanthus savala) [7]. The genetic structure of numerous marine organisms may be influenced by larval dispersal, which is believed to be affected by the ocean current system in the northwest Pacific. In general, ocean currents can facilitate the migration of marine organisms over long distances and enhance population connectivity. Conversely, convergent ocean currents can create an effective barrier to gene flow to a certain degree [8]. Prominent ocean currents in the NWP and its marginal seas include the cyclonic circulation of the Kuroshio Current (KC) and coastal currents. The Kuroshio Current is characterized by a warm current that originates from the Philippines and flows northeastward to the Pacific coast of southern mainland Japan via Taiwan and the Ryukyu Islands. The NWP marginal seas are characterized by complex coastal currents, such as the China Coastal Current (CCC), Tsushima Current (TC), and Korean Coastal Current (KCC). By altering the dynamics of larval supply, enhancing population connectivity, or producing different outcomes from the sea surface temperature gradient, ocean currents play an important role in shaping genetic structure. Meanwhile, some phylogeographic studies in marine organisms have shown a lack of genetic structuring among distant populations with continuous gene flow via current-driven larval dispersal, such as in bigeye tuna (Thunnus obesus) [9], Chinese beard eel (Cirrhimuraena chinensis) [10], and the large yellow croaker (Larimichthys crocea) [11].
The marine trochid gastropod Monodonta belongs to the family Trochidae and is commonly found in the mid-intertidal zone. It exhibits a wide distribution in the Pacific Northwest. This herbivorous snail is known to feed on microalgae [12]. The genus Monodonta consists of eleven species and one subspecies, of which eight species and one subspecies are distributed in East Asia: M. australis, M. canalifera, M. confusa, M. glabrata, M. labio, M. viridis, M. neritoides, M. perplexa perplexa, and M. perplexa boninensis (listed in Donald [12]), including two species (M. confusa and M. perplexa) that undergo a planktonic stage lasting 3 days in their life cycle [13]. Biological traits such as effective dispersal capability through a planktonic larval stage and life cycle can strongly influence the genetic differentiation and speciation of marine organisms [14]. In general, the presence of a long-lived planktonic stage is often considered important for dispersal and lower geographic variations, such as Sicyopterus japonicus [15] and Thais clavigera [16]. In contrast, a short planktonic stage in marine organisms has revealed a clear genetic structure, such as in the Japanese turban shell (Turbo (Batillus) cornutus) [17] and the Moon Turban Snail (Lunella granulate) [18]. Previous studies of M. labio detected intraspecific variations in the shell morphology and genetic structure in NWP [19]. According to Zhao et al. [20], M. labio can be divided into five primary clades distributed along the Chinese coastline, and the Dongshan land bridge served as a separation between the ancient East China Sea and the ancient South China Sea, forming two possible refugia. Another similar study proposed that M. labio and M. confusa were primarily separated into different lineages with a boundary between the Japanese mainland and the Ryukyu Islands [21]. The genetic structure of M. labio and M. confusa is also influenced by ecological factors such as habitat specificity and adaptability. Molecular studies have indicated that M. labio displays lower genetic diversity and higher genetic differentiation compared to M. confusa, potentially influenced by specific ecological factors [22]. This pattern is known as the “specialist-generalist variation hypothesis” (SGVH) [22]. Monodonta labio and M. confusa have been regarded as excellent model systems in the NWP for examining the impact of geographical and ecological isolation on genetic divergence. Previous studies of genetic structure in M. labio and M. confusa focused on mainland China, the Ryukyu Islands, and Japan, but provided little information on Taiwan. Furthermore, the phylogenetic processes of M. labio and M. confusa on mainland China, Taiwan, the Ryukyu Islands, and Japan have not been thoroughly or convincingly elucidated, and there have been no investigations into the population relationships between Taiwan and other regions. Statistical phylogeography is a powerful approach to understanding the processes that have shaped the spatio-temporal demographic and distributional changes in M. labio and M. confusa. Relying on the former molecular results, our objective was to use ABC simulations, which offer a framework for testing competing hypotheses to determine the best model of population divergence and demographics that fits the patterns of key demographic parameters, such as changes in effective population size through time. This study aimed to compare the genetic variation between two co-occurring species of intertidal snails, M. labio and M. confusa, from the genus Monodonta, found on the Northwestern Pacific Coast. To achieve this, we utilized the mitochondrial COI + 16S genes to establish the phylogeographic patterns of these species across mainland China, Taiwan, the Ryukyu Islands, and Japan. The three main objectives were to: (1) determine the genetic diversity and structure of M. labio and M. confusa; (2) reconstruct the phylogeographic patterns of M. labio and M. confusa, understand the influence of paleoclimate oscillations on their current distribution; and (3) examine the demographic response of M. labio and M. confusa to repeated climate changes during the Quaternary period.

2. Materials and Methods

2.1. Sample Collection and Sequencing

The mitochondrial COI + 16S gene, which was 1340 bp long, was sequenced from 85 specimens of M. labio and 41 specimens of M. confusa collected from thirteen localities in Taiwan, mainland China, Jeju Island in Korea [was missing], and representing the Ryukyu Islands and more northern localities within Japan (Table 1 and Figure 1). Based on their geographical distribution, the thirteen populations were classified into four regions: Taiwan Island (Audi, AD; Miaoli, WW; Tainan, CG; Taitung, TD); Ryukyu Islands (Okinawa, OK; Ishigaki, SY); Japan (Inubosaki, CF; Ine, YG; Kagoshima, SI); Mainland China (Kinmen, KM; Matsu, MZ); and Korea (Jeju, GJ) (Table 1 and Figure 1). Table 1 provides details on the location and sample sizes.
Samples were manually collected from the field sites, fixed, and stored in 95% ethanol. All animal experiments were conducted in the laboratory of Yuh-Wen Chiu at National Chiayi University, following the guidelines and approval of the Animal Research and Ethics Committee of the same university. The use of gastropods as experimental animals in Taiwan does not require a license. Muscle tissue was used to extract genomic DNA using a tissue and cell genomic DNA purification kit ((GeneMark Technology Co., Tainan, Taiwan)). Fragments of the mitochondrial cytochrome subunit I gene (COI) and 16S rRNA region (16S) were amplified by polymerase chain reaction (PCR) using the primers LCO-1490 (5′- GGTCAACAAATCATAAAGATATTGG-3′) and HCO-2198 (5′-TAAACTTAGGGTGACCAAAAATCA-3′) [23] for COI and 16Sar (5′-CGCCTGTTTATCA AAAACAT-3′) and 16Sbr (5′-CCGGTCTGAACTCAGATCACGT-3′) [24] for 16S. The polymerase chain reaction (PCR) was performed in a 25 μL reaction volume according to the conditions specified in Chiu et al. (2013) [18]. Subsequently, the PCR products were run on an ABI 377 automated sequencer (Applied Biosystems, Foster City, CA, USA). All sequences have been deposited in GenBank under the following accession numbers: OQ880826-OQ 880919. Other sequences of M. labio and M. confusa are available on Genbank [20,21,22], with additional sequences of M. perplexa selected as the outgroup.

2.2. Sequence Alignment and Data Analysis

We aligned the COI + 16S genes of the mtDNA sequences using CLUSTAL X v. 2.0 software [25] and manually adjusted the alignment. To calculate the genetic diversity of each population, we used DnaSP v5.0 software [26] and determined the number of haplotypes (N), haplotype diversity (h), private haplotypes (ph), nucleotide diversity (π), and nucleotide diversity (θ) [27]. To examine the spatial partitioning of genetic variation among populations, pairwise FST values and analysis of molecular variance (AMOVA) were computed using Arlequin v3.5 [28]. For the hierarchical analysis, populations were grouped and defined based on geographical factors: (1) two geographical groups, Pacific coast and other, were primarily divided by the Kuroshio Current; (2) three geographical groups, Ryukyu Islands, Taiwan, mainland China (Taiwan Strait), and Japan, were primarily divided by the ocean region; and (3) five geographical groups, Ryukyu Islands, Taiwan Island, mainland China, Japan, and Korea, were primarily divided by the sampling region.
Phylogenetic relationships were reconstructed using neighbor-joining (NJ), maximum likelihood (ML), and Bayesian inference (BI) methods, implemented in MEGA-X [29], PhyML 3.0 [30], and MrBayes v3.1.2 [31], respectively. The best-fit nucleotide substitution model was determined using Smart Model Selection with the Akaike Information Criterion (AIC) [32] in PhyML 3.0, and the GTR + G + I model [33] was selected as the best nucleotide replacement model. We conducted BI analyses with two concurrent runs of approximately 1,500,000 iterations, with samples recorded every 5000 iterations. The first 25% of trees were discarded as burn-in, and the remaining trees were used to generate a 50% majority rule consensus. For the ML analysis, we used the GTR + G + I substitution model for the nucleotide dataset, and bootstrap resampling was performed using the rapid option with 1000 iterations. Branch support was evaluated with 100 pseudoreplicates. For NJ analysis, we applied the Kimura 2-parameter model, and bootstrapping was performed with 1000 replications. Haplotype networks were constructed using the MINSPNET algorithm implemented in Popart [34].

2.3. Historical Demography, Divergence Time Estimation and Biogeographic Analysis

To test the assumption of neutrality, we analyzed historical demographic scenarios using Tajima’s D test [35] and Fu’s Fs test [36], and computed mismatch distributions in DnaSP v5.0 [26]. Bayesian skyline plot (BSP) analyses were performed in BEAST v1.8.2 for each lineage in M. labio [37] to determine effective population size changes over time. To ensure convergence of all parameters (ESSs > 200), we ran >200,000,000 MCMC iterations and discarded the first 10% of samples for each chain as burn-in. Plots for each analysis were generated using Tracer v1.6 [38]. Additionally, we used BEAST v1.8.2 software [37] to calculate the time to the most recent common ancestor (TMRCA). In this study, we adopted a mutation rate of 2.4% per million years calculated from two Tegula species separated by the Isthmus of Panama and calibrated for the mtDNA COI genes for population expansion [17,39]. We utilized MIGRATE v3.2.6 to ascertain the historical migration rate direction (M = m/μ, where m = migration rate and μ = the per-locus mutation rate) between Taiwan, the Ryukyu Islands, and Japan (including GJ) [40]. Maximum likelihood inference was performed in MIGRATE, comparing a full migration model (where θ and M were estimated jointly from the data) to a restricted model (where θ was averaged and M was assumed to be symmetrical among populations). Each run included five replicates of ten short chains (with 10,000 genealogies sampled) and three long chains (with 100,000 genealogies sampled), discarding the first 10,000 genealogies. The model was run five times to ensure parameter estimates’ convergence, assuming a migration matrix model with unequal population sizes and different migration rates [40].
We employed approximate Bayesian Computation (ABC) methods with DIYABC v.2.0.4 [41] to investigate the potential historical demographic scenarios of M. labio and M. confusa. We defined five demographic scenarios, as suggested by Cabrera and Palsboll [42] (Figure 2), to explore the demographic histories of M. labio and M. confusa. The five scenarios were as follows: (1) a constant size model (CON); (2) a bottleneck event model (DEC); (3) a recent expansion model (INC); (4) an expansion followed by a decrease model (INCDEC); and (5) a decrease followed by an expansion model (DECINC). We generated three million datasets for each scenario and calculated summary statistics using DIYABC [40]. The posterior probability of each scenario was obtained using both the direct approach and the logistic regression approach implemented in DIYABC [40].

3. Results

3.1. Genetic Diversity

A total of 1340 base pairs (bp) of mtDNA COI + 16S gene sequences from 85 specimens were analyzed and found to contain 25.0% adenine, 32.1% thymine, 20.7% cytosine, and 22.1% guanine (42.8% GC content) in M. labio. Out of the 41 specimens analyzed in M. confusa, the DNA composition revealed 25.2% adenine, 32.4% thymine, 20.3% cytosine, and 22.1% guanine, resulting in a GC content of 42.4%. The lengths of the COI and 16S markers were 769 and 571 bp, respectively. A total of 61 haplotypes were identified, and the average haplotype diversity was high (0.987), ranging from 0.901 (SY) to 1.000 (AD, WW, MZ, and GJ) in M. labio. In M. confusa, we identified a total of 33 haplotypes, and the average haplotype diversity was high at 0.983, ranging from 0.867 (CF) to 1.000 (AD, TD, and YG). The average nucleotide diversity (θπ) was low (0.0131), ranging from 0.0012 (YG) to 0.0042 (GY) in M. labio (Table 1). In M. confusa, the average nucleotide diversity (θπ) was found to be low at 0.0082, with values ranging from 0.0031 (CF) to 0.0093 (AD), as shown in Table 1. There were only six haplotypes that occurred in more than one population, with one haplotype (H29) observed in five populations and the others only shared by specimens from two populations of M. labio (Table S1). In M. confusa, only four haplotypes were found to occur in more than one population, being shared by specimens from two populations (Table S1). Only two populations (OK and SY) in the Ryukyu Islands do not share any haplotypes with other populations of M. labio.

3.2. Phylogenetic Reconstruction and Genetic Structure

The significant correlation between phylogeny and geography was evident through the considerably larger NST values compared to GST (0.914 and 0.035 in M. labio, 0.266 and 0.026 in M. confusa, respectively), consistent with findings from Pons and Petit [43]. Using COI + 16S for phylogenetic analysis, we identified three main lineages (IA, IB, and IC) in M. labio. Lineage IA was exclusively found in populations OK and SY in the Ryukyu Islands. Lineage IB consisted of individuals from populations in Taiwan and mainland China, while lineage IC was distributed among populations from Taiwan Island, mainland China, Japan, and one population (GJ) from Korea. Based on the phylogenetic tree, three main lineages (IIA, IIB, and IIC) were identified in M. confusa. Lineages IIA and IIB were present in populations from Taiwan Island and mainland China, while lineage IIC exclusively occurred in populations from Japan. The NJ, ML, and BI methods yielded the same tree topologies in the phylogenetic analysis (Figure 3), which supported the haplotype minimum-spanning network topology and showed no clustering corresponding to the sampling site (Figure 4). The topological relationships obtained from the phylogenetic analysis with other sequences of M. labio and M. confusa from GenBank also support the formation of two major lineages (Figure S1).
Diversity 15 01021 g003
Figure 3. Reconstructed phylogeny of Monodonta labio and M. confusa using neighbor-joining (NJ) analyses based on the mitochondrial COI + 16S gene among 13 populations using 94 haplotypes The numbers above each branch are the posterior probabilities for bootstrap values for the NJ, ML, and Bayesian analyses.
Figure 3. Reconstructed phylogeny of Monodonta labio and M. confusa using neighbor-joining (NJ) analyses based on the mitochondrial COI + 16S gene among 13 populations using 94 haplotypes The numbers above each branch are the posterior probabilities for bootstrap values for the NJ, ML, and Bayesian analyses.
Diversity 15 01021 g003
The range of pairwise FST values between populations ranged from −0.047 (between MZ and WW) to 0.985 (between SI and PH), with a mean value of 0.911 in M. labio. In M. confusa, the pairwise FST values between populations varied from −0.138 (between MZ and TD) to 0.605 (between CF and PH), with an average value of 0.265. Moreover, the pairwise FST values between populations from different regions were relatively high and significant, indicating that high differentiation existed among these populations (Table 2). The AMOVA results indicate significant genetic structure at scenarios II and III in M. labio and M. confusa (Table 3). The results from the five geographic regions (Scenario III: Taiwan, mainland China, Ryukyu Islands, Japan, and Korea) demonstrated that significant spatial genetic structuring among groups was 59.93% (FCT = 0.599, p < 0.001) but was 31.74% (FSC = 0.792, p < 0.001) among populations within groups and 8.31% (FST = 0.916, p < 0.001) within populations in M. labio (Table 3). Furthermore, the findings from the five geographic regions (Scenario III) revealed significant spatial genetic structuring among groups at 33.51% (FCT = 0.335, p < 0.001). The genetic structuring among populations within groups was minimal at −0.27% (FSC = −0.004, p = 0.546), with the majority of the variation, 66.75%, being attributed to within-population differences (FST = 0.332, p < 0.001) in M. confusa (Table 3). However, AMOVA from the two geographic regions (Pacific Coast and others) identified only −5.50% and −10.46% of the variants to be present among groups, 72.41% and 35.55% of the variation among populations within groups, and 33.10% and 74.90% of the variation within populations in M. labio and M. confusa, respectively (Table 3).
Estimates of gene flow calculated with the MIGRATE program indicated the possibility of genetic exchange among the southern region (Taiwan and mainland China group), the central region (Ryukyu Islands group), and the northern region (Japan and Korea group). Asymmetrical gene flow among the three regions was found, indicating that individuals from the southern region may have colonized the northern region. The estimate of historical gene flow from the northern region to the southern region and the southern region to the central region was larger in M. labio (M S→N = 215.1, 97.5% CI: 78~328) and M C→N = 101.8, 97.5% CI: 25.3~184.0, where the suffixes S, C, and N indicate the southern region, the central region, and the northern region, respectively. In M. confusa, the estimation of historical gene flow indicated similarity in gene flow from the northern region to the southern region and from the southern region to the northern region.

3.3. Molecular Dating and Historical Demographic Expansion

The estimated Tajima’s D and Fu’s Fs indexes were not consistent in the overall population of M. labio, with positive Tajima’s D (1.51946) and significantly negative Fu’s Fs (−1.931). In M. confusa, the estimated Tajima’s D and Fu’s Fs indexes showed inconsistency. Tajima’s D was nonsignificantly negative (−1.16158), while Fu’s Fs was significantly negative (−15.823). Multiple-modal mismatch distributions were observed in the total population, indicating no recent demographic expansion for the populations of M. labio and M. confusa. The Bayesian skyline plot analysis, which assumed a substitution rate of 2.4% per million years, revealed that M. labio populations underwent a population decline followed by rapid population growth around 0.1 million years ago, as shown in Figure 5a. Figure 5b indicates that populations of M. confusa experienced a population growth event around 0.075 million years ago. These results showed that the history of population size changes in M. labio and M. confusa is complex. We used the five scenarios of ABC analysis to determine the possible demographic history of M. labio. In the ABC modeling, Scenario 5 (DECINC model) was highly favored (posterior probability = 1.0000 [1.0000, 1.0000]) compared with the other scenarios in M. labio. In M. confusa, Scenario 1 (CON) was highly favored (posterior probability = 0.9698 [0.9613, 0.9783]) compared with the other scenarios. This scenario suggested that M. labio experienced a reduction in the effective population size in the past, followed by a single instantaneous increase in population size. The estimated time to coalescence in M. labio was estimated at 2.885 ± 0.441 myr based on results obtained with BEAST. The estimated time to the most recent common ancestor (TMRCA) of M. confusa is approximately 0.286 ± 0.103 myr.

4. Discussion

4.1. Genetic Diversity

Maintaining genetic diversity is crucial for populations to persist and enable them to adapt and evolve in response to environmental changes. Several factors can affect the genetic diversity of marine organisms, including overfishing, selection, pollution, the presence of dispersal barriers, dispersal capacity, genetic drift within the population, and inbreeding. In this study, analyzing the genetic diversity of M. labio and M. confusa populations revealed higher haplotype diversity (0.987, 0.983) and nucleotide diversity (0.0131, 0.0082) than previous studies conducted in mainland China and Japan (e.g., [20,21,22]). In our results, the remarkably lower nucleotide diversity may be due to the lower effective population size of M. labio and M. confusa in the Ryukyu Islands and Japan, as compared to that of Taiwan and mainland China (Table 1). Monodonta labio and M. confusa are widely distributed along the coast of East Asia. In our studies, populations in Taiwan and mainland China at the center of the distribution are likely to have high genetic diversity due to gene flow from neighboring populations in the Ryukyu Islands and Japan that are better adapted to different conditions. The combination of low nucleotide diversity and high haplotype diversity has often been linked to founder effects or population bottlenecks [43]. In previous studies, M. labio was identified as a habitat specialist, whereas M. confusa showed no habitat specificity. Additionally, M. labio exhibited lower genetic diversity and higher genetic differentiation among populations when compared to M. confusa. In our studies, we observed a comparable pattern of nucleotide diversity in M. labio, which showed higher nucleotide diversity than M. confusa. This suggests that our study had a larger sampling range compared to previous studies [22]. However, the abundance and uniform distribution of private haplotypes within M. labio and M. confusa populations indicate a large effective population size, likely due to their wide distribution throughout East Asia.

4.2. Population Structure and Demographic History

Pairwise FST values suggested genetic differentiation was substantially greater among populations between regions than among populations within regions. Furthermore, our AMOVA and IBD analyses among populations of M. labio and M. confusa suggested lower proportions of genetic admixture. The AMOVA results revealed that the genetic differentiation of the total genetic variation among groups in M. labio and M. confusa from the three geographic regions (Scenario II) accounted for 68.40% (FCT = 0.683, p = 0.016) and 43.25% (FCT = 0.432, p = 0.032), respectively. However, the variation among populations within groups was 24.36% (FSC = 0.770, p < 0.000) for M. labio and −2.04% (FSC = 0.035, p = 0.686) for M. confusa, while the variation within populations was 7.24% (FST = 0.927, p < 0.001) for M. labio and 66.75% (FST = 0.332, p < 0.000) for M. confusa (Table 2). These results were consistent with our IBD results. The life cycle of M. labio includes a relatively short larval stage, which might result in genetic differentiation [22]. The findings of our study align with previous research, indicating that M. labio displays higher genetic differentiation among populations than M. confusa [22]. The higher genetic differentiation observed in M. labio is consistent with previous studies that reported significant differences between populations in Japan and mainland China [20,21,22]. Yamazaki et al. [22] indicated that the genetic structure of M. labio and M. confusa seems to be strongly influenced by geographical factors, such as the dispersal barrier between the Japanese Islands and the Ryukyu Islands. In our study, lineage IA was a monophyletic group only distributed in the Ryukyu Islands. We suggested that M. labio has relatively higher genetic differentiation and lower nucleotide diversity in the Ryukyu Islands, which have experienced significant founder effects and bottlenecks. Within M. confusa, lineage IIC constituted a monophyletic group solely distributed in the Japan Islands, implying relatively higher genetic differentiation between Taiwan Island, mainland China, and Japan.
In general, genetic differentiation among populations is influenced by their dispersal abilities in marine invertebrate populations [44]. The results from MIGRATE suggested asymmetrical gene flow, with higher migration rates from the southern and central regions to the northern region. Such a scenario would occur in the context of the biological characteristics and population demographics of M. labio. In the spring, M. labio lays its eggs in shallow coastal waters and exhibits a planktonic larval stage that facilitates dispersal, following the coastal current that flows northeastward through the NWP marginal seas in spring and summer [45]. Kuroshio Current systems probably played a more predominant role in shaping the genetic homogeneity among populations of M. labio in NWP marginal seas [22]. The influence of combined southward-running currents (the China Coastal Current (CCC)) on the population genetic connectivity of M. labio is presumably weaker than the solely northward-flowing Kuroshio Current. The results confirm that the dynamic system associated with the variation in dispersal ability may be affected by the dominant Kuroshio current systems of the Pacific Ocean. Such a scenario might apply to other marine species in NWP marginal seas, such as Sicyopterus japonicas [15] and Sargassum fusiforme [46]. In ecological terms, M. confusa demonstrates a broad tolerance for different wave exposures and exhibits no specific habitat preference. Regarding population differentiation, M. confusa shows lower levels compared to M. labio. Based on the results from the MIGRATE analysis, the historical gene flow estimation indicates similarity in gene flow between populations. Therefore, we conclude that the complex currents in the Western Pacific consistently influence the dispersal of M. confusa.
Neutrality indices by Tajima’s D tests and multi-modal mismatch distribution analyses suggested that M. labio and M. confusa populations indicate no recent demographic expansion, but the significantly negative value of Fu’s Fs tests indicated recent demographic expansion. In general, Fu’s Fs tests are more sensitive at detecting demographic events than Tajima’s D tests. However, high haplotype diversity and low nucleotide diversity were found in the population, and the results of Bayesian skyline plot analysis were consistent with a recent population expansion [35,36,47]. The Bayesian skyline plot analyses unveiled a late Pleistocene demographic expansion for M. labio, which commenced approximately 0.1 million years ago. Additionally, M. confusa populations exhibited a population growth event around 0.075 million years ago. Similar results were observed in previous studies of M. labio and M. confusa in mainland China, as well as in other marine organisms such as S. macrorhynchos [48] and Nemipterus bathybius [49]. The demographic history of M. labio and M. confusa is more intricate, as evidenced by our approximate Bayesian Computation (ABC) analysis. The results suggested that all populations of M. labio experienced a decline in population size after a recent expansion, while M. confusa populations remained constant in size.

4.3. Phylogeography of Monodonta labio

Yamazaki et al. [21] noted that the genus Monodonta undergoes a planktonic stage in their lifecycle, and the planktonic period is reported as three days, revealing that species were likely affected by geographic isolation with regard to their intraspecific divergence, population structure, and geographical distribution. Tempo-spatial genetic architecture is observed between two mitochondrial lineages based on the haplotype network and phylogenetic analysis of M. labio and M. confusa. We incorporated sequences published by other studies into our analysis and found that they can be divided into two groups, lineages I (M. labio) and II (M. confusa). Among the individuals in the lineage II group, Zhao et al. [20] identified them as M. labio (refer to “clade D” in Zhao et al. [20]), while Yamazaki et al. [20,21] identified them as M. confusa (refer to “clade D” in Yamazaki et al. [20,21]). Due to the wide distribution range of this group, in our study, to avoid confusion in the classification system and to focus on the population’s historical dynamics, we agree with the viewpoint of Yamazaki et al. [21] and consider the individuals in lineage II as M. confusa (Figure S1, Table S2). Historically, many phylogeographic scenarios for marine organisms have responded to sea-level fluctuations during Pleistocene glaciations. The Taiwan Strait acted as a land bridge due to the fact that the sea level fell to 130 m below the present sea level (generally shallower than 60 m) within glacial periods [50]. The estimated time to the most recent common ancestor (TMRCA) of lineages I (M. labio) and II (M. confusa) is estimated at 2.885 myr and 0.441 myr, respectively; however, no direct evidence demonstrates that the Taiwan Strait Land Bridge, which acted as a physical barrier for the M. labio and M. confusa populations, existed during the period [20]. Many previous studies of some marine organisms suggested that two distinct genetic lineages were consistent with our findings in the Northwest Pacific Ocean, such as Macrobrachium nipponense [51], Terapon jarbua [52], Chinese black sleeper (Bostrychus sinensis) [53], and spotted scat (S. argus) [6]. The lineage I (M. labio) was subdivided into three sub-lineages (IA, IB, and IC): the Ryukyu Islands and the populations in Taiwan, mainland China, and Japan. The same pattern of genetic differentiation in the Ryukyu Islands has previously been observed in the same species [20,21]. In other studies, individuals of sub-lineage IA were also found to be distributed in Vietnam and the Gulf of Tonkin region (Figure S1) [20]. Sub-lineages IB are mainly distributed in Taiwan, Penghu, and mainland China, while sub-lineages IC are mainly distributed in Japan and South Korea [20,21,22]. Lineages II (M. confusa) have the widest distribution in the Northwest Pacific region, with distribution across various regions including Japan, South Korea, mainland China, and Taiwan, except for the Ryukyu Islands (Figure S1, Table S2) [20,21,22]. In previous studies of M. confusa, the research locations were primarily focused on Japan and the Ryukyu Islands, without including Taiwan and mainland China. The genetic differentiation between Japan and the Ryukyu Islands was not significant. In our study, among the three sub-lineages, two (IIA and IIB) are mainly distributed in Taiwan and mainland China, while the IIC sub-lineage is primarily distributed in Japan. This indicates a certain level of differentiation among the sub-lineages, which we believe is a result of distance and ocean currents. We suggested that the specialist-generalist variation hypothesis (SGVH) is the most important factor for genetic differentiation and that the Kuroshio Current is an effective influencer of the genetic differentiation of M. labio, supporting previous studies [21,54].

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/d15091021/s1. Table S1: The distribution information of the shared haplotypes (Hm1–Hm 94) (COI + 16S) p indicates the number of private haplotypes in each population. Table S2: The haplotype information downloaded from NCBI. Figure S1: The phylogeny of Monodonta labio (Lineage I) and M. confusa (Lineage II) was reconstructed using neighbor-joining (NJ) analyses based on the mitochondrial COI + 16S gene. This analysis involved a total of \mbox{126 haplotypes}, including 68 sequences obtained from NCBI and 94 haplotypes derived from our study. The numbers above each branch are the posterior probabilities for bootstrap values for the NJ, ML, and Bayesian analyses. The haplotypes marked in black on the phylogenetic tree are haplotypes downloaded from NCBI.

Author Contributions

Conceptualization, Y.-W.C. and B.-S.S.; methodology, H.B.; software, H.B. and J.-X.W.; validation, H.B., B.-S.S. and J.-X.W.; investigation, Y.-W.C. and H.B.; resources, Y.-W.C.; data curation, Y.-W.C., H.B. and H.-D.L.; writing—original draft preparation, Y.-W.C. and H.-D.L.; writing—review and editing, Y.-W.C. and H.-D.L. project administration, Y.-W.C.; funding acquisition, Y.-W.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by National Science and Technology Council of Taiwan, MOST 110-2621-B-415-004 and MOST 105-2621-B-006-005.

Institutional Review Board Statement

All animal experiments were conducted in the laboratory of Yuh-Wen Chiu at Nation-al Chiayi University, following the guidelines and approval of the Animal Research and Ethics Committee of the same university.

Data Availability Statement

All sequences have been deposited in GenBank under the following accession numbers: OQ880826-OQ 880919.

Acknowledgments

We thank Da-Ji Huang for his assistance in the field and help with the sample collection. Thanks also to Shih-Hsiung Liang for their invaluable review of and improvements to our manuscript. We also thank the anonymous reviewers and associate editor for their helpful comments.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hewitt, G. The genetic legacy of the Quaternary ice ages. Nature 2000, 405, 907–913. [Google Scholar] [CrossRef] [PubMed]
  2. Galindo, H.M.; Olson, D.B.; Palumbi, S.R. Seascape genetics: A coupled oceanographic-genetic model predicts population structure of Caribbean corals. Curr. Biol. 2006, 16, 1622–1626. [Google Scholar] [CrossRef] [PubMed]
  3. Briggs, J.C.; Bowen, B.W. A realignment of marine biogeographic provinces with particular reference to fish distributions. J. Biogeogr. 2012, 39, 12–30. [Google Scholar] [CrossRef]
  4. Yang, Y.; Qi, L.; Kong, L.; Li, Q. Phylogeography of the marine gastropod Reticunassa festiva complex (Nassariidae) in the coast of China. J. Shellfish Res. 2020, 39, 419–431. [Google Scholar] [CrossRef]
  5. Shen, K.N.; Jamandre, B.W.; Hsu, C.C.; Tzeng, W.N.; Durand, J.D. Plio-Pleistocene sea level and temperature fluctuations in the northwestern Pacific promoted speciation in the globally-distributed flathead mullet Mugil cephalus. BMC Evol. Biol. 2011, 11, 83. [Google Scholar] [CrossRef]
  6. Yan, Y.R.; Hsu, K.C.; Yi, M.R.; Li, B.; Wang, W.K.; Kang, B.; Lin, H.D. Cryptic diversity of the spotted scat Scatophagus argus (Perciformes: Scatophagidae) in the South China Sea: Pre-or post-production isolation. Mar. Freshw. Res. 2020, 71, 1640–1650. [Google Scholar] [CrossRef]
  7. Gu, S.; Yi, M.R.; He, X.B.; Lin, P.S.; Liu, W.H.; Luo, Z.S.; Lin, H.D.; Yan, Y.R. Genetic diversity and population structure of cutlassfish (Lepturacanthus savala) along the coast of mainland China, as inferred by mitochondrial and microsatellite DNA markers. Reg. Stud. Mar. Sci. 2021, 43, 101702. [Google Scholar] [CrossRef]
  8. Wang, L.; Wan, Z.Y.; Lim, H.S.; Yue, G.H. Genetic variability, local selection and demographic history: Genomic evidence of evolving towards allopatric speciation in Asian seabass. Mol. Ecol. 2016, 25, 3605–3621. [Google Scholar] [CrossRef]
  9. Chiang, H.C.; Hsu, C.C.; Lin, H.D.; Ma, G.C.; Chiang, T.Y.; Yang, H.Y. Population structure of bigeye tuna (Thunnus obesus) in the South China Sea, Philippine Sea and western Pacific Ocean inferred from mitochondrial DNA. Fish. Res. 2006, 79, 219–225. [Google Scholar] [CrossRef]
  10. Li, H.; Lin, H.; Li, J.; Ding, S. Phylogeography of the Chinese beard eel, Cirrhimuraena chinensis Kaup, inferred from mitochondrial DNA: A range expansion after the last glacial maximum. Int. J. Mol. Sci. 2014, 15, 13564–13577. [Google Scholar] [CrossRef]
  11. Liu, Q.; Lin, H.; Chen, J.; Ma, J.; Liu, R.; Ding, S. Genetic variation and population genetic structure of the large yellow croaker (Larimichthys crocea) based on genome-wide single nucleotide polymorphisms in farmed and wild populations. Fish. Res. 2020, 232, 105718. [Google Scholar] [CrossRef]
  12. Donald, K.M.; Kennedy, M.; Spencer, H.G. The phylogeny and taxonomy of austral monodontine topshells (Mollusca: Gastropoda: Trochidae), inferred from DNA sequences. Mol. Phylogenet. Evol. 2005, 37, 474–483. [Google Scholar] [CrossRef] [PubMed]
  13. Sasaki, R. Assessment and appearance of planktonic larvae of Haliotis (Nordotis) discus hannai around Kesennuma Bay. Aquaculture 1985, 32, 199. [Google Scholar]
  14. Palumbi, S.R. Genetic divergence, reproductive isolation, and marine speciation. Annu. Rev. Ecol. Evol. Syst. 1994, 25, 547–572. [Google Scholar] [CrossRef]
  15. Ju, Y.M.; Hsu, C.H.; Fang, L.S.; Lin, H.D.; Wu, J.H.; Han, C.C.; Chen, I.S.; Chiang, T.Y. Population structure and demographic history of Sicyopterus japonicus (Perciformes; Gobiidae) in Taiwan inferred from mitochondrial control region sequences. Genet. Mol. Res. 2013, 12, 4046–4059. [Google Scholar] [CrossRef]
  16. Guo, X.; Zhao, D.; Jung, D.; Li, Q.; Kong, L.F.; Ni, G.; Nakano, T.; Matsukuma, A.; Kim, S.; Park, C.; et al. Phylogeography of the rock shell Thais clavigera (Mollusca): Evidence for long-distance dispersal in the northwestern Pacific. PLoS ONE 2015, 10, e0129715. [Google Scholar] [CrossRef]
  17. Kojima, S.; Segawa, R.; Hayashi, I. Genetic differentiation among populations of the Japanese turban shell Turbo (Batillus) cornutus corresponding to warm currents. Mar. Ecol. Prog. Ser. 1997, 150, 149–155. [Google Scholar] [CrossRef]
  18. Chiu, Y.W.; Bor, H.; Tan, M.S.; Lin, H.D.; Jean, C.T. Phylogeography and genetic differentiation among populations of the moon turban snail Lunella granulata Gmelin, 1791 (Gastropoda: Turbinidae). Int. J. Mol. Sci. 2013, 14, 9062–9079. [Google Scholar] [CrossRef]
  19. Zhao, D.; Kong, L.F.; Sasaki, T.; Li, Q. Shell variations in the gastropod, Monodonta labio, in the North-western Pacific: The important role of temperature in the evolution process. J. Mar. Biolog. Assoc. UK 2019, 99, 1591–1599. [Google Scholar] [CrossRef]
  20. Zhao, D.; Li, Q.; Kong, L.; Yu, H. Cryptic diversity of marine gastropod Monodonta labio (Trochidae): Did the early Pleistocene glacial isolation and sea surface temperature gradient jointly drive diversification of sister species and/or subspecies in the Northwestern Pacific? Mar. Ecol. 2017, 38, e12443. [Google Scholar] [CrossRef]
  21. Yamazaki, D.; Miura, O.; Ikeda, M.; Kijima, A.; Van Tu, D.; Sasaki, T.; Chiba, S. Genetic diversification of intertidal gastropoda in an archipelago: The effects of islands, oceanic currents, and ecology. Mar. Biol. 2017, 164, 184. [Google Scholar] [CrossRef]
  22. Yamazaki, D.; Miura, O.; Uchida, S.; Ikeda, M.; Chiba, S. Comparative seascape genetics of co-distributed intertidal snails Monodonta spp. in the Japanese and Ryukyu archipelagoes. Mar. Ecol. Prog. Ser. 2021, 657, 135–146. [Google Scholar] [CrossRef]
  23. Folmer, O.; Black, M.; Hoeh, W.; Lutz, R.; Vrijenhoek, R. DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates. Mol. Mar. Biol. Biotechnol. 1994, 3, 294–297. [Google Scholar] [PubMed]
  24. Simon, C.; Franke, A.; Martin, A. The polymerase chain reaction: DNA extraction and amplification. Mol. Tech. Taxon. 1991, 57, 329–355. [Google Scholar] [CrossRef]
  25. Thompson, J.D.; Gibson, T.J.; Plewniak, F.; Jeanmougin, F.; Higgins, D.G. The CLUSTAL X windows interface: Flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res. 1997, 25, 4876–4882. [Google Scholar] [CrossRef]
  26. Librado, P.; Rozas, J. DnaSP v5: A software for comprehensive analysis of DNA polymorphism data. Bioinformatics 2009, 25, 1451–1452. [Google Scholar] [CrossRef]
  27. Jukes, T.H.; Cantor, C.R. Evolution of protein molecules. Mamm. Protein Metab. 1969, 3, 21–132. [Google Scholar]
  28. Excoffier, L.; Lischer, H.E.L. Arlequin suite ver 3.5: A new series of programs to perform population genetics analyses under Linux and Windows. Mol. Ecol. Resour. 2010, 10, 564–567. [Google Scholar] [CrossRef]
  29. Kumar, S.; Stecher, G.; Li, M.; Knyaz, C.; Tamura, K. MEGA X: Molecular evolutionary genetics analysis across computing platforms. Mol. Biol. Evol. 2018, 35, 1547. [Google Scholar] [CrossRef]
  30. Guindon, S.; Dufayard, J.F.; Lefort, V.; Anisimova, M.; Hordijk, W.; Gascuel, O. New algorithms and methods to estimate maximum-likelihood phylogenies: Assessing the performance of PhyML 3.0. Syst. Biol. 2010, 59, 307–321. [Google Scholar] [CrossRef]
  31. Huelsenbeck, J.P.; Ronquist, F. MRBAYES: Bayesian inference of phylogenetic trees. Bioinformatics 2001, 17, 754–755. [Google Scholar] [CrossRef]
  32. Posada, D. jModelTest: Phylogenetic model averaging. Mol. Biol. Evol. 2008, 25, 1253–1256. [Google Scholar] [CrossRef]
  33. Lefort, V.; Longueville, J.E.; Gascuel, O. SMS: Smart model selection in PhyML. Mol. Biol. Evol. 2017, 34, 2422–2424. [Google Scholar] [CrossRef]
  34. Leigh, J.W.; Bryant, D. POPART: Full-feature software for haplotype network construction. Methods Ecol. Evol. 2015, 6, 1110–1116. [Google Scholar] [CrossRef]
  35. Tajima, F. The effect of change in population size on DNA polymorphism. Genetics 1989, 123, 597–601. [Google Scholar] [CrossRef] [PubMed]
  36. Fu, Y.X. Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection. Genetics 1997, 147, 915–925. [Google Scholar] [CrossRef]
  37. Drummond, A.J.; Suchard, M.A.; Xie, D.; Rambaut, A. Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol. Biol. Evol. 2012, 29, 1969–1973. [Google Scholar] [CrossRef] [PubMed]
  38. Rambaut, A.; Suchard, M.A.; Xie, D.; Drummond, A.J. Tracer v1.6. Available online: http://beast.bio.ed.ac.uk/Tracer (accessed on 26 March 2021).
  39. Hellberg, M.E.; Vacquier, V.D. Rapid evolution of fertilization selectivity and lysin cDNA sequences in teguline gastropods. Mol. Biol. Evol. 1999, 16, 839–848. [Google Scholar] [CrossRef]
  40. Beerli, P.; Felsenstein, J. Maximum likelihood estimation of a migration matrix and effective population sizes in n subpopulations by using a coalescent approach. Proc. Natl. Acad. Sci. USA 2001, 98, 4563–4568. [Google Scholar] [CrossRef]
  41. Cornuet, J.M.; Pudlo, P.; Veyssier, J.; Dehne-Garcia, A.; Gautier, M.; Leblois, R.; Marin, J.M.; Estoup, A. DIYABC v2.0: A software to make approximate Bayesian computation inferences about population history using single nucleotide polymorphism, DNA sequence and microsatellite data. Bioinformatics 2014, 30, 1187–1189. [Google Scholar] [CrossRef]
  42. Cabrera, A.A.; Palsbøll, P.J. Inferring past demographic changes from contemporary genetic data: A simulation-based evaluation of the ABC methods implemented in diyabc. Mol. Ecol. Resour. 2017, 17, e94–e110. [Google Scholar] [CrossRef] [PubMed]
  43. Pons, O.; Petit, R.J. Measuring and testing genetic differentiation with ordered versus unordered alleles. Genetics 1996, 144, 1237–1245. [Google Scholar] [CrossRef] [PubMed]
  44. Benzie, J.A.H.; Ballment, E.; Forbes, A.T.; Demetriades, N.T.; Sugama, K.; Moria, S. Mitochondrial DNA variation in Indo-Pacific populations of the giant tiger prawn, Penaeus monodon. Mol. Ecol. 2002, 11, 2553–2569. [Google Scholar] [CrossRef] [PubMed]
  45. Kojima, Y. On the spawning of a top-shell, Monodonta labio (Linne). Venus (Jpn. J. Malacol.) 1962, 22, 200–203. [Google Scholar]
  46. Hu, Z.M.; Zhang, J.; Lopez-Bautista, J.; Duan, D.L. Asymmetric genetic exchange in the brown seaweed Sargassum fusiforme (Phaeophyceae) driven by oceanic currents. Mar. Biol. 2013, 160, 1407–1414. [Google Scholar] [CrossRef]
  47. Grant, W.A.S.; Bowen, B.W. Shallow population histories in deep evolutionary lineages of marine fishes: Insights from sardines and anchovies and lessons for conservation. J. Hered. 1998, 89, 415–426. [Google Scholar] [CrossRef]
  48. Chen, X.; Wang, J.J.; Ai, W.M.; Chen, H.; Lin, H.D. Phylogeography and genetic population structure of the spadenose shark (Scoliodon macrorhynchos) from the Chinese coast. Mitochondrial DNA Part A 2018, 29, 1100–1107. [Google Scholar] [CrossRef] [PubMed]
  49. Yi, M.R.; Hsu, K.C.; Wang, J.X.; Feng, B.; Lin, H.D.; Yan, Y.R. Genetic structure and diversity of the yellowbelly threadfin bream Nemipterus bathybius in the Northern South China Sea. Diversity 2021, 13, 324. [Google Scholar] [CrossRef]
  50. Gu, S.; Yan, Y.R.; Yi, M.R.; Luo, Z.S.; Wen, H.; Jiang, C.P.; Lin, H.D.; He, X.B. Genetic pattern and demographic history of cutlassfish (Trichiurus nanhaiensis) in South China Sea by the influence of Pleistocene climatic oscillations. Sci. Rep. 2022, 12, 14716. [Google Scholar] [CrossRef] [PubMed]
  51. Chen, P.C.; Shih, C.H.; Chu, T.J.; Wang, D.; Lee, Y.C.; Tzeng, T.D. Population structure and historical demography of the oriental river prawn (Macrobrachium nipponense) in Taiwan. PLoS ONE 2015, 10, e0145927. [Google Scholar] [CrossRef] [PubMed]
  52. Liu, S.Y.V.; Huang, I.H.; Liu, M.Y.; Lin, H.D.; Wang, F.Y.; Liao, T.Y. Genetic stock structure of Terapon jarbua in Taiwanese waters. Mar. Coast. Fish. 2015, 7, 464–473. [Google Scholar] [CrossRef]
  53. Ding, S.; Mishra, M.; Wu, H.; Liang, S.; Miyamoto, M.M. Characterization of hybridization within a secondary contact region of the inshore fish, Bostrychus sinensis, in the East China Sea. Heredity 2018, 120, 51–62. [Google Scholar] [CrossRef] [PubMed]
  54. Hirase, S.; Kanno, M.; Ikeda, M.; Kijima, A. Evidence of the restricted gene flow within a small spatial scale in the Japanese common intertidal goby Chaenogobius annularis. Mar. Ecol. 2012, 33, 481–489. [Google Scholar] [CrossRef]
Diversity 15 01021 g001
Figure 1. Map showing the 13 sampling localities of Monodonta labio (Lineage I) and M. confusa (Lineage II) examined in this study. Collection sites (circule) correspond to locations given in the text and Table 1. The frequencies of the haplogroups (Figure 3) in each population are displayed on the map. The areas shaded with diagonal lines in the figure represent the regions that were exposed at sea level during the glacial period, indicating the coastline at that time.
Figure 1. Map showing the 13 sampling localities of Monodonta labio (Lineage I) and M. confusa (Lineage II) examined in this study. Collection sites (circule) correspond to locations given in the text and Table 1. The frequencies of the haplogroups (Figure 3) in each population are displayed on the map. The areas shaded with diagonal lines in the figure represent the regions that were exposed at sea level during the glacial period, indicating the coastline at that time.
Diversity 15 01021 g001
Diversity 15 01021 g002
Figure 2. Schematic representation of five demographic scenarios for Monodonta labio and M. confusa tested by approximate Bayesian computation (ABC) Time (t1–t3) and effective population size (N–N5) are not scaled.
Figure 2. Schematic representation of five demographic scenarios for Monodonta labio and M. confusa tested by approximate Bayesian computation (ABC) Time (t1–t3) and effective population size (N–N5) are not scaled.
Diversity 15 01021 g002
Diversity 15 01021 g004
Figure 4. The median-joining network of 94 haplotypes of Monodonta labio and M. confusa. Each circle represents a haplotype, and the circle size is proportional to the haplotype frequency. The bars across the branches are the number of mutational steps. Black circles represent missing haplotypes. The sub-figures (I, II, A, B, C) referred to in this figure correspond to the clusters in Figure 3.
Figure 4. The median-joining network of 94 haplotypes of Monodonta labio and M. confusa. Each circle represents a haplotype, and the circle size is proportional to the haplotype frequency. The bars across the branches are the number of mutational steps. Black circles represent missing haplotypes. The sub-figures (I, II, A, B, C) referred to in this figure correspond to the clusters in Figure 3.
Diversity 15 01021 g004
Diversity 15 01021 g005
Figure 5. Bayesian skyline plot of the effective population sizes over time for Monodonta labio (a) and M. confusa (b) The units of the X axis are millions of years, and the units of the Y axis are Ne (the log-transformed product of effective population size and generation length in years).
Figure 5. Bayesian skyline plot of the effective population sizes over time for Monodonta labio (a) and M. confusa (b) The units of the X axis are millions of years, and the units of the Y axis are Ne (the log-transformed product of effective population size and generation length in years).
Diversity 15 01021 g005
Table 1. List of sampling locations and the sample sizes, the number of haplotypes (N), haplotype diversities (h), and nucleotide diversities of the mtDNA COI + 16S region sequences for each location.
Table 1. List of sampling locations and the sample sizes, the number of haplotypes (N), haplotype diversities (h), and nucleotide diversities of the mtDNA COI + 16S region sequences for each location.
LocationsSampling LocationCodeCoordinatesNNhhθπθω
Monodonta labio85610.9870.01310.0151
Taiwan Island 22200.9910.0.04480.0376
AudiAD121.92 25.08551.0000.00320.0035
MiaoliWW120.73 24.61661.0000.00410.0052
TainanCG120.10 23.1211100.9820.00190.0028
Penghu 970.9170.00170.0027
PenghuPH119.58 23.58970.9170.00170.0027
Coast of mainland China 15120.9620.01790.0387
KinmenKM118.35 24.43870.9640.03150.0480
MatsuMZ119.98 26.22771.0000.00270.0030
Japan 15100.8950.00150.0029
IneYG135.29 35.70970.9170.00120.0019
KagoshimaSI130.46 32.51650.9330.00180.0022
Ryukyu 24170.9640.00370.0056
OkinawaOK127.79 24.251090.9780.00350.0044
IshigakiSY124.11 24.361480.9010.00320.0035
Monodonta confusa41330.9830.00820.0115
Taiwan Island 12121.0000.0.00860.0096
AudiAD121.92 25.08771.0000.00930.0097
TaitungTD121.19 22.82551.0000.00890.0086
Penghu 551.0000.00930.0093
PenghuPH119.58 23.58551.0000.00930.0093
Coast of mainland China 760.9520.00500.0067
MatsuMZ119.98 26.22760.9520.00500.0067
Japan 13100.9230.00370.0050
InubosakiCF140.87 25.701070.8670.00310.0034
IneYG135.29 35.70331.0000.00490.0049
Korea GY 441.0000.00420.0044
Total 126940.9920.08250.0441
N: sample size π: nucleotide diversity, Nh: number of haplotypes h: haplotype diversity.
Table 2. (A). Matrix of pairwise FST among thirteen populations based on mitochondrial COI + 16S gene (below diagonal) and p values (above diagonal) in Monodonta labio. (B). Matrix of pairwise FST among thirteen populations based on mitochondrial COI + 16S gene (below diagonal) and p values (above diagonal) in Monodonta confusa.
Table 2. (A). Matrix of pairwise FST among thirteen populations based on mitochondrial COI + 16S gene (below diagonal) and p values (above diagonal) in Monodonta labio. (B). Matrix of pairwise FST among thirteen populations based on mitochondrial COI + 16S gene (below diagonal) and p values (above diagonal) in Monodonta confusa.
(A)
ADWWCGPHKMMZYGSIOKSY
AD 0.0000.0000.0000.0000.0000.0000.2700.0000.000
WW0.968 0.4140.0000.8910.7740.0000.0000.0000.000
CG0.9800.000 0.0000.2520.5130.0000.0000.0000.000
PH0.9810.4550.594 0.0000.0000.0000.0000.0000.000
KM0.799−0.0410.0340.145 0.9810.0000.0000.0000.000
MZ0.975−0.047−0.0080.519−0.030 0.0000.0000.0000.000
YG0.1520.9790.9860.9870.8490.984 0.5580.0000.000
SI0.0290.9740.9830.9850.8160.980−0.015 0.0000.000
OK0.9730.7950.8490.8500.4830.8190.9810.977 0.000
SY0.9750.7970.8430.8440.5220.8170.9810.9780.154
(B)
ADTDPHMZCFYGGJ
AD 0.9720.3960.9360.0000.0360.027
TD−0.133 0.5940.9900.0000.0090.072
PH 0.013−0.042 0.3420.0000.0000.009
MZ −0.137−0.138−0.002 0.0000.0360.036
CF 0.3950.4420.6050.459 0.0360.279
YG 0.3320.3510.0530.3560.193 0.576
GJ0.2600.2840.5280.3060.018−0.049　
Table 3. Analysis of molecular variance (AMOVA) for Monodonta labio and M. confusa populations based on mitochondrial COI + 16S genes.
Table 3. Analysis of molecular variance (AMOVA) for Monodonta labio and M. confusa populations based on mitochondrial COI + 16S genes.
Source of VariationPercentage Variation p
Monodonta labio
1. Scenario I: two groups (WW, CG, PH, KM, MZ, YG) (AD, OK, SY, SI)
Among groups−5.50FCT = −0.0550.461
Among populations within groups72.41FSC = 0.6860.000
Among individuals within populations33.10FST = 0.6690.000
2. Scenario II: three groups (AD, WW, CG, PH, KM, MZ) (OK, SY) (YG, SI)
Among groups68.40FCT = 0.6830.016
Among populations within groups24.36FSC = 0.7700.000
Among individuals within populations7.24FST = 0.9270.000
3. Scenario III: five groups (AD, WW, CG, PH) (KM, MZ) (OK, SY) (YG, SI)
Among groups59.93FCT = 0.5990.000
Among populations within groups31.74FSC = 0.7920.000
Among individuals within populations8.31FST = 0.9160.000
Monodonta confusa
1. Scenario I: two groups (PH, MZ, YG, GJ) (AD, TD, CF)
Among groups−10.46FCT = −0.1041.000
Among populations within groups35.55FSC = 0.3210.000
Among individuals within populations74.90FST = 0.2500.000
2. Scenario II: three groups (AD, TD, PH, MZ) (CF, YG, GJ)
Among groups43.25FCT = 0.4320.032
Among populations within groups−2.04FSC = −0.0350.686
Among individuals within populations58.78FST =0.4120.000
3. Scenario III: five groups (AD, TD, PH) (MZ) (CF, YG) (GJ)
Among groups33.51FCT = 0.3350.000
Among populations within groups−0.27FSC = −0.0040.546
Among individuals within populations66.75FST = 0.3320.000
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Chiu, Y.-W.; Bor, H.; Wu, J.-X.; Shieh, B.-S.; Lin, H.-D. Population Structure and Phylogeography of Marine Gastropods Monodonta labio and M. confusa (Trochidae) along the Northwestern Pacific Coast. Diversity 2023, 15, 1021. https://doi.org/10.3390/d15091021

AMA Style

Chiu Y-W, Bor H, Wu J-X, Shieh B-S, Lin H-D. Population Structure and Phylogeography of Marine Gastropods Monodonta labio and M. confusa (Trochidae) along the Northwestern Pacific Coast. Diversity. 2023; 15(9):1021. https://doi.org/10.3390/d15091021

Chicago/Turabian Style

Chiu, Yuh-Wen, Hor Bor, Jin-Xian Wu, Bao-Sen Shieh, and Hung-Du Lin. 2023. "Population Structure and Phylogeography of Marine Gastropods Monodonta labio and M. confusa (Trochidae) along the Northwestern Pacific Coast" Diversity 15, no. 9: 1021. https://doi.org/10.3390/d15091021

APA Style

Chiu, Y. -W., Bor, H., Wu, J. -X., Shieh, B. -S., & Lin, H. -D. (2023). Population Structure and Phylogeography of Marine Gastropods Monodonta labio and M. confusa (Trochidae) along the Northwestern Pacific Coast. Diversity, 15(9), 1021. https://doi.org/10.3390/d15091021

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop