The phylogeny and evolution of host choice in the Hippoboscoidea (Diptera) as reconstructed using four molecular markers

Hippoboscoidea is a superfamily of Diptera that contains the Glossinidae or tsetse flies, the Hippoboscidae or louse flies, and two families of bat flies, the Streblidae and the Nycteribiidae. We reconstruct the phylogenetic relationships within Hippoboscoidea using maximum parsimony and Bayesian methods based on nucleotide sequences from fragments of four genes: nuclear 28S ribosomal DNA and the CPSase domain of CAD, and mitochondrial 16S rDNA and cytochrome oxidase I. We recover monophyly for most of the presently recognized groups within Hippoboscoidea including the superfamily as a whole, the Hippoboscidae, the Nycteribiidae, the bat flies, and the Pupipara (=Hippoboscidae+Nycteribiidae+Streblidae), as well as several subfamilies within the constituent families. Streblidae appear to be paraphyletic. Our phylogenetic hypothesis is well supported and decisive in that most competing topological hypotheses for the Hippoboscoidea require significantly longer trees. We confirm a single shift from a free-living fly to a blood-feeding ectoparasite of vertebrates and demonstrate that at least two host shifts from mammals to birds have occurred. Wings have been repeatedly lost, but never regained. The hippoboscoid ancestor also evolved adenotrophic viviparity and our cladogram is consistent with a gradual reduction in the motility of the deposited final instar larvae from active burrowing in the soil to true pupiparity where adult females glue the puparium within the confines of bat roosts.     

Higher level phylogeny and the first divergence time estimation of Heteroptera (Insecta: Hemiptera) based on multiple genes

Heteroptera, or true bugs, are the largest, morphologically diverse and economically important group of insects with incomplete metamorphosis. However, the phylogenetic relationships within Heteroptera are still in dispute and most of the previous studies were based on morphological characters or with single gene (partial or whole 18S rDNA). Besides, so far, divergence time estimates for Heteroptera totally rely on the fossil record, while no studies have been performed on molecular divergence rates. Here, for the first time, we used maximum parsimony (MP), maximum likelihood (ML) and Bayesian inference (BI) with multiple genes (18S rDNA, 28S rDNA, 16S rDNA and COI) to estimate phylogenetic relationships among the infraorders, and meanwhile, the Penalized Likelihood (r8s) and Bayesian (BEAST) molecular dating methods were employed to estimate divergence time of higher taxa of this suborder. Major results of the present study included: Nepomorpha was placed as the most basal clade in all six trees (MP trees, ML trees and Bayesian trees of nuclear gene data and four-gene combined data, respectively) with full support values. The sister-group relationship of Cimicomorpha and Pentatomomorpha was also strongly supported. Nepomorpha originated in early Triassic and the other six infraorders originated in a very short period of time in middle Triassic. Cimicomorpha and Pentatomomorpha underwent a radiation at family level in Cretaceous, paralleling the proliferation of the flowering plants. Our results indicated that the higher-group radiations within hemimetabolous Heteroptera were simultaneously with those of holometabolous Coleoptera and Diptera which took place in the Triassic. While the aquatic habitat was colonized by Nepomorpha already in the Triassic, the Gerromorpha independently adapted to the semi-aquatic habitat in the Early Jurassic.     

Influence of data conflict and molecular phylogeny of major clades in Cimicomorphan true bugs (Insecta: Hemiptera: Heteroptera)

