Species delimitation using a combined coalescent and information-theoretic approach: an example from North American Myotis bats

Coalescent model-based methods for phylogeny estimation force systematists to confront issues related to the identification of species boundaries. Unlike conventional phylogenetic analysis, where species membership can be assessed qualitatively after the phylogeny is estimated, the phylogenies that are estimated under a coalescent model treat aggregates of individuals as the operational taxonomic units and thus require a priori definition of these sets because the models assume that the alleles in a given lineage are sampled from a single panmictic population. Fortunately, the use of coalescent model-based approaches allows systematists to conduct probabilistic tests of species limits by calculating the probability of competing models of lineage composition. Here, we conduct the first exploration of the issues related to applying such tests to a complex empirical system. Sequence data from multiple loci were used to assess species limits and phylogeny in a clade of North American Myotis bats. After estimating gene trees at each locus, the likelihood of models representing all hierarchical permutations of lineage composition was calculated and Akaike information criterion scores were computed. Metrics borrowed from information theory suggest that there is strong support for several models that include multiple evolutionary lineages within the currently described species Myotis lucifugus and M. evotis. Although these results are preliminary, they illustrate the practical importance of coupled species delimitation and phylogeny estimation.

Species delineation using Bayesian model-based assignment tests: a case study using Chinese toad-headed agamas (genus Phrynocephalus)

Background

Species are fundamental units in biology, yet much debate exists surrounding how we should delineate species in nature. Species discovery now requires the use of separate, corroborating datasets to quantify independently evolving lineages and test species criteria. However, the complexity of the speciation process has ushered in a need to infuse studies with new tools capable of aiding in species delineation. We suggest that model-based assignment tests are one such tool. This method circumvents constraints with traditional population genetic analyses and provides a novel means of describing cryptic and complex diversity in natural systems. Using toad-headed agamas of the Phrynocephalus vlangalii complex as a case study, we apply model-based assignment tests to microsatellite DNA data to test whether P. putjatia, a controversial species that closely resembles P. vlangalii morphologically, represents a valid species. Mitochondrial DNA and geographic data are also included to corroborate the assignment test results.

Results

Assignment tests revealed two distinct nuclear DNA clusters with 95% (230/243) of the individuals being assigned to one of the clusters with > 90% probability. The nuclear genomes of the two clusters remained distinct in sympatry, particularly at three syntopic sites, suggesting the existence of reproductive isolation between the identified clusters. In addition, a mitochondrial ND2 gene tree revealed two deeply diverged clades, which were largely congruent with the two nuclear DNA clusters, with a few exceptions. Historical mitochondrial introgression events between the two groups might explain the disagreement between the mitochondrial and nuclear DNA data. The nuclear DNA clusters and mitochondrial clades corresponded nicely to the hypothesized distributions of P. vlangalii and P. putjatia.

Conclusions

These results demonstrate that assignment tests based on microsatellite DNA data can be powerful tools for distinguishing closely related species and support the validity of P. putjatia. Assignment tests have the potential to play a significant role in elucidating biodiversity in the era of DNA data. Nonetheless, important limitations do exist and multiple independent datasets should be used to corroborate results from assignment programs.

New insights into diversity and evolution of deep-sea Mytilidae (Mollusca: Bivalvia)

Bathymodiolinae mussels have been used as a biological model to better understand the evolutionary origin of faunas associated with deep-sea hydrothermal vents and cold seeps. Most studies to date, however, have sampled with a strong bias towards vent and seep species, mainly because of a lack of knowledge of closely related species from organic falls. Here we reassess the species diversity of deep-sea mussels using two genes and a large taxon sample from the South-Western Pacific. This new taxonomic framework serves as a basis for a phylogenetic investigation of their evolutionary history. We first highlight an unexpected allopatric pattern and suggest that mussels usually reported from organic falls are in fact poorly specialized with regard to their environment. This challenges the adaptive scenarios proposed to explain the diversification of the group. Second, we confirm that deep-sea mussels arose from organic falls and then colonized hydrothermal vents and cold seeps in multiple events. Overall, this study constitutes a new basis for further phylogenetic investigations and a global systematic revision of deep-sea mussels.

Bayesian species delimitation in West African forest geckos (Hemidactylus fasciatus)

