Loss of all plastid ndh genes in Gnetales and conifers: extent and evolutionary significance for the seed plant phylogeny

Alternative electron transport pathways in photosynthesis: a confluence of regulation

Photosynthetic reactions proceed along a linear electron transfer chain linking water oxidation at photosystem II (PSII) to CO₂ reduction in the Calvin-Benson-Bassham cycle. Alternative pathways poise the electron carriers along the chain in response to changing light, temperature and CO₂ inputs, under prolonged hydration stress and during development. We describe recent literature that reports the physiological functions of new molecular players. Such highlights include the flavodiiron proteins and their important role in the green lineage. The parsing of the proton-motive force between ΔpH and Δψ, regulated in many different ways (cyclic electron flow, ATPsynthase conductivity, ion/H⁺ transporters), is comprehensively reported. This review focuses on an integrated description of alternative electron transfer pathways and how they contribute to photosynthetic productivity in the context of plant fitness to the environment.

Concomitant loss of NDH complex-related genes within chloroplast and nuclear genomes in some orchids

The chloroplast NAD(P)H dehydrogenase-like (NDH) complex consists of about 30 subunits from both the nuclear and chloroplast genomes and is ubiquitous across most land plants. In some orchids, such as Phalaenopsis equestris, Dendrobium officinale and Dendrobium catenatum, most of the 11 chloroplast genome-encoded ndh genes (cp-ndh) have been lost. Here we investigated whether functional cp-ndh genes have been completely lost in these orchids or whether they have been transferred and retained in the nuclear genome. Further, we assessed whether both cp-ndh genes and nucleus-encoded NDH-related genes can be lost, resulting in the absence of the NDH complex. Comparative analyses of the genome of Apostasia odorata, an orchid species with a complete complement of cp-ndh genes which represents the sister lineage to all other orchids, and three published orchid genome sequences for P. equestris, D. officinale and D. catenatum, which are all missing cp-ndh genes, indicated that copies of cp-ndh genes are not present in any of these four nuclear genomes. This observation suggests that the NDH complex is not necessary for some plants. Comparative genomic/transcriptomic analyses of currently available plastid genome sequences and nuclear transcriptome data showed that 47 out of 660 photoautotrophic plants and all the heterotrophic plants are missing plastid-encoded cp-ndh genes and exhibit no evidence for maintenance of a functional NDH complex. Our data indicate that the NDH complex can be lost in photoautotrophic plant species. Further, the loss of the NDH complex may increase the probability of transition from a photoautotrophic to a heterotrophic life history.

Comparative analysis of Dendrobium plastomes and utility of plastomic mutational hotspots

Dendrobium is one of the largest genera in Orchidaceae, comprising about 800-1500 species mainly distributed in tropical Asia, Australasia, and Australia. There are 74 species and two varieties of this genus in China. Because of their ornamental and commercial value, Dendrobium orchids have been studied at low taxonomic levels. However, structural changes and effective mutational hotspots of Dendrobium plastomes have rarely been documented. Here, 30 Dendrobium plastomes were compared, comprising 25 newly sequenced in this study and five previously published. Except for their differences in NDH genes, these plastomes shared identical gene content and order. Comparative analyses revealed that the variation in size of Dendroubium plastomes was associated with dramatically changed length of InDels. Furthermore, ten loci were identified as the top-ten mutational hotspots, whose sequence variability was almost unchanged with more than 10 plastomes sampled, suggesting that they may be powerful markers for Dendrobium species. In addition, primer pairs of 47 polymorphic microsatellites were developed. After assessing the mean BS values of all combinations derived from the top-ten hotspots, we recommend that the combination of five hotspots-trnT-trnL, rpl32-trnL, clpP-psbB, trnL intron, and rps16-trnQ-should be used in the phylogenetic and identification studies of Dendrobium.

Comparative analysis of Dendrobium plastomes and utility of plastomic mutational hotspots

Dendrobium is one of the largest genera in Orchidaceae, comprising about 800-1500 species mainly distributed in tropical Asia, Australasia, and Australia. There are 74 species and two varieties of this genus in China. Because of their ornamental and commercial value, Dendrobium orchids have been studied at low taxonomic levels. However, structural changes and effective mutational hotspots of Dendrobium plastomes have rarely been documented. Here, 30 Dendrobium plastomes were compared, comprising 25 newly sequenced in this study and five previously published. Except for their differences in NDH genes, these plastomes shared identical gene content and order. Comparative analyses revealed that the variation in size of Dendroubium plastomes was associated with dramatically changed length of InDels. Furthermore, ten loci were identified as the top-ten mutational hotspots, whose sequence variability was almost unchanged with more than 10 plastomes sampled, suggesting that they may be powerful markers for Dendrobium species. In addition, primer pairs of 47 polymorphic microsatellites were developed. After assessing the mean BS values of all combinations derived from the top-ten hotspots, we recommend that the combination of five hotspots-trnT-trnL, rpl32-trnL, clpP-psbB, trnL intron, and rps16-trnQ-should be used in the phylogenetic and identification studies of Dendrobium.

The Complete Plastome Sequences of Four Orchid Species: Insights into the Evolution of the Orchidaceae and the Utility of Plastomic Mutational Hotspots