Cimicomorpha, which consists of 16 families representing more than 19,400 species, is the largest infraorder in Heteroptera, Insecta. We present the first molecular phylogenetic investigation of family relationships of Cimicomorpha, including 46 taxa from 12 of 16 Cimicomorphan families. Three genes, with a total of 3277 bp of sequence data (nuclear 18S rDNA: 2022 bp, 28S rDNA: 755 bp, and mitochondrial 16S rDNA: 498 bp) were analyzed. Data partitions were analyzed separately and in combination, by employing ML (maximum likelihood), MP (maximum parsimony), and Bayesian methods. As saturation was detected in substitutions of 16S rDNA, influence of data conflict in combined analyses was further explored by three methods: the incongruence length difference (ILD) test, the partitioned Bremer support (PBS), and the partition addition bootstrap alteration approach (PABA). PBS and PABA approaches suggested that 16S rDNA was not very suitable for addressing relationships at this level in Cimicomorpha. Our results also supported the nabid-cimicoid lineage for Cimicoidea proposed by Schuh and Stys [Schuh, R.T., Stys, P., 1991. Phylogenetic analysis of Cimicomorphan family relationships (Heteroptera). J. NY Entomol. Soc. 99 (3), 298-350]. Data incongruence and the utility of the three genes were briefly discussed.     

Distributional Evidence for Cospeciation between Neotropical Bats and their Bat Fly Ectoparasites

Distributional evidence bearing on cospeciation of New World bats (superfamily Noctilionoidea) and associated bat flies (Diptera: Streblidae) was analyzed and revealed substantial phylogenetic structure, especially when we focused on "normal" host associations and encoded the data set hierarchically. Brooks' parsimony analyses of bat fly occurrences among hosts recovered generally accepted sets of host relationships. Putatively monophyletic host taxa provide a necessary but insufficient condition for coevolution via cospeciation.     

Identification and phylogenetic analysis of Arsenophonus- and Photorhabdus-type bacteria from adult Hippoboscidae and Streblidae (Hippoboscoidea)

This is the first report of Arsenophonus- and Photorhabdus-type bacteria from Streblidae (bat flies) and Hippoboscidae (louse flies, keds). Strains were detected by means of polymerase chain reaction of 16S rDNA, and phylogenetic analysis determined the relationship of the obtained sequences to previously reported sequences in GenBank. Phylogenetic analysis by means of maximum parsimony revealed that all isolated Arsenophonus spp. 16S rDNA sequences formed a monophyletic sub-clade within other insect Arsenophonus spp., while the Photorhabdus spp. sequences are part of a monophyletic clade including Photorhabdus spp., Xenorhabdus spp. and Proteus spp.     

Phylogenetics of barbets (Aves: Piciformes) based on nuclear and mitochondrial DNA sequence data

A combination of nuclear (beta-fibrinogen intron 7; 938 bp) and mitochondrial (cytochrome b; 1045 bp) DNA sequence data analyzed with model-based phylogenetic methods yields a hypothesis of barbet and toucan relationships supported by high Bayesian posterior probabilities and several synapomorphic indels in the nuclear intron data. The basal topology differs from previous morphology and mitochondrial DNA sequence based hypotheses, but is consistent with DNA-DNA hybridization results. The mitochondrial DNA sequence data provide resolution at the tips of the tree, but when analyzed alone, yield a different basal topology than the combined data. However, the basal nodes in the mitochondrial-based tree have little statistical support. Combined data analysis produced strong statistical support at basal nodes and a relatively simple geographic structure. Barbets from each of the three tropical regions are monophyletic, but the Old World barbets are paraphyletic. The African and New World clades are sister taxa, and the Asian clade is basal. This geographic structure indicates that similarities in plumage, voice, and behavior between Gymnobucco and Calorhamphus are convergent. The data are not conclusive, but suggest paraphyly of the New World barbets. Only 4% of the Bayesian posterior distribution unites Semnornis with the other New World barbets.     

Phylogeny of Anophelinae (Diptera: Culicidae) based on nuclear ribosomal and mitochondrial DNA sequences