Genealogical data are an important source of evidence for delimiting species, yet few statistical methods are available for calculating the probabilities associated with different species delimitations. Bayesian species delimitation uses reversible-jump Markov chain Monte Carlo (rjMCMC) in conjunction with a user-specified guide tree to estimate the posterior distribution for species delimitation models containing different numbers of species. We apply Bayesian species delimitation to investigate the speciation history of forest geckos (Hemidactylus fasciatus) from tropical West Africa using five nuclear loci (and mtDNA) for 51 specimens representing 10 populations. We find that species diversity in H. fasciatus is currently underestimated, and describe three new species to reflect the most conservative estimate for the number of species in this complex. We examine the impact of the guide tree, and the prior distributions on ancestral population sizes () and root age (τ(0)), on the posterior probabilities for species delimitation. Mis-specification of the guide tree or the prior distribution for can result in strong support for models containing more species. We describe a new statistic for summarizing the posterior distribution of species delimitation models, called speciation probabilities, which summarize the posterior support for each speciation event on the starting guide tree.

Phylogeny and species boundaries of the endemic species complex, Eliurus antsingy and E. carletoni (Rodentia: Muroidea: Nesomyidae), in Madagascar using mitochondrial and nuclear DNA sequence data

Current understanding of the taxonomy and species relationships within the genus Eliurus, a member of the endemic Malagasy rodent subfamily Nesomyinae, is largely based on morphological characters. Little has been published about the historical biogeography and diversification factors in members of this genus. Here we examined, using molecular data, the evolutionary relationships, and species boundaries within a species complex of Eliurus and the role of rivers in their evolution. Specifically, we tested the hypothesized monophyly of E. antsingy (west) and the recently described E. carletoni (north) using two mitochondrial genes, cytochrome b, the control region, and a nuclear marker, the Interstitial Retinol Binding Protein (IRBP). We also assessed the species boundaries of E. carletoni and examined phylogeographic patterns using two mitochondrial markers and an extensive geographic sampling. Molecular data show that members of the E. antsingy and E. carletoni complex together form a monophyletic group. Although samples provisionally assigned to E. carletoni based on morphological characters from different sites in northern Madagascar cluster in three separate phylogroups, they are monophyletic and assignable to this species. Rivers do not act as strict barriers to dispersal in E. carletoni. The demonstrated monophyly of both E. antsingy and E. carletoni contributes to the amounting evidence for a biogeographic link between the north and western regions of the island.

Integrative taxonomy: a multisource approach to exploring biodiversity

Good alpha taxonomy is central to biology. On the basis of a survey of arthropod studies that used multiple disciplines for species delimitation, we evaluated the performance of single disciplines. All included disciplines had a considerable failure rate. Rigor in species delimitation can thus be increased when several disciplines chosen for complementarity are used. We present a flexible procedure and stopping rule for integrative taxonomy that uses the information from different disciplines separately. Disagreement among disciplines over the number and demarcation of species is resolved by elucidating and invoking evolutionary explanations for disagreement. With the identification of further promising study organisms and of new questions for in-depth analysis, evolutionary biology should profit from integrative taxonomy. An important rationale is clarity in researcher bias in the decision-making process. The success of integrative taxonomy will further increase through methodological progress, taxonomic training of evolutionary biologists, and balanced resource allocation.

Species delimitation: a case study in a problematic ant taxon

Species delimitation has been invigorated as a discipline in systematics by an influx of new character sets, analytical methods, and conceptual advances. We use genetic data from 68 markers, combined with distributional, bioclimatic, and coloration information, to hypothesize boundaries of evolutionarily independent lineages (species) within the widespread and highly variable nominal fire ant species Solenopsis saevissima, a member of a species group containing invasive pests as well as species that are models for ecological and evolutionary research. Our integrated approach uses diverse methods of analysis to sequentially test whether populations meet specific operational criteria (contingent properties) for candidacy as morphologically cryptic species, including genetic clustering, monophyly, reproductive isolation, and occupation of distinctive niche space. We hypothesize that nominal S. saevissima comprises at least 4-6 previously unrecognized species, including several pairs whose parapatric distributions implicate the development of intrinsic premating or postmating barriers to gene flow. Our genetic data further suggest that regional genetic differentiation in S. saevissima has been influenced by hybridization with other nominal species occurring in sympatry or parapatry, including the quite distantly related Solenopsis geminata. The results of this study illustrate the importance of employing different classes of genetic data (coding and noncoding regions and nuclear and mitochondrial DNA [mtDNA] markers), different methods of genetic data analysis (tree-based and non-tree based methods), and different sources of data (genetic, morphological, and ecological data) to explicitly test various operational criteria for species boundaries in clades of recently diverged lineages, while warning against over reliance on any single data type (e.g., mtDNA sequence variation) when drawing inferences.

