The problem of rooting rapid radiations

There are many examples of groups (such as birds, bees, mammals, multicellular animals, and flowering plants) that have undergone a rapid radiation. In such cases, where there is a combination of short internal and long external branches, correctly estimating and rooting phylogenetic trees is known to be a difficult problem. In this simulation study, we tested the performances of different phylogenetic methods at estimating a tree that models a rapid radiation. We found that maximum likelihood, corrected and uncorrected neighbor-joining, and corrected and uncorrected parsimony, all suffer from biases toward specific tree topologies. In addition, we found that using a single-taxon outgroup to root a tree frequently disrupts an otherwise correct ingroup phylogeny. Moreover, for uncorrected parsimony, we found cases where several individual trees (in which the outgroup was placed incorrectly) were selected more frequently than the correct tree. Even for parameter settings where the correct tree was selected most frequently when using extremely long sequences, for sequences of up to 60,000 nucleotides the incorrectly rooted trees were each selected more frequently than the correct tree. For all the cases tested here, tree estimation using a two taxon outgroup was more accurate than when using a single-taxon outgroup. However, the ingroup was most accurately recovered when no outgroup was used.     

An empirical test of the midpoint rooting method

The outgroup method is widely used to root phylogenetic trees. An accurate root indication, however, strongly depends on the availability of a proper outgroup. An alternate rooting method is the midpoint rooting (MPR). In this case, the root is set at the midpoint between the two most divergent operational taxonomic units. Although the midpoint rooting algorithm has been extensively used, the efficiency of this method in retrieving the correct root remains untested. In the present study, we empirically tested the success rate of the MPR in obtaining the outgroup root for a given phylogenetic tree. This was carried out by eliminating outgroups in 50 selected data sets from 33 papers and rooting the trees with the midpoint method. We were thus able to compare the root position retrieved by each method. Data sets were separated into three categories with different root consistencies: data sets with a single outgroup taxon (54% success rate for MPR), data sets with multiple outgroup taxa that showed inconsistency in root position (82% success rate), and data sets with multiple outgroup taxa in which root position was consistent (94% success rate). Interestingly, the more consistent the outgroup root is, the more successful MPR appears to be. This is a strong indication that the MPR method is valuable, particularly for cases where a proper outgroup is unavailable.  © 2007 The Linnean Society of London, Biological Journal of the Linnean Society , 2007, 92 , 669–674.     

Different gymnosperm outgroups have (mostly) congruent signal regarding the root of flowering plant phylogeny

We examined multiple plastid genes from a diversity of gymnosperm lineages to explore the consistency of signal among different outgroups for rooting flowering plant phylogeny. For maximum parsimony (MP), most outgroups attach on a branch of the underlying ingroup tree that leads to Amborella. Maximum likelihood (ML) analyses either root angiosperms on a nearby branch or find split support for these neighboring root placements, depending on the outgroup. The inclusion of two species of Hydatellaceae, recently recognized as an ancient line of angiosperms, does not aid in inference of the root. Cost profiles for placing the root in suboptimal locations are highly correlated across most outgroup comparisons, even comparing MP and ML profiles. Those for Gnetales are the most deviant of all those considered. This divergent outgroup either attaches on a long eudicot branch with moderate bootstrap support in MP analyses or supports no particular root location in ML analysis. Removing the most rapidly evolving sites in rate classifications based on two divergent angiosperm root placements with Gnetales yields strongly conflicting root placements in MP analysis, despite substantial overlap in the estimated sets of conservative sites. However, the generally high consistency in rooting signal among distantly related gymnosperm clades suggests that the long branch connecting angiosperms to their extant relatives may not interfere substantially with inference of the angiosperm root.     

Molecular phylogeny of some European heteronemertean (Nemertea) species and the monophyletic status of Riseriellus, Lineus, and Micrura

The 16S rRNA mitochondrial gene was used to reconstruct the relationships among 10 heteronemertean species (subclass Heteronemertea, phylum Nemertea); Lineus ruber and L. viridis are represented by more than one specimen to assess intraspecific variation in these enigmatic species, and the analysis includes in total 14 terminal taxa incorporating one palaeonemertean species (Tubulanus annulatus) for outgroup rooting. The aligned sequences were subjected to maximum parsimony, maximum-likelihood, and neighbor-joining analyses to estimate the phylogenetic relationship of the species. The results were concordant from all analyses and indicate that neither Lineus nor Micrura are monophyletic taxa, and that there is no support from a phylogenetic point of view to establish the monotypic genus Riseriellus.     

