Perfectly misleading distances from ternary characters

D. Huson and M. Steel showed that for any two binary phylogenetic trees on the same set of n taxa, there exists a sequence of multistate characters that is homoplasy-free only on the first tree but perfectly additive only on the second one. The original construction of such a sequence required n - 1 character states and it remained an open question whether a sequence using fewer character states can always be found. In the present note we will answer this question by showing that three character states suffice to construct such misleading sequences--even if we insist that they conform to an ultrametric (i.e., fit a molecular clock).     

Reconstructing an ultrametric galled phylogenetic network from a distance matrix

Given a distance matrix M that specifies the pairwise evolutionary distances between n species, the phylogenetic tree reconstruction problem asks for an edge-weighted phylogenetic tree that satisfies M, if one exists. We study some extensions of this problem to rooted phylogenetic networks. Our main result is an O(n(2) log n)-time algorithm for determining whether there is an ultrametric galled network that satisfies M, and if so, constructing one. In fact, if such an ultrametric galled network exists, our algorithm is guaranteed to construct one containing the minimum possible number of nodes with more than one parent (hybrid nodes). We also prove that finding a largest possible submatrix M' of M such that there exists an ultrametric galled network that satisfies M' is NP-hard. Furthermore, we show that given an incomplete distance matrix (i.e. where some matrix entries are missing), it is also NP-hard to determine whether there exists an ultrametric galled network which satisfies it.     

Is cladogenesis heritable?

The heritability of speciation rates and extinction risks is a crucial parameter in models of macroevolution, but little direct evidence has been available to assess the occurrence, strength, or generality of this heritability. We tested for heritability using correlations between ancestral and descendent branch lengths in phylogenetic trees, an approach first applied to a bird phylogeny by Harvey et al. (1991, pages 123-137 in Genes in ecology [R. J. Berry et al., eds.], Blackwell Scientific, Oxford). We applied Harvey et al.'s test to some of the largest DNA sequence-based phylogenetic analyses published to date for plants, insects, fungi, and bacteria. If one of two parent lineages splits first and if this is the case for any heritable reason, then on average we expect its daughter lineages to also split first. We also used a randomization procedure to assess significance of branch length heritability. Using maximum parsimony and maximum likelihood branch lengths and trees made ultrametric after nonparametric rate smoothing or by enforcing a molecular clock, we found a pattern for most clades consistent with heritable net cladogenesis. Heritability of cladogenesis may be a general phenomenon, detectable across a large number of lineages and a broad range of taxa.     

Molecular clock fork phylogenies: closed form analytic maximum likelihood solutions

Maximum likelihood (ML) is increasingly used as an optimality criterion for selecting evolutionary trees, but finding the global optimum is a hard computational task. Because no general analytic solution is known, numeric techniques such as hill climbing or expectation maximization (EM) are used in order to find optimal parameters for a given tree. So far, analytic solutions were derived only for the simplest model-three-taxa, two-state characters, under a molecular clock. Quoting Ziheng Yang, who initiated the analytic approach,"this seems to be the simplest case, but has many of the conceptual and statistical complexities involved in phylogenetic estimation."In this work, we give general analytic solutions for a family of trees with four-taxa, two-state characters, under a molecular clock. The change from three to four taxa incurs a major increase in the complexity of the underlying algebraic system, and requires novel techniques and approaches. We start by presenting the general maximum likelihood problem on phylogenetic trees as a constrained optimization problem, and the resulting system of polynomial equations. In full generality, it is infeasible to solve this system, therefore specialized tools for the molecular clock case are developed. Four-taxa rooted trees have two topologies-the fork (two subtrees with two leaves each) and the comb (one subtree with three leaves, the other with a single leaf). We combine the ultrametric properties of molecular clock fork trees with the Hadamard conjugation to derive a number of topology dependent identities. Employing these identities, we substantially simplify the system of polynomial equations for the fork. We finally employ symbolic algebra software to obtain closed formanalytic solutions (expressed parametrically in the input data). In general, four-taxa trees can have multiple ML points. In contrast, we can now prove that each fork topology has a unique(local and global) ML point.     

Testing a molecular clock without an outgroup: derivations of induced priors on branch-length restrictions in a Bayesian framework

