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

The Systematics and Evolution of Gymnosperms with an Emphasis on a Few Problematic Taxa

Gymnosperms originated in the Middle Devonian and have experienced a long evolutionary history with pulses of speciation and extinction, which resulted in the four morphologically distinct extant groups, i.e., cycads, Ginkgo, conifers and gnetophytes. For over a century, the systematic relationships within the extant gymnosperms have been debated because different authors emphasized different characters. Recent phylogenomic studies of gymnosperms have given a consistent topology, which aligns well with extant gymnosperms classified into three classes, five subclasses, eight orders, and 13 families. Here, we review the historical opinions of systematics of gymnosperms with special reference to several problematic taxa and reconsider the evolution of some key morphological characters previously emphasized by taxonomists within a phylogenomic context. We conclude that (1) cycads contain two families, i.e., the Cycadaceae and the Zamiaceae; (2) Ginkgo is sister to cycads but not to conifers, with the similarities between Ginkgo and conifers being the result of parallel evolution including a monopodial growth pattern, pycnoxylic wood in long shoots, and the compound female cones, and the reproductive similarities between Ginkgo and cycads are either synapomorphic or plesiomorphic, e.g., the boat-shaped pollen, the branched pollen tube, and the flagellate sperms; (3) conifers are paraphyletic with gnetophytes nested within them, thus gnetophytes are derived conifers, and our newly delimited coniferophytes are equivalent to the Pinopsida and include three subclasses, i.e., Pinidae, Gnetidae, and Cupressidae; (4) fleshy cones of conifers originated multiple times, the Podocarpaceae are sister to the Araucariaceae, the Cephalotaxaceae and the Taxaceae comprise a small clade, which is sister to the Cupressaceae; (5) the Cephalotaxaceae are distinct from the Taxaceae, because the former family possesses typical female cones and the fleshy part of the seed is derived from the fleshiness of integument, while the latter family has reduced female cones and preserves no traces of the seed scale complexes.

Assembly and Annotation of Red Spruce (Picea rubens) Chloroplast Genome, Identification of Simple Sequence Repeats, and Phylogenetic Analysis in Picea

We have sequenced the chloroplast genome of red spruce (Picea rubens) for the first time using the single-end, short-reads (44 bp) Illumina sequences, assembled and functionally annotated it, and identified simple sequence repeats (SSRs). The contigs were assembled using SOAPdenovo2 following the retrieval of chloroplast genome sequences using the black spruce (Picea mariana) chloroplast genome as the reference. The assembled genome length was 122,115 bp (gaps included). Comparatively, the P. rubens chloroplast genome reported here may be considered a near-complete draft. Global genome alignment and phylogenetic analysis based on the whole chloroplast genome sequences of Picea rubens and 10 other Picea species revealed high sequence synteny and conservation among 11 Picea species and phylogenetic relationships consistent with their known classical interrelationships and published molecular phylogeny. The P. rubens chloroplast genome sequence showed the highest similarity with that of P. mariana and the lowest with that of P. sitchensis. We have annotated 107 genes including 69 protein-coding genes, 28 tRNAs, 4 rRNAs, few pseudogenes, identified 42 SSRs, and successfully designed primers for 26 SSRs. Mononucleotide A/T repeats were the most common followed by dinucleotide AT repeats. A similar pattern of microsatellite repeats occurrence was found in the chloroplast genomes of 11 Picea species.

Recent advances on phylogenomics of gymnosperms and a new classification

Living gymnosperms comprise four major groups: cycads, Ginkgo, conifers, and gnetophytes. Relationships among/within these lineages have not been fully resolved. Next generation sequencing has made available a large number of sequences, including both plastomes and single-copy nuclear genes, for reconstruction of solid phylogenetic trees. Recent advances in gymnosperm phylogenomic studies have updated our knowledge of gymnosperm systematics. Here, we review major advances of gymnosperm phylogeny over the past 10 years and propose an updated classification of extant gymnosperms. This new classification includes three classes (Cycadopsida, Ginkgoopsida, and Pinopsida), five subclasses (Cycadidae, Ginkgoidae, Cupressidae, Pinidae, and Gnetidae), eight orders (Cycadales, Ginkgoales, Araucariales, Cupressales, Pinales, Ephedrales, Gnetales, and Welwitschiales), 13 families, and 86 genera. We also described six new tribes including Acmopyleae Y. Yang, Austrocedreae Y. Yang, Chamaecyparideae Y. Yang, Microcachrydeae Y. Yang, Papuacedreae Y. Yang, and Prumnopityeae Y. Yang, and made 27 new combinations in the genus Sabina.

Comparative Wood Anatomy in Pinaceae with Reference to Its Systematic Position