Abstract Phylogenetic relationships among thirty-two species of mosquitoes in subfamily Anophelinae are inferred from portions of the mitochondrial genes COI and COII, the nuclear 18S small subunit rRNA gene and the expansion D2 region of the nuclear large subunit 28S rRNA gene. Sequences were obtained from the genera Anopheles , Bironella and Chagasia . Representatives of all six subgenera of Anopheles were included: Anopheles , Cellia , Kerteszia , Lophopodomyia , Nyssorhynchus and Stethomyia. Using parsimony and maximum likelihood methods, various combinations of these DNA sequence data were analysed separately: 18S, 28S, combined 18S and 28S, combined COI and COII, and combined 18S, 28S, COI and COII ('total evidence'). The combined rDNA data contain strong phylogenetic signal, moderately to strongly supporting most clades in MP and ML analyses; however, the mtDNA data (analysed as either nucleotide or amino acid sequences) contain little phylogenetic signal, except for relationships of very recently derived groups of species and, at the deepest level, for the monophyly of Anophelinae. The paraphyly of Anopheles relative to Bironella is confirmed by most analyses and statistical tests. Support for the monophyly of subgenera Anopheles , Cellia , Kerteszia and Nyssorhynchus is indicated by most analyses. Subgenus Lophopodomyia is reconstructed as the sister to Bironella , nested within a clade also containing Nyssorhynchus and Kerteszia . The most basal relationships within genus Anopheles are not well resolved by any of the data partitions, although the results of statistical analyses of the rDNA data (S-H-tests, likelihood ratio tests for monophyly and Bayesian MCMC analyses) suggest that the clade consisting of Bironella , Lophopodomyia , Nyssorhynchus and Kerteszia is the sister to the clade containing Cellia and Anopheles .     

Phylogenetic relationships of land plants using mitochondrial small-subunit rDNA sequences

Phylogenetic relationships among embryophytes (tracheophytes, mosses, liverworts, and hornworts) were examined using 21 newly generated mitochondrial small-subunit (19S) rDNA sequences. The "core" 19S rDNA contained more phylogenetically informative sites and lower homoplasy than either nuclear 18S or plastid 16S rDNA. Results of phylogenetic analyses using parsimony (MP) and likelihood (ML) were generally congruent. Using MP, two trees were obtained that resolved either liverworts or hornworts as the basal land plant clade. The optimal ML tree showed hornworts as basal. That topology was not statistically different from the two MP trees, thus both appear to be equally viable evolutionary hypotheses. High bootstrap support was obtained for the majority of higher level embryophyte clades named in a recent morphologically based classification, e.g., Tracheophyta, Euphyllophytina, Lycophytina, and Spermatophytata. Strong support was also obtained for the following monophyletic groups: hornworts, liverworts, mosses, lycopsids, leptosporangiate and eusporangiate ferns, gymnosperms and angiosperms. This molecular analysis supported a sister relationship between Equisetum and leptosporangiate ferns and a monophyletic gymnosperms sister to angiosperms. The topologies of deeper clades were affected by taxon inclusion (particularly hornworts) as demonstrated by jackknife analyses. This study represents the first use of mitochondrial 19S rDNA for phylogenetic purposes and it appears well-suited for examining intermediate to deep evolutionary relationships among embryophytes.     

Molecular systematics and biogeography of the southern South american freshwater "crabs" Aegla (decapoda: Anomura: Aeglidae) using multiple heuristic tree search approaches

