Relationships among seed plants inferred from highly conserved genes: sorting conflicting phylogenetic signals among ancient lineages

Model Species to Investigate the Origin of Flowers

The angiosperms, or flowering plants, arose at least 135 million years ago (Ma) and rapidly diversified to form over 300,000 species alive today. This group appears, however, to have separated from its closest living relatives, the extant gymnosperms, much earlier: over 300 Ma. Representatives of basally-diverging angiosperm lineages are of key importance to studies aimed at reconstructing the most recent common ancestor of living angiosperms, including its morphological, anatomical, eco-physiological and molecular aspects. Furthermore, evo-devo comparisons of angiosperms with living gymnosperms may help to determine how the many novel aspects of angiosperms, including those of the flower, first came about. This chapter reviews literature on the origin of angiosperms and focusses on basally-diverging angiosperms and gymnosperms that show advantages as potential experimental models, reviewing information and protocols for the use of these species in an evo-devo context. The final section suggests a means by which data from living and fossil groups could be integrated to better elucidate evolutionary events that took place on the long stem-lineage that apparently preceded the radiation of living angiosperms.

Cutting the long branches: Consilience as a path to unearth the evolutionary history of Gnetales

The Gnetales are one of the most fascinating groups within seed plants. Although the advent of molecular phylogenetics has generated some confidence in their phylogenetic placement of Gnetales within seed plants, their macroevolutionary history still presents many unknowns. Here, we review the reasons for such unknowns, and we focus the discussion on the presence of “long branches” both in their molecular and morphological history. The increased rate of molecular evolution and genome instability as well as the numerous unique traits (both reproductive and vegetative) in the Gnetales have been obstacles to a better understanding of their evolution. Moreover, the fossil record of the Gnetales, though relatively rich, has not yet been properly reviewed and investigated using a phylogenetic framework. Despite these apparent blocks to progress we identify new avenues to enable us to move forward. We suggest that a consilience approach, involving different disciplines such as developmental genetics, paleobotany, molecular phylogenetics, and traditional anatomy and morphology might help to “break” these long branches, leading to a deeper understanding of this mysterious group of plants.

Evolution of the coniferous seed scale

BACKGROUND

The Florin model is the commonly accepted theory of coniferous seed scale evolution. It describes the derivation of extant seed scale morphology from the morphology of fossil conifers via the reduction of complex to simple axillary structures. In this framework the seed scale is composed of a reduced lateral shoot with fertile and sterile appendages which are interpreted as leaf homologues.

SCOPE

The Florin model has three crucial problems that we address here: (1) the original derivation series does not take the ontogeny of extant conifers into account, (2) it cannot explain the morphology of all extant conifers and (3) Taxaceae were originally excluded. Examination of seed cones of extant conifers shows that ovules occur in three different positions in the cone: in an axillary position, replacing a leaf or terminating the cone axis. By interpreting the fertile appendage or seed-bearing structure as a leaf, not all positions are possible. The exclusion of Taxaceae from conifers is in stark contrast to recent molecular phylogenetic studies, which include Taxaceae in conifers as sister to Cupressaceae. Therefore, the Florin model does not offer an adequate explanation for taxaceous morphology.

CONCLUSION

We conclude that the seed-bearing structure of conifers cannot be interpreted as homologous to a leaf. In the interpretation we present here, the seed-bearing structure is the modified funiculus of the ovule, multiples of which laterally fuse to form the seed scale. The seed scales of all extant conifers can be derived from a Cunninghamia-like morphology via fusion and reduction of individual funiculi.

Plastid phylogenomic analysis of green plants: A billion years of evolutionary history

PREMISE OF THE STUDY

For the past one billion years, green plants (Viridiplantae) have dominated global ecosystems, yet many key branches in their evolutionary history remain poorly resolved. Using the largest analysis of Viridiplantae based on plastid genome sequences to date, we examined the phylogeny and implications for morphological evolution at key nodes.

METHODS

We analyzed amino acid sequences from protein-coding genes from complete (or nearly complete) plastomes for 1879 taxa, including representatives across all major clades of Viridiplantae. Much of the data used was derived from transcriptomes from the One Thousand Plants Project (1KP); other data were taken from GenBank.

KEY RESULTS

Our results largely agree with previous plastid-based analyses. Noteworthy results include (1) the position of Zygnematophyceae as sister to land plants (Embryophyta), (2) a bryophyte clade (hornworts, mosses + liverworts), (3) Equisetum + Psilotaceae as sister to Marattiales + leptosporangiate ferns, (4) cycads + Ginkgo as sister to the remaining extant gymnosperms, within which Gnetophyta are placed within conifers as sister to non-Pinaceae (Gne-Cup hypothesis), and (5) Amborella, followed by water lilies (Nymphaeales), as successive sisters to all other extant angiosperms. Within angiosperms, there is support for Mesangiospermae, a clade that comprises magnoliids, Chloranthales, monocots, Ceratophyllum, and eudicots. The placements of Ceratophyllum and Dilleniaceae remain problematic. Within Pentapetalae, two major clades (superasterids and superrosids) are recovered.

CONCLUSIONS

This plastid data set provides an important resource for elucidating morphological evolution, dating divergence times in Viridiplantae, comparisons with emerging nuclear phylogenies, and analyses of molecular evolutionary patterns and dynamics of the plastid genome.

Fossil record of Ephedra in the Lower Cretaceous (Aptian), Argentina

Fossil plants from the Lower Cretaceous (upper Aptian) of the La Cantera Formation, Argentina, are described. The fossils studied represent a leafy shooting system with several orders of articulated and striated axes and attached leaves with unequivocal ephedroid affinity. We also found associated remains of ovulate cones with four whorls of sterile bracts, which contain two female reproductive units (FRU). Ovulate cone characters fit well within the genus Ephedra. Special characters in the ovulate cones including an outer seed envelope with two types of trichomes, allowed us to consider our remains as a new Ephedra species. Abundant dispersed ephedroid pollen obtained from the macrofossil-bearing strata also confirms the abundance of Ephedraceae in the basin. The co-occurrence of abundant fossil of Ephedra (adapted to dry habitats) associated with thermophilic cheirolepideacean conifer pollen (Classopollis) in the unit would suggest marked seasonality at the locality during the Early Cretaceous. Furthermore, the floristic association is linked to dry sensitive rocks in the entire section. The macro- and microflora from San Luis Basin are similar in composition to several Early Cretaceous floras from the Northern Gondwana floristic province, but it may represent one of the southernmost records of an arid biome in South America.

