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

Characterizing conflict and congruence of molecular evolution across organellar genome sequences for phylogenetics in land plants

Chloroplasts and mitochondria each contain their own genomes, which have historically been and continue to be important sources of information for inferring the phylogenetic relationships among land plants. The organelles are predominantly inherited from the same parent, and therefore should exhibit phylogenetic concordance. In this study, we examine the mitochondrion and chloroplast genomes of 226 land plants to infer the degree of similarity between the organelles' evolutionary histories. Our results show largely concordant topologies are inferred between the organelles, aside from four well-supported conflicting relationships that warrant further investigation. Despite broad patterns of topological concordance, our findings suggest that the chloroplast and mitochondrial genomes evolved with significant differences in molecular evolution. The differences result in the genes from the chloroplast and the mitochondrion preferentially clustering with other genes from their respective organelles by a program that automates selection of evolutionary model partitions for sequence alignments. Further investigation showed that changes in compositional heterogeneity are not always uniform across divergences in the land plant tree of life. These results indicate that although the chloroplast and mitochondrial genomes have coexisted for over 1 billion years, phylogenetically, they are still evolving sufficiently independently to warrant separate models of evolution. As genome sequencing becomes more accessible, research into these organelles' evolution will continue revealing insight into the ancient cellular events that shaped not only their history, but the history of plants as a whole.

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.

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.

Evolution of mitochondrial RNA editing in extant gymnosperms

To unveil the evolution of mitochondrial RNA editing in gymnosperms, we characterized mitochondrial genomes (mitogenomes), plastid genomes, RNA editing sites, and pentatricopeptide repeat (PPR) proteins from 10 key taxa representing four of the five extant gymnosperm clades. The assembled mitogenomes vary in gene content due to massive gene losses in Gnetum and Conifer II clades. Mitochondrial gene expression levels also vary according to protein function, with the most highly expressed genes involved in the respiratory complex. We identified 9132 mitochondrial C-to-U editing sites, as well as 2846 P-class and 8530 PLS-class PPR proteins. Regains of editing sites were demonstrated in Conifer II rps3 transcripts whose corresponding mitogenomic sequences lack introns due to retroprocessing. Our analyses reveal that non-synonymous editing is efficient and results in more codons encoding hydrophobic amino acids. In contrast, synonymous editing, although performed with variable efficiency, can increase the number of U-ending codons that are preferentially utilized in gymnosperm mitochondria. The inferred loss-to-gain ratio of mitochondrial editing sites in gymnosperms is 2.1:1, of which losses of non-synonymous editing are mainly due to genomic C-to-T substitutions. However, such substitutions only explain a small fraction of synonymous editing site losses, indicating distinct evolutionary mechanisms. We show that gymnosperms have experienced multiple lineage-specific duplications in PLS-class PPR proteins. These duplications likely contribute to accumulated RNA editing sites, as a mechanistic correlation between RNA editing and PLS-class PPR proteins is statistically supported.

Deep evaluation of the evolutionary history of the Heat Shock Factor (HSF) gene family and its expansion pattern in seed plants

Heat shock factor (HSF) genes are essential in some of the basic developmental pathways in plants. Despite extensive studies on the structure, functional diversification, and evolution of HSF genes, their divergence history and gene duplication pattern remain unknown. To further illustrate the probable divergence patterns in these subfamilies, we analyzed the evolutionary history of HSF genes using phylogenetic reconstruction and genomic syntenic analyses, taking advantage of the increased sampling of genomic data from pteridophytes, gymnosperms and basal angiosperms. We identified a novel clade that includes HSFA2, HSFA6, HSFA7, and HSFA9 with a complex relationship, which is very likely due to orthologous or paralogous genes retained after frequent gene duplication events. We hypothesized that HSFA9 derives from HSFA2 through gene duplication in eudicots at the ancestral state, and then expanded in a lineage-specific way. Our findings indicate that HSFB3 and HSFB5 emerged before the divergence of ancestral angiosperms, but were lost in the most recent common ancestors of monocots. We also presumed that HSFC2 derives from HSFC1 in ancestral monocots. This work proposes that during the radiation of flowering plants, an era during which there was a differentiation of angiosperms, the size of the HSF gene family was also being adjusted with considerable sub- or neo-functionalization. The independent evolution of HSFs in eudicots and monocots, including lineage-specific gene duplication, gave rise to a new gene in ancestral eudicots and monocots, and lineage-specific gene loss in ancestral monocots. Our analyses provide essential insights for studying the evolutionary history of this multigene family.

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.

Insights from an ancient gymnosperm lineage: ambient temperature and light and the timing of thermogenesis in cycad cones

PREMISE

Although maintaining the appropriate mid-day timing of the diel thermogenic events of cones of the dioecious cycads Macrozamia lucida and M. macleayi is central to the survival of both plant and pollinator in this obligate pollination mutualism, the nature of the underlying mechanism remains obscure. We investigated whether it is under circadian control. Circadian mechanisms control the timing of many ecologically important processes in angiosperms, yet only a few gymnosperms have been studied in this regard.

METHODS

We subjected cones to different ambient temperature and lighting regimens (constant temperature and darkness; stepwise cool/warm ambient temperatures in constant darkness; stepwise dark/light exposures at constant temperature) to determine whether the resulting timing of their thermogenic events was consistent with circadian control.

RESULTS

Cones exposed to constant ambient temperature and darkness generated multiple temperature peaks endogenously, with an average interpeak-temperature period of 20.7 (±0.20) h that is temperature-compensated (Q₁₀  = 1.02). Exposure to 24-h ambient temperature cycles (12 h cool/12 h warm, constant darkness) yielded an interpeak-temperature period of 24.0 (±0.05) h, accurately and precisely replicating the ambient temperature period. Exposure to 24-h photo-cycles (12 h light/12 h dark, constant ambient temperature) yielded a shorter, more variable interpeak-temperature period of 23 (±0.23) h.