Orchidaceae (orchids) is the largest family in the monocots, including about 25,000 species in 880 genera and five subfamilies. Many orchids are highly valued for their beautiful and long-lasting flowers. However, the phylogenetic relationships among the five orchid subfamilies remain unresolved. The major dispute centers on whether the three one-stamened subfamilies, Epidendroideae, Orchidoideae, and Vanilloideae, are monophyletic or paraphyletic. Moreover, structural changes in the plastid genome (plastome) and the effective genetic loci at the species-level phylogenetics of orchids have rarely been documented. In this study, we compared 53 orchid plastomes, including four newly sequenced ones, that represent four remote genera: Dendrobium, Goodyera, Paphiopedilum, and Vanilla. These differ from one another not only in their lengths of inverted repeats and small single copy regions but also in their retention of ndh genes. Comparative analyses of the plastomes revealed that the expansion of inverted repeats in Paphiopedilum and Vanilla is associated with a loss of ndh genes. In orchid plastomes, mutational hotspots are genus specific. After having carefully examined the data, we propose that the three loci 5'trnK-rps16, trnS-trnG, and rps16-trnQ might be powerful markers for genera within Epidendroideae, and clpP-psbB and rps16-trnQ might be markers for genera within Cypripedioideae. After analyses of a partitioned dataset, we found that our plastid phylogenomic trees were congruent in a topology where two one-stamened subfamilies (i.e., Epidendroideae and Orchidoideae) were sisters to a multi-stamened subfamily (i.e., Cypripedioideae) rather than to the other one-stamened subfamily (Vanilloideae), suggesting that the living one-stamened orchids are paraphyletic.

Large-Scale Comparative Analysis Reveals the Mechanisms Driving Plastomic Compaction, Reduction, and Inversions in Conifers II (Cupressophytes)

Conifers II (cupressophytes), comprising about 400 tree species in five families, are the most diverse group of living gymnosperms. Their plastid genomes (plastomes) are highly variable in size and organization, but such variation has never been systematically studied. In this study, we assessed the potential mechanisms underlying the evolution of cupressophyte plastomes. We analyzed the plastomes of 24 representative genera in all of the five cupressophyte families, focusing on their variation in size, noncoding DNA content, and nucleotide substitution rates. Using a tree-based method, we further inferred the ancestral plastomic organizations of internal nodes and evaluated the inversions across the evolutionary history of cupressophytes. Our data showed that variation in plastome size is statistically associated with the dynamics of noncoding DNA content, which results in different degrees of plastomic compactness among the cupressophyte families. The degrees of plastomic inversions also vary among the families, with the number of inversions per genus ranging from 0 in Araucariaceae to 1.27 in Cupressaceae. In addition, we demonstrated that synonymous substitution rates are significantly correlated with plastome size as well as degree of inversions. These data suggest that in cupressophytes, mutation rates play a critical role in driving the evolution of plastomic size while plastomic inversions evolve in a neutral manner.

Development of the photosynthetic apparatus of Cunninghamia lanceolata in light and darkness

Here, we compared the development of dark- and light-grown Chinese fir (Cunninghamia lanceolata) cotyledons, which synthesize chlorophyll in the dark, representing a different phenomenon from angiosperm model plants. We determined that the grana lamellar membranes were well developed in both chloroplasts and etiochloroplasts. The accumulation of thylakoid membrane protein complexes was similar between chloroplasts and etiochloroplasts. Measurement of chlorophyll fluorescence parameters indicated that photosystem II (PSII) had low photosynthetic activities, whereas the photosystem I (PSI)-driven cyclic electron flow (CEF) rate exceeded the rate of PSII-mediated photon harvesting in etiochloroplasts. Analysis of the protein contents in etiochloroplasts indicated that the light-harvesting complex II remained mostly in its monomeric conformation. The ferredoxin NADP⁺ oxidoreductase and NADH dehydrogenase-like complexes were relatively abundantly expressed in etiochloroplasts for Chinese fir. Our transcriptome analysis contributes a global expression database for Chinese fir cotyledons, providing background information on the regulatory mechanisms of different genes involved in the development of dark- and light-grown cotyledons. In conclusion, we provide a novel description of the early developmental status of the light-dependent and light-independent photosynthetic apparatuses in gymnosperms.

The Chloroplast Genome of Utricularia reniformis Sheds Light on the Evolution of the ndh Gene Complex of Terrestrial Carnivorous Plants from the Lentibulariaceae Family

Lentibulariaceae is the richest family of carnivorous plants spanning three genera including Pinguicula, Genlisea, and Utricularia. Utricularia is globally distributed, and, unlike Pinguicula and Genlisea, has both aquatic and terrestrial forms. In this study we present the analysis of the chloroplast (cp) genome of the terrestrial Utricularia reniformis. U. reniformis has a standard cp genome of 139,725bp, encoding a gene repertoire similar to essentially all photosynthetic organisms. However, an exclusive combination of losses and pseudogenization of the plastid NAD(P)H-dehydrogenase (ndh) gene complex were observed. Comparisons among aquatic and terrestrial forms of Pinguicula, Genlisea, and Utricularia indicate that, whereas the aquatic forms retained functional copies of the eleven ndh genes, these have been lost or truncated in terrestrial forms, suggesting that the ndh function may be dispensable in terrestrial Lentibulariaceae. Phylogenetic scenarios of the ndh gene loss and recovery among Pinguicula, Genlisea, and Utricularia to the ancestral Lentibulariaceae cladeare proposed. Interestingly, RNAseq analysis evidenced that U. reniformis cp genes are transcribed, including the truncated ndh genes, suggesting that these are not completely inactivated. In addition, potential novel RNA-editing sites were identified in at least six U. reniformis cp genes, while none were identified in the truncated ndh genes. Moreover, phylogenomic analyses support that Lentibulariaceae is monophyletic, belonging to the higher core Lamiales clade, corroborating the hypothesis that the first Utricularia lineage emerged in terrestrial habitats and then evolved to epiphytic and aquatic forms. Furthermore, several truncated cp genes were found interspersed with U. reniformis mitochondrial and nuclear genome scaffolds, indicating that as observed in other smaller plant genomes, such as Arabidopsis thaliana, and the related and carnivorous Genlisea nigrocaulis and G. hispidula, the endosymbiotic gene transfer may also shape the U. reniformis genome in a similar fashion. Overall the comparative analysis of the U. reniformis cp genome provides new insight into the ndh genes and cp genome evolution of carnivorous plants from Lentibulariaceae family.

