Implied weighting and its utility in palaeontological datasets: a study using modelled phylogenetic matrices

Implied weighting, a method for phylogenetic inference that actively seeks to downweight supposed homoplasy, has in recent years begun to be widely utilized in palaeontological datasets. Given the method's purported ability at handling widespread homoplasy/convergence, we investigate the effects of implied weighting on modelled phylogenetic data. We generated 100 character matrices consisting of 55 characters each using a Markov Chain morphology model of evolution based on a known phylogenetic tree. Rates of character evolution in these datasets were variable and generated by pulling from a gamma distribution for each character in the matrix. These matrices were then analysed under equal weighting and four settings of implied weights (k = 1, 3, 5, and 10). Our results show that implied weighting is inconsistent in its ability to retrieve a known phylogenetic tree. Equally weighted analyses are found to generally be more conservative, retrieving higher frequency of polytomies but being less likely to generate erroneous topologies. Implied weighting is found to generally resolve polytomies while also propagating errors, resulting in an increase in both correctly and incorrectly resolved nodes with a tendency towards higher rates of error compared to equal weighting. Our results suggest that equal weights may be a preferable method for parsimony analysis.     

Uninode coding vs gene tree parsimony for phylogenetic reconstruction using duplicate genes

Two different methods of using paralogous genes for phylogenetic inference have been proposed: reconciled trees (or gene tree parsimony) and uninode coding. Gene tree parsimony suffers from 10 serious problems, including differential weighting of nucleotide and gap characters, undersampling which can be misinterpreted as synapomorphy, all of the characters not being allowed to interact, and conflict between gene trees being given equal weight, regardless of branch support. These problems are largely avoided by using uninode coding. The uninode coding method is elaborated to address multiple gene duplications within a single gene tree family and handle problems caused by lack of gene tree resolution. An example of vertebrate phylogeny inferred from nine genes is reanalyzed using uninode coding. We suggest that uninode coding be used instead of gene tree parsimony for phylogenetic inference from paralogous genes.     

THE PHYLOGENY OFLECANACTIS(OPEGRAPHACEAE)

The genusLecanactis, with 24 species, has been phylogenetically analysed using cladistic parsimony methods and support tests. Morphological, anatomical and chemical data were used, comprising 38 characters. Twelve equally most parsimonious trees were obtained. The successive approximations character weighting method gave one most parsimonious tree. The ingroup,Lecanactis, is supported as monophyletic. Although parsimony jackknifing and Bremer support indicate that the trees are poorly supported, some groups are wholly or partly distinguished in both the strict consensus tree, the successive weighting tree and the Jac tree.     

Probabilistic methods outperform parsimony in the phylogenetic analysis of data simulated without a probabilistic model

To understand patterns and processes of the diversification of life, we require an accurate understanding of taxon interrelationships. Recent studies have suggested that analyses of morphological character data using the Bayesian and maximum likelihood Mk model provide phylogenies of higher accuracy compared to parsimony methods. This has proved controversial, particularly studies simulating morphology‐data under Markov models that assume shared branch lengths for characters, as it is claimed this leads to bias favouring the Bayesian or maximum likelihood Mk model over parsimony models which do not explicitly make this assumption. We avoid these potential issues by employing a simulation protocol in which character states are randomly assigned to tips, but datasets are constrained to an empirically realistic distribution of homoplasy as measured by the consistency index. Datasets were analysed with equal weights and implied weights parsimony, and the maximum likelihood and Bayesian Mk model. We find that consistent (low homoplasy) datasets render method choice largely irrelevant, as all methods perform well with high consistency (low homoplasy) datasets, but the largest discrepancies in accuracy occur with low consistency datasets (high homoplasy). In such cases, the Bayesian Mk model is significantly more accurate than alternative models and implied weights parsimony never significantly outperforms the Bayesian Mk model. When poorly supported branches are collapsed, the Bayesian Mk model recovers trees with higher resolution compared to other methods. As it is not possible to assess homoplasy independently of a tree estimate, the Bayesian Mk model emerges as the most reliable approach for categorical morphological analyses.     

Conflict and congruence in a combined DNA–morphology analysis of megachiropteran bat relationships (Mammalia: Chiroptera: Pteropodidae)

