SKEWNESS AND PERMUTATION

Maximum Parsimony and the Skewness Test: A Simulation Study of the Limits of Applicability

The maximum parsimony (MP) method for inferring phylogenies is widely used, but little is known about its limitations in non-asymptotic situations. This study employs large-scale computations with simulated phylogenetic data to estimate the probability that MP succeeds in finding the true phylogeny for up to twelve taxa and 256 characters. The set of candidate phylogenies are taken to be unrooted binary trees; for each simulated data set, the tree lengths of all (2n - 5)!! candidates are computed to evaluate quantities related to the performance of MP, such as the probability of finding the true phylogeny, the probability that the tree with the shortest length is unique, the probability that the true phylogeny has the shortest tree length, and the expected inverse of the number of trees sharing the shortest length. The tree length distributions are also used to evaluate and extend the skewness test of Hillis for distinguishing between random and phylogenetic data. The results indicate, for example, that the critical point after which MP achieves a success probability of at least 0.9 is roughly around 128 characters. The skewness test is found to perform well on simulated data and the study extends its scope to up to twelve taxa.

Bush Blitz aids description of three new species and a new genus of Australian beeflies (Diptera, Bombyliidae, Exoprosopini)

Bush Blitz is a three-year multimillion dollar program to document the plants and animals in hundreds of properties across Australia's National Reserve System. The core focus is on nature discovery - identifying and describing new species of plants and animals. The Bush Blitz program has enabled the collection and description of beeflies (Diptera, Bombyliidae) from surveys in Western Australia and Queensland. Three new species of Australian beeflies belonging to the Exoprosopini are described; Palirika mackenziei Lambkin sp. n., Palirika culgoafloodplainensis Lambkin sp. n., and Larrpana bushblitz Lambkin sp. n. Phylogenetic analysis of 40 Australian exoprosopine species belonging to the Balaana generic-group Lambkin & Yeates 2003 supports the placement of the three new species into existing genera, and the erection and description of the new genus Ngalki Lambkin gen. n. for Ngalki trigonium (Lambkin & Yeates 2003) comb. n. Revised keys are provided for the genera of the Australian Balaana genus-group and the species of Palirika Lambkin & Yeates, 2003 and Larrpana Lambkin & Yeates 2003. With the description of the three new species and the transferral of Munjua trigona Lambkin & Yeates 2003 into the new genus Ngalki Lambkin gen. n., three genera are rediagnosed; Munjua Lambkin & Yeates 2003, Palirika and Larrpana.

The monogenean parasite fauna of cichlids: a potential tool for host biogeography

We discuss geographical distribution and phylogeny of Dactylogyridea (Monogenea) parasitizing Cichlidae to elucidate their hosts' history. Although mesoparasitic Monogenea (Enterogyrus spp.) show typical vicariant distribution, ectoparasitic representatives from different continents are not considered sister taxa, hence their distribution cannot result from vicariance alone. Because of the close host-parasite relationship, this might indicate that present-day cichlid distribution may also reflect dispersal through coastal or brackish waters. Loss of ectoparasites during transoceanic migration, followed by lateral transfer from other fish families might explain extant host-parasite associations. Because of its mesoparasitic nature, hence not subject to salinity variations of the host's environment, Enterogyrus could have survived marine migrations, intolerable for ectoparasites. Host-switches and salinity transitions may be invoked to explain the pattern revealed by a preliminary morphological phylogeny of monogenean genera from Cichlidae and other selected Monogenea genera, rendering the parasite distribution explicable under both vicariance and dispersal. Testable hypotheses are put forward in this parasitological approach to cichlid biogeography. Along with more comprehensive in-depth morphological phylogeny, comparison with molecular data, clarifying dactylogyridean evolution on different continents and from various fish families, and providing temporal information on host-parasite history, are needed to discriminate between the possible scenarios.

Species delimitation in downy mildews: the case of Hyaloperonospora in the light of nuclear ribosomal ITS and LSU sequences

