Causations of phylogeographic barrier of some rocky shore species along the Chinese coastline
© Wang et al. 2015
Received: 25 September 2014
Accepted: 22 May 2015
Published: 15 June 2015
Substrate, ocean current and freshwater discharge are recognized as important factors that control the larval dispersal and recruitment of intertidal species. Life history traits of individual species will determine the differential responses to these physical factors, and hence resulting in contrasting phylogeography across the same biogeographic barrier. To determine how these factors affect genetic structure of rocky shore species along the China coast, a comparative phylogeographic study of four intertidal and subtidal species was conducted using mitochondrial and nuclear DNA by combining new sequences from Siphonaria japonica with previously published sequences from three species (Cellana toreuma, Sargassum horneri and Atrina pectinata).
Analysis of molecular variance and pairwise ΦST revealed significant genetic differences between the Yellow Sea (YS) and the other two marginal seas (East China Sea, ECS and South China Sea, SCS) for rocky-shore species (S. japonica, C. toreuma, S. horneri), but not for muddy-shore species Atrina pectinata. Demographic history analysis proved that the population size of all these four species were persistent though the Last Glacial Maximum (LGM, ~20 ka BP). Migration analysis revealed that gene flow differentiated northward and southward migration for these four species. However, the inferred direction of gene flow using alternatively mitochondrial or nuclear markers was contradictory in S. japonica.
It is concluded that there is a phylogeographical break at the Yangtze River estuary for the rocky shore species and the causation of the barrier is mainly due to the unsuitable substratum and freshwater discharge. All four intertidal and subtidal species appear to have persisted through the LGM in China, indicating the lower impact of LGM on intertidal and subtidal species than generally anticipated. The imbalanced gene flow between YS and ESCS groups for these four species could be explained by historical refugia. The discordance between mitochondrial and nuclear markers in the MIGRATE analysis of S. japonica prove the importance of employing multi-locus data in biogeographic study. Climate change, land reclamation and dam construction, which are changing substrate and hydrological conditions around Yangtze River estuary, will consequently affect the biogeographic pattern of intertidal species.
Phylogeographical patterns of marine faunas are complex and affected by multiple biotic and abiotic factors, and it is in the long-term interest of marine ecologists to understand the roles of these factors in determining the distribution and genetic structuring of species. Glacial-interglacial climate fluctuations during the Pleistocene led to changes of sea level [1, 2] and then caused habitat contractions or expansions , which impeded gene flow of marine species and resulted in genetic divergence. Following the glacial retreat, demographic expansion occurred in most marine taxa . Postglacial exchanges of propagules may erase the signals of historic isolation. Accordingly, the biology of the species (e.g. dispersal capacity), availability of suitable habitat and ocean current regimes, which determine the contemporary level of gene flow, would significantly contribute to promote, maintain or homogenize the genetic divergence created by glacial periods. Suitable habitat is assumed to have a significant impact on the phylogeographical distribution of coastal species [5–10] .
The impacts of ocean currents on larval dispersal are variable. Sometimes oceanic currents can promote larval dispersal and enrich population connectivity [11–14]. However, converging ocean currents may also pose a potential barrier to gene flow to some extents . Additionally, a river outflow carrying large amount of sediments and freshwater discharge can influence physical and chemical characteristics (e.g. geomorphology, turbidity, salinity, nutrients and dissolved oxygen etc.) of estuarine and coastal water, which control the biotic dynamics in estuary and coastal ecosystems . This has significant impacts on the dispersal of marine taxa [5, 15, 17–20].
The marginal seas in the northwestern Pacific have changed dramatically in area and configuration, particularly during the Pleistocene glacial-interglacial cycles . Three of the marginal seas, namely the South China Sea (SCS), the Yellow Sea (YS) plus the East China Sea (ECS), and the Sea of Japan (SOJ), served as three independent glacial refugia and resulted in vicariant divergence in marine fauna [22–24]. On the other hand, the lack of genetic structuring in some sympatrically distributed taxa was attributed to their postglacial colonization of these regions . Yet, these hypotheses were not rigorously tested using comprehensive phylogenetic data because most of the studies performed in the region were based on a single marker (in most case the mitochondrial DNA) and largely biased to commercial species, of which the effect of anthropogenic introduction for aquaculture purposes was unknown. Hence, additional investigation using both mitochondrial and nuclear markers was advocated .
Summary of life history features of four intertidal and subtidal species. The life history features include period of spawning, larval type and larval duration for four intertidal species
Intertidal rocky shore
Egg ribbons + Pelagic larvae
Intertidal rocky shore
Intertidal and subtidal rocky shore
Muddy and sandy shore
Sampling sites and summary of diversity indices of Siphonaria japonica. For each locality, individual numbers (n), haplotype diversity (h) and nucleotide diversity (π) for both mitochondrial and nuclear DNA are listed
h (Mean ± S.D.)
π (Mean ± S.D.)
h (Mean ± S.D.)
π (Mean ± S.D.)
0.963 ± 0.029
0.0085 ± 0.0047
0.986 ± 0.016
0.0067 ± 0.0036
0.989 ± 0.019
0.0092 ± 0.0051
1.000 ± 0.016
0.0094 ± 0.0050
0.984 ± 0.021
0.0077 ± 0.0044
0.996 ± 0.014
0.0088 ± 0.0046
0.943 ± 0.037
0.0098 ± 0.0054
0.988 ± 0.016
0.0075 ± 0.0040
East China Sea
0.944 ± 0.035
0.0079 ± 0.0044
1.000 ± 0.009
0.0073 ± 0.0039
0.902 ± 0.052
0.0050 ± 0.0030
1.000 ± 0.013
0.0095 ± 0.0050
0.944 ± 0.037
0.0080 ± 0.0044
1.000 ± 0.010
0.0065 ± 0.0035
South China Sea
0.929 ± 0.036
0.0075 ± 0.0042
0.998 ± 0.009
0.0072 ± 0.0038
0.942 ± 0.037
0.0082 ± 0.0045
1.000 ± 0.011
0.0069 ± 0.0037
Cellana toreuma, Sargassum horneri and Atrina pectinata
Dong et al.  had suggested that the two populations of C. toreuma from YS have significant higher haplotype and nucleotide diversity as compared to the populations from ECS and SCS. Re-analysis of published data of S. horneri revealed that the difference between populations from YS and populations from ECS/SCS was significant in nucleotide diversity (P = 0.034), but not in haplotype diversity (P = 0.187). For A. pectinata, however, there were no significant difference between populations from YS and populations from ECS in both haplotype diversity and nucleotide diversity.
Pairwise genetic distances (ΦST) among locations of Siphonaria japonica. Pairwise ΦST values for nuclear sequence ITS and mitochondrial sequence COI of Siphonaria japonica are given in the upper and lower diagonals, respectively. *P < 0.05; **P < 0.01, ***P < 0.001
Pairwise genetic distances (ΦST) among groups for each species. The locations of each intertidal and subtidal species were divided into Yellow Sea (YS), East China Sea (ECS) and South China Sea (SCS) along the China Coast. Reference: Cellana toreuma, Dong et al. ; Sargassum horneri, Hu et al. ; Atrina pectinata, Liu et al. .*P < 0.05; ***P < 0.001
Rocky shore species
Muddy shore species
Results from analysis of molecular variance (AMOVA) in four intertidal and subtidal species. For each species, locations were divided into two groups (Yellow Sea group and East plus South China Seas group) according to the Yangtze River Estuary. Reference: Cellana toreuma, Dong et al. ; Sargassum horneri, Hu et al. ; Atrina pectinata, Liu et al. . **P < 0.01, ***P < 0.001
Among locations within groups
Cellana toreuma, Sargassum horneri and Atrina pectinata
One-group AMOVA indicated that the genetic variation among all samples was negative and insignificant in A. pectinata (ΦST = −0.01), but positive and significant in C. toreuma and S. horneri (Additional file 1: Table S1). At location level, re-analysis of pairwise ΦST showed that there were no significant differences between all eight locations in A. pectinata (Additional file 2: Table S2). However, there existed significant difference between locations from YS and locations from the ECS and SCS (except for FJ location) (Additional file 3: Table S3) in S. horneri. At group level, the YS group was significantly different from the ECS group and the SCS group and the difference between the ECS group and the SCS group was low and non-significant in S. horneri (Table 4). However, there was no significant difference between the YS group and the ECS group for A. pectinata (ΦST = −0.0001, P = 0.436, Table 4).
Under the grouping criteria mentioned above, AMOVA analysis for COI of C. toreuma showed that there were significant genetic variations among groups (19.39 %, P < 0.001) and within locations (80.82 %, P < 0.001) (Table 5). AMOVA analysis for COIII of S. horneri showed significantly high variations among locations within groups (15.14 %, P < 0.001) and variation within locations (85.49 %, P < 0.001). The variance among groups was relatively low and insignificant (0.37 %, P = 0.331; Table 5). For COI sequences of A. pectinata, variance components at all three levels were non-significant (Table 5).
Cellana toreuma, Sargassum horneri and Atrina pectinata
The haplotype networks of other three species (C. toreuma, S. horneri and A. pectinata) showed a pattern which did not exhibit obvious subdivision according to geographical locations (Fig. 2c, d and e).
Results of population expansion tests on the four species. Fu’s F S , Tajima’s D and mismatch distribution estimate (τ) and the real expansion time (t) with 90 % credibility intervals in parentheses for each species are included. All COI sequences of S. japonica were divided into two groups: northern group and southern group, according to the phylogenetic analyses (See the results in the phylogenetic analysis). Reference: Cellana toreuma, Dong et al. ; Sargassum horneri, Hu et al. ; Atrina pectinata, Liu et al. . *P < 0.05; **P < 0.01, ***P < 0.001
Fu’s F S
A posterior simulation-based analogue of Akaike’s information reiteration through MCMC (AICM) test for northern group revealed that the expansion model favored over the other two models (Additional file 6: Table S4), and the BSP showed a demographic expansion. The constant size model was the best-fitting demographic model for southern group (Additional file 6: Table S4) which exhibited a similar pattern to the BSP.
Cellana toreuma, Sargassum horneri and Atrina pectinata
Although mitochondrial sequences in both C. toreuma and S. horneri showed a significant value of ΦST, we reconstructed their demographic history with mismatch distribution and BSP because, no spatial subdivision of haplotypes of mitochondrial sequence existed in the network (Fig. 2c, d). The BSP analysis showed that the population size of C. toreuma had an extremely gentle increase since 260 ka BP (Additional file 7: Figure S3A), and the population size of S. horneri experienced a very slow increase since 194 ka BP with the mutation rate of 1.675 % MY−1 (Additional file 7: Figure S3B). The population size of A. pectinata increased significantly since 232 ka BP with the mutation rate of 0.775 % MY−1 (Additional file 7: Figure S3C). Mismatch distribution analysis showed that the shapes of mismatch distribution were unimodal (Fig. 3c, d, e) and the time of demographic expansion for C. toreuma, S. horneri and A. pectinata were 502.5 ka (90 % CI: 70.7-538.2 ka) BP, 119.3 ka (90 % CI: 47.7-195.5 ka) BP and 170.1 ka (90 % CI: 133.7-421.7 ka) BP respectively (Table 6). Therefore, the demographic expansion of the two groups of S. japonica and other three species significantly predated the Last Glacial Maxima (~20 ka BP).
AICM tests for C. toreuma and S. horneri showed that the constant size model was a better fit to the data than the BSP and expansion models (Additional file 6: Table S4). It was not surprising for these two rocky-shore specie as they both showed relatively flat BSPs. The expansion model was strongly favored for A. pectinata that showed a significant demographic expansion in BSP.
Because the present study mainly focused on the gene flow across the Yangtze River Estuary, locations of each species were divided into YS group and ESCS group to detect the possible phylogeographical break. In the present studies of S. japonica, contrasting patters of gene flow between mitochondrial and nuclear markers were suggested (Fig. 2a and b, Additional file 8: Table S5). For the mitochondrial sequences, there is southward gene flow from YS to ESCS (Nm (YS-ESCS) = 5.79-24.93), but little gene flow northwards from ESCS to YS (Nm (ESCS-YS) = 0.00-8.14). The ITS analyses, however, indicated there was almost no gene flow in a southward direction from YS to ESCS (Nm (YS-ESCS) = 0.00-4.55), compared with a high gene flow from ESCS to YS (Nm (ESCS-YS) = 18.37-138.29).
Cellana toreuma, Sargassum horneri and Atrina pectinata
Re-analysis of published data showed that significant higher southward gene flow could be observed in C. toreuma (Nm (YS-ESCS) = 43.41-359.13, Fig. 2c) and A. pectinata (Nm (YS-ESCS) = 34.87-3620.72, Fig. 2e). On the contrary, there is almost no gene flow in a southward direction (Nm (YS-ESCS) = 0.00-31.57, Fig. 2d) in S. horneri, while large gene flow can be observed from ESCS to YS (Nm (ESCS-YS) = 59.93-3629.67).
Phylogeography in intertidal and subtidal species along the China coast: northern and southern differentiation
Significant population differentiation existed among northern (Yellow Sea) and southern groups (East China Sea and South China Sea) of Siphonaria japonica, Cellana toreuma and Sargassum horneri respectively, indicating that there is a genetic break between northern and southern populations of rocky-shore species. For mitochondrial and nuclear sequences of S. japonica, pairwise ΦST analyses among different locations indicated that Yellow Sea locations were significantly different from locations from the East China Sea and the South China Sea (Table 3). When these locations were divided into three groups (YS, ECS and SCS) based on their geographic locations, the ΦST values between the YS group and the other two groups (ECS and SCS) were significant for limpet S. japonica and macroalga S. horneri. This result is similar to results for the limpet C. toreuma, and suggests that the Yellow Sea group is relatively isolated from the other two groups. AMOVA analyses for COI and ITS of S. japonica and for COI of C. toreuma showed that genetic differentiation among YS and ESCS groups accounted for a high proportion of the molecular variance among groups (Table 5), which suggested that populations were mainly grouped according to the Yangtze River Estuary. Phylogenetic analysis showed that two putative groups existed in COI of S. japonica, which also demonstrated a genetic break around the Yangtze River Estuary. A geographic genetic break can also be detected by changes in gene diversity between populations. Significant difference in genetic diversity between locations from YS and locations from ECS/SCS for these three rocky-shore species supported the suggestion that a genetic break existed between YS and ECS/SCS groups.
In contrast to the rocky-shore species, no obvious population structure for the muddy-shore species (Atrina pectinata) was found based on the results from AMOVA analysis and pairwise ΦST values. The absence of genetic difference between YS and ECS groups as detected for the muddy-shore species was similar to previous studies of Rapana venosa , Tegillarca granosa  and Cyclina sinensis .
During the Pleistocene glacial-interglacial cycles, areas and configurations of the marginal seas have changed dramatically in the West Pacific [21, 47]. When the sea level fell about 120-140 m over the past ~ 800 kyr , the Yellow Sea (YS) and the East China Sea (ECS) were reduced to an elongated trough, the Okinawa Trough, and the South China Sea (SCS) became a semi-enclosed gulf . ECS and SCS were separated by a land bridge which connected Taiwan and the continent . Previous comparative phylogeographical studies in marginal seas of the northwestern Pacific have suggested that historical isolation between SCS and ECS plays important roles for the present-day distribution of genetic variation of coastal species such as some fishes, crustaceans and muddy/sandy shore mollusks .
Multi-factors controlling the population genetic differentiation of rocky-shore species
The contrasting phylogeographical patterns between rocky and muddy intertidal species indicates that substrate plays important roles in the phylogeographical patterns of intertidal species. From Lianyungang, Jiangsu Province to Shaoxing, Zhejiang Province, the absence of appropriate habitats (~600 km salt marsh shore) could hamper the settlement of rocky intertidal species and consequently genetic exchange between Yellow Sea and East China Sea populations. The negative impacts of a lack of appropriate substrate on genetic connectivity have been reported in other intertidal fauna. For instance, the existence of long stretches of sandy beach serves as a barrier to dispersal between Cape St Lucia and Zinkwazi Beach in the limpet S. nigerrima in southeast Africa . Recently, a genetic analysis of S. japonica specimens collected from Yangguang Island (32°36′N, 121°08′E), an artificial island between Zhoushan (30°01′N, 122°06′E) and Lianyungang (34°36′N, 119°13′E), was carried out and revealed that this location is genetically similar to the location in Zhoushan. These results indicate that colonization of some rocky intertidal species can happen across the barrier if suitable habitat is provided (Huang XW, Wang W and Dong YW, unpublished observations). In contrast, as seen above, the species living in muddy substrate appear unaffected by the salt marsh around at the Yangtze River estuary.
Freshwater discharge influences hydrological condition nearby estuary and may have some impacts on gene flow between the Yellow Sea and East China Sea. During spring and early summer, the spawning season of S. japonica, the Taiwan Warm Current (TWC) and the China Coast Current (CCC) in ECS and SCS flowing northward  transport pelagic larvae from SCS to ECS. However, the freshwater discharge from the Yangtze River will cause deflection of the East China Sea coastal current. The size and distance of deflection is more prominent in spring and summer due to the increased amount of surface runoff in the rainy season . In spring and summer, the plume of water from the Yangtze River discharge shifts in a northerly direction in parallel with the TWC with a clockwise deflection (Fig. 1b). When it reaches its maximum flow, the Yangtze River discharge affects surrounding hydrological conditions [26, 28] and causes a decline in the salinity of the upper layer of the Kuroshio Current , which could influence the northward transport of larvae into the Yellow Sea. Dong et al.  also suggested that unique haplotype and higher genetic diversity in YS group were mainly contributed to the ocean current and freshwater discharge during the spawning season of C. toreuma. The impact of freshwater discharge on the phylogenetic distribution of marine species has been widely observed in previous studies. For example, the outflow from the Amazon River has been invoked as a major factor for the biogeographic break between Brazilian and Caribbean faunas (e.g. barnacle Chthamalus proteus  and surgeonfishes Acanthuridae sp ).
Although the same phylogeographic break around the Yangtze River estuary was observed for two rocky intertidal limpets (C. toreuma ; S. japonica, the present study), haplotype networks of these two limpets were different. The haplotype network of C. toreuma is a single star-like network  and S. japonica presents a relatively complex network with two putative groups (Fig. 2). C. toreuma lays its eggs into sea water directly and the length of the pelagic larval stage varies between 4–18 days in congenerics . S. japonica deposits gelatinous egg ribbons containing numerous small eggs on the rocky shore, and then the eggs develop into planktotrophic veliger larvae within ~15 days. The veliger larvae of S. japonica could at least maintain in the plankton about seven days (Wang W and Dong YW, unpublished observations). Therefore, the discrepancy of phylogeographic patterns between S. japonica and C. toreuma could be associated with their different reproductive modes and larval dispersal capabilities. The high dispersal capability of the brown macroalga S. horneri could be partly attributed to low variance among groups based on the AMOVA analysis. S. horneri can breed by both sexual reproduction and asexual reproduction, and can form floating mats which can drift about 1–5 months after being detached from the substratum [55–57].
Historical demography of intertidal and subtidal species along the China coast
Even though the expansion model and the constant size model received strong support over the BSP for the putative northern and southern groups in S. japonica, respectively, evidences from the mismatch distribution analysis and the BSP suggested northern group and southern group had different timings of population expansion events. The expansion time of northern group (313.9 ~ 682.4 ka BP) was earlier than that of southern group (104.4 ~ 259.3 ka BP). This corroborates recent evidence of the existence of a northern refugium in the Northwestern Pacific observed in other marine organisms (e.g. seaweed Ishige okamurae , limpet C. toreuma  and barnacle Chthamalus challengeri ). In northwestern Pacific, South China Sea (SCS) and East China Sea (ECS) are widely accepted as southern glacial refugia [22–24]. Thus, new evidences of the existence of a northern refugium in this area may be helpful to illustrate the influence both of past populations in glacial refugia and of contemporary gene flow in shaping current phylogeographical patterns in the future studies. Furthermore, the population sizes of these four species were persistent through the Last Glacial Maximum (LGM; ~20 ka BP). Recent meta-analyses of demographic history of intertidal rocky organisms in the northeastern Pacific  and northwestern Pacific  also converged on similar finding that majority of the species were not extirpated entirely by the LGM and regional persistence maybe more prevalent.
At the contemporary level, MIGRATE analyses revealed that migration rates were unbalanced between ESCS and YS groups of S. japonica. However, the inferred direction of gene flow using alternatively mitochondrial or nuclear markers was contradictory. COI analyses suggested that southward migration rate was almost ten times the northward migration rate. ITS sequences suggested the northward migration rate was significantly higher (Fig. 2b). Discordances between nuclear and mitochondrial data in animal biogeographic studies are not uncommonly reported . Adaptive introgression, demographic disparity and sex biased migration are commonly invoked as potential explanation of discordance observed . Sex biased gene flow appears unlikely for this hermaphroditic limpet with planktonic larvae transported by ocean current while adaptive introgression is difficult to test based on the present data. The mitochondrial genes have a higher mutation rates and smaller effective population size than the nuclear genes, and hence considered to be more informative in shallow relationships. Besides, the high imbalance between immigration and emigration rates in other three species (Fig. 2c, d and e) could be associated with historical refugia and suggested the likely existence of sources and sinks at metapopulation level. Recent studies emphasize the need of multi-locus data for accurate estimation of various population parameters [3, 61]. As a result, more data are required to determine the underlying mechanism leading to the discordance, and our results emphasize the importance of employing multi-locus data in biogeographic study.
Biogeography of intertidal species in China
The phylogeographic barrier of rocky intertidal species around Yangtze River estuary will possibly disappear under the coupled impacts of climate change and human activities. Firstly, northward shift of rocky intertidal species is ubiquitous in the scenario of climate change [62–65]. The coastal sea surface temperature (SST) and extreme hot days have continued to increase from 1982 to 2010 [66, 67]. The increasing temperature will potentially force the northward shift of intertidal species along China coast; Secondly, land reclamation has resulted in about 55 % loss of coastal wetland in China from 1949 to 2002 , and large numbers of artificial structures provide hard substrates in areas where these are generally absent and act as stepping stones to connect the Yellow Sea populations and East China Sea/South China Sea populations; Finally, about 50,000 dams constructed in the Yangtze River catchment from the late 1950 to 2003 have a storage capacity of 22 % of the annual water discharge (200 × 109 m3 in 2003) and result in a strong decrease of sediment and freshwater discharge in spring and summer . The changing hydrological condition around the Yangtze River estuary could enhance the possibility of larval dispersal across it. Overall, the potential northward shift of rocky intertidal species will change the biogeographic pattern along the Chinese coast and adaptive management should be considered for future management of rocky intertidal ecosystem in China.
A significant phylogeographic break, occurring around the Yangtze River Estuary, was observed for populations of rocky-shore species along the China coast. Substrate, ocean current and freshwater discharge are suggested as major factors determining the contemporary structure of rocky-shore species by limiting north–south dispersal of planktonic larvae. In addition, historical events and life history characteristics can also influence the contemporary phylogeographic distribution of intertidal and subtidal species. However, human activities are changing the habitat of intertidal species. Large-scale land reclamation activities can change the sedimentary substrate, and so provide suitable hard substrates for colonization by rocky-shore species. On the other hand, numerous dams constructed in the Yangtze River catchment decrease the riverine sediment supply to the sea and will impact on the environment of the Yangtze Delta and the nearby coastal ocean.
Sampling and sequencing
A DNA sequence database consisted of new sequences (cytochrome c oxidase subunit I gene, COI; internal transcribed spacer, ITS) from Siphonaria japonica (specimens were collected from nine rocky shore localities along the Chinese coastline between May 2012 and January 2013) and previously published data (mitochondrial gene) from three other intertidal and subtidal species in the Pacific Northwest (Fig. 2).
Partial sequences of COI and ITS in S. japonica were amplified by polymerase chain reaction (PCR). COI sequences were amplified with universal primers LCO1490 and HCO2198  and ITS sequences were amplified using primers its-1d and its-4r . PCRs were conducted in a 25-μL reaction volume containing 2.5 μL of 10 × buffer (Mg2+ Plus), 2 μL of 2.5 mM dNTPs, 1 μL of each 10 mM primers, 0.25 μL (1.25 U) of Taq DNA polymerase and 200 ng DNA template. Amplification was initiated with denaturing at 95 °C for 3 min, followed by 35 cycles of 95 °C for 1 min, annealing at 40 °C for COI and 54 °C for ITS for 1 min and 72 °C for 1 min and then a final extension at 72 °C for 10 min. After visualizing the target amplicon in 1.5 % agarose gels, PCR products of COI gene were sent to a commercial company for sequencing (Invitrogen Biotechnology Co., Ltd., Shanghai, China). The ITS products were purified using the PCR Purification Kit (Aidlab Biotechnologies Co., Ltd., Beijing, China), ligated to pMD19-T Vector (TaKaRa Biotechnology, Dalian, China) and then transformed into competent cell of Escherichia coli DH5α (TaKaRa Biotechnology, Dalian, China). Finally, a positive clone per individual was sequenced in both directions using M13 and M17 primers by Invitrogen Biotechnology Co., Ltd. (Shanghai China).
Sequence variation and population genetic analysis
All sequences in Siphonaria japonica were edited by comparing both strands using DNAMAN 7 software (LynnonBioSoft, Quebec, Canada), and then aligned with MUSCLE  using MEGA 5  with default settings. Standard molecular diversity indices including haplotype diversity (h) and nucleotide diversity (π), and neutrality test including Fu’s F s  and Tajima’s D  were calculated using ARLEQUIN 3.5 .
jModelTest 2.1.1  was utilized to estimate the best-fitting substitution model and substitution parameters with the Bayesian information criterion (BIC) for each species (Table 3). For sequences from S. japonica, phylogenetic analyses were performed based on the neighbour-joining (NJ) approach using MEGA incorporating the Tamura-Nei model (TrN)  with corresponding gamma correction for COI and ITS (the closest model in MEGA to TVM + G is the TrN + G model for NJ tree building). 1 000 bootstrap replicates were carried out to assess the clade credibility of the NJ phylogram. A haplotype network illustrating genealogical relationships between haplotypes was constructed for each species using the median-joining (MJ) method with Network 4.6 .
Pairwise ΦST measures were calculated to evaluate the levels of genetic differentiation and significance was estimated with 10 000 permutations using ARLEQUIN. Furthermore, hierarchical analysis of molecular variance (AMOVA)  was performed to test for possible phylogeographic separation. In the present study, genetic differentiation among locations was only analyzed in S. japonica, because the published data from other three species have been analyzed. Because this study focused on the gene flow and phylogeographic break, the locations of each species were divided into three groups: the Yellow Sea (YS) group, the East China Sea (ECS) group and the South China Sea (SCS) group based on the geographical locations (refer to Table 4 for the allocation of sites). Because no significantly genetic differentiation was observed between ECS and SCS (see the Results Table 3), the ECS and SCS were combined as the ESCS group to compare with the YS group to test for the hypothesis of reduced gene flow across the Yangtze River outflow. On the other hand, samples from all locations for each species were considered as a single group to verify the significance of partitioning of genetic variance among all samples.
The coalescent-based approach is widely used to reconstruct the demographic history. In the present study, the Bayesian skyline plot (BSP) was generated in BEAST v1.7.4  to estimate the change in effective population size over time. The analyses were only performed with the mitochondrial dataset because no distinct clade was revealed in ITS and the mutation rate is unknown for the marker, hampering the estimate of temporal scale. Analyses were conducted under corresponding model suggested by jModelTest for each species with constant Bayesian skyline tree priors with 10 groups under a strict clock model. Default priors were used for parameter settings. Three independent MCMCMC searches were run for 100 million generations and parameters were recorded every 10 000 generations with the first 10 million generations discarded as burn-in. For each run, the effective sample sizes (ESS) of important parameters sampled from the MCMCMC were >1000 in all three replicate runs calculated by TRACER 1.4 . Results from three runs were then pooled with LOGCOMBINER version 1.7  for final reconstruction of the BSP. Pairwise mismatch distribution analysis was also performed using the mitochondrial dataset of each species with ARLEQUIN 3.5  to validate the result from BSP. Population expansion parameters (τ) were transformed to estimates of real-time since expansion (t) with the formula τ = 2μkt, where μ is the mutation rate and k is the sequence length. As two geographically-partitioned groups were revealed in phylogenetic analysis of COI in S. japonica (Fig. 2a), the demographic history of each putative group was reconstructed using both methods above.
Because of the absence of clear fossil or geological record, Colgan & Costa  and Teske et al.  had used a mutation rate of 1 % per million years to estimate divergence times in Siphonaria genus, when referring to available calibrated fossil data for marine gastropod. Thus, in this study, such a mutation rate (1 % Myr−1) with a generation of 1 year was applied for COI of S. japonica. A divergence rate of 0.85-1.15 % Myr−1 for COI was set in C. toreuma, as used for estimation times of the C. nigrolineata . The divergence rate was set at 2.6-4.1 % Myr−1 for COIII of S. horneri  and 0.7-2.4 % Myr−1 for COI of A. pectinata , respectively, which were proposed in the original papers. In addition a generation of 1 year for C. toreuma  and S. horneri  and a generation of 2 years for A. pectinata  were assumed. In the present study, the mutation rate for these three species (C. toreuma, 0.5 % Myr−1; S. horneri, 1.675 % Myr−1; A. pectinata, 0.775 % Myr−1) was obtained by averaging the divergence rates used above and then dividing by two.
To evaluate whether the BSP was the best model in reconstructing demographic histories, BSP model was compared with two simple moles in Tracer V 1.6: constant population size and expansion growth. The two alternative models were run in BEAST using the same way as mentioned above for BSP. A posterior simulation-based analogue of Akaike’s information criterion through MCMC (AICM)  was used to compare all three models, which measured AIC from the posterior of each model, with score > 10 as strong evidence in favor of one model over the others .
Gene flow analysis
To examine the level and direction of contemporary gene flow across the phylogeographic break, we adopted the coalescent-based approach implemented in MIGRATE-N 3.5.1  for estimating the migration rates between groups within the four studied species. Locations were grouped into a “YS” group and an “ESCS” group according to their geographical distribution. Random sub-samples were performed to allow the two groups to contain comparable number of individuals. The Bayesian approach was utilized to infer the mutation-scaled effective population size (Θ = 2Nμ, with N = effective population size and μ = mutation rate) and the mutation-scaled effective immigration rate (Μ = m/μ, with m = immigration rate) . Analyses were conducted with a full migration matrix model (Θ and Μ were estimated jointly from the data). The effect number of migrates per generation (Nm) among groups can be calculated by multiplying Θ and Μ together. For the Bayesian approach, a single long chain with slice sampling for the proposal distribution was used according to the recommendation of the author . We performed Migrate with the DNA sequence model. Initially, short runs were performed to estimate Θ and Μ with F ST with a uniform prior for parameters, followed by subsequent runs set using parameters estimated in the short runs. Five independent sets of runs were conducted, each containing one long chain of 5 000 000 steps with a burn-in time of 500 000, a sampling increment of 1000, and an adaptive heating scheme with four chains and temperatures of 1.0, 1.5, 3.0, and 10 000. The results of the independent runs were congruent and the last run was chosen for interpretations.
Availability of supporting data
The newly obtained DNA sequences: GenBank accessions: KF716505-KF716747 mtDNA; KF716748-KF716980 nuclear gene for Siphonaria japonica. The reanalyzed data: Cellana toreuma (COI, GenBank accession number, JQ313140-JQ313557; ), Sargassum horneri (COIII, GenBank accession number, JF461002-JF461052; ), Atrina pectinata (COI, GenBank accession number, HQ449254-HQ449388; ). All the alignments data supporting the results of this article are available in the Dyad data repository under doi: 10.5061/dryad.2mv73 .
We would like to thank Guo-dong Han for his assistance in the laboratory, Dr. Colin Little for his great helps in preparing the manuscript. This work was supported by grants from National Natural Science Foundation of China (41476115 41276126), National Basic Research Program of China (2013CB956504), Nature Science funds for Distinguished Young Scholars of Fujian Province, China (2011 J06017), the Fundamental Research Funds for the Central Universities, Program for New Century Excellent Talents in University of Fujian Province and Program for New Century Excellent Talents in University.
- Imbrie J, Boyle EA, Clemens SC, Duffy A, Howard WR, Kukla G, et al. On the structure and origin of major glaciation cycles 1. Linear responses to milankovitch forcing. Paleoceanography. 1992;7:701–38.View ArticleGoogle Scholar
- Imbrie J, Berger A, Boyle EA, Clemens SC, Duffy A, Howard WR, et al. On the structure and origin of major glaciation cycles 2. The 100, 000 Year cycle. Paleoceanography. 1993;8:699–735.View ArticleGoogle Scholar
- Marko PB, Hoffman JM, Emme SA, McGovern TM, Keever CC, Nicole Cox L. The ‘Expansion-Contraction’ model of Pleistocene biogeography: rocky shores suffer a sea change? Mol Ecol. 2010;19:146–69.PubMedView ArticleGoogle Scholar
- Hellberg ME. Gene flow and isolation among populations of marine animals. Annu Rev Ecol Evol S. 2009;40:291–310.View ArticleGoogle Scholar
- Rocha LA, Bass AL, Robertson G, Bowen BW: Adult habitat preferences, larval dispersal, and the comparative phylogeography of three Atlantic surgeonfishes (Teleoste: Acanthuridae). Mol Ecol Notes 2002, 11:243-251.Google Scholar
- Marko P. ‘What's larvae got to do with it?’Disparate patterns of post‐glacial population structure in two benthic marine gastropods with identical dispersal potential. Mol Ecol. 2004;13:597–611.PubMedView ArticleGoogle Scholar
- Lourie S, Green D, Vincent A. Dispersal, habitat differences, and comparative phylogeography of Southeast Asian seahorses (Syngnathidae: Hippocampus). Mol Ecol. 2005;14:1073–94.PubMedView ArticleGoogle Scholar
- Nakano T, Sasaki T, Kase T. Color polymorphism and historical biogeography in the Japanese patellogastropod limpet Cellana nigrolineata (Reeve) (Patellogastropoda: Nacellidae). Zoolog Sci. 2010;27:811–20.PubMedView ArticleGoogle Scholar
- Teske PR, Papadopoulos I, Mmonwa KL, Matumba T, McQUAID CD, Barker NP, et al. Climate‐driven genetic divergence of limpets with different life histories across a southeast African marine biogeographic disjunction: different processes, same outcome. Mol Ecol. 2011;20:5025–41.PubMedView ArticleGoogle Scholar
- Li W, Zhong J, Wang Y. Genetic diversity and population structure of two lancelets along the coast of China. Zoolog Sci. 2013;30:83–91.PubMedView ArticleGoogle Scholar
- Scheltema RS. Larval dispersal as a means of genetic exchange between geographically separated populations of shallow-water benthic marine gastropods. BioBull. 1971;140:284–322.View ArticleGoogle Scholar
- Kyle CJ, Boulding EG. Comparative population genetic structure of marine gastropods (Littorina spp) with and without pelagic larval dispersal. Mar Biol. 2000;127:835–45.View ArticleGoogle Scholar
- Sponaugle S, Cowen RK, Shanks A, Morgan SG, Leis JM, Pineda J, et al. Predicition self-recruitment in marine populations: biophysical correlates and mechanisms. B Mar Sci. 2002;70:341–75.Google Scholar
- Cowen RK, Sponaugle S. Larval dispersal and marine population connectivity. Ann Rev Mar Sci. 2009;1:443–66.PubMedView ArticleGoogle Scholar
- Miller AD, Versace VL, Matthews TG, Montgomery S, Bowie KC. Ocean currents influence the genetic structure of an intertidal mollusc in southeastern Australia–implications for predicting the movement of passive dispersers across a marine biogeographic barrier. Ecol Evol. 2013;3:1248–61.PubMed CentralPubMedView ArticleGoogle Scholar
- Whitfield AK, Wooldridge TH: Changes in freshwater supplies to southern African estuaries: some theoretical and practical considerations. In: Changes in fluxes in estuaries: implications from science to management. Dyer KR, Orth RJ, editors. Olsen & Olsen, Fredensborg, Denmark. 1994: p.41-50.Google Scholar
- Shanks AL, Eckert GL. Population persistence of California current fishes and benthic crustaceans: a marine drift paradox. Ecol Monogr. 2005;75:505–24.View ArticleGoogle Scholar
- Cheang CC, Chu KH, Ang Jr PO. Phylogeography of the marine macroalga Sargassum hemiphyllum (Phaeophyceae, Heterokontophyta) in northwestern Pacific. Mol Ecol. 2010;19:2933–48.PubMedView ArticleGoogle Scholar
- Kelly RP, Palumbi SR. Genetic structure among 50 species of the northeastern pacific rocky intertidal community. PLoS One. 2010;5:e8594.PubMed CentralPubMedView ArticleGoogle Scholar
- Dong YW, Wang HS, Han GD, Ke CH, Zhan X, Nakano T, et al. The impact of Yangtze River discharge, ocean currents and historical events on the biogeographic pattern of Cellana toreuma along the China coast. PLoS One. 2012;7:e36178.PubMed CentralPubMedView ArticleGoogle Scholar
- Wang PX. Response of western Pacific marginal seas to glacial cycles: paleoceanographic and sedimentological features. Mar Geol. 1999;156:5–39.View ArticleGoogle Scholar
- Li YL, Kong XY, Yu ZN, Kong J, Ma S, Chen LM. Genetic diversity and historical demography of Chinese shrimp feneropenaeus chinensis in yellow Sea and bohai Sea based on mitochondrial DNA analysis. Afr J Biotechnol. 2009;8:1193–202.Google Scholar
- Xu J, Chan TY, Tsang LM, Chu KH. Phylogeography of the mitten crab Eriocheir sensu stricto in East Asia: Pleistocene isolation, population expansion and secondary contact. Mol Phylogenet Evol. 2009;52:45–56.PubMedView ArticleGoogle Scholar
- Ni G, Li Q, Kong L, Yu H. Comparative phylogeography in marginal seas of the northwestern Pacific. Mol Ecol. 2014;23:534–48.PubMedView ArticleGoogle Scholar
- Fang J, Li Y, Sun S, Deng J. Analysis of runoff change characteristics at Datong station of Yangtze River. Water Res Power. 2011;29:9–12.Google Scholar
- Beardsley RC, Limeburner R, Yu H, Cannon GA. Discharge of the Changjiang (Yangtze River) into the East China Sea. Con Shelf Res. 1985;4:57.View ArticleGoogle Scholar
- Senjyu T, Enomoto H, Matsuno T, Matsui S. Interannual salinity variations in the Tsushima Strait and its relation to the Changjiang discharge. J Oceanogr. 2006;62:681–92.View ArticleGoogle Scholar
- Chu P, Chen YC, Kuninaka A. Seasonal variability of the Yellow Sea/East China Sea surface fluxes and thermohaline structure. Adv Atmos Sci. 2005;22:1–20.View ArticleGoogle Scholar
- Wu Z-Y, Chen H-X, Liu N. Relationship between east China Sea Kuroshio and climatic elements in East China. B Mar Sci. 2010;12:1–9.Google Scholar
- Feng S, Li F, Li S: An introduction to marine science. China Higher Education Press, Beijing, China; 1999.Google Scholar
- Coleman JM, Roberts HH. Deltaic coastal wetlands. Geol Mijnbouw. 1989;68:1–24.Google Scholar
- Orton GJ, Reading HG. Variability of deltaic processes in terms of sediment supply, with particular emphasis on grain-size. Sedimentology. 1993;40:475–512.View ArticleGoogle Scholar
- Zhu ZW, Gao MS, Zhu YF. Quantitative analysis of basic types of coast zones and their distribution. Front Earth Sci. 2008;15:315–21.Google Scholar
- Zhao Y, Li Q, Kong L, Mao Y. Genetic and morphological variation in the venus clam Cyclina sinensis along the coast of China. Hydrobiologia. 2009;635:227–35.View ArticleGoogle Scholar
- Cheang CC, Chu KH, Ang Jr PO. Morphological and genetic variation in the populations of Sargassum Hemiphyllum (Phaeophyceae) in the northwestern Pacific. J Phycol. 2008;44:855–65.View ArticleGoogle Scholar
- Hu ZM, Uwai S, Yu SH, Komatsu T, Ajisaka T, Duan DL. Phylogeographic heterogeneity of the brown macroalga Sargassum horneri (Fucaceae) in the northwestern Pacific in relation to late Pleistocene glaciation and tectonic configurations. Mol Ecol. 2011;20:3894–909.PubMedView ArticleGoogle Scholar
- Liu J, Li Q, Kong L, Zheng X. Cryptic diversity in the pen shell Atrina pectinata (Bivalvia: Pinnidae): high divergence and hybridization revealed by molecular and morphological data. Mol Ecol. 2011;20:4332–45.PubMedView ArticleGoogle Scholar
- Hirano Y. Siphonaria (pulmonate limpet) survey of Japan. II. Periodicity of spawning activity in Siphonaria japonica. Publications of the Seto Marine Biological Laboratory. 1980;25:335–42.Google Scholar
- Wang Z, Wu C. Study on the age and growth of Cellauma toreuma (Reeve) in north Zhejiang coastal area. J Zhejiang Ocean University (Nat Sci). 1999;19:316–23.Google Scholar
- Pang SJ, Liu F, Shan TF, Gao SQ, Zhang ZH. Cultivation of the brown alga Sargassum horneri: sexual reproduction and seedling production in tank culture under reduced solar irradiance in ambient temperature. J Appl Phycol. 2009;21:413–22.View ArticleGoogle Scholar
- Colgan DJ, da Costa P. Possible drivers of biodiversity generation in the Siphonaria of southeastern Australia. Marine Biodiversity. 2012;43:73–85.View ArticleGoogle Scholar
- Lisiecki LE, Raymo ME: A Pliocene‐Pleistocene stack of 57 globally distributed benthic δ18O records. Paleoceanography 2005, 20:1-17.Google Scholar
- Domitsu H, Uchida J, Ogane K, Dobuchi N, Sato T, Ikehara M, et al. Stratigraphic relationships between the last occurrence of Neogloboquadrina inglei and marine isotope stages in the northwest Pacific, D/V Chikyu Expedition 902, Hole C9001C. Newsl Stratigr. 2011;44:113–22.View ArticleGoogle Scholar
- Yang J, Li Q, Kong L, Zheng X, Wang R. Genetic structure of the veined rapa whelk (Rapana venosa) populations along the coast of China. Biochem Genet. 2008;46:539–48.PubMedView ArticleGoogle Scholar
- Zheng WJ, Zhu SH, Shen XQ, Liu BQ, Pan ZC, Ye YF. Genetic differentiation of Tegillarca granosa based on mitochondrial COI gene sequences. Zoological Res. 2009;30:17–23.View ArticleGoogle Scholar
- Ni G, Li Q, Kong L, Zheng X. Phylogeography of bivalve Cyclina sinensis: testing the historical glaciations and Changjiang River outflow hypotheses in northwestern Pacific. PLoS One. 2012;7:e49487.PubMed CentralPubMedView ArticleGoogle Scholar
- Voris HK. Maps of Pleistocene sea levels in Southeast Asia: shorelines, river systems and time durations. J Biogeogr. 2000;27:1153–67.View ArticleGoogle Scholar
- Lambeck K, Esat TM, Potter EK. Links between climate and sea levels for the past three million years. Nature. 2002;419:199–206.PubMedView ArticleGoogle Scholar
- Wang PX, Sun XJ. Last glacial maximum in China: comparison between land and sea. Catena. 1994;23:341–53.View ArticleGoogle Scholar
- Kimura M. Paleogeography of the Ryukyu Islands. Tropics. 2000;10:5–24.View ArticleGoogle Scholar
- Qiao FL: Regional oceanography of China seas: physical oceanography. China Ocean Press, Beijing, China; 2012.Google Scholar
- Liu JP, Xu KH, Li AC, Milliman JD, Velozzi DM, Xiao SB, et al. Flux and fate of Yangtze river sediment delivered to the East China Sea. Geomorphology. 2007;85:208–24.View ArticleGoogle Scholar
- Zardus JD, Hadfield MG. Multiple origins and incursions of the Atlantic barnacle Chthamalus proteus in the Pacific. Mol Ecol. 2005;14:3719–33.PubMedView ArticleGoogle Scholar
- Bird CE, Holland BS, Bowen BW, Toonen RJ. Contrasting phylogeography in three endemic Hawaiian limpets (Cellana spp.) with similar life histories. Mol Ecol. 2007;16:3173–86.PubMedView ArticleGoogle Scholar
- Yoshida T. Studies on the distribution and drift of the floating seaweeds. Bull Tohoku Reg Fish Res Lab. 1963;23:141–86.Google Scholar
- Ohno M: Observation on the floating seaweeds of near-shore waters of southern Japan. Hydrobiologia. 1984;116:408-412.Google Scholar
- Komatsu T, Tatsukawa K, Filippi JB, Sagawa T, Matsunaga D, Mikami A, et al. Distribution of drifting seaweeds in eastern East China Sea. J Mar Syst. 2007;67:245–52.View ArticleGoogle Scholar
- Lee KM, Yang EC, Coyer JA, Zuccarello GC, Wang WL, Choi CG, et al. Phylogeography of the seaweed Ishige okamurae (Phaeophyceae): evidence for glacial refugia in the northwest Pacific region. Mar Biol. 2012;159:1021–8.View ArticleGoogle Scholar
- Cheang CC, Tsang LM, Ng WC, Williams GA, Chu KH, Chan BK. Phylogeography of the cold water barnacle Chthamalus challengeri in the north western Pacific: effect of past population expansion and contemporary gene flow. J Biogeogr. 2012;39:1819–35.View ArticleGoogle Scholar
- Toews DP, Brelsford A. The biogeography of mitochondrial and nuclear discordance in animals. Mol Ecol. 2012;21:3907–30.PubMedView ArticleGoogle Scholar
- Beerli P: Migrate documentation, Version 3.2.1, Available from http://popgen.sc.fsu.edu/Migrate/Download.html; 2012.Google Scholar
- Barry JP, Baxter CH, Sagarin RD, Gilman SE. Climate-related, long-term faunal changes in a California rocky intertidal community. Science(Washington). 1995;267:672–5.View ArticleGoogle Scholar
- Hiscock K, Southward A, Tittley I, Hawkins S. Effects of changing temperature on benthic marine life in Britain and Ireland. Aquat Conserv. 2004;14:333–62.View ArticleGoogle Scholar
- Mieszkowska N, Kendall M, Hawkins S, Leaper R, Williamson P, Hardman-Mountford N, et al. Changes in the range of some common rocky shore species in Britain–a response to climate change? Hydrobiologia. 2006;555:241–51.View ArticleGoogle Scholar
- Hawkins SJ, Moore PJ, Burrows MT, Poloczanska E, Mieszkowska N, Herbert RJH,J. Complex interactions in a rapidly changing world: responses of rocky shore communities to recent climate change. Clim Res. 2008;37:123–33.View ArticleGoogle Scholar
- Lima FP, Wethey DS. Three decades of high-resolution coastal sea surface temperatures reveal more than warming. Nat Commun. 2012;3:307.View ArticleGoogle Scholar
- Bao B, Ren G. Climatological characteristics and long-term change of SST over the marginal seas of China. Cont Shelf Res. 2014;77:96–106.View ArticleGoogle Scholar
- Sun S. A discussion on management of coastl reclamation. Ocean Develop Manage. 2004;6:21–3.Google Scholar
- Yang SL, Zhang J, Zhu J, Smith JP, Dai SB, Gao A, et al.: Impact of dams on Yangtze River sediment supply to the sea and delta intertidal wetland response. J Geophys Res Earth. 2005;110:1-12.Google Scholar
- 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–9.PubMedGoogle Scholar
- Oliverio M, Mariottini P. Contrasting morphological and molecular variation in Coralliophila meyendorffii (Muricidae, Coralliophilinae). J Mollus Stud. 2001;67:243–5.View ArticleGoogle Scholar
- Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004;32:1792–7.PubMed CentralPubMedView ArticleGoogle Scholar
- Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S. MEGA5: Molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011;28:2731–9.PubMed CentralPubMedView ArticleGoogle Scholar
- Fu YX. Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection. Genetics. 1997;147:915–25.PubMed CentralPubMedGoogle Scholar
- Tajima F. Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics. 1989;123:585–95.PubMed CentralPubMedGoogle Scholar
- Excoffier L, Lischer HE. 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–7.PubMedView ArticleGoogle Scholar
- Darriba D, Taboada GL, Doallo R, Posada D. jModelTest 2: more models, new heuristics and parallel computing. Nat Methods. 2012;9:772.PubMedView ArticleGoogle Scholar
- Tamura K, Nei M. Estimation of the number of nucleotide substitutions in the control region of mitochondrial DNA in humans and chimpanzees. Mol Biol Evol. 1993;10:512–26.PubMedGoogle Scholar
- Bandelt HJ, Forster P, Rohl A. Median-joining networks for inferring intraspecific phylogenies. Mol Biol Evol. 1999;16:37–48.PubMedView ArticleGoogle Scholar
- Excoffier L, Smouse PE, Quattro JM. Analysis of molecular variance inferred from metric distances among DNA haplotypes: application to human mitochondrial DNA restriction data. Genetics. 1992;131:479–91.PubMed CentralPubMedGoogle Scholar
- Drummond AJ, Suchard MA, Xie D, Rambaut A. Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol Biol Evol. 2012;29:1969–73.PubMed CentralPubMedView ArticleGoogle Scholar
- Ajisaka T. On the fluctuation in maturation period of Sargassum population in tropical area. Gekkan Kaiyo. 1997;29:462–6.Google Scholar
- Wang ZR: Fauna sinica: Phylum mollusca. Order mytiloida. Chinese Science Press, Beijing, China; 1997.Google Scholar
- Raftery A, Newton M, Satagopan J, Krivitsky P. Estimating the integrated likelihood via posterior simulation using the harmonic mean identity. In: Bayesian statistics. Bernardo JM, Bayarri MJ, Berger JO, editors. Oxford: Oxford University Press; 2007. p. 1–45.Google Scholar
- Baele G, Lemey P, Bedford T, Rambaut A, Suchard MA, Alekseyenko AV. Improving the accuracy of demographic and molecular clock model comparison while accommodating phylogenetic uncertainty. Mol Biol Evol. 2012;29:2157–67.PubMed CentralPubMedView ArticleGoogle Scholar
- Wang J, Tsang LM, Dong YW. Causations of phylogeographic barrier of some rocky shore species along the Chinese coastline. 2015. Dryad. http://dx.doi.org/10.5061/dryad.2mv73.
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.