Bootstrap Proportions Research Articles

Two New Edible Lyophyllum Species from Tibetan Areas, China

Two new species, Lyophyllum yiqunyang and L. heimogu, that belong to the section Difformia of the genus Lyophyllum, are described based on collections from Tibetan areas, China. The two species are delicious edible low-temperature mushrooms and are widely collected and eaten by local people. Lyophyllum yiqunyang sp. nov. is saprotrophic and has medium-sized basidiomata, olive-grey pileus, cheilocystidia, absent pleurocystidia, globose to subglobose basidiospores (6.12–6.31 × 6.02–6.23 μm) and clamp connections at the pileus context, hymenophoral trama and stipe. Lyophyllum heimogu sp. nov. is saprotrophic and has a dark grey to olive pileus, medium-sized basidiomata and globose to subglobose basidiospores (5.31–5.63 × 5.22–5.41 μm). In the phylogenetic analyses, our two new species formed distinct clades that are well supported by posterior probabilities and bootstrap proportions. Detailed descriptions, colour photos, illustrations and a phylogenetic tree to show the positions of the two new species are presented.

Robustness of Felsenstein's Versus Transfer Bootstrap Supports With Respect to Taxon Sampling.

The bootstrap method is based on resampling sequence alignments and re-estimating trees. Felsenstein's bootstrap proportions (FBP) are the most common approach to assess the reliability and robustness of sequence-based phylogenies. However, when increasing taxon sampling (i.e., the number of sequences) to hundreds or thousands of taxa, FBP tend to return low support for deep branches. The transfer bootstrap expectation (TBE) has been recently suggested as an alternative to FBP. TBE is measured using a continuous transfer index in [0,1] for each bootstrap tree, instead of the binary {0,1} index used in FBP to measure the presence/absence of the branch of interest. TBE has been shown to yield higher and more informative supports while inducing a very low number of falsely supported branches. Nonetheless, it has been argued that TBE must be used with care due to sampling issues, especially in datasets with a high number of closely related taxa. In this study, we conduct multiple experiments by varying taxon sampling and comparing FBP and TBE support values on different phylogenetic depths, using empirical datasets. Our results show that the main critique of TBE stands in extreme cases with shallow branches and highly unbalanced sampling among clades, but that TBE is still robust in most cases, while FBP is inescapably negatively impacted by high taxon sampling. We suggest guidelines and good practices in TBE (and FBP) computing and interpretation.

Introduction of two novel species of Hymenopellis (Agaricales, Physalacriaceae) from Thailand.

Hymenopellis is the most diverse genus in the group of oudemansielloid/xeruloid taxa (Physalacriaceae). This genus has a worldwide distribution with records mostly from Europe and America. Asian taxa are least represented. In this paper on Hymenopellis from Thailand, two novel species are introduced, and a Hymenopellis collection affine to H.orientalis is described. Macro and micromorphological characters are described. Maximum likelihood and Bayesian phylogenetic analyses were performed on combined ITS and nrLSU regions to confirm taxonomical placement and infer the phylogenetic affinities of the studied species. Hymenopellisstraminea sp. nov. is straw-yellow, with medium-sized basidiomata, abundant and diverse in form cheilocystidia, few, narrowly lageniform to fusiform pleurocystidia, and clamp connections at the lower part of the stipe. Hymenopellisutriformis sp. nov. has mostly utriform pleurocystidia and 2-spored basidia. In the inferred phylogenies, the new species from this study formed distinct clades well supported by bootstrap proportions and posterior probabilities. The studied specimen affine to H.orientalis produced 2-spored basidia whereas published descriptions of other specimens mention 4-spored basidia. Moreover, the genetic distance between ITS sequences of this specimen and that of a Hymenopellisorientalis specimen from GenBank was 1.30-2.57%. Therefore, the conspecificity of our specimen with H.orientalis is uncertain, and additional specimens are needed to fully confirm its identity.

Loop-mediated isothermal amplification assay for the rapid detection of Anaplasma marginale in cattle based on major surface protein 5 gene.