The phylogeny of megabats (Mammalia: Chiroptera: Megachiroptera) has been addressed only on molecular grounds, as little effort has previously been made to describe the impressive morphological variation of the group in terms of phylogenetically informative characters. Here we provide a morphological matrix of 236 characters from the integument, dentition, cranial and post‐cranial skeleton, digestive apparatus and urogenital system. This data set covers most characters discussed previously in more restricted taxonomic contexts, as well a large number of new characters. Our aim was to generate a phylogenetic hypothesis for megabats based on a combined analysis of morphological characters and available gene sequence data from four mitochondrial and one nuclear loci. We used direct optimization under conventional equal costs, as well as under a cost ratio that maximizes homology when inapplicables (gaps) are present. Our results contradict the allegedly high level of conflict between the molecular and morphological partitions. We found that, although morphology alone recovered trees different and to some extent incompatible with those from previous molecular analyses, the combination of the two sources of evidence easily accommodated the morphological and molecular signals, yielding a resolved and relatively well‐supported phylogeny of Megachiroptera that is in reasonable agreement with the current morphology‐based taxonomy of the group. Overall congruence favored the maximization of homology by a narrow margin. In addition, partial analyses showed that implied weighting of morphology performed slightly better than equal weighting with respect to the combined analyses. © The Willi Hennig Society 2005.     

Phylogeny of the parasitic wasp subfamily Charipinae (Hymenoptera, Cynipoidea, Figitidae)

The Charipinae are a major group of hyperparasitoids of Hemiptera. Here, we present the first cladistic analysis of this subfamily's internal relationships, based on 96 morphological characters of adults. The data matrix was analysed using uniformly weighted parsimony. The effects of using alternative weighting schemes were explored by performing additional searches employing implied weights criteria. One of the caveats of implied weights analysis is that it lacks an objective criterion for selecting the value of the concavity function. In the present study, differential weighting was used to explore the sensitivity of our results to the alternative assumptions made in the analysis and to select one of the most parsimonious trees under equal weights, which we regard as being the hypothesis that minimizes the amount of ad hoc assumptions. The validity of the two existing tribes and the monophyly of all the genera of Charipinae were tested, in particular the cosmopolitan and highly species-rich Alloxysta and Phaenoglyphis , which appear repeatedly in ecological and biochemical studies of host–parasitoid associations. The evolution of several major characters and the relationships between genera are discussed. On the basis of the phylogenetic results, we discuss a number of taxonomic issues. A new classification of the subfamily is proposed in which no tribes are maintained, Carvercharips is synonymyzed with Alloxysta , and the creation of a new genus from Nepal is justified. Our analysis points to the need for a world revision of the basal genus Phaenoglyphis , which is shown as paraphyletic.     

A morphological phylogeny of Agrotis Ochsenheimer (Lepidoptera,Noctuidae), with emphasis on the South American species

The monophyly of Agrotis Ochsenheimer in addition to the relationships among the South American species of the genus and the relationship between Agrotis and Feltia Walker are examined. A cladistic analysis was performed based on 45 terminals, with 28 ingroup and 17 outgroup taxa, and 126 characters (two continuous and 124 discrete) from adult morphology, including male and female genitalia. Parsimony analyses were performed under equal and implied weighting. Results support Agrotis as a monophyletic group, sister to the South American species of Feltia (Feltia tent.), and the latter being more closely related to Agrotis than to the ‘true’ Feltia. Species of Agrotis that branched off first (including the type species) have characters shared with both Feltia and Agrotis. South American species of Agrotis (with some proviso) form a clade that branched off later, with some weakly supported species groups that had been proposed in previous works. Biogeographical aspects of the group are discussed after optimizing species distributions on the tree. Trees obtained after analyses using implied weights show similar relationships to those under equal weights, particularly regarding Feltia tent., showing Feltia tent. plus Agrotis forming a monophyletic group, sister to the ‘true’ Feltia.     

Phylogeny of the Hymenoptera: a reanalysis of the Ronquist et al. (1999) reanalysis, emphasizing wing venation and apocritan relationships

Ronquist et al . (1999) recoded Rasnitsyn's (1988) analysis to effect an explicit numerical cladistic analysis of his data. Here we examine their analysis and reveal that much of the resolution obtained for apocritan relationships is dependant on reductional wing characters. The wing characters in their matrix are replaced with a revised set of wing characters and reanalysed using strict parsimony. For apocritan taxa the resulting strict consensus tree is considerably less resolved, but perhaps preferable as a more conservative starting point in the continuing investigation of higher level hymenopteran relationships.     

