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

An evolutionary perspective on leaf economics: phylogenetics of leaf mass per area in vascular plants

In plant leaves, resource use follows a trade-off between rapid resource capture and conservative storage. This "worldwide leaf economics spectrum" consists of a suite of intercorrelated leaf traits, among which leaf mass per area, LMA, is one of the most fundamental as it indicates the cost of leaf construction and light-interception borne by plants. We conducted a broad-scale analysis of the evolutionary history of LMA across a large dataset of 5401 vascular plant species. The phylogenetic signal in LMA displayed low but significant conservatism, that is, leaf economics tended to be more similar among close relatives than expected by chance alone. Models of trait evolution indicated that LMA evolved under weak stabilizing selection. Moreover, results suggest that different optimal phenotypes evolved among large clades within which extremes tended to be selected against. Conservatism in LMA was strongly related to growth form, as were selection intensity and phenotypic evolutionary rates: woody plants showed higher conservatism in relation to stronger stabilizing selection and lower evolutionary rates compared to herbaceous taxa. The evolutionary history of LMA thus paints different evolutionary trajectories of vascular plant species across clades, revealing the coordination of leaf trait evolution with growth forms in response to varying selection regimes.

Leaf anatomy and its implications for phylogenetic relationships in Taxaceae s. l

The comparative study on leaf anatomy and stomata structures of six genera of Taxaceae s. l. was conducted. Leaf anatomical structures were very comparable to each other in tissue shape and their arrangements. Taxus, Austrotaxus, and Pseudotaxus have no foliar resin canal, whereas Amentotaxus, Cephalotaxus, and Torreya have a single resin canal located below the vascular bundle. Among them, Torreya was unique with thick-walled, almost round sclerenchymatous epidermal cells. In addition, Amentotaxus and Torreya were comprised of some fiber cells around the vascular bundle. Also, Amentotaxus resembled Cephalotaxus harringtonia and its var. nana because they have discontinuous fibrous hypodermis. However, C. fortunei lacked the same kind of cells. Stomata were arranged in two stomatal bands separated by a mid-vein. The most unique stomatal structure was of Taxus with papillose accessory cells forming stomatal apparatus and of Torreya with deeply seated stomata covered with a special filament structure. Some morphological and molecular studies have already been discussed for the alternative classification of taxad genera into different minor families. The present study is also similar to these hypotheses because each genus has their own individuality in anatomical structure and stomata morphology. In conclusion, these differences in leaf and stomata morphology neither strongly support the two tribes in Taxaceae nor fairly recognize the monogeneric family, Cephalotaxaceae. Rather, it might support an alternative classification of taxad genera in different minor families or a single family Taxaceae including Cephalotaxus. In this study, we would prefer the latter one because there is no clear reason to separate Cephalotaxus from the rest genera of Taxaceae. Therefore, Taxaceae should be redefined with broad circumscriptions including Cephalotaxus.

Highly rearranged and size-variable chloroplast genomes in conifers II clade (cupressophytes): evolution towards shorter intergenic spacers

Although conifers are of immense ecological and economic value, bioengineering of their chloroplasts remains undeveloped. Understanding the chloroplast genomic organization of conifers can facilitate their bioengineering. Members of the conifer II clade (or cupressophytes) are highly diverse in both morphologic features and chloroplast genomic organization. We compared six cupressophyte chloroplast genomes (cpDNAs) that represent four of the five cupressophyte families, including three genomes that are first reported here (Agathis dammara, Calocedrus formosana and Nageia nagi). The six cupressophyte cpDNAs have lost a pair of large inverted repeats (IRs) and vary greatly in size, organization and tRNA copies. We demonstrate that cupressophyte cpDNAs have evolved towards reduced size, largely due to shrunken intergenic spacers. In cupressophytes, cpDNA rearrangements are capable of extending intergenic spacers, and synonymous mutations are negatively associated with the size and frequency of rearrangements. The variable cpDNA sizes of cupressophytes may have been shaped by mutational burden and genomic rearrangements. On the basis of cpDNA organization, our analyses revealed that in gymnosperms, cpDNA rearrangements are phylogenetically informative, which supports the 'gnepines' clade. In addition, removal of a specific IR influences the minimal rearrangements required for the gnepines and cupressophyte clades, whereby Pinaceae favours the removal of IRB but cupressophytes exclusion of IRA. This result strongly suggests that different IR copies have been lost from conifers I and II. Our data help understand the complexity and evolution of cupressophyte cpDNAs.

From algae to angiosperms-inferring the phylogeny of green plants (Viridiplantae) from 360 plastid genomes

BACKGROUND

Next-generation sequencing has provided a wealth of plastid genome sequence data from an increasingly diverse set of green plants (Viridiplantae). Although these data have helped resolve the phylogeny of numerous clades (e.g., green algae, angiosperms, and gymnosperms), their utility for inferring relationships across all green plants is uncertain. Viridiplantae originated 700-1500 million years ago and may comprise as many as 500,000 species. This clade represents a major source of photosynthetic carbon and contains an immense diversity of life forms, including some of the smallest and largest eukaryotes. Here we explore the limits and challenges of inferring a comprehensive green plant phylogeny from available complete or nearly complete plastid genome sequence data.

RESULTS

We assembled protein-coding sequence data for 78 genes from 360 diverse green plant taxa with complete or nearly complete plastid genome sequences available from GenBank. Phylogenetic analyses of the plastid data recovered well-supported backbone relationships and strong support for relationships that were not observed in previous analyses of major subclades within Viridiplantae. However, there also is evidence of systematic error in some analyses. In several instances we obtained strongly supported but conflicting topologies from analyses of nucleotides versus amino acid characters, and the considerable variation in GC content among lineages and within single genomes affected the phylogenetic placement of several taxa.

CONCLUSIONS

Analyses of the plastid sequence data recovered a strongly supported framework of relationships for green plants. This framework includes: i) the placement of Zygnematophyceace as sister to land plants (Embryophyta), ii) a clade of extant gymnosperms (Acrogymnospermae) with cycads + Ginkgo sister to remaining extant gymnosperms and with gnetophytes (Gnetophyta) sister to non-Pinaceae conifers (Gnecup trees), and iii) within the monilophyte clade (Monilophyta), Equisetales + Psilotales are sister to Marattiales + leptosporangiate ferns. Our analyses also highlight the challenges of using plastid genome sequences in deep-level phylogenomic analyses, and we provide suggestions for future analyses that will likely incorporate plastid genome sequence data for thousands of species. We particularly emphasize the importance of exploring the effects of different partitioning and character coding strategies.