Evidence for an ancient whole genome duplication in the cycad lineage

Contrary to the many whole genome duplication events recorded for angiosperms (flowering plants), whole genome duplications in gymnosperms (non-flowering seed plants) seem to be much rarer. Although ancient whole genome duplications have been reported for most gymnosperm lineages as well, some are still contested and need to be confirmed. For instance, data for ginkgo, but particularly cycads have remained inconclusive so far, likely due to the quality of the data available and flaws in the analysis. We extracted and sequenced RNA from both the cycad Encephalartos natalensis and Ginkgo biloba. This was followed by transcriptome assembly, after which these data were used to build paralog age distributions. Based on these distributions, we identified remnants of an ancient whole genome duplication in both cycads and ginkgo. The most parsimonious explanation would be that this whole genome duplication event was shared between both species and had occurred prior to their divergence, about 300 million years ago.

Specialization and generalization in the diversification of phytophagous insects: tests of the musical chairs and oscillation hypotheses

Evolutionary biologists have often assumed that ecological generalism comes at the expense of less intense exploitation of specific resources and that this trade-off will promote the evolution of ecologically specialized daughter species. Using a phylogenetic comparative approach with butterflies as a model system, we test hypotheses that incorporate changes in niche breadth and location into explanations of the taxonomic diversification of insect herbivores. Specifically, we compare the oscillation hypothesis, where speciation is driven by host-plant generalists giving rise to specialist daughter species, to the musical chairs hypothesis, where speciation is driven by host-plant switching, without changes in niche breadth. Contrary to the predictions of the oscillation hypothesis, we recover a negative relationship between host-plant breadth and diversification rate and find that changes in host breadth are seldom coupled to speciation events. By contrast, we present evidence for a positive relationship between rates of host switching and butterfly diversification, consonant with the musical chairs hypothesis. These results suggest that the costs of trophic generalism in plant-feeding insects may have been overvalued and that transitions from generalists to ecological specialists may not be an important driver of speciation in general.

Impact of gene molecular evolution on phylogenetic reconstruction: a case study in the rosids (Superorder Rosanae, Angiosperms)

Rate of substitution of genomic regions is among the most debated intrinsic features that impact phylogenetic informativeness. However, this variable is also coupled with rates of nonsynonymous substitutions that underscore the nature and degree of selection on the selected genes. To empirically address these variables, we constructed four completely overlapping data sets of plastid matK, atpB, rbcL, and mitochondrial matR genes and used the rosid lineage (angiosperms) as a working platform. The genes differ in combinations of overall rates of nucleotide and amino acid substitutions. Tree robustness, homoplasy, accuracy in contrast to a reference tree, and phylogenetic informativeness are evaluated. The rapidly evolving/unconstrained matK faired best, whereas remaining genes varied in degrees of contribution to rosid phylogenetics across the lineage's 108 million years evolutionary history. Phylogenetic accuracy was low with the slowly evolving/unconstrained matR despite least amount of homoplasy. Third codon positions contributed the highest amount of parsimony informative sites, resolution and informativeness, but magnitude varied with gene mode of evolution. These findings are in clear contrast with the views that rapidly evolving regions and the 3rd codon position have inevitable negative impact on phylogenetic reconstruction at deep historic level due to accumulation of multiple hits and subsequent elevation in homoplasy and saturation. Relaxed evolutionary constraint in rapidly evolving genes distributes substitutions across codon positions, an evolutionary mode expected to reduce the frequency of multiple hits. These findings should be tested at deeper evolutionary histories.

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.

Land plant evolutionary timeline: gene effects are secondary to fossil constraints in relaxed clock estimation of age and substitution rates

PREMISE OF THE STUDY

Land plants play an essential role in the evolution of terrestrial life. Their time of origin and diversification is fundamental to understanding the evolution of life on land. We investigated the timing and the rate of molecular evolution of land plants, evaluating the effects of different types of molecular data, including temporal information from fossils, and using different molecular clock methods. •

METHODS

Ages and absolute rates were estimated independently with two substitutionally different data sets: a highly conserved 4-gene data set and matK, a fast-evolving gene. The vascular plant backbone and the crown nodes of all major lineages were calibrated with fossil-derived ages. Dates and absolute rates were estimated while including or excluding the calibrations and using two relaxed clocks that differ in their implementation of temporal autocorrelation. •

KEY RESULTS

Land plants diverged from streptophyte alga 912 (870-962) million years ago (Mya) but diversified into living lineages 475 (471-480) Mya. Ages estimated for all major land-plant lineages agree with their fossil record, except for angiosperms. Different genes estimated very similar ages and correlated absolute rates across the tree. Excluding calibrations resulted in the greatest age differences. Different relaxed clocks provided similar ages, but different and uncorrelated absolute rates. •

CONCLUSIONS

Whole-genome rate accelerations or decelerations may underlie the similar ages and correlated absolute rates estimated with different genes. We suggest that pronounced substitution rate changes around the angiosperm crown node may represent a challenge for relaxed clocks to model adequately.

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.

Molecular Evolution of Plant AAP and LHT Amino Acid Transporters

Nitrogen is an essential mineral nutrient and it is often transported within living organisms in its reduced form, as amino acids. Transport of amino acids across cellular membranes requires proteins, and here we report the phylogenetic analysis across taxa of two amino acid transporter families, the amino acid permeases (AAPs) and the lysine-histidine-like transporters (LHTs). We found that the two transporter families form two distinct groups in plants supporting the concept that both are essential. AAP transporters seem to be restricted to land plants. They were found in Selaginella moellendorffii and Physcomitrella patens but not in Chlorophyte, Charophyte, or Rhodophyte algae. AAPs were strongly represented in vascular plants, consistent with their major function in phloem (vascular tissue) loading of amino acids for sink nitrogen supply. LHTs on the other hand appeared prior to land plants. LHTs were not found in chlorophyte algae Chlamydomonas reinhardtii and Volvox carterii. However, the characean alga Klebsormidium flaccidum encodes KfLHT13 and phylogenetic analysis indicates that it is basal to land plant LHTs. This is consistent with the hypothesis that characean algae are ancestral to land plants. LHTs were also found in both S. moellendorffii and P. patens as well as in monocots and eudicots. To date, AAPs and LHTs have mainly been characterized in Arabidopsis (eudicots) and these studies provide clues to the functions of the newly identified homologs.