The wood anatomy of 132 species of the genera Abies, Cathaya, Cedrus, Keteleeria, Larix, Nothotsuga, Picea, Pinus, Pseudolarix, Pseudotsuga and Tsuga was studied to determine the elements that characterise the xylem of each genus and discuss possible groupings by wood anatomy for comparison with clades established by molecular phylogeny. The presence of resin canals and ray tracheids supports the family Pinaceae, although the absence of ray tracheids in Keteleeria and their occasional presence in Abies and Pseudolarix weakens it. Based on wood structure, Pinaceae clearly supports division into two groups, coinciding with molecular phylogeny: Pinoideae (Cathaya-Larix-Picea-Pinus-Pseudotsuga) and Abietoideae (Abies-Cedrus-Keteleeria-Nothotsuga-Pseudolarix-Tsuga). Although differences between genera are slight in Pinoideae, the Abietoideae group presents problems such as the presence of only axial resin canals in Keteleeria and Nothotsuga, absence of ray tracheids in Keteleeria and presence of traumatic radial resin canals in Cedrus. However, other features such as pitted horizontal walls and nodular end walls of ray parenchyma cells, indentures, scarce marginal axial parenchyma and presence of crystals in ray parenchyma strengthen the Abietoideae group.

A new macrofossil ephedroid plant with unusual bract morphology from the Lower Cretaceous Jiufotang Formation of northeastern China

BACKGROUND

The evolution of the Jehol Biota of western Liaoning in China includes three phases, initiation in the Dabeigou phase, radiation in the Yixian phase, and decline in the Jiufotang phase. Numerous ephedroid macrofossils were reported from the Lower Cretaceous Yixian Formation. However, so far none has been found in the younger Jiufotang Formation (ca. 120.3 Ma) of western Liaoning.

RESULTS

Here we report a new species Jianchangia verticillata gen. et sp. nov. with unusual morphology from the Lower Cretaceous of the Jiufotang Formation, Lamadong Village, Jianchang County, Liaoning. This species is the first record of gnetophytes from the Jiufotang Formation. It is similar to other ephedroid species from the Yixian Formation in possessing linear leaves with parallel veins, jointed shoots with swollen nodes and longitudinally furrowed internodes, and ovulate cones possessing two whorls of bracts enclosing two chlamydosperms, but differs from all known species by the ovulate cone having multiple fine linear verticillate bracts.

CONCLUSIONS

This study expands our knowledge about the diversity of early gnetophytes in the Lower Cretaceous, and demonstrates the lineage continuity of gnetophytes from the Yixian Formation to the younger Jiufotang Formation.

The complete mitochondrial genome of Taxus cuspidata (Taxaceae): eight protein-coding genes have transferred to the nuclear genome

BACKGROUND

Gymnosperms represent five of the six lineages of seed plants. However, most sequenced plant mitochondrial genomes (mitogenomes) have been generated for angiosperms, whereas mitogenomic sequences have been generated for only six gymnosperms. In particular, complete mitogenomes are available for all major seed plant lineages except Conifer II (non-Pinaceae conifers or Cupressophyta), an important lineage including six families, which impedes a comprehensive understanding of the mitogenomic diversity and evolution in gymnosperms.

RESULTS

Here, we report the complete mitogenome of Taxus cuspidata in Conifer II. In comparison with previously released gymnosperm mitogenomes, we found that the mitogenomes of Taxus and Welwitschia have lost many genes individually, whereas all genes were identified in the mitogenomes of Cycas, Ginkgo and Pinaceae. Multiple tRNA genes and introns also have been lost in some lineages of gymnosperms, similar to the pattern observed in angiosperms. In general, gene clusters could be less conserved in gymnosperms than in angiosperms. Moreover, fewer RNA editing sites were identified in the Taxus and Welwitschia mitogenomes than in other mitogenomes, which could be correlated with fewer introns and frequent gene losses in these two species.

CONCLUSIONS

We have sequenced the Taxus cuspidata mitogenome, and compared it with mitogenomes from the other four gymnosperm lineages. The results revealed the diversity in size, structure, gene and intron contents, foreign sequences, and mutation rates of gymnosperm mitogenomes, which are different from angiosperm mitogenomes.

Phylogenetically distant group of terpene synthases participates in cadinene and cedrane-type sesquiterpenes accumulation in Taiwania cryptomerioides

Along with the species evolution, plants have evolved ways to produce a different collection of terpenoids to accommodate its biotic and abiotic environment, and terpene synthase (TPS) is one of the major contributors to various terpene compounds. The timber of a monotypic and relictual conifer species of Cupressace, Taiwania cryptomerioides, has excellent durability, and one of the essential factors for Taiwania to resist decay and insect pests is sesquiterpene. Compared to other conifers, Taiwania has much higher abundance of cadinene-type sesquiterpenes, and the presence of cedrene-type sesquiterpenes. To understand sesquiterpene biosynthesis in Taiwania, we functionally characterized 10 T. cryptomerioides TPSs (TcTPSs) in vivo or in planta, which could catalyze sesquiterpene formation and potentially are involved in biosynthesis of diverse sesquiterpenoids in Taiwania. The distant phylogenetic relationship and the intron loss event of TcTPSs correlate to the differentiation of chemical profile Taiwania compared to other conifers. Furthermore, we identified TcTPS3 and TcTPS12 as δ-cadinene synthase, and TcTPS6 as cedrol synthase, which demonstrates the important contributions of dynamic evolution in TPSs to the chemical diversity in plants. Combining with functional characterization and comparison of catalytic residues, we conclude at least three catalytic routes for sesquiterpene biosynthesis in this species, and the skeleton diversity has been expended in T. cryptomeriodes.

