Markus Pfenninger

Cryptic animal species are homogeneously distributed among taxa and biogeographical regions.

Cryptic species are two or more distinct but morphologically similar species that were classified as a single species. During the past two decades we observed an exponential growth of publications on cryptic species. Recently published reviews have demonstrated cryptic species have profound consequences on many biological disciplines. It has been proposed that their distribution is non-random across taxa and biomes. We analysed a literature database for the taxonomic and biogeographical distribution of cryptic animal species reports. Results from regression analysis indicate that cryptic species are almost evenly distributed among major metazoan taxa and biogeographical regions when corrected for species richness and study intensity. This indicates that morphological stasis represents an evolutionary constant and that cryptic metazoan diversity does predictably affect estimates of earths animal diversity. Our findings have direct theoretical and practical consequences for a number of prevailing biological questions with regard to global biodiversity estimates, conservation efforts and global taxonomic initiatives.BackgroundCryptic species are two or more distinct but morphologically similar species that were classified as a single species. During the past two decades we observed an exponential growth of publications on cryptic species. Recently published reviews have demonstrated cryptic species have profound consequences on many biological disciplines. It has been proposed that their distribution is non-random across taxa and biomes.ResultsWe analysed a literature database for the taxonomic and biogeographical distribution of cryptic animal species reports. Results from regression analysis indicate that cryptic species are almost evenly distributed among major metazoan taxa and biogeographical regions when corrected for species richness and study intensity.ConclusionThis indicates that morphological stasis represents an evolutionary constant and that cryptic metazoan diversity does predictably affect estimates of earths animal diversity. Our findings have direct theoretical and practical consequences for a number of prevailing biological questions with regard to global biodiversity estimates, conservation efforts and global taxonomic initiatives.

Phylogeographic history of the land snail Candidula unifasciata (Helicellinae, Stylommatophora): fragmentation, corridor migration, and secondary contact.

Abstract We studied sequence variation in 16S rDNA in 204 individuals from 37 populations of the land snail Candidula unifasciata (Poiret 1801) across the core species range in France, Switzerland, and Germany. Phylogeographic, nested clade, and coalescence analyses were used to elucidate the species evolutionary history. The study revealed the presence of two major evolutionary lineages that evolved in separate refuges in southeast France as result of previous fragmentation during the Pleistocene. Applying a recent extension of the nested clade analysis (Templeton 2001), we inferred that range expansions along river valleys in independent corridors to the north led eventually to a secondary contact zone of the major clades around the Geneva Basin. There is evidence supporting the idea that the formation of the secondary contact zone and the colonization of Germany might be postglacial events. The phylogeographic history inferred for C. unifasciata differs from general biogeographic patterns of postglacial colonization previously identified for other taxa, and it might represent a common model for species with restricted dispersal.

Short read Illumina data for the de novo assembly of a non-model snail species transcriptome (Radix balthica, Basommatophora, Pulmonata), and a comparison of assembler performance

BackgroundUntil recently, read lengths on the Solexa/Illumina system were too short to reliably assemble transcriptomes without a reference sequence, especially for non-model organisms. However, with read lengths up to 100 nucleotides available in the current version, an assembly without reference genome should be possible. For this study we created an EST data set for the common pond snail Radix balthica by Illumina sequencing of a normalized transcriptome. Performance of three different short read assemblers was compared with respect to: the number of contigs, their length, depth of coverage, their quality in various BLAST searches and the alignment to mitochondrial genes.ResultsA single sequencing run of a normalized RNA pool resulted in 16,923,850 paired end reads with median read length of 61 bases. The assemblies generated by VELVET, OASES, and SeqMan NGEN differed in the total number of contigs, contig length, the number and quality of gene hits obtained by BLAST searches against various databases, and contig performance in the mt genome comparison. While VELVET produced the highest overall number of contigs, a large fraction of these were of small size (< 200bp), and gave redundant hits in BLAST searches and the mt genome alignment. The best overall contig performance resulted from the NGEN assembly. It produced the second largest number of contigs, which on average were comparable to the OASES contigs but gave the highest number of gene hits in two out of four BLAST searches against different reference databases. A subsequent meta-assembly of the four contig sets resulted in larger contigs, less redundancy and a higher number of BLAST hits.ConclusionOur results document the first de novo transcriptome assembly of a non-model species using Illumina sequencing data. We show that de novo transcriptome assembly using this approach yields results useful for downstream applications, in particular if a meta-assembly of contig sets is used to increase contig quality. These results highlight the ongoing need for improvements in assembly methodology.

Utility of DNA taxonomy and barcoding for the inference of larval community structure in morphologically cryptic Chironomus (Diptera) species