CONCLUSIONS

Our results indicate that cycad cone thermogenesis is under circadian clock control and differentially affected by ambient temperature and light cycles. Our data from cycads (an ancient gymnosperm lineage) adds to what little is known about circadian timing in gymnosperms, which have rarely been studied from the circadian perspective.

The dynamic history of gymnosperm plastomes: Insights from structural characterization, comparative analysis, phylogenomics, and time divergence

Gymnosperms are among the most endangered groups of plant species; they include ginkgo, pines (Conifers I), cupressophytes (Conifers II), cycads, and gnetophytes. The relationships among the five extant gymnosperm groups remain equivocal. We analyzed 167 available gymnosperm plastomes and investigated their diversity and phylogeny. We found that plastome size, structure, and gene order were highly variable in the five gymnosperm groups, of which Parasitaxus usta (Vieill.) de Laub. and Macrozamia mountperriensis F.M.Bailey had the smallest and largest plastomes, respectively. The inverted repeats (IRs) of the five groups were shown to have evolved through distinctive evolutionary scenarios. The IRs have been lost in all conifers but retained in cycads and gnetophytes. A positive association between simple sequence repeat (SSR) abundance and plastome size was observed, and the SSRs with the most variation were found in Pinaceae. Furthermore, the number of repeats was negatively correlated with IR length; thus, the highest number of repeats was detected in Conifers I and II, in which the IRs had been lost. We constructed a phylogeny based on 29 shared genes from 167 plastomes. With the plastome tree and 13 calibrations, we estimated the tree height between present-day angiosperms and gymnosperms to be ∼380 million years ago (mya). The placement of Gnetales in the tree agreed with the Gnetales-other gymnosperms hypothesis. The divergence between Ginkgo and cycads was estimated as ∼284 mya; the crown age of the cycads was 251 mya. Our time-calibrated plastid-based phylogenomic tree provides a framework for comparative studies of gymnosperm evolution.

Endosperm of Angiosperms and Genomic Imprinting

Modern ideas about the role of epigenetic systems in the regulation of gene expression allow us to understand the mechanisms of vital activities in plants, such as genomic imprinting. It is important that genomic imprinting is known first and foremost for the endosperm, which not only provides an embryo with necessary nutrients, but also plays a special biological role in the formation of seeds and fruits. Available data on genomic imprinting in the endosperm have been obtained only for the triploid endosperm in model plants, which develops after double fertilization in a Polygonum-type embryo sac, the most common type among angiosperms. Here we provide a brief overview of a wide diversity of embryo sacs and endosperm types and ploidy levels, as well as their distribution in the angiosperm families, positioned according to the Angiosperm Phylogeny Group IV (APG IV) phylogenetic classification. Addition of the new, non-model taxa to study gene imprinting in seed development will extend our knowledge about the epigenetic mechanisms underlying angiosperm fertility.

An ancient push-pull pollination mechanism in cycads

Most cycads engage in brood-site pollination mutualisms, yet the mechanism by which the Cycadales entice pollination services from diverse insect mutualists remains unknown. Here, we characterize a push-pull pollination mechanism between a New World cycad and its weevil pollinators that mirrors the mechanism between a distantly related Old World cycad and its thrips pollinators. The behavioral convergence between weevils and thrips, combined with molecular phylogenetic dating and a meta-analysis of thermogenesis and coordinated patterns of volatile attraction and repulsion suggest that a push-pull pollination mutualism strategy is ancestral in this ancient, dioecious plant group. Hence, it may represent one of the earliest insect/plant pollination mechanisms, arising long before the evolution of visual floral signaling commonly used by flowering plants.

Reconstructing Krassilovia mongolica supports recognition of a new and unusual group of Mesozoic conifers

Previously unrecognized anatomical features of the cone scales of the enigmatic Early Cretaceous conifer Krassilovia mongolica include the presence of transversely oriented paracytic stomata, which is unusual for all other extinct and extant conifers. Identical stomata are present on co-occurring broad, linear, multiveined leaves assigned to Podozamites harrisii, providing evidence that K. mongolica and P. harrisii are the seed cones and leaves of the same extinct plant. Phylogenetic analyses of the relationships of the reconstructed Krassilovia plant place it in an informal clade that we name the Krassilovia Clade, which also includes Swedenborgia cryptomerioides-Podozamites schenkii, and Cycadocarpidium erdmanni-Podozamites schenkii. All three of these plants have linear leaves that are relatively broad compared to most living conifers, and that are also multiveined with transversely oriented paracytic stomata. We propose that these may be general features of the Krassilovia Clade. Paracytic stomata, and other features of this new group, recall features of extant and fossil Gnetales, raising questions about the phylogenetic homogeneity of the conifer clade similar to those raised by phylogenetic analyses of molecular data.

Phylogenetic Conflicts, Combinability, and Deep Phylogenomics in Plants

Studies have demonstrated that pervasive gene tree conflict underlies several important phylogenetic relationships where different species tree methods produce conflicting results. Here, we present a means of dissecting the phylogenetic signal for alternative resolutions within a data set in order to resolve recalcitrant relationships and, importantly, identify what the data set is unable to resolve. These procedures extend upon methods for isolating conflict and concordance involving specific candidate relationships and can be used to identify systematic error and disambiguate sources of conflict among species tree inference methods. We demonstrate these on a large phylogenomic plant data set. Our results support the placement of Amborella as sister to the remaining extant angiosperms, Gnetales as sister to pines, and the monophyly of extant gymnosperms. Several other contentious relationships, including the resolution of relationships within the bryophytes and the eudicots, remain uncertain given the low number of supporting gene trees. To address whether concatenation of filtered genes amplified phylogenetic signal for relationships, we implemented a combinatorial heuristic to test combinability of genes. We found that nested conflicts limited the ability of data filtering methods to fully ameliorate conflicting signal amongst gene trees. These analyses confirmed that the underlying conflicting signal does not support broad concatenation of genes. Our approach provides a means of dissecting a specific data set to address deep phylogenetic relationships while also identifying the inferential boundaries of the data set. [Angiosperms; coalescent; gene-tree conflict; genomics; phylogenetics; phylogenomics.]

