Molecular data from the chloroplast rpoC1 gene suggest a deep and distinct dichotomy of contemporary spermatophytes into two monophyla: gymnosperms (including Gnetales) and angiosperms

Development and evolution of the female gametophyte and fertilization process in Welwitschia mirabilis (Welwitschiaceae)

PREMISE OF THE STUDY

The female gametophyte of Welwitschia has long been viewed as highly divergent from other members of the Gnetales and, indeed, all other seed plants. However, the formation of female gametes and the process of fertilization have never been observed.

METHODS

Standard histological techniques were applied to study gametophyte development and the fertilization process in Welwitschia.

KEY RESULTS

In Welwitschia, fertilization events occur when pollen tubes with binucleate sperm cells grow down through the nucellus and encounter prothallial tubes, free nuclear tubular extensions of the micropylar end of the female gametophyte that grow up through the nucellus. Entry of a binucleate sperm cell into a vacuolate prothallial tube appears to stimulate the rapid coagulation of cytoplasm around a single female nucleus, which differentiates into an egg cell. One sperm nucleus enters the female gamete, while the second sperm nucleus remains outside and ultimately degenerates. Only a single fertilization event occurs per mating pair of pollen tube and prothallial tube.

CONCLUSIONS

Welwitschia lacks the gnetalean pattern of regular double fertilization, as found in Ephedra and Gnetum, involving sperm from a single pollen tube to yield two zygotes. Moreover, an analysis of character evolution indicates that the female gametophyte of Welwitschia is highly apomorphic both among seed plants, and specifically within Gnetales, but also shares several key synapomorphies with its sister taxon Gnetum. Finally, the biological role of prothallial tubes in Welwitschia is examined from the perspectives of gamete competition and kin conflict.

Genetic diversity, genetic structure and demographic history of Cycas simplicipinna (Cycadaceae) assessed by DNA sequences and SSR markers

BACKGROUND

Cycas simplicipinna (T. Smitinand) K. Hill. (Cycadaceae) is an endangered species in China. There were seven populations and 118 individuals that we could collect were genotyped in this study. Here, we assessed the genetic diversity, genetic structure and demographic history of this species.

RESULTS

Analyses of data of DNA sequences (two maternally inherited intergenic spacers of chloroplast, cpDNA and one biparentally inherited internal transcribed spacer region ITS4-ITS5, nrDNA) and sixteen microsatellite loci (SSR) were conducted in the species. Of the 118 samples, 86 individuals from the seven populations were used for DNA sequencing and 115 individuals from six populations were used for the microsatellite study. We found high genetic diversity at the species level, low genetic diversity within each of the seven populations and high genetic differentiation among the populations. There was a clear genetic structure within populations of C. simplicipinna. A demographic history inferred from DNA sequencing data indicates that C. simplicipinna experienced a recent population contraction without retreating to a common refugium during the last glacial period. The results derived from SSR data also showed that C. simplicipinna underwent past effective population contraction, likely during the Pleistocene.

CONCLUSIONS

Some genetic features of C. simplicipinna such as having high genetic differentiation among the populations, a clear genetic structure and a recent population contraction could provide guidelines for protecting this endangered species from extinction. Furthermore, the genetic features with population dynamics of the species in our study would help provide insights and guidelines for protecting other endangered species effectively.

A comparison of 16 DNA regions for use as phylogenetic markers in the pleurocarpous moss genus Plagiothecium (Hypnales)

PREMISE OF THE STUDY

Within the Hypnales-the most derived and species-rich order of pleurocarpous mosses - phylogenies at or below the family level often show poor resolution. In preparation for a phylogeny of the genus Plagiothecium, we wished to identify the DNA markers best suited for evolutionary reconstruction in this group of hypnalean pleurocarps.

METHODS

For each of 25 collections of Plagiothecium and associated taxa, 16 DNA regions were sequenced: nuclear ITS and 26S, and plastid rps4, rps4-trnL, trnL-F, trnK (matK)-psbA, psbA-trnH, trnM-V, trnD-T, rbcL, atpB-rbcL, psbT-H, rpoC1 exon 2 (partial), the trnG intron, the rpl16 intron and the plastid ribosomal spacer DNA (cpITS). Each region was evaluated on the basis of its ability to resolve clades, the amount of homoplasy present in the data set, and the relative ease of obtaining the data. Descriptive statistics for each region are given.