Biodiversity studies require species level analyses for the accurate assessment of community structures. However, while specialized taxonomic knowledge is only rarely available for routine identifications, DNA taxonomy and DNA barcoding could provide the taxonomic basis for ecological inferences. In this study, we assessed the community structure of sediment dwelling, morphologically cryptic Chironomus larvae in the Rhine‐valley plain/Germany, comparing larval type classification, cytotaxonomy, DNA taxonomy and barcoding. While larval type classification performed poorly, cytotaxonomy and DNA‐based methods yielded comparable results: detrended correspondence analysis and permutation analyses indicated that the assemblages are not randomly but competitively structured. However, DNA taxonomy identified an additional species that could not be resolved by the traditional method. We argue that DNA‐based identification methods such as DNA barcoding can be a valuable tool to increase accuracy, objectivity and comparability of the taxonomic assessment in biodiversity and community ecology studies.

Inferring the past to predict the future: climate modelling predictions and phylogeography for the freshwater gastropod Radix balthica (Pulmonata, Basommatophora)

Understanding the impact of past climatic events on species may facilitate predictions of how species will respond to future climate change. To this end, we sampled populations of the common pond snail Radix balthica over the entire species range (northwestern Europe). Using a recently developed analytical framework that employs ecological niche modelling to obtain hypotheses that are subsequently tested with statistical phylogeography, we inferred the range dynamics of R. balthica over time. A Maxent modelling for present‐day conditions was performed to infer the climate envelope for the species, and the modelled niche was used to hindcast climatically suitable range at the last glacial maximum (LGM) c. 21 000 years ago. Ecological niche modelling predicted two suitable areas at the LGM within the present species range. Phylogeographic model selection on a COI mitochondrial DNA data set confirmed that R. balthica most likely spread from these two disjunct refuges after the LGM. The match observed between the potential range of the species at the LGM given its present climatic requirements and the phylogeographically inferred refugial areas was a clear argument in favour of niche conservatism in R. balthica, thus allowing to predict the future range. The subsequent projection of the potential range under a global change scenario predicts a moderate pole‐ward shift of the northern range limits, but a dramatic loss of areas currently occupied in France, western Great Britain and southern Germany.

Comparing the efficacy of morphologic and DNA-based taxonomy in the freshwater gastropod genus Radix (Basommatophora, Pulmonata)

BackgroundReliable taxonomic identification at the species level is the basis for many biological disciplines. In order to distinguish species, it is necessary that taxonomic characters allow for the separation of individuals into recognisable, homogeneous groups that differ from other such groups in a consistent way. We compared here the suitability and efficacy of traditionally used shell morphology and DNA-based methods to distinguish among species of the freshwater snail genus Radix (Basommatophora, Pulmonata).ResultsMorphometric analysis showed that shell shape was unsuitable to define homogeneous, recognisable entities, because the variation was continuous. On the other hand, the Molecularly defined Operational Taxonomic Units (MOTU), inferred from mitochondrial COI sequence variation, proved to be congruent with biological species, inferred from geographic distribution patterns, congruence with nuclear markers and crossing experiments. Moreover, it could be shown that the phenotypically plastic shell variation is mostly determined by the environmental conditions experienced.ConclusionContrary to DNA-taxonomy, shell morphology was not suitable for delimiting and recognising species in Radix. As the situation encountered here seems to be widespread in invertebrates, we propose DNA-taxonomy as a reliable, comparable, and objective means for species identification in biological research.

Rangewide phylogeography of a terrestrial slug in Europe: evidence for Alpine refugia and rapid colonization after the Pleistocene glaciations

Intraspecific phylogeographical patterns largely depend on the life history traits of a species. Especially species with a high degree of cold tolerance, limited requirements towards habitat preferences, and relatively low active dispersal capacities may have responded in a different way to the Pleistocene climatological fluctuations than the majority of taxa studied so far. To evaluate this possibility, we studied Arion fuscus (Müller, 1774), a common and widespread European terrestrial slug, from 88 locations (N = 964). Sequence variation was assessed for fragments of the mitochondrial 16S rDNA and COI genes by means of single‐strand conformation polymorphisms (SSCP) and subsequent DNA sequencing. Additionally, eight allozyme loci were scored in 843 individuals. Phylogenetic analysis revealed the presence of two major evolutionary lineages, one in the Balkan region and another in the Alps and the rest of Europe. The sequence divergence between the two lineages was limited (3.3%), but gene flow between the regions was absent, suggesting that the two regions have been isolated since the late Pliocene or early Pleistocene. Allozyme differentiation among geographical regions and mitochondrial DNA (mtDNA) lineages was low. The geographical patterns observed in our data showed that (i) haplotype and nucleotide diversities are very low in northern Europe, suggesting that single haplotypes rapidly colonized large areas; (ii) recently expanded haplotype clades have restricted distribution ranges, suggesting that current gene flow is low; and (iii) genetic diversity in the Alps is much higher than in other regions and estimated past gene flow from the Eastern Alps to other regions was high, suggesting that this was a refugial zone during the Pleistocene. This full‐range phylogeography suggests the existence of an alternative refugial zone, situated north of the refugial areas currently recognized in most other taxa.

