Should we be worried about long-branch attraction in real data sets? Investigations using metazoan 18S rDNA

Although long-branch attraction (LBA) is frequently cited as the cause of anomalous phylogenetic groupings, few examples of LBA involving real sequence data are known. We have found several cases of probable LBA by analyzing subsamples from an alignment of 18S rDNA sequences for 133 metazoans. In one example, maximum parsimony analysis of sequences from two rotifers, a ctenophore, and a polychaete annelid resulted in strong support for a tree grouping two "long-branch taxa" (a rotifer and the ctenophore). Maximum-likelihood analysis of the same sequences yielded strong support for a more biologically reasonable "rotifer monophyly" tree. Attempts to break up long branches for problematic subsamples through increased taxon sampling reduced, but did not eliminate, LBA problems. Exhaustive analyses of all quartets for a subset of 50 sequences were performed in order to compare the performance of maximum likelihood, equal-weights parsimony, and two additional variants of parsimony; these methods do differ substantially in their rates of failure to recover trees consistent with well established, but highly unresolved phylogenies. Power analyses using simulations suggest that some incorrect inferences by maximum parsimony are due to statistical inconsistency and that when estimates of central branch lengths for certain quartets are very low, maximum-likelihood analyses have difficulty recovering accepted phylogenies even with large amounts of data. These examples demonstrate that LBA problems can occur in real data sets, and they provide an opportunity to investigate causes of incorrect inferences.     

Evolutionary history of "early-diverging" eukaryotes: the excavate taxon Carpediemonas is a close relative of Giardia

Diplomonads, such as Giardia, and their close relatives retortamonads have been proposed as early-branching eukaryotes that diverged before the acquisition-retention of mitochondria, and they have become key organisms in attempts to understand the evolution of eukaryotic cells. In this phylogenetic study we focus on a series of eukaryotes suggested to be relatives of diplomonads on morphological grounds, the "excavate taxa". Phylogenies of small subunit ribosomal RNA (SSU rRNA) genes, alpha-tubulin, beta-tubulin, and combined alpha- + beta-tubulin all scatter the various excavate taxa across the diversity of eukaryotes. But all phylogenies place the excavate taxon Carpediemonas as the closest relative of diplomonads (and, where data are available, retortamonads). This novel relationship is recovered across phylogenetic methods and across various taxon-deletion experiments. Statistical support is strongest under maximum-likelihood (ML) (when among-site rate variation is modeled) and when the most divergent diplomonad sequences are excluded, suggesting a true relationship rather than an artifact of long-branch attraction. When all diplomonads are excluded, our ML SSU rRNA tree actually places retortamonads and Carpediemonas away from the base of the eukaryotes. The branches separating excavate taxa are mostly not well supported (especially in analyses of SSU rRNA data). Statistical tests of the SSU rRNA data, including an "expected likelihood weights" approach, do not reject trees where excavate taxa are constrained to be a clade (with or without parabasalids and Euglenozoa). Although diplomonads and retortamonads lack any mitochondria-like organelle, Carpediemonas contains double membrane-bounded structures physically resembling hydrogenosomes. The phylogenetic position of Carpediemonas suggests that it will be valuable in interpreting the evolutionary significance of many molecular and cellular peculiarities of diplomonads.     

Molecular Phylogeny of Conjugating Green Algae (Zygnemophyceae,Streptophyta) Inferred from SSU rDNA Sequence Comparisons

