Morphological Species Delimitation in The Western Pond Turtle (Actinemys): Can Machine Learning Methods Aid in Cryptic Species Identification?

As the discovery of cryptic species has increased in frequency, there has been an interest in whether geometric morphometric data can detect fine-scale patterns of variation that can be used to morphologically diagnose such species. We used a combination of geometric morphometric data and an ensemble of five supervised machine learning methods (MLMs) to investigate whether plastron shape can differentiate two putative cryptic turtle species, Actinemys marmorata and Actinemys pallida. Actinemys has been the focus of considerable research due to its biogeographic distribution and conservation status. Despite this work, reliable morphological diagnoses for its two species are still lacking. We validated our approach on two datasets, one consisting of eight morphologically disparate emydid species, the other consisting of two subspecies of Trachemys (T. scripta scripta, T. scripta elegans). The validation tests returned near-perfect classification rates, demonstrating that plastron shape is an effective means for distinguishing taxonomic groups of emydids via MLMs. In contrast, the same methods did not return high classification rates for a set of alternative phylogeographic and morphological binning schemes in Actinemys. All classification hypotheses performed poorly relative to the validation datasets and no single hypothesis was unequivocally supported for Actinemys. Two hypotheses had machine learning performance that was marginally better than our remaining hypotheses. In both cases, those hypotheses favored a two-species split between A. marmorata and A. pallida specimens, lending tentative morphological support to the hypothesis of two Actinemys species. However, the machine learning results also underscore that Actinemys as a whole has lower levels of plastral variation than other turtles within Emydidae, but the reason for this morphological conservatism is unclear.

Coalescent-Based Species Delimitation in Herbaceous Bamboos (Bambusoideae, Olyreae) from Eastern Brazil: Implications for Taxonomy and Conservation in a Group with Weak Morphological Divergence Coupled with Low Genetic Diversity

Species delimitation in herbaceous bamboos has been complex and, in some genera, a great part of its diversity has been confirmed only based on genetic information, as is the case of the genus Raddia. It includes nine species, all occurring in Brazil, but only R. portoi predominates in dry forests of the Northeast associated with the Caatinga phytogeographic domain. This species is morphologically close to R. angustifolia, which is known for a single location in the Atlantic Forest in Southern Bahia, and is considered to be threatened by extinction. Besides problems with taxonomic focus, actions for its conservation are complicated because it is not certain if it must be considered an independent species or included in the more widespread R. portoi. In this study, we used coalescent multispecies (MSC) theory approaches combined with genetic structure analyses in an attempt to delimit these two species. Different analyses were congruent and the species delimitation using MSC inferred distinct lineages supporting their recognition as two species. These results solved the taxonomic doubts and also showed the power of these approaches to delimit species as lineages, even in groups with weak morphological divergence and low genetic variability, and also impacting our knowledge for conservation purposes.

The Complete Chloroplast Genome of Endangered Species Stemona parviflora: Insight into the Phylogenetic Relationship and Conservation Implications

Stemona parviflora is an endangered species, narrowly endemic to Hainan and Southwest Guangdong. The taxonomic classification of S. parviflora remains controversial. Moreover, studying endangered species is helpful for current management and conservation. In this study, the first complete chloroplast genome of S. parviflora was assembled and compared with other Stemona species. The chloroplast genome size of S. parviflora was 154,552 bp, consisting of 87 protein-coding genes, 38 tRNA genes, 8 rRNA genes, and one pseudogene. The ψycf1 gene was lost in the cp genome of S. sessilifolia, but it was detected in four other species of Stemona. The inverted repeats (IR) regions have a relatively lower length variation compared with the large single copy (LSC) and small single copy (SSC) regions. Long repeat sequences and simple sequence repeat (SSR) were detected, and most SSR were distributed in the LSC region. Codon usage bias analyses revealed that the RSCU value of the genus Stemona has almost no difference. As with most angiosperm chloroplast genomes, protein-coding regions were more conservative than the inter-gene spacer. Seven genes (atpI, ccsA, cemA, matK, ndhA, petA, and rpoC1) were detected under positive selection in different Stemona species, which may result from adaptive evolution to different habitats. Phylogenetic analyses show the Stemona cluster in two main groups; S. parviflora were closest to S. tuberosa. A highly suitable region of S. parviflora was simulated by Maxent in this study; it is worth noting that the whole territory of Taiwan has changed to a low fitness area and below in the 2050 s, which may not be suitable for the introduction and cultivation of S. parviflora. In addition, limited by the dispersal capacity of S. parviflora, it is necessary to carry out artificial grafts to expand the survival areas of S. parviflora. Our results provide valuable information on characteristics of the chloroplast genome, phylogenetic relationships, and potential distribution range of the endangered species S. parviflora.

Diversity investigation by application of DNA barcoding: A case study of lepidopteran insects in Xinjiang wild fruit forests, China

