Next Article in Journal
Assessment of the Centrality of the Cruise Ship Navigation Networks in Southern Europe through the Application of Social Network Analysis
Next Article in Special Issue
Seasonal Variations in the Fish Species Composition and Community Structures on the Eastern Coast of South Korea
Previous Article in Journal
Improved Robust High-Degree Cubature Kalman Filter Based on Novel Cubature Formula and Maximum Correntropy Criterion with Application to Surface Target Tracking
Previous Article in Special Issue
Diet of Adult Sardine Sardina pilchardus in the Gulf of Trieste, Northern Adriatic Sea
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genetic Variations and Expansion of the Blue Swimmer Crab (Portunus pelagicus) in Southeast Asia

1
Institute of Fisheries Science, National Taiwan University, Taipei City 10617, Taiwan
2
Department of Leisure & Tourism Management, Shu-Te University, Kaohsiung 82445, Taiwan
3
Fisheries College, Jimei University, Xiamen 361021, China
4
Fujian Provincial Key Laboratory of Marine Fishery Resources and Eco-Environment, Xiamen 361021, China
5
Third Institute of Oceanography, Ministry of Natural Resources, Xiamen 361005, China
6
Department of Hospitality and Baking Management, Shu-Te University, Kaohsiung 82445, Taiwan
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
J. Mar. Sci. Eng. 2022, 10(8), 1071; https://doi.org/10.3390/jmse10081071
Submission received: 25 May 2022 / Revised: 1 August 2022 / Accepted: 2 August 2022 / Published: 4 August 2022
(This article belongs to the Special Issue Marine Fisheries Management)

Abstract

:
Blue swimmer crabs (Portunus pelagicus) primarily inhabit the coastal regions of Southeast Asia. This study aimed to examine their population genetic structure, phylogeography, and historical demography by applying partial sequences to the control region in the mitochondrial DNA. The analysis focused on 6 populations comprising 166 individual samples collected from Taiwan (Kezailiao Kaohsiung (KK) and Penghu (PH)), mainland China (Xiamen (XM) and Hong Kong (HK)), Vietnam (Hanoi (VN)), and Singapore (SGP). The estimated nucleotide diversity (π) for all of the samples was 0.062, with values ranging between 0.018 (PH) and 0.045 (HK). Our study collected 116 haplotypes and classified them into lineages A and B. The XM, HK, VN, and SGP populations made up lineage A, and the PH and KK populations comprised lineage B. For each lineage, the results highlighted indications of demographic expansion. Lineage A can be traced back to the second-to-last interglacial period, while lineage B possibly originated in the last glacial period. There were significant differences in the FST values among the six populations, except for the pairs HK–XM and PH–KK. In conclusion, the phylogeography and genetic structure of P. pelagicus in Southeast Asia were potentially affected by the Pleistocene glacial cycles and human introduction. These findings further highlight the possible dispersal routes of P. pelagicus throughout Southeast Asia.

1. Introduction

Blue swimmer crabs (Portunus pelagicus (Linnaeus)), are found in the Indo-Pacific, Southeast and East Asia (i.e., Japan, Korea, China, Taiwan, and the Bay of Bengal), and the west, north, and east of Australia [1,2]. They tend to be dispersed across several kilometres by ocean currents, are generally not limited by demographic boundaries, and have high population exchange levels [3,4]. The species is particularly valued in the commercial fisheries of Southeast Asia, and a key factor in the successful and sustainable management of fishery resources is population genetic structure [5]. Research on population genetics [6] and historical demography [7,8] has offered valuable data on catches and age composition [9]. Studies on the population genetic structure of P. pelagicus have demonstrated that highly dispersed organisms tend to show low genetic distinction across a large geographic scale [10,11,12]. Moreover, the gene flow rate is high among marine populations [13]. However, given their high dispersal potential and the lack of obvious barriers to dispersal in the ocean, marine organisms are more genetically structured than expected, in contrast with the findings of previous studies on genetic structure [14,15,16]. However, habitats, local oceanic conditions, and historical events impede the dispersal of marine organisms [17], creating the potential for genetic distinction [18].
Marine organism research has applied a wide variety of methodologies to gain an understanding of population genetic structures [18]. Some studies have analysed the distribution and abundance of life-history stages, while others have explored marking and tagging, morphological characteristics, allozymes, or DNA markers [1,2]. A multipronged approach that integrates the outcomes of several methods could maximise the probability of accurately defining stocks [19]. Mitochondrial (mt) DNA is a key component in population genetic studies given its rapid evolution rate, lack of recombination, and maternal inheritance [14,15,20]; mtDNA—specifically, its control region—is an ideal marker of geographical patterns in genetic variations within and across crab populations [20,21,22,23,24,25,26,27,28,29,30,31]. However, research findings on the phylogeography and genetic structure of P. pelagicus in Southeast Asia have remained inconsistent.
Geological evidence indicates that land bridges connected the island of Taiwan to the Asian continent 2–3 times during the Pleistocene [20]. During the Late Pleistocene glaciation, sea levels in the area were about 140 m lower than at present [20]. As a result, the floor of the Taiwan Strait was exposed as a broad land bridge that was crossed by mainland fauna until the beginning of the Holocene 10,000 years ago [20].The changing sea levels and temperature during the glacial period critically altered the marine environment and, with it, the distribution, connectivity, and genetic structure of marine organisms [32]. Paleoclimate changes influenced the modern biodiversity patterns of the continent and these islands [33]. The Pleistocene glacial period resulted in fluctuations in East Asia’s sea levels which, in turn, repeatedly altered the connections between these islands or with the mainland [34,35]. More specifically, the Taiwan Strait was almost separated from the Pacific Ocean, the South China Sea was completely exposed, and the major islands of Taiwan were connected to the continents [31]. Studies have reported that diverse and unique biota emerging from complex geological and climatic factors contributed to ancient substitution and diffusion events [36]. After the end of the glacial period, the islands were once again isolated
The abovementioned historical events critically influenced the genetic isolation, distribution, and structure of East Asian organisms [37]. Other significant factors included ocean currents [38] and past and current land barriers [39,40]. A genetic study using gene markers in the mtDNA cytochrome oxidase (COI) to determine the genetic structure of P. pelagicus around the South China Sea highlighted historical vicariance as an influencing factor [41]. The mtDNA D-loop has high variability, a high level of polymorphism in other crustacean species, and can be used to distinguish between populations [42,43,44]. Approximately 0.01–0.11 million years ago (Mya), the South China Sea was 100–120 m lower owing to global glaciation, and the Pacific and Indian oceans were separated by the Sundaland—an exposed land ridge [11,12,45]. This study explores the impact of these areas on the population structure of the subsamples of species distribution. This study aims to explore the population genetic structure of blue swimmer crabs in Southeast Asia by conducting sequence analyses on the mtDNA control region in the species.

