Angiosperm phylogeny based on matK sequence information

Plastid matK gene sequences for 374 genera representing all angiosperm orders and 12 genera of gymnosperms were analyzed using parsimony (MP) and Bayesian inference (BI) approaches. Traditionally, slowly evolving genomic regions have been preferred for deep-level phylogenetic inference in angiosperms. The matK gene evolves approximately three times faster than the widely used plastid genes rbcL and atpB. The MP and BI trees are highly congruent. The robustness of the strict consensus tree supercedes all individual gene analyses and is comparable only to multigene-based phylogenies. Of the 385 nodes resolved, 79% are supported by high jackknife values, averaging 88%. Amborella is sister to the remaining angiosperms, followed by a grade of Nymphaeaceae and Austrobaileyales. Bayesian inference resolves Amborella + Nymphaeaceae as sister to the rest, but with weak (0.42) posterior probability. The MP analysis shows a trichotomy sister to the Austrobaileyales representing eumagnoliids, monocots + Chloranthales, and Ceratophyllum + eudicots. The matK gene produces the highest internal support yet for basal eudicots and, within core eudicots, resolves a crown group comprising Berberidopsidaceae/Aextoxicaceae, Santalales, and Caryophyllales + asterids. Moreover, matK sequences provide good resolution within many angiosperm orders. Combined analyses of matK and other rapidly evolving DNA regions with available multigene data sets have strong potential to enhance resolution and internal support in deep level angiosperm phylogenetics and provide additional insights into angiosperm evolution.     

Phylogenetics of early branching eudicots:Comparing phylogenetic signal across plastid introns,spacers,and genes

<正>Recent phylogenetic analyses revealed a grade with Ranunculales,Sabiales,Proteales,Trochodendrales, and Buxales as first branching eudicots,with the respective positions of Proteales and Sabiales still lacking statistical confidence.As previous analyses of conserved plastid genes remain inconclusive,we aimed to use and evaluate a representative set of plastid introns(group I:trnL;group II:petD,rpll6,trnK) and intergenic spacers(trnL-F,petB-petD, atpB-rbcL,rps3-rpl16) in comparison to the rapidly evolving matK and slowly evolving atpB and rbcL genes. Overall patterns of micro structural mutations converged across genomic regions,underscoring the existence of a general mutational pattern throughout the plastid genome.Phylogenetic signal differed strongly between functionally and structurally different genomic regions and was highest in matK,followed by spacers,then group II and group I introns.The more conserved atpB and rbcL coding regions showed distinctly lower phylogenetic information content.Parsimony,maximum likelihood,and Bayesian phylogenetic analyses based on the combined dataset of non-coding and rapidly evolving regions(14 000 aligned characters) converged to a backbone topology of eudicots with Ranunculales branching first,a Proteales-Sabiales clade second,followed by Trochodendrales and Buxales. Gunnerales generally appeared as sister to all remaining core eudicots with maximum support.Our results show that a small number of intron and spacer sequences allow similar insights into phylogenetic relationships of eudicots compared to datasets of many combined genes.The non-coding proportion of the plastid genome thus can be considered an important information source for plastid phylogenomics.     

Phylogenetic systematics and character evolution in the angiosperm family Haloragaceae

The poorly known Haloragaceae R. Br. (Saxifragales) are highly diverse in habit (small trees to submerged aquatics) and labile in floral merosity (2-4), both uncommon among the core eudicots. This family has a cosmopolitan distribution, but taxonomic diversity is concentrated in Australia. An explicit phylogenetic approach has not previously been utilized to examine relationships or character evolution in this family. We used molecular evidence from nrDNA ITS and cpDNA trnK and matK regions under both Bayesian and parsimony analyses to address phylogenetic relationships. Combined molecular analyses defined a monophyletic Haloragaceae with the woody genera (Haloragodendron, Glischrocaryon) sister to the rest. Relationships among many genera were well resolved, with genera as currently delimited generally well supported, although there were notable exceptions; a new genus (Trihaloragis) is recognized, and the aquatic genus Meionectes is again distinct from Haloragis. Three new species combinations are also recognized. There are multiple (two or three) origins of the submerged aquatic habit in the family and potentially an intermediate reversal to the terrestrial habit, neither previously demonstrated in a core eudicot family using an explicit phylogenetic hypothesis. Ancestral character analyses suggest two origins of trimerous flowers and multiple reductions to dimerous flowers throughout Haloragaceae.     