To investigate the species diversity of lepidopteran insects in Xinjiang wild fruit forests, establish insect community monitoring systems, and determine the local species pool, we test the applicability of DNA barcoding based on cytochrome c oxidase subunit I (COI) gene for accurate and rapid identification of insect species. From 2017 to 2019, a total of 212 samples with ambiguous morphological identification were selected for DNA barcoding analysis. Five sequence-based methods for species delimitation (ABGD, BINs, GMYC, jMOTU, and bPTP) were conducted for comparison to traditional morphology-based identification. In total, 2,422 samples were recorded, representing 143 species of 110 genera in 17 families in Lepidoptera. The diversity analysis showed that the richness indices for Noctuidae was the highest (54 species), and for Pterophoridae, Cossidae, Limacodidae, Lasiocampidae, Pieridae, and Lycaenidae were the lowest (all with 1 species). The Shannon-Wiener species diversity index (H') and Pielou's evenness (J') of lepidopteran insects first increased and then decreased across these 3 years, while the Simpson diversity index showed a trend of subtracted then added. For molecular-based identification, 67 lepidopteran species within 61 genera in 14 families were identified through DNA barcoding. Neighbor-joining (NJ) analysis showed that conspecific individuals were clustered together and formed monophyletic groups with a high support value, except for Lacanobia contigua (Denis & Schiffermüller, 1775) (Noctuidae: Hadeninae). Sixty-seven morphospecies were classified into various numbers of MOTUs based on ABGD, BINs, GMYC, jMOTU, and bPTP (70, 96, 2, 71, and 71, respectively). In Xinjiang wild fruit forests, the family with the largest number of species is Noctuidae, followed by Geometridae, Crambidae, and the remaining families. The highest Shannon diversity index is observed for the family Noctuidae. Our results indicate that the distance-based methods (ABGD and jMOTU) and character-based method (bPTP) outperform GMYC. BINs is inclined to overestimate species diversity compared to other methods.

An integrative genomic and phenomic analysis to investigate the nature of plant species in Escallonia (Escalloniaceae)

What we mean by species and whether they have any biological reality has been debated since the early days of evolutionary biology. Some biologists even suggest that plant species are created by taxonomists as a subjective, artificial division of nature. However, the nature of plant species has been rarely tested critically with data while ignoring taxonomy. We integrate phenomic and genomic data collected across hundreds of individuals at a continental scale to investigate this question in Escallonia (Escalloniaceae), a group of plants which includes 40 taxonomic species (the species proposed by taxonomists). We first show that taxonomic species may be questionable as they match poorly to patterns of phenotypic and genetic variation displayed by individuals collected in nature. We then use explicit statistical methods for species delimitation designed for phenotypic and genomic data, and show that plant species do exist in Escallonia as an objective, discrete property of nature independent of taxonomy. We show that such species correspond poorly to current taxonomic species (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$< 20\%$$\end{document}<20%) and that phenomic and genomic data seldom delimit congruent entities (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$< 20\%$$\end{document}<20%). These discrepancies suggest that evolutionary forces additional to gene flow can maintain the cohesion of species. We propose that phenomic and genomic data analyzed on an equal footing build a broader perspective on the nature of plant species by helping delineate different ‘types of species’. Our results caution studies which take the accuracy of taxonomic species for granted and challenge the notion of plant species without empirical evidence. Note: A version of the complete manuscript in Spanish is available in the Supplemental Materials.

How to Render Species Comparable Taxonomic Units Through Deep Time: a Case Study on Intraspecific Osteological Variability in Extant and Extinct Lacertid Lizards

Generally, the species is considered to be the only naturally occurring taxon. However, species recognised and defined using different species delimitation criteria cannot readily be compared, impacting studies of biodiversity through Deep Time. This comparability issue is particularly marked when comparing extant with extinct species, because the only available data for species delimitation in fossils is derived from their preserved morphology, which is generally restricted to osteology in vertebrates. Here, we quantify intraspecific, intrageneric, and intergeneric osteological variability in extant species of lacertid lizards using pairwise dissimilarity scores based on a dataset of 253 discrete osteological characters for 99 specimens referred to 24 species. Variability is always significantly lower intraspecifically than between individuals belonging to distinct species of a single genus, which is in turn significantly lower than intergeneric variability. Average values of intraspecific variability and associated standard deviations are consistent (with few exceptions), with an overall average within a species of 0.208 changes per character scored. Application of the same methods to six extinct lacertid species (represented by 40 fossil specimens) revealed that intraspecific osteological variability is inconsistent, which can at least in part be attributed to different researchers having unequal expectations of the skeletal dissimilarity within species units. Such a divergent interpretation of intraspecific and interspecific variability among extant and extinct species reinforces the incomparability of the species unit. Lacertidae is an example where extant species recognised and defined based on a number of delimitation criteria show comparable and consistent intraspecific osteological variability. Here, as well as in equivalent cases, application of those skeletal dissimilarity values to palaeontological species delimitation potentially provides a way to ameliorate inconsistencies created by the use of morphology to define species.

The Exploitation of Distance Distributions for Clustering

Although distance measures are used in many machine learning algorithms, the literature on the context-independent selection and evaluation of distance measures is limited in the sense that prior knowledge is used. In cluster analysis, current studies evaluate the choice of distance measure after applying unsupervised methods based on error probabilities, implicitly setting the goal of reproducing predefined partitions in data. Such studies use clusters of data that are often based on the context of the data as well as the custom goal of the specific study. Depending on the data context, different properties for distance distributions are judged to be relevant for appropriate distance selection. However, if cluster analysis is based on the task of finding similar partitions of data, then the intrapartition distances should be smaller than the interpartition distances. By systematically investigating this specification using distribution analysis through the mirrored-density (MD plot), it is shown that multimodal distance distributions are preferable in cluster analysis. As a consequence, it is advantageous to model distance distributions with Gaussian mixtures prior to the evaluation phase of unsupervised methods. Experiments are performed on several artificial datasets and natural datasets for the task of clustering.

Evaluating Species Delimitation Methods in Radiations: The Land Snail Albinaria cretensis Complex on Crete

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.

Reconstruction of forest dynamics in the Western Palaearctic based on phylogeographic analysis of the ringlet butterfly Erebia aethiops

Glacial refugia are centers of high biodiversity. Therefore, knowledge on their locations and reactions of associated populations and landscapes to climatic changes is crucial for conservation management. We here investigated the biogeography of a butterfly species linked to open forest habitats. Using mitochondrial and nuclear markers in combination with Bayesian simulations, we analyzed the location and age of potential glacial refugia of the species. We identified five putative refugia in Europe. Considering the ecological needs of our study species, tree density within these refugial areas, in contrast to earlier assumptions, must have exceeded the level of individually scattered trees. Our results also provide evidence that especially the refuge areas in the Carpathians were previously underestimated regarding their age: the refugia in the Southern Carpathians presented suitable conditions throughout several glacial cycles, probably since the Mindel or Riss cycles. Additionally, our analyses provided support for a forest refugium near the Tatra Mountains persisting the last glacial maximum. Our results underline the usefulness of this and probably other butterfly species as indicators of forest refugia.

Evolutionary systematics of the viviparous gastropod Sermyla (Gastropoda: Cerithioidea: Thiaridae), with the description of a new species

While most Cerithioidea are marine, some occur in brackish and freshwater habitats. Many members are systematically problematic due to variability or homoplasy in conchological characters, which has led to taxonomic redundancy, but also because of discrepancies between phylogenetic trees and morphologically distinguishable units as revealed in recent molecular genetic studies. We have chosen an evolutionary systematic approach and combine analyses of shell biometry and geometric morphometrics with the analyses of reproductive traits and molecular genetics based on mtDNA and AFLP markers in order to resolve the relationships among species of the genus Sermyla. We describe a new species from Sulawesi, Sermyla kupaensis sp. nov., which is characterized by a distinct reproductive strategy. This unique reproductive strategy corresponds with its distinct molecular genetic signal. However, it is not possible to distinguish S. kupaensis from S. riquetii based on shell morphology alone. We also provide data on the population structure of the endemic Australian species Sermyla carbonata, for which we found a drainage-based population structure. Overall, we present a new concept of the relationships among the species within the genus Sermyla based on morphological and genetic data.

Clustering benchmark datasets exploiting the fundamental clustering problems

The Fundamental Clustering Problems Suite (FCPS) offers a variety of clustering challenges that any algorithm should be able to handle given real-world data. The FCPS consists of datasets with known a priori classifications that are to be reproduced by the algorithm. The datasets are intentionally created to be visualized in two or three dimensions under the hypothesis that objects can be grouped unambiguously by the human eye. Each dataset represents a certain problem that can be solved by known clustering algorithms with varying success. In the R package “Fundamental Clustering Problems Suite” on CRAN, user-defined sample sizes can be drawn for the FCPS. Additionally, the distances of two high-dimensional datasets called Leukemia and Tetragonula are provided here. This collection is useful for investigating the shortcomings of clustering algorithms and the limitations of dimensionality reduction methods in the case of three-dimensional or higher datasets. This article is a simultaneous co-submission with Swarm Intelligence for Self-Organized Clustering [1].

Diversification and independent polyploid origins in the disjunct species Alyssum repens from the Southeastern Alps and the Carpathians

PREMISE

Disjunct distributions have been commonly observed in mountain plant species and have stimulated phylogeographic and phylogenetic research. Here we studied Alyssum repens, a member of the polyploid species complex A. montanum-A. repens, which exhibits SE Alpine-Carpathian disjunctions with a large elevational span and consists of diploid and tetraploid populations. We aimed to investigate the species' genetic and cytotype structure in the context of its distribution patterns, to elucidate the polyploid origins and to propose an appropriate taxonomic treatment.

METHODS

We combined AFLP fingerprinting markers, sequence variation of the highly repetitive ITS region of rDNA and the low-copy DET1 nuclear gene, genome size, and morphometric data.

RESULTS

We identified four geographically structured genetic lineages. One consisted of diploid populations from the foothills of the Southeastern Alps and neighboring regions, and the three others were allopatric montane to alpine groups comprising diploids and tetraploids growing in the Southeastern Carpathians and the Apuseni Mts. in Romania.

CONCLUSIONS

We inferred a vicariance scenario associated with Quaternary climatic oscillations, accompanied by one auto- and two allopolyploidization events most likely involving a northern Balkan relative. Whereas genetic differentiation and allopatric distribution would favor the taxonomic splitting of this species, the genetic lineages largely lack morphological distinguishability, and their ecological, cytotype and genome size divergence is only partial. Even though we probably face here a case of incipient speciation, we propose to maintain the current taxonomic treatment of Alyssum repens as a single, albeit variable, species.

Integrating multilocus DNA data and 3D geometric morphometrics to elucidate species boundaries in the case of Pyrenaearia (Pulmonata: Hygromiidae)

To accurately delimit species the use of multiple character types is essential as all speciation processes are not equally reflected in different data (e.g. morphological, molecular or ecological characters). With the introduction of geometric morphometrics methods and advances in 3D technology, a comprehensive combination of molecular and morphological data has been enabled in groups where exhaustively quantifying and measuring morphological shape change was not possible before such as gastropod shells. In this study, we combined multilocus coalescent species delimitation methods with 3D geometric morphometrics of shell shape to delimit species within the land snail genus Pyrenaearia. A new taxonomic scheme was constructed for the genus identifying ten species. Two nominal species were synonymized and a hitherto unrecognized cryptic species was identified. Our findings support the importance of combining multiple lines of evidence as molecular and morphological data on their own do not yield the same information. Further, the integration of morphological and molecular data shows the importance of allometry in shell shape and suggests a combined effect of population history and selection in different environments on shells morphological variation. Our new taxonomy and phylogenetic reconstruction suggest that, besides the glacial cycles of the Pleistocene, passive dispersal and rock substrate complexity could also have been involved in the speciation of the genus.

Issues and Perspectives in Species Delimitation using Phenotypic Data: Atlantean Evolution in Darwin's Finches

Progress in the development and use of methods for species delimitation employing phenotypic data lags behind conceptual and practical advances in molecular genetic approaches. The basic evolutionary model underlying the use of phenotypic data to delimit species assumes random mating and quantitative polygenic traits, so that phenotypic distributions within a species should be approximately normal for individuals of the same sex and age. Accordingly, two or more distinct normal distributions of phenotypic traits suggest the existence of multiple species. In light of this model, we show that analytical approaches employed in taxonomic studies using phenotypic data are often compromised by three issues: 1) reliance on graphical analyses that convey little information on phenotype frequencies; 2) exclusion of characters potentially important for species delimitation following reduction of data dimensionality; and 3) use of measures of central tendency to evaluate phenotypic distinctiveness. We outline approaches to overcome these issues based on statistical developments related to normal mixture models (NMMs) and illustrate them empirically with a reanalysis of morphological data recently used to claim that there are no morphologically distinct species of Darwin's ground-finches (Geospiza). We found negligible support for this claim relative to taxonomic hypotheses recognizing multiple species. Although species limits among ground-finches merit further assessments using additional sources of information, our results bear implications for other areas of inquiry including speciation research: because ground-finches have likely speciated and are not trapped in a process of "Sisyphean" evolution as recently argued, they remain useful models to understand the evolutionary forces involved in speciation. Our work underscores the importance of statistical approaches grounded on appropriate evolutionary models for species delimitation. We discuss how NMMs offer new perspectives in the kind of inferences available to systematists, with significant repercussions on ideas about the phenotypic structure of biodiversity.