Revisiting the Plastid Phylogenomics of Pinaceae with Two Complete Plastomes of Pseudolarix and Tsuga

Phylogeny of the ten Pinaceous genera has long been contentious. Plastid genomes (plastomes) provide an opportunity to resolve this problem because they contain rich evolutionary information. To comprehend the plastid phylogenomics of all ten Pinaceous genera, we sequenced the plastomes of two previously unavailable genera, Pseudolarix amabilis (122,234 bp) and Tsuga chinensis (120,859 bp). Both plastomes share similar gene repertoire and order. Here for the first time we report a unique insertion of tandem repeats in accD of T. chinensis From the 65 plastid protein-coding genes common to all Pinaceous genera, we re-examined the phylogenetic relationship among all Pinaceous genera. Our two phylogenetic trees are congruent in an identical tree topology, with the five genera of the Abietoideae subfamily constituting a monophyletic clade separate from the other three subfamilies: Pinoideae, Piceoideae, and Laricoideae. The five genera of Abietoideae were grouped into two sister clades consisting of (1) Cedrus alone and (2) two sister subclades of Pseudolarix-Tsuga and Abies-Keteleeria, with the former uniquely losing the gene psaM and the latter specifically excluding the 3 psbA from the residual inverted repeat.

Chloroplast genomes: diversity, evolution, and applications in genetic engineering

Chloroplasts play a crucial role in sustaining life on earth. The availability of over 800 sequenced chloroplast genomes from a variety of land plants has enhanced our understanding of chloroplast biology, intracellular gene transfer, conservation, diversity, and the genetic basis by which chloroplast transgenes can be engineered to enhance plant agronomic traits or to produce high-value agricultural or biomedical products. In this review, we discuss the impact of chloroplast genome sequences on understanding the origins of economically important cultivated species and changes that have taken place during domestication. We also discuss the potential biotechnological applications of chloroplast genomes.

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.

NDH-1 and NDH-2 Plastoquinone Reductases in Oxygenic Photosynthesis

Oxygenic photosynthesis converts solar energy into chemical energy in the chloroplasts of plants and microalgae as well as in prokaryotic cyanobacteria using a complex machinery composed of two photosystems and both membrane-bound and soluble electron carriers. In addition to the major photosynthetic complexes photosystem II (PSII), cytochrome b6f, and photosystem I (PSI), chloroplasts also contain minor components, including a well-conserved type I NADH dehydrogenase (NDH-1) complex that functions in close relationship with photosynthesis and likewise originated from the endosymbiotic cyanobacterial ancestor. Some plants and many microalgal species have lost plastidial ndh genes and a functional NDH-1 complex during evolution, and studies have suggested that a plastidial type II NADH dehydrogenase (NDH-2) complex substitutes for the electron transport activity of NDH-1. However, although NDH-1 was initially thought to use NAD(P)H as an electron donor, recent research has demonstrated that both chloroplast and cyanobacterial NDH-1s oxidize reduced ferredoxin. We discuss more recent findings related to the biochemical composition and activity of NDH-1 and NDH-2 in relation to the physiology and regulation of photosynthesis, particularly focusing on their roles in cyclic electron flow around PSI, chlororespiration, and acclimation to changing environments.

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.

Plastid genomes reveal support for deep phylogenetic relationships and extensive rate variation among palms and other commelinid monocots

Despite progress based on multilocus, phylogenetic studies of the palms (order Arecales, family Arecaceae), uncertainty remains in resolution/support among major clades and for the placement of the palms among the commelinid monocots. Palms and related commelinids represent a classic case of substitution rate heterogeneity that has not been investigated in the genomic era. To address questions of relationships, support and rate variation among palms and commelinid relatives, 39 plastomes representing the palms and related family Dasypogonaceae were generated via genome skimming and integrated within a monocot-wide matrix for phylogenetic and molecular evolutionary analyses. Support was strong for 'deep' relationships among the commelinid orders, among the five palm subfamilies, and among tribes of the subfamily Coryphoideae. Additionally, there was extreme heterogeneity in the plastid substitution rates across the commelinid orders indicated by model based analyses, with c. 22 rate shifts, and significant departure from a global clock. To date, this study represents the most comprehensively sampled matrix of plastomes assembled for monocot angiosperms, providing genome-scale support for phylogenetic relationships of monocot angiosperms, and lays the phylogenetic groundwork for comparative analyses of the drivers and correlates of such drastic differences in substitution rates across a diverse and significant clade.

Seven New Complete Plastome Sequences Reveal Rampant Independent Loss of the ndh Gene Family across Orchids and Associated Instability of the Inverted Repeat/Small Single-Copy Region Boundaries

Earlier research has revealed that the ndh loci have been pseudogenized, truncated, or deleted from most orchid plastomes sequenced to date, including in all available plastomes of the two most species-rich subfamilies, Orchidoideae and Epidendroideae. This study sought to resolve deeper-level phylogenetic relationships among major orchid groups and to refine the history of gene loss in the ndh loci across orchids. The complete plastomes of seven orchids, Oncidium sphacelatum (Epidendroideae), Masdevallia coccinea (Epidendroideae), Sobralia callosa (Epidendroideae), Sobralia aff. bouchei (Epidendroideae), Elleanthus sodiroi (Epidendroideae), Paphiopedilum armeniacum (Cypripedioideae), and Phragmipedium longifolium (Cypripedioideae) were sequenced and analyzed in conjunction with all other available orchid and monocot plastomes. Most ndh loci were found to be pseudogenized or lost in Oncidium, Paphiopedilum and Phragmipedium, but surprisingly, all ndh loci were found to retain full, intact reading frames in Sobralia, Elleanthus and Masdevallia. Character mapping suggests that the ndh genes were present in the common ancestor of orchids but have experienced independent, significant losses at least eight times across four subfamilies. In addition, ndhF gene loss was correlated with shifts in the position of the junction of the inverted repeat (IR) and small single-copy (SSC) regions. The Orchidaceae have unprecedented levels of homoplasy in ndh gene presence/absence, which may be correlated in part with the unusual life history of orchids. These results also suggest that ndhF plays a role in IR/SSC junction stability.