We propose a Bayesian method for testing molecular clock hypotheses for use with aligned sequence data from multiple taxa. Our method utilizes a nonreversible nucleotide substitution model to avoid the necessity of specifying either a known tree relating the taxa or an outgroup for rooting the tree. We employ reversible jump Markov chain Monte Carlo to sample from the posterior distribution of the phylogenetic model parameters and conduct hypothesis testing using Bayes factors, the ratio of the posterior to prior odds of competing models. Here, the Bayes factors reflect the relative support of the sequence data for equal rates of evolutionary change between taxa versus unequal rates, averaged over all possible phylogenetic parameters, including the tree and root position. As the molecular clock model is a restriction of the more general unequal rates model, we use the Savage-Dickey ratio to estimate the Bayes factors. The Savage-Dickey ratio provides a convenient approach to calculating Bayes factors in favor of sharp hypotheses. Critical to calculating the Savage-Dickey ratio is a determination of the prior induced on the modeling restrictions. We demonstrate our method on a well-studied mtDNA sequence data set consisting of nine primates. We find strong support against a global molecular clock, but do find support for a local clock among the anthropoids. We provide mathematical derivations of the induced priors on branch length restrictions assuming equally likely trees. These derivations also have more general applicability to the examination of prior assumptions in Bayesian phylogenetics.     

Maximum likelihood molecular clock comb: analytic solutions.

Maximum likelihood (ML) is increasingly used as an optimality criterion for selecting evolutionary trees, but finding the global optimum is a hard computational task. Because no general analytic solution is known, numeric techniques such as hill climbing or expectation maximization (EM), are used in order to find optimal parameters for a given tree. So far, analytic solutions were derived only for the simplest model--three taxa, two state characters, under a molecular clock. Four taxa rooted trees have two topologies--the fork (two subtrees with two leaves each) and the comb (one subtree with three leaves, the other with a single leaf). In a previous work, we devised a closed form analytic solution for the ML molecular clock fork. In this work, we extend the state of the art in the area of analytic solutions ML trees to the family of all four taxa trees under the molecular clock assumption. The change from the fork topology to the comb incurs a major increase in the complexity of the underlying algebraic system and requires novel techniques and approaches. We combine the ultrametric properties of molecular clock trees with the Hadamard conjugation to derive a number of topology dependent identities. Employing these identities, we substantially simplify the system of polynomial equations. We finally use tools from algebraic geometry (e.g., Gr?bner bases, ideal saturation, resultants) and employ symbolic algebra software to obtain analytic solutions for the comb. We show that in contrast to the fork, the comb has no closed form solutions (expressed by radicals in the input data). In general, four taxa trees can have multiple ML points. In contrast, we can now prove that under the molecular clock assumption, the comb has a unique (local and global) ML point. (Such uniqueness was previously shown for the fork.).     

A mitogenomic timescale for birds detects variable phylogenetic rates of molecular evolution and refutes the standard molecular clock

Current understanding of the diversification of birds is hindered by their incomplete fossil record and uncertainty in phylogenetic relationships and phylogenetic rates of molecular evolution. Here we performed the first comprehensive analysis of mitogenomic data of 48 vertebrates, including 35 birds, to derive a Bayesian timescale for avian evolution and to estimate rates of DNA evolution. Our approach used multiple fossil time constraints scattered throughout the phylogenetic tree and accounts for uncertainties in time constraints, branch lengths, and heterogeneity of rates of DNA evolution. We estimated that the major vertebrate lineages originated in the Permian; the 95% credible intervals of our estimated ages of the origin of archosaurs (258 MYA), the amniote-amphibian split (356 MYA), and the archosaur-lizard divergence (278 MYA) bracket estimates from the fossil record. The origin of modern orders of birds was estimated to have occurred throughout the Cretaceous beginning about 139 MYA, arguing against a cataclysmic extinction of lineages at the Cretaceous/Tertiary boundary. We identified fossils that are useful as time constraints within vertebrates. Our timescale reveals that rates of molecular evolution vary across genes and among taxa through time, thereby refuting the widely used mitogenomic or cytochrome b molecular clock in birds. Moreover, the 5-Myr divergence time assumed between 2 genera of geese (Branta and Anser) to originally calibrate the standard mitochondrial clock rate of 0.01 substitutions per site per lineage per Myr (s/s/l/Myr) in birds was shown to be underestimated by about 9.5 Myr. Phylogenetic rates in birds vary between 0.0009 and 0.012 s/s/l/Myr, indicating that many phylogenetic splits among avian taxa also have been underestimated and need to be revised. We found no support for the hypothesis that the molecular clock in birds "ticks" according to a constant rate of substitution per unit of mass-specific metabolic energy rather than per unit of time, as recently suggested. Our analysis advances knowledge of rates of DNA evolution across birds and other vertebrates and will, therefore, aid comparative biology studies that seek to infer the origin and timing of major adaptive shifts in vertebrates.     

