Multispecies coalescent delimits structure, not species

Sequence clustering threshold has little effect on the recovery of microbial community structure

Analysis of microbial community structure by multivariate ordination methods, using data obtained by high-throughput sequencing of amplified markers (i.e., DNA metabarcoding), often requires clustering of DNA sequences into operational taxonomic units (OTUs). Parameters for the clustering procedure tend not to be justified but are set by tradition rather than being based on explicit knowledge. In this study, we explore the extent to which ordination results are affected by variation in parameter settings for the clustering procedure. Amplicon sequence data from nine microbial community studies, representing different sampling designs, spatial scales and ecosystems, were subjected to clustering into OTUs at seven different similarity thresholds (clustering thresholds) ranging from 87% to 99% sequence similarity. The 63 data sets thus obtained were subjected to parallel DCA and GNMDS ordinations. The resulting community structures were highly similar across all clustering thresholds. We explain this pattern by the existence of strong ecological structuring gradients and phylogenetically diverse sets of abundant OTUs that are highly stable across clustering thresholds. Removing low-abundance, rare OTUs had negligible effects on community patterns. Our results indicate that microbial data sets with a clear gradient structure are highly robust to choice of sequence clustering threshold.

CLADES: A classification-based machine learning method for species delimitation from population genetic data

Species are considered to be the basic unit of ecological and evolutionary studies. As multilocus genomic data are increasingly available, there have been considerable interests in the use of DNA sequence data to delimit species. In this study, we show that machine learning can be used for species delimitation. Our method treats the species delimitation problem as a classification problem for identifying the category of a new observation on the basis of training data. Extensive simulation is first conducted over a broad range of evolutionary parameters for training purposes. Each pair of known populations is combined to form training samples with a label of "same species" or "different species". We use support vector machine (SVM) to train a classifier using a set of summary statistics computed from training samples as features. The trained classifier can classify a test sample to two outcomes: "same species" or "different species". Given multilocus genomic data of multiple related organisms or populations, our method (called CLADES) performs species delimitation by first classifying pairs of populations. CLADES then delimits species by maximizing the likelihood of species assignment for multiple populations. CLADES is evaluated through extensive simulation and also tested on real genetic data. We show that CLADES is both accurate and efficient for species delimitation when compared with existing methods. CLADES can be useful especially when existing methods have difficulty in delimitation, for example with short species divergence time and gene flow.

Phylogenomics and barcoding of Panax: toward the identification of ginseng species

BACKGROUND

The economic value of ginseng in the global medicinal plant trade is estimated to be in excess of US$2.1 billion. At the same time, the evolutionary placement of ginseng (Panax ginseng) and the complex evolutionary history of the genus is poorly understood despite several molecular phylogenetic studies. In this study, we use a full plastome phylogenomic framework to resolve relationships in Panax and to identify molecular markers for species discrimination.

RESULTS

We used high-throughput sequencing of MBD2-Fc fractionated Panax DNA to supplement publicly available plastid genomes to create a phylogeny based on fully assembled and annotated plastid genomes from 60 accessions of 8 species. The plastome phylogeny based on a 163 kbp matrix resolves the sister relationship of Panax ginseng with P. quinquefolius. The closely related species P. vietnamensis is supported as sister of P. japonicus. The plastome matrix also shows that the markers trnC-rps16, trnS-trnG, and trnE-trnM could be used for unambiguous molecular identification of all the represented species in the genus.

CONCLUSIONS

MBD2 depletion reduces the cost of plastome sequencing, which makes it a cost-effective alternative to Sanger sequencing based DNA barcoding for molecular identification. The plastome phylogeny provides a robust framework that can be used to study the evolution of morphological characters and biosynthesis pathways of ginsengosides for phylogenetic bioprospecting. Molecular identification of ginseng species is essential for authenticating ginseng in international trade and it provides an incentive for manufacturers to create authentic products with verified ingredients.

Genetic connectivity from the Arctic to the Antarctic: Sclerolinum contortum and Nicomache lokii (Annelida) are both widespread in reducing environments

