Mitochondrial <Emphasis Type="Italic">matR</Emphasis> sequences help to resolve deep phylogenetic relationships in rosids

Background  

Rosids are a major clade in the angiosperms containing 13 orders and about one-third of angiosperm species. Recent molecular analyses recognized two major groups (i.e., fabids with seven orders and malvids with three orders). However, phylogenetic relationships within the two groups and among fabids, malvids, and potentially basal rosids including Geraniales, Myrtales, and Crossosomatales remain to be resolved with more data and a broader taxon sampling. In this study, we obtained DNA sequences of the mitochondrial matR gene from 174 species representing 72 families of putative rosids and examined phylogenetic relationships and phylogenetic utility of matR in rosids. We also inferred phylogenetic relationships within the "rosid clade" based on a combined data set of 91 taxa and four genes including matR, two plastid genes (rbcL, atpB), and one nuclear gene (18S rDNA).     

Persistent nuclear ribosomal DNA sequence polymorphism in the Amelanchier agamic complex (Rosaceae)

Individual plants of several Amelanchier taxa contain many polymorphic nucleotide sites in the internal transcribed spacers (ITS) of nuclear ribosomal DNA (nrDNA). This polymorphism is unusual because it is not recent in origin and thus has resisted homogenization by concerted evolution. Amelanchier ITS sequence polymorphism is hypothesized to be the result of gene flow between two major North American clades resolved by phylogenetic analysis of ITS sequences. Western North American species plus A. humilis and A. sanguinea of eastern North America form one clade (A), and the remaining eastern North American Amelanchier make up clade B. Five eastern North American taxa are polymorphic at many of the nucleotide sites where clades A and B have diverged and are thought to be of hybrid origin, with A. humilis or A. sanguinea as one parent and various members of clade B as the other parent. Morphological evidence suggests that A. humilis is one of the parents of one of the polymorphic taxa, a microspecies that we refer to informally as A. "erecta." Sequences of 21 cloned copies of the ITS1- 5.8S gene-ITS2 region from one A. "erecta" individual are identical to A. humilis sequence or to the clade B consensus sequence, or they are apparent recombinants of A. humilis and clade B ITS repeats. Amelanchier "erecta" and another polymorphic taxon are suspected to be relatively old because both grow several hundred kilometers beyond the range of one of their parents. ITS sequence polymorphisms have apparently persisted in these two taxa perhaps because of polyploidy and/or agamospermy (asexual seed production), which are prevalent in the genus.      

A decade of progress in plant molecular phylogenetics

Over the past decade, botanists have produced several thousand phylogenetic analyses based on molecular data, with particular emphasis on sequencing rbcL, the plastid gene encoding the large subunit of Rubisco (ribulose bisphosphate carboxylase). Because phylogenetic trees retrieved from the three plant genomes (plastid, nuclear and mitochondrial) have been highly congruent, the ‘Angiosperm Phylogeny Group’ has used these DNA-based phylogenetic trees to reclassify all families of flowering plants. However, in addition to taxonomy, these major phylogenetic efforts have also helped to define strategies to reconstruct the ‘tree of life’, and have revealed the size of the ancestral plant genome, uncovered potential candidates for the ancestral flower, identified molecular living fossils, and linked the rate of neutral substitutions with species diversity. With an increased interest in DNA sequencing programmes in non-model organisms, the next decade will hopefully see these phylogenetic findings integrated into new genetic syntheses, from genomes to taxa.     

Piecing together the "new" Plantaginaceae

Scrophulariaceae is one of the families that has been divided extensively due to the results of DNA sequence studies. One of its segregates is a vastly enlarged Plantaginaceae. In a phylogenetic study of 47 members of Plantaginaceae and seven outgroups based on 3561 aligned characters from four DNA regions (the nuclear ribosomal ITS region and the plastid trnL-F, rps16 intron, and matK-trnK intron regions), the relationships within this clade were analyzed. The results from parsimony and Bayesian analyses support the removal of the Lindernieae from Gratioleae to a position outside Plantaginaceae. A group of mainly New World genera is paraphyletic with respect to a clade of Old World genera. Among the New World taxa, those offering oil as a pollinator reward cluster together. Ourisia is sister to this clade. Gratioleae consist of Gratiola, Otacanthus, Bacopa, Stemodia, Scoparia, and Mecardonia. Cheloneae plus Russelia and Tetranema together constitute the sister group to a clade predominantly composed of Old World taxa. Among the Old World clade, Ellisiophyllum and Lafuentea have been analyzed for the first time in a molecular phylogenetic analysis. The former genus is sister to Sibthorpia and the latter is surprisingly the sister to Antirrhineae.     

