Harilaos A. Lessios

POPULATION STRUCTURE AND SPECIATION IN TROPICAL SEAS: GLOBAL PHYLOGEOGRAPHY OF THE SEA URCHIN DIADEMA

Abstract.— The causes of speciation in the sea are rarely obvious, because geographical barriers are not conspicuous and dispersal abilities or marine organisms, particularly those of species with planktonic larvae, are hard to determine. The phylogenetic relations of species in cosmopolitan genera can provide information on the likely mode of their formation. We reconstructed the phylogeny of the pantropical and subtropical sea urchin genus Diadema, using sequences of mitochondrial DNA from 482 individuals collected around the world, to determine the efficacy of barriers to gene flow and to ascertain the history of possible dispersal and vicariance events that led to speciation. We also compared 22 isozyme loci between all described species except D. palmeri. The mitochondrial DNA data show that the two deepest lineages are found in the Indian and West Pacific Oceans. (Indo‐Pacific) Diadema setosum diverged first from all other extant Diadema, probably during the initiation of wide fluctuations in global sea levels in the Miocene. The D. setosum clade then split 3‐5 million years ago into two clades, one found around the Arabian Peninsula and the other in the Indo‐West Pacific. On the lineage leading to the other species of Diadema, the deepest branch is composed of D. palmeri, apparently separated when the climate of New Zealand became colder and other tropical echinoids at these islands went extinct. The next lineage to separate is composed of a currently unrecognized species of Diadema that is found at Japan and the Marshall Islands. Diadema mexicanum in the eastern Pacific separated next, whereas D. paucispinum, D. savignyi, and D. antillarum from the western and central Atlantic, and (as a separate clade) D. antillarum from the eastern Atlantic form a shallow polytomy. Apparently, Indo‐Pacific populations of Diadema maintained genetic contact with Atlantic ones around the southern tip of Africa for some time after the Isthmus of Panama was complete. Diadema paucispinum contains two lineages: D. paucispinum sensu stricto is not limited to Hawaii as previously thought, but extends to Easter Island, Pitcairn, and Okinawa; A second mitochondrial clade of D. paucispinum extends from East Africa and Arabia to the Philippines and New Guinea. A more recent separation between West Indian Ocean and West Pacific populations was detected in D. setosum. Presumably, these genetic discontinuities are the result of water flow restrictions in the straits between northern Australia and Southeast Asia during Pleistocene episodes of low sea level. Diadema savignyi is characterized by high rates of gene flow from Kiribati in the central Pacific all the way to the East African Coast. In the Atlantic, there is a biogeographic barrier between the Caribbean and Brazil, possibly caused by fresh water outflow from the Amazon and the Orinoco Rivers. Diadema antillarum populations of the central Atlantic islands of Ascension and St. Helena are genetically isolated and phylogenetically derived from Brazil. Except for its genetic separation by the mid‐Atlantic barrier, Diadema seems to have maintained connections through potential barriers to dispersal (including the Isthmus of Panama) more recently than did Eucidaris or Echinometra, two other genera of sea urchins in which phylogeography has been studied. Nevertheless, the mtDNA phylogeography of Diadema includes all stages expected from models of allopatric differentiation. There are anciently separated clades that now overlap in their geographic distribution, clades isolated in the periphery of the genus range that have remained in the periphery, clades that may have been isolated in the periphery but have since spread towards the center, closely related clades on either side of an existing barrier, and closely related monophyletic entities on either side of an historical barrier that have crossed the former barrier line, but have not attained genetic equilibrium. Except for D. paucispinum and D. savignyi, in which known hybridization may have lodged mtDNA from one species into the genome of the other, closely related clades are always allopatric, and only distantly related ones overlap geographically. Thus, the phylogenetic history and distribution of extant species of Diadema is by and large consistent with allopatric speciation.

PHYLOGEOGRAPHY OF THE PANTROPICAL SEA URCHIN TRIPNEUSTES: CONTRASTING PATTERNS OF POPULATION STRUCTURE BETWEEN OCEANS