Timely intervention is needed to minimize the economic losses of vector-borne bovine anaplasmosis which can be possible by the isothermal amplification assay. Anaplasma marginale in the cattle of south Gujarat, India was detected in the PCR and LAMP by amplifying the fragment of msp5 gene. The PCR product was digested with EcoRI, and sequenced to confirm its pathogen specific detection. Species specific PCR observed a band of 457 bp of msp5 DNA following 1% agarose gel electrophoresis. Positive LAMP reaction turned into yellow colour while negative sample depicted original pink colour. A detection limit of PCR and LAMP was up to 10-6 and 10-8 of the original genomic DNA of A. marginale, respectively. A single cut site of EcoRI was observed in the PCR product. Current msp5 DNA sequences of A. marginale (MW538962 and MW538961) showed 100% homology with the published sequences. Monophyletic lineage type relationship was observed with high bootstrap proportion among the msp5 DNA sequences of A. marginale in the phylogram. Prevalence rate of A. marginale was significantly higher (p<0.05) in the PCR [43/280 (15.36%)] and LAMP [62/280 (22.14%)] than the microscopic technique [17/280 (6.07%)]. Diagnostic sensitivity, specificity, positive and negative predictive values at 95% CI for LAMP assay with respect to PCR were 93.02%, 90.72%, 64.52% and 98.62%, respectively. Thus LAMP can be a practical alternative to the PCR for the diagnosis of A. marginale infection in the cattle even in field condition.

Quantifying the Effectiveness of Advertising: A Bootstrap Proportion Test for Brand Lift Testing

Quantifying the Effectiveness of Advertising: A Bootstrap Proportion Test for Brand Lift Testing

On the Need for New Measures of Phylogenomic Support.

The scale of data sets used to infer phylogenies has grown dramatically in the last decades, providing researchers with an enormous amount of information with which to draw inferences about evolutionary history. However, standard approaches to assessing confidence in those inferences (e.g., nonparametric bootstrap proportions [BP] and Bayesian posterior probabilities [PPs]) are still deeply influenced by statistical procedures and frameworks that were developed when information was much more limited. These approaches largely quantify uncertainty caused by limited amounts of data, which is often vanishingly small with modern, genome-scale sequence data sets. As a consequence, today’s phylogenomic studies routinely report near-complete confidence in their inferences, even when different studies reach strongly conflicting conclusions and the sites and loci in a single data set contain much more heterogeneity than our methods assume or can accommodate. Therefore, we argue that BPs and marginal PPs of bipartitions have outlived their utility as the primary means of measuring phylogenetic support for modern phylogenomic data sets with large numbers of sites relative to the number of taxa. Continuing to rely on these measures will hinder progress towards understanding remaining sources of uncertainty in the most challenging portions of the Tree of Life. Instead, we encourage researchers to examine the ideas and methods presented in this special issue of Systematic Biology and to explore the area further in their own work. The papers in this special issue outline strategies for assessing confidence and uncertainty in phylogenomic data sets that move beyond stochastic error due to limited data and offer promise for more productive dialogue about the challenges that we face in reaching our shared goal of understanding the history of life on Earth.[Big data; gene tree variation; genomic era; statistical bias.]

Propagating clade and model uncertainty to confidence intervals of divergence times and branch lengths

Confidence intervals of divergence times and branch lengths do not reflect uncertainty about their clades or about the prior distributions and other model assumptions on which they are based. Uncertainty about the clade may be propagated to a confidence interval by multiplying its confidence level by the bootstrap proportion of its clade or by another probability that the clade is correct. (If the confidence level is 95% and the bootstrap proportion is 90%, then the uncertainty-adjusted confidence level is (0.95)(0.90) = 86%.) Uncertainty about the model can be propagated to the confidence interval by reporting the union of the confidence intervals from all the plausible models. Unless there is no overlap between the confidence intervals, that results in an uncertainty-adjusted interval that has as its lower and upper limits the most extreme limits of the models. The proposed methods of uncertainty quantification may be used together.

Anchored hybrid enrichment challenges the traditional classification of flesh flies (Diptera: Sarcophagidae)