Recently new heuristic genetic algorithms such as Treefinder and MetaGA have been developed to search for optimal trees in a maximum likelihood (ML) framework. In this study we combined these methods with other standard heuristic approaches such as ML and maximum parsimony hill-climbing searches and Bayesian inference coupled with Markov chain Monte Carlo techniques under homogeneous and mixed models of evolution to conduct an extensive phylogenetic analysis of the most abundant and widely distributed southern South American freshwater"crab,"the Aegla(Anomura: Aeglidae). A total of 167 samples representing 64 Aegla species and subspecies were sequenced for one nuclear (28S rDNA) and four mitochondrial (12S and 16S rDNA, COI, and COII) genes (5352 bp total). Additionally, six other anomuran species from the genera Munida,Pachycheles, and Uroptychus(Galatheoidea), Lithodes(Paguroidea), and Lomis(Lomisoidea) and the nuclear 18S rDNA gene (1964 bp) were included in preliminary analyses for rooting the Aegla tree. Nonsignificantly different phylogenetic hypotheses resulted from all the different heuristic methods used here, although the best scored topologies found under the ML hill-climbing, Bayesian, and MetaGA approaches showed considerably better likelihood scores (Delta> 54) than those found under the MP and Treefinder approaches. Our trees provided strong support for most of the recognized Aegla species except for A. cholchol,A. jarai,A. parana,A. marginata, A. platensis, and A. franciscana, which may actually represent multiple species. Geographically, the Aegla group was divided into a basal western clade (21 species and subspecies) composed of two subclades with overlapping distributions, and a more recent central-eastern clade (43 species) composed of three subclades with fairly well-recognized distributions. This result supports the Pacific-Origin Hypothesis postulated for the group; alternative hypotheses of Atlantic or multiple origins were significantly rejected by our analyses. Finally, we combined our phylogenetic results with previous hypotheses of South American paleodrainages since the Jurassic to propose a biogeographical framework of the Aegla radiation.     

A phylogeny of the "evil tribe" (Vernonieae: Compositae) reveals Old/New World long distance dispersal: support from separate and combined congruent datasets (trnL-F, ndhF, ITS)

The Vernonieae is one of the major tribes of the largest family of flowering plants, the sunflower family (Compositae or Asteraceae), with ca. 25,000 species. While the family's basal members (the Barnadesioideae) are found in South America, the tribe Vernonieae originated in the area of southern Africa/Madagascar. Its sister tribe, the Liabeae, is New World, however. This is the only such New/Old World sister tribe pairing anywhere in the family. The Vernonieae is now found on islands and continents worldwide and includes more than 1500 taxa. The Vernonieae has been called the "evil tribe" because overlapping character states make taxonomic delimitations difficult at all levels from the species to the subtribe for the majority of taxa. Juxtaposed with these difficult-to-separate entities are monotypic genera with highly distinctive morphologies and no obvious affinities to any other members of the tribe. The taxonomic frustration generated by these contrary circumstances has resulted in a lack of any phylogeny for the tribe until now. A combined approach using DNA sequence data from two chloroplast regions, the ndhF gene and the noncoding spacer trnL-F, and from the nuclear rDNA ITS region for 90 taxa from throughout the world was used to reconstruct the evolutionary history of the tribe. The data were analyzed separately and in combination using maximum parsimony (MP), minimum evolution neighbor-joining (NJ), and Bayesian analysis, the latter producing the best resolved and most strongly supported tree. In general, the phylogeny shows Old World taxa to be basal and New World taxa to be derived, but this is not always the case. Old and New World species are found together in two separate and only distantly related clades. This is best explained by long-distance dispersal with a minimum of two trans-oceanic exchanges. Meso/Central America has had an important role in ancient dispersals between the Old and New World and more recent movements from South to North America in the New World.     

Phylogeny and zoogeography of six squirrel species of the genus sciurus (mammalia, rodentia), inferred from cytochrome B gene sequences

To investigate the phylogenetic relationships between the New World Sciurus and the Old World Sciurus and their biogeographic history, the partial mitochondrial cytochrome b gene sequences (1,040 base pairs) were analyzed on six Sciurus species: S. aberti, S. carolinensis, S. lis, S. niger, S. stramineus, and S. vulgaris. Phylogenetic trees (maximum parsimony, neighbor-joining, and maximum likelihood methods) commonly showed two groups with high bootstrap values (73-100%): one consisting of the New World Sciurus and the other consisting of the Old World Sciurus. Genetic distances among the New World Sciurus species were remarkably larger than that between two Sciurus species of the Old World, suggesting the earlier radiation of the New World Sciurus than the Old World Sciurus.     

Bayesian mixed models and the phylogeny of pitvipers (Viperidae: Serpentes)