The paradigm of large geographic ranges in the deep sea has been challenged by genetic studies, which often reveal putatively widespread species to be several taxa with more restricted ranges. Recently, a phylogeographic study revealed that the tubeworm Sclerolinum contortum (Siboglinidae) inhabits vents and seeps from the Arctic to the Antarctic. Here, we further test the conspecificity of the same populations of S. contortum with additional mitochondrial and nuclear markers. We also investigate the genetic connectivity of another species with putatively the same wide geographic range - Nicomache lokii (Maldanidae). Our results support the present range of S. contortum, and the range of N. lokii is extended from vents and seeps in the Nordic Seas to mud volcanoes in the Barbados Trench and Antarctic vents. Sclerolinum contortum shows more pronounced geographic structure than N. lokii, but whether this is due to different dispersal capacities or reflects the geographic isolation of the sampled localities is unclear. Two distinct mitochondrial lineages of N. lokii are present in the Antarctic, which may result from two independent colonization events. The environmental conditions inhabited by the two species and implications for their distinct habitat preference is discussed.

Impact of Model Violations on the Inference of Species Boundaries Under the Multispecies Coalescent

The use of genetic data for identifying species-level lineages across the tree of life has received increasing attention in the field of systematics over the past decade. The multispecies coalescent model provides a framework for understanding the process of lineage divergence and has become widely adopted for delimiting species. However, because these studies lack an explicit assessment of model fit, in many cases, the accuracy of the inferred species boundaries are unknown. This is concerning given the large amount of empirical data and theory that highlight the complexity of the speciation process. Here, we seek to fill this gap by using simulation to characterize the sensitivity of inference under the multispecies coalescent (MSC) to several violations of model assumptions thought to be common in empirical data. We also assess the fit of the MSC model to empirical data in the context of species delimitation. Our results show substantial variation in model fit across data sets. Posterior predictive tests find the poorest model performance in data sets that were hypothesized to be impacted by model violations. We also show that while the inferences assuming the MSC are robust to minor model violations, such inferences can be biased under some biologically plausible scenarios. Taken together, these results suggest that researchers can identify individual data sets in which species delimitation under the MSC is likely to be problematic, thereby highlighting the cases where additional lines of evidence to identify species boundaries are particularly important to collect. Our study supports a growing body of work highlighting the importance of model checking in phylogenetics, and the usefulness of tailoring tests of model fit to assess the reliability of particular inferences. [Populations structure, gene flow, demographic changes, posterior prediction, simulation, genetics.].

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.

Species delimitation and biogeography of the gnatcatchers and gnatwrens (Aves: Polioptilidae)

The New World avian family Polioptilidae (gnatcatchers and gnatwrens) is distributed from Argentina to Canada and includes 15 species and more than 60 subspecies. No study to date has evaluated phylogenetic relationships within this family and the historical pattern of diversification within the group remains unknown. Moreover, species limits, particularly in widespread taxa that show geographic variation, remain unclear. In this study, we delimited species and estimated phylogenetic relationships using multilocus data for the entire family. We then used the inferred diversity along with alternative taxonomic classification schemes to evaluate how lumping and splitting of both taxa and geographical areas influenced biogeographic inference. Species-tree analyses grouped Polioptilidae into four main clades: Microbates, Ramphocaenus, a Polioptila guianensis complex, and the remaining members of Polioptila. Ramphocaenus melanurus was sister to the clade containing M. cinereiventris and M. collaris, which formed a clade sister to all species within Polioptila. Polioptila was composed of two clades, the first of which included the P. guianensis complex; the other contained all remaining species in the genus. Using multispecies coalescent modeling, we inferred a more than 3-fold increase in species diversity, of which 87% represent currently recognized species or subspecies. Much of this diversity corresponded to subspecies that occur in the Neotropics. We identified three polyphyletic species, and delimited 4-6 previously undescribed candidate taxa. Probabilistic modeling of geographic ranges on the species tree indicated that the family likely had an ancestral origin in South America, with all three genera independently colonizing North America. Support for this hypothesis, however, was sensitive to the taxonomic classification scheme used and the number of geographical areas allowed. Our study proposes the first phylogenetic hypothesis for Polioptilidae and provides genealogical support for the reclassification of species limits. Species limits and the resolution of geographical areas that taxa inhabit influence the inferred spatial diversification history.

A biting commentary: Integrating tooth characters with molecular data doubles known species diversity in a lineage of sea slugs that consume "killer algae"

