Phylogeography of Ceriops tagal (Rhizophoraceae) in Southeast Asia: the land barrier of the Malay Peninsula has caused population differentiation between the Indian Ocean and South China Sea

Genomic evidence for demographic fluctuations, genetic burdens and adaptive divergence in fourfinger threadfin Eleutheronema rhadinum

Declining populations and bottlenecks lead to the accumulation of deleterious mutations in fish populations. These processes also trigger genetic purging, which is a key genetic factor in reducing the deleterious burdens and increasing population viability. However, there is a lack of empirical evidence on the interaction between demographic history and the genome-wide pattern of deleterious variations. Here, we generated genome resequencing data of Eleutheronema rhadinum from China and Thailand, representing the major distribution of the species' southern regions. E. rhadinum had exceptionally low genome-wide variability and experienced dramatic population expansions followed by continuous declines. The geographical divergence, which occurred ~ 23,000 years ago, shaped different demographic trajectories and generated different regional patterns of deleterious mutations in China and Thailand populations. Several lines of evidence revealed that this geographical pattern of deleterious mutation was driven by the purging of highly deleterious mutations. We showed that purifying selection had inbreeding-associated fitness costs and was more efficient against missense mutations in the Thailand population, which had the lowest genetic burden of homozygous deleterious mutations. Multiple evolutionarily conserved protein domains were disrupted by the loss-of-function mutations, posing a high probability of gene functionality elimination. Moreover, thermal and salinity genes (Trpm3, Nek4, Gtf2f2, Cldn14) were identified in genomic divergence regions of E. rhadinum among China and Thailand populations. Our findings highlight the importance of demographic history factors shaping the geographical patterns of deleterious mutations. The results serve to deepen our understanding of the adaptive evolution and divergence of E. rhadinum with implications for other marine fish.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42995-024-00276-4.

Assessment of the Genetic Diversity and Population Structure of Rhizophora mucronata along Coastal Areas in Thailand

Unique and biodiverse, mangrove ecosystems provide humans with benefits and contribute to coastal protection. Rhizophora mucronata, a member of the Rhizophoraceae family, is prevalent in the mangrove forests of Thailand. R. mucronata's population structure and genetic diversity have received scant attention. Here, we sequenced the entire genome of R. mucronata using 10× Genomics technology and obtained an assembly size of 219 Mb with the N50 length of 542,540 bases. Using 2857 single nucleotide polymorphism (SNP) markers, this study investigated the genetic diversity and population structure of 80 R. mucronata accessions obtained from the mangrove forests in Thailand. The genetic diversity of R. mucronata was moderate (I = 0.573, Ho = 0.619, He = 0.391). Two subpopulations were observed and confirmed from both population structure and principal component analysis (PCA). Analysis of molecular variance (AMOVA) showed that there was more variation within populations than between them. Mean pairwise genetic differentiation (FST = 0.09) showed that there was not much genetic difference between populations. Intriguingly, the predominant clustering pattern in the R. mucronata population did not correspond to the Gulf of Thailand and the Andaman Sea, which are separated by the Malay Peninsula. Several factors could have influenced the R. mucronata genetic pattern, such as hybridization and anthropogenic factors. This research will provide important information for the future conservation and management of R. mucronata in Thailand.

Barrier to Gene Flow of Grey Mangrove Avicennia marina Populations in the Malay Peninsula as Revealed From Nuclear Microsatellites and Chloroplast Haplotypes

