Mitochondrial DNA CoI haplotype variation in sibling species of rough periwinkles

Two new pseudocryptic species in the medium-sized common European land snails, Fruticicola Held, 1838; as a result of phylogeographic analysis of Fruticicola fruticum (O. F. Müller, 1774) (Gastropoda: Helicoidea: Camaenidae)

Fruticicola fruticum (O. F. Müller, 1774), a medium-sized helicoid snail in the Bradybaenidae, has a wide range in Europe, reaching from the Urals and the Caucasus to the Balkans, and from the southern part of Scandinavia, through Central Europe to eastern and central France and northern Italy. There are numerous studies on its distribution, biology, life cycle, etc., but little is known about the genetic diversity of this taxon. Here, we studied the phylogeny and phylogeography of F. fruticum using two mitochondrial markers: cytochrome oxidase subunit I (COI) and 16S ribosomal RNA (16S); and four nuclear markers: 18S ribosomal RNA (18S), 28S ribosomal RNA (28S), internal transcribed spacer (ITS-2), and histone 3 (H3). The study was based on 59 populations sampled across the range. Whereas nuclear markers showed little differentiation, phylogenetic analysis of COI sequences clearly confirmed the distinctness of the European Fruticicola and Asian Bradybaena (p-distance 0.229). Within Fruticicola 54 haplotypes were detected, haplotype diversity (Hd) = 0.973±0.006, nucleotide diversity (π) = 0.137±0.005. ABGD and PTP delimitation analyzes distinguished eight mOTUs. Two sequences (our mOTU C) from Russia were published in the GenBank as two distinct species: F. schrenckii and F. transbaicalia. Seven further mOTUs identified in our study formed three distinct lineages, regarded as species. The first (mOTU A and mOTU B), represented by 40 populations, occupies a wide range across northern and central Europe, extending east to Ukraine and south to northern Croatia (mOTU B). It encompasses the type locality of F. fruticum, and can be recognized as F. fruticum sensu stricto. Another lineage (mOTU D and mOTU E), represented by six populations in central Romania, appears to form another species. Both mOTUs were found together in one population. A third lineage, containing mOTUs F, G and H, represented by 14 populations, was distributed across the Balkans from N.E. Croatia to Bulgaria. p-distances between the three species ranged from 0.172 to 0.219, and between all the mOTUs, pooled together, from 0.172 to 0.258. The highest genetic diversity was found in species 3 (0.112) and the lowest in species 1 (0.025), despite its largest geographic distribution. Pairwise p-distances, Tamura 3-parameter distances, composite likelihood distances, as well as the coancestry coefficient F_ST, calculated for all populations pooled together were significantly associated with geographic distance, but this was not the case within each of these three species. The significant association for all populations reflected high diversity between the species coupled with high geographic distances between their populations, not the character of intraspecies diversity. With a few exceptions, there hold a rather infinite island model with low migration. AMOVA detected 78% of the variance between the three species, 18% among populations within the species, and only 3.6% within the populations. The low genetic diversity of widespread F. fruticum s. stricto, compared with much higher diversity of two narrowly distributed newly found species of Fruticicola, may reflect the rapid spread of the former into previously uninhabitable regions, while the latter were able to maintain populations in glacial refugia. The estimated time of divergence between the three species, 1.7-2.19 mya, suggests their ancestors' isolation in southern European refugia during the lower Pleistocene, the Gelasian/Calabrian. There was no clear association of variation in shell morphology and lineage or mOTU identity; on external characters, these species are semicryptic, subtle differences in reproductive anatomy among them were found.

Proteomic similarity of the Littorinid snails in the evolutionary context

BACKGROUND

The introduction of DNA-based molecular markers made a revolution in biological systematics. However, in cases of very recent divergence events, the neutral divergence may be too slow, and the analysis of adaptive part of the genome is more informative to reconstruct the recent evolutionary history of young species. The advantage of proteomics is its ability to reflect the biochemical machinery of life. It may help both to identify rapidly evolving genes and to interpret their functions.

METHODS