AbstractSarcophagidae is one of the most species‐rich families within the superfamily Oestroidea. This diversity is usually represented by three lineages: Miltogramminae, Paramacronychiinae and Sarcophaginae. Historically, the phylogenetic relationships among these lineages have been elusive, due to poorly supported hypotheses or small taxon sets, or both. This study provides a dramatic increase in molecular data, more balanced sampling of all three lineages from all biogeographical regions and a reassessment of morphological characters using scanning electron microscopy in the most comprehensive assessment of subfamily‐level phylogeny in Sarcophagidae to date. This analysis of the largest molecular dataset ever produced for a phylogenetic analysis of a fly lineage, with 950 loci from anchored hybrid enrichment comprising 435 930 bp from 101 species, revealed Paramacronychiinae as sister to Miltogramminae, not to Sarcophaginae, as suggested by adult morphology. Maximum likelihood analysis produced a well‐supported topology, with 91% of the nodes receiving strong bootstrap proportions (&gt; 97%). In contrast to the molecular data, three out of nine morphological characters studied point to a sister‐group relationship of (Sarcophaginae + Paramacronychiinae) and the remaining six characters are either silent on subfamily relationships or in need of further study. Re‐examination of morphological structures provides new insights into the evolution of male genitalic traits within Sarcophagidae and highlights their convergence producing conflicting phylogenetic signal. Our phylogeny reconciles older and widely used systems of classification with tree‐based thinking and sets up a classification of flesh flies that is more aligned with their evolutionary history.

Molecular Assessment of Anaplasma marginale in Bovine and Rhipicephalus (Boophilus) microplus Tick of Endemic Tribal Belt of Coastal South Gujarat, India.

This study described the detection, prevalence and phylogeny of Anaplasma marginale in the bovine (cattle and buffaloes) and Rhipicephalus (Boophilus) microplus tick belonged to the tribal area of coastal South Gujarat, India, by amplifying 576bp of major surface protein (msp) 5 gene using custom designed primers in the polymerase chain reaction (PCR). The PCR detection limit was up to 20 parasites/µl of blood in sensitivity experiment, and observed 100% specificity against Trypanosoma evansi, Babesia bigemina and Theileria annulata. Prevalence rate of the A. marginale in the bovine (n = 211)) was 18.48% and 6.64% (p < 0.05) as per the PCR and Giemsa stained blood smear, respectively. Febrile animals (35%) observed significantly (p < 0.05) higher incidence rate than the non-febrile (14.62%). The amplified msp5 had single cut site for the EcoR1 enzyme, upon digestion yielded two fragments of 365 and 211bp on 1.0% agarose gel. The current sequence (KC811329) showed 100% homology and 1064 total score with the published nucleotide sequences of msp5 of A. marginale in the NCBI-BLAST study. Monophyletic relationship was observed with high bootstrap proportion (>76% in Neighbor-Joining/ Maximum Likelihood) between the current and published nucleotide sequences in the phylogeny. Twenty out of 39 A. marginale infected bovine recorded R. (B.) microplus on their body surface, out of which 18 had detected the infection. The rickettsia was in 55%, 65% and 25% of anterior half, posterior half and egg of tick, respectively. The test detected A. marginale in a carrier, pre-symptomatic and symptomatic vertebrate hosts (cattle and buffalo) and different body parts of the starved R. (B.) microplus including its egg. The current genotype could be an explanation for the frequent outbreaks of bovine anaplasmosis in the targeted areas.

Bootstrap and Rogue Identification Tests for Phylogenetic Analyses.

Most phylogenetic tree-generating programs produce a fully dichotomous phylogenetic tree. However, as different markers may produce distinct topologies for the same set of organisms, topological tests are used to estimate the statistical reliability of the clades. In this protocol, we provide step-by-step instructions on how to perform the widely used bootstrap test using MEGA. However, a single unstable lineage, also known as a rogue lineage, may decrease the bootstrap proportions in many branches of the tree. This occurs because rogue taxa tend to bounce between clades from one pseudo-replicate to the next, lowering bootstrap proportions for many correct clades. Thus, it is important to identify and exclude rogue taxa before initiating a final phylogenetic analysis; here, we provide this protocol using the RogueNaRok platform.

