Homogeneous Population of the Brown Alga Sargassum polycystum in Southeast Asia: Possible Role of Recent Expansion and Asexual Propagation

Sze Wai Chan; Chi Chiu Cheang; Anong Chirapart; Grevo Gerung; Chea Tharith; Put Ang

doi:10.1371/journal.pone.0077662

Abstract

Southeast Asia has been known as one of the biodiversity hotspots in the world. Repeated glacial cycles during Pleistocene were believed to cause isolation of marine taxa in refugia, resulting in diversification among lineages. Recently, ocean current was also found to be another factor affecting gene flow by restricting larval dispersal in animals. Macroalgae are unique in having mode of reproduction that differs from that of animals. Our study on the phylogeographical pattern of the brown macroalga Sargassum polycystum using nuclear Internal Transcribed Spacer 2 (ITS2), plastidal RuBisCO spacer (Rub spacer) and mitochondrial cytochrome oxidase subunit-III (Cox3) as molecular markers revealed genetic homogeneity across 27 sites in Southeast Asia and western Pacific, in sharp contrast to that revealed from most animal studies. Our data suggested that S. polycystum persisted in single refugium during Pleistocene in a panmixia pattern. Expansion occurred more recently after the Last Glacial Maximum and recolonization of the newly flooded Sunda Shelf could have involved asexual propagation of the species. High dispersal ability through floating fronds carrying developing germlings may also contribute to the low genetic diversity of the species.

Citation: Chan SW, Cheang CC, Chirapart A, Gerung G, Tharith C, Ang P (2013) Homogeneous Population of the Brown Alga Sargassum polycystum in Southeast Asia: Possible Role of Recent Expansion and Asexual Propagation. PLoS ONE 8(10): e77662. https://doi.org/10.1371/journal.pone.0077662

Editor: Valerio Ketmaier, Institute of Biochemistry and Biology, Germany

Received: July 17, 2013; Accepted: September 12, 2013; Published: October 17, 2013

Copyright: © 2013 chan et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Funding: This work was supported by The Research Grant Council of Hong Kong. The funder had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Competing interests: The authors have declared that no competing interests exist.

Introduction

Southeast Asia is well known to be a hotspot for biodiversity. This high biodiversity is believed to have been contributed by its complex geological history. The repeated glacial cycles during Pleistocene caused fluctuation of sea level. This resulted in the formation of exposed land (e.g. Sunda Shelf) in relatively shallow sea that connected adjacent islands to form land bridges [1]. Pacific and Indian Oceans were separated and this separation was suggested to be the main event causing species differentiation between the two Oceans. Meanwhile, isolated basins such as the South China Sea, Sulu Sea, Celebes Sea and Flores Sea served as refugia for marine taxa. Their isolation led to the diversification of lineages among populations [2–4]. Recent oceanic currents, such as the Indonesian Thoughflow, is found to be restricting the larval dispersal of marine animals, such as the stomatopods, between Java Sea and Flores Sea [3]. The high genetic diversity in marine taxa in Southeast Asia is thus caused not only by past isolations, but also by restricted larval dispersal by recent oceanic conditions.

Most of the marine phylogeographical studies in Southeast Asia, however, were on animals. Seaweeds, which exhibit a largely different dispersal mode than the animals, may provide new insights in understanding the phylogeographical processes in this region. Oceanic circulation has a profound effect on much of the population structure of marine animals as genetic connectively is highly associated with their pelagic larval duration. Unlike animals, seaweeds produce immobile germlings after fertilization and dispersal distance is expected to be much shorter. This may lead to higher genetic differentiation among seaweed populations. In North Atlantic, extensive phylogeographical studies have been carried out using multilocus molecular markers on brown algae [5,6]. Not much can be said of such studies in Southeast Asia except for a recent published research on the red alga Gracilaria changii in Peninsular Malaysia [7].

There are currently more than 300 species of Sargassum [8] reported from Southeast Asia. This makes Sargassum an ideal tool for biogeographical and phylogeographical studies because of its wide distribution range and high species diversity. The dispersal mode of the genus Sargassum is also unique among macroalgae. While they have short dispersal distance (<1m) of their germlings [9], they can be dispersed in the form of drifting fronds carrying with them developing propagules [10]. With the frequent isolation of many local basins in Southeast Asia in the geological past and the recent existence of complex oceanic circulation that is associated with monsoon and directional current flows (e.g. Indonesian Thoughflow), it is worth investigating the phylogeography of this region using this unique group of marine organisms. The target species in this study, Sargassum polycystum C. Agardh has a wide distribution range over the Indo-Pacific region, and a large distribution record in Southeast Asia [11,12]. This study has therefore the following objectives: (1) to investigate the genetic diversity and population structure of S. polycystum, with focus in Southeast Asia; and (2), to evaluate the potential mechanisms involved in bringing about any genetic structure observed. Three molecular markers were used to achieve these objectives, including the nuclear Internal Transcribed Spacer 2 (ITS2), plastidal RuBisCo spacer (Rub spacer) and mitochondrial cytochrome oxidase subunit-III (Cox3).

Materials and Methods

Sample collection and DNA extraction

Field studies did not involve any endangered or protected species. This study also involved collaborative works between institutions so no special permits were needed for sampling. Samples were collected from 22 sampling sites either by snorkeling or sampling during low tide (Figure 1; see also Table S1 in File S1). Individuals were collected more than 1m apart to avoid having the same mother plant as Sargassum germlings have short dispersal distance [9]. Leaf tips of 3-5 cm were dried and stored in silica gel for molecular analysis. The voucher plants were air-dried and deposited at the Simon F.S. Li Marine Science Laboratory Herbarium, The Chinese University of Hong Kong. Genomic DNA was extracted by modified cetyltrimethylammonium bromide (CTAB) method [2]. It was further purified by GENECLEAN II kit (Obiogene Inc.), following the manufacturer’s instructions.

Download:

Figure 1. Haplotype distribution and haplotype network of Sargassum polycystum for ITS2, Rub spacer and Cox3.

Pie chart size is proportional to sample size. Abbreviations for sample sites are given in Table S1 in File S1. Dominant currents are shown in solid lines [13] and seasonally reversing currents in dashed lines [14]. Light gray area shows the coastal outline during Pleistocene maximum low sea level of 120m [1]. SEC: South Equatorial Current, NEC: North Equatorial Current, NECC: North Equatorial Counter Current, NGCC: New Guinea Coastal Current, ME: Mindanao Eddy, HE: Halmahera Eddy.

https://doi.org/10.1371/journal.pone.0077662.g001

Polymerase Chain Reaction (PCR) and Sequencing

Primers for the three markers and PCR profiles were listed in [15]. All PCR products were sequenced in both directions with the same primers by Macrogen (Macrogen Inc., Seoul, Korea) employing the BigDye^TM terminator method. The presence of multiple peaks in chromatograms was checked in both directions to see whether there were intragenomic variations. Sequences were aligned and edited together with the sequences of S. polycystum available from GenBank by naked eyes using MEGA ver. 5 [16].

Phylogeny and haplotype network construction

To ensure all the samples were S. polycystum with no misidentification before subsequent analysis, representative sequences from each haplotype of combined ITS2, Rub spacer and Cox3 sequences were initially screened to ensure they formed a monophyletic lineage and single taxonomic unit. The sequences were compared with those of other species in the subgenus Sargassum [15] using Bayesian Inference (BI) and Maximum Likelihood (ML) methods by MRBAYES 3.2.1 [17] and PhyML 3.0 [18] respectively. The Corrected Akaike information criterion (AICc) implemented in jMODELTEST 2.1.1 [19] was used to determine the best-fit substitution model. Turbinaria ornata, as suggested by Stiger and Horiguchi [20] and Phillips et al. [21], was used as outgroup. Any sequence that was not monophyletic was discarded to ensure consistency in subsequent analysis.