2. Materials and Methods

2.1. Sample Collection

This study is based on 166 individual specimens of P. pelagicus collected from 6 regions: Kezailiao Kaohsiung (KK) and Penghu (PH) in Taiwan, Xiamen (XM) and Hong Kong (HK) in mainland China, Hanoi (VN) in Vietnam, and Singapore (SGP). The specimens were collected using commercial crab traps or gillnets (Table 1) and immediately stored at −75 °C before DNA was extracted. This study was conducted under permits granted by wildlife conservation authorities in Vietnam, Singapore, China, and Taiwan (Figure 1).

2.2. DNA Extraction, PCR, and Sequencing

Genomic DNA was extracted from the pereopod muscle using a QIAamp DNA Mini Kit [46]. A fragment of the mtDNA control region was amplified and sequenced using YEN-F (5′- GCA AAT ACA CGC AAT AAC TCT CAT AC -3′) and YEN-B (5′- TGT AAA TCC GTT ACG AAT AAT ATA GG -3′) primers. Thermal cycling was performed using a GeneAmp 2400 thermal cycler (PerkinElmer, Norwalk, CT, USA). The PCR conditions were as follows: 39 cycles of denaturation at 95 °C for 50 s, annealing at 50 °C for 1 min, and an extension at 72 °C for 1.5 min. The first cycle included an initial denaturation step at 95 °C for 5 min, and the last cycle had a final extension at 72 °C for 10 min. The PCR product was separated on 1.5% agarose gels using electrophoresis. Purification was conducted using the Gene Clean II Kit (Bio101, Vista, CA, USA). Sequencing was peformed using the ABI 377 DNA sequencer (Applied Biosystems, Inc., Foster City, CA, USA).

2.3. Analyses of Sequences

MegAlign (DNASTAR, Lasergene, Madison, WI, USA) was used to align all of the sequences, and DnaSP version 5.00 was used to estimate the numbers of variable and parsimony-informative sites, base composition, and haplotype and nucleotide diversities [47,48,49,50].
The neighbour-joining (NJ) and maximum likelihood (ML) methods were employed to conduct phylogeographic analyses of the control region using MEGA X [31,49,50]. Bootstrap analyses were conducted on 1000 replicates to evaluate the phylogenetic relationships of all haplotypes. MEGA was used to determine the optimal substitution model. The median-joining method with Network version 4.6.1.3 (www.fluxus-engineering.com [accessed 1 February 2021]) was employed to construct a network of haplotypes [50]. Insight into the expansion of the historical demography was derived by using DnaSP to examine the frequency distributions of pairwise differences between sequences (mismatched distribution) [31]. Rough dates of population expansion were estimated using the formula T = τ/2 u [7,51], where T = time since expansion, τ = expansion time, and 2 u = μ (mutation rate) × generation time × number of bases sequenced. The average divergence rate was 1.72–2.91% per million years, with a generation time of one year.
Pairwise FST statistics were estimated using ARLEQUIN version 3.5 (Excoffier and Lischer 2010) to examine genetic differentiation between sites [51,52]. Next, on the basis of the FST values, a dendrogram of the six sampling sites was constructed using the unweighted pair–group method with arithmetic means (UPGMA). Population structure was assessed by conducting an analysis of molecular variance (AMOVA [6]) using ARLEQUIN. The UPGMA tree highlighted various groupings for the six populations, of which one grouping reported the maximal value for ΦCT. The groupings significantly differed from a random organisation of similar groupings. Furthermore, they indicated geographic subdivisions that were the most probable [53,54]. The statistical results were tested for significance using 10,000 random permutations.
Population expansion or deviation from neutrality was analysed using three methods. Furthermore, to determine whether population expansion occurred in the past, deviations from neutrality were examined by estimating Tajima’s D [55]. In addition, evidence of population expansion was analysed by conducting Fu’s Fs test [8] using DnaSP. A mismatch analysis was also performed using DnaSP to check for population expansion and the frequency distributions in the nucleotide differences as a function of frequency [56].

3. Results

