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.

Three-dimensional functional gradients direct stem curling in the resurrection plant Selaginella lepidophylla

Upon hydration and dehydration, the vegetative tissue of Selaginella lepidophylla can reversibly swell and shrink to generate complex morphological transformations. Here, we investigate how structural and compositional properties at tissue and cell wall levels in S. lepidophylla lead to different stem curling profiles between inner and outer stems. Our results show that directional bending in both stem types is associated with cross-sectional gradients of tissue density, cell orientation and secondary cell wall composition between adaxial and abaxial stem sides. In inner stems, longitudinal gradients of cell wall thickness and composition affect tip-to-base tissue swelling and shrinking, allowing for more complex curling as compared to outer stems. Together, these features yield three-dimensional functional gradients that allow the plant to reproducibly deform in predetermined patterns that vary depending on the stem type. This study is the first to demonstrate functional gradients at different hierarchical levels combining to operate in a three-dimensional context.

Functional study of the brassinosteroid biosynthetic genes from Selagnella moellendorfii in Arabidopsis

Brassinosteroids (BRs) are essential hormones for plant growth and development. Enzymes DET2 and CYP90 family are responsible for BR biosynthesis in seed plants. Yet, their roles in non-seed plants are unknown. Here, we report the first functional study of DET2 and all 4 CYP90 genes isolated from Selaginella moellendorfii. Sm89026 (SmCPD) belonged to a clade with CYP90A1 (CPD) and CYP90B1 (DWF4) while Sm182839, Sm233379 and Sm157387 formed a distinct clade with CYP90C1 (ROT3) and CYP90D1. SmDET2, SmCPD and Sm157387 were highly expressed in both leaves and strobili while Sm233379 was only highly expressed in the leaves but not strobili, implying their differential functions in a tissue-specific manner in S. moellendorfii. We showed that only SmDET2 and SmCPD completely rescued Arabidopsis det2 and cpd mutant phenotypes, respectively, suggestive of their conserved BR biosynthetic functions. However, neither SmCPD nor other CYP90 genes rescued any other cyp90 mutants. Yet overexpression of Sm233379 altered plant fertility and BR response, which means that Sm233379 is not an ortholog of any CYP90 genes in Arabidopsis but appears to have a BR function in the S. moellendorfii leaves. This function is likely turned off during the development of the strobili. Our results suggest a dramatic functional divergence of CYP90 family in the non-seed plants. While some of them are functionally similar to that of seed plants, the others may be functionally distinct from that of seed plants, shedding light for future exploration.

cis-Cinnamic Acid Is a Novel, Natural Auxin Efflux Inhibitor That Promotes Lateral Root Formation

Auxin steers numerous physiological processes in plants, making the tight control of its endogenous levels and spatiotemporal distribution a necessity. This regulation is achieved by different mechanisms, including auxin biosynthesis, metabolic conversions, degradation, and transport. Here, we introduce cis-cinnamic acid (c-CA) as a novel and unique addition to a small group of endogenous molecules affecting in planta auxin concentrations. c-CA is the photo-isomerization product of the phenylpropanoid pathway intermediate trans-CA (t-CA). When grown on c-CA-containing medium, an evolutionary diverse set of plant species were shown to exhibit phenotypes characteristic for high auxin levels, including inhibition of primary root growth, induction of root hairs, and promotion of adventitious and lateral rooting. By molecular docking and receptor binding assays, we showed that c-CA itself is neither an auxin nor an anti-auxin, and auxin profiling data revealed that c-CA does not significantly interfere with auxin biosynthesis. Single cell-based auxin accumulation assays showed that c-CA, and not t-CA, is a potent inhibitor of auxin efflux. Auxin signaling reporters detected changes in spatiotemporal distribution of the auxin response along the root of c-CA-treated plants, and long-distance auxin transport assays showed no inhibition of rootward auxin transport. Overall, these results suggest that the phenotypes of c-CA-treated plants are the consequence of a local change in auxin accumulation, induced by the inhibition of auxin efflux. This work reveals a novel mechanism how plants may regulate auxin levels and adds a novel, naturally occurring molecule to the chemical toolbox for the studies of auxin homeostasis.

Conserved Gene Expression Programs in Developing Roots from Diverse Plants

