BRANCH SUPPORT AND TREE STABILITY

Abstract— Branch support is quantified as the extra length needed to lose a branch in the consensus of near-most-parsimonious trees. This approach is based solely on the original data, as opposed to the data perturbation used in the bootstrap procedure. If trees have been generated by Farris's successive approximations approach to character weighting, branch support should be examined in terms of weighted extra length needed to lose a branch. The sum of all branch support values over the tree divided by the length of the most parsimonious tree[s] provides a new index, the total support index. This index is a measure of tree stability in terms of supported resolutions, which is of prime importance in cladistic analysis.     

Generic interrelationships within the Spionidae (Annelida: Polychaeta)

The phylogenetic relationships of spionid genera are estimated from parsimony analyses of morphological characters, with Trochochaetidae, Poecilochaetidae and Uncispionidae as outgroups. A first analysis of currendy recognised genera proved inconclusive and even exclusion of six of the most polymorphic genera resulted in 13 305 equally parsimonious trees and a fully collapsed consensus tree. A second analysis using only the type species of each genus, yielded four equally parsimonious trees; reduced to two after successive weighting. The topologies of these two trees indicated division of the family into four main groups: (1) Aonidella and Xandaros; (2) Prionospio (sensu fato)-complex, Laonice, Spiophanes and Aonides; (3) a large assemblage of genera, including Polydora-{senm late), Scolelepis, Malacoceros and Spio; (4) Atherospio, Pseudatherospio and Pygospiopsis. Earlier literature classifications of the group are evaluated and compared with die new results.     

CLADISTIC AND BIOGEOGRAPHIC ANALYSES OF HAWAIIAN PIPUNCULIDAE (DIPTERA) REVISITED

Abstract — Pipunculidae (Diptera) of the Hawaiian islands belong to the endemic hawaiiensis subgroup of the subgenus Cephalops ( Semicephalops ). In total 36 species are known from the Hawaiian islands. Cladistic analysis of 31 species, using 21 morphological characters from the male terminalia, resulted in 480 equally parsimonious trees. Three rounds of successive weighting resulted in 610 equally parsimonious trees. A biogeographic analysis was carried out, using the three-area statement technique and based on the strict consensus tree produced from the tree set obtained by successive weighting. The analysis puts the sequence of the island groups in congruence with the geological history of the islands. The results of the cladistic and biogeographic analyses were compared with earlier similar analyses.     

Chloroplast DNA restriction site variation and phylogeny of the Berberidaceae

Comparative restriction site mapping of the chloroplast genome was performed to examine phylogenetic relationships among 27 species representing 16 genera of the Berberidaceae and two outgroups. Chloroplast genomes of the species included in this study showed no major structural rearrangements (i.e., they are collinear to tobacco cpDNA) except for the extension of the inverted repeat in species of Berberis and Mahonia. Excluding several regions that exhibited severe length variation, a total of 501 phylogenetically informative sites was mapped for ten restriction enzymes. The strict consensus tree of 14 equally parsimonious trees indicated that some berberidaceous genera (Berberis, Mahonia, Diphylleia) are not monophyletic. To explore phylogenetic utility of different parsimony methods phylogenetic trees were generated using Wagner, Dollo, and weighted parsimony for a reduced data set that included 18 species. One of the most significant results was the recognition of the four chromosomal groups, which were strongly supported regardless of the parsimony method used. The most notable difference among the trees produced by the three parsimony methods was the relationships among the four chromosomal groups. The cpDNA trees also strongly supported a close relationship of several generic pairs (e.g., Berberis-Mahonia, Epimedium-Vancouveria, etc.). Maximum likelihood values were computed for the four different tree topologies of the chromosomal groups, two Wagner, one Dollo, and one weighted topology. The results indicate that the weighted tree has the highest likelihood value. The lowest likelihood value was obtained for the Dollo tree, which had the highest bootstrap and decay values. Separate analyses using only the Inverted Repeat (IR) region resulted in a tree that is identical to the weighted tree. Poor resolution and/or support for the relationships among the four chromosomal lineages of the Berberidaceae indicate that they may have radiated from an ancestral stock in a relatively short evolutionary time.     

Problems in using Robertsonian rearrangements in determining monophyly: examples from the genera Tatera and Gerbillurus

