Nicolas Bierne

Introgression patterns in the mosaic hybrid zone between Mytilus edulis and M. galloprovincialis

Hybrid zones are fascinating systems to investigate the structure of genetic barriers. Marine hybrid zones deserve more investigation because of the generally high dispersion potential of planktonic larvae which allows migration on scales unrivalled by terrestrial species. Here we analyse the genetic structure of the mosaic hybrid zone between the marine mussels Mytilus edulis and M. galloprovincialis, using three length‐polymorphic PCR loci as neutral and diagnostic markers on 32 samples along the Atlantic coast of Europe. Instead of a single genetic gradient from M. galloprovincialis on the Iberian Peninsula to M. edulis populations in the North Sea, three successive transitions were observed in France. From South to North, the frequency of alleles typical of M. galloprovincialis first decreases in the southern Bay of Biscay, remains low in Charente, then increases in South Brittany, remains high in most of Brittany, and finally decreases again in South Normandy. The two enclosed patches observed in the midst of the mosaic hybrid zone in Charente and Brittany, although predominantly M. edulis‐like and M. galloprovincialis‐like, respectively, are genetically original in two respects. First, considering only the various alleles typical of one species, the patches show differentiated frequencies compared to the reference external populations. Second, each patch is partly introgressed by alleles of the other species. When introgression is taken into account, linkage disequilibria appear close to their maximum possible values, indicating a strong genetic barrier within all transition zones. Some pre‐ or postzygotic isolation mechanisms (habitat specialization, spawning asynchrony, assortative fertilization and hybrid depression) have been documented in previous studies, although their relative importance remains to be evaluated. We also provided evidence for a recent migratory ‘short‐cut’ connecting M. edulis‐like populations of the Charente patch to an external M. edulis population in Normandy and thought to reflect artificial transfer of spat for aquaculture.

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.

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

Abstract Assortative mating (prezygotic isolation) and reduced hybrid fitness (postzygotic isolation) are typically invoked to explain the stability of hybrid zones. In the tension zone model, these factors work in opposition to migration, which promotes genetic homogeneity. Many marine animals migrate over long distances through a planktonic larval stage. Therefore, strong reproductive isolation is needed to maintain stable marine hybrid zones. However, surprisingly little is known about mating preferences and hybrid fitness in marine organisms. Smooth‐shelled mussels (Mytilus spp.) form a well‐known species complex, with hybridization over extensive areas such as the contact zone of M. edulis and M. galloprovincialis around European Atlantic coasts. This paper reports direct experimental evidence of assortative fertilization, hybrid larval inviability, and early heterosis for growth rate in M. edulis and M. galloprovincialis. Four crosses between pure M. edulis and M. galloprovincialis were analyzed with a new polymerasechain‐reaction‐based diagnostic marker. Gamete competition between taxa was allowed in two out of the four crosses. Genotype frequencies observed at an early stage (36 h after fertilization) unambiguously revealed assortative fertilization when gamete competition was allowed. A significant reduction in hybrid viability was subsequently observed during the larval stage. At the same stage an antagonistic effect, heterosis, was observed on growth rate. However, even if heterosis is observed in the F1, it is expected to vanish in subsequent hybrid generations. Although specialization for different habitats and asynchronous spawning have been mentioned as factors contributing to the maintenance of the blue mussel hybrid zone in Europe, we argue that assortative fertilization and reduced hybrid fitness are important factors that also contribute to the stabilization of this zone. These results emphasize that multiple factors may act concomitantly in a barrier to gene flow, especially in complex life cycles. Furthermore, they show that assortative mating through gamete preference, as already demonstrated for sea urchins, may play a role in speciation processes taking place in the sea.

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.

Deleterious mutations in a hybrid zone: can mutational load decrease the barrier to gene flow?