Resolution of a concatenation/coalescence kerfuffle: partitioned coalescence support and a robust family‐level tree for Mammalia

Recent phylogenetic analyses of a large dataset for mammalian families (169 taxa, 26 loci) portray contrasting results. Supermatrix (concatenation) methods support a generally robust tree with only a few inconsistently resolved polytomies, whereas MP‐EST coalescence analysis of the same dataset yields a weakly supported tree that conflicts with many traditionally recognized clades. Here, we evaluate this discrepancy via improved coalescence analyses with reference to the rich history of phylogenetic studies on mammals. This integration clearly demonstrates that both supermatrix and coalescence analyses of just 26 loci yield a congruent, well‐supported phylogenetic hypothesis for Mammalia. Discrepancies between published studies are explained by implementation of overly simple DNA substitution models, inadequate tree‐search routines and limitations of the MP‐EST method. We develop a simple measure, partitioned coalescence support (PCS), which summarizes the distribution of support and conflict among gene trees for a given clade. Extremely high PCS scores for outlier gene trees at two nodes in the mammalian tree indicate a troubling bias in the MP‐EST method. We conclude that in this age of phylogenomics, a solid understanding of systematics fundamentals, choice of valid methodology and a broad knowledge of a clade's taxonomic history are still required to yield coherent phylogenetic inferences.     

Tree of life for the genera of Chinese vascular plants