Delimiting species boundaries within the Neotropical bamboo Otatea (Poaceae: Bambusoideae) using molecular, morphological and ecological data

Species delimitation is a task that has engaged taxonomists for more than two centuries. Recently, it has been demonstrated that molecular data and ecological niche modeling are useful in species delimitation. In this paper multiple data sets (molecular, morphological, ecological) were utilized to set limits for the species belonging to the Neotropical bamboo Otatea, because there is disagreement about species circumscriptions and also because the genus has an interesting distribution, with most of its populations in Mexico and a single disjunct population in Colombia. Molecular and morphological phylogenetic analyses recovered trees with conflicting topologies. Tree-based morphological and character-based analyses recognized the same entities. Ecological niche models and PCA/MANOVAS agreed with the recognition of the same entities that resulted from the morphological analyses. Morphological analyses retrieved clades supported by diagnostic characters and coherent geographical distributions. Based on these results seven entities should be recognized in Otatea, instead of the three previously described species.

Frequent mitochondrial gene introgression among high elevation Tibetan megophryid frogs revealed by conflicting gene genealogies

Historical mitochondrial introgression causes differences between a species' mitochondrial gene genealogy and its nuclear gene genealogy, making tree-based species delineation ambiguous. Using sequence data from one mitochondrial gene (cytochrome b) and three nuclear genes (introns), we examined the evolutionary history of four high elevation Tibetan megophryid frog species, Scutiger boulengeri, Scutiger glandulatus, Scutiger mammatus and Scutiger tuberculatus. The three nuclear genes shared a similar history but the mitochondrial gene tree suggested a drastically different evolutionary scenario. The conflicts between them were explained by multiple episodes of mitochondrial introgression events via historical interspecific hybridization. 'Foreign' mitochondrial genomes might have been fixed in populations and extended through a large portion of the species' distribution. Some hybridization events were probably as old as 10 Myr, while others were recent. An F(1) hybrid was also identified. Historical hybridization events among the four species appeared to be persistent and were not restricted to the period of Pleistocene glaciation, as in several other well-studied cases. Furthermore, hybridization involved several species and occurred in multiple directions, and there was no indication of one mitochondrial genome being superior to others. In addition, incomplete lineage sorting resulting from budding speciation may have also explained some discrepancies between the mitochondrial DNA and nuclear gene trees. Combining all evidences, the former 'Scutiger mammatus' appeared to be two species, including a new species. With the availability of a wide range of highly variable nuclear gene markers, we recommend using a combination of mitochondrial gene and multiple nuclear genes to reveal a complete species history.

Vast underestimation of Madagascar's biodiversity evidenced by an integrative amphibian inventory

Amphibians are in decline worldwide. However, their patterns of diversity, especially in the tropics, are not well understood, mainly because of incomplete information on taxonomy and distribution. We assess morphological, bioacoustic, and genetic variation of Madagascar's amphibians, one of the first near-complete taxon samplings from a biodiversity hotspot. Based on DNA sequences of 2,850 specimens sampled from over 170 localities, our analyses reveal an extreme proportion of amphibian diversity, projecting an almost 2-fold increase in species numbers from the currently described 244 species to a minimum of 373 and up to 465. This diversity is widespread geographically and across most major phylogenetic lineages except in a few previously well-studied genera, and is not restricted to morphologically cryptic clades. We classify the genealogical lineages in confirmed and unconfirmed candidate species or deeply divergent conspecific lineages based on concordance of genetic divergences with other characters. This integrative approach may be widely applicable to improve estimates of organismal diversity. Our results suggest that in Madagascar the spatial pattern of amphibian richness and endemism must be revisited, and current habitat destruction may be affecting more species than previously thought, in amphibians as well as in other animal groups. This case study suggests that worldwide tropical amphibian diversity is probably underestimated at an unprecedented level and stresses the need for integrated taxonomic surveys as a basis for prioritizing conservation efforts within biodiversity hotspots.