Chromosomal banding data on three species of Tatera from Kenya significantly alter the previous hypothesis of relationships between and within the genera Tatera and Gerbillurus based on cladistic analyses and the rule of parsimony (Qumsiyeh, 1986b). Of the many possible hypothetical relationships, the most parsimonious tree showed three homoplasies and allowed the genus Gerbillurus to be paraphyletic. The alternative trees, depicting larger number of homoplasies but with homoplasies restricted to fusion or fission events, were compatible with the morphological data in supporting the monophyly of Gerbillurus. To choose between the hypothesis based on the most parsimonious chromosomal tree and those supported by both morphology and slightly less parsimonious chromosomal trees, we performed an electrophoretic study on this group of gerbils. The conclusions are that the genus Gerbillurus is monophyletic and represents a branch that is closely related to the T. robusta group of Taterillini. The study documents that fissions and fusions must have occurred frequently and that in some cases the same fusions were acquired in two independent lineages in numbers exceeding those that are predicted by strict parsimony. The results raise questions about the validity of systematic conclusions based solely on fusion/fission data and utilizing the parsimony criterion.     

CHLOROPLAST DNA RESTRICTION SITE VARIATION IN THE FERN GENUS PELLAEA: PHYLOGENETIC RELATIONSHIPS OF THE PELLAEA GLABELLA COMPLEX

The cheilanthoid ferns have long resisted efforts to circumscribe well-defined, phylogenetically natural generic and infrageneric groups, presumably because of homoplastic morphologies associated with their xeric habitats. This cladistic analysis of phylogenetically informative chloroplast DNA restriction site data from 14 enzymes and seven taxa in the cheilanthoid genus Pellaea provides new insights into the phylogenetic relationships of the P. glabella complex. It also assesses the congruence of results based on restriction site data at inter- and intraspecific levels in these sexually and apogamously reproducing ferns with those of earlier morphological and isozyme analyses of the same group. Wagner parsimony yielded a single most parsimonious tree of 187 steps and 11% homoplasy, based on a data matrix of 166 restriction sites of which 66 were phylogenetically informative. Phylogenetic analysis based on user-defined stepmatrix character-state weighting of site gains over losses produced an identical single most parsimonious tree. Dollo parsimony yielded two most parsimonious trees, one of which was topologically identical to the Wagner tree. Specific and infraspecific relationships in the P. glabella complex determined by the completely independent restriction site and isozyme data sets are identical. This lends confidence to the ongoing use of restriction site data in a broader study of Pellaea and other cheilanthoid taxa and to the present conclusions that P. atropurpurea is sister to the P. glabella complex, whereas P. breweri, previously considered the closest relative of this complex, is actually more distantly related to it than are the other taxa in this study.     

SYSTEMATICS AND BIOGEOGRAPHY OF THE AUSTRALIAN "GREEN ASH" EUCALYPTS (MONOCALYPTUS)

Abstract— A cladistic analysis of the "green ash" eucalypts, informal subgenus " Monocalyptus ", is presented- As a first step, ordination methods of principal coordinates analysis and multidimensional scaling delineated some terminal taxa. The cladistic analysis, applying parsimony methods to the unweighted data set, yielded 25 equally parsimonious trees, each with a consistency index of 0.57.
Farris' successive approximations approach to character weighting produced one tree with a consistency index of 0.74.
An informal classification of the group, superseries Eucalyptus , is based on that cladogram. The biogeographic history of superseries Eucalyptus is interpreted from the cladogram as having been caused by lour vicariant events in southeastern Australia, in combination with a suite of ecological features that overlie the biogeographic area-pattern.     

红豆杉科及相关类群rbcL基因PCR—RFLP分析

运用RFLP方法对红豆杉科及相关类群14种植物叶绿体rbcL基因PCR产物进行限制酶酶切分析,共获29个酶切变异位点。采用PHYLIP软件包对限制位点变异数据进行极大简约法分析得到18个步长为6的最简约树并求得一致树,结果显示：⑴红豆杉科和三尖杉科属单系群;⑵穗花杉属Amentotaxus以置于红豆杉科内为宜,不支持将穗花杉属独立成科的处理方式;⑶白豆杉应为红豆杉科内一个属Pseudotaxus;⑷三尖杉属内篦子三尖杉地位特殊,可设篦子三尖杉组;⑸不赞同将竹柏类从罗汉松属中分离出去成立新科;⑹红豆杉科、三尖杉科和罗汉松科三者间,前两者的关系更为接近。     

Morphological phylogeny of gall‐forming aphids of the tribe Eriosomatini (Aphididae: Eriosomatinae)