The aim of this paper is to investigate the effect of deleterious mutations in a hybrid zone maintained by selection against hybrids. In such zones, linkage disequilibria among hybrid depression loci, resulting from a balance between migration and selection, are crucial in maintaining the barrier because they allow each locus, in addition to its own selection coefficient, to cumulate indirect selective effects from other loci. Deleterious alleles produce heterosis and increase by this means the effective migration rate in structured populations. In a hybrid zone, they therefore contribute to decrease linkage disequilibria as well as the barrier to gene flow imposed by hybrid depression. However, deleterious mutations have no effect: (i) when selection against hybrids is weak, because linkage disequilibria are small even without heterosis in this case, or (ii) when selection against hybrids is so strong that it overwhelms heterosis. On the other hand, with moderate selection against hybrids, the decrease in the strength of the barrier due to heterosis may reach detectable levels, although it requires relatively small population sizes and/or migration rates. The effect is expected to be small and only within small genomes where loci are tightly linked can it become strong. Nevertheless, neglecting mutational load may to some extent obscure the estimations of selective parameters based either on artificial F1 crosses or on cline characteristics.

Codon Usage Bias in Animals: Disentangling the Effects of Natural Selection, Effective Population Size, and GC-Biased Gene Conversion

Abstract Selection on codon usage bias is well documented in a number of microorganisms. Whether codon usage is also generally shaped by natural selection in large organisms, despite their relatively small effective population size (Ne), is unclear. In animals, the population genetics of codon usage bias has only been studied in a handful of model organisms so far, and can be affected by confounding, nonadaptive processes such as GC‐biased gene conversion and experimental artefacts. Using population transcriptomics data, we analyzed the relationship between codon usage, gene expression, allele frequency distribution, and recombination rate in 30 nonmodel species of animals, each from a different family, covering a wide range of effective population sizes. We disentangled the effects of translational selection and GC‐biased gene conversion on codon usage by separately analyzing GC‐conservative and GC‐changing mutations. We report evidence for effective translational selection on codon usage in large‐Ne species of animals, but not in small‐Ne ones, in agreement with the nearly neutral theory of molecular evolution. C‐ and T‐ending codons tend to be preferred over synonymous G‐ and A‐ending ones, for reasons that remain to be determined. In contrast, we uncovered a conspicuous effect of GC‐biased gene conversion, which is widespread in animals and the main force determining the fate of AT↔GC mutations. Intriguingly, the strength of its effect was uncorrelated with Ne.

Analysis of genome-wide differentiation between native and introduced populations of the cupped oysters Crassostrea gigas and Crassostrea angulata

Abstract The Pacific cupped oyster is genetically subdivided into two sister taxa, Crassostrea gigas and Crassostrea angulata, which are in contact in the north-western Pacific. The nature and origin of their genetic and taxonomic differentiation remains controversial due the lack of known reproductive barriers and the high degree of morphologic similarity. In particular, whether the presence of ecological and/or intrinsic isolating mechanisms contributes to species divergence is unknown. The recent co-introduction of both taxa into Europe offers a unique opportunity to test how genetic differentiation is maintained under new environmental and demographic conditions. We generated a pseudochromosome assembly of the Pacific oyster genome using a combination of BAC-end sequencing and scaffold anchoring to a new high-density linkage map. We characterized genome-wide differentiation between C. angulata and C. gigas in both their native and introduced ranges, and showed that gene flow between species has been facilitated by their recent co-introductions in Europe. Nevertheless, patterns of genomic divergence between species remain highly similar in Asia and Europe, suggesting that the environmental transition caused by the co-introduction of the two species did not affect the genomic architecture of their partial reproductive isolation. Increased genetic differentiation was preferentially found in regions of low recombination. Using historical demographic inference, we show that the heterogeneity of differentiation across the genome is well explained by a scenario whereby recent gene flow has eroded past differentiation at different rates across the genome after a period of geographical isolation. Our results thus support the view that low-recombining regions help in maintaining intrinsic genetic differences between the two species.