Suzanne T. Williams

SPECIATION AND DIVERSITY ON TROPICAL ROCKY SHORES: A GLOBAL PHYLOGENY OF SNAILS OF THE GENUS ECHINOLITTORINA

Abstract A phylogenetic approach to the origin and maintenance of species diversity ideally requires the sampling of all species within a clade, confirmation that they are evolutionarily distinct entities, and knowledge of their geographical distributions. In the marine tropics such studies have mostly been of fish and reef‐associated organisms, usually with high dispersal. In contrast, snails of the genus Echinolittorina (Littorinidae) are restricted to rocky shores, have a four‐week pelagic development (and recorded dispersal up to 1400 km), and show different evolutionary patterns. We present a complete molecular phylogeny of Echinolittorina, derived from Bayesian analysis of sequences from nuclear 28S rRNA and mitochondrial 12S rRNA and COI genes (nodal support indicated by posterior probabilities, maximum likelihood, and neighbor‐joining bootstrap). This consists of 59 evolutionarily significant units (ESUs), including all 50 known taxonomic species. The 26 ESUs found in the Indo‐West Pacific region form a single clade, whereas the eastern Pacific and Atlantic species are basal. The earliest fossil occurred in the Tethys during the middle Eocene and we suggest that the Indo‐West Pacific clade has been isolated since closure of the Tethyan seaway in the early Miocene. The geographical distributions of all species (based on more than 3700 locality records) appear to be circumscribed by barriers of low temperature, unsuitable sedimentary habitat, stretches of open water exceeding about 1400 km, and differences in oceanographic conditions on the continuum between oceanic and continental. The geographical ranges of sister species show little or no overlap, indicating that the speciation mode is predominantly allopatric. Furthermore, range expansion following speciation appears to have been limited, because a high degree of allopatry is maintained through three to five branching points of the phylogeny. This may be explained by infrequent long‐distance colonization, habitat specialization on the oceanic/continental gradient, and perhaps by interspecific competition. In the eastern Pacific plus Atlantic we identify five cases of divergence on either side of the Isthmus of Panama, but our estimates of their ages pre‐date the emergence of the Isthmus. There are three examples of sister relationships between species in the western Atlantic and eastern Atlantic, all resulting from dispersal to the east. Within the Indo‐West Pacific, we find no geographical pattern of speciation events; narrowly endemic species of recent origin are present in both peripheral and central parts of the region. Evidence from estimated divergence times of sister species, and from a plot of the number of lineages over time, suggest that there has been no acceleration of diversification during the glacioeustatic cycles of the Plio‐Pleistocene. In comparison with reefal organisms, species of Echinolittorina on rocky shores may be less susceptible to extinction or isolation during sea‐level fluctuations. The species richness of Echinolittorina in the classical biogeographic provinces conforms to the common pattern of highest diversity (11 species) in the central “East Indies Triangle of the Indo‐West Pacific, with a subsidiary focus in the eastern Pacific and western Atlantic, and lowest diversity in the eastern Atlantic. The diversity focus in the East Indies Triangle is produced by a mosaic of restricted allopatric species and overlap of a few widespread ones, and is the result of habitat specialization rather than historical vicariance. This study emphasizes the plurality of biogeographic histories and speciation patterns in the marine tropics.

GENETIC STRUCTURE OF GIANT CLAM (TRIDACNA MAXIMA) POPULATIONS IN THE WEST PACIFIC IS NOT CONSISTENT WITH DISPERSAL BY PRESENT-DAY OCEAN CURRENTS

The Pacific marine biota, particularly species with long planktonic larval stages, are thought to disperse widely throughout the Pacific via ocean currents. The little genetic data available to date has supported this view in that little or no significant regional differentiation of populations has been found over large geographical distances. However, recent data from giant clams has demonstrated not only significant regional differentiation of populations, but routes of gene flow that run perpendicular to the main present‐day ocean currents. Extensive surveys of genetic variation at eight polymorphic loci in 19 populations of the giant clam Tridacna maxima, sampled throughout the West and Central Pacific, confirmed that the patterns of variation seen so far in T. gigas were not unique to that species, and may reflect a fundamental genetic structuring of shallow‐water marine taxa. Populations of T. maxima within highly connected reef systems like the Great Barrier Reef were panmictic (average FST < 0.003), but highly significant genetic differences between reef groups on different archipelagos (average FST = 0.084) and between West and Central Pacific regions (average FST = 0.156) were found. Inferred gene flow was high (Nem usually > 5) between the Philippines and the Great Barrier Reef, between the Philippines and Melanesia (the Solomon Islands and Fiji), and between the Philippines and the Central Pacific island groups (Marshall Islands, Kiribati, Tuvalu and Cook Islands). Gene flow was low between these three sets of island chains (Nem < 2). These routes of gene flow are perpendicular to present‐day ocean currents. It is suggested that the spatial patterns of gene frequencies reflect past episodes of dispersal at times of lower sea levels which have not been erased by subsequent dispersal by present‐day circulation. The patterns are consistent with extensive dispersal of marine species in the Pacific, and with traditional views of dispersal from the Indo‐Malay region. However, they demonstrate that dispersal along present‐day ocean surface currents cannot be assumed, that other mechanisms may operate today or that major dispersal events are intermittent (perhaps separated by several thousands of years), and that the nature and timing of dispersal of Pacific marine species is more complex than has been thought.