The sequencing analyses of the 166 specimens revealed 92 variable sites and 199 parsimony-informative sites, thus detecting a total of 116 unique haplotypes (S1 Dataset). The D-loop fragment was estimated at 715 bp. All of the sequences were deposited in GenBank (accession numbers MZ512410–MZ512575), and no gap was detected. The nucleotide composition frequency indicated an AT bias (G + C content: 46.0%).
Table 1 presents the estimates for the haplotype (h) and nucleotide (π) diversities of all populations. The average value for haplotype diversity was between 0.949 and 0.995, with values ranging between 0.949 (VN) and 0.995 (HK), while that for nucleotide diversity was 0.064, with values ranging from 0.020 (KK) to 0.045 (HK).
The study’s data were best explained by the T92 model, which was subsequently applied to NJ and ML reconstructions and AMOVA analyses. Figure 1 illustrates the phylogenetic tree for all haplotypes. The tree and the network highlight two distinct lineages (A and B) (Figure 1 and Figure 2). The individual specimens in KK belonged only to lineage B; those in HK, XM, VN, and SGP were only in lineage A; and specimens in PH were in either lineage A or B. For lineages A and B, the estimated values for haplotype diversity (h) were 0.993 and 0.990, respectively, while those for nucleotide diversity (π) were 0.040 and 0.031, respectively (Table 1). The τ values for lineages A and B were 19.439/2µ and 11.712/2µ generations, respectively.
Table 2 reports the FST results. The tests showed significant values for all six populations, except for the relationships between HK and XM, and between PH and KK. Figure 3 illustrates the UPGMA tree for the six sampling areas. According to the intergroup variance component test results, five of the six designated grouping schemes (G1, G2, G3, G4, and G5) were statistically significant. The selection of the G2, G3, and G4 grouping schemes as the best grouping schemes was based on the statistics of ΦCT values (Figure 4).
The Tajima’s D values were not significant for either lineage, for the six individual populations, or for the total population (Table 1). The Fu’s Fs test results were significant for all populations (Table 1), both lineages, and for the total population. Statistically, a unimodal distribution was obtained from lineage A or B when measured by squared deviation (p > 0.05). The distribution did not significantly differ from that estimated by the growth expansion model (Figure 5).
Table 3 shows the AMOVA results, and for all six populations, they reported a significant FST value of 0.41978. Specifically, there was significant heterogeneity in at least one of the pairwise populations. ΦCT values were significant for all groupings, of which group 2 reported the highest ΦCT value of 0.46573. These results support the conclusion that the six populations can be divided into two main groups: the first group comprises the HK, XM, VN, and SGP populations, while the second group comprises the PH and KK populations. Furthermore, the findings suggest the possibility of additional genetic discontinuity between populations, given the significant ΦCT values for different groupings.

4. Discussion

This study analysed two lineages (A and B) of blue swimmer crabs (P. pelagicus) inhabiting Southeast Asia. The findings are consistent with those of previous research [57,58,59,60]. Greater nucleotide and haplotype diversities in populations with ancestral genotypes are an outcome of the long-term accumulation of mutations [61,62,63]. The study reported significantly higher nucleotide and haplotype diversities in lineage A (π = 0.044; h = 0.992) than in lineage B (π = 0.031; h = 0.990), indicating that lineage A was older than lineage B. Population age was estimated on the basis of an average mutation rate of 2.3%/My. The time of expansion was 597,718 and 360,125 years for lineages A and B, respectively. It appeared that lineage B had a smaller original population prior to an expansion or bottleneck given its steeper peak in the mismatch distributions (Figure 5) [7]. Furthermore, the pairwise distribution was more clearly displaced to the right of the distribution pattern, suggesting that lineage B’s expansion could have been more recent [58,59,60].
As swimmer crabs progress through their life stages, they tend to migrate from inshore to offshore. However, given the limited distance of migration [38], larval dispersal remains the primary source of gene flow, with ocean currents playing an important role. This explains why the FST values between the KK and PH populations were at least twice as high as those obtained among the other populations. The mixing of populations and reduced genetic variation can also be attributed to disappearing habitats [17]. Since the Last Glacial Maximum, there has been a gradual expansion in blue swimmer crabs, depending on the sea level of the East China Sea.
While this study found two lineages (Figure 2 and Figure 3), samples from the six populations did not always belong to both lineages. The XM, HK, VN, and SGP populations were only in lineage A, whereas the PH and KK populations comprised lineage B. Geographical conditions such as coastal currents traveling from mainland China to the South China Sea, the Chinese coastal stream flowing along Fujian to Vietnam [64], and the thermohaline circulation in Singapore’s waters played a crucial role in the population structure. However, there was an isolated group in lineage B given the lack of obvious interactions between the KK samples in Taiwan and the coastal current in China for lineage A. P. pelagicus has an extended planktonic larval stage with potentially high larval dispersal, causing extensive gene flow between conspecific samples within a geographic mesoscale of 10–100 km [65]. In China, however, the population genetic structure had a limited occurrence and high gene flow along the distribution areas [41]. In Malaysia, microsatellite analyses have indicated low genetic differentiation between P. pelagicus populations [11], which is consistent with findings from Moreton Bay, Queensland, where the species was found to travel as far as 20 km in one day [66,67]. This study found that blue swimmer crabs tend to be mixed within a mesoscale region (e.g., the Java Sea), but differ genetically at larger scales [68]. This finding is consistent with those of previous studies that reported that adult P. pelagicus travel within a geographic mesoscale on a daily basis, making them vagile [69].
Past and present habitat conformations have considerably influenced the genetic patterns of P. pelagicus. The exposed land mass of China and Southeast Asia contributed to lineage A, which has an older evolutionary age than lineage B. During the Last Glacial Maximum, the species in lineage A may have dispersed along the coast from China to Southeast Asia. As the sea levels declined, they first colonised the HK area, and subsequently spread to the SGP coasts with the rise in sea levels and the withdrawal of shorelines. The rising sea immersed the coastline, giving way to the PH and KK populations in lineage B. The warm currents of the South China Sea possibly dispersed lineage B to Penghu’s waters.
The differing ages of the two lineages can be partially ascribed to the varying dispersal routes and times of origin, approximately one million years ago [1,41,43]. Lineage A—specifically the XM, HK, VN, PH, and SGP populations—carried the most prevalent allele. In other words, these populations likely share a common ancestor. For populations with ancestral genotypes, the long-term accumulation of mutations results in higher nucleotide and haplotype diversity [44,45,46]. The HK population may be the oldest population, because it showed significantly higher nucleotide (π = 0.046) and haplotype (h = 0.995) diversity compared to the other populations (Table 1). The KK population, on the other hand, may be the youngest. Our study also highlights that P. pelagicus originated in mainland China [16].
From the discussions thus far, it appears that the dispersal route for lineage B originated from the south of China (XM) to the north of areas at the same latitude (PH) and further toward southern Taiwan (KK). The SGP and VN populations in lineage A are to the extreme south of the sampling areas. The XM population in lineage B showed the lowest nucleotide and haplotype diversity. Furthermore, the relationship between the HK and SGP populations was not as close as that between the XM and HK populations (Figure 2). Thus, it is possible that lineage A had two distinct dispersal routes: one from the south of China to the south of the low-latitude areas (VN and SGP) and further north of Xiamen (XM), and the other one the south of China to the north of areas at the same latitude (PH and KK). Both lineages demonstrated demographic expansions, with lineage A starting in the second-to-latest interglacial period (597,718 years ago) and lineage B in the Last Glacial Maximum (360,125 years ago).
The distribution limit of lineage B probably extended to the Indian Ocean through the Strait of Malacca, given that the Sunda Shelf was blocked during the Last Glacial Maximum [17]. Thus, P. pelagicus found in the Indian Ocean possibly belong to another species, as indicated by Lai et al. [1]. More specifically, there was a sharp genetic break on the west side of the strait. Thus, the distribution of lineage B was limited to the South China Sea, and did not extend to the West Pacific (lineage A).
Studies have highlighted chaotic mismatch distributions only for populations in a long-term stable demographic equilibrium. However, recent rapid expansions or bottlenecks have resulted in unimodal (approximately Poisson) profiles with a steeper wave, indicating a smaller initial population prior to expansion [1,41]. This study revealed a steeper wave for lineage B (Figure 4).

