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

Artificial intelligence enables unified analysis of historical and landscape influences on genetic diversity

While genetic variation in any species is potentially shaped by a range of processes, phylogeography and landscape genetics are largely concerned with inferring how environmental conditions and landscape features impact neutral intraspecific diversity. However, even as both disciplines have come to utilize SNP data over the last decades, analytical approaches have remained for the most part focused on either broad-scale inferences of historical processes (phylogeography) or on more localized inferences about environmental and/or landscape features (landscape genetics). Here we demonstrate that an artificial intelligence model-based analytical framework can consider both deeper historical factors and landscape-level processes in an integrated analysis. We implement this framework using data collected from two Brazilian anurans, the Brazilian sibilator frog (Leptodactylus troglodytes) and granular toad (Rhinella granulosa). Our results indicate that historical demographic processes shape most the genetic variation in the sibulator frog, while landscape processes primarily influence variation in the granular toad. The machine learning framework used here allows both historical and landscape processes to be considered equally, rather than requiring researchers to make an a priori decision about which factors are important.

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.

Likelihood-Based Tests of Species Tree Hypotheses

Likelihood-based tests of phylogenetic trees are a foundation of modern systematics. Over the past decade, an enormous wealth and diversity of model-based approaches have been developed for phylogenetic inference of both gene trees and species trees. However, while many techniques exist for conducting formal likelihood-based tests of gene trees, such frameworks are comparatively underdeveloped and underutilized for testing species tree hypotheses. To date, widely used tests of tree topology are designed to assess the fit of classical models of molecular sequence data and individual gene trees and thus are not readily applicable to the problem of species tree inference. To address this issue, we derive several analogous likelihood-based approaches for testing topologies using modern species tree models and heuristic algorithms that use gene tree topologies as input for maximum likelihood estimation under the multispecies coalescent. For the purpose of comparing support for species trees, these tests leverage the statistical procedures of their original gene tree-based counterparts that have an extended history for testing phylogenetic hypotheses at a single locus. We discuss and demonstrate a number of applications, limitations, and important considerations of these tests using simulated and empirical phylogenomic data sets that include both bifurcating topologies and reticulate network models of species relationships. Finally, we introduce the open-source R package SpeciesTopoTestR (SpeciesTopology Tests in R) that includes a suite of functions for conducting formal likelihood-based tests of species topologies given a set of input gene tree topologies.

MtDNA species-level phylogeny and delimitation support significantly underestimated diversity and endemism in the largest Neotropical cichlid genus (Cichlidae: Crenicichla)

Crenicichla is the largest and most widely distributed genus of Neotropical cichlids. Here, we analyze a mtDNA dataset comprising 681 specimens (including Teleocichla, a putative ingroup of Crenicichla) and 77 out of 105 presently recognized valid species (plus 10 out of 36 nominal synonyms plus over 50 putatively new species) from 129 locations in 31 major river drainages throughout the whole distribution of the genus in South America. Based on these data we make an inventory of diversity and highlight taxa and biogeographic areas worthy of further sampling effort and conservation protection. Using three methods of molecular species delimitation, we find between 126 and 168 species-like clusters, i.e., an average increase of species diversity of 65–121% with a range of increase between species groups. The increase ranges from 0% in the Missioneira and Macrophthama groups, through 25–40% (Lacustris group), 50–87% (Reticulata group, Teleocichla), 68–168% (Saxatilis group), 125–200% (Wallacii group), and 158–241% in the Lugubris group. We found a high degree of congruence between clusters derived from the three used methods of species delimitation. Overall, our results recognize substantially underestimated diversity in Crenicichla including Teleocichla. Most of the newly delimited putative species are from the Amazon-Orinoco-Guiana (AOG) core area (Greater Amazonia) of the Neotropical region, especially from the Brazilian and Guiana shield areas of which the former is under the largest threat and largest degree of environmental degradation of all the Amazon.

Reply to the comment by Morales et al. on “Population genetics reveal Myotis keenii (Keen’s myotis) and Myotis evotis (long-eared myotis) to be a single species”

A.E. Morales et al. (2021. Can. J. Zool. 99(5): 415–422) provided no new evidence to alter the conclusions of C.L. Lausen et al. (2019. Can. J. Zool. 97(3): 267–279). We present background information, relevant comparisons, and clarification of analyses to further strengthen our conclusions. The genesis of the original “evotis–keenii” study in British Columbia (Canada) was to differentiate Myotis keenii (Merriam, 1895) (Keen’s myotis), with one of the smallest North American bat distributions, from sympatric Myotis evotis (H. Allen, 1864) (long-eared myotis), using something other than the suggested post-mortem skull size comparison, but no differentiating trait could be found, leading to the molecular genetics examination of C.L. Lausen et al. (2019). We present cumulative data that rejects the 1979 hypothesis of M. keenii as a distinct species. A.E. Morales et al. (2021) inaccurately portray C.L. Lausen et al.’s (2019) question and results; present inaccurate morphological and outdated distribution data; overstate the impact of homoplasy without supporting evidence; and misinterpret evidence of population structure.