KEY RESULTS

Under-utilized plastid markers for bryophytes such as trnK-psbA, rps4-trnL, and trnD-T outperformed more traditional markers such as trnL-F and rps4. Individual plastid topologies were similar, suggesting that only a limited amount of plastid data are needed to recover a backbone phylogeny. Adding a small amount of nuclear ribosomal data to a large plastid matrix restructured the recovered topology, emphasizing the importance of sampling multiple genomes and the need for new low-copy nuclear markers in bryophyte systematics.

CONCLUSIONS

Future genus-level phylogenies of pleurocarpous mosses should target under-utilized plastid markers such as trnK-psbA and rps4-trnL in conjunction with low-copy nuclear markers.

Phylogenomics and coalescent analyses resolve extant seed plant relationships

The extant seed plants include more than 260,000 species that belong to five main lineages: angiosperms, conifers, cycads, Ginkgo, and gnetophytes. Despite tremendous effort using molecular data, phylogenetic relationships among these five lineages remain uncertain. Here, we provide the first broad coalescent-based species tree estimation of seed plants using genome-scale nuclear and plastid data By incorporating 305 nuclear genes and 47 plastid genes from 14 species, we identify that i) extant gymnosperms (i.e., conifers, cycads, Ginkgo, and gnetophytes) are monophyletic, ii) gnetophytes exhibit discordant placements within conifers between their nuclear and plastid genomes, and iii) cycads plus Ginkgo form a clade that is sister to all remaining extant gymnosperms. We additionally observe that the placement of Ginkgo inferred from coalescent analyses is congruent across different nucleotide rate partitions. In contrast, the standard concatenation method produces strongly supported, but incongruent placements of Ginkgo between slow- and fast-evolving sites. Specifically, fast-evolving sites yield relationships in conflict with coalescent analyses. We hypothesize that this incongruence may be related to the way in which concatenation methods treat sites with elevated nucleotide substitution rates. More empirical and simulation investigations are needed to understand this potential weakness of concatenation methods.

A functional phylogenomic view of the seed plants

A novel result of the current research is the development and implementation of a unique functional phylogenomic approach that explores the genomic origins of seed plant diversification. We first use 22,833 sets of orthologs from the nuclear genomes of 101 genera across land plants to reconstruct their phylogenetic relationships. One of the more salient results is the resolution of some enigmatic relationships in seed plant phylogeny, such as the placement of Gnetales as sister to the rest of the gymnosperms. In using this novel phylogenomic approach, we were also able to identify overrepresented functional gene ontology categories in genes that provide positive branch support for major nodes prompting new hypotheses for genes associated with the diversification of angiosperms. For example, RNA interference (RNAi) has played a significant role in the divergence of monocots from other angiosperms, which has experimental support in Arabidopsis and rice. This analysis also implied that the second largest subunit of RNA polymerase IV and V (NRPD2) played a prominent role in the divergence of gymnosperms. This hypothesis is supported by the lack of 24nt siRNA in conifers, the maternal control of small RNA in the seeds of flowering plants, and the emergence of double fertilization in angiosperms. Our approach takes advantage of genomic data to define orthologs, reconstruct relationships, and narrow down candidate genes involved in plant evolution within a phylogenomic view of species' diversification.

The impact of outgroup choice and missing data on major seed plant phylogenetics using genome-wide EST data

BACKGROUND

Genome level analyses have enhanced our view of phylogenetics in many areas of the tree of life. With the production of whole genome DNA sequences of hundreds of organisms and large-scale EST databases a large number of candidate genes for inclusion into phylogenetic analysis have become available. In this work, we exploit the burgeoning genomic data being generated for plant genomes to address one of the more important plant phylogenetic questions concerning the hierarchical relationships of the several major seed plant lineages (angiosperms, Cycadales, Gingkoales, Gnetales, and Coniferales), which continues to be a work in progress, despite numerous studies using single, few or several genes and morphology datasets. Although most recent studies support the notion that gymnosperms and angiosperms are monophyletic and sister groups, they differ on the topological arrangements within each major group.