Plastome Evolution in Hemiparasitic Mistletoes

Santalales is an order of plants consisting almost entirely of parasites. Some, such as Osyris, are facultative root parasites whereas others, such as Viscum, are obligate stem parasitic mistletoes. Here, we report the complete plastome sequences of one species of Osyris and three species of Viscum, and we investigate the evolutionary aspects of structural changes and changes in gene content in relation to parasitism. Compared with typical angiosperms plastomes, the four Santalales plastomes are all reduced in size (10–22% compared with Vitis), and they have experienced rearrangements, mostly but not exclusively in the border areas of the inverted repeats. Additionally, a number of protein-coding genes (matK, infA, ccsA, rpl33, and all 11 ndh genes) as well as two transfer RNA genes (trnG-UCC and trnV-UAC) have been pseudogenized or completely lost. Most of the remaining plastid genes have a significantly changed selection pattern compared with other dicots, and the relaxed selection of photosynthesis genes is noteworthy. Although gene loss obviously reduces plastome size, intergenic regions were also shortened. As plastome modifications are generally most prominent in Viscum, they are most likely correlated with the increased nutritional dependence on the host compared with Osyris.

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.

The location and translocation of ndh genes of chloroplast origin in the Orchidaceae family

The NAD(P)H dehydrogenase complex is encoded by 11 ndh genes in plant chloroplast (cp) genomes. However, ndh genes are truncated or deleted in some autotrophic Epidendroideae orchid cp genomes. To determine the evolutionary timing of the gene deletions and the genomic locations of the various ndh genes in orchids, the cp genomes of Vanilla planifolia, Paphiopedilum armeniacum, Paphiopedilum niveum, Cypripedium formosanum, Habenaria longidenticulata, Goodyera fumata and Masdevallia picturata were sequenced; these genomes represent Vanilloideae, Cypripedioideae, Orchidoideae and Epidendroideae subfamilies. Four orchid cp genome sequences were found to contain a complete set of ndh genes. In other genomes, ndh deletions did not correlate to known taxonomic or evolutionary relationships and deletions occurred independently after the orchid family split into different subfamilies. In orchids lacking cp encoded ndh genes, non cp localized ndh sequences were identified. In Erycina pusilla, at least 10 truncated ndh gene fragments were found transferred to the mitochondrial (mt) genome. The phenomenon of orchid ndh transfer to the mt genome existed in ndh-deleted orchids and also in ndh containing species.

The rise of the photosynthetic rate when light intensity increases is delayed in ndh gene-defective tobacco at high but not at low CO2 concentrations

The 11 plastid ndh genes have hovered frequently on the edge of dispensability, being absent in the plastid DNA of many algae and certain higher plants. We have compared the photosynthetic activity of tobacco (Nicotiana tabacum, cv. Petit Havana) with five transgenic lines (ΔndhF, pr-ΔndhF, T181D, T181A, and ndhF FC) and found that photosynthetic performance is impaired in transgenic ndhF-defective tobacco plants at rapidly fluctuating light intensities and higher than ambient CO2 concentrations. In contrast to wild type and ndhF FC, which reach the maximum photosynthetic rate in less than 1 min when light intensity suddenly increases, ndh defective plants (ΔndhF and T181A) show up to a 5 min delay in reaching the maximum photosynthetic rate at CO2 concentrations higher than the ambient 360 ppm. Net photosynthesis was determined at different CO2 concentrations when sequences of 130, 870, 61, 870, and 130 μmol m(-2) s(-1) PAR sudden light changes were applied to leaves and photosynthetic efficiency and entropy production (Sg) were determined as indicators of photosynthesis performance. The two ndh-defective plants, ΔndhF and T181A, had lower photosynthetic efficiency and higher Sg than wt, ndhF FC and T181D tobacco plants, containing full functional ndh genes, at CO2 concentrations above 400 ppm. We propose that the Ndh complex improves cyclic electron transport by adjusting the redox level of transporters during the low light intensity stage. In ndhF-defective strains, the supply of electrons through the Ndh complex fails, transporters remain over-oxidized (specially at high CO2 concentrations) and the rate of cyclic electron transport is low, impairing the ATP level required to rapidly reach high CO2 fixation rates in the following high light phase. Hence, ndh genes could be dispensable at low but not at high atmospheric concentrations of CO2.

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.

Investigating the path of plastid genome degradation in an early-transitional clade of heterotrophic orchids, and implications for heterotrophic angiosperms

