Non-seed plants are emerging gene sources for agriculture and insect control proteins.

The non-seed plants (e.g., charophyte algae, bryophytes, and ferns) have multiple human uses, but their contributions to agriculture and research have lagged behind seed plants. While sharing broadly conserved biology with seed plants and the major crops, non-seed plants sometimes possess alternative molecular and physiological adaptations. These adaptations may guide crop improvements. One such area is the presence of multiple classes of insecticidal proteins found in non-seed plant genomes which are either absent or widely diverged in seed plants. There are documented uses of non-seed plants, and ferns for example have been used in human diets. Among the occasional identifiable toxins or antinutritive components present in non-seed plants, none include these insecticidal proteins. Apart from these discrete risk factors which can be addressed in the safety assessment, there should be no general safety concern about sourcing genes from non-seed plant species.

The first homosporous lycophyte genome revealed the association between the recent dynamic accumulation of LTR-RTs and genome size variation.

The contrasting genome size between homosporous and heterosporous plants is fascinating. Different from the heterosporous seed plants and mainly homosporous ferns, the lycophytes are either heterosporous (Isoetales and Selaginellales) or homosporous (Lycopodiales). Many lycophytes are the resource plants of Huperzine A (HupA) which is invaluable for treating Alzheimer's disease. For the seed-free vascular plants, several high-quality genomes of heterosporous Selaginella, homosporous ferns (maidenhair fern, monkey spider tree fern), and heterosporous ferns (Azolla) have been published and provided important insights into the origin and evolution of early land plants. However, the homosporous lycophyte genome has not been decoded. Here, we assembled the first homosporous lycophyte genome and conducted comparative genomic analyses by applying a reformed pipeline for filtering out non-plant sequences. The obtained genome size of Lycopodium clavatum is 2.30 Gb, distinguished in more than 85% repetitive elements of which 62% is long terminal repeat (LTR). This study disclosed a high birth rate and a low death rate of the LTR-RTs in homosporous lycophytes, but the opposite occurs in heterosporous lycophytes. we propose that the recent activity of LTR-RT is responsible for the immense genome size variation between homosporous and heterosporous lycophytes. By combing Ks analysis with a phylogenetic approach, we discovered two whole genome duplications (WGD). Morover, we identified all the five recognized key enzymes for the HupA biosynthetic pathway in the L. clavatum genome, but found this pathway incomplete in other major lineages of land plants. Overall, this study is of great importance for the medicinal utilization of lycophytes and the decoded genome data will be a key cornerstone to elucidate the evolution and biology of early vascular land plants.

Evapotranspiration-linked silica deposition in a basal tracheophyte plant (Lycopodiaceae: Lycopodiella alopecuroides): implications for the evolutionary origins of phytoliths.

Phytoliths, microscopic deposits of hydrated silica within plants, play a myriad of functional roles in extant tracheophytes - yet their evolutionary origins and the original selective pressures leading to their deposition remain poorly understood. To gain new insights into the ancestral condition of tracheophyte phytolith production and function, phytolith content was intensively assayed in a basal, morphologically conserved tracheophyte: the foxtail clubmoss Lycopodiella alopecuroides. Wet ashing was employed to perform phytolith extractions from every major anatomical region of L. alopecuroides. Phytolith occurrence was recorded, alongside abundance, morphometric information, and morphological descriptions. Phytoliths were recovered exclusively from the microphylls, which were apicodistally silicified into multiphytolith aggregates. Phytolith aggregates were larger and more numerous in anatomical regions engaging in greater evapotranspirational activity. The tissue distribution of L. alopecuroides phytoliths is inconsistent with the expectations of proposed adaptive hypotheses of phytolith evolutionary origin. Instead, it is hypothesized that phytoliths may have arisen incidentally in the L. alopecuroides-like ancestral plant, polymerizing from intraplant silicon accumulations arising via bulk flow and 'leaky' cellular micronutrient channels. This basal, nonadaptive phytolith formation model would provide the evolutionary 'raw material' for later modification into the useful, adaptative, phytolith deposits seen in later-diverging plant clades.