Two Independent Plastid accD Transfers to the Nuclear Genome of Gnetum and Other Insights on Acetyl-CoA Carboxylase Evolution in Gymnosperms

Acetyl-CoA carboxylase (ACCase) is the key regulator of fatty acid biosynthesis. In most plants, ACCase exists in two locations (cytosol and plastids) and in two forms (homomeric and heteromeric). Heteromeric ACCase comprises four subunits, three of them (ACCA-C) are nuclear encoded (nr) and the fourth (ACCD) is usually plastid encoded. Homomeric ACCase is encoded by a single nr-gene (ACC). We investigated the ACCase gene evolution in gymnosperms by examining the transcriptomes of newly sequenced Gnetum ula, combined with 75 transcriptomes and 110 plastomes of other gymnosperms. AccD-coding sequences are elongated through the insertion of repetitive DNA in four out of five cupressophyte families (except Sciadopityaceae) and were functionally transferred to the nucleus of gnetophytes and Sciadopitys. We discovered that, among the three genera of gnetophytes, only Gnetum has two copies of nr-accD. Furthermore, using protoplast transient expression assays, we experimentally verified that the nr-accD precursor proteins in Gnetum and Sciadopitys can be delivered to the plastids. Of the two nr-accD copies of Gnetum, one dually targets plastids and mitochondria, whereas the other potentially targets plastoglobuli. The distinct transit peptides, gene architectures, and flanking sequences between the two Gnetum accDs suggest that they have independent origins. Our findings are the first account of two distinctly targeted nr-accDs of any green plants and the most comprehensive analyses of ACCase evolution in gymnosperms to date.

One thousand plant transcriptomes and the phylogenomics of green plants

Green plants (Viridiplantae) include around 450,000-500,000 species1,2 of great diversity and have important roles in terrestrial and aquatic ecosystems. Here, as part of the One Thousand Plant Transcriptomes Initiative, we sequenced the vegetative transcriptomes of 1,124 species that span the diversity of plants in a broad sense (Archaeplastida), including green plants (Viridiplantae), glaucophytes (Glaucophyta) and red algae (Rhodophyta). Our analysis provides a robust phylogenomic framework for examining the evolution of green plants. Most inferred species relationships are well supported across multiple species tree and supermatrix analyses, but discordance among plastid and nuclear gene trees at a few important nodes highlights the complexity of plant genome evolution, including polyploidy, periods of rapid speciation, and extinction. Incomplete sorting of ancestral variation, polyploidization and massive expansions of gene families punctuate the evolutionary history of green plants. Notably, we find that large expansions of gene families preceded the origins of green plants, land plants and vascular plants, whereas whole-genome duplications are inferred to have occurred repeatedly throughout the evolution of flowering plants and ferns. The increasing availability of high-quality plant genome sequences and advances in functional genomics are enabling research on genome evolution across the green tree of life.

Plant biodiversity and regulation of photosynthesis in the natural environment

Investigation of photosynthesis regulation in different plant groups exposed to variable conditions showed that all species have similar photosynthetic electron transport modulation while excess energy dissipation is species specific. Photosynthesis is regulated in response to dynamic environmental conditions to satisfy plant metabolic demands while also avoiding possible over-excitation of the electron transport chain and the generation of harmful reactive oxygen species. Photosynthetic organisms evolved several mechanisms to modulate light harvesting and electron transport efficiency to respond to conditions changing at different timescales, going from fast sun flecks to slow seasonal variations. These regulatory mechanisms changed during evolution of photosynthetic organisms, also adapting to various ecological niches, making the investigation of plant biodiversity highly valuable to uncover conserved traits and plasticity of photosynthetic regulation and complement studies on model species. In this work, a set of plants belonging to different genera of angiosperms, gymnosperms, ferns and lycophytes were investigated by monitoring their photosynthetic parameters in different seasons looking for common trends and differences. In all plants, analysed photosynthetic electron transport rate was found to be modulated by growth light intensity, ensuring a balance between available energy and photochemical capacity. Growth light also influenced the threshold where heat dissipation of excitation energy, a mechanism called non-photochemical quenching (NPQ), was activated. On the contrary, NPQ amplitude did not correlate with light intensity experienced by the plants but was a species-specific feature. The zeaxanthin-dependent component of NPQ, qZ, was found to be the most variable in different plants and its modulation influenced the intensity and the kinetic properties of the response.

Origin of the Taxaceae aril: evolutionary implications of seed-cone teratologies in Pseudotaxus chienii

BACKGROUND AND AIMS

Fleshy structures that promote biotic dispersal by ingestion have evolved many times in seed plants. Within the yew family Taxaceae sensu lato (six genera, including Cephalotaxus), it remains controversial whether the characteristic fleshy structure surrounding the seed is interpreted as a novel outgrowth of the base of the ovule (i.e. an aril) or a fleshy seed coat that is entirely derived from the integument (i.e. a sarcotesta). This paper presents a detailed study of both wild-type and teratological seed cones of Pseudotaxus chienii, including morphology, anatomy and ontogeny.

METHODS

Wild-type and teratological seed cones were investigated with the classical paraffin technique and subsequent astrablue/safranin staining and scanning electron microscopy.

KEY RESULTS

The wild-type seed cone of Pseudotaxus possesses a fleshy white aril that is cup-like and not entirely fused to the seed. In the teratological seed cones investigated, the aril was bilobed and consisted of two free halves. In both wild-type and teratological cones, the aril was initiated as two lateral primordia in a transverse plane, but in wild-type cones the two primordia became extended into a ring primordium, which grew apically, leading to the cup-like shape. The teratological seed cones lacked a ring primordium and the two lateral aril lobes remained free throughout their entire ontogeny, alternating with the scale-like leaves inserted below them on the same branch; in some cases, these leaves also became fleshy.

