Phylotranscriptomic Analyses of Mycoheterotrophic Monocots Show a Continuum of Convergent Evolutionary Changes in Expressed Nuclear Genes From Three Independent Nonphotosynthetic Lineages.

Mycoheterotrophy is an alternative nutritional strategy whereby plants obtain sugars and other nutrients from soil fungi. Mycoheterotrophy and associated loss of photosynthesis have evolved repeatedly in plants, particularly in monocots. Although reductive evolution of plastomes in mycoheterotrophs is well documented, the dynamics of nuclear genome evolution remains largely unknown. Transcriptome datasets were generated from four mycoheterotrophs in three families (Orchidaceae, Burmanniaceae, Triuridaceae) and related green plants and used for phylogenomic analyses to resolve relationships among the mycoheterotrophs, their relatives, and representatives across the monocots. Phylogenetic trees based on 602 genes were mostly congruent with plastome phylogenies, except for an Asparagales + Liliales clade inferred in the nuclear trees. Reduction and loss of chlorophyll synthesis and photosynthetic gene expression and relaxation of purifying selection on retained genes were progressive, with greater loss in older nonphotosynthetic lineages. One hundred seventy-four of 1375 plant benchmark universally conserved orthologous genes were undetected in any mycoheterotroph transcriptome or the genome of the mycoheterotrophic orchid Gastrodia but were expressed in green relatives, providing evidence for massively convergent gene loss in nonphotosynthetic lineages. We designate this set of deleted or undetected genes Missing in Mycoheterotrophs (MIM). MIM genes encode not only mainly photosynthetic or plastid membrane proteins but also a diverse set of plastid processes, genes of unknown function, mitochondrial, and cellular processes. Transcription of a photosystem II gene (psb29) in all lineages implies a nonphotosynthetic function for this and other genes retained in mycoheterotrophs. Nonphotosynthetic plants enable novel insights into gene function as well as gene expression shifts, gene loss, and convergence in nuclear genomes.

Phylogenomic resolution of order- and family-level monocot relationships using 602 single-copy nuclear genes and 1375 BUSCO genes.

We assess relationships among 192 species in all 12 monocot orders and 72 of 77 families, using 602 conserved single-copy (CSC) genes and 1375 benchmarking single-copy ortholog (BUSCO) genes extracted from genomic and transcriptomic datasets. Phylogenomic inferences based on these data, using both coalescent-based and supermatrix analyses, are largely congruent with the most comprehensive plastome-based analysis, and nuclear-gene phylogenomic analyses with less comprehensive taxon sampling. The strongest discordance between the plastome and nuclear gene analyses is the monophyly of a clade comprising Asparagales and Liliales in our nuclear gene analyses, versus the placement of Asparagales and Liliales as successive sister clades to the commelinids in the plastome tree. Within orders, around six of 72 families shifted positions relative to the recent plastome analysis, but four of these involve poorly supported inferred relationships in the plastome-based tree. In Poales, the nuclear data place a clade comprising Ecdeiocoleaceae+Joinvilleaceae as sister to the grasses (Poaceae); Typhaceae, (rather than Bromeliaceae) are resolved as sister to all other Poales. In Commelinales, nuclear data place Philydraceae sister to all other families rather than to a clade comprising Haemodoraceae+Pontederiaceae as seen in the plastome tree. In Liliales, nuclear data place Liliaceae sister to Smilacaceae, and Melanthiaceae are placed sister to all other Liliales except Campynemataceae. Finally, in Alismatales, nuclear data strongly place Tofieldiaceae, rather than Araceae, as sister to all the other families, providing an alternative resolution of what has been the most problematic node to resolve using plastid data, outside of those involving achlorophyllous mycoheterotrophs. As seen in numerous prior studies, the placement of orders Acorales and Alismatales as successive sister lineages to all other extant monocots. Only 21.2% of BUSCO genes were demonstrably single-copy, yet phylogenomic inferences based on BUSCO and CSC genes did not differ, and overall functional annotations of the two sets were very similar. Our analyses also reveal significant gene tree-species tree discordance despite high support values, as expected given incomplete lineage sorting (ILS) related to rapid diversification. Our study advances understanding of monocot relationships and the robustness of phylogenetic inferences based on large numbers of nuclear single-copy genes that can be obtained from transcriptomes and genomes.