To investigate the phylogenetic relationship within S. polycystum, Maximum Likelihood (ML) using PhyML 3.0 and Bayesian Inference (BI) using MRBAYES 3.2.1were applied. Best-fit substitution models were determined by jMODELTEST 2.1.1. Sargassum ilicifolium was used as outgroup as it is closely related within the same subgenus Sargassum [15]. Significance of the branching was evaluated by 1000 bootstrap in ML and posterior probability in BI analysis. For BI, two independent Markov Chain Monte Carlo (MCMC) searches were conducted until the divergence became small and stationary. Trees were sampled every 100 cycles in 1 500 000 generations and first 25% of the trees were discarded as burn-in. Haplotype networks for the three markers were generated to estimate gene genealogies using TCS 1.21 [22].

Population structure analysis

Population structure of S. polycystum was analyzed by ARLEQUIN 3.5.1.2 [23]. Haplotype (h) and nucleotide (π) diversity [24,25] were calculated for each site. Tajima’s D [26] and Fu’s Fs [27] were calculated to test for the selective neutrality of the markers. Tajima’s D can distinguish between a recent population bottleneck and population expansion. Fu’s Fs is more sensitive to recent population expansion than Tajima’s D which gives large Fs negative values for excess number of alleles. Pairwise Φ_ST was calculated to determine the level of gene flow between populations. Hierarchical Analysis of Molecular Variance (AMOVA) with 10,000 permutations was conducted based on Φ_ST to find out the best spatial groupings of populations. Several groupings were tested iteratively, based on geographical proximity of the collecting sites, possible connectivity by current flow and haplotype similarity, to investigate whether there were any biogeographical barriers in between groups or localized haplotypes within groups that shaped the population structure. Each iteration was repeatedly tested by AMOVA until the highest percentages of variances between groups could be obtained. In addition, sampling sites were divided into smaller subregions corresponding to known marine basins or putative refugia, in order to investigate the direction and amount of gene flow between them. These subregions included South China Sea and Gulf of Thailand (BL, PP, PA, HN, NT, OB, SS, MS; see Figure 1 for site location), West Coast of Malay Peninsula (PC, PD, PS, PDo), West Java (WJ, SE, MB, PK, TE, PR), Celebes Sea and Flores Sea (IC, MA, KE), Guam (CC, PG) and Pacific Islands (SO, VA, FI, WF). All haplotypes in the populations were analyzed using the software MIGRATE 3.6 [28,29]. For each molecular marker, Bayesian search strategy was used in single 1,000,000 step chain with sampling increment of 100. Default setting was used, assuming unrestricted dispersal, to estimate number of immigrants per generation as calculated by multiplying θ (mutation-scaled effective population size) with M (mutation-scaled effective immigration rate) [28,29]. Samples of initial 25% were discarded as burn-in. Two independent runs were conducted and the results combined to represent >10,000 effective population size.