Species definitions for plant pathogens have considerable practical impact for measures such as plant protection or biological control, and are also important for comparative studies involving model organisms. However, in many groups, the delimitation of species is a notoriously difficult taxonomic problem. This is particularly evident in the obligate biotrophic downy mildew genera (Peronosporaceae, Peronosporales, Oomycetes), which display a considerable diversity with respect to genetic distances and host plants, but are, for the most part, morphologically rather uniform. The recently established genus Hyaloperonospora is of particular biological interest because it shows an impressive radiation on virtually a single host family, Brassicaceae, and it contains the downy mildew parasite, Arabidopsis thaliana, of importance as a model organism. Based on the most comprehensive molecular sampling of specimens from a downy mildew genus to date, including various collections from different host species and geographic locations, we investigate the phylogenetic relationships of Hyaloperonospora by molecular analysis of the nuclear ribosomal ITS and LSU sequences. Phylogenetic trees were inferred with ML and MP from the combined dataset; partitioned Bremer support (PBrS) was used to assess potential conflict between data partitions. As in other downy mildew groups, the molecular data clearly corroborate earlier results that supported the use of narrow species delimitations and host ranges as taxonomic markers. With few exceptions, suggested species boundaries are supported without conflict between different data partitions. The results indicate that a combination of molecular and host features is a reliable means to discriminate downy mildew species for which morphological differences are unknown.

Higher-order phylogeny of modern birds (Theropoda, Aves: Neornithes) based on comparative anatomy. II. Analysis and discussion