Origin and evolution of the chloroplast trnK (matK) intron: a model for evolution of group II intron RNA structures

The trnK intron of plants encodes the matK open reading frame (ORF), which has been used extensively as a phylogenetic marker for classification of plants. Here we examined the evolution of the trnK intron itself as a model for group II intron evolution in plants. Representative trnK intron sequences were compiled from species spanning algae to angiosperms, and four introns were newly sequenced. Phylogenetic analyses showed that the matK ORFs belong to the ML (mitochondrial-like) subclass of group II intron ORFs, indicating that they were derived from a mobile group II intron of the class. RNA structures of the introns were folded and analyzed, which revealed progressive RNA structural deviations and degenerations throughout plant evolution. The data support a model in which plant organellar group II introns were derived from bacterial-like introns that had "standard" RNA structures and were competent for self-splicing and mobility and that subsequently the ribozyme structures degenerated to ultimately become dependent upon host-splicing factors. We propose that the patterns of RNA structure evolution seen for the trnK intron will apply to the other group II introns in plants.     

Occurrence of matK in a trnK group II intron in charophyte green algae and phylogeny of the Characeae

A group II intron containing the matK gene, which encodes a splicing-associated maturase, was found in the trnK (lysine tRNA) exon in the chloroplast genome of the six extant genera of green algae in the family Characeae, which among green algae are the sister group to embryophytes (land plants). The characean trnK intron (～2.5 kilobases [kb]) and matK ORF (～1.5 kb) are comparable in size to the intron and ORF of land plants, in which they are similarly found inserted in the trnK exon. Domain X, a sequence of conserved amino acid residues within matK, occurs in the Characeae. Phylogenetic analysis using maximum likelihood (GTR + I + gamma likelihood model) and parsimony (branch and bound search) yielded one tree with high bootstrap support for all branches. The matK tree was congruent with the rbcL tree for the same taxa. The number and proportion of informative sites was higher in matK (501, 31% of matK sequence) compared to rbcL (122, 10%). Characeae branch lengths were on average more than five times longer for matK compared to rbcL and provided better resolution within the Characeae. These findings along with recent genomic analyses demonstrate that the intron and matK invaded the chloroplast genome of green algae prior to the evolution of land plants.     

Group II introns as phylogenetic tools: structure, function, and evolutionary constraints

Group II introns comprise the majority of noncoding DNA in many plant chloroplast genomes and include the commonly sequenced regions trnK/matK, the rps16 intron, and the rpl16 intron. As demand increases for nucleotide characters at lower taxonomic levels, chloroplast introns may come to provide the bulk of plastome sequence data for assessment of evolutionary relationships in infrageneric, intergeneric, and interfamilial studies. Group II introns have many attractive properties for the molecular systematist: they are confined to organellar genomes in eukaryotes and the majority are single-copy; they share a well-defined and empirically tested secondary and tertiary structure; and many are easily amplified due to highly conserved sequence in flanking exons. However, structure-linked mutation patterns in group II intron sequences are more complex than generally supposed and have important implications for aligning nucleotides, assessing mutational biases in the data, and selecting appropriate models of character evolution for phylogenetic analysis. This paper presents a summary of group II intron function and structure, reviews the link between that structure and specific mutational constraints in group II intron sequences, and discusses strategies for accommodating the resulting complex mutational patterns in subsequent phylogenetic analyses.     

Angiosperm phylogeny inferred from sequences of four mitochondrial genes