Pinaceae-like reproductive morphology in Schizolepidopsis canicularis sp. nov. from the Early Cretaceous (Aptian-Albian) of Mongolia

PREMISE OF THE STUDY

Seed cone scales assigned to the genus Schizolepidopsis are widespread in Late Triassic to Cretaceous Eurasian deposits. They have been linked to the conifer family Pinaceae based on associated vegetative remains, but their exact affinities are uncertain. Recently discovered material from the Early Cretaceous of Mongolia reveals important new information concerning Schizolepidopsis cone scales and seeds, and provides support for a relationship between the genus and extant Pinaceae.

METHODS

Specimens were collected from Early Cretaceous (probable Aptian-Albian) lignite deposits in central Mongolia. Lignite samples were disaggregated, cleaned in hydrofluoric acid, and washed in water. Specimens were selected for further study using light and electron microscopy.

KEY RESULTS

Schizolepidopsis canicularis seed cones consist of loosely arranged, bilobed ovulate scales subtended by a small bract. A single inverted seed with an elongate micropyle is borne on each lobe of the ovulate scale. Each seed has a wing formed by the separation of the adaxial surface of the ovulate scale.

CONCLUSIONS

Schizolepidopsis canicularis produced winged seeds that formed in a manner that is unique to Pinaceae among extant conifers. We do not definitively place this species in Pinaceae pending more complete information concerning its pollen cones and vegetative remains. Nevertheless, this material suggests that Schizolepidopsis may be important for understanding the early evolution of Pinaceae, and may potentially help reconcile the appearance of the family in the fossil record with results based on phylogenetic analyses of molecular data.

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.

The complete chloroplast genome sequence of Mahonia bealei (Berberidaceae) reveals a significant expansion of the inverted repeat and phylogenetic relationship with other angiosperms

Mahonia bealei (Berberidaceae) is a frequently-used traditional Chinese medicinal plant with efficient anti-inflammatory ability. This plant is one of the sources of berberine, a new cholesterol-lowering drug with anti-diabetic activity. We have sequenced the complete nucleotide sequence of the chloroplast (cp) genome of M. bealei. The complete cp genome of M. bealei is 164,792 bp in length, and has a typical structure with large (LSC 73,052 bp) and small (SSC 18,591 bp) single-copy regions separated by a pair of inverted repeats (IRs 36,501 bp) of large size. The Mahonia cp genome contains 111 unique genes and 39 genes are duplicated in the IR regions. The gene order and content of M. bealei are almost unarranged which is consistent with the hypothesis that large IRs stabilize cp genome and reduce gene loss-and-gain probabilities during evolutionary process. A large IR expansion of over 12 kb has occurred in M. bealei, 15 genes (rps19, rpl22, rps3, rpl16, rpl14, rps8, infA, rpl36, rps11, petD, petB, psbH, psbN, psbT and psbB) have expanded to have an additional copy in the IRs. The IR expansion rearrangement occurred via a double-strand DNA break and subsequence repair, which is different from the ordinary gene conversion mechanism. Repeat analysis identified 39 direct/inverted repeats 30 bp or longer with a sequence identity ≥ 90%. Analysis also revealed 75 simple sequence repeat (SSR) loci and almost all are composed of A or T, contributing to a distinct bias in base composition. Comparison of protein-coding sequences with ESTs reveals 9 putative RNA edits and 5 of them resulted in non-synonymous modifications in rpoC1, rps2, rps19 and ycf1. Phylogenetic analysis using maximum parsimony (MP) and maximum likelihood (ML) was performed on a dataset composed of 65 protein-coding genes from 25 taxa, which yields an identical tree topology as previous plastid-based trees, and provides strong support for the sister relationship between Ranunculaceae and Berberidaceae. Molecular dating analyses suggest that Ranunculaceae and Berberidaceae diverged between 90 and 84 mya, which is congruent with the fossil records and with recent estimates of the divergence time of these two taxa.

Gymnosperms have increased capacity for electron leakage to oxygen (Mehler and PTOX reactions) in photosynthesis compared with angiosperms

Oxygen plays an important role in photosynthesis by participating in a number of O2-consuming reactions. O2 inhibits CO2 fixation by stimulating photorespiration, thus reducing plant production. O2 interacts with photosynthetic electron transport in the chloroplasts' thylakoids in two main ways: by accepting electrons from PSI (Mehler reaction); and by accepting electrons from reduced plastoquinone (PQ) mediated by the plastid terminal oxidase (PTOX). In this study, we show, using 101 plant species, that there is a difference in the potential for photosynthetic electron flow to O2 between angiosperms and gymnosperms. We found, from measurements of Chl fluorescence and leaf absorbance at 830 nm, (i) that electron outflow from PSII, as determined by decay kinetics of Chl fluorescence after application of a saturating light pulse, is more rapid in gymnosperms than in angiosperms; (ii) that the reaction center Chl of PSI (P700) is rapidly and highly oxidized in gymnosperms during induction of photosynthesis; and (iii) that these differences are dependent on oxygen. Finally, rates of O2 uptake measured by mass spectrometry in the absence of photorespiration were significantly promoted by illumination in dark-adapted leaves of gymnosperms, but not in those of angiosperms. The light-stimulated O2 uptake was around 10% of the maximum O2 evolution in gymnosperms and 1% in angiosperms. These results suggest that gymnosperms have increased capacity for electron leakage to oxygen in photosynthesis compared with angiosperms. The involvement of the Mehler reaction and PTOX in the electron flow to O2 is discussed.

Chloroplast phylogenomics indicates that Ginkgo biloba is sister to cycads

Molecular phylogenetic studies have not yet reached a consensus on the placement of Ginkgoales, which is represented by the only living species, Ginkgo biloba (common name: ginkgo). At least six discrepant placements of ginkgo have been proposed. This study aimed to use the chloroplast phylogenomic approach to examine possible factors that lead to such disagreeing placements. We found the sequence types used in the analyses as the most critical factor in the conflicting placements of ginkgo. In addition, the placement of ginkgo varied in the trees inferred from nucleotide (NU) sequences, which notably depended on breadth of taxon sampling, tree-building methods, codon positions, positions of Gnetopsida (common name: gnetophytes), and including or excluding gnetophytes in data sets. In contrast, the trees inferred from amino acid (AA) sequences congruently supported the monophyly of a ginkgo and Cycadales (common name: cycads) clade, regardless of which factors were examined. Our site-stripping analysis further revealed that the high substitution saturation of NU sequences mainly derived from the third codon positions and contributed to the variable placements of ginkgo. In summary, the factors we surveyed did not affect results inferred from analyses of AA sequences. Congruent topologies in our AA trees give more confidence in supporting the ginkgo-cycad sister-group hypothesis.