Impacts of methyl jasmonate on Selaginella martensii: volatiles, transcriptomics, phytohormones, and gas exchange.

Methyl jasmonate (MeJA) induces various defence responses in seed plants, but for early plant lineages, information on the potential of jasmonates to elicit stress signalling and trigger physiological modifications is limited. The spikemoss Selaginella martensii was exposed to a range of MeJA concentrations (0, 10, 25, and 50 mM), and biogenic volatile organic compound (BVOC) emissions, photosynthetic rate (A), and stomatal conductance (gs) were continuously measured. In addition, changes in phytohormone concentrations and gene expression were studied. Enhancement of methanol, lipoxygenase pathway volatiles and linalool emissions, and reductions in A and gs, were MeJA dose-dependent. Before MeJA treatment, the concentration of 12-oxo-phytodienoic acid (OPDA) was 7-fold higher than jasmonic acid (JA). MeJA treatment rapidly increased OPDA and JA concentrations (within 30 min), with the latter more responsive. Some genes involved in BVOC biosynthesis and OPDA-specific response were up-regulated at 30 min after MeJA spraying, whereas those in the JA signalling pathway were not affected. Although JA was synthesized in S. martensii, OPDA was prioritized as a signalling molecule upon MeJA application. MeJA inhibited primary and enhanced secondary metabolism; we propose that fast-emitted linalool could serve as a marker of elicitation of stress-induced metabolism in lycophytes.

Plastome structure, evolution, and phylogeny of Selaginella.

As one of the earliest land plant lineages, Selaginella is important for studying land plant evolution. It is the largest genus of lycophytes containing 700-800 species. Some unique characters of Selaginella plastomes have been reported, but based only on 20 species. There have been no plastome phylogenies of Selaginella based on a relatively large sampling, and no efforts have been made to resolve the phylogeny of the enigmatic Sinensis group whose relationships have been unclear based on small datasets. Here we investigated the structures of 59 plastomes representing 51 species covering all six subgenera and 18 sections of Selaginella except two sections and including the intriguing Sinensis group for the first time. Our major results include: (1) the plastome size of Selaginella ranges tremendously from 78,492 bp to 187,632 bp; (2) there are numerous gene losses in Selaginella comparing with other lycophytes, Isoëtaceae and Lycopodiaceae; (3) the gene contents and plastome structures in Selaginella vary lineage-specifically and all infrageneric taxa are well supported in the plastome phylogeny; (4) the ndh gene family tends to lose or pseudogenize in those species with DR structure and without other short or medium repeats; (5) the short and medium repeat regions in SC mediate many conformations causing diverse and complex plastome structures, and six new conformations are discovered; (6) forty-eight species sampled have high GC content (>50%) but three species in the Sinensis group have âˆ¼ 30% GC content in plastomes, similar to most vascular plants; (7) the Sinensis group is monophyletic, includes at least two subgroups, and has the smallest plastomes in land plants except some parasitic plants, and their plastomes do not contain any tRNAs; (8) the younger lineages in Selaginella tend to have higher GC content, whereas the older lineages tend to have lower GC content; and (9) because of incomplete genomic data and abnormal structures or some unknown reasons, even the concatenated plastomes could not well resolve the phylogenetic relationships in Selaginella with confidence, highlighting the difficulty in resolving the phylogeny and evolution of this particularly important land plant lineage.

Biochemical Characterization and Function of Eight Microbial Type Terpene Synthases from Lycophyte Selaginella moellendorffii.

