Assortative fertilization and selection at larval stage in the mussels Mytilus edulis and M. galloprovincialis

Patterns of admixture and introgression in a mosaic Mytilus galloprovincialis and Mytilus edulis hybrid zone in SW England

The study of hybrid zones offers important insights into speciation. Earlier studies on hybrid populations of the marine mussel species Mytilus edulis and Mytilus galloprovincialis in SW England provided evidence of admixture but were constrained by the limited number of molecular markers available. We use 57 ancestry-informative SNPs, most of which have been mapped genetically, to provide evidence of distinctive differences between admixed populations in SW England and asymmetrical introgression from M. edulis to M. galloprovincialis. We combine the genetic study with analysis of phenotypic traits of potential ecological and adaptive significance. We demonstrate that hybrid individuals have brown mantle edges unlike the white or purple in the parental species, suggesting allelic or non-allelic genomic interactions. We report differences in gonad development stage between the species consistent with a prezygotic barrier between the species. By incorporating results from publications dating back to 1980, we confirm the long-term stability of the hybrid zone despite higher viability of M. galloprovincialis. This stability coincides with a dramatic change in temperature of UK coastal waters and suggests that these hybrid populations might be resisting the effects of global warming. However, a single SNP locus associated with the Notch transmembrane signalling protein shows a markedly different pattern of variation to the others and might be associated with adaptation of M. galloprovincialis to colder northern temperatures.

Pre-introduction introgression contributes to parallel differentiation and contrasting hybridization outcomes between invasive and native marine mussels

Non-native species experience novel selection pressures in introduced environments and may interbreed with native lineages. Species introductions therefore provide opportunities to investigate repeated patterns of adaptation and introgression across replicated contact zones. Here, we investigate genetic parallelism between multiple introduced populations of the invasive marine mussel, Mytilus galloprovincialis, in the absence (South Africa and California) and presence of hybridization with a native congener (Mytilus planulatus in Batemans Bay and Sydney Harbour, Australia). Repeatability in post-introduction differentiation from native-range populations varied between genetically distinct Atlantic and Mediterranean lineages, with Atlantic-derived introductions displaying high differentiation (maxF_ST  > 0.4) and parallelism at outlier loci. Identification of long noncoding RNA transcripts (lncRNA) additionally allowed us to clarify that parallel responses are largely limited to protein-coding loci, with lncRNAs likely evolving under evolutionary constraints. Comparisons of independent hybrid zones revealed differential introgression most strongly in Batemans Bay, with an excess of M. galloprovincialis ancestry and resistance to introgression at loci differentiating parental lineages (M. planulatus and Atlantic M. galloprovincialis). Additionally, contigs putatively introgressed with divergent alleles from a closely related species, Mytilus edulis, showed stronger introgression asymmetries compared with genome-wide trends and also diverged in parallel in both Atlantic-derived introductions. These results suggest that divergent demographic histories experienced by introduced lineages, including pre-introduction introgression, influence contemporary admixture dynamics. Our findings build on previous investigations reporting contributions of historical introgression to intrinsic reproductive architectures shared between marine lineages and illustrate that interspecific introgression history can shape differentiation between colonizing populations and their hybridization with native congeners.

Massive gene presence-absence variation shapes an open pan-genome in the Mediterranean mussel

BACKGROUND

The Mediterranean mussel Mytilus galloprovincialis is an ecologically and economically relevant edible marine bivalve, highly invasive and resilient to biotic and abiotic stressors causing recurrent massive mortalities in other bivalves. Although these traits have been recently linked with the maintenance of a high genetic variation within natural populations, the factors underlying the evolutionary success of this species remain unclear.

RESULTS

Here, after the assembly of a 1.28-Gb reference genome and the resequencing of 14 individuals from two independent populations, we reveal a complex pan-genomic architecture in M. galloprovincialis, with a core set of 45,000 genes plus a strikingly high number of dispensable genes (20,000) subject to presence-absence variation, which may be entirely missing in several individuals. We show that dispensable genes are associated with hemizygous genomic regions affected by structural variants, which overall account for nearly 580 Mb of DNA sequence not included in the reference genome assembly. As such, this is the first study to report the widespread occurrence of gene presence-absence variation at a whole-genome scale in the animal kingdom.

CONCLUSIONS

Dispensable genes usually belong to young and recently expanded gene families enriched in survival functions, which might be the key to explain the resilience and invasiveness of this species. This unique pan-genome architecture is characterized by dispensable genes in accessory genomic regions that exceed by orders of magnitude those observed in other metazoans, including humans, and closely mirror the open pan-genomes found in prokaryotes and in a few non-metazoan eukaryotes.

How do species barriers decay? Concordance and local introgression in mosaic hybrid zones of mussels

The Mytilus complex of marine mussel species forms a mosaic of hybrid zones, found across temperate regions of the globe. This allows us to study 'replicated' instances of secondary contact between closely related species. Previous work on this complex has shown that local introgression is both widespread and highly heterogeneous, and has identified SNPs that are outliers of differentiation between lineages. Here, we developed an ancestry-informative panel of such SNPs. We then compared their frequencies in newly sampled populations, including samples from within the hybrid zones, and parental populations at different distances from the contact. Results show that close to the hybrid zones, some outlier loci are near to fixation for the heterospecific allele, suggesting enhanced local introgression, or the local sweep of a shared ancestral allele. Conversely, genomic cline analyses, treating local parental populations as the reference, reveal a globally high concordance among loci, albeit with a few signals of asymmetric introgression. Enhanced local introgression at specific loci is consistent with the early transfer of adaptive variants after contact, possibly including asymmetric bi-stable variants (Dobzhansky-Muller incompatibilities), or haplotypes loaded with fewer deleterious mutations. Having escaped one barrier, however, these variants can be trapped or delayed at the next barrier, confining the introgression locally. These results shed light on the decay of species barriers during phases of contact.