Bias and conflict in phylogenetic inference of myco-heterotrophic plants: a case study in Thismiaceae

Due to morphological reduction and absence of amplifiable plastid genes, the identification of photosynthetic relatives of heterotrophic plants is problematic. Although nuclear and mitochondrial gene sequences may offer a welcome alternative source of phylogenetic markers, the presence of rate heterogeneity in these genes may introduce bias/systematic error in phylogenetic analyses. We examine the phylogenetic position of Thismiaceae based on nuclear 18S rDNA and mitochondrial atpA DNA sequence data, as well as using parsimony, likelihood and Bayesian inference methods. Significant differences in evolutionary rates of these genes between closely related taxa lead to conflicting results: while parsimony analyses of 18S rDNA and combined data strongly support the monophyly of Thismiaceae, Bayesian inference, with and without a relaxed molecular clock, as well as the Swofford–Olsen–Waddell–Hillis (SOWH) test confidently reject this hypothesis. We show that rate heterogeneity in our data leads to long-branch attraction artifacts in parsimony analysis. However, using model-based inference methods the question of whether Thismiaceae are monophyletic remains elusive. On the one hand maximum likelihood nonparametric bootstrapping and parametric hypothesis tests fail to support a paraphyletic Thismiaceae, on the other hand Bayesian inference methods (both without and with a relaxed clock) significantly reject a monophyletic Thismiaceae. These results show that an adequate sampling, the use of rate homogeneous data, and the application of different inference methods are important factors for developing phylogenetic hypotheses of myco-heterotrophic plants. © The Willi Hennig Society 2009.     

Phylogeny and evolution of Cervidae based on complete mitochondrial genomes

Mitochondrial DNA sequences can be used to estimate phylogenetic relationships among animal taxa and for molecular phylogenetic evolution analysis. With the development of sequencing technology, more and more mitochondrial sequences have been made available in public databases, including whole mitochondrial DNA sequences. These data have been used for phylogenetic analysis of animal species, and for studies of evolutionary processes. We made phylogenetic analyses of 19 species of Cervidae, with Bos taurus as the outgroup. We used neighbor joining, maximum likelihood, maximum parsimony, and Bayesian inference methods on whole mitochondrial genome sequences. The consensus phylogenetic trees supported monophyly of the family Cervidae; it was divided into two subfamilies, Plesiometacarpalia and Telemetacarpalia, and four tribes, Cervinae, Muntiacinae, Hydropotinae, and Odocoileinae. The divergence times in these families were estimated by phylogenetic analysis using the Bayesian method with a relaxed molecular clock method; the results were consistent with those of previous studies. We concluded that the evolutionary structure of the family Cervidae can be reconstructed by phylogenetic analysis based on whole mitochondrial genomes; this method could be used broadly in phylogenetic evolutionary analysis of animal taxa.     

Disentangling causes of disjunction on the South Island of New Zealand: the Alpine fault hypothesis of vicariance revisited

Many elements of the flora and fauna of New Zealand's South Island show disjunct distributions with conspecific populations or closely-related species that occur in the north-west and south separated by a central gap. Three events have been implicated to account for this pattern: Pleistocene glaciations, Pliocene mountain building, or displacement along the Alpine fault, the border of the Pacific and Australian plates stretching diagonally across the South Island from south-west to north-east that formed during the Miocene. Disjunct distributions of species level taxa are probably too young to be due to Alpine fault vicariance. It has therefore been suggested that the biogeographical impact of the Alpine fault, if any, should be apparent on deeper phylogenetic levels. We tested this hypothesis by reconstructing the phylogenetic relationships of the hydrobiid gastropods of New Zealand based on mitochondrial DNA fragments of cytochrome oxidase subunit I (CO I ) and 16S rDNA. The creno- and stygobiont species of this family are typically poor dispersers. Therefore, ancient patterns of distribution may be conserved. The phylogenetic reconstructions were in accordance with the Alpine fault hypothesis uniting genera occurring on either side of the fault. Divergence estimates based on a molecular clock of CO I indicated splits predating the Pliocene uplift of the Alps.  © 2007 The Linnean Society of London, Biological Journal of the Linnean Society , 2007, 91 , 361–374.     

