Brooks Parsimony Analysis: a valiant failure

A recent comparison of two methods for examining correlated host and parasite phylogenies, namely TreeMap 1.0 and Brooks Parsimony Analysis concluded that the latter method performed better and is to be preferred. Reevaluation of the examples contrived for that study demonstrates that the two methods were only compared on one kind of problem (widespread parasite) for which there is an easy fix for TreeMap 1.0. Other kinds of problems like host‐switching among sister taxa or host‐switching over great distances across a host tree befuddle BPA even as they are readily resolved parsimoniously by TreeMap 1.0. These difficulties, compounded with inaccurate counting of ad hoc hypotheses required by its solutions render BPA unsuitable for comparison of host and parasite phylogenies.     

The Parsimony Ratchet, a New Method for Rapid Parsimony Analysis

The Parsimony Ratchet ¹ 1 This method, the Parsimony Ratchet, was originally presented at the Numerical Cladistics Symposium at the American Museum of Natural History, New York, in May 1998 (see Horovitz, 1999) and at the Meeting of the Willi Hennig Society (Hennig XVII) in September 1998 in São Paulo, Brazil.
is presented as a new method for analysis of large data sets. The method can be easily implemented with existing phylogenetic software by generating batch command files. Such an approach has been implemented in the programs DADA (Nixon, 1998) and Winclada (Nixon, 1999). The Parsimony Ratchet has also been implemented in the most recent versions of NONA (Goloboff, 1998). These implementations of the ratchet use the following steps: (1) Generate a starting tree (e.g., a “Wagner” tree followed by some level of branch swapping or not). (2) Randomly select a subset of characters, each of which is given additional weight (e.g., add 1 to the weight of each selected character). (3) Perform branch swapping (e.g., “branch-breaking” or TBR) on the current tree using the reweighted matrix, keeping only one (or few) tree. (4) Set all weights for the characters to the “original” weights (typically, equal weights). (5) Perform branch swapping (e.g., branch-breaking or TBR) on the current tree (from step 3) holding one (or few) tree. (6) Return to step 2. Steps 2–6 are considered to be one iteration, and typically, 50–200 or more iterations are performed. The number of characters to be sampled for reweighting in step 2 is determined by the user; I have found that between 5 and 25% of the characters provide good results in most cases. The performance of the ratchet for large data sets is outstanding, and the results of analyses of the 500 taxon seed plant rbcL data set (Chase et al., 1993) are presented here. A separate analysis of a three-gene data set for 567 taxa will be presented elsewhere (Soltis et al., in preparation) demonstrating the same extraordinary power. With the 500-taxon data set, shortest trees are typically found within 22 h (four runs of 200 iterations) on a 200-MHz Pentium Pro. These analyses indicate efficiency increases of 20×–80× over “traditional methods” such as varying taxon order randomly and holding few trees, followed by more complete analyses of the best trees found, and thousands of times faster than nonstrategic searches with PAUP. Because the ratchet samples many tree islands with fewer trees from each island, it provides much more accurate estimates of the “true” consensus than collecting many trees from few islands. With the ratchet, Goloboff's NONA, and existing computer hardware, data sets that were previously intractable or required months or years of analysis with PAUP* can now be adequately analyzed in a few hours or days.     

On the Three-Taxon Approach to Parsimony Analysis

The following three basic defects for which three-taxon analysis has been rejected as a method for biological systematics are reviewed: (1) character evolution is a priori assumed to be irreversible; (2) basic statements that are not logically independent are treated as if they are; (3) three-taxon statements that are considered as independent support for a given tree may be mutually exclusive on that tree. It is argued that these criticisms only relate to the particular way the three-taxon approach was originally implemented. Four-taxon analysis, an alternative implementation that circumvents these problems, is derived. Four-taxon analysis is identical to standard parsimony analysis except for an unnatural restriction on the maximum amount of homoplasy that may be concentrated in a single character state. This restriction follows directly from the basic tenet of the three-taxon approach, that character state distributions should be decomposed into basic statements that are, in themselves, still informative with respect to relationships. A reconsideration of what constitutes an elementary relevant statement in systematics leads to a reformulation of standard parsimony as two-taxon analysis and to a rejection of four-taxon analysis as a method for biological systematics.     