Selaginella moellendorffii is a lycophyte, a member of an ancient vascular plant lineage. Two distinct types of terpene synthase (TPS) genes were identified from this species, including S. moellendorffii TPS genes (SmTPSs) and S. moellendorffii microbial TPS-like genes (SmMTPSLs). The goal of this study was to investigate the biochemical functions of SmMTPSLs. Here, eight full-length SmMTPSL genes (SmMTPSL5, -15, -19, -23, -33, -37, -46, and -47) were functionally characterized from S. moellendorffii. Escherichia coli-expressed recombinant SmMTPSLs were tested for monoterpenes synthase and sesquiterpenes synthase activities. These enzymatic products were typical monoterpenes and sesquiterpenes that have been previous shown to be generated by typical plant TPSs when provided with geranyl diphosphate (GPP) and farnesyl diphosphate (FPP) as the substrates. Meanwhile, SmMTPSL23, -33, and -37 were up-regulated when induced by alamethicin (ALA) and methyl jasmonate (MeJA), suggesting a role for these genes in plants response to abiotic stresses. Furthermore, this study pointed out that the terpenoids products of SmMTPSL23, -33, and -37 have an antibacterial effect on Pseudomonas syringae pv. tomato DC3000 and Staphylococcus aureus. Taken together, these results provide more information about the catalytic and biochemical function of SmMTPSLs in S. moellendorffii plants.

An Evolutionarily Primitive and Distinct Auxin Metabolism in the Lycophyte Selaginella moellendorffii.

Auxin is a key regulator of plant growth and development. Indole-3-acetic acid (IAA), a plant auxin, is mainly produced from tryptophan via indole-3-pyruvate (IPA) in both bryophytes and angiosperms. Angiosperms have multiple, well-documented IAA inactivation pathways, involving conjugation to IAA-aspartate (IAA-Asp)/glutamate by the GH3 auxin-amido synthetases, and oxidation to 2-oxindole-3-acetic acid (oxIAA) by the DAO proteins. However, IAA biosynthesis and inactivation processes remain elusive in lycophytes, an early lineage of spore-producing vascular plants. In this article, we studied IAA biosynthesis and inactivation in the lycophyte Selaginella moellendorffii. We demonstrate that S. moellendorffii mainly produces IAA from the IPA pathway for the regulation of root growth and response to high temperature, similar to the angiosperm Arabidopsis. However, S. moellendorffii exhibits a unique IAA metabolite profile with high IAA-Asp and low oxIAA levels, distinct from Arabidopsis and the bryophyte Marchantia polymorpha, suggesting that the GH3 family is integral for IAA homeostasis in the lycophytes. The DAO homologs in S. moellendorffii share only limited similarity to the well-characterized rice and Arabidopsis DAO proteins. We therefore suggest that these enzymes may have a limited role in IAA homeostasis in S. moellendorffii compared to angiosperms. We provide new insights into the functional diversification of auxin metabolic genes in the evolution of land plants.

Drought-induced lacuna formation in the stem causes hydraulic conductance to decline before xylem embolism in Selaginella.

Lycophytes are the earliest diverging extant lineage of vascular plants, sister to all other vascular plants. Given that most species are adapted to ever-wet environments, it has been hypothesized that lycophytes, and by extension the common ancestor of all vascular plants, have few adaptations to drought. We investigated the responses to drought of key fitness-related traits such as stomatal regulation, shoot hydraulic conductance (Kshoot ) and stem xylem embolism resistance in Selaginella haematodes and S. pulcherrima, both native to tropical understory. During drought stomata in both species were found to close before declines in Kshoot , with a 50% loss of Kshoot occurring at -1.7 and -2.5 MPa in S. haematodes and S. pulcherrima, respectively. Direct observational methods revealed that the xylem of both species was resistant to embolism formation, with 50% of embolized xylem area occurring at -3.0 and -4.6 MPa in S. haematodes and S. pulcherrima, respectively. X-ray microcomputed tomography images of stems revealed that the decline in Kshoot occurred with the formation of an air-filled lacuna, disconnecting the central vascular cylinder from the cortex. We propose that embolism-resistant xylem and large capacitance, provided by collapsing inner cortical cells, is essential for Selaginella survival during water deficit.

The enigma of sex allocation in Selaginella.

