CORRELATED EVOLUTIONARY DIVERGENCE OF EGG SIZE AND A MITOCHONDRIAL PROTEIN ACROSS THE ISTHMUS OF PANAMA

Recalibrating the molecular clock for Arctic marine invertebrates based on DNA barcodes

Divergence times for species assemblages of Arctic marine invertebrates have often been estimated using a standard rate (1.4%/MY) of molecular evolution calibrated using a single sister pair of tropical crustaceans. Because rates of molecular evolution vary among taxa and environments, it is essential to obtain clock calibrations from northern lineages. The recurrent opening and closure of the Bering Strait provide an exceptional opportunity for clock calibration. Here, we apply the iterative calibration approach to investigate patterns of molecular divergence among lineages of northern marine molluscs and arthropods using publicly available sequences of the cytochrome c oxidase subunit I (COI) gene and compare these results with previous estimates of trans-Bering divergences for echinoderms and polychaetes. The wide range of Kimura two-parameter (K2P) divergences among 73 trans-Bering sister pairs (0.12%–16.89%) supports multiple pulses of migration through the Strait. Overall, the results indicate a rate of K2P divergence of 3.2%/MY in molluscs, 5%–5.2%/MY in arthropods, and 3.5%–4.7%/MY in polychaetes. While these rates are considerably higher than the often-adopted 1.4%/MY rate, they are similar to calibrations (3%–5%/MY) in several other studies of marine invertebrates. This upward revision in rates means there is a need both to reevaluate the evolutionary history of marine lineages and to reexamine the impact of prior climatic changes upon the diversification of marine life.

EGG SIZE EVOLUTION IN TROPICAL AMERICAN ARCID BIVALVES: THE COMPARATIVE METHOD AND THE FOSSIL RECORD

Marine organisms exhibit a wide range of egg sizes, even among closely related taxa, and egg size is widely considered to be one of the most important components of the life histories of marine species. The nature of the trade-off between egg size and number and the consequences of variation in egg size for offspring growth and survivorship have been extensively modeled. Yet, there is little empirical evidence that supports the relative importance of particular environmental parameters in engendering the tremendous variation in egg size seen in marine organisms. This study compares egg sizes between six geminate species pairs of bivalves in the family Arcidae to determine whether egg size differs in predictable directions between geminate species in the two oceans separated by the Central American isthmus, and whether the direction and timing of egg size evolution among geminates in this family is correlated with both modern and paleoceanographic patterns of oceanic productivity. In all modern members of six geminate pairs, egg size was larger in the species in the western Atlantic than in its sister species the eastern Pacific. This pattern supports the hypothesis that optimal egg size differs in the two oceans due to the low productivity and poor larval feeding environment in the western Atlantic relative to the eastern Pacific. The fossil record of one geminate pair shows that egg size has remained consistently large in the western Atlantic from the Miocene to the Recent, while egg size in the eastern Pacific has decreased to the current small size in less than 2 million years; this suggests that modern-day differences between egg sizes in the western Atlantic and eastern Pacific are due to either an increase in productivity in the eastern Pacific and subsequent selection for smaller eggs in that ocean, or differential patterns of extinction that occurred well after the rise of the isthmus. These results agree with ancestral character state reconstruction using linear parsimony, but differ from squared-change parsimony reconstructions.

CALIBRATING A MOLECULAR CLOCK FROM PHYLOGEOGRAPHIC DATA: MOMENTS AND LIKELIHOOD ESTIMATORS

We present moments and likelihood methods that estimate a DNA substitution rate from a group of closely related sister species pairs separated at an assumed time, and we test these methods with simulations. The methods also estimate ancestral population size and can test whether there is a significant difference among the ancestral population sizes of the sister species pairs. Estimates presented in the literature often ignore the ancestral coalescent prior to speciation and therefore should be biased upward. The simulations show that both methods yield accurate estimates given sample sizes of five or more species pairs and that better likelihood estimates are obtained if there is no significant difference among ancestral population sizes. The model presented here indicates that the larger than expected variation found in multitaxa datasets can be explained by variation in the ancestral coalescence and the Poisson mutation process. In this context, observed variation can often be accounted for by variation in ancestral population sizes rather than invoking variation in other parameters, such as divergence time or mutation rate. The methods are applied to data from two groups of species pairs (sea urchins and Alpheus snapping shrimp) that are thought to have separated by the rise of Panama three million years ago.