Information-theoretic indices and an approximate significance test for testing the molecular clock hypothesis with genetic distances

Distance-based phylogenetic methods are widely used in biomedical research. However, distance-based dating of speciation events and the test of the molecular clock hypothesis are relatively underdeveloped. Here I develop an approximate test of the molecular clock hypothesis for distance-based trees, as well as information-theoretic indices that have been used frequently in model selection, for use with distance matrices. The results are in good agreement with the conventional sequence-based likelihood ratio test. Among the information-theoretic indices, AICu is the most consistent with the sequence-based likelihood ratio test. The confidence in model selection by the indices can be evaluated by bootstrapping. I illustrate the usage of the indices and the approximate significance test with both empirical and simulated sequences. The tests show that distance matrices from protein gel electrophoresis and from genome rearrangement events do not violate the molecular clock hypothesis, and that the evolution of the third codon position conforms to the molecular clock hypothesis better than the second codon position in vertebrate mitochondrial genes. I outlined evolutionary distances that are appropriate for phylogenetic reconstruction and dating.     

Caveats on the use of fossil calibrations for molecular dating: a comment on Near et al

We reassess a study on a fossil-calibrated molecular clock that provides a new method for evaluating the accuracy of calibration points. We address several pitfalls that molecular systematists should be aware of when calculating rates of molecular evolution based on fossil calibrations. These caveats involve the substantiation and accurate use of geologic dates, the inappropriate use of fixed calibration points, and the explicit and objective phylogenetic placement of fossil taxa. Paleontological data, like molecular data, should be treated with the utmost rigor.     

Phylogenetic analysis and trait evolution in Australian lineages of drywood termites (Isoptera, Kalotermitidae)

A phylogenetic analysis of Australian drywood termites (Isoptera, Kalotermitidae) based on partial sequence from the cytochrome oxidase II (COII) and cytochrome b genes is presented. In addition to providing new information on the evolutionary relationships among 25 species from seven genera, we evaluate the relative likelihoods of alternative topological hypotheses, including those derived from morphology-based classifications. We also test the applicability of a molecular clock for estimating the age of the Kalotermitidae and infer the evolution of species-specific variation for habitat type and soldier caste phragmosis by mapping this information onto the independently derived phylogeny. Maximum-likelihood analysis of both nucleotide and protein sequences from a multigene data set jointly support a single topology, which is shown to be the best estimate of the true phylogeny among the alternatives tested. Our results support the monophyly of all genera but question the discrimination between Procryptotermes and Cryptotermes. A basal dichotomy among generic groups suggests two principle lines of divergence within the family. Intergeneric relationships show mixed congruence to previous proposals, resulting in one morphology-based classification being rejected. A molecular clock hypothesis is not supported due to significant among-lineage rate heterogeneity in the COII gene. Patterns revealed through trait mapping suggest that the most recently diverged taxa tend to occupy the driest habitats and that these same taxa reflect a defensive transition away from large mandibulate soldiers toward small phragmotic soldiers. The association between habitat and defensibility supports the hypothesis that these two characters have been tightly linked throughout the social diversification of termites.     

Molecular phylogenetics of porcini mushrooms (Boletus section Boletus)

Porcini (Boletus section Boletus: Boletaceae: Boletineae: Boletales) are a conspicuous group of wild, edible mushrooms characterized by fleshy fruiting bodies with a poroid hymenophore that is "stuffed" with white hyphae when young. Their reported distribution is with ectomycorrhizal plants throughout the Northern Hemisphere. Little progress has been made on the systematics of this group using modern molecular phylogenetic tools because sampling has been limited primarily to European species and the genes employed were insufficient to resolve the phylogeny. We examined the evolutionary history of porcini by using a global geographic sampling of most known species, new discoveries from little explored areas, and multiple genes. We used 78 sequences from the fast-evolving nuclear internal transcribed spacers and are able to recognize 18 reciprocally monophyletic species. To address whether or not porcini form a monophyletic group, we compiled a broadly sampled dataset of 41 taxa, including other members of the Boletineae, and used separate and combined phylogenetic analysis of sequences from the nuclear large subunit ribosomal DNA, the largest subunit of RNA polymerase II, and the mitochondrial ATPase subunit six gene. Contrary to previous studies, our separate and combined phylogenetic analyses support the monophyly of porcini. We also report the discovery of two taxa that expand the known distribution of porcini to Australia and Thailand and have ancient phylogenetic connections to the rest of the group. A relaxed molecular clock analysis with these new taxa dates the origin of porcini to between 42 and 54 million years ago, coinciding with the initial diversification of angiosperms, during the Eocene epoch when the climate was warm and humid. These results reveal an unexpected diversity, distribution, and ancient origin of a group of commercially valuable mushrooms that may provide an economic incentive for conservation and support the hypothesis of a tropical origin of the ectomycorrhizal symbiosis.     