Phylogeny of the genera of the parasitic wasps subfamily Doryctinae (Hymenoptera: Braconidae) based on morphological evidence

The evolutionary relationships among most (143 genera) of the currently recognized genera of the braconid wasp subfamily Doryctinae were investigated using maximum parsimony analysis, employing 100 characters from external morphology and four additional, less well‐studied character systems (male genitalia, ovipositor structure, venom apparatus and larval cephalic structure). We investigated the ‘performance’ of characters from external morphology and the other character systems and the effects of abundant missing entries by comparing the data decisiveness, retention and consistency indices of four different character partitions. The results indicate that the performances of the different partitions are not related to the proportions of missing entries, but instead are negatively correlated to their proportion of informative characters, suggesting that the morphological information in this group is subject to high levels of homoplasy. The external morphological partition is significantly incongruent with respect to a data set comprising the other character systems based on the ILD test. Analyses supported neither the monophyly of the large tribes Doryctini and Hecabolini, nor the monophyly of the Spathiini and Heterospilini. Relationships obtained from successive approximation weighting analysis for the complete data differ considerably from the currently accepted tribal and subtribal classifications. The only exceptions were the Ypsistocerini and the Ecphylini, whose recognized members were recovered in single clades. A close relationship between the Binaerini and Holcobraconini, and also Monarea, is consistently supported by venom apparatus and ovipositor structure characters but is not indicated by external morphological data. Low bootstrap values obtained for most of the recovered clades in all analyses do not allow us to propose a meaningful reclassification for the group at this time. A complete list of the recognized genera and their synonymies is given. © 2004 The Linnean Society of London, Zoological Journal of the Linnean Society, 2004, 142 , 369–404.     

Extended implied weighting

Several extensions to implied weighting, recently implemented in TNT, allow a better treatment of data sets combining morphological and molecular data sets, as well as those comprising large numbers of missing entries (e.g. palaeontological matrices, or combined matrices with some genes sequenced for few taxa). As there have been recent suggestions that molecular matrices may be better analysed using equal weights (rather than implied weighting), a simple way to apply implied weighting to only some characters (e.g. morphology), leaving other characters with a constant weight (e.g. molecules), is proposed. The new methods also allow weighting entire partitions according to their average homoplasy, giving each of the characters in the partition the same weight (this can be used for dynamically weighting, e.g. entire genes, or first, second, and third positions collectively). Such an approach is easily implemented in schemes like successive weighting, but in the case of implied weighting poses some particular problems. The approach has the peculiar implication that the inclusion of uninformative characters influences the results (by influencing the implied weights for the partitions). Last, the concern that characters with many missing entries may receive artificially inflated weights (because they necessarily display less homoplasy) can be solved by allowing the use of different weighting functions for different characters, in such a way that the cost of additional transformations decreases more rapidly for characters with more missing entries (thus effectively assuming that the unobserved entries are likely to also display some unobserved homoplasy). The conceptual and practical aspects of all these problems, as well as details of the implementation in TNT, are discussed.     

Phylogeny of snakes (Serpentes): Combining morphological and molecular data in likelihood,Bayesian and parsimony analyses

Abstract

The phylogeny of living and fossil snakes is assessed using likelihood and parsimony approaches and a dataset combining 263 morphological characters with mitochondrial (2693 bp) and nuclear (1092 bp) gene sequences. The ‘no common mechanism’ (NCMr) and ‘Markovian’ (Mkv) models were employed for the morphological partition in likelihood analyses; likelihood scores in the NCMr model were more closely correlated with parsimony tree lengths. Both models accorded relatively less weight to the molecular data than did parsimony, with the effect being milder in the NCMr model. Partitioned branch and likelihood support values indicate that the mtDNA and nuclear gene partitions agree more closely with each other than with morphology. Despite differences between data partitions in phylogenetic signal, analytic models, and relative weighting, the parsimony and likelihood analyses all retrieved the following widely accepted groups: scolecophidians, alethinophidians, cylindrophiines, macrostomatans (sensu lato) and caenophidians. Anilius alone emerged as the most basal alethinophidian; the combined analyses resulted in a novel and stable position of uropeltines and cylindrophiines as the second‐most basal clade of alethinophidians. The limbed marine pachyophiids, along with Dinilysia and Wonambi, were always basal to all living snakes. Other results stable in all combined analyses include: Xenopeltis and Loxocemus were sister taxa (fide morphology) but clustered with pythonines (fide molecules), and Ungaliophis clustered with a boine‐erycine clade (fide molecules). Tropidophis remains enigmatic; it emerges as a basal alethinophidian in the parsimony analyses (fide molecules) but a derived form in the likelihood analyses (fide morphology), largely due to the different relative weighting accorded to data partitions.     