CONCLUSIONS

Based on the ontogeny and arrangement of the two fleshy aril lobes in the teratological seed cones of Pseudotaxus, we suggest that the typical aril of Taxaceae could be readily interpreted as a fused pair of strongly swollen leaves rather than a modified integument. Our investigations of the cup-like aril of Pseudotaxus demonstrate a similarity not only with other Taxaceae but also with relatively distantly related conifers such as Phyllocladus (Podocarpaceae).

Insights on the persistence of pines (Pinus species) in the Late Cretaceous and their increasing dominance in the Anthropocene

Although gymnosperms were nearly swept away by the rise of the angiosperms in the Late Cretaceous, conifers, and pines (Pinus species) in particular, survived and regained their dominance in some habitats. Diversification of pines into fire-avoiding (subgenus Haploxylon) and fire-adapted (subgenus Diploxylon) species occurred in response to abiotic and biotic factors in the Late Cretaceous such as competition with emerging angiosperms and changing fire regimes. Adaptations/traits that evolved in response to angiosperm-fuelled fire regimes and stressful environments in the Late Cretaceous were key to pine success and are also contributing to a new "pine rise" in some areas in the Anthropocene. Human-mediated activities exert both positive and negative impacts of range size and expansion and invasions of pines. Large-scale afforestation with pines, human-mediated changes to fire regimes, and other ecosystem processes are other contributing factors. We discuss traits that evolved in response to angiosperm-mediated fires and stressful environments in the Cretaceous and that continue to contribute to pine persistence and dominance and the numerous ways in which human activities favor pines.

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.

Incorporating fossils into the Pinaceae tree of life

PREMISE OF THE STUDY

Pinaceae have a rich but enigmatic early fossil record, much of which is represented by permineralized seed cones. Our incomplete knowledge of morphology and anatomy in living and extinct species poses an important barrier to understanding their phylogenetic relationships and timing of diversification.

METHODS

We expanded a morphology matrix to 46 fossil and 31 extant Pinaceae species, mainly adding characters from stem and leaf anatomy and seed cones. Using parsimony and Bayesian inference, we compared phylogenetic relationships for extant taxa with and without fossils from the morphology matrix combined with an alignment of plastid gene sequences.

KEY RESULTS

Combined analysis of morphological and molecular characters resulted in a phylogeny of extant Pinaceae that was robust at all nodes except those relating to the interrelationships of Pinus, Picea, and Cathaya and the position of Cedrus. Simultaneous analysis of all fossil and extant species did not result in changes in the relationships among the extant species but did greatly reduce branch support. We found that the placement of most fossils was sensitive to the method of phylogenetic reconstruction when analyzing them singly with the extant species.

CONCLUSIONS

A robust phylogenetic hypothesis for the main lineages of Pinaceae is emerging. Most Early Cretaceous fossils are stem or crown lineages of Pinus, but close relationships also were found between fossils and several other extant genera. The phylogenetic position of fossils broadly supports the existence of extant genera in the Lower Cretaceous.

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.

Single-Copy Genes as Molecular Markers for Phylogenomic Studies in Seed Plants

Phylogenetic relationships among seed plant taxa, especially within the gymnosperms, remain contested. In contrast to angiosperms, for which several genomic, transcriptomic and phylogenetic resources are available, there are few, if any, molecular markers that allow broad comparisons among gymnosperm species. With few gymnosperm genomes available, recently obtained transcriptomes in gymnosperms are a great addition to identifying single-copy gene families as molecular markers for phylogenomic analysis in seed plants. Taking advantage of an increasing number of available genomes and transcriptomes, we identified single-copy genes in a broad collection of seed plants and used these to infer phylogenetic relationships between major seed plant taxa. This study aims at extending the current phylogenetic toolkit for seed plants, assessing its ability for resolving seed plant phylogeny, and discussing potential factors affecting phylogenetic reconstruction. In total, we identified 3,072 single-copy genes in 31 gymnosperms and 2,156 single-copy genes in 34 angiosperms. All studied seed plants shared 1,469 single-copy genes, which are generally involved in functions like DNA metabolism, cell cycle, and photosynthesis. A selected set of 106 single-copy genes provided good resolution for the seed plant phylogeny except for gnetophytes. Although some of our analyses support a sister relationship between gnetophytes and other gymnosperms, phylogenetic trees from concatenated alignments without 3rd codon positions and amino acid alignments under the CAT + GTR model, support gnetophytes as a sister group to Pinaceae. Our phylogenomic analyses demonstrate that, in general, single-copy genes can uncover both recent and deep divergences of seed plant phylogeny.

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.

Multiple measures could alleviate long-branch attraction in phylogenomic reconstruction of Cupressoideae (Cupressaceae)

Long-branch attraction (LBA) is a major obstacle in phylogenetic reconstruction. The phylogenetic relationships among Juniperus (J), Cupressus (C) and the Hesperocyparis-Callitropsis-Xanthocyparis (HCX) subclades of Cupressoideae are controversial. Our initial analyses of plastid protein-coding gene matrix revealed both J and C with much longer stem branches than those of HCX, so their sister relationships may be attributed to LBA. We used multiple measures including data filtering and modifying, evolutionary model selection and coalescent phylogenetic reconstruction to alleviate the LBA artifact. Data filtering by strictly removing unreliable aligned regions and removing substitution saturation genes and rapidly evolving sites could significantly reduce branch lengths of subclades J and C and recovered a relationship of J (C, HCX). In addition, using coalescent phylogenetic reconstruction could elucidate the LBA artifact and recovered J (C, HCX). However, some valid methods for other taxa were inefficient in alleviating the LBA artifact in J-C-HCX. Different strategies should be carefully considered and justified to reduce LBA in phylogenetic reconstruction of different groups. Three subclades of J-C-HCX were estimated to have experienced ancient rapid divergence within a short period, which could be another major obstacle in resolving relationships. Furthermore, our plastid phylogenomic analyses fully resolved the intergeneric relationships of Cupressoideae.

