Genetic Variation in the Asian Shore Crab <italic>Hemigrapsus sanguineus</italic> in Korean Coastal Waters as Inferred from Mitochondrial DNA Sequences
Genetic Variation in the Asian Shore Crab Hemigrapsus sanguineus in Korean Coastal Waters as Inferred from Mitochondrial DNA Sequences
Fisheries and aquatic sciences. 2012. Mar, 15(1): 49-56
Copyright ©2012, The Korean Society of Fisheries and Aquatic Science
This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License ( which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.
  • Received : October 10, 2011
  • Accepted : February 02, 2012
  • Published : March 30, 2012
Export by style
Cited by
About the Authors
Sung-Eic Hong
Department of Marine Bioscience, Gangneung-Wonju National University, Gangneung 210-702, Korea
Jin-Koo Kim
Department of Marine Biology, Pukyong National University, Busan 608-737, Korea
Jeong-Nam Yu
Environmental Research Complex, Incheon 404-708, Korea
Keun-Yong Kim
Department of Marine Biotechnology, Soonchunhyang University, Asan 336-745, Korea
Chung Il Lee
Department of Marine Bioscience, Gangneung-Wonju National University, Gangneung 210-702, Korea
Kwan Eui Hong
Yangyang Salmon Station, Korea Fisheries Resources Agency, Yangyang 215-851, Korea
Kie Young Park
Department of Marine Bioscience, Gangneung-Wonju National University, Gangneung 210-702, Korea
Genetic variation in the Asian shore crab Hemigrapsus sanguineus was determined from partial mitochondrial DNA (mtDNA) sequences of the cytochrome b ( Cyt b) gene. Samples included 143 crabs from six localities along three coastlines in South Korea. A nucleotide sequence analysis revealed 38 variable sites in a 470-bp sequence, which defined 37 haplotypes. The haplotypes were not associated geographically and had a shallow genealogy. Pairwise F ST tests and a two-dimensional scaling analysis revealed no significant genetic differentiation among most of the populations. The low pairwise comparison values, but significant genetic differentiation of a northeastern population from all other populations, might have been influenced by a restriction in gene flow caused by hydrographic conditions such as ocean boundaries. The high haplotype diversity, low nucleotide diversity, and time since H. sanguineus expansion in Korean coastal waters indicate rapid population growth and a recent, sudden expansion in the Late Pleistocene.
The Asian shore crab Hemigrapsus sanguineus is distrib-uted widely in the western Pacific, including Hong Kong, Taiwan, Korean, Chinese, and Japanese coastal waters, and as far north as Sakhalin Island, Russia (Sakai, 1976; Fukui, 1988; Dai and Yang, 1991; Hwang et al., 1993). It is abundant in rocky intertidal habitats and is an ecologically important predator in coastal ecosystems, as are other shore crab species (Kikuchi et al., 1981; Takada and Kikuchi, 1991; McDermott, 1998a, 1998b; Lohrer et al., 2000). In Japan, the breeding season of H. sanguineus is March-October, with a main peak May-June (Fukui, 1988). It has a planktonic larval stage of more than 1 month before developing into the juvenile crab (Fukui, 1988). The larval dispersal pattern and preferred habi-tat might have caused geographically distinct regional popula-tions to become homogeneous.
Estimating genetic structure among populations using mo-lecular markers has become a common approach to determin-ing sustainable yields and genetic diversity (Dunham, 2004). The population genetic structures of some marine species are influenced by their larval dispersal pattern and behavior after spawning, which are determined by oceanographic features including sea currents, hydrological conditions, and physical barriers (Doyle et al., 1993; Hsieh et al., 2010). In general, most marine species have limited population substructures and high levels of gene flow because of the effect of sea cur-rents. Ocean structure and dynamics, including current bound-aries and hydrographic conditions, have caused reproductive and partial genetic isolation in geographically distinct regional populations (Wares et al., 2001; Bilton et al., 2002).
Although DNA markers are expected to overcome defi-ciencies in allozyme analysis by increasing the accuracy and resolution of population-structure assessments in crab species (McMillen-Jackson et al., 1994; Creasey et al., 1997), there are few reports on the use of mitochondrial DNA (mtDNA) to measure genetic variation in shore crabs (Cassone and Bould-ing, 2006). Maternally inherited mtDNA has greater sequence variability than do most nuclear genes (Brown et al., 1979). Moreover, it has a compact genome size in both conserved and
PPT Slide
Lager Image
Sampling locations of six Asian shore crab Hemigrapsus sanguineus populations and ocean current flow pattern in the Korean waters. BUA, Buan; GUJ, Geoje; GUR, Guryonpo; HAE, Haenam; JUM, Jumunjin; TAE, Taean.
Sampling sites, and dates, geographical coordinates, number of individuals examined (n), haplotype and nucleotide diversity (π) of Asian shore crab Hemigrapsus sanguineus populations
PPT Slide
Lager Image
Sampling sites, and dates, geographical coordinates, number of individuals examined (n), haplotype and nucleotide diversity (π) of Asian shore crab Hemigrapsus sanguineus populations
variable regions and a conserved gene content and arrange-ment (Anderson et al., 1981). Therefore, mtDNA analysis has become a key method in evolutionary and ecological studies of crabs (Pfeiler et al., 2005; Cassone and Boulding, 2006; Azuma et al., 2008; Wang et al., 2008).
This study investigated the genetic variation and popula-tion structure of H. sanguineus in Korean coastal waters using mtDNA cytochrome b ( Cyt b) sequences to assess phylogeo-graphic and demographic patterns.
Materials and Methods
- Sampling
Muscle samples were taken from 143 live crabs collected from six localities along three coastlines of Korea from 2009 to 2010 ( Table 1 , Fig. 1 ). The collected samples were stored at -20℃ or kept in 100% ethanol at room temperature until used.
- PCR amplification and sequence analysis
Genomic DNA was extracted from about 20 mg of each specimen by a PUREGENE DNA isolation kit (Gentra Sys-tems, Minneapolis, MN, USA) following the manufacturer’s instructions. The purified DNA was dried at room temperature and dissolved in TE buffer (10 mM Tris-HCl, 1 mM EDTA; pH 8.0). To amplify the Cyt b gene, a pair of degenerate primers (Cyt9237 5′-GTWGCHCAYATTTGYCGAGA-3′; Cyt10050 5′-ACWGGKCGWGCWCCAATTCA-3′) with expected amplicon sizes of 850 bp was designed based on mtDNA sequences of the closely related crab species available in GenBank (AB093006, AY659990, AY562127, FJ797435, FJ827758-827761). The PCR amplification was performed with a thermocycler DNA Engine (MJ Research, Tokyo, Ja-pan). in a 20-μL reaction volume containing 1-2 μL of genom-ic DNA, 2 μM of each primer, 0.25 mM of each dNTP, 1 unit of Takara LA Taq DNA polymerase (Takara Shuzo, Kyoto, Japan), and 2 μL of 10× LA Taq reaction buffer (Takara Shuzo). The PCR conditions consisted of preheating at 94℃ for 5 min, followed by 35 cycles of 94℃ for 30 s, 55℃ for 30 s, and 72℃ for 30 s, with final extensions at 72℃ for 5 min. The fragment size of the PCR product was verified using 2% aga-rose gel electrophoresis after ethidium bromide staining. The PCR product was purified with the AccuPrep PCR Purification Kit (Bioneer, Daejeon, Korea). After cycle sequencing with the ABI PRISM BigDye Terminator v3.1 Cycle Sequencing Ready Reaction Kit (Applied Biosystems, Foster City, CA, USA), the purified PCR product was sequenced directly on an ABI 3730xl DNA Analyzer (Applied Biosystems) with two newly designed internal primers, HsCytF (5′-GGGGT-CAAATATCATTCTGG-3′) and HsCytR (5′-GCCTTTG-GAATTTTGAAGAG-3′).
- Data analysis
The sequence fragments obtained in this study were aligned with GENETIX-WIN ver. 4.0.1 (Software Development, To-kyo, Japan) to identify sequence variants. The integrated soft-ware package DnaSP version 5 (Librado and Rozas, 2009) was used to determine the genotypes or haplotypes. A parsi-mony network connecting the observed haplotypes to resolve the genealogy was plotted with TCS version 1.21 (Clement et al., 2000). Genetic variation within the populations, expressed as haplotype diversity ( h ) and nucleotide diversity (π), was estimated according to Nei (1987) based on Kimura’s two-parameter distance method using K and DA in the REAP pro-gram (McElroy et al., 1992). Pairwise population F ST values were calculated to estimate the genetic differentiation among the populations according to Slatkin and Hudson (1991) us-ing ARLEQUIN version 3.1 (Excoffier et al., 2005). The sig-nificance of each F ST value was tested using 10,000 random permutations. Gene flow among populations was estimated with N e m, the number of migrations per generation between population pairs (Slatkin, 1993) using the following equation:
  • Nem = (1/FST-1)/4,
where N is the effective population size, m is the effective proportion of immigrants, and F ST is the fixation index. To ex-amine the relationships among populations visually, we used a two-dimensional scaling (TDS) analysis based on the matrix of pairwise F ST values from the Cyt b sequence data calculated by SPSS version 14.0K (SPSS, Inc., Chicago, IL, USA).
Neutral expectation and historic demographic expansions were investigated by examining Fu’s F S and Tajima’s D and mismatch distributions with the sudden expansion model (Rogers and Harpending, 1992). A goodness-of-fit test was used to test the validity of the sudden-expansion model using a parametric bootstrap approach based on the sum of squared deviations (SSD) to compare the observed and estimated mis-match distributions (Schneider and Excoffier, 1999). Both the neutrality test and mismatch distribution analysis were per-formed in ARLEQUIN version 2000 (Schneider et al., 2000). Since the mutation rate of the H. sanguineus Cyt b gene over the estimated time since expansion was unknown, a molecu-lar clock was calculated using the sequence-divergence rates of mitochondrial protein coding regions for other marine crustaceans; rates ranged from 2.2 to 3.1% per million years (Knowlton and Weigt, 1998; Schubart et al., 1998).
- Cytb variation in H. sanguineus
The degenerate primers newly designed in this study (Cyt9237 and Cyt10050) successfully amplified the mtDNA Cyt b region of the 143 H. sanguineus individuals. Direct se-quencing of the PCR products with two internal primers (Hs-CytF and HsCytR) yielded a fragment with an amplicon size of 470 bp, which revealed 38 variable nucleotide sites defin-ing 37 haplotypes among the populations ( Table 2 ). The vari-able nucleotide sites observed consisted of 33 transitions and four transversions. All substitutions were biallelic except one, which was triallelic, suggesting the occurrence of a single base substitution among sequences. The nucleotide sequences of the 37 haplotypes have been deposited in the DDBJ/GenBank database under accession numbers AB570203-AB570239.
The parsimony network of the Cyt b haplotypes in H. san-guineus did not provide evidence of geographical associa-tion ( Fig. 2 ). In the network, two focal haplotypes, HeS1 and HeS12, were abundant, whereas the others, including single-tons, were rare and radiated from these focal haplotypes. The distribution of the 37 haplotypes among the six H. sanguin-eus populations is presented in Table 3 . Among them, 28 were found at single localities, and nine (HeS1, HeS4, HeS9,
PPT Slide
Lager Image
A single minimum spanning tree for the 38 mitochondrial cytochrome b region haplotypes of Asian shore crab Hemigrapsus sanguineus (Table 2). Circle sizes reflect haplotype abundance..
Variable nucleotide sites in 470 bp sequences of mitochondrial cytochrome b gene partial sequences in Asian shore crab Hemigrapsus sanguineus
PPT Slide
Lager Image
Variable nucleotide sites in 470 bp sequences of mitochondrial cytochrome b gene partial sequences in Asian shore crab Hemigrapsus sanguineus
Distribution of mitochondrial cytochrome b haplotypes among 6 populations of Asian shore crab Hemigrapsus sanguineus. Sam-pling site abbreviations (first row) are listed inTable 1TAE, Taean; BUA, Buan; HAE, Haenam; GUJ, Geoje; GUR, Guryonpo; JUM, Jumunjin.
PPT Slide
Lager Image
Distribution of mitochondrial cytochrome b haplotypes among 6 populations of Asian shore crab Hemigrapsus sanguineus. Sam-pling site abbreviations (first row) are listed in Table 1 TAE, Taean; BUA, Buan; HAE, Haenam; GUJ, Geoje; GUR, Guryonpo; JUM, Jumunjin.
HeS12, HeS21, HeS26, HeS30, HeS32, and HeS35) were ob-served in two or more localities. A number of individuals from the populations examined had the HeS12 and HeS26 haplo-types, although these were not found in JUM. HeS1 occurred in all populations.
The haplotype and nucleotide diversities of H. sanguineus demonstrated varying levels of genetic variation among populations. Variation was higher in the Guryonpo (GUR), Haenam (HAE), and Jumunjin (JUM) populations ( Table 1 ). Reduced haplotype diversity was observed in the Taean (TAE), Buan (BUA), and Geoje (GUJ) populations.
- Population genetic analysis
The pairwise population F ST estimates and migration rate based on the Cyt b sequences are presented in Table 4 . Most pairwise values were not significantly different among popu-lations. Although the pairwise F ST estimates showed no clear genetic differentiation (low values among populations), the
PPT Slide
Lager Image
Two-dimensional scaling analysis based on the matrix of pairwise FST values from the mitochondrial cytochrome b sequence data for Asian shore crab Hemigrapsus sanguineus. BUA, Buan; GUJ, Geoje; GUR, Guryonpo; HAE, Haenam; JUM, Jumunjin; TAE, Taean.
PPT Slide
Lager Image
Mismatch distribution constructed using pairwise differences among mitochondrial DNA (mtDNA) haplotypes of Asian shore crab Hemigrapsus sanguineus from the pooled populations and from a single northeastern population (Jumunjin [JUM]). Solid lines, observed frequency; dashed lines, frequency distribution expected from a sudden expansion model.
pairwise population F ST estimates were large when JUM was compared with the other populations. These findings were evi-dent in the TDS analysis and N e m values ( Fig. 3 ). Overall, the TAE, BUA, HAE, GUJ, and GUR populations were geneti-cally close to one another, but distinct from JUM.
Mismatch distributions for all populations pooled and the northernmost population (JUM) are shown in Fig. 4 . The pooled populations and JUM had no additional distribution peaks; the highest frequency occurred at one difference. The D , Fu’s F S , and SSD ( Table 5 ) values indicated expansion in these populations. Tajima’s D and Fu’s F S were significantly negative, with markedly reduced SSD in the pooled popula-
Pairwise FST(below the diagonal) and Nem (above the diago-nal) estimates among the six populations of Asian shore crab Hemigrapsus sanguineus.TAE, Taean; BUA, Buan; HAE, Haenam; GUJ, Geoje; GUR, Guryonpo; JUM, Jumunjin.*Significant differentiation (P< 0.05) on the exact test (Raymond and Rousset, 1995).
PPT Slide
Lager Image
Pairwise FST (below the diagonal) and Nem (above the diago-nal) estimates among the six populations of Asian shore crab Hemigrapsus sanguineus. TAE, Taean; BUA, Buan; HAE, Haenam; GUJ, Geoje; GUR, Guryonpo; JUM, Jumunjin. *Significant differentiation (P < 0.05) on the exact test (Raymond and Rousset, 1995).
tion group and JUM population, a parameter combination that strongly supported sudden expansion. Sudden expansion of the pooled population group was estimated to have occurred 0.042-0.060 million years ago (Ma) and that in the JUM popu-lation, 0.049-0.069 Ma ( Table 5 ).
The genetic variation analysis of the mitochondrial Cyt b gene revealed the followings: (i) haplotype and nucleotide diversity occurred at various levels of genetic variation and (ii) there was a low level of genetic differentiation in coastal areas of Korea, but there was some genetic differentiation of the northernmost population (JUM) in the East Sea from the other populations. The observed haplotypes from the Cyt b region in H. sanguineus were arrayed in star-like genealogies, each with several closely related, low-frequency haplotypes around a central high-fre-quency haplotype, suggesting a shallow haplotype genealogy. The star-like pattern and shallow genealogy indicate recent appearance and rapid population growth (Slatkin and Hudson, 1991; Rogers and Harpending, 1992). Our estimates showed that H. sanguineus began expanding 42,862-60,396 years ago ( Table 5 ). High haplotype diversity (0.616-0.880) but low nucle-otide diversity (0.0017-0.0034) within the populations indicates that the populations might have experienced historically rapid population growth from an ancestral population with a small effective population size in the Late Pleistocene (Avise, 2000).
Asian shore crab Hemigrapsus sanguineus: parameters of the sudden expansion model and estimated time since expansionD, Tajima’s D;FS, Fu’sFS; SSD, sum of squared deviations in the goodness of-fit-test; τ , time since expansion measured in mutational time units; CI, confi-dence interval; θ0and θ1, population sizes scaled by mutation rate before and after expansion, respectively.
PPT Slide
Lager Image
Asian shore crab Hemigrapsus sanguineus: parameters of the sudden expansion model and estimated time since expansion D, Tajima’s D; FS, Fu’s FS; SSD, sum of squared deviations in the goodness of-fit-test; τ , time since expansion measured in mutational time units; CI, confi-dence interval; θ0 and θ1, population sizes scaled by mutation rate before and after expansion, respectively.
A reduced genetic diversity in the TAE and BUA populations from the west coast in the Yellow Sea compared with the other regional populations, except the GUJ population, was observed from the haplotype and nucleotide diversities. The western Ko-rean populations showed no genetic differentiation from the southern populations as inferred from the pairwise F ST values. Therefore, crabs might have been introduced continuously with high rates of gene flow to the TAE and BUA populations from other sources following a reduction in effective population size.
The estimated pairwise F ST values, N e m rate, and TDS analy-sis indicate that substantial gene flow has occurred among the populations, suggesting that larval behavior and sea currents are responsible for the high rates of gene flow. The passive dispersal of planktonic larvae and sedentary lifestyle of adult marine in-vertebrates limit the formation of population substructure (Les-sios et al., 2003; Waters and Roy, 2004). Many other marine crab species reportedly have high levels of gene flow between populations or regions (Beckwitt, 1985; Merkouris et al., 1998; Azuma et al., 2008). Hwang et al. (1993) reported that the plank-tonic larval stage of H. sanguineus persists for 25 days until metamorphosis on the sea floor to form the first crab stage. In general, crustacean species with similar persistent pelagic larval stages have a high dispersal potential that produces genetic ho-mogeneity among local populations (Palumbi, 1994). The Tsu-shima Warm Current (TC) branches off the Kuroshio Current, with part of the TC running into the Yellow Sea and the main part entering the East Sea along the Korean Peninsula ( Fig. 1 ). Therefore, the TC might transport H. sanguineus larvae to the western coast of Korea.
Despite the lack of geographically associated haplotypes and genetic structure within and among populations, our analyses indicated a degree of genetic differentiation between the north-ernmost population (JUM) in the East Sea and the other popula-tions. The sub-polar front in the East Sea is similar to the west-ern boundary current in that a polar front forms at the boundary between the low-temperature, low-salinity waters of the north-ern region and the high-temperature, high-salinity waters of the southern region (Rhein et al., 1995; Pickart et al., 1997). The sub-polar front, which extends along the coast of Japan before turning abruptly at the Noto Peninsula frontal region toward the center of the East Sea, has a close relationship to the TC and cold-water currents, including the Liman Current (LC) (Senjyu, 1999; Ichikawa and Beardsley, 2002) ( Fig. 1 ). Restricted or se-lective gene flow from main distributions can result from such physical barriers, which restrict larval transport (Hedgecock, 1986; Scheltema, 1986; Bowen and Avise, 1990; Palumbi, 1994; Burton, 1998). The JUM population is located at the sub-polar front, where the cold and warm currents of the East Sea meet. Therefore, the balanced effects of the cold (LC) and warm (TC) water currents in the East Sea might explain the genetic differen-tiation between the JUM population and other wild populations.
The mismatch distribution in our data was unimodal, and the neutrality test gave a significantly negative value, suggesting re-cent population expansion of H. sanguineus in Korea. Estimates of other mismatch distribution parameters corroborate this evi-dence. Therefore, H. sanguineus population expansion in Ko-rean coastal waters resulted from rapid population growth and recent, sudden expansion in the Late Pleistocene. This was cor-roborated by the star-like genealogy of the haplotypes, high hap-lotype diversities, close genetic similarities among haplotypes, mismatch distribution pattern supporting a sudden expansion model, and estimated expansion time. However, this perspective remains ambiguous, and further extended sampling from Hong Kong through Japan to Russia is necessary to clarify the histori-cal influences.
Based on our results, the H. sanguineus populations might be one large panmictic population in Korean coastal waters, with the exception of the northeastern population. This genetic information about the current condition of H. sanguineus will be useful in developing a conservation strategy and subsequent ecological monitoring. H. sanguineus has colonized the eastern coast of the USA, Atlantic France, and the Netherlands in recent years. Comparative studies of genetic variation and structure in donor and invader populations using the molecular procedures employed in this study are essential for developing a greater un-derstanding of the spatiotemporal invasion-pattern mechanism. Finally, we demonstrated that genetic variation analysis using mtDNA Cyt b sequences is a useful model for research on popu-lation-level studies in closely related species.
This research was a part of the project titled “East Asian Seas Time Series - I (EAST-I)” and “Long-term change of structure and function in marine ecosystems of Korea” fund-ed by the Ministry of Land, Transport and Maritime Affairs, Korea.
Anderson S , Bankier AT , Barrell BG , de Bruijn MHL , Coulson AR , Drouin J , Eperon IC , Nierlich DP , Roe BA , Sanger F , Schreier PH , Smith AJH , Staden R , Young IG 1981 Sequence and organi-zation of the human mitochondrial genome. Nature 290 457 - 465
Avise JC 2000 Phylogeography: the History and Formation of Species. Harvard University Press Cambridge, MA, US.
Azuma N , Kunihiro Y , Sasaki J , Mihara E , Mihara Y , Yasunaga T , Jin DH , Abe S 2008 Genetic variation and population structure of hair crab (Erimacrus isenbeckii) in Japan inferred from mitochon-drial DNA sequence analysis. Mar Biotechnol 10 39 - 48
Beckwitt R 1985 Population genetics of the sand crab, Emerita analoga Stimpson, in southern California. J Exp Mar Biol Ecol 91 45 - 52
Bilton DT , Paula J , Bishop JDD 2002 Dispersal genetic differen-tiation and speciation in estuarine organisms. Estuar Coast Shelf Sci 55 937 - 952
Bowen BW , Avise JC 1990 Genetic structure of Atlantic and Gulf of Mexico populations of sea bass, menhaden, and sturgeon influ-ence of zoogeographic factors and life-history patterns. Mar Biol 107 371 - 381
Brown WM , George M Jr , Wilson AC 1979 Rapid evolution of ani-mal mitochondrial DNA. Proc Natl Acad Sci U S A 76 1967 - 1971
Burton RS 1998 Intraspecific phylogeography across the Point Con-ception biogeographic boundary. Evolution 52 734 - 745
Cassone BJ , Boulding EG 2006 Genetic structure and phylogeog-raphy of the lined shore crab, Pachygrapsus crassipes, along the northeastern and western Pacific coasts. Mar Biol 149 213 - 226
Clement M , Posada D , Crandall KA 2000 TCS: a computer pro-gram to estimate gene genealogies. Mol Ecol 9 1657 - 1659
Creasey S , Rogers AD , Tyler P , Young C , Gage J 1997 The popula-tion biology and genetics of the deep-sea spider crab, Encephaloi-des armstrongi Wood-Mason 1891 (Decapoda: Majidae). Philos Trans R Soc Lond B Biol Sci 352 365 - 379
Dai AY , Yang SL 1991 Crabs of the China Seas. China Ocean Press Beijing, CN.
Doyle MJ , Morse WW , Kendall AW Jr 1993 A comparison of lar-val fish assemblages in the temperate zone of the northeast Pacific and northwest Atlantic Oceans. Bull Mar Sci 53 588 - 644
Dunham RA 2004 Aquaculture and Fisheries Biotechnology. CABI Publishing Press Wallingford, GB.
Excoffier L , Laval G , Schneider S 2005 Arlequin (version 3.0): an integrated software package for population genetics data analysis. Evol Bioinform Online 1 47 - 50
Fukui Y 1988 Comparative studies on the life history of the grapsid crabs (Crustacea Brachyura) inhabiting intertidal cobble and boul-der shores. Publ Seto Mar Biol Lab 33 121 - 162
Hedgecock D 1986 Is gene flow from pelagic larval dispersal important in the adaptation and evolution of marine invertebrates? Bull Mar Sci 39 550 - 564
Hsieh HY , Lo WT , Liu DC , Su CW 2010 Influence of hydrographic features on larval fish distribution during the south-westerly mon-soon in the waters of Taiwan, western North Pacific Ocean. J Fish Biol 76 2521 - 2539
Hwang SG , Lee C , Kim CH 1993 Complete larval development of Hemigrapsus sanguineus (Decapoda, Brachyura, Grapsidae) reared in the laboratory. Korean J Syst Zool 9 69 - 86
Ichikawa H , Beardsley RC 2002 The current system in the Yellow and East China Seas. J Oceanogr 58 77 - 92
Kikuchi T , Tanaka M , Nojima S , Takahashi T 1981 Ecological studies on the pebble crab, Gaetice depressus (de Haan). I. Eco-logical distribution of the crab and environmental conditions. Publ Amakusa Mar Biol Lab 6 23 - 34
Knowlton N , Weigt LA 1998 New dates and new rates for diver-gence across the Isthmus of Panama. Proc R Soc Lond B Biol Sci 265 2257 - 2263
Lessios HA , Kane J , Robertson DR 2003 Phylogeography of the pantropical sea urchin Tripneustes: Contrasting patterns of popula-tion structure between oceans. Evolution 57 2026 - 2036
Librado P , Rozas J 2009 DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics 25 1451 - 1452
Lohrer AM , Fukui Y , Wada K , Whitlatch RB 2000 Structural com-plexity and vertical zonation of intertidal crabs, with focus on habi-tat requirements of the invasive Asian shore crab, Hemigrapsus sanguineus (de Haan). J Exp Mar Biol Ecol 244 203 - 217
McDermott JJ 1998a The western Pacific brachyuran Hemigrapsus sanguineus (Grapsidae) in its new habitat along the Atlantic coast of the Untied States: reproduction. J Crustac Biol 18 308 - 316
McDermott JJ 1998b The western Pacific brachyuran (Hemigrapsus sanguineus: Grapsidae), in its new habitat along the Atlantic coast of the Untied States: geographic distribution and ecology. ICES J Mar Sci 55 289 - 298
McElroy D , Moran P , Bermingham E , Kornfield I 1992 REAP: An integrated environment for the manipulation and phylogenetic analysis of restriction data. J Hered 83 157 - 158
McMillen-Jackson AL , Bert TM , Steele P 1994 Population genetics of the blue crab Callinectes sapidus: modest population structuring in a background of high gene flow. Mar Biol 118 53 - 65
Merkouris SE , Seeb LW , Murphy MC 1998 Low levels of genetic diversity in highly exploited populations of Alaskan tanner crabs, Chionoecetes bairdi, and Alaskan and Atlantic snow crabs, C. opilio. Fish Bull 96 525 - 537
Nei M 1987 Molecular Evolutionary Genetics. Columbia University Press New York, US.
Palumbi SR 1994 Genetic divergence, reproductive isolation, and ma-rine speciation. Annu Rev Ecol Syst 25 547 - 572
Pfeiler E , Hurtado LA , Knowles LL , Torre-Cosío J , Bourillón-Moreno L , Márquez-Farías JF , Montemayor-López G 2005 Population genetics of the swimming crab Callinectes bellicosus (Brachyura: Portunidae) from the eastern Pacific Ocean. Mar Biol 146 559 - 569
Pickart RS , Spall MA , Lazier JRN 1997 Mid-depth ventilation in the western boundary current system of the sub-polar gyre. Deep Sea Res Part I Oceanogr Res Pap 44 1025 - 1054
Raymond M , Rousset F 1995 An exact test for population differen-tiation. Evolution 49 1280 - 1283
Rhein M , Stramma L , Send U 1995 The Atlantic deep western boundary current: water masses and transports near the equator. J Geophys Res 100 2441 - 2457
Rogers AR , Harpending H 1992 Population growth makes waves in the distribution of pairwise genetic differences. Mol Biol Evol 9 552 - 569
Sakai T 1976 Crabs of Japan and the Adjacent Seas. Kodansha Ltd Tokyo, JP.
Scheltema RS 1986 On dispersal and planktonic larvae of benthic in-vertebrates: an eclectic overview and summary of problems. Bull Mar Sci 39 290 - 322
Schneider S , Excoffier L 1999 Estimation of past demographic parameters from the distribution of pairwise differences when the mutation rates vary among sites: application to human mitochon-drial DNA. Genetics 152 1079 - 1089
Schneider S , Roessli D , Excoffier L 2000 Arlequin: a Software for Population Genetics Data Analysis. User Manual Ver. 2.000. Genetics and Biometry Laboratory, Department of Anthropology University of Geneva Geneva, CH.
Schubart CD , Diesel R , Hedges SB 1998 Rapid evolution to ter-restrial life in Jamaican crabs. Nature 393 363 - 365
Senjyu T 1999 The Japan Sea intermediate water: its characteristics and circulation. J Oceanogr 55 111 - 122
Slatkin M 1993 Isolation by distance in equilibrium and non-equilibri-um populations. Evolution 47 264 - 279
Slatkin M , Hudson RR 1991 Pairwise comparisons of mitochon-drial DNA sequences in stable and exponentially growing popula-tions. Genetics 129 555 - 562
Takada Y , Kikuchi T 1991 Seasonal and vertical variation of the boulder shore fauna in Amakusa. Publ Amakusa Mar Biol Lab 11 1 - 17
Wang C , Li C , Li S 2008 Mitochondrial DNA-inferred population structure and demographic history of the mitten crab (Eriocheir sensu stricto) found along the coast of mainland China. Mol Ecol 17 3515 - 3527
Wares JP , Gaines SD , Cunningham CW 2001 A comparative study of asymmetric migration events across a marine biogeographic boundary. Evolution 55 295 - 306
Waters JM , Roy MS 2004 Phylogeography of a high-dispersal New Zealand sea-star: does upwelling block gene-flow? Mol Ecol 13 2797 - 2806