Identification of quantitative trait loci controlling fibre length and lignin content in Arabidopsis thaliana stems

Fibre properties and the biochemical composition of cell walls are important traits in many applications. For example, the lengths of fibres define the strength and quality of paper, and lignin content is a critical parameter for the use of biomass in biofuel production. Identifying genes controlling these traits is comparatively difficult in woody species, because of long generation times and limited amenability to high-resolution genetic mapping. To address this problem, this study mapped quantitative trait loci (QTLs) defining fibre length and lignin content in the Arabidopsis recombinant inbred line population Col-4 × Ler-0. Adapting high-throughput phenotyping techniques for both traits for measurements in Arabidopsis inflorescence stems identified significant QTLs for fibre length on chromosomes 2 and 5, as well as one significant QTL affecting lignin content on chromosome 2. For fibre length, total variation within the population was 208% higher than between parental lines and the identified QTLs explained 50.58% of the observed variation. For lignin content, the values were 261 and 26.51%, respectively. Bioinformatics analysis of the associated intervals identified a number of candidate genes for fibre length and lignin content. This study demonstrates that molecular mapping of QTLs pertaining to wood and fibre properties is possible in Arabidopsis, which substantially broadens the use of Arabidopsis as a model species for the functional characterization of plant genes.

Using multiple analytical methods to improve phylogenetic hypotheses in Minaria (Apocynaceae)

Metastelmatinae is a neotropical subtribe of Asclepiadoideae (Apocynaceae), comprising 13 genera and around 260 species whose phylogenetic relationships are often unresolved or incongruent between plastid and nuclear datasets. The genus Minaria is one of the first lineages to emerge in the Metastelmatinae and is highly supported based on plastid markers. It comprises 21 species, most of which are endemic to small areas with open vegetation in the Espinhaço Range, Brazil. In the work presented here, we use plastid (rps16, trnH-psbA, trnS-trnG, and trnD-trnT) and nuclear (ITS and ETS) datasets to investigate the relationships within Minaria. We show that the three methods mostly used in phylogenetic studies, namely, maximum parsimony, maximum likelihood, and Bayesian Inference, have different performances and that a pluralistic analytical approach combining results from them can increase tree resolution and clade confidence, providing valuable phylogenetic information.

Biosynthesis of amphetamine analogs in plants

Amphetamine analogs are produced by plants in the genus Ephedra and by Catha edulis, and include the widely used decongestants and appetite suppressants pseudoephedrine and ephedrine. A combination of yeast (Candida utilis or Saccharomyces cerevisiae) fermentation and subsequent chemical modification is used for the commercial production of these compounds. The availability of certain plant biosynthetic genes would facilitate the engineering of yeast strains capable of de novo pseudoephedrine and ephedrine biosynthesis. Chemical synthesis has yielded amphetamine analogs with myriad functional group substitutions and diverse pharmacological properties. The isolation of enzymes with the serendipitous capacity to accept novel substrates could allow the production of substituted amphetamines in synthetic biosystems. Here, we review the biology, biochemistry and biotechnological potential of amphetamine analogs in plants.

Developmental pattern of Ginkgo biloba levopimaradiene synthase (GbLPS) as probed by promoter analysis in Arabidopsis thaliana

Levopimaradiene synthase (GbLPS) of Ginkgo biloba catalyzes the first committed step in ginkgolide biosynthesis by converting geranylgeranyl diphosphate into levopimaradiene, which subsequently undergoes complex oxidation step and rearrangement of carbon skeleton, leading to formation of ginkgolides. To assess the organ-specificity and developmental characteristics of GbLPS expression, the GbLPS promoter-driven GUS expression in transgenic Arabidopsis was studied. Histological analysis of the transgenic Arabidopsis plant showed that the GUS accumulation was mainly localized in the epidermis of leaves, phloem of the shoots, ovaries and stamens of flowers, and vasculature of roots. These observations correlate with the occurrence of LPS transcripts in roots and male strobili of G. biloba. Treatment of methyl jasmonate on the transformant exhibited significant upregulation of the reporter gene in the roots with little change in leaves and flowers. The present findings support biosynthesis of ginkgolide in the roots of Ginkgo plant and suggest translocation occurs through the phloem.

Morphological evolution in land plants: new designs with old genes

The colonization and radiation of multicellular plants on land that started over 470 Ma was one of the defining events in the history of this planet. For the first time, large amounts of primary productivity occurred on the continental surface, paving the way for the evolution of complex terrestrial ecosystems and altering global biogeochemical cycles; increased weathering of continental silicates and organic carbon burial resulted in a 90 per cent reduction in atmospheric carbon dioxide levels. The evolution of plants on land was itself characterized by a series of radical transformations of their body plans that included the formation of three-dimensional tissues, de novo evolution of a multicellular diploid sporophyte generation, evolution of multicellular meristems, and the development of specialized tissues and organ systems such as vasculature, roots, leaves, seeds and flowers. In this review, we discuss the evolution of the genes and developmental mechanisms that drove the explosion of plant morphologies on land. Recent studies indicate that many of the gene families which control development in extant plants were already present in the earliest land plants. This suggests that the evolution of novel morphologies was to a large degree driven by the reassembly and reuse of pre-existing genetic mechanisms.

Extended linkage disequilibrium in noncoding regions in a conifer, Cryptomeria japonica

We measured linkage disequilibrium in mostly noncoding regions of Cryptomeria japonica, a conifer belonging to Cupressaceae. Linkage disequilibrium was extensive and did not decay even at a distance of 100 kb. The average estimate of the population recombination rate per base pair was 1.55 × 10(-5) and was <1/70 of that in the coding regions. We discuss the impact of low recombination rates in a large part of the genome on association studies.

Systematic error in seed plant phylogenomics