5. Conclusions

This study found a strong genetic structure in Southeast Asia’s blue swimmer crabs. It identified two main lineages (A and B), of which lineage A was the older. Both lineages had experienced population expansion since the Pleistocene glacial cycles in Southeast Asia (approximately 597,718 and 360,125 years ago for lineages A and B, respectively). Our research also examined the possible dispersal routes of the species, the findings of which indicated that lineage B could have experienced expansion more recently than lineage A. Finally, land bridges formed during Pleistocene glacial maxima may have influenced the phylogeography and genetic structure of the species.

Author Contributions

Conceptualisation, Y.-M.L., C.-H.S., P.-C.C. and T.-D.T.; methodology, Y.-M.L., C.-H.S., P.-C.C. and T.-D.T.; software, Y.-M.L., C.-H.S. and P.-C.C.; validation, Y.-M.L., C.-H.S., P.-C.C., W.-C.K., Y.-C.L. and T.-D.T.; formal analysis, Y.-M.L., C.-H.S., P.-C.C., W.-C.K., Y.-C.L. and T.-D.T.; investigation, Y.-M.L., C.-H.S., P.-C.C., W.-C.K., Y.-C.L. and T.-D.T.; resources, Y.-M.L., C.-H.S. and P.-C.C.; data curation, Y.-M.L., C.-H.S. and P.-C.C.; writing—original draft preparation, Y.-M.L., C.-H.S. and P.-C.C.; writing—review and editing, Y.-M.L., C.-H.S., Y.-S.H. and P.-C.C.; visualisation, Y.-M.L., C.-H.S. and P.-C.C.; supervision, Y.-M.L., C.-H.S., P.-C.C. and T.-D.T.; project administration, Y.-M.L., C.-H.S., P.-C.C. and T.-D.T.; funding acquisition, Y.-M.L., C.-H.S., P.-C.C. and T.-D.T. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Science and Technology Council, grant number MOST 110-2121-M-366-001. The funders had no role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We thank the Global Institute for Green Tourism, University of California, Berkeley, for supporting the cruises of the biological survey. We would also like to thank the anonymous reviewers, whose useful suggestions were incorporated into the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lai, J.C.Y.; Ng, P.K.L.; Davie, P.J.F. A revision of the Portunus pelagicus (Linnaeus, 1758) species complex (Crustacea: Brachyura: Portunidae), with the recognition of four species. Raffles Bull. Zool. 2010, 58, 199–237. [Google Scholar]
  2. Huang, Y.H.; Shih, H.T. Diversity in the Taiwanese swimming crabs (Crustacea: Brachyura: Portunidae) estimated through DNA barcodes, with descriptions of 14 new records. Zool Stud. 2021, 60, 60. [Google Scholar]
  3. Becker, B.J.; Levin, L.A.; Fodrie, F.J.; McMillan, P.A. Complex larval connectivity patterns among marine invertebrate populations. Proc. Natl. Acad. Sci. USA 2007, 104, 3267–3272. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Zhou, S.; Smith, A.D.; Knudsen, E.E. Ending overfishing while catching more fish. Fish Fish. 2015, 16, 716–722. [Google Scholar] [CrossRef]
  5. Hillis, D.M.; Moritz, C.; Mable, B.K. Molecular Systematics, 1st ed.; Sinauer Associates Publishers Inc.: Sunderland, MA, USA, 1996. [Google Scholar]
  6. Excoffier, L.; Smouse, P.E.; Quattro, J.M. Analysis of molecular variance inferred from metric distances among DNA haplotypes: Application to human mitochondrial DNA restriction data. Genetics 1992, 131, 479–491. [Google Scholar] [CrossRef] [PubMed]
  7. Rogers, A.R.; Harpending, H. Population growth makes waves in the distribution of pairwise genetic differences. Mol. Biol. Evol. 1992, 9, 552–569. [Google Scholar]
  8. Fu, Y.X. Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection. Genetics 1997, 147, 915–925. [Google Scholar] [CrossRef]
  9. Pauly, D.; Christensen, V.; Guénette, S.; Pitcher, T.J.; Sumaila, U.R.; Walter, C.J.; Watson, R.; Zeller, D. Towards Sustainability in World Fisheries. Nature 2002, 418, 689–695. [Google Scholar] [CrossRef]
  10. Hellberg, M.E. Dependence of gene flow on geographic distance in two solitary corals with different larval dispersal capabilities. Evolution 1996, 50, 1167–1175. [Google Scholar] [CrossRef]
  11. Chai, C.J.; Bin Esa, Y.; Ismail, S.; Kamarudin, M.S. Population Structure of the Blue Swimmer Crab Portunus pelagicus in Coastal Areas of Malaysia Inferred from Microsatellites. Zool Stud. 2017, 16, 56.e26. [Google Scholar]
  12. Madduppa, H.; Martaulina, R.; Zairion, Z.; Renjani, R.M.; Kawaroe, M.; Anggraini, N.P.; Subhan, B.; Verawati, I.; Sani, L.M. Genetic population subdivision of the blue swimming crab (Portunus pelagicus) across Indonesia inferred from mitochondrial DNA: Implication to sustainable fishery. PLoS ONE 2021, 16, e0240951. [Google Scholar] [CrossRef]
  13. Haye, P.A.; Segovia, N.I.; Munoz-Herrera, N.C.; Galvez, F.E.; Martinez, A.; Meynard, A.; Pardo-Gandarillas, M.C.; Poulin, E.; Faugeron, S. Phylogeographic structure in benthic marine invertebrates of the southeast Pacific coast of Chile with differing dispersal potential. PLoS ONE 2014, 9, e88613. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Palumbi, S.R. Molecular biogeography of the Pacific. Coral Reefs 1997, 16, S47–S52. [Google Scholar] [CrossRef]
  15. Benzie, J.A.H. Genetic structure of coral reef organisms: Eghosts of dispersal past. Am. Zool. 1999, 39, 131–145. [Google Scholar] [CrossRef] [Green Version]
  16. Briggs, J.S. Coincident biogeographic patterns: Indo-West Pacific Ocean. Evolution 1999, 53, 326–335. [Google Scholar] [CrossRef]
  17. Benzie, J.A.H.; Williams, S.T. Genetic structure of giant clam (Tridacna maxima) populations in the West Pacific is not consistent with dispersal by present-day ocean currents. Evolution 1997, 51, 768–783. [Google Scholar]
  18. Palumbi, S.R. Genetic divergence, reproductive isolation, and marine speciation. Annu. Rev. Ecol. Syst. 1994, 25, 547–572. [Google Scholar] [CrossRef]
  19. Pawson, M.G.; Jennings, S. A critique of methods for stock identification in marine captured fisheries. Fish. Res. 1996, 25, 203–217. [Google Scholar] [CrossRef]
  20. Chen, P.C.; Shih, C.H.; Chu, T.J.; Lee, Y.C.; Tzeng, T.D. Phylogeography and genetic structure of the oriental river prawn Macrobrachium nipponense (Crustacea: Decapoda: Palaemonidae) in East Asia. PLoS ONE 2017, 12, e0173490. [Google Scholar] [CrossRef] [Green Version]
  21. Held, C.; Koenemann, S.; Schubart, C.D. Phylogeography and Population Genetics in Crustacea; CRC Publishers Inc.: Boca Raton, FL, USA, 2016. [Google Scholar]
  22. Yang, C.P.; Li, H.X.; Li, L.; Yan, Y. Occurrence and effects of the rhizocephalan parasite Diplothylacus sinensis (Cirripedia: Rhizocephala: Thomsoniidae) in the swimming crab Portunus sanguinolentus (Decapoda: Portunidae) in Honghai Bay, South China Sea. J. Crustacean Biol. 2014, 34, 573–580. [Google Scholar] [CrossRef] [Green Version]
  23. Tzeng, T.D. Stock identification of sword prawn (Parapenaeopsis hardwickii) in the East China Sea and Taiwan Strait inferred by morphometric variation. Fish Sci. 2004, 70, 758–764. [Google Scholar] [CrossRef]
  24. Tzeng, T.D. Population structure of the sword prawn (Parapenaeopsis hardwickii) (Decapoda: Penaeidae) in the East China Sea and waters adjacent to Taiwan inferred from the mitochondrial control region. Zool. Stud. 2007, 46, 561–568. [Google Scholar]
  25. Wu, C.; Kraft, P.; Zhai, K.; Chang, J.; Wang, Z.; Li, Y. Genome-wide association analyses of esophageal squamous cell carcinoma in Chinese identify multiple susceptibility loci and gene-environment interactions. Nat. Genet. 2012, 44, 1090. [Google Scholar] [CrossRef] [PubMed]
  26. Fratini, S.; Ragionieri, L.; Cannicci, S. Stock structure and demographic history of the Indo-West Pacific mud crab Scylla serrata. Estuarine Coast. Shelf Sci. 2010, 86, 51–61. [Google Scholar] [CrossRef]
  27. Chu, T.J.; Wang, D.; Haung, H.L.; Lin, F.J.; Tzeng, T.D. Genetic variations and expansion of whiskered velvet shrimp (Metapenaeopsis barbata) off China and Taiwan inferred from intron sequence. Biochem. Syst. Ecol. 2011, 39, 520–525. [Google Scholar] [CrossRef]
  28. Chu, T.J.; Wang, D.; Huang, H.L.; Lin, F.J.; Tzeng, T.D. Population structure and historical demography of the whiskered velvet shrimp (Metapenaeopsis barbata) off China and Taiwan inferred from the mitochondrial control region. Zool. Stud. 2012, 51, 99–107. [Google Scholar]
  29. Ma, H.Y.; Ma, C.Y.; Ma, L.B. Population genetic diversity of mud crab (Scylla paramamosain) in Hainan Island of China based on mitochondrial DNA. Biochem. Syst. Ecol. 2011, 39, 434–440. [Google Scholar] [CrossRef]
  30. Han, Z.Q.; Zhu, W.B.; Zheng, W.; Li, P.F.; Shui, B.O. Significant genetic differentiation between the Yellow Sea and East China Sea populations of cocktail shrimp Trachypenaeus curvirostris revealed by the mitochondrial DNA COI gene. Biochem. Syst. Ecol. 2015, 59, 78–84. [Google Scholar] [CrossRef]
  31. 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] [Green Version]
  32. Grosberg, R.K.; Cunningham, C.W. Genetic Structure in the Sea: From Populations to Communities. In Marine Community Ecology; Bertness, M.D., Gaines, S.D., Hay, M.E., Eds.; Sinauer Associates Publishers Inc.: Sunderland, MA, USA, 2001; pp. 61–84. [Google Scholar]
  33. Gallagher, S.J.; Kitamura, A.; Iryu, Y.; Itaki, T.; Koizumi, I.; Hoiles, P.W. The Pliocene to recent history of the Kuroshio and Tsushima Currents: A multi-proxy approach. Prog. Earth Planet. Sci. 2015, 2, 2–17. [Google Scholar] [CrossRef] [Green Version]
  34. Wang, L.; Sarnthein, M.; Erlenkeuser, H.; Grimalt, J.; Grootes, P.; Heilig, S.; Ivanova, E.; Kienast, M.; Pelejero, C.; Pflaumann, U. East Asian monsoon climate during the late Pleistocene: High-resolution sediment records from the South China Sea. Mar. Geol. 1999, 156, 245–284. [Google Scholar] [CrossRef] [Green Version]
  35. Hansen, J.E.; Sato, M. Paleoclimate implications for human-made climate change. In Climate Change: Inferences from Paleoclimate and Regional Aspects; Berger, A., Mesinger, F., Sijacki, D., Eds.; Springer Publishers Inc.: New York, NY, USA, 2012. [Google Scholar]
  36. Li, Y.; Chen, J.; Jiang, L.; Qiao, G. Islands conserve high species richness and areas of endemism of Hormaphidinae aphids. Curr. Zool. 2017, 63, 623–632. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. 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 fathead mullet Mugil cephalus. BMC Evol. Biol. 2011, 11, 83. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Dall, W.; Hill, B.J.; Rothlisberg, P.C.; Staples, D.J. The biology of the Penaeidae. In Advances in Marine Biology; Blaxter, J.H.S., Southward, A.J., Eds.; Academic Press: San Diego, CA, USA, 1990. [Google Scholar]
  39. Dodd, R.S.; Afzal-Rafii, Z. Evolutionary genetics of mangroves: Continental drift to recent climate change. Trees 2002, 16, 80–86. [Google Scholar] [CrossRef]
  40. Duke, N.C.; Lo, E.; Sun, M. Global distribution and genetic discontinuities of mangroves—Emerging patterns in the evolution of Rhizophora. Trees 2002, 16, 65–79. [Google Scholar] [CrossRef]
  41. Ren, G.; Miao, G.; Ma, C.; Lu, J.; Yang, X.; Ma, H. Genetic structure and historical demography of the blue swimming crab (Portunus pelagicus) from southeastern sea of China based on mitochondrial COI gene. Mitochondrial DNA Part A 2018, 29, 192–198. [Google Scholar] [CrossRef]
  42. McMillen-Jackson, A.L.; Bert, T.M. Disparate patterns of population genetic structure and population history in two sympatric penaeid shrimp species (Farfantepenaeus aztecus and Litopenaeus setiferus) in the eastern United States. Mol. Ecol. 2003, 12, 2895–2905. [Google Scholar] [CrossRef]
  43. McMillen-Jackson, A.; Bert, T. Mitochondrial DNA variation and population genetic structure of the blue crab Callinectes sapidus in the eastern United States. Mar. Biol. 2004, 145, 769–777. [Google Scholar] [CrossRef]
  44. Gallagher, J.; Finarelli, J.A.; Jonasson, J.P.; Carlsson, J. Mitochondrial D-loop DNA analyses of Norway lobster (Nephrops norvegicus) reveals genetic isolation between Atlantic and East Mediterranean populations. J. Mar. Biol. Assoc. UK 2019, 99, 933–940. [Google Scholar] [CrossRef]
  45. Voris, H.K. Maps of Pleistocene sea levels in Southeast Asia: Shorelines, river systems and time durations. J. Biogeogr. 2000, 27, 1153–1167. [Google Scholar] [CrossRef] [Green Version]
  46. Sambrook, J.; Russell, D.W. Molecular Cloning: A Laboratory Manual, 3rd ed.; Cold Spring Harbor Publishers Inc.: New York, NY, USA, 2001. [Google Scholar]
  47. Rozas, J.; Sanchez-Delbarrio, J.C.; Messeguer, X.; Rozas, R. DnaSP, DNA polymorphism analyses by the coalescent and other methods. Bioinformatics 2003, 19, 2496–2497. [Google Scholar] [CrossRef] [PubMed]
  48. Nei, M. Molecular Evolutionary Genetics, 1st ed.; Columbia Publishers Inc.: New York, NY, USA, 1987. [Google Scholar]
  49. Tamura, K.; Peterson, D.; Peterson, N.; Stecher, G.; Nei, M.; Kumar, S. MEGA6: Molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol. Biol. Evol. 2011, 28, 2731–2739. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Bandelt, H.J.; Forster, P.; Rohl, A. Median-joining networks for inferring intraspecific phylogenies. Mol. Biol. Evol. 1999, 16, 37–48. [Google Scholar] [CrossRef]
  51. Nei, M.; Tajima, F. DNA polymorphism detectable by restriction endonucleases. Genetics 1981, 97, 145–163. [Google Scholar] [CrossRef] [PubMed]
  52. Excoffier, L.; Lischer, H.E. 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] [PubMed]
  53. Tzeng, T.D.; Yeh, S.Y.; Hui, C.F. Population genetic structure of the kuruma prawn (Penaeus japonicus) in East Asia inferred from mitochondrial DNA sequences. ICES J. Mar. Sci. 2004, 61, 913–920. [Google Scholar] [CrossRef]
  54. Stanley, H.F.; Casey, S.; Carnahan, J.M.; Goodman, S.; Harwood, J.; Wayne, K. Worldwide patterns of mitochondrial DNA differentiation in the harbor seal (Phoca vitulina). Mol. Biol. Evol. 1996, 13, 368–382. [Google Scholar] [CrossRef] [Green Version]
  55. Tajima, F. Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics 1989, 123, 585–595. [Google Scholar] [CrossRef]
  56. Tzeng, T.D.; Chu, T.J.; Wang, D.; Haung, H.L.; Yeh, S.Y. Population Structure in the Sword Prawn (Parapenaeopsis Hardwickii) from the East China Sea and Taiwan Strait Inferred from Intron Sequences. J. Crust. Biol. 2008, 28, 234–239. [Google Scholar] [CrossRef]
  57. Sezmiş, E. The Population Genetic Structure of Portunus pelagicus in Australian Waters; Murdoch University Publishers Inc.: Filosopi, Australia, 2004. [Google Scholar]
  58. Chan, B.K.K.; Tsang, L.M.; Chu, K.H. Morphological and genetic differentiation of the acorn barnacle Tetraclita squamosa (Crustacea, Cirripedia) in East Asia and description of a new species of Tetraclita. Zool. Scr. 2007, 36, 79–91. [Google Scholar] [CrossRef]
  59. Wong, K.J.H.; Chan, B.K.K.; Shih, H.T. Taxonomy of the sand bubbler crabs Scopimera globosa De Haan, 1835, and S. tuberculata Stimpson, 1858 (Crustacea: Decapoda: Dotillidae) in East Asia with description of a new species from the Ryukyus, Japan. Zootaxa 2010, 2345, 43–59. [Google Scholar] [CrossRef] [Green Version]
  60. Tsang, L.M.; Wu, T.H.; Shih, H.T.; Williams, G.A.; Chu, K.H.; Chan, B.K.K. Genetic and morphological differentiation of the Indo-West Pacific intertidal barnacle Chthamalus malayensis. Integr. Comp. Biol. 2012, 52, 388–409. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  61. Crandall, K.A.; Templeton, A.R. Empirical tests of some predictions from coalescent theory with applications to intraspecific phylogeny reconstruction. Genetics 1993, 134, 959–969. [Google Scholar] [CrossRef] [PubMed]
  62. Chiang, T.Y.; Schaal, B.A. Phylogeography of North American populations of the moss species Hylocomium splendens based on the nucleotide sequence of internal transcribed spacer 2 of nuclear ribosomal DNA. Mol. Ecol. 1999, 8, 1037–1042. [Google Scholar] [CrossRef]
  63. Wang, J.P.; Hsu, K.C.; Chiang, T.Y. Mitochondrial DNA phylogeography of Acrossocheilus paradoxus (Cyprinidae) in Taiwan. Mol. Ecol. 2000, 9, 1483–1494. [Google Scholar] [CrossRef]
  64. Hu, J.; Kawamura, H.; Hong, H.; Qi, Y. A Review on the currents in the South China Sea: Seasonal circulation, South China Sea warm current and Kuroshio intrusion. J. Oceanogr. 2000, 56, 607–624. [Google Scholar] [CrossRef]
  65. Yap, E.S.; Sezmis, E.; Chapin, J.A.; Potter, I.C.; Spencer, P.B.S. Isolation and characterization of micro satellite loci in Portunus pelagicus (Crustacea: Portunidae). Mol. Ecol. Notes. 2002, 2, 30–32. [Google Scholar] [CrossRef]
  66. Smith, G.S.; Sumpton, W.D. Behaviour of the commercial sand crab Portunus pelagicus (L.) at trap entrances. Asian Fish Sci. 1989, 3, 101–113. [Google Scholar]
  67. Sumpton, W.D.; Potter, M.A.; Smith, G.S. Reproduction and growth of the Commercial Sand Crab, Portunus pelagicus (L.) in Moreton Bay, Queensland. Asian Fish Sci. 1994, 7, 103–113. [Google Scholar]
  68. Kangas, M.I. Synopsis of the Biology and Exploitation of the Blue Swimmer Crab, Portunus pelagicus Linnaeus, in Western Australia; Fisheries Research Report No. 121; Fisheries Western Australia: Perth, Australia, 2000; pp. 1–22. [Google Scholar]
  69. Abdul Halim, S.A.A.; Othman, A.S.; Akib, N.A.M.; Jamaludin, N.A.; Esa, Y.; Nor, S.A.M. Mitochondrial markers identify a genetic boundary of the Green Tiger Prawn (Penaeus semisulcatus) in the Indo-Pacific Ocean. Zool. Stud. 2021, 60, e8. [Google Scholar]