Relative strength of fine-scale spatial genetic structure in paternally vs biparentally inherited DNA in a dioecious plant depends on both sex proportions and pollen-to-seed dispersal ratio

In plants, the spatial genetic structure (SGS) is shaped mainly by gene dispersal and effective population density. Among additional factors, the mode of DNA inheritance and dioecy influence SGS. However, their joint impact on SGS remains unclear, especially in the case of paternally inherited DNA. Using theoretical approximations and computer simulations, here we showed that the relative intensity of SGS measured in paternally and biparentally inherited DNA in a dioecious plant population depends on both the proportion of males and the pollen-to-seed dispersal ratio. As long as males do not prevail in a population, SGS is more intense in paternally than biparentally inherited DNA. When males prevail, the intensity of SGS in paternally vs biparentally inherited DNA depends on the compound effect of sex proportions and the pollen-to-seed dispersal ratio. To empirically validate our predictions, we used the case of Taxus baccata, a dioecious European tree. First, we showed that mitochondrial DNA (mtDNA) in T. baccata is predominantly (98%) paternally inherited. Subsequently, using nuclear DNA (nuDNA) and mitochondrial microsatellite data, we compared the fine-scale SGS intensity at both marker types in two natural populations. The population with equal sex proportions showed stronger SGS in mtDNA than in nuDNA. On the other hand, we found lower SGS intensity in mtDNA than in nuDNA in the population with 67% males. Thus, the empirical results provided good support for the theoretical predictions, suggesting that knowledge about SGS in paternally vs biparentally inherited DNA may provide insight into effective sex proportions within dioecious populations.

Birth of Four Chimeric Plastid Gene Clusters in Japanese Umbrella Pine

Many genes in the plastid genomes (plastomes) of plants are organized as gene clusters, in which genes are co-transcribed, resembling bacterial operons. These plastid operons are highly conserved, even among conifers, whose plastomes are highly rearranged relative to other seed plants. We have determined the complete plastome sequence of Sciadopitys verticillata (Japanese umbrella pine), the sole member of Sciadopityaceae. The Sciadopitys plastome is characterized by extensive inversions, pseudogenization of four tRNA genes after tandem duplications, and a unique pair of 370-bp inverted repeats involved in the formation of isomeric plastomes. We showed that plastomic inversions in Sciadopitys have led to shuffling of the remote conserved operons, resulting in the birth of four chimeric gene clusters. Our data also demonstrated that the relocated genes can be co-transcribed in these chimeric gene clusters. The plastome of Sciadopitys advances our current understanding of how the conifer plastomes have evolved toward increased diversity and complexity.

Multiplexed shotgun sequencing reveals congruent three-genome phylogenetic signals for four botanical sections of the flax genus Linum

A genome-wide detection of phylogenetic signals by next generation sequencing (NGS) has recently emerged as a promising genomic approach for phylogenetic analysis of non-model organisms. Here we explored the use of a multiplexed shotgun sequencing method to assess the phylogenetic relationships of 18 Linum samples representing 16 species within four botanical sections of the flax genus Linum. The whole genome DNAs of 18 Linum samples were fragmented, tagged, and sequenced using an Illumina MiSeq. Acquired sequencing reads per sample were further separated into chloroplast, mitochondrial and nuclear sequence reads. SNP calls upon genome-specific sequence data sets revealed 6143 chloroplast, 2673 mitochondrial, and 19,562 nuclear SNPs. Phylogenetic analyses based on three-genome SNP data sets with and without missing observations showed congruent three-genome phylogenetic signals for four botanical sections of the Linum genus. Specifically, two major lineages showing a separation of Linum-Dasylinum sections and Linastrum-Syllinum sections were confirmed. The Linum section displayed three major branches representing two major evolutionary stages leading to cultivated flax. Cultivated flax and its immediate progenitor were formed as its own branch, genetically more closely related to L. decumbens and L. grandiflorum with chromosome count of eight, and distantly apart from six other species with chromosome count of nine. Five species of the Linastrum and Syllinum sections were genetically more distant from cultivated flax, but they appeared to be more closely related to each other, even with variable chromosome counts. These findings not only provide the first evidence of congruent three-genome phylogenetic pathways within the Linum genus, but also demonstrate the utility of the multiplexed shotgun sequencing in acquisition of three-genome phylogenetic signals of non-model organisms.

An Ancient Transkingdom Horizontal Transfer of Penelope-Like Retroelements from Arthropods to Conifers

Comparative genomics analyses empowered by the wealth of sequenced genomes have revealed numerous instances of horizontal DNA transfers between distantly related species. In eukaryotes, repetitive DNA sequences known as transposable elements (TEs) are especially prone to move across species boundaries. Such horizontal transposon transfers, or HTTs, are relatively common within major eukaryotic kingdoms, including animals, plants, and fungi, while rarely occurring across these kingdoms. Here, we describe the first case of HTT from animals to plants, involving TEs known as Penelope-like elements, or PLEs, a group of retrotransposons closely related to eukaryotic telomerases. Using a combination of in situ hybridization on chromosomes, polymerase chain reaction experiments, and computational analyses we show that the predominant PLE lineage, EN(+)PLEs, is highly diversified in loblolly pine and other conifers, but appears to be absent in other gymnosperms. Phylogenetic analyses of both protein and DNA sequences reveal that conifers EN(+)PLEs, or Dryads, form a monophyletic group clustering within a clade of primarily arthropod elements. Additionally, no EN(+)PLEs were detected in 1,928 genome assemblies from 1,029 nonmetazoan and nonconifer genomes from 14 major eukaryotic lineages. These findings indicate that Dryads emerged following an ancient horizontal transfer of EN(+)PLEs from arthropods to a common ancestor of conifers approximately 340 Ma. This represents one of the oldest known interspecific transmissions of TEs, and the most conspicuous case of DNA transfer between animals and plants.