We reconstructed a phylogenetic tree of Chinese vascular plants (Tracheophyta) using sequences of the chloroplast genes atpB, matK, ndhF, and rbcL and mitochondrial matR. We produced a matrix comprising 6098 species and including 13?695 DNA sequences, of which 1803 were newly generated. Our taxonomic sampling spanned 3114 genera representing 323 families of Chinese vascular plants, covering more than 93% of all genera known from China. The comprehensive large phylogeny supports most relationships among and within families recognized by recent molecular phylogenetic studies for lycophytes, ferns (monilophytes), gymnosperms, and angiosperms. For angiosperms, most families in Angiosperm Phylogeny Group IV are supported as monophyletic, except for a paraphyletic Dipterocarpaceae and Santalaceae. The infrafamilial relationships of several large families and monophyly of some large genera are well supported by our dense taxonomic sampling. Our results showed that two species of Eberhardtia are sister to a clade formed by all other taxa of Sapotaceae, except Sarcosperma. We have made our phylogeny of Chinese vascular plants publically available for the creation of subtrees via SoTree (http://www.darwintree.cn/flora/index.shtml), an automated phylogeny assembly tool for ecologists.     

Phylogeny,divergence times and biogeography of window flies (Scenopinidae) and the therevoid clade (Diptera: Asiloidea)

The evolution of the ‘therevoid’ clade, with an emphasis on window flies (Scenopinidae), is presented by combining DNA sequence data with morphological characters for living and fossil species. The therevoid clade represents a group of four families (Apsilocephalidae, Evocoidae, Scenopinidae and Therevidae) of lower brachyceran Diptera in the superfamily Asiloidea. A comprehensive phylogenetic analysis using parsimony and likelihood methods was undertaken using extensive taxon sampling from all families and subfamilies, and compared with outgroup taxa sampled from the related families Asilidae, Mydidae, Apioceridae and Empididae. Fifty‐nine morphological characters (adult, larval and pupal) were combined with 6.4 kb of DNA sequences for two ribosomal genes (16S and 18S ribosomal DNA) and three protein‐encoding genes [cytochrome oxidase I (COI), triose phosphate isomerase (TPI) and the CPSase region of carbamoyl‐phosphate synthase‐aspartate transcarbamoylase‐dihydroorotase (CAD)]. Results from combined analyses of morphological and molecular data for 78 taxa representing all families of the therevoid clade are presented. Specific hypotheses of the relationship between respective families and subfamilies were tested statistically using four‐cluster likelihood mapping. The therevoid clade is a well‐supported monophyletic group within Asiloidea, with Evocoidae sister to Apsilocephalidae and Therevidae sister to Scenopinidae. Temporal and zoogeographical aspects of therevoid clade evolution were investigated using Bayesian divergence time estimates and Lagrange ancestral range scenarios. The effect of inclusion of fossils as terminal taxa on phylogenetic and divergence time estimation was investigated, with morphological scoring for fossil representatives included in the analyses rather than used simply as minimum age constraints. In each analysis there was either improvement in estimation, or only marginal and localized loss in tree resolution, and with younger estimates of divergence time across the tree. The historical biogeography of the therevoid clade was examined with multiple trans‐Antarctic vicariance events between Australasia and South America evident during the Late Cretaceous to early Palaeogene. Scenopininae is newly subdivided into two tribes, Metatrichini trib.n. and Scenopinini Fallén stat.r. This published work has been registered in ZooBank, http://zoobank.org/urn:lsid:zoobank.org:pub:4974EBF8‐3117‐4189‐B6DE‐7D5BF9B23E53 .     

Phylogenetic Placement of Pentatomid Stink Bug Gut Symbionts

Insect bacterial symbionts are ubiquitous, however, only a few groups of host families have been well studied in relation to their associations with microbes. The determination of the phylogenetic relationships among bacteria associated with different species within an insect family can provide insights into the biology and evolution of these interactions. We studied the phylogenetic placement of vertically transmitted bacterial symbionts associated with the posterior midgut (crypt-bearing) region of pentatomid stink bugs (Hemiptera, Pentatomidae). Our results demonstrate that different host species carried one major bacterium in their midgut. Phylogenetic analyses of the 16S rRNA gene sequences obtained from the midgut of stink bugs placed all symbionts in a clade with Erwinia and Pantoea species, both plant-associated bacteria. Results indicate that symbiont monophyly occurs among recently diverged taxa (e.g., within a genus) but does not occur in the Pentatomidae. Results suggest that these vertically transmitted symbionts are occasionally replaced by other taxonomically similar bacteria over evolutionary time. Our findings highlight how the evolutionary history of hemipteran symbionts in unexplored host families may have unpredictable levels of complexity.     

Inter-generic relationships of the crows, jays, magpies and allied groups (Aves: Corvidae) based on nucleotide sequence data

Phylogenetic relationships were studied based on DNA sequences obtained from all recognized genera of the family Corvidae sensu stricto . The aligned data set consists 2589 bp obtained from one mitochondrial and two nuclear genes. Maximum parsimony, maximum-likelihood, and Bayesian inference analyses were used to estimate phylogenetic relationships. The analyses were done for each gene separately, as well as for all genes combined. An analysis of a taxonomically expanded data set of cytochrome b sequences was performed in order to infer the phylogenetic positions of six genera for which nuclear genes could not be obtained. Monophyly of the Corvidae is supported by all analyses, as well as by the occurrence of a deletion of 16 bp in the β-fibrinogen intron in all ingroup taxa. Temnurus and Pyrrhocorax are placed as the sister group to all other corvids, while Cissa and Urocissa appear as the next clade inside them. Further up in the tree, two larger and well-supported clades of genera were recovered by the analyses. One has an entirely New World distribution (the New World jays), while the other includes mostly Eurasian (and one African) taxa. Outside these two major clades are Cyanopica and Perisoreus whose phylogenetic positions could not be determined by the present data. A biogeographic analysis of our data suggests that the Corvidae underwent an initial radiation in Southeast Asia. This is consistent with the observation that almost all basal clades in the phylogenetic tree consist of species adapted to tropical and subtropical forest habitats.     

VITCOMIC2: visualization tool for the phylogenetic composition of microbial communities based on 16S rRNA gene amplicons and metagenomic shotgun sequencing

Background

The 16S rRNA gene-based amplicon sequencing analysis is widely used to determine the taxonomic composition of microbial communities. Once the taxonomic composition of each community is obtained, evolutionary relationships among taxa are inferred by a phylogenetic tree. Thus, the combined representation of taxonomic composition and phylogenetic relationships among taxa is a powerful method for understanding microbial community structure; however, applying phylogenetic tree-based representation with information on the abundance of thousands or more taxa in each community is a difficult task. For this purpose, we previously developed the tool VITCOMIC (VIsualization tool for Taxonomic COmpositions of MIcrobial Community), which is based on the genome-sequenced microbes’ phylogenetic information. Here, we introduce VITCOMIC2, which incorporates substantive improvements over VITCOMIC that were necessary to address several issues associated with 16S rRNA gene-based analysis of microbial communities.

Results

We developed VITCOMIC2 to provide (i) sequence identity searches against broad reference taxa including uncultured taxa; (ii) normalization of 16S rRNA gene copy number differences among taxa; (iii) rapid sequence identity searches by applying the graphics processing unit-based sequence identity search tool CLAST; (iv) accurate taxonomic composition inference and nearly full-length 16S rRNA gene sequence reconstructions for metagenomic shotgun sequencing; and (v) an interactive user interface for simultaneous representation of the taxonomic composition of microbial communities and phylogenetic relationships among taxa. We validated the accuracy of processes (ii) and (iv) by using metagenomic shotgun sequencing data from a mock microbial community.

Conclusions

The improvements incorporated into VITCOMIC2 enable users to acquire an intuitive understanding of microbial community composition based on the 16S rRNA gene sequence data obtained from both metagenomic shotgun and amplicon sequencing.

     

Phylogenetic utility of the nuclear gene arginine decarboxylase: an example from Brassicaceae

Arginine decarboxylase (ADC) is an important enzyme in the production of putrescine and polyamines in plants. It is encoded by a single or low-copy nuclear gene that lacks introns in sequences studied to date. The rate of Adc amino acid sequence evolution is similar to that of ndhF for the angiosperm family studied. Highly conserved regions provide several target sites for PCR priming and sequencing and aid in nucleotide and amino acid sequence alignment across a range of taxonomic levels, while a variable region provides an increased number of potentially informative characters relative to ndhF for the taxa surveyed. The utility of the Adc gene in plant molecular systematic studies is demonstrated by analysis of its partial nucleotide sequences obtained from 13 representatives of Brassicaceae and 3 outgroup taxa, 2 from the mustard oil clade (order Capparales) and 1 from the related order Malvales. Two copies of the Adc gene, Adc1 and Adc2, are found in all members of the Brassicaceae studied to data except the basal genus Aethionema. The resulting Adc gene tree provides robust phylogenetic data regarding relationships within the complex mustard family, as well as independent support for proposed tribal realignments based on other molecular data sets such as those from chloroplast DNA.      

Taxon sampling to address an ancient rapid radiation: a supermatrix phylogeny of early brachyceran flies (Diptera)

Early diverging brachyceran fly lineages underwent a rapid radiation approximately 180 Ma, coincident in part with the origin of flowering plants. This region of the fly tree includes 25 000 described extant species with diverse ecological roles such as blood‐feeding (haematophagy), parasitoidism, predation, pollination and wood‐feeding (xylophagy). Early diverging brachyceran lineages were once considered a monophyletic group of families called Orthorrhapha, based on the shared character of a longitudinal break in the pupal skin made during the emergence of the adult. Yet other morphological and molecular evidence generally supports a paraphyletic arrangement of ‘Orthorrhapha’, with strong support for one orthorrhaphan lineage – dance flies and relatives – as the closest relative to all higher flies (Cyclorrhapha), together called Eremoneura. In order to establish a comprehensive estimate of the relationships among orthorrhaphan lineages using a thorough sample of publicly available data, we compiled and analysed a dataset including 1217 taxa representing major lineages and 20 molecular markers. Our analyses suggest that ‘Orthorrhapha’ excluding Eremoneura is not monophyletic; instead, we recover two main lineages of early brachyceran flies: Homeodactyla and Heterodactyla. Homeodactyla includes Nemestrinoidea (uniting two parasitic families Acroceridae + Nemestrinidae) as the closest relatives to the large SXT clade, comprising Stratiomyomorpha, Xylophagidae and Tabanomorpha. Heterodactyla includes Bombyliidae with a monophyletic Asiloidea (exclusive of Bombyliidae) as the closest relatives to Eremoneura. Reducing missing data, modifying the distribution of genes across taxa, and, in particular, removing rogue taxa significantly improved tree resolution and statistical support. Although our analyses rely on dense taxonomic sampling and substantial gene coverage, our results pinpoint the limited resolving power of Sanger sequencing‐era molecular phylogenetic datasets with respect to ancient, hyperdiverse radiations.     

The Evolutionary History of Sarco(endo)plasmic Calcium ATPase (SERCA)

Investigating the phylogenetic relationships within physiologically essential gene families across a broad range of taxa can reveal the key gene duplication events underlying their family expansion and is thus important to functional genomics studies. P-Type II ATPases represent a large family of ATP powered transporters that move ions across cellular membranes and includes Na⁺/K⁺ transporters, H⁺/K⁺ transporters, and plasma membrane Ca²⁺ pumps. Here, we examine the evolutionary history of one such transporter, the Sarco(endo)plasmic reticulum calcium ATPase (SERCA), which maintains calcium homeostasis in the cell by actively pumping Ca²⁺ into the sarco(endo)plasmic reticulum. Our protein-based phylogenetic analyses across Eukaryotes revealed two monophyletic clades of SERCA proteins, one containing animals, fungi, and plants, and the other consisting of plants and protists. Our analyses suggest that the three known SERCA proteins in vertebrates arose through two major gene duplication events after the divergence from tunicates, but before the separation of fishes and tetrapods. In plants, we recovered two SERCA clades, one being the sister group to Metazoa and the other to Apicomplexa clade, suggesting an ancient duplication in an early eukaryotic ancestor, followed by subsequent loss of one copy in Opisthokonta, the other in protists, and retention of both in plants. We also report relatively recent and independent gene duplication events within invertebrate taxa including tunicates and the leech Helobdella robusta. Thus, it appears that both ancient and recent gene duplication events have played an important role in the evolution of this ubiquitous gene family across the eukaryotic domain.     

Nuclear and mitochondrial sequence data reveal the major lineages of starlings, mynas and related taxa

We investigated the phylogenetic relationships among the major lineages of the avian family Sturnidae and their placement within the Muscicapoidea clade using two nuclear (RAG-1 and myoglobin) and one mitochondrial gene (ND2). Among Muscicapoidea, we recovered three clades corresponding to the families Cinclidae, Muscicapidae and Sturnidae (sensu [Sibley, C.G., Monroe Jr., B.L., 1990. Distribution and Taxonomy of Birds of the World. Yale University Press, New Haven, CT]). Within the sturnoid lineage Mimini and Sturnini are sister groups, with Buphagus basal to them. We identified three major lineages of starlings: the Philippine endemic genus Rhabdornis, an Oriental-Australasian clade (genera Scissirostrum, Gracula, Mino, Ampeliceps, Sarcops, Aplonis), and an Afrotropical-Palaearctic clade (all African taxa, Sturnus and Acridotheres). We discuss the biogeographic implications of our findings and suggest an Asiatic origin for this family. The congruence between the age of major clades, estimated by NPRS, and palaeoclimatic data present evidence for the role of climatic changes in shaping present day distribution of the group.     

Exclusive conservation of mitochondrial group II intron nad4i548 among liverworts and its use for phylogenetic studies in this ancient plant clade

Liverworts occupy a pivotal position in land plant (embryophyte) phylogeny as the presumed earliest-branching major clade, sister to all other land plants, including the mosses, hornworts, lycophytes, monilophytes and seed plants. Molecular support for this earliest dichotomy in land plant phylogeny comes from strikingly different occurrences of introns in mitochondrial genes distinguishing liverworts from all other embryophytes. Exceptionally, however, the nad5 gene--the mitochondrial locus hitherto used most widely to elucidate early land plant phylogeny--carries a group I type intron that is shared between liverworts and mosses. We here explored whether a group II intron, the other major type of organellar intron, would similarly be conserved in position across the entire diversity of extant liverworts and could be of use for phylogenetic analyses in this supposedly most ancient embryophyte clade. To this end, we investigated the nad4 gene as a candidate locus possibly featuring different introns in liverworts as opposed to the non-liverwort embryophyte (NLE) lineage. We indeed found group II intron nad4i548 universally conserved in a wide phylogenetic sampling of 55 liverwort taxa, confirming clade specificity and surprising evolutionary stability of plant mitochondrial introns. As expected, intron nad4i548g2 carries phylogenetic information in its variable sequences, which confirms and extends previous cladistic insights on liverwort evolution. We integrate the new nad4 data with those of the previously established mitochondrial nad5 and the chloroplast rbcL and rps4 genes and present a phylogeny based on the fused datasets. Notably, the phylogenetic analyses suggest a reconsideration of previous phylogenetic and taxonomic assignments for the genera Calycularia and Mylia and resolve a sister group relationship of Ptilidiales and Porellales.     

Fingerprinting the Asterid species using subtracted diversity array reveals novel species-specific sequences

Background

Asterids is one of the major plant clades comprising of many commercially important medicinal species. One of the major concerns in medicinal plant industry is adulteration/contamination resulting from misidentification of herbal plants. This study reports the construction and validation of a microarray capable of fingerprinting medicinally important species from the Asterids clade.

Methodology/Principal Findings

Pooled genomic DNA of 104 non-asterid angiosperm and non-angiosperm species was subtracted from pooled genomic DNA of 67 asterid species. Subsequently, 283 subtracted DNA fragments were used to construct an Asterid-specific array. The validation of Asterid-specific array revealed a high (99.5%) subtraction efficiency. Twenty-five Asterid species (mostly medicinal) representing 20 families and 9 orders within the clade were hybridized onto the array to reveal its level of species discrimination. All these species could be successfully differentiated using their hybridization patterns. A number of species-specific probes were identified for commercially important species like tea, coffee, dandelion, yarrow, motherwort, Japanese honeysuckle, valerian, wild celery, and yerba mate. Thirty-seven polymorphic probes were characterized by sequencing. A large number of probes were novel species-specific probes whilst some of them were from chloroplast region including genes like atpB, rpoB, and ndh that have extensively been used for fingerprinting and phylogenetic analysis of plants.

Conclusions/Significance

Subtracted Diversity Array technique is highly efficient in fingerprinting species with little or no genomic information. The Asterid-specific array could fingerprint all 25 species assessed including three species that were not used in constructing the array. This study validates the use of chloroplast genes for bar-coding (fingerprinting) plant species. In addition, this method allowed detection of several new loci that can be explored to solve existing discrepancies in phylogenetics and fingerprinting of plants.     

Phylogeny of Myrmeleontiformia based on larval morphology (Neuropterida: Neuroptera)

The suborder Myrmeleontiformia is a derived lineage of lacewings (Insecta: Neuroptera) including the families Psychopsidae, Nemopteridae, Nymphidae, Ascalaphidae and Myrmeleontidae. In particular, Myrmeleontidae (antlions) are the most diverse neuropteran family, representing a conspicuous component of the insect fauna of xeric environments. We present the first detailed quantitative phylogenetic analysis of Myrmeleontiformia, based on 107 larval morphological and behavioural characters for 36 genera whose larvae are known (including at least one representative of all the subfamilies of the suborder). Four related families were used as outgroups to polarize character states. Phylogenetic analyses were conducted using both parsimony and Bayesian methods. The reconstructions resulting from our analyses corroborate the monophyly of Myrmeleontiformia. Within this clade, Psychopsidae are recovered as the sister family to all the remaining taxa. Nemopteridae (including both subfamilies Nemopterinae and Crocinae) are recovered as monophyletic and sister to the clade comprising Nymphidae + (Myrmeleontidae + Ascalaphidae). Nymphidae consist of two well‐supported clades corresponding to the subfamilies Nymphinae and Myiodactylinae. Our results suggest that Ascalaphidae may not be monophyletic, as they collapse into an unresolved polytomy under the Bayesian analysis. In addition, the recovered phylogenetic relationships diverge from the traditional classification scheme for ascalaphids. Myrmeleontidae are reconstructed as monophyletic, with the subfamilies Stilbopteryginae, Palparinae and Myrmeleontinae. We retrieved a strongly supported clade comprising taxa with a fossorial habit of the preimaginal instars, which represents a major antlion radiation, also including the monophyletic pit‐trap building species.     

Supermatrix data highlight the phylogenetic relationships of photosynthetic stramenopiles

Molecular data had consistently recovered monophyletic classes for the heterokont algae, however, the relationships among the classes had remained only partially resolved. Furthermore, earlier studies did not include representatives from all taxonomic classes. We used a five-gene (nuclear encoded SSU rRNA; plastid encoded rbcL, psaA, psbA, psbC) analysis with a subset of 89 taxa representing all 16 heterokont classes to infer a phylogenetic tree. There were three major clades. The Aurearenophyceae, Chrysomerophyceae, Phaeophyceae, Phaeothamniophyceae, Raphidophyceae, Schizocladiophyceae and Xanthophyceae formed the SI clade. The Chrysophyceae, Eustigmatophyceae, Pinguiophyceae, Synchromophyceae and Synurophyceae formed the SII clade. The Bacillariophyceae, Bolidophyceae, Dictyochophyceae and Pelagophyceae formed the SIII clade. These three clades were also found in a ten-gene analysis. The approximately unbiased test rejected alternative hypotheses that forced each class into either of the other two clades. Morphological and biochemical data were not available for all 89 taxa, however, existing data were consistent with the molecular phylogenetic tree, especially for the SIII clade.