Complete mitochondrial genome of the blind vent crab Gandalfus puia (Crustacea: Bythograeidae) from the Tonga Arc

The brachyuran crab Gandalfus puia is a species endemic to the hydrothermal vent fields in the Tonga-Kermadec Arc. In order to understand G. puia at the genomic level, we sequenced its mitochondrial genome (mitogenome) and then compared to other bythograeids. The mitogenome is 15,548 bp in length and exhibits brachyuran-typical gene arrangement. Its protein-coding genes were very similar to other bythograeid species with respect to length, AT content and start and stop codons. Additionally, we compared the mitogenomes of Gandalfus and the closely related Austinograea. The inter-specific nucleotide divergence was 13.4% in Gandalfus and 13.7–14.0% in Austinograea. The inter-generic nucleotide divergence between Gandalfus and Austinograea was 16.3–19.7%. Based on the phylogenetic tree constructed with maximum likelihood and Bayesian inference, bythograeid crabs were recognized as the monophyletic taxon with the high supporting values (100% bootstrap proportions and 1.00 posterior probabilities). These results are useful for understanding the phylogenetic relationships and evolution of bythograeid crabs.

Impacts of Terraces on Phylogenetic Inference.

Terraces are sets of trees with precisely the same likelihood or parsimony score, which can be induced by missing sequences in partitioned multi-locus phylogenetic data matrices. The potentially large set of trees on a terrace can be characterized by enumeration algorithms or consensus methods that exploit the pattern of partial taxon coverage in the data, independent of the sequence data themselves. Terraces can add ambiguity and complexity to phylogenetic inference, particularly in settings where inference is already challenging: data sets with many taxa and relatively few loci. In this article we present five new findings about terraces and their impacts on phylogenetic inference. First, we clarify assumptions about partitioning scheme model parameters that are necessary for the existence of terraces. Second, we explore the dependence of terrace size on partitioning scheme and indicate how to find the partitioning scheme associated with the largest terrace containing a given tree. Third, we highlight the impact of terrace size on bootstrap estimates of confidence limits in clades, and characterize the surprising result that the bootstrap proportion for a clade, as it is usually calculated, can be entirely determined by the frequency of bipartitions on a terrace, with some bipartitions receiving high support even when incorrect. Fourth, we dissect some effects of prior distributions of edge lengths on the computed posterior probabilities of clades on terraces, to understand an example in which long edges "attract" each other in Bayesian inference. Fifth, we describe how assuming relationships between edge-lengths of different loci, as an attempt to avoid terraces, can also be problematic when taxon coverage is partial, specifically when heterotachy is present. Finally, we discuss strategies for remediation of some of these problems. One promising approach finds a minimal set of taxa which, when deleted from the data matrix, reduces the size of a terrace to a single tree.

Concatenated alignments and the case of the disappearing tree.

BackgroundAnalyzed individually, gene trees for a given taxon set tend to harbour incongruent or conflicting signals. One popular approach to deal with this circumstance is to use concatenated data. But especially in prokaryotes, where lateral gene transfer (LGT) is a natural mechanism of generating genetic diversity, there are open questions as to whether concatenation amplifies or averages phylogenetic signals residing in individual genes. Here we investigate concatenations of prokaryotic and eukaryotic datasets to investigate possible sources of incongruence in phylogenetic trees and to examine the level of overlap between individual and concatenated alignments.ResultsWe analyzed prokaryotic datasets comprising 248 invidual gene trees from 315 genomes at three taxonomic depths spanning gammaproteobacteria, proteobacteria, and prokaryotes (bacteria plus archaea), and eukaryotic datasets comprising 279 invidual gene trees from 85 genomes at two taxonomic depths: across plants-animals-fungi and within fungi. Consistent with previous findings, the branches in trees made from concatenated alignments are, in general, not supported by any of their underlying individual gene trees, even though the concatenation trees tend to possess high bootstrap proportions values. For the prokaryote data, this observation is independent of phylogenetic depth and sequence conservation. The eukaryotic data show much better agreement between concatenation and single gene trees. LGT frequencies in trees were estimated using established methods. Sequence length in individual alignments, but not sequence divergence, was found to correlate with the generation of branches that correspond to the concatenated tree.ConclusionsThe weak correspondence of concatenation trees with single gene trees gives rise to the question where the phylogenetic signal in concatenated trees is coming from. The eukaryote data reveals a better correspondence between individual and concatenation trees than the prokaryote data. The question of whether the lack of correspondence between individual genes and the concatenation tree in the prokaryotic data is due to LGT or phylogenetic artefacts remains unanswered. If LGT is the cause of incongruence between concatenation and individual trees, we would have expected to see greater degrees of incongruence for more divergent prokaryotic data sets, which was not observed, although estimated rates of LGT suggest that LGT is responsible for at least some of the observed incongruence.Electronic supplementary materialThe online version of this article (doi:10.1186/s12862-014-0266-0) contains supplementary material, which is available to authorized users.