Nuclear-encoded SSU rDNA sequences have been obtained from 64 strains of conjugating green algae (Zygnemophyceae, Streptophyta, Viridiplantae). Molecular phylogenetic analyses of 90 SSU rDNA sequences of Viridiplantae (inciuding 78 from the Zygnemophyceae) were performed using complex evolutionary models and maximum likelihood, distance, and maximum parsimony methods. The significance of the results was tested by bootstrap analyses, deletion of long-branch taxa, relative rate tests, and Kishino-Hasegawa tests with user-defined trees. All results support the monophyly of the class Zygnemophyceae and of the order Desmidiales. The second order, Zygnematales, forms a series of early-branching clades in paraphyletic succession, with the two traditional families Mesotaeniaceae and Zygnemataceae not recovered as lineages. Instead, a long-branch Spirogyra/Sirogonium clade and the later-diverging Netrium and Roya clades represent independent clades. Within the order Desmidiales, the families Gonatozygaceae and Closteriaceae are monophyletic, whereas the Peniaceae (represented only by Penium margaritaceum) and the Desmidiaceae represent a single weakly supported lineage. Within the Desmidiaceae short internal branches and varying rates of sequence evolution among taxa reduce the phylogenetic resolution significantly. The SSU rDNA-based phylogeny is largely congruent with a published analysis of the rbcL phylogeny of the Zygnemophyceae (McCourt et al. 2000) and is also in general agreement with classification schemes based on cell wall ultrastructure. The extended taxon sampling at the subgenus level provides solid evidence that many genera in the Zygnemophyceae are not monophyletic and that the genus concept in the group needs to be revised.     

Revised small subunit rRNA analysis provides further evidence that Foraminifera are related to Cercozoa

There is accumulating evidence that the general shape of the ribosomal DNA-based phylogeny of Eukaryotes is strongly biased by the long-branch attraction phenomenon, leading to an artifactual basal clustering of groups that are probably highly derived. Among these groups, Foraminifera are of particular interest, because their deep phylogenetic position in ribosomal trees contrasts with their Cambrian appearance in the fossil record. A recent actin-based phylogeny of Eukaryotes has proposed that Foraminifera might be closely related to Cercozoa and, thus, branch among the so-called crown of Eukaryotes. Here, we reanalyze the small-subunit ribosomal RNA gene (SSU rDNA) phylogeny by removing all long-branching lineages that could artifactually attract foraminiferan sequences to the base of the tree. Our analyses reveal that Foraminifera branch together with the marine testate filosean Gromia oviformis as a sister group to Cercozoa, in agreement with actin phylogeny. Our study confirms the utility of SSU rDNA as a phylogenetic marker of megaevolutionary history, provided that the artifacts due to the heterogeneity of substitution rates in ribosomal genes are circumvented.     

Global eukaryote phylogeny: Combined small- and large-subunit ribosomal DNA trees support monophyly of Rhizaria, Retaria and Excavata

Resolution of the phylogenetic relationships among the major eukaryotic groups is one of the most important problems in evolutionary biology that is still only partially solved. This task was initially addressed using a single marker, the small-subunit ribosomal DNA (SSU rDNA), although in recent years it has been shown that it does not contain enough phylogenetic information to robustly resolve global eukaryotic phylogeny. This has prompted the use of multi-gene analyses, especially in the form of long concatenations of numerous conserved protein sequences. However, this approach is severely limited by the small number of taxa for which such a large number of protein sequences is available today. We have explored the alternative approach of using only two markers but a large taxonomic sampling, by analysing a combination of SSU and large-subunit (LSU) rDNA sequences. This strategy allows also the incorporation of sequences from non-cultivated protists, e.g., Radiozoa (=radiolaria minus Phaeodarea). We provide the first LSU rRNA sequences for Heliozoa, Apusozoa (both Apusomonadida and Ancyromonadida), Cercozoa and Radiozoa. Our Bayesian and maximum likelihood analyses for 91 eukaryotic combined SSU+LSU sequences yielded much stronger support than hitherto for the supergroup Rhizaria (Cercozoa plus Radiozoa plus Foraminifera) and several well-recognised groups and also for other problematic clades, such as the Retaria (Radiozoa plus Foraminifera) and, with more moderate support, the Excavata. Within opisthokonts, the combined tree strongly confirms that the filose amoebae Nuclearia are sisters to Fungi whereas other Choanozoa are sisters to animals. The position of some bikont taxa, notably Heliozoa and Apusozoa, remains unresolved. However, our combined trees suggest a more deeply diverging position for Ancyromonas, and perhaps also Apusomonas, than for other bikonts, suggesting that apusozoan zooflagellates may be central for understanding the early evolution of this huge eukaryotic group. Multiple protein sequences will be needed fully to resolve basal bikont phylogeny. Nonetheless, our results suggest that combined SSU+LSU rDNA phylogenies can help to resolve several ambiguous regions of the eukaryotic tree and identify key taxa for subsequent multi-gene analyses.     