In recent years, avian systematics has been characterized by a diminished reliance on morphological cladistics of modern taxa, intensive palaeornithogical research stimulated by new discoveries and an inundation by analyses based on DNA sequences. Unfortunately, in contrast to significant insights into basal origins, the broad picture of neornithine phylogeny remains largely unresolved. Morphological studies have emphasized characters of use in palaeontological contexts. Molecular studies, following disillusionment with the pioneering, but non-cladistic, work of Sibley and Ahlquist, have differed markedly from each other and from morphological works in both methods and findings. Consequently, at the turn of the millennium, points of robust agreement among schools concerning higher-order neornithine phylogeny have been limited to the two basalmost and several mid-level, primary groups. This paper describes a phylogenetic (cladistic) analysis of 150 taxa of Neornithes, including exemplars from all non-passeriform families, and subordinal representatives of Passeriformes. Thirty-five outgroup taxa encompassing Crocodylia, predominately theropod Dinosauria, and selected Mesozoic birds were used to root the trees. Based on study of specimens and the literature, 2954 morphological characters were defined; these characters have been described in a companion work, approximately one-third of which were multistate (i.e. comprised at least three states), and states within more than one-half of these multistate characters were ordered for analysis. Complete heuristic searches using 10 000 random-addition replicates recovered a total solution set of 97 well-resolved, most-parsimonious trees (MPTs). The set of MPTs was confirmed by an expanded heuristic search based on 10 000 random-addition replicates and a full ratchet-augmented exploration to ascertain global optima. A strict consensus tree of MPTs included only six trichotomies, i.e. nodes differing topologically among MPTs. Bootstrapping (based on 10 000 replicates) percentages and ratchet-minimized support (Bremer) indices indicated most nodes to be robust. Several fossil Neornithes (e.g. Dinornithiformes, Aepyornithiformes) were placed within the ingroup a posteriori either through unconstrained, heursitic searches based on the complete matrix augmented by these taxa separately or using backbone-constraints. Analysis confirmed the topology among outgroup Theropoda and achieved robust resolution at virtually all levels of the Neornithes. Findings included monophyly of the palaeognathous birds, comprising the sister taxa Tinamiformes and ratites, respectively, and the Anseriformes and Galliformes as monophyletic sister-groups, together forming the sister-group to other Neornithes exclusive of the Palaeognathae (Neoaves). Noteworthy inferences include: (i) the sister-group to remaining Neoaves comprises a diversity of marine and wading birds; (ii) Podicipedidae are the sister-group of Gaviidae, and not closely related to the Phoenicopteridae, as recently suggested; (iii) the traditional Pelecaniformes, including the shoebill (Balaeniceps rex) as sister-taxon to other members, are monophyletic; (iv) traditional Ciconiiformes are monophyletic; (v) Strigiformes and Falconiformes are sister-groups; (vi) Cathartidae is the sister-group of the remaining Falconiformes; (vii) Ralliformes (Rallidae and Heliornithidae) are the sister-group to the monophyletic Charadriiformes, with the traditionally composed Gruiformes and Turniciformes (Turnicidae and Mesitornithidae) sequentially paraphyletic to the entire foregoing clade; (viii) Opisthocomus hoazin is the sister-taxon to the Cuculiformes (including the Musophagidae); (ix) traditional Caprimulgiformes are monophyletic and the sister-group of the Apodiformes; (x) Trogoniformes are the sister-group of Coliiformes; (xi) Coraciiformes, Piciformes and Passeriformes are mutually monophyletic and closely related; and (xii) the Galbulae are retained within the Piciformes. Unresolved portions of the Neornithes (nodes having more than one most-parsimonious solution) comprised three parts of the tree: (a) several interfamilial nodes within the Charadriiformes; (b) a trichotomy comprising the (i) Psittaciformes, (ii) Columbiformes and (iii) Trogonomorphae (Trogoniformes, Coliiformes) + Passerimorphae (Coraciiformes, Piciformes, Passeriformes); and (c) a trichotomy comprising the Coraciiformes, Piciformes and Passeriformes. The remaining polytomies were among outgroups, although several of the highest-order nodes were only marginally supported; however, the majority of nodes were resolved and met or surpassed conventional standards of support. Quantitative comparisons with alternative hypotheses, examination of highly supportive and diagnostic characters for higher taxa, correspondences with prior studies, complementarity and philosophical differences with palaeontological phylogenetics, promises and challenges of palaeogeography and calibration of evolutionary rates of birds, and classes of promising evidence and future directions of study are reviewed. Homology, as applied to avian examples of apparent homologues, is considered in terms of recent theory, and a revised annotated classification of higher-order taxa of Neornithes and other closely related Theropoda is proposed. (c) 2007 The Linnean Society of London, Zoological Journal of the Linnean Society, 2007, 149, 1-95.

Probing the relationships of the branchiopod crustaceans

The Branchiopoda display extraordinary variation in body form, even within the morphologically diverse crustaceans. To fully understand the origin and evolution of these morphological reconfigurations, a robust phylogeny of the group is essential. To infer the affinities among branchiopods, we employed two approaches to taxon and gene sampling, presented new sequence data from three genes, incorporated previously published sequence data from three additional genes, and utilized comprehensive techniques of phylogeny reconstruction. The results provided support for a number of longstanding hypotheses concerning the relationships among the orders. For example, we obtained support for the Cladoceromorpha and Gymnomera, and favoured a unique arrangement of the cladoceran orders. A few affinities remain to be resolved, particularly at the base of the Phyllopoda and within the Anomopoda. However, the results suggest that increased gene sampling is recommended for future investigations of branchiopod systematics.

Support versus corroboration

Numerous metrics have been developed that attempt to assess the reliability of phylogenetic trees. Several of these commonly used measures of tree and tree branch support are described and discussed in the context of their relationship to Popperian corroboration. Claims that measures of support indicate the accuracy of phylogenetic trees or provide information for tree choice are rebutted. Measures of support are viewed as being of heuristic value within a given phylogenetic framework for describing the precision of the data based on perturbations to the data. However, no direct link is observed between the calculation of measures of support and corroboration. Direct measures of support, but not re-sampling or randomization methods, may play a more specific role in phylogenetic inference by providing the tools to search for falsifiers that could be the subject of future rounds of hypothesis testing.