An overview of extant conifer evolution from the perspective of the fossil record

PREMISE OF THE STUDY

Conifers are an important living seed plant lineage with an extensive fossil record spanning more than 300 million years. The group therefore provides an excellent opportunity to explore congruence and conflict between dated molecular phylogenies and the fossil record.

METHODS

We surveyed the current state of knowledge in conifer phylogenetics to present a new time-calibrated molecular tree that samples ~90% of extant species diversity. We compared phylogenetic relationships and estimated divergence ages in this new phylogeny with the paleobotanical record, focusing on clades that are species-rich and well known from fossils.

KEY RESULTS

Molecular topologies and estimated divergence ages largely agree with the fossil record in Cupressaceae, conflict with it in Araucariaceae, and are ambiguous in Pinaceae and Podocarpaceae. Molecular phylogenies provide insights into some fundamental questions in conifer evolution, such as the origin of their seed cones, but using them to reconstruct the evolutionary history of specific traits can be challenging.

CONCLUSIONS

Molecular phylogenies are useful for answering deep questions in conifer evolution if they depend on understanding relationships among extant lineages. Because of extinction, however, molecular datasets poorly sample diversity from periods much earlier than the Late Cretaceous. This fundamentally limits their utility for understanding deep patterns of character evolution and resolving the overall pattern of conifer phylogeny.

Macrofossil evidence unveiling evolution of male cones in Ephedraceae (Gnetidae)

Background

Male cones of modern Ephedraceae are compound and compact. No fossil evidence has so far been found to support an origin of the compact compound male cone from a hypothetical loosely-arranged shoot system.

Results

Here we describe a new macrofossil taxon, Eamesia chinensis Yang, Lin, Ferguson et Wang, gen. et sp. nov., from the Early Cretaceous of western Liaoning, northeastern China. It was an ephedroid shrub bearing male spikes terminal to twigs, but differs from modern Ephedraceae by its loosely-arranged male cones, the axillary male shoot consisting of an elongated synangiophore on which leaf-like foliar organs were inserted, and four sessile synangia terminal to the apex.

Conclusions

The morphology of this fossil suggests that the modern compact male cone of Ephedra was indeed derived from a once loosely-arranged shoot system, and the male reproductive unit originated from a once elongated axillary male shoot. This new fossil species thus provides a transitional link from the hypothetical ancestral shoot system to the modern compact morphology. Changes of habitat from closed humid forests to open dry deserts and shifts of the pollination syndrome may have acted as the driving forces behind this morphological evolution.

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.

Localized Retroprocessing as a Model of Intron Loss in the Plant Mitochondrial Genome

Loss of introns in plant mitochondrial genes is commonly explained by retroprocessing. Under this model, an mRNA is reverse transcribed and integrated back into the genome, simultaneously affecting the contents of introns and edited sites. To evaluate the extent to which retroprocessing explains intron loss, we analyzed patterns of intron content and predicted RNA editing for whole mitochondrial genomes of 30 species in the monocot order Alismatales. In this group, we found an unusually high degree of variation in the intron content, even expanding the hitherto known variation among angiosperms. Some species have lost some two-third of the cis-spliced introns. We found a strong correlation between intron content and editing frequency, and detected 27 events in which intron loss is consistent with the presence of nucleotides in an edited state, supporting retroprocessing. However, we also detected seven cases of intron loss not readily being explained by retroprocession. Our analyses are also not consistent with the entire length of a fully processed cDNA copy being integrated into the genome, but instead indicate that retroprocessing usually occurs for only part of the gene. In some cases, several rounds of retroprocessing may explain intron loss in genes completely devoid of introns. A number of taxa retroprocessing seem to be very common and a possibly ongoing process. It affects the entire mitochondrial genome.

Evidence of intense chromosomal shuffling during conifer evolution

Although recent advances have been gained on genome evolution in angiosperm lineages, virtually nothing is known about karyotype evolution in the other group of seed plants, the gymnosperms. Here, we used high-density gene-based linkage mapping to compare the karyotype structure of two families of conifers (the most abundant group of gymnosperms) separated around 290 Ma: Pinaceae and Cupressaceae. We propose for the first time a model based on the fusion of 20 ancestral chromosomal blocks that may have shaped the modern karyotpes of Pinaceae (with n = 12) and Cupressaceae (with n = 11). The considerable difference in modern genome organization between these two lineages contrasts strongly with the remarkable level of synteny already reported within the Pinaceae. It also suggests a convergent evolutionary mechanism of chromosomal block shuffling that has shaped the genomes of the spermatophytes.