An ILP solution for the gene duplication problem

Background

The gene duplication (GD) problem seeks a species tree that implies the fewest gene duplication events across a given collection of gene trees. Solving this problem makes it possible to use large gene families with complex histories of duplication and loss to infer phylogenetic trees. However, the GD problem is NP-hard, and therefore, most analyses use heuristics that lack any performance guarantee.

Results

We describe the first integer linear programming (ILP) formulation to solve instances of the gene duplication problem exactly. With simulations, we demonstrate that the ILP solution can solve problem instances with up to 14 taxa. Furthermore, we apply the new ILP solution to solve the gene duplication problem for the seed plant phylogeny using a 12-taxon, 6, 084-gene data set. The unique, optimal solution, which places Gnetales sister to the conifers, represents a new, large-scale genomic perspective on one of the most puzzling questions in plant systematics.

Conclusions

Although the GD problem is NP-hard, our novel ILP solution for it can solve instances with data sets consisting of as many as 14 taxa and 1, 000 genes in a few hours. These are the largest instances that have been solved to optimally to date. Thus, this work can provide large-scale genomic perspectives on phylogenetic questions that previously could only be addressed by heuristic estimates.

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.

Nectar and pollination drops: how different are they?

BACKGROUND

Pollination drops and nectars (floral nectars) are secretions related to plant reproduction. The pollination drop is the landing site for the majority of gymnosperm pollen, whereas nectar of angiosperm flowers represents a common nutritional resource for a large variety of pollinators. Extrafloral nectars also are known from all vascular plants, although among the gymnosperms they are restricted to the Gnetales. Extrafloral nectars are not generally involved in reproduction but serve as 'reward' for ants defending plants against herbivores (indirect defence).

SCOPE

Although very different in their task, nectars and pollination drops share some features, e.g. basic chemical composition and eventual consumption by animals. This has led some authors to call these secretions collectively nectar. Modern techniques that permit chemical analysis and protein characterization have very recently added important information about these sugary secretions that appear to be much more than a 'reward' for pollinating (floral nectar) and defending animals (extrafloral nectar) or a landing site for pollen (pollination drop).

CONCLUSIONS

Nectar and pollination drops contain sugars as the main components, but the total concentration and the relative proportions are different. They also contain amino acids, of which proline is frequently the most abundant. Proteomic studies have revealed the presence of common functional classes of proteins such as invertases and defence-related proteins in nectar (floral and extrafloral) and pollination drops. Invertases allow for dynamic rearrangement of sugar composition following secretion. Defence-related proteins provide protection from invasion by fungi and bacteria. Currently, only few species have been studied in any depth. The chemical composition of the pollination drop must be investigated in a larger number of species if eventual phylogenetic relationships are to be revealed. Much more information can be provided from further proteomic studies of both nectar and pollination drop that will contribute to the study of plant reproduction and evolution.

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.

The complete plastid genome sequence of Welwitschia mirabilis: an unusually compact plastome with accelerated divergence rates

Background

Welwitschia mirabilis is the only extant member of the family Welwitschiaceae, one of three lineages of gnetophytes, an enigmatic group of gymnosperms variously allied with flowering plants or conifers. Limited sequence data and rapid divergence rates have precluded consensus on the evolutionary placement of gnetophytes based on molecular characters. Here we report on the first complete gnetophyte chloroplast genome sequence, from Welwitschia mirabilis, as well as analyses on divergence rates of protein-coding genes, comparisons of gene content and order, and phylogenetic implications.

Results

The chloroplast genome of Welwitschia mirabilis [GenBank: EU342371] is comprised of 119,726 base pairs and exhibits large and small single copy regions and two copies of the large inverted repeat (IR). Only 101 unique gene species are encoded. The Welwitschia plastome is the most compact photosynthetic land plant plastome sequenced to date; 66% of the sequence codes for product. The genome also exhibits a slightly expanded IR, a minimum of 9 inversions that modify gene order, and 19 genes that are lost or present as pseudogenes. Phylogenetic analyses, including one representative of each extant seed plant lineage and based on 57 concatenated protein-coding sequences, place Welwitschia at the base of all seed plants (distance, maximum parsimony) or as the sister to Pinus (the only conifer representative) in a monophyletic gymnosperm clade (maximum likelihood, bayesian). Relative rate tests on these gene sequences show the Welwitschia sequences to be evolving at faster rates than other seed plants. For these genes individually, a comparison of average pairwise distances indicates that relative divergence in Welwitschia ranges from amounts about equal to other seed plants to amounts almost three times greater than the average for non-gnetophyte seed plants.

Conclusion

Although the basic organization of the Welwitschia plastome is typical, its compactness, gene content and high nucleotide divergence rates are atypical. The current lack of additional conifer plastome sequences precludes any discrimination between the gnetifer and gnepine hypotheses of seed plant relationships. However, both phylogenetic analyses and shared genome features identified here are consistent with either of the hypotheses that link gnetophytes with conifers, but are inconsistent with the anthophyte hypothesis.

Widespread Genealogical Nonmonophyly in Species of Pinus Subgenus Strobus

Phylogenetic relationships among Pinus species from subgenus Strobus remain unresolved despite combined efforts based on nrITS and cpDNA. To provide greater resolution among these taxa, a 900-bp intron from a late embryogenesis abundant (LEA)-like gene (IFG8612)was sequenced from 39 pine species, with two or more alleles representing 33 species. Nineteen of 33 species exhibited allelic nonmonphyly in the strict consensus tree, and 10 deviated significantly from allelic monophyly based on topology incongruence tests. Intraspecific nucleotide diversity ranged from 0.0 to 0.0211, and analysis of variance shows that nucleotide diversity was strongly associated (P < 0.0001)with the degree of species monophyly. Although species nonmonophyly complicates phylogenetic interpretations, this nuclear locus offers greater topological support than previously observed for cpDNA or nrITS. Lacking evidence for hybridization, recombination, or imperfect taxonomy, we feel that incomplete lineage sorting remains the best explanation for the polymorphisms shared among species. Depending on the species, coalescent expectations indicate that reciprocal monophyly will be more likely than paraphyly in 1.71 to 24.0 x 10(6) years, and that complete genome-wide coalescence in these species may require up to 76.3 x 10(6) years. The absence of allelic coalescence is a severe constraint in the application of phylogenetic methods in Pinus, and taxa sharing similar life history traits with Pinus are likely to show species nonmonophyly using nuclear markers.

