The molecular background of the aspartate aminotransferase polymorphism in Littorina snails maintained by strong selection on small spatial scales.

Allozymes present several classical examples of divergent selection, including the variation in the cytosolic aspartate aminotransferase (AAT) in the intertidal snails Littorina saxatilis. AAT is a part of the asparate-malate shuttle, in the interidal molluscs involved in the anaerobic respiration during desiccation. Previous allozyme studies reported the sharp gradient in the frequencies of the AAT100and the AAT120 alleles between the low and high shores in the Northern Europe and the differences in their enzymatic activity, supporting the role of AAT in adaptation to desiccation. However, the populations in the Iberian Peninsula showed the opposite allele cline. Using the mRNA sequencing and the genome pool-seq analyses we characterize DNA sequences of the different AAT alleles, report the amino acid replacements behind the allozyme variation and show that same allozyme alleles in Northern and Southern populations have different protein sequences. Gene phylogeny reveals that the AAT100 and the northern AAT120 alleles represent the old polymorphism, shared among the closely related species of Littorina, while the southern AAT120 allele is more recently derived from AAT100. Further, we show that the Aat gene is expressed at constitutive level in different genotypes and conditions, supporting the role of structural variation in regulation of enzyme activity. Finally, we report the location and the structure of the gene in the L. saxatilis genome and the presence of two additional non-functional gene copies. Altogether, we provide a missing link between the classical allozyme studies and the genome scans and bring together the results produced over decades of the genetic research.

An allozyme polymorphism is associated with a large chromosomal inversion in the marine snail Littorina fabalis.

Understanding the genetic targets of natural selection is one of the most challenging goals of population genetics. Some of the earliest candidate genes were identified from associations between allozyme allele frequencies and environmental variation. One such example is the clinal polymorphism in the arginine kinase (Ak) gene in the marine snail Littorina fabalis. While other enzyme loci do not show differences in allozyme frequencies among populations, the Ak alleles are near differential fixation across repeated wave exposure gradients in Europe. Here, we use this case to illustrate how a new sequencing toolbox can be employed to characterize the genomic architecture associated with historical candidate genes. We found that the Ak alleles differ by nine nonsynonymous substitutions, which perfectly explain the different migration patterns of the allozymes during electrophoresis. Moreover, by exploring the genomic context of the Ak gene, we found that the three main Ak alleles are located on different arrangements of a putative chromosomal inversion that reaches near fixation at the opposing ends of two transects covering a wave exposure gradient. This shows Ak is part of a large (3/4 of the chromosome) genomic block of differentiation, in which Ak is unlikely to be the only target of divergent selection. Nevertheless, the nonsynonymous substitutions among Ak alleles and the complete association of one allele with one inversion arrangement suggest that the Ak gene is a strong candidate to contribute to the adaptive significance of the inversion.

Genetic and morphological divergence between Littorina fabalis ecotypes in Northern Europe.

Low dispersal marine intertidal species facing strong divergent selective pressures associated with steep environmental gradients have a great potential to inform us about local adaptation and reproductive isolation. Among these, gastropods of the genus Littorina offer a unique system to study parallel phenotypic divergence resulting from adaptation to different habitats related with wave exposure. In this study, we focused on two Littorina fabalis ecotypes from Northern European shores and compared patterns of habitat-related phenotypic and genetic divergence across three different geographic levels (local, regional and global). Geometric morphometric analyses revealed that individuals from habitats moderately exposed to waves usually present a larger shell size with a wider aperture than those from sheltered habitats. The phenotypic clustering of L. fabalis by habitat across most locations (mainly in terms of shell size) support an important role of ecology in morphological divergence. A genome scan based on amplified fragment length polymorphisms (AFLPs) revealed a heterogeneous pattern of differentiation across the genome between populations from the two different habitats, suggesting ecotype divergence in the presence of gene flow. The contrasting patterns of genetic structure between nonoutlier and outlier loci, and the decreased sharing of outlier loci with geographic distance among locations are compatible with parallel evolution of phenotypic divergence, with an important contribution of gene flow and/or ancestral variation. In the future, model-based inference studies based on sequence data across the entire genome will help unravelling these evolutionary hypotheses, improving our knowledge about adaptation and its influence on diversification within the marine realm.

Phylogeographic history of flat periwinkles, Littorina fabalis and L. obtusata.

BACKGROUND: The flat periwinkles, Littorina fabalis and L. obtusata, are two sister species widely distributed throughout the Northern Atlantic shores with high potential to inform us about the process of ecological speciation in the intertidal. However, whether gene flow has occurred during their divergence is still a matter of debate. A comprehensive assessment of the genetic diversity of these species is also lacking and their main glacial refugia and dispersal barriers remain largely unknown. In order to fill these gaps, we sequenced two mitochondrial genes and two nuclear fragments to perform a phylogeographic analysis of flat periwinkles across their distribution range. RESULTS: We identified two main clades largely composed by species-specific haplotypes corresponding to L. obtusata and L. fabalis, with moderate to strong support, respectively. Importantly, a model of divergence with gene flow between the two species (from L. obtusata to L. fabalis) was better supported, both in Iberia and in northern-central Europe. Three mitochondrial clades were detected within L. fabalis and two within L. obtusata, with strong divergence between Iberia and the remaining populations. The largest component of the genetic variance within each species was explained by differences between geographic regions associated with these clades. Our data suggests that overall intraspecific genetic diversity is similar between the two flat periwinkle species and that populations from Iberia tend to be less diverse than populations from northern-central Europe. CONCLUSIONS: The phylogeographic analysis of this sister-species pair supports divergence with gene flow. This system thus provides us with the opportunity to study the contribution of gene flow and natural selection during diversification. The distribution of the different clades suggests the existence of glacial refugia in Iberia and northern-central Europe for both species, with a main phylogeographic break between these regions. Although the genetic diversity results are not fully conclusive, the lower diversity observed in Iberia could reflect marginal conditions at the southern limit of their distribution range during the current interglacial period.