Figure 1. Sampling areas and haplotype frequencies of P. pelagicus in Southeast Asia. Table 1 presents the numbers of lineages A and B at each sampling site.
Figure 1. Sampling areas and haplotype frequencies of P. pelagicus in Southeast Asia. Table 1 presents the numbers of lineages A and B at each sampling site.
Jmse 10 01071 g001
Figure 2. Neighbour-joining (NJ) tree based on mtDNA D-loop sequences with bootstrap values (NJ/ML, respectively) shown adjacent to the two corresponding lineages for Portunus pelagicus. The numbers at the nodes indicate bootstrap values (expressed as percentages) with 1000 replicates.
Figure 2. Neighbour-joining (NJ) tree based on mtDNA D-loop sequences with bootstrap values (NJ/ML, respectively) shown adjacent to the two corresponding lineages for Portunus pelagicus. The numbers at the nodes indicate bootstrap values (expressed as percentages) with 1000 replicates.
Jmse 10 01071 g002
Figure 3. The haplotype network of Portunus pelagicus in all sampling sites.
Figure 3. The haplotype network of Portunus pelagicus in all sampling sites.
Jmse 10 01071 g003
Figure 4. UPGMA tree showing relationships between the 6 sampling sites.
Figure 4. UPGMA tree showing relationships between the 6 sampling sites.
Jmse 10 01071 g004
Figure 5. The observed pairwise differences and the expected mismatch distributions under the sudden expansion model of swimming crabs: all populations; lineage A; lineage B.
Figure 5. The observed pairwise differences and the expected mismatch distributions under the sudden expansion model of swimming crabs: all populations; lineage A; lineage B.
Jmse 10 01071 g005
Table 1. Codes of sampling sites, sample sizes (n), number of haplotypes (nh), gene diversity (h), nucleotide diversity (π), Tajima’s D, and Fu’s Fs statistics for 6 populations of Portunus pelagicus in Southeast Asia.
Table 1. Codes of sampling sites, sample sizes (n), number of haplotypes (nh), gene diversity (h), nucleotide diversity (π), Tajima’s D, and Fu’s Fs statistics for 6 populations of Portunus pelagicus in Southeast Asia.
PopulationnnhLineage ALineage Bh Diversityπ DiversityTajima’s DFu’s Fs
XM22162200.970 ± 0.0220.045 ± 0.00462−0.871971.271
PH30201290.961 ± 0.0210.023 ± 0.00431−1.53515−1.445
KK36230360.962 ± 0.0190.020 ± 0.00114−1.03393−2.597
HK28262800.995 ± 0.0110.046 ± 0.00252−1.15632−4.889
VN27152700.949 ± 0.0220.025 ± 0.001940.56381.662
SGP23172300.972 ± 0.0200.032 ± 0.001520.03936−0.197
Lineage A10175 0.993 ± 0.0030.042 ± 0.00186−1.59138 *−24.23 *
Lineage B6541 0.979 ± 0.0070.020 ± 0.00097−1.30135−11.305 *
Total166116 0.994 ± 0.0020.064 ± 0.00120−1.0188−31.544 *
* p < 0.05. XM: Xiamen, PH: Penghu, KK: Kezailiao Kaohsiung, HK: Hong Kong, VN: Vietnam, SGP: Singapore.
Table 2. Matrix of pairwise FST (below diagonal) and p-values (above diagonal) among 6 populations of Portunus pelagicus in Southeast Asia.
Table 2. Matrix of pairwise FST (below diagonal) and p-values (above diagonal) among 6 populations of Portunus pelagicus in Southeast Asia.
HKXMVNSGPPHKK
HK-0.0830.0000.0000.0000.000
XM0.037-0.0000.0000.0000.000
VN0.1220.139-0.0000.0000.000
SGP0.1760.1620.250-0.0000.000
PH0.5880.5710.7030.447-0.090
KK0.5320.5140.6290.3890.035-
XM: Xiamen, PH: Penghu, KK: Kezailiao Kaohsiung, HK: Hong Kong, VN: Vietnam, SGP: Singapore.
Table 3. AMOVA results for 6 populations of Portunus pelagicus in Southeast Asia.
Table 3. AMOVA results for 6 populations of Portunus pelagicus in Southeast Asia.
PopulationGroupingSource of Variation Percentage of VariationΦ Statisticsp (More Extreme Value)
One groupGroup 1 (HK, XM, VN, SGP, PH, KK)AP41.98ΦST = 0.41978***
WP58.02
Two groupsGroup 1 (HK, XM, VN, SGP)AG46.57ΦCT = 0.46573
Group 2 (PH, KK)AP/WG7.07ΦSC = 0.13235***
WP46.36ΦST = 0.53644***
Three groupsGroup 1 (HK, XM, VN)AG44.45ΦCT = 0.44445*
Group 2 (SGP)AP/WG4ΦSC = 0.07201***
Group 3 (PH, KK)WP51.55ΦST = 0.48446***
Four groupsGroup 1 (HK, XM)AG43.42ΦCT = 0.43419*
Group 2 (VN)AP/WG2.13ΦSC = 0.03759***
Group 3 (SGP)WP54.45ΦST = 0.45546***
Group 4 (PH, KK)
AG is the intergroup component of variance, AP/WG is the inter-population/intragroup component of variance, and WP is the intra-population component of variance. * Significant at * p < 0.05, *** p < 0.001, based on the permutation test.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lu, Y.-M.; Shih, C.-H.; Chen, P.-C.; Kao, W.-C.; Lee, Y.-C.; Han, Y.-S.; Tzeng, T.-D. Genetic Variations and Expansion of the Blue Swimmer Crab (Portunus pelagicus) in Southeast Asia. J. Mar. Sci. Eng. 2022, 10, 1071. https://doi.org/10.3390/jmse10081071

AMA Style

Lu Y-M, Shih C-H, Chen P-C, Kao W-C, Lee Y-C, Han Y-S, Tzeng T-D. Genetic Variations and Expansion of the Blue Swimmer Crab (Portunus pelagicus) in Southeast Asia. Journal of Marine Science and Engineering. 2022; 10(8):1071. https://doi.org/10.3390/jmse10081071

Chicago/Turabian Style

Lu, Yu-Ming, Chun-Han Shih, Po-Cheng Chen, Wei-Chieh Kao, Ying-Chou Lee, Yu-San Han, and Tzong-Der Tzeng. 2022. "Genetic Variations and Expansion of the Blue Swimmer Crab (Portunus pelagicus) in Southeast Asia" Journal of Marine Science and Engineering 10, no. 8: 1071. https://doi.org/10.3390/jmse10081071

APA Style

Lu, Y. -M., Shih, C. -H., Chen, P. -C., Kao, W. -C., Lee, Y. -C., Han, Y. -S., & Tzeng, T. -D. (2022). Genetic Variations and Expansion of the Blue Swimmer Crab (Portunus pelagicus) in Southeast Asia. Journal of Marine Science and Engineering, 10(8), 1071. https://doi.org/10.3390/jmse10081071

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