Aphids of the tribe Eriosomatini are typically associated with the tree genera Ulmus and Zelkova as the primary host, on which they induce several types of leaf gall. To elucidate evolutionary changes in the aphid–host associations and gall morphology, phylogenetic relationships were inferred using 36 species (28 in the ingroup) and based on 52 morphological characters. Phylogenetic analysis with equal character weighting led to hundreds of most‐parsimonious trees, and the strict consensus of these was poorly resolved. However, the successive weighting of characters yielded three most‐parsimonious trees. The strict consensus of these supported the monophyly of the Eriosomatini and the monophyly of most genera. Reconstruction of the aphid–host associations on the consensus tree indicated that the ancestral Eriosomatini were associated with Zelkova and that the association with Ulmus evolved twice independently. Ancestral reconstruction suggests that galls of the leaf‐roll type are ancestral to those of the completely and incompletely closed pouch type, and that each type of pouch galls evolved independently. It is suggested that despite early diversification of the Eriosomatini on Zelkova species, Zelkova‐associated eriosomatines had become extinct or survived as relict parthenogens on the secondary host due to the elimination of Zelkova in most regions since the late Tertiary. In contrast, two large genera in the Eriosomatini, Eriosoma and Tetraneura, are associated with the largest elm section Ulmus whose elements are distributed widely in Eurasia, including boreal regions. Therefore, the available evidence suggests that the present species diversity and distribution patterns of the eriosomatines have been largely affected by the diversification and extinction of their host plants during the late Tertiary and Quaternary.     

Phylogeny of the bee genus Halictus (Hymenoptera: halictidae) based on parsimony and likelihood analyses of nuclear EF-1alpha sequence data

We investigated higher-level phylogenetic relationships within the genus Halictus based on parsimony and maximum likelihood (ML) analysis of elongation factor-1alpha DNA sequence data. Our data set includes 41 OTUs representing 35 species of halictine bees from a diverse sample of outgroup genera and from the three widely recognized subgenera of Halictus (Halictus s.s., Seladonia, and Vestitohalictus). We analyzed 1513 total aligned nucleotide sites spanning three exons and two introns. Equal-weights parsimony analysis of the overall data set yielded 144 equally parsimonious trees. Major conclusions supported in this analysis (and in all subsequent analyses) included the following: (1) Thrincohalictus is the sister group to Halictus s.l., (2) Halictus s.l. is monophyletic, (3) Vestitohalictus renders Seladonia paraphyletic but together Seladonia + Vestitohalictus is monophyletic, (4) Michener's Groups 1 and 3 are monophyletic, and (5) Michener's Group 1 renders Group 2 paraphyletic. In order to resolve basal relationships within Halictus we applied various weighting schemes under parsimony (successive approximations character weighting and implied weights) and employed ML under 17 models of sequence evolution. Weighted parsimony yielded conflicting results but, in general, supported the hypothesis that Seladonia + Vestitohalictus is sister to Michener's Group 3 and renders Halictus s.s. paraphyletic. ML analyses using the GTR model with site-specific rates supported an alternative hypothesis: Seladonia + Vestitohalictus is sister to Halictus s.s. We mapped social behavior onto trees obtained under ML and parsimony in order to reconstruct the likely historical pattern of social evolution. Our results are unambiguous: the ancestral state for the genus Halictus is eusociality. Reversal to solitary behavior has occurred at least four times among the species included in our analysis.     

Phylogeny of the Bee Genus Halictus (Hymenoptera: Halictidae) Based on Parsimony and Likelihood Analyses of Nuclear EF-1α Sequence Data

We investigated higher-level phylogenetic relationships within the genus Halictus based on parsimony and maximum likelihood (ML) analysis of elongation factor-1α DNA sequence data. Our data set includes 41 OTUs representing 35 species of halictine bees from a diverse sample of outgroup genera and from the three widely recognized subgenera of Halictus (Halictus s.s., Seladonia, and Vestitohalictus). We analyzed 1513 total aligned nucleotide sites spanning three exons and two introns. Equal-weights parsimony analysis of the overall data set yielded 144 equally parsimonious trees. Major conclusions supported in this analysis (and in all subsequent analyses) included the following: (1) Thrincohalictus is the sister group to Halictus s.l., (2) Halictus s.l. is monophyletic, (3) Vestitohalictus renders Seladonia paraphyletic but together Seladonia + Vestitohalictus is monophyletic, (4) Michener's Groups 1 and 3 are monophyletic, and (5) Michener's Group 1 renders Group 2 paraphyletic. In order to resolve basal relationships within Halictus we applied various weighting schemes under parsimony (successive approximations character weighting and implied weights) and employed ML under 17 models of sequence evolution. Weighted parsimony yielded conflicting results but, in general, supported the hypothesis that Seladonia + Vestitohalictus is sister to Michener's Group 3 and renders Halictus s.s. paraphyletic. ML analyses using the GTR model with site-specific rates supported an alternative hypothesis: Seladonia + Vestitohalictus is sister to Halictus s.s. We mapped social behavior onto trees obtained under ML and parsimony in order to reconstruct the likely historical pattern of social evolution. Our results are unambiguous: the ancestral state for the genus Halictus is eusociality. Reversal to solitary behavior has occurred at least four times among the species included in our analysis.     