The subfamily Crotalinae (pitvipers) contains over 190 species of venomous snakes distributed in both the Old and New World. We incorporated an extensive sampling of taxa (including 28 of 29 genera), and sequences of four mitochondrial gene fragments (2.3kb) per individual, to estimate the phylogeny of pitvipers based on maximum parsimony and Bayesian phylogenetic methods. Our Bayesian analyses incorporated complex mixed models of nucleotide evolution that allocated independent models to various partitions of the dataset within combined analyses. We compared results of unpartitioned versus partitioned Bayesian analyses to investigate how much unpartitioned (versus partitioned) models were forced to compromise estimates of model parameters, and whether complex models substantially alter phylogenetic conclusions to the extent that they appear to extract more phylogenetic signal than simple models. Our results indicate that complex models do extract more phylogenetic signal from the data. We also address how differences in phylogenetic results (e.g., bipartition posterior probabilities) obtained from simple versus complex models may be interpreted in terms of relative credibility. Our estimates of pitviper phylogeny suggest that nearly all recently proposed generic reallocations appear valid, although certain Old and New World genera (Ovophis, Trimeresurus, and Bothrops) remain poly- or paraphyletic and require further taxonomic revision. While a majority of nodes were resolved, we could not confidently estimate the basal relationships among New World genera and which lineage of Old World species is most closely related to this New World group.     

Is 16S rDNA a Reliable Phylogenetic Marker to Characterize Relationships Below the Family Level in the Enterobacteriaceae?

The phylogenetic relationships of multiple enterobacterial species were reconstructed based on 16S rDNA gene sequences to evaluate the robustness of this housekeeping gene in the taxonomic placement of the enteric plant pathogens Erwinia, Brenneria, Pectobacterium, and Pantoea. Four data sets were compiled, two of which consisted of previously published data. The data sets were designed in order to evaluate how 16S rDNA gene phylogenies are affected by the use of different plant pathogen accessions and varying numbers of animal pathogen and outgroup sequences. DNA data matrices were analyzed using maximum likelihood (ML) algorithms, and character support was determined by ML bootstrap and Bayesian analyses. As additional animal pathogen sequences were added to the phylogenetic analyses, taxon placement changed. Further, the phylogenies varied in their placement of the plant pathogen species, and only the genus Pantoea was monophyletic in all four trees. Finally, bootstrap and Bayesian support values were low for most of the nodes, and all nonterminal branches collapsed in strict consensus trees. Inspection of 16S rDNA nucleotide alignments revealed several highly variable blocks punctuated by regions of conserved sequence. These data suggest that 16S rDNA, while effective for both species-level and family-level phylogenetic reconstruction, may underperform for genus-level phylogenetic analyses in the Enterobacteriaceae.     

Higher level phylogenetic relationships of euascomycetes (Pezizomycotina) inferred from a combined analysis of nuclear and mitochondrial sequence data

A combined data set of nuclear SSU rDNA, LSU rDNA, and mitochondrial SSU rDNA sequences was analyzed in order to examine the relationships of the major clades of euascomycetes. Partial sequences of 14 ascomycetes were determined and aligned with the corresponding sequences of 16 other ascomycetes retrieved from Genbank. The alignment was analyzed using maximum parsimony (MP) and a Bayesian analysis with Markov chain Monte Carlo (B/MCMC). The classification based on single-gene studies is supported, but the confidence is enhanced in the concatenated analysis. The monophyly of the superclass Leotiomyceta, which includes all euascomycetes with inoperculate asci, is strongly supported. The polyphyly of ascolocularous fungi is supported. The group is divided into two groups: the Dothideomycetes basal to all other Leotiomyceta and the Chaetothyriomycetes as sister-group to Eurotiomycetes. The Lecanoromycetes appear as a monophyletic group with strong support and form a sister-group to the Chaetothyriomycetes/Eurotiomycetes clade, but this lacks support. The Leotiomycetes and Sordariomycetes form a strongly supported sister-group. Alternative topologies are tested using parametric bootstrapping; a basal position of the Eurotiomycetes and Leotiomycetes in the Leotiomyceta cannot be rejected, while such a position can be rejected for Chaetothyriomycetes, Lecanoromycetes and Sordariomycetes. The character evolution with regard to ascoma type, ascus type and ascoma-ontogeny is examined using MP and maximum likelihood (ML). While it appears most likely that the ancestor of the inoperculate ascomycetes had apothecia and an ascohymenial ascoma-ontogeny using MP methods, the ML approach shows that there is some uncertainty at the current state of knowledge. The improvement of confidence of the combined data set in comparison with single-gene studies makes us confident that analyses with additional data sets will further improve the confidence and eventually uncover the branching order of euascomycetes.     