Replicated anthropogenic hybridisations reveal parallel patterns of admixture in marine mussels

Human-mediated transport creates secondary contacts between genetically differentiated lineages, bringing new opportunities for gene exchange. When similar introductions occur in different places, they provide informally replicated experiments for studying hybridisation. We here examined 4,279 Mytilus mussels, sampled in Europe and genotyped with 77 ancestry-informative markers. We identified a type of introduced mussels, called "dock mussels," associated with port habitats and displaying a particular genetic signal of admixture between M. edulis and the Mediterranean lineage of M. galloprovincialis. These mussels exhibit similarities in their ancestry compositions, regardless of the local native genetic backgrounds and the distance separating colonised ports. We observed fine-scale genetic shifts at the port entrance, at scales below natural dispersal distance. Such sharp clines do not fit with migration-selection tension zone models, and instead suggest habitat choice and early-stage adaptation to the port environment, possibly coupled with connectivity barriers. Variations in the spread and admixture patterns of dock mussels seem to be influenced by the local native genetic backgrounds encountered. We next examined departures from the average admixture rate at different loci, and compared human-mediated admixture events, to naturally admixed populations and experimental crosses. When the same M. galloprovincialis background was involved, positive correlations in the departures of loci across locations were found; but when different backgrounds were involved, no or negative correlations were observed. While some observed positive correlations might be best explained by a shared history and saltatory colonisation, others are likely produced by parallel selective events. Altogether, genome-wide effect of admixture seems repeatable and more dependent on genetic background than environmental context. Our results pave the way towards further genomic analyses of admixture, and monitoring of the spread of dock mussels both at large and at fine spacial scales.

Comparison between single and multi-locus approaches for specimen identification in Mytilus mussels

Mytilus mussels have been the object of much research given their sentinel role in coastal ecosystems and significant value as an aquaculture resource appreciated for both, its flavour and nutritional content. Some of the most-studied Mytilus species are M. edulis, M. galloprovincialis, M. chilensis and M. trossulus. As species identification based on morphological characteristics of Mytilus specimens is difficult, molecular markers are often used. Single-locus markers can give conflicting results when used independently; not all markers differentiate among all species, and the markers target genomic regions with different evolutionary histories. We evaluated the concordance between the PCR-RFLP markers most commonly-used for species identification in mussels within the Mytilus genus (Me15-16, ITS, mac-1, 16S rRNA and COI) when used alone (mono-locus approach) or together (multi-locus approach). In this study, multi-locus strategy outperformed the mono-locus methods, clearly identifying all four species and also showed similar specimen identification performance than a 49 SNPs panel. We hope that these findings will contribute to a better understanding of DNA marker-based analysis of Mytilus taxa. These results support the use of a multi-locus approach when studying this important marine resource, including research on food quality and safety, sustainable production and conservation.

The hidden side of a major marine biogeographic boundary: a wide mosaic hybrid zone at the Atlantic-Mediterranean divide reveals the complex interaction between natural and genetic barriers in mussels

The Almeria-Oran Front (AOF) is a recognised hotspot of genetic differentiation in the sea, with genetic discontinuities reported in more than 50 species. The AOF is a barrier to dispersal and an ecological boundary; both can determine the position of these genetic breaks. However, the maintenance of genetic differentiation is likely reinforced by genetic barriers. A general drawback of previous studies is an insufficient density of sampling sites at the transition zone, with a conspicuous lack of samples from the southern coastline. We analysed the fine-scale genetic structure in the mussel Mytilus galloprovincialis using a few ancestry-informative loci previously identified from genome scans. We discovered a 600-km-wide mosaic hybrid zone eastward of the AOF along the Algerian coasts. This mosaic zone provides a new twist to our understanding of the Atlantic-Mediterranean transition because it demonstrates that the two lineages can live in sympatry with ample opportunities to interbreed in a large area, but they hardly do so. This implies that some form of reproductive isolation must exist to maintain the two genetic backgrounds locally cohesive. The mosaic zone ends with an abrupt genetic shift at a barrier to dispersal in the Gulf of Bejaia, Eastern Algeria. Simulations of endogenous or exogenous selection in models that account for the geography and hydrodynamic features of the region support the hypothesis that sister hybrid zones could have been differentially trapped at two alternative barriers to dispersal and/or environmental boundaries, at Almeria in the north and Bejaia in the south. A preponderantly unidirectional north-south gene flow next to the AOF can also maintain a patch of intrinsically maintained genetic background in the south and the mosaic structure, even in the absence of local adaptation. Our results concur with the coupling hypothesis that suggests that natural barriers can explain the position of genetic breaks, while their maintenance depends on genetic barriers.

Connectivity of larval stages of sedentary marine communities between hard substrates and offshore structures in the North Sea

Man-made structures including rigs, pipelines, cables, renewable energy devices, and ship wrecks, offer hard substrate in the largely soft-sediment environment of the North Sea. These structures become colonised by sedentary organisms and non-migratory reef fish, and form local ecosystems that attract larger predators including seals, birds, and fish. It is possible that these structures form a system of interconnected reef environments through the planktonic dispersal of the pelagic stages of organisms by ocean currents. Changes to the overall arrangement of hard substrate areas through removal or addition of individual man-made structures will affect the interconnectivity and could impact on the ecosystem. Here, we assessed the connectivity of sectors with oil and gas structures, wind farms, wrecks, and natural hard substrate, using a model that simulates the drift of planktonic stages of seven organisms with sedentary adult stages associated with hard substrate, applied to the period 2001-2010. Connectivity was assessed using a classification system designed to address the function of sectors in the network. Results showed a relatively stable overall spatial distribution of sector function but with distinct variations between species and years. The results are discussed in the context of decommissioning of oil and gas infrastructure in the North Sea.