An angiosperm phylogeny was reconstructed in a maximum likelihood analysis of sequences of four mitochondrial genes, atpl, matR, had5, and rps3, from 380 species that represent 376 genera and 296 families of seed plants. It is largely congruent with the phylogeny of angiosperms reconstructed from chloroplast genes atpB, matK, and rbcL, and nuclear 18S rDNA. The basalmost lineage consists of Amborella and Nymphaeales (including Hydatellaceae). Austrobaileyales follow this clade and are sister to the mesangiosperms, which include Chloranthaceae, Ceratophyllum, magnoliids, monocots, and eudicots. With the exception of Chloranthaceae being sister to Ceratophyllum, relationships among these five lineages are not well supported. In eudicots, Ranunculales, Sabiales, Proteales, Trochodendrales, Buxales, Gunnerales, Saxifragales, Vitales, Berberidopsidales, and Dilleniales form a basal grade of lines that diverged before the diversification of rosids and asterids. Within rosids, the COM (Celastrales-Oxalidales-Malpighiales) clade is sister to malvids (or rosid Ⅱ), instead of to the nitrogen-fixing clade as found in all previous large-scale molecular analyses of angiosperms. Santalales and Caryophyllales are members of an expanded asterid clade. This study shows that the mitochondrial genes are informative markers for resolving relationships among genera, families, or higher rank taxa across angiosperms. The low substitution rates and low homoplasy levels of the mitochondrial genes relative to the chloroplast genes, as found in this study, make them particularly useful for reconstructing ancient phylogenetic relationships. A mitochondrial gene-based angiosperm phylogeny provides an independent and essential reference for comparison with hypotheses of angiosperm phylogeny based on chloroplast genes, nuclear genes, and non-molecular data to reconstruct the underlying organismal phylogeny.     

连香树科及其近缘植物matK序列分析和系统学意义

测定和分析了连香树科(Cercidiphyllaeeae)、交让木科(Daplmiphyllaceac)、金缕梅科(Hamamelidaceae)代表植物的叶绿体marK序列(5′端31bps除外),以木兰属作为外类群,应用邻接法构建分子系统树,结果表明：连香树科与水青树科的亲缘关系较远。连香树科、交让木科和金缕梅科形成了一个自展数据支持率(bootstrap)为100％的单系类群,其中金缕梅科枫香属(Liquidambar)、红花荷属(Rhodoleia)和金缕梅属(Hamamelis)虽构成了一个单系类群,但自展数据支持率仅为68％;连香树科与交让木科构成的单系分支自展数据支持率仅为53％。由于连香树科、交让木科、金缕梅科之间的进化距离相当短,表明这3个科之间亲缘关系密切,内部分支的自展数据支持率不高,表明它们之间准确的亲缘关系有待进一步研究。本研究结果与rbcL、aptB、18S rDNA序列分析结果相似,但自展数据支持率更高,表明marK序列分析可应用于较高等级分类群系统发育关系的研究。     

Sequences of 72 plastid genes support the early divergence of Cornales in the asterids

Recent phylogenetic analyses of molecular data have supported different hypotheses of relationships among Cornales,Ericales,and core asterids.Such hypotheses have implications for the evolution of important morphological and embryological features of asterids.In this study we generated plastid genome-scale data of Davidia (Cornales) and Franklinia (Ericales) and combined them with published sequence data of eudicots.Our maximum parsimony,maximum likelihood,and Bayesian analyses generated strongly supported and congruent phylogenetic relationships among the three major lineages of the asterids.Cornales diverges first in asterids; Ericales is sister to the core asterids.Adding two more taxa helps mitigate long branch attraction in parsimony analyses.Sampling 26-28 plastid protein-coding genes may provide satisfactory resolution and support for relationships of eudicots including basal lineages of asterids.     

A Molecular Phylogeny of Costaceae (Zingiberales)