Identifying cryptic diversity with predictive phylogeography

Identifying units of biological diversity is a major goal of organismal biology. An increasing literature has focused on the importance of cryptic diversity, defined as the presence of deeply diverged lineages within a single species. While most discoveries of cryptic lineages proceed on a taxon-by-taxon basis, rapid assessments of biodiversity are needed to inform conservation policy and decision-making. Here, we introduce a predictive framework for phylogeography that allows rapidly identifying cryptic diversity. Our approach proceeds by collecting environmental, taxonomic and genetic data from codistributed taxa with known phylogeographic histories. We define these taxa as a reference set, and categorize them as either harbouring or lacking cryptic diversity. We then build a random forest classifier that allows us to predict which other taxa endemic to the same biome are likely to contain cryptic diversity. We apply this framework to data from two sets of disjunct ecosystems known to harbour taxa with cryptic diversity: the mesic temperate forests of the Pacific Northwest of North America and the arid lands of Southwestern North America. The predictive approach presented here is accurate, with prediction accuracies placed between 65% and 98.79% depending of the ecosystem. This seems to indicate that our method can be successfully used to address ecosystem-level questions about cryptic diversity. Further, our application for the prediction of the cryptic/non-cryptic nature of unknown species is easily applicable and provides results that agree with recent discoveries from those systems. Our results demonstrate that the transition of phylogeography from a descriptive to a predictive discipline is possible and effective.