Comment on “Population genetics reveal Myotis keenii (Keen’s myotis) and Myotis evotis (long-eared myotis) to be a single species”

Genetic exchange and hybridization appear common among the western long-eared bats from North America. Multiple sources of evidence indicate that lineages within this group are evolving independently, despite genetic exchange. However, evidence of gene flow raises questions about the species-level status of some lineages. C.L. Lausen et al. (2019. Can. J. Zool. 97(3): 267–279) proposed that Myotis evotis (H. Allen, 1864) (long-eared myotis) and Myotis keenii (Merriam, 1895) (Keen’s myotis) are one species, not two. This conclusion is based on analyses of cytochrome b and microsatellite data suggesting gene flow between these taxa. Microsatellites are not reliable markers for identifying species because homoplasy can be a major confounding factor, which appears to be true in this case. We reanalyzed the dataset of C.L. Lausen et al. (2019) and show that it is not reliable to distinguish between gene flow or homoplasy, and that these data do not support the conclusion that M. evotis and M. keenii represent a single species. Previous morphological and genomic studies indicate that these are separate species despite previous genetic exchange between them. Failing to recognize that gene flow can occur between independently evolving lineages is counterproductive for conservation because it can lead to neglect of important independent lineages, and likewise failing to use proper tools to delimit species is counterproductive to efforts to quantify biodiversity and design conservation strategies.

Phylogeographic analysis and species distribution modelling of the wood frog Batrachyla leptopus (Batrachylidae) reveal interglacial diversification in south western Patagonia

Background

The evolutionary history of southern South American organisms has been strongly influenced by Pleistocene climate oscillations. Amphibians are good models to evaluate hypotheses about the influence of these climate cycles on population structure and diversification of the biota, because they are sensitive to environmental changes and have restricted dispersal capabilities. We test hypotheses regarding putative forest refugia and expansion events associated with past climatic changes in the wood frog Batrachyla leptopus distributed along ∼1,000 km of length including glaciated and non-glaciated areas in southwestern Patagonia.

Methods

Using three mitochondrial regions (D-loop, cyt b, and coI) and two nuclear loci (pomc and crybA1), we conducted multilocus phylogeographic analyses and species distribution modelling to gain insights of the evolutionary history of this species. Intraspecific genealogy was explored with maximum likelihood, Bayesian, and phylogenetic network approaches. Diversification time was assessed using molecular clock models in a Bayesian framework, and demographic scenarios were evaluated using approximate Bayesian computation (ABC) and extended Bayesian skyline plot (EBSP). Species distribution models (SDM) were reconstructed using climatic and geographic data.

Results

Population structure and genealogical analyses support the existence of four lineages distributed north to south, with moderate to high phylogenetic support (Bootstrap > 70%; BPP > 0.92). The diversification time of B. leptopus’ populations began at ∼0.107 mya. The divergence between A and B lineages would have occurred by the late Pleistocene, approximately 0.068 mya, and divergence between C and D lineages was approximately 0.065 mya. The ABC simulations indicate that lineages coalesced at two different time periods, suggesting the presence of at least two glacial refugia and a postglacial colonization route that may have generated two southern lineages (p = 0.93, type I error: <0.094, type II error: 0.134). EBSP, mismatch distribution and neutrality indexes suggest sudden population expansion at ∼0.02 mya for all lineages. SDM infers fragmented distributions of B. leptopus associated with Pleistocene glaciations. Although the present populations of B. leptopus are found in zones affected by the last glacial maximum (∼0.023 mya), our analyses recover an older history of interglacial diversification (0.107–0.019 mya). In addition, we hypothesize two glacial refugia and three interglacial colonization routes, one of which gave rise to two expanding lineages in the south.

Cryptic Patterns of Speciation in Cryptic Primates: Microendemic Mouse Lemurs and the Multispecies Coalescent

Mouse lemurs (Microcebus) are a radiation of morphologically cryptic primates distributed throughout Madagascar for which the number of recognized species has exploded in the past two decades. This taxonomic revision has prompted understandable concern that there has been substantial oversplitting in the mouse lemur clade. Here, we investigate mouse lemur diversity in a region in northeastern Madagascar with high levels of microendemism and predicted habitat loss. We analyzed RADseq data with multispecies coalescent (MSC) species delimitation methods for two pairs of sister lineages that include three named species and an undescribed lineage previously identified to have divergent mtDNA. Marked differences in effective population sizes, levels of gene flow, patterns of isolation-by-distance, and species delimitation results were found among the two pairs of lineages. Whereas all tests support the recognition of the presently undescribed lineage as a separate species, the species-level distinction of two previously described species, M. mittermeieri and M. lehilahytsara is not supported-a result that is particularly striking when using the genealogical discordance index (gdi). Nonsister lineages occur sympatrically in two of the localities sampled for this study, despite an estimated divergence time of less than 1 Ma. This suggests rapid evolution of reproductive isolation in the focal lineages and in the mouse lemur clade generally. The divergence time estimates reported here are based on the MSC calibrated with pedigree-based mutation rates and are considerably more recent than previously published fossil-calibrated relaxed-clock estimates. We discuss the possible explanations for this discrepancy, noting that there are theoretical justifications for preferring the MSC estimates in this case. [Cryptic species; effective population size; microendemism; multispecies coalescent; speciation; species delimitation.].