The molecular basis for the origin and diversification of morphological adaptations is a central issue in evolutionary developmental biology. Here, we defined temporal transcript accumulation in developing roots from seven vascular plants, permitting a genome-wide comparative analysis of the molecular programs used by a single organ across diverse species. The resulting gene expression maps uncover significant similarity in the genes employed in roots and their developmental expression profiles. The detailed analysis of a subset of 133 genes known to be associated with root development in Arabidopsis thaliana indicates that most of these are used in all plant species. Strikingly, this was also true for root development in a lycophyte (Selaginella moellendorffii), which forms morphologically different roots and is thought to have evolved roots independently. Thus, despite vast differences in size and anatomy of roots from diverse plants, the basic molecular mechanisms employed during root formation appear to be conserved. This suggests that roots evolved in the two major vascular plant lineages either by parallel recruitment of largely the same developmental program or by elaboration of an existing root program in the common ancestor of vascular plants.

Insights into the evolution and diversification of the AT-hook Motif Nuclear Localized gene family in land plants

BACKGROUND

Members of the ancient land-plant-specific transcription factor AT-Hook Motif Nuclear Localized (AHL) gene family regulate various biological processes. However, the relationships among the AHL genes, as well as their evolutionary history, still remain unexplored.

RESULTS

We analyzed over 500 AHL genes from 19 land plant species, ranging from the early diverging Physcomitrella patens and Selaginella to a variety of monocot and dicot flowering plants. We classified the AHL proteins into three types (Type-I/-II/-III) based on the number and composition of their functional domains, the AT-hook motif(s) and PPC domain. We further inferred their phylogenies via Bayesian inference analysis and predicted gene gain/loss events throughout their diversification. Our analyses suggested that the AHL gene family emerged in embryophytes and further evolved into two distinct clades, with Type-I AHLs forming one clade (Clade-A), and the other two types together diversifying in another (Clade-B). The two AHL clades likely diverged before the separation of Physcomitrella patens from the vascular plant lineage. In angiosperms, Clade-A AHLs expanded into 5 subfamilies; while, the ones in Clade-B expanded into 4 subfamilies. Examination of their expression patterns suggests that the AHLs within each clade share similar expression patterns with each other; however, AHLs in one monophyletic clade exhibit distinct expression patterns from the ones in the other clade. Over-expression of a Glycine max AHL PPC domain in Arabidopsis thaliana recapitulates the phenotype observed when over-expressing its Arabidopsis thaliana counterpart. This result suggests that the AHL genes from different land plant species may share conserved functions in regulating plant growth and development. Our study further suggests that such functional conservation may be due to conserved physical interactions among the PPC domains of AHL proteins.

CONCLUSIONS

Our analyses reveal a possible evolutionary scenario for the AHL gene family in land plants, which will facilitate the design of new studies probing their biological functions. Manipulating the AHL genes has been suggested to have tremendous effects in agriculture through increased seedling establishment, enhanced plant biomass and improved plant immunity. The information gleaned from this study, in turn, has the potential to be utilized to further improve crop production.

Brassinosteroids regulate plant growth through distinct signaling pathways in Selaginella and Arabidopsis

Brassinosteroids (BRs) are growth-promoting steroid hormones that regulate diverse physiological processes in plants. Most BR biosynthetic enzymes belong to the cytochrome P450 (CYP) family. The gene encoding the ultimate step of BR biosynthesis in Arabidopsis likely evolved by gene duplication followed by functional specialization in a dicotyledonous plant-specific manner. To gain insight into the evolution of BRs, we performed a genomic reconstitution of Arabidopsis BR biosynthetic genes in an ancestral vascular plant, the lycophyte Selaginella moellendorffii. Selaginella contains four members of the CYP90 family that cluster together in the CYP85 clan. Similar to known BR biosynthetic genes, the Selaginella CYP90s exhibit eight or ten exons and Selaginella produces a putative BR biosynthetic intermediate. Therefore, we hypothesized that Selaginella CYP90 genes encode BR biosynthetic enzymes. In contrast to typical CYPs in Arabidopsis, Selaginella CYP90E2 and CYP90F1 do not possess amino-terminal signal peptides, suggesting that they do not localize to the endoplasmic reticulum. In addition, one of the three putative CYP reductases (CPRs) that is required for CYP enzyme function co-localized with CYP90E2 and CYP90F1. Treatments with a BR biosynthetic inhibitor, propiconazole, and epi-brassinolide resulted in greatly retarded and increased growth, respectively. This suggests that BRs promote growth in Selaginella, as they do in Arabidopsis. However, BR signaling occurs through different pathways than in Arabidopsis. A sequence homologous to the Arabidopsis BR receptor BRI1 was absent in Selaginella, but downstream components, including BIN2, BSU1, and BZR1, were present. Thus, the mechanism that initiates BR signaling in Selaginella seems to differ from that in Arabidopsis. Our findings suggest that the basic physiological roles of BRs as growth-promoting hormones are conserved in both lycophytes and Arabidopsis; however, different BR molecules and BRI1-based membrane receptor complexes evolved in these plants.

