Peter B. Marko

The ‘Expansion–Contraction’ model of Pleistocene biogeography: rocky shores suffer a sea change?

Approximately 20 000 years ago the last glacial maximum (LGM) radically altered the distributions of many Northern Hemisphere terrestrial organisms. Fewer studies describing the biogeographic responses of marine species to the LGM have been conducted, but existing genetic data from coastal marine species indicate that fewer taxa show clear signatures of post‐LGM recolonization. We have assembled a mitochondrial DNA (mtDNA) data set for 14 co‐distributed northeastern Pacific rocky‐shore species from four phyla by combining new sequences from ten species with previously published sequences from eight species. Nuclear sequences from four species were retrieved from GenBank, plus we gathered new elongation factor 1‐α sequences from the barnacle Balanus glandula. Results from demographic analyses of mtDNA for five (36%) species (Evasterias troschelii, Pisaster ochraceus, Littorina sitkana, L. scutulata, Xiphister mucosus) were consistent with large population expansions occurring near the LGM, a pattern expected if these species recently recolonized the region. However, seven (50%) species (Mytilus trossulus, M. californianus, B. glandula, S. cariosus, Patiria miniata, Katharina tunicata, X. atropurpureus) exhibited histories consistent with long‐term stability in effective population size, a pattern indicative of regional persistence during the LGM. Two species of Nucella with significant mtDNA genetic structure showed spatially variable demographic histories. Multilocus analyses for five species were largely consistent with mtDNA: the majority of multilocus interpopulation divergence times significantly exceeded the LGM. Our results indicate that the LGM did not extirpate the majority of species in the northeastern Pacific; instead, regional persistence during the LGM appears a common biogeographic history for rocky‐shore organisms in this region.

The complex analytical landscape of gene flow inference

Gene flow estimation is essential for characterizing local adaptation, speciation potential and connectivity among threatened populations. New model-based population genetic methods can resolve complex demographic histories, but many studies in fields such as landscape genetics continue to rely on simple rules of thumb focused on gene flow to explain patterns of spatial differentiation. Here, we show how methods that use gene genealogies can reveal cryptic demographic histories and provide better estimates of gene flow with other parameters that contribute to genetic variation across landscapes and seascapes. We advocate for the expanded use and development of methods that consider spatial differentiation as the product of multiple forces interacting over time, and caution against a routine reliance on post-hoc gene flow interpretations.

The cytochrome P450 (CYP) gene superfamily in Daphnia pulex.

BackgroundCytochrome P450s (CYPs) in animals fall into two categories: those that synthesize or metabolize endogenous molecules and those that interact with exogenous chemicals from the diet or the environment. The latter form a critical component of detoxification systems.ResultsData mining and manual curation of the Daphnia pulex genome identified 75 functional CYP genes, and three CYP pseudogenes. These CYPs belong to 4 clans, 13 families, and 19 subfamilies. The CYP 2, 3, 4, and mitochondrial clans are the same four clans found in other sequenced protostome genomes. Comparison of the CYPs from D. pulex to the CYPs from insects, vertebrates and sea anemone (Nematostella vectensis) show that the CYP2 clan, and to a lesser degree, the CYP4 clan has expanded in Daphnia pulex, whereas the CYP3 clan has expanded in insects. However, the expansion of the Daphnia CYP2 clan is not as great as the expansion observed in deuterostomes and the nematode C. elegans. Mapping of CYP tandem repeat regions demonstrated the unusual expansion of the CYP370 family of the CYP2 clan. The CYP370s are similar to the CYP15s and CYP303s that occur as solo genes in insects, but the CYP370s constitute ~20% of all the CYP genes in Daphnia pulex. Lastly, our phylogenetic comparisons provide new insights into the potential origins of otherwise mysterious CYPs such as CYP46 and CYP19 (aromatase).ConclusionOverall, the cladoceran, D. pulex has a wide range of CYPs with the same clans as insects and nematodes, but with distinct changes in the size and composition of each clan.

Divergence genetics analysis reveals historical population genetic processes leading to contrasting phylogeographic patterns in co-distributed species