A novel method to infer the origin of polyploids from Amplified Fragment Length Polymorphism data reveals that the alpine polyploid complex of Senecio carniolicus (Asteraceae) evolved mainly via autopolyploidy

Despite its evolutionary and ecological relevance, the mode of polyploid origin has been notoriously difficult to be reconstructed from molecular data. Here, we present a method to identify the putative parents of polyploids and thus to infer the mode of their origin (auto- vs. allopolyploidy) from Amplified Fragment Length Polymorphism (AFLP) data. To this end, we use Cohen's d of distances between in silico polyploids, generated within a priori defined scenarios of origin from a priori delimited putative parental entities (e.g. taxa, genetic lineages), and natural polyploids. Simulations show that the discriminatory power of the proposed method increases mainly with increasing divergence between the lower-ploid putative ancestors and less so with increasing delay of polyploidization relative to the time of divergence. We apply the new method to the Senecio carniolicus aggregate, distributed in the European Alps and comprising two diploid, one tetraploid and one hexaploid species. In the eastern part of its distribution, the S. carniolicus aggregate was inferred to comprise an autopolyploid series, whereas for western populations of the tetraploid species, an allopolyploid origin involving the two diploid species was the most likely scenario. Although this suggests that the tetraploid species has two independent origins, other evidence (ribotype distribution, morphology) is consistent with the hypothesis of an autopolyploid origin with subsequent introgression by the second diploid species. Altogether, identifying the best among alternative scenarios using Cohen's d can be straightforward, but particular scenarios, such as allopolyploid origin vs. autopolyploid origin with subsequent introgression, remain difficult to be distinguished.