Evolution of the bHLH genes involved in stomatal development: implications for the expansion of developmental complexity of stomata in land plants

Stomata play significant roles in plant evolution. A trio of closely related basic Helix-Loop-Helix (bHLH) subgroup Ia genes, SPCH, MUTE and FAMA, mediate sequential steps of stomatal development, and their functions may be conserved in land plants. However, the evolutionary history of the putative SPCH/MUTE/FAMA genes is still greatly controversial, especially the phylogenetic positions of the bHLH Ia members from basal land plants. To better understand the evolutionary pattern and functional diversity of the bHLH genes involved in stomatal development, we made a comprehensive evolutionary analysis of the homologous genes from 54 species representing the major lineages of green plants. The phylogenetic analysis indicated: (1) All bHLH Ia genes from the two basal land plants Physcomitrella and Selaginella were closely related to the FAMA genes of seed plants; and (2) the gymnosperm 'SPCH' genes were sister to a clade comprising the angiosperm SPCH and MUTE genes, while the FAMA genes of gymnosperms and angiosperms had a sister relationship. The revealed phylogenetic relationships are also supported by the distribution of gene structures and previous functional studies. Therefore, we deduce that the function of FAMA might be ancestral in the bHLH Ia subgroup. In addition, the gymnosperm "SPCH" genes may represent an ancestral state and have a dual function of SPCH and MUTE, two genes that could have originated from a duplication event in the common ancestor of angiosperms. Moreover, in angiosperms, SPCHs have experienced more duplications and harbor more copies than MUTEs and FAMAs, which, together with variation of the stomatal development in the entry division, implies that SPCH might have contributed greatly to the diversity of stomatal development. Based on the above, we proposed a model for the correlation between the evolution of stomatal development and the genes involved in this developmental process in land plants.

Conserved transport mechanisms but distinct auxin responses govern shoot patterning in Selaginella kraussiana

To provide a comparative framework to understand the evolution of auxin regulation in vascular plants, the effect of perturbed auxin homeostasis was examined in the lycophyte Selaginella kraussiana. Polar auxin transport was measured by tracing tritiated IAA in excised shoots. Shoots were cultured in the presence of auxin efflux inhibitors and exogenous auxin, and developmental abnormalities were documented. Auxin transport in Selaginella shoots is exclusively basipetal, as in angiosperms. Perturbed auxin transport results in the loss of meristem maintenance and abnormal shoot architecture. Dichotomous root branching in Selaginella appears to be regulated by an antagonistic relationship between auxin and cytokinin. The results suggest that basipetal polar auxin transport occurred in the common ancestor of lycophytes and euphyllophytes. Although the mechanisms of auxin transport appear to be conserved across all vascular plants, distinct auxin responses govern shoot growth and development in lycophytes and euphyllophytes.

ARABIDILLO gene homologues in basal land plants: species-specific gene duplication and likely functional redundancy

ARABIDILLO proteins regulate multicellular root development in Arabidopsis thaliana. Conserved ARABIDILLO homologues are present throughout land plants, even in early-evolving plants that do not possess complex root architecture, suggesting that ARABIDILLO genes have additional functions. Here, we have cloned and characterised ARABIDILLO gene homologues from two early-evolving land plants, the bryophyte Physcomitrella patens and the lycophyte Selaginella moellendorffii. We show that two of the PHYSCODILLO genes (PHYSCODILLO1A and -1B) exist as a tail-to-tail tandem array of two almost identical 12 kb sequences, while a third related gene (PHYSCODILLO2) is located elsewhere in the Physcomitrella genome. Physcomitrella possesses a very low percentage of tandemly arrayed genes compared with the later-evolving plants whose genomes have been sequenced to date. Thus, PHYSCODILLO1A and -1B genes represent a relatively unusual gene arrangement. PHYSCODILLO promoters are active largely in the haploid gametophyte, with additional activity at the foot of the sporophyte. The pattern of promoter activity is uniform in filamentous and leafy tissues, suggesting pleiotropic gene functions and likely functional redundancy: the latter possibility is confirmed by the lack of discernible phenotype in a physcodillo2 deletion mutant. Interestingly, the pattern of PHYSCODILLO promoter activity in female reproductive organs is strikingly similar to that of an Arabidopsis homologue, suggesting co-option of some PHYSCODILLO functions or regulation into both the sporophyte and gametophyte. In conclusion, our work identifies and characterises some of the earliest-evolving land plant ARABIDILLO homologues. We confirm that all land plant ARABIDILLO genes arose from a single common ancestor and suggest that PHYSCODILLO proteins have novel and pleiotropic functions, some of which may be conserved in later-evolving plants.