Background and Aims: The division of resource investment between male and female functions is poorly known for land plants other than angiosperms. The ancient lycophyte genus Selaginella is similar in some ways to angiosperms (in heterospory and in having sex allocation occur in the sporophyte generation, for example) but lacks the post-fertilization maternal investments that angiosperms make via fruit and seed tissues. One would therefore expect Selaginella to have sex allocation values less female-biased than in flowering plants and closer to the theoretical prediction of equal investment in male and female functions. Nothing is currently known of sex allocation in the genus, so even the simplest predictions have not been tested. Methods: Volumetric measurements of microsporangial and megasporangial investment were made in 14 species of Selaginella from four continents. In five of these species the length of the main above-ground axis of each plant was measured to determine whether sex allocation is related to plant size. Key Results: Of the 14 species, 13 showed male-biased allocations, often extreme, in population means and among the great majority of individual plants. There was some indication from the five species with axis length measurements that relative male allocation might be related to the release height of spores, but this evidence is preliminary. Conclusions: Sex allocation in Selaginella provides a phylogenetic touchstone showing how the innovations of fruit and seed investment in the angiosperm life cycle lead to typically female-biased allocations in that lineage. Moreover, the male bias we found in Selaginella requires an evolutionary explanation. The bias was often greater than what would occur from the mere absence of seed and fruit investments, and thus poses a challenge to sex allocation theory. It is possible that differences between microspores and megaspores in their dispersal ecology create selective effects that favour male-biased sexual allocation. This hypothesis remains tentative.

Purification, cDNA cloning, and characterization of plant chitinase with a novel domain combination from lycophyte Selaginella doederleinii.

Chitinase-A from a lycophyte Selaginella doederleinii (SdChiA), having molecular mass of 53 kDa, was purified to homogeneity by column chromatography. The cDNA encoding SdChiA was cloned by rapid amplification of cDNA ends and polymerase chain reaction. It consisted of 1477 nucleotides and its open reading frame encoded a polypeptide of 467 amino acid residues. The deduced amino acid sequence indicated that SdChiA consisted of two N-terminal chitin-binding domains and a C-terminal plant class V chitinase catalytic domain, belonging to the carbohydrate-binding module family 18 (CBM18) and glycoside hydrolase family 18 (GH18), respectively. SdChiA had chitin-binding ability. The time-dependent cleavage pattern of (GlcNAc)4 by SdChiA showed that SdChiA specifically recognizes the ß-anomer in the + 2 subsite of the substrate (GlcNAc)4 and cleaves the glycoside bond at the center of the substrate. This is the first report of the occurrence of a family 18 chitinase containing CBM18 chitin-binding domains. ABBREVIATIONS: AtChiC: Arabidopsis thaliana class V chitinase; CBB: Coomassie brilliant blue R250; CBM: carbohydrate binding module family; CrChi-A: Cycas revolute chitinase-A; EaChiA: Equisetum arvense chitinase-A; GH: glycoside hydrolase family, GlxChi-B: gazyumaru latex chitinase-B; GlcNAc: N-acetylglucosamine; HPLC: high performance liquid chromatography; LysM; lysin motif; MtNFH1: Medicago truncatula ecotypes R108-1 chitinase; NCBI: national center for biotechnology information; NF: nodulation factor; NtChiV: Nicotiana tabacum class V chitinase; PCR: polymerase chain reaction; PrChi-A: Pteris ryukyuensis chitinase-A; RACE: rapid amplification of cDNA ends; SDS-PAGE: sodium dodecyl sulfate-polyacrylamide gel electrophoresis; SdChiA: Selaginella doederleinii chitinase-A.

Evolution of allosteric regulation in chorismate mutases from early plants.