Contemporary mangrove forest areas took shape historically and their genetic connectivity depends on sea-faring propagules, subsequent settlement, and persistence in suitable environments. Mangrove species world-wide may experience genetic breaks caused by major land barriers or opposing ocean currents influencing their population genetic structure. For Malay Peninsula, several aquatic species showed strong genetic differentiation between East and West coast regions due to the Sunda shelf flooding since the Last Glacial Maximum. In this study genetic diversity and structure of Avicennia marina populations in Malay Peninsula were assessed using nuclear microsatellite markers and chloroplast sequences. Even though all populations showed identical morphological features of A. marina, three evolutionary significant units were obtained with nuclear and cytoplasmic markers. Avicennia marina along a 586 km stretch of the West coast differed strongly from populations along an 80 km stretch of the East coast featuring chloroplast capture of Avicennia alba in an introgressive A. marina. Over and above this expected East-West division, an intra-regional subdivision was detected among A. marina populations in the narrowest region of the Strait of Malacca. The latter genetic break was supported by an amova, structure, and barrier analysis whereas RST > FST indicated an evolutionary signal of long-lasting divergence. Two different haplotypes along the Western coast showed phylogeographic relationship with either a northern or a putative southern lineage, thereby assuming two Avicennia sources facing each other during Holocene occupation with prolonged separation in the Strait of Malacca. Migrate-n model testing supported a northward unidirectional stepping-stone migration route, although with an unclear directionality at the genetic break position, most likely due to weak oceanic currents. Low levels of genetic diversity and southward connectivity was detected for East coast Avicennia populations. We compared the fine-scale spatial genetic structure (FSGS) of Avicennia populations along the exposed coast in the East vs. the sheltered coast in the West. A majority of transects from both coastlines revealed no within-site kinship-based FSGS, although the remoteness of the open sea is important for Avicennia patches to maintain a neighborhood. The results provide new insights for mangrove researchers and managers for future in-depth ecological-genetic-based species conservation efforts in Malay Peninsula.

Strong genetic differentiation in tropical seagrass Enhalus acoroides (Hydrocharitaceae) at the Indo-Malay Archipelago revealed by microsatellite DNA

The Indo-Malay Archipelago is regarded as a barrier that separates organisms of the Indian and Pacific Oceans. Previous studies of marine biota from this region have found a variety of biogeographic barriers, seemingly dependent on taxon and methodology. Several hypotheses, such as emergence of the Sunda Shelf and recent physical oceanography, have been proposed to account for the genetic structuring of marine organisms in this region. Here, we used six microsatellite loci to infer genetic diversity, population differentiation and phylogeographic patterns of Enhalus acoroides across the Indo-Malay Archipelago. Heterozygosities were consistently high, and significant isolation-by-distance, consistent with restricted gene flow, was observed. Both a neighbour joining tree based on D_A distance and Bayesian clustering revealed three major clusters of E. acoroides. Our results indicate that phylogeographic patterns of E. acoroides have possibly been influenced by glaciation and deglaciation during the Pleistocene. Recent physical oceanography such as the South Java Current and the Seasonally Reversing Current may also play a role in shaping the genetic patterns of E. acoroides.

Effects of Pleistocene sea-level fluctuations on mangrove population dynamics: a lesson from Sonneratia alba

Background

A large-scale systematical investigation of the influence of Pleistocene climate oscillation on mangrove population dynamics could enrich our knowledge about the evolutionary history during times of historical climate change, which in turn may provide important information for their conservation.

Results

In this study, phylogeography of a mangrove tree Sonneratia alba was studied by sequencing three chloroplast fragments and seven nuclear genes. A low level of genetic diversity at the population level was detected across its range, especially at the range margins, which was mainly attributed to the steep sea-level drop and associated climate fluctuations during the Pleistocene glacial periods. Extremely small effective population size (Ne) was inferred in populations from both eastern and western Malay Peninsula (44 and 396, respectively), mirroring the fragility of mangrove plants and their paucity of robustness against future climate perturbations and human activity. Two major genetic lineages of high divergence were identified in the two mangrove biodiversity centres: the Indo-Malesia and Australasia regions. The estimated splitting time between these two lineages was 3.153 million year ago (MYA), suggesting a role for pre-Pleistocene events in shaping the major diversity patterns of mangrove species. Within the Indo-Malesia region, a subdivision was implicated between the South China Sea (SCS) and the remaining area with a divergence time of 1.874 MYA, corresponding to glacial vicariance when the emerged Sunda Shelf halted genetic exchange between the western and eastern coasts of the Malay Peninsula during Pleistocene sea-level drops. Notably, genetic admixture was observed in populations at the boundary regions, especially in the two populations near the Malacca Strait, indicating secondary contact between divergent lineages during interglacial periods. These interregional genetic exchanges provided ample opportunity for the re-use of standing genetic variation, which could facilitate mangrove establishment and adaptation in new habitats, especially in the context of global climate changes.

Conclusion

Phylogeogrpahic analysis in this study reveal that Pleistocene sea-level fluctuations had profound influence on population differentiation of the mangrove tree S. alba. Our study highlights the fragility of mangrove plants and offers a guide for the conservation of coastal mangrove communities experiencing ongoing changes in sea-level.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0849-z) contains supplementary material, which is available to authorized users.