The Selaginella genome identifies genetic changes associated with the evolution of vascular plants

Vascular plants appeared ~410 million years ago, then diverged into several lineages of which only two survive: the euphyllophytes (ferns and seed plants) and the lycophytes. We report here the genome sequence of the lycophyte Selaginella moellendorffii (Selaginella), the first nonseed vascular plant genome reported. By comparing gene content in evolutionarily diverse taxa, we found that the transition from a gametophyte- to a sporophyte-dominated life cycle required far fewer new genes than the transition from a nonseed vascular to a flowering plant, whereas secondary metabolic genes expanded extensively and in parallel in the lycophyte and angiosperm lineages. Selaginella differs in posttranscriptional gene regulation, including small RNA regulation of repetitive elements, an absence of the trans-acting small interfering RNA pathway, and extensive RNA editing of organellar genes.

Step-by-step acquisition of the gibberellin-DELLA growth-regulatory mechanism during land-plant evolution

Angiosperms (flowering plants) evolved relatively recently and are substantially diverged from early land plants (bryophytes, lycophytes, and others [1]). The phytohormone gibberellin (GA) adaptively regulates angiosperm growth via the GA-DELLA signaling mechanism [2-7]. GA binds to GA receptors (GID1s), thus stimulating interactions between GID1s and the growth-repressing DELLAs [8-12]. Subsequent 26S proteasome-mediated destruction of the DELLAs promotes growth [13-17]. Here we outline the evolution of the GA-DELLA mechanism. We show that the interaction between GID1 and DELLA components from Selaginella kraussiana (a lycophyte) is GA stimulated. In contrast, GID1-like (GLP1) and DELLA components from Physcomitrella patens (a bryophyte) do not interact, suggesting that GA-stimulated GID1-DELLA interactions arose in the land-plant lineage after the bryophyte divergence ( approximately 430 million years ago [1]). We further show that a DELLA-deficient P. patens mutant strain lacks the derepressed growth characteristic of DELLA-deficient angiosperms, and that both S. kraussiana and P. patens lack detectable growth responses to GA. These observations indicate that early land-plant DELLAs do not repress growth in situ. However, S. kraussiana and P. patens DELLAs function as growth-repressors when expressed in the angiosperm Arabidopsis thaliana. We conclude that the GA-DELLA growth-regulatory mechanism arose during land-plant evolution and via independent stepwise recruitment of GA-stimulated GID1-DELLA interaction and DELLA growth-repression functions.

Growth from two transient apical initials in the meristem of Selaginella kraussiana

A major transition in land plant evolution was from growth in water to growth on land. This transition necessitated major morphological innovations that were accompanied by the development of three-dimensional apical growth. In extant land plants, shoot growth occurs from groups of cells at the apex known as meristems. In different land plant lineages, meristems function in different ways to produce distinct plant morphologies, yet our understanding of the developmental basis of meristem function is limited to the most recently diverged angiosperms. To redress this balance, we have examined meristem function in the lycophyte Selaginella kraussiana. Using a clonal analysis, we show that S. kraussiana shoots are derived from the activity of two short-lived apical initials that facilitate the formation of four axes of symmetry in the shoot. Leaves are initiated from just two epidermal cells, and the mediolateral leaf axis is the first to be established. This pattern of development differs from that seen in flowering plants. These differences are discussed in the context of the development and evolution of diverse land plant forms.

Distinct developmental mechanisms reflect the independent origins of leaves in vascular plants