Comparative Phylogeography of Ethiopian anurans: impact of the Great Rift Valley and Pleistocene climate change

Background

The Ethiopian highlands are a biodiversity hotspot, split by the Great Rift Valley into two distinct systems of plateaus and mountains. The Rift Valley is currently hot and dry and acts as a barrier to gene flow for highland-adapted species. It is however unlikely that the conditions in the Rift were inhospitable to highland species during the entire Pleistocene. To assess the significance of the Ethiopian Rift as a biogeographic barrier as well as the impact Pleistocene climatic changes have had on the evolution of Ethiopian organisms, we performed phylogeographic analyses and developed present and past niche models on seven anuran species with different elevational and ecological preferences.

Results

We found that highland species on the east and the west sides of the Rift are genetically differentiated and have not experienced any detectable gene flow for at least 0.4 my. In contrast, species found at elevations lower than 2500 m do not show any population structure. We also determined that highland species have lower effective population sizes than lowland species, which have experienced a large, yet gradual, demographic expansion, starting approximately half a million year ago.

Conclusions

The pattern we report here is consistent with the increasingly warmer and drier conditions of the Pleistocene in East Africa, which resulted in the expansion of savanna, the fragmentation of forests and the shrinking of highland habitats. Climatic niche models indicated that the Rift is currently non suitable for most of the studied species, but it could have been a more permeable barrier during the Last Glacial Maximum. However, considering the strong genetic structure of highland species, we hypothesize that the barrier mechanisms at the Rift are not only climatic but also topographical.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0774-1) contains supplementary material, which is available to authorized users.

Population structure and connectivity in the Mediterranean sponge Ircinia fasciculata are affected by mass mortalities and hybridization

Recent episodes of mass mortalities in the Mediterranean Sea have been reported for the closely related marine sponges Ircinia fasciculata and Ircinia variabilis that live in sympatry. In this context, the assessment of the genetic diversity, bottlenecks and connectivity of these sponges has become urgent in order to evaluate the potential effects of mass mortalities on their latitudinal range. Our study aims to establish (1) the genetic structure, connectivity and signs of bottlenecks across the populations of I. fasciculata and (2) the hybridization levels between I. fasciculata and I. variabilis. To accomplish the first objective, 194 individuals of I. fasciculata from 12 locations across the Mediterranean were genotyped at 14 microsatellite loci. For the second objective, mitochondrial cytochrome c oxidase subunit I sequences of 16 individuals from both species were analyzed along with genotypes at 12 microsatellite loci of 40 individuals coexisting in 3 Mediterranean populations. We detected strong genetic structure along the Mediterranean for I. fasciculata, with high levels of inbreeding in all locations and bottleneck signs in most locations. Oceanographic barriers like the Almeria-Oran front, North-Balearic front and the Ligurian-Thyrrenian barrier seem to be impeding gene flow for I. fasciculata, adding population divergence to the pattern of isolation by distance derived from the low dispersal abilities of sponge larvae. Hybridization between both species occurred in some populations that might be increasing genetic diversity and somewhat palliating the genetic loss caused by population decimation in I. fasciculata.

Investigation of population structure in Gulf of Mexico Seepiophila jonesi (Polychaeta, Siboglinidae) using cross-amplified microsatellite loci

BACKGROUND

Vestimentiferan tubeworms are some of the most recognizable fauna found at deep-sea cold seeps, isolated environments where hydrocarbon rich fluids fuel biological communities. Several studies have investigated tubeworm population structure; however, much is still unknown about larval dispersal patterns at Gulf of Mexico (GoM) seeps. As such, researchers have applied microsatellite markers as a measure for documenting the transport of vestimentiferan individuals. In the present study, we investigate the utility of microsatellites to be cross-amplified within the escarpiid clade of seep vestimentiferans, by determining if loci originally developed for Escarpia spp. could be amplified in the GoM seep tubeworm, Seepiophila jonesi. Additionally, we determine if cross-amplified loci can reliably uncover the same signatures of high gene flow seen in a previous investigation of S. jonesi.

METHODS

Seventy-seven S. jonesi individuals were collected from eight seep sites across the upper Louisiana slope (<1,000 m) in the GoM. Forty-eight microsatellite loci that were originally developed for Escarpia laminata (18 loci) and Escarpia southwardae (30 loci) were tested to determine if they were homologous and polymorphic in S. jonesi. Loci found to be both polymorphic and of high quality were used to test for significant population structuring in S. jonesi.

RESULTS

Microsatellite pre-screening identified 13 (27%) of the Escarpia loci were homologous and polymorphic in S. jonesi, revealing that microsatellites can be amplified within the escarpiid clade of vestimentiferans. Our findings uncovered low levels of heterozygosity and a lack of genetic differentiation amongst S. jonesi from various sites and regions, in line with previous investigations that employed species-specific polymorphic loci on S. jonesi individuals retrieved from both the same and different seep sites. The lack of genetic structure identified from these populations supports the presence of significant gene flow via larval dispersal in mixed oceanic currents.

DISCUSSION

The ability to develop "universal" microsatellites reduces the costs associated with these analyses and allows researchers to track and investigate a wider array of taxa, which is particularly useful for organisms living at inaccessible locations such as the deep sea. Our study highlights that non-species specific microsatellites can be amplified across large evolutionary distances and still yield similar findings as species-specific loci. Further, these results show that S. jonesi collected from various localities in the GoM represents a single panmictic population, suggesting that dispersal of lecithotrophic larvae by deep sea currents is sufficient to homogenize populations. These data are consistent with the high levels of gene flow seen in Escarpia spp., which advocates that differences in microhabitats of seep localities lead to variation in biogeography of separate species.