Testing the monophyly of the New Zealand and Australian endemic family Conoesucidae Ross based on combined molecular and morphological data (Insecta: Trichoptera: Sericostomatoidea)

Conoesucidae (Trichoptera, Insecta) are restricted to SE Australia, Tasmania and New Zealand. The family includes 42 described species in 12 genera, and each genus is endemic to either New Zealand or Australia. Although monophyly has been previously assumed, no morphological characters have been proposed to represent synapomorphies for the group. We collected molecular data from two mitochondrial genes (16S and cytochrome oxidase I), one nuclear gene (elongation factor 1-α) (2237–2277 bp in total), and 12 morphological characters to produce the first phylogeny of the family. We combined the molecular and morphological characters and performed both a maximum parsimony analysis and a Bayesian analysis to test the monophyly of the family, and to hypothesize the phylogeny among its genera. The parsimony analysis revealed a single most parsimonious tree with Conoesucidae being a monophyletic taxon and sistergroup to the Calocidae. The Bayesian inference produced a distribution of trees, the consensus of which is supported with posterior probabilities of 100% for 15 out of 22 possible ingroup clades including the most basal branch of the family, indicating strong support for a monophyletic Conoesucidae. The most parsimonious tree and the tree from the Bayesian analysis were identical except that the ingroup genus Pycnocentria changed position by jumping to a neighbouring clade. Based on the assumption that the ancestral conoesucid species was present on both New Zealand and Australia, a biogeographical analysis using the dispersal-vicariance criteria demonstrated that one or two (depending on which of the two phylogenetic reconstructions were applied) sympatric speciation events took place on New Zealand prior to a single, late dispersal from New Zealand to Australia.     

Plethodontid salamander mitochondrial genomics: A parsimony evaluation of character conflict and implications for historical biogeography

A new parsimony analysis of 27 complete mitochondrial genomic sequences is conducted to investigate the phylogenetic relationships of plethodontid salamanders. This analysis focuses on the amount of character conflict between phylogenetic trees recovered from newly conducted parsimony searches and the Bayesian and maximum likelihood topology reported by Mueller et al. (2004 ; PNAS, 101, 13820–13825). Strong support for Hemidactylium as the sister taxon to all other plethodontids is recovered from parsimony analyses. Plotting area relationships on the most parsimonious phylogenetic tree suggests that eastern North America is the origin of the family Plethodontidae supporting the “Out of Appalachia” hypothesis. A new taxonomy that recognizes clades recovered from phylogenetic analyses is proposed. © The Willi Hennig Society 2005.     

红豆杉科及相关类群rbc L基因PCR-RFLP分析

运用RFLP方法对红豆杉科及相关类群14种植物叶绿体rbcL基因PCR产物进行限制酶酶切分析,共获29个酶切变异位点。采用PHYLIP 软件包对限制位点变异数据进行极大简约法分析得到18个步长为6的最简约树并求得一致树,结果显示：（1）红豆杉科和三尖杉科属单系群;(2) 穗花杉属Amentotaxus以置于红豆杉科内为宜,不支持将穗花杉属独立成科的处理方式;(3)白豆杉应为红豆杉科内一个属Pseudotaxus;(4) 三尖杉属内篦子三尖杉地位特殊,可设篦子三尖杉组;(5) 不赞同将竹柏类从罗汉松属中分离出去成立新科;(6) 红豆杉科、三尖杉科和罗汉松科三者间,前两者的关系更为接近。     

Parsimony-based test for identifying changes in evolutionary trends for quantitative characters: implications for the origin of the amniotic egg

The origin of the amniotic egg was a major event in vertebrate evolution and is thought to have contributed to the spectacular evolutionary radiation of amniotes. We test one of the most popular scenarios proposed by Carroll in 1970 to explain the origin of the amniotic egg using a novel method based on an asymmetric version of linear parsimony (aka Wagner parsimony) for identifying the most parsimonious split of a tree into two parts between which the evolution of the character is allowed to differ. The new method evaluates the cost of splitting a phylogenetic tree at a given node as the integral, over all pairs of asymmetry parameters, of the most parsimonious costs that can be achieved by using the first parameter on the subtree pending from this node and the second parameter elsewhere. By testing all the nodes, we then obtain the most parsimonious split of a tree with regard to the character values at its tips. Among the nine trees and two characters tested, our method yields a total of 517 parsimonious trend changes in Permo-Carboniferous stegocephalians, a single one of which occurs in a part of the tree (among stem-amniotes) where Carroll's scenario predicts that there should have been distinct changes in body size evolutionary trends. This refutes the scenario because the amniote stem does not appear to have elevated rates of evolutionary trend shifts. Our nodal body size estimates offer less discriminating power, but they likewise fail to find strong support for Carroll's scenario.     