Here we applied a comparative gel-based proteomic analysis to several species from the gastropod family Littorinidae. Proteomes were clustered to assess differences related to species, geographic location, sex and body part, using data on presence/absence of proteins in samples and data on protein occurrence frequency in samples of different species. Cluster support was assessed using multiscale bootstrap resampling and the stability of clustering-using cluster-wise index of cluster stability. Taxon-specific protein markers were derived using IndVal method. Proteomic trees were compared to consensus phylogenetic tree (based on neutral genetic markers) using estimates of the Robinson-Foulds distance, the Fowlkes-Mallows index and cophenetic correlation.

RESULTS

Overall, the DNA-based phylogenetic tree and the proteomic similarity tree had consistent topologies. Further, we observed some interesting deviations of the proteomic littorinid tree from the neutral expectations. (1) There were signs of molecular parallelism in two Littoraria species that phylogenetically are quite distant, but live in similar habitats. (2) Proteome divergence was unexpectedly high between very closely related Littorina fabalis and L. obtusata, possibly reflecting their ecology-driven divergence. (3) Conservative house-keeping proteins were usually identified as markers for cryptic species groups ("saxatilis" and "obtusata" groups in the Littorina genus) and for genera (Littoraria and Echinolittorina species pairs), while metabolic enzymes and stress-related proteins (both potentially adaptively important) were often identified as markers supporting species branches. (4) In all five Littorina species British populations were separated from the European mainland populations, possibly reflecting their recent phylogeographic history. Altogether our study shows that proteomic data, when interpreted in the context of DNA-based phylogeny, can bring additional information on the evolutionary history of species.

DNA barcoding of crickets, katydids and grasshoppers (Orthoptera) from Central Europe with focus on Austria, Germany and Switzerland

We present a DNA barcoding study on the insect order Orthoptera that was generated in collaboration between four barcoding projects in three countries, viz. Barcoding Fauna Bavarica (Germany), German Barcode of Life, Austrian Barcode of Life and Swiss Barcode of Life. Our data set includes 748 COI sequences from 127 of the 162 taxa (78.4%) recorded in the three countries involved. Ninety-three of these 122 species (76.2%, including all Ensifera) can be reliably identified using DNA barcodes. The remaining 26 caeliferan species (families Acrididae and Tetrigidae) form ten clusters that share barcodes among up to five species, in three cases even across different genera, and in six cases even sharing individual barcodes. We discuss incomplete lineage sorting and hybridization as most likely causes of this phenomenon, as the species concerned are phylogenetically young and hybridization has been previously observed. We also highlight the problem of nuclear mitochondrial pseudogenes (numts), a known problem in the barcoding of orthopteran species, and the possibility of Wolbachia infections. Finally, we discuss the possible taxonomic implications of our barcoding results and point out future research directions.

Species and gene divergence in Littorina snails detected by array comparative genomic hybridization

Background

Array comparative genomic hybridization (aCGH) is commonly used to screen different types of genetic variation in humans and model species. Here, we performed aCGH using an oligonucleotide gene-expression array for a non-model species, the intertidal snail Littorina saxatilis. First, we tested what types of genetic variation can be detected by this method using direct re-sequencing and comparison to the Littorina genome draft. Secondly, we performed a genome-wide comparison of four closely related Littorina species: L. fabalis, L. compressa, L. arcana and L. saxatilis and of populations of L. saxatilis found in Spain, Britain and Sweden. Finally, we tested whether we could identify genetic variation underlying “Crab” and “Wave” ecotypes of L. saxatilis.

Results

We could reliably detect copy number variations, deletions and high sequence divergence (i.e. above 3%), but not single nucleotide polymorphisms. The overall hybridization pattern and number of significantly diverged genes were in close agreement with earlier phylogenetic reconstructions based on single genes. The trichotomy of L. arcana, L. compressa and L. saxatilis could not be resolved and we argue that these divergence events have occurred recently and very close in time. We found evidence for high levels of segmental duplication in the Littorina genome (10% of the transcripts represented on the array and up to 23% of the analyzed genomic fragments); duplicated genes and regions were mostly the same in all analyzed species. Finally, this method discriminated geographically distant populations of L. saxatilis, but we did not detect any significant genome divergence associated with ecotypes of L. saxatilis.

Conclusions

The present study provides new information on the sensitivity and the potential use of oligonucleotide arrays for genotyping of non-model organisms. Applying this method to Littorina species yields insights into genome evolution following the recent species radiation and supports earlier single-gene based phylogenies. Genetic differentiation of L. saxatilis ecotypes was not detected in this study, despite pronounced innate phenotypic differences. The reason may be that these differences are due to single-nucleotide polymorphisms.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-687) contains supplementary material, which is available to authorized users.