The efficiency of different search strategies in estimating parsimony jackknife, bootstrap, and Bremer support

BACKGROUND

For parsimony analyses, the most common way to estimate confidence is by resampling plans (nonparametric bootstrap, jackknife), and Bremer support (Decay indices). The recent literature reveals that parameter settings that are quite commonly employed are not those that are recommended by theoretical considerations and by previous empirical studies. The optimal search strategy to be applied during resampling was previously addressed solely via standard search strategies available in PAUP*. The question of a compromise between search extensiveness and improved support accuracy for Bremer support received even less attention. A set of experiments was conducted on different datasets to find an empirical cut-off point at which increased search extensiveness does not significantly change Bremer support and jackknife or bootstrap proportions any more.

RESULTS

For the number of replicates needed for accurate estimates of support in resampling plans, a diagram is provided that helps to address the question whether apparently different support values really differ significantly. It is shown that the use of random addition cycles and parsimony ratchet iterations during bootstrapping does not translate into higher support, nor does any extension of the search extensiveness beyond the rather moderate effort of TBR (tree bisection and reconnection branch swapping) plus saving one tree per replicate. Instead, in case of very large matrices, saving more than one shortest tree per iteration and using a strict consensus tree of these yields decreased support compared to saving only one tree. This can be interpreted as a small risk of overestimating support but should be more than compensated by other factors that counteract an enhanced type I error. With regard to Bremer support, a rule of thumb can be derived stating that not much is gained relative to the surplus computational effort when searches are extended beyond 20 ratchet iterations per constrained node, at least not for datasets that fall within the size range found in the current literature.

CONCLUSION

In view of these results, calculating bootstrap or jackknife proportions with narrow confidence intervals even for very large datasets can be achieved with less expense than often thought. In particular, iterated bootstrap methods that aim at reducing statistical bias inherent to these proportions are more feasible when the individual bootstrap searches require less time.

A complete estimate of the phylogenetic relationships in Ruminantia: a dated species-level supertree of the extant ruminants

This paper presents the first complete estimate of the phylogenetic relationships among all 197 species of extant and recently extinct ruminants combining morphological, ethological and molecular information. The composite tree is derived by applying matrix representation using parsimony analysis to 164 previous partial estimates, and is remarkably well resolved, containing 159 nodes (> 80 % of the potential nodes in the completely resolved phylogeny). Bremer decay index has been used to indicate the degree of certainty associated with each clade. The ages of over 80% of the clades in the tree have been estimated from information in the literature. The supertree for Ruminantia illustrates which areas of ruminant phylogeny are still only roughly known because of taxa with controversial relationships (e.g. Odocoileini, Antilopinae) or not studied in great detail (e.g. Muntiacus). It supports the monophyly of the ruminant families and Pecora. According to this analysis Antilocapridae and Giraffidae constitute the superfamily Giraffoidea, which is the sister group of a clade clustering Bovoidea and Cervoidea. The position of several taxa whose systematic positions have remained controversial in the past (Saiga, Pelea, Aepycerus, Pantholops, Ammotragus, Pseudois) is unambiguously established. Nevertheless, the position of Neotragus and Oreotragus within the original radiation of the non-bovine bovids remains unresolved in the present analysis. It also shows that six successive rapid cladogenesis events occurred within the infraorder Pecora during the Oligocene to middle Pliocene, which coincided with periods of global climatic change. Finally, the presented supertree will be a useful framework for comparative and evolutionary biologists interested in studies involving the ruminants.

Phylogenetic signal in AFLP data sets