Re-Evaluation of Phylogenetic Relationships among Species of the Mangrove Genus Avicennia from Indo-West Pacific Based on Multilocus Analyses

Avicennia L. (Avicenniaceae), one of the most diverse mangrove genera, is distributed widely in tropical and subtropical intertidal zones worldwide. Five species of Avicennia in the Indo-West Pacific region have been previously described. However, their phylogenetic relationships were determined based on morphological and allozyme data. To enhance our understanding of evolutionary patterns in the clade, we carried out a molecular phylogenetic study using wide sampling and multiple loci. Our results support two monophyletic clades across all species worldwide in Avicennia: an Atlantic-East Pacific (AEP) lineage and an Indo-West Pacific (IWP) lineage. This split is in line with biogeographic distribution of the clade. Focusing on the IWP branch, we reconstructed a detailed phylogenetic tree based on sequences from 25 nuclear genes. The results identified three distinct subclades, (1) A. rumphiana and A. alba, (2) A. officinalis and A. integra, and (3) the A. marina complex, with high bootstrap support. The results strongly corresponded to two morphological traits in floral structure: stigma position in relation to the anthers and style length. Using Bayesian dating methods we estimated diversification of the IWP lineage was dated to late Miocene (c. 6.0 million years ago) and may have been driven largely by the fluctuating sea levels since that time.

Pronounced genetic differentiation and recent secondary contact in the mangrove tree Lumnitzera racemosa revealed by population genomic analyses

Systematically investigating the impacts of Pleistocene sea-level fluctuations on mangrove plants may provide a better understanding of their demographic history and useful information for their conservation. Therefore, we conducted population genomic analyses of 88 nuclear genes to explore the population dynamics of a mangrove tree Lumnitzera racemosa across the Indo-West Pacific region. Our results revealed pronounced genetic differentiation in this species between the populations from the Indian Ocean and the Pacific Ocean, which may be attributable to the long-term isolation between the western and eastern coasts of the Malay Peninsula during sea-level drops in the Pleistocene glacial periods. The mixing of haplotypes from the two highly divergent groups was identified in a Cambodian population at almost all 88 nuclear genes, suggesting genetic admixture of the two lineages at the boundary region. Similar genetic admixture was also found in other populations from Southeast Asia based on the Bayesian clustering analysis of six nuclear genes, which suggests extensive and recent secondary contact of the two divergent lineages in Southeast Asia. Computer simulations indicated substantial migration from the Indian Ocean towards the South China Sea, which likely results in the genetic admixture in Southeast Asia.

Comparative Analysis of the Pattern of Population Genetic Diversity in Three Indo-West Pacific Rhizophora Mangrove Species

Rhizophora species are the most widely distributed mangrove trees in the Indo-West Pacific (IWP) region. Comparative studies of these species with shared life history traits can help identify evolutionary factors that have played most important roles in determining genetic diversity within and between populations in ocean-current dispersed mangrove tree species. We sampled 935 individuals from 54 natural populations for genotyping with 13 microsatellite markers to investigate the level of genetic variation, population structure, and gene flow on a broad geographic scale in Rhizophora apiculata, Rhizophora mucronata, and Rhizophora stylosa across the IWP region. In contrast to the pattern expected of long-lived woody plants with predominant wind-pollination, water-dispersed seeds and wide geographic range, genetic variation within populations was generally low in all the three species, especially in those peripheral populations from geographic range limits. Although the large water-buoyant propagules of Rhizophora have capacity for long distance dispersal, such events might be rare in reality, as reflected by the low level of gene flow and high genetic differentiation between most of population pairs within each species. Phylogeographic separation of Australian and Pacific island populations from SE Asian lineages previously revealed with DNA sequence data was still detectable in R. apiculata based on genetic distances, but this pattern of disjunction was not always evident in R. mucronata and R. stylosa, suggesting that fast-evolving molecular markers could be more suitable for detecting contemporary genetic structure but not deep evolutionary divergence caused by historical vicariance. Given that mangrove species generally have small effective population sizes, we conclude that genetic drift coupled with limited gene flow have played a dominant role in producing the current pattern of population genetic diversity in the IWP Rhizophora species, overshadowing the effects of their life history traits. Recent population fragmentation and disturbances arising from human activities could further endanger genetic diversity in mangrove trees.

Ancient Geographical Barriers Drive Differentiation among Sonneratia caseolaris Populations and Recent Divergence from S. lanceolata