Coadapted genomes and selection on hybrids: Fisher's geometric model explains a variety of empirical patterns

Natural selection plays a variety of roles in hybridization, speciation, and admixture. Most research has focused on two extreme cases: crosses between closely related inbred lines, where hybrids are fitter than their parents, or crosses between effectively isolated species, where hybrids suffer severe breakdown. But many natural populations must fall into intermediate regimes, with multiple types of gene interaction, and these are more difficult to study. Here, we develop a simple fitness landscape model, and show that it naturally interpolates between previous modeling approaches, which were designed for the extreme cases, and invoke either mildly deleterious recessives, or discrete hybrid incompatibilities. Our model yields several new predictions, which we test with genomic data from Mytilus mussels, and published data from plants (Zea, Populus, and Senecio) and animals (Mus, Teleogryllus, and Drosophila). The predictions are generally supported, and the model explains a number of surprising empirical patterns. Our approach enables novel and complementary uses of genome-wide datasets, which do not depend on identifying outlier loci, or "speciation genes" with anomalous effects. Given its simplicity and flexibility, and its predictive successes with a wide range of data, the approach should be readily extendable to other outstanding questions in the study of hybridization.

RNA-seq coupled to proteomic analysis reveals high sperm proteome variation between two closely related marine mussel species

Speciation mechanisms in marine organisms have attracted great interest because of the apparent lack of substantial barriers to genetic exchange in marine ecosystems. Marine mussels of the Mytilus edulis species complex provide a good model to study mechanisms underlying species formation. They hybridise extensively at many localities and both pre- and postzygotic isolating mechanisms may be operating. Mussels have external fertilisation and sperm cells should show specific adaptations for survival and successful fertilisation. Sperm thus represent key targets in investigations of the molecular mechanisms underlying reproductive isolation. We undertook a deep transcriptome sequencing (RNA-seq) of mature male gonads and a 2DE/MS-based proteome analysis of sperm from Mytilus edulis and M. galloprovincialis raised in a common environment. We provide evidence of extensive expression differences between the two mussel species, and general agreement between the transcriptomic and proteomic results in the direction of expression differences between species. Differential expression is marked for mitochondrial genes and for those involved in spermatogenesis, sperm motility, sperm-egg interactions, the acrosome reaction, sperm capacitation, ATP reserves and ROS production. Proteins and their corresponding genes might thus be good targets in further genomic analysis of reproductive barriers between these closely related species. SIGNIFICANCE: Model systems for the study of fertilization include marine invertebrates with external fertilisation, such as abalones, sea urchins and mussels, because of the ease with which large quantities of gametes released into seawater can be collected after induced spawning. Unlike abalones and sea urchins, hybridisation has been reported between mussels of different Mytilus spp., which thus makes them very appealing for the study of reproductive isolation at both pre- and postzygotic levels. There is a lack of empirical proteomic studies on sperm samples comparing different Mytilus species, which could help to advance this study. A comparative analysis of sperm proteomes across different taxa may provide important insights into the fundamental molecular processes and mechanisms involved in reproductive isolation. It might also contribute to a better understanding of sperm function and of the adaptive evolution of sperm proteins in different taxa. There is now growing evidence from genomics studies that multiple protein complexes and many individual proteins might have important functions in sperm biology and the fertilisation process. From an applied perspective, the identification of sperm-specific proteins could also contribute to the improved understanding of fertility problems and as targets for fertility control.

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.

Adaptive responses along a depth and a latitudinal gradient in the endemic seagrass Posidonia oceanica

Seagrass meadows provide important ecosystem services and are critical for the survival of the associated invertebrate community. However, they are threatened worldwide by human-driven environmental change. Understanding the seagrasses' potential for adaptation is critical to assess not only their ability to persist under future global change scenarios, but also to assess the persistence of the associated communities. Here we screened a wild population of Posidonia oceanica, an endemic long-lived seagrass in the Mediterranean Sea, for genes that may be target of environmental selection, using an outlier and a genome-wide transcriptome analysis. We identified loci where polymorphism or differential expression was associated with either a latitudinal or a bathymetric gradient, as well as with both gradients in an effort to identify loci associated with temperature and light. We found the candidate genes underlying growth and immunity to be divergent between populations adapted to different latitudes and/or depths, providing evidence for local adaptation. Furthermore, we found evidence of reduced gene flow among populations including adjacent populations. Reduced gene flow, combined with low sexual recombination, small effective population size, and long generation time of P. oceanica raises concerns for the long-term persistence of this species, especially in the face of rapid environmental change driven by human activities.

Do Sperm Really Compete and Do Eggs Ever Have a Choice? Adult Distribution and Gamete Mixing Influence Sexual Selection, Sexual Conflict, and the Evolution of Gamete Recognition Proteins in the Sea

The evolution of gametic compatibility and the effectiveness of compatibility, within and across species, depend on whether sperm from different males directly compete for an egg and whether eggs ever have a choice. Direct sperm competition and egg choice depend on whether sperm from different males arrive at an egg in the brief interval between first sperm contact and fertilization. Although this process may be relevant for all sexually reproducing organisms, it is most easily examined in aquatic external fertilizers. When sperm are released into the sea, packets of seawater at the spatial scale relevant to single eggs might contain sperm from only one male, eliminating the potential for direct sperm competition and egg choice. Field experiments and a simple heuristic model examining the degree of sperm mixing for the sea urchin Strongylocentrotus franciscanus indicate that degree of competitive fertilization depends on density and distribution of competing males and that the nature of this competition influences whether males with high- or low-affinity gamete recognition protein genotypes have higher reproductive success. These results provide a potential explanation for the generation and maintenance of variation in gamete recognition proteins and why effectiveness of conspecific sperm precedence can be density dependent.

Metapopulation patterns of additive and nonadditive genetic variance in the sea bass (Dicentrarchus labrax)