Warneckea parvifolia (Melastomataceae–Olisbeoideae), a new “sand-forest” endemic from northeastern KwaZulu-Natal (South Africa) and southernmost Mozambique, and a phylogenetic analysis of eastern and southern African representatives of W. section Warneckea

Warneckea populations from “sand-forest” or “sand-thicket” habitats in Tembe Elephant Park, South Africa, and Licuati Forest Reserve in adjacent southern Mozambique were previously thought to be a small-leaved form of W. sousae, which typically includes larger-leaved plants ranging from central Mozambique northward to Tanzania. We examine this hypothesis using molecular and morphological evidence. Maximum-likelihood phylogenetic analysis of combined nrDNA ETS and ITS sequence data failed to resolve W. sousae and the Maputaland populations as an exclusively monophyletic group. Instead, the Kenyan endemic W. mouririifolia was strongly supported as the sister species of W. sousae, and the Maputaland plants were resolved in a separate, strongly supported clade together with populations of an as-yet undetermined Warneckea species from northern Mozambique. A hypothesis of exclusive monophyly for the plants from Tembe and Licuati had moderate support in separate ETS and ITS1 analyses (bootstrap proportions of 88% and 81%, respectively). Statistically significant differences in leaf dimensions and internode length were found between the Maputaland plants and typical W. sousae. We conclude that the populations from Tembe and Licuati represent a distinct species, which we describe as W. parvifolia. The species differs from W. sousae in having shorter internodes (mostly 5–25mm not 10–60mm long), smaller leaves (mostly 14–32×8–19mm not 40–76×22–52mm), shorter petioles (mostly 1–1.5mm not 1.5–6mm long), smaller flowers (hypanthium 1×1.5–1.75mm not 1.5–2×2mm; calyx lobes 0.5mm not 0.75mm long; staminal filaments 3–4mm not 5mm long; style 4–5mm not 9mm long), and globose fruit (not obovoid). An IUCN conservation status of Endangered (EN) B1a, b(ii, iii) is indicated for W. parvifolia, due to its limited distribution and projected declines in its habitat quality and area of occupancy.

How Many SNPs Should Be Used for the Human Phylogeny of Highly Related Ethnicities? A Case of Pan Asian 63 Ethnicities

In planning a model-based phylogenic study for highly related ethnic data, the SNP marker number is an important factor to determine for relationship inferences. Genotype frequency data, utilizing a sub sampling method, from 63 Pan Asian ethnic groups was used for determining the minimum SNP number required to establish such relationships. Bootstrap random sub-samplings were done from 5.6K PASNPi SNP data. DA distance was calculated and neighbour-joining trees were drawn with every re-sampling data set. Consensus trees were made with the same 100 sub-samples and bootstrap proportions were calculated. The tree consistency to the one obtained from the whole marker set, improved with increasing marker numbers. The bootstrap proportions became reliable when more than 7,000 SNPs were used at a time. Within highly related ethnic groups, the minimum SNPs number for a robust neighbor-joining tree inference was about 7,000 for a 95% bootstrap support.