Predicting biotic resistance to highly invasive strains of "killer algae" (Caulerpa spp.) requires understanding the diversity and feeding preferences of native consumers, including sea slugs in family Oxynoidae. Past studies reported low algal host specificity for Oxynoe (6 spp.) and Lobiger (4 spp.), but these taxonomically challenging slugs may represent species complexes of unrecognized specialists that prefer different Caulerpa spp. Here, we assess global diversity of these genera by integrating gene sequences with morphological data from microscopic teeth and internal shells, the only hard parts in these soft-bodied invertebrates. Four delimitation methods applied to datasets comprising mtDNA and/or nuclear alleles yielded up to 16 species hypotheses for samples comprising five nominal taxa, including five highly divergent species in Lobiger and five in Oxynoe. Depending on the analysis, a further four to six species were recovered in the O. antillarum-viridis complex, a clade in which mitochondrial divergence was low and nuclear alleles were shared among lineages. Bayesian species delimitation using only morphological data supported most candidate species, however, and integrative analyses combining morphological and genetic data fully supported all complex members. Collectively, our findings double the recognized biodiversity in Oxynoidae, and illustrate the value of including data from traits that mediate fast-evolving ecological interactions during species delimitation. Preference for Caulerpa spp. and radular tooth characteristics covaried among newly delimited species, highlighting an unappreciated degree of host specialization and coevolution in these taxa that may help predict their role in containing outbreaks of invasive algae.

Revisiting the phylogeography, demography and taxonomy of the frog genus Ptychadena in the Ethiopian highlands with the use of genome-wide SNP data

Understanding the diversification of biological lineages is central to evolutionary studies. To properly study the process of speciation, it is necessary to link micro-evolutionary studies with macro-evolutionary mechanisms. Micro-evolutionary studies require proper sampling across a taxon's range to adequately infer genetic diversity. Here we use the grass frogs of the genus Ptychadena from the Ethiopian highlands as a model to study the process of lineage diversification in this unique biodiversity hotspot. We used thousands of genome-wide SNPs obtained from double digest restriction site associated DNA sequencing (ddRAD-seq) in populations of the Ptychadena neumanni species complex from the Ethiopian highlands in order to infer their phylogenetic relationships and genetic structure, as well as to study their demographic history. Our genome-wide phylogenetic study supports the existence of approximately 13 lineages clustered into 3 species groups. Our phylogenetic and phylogeographic reconstructions suggest that those endemic lineages diversified in allopatry, and subsequently specialized to different habitats and elevations. Demographic analyses point to a continuous decrease in the population size across the majority of lineages and populations during the Pleistocene, which is consistent with a continuous period of aridification that East Africa experienced since the Pliocene. We discuss the taxonomic implications of our analyses and, in particular, we warn against the recent practice to solely use Bayesian species delimitation methods when proposing taxonomic changes.

Rhachotropis (Eusiroidea, Amphipoda) from the North East Atlantic

The genus Rhachotropis has the widest geographic and bathymetric distribution of all amphipod genera worldwide. Molecular and morphological investigations of specimens sampled around Iceland and off the Norwegian coast allow the first insights into the relationships of North East Atlantic Rhachotropis. The 31 cytochrome oxidase subunit I (COI) sequences generated for this study were assigned 13 Barcode Index Numbers (BINs) in the Barcode of Life database (BOLD), of which 12 are new to the database. Molecular analyses of COI and 16S sequences could not confirm a theory that depth has a greater influence on the phylogeny of Rhachotropis than geographic distance. Although the North East Atlantic is a well-studied area, our molecular investigations revealed the genus Rhachotropis may contain cryptic species, which indicates a higher biodiversity than currently known. For example, the specimens which key to Rhachotropis helleri is a complex of three COI clades, two of which cannot be identified with morphological traits. One specimen of each of the clades in the cladogram was documented by high definition photographs. A special focus was on the visual morphology of the eyes, as this character shows interspecific differences within the genus Rhachotropis in response to fixation in ethanol. Detailed morphological investigation showed that some clades thought to be indistinguishable can be separated by minute but consistent morphological characters. Datamining Genbank to examine all registered COI-sequences of R. aculeata, the only previously known Rhachotropis BIN in the North Atlantic and sub-Arctic, showed R. aculeata to be subdivided by an Arctic and a North Atlantic population.

Genome Sequences Reveal Cryptic Speciation in the Human Pathogen Histoplasma capsulatum