Describing and explaining the geographic within‐species variation in phenotypes (“phenogeography”) in the sea over a species distribution range is central to our understanding of a variety of eco‐evolutionary topics. However, phenogeographic studies that have a large potential to investigate adaptive variation are overcome by phylogeographic studies, still mainly focusing on neutral markers. How genotypic and phenotypic data could covary over large geographic scales remains poorly understood in marine species. We crossed 75 noninbred sires (five origins) and 26 dams (two origins; each side of a hybrid zone) in a factorial diallel cross in order to investigate geographic variation for early survival and sex ratio in the metapopulation of the European sea bass (Dicentrarchus labrax), a highly prized marine fish species. Full‐sib families (N = 1,950) were produced and reared in a common environment. Parentage assignment of 7,200 individuals was performed with seven microsatellite markers. Generalized linear models showed significant additive effects for both traits and pleiotropy between traits. A significant nonadditive genetic effect was detected. Different expression of traits and distinct relative performances were found for reciprocal crosses involving populations located on each side of the main hybrid zone located at the Almeria‐Oran front, illustrating asymmetric reproductive isolation. The poor fitness performance observed for the Western Mediterranean population of sea bass is discussed as it represents the main source of seed hatchery production, but also because it potentially illustrates nonadaptive introgression and maladaptation.

A First Insight into the Genome of the Filter-Feeder Mussel Mytilus galloprovincialis

Mussels belong to the phylum Mollusca, one of the largest and most diverse taxa in the animal kingdom. Despite their importance in aquaculture and in biology in general, genomic resources from mussels are still scarce. To broaden and increase the genomic knowledge in this family, we carried out a whole-genome sequencing study of the cosmopolitan Mediterranean mussel (Mytilus galloprovincialis). We sequenced its genome (32X depth of coverage) on the Illumina platform using three pair-end libraries with different insert sizes. The large number of contigs obtained pointed out a highly complex genome of 1.6 Gb where repeated elements seem to be widespread (~30% of the genome), a feature that is also shared with other marine molluscs. Notwithstanding the limitations of our genome sequencing, we were able to reconstruct two mitochondrial genomes and predict 10,891 putative genes. A comparative analysis with other molluscs revealed a gene enrichment of gene ontology categories related to multixenobiotic resistance, glutamate biosynthetic process, and the maintenance of ciliary structures.

Low reproductive isolation and highly variable levels of gene flow reveal limited progress towards speciation between European river and brook lampreys

Ecologically based divergent selection is a factor that could drive reproductive isolation even in the presence of gene flow. Population pairs arrayed along a continuum of divergence provide a good opportunity to address this issue. Here, we used a combination of mating trials, experimental crosses and population genetic analyses to investigate the evolution of reproductive isolation between two closely related species of lampreys with distinct life histories. We used microsatellite markers to genotype over 1000 individuals of the migratory parasitic river lamprey (Lampetra fluviatilis) and freshwater-resident nonparasitic brook lamprey (Lampetra planeri) distributed in 10 sympatric and parapatric population pairs in France. Mating trials, parentage analyses and artificial fertilizations demonstrated a low level of reproductive isolation between species even though size-assortative mating may contribute to isolation. Most parapatric population pairs were strongly differentiated due to the joint effects of geographic distance and barriers to migration. In contrast, we found variable levels of gene flow between sympatric populations ranging from panmixia to moderate differentiation, which indicates a gradient of divergence with some population pairs that may correspond to alternative morphs or ecotypes of a single species and others that remain partially isolated. Ecologically based divergent selection may explain these variable levels of divergence among sympatric population pairs, but incomplete genome swamping following secondary contact could have also played a role. Overall, this study illustrates how highly differentiated phenotypes can be maintained despite high levels of gene flow that limit the progress towards speciation.

Does natural selection explain the fine scale genetic structure at the nuclear exon Glu-5' in blue mussels from Kerguelen?