Parasitic organisms exemplify morphological and genomic reduction. Some heterotrophic, parasitic plants harbor drastically reduced and degraded plastid genomes resulting from relaxed selective pressure on photosynthetic function. However, few studies have addressed the initial stages of plastome degradation in groups containing both photosynthetic and nonphotosynthetic species. Corallorhiza is a genus of leafless, heterotrophic orchids that contains both green, photosynthetic species and nongreen, putatively nonphotosynthetic species, and represents an ideal system in which to assess the beginning of the transition to a "minimal plastome." Complete plastomes were generated for nine taxa of Corallorhiza using Illumina paired-end sequencing of genomic DNA to assess the degree of degradation among taxa, and for comparison with a general model of degradation among angiosperms. Quantification of total chlorophyll suggests that nongreen Corallorhiza still produce chlorophyll, but at 10-fold lower concentrations than green congeners. Complete plastomes and partial nuclear rDNA cistrons yielded a fully resolved tree for Corallorhiza, with at least two independent losses of photosynthesis, evidenced by gene deletions and pseudogenes in Co. striata and nongreen Co. maculata. All Corallorhiza show some evidence of degradation in genes of the NAD(P)H dehydrogenase complex. Among genes with open reading frames, photosynthesis-related genes displayed evidence of neutral evolution in nongreen Corallorhiza, whereas genes of the ATP synthase complex displayed some evidence of positive selection in these same groups, though for reasons unknown. Corallorhiza spans the early stages of a general model of plastome degradation and has added critical insight for understanding the process of plastome evolution in heterotrophic angiosperms.

Disproportional plastome-wide increase of substitution rates and relaxed purifying selection in genes of carnivorous Lentibulariaceae

Carnivorous Lentibulariaceae exhibit the most sophisticated implementation of the carnivorous syndrome in plants. Their unusual lifestyle coincides with distinct genomic peculiarities such as the smallest angiosperm nuclear genomes and extremely high nucleotide substitution rates across all genomic compartments. Here, we report the complete plastid genomes from each of the three genera Pinguicula, Utricularia, and Genlisea, and investigate plastome-wide changes in their molecular evolution as the carnivorous syndrome unfolds. We observe a size reduction by up to 9% mostly due to the independent loss of genes for the plastid NAD(P)H dehydrogenase and altered proportions of plastid repeat DNA, as well as a significant plastome-wide increase of substitution rates and microstructural changes. Protein-coding genes across all gene classes show a disproportional elevation of nonsynonymous substitutions, particularly in Utricularia and Genlisea. Significant relaxation of purifying selection relative to noncarnivores occurs in the plastid-encoded fraction of the photosynthesis ATP synthase complex, the photosystem I, and in several other photosynthesis and metabolic genes. Shifts in selective regimes also affect housekeeping genes including the plastid-encoded polymerase, for which evidence for relaxed purifying selection was found once during the transition to carnivory, and a second time during the diversification of the family. Lentibulariaceae significantly exhibit enhanced rates of nucleotide substitution in most of the 130 noncoding regions. Various factors may underlie the observed patterns of relaxation of purifying selection and substitution rate increases, such as reduced net photosynthesis rates, alternative paths of nutrient uptake (including organic carbon), and impaired DNA repair mechanisms.

The plastid genome of Najas flexilis: adaptation to submersed environments is accompanied by the complete loss of the NDH complex in an aquatic angiosperm

The re-colonization of aquatic habitats by angiosperms has presented a difficult challenge to plants whose long evolutionary history primarily reflects adaptations to terrestrial conditions. Many aquatics must complete vital stages of their life cycle on the water surface by means of floating or emergent leaves and flowers. Only a few species, mainly within the order Alismatales, are able to complete all aspects of their life cycle including pollination, entirely underwater. Water-pollinated Alismatales include seagrasses and water nymphs (Najas), the latter being the only freshwater genus in the family Hydrocharitaceae with subsurface water-pollination. We have determined the complete nucleotide sequence of the plastid genome of Najas flexilis. The plastid genome of N. flexilis is a circular AT-rich DNA molecule of 156 kb, which displays a quadripartite structure with two inverted repeats (IR) separating the large single copy (LSC) from the small single copy (SSC) regions. In N. flexilis, as in other Alismatales, the rps19 and trnH genes are localized in the LSC region instead of within the IR regions as in other monocots. However, the N. flexilis plastid genome presents some anomalous modifications. The size of the SSC region is only one third of that reported for closely related species. The number of genes in the plastid is considerably less. Both features are due to loss of the eleven ndh genes in the Najas flexilis plastid. In angiosperms, the absence of ndh genes has been related mainly to the loss of photosynthetic function in parasitic plants. The ndh genes encode the NAD(P)H dehydrogenase complex, believed essential in terrestrial environments, where it increases photosynthetic efficiency in variable light intensities. The modified structure of the N. flexilis plastid genome suggests that adaptation to submersed environments, where light is scarce, has involved the loss of the NDH complex in at least some photosynthetic angiosperms.

Chengia laxispicata gen. et sp. nov., a new ephedroid plant from the Early Cretaceous Yixian Formation of western Liaoning, Northeast China: evolutionary, taxonomic, and biogeographic implications

Background

The extant Gnetales include three monotypic families, namely, Ephedraceae (Ephedra), Gnetaceae (Gnetum), and Welwitschiaceae (Welwitschia), all of which possess compound female cones that comprise a main axis and 1 to multiple pairs/whorls of bracts subtending a female reproductive unit or having lower pairs/whorls of bracts sterile. However, the evolutionary origin of such a reproductive architecture in Gnetales is controversial in the light of the competing anthophyte versus gnetifer hypotheses of seed plant relationships. Hence, macrofossils demonstrating the structure of compound female cones of the Gnetales should be important to decipher the early evolution of the order.

Results

A new ephedroid plant Chengia laxispicata gen. et sp. nov. is described from the Early Cretaceous Yixian Formation of western Liaoning, Northeast China. The fossil represents a part of a leafy shooting system with reproductive organs attached. The main shoot bears internodes and swollen nodes, from which lateral branches arise oppositely. Reproductive organs consist of female spikes terminal to twigs or axillary to linear leaves. Spikes are loosely arranged, having prominent nodes and internodes. Bracts of the spikes are decussately opposite and comprise 4—8 pairs of bracts. Each bract subtends an ellipsoid seed. Seeds are sessile, with a thin outer envelope and a distal micropylar tube.