Evolutionary loss of light-harvesting proteins Lhcb6 and Lhcb3 in major land plant groups--break-up of current dogma

Photosynthesis in plants and algae relies on the coordinated function of photosystems (PS) I and II. Their efficiency is augmented by finely-tuned light-harvesting proteins (Lhcs) connected to them. The most recent Lhcs (in evolutionary terms), Lhcb6 and Lhcb3, evolved during the transition of plants from water to land and have so far been considered to be an essential characteristic of land plants. We used single particle electron microscopy and sequence analysis to study architecture and composition of PSII supercomplex from Norway spruce and related species. We have found that there are major land plant families that lack functional lhcb6 and lhcb3 genes, which notably changes the organization of PSII supercomplexes. The Lhcb6 and Lhcb3 proteins have been lost in the gymnosperm genera Picea and Pinus (family Pinaceae) and Gnetum (Gnetales). We also revealed that the absence of these proteins in Norway spruce modifies the PSII supercomplex in such a way that it resembles its counterpart in the alga Chlamydomonas reinhardtii, an evolutionarily older organism. Our results break a deep-rooted concept of Lhcb6 and Lhcb3 proteins being the essential characteristic of land plants, and beg the question of what the evolutionary benefit of their loss could be.

Specific Microbial Communities Associate with the Rhizosphere of Welwitschia mirabilis, a Living Fossil

Welwitschia mirabilis is an ancient and rare plant distributed along the western coast of Namibia and Angola. Several aspects of Welwitschia biology and ecology have been investigated, but very little is known about the microbial communities associated with this plant. This study reports on the bacterial and fungal communities inhabiting the rhizosphere of W. mirabilis and the surrounding bulk soil. Rhizosphere communities were dominated by sequences of Alphaproteobacteria and Euromycetes, while Actinobacteria, Alphaproteobacteria, and fungi of the class Dothideomycetes jointly dominated bulk soil communities. Although microbial communities within the rhizosphere and soil samples were highly variable, very few “species” (OTUs defined at a 97% identity cut-off) were shared between these two environments. There was a small ‘core’ rhizosphere bacterial community (formed by Nitratireductor, Steroidobacter, Pseudonocardia and three Phylobacteriaceae) that together with Rhizophagus, an arbuscular mycorrhizal fungus, and other putative plant growth-promoting microbes may interact synergistically to promote Welwitschia growth.

Functional characterization of blue-light-induced responses and PHOTOTROPIN 1 gene in Welwitschia mirabilis

The blue light (BL) receptor phototropin (phot) is specifically found in green plants; it regulates various BL-induced responses such as phototropism, chloroplast movement, stomatal opening, and leaf flattening. In Arabidopsis thaliana, two phototropins--phot1 and phot2--respond to blue light in overlapping but distinct ways. These BL-receptor-mediated responses enhance the photosynthetic activity of plants under weak light and minimize photodamage under strong light conditions. Welwitschia mirabilis Hook.f. found in the Namib Desert, and it has adapted to severe environmental stresses such as limiting water and strong sunlight. Although the plant has physiologically and ecologically unique features, it is unknown whether phototropin is functional in this plant. In this study, we assessed the functioning of phot-mediated BL responses in W. mirabilis. BL-dependent phototropism and stomatal opening was observed but light-dependent chloroplast movement was not detected. We performed a functional analysis of the PHOT1 gene of W. mirabilis, WmPHOT1, in Arabidopsis thaliana. We generated transgenic A. thaliana lines expressing WmPHOT1 in a phot1 phot2 double mutant background. Several Wmphot1 transgenic plants showed normal growth, although phot1 phot2 double mutant plants showed stunted growth. Furthermore, Wmphot1 transgenic plants showed normal phot1-mediated responses including phototropism, chloroplast accumulation, stomatal opening, and leaf flattening, but lacked the chloroplast avoidance response that is specifically mediated by phot2. Thus, our findings indicate that W. mirabilis possesses typical phot-mediated BL responses that were at least partially mediated by functional phototropin 1, an ortholog of Atphot1.

The Diversity of the Pollen Tube Pathway in Plants: Toward an Increasing Control by the Sporophyte

Plants, unlike animals, alternate multicellular diploid, and haploid generations in their life cycle. While this is widespread all along the plant kingdom, the size and autonomy of the diploid sporophyte and the haploid gametophyte generations vary along evolution. Vascular plants show an evolutionary trend toward a reduction of the gametophyte, reflected both in size and lifespan, together with an increasing dependence from the sporophyte. This has resulted in an overlooking of the importance of the gametophytic phase in the evolution of higher plants. This reliance on the sporophyte is most notorious along the pollen tube journey, where the male gametophytes have to travel a long way inside the sporophyte to reach the female gametophyte. Along evolution, there is a change in the scenery of the pollen tube pathway that favors pollen competition and selection. This trend, toward apparently making complicated what could be simple, appears to be related to an increasing control of the sporophyte over the gametophyte with implications for understanding plant evolution.

Non-functional plastid ndh gene fragments are present in the nuclear genome of Norway spruce (Picea abies L. Karsch): insights from in silico analysis of nuclear and organellar genomes