Plants, fungi, and bacteria synthesize the aromatic amino acids: l-phenylalanine, l-tyrosine, and l-tryptophan. Chorismate mutase catalyzes the branch point reaction of phenylalanine and tyrosine biosynthesis to generate prephenate. In Arabidopsis thaliana, there are two plastid-localized chorismate mutases that are allosterically regulated (AtCM1 and AtCM3) and one cytosolic isoform (AtCM2) that is unregulated. Previous analysis of plant chorismate mutases suggested that the enzymes from early plants (i.e. bryophytes/moss, lycophytes, and basal angiosperms) formed a clade distinct from the isoforms found in flowering plants; however, no biochemical information on these enzymes is available. To understand the evolution of allosteric regulation in plant chorismate mutases, we analyzed a basal lineage of plant enzymes homologous to AtCM1 based on sequence similarity. The chorismate mutases from the moss/bryophyte Physcomitrella patens (PpCM1 and PpCM2), the lycophyte Selaginella moellendorffii (SmCM), and the basal angiosperm Amborella trichopoda (AmtCM1 and AmtCM2) were characterized biochemically. Tryptophan was a positive effector for each of the five enzymes examined. Histidine was a weak positive effector for PpCM1 and AmtCM1. Neither tyrosine nor phenylalanine altered the activity of SmCM; however, tyrosine was a negative regulator of the other four enzymes. Phenylalanine down-regulates both moss enzymes and AmtCM2. The 2.0 ŠX-ray crystal structure of PpCM1 in complex with the tryptophan identified the allosteric effector site and reveals structural differences between the R- (more active) and T-state (less active) forms of plant chorismate mutases. Molecular insight into the basal plant chorismate mutases guides our understanding of the evolution of allosteric regulation in these enzymes.

Four hundred million years of silica biomineralization in land plants.

Biomineralization plays a fundamental role in the global silicon cycle. Grasses are known to mobilize significant quantities of Si in the form of silica biominerals and dominate the terrestrial realm today, but they have relatively recent origins and only rose to taxonomic and ecological prominence within the Cenozoic Era. This raises questions regarding when and how the biological silica cycle evolved. To address these questions, we examined silica abundances of extant members of early-diverging land plant clades, which show that silica biomineralization is widespread across terrestrial plant linages. Particularly high silica abundances are observed in lycophytes and early-diverging ferns. However, silica biomineralization is rare within later-evolving gymnosperms, implying a complex evolutionary history within the seed plants. Electron microscopy and X-ray spectroscopy show that the most common silica-mineralized tissues include the vascular system, epidermal cells, and stomata, which is consistent with the hypothesis that biomineralization in plants is frequently coupled to transpiration. Furthermore, sequence, phylogenetic, and structural analysis of nodulin 26-like intrinsic proteins from diverse plant genomes points to a plastic and ancient capacity for silica accumulation within terrestrial plants. The integration of these two comparative biology approaches demonstrates that silica biomineralization has been an important process for land plants over the course of their >400 My evolutionary history.

Global transcriptome analysis reveals extensive gene remodeling, alternative splicing and differential transcription profiles in non-seed vascular plant Selaginella moellendorffii.