Glacial vicariance is thought to influence population dynamics and speciation of many marine organisms. Mangroves, a plant group inhabiting intertidal zones, were also profoundly influenced by Pleistocene glaciations. In this study, we investigated phylogeographic patterns of a widespread mangrove species Sonneratia caseolaris and a narrowly distributed, closely related species S. lanceolata to infer their divergence histories and related it to historical geological events. We sequenced two chloroplast fragments and five nuclear genes for one population of S. lanceolata and 12 populations of S. caseolaris across the Indo-West Pacific (IWP) region to evaluate genetic differentiation and divergence time among them. Phylogenetic analysis based on sequences of nuclear ribosomal internal transcribed spacer and a nuclear gene rpl9 for all Sonneratia species indicate that S. lanceolata individuals are nested within S. caseolaris. We found strong genetic structure among geographic regions (South China Sea, the Indian Ocean, and eastern Australia) inhabited by S. caseolaris. We estimated that divergence between the Indo-Malesia and Australasia populations occurred 4.035 million years ago (MYA), prior to the onset of Pleistocene. BARRIERS analysis suggested that complex geographic features in the IWP region had largely shaped the phylogeographic patterns of S. caseolaris. Furthermore, haplotype analyses provided convincing evidence for secondary contact of the South China Sea and the Indian Ocean lineages at the Indo-Pacific boundary. Demographic history inference under isolation and migration (IM) model detected substantial gene flow from the Sri Lanka populations to the populations in the Java Island. Moreover, multi-locus sequence analysis indicated that S. lanceolata was most closely related to the Indian Ocean populations of S. caseolaris and the divergence time between them was 2.057 MYA, coinciding with the onset of the Pleistocene glaciation. Our results suggest that geographic isolation driven by the Pleistocene ice age resulted in the recent origin of S. lanceolata.

Phylogeography of the finless porpoise (genus Neophocaena): testing the stepwise divergence hypothesis in the northwestern Pacific

We used 344 mitochondrial control region (717 bp) sequences from the finless porpoise (genus Neophocaena) from the northwestern Pacific to investigate the extent and manner in which past climatic oscillations may have shaped patterns of genetic diversity for this marine mammal. Both SplitsTree and Analysis of Molecular Variance (AMOVA) revealed the presence of a deep divergence among N. phocaenoides in subtropical waters compared with N. asiaeorientalis in temperate waters. Results from Migrate-n indicated that migration increased along the continent during the early Pleistocene period. Migration increased, although to a lesser extent than that during the Pleistocene, along the marginal shelf in the Yellow/Bohai Sea during the Last Glacial Maximum (LGM) due to a shortening coastline. Our results suggest that the current patterns of genetic diversity of Neophocaena vary at a hierarchy on a temporal and spatial scale, and phylogeographic history should be taken into account when examining species population structure and taxonomy.

Inferences of evolutionary history of a widely distributed mangrove species, Bruguiera gymnorrhiza, in the Indo-West Pacific region

Inference of genetic structure and demographic history is fundamental issue in evolutionary biology. We examined the levels and patterns of genetic variation of a widespread mangrove species in the Indo-West Pacific region, Bruguiera gymnorrhiza, using ten nuclear gene regions. Genetic variation of individual populations covering its distribution range was low, but as the entire species it was comparable to other plant species. Genetic differentiation among the investigated populations was high. They could be divided into two genetic clusters: the West and East clusters of the Malay Peninsula. Our results indicated that these two genetic clusters derived from their ancestral population whose effective size of which was much larger compared to the two extant clusters. The point estimate of speciation time between B. gymnorrhiza and Bruguiera sexangula was two times older than that of divergence time between the two clusters. Migration from the West cluster to the East cluster was much higher than the opposite direction but both estimated migration rates were low. The past Sundaland and/or the present Malay Peninsula are likely to prevent gene flow between the West and East clusters and function as a geographical or land barrier.

Recruitment constraints in Singapore's fluted giant clam (Tridacna squamosa) population--a dispersal model approach