Morphology and phylogenetic relationships of a remarkable new genus and two new species of Neotropical freshwater stingrays from the Amazon basin (Chondrichthyes: Potamotrygonidae)

The morphology and phylogenetic relationships of a new genus and two new species of Neotropical freshwater stingrays, family Potamotrygonidae, are investigated and described in detail. The new genus, Heliotrygon, n. gen., and its two new species, Heliotrygon gomesi, n. sp. (type-species) and Heliotrygon rosai, n. sp., are compared to all genera and species of potamotrygonids, based on revisions in progress. Some of the derived features of Heliotrygon include its unique disc proportions (disc highly circular, convex anteriorly at snout region, its width and length very similar), extreme subdivision of suborbital canal (forming a complex honeycomb-like pattern anterolaterally on disc), stout and triangular pelvic girdle, extremely reduced caudal sting, basibranchial copula with very slender and acute anterior extension, and precerebral and frontoparietal fontanellae of about equal width, tapering very little posteriorly. Both new species can be distinguished by their unique color patterns: Heliotrygon gomesi is uniform gray to light tan or brownish dorsally, without distinct patterns, whereas Heliotrygon rosai is characterized by numerous white to creamy-white vermiculate markings over a light brown, tan or gray background color. Additional proportional characters that may further distinguish both species are also discussed. Morphological descriptions are provided for dermal denticles, ventral lateral-line canals, skeleton, and cranial, hyoid and mandibular muscles of Heliotrygon, which clearly corroborate it as the sister group of Paratrygon. Both genera share numerous derived features of the ventral lateral-line canals, neurocranium, scapulocoracoid, pectoral basals, clasper morphology, and specific patterns of the adductor mandibulae and spiracularis medialis muscles. Potamotrygon and Plesiotrygon are demonstrated to share derived characters of their ventral lateral-line canals, in addition to the presence of angular cartilages. Our morphological phylogeny is further corroborated by a molecular phylogenetic analysis of cytochrome b based on four sequences (637 base pairs in length), representing two distinct haplotypes for Heliotrygon gomesi. Parsimony analysis produced a single most parsimonious tree revealing Heliotrygon and Paratrygon as sister taxa (bootstrap proportion of 70%), which together are the sister group to a clade including Plesiotrygon and species of Potamotrygon. These unusual stingrays highlight that potamotrygonid diversity, both in terms of species composition and undetected morphological and molecular patterns, is still poorly known.

LifePrint: a novel k-tuple distance method for construction of phylogenetic trees

PurposeHere we describe LifePrint, a sequence alignment-independent k-tuple distance method to estimate relatedness between complete genomes.MethodsWe designed a representative sample of all possible DNA tuples of length 9 (9-tuples). The final sample comprises 1878 tuples (called the LifePrint set of 9-tuples; LPS9) that are distinct from each other by at least two internal and noncontiguous nucleotide differences. For validation of our k-tuple distance method, we analyzed several real and simulated viroid genomes. Using different distance metrics, we scrutinized diverse viroid genomes to estimate the k-tuple distances between these genomic sequences. Then we used the estimated genomic k-tuple distances to construct phylogenetic trees using the neighbor-joining algorithm. A comparison of the accuracy of LPS9 and the previously reported 5-tuple method was made using symmetric differences between the trees estimated from each method and a simulated “true” phylogenetic tree.ResultsThe identified optimal search scheme for LPS9 allows only up to two nucleotide differences between each 9-tuple and the scrutinized genome. Similarity search results of simulated viroid genomes indicate that, in most cases, LPS9 is able to detect single-base substitutions between genomes efficiently. Analysis of simulated genomic variants with a high proportion of base substitutions indicates that LPS9 is able to discern relationships between genomic variants with up to 40% of nucleotide substitution.ConclusionOur LPS9 method generates more accurate phylogenetic reconstructions than the previously proposed 5-tuples strategy. LPS9-reconstructed trees show higher bootstrap proportion values than distance trees derived from the 5-tuple method.

Looks can deceive: Molecular phylogeny of a family of flatworm ectoparasites (Monogenea: Capsalidae) does not reflect current morphological classification