A Mantel test [30] with 10,000 permutations was applied to test for the isolation by distance model by comparing the pairwise Φ_ST and the matrix of geographical distances among sampling sites. Geographical distances were measured using the shortest distance between two sites by Google Earth (http://www.google.com/earth/index.html).

Demographic history

Demographic history was only conducted for Cox3 data as there was no molecular clock available for ITS2 and Rub spacer. Mismatch distribution was used to detect any recent population expansion which is indicated with a unimodel distribution [31] using ARLEQUIN 3.5.1.2. Time of expansion (in generation), t, can be estimated by τ= 2μt. τ is the crest of mismatch distribution, μ is the mutation rate of the marker per generation. In addition, mismatch distribution was also conducted in ITS2 and Rub spacer without calculating the time of expansion. A Bayesian Skyline Plot of Cox3 was generated with BEAST 1.5.3 [32] using Markov Chain Monte Carlo (MCMC) sampling procedures with 10⁸ steps, sampled at every 1000 steps. Runs were repeated until effective sample size > 200 was reached in all parameters, as recommended in the user manual [33]. All runs were pooled in LogCombiner 1.7.3 with first 10% of generations discarded as burn-in. The Bayesian Skyline Plot was generated and visualized in Tracer 1.5.

The mutation rate of Cox3 of Sargassum was calculated based on the psbA clock (0.08%-0.12%/Myr) calibrated for Fucus [5]. PsbA sequences of S. ilicifolium were used as this species is closely related to S. polycystum within the same subgenus [15]. Divergence of psbA sequences between S. ilicifolium and F. vesiculosus (Genbank accession no: DQ307679) calculated was 5.67%. Therefore the estimated divergence time was 47.25 to 79.88 Myr (5.67/0.12 and 5.67/0.0.8 respectively). Divergence of Cox3 sequences between S. ilicifolium and F. vesiculosus (AY494079) was 14.7% and the estimated divergence rate of Cox3 was 0.207 to 0.311%/Myr (14.7/79.88 and 14.7/47.25 respectively). As the mutation rate should be half of the divergence rate, the estimated mutation rate was 1.035×10^-9 to 1.555×10^-9 substitutions per site per year. The generation time of 1 year was used as most Sargassum species have an annual cycle of growth and regeneration [34] and S. polycystum was shown to have a unimodal growth cycle [35].

Results

Genetic diversity

A total of 261, 191 and 185 sequences were obtained from 27 sampling sites and Genbank for ITS2 (294 bp), Rub spacer (163 bp) and Cox3 (379 bp) respectively (Table S1 in File S1). Cox3 showed the highest genetic variability with highest haplotype and nucleotide diversity (H = 0.40 ± 0.04, π = 0.089 ± 0.082). Rub spacer had the lowest variability (H = 0.18 ± 0.04, π = 0.079 ± 0.101). There were 9, 4, 5 haplotypes for ITS2, Rub spacer and Cox3 respectively with 13, 3, 6 polymorphic sites. No intragenomic variations were observed for all the markers as all the chromatographs showed clear and consistent peaks in both sequencing directions.

Interspecific and intraspecific phylogenetic relationship

HKY+G model determined by jMODELTEST was used to analyze the taxonomic status of haplotypes from combined sequences of ITS2, Rub spacer and Cox3. Phylogenetic trees inferred from Bayesian Inference and Maximum Likelihood methods with other species in the subgenus Sargassum showed that sequences from this study were monophyletic (see Figure S1). All sequences thus were used for subsequent analysis. HKY, F81+I and HKY models were found to be optimum for ITS2, Rub spacer and Cox3 respectively. ML and BI analyses consistently gave same tree topology for all markers with a single clade revealed (see Figures S2-S4).

Population structure

Tajima’s D and Fu’s Fs gave significant negative value for ITS2 and insignificant negative value for Rub spacer and Cox3 (Table S1 in File S1). Neutrality tests indicated there may have been recent population expansion for S. polycystum. Most of the populations showed non-significant pairwise Φ_ST values in all markers indicating shallow population structure (Tables S2-S4 in File S1). Result of hierarchical AMOVA is summarized in Table S5 in File S1. No groupings give significant population structure in ITS2 and Rub spacer. For Cox3, highest percentage of variance (85.38%) was found with the groupings: (1) Southeast Asia; (2) Guam; (3) Pacific Islands. MIGRATE analysis for all markers showed comparable rates of import and export of migrants between subregions (Figure 2 , Tables S6-8 in File S1), indicating a balanced gene flow among these subregions. Mantel test for ITS2, Rub spacer and Cox3 showed low regression coefficient (r = -3.8 × 10^-5; r = -2.7× 10^-5, r = -1.9 × 10^-5, p > 0.05 respectively) that did not fit the isolation by distance model.

Download:

Figure 2. Subregions of the study area defined for MIGRATE analysis.

1: South China Sea and Gulf of Thailand; 2: West Coast of Malay Peninsula; 3: West Java; 4: Celebes Sea and Flores Sea; 5: Guam; 6: Pacific Islands. Directionality of gene flows between selected subregions based on Cox3 is shown with N_e, the mean number of effective migrants per generation. All values are within the range of 440 to 665. Other details of gene flow among subregions based on all three markers are listed in Tables S6-S8 in File S1 respectively.

https://doi.org/10.1371/journal.pone.0077662.g002

Demographic history

Mismatch distribution of all markers fit the sudden expansion model (Figure 3, ITS2: Sum of squared deviation = 0.00260, p>0.05; Rub spacer: Sum of squared deviation = 0.000894, p>0.05; Cox3: Sum of squared deviation = 0.00814, p>0.05). Estimated expansion time calculated from Cox3 was 0.62- 0.41Mya, middle Pleistocene (τ = 0.484, mutation rate = 1.035×10^-9 to 1.555×10^-9 substitutions per site per year, length of Cox3 sequence =379). Bayesian Skyline Plot did not show any recent population expansion reflected by a flat curve (Figure S5).

Download:

Figure 3. Mismatch distribution of Sargassum polycystum using ITS2, Rub spacer and Cox3.

Bar chart and line graph indicate the observed and expected frequency respectively.

https://doi.org/10.1371/journal.pone.0077662.g003

Discussion

Genetic variability of S. polycystum

The overall genetic diversity of S. polycystum is relatively low compared to that of other brown macroalgae. For examples, 51 haplotypes of cox3 were revealed in Sargassum horneri [36], 17 were found in Ishige okamurae [37] and nine were found in Undaria pinnatifida [38] in Northwestern Pacific. In Southeast Asia, study on Gracilaria changii using Cox1 gene revealed six haplotypes in Peninsular Malaysia [7]. The low genetic diversity of S. polycystum is unexpected as the sampling area in the present study is much larger than those in previous studies. Study on S. ilicifolium, having similar sampling sites as the present study, revealed 16, 8, 17 haplotypes for ITS2, Rub spacer and Cox3 respectively (Chan et al., unpublished data). This indicated that the molecular markers used in this study were variable enough to give significant population structures in other closely related Sargassum species and homogeneity found in S. polycystum is not an artifact of the markers used. The difference in genetic diversity may reflect a difference in evolutionary history experienced by S. polycystum and/or the physical environment it is subjected to.

Among the three molecular markers used, Rub spacer showed lowest genetic variability that is comparable to that in other studies on Sargassum [2,36,39]. This showed that Rub spacer has low mutation rate and may not be suitable for phylogeographical studies of Sargassum. On the other hand, genetic variabilities of ITS2 and Cox3 found in the present study are similar (H = 0.28 ± 0.04, π = 0.030 ± 0.032; H = 0.39 ± 0.03, π =0.079 ± 0.076).

Population Structure

Overall, the population structure of S. polycystum is relatively homogeneous with few haplotypes derived from the main centre haplotype (Figure 1). In contrast, for animals, high genetic diversity is observed in giant clam [40,41], stomatopods [3], clownfish [4] and coral [42] within the same region due to past isolation in refugia and recent oceanographic currents which limited gene flow. For example, exposed Sunda Shelf during Last Glacial Maximum caused isolation of clownfish populations between Sulawesi and Indo-Malay Archipelago [4] and downward directional flow of Indonesian Throughflow maintained the genetic connectivity of stomatopod populations along the west coast of Sulawesi and East Java [3]. Phylogeographical analysis of three Sargassum spp. in northwestern Pacific [2,36,43] also showed Sargassum populations could both be structured by Pleistocene isolation and postglacial recolonization associated with current flows (e.g. Kuroshio Current). These factors, however, seem to have little effect on S. polycystum as no significant population structure can be observed within the Southeast Asian region. No population isolation is observed in Sulawesi and other identified refugia (e.g. South China Sea and Sulu Sea) and the Indonesian Throughflow did not appear to impose strong restricted dispersal between Java Sea and Flores Sea. This lack of population structure is further supported by both AMOVA and MIGRATE results. AMOVA did not identify clear groupings of S. polycyctum populations within Southeast Asian region so AMOVA results could not be used as a reference to estimate the gene flows between groups. Nonetheless, subregions based on putative marine basins or refugia were used in MIGRATE analysis and the results also revealed relatively balanced gene flows between subregions.

Different from the marine animals studied, S. polycystum is capable of both sexual and asexual reproductions. It was reported that S. polycystum can reproduce asexually with new main axes growing out from its ramifying holdfast (stolons) [44]. Indeed, our own experiments also indicated that fragments of this ramifying holdfast (stolons) can develop new shoots and new axes, which subsequently could reattach to the substrata (Ang et al, unpublished data). This feature is unique and uncommon in the genus Sargassum [11]. Interestingly, S. hemiphyllum var. chinensis, which also has a ramifying holdfast, also showed a low genetic diversity within its range of distribution in southern China [2]. Asexual propagation of S. polycystum could thus be suggested to be involved in maintaining the genetic homogeneity of the species. Similar process has also been suggested for the brown alga Fucus vesiculosus [45]. However, the ratio of sexual and asexual reproductions in the populations of S. polycystum remains unknown. Calculation of incidence of asexual vs. sexual dispersal is difficult as comparison with known asexually vs. sexually propagated populations of other Sargassum species within the same distribution range is needed. To date, however, only S. polycystum is putatively known to propagate asexually. Further investigation using multilocus genetic markers would also be needed to further confirm the roles of asexual propagation.

Asexual propagation alone would not be sufficient to explain the lack of genetic structure in S. polycystum populations across such a big geographical region covered in the present study. Homogeneity of the populations could also indicate high dispersal ability of the species as no isolation by distance was observed in S. polycystum nor was there distinct directional gene flow revealed by MIGRATE analysis. For animal species, population structures are determined by larval duration and current flow. In contrast, Sargassum is capable of dispersing for long distance in the form of drifting fronds with germlings [10] and rapid recolonzation of a site could be facilitated. Homogeneity was also observed in other rafting seaweeds. For example, samples from subantarctic islands of Southern bull kelp Durvillaea antarctica showed homogeneity dominated by one single haplotype over a distance of 10,000 km [46]. High genetic connectivity was also observed in the giant kelp Macrocystis pyrifera along the subantarctic islands and southern-central Chile [47] with high degree of shared haplotypes. All these can be explained by high dispersal ability of kelp rafts [47]. In northwestern Pacific, drifting of Sargassum horneri with coastal currents was found to contribute genetic connectivity in East China Sea and allowed secondary contacts between the Chinese and Japanese populations [36]. In contrast, study in central Philippines revealed 21 cryptic species of the red alga Portieria, suggesting that speciation may occur within < 100 km due to the absence of propagules, which in turn limits the dispersal ability of the populations [48]. High endemism was thus observed with no or limited overlaps of species even though islands in central Philippines are highly connected by ocean currents. The dispersal mode of Sargassum, having drifting fronds with germlings, suggests that it has a much higher dispersal potential than Portieria. The long persistence of drifting fronds of Sargassum allowed eventual mixing of isolated populations (e.g. by Indonesian Throughflow) and explains the homogeneity across long distance, i.e. from Pacific Islands to Southeast Asia.

Nonetheless, despite an overall homogeneous pattern of S. polycystum populations observed, small discrepancies could still be found in the Philippine, Indonesian and Thai populations based on ITS2 and Guam and eastern Indonesian populations based on Cox3. These discrepancies were not consistent between the patterns generated from ITS2 and Rub spacer and those from Cox3 (Figure 1). In either situation, some isolated cases of mutation had occurred in these populations but their dispersal was restricted. For example, there appears to have a single base pair mutation in the Guam population, as suggested by Cox3. By and large, the Guam population could have been restricted in its dispersal by the directional flow of the North Equatorial Current (NEC in Figure 1). However, some individuals could have gotten into the Indonesia Throughflow and eventually reached Java (TF and PK). The frequency of the Guam haplotype in Javan populations is very low and the absence of this haplotype in any sites between Guam and Java may simply be a sampling problem.

Demographic history

Star-shaped haplotype network (Figure 1) and mismatch distribution of Cox3 (Figure 3) indicated demographic expansion of S. polycystum populations in mid-Pleistocene (0.62- 0.41Mya) that remained constant recently, as inferred from Bayesian Skyline Plot (Figure S5). The homogeneous pattern resembles null model I of panmixia proposed by Maggs et al. [49]. This model suggested no significant population differentiation but with only shallow mixture of derived haplotypes. Panmixia was also found in the European eel, Anguilla anguilla L. with low genetic differentiation between populations (FST = 0.0014; p<0.01) [50]. It can be suggested that population of S. polycystum retracted during mid-Pleistocene and persisted in a single refugium which would have included the Philippines and western Pacific islands. After the Last Glacial Maximum, recolonization of the rest of the slowly flooded Sunda Shelf was likely achieved mainly through asexual propagation and rapid dispersal in the form of drifting fronds with germlings. Mutation of the species was slow resulting in homogeneity of the populations across the whole Southeast Asian region.

Conclusion

The absence of significant population structure in S. polycystum over a wide area in Southeast Asia and western Pacific suggests high genetic connectivity within this region. Asexual propagation and high dispersal ability are proposed to be the mechanisms that contributed to this observed phenomenon. There appears to be no significant isolation of S. polycystum populations during the lowering of sea levels and demographic expansion may have occurred more recently, after the Last Glacial Maximum. The species could have persisted in a single refugium before recolonization occurred. To further reveal the evolutionary history of this species in Southeast Asia, physiological experiment on the ratio of sexual and asexual reproduction would be needed. Multilocus molecular markers with higher variability (e.g. microsatellites) are suggested to reveal the hidden diversity among the populations. More sampling sites, especially in Indian Ocean, may also help to elucidate the extent of population homogeneity in this species.

Supporting Information

File S1.

Combined file containing all supporting tables. Table S1: Sample localities, diversity indices and neutrality tests of Sargassum polycystum using ITS2, Rub spacer and Cox3 (from top to bottom) including sequences from Genbank. Table S2: Pairwise Φ_ST among populations of Sargassum polycystum based on ITS2. Table S3: Pairwise Φ_ST among populations of Sargassum polycystum based on Rub spacer. Table S4: Pairwise Φ_ST among populations of Sargassum polycystum based on Cox3. Table S5: Hierarchical Analysis of Molecular Variance (AMOVA) based on the three markers ITS2, Rub spacer and Cox3 from Sargassum polycystum. Table S6: Migration estimates of gene flow between subregions 1 to 6 and 95% confidence values as estimated by MIGRATE using ITS2. Table S7: Migration estimates of gene flow between subregions 1 to 6 and 95% confidence values as estimated by MIGRATE using Rub spacer. Table S8: Migration estimates of gene flow between subregions 1 to 6 and 95% confidence values as estimated by MIGRATE using Cox3.

https://doi.org/10.1371/journal.pone.0077662.s001

(DOCX)

Figure S1.

Phylogenetic tree of representative sequences from each haplotype of combined ITS2, Rub spacer and Cox3 sequences from Sargassum polycystum with other Sargassum spp. in subgenus Sargassum.

Voucher numbers of sequences from Genbank are presented after species name. Representative sequences from each haplotype of S. polycystum are labeled in black bar. Turbinaria ornata is used as outgroup. Posterior probabilities of Bayesian Inference and boostrap value of Maximum Likelihood are shown. Tables of sample localities, diversity indices, neutrality tests; Pairwise <PHI>ST; Hierarchical Analysis of Molecular Variance (AMOVA) and Migration estimates of gene flow between subregions of three genes of S. polycystum.

https://doi.org/10.1371/journal.pone.0077662.s002

(TIF)

Figure S2.

Phylogenetic tree inferred from ITS2 using Sargassum ilicifolium as outgroup. Posterior probabilities of Bayesian Inference and boostrap values of Maximum likelihood are shown.

https://doi.org/10.1371/journal.pone.0077662.s003

(TIF)

Figure S3.

Phylogenetic tree inferred from Rub spacer using Sargassum ilicifolium as outgroup. Posterior probabilities of Bayesian Inference and boostrap values of Maximum likelihood are shown.

https://doi.org/10.1371/journal.pone.0077662.s004

(TIF)

Figure S4.

Phylogenetic tree inferred from Cox3 using Sargassum ilicifolium as outgroup. Posterior probabilities of Bayesian Inference and boostrap values of Maximum likelihood are shown.

https://doi.org/10.1371/journal.pone.0077662.s005

(TIF)

Figure S5.

Bayesian Skyline Plot of Cox3 in effective population size with function of time (year before present). 95% confidence interval is shown in blue.

https://doi.org/10.1371/journal.pone.0077662.s006

(TIF)

Acknowledgments

We thank Phang Siew Moi (Malaysia), Floredel Dangan, Edna Ganzon-Fortes, Lea Jimenez, Danilo Largo, Filemon Romero (Philippines), Jaruwan Mayakun, Achana Prathep (Thailand), Dang Diem Hong, Dam Duc Tien (Vietnam), for their hospitality during our collecting trips or for help with the collection of samples. Yeung Chung Wing assisted in collection especially in sites in Indonesia. We also thank Jeanine Olsen (Netherland) for a fruitful discussion on asexual propagation in Sargassum. Comments from two anonymous reviewers were also very helpful.

Author Contributions

Conceived and designed the experiments: PA SWC. Performed the experiments: SWC. Analyzed the data: SWC CCC PA AC. Contributed reagents/materials/analysis tools: SWC PA CCC AC GG CT. Wrote the manuscript: SWC PA CCC. Served as contact person / collaborator in Indonesia so no specific permission for collecting in Indonesia was required: GG.

References

1. Voris HK (2000) Maps of Pleistocene sea levels in Southeast Asia: shorelines, river systems and time durations. J Biogeogr 27: 1153–1167. doi:https://doi.org/10.1046/j.1365-2699.2000.00489.x.
- View Article
- Google Scholar
2. Cheang CC, Chu KH, Ang PO (2010) Phylogeography of the marine macroalga Sargassum hemiphyllum (Phaeophyceae, Heterokontophyta) in northwestern Pacific. Mol Ecol 19: 2933–2948. doi:https://doi.org/10.1111/j.1365-294X.2010.04685.x. PubMed: 20584136.
- View Article
- Google Scholar
3. Barber PH, Erdmann MV, Palumbi SR (2006) Comparative phylogeography of three codistributed stomatopods: origins and timing of regional lineage diversification in the Coral Triangle. Evolution 60: 1825–1839. doi:https://doi.org/10.1111/j.0014-3820.2006.tb00526.x. PubMed: 17089967.
- View Article
- Google Scholar
4. Timm J, Kochzius M (2008) Geological history and oceanography of the Indo-Malay Archipelago shape the genetic population structure in the false clown anemonefish (Amphiprion ocellaris). Mol Ecol 17: 3999–4014. doi:https://doi.org/10.1111/j.1365-294X.2008.03881.x. PubMed: 19238702.
- View Article
- Google Scholar
5. Hoarau G, Coyer JA, Veldsink JH, Stam WT, Olsen JL (2007) Glacial refugia and recolonization pathways in the brown seaweed Fucus serratus. Mol Ecol 16: 3606–3616. doi:https://doi.org/10.1111/j.1365-294X.2007.03408.x. PubMed: 17845434.
- View Article
- Google Scholar
6. Andreakis N, Kooistra WHCF, Procaccini G (2009) High genetic diversity and connectivity in the polyploid invasive seaweed Asparagopsis taxiformis (Bonnemaisoniales) in the Mediterranean, explored with microsatellite alleles and multilocus genotypes. Mol Ecol 18: 212–226. doi:https://doi.org/10.1111/j.1365-294X.2008.04022.x. PubMed: 19192176.
- View Article
- Google Scholar
7. Yow YY, Lim PE, Phang SM (2010) Genetic diversity of Gracilaria changii (Gracilariaceae, Rhodophyta) from west coast, Peninsular Malaysia based on mitochondrial cox1 gene analysis. J Appl Phycol 23: 219–226.
- View Article
- Google Scholar
8. Guiry MD, Guiry GM (2013) AlgaeBase. World-wide electronic publication, National University of Ireland, Galway. Available: http://www.algaebase.org. Accessed 10 June 2013.
9. Kendrick GA, Walker DI (1995) Dispersal of propagules of Sargassum spp. (Sargassaceae: Phaeophyta): Observations of local patterns of dispersal and consequences for recruitment and population structure. J Exp Mar Biol Ecol 192: 273–288. doi:https://doi.org/10.1016/0022-0981(95)00076-4.
- View Article
- Google Scholar
10. Komatsu T, Tatsukawa K, Filippi JB, Sagawa T, Matsunaga D et al. (2007) Distribution of drifting seaweeds in eastern East China Sea. J Mar Syst 67: 245–252. doi:https://doi.org/10.1016/j.jmarsys.2006.05.018.
- View Article
- Google Scholar
11. Tseng CK, Baoren Lu (1992) Studies on Chinese species of zygocarpic Sargassum. In IA Abbott. Taxonomy of Economic Seaweeds with reference to some Pacific and Caribbean species. Vol.2. California Sea Grant College Program. La Jolla: University of California. pp. 23-54.
12. Chiang YM, Yoshida T, Ajisaka T, Trono G Jr, Tseng CK et al. (1992) Distribution and variation in Sargassum polycystum C.A. Agardh (Fucales, Phaeophyta). In IA Abbott. Taxonomy of Economic Seaweeds with reference to some Pacific and Western Atlantic species. Vol.3. California Sea Grant College Program. La Jolla: University of California. pp. 35-42.
13. Gordon A, Fine R (1996) Pathways of water between the Pacific and Indian oceans in the Indonesian seas. Nature 379: 146–149. doi:https://doi.org/10.1038/379146a0.
- View Article
- Google Scholar
14. Wyrtki K (1961) Scientific results of marine investigations of the South China Sea and the Gulf of Thailand 1959-1961. NAGA Report 2. : 1–195.
15. Mattio L, Payri CE, Verlaque M (2009) Taxonomic revision and geographic distribution of the subgenus Sargassum (Fucales, Phaeophyceae) in the western and central Pacific islands based on morphological and molecular analyses. J Phycol 45: 1213–1227. doi:https://doi.org/10.1111/j.1529-8817.2009.00737.x.
- View Article
- Google Scholar
16. Tamura K, Peterson D, Peterson N, Stecher G, Nei M et al. (2011) MEGA5: Molecular Evolutionary Genetics Analysis using Maximum Likelihood, Evolutionary Distance, and Maximum Parsimony Methods. Mol Biol Evol 28: 2731–2739. doi:https://doi.org/10.1093/molbev/msr121. PubMed: 21546353.
- View Article
- Google Scholar
17. Ronquist F, Teslenko M, van der Mark P, Ayres DL, Darling A et al. (2012) MrBayes 3.2: efficient Bayesian phylogenetic inference and model choice across a large model space. Syst Biol 61: 539–542. doi:https://doi.org/10.1093/sysbio/sys029. PubMed: 22357727.
- View Article
- Google Scholar
18. Guindon S, Dufayard J-F, Lefort V, Anisimova M, Hordijk W et al. (2010) New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML. Syst Biol 3.0: 59: 307–321.
- View Article
- Google Scholar
19. Darriba D, Taboada GL, Doallo R, Posada D (2012) jModelTest 2: more models, new heuristics and parallel computing. Nat Methods 9: 772. doi:https://doi.org/10.1038/nmeth.2111.
- View Article
- Google Scholar
20. Stiger V, Horiguchi T (2003) Phylogenetic relationships within the genus Sargassum ( Fucales , Phaeophyceae ), inferred from ITS-2 nrDNA , with an emphasis on the taxonomic subdivision of the genus. Phycol Res 51: 1–10. doi:https://doi.org/10.1111/j.1440-1835.2003.tb00164.x.
- View Article
- Google Scholar
21. Phillips N, Smith C, Morden C (2005) Testing systematic concepts of Sargassum (Fucales, Phaeophyceae) using portions of the rbcLS operon. Phycol Res 53: 1–10. doi:https://doi.org/10.1111/j.1440-1835.2005.tb00353.x.
- View Article
- Google Scholar
22. Clement M, Posada D, Crandall KA (2000) TCS: a computer program to estimate gene genealogies. Mol Ecol 9: 1657–1659. doi:https://doi.org/10.1046/j.1365-294x.2000.01020.x. PubMed: 11050560.
- View Article
- Google Scholar
23. Excoffier L, Lischer HEL (2010) Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows. Mol Ecol Resour 10: 564–567. doi:https://doi.org/10.1111/j.1755-0998.2010.02847.x. PubMed: 21565059.
- View Article
- Google Scholar
24. Nei M (1987) Molecular Evolutionary Genetics. New York: Columbia University Press.
25. Nei M, Li J (1989) Variances of the average numbers of nucleotide substitutions within and between populations. Mol Biol Evol 6: 290–300. PubMed: 2576093.
- View Article
- Google Scholar
26. Tajima F (1989) Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics 123: 585–595. PubMed: 2513255.
- View Article
- Google Scholar
27. Fu YX (1997) Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection. Genetics 147: 915–925. PubMed: 9335623.
- View Article
- Google Scholar
28. Beerli P (2006) Comparison of Bayesian and maximum likelihood inference of population genetic parameters. Bioinformatics 22 : 341–345. doi:https://doi.org/10.1093/bioinformatics/bti803. PubMed: 16317072.
- View Article
- Google Scholar
29. Beerli P (2009) How to use migrate or why are markov chain monte carlo programs difficult to use? In: G. BertorelleMW BrufordHC HauffeA. RizzoliC. Vernesi. Population Genetics for Animal Conservation, volume 17 of Conservation Biology. Cambridge, UK: Cambridge University Press. pp, 42–79.
30. Rousset F (1997) Genetic differentiation and estimation of gene flow from F-statistics under isolation by distance. Genetics 145: 1219–1228. PubMed: 9093870.
- View Article
- Google Scholar
31. Rogers AR, Harpending H (1992) Population growth makes waves in the distribution of pairwise genetic differences. Mol Biol Evol 9: 552–569. PubMed: 1316531.
- View Article
- Google Scholar
32. Drummond AJ, Rambaut A, Shapiro B, Pybus OG (2005) Bayesian coalescent inference of past population dynamics from molecular sequences. Mol Biol Evol 22: 1185–1192. doi:https://doi.org/10.1093/molbev/msi103. PubMed: 15703244.
- View Article
- Google Scholar
33. Drummond AJ, Rambaut A (2007) BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol Biol 7: 214. doi:https://doi.org/10.1186/1471-2148-7-214. PubMed: 17996036.
- View Article
- Google Scholar
34. Trono GC, Lluisma AO (1990) Seasonality of standing crop of a Sargassum (Fucales, Phaeophyta) bed in Bolinao, Pangasinan, Philippines. Hydrobiologia 204/205: 331–338. doi:https://doi.org/10.1007/BF00040253.
- View Article
- Google Scholar
35. Rao A, Rao M (2002) Seasonal growth pattern in Sargassum polycystum C Agardh (Phaeophyta, Fucales) occurring at Visakhapatnam, east coast of India. Indian J Mar Sci 31: 26–32.
- View Article
- Google Scholar
36. Hu ZM, Uwai S, Yu SH, Komatsu T, Ajisaka T et al. (2011) Phylogeographic heterogeneity of the brown macroalga Sargassum horneri (Fucaceae) in the northwestern Pacific in relation to late Pleistocene glaciation and tectonic configurations. Mol Ecol 20: 3894–3909. doi:https://doi.org/10.1111/j.1365-294X.2011.05220.x. PubMed: 21851438.
- View Article
- Google Scholar
37. Lee KM, Yang EC, Coyer JA, Zuccarello GC, Wang WL et al. (2012) Phylogeography of the seaweed Ishige okamurae (Phaeophyceae): evidence for glacial refugia in the northwest Pacific region. Mar Biol 159: 1021–1028. doi:https://doi.org/10.1007/s00227-012-1882-0.
- View Article
- Google Scholar
38. Uwai S, Yotsukura N, Serisawa Y, Muraoka D, Hiraoka M et al. (2006 ) Intraspecific Genetic Diversity of Undaria pinnatifida in Japan, based on the Mitochondrial cox3 Gene and the ITS1 of nrDNA. Hydrobiologia 553: 345–356. doi:https://doi.org/10.1007/s10750-005-0883-0.
- View Article
- Google Scholar
39. Cheang CC, Chu KH, Fujita D, Yoshida G, Hiraoka M et al. (2010) Low genetic variability of Sargassum muticum (Phaeophyceae) revealed by a global analysis of native and introduced populations. J Phycol 46: 1063–1074. doi:https://doi.org/10.1111/j.1529-8817.2010.00901.x.
- View Article
- Google Scholar
40. Deboer TS, Subia MD, Erdmann MV, Kovitvongsa K, Barber PH (2008) Phylogeography and limited genetic connectivity in the endangered boring giant clam across the Coral Triangle. Conserv Biol 22: 1255–1266. doi:https://doi.org/10.1111/j.1523-1739.2008.00983.x. PubMed: 18637905.
- View Article
- Google Scholar
41. Kochzius M, Nuryanto A (2008) Strong genetic population structure in the boring giant clam, Tridacna crocea, across the Indo-Malay Archipelago: implications related to evolutionary processes and connectivity. Mol Ecol 17: 3775–3787. doi:https://doi.org/10.1111/j.1365-294X.2008.03803.x. PubMed: 18662232.
- View Article
- Google Scholar
42. Knittweis L, Kraemer WE, Timm J, Kochzius M (2008) Genetic structure of Heliofungia actiniformis (Scleractinia: Fungiidae) populations in the Indo-Malay Archipelago: implications for live coral trade management efforts. Conserv Genet 10: 241–249.
- View Article
- Google Scholar
43. Hu ZM, Zhang J, Lopez-Bautista J, Duan DL (2013) Asymmetric genetic exchange in the brown seaweed Sargassum fusiforme (Phaeophyceae) driven by oceanic currents. Mar Biol 160 : 1407–1414. doi:https://doi.org/10.1007/s00227-013-2192-x.
- View Article
- Google Scholar
44. Yeong MLB, Wong CL (2012) Seasonal growth rate of Sargassum species at Teluk Kemang, Port Dickson, Malaysia. J Appl Phycol 25: 805–814.
- View Article
- Google Scholar
45. Tatarenkov A, Bergström L, Jönsson RB, Serrão EA, Kautsky L et al. (2005) Intriguing asexual life in marginal populations of the brown seaweed Fucus vesiculosus. Mol Ecol 14: 647–651. doi:https://doi.org/10.1111/j.1365-294X.2005.02425.x. PubMed: 15660953.
- View Article
- Google Scholar
46. Fraser CI, Spencer HG, Waters JM (2009) Glacial oceanographic contrasts explain phylogeography of Australian bull kelp. Mol Ecol 18: 2287–2296. doi:https://doi.org/10.1111/j.1365-294X.2009.04201.x. PubMed: 19389161.
- View Article
- Google Scholar
47. Macaya EC, Zuccarello G (2010) Genetic structure of the giant kelp Macrocystis pyrifera along the southeastern Pacific. Mar Ecol-Prog Ser 420: 103–112. doi:https://doi.org/10.3354/meps08893.
- View Article
- Google Scholar
48. Payo DA, Leliaert F, Verbruggen H, D’hondt S, Calumpong HP et al. (2013) Extensive cryptic species diversity and fine-scale endemism in the marine red alga Portieria in the Philippines. Proc R Soc Lond B Biol Sci 280: 1753–1760. PubMed: 23269854.
- View Article
- Google Scholar
49. Maggs CA, Castilho R, Foltz D, Henzler C, Jolly MT et al. (2008) Evaluating signatures of glacial refugia for North Atlantic benthic marine taxa. Ecology 89: S108–S122. doi:https://doi.org/10.1890/08-0257.1. PubMed: 19097488.
- View Article
- Google Scholar
50. Dannewitz J (2005) Panmixia in the European eel: a matter of time. Proc R Soc Lond B Biol Sci 272: 1129–1137. doi:https://doi.org/10.1098/rspb.2005.3064. PubMed: 16024374.
- View Article
- Google Scholar

[ref1] 1. Voris HK (2000) Maps of Pleistocene sea levels in Southeast Asia: shorelines, river systems and time durations. J Biogeogr 27: 1153–1167. doi:https://doi.org/10.1046/j.1365-2699.2000.00489.x.
View Article
Google Scholar

[2] View Article

[3] Google Scholar

[ref2] 2. Cheang CC, Chu KH, Ang PO (2010) Phylogeography of the marine macroalga Sargassum hemiphyllum (Phaeophyceae, Heterokontophyta) in northwestern Pacific. Mol Ecol 19: 2933–2948. doi:https://doi.org/10.1111/j.1365-294X.2010.04685.x. PubMed: 20584136.
View Article
Google Scholar

[5] View Article

[6] Google Scholar

[ref3] 3. Barber PH, Erdmann MV, Palumbi SR (2006) Comparative phylogeography of three codistributed stomatopods: origins and timing of regional lineage diversification in the Coral Triangle. Evolution 60: 1825–1839. doi:https://doi.org/10.1111/j.0014-3820.2006.tb00526.x. PubMed: 17089967.
View Article
Google Scholar

[8] View Article

[9] Google Scholar

[ref4] 4. Timm J, Kochzius M (2008) Geological history and oceanography of the Indo-Malay Archipelago shape the genetic population structure in the false clown anemonefish (Amphiprion ocellaris). Mol Ecol 17: 3999–4014. doi:https://doi.org/10.1111/j.1365-294X.2008.03881.x. PubMed: 19238702.
View Article
Google Scholar

[11] View Article

[12] Google Scholar

[ref5] 5. Hoarau G, Coyer JA, Veldsink JH, Stam WT, Olsen JL (2007) Glacial refugia and recolonization pathways in the brown seaweed Fucus serratus. Mol Ecol 16: 3606–3616. doi:https://doi.org/10.1111/j.1365-294X.2007.03408.x. PubMed: 17845434.
View Article
Google Scholar

[14] View Article

[15] Google Scholar

[ref6] 6. Andreakis N, Kooistra WHCF, Procaccini G (2009) High genetic diversity and connectivity in the polyploid invasive seaweed Asparagopsis taxiformis (Bonnemaisoniales) in the Mediterranean, explored with microsatellite alleles and multilocus genotypes. Mol Ecol 18: 212–226. doi:https://doi.org/10.1111/j.1365-294X.2008.04022.x. PubMed: 19192176.
View Article
Google Scholar

[17] View Article

[18] Google Scholar

[ref7] 7. Yow YY, Lim PE, Phang SM (2010) Genetic diversity of Gracilaria changii (Gracilariaceae, Rhodophyta) from west coast, Peninsular Malaysia based on mitochondrial cox1 gene analysis. J Appl Phycol 23: 219–226.
View Article
Google Scholar

[20] View Article

[21] Google Scholar

[ref8] 8. Guiry MD, Guiry GM (2013) AlgaeBase. World-wide electronic publication, National University of Ireland, Galway. Available: http://www.algaebase.org. Accessed 10 June 2013.

[ref9] 9. Kendrick GA, Walker DI (1995) Dispersal of propagules of Sargassum spp. (Sargassaceae: Phaeophyta): Observations of local patterns of dispersal and consequences for recruitment and population structure. J Exp Mar Biol Ecol 192: 273–288. doi:https://doi.org/10.1016/0022-0981(95)00076-4.
View Article
Google Scholar

[24] View Article

[25] Google Scholar

[ref10] 10. Komatsu T, Tatsukawa K, Filippi JB, Sagawa T, Matsunaga D et al. (2007) Distribution of drifting seaweeds in eastern East China Sea. J Mar Syst 67: 245–252. doi:https://doi.org/10.1016/j.jmarsys.2006.05.018.
View Article
Google Scholar

[27] View Article

[28] Google Scholar

[ref11] 11. Tseng CK, Baoren Lu (1992) Studies on Chinese species of zygocarpic Sargassum. In IA Abbott. Taxonomy of Economic Seaweeds with reference to some Pacific and Caribbean species. Vol.2. California Sea Grant College Program. La Jolla: University of California. pp. 23-54.

[ref12] 12. Chiang YM, Yoshida T, Ajisaka T, Trono G Jr, Tseng CK et al. (1992) Distribution and variation in Sargassum polycystum C.A. Agardh (Fucales, Phaeophyta). In IA Abbott. Taxonomy of Economic Seaweeds with reference to some Pacific and Western Atlantic species. Vol.3. California Sea Grant College Program. La Jolla: University of California. pp. 35-42.

[ref13] 13. Gordon A, Fine R (1996) Pathways of water between the Pacific and Indian oceans in the Indonesian seas. Nature 379: 146–149. doi:https://doi.org/10.1038/379146a0.
View Article
Google Scholar

[32] View Article

[33] Google Scholar

[ref14] 14. Wyrtki K (1961) Scientific results of marine investigations of the South China Sea and the Gulf of Thailand 1959-1961. NAGA Report 2. : 1–195.

[ref15] 15. Mattio L, Payri CE, Verlaque M (2009) Taxonomic revision and geographic distribution of the subgenus Sargassum (Fucales, Phaeophyceae) in the western and central Pacific islands based on morphological and molecular analyses. J Phycol 45: 1213–1227. doi:https://doi.org/10.1111/j.1529-8817.2009.00737.x.
View Article
Google Scholar

[36] View Article

[37] Google Scholar

[ref16] 16. Tamura K, Peterson D, Peterson N, Stecher G, Nei M et al. (2011) MEGA5: Molecular Evolutionary Genetics Analysis using Maximum Likelihood, Evolutionary Distance, and Maximum Parsimony Methods. Mol Biol Evol 28: 2731–2739. doi:https://doi.org/10.1093/molbev/msr121. PubMed: 21546353.
View Article
Google Scholar

[39] View Article

[40] Google Scholar

[ref17] 17. Ronquist F, Teslenko M, van der Mark P, Ayres DL, Darling A et al. (2012) MrBayes 3.2: efficient Bayesian phylogenetic inference and model choice across a large model space. Syst Biol 61: 539–542. doi:https://doi.org/10.1093/sysbio/sys029. PubMed: 22357727.
View Article
Google Scholar

[42] View Article

[43] Google Scholar

[ref18] 18. Guindon S, Dufayard J-F, Lefort V, Anisimova M, Hordijk W et al. (2010) New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML. Syst Biol 3.0: 59: 307–321.
View Article
Google Scholar

[45] View Article

[46] Google Scholar

[ref19] 19. Darriba D, Taboada GL, Doallo R, Posada D (2012) jModelTest 2: more models, new heuristics and parallel computing. Nat Methods 9: 772. doi:https://doi.org/10.1038/nmeth.2111.
View Article
Google Scholar

[48] View Article

[49] Google Scholar

[ref20] 20. Stiger V, Horiguchi T (2003) Phylogenetic relationships within the genus Sargassum ( Fucales , Phaeophyceae ), inferred from ITS-2 nrDNA , with an emphasis on the taxonomic subdivision of the genus. Phycol Res 51: 1–10. doi:https://doi.org/10.1111/j.1440-1835.2003.tb00164.x.
View Article
Google Scholar

[51] View Article

[52] Google Scholar

[ref21] 21. Phillips N, Smith C, Morden C (2005) Testing systematic concepts of Sargassum (Fucales, Phaeophyceae) using portions of the rbcLS operon. Phycol Res 53: 1–10. doi:https://doi.org/10.1111/j.1440-1835.2005.tb00353.x.
View Article
Google Scholar

[54] View Article

[55] Google Scholar

[ref22] 22. Clement M, Posada D, Crandall KA (2000) TCS: a computer program to estimate gene genealogies. Mol Ecol 9: 1657–1659. doi:https://doi.org/10.1046/j.1365-294x.2000.01020.x. PubMed: 11050560.
View Article
Google Scholar

[57] View Article

[58] Google Scholar

[ref23] 23. Excoffier L, Lischer HEL (2010) Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows. Mol Ecol Resour 10: 564–567. doi:https://doi.org/10.1111/j.1755-0998.2010.02847.x. PubMed: 21565059.
View Article
Google Scholar

[60] View Article

[61] Google Scholar

[ref24] 24. Nei M (1987) Molecular Evolutionary Genetics. New York: Columbia University Press.

[ref25] 25. Nei M, Li J (1989) Variances of the average numbers of nucleotide substitutions within and between populations. Mol Biol Evol 6: 290–300. PubMed: 2576093.
View Article
Google Scholar

[64] View Article

[65] Google Scholar

[ref26] 26. Tajima F (1989) Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics 123: 585–595. PubMed: 2513255.
View Article
Google Scholar

[67] View Article

[68] Google Scholar

[ref27] 27. Fu YX (1997) Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection. Genetics 147: 915–925. PubMed: 9335623.
View Article
Google Scholar

[70] View Article

[71] Google Scholar

[ref28] 28. Beerli P (2006) Comparison of Bayesian and maximum likelihood inference of population genetic parameters. Bioinformatics 22 : 341–345. doi:https://doi.org/10.1093/bioinformatics/bti803. PubMed: 16317072.
View Article
Google Scholar

[73] View Article

[74] Google Scholar

[ref29] 29. Beerli P (2009) How to use migrate or why are markov chain monte carlo programs difficult to use? In: G. BertorelleMW BrufordHC HauffeA. RizzoliC. Vernesi. Population Genetics for Animal Conservation, volume 17 of Conservation Biology. Cambridge, UK: Cambridge University Press. pp, 42–79.

[ref30] 30. Rousset F (1997) Genetic differentiation and estimation of gene flow from F-statistics under isolation by distance. Genetics 145: 1219–1228. PubMed: 9093870.
View Article
Google Scholar

[77] View Article

[78] Google Scholar

[ref31] 31. Rogers AR, Harpending H (1992) Population growth makes waves in the distribution of pairwise genetic differences. Mol Biol Evol 9: 552–569. PubMed: 1316531.
View Article
Google Scholar

[80] View Article

[81] Google Scholar

[ref32] 32. Drummond AJ, Rambaut A, Shapiro B, Pybus OG (2005) Bayesian coalescent inference of past population dynamics from molecular sequences. Mol Biol Evol 22: 1185–1192. doi:https://doi.org/10.1093/molbev/msi103. PubMed: 15703244.
View Article
Google Scholar

[83] View Article

[84] Google Scholar

[ref33] 33. Drummond AJ, Rambaut A (2007) BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol Biol 7: 214. doi:https://doi.org/10.1186/1471-2148-7-214. PubMed: 17996036.
View Article
Google Scholar

[86] View Article

[87] Google Scholar

[ref34] 34. Trono GC, Lluisma AO (1990) Seasonality of standing crop of a Sargassum (Fucales, Phaeophyta) bed in Bolinao, Pangasinan, Philippines. Hydrobiologia 204/205: 331–338. doi:https://doi.org/10.1007/BF00040253.
View Article
Google Scholar

[89] View Article

[90] Google Scholar

[ref35] 35. Rao A, Rao M (2002) Seasonal growth pattern in Sargassum polycystum C Agardh (Phaeophyta, Fucales) occurring at Visakhapatnam, east coast of India. Indian J Mar Sci 31: 26–32.
View Article
Google Scholar

[92] View Article

[93] Google Scholar

[ref36] 36. Hu ZM, Uwai S, Yu SH, Komatsu T, Ajisaka T et al. (2011) Phylogeographic heterogeneity of the brown macroalga Sargassum horneri (Fucaceae) in the northwestern Pacific in relation to late Pleistocene glaciation and tectonic configurations. Mol Ecol 20: 3894–3909. doi:https://doi.org/10.1111/j.1365-294X.2011.05220.x. PubMed: 21851438.
View Article
Google Scholar

[95] View Article

[96] Google Scholar

[ref37] 37. Lee KM, Yang EC, Coyer JA, Zuccarello GC, Wang WL et al. (2012) Phylogeography of the seaweed Ishige okamurae (Phaeophyceae): evidence for glacial refugia in the northwest Pacific region. Mar Biol 159: 1021–1028. doi:https://doi.org/10.1007/s00227-012-1882-0.
View Article
Google Scholar

[98] View Article

[99] Google Scholar

[ref38] 38. Uwai S, Yotsukura N, Serisawa Y, Muraoka D, Hiraoka M et al. (2006 ) Intraspecific Genetic Diversity of Undaria pinnatifida in Japan, based on the Mitochondrial cox3 Gene and the ITS1 of nrDNA. Hydrobiologia 553: 345–356. doi:https://doi.org/10.1007/s10750-005-0883-0.
View Article
Google Scholar

[101] View Article

[102] Google Scholar

[ref39] 39. Cheang CC, Chu KH, Fujita D, Yoshida G, Hiraoka M et al. (2010) Low genetic variability of Sargassum muticum (Phaeophyceae) revealed by a global analysis of native and introduced populations. J Phycol 46: 1063–1074. doi:https://doi.org/10.1111/j.1529-8817.2010.00901.x.
View Article
Google Scholar

[104] View Article

[105] Google Scholar

[ref40] 40. Deboer TS, Subia MD, Erdmann MV, Kovitvongsa K, Barber PH (2008) Phylogeography and limited genetic connectivity in the endangered boring giant clam across the Coral Triangle. Conserv Biol 22: 1255–1266. doi:https://doi.org/10.1111/j.1523-1739.2008.00983.x. PubMed: 18637905.
View Article
Google Scholar

[107] View Article

[108] Google Scholar

[ref41] 41. Kochzius M, Nuryanto A (2008) Strong genetic population structure in the boring giant clam, Tridacna crocea, across the Indo-Malay Archipelago: implications related to evolutionary processes and connectivity. Mol Ecol 17: 3775–3787. doi:https://doi.org/10.1111/j.1365-294X.2008.03803.x. PubMed: 18662232.
View Article
Google Scholar

[110] View Article

[111] Google Scholar

[ref42] 42. Knittweis L, Kraemer WE, Timm J, Kochzius M (2008) Genetic structure of Heliofungia actiniformis (Scleractinia: Fungiidae) populations in the Indo-Malay Archipelago: implications for live coral trade management efforts. Conserv Genet 10: 241–249.
View Article
Google Scholar

[113] View Article

[114] Google Scholar

[ref43] 43. Hu ZM, Zhang J, Lopez-Bautista J, Duan DL (2013) Asymmetric genetic exchange in the brown seaweed Sargassum fusiforme (Phaeophyceae) driven by oceanic currents. Mar Biol 160 : 1407–1414. doi:https://doi.org/10.1007/s00227-013-2192-x.
View Article
Google Scholar

[116] View Article

[117] Google Scholar

[ref44] 44. Yeong MLB, Wong CL (2012) Seasonal growth rate of Sargassum species at Teluk Kemang, Port Dickson, Malaysia. J Appl Phycol 25: 805–814.
View Article
Google Scholar

[119] View Article

[120] Google Scholar

[ref45] 45. Tatarenkov A, Bergström L, Jönsson RB, Serrão EA, Kautsky L et al. (2005) Intriguing asexual life in marginal populations of the brown seaweed Fucus vesiculosus. Mol Ecol 14: 647–651. doi:https://doi.org/10.1111/j.1365-294X.2005.02425.x. PubMed: 15660953.
View Article
Google Scholar

[122] View Article

[123] Google Scholar

[ref46] 46. Fraser CI, Spencer HG, Waters JM (2009) Glacial oceanographic contrasts explain phylogeography of Australian bull kelp. Mol Ecol 18: 2287–2296. doi:https://doi.org/10.1111/j.1365-294X.2009.04201.x. PubMed: 19389161.
View Article
Google Scholar

[125] View Article

[126] Google Scholar

[ref47] 47. Macaya EC, Zuccarello G (2010) Genetic structure of the giant kelp Macrocystis pyrifera along the southeastern Pacific. Mar Ecol-Prog Ser 420: 103–112. doi:https://doi.org/10.3354/meps08893.
View Article
Google Scholar

[128] View Article

[129] Google Scholar

[ref48] 48. Payo DA, Leliaert F, Verbruggen H, D’hondt S, Calumpong HP et al. (2013) Extensive cryptic species diversity and fine-scale endemism in the marine red alga Portieria in the Philippines. Proc R Soc Lond B Biol Sci 280: 1753–1760. PubMed: 23269854.
View Article
Google Scholar

[131] View Article

[132] Google Scholar

[ref49] 49. Maggs CA, Castilho R, Foltz D, Henzler C, Jolly MT et al. (2008) Evaluating signatures of glacial refugia for North Atlantic benthic marine taxa. Ecology 89: S108–S122. doi:https://doi.org/10.1890/08-0257.1. PubMed: 19097488.
View Article
Google Scholar

[134] View Article

[135] Google Scholar

[ref50] 50. Dannewitz J (2005) Panmixia in the European eel: a matter of time. Proc R Soc Lond B Biol Sci 272: 1129–1137. doi:https://doi.org/10.1098/rspb.2005.3064. PubMed: 16024374.
View Article
Google Scholar

[137] View Article

[138] Google Scholar