Nuclear and Nucleomorph SSU rDNA Phylogeny in the Cryptophyta and the Evolution of Cryptophyte Diversity

The plastid-bearing members of the Cryptophyta contain two functional eukaryotic genomes of different phylogenetic origin, residing in the nucleus and in the nucleomorph, respectively. These widespread and diverse protists thus offer a unique opportunity to study the coevolution of two different eukaryotic genomes within one group of organisms. In this study, the SSU rRNA genes of both genomes were PCR-amplified with specific primers and phylogenetic analyses were performed on different data sets using different evolutionary models. The results show that the composition of the principal clades obtained from the phylogenetic analyses of both genes was largely congruent, but striking differences in evolutionary rates were observed. These affected the topologies of the nuclear and nucleomorph phylogenies differently, resulting in long-branch attraction artifacts when simple evolutionary models were applied. Deletion of long-branch taxa stabilized the internal branching order in both phylogenies and resulted in a completely resolved topology in the nucleomorph phylogeny. A comparison of the tree topologies derived from SSU rDNA sequences with characters previously used in cryptophyte systematics revealed that the biliprotein type was congruent, but the type of inner periplast component incongruent, with the molecular trees. The latter is indicative of a hidden cellular dimorphism (cells with two periplast types present in a single clonal strain) of presumably widespread occurrence throughout cryptophyte diversity, which, in consequence, has far-reaching implications for cryptophyte systematics as it is practiced today.     

Dinoflagellate phylogeny as inferred from heat shock protein 90 and ribosomal gene sequences

Background

Interrelationships among dinoflagellates in molecular phylogenies are largely unresolved, especially in the deepest branches. Ribosomal DNA (rDNA) sequences provide phylogenetic signals only at the tips of the dinoflagellate tree. Two reasons for the poor resolution of deep dinoflagellate relationships using rDNA sequences are (1) most sites are relatively conserved and (2) there are different evolutionary rates among sites in different lineages. Therefore, alternative molecular markers are required to address the deeper phylogenetic relationships among dinoflagellates. Preliminary evidence indicates that the heat shock protein 90 gene (Hsp90) will provide an informative marker, mainly because this gene is relatively long and appears to have relatively uniform rates of evolution in different lineages.

Methodology/Principal Findings

We more than doubled the previous dataset of Hsp90 sequences from dinoflagellates by generating additional sequences from 17 different species, representing seven different orders. In order to concatenate the Hsp90 data with rDNA sequences, we supplemented the Hsp90 sequences with three new SSU rDNA sequences and five new LSU rDNA sequences. The new Hsp90 sequences were generated, in part, from four additional heterotrophic dinoflagellates and the type species for six different genera. Molecular phylogenetic analyses resulted in a paraphyletic assemblage near the base of the dinoflagellate tree consisting of only athecate species. However, Noctiluca was never part of this assemblage and branched in a position that was nested within other lineages of dinokaryotes. The phylogenetic trees inferred from Hsp90 sequences were consistent with trees inferred from rDNA sequences in that the backbone of the dinoflagellate clade was largely unresolved.

Conclusions/Significance

The sequence conservation in both Hsp90 and rDNA sequences and the poor resolution of the deepest nodes suggests that dinoflagellates reflect an explosive radiation in morphological diversity in their recent evolutionary past. Nonetheless, the more comprehensive analysis of Hsp90 sequences enabled us to infer phylogenetic interrelationships of dinoflagellates more rigorously. For instance, the phylogenetic position of Noctiluca, which possesses several unusual features, was incongruent with previous phylogenetic studies. Therefore, the generation of additional dinoflagellate Hsp90 sequences is expected to refine the stem group of athecate species observed here and contribute to future multi-gene analyses of dinoflagellate interrelationships.     