Ending a decade of deception: a valiant failure, a not-so-valiant failure, and a success story

Prior studies involving two methods, Brooks Parsimony Analysis (BPA) and TreeMap, have found BPA to be the more reliable method. Recent criticisms leveled at these studies argue that the tests were unfairly created and biased in favor of BPA. The authors of a recent critique offered new exemplars to demonstrate flaws in BPA, plus a simple fix to correct the flaws found in TreeMap. A re‐evaluation of their exemplars clearly shows that the authors' calculations are incorrect, their understanding of the methods is lacking, and that their simple fix does not work. Additional analyses using TreeMap 2.02 are run to show that TreeMap 2.02, like TreeMap 1.0, cannot adequately deal with widespread parasites, contrary to the claims of its supporters. Furthermore, the exemplars corroborate previous findings that BPA, when calculated correctly, is more reliable than TreeMap1.0 and TreeMap 2.02 and therefore the method of choice in coevolutionary and biogeographic studies.     

Innovative combination strategy to enhance effect and diminish adverse effects of glucocorticoids: another promise?

In a paper by Zimmermann and colleagues in this issue of Arthritis Research & Therapy, results of extended laboratory research with the drug combination of prednisolone and dipyridamole are reported. There seems to be a boost and extension of the glucocorticoid effect by the combination, without a clear increase of adverse effects, potentially allowing the application of lower dosages. However, laboratory models are not patients and the glucocorticoid mechanisms leading to effects and adverse effects are manifold. The next required step will be to demonstrate the improved therapeutic window in patients in adequate comparative clinical trials, assessing predefined beneficial effects and adverse effects in a standardized way.     

Parsimony Analysis as a Specific Kind of Homology Estimation and the Implications for Character Weighting

“Remane-Hennigian systematists” still reject parsimony analysis for phylogenetics, because homology or apomorphy analyses are not included. In contrast, “pattern cladists” regard homology as a deductive concept after applying a parsimony test of character congruence. However, as in molecular phylogeny, selection of “good” characters is always done on the basis of ana priorihomology analysis. The distribution criterion of homology—“homologous characters have identical or hierarchical distribution”—is the basis of parsimony analysis. Because this criterion also might fail in cases of genealogical reticulation or concerted homoplasy, character congruence is not a strict test but another probabilistic criterion of homology. A synthetic approach is proposed for phenotypic analysis with application ofa prioricriteria of homology. The resultinga prioriprobabilities of homology serve as criteria for selection and weighting of characters (very low = not selected/poor/mediocre/good/Dollo characters). After application of a parsimony algorithm the final cladogram decides homology estimations.     

Analysis of the Distribution of Bootstrap Tree Lengths Using the Maximum Parsimony Method

The bootstrap is an important tool for estimating the confidence interval of monophyletic groups within phylogenies. Although bootstrap analyses are used in most evolutionary studies, there is no clear consensus as how best to interpret bootstrap probability values. To study further the bootstrap method, nine small subunit ribosomal DNA (SSU rDNA) data sets were submitted to bootstrapped maximum parsimony (MP) analyses using unweighted and weighted sequence positions. Analyses of the lengths (i.e., parsimony steps) of the bootstrap trees show that the shape and mean of the bootstrap tree distribution may provide important insights into the evolutionary signal within the sequence data. With complex phylogenies containing nodes defined by short internal branches (multifurcations), the mean of the bootstrap tree distribution may differ by 2 standard deviations from the length of the best tree found from the original data set. Weighting sequence positions significantly increases the bootstrap values at internal nodes. There may, however, be strong bootstrap support for conflicting species groupings among different data sets. This phenomenon appears to result from a correlation between the topology of the tree used to create the weights and the topology of the bootstrap consensus tree inferred from the MP analysis of these weighted data. The analyses also show that characteristics of the bootstrap tree distribution (e.g., skewness) may be used to choose between alternative weighting schemes for phylogenetic analyses.     