Abstract To understand how allopatric speciation proceeds, we need information on barriers to gene flow, their antiquity, and their efficacy. For marine organisms with planktonic larvae, much of this information can only be obtained through the determination of divergence between populations. We evaluated the importance of ocean barriers by studying the mitochondrial DNA phylogeography of Tripneustes, a pantropical genus of shallow water sea urchin. A region of cytochrome oxidase I (COI) was sequenced in 187 individuals from locations around the globe. The COI phylogeny agreed with a previously published phylogeny of bindin that barriers important to the evolution of Tripneustes are: (1) the cold water upwelling close to the tip of South Africa, (2) the Isthmus of Panama, (3) the long stretch of deep water separating the eastern from the western Atlantic, and (4) the freshwater plume of the Orinoco and the Amazon rivers between the Caribbean and the coast of Brazil. These barriers have previously been shown to be important in at least a subset of the shallow water marine organisms in which phylogeography has been studied. In contrast, the Eastern Pacific Barrier, 5000 km of deep water between the central and the eastern Pacific that has caused the deepest splits in other genera of sea urchins, is remarkably unimportant as a cause of genetic subdivision in Tripneustes. There is also no discernible subdivision between the Pacific and Indian Ocean populations of this genus. The most common COI haplotype is found in the eastern, central, and western Pacific as well as the Indian Ocean. Morphology, COI, and bindin data agree that T. depressus from the eastern Pacific and T. gratilla from the western Pacific are, in fact, the same species. The distribution of haplotype differences in the Indo‐Pacific exhibits characteristics expected from a sea urchin genus with ephemeral local populations, but with high fecundity, dispersal, and growth: there is little phylogenetic structure, and mismatch distributions conform to models of recent population expansion on a nearly global scale. Yet, comparisons between local populations produce large and significant FST values, indicating nonrandom haplotype distribution. This apparent local differentiation is only weakly reflected in regional divergence, and there is no evidence of isolation by distance in correlations between FST values and either geographical or current distance. Thus, Tripneustes in the Indo‐Pacific (but not in the Atlantic) seems to be one large metapopulation spanning two oceans and containing chaotic, nonequilibrium local variation, produced by the haphazard arrival of larvae or by unpredictable local extinction.

Phylogeography of the Pantropical Sea Urchin Eucidaris in Relation to Land Barriers and Ocean Currents

The pantropical sea urchin genus Eucidaris contains four currently recognized species, all of them allopatric: E. metularia in the Indo‐West Pacific, E. thouarsi in the eastern Pacific, E. tribuloides in both the western and eastern Atlantic, and E. clavata at the central Atlantic islands of Ascension and St. Helena. We sequenced a 640‐bp region of the cytochrome oxidase I (COI) gene of mitochondrial DNA to determine whether this division of the genus into species was confirmed by molecular markers, to ascertain their phylogenetic relations, and to reconstruct the history of possible dispersal and vicariance events that led to present‐day patterns of species distribution. We found that E. metularia split first from the rest of the extant species of the genus. If COI divergence is calibrated by the emergence of the Isthmus of Panama, the estimated date of the separation of the Indo‐West Pacific species is 4.7–6.4 million years ago. This date suggests that the last available route of genetic contact between the Indo‐Pacific and the rest of the tropics was from west to east through the Eastern Pacific Barrier, rather than through the Tethyan Sea or around the southern tip of Africa. The second cladogenic event was the separation of eastern Pacific and Atlantic populations by the Isthmus of Panama. Eucidaris at the outer eastern Pacific islands (Galapagos, Isla del Coco, Clipperton Atoll) belong to a separate clade, so distinct from mainland E. thouarsi as to suggest that this is a different species, for which the name E. galapagensis is revived from the older taxonomic literature. Complete lack of shared alleles in three allozyme loci between island and mainland populations support their separate specific status. Eucidaris galapagensis and E. thouarsi are estimated from their COI divergence to have split at about the same time that E. thouarsi and E. tribuloides were being separated by the Isthmus of Panama. Even though currents could easily convey larvae between the eastern Pacific islands and the American mainland, the two species do not appear to have invaded each others ranges. Conversely, the central Atlantic E. clavata at St. Helena and Ascension is genetically similar to E. tribuloides from the American and African coasts. Populations on these islands are either genetically connected to the coasts of the Atlantic or have been colonized by extant mitochondrial DNA lineages of Eucidaris within the last 200,000 years. Although it is hard to explain how larvae can cross the entire width of the Atlantic within their competent lifetimes, COI sequences of Eucidaris from the west coast of Africa are very similar to those of E. tribuloides from the Caribbean. FST statistics indicate that gene flow between E. metularia from the Indian Ocean and from the western and central Pacific is restricted. Low gene flow is also evident between populations of E. clavata from Ascension and St. Helena. Rates of intraspecific exchange of genes in E. thouarsi, E. galapagensis, and E. tribuloides, on the other hand, are high. The phylogeny of Eucidaris confirms Ernst Mayrs conclusions that major barriers to the dispersal of tropical echinoids have been the wide stretch of deep water between central and eastern Pacific, the cold water off the southwest coast of Africa, and the Isthmus of Panama. It also suggests that a colonization event in the eastern Pacific has led to speciation between mainland and island populations.