METHODOLOGY

We exploited the EST database to construct a supermatrix of DNA sequences (over 1,200 concatenated orthologous gene partitions for 17 taxa) to examine non-flowering seed plant relationships. This analysis employed programs that offer rapid and robust orthology determination of novel, short sequences from plant ESTs based on reference seed plant genomes. Our phylogenetic analysis retrieved an unbiased (with respect to gene choice), well-resolved and highly supported phylogenetic hypothesis that was robust to various outgroup combinations.

CONCLUSIONS

We evaluated character support and the relative contribution of numerous variables (e.g. gene number, missing data, partitioning schemes, taxon sampling and outgroup choice) on tree topology, stability and support metrics. Our results indicate that while missing characters and order of addition of genes to an analysis do not influence branch support, inadequate taxon sampling and limited choice of outgroup(s) can lead to spurious inference of phylogeny when dealing with phylogenomic scale data sets. As expected, support and resolution increases significantly as more informative characters are added, until reaching a threshold, beyond which support metrics stabilize, and the effect of adding conflicting characters is minimized.

Phylogeny of flowering plants by the chloroplast genome sequences: in search of a "lucky gene"

One of the most complicated remaining problems of molecular-phylogenetic analysis is choosing an appropriate genome region. In an ideal case, such a region should have two specific properties: (i) results of analysis using this region should be similar to the results of multigene analysis using the maximal number of regions; (ii) this region should be arranged compactly and be significantly shorter than the multigene set. The second condition is necessary to facilitate sequencing and extension of taxons under analysis, the number of which is also crucial for molecular phylogenetic analysis. Such regions have been revealed for some groups of animals and have been designated as "lucky genes". We have carried out a computational experiment on analysis of 41 complete chloroplast genomes of flowering plants aimed at searching for a "lucky gene" for reconstruction of their phylogeny. It is shown that the phylogenetic tree inferred from a combination of translated nucleotide sequences of genes encoding subunits of plastid RNA polymerase is closest to the tree constructed using all protein coding sites of the chloroplast genome. The only node for which a contradiction is observed is unstable according to the different type analyses. For all the other genes or their combinations, the coincidence is significantly worse. The RNA polymerase genes are compactly arranged in the genome and are fourfold shorter than the total length of protein coding genes used for phylogenetic analysis. The combination of all necessary features makes this group of genes main candidates for the role of "lucky gene" in studying phylogeny of flowering plants.

Phylogenetic signal in nucleotide data from seed plants: implications for resolving the seed plant tree of life

Effects of taxonomic sampling and conflicting signal on the inference of seed plant trees supported in previous molecular analyses were explored using 13 single-locus data sets. Changing the number of taxa in single-locus analyses had limited effects on log likelihood differences between the gnepine (Gnetales plus Pinaceae) and gnetifer (Gnetales plus conifers) trees. Distinguishing among these trees also was little affected by the use of different substitution parameters. The 13-locus combined data set was partitioned into nine classes based on substitution rates. Sites evolving at intermediate rates had the best likelihood and parsimony scores on gnepine trees, and those evolving at the fastest rates had the best parsimony scores on Gnetales-sister trees (Gnetales plus other seed plants). When the fastest evolving sites were excluded from parsimony analyses, well-supported gnepine trees were inferred from the combined data and from each genomic partition. When all sites were included, Gnetales-sister trees were inferred from the combined data, whereas a different tree was inferred from each genomic partition. Maximum likelihood trees from the combined data and from each genomic partition were well-supported gnepine trees. A preliminary stratigraphic test highlights the poor fit of Gnetales-sister trees to the fossil data.

The evolution of plant development

The last decade has witnessed a resurgence in the study of the evolution of plant development, combining investigations in systematics, developmental morphology, molecular developmental genetics, and molecular evolution. The integration of phylogenetic studies, structural analyses of fossil and extant taxa, and molecular developmental genetic information allows the formulation of explicit and testable hypotheses for the evolution of morphological characters. These comprehensive approaches provide opportunities to dissect the evolution of major developmental transitions among land plants, including those associated with apical meristems, the origins of the root/shoot dichotomy, diversification of leaves, and origin and subsequent modification of flower structure. The evolution of these major developmental innovations is discussed within both phylogenetic and molecular genetic contexts. We conclude that it is the combination of these approaches that will lead to the greatest understanding of the evolution of plant development.