Sociality and the rate of molecular evolution

The molecular clock does not tick at a uniform rate in all taxa but may be influenced by species characteristics. Eusocial species (those with reproductive division of labor) have been predicted to have faster rates of molecular evolution than their nonsocial relatives because of greatly reduced effective population size; if most individuals in a population are nonreproductive and only one or few queens produce all the offspring, then eusocial animals could have much lower effective population sizes than their solitary relatives, which should increase the rate of substitution of "nearly neutral" mutations. An earlier study reported faster rates in eusocial honeybees and vespid wasps but failed to correct for phylogenetic nonindependence or to distinguish between potential causes of rate variation. Because sociality has evolved independently in many different lineages, it is possible to conduct a more wide-ranging study to test the generality of the relationship. We have conducted a comparative analysis of 25 phylogenetically independent pairs of social lineages and their nonsocial relatives, including bees, wasps, ants, termites, shrimps, and mole rats, using a range of available DNA sequences (mitochondrial and nuclear DNA coding for proteins and RNAs, and nontranslated sequences). By including a wide range of social taxa, we were able to test whether there is a general influence of sociality on rates of molecular evolution and to test specific predictions of the hypothesis: (1) that social species have faster rates because they have reduced effective population sizes; (2) that mitochondrial genes would show a greater effect of sociality than nuclear genes; and (3) that rates of molecular evolution should be correlated with the degree of sociality. We find no consistent pattern in rates of molecular evolution between social and nonsocial lineages and no evidence that mitochondrial genes show faster rates in social taxa. However, we show that the most highly eusocial Hymenoptera do have faster rates than their nonsocial relatives. We also find that social parasites (that utilize the workers from related species to produce their own offspring) have faster rates than their social relatives, which is consistent with an effect of lower effective population size on rate of molecular evolution. Our results illustrate the importance of allowing for phylogenetic nonindependence when conducting investigations of determinants of variation in rate of molecular evolution.     

Testing biogeographical hypotheses in space and time: faunal relationships of the putative ancient Lake Egˇirdir in Asia Minor

Aim The aims of this study are to establish a multi‐locus phylogeny‐based hypothesis for the biogeographical relationship of gastropods from the putative ancient Lake Egˇirdir, to test the respective null hypothesis, to estimate the timing of biogeographical events based on independent molecular clock approaches, and to interpret the data with respect to the putative ancient character of Lake Egˇirdir. Location Lake Egˇirdir, western Taurus Lake District, Turkey. Methods DNA sequences from the putatively only extant endemic taxon of Lake Egˇirdir, Falsipyrgula pfeiferi, as well as representatives of other pyrgulinid genera from Europe and western Asia are used for phylogenetic analyses based on Bayesian inference. The respective null hypothesis is tested utilizing parametric bootstrapping. The timing of evolutionary events is estimated based on two independent molecular clock approaches, which involve the modelling of judicious errors associated with branch‐length calculations and calibration points. Results Bayesian inference indicates a very close relationship between the Lake Egˇirdir and Ponto‐Caspian taxa. Parametric bootstrapping rejects the null hypothesis that these taxa are not monophyletic (P ≤ 0.01). The alternative hypothesis, namely monophyly of the Ponto‐Caspian and Egˇirdir taxa, can therefore be accepted. The two independent molecular clock approaches show diversion times for the Ponto‐Caspian/Egˇirdir taxa of 0.42 ± 0.18 and 0.43 ± 0.63 Ma. Main conclusions The present study shows that there is no close biogeographical affiliation between the probably only remaining endemic taxon of Lake Egˇirdir and taxa from central Europe or the Balkan region. Instead, there is a very close and relatively young (i.e. late Pleistocene) biogeographical relationship with the Ponto‐Caspian pyrgulinids. However, fossil and comparative data from other invertebrates indicate that biogeographical connections between Lake Egˇirdir and the Ponto‐Caspian region existed during various time periods, i.e. the Miocene/Pliocene, early Pleistocene, and late Pleistocene. Acknowledging the still‐restricted knowledge of the evolutionary history of the lake, the data presented here do not reject the putative ancient status of Lake Egˇirdir. Future studies utilizing endemic taxa of other lakes in the region need to show whether the western Taurus Lake District represents a melting pot of Pleistocene refuge biodiversity from different regions, and whether the admixture of divergent lineages has created a genetically distinct set of taxa that would warrant the designation of the area as a unique biogeographical subregion.     