Crossing the impassable: genetic connections in 20 reef fishes across the eastern Pacific barrier

The ‘impassable’ Eastern Pacific Barrier (EPB), ca 5000 km of deep water separating the eastern from the central Pacific, is the Worlds widest marine biogeographic barrier. Sequencing of mitochondrial DNA in 20 reef fish morphospecies encountered on both sides of the barrier revealed cryptic speciation in two. Among the other 18 species only two showed significant differentiation (as revealed by haplotype networks and FST statistics) between the eastern and the central Pacific. Coalescence analyses indicated that genetic similarity in the 18 truly transpacific species resulted from different combinations of ages of most recent invasion and of levels of recurrent gene flow, with estimated times of initial separation ranging from approximately 30 000 to 1 Myr (ago). There is no suggestion of simultaneous interruptions of gene flow among the species. Migration across the EPB was previously thought to be exclusively eastward, but our evidence showed two invasions from east to west and eight cases in which subsequent gene flow possibly proceeded in the same direction. Thus, the EPB is sporadically permeable to propagules originating on either side.

Dispersal barriers in tropical oceans and speciation in Atlantic and eastern Pacific sea urchins of the genus Echinometra

Echinometra is a pantropical sea urchin made famous through studies of phylogeny, speciation, and genetic structure of the Indo‐West Pacific (IWP) species. We sequenced 630 bp of the cytochrome oxidase I (COI) mitochondrial gene to provide comparable information on the eastern Pacific and Atlantic species, using divergence between those separated by closure of the Isthmus of Panama 3.1 million years ago (Ma) to estimate dates for cladogenic events. Most recently (1.27–1.62 Ma), the Atlantic species E. lucunter and E. viridis diverged from each other, at a time in the Pleistocene that sea levels fell and Caribbean coral speciation and extinction rates were high. An earlier split, assumed to have been coincident with the completion of the Isthmus of Panama, separated the eastern Pacific E. vanbrunti from the Atlantic common ancestor. Transisthmian COI divergence similar to that in the sea urchin genus Eucidaris supports this assumption. The most ancient split in Echinometra occurred between the IWP and the neotropical clades, due to cessation of larval exchange around South Africa or across the Eastern Pacific Barrier. Gene flow within species is generally high; however, there are restrictions to genetic exchange between E. lucunter populations from the Caribbean and those from the rest of the Atlantic. Correlation between cladogenic and vicariant events supports E. Mayr’s contention that marine species, despite their high dispersal potential, form by means of geographical separation. That sympatric, nonhybridizing E. lucunter and E. viridis were split so recently suggests, however, that perfection of reproductive barriers between marine species with large populations can occur in less than 1.6 million years (Myr).

Massive gene flow across the world's most potent marine biogeographic barrier

The ‘Eastern Pacific Barrier’ (EPB), 5400 km of uninterrupted deep water between the central and eastern Pacific, constitutes the greatest marine obstacle to the dispersal of shallow–water organisms. However, some species are found on both sides of the EPB. These ‘transpacific’ species are considered by ‘dispersal’ biogeographers as evidence of invasions through the barrier. ‘Vicariance’ biogeographers, on the other hand, think that transpacific species are morphologically conservative remnants of previously continuous distributions. We compared nucleotide sequences in a 642 bp region of mitochondrial DNA, and electrophoretically detected alleles in 17 enzymatic loci of central and eastern Pacific populations of Echinothrix diadema, an Indo–Pacific sea urchin recently reported from the eastern Pacific. Both types of molecules produced clear evidence of massive, recent gene flow across the EPB. Thus, rather than being isolated relicts of Tethyan distributions, conspecific populations from the eastern and central Pacific are genetically connected. Though the EPB is biogeographically important as a cause of speciation in many groups, it allows genetic connections in others, possibly through larval transport during El Niño events.

TEST FOR SIMULTANEOUS DIVERGENCE USING APPROXIMATE BAYESIAN COMPUTATION