Random roots and lineage sorting

Lineage sorting has been suggested as a major force in generating incongruent phylogenetic signal when multiple gene partitions are examined. The degree of lineage sorting can be estimated using the coalescent process and simulation studies have also pointed to a major role for incomplete lineage sorting as a factor in phylogenetic inference. Some recent empirical studies point to an extreme role for this phenomenon with up to 50-60% of all informative genes showing incongruence as a result of lineage sorting. Here, we examine seven large multi-partition genome level data sets over a large range of taxonomic representation. We took the approach of examining outgroup choice and its impact on tree topology, by swapping outgroups into analyses with successively larger genetics distances to the ingroup. Our results indicate a linear relationship of outgroup distance with incongruence in the data sets we examined suggesting a strong random rooting effect. In addition, we attempted to estimate the degree of lineage sorting in several large genome level data sets by examining triads of very closely related taxa. This exercise resulted in much lower estimates of incongruent genes that could be the result of lineage sorting, with an overall estimate of around 10% of the total number of genes in a genome showing incongruence as a result of true lineage sorting. Finally we examined the behavior of likelihood and parsimony approaches on the random rooting phenomenon. Likelihood tends to stabilize incongruence as outgroups get further and further away from the ingroup. In one extreme case, likelihood overcompensates for sequence divergence but increases random rooting causing long branch repulsion.     

Outgroup misplacement and phylogenetic inaccuracy under a molecular clock--a simulation study

We conducted a simulation study of the phylogenetic methods UPGMA, neighbor joining, maximum parsimony, and maximum likelihood for a five-taxon tree under a molecular clock. The parameter space included a small region where maximum parsimony is inconsistent, so we tested inconsistency correction for parsimony and distance correction for neighbor joining. As expected, corrected parsimony was consistent. For these data, maximum likelihood with the clock assumption outperformed each of the other methods tested. The distance-based methods performed marginally better than did maximum parsimony and maximum likelihood without the clock assumption. Data correction was generally detrimental to accuracy, especially for short sequence lengths. We identified another region of the parameter space where, although consistent for a given method, some incorrect trees were each selected with up to twice the frequency of the correct (generating) tree for sequences of bounded length. These incorrect trees are those where the outgroup has been incorrectly placed. In addition to this problem, the placement of the outgroup sequence can have a confounding effect on the ingroup tree, whereby the ingroup is correct when using the ingroup sequences alone, but with the inclusion of the outgroup the ingroup tree becomes incorrect.     

Evolutionary relationships within the protostome phylum Sipuncula: a molecular analysis of ribosomal genes and histone H3 sequence data

The phylogenetic relationships of the members of the phylum Sipuncula are investigated by means of DNA sequence data from three nuclear markers, two ribosomal genes (18S rRNA and the D3 expansion fragment of 28S rRNA), and one protein-coding gene, histone H3. Phylogenetic analysis via direct optimization of DNA sequence data using parsimony as optimality criterion is executed for 12 combinations of parameter sets accounting for different indel costs and transversion/transition cost ratios in a sensitivity analysis framework. Alternative outgroup analyses are also performed to test whether they affected rooting of the sipunculan topology. Nodal support is measured by parsimony jackknifing and Bremer support values. Results from the different partitions are highly congruent, and the combined analysis for the parameter set that minimizes overall incongruence supports monophyly of Sipuncula, but nonmonophyly of several higher taxa recognized for the phylum. Mostly responsible for this is the split of the family Sipunculidae in three main lineages, with the genus Sipunculus being the sister group to the remaining sipunculans, the genus Phascolopsis nesting within the Golfingiiformes, and the genus Siphonosoma being associated to the Phascolosomatidea. Other interesting results are the position of Phascolion within Golfingiidae and the position of Antillesoma within Aspidosiphonidae. These results are not affected by the loci selected or by the outgroup chosen. The position of Apionsoma is discussed, although more data would be needed to better ascertain its phylogenetic affinities. Monophyly of the genera with multiple representatives (Themiste, Aspidosiphon, and Phascolosoma) is well supported, but not the monophyly of the genera Nephasoma or Golfingia. Interesting phylogeographic questions arise from analysis of multiple representatives of a few species.     