Resolving the closest relatives of Gnetales has been an enigmatic problem in seed plant phylogeny. The problem is known to be difficult because of the extent of divergence between this diverse group of gymnosperms and their closest phylogenetic relatives. Here, we investigate the evolutionary properties of conifer chloroplast DNA sequences. To improve taxon sampling of Cupressophyta (non-Pinaceae conifers), we report sequences from three new chloroplast (cp) genomes of Southern Hemisphere conifers. We have applied a site pattern sorting criterion to study compositional heterogeneity, heterotachy, and the fit of conifer chloroplast genome sequences to a general time reversible + G substitution model. We show that non-time reversible properties of aligned sequence positions in the chloroplast genomes of Gnetales mislead phylogenetic reconstruction of these seed plants. When 2,250 of the most varied sites in our concatenated alignment are excluded, phylogenetic analyses favor a close evolutionary relationship between the Gnetales and Pinaceae—the Gnepine hypothesis. Our analytical protocol provides a useful approach for evaluating the robustness of phylogenomic inferences. Our findings highlight the importance of goodness of fit between substitution model and data for understanding seed plant phylogeny.

Loss of different inverted repeat copies from the chloroplast genomes of Pinaceae and cupressophytes and influence of heterotachy on the evaluation of gymnosperm phylogeny

The relationships among the extant five gymnosperm groups--gnetophytes, Pinaceae, non-Pinaceae conifers (cupressophytes), Ginkgo, and cycads--remain equivocal. To clarify this issue, we sequenced the chloroplast genomes (cpDNAs) from two cupressophytes, Cephalotaxus wilsoniana and Taiwania cryptomerioides, and 53 common chloroplast protein-coding genes from another three cupressophytes, Agathis dammara, Nageia nagi, and Sciadopitys verticillata, and a non-Cycadaceae cycad, Bowenia serrulata. Comparative analyses of 11 conifer cpDNAs revealed that Pinaceae and cupressophytes each lost a different copy of inverted repeats (IRs), which contrasts with the view that the same IR has been lost in all conifers. Based on our structural finding, the character of an IR loss no longer conflicts with the "gnepines" hypothesis (gnetophytes sister to Pinaceae). Chloroplast phylogenomic analyses of amino acid sequences recovered incongruent topologies using different tree-building methods; however, we demonstrated that high heterotachous genes (genes that have highly different rates in different lineages) contributed to the long-branch attraction (LBA) artifact, resulting in incongruence of phylogenomic estimates. Additionally, amino acid compositions appear more heterogeneous in high than low heterotachous genes among the five gymnosperm groups. Removal of high heterotachous genes alleviated the LBA artifact and yielded congruent and robust tree topologies in which gnetophytes and Pinaceae formed a sister clade to cupressophytes (the gnepines hypothesis) and Ginkgo clustered with cycads. Adding more cupressophyte taxa could not improve the accuracy of chloroplast phylogenomics for the five gymnosperm groups. In contrast, removal of high heterotachous genes from data sets is simple and can increase confidence in evaluating the phylogeny of gymnosperms.

Morphological and functional stasis in mycorrhizal root nodules as exhibited by a Triassic conifer

Mycorrhizal root nodules occur in the conifer families Araucariaceae, Podocarpaceae, and Sciadopityaceae. Although the fossil record of these families can be traced back into the early Mesozoic, the oldest fossil evidence of root nodules previously came from the Cretaceous. Here we report on cellularly preserved root nodules of the early conifer Notophytum from Middle Triassic permineralized peat of Antarctica. These fossil root nodules contain fungal arbuscules, hyphal coils, and vesicles in their cortex. Numerous glomoid-type spores are found in the peat matrix surrounding the nodules. This discovery indicates that mutualistic associations between conifer root nodules and arbuscular mycorrhizal fungi date back to at least the early Mesozoic, the period during which most of the modern conifer families first appeared. Notophytum root nodules predate the next known appearance of this association by 100 million years, indicating that this specialized form of mycorrhizal symbiosis has ancient origins.

Chromosome number variation in two antipodean floras

BACKGROUND AND AIMS

We compared chromosome number (CN) variation in the nearly antipodean Italian and New Zealand floras to verify (i) whether patterns of variation reflect their similar latitudinal ranges or their different biogeographic/taxonomic contexts, (ii) if any differences are equally distributed across major taxa/lineages and (iii) if the frequency, number and taxonomic distribution of B-chromosomes differ between the two countries.

METHODOLOGY

We compared two datasets comprising 3426 (Italy) and 2525 (New Zealand) distinct cytotypes. We also compared a subset based on taxonomic orders and superimposed them onto a phylogeny of vascular plants. We used standard statistics, histograms, and either analysis of variance or Kruskal-Wallis tests to analyse the data.

PRINCIPAL RESULTS

Mean CN of the vascular New Zealand flora is about twice that of Italy. For most orders, mean CN values for New Zealand are higher than those of the Italian flora and the differences are statistically significant. Further differences in CN variation among the orders and main clades that we studied, irrespective of geographical distinctions, are revealed. No correlation was found between chromosome and B-chromosome number.

CONCLUSIONS

Mean CN of the whole New Zealand dataset is about twice that of the Italian flora. This suggests that extensive polyploidization played a major role in the evolution of the New Zealand vascular flora that is characterized by a rate of high endemism. Our results show that the hypothesis of a polyploid increase proportional to distance from the Equator cannot be applied to territories with the same latitudinal ranges but placed in different hemispheres. We suggest that bioclimatic gradients, rather than or in addition to latitudinal gradients, might account for a polyploidy increase. Our data also suggest that any adaptive role of B-chromosomes at geographic scale may be sought in their frequency rather than in their number.

Recent loss of plastid-encoded ndh genes within Erodium (Geraniaceae)

Plastid genomes in the flowering plant family Geraniaceae are known to be highly rearranged based on complete sequences representing the four major genera Erodium, Geranium, Monsonia, and Pelargonium. In this paper we report on the genome sequence of a second species of Erodium, E. carvifolium, representing the second major clade (clade II) in the phylogeny of this genus. Comparison of this genome sequence to the previously published sequence of E. texanum from clade I demonstrates that the plastid genomes of these two species encode the same number of proteins but differ greatly in their relative degree of rearrangement; 14 kb of additional sequence in E. texanum contains complex repeats associated with rearrangement endpoints, whereas the plastid genome of E. carvifolium is streamlined at 116 kb and displays no unique alterations in gene order. Furthermore, these species from both major clades of Erodium contain intact NADH dehydrogenase (ndh) genes, but the 11 ndh genes are represented as pseudogenes in a small clade of 13 species. It is unclear whether plastid-encoded ndh genes have been lost entirely or functionally transferred to the nucleus. This is the third report of the absence of functional ndh genes, and the current study describes the most recent loss of these genes among photosynthetic seed plants and the second such loss among angiosperms. The other ndh losses from Pinaceae/Gnetales and Orchidaceae are much more ancient. Comparative biochemistry between Erodium species with and without plastid-encoded ndh genes may elucidate changes in photosynthetic function and the role of the Ndh complex.