Molecular Phylogeny of Conjugating Green Algae (Zygnemophyceae, Streptophyta) Inferred from SSU rDNA Sequence Comparisons

Abstract Nuclear-encoded SSU rDNA sequences have been obtained from 64 strains of conjugating green algae (Zygnemophyceae, Streptophyta, Viridiplantae). Molecular phylogenetic analyses of 90 SSU rDNA sequences of Viridiplantae (inciuding 78 from the Zygnemophyceae) were performed using complex evolutionary models and maximum likelihood, distance, and maximum parsimony methods. The significance of the results was tested by bootstrap analyses, deletion of long-branch taxa, relative rate tests, and Kishino–Hasegawa tests with user-defined trees. All results support the monophyly of the class Zygnemophyceae and of the order Desmidiales. The second order, Zygnematales, forms a series of early-branching clades in paraphyletic succession, with the two traditional families Mesotaeniaceae and Zygnemataceae not recovered as lineages. Instead, a long-branch Spirogyra/Sirogonium clade and the later-diverging Netrium and Roya clades represent independent clades. Within the order Desmidiales, the families Gonatozygaceae and Closteriaceae are monophyletic, whereas the Peniaceae (represented only by Penium margaritaceum) and the Desmidiaceae represent a single weakly supported lineage. Within the Desmidiaceae short internal branches and varying rates of sequence evolution among taxa reduce the phylogenetic resolution significantly. The SSU rDNA-based phylogeny is largely congruent with a published analysis of the rbcL phylogeny of the Zygnemophyceae (McCourt et al. 2000) and is also in general agreement with classification schemes based on cell wall ultrastructure. The extended taxon sampling at the subgenus level provides solid evidence that many genera in the Zygnemophyceae are not monophyletic and that the genus concept in the group needs to be revised.     

Intragenomic variation in the ITS rDNA region obscures phylogenetic relationships and inflates estimates of operational taxonomic units in genus Laetiporus

Regions of rDNA are commonly used to infer phylogenetic relationships among fungal species and as DNA barcodes for identification. These regions occur in large tandem arrays, and concerted evolution is believed to reduce intragenomic variation among copies within these arrays, although some variation still might exist. Phylogenetic studies typically use consensus sequencing, which effectively conceals most intragenomic variation, but cloned sequences containing intragenomic variation are becoming prevalent in DNA databases. To understand effects of using cloned rDNA sequences in phylogenetic analyses we amplified and cloned the ITS region from pure cultures of six Laetiporus species and one Wolfiporia species (Basidiomycota, Polyporales). An average of 66 clones were selected randomly and sequenced from 21 cultures, producing a total of 1399 interpretable sequences. Significant variation (≥ 5% variation in sequence similarity) was observed among ITS copies within six cultures from three species clades (L. cincinnatus, L. sp. clade J, and Wolfiporia dilatohypha) and phylogenetic analyses with the cloned sequences produced different trees relative to analyses with consensus sequences. Cloned sequences from L. cincinnatus fell into more than one species clade and numerous cloned L. cincinnatus sequences fell into entirely new clades, which if analyzed on their own most likely would be recognized as "undescribed" or "novel" taxa. The use of a 95% cut off for defining operational taxonomic units (OTUs) produced seven Laetiporus OTUs with consensus ITS sequences and 20 OTUs with cloned ITS sequences. The use of cloned rDNA sequences might be problematic in fungal phylogenetic analyses, as well as in fungal bar-coding initiatives and efforts to detect fungal pathogens in environmental samples.     

Turning the crown upside down: gene tree parsimony roots the eukaryotic tree of life