Bayesian species delimitation in Pleophylla chafers (Coleoptera) - the importance of prior choice and morphology

BACKGROUND

Defining species units can be challenging, especially during the earliest stages of speciation, when phylogenetic inference and delimitation methods may be compromised by incomplete lineage sorting (ILS) or secondary gene flow. Integrative approaches to taxonomy, which combine molecular and morphological evidence, have the potential to be valuable in such cases. In this study we investigated the South African scarab beetle genus Pleophylla using data collected from 110 individuals of eight putative morphospecies. The dataset included four molecular markers (cox1, 16S, rrnL, ITS1) and morphometric data based on male genital morphology. We applied a suite of molecular and morphological approaches to species delimitation, and implemented a novel Bayesian approach in the software iBPP, which enables continuous morphological trait and molecular data to be combined.

RESULTS

Traditional morphology-based species assignments were supported quantitatively by morphometric analyses of the male genitalia (eigenshape analysis, CVA, LDA). While the ITS1-based delineation was also broadly congruent with the morphospecies, the cox1 data resulted in over-splitting (GMYC modelling, haplotype networks, PTP, ABGD). In the most extreme case morphospecies shared identical haplotypes, which may be attributable to ILS based on statistical tests performed using the software JML. We found the strongest support for putative morphospecies based on phylogenetic evidence using the combined approach implemented in iBPP. However, support for putative species was sensitive to the use of alternative guide trees and alternative combinations of priors on the population size (θ) and rootage (τ 0 ) parameters, especially when the analysis was based on molecular or morphological data alone.

CONCLUSIONS

We demonstrate that continuous morphological trait data can be extremely valuable in assessing competing hypotheses to species delimitation. In particular, we show that the inclusion of morphological data in an integrative Bayesian framework can improve the resolution of inferred species units. However, we also demonstrate that this approach is extremely sensitive to guide tree and prior parameter choice. These parameters should be chosen with caution - if possible - based on independent empirical evidence, or careful sensitivity analyses should be performed to assess the robustness of results. Young species provide exemplars for investigating the mechanisms of speciation and for assessing the performance of tools used to delimit species on the basis of molecular and/or morphological evidence.

Rarity and Incomplete Sampling in DNA-Based Species Delimitation

DNA-based species delimitation may be compromised by limited sampling effort and species rarity, including "singleton" representatives of species, which hampers estimates of intra- versus interspecies evolutionary processes. In a case study of southern African chafers (beetles in the family Scarabaeidae), many species and subclades were poorly represented and 48.5% of species were singletons. Using cox1 sequences from >500 specimens and ∼100 species, the Generalized Mixed Yule Coalescent (GMYC) analysis as well as various other approaches for DNA-based species delimitation (Automatic Barcode Gap Discovery (ABGD), Poisson tree processes (PTP), Species Identifier, Statistical Parsimony), frequently produced poor results if analyzing a narrow target group only, but the performance improved when several subclades were combined. Hence, low sampling may be compensated for by "clade addition" of lineages outside of the focal group. Similar findings were obtained in reanalysis of published data sets of taxonomically poorly known species assemblages of insects from Madagascar. The low performance of undersampled trees is not due to high proportions of singletons per se, as shown in simulations (with 13%, 40% and 52% singletons). However, the GMYC method was highly sensitive to variable effective population size ([Formula: see text]), which was exacerbated by variable species abundances in the simulations. Hence, low sampling success and rarity of species affect the power of the GMYC method only if they reflect great differences in [Formula: see text] among species. Potential negative effects of skewed species abundances and prevalence of singletons are ultimately an issue about the variation in [Formula: see text] and the degree to which this is correlated with the census population size and sampling success. Clade addition beyond a limited study group can overcome poor sampling for the GMYC method in particular under variable [Formula: see text] This effect was less pronounced for methods of species delimitation not based on coalescent models.

Evolution along the Great Rift Valley: phenotypic and genetic differentiation of East African white-eyes (Aves, Zosteropidae)

The moist and cool cloud forests of East Africa represent a network of isolated habitats that are separated by dry and warm lowland savannah, offering an opportunity to investigate how strikingly different selective regimes affect species diversification. Here, we used the passerine genus Zosterops (white-eyes) from this region as our model system. Species of the genus occur in contrasting distribution settings, with geographical mountain isolation driving diversification, and savannah interconnectivity preventing differentiation. We analyze (1) patterns of phenotypic and genetic differentiation in high- and lowland species (different distribution settings), (2) investigate the potential effects of natural selection and temporal and spatial isolation (evolutionary drivers), and (3) critically review the taxonomy of this species complex. We found strong phenotypic and genetic differentiation among and within the three focal species, both in the highland species complex and in the lowland taxa. Altitude was a stronger predictor of phenotypic patterns than the current taxonomic classification. We found longitudinal and latitudinal phenotypic gradients for all three species. Furthermore, wing length and body weight were significantly correlated with altitude and habitat type in the highland species Z. poliogaster. Genetic and phenotypic divergence showed contrasting inter- and intraspecific structures. We suggest that the evolution of phenotypic characters is mainly driven by natural selection due to differences in the two macro-habitats, cloud forest and savannah. In contrast, patterns of neutral genetic variation appear to be rather driven by geographical isolation of the respective mountain massifs. Populations of the Z. poliogaster complex, as well as Z. senegalensis and Z. abyssinicus, are not monophyletic based on microsatellite data and have higher levels of intraspecific differentiation compared to the currently accepted species.