Molecular clocks, biogeography and species diversity in Herichthys with evaluation of the role of Punta del Morro as a vicariant brake along the Mexican Transition Zone in the context of local and global time frame of cichlid diversification

Using molecular dated phylogenies and biogeographic reconstructions, the species diversity, biogeography and time frame of evolution of the genus Herichthys were evaluated. In particular, we test the role of Punta del Morro (PdM) as a vicariant brake along the Mexican Transition Zone in the context of local and global time frame of cichlid diversification using several sets of calibrations. Species diversity in Herichthys is complex and the here employed dating methods suggest young age and rapid divergence for many species while species delimitation methods did not resolve these young species including both sympatric species pairs. Based on our molecular clock dating analyses, Herichthys has colonized its present distribution area significantly prior to the suggested vicariance by PdM (10-17.1 Ma vs. 5 to 7.5 Ma). The PdM constraint is in conflict with all other paleogeographic and fossil constraints including novel ones introduced in this study that are, however, congruent among each other. Our study demonstrates that any cichlid datings significantly older or younger than the bounds presented by our analyses and discussion have to be taken as highly questionable from the point of view of Middle American paleogeography and cichlid biogeography unless we allow the option that cichlid biogeography is completely independent from ecological and geological constraints.

Coalescent-based delimitation and species-tree estimations reveal Appalachian origin and Neogene diversification in Russula subsection Roseinae

Numerous lineages of mushroom-forming fungi have been subject to bursts of diversification throughout their evolutionary history, events that can impact our ability to infer well-resolved phylogenies. However, groups that have undergone quick genetic change may have the highest adaptive potential. As the second largest genus of mushroom-forming fungi, Russula provides an excellent model for studying hyper-diversification and processes in evolution that drives it. This study focuses on the morphologically defined group - Russula subsection Roseinae. Species hypotheses based on morphological differentiation and multi-locus phylogenetic analyses are tested in the Roseinae using different applications of the multi-species coalescent model. Based on this combined approach, we recognize fourteen species in Roseinae including the Albida and wholly novel Magnarosea clades. Reconstruction of biogeographic and host association history suggest that parapatric speciation in refugia during glacial cycles of the Pleistocene drove diversification within the Roseinae, which is found to have a Laurasian distribution with an evolutionary origin in the Appalachian Mountains of eastern North America. Finally, we detect jump dispersal at a continental scale that has driven diversification since the most recent glacial cycles.

Population genetics revealMyotis keenii(Keen’s myotis) andMyotis evotis(long-eared myotis) to be a single species

Recognizing delineations of gene flow among groups of animals can be challenging but is necessary for conservation and management. Of particular importance is the identification of species boundaries. Several physical and genetic traits have been used with mixed success to distinguish Myotis keenii (Merriam, 1895) (Keen’s myotis) and Myotis evotis (H. Allen, 1864) (long-eared myotis), but it is unclear whether species distinction is biologically warranted. We generated 12–14 microsatellite locus genotypes for 275 long-eared Myotis representing four species — M. keenii, M. evotis, Myotis septentrionalis (Trouessart, 1897) (northern myotis), and Myotis thysanodes Miller, 1897 (fringed myotis) — from across northwestern North America and 23 Myotis lucifugus (Le Conte, 1831) (little brown myotis) as the outgroup. Population genetic analyses revealed four well-defined groups (species): M. septentrionalis, M. thysanodes, M. lucifugus, and a single group comprising M. keenii and M. evotis. We document high rates of gene flow within M. evotis/M. keenii. Cytochrome b gene (mtDNA) sequencing failed to resolve morphologically identifiable species. We highlight the importance of geographically thorough investigation of genetic connectivity (nuclear markers) when assessing taxonomic status of closely related groups. We document a morphometric cline within M. evotis/M. keenii that may in part explain earlier analyses that led to the description of the smaller bodied M. keenii (type locality: Haida Gwaii). We conclude that M. keenii does not qualify as a genetic or biological species.

Description of a new species of the Neotropical cichlid genus Gymnogeophagus Miranda Ribeiro, 1918 (Teleostei: Cichliformes) from the Middle Paraná basin, Misiones, Argentina