Comparative variation of morphological and molecular evolution through geologic time: 28S ribosomal RNA versus morphology in echinoids.

The comparatively good fossil record of post-Palaeozoic echinoids allows rates of morphological change to be estimated over the past 260 million years and compared with rates of molecular evolution. Parsimony analysis of morphological data, based predominantly on skeletal characteristics, and parsimony, distance and maximum likelihood analyses of molecular data, from the first 380 bases from the 5' end of the 28S rRNA molecule, for 10 species of echinoid produce congruent phylogenies. The molecular sequence chosen is demonstrably far from saturation and sister groups have divergence times ranging from about 15 to 260 Ma. Parsimony analysis allows the great majority of molecular and morphological apomorphies to be placed in one of 18 independent geological time intervals, providing a direct measure of rates of evolution for periods in the geological past. Because most molecular fixed point mutations in our sequences cannot be polarized unambiguously by outgroup comparison (making the outgroup states effectively random), distance and parsimony analyses both tend spuriously to root the echinoid tree on the longest internal branch. A topology identical to that derived from morphological data is, however, obtained using Maximum Likelihood and also parsimony analysis where outgroup rooting is restricted to more conserved regions. This is taken as the correct topology for assessing rates of evolution. Overall, both morphological and molecular changes show a moderately strong correlation with time elapsed, but a weaker correlation with one another. Statistically significant differences in evolutionary rate are found between some, but not all, pair-wise comparisons of sister lineages for both molecular and morphological data. The molecular clock rate for echinaceans is three times faster than that for cidaroids and irregular echinoids. Spearman's rank correlation test, which requires only relative magnitude of changes to be known, suggests that morphological change has a slightly better correlation with time than does molecular change, averaged over all ten species. However, when just echinaceans are considered an extremely good correlation is found between the number of molecular changes and time elapsed, whereas morphological change remains poorly correlated. Thus, molecular rates approximate to a clocklike model within restricted echinoid clades, but vary significantly between clades. Averaging results over all echinoids produces a correlation that is no better than the correlation between morphological change and time elapsed.     

Molecular Phylogeny of Some European Heteronemertean (Nemertea) Species and the Monophyletic Status ofRiseriellus, Lineus,andMicrura

The 16S rRNA mitochondrial gene was used to reconstruct the relationships among 10 heteronemertean species (subclass Heteronemertea, phylum Nemertea);Lineus ruberandL. viridisare represented by more than one specimen to assess intraspecific variation in these enigmatic species, and the analysis includes in total 14 terminal taxa incorporating one palaeonemertean species (Tubulanus annulatus) for outgroup rooting. The aligned sequences were subjected to maximum parsimony, maximum-likelihood, and neighbor-joining analyses to estimate the phylogenetic relationship of the species. The results were concordant from all analyses and indicate that neitherLineusnorMicruraare monophyletic taxa, and that there is no support from a phylogenetic point of view to establish the monotypic genusRiseriellus.     

The phylogeny and classification of Scaphopoda (Mollusca): an assessment of current resolution and cladistic reanalysis

The first cladistic analysis of phylogeny in the class Scaphopoda (Steiner 1992a,1996) examined relationships among family and selected sub-family taxa using morphological data. A preferred/ consensus tree of relationships illustrated monophyly of the orders Dentaliida and Gadilida, partial resolution among dentaliid families, and complete resolution among gadilid taxa. However, several alternative replications of the analysis, including use of a revised data matrix, did not produce the reported tree number or level of resolution; in all cases, monophyly of the Dentaliida was not supported by strict consensus of resultant parsimonious trees. Reanalysis, using unordered characters and outgroup rooting, only clearly resolves monophyly of the Gadilida and the sister relationship of the Entalinidae with the remaining gadilid families. These analyses emphasize the need for more comparative data and thorough parsimony analysis in scaphopod cladistic phylogenetics, as relationships in this class are still some way from resolution.     

Molecular phylogeny of the Siphonocladales (Chlorophyta: Cladophorophyceae)