Achievement of genetics in plant reproduction research: the past decade for the coming decade

In the last decade, a variety of innovations of emerging technologies in science have been accomplished. Advanced research environment in plant science has made it possible to obtain whole genome sequence in plant species. But now we recognize this by itself is not sufficient to understand the overall biological significance. Since Gregor Mendel established a principle of genetics, known as Mendel's Laws of Inheritance, genetics plays a prominent role in life science, and this aspect is indispensable even in modern plant biology. In this review, we focus on achievements of genetics on plant sexual reproduction research in the last decade and discuss the role of genetics for the coming decade. It is our hope that this will shed light on the importance of genetics in plant biology and provide valuable information to plant biologists.     

Asymmetry of the os pubis: Implications for the Suchey‐Brooks method

Studies of skeletal development frequently document populational incidences of bilateral asymmetry. Degenerative morphological skeletal changes, attributed to age related and irregular ossification, may also progress asymmetrically, either as the result of asymmetric biomechanical factors expressed over the lifespan, asymmetric expression of physiological processes, or progressive magnification of asymmetry acquired previously during development. This study illustrates the effects of bilateral asymmetry on age at death estimates obtained from human skeletal remains. The Suchey‐Brooks method, which uses the pubic symphyseal face for age estimation (Katz and Suchey, Am J Phys Anthropol 69 1986 427–435), was selected for the study based on its widespread use. Asymmetry in the Suchey‐Brooks symphyseal age phases was found in over 60% of a sample composed of 20th century White male individuals from 18 to 86 years of age (N = 130). However, results suggest that the presence of asymmetry does not compromise the accuracy of the Suchey‐Brooks method if the morphologically older symphyseal face of an asymmetric individual is used to estimate age at death. In addition, weak directional asymmetry and a correlation between age and asymmetry were found. This suggests that a comparison of asymmetry in this area with that in other skeletal areas, where the factors originating and influencing asymmetry are better understood, may be useful in better understanding the biological processes which underlie the age markers used in the Suchey‐Brooks method. Am J Phys Anthropol 2009. © 2009 Wiley‐Liss, Inc.     

Parsimony analysis of unaligned sequence data: some clarifications

De Laet (2015) claimed that minimization of ad hoc hypotheses of homoplasy does not lead to a preference for trivial optimizations when analysing unaligned sequence data, as claimed by Wheeler (2012; see also Kluge and Grant, 2006). In addition, De Laet asserted that Kluge and Grant's (2006) parsimony rationale is internally inconsistent in terms of Baker's (2003) theoretical framework. We argue that De Laet used extraneous presuppositions to critique Wheeler's position and, as such, his criticism should be considered cautiously in terms of its scope. Finally, we demonstrate that considering Kluge and Grant's parsimony rationale as inconsistent rests on De Laet's misunderstanding of the ideographic character concept and the consequences of relating it to Baker's rationale.     

Harvesting the genomic promise: recombineering sequences for phenotypes

The past decade has witnessed the construction of linkage and physical maps defining quantitative trait loci (QTL) in various domesticated species. Targeted chromosomal regions are being further characterized through the construction of bacterial artificial chromosome (BAC) contigs in order to isolate and characterize genes contributing towards phenotypic variation. Whole-genome BAC contigs are also being constructed that will serve as the tiling path for genomic sequencing. Harvesting this genetic information for biological gain requires either genetic selection or the production of genetically modified animals. This later approach when coupled with nuclear transfer technology (NT) provides "clones" of genetically modified animals. However, to date, the production of genetically modified animals has been limited to either microinjection of small gene constructs into embryos with random insertion or complex gene constructs designed to knock-out targeted gene expression. Neither of these approaches provides for introducing directed genetic manipulation allowing for allelic substitution [knock-in], subsequent analyses of gene expression, and cloning. An alternative approach utilizing genomic sequence information and recombineering to direct gene targeting of specific porcine BACs is presented here.     