The first analyses of gene sequence data indicated that the eukaryotic tree of life consisted of a long stem of microbial groups "topped" by a crown-containing plants, animals, and fungi and their microbial relatives. Although more recent multigene concatenated analyses have refined the relationships among the many branches of eukaryotes, the root of the eukaryotic tree of life has remained elusive. Inferring the root of extant eukaryotes is challenging because of the age of the group (～1.7-2.1 billion years old), tremendous heterogeneity in rates of evolution among lineages, and lack of obvious outgroups for many genes. Here, we reconstruct a rooted phylogeny of extant eukaryotes based on minimizing the number of duplications and losses among a collection of gene trees. This approach does not require outgroup sequences or assumptions of orthology among sequences. We also explore the impact of taxon and gene sampling and assess support for alternative hypotheses for the root. Using 20 gene trees from 84 diverse eukaryotic lineages, this approach recovers robust eukaryotic clades and reveals evidence for a eukaryotic root that lies between the Opisthokonta (animals, fungi and their microbial relatives) and all remaining eukaryotes.     

Unusually high evolutionary rate of the elongation factor 1 alpha genes from the Ciliophora and its impact on the phylogeny of eukaryotes

The elongation factor 1 alpha (EF-1 alpha) has become widely employed as a phylogenetic marker for studying eukaryotic evolution. However, a disturbing problem, the artifactual polyphyly of ciliates, is always observed. It has been suggested that the addition of new sequences will help to circumvent this problem. Thus, we have determined 15 new ciliate EF-1 alpha sequences, providing for a more comprehensive taxonomic sampling of this phylum. These sequences have been analyzed together with a representation of eukaryotic sequences using distance-, parsimony-, and likelihood-based phylogenetic methods. Such analyses again failed to recover the monophyly of Ciliophora. A study of the substitution rate showed that ciliate EF-1 alpha genes exhibit a high evolutionary rate, produced in part by an increased number of variable positions. This acceleration could be related to alterations of the accessory functions acquired by this protein, likely to those involving interactions with the cytoskeleton, which is very modified in the Ciliophora. The high evolutionary rate of these sequences leads to an artificial basal emergence of some ciliates in the eukaryotic tree by effecting a long-branch attraction artifact that produces an asymmetric topology for the basal region of the tree. The use of a maximum-likelihood phylogenetic method (which is less sensitive to long-branch attraction) and the addition of sequences to break long branches allow retrieval of more symmetric topologies, which suggests that the asymmetric part of the tree is most likely artifactual. Therefore, the sole reliable part of the tree appears to correspond to the apical symmetric region. These kinds of observations suggest that the general eukaryotic evolution might have consisted of a massive radiation followed by an increase in the evolutionary rates of certain groups that emerge artificially as early branches in the asymmetric base of the tree. Ciliates in the case of the EF-1 alpha genes would offer clear evidence for this hypothesis.     

Is the Felsenstein zone a fly trap?

Although long-branch attraction, the incorrect grouping of long lineages in a phylogeny because of systematic error, has been identified as a potential source of error in phylogenetic analysis for almost two decades, no empirical examples of the phenomenon exist. Here, I outline several criteria for identifying long-branch attraction and apply these criteria to 18S ribosomal DNA (rDNA) sequence data for 13 insects. Parsimony and minimum evolution with p distances group the two longest branches together (those leading to Strepsiptera and Diptera). Simulation studies show that the long branches are long enough to attract. When a tree is assumed in which Strepsiptera and Diptera are separated and many data sets are simulated for that tree (using the parameter estimates for that tree for the original data), parsimony analysis of the simulated data consistently groups Strepsiptera and Diptera. Analyses of the 18S rDNA sequences using methods that are less sensitive to the problem of long-branch attraction estimate trees in which the long branches are separate.     

A molecular phylogeny of annelids