Histoplasma capsulatum is a pathogenic fungus that causes life-threatening lung infections. About 500,000 people are exposed to H. capsulatum each year in the United States, and over 60% of the U.S. population has been exposed to the fungus at some point in their life. We performed genome-wide population genetics and phylogenetic analyses with 30 Histoplasma isolates representing four recognized areas where histoplasmosis is endemic and show that the Histoplasma genus is composed of at least four species that are genetically isolated and rarely interbreed. Therefore, we propose a taxonomic rearrangement of the genus.IMPORTANCE The evolutionary processes that give rise to new pathogen lineages are critical to our understanding of how they adapt to new environments and how frequently they exchange genes with each other. The fungal pathogen Histoplasma capsulatum provides opportunities to precisely test hypotheses about the origin of new genetic variation. We find that H. capsulatum is composed of at least four different cryptic species that differ genetically and also in virulence. These results have implications for the epidemiology of histoplasmosis because not all Histoplasma species are equivalent in their geographic range and ability to cause disease.

Multiple markers, niche modelling, and bioregions analyses to evaluate the genetic diversity of a plant species complex

Background

The classification of closely related plants is not straightforward. These morphologically similar taxa frequently maintain their inter-hybridization potential and share ancestral polymorphisms as a consequence of their recent divergence. Under the biological species concept, they may thus not be considered separate species. The Petunia integrifolia complex is especially interesting because, in addition to the features mentioned above, its taxa share a pollinator, and their geographical ranges show multiple overlaps. Here, we combined plastid genome sequences, nuclear microsatellites, AFLP markers, ecological niche modelling, and bioregions analysis to investigate the genetic variability between the different taxa of the P. integrifolia complex in a comprehensive sample covering the entire geographical range of the complex.

Results

Results from molecular markers did not fully align with the current taxonomic classification. Niche modelling and bioregions analyses revealed that taxa were associated with different ecological constraints, indicating that the habitat plays an important role in preserving species boundaries. For three taxa, our analyses showed a mostly conserved, non-overlapping geographical distribution over time. However, for two taxa, niche modelling found an overlapping distribution over time; these taxa were also associated with the same bioregions.

Conclusions

cpDNA markers were better able to discriminate between Petunia taxa than SSRs and AFLPs. Overall, our results suggest that the P. integrifolia complex represents a continuum of individuals from distant and historically isolated populations, which share some morphological traits, but are established in four different evolutionary lineages.

Electronic supplementary material

The online version of this article (10.1186/s12862-017-1084-y) contains supplementary material, which is available to authorized users.

Coalescent-based delimitation outperforms distance-based methods for delineating less divergent species: the case of Kurixalus odontotarsus species group

Few empirical studies have compared coalescent-based methods to distance-based methods for delimitation of less divergent species. In this study, we used two coalescent-based (BFD and BPP) and two distance-based barcoding (ABGD and jMOTU) methods to delimit closely related species in the Kurixalus odontotarsus species group. Phylogenetic analyses revealed that the K. odontotarsus species group comprises 11 distinct maternal clades with strong support values. Based on the genetic and morphological evidences, we consider that species diversity in the K. odontotarsus species group was underestimated and the 11 clades represent 11 species, of which six are unnamed. The coalescent-based delimitations decisively supported the scenario of 11-species corresponding to the 11 clades. However, the distance-based ABGD only obtained 3-6 candidate species, which is not consistent with morphological evidence. These results indicate that BFD and BPP are more conservative than ABGD to false negatives (lumping). Method of fixed threshold (jMOTU) may obtain a resolution similar to that inferred by BFD and BPP, but it severely relies on subjective choice of the threshold and lacks statistical support. We consider that coalescent-based BFD and BPP approaches outperform distance-based methods for delineation of less divergent species.

Maximizing Power in Phylogenetics and Phylogenomics: A Perspective Illuminated by Fungal Big Data

Since its original inception over 150 years ago by Darwin, we have made tremendous progress toward the reconstruction of the Tree of Life. In particular, the transition from analyzing datasets comprised of small numbers of loci to those comprised of hundreds of loci, if not entire genomes, has aided in resolving some of the most vexing of evolutionary problems while giving us a new perspective on biodiversity. Correspondingly, phylogenetic trees have taken a central role in fields that span ecology, conservation, and medicine. However, the rise of big data has also presented phylogenomicists with a new set of challenges to experimental design, quantitative analyses, and computation. The sequencing of a number of very first genomes presented significant challenges to phylogenetic inference, leading fungal phylogenomicists to begin addressing pitfalls and postulating solutions to the issues that arise from genome-scale analyses relevant to any lineage across the Tree of Life. Here we highlight insights from fungal phylogenomics for topics including systematics and species delimitation, ecological and phenotypic diversification, and biogeography while providing an overview of progress made on the reconstruction of the fungal Tree of Life. Finally, we provide a review of considerations to phylogenomic experimental design for robust tree inference. We hope that this special issue of Advances in Genetics not only excites the continued progress of fungal evolutionary biology but also motivates the interdisciplinary development of new theory and methods designed to maximize the power of genomic scale data in phylogenetic analyses.