Micro-spatial distribution of two sibling periwinkle species across the intertidal indicates hybrdization

Populations of periwinkles Littorina saxatilis (Olivi 1792) and L. arcana Hannaford Ellis, 1978 are well suited for microevolutionary studies, being at the same time closely related and intraspecifically diverse. The divergence between these two sibling species, sympatric over large parts of their distribution areas, is small, the only morphological difference being the pallial gland complex structure in females. Molecular identification is possible with the use of a RAPD nuclear marker (cloned A2.8 DNA fragment) typical for L. arcana. However, in some individuals from sympatric populations molecular and morphological criteria suggest conflicting species affiliation, which may be explained either by hybridization or by shared ancestral polymorphism. We tested the hybridization hypotheses examining the micro-spatial distribution of these two species across the intertidal zone in two distant sites at the Barents Sea. We found that (a) the frequency of putative hybrids in sympatric populations was proportional to the frequency of L. arcana; (b) L. saxatilis bearing A2.8 DNA fragment were almost absent in the lower part of the intertidal zone, where L. arcana was absent too; (c) there was a close positive correlation between the distribution of potential parent molluscs and putative hybrids. Moreover, logistic regression models showed a good agreement between the distribution of putative hybrid frequencies and that of parental species frequencies. All our observations taken together support the hypothesis of hybridization between L. saxatilis and L. arcana. Elucidating the mechanisms that support the species status of these sympatric populations is necessary.

Phylogeographic analysis reveals a deep lineage split within North Atlantic Littorina saxatilis

Phylogeographic studies provide critical insight into the evolutionary histories of model organisms; yet, to date, range-wide data are lacking for the rough periwinkle Littorina saxatilis, a classic example of marine sympatric speciation. Here, we use mitochondrial DNA (mtDNA) sequence data to demonstrate that L. saxatilis is not monophyletic for this marker, but is composed of two distinct mtDNA lineages (I and II) that are shared with sister species Littorina arcana and Littorina compressa. Bayesian coalescent dating and phylogeographic patterns indicate that both L. saxatilis lineages originated in the eastern North Atlantic, around the British Isles, at approximately 0.64 Ma. Both lineages are now distributed broadly across the eastern, central and western North Atlantic, and show strong phylogeographic structure among regions. The Iberian Peninsula is genetically distinct, suggesting prolonged isolation from northeastern North Atlantic populations. Western North Atlantic populations of L. saxatilis lineages I and II predate the last glacial maximum and have been isolated from eastern North Atlantic populations since that time. This identification of two distinct, broadly distributed mtDNA lineages further complicates observed patterns of repeated incipient ecological speciation in L. saxatilis, because the sympatric origins of distinct ecotype pairs on eastern North Atlantic shores may be confounded by admixture of divergent lineages.

A potential species-specific molecular marker suggests interspecific hybridization between sibling species Littorina arcana and L. saxatilis (Mollusca, Caenogastropoda) in natural populations

Three sister species of rough periwinkles, viz. Littorina saxatilis (Olivi 1792), L. arcana (Hannaford Ellis 1978) and L. compressa (Jeffreys 1865) from the Barents Sea (Russia), the White Sea (Russia) and the Norwegian Sea (Norway) were studied. The identification of two sibling species L. saxatilis and L. arcana is often difficult as both species have extremely similar shell morphology and reproductive systems. Only mature females can be unambiguously distinguished, with a jelly gland present in female L. arcana, but which is replaced by a brood pouch containing developing embryos in L. saxatilis. No clear-cut diagnostic features have been found to discriminate between males or juveniles of the two species. The very first diagnostic DNA marker (DNA fragment A2.8, 271 bp length) for L. arcana and L. saxatilis separation was developed. The marker was derived from apparently species-specific L. arcana DNA fragments obtained via Random Amplified Polymorphic DNA (RAPD) analysis. This fragment was cloned and sequenced, whereupon specific primers were designed and the amplification was surveyed in a large number of morphologically well-identified females of both species. Subsequently, the specific DNA marker was used for the identification of male L. arcana and partners in copulating pairs. In this way, we obtained evidence of possible interspecific hybridization between the sibling species L. arcana and L. saxatilis living in sympatry in natural populations: the presence of A2.8 fragment in 12% of morphologically well identified L. saxatilis females and its absence in 14% of morphologically well identified L. arcana females. The A2.8 fragment never amplified in L. saxatilis from sites without L. arcana. The A2.8 fragment did not amplify in L. compressa, not even in microsympatric populations, and we did not observe interspecific copulations between L. arcana and L. compressa.