Recruitment constraints on Singapore's dwindling fluted giant clam, Tridacna squamosa, population were studied by modelling fertilisation, larval transport, and settlement using real-time hydrodynamic forcing combined with knowledge of spawning characteristics, larval development, behaviour, and settlement cues. Larval transport was simulated using a finite-volume advection-diffusion model coupled to a three-dimensional hydrodynamic model. Three recruitment constraint hypotheses were tested: 1) there is limited connectivity between Singapore's reefs and other reefs in the region, 2) there is limited exchange within Singapore's Southern Islands, and 3) there exist low-density constraints to fertilisation efficacy (component Allee effects). Results showed that connectivity among giant clam populations was primarily determined by residual hydrodynamic flows and spawning time, with greatest chances of successful settlement occurring when spawning and subsequent larval dispersal coincided with the period of lowest residual flow. Simulations suggested poor larval transport from reefs located along the Peninsular Malaysia to Singapore, probably due to strong surface currents between the Andaman Sea and South China Sea combined with a major land barrier disrupting larval movement among reefs. The model, however, predicted offshore coral reefs to the southeast of Singapore (Bintan and Batam) may represent a significant source of larvae. Larval exchange within Singapore's Southern Islands varied substantially depending on the locations of source and sink reefs as well as spawning time; but all simulations resulted in low settler densities (2.1-68.6 settled individuals per 10,000 m(2)). Poor fertilisation rates predicted by the model indicate that the low density and scattered distribution of the remaining T. squamosa in Singapore are likely to significantly inhibit any natural recovery of local stocks.

Molecular data reveal isolation by distance and past population expansion for the shea tree (Vitellaria paradoxa C.F. Gaertn) in West Africa

While the genetic structure of many tree species in temperate, American and Asian regions is largely explained by climatic oscillations and subsequent habitat contractions and expansions, little is known about Africa. We investigated the genetic diversity and structure of shea tree (Vitellaria paradoxa,) in Western Africa, an economically important tree species in the Sudano-Sahelian zone. Eleven nuclear microsatellites (nuc) were used to genotype 673 trees selected in 38 populations. They revealed moderate to high within-population diversity: allelic richness ranged from R(nuc) = 3.99 to 5.63. This diversity was evenly distributed across West Africa. Populations were weakly differentiated (F(STnuc) = 0.085; P < 0.0001) and a pattern of isolation by distance was noted. No phylogeographic signal could be detected across the studied sample. Additionally, two chloroplast microsatellite loci, leading to 11 chlorotypes, were used to analyse a sub-set of 370 individuals. Some variation in chloroplast allelic richness among populations could be detected (R(cp) = 0.00 to 4.36), but these differences were not significant. No trend with latitude and longitude were observed. Differentiation was marked (G(STcp) = 0.553; P < 0.0001), but without a significant phylogeographical signal. Population expansion was detected considering the total population using approximate Bayesian computation (nuclear microsatellites) and mismatch distribution (chloroplast microsatellites) methods. This expansion signal and the isolation by distance pattern could be linked to the past climatic conditions in West Africa during the Pleistocene and Holocene which should have been favourable to shea tree development. In addition, human activities through agroforestry and domestication (started 10,000 bp) have probably enhanced gene flow and population expansion.

Levels and patterns of DNA variation in two sympatric mangrove species, Rhizophora apiculata and R. mucronata from Thailand

In mangrove species the past geomorphic changes in coastal regions and reproductive systems are important factors of their distribution and genetic structure of populations. However, very little is known about the levels of genetic variation of Rhiozophora species in Southeast Asia. In this study, we surveyed levels and patterns of genetic variation as well as population structure of two sympatric mangrove species, Rhizophora apiculata and R. mucronata in Thailand, using five nuclear genes and two cpDNA regions. In all investigated DNA regions, nucleotide variation within species was low, while nucleotide divergence between the two species was considerable. The nuclear genes evolved 10 times faster than the cpDNA regions. In both R. apiculata and R. mucronata, significant positive F(IS) values were found, indicating deviation from Hardy-Weinberg proportions and a deficiency of heterozygotes. In both species, we found significant genetic differentiation between populations. However, the pattern of population differentiation (F(ST)) of R. apiculata differed from that of R. mucronata. Our results suggest that the two investigated species have different demographic history, even though they are sympatric and have similar reproductive systems.

Global diversification of mangrove fauna: a molecular phylogeny of Littoraria (Gastropoda: Littorinidae)