Do ampharetids take sedimented steps between vents and seeps? Phylogeny and habitat-use of Ampharetidae (Annelida, Terebelliformia) in chemosynthesis-based ecosystems

BACKGROUND

A range of higher animal taxa are shared across various chemosynthesis-based ecosystems (CBEs), which demonstrates the evolutionary link between these habitats, but on a global scale the number of species inhabiting multiple CBEs is low. The factors shaping the distributions and habitat specificity of animals within CBEs are poorly understood, but geographic proximity of habitats, depth and substratum have been suggested as important. Biogeographic studies have indicated that intermediate habitats such as sedimented vents play an important part in the diversification of taxa within CBEs, but this has not been assessed in a phylogenetic framework. Ampharetid annelids are one of the most commonly encountered animal groups in CBEs, making them a good model taxon to study the evolution of habitat use in heterotrophic animals. Here we present a review of the habitat use of ampharetid species in CBEs, and a multi-gene phylogeny of Ampharetidae, with increased taxon sampling compared to previous studies.

RESULTS

The review of microhabitats showed that many ampharetid species have a wide niche in terms of temperature and substratum. Depth may be limiting some species to a certain habitat, and trophic ecology and/or competition are identified as other potentially relevant factors. The phylogeny revealed that ampharetids have adapted into CBEs at least four times independently, with subsequent diversification, and shifts between ecosystems have happened in each of these clades. Evolutionary transitions are found to occur both from seep to vent and vent to seep, and the results indicate a role of sedimented vents in the transition between bare-rock vents and seeps.

CONCLUSION

The high number of ampharetid species recently described from CBEs, and the putative new species included in the present phylogeny, indicates that there is considerable diversity still to be discovered. This study provides a molecular framework for future studies to build upon and identifies some ecological and evolutionary hypotheses to be tested as new data is produced.

Molecular and morphological data reveal non-monophyly and speciation in imperiled freshwater mussels (Anodontoides and Strophitus)

Accurate taxonomic placement is vital to conservation efforts considering many intrinsic biological characteristics of understudied species are inferred from closely related taxa. The rayed creekshell, Anodontoides radiatus (Conrad, 1834), exists in the Gulf of Mexico drainages from western Florida to Louisiana and has been petitioned for listing under the Endangered Species Act. We set out to resolve the evolutionary history of A. radiatus, primarily generic placement and species boundaries, using phylogenetic, morphometric, and geographic information. Our molecular matrix contained 3 loci: cytochrome c oxidase subunit I, NADH dehydrogenase subunit I, and the nuclear-encoded ribosomal internal transcribed spacer I. We employed maximum likelihood and Bayesian inference to estimate a phylogeny and test the monophyly of Anodontoides and Strophitus. We implemented two coalescent-based species delimitation models to test seven species models and evaluate species boundaries within A. radiatus. Concomitant to molecular data, we also employed linear morphometrics and geographic information to further evaluate species boundaries. Molecular and morphological evidence supports the inclusion of A. radiatus in the genus Strophitus, and we resurrect the binomial Strophitus radiatus to reflect their shared common ancestry. We also found strong support for polyphyly in Strophitus and advocate the resurrection of the genus Pseudodontoideus to represent 'Strophitus' connasaugaensis and 'Strophitus' subvexus. Strophitus radiatus exists in six well-supported clades that were distinguished as evolutionary independent lineages using Bayesian inference, maximum likelihood, and coalescent-based species delimitation models. Our integrative approach found evidence for as many as 4 evolutionary divergent clades within S. radiatus. Therefore, we formally describe two new species from the S. radiatus species complex (Strophitus williamsi and Strophitus pascagoulaensis) and recognize the potential for a third putative species (Strophitus sp. cf. pascagoulaensis). Our findings aid stakeholders in establishing conservation and management strategies for the members of Anodontoides, Strophitus, and Pseudodontoideus.