The Siphonocladales are tropical to warm-temperate, marine green macro-algae characterized by a wide variety of thallus morphologies, ranging from branched filaments to pseudo-parenchymatous plants. Phylogenetic analyses of partial large subunit (LSU) rDNA sequences sampled from 166 isolates revealed nine well-supported siphonocladalean clades. Analyses of a concatenated dataset of small subunit (SSU) and partial LSU rDNA sequences greatly clarified the phylogeny of the Siphonocladales. However, the position of the root of the Siphonocladales could not be determined unambiguously, as outgroup rooting and molecular clock rooting resulted in a different root placement. Different phylogenetic methods (likelihood, parsimony and distance) yielded similar tree topologies with comparable internal node resolution. Likewise, analyses under more realistic models of sequence evolution, taking into account differences in evolution between stem and loop regions of rRNA, did not differ markedly from analyses using standard four-state models. The molecular phylogeny revealed that all siphonocladalean architectures may be derived from a single Cladophora-like ancestor. Parallel and convergent evolution of various morphological characters (including those traditionally employed to circumscribe the families and genera) have occurred in the Siphonocladales. Consequently, incongruence with traditional classifications, including non-monophyly in all families and most genera, was shown.     

Rooting phylogenies using gene duplications: an empirical example from the bees (Apoidea)

The placement of the root node in a phylogeny is fundamental to characterizing evolutionary relationships. The root node of bee phylogeny remains unclear despite considerable previous attention. In order to test alternative hypotheses for the location of the root node in bees, we used the F1 and F2 paralogs of elongation factor 1-alpha (EF-1α) to compare the tree topologies that result when using outgroup versus paralogous rooting. Fifty-two taxa representing each of the seven bee families were sequenced for both copies of EF-1α. Two datasets were analyzed. In the first (the "concatenated" dataset), the F1 and F2 copies for each species were concatenated and the tree was rooted using appropriate outgroups (sphecid and crabronid wasps). In the second dataset (the "duplicated" dataset), the F1 and F2 copies were aligned to each another and each copy for all taxa were treated as separate terminals. In this dataset, the root was placed between the F1 and F2 copies (e.g., paralog rooting). Bayesian analyses demonstrate that the outgroup rooting approach outperforms paralog rooting, recovering deeper clades and showing stronger support for groups well established by both morphological and other molecular data. Sequence characteristics of the two copies were compared at the amino acid level, but little evidence was found to suggest that one copy is more functionally conserved. Although neither approach yields an unambiguous root to the tree, both approaches strongly indicate that the root of bee phylogeny does not fall near Colletidae, as has been previously proposed. We discuss paralog rooting as a general strategy and why this approach performs relatively poorly with our particular dataset.     

Phylogenetic relationships ofCrepidiastrum (Asteraceae) of the Bonin (Ogasawara) Islands based on chloroplast DNA restriction site variation

RFLPs of cpDNA were investigated for seven species ofCrepidiastrum, of which three are endemic to the Bonin Islands. As an outgroup for rooting the tree,Paraixeris denticulata was also examined. Approximately 350 restriction sites were surveyed using 16 restriction enzymes. A total of 26 restriction site mutations were detected, and seven of them were phylogenetically informative. Wagner parsimony analysis resulted in four most parsimonious trees. In the tree obtained, the Bonin endemics are monophyletic with four mutations and a bootstrap value of 0.98 for the branch. This result agrees with that obtained from a previous electrophoretic analysis (Ito and Ono 1990), and supports the hypothesis that the three Bonin endemics have been evolved from a common ancestor. The Bonin endemics cluster together withC. lanceolatum, suggesting that the endemics andC. lanceolatum share a common ancestor, although the bootstrap probability is not very high (0.78) and thus other possibilities cannot be rejected.     

NUCLEIC ACID SEQUENCE PHYLOGENY AND RANDOM OUTGROUPS

Abstract— When divergent taxa are used to root networks, it is assumed that the character stales in the outgroup have historical similarity to those in the ingroup. Yet, if the data are nucleic acid sequences, the character stales shared by a divergent outgroup may be based not on history but on random similarity. A simple procedure is proposed to test this possibility. In the absence of an appropriate outgroup, root position can be estimated with the use of an asymmetrical character transformation matrix. If the matrix is sufficiently biased, it can supply the polarity information usually derived from an outgroup. This outgroup test and rooting procedure are demonstrated with ADH sequences from the genus Drosophila .     