Multigene phylogeny of the green lineage reveals the origin and diversification of land plants

The Viridiplantae (green plants) include land plants as well as the two distinct lineages of green algae, chlorophytes and charophytes. Despite their critical importance for identifying the closest living relatives of land plants, phylogenetic studies of charophytes have provided equivocal results [1-5]. In addition, many relationships remain unresolved among the land plants, such as the position of mosses, liverworts, and the enigmatic Gnetales. Phylogenomics has proven to be an insightful approach for resolving challenging phylogenetic issues, particularly concerning deep nodes [6-8]. Here we extend this approach to the green lineage by assembling a multilocus data set of 77 nuclear genes (12,149 unambiguously aligned amino acid positions) from 77 taxa of plants. We therefore provide the first multigene phylogenetic evidence that Coleochaetales represent the closest living relatives of land plants. Moreover, our data reinforce the early divergence of liverworts and the close relationship between Gnetales and Pinaceae. These results provide a new phylogenetic framework and represent a key step in the evolutionary interpretation of developmental and genomic characters in green plants.

Organization, anatomy, and fungal endophytes of a Triassic conifer embryo

PREMISE OF THE STUDY

Much is known about the Paleozoic conifers and the conifers assignable to extant families that appear in the Jurassic; however, relatively little is known about the transitional conifers of the early Mesozoic, especially the Voltziales. To better understand the evolution of this group, we aim to increase knowledge of voltzialean anatomy and morphology. A fossil embryo of the group is described in this study. •

METHODS

Several permineralized seeds, one containing a well-preserved embryo, were collected from the early Middle Triassic Fremouw Peak locality in the Central Transantarctic Mountains. Samples were prepared using the standard acetate peel technique and studied in transmitted light. •

KEY RESULTS

The embryo belongs to Parasciadopitys aequata, a member of the Voltziales. The embryo and megagametophyte tissues are exquisitely preserved. The embryo is colonized by two distinct fungi. Sporocarps of the fungi are found in the megagametophyte and in the space between the megagametophyte and nucellus. The additions of embryo and megagametophyte characters to the description of P. aequata have made it one of the most completely known fossil taxa reproductively. •

CONCLUSIONS

Preservation of fossil embryos, although rare, can expand the array of characters available in tracing the evolutionary history of plants. The embryo of P. aequata shares similarities with embryos of other extinct and extant conifer families. The association of the embryo with Combresomyces cornifer and Mycocarpon asterineum increases our understanding of the roles of microorganisms in Triassic ecosystems.

Comparative expression pattern analysis of WUSCHEL-related homeobox 2 (WOX2) and WOX8/9 in developing seeds and somatic embryos of the gymnosperm Picea abies

• In seed plants, current knowledge concerning embryonic pattern formation by polar auxin transport (PAT) and WUSCHEL-related homeobox (WOX) gene activity is primarily derived from studies on angiosperms, while less is known about these processes in gymnosperms. In view of the differences in their embryogeny, and the fact that somatic embryogenesis is used for mass propagation of conifers, a better understanding of embryo development is vital. • The expression patterns of PaWOX2 and PaWOX8/9 were followed with quantitative reverse transcription-polymerase chain reaction (qRT-PCR) and in situ hybridization (ISH) during seed and somatic embryo development in Norway spruce (Picea abies), and in somatic embryos treated with the PAT inhibitor N-1-naphthylphthalamic acid (NPA). • Both PaWOX2 and PaWOX8/9 were highly expressed at the early growth stages of zygotic and somatic embryos, and shared a similar expression pattern over the entire embryo. At later embryo stages, high expression of PaWOX8/9 became restricted to cotyledon primordia, epidermis, procambium and root apical meristem (RAM), which became most evident in NPA-treated somatic embryos, while expression of PaWOX2 was much lower. • Our results suggest an ancestral role of WOX in seed plant embryo development, and strengthen the proposed connection between PAT, PIN-FORMED (PIN) and WOX in the regulation of embryo patterning in seed plants.

Molecular interactions of orthologues of floral homeotic proteins from the gymnosperm Gnetum gnemon provide a clue to the evolutionary origin of 'floral quartets'

Several lines of evidence suggest that the identity of floral organs in angiosperms is specified by multimeric transcription factor complexes composed of MADS-domain proteins. These bind to specific cis-regulatory elements ('CArG-boxes') of their target genes involving DNA-loop formation, thus constituting 'floral quartets'. Gymnosperms, angiosperms' closest relatives, contain orthologues of floral homeotic genes, but when and how the interactions constituting floral quartets were established during evolution has remained unknown. We have comprehensively studied the dimerization and DNA-binding of several classes of MADS-domain proteins from the gymnosperm Gnetum gnemon. Determination of protein-protein and protein-DNA interactions by yeast two-hybrid, in vitro pull-down and electrophoretic mobility shift assays revealed complex patterns of homo- and heterodimerization among orthologues of floral homeotic class B, class C and class E proteins and B(sister) proteins. Using DNase I footprint assays we demonstrate that both orthologues of class B with C proteins, and orthologues of class C proteins alone, but not orthologues of class B proteins alone can loop DNA in floral quartet-like complexes. This is in contrast to class B and class C proteins from angiosperms, which require other factors such as class E floral homeotic proteins to 'glue' them together in multimeric complexes. Our findings suggest that the evolutionary origin of floral quartet formation is based on the interaction of different DNA-bound homodimers, does not depend on class E proteins, and predates the origin of angiosperms.

Plastid ndh genes in plant evolution

The plastid ndh genes encode components of the thylakoid Ndh complex which purportedly acts as an electron feeding valve to adjust the redox level of the cyclic photosynthetic electron transporters. During the process of evolution from endosymbiosis to modern chloroplast, most cyanobacterial genes were lost or transferred to nucleus. Eleven ndh genes are among the 150-200 genes remaining in higher plant chloroplast DNA, out of some 3000 genes in the original prokaryotic Cyanobacteria in which homologues to ndh genes encode components of the respiratory Complex I and probably other complexes. The ndh genes are absent in all sequenced plastid DNAs of algae except for the Charophyceae and some Prasinophyceae. With the possible exclusion of some Conifers and Gnetales, the plastid DNA of all photosynthetic land plants contains the ndh genes, whereas they are absent in epiphytic plants that have also lost genes for the photosynthetic machinery. Therefore, the functional role of the ndh genes seems closely related to the land adaptation of photosynthesis. Transcripts of several plastid genes require C to U editing. The ndh genes concentrate about 50% of the editing sites of angiosperm plastid transcripts. Editing sites may be remnants from an ancestor in which a number of T to C inactivating mutations took place in the ndh genes which, during evolution, are being corrected back to T. The comparison of homologous editing sites in the mRNAs of angiosperm ndh genes provides a tool to investigate selective and permissive environmental conditions of past evolutionary events.