Chloroplast genome (cpDNA) of Cycas taitungensis and 56 cp protein-coding genes of Gnetum parvifolium: insights into cpDNA evolution and phylogeny of extant seed plants

Phylogenetic relationships among the 5 groups of extant seed plants are presently unsettled. To reexamine this long-standing debate, we determine the complete chloroplast genome (cpDNA) of Cycas taitungensis and 56 protein-coding genes encoded in the cpDNA of Gnetum parvifolium. The cpDNA of Cycas is a circular molecule of 163,403 bp with 2 typical large inverted repeats (IRs) of 25,074 bp each. We inferred phylogenetic relationships among major seed plant lineages using concatenated 56 protein-coding genes in 37 land plants. Phylogenies, generated by the use of 3 independent methods, provide concordant and robust support for the monophylies of extant seed plants, gymnosperms, and angiosperms. Within the modern gymnosperms are 2 highly supported sister clades: Cycas-Ginkgo and Gnetum-Pinus. This result agrees with both the "gnetifer" and "gnepines" hypotheses. The sister relationships in Cycas-Ginkgo and Gnetum-Pinus clades are further reinforced by cpDNA structural evidence. Branch lengths of Cycas-Ginkgo and Gnetum were consistently the shortest and the longest, respectively, in all separate analyses. However, the Gnetum relative rate test revealed this tendency only for the 3rd codon positions and the transversional sites of the first 2 codon positions. A PsitufA located between psbE and petL genes is here first detected in Anthoceros (a hornwort), cycads, and Ginkgo. We demonstrate that the PsitufA is a footprint descended from the chloroplast tufA of green algae. The duplication of ycf2 genes and their shift into IRs should have taken place at least in the common ancestor of seed plants more than 300 MYA, and the tRNAPro-GGG gene was lost from the angiosperm lineage at least 150 MYA. Additionally, from cpDNA structural comparison, we propose an alternative model for the loss of large IR regions in black pine. More cpDNA data from non-Pinaceae conifers are necessary to justify whether the gnetifer or gnepines hypothesis is valid and to generate solid structural evidence for the monophyly of extant gymnosperms.

Phytochrome-mediated development in land plants: red light sensing evolves to meet the challenges of changing light environments

Phytochromes are photoreceptors that provide plants with circadian, seasonal, and positional information critical for the control of germination, seedling development, shade avoidance, reproduction, dormancy, and sleep movements. Phytochromes are unique among photoreceptors in their capacity to interconvert between a red-absorbing form (absorption maximum of approximately 660 nm) and a far-red absorbing form (absorption maximum of approximately 730 nm), which occur in a dynamic equilibrium within plant cells, corresponding to the proportions of red and far-red energy in ambient light. Because pigments in stems and leaves absorb wavelengths below about 700 nm, this provides plants with an elegant system for detecting their position relative to other plants, with which the plants compete for light. Certain aspects of phytochrome-mediated development outside of flowering plants are strikingly similar to those that have been characterized in Arabidopsis thaliana and other angiosperms. However, early diverging land plants have fewer distinct phytochrome gene lineages, suggesting that both diversification and subfunctionalization have been important in the evolution of the phytochrome gene family. There is evidence that subfunctionalization proceeded by the partitioning among paralogues of photosensory specificity, physiological response modes, and light-regulated gene expression and protein stability. Parallel events of duplication and functional divergence may have coincided with the evolution of canopy shade and the increasing complexity of the light environment. Within angiosperms, patterns of functional divergence are clade-specific and the roles of phytochromes in A. thaliana change across environments, attesting to the evolutionary flexibility and contemporaneous plasticity of phytochrome signalling in the control of development.

Dating dispersal and radiation in the gymnosperm Gnetum (Gnetales)--clock calibration when outgroup relationships are uncertain

Most implementations of molecular clocks require resolved topologies. However, one of the Bayesian relaxed clock approaches accepts input topologies that include polytomies. We explored the effects of resolved and polytomous input topologies in a rate-heterogeneous sequence data set for Gnetum, a member of the seed plant lineage Gnetales. Gnetum has 10 species in South America, 1 in tropical West Africa, and 20 to 25 in tropical Asia, and explanations for the ages of these disjunctions involve long-distance dispersal and/or the breakup of Gondwana. To resolve relationships within Gnetum, we sequenced most of its species for six loci from the chloroplast (rbcL, matK, and the trnT-trnF region), the nucleus (rITS/5.8S and the LEAFY gene second intron), and the mitochondrion (nad1 gene second intron). Because Gnetum has no fossil record, we relied on fossils from other Gnetales and from the seed plant lineages conifers, Ginkgo, cycads, and angiosperms to constrain a molecular clock and obtain absolute times for within-Gnetum divergence events. Relationships among Gnetales and the other seed plant lineages are still unresolved, and we therefore used differently resolved topologies, including one that contained a basal polytomy among gymnosperms. For a small set of Gnetales exemplars (n = 13) in which rbcL and matK satisfied the clock assumption, we also obtained time estimates from a strict clock, calibrated with one outgroup fossil. The changing hierarchical relationships among seed plants (and accordingly changing placements of distant fossils) resulted in small changes of within-Gnetum estimates because topologically closest constraints overrode more distant constraints. Regardless of the seed plant topology assumed, relaxed clock estimates suggest that the extant clades of Gnetum began diverging from each other during the Upper Oligocene. Strict clock estimates imply a mid-Miocene divergence. These estimates, together with the phylogeny for Gnetum from the six combined data sets, imply that the single African species of Gnetum is not a remnant of a once Gondwanan distribution. Miocene and Pliocene range expansions are inferred for the Asian subclades of Gnetum, which stem from an ancestor that arrived from Africa. These findings fit with seed dispersal by water in several species of Gnetum, morphological similarities among apparently young species, and incomplete concerted evolution in the nuclear ITS region.

Evolution of host specialization in the Adelgidae (Insecta: Hemiptera) inferred from molecular phylogenetics