Gymnogeophagus jaryi, new species, is described from Southern tributaries of the Middle Paraná basin in Misiones. It can be distinguished from all other members of the genus, except from G. australis and G. caaguazuensis, by the presence of a hyaline to grey anterior portion of the dorsal fin. Gymnogeophagus jaryi differs from G. caaguazuensis by a longer caudal peduncle, caudal fin not lyrate, central portion of scales on dorsal portion of trunk light iridescent blue and by white spots in soft portion of dorsal fin in adult males, and from G. australis by the light iridescent blue coloration of central portion of scales on the dorsal portion of trunk and tail, and by the lack of scales on the soft portion of the dorsal fin. Additionally, it can be diagnosed by the following unique combination of characters: 10-11 dorsal-fin branched rays, 27-30 E1 scales, absence of lips thickening, and, in males, by the possession of a hump in adults, caudal fin not lyrate, presence of large white spots forming transversal stripes distally and in anterior area of the dorsal fin's soft portion, central area of scales on the dorsal portion of the trunk light iridescent blue, lack of scales on the base of the dorsal fin's soft portion, absence of a conspicuous and oblique dark band from the eye to the anterior border of the head, anterior portion of dorsal fin hyaline to grey, scales of the midlateral spot each bearing a semicircular light blue blotch, head hump starting at the horizontal through the eyes, concave anterior profile in lateral view, base of unpaired fins yellow, and whitish hyaline spots on caudal fin. The new species, based on mtDNA phylogeny, is the sister species of G. caaguazuensis from the Paraguay basin and is closely related to G. australis.

Modern Phylogenomics: Building Phylogenetic Trees Using the Multispecies Coalescent Model

The multispecies coalescent (MSC) model provides a compelling framework for building phylogenetic trees from multilocus DNA sequence data. The pure MSC is best thought of as a special case of so-called "multispecies network coalescent" models, in which gene flow is allowed among branches of the tree, whereas MSC methods assume there is no gene flow between diverging species. Early implementations of the MSC, such as "parsimony" or "democratic vote" approaches to combining information from multiple gene trees, as well as concatenation, in which DNA sequences from multiple gene trees are combined into a single "supergene," were quickly shown to be inconsistent in some regions of tree space, in so far as they converged on the incorrect species tree as more gene trees and sequence data were accumulated. The anomaly zone, a region of tree space in which the most frequent gene tree is different from the species tree, is one such region where many so-called "coalescent" methods are inconsistent. Second-generation implementations of the MSC employed Bayesian or likelihood models; these are consistent in all regions of gene tree space, but Bayesian methods in particular are incapable of handling the large phylogenomic data sets currently available. Two-step methods, such as MP-EST and ASTRAL, in which gene trees are first estimated and then combined to estimate an overarching species tree, are currently popular in part because they can handle large phylogenomic data sets. These methods are consistent in the anomaly zone but can sometimes provide inappropriate measures of tree support or apportion error and signal in the data inappropriately. MP-EST in particular employs a likelihood model which can be conveniently manipulated to perform statistical tests of competing species trees, incorporating the likelihood of the collected gene trees on each species tree in a likelihood ratio test. Such tests provide a useful alternative to the multilocus bootstrap, which only indirectly tests the appropriateness of competing species trees. We illustrate these tests and implementations of the MSC with examples and suggest that MSC methods are a useful class of models effectively using information from multiple loci to build phylogenetic trees.

Speciation in a biodiversity hotspot: Phylogenetic relationships, species delimitation, and divergence times of Patagonian ground frogs from the Eupsophus roseus group (Alsodidae)

The alsodid ground frogs of the Eupsophus genus are divided into two groups, the roseus (2n = 30) and vertebralis (2n = 28), which are distributed throughout the temperate Nothofagus forests of South America. Currently, the roseus group is composed by four species, while the vertebralis group consists of two. Phylogenetic relationships and species delimitation within each group are controversial. In fact, previous analyses considered that the roseus group was composed of between four to nine species. In this work, we evaluated phylogenetic relationships, diversification times, and species delimitation within the roseus group using a multi-locus dataset. For this purpose, mitochondrial (D-loop, Cyt b, and COI) and nuclear (POMC and CRYBA1) partial sequences from 164 individuals were amplified, representing all species. Maximum Likelihood (ML) and Bayesian approaches were used to reconstruct phylogenetic relationships. Species tree was estimated using BEAST and singular value decomposition scores for species quartets (SVDquartets). Species limits were evaluated with six coalescent approaches. Diversification times were estimated using mitochondrial and nuclear rates with LogNormal relaxed clock in BEAST. Nine well-supported monophyletic lineages were recovered in Bayesian, ML, and SVDquartets, including eight named species and a lineage composed by specimens from the Villarrica population (Bootstrap:>70, PP:> 0.99). Single-locus species delimitation analyses overestimated the species number in E. migueli, E. calcaratus, and E. roseus lineages, while multi-locus analyses recovered as species the nine lineages observed in phylogenetic analyses (Ctax = 0.69). It is hypothesized that Eupsophus diversification occurred during Mid-Pleistocene (0.42-0.14 Mya), with most species having originated after the Last Southern Patagonian Glaciation (0.18 Mya). Our results revitalize the hypothesis that the E. roseus group is composed of eight species and support the Villarrica lineage as a new putative species.

Integrating coalescent species delimitation with analysis of host specificity reveals extensive cryptic diversity despite minimal mitochondrial divergence in the malaria parasite genus Leucocytozoon

BACKGROUND

Coalescent methods that use multi-locus sequence data are powerful tools for identifying putatively reproductively isolated lineages, though this approach has rarely been used for the study of microbial groups that are likely to harbor many unrecognized species. Among microbial symbionts, integrating genetic species delimitation methods with trait data that could indicate reproductive isolation, such as host specificity data, has rarely been used despite its potential to inform species limits. Here we test the ability of an integrative approach combining genetic and host specificity data to delimit species within the avian malaria parasite genus Leucocytozoon in central Alaska.