The Kerguelen archipelago, isolated in the Southern Ocean, shelters a blue mussel Mytilus metapopulation far from any influence of continental populations or any known hybrid zone. The finely carved coast leads to a highly heterogeneous habitat. We investigated the impact of the environment on the genetic structure in those Kerguelen blue mussels by relating allele frequencies to habitat descriptors. A total sample comprising up to 2248 individuals from 35 locations was characterized using two nuclear markers, mac-1 and Glu-5', and a mitochondrial marker (COI). The frequency data from 9 allozyme loci in 9 of these locations were also reanalyzed. Two other nuclear markers (EFbis and EFprem's) were monomorphic. Compared to Northern Hemisphere populations, polymorphism in Kerguelen blue mussels was lower for all markers except for the exon Glu-5'. At Glu-5', genetic differences were observed between samples from distinct regions (F CT = 0.077), as well as within two regions, including between samples separated by <500 m. No significant differentiation was observed in the AMOVA analyses at the two other markers (mac-1 and COI). Like mac-1, all allozyme loci genotyped in a previous publication, displayed lower differentiation (Jost's D) and F ST values than Glu-5'. Power simulations and confidence intervals support that Glu-5' displays significantly higher differentiation than the other loci (except a single allozyme for which confidence intervals overlap). AMOVA analyses revealed significant effects of the giant kelp Macrocystis and wave exposure on this marker. We discuss the influence of hydrological conditions on the genetic differentiation among regions. In marine organisms with high fecundity and high dispersal potential, gene flow tends to erase differentiation, but this study showed significant differentiation at very small distance. This may be explained by the particular hydrology and the carved coastline of the Kerguelen archipelago, together with spatially variable selection at Glu-5'.

Conspecific sperm precedence is a reproductive barrier between free-spawning marine mussels in the northwest atlantic mytilus hybrid zone

Reproductive isolation at the gamete stage has become a focus of speciation research because of its potential to evolve rapidly between closely related species. Conspecific sperm precedence (CSP), a type of gametic isolation, has been demonstrated in a number of taxa, both marine and terrestrial, with the potential to play an important role in speciation. Free-spawning marine invertebrates are ideal subjects for the study of CSP because of a likely central role for gametic barriers in reproductive isolation. The western Atlantic Mytilus blue mussel hybrid zone, ranging from the Atlantic Canada to eastern Maine, exhibits characteristics conducive to the study of CSP. Previous studies have shown that gametic incompatibility is incomplete, variable in strength and the genotype distribution is bimodal-dominated by the parental species, with a low frequency of hybrids. We conducted gamete crossing experiments using M. trossulus and M. edulis individuals collected from natural populations during the spring spawning season in order to detect the presence or absence of CSP within this hybrid zone. We detected CSP, defined here as a reduction in heterospecific offspring from competitive fertilizations in vitro compared to that seen in non-competitive fertilizations, in five of the twelve crosses in which conspecific crosses were detectable. This is the first finding of CSP in a naturally hybridizing population of a free-spawning marine invertebrate. Our findings support earlier predictions that CSP can promote assortative fertilization in bimodal hybrid zones, further advancing their hypothesized progression towards full speciation. Despite strong CSP numerous heterospecific fertilizations remain, reinforcing the hypothesis that compatible females are a source of hybrid offspring in mixed natural spawns.

Gene-flow in a mosaic hybrid zone: is local introgression adaptive?

Genome-wide scans of genetic differentiation between hybridizing taxa can identify genome regions with unusual rates of introgression. Regions of high differentiation might represent barriers to gene flow, while regions of low differentiation might indicate adaptive introgression-the spread of selectively beneficial alleles between reproductively isolated genetic backgrounds. Here we conduct a scan for unusual patterns of differentiation in a mosaic hybrid zone between two mussel species, Mytilus edulis and M. galloprovincialis. One outlying locus, mac-1, showed a characteristic footprint of local introgression, with abnormally high frequency of edulis-derived alleles in a patch of M. galloprovincialis enclosed within the mosaic zone, but low frequencies outside of the zone. Further analysis of DNA sequences showed that almost all of the edulis allelic diversity had introgressed into the M. galloprovincialis background in this patch. We then used a variety of approaches to test the hypothesis that there had been adaptive introgression at mac-1. Simulations and model fitting with maximum-likelihood and approximate Bayesian computation approaches suggested that adaptive introgression could generate a "soft sweep," which was qualitatively consistent with our data. Although the migration rate required was high, it was compatible with the functioning of an effective barrier to gene flow as revealed by demographic inferences. As such, adaptive introgression could explain both the reduced intraspecific differentiation around mac-1 and the high diversity of introgressed alleles, although a localized change in barrier strength may also be invoked. Together, our results emphasize the need to account for the complex history of secondary contacts in interpreting outlier loci.

Asymmetric conspecific sperm precedence in relation to spawning times in the Montastraea annularis species complex (Cnidaria: Scleractinia)

In broadcast spawners, prezygotic reproductive isolation depends on differences in the spatial and temporal patterns of gamete release and gametic incompatibility. Typically, gametic incompatibility is measured in no-choice crosses, but conspecific sperm precedence (CSP) can prevent hybridization in gametes that are compatible in the absence of sperm competition. Broadcast spawning corals in the Montastraea annularis species complex spawn annually on the same few evenings. Montastraea franksi spawns an average of 110 min before M. annularis, with a minimum gap of approximately 40 min. Gametes are compatible in no-choice heterospecific assays, but it is unknown whether eggs exhibit choice when in competition. Hybridization depends on either M. franksi eggs remaining unfertilized and in proximity to M. annularis when the latter species spawns or M. franksi sperm remaining in sufficient viable concentrations when M. annularis spawns. We found that the eggs of the early spawning M. franksi demonstrate strong CSP, whereas CSP appears to be lacking for M. annularis eggs. This study provides evidence of diverging gamete affinities between these recently separated species and suggests for the first time that selection may favour CSP in earlier spawning species when conspecific sperm is diluted and aged and is otherwise at a numeric and viability disadvantage with heterospecific sperm.

Assessing the potential for egg chemoattractants to mediate sexual selection in a broadcast spawning marine invertebrate

In numerous species, egg chemoattractants play a critical role in guiding sperm towards unfertilized eggs (sperm chemotaxis). Until now, the known functions of sperm chemotaxis include increasing the effective target size of eggs, thereby promoting sperm-egg encounters, and facilitating species recognition. Here, we report that in the broadcast spawning mussel, Mytilus galloprovincialis, egg chemoattractants may play an unforeseen role in sexual selection by enabling sperm to effectively 'choose' between the eggs of different conspecific females. In an initial experiment, we confirmed that sperm chemotaxis occurs in M. galloprovincialis by showing that sperm are attracted towards unfertilized eggs when given the choice of eggs or no eggs in a dichotomous chamber. We then conducted two cross-classified mating experiments, each comprising the same individual males and females crossed in identical male × female combinations, but under experimental conditions that offered sperm 'no-choice' (each fertilization trial took place in a Petri dish and involved a single male and female) or a 'choice' of a female's eggs (sperm were placed in the centre of a dichotomous choice chamber and allowed to choose eggs from different females). We show that male-by-female interactions characterized fertilization rates in both experiments, and that there was remarkable consistency between patterns of sperm migration in the egg-choice experiment and fertilization rates in the no-choice experiment. Thus, sperm appear to exploit chemical cues to preferentially swim towards eggs with which they are most compatible during direct sperm-to-egg encounters. These results reveal that sperm differentially select eggs on the basis of chemical cues, thus exposing the potential for egg chemoattractants to mediate mate choice for genetically compatible partners. Given the prevalence of sperm chemotaxis across diverse taxa, our findings may have broad implications for sexual selection in other mating systems.

Differential fertilization success between two populations of eastern oyster, Crassostrea virginica

Identification of mechanisms promoting prezygotic reproductive isolation and their prevalence are key goals in evolutionary biology because of their potential role in speciation. In marine broadcast-spawning species, molecular interactions between gamete surface proteins are more important than mating behavior for determining reproductive compatibility. Evidence for differential fertilization capacity has been reported from experiments utilizing competing sperm from two males sampled within populations and between species, but to our knowledge conspecific populations that might have diverged in allopatry have never been tested on the basis of sperm competition. In the present study, the gametic compatibility and embryo survivorship from matings between two allopatric populations of Crassostrea virginica, the eastern oyster, on either side of a genetic step cline were investigated. Fertilization success, embryo survival, and paternity data all indicated an absence of strong reproductive barriers between the two oyster populations, implicating other mechanisms for maintenance of the cline step. Sperm from northern male oysters showed a tendency to produce more larvae than expected when competing with sperm from southern male oysters. Although the northern male advantage was not strong, the trend implies that long-distance dispersal across the step cline might more successfully result in north-to-south gene flow than the reverse, providing a mechanistic hypothesis explaining the asymmetric cline shape.

Incomplete reproductive isolation in the blue mussel (Mytilus edulis and M. trossulus) hybrid zone in the Northwest Atlantic: role of gamete interactions and larval viability

Reproductive isolation in free-spawning organisms may involve only small changes in the gamete surface molecules that control fertilization, linking gamete incompatibility and speciation. Most studies have focused on species in which natural hybrids are absent and reproductive isolation is complete, but how gamete incompatibility evolves remains unclear. Reproductive isolation is incomplete between two sympatric mussel species (Mytilus edulis, M. trossulus) that hybridize in nature. In this study prezygotic and postzygotic components of reproductive incompatibility were examined in laboratory crosses. Conspecific crosses showed significantly greater rates of fertilization than heterospecific crosses, although some females of both species showed heterospecific gamete compatibility. The proportion of fertilized eggs developing into normal larvae was not significantly different between conspecific and heterospecific crosses, but survival of normal larvae was greater for conspecific crosses. Mixed-species sperm experiments suggested that conspecific sperm preference may further limit hybridization. The different components of reproductive incompatibility and total incompatibility varied among females of both species. Although our study has shown that partial reproductive isolation between M. edulis and M. trossulus involves both prezygotic gamete interactions and postzygotic larval survival, further research is required to determine the potential role of gamete incompatibility in the evolution of complete reproductive isolation between these species.

Assortative mating and spatial structure in hybrid zones

The spatial genetic composition of hybrid zones exhibits a range of possible patterns, with many characterized by patchy distributions. While several hypothetical explanations exist for the maintenance of these "mosaic" hybrid zones, they remain virtually unexplored theoretically. Using computer simulations we investigate the roles of dispersal and assortative mating in the formation and persistence of hybrid zone structure. To quantify mosaic structure we develop a likelihood method, which we apply to simulation and empirical data. We find that long distance dispersal can lead to a patchy distribution that assortative mating can then reinforce, ultimately producing a mosaic capable of persisting over evolutionarily significant periods of time. By reducing the mating success of rare males, assortative mating creates a positive within-patch frequency-dependent selective pressure. Selection against heterozygotes can similarly create a rare-type disadvantage and we show that it can also preserve structure. We find that mosaic structure is maintained across a range of assumptions regarding the form and strength of assortative mating. Interestingly, we find that higher levels of mosaic structure are sometimes observed for intermediate assortment strengths. The high incidence of assortment documented in hybrid zones suggests that it may play a key role in stabilizing their form and structure.

The flow of antimicrobial peptide genes through a genetic barrier between Mytilus edulis and M. galloprovincialis

We studied the population genetics of two antimicrobial peptide (AMP) loci, called Mytilin B and Mytilus galloprovincialis defensin 2 (MGD2), in the secondary contact mosaic hybrid zone between Mytilus edulis and M. galloprovincialis. The isolation period between the two species was estimated to be approximately 1 million years (range, 0.5 million to 2 million years) long. During this period, coevolution between microbes and the immune system has likely occurred. The secondary contact, which would date back to approximately 25,000 (0-200,000) years, recently allowed these coadaptations to be rearranged through hybridization. Distinctive polymorphisms were uncovered in coding sequences of the two AMP loci such as insertion/deletion of codons or bisubstituted codons. Very low levels of differentiation were observed between populations of the two species at both loci, while other nuclear loci often showed marked structure among the same samples. The absence of population differentiation proved to be the consequence of secondary introgression of highly divergent alleles. While only a few recombinants were observed at the Mytilin B locus, the MGD2 locus showed a high intragenic recombination rate, which increased in the exon coding for the mature peptide. In addition, standard neutrality tests revealed significant deviations from the mutation-drift equilibrium at both loci. These results suggest that either balancing or directional selection is likely to play a role in the evolution of the two AMPs and introgression would be adaptive. However, evidence accumulated at the Mytilin B locus allows neither for identification of the direction of selection nor for any conclusions on whether selection acted directly on the antimicrobial peptide itself. At the MGD2 locus, a spatial variation of polymorphism patterns along the sequence suggests that selection was direct, although the precise nature of the selection (directional vs. balancing) remains unclear. This study concurs with previous reports of an effect of slight selection on AMP genes evolution in other invertebrates, although selection does not necessarily act on the mature peptides.

Genetic hitchhiking in a subdivided population of Mytilus edulis

BACKGROUND

Few models of genetic hitchhiking in subdivided populations have been developed and the rarity of empirical examples is even more striking. We here provide evidences of genetic hitchhiking in a subdivided population of the marine mussel Mytilus edulis. In the Bay of Biscay (France), a patch of M. edulis populations happens to be separated from its North Sea conspecifics by a wide region occupied only by the sister species M. galloprovincialis. Although genetic differentiation between the two M. edulis regions is largely non-significant at ten marker loci (average FST~0.007), a strong genetic differentiation is observed at a single locus (FST = 0.25). We validated the outlier status of this locus, and analysed DNA sequence polymorphism in order to identify the nature of the selection responsible for the unusual differentiation.

RESULTS

We first showed that introgression of M. galloprovincialis alleles was very weak in both populations and did not significantly affect their differentiation. Secondly, we observed the genetic signature of a selective sweep within both M. edulis populations in the form of a star-shaped clade of alleles. This clade was nearly fixed in the North Sea and was segregating at a moderate frequency in the Bay of Biscay, explaining their genetic differentiation. Incomplete fixation reveals that selection was not direct on the locus but that the studied sequence recombined with a positively selected allele at a linked locus while it was on its way to fixation. Finally, using a deterministic model we showed that the wave of advance of a favourable allele at a linked locus, when crossing a strong enough barrier to gene flow, generates a step in neutral allele frequencies comparable to the step observed between the two M. edulis populations at the outlier locus. In our case, the position of the barrier is now materialised by a large patch of heterospecific M. galloprovincialis populations.

CONCLUSION

High FST outlier loci are usually interpreted as being the consequence of ongoing divergent local adaptation. Combining models and data we show that among-population differentiation can also dramatically increase following a selective sweep in a structured population. Our study illustrates how a striking geographical pattern of neutral diversity can emerge from past indirect hitchhiking selection in a structured population. NOTE: Nucleotide sequences reported in this paper are available in the GenBanktrade mark database under the accession numbers EU684165 - EU684228.

Replaying the tape: recurring biogeographical patterns in Cape Verde Conus after 12 million years

Isolated oceanic islands are excellent natural laboratories to test the relative role of historical contingency and determinism in evolutionary diversification. Endemics of the marine venomous snail Conus in the Cape Verde archipelago were originated from at least two independent colonizations of 'small' and 'large' shelled species separated by 12 million years. In this study, we have reconstructed phylogenetic relationships within large-shelled Conus (C. ateralbus, C. pseudonivifer, C. trochulus, and C. venulatus) based on mitochondrial cox1 and nad4 haplotype sequences. The reconstructed molecular phylogeny revealed three well-supported and relatively divergent clades (A, B, and C) that do not correspond to current species classification based on shell colour and banding patterns. Clade A grouped specimens assigned either to C. pseudonivifer or C. trochulus, clade B is composed of specimens assigned to C. venulatus, and clade C comprises specimens assigned either to C. venulatus or C. ateralbus. Geometric morphometric analyses found significant differences between the radular teeth shape of C. pseudonivifer/C. trochulus and C. venulatus/C. ateralbus. In clades A and B, northwestern Boavista and Maio specimens cluster together to the exclusion of eastern Boavista samples. In Sal, populations form a monophyletic island assemblage (clade C). The large-shelled Conus have remarkably replicated biogeographical patterns of diversification of small-shelled Conus. Similar selective forces (i.e. nonplanktonic lecithotrophy with limited larval dispersal and allopatric diversification) together with repeated instances of low sea level stands during glacial maxima that allowed connection between islands, have overcome the effect of historical contingency, and explain the observed recurring biogeographical patterns.

Evolution of 2-DE protein patterns in a mussel hybrid zone

Protein expression patterns were compared in mussels from a hybrid zone between the two species Mytilus edulis and M. galloprovincialis using 2-DE. Significant differences in expression pattern were observed between species and between mussels within species. Hybrid mussels had more variable protein expression patterns than mussels of each species. This could be due to segregation at expression modifier loci in the hybrids. It is proposed that unusual hybrid expression patterns might contribute to postzygotic isolation, if such patterns reduce fitness. The use of proteomic data for testing evolutionary models, as in many transcriptomics studies, is explored with results consistent with the expectations of neutral theory.

Mitochondrial DNA phylogeography and mating compatibility reveal marked genetic structuring and speciation in the NE Atlantic bryozoan Celleporella hyalina

The marine bryozoan Celleporella hyalina is a species complex composed of many highly divergent and mostly allopatric genetic lineages that are reproductively isolated but share a remarkably similar morphology. One such lineage commonly encrusts macroalgae throughout the NE Atlantic coast. To explore the processes leading to geographical diversification, reproductive isolation and speciation in this taxon, we (i) investigated NE Atlantic C. hyalina mitochondrial DNA phylogeography, and (ii) used breeding trials between geographical isolates to ascertain reproductive isolation. We find that haplotype diversity is geographically variable and there is a strong population structure, with significant isolation by distance. NE Atlantic C. hyalina is structured into two main parapatric lineages that appear to have had independent Pleistocene histories. Range expansions have resulted in two contact zones in Spain and W Ireland. Lineage 1 is found from Ireland to Spain and has low haplotype diversity, with closely related haplotypes, suggesting a recent population expansion into the Irish Sea, S Ireland, S England and Spain. Lineage 2 is found from Iceland to Spain and has high haplotype diversity. Complete reproductive isolation was found between some geographical isolates representing both lineages, whereas it was incomplete or asymmetric between others, suggesting these latter phylogeographical groups probably represent incipient species. The phylogeographical distribution of NE Atlantic C. hyalina does not fall easily into a pattern of southern refugia, and we discuss likely differences between terrestrial and marine system responses to Pleistocene glacial cycles.

Geographic Variation and Positive Selection on M7 Lysin, an Acrosomal Sperm Protein in Mussels (Mytilus spp.)

Successful fertilization in free-spawning marine organisms depends on the interactions between genes expressed on the surfaces of eggs and sperm. Positive selection frequently characterizes the molecular evolution of such genes, raising the possibility that some common deterministic process drives the evolution of gamete recognition genes and may even be important for understanding the evolution of prezygotic isolation and speciation in the marine realm. One hypothesis is that gamete recognition genes are subject to selection for prezygotic isolation, namely, reinforcement. In a previous study, positive selection on the gene coding for the acrosomal sperm protein M7 lysin was demonstrated among allopatric populations of mussels in the Mytilus edulis species group (M. edulis, Mytilus galloprovincialis, and Mytilus trossulus). Here, we expand sampling to include M7 lysin haplotypes from populations where mussel species are sympatric and hybridize to determine whether there is a pattern of reproductive character displacement (RCD), which would be consistent with reinforcement driving selection on this gene. We do not detect a strong pattern of RCD; neither are there unique haplotypes in sympatry nor is there consistently greater population structure in comparisons involving sympatric populations. One distinct group of haplotypes, however, is strongly affected by natural selection, and this group of haplotypes is found within M. galloprovincialis populations throughout the Northern Hemisphere concurrent with haplotypes common to M. galloprovincialis and M. edulis. We suggest that balancing selection, perhaps resulting from sexual conflicts between sperm and eggs, maintains old allelic diversity within M. galloprovincialis.

Fitness landscapes support the dominance theory of post-zygotic isolation in the mussels Mytilus edulis and M. galloprovincialis

We studied the genetic basis of post-zygotic isolation in the marine mussels Mytilus edulis and Mytilus galloprovincialis. Evidence was obtained for a high number of recessive Dobzhansky-Muller substitutions in the genome of these two mussel taxa. We analysed the segregation of unlinked diagnostic markers in the progeny of two backcrosses and an F2 cross, 36 h and 200 days after fertilization. Directional selection favouring M. galloprovincialis genotypes was observed in both kinds of cross. In the F2, epistatic interactions between each pair of chromosome fragments mapped by the markers were identified in addition. Our results imply that homozygous-homozygous interactions are required for breakdown of coadaptation, in accordance with the dominance theory of post-zygotic isolation. Endogenous post-zygotic selection distributed over many loci throughout the genome provides the missing factor explaining the astonishing persistence and strength of barriers to neutral introgression in such a dispersive taxon as Mytilus.

Gamete compatibility and sperm competition affect paternity and hybridization between sympatric Asterias sea stars

Gamete interactions may strongly influence speciation and hybridization in sympatric broadcast-spawning marine invertebrates. We examined the role of gamete compatibility in species integrity using cross-fertilization studies between sympatric Asterias sea stars from a secondary contact zone in the northwest Atlantic. In crosses between single males and single females, gametes of both species were compatible and produced viable, fertile hybrid offspring, but with considerable variation in the receptivity of eggs to heterospecific sperm. Differential compatibility of heterospecific gametes was detected in sperm competition studies in which we used a nuclear DNA marker to assign paternity to larval offspring. Several families showed conspecific sperm precedence in A. forbesi eggs, and one family showed competitive superiority of A. forbesi sperm fertilizing A. rubens eggs. Gametic interactions are an important component of prezygotic reproductive isolation in sympatric Asterias. The interaction between gametes of these closely related sea stars is consistent with the function of gamete recognition systems that are known to mediate fertilization success and speciation in other marine invertebrates.

Disruption of doubly uniparental inheritance of mitochondrial DNA in hybrid mussels (Mytilus edulis x M. galloprovincialis)

Blue mussels of the genus Mytilus have an unusual mode of mitochondrial DNA inheritance termed doubly uniparental inheritance (DUI). Females are homoplasmic for the F mitotype which is inherited maternally, whereas males are heteroplasmic for this and the paternally inherited M mitotype. In areas where species distributions overlap a varying degree of hybridization occurs; yet genetic differences between allopatric populations are maintained. Observations from natural populations and previous laboratory experiments suggest that DUI may be disrupted by hybridization, giving rise to heteroplasmic females and homoplasmic males. We carried out controlled laboratory crosses between Mytilus edulis and M. galloprovincialis to produce pure species and hybrid larvae of known parentage. DNA markers were used to follow the fate of the F and M mitotypes through larval development. Disruption of the mechanism which determines whether the M mitotype is retained or eliminated occurred in an estimated 38% of M. edulis x M. galloprovincialis hybrid larvae, a level double that previously observed in adult mussels from a natural M. edulis x M. galloprovincialis hybrid population. Furthermore, reciprocal hybrid crosses exhibited contrasting types of DUI disruption. The results indicate that disruption of DUI in hybrid mussels may be associated with increased mortality and hence could be a factor in the maintenance of genetic integrity for each species.

Habitat preference and the marine-speciation paradox

Marine organisms challenge the classical theories of local adaptation and speciation because their planktonic larvae have the potential to maintain high gene flow. The marine-speciation paradox is illustrated by contact zones between incipient species that are so large that allopatric divergence seems unlikely. For this reason any mechanism preventing sympatric larvae of two incipient species from coexisting in the same habitats can be a powerful promoter of speciation. The contact zone between two hybridizing taxa of mussel, Mytilus edulis and M. galloprovincialis, in Europe provides an excellent example. Although the zone itself extends over thousands of kilometres, the opportunities for interbreeding are considerably reduced by the small-scale mosaic structure of the zone, where local patches of each taxon alternate at scales of kilometres or less, in response to locally variable ecological factors. Habitat choice by settling larvae would be a less costly mechanism than post-settlement selection to maintain such a mosaic structure. Unfortunately the role of selective settlement has remained hypothetical because larvae could not be scored by classical genetic markers. PCR markers allowed us to study larvae and settlement in ecologically contrasting sites within the zone. We show that only a subset of the genotypes present in the plankton settle in some sites, and that the adults on these sites show the same genetic bias. Genetically based variation in pre-settlement processes therefore accounts for the ecological segregation observed, though it is not the only factor involved in limiting successful interbreeding. The present dataset also supports previous reports of partial spawning asynchrony.