Conclusions

Chengia laxispicata gen. et sp. nov. provides a missing link between archetypal fertile organs in the crown lineage of the Gnetales and compound female cones of the extant Ephedraceae. Combined with a wealth of Ephedra and ephedroid macrofossils from the Early Cretaceous, we propose a reduction and sterilization hypothesis that the female cone of the extant Ephedraceae may have stemmed from archetypal fertile organs in the crown lineage of the Gnetales. These have undergone sequentially intermediate links similar to female cones of Cretaceous Siphonospermum, Chengia, and Liaoxia by reduction and sterilization of the lower fertile bracts, shortenings of internodes and peduncles as well as loss of reproductive units in all inferior bracts. The basal family Ephedraceae including Ephedra of the extant Gnetales was demonstrated to have considerable diversity by the Early Cretaceous, so an emended familial diagnosis is given here. The Jehol Biota in Northeast China and adjacent areas contains a plethora of well-preserved macrofossils of Ephedra and ephedroids that show different evolutionary stages including primitive and derived characters of Ephedraceae, so Northeast China and adjacent areas may represent either the centre of origination or one of the centres for early diversification of the family.

Phylogenetic viewpoints on regulation of light harvesting and electron transport in eukaryotic photosynthetic organisms

The comparative study of photosynthetic regulation in the thylakoid membrane of different phylogenetic groups can yield valuable insights into mechanisms, genetic requirements and redundancy of regulatory processes. This review offers a brief summary on the current understanding of light harvesting and photosynthetic electron transport regulation in different photosynthetic eukaryotes, with a special focus on the comparison between higher plants and unicellular algae of secondary endosymbiotic origin. The foundations of thylakoid structure, light harvesting, reversible protein phosphorylation and PSI-mediated cyclic electron transport are traced not only from green algae to vascular plants but also at the branching point between the "green" and the "red" lineage of photosynthetic organisms. This approach was particularly valuable in revealing processes that (1) are highly conserved between phylogenetic groups, (2) serve a common physiological role but nevertheless originate in divergent genetic backgrounds or (3) are missing in one phylogenetic branch despite their unequivocal importance in another, necessitating a search for alternative regulatory mechanisms and interactions.

Plastid genome evolution across the genus Cuscuta (Convolvulaceae): two clades within subgenus Grammica exhibit extensive gene loss

The genus Cuscuta (Convolvulaceae, the morning glory family) is one of the most intensely studied lineages of parasitic plants. Whole plastome sequencing of four Cuscuta species has demonstrated changes to both plastid gene content and structure. The presence of photosynthetic genes under purifying selection indicates that Cuscuta is cryptically photosynthetic. However, the tempo and mode of plastid genome evolution across the diversity of this group (~200 species) remain largely unknown. A comparative investigation of plastid genome content, grounded within a phylogenetic framework, was conducted using a slot-blot Southern hybridization approach. Cuscuta was extensively sampled (~56% of species), including groups previously suggested to possess more altered plastomes compared with other members of this genus. A total of 56 probes derived from all categories of protein-coding genes, typically found within the plastomes of flowering plants, were used. The results indicate that two clades within subgenus Grammica (clades 'O' and 'K') exhibit substantially more plastid gene loss relative to other members of Cuscuta. All surveyed members of the 'O' clade show extensive losses of plastid genes from every category of genes typically found in the plastome, including otherwise highly conserved small and large ribosomal subunits. The extent of plastid gene losses within this clade is similar in magnitude to that observed previously in some non-asterid holoparasites, in which the very presence of a plastome has been questioned. The 'K' clade also exhibits considerable loss of plastid genes. Unlike in the 'O' clade, in which all species seem to be affected, the losses in clade 'K' progress phylogenetically, following a pattern consistent with the Evolutionary Transition Series hypothesis. This clade presents an ideal opportunity to study the reduction of the plastome of parasites 'in action'. The widespread plastid gene loss in these two clades is hypothesized to be a consequence of the complete loss of photosynthesis. Additionally, taxa that would be the best candidates for entire plastome sequencing are identified in order to investigate further the loss of photosynthesis and reduction of the plastome within Cuscuta.

Complete plastid genomes from Ophioglossum californicum, Psilotum nudum, and Equisetum hyemale reveal an ancestral land plant genome structure and resolve the position of Equisetales among monilophytes

BACKGROUND

Plastid genome structure and content is remarkably conserved in land plants. This widespread conservation has facilitated taxon-rich phylogenetic analyses that have resolved organismal relationships among many land plant groups. However, the relationships among major fern lineages, especially the placement of Equisetales, remain enigmatic.

RESULTS

In order to understand the evolution of plastid genomes and to establish phylogenetic relationships among ferns, we sequenced the plastid genomes from three early diverging species: Equisetum hyemale (Equisetales), Ophioglossum californicum (Ophioglossales), and Psilotum nudum (Psilotales). A comparison of fern plastid genomes showed that some lineages have retained inverted repeat (IR) boundaries originating from the common ancestor of land plants, while other lineages have experienced multiple IR changes including expansions and inversions. Genome content has remained stable throughout ferns, except for a few lineage-specific losses of genes and introns. Notably, the losses of the rps16 gene and the rps12i346 intron are shared among Psilotales, Ophioglossales, and Equisetales, while the gain of a mitochondrial atp1 intron is shared between Marattiales and Polypodiopsida. These genomic structural changes support the placement of Equisetales as sister to Ophioglossales + Psilotales and Marattiales as sister to Polypodiopsida. This result is augmented by some molecular phylogenetic analyses that recover the same relationships, whereas others suggest a relationship between Equisetales and Polypodiopsida.