RESULTS

We sequenced seven nuclear loci for 69 Leucocytozoon samples and used multiple species delimitation methods (GMYC and BPP models), tested for differences in host infection patterns among putative species based on 406 individual infections, and characterized parasite morphology. We found that cryptic morphology has masked a highly diverse Leucocytozoon assemblage, with most species delimitation methods recovering support for at least 21 separate species that occur sympatrically and have divergent host infection patterns. Reproductive isolation among putative species appears to have evolved despite low mtDNA divergence, and in one instance two Leucocytozoon cytb haplotypes that differed by a single base pair (~ 0.2% divergence) were supported as separate species. However, there was no consistent association between mtDNA divergence and species limits. Among cytb haplotypes that differed by one to three base pairs we observed idiosyncratic patterns of nuclear and ecological divergence, with cytb haplotype pairs found to be either conspecific, reproductively isolated with no divergence in host specificity, or reproductively isolated with divergent patterns of host specialization.

CONCLUSION

Integrating multi-locus genetic species delimitation methods and non-traditional ecological data types such as host specificity provide a novel view of the diversity of avian malaria parasites that has been missed previously using morphology and mtDNA barcodes. Species delimitation methods show that Leucocytozoon is highly species-rich in Alaska, and the genus is likely to harbor extraordinary species-level diversity worldwide. Integrating genetic and ecological data will be an important approach for understanding the diversity and evolutionary history of microbial symbionts moving forward.

A hyper-diverse genus of acanthocephalans revealed by tree-based and non-tree-based species delimitation methods: Ten cryptic species of Neoechinorhynchus in Middle American freshwater fishes

The genus Neoechinorhynchus represents a hyper-diverse group of acanthocephalans, parasites of fresh and brackish water fish and freshwater turtles, with approximately 116 species described worldwide. Forty-nine species have been recorded in the Americas, nine of them in Middle America. Even though species delimitation methods using DNA sequences have been rarely used for parasitic helminths, the genetic library for species of Neoechinorhynchus has grown in the past few years, enhancing the possibility of using these methods for inferring evolutionary relationships and for establishing more robust species boundaries. In this study, we used non-tree-based and tree-based methods through a coalescent approach to explore the species limits of specimens of Neoechinorhynchus collected in 57 localities across Middle America. We sequenced a large number of individuals to build a comprehensive dataset for three genes: the mitochondrial cytochrome c oxidase subunit I (352 individuals), the internal transcribed spacers (330 individuals), and the D2 + D3 domains of the large subunit (278 individuals). Several species delimitation methods were implemented, i.e., Automatic Barcode Gap Discovery (ABGD), General Mixed Yule-Coalescent Model (GMYC), Bayesian species delimitation (BPP) and species tree (∗BEAST). Additionally, we conducted a detailed morphological study of the diagnostic traits associated with the proboscis of 184 males and 169 females. Overall, our analyses allowed us to validate nine nominal species of Neoechinorhynchus and to identify 10 additional genetic lineages herein regarded as candidate species. This unexpected genetic diversity and the lack of reliable morphological traits show that the genus Neoechinorhynchus includes a group of cryptic species, at least in Middle America.

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.].

Conflicting Evolutionary Histories of the Mitochondrial and Nuclear Genomes in New World Myotis Bats

The rapid diversification of Myotis bats into more than 100 species is one of the most extensive mammalian radiations available for study. Efforts to understand relationships within Myotis have primarily utilized mitochondrial markers and trees inferred from nuclear markers lacked resolution. Our current understanding of relationships within Myotis is therefore biased towards a set of phylogenetic markers that may not reflect the history of the nuclear genome. To resolve this, we sequenced the full mitochondrial genomes of 37 representative Myotis, primarily from the New World, in conjunction with targeted sequencing of 3648 ultraconserved elements (UCEs). We inferred the phylogeny and explored the effects of concatenation and summary phylogenetic methods, as well as combinations of markers based on informativeness or levels of missing data, on our results. Of the 294 phylogenies generated from the nuclear UCE data, all are significantly different from phylogenies inferred using mitochondrial genomes. Even within the nuclear data, quartet frequencies indicate that around half of all UCE loci conflict with the estimated species tree. Several factors can drive such conflict, including incomplete lineage sorting, introgressive hybridization, or even phylogenetic error. Despite the degree of discordance between nuclear UCE loci and the mitochondrial genome and among UCE loci themselves, the most common nuclear topology is recovered in one quarter of all analyses with strong nodal support. Based on these results, we re-examine the evolutionary history of Myotis to better understand the phenomena driving their unique nuclear, mitochondrial, and biogeographic histories.

Speciation with Gene Flow in North American Myotis Bats