The sequence of the largest subunit of RNA polymerase II is a useful marker for inferring seed plant phylogeny

We used RT-PCR to sequence approximately 3 kb of the gene coding for the largest subunit of RNA polymerase II (rpb1) from nine land plants. Our results show that plant rpb1 genes all have a similar GC-content and that their amino acid sequences evolve at a similar rate in most species we examined, except for the Arabidopsis thaliana and rice sequences which evolve faster. This gene also exists as a single copy in most species and contains enough phylogenetically informative sites to resolve the evolutionary relationships among seed plants. Protein maximum parsimony, as well as neighbor-joining and maximum likelihood analyses of DNA and protein sequences, all generated identical tree topologies with similar strong support values at each node. The angiosperms are a clade comprising Amborella as a sister group to all other angiosperms, followed by Nymphaea, Magnolia, Arabidopsis, and a monocot clade containing maize and rice. The gymnosperms also form a monophyletic clade with Welwitschia and pine grouped together and sister to a Cycas and Zamia clade. These findings concur with recent studies that refute the Anthophyte Hypothesis and place Amborella at the base of the angiosperm tree. These rpb1 sequences also give a more consistent picture of seed plant relationships than similar analyses performed on data sets made of 18S rDNA, atpB, and rbcL sequences from the same species. These sequences therefore show great promise to help further resolve the phylogenetic relationships of seed plants.

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 evolution of embryo size in angiosperms and other seed plants: implications for the evolution of seed dormancy

Seed dormancy plays an important role in germination ecology and seed plant evolution. Morphological seed dormancy is caused by an underdeveloped embryo that must mature prior to germination. It has been suggested that the presence of an underdeveloped embryo is plesiomorphic among seed plants and that parallel directional change in embryo morphology has occurred separately in gymnosperms and in angiosperms. We test these hypotheses using original data on embryo morphology of key basal taxa, a published dataset, and the generalized least squares (GLS) method of ancestral character state reconstruction. Reconstructions for embryo to seed ratio (E:S) using family means for 179 families showed that E:S has increased between the ancestral angiosperm and almost all extant angiosperm taxa. Species in the rosid clade have particularly large embryos relative to the angiosperm ancestor. Results for the gymnosperms show a similar but smaller increase. There were no statistically significant differences in E:S between basal taxa and any derived group due to extremely large standard errors produced by GLS models. However, differences between reconstructed values for the angiosperm ancestor and more highly nested nodes are large and these results are robust to topological and branch-length manipulations. Our analysis supports the idea that the underdeveloped embryo is primitive among seed plants and that there has been a directional change in E:S within both angiosperms and gymnosperms. Our analysis suggests that dormancy enforced by an underdeveloped embryo is plesiomorphic among angiosperms and that nondormancy and other dormancy types probably evolved within the angiosperms. The shift in E:S was likely a heterochronic change, and has important implications for the life history of seed plants.

On the origin of class B floral homeotic genes: functional substitution and dominant inhibition in Arabidopsis by expression of an orthologue from the gymnosperm Gnetum