Lineage-specific group II intron gains and losses of the mitochondrial rps3 gene in gymnosperms

According to PCR assays and sequencing, we now report the shared presence of two rps3 introns, namely the rps3i74 and the rps3i249, in the mitochondria of all the classes representing the surviving lineages of gymnosperms, and unveil several lineages experiencing intron loss. Interestingly, the rps3 intron gains and losses within the four groups of gymnosperms let us sort out the Pinaceae and the non-Pinaceae into intron (+)- and intron (-)-lineages, respectively. Worthy of mention is also the finding that only Gnetum within the Gnetales harbours both the rps3 introns. This intron distribution pattern is consistent with the hypothesis that the two rps3 introns were likely present in the common ancestor of the seed plants and, then, independently lost in the non-Pinaceae during gymnosperm evolution. The derived secondary structural model of the novel group IIA intron improves our understanding of the significance and origin of the extraordinary length polymorphisms observed among rps3i249 orthologs. Despite the remarkable structural plasticity to adopt and reject introns, the rps3 mRNAs undergo accurate processing by splicing and extensive editing in gymnosperm mitochondria. This study provides additional insights into the evolutionarily high dynamics of mitochondrial introns which may come and go in closely related plant species. The turnover of the mitochondrial rps3 group II introns seen among lineages of seed plants further suggests that these introns might be an additional signature to discriminate between particularly cryptical taxonomic groups for which there is a need of a further evaluation of their evolutionary affiliation.

The position of gnetales among seed plants: overcoming pitfalls of chloroplast phylogenomics

The phylogenetic position of Gnetales is one of the most contentious issues in the seed plant systematics. To elucidate the Gnetales position, an improved amino acid substitution matrix was estimated based on 64 chloroplast (cp) genomes and was applied to cp genome data including all three lineages of Gnetales in maximum likelihood analyses of proteins. Although the initial analysis strongly supported the sister relation of Gnetales with Cryptomeria (Cupressophyta or non-Pinaceae conifers) (the "Gnecup" hypothesis), the support seems to be caused by a long-branch attraction (LBA) artifact. Indeed, by removing fastest evolving proteins that are most likely associated with the LBA, the support drastically declined. Furthermore, another analysis of partial genome data with dense taxon sampling of conifers showed that, in psbC, rpl2, and rps7 proteins, there are many parallel amino acid substitutions between the lineages leading to Gnetales and to Cryptomeria, and by further excluding these three genes, the sister relation of Gnetales with Pinaceae (the "Gnepine" hypothesis) became supported. Overall, our analyses indicate that the LBA and parallel substitutions cause a seriously biased inference of phylogenetic position of Gnetales with the cp genome data.

Expression of a gymnosperm PIN homologous gene correlates with auxin immunolocalization pattern at cotyledon formation and in demarcation of the procambium during Picea abies somatic embryo development and in seedling tissues

In seed plants, the body organization is established during embryogenesis and is uniform across gymnosperms and angiosperms, despite differences during early embryogeny. Evidence from angiosperms implicates the plant hormone auxin and its polar transport, mainly established by the PIN family of auxin efflux transporters, in the patterning of embryos. Here, PaPIN1 from Norway spruce (Picea abies [L.] Karst.), a gene widely expressed in conifer tissues and organs, was characterized and its expression and localization patterns were determined with reverse transcription polymerase chain reaction and in situ hybridization during somatic embryo development and in seedlings. PaPIN1 shares the predicted structure of other PIN proteins, but its central hydrophilic loop is longer than most PINs. In phylogenetic analyses, PaPIN1 clusters with Arabidopsis thaliana (L.) Heynh. PIN3, PIN4 and PIN7, but its expression pattern also suggests similarity to PIN1. The PaPIN1 expression signal was high in the protoderm of pre-cotyledonary embryos, but not if embryos were pre-treated with the auxin transport inhibitor N-1-naphthylphthalamic acid (NPA). This, together with a high auxin immunolocalization signal in this cell layer, suggests a role of PaPIN1 during cotyledon formation. At later stages, high PaPIN1 expression was observed in differentiating procambium, running from the tip of incipient cotyledons down through the embryo axis and to the root apical meristem (RAM), although the mode of RAM specification in conifer embryos differs from that of most angiosperms. Also, the PaPIN1 in situ signal was high in seedling root tips including root cap columella cells. The results thus suggest that PaPIN1 provides an ancient function associated with auxin transport and embryo pattern formation prior to the separation of angiosperms and gymnosperms, in spite of some morphological differences.

A duplicate gene rooting of seed plants and the phylogenetic position of flowering plants

Flowering plants represent the most significant branch in the tree of land plants, with respect to the number of extant species, their impact on the shaping of modern ecosystems and their economic importance. However, unlike so many persistent phylogenetic problems that have yielded to insights from DNA sequence data, the mystery surrounding the origin of angiosperms has deepened with the advent and advance of molecular systematics. Strong statistical support for competing hypotheses and recent novel trees from molecular data suggest that the accuracy of current molecular trees requires further testing. Analyses of phytochrome amino acids using a duplicate gene-rooting approach yield trees that unite cycads and angiosperms in a clade that is sister to a clade in which Gingko and Cupressophyta are successive sister taxa to gnetophytes plus Pinaceae. Application of a cycads + angiosperms backbone constraint in analyses of a morphological dataset yields better resolved trees than do analyses in which extant gymnosperms are forced to be monophyletic. The results have implications both for our assessment of uncertainty in trees from sequence data and for our use of molecular constraints as a way to integrate insights from morphological and molecular evidence.

The female reproductive unit of ephedra (Gnetales): comparative morphology and evolutionary perspectives