Species boundaries in the human pathogen Paracoccidioides

The use of molecular taxonomy for identifying recently diverged species has transformed the study of speciation in fungi. The pathogenic fungus Paracoccidioides spp has been hypothesized to be composed of five phylogenetic species, four of which compose the brasiliensis species complex. Nuclear gene genealogies support this divergence scenario, but mitochondrial loci do not; while all species from the brasiliensis complex are differentiated at nuclear coding loci, they are not at mitochondrial loci. We addressed the source of this incongruity using 11 previously published gene fragments, 10 newly-sequenced nuclear non-coding loci, and 10 microsatellites. We hypothesized and further demonstrated that the mito-nuclear incongruence in the brasiliensis species complex results from interspecific hybridization and mitochondrial introgression, a common phenomenon in eukaryotes. Additional population genetic analyses revealed possible nuclear introgression but much less than that seen in the mitochondrion. Our results are consistent with a divergence scenario of secondary contact and subsequent mitochondrial introgression despite the continued persistence of species boundaries. We also suggest that yeast morphology slightly-but significantly-differs across all five Paracoccidioides species and propose to elevate four of these phylogenetic species to formally described taxonomic species.

Validation and description of two new north-western Australian Rainbow skinks with multispecies coalescent methods and morphology

While methods for genetic species delimitation have noticeably improved in the last decade, this remains a work in progress. Ideally, model based approaches should be applied and considered jointly with other lines of evidence, primarily morphology and geography, in an integrative taxonomy framework. Deep phylogeographic divergences have been reported for several species of Carlia skinks, but only for some eastern taxa have species boundaries been formally tested. The present study does this and revises the taxonomy for two species from northern Australia, Carlia johnstonei and C. triacantha. We introduce an approach that is based on the recently published method StarBEAST2, which uses multilocus data to explore the support for alternative species delimitation hypotheses using Bayes Factors (BFD). We apply this method, jointly with two other multispecies coalescent methods, using an extensive (from 2,163 exons) data set along with measures of 11 morphological characters. We use this integrated approach to evaluate two new candidate species previously revealed in phylogeographic analyses of rainbow skinks (genus Carlia) in Western Australia. The results based on BFD StarBEAST2, BFD* SNAPP and BPP genetic delimitation, together with morphology, support each of the four recently identified Carlia lineages as separate species. The BFD StarBEAST2 approach yielded results highly congruent with those from BFD* SNAPP and BPP. This supports use of the robust multilocus multispecies coalescent StarBEAST2 method for species delimitation, which does not require a priori resolved species or gene trees. Compared to the situation in C. triacantha, morphological divergence was greater between the two lineages within Kimberley endemic C. johnstonei, which also had deeper divergent histories. This congruence supports recognition of two species within C. johnstonei. Nevertheless, the combined evidence also supports recognition of two taxa within the more widespread C. triacantha. With this work, we describe two new species, Carlia insularis sp. nov and Carlia isostriacantha sp. nov. in the northwest of Australia. This contributes to increasing recognition that this region of tropical Australia has a rich and unique fauna.

The complex evolutionary history and phylogeography of Caridina typus (Crustacea: Decapoda): long-distance dispersal and cryptic allopatric species

The evolutionary history of the old, diverse freshwater shrimp genus Caridina is still poorly understood, despite its vast distribution – from Africa to Polynesia. Here, we used nuclear and mitochondrial DNA to infer the phylogeographic and evolutionary history of C. typus, which is one of only four species distributed across the entire range of the genus. Despite this species’ potential for high levels of gene flow, questions have been raised regarding its phylogeographic structure and taxonomic status. We identified three distinct lineages that likely diverged in the Miocene. Molecular dating and ancestral range reconstructions are congruent with C. typus’ early dispersal to Africa, possibly mediated by the Miocene Indian Ocean Equatorial Jet, followed by back dispersal to Australasia after the Jet’s closure. Furthermore, several different species delimitation methods indicate each lineage represents a distinct (cryptic) species, contradicting current morphospecies delimitation of a single C. typus taxon. The evolutionary history of C. typus lineages is complex, in which ancient oceanic current systems and (currently unrecognised) speciation events preceded secondary sympatry of these cryptic species.

Hybridization patterns in two contact zones of grass snakes reveal a new Central European snake species