Coalescent samplers are computational time machines for inferring the historical demographic genetic processes that have given rise to observable patterns of spatial genetic variation among contemporary populations. We have used traditional characterizations of population structure and coalescent‐based inferences about demographic processes to reconstruct the population histories of two co‐distributed marine species, the frilled dog whelk, Nucella lamellosa, and the bat star, Patiria miniata. Analyses of population structure were consistent with previous work in both species except that additional samples of N. lamellosa showed a larger regional genetic break on Vancouver Island (VI) rather than between the southern Alexander Archipelago as in P. miniata. Our understanding of the causes, rather than just the patterns, of spatial genetic variation was dramatically improved by coalescent analyses that emphasized variation in population divergence times. Overall, gene flow was greater in bat stars (planktonic development) than snails (benthic development) but spatially homogeneous within species. In both species, these large phylogeographic breaks corresponded to relatively ancient divergence times between populations rather than regionally restricted gene flow. Although only N. lamellosa shows a large break on VI, population separation times on VI are congruent between species, suggesting a similar response to late Pleistocene ice sheet expansion. The absence of a phylogeographic break in P. miniata on VI can be attributed to greater gene flow and larger effective population size in this species. Such insights put the relative significance of gene flow into a more comprehensive historical biogeographic context and have important implications for conservation and landscape genetic studies that emphasize the role of contemporary gene flow and connectivity in shaping patterns of population differentiation.

It’s About Time: Divergence, Demography, and the Evolution of Developmental Modes in Marine Invertebrates

Differences in larval developmental mode are predicted to affect ecological and evolutionary processes ranging from gene flow and population bottlenecks to rates of population recovery from anthropogenic disturbance and capacity for local adaptation. The most powerful tests of these predictions use comparisons among species to ask how phylogeographic patterns are correlated with the evolution and loss of prolonged planktonic larval development. An important and largely untested assumption of these studies is that interspecific differences in population genetic structure are mainly caused by differences in dispersal and gene flow (rather than by differences in divergence times among populations or changes in effective population sizes), and that species with similar patterns of spatial genetic variation have similar underlying temporal demographic histories. Teasing apart these temporal and spatial patterns is important for understanding the causes and consequences of evolutionary changes in larval developmental mode. New analytical methods that use the coalescent history of allelic diversity can reveal these temporal patterns, test the strength of traditional population-genetic explanations for variation in spatial structure based on differences in dispersal, and identify strongly supported alternative explanations for spatial structure based on demographic history rather than on gene flow alone. We briefly review some of these recent analytical developments, and show their potential for refining ideas about the correspondence between the evolution of larval developmental mode, population demographic history, and spatial genetic variation.

A SIMPLE TECHNIQUE FOR PHYSICAL MARKING OF LARVAE OF MARINE BIVALVES

Abstract The identification of effective, nontoxic means for physically marking and tracking marine invertebrate larvae is a necessary step towards meeting a major goal of modern marine population biology, the direct measurement of larval dispersal. An inexpensive, rapid and effective means for marking bivalve larvae would be particularly useful because, as a taxonomic group, bivalves contain many commercially important and exploited species. Likewise, bivalves produce large numbers of propagules for experimental procedures and, for many species, methods for rearing larvae have been well established. Calcein has been used as a marker in numerous studies of adults and juveniles of calcium-carbonate-containing marine organisms, but its effects on small and sensitive life history stages such as embryos and larvae can be detrimental. We show that calcein can be used to rapidly and effectively mark large numbers of larvae from two bivalve species, Argopecten irradians concentricus (Say, 1822) (the Bay Scallop) and Mytilus trossulus Gould, 1850 (the Bay Mussel). Calcein had no detectable negative effects on growth or survivorship of larvae of either species; therefore, this fluorescent mark should serve as a useful tool for directly tracking dispersal of these species in the field. Our marking method is simple and inexpensive and can easily be used to determine the effectiveness and potential toxicity of the calcein mark for other bivalves.

Retrospective coalescent methods and the reconstruction of metapopulation histories in the sea