The phylogenetic relationships of Morgan's Sphinx, Xanthopan morganii (Walker), the tribe Acherontiini, and allied long-tongued hawkmoths (Lepidoptera: Sphingidae, Sphinginae)

A cladistic analysis is presented of the hawkmoths of the tribe Acherontiini, Morgan's Sphinx (Xanthopan morganii (Walker)), and related genera. The study aims to test the monophyly of tribe Acherontiini; the hypothesis that all taxa with extremely long probosces (some Acherontiini, Meganoton rubescens, Neococytius, Xanthopan) form a monophyletic group, or at least fall within a single reasonably compact clade; and, within this group, to determine whether Xanthopan is more closely related to Acherontiini or to Cocytius and Neococytius. The data set comprises 109 characters derived from adult and immature stage morphology, biology and behaviour. These data were analysed using equal weighting, successive approximations character weighting (SACW) and implied weighting. All weighting schemes agreed on the monophyly of Acherontiini and of a group of genera comprising Amphimoea, Cocytius and Neococytius (the Cocytius group). Several other generic and suprageneric clades were also consistently recovered. However, those hawkmoths with extremely long probosces were never recovered as a monophyletic group. The relationships of Xanthopan were also ambiguous. Equal weighting and SACW placed Xanthopan + Meganoton rubescens (Butler) as sister to the Cocytius group, while implied weighting placed Xanthopan as sister to Acherontiini. This latter relationship is based primarily on shared possession of a pilifer/palp hearing organ. Further analyses suggested the two components of this organ were not biologically independent. Downweighting this feature accordingly resulted in all weighting schemes converging on the topology found by equal weighting. Exclusion of the incomplete subset of immature stage data had no effect under implied weighting but equal weighting and SACW now recovered a Neotropical clade comprising Manduca and the Cocytius group, while Xanthopan was placed with M. rubescens and Panogena. Downweighting the pilifer/palp hearing organ under implied weighting again caused convergence with the equal weighting/SACW results. Thus, the relationships of Xanthopan remain equivocal and further data, particularly from the immature stages, will be required to elucidate its phylogenetic position further. © 2002 The Linnean Society of London. Zoological Journal of the Linnean Society, 2002, 135 , 471–527.     

Combined phylogenetic analysis of the subclass Marchantiidae (Marchantiophyta): towards a robustly diagnosed classification

The most extensive combined phylogenetic analyses of the subclass Marchantiidae yet undertaken was conducted on the basis of morphological and molecular data. The morphological data comprised 126 characters and 56 species. Taxonomic sampling included 35 ingroup species with all genera and orders of Marchantiidae sampled, and 21 outgroup species with two genera of Blasiidae (Marchantiopsida), 15 species of Jungermanniopsida (the three subclasses represented) and the three genera of Haplomitriopsida. Takakia ceratophylla (Bryophyta) was employed to root the trees. Character sampling involved 92 gametophytic and 34 sporophytic traits, supplemented with ten continuous characters. Molecular data included 11 molecular markers: one nuclear ribosomal (26S), three mitochondrial genes (nad1, nad5, rps3) and seven chloroplast regions (atpB, psbT‐psbH, rbcL, ITS, rpoC1, rps4, psbA). Searches were performed under extended implied weighting, weighting the character blocks against the average homoplasy. Clade stability was assessed across three additional weighting schemes (implied weighting corrected for missing entries, standard implied weighting and equal weighting) in three datasets (molecular, morphological and combined). The contribution from different biological phases regarding node recovery and diagnosis was evaluated. Our results agree with many of the previous studies but cast doubt on some relationships, mainly at the family and interfamily level. The combined analyses underlined the fact that, by combining data, taxonomic enhancements could be achieved regarding taxon delimitation and quality of diagnosis. Support values for many clades of previous molecular studies were improved by the addition of morphological data. The long‐held assumption that morphology may render spurious or low‐quality results in this taxonomic group is challenged. The morphological trends previously proposed are re‐evaluated in light of the new phylogenetic scheme.     

Phylogeny of the death's head hawkmoths, Acherontia[Laspeyres], and related genera (Lepidoptera: Sphingidae: Sphinginae: Acherontiini)