MADS goes genomic in conifers: towards determining the ancestral set of MADS-box genes in seed plants

BACKGROUND AND AIMS

MADS-box genes comprise a gene family coding for transcription factors. This gene family expanded greatly during land plant evolution such that the number of MADS-box genes ranges from one or two in green algae to around 100 in angiosperms. Given the crucial functions of MADS-box genes for nearly all aspects of plant development, the expansion of this gene family probably contributed to the increasing complexity of plants. However, the expansion of MADS-box genes during one important step of land plant evolution, namely the origin of seed plants, remains poorly understood due to the previous lack of whole-genome data for gymnosperms.

METHODS

The newly available genome sequences of Picea abies, Picea glauca and Pinus taeda were used to identify the complete set of MADS-box genes in these conifers. In addition, MADS-box genes were identified in the growing number of transcriptomes available for gymnosperms. With these datasets, phylogenies were constructed to determine the ancestral set of MADS-box genes of seed plants and to infer the ancestral functions of these genes.

KEY RESULTS

Type I MADS-box genes are under-represented in gymnosperms and only a minimum of two Type I MADS-box genes have been present in the most recent common ancestor (MRCA) of seed plants. In contrast, a large number of Type II MADS-box genes were found in gymnosperms. The MRCA of extant seed plants probably possessed at least 11-14 Type II MADS-box genes. In gymnosperms two duplications of Type II MADS-box genes were found, such that the MRCA of extant gymnosperms had at least 14-16 Type II MADS-box genes.

CONCLUSIONS

The implied ancestral set of MADS-box genes for seed plants shows simplicity for Type I MADS-box genes and remarkable complexity for Type II MADS-box genes in terms of phylogeny and putative functions. The analysis of transcriptome data reveals that gymnosperm MADS-box genes are expressed in a great variety of tissues, indicating diverse roles of MADS-box genes for the development of gymnosperms. This study is the first that provides a comprehensive overview of MADS-box genes in conifers and thus will provide a framework for future work on MADS-box genes in seed plants.

Phylogenomics and coalescent analyses resolve extant seed plant relationships

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

A new phylogeny for the genus Picea from plastid, mitochondrial, and nuclear sequences

Studies over the past ten years have shown that the crown groups of most conifer genera are only about 15-25 Ma old. The genus Picea (spruces, Pinaceae), with around 35 species, appears to be no exception. In addition, molecular studies of co-existing spruce species have demonstrated frequent introgression. Perhaps not surprisingly therefore previous phylogenetic studies of species relationships in Picea, based mostly on plastid sequences, suffered from poor statistical support. We therefore generated mitochondrial, nuclear, and further plastid DNA sequences from carefully sourced material, striking a balance between alignability with outgroups and phylogenetic signal content. Motif duplications in mitochondrial introns were treated as characters in a stochastic Dollo model; molecular clock models were calibrated with fossils; and ancestral ranges were inferred under maximum likelihood. In agreement with previous findings, Picea diverged from its sister clade 180 million years ago (Ma), and the most recent common ancestor of today's spruces dates to 28 Ma. Different from previous analyses though, we find a large Asian clade, an American clade, and a Eurasian clade. Two expansions occurred from Asia to North America and several between Asia and Europe. Chinese P. brachytyla, American P. engelmannii, and Norway spruce, P. abies, are not monophyletic, and North America has ten, not eight species. Divergence times imply that Pleistocene refugia are unlikely to be the full explanation for the relationships between the European species and their East Asian relatives. Thus, northern Norway spruce may be part of an Asian species complex that diverged from the southern Norway spruce lineage in the Upper Miocene, some 6 Ma, which can explain the deep genetic gap noted in phylogeographic studies of Norway spruce. The large effective population sizes of spruces, and incomplete lineage sorting during speciation, mean that the interspecific relationships within each of the geographic clades require further studies, especially based on genomic information and population genetic data.

Natural stilbenoids: distribution in the plant kingdom and chemotaxonomic interest in Vitaceae

Stilbenoids, a family of polyphenols known for the complexity of their structure and for their diverse biological activities, occur with a limited but heterogeneous distribution in the plant kingdom. The most prominent stilbene containing plant family, the Vitaceae, represented by the famous wine producing grape vines Vitis vinifera L., is one of the richest sources of novel stilbenes currently known, together with other families, such as Dipterocarpaceae, Gnetaceae and Fabaceae. This review focuses on the distribution of stilbenes and 2-arylbenzofuran derivatives in the plant kingdom, the chemical structure of stilbenes in the Vitaceae family and their taxonomic implication.

Genes and processed paralogs co-exist in plant mitochondria