Dynamic genome evolution in a model fern.

The large size and complexity of most fern genomes have hampered efforts to elucidate fundamental aspects of fern biology and land plant evolution through genome-enabled research. Here we present a chromosomal genome assembly and associated methylome, transcriptome and metabolome analyses for the model fern species Ceratopteris richardii. The assembly reveals a history of remarkably dynamic genome evolution including rapid changes in genome content and structure following the most recent whole-genome duplication approximately 60 million years ago. These changes include massive gene loss, rampant tandem duplications and multiple horizontal gene transfers from bacteria, contributing to the diversification of defence-related gene families. The insertion of transposable elements into introns has led to the large size of the Ceratopteris genome and to exceptionally long genes relative to other plants. Gene family analyses indicate that genes directing seed development were co-opted from those controlling the development of fern sporangia, providing insights into seed plant evolution. Our findings and annotated genome assembly extend the utility of Ceratopteris as a model for investigating and teaching plant biology.

Plant ecological genomics at the limits of life in the Atacama Desert.

The Atacama Desert in Chile-hyperarid and with high-ultraviolet irradiance levels-is one of the harshest environments on Earth. Yet, dozens of species grow there, including Atacama-endemic plants. Herein, we establish the Talabre-Lejía transect (TLT) in the Atacama as an unparalleled natural laboratory to study plant adaptation to extreme environmental conditions. We characterized climate, soil, plant, and soil-microbe diversity at 22 sites (every 100 m of altitude) along the TLT over a 10-y period. We quantified drought, nutrient deficiencies, large diurnal temperature oscillations, and pH gradients that define three distinct vegetational belts along the altitudinal cline. We deep-sequenced transcriptomes of 32 dominant plant species spanning the major plant clades, and assessed soil microbes by metabarcoding sequencing. The top-expressed genes in the 32 Atacama species are enriched in stress responses, metabolism, and energy production. Moreover, their root-associated soils are enriched in growth-promoting bacteria, including nitrogen fixers. To identify genes associated with plant adaptation to harsh environments, we compared 32 Atacama species with the 32 closest sequenced species, comprising 70 taxa and 1,686,950 proteins. To perform phylogenomic reconstruction, we concatenated 15,972 ortholog groups into a supermatrix of 8,599,764 amino acids. Using two codon-based methods, we identified 265 candidate positively selected genes (PSGs) in the Atacama plants, 64% of which are located in Pfam domains, supporting their functional relevance. For 59/184 PSGs with an Arabidopsis ortholog, we uncovered functional evidence linking them to plant resilience. As some Atacama plants are closely related to staple crops, these candidate PSGs are a "genetic goldmine" to engineer crop resilience to face climate change.

Cycad-Weevil Pollination Symbiosis Is Characterized by Rapidly Evolving and Highly Specific Plant-Insect Chemical Communication.

Coevolution between plants and insects is thought to be responsible for generating biodiversity. Extensive research has focused largely on antagonistic herbivorous relationships, but mutualistic pollination systems also likely contribute to diversification. Here we describe an example of chemically-mediated mutualistic species interactions affecting trait evolution and lineage diversification. We show that volatile compounds produced by closely related species of Zamia cycads are more strikingly different from each other than are other phenotypic characters, and that two distantly related pollinating weevil species have specialized responses only to volatiles from their specific host Zamia species. Plant transcriptomes show that approximately a fifth of genes related to volatile production are evolving under positive selection, but we find no differences in the relative proportion of genes under positive selection in different categories. The importance of phenotypic divergence coupled with chemical communication for the maintenance of this obligate mutualism highlights chemical signaling as a key mechanism of coevolution between cycads and their weevil pollinators.