The Adelgidae form a small group of insects in the Aphidoidea. They are cyclically parthenogenetic with host alternating, multiple-generation complex life cycles and are restricted to certain host genera in the Pinaceae. Species that host alternate always have Picea as the primary host where sexual reproduction and gall formation occur, and another genus in the Pinaceae as the secondary host where a series of parthenogenetic generations are produced. Other species that do not host alternate complete their entire life cycle on one host and only reproduce parthenogenetically. We studied relationships within Adelgidae using DNA sequences from the mitochondrial COI, COII, and cytb genes, and the nuclear EF1alpha gene. Analysis of the combined data resulted in a well-resolved phylogeny in which the major adelgid clades correspond neatly to their association with secondary host genera. Specialization on each secondary host genus occurred only once and was followed by diversification on the host genus. Molecular dating of divergence times in the Adelgidae suggest that diversification among host genera occurred in the Late Cretaceous and Early Tertiary when the Pinaceae genera were diverging. It is not clear, however, whether the Adelgidae and Pinaceae co-diversified because the relationships among the Pinaceae genera are not fully resolved. We discuss implications for adelgid taxonomy, life cycle evolution, and evolution of the interaction between adelgids and their host plants.

Former diversity of Ephedra (Gnetales): evidence from Early Cretaceous seeds from Portugal and North America

BACKGROUND AND AIMS

The extant species of the seed plant group Gnetales (Ephedra, Gnetum and Welwitschia) have been considered a remnant of a much greater, now extinct, diversity due to the pronounced differences in form and ecology among the genera. Until recently, this hypothesis has not been supported by evidence from the fossil record. This paper adds to the expanding information on Gnetales from the Early Cretaceous and describes coalified seeds from Barremian-Albian localities in Portugal and USA.

METHODS

The fossils were extracted from sediment samples by sieving in water. Adhering mineral matrix was removed by chemical treatment. Seeds were investigated using light and scanning electron microscopy. Morphology and anatomy of the seeds were documented and compared with those of extant species.

KEY RESULTS

The fossils share characters with extant Ephedra, for example papillae on the inner surface of the seed envelope and in situ polyplicate pollen grains that shed the exine during germination. They differ from extant Ephedra seeds in morphological and anatomical details as well as in their smaller size. Two new species of Ephedra are described together with one species assigned to a new genus of Gnetales. Other Ephedra-like seeds, for which pollen and critical morphological details are currently unknown, are also present in the samples.

CONCLUSIONS

These Cretaceous seeds document that key reproductive characters and pollen germination processes have remained unchanged within Ephedra for about 120 million years or more. There is sufficient variety in details of morphology to suggest that a diversity of Ephedra and Ephedra-like species were present in the Early Cretaceous flora. Their presence in Portugal and eastern North America indicates that they were widespread on the Laurasian continent. The fossil seeds are similar to seeds of Erdtmanithecales and this supports the previously suggested relationship between Erdtmanithecales and Gnetales.

Seed plant phylogeny: gnetophytes are derived conifers and a sister group to Pinaceae

The phylogenetic position of gnetophytes has long been controversial. We sequenced parts of the genes coding for the largest subunit of nuclear RNA polymerase I, II, and III and combined these sequences with those of four chloroplast genes, two mitochondrial genes, and 18S rRNA genes to address this issue. Both maximum likelihood and maximum parsimony analyses of the sites not affected by high substitution levels strongly support a phylogeny where gymnosperms and angiosperms are monophyletic, where cycads are at the base of gymnosperm tree and are followed by ginkgos, and where gnetophytes are grouped within conifers as the sister group of pines. The evolution of several morphological and molecular characters of gnetophytes and conifers will therefore need to be reinterpreted.

EST analysis in Ginkgo biloba: an assessment of conserved developmental regulators and gymnosperm specific genes

Background

Ginkgo biloba L. is the only surviving member of one of the oldest living seed plant groups with medicinal, spiritual and horticultural importance worldwide. As an evolutionary relic, it displays many characters found in the early, extinct seed plants and extant cycads. To establish a molecular base to understand the evolution of seeds and pollen, we created a cDNA library and EST dataset from the reproductive structures of male (microsporangiate), female (megasporangiate), and vegetative organs (leaves) of Ginkgo biloba.

Results

RNA from newly emerged male and female reproductive organs and immature leaves was used to create three distinct cDNA libraries from which 6,434 ESTs were generated. These 6,434 ESTs from Ginkgo biloba were clustered into 3,830 unigenes. A comparison of our Ginkgo unigene set against the fully annotated genomes of rice and Arabidopsis, and all available ESTs in Genbank revealed that 256 Ginkgo unigenes match only genes among the gymnosperms and non-seed plants – many with multiple matches to genes in non-angiosperm plants. Conversely, another group of unigenes in Gingko had highly significant homology to transcription factors in angiosperms involved in development, including MADS box genes as well as post-transcriptional regulators. Several of the conserved developmental genes found in Ginkgo had top BLAST homology to cycad genes. We also note here the presence of ESTs in G. biloba similar to genes that to date have only been found in gymnosperms and an additional 22 Ginkgo genes common only to genes from cycads.

Conclusion

Our analysis of an EST dataset from G. biloba revealed genes potentially unique to gymnosperms. Many of these genes showed homology to fully sequenced clones from our cycad EST dataset found in common only with gymnosperms. Other Ginkgo ESTs are similar to developmental regulators in higher plants. This work sets the stage for future studies on Ginkgo to better understand seed and pollen evolution, and to resolve the ambiguous phylogenetic relationship of G. biloba among the gymnosperms.

EST analysis in Ginkgo biloba: an assessment of conserved developmental regulators and gymnosperm specific genes

BACKGROUND

Ginkgo biloba L. is the only surviving member of one of the oldest living seed plant groups with medicinal, spiritual and horticultural importance worldwide. As an evolutionary relic, it displays many characters found in the early, extinct seed plants and extant cycads. To establish a molecular base to understand the evolution of seeds and pollen, we created a cDNA library and EST dataset from the reproductive structures of male (microsporangiate), female (megasporangiate), and vegetative organs (leaves) of Ginkgo biloba.

RESULTS

RNA from newly emerged male and female reproductive organs and immature leaves was used to create three distinct cDNA libraries from which 6,434 ESTs were generated. These 6,434 ESTs from Ginkgo biloba were clustered into 3,830 unigenes. A comparison of our Ginkgo unigene set against the fully annotated genomes of rice and Arabidopsis, and all available ESTs in Genbank revealed that 256 Ginkgo unigenes match only genes among the gymnosperms and non-seed plants--many with multiple matches to genes in non-angiosperm plants. Conversely, another group of unigenes in Gingko had highly significant homology to transcription factors in angiosperms involved in development, including MADS box genes as well as post-transcriptional regulators. Several of the conserved developmental genes found in Ginkgo had top BLAST homology to cycad genes. We also note here the presence of ESTs in G. biloba similar to genes that to date have only been found in gymnosperms and an additional 22 Ginkgo genes common only to genes from cycads.