Many genes have been lost from the prokaryote plastidial genome during the early events of endosymbiosis in eukaryotes. Some of them were definitively lost, but others were relocated and functionally integrated to the host nuclear genomes through serial events of gene transfer during plant evolution. In gymnosperms, plastid genome sequencing has revealed the loss of ndh genes from several species of Gnetales and Pinaceae, including Norway spruce (Picea abies). This study aims to trace the ndh genes in the nuclear and organellar Norway spruce genomes. The plastid genomes of higher plants contain 11 ndh genes which are homologues of mitochondrial genes encoding subunits of the proton-pumping NADH-dehydrogenase (nicotinamide adenine dinucleotide dehydrogenase) or complex I (electron transport chain). Ndh genes encode 11 NDH polypeptides forming the Ndh complex (analogous to complex I) which seems to be primarily involved in chloro-respiration processes. We considered ndh genes from the plastidial genome of four gymnosperms (Cryptomeria japonica, Cycas revoluta, Ginkgo biloba, Podocarpus totara) and a single angiosperm species (Arabidopsis thaliana) to trace putative homologs in the nuclear and organellar Norway spruce genomes using tBLASTn to assess the evolutionary fate of ndh genes in Norway spruce and to address their genomic location(s), structure, integrity and functionality. The results obtained from tBLASTn were subsequently analyzed by performing homology search for finding ndh specific conserved domains using conserved domain search. We report the presence of non-functional plastid ndh gene fragments, excepting ndhE and ndhG genes, in the nuclear genome of Norway spruce. Regulatory transcriptional elements like promoters, TATA boxes and enhancers were detected in the upstream regions of some ndh fragments. We also found transposable elements in the flanking regions of few ndh fragments suggesting nuclear rearrangements in those regions. These evidences support the hypothesis that, at least in Picea, ndh translocations from the plastid to the nuclear genome have occurred, and that there might have been a functional machinery at some time during evolution to accommodate them within a nuclear-encoded environment, or attempts to form it.

Astonishing 35S rDNA diversity in the gymnosperm species Cycas revoluta Thunb

In all eukaryotes, the highly repeated 35S ribosomal DNA (rDNA) sequences encoding 18S-5.8S-26S ribosomal RNA (rRNA) typically show high levels of intragenomic uniformity due to homogenisation processes, leading to concerted evolution of 35S rDNA repeats. Here, we compared 35S rDNA divergence in several seed plants using next generation sequencing and a range of molecular and cytogenetic approaches. Most species showed similar 35S rDNA homogeneity indicating concerted evolution. However, Cycas revoluta exhibits an extraordinary diversity of rDNA repeats (nucleotide sequence divergence of different copies averaging 12 %), influencing both the coding and non-coding rDNA regions nearly equally. In contrast, its rRNA transcriptome was highly homogeneous suggesting that only a minority of genes (<20 %) encode functional rRNA. The most common SNPs were C > T substitutions located in symmetrical CG and CHG contexts which were also highly methylated. Both functional genes and pseudogenes appear to cluster on chromosomes. The extraordinary high levels of 35S rDNA diversity in C. revoluta, and probably other species of cycads, indicate that the frequency of repeat homogenisation has been much lower in this lineage, compared with all other land plant lineages studied. This has led to the accumulation of methylation-driven mutations and pseudogenisation. Potentially, the reduced homology between paralogs prevented their elimination by homologous recombination, resulting in long-term retention of rDNA pseudogenes in the genome.

A protocol for targeted enrichment of intron-containing sequence markers for recent radiations: A phylogenomic example from Heuchera (Saxifragaceae)

PREMISE OF THE STUDY

Phylogenetic inference is moving to large multilocus data sets, yet there remains uncertainty in the choice of marker and sequencing method at low taxonomic levels. To address this gap, we present a method for enriching long loci spanning intron-exon boundaries in the genus Heuchera.

METHODS

Two hundred seventy-eight loci were designed using a splice-site prediction method combining transcriptomic and genomic data. Biotinylated probes were designed for enrichment of these loci. Reference-based assembly was performed using genomic references; additionally, chloroplast and mitochondrial genomes were used as references for off-target reads. The data were aligned and subjected to coalescent and concatenated phylogenetic analyses to demonstrate support for major relationships.

RESULTS

Complete or nearly complete (>99%) sequences were assembled from essentially all loci from all taxa. Aligned introns showed a fourfold increase in divergence as opposed to exons. Concatenated analysis gave decisive support to all nodes, and support was also high and relationships mostly similar in the coalescent analysis. Organellar phylogenies were also well-supported and conflicted with the nuclear signal.

DISCUSSION

Our approach shows promise for resolving a recent radiation. Enrichment for introns is highly successful with little or no sequencing dropout at low taxonomic levels despite higher substitution and indel frequencies, and should be exploited in studies of species complexes.

A protocol for targeted enrichment of intron-containing sequence markers for recent radiations: A phylogenomic example from Heuchera (Saxifragaceae)

PREMISE OF THE STUDY

Phylogenetic inference is moving to large multilocus data sets, yet there remains uncertainty in the choice of marker and sequencing method at low taxonomic levels. To address this gap, we present a method for enriching long loci spanning intron-exon boundaries in the genus Heuchera.

METHODS

Two hundred seventy-eight loci were designed using a splice-site prediction method combining transcriptomic and genomic data. Biotinylated probes were designed for enrichment of these loci. Reference-based assembly was performed using genomic references; additionally, chloroplast and mitochondrial genomes were used as references for off-target reads. The data were aligned and subjected to coalescent and concatenated phylogenetic analyses to demonstrate support for major relationships.

RESULTS

Complete or nearly complete (>99%) sequences were assembled from essentially all loci from all taxa. Aligned introns showed a fourfold increase in divergence as opposed to exons. Concatenated analysis gave decisive support to all nodes, and support was also high and relationships mostly similar in the coalescent analysis. Organellar phylogenies were also well-supported and conflicted with the nuclear signal.

DISCUSSION

Our approach shows promise for resolving a recent radiation. Enrichment for introns is highly successful with little or no sequencing dropout at low taxonomic levels despite higher substitution and indel frequencies, and should be exploited in studies of species complexes.

Catharanthus mosaic virus: A potyvirus from a gymnosperm, Welwitschia mirabilis