RNA-mediated gene duplication has been proposed to create processed paralogs in the plant mitochondrial genome. A processed paralog may retain signatures left by the maturation process of its RNA precursor, such as intron removal and no need of RNA editing. Whereas it is well documented that an RNA intermediary is involved in the transfer of mitochondrial genes to the nucleus, no direct evidence exists for insertion of processed paralogs in the mitochondria (i.e., processed and un-processed genes have never been found simultaneously in the mitochondrial genome). In this study, we sequenced a region of the mitochondrial gene nad1, and identified a number of taxa were two different copies of the region co-occur in the mitochondria. The two nad1 paralogs differed in their (a) presence or absence of a group II intron, and (b) number of edited sites. Thus, this work provides the first evidence of co-existence of processed paralogs and their precursors within the plant mitochondrial genome. In addition, mapping the presence/absence of the paralogs provides indirect evidence of RNA-mediated gene duplication as an essential process shaping the mitochondrial genome in plants.

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.

Are substitution rates and RNA editing correlated?

Background

RNA editing is a post-transcriptional process that, in seed plants, involves a cytosine to uracil change in messenger RNA, causing the translated protein to differ from that predicted by the DNA sequence. RNA editing occurs extensively in plant mitochondria, but large differences in editing frequencies are found in some groups. The underlying processes responsible for the distribution of edited sites are largely unknown, but gene function, substitution rate, and gene conversion have been proposed to influence editing frequencies.

Results

We studied five mitochondrial genes in the monocot order Alismatales, all showing marked differences in editing frequencies among taxa. A general tendency to lose edited sites was observed in all taxa, but this tendency was particularly strong in two clades, with most of the edited sites lost in parallel in two different areas of the phylogeny. This pattern is observed in at least four of the five genes analyzed. Except in the groups that show an unusually low editing frequency, the rate of C-to-T changes in edited sites was not significantly higher that in non-edited 3^rd codon positions. This may indicate that selection is not actively removing edited sites in nine of the 12 families of the core Alismatales. In all genes but ccmB, a significant correlation was found between frequency of change in edited sites and synonymous substitution rate. In general, taxa with higher substitution rates tend to have fewer edited sites, as indicated by the phylogenetically independent correlation analyses. The elimination of edited sites in groups that lack or have reduced levels of editing could be a result of gene conversion involving a cDNA copy (retroprocessing). If so, this phenomenon could be relatively common in the Alismatales, and may have affected some groups recurrently. Indirect evidence of retroprocessing without a necessary correlation with substitution rate was found mostly in families Alismataceae and Hydrocharitaceae (e.g., groups that suffered a rapid elimination of all their edited sites, without a change in substitution rate).

Conclusions

The effects of substitution rate, selection, and/or gene conversion on the dynamics of edited sites in plant mitochondria remain poorly understood. Although we found an inverse correlation between substitution rate and editing frequency, this correlation is partially obscured by gene retroprocessing in lineages that have lost most of their edited sites. The presence of processed paralogs in plant mitochondria deserves further study, since most evidence of their occurrence is circumstantial.

Comparative chloroplast genomics reveals the evolution of Pinaceae genera and subfamilies

As the largest and the basal-most family of conifers, Pinaceae provides key insights into the evolutionary history of conifers. We present comparative chloroplast genomics and analysis of concatenated 49 chloroplast protein-coding genes common to 19 gymnosperms, including 15 species from 8 Pinaceous genera, to address the long-standing controversy about Pinaceae phylogeny. The complete cpDNAs of Cathaya argyrophylla and Cedrus deodara (Abitoideae) and draft cpDNAs of Larix decidua, Picea morrisonicola, and Pseudotsuga wilsoniana are reported. We found 21- and 42-kb inversions in congeneric species and different populations of Pinaceous species, which indicates that structural polymorphics may be common and ancient in Pinaceae. Our phylogenetic analyses reveal that Cedrus is clustered with Abies-Keteleeria rather than the basal-most genus of Pinaceae and that Cathaya is closer to Pinus than to Picea or Larix-Pseudotsuga. Topology and structural change tests and indel-distribution comparisons lend further evidence to our phylogenetic finding. Our molecular datings suggest that Pinaceae first evolved during Early Jurassic, and diversification of Pinaceous subfamilies and genera took place during Mid-Jurassic and Lower Cretaceous, respectively. Using different maximum-likelihood divergences as thresholds, we conclude that 2 (Abietoideae and Larix-Pseudotsuga-Piceae-Cathaya-Pinus), 4 (Cedrus, non-Cedrus Abietoideae, Larix-Pseudotsuga, and Piceae-Cathaya-Pinus), or 5 (Cedrus, non-Cedrus Abietoideae, Larix-Pseudotsuga, Picea, and Cathaya-Pinus) groups/subfamilies are more reasonable delimitations for Pinaceae. Specifically, our views on subfamilial classifications differ from previous studies in terms of the rank of Cedrus and with recognition of more than two subfamilies.

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.

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.