Flower and Spikelet Construction in Rapateaceae (Poales).

The family Rapateaceae represents an early-divergent lineage of Poales with biotically pollinated showy flowers. We investigate developmental morphology and anatomy in all three subfamilies and five tribes of Rapateaceae to distinguish between contrasting hypotheses on spikelet morphology and to address questions on the presence of nectaries and gynoecium structure. We support an interpretation of the partial inflorescence (commonly termed spikelet), as a uniaxial system composed of a terminal flower and numerous empty phyllomes. A terminal flower in an inflorescence unit is an autapomorphic feature of Rapateaceae. The gynoecium consists of synascidiate, symplicate, and usually asymplicate zones, with gynoecium formation encompassing congenital and often also postgenital fusions between carpels. Species of Rapateaceae differ in the relative lengths of the gynoecial zones, the presence or absence of postgenital fusion between the carpels and placentation in the ascidiate or plicate carpel zones. In contrast with previous reports, septal nectaries are lacking in all species. The bird-pollinated tribe Schoenocephalieae is characterized by congenital syncarpy; it displays an unusual type of gynoecial (non-septal) nectary represented by a secretory epidermis at the gynoecium base.

Organellomic data sets confirm a cryptic consensus on (unrooted) land-plant relationships and provide new insights into bryophyte molecular evolution.

PREMISE: Phylogenetic trees of bryophytes provide important evolutionary context for land plants. However, published inferences of overall embryophyte relationships vary considerably. We performed phylogenomic analyses of bryophytes and relatives using both mitochondrial and plastid gene sets, and investigated bryophyte plastome evolution. METHODS: We employed diverse likelihood-based analyses to infer large-scale bryophyte phylogeny for mitochondrial and plastid data sets. We tested for changes in purifying selection in plastid genes of a mycoheterotrophic liverwort (Aneura mirabilis) and a putatively mycoheterotrophic moss (Buxbaumia), and compared 15 bryophyte plastomes for major structural rearrangements. RESULTS: Overall land-plant relationships conflict across analyses, generally weakly. However, an underlying (unrooted) four-taxon tree is consistent across most analyses and published studies. Despite gene coverage patchiness, relationships within mosses, liverworts, and hornworts are largely congruent with previous studies, with plastid results generally better supported. Exclusion of RNA edit sites restores cases of unexpected non-monophyly to monophyly for Takakia and two hornwort genera. Relaxed purifying selection affects multiple plastid genes in mycoheterotrophic Aneura but not Buxbaumia. Plastid genome structure is nearly invariant across bryophytes, but the tufA locus, presumed lost in embryophytes, is unexpectedly retained in several mosses. CONCLUSIONS: A common unrooted tree underlies embryophyte phylogeny, [(liverworts, mosses), (hornworts, vascular plants)]; rooting inconsistency across studies likely reflects substantial distance to algal outgroups. Analyses combining genomic and transcriptomic data may be misled locally for heavily RNA-edited taxa. The Buxbaumia plastome lacks hallmarks of relaxed selection found in mycoheterotrophic Aneura. Autotrophic bryophyte plastomes, including Buxbaumia, hardly vary in overall structure.

The Plant Ontology Facilitates Comparisons of Plant Development Stages Across Species.