We present parsimony analyses of annelids based on the largest taxon sample and most extensive molecular data set yet assembled, with two nuclear ribosomal genes (18S rDNA and the D1 region of 28S rDNA), one nuclear protein coding‐gene (Histone H3) and one mitochondrial ribosomal gene (16S rDNA) from 217 terminal taxa. Of these, 267 sequences are newly sequenced, and the remaining were obtained from GenBank. The included taxa are based on the criteria that the taxon must have 18S rDNA or at least two other loci. Our analyses show that 68% of annelid family ranked taxa represented by more than one taxon in our study are supported by a jackknife value > 50%. In spite of the size of our data set, the phylogenetic signal in the deepest part of the tree remains weak and the majority of the currently recognized major polychaete clades (except Amphinomida and Aphroditiformia) could not be recovered. Terbelliformia is monophyletic (with the exclusion of Pectinariidae, for which only 18S data were available), whereas members of taxa such as Phyllodocida, Cirratuliformia, Sabellida and Scolecida are scattered over the trees. Clitellata is monophyletic, although Dinophilidae should possibly be included, and Clitellata has a sister group within the polychaetes. One major problem is the current lack of knowledge on the closest relatives to annelids and the position of the annelid root. We suggest that the poor resolution in the basal parts of the trees presented here may be due to lack of signal connected to incomplete data sets both in terms of terminal and gene sampling, rapid radiation events and/or uneven evolutionary rates and long‐branch attraction. © The Willi Hennig Society 2006.     

Early diverging Ascomycota: phylogenetic divergence and related evolutionary enigmas

The early diverging Ascomycota lineage, detected primarily from nSSU rDNA sequence-based phylogenetic analyses, includes enigmatic key taxa important to an understanding of the phylogeny and evolution of higher fungi. At the moment six representative genera of early diverging ascomycetes (i.e. Taphrina, Protomyces, Saitoella, Schizosaccharomyces, Pneumocystis and Neolecta) have been assigned to "Archiascomycetes" sensu Nishida and Sugi ama (1994) or the subphylum "Taphrinomycotina" sensu Eriksson and Winka (1997). The group includes fungi that are ecologically and morphologically diverse, and it is difficult therefore to define the group based on common phenotypic characters. Bayesian analyses of nSSU rDNA or combined nSSU and nLSU rDNA sequences supported previously published Ascomycota frameworks that consist of three major lineages (i.e. a group of early diverging Ascomycota. [Taphrinomycotina], Saccharomycotina and Pezizomycotina); Taphrinomycotina is the sister group of Saccharomycotina and Pezizomycotina. The 50% majority rule consensus of 18000 Bayesian MCMCMC-generated trees from multilocus gene sequences of nSSU rDNA, nLSU rDNA (D1/D2), RPB2 and beta-tubulin also showed the monophyly of the three subphyla and the basal position of Taphrinomycotina in Ascomycota with significantly higher statistical support. However to answer controversial questions on the origin, monophyly and evolution of the Taphrinomycotina, additional integrated phylogenetic analyses might be necessary using sequences of more genes with broader taxon sampling from the early diverging Ascomycota.     

Long-Branch Abstractions

Recent attention has been focused on the sensitivities of various tree reconstructing algorithms to sequence rate heterogeneity (long-branch attraction). Phylogenetic conclusions from two recent empirical studies have been indicted as artifacts attributable to long-branch attraction. Siddall et al. (1995) concluded that Myxozoa are cnidarians and sister group to Polypodium based on 18S rDNA and morphology. Hanelt et al. (1996) argued that this result is due to long-branch attraction. Whiting et al. (1997) concluded that the Strepsiptera are sister group to Diptera based on parsimony analysis of 18S rDNA, 28S rDNA, and morphology. Huelsenbeck (1997) argued that this result also is attributable to long-branch attraction. We demonstrate that the analyses and arguments dismissing these results as the effects of long-branch attraction are fundamentally flawed. The criteria employed by these authors were applied arbitrarily by them to the groups that they did not want, and yet using those same criteria, there is more reason to exclude other taxa besides Polypodium and there is more reason to disbelieve monophyly of Diptera than monophyly of Strepsiptera with Diptera. Moreover, it is asserted, long-branch attraction cannot explain the presence of nematocysts in Myxozoa and halteres in Strepsiptera. For these reasons, and in light of the demonstration that long branches cannot attract each other in their mutual absence, we conclude that the monophyly of Myxozoa + Polypodium and Strepsiptera + Diptera is not due to long-branch attraction. We suggest that maximum likelihood methods are extremely sensitive to taxon and character sampling and that these data sets are demonstrative of the long-branch repulsion problem.     