Next-generation sampling: Pairing genomics with herbarium specimens provides species-level signal in Solidago (Asteraceae)

PREMISE OF THE STUDY

The ability to conduct species delimitation and phylogeny reconstruction with genomic data sets obtained exclusively from herbarium specimens would rapidly enhance our knowledge of large, taxonomically contentious plant genera. In this study, the utility of genotyping by sequencing is assessed in the notoriously difficult genus Solidago (Asteraceae) by attempting to obtain an informative single-nucleotide polymorphism data set from a set of specimens collected between 1970 and 2010.

METHODS

Reduced representation libraries were prepared and Illumina-sequenced from 95 Solidago herbarium specimen DNAs, and resulting reads were processed with the nonreference Universal Network-Enabled Analysis Kit (UNEAK) pipeline. Multidimensional clustering was used to assess the correspondence between genetic groups and morphologically defined species.

RESULTS

Library construction and sequencing were successful in 93 of 95 samples. The UNEAK pipeline identified 8470 single-nucleotide polymorphisms, and a filtered data set was analyzed for each of three Solidago subsections. Although results varied, clustering identified genomic groups that often corresponded to currently recognized species or groups of closely related species.

DISCUSSION

These results suggest that genotyping by sequencing is broadly applicable to DNAs obtained from herbarium specimens. The data obtained and their biological signal suggest that pairing genomics with large-scale herbarium sampling is a promising strategy in species-rich plant groups.

Polyploidisation and geographic differentiation drive diversification in a European High Mountain Plant Group (Doronicum clusii Aggregate, Asteraceae)

Range shifts (especially during the Pleistocene), polyploidisation and hybridization are major factors affecting high-mountain biodiversity. A good system to study their role in the European high mountains is the Doronicum clusii aggregate (Asteraceae), whose four taxa (D. clusii s.s., D. stiriacum, D. glaciale subsp. glaciale and D. glaciale subsp. calcareum) are differentiated geographically, ecologically (basiphilous versus silicicolous) and/or via their ploidy levels (diploid versus tetraploid). Here, we use DNA sequences (three plastid and one nuclear spacer) and AFLP fingerprinting data generated for 58 populations to infer phylogenetic relationships, origin of polyploids-whose ploidy level was confirmed by chromosomally calibrated DNA ploidy level estimates-and phylogeographic history. Taxonomic conclusions were informed, among others, by a Gaussian clustering method for species delimitation using dominant multilocus data. Based on molecular data we identified three lineages: (i) silicicolous diploid D. clusii s.s. in the Alps, (ii) silicicolous tetraploid D. stiriacum in the eastern Alps (outside the range of D. clusii s.s.) and the Carpathians and (iii) the basiphilous diploids D. glaciale subsp. glaciale (eastern Alps) and D. glaciale subsp. calcareum (northeastern Alps); each taxon was identified as distinct by the Gaussian clustering, but the separation of D. glaciale subsp. calcareum and D. glaciale subsp. glaciale was not stable, supporting their taxonomic treatment as subspecies. Carpathian and Alpine populations of D. stiriacum were genetically differentiated suggesting phases of vicariance, probably during the Pleistocene. The origin (autopolyploid versus allopolyploid) of D. stiriacum remained unclear. Doronicum glaciale subsp. calcareum was genetically and morphologically weakly separated from D. glaciale subsp. glaciale but exhibited significantly higher genetic diversity and rarity. This suggests that the more widespread D. glaciale subsp. glaciale originated from D. glaciale subsp. calcareum, which is restricted to a prominent Pleistocene refugium previously identified in other alpine plant species.

Population signatures of large-scale, long-term disjunction and small-scale, short-term habitat fragmentation in an Afromontane forest bird

The Eastern Afromontane cloud forests occur as geographically distinct mountain exclaves. The conditions of these forests range from large to small and from fairly intact to strongly degraded. For this study, we sampled individuals of the forest bird species, the Montane White-eye Zosterops poliogaster from 16 sites and four mountain archipelagos. We analysed 12 polymorphic microsatellites and three phenotypic traits, and calculated Species Distribution Models (SDMs) to project past distributions and predict potential future range shifts under a scenario of climate warming. We found well-supported genetic and morphologic clusters corresponding to the mountain ranges where populations were sampled, with 43% of all alleles being restricted to single mountains. Our data suggest that large-scale and long-term geographic isolation on mountain islands caused genetically and morphologically distinct population clusters in Z. poliogaster. However, major genetic and biometric splits were not correlated to the geographic distances among populations. This heterogeneous pattern can be explained by past climatic shifts, as highlighted by our SDM projections. Anthropogenically fragmented populations showed lower genetic diversity and a lower mean body mass, possibly in response to suboptimal habitat conditions. On the basis of these findings and the results from our SDM analysis we predict further loss of genotypic and phenotypic uniqueness in the wake of climate change, due to the contraction of the species' climatic niche and subsequent decline in population size.

The genetic integrity of the ex situ population of the European wildcat (Felis silvestris silvestris) is seriously threatened by introgression from domestic cats (Felis silvestris catus)