CONCLUSION

Our analysis of an EST dataset from G. biloba revealed genes potentially unique to gymnosperms. Many of these genes showed homology to fully sequenced clones from our cycad EST dataset found in common only with gymnosperms. Other Ginkgo ESTs are similar to developmental regulators in higher plants. This work sets the stage for future studies on Ginkgo to better understand seed and pollen evolution, and to resolve the ambiguous phylogenetic relationship of G. biloba among the gymnosperms.

Phylogenetic relationships in Ephedra (Ephedraceae) inferred from chloroplast and nuclear DNA sequences

Sequences of the nuclear ribosomal DNA internal transcribed spacer region 1 and the chloroplast-encoded genes maturase K and ribulose-1,5 biphosphate carboxylase large subunit were obtained from species of Ephedra (Ephedraceae) representing the geographic range and morphological diversity of the genus. Phylogenetic analyses of the DNA data indicate that relationships within the genus are better predicted by geographic region of origin than by ovulate cone characters. The sampled species with dry, winged (versus fleshy) ovulate cone bracts or single-seeded cones do not form monophyletic groups and therefore the previous classification systems of Ephedra based on these aspects of bract morphology appear to be largely unnatural. Three groups were identified among the Old World species studied, one comprising European and Mediterranean species and two including only Asian species. The sequence data suggest a possible early divergence of a New World clade of Ephedra from among the Old World groups. The South American species form a distinct clade apparently related to one of two groups of North American species, which accords with a frequent floristic pattern of close relationships between species groups in western South America and southwestern North America.

Missing the Forest for the Trees: Phylogenetic Compression and Its Implications for Inferring Complex Evolutionary Histories

Phylogenetic tree reconstruction is difficult in the presence of lateral gene transfer and other processes generating conflicting signals. We develop a new approach to this problem using ideas borrowed from algorithmic information theory. It selects the hypothesis that simultaneously minimizes the descriptive complexity of the tree(s) plus the data when encoded using those tree(s). In practice this is the hypothesis that can compress the data the most. We show not only that phylogenetic compression is an efficient method for encoding most phylogenetic data sets and is more efficient than compression schemes designed for single sequences, but also that it provides a clear information theoretic rule for determining when a collection of conflicting trees is a better explanation of the data than a single tree. By casting the parsimony problem in this more general framework, we also conclude that the so-called total-evidence tree--the tree constructed from all the data simultaneously--is not always the most economical explanation of the data.

Trans-species shared polymorphisms at orthologous nuclear gene loci among distant species in the conifer Picea (Pinaceae): implications for the long-term maintenance of genetic diversity in trees

For each of three nuclear gene loci, intraspecific- as well as trans-specific shared polymorphisms were detected in DNA among three distantly related species in the genus Picea. Few fixed interspecific polymorphisms were observed. Allele genealogies did not match species phylogenies, and species lineages were not reciprocally monophyletic. Based on molecular clocks and morphological evidence from the fossil record, the divergence time between species was estimated at 13-20 million years (my), and a mutation rate of 2.23 × 10(-10) to 3.42 × 10(-10) per site per year was estimated. Large historical population sizes in excess of 100 000 were inferred, which would have delayed the fixation of polymorphisms. These numbers translated into allele coalescence times in the order of 10 to 18 my, which implies the sharing of polymorphisms since common ancestry. These results suggest that trans-species shared polymorphisms might be frequent at plant nuclear gene loci, leading to high allelic diversity. Such a trend is more likely in trees and plants characterized by ecological and life-history determinants favoring large population sizes such as an outcrossing mating system, wind pollination, and a dominant position in ecosystem. These polymorphisms also call for caution in estimating congeneric species phylogenies from nuclear gene sequences in such plant groups.

ANGIOSPERM DIVERGENCE TIMES: THE EFFECT OF GENES, CODON POSITIONS, AND TIME CONSTRAINTS

An understanding of the evolution of modern terrestrial ecosystems requires an understanding of the dynamics associated with angiosperm evolution, including the timing of their origin and diversification into their extraordinary present-day diversity. Molecular estimates of angiosperm age have varied widely, and many substantially predate the Early Cretaceous fossil appearance of the group. In this study, the effect of different genes, codon positions, and chronological constraints on node ages are examined on divergence time estimates across seed plants, with a special focus on angiosperms. Penalized likelihood was used to estimate divergence times on a phylogenetic hypothesis for seed plants derived from Bayesian analysis, with branch lengths estimated with maximum likelihood. The plastid genes atpB, psaA, psbB, and rbcL were used individually and in combination, using first and second, third, and the three codon positions, including and excluding age constraints on 20 nodes derived from a critical examination of the land-plant fossil record. The optimal level of rate smoothing according to each unconstrained and constrained dataset was obtained with penalized likelihood. Tests for a molecular clock revealed significantly unclocklike rates in all datasets. Addition of fossil constraints resulted in even greater departures from constancy. Consistently with significant deviations from a clock, estimated optimal smoothing values were low, but a strict correlation between rate heterogeneity and optimal smoothing value was not found. Age estimates for nodes across the phylogeny varied, sometimes substantially, with gene and codon position. Nevertheless, estimates based on the four concatenated genes are very similar to the mean of the four individual gene estimates. For any given node, unconstrained age estimates are more variable than constrained estimates and are frequently younger than well-substantiated fossil members of the clade. Constrained estimates of ages of clades are older than unconstrained estimates and oldest fossil representatives, sometimes substantially so. Angiosperm age estimates decreased as rate smoothing increased. Whereas the range of unconstrained angiosperm age estimates spans the fossil age of the clade, the range of constrained estimates is narrower (and older) than the earliest angiosperm fossils. Results unambiguously indicate the relevance of constraints in reducing the variability of ages derived from different partitions of the data and diminishing the effect of the smoothing parameter. Constrained optimizations of divergence times and substitution rates across the phylogeny suggest appreciably different evolutionary dynamics for angiosperms and for gymnosperms. Whereas the gymnosperm crown group originated shortly after the origin of seed plants, a long time elapsed before the origin of crown group angiosperms. Although absolute age estimates of angiosperms and angiosperm clades are older than their earliest fossils, the estimated pace of phylogenetic diversification largely agrees with the rapid appearance of angiosperm lineages in stratigraphic sequences.