Growing evidence supports the idea that species can diverge in the presence of gene flow. However, most methods of phylogeny estimation do not consider this process, despite the fact that ignoring gene flow is known to bias phylogenetic inference. Furthermore, studies that do consider divergence-with-gene-flow typically do so by estimating rates of gene flow using a isolation-with-migration model (IM), rather than evaluating scenarios of gene flow (such as divergence-with-gene flow or secondary contact) that represent very different types of diversification. In this investigation, we aim to infer the recent phylogenetic history of a clade of western long-eared bats while evaluating a number of different models that parameterize gene flow in a variety of ways. We utilize PHRAPL, a new tool for phylogeographic model selection, to compare the fit of a broad set of demographic models that include divergence, migration, or both among Myotis evotis, $M$. thysanodes and M. keenii. A genomic data set consisting of 808 loci of ultraconserved elements was used to explore such models in three steps using an incremental design where each successive set was informed by, and thus more focused than, the previous set of models. Specifically, the three steps were to (i) assess whether gene flow should be modeled and identify the best topologies, (ii) infer directionality of migration using the best topologies, and (iii) estimate the timing of gene flow. The best model (AIC model weight ${\sim}0.98$) included two divergence events (($M$. evotis, $M$. thysanodes), M. keenii) accompanied by gene flow at the initial stages of divergence. These results provide a striking example of speciation-with-gene-flow in an evolutionary lineage. [Myotis bats; PHRAPL; P2C2M; phylogeographic model selection; speciation with gene flow.].

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.

The phylogeography of Myotis bat-associated rabies viruses across Canada

As rabies in carnivores is increasingly controlled throughout much of the Americas, bats are emerging as a significant source of rabies virus infection of humans and domestic animals. Knowledge of the bat species that maintain rabies is a crucial first step in reducing this public health problem. In North America, several bat species are known to be rabies virus reservoirs but the role of bats of the Myotis genus has been unclear due to the scarcity of laboratory confirmed cases and the challenges encountered in species identification of poorly preserved diagnostic submissions by morphological traits alone. This study has employed a collection of rabid bat specimens collected across Canada over a 25 year period to clearly define the role of particular Myotis species as rabies virus reservoirs. The virus was characterised by partial genome sequencing and host genetic barcoding, used to confirm species assignment of specimens, proved crucial to the identification of certain bat species as disease reservoirs. Several variants were associated with Myotis species limited in their Canadian range to the westernmost province of British Columbia while others were harboured by Myotis species that circulate across much of eastern and central Canada. All of these Myotis-associated viral variants, except for one, clustered as a monophyletic MYCAN clade, which has emerged from a lineage more broadly distributed across North America; in contrast one distinct variant, associated with the long-legged bat in Canada, represents a relatively recent host jump from a big brown bat reservoir. Together with evidence from South America, these findings demonstrate that rabies virus has emerged in the Myotis genus independently on multiple occasions and highlights the potential for emergence of new viral-host associations within this genus.

Reducing the public health burden of rabies is most effectively achieved by elimination of the disease from its hosts. While the role of dogs and many wild carnivore species in maintaining and transmitting rabies virus is well established, our understanding of the role of many bat species in this regard is still incomplete. Several North American bat species are known to be rabies virus reservoirs, but the role of bats of the Myotis genus has been unclear due to the very limited number of laboratory confirmed rabies cases detected in these bats and the challenge in assigning species to poorly preserved bat carcasses. Our study utilised a collection of rabid bats collected across Canada over 25 years to address this issue. Genetic barcoding was used to identify the specimens to species and the virus was characterised by partial genome sequencing. Host barcoding proved to be crucial for correct species assignment to many specimens and allowed the identification of certain Myotis species as reservoirs for several genetically distinct variants of the rabies virus. The geographical distribution of these viral variants, and its correlation with the ecological and taxonomic properties of the hosts, are described. While most of these variants belong to a single rabies virus lineage, one variant associated with the long-legged bat appears to have emerged relatively recently following transmission from a big brown bat reservoir. Consistent with this observation, comparison of the rabies viruses associated with Myotis bats in North and South America clearly shows that rabies has emerged in members of this genus on several separate occasions.

Multilocus approaches reveal underestimated species diversity and inter-specific gene flow in pikas (Ochotona) from southwestern China

The phylogeny of living pikas (Ochotonidae, Ochotona) remains obscure, and pika species diversity in southwestern China has never been well explored. In this study, 96 tissue samples from 11 valid species in three classified subgenera (Pika, Ochotona and Conothoa) from 23 locations were characterized using multilocus sequences of 7031bp. Two mitochondrial (CYT B and COI) and five nuclear gene segments (RAG1, RAG2, TTN, OXAIL and IL1RAPL1) were sequenced. We analysed evolutionary histories using maximum likelihood (RAxML) and Bayesian analyses (BEAST), and we also used molecular species delimitation analyses (BPP) to explore species diversity. Our study supported O. syrinx (O. huangensis) as a distinct clade from all named subgenera. Relationships among subgenera were not fully resolved, which may be due to a rapid diversification in the middle Miocene (∼13.90Ma). Conflicting gene trees implied mitochondrial introgression from O. cansus to O. curzoniae. We uncovered three cryptic species from Shaanxi, Sichuan and Yunnan with strong support, suggesting an underestimation of species diversity in the "sky-island" mountains of southwest China.