Inferring angiosperm phylogeny from EST data with widespread gene duplication

Background

Most studies inferring species phylogenies use sequences from single copy genes or sets of orthologs culled from gene families. For taxa such as plants, with very high levels of gene duplication in their nuclear genomes, this has limited the exploitation of nuclear sequences for phylogenetic studies, such as those available in large EST libraries. One rarely used method of inference, gene tree parsimony, can infer species trees from gene families undergoing duplication and loss, but its performance has not been evaluated at a phylogenomic scale for EST data in plants.

Results

A gene tree parsimony analysis based on EST data was undertaken for six angiosperm model species and Pinus, an outgroup. Although a large fraction of the tentative consensus sequences obtained from the TIGR database of ESTs was assembled into homologous clusters too small to be phylogenetically informative, some 557 clusters contained promising levels of information. Based on maximum likelihood estimates of the gene trees obtained from these clusters, gene tree parsimony correctly inferred the accepted species tree with strong statistical support. A slight variant of this species tree was obtained when maximum parsimony was used to infer the individual gene trees instead.

Conclusion

Despite the complexity of the EST data and the relatively small fraction eventually used in inferring a species tree, the gene tree parsimony method performed well in the face of very high apparent rates of duplication.

     

Molecular Phylogeny of Neotropical Anopheles (Nyssorhynchus) Albitarsis Species Complex (Diptera: Culicidae)

A phylogeny was reconstructed for four species belonging to the Neotropical Anopheles (Nyssorhynchus) albitarsis complex using partial sequences from the mitochondrial cytochrome oxidase I (COI) and NADH dehydrogenase 4 (ND4) genes and the ribosomal DNA ITS2 and D2 expansion region of the 28S subunit. The basis for initial characterization of each member of the complex was by correlated random amplification of polymorphic DNA-polymerase chain reaction (RAPD-PCR) markers. Analyses were carried out with and without an outgroup (An.(Nys.) argyritarsis Robineau-Desvoidy) by using maximum parsimony, maximum likelihood, and Bayesian methods. A total evidence approach without the outgroup, using separate models for "fast" (COI and ND4 position 3) and "slow" (rDNA ITS2 and D2, and COI and ND4 position 1) partitions, gave the best supported topology, showing close relationships of An. albitarsis Lynch-Arribálzaga to An. albitarsis B and An. marajoara Galv?o & Damasceno to An. deaneorum Rosa-Freitas. Analyses with the outgroup included showed poorer support, possibly because of a long branch attraction effect caused by a divergent outgroup, which caused one of the An. marajoara specimens to cluster with An. deaneorum in some analyses. The relationship of the above-mentioned result to a separately proposed hypothesis suggesting a fifth species in the complex is discussed.     

Hypothetical Ancestors and Rooting in Cladistic Analysis

Most hypothetical ancestors that are used to root trees in cladistic analyses summarize character-state information in one or more outgroup taxa. Nonetheless, hypothetical ancestors also provide a means of rooting trees using the ontogenetic and paleontological methods of polarizing character transformations, and for incorporating the inferences of more than one of these methods into a single analysis. However, the use of one hypothetical ancestor that combines inferences based on outgroup comparison with those based on other methods of polarizing character transformations to root a cladogram is invalid. Inferences regarding plesiomorphic character states based on outgroup comparison apply to the outgroup node, whereas inferences based on either the ontogenetic or paleontological method apply to the ingroup node. These inferences cannot be combined into a single hypothetical construct. A hypothetical ancestor based on outgroup information is included in the data matrix and used to root the resulting network; however, because this ancestor places potentially problematic constraints on the analysis, the use of actual outgroup taxa is preferable in most instances. Correct use of a hypothetical ancestor inferred with the ontogenetic and paleontological methods involves the Lundberg method in which the shortest ingroup network is rooted at the internode to which the hypothetical ancestor attaches most parsimoniously. Because inferences of polarity based on outgroup comparison cannot be combined directly with those based on other polarization methods, the synthesis of information from all three methods in a single tree must involve taxonomic congruence.     

Molecular phylogenetics of Fouquieriaceae: evidence from nuclear rDNA ITS studies