rRNA及MaSp1基因证据不支持圆网蛛类(蜘蛛目:妖面蛛总科+园蛛总科)的单系起源

圆网蛛类(妖面蛛总科 园蛛总科)是否为单系,圆网究竟经历一次进化还是多次进化,这是多年来有争论的、悬而未决的蛛形学难题之一。本文测定了包括妖面蛛总科、园蛛总科和非圆网蛛类等类群在内的9科10种蜘蛛线粒体12S rDNA、16S rDNA及核18S rDNA、28S rDNA等4个基因片段序列,并基于4个基因序列的整合数据,分别通过邻接(NJ)法、最大简约(MP)法、最大似然(ML)法和贝叶斯法(Bayesian)分析,对园蛛总科和妖面蛛总科蜘蛛之间的分子系统关系进行了探讨。系统发生结果表明:1)园蛛总科和妖面蛛总科蜘蛛不是姊妹群,从而支持这两个类群的圆网是平行演化而非同源演化的观点;2)筛器类蜘蛛并非单系发生而为多系发生。另外,依据编码大壶状腺丝蛋白-1(MaSp1)C末端非重复氨基酸序列区段的核酸序列重建的系统发生树也证实圆网蛛类并非单系发生。     

Inverse dispersal patterns in a group of ant parasitoids (Hymenoptera: Eucharitidae: Oraseminae) and their ant hosts

When postulating evolutionary hypotheses for diverse groups of taxa using molecular data, there is a tradeoff between sampling large numbers of taxa with a few Sanger-sequenced genes or sampling fewer taxa with hundreds to thousands of next-generation-sequenced genes. High taxon sampling enables the testing of evolutionary hypotheses that are sensitive to sampling bias (i.e. dating, biogeography and diversification analyses), whereas high character sampling improves resolution of critical nodes. In a group of ant parasitoids (Hymenoptera: Eucharitidae: Oraseminae), we analyse both of these types of datasets independently (203 taxa with five Sanger loci, 92 taxa with 348 anchored hybrid enrichment loci) and in combination (229 taxa, 353 loci) to explore divergence dating, biogeography, host relationships and differential rates of diversification. Oraseminae specialize as parasitoids of the immature stages of ants in the subfamily Myrmicinae (Hymenoptera: Formicidae), with ants in the genus Pheidole being their most common and presumed ancestral host. A general assumption is that the distribution of the parasite must be limited by any range contraction or expansion of its host. Recent studies support a single New World to Old World dispersal pattern for Pheidole at c. 11–27 Ma. Using multiple phylogenetic inference methods (parsimony, maximum likelihood, dated Bayesian and coalescent analyses), we provide a robust phylogeny showing that Oraseminae dispersed in the opposite direction, from Old World to New World, c. 24–33 Ma, which implies that they existed in the Old World before their presumed ancestral hosts. Their dispersal into the New World appears to have promoted an increased diversification rate. Both the host and parasitoid show single unidirectional dispersals in accordance with the presence of the Beringian Land Bridge during the Oligocene, a time when the changing northern climate probably limited the dispersal ability of such tropically adapted groups.     

从核rDNA序列探讨隙蛛亚科Coelotinae的分类地位