Recent studies found major conflicts between traditional taxonomy and genetic differentiation of grass snakes and identified previously unknown secondary contact zones. Until now, little is known about gene flow across these contact zones. Using two mitochondrial markers and 13 microsatellite loci, we examined two contact zones. One, largely corresponding to the Rhine region, involves the western subspecies Natrix natrix helvetica and the eastern subspecies N. n. natrix, whereas in the other, more easterly, contact zone two lineages meet that are currently identified with N. n. natrix and N. n. persa. This second contact zone runs across Central Europe to the southern Balkans. Our analyses reveal that the western contact zone is narrow, with parapatrically distributed mitochondrial lineages and limited, largely unidirectional nuclear gene flow. In contrast, the eastern contact zone is very wide, with massive nuclear admixture and broadly overlapping mitochondrial lineages. In combination with additional lines of evidence (morphology, phylogeny, divergence times), we conclude that these differences reflect different stages in the speciation process and that Natrix helvetica should be regarded as a distinct species. We suggest a nomenclatural framework for presently recognized grass snake taxa and highlight the need for reconciling the conflicts between genetics and taxonomy.

Genomic signatures of paleodrainages in a freshwater fish along the southeastern coast of Brazil: genetic structure reflects past riverine properties

Past shifts in connectivity in riverine environments (for example, sea-level changes) and the properties of current drainages can act as drivers of genetic structure and demographic processes in riverine population of fishes. However, it is unclear whether the same river properties that structure variation on recent timescales will also leave similar genomic signatures that reflect paleodrainage properties. By characterizing genetic structure in a freshwater fish species (Hollandichthys multifasciatus) from a system of basins along the Atlantic coast of Brazil we test for the effects of paleodrainages caused by sea-level changes during the Pleistocene. Given that the paleodrainage properties differ along the Brazilian coast, we also evaluate whether estimated genetic diversity within paleodrainages can be explained by past riverine properties (i.e., area and number of rivers in a paleodrainage). Our results demonstrate that genetic structure between populations is not just highly concordant with paleodrainages, but that differences in the genetic diversity among paleodrainages correspond to the joint effect of differences in the area encompassed by, and the number of rivers, within a paleodrainage. Our findings extend the influence of current riverine properties on genetic diversity to those associated with past paleodrainage properties. We discuss how these findings may explain the inconsistent support for paleodrainages in structuring divergence from different global regions and the importance of taking into account past conditions for understanding the high species diversity of freshwater fish that we currently observe in the world, and especially in the Neotropics.

Integrative species delimitation of the widespread North American jumping mice (Zapodinae)

Delimiting species can be challenging, but is a key step for the critical examination of evolutionary history and for prioritizing conservation efforts. Because systematic relationships are often determined iteratively using tests based on taxonomy, such methods can fail to detect cryptic variation and result in biased conclusions. Conversely, discovery-based approaches provide a powerful way to define operational taxonomic units and test species boundaries. We compare both approaches (taxonomy-based delimitation - TBD and discovery-based delimitation - DBD) within North American jumping mice (Zapodinae) using broad sampling, multilocus analyses, and ecological tests. This group diversified through the dynamic glacial-interglacial periods of the Quaternary and phylogeographic tests reveal 28 lineages that correspond poorly with current taxonomy (4 species, 32 nominal subspecies). However, neither the 4-species or 28-lineage hypotheses are optimal for species-level classification. Rather, information theoretic approaches (Bayes Factors) indicate a 15-species hypothesis is best for characterizing genetic variation in this group, with subsequent iterative pairwise ecological tests failing to confirm four species pairs. Taken together, evolutionary and ecological tests capture divergence among 11 putative species that, if upheld by additional tests, will lead to taxonomic revision and reevaluation of conservation plans.

Revision of the genus Reinmara Schaus, 1928 (Lepidoptera, Mimallonoidea, Mimallonidae) with the descriptions of four new species from South America

The mimallonid genus Reinmara Schaus, 1928 is revised. The three previously described species, R. enthona (Schaus, 1905), R. minasa Schaus, 1928, and R. wolfei Herbin & C. Mielke, 2014 are redescribed and the females of each are described and figured for the first time. Additionally, we describe four new species, two Andean: R. andensissp. n. and R. occidentalissp. n., and two Brazilian: R. atlanticasp. n. and R. igneasp. n.. The new species R. ignea and R. atlantica are likely of conservation concern due to their rarity in collections and their apparent endemism to an endangered biome, the Brazilian Atlantic Forest.