Relationships among seed plants inferred from highly conserved genes: sorting conflicting phylogenetic signals among ancient lineages

Phylogenetic studies based on different types and treatment of data provide substantially conflicting hypotheses of relationships among seed plants. We conducted phylogenetic analyses of sequences of two highly conserved chloroplast genes, psaA and psbB, for a comprehensive taxonomic sample of seed plants and land plants. Parsimony analyses of two different codon position partitions resulted in well-supported, but significantly conflicting, phylogenetic trees. First and second codon positions place angiosperms and gymnosperms as sister clades and Gnetales as sister to Pinaceae. Third positions place Gnetales as sister to all other seed plants. Maximum likelihood trees for the two partitions are also in conflict. Relationships among the main seed plant clades according to first and second positions are similar to those found in parsimony analysis for the same data, but the third position maximum likelihood tree is substantially different from the corresponding parsimony tree, although it agrees partially with the first and second position trees in placing Gnetales as the sister group of Pinaceae. Our results document high rate heterogeneity among lineages, which, together with the greater average rate of substitution for third positions, may reduce phylogenetic signal due to long-branch attraction in parsimony reconstructions. Whereas resolution of relationships among major seed plant clades remains pending, this study provides increased support for relationships within major seed plant clades.     

The phylogeny of land plants inferred from 18S rDNA sequences: pushing the limits of rDNA signal?

Previous studies of the phylogeny of land plants based on analysis of 18S ribosomal DNA (rDNA) sequences have generally found weak support for the relationships recovered and at least some obviously spurious relationships, resulting in equivocal inferences of land plant phylogeny. We hypothesized that greater sampling of both characters and taxa would improve inferences of land plant phylogeny based on 18S rDNA sequences. We therefore conducted a phylogenetic analysis of complete (or nearly complete) 18S rDNA sequences for 93 species of land plants and 7 green algal relatives. Parsimony analyses with equal weighting of characters and characters state changes and parsimony analyses weighting (1) stem bases half as much as loop bases and (2) transitions half as much as transversions did not produce substantially different topologies. Although the general structure of the shortest trees is consistent with most hypotheses of land plant phylogeny, several relationships, particularly among major groups of land plants, appear spurious. Increased character and taxon sampling did not substantially improve the performance of 18S rDNA in phylogenetic analyses of land plants, nor did analyses designed to accommodate variation in evolutionary rates among sites. The rate and pattern of 18S rDNA evolution across land plants may limit the usefulness of this gene for phylogeny reconstruction at deep levels of plant phylogeny. We conclude that the mosaic structure of 18S rDNA, consisting of highly conserved and highly variable regions, may contain historical signal at two levels. Rapidly evolving regions are informative for relatively recent divergences (e.g., within angiosperms, seed plants, and ferns), but homoplasy at these sites makes it difficult to resolve relationships among these groups. At deeper levels, changes in the highly conserved regions of small-subunit rDNAs provide signal across all of life. Because constraints imposed by the secondary structure of the rRNA may affect the phylogenetic information content of 18S rDNA, we suggest that 18S rDNA sequences be combined with other data and that methods of analysis be employed to accommodate these differences in evolutionary patterns, particularly across deep divergences in the tree of life.     

Phylogenetic inference regarding Parergodrilidae and Hrabeiella periglandulata ('Polychaeta', Annelida) based on 18S rDNA, 28S rDNA and COI sequences