Multiple introductions promote range expansion of the mollusc Cyclope neritea (Nassariidae) in France: evidence from mitochondrial sequence data

Since the 1970s, the nassariid gastropod Cyclope neritea has been extending its range north along the French Atlantic coasts from the Iberian Peninsula. This may be due to natural spread because of the recent warming of the northeastern Atlantic. However, human-mediated introductions related to shellfish culture may also be a probable explanation for this sudden range expansion. To examine these two hypotheses, we carried out a comprehensive study based on mitochondrial gene sequences (cytochrome oxidase I) of the five recently colonized French bays as well as 14 populations located in the recognized native range of the species. From a total of 594 individuals, we observed 29 haplotypes to split into three divergent clades. In the native range, we observed a low molecular diversity, strong genetic structure and agreement between geography and gene genealogies. Along the French coasts, we observed the opposite: high genetic diversity and low genetic structure. Our results show that recurrent human-mediated introductions from several geographical areas in the native range may be a source for the French Atlantic populations. However, despite the low dispersal ability of C. neritea, the isolation-by-distance pattern in France suggested that this gastropod may have been present (although unnoticed) on the French Atlantic coasts before the 1970s. As C. neritea shows characteristics of a cryptogenic species, the classification of Atlantic populations as either native or introduced is not straightforward. Cryptogenic species should be studied further to determine the status of new populations close to their recognized native range.

Invertebrate species with nonpelagic larvae have elevated levels of nonsynonymous substitutions and reduced nucleotide diversities

Under a nearly neutral model in which most amino acid substitutions are slightly deleterious, variation in demography, population structure, and other ecological factors among closely related species can potentially modify the effective population size or the selective regime, leading to differences in the rate of nonsynonymous substitution. Ratios of nonsynonymous to synonymous substitutions (d(N)/d(S)) between species were analyzed in a sea star genus (Patiriella) and a molluscan genus (Littorina), each with diverse modes of reproduction, including multiple lineages with pelagic and nonpelagic larvae. In both genera, lineages with nonpelagic larvae had significantly higher d(N)/d(S) ratios than lineages with pelagic larvae. The hypothesis that the elevated d(N)/d(S) ratios in species with nonpelagic larvae was due to reduced effective population size was tested by comparing nucleotide diversities in three genera of gastropod mollusks (Littorina, Crepidula, and Hydrobia), each with several modes of reproduction. Overall, there was a significant (p < 0.05) reduction in nucleotide diversity in species with nonpelagic larvae compared to species with pelagic larvae.

Nuclear DNA restriction site polymorphisms and the phylogeny and population structure of an intertidal snail species complex (Littorina)

Primers for amplification of four novel, unlinked nuclear DNA loci, the first reported for the rough periwinkles of the genus Littorina, are described. Patterns of restriction site polymorphism for these loci are detailed within the rough periwinkles. RFLPs are not found to be diagnostic for any of the currently accepted species within this group, nor for any of the contentious subspecies, or forms, whose taxonomic status is uncertain. However, there are important differences in allele frequencies between these taxa and certain of these mirror differences detected in a previous study of the mitochondrial DNA. These allele frequency data are used to construct a phylogeny in which groupings of the three recognised species are obvious when either Nei's genetic distances or Reynold's distances are clustered. Contentious forms (L. neglecta, L. saxatilis 'b' and L. tenebrosa) do not cluster as distinct taxa, although populations of L. neglecta have important allele frequency differences from L. saxatilis. These four loci have confirmed the consensus view of Littorina phylogeny and provided important information on population structure-however four loci is insufficient for reaching definitive conclusions. Since analysis of nuclear DNA polymorphisms such as these is invaluable for analysis of phylogeny, population structure and phylogeography, identification of additional loci is considered imperative.