Abstract. Adult death's head hawkmoths (Acherontia species) have a unique feeding biology as cleptoparasites of honeybees, stealing honey from the combs, rather than imbibing nectar from flowers. The moths have a range of features, both morphological and behavioural, that enable them to successfully enter, feed and escape from the colonies. These adaptations may vary among the three Acherontia species and allow them each to target different species of honeybee. A cladistic analysis is presented of the hawkmoths of tribe Acherontiini. The study aims to resolve the relationships of the genera and species of Acherontiini, with a particular focus on the three species of Acherontia. The dataset comprises sixty‐five characters derived from adult, larval and pupal morphology, and larva host‐plant biology. These data are analysed using equal weighting and implied weighting. Acherontiini and each constituent genus are recovered as monophyletic. However, within Coelonia, there is ambiguity in that the sister‐species relationships C. brevis+C. fulvinotata and C. fulvinotata + C. solani are equally parsimonious under both weighting schemes. Furthermore, under equal weighting Agrius is placed equally parsimoniously as the sister group of either Acherontia + Coelonia or Callosphingia. Under implied weighting, however, only the latter relationship is most parsimonious (fit). Within Acherontia, A. atropos and A. styx are always recovered as sister species to the exclusion of A. lachesis. The results of the phylogenetic analysis provide an objective basis for future studies of the unique cleptoparasitic association of these moths.     

Noise

The proliferation of DNA sequence data has generated a concern about the effects of "noise" on phylogeny reconstruction. This concern has led to various recommendations for weighting schemes and for separating data types prior to analysis. A new technique is explored to examine directly how noise influences the stability of parsimony reconstruction. By appending purely random characters onto a matrix of pure signal, or by replacing characters in a matrix of signal by random states, one can measure the degree to which a matrix is robust against noise. Reconstructions were sensitive to tree topology and clade size when noise was added, but were less so when character states were replaced with noise. When a signal matrix is complemented with a noise matrix of equal size, parsimony will trace the original signal about half the time when there is only one synapomorphy per node, and about 90% of the time when there are three synapomorphies per node. Similar results obtain when 20% of a matrix is replaced by noise. Successive weighting does not improve performance. Adding noise to only some taxa is more damaging, but replacing characters in only some taxa is less so. The bootstrap and g1 (tree skewness) statistics are shown to be uninterpretable measures of noise or departures from randomness. Empirical data sets illustrate that commonly recommended schemes of differential weighting (e.g. downweighting third positions) are not well supported from the point of view of reducing the influence of noise nor are more noisy data sets likely to degrade signal found in less noisy data sets.     

Morphology,molecules and the phylogeny of Characidae (Teleostei,Characiformes)

This is the most comprehensive phylogenetic analysis of the Characidae to date and the first large-scale hypothesis of the family, combining myriad morphological data with molecular information. A total of 520 morphological characters were analysed herein, of which 98 are newly defined. Among the analysed taxa, 259 species were coded by examining specimens, three fossil species were coded from the literature, one species was coded almost completely from published figures, 122 were partially coded from the literature, and 88 were analysed exclusively from molecular data. The total number of species in the analysed dataset is 473. Analyses were made by parsimony under equal and extended implied weighting with a broad range of parameters. The final hypothesis was selected using a stability criterion that chooses among the most parsimonious trees of all searches. It was found by weighting molecular characters with the average homoplasy of entire partitions (markers). The resulting hypothesis is congruent with previous molecular-based phylogenies of the family. The Characidae are monophyletic, with four main clades: the Spintherobolinae new subfamily; an expanded Stethaprioninae including the Grundulini, Gymnocharacini, Rhoadsiini and Stethaprionini; the Stevardiinae; and a clade composed of the Aphyocharacinae, Characinae, Cheirodontinae, Exodontinae and Tetragonopterinae. Also, a stem Characidae was found, as formed by the Eocene–Oligocene genera †Bryconetes and †Paleotetra as successive sister groups of extant members of the family. A subfamilial classification is proposed, but deep changes in the systematics that are beyond the scope of this study are still needed to classify the Characidae into monophyletic genera.     

Phylogenetic inference using discrete characters: performance of ordered and unordered parsimony and of three‐item statements