Abstract Comparative phylogeographic studies often reveal disparate levels of sequence divergence between lineages spanning a common geographic barrier, leading to the conclusion that isolation was nonsynchronous. However, only rarely do researchers account for the expected variance associated with ancestral coalescence and among-taxon variation in demographic history. We introduce a flexible approximate Bayesian computational (ABC) framework that can test for simultaneous divergence (TSD) using a hierarchical model that incorporates idiosyncratic differences in demographic history across taxon pairs. The method is tested across a range of conditions and is shown to be accurate even with single-locus mitochondrial DNA (mtDNA) data. We apply this method to a landmark dataset of putative simultaneous vicariance, eight geminate echinoid taxon pairs thought to have been split by the Isthmus of Panama 3.1 million years ago. The ABC posterior estimates are not consistent with a history of simultaneous vicariance given these data. Subsequent ABC estimates under a constrained model that assumes two divergence times across the eight taxon pairs suggests simultaneous divergence 3.1 million years ago in seven of the taxon pairs and a more recent divergence in the remaining taxon pair. These ABC estimates on the simultaneous divergence of the seven taxon pairs correspond to a DNA substitution rate of approximately 1.59% per lineage per million years at the mtDNA cytochrome oxidase I gene. This ABC framework can easily be modified to analyze single taxon-pair datasets and/or be expanded to include multiple loci, migration, recombination, and other idiosyncratic demographic histories. The flexible aspect of ABC and its built-in evaluation of estimator bias and statistical power has the potential to greatly enhance statistical rigor in phylogeographic studies.

Formation of the Isthmus of Panama

Independent evidence from rocks, fossils, and genes converge on a cohesive narrative of isthmus formation in the Pliocene. The formation of the Isthmus of Panama stands as one of the greatest natural events of the Cenozoic, driving profound biotic transformations on land and in the oceans. Some recent studies suggest that the Isthmus formed many millions of years earlier than the widely recognized age of approximately 3 million years ago (Ma), a result that if true would revolutionize our understanding of environmental, ecological, and evolutionary change across the Americas. To bring clarity to the question of when the Isthmus of Panama formed, we provide an exhaustive review and reanalysis of geological, paleontological, and molecular records. These independent lines of evidence converge upon a cohesive narrative of gradually emerging land and constricting seaways, with formation of the Isthmus of Panama sensu stricto around 2.8 Ma. The evidence used to support an older isthmus is inconclusive, and we caution against the uncritical acceptance of an isthmus before the Pliocene.

Mass mortality of Diadema antillarum on the Caribbean coast of Panama

The ecologically important sea urchin Diadema antillarum suffered mass mortalities in 1983, first noted in Panama and then reported from the rest of the Caribbean. We documented the effects of this mortality at two localities on the Atlantic coast of Panama, Punta Galeta and the San Blas Archipelago. At Punta Galeta, affected by the mortality in January 1983, the numbers of D. antillarum changed from an estimated 14,000 per ha in June 1982 to 0.5 per ha in May 1983; by February 1984 they had increased to 38 per ha. In the San Blas, where mass mortality started in April 1983, the number of D. antillarum in permanent quadrats on 8 reefs was reduced by an average of 94.2%. The average reduction in population density measured in transects on nine reefs was 98.9%. Data taken in permanent quadrats on four reefs in 1978, 1979 and 1980 indicate that population fluctuations of D. antillarum are normally much smaller, justifying the labeling of the 1983 event as “mass mortality”. Size structure of the San Blas populations was also affected; mean test diameter of D. antillarum on four reefs was reduced from 48.6 mm to 25.0 mm. Other echinoids (Echinometra viridis, E. lucunter, Lytechinus variegatus, L. williamsi, Eucidaris tribuloides, Tripneustes ventricosus, Clypeaster rosaceus and Echinoneus cyclostomus) suffered no ill effects at either Galeta or the San Blas; their population densities remained stable or increased. Density determinations of Diadema mexicanum at the island of Taboguilla on the Pacific side of Panama indicate that Diadema mass mortality did not extend to the eastern Pacific. Sea surface temperatures, tidal levels, rainfall and salinity showed no abnormal fluctuations during the time of D. antillarum mass mortality at Galeta, suggesting that mortality was not due to physical stress. The wide geographical spread and species-specificity of the mortality suggest a water-borne pathogen as the most likely causative agent. Recovery of D. antillarum populations is likely to be slow because there are few, if any, unaffected populations in the Caribbean to contribute larvae for the recolonization of depleted areas. The absence of D. antillarum will probably be reflected by changes in the algal, coral and echinoid communities, and by altered patterns of bioerosion.

A threat to coral reefs multiplied? Four species of crown-of-thorns starfish

In the face of ever-increasing threats to coral reef ecosystems, it is essential to understand the impact of natural predators in order to devise appropriate management strategies. Destructive population explosions of the crown-of-thorns starfish Acanthaster planci have devastated coral reefs throughout the Indo-Pacific for decades. But despite extensive research, the causes of outbreaks are still unclear. An important consideration in this research is that A. planci has been regarded as a single taxonomic entity. Using molecular data from its entire distribution, we find that A. planci is in fact a species complex. This discovery has important consequences for future coral reef research, and might prove critical for successful reef conservation management.