隙蛛亚科Coelotinae主要分布于东亚地区,其中我国的已有种类占到全世界种数的一半以上,因此对于我国隙蛛类蜘蛛的研究已经成为世界暗蛛科研究的重点之一。隙蛛亚科属于无筛器类群,于1893年,由Cambridge以隙蛛属为模式属而建立,归属于无筛器的漏斗蛛科。之后,虽然经历了数次修订     

Phylogenetic relationships of leptophlebiid mayflies as inferred by histone H3 and 28S ribosomal DNA

Abstract Leptophlebiidae is among the largest and most diverse groups of extant mayflies (Ephemeroptera), but little is known of family‐level phylogenetic relationships. Using two nuclear genes (the D2 + D3 region of 28S ribosomal DNA and histone H3) and maximum parsimony (MP), maximum likelihood (ML) and Bayesian inference (BI), we inferred the evolutionary relationships of 69 leptophlebiids sampled from six continents and representing 30 genera plus 11 taxa of uncertain taxonomic rank from Madagascar and Papua New Guinea. Although we did not recover monophyly of the Leptophlebiidae, monophyly of two of the three leptophlebiid subfamilies, Habrophlebiinae and Leptophlebiinae, was recovered with moderate to strong support in most analyses. The Atalophlebiinae was rendered paraphyletic as a result of the inclusion of members of Ephemerellidae or the Leptophlebiinae clade. For the species‐rich Atalophlebiinae, four groups of taxa were recovered with moderate to strong branch support: (i) an endemic Malagasy clade, (ii) a Paleoaustral group, a pan‐continental cluster with members drawn from across the southern hemisphere, (iii) a group, uniting fauna from North America, southeast Asia and Madagascar, which we call the Choroterpes group and (iv) a group uniting three New World genera, Thraulodes, Farrodes and Traverella. Knowledge of the phylogenetic relationships of the leptophlebiids will aid in future studies of morphological evolution and biogeographical patterns in this highly diverse and speciose family of mayflies.     

Phylogeny of seed plants based on evidence from eight genes

Relationships among the five groups of extant seed plants (cycads, Ginkgo, conifers, Gnetales, and angiosperms) remain uncertain. To explore relationships among groups of extant seed plants further and to attempt to explain the conflict among molecular data sets, we assembled a data set of four plastid (cpDNA) genes (rbcL, atpB, psaA, and psbB), three mitochondrial (mtDNA) genes (mtSSU, coxI, and atpA), and one nuclear gene (18S rDNA) for 19 exemplars representing the five groups of living seed plants. Analyses of the combined eight-gene data set (15?772 base pairs/taxon) with maximum parsimony (MP), maximum likelihood (ML), and Bayesian approaches reveal a gymnosperm clade that is sister to angiosperms. Within the gymnosperms, a conifer clade includes Gnetales as sister to Pinaceae. Cycads and Ginkgo are either successive sisters to this conifer clade (including Gnetales) or a clade that is sister to conifers and Gnetales. All analyses of the mtDNA partition and ML analyses of the nuclear partition yield very similar topologies. However, MP analyses of the combined cpDNA genes place Gnetales as sister to all other seed plants with strong bootstrap support, whereas ML and Bayesian analyses of the cpDNA data set place Gnetales as sister to Pinaceae. Maximum parsimony and ML analyses of first and second codon positions of the cpDNA partiation also place Gnetales as sister to Pinaceae. In contrast, MP analyses of third codon positions place Gnetales as sister to other seed plants, although ML analyses of third codon positions place Gnetales with Pinaceae. Thus, most of the discrepancies in seed plant topologies involve third codon positions of cpDNA genes. The likelihood ratio (LR) and Shimodaira-Hasegasa (SH) tests were applied to the cpDNA data. The preferred topology based on the LR test is that Gnetales are sister to Pseudotsuga. The SH test based on first and second codon and all three codon positions indicated that there is no significant difference between the best topology (Gnetales sister to Pseudotsuga) and Gnetales sister to a conifer clade. However, there is a significant difference between the best topology and topologies in which Gnetales are sister to the rest of the seed plants or Gnetales sister to angiosperms.