Phylogeny of nereidids (Polychaeta, Nereididae) with paragnaths

A phylogenetic analysis was conducted of the Nereidinae — those members of the Nereididae (Polychaeta) with pharyngeal paragnaths. We had two objectives: to test the monophyly of currently accepted genera, subgenera and informal subgeneric groupings within the Nereidinae, and, if warranted, to propose a more natural classification of the Nereidinae. Parsimony analyses were undertaken, including 52 terminal taxa from all genera and informal groupings from the large heterogeneous genera Nereis , Ceratonereis , Neanthes and Perinereis . Analyses of a character set of 52 informative characters yielded more than 10 000 equally parsimonious trees with a length of 176 steps (consistency index [CI] = 0.34, retention index [RI] = 0.66). Reweighting three times resulted in 445 most parsimonious trees with length 54.62 (CI = 0.59, RI = 0.79). Many characters widely used in nereidid systematics were found to exhibit high levels of homoplasy. The most parsimonious trees could not be rooted such that the selected ingroup, 'Nereididae with paragnaths', was monophyletic, causing us to reject the monophyly of the Nereidinae as currently defined. The following genera were well supported by the parsimony analyses and are newly diagnosed: Alitta , Ceratonereis , Pseudonereis , Simplisetia , Solomononereis and Unanereis . Alitta succinea , Pseudonereis cortezi , Pseudonereis noodti and Pseudonereis pseudonoodti are proposed as new combinations. The parsimony analysis supported the monophyly of neither Composetia , Neanthes , Nereis and Perinereis nor of any new groupings of remaining species presently placed in those genera. It is these poorly supported genera that comprise most species of Nereididae.     

An investigation into the application of the Wagner parsimony method in synecology

Synecological analyses are usually based on typological, phenetic and cladistic methods. The disadvantages of these techniques are shown. The application of the Wagner parsimony method to synecology is considered. All the methods need some prerequisites, viz. definitions of localities and characters (the most simple one being the presence/absence of taxa); the choice of taxonomic level of taxa; their autochthony. The application of Wagner parsimony needs a new terminology. The congruence of any environmental condition, including freshwater monitoring indices, can be tested on parsimonious trees. The Wagner parsimony method not only provides various indices (tree length, CI, HI, RC, RI) which allow the comparison of trees but also minimal trees which are direct tools in synecology.     

Maximum Parsimony on Phylogenetic networks

Background

Phylogenetic networks are generalizations of phylogenetic trees, that are used to model evolutionary events in various contexts. Several different methods and criteria have been introduced for reconstructing phylogenetic trees. Maximum Parsimony is a character-based approach that infers a phylogenetic tree by minimizing the total number of evolutionary steps required to explain a given set of data assigned on the leaves. Exact solutions for optimizing parsimony scores on phylogenetic trees have been introduced in the past.

Results

In this paper, we define the parsimony score on networks as the sum of the substitution costs along all the edges of the network; and show that certain well-known algorithms that calculate the optimum parsimony score on trees, such as Sankoff and Fitch algorithms extend naturally for networks, barring conflicting assignments at the reticulate vertices. We provide heuristics for finding the optimum parsimony scores on networks. Our algorithms can be applied for any cost matrix that may contain unequal substitution costs of transforming between different characters along different edges of the network. We analyzed this for experimental data on 10 leaves or fewer with at most 2 reticulations and found that for almost all networks, the bounds returned by the heuristics matched with the exhaustively determined optimum parsimony scores.

Conclusion

The parsimony score we define here does not directly reflect the cost of the best tree in the network that displays the evolution of the character. However, when searching for the most parsimonious network that describes a collection of characters, it becomes necessary to add additional cost considerations to prefer simpler structures, such as trees over networks. The parsimony score on a network that we describe here takes into account the substitution costs along the additional edges incident on each reticulate vertex, in addition to the substitution costs along the other edges which are common to all the branching patterns introduced by the reticulate vertices. Thus the score contains an in-built cost for the number of reticulate vertices in the network, and would provide a criterion that is comparable among all networks. Although the problem of finding the parsimony score on the network is believed to be computationally hard to solve, heuristics such as the ones described here would be beneficial in our efforts to find a most parsimonious network.