Using Genes as Characters and a Parsimony Analysis to Explore the Phylogenetic Position of Turtles

The phylogenetic position of turtles within the vertebrate tree of life remains controversial. Conflicting conclusions from different studies are likely a consequence of systematic error in the tree construction process, rather than random error from small amounts of data. Using genomic data, we evaluate the phylogenetic position of turtles with both conventional concatenated data analysis and a “genes as characters” approach. Two datasets were constructed, one with seven species (human, opossum, zebra finch, chicken, green anole, Chinese pond turtle, and western clawed frog) and 4584 orthologous genes, and the second with four additional species (soft-shelled turtle, Nile crocodile, royal python, and tuatara) but only 1638 genes. Our concatenated data analysis strongly supported turtle as the sister-group to archosaurs (the archosaur hypothesis), similar to several recent genomic data based studies using similar methods. When using genes as characters and gene trees as character-state trees with equal weighting for each gene, however, our parsimony analysis suggested that turtles are possibly sister-group to diapsids, archosaurs, or lepidosaurs. None of these resolutions were strongly supported by bootstraps. Furthermore, our incongruence analysis clearly demonstrated that there is a large amount of inconsistency among genes and most of the conflict relates to the placement of turtles. We conclude that the uncertain placement of turtles is a reflection of the true state of nature. Concatenated data analysis of large and heterogeneous datasets likely suffers from systematic error and over-estimates of confidence as a consequence of a large number of characters. Using genes as characters offers an alternative for phylogenomic analysis. It has potential to reduce systematic error, such as data heterogeneity and long-branch attraction, and it can also avoid problems associated with computation time and model selection. Finally, treating genes as characters provides a convenient method for examining gene and genome evolution.     

Systems vaccinology: a promise for the young and the poor

As a child, the risk of suffering and dying from infection is higher the younger you are; and higher, the less developed a region you are born in. Childhood vaccination programmes have greatly reduced mortality around the world, but least so for the very young among the very poor of the world. This appears partly owing to suboptimal vaccine effectiveness. Unfortunately, although most vaccines are administered to the newborn and very young infant (less than or equal to two months), we know the least about their host response to vaccination. We thus currently lack the knowledge to guide efforts aimed at improving vaccine effectiveness in this vulnerable population. Systems vaccinology, the study of molecular networks activated by immunization, has begun to provide unprecedented insights into mechanisms leading to vaccine-induced protection from infection or disease. However, all published reports of systems vaccinology have focused on either adults or at most children and older infants, not those most in need, i.e. newborns and very young infants. Given that the tools of systems vaccinology work perfectly well with very small sample volumes, it is time we deliver the promise that systems vaccinology holds for those most in need of vaccine-mediated protection from infection.     

The Wagner Parsimony using morphological characters: a new method for palaeosynecological studies

The limits and difficulties related to the tools currently in use for palaeosynecological comparisons of faunas or floras of different geological periods are discussed. The new method of the Wagner parsimony Applied to Palaeosynecology Using Morphology (WAPUM method), is defined and tested on morphological characters gathered from two insect groups Odonatoptera and Thripida. The difficulties related to the monophyly of the taxonomic groups used in the more traditional approaches are no longer a problem when using the WAPUM method. In the WAPUM a character is ‘presence versus absence of species bearing a morphological structure’. The results obtained from use of the WAPUM minimize the number of changes among character states. Application of the WAPUM could reveal signals to confirm or object the currently available scenarios for the global changes in the evolution of past diversity and disparity of organisms (major changes or global crises of diversity).