CONCLUSIONS

Although molecular analyses were inconsistent with respect to the position of Marattiales and Equisetales, several genomic structural changes have for the first time provided a clear placement of these lineages within the ferns. These results further demonstrate the power of using rare genomic structural changes in cases where molecular data fail to provide strong phylogenetic resolution.

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.

Does the mode of plastid inheritance influence plastid genome architecture?

Plastid genomes show an impressive array of sizes and compactnesses, but the forces responsible for this variation are unknown. It has been argued that species with small effective genetic population sizes are less efficient at purging excess DNA from their genomes than those with large effective population sizes. If true, one may expect the primary mode of plastid inheritance to influence plastid DNA (ptDNA) architecture. All else being equal, biparentally inherited ptDNAs should have a two-fold greater effective population size than those that are uniparentally inherited, and thus should also be more compact. Here, we explore the relationship between plastid inheritance pattern and ptDNA architecture, and consider the role of phylogeny in shaping our observations. Contrary to our expectations, we found no significant difference in plastid genome size or compactness between ptDNAs that are biparentally inherited relative to those that are uniparentally inherited. However, we also found that there was significant phylogenetic signal for the trait of mode of plastid inheritance. We also found that paternally inherited ptDNAs are significantly smaller (n = 19, p = 0.000001) than those that are maternally, uniparentally (when isogamous), or biparentally inherited. Potential explanations for this observation are discussed.

Activity, polypeptide and gene identification of thylakoid Ndh complex in trees: potential physiological relevance of fluorescence assays

Three evergreen (Laurus nobilis, Viburnum tinus and Thuja plicata) and two autumnal abscission deciduous trees (Cydonia oblonga and Prunus domestica) have been investigated for the presence (zymogram and immunodetection) and functionality (post-illumination chlorophyll fluorescence) of the thylakoid Ndh complex. The presence of encoding ndh genes has also been investigated in T. plicata. Western assays allowed tentative identification of zymogram NADH dehydrogenase bands corresponding to the Ndh complex after native electrophoresis of solubilized fractions from L. nobilis, V. tinus, C. oblonga and P. domestica leaves, but not in those of T. plicata. However, Ndh subunits were detected after SDS-PAGE of thylakoid solubilized proteins of T. plicata. The leaves of the five plants showed the post-illumination chlorophyll fluorescence increase dependent on the presence of active Ndh complex. The fluorescence increase was higher in autumn in deciduous, but not in evergreen trees, which suggests that the thylakoid Ndh complex could be involved in autumnal leaf senescence. Two ndhB genes were sequenced from T. plicata that differ at the 350 bp 3' end sequence. Comparison with the mRNA revealed that ndhB genes have a 707-bp type II intron between exons 1 (723 bp) and 2 (729 bp) and that the UCA 259th codon is edited to UUA in mRNA. Phylogenetically, the ndhB genes of T. plicata group close to those of Metasequoia, Cryptomeria, Taxodium, Juniperus and Widdringtonia in the cupresaceae branch and are 5' end shortened by 18 codons with respect to that of angiosperms.

The plastid genome of the mycoheterotrophic Corallorhiza striata (Orchidaceae) is in the relatively early stages of degradation

PREMISE OF THE STUDY

Plastid genomes of nonphotosynthetic, mycoheterotrophic plants represent apt systems in which to study effects of relaxed evolutionary constraints. The few mycoheterotrophic angiosperm plastomes sequenced to date display drastic patterns of degradation/reduction relative to those of photosynthetic relatives. The goal of this study was to focus on a mycoheterotrophic orchid hypothesized to be in the "early" stages of plastome degradation, to provide perspective on this process.

METHODS

Short-read sequencing was used to generate a complete plastome sequence for Corallorhiza striata var. vreelandii, a mycoheterotrophic orchid, to investigate the extent of plastome degradation. Patterns of nonsynonymous/synonymous mutations were also assessed, and comparisons were made between Corallorhiza and other heterotrophic plant lineages.

KEY RESULTS

Corallorhiza yielded a plastome of 137505 bp, with several photosynthesis-related genes either lost or pseudogenized. Members of all major photosynthesis complexes, except ATP-synthase genes, were affected. "Housekeeping" genes were intact, despite the loss of a single tRNA. Intact photosynthesis genes (excluding atp genes) together displayed elevated nonsynonymous changes, while housekeeping genes did not.

CONCLUSIONS

The Corallorhiza plastome is not drastically reduced in overall size (∼6% reduction relative to that of photosynthetic Oncidium), but displays a pattern congruent with a loss of photosynthetic function. Comparing Corallorhiza with other heterotrophs allows some emergent evolutionary patterns to be inferred, but these remain as hypotheses to be tested, especially at lower taxonomic levels, and in lineages illustrating transitions from autotrophy to heterotrophy. The independent, unique processes of plastome modification among mycoheterotrophic lineages illustrate the urgency of their conservation.

The evolutionary root of flowering plants

Correct rooting of the angiosperm radiation is both challenging and necessary for understanding the origins and evolution of physiological and phenotypic traits in flowering plants. The problem is known to be difficult due to the large genetic distance separating flowering plants from other seed plants and the sparse taxon sampling among basal angiosperms. Here, we provide further evidence for concern over substitution model misspecification in analyses of chloroplast DNA sequences. We show that support for Amborella as the sole representative of the most basal angiosperm lineage is founded on sequence site patterns poorly described by time-reversible substitution models. Improving the fit between sequence data and substitution model identifies Trithuria, Nymphaeaceae, and Amborella as surviving relatives of the most basal lineage of flowering plants. This finding indicates that aquatic and herbaceous species dominate the earliest extant lineage of flowering plants. [; ; ; ; ; .].