On the evolutionary history of Ephedra: Cretaceous fossils and extant molecules

Gnetales comprise three unusual genera of seed plants, Ephedra, Gnetum, and Welwitschia. Their extraordinary morphological diversity suggests that they are survivors of an ancient, more diverse group. Gnetalean antiquity is also supported by fossil data. Dispersed "ephedroid" (polyplicate) pollen first appeared in the Permian >250 million years ago (Myr), and a few megafossils document the presence of gnetalean features in the early Cretaceous. The Cretaceous welwitschioid seedling Cratonia cotyledon dates the split between Gnetum and Welwitschia to before 110 Myr. Ages and character evolution of modern diversity are, however, controversial, and, based on molecular data, it has recently been suggested that Ephedra is very young, only 8-32 Myr. Here, we present data on the evolutionary history of Ephedra. Fossil seeds from Buarcos, Portugal, unequivocally link one type of Cretaceous polyplicate pollen to Ephedra and document that plants with unique characters, including the peculiar naked male gametophyte, were established already in the Early Cretaceous. Clades in our molecular phylogeny of extant species correspond to geographical regions, with African species in a basal grade/clade. The study demonstrates extremely low divergence in both molecular and morphological characters in Ephedra. Features observed in the fossils are present in all major extant clades, showing that modern species have retained unique reproductive characters for >110 million years. A recent origin of modern species of Ephedra would imply that the Cretaceous Ephedra fossils discussed here were members of widespread, now extinct sister lineage(s), and that no morphological innovations characterized the second diversification.

Phylogeny and evolution of ferns (monilophytes) with a focus on the early leptosporangiate divergences

The phylogenetic structure of ferns (= monilophytes) is explored here, with a special focus on the early divergences among leptosporangiate lineages. Despite considerable progress in our understanding of fern relationships, a rigorous and comprehensive analysis of the early leptosporangiate divergences was lacking. Therefore, a data set was designed here to include critical taxa that were not included in earlier studies. More than 5000 bp from the plastid (rbcL, atpB, rps4) and the nuclear (18S rDNA) genomes were sequenced for 62 taxa. Phylogenetic analyses of these data (1) confirm that Osmundaceae are sister to the rest of the leptosporangiates, (2) resolve a diverse set of ferns formerly thought to be a subsequent grade as possibly monophyletic (((Dipteridaceae, Matoniaceae), Gleicheniaceae), Hymenophyllaceae), and (3) place schizaeoid ferns as sister to a large clade of "core leptosporangiates" that includes heterosporous ferns, tree ferns, and polypods. Divergence time estimates for ferns are reported from penalized likelihood analyses of our molecular data, with constraints from a reassessment of the fossil record.

Genome-scale data, angiosperm relationships, and "ending incongruence": a cautionary tale in phylogenetics

As systematists grapple with assembling the Tree of Life, recent studies have encouraged a genomic-scale approach, obtaining DNA sequence data for entire nuclear, plastid or mitochondrial genomes for a few exemplar taxa. Some have proclaimed that this comparative genomic strategy heralds the end of incongruence in phylogeny reconstruction. Although we applaud the use of many genes to resolve phylogenetic patterns, there is a significant caveat. In spite of, or even because of, the abundant data per taxon, whole-genome sequencing for a few exemplars can provide completely resolved and strongly supported, but incorrect, evolutionary reconstructions. We provide a conspicuous example that includes Amborella, the putative sister of all other extant angiosperms, highlighting the limits of phylogenetics when whole genomes are used but taxon sampling is poor.

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.

Fossils and plant phylogeny

Developing a detailed estimate of plant phylogeny is the key first step toward a more sophisticated and particularized understanding of plant evolution. At many levels in the hierarchy of plant life, it will be impossible to develop an adequate understanding of plant phylogeny without taking into account the additional diversity provided by fossil plants. This is especially the case for relatively deep divergences among extant lineages that have a long evolutionary history and in which much of the relevant diversity has been lost by extinction. In such circumstances, attempts to integrate data and interpretations from extant and fossil plants stand the best chance of success. For this to be possible, what will be required is meticulous and thorough descriptions of fossil material, thoughtful and rigorous analysis of characters, and careful comparison of extant and fossil taxa, as a basis for determining their systematic relationships.

Fossil evidence and phylogeny: the age of major angiosperm clades based on mesofossil and macrofossil evidence from Cretaceous deposits

The fossil record has played an important role in the history of evolutionary thought, has aided the determination of key relationships through mosaics, and has allowed an assessment of a number of ecological hypotheses. Nonetheless, expectations that it might accurately and precisely mirror the progression of taxa through time seem optimistic in light of the many factors potentially interfering with uniform preservation. In view of these limitations, attempts to use the fossil record to corroborate phylogenetic hypotheses based on extensive comparisons among extant taxa may be misplaced. Instead we suggest a method-minimum age node mapping-for combining reliable fossil evidence with hypotheses of phylogeny. We use this methodology in conjunction with a phylogeny for angiosperms to assess timing in the history of major angiosperm clades. This method places many clades both with and without fossil records in temporal perspective, reveals discrepancies among clades in propensities for preservation, and raises some interesting questions about angiosperm evolution. By providing a context for understanding the gaps in the angiosperm fossil record this technique lends credibility and support to the remainder of the angiosperm record and to its applications in understanding a variety of aspects of angiosperm history. In effect, this methodology empowers the fossil record.

Long branch attraction, taxon sampling, and the earliest angiosperms: Amborella or monocots?

Background

Numerous studies, using in aggregate some 28 genes, have achieved a consensus in recognizing three groups of plants, including Amborella, as comprising the basal-most grade of all other angiosperms. A major exception is the recent study by Goremykin et al. (2003; Mol. Biol. Evol. 20:1499–1505), whose analyses of 61 genes from 13 sequenced chloroplast genomes of land plants nearly always found 100% support for monocots as the deepest angiosperms relative to Amborella, Calycanthus, and eudicots. We hypothesized that this conflict reflects a misrooting of angiosperms resulting from inadequate taxon sampling, inappropriate phylogenetic methodology, and rapid evolution in the grass lineage used to represent monocots.

Results