The phylogenetic relationships of Costaceae, a tropical monocotyledonous family sister to the gingers (Zingiberaceae), were investigated with a combination of two chloroplast loci (the trnL-F locus, including the trnL intron, the 3'trnL exon, and the trnL-F intergenic spacer, and the trnK locus, including the trnK intron and the matK coding region) and one nuclear locus (ITS1-5.8s-ITS2). The resulting parsimony analysis of selected taxa that demonstrate the range of floral morphological variation in the family shows that the Cadalvena-type [corrected] floral morphology is ancestral to the group and that both Tapeinochilos species and a Monocostus + Dimerocostus clade represent recent divergences. The genus Costus is broadly paraphyletic but Costus subgenus Eucostus K. Schum. represents a large monophyletic radiation that is poorly resolved. Within this clade, secondary analyses suggest that pollination syndrome, traditionally used for taxonomic and classification purposes within the genus Costus, is a relatively plastic trait of limited phylogenetic utility. This represents the first detailed investigation into intrageneric and interspecific evolutionary relationships within the family Costaceae and presents some novel evolutionary trends with respect to floral morphology and biogeography.     

Molecular phylogenetic evidence for the monophyly of Fritillaria and Lilium (Liliaceae; Liliales) and the infrageneric classification of Fritillaria

We present phylogenetic analyses of 37 taxa of Fritillaria (Liliaceae), 15 species of Lilium, and several outgroup taxa from Liliaceae s.s. to investigate the generic delimitation of Fritillaria in relation to Lilium as well as infrageneric relationships within Fritillaria. We used DNA sequences from the maturase-coding plastid matK gene and the trnK intron, the intron of the ribosomal protein-coding rpl16 plastid gene, and the nuclear ribosomal internal transcribed spacers (ITS). Phylogenetic analysis using maximum parsimony defined Fritillaria and Lilium (the latter including Nomocharis) as sister taxa. Fritillaria sections Fritillaria and Liliorhiza are supported in part, and some of the most enigmatic species usually included in Fritillaria (sections Petilium and Theresia and the monotypic genus Korolkowia) are closely related. The results support the new classification of Fritillaria proposed by Rix. We postulate independent origins of the underground bulbils found in Fritillaria davidii and the remainder of subgenus Liliorhiza.     

Nucleotide sequence variation of the chloroplast trnK/matK region in two wild Fagopyrum (Polygonaceae) species, F. leptopodum and F. statice

Nucleotide sequence polymorphisms of the intron of the chloroplast trnK (UUU) gene, including a matK gene, were investigated within two wild Fagopyrum species, F. leptopodum and F. statice, to assess the degree and pattern of the inter- and intraspecific differences in coding and noncoding chloroplast DNA regions in higher plants. Ten and five accessions were used for F. leptopodum and F. statice, respectively. The length of the trnK intron region in these species ranged from 2494 to 2506 bp. In the trnK intron, the net nucleotide substitution number per site (Da) between the two species was 0.00109, lower than the nucleotide diversity (pi), 0.00195 for F. leptopodum and 0.00144 for F. statice, suggesting a low level of interspecific divergence. This result seems to be due to the phylogenetic pattern that both species are interspersed with each other, which was revealed by the phylogenetic analyses based on the nucleotide substitutions and indels. In the matK gene region (1524 bp), seven and two nucleotide substitutions were found within F. leptopodum and F. statice, respectively. All of the nine nucleotide substitutions (eight of which were nonsynonymous) within and between F. leptopodum and F. statice were clustered in the 5' part of the matK gene region, and no variation was found in the 3' part. This suggests that most of the 3' part is occupied by the conserved domains that are important for the binding activity of the gene product to the precursor mRNA, and therefore implies that the 3' part is more functionally constrained than the 5' part.     

The tortoise and the hare II: relative utility of 21 noncoding chloroplast DNA sequences for phylogenetic analysis

Chloroplast DNA sequences are a primary source of data for plant molecular systematic studies. A few key papers have provided the molecular systematics community with universal primer pairs for noncoding regions that have dominated the field, namely trnL-trnF and trnK/matK. These two regions have provided adequate information to resolve species relationships in some taxa, but often provide little resolution at low taxonomic levels. To obtain better phylogenetic resolution, sequence data from these regions are often coupled with other sequence data. Choosing an appropriate cpDNA region for phylogenetic investigation is difficult because of the scarcity of information about the tempo of evolutionary rates among different noncoding cpDNA regions. The focus of this investigation was to determine whether there is any predictable rate heterogeneity among 21 noncoding cpDNA regions identified as phylogenetically useful at low levels. To test for rate heterogeneity among the different cpDNA regions, we used three species from each of 10 groups representing eight major phylogenetic lineages of phanerogams. The results of this study clearly show that a survey using as few as three representative taxa can be predictive of the amount of phylogenetic information offered by a cpDNA region and that rate heterogeneity exists among noncoding cpDNA regions.     

Revisiting the phylogeny of Ranunculeae: Implications for divergence time estimation and historical biogeography

Previous studies based on different molecular datasets have generated conflicting topologies for Ranunculeae. Here, we revisit the phylogeny of Ranunculeae by analyzing the individual matK/trnK, psbJ-petA, and internal transcribed spacer (ITS) data, the combined matK/trnK, psbJ-petA, and ITS dataset, and the combinedrbcL, trnL-F, matK/trnK, psbJ-petA, and ITS dataset. Based on the tree-based comparisons, blast searches against NCBI of the sequences, and close examination of the alignment, we found that 10 psbJ-petA sequences previously used were questionable (erroneous or problematic) and responsible for previous conflicting topologies. After omitting these questionable sequences, we provide a new phylogeny for Ranunculeae, in which Beckwithia–Cyrtorhyncha, Kumlienia, andPeltocalathos were replaced. These new replacements are supported by corresponding morphological characters. Moreover, three previously proposed intercontinental disjunct distributions within Ranunculus were also refuted. In our framework, our divergence time and biogeographic analyses indicate that divergence time estimates and the ancestral areas reconstructed for 10 of the 15 nodes in the genus-level phylogeny were influenced by elimination of the questionable sequences. The most recent common ancestor of Ranunculeae was inferred to be present in Europe and North America during the late Eocene. Clades I and II began to diversify in Europe and North America, respectively, and subsequently migrated to other continents. This study shows that it is necessary to analyze individual chloroplast DNA region datasets separately to detect questionable sequences early in the study. The combined dataset including the questionable sequences resulted in an erroneous phylogenetic tree, and the use of this tree subsequently affected age estimates and biogeographic analyses.     

The mitochondrial nad2 gene as a novel marker locus for phylogenetic analysis of early land plants: a comparative analysis in mosses

The mitochondrial nad2 gene is established as a novel marker locus for phylogenetic analyses among early land plants. The potential of this gene for phylogenetic resolution was checked with a broad taxon sampling of 42 mosses (Bryopsida, including the enigmatic genus Takakia) to allow both a comparative analysis with the recently explored nad5 gene and the fusion of independent data sets. The mitochondrial gene sequences provide valuable phylogenetic information on the relationships of classically defined orders and their respective monophylies. The more rapidly diverging sequences of a group I intron in nad5 and of a group II intron in nad2 add information for fine resolution. Although both genes provide phylogenetic information in the same taxonomic range (above family level), the combined sequence alignment results in an approximate doubling in the number of nodes with significant bootstrap support (>90). According to our data, Buxbaumiales are a paraphyletic taxon in a key position between the earliest branching taxa (Sphagnales, Takakiales, Andreaeales, Polytrichales, and Tetraphidales) and all other orders, possibly to be placed in the subclass Bryidae. A dichotomy in the latter recalls two previously suggested superorders Hypnanae and Dicrananae. Both genes independently question the monophyly of the orders Dicranales and Neckerales and reject the inclusion of the genera Schistostega, Timmia, and Encalypta among Eubryales.     

Phylogeny of basal eudicots: Insights from non-coding and rapidly evolving DNA

Sequence data of the trnL group I intron, the petD group II intron, the trnL-F and petB-D spacers, and the rapidly evolving matK gene were analysed from all families of the basal eudicot grade and from representatives of 19 core eudicot orders. The dataset comprised 5654 positions of aligned sequence plus a matrix of 1087 binary indel characters. Mutational hotspots correspond in number and extension to hotspots already known from basal angiosperms and, with respect to secondary structure, are generally located in terminal parts of stem-loop regions. Parsimony, Bayesian, and likelihood analyses depict Ranunculales as sister to all remaining eudicots with maximum support. The branching order in the basal eudicot grade is further resolved as Sabiales, Proteales, Trochodendrales, and Buxales. Nearly all of the backbone nodes gain high confidence, except for the node showing Proteales diverging before Trochodendrales, which is only moderately supported (83% JK). In Ranunculales, the woody Eupteleaceae are first-branching, with Papaveraceae plus Fumariaceae coming next. Within Proteales, Nelumbo is clearly resolved as sister to a Platanaceae–Proteaceae clade. Gunnerales are found as the first branch in core eudicots, with maximum support in the combined analysis. This node is also resolved with matK alone, but unsupported. It appears that the combined analysis of sequence data from rapidly evolving and non-coding genomic regions leads to significantly improved statistical support values in comparison to earlier studies of basal eudicots using multiple conserved genes.See also Electronic Supplement at doi:10.1016/j.ode.2006.08.001     

Nuclear and cpDNA sequences combined provide strong inference of higher phylogenetic relationships in the phlox family (Polemoniaceae)

Members of the phlox family (Polemoniaceae) serve as useful models for studying various evolutionary and biological processes. Despite its biological importance, no family-wide phylogenetic estimate based on multiple DNA regions with complete generic sampling is available. Here, we analyze one nuclear and five chloroplast DNA sequence regions (nuclear ITS, chloroplast matK, trnL intron plus trnL-trnF intergeneric spacer, and the trnS-trnG, trnD-trnT, and psbM-trnD intergenic spacers) using parsimony and Bayesian methods, as well as assessments of congruence and long branch attraction, to explore phylogenetic relationships among 84 ingroup species representing all currently recognized Polemoniaceae genera. Relationships inferred from the ITS and concatenated chloroplast regions are similar overall. A combined analysis provides strong support for the monophyly of Polemoniaceae and subfamilies Acanthogilioideae, Cobaeoideae, and Polemonioideae. Relationships among subfamilies, and thus for the precise root of Polemoniaceae, remain poorly supported. Within the largest subfamily, Polemonioideae, four clades corresponding to tribes Polemonieae, Phlocideae, Gilieae, and Loeselieae receive strong support. The monogeneric Polemonieae appears sister to Phlocideae. Relationships within Polemonieae, Phlocideae, and Gilieae are mostly consistent between analyses and data permutations. Many relationships within Loeselieae remain uncertain. Overall, inferred phylogenetic relationships support a higher-level classification for Polemoniaceae proposed in 2000.     

First insights into fern matK phylogeny

MatK, the only maturase gene in the land plant plastid genome, is a very popular phylogenetic marker that has been extensively applied in reconstructing angiosperm phylogeny. However, the use of matK in fern phylogeny is largely unknown, due to difficulties with amplification: ferns have lost the flanking trnK exons, typically the region used for designing stable priming sites. We developed primers that are either universal or lineage-specific that successfully amplify matK across all fern families. To evaluate whether matK is as powerful a phylogenetic marker in ferns as in angiosperms, we compared its sequence characteristics and phylogenetic performance to those of rbcL and atpA. Among these three genes, matK has the highest variability and substitution evenness, yet shows the least homoplasy. Most importantly, applying matK in fern phylogenetics better resolved relationships among families, especially within eupolypods I and II. Here we demonstrate the power of matK for fern phylogenetic reconstruction, as well as provide primers and extensive sequence data that will greatly facilitate future evolutionary studies of ferns.