AFLP markers provide a potential source of phylogenetic information for molecular systematic studies. However, there are properties of restriction fragment data that limit phylogenetic interpretation of AFLPs. These are (a) possible nonindependence of fragments, (b) problems of homology assignment of fragments, (c) asymmetry in the probability of losing and gaining fragments, and (d) problems in distinguishing heterozygote from homozygote bands. In the present study, AFLP data sets of Lactuca s.l. were examined for the presence of phylogenetic signal. An indication of this signal was provided by carrying out tree length distribution skewness (g1) tests, permutation tail probability (PTP) tests, and relative apparent synapomorphy analysis (RASA). A measure of the support for internal branches in the optimal parsimony tree (MPT) was made using bootstrap, jackknife, and decay analysis. Finally, the extent of congruence in MPTs for AFLP and internal transcribed spacer (ITS)-1 data sets for the same taxa was made using the partition homogeneity test (PHT) and the Templeton test. These analytical studies suggested the presence of phylogenetic signal in the AFLP data sets, although some incongruence was found between AFLP and ITS MPTs. An extensive literature survey undertaken indicated that authors report a general congruence of AFLP and ITS tree topologies across a wide range of taxonomic groups, suggesting that the present results and conclusions have a general bearing. In these earlier studies and those for Lactuca s.l., AFLP markers have been found to be informative at somewhat lower taxonomic levels than ITS sequences. Tentative estimates are suggested for the levels of ITS sequence divergence over which AFLP profiles are likely to be phylogenetically informative.

Inferring hominoid and early hominid phylogeny using craniodental characters: the role of fossil taxa

Recent discoveries of new fossil hominid species have been accompanied by several phylogenetic hypotheses. All of these hypotheses are based on a consideration of hominid craniodental morphology. However, Collard and Wood (2000) suggested that cladograms derived from craniodental data are inconsistent with the prevailing hypothesis of ape phylogeny based on molecular data. The implication of their study is that craniodental characters are unreliable indicators of phylogeny in hominoids and fossil hominids but, notably, their analysis did not include extinct species. We report here on a cladistic analysis designed to test whether the inclusion of fossil taxa affects the ability of morphological characters to recover the molecular ape phylogeny. In the process of doing so, the study tests both Collard and Wood's (2000) hypothesis of character reliability, and the several recently proposed hypotheses of early hominid phylogeny. One hundred and ninety-eight craniodental characters were examined, including 109 traits that traditionally have been of interest in prior studies of hominoid and early hominid phylogeny, and 89 craniometric traits that represent size-corrected linear dimensions measured between standard cranial landmarks. The characters were partitioned into two data sets. One set contained all of the characters, and the other omitted the craniometric characters. Six parsimony analyses were performed; each data set was analyzed three times, once using an ingroup that consisted only of extant hominoids, a second time using an ingroup of extant hominoids and extinct early hominids, and a third time excluding Kenyanthropus platyops. Results suggest that the inclusion of fossil taxa can play a significant role in phylogenetic analysis. Analyses that examined only extant taxa produced most parsimonious cladograms that were inconsistent with the ape molecular tree. In contrast, analyses that included fossil hominids were consistent with that tree. This consistency refutes the basis for the hypothesis that craniodental characters are unreliable for reconstructing phylogenetic relationships. Regarding early hominids, the relationships of Sahelanthropus tchadensis and Ardipithecus ramidus were relatively unstable. However, there is tentative support for the hypotheses that S. tchadensis is the sister taxon of all other hominids. There is support for the hypothesis that A. anamensis is the sister taxon of all hominids except S. tchadensis and Ar. ramidus. There is no compelling support for the hypothesis that Kenyanthropus platyops shares especially close affinities with Homo rudolfensis. Rather, K. platyops is nested within the Homo + Paranthropus + Australopithecus africanus clade. If K. platyops is a valid species, these relationships suggest that Homo and Paranthropus are likely to have diverged from other hominids much earlier than previously supposed. There is no support for the hypothesis that A. garhi is either the sister taxon or direct ancestor of the genus Homo. Phylogenetic relationships indicate that Australopithecus is paraphyletic. Thus, A. anamensis and A. garhi should be allocated to new genera.