We used two main approaches to test this hypothesis. First, we sequenced a large number of chloroplast genes from the monocot Acorus and added these plus previously sequenced Acorus genes to the Goremykin et al. (2003) dataset in order to explore the effects of altered monocot sampling under the same analytical conditions used in their study. With Acorus alone representing monocots, strongly supported Amborella-sister trees were obtained in all maximum likelihood and parsimony analyses, and in some distance-based analyses. Trees with both Acorus and grasses gave either a well-supported Amborella-sister topology or else a highly unlikely topology with 100% support for grasses-sister and paraphyly of monocots (i.e., Acorus sister to "dicots" rather than to grasses). Second, we reanalyzed the Goremykin et al. (2003) dataset focusing on methods designed to account for rate heterogeneity. These analyses supported an Amborella-sister hypothesis, with bootstrap support values often conflicting strongly with cognate analyses performed without allowing for rate heterogeneity. In addition, we carried out a limited set of analyses that included the chloroplast genome of Nymphaea, whose position as a basal angiosperm was also, and very recently, challenged.

Conclusions

These analyses show that Amborella (or Amborella plus Nymphaea), but not monocots, is the sister group of all other angiosperms among this limited set of taxa and that the grasses-sister topology is a long-branch-attraction artifact leading to incorrect rooting of angiosperms. These results highlight the danger of having lots of characters but too few and, especially, molecularly divergent taxa, a situation long recognized as potentially producing strongly misleading molecular trees. They also emphasize the importance in phylogenetic analysis of using appropriate evolutionary models.

The mitochondrial DNA of land plants: peculiarities in phylogenetic perspective

Land plants exhibit a significant evolutionary plasticity in their mitochondrial DNA (mtDNA), which contrasts with the more conservative evolution of their chloroplast genomes. Frequent genomic rearrangements, the incorporation of foreign DNA from the nuclear and chloroplast genomes, an ongoing transfer of genes to the nucleus in recent evolutionary times and the disruption of gene continuity in introns or exons are the hallmarks of plant mtDNA, at least in flowering plants. Peculiarities of gene expression, most notably RNA editing and trans-splicing, are significantly more pronounced in land plant mitochondria than in chloroplasts. At the same time, mtDNA is generally the most slowly evolving of the three plant cell genomes on the sequence level, with unique exceptions in only some plant lineages. The slow sequence evolution and a variable occurrence of introns in plant mtDNA provide an attractive reservoir of phylogenetic information to trace the phylogeny of older land plant clades, which is as yet not fully resolved. This review attempts to summarize the unique aspects of land plant mitochondrial evolution from a phylogenetic perspective.

Amborella not a "basal angiosperm"? Not so fast

The sequence of the plastid genome of Amborella trichopoda, the putative sister to all other extant angiosperms, was recently reported (Molecular Biology and Evolution 20: 1499-1505). Goremykin et al. used sequence data for 61 plastid genes from Amborella and 12 other embryophytes in phylogenetic analyses and concluded that Amborella is not the sister to the remaining flowering plants; the monocots instead occupy this position. The authors attributed their results, which differ substantially from all recent phylogenetic analyses of angiosperms, to the increased character sampling (30 017 nucleotides in their aligned matrix) in their analysis relative to published studies that included fewer genes but more taxa. We hypothesized that the difference in topology is not due to limited character sampling in previous studies but to limited taxon sampling in the analysis by Goremykin et al. To test this, we conducted a series of phylogenetic analyses using a three-gene, 12 (or more)-taxon data set to evaluate the topological effects of (i) including three vs. 61 genes for (nearly) the same set of taxa, (ii) analyzing different codon positions, (iii) substituting representatives of other basal lineages for Amborella, (iv) replacing the grasses used to represent the monocots with other monocots, selected either for their phylogenetic position or randomly, and (v) adding other basal taxa-Nymphaea, Austrobaileya, magnoliids, and monocots-to the 12-taxon data set. Our results demonstrate that the "monocots basal" topology obtained by Goremykin et al. is not due to increased character sampling of the plastid genome; their topology was obtained using only two plastid genes or two plastid genes and one nuclear gene. This topology was also retained when either Nymphaea or Austrobaileya was substituted for Amborella, demonstrating that any of the three basal lineages will attach to Calycanthus for lack of any other close branch. Furthermore, the "monocots basal" topology is not robust to changes in sampling of monocots. Simply adding Oncidium, for example, places Amborella sister to the other angiosperms. Thus, limited taxon sampling, focusing on organisms with complete genome sequences, can lead to artifactual results.

Ontogeny of triovulate cones of Ephedra intermedia and origin of the outer envelope ofovules of Ephedraceae

The origin of the two envelopes of ovules of extant species of Ephedraceae has been discussed previously in different ways based only on investigation of bimerous female cones. To generalize the characteristics of the two envelopes, ontogenetic and metamorphic patterns of the two envelopes of ovules of trimerous female cones of Ephedra intermedia Schrenk et C. A. Meyer were studied using scanning electron microscopy. The outer envelope of the ovules of triovulate female cones of Ephedra intermedia is initiated as two adaxial lateral protuberances. Thus, the outer envelope of extant Ephedraceae may have been derived from three foliar components of the second proximal whorl of the secondary reproductive shoot of the ancestral compound female cone. Two adaxial lateral foliar components likely gave rise to the outer envelope by fusion, while the third abaxial foliar component was reduced. These new findings suggest that the Gnetales derived their female cones from trimerous multiaxial female reproductive organs of their progenitors. Some evolutionary aspects of female cones of Ephedra are also discussed.

Horizontal gene transfer from flowering plants to Gnetum

Although horizontal gene transfer is well documented in microbial genomes, no case has been reported in higher plants. We discovered horizontal transfer of the mitochondrial nad1 intron 2 and adjacent exons b and c from an asterid to Gnetum (Gnetales, gymnosperms). Gnetum has two copies of intron 2, a group II intron, that differ in their exons, nucleotide composition, domain lengths, and structural characteristics. One of the copies, limited to an Asian clade of Gnetum, is almost identical to the homologous locus in angiosperms, and partial sequences of its exons b and c show characteristic substitutions unique to angiosperms. Analyses of 70 seed plant nad1 exons b and c and intron 2 sequences, including representatives of all angiosperm clades, support that this copy originated from a euasterid and was horizontally transferred to Gnetum. Molecular clock dating, using calibrations provided by gnetalean macrofossils, suggests an age of 5 to 2 million years for the Asian clade that received the horizontal transfer.