Did tectonic activity stimulate Oligo-Miocene speciation in the Indo-West Pacific?

Abstract Analyses of molecular phylogenies of three unrelated tropical marine gastropod genera, Turbo, Echinolittorina, and Conus, reveal an increase in the rate of cladogenesis of some Indo-West Pacific (IWP) clades beginning in the Late Oligocene or Early Miocene between 23.7 and 21.0 million years ago. In all three genera, clades with an increased rate of diversification reach a maximum of diversity, in terms of species richness, in the central IWP. Congruence in both the geographical location and the narrow interval of timing suggests a common cause. The collision of the Australia and New Guinea plate with the southeast extremity of the Eurasian plate approximately 25 Mya resulted in geological changes to the central IWP, including an increase in shallow-water areas and length of coastline, and the creation of a mosaic of distinct habitats. This was followed by a period of rapid diversification of zooxanthellate corals between 20 and 25 Mya. The findings reported here provide the first molecular evidence from multiple groups that part of the present-day diversity of shallow-water gastropods in the IWP arose from a rapid pulse of speciation when new habitats became available in the Late Oligocene to Early Miocene. After the new habitats were filled, the rate of speciation likely decreased and this combined with high levels of extinction (in some groups), resulted in a slow down in the rate of diversification in the genera examined.

A molecular phylogeny of heterodont bivalves (Mollusca: Bivalvia: Heterodonta): new analyses of 18S and 28S rRNA genes

A new molecular phylogeny is presented for the highly diverse, bivalve molluscan subclass Heterodonta. The study, the most comprehensive for heterodonts to date, used new sequences of 18S and 28S rRNA genes for 103 species from 49 family groups with species of Palaeoheterodonta (Trigoniidae, Margaritiferidae and Unionidae) as outgroups. Results confirm previous analyses that the Carditidae/Astartidae/Crassatellidae clade is basal to all other heterodonts including Anomalodesmata (often classified as a separate subclass or order). Thyasiroidea occupy a near basal position between the Crassatelloidea and Anomalodesmata. Lucinidae form a well‐supported monophyletic group distinct from Thyasiridae and Ungulinidae. The Solenoidea and Hiatelloidea link as sister groups distant from the Tellinoidea and Myoidea, respectively, where they had been previously associated. The position of the Gastrochaenidae is unstable but does not group with myoidean taxa. Species of four families of Galeommatoidea form a clade that also includes Sportellidae of the Cyamioidea. The Cardioidea and Tellinoidea form highly supported, long branched, individual clades but group as sister taxa. A major clade including Veneroidea, Mactroidea, Myoidea and other families is given the unranked name Neoheterodontei. There is no support for a separate order Myoida (Myoidea and Pholadoidea). Dreissenidae group within the clade including Myidae, Corbulidae, Pholadidae and Teredinidae. The Corbiculoidea is confirmed as polyphyletic with the Sphaeriidae and Corbiculidae forming separate clades within the Neoheterodontei; Corbiculidae grouping with the Glauconomidae. Hemidonacidae are unrelated to the Cardiidae, as previously proposed, but nest within the Neoheterodontei. The Gaimardiidae group near to the Ungulinidae and not with Cyamioidea where most recently classified. The family Ungulinidae, previously classified in the Lucinoidea, forms a well‐supported clade within the Neoheterodontei and is elevated to superfamily rank — Ungulinoidea. The monophyletic status of Glossoidea, Arcticoidea and Veneroidea is unconfirmed. A brief review of the fossil record of the heterodonts indicates that the basal clades of Crassatelloidea, Anomalodesmata and Lucinoidea diverged very early in the Lower Palaeozoic. Other groups such as the Hiatelloidea, Solenoidea, Gastrochaenidae probably were of late Palaeozoic origins. The Cardioidea and Tellinoidea originated in the Triassic while major groups of Neoheterodontei radiated in the Late Mesozoic. The phylogenetic position of the Thyasiroidea and Galeommatoidea suggests a longer fossil history than has so far been recognized.