Class B floral homeotic genes are involved in specifying stamen and petal identity in angiosperms (flowering plants). Here we report that gymnosperms, the closest relatives of the angiosperms, contain at least two different clades representing putative orthologues of class B genes, termed GGM2-like and DAL12-like genes. To obtain information about the functional conservation of the class B genes in seed plants, the representative of one of these clades from Gnetum, termed GGM2, was expressed under the control of the CaMV 35S promoter in Arabidopsis wild-type plants and in different class B mutants. In wild-type plants and in a conditional mutant grown at a permissive temperature, gain-of-function phenotypes were obtained in whorls 1 and 4, where class B genes are usually not expressed. In contrast, loss-of-function phenotypes were observed in whorls 2 and 3, where class B genes are expressed. In different class B gene null mutants of Arabidopsis, and in the conditional B mutant grown at the non-permissive temperature, a partial complementation of the mutant phenotype was obtained. In situ hybridization studies and class B gene promoter test fusion experiments demonstrated that the gain-of-function phenotypes are not due to an upregulation of the endogenous B genes from Arabidopsis, and hence probably involve interactions between GGM2 protein homodimers and class B protein target genes other than the Arabidopsis class B genes itself. To our knowledge, this is the first time that partial complementation of a homeotic mutant by an orthologous gene from a distantly related species has been reported. These data suggest that GGM2 has a function in the gymnosperm Gnetum which is related to that of class B floral organ identity genes of angiosperms. That function may be in the specification of male reproductive organ identity, and in distinguishing male from female reproductive organs.

Anatomy of umkomasia (corystospermales) from the triassic of antarctica

The permineralized, corystosperm, cupulate, ovule-bearing organ Umkomasia resinosa is described from the early Middle Triassic of Antarctica. This is the first description of anatomically preserved Umkomasia, which consists of a determinate cupulate branch with helically arranged, recurved, pedicellate megasporophylls, each of which bears one or two abaxially attached unitegmic ovules. Cupules are ovoid, bilobed with elongate ventral and dorsal openings or unlobed with a single ventral opening, and have a two-zoned parenchymatous cortex and abundant secretory cavities. Ovules are small, orthotropous, and possess a thin integument that contains numerous secretory cavities. The ovules are broadly attached at the base, with a bifid integumentary apex that extends past the cupule lobes. The cupulate branch displays stem-like anatomy, producing paired traces into each cupule stalk. These structurally preserved ovulate organs can be related to other corystosperm organs from Antarctica, particularly the pollen-organ Pteruchus fremouwensis. Both anatomical and morphological features support interpretation of corystosperm reproductive structures as branching systems rather than as compound sporophylls. As a result of an increased understanding of the organization of Umkomasia, it appears doubtful that any direct relationship exists between the corystosperm and angiosperm lineages.

The evolutionary split of Pinaceae from other conifers: evidence from an intron loss and a multigene phylogeny

The second intron in the mitochondrial gene nad1 was surveyed using PCR, DNA sequencing, or Southern hybridization in 323 species (313 genera, 212 families) of seed plants. The intron was absent in all 22 species (22 genera, 8 families) of non-Pinaceae conifers studied, in Welwitschia mirabilis, and in seven angiosperms. Whereas absence of the intron in seven angiosperms and Welwitschia is likely due to seven independent losses when evaluated against the recently published multigene phylogenies, the lack of the intron in all non-Pinaceae conifers can be best explained by a single loss. These data suggest that the non-Pinaceae conifers represent a monophyletic group. We also conducted a phylogenetic analysis of seed plants using a combined data set of the partial exon and intron sequences of nad1 generated from this study and published sequences of mitochondrial cox1 and small subunit (SSU) rDNA, chloroplast rbcL, and nuclear 18S rDNA. The results supported the split of conifers into two groups: Pinaceae and non-Pinaceae conifers. The Gnetales were sister to Pinaceae, in agreement with the conclusion from other recent molecular phylogenetic studies that refute the anthophyte hypothesis.

Developmental and evolutionary hypotheses for the origin of double fertilization and endosperm

The discovery of a second fertilization event that initiates endosperm in flowering plants, just over a century ago, stimulated intense interest in the evolutionary history and homology of endosperm, the genetically biparental embryo-nourishing tissue that is found only in angiosperms. Two alternative hypotheses for the origin of double fertilization and endosperm have been invoked to explain the origin of the angiosperm reproductive syndrome from a typical non-flowering seed plant reproductive syndrome. Endosperm may have arisen from a developmental transformation of a supernumerary embryo derived from a rudimentary second fertilization event that first evolved in the ancestors of angiosperms (endosperm homologous with an embryo). Conversely, endosperm may represent the developmental transformation of the cellular phase of non-flowering seed plant female gametophyte ontogeny that was later sexualized by the addition of a second fertilization event in a strongly progenetic female gametophyte (endosperm homologous with a female gametophyte). For the first time, explicit developmental and evolutionary transitions for both of these hypotheses are examined and compared. In addition, current data that may be congruent with either of these hypotheses are discussed. It is clear that much remains to be accomplished if the evolutionary significance of the process of double fertilization and the formation of endosperm is to be fully understood.