The morphological based taxonomy of highly derived parasite groups is likely to poorly reflect their evolutionary relationships. The taxonomy of the monogenean family Capsalidae, which comprises approximately 180 species of flatworm parasites that predominantly attach to external surfaces of chondrichthyan and teleost fishes, is based mainly on six morphological characters. The phylogenetic history of the family is largely unknown. We reconstructed the phylogenetic relationships of 47 species in 20 genera from eight of the nine subfamilies, from nucleotide sequences of three unlinked nuclear genes, 28S ribosomal RNA, Histone 3 and Elongation Factor 1 α. Our phylogeny was well corroborated, with 75% of branches receiving strong support from both Bayesian posterior probabilities and maximum likelihood bootstrap proportions and all nodes showed positive partitioned likelihood support for each of the three genes. We found that the family was monophyletic, with the Gyrodactylidae and Udonellidae forming the sister group. The Capsalinae was monophyletic, however, our data do not support monophyly for the Benedeniinae, Entobdellinae and Trochopodinae. Monophyly was supported for Capsala, Entobdella, Listrocephalos, Neobenedenia and Tristoma, but Benedenia and Neoentobdella were polyphyletic. Comparisons of the distribution of character states for the small number of morphological characters on the molecular phylogeny show a high frequency of apparent homoplasy. Consequently the current morphological classification shows little correspondence with the phylogenetic relationships within the family.

Homoplasy and Clade Support

Distinguishing phylogenetic signal from homoplasy (shared similarities among taxa that do not arise by common ancestry) is an implicit goal of any phylogenetic study. Large amounts of homoplasy can interfere with accurate tree inference, and it is expected that common measures of clade support, including bootstrap proportions and Bayesian posterior probabilities, should also be impacted to some degree by homoplasy. Through data simulation and analysis of 38 empirical data sets, we show that high amounts of homoplasy will affect all measures of clade support in a manner that is dependent on clade size. More specifically, the smallest taxon bipartitions in an unrooted tree topology will receive higher support relative to clades of intermediate sizes, even when all clades are supported by the same amount of data. We determine that the ultimate causes of this effect are the inclusion of random trees (due to homoplasy) during bootstrap resampling and Markov chain Monte Carlo (MCMC) topology searching and the higher relative proportion of small taxon bipartitions (i.e., 2 or 3 taxa) to larger sized bipartitions. However, the use of explicit model-based methods, especially Bayesian MCMC methods, effectively overcomes this clade size effect even when very small amounts of phylogenetic signal are present. We develop a post hoc statistic, the clade disparity index (CDI), to measure both the relative magnitude of the clade size effect and its statistical significance. In analyses of both simulated and empirical data, CDI values indicate that Bayesian MCMC analyses are substantially more likely to estimate clade support values that are uncorrelated with clade size than are maximum parsimony and maximum likelihood bootstrap analyses and thus less affected by homoplasy. These results may be especially relevant to "deep" phylogenetic problems, such as reconstructing the tree of life, as they represent the largest possible extremes of time and evolutionary rates, 2 factors that cause homoplasy.

Formalizing reliability in the taxonomic congruence approach

In the ‘total evidence’ approach to phylogenetics, the reliability of a clade is implicitly measured by its degree of support, often embodied in a robustness index such as a bootstrap proportion. In the taxonomic congruence approach, the measurement of reliability has been implemented by various consensus or supertree methods, but was seldom explicitly discussed as such. We explore a reliability index for clades using their repetition across independent data sets. All possible combinations of the elementary data sets are used to compose the sets of independent data sets, across which the repetitions are counted. The more a clade occurs across such independent combinations, the higher its index. However, if other repeated clades occur that are incompatible with that clade, its index is decreased to take into account the uncertainty resulting from conflicting hypotheses. Results can be summarized through a greedy consensus tree in which clades appear according to their repetition indices. This index is tested on a 73 acanthomorph taxa data set composed of five independent molecular markers and multiple combinations of them. On this particular application, we confirm that reliability as defined here and robustness (estimated by bootstrap proportions obtained from a ‘total evidence’ approach) should be clearly distinguished.