Phylogeographic analyses are a key interface between ecological and evolutionary ways of knowing because such analyses integrate the cumulative effects of demographic (ecological) processes over geological (evolutionary) time scales. Newly developed coalescent methods allow evolutionary ecologists to overcome some limitations associated with inferring population history from classic methods such as Wright’s FST. Here we briefly contrast classic and coalescent methods for looking backward in time through a population genetic lens, focusing on the key advantages of the isolation-with-migration (IM) class of coalescent methods for distinguishing ancient connectedness from actual recurrent contemporary gene flow as causes of genetic similarity or differentiation among populations. Making this critical distinction can lead to the discovery of otherwise obscured histories underlying conventional patterns of spatial variation. We illustrate the importance of these insights using analyses of Pacific fishes, snails, and sea stars in which population sizes and divergence times are more important than rates of contemporary gene flow as determinants of population genetic differentiation. We then extend the IM method to genetic data from two model metapopulation species (California abalone, Australian damselfish). The analyses show the potential use of non-equilibrium IM methods for differentiating among metapopulation models that make different predictions about population parameters and have different implications for the design of marine protected areas and other conservation goals. At face value, the results largely rule out classic metapopulation dynamics (dominated by extinction and colonization rather than connectivity via ongoing recurrent gene flow) but, at the same time, do not strongly support a modern marine metapopulation dynamic (ecologically significant connectivity between demes). However, the results also highlight the need for much more data (i.e., loci) sampled on different spatial scales in order to determine whether metapopulation dynamics might exist on smaller scales than are typically sampled by most phylogeographers and landscape geneticists.

Demographic Processes Underlying Subtle Patterns of Population Structure in the Scalloped Hammerhead Shark, Sphyrna lewini

Genetic diversity (θ), effective population size (Ne), and contemporary levels of gene flow are important parameters to estimate for species of conservation concern, such as the globally endangered scalloped hammerhead shark, Sphyrna lewini. Therefore, we have reconstructed the demographic history of S. lewini across its Eastern Pacific (EP) range by applying classical and coalescent population genetic methods to a combination of 15 microsatellite loci and mtDNA control region sequences. In addition to significant population genetic structure and isolation-by-distance among seven coastal sites between central Mexico and Ecuador, the analyses revealed that all populations have experienced a bottleneck and that all current values of θ are at least an order of magnitude smaller than ancestral θ, indicating large decreases in Ne (θ = 4Neμ), where μ is the mutation rate. Application of the isolation-with-migration (IM) model showed modest but significant genetic connectivity between most sampled sites (point estimates of Nm = 0.1–16.7), with divergence times (t) among all populations significantly greater than zero. Using a conservative (i.e., slow) fossil-based taxon-specific phylogenetic calibration for mtDNA mutation rates, posterior probability distributions (PPDs) for the onset of the decline in Ne predate modern fishing in this region. The cause of decline over the last several thousand years is unknown but is highly atypical as a post-glacial demographic history. Regardless of the cause, our data and analyses suggest that S. lewini was far more abundant throughout the EP in the past than at present.

New microsatellite loci for the endangered scalloped hammerhead shark, Sphyrna lewini.

We isolated 15 microsatellite markers for the scalloped hammerhead shark, Sphyrna lewini. Loci were tested on 80 specimens of S. lewini from four Eastern Pacific samples. The number of alleles per locus ranged from 6 to 31 (mean = 14). Observed and expected levels of heterozygosity per locus ranged from 0.39 to 0.91 (mean = 0.70) and from 0.54 to 0.90 (mean = 0.76), respectively. No pairs of loci were in gametic disequilibrium after Bonferroni correction of α. One locus showed significantly lower heterozygosity than expected under Hardy–Weinberg proportions in two populations, possibly caused by null alleles.

Selection and demographic history shape the molecular evolution of the gamete compatibility protein bindin in Pisaster sea stars

Reproductive compatibility proteins have been shown to evolve rapidly under positive selection leading to reproductive isolation, despite the potential homogenizing effects of gene flow. This process has been implicated in both primary divergence among conspecific populations and reinforcement during secondary contact; however, these two selective regimes can be difficult to discriminate from each other. Here, we describe the gene that encodes the gamete compatibility protein bindin for three sea star species in the genus Pisaster. First, we compare the full-length bindin-coding sequence among all three species and analyze the evolutionary relationships between the repetitive domains of the variable second bindin exon. The comparison suggests that concerted evolution of repetitive domains has an effect on bindin divergence among species and bindin variation within species. Second, we characterize population variation in the second bindin exon of two species: We show that positive selection acts on bindin variation in Pisaster ochraceus but not in Pisaster brevispinus, which is consistent with higher polyspermy risk in P. ochraceus. Third, we show that there is no significant genetic differentiation among populations and no apparent effect of sympatry with congeners that would suggest selection based on reinforcement. Fourth, we combine bindin and cytochrome c oxidase 1 data in isolation-with-migration models to estimate gene flow parameter values and explore the historical demographic context of our positive selection results. Our findings suggest that positive selection on bindin divergence among P. ochraceus alleles can be accounted for in part by relatively recent northward population expansions that may be coupled with the potential homogenizing effects of concerted evolution.