Flies on thistles: support for synchronous speciation?

Synchronous speciation of hosts and herbivorous insects predicts a congruent topology of host and insect phylogenies and similar evolutionary ages of host and insect taxa. To test these predictions for the specialized herbivorous fly genus Urophora (Diptera: Tephritidae), we used three different approaches. (i) We generated a phylogenetic tree of 11 European Urophora species from allozyme data and constructed a phylogeny of their hosts from published sources. Superimposing the Urophora tree on the host-plant tree we found no evidence for general congruence. (ii) We correlated genetic distances (Nei distances) of the host plants vs. the genetic distances of associated Urophora species. Overall, the relationship was not positive. Nevertheless, for some pairs of Urophora species and host plants genetic distances were in the same order of magnitude. (iii) We collected allozyme data for pairs of thistle taxa and pairs of herbivores on thistles together with independent time estimates. With these data we calibrated a molecular clock. There was a non-linear relationship between phylogenetic age and genetic distance, rendering the dating of deep events in thistle–insect evolution difficult. Nevertheless the derived molecular clock showed that the split of insect taxa lagged behind the split of hosts.  © 2005 The Linnean Society of London, Biological Journal of the Linnean Society , 2005, 84 , 775–783.     

Metrics for Phylogenetic Networks II: Nodal and Triplets Metrics

The assessment of phylogenetic network reconstruction methods requires the ability to compare phylogenetic networks. This is the second in a series of papers devoted to the analysis and comparison of metrics for tree-child time consistent phylogenetic networks on the same set of taxa. In this paper, we generalize to phylogenetic networks two metrics that have already been introduced in the literature for phylogenetic trees: the nodal distance and the triplets distance. We prove that they are metrics on any class of tree-child time consistent phylogenetic networks on the same set of taxa, as well as some basic properties for them. To prove these results, we introduce a reduction/expansion procedure that can be used not only to establish properties of tree-child time consistent phylogenetic networks by induction, but also to generate all tree-child time consistent phylogenetic networks with a given number of leaves.     

Species identification in the taxonomically neglected, highly diverse, neotropical parasitoid wasp genus Notiospathius (Braconidae: Doryctinae) based on an integrative molecular and morphological approach

Various DNA sequence-based methods for species delineation have recently been developed to assess the species-richness of highly diverse, neglected invertebrate taxa. These methods, however, need to be tested under a variety of conditions, including the use of different markers and parameters. Here, we explored the species diversity of a species-rich group of braconid parasitoid wasps, the Neotropical genus Notiospathius, including 233 specimens from 10 different countries. We examined sequences of two mitochondrial (mt) (COI, cyt b) and one nuclear (wg) gene fragments. We analysed them separately as well as concatenating the mt data with the general mixed Yule-coalescent (GMYC) model for species delineation using different tree-building methods and parameters for reconstructing ultrametric trees. We evaluated the performance of GMYC analyses by comparing their species delineations with our morphospecies identifications. Reconstructing ultrametric trees with a relaxed lognormal clock rate using the program BEAST gave the most congruent results with morphology for the two mt markers. A tree obtained with wg using the programs MrBayes+Pathd8 had the fewest cases of incongruence with morphology, though the performance of this nuclear marker was considerably lower than that of COI and cyt b. Species delimitation using the coalescent prior to obtain ultrametric trees was morphologically more congruent with COI, whereas the Yule prior was more congruent with cyt b. The analyses concatenating the mt datasets failed to recover some species supported both by morphology and the separate analyses of the mt markers. The highest morphological congruence was obtained with the GMYC analysis on an ultrametric tree reconstructed with cyt b using the relaxed lognormal clock rate and the Yule prior, thus supporting the importance of using alternative markers when the information of the barcoding locus (COI) is not concordant with morphological evidence. Seventy-one species were delimited based on the congruence found among COI, cyt b and morphology. Both mt markers also revealed the existence of seven potential cryptic species. This high species richness from a scattered geographical sampling indicates that there is a remarkable number of Notiospathius species that remains undiscovered.