The genus Littoraria is one of very few molluscan groups that are closely associated with mangroves. We document its global evolutionary radiation and compare biogeographic patterns with those of mangrove plants, based on phylogenetic and fossil evidence. Using sequences from three genes (nuclear 28S rRNA, mitochondrial 12S rRNA and COI) we reconstruct a phylogeny of 37 of the 39 living morphospecies. Six monophyletic subgenera are defined (Bulimilittorina, Lamellilitorina, Littoraria, Palustorina, Protolittoraria, Littorinopsis) and we synonymize L. coccinea and L. glabrata. A deep division between Palustorina from the Indo-West Pacific and Littoraria from the Atlantic and Eastern Pacific is estimated by a Bayesian relaxed-clock method to be of Middle Eocene to Palaeocene age (43.2-62.7 Ma), which far predates the Early Miocene (18 Ma) closure of the Tethyan Seaway; this, as in mangrove plants, may reflect vicariance by climatic cooling, rather than tectonic processes. The age of Littoraria angulifera in the Atlantic is, however, consistent with Early Miocene vicariance of a Tethyan ancestor. We infer that speciation events are mainly of Miocene or older age, and that diversification has not been driven by depletion of mangrove habitats during recent glacial intervals. Parsimonious reconstruction of ancestral habitats suggests that the genus has inhabited mangrove or wood substrates since its origin, while the rock-dwelling habit of the four members of Protolittoraria is derived. Three species span the Eastern Pacific Barrier, and one is amphi-Atlantic, consistent with a long larval phase of up to 10 weeks. Allopatric speciation is inferred, but usually with subsequent range overlap. Ovoviviparity (interpreted as an adaptation to life in mangroves) has arisen twice.

Mating system, pollen and propagule dispersal, and spatial genetic structure in a high-density population of the mangrove tree Kandelia candel

Mangrove tree species form ecologically and economically important forests along the tropical and subtropical coastlines of the world. Although low intrapopulation genetic diversity and high interpopulation genetic differentiation have been detected in most mangrove tree species, no direct investigation of pollen and propagule dispersal through paternity and/or parentage analysis and spatial genetic structure within populations has been conducted. We surveyed the mating system, pollen and propagule dispersal, and spatial genetic structure in a natural population of Kandelia candel, one of the typical viviparous mangrove tree species, using nuclear and chloroplast microsatellite markers. High diversity and outcrossing rates were observed. Paternity and parentage analysis and modelling estimations revealed the presence of an extremely short-distance component of pollen and propagule dispersal (pollen: 15.2+/-14.9 m (SD) by paternity analysis and 34.4 m by modelling; propagule: 9.4+/-13.8 m (SD) by parentage analysis, and 18.6 m by modelling). Genetic structure was significant at short distances, and a clumped distribution of chloroplast microsatellite genotypes was seen in K. candel adults. We conclude that the K. candel population was initiated by limited propagule founders from outside by long-distance dispersal followed by limited propagule dispersal from the founders, resulting in a half-sib family structure.

Population genetic structure of three tree species in the mangrove genus Ceriops (Rhizophoraceae) from the Indo West Pacific

Ceriops is a viviparous mangrove with widespread species Ceriops decandra and C. tagal, and an endemic species C. australis. Genetic diversity of the three species was screened in 30 populations collected from 23 locations in the Indo West Pacific (IWP) using Inter-simple sequence repeats (ISSR) and sequences of partial nuclear gene (G3pdh) and chloroplast DNA (trnV-trnM). At the species level, the total gene diversity (Ht) revealed by ISSRs was 0.270, 0.118, and 0.089 in C. decandra, C. tagal, and C. australis, respectively. A total of six haplotypes of G3pdh and five haplotypes of trnV-trnM were recognized among the three species. Only C. decandra was detected containing more than one haplotype from each sequence data set (four G3pdh haplotypes and three trnV-trnM haplotypes). At the population level, genetic diversity of Ceriops was relatively low inferred from ISSRs (He = 0.028, 0.023, and 0.053 in C. decandra, C. tagal, and C. australis, respectively). No haplotype diversity within population was detected from any of the three species. Cluster analysis based on ISSRs identified three major geographical groups in correspond to the East Indian Ocean (EIO), South China Sea (SCS), and North Australia (NA) in both C. decandra and C. tagal. The cladogram from DNA sequences also detected the same three geographical groups in C. decandra. Analysis of molecular variance (AMOVA) revealed that most of the total variation was accounted for by differentiation between the three major geographical regions of both C. decandra and C. tagal. The significant genetic structure may result from the geological events in these regions during the recent Pleistocene glaciations. This study also provided insights into the phylogenetics of Ceriops.