Indo-West Pacific patterns of genetic differentiation in the high-dispersal starfish Linckia laevigata

Genetic variation in four natural populations of the starfish Linckia laevigata from the Indo‐West Pacific was examined using restriction fragment analysis of a portion of the mtDNA including the control region. Digestion with seven restriction enzymes identified 47 haplotypes in a sample of 326 individuals. Samples collected from reef sites within each location were not significantly differentiated based on ΦST or spatial distribution of haplotypes, indicating that dispersal is high over short to moderate distances. Evidence of gene flow is further supported by the low divergence among haplotypes and the lack of any clear geographical structuring among different haplotypes in the gene phylogeny. However, analysis of molecular variance (AMOVA), ΦST and contingency χ2 analyses of the spatial distribution of haplotypes demonstrate the presence of significant broad scale population genetic structure among the four widespread locations examined. RFLP data are consistent with high gene flow between the Philippines and Western Australia and moderate gene flow between the Great Barrier Reef (GBR) and Fiji, but only limited gene flow between either the Philippines or Western Australia and either the GBR or Fiji. The presence of mtDNA structure contrasts with previous allozyme data which suggest that dispersal among widely separated locations is equivalent to dispersal among populations within the highly connected GBR studies. This discordance between patterns of gene flow inferred from these two markers cannot be fully accounted for by differences in effective population size for mtDNA. This might suggest that while mtDNA variation may represent contemporary patterns of gene flow, allozyme variation among populations is yet to reach equilibrium between drift and migration over the range surveyed.

The Marine Indo-West Pacific Break: Contrasting the Resolving Power of Mitochondrial and Nuclear Genes

Abstract Simultaneous studies of both nuclear and mitochondrial markers were undertaken in two widespread Indo-West Pacific (IWP) marine invertebrates to compare and contrast the ability of these markers to resolve genetic structure. In particular, we were interested in the resolution of a genetic break between the Indian and Pacific Oceans due to historical isolation. Sequence variation from the nuclear gene encoding myosin heavy chain (MyHC) and the mitochondrial gene cytochrome oxidase I (COI) were examined for the snapping shrimp Alpheus lottini from wide-ranging populations throughout the Indian and Pacific Oceans. A previously identified genetic break between oceans based on COI sequences appears to have been an artifact caused by the inadvertent inclusion of pseudogene sequences; our new COI data provide evidence only of a break between IWP and East Pacific populations. Distribution of a single nucleotide polymorphism in MyHC, on the other hand, shows evidence of a cline between Indian and Pacific Oceans. New allozyme and mtDNA sequence data were also obtained for the starfish Linckia laevigata. Allozyme data show a clear genetic break between Indian Ocean populations and Pacific (including western Australian) populations, whereas the distribution of mtDNA haplotypes shows a region of overlap in the central IWP. Comparisons of our data for both Alpheus and Linckia with data from other population genetic studies in the IWP suggest that nuclear markers (allozymes, sequence data and morphological characters) may in some instances reveal historical patterns of genetic population structure whereas mtDNA variation better reflects present day patterns of gene flow.

Genetic consequences of long larval life in the starfish Linckia laevigata (Echinodermata: Asteroidea) on the great barrier reef

Gene flow between populations of the asteroid Linckia laevigata (Linnaeus) was investigated by examining over 1000 individuals collected from ten reefs throughout the Great Barrier Reef (GBR), Australia, for genetic variation at seven polymorphic enzyme loci. Despite geographic separations in excess of 1000 km, Neis unbiased genetic distance (0 to 0.003) and standardised genetic variation between populations (FST) values (mean 0.0011) were small and not significant. Genetic homogeneity among L. laevigata populations is consistent with the long-distance dispersal capability of its 28 d planktonic larval phase, and is greater than that observed for other asteroid species, including another high-dispersal species, Acanthaster planci, which has a 14 d larval phase. Variation within populations was also higher than previously recorded for asteroids (mean heterozygosity=0.384; number of alleles per locus ranged from 5.1 to 6.0 in each population). Among asteroids, dispersal ability is positively correlated with gene flow and levels of variation, and negatively correlated with levels of differentiation.