Morphological variation in Ephedra (Gnetales) is limited and confusing from an evolutionary perspective, with parallelisms and intraspecific variation. However, recent analyses of molecular data provide a phylogenetic framework for investigations of morphological traits, albeit with few informative characters in the investigated gene regions. We document morphological, anatomical and histological variation patterns in the female reproductive unit and test the hypothesis that some Early Cretaceous fossils, which share synapomorphies with Ephedra, are members of the extant clade. Results indicate that some morphological features are evolutionarily informative although intraspecific variation is evident. Histology and anatomy of cone bracts and seed envelopes show clade-specific variation patterns. There is little evidence for an inclusion of the Cretaceous fossils in the extant clade. Rather, a hypothesized general pattern of reduction of the vasculature in the ephedran seed envelope, probably from four vascular bundles in the fossils, to ancestrally three in the living clade, and later to two, is consistent with phylogenetic and temporal analyses, which indicate that extant diversity evolved after the Cretaceous-Tertiary boundary. Notwithstanding striking similarities between living and Cretaceous Ephedra, available data indicate that the Mesozoic diversity went almost entirely extinct in the late Cretaceous causing a bottleneck effect in Ephedra, still reflected today by an extraordinarily low level of genetic and structural diversity.

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.

Analysis of the soluble cell wall proteome of gymnosperms

We analyzed the cell wall proteome of lignifying suspension cell cultures (SCCs) from four gymnosperms that differ in evolution degree. This analysis showed the presence of "peptide sequence tags" (PSTs) corresponding to glucan endo-1,3-beta-D-glucosidase, xyloglucan-endotrans-glucosylase/hydrolase, chitinases, thaumatin-like proteins and proteins involved in lignin/lignan biosynthesis, such as dirigent-like proteins and peroxidases. Surprisingly, and given the abundance of peroxidases in the cell wall proteome of these gymnosperms, PSTs corresponding to peroxidases were only detected in tryptic fragments of the cell wall proteome of Cycas revoluta. The current lack of knowledge regarding C. revoluta peroxidases led us to purify, characterize and partially sequence the peroxidases responsible for lignin biosynthesis in this species. This yielded three peroxidase-enriched fractions: CrPrx 1, CrPrx 2 and CrPrx 3. Analyses of tryptic peptides of CrPrx 2 (32kDa) and CrPrx 3 (26kDa) suggest that CrPrx 3 arises from CrPrx 2 by protein truncation, and that CrPrx 3 apparently constitutes a post-translational modification of CrPrx 2. That CrPrx 2 and CrPrx 3 are apparently the same enzyme was also deduced from the similarity between the k(cat) shown by both peroxidases for the three monolignols. These results emphasize the analogies between the cell wall proteome of gymnosperms and angiosperms, the complexity of the peroxidase proteome, and the difficulties involved in establishing fine structure-function relationships.

Distribution and evolutionary trends of photoprotective isoprenoids (xanthophylls and tocopherols) within the plant kingdom

The earliest land photosynthesis would have increased the risk of photo-oxidations and the demand of anti-oxidative protection. In this work, we aimed to determine the evolutionary trends in photoprotection across a wide representation of the plant kingdom and to verify whether the non-ubiquitous lutein-epoxide (Lx) cycle is a polyphyletic or an ancient character. Carotenoids and alpha-tocopherol (alpha-toc) were analysed by HPLC in 266 species. Phylogenetic analyses of the presence of photoprotective compounds and zeaxanthin-epoxidase (ZE) sequences were performed. Violaxanthin-cycle pigments (VAZ) and alpha-toc were taxonomically ubiquitous. Ancient groups showed higher contents of VAZ than vascular plants, while alpha-toc showed the opposite pattern. Lutein-epoxide was present in 45% of the species. It showed a remarkable variation across groups but with a clear increasing trend from algae to basal angiosperms. Lutein-epoxide was also related to woody trait and leaf longevity. No correlation between the presence of Lx and recurrent mutations in ZE sequences, including the duplications, was found. Thus, there is an evolutionary trend to increase the content of alpha-toc and to decrease the total amount of VAZ pigments. Absence of Lx in algae discards an ancestral origin. Present results are also inconsistent with a polyphyletic origin of Lx in angiosperms.

The relative ages of ectomycorrhizal mushrooms and their plant hosts estimated using Bayesian relaxed molecular clock analyses

BACKGROUND

Ectomycorrhizae (ECM) are symbioses formed by polyphyletic assemblages of fungi (mostly Agaricomycetes) and plants (mostly Pinaceae and angiosperms in the rosid clade). Efforts to reconstruct the evolution of the ECM habit in Agaricomycetes have yielded vastly different results, ranging from scenarios with many relatively recent origins of the symbiosis and no reversals to the free-living condition; a single ancient origin of ECM and many subsequent transitions to the free-living condition; or multiple gains and losses of the association. To test the plausibility of these scenarios, we performed Bayesian relaxed molecular clock analyses including fungi, plants, and other eukaryotes, based on the principle that a symbiosis cannot evolve prior to the origin of both partners. As we were primarily interested in the relative ages of the plants and fungi, we did not attempt to calibrate the molecular clock using the very limited fossil record of Agaricomycetes.

RESULTS

Topologically constrained and unconstrained analyses suggest that the root node of the Agaricomycetes is much older than either the rosids or Pinaceae. The Agaricomycetidae, a large clade containing the Agaricales and Boletales (collectively representing 70% of Agaricomycetes), is also significantly older than the rosids. The relative age of Agaricomycetidae and Pinaceae, however, is sensitive to tree topology, and the inclusion or exclusion of the gnetophyte Welwitschia mirabilis.

CONCLUSION

The ancestor of the Agaricomycetes could not have been an ECM species because it existed long before any of its potential hosts. Within more derived clades of Agaricomycetes, there have been at least eight independent origins of ECM associations involving angiosperms, and at least six to eight origins of associations with gymnosperms. The first ECM symbioses may have involved Pinaceae, which are older than rosids, but several major clades of Agaricomycetes, such as the Boletales and Russulales, are young enough to have been plesiomorphically associated with either rosids or Pinaceae, suggesting that some contemporary ECM partnerships could be of very ancient origin.

The presence of sinapyl lignin in Ginkgo biloba cell cultures changes our views of the evolution of lignin biosynthesis

Suspension cell cultures (SCCs) from one of the oldest seed plants, Ginkgo biloba, show unpredictable alterations in the nature of the lignins, such as is the recruitment of sinapyl alcohol for lignin biosynthesis, compared with the woody tissues of the same species, which lack syringyl (S) lignins. These results show that, in this gymnosperm, the genes involved in sinapyl alcohol biosynthesis are latent and that their regulatory regions respond, by initiating gene expression, to the developmental signals and the environmental clues, which condition its in vitro culture. G. biloba SCCs not only synthesize S lignins but also their extracellular proteome contains both class III peroxidases capable of oxidizing sinapyl alcohol and enzymes involved in H2O2 production, observation which suggests that the peroxidase branch for the oxidative coupling of sinapyl alcohol units into lignins is operative. The incomplete knowledge of the G. biloba peroxidase-encoding genes led us to purify, characterize and partially sequence the peroxidase responsible for monolignol oxidation. When the major peroxidase from G. biloba SCCs (GbPrx) was purified to homogeneity, it showed absorption maxima in the visible region at 414 (Soret band), and at 543 and 570 nm, which calls to mind those shown by low-spin ferric peroxidases. However, the results also showed that the paraperoxidase-like character of GbPrx is not an obstacle for oxidizing the three monolignols compared with high-spin ferric peroxidases. Taken together, these results mean that the time at which the evolutionary gain of the segment of the route that leads to the biosynthesis of S lignins took place in seed plants needs to be revised.

