Abstract
Widely distributed, broadcast-spawning Diadema sea urchins have been used as model in-vertebrate species for studying the zoogeography of the tropical Indo-Pacific. So far, the Indo-Malay archipelago, a wide and geographically complex maritime region extending from the eastern Indian Ocean to the western Pacific Ocean, has been under-sampled. This study aims to fill this sampling gap and uncover the phylogeographic structure of the long-spined sea-urchin D. setosum in the central Indo-West pacific region. D. setosum samples (total N = 718) were collected in 13 sites throughout the Indo-Malay archipelago. We se-quenced 1157 base pairs of the COI gene. The phylogeographic structure was derived from pairwise Фst estimates using multidimensional scaling and hierarchical clustering analysis; biogeographic hypotheses were tested by AMOVA; genetic relationships between haplo-types were summarised in the form of a minimum-spanning network; pairwise mismatch distributions were compared to the expectations from demographic and spatial expansion models. All samples from the Indo-West Pacific were of the previously uncovered D. se-tosum-a lineage. Phylogeographic structure was evident: the Andaman Sea population and the northern New Guinea population were genetically distinct. Subtler but significant hap-lotype-frequency differences distinguished two populations within the Indonesian seas, distributed in a parapatric-like fashion. The phylogeographic partition observed was insuffi-ciently explained by previous biogeographic hypotheses. The haplotype network showed a series of closely related star-shaped haplogroups with a high proportion of singletons. Nu-cleotide-pairwise mismatch patterns in the two populations from the Indonesian seas were consistent with both demographic and spatial expansion models. While geographic barriers to gene flow were inferred at the western and eastern extremities of the Indo-Malay archi-pelago, the subtler parapatric pattern observed within the Indonesian seas indicated re-striction in gene flow, in a fashion that can hardly be explained by geographic isolation giv-en the dynamic current systems that cross this region. Our results thus raise the hypothesis of subtle reproductive isolation between ecologically incompatible populations. While the coalescence pattern of the Andaman-Sea population suggested demographic stability over evolutionary timescales, that of the two populations from the Indonesian seas indicated recent population expansion, possibly linked to the rapid changes in available D. setosum habitat caused by sea-level oscillations in the late Pleistocene. The phylogeographic pat-terns observed in this study point to likely allopatric differentiation in the central Indo-West Pacific region. Genetic differences between populations were likely reinforced during inter-glacials by some form of reproductive isolation.