Perspective: the origin of flowering plants and their reproductive biology--a tale of two phylogenies

Recently, two areas of plant phylogeny have developed in ways that could not have been anticipated, even a few years ago. Among extant seed plants, new phylogenetic hypotheses suggest that Gnetales, a group of nonflowering seed plants widely hypothesized to be the closest extant relatives of angiosperms, may be less closely related to angiosperms than was believed. In addition, recent phylogenetic analyses of angiosperms have, for the first time, clearly identified the earliest lineages of flowering plants: Amborella, Nymphaeales, and a clade that includes Illiciales/ Trimeniaceae/Austrobaileyaceae. Together, the new seed plant and angiosperm phylogenetic hypotheses have major implications for interpretation of homology and character evolution associated with the origin and early history of flowering plants. As an example of the complex and often unpredictable interplay of phylogenetic and comparative biology, we analyze the evolution of double fertilization, a process that forms a diploid embryo and a triploid endosperm, the embryo-nourishing tissue unique to flowering plants. We demonstrate how the new phylogenetic hypotheses for seed plants and angiosperms can significantly alter previous interpretations of evolutionary homology and firmly entrenched assumptions about what is synapomorphic of flowering plants. In the case of endosperm, a solution to the century-old question of its potential homology with an embryo or a female gametophyte (the haploid egg-producing generation within the life cycle of a seed plant) remains complex and elusive. Too little is known of the comparative reproductive biology of extant nonflowering seed plants (Gnetales, conifers, cycads, and Ginkgo) to analyze definitively the potential homology of endosperm with antecedent structures. Remarkably, the new angiosperm phylogenies reveal that a second fertilization event to yield a biparental endosperm, long assumed to be an important synapomorphy of flowering plants, cannot be conclusively resolved as ancestral for flowering plants. Although substantive progress has been made in the analysis of phylogenetic relationships of seed plants and angiosperms, these efforts have not been matched by comparable levels of activity in comparative biology. The consequence of inadequate comparative biological information in an age of phylogenetic biology is a severe limitation on the potential to reconstruct key evolutionary historical events.

Ferulic acid is bound to the primary cell walls of all gymnosperm families

Unlignified primary cell walls containing ester-linked ferulic acid fluoresce blue in ultraviolet radiation which changes to green with increased intensity on treatment with ammonium hydroxide. Using this fluorescence behaviour, we detected ester-linked ferulic acid in the primary cell walls of all 41 species of gymnosperms we examined. These species were in 17 families representing all four extant classes of gymnosperms. In addition, we obtained cell-wall preparations containing >95% primary cell walls from nine gymnosperm species in nine families, representing all four extant classes. These preparations were analysed for ester-linked monomeric phenolic acids. We found ferulic acid (mostly trans) (88-1,561µg/g cell walls) in all of the preparations and p-coumaric acid (mostly trans) (0-106µg/g cell walls) in all except one of them. Ferulic acid ester-linked to primary cell walls has previously been found in angiosperms: in the commelinoid monocotyledons and in the dicotyledon order Caryophyllales, both monophyletic groups. From the present results, we postulate that the extant classes of gymnosperms are monophyletic and no class is sister to the angiosperms.

MADS-Box gene diversity in seed plants 300 million years ago