Evolution of reduced and compact chloroplast genomes (cpDNAs) in gnetophytes: selection toward a lower-cost strategy

The cpDNA of Welwitschia mirabilis (the only species of Welwitschiales) was recently reported to be the most reduced and compact among photosynthetic land plants. However, cpDNAs of the other two gnetophyte lineages (viz. Ephedrales and Gnetales) have not yet been studied. It remains unclear what underlining mechanisms have downsized the cpDNA. To pin down major factors for cpDNA reduction and compaction in gnetophytes, we have determined 4 complete cpDNAs, including one from each of the 3 gnetophyte orders, Ephedra equisetina, Gnetum parvifolium, and W. mirabilis, and one from the non-Pinus Pinaceae, Keteleeria davidiana. We report that the cpDNAs of E. equisetina (109,518bp) and G.parvifolium (114,914bp) are not only smaller but more compact than that of W. mirabilis (118,919bp). The gnetophyte cpDNAs have commonly lost at least 18 genes that are retained in other seed plants. Furthermore, they have significantly biased usages of AT-rich codons and shorter introns and intergenic spaces, which are largely due to more deletions at inter-operon than intra-operon spaces and removal of segment sequences rather than single-nucleotides. We show that the reduced gnetophyte cpDNAs clearly resulted from selection for economy by deletions of genes and non-coding sequences, which then led to the compactness and the accelerated substitution rates. The smallest C-values in gnetophyte nuclear DNAs and the competitive or resource-poor situations encountered by gnetophytes further suggest a critical need for an economic strategy.

Is the anthophyte hypothesis alive and well? New evidence from the reproductive structures of Bennettitales

Bennettitales is an extinct group of seed plants with reproductive structures that are similar in some respects to both Gnetales and angiosperms, but systematic relationships among the three clades remain controversial. This study summarizes characters of bennettitalean plants and presents new evidence for the structure of cones and seeds that help clarify relationships of Bennettitales to flowering plants, Gnetales, and other potential angiosperm sister groups. Bennettitales have simple mono- or bisporangiate cones. Seeds are borne terminally on sporophylls. They have a unique structure that includes a nucellus with a solid apex, no pollen chamber, and a single integument, and they are clearly not enclosed by a cupule or other specialized structures. Such features differ substantially from Gnetales, flowering plants, and the seed fern Caytonia, providing no compelling evidence for the origin of the angiospermous carpel. Cladistic tests were performed to assess the strength of the "anthophyte hypothesis" and possible relationships of Bennettitales, Gnetales, and Caytonia to flowering plants. Our results do not support the anthophyte hypothesis for the origin of angiosperms by a transformation of fertile organs that were already aggregated into a cone or flower-like structure. However, the anthophyte topology of the seed plant tree continues to be supported by morphological analyses of living and extinct taxa.

Phylogenetic relationships among seed plants: Persistent questions and the limits of molecular data

Trees inferred from DNA sequence data provide only limited insight into the phylogeny of seed plants because the living lineages (cycads, Ginkgo, conifers, gnetophytes, and angiosperms) represent fewer than half of the major lineages that have been detected in the fossil record. Nevertheless, phylogenetic trees of living seed plants inferred from sequence data can provide a test of relationships inferred in analyses that include fossils. So far, however, significant uncertainty persists because nucleotide data support several conflicting hypotheses. It is likely that improved sampling of gymnosperm diversity in nucleotide data sets will help alleviate some of the analytical issues encountered in the estimation of seed plant phylogeny, providing a more definitive test of morphological trees. Still, rigorous morphological analyses will be required to answer certain fundamental questions, such as the identity of the angiosperm sister group and the rooting of crown seed plants. Moreover, it will be important to identify approaches for incorporating insights from data that may be accurate but less likely than sequence data to generate results supported by high bootstrap values. How best to weigh evidence and distinguish among hypotheses when some types of data give high support values and others do not remains an important problem.

The meaning of Darwin's 'abominable mystery'

Charles Darwin's "abominable mystery" has come to symbolize just about all aspects of the origin and early evolution of flowering plants. Yet, there has never been an analysis of precisely what Darwin thought was so abominably mysterious. Here I explicate Darwin's thoughts and frustrations with the fossil record of flowering plants as revealed in correspondence with Joseph Hooker, Gaston de Saporta, and Oswald Heer between 1875 and 1881. I also examine the essay by John Ball that prompted Darwin to write his "abominable mystery" letter to Hooker in July of 1879. Contrary to what is generally believed, Darwin's abominable mystery has little if anything to do with the fossil prehistory of angiosperms, identification of the closest relatives of flowering plants, questions of the homologies (and character transformations) of defining features of flowering plants, or the phylogeny of flowering plants themselves. Darwin's abominable mystery and his abiding interest in the radiation of angiosperms were never driven primarily by a need to understand the literal text of the evolutionary history of flowering plants. Rather, Darwin was deeply bothered by what he perceived to be an abrupt origin and highly accelerated rate of diversification of flowering plants in the mid-Cretaceous. This led Darwin to create speculative arguments for a long, gradual, and undiscovered pre-Cretaceous history of flowering plants on a lost island or continent. Darwin also took refuge in the possibility that a rapid diversification of flowering plants in the mid-Cretaceous might, if real, have a biological explanation involving coevolutionary interactions between pollinating insects and angiosperms. Nevertheless, although generations of plant biologists have seized upon Darwin's abominable mystery as a metaphor for their struggle to understand angiosperm history, the evidence strongly suggests that the abominable mystery is not about angiosperms per se. On the contrary, Darwin's abominable mystery is about his abhorrence that evolution could be both rapid and potentially even saltational. Throughout the last years of his life, it just so happens that flowering plants, among all groups of organisms, presented Darwin with the most extreme exception to his strongly held notion natura non facit saltum, nature does not make a leap.