The comparative chloroplast genomic analysis of photosynthetic orchids and developing DNA markers to distinguish Phalaenopsis orchids

The chloroplast genome of Phalaenopsis equestris was determined and compared to those of Phalaenopsis aphrodite and Oncidium Gower Ramsey in Orchidaceae. The chloroplast genome of P. equestris is 148,959 bp, and a pair of inverted repeats (25,846 bp) separates the genome into large single-copy (85,967 bp) and small single-copy (11,300 bp) regions. The genome encodes 109 genes, including 4 rRNA, 30 tRNA and 75 protein-coding genes, but loses four ndh genes (ndhA, E, F and H) and seven other ndh genes are pseudogenes. The rate of inter-species variation between the two moth orchids was 0.74% (1107 sites) for single nucleotide substitution and 0.24% for insertions (161 sites; 1388 bp) and deletions (189 sites; 1393 bp). The IR regions have a lower rate of nucleotide substitution (3.5-5.8-fold) and indels (4.3-7.1-fold) than single-copy regions. The intergenic spacers are the most divergent, and based on the length variation of the three intergenic spacers, 11 native Phalaenopsis orchids could be successfully distinguished. The coding genes, IR junction and RNA editing sites are relatively more conserved between the two moth orchids than between those of Phalaenopsis and Oncidium spp.

Plastid genome evolution in mycoheterotrophic Ericaceae

Unlike parasitic plants, which are linked to their hosts directly through haustoria, mycoheterotrophic (MHT) plants derive all or part of their water and nutrients from autothrophs via fungal mycorrhizal intermediaries. Ericaceae, the heather family, are a large and diverse group of plants known to form elaborate symbiotic relationships with mycorrhizal fungi. Using PHYA sequence data, we first investigated relationships among mycoheterotrophic Ericaceae and their close autotrophic relatives. Phylogenetic results suggest a minimum of two independent origins of MHT within this family. Additionally, a comparative investigation of plastid genomes (plastomes) grounded within this phylogenetic framework was conducted using a slot-blot Southern hybridization approach. This survey encompassed numerous lineages of Ericaceae with different life histories and trophic levels, including multiple representatives from mixotrophic Pyroleae and fully heterotrophic Monotropeae and Pterosporeae. Fifty-four probes derived from all categories of protein coding genes typically found within the plastomes of flowering plants were used. Our results indicate that the holo-mycoheterotrophic Ericaceae exhibit extensive loss of genes relating to photosynthetic function and expression of the plastome but retain genes with possible functions outside photosynthesis. Mixotrophic taxa tend to retain most genes relating to photosynthetic functions but are varied regarding the plastid ndh gene content. This investigation extends previous inferences that the loss of the NDH complex occurs prior to becoming holo-heterotrophic and it shows that the pattern of gene losses among mycoheterotrophic Ericaceae is similar to that of haustorial parasites. Additionally, we identify the most desirable candidate species for entire plastome sequencing.

A functional phylogenomic view of the seed plants

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

Sequencing and analysis of plastid genome in mycoheterotrophic orchid Neottia nidus-avis

Plastids are the semiautonomous organelles that possess their own genome inherited from the cyanobacterial ancestor. The primary function of plastids is photosynthesis so the structure and evolution of plastid genomes are extensively studied in photosynthetic plants. In contrast, little is known about the plastomes of nonphotosynthetic species. In higher plants, plastid genome sequences are available for only three strictly nonphotosynthetic species, the liverwort Aneura mirabilis and two flowering plants, Epifagus virginiana and Rhizanthella gardneri. We report here the complete sequence of a plastid genome of nonphotosynthetic mycoheterotrophic orchid Neottia nidus-avis, determined using 454 pyrosequencing technology. It was found to be reduced in both genome size and gene content; this reduction is however not as drastic as in the other nonphotosynthetic orchid, R. gardneri. Neottia plastome lacks all genes encoding photosynthetic proteins, RNA polymerase subunits but retains most genes of translational apparatus. Those genes that are retained have an increased rate of both synonymous and nonsynonymous substitutions but do not exhibit relaxation of purifying selection either in Neottia or in Rhizanthella.

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

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

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

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

Southward Pleistocene migration of Douglas-fir into Mexico: phylogeography, ecological niche modeling, and conservation of 'rear edge' populations

• Poleward Pleistocene plant migration has been an important process structuring modern temperate and boreal plant communities, but the contribution of equatorward migration remains poorly understood. Paleobotanical evidence suggests Miocene or Pleistocene origin for temperate 'sky island' plant taxa in Mexico. These 'rear edge' populations situated in a biodiversity hotspot may be an important reserve of genetic diversity in changing climates. • We used mtDNA sequences, cpDNA sequences and chloroplast microsatellites to test hypotheses of Miocene vs Pleistocene colonization of temperate Douglas-fir in Mexico, explore geographic patterns of molecular variation in relation to Pleistocene climate history using ecological niche models, and assess the taxonomic and conservation implications. • We found strong evidence for Pleistocene divergence of Douglas-fir in Mexico (958 thousand yr before present (ka) with the 90% highest posterior density interval ranging from 1.6 million yr before present (Ma) to 491 ka), consistent with the southward Pleistocene migration hypothesis. Genetic diversity was high and strongly partitioned among populations. Spatial patterns of molecular variation and ecological niche models suggest a complex late Pleistocene history involving periods of isolation and expansion along mountain corridors. • These results highlight the importance of southward Pleistocene migration in establishing modern high-diversity plant communities and provide critical insights into proposals to conserve the unique biodiversity of Mexican Douglas-fir and associated taxa.

An ILP solution for the gene duplication problem

Background

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

Results

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

Conclusions

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

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.

Plastid ndh genes in plant evolution

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

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

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

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

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