Mitochondrial phylogeography of baboons (Papio spp.): indication for introgressive hybridization?

Background

Baboons of the genus Papio are distributed over wide ranges of Africa and even colonized parts of the Arabian Peninsula. Traditionally, five phenotypically distinct species are recognized, but recent molecular studies were not able to resolve their phylogenetic relationships. Moreover, these studies revealed para- and polyphyletic (hereafter paraphyletic) mitochondrial clades for baboons from eastern Africa, and it was hypothesized that introgressive hybridization might have contributed substantially to their evolutionary history. To further elucidate the phylogenetic relationships among baboons, we extended earlier studies by analysing the complete mitochondrial cytochrome b gene and the 'Brown region' from 67 specimens collected at 53 sites, which represent all species and which cover most of the baboons' range.

Results

Based on phylogenetic tree reconstructions seven well supported major haplogroups were detected, which reflect geographic populations and discordance between mitochondrial phylogeny and baboon morphology. Our divergence age estimates indicate an initial separation into southern and northern baboon clades 2.09 (1.54–2.71) million years ago (mya). We found deep divergences between haplogroups within several species (~2 mya, northern and southern yellow baboons, western and eastern olive baboons and northern and southern chacma baboons), but also recent divergence ages among species (< 0.7 mya, yellow, olive and hamadryas baboons in eastern Africa).

Conclusion

Our study confirms earlier findings for eastern Africa, but shows that baboon species from other parts of the continent are also mitochondrially paraphyletic. The phylogenetic patterns suggest a complex evolutionary history with multiple phases of isolation and reconnection of populations. Most likely all these biogeographic events were triggered by multiple cycles of expansion and retreat of savannah biomes during Pleistocene glacial and inter-glacial periods. During contact phases of populations reticulate events (i.e. introgressive hybridization) were highly likely, similar to ongoing hybridization, which is observed between East African baboon populations. Defining the extent of the introgressive hybridization will require further molecular studies that incorporate additional sampling sites and nuclear loci.

Cryptic diversity in vertebrates: molecular data double estimates of species diversity in a radiation of Australian lizards (Diplodactylus, Gekkota)

A major problem for biodiversity conservation and management is that a significant portion of species diversity remains undocumented (the 'taxonomic impediment'). This problem is widely acknowledged to be dire among invertebrates and in developing countries; here, we demonstrate that it can be acute even in conspicuous animals (reptiles) and in developed nations (Australia). A survey of mtDNA, allozyme and chromosomal variation in the Australian gecko, genus Diplodactylus, increases overall species diversity estimates from 13 to 29. Four nominal species each actually represent multi-species complexes; three of these species complexes are not even monophyletic. The high proportion of cryptic species discovered emphasizes the importance of continuing detailed assessments of species diversity, even in apparently well-known taxa from industrialized countries.

Molecular and morphological characterization of heterorhabditid entomopathogenic nematodes from the tropical rainforest in Brazil

Despite massive losses of primary forest, the Amazonian rainforest remains an extremely rich source of biodiversity. In recent years, entomopathogenic nematodes (EPNs) have been isolated from soil in various parts of the world and used successfully as biological control agents against numerous insect pests. Therefore, a sampling in the rainforest of Monte Negro, Rondônia, Brazil was conducted with the aim of discovering new strains and/or species of EPNs for future development as biological control agents. From 156 soil samples taken at nine collecting sites, 19 isolates were obtained, all of them belonging to the genus Heterorhabditis. Four strains were subjected to detailed morphological and molecular evaluation. Based on morphometrics and internal transcribed spacer (ITS) sequence data, the strains LPP1, LPP2 and LPP4 were identified as Heterorhabditis indica, whereas LPP7 was considered Heterorhabditis baujardi. Comparative analysis of the ITS1 sequence of H. indica and H. baujardi isolates showed a polymorphic site for the restriction enzyme Tth 111 that could be used to distinguish the two species. Consequently, strains LPP1, LPP2, LPP3, LPP4, and LPP9 were identified as H. indica, whereas LPP5, LPP7, LPP8 and LPP10 were identified as H. baujardi.

Speciation and DNA barcodes: testing the effects of dispersal on the formation of discrete sequence clusters