Frequent, phylogenetically local horizontal transfer of the cox1 group I Intron in flowering plant mitochondria

Horizontal gene transfer is surprisingly common among plant mitochondrial genomes. The first well-established case involves a homing group I intron in the mitochondrial cox1 gene shown to have been frequently acquired via horizontal transfer in angiosperms. Here, we report extensive additional sampling of angiosperms, including 85 newly sequenced introns from 30 families. Analysis of all available data leads us to conclude that, among the 640 angiosperms (from 212 families) whose cox1 intron status has been characterized thus far, the intron has been acquired via roughly 70 separate horizontal transfer events. We propose that the intron was originally seeded into angiosperms by a single transfer from fungi, with all subsequent inferred transfers occurring from one angiosperm to another. The pattern of angiosperm-to-angiosperm transfer is biased toward exchanges between plants belonging to the same family. Illegitimate pollination is proposed as one potential factor responsible for this pattern, given that aberrant, cross-species pollination is more likely between close relatives. Other potential factors include shared vectoring agents or common geographic locations. We report the first apparent cases of loss of the cox1 intron; losses are accompanied by retention of the exonic coconversion tract, which is located immediately downstream of the intron and which is a product of the intron's self-insertion mechanism. We discuss the many reasons why the cox1 intron is so frequently and detectably transferred, and rarely lost, and conclude that it should be regarded as the "canary in the coal mine" with respect to horizontal transfer in angiosperm mitochondria.

The relevance of phylogeny to studies of global change

Phylogenetic thinking has infiltrated many areas of biological research, but has had little impact on studies of global ecology or climate change. Here, we illustrate how phylogenetic information can be relevant to understanding vegetation-atmosphere dynamics at ecosystem or global scales by re-analyzing a data set of carbonic anhydrase (CA) activity in leaves that was used to estimate terrestrial gross primary productivity. The original calculations relied on what appeared to be low CA activity exclusively in C4 grasses, but our analyses indicate that such activity might instead characterize the PACCAD grass lineage, which includes many widespread C3 species. We outline how phylogenetics can guide better taxon sampling of key physiological traits, and discuss how the emerging field of phyloinformatics presents a promising new framework for scaling from organism physiology to global processes.

Characterization of genes for a pollen allergen, Cry j 2, of Cryptomeria japonica

BACKGROUND

Cry j 2 is one of the major pollen allergens of Cryptomeria japonica. The polymorphism of Cry j 2 isoforms and the conservation of the structure of Cry j 2 in coniferous species remain to be analyzed.

METHODS

A cDNA library derived from the pollen of C. japonica was screened using a fragment of Cry j 2 cDNA. Restriction fragment length polymorphism analysis was performed to examine the diversity of Cry j 2 genes. The promoters of Cry j 2 genes were isolated with a commercially available cloning kit. Clonal variations in the expression of Cry j 2 in pollen were examined by RNA gel blot analysis, and the conservation of the structure of the Cry j 2 gene in coniferous species was evaluated by DNA gel blot analysis.

RESULTS

We isolated three cDNA clones encoding novel isoforms of Cry j 2. We also sequenced a total of 16 promoter regions from 10 specimens. The sequences of promoter regions of Cry j 2 genes were highly divergent. The amount of Cry j 2 mRNA also varied considerably. The Cry j 2 gene was found to be conserved among species belonging to Taxodiaceae and Cupressaceae but to vary between Taxodiaceae and Pinaceae.

CONCLUSIONS

The coding and promoter regions of Cry j 2 genes contain large numbers of polymorphisms. Our analysis revealed large variations in the expression of Cry j 2 at the transcriptional level, and we suggest that conserved homologs of Cry j 2 confer cross-allergenicity among Taxodiaceae and Cupressaceae.

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 class III homeodomain-leucine zipper genes in streptophytes

Land plants underwent tremendous evolutionary change following the divergence of the ancestral lineage from algal relatives. Several important developmental innovations appeared as the embryophyte clade diversified, leading to the appearance of new organs and tissue types. To understand how these changes came about, we need to identify the fundamental genetic developmental programs that are responsible for growth, patterning, and differentiation and describe how these programs were modified and elaborated through time to produce novel morphologies. Class III homeodomain-leucine zipper (class III HD-Zip) genes, identified in the model plant Arabidopsis thaliana, provide good candidates for basic land plant patterning genes. We show that these genes may have evolved in a common ancestor of land plants and their algal sister group and that the gene family has diversified as land plant lineages have diversified. Phylogenetic analysis, expression data from nonflowering lineages, and evidence from Arabidopsis and other flowering plants indicate that class III HD-Zip genes acquired new functions in sporophyte apical growth, vascular patterning and differentiation, and leaf development. Modification of expression patterns that accompanied diversification of class III HD-Zip genes likely played an important role in the evolution of land plant form.

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.

Molecular evolution of the AP2 subfamily