A virus from a symptomatic plant of the gymnosperm Welwitschia mirabilis Hook. growing as an ornamental plant in a domestic garden in Western Australia was inoculated to a plant of Nicotiana benthamiana where it established a systemic infection. The complete genome sequence of 9636 nucleotides was determined using high-throughput and Sanger sequencing technologies. The genome sequence shared greatest identity (83% nucleotides and 91% amino acids) with available partial sequences of catharanthus mosaic virus, indicating that the new isolate belonged to that taxon. Analysis of the phylogeny of the complete virus sequence placed it in a monotypic group in the genus Potyvirus. This is the first record of a virus from W. mirabilis, the first complete genome sequence of catharanthus mosaic virus determined, and the first record from Australia. This finding illustrates the risk to natural and managed systems posed by the international trade in live plants and propagules, which enables viruses to establish in new regions and infect new hosts.

NDH expression marks major transitions in plant evolution and reveals coordinate intracellular gene loss

BACKGROUND

Key innovations have facilitated novel niche utilization, such as the movement of the algal predecessors of land plants into terrestrial habitats where drastic fluctuations in light intensity, ultraviolet radiation and water limitation required a number of adaptations. The NDH (NADH dehydrogenase-like) complex of Viridiplantae plastids participates in adapting the photosynthetic response to environmental stress, suggesting its involvement in the transition to terrestrial habitats. Although relatively rare, the loss or pseudogenization of plastid NDH genes is widely distributed across diverse lineages of photoautotrophic seed plants and mutants/transgenics lacking NDH function demonstrate little difference from wild type under non-stressed conditions. This study analyzes large transcriptomic and genomic datasets to evaluate the persistence and loss of NDH expression across plants.

RESULTS

Nuclear expression profiles showed accretion of the NDH gene complement at key transitions in land plant evolution, such as the transition to land and at the base of the angiosperm lineage. While detection of transcripts for a selection of non-NDH, photosynthesis related proteins was independent of the state of NDH, coordinate, lineage-specific loss of plastid NDH genes and expression of nuclear-encoded NDH subunits was documented in Pinaceae, gnetophytes, Orchidaceae and Geraniales confirming the independent and complete loss of NDH in these diverse seed plant taxa.

CONCLUSION

The broad phylogenetic distribution of NDH loss and the subtle phenotypes of mutants suggest that the NDH complex is of limited biological significance in contemporary plants. While NDH activity appears dispensable under favorable conditions, there were likely sufficiently frequent episodes of abiotic stress affecting terrestrial habitats to allow the retention of NDH activity. These findings reveal genetic factors influencing plant/environment interactions in a changing climate through 450 million years of land plant evolution.

Transcriptome profiling of khat (Catha edulis) and Ephedra sinica reveals gene candidates potentially involved in amphetamine-type alkaloid biosynthesis

Amphetamine analogues are produced by plants in the genus Ephedra and by khat (Catha edulis), and include the widely used decongestants and appetite suppressants (1S,2S)-pseudoephedrine and (1R,2S)-ephedrine. The production of these metabolites, which derive from L-phenylalanine, involves a multi-step pathway partially mapped out at the biochemical level using knowledge of benzoic acid metabolism established in other plants, and direct evidence using khat and Ephedra species as model systems. Despite the commercial importance of amphetamine-type alkaloids, only a single step in their biosynthesis has been elucidated at the molecular level. We have employed Illumina next-generation sequencing technology, paired with Trinity and Velvet-Oases assembly platforms, to establish data-mining frameworks for Ephedra sinica and khat plants. Sequence libraries representing a combined 200,000 unigenes were subjected to an annotation pipeline involving direct searches against public databases. Annotations included the assignment of Gene Ontology (GO) terms used to allocate unigenes to functional categories. As part of our functional genomics program aimed at novel gene discovery, the databases were mined for enzyme candidates putatively involved in alkaloid biosynthesis. Queries used for mining included enzymes with established roles in benzoic acid metabolism, as well as enzymes catalyzing reactions similar to those predicted for amphetamine alkaloid metabolism. Gene candidates were evaluated based on phylogenetic relationships, FPKM-based expression data, and mechanistic considerations. Establishment of expansive sequence resources is a critical step toward pathway characterization, a goal with both academic and industrial implications.

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

PREMISE OF THE STUDY

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

METHODS

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

KEY RESULTS

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

CONCLUSIONS

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

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.

Phylogeny and divergence times of gymnosperms inferred from single-copy nuclear genes

Phylogenetic reconstruction is fundamental to study evolutionary biology and historical biogeography. However, there was not a molecular phylogeny of gymnosperms represented by extensive sampling at the genus level, and most published phylogenies of this group were constructed based on cytoplasmic DNA markers and/or the multi-copy nuclear ribosomal DNA. In this study, we use LFY and NLY, two single-copy nuclear genes that originated from an ancient gene duplication in the ancestor of seed plants, to reconstruct the phylogeny and estimate divergence times of gymnosperms based on a complete sampling of extant genera. The results indicate that the combined LFY and NLY coding sequences can resolve interfamilial relationships of gymnosperms and intergeneric relationships of most families. Moreover, the addition of intron sequences can improve the resolution in Podocarpaceae but not in cycads, although divergence times of the cycad genera are similar to or longer than those of the Podocarpaceae genera. Our study strongly supports cycads as the basal-most lineage of gymnosperms rather than sister to Ginkgoaceae, and a sister relationship between Podocarpaceae and Araucariaceae and between Cephalotaxaceae-Taxaceae and Cupressaceae. In addition, intergeneric relationships of some families that were controversial, and the relationships between Taxaceae and Cephalotaxaceae and between conifers and Gnetales are discussed based on the nuclear gene evidence. The molecular dating analysis suggests that drastic extinctions occurred in the early evolution of gymnosperms, and extant coniferous genera in the Northern Hemisphere are older than those in the Southern Hemisphere on average. This study provides an evolutionary framework for future studies on gymnosperms.