Vascular plants diverged more than 400 million years ago into two lineages, the lycophytes and the euphyllophytes . Leaf-like organs evolved independently in these two groups . Microphylls in lycophytes are hypothesized to have originated as lateral outgrowths of tissue that later became vascularized (the enation theory) or through the sterilization of sporangia (the sterilization hypothesis) . Megaphylls in euphyllophytes are thought to represent modified lateral branches . The fossil record also indicates that the seed plant megaphyll evolved uniquely in the ancestor of seed plants, independent of megaphylls in ferns, because seed plants evolved from leafless progymnosperm ancestors . Surprisingly, a recent study of KNOX and ARP gene expression in a lycophyte was reported to indicate recruitment of a similar mechanism for determinacy in both types of leaves . We examined the expression of Class III HD-Zip genes in the lycophyte Selaginella kraussiana and in two gymnosperms, Ginkgo and Pseudotsuga. Our data indicate that mechanisms promoting leaf initiation, vascularization, and polarity are quite different in lycophytes and seed plants, consistent with the hypotheses that megaphylls originated as lateral branches whereas microphylls originated as tissue outgrowths.

Deterministic cellular descendance and its relationship to the branching of plant organ axes

A double-wall map L-system, designated as S(5-5), was developed to simulate the cellular pattern found at the summit of shoot apices of Psilotum nudum. Commencing from a 3-sided autoreproductive founder cell, fives steps of simulation established a basic set of ten different cell types. Continuing the simulation beyond the fifth step revealed that, in addition to the regular production of new 3-sided cells, a group of autoreproductive 5-sided cells came into being. A close correspondence exists between the cells of the two-dimensional simulation and the two-dimensional cellular patterns found on the epidermis of the apices of Psilotum species. The 3-sided cells produced during the simulation correspond to the potentially organogenetic 3-sided cells that can be seen upon the apical surfaces. Successive generations of these 3-sided apical cells (which are actually 4-sided tetrahedral cells when viewed in three dimensions) and their immediate descendants are thought to be selected to organise the successive pairs of apices that bring about the repeated bifurcation of the Psilotum shoots. The 5-sided cells contribute to the cellular "pavements" which separate these pairs of organogenetic centres, each with their 3-sided apical cells. The cellular patterns simulated by the S(5-5) system may also correspond to the cellular patterns found on the surfaces of some other pteridophyte apices, including that of the rhizophores of Selaginella species. Tritiated-thymidine labelling of rhizophore apices revealed a group of nonproliferating cells that was associated with rhizophore bifurcation and which may correspond to a group of pavement cells. Nonproliferating cells, by regulating the siting of new organogenetic centres, may have evolved as an accompaniment to branching events such as the bifurcation of root and organ axes.

Stages in development of Selaginella diffusa megaspores

In mature megaspores of Selaginella diffusa (C. Presl) Spring the units of the exospore are ordered and become unordered toward the outer and inner surfaces. The exospore surface is coated with silica at maturity. The insertion of the future gap begins in early stages with formation of many minigaps within the inner part of the exospore distally. The mesospore, like the exospore, is resistant to the acetolysis reaction and can, thus, provisionally be considered to consist of sporopollenin. Unit structures within the outer part of the mesospore are unordered, but become ordered in the middle and inner parts. The inner surface of the mesospore appears verrucate. In maturing megaspores, the mesospore is mostly disintegrated and the inner exospore, which encapsulated the mesospore, remains as a somewhat isolated structure, and is again near the outer exospore. There are connecting strands across the gap between the inner surface of the outer exospore and the surface of the inner exospore. There are also spheres on the outer surface of the inner exospore.

Inhibition of plastid division by ampicillin in the pteridophyte Selaginella nipponica Fr. et Sav

We investigated the effect of the beta-lactam antibiotic, ampicillin, on plastid division in the pteridophyte Selaginella nipponica. Guard cells of plantlets treated with 1 mM ampicillin only often had one plastid, whereas guard cells of untreated plantlets had two to four plastids. We generated a S. nipponica cell culture system and used it to investigate the effects of ampicillin. Treatment with 1 mM ampicillin had no effect on cell division in culture. We classified cultured cells into four types based on the number of plastids they contained: one (Type I), two (Type II), three or four (Type III) and more than five (Type IV). After 3 d in culture, the percentage of each cell type (I-IV) was 29.5, 46.7, 20.9, and 1.9%, respectively. Subsequently, the percentage of Types III and IV increased gradually, reaching 61.9 and 11.4%, respectively, after 15 d in culture in the absence of ampicillin. When 1 mM ampicillin was added, there was a minimal increase in the number of Type III and IV cells, with high percentages of Type I and II cells (32.4 and 45.7%, respectively) after 15 d. These results suggest that ampicillin inhibits plastid division in S. nipponica.