Postglacial recolonization at a snail's pace (Trochulus villosus): confronting competing refugia hypotheses using model selection

The localization of Last Glacial Maximum (LGM) refugia is crucial information to understand a species’ history and predict its reaction to future climate changes. However, many phylogeographical studies often lack sampling designs intensive enough to precisely localize these refugia. The hairy land snail Trochulus villosus has a small range centred on Switzerland, which could be intensively covered by sampling 455 individuals from 52 populations. Based on mitochondrial DNA sequences (COI and 16S), we identified two divergent lineages with distinct geographical distributions. Bayesian skyline plots suggested that both lineages expanded at the end of the LGM. To find where the origin populations were located, we applied the principles of ancestral character reconstruction and identified a candidate refugium for each mtDNA lineage: the French Jura and Central Switzerland, both ice‐free during the LGM. Additional refugia, however, could not be excluded, as suggested by the microsatellite analysis of a population subset. Modelling the LGM niche of T. villosus, we showed that suitable climatic conditions were expected in the inferred refugia, but potentially also in the nunataks of the alpine ice shield. In a model selection approach, we compared several alternative recolonization scenarios by estimating the Akaike information criterion for their respective maximum‐likelihood migration rates. The ‘two refugia’ scenario received by far the best support given the distribution of genetic diversity in T. villosus populations. Provided that fine‐scale sampling designs and various analytical approaches are combined, it is possible to refine our necessary understanding of species responses to environmental changes.

Isolation by Distance in a Population of a Small Land Snail Trochoidea geyeri: Evidence from Direct and Indirect Methods

Population structure was estimated in a continuous population of a small land snail (Trochoidea geyeri). Mark–recapture experiments and randomly amplified polymorphic DNA analyses indicate that the population structure can be described by the isolation by distance model of Wright (1946). Estimates of density and dispersal suggest a neighbourhood size of 70–208 individuals on an area of 13–21 m2. A principal component analysis of the randomly amplified polymorphic DNA data reveals clinal variation of genetic composition across the population, as predicted by the neighbourhood concept. An analysis of molecular variance indicates substantial genetic structuring. Comparisons of the genetic distances, expressed as euclidean distances among individuals, versus the geographic distance between sampling sites yield a highly significant positive correlation (Mantel test: r = 0.567, p < 0.0001). The revealed pattern of populational subdivision on a microgeographic scale seems to be one of the principal processes generating and maintaining genetic diversity within populations of small land gastropods.

Separating historic events from recurrent processes in cryptic species: phylogeography of mud snails (Hydrobia spp.)

The present study combines methods that were designed to infer intraspecific relationships (e.g. nested‐clade analysis (NCA), mismatch distributions and maximum likelihood gene flow analysis) to analyse historic events and recurrent processes in the cryptic mud snail species Hydrobia acuta and H. glyca. Specifically, we test the proposed allopatry of cryptic species and whether the peculiar range‐subdivision of the putative subspecies H. a. acuta and H. a. neglecta is a result of long‐distance dispersal or continuous range expansion. The NCA indicates a past fragmentation of the two H. acuta subspecies as well as past fragmentations within H. glyca. Gene‐flow analyses show extensive gene flow in an E–W direction (towards the Atlantic) in the Mediterranean H. a. acuta, generally low gene flow in a W–E direction in the Atlantic H. a. neglecta and complex gene‐flow pattern in a N–S but also in a S–N direction (against the Gulf Stream) in H. glyca. Based on these data and supportive ecological and oceanographical data, we hypothesize that the separation of the two H. acuta subspecies was not caused by long‐distance dispersal but by a range shift and/or range expansion of the closely related competitor H. glyca as a result of an interglacial warming with a subsequent range shift in H. acuta. Moreover, our data do not show evidence for a long‐term, stable sympatry of Hydrobia species, supporting the concept of allopatric relationships within cryptic radiations. NCA and gene‐flow analyses indicate that the only sympatric population found in our study is the result of a recent dispersal event from the nearby Mediterranean. It is assumed that allopatric relationships in ephemeral Hydrobia populations constitute an evolutionary advantage relative to competition, recruitment and re‐establishment of habitats. Mechanisms that could be of relevance for maintaining allopatry are discussed.