The AP2 (APETALA2)/EREBP (Ethylene Responsive Element Binding Protein) multigene family includes developmentally and physiologically important transcription factors. AP2/EREBP genes are divided into two subfamilies: AP2 genes with two AP2 domains and EREBP genes with a single AP2/ERF (Ethylene Responsive Element Binding Factor) domain. Based on previous phylogenetic analyses, AP2 genes can be divided into two clades, AP2 and ANT groups. To clarify the molecular evolution of the AP2 subfamily, we isolated and sequenced genes with two AP2 domains from three gymnosperms, Cycas revoluta, Ginkgo biloba, and Gnetum parvifolium,as well as from the moss Physcomitrella patens. Expressions of AP2-like genes, including AP2, in Arabidopsis thaliana are regulated by the microRNA miR172. We found that the target site of miR172 is significantly conserved in gymnosperm AP2 homologs, suggesting that regulatory mechanisms of gene expression using microRNA have been conserved over the three hundred million years since the divergence of gymnosperm and flowering plant lineages. We inferred a phylogenetic relationship of these genes with the green alga Chlamydomonas reinhardtii and seed-plant genes available in public DNA databases. The phylogenetic tree showed that the AP2 subfamily diverged into the AP2 and ANT groups before the last common ancestor of land plants and after C. reinhardtii diverged from the land-plant lineage. The tree also indicated that each AP2 and ANT group further diverged into several clades through gene duplications prior to the divergence of gymnosperms and angiosperms.

Comparative analysis of expressed sequence tags of conifers and angiosperms reveals sequences specifically conserved in conifers

To identify and characterize lineage-specific genes of conifers, two sets of ESTs (with 12791 and 5902 ESTs, representing 5373 and 3018 gene transcripts, respectively) were generated from the Cupressaceae species Cryptomeria japonica and Chamaecyparis obtusa. These transcripts were compared with non-redundant sets of genes generated from Pinaceae species, other gymnosperms and angiosperms. About 6% of tentative unique genes (Unigenes) of C. japonica and C. obtusa had homologs in other conifers but not angiosperms, and about 70% had apparent homologs in angiosperms. The calculated GC contents of orthologous genes showed that GC contents of coniferous genes are likely to be lower than those of angiosperms. Comparisons of the numbers of homologous genes in each species suggest that copy numbers of genes may be correlated between diverse seed plants. This correlation suggests that the multiplicity of such genes may have arisen before the divergence of gymnosperms and angiosperms.

Resolution of a deep animal divergence by the pattern of intron conservation

The relationship between three biologically important groups, arthropods, nematodes, and deuterostomes, remains unresolved. It is unknown whether arthropods are more closely related to nematodes (consistent with the "ecdysozoa" hypothesis) or to deuterostomes (consistent with "coelomata"). We present a method in which we use the pattern of spliceosomal intron conservation to develop a series of inequalities that characterize each possible relationship. We find that only the ecdysozoa grouping satisfies these predictions, with P < 10(-6). Simulations show that our method, unlike some previous methods, is largely insensitive to rate variation between branches.

Many independent origins of trans splicing of a plant mitochondrial group II intron

We examined the cis- vs. trans-splicing status of the mitochondrial group II intron nad1i728 in 439 species (427 genera) of land plants, using both Southern hybridization results (for 416 species) and intron sequence data from the literature. A total of 164 species (157 genera), all angiosperms, was found to have a trans-spliced form of the intron. Using a multigene land plant phylogeny, we infer that the intron underwent a transition from cis to trans splicing 15 times among the sampled angiosperms. In 10 cases, the intron was fractured between its 5' end and the intron-encoded matR gene, while in the other 5 cases the fracture occurred between matR and the 3' end of the intron. The 15 intron fractures took place at different time depths during the evolution of angiosperms, with those in Nymphaeales, Austrobaileyales, Chloranthaceae, and eumonocots occurring early in angiosperm evolution and those in Syringodium filiforme, Hydrocharis morsus- ranae, Najas, and Erodium relatively recently. The trans-splicing events uncovered in Austrobaileyales, eumonocots, Polygonales, Caryophyllales, Sapindales, and core Rosales reinforce the naturalness of these major clades of angiosperms, some of which have been identified solely on the basis of recent DNA sequence analyses.

Morphology and affinities of an Early Cretaceous Ephedra (Ephedraceae) from China

Detailed investigations on Lower Cretaceous Ephedra L. fossils (Gnetopsida) reveal morphological characters similar to those of extant Ephedra rhytidosperma Pachomova, including articulate branches with many fine longitudinal striations, a dichasial branching pattern, uni- or bi-ovulate cones with paired bracts, cones terminal on branchlets, and seeds with a short, straight micropylar tubes, covered by numerous regular and prominent transverse laminar protuberances. Fossils are similar to extant E. rhytidosperma reproductive organs but differ in some vegetative structures and are described and discussed here as Ephedra archaeorhytidosperma Y. Yang et al. Because E. rhytidosperma is currently considered one of the most specialized members in Ephedra L. section Pseudobaccatae Stapf, the occurrence of E. archaeorhytidosperma in the Yixian Formation suggests that Ephedra L. was perhaps a more diverse genus in the Lower Cretaceous. Perhaps the evolution and diversity of Ephedra L. was already in place by the Lower Cretaceous and certainly before the end of the Mesozoic.