The cladistic literature does not always specify the kind of multistate character treatment that is applied for an analysis. Characters can be treated either as unordered transformation series or as rooted [three‐item analysis (3ia)] or unrooted state trees (ordered characters). We aimed to measure the impact of these character treatments on phylogenetic inference. Discrete characters can be represented either as rows or columns in matrices (e.g. for parsimony) or as hierarchies for 3ia. In the present study, we use simulated and empirical examples to assess the relative merits of each method considering both the character treatment and representation. We measure two parameters (resolving power and artefactual resolution) using a new tree comparison metric, ITRI (inter‐tree retention index). Our results suggest that the hierarchical character representation not only results (with our simulation settings) in the greatest resolving power, but also in the highest artefactual resolution. Our empirical examples provide equivocal results. Parsimony unordered states yield less resolving power and more artefactual resolutions than parsimony ordered states, both with our simulated and empirical data. Relationships between three operational taxonomic units (OTUs), irrespective of their relationships with other OTUs, are called three‐item statements (3is). We compare the intersection tree (which reconstructs a single tree from all of the common 3is of source trees) with the traditional strict consensus and show that the intersection tree retains more of the information contained in the source trees. © 2013 The Linnean Society of London, Biological Journal of the Linnean Society, 2013, 110 , 914–930.     

The phylogeny of elateroidea (Coleoptera), or which tree is best today?

Abstract — The relationships of the clicking Elateroidea beetles were studied with the help of parsimony analysis using Hennig86. The character matrix included 70 characters and 27 taxa. The results demonstrate the monophyly of the group Throscidae sensu Crowson, contrary to views presented in other papers. Methods for solving this problem were sought. When several minimum length solutions were obtained, successive weighting and a search for a strict consensus tree identical with one of the original trees appeared to be acceptable ways for trying to identify the preferred solution. When conflicting trees from separate data sets were compared, a combined global analysis turned out to be impossible to perform because the data sets used different terminal taxa. In this case, the incongruence and total support tests provided by Farris' programs RNA and KON proved indispensable. The conflict found between the results obtained here and those presented by other workers using a large suite of larval characters were shown to be caused by an incongruent data matrix used in the latter study—the larval data set resulted in a polyphyletic ingroup and suggests relationships quite different from adult data alone. Directed large scale homoplasy due to repeated re-invasion of two major habitats by separate clades may be the factor causing difficulties in coding the larval characters.     

First phylogenetic analysis of the family Neriidae (Diptera), with a study on the issue of scaling continuous characters

Neriidae are a small family of acalyptratae flies, mostly distributed in the tropics. Very little is known about their biology, and the evolutionary relationships among species have never been evaluated. We perform the first comprehensive phylogenetic analysis of the family, including 48 species from all biogeographic regions inhabited, as well as five species of Micropezidae and one Cypselosomatidae as outgroups. We build a morphological data matrix of 194 characters, including 72 continuous characters. We first explore ways to deal with the issue of scaling continuous characters, including rescaling ranges to unity and using implied weighting. We find that both strategies result in very different phylogenetic hypotheses, and that implied weighting reduces the problem of scaling, but only partially. Furthermore, using implied weighting after rescaling characters improves the congruence between partitions and results in higher values of group support. With respect to the Neriidae, we confirm the monophyly of the family and of most its genera, although we do not obtain any of the currently accepted suprageneric groups. We propose to restrict the Eoneria and Nerius groups exclusively to the Neotropical fauna, and synonymize Glyphidops subgenus Oncopsia Enderlein with Glyphidops subgenus Glyphidops Enderlein, eliminating the subgeneric divisions. This revised phylogeny presents a striking biogeographic consistency, and shows that previous main divisions of the family were based on events of convergence.     

单子叶植物科级分类阶元形态性状分支分析

基于93个形态形状,采用13个被子植物基部类群做为外类群,对49个单子叶植物科级分类阶元进行了分支系统学分析。经过简约性分析,得到了1684棵同等最大简约分支树。严格一致树的分支结构图表明：1）古草本类植物和单子叶植物是姐妹群关系;2）具有网状脉的类群,薯蓣科,菝葜科,百部科是单子叶植物的最基部类群。由于性状状态间存在着较多的平行和逆转进化,这在一定程度上影响了系统发育重建的准确性;所选择的性状状态之间的演化很可能是平行的、多次的或者是特化的状态,因此这样复杂的演化关系的探索关键在于找到一些能确切反映其系统演化关系的形态性状。目前很难通过简约化的形态分支分析来解开整个单子叶植物的起源和演化之谜。为了避开对系统学分析造成干扰的误导性状,形态数据结合DNA序列分析很可能是必需的。