Large-scale sequencing of short mtDNA fragments for biodiversity inventories ('DNA barcoding') indicates that sequence variation in animal mtDNA is highly structured and partitioned into discrete genetic clusters that correspond broadly to species-level entities. Here we explore how the migration rate, an important demographic parameter that is directly related to population isolation, might affect variation in the strength of mtDNA clustering among taxa. Patterns of mtDNA variation were investigated in two groups of beetles that both contain lineages occupying habitats predicted to select for different dispersal abilities: predacious diving beetles (Dytiscidae) in the genus Bidessus from lotic and lentic habitats across Europe and darkling beetles (Tenebrionidae) in the genus Eutagenia from sand and other soil types in the Aegean Islands. The degree of genetic clustering was determined using the recently developed 'mixed Yule coalescent' (MYC) model that detects the transition from between-species to within-population branching patterns. Lineages from presumed stable habitats, and therefore displaying lower dispersal ability and migration rates, showed greater levels of mtDNA clustering and geographical subdivision than their close relatives inhabiting ephemeral habitats. Simulations of expected patterns of mtDNA variation under island models showed that MYC clusters are only detected when the migration rates are much lower than the value of Nm=1 typically used to define the threshold for neutral genetic divergence. Therefore, discrete mtDNA clusters provide strong evidence for independently evolving populations or species, but their formation is suppressed even under very low levels of dispersal.

The role of linkage disequilibrium in the evolution of premating isolation

The suggestion that speciation may often occur, or be completed, in the presence of gene flow has long been contentious, due to an appreciation of the challenges to maintaining population- or species-specific gene combinations when gene flow is occurring. Linkage disequilibrium between loci involved in postzygotic and premating isolation must often be built and maintained as the source of these species-specific genotypes. Here, I discuss proposed solutions to facilitate the establishment and maintenance of this linkage disequilibrium. I concentrate primarily on two such factors: one-allele versus two-allele mechanisms of premating isolation, and the form of selection against hybrids as it relates to its effect on the pathway between postzygotic and prezygotic isolation. The goal of this discussion is not to thoroughly review these factors, but instead to concentrate on aspects and implications of these solutions that are currently underemphasized in the speciation literature.

A genealogical approach to quantifying lineage divergence

We introduce a statistic, the genealogical sorting index (gsi), for quantifying the degree of exclusive ancestry of labeled groups on a rooted genealogy and demonstrate its application. The statistic is simple, intuitive, and easily calculated. It has a normalized range to facilitate comparisons among different groups, trees, or studies and it provides information on individual groups rather than a composite measure for all groups. It naturally handles polytomies and accommodates measures of uncertainty in phylogenetic relationships. We use coalescent simulations to explore the behavior of the gsi across a range of divergence times, with the mean value increasing to 1, the maximum value when exclusivity within a group reached monophyly. Simulations also demonstrate that the power to reject the null hypothesis of mixed genealogical ancestry increased markedly as sample size increased, and that the gsi provides a statistically more powerful measure of divergence than FST. Applications to data from published studies demonstrated that the gsi provides a useful way to detect significant exclusivity even when groups are not monophyletic. Although we describe this statistic in the context of divergence, it is more broadly applicable to quantify and assess the significance of clustering of observations in labeled groups on any tree.

Patterns and processes underlying evolutionary significant units in the Platypleura stridula L. species complex (Hemiptera: Cicadidae) in the Cape Floristic Region, South Africa

Cicadas have been shown to be useful organisms for examining the effects of distribution, plant association and geographical barriers on gene flow between populations. The cicadas of the Platypleura stridula species complex are restricted to the biologically diverse Cape Floristic Region (CFR) of South Africa. They are thus an excellent study group for elucidating the mechanisms by which hemipteran diversity is generated and maintained in the CFR. Phylogeographical analysis of this species complex using mitochondrial DNA Cytochrome Oxidase I (COI) and ribosomal 16S sequence data, coupled with preliminary morphological and acoustic data, resolves six clades, each of which has specific host-plant associations and distinct geographical ranges. The phylogeographical structure implies simultaneous or near-simultaneous radiation events, coupled with shifts in host-plant associations. When calibrated using published COI and 16S substitution rates typical for related insects, these lineages date back to the late Pliocene - early Pleistocene, coincident with vegetation change, altered drainage patterns and accelerated erosion in response to neotectonic crustal uplift and cyclic Pleistocene climate change, and glaciation-associated changes in climate and sea level.