AimTaxon co-occurrence analysis is commonly used in ecology, but it has not been applied to range-wide distribution data of partly allopatric taxa because existing methods cannot differentiate between distribution-related effects and taxon interactions. Our first aim was to develop a taxon co-occurrence analysis method that is also capable of taking into account the effect of species ranges and can handle faunistic records from museum databases or biodiversity inventories. Our second aim was to test the independence of taxon co-occurrences of rock-dwelling gastropods at different taxonomic levels, with a special focus on the Clausiliidae subfamily Alopiinae, and in particular the genus Montenegrina.LocationBalkan Peninsula in south-eastern Europe (46N–36N, 13.5E–28E).MethodsWe introduced a taxon-specific metric that characterizes the occurrence probability at a given location. This probability was calculated as a distance-weighted mean of the taxon’s presence and absence records at all sites. We applied corrections to account for the biases introduced by varying sampling intensity in our dataset. Then we used probabilistic null-models to simulate taxon distributions under the null hypothesis of no taxon interactions and calculated pairwise and cumulated co-occurrences. Independence of taxon occurrences was tested by comparing observed co-occurrences to simulated values.ResultsWe observed significantly fewer co-occurrences among species and intra-generic lineages of Montenegrina than expected under the assumption of no taxon interaction.Main conclusionsFewer than expected co-occurrences among species and intrageneric clades indicate that species divergence preceded niche partitioning. This suggests a primary role of non-adaptive processes in the speciation of rock-dwelling gastropods. The method can account for the effects of distributional constraints in range-wide datasets, making it suitable for testing ecological, biogeographical, or evolutionary hypotheses where interactions of partly allopatric taxa are in question.

No abstract

Delimiting species in radiations is notoriously difficult because of the small differences between the incipient species, the star-like tree with short branches between species, incomplete lineage sorting, and the possibility of introgression between several of the incipient species. Next generation sequencing data may help to overcome some of these problems. We evaluated methods for species delimitation based on genome-wide markers in a land snail radiation on Crete. Species delimitation in the Albinaria cretensis group was based exclusively on shell characters until now and resulted in classifications distinguishing 3–9 species. We generated sequences of 4270 loci for 140 specimens of the Albinaria cretensis group from 48 populations by double-digest restriction site-associated DNA sequencing. We evaluated three methods for species discovery. The multispecies coalescent approach implemented in the program Bayesian Phylogenetics and Phylogeography resulted in a drastic overestimating of the number of species, whereas Gaussian clustering resulted in an overlumping. Primary species hypotheses based on the maximum percentage of the genome of the individuals derived from ancestral populations as estimated with the program ADMIXTURE moderately overestimated the number of species, but this was the only approach that provided information about gene flow between groups. Two of the methods for species validation that we applied, BFD* and delimitR, resulted in an acceptance of almost all primary species hypotheses, even such based on arbitrary subdivisions of hypotheses based on ADMIXTURE. In contrast, secondary species hypotheses, resulting from an evaluation of primary species hypotheses based on ADMIXTURE with isolation by distance tests, approached the morphological classification, but also uncovered two cryptic species and indicated that some of the previously delimited units should be combined. Thus, we recommend this combination of approaches that provided more detailed insights in the distinctness of barriers between the taxa of a species complex and the spatial distribution of admixture between them than the other methods. The recognition and delimitation of undersampled species remained a major challenge.

The Austrian endemic land snail species Noricella oreinos (formerly Trochulus oreinos) occurs in the Northeastern Calcareous Alps at high elevations. Two morphologically highly similar subspecies N. o. oreinos and N. o. scheerpeltzi have been described. First analyses of mitochondrial and nuclear marker sequences indicated a high genetic divergence between them. In the present study, we aimed to assess gene flow between the two subspecies which should help to re‐evaluate their taxonomic status. Sequence data and amplified fragment length polymorphism (AFLP) markers of 255 Noricella specimens covering the whole distribution range were analyzed. A clear geographic separation was found within the potential contact zone, the Haller Mauern mountain range. Samples of all western sites were part of the clade representing N. o. scheerpeltzi and almost all samples from the eastern sites clustered with N. o. oreinos. However, within two sampling sites of the eastern Haller Mauern, a few individuals possessed a COI sequence matching the N. o. oreinos clade whereas at the ITS2 locus they were heterozygous possessing the alleles of both taxa. Contrary to the ITS2 results indicating historical and/or ongoing hybridization, AFLP analyses of 202 individuals confirmed a clear separation of the two taxa congruent with the mitochondrial data. Although they occur on the same mountain range without any physical barrier, no indication of ongoing gene flow between the two taxa was found. Thus, we conclude that the two taxa are separate species N. oreinos and N. scheerpeltzi.