BACKGROUND: Selaginella moellendorffii, a lycophyte, is a model plant to study the early evolution and development of vascular plants. As the first and only sequenced lycophyte to date, the genome of S. moellendorffii revealed many conserved genes and pathways, as well as specialized genes different from flowering plants. Despite the progress made, little is known about long noncoding RNAs (lncRNA) and the alternative splicing (AS) of coding genes in S. moellendorffii. Its coding gene models have not been fully validated with transcriptome data. Furthermore, it remains important to understand whether the regulatory mechanisms similar to flowering plants are used, and how they operate in a non-seed primitive vascular plant. RESULTS: RNA-sequencing (RNA-seq) was performed for three S. moellendorffii tissues, root, stem, and leaf, by constructing strand-specific RNA-seq libraries from RNA purified using RiboMinus isolation protocol. A total of 176 million reads (44 Gbp) were obtained from three tissue types, and were mapped to S. moellendorffii genome. By comparing with 22,285 existing gene models of S. moellendorffii, we identified 7930 high-confidence novel coding genes (a 35.6% increase), and for the first time reported 4422 lncRNAs in a lycophyte. Further, we refined 2461 (11.0%) of existing gene models, and identified 11,030 AS events (for 5957 coding genes) revealed for the first time for lycophytes. Tissue-specific gene expression with functional implication was analyzed, and 1031, 554, and 269 coding genes, and 174, 39, and 17 lncRNAs were identified in root, stem, and leaf tissues, respectively. The expression of critical genes for vascular development stages, i.e. formation of provascular cells, xylem specification and differentiation, and phloem specification and differentiation, was compared in S. moellendorffii tissues, indicating a less complex regulatory mechanism in lycophytes than in flowering plants. The results were further strengthened by the evolutionary trend of seven transcription factor families related to vascular development, which was observed among four representative species of seed and non-seed vascular plants, and nonvascular land and aquatic plants. CONCLUSIONS: The deep RNA-seq study of S. moellendorffii discovered extensive new gene contents, including novel coding genes, lncRNAs, AS events, and refined gene models. Compared to flowering vascular plants, S. moellendorffii displayed a less complexity in both gene structure, alternative splicing, and regulatory elements of vascular development. The study offered important insight into the evolution of vascular plants, and the regulation mechanism of vascular development in a non-seed plant.

The evolution of lycopsid rooting structures: conservatism and disparity.

Contents 538 I. 538 II. 539 III. 541 IV. 542 543 References 543 SUMMARY: The evolution of rooting structures was a crucial event in Earth's history, increasing the ability of plants to extract water, mine for nutrients and anchor above-ground shoot systems. Fossil evidence indicates that roots evolved at least twice among vascular plants, in the euphyllophytes and independently in the lycophytes. Here, we review the anatomy and evolution of lycopsid rooting structures. Highlighting recent discoveries made with fossils we suggest that the evolution of lycopsid rooting structures displays two contrasting patterns - conservatism and disparity. The structures termed roots have remained structurally similar despite hundreds of millions of years of evolution - an example of remarkable conservatism. By contrast, and over the same time period, the organs that give rise to roots have diversified, resulting in the evolution of numerous novel and disparate organs.

Vascularization of the Selaginella rhizophore: anatomical fingerprints of polar auxin transport with implications for the deep fossil record.

The Selaginella rhizophore is a unique and enigmatic organ whose homology with roots, shoots, or neither of the two remains unresolved. Nevertheless, rhizophore-like organs have been documented in several fossil lycophytes. Here we test the homology of these organs through comparisons with the architecture of rhizophore vascularization in Selaginella. We document rhizophore vascularization in nine Selaginella species using cleared whole-mounts and histological sectioning combined with three-dimensional reconstruction. Three patterns of rhizophore vascularization are present in Selaginella and each is comparable to those observed in rhizophore-like organs of fossil lycophytes. More compellingly, we found that all Selaginella species sampled exhibit tracheids that arc backward from the stem and side branch into the rhizophore base. This tracheid curvature is consistent with acropetal auxin transport previously documented in the rhizophore and is indicative of the redirection of basipetal auxin from the shoot into the rhizophore during development. The tracheid curvature observed in Selaginella rhizophores provides an anatomical fingerprint for the patterns of auxin flow that underpin rhizophore development. Similar tracheid geometry may be present and should be searched for in fossils to address rhizophore homology and the conservation of auxin-related developmental mechanisms from early stages of lycophyte evolution.

Evaluating the spore genome sizes of ferns and lycophytes: a flow cytometry approach.