Studies on the genetic diversity and relatedness of zoo populations are crucial for implementing successful breeding programmes. The European wildcat, Felis s. silvestris, is subject to intensive conservation measures, including captive breeding and reintroduction. We here present the first systematic genetic analysis of the captive population of Felis s. silvestris in comparison with a natural wild population. We used microsatellites and mtDNA sequencing to assess genetic diversity, structure and integrity of the ex situ population. Our results show that the ex situ population of the European wildcat is highly structured and that it has a higher genetic diversity than the studied wild population. Some genetic clusters matched the breeding lines of certain zoos or groups of zoos that often exchanged individuals. Two mitochondrial haplotype groups were detected in the in situ populations, one of which was closely related to the most common haplotype found in domestic cats, suggesting past introgression in the wild. Although native haplotypes were also found in the captive population, the majority (68%) of captive individuals shared a common mtDNA haplotype with the domestic cat (Felis s. catus). Only six captive individuals (7.7%) were assigned as wildcats in the STRUCTURE analysis (at K = 2), two of which had domestic cat mtDNA haplotypes and only two captive individuals were assigned as purebred wildcats by NewHybrids. These results suggest that the high genetic diversity of the captive population has been caused by admixture with domestic cats. Therefore, the captive population cannot be recommended for further breeding and reintroduction.

Effects of recent and past climatic shifts on the genetic structure of the high mountain yellow-spotted ringlet butterfly Erebia manto (Lepidoptera, Satyrinae): a conservation problem

Mountain species have evolved important genetic differentiation due to past climatic fluctuations. The genetic uniqueness of many of these lineages is now at risk due to global warming. Here, we analyse allozyme polymorphisms of 1306 individuals (36 populations) of the mountain butterfly Erebia manto and perform Species Distribution Models (SDMs). As a consensus of analyses, we obtained six most likely genetic clusters: (i) Pyrenees with Massif Central; (ii) Vosges; (iii-v) Alps including the Slovakian Carpathians; (vi) southern Carpathians. The Vosges population showed the strongest genetic split from all other populations, being almost as strong as the split between E. manto and its sister species Erebia eriphyle. The distinctiveness of the Pyrenees-Massif Central group and of the southern Carpathians group from all other groups is also quite high. All three groups are assumed to have survived more than one full glacial-interglacial cycle close to their current distributions with up-hill and down-slope shifts conforming climatic conditions. In contrast with these well-differentiated groups, the three groups present in the Alps and the Slovakian Carpathians show a much shallower genetic structure and thus also should be of a more recent origin. As predicted by our SDM projections, rising temperatures will strongly impact the distribution of E. manto. While the populations in the Alps are predicted to shrink, the survival of the three lineages present here should not be at risk. The situation of the three other lineages is quite different. All models predict the extinction of the Vosges lineage in the wake of global warming, and also the southern Carpathians and Pyrenees-Massif Central lineages might be at high risk to disappear. Thus, albeit global warming will therefore be unlikely to threaten E. manto as a species, an important proportion of the species' intraspecific differentiation and thus uniqueness might be lost.

Integrating coalescent and phylogenetic approaches to delimit species in the lichen photobiont Trebouxia

The accurate assessment of species boundaries in symbiotic systems is a prerequisite for the study of speciation, co-evolution and selectivity. Many studies have shown the high genetic diversity of green algae from the genus Trebouxia, the most common photobiont of lichen-forming fungi. However, the phylogenetic relationships, and the amount of cryptic diversity of these algae are still poorly understood, and an adequate species concept for trebouxiophycean algae is still missing. In this study we used a multifaceted approach based on coalescence (GMYC, STEM) and phylogenetic relationships to assess species boundaries in the trebouxioid photobionts of the lichen-forming fungus Lasallia pustulata. We further investigated whether putative species of Trebouxia found in L. pustulata are shared with other lichen-forming fungi. We found that L. pustulata is associated with at least five species of Trebouxia and most of them are shared with other lichen-forming fungi, showing different patterns of species-to-species and species-to-community interactions. We also show that one of the putative Trebouxia species is found exclusively in association with L. pustulata and is restricted to thalli from localities with Mediterranean microclimate. We suggest that the species delimitation method presented in this study is a promising tool to address species boundaries within the heterogeneous genus Trebouxia.

Species detection and individual assignment in species delimitation: can integrative data increase efficacy?

Statistical species delimitation usually relies on singular data, primarily genetic, for detecting putative species and individual assignment to putative species. Given the variety of speciation mechanisms, singular data may not adequately represent the genetic, morphological and ecological diversity relevant to species delimitation. We describe a methodological framework combining multivariate and clustering techniques that uses genetic, morphological and ecological data to detect and assign individuals to putative species. Our approach recovers a similar number of species recognized using traditional, qualitative taxonomic approaches that are not detected when using purely genetic methods. Furthermore, our approach detects groupings that traditional, qualitative taxonomic approaches do not. This empirical test suggests that our approach to detecting and assigning individuals to putative species could be useful in species delimitation despite varying levels of differentiation across genetic, phenotypic and ecological axes. This work highlights a critical, and often overlooked, aspect of the process of statistical species delimitation-species detection and individual assignment. Irrespective of the species delimitation approach used, all downstream processing relies on how individuals are initially assigned, and the practices and statistical issues surrounding individual assignment warrant careful consideration.