Parergodrilidae and Hrabeiella periglandulata are Annelida showing different combinations of clitellate-like and aclitellate characters. Similarities between both of these taxa and Clitellata have widely been regarded as the result of convergent evolution due to similar selection pressures. The position of the three taxa in the phylogenetic system of Annelida is still in debate. However, in analyses based on 18S rDNA sequences a close relationship of Parergodrilidae with Orbiniidae and Questidae was suggested. To infer their phylogeny the sequences of the 28S rDNA and of the cytochrome oxidase I (COI) gene of Stygocapitella subterranea , Parergodrilus heideri and H. periglandulata were determined. The data were extended by sequences of various species including species from Clitellata and Orbiniidae. Prior to tree reconstruction the dataset was analysed in detail for phylogenetic content by applying a sliding window analysis, a likelihood mapping and Modeltest V.3.04. Subsequently, generalized parsimony and maximum likelihood methods were employed. Clade robustness was estimated by bootstrapping. In addition, combined analyses of the sequences of 18S rDNA and 28S rDNA as well as of 18S rDNA, 28S rDNA and COI were performed. The combination of the data of the two structure genes and a mitochondrial gene improved the resolution obtained with the single datasets slightly. These analyses support a close relationship of Parergodrilidae and Orbiniidae but cannot resolve the position of H. periglandulata . In every analysis Clitellata cluster within 'Polychaeta', confirming previous investigations.     

Escaping from the Felsenstein Zone by Detecting Long Branches in Phylogenetic Data

Long branches in a true phylogeny tend to disrupt hierarchical character covariation (phylogenetic signal) in the distribution of traits among organisms. The distortion of hierarchical structure in character-state matrices can lead to errors in the estimation of phylogenetic relationships and inconsistency of methods of phylogenetic inference. Examination of trees distorted by long-branch attraction will not reveal the identities of problematic taxa, in part because the distortion can mask long branches by reducing inferred branch lengths and through errors in branching order. Here we present a simple method for the detection of taxa whose placement in evolutionary trees is made difficult by the effects of long-branch attraction. The method is an extension of a tree-independent conceptual framework of phylogenetic data exploration (RASA). Taxa that are likely to attract are revealed because long branches leave distinct footprints in the distribution of character states among taxa, and these traces can be directly observed in the error structure of the RASA regression. Problematic taxa are identified using a new diagnostic plot called the taxon variance plot, in which the apparent cladistic and phenetic variances contributed by individual taxa are compared. The procedure for identifying long edges employs algorithms solved in polynomial time and can be applied to morphological, molecular, and mixed characters. The efficacy of the method is demonstrated using simulated evolution and empirical evidence of long branches in a set of recently published sequences. We show that the accuracy of evolutionary trees can be improved by detecting and combating the potentially misleading influences of long-branch taxa.     

Covarion shifts cause a long-branch attraction artifact that unites microsporidia and archaebacteria in EF-1alpha phylogenies

Microsporidia branch at the base of eukaryotic phylogenies inferred from translation elongation factor 1alpha (EF-1alpha) sequences. Because these parasitic eukaryotes are fungi (or close relatives of fungi), it is widely accepted that fast-evolving microsporidian sequences are artifactually "attracted" to the long branch leading to the archaebacterial (outgroup) sequences ("long-branch attraction," or "LBA"). However, no previous studies have explicitly determined the reason(s) why the artifactual allegiance of microsporidia and archaebacteria ("M + A") is recovered by all phylogenetic methods, including maximum likelihood, a method that is supposed to be resistant to classical LBA. Here we show that the M + A affinity can be attributed to those alignment sites associated with large differences in evolutionary site rates between the eukaryotic and archaebacterial subtrees. Therefore, failure to model the significant evolutionary rate distribution differences (covarion shifts) between the ingroup and outgroup sequences is apparently responsible for the artifactual basal position of microsporidia in phylogenetic analyses of EF-1alpha sequences. Currently, no evolutionary model that accounts for discrete changes in the site rate distribution on particular branches is available for either protein or nucleotide level phylogenetic analysis, so the same artifacts may affect many other "deep" phylogenies. Furthermore, given the relative similarity of the site rate patterns of microsporidian and archaebacterial EF-1alpha proteins ("parallel site rate variation"), we suggest that the microsporidian orthologs may have lost some eukaryotic EF-1alpha-specific nontranslational functions, exemplifying the extreme degree of reduction in this parasitic lineage.