ASSORTATIVE FERTILIZATION AND SELECTION AT LARVAL STAGE IN THE MUSSELS MYTILUS EDULIS AND M. GALLOPROVINCIALIS

The divergence history of European blue mussel species reconstructed from Approximate Bayesian Computation: the effects of sequencing techniques and sampling strategies

Genome-scale diversity data are increasingly available in a variety of biological systems, and can be used to reconstruct the past evolutionary history of species divergence. However, extracting the full demographic information from these data is not trivial, and requires inferential methods that account for the diversity of coalescent histories throughout the genome. Here, we evaluate the potential and limitations of one such approach. We reexamine a well-known system of mussel sister species, using the joint site frequency spectrum (jSFS) of synonymous mutations computed either from exome capture or RNA-seq, in an Approximate Bayesian Computation (ABC) framework. We first assess the best sampling strategy (number of: individuals, loci, and bins in the jSFS), and show that model selection is robust to variation in the number of individuals and loci. In contrast, different binning choices when summarizing the jSFS, strongly affect the results: including classes of low and high frequency shared polymorphisms can more effectively reveal recent migration events. We then take advantage of the flexibility of ABC to compare more realistic models of speciation, including variation in migration rates through time (i.e., periodic connectivity) and across genes (i.e., genome-wide heterogeneity in migration rates). We show that these models were consistently selected as the most probable, suggesting that mussels have experienced a complex history of gene flow during divergence and that the species boundary is semi-permeable. Our work provides a comprehensive evaluation of ABC demographic inference in mussels based on the coding jSFS, and supplies guidelines for employing different sequencing techniques and sampling strategies. We emphasize, perhaps surprisingly, that inferences are less limited by the volume of data, than by the way in which they are analyzed.

SeaBase: a multispecies transcriptomic resource and platform for gene network inference

Marine and aquatic animals are extraordinarily useful as models for identifying mechanisms of development and evolution, regeneration, resistance to cancer, longevity and symbiosis, among many other areas of research. This is due to the great diversity of these organisms and their wide-ranging capabilities. Genomics tools are essential for taking advantage of these "free lessons" of nature. However, genomics and transcriptomics are challenging in emerging model systems. Here, we present SeaBase, a tool for helping to meet these needs. Specifically, SeaBase provides a platform for sharing and searching transcriptome data. More importantly, SeaBase will support a growing number of tools for inferring gene network mechanisms. The first dataset available on SeaBase is a developmental transcriptomic profile of the sea anemone Nematostella vectensis (Anthozoa, Cnidaria). Additional datasets are currently being prepared and we are aiming to expand SeaBase to include user-supplied data for any number of marine and aquatic organisms, thereby supporting many potentially new models for gene network studies. SeaBase can be accessed online at: http://seabase.core.cli.mbl.edu.

Adaptive gamete-recognition divergence in a hybridizing Mytilus population

Gamete-recognition proteins often evolve rapidly, but it is not known if their divergence occurs within species and corresponds with the evolution of reproductive isolation, or if divergence typically accumulates between already isolated lineages. We examined the evolution of a candidate gamete-recognition protein in several sympatric and allopatric populations of Mytilus blue mussels, species that hybridize in nature. Within a single species, Mytilus galloprovincialis, we found adaptive divergence of Lysin-M7, a sperm acrosomal protein that dissolves the egg vitelline envelope during fertilization. Mytilus galloprovincialis Lysin-M7 alleles group into two distinct clades (termed G and G(D)), and individual alleles in these clades are separated from each other by at least three and up to eleven amino-acid substitutions. Maximum-likelihood estimates of selective pressure (dN/dS =omega) implicate selection in the divergence between M. galloprovincialis Lysin-M7 clades, and within the G(D) clade. Exact tests of population differentiation indicate that the relative frequency of G and G(D) Lysin-M7 alleles differs significantly among M. galloprovincialis populations. Compared with allopatric Mediterranean samples, Lysin-M7 alleles in the G(D) clade are found at elevated frequency in samples from the East Atlantic and California, areas of secondary contact and hybridization between Mytilus species, and Australia, an area of unknown species composition. Adaptive divergence between the alleles most common in allopatry and those found at elevated frequency in samples from sympatry suggests that selection pressures acting in hybridizing populations, likely following Pleistocene secondary contact with M. edulis in the East Atlantic, drove the divergence of Lysin-M7 in M. galloprovincialis.

Patterns of larval dispersal and their effect on the maintenance of a blue mussel hybrid zone in southwestern England

The blue mussels Mytilus edulis and M. galloprovincialis hybridize in southwestern England. Within this hybrid zone environmentally based directional selection favors individuals with alleles specific to M. galloprovincialis. What forces are countering this directional selection and allowing for the maintenance of a stable hybrid population are unknown. We used both the genetics of recently settled larvae and a fine-scale model of the physical oceanography of the region to determine the patterns of larval dispersal throughout the hybrid zone and the bordering parental populations. Evidence from both the model and the genetics suggests that the hybrid zone lies between two barriers to dispersal. Start Point separates the M. edulis population from the hybrid zone and allows minimal dispersal from the hybrid zone into the M. edulis population, but none in the other direction. Likewise, the M. galloprovincialis populations along the northern coast of Cornwall regularly receive immigrating larvae from the hybrid zone, but larvae from the M. galloprovincialis population do not enter the hybrid zone. However, larvae settling at hybrid zone sites have high frequencies of alleles specific to M. edulis, suggesting that reproductive barriers, selection in the larval stage, or gene flow from an undetermined source is effectively balancing the directional selection observed in the adults.