Molecular sequence data from the 18S-26S rDNA internal transcribed spacer (ITS) region support the monophyly of Fouquieria sensu lato (Fouquieriaceae) and the three subgenera (subg. Fouquieria, subg. Bronnia, subg. Idria) previously recognized within it. Resolution within subg. Fouquieria differs somewhat between parsimony and maximum likelihood (ML) trees. Section Fouquieria and sect. Ocotilla within subg. Fouquieria are not well supported as monophyletic groups. Uncertainty regarding placement of the root within Fouquieriaceae makes discussion of character evolution within the family difficult. Three root positions are consistent with rate-constant evolution of ITS sequences: (1) along the branch to subg. Idria, (2) along the branch to subg. Bronnia, and (3) along the branch to subg. Fouquieria. The first root position listed is equivalent to an outgroup rooting. The third root position listed is equivalent to a midpoint rooting. Of the three root positions above, only the third is along a branch that may be sufficiently long to act as a long-branch attractor. The first two root positions would result in character reconstruction suggesting that succulent growth forms and white floral pigmentation are ancestral within the family, with shifts to woody growth forms and to red floral pigmentation. The third root position results in equivocal reconstruction of the ancestral growth form, equivocal reconstruction of ancestral floral pigmentation in parsimony trees, and a suggestion of white floral pigmentation as ancestral in ML trees. Two previous hypotheses of polyploid origins are compatible with the molecular data presented here: (1) origin of the tetraploid F. diguetii from F. macdougalii, and (2) allopolyploid origin of the hexaploid F. burragei from the tetraploid F. diguetii and a diploid species similar to F. splendens. Direct descent of the hexaploid F. columnaris from the subg. Bronnia lineage is not supported by our data. An amphiploid origin of F. columnaris involving a member of the subg. Bronnia lineage and an extinct taxon outside subg. Bronnia, however, cannot be ruled out.     

Optimal outgroup analysis

We present and critically examine a statistical criterion for the selection of outgroup taxa for rooting evolutionary trees. The criterion is the amount of phylogenetic signal for the ingroup when the states of the candidate outgroup taxa are assumed to be plesiomorphic relative to the ingroup for the purpose of measuring plesiomorphy content of the outgroup taxon. A statistical measure of rooted, ingroup signal was subjected to a suite of critical tests which indicate that it provides a proxy measure of plesiomorphy content. As the evolutionary distance between the ingroup ancestral node and outgroup taxa increases, the tree-independent measure of signal decreases, tracking the decay in plesiomorphy content and the increase in convergence to the ingroup states. We show that a priori generalizations about optimal outgroup taxon sampling strategies are likely to be misleading, and that testing for the suitability of available outgroup taxon sampling in specific instances is warranted. Software for optimal outgroup analysis is available.     

ONTOGENY AND CHARACTER PHYLOGENY

Abstract— The doctrine of recapitulation was long ago exchanged for a view that ontogeny is orderly and recapitulates ancestral ontogenies. Only recently has the pattern of ontogeny, as a means for determining individual character phylogenies (polarities), been explicitly explored by Gareth Nelson and other cladists. Constraints on outgroup comparisons are suggested to apply to Nelson's rule equally, and the supposed distinction between "direct" and "indirect" methods is suggested to be nonexistent. Nelson's rule is concluded to measure character adjacency directly, but polarity indirectly. Nelson's (ontogenetic) rule is compared to the Outgroup Rule of Watrous and Wheeler, based on 60 postembryonic larval characters of three species of slime-mold beetles of the genus Agathidium ( onisemdes Palisot de Beauvois, pulchrum LeConte, aristerium Wheeler), and reference to larvae of Amsotoma basalts (Le Conte) as an outgroup. Polarities were hypothesized based on each rule and cladograms constructed with the microcomputer programs Henning86 (J. S. Farris) and CLADOS (K. C. Nixon). With levels of analytical error, or homoplasy, measured as the consistency index, the results were compared and the least homoplastic solutions preferred. It is concluded that the outgroup rule and Nelson's rule arc of about the same efficacy as criteria for polarity, and that each is ultimately justified on the basis of parsimony. The following hypothesis of relationships is accepted for the species studied: ( Anisotoma basalts + ( Agathidium oniscoides + ( Agathidium pulchrum + Agathidium arislerium ))).