MADS-box genes encode a family of transcription factors which control diverse developmental processes in flowering plants ranging from root development to flower and fruit development. Through phylogeny reconstructions, most of these genes can be subdivided into defined monophyletic gene clades whose members share similar expression patterns and functions. Therefore, the establishment of the diversity of gene clades was probably an important event in land plant evolution. In order to determine when these clades originated, we isolated cDNAs of 19 different MADS-box genes from Gnetum gnemon, a gymnosperm model species and thus a representative of the sister group of the angiosperms. Phylogeny reconstructions involving all published MADS-box genes were then used to identify gene clades containing putative orthologs from both angiosperm and gymnosperm lineages. Thus, the minimal number of MADS-box genes that were already present in the last common ancestor of extant gymnosperms and angiosperms was determined. Comparative expression studies involving pairs of putatively orthologous genes revealed a diversity of patterns that has been largely conserved since the time when the angiosperm and gymnosperm lineages separated. Taken together, our data suggest that there were already at least seven different MADS-box genes present at the base of extant seed plants about 300 MYA. These genes were probably already quite diverse in terms of both sequence and function. In addition, our data demonstrate that the MADS-box gene families of extant gymnosperms and angiosperms are of similar complexities.

Phylogeny of seed plants based on all three genomic compartments: extant gymnosperms are monophyletic and Gnetales' closest relatives are conifers

Efforts to resolve Darwin's "abominable mystery"-the origin of angiosperms-have led to the conclusion that Gnetales and various fossil groups are sister to angiosperms, forming the "anthophytes." Morphological homologies, however, are difficult to interpret, and molecular data have not provided clear resolution of relationships among major groups of seed plants. We introduce two sequence data sets from slowly evolving mitochondrial genes, cox1 and atpA, which unambiguously reject the anthophyte hypothesis, favoring instead a close relationship between Gnetales and conifers. Parsimony- and likelihood-based analyses of plastid rbcL and nuclear 18S rDNA alone and with cox1 and atpA also strongly support a gnetophyte-conifer grouping. Surprisingly, three of four genes (all but nuclear rDNA) and combined three-genome analyses also suggest or strongly support Gnetales as derived conifers, sister to Pinaceae. Analyses with outgroups screened to avoid long branches consistently identify all gymnosperms as a monophyletic sister group to angiosperms. Combined three- and four-gene rooted analyses resolve the branching order for the remaining major groups-cycads separate from other gymnosperms first, followed by Ginkgo and then (Gnetales + Pinaceae) sister to a monophyletic group with all other conifer families. The molecular phylogeny strongly conflicts with current interpretations of seed plant morphology, and implies that many similarities between gnetophytes and angiosperms, such as "flower-like" reproductive structures and double fertilization, were independently derived, whereas other characters could emerge as synapomorphies for an expanded conifer group including Gnetales. An initial angiosperm-gymnosperm split implies a long stem lineage preceding the explosive Mesozoic radiation of flowering plants and suggests that angiosperm origins and homologies should be sought among extinct seed plant groups.

Seed plant phylogeny inferred from all three plant genomes: monophyly of extant gymnosperms and origin of Gnetales from conifers

Phylogenetic relationships among the five groups of extant seed plants are presently quite unclear. For example, morphological studies consistently identify the Gnetales as the extant sister group to angiosperms (the so-called "anthophyte" hypothesis), whereas a number of molecular studies recover gymnosperm monophyly, and few agree with the morphology-based placement of Gnetales. To better resolve these and other unsettled issues, we have generated a new molecular data set of mitochondrial small subunit rRNA sequences, and have analyzed these data together with comparable data sets for the nuclear small subunit rRNA gene and the chloroplast rbcL gene. All nuclear analyses strongly ally Gnetales with a monophyletic conifers, whereas all mitochondrial analyses and those chloroplast analyses that take into account saturation of third-codon position transitions actually place Gnetales within conifers, as the sister group to the Pinaceae. Combined analyses of all three genes strongly support this latter relationship, which to our knowledge has never been suggested before. The combined analyses also strongly support monophyly of extant gymnosperms, with cycads identified as the basal-most group of gymnosperms, Ginkgo as the next basal, and all conifers except for Pinaceae as sister to the Gnetales + Pinaceae clade. According to these findings, the Gnetales may be viewed as extremely divergent conifers, and the many morphological similarities between angiosperms and Gnetales (e.g., double fertilization and flower-like reproductive structures) arose independently.

Seed plant phylogeny: Demise of the anthophyte hypothesis?

Recent molecular phylogenetic studies indicate, surprisingly, that Gnetales are related to conifers, or even derived from them, and that no other extant seed plants are closely related to angiosperms. Are these results believable? Is this a clash between molecules and morphology?