Molecular systematics of Vetigastropoda: Trochidae, Turbinidae and Trochoidea redefined

Trochoidea are a large superfamily of morphologically and ecologically diverse marine gastropods. We present here an appraisal of the composition and relationships among trochoidean families based on molecular data, with an especial focus on the family Trochidae. Bayesian analyses of sequences from three genes (18S rRNA, 28S rRNA and COI) including data from 162 vetigastropod species show that the gastropod family Trochidae (sensu Hickman & McLean (1990 ), Natural History Museum Los Angeles County Science Series, 35, 1–169) is not monophyletic. Recognition of Chilodontidae, Solariellidae and Calliostomatidae at the family level is supported. Our new, more limited, definition of Trochidae includes the subfamilies Stomatellinae, Lirulariinae and Umboniinae and redefined Trochinae, Cantharidinae and Monodontinae. Halistylinae are provisionally retained in the Trochidae based on previous morphological studies. As redefined, Trochidae are a predominantly shallow‐water radiation in the tropics and subtropics. Some subfamilies and genera previously included in Trochidae have been moved to an enlarged family Turbinidae. The family Turbinidae has been redefined to include Turbininae, Skeneinae, Margaritinae, Tegulinae, Prisogasterinae and most surprisingly the commercially important genus Tectus Montfort, 1810. The new definition of Turbinidae means that the family includes both predominantly shallow and deep‐water clades as well as genera that are distributed across the globe from the poles to the tropics. A greater range of habitat is now seen in Turbinidae than in Trochidae. The redefined Trochidae and Turbinidae, together with Solariellidae, Calliostomatidae and Liotiidae, make up the superfamily Trochoidea. Phasianellidae and Colloniidae are recognized as belonging in a new superfamily, Phasianelloidea, and Angaria Röding, 1798 is recognized as belonging in a new superfamily, Angarioidea. Placement of Areneidae into a superfamily awaits further work.

The complete mitochondrial genome of a turbinid vetigastropod from MiSeq Illumina sequencing of genomic DNA and steps towards a resolved gastropod phylogeny.

A need to increase sampling of mitochondrial genomes for Vetigastropoda has been identified as an important step towards resolving relationships within the Gastropoda. We used shotgun sequencing of genomic DNA, using an Illumina MiSeq, to obtain the first mitochondrial genome for the vetigastropod family Turbinidae, doubling the number of genomes for the species-rich superfamily Trochoidea. This method avoids the necessity of finding suitable primers for long PCRs or primer-walking amplicons, resulting in a timely and cost-effective method for obtaining whole mitochondrial genomes from ethanol-preserved tissue samples. Bayesian analysis of amino acid variation for all available gastropod genomes including the new turbinid mtgenome produced a well resolved tree with high nodal support for most nodes. Major clades within Gastropoda were recovered with strong support, with the exception of Littorinimorpha, which was polyphyletic. We confirm here that mitogenomics is a useful tool for molluscan phylogenetics, especially when using powerful new models of amino acid evolution, but recognise that increased taxon sampling is still required to resolve existing differences between nuclear and mitochondrial gene trees.

A molecular phylogenetic framework for the Muricidae, a diverse family of carnivorous gastropods.

With over 1600 extant described species, the Muricidae are one of the most species-rich and morphologically diverse families of molluscs. As predators of molluscs, polychaetes, anthozoans barnacles and other invertebrates, they form an important component of many benthic communities. Traditionally, the classification of muricids at specific and generic levels has been based primarily on shells, while subfamilies have been defined largely by radular morphology, although the composition and relationships of suprageneric groups have never been studied exhaustively. Here we present the phylogenetic relationships of 77 muricid species belonging to nine of the ten currently recognized subfamilies, based on Bayesian inference and Maximum Likelihood analyses of partial sequences of three mitochondrial (12S, 16S and COI) and one nuclear (28S) genes. The resulting topologies are discussed with respect to traditional subfamilial arrangements, and previous anatomical and molecular findings. We confirm monophyly of each of the subfamilies Ergalataxinae, Rapaninae, Coralliophilinae, Haustrinae, Ocenebrinae and Typhinae as previously defined, but earlier concepts of Muricinae, Trophoninae and Muricopsinae are shown to be polyphyletic. Based on our phylogenetic hypothesis, a new arrangement of these subfamilies is proposed.