A Rapid and Scalable Method for Multilocus Species Delimitation Using Bayesian Model Comparison and Rooted Triplets

Multilocus sequence data provide far greater power to resolve species limits than the single locus data typically used for broad surveys of clades. However, current statistical methods based on a multispecies coalescent framework are computationally demanding, because of the number of possible delimitations that must be compared and time-consuming likelihood calculations. New methods are therefore needed to open up the power of multilocus approaches to larger systematic surveys. Here, we present a rapid and scalable method that introduces 2 new innovations. First, the method reduces the complexity of likelihood calculations by decomposing the tree into rooted triplets. The distribution of topologies for a triplet across multiple loci has a uniform trinomial distribution when the 3 individuals belong to the same species, but a skewed distribution if they belong to separate species with a form that is specified by the multispecies coalescent. A Bayesian model comparison framework was developed and the best delimitation found by comparing the product of posterior probabilities of all triplets. The second innovation is a new dynamic programming algorithm for finding the optimum delimitation from all those compatible with a guide tree by successively analyzing subtrees defined by each node. This algorithm removes the need for heuristic searches used by current methods, and guarantees that the best solution is found and potentially could be used in other systematic applications. We assessed the performance of the method with simulated, published, and newly generated data. Analyses of simulated data demonstrate that the combined method has favorable statistical properties and scalability with increasing sample sizes. Analyses of empirical data from both eukaryotes and prokaryotes demonstrate its potential for delimiting species in real cases.

Inferring phylogeny and speciation of Gymnosporangium species, and their coevolution with host plants

Gymnosporangium species (Pucciniaceae, Pucciniales) cause serious diseases and significant economic losses to apple cultivars. Most of the reported species are heteroecious and complete their life cycles on two different plant hosts belonging to two unrelated genera, i.e. Juniperus and Malus. However, the phylogenetic relationships among Gymnosporangium species and the evolutionary history of Gymnosporangium on its aecial and telial hosts were still undetermined. In this study, we recognized species based on rDNA sequence data by using coalescent method of generalized mixed Yule-coalescent (GMYC) and Poisson Tree Processes (PTP) models. The evolutionary relationships of Gymnosporangium species and their hosts were investigated by comparing the cophylogenetic analyses of Gymnosporangium species with Malus species and Juniperus species, respectively. The concordant results of GMYC and PTP analyses recognized 14 species including 12 known species and two undescribed species. In addition, host alternations of 10 Gymnosporangium species were uncovered by linking the derived sequences between their aecial and telial stages. This study revealed the evolutionary process of Gymnosporangium species, and clarified that the aecial hosts played more important roles than telial hosts in the speciation of Gymnosporangium species. Host switch, losses, duplication and failure to divergence all contributed to the speciation of Gymnosporangium 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.

An Integrative Approach for Understanding Diversity in the Punctelia rudecta Species Complex (Parmeliaceae, Ascomycota)

High levels of cryptic diversity have been documented in lichenized fungi, especially in Parmeliaceae, and integrating various lines of evidence, including coalescent-based species delimitation approaches, help establish more robust species circumscriptions. In this study, we used an integrative taxonomic approach to delimit species in the lichen-forming fungal genus Punctelia (Parmeliaceae), with a particular focus on the cosmopolitan species P. rudecta. Nuclear, mitochondrial ribosomal DNA and protein-coding DNA sequences were analyzed in phylogenetic and coalescence-based frameworks. Additionally, morphological, ecological and geographical features of the sampled specimens were evaluated. Five major strongly supported monophyletic clades were recognized in the genus Punctelia, and each clade could be characterized by distinct patterns in medullary chemistry. Punctelia rudecta as currently circumscribed was shown to be polyphyletic. A variety of empirical species delimitation methods provide evidence for a minimum of four geographically isolated species within the nominal taxon Punctelia rudecta, including a newly described saxicolous species, P. guanchica, and three corticolous species. In order to facilitate reliable sample identification for biodiversity, conservation, and air quality bio-monitoring research, these three species have been epitypified, in addition to the description of a new species.

Multilocus coalescent species delimitation reveals widespread cryptic differentiation among Drakensberg mountain-living freshwater crabs (Decapoda : Potamonautes)

Cryptic lineages present major challenges for evolutionary and conservation studies, particularly where these lineages remain undiscovered. Freshwater crabs are known to harbour cryptic diversity, in most cases with limited morphological differences. During the present study, we used a multilocus (12S rRNA, 16S rRNA, COI, 28S rRNA, DecapANT and PEPCK) Bayesian species delimitation to examine cryptic diversity within a freshwater crab species complex (Potamonautes clarus/P. depressus). We sampled 25 highland rivers in the Tugela and uMkomazi River drainage systems of the Drakensberg Mountain range, in the KwaZulu–Natal province of South Africa. Our results showed there to be at least eight lineages: six novel potamonautid freshwater crabs, and two described taxa P. clarus and P. depressus. Divergence from the most recent common ancestor occurred between the mid- and late Miocene (12.1 Mya), while divergence within the species complex occurred ~10.3 Mya up until the Holocene (0.11 Mya). The discovery of six novel lineages of freshwater crabs from a seemingly restricted distribution range has conservation implications, but to date most conservation planning strategies have focussed on freshwater vertebrates. By conducting a fine-scale phylogenetic survey using invertebrates, this study provides a platform for the inclusion of freshwater invertebrates in future conservation assessments.