The Plant Ontology (PO) is a community resource consisting of standardized terms, definitions, and logical relations describing plant structures and development stages, augmented by a large database of annotations from genomic and phenomic studies. This paper describes the structure of the ontology and the design principles we used in constructing PO terms for plant development stages. It also provides details of the methodology and rationale behind our revision and expansion of the PO to cover development stages for all plants, particularly the land plants (bryophytes through angiosperms). As a case study to illustrate the general approach, we examine variation in gene expression across embryo development stages in Arabidopsis and maize, demonstrating how the PO can be used to compare patterns of expression across stages and in developmentally different species. Although many genes appear to be active throughout embryo development, we identified a small set of uniquely expressed genes for each stage of embryo development and also between the two species. Evaluating the different sets of genes expressed during embryo development in Arabidopsis or maize may inform future studies of the divergent developmental pathways observed in monocotyledonous versus dicotyledonous species. The PO and its annotation database (http://www.planteome.org) make plant data for any species more discoverable and accessible through common formats, thus providing support for applications in plant pathology, image analysis, and comparative development and evolution.

Complex reproductive secretions occur in all extant gymnosperm lineages: a proteomic survey of gymnosperm pollination drops.

KEY MESSAGE: Complex protein-containing reproductive secretions are a conserved trait amongst all extant gymnosperms; the pollination drops of most groups include carbohydrate-modifying enzymes and defence proteins. Pollination drops are aqueous secretions that receive pollen and transport it to the ovule interior in gymnosperms (Coniferales, Cycadales, Ginkgoales, Gnetales). Proteins are well established as components of pollination drops in conifers (Coniferales) and Ephedra spp. (Gnetales), but it is unknown whether proteins are also present in the pollination drops of cycads (Cycadales), Ginkgo (Ginkgoales), Gnetum (Gnetales), or in the pollination drops produced by sterile ovules occurring on pollen plants in the Gnetales. We used liquid chromatography-tandem mass spectrometry followed by database-derived protein identification to conduct proteomic surveys of pollination drops collected from: Ceratozamia hildae, Zamia furfuracea and Cycas rumphii (Cycadales); Ginkgo biloba (Ginkgoales); Gnetum gnemon and Welwitschia mirabilis, including pollination drops from both microsporangiate and ovulate plants (Gnetales). We identified proteins in all samples: C. hildae (61), Z. furfuracea (40), C. rumphii (9), G. biloba (57), G. gnemon ovulate (17) and sterile ovules from microsporangiate plants (25) and W. mirabilis fertile ovules (1) and sterile ovules from microsporangiate plants (138). Proteins involved in defence and carbohydrate modification occurred in the drops of most groups, indicating conserved functions for proteins in pollination drops. Our study demonstrates that all extant gymnosperm groups produce complex reproductive secretions containing proteins, an ancient trait that likely contributed to the evolutionary success of seed plants.

Monocot plastid phylogenomics, timeline, net rates of species diversification, the power of multi-gene analyses, and a functional model for the origin of monocots.

PREMISE OF THE STUDY: We present the first plastome phylogeny encompassing all 77 monocot families, estimate branch support, and infer monocot-wide divergence times and rates of species diversification. METHODS: We conducted maximum likelihood analyses of phylogeny and BAMM studies of diversification rates based on 77 plastid genes across 545 monocots and 22 outgroups. We quantified how branch support and ascertainment vary with gene number, branch length, and branch depth. KEY RESULTS: Phylogenomic analyses shift the placement of 16 families in relation to earlier studies based on four plastid genes, add seven families, date the divergence between monocots and eudicots+Ceratophyllum at 136 Mya, successfully place all mycoheterotrophic taxa examined, and support recognizing Taccaceae and Thismiaceae as separate families and Arecales and Dasypogonales as separate orders. Only 45% of interfamilial divergences occurred after the Cretaceous. Net species diversification underwent four large-scale accelerations in PACMAD-BOP Poaceae, Asparagales sister to Doryanthaceae, Orchidoideae-Epidendroideae, and Araceae sister to Lemnoideae, each associated with specific ecological/morphological shifts. Branch ascertainment and support across monocots increase with gene number and branch length, and decrease with relative branch depth. Analysis of entire plastomes in Zingiberales quantifies the importance of non-coding regions in identifying and supporting short, deep branches. CONCLUSIONS: We provide the first resolved, well-supported monocot phylogeny and timeline spanning all families, and quantify the significant contribution of plastome-scale data to resolving short, deep branches. We outline a new functional model for the evolution of monocots and their diagnostic morphological traits from submersed aquatic ancestors, supported by convergent evolution of many of these traits in aquatic Hydatellaceae (Nymphaeales).

Gymnosperms on the EDGE.

Driven by limited resources and a sense of urgency, the prioritization of species for conservation has been a persistent concern in conservation science. Gymnosperms (comprising ginkgo, conifers, cycads, and gnetophytes) are one of the most threatened groups of living organisms, with 40% of the species at high risk of extinction, about twice as many as the most recent estimates for all plants (i.e. 21.4%). This high proportion of species facing extinction highlights the urgent action required to secure their future through an objective prioritization approach. The Evolutionary Distinct and Globally Endangered (EDGE) method rapidly ranks species based on their evolutionary distinctiveness and the extinction risks they face. EDGE is applied to gymnosperms using a phylogenetic tree comprising DNA sequence data for 85% of gymnosperm species (923 out of 1090 species), to which the 167 missing species were added, and IUCN Red List assessments available for 92% of species. The effect of different extinction probability transformations and the handling of IUCN data deficient species on the resulting rankings is investigated. Although top entries in our ranking comprise species that were expected to score well (e.g. Wollemia nobilis, Ginkgo biloba), many were unexpected (e.g. Araucaria araucana). These results highlight the necessity of using approaches that integrate evolutionary information in conservation science.

Crowds Replicate Performance of Scientific Experts Scoring Phylogenetic Matrices of Phenotypes.

Scientists building the Tree of Life face an overwhelming challenge to categorize phenotypes (e.g., anatomy, physiology) from millions of living and fossil species. This biodiversity challenge far outstrips the capacities of trained scientific experts. Here we explore whether crowdsourcing can be used to collect matrix data on a large scale with the participation of nonexpert students, or "citizen scientists." Crowdsourcing, or data collection by nonexperts, frequently via the internet, has enabled scientists to tackle some large-scale data collection challenges too massive for individuals or scientific teams alone. The quality of work by nonexpert crowds is, however, often questioned and little data have been collected on how such crowds perform on complex tasks such as phylogenetic character coding. We studied a crowd of over 600 nonexperts and found that they could use images to identify anatomical similarity (hypotheses of homology) with an average accuracy of 82% compared with scores provided by experts in the field. This performance pattern held across the Tree of Life, from protists to vertebrates. We introduce a procedure that predicts the difficulty of each character and that can be used to assign harder characters to experts and easier characters to a nonexpert crowd for scoring. We test this procedure in a controlled experiment comparing crowd scores to those of experts and show that crowds can produce matrices with over 90% of cells scored correctly while reducing the number of cells to be scored by experts by 50%. Preparation time, including image collection and processing, for a crowdsourcing experiment is significant, and does not currently save time of scientific experts overall. However, if innovations in automation or robotics can reduce such effort, then large-scale implementation of our method could greatly increase the collective scientific knowledge of species phenotypes for phylogenetic tree building. For the field of crowdsourcing, we provide a rare study with ground truth, or an experimental control that many studies lack, and contribute new methods on how to coordinate the work of experts and nonexperts. We show that there are important instances in which crowd consensus is not a good proxy for correctness.

The Physcomitrella patens chromosome-scale assembly reveals moss genome structure and evolution.

The draft genome of the moss model, Physcomitrella patens, comprised approximately 2000 unordered scaffolds. In order to enable analyses of genome structure and evolution we generated a chromosome-scale genome assembly using genetic linkage as well as (end) sequencing of long DNA fragments. We find that 57% of the genome comprises transposable elements (TEs), some of which may be actively transposing during the life cycle. Unlike in flowering plant genomes, gene- and TE-rich regions show an overall even distribution along the chromosomes. However, the chromosomes are mono-centric with peaks of a class of Copia elements potentially coinciding with centromeres. Gene body methylation is evident in 5.7% of the protein-coding genes, typically coinciding with low GC and low expression. Some giant virus insertions are transcriptionally active and might protect gametes from viral infection via siRNA mediated silencing. Structure-based detection methods show that the genome evolved via two rounds of whole genome duplications (WGDs), apparently common in mosses but not in liverworts and hornworts. Several hundred genes are present in colinear regions conserved since the last common ancestor of plants. These syntenic regions are enriched for functions related to plant-specific cell growth and tissue organization. The P. patens genome lacks the TE-rich pericentromeric and gene-rich distal regions typical for most flowering plant genomes. More non-seed plant genomes are needed to unravel how plant genomes evolve, and to understand whether the P. patens genome structure is typical for mosses or bryophytes.

The Planteome database: an integrated resource for reference ontologies, plant genomics and phenomics.

The Planteome project (http://www.planteome.org) provides a suite of reference and species-specific ontologies for plants and annotations to genes and phenotypes. Ontologies serve as common standards for semantic integration of a large and growing corpus of plant genomics, phenomics and genetics data. The reference ontologies include the Plant Ontology, Plant Trait Ontology and the Plant Experimental Conditions Ontology developed by the Planteome project, along with the Gene Ontology, Chemical Entities of Biological Interest, Phenotype and Attribute Ontology, and others. The project also provides access to species-specific Crop Ontologies developed by various plant breeding and research communities from around the world. We provide integrated data on plant traits, phenotypes, and gene function and expression from 95 plant taxa, annotated with reference ontology terms. The Planteome project is developing a plant gene annotation platform; Planteome Noctua, to facilitate community engagement. All the Planteome ontologies are publicly available and are maintained at the Planteome GitHub site (https://github.com/Planteome) for sharing, tracking revisions and new requests. The annotated data are freely accessible from the ontology browser (http://browser.planteome.org/amigo) and our data repository.

Evolution of strigolactone receptors by gradual neo-functionalization of KAI2 paralogues.

BACKGROUND: Strigolactones (SLs) are a class of plant hormones that control many aspects of plant growth. The SL signalling mechanism is homologous to that of karrikins (KARs), smoke-derived compounds that stimulate seed germination. In angiosperms, the SL receptor is an α/ß-hydrolase known as DWARF14 (D14); its close homologue, KARRIKIN INSENSITIVE2 (KAI2), functions as a KAR receptor and likely recognizes an uncharacterized, endogenous signal ('KL'). Previous phylogenetic analyses have suggested that the KAI2 lineage is ancestral in land plants, and that canonical D14-type SL receptors only arose in seed plants; this is paradoxical, however, as non-vascular plants synthesize and respond to SLs. RESULTS: We have used a combination of phylogenetic and structural approaches to re-assess the evolution of the D14/KAI2 family in land plants. We analysed 339 members of the D14/KAI2 family from land plants and charophyte algae. Our phylogenetic analyses show that the divergence between the eu-KAI2 lineage and the DDK (D14/DLK2/KAI2) lineage that includes D14 occurred very early in land plant evolution. We show that eu-KAI2 proteins are highly conserved, and have unique features not found in DDK proteins. Conversely, we show that DDK proteins show considerable sequence and structural variation to each other, and lack clearly definable characteristics. We use homology modelling to show that the earliest members of the DDK lineage structurally resemble KAI2 and that SL receptors in non-seed plants likely do not have D14-like structure. We also show that certain groups of DDK proteins lack the otherwise conserved MORE AXILLARY GROWTH2 (MAX2) interface, and may thus function independently of MAX2, which we show is highly conserved throughout land plant evolution. CONCLUSIONS: Our results suggest that D14-like structure is not required for SL perception, and that SL perception has relatively relaxed structural requirements compared to KAI2-mediated signalling. We suggest that SL perception gradually evolved by neo-functionalization within the DDK lineage, and that the transition from KAI2-like to D14-like protein may have been driven by interactions with protein partners, rather than being required for SL perception per se.

Insights into the Evolution of Hydroxyproline-Rich Glycoproteins from 1000 Plant Transcriptomes.

The carbohydrate-rich cell walls of land plants and algae have been the focus of much interest given the value of cell wall-based products to our current and future economies. Hydroxyproline-rich glycoproteins (HRGPs), a major group of wall glycoproteins, play important roles in plant growth and development, yet little is known about how they have evolved in parallel with the polysaccharide components of walls. We investigate the origins and evolution of the HRGP superfamily, which is commonly divided into three major multigene families: the arabinogalactan proteins (AGPs), extensins (EXTs), and proline-rich proteins. Using motif and amino acid bias, a newly developed bioinformatics pipeline, we identified HRGPs in sequences from the 1000 Plants transcriptome project (www.onekp.com). Our analyses provide new insights into the evolution of HRGPs across major evolutionary milestones, including the transition to land and the early radiation of angiosperms. Significantly, data mining reveals the origin of glycosylphosphatidylinositol (GPI)-anchored AGPs in green algae and a 3- to 4-fold increase in GPI-AGPs in liverworts and mosses. The first detection of cross-linking (CL)-EXTs is observed in bryophytes, which suggests that CL-EXTs arose though the juxtaposition of preexisting SPn EXT glycomotifs with refined Y-based motifs. We also detected the loss of CL-EXT in a few lineages, including the grass family (Poaceae), that have a cell wall composition distinct from other monocots and eudicots. A key challenge in HRGP research is tracking individual HRGPs throughout evolution. Using the 1000 Plants output, we were able to find putative orthologs of Arabidopsis pollen-specific GPI-AGPs in basal eudicots.

Microbial-type terpene synthase genes occur widely in nonseed land plants, but not in seed plants.

The vast abundance of terpene natural products in nature is due to enzymes known as terpene synthases (TPSs) that convert acyclic prenyl diphosphate precursors into a multitude of cyclic and acyclic carbon skeletons. Yet the evolution of TPSs is not well understood at higher levels of classification. Microbial TPSs from bacteria and fungi are only distantly related to typical plant TPSs, whereas genes similar to microbial TPS genes have been recently identified in the lycophyte Selaginella moellendorffii The goal of this study was to investigate the distribution, evolution, and biochemical functions of microbial terpene synthase-like (MTPSL) genes in other plants. By analyzing the transcriptomes of 1,103 plant species ranging from green algae to flowering plants, putative MTPSL genes were identified predominantly from nonseed plants, including liverworts, mosses, hornworts, lycophytes, and monilophytes. Directed searching for MTPSL genes in the sequenced genomes of a wide range of seed plants confirmed their general absence in this group. Among themselves, MTPSL proteins from nonseed plants form four major groups, with two of these more closely related to bacterial TPSs and the other two to fungal TPSs. Two of the four groups contain a canonical aspartate-rich "DDxxD" motif. The third group has a "DDxxxD" motif, and the fourth group has only the first two "DD" conserved in this motif. Upon heterologous expression, representative members from each of the four groups displayed diverse catalytic functions as monoterpene and sesquiterpene synthases, suggesting these are important for terpene formation in nonseed plants.

Challenging the paradigms of leaf evolution: Class III HD-Zips in ferns and lycophytes.

Despite the extraordinary significance leaves have for life on Earth, their origin and development remain vigorously debated. More than a century of paleobotanical, morphological, and phylogenetic research has still not resolved fundamental questions about leaves. Developmental genetic data are sparse in ferns, and comparative studies of lycophytes and seed plants have reached opposing conclusions on the conservation of a leaf developmental program. We performed phylogenetic and expression analyses of a leaf developmental regulator (Class III HD-Zip genes; C3HDZs) spanning lycophytes and ferns. We show that a duplication and neofunctionalization of C3HDZs probably occurred in the ancestor of euphyllophytes, and that there is a common leaf developmental mechanism conserved between ferns and seed plants. We show C3HDZ expression in lycophyte and fern sporangia and show that C3HDZs have conserved expression patterns during initiation of lateral primordia (leaves or sporangia). This expression is maintained throughout sporangium development in lycophytes and ferns and indicates an ancestral role of C3HDZs in sporangium development. We hypothesize that there is a deep homology of all leaves and that a sporangium-specific developmental program was coopted independently for the development of lycophyte and euphyllophyte leaves. This provides molecular genetic support for a paradigm shift in theories of lycophyte leaf evolution.

Molecular phylogenetics of the Ronnbergia Alliance (Bromeliaceae, Bromelioideae) and insights into their morphological evolution.

The tank-epiphytic clade of berry-fruited bromeliads, also known as the Core Bromelioideae, represents a remarkable event of adaptive radiation within the Bromeliaceae; however, the details of this radiation have been difficult to study because this lineage is plagued with generic delimitation problems. In this study, we used a phylogenetic approach to investigate a well supported, albeit poorly understood, lineage nested within the Core Bromelioideae, here called the "Ronnbergia Alliance". In order to assess the monophyly and phylogenetic relationships of this group, we used three plastid and three nuclear DNA sequence markers combined with a broad sampling across three taxonomic groups and allied species of Aechmea expected to comprise the Ronnbergia Alliance. We combined the datasets to produce a well-supported and resolved phylogenetic hypothesis. Our main results indicated that the Ronnbergia Alliance was a well-supported monophyletic group, sister to the remaining Core Bromelioideae, and it was composed by species of the polyphyletic genera Aechmea, Hohenbergia and Ronnbergia. We identified two major internal lineages with high geographic structure within the Ronnbergia Alliance. The first of these lineages, called the Pacific Clade, contained species of Aechmea and Ronnbergia that occur exclusively from southern Central America to northwestern South America. The second clade, called the Atlantic Clade, contained species of Aechmea, Hohenbergia and Ronnbergia mostly limited to the Atlantic Forest and the Caribbean. We also explored the diagnostic and evolutionary importance of 13 selected characters using ancestral character reconstructions on the phylogenetic hypothesis. We found that the combination of tubular corollas apically spreading and unappendaged ovules had diagnostic value for the Ronnbergia Alliance, whereas flower size, length of the corolla tube, and petal pigmentation and apex were important characters to differentiate the Pacific and Atlantic clades. This study opens new perspectives for future taxonomic reorganizations and provides a framework for evolutionary and biogeographic studies.

Phylogenomics and historical biogeography of the monocot order Liliales: out of Australia and through Antarctica.

We present the first phylogenomic analysis of relationships among all ten families of Liliales, based on 75 plastid genes from 35 species in 29 genera, and 97 additional plastomes stratified across angiosperm lineages. We used a supermatrix approach to extend our analysis to 58 of 64 genera of Liliales, and calibrated the resulting phylogeny against 17 fossil dates to produce a new timeline for monocot evolution. Liliales diverged from other monocots 124 Mya and began splitting into separate families 113 Mya. Our data support an Australian origin for Liliales, with close relationships between three pairs of lineages (Corsiaceae/Campynemataceae, Philesiaceae/Ripogonaceae, tribes Alstroemerieae/Luzuriageae) in South America and Australia or New Zealand reflecting teleconnections of these areas via Antarctica. Long-distance dispersal (LDD) across the Pacific and Tasman Sea led to re-invasion of New Zealand by two lineages (Luzuriaga, Ripogonum); LDD allowed Campynemanthe to colonize New Caledonia after its submergence until 37 Mya. LDD permitted Colchicaceae to invade East Asia and Africa from Australia, and re-invade Africa from Australia. Periodic desert greening permitted Gloriosa and Iphigenia to colonize Southeast Asia overland from Africa, and Androcymbium-Colchicum to invade the Mediterranean from South Africa. Melanthiaceae and Liliaceae crossed the Bering land-bridge several times from the Miocene to the Pleistocene.