Distribution of introns in the mitochondrial gene nad1 in land plants: phylogenetic and molecular evolutionary implications

Forty-six species of diverse land plants were investigated by sequencing for their intron content in the mitochondrial gene nad1. A total of seven introns, all belonging to group II, were found, and two were newly discovered in this study. All 13 liverworts examined contain no intron, the same condition as in green algae. Mosses and hornworts, however, share one intron by themselves and another one with vascular plants. These intron distribution patterns are consistent with the hypothesis that liverworts represent the basal-most land plants and that the two introns were gained in the common ancestor of mosses-hornworts-vascular plants after liverworts had diverged. Hornworts also possess a unique intron of their own. A fourth intron was found only in Equisetum L., Marattiaceae, Ophioglossum L., Osmunda L., Asplenium L., and Adiantum L., and was likely acquired in their common ancestor, which supports the monophyly of moniliformopses. Three introns that were previously characterized in angiosperms and a few pteridophytes are now all extended to lycopods, and were likely gained in the common ancestor of vascular plants. Phylogenetic analyses of the intron sequences recovered topologies mirroring those of the plants, suggesting that the introns have all been vertically inherited. All seven nad1 group II introns show broad phylogenetic distribution patterns, with the narrowest being in moniliformopses and hornworts, lineages that date back to at least the Devonian (345 million years ago) and Silurian (435 million years ago), respectively. Hence, these introns must have invaded the genes via ancient transpositional events during the early stage of land plant evolution. Potentially heavy RNA editing was observed in nad1 of Haplomitrium Dedecek, Takakia Hatt. & Inoue, hornworts, Isoetes L., Ophioglossum, and Asplenium. A new nomenclature is proposed for group II introns.

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.

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

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

The 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.

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.

Cross-species amplification of mitochondrial DNA sequence-tagged-site markers in conifers: the nature of polymorphism and variation within and among species in Picea

Primers previously developed to amplify specific non-coding regions of the mitochondrial genome in Angiosperms, and new primers for additional non-coding mtDNA regions, were tested for their ability to direct DNA amplification in 12 conifer taxa and to detect sequence-tagged-site (STS) polymorphisms within and among eight species in Picea. Out of 12 primer pairs, nine were successful at amplifying mtDNA in most of the taxa surveyed. In conifers, indels and substitutions were observed for several loci, allowing them to distinguish between families, genera and, in some cases, between species within genera. In Picea, interspecific polymorphism was detected for four loci, while intraspecific variation was observed for three of the mtDNA regions studied. One of these (SSU rRNA V1 region) exhibited indel polymorphisms, and the two others ( nad1 intron b/c and nad5 intron1) revealed restriction differences after digestion with Sau3AI (PCR-RFLP). A fourth locus, the nad4L- orf25 intergenic region, showed a multibanding pattern for most of the spruce species, suggesting a possible gene duplication. Maternal inheritance, expected for mtDNA in conifers, was observed for all polymorphic markers except the intergenic region nad4L- orf25. Pooling of the variation observed with the remaining three markers resulted in two to six different mtDNA haplotypes within the different species of Picea. Evidence for intra-genomic recombination was observed in at least two taxa. Thus, these mitotypes are likely to be more informative than single-locus haplotypes. They should be particularly useful for the study of biogeography and the dynamics of hybrid zones.

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

Phylogenetic studies based on different types and treatment of data provide substantially conflicting hypotheses of relationships among seed plants. We conducted phylogenetic analyses of sequences of two highly conserved chloroplast genes, psaA and psbB, for a comprehensive taxonomic sample of seed plants and land plants. Parsimony analyses of two different codon position partitions resulted in well-supported, but significantly conflicting, phylogenetic trees. First and second codon positions place angiosperms and gymnosperms as sister clades and Gnetales as sister to Pinaceae. Third positions place Gnetales as sister to all other seed plants. Maximum likelihood trees for the two partitions are also in conflict. Relationships among the main seed plant clades according to first and second positions are similar to those found in parsimony analysis for the same data, but the third position maximum likelihood tree is substantially different from the corresponding parsimony tree, although it agrees partially with the first and second position trees in placing Gnetales as the sister group of Pinaceae. Our results document high rate heterogeneity among lineages, which, together with the greater average rate of substitution for third positions, may reduce phylogenetic signal due to long-branch attraction in parsimony reconstructions. Whereas resolution of relationships among major seed plant clades remains pending, this study provides increased support for relationships within major seed plant clades.