Ferns and lycophytes produce spores to initiate the gametophyte stage for sexual reproduction. Approximately 10% of these seedless vascular plants are apomictic, and produce genomic unreduced spores. Genome size comparisons between spores and leaves are a reliable, and potentially easier way to determine their reproductive mode compared to traditional approaches. However, estimation of the spore genome sizes of these plants has not been attempted. We attempted to evaluate the spore genome sizes of ferns and lycophytes using flow cytometry, collected spores from selected species representing different spore physical properties and taxonomic groups, and sought to optimize bead-vortexing conditions. By evaluating the spore and sporophyte genome sizes, we examined whether reproductive modes could be ascertained from these flow cytometry results. We proposed two separate sets of optimized bead-vortexing conditions for the nuclear extraction of green and nongreen spores. We further successfully extracted spore nuclei of 19 families covering most orders, and the qualities and quantities of these extractions satisfied the C-value criteria. These evaluated genome sizes further supported the reproductive modes reported previously. In the current study, flow cytometry was used for the first time to evaluate the spore genome sizes of ferns and lycophytes. This use of spore flow cytometry provides a new, efficient approach to ascertaining the reproductive modes of these plants.

Chloroplast RNA editing going extreme: more than 3400 events of C-to-U editing in the chloroplast transcriptome of the lycophyte Selaginella uncinata.

RNA editing in chloroplasts and mitochondria of land plants differs significantly in abundance. For example, some 200-500 sites of cytidine-to-uridine RNA editing exist in flowering plant mitochondria as opposed to only 30-50 such C-to-U editing events in their chloroplasts. In contrast, we predicted significantly more chloroplast RNA editing for the protein-coding genes in the available complete plastome sequences of two species of the spike moss genus Selaginella (Lycopodiophyta). To evaluate these predictions we investigated the Selaginella uncinata chloroplast transcriptome. Our exhaustive cDNA studies identified the extraordinary number of 3415 RNA-editing events, exclusively of the C-to-U type, in the 74 mRNAs encoding intact reading frames in the S. uncinata chloroplast. We find the overwhelming majority (61%) of the 428 silent editing events leaving codon meanings unaltered directly neighboring other editing events, possibly suggesting a sterically more flexible RNA-editing deaminase activity in Selaginella. No evidence of RNA editing was found for tRNAs or rRNAs but we identified a total of 74 editing sites in cDNA sequences of four group II introns (petBi6g2, petDi8g2, ycf3i124g2, and ycf3i354g2) retained in partially matured transcripts, which strongly contribute to improved base-pairing in the intron secondary structures as a likely prerequisite for their splicing.

Root evolution at the base of the lycophyte clade: insights from an Early Devonian lycophyte.

BACKGROUND AND AIMS: The evolution of complex rooting systems during the Devonian had significant impacts on global terrestrial ecosystems and the evolution of plant body plans. However, detailed understanding of the pathways of root evolution and the architecture of early rooting systems is currently lacking. We describe the architecture and resolve the structural homology of the rooting system of an Early Devonian basal lycophyte. Insights gained from these fossils are used to address lycophyte root evolution and homology. METHODS: Plant fossils are preserved as carbonaceous compressions at Cottonwood Canyon (Wyoming), in the Lochkovian-Pragian (∼411 Ma; Early Devonian) Beartooth Butte Formation. We analysed 177 rock specimens and documented morphology, cuticular anatomy and structural relationships, as well as stratigraphic position and taphonomic conditions. KEY RESULTS: The rooting system of the Cottonwood Canyon lycophyte is composed of modified stems that bear fine, dichotomously branching lateral roots. These modified stems, referred to as root-bearing axes, are produced at branching points of the above-ground shoot system. Root-bearing axes preserved in growth position exhibit evidence of positive gravitropism, whereas the lateral roots extend horizontally. Consistent recurrence of these features in successive populations of the plant preserved in situ demonstrates that they represent constitutive structural traits and not opportunistic responses of a flexible developmental programme. CONCLUSIONS: This is the oldest direct evidence for a rooting system preserved in growth position. These rooting systems, which can be traced to a parent plant, include some of the earliest roots known to date and demonstrate that substantial plant-substrate interactions were under way by Early Devonian time. The morphological relationships between stems, root-bearing axes and roots corroborate evidence that positive gravitropism and root identity were evolutionarily uncoupled in lycophytes, and challenge the hypothesis that roots evolved from branches of the above-ground axial system, suggesting instead that lycophyte roots arose as a novel organ.