Integrating a Numerical Taxonomic Method and Molecular Phylogeny for Species Delimitation of Melampsora Species (Melampsoraceae, Pucciniales) on Willows in China

The species in genus Melampsora are the causal agents of leaf rust diseases on willows in natural habitats and plantations. However, the classification and recognition of species diversity are challenging because morphological characteristics are scant and morphological variation in Melampsora on willows has not been thoroughly evaluated. Thus, the taxonomy of Melampsora species on willows remains confused, especially in China where 31 species were reported based on either European or Japanese taxonomic systems. To clarify the species boundaries of Melampsora species on willows in China, we tested two approaches for species delimitation inferred from morphological and molecular variations. Morphological species boundaries were determined based on numerical taxonomic analyses of morphological characteristics in the uredinial and telial stages by cluster analysis and one-way analysis of variance. Phylogenetic species boundaries were delineated based on the generalized mixed Yule-coalescent (GMYC) model analysis of the sequences of the internal transcribed spacer (ITS1 and ITS2) regions including the 5.8S and D1/D2 regions of the large nuclear subunit of the ribosomal RNA gene. Numerical taxonomic analyses of 14 morphological characteristics recognized in the uredinial-telial stages revealed 22 morphological species, whereas the GMYC results recovered 29 phylogenetic species. In total, 17 morphological species were in concordance with the phylogenetic species and 5 morphological species were in concordance with 12 phylogenetic species. Both the morphological and molecular data supported 14 morphological characteristics, including 5 newly recognized characteristics and 9 traditionally emphasized characteristics, as effective for the differentiation of Melampsora species on willows in China. Based on the concordance and discordance of the two species delimitation approaches, we concluded that integrative taxonomy by using both morphological and molecular variations was an effective approach for delimitating Melampsora species on willows in China.

Speciation with gene flow in whiptail lizards from a Neotropical xeric biome

Two main hypotheses have been proposed to explain the diversification of the Caatinga biota. The riverine barrier hypothesis (RBH) claims that the São Francisco River (SFR) is a major biogeographic barrier to gene flow. The Pleistocene climatic fluctuation hypothesis (PCH) states that gene flow, geographic genetic structure and demographic signatures on endemic Caatinga taxa were influenced by Quaternary climate fluctuation cycles. Herein, we analyse genetic diversity and structure, phylogeographic history, and diversification of a widespread Caatinga lizard (Cnemidophorus ocellifer) based on large geographical sampling for multiple loci to test the predictions derived from the RBH and PCH. We inferred two well-delimited lineages (Northeast and Southwest) that have diverged along the Cerrado-Caatinga border during the Mid-Late Miocene (6-14 Ma) despite the presence of gene flow. We reject both major hypotheses proposed to explain diversification in the Caatinga. Surprisingly, our results revealed a striking complex diversification pattern where the Northeast lineage originated as a founder effect from a few individuals located along the edge of the Southwest lineage that eventually expanded throughout the Caatinga. The Southwest lineage is more diverse, older and associated with the Cerrado-Caatinga boundaries. Finally, we suggest that C. ocellifer from the Caatinga is composed of two distinct species. Our data support speciation in the presence of gene flow and highlight the role of environmental gradients in the diversification process.

Population genetic structure of a common host predicts the spread of white-nose syndrome, an emerging infectious disease in bats

Landscape complexity influences patterns of animal dispersal, which in turn may affect both gene flow and the spread of pathogens. White-nose syndrome (WNS) is an introduced fungal disease that has spread rapidly throughout eastern North America, causing massive mortality in bat populations. We tested for a relationship between the population genetic structure of the most common host, the little brown myotis (Myotis lucifugus), and the geographic spread of WNS to date by evaluating logistic regression models of WNS risk among hibernating colonies in eastern North America. We hypothesized that risk of WNS to susceptible host colonies should increase with both geographic proximity and genetic similarity, reflecting historical connectivity, to infected colonies. Consistent with this hypothesis, inclusion of genetic distance between infected and susceptible colonies significantly improved models of disease spread, capturing heterogeneity in the spatial expansion of WNS despite low levels of genetic differentiation among eastern populations. Expanding our genetic analysis to the continental range of little brown myotis reveals strongly contrasting patterns of population structure between eastern and western North America. Genetic structure increases markedly moving westward into the northern Great Plains, beyond the current distribution of WNS. In western North America, genetic differentiation of geographically proximate populations often exceeds levels observed across the entire eastern region, suggesting infrequent and/or locally restricted dispersal, and thus relatively limited opportunities for pathogen introduction in western North America. Taken together, our analyses suggest